Comparative transcriptome analysis of papilla and skin in the sea cucumber, Apostichopus japonicus

Fine-mapping and association analysis of candidate genes for papilla number in sea cucumber, Apostichopus japonicus

The papilla number is one of the most economically important traits of sea cucumber in the China marketing trade. However, the genetic basis for papilla number diversity in holothurians is still scarce. In the present study, we conducted genome-wide association studies (GWAS) for the trait papilla number of sea cucumbers utilizing a set of 400,186 high-quality SNPs derived from 200 sea cucumbers. Two significant trait-associated SNPs that passed Bonferroni correction (P < 1.25E-7) were located in the intergenic region near PATS1 and the genic region of EIF4G, which were reported to play a pivotal role in cell growth and proliferation. The fine-mapping regions around the top two lead SNPs provided precise causative loci/genes related to papilla formation and cellular activity, including PPP2R3C, GBP1, and BCAS3. Potential SNPs with P < 1E-4 were acquired for the following GO and KEGG enrichment analysis. Moreover, the two lead SNPs were verified in another population of sea cucumber, and the expressive detection of three potential candidate genes PATS1, PPP2R3C, and EIF4G that near or cover the two lead SNPs was conducted in papilla tissue of TG (Top papilla number group) and BG (Bottom papilla number group) by qRT-PCR. We found the significantly higher expression profile of PATS1 (3.34-fold), PPP2R3C (4.90-fold), and EIF4G (4.23-fold) in TG, implying their potential function in papilla polymorphism. The present results provide valuable information to decipher the phenotype differences of the papilla trait and will provide a scientific basis for selective breeding in sea cucumbers.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-022-00139-w.

Analysis of the gene transcription patterns and DNA methylation characteristics of triploid sea cucumbers (Apostichopus japonicus)

Breeding of polyploid aquatic animals is still an important approach and research hotspot for realizing the economic benefits afforded by the improvement of aquatic animal germplasm. To better understand the molecular mechanisms of the growth of triploid sea cucumbers, we performed gene expression and genome-wide comparisons of DNA methylation using the body wall tissue of triploid sea cucumbers using RNA-seq and MethylRAD-seq technologies. We clarified the expression pattern of triploid sea cucumbers and found no dosage effect. DEGs were significantly enriched in the pathways of nucleic acid and protein synthesis, cell growth, cell division, and other pathways. Moreover, we characterized the methylation pattern changes and found 615 differentially methylated genes at CCGG sites and 447 differentially methylated genes at CCWGG sites. Integrative analysis identified 23 genes (such as Guf1, SGT, Col5a1, HAL, HPS1, etc.) that exhibited correlations between promoter methylation and expression. Altered DNA methylation and expression of various genes suggested their roles and potential functional interactions in the growth of triploid sea cucumbers. Our data provide new insights into the epigenetic and transcriptomic alterations of the body wall tissue of triploid sea cucumbers and preliminarily elucidate the molecular mechanism of their growth, which is of great significance for the breeding of fine varieties of sea cucumbers.

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.

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.

Reference assembly and gene expression analysis of Apostichopus japonicus larval development

The transcriptome of the holothurian Apostichopus japonicus was sequenced at four developmental stages—blastula, gastrula, auricularia, pentactula—on an Illumina sequencer. Based on our RNA-seq data and the paired-end reads from 16 libraries obtained by other researchers earlier, we have achieved the currently most complete transcriptome assembly for A. japonicus with the best basic statistical parameters. An analysis of the obtained transcriptome has revealed 174 differentially expressed transcription factors, as well as stage-specific transcription factors that are most promising for further study. In addition, a total of 1,174,999 high-quality single nucleotide polymorphisms have been identified, including 58,932 indels. A GO enrichment analysis of contigs containing polymorphic loci shows the predominance of GO terms associated with immune response. The data obtained by us provide an additional basis for a deeper study of the mechanisms of the planktotrophic-type development in holothurians and can be used in commercial sea cucumber breeding programs.

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.