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

Sea cucumbers: an emerging system in evo-devo

A challenge for evolutionary developmental (evo-devo) biology is to expand the breadth of research organisms used to investigate how animal diversity has evolved through changes in embryonic development. New experimental systems should couple a relevant phylogenetic position with available molecular tools and genomic resources. As a phylum of the sister group to chordates, echinoderms extensively contributed to our knowledge of embryonic patterning, organ development and cell-type evolution. Echinoderms display a variety of larval forms with diverse shapes, making them a suitable group to compare the evolution of embryonic developmental strategies. However, because of the laboratory accessibility and the already available techniques, most studies focus on sea urchins and sea stars mainly. As a comparative approach, the field would benefit from including information on other members of this group, like the sea cucumbers (holothuroids), for which little is known on the molecular basis of their development. Here, we review the spawning and culture methods, the available morphological and molecular information, and the current state of genomic and transcriptomic resources on sea cucumbers. With the goal of making this system accessible to the broader community, we discuss how sea cucumber embryos and larvae can be a powerful system to address the open questions in evo-devo, including understanding the origins of bilaterian structures.

Integration of mRNA and miRNA Profiling Reveals Heterosis in Oreochromis niloticus × O. aureus Hybrid Tilapia

Heterosis is a widespread biological phenomenon in fishes, in which hybrids have superior traits to parents. However, the underlying molecular basis for heterosis remains uncertain. Heterosis in growth and survival rates is apparent in hybrid tilapia (Oreochromis niloticus ♀ × O. aureus ♂). Comparisons of growth and hematological biochemical characteristics and mRNA and miRNA transcriptional analyses were performed in hybrid and parents tilapia stocks to investigate the underlying molecular basis for heterosis. Growth characteristics and hematological glucose and cholesterol parameters were significantly improved in hybrids. Of 3097 differentially expressed genes (DEGs) and 120 differentially expressed miRNAs (DEMs) identified among three stocks (O. niloticus, O. aureus, and hybrids), 1598 DEGs and 62 DEMs were non-additively expressed in hybrids. Both expression level dominance and overdominance patterns occurred for DEGs and DEMs, indicating that dominance and overdominance models are widespread in the transcriptional and post-transcriptional regulation of genes involved in growth, metabolism, immunity, and antioxidant capacity in hybrid tilapia. Moreover, potential negative regulation networks between DEMs and predicted target DEGs revealed that most DEGs from miRNA-mRNA pairs are up-regulated. Dominance and overdominance models in levels of transcriptome and miRNAome facilitate the integration of advantageous parental alleles into hybrids, contributing to heterosis of growth and improved survival. The present study provides new insights into molecular heterosis in hybrid tilapia, advancing our understanding of the complex mechanisms involved in this phenomenon in aquatic animals.

Transcriptome Profiling Revealed Basis for Growth Heterosis in Hybrid Tilapia (Oreochromis niloticus ♀ × O. aureus ♂)

Hybrid tilapia were produced from hybridization of Nile tilapia (Oreochromis niloticus) and blue tilapia (O. aureus). Comparative transcriptome analysis was carried out on the liver of hybrid tilapia and their parents by RNA sequencing. A total of 2319 differentially expressed genes (DEGs) were identified. Trend co-expression analysis showed that non-additive gene expression accounted for 67.1% of all DEGs. Gene Ontology and KEGG enrichment analyses classified the respective DEGs. Gene functional enrichment analysis indicated that most up-regulated genes, such as FASN, ACSL1, ACSL3, ACSL6, ACACA, ELOVL6, G6PD, ENO1, GATM, and ME3, were involved in metabolism, including fatty acid biosynthesis, unsaturated fatty acid biosynthesis, glycolysis, pentose phosphate pathway, amino acid metabolism, pyruvate metabolism, and the tricarboxylic acid cycle. The expression levels of a gene related to ribosomal biosynthesis in eukaryotes, GSH-Px, and those associated with heat shock proteins (HSPs), such as HSPA5 and HSP70, were significantly down-regulated compared with the parent tilapia lineages. The results revealed that the metabolic pathway in hybrid tilapia was up-regulated, with significantly improved fatty acid metabolism and carbon metabolism, whereas ribosome biosynthesis in eukaryotes and basal defense response were significantly down-regulated. These findings provide new insights into our understanding of growth heterosis in hybrid tilapia.

Transcriptome analysis reveals the molecular mechanisms of heterosis on thermal resistance in hybrid abalone

BACKGROUND

Heterosis has been exploited for decades in different animals and crops due to it resulting in dramatic increases in yield and adaptability. Hybridization is a classical breeding method that can effectively improve the genetic characteristics of organisms through heterosis. Abalone has become an increasingly economically important aquaculture resource with high commercial value. However, due to changing climate, abalone is now facing serious threats of high temperature in summer. Interspecific hybrid abalone (Haliotis gigantea ♀ × H. discus hannai ♂, SD) has been cultured at large scale in southern China and has been shown high survival rates under heat stress in summer. Therefore, SD has become a good model material for heterosis research, but the molecular basis of heterosis remains elusive.

RESULTS

Heterosis in thermal tolerance of SD was verified through Arrhenius break temperatures (ABT) of cardiac performance in this study. Then RNA-Sequencing was conducted to obtain gene expression patterns and alternative splicing events at control temperature (20 °C) and heat stress temperature (30 °C). A total of 356 (317 genes), 476 (435genes), and 876 (726 genes) significantly diverged alternative splicing events were identified in H. discus hannai (DD), H. gigantea (SS), and SD in response to heat stress, respectively. In the heat stress groups, 93.37% (20,512 of 21,969) of the expressed genes showed non-additive expression patterns, and over-dominance expression patterns of genes account for the highest proportion (40.15%). KEGG pathway enrichment analysis showed that the overlapping genes among common DEGs and NAGs were significantly enriched in protein processing in the endoplasmic reticulum, mitophagy, and NF-κB signaling pathway. In addition, we found that among these overlap genes, 39 genes had undergone alternative splicing events in SD. These pathways and genes may play an important role in the thermal resistance of hybrid abalone.

CONCLUSION

More alternative splicing events and non-additive expressed genes were detected in hybrid under heat stress and this may contribute to its thermal heterosis. These results might provide clues as to how hybrid abalone has a better physiological regulation ability than its parents under heat stress, to increase our understanding of heterosis in abalone.

Distant hybrids of Heliocidaris crassispina (♀) and Strongylocentrotus intermedius (♂): identification and mtDNA heteroplasmy analysis

Background

Distant hybridization between the sea urchin Heliocidaris crassispina (♀) and the sea urchin Strongylocentrotus intermedius (♂) was successfully performed under laboratory conditions. A new variety of hybrid sea urchin (HS hybrid) was obtained. However, the early-development success rates for the HS hybrids were significantly lower than those of purebred H. crassispina or S. intermedius offspring. In addition, it was difficult to distinguish the HS-hybrid adults from the pure H. crassispina adults, which might lead to confusion in subsequent breeding attempts. In this study, we attempted to develop a method to quickly and effectively identify HS hybrids, and to preliminarily investigate the molecular mechanisms underlying the poor early-development success rates in the HS hybrids.

Results

The hybrid sea urchins (HS hybrids) were identified both morphologically and molecularly. There were no significant differences in the test height to test diameter ratios between the HS hybrids and the parents. The number and arrangement of ambulacral pore pairs in the HS hybrids differed from those of the parental lines, which might serve as a useful morphological character for the identification of the HS hybrids. A primer pair that identified the HS hybrids was screened by comparing the mitochondrial genomes of the parental lines. Moreover, paternal leakage induced mitochondrial DNA heteroplasmy in the HS hybrids, which might explain the low rates of early development success in these hybrids.

Conclusions

The distant-hybrid sea urchins were accurately identified using comparative morphological and molecular genetic methods. The first evidence of mtDNA heteroplasmy after the distant hybridization of an echinoderm was also provided.

Identification of molecular markers for superior quantitative traits in a novel sea cucumber strain by comparative microRNA-mRNA expression profiling

To investigate the adaptability of Apostichopus japonicus (A. japonicus) strain "Anyuan No. 1" in the South China Sea, field monitoring and microRNA-mRNA integrated analyses were conducted between "Anyuan No. 1" and a regular A. japonicus population from Wendeng (Shandong Province, as a control) in the Xiapu farming area in Fujian Province, China. The results showed that "Anyuan No. 1" exhibited greater body weight increase and a higher number of papillae compared to the control during two and a half months of field monitoring. Comparative microRNA (miRNA) and mRNA transcriptome analyses identified 12 differentially expressed miRNAs (DEMs) and 165 differentially expressed genes (DEGs) in "Anyuan No. 1" compared to the control. Long-chain specific acyl-CoA dehydrogenase (ACADL), transmembrane protein 251 (TMEM251), dehydrogenase/reductase SDR family protein 7-like (Dhrs7), insulin-like growth factor-binding protein 7 (IGFBP-7), CDK5 regulatory subunit-associated protein 1 (CDK5RAP1), visual pigment-like receptor peropsin, 39S ribosomal protein, miR-10, miR-153, miR-7, and miR-3529 were identified as gene and miRNA candidates correlated with superior economic traits in "Anyuan No. 1". Collectively, "Anyuan No. 1" is suitable for large-scale cultivation extension due to its better adaptability to the South China Sea area. Furthermore, we identified "miR10-ACADL" as a potential module for further molecular marker-assisted selective breeding of A. japonicus.

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.

Transcriptome analysis of body wall reveals growth difference between the largest and smallest individuals in the pure and hybrid populations of Apostichopus japonicus

Long-term inbreeding of sea cucumber has resulted in a decrease in its growth rate, which has severely affected yield and economic efficiency. In this study, three Apostichopus japonicus families were constructed and screened into the weight of smallest and largest, which included Russian, Chinese, and their hybrids (RC). We examined the transcriptional profiles of hybrid (RC) and purebred (CC and RR). A total of 49.69 Gb clean reads were obtained, and the Q30 base percentage was above 90.47%. A total of 5191 novel genes were discovered, of which 2592 genes were annotated. Differentially expressed genes (DEGs) were identified, and functional annotation and enrichment analysis were performed. Approximately 1874 DEGs were screened in the Chinese sea cucumber (CC) difference group; 2591 DEGs were obtained in the hybrid sea cucumber difference group (RC), and 3006 DEGs were obtained in the Russian sea cucumber difference group (RR). In Gene Ontology (GO) analysis, highest DEG enrichment was observed for the functional categories of cellular process and metabolic process. In terms of cellular components, DEG enrichment was observed in cell part, cell; for molecular function, DEG enrichment was detected in catalytic activity, binding, hydrolase activity, transferase activity. According to the differential expression analysis, we found that 15 heat shock protein (HSP) genes that have the same expression trends, which were upregulated in the smallest weight of three sea cucumber lines. In addition, COG analysis of defense genes was conducted. All defense genes (ATP-binding cassette transporters (ABCs), multidrug resistance protein (MRPs), and beta-lactamase) showed the same expression trend, which was significantly upregulated in smallest individuals compared to that of largest individuals in RC lines, which implied the smallest individuals are exposed to more pressure during growth. These results may lead to the smallest individuals showing slow growth. Additionally, we selected 12 DEGs to validate the result by qPCR. Those DEGs were included in growth-related and resistance genes. Sequencing of the A. japonicus transcriptome improves our understanding of the transcriptional regulatory apparatus that controls individual development and growth.