Genome-Wide Analysis of Alternative Splicing Provides Insights into Stress Adaptation of the Pacific Oyster

Alternative splicing plays a nonredundant role in greater amberjack (Seriola dumerili) in acclimation to ambient salinity fluctuations

Alternative splicing (AS) is an important post-transcriptional mechanism for adaptation of fish to environmental stress. Here, we performed a genome-wide investigation to AS dynamics in greater amberjack (Seriola dumerili), an economical marine teleost, in response to hypo- (10 ppt) and hyper-salinity (40 ppt) stresses. Totally, 2267-2611 differentially spliced events were identified in gills and kidney upon the exposure to undesired salinity regimes. In gills, genes involved in energy metabolism, stimulus response and epithelial cell differentiation were differentially spliced in response to salinity variation, while sodium ion transport and cellular amide metabolism were enhanced in kidney to combat the adverse impacts of salinity changes. Most of these differentially spliced genes were not differentially expressed, and AS was found to regulate different biological processes from differential gene expression, indicative of the functionally nonredundant role of AS in modulating salinity acclimation in greater amberjack. Together, our study highlights the important contribution of post-transcriptional mechanisms to the adaptation of fish to ambient salinity fluctuations and provides theoretical guidance for the conservation of marine fishery resources against increasingly environmental challenges.

Alternative Splicing Reveals Acute Stress Response of Litopenaeus vannamei at High Alkalinity

Alkalinity is regarded as one of the primary stressors for aquatic animals in saline-alkaline water. Alternative splicing (AS) can significantly increase the diversity of transcripts and play key roles in stress response; however, the studies on AS under alkalinity stress of crustaceans are still limited. In the present study, we devoted ourselves to the study of AS under acute alkalinity stress at control (50 mg/L) and treatment groups (350 mg/L) by RNA-seq in pacific white shrimp (Litopenaeus vannamei). We identified a total of 10,556 AS events from 4865 genes and 619 differential AS (DAS) events from 519 DAS genes in pacific white shrimp. Functional annotation showed that the DAS genes primarily involved in spliceosome. Five splicing factors (SFs), U2AF1, PUF60, CHERP, SR140 and SRSF2 were significantly up-regulated and promoted AS. Furthermore, alkalinity activated the Leukocyte transendothelial migration, mTOR signaling pathway and AMPK signaling pathway, which regulated MAPK1, EIF3B and IGFP-RP1 associated with these pathways. We also studied three SFs (HSFP1, SRSF2 and NHE-RF1), which underwent AS to form different transcript isoforms. The above results demonstrated that AS was a regulatory mechanism in pacific white shrimp in response to acute alkalinity stress. SFs played vital roles in AS of pacific white shrimp, such as HSFP1, SRSF2 and NHE-RF1. DAS genes were significantly modified in immunity of pacific white shrimp to cope with alkalinity stress. This is the first study on the response of AS to acute alkalinity stress, which provided scientific basis for AS mechanism of crustaceans response to alkalinity stress.

Identification, diversity, and evolution analysis of thioester-containing protein family in Pacific oyster (Crassostrea gigas) and immune response to biotic and abiotic stresses

Thioester-containing proteins (TEPs) play a vital role in the innate immune response to biotic and abiotic stresses. In this study, the TEPs in C. gigas were identified, and their gene structure, phylogenetic relationships, collinearity relationships, expression profiles, sequence diversity, and alternative splicing were analyzed. Eight Tep genes were identified in C. gigas genome. Functional analysis and evolutionary relationships indicated a high level of homology to other mollusks TEPs. The transcriptome quantitative analysis results showed that the Tep genes in C. gigas respond to heat stress and Vibrio stress. Alternative splicing analysis revealed four Tep genes (designated A2M_1, CD109_3, CD109_5, complement C3) encode multiple alternative splice variants. Analysis of gene structure and multiple alignments revealed that seven CD109_5 variants are produced through the alternative splicing of the 19th exon, which encodes the highly variable central region. Sequence diversity analysis revealed thirteen missense variants within the 19th exon region of these seven CD109_5 alternative splice variants. Furthermore, the differential alternative splicing analysis showed significant induction of CD109_5, A2M_1 and A2M_2 variants after infection with V. parahaemolyticus. This study explores the Tep genes of C. gigas, providing insights into the molecular mechanisms underlying the involvement of C. gigas TEPs in innate immunity.

Ancestors' Gift: Parental Early Exposure to the Environmentally Realistic Pesticide Mixture Drives Offspring Phenotype in a Larger Extent Than Direct Exposure in the Pacific Oyster, Crassostrea gigas

Marine organisms are threatened by the presence of pesticides in coastal waters. Among them, the Pacific oyster is one of the most studied invertebrates in marine ecotoxicology where numerous studies highlighted the multiscale impacts of pesticides. In the past few years, a growing body of literature has reported the epigenetic outcomes of xenobiotics. Because DNA methylation is an epigenetic mark implicated in organism development and is meiotically heritable, it raises the question of the multigenerational implications of xenobiotic-induced epigenetic alterations. Therefore, we performed a multigenerational exposure to an environmentally relevant mixture of 18 pesticides (nominal sum concentration: 2.85 μg·L-1) during embryo-larval stages (0-48 hpf) of a second generation (F1) for which parents where already exposed or not in F0. Gene expression, DNA methylation, and physiological end points were assessed throughout the life cycle of individuals. Overall, the multigenerational effect has a greater influence on the phenotype than the exposure itself. Thus, multigenerational phenotypic effects were observed: individuals descending from exposed parents exhibited lower epinephrine-induced metamorphosis and field survival rates. At the molecular level, RNA-seq and Methyl-seq data analyses performed in gastrula embryos and metamorphosis-competent pediveliger (MCP) larvae revealed a clear F0 treatment-dependent discrimination. Some genes implicated into shell secretion and immunity exhibited F1:F0 treatment interaction patterns (e.g., Calm and Myd88). Those results suggest that low chronic environmental pesticide contamination can alter organisms beyond the individual scale level and have long-term adaptive implications.

The Responses of Alternative Splicing during Heat Stress in the Pacific White Shrimp Litopenaeus vannamei

Heat tolerance is increasingly becoming a crucial trait for aquaculture species in the face of rapidly changing climate conditions. Alternative splicing (AS) is a vital mechanism within cells that modulates gene abundance and functional diversity, enabling organisms to effectively respond to diverse stressful conditions, including thermal stress. However, it is still uncertain whether AS contributes to heat tolerance in shrimp. In this study, we conducted an extensive transcriptome analysis on the Pacific white shrimp, Litopenaeus vannamei, revealing a total of 1267, 987, and 130 differential AS events (DAS) in the gill, hepatopancreas, and muscle, respectively, following exposure to heat stress. Among all of the DAS events, exon skipping (ES) was the predominant form of splicing modification observed. Interestingly, a minor portion of DAS genes exhibited overlap across the three tissues, implying that heat stress exerts unique effects on various tissue types. Moreover, the functional enrichment analysis demonstrated that commonly identified DAS genes were primarily associated with the "spliceosome" pathway, indicating that the AS of splicing-related genes played a crucial role in the response to heat stress. Our findings also revealed that heat stress tended to induce longer mRNA isoforms through differential alternative 3' splice site (A3SS) events. Notably, A3SS events exhibited the highest proportion of maintained open reading frames (ORFs) under heat stress. Interestingly, we observed a limited overlap between the genes exhibiting DAS and those showing differential gene expression (DEG), indicating that AS may function as a distinct regulatory mechanism independent of transcriptional regulation in response to heat stress. This is the first comprehensive study on AS in crustacea species under heat stress, which broadens our understanding of the regulatory mechanisms governing the crustaceans' response to environmental stress, providing valuable insights for the aquaculture breeding of shrimp and other aquatic animals.

Identification of sex-specific splicing via comparative transcriptome analysis of gonads from sea cucumber Apostichopus japonicus

Alternative splicing (AS) is an essential post-transcriptional regulation mechanism for sex differentiation and gonadal development, which has rarely been reported in marine invertebrates. Sea cucumber (Apostichopus japonicus) is an economically important marine benthic echinoderm with a potential XX/XY sex determination mechanism, whose molecular mechanism in the gonadal differentiation has not been clearly understood. In this study, we analyzed available RNA-seq datasets of male and female gonads to explore if AS mechanism exerts an essential function in sex differentiation and gonadal development of A. japonicus. In our results, a total of 20,338 AS events from 7219 alternatively spliced genes, and 189 sexually differential alternative splicing (DAS) events from 156 genes were identified in gonadal transcriptome of sea cucumber. Gene Ontology analysis indicated that these DAS genes were significantly enriched in spermatogenesis-related GO terms. Maximal Clique Centrality (MCC) was then applied for protein-protein interaction (PPI) analysis to search for protein interactions and hub DAS gene. Among all DAS genes, we identified 10 DAS genes closely related to spermatogenesis and (or) sperm motility and a hub gene dnah1. Thus, this study revealed that alternative isoforms were generated from certain genes in female and male gonads through alternative splicing, which may provide direct evidence that alternative splicing mechanisms participate in female and male gonads. These results suggested a novel perspective for explaining the molecular mechanisms underlying gonadal differentiation between male and female sea cucumbers.

RNA-seq Analysis Reveals Alternative Splicing Under Heat Stress in Rainbow Trout (Oncorhynchus mykiss)

Rainbow trout (Oncorhynchus mykiss) is one of the most economically important cold-water farmed species in the world, and transcriptomic studies in response to heat stress have been conducted and will be studied in depth. Alternative splicing (AS), a post-transcriptional regulatory process that regulates gene expression and increases proteomic diversity, is still poorly understood in rainbow trout under heat stress. In the present study, 18,623 alternative splicing events were identified from 9936 genes using RNA transcriptome sequencing technology (RNA-Seq) and genomic information. A total of 2731 differential alternative splicing (DAS) events were found among 2179 differentially expressed genes (DEGs). Gene ontology analysis revealed that the DEGs were mainly enriched in cellular metabolic process, cell part, and organic cyclic compound binding under heat stress. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis displayed that the DEGs were enriched for 39 pathways, and some key pathways, such as lysine degradation, are involved in the regulation of heat stress in liver tissues of rainbow trout. The results were validated by qRT-PCR, confirming reliability of our bioinformatics analysis.

Population Genomics, Transcriptional Response to Heat Shock, and Gut Microbiota of the Hong Kong Oyster Magallana hongkongensis

The Hong Kong oyster Magallana hongkongensis, previously known as Crassostrea hongkongensis, is a true oyster species native to the estuarine-coast of the Pearl River Delta in southern China. The species—with scientific, ecological, cultural, and nutritional importance—has been farmed for hundreds of years. However, there is only limited information on its genetics, stress adaptation mechanisms, and gut microbiota, restricting the sustainable production and use of oyster resources. Here, we present population structure analysis on M. hongkongensis oysters collected from Deep Bay and Lantau Island in Hong Kong, as well as transcriptome analysis on heat shock responses and the gut microbiota profile of M. hongkongensis oysters collected from Deep Bay. Single nucleotide polymorphisms (SNPs), including those on the homeobox genes and heat shock protein genes, were revealed by the whole genome resequencing. Transcriptomes of oysters incubated at 25 °C and 32 °C for 24 h were sequenced which revealed the heat-induced regulation of heat shock protein pathway genes. Furthermore, the gut microbe community was detected by 16S rRNA sequencing which identified Cyanobacteria, Proteobacteria and Spirochaetes as the most abundant phyla. This study reveals the molecular basis for the adaptation of the oyster M. hongkongensis to environmental conditions.

Nanopore long-read RNAseq reveals regulatory mechanisms of thermally variable reef environments promoting heat tolerance of scleractinian coral Pocillopora damicornis

Some scleractinian corals exhibit high thermal adaptability to climate changes, although the mechanism of their adaptation is unclear. This study investigated the adaptability of scleractinian coral Pocillopora damicornis to thermally variable reef environments by applying a nanopore-based RNA sequencing method to characterize different transcription responses that promote heat tolerance of P. damicornis. We identified 1414 novel genes and optimized 6256 mis-annotated loci. Based on full-length transcriptome data, we identified complex alternative polyadenylation and alternative splicing events, which can improve our understanding of the genome annotation and gene structures of P. damicornis. Furthermore, we constructed differentially expressed lncRNA-mRNA co-expression networks, which may play a crucial role in the P. damicornis thermal adaptive response. KEGG function enrichment analysis revealed that P. damicornis from the high-temperature pool had a lower metabolic rate than that from the low-temperature pool. We hypothesize that metabolic readjustment, in the form of a lower metabolic rate, positively correlated with increased heat tolerance in P. damicornis in thermally variable reef environments. Our study provides novel insights into lncRNAs that promote thermally tolerance of scleractinian corals in the thermally variable reef environment, suggesting potential mechanisms for their adaptation to global warming in the future.

Alternative splicing is highly variable among Daphnia pulex lineages in response to acute copper exposure

BACKGROUND

Despite being one of the primary mechanisms of gene expression regulation in eukaryotes, alternative splicing is often overlooked in ecotoxicogenomic studies. The process of alternative splicing facilitates the production of multiple mRNA isoforms from a single gene thereby greatly increasing the diversity of the transcriptome and proteome. This process can be important in enabling the organism to cope with stressful conditions. Accurate identification of splice sites using RNA sequencing requires alignment to independent exonic positions within the genome, presenting bioinformatic challenges, particularly when using short read data. Although technological advances allow for the detection of splicing patterns on a genome-wide scale, very little is known about the extent of intraspecies variation in splicing patterns, particularly in response to environmental stressors. In this study, we used RNA-sequencing to study the molecular responses to acute copper exposure in three lineages of Daphnia pulex by focusing on the contribution of alternative splicing in addition to gene expression responses.

RESULTS

By comparing the overall gene expression and splicing patterns among all 15 copper-exposed samples and 6 controls, we identified 588 differentially expressed (DE) genes and 16 differentially spliced (DS) genes. Most of the DS genes (13) were not found to be DE, suggesting unique transcriptional regulation in response to copper that went unnoticed with conventional DE analysis. To understand the influence of genetic background on gene expression and alternative splicing responses to Cu, each of the three lineages was analyzed separately. In contrast to the overall analysis, each lineage had a higher proportion of unique DS genes than DE genes suggesting that genetic background has a larger influence on DS than on DE. Gene Ontology analysis revealed that some pathways involved in stress response were jointly regulated by DS and DE genes while others were regulated by only transcription or only splicing.

CONCLUSIONS

Our findings suggest an important role for alternative splicing in shaping transcriptome diversity in response to metal exposure in Daphnia, highlighting the importance of integrating splicing analyses with gene expression surveys to characterize molecular pathways in evolutionary and environmental studies.

Evidence from ileum and liver transcriptomes of resistance to high-salt and water-deprivation conditions in camel

Camels have evolved various resistance characteristics adaptive to their desert habitats. In the present study, we used high-throughput sequencing to investigate stress-induced alternative splicing events as well as different genes involved in resistance to water deprivation and salt absorption in the ileum and liver in Camelus bactrianus. Through association analyses of mRNA, miRNA and lncRNA, we sought to explicate how camels respond to high salt and water scarcity conditions. There were two modes by which genes driven by alternative splicing were enriched to molecular functions, invoking of which was potentially fixed by organ and stress types. With qRT-PCR detection, the differentially expressed MUC6, AQP5, LOC105076960, PKP4, CDH11, TENM1, SDS, LOC105061856, PLIN2 and UPP2 were screened as functionally important genes, along with miR-29b, miR-484, miR-362-5p, miR-96, miR-195, miR-128 and miR-148a. These genes contributed to cellular stress resistance, for instance by reducing water loss, inhibiting excessive import of sodium, improving protective barriers and sodium ion homeostasis, and maintaining uridine content. The underlying competing endogenous RNAs referred to LNC001664, let-7e and LOC105076960 mRNA in ileum, and LNC001438, LNC003417, LNC001770, miR-199c and TENM1 mRNA in liver. Besides competent interpretation to resistance, there may be inspirations for curing human diseases triggered by high-salt intake.

Genome-Wide Characterization of Alternative Splicing Events and Their Responses to Cold Stress in Tilapia

Alternative splicing (AS) is an important post-transcriptional regulatory mechanism for cells to generate transcript variability and proteome diversity. No systematic investigation of AS events among different tissues in response to stressors is available for tilapia currently. In this study, AS among different tissues was identified and the cold stress-related AS events were explored in a Nile tilapia (Oreochromis niloticus) line based on 42 RNA-seq datasets using a bioinformatics pipeline. 14,796 (82.76%; SD = 2,840) of the expression genes showed AS events. The two most abundant AS types were alternative transcription start site (TSS) and terminal site (TTS) in tilapia. Testis, brain and kidney possess the most abundant AS gene number, while the blood, muscle and liver possess the least number in each tissue. Furthermore, 208 differentially alternative splicing (DAS) genes in heart and 483 DAS in brain in response to cold stress. The number of AS types for alternative exon end, exon skipping and retention of single intron increased significantly under cold stress. GO enrichment and pathway overrepresentation analysis indicated that many DAS genes, e.g., genes in circadian clock pathway, may influence expression of downstream genes under cold stress. Our study revealed that AS exists extensively in tilapia and plays an important role in cold adaption.

Characterization of Full-Length Transcriptome Sequences and Splice Variants of Lateolabrax maculatus by Single-Molecule Long-Read Sequencing and Their Involvement in Salinity Regulation

Transcriptome complexity plays crucial roles in regulating the biological functions of eukaryotes. Except for functional genes, alternative splicing and fusion transcripts produce a vast expansion of transcriptome diversity. In this study, we applied PacBio single-molecule long-read sequencing technology to unveil the whole transcriptome landscape of Lateolabrax maculatus. We obtained 28,809 high-quality non-redundant transcripts, including 18,280 novel isoforms covering 8,961 annotated gene loci within the current reference genome and 3,172 novel isoforms. A total of 10,249 AS events were detected, and intron retention was the predominant AS event. In addition, 1,359 alternative polyadenylation events, 3,112 lncRNAs, 29,609 SSRs, 365 fusion transcripts, and 1,194 transcription factors were identified in this study. Furthermore, we performed RNA-Seq analysis combined with Iso-Seq results to investigate salinity regulation mechanism at the transcripts level. A total of 518 transcripts were differentially expressed, which were further divided into 8 functional groups. Notably, transcripts from the same genes exhibited similar or opposite expression patterns. Our study provides a comprehensive view of the transcriptome complexity in L. maculatus, which significantly improves current gene models. Moreover, the diversity of the expression patterns of transcripts may enhance the understanding of salinity regulatory mechanism in L. maculatus and other euryhaline teleosts.

Transcriptome and Gene Coexpression Network Analyses of Two Wild Populations Provides Insight into the High-Salinity Adaptation Mechanisms of Crassostrea ariakensis

Crassostrea ariakensis naturally distributes in the intertidal and estuary region with relative low salinity ranging from 10 to 25‰. To understand the adaptive capacity of oysters to salinity stress, we conducted transcriptome analysis to investigate the metabolic pathways of salinity stress effectors in oysters from two different geographical sites, namely at salinities of 16, 23, and 30‰. We completed transcriptome sequencing of 18 samples and a total of 52,392 unigenes were obtained after assembly. Differentially expressed gene (DEG) analysis and weighted gene correlation network analysis (WGCNA) were performed using RNA-Seq transcriptomic data from eye-spot larvae at different salinities and from different populations. The results showed that at moderately high salinities (23 and 30‰), genes related to osmotic agents, oxidation-reduction processes, and related regulatory networks of complex transcriptional regulation and signal transduction pathways dominated to counteract the salinity stress. Moreover, there were adaptive differences in salinity response mechanisms, especially at high salinity, in oyster larvae from different populations. These results provide a framework for understanding the interactions of multiple pathways at the system level and for elucidating the complex cellular processes involved in responding to osmotic stress and maintaining growth. Furthermore, the results facilitate further research into the biological processes underlying physiological adaptations to hypertonic stress in marine invertebrates and provide a molecular basis for our subsequent search for high salinity-tolerant populations.

Genome-Wide Analysis of Alternative Splicing Provides Insights Into Stress Response of the Pacific White Shrimp Litopenaeus vanname

Alternative splicing (AS) can enhance transcript diversity dramatically and play an important role in stress adaptation. Limited researches of AS have been reported in the Pacific white shrimp (Litopenaeus vannamei), which is an important aquaculture species in the world. Here, we performed a genome-wide identification of AS events in L. vannamei based on eight transcriptomes. We identified 38,781 AS events in the shrimp genome, and some of them were validated by polymerase chain reaction experiments. These AS events correspond to 9,209 genes, accounting for 36% of protein-coding genes in the shrimp genome. The number of AS events increased after virus or bacteria infection and low salinity stress. Type 1 AS genes (AS was initially activated) were mainly enriched in substance and energy metabolism, such as carbon metabolism and amino metabolism. However, type 2 AS genes (AS events changed) displayed specific enrichment under different stress challenges. Specifically, type 2 AS genes under biotic stresses were mainly enriched in the pathogenic pathway and immune network, and the AS genes under low salinity stress were significantly enriched for betalain biosynthesis. In summary, our study indicates that AS events are complex in shrimp and may be related to stress adaptation. These results will provide valuable resource for functional genomic studies on crustaceans.

Patterns of alternative splicing in response to cold acclimation in fish

Phenotypic plasticity is an important aspect of an organism's response to environmental change that often requires the modulation of gene expression. These changes in gene expression can be quantitative, as a result of increases or decreases in the amounts of specific transcripts, or qualitative, as a result of the expression of alternative transcripts from the same gene (e.g. via alternative splicing of pre-mRNAs). Although the role of quantitative changes in gene expression in phenotypic plasticity is well known, relatively few studies have examined the role of qualitative changes. Here, we use skeletal muscle RNA-seq data from Atlantic killifish (Fundulus heteroclitus), threespine stickleback (Gasterosteus aculeatus) and zebrafish (Danio rerio) to investigate the extent of qualitative changes in gene expression in response to cold acclimation. Fewer genes demonstrated alternative splicing than differential expression as a result of cold acclimation; however, differences in splicing were detected for 426 to 866 genes depending on species, indicating that large numbers of qualitative changes in gene expression are associated with cold acclimation. Many of these alternatively spliced genes were also differentially expressed, and there was functional enrichment for involvement in muscle contraction among the genes demonstrating qualitative changes in response to cold acclimation. Additionally, there was a common group of 29 genes with cold-acclimation-mediated changes in splicing in all three species, suggesting that there may be a set of genes with expression patterns that respond qualitatively to prolonged exposure to cold temperatures across fishes.

Effects of cold acclimation and dsRNA injections on Gs1l gene splicing in Drosophila montana

Alternative splicing, in which one gene produce multiple transcripts, may influence how adaptive genes respond to specific environments. A newly produced transcriptome of Drosophila montana shows the Gs1-like (Gs1l) gene to express multiple splice variants and to be down regulated in cold acclimated flies with increased cold tolerance. Gs1l’s effect on cold tolerance was further tested by injecting cold acclimated and non-acclimated flies from two distantly located northern and southern fly populations with double stranded RNA (dsRNA) targeting Gs1l. While both populations had similar cold acclimation responses, dsRNA injections only effected the northern population. The nature of splicing expression was then investigated in the northern population by confirming which Gs1l variants are present, by comparing the expression of different gene regions and by predicting the protein structures of splices using homology modelling. We find different splices of Gs1l not only appear to have independent impacts on cold acclimation but also elicit different effects in populations originating from two very different environments. Also, at the protein level, Gs1l appears homologous to the human HDHD1A protein and some splices might produce functionally different proteins though this needs to be verified in future studies by measuring the particular protein levels. Taken together, Gs1l appears to be an interesting new candidate to test how splicing influences adaptations.

Comparative Transcriptome Analysis Reveals Growth-Related Genes in Juvenile Chinese Sea Cucumber, Russian Sea Cucumber, and Their Hybrids

Heterosis is important for sea cucumber breeding, but its molecular mechanism remains largely unexplored. In this study, parental lines of Apostichopus japonicus from Russia (R) and China (C) were used to construct hybrids (CR and RC) by reciprocal crossing. We examined the transcriptional profiles of the hybrids (CR and RC) and the purebreds (CC and RR) at different developmental times. A total of 60.27 Gb of clean data was obtained, and 176,649 unigenes were identified, of which 50,312 unigenes were annotated. A total of 414,536 SNPs were identified. A total of 7011 differentially expressed genes (DEGs) were obtained between the purebreds and hybrids at 45 days after fertilization (DAF), and a total of 8218 DEGs were obtained between the purebreds and hybrids at 75 DAF. In addition, a total of 7652 DEGs were obtained between 45 DAF and 75 DAF. The significant DEGs were mainly involved in the MAPK and FOXO signaling pathways, especially in the Ras-Raf-MEK1/2-ERK module, which may be a key regulator of development and growth in juvenile A. japonicus. In addition, we also identified key growth-related genes, such as fgfs, igfs, megfs and hgfs, which were upregulated in the hybrids (RC and CR); these genes may play important roles in heterosis in A. japonicus. Our study provides fundamental information on the molecular mechanisms underlying heterosis in sea cucumber and might suggest strategies for the selection of rapidly growing strains of sea cucumber in aquaculture.

A novel and stress adaptive alternative oxidase derived from alternative splicing of duplicated exon in oyster Crassostrea virginica

Alternative oxidase (AOX) is a mitochondrial inner-membrane oxidase that accepts electrons directly from ubiquinol and reduces oxygen to water without involving cytochrome-linked electron transport chain. It is highly conserved in many non-vertebrate taxa and may protect cells against hypoxia and oxidative stress. We identified two AOX mRNAs in eastern oyster Crassostrea virginica, CvAOXA and CvAOXB, which differ by 170 bp but encode AOXs of the same size. Sequence analyses indicate that CvAOX has 10 exons with a tandem duplication of exon 10, and 3' alternative splicing using either the first or second exon 10 produces the two variants CvAOXB or CvAOXA, respectively. The second exon 10 in CvAOXA is more conserved across taxa, while the first exon 10 in CvAOXB contains novel mutations surrounding key functional sites. Both variants are expressed in all organs with the expression of CvAOXA higher than that of CvAOXB under normal condition. Under stress by air exposure, CvAOXB showed significantly higher expression than CvAOXA and became the dominant variant. This is the first case of alternative splicing of duplicated exon in a mollusc that produces a novel variant adaptive to stress, highlighting genome's versatility in generating diversity and phenotypic plasticity.