Transcriptome profiling analysis on whole bodies of microbial challenged Eriocheir sinensis larvae for immune gene identification and SNP development

Comparative transcriptome, digital gene expression and proteome profiling analyses provide insights into the brachyurization from the megalopa to the first juvenile in Eriocheir sinensis

Eriocheir sinensis larva normally experiences 11 stages. The reduced abdomen folded beneath the thorax is the most prominent characteristic of morphological alteration from megalopa to juvenile crab. Up to date, the molecular mechanisms of brachyurization remain a mystery. Here, transcriptome library, digital gene expression (DGE) libraries and proteome libraries at two developmental stages [the megalopa stage of E. sinensis (stage M) and the first stage of juvenile crab (stage J1)] of the Chinese mitten crab larva were constructed for RNA sequencing and iTRAQ approaches followed by bioinformatics analysis, respectively. In total, 1106 genes and 871 proteins were differentially expressed between the stage M and stage J1. Moreover, several important pathways were identified, including biosynthesis of secondary metabolites, metabolic pathways, focal adhesion, and some disease pathways. Besides, muscle contraction, oxidative phosphorylation, calcium signaling, PI3K-Akt, DNA replication pathway, and integrin signaling pathway also had important functions in brachyurization process. Furthermore, the components, actin, actin-related protein, collagens, filamin-A/B, laminin, integrins, paxillin, and fibronectin had up-regulated expression levels in M stage compared to J1 stage.

Toxic effects of metal copper stress on immunity, metabolism and pathologic changes in Chinese mitten crab (Eriocheir japonica sinensis)

Copper (Cu²⁺), which represents a major physiological challenge for crab culture, is ubiquitous in the aquatic culture environment, and gills are the first organs that come into direct contact with the environment. However, the molecular basis of the response of crabs to Cu²⁺ stress remains unclear. Here, we conducted a transcriptome and differential expression analysis on the gills from Chinese mitten crab unexposed and exposed to Cu²⁺ for 24 h. The comparative transcriptome analysis identified 2486 differentially expressed genes (DEGs). GO functional analysis and KEGG pathway analysis revealed some DEGs, which were mostly related to immunity, metabolism, osmotic regulation, Cu²⁺ homeostasis regulation, antioxidant activity, and detoxification process. Some pathways related to humoral and cellular immunity, such as phagosome, peroxisome, lysosome, mTOR signaling pathway, PI3K-Akt signaling pathway, Toll-like receptor signaling pathway, and T cell receptor signaling pathway were enhanced under Cu²⁺ stress. In addition, Cu²⁺ stress altered the expression patterns of key phagocytosis and apoptosis genes (lectin, cathepsin L, Rab7, and HSP70), confirming that Cu²⁺ can induce oxidative stress and eventually even apoptosis. Histological analysis revealed that the copper can induce damage at the cellular level. This comparative transcriptome analysis provides valuable molecular information to aid future study of the immune mechanism of Chinese mitten crab in response to Cu²⁺ stress and provides a foundation for further understanding of the effects of metal toxicity.

High efficiency CRISPR/Cas9 genome editing system with an eliminable episomal sgRNA plasmid in Pichia pastoris

Pichia pastoris is a methylotrophic yeast in which host heterologous expression of proteins has been developed owing to the strong inducible alcohol oxidase promoter (PAOX1). However, it is difficult to manipulate the genome in P. pastoris. Based on previous attempts to apply the CRISPR/Cas9 system in P. pastoris, a CRISPR/Cas9 system with episomal sgRNA plasmid was developed and 100 % genome editing efficiency, high multicopy gene editing and stable multigene editing were obtained without a sharp decline caused by multi-sgRNA. And 28/34 (∼82 %) sgRNAs tested were effective. The CGG may have a slightly higher and more stable cleavage efficiency than the other three NGG motifs, and a low GC content may be preferable for higher cleavage efficiency. This provides researchers with a stable genome editing tool that shows a high editing efficiency, shortening the experimentation period. Furthermore, we introduced dCas9 into P. pastoris and achieved target gene interference, expanding the CRISPR/Cas9 toolbox in P. pastoris.

The Single-molecule long-read sequencing of Scylla paramamosain

Scylla paramamosain is an important aquaculture crab, which has great economical and nutritional value. To the best of our knowledge, few full-length crab transcriptomes are available. In this study, a library composed of 12 different tissues including gill, hepatopancreas, muscle, cerebral ganglion, eyestalk, thoracic ganglia, intestine, heart, testis, ovary, sperm reservoir, and hemocyte was constructed and sequenced using Pacific Biosciences single-molecule real-time (SMRT) long-read sequencing technology. A total of 284803 full-length non-chimeric reads were obtained, from which 79005 high-quality unique transcripts were obtained after error correction and sequence clustering and redundant. Additionally, a total of 52544 transcripts were annotated against protein database (NCBI nonredundant, Swiss-Prot, KOG, and KEGG database). A total of 23644 long non-coding RNAs (lncRNAs) and 131561 simple sequence repeats (SSRs) were identified. Meanwhile, the isoforms of many genes were also identified in this study. Our study provides a rich set of full-length cDNA sequences for S. paramamosain, which will greatly facilitate S. paramamosain research.

Transcriptome Profile Analysis on Ovarian Tissues of Autotetraploid Fish and Diploid Red Crucian Carp

Polyploidization can significantly alter the size of animal gametes. Autotetraploid fish (RRRR, 4nRR = 200) (4nRR) possessing four sets of chromosomes were derived from whole-genome duplication in red crucian carp (RR, 2n = 100) (RCC). The diploid eggs of the 4nRR fish were significantly larger than the eggs of RCC. To explore the differences between the ovaries of these two ploidies of fishes at the molecular level, we compared the ovary transcriptome profiles of 4nRR fish and RCC and identified differentially expressed genes (DEGs). A total of 19,015 unigenes were differentially expressed between 4nRR fish and RCC, including 12,591 upregulated and 6,424 downregulated unigenes in 4nRR fish. Functional analyses revealed that eight genes (CDKL1, AHCY, ARHGEF3, TGFβ, WNT11, CYP27A, GDF7, and CKB) were involved in the regulation of cell proliferation, cell division, gene transcription, ovary development and energy metabolism, suggesting that these eight genes were related to egg size in 4nRR fish and RCC. We validated the expression levels of these eight DEGs in 4nRR fish and RCC using quantitative PCR. The study results provided insights into the regulatory mechanisms underlying the differences in crucian carp egg sizes.

Comparative transcriptome analysis explores maternal to zygotic transition during Eriocheir sinensis early embryogenesis

The maternal genome directs almost all aspects of early animal development. As development proceeds, the elimination of maternal gene products and zygotic genome activation (ZGA) occur during the maternal to zygotic transition (MZT). To study the molecular mechanisms regulating this developmental event in Eriocheir sinensis, RNA-Seq technology was applied to generate comprehensive information on transcriptome dynamics during early embryonic stages. In total, 32,088 annotated unigenes were obtained from the transcriptomes of fertilized eggs and embryos at the cleavage (2-4 cell) and blastula stage. A total of 566 maternal genes and 1165 zygotic genes were isolated, among which 103 and 266 genes were predicted conserved maternal transcripts (COMATs) and conserved zygotic transcripts (COZYTs), respectively. The COMATs performed housekeeping gene functions and may be essential for initiating early embryogenesis of the Bilateria. Furthermore, 87, 76 and 117 differentially expressed genes associated with the MZT, morphogenesis and immunity were identified when compared the three transcriptomic datasets. We also unmask that the MZT takes place around the cleavage stage, when the genes involved in the clearance of maternal gene products and the ZGA were significantly up-regulated. Taken together, these datasets provide a valuable resource for understanding the mechanisms of early developmental events in E. sinensis, and facilitate further studies on molecular mechanisms of asynchronous development in crabs.

Transcriptome analysis of mantle tissues reveals potential biomineralization-related genes in Tectus pyramis Born

The marine mollusk Tectus pyramis is a valuable shellfish primarily distributed in the tropical waters of the South China Sea, as well as in the Indo-Pacific Ocean and areas near the southern portion of the Japanese Peninsula. Despite major economic interest in this mollusk, limited genomic resources are available for this species, which has prevented studies of the molecular mechanism, such as biomineralization. Here, we report the first comprehensive transcript dataset of T. pyramis mantle tissue. From a total of 16,801,141 reads, 173,671 unique transcripts were assembled, which provides new genomic resources for the understanding of biomineralization in T. pyramis. The most abundant unique sequences of the top 30 most highly expressed genes were annotated as shematrin, while other highly expressed genes included glycine-rich protein and shematrin-1. Based on transcriptome annotation and Gene Ontology classification, 130 biomineralization-related genes were found including members of the BMP (bone morphogenetic proteins), calmodulin, perlucin, and shematrin families, as well as mantle genes, nacrein, and MSI60. The results of qPCR showed that 14 of 24 examined genes were highly expressed in the mantle. A phylogenetic tree of BMP, perlucin, shematrin proteins revealed conservation of their structure and functions and indicated that some members participated in biomineralization in T. pyramis. Taken together, the results presented herein will be useful in studies of molecular mechanisms and pathways of biomineralization in T. pyramis, as well as provide new insight into the mechanisms of biomineralization in gastropods.

Transcriptome analysis of the Yesso scallop, Patinopecten yessoensis gills in response to water temperature fluctuations

Water temperature fluctuations are considered to be a major factor affecting the immune functions and metabolic processes of scallops. To better understand the immune defense mechanisms of Yesso scallop, Patinopecten yessoensis following exposure to water temperature fluctuations, transcriptomic profiles in the gills from high-frequency fluctuations (HF_G), low-frequency fluctuations (LF_G), and no fluctuations (NF_G) groups were obtained using HiSeq™ 2500 (Illumina). For HF_G, scallops were transferred directly between 18 and 8 °C every 4 h and for 10 fluctuations, while scallops in LF_G were transferred between 18 and 13 °C every 12 h, for a total of 4 fluctuations. A total of 442,922,590 clean reads were generated in 9 libraries and then assembled into 210,780 unigenes with an average length of 705 bp and an N50 of 1253 bp. Based on sequence similarity, 54,529 unigenes (25.87%) were annotated in at least one database. Comparative analysis revealed that 696 unigenes differentially expressed in temperature stressed groups compared with the control, including 229 unigenes between HF_G and NF_G, and 548 unigenes between LF_G and NF_G, respectively. Additionally, among these differentially expressed genes (DEGs), there were 41 immune-related unigenes and 16 protein metabolism-related unigenes. These results provide fundamental information on the molecular defense mechanisms in the Yesso scallop gills after exposure to water temperature fluctuations.

Transcriptome-seq provides insights into sex-preference pattern of gene expression between testis and ovary of the crucifix crab (Charybdis feriatus)

The crucifix crab, Charybdis feriatus, which mainly inhabits Indo-Pacific region, is regarded as one of the most high-potential species for domestication and incorporation into the aquaculture sector. However, the regulatory mechanisms of sex determination and differentiation of this species remain unclear. To identify candidate genes involved in sex determination and differentiation, high throughput sequencing of transcriptome from the testis and ovary of C. feriatus was performed by the Illumina platform. After removing adaptor primers, low-quality sequences and very short (<50 nt) reads, we obtained 80.9 million and 66.2 million clean reads from testis and ovary, respectively. A total of 86,433 unigenes were assembled, and ~43% (37,500 unigenes) were successfully annotated to the NR, NT, Swiss-Prot, KEGG, COG, GO databases. By comparing the testis and ovary libraries, we obtained 27,636 differentially expressed genes. Some candidate genes involved in the sex determination and differentiation of C. feriatus were identified, such as vasa, pgds, vgr, hsp90, dsx-f, fem-1, and gpr. In addition, 88,608 simple sequence repeats were obtained, and 61,929 and 77,473 single nucleotide polymorphisms from testis and ovary were detected, respectively. The transcriptome profiling was validated by quantitative real-time PCR in 30 selected genes, which showed a good consistency. The present study is the first high-throughput transcriptome sequencing of C. feriatus. These findings will be useful for future functional analysis of sex-associated genes and molecular marker-assisted selections in C. feriatus.

Transcriptome Analysis and Differential Gene Expression on the Testis of Orange Mud Crab, Scylla olivacea, during Sexual Maturation

Adequate genetic information is essential for sustainable crustacean fisheries and aquaculture management. The commercially important orange mud crab, Scylla olivacea, is prevalent in Southeast Asia region and is highly sought after. Although it is a suitable aquaculture candidate, full domestication of this species is hampered by the lack of knowledge about the sexual maturation process and the molecular mechanisms behind it, especially in males. To date, data on its whole genome is yet to be reported for S. olivacea. The available transcriptome data published previously on this species focus primarily on females and the role of central nervous system in reproductive development. De novo transcriptome sequencing for the testes of S. olivacea from immature, maturing and mature stages were performed. A total of approximately 144 million high-quality reads were generated and de novo assembled into 160,569 transcripts with a total length of 142.2 Mb. Approximately 15–23% of the total assembled transcripts were annotated when compared to public protein sequence databases (i.e. UniProt database, Interpro database, Pfam database and Drosophila melanogaster protein database), and GO-categorised with GO Ontology terms. A total of 156,181 high-quality Single-Nucleotide Polymorphisms (SNPs) were mined from the transcriptome data of present study. Transcriptome comparison among the testes of different maturation stages revealed one gene (beta crystallin like gene) with the most significant differential expression—up-regulated in immature stage and down-regulated in maturing and mature stages. This was further validated by qRT-PCR. In conclusion, a comprehensive transcriptome of the testis of orange mud crabs from different maturation stages were obtained. This report provides an invaluable resource for enhancing our understanding of this species’ genome structure and biology, as expressed and controlled by their gonads.

Transcriptome Analysis for Identification of Genes Related to Gonad Differentiation, Growth, Immune Response and Marker Discovery in The Turbot (Scophthalmus maximus)

Background

Turbot Scophthalmus maximus is an economically important species extensively aquacultured in China. The genetic selection program is necessary and urgent for the sustainable development of this industry, requiring more and more genome background knowledge. Transcriptome sequencing is an excellent alternative way to identify transcripts involved in specific biological processes and exploit a considerable quantity of molecular makers when no genome sequences are available. In this study, a comprehensive transcript dataset for major tissues of S. maximus was produced on basis of an Illumina platform.

Results

Total RNA was isolated from liver, spleen, kidney, cerebrum, gonad (testis and ovary) and muscle. Equal quantities of RNA from each type of tissues were pooled to construct two cDNA libraries (male and female). Using the Illumina paired-end sequencing technology, nearly 44.22 million clean reads in length of 100 bp were generated and then assembled into 106,643 contigs, of which 71,107 were named unigenes with an average length of 892 bp after the elimination of redundancies. Of these, 24,052 unigenes (33.83% of the total) were successfully annotated. GO, KEGG pathway mapping and COG analysis were performed to predict potential genes and their functions. Based on our sequence analysis and published documents, many candidate genes with fundamental roles in sex determination and gonad differentiation (dmrt1), growth (ghrh, myf5, prl/prlr) and immune response (TLR1/TLR21/TLR22, IL-15/IL-34), were identified for the first time in this species. In addition, a large number of credible genetic markers, including 21,192 SSRs and 8,642 SNPs, were identified in the present dataset.

Conclusion

This informative transcriptome provides valuable new data to increase genomic resources of Scophthalmus maximus. The future studies of corresponding gene functions will be very useful for the management of reproduction, growth and disease control in turbot aquaculture breeding programs. The molecular markers identified in this database will aid in genetic linkage analyses, mapping of quantitative trait loci, and acceleration of marker assisted selection programs.

High-density linkage mapping aided by transcriptomics documents ZW sex determination system in the Chinese mitten crab Eriocheir sinensis

The sex determination system in crabs is believed to be XY-XX from karyotypy, but centromeres could not be identified in some chromosomes and their morphology is not completely clear. Using quantitative trait locus mapping of the gender phenotype, we revealed a ZW-ZZ sex determination system in Eriocheir sinensis and presented a high-density linkage map covering ~98.5% of the genome, with 73 linkage groups corresponding to the haploid chromosome number. All sex-linked markers in the family we used were located on a single linkage group, LG60, and sex linkage was confirmed by genome-wide association studies (GWAS). Forty-six markers detected by GWAS were heterozygous and segregated only in the female parent. The female LG60 was thus the putative W chromosome, with the homologous male LG60 as the Z chromosome. The putative Z and W sex chromosomes were identical in size and carried many homologous loci. Sex ratio (5:1) skewing towards females in induced triploids using unrelated animals also supported a ZW-ZZ system. Transcriptome data were used to search for candidate sex-determining loci, but only one LG60 gene was identified as an ankyrin-2 gene. Double sex- and mab3-related transcription factor 1 (Dmrt1), a Z-linked gene in birds, was located on a putative autosome. With complete genome sequencing and transcriptomic data, more genes on putative sex chromosomes will be characterised, thus leading towards a comprehensive understanding of the sex determination and differentiation mechanisms of E. sinensis, and decapod crustaceans in general.

Characterization of the mantle transcriptome of yesso scallop (Patinopecten yessoensis): identification of genes potentially involved in biomineralization and pigmentation

The Yesso scallop Patinopecten yessoensis is an economically important marine bivalve species in aquaculture and fishery in Asian countries. However, limited genomic resources are available for this scallop, which hampers investigations into molecular mechanisms underlying their unique biological characteristics, such as shell formation and pigmentation. Mantle is the special tissue of P. yessoensis that secretes biomineralization proteins inducing shell deposition as well as pigmentation on the shells. However, a current deficiency of transcriptome information limits insight into mechanisms of shell formation and pigmentation in this species. In this study, the transcriptome of the mantle of P. yessoensis was deeply sequenced and characterized using Illumina RNA-seq technology. A total of 86,521 unique transcripts are assembled from 55,884,122 reads that passed quality filters, and annotated, using Gene Ontology classification. A total of 259 pathways are identified in the mantle transcriptome, including the calcium signaling and melanogenesis pathways. A total of 237 unigenes that are homologous to 102 reported biomineralization genes are identified, and 121 unigenes that are homologous to 93 known proteins related to melanin biosynthesis are found. Twenty-three annotated unigenes, which are mainly homologous to calmodulin and related proteins, Ca2+/calmodulin-dependent protein kinase, adenylate/guanylate cyclase, and tyrosinase family are potentially involved in both biomineralization and melanin biosynthesis. It is suggested that these genes are probably not limited in function to induce shell deposition by calcium metabolism, but may also be involved in pigmentation of the shells of the scallop. This potentially supports the idea that there might be a link between calcium metabolism and melanin biosynthesis, which was previously found in vertebrates. The findings presented here will notably advance the understanding of the sophisticated processes of shell formation as well as shell pigmentation in P. yessoensis and other bivalve species, and also provide new evidence on gene expression for the understanding of pigmentation and biomineralization not only in invertebrates but also probably in vertebrates.

Calmodulin is a stress and immune response gene in Chinese mitten crab Eriocheir sinensis

Calmodulin (CaM) is a multifunctional calcium sensor protein that participates in various cellular processes under normal, stress and pathological conditions. In crabs, however, the involvement of CaM in response to environmental stress and immune challenges has not been revealed yet. In the present study, a CaM cDNA (EsCaM) was identified from Chinese mitten crab Eriocheir sinensis and its mRNA expression patterns in response to ambient (salinity and pH) stress and immune challenges was examined. EsCaM encodes a 149-amino-acid protein with a calculated molecular mass of 16.8 kDa and an isoelectric point of 4.09. In unstimulated healthy E. sinensis, EsCaM mRNA transcript was detected in all tested tissues with predominant expression in hepatopancreas and the lowest expression in haemocytes. Ambient salinity (15‰ and 30‰ salinities) and pH (pH 6 and 8.5) stress significantly altered EsCaM mRNA expression in gill, hepatopancreas, haemocytes, intestine and muscle in Chinese mitten crab. In addition, EsCaM gene expression was significantly and rapidly induced as early as 2 h after LPS and Poly(I:C) immune stimulations in haemocytes in vitro. Furthermore, EsCaM expression was significantly up-regulated in E. sinensis haemocytes, gill, hepatopancreas, intestine and muscle in response to Edwardsiella tarda and Vibrio anguillarum challenges. Collectively, our findings suggest that EsCaM is an important stress and immune response gene in Chinese mitten crab.