1. Introduction
Salinity is one of the most important environmental factors, and it greatly affects the survival, growth, development, reproduction, metabolism and other physiological activities of teleosts [1,2]. To cope with salinity changes, many euryhaline teleosts have evolved complex physiological strategies that involve a series of various biological processes [3,4,5]. With the development of high-throughput sequencing technology, RNA-Seq has been regarded as a powerful tool to interpret the functional genomic elements and uncover the molecular mechanisms in organisms [6,7].
Over the past few years, RNA-Seq analyses associated with salinity adaptation have been realized in some teleosts, such as striped catfish (Pangasianodon hypophthalmus) [7], tongue sole (Cynoglossus semilaevis) [8], spotted sea bass (Lateolabrax maculatus) [9], steelhead trout (Oncorhynchus mykiss) [10] and turbot (Scophthalmus maximus) [11]. However, most of these studies were only focused on the differentially expressed transcripts or genes after salinity changes. The knowledge regarding roles of post-transcriptional regulation mechanism in salinity adaptation has been quite limited until now.
Alternative splicing (AS) is a key post-transcriptional regulatory mechanism that produces various splicing variants from same pre-mRNA, making great contributions to increase the transcriptome complexity and enhance proteome diversity in eukaryotes [12,13]. Genome-wide studies have reported that ~95% of multi-exon genes in human (Homo sapiens) undergo AS events [14] compared with ~63% of genes in mouse (Mus musculus) [15] and ~61% in thale cress (Arabidopsis thaliana) [16]. Of them, some genes could even produce hundreds of spliced variants. For instance, there are 247 unique spliced variants derived from the neurexin-1-alpha gene in mouse [17].
Based on different splicing patterns, AS events can be roughly classified into five basic types, including exon skipping (ES), intron retention (IR), mutually exclusive exon (MXE), alternative 5′ splice site (A5SS) and alternative 3′ splice site (A3SS). In general, ES is the most frequent AS event in animals, while IR is the most prevalent type of AS events in plants [13]. The insertion or deletion of alternative fragments in mRNAs is likely to directly change the protein products and heavily affect the biological characteristics, such as the structure, function, stability and subcellular localization [13,18].
The two AS isoforms from the Bcl-x gene in fruit flies (Drosophila melanogaster) have been proven to participate in activating and inhibiting apoptosis [19]. AS could also modulate the abundance of mRNA via introduction of premature termination codons (PTCs). The mRNA with PTCs would be recognized by specific proteins and degraded by nonsense-mediated decay pathways in cells [20,21,22].
Mounting evidence has revealed that AS acts an important regulator in response to environmental stresses in teleosts. It was reported that the exon 9a of heat shock transcription factor 1 gene is differentially alternative spliced between control and heat-intolerant groups of channel catfish (Ictalurus punctatus) under heat stress, suggesting its potential roles in heat tolerance [23]. In rainbow trout (Oncorhynchus mykiss), differentially alternative splicing (DAS) events are significantly increased following heat exposure [24]. Hundreds of DAS events in the heart and brain of Nile tilapia (Oreochromis niloticus), enriched in the circadian clock pathway, are proven to be correlated with cold adaptation [25]. Similarly, differences in splicing patterns are observed for genes involved in striated muscle contraction and the mus cle myosin complex, which are thought to be key in the cold adaptation in skeletal muscle of zebrafish (Danio rerio), Atlantic killifish (Fundulus heteroclitus) and threespine stickleback (Gasterosteus aculeatus) (Healy & Schulte, 2019). Under acute hypoxia environment, elongation factor 2 kinase, NADH dehydrogenase (ubiquinone) Fe-S protein 1, family with sequence similarity 162 member A and NDRG family member 2 genes undergo differential usage of exons, which have been demonstrated to be tightly related to hypoxia tolerance in heart and gill of Nile tilapia [26,27].
Numerous AS events are greatly induced (over a 20% increase) by high-salinity environment in both the liver and gills of spotted sea bass [28]. This finding suggests that AS may be an important regulator in response to salinity changes in spotted sea bass. Nevertheless, the common mechanism of AS in salinity adaptation across teleosts has remained unknown. Turbot and tongue sole are typical euryhaline teleosts with the remarkable ability to survive in a wide range of salinity environments from 0 to 50 ppt [8,11,29,30,31,32]. Steelhead trout is also capable of adapting and of living in 0-32 ppt salinity environments [10,33]. Therefore, they provide excellent models to investigate molecular mechanism of salinity adaptation in teleosts. In the present study, AS analysis based on RNA-Seq datasets were performed to investigate the potential roles of various splicing patterns in salinity adaptation in the livers of turbot, tongue sole and steelhead trout. The study not only provides a new insight into uncovering the molecular mechanism of salinity adaptation in teleosts, but also serves as a valuable reference for related studies in the future.
2. Materials and Methods
2.1. Data Acquisition and Experimental Design
A total of 18 RNA-Seq datasets of livers in turbot, tongue sole and steelhead trout were derived from previously published studies, the objects of which were to investigate quantitative changes in gene expression in response to different salinity environments [8,10,11]. The salinity experiments of each study were briefly described as followed: healthy turbot were obtained and placed in six identical tanks (4000 L) full of seawater (salinity: 30 ± 0.2 ppt) for one-week acclimatization. Then, individuals from three tanks were directly exposed to dechlorinated tap water and regarded as the low-salinity group. The others were maintained in seawater and labeled as the high-salinity group. For both low- and high-salinity group, six individuals per tank were randomly collected and quickly anesthetized with MS-222 at 24 h after salinity changes. The livers were obtained and frozen in liquid nitrogen for following RNA extraction and transcriptome sequencing. For tongue sole, female individuals with similar sizes were selected and acclimatized in six experimental tanks at an optimal salinity of 30 ppt for one week. After that, six tanks were equally divided into low- and high-salinity groups. The salinity of high-salinity group were kept at 30 ppt, while the salinity of low-salinity group was gradually decreased from 30 to 15 ppt at a constant rate of 5 ppt per day for 3 days. After 60-day culture, three individuals per tank were euthanized and anesthetized with MS-222. The livers were immediately removed and frozen in liquid nitrogen for RNA extraction and transcriptome sequencing. For steelhead trout, healthy individuals were collected and reared in six cylindrical tanks under freshwater for two weeks prior to the salinity experiment. After that, low- and high-salinity groups were established in three replicated tanks. The low-salinity group was maintained in freshwater. Simultaneously, the salinity of high-salinity group was directly increased to 14 ppt, followed by a daily increase of 2 ppt until 30 ppt. The whole experiment lasted for 40 days. At the end, three individuals per tank were anesthetized with MS-222. Following anesthesia, livers were quickly dissected and immersed into liquid nitrogen for high-throughput sequencing.
2.2. Identification of AS and Differentially Alternative Splicing (DAS) Events
To construct the AS landscapes in the livers of turbot, tongue sole and steelhead trout, a total of 18 RNA-Seq datasets (six datasets from each fish) were selected and analyzed in the present study. The raw data of RNA-Seq derived from NCBI SRA database were converted to FASTQ format using SRA Toolkit v 2.9.0. Quality trimming and adapter clipping were performed to provide clean reads using fastp v0.23.2 with default parameters. The qualities of clean reads was further assessed using FastQC v0.11.9. The reference genomes and their corresponding annotation files of turbot, tongue sole and steelhead trout were also obtained from NCBI database with the assembly number of GCA_013347765.1, GCA_000523025.1 and GCA_002163495.1, respectively. High-quality clean reads of turbot, tongue sole and steelhead trout were then aligned to their corresponding reference genomes using Hisat2 v2.2.1 with default parameters. The aligned results were strictly filtered, and only unique alignments were retained for subsequent analysis. Reference-based transcript assembly was performed using StringTie v2.1.7. AStalavista v4.0 was utilized to automatically characterize the complete landscapes of AS events in the livers of the three species. The five typical AS events were defined by specific codes, including “0, 1-2^” for ES, “0, 1^2-” for IR, “1-2^, 3-4^” for MXE, “1^, 2^” for A5SS and “1-, 2-” for A3SS.
In addition, rMATs v.4.0.1 was also used to detect the AS events in each fish. Only AS events detected by both AStalavista and rMATs were considered stable for following analysis. Moreover, salt-responsive AS events were also determined using rMATs software by computing the inclusion level from two-group RNA-Seq data with replicates and constructing a hierarchical framework to model the variability among replicates. The threshold of DAS events was set as FDR-adjusted p-value < 0.05.
2.3. Identification of AS and Differentially Expressed Genes (DEGs) Events
Read matrix were constructed based on filtered alignments using featureCounts v 2.0.1. Then, read counts were transformed into fragments per kilobase of transcript per million fragments mapped (FPKM) for the normalization of gene expression levels. DESeq2 v1.34.0 R package was employed to determine the differentially expressed genes (DEGs) in the livers of all three species under different salinity environments. The significant thresholds of DEGs were set as |log2 (fold change)| ≥ 1 and p-value < 0.05.
2.4. Functional Annotation and Enrichment Analysis
The functional annotations of protein-coding genes in each fish were conducted using eggNOG-mapper v5.0 with diamond search mode. The functional categories of Gene Ontology (GO) terms, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were extracted and used to construct specific Org.db databases using AnnotationForge R package v3.14. Based on these Org.db databases, GO and KEGG enrichment analyses of turbot, tongue sole and steelhead trout were implemented using clusterProfiler v4.2.0 R package, respectively. Additionally, enrichment results were further visualized using the ggplot2 R package.
2.5. Orthologues Gene Identification
To obtain the putative orthologs among turbot, tongue sole and steelhead trout, a reciprocal best hit (RHB) strategy was performed using the BLAST v2.12.0 program. Based on the information of annotated gene model, only the longest transcripts of each gene were extracted and translated into amino acid sequences. Then, the predicated amino acid sequences of turbot were used as queries for the BLAST program with the parameters of e-value < 1 × 10−5. The RHB were performed for each pair of genes between tongue sole and turbot, steelhead trout and turbot. The putative orthologs were identified based on the top BLAST hits in each reciprocal BLAST pair.
3. Results
3.1. AS Landscapes in Livers of Turbot, Tongue Sole and Steelhead Trout
In total, 10,826, 10,741 and 10,112 AS events were identified in the livers of turbot, tongue sole and steelhead trout, respectively (Supplementary Table S1). These AS events were classified into five basic types, including ES, IR, MXE, A5SS and A3SS (Figure 1A). Among them, ES was the most abundant type, accounting for more than half of AS events in the livers of all three species, followed by A3SS, A5SS, IR and MXE (Figure 1B).
Additionally, we found that the 10,826, 10,741 and 10,112 AS events were derived from 5281, 5516 and 5614 functional genes, demonstrating that numerous genes may undergo at least two different AS events. The intersections between genes and their AS events in three species were intuitively visualized through the construction of UpSet plots, as shown in Figure 2A–C. The results clearly showed that most genes in turbot, tongue sole and steelhead trout only produced single AS events. Nevertheless, 64 genes in turbot, 58 genes in tongue sole and 42 genes in steelhead trout had the ability to produce greater than three different AS events.
These functional genes were closely associated with growth, signal transduction, the cell cycle and immune response. Strikingly, five, six and two genes even harbored all five different AS events in three species, including vegfa, pfkfb2b, nlrc3, p4ha2 and uncharacterized LOC118318343 in turbot; sox11, hcf1, rbms3, spo11, uncharacterized LOC103399889 and uncharacterized LOC112487504 in tongue sole; and sec31 and pps1 in steelhead trout. Their detailed information is summarized in Supplementary Table S2.
3.2. Distribution of AS Events across the Genomes of Turbot, Tongue Sole and Steelhead Trout
We investigated the distribution of AS events across the whole genomes of turbot, tongue sole and steelhead trout, which were further visualized using Circos plots (Figure 3A–C). To well characterize the distribution of AS events, the density of AS events (AS event number/gene number) was calculated within non-overlap 200 kb bins on the genomes of three species. The average densities of AS events in turbot, tongue sole and steelhead trout were estimated as 2.08, 1.94 and 1.79, respectively.
We also found that the AS event densities varied dramatically across the same chromosome, organized into numerous hotspots with the density of AS events greater than 5. Of the identified hotspots, each gene locating in 15.40–15.60 Mb on Chr12 and 18.60–18.80 Mb on Chr16 of turbot could generate 17 AS events on average. Similar observations were also detected in 13.60–13.80 Mb on Chr18 of tongue sole and 55.40–54.60 Mb on Chr25 of steelhead trout, the AS event densities of which were 12.66 and 11.67, respectively.
There were 109, 103 and 107 functional genes located in the AS hotspots of turbot, tongue sole and steelhead trout. KEGG enrichment analysis was carried out to explore the potential functions of these genes. The results indicated that genes in turbot were significantly enriched in apoptosis, folate biosynthesis and the complement and coagulation cascades pathways. The genes in tongue sole were closely associated with complement and coagulation cascades, and some disease-related pathways, such as staphylococcus aureus infection and legionellosis.
Genes in steelhead trout were tightly related to the complement and coagulation cascades, alpha-Linolenic acid metabolism and biosynthesis of unsaturated fatty acids. It is noted that genes located in AS hotspots of all three species were commonly enriched in complement and coagulation cascades, particularly complement cascades. The related elements in complement and coagulation cascades and functional gene information located in AS hotspots are visualized and given in Supplementary Figure S1 and Supplementary Table S3, respectively.
3.3. Identification of DAS Events under Different Salinity Environments
DAS events were determined by comparing the RNA-Seq datasets between the low- and high-salinity groups in the livers of all three species. As results, a total of 940, 590 and 553 DAS events were identified in turbot, tongue sole and steelhead trout, corresponding to 769, 482 and 467 functional genes, respectively (Table 1). Among them, 37, 56 and 23 genes showed differentially expression levels and differentially alternative splices in the livers of turbot, tongue sole and steelhead trout (Supplementary Figure S2; Supplementary Table S4).
In agreement with the classification patterns of all AS events, most DAS events in all three species belonged to the ES type, followed by A3SS, A5SS, IR and MXE (Table 1; Supplementary Figure S3). In addition, a pair of shorter and longer isoforms in AS events were defined as inclusion and exclusion by rMATS, respectively (Figure 1A). In the present study, exclusion and inclusion isoforms in each type of DAS events were systematically characterized according to the values of the relative inclusion levels (Figure 4). The results clearly showed that different salinity environments resulted in diverse patterns of exclusion and inclusion isoforms in most AS types, except for IR.
Global increased inclusion levels and decreased proportions of exclusion isoforms in IR were detected in the livers of all three species under the high-salinity environment (Figure 4). It was indicated that the high-salinity environment could elevate the intron retention levels and promote the selection of longer introns in some differentially spliced genes. Splicing models of several genes were constructed as shown in Supplementary Figure S4. The functional analysis showed that these genes with different IRs were tightly related to the regulation of mRNA 3′-end processing, cap mRNA methylation and protein dilapidation (Supplementary Figure S3).
In order to detect the common DAS genes across species under different salinity environments, the putative orthologs among turbot, tongue sole and steelhead trout were obtained using the RHB strategy. As result, a total of 10 common DAS genes were detected in all three species (Supplementary Figures S5–S7). The functional annotation revealed that these genes were involved in transcription regulation, RNA splicing, cytoskeleton regulation, signal transduction, neuronal development, protein transport, endoplasmic reticulum formation and hyaluronan metabolism (Supplementary Table S5).
3.4. Enrichment Analysis of DAS Genes under Different Salinity Environments
Enrichment analysis was performed to explore the potential functions of all the DAS genes in the livers of all three species. These GO terms and KEGG pathways, commonly enriched in at least two fish and ranked in top 30, were displayed in Figure 5. The results clearly indicated that many DAS genes were primarily enriched in some GO terms associated with the biological processes of RNA processing, such as mRNA binding, mRNA processing, mRNA splicing, via spliceosome, RNA splicing, regulation of mRNA metabolic, regulation of RNA splicing and so on (Figure 5A).
The KEGG enrichment analysis showed that complement and coagulation cascades, glycerophospholipid metabolism and herpes simplex virus 1 infection were also enriched in the livers of turbot, tongue sole and steelhead trout. Moreover, we found that these DAS genes were tightly related to the spliceosome pathway and were in parallel with the results of the GO enrichment analysis (Figure 5B). The detailed information of these enriched GO terms and KEGG pathways are given in Supplementary Tables S6 and S7.
3.5. Function Categories of Common DAS Genes Associated with RNA Processing
The GO and KEGG enrichment analyses provided considerable evidences that DAS genes associated with RNA processing were crucial in response to different salinity environments. Based on the orthologous relationships, these DAS genes were associated with RNA processing and existed in at least two fish. Their specific biology functions were further investigated by reviewing published documents. As shown in Figure 6, most DAS genes associated with RNA processing appear to be species-specific, while a total of common 12 DAS genes were detected in at least two fish, suggesting important roles in salinity regulation in teleosts.
A putative model of functional interaction of these DAS genes was proposed as in Figure 7. Of them, most DAS genes were involved in RNA splicing, including celf1, rmb39, ptbp1, ptbp2, hnrnpm, u2surp, srsf2a, srsf3a and srsf7a. Additionally, ptbp2 and celf1 were also implicated in translation regulation. Research reported that fubp3 could work for transcriptional regulation and elavl1 for mRNA stability. The specific biology functions of these common DAS genes would be discussed in the next section.
4. Discussion
Salinity adaptation is a complicated process that involves multiple physiological responses. In the present study, we conducted a comprehensive study to identify and characterize the AS events associated with salinity adaptation in the livers of turbot, tongue sole and steelhead trout, providing valuable information to better understand the underlying molecular mechanisms of salinity adaptation in teleosts. In fish, the liver is generally considered as an essential metabolic organ, supplying energy to other osmoregulatory organs [1,34]. Furthermore, it has been reported that the liver is also an important antioxidant defense system in response to environmental stresses [35].
The present study explored genome-wide AS events in the livers of three species based on the 18 RNA-Seq datasets. Thousands of AS events were identified and characterized, suggesting that AS could be a widespread phenomenon in the livers of the three species. Among these identified AS events, ES was the most frequent type, which is largely consistent with previous observation in teleosts, such as channel catfish [23], rainbow trout [24] and spotted sea bass [28]. As an important post-transcriptional regulation mechanism, AS events could generate different transcripts from the same pre-mRNA. Genes experiencing different AS events can create a much greater extent of transcript diversity [36]. These transcripts may be degraded by nonsense-mediated decay pathways to regulate the transcript abundance or encoding proteins with subtle functional differences that can have profound biological consequences [37]. In the present study, we found that several genes could produce all five AS events in the livers of all three fish, such as pfkfb2b in turbot, sox11 in tongue sole and sec31 in steelhead trout. As is known, pfkfb2b could work not only as a key regulator of glycolysis through catalyzing the synthesis and degradation of fructose-2, 6-bisphosphatase but also could affect the biological process of liver regeneration [38,39]. Sox11 generally acts as a transcriptional activator to enhance the transcriptional activities of some functional genes [40]. A recent study demonstrated that sox transcription factors are also implicated in liver development [41]. In addition, sec31 could function in vesicle budding and cargo export from the endoplasmic reticulum, participating in protein transport [42]. The diverse mRNA of these genes may coordinate physiologically meaningful changes in protein isoforms to maintain the normal biological functions of the liver in these fish, such as liver regeneration, development and protein metabolism. However, their specific molecular mechanisms required further study in the future. Additionally, functional genes with diverse AS events were also detected in other organisms, such as the Na+/K+/2Cl− cotransporter in humans [43], signal transducer and activator of transcription 3 in zebrafish [44], immunoglobulin heavy chain in spotted sea bass [9] and acetyl-CoA carboxylase alpha in swamp eels (Monopterus albus) [45].
Additionally, many hotspots of AS events were also detected across the genomes of turbot, tongue sole and steelhead trout. In all three species, the genes located in AS hotspot regions were commonly enriched in complement and coagulation cascades pathways, particularly the complement cascade system. The complement and coagulation systems are two distinct multi-component protein networks, both of which work as important innate defenses against external threats [46,47,48]. The activated complement system could elicit a pro-inflammatory response and recruit immune cells from both innate and adaptive branches of the immune system, playing crucial roles in the defense against pathogens [49,50,51]. Whereas coagulation is able to culminate the formation of thrombin and convert soluble fibrinogen to insoluble fibrin clots, acting as a main column in hemostasis in organisms [47,52]. Our results have revealed that many key elements, such as C2, C3, C5, C8, C1INH, FH, PLG and CR1, were located in the AS hotspot regions, suggesting that AS events frequently appear in these genes and produce a large number of alternative spliced transcripts. Furthermore, we found that the most abundant type of AS events in these elements is ES. Through changing the array of exons, AS events may affect the encoding sequences of transcripts and increase their proteomic diversity to modulate the functionality of complement and coagulation cascades.
Mounting evidence indicates that the AS mechanism is closely implicated in the response to environmental stimuli in teleosts, such as the heat tolerance of rainbow trout [24] and hypoxia adaptation of Nile tilapia [26]. In the present study, there were substantial changes in the splicing patterns of functional genes in the livers of all three species between low- and high-salinity environments, resulting in numerous DAS events. Of these DAS events, ES is the most common type as expected. IR is greatly induced by high-salinity environment in the livers of all three species. It is known that ES is the most-prevalent type of alternative splicing in animals, which could work as a major contributor to increase the complexity of proteomes [53,54]. Whereas IR is generally considered an aberrant splicing event with little functional consequence [55]. Through coupling to nonsense-mediated decay, IR can be used to govern gene expression levels [56]. Based on the functional analysis, ES and IR events may have affected these genes associated with mRNA binding, cytoskeleton organization, site selection of RNA transcription, mRNA methylation and protein dilapidation in the livers of three species under different salinity environments. Although the functional differences among alternative spliced transcripts remained largely unexplored, the models of the fubp1, srsf2a, mark3, cmtr1, fus, atg4b, ewsr1b and tom1l2 genes clearly showed ES and IR caused the frame shift mutations or inserted premature termination codes. This would result in different protein products or nonsense-mediated decay to influence gene functions in response to different salinity environments in the livers of three species.
The transcribed pre-mRNAs have to undergo intron removal and exon joints to form mature mRNA. Pre-mRNA splicing takes place in the nucleus carried out by spliceosome, a large molecular complex [13]. The spliceosome core is composed of five small nuclear ribonucleoproteins and many splice-related functional proteins, which assemble at introns in a precise order [57,58]. In addition to the core components of the spliceosome, some RNA-binding proteins also play important roles in the determination of splice-site selection in pre-mRNA, such as serine/arginine-rich (SR) splicing factors and heterogeneous nuclear ribonucleoproteins (hnRNPs) [59,60]. In the present study, we found that the many DAS genes in the livers of all three species under different salinity environments were tightly related to RNA processing, especially RNA splicing. This result is largely consistent with previous observations in the salinity adaptation of spotted sea bass [28], heat tolerance of rainbow trout [24] and bacteria defense of channel catfish [61]. There were two common RNA splicing regulators (celf1 and hnrnpm) with DAS events in the livers of all three species under different salinity environments. The functional protein, encoded by celf1, is implicated in mediating the inclusion and exclusion of some specific exons in functional genes [62,63]. In addition, hnrnpm is reported to be involved in the correct splicing of a subset of pre-mRNA [64]. The other common DAS genes identified in at least two species include ptbp1, ptbp2, rmb39, u2surp, srsf2a, srsf3a, srsf7a, fubp3 and elval1, which also act as important splicing factors and take part in the regulation of RNA splicing [63,64,65,66,67,68,69,70,71]. Different salinity environments changed the splicing patterns and resulted in different splicing transcripts of these RNA splicing regulators, which are likely to alter their functions in pre-mRNA splicing and affect the splicing decisions of downstream target genes in response to salinity changes. Although not enriched in GO terms or KEGG enrichments, many DAS genes associated with salinity adaptation were accompanied with differentially expressed levels under different salinity environments, such as sodium- and chloride-dependent taurine transporter and solute carrier family 12 member 7 for ion and amino transport, fatty acid desaturase 2 and acyl-coenzyme A thioesterase 1 for the energy metabolism. This indicates that RNA splicing regulators may regulate these functional genes through affecting their transcript abundance. The involvement of these genes indicates their importance in response to salinity changes.
5. Conclusions
In the present study, a total of 18 RNA-Seq datasets were used to investigate the potential roles of AS events in salinity adaptation in the livers of turbot, tongue sole and steelhead trout. A total of 10,826, 10,741 and 10,112 AS events were identified in the livers of turbot, tongue sole and steelhead trout, respectively. The characteristics of these AS events were systematically investigated. A total of 940, 590 and 553 DAS events were identified between low- and high-salinity environments in the livers of three species. The results of GO and KEGG enrichment analysis indicated that these DAS genes were tightly related to RNA processing, particularly RNA splicing. The common DAS genes in three species worked as RNA-binding proteins, thereby, playing crucial roles in the regulation of RNA splicing. This finding suggests that different salinity environments significantly changed the splicing patterns of these RNA processing regulators, which, in turn, may affect the AS of downstream target genes in response to salinity changes in teleosts. Taken together, this study provides preliminary evidence for the important roles of AS events in salinity adaptation in teleosts. Further functional studies are needed to elucidate the underlying mechanisms.
Conceptualization, Q.G. and Y.T.; software, H.Y. and D.L.; supervision, S.D.; data curation, Y.Z.; writing—original draft preparation, Y.T.; writing—review and editing, Q.G.; visualization, W.Y. All authors have read and agreed to the published version of the manuscript.
This study was funded by the National Key Research and Development Program of China (2019YFD0901000), the National Natural Science Foundation of China (31872575) and the Key Research and Development Program of Shandong Province (2021SFGC0701).
All animal experiments in the present study were conducted in accordance with the guidelines and approval of the respective Animal Research and Ethics Committees of Ocean University of China (Permit Number: 20141201). The field studies did not involve endangered or protected species.
Not applicable.
RNA-Seq data of turbot, tongue sole and steelhead trout can be derived from NCBI SRA database with BioProject accession number of PRJNA656408, PRJNA315357 and PRJNA797919, respectively.
The authors declare no conflict of interest.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Enlarge this image.Figure 1. Landscapes of AS events. (A) Schematic diagram of five typical AS events. The black solid boxes are used to represent the constitutive exons, while e1 and e2 indicate the alternative spliced fragments. The “I” and “E” letters represent the inclusion and exclusion isoforms, and their reads are marked by carmine and azure lines, respectively. (B) The statistics of alternative splicing events in the livers of turbot, tongue sole and steelhead trout.
Enlarge this image.Figure 2. UpSet plot showing the intersections between genes and their AS events in the livers of turbot (A), tongue sole (B) and steelhead trout (C). The upper and left bars were used to represent the number of genes and their AS events. Genes with four and five different typical AS events were marked with azure and red colors.
Enlarge this image.Figure 3. The distribution of protein-coding genes and AS events across the whole genomes of turbot (A), tongue sole (B) and steelhead trout (C). The bar plots from the outer to inner were used to represent the: (a) genome karyotype; (b) distribution of all the protein-coding genes; (c) distribution of protein-coding genes with AS events; (d) distribution of all the AS events; (e) density of AS events (AS event number/gene number); and (f) AS hotspot distribution. The distribution of both genes and AS events were calculated within non-overlap 200 kb bins.
Enlarge this image.Figure 4. Inclusion level differences of DAS events in the livers of turbot (A), tongue sole (B) and steelhead trout (C) between low- and high-salinity groups. These differences indicate the inclusion levels of low-salinity relative to high-salinity groups.
Enlarge this image.Figure 5. GO (A) and KEGG (B) enrichment analysis of identified DAS genes in the livers of turbot, tongue sole and steelhead trout under different salinity environments.
Enlarge this image.Figure 6. Common DAS genes associated with RNA processing in the livers of turbot, tongue sole and steelhead trout under different salinity environments. (A) Venn diagram showing the numbers of common DAS genes among three species. (B) GOChord plot of the relationship between the list of common DAS genes and their corresponding species classification. The full names of common DAS genes are given in Supplementary Table S8.
Enlarge this image.Figure 7. Schematic diagram of common DAS genes associated with RNA processing in the livers of turbot, tongue sole and steelhead trout under different salinity environments. Red and green arrows are used to represent positive and negative regulation function, respectively. The proteins, coded by common DAS genes, are marked by colorful modules. The full names of common DAS genes are given in Supplementary Table S8.
DAS events in the livers of turbot, tongue sole and steelhead trout under different salinity environments. The DAS events were determined using calculating the inclusion levels of low-salinity groups in relative to high-salinity groups in three species. The gene number of each AS type are listed in the brackets.
| AS Events | Turbot | Tongue Sole | Steelhead Trout |
|---|---|---|---|
| ES | 544 (469) | 404 (344) | 222 (201) |
| IR | 80 (75) | 45 (45) | 84 (75) |
| MXE | 47 (38) | 32 (26) | 40 (33) |
| A5SS | 115 (108) | 39 (38) | 84 (81) |
| A3SS | 154 (141) | 70 (65) | 123 (118) |
| Total | 940 (769) | 590 (482) | 553 (467) |
Supplementary Materials
The following supporting information can be downloaded at:
References
1. Tseng, Y.C.; Hwang, P.P. Some insights into energy metabolism for osmoregulation in fish. Comp. Biochem. Physiol. C Toxicol. Harmacol.; 2008; 148, pp. 419-429. [DOI: https://dx.doi.org/10.1016/j.cbpc.2008.04.009] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18539088]
2. Kültz, D. Physiological mechanisms used by fish to cope with salinity stress. J. Exp. Biol.; 2015; 218, Pt 12, pp. 1907-1914. [DOI: https://dx.doi.org/10.1242/jeb.118695] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26085667]
3. Tine, M.; Kuhl, H.; Gagnaire, P.A.; Louro, B.; Desmarais, E.; Martins, R.S.; Hecht, J.; Knaust, F.; Belkhir, K.; Klages, S. et al. European sea bass genome and its variation provide insights into adaptation to euryhalinity and speciation. Nat. Commun.; 2014; 5, 5770. [DOI: https://dx.doi.org/10.1038/ncomms6770] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25534655]
4. McCormick, S.D. Endocrine Control of Osmoregulation in Teleost Fish. Am. Zool.; 2001; 41, pp. 781-794. [DOI: https://dx.doi.org/10.1093/icb/41.4.781]
5. Zhang, H.; Hou, J.; Liu, H.; Zhu, H.; Xu, G.; Xu, J. Adaptive evolution of low-salinity tolerance and hypoosmotic regulation in a euryhaline teleost, Takifugu obscurus. Mar. Biol.; 2020; 167, 90. [DOI: https://dx.doi.org/10.1007/s00227-020-03705-x]
6. Wang, Z.; Gerstein, M.; Snyder, M. RNA-Seq: A revolutionary tool for transcriptomics. Nat. Rev. Genet.; 2009; 10, pp. 57-63. [DOI: https://dx.doi.org/10.1038/nrg2484]
7. Nguyen, T.V.; Jung, H.; Nguyen, T.M.; Hurwood, D.; Mather, P. Evaluation of potential candidate genes involved in salinity tolerance in striped catfish (Pangasianodon hypophthalmus) using an RNA-Seq approach. Mar. Genom.; 2016; 25, pp. 75-88. [DOI: https://dx.doi.org/10.1016/j.margen.2015.11.010]
8. Si, Y.; Wen, H.; Li, Y.; He, F.; Li, J.; Li, S.; He, H. Liver transcriptome analysis reveals extensive transcriptional plasticity during acclimation to low salinity in Cynoglossus semilaevis. BMC Genom.; 2018; 19, 464. [DOI: https://dx.doi.org/10.1186/s12864-018-4825-4]
9. Tian, Y.; Wen, H.; Qi, X.; Zhang, X.; Liu, S.; Li, B.; Sun, Y.; Li, J.; He, F.; Yang, W. et al. Characterization of Full-Length Transcriptome Sequences and Splice Variants of Lateolabrax maculatus by Single-Molecule Long-Read Sequencing and Their Involvement in Salinity Regulation. Front. Genet.; 2019; 10, 1126. [DOI: https://dx.doi.org/10.3389/fgene.2019.01126]
10. Xiong, Y.; Dong, S.; Huang, M.; Li, Y.; Wang, X.; Wang, F.; Ma, S.; Zhou, Y. Growth, osmoregulatory response, adenine nucleotide contents, and liver transcriptome analysis of steelhead trout (Oncorhynchus mykiss) under different salinity acclimation methods. Aquaculture; 2020; 520, 734937. [DOI: https://dx.doi.org/10.1016/j.aquaculture.2020.734937]
11. Liu, Z.; Ma, A.; Yuan, C.; Zhao, T.; Chang, H.; Zhang, J. Transcriptome analysis of liver lipid metabolism disorders of the turbot Scophthalmus maximus in response to low salinity stress. Aquaculture; 2021; 534, 736273. [DOI: https://dx.doi.org/10.1016/j.aquaculture.2020.736273]
12. Baralle, F.E.; Giudice, J. Alternative splicing as a regulator of development and tissue identity. Nat. Rev. Mol. Cell Biol.; 2017; 18, pp. 437-451. [DOI: https://dx.doi.org/10.1038/nrm.2017.27] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28488700]
13. Laloum, T.; Martin, G.; Duque, P. Alternative Splicing Control of Abiotic Stress Responses. Trends Plant Sci.; 2018; 23, pp. 140-150. [DOI: https://dx.doi.org/10.1016/j.tplants.2017.09.019] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29074233]
14. Pan, Q.; Shai, O.; Lee, L.J.; Frey, B.J.; Blencowe, B.J. Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing. Nat. Genet.; 2008; 40, pp. 1413-1415. [DOI: https://dx.doi.org/10.1038/ng.259] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18978789]
15. Lee, Y.; Rio, D.C. Mechanisms and Regulation of Alternative Pre-mRNA Splicing. Annu. Rev. Biochem.; 2015; 84, pp. 291-323. [DOI: https://dx.doi.org/10.1146/annurev-biochem-060614-034316] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25784052]
16. Marquez, Y.; Brown, J.W.; Simpson, C.; Barta, A.; Kalyna, M. Transcriptome survey reveals increased complexity of the alternative splicing landscape in Arabidopsis. Genome Res.; 2012; 22, pp. 1184-1195. [DOI: https://dx.doi.org/10.1101/gr.134106.111]
17. Treutlein, B.; Gokce, O.; Quake, S.R.; Sudhof, T.C. Cartography of neurexin alternative splicing mapped by single-molecule long-read mRNA sequencing. Proc. Natl. Acad. Sci. USA; 2014; 111, pp. E1291-E1299. [DOI: https://dx.doi.org/10.1073/pnas.1403244111]
18. Nilsen, T.W.; Graveley, B.R. Expansion of the eukaryotic proteome by alternative splicing. Nature; 2010; 463, pp. 457-463. [DOI: https://dx.doi.org/10.1038/nature08909] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20110989]
19. Li, C.Y.; Chu, J.Y.; Yu, J.K.; Huang, X.Q.; Liu, X.J.; Li, S.; Che, Y.C.; Xie, J.Y. Regulation of alternative splicing of Bcl-x by IL-6, GM-CSF and TPA. Cell Res.; 2004; 14, pp. 473-479. [DOI: https://dx.doi.org/10.1038/sj.cr.7290250]
20. Ni, J.Z.; Grate, L.; Donohue, J.P.; Preston, C.; Nobida, N.; O’Brien, G.; Shiue, L.; Clark, T.A.; Blume, J.E.; Ares, M., Jr. Ultraconserved elements are associated with homeostatic control of splicing regulators by alternative splicing and nonsense-mediated decay. Genes Dev.; 2007; 21, pp. 708-718. [DOI: https://dx.doi.org/10.1101/gad.1525507]
21. Kalyna, M.; Simpson, C.G.; Syed, N.H.; Lewandowska, D.; Marquez, Y.; Kusenda, B.; Marshall, J.; Fuller, J.; Cardle, L.; McNicol, J. et al. Alternative splicing and nonsense-mediated decay modulate expression of important regulatory genes in Arabidopsis. Nucleic Acids Res.; 2012; 40, pp. 2454-2469. [DOI: https://dx.doi.org/10.1093/nar/gkr932]
22. Hamid, F.M.; Makeyev, E.V. Emerging functions of alternative splicing coupled with nonsense-mediated decay. Biochem. Soc. Trans.; 2014; 42, pp. 1168-1173. [DOI: https://dx.doi.org/10.1042/BST20140066]
23. Tan, S.; Wang, W.; Tian, C.; Niu, D.; Zhou, T.; Jin, Y.; Yang, Y.; Gao, D.; Dunham, R.; Liu, Z. Heat stress induced alternative splicing in catfish as determined by transcriptome analysis. Comp. Biochem. Physiol. Part D Genom. Proteomics.; 2019; 29, pp. 166-172. [DOI: https://dx.doi.org/10.1016/j.cbd.2018.11.008]
24. Sun, J.; Liu, Z.; Quan, J.; Li, L.; Zhao, G.; Lu, J. RNA-seq Analysis Reveals Alternative Splicing Under Heat Stress in Rainbow Trout (Oncorhynchus mykiss). Mar. Biotechnol.; 2021; [DOI: https://dx.doi.org/10.1007/s10126-021-10082-z]
25. Li, B.J.; Zhu, Z.X.; Qin, H.; Meng, Z.N.; Lin, H.R.; Xia, J.H. Genome-Wide Characterization of Alternative Splicing Events and Their Responses to Cold Stress in Tilapia. Front. Genet.; 2020; 11, 244. [DOI: https://dx.doi.org/10.3389/fgene.2020.00244]
26. Li, H.L.; Lin, H.R.; Xia, J.H. Differential Gene Expression Profiles and Alternative Isoform Regulations in Gill of Nile Tilapia in Response to Acute Hypoxia. Mar. Biotechnol.; 2017; 19, pp. 551-562. [DOI: https://dx.doi.org/10.1007/s10126-017-9774-4]
27. Xia, J.H.; Li, H.L.; Li, B.J.; Gu, X.H.; Lin, H.R. Acute hypoxia stress induced abundant differential expression genes and alternative splicing events in heart of tilapia. Gene; 2018; 639, pp. 52-61. [DOI: https://dx.doi.org/10.1016/j.gene.2017.10.002]
28. Tian, Y.; Wen, H.; Qi, X.; Zhang, X.; Sun, Y.; Li, J.; He, F.; Zhang, M.; Zhang, K.; Yang, W. et al. Alternative splicing (AS) mechanism plays important roles in response to different salinity environments in spotted sea bass. Int. J. Biol. Macromol.; 2020; 155, pp. 50-60. [DOI: https://dx.doi.org/10.1016/j.ijbiomac.2020.03.178]
29. Gonzalez, R.J. The physiology of hyper-salinity tolerance in teleost fish: A review. J. Comp. Physiol. B; 2012; 182, pp. 321-329. [DOI: https://dx.doi.org/10.1007/s00360-011-0624-9]
30. Fiol, D.F.; Kültz, D. Osmotic stress sensing and signaling in fishes. FEBS J.; 2007; 274, pp. 5790-5798. [DOI: https://dx.doi.org/10.1111/j.1742-4658.2007.06099.x]
31. Wang, G.; Xu, K.; Tian, X.; Dong, S.; Fang, Z. Changes in plasma osmolality, cortisol and amino acid levels of tongue sole (Cynoglossus semilaevis) at different salinities. J. Ocean Univ. China; 2015; 14, pp. 881-887. [DOI: https://dx.doi.org/10.1007/s11802-015-2598-9]
32. Cui, W.; Ma, A.; Huang, Z.; Wang, X.; Liu, Z.; Xia, D.; Yang, S.; Zhao, T. Comparative transcriptomic analysis reveals mechanisms of divergence in osmotic regulation of the turbot Scophthalmus maximus. Fish Physiol. Biochem.; 2020; 46, pp. 1519-1536. [DOI: https://dx.doi.org/10.1007/s10695-020-00808-6]
33. Lee, M.; Lee, B.; Kim, K.; Yoon, M.; Lee, J.W. Differential effects of two seawater transfer regimes on the hypoosmoregulatory adaptation, hormonal response, feed efficiency, and growth performance of juvenile steelhead trout. Aquac. Rep.; 2022; 22, 101004. [DOI: https://dx.doi.org/10.1016/j.aqrep.2021.101004]
34. Aluru, N.; Vijayan, M.M. Hepatic transcriptome response to glucocorticoid receptor activation in rainbow trout. Physiol. Genom.; 2007; 31, pp. 483-491. [DOI: https://dx.doi.org/10.1152/physiolgenomics.00118.2007]
35. Chang, C.H.; Wang, Y.C.; Lee, T.H. Hypothermal stress-induced salinity-dependent oxidative stress and apoptosis in the livers of euryhaline milkfish, Chanos chanos. Aquaculture; 2021; 534, 736280. [DOI: https://dx.doi.org/10.1016/j.aquaculture.2020.736280]
36. Graveley, B.R. Alternative splicing: Increasing diversity in the proteomic world. Trends Genet.; 2001; 17, pp. 100-107. [DOI: https://dx.doi.org/10.1016/S0168-9525(00)02176-4]
37. Stamm, S.; Ben-Ari, S.; Rafalska, I.; Tang, Y.; Zhang, Z.; Toiber, D.; Thanaraj, T.A.; Soreq, H. Function of alternative splicing. Gene; 2005; 344, pp. 1-30. [DOI: https://dx.doi.org/10.1016/j.gene.2004.10.022]
38. Heine-Suñer, D.; Díaz-Guillén, M.A.; Lange, A.J.; Rodríguez de Córdoba, S. Sequence and structure of the human 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase heart isoform gene (PFKFB2). Eur. J. Biochem.; 1998; 254, pp. 103-110. [DOI: https://dx.doi.org/10.1046/j.1432-1327.1998.2540103.x]
39. Tang, N.; Zhang, J.; Fu, X.; Xie, W.; Qiu, Y. PP2Acα inhibits PFKFB2-induced glycolysis to promote termination of liver regeneration. Biochem. Biophys. Res. Commun.; 2020; 526, pp. 1-7. [DOI: https://dx.doi.org/10.1016/j.bbrc.2020.03.002]
40. Jankowski, M.P.; Cornuet, P.K.; McIlwrath, S.; Koerber, H.R.; Albers, K.M. SRY-box containing gene 11 (Sox11) transcription factor is required for neuron survival and neurite growth. Neuroscience; 2006; 143, pp. 501-514. [DOI: https://dx.doi.org/10.1016/j.neuroscience.2006.09.010]
41. Yin, C. Molecular mechanisms of Sox transcription factors during the development of liver, bile duct, and pancreas. Semin. Cell Dev. Biol.; 2017; 63, pp. 68-78. [DOI: https://dx.doi.org/10.1016/j.semcdb.2016.08.015]
42. Ong, Y.S.; Tang, B.L.; Loo, L.S.; Hong, W. p125A exists as part of the mammalian Sec13/Sec31 COPII subcomplex to facilitate ER-Golgi transport. J. Cell Biol.; 2010; 190, pp. 331-345. [DOI: https://dx.doi.org/10.1083/jcb.201003005]
43. Schie, I.M.; Castrop, H. Regulation of NKCC2 splicing and phosphorylation. Curr. Opin. Nephrol. Hypertens.; 2015; 24, pp. 457-462. [DOI: https://dx.doi.org/10.1097/MNH.0000000000000150]
44. Wang, X.; Du, L.; Wei, H.; Zhang, A.; Yang, K.; Zhou, H. Identification of two Stat3 variants lacking a transactivation domain in grass carp: New insights into alternative splicing in the modification of teleost Stat3 signaling. Fish Shellfish Immunol.; 2018; 77, pp. 13-21. [DOI: https://dx.doi.org/10.1016/j.fsi.2018.03.022]
45. Liu, L.; Zhou, Q.; Lin, C.; He, L.; Wei, L. Integrative analyses of gene expression and alternative splicing to gain insights into the effects of copper on hepatic lipid metabolism in swamp eel (Monopterus albus). Aquaculture; 2022; 546, 737367. [DOI: https://dx.doi.org/10.1016/j.aquaculture.2021.737367]
46. Amara, U.; Flierl, M.A.; Rittirsch, D.; Klos, A.; Chen, H.; Acker, B.; Brückner, U.B.; Nilsson, B.; Gebhard, F.; Lambris, J.D. et al. Molecular intercommunication between the complement and coagulation systems. J. Immunol.; 2010; 185, pp. 5628-5636. [DOI: https://dx.doi.org/10.4049/jimmunol.0903678]
47. Oikonomopoulou, K.; Ricklin, D.; Ward, P.A.; Lambris, J.D. Interactions between coagulation and complement--their role in inflammation. Semin. Immunopathol.; 2012; 34, pp. 151-165. [DOI: https://dx.doi.org/10.1007/s00281-011-0280-x]
48. Dzik, S. Complement and Coagulation: Cross Talk Through Time. Transfus. Med. Rev.; 2019; 33, pp. 199-206. [DOI: https://dx.doi.org/10.1016/j.tmrv.2019.08.004]
49. Noris, M.; Remuzzi, G. Overview of complement activation and regulation. Semin. Nephrol.; 2013; 33, pp. 479-492. [DOI: https://dx.doi.org/10.1016/j.semnephrol.2013.08.001]
50. Bajic, G.; Degn, S.E.; Thiel, S.; Andersen, G.R. Complement activation, regulation, and molecular basis for complement-related diseases. EMBO J.; 2015; 34, pp. 2735-2757. [DOI: https://dx.doi.org/10.15252/embj.201591881]
51. Jiang, C.; Zhang, J.; Yao, J.; Liu, S.; Li, Y.; Song, L.; Li, C.; Wang, X.; Liu, Z. Complement regulatory protein genes in channel catfish and their involvement in disease defense response. Dev. Comp. Immunol.; 2015; 53, pp. 33-41. [DOI: https://dx.doi.org/10.1016/j.dci.2015.06.002]
52. Keihanian, F.; Saeidinia, A.; Bagheri, R.K.; Johnston, T.P.; Sahebkar, A. Curcumin, hemostasis, thrombosis, and coagulation. J. Cell. Physiol.; 2018; 233, pp. 4497-4511. [DOI: https://dx.doi.org/10.1002/jcp.26249]
53. Miriami, E.; Margalit, H.; Sperling, R. Conserved sequence elements associated with exon skipping. Nucleic Acids Res.; 2003; 31, pp. 1974-1983. [DOI: https://dx.doi.org/10.1093/nar/gkg279]
54. Grau-Bové, X.; Ruiz-Trillo, I.; Irimia, M. Origin of exon skipping-rich transcriptomes in animals driven by evolution of gene architecture. Genome Biol.; 2018; 19, 135. [DOI: https://dx.doi.org/10.1186/s13059-018-1499-9]
55. Ge, Y.; Porse, B.T. The functional consequences of intron retention: Alternative splicing coupled to NMD as a regulator of gene expression. Bioessays; 2014; 36, pp. 236-243. [DOI: https://dx.doi.org/10.1002/bies.201300156]
56. Grabski, D.F.; Broseus, L.; Kumari, B.; Rekosh, D.; Hammarskjold, M.L.; Ritchie, W. Intron retention and its impact on gene expression and protein diversity: A review and a practical guide. Wiley Interdiscip Rev. RNA; 2021; 12, e1631. [DOI: https://dx.doi.org/10.1002/wrna.1631]
57. Koncz, C.; Dejong, F.; Villacorta, N.; Szakonyi, D.; Koncz, Z. The spliceosome-activating complex: Molecular mechanisms underlying the function of a pleiotropic regulator. Front. Plant Sci.; 2012; 3, 9. [DOI: https://dx.doi.org/10.3389/fpls.2012.00009]
58. Wilkinson, M.E.; Charenton, C.; Nagai, K. RNA Splicing by the Spliceosome. Annu. Rev. Biochem.; 2020; 89, pp. 359-388. [DOI: https://dx.doi.org/10.1146/annurev-biochem-091719-064225]
59. Martinez-Contreras, R.; Cloutier, P.; Shkreta, L.; Fisette, J.F.; Revil, T.; Chabot, B. hnRNP proteins and splicing control. Adv. Exp. Med. Biol.; 2007; 623, pp. 123-147. [DOI: https://dx.doi.org/10.1007/978-0-387-77374-2_8]
60. Morton, M.; AlTamimi, N.; Butt, H.; Reddy, A.S.N.; Mahfouz, M. Serine/Arginine-rich protein family of splicing regulators: New approaches to study splice isoform functions. Plant Sci.; 2019; 283, pp. 127-134. [DOI: https://dx.doi.org/10.1016/j.plantsci.2019.02.017]
61. Tan, S.; Wang, W.; Tian, C.; Niu, D.; Zhou, T.; Yang, Y.; Gao, D.; Liu, Z. Post-transcriptional regulation through alternative splicing after infection with Flavobacterium columnare in channel catfish (Ictalurus punctatus). Fish Shellfish Immunol.; 2019; 91, pp. 188-193. [DOI: https://dx.doi.org/10.1016/j.fsi.2019.05.008]
62. Ladd, A.N.; Charlet, N.; Cooper, T.A. The CELF family of RNA binding proteins is implicated in cell-specific and developmentally regulated alternative splicing. Mol. Cell. Biol.; 2001; 21, pp. 1285-1296. [DOI: https://dx.doi.org/10.1128/MCB.21.4.1285-1296.2001]
63. Gromak, N.; Matlin, A.J.; Cooper, T.A.; Smith, C.W. Antagonistic regulation of alpha-actinin alternative splicing by CELF proteins and polypyrimidine tract binding protein. RNA; 2003; 9, pp. 443-456. [DOI: https://dx.doi.org/10.1261/rna.2191903]
64. Hase, M.E.; Yalamanchili, P.; Visa, N. The Drosophila heterogeneous nuclear ribonucleoprotein M protein, HRP59, regulates alternative splicing and controls the production of its own mRNA. J. Biol. Chem.; 2006; 281, pp. 39135-39141. [DOI: https://dx.doi.org/10.1074/jbc.M604235200]
65. Markovtsov, V.; Nikolic, J.M.; Goldman, J.A.; Turck, C.W.; Chou, M.Y.; Black, D.L. Cooperative assembly of an hnRNP complex induced by a tissue-specific homolog of polypyrimidine tract binding protein. Mol. Cell. Biol.; 2000; 20, pp. 7463-7479. [DOI: https://dx.doi.org/10.1128/MCB.20.20.7463-7479.2000]
66. De Maio, A.; Yalamanchili, H.K.; Adamski, C.J.; Gennarino, V.A.; Liu, Z.; Qin, J.; Jung, S.Y.; Richman, R.; Orr, H.; Zoghbi, H.Y. RBM17 Interacts with U2SURP and CHERP to Regulate Expression and Splicing of RNA-Processing Proteins. Cell Rep.; 2018; 25, pp. 726-736.e7. [DOI: https://dx.doi.org/10.1016/j.celrep.2018.09.041]
67. Xu, Y.; Nijhuis, A.; Keun, H.C. RNA-binding motif protein 39 (RBM39): An emerging cancer target. Br. J. Pharmacol.; 2020; [DOI: https://dx.doi.org/10.1111/bph.15331]
68. Zheng, X.; Peng, Q.; Wang, L.; Zhang, X.; Huang, L.; Wang, J.; Qin, Z. Serine/arginine-rich splicing factors: The bridge linking alternative splicing and cancer. Int. J. Biol. Sci.; 2020; 16, pp. 2442-2453. [DOI: https://dx.doi.org/10.7150/ijbs.46751]
69. Iakova, P.; Wang, G.L.; Timchenko, L.; Michalak, M.; Pereira-Smith, O.M.; Smith, J.R.; Timchenko, N.A. Competition of CUGBP1 and calreticulin for the regulation of p21 translation determines cell fate. EMBO J.; 2004; 23, pp. 406-417. [DOI: https://dx.doi.org/10.1038/sj.emboj.7600052]
70. Zhou, W.; Chung, Y.J.; Parrilla Castellar, E.R.; Zheng, Y.; Chung, H.J.; Bandle, R.; Liu, J.; Tessarollo, L.; Batchelor, E.; Aplan, P.D. et al. Far Upstream Element Binding Protein Plays a Crucial Role in Embryonic Development, Hematopoiesis, and Stabilizing Myc Expression Levels. Am. J. Pathol.; 2016; 186, pp. 701-715. [DOI: https://dx.doi.org/10.1016/j.ajpath.2015.10.028]
71. Zhang, Z.; Yao, Z.; Wang, L.; Ding, H.; Shao, J.; Chen, A.; Zhang, F.; Zheng, S. Activation of ferritinophagy is required for the RNA-binding protein ELAVL1/HuR to regulate ferroptosis in hepatic stellate cells. Autophagy; 2018; 14, pp. 2083-2103. [DOI: https://dx.doi.org/10.1080/15548627.2018.1503146] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30081711]
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Gao, Qinfeng 2 ; Dong, Shuanglin 2 ; Zhou, Yangen 2 ; Han, Yu 1 ; Liu, Dazhi 1 ; Yang, Wenzhao 1 1 Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, Qingdao 266003, China;
2 Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, Qingdao 266003, China;