De novo transcriptome sequencing and analysis of male and female swimming crab (Portunus trituberculatus) reproductive systems during mating embrace (stage II)

Overcoming research challenges: In vitro cultivation of Ameson portunus (Phylum Microsporidia)

Ameson portunus is an intracellular pathogen that infects marine crabs Portunus trituberculatus and Scylla paramamosain, causing significant economic losses. However, research into this important parasite has been limited due to the absence of an in vitro culture system. To address this challenge, we developed an in vitro cultivation model of A. portunus using RK13 cell line in this study. The fluorescent labeling assay indicated a high infection rate (∼60 %) on the first day post-infection and quantitative PCR (qPCR) detection demonstrated successful infection as early as six hours post-inoculation. Fluorescence in situ hybridization (FISH) and qPCR were used for the detection of A. portunus infected cells. The FISH probe we designed allowed detection of A. portunus in infected cells and qPCR assay provided accurate quantification of A. portunus in the samples. Transmission electron microscopy (TEM) images revealed that A. portunus could complete its entire life cycle and produce mature spores in RK13 cells. Additionally, we have identified novel life cycle characteristics during the development of A. portunus in RK 13 cells using TEM. These findings contribute to our understanding of new life cycle pathways of A. portunus. The establishment of an in vitro culture model for A. portunus is critical as it provides a valuable tool for understanding the molecular and immunological events that occur during infection. Furthermore, it will facilitate the development of effective treatment strategies for this intracellular pathogen.

Phenotypic sorting of individual male and female intersex Cherax quadricarinatus and analysis of molecular differences in the gonadal transcriptome

Cherax quadricarinatus exhibit sexual dimorphism, with males outpacing females in size specification and growth rate. However, there is limited understanding of the molecular mechanisms underlying sex determination and sex differentiation in crustaceans. To study the differences between intersex individuals and normal individuals, this study counted the proportion of intersex individuals in the natural population, collected the proportion of 7 different phenotypes in 200 intersex individuals, and observed the differences in tissue sections. RNA-seq was used to study the different changes in the transcriptome of normal and intersex gonads. The results showed that: the percentage of intersex in the natural population was 1.5 %, and the percentage of different types of intersex ranged from 0.5 % to 22.5 %; the sections revealed that the development of normal ovaries was stagnant at the primary oocyte stage when intersex individuals with ovaries were present; We screened for pathways and genes that may be associated with gonadal development and sex, including ovarian steroid synthesis, estrogen signaling pathway, oocyte meiosis, progesterone-mediated oocyte maturation, etc. Relevant genes including tra2a, dmrta2, ccnb2, foxl2, and smad4. This study provides an important molecular basis for sex determination, sex-controlled breeding, and unisex breeding in red crayfish.

Identifying sex-differential gene expression in the antennal gland of the swimming crab by transcriptomic analysis

The antennal glands (AnGs) are recognized as an important organ that functions in ion regulation and excretion in decapods. Previously, many studies had explored this organ at the biochemical, physiological, and ultrastructural levels but had few molecular resources. In this study, the transcriptomes of the male and female AnGs of Portunus trituberculatus were sequenced using RNA sequencing (RNA-Seq) technology. Genes involved in osmoregulation and organic/inorganic solute transport were identified. This suggests that AnGs might be involved in these physiological functions as versatile organs. A total of 469 differentially expressed genes (DEGs) were further identified between male and female transcriptomes and found to be male-biased. Enrichment analysis showed that females were enriched in amino acid metabolism and males were enriched in nucleic acid metabolism. These results suggested differences in possible metabolic patterns between males and females. Furthermore, two transcription factors related to reproduction, namely AF4/FMR2 family members Lilli (Lilli) and Virilizer (Vir), were identified in DEGs. Lilli was found to be specifically expressed in the male AnGs, whereas Vir showed high expression levels in the female AnGs. The expression of up-regulated metabolism and sexual development-related genes in three males and six females was verified by qRT-PCR and the pattern was found to be consistent with the transcriptome expression pattern. Our results suggest that although the AnG is a unified somatic tissue composed of individual cells, it still demonstrates distinct sex-specific expression patterns. These results provide foundational knowledge of the function and differences between male and female AnGs in P. trituberculatus.

Development of Microsatellite Markers Based on Transcriptome Sequencing and Evaluation of Genetic Diversity in Swimming Crab (Portunus trituberculatus)

P. trituberculatus is an economically important mariculture species in China. Evaluating its genetic diversity and population structure can contribute to the exploration of germplasm resources and promote sustainable aquaculture production. In this study, a total of 246,243 SSRs were generated by transcriptome sequencing of P. trituberculatus. Among the examined 254,746 unigenes, 66,331 had more than one SSR. Among the different SSR motif types, dinucleotide repeats (110,758, 44.98%) were the most abundant. In 173 different base repeats, A/T (96.86%), AC/GT (51.46%), and ACC/GGT (26.20%) were dominant in mono-, di-, and trinucleotide, respectively. GO annotations showed 87,079 unigenes in 57 GO terms. Cellular process, cell, and binding were the most abundant terms in biological process, cellular component, and molecular function categories separately. A total of 34,406 annotated unigenes were classified into 26 functional categories according to the functional annotation analysis of KOG, of which “general function prediction only” was the biggest category (6,028 unigenes, 17.52%). KEGG pathway annotations revealed the clustering of 34,715 unigenes into 32 different pathways. Nineteen SSRs were identified as polymorphic and, thus, used to assess the genetic diversity and structure of 240 P. trituberculatus individuals from four populations in the Bohai Sea. Genetic parameter analysis showed a similar level of genetic diversity within wild populations, and the cultured population indicated a reduction in genetic diversity compared with wild populations. The pairwise F_ST values were between 0.001 and 0.04 with an average of 0.0205 (p < 0.05), suggesting a low but significant level of genetic differentiation among the four populations. Structure analysis demonstrated that the four populations were classified into two groups including the cultured group and other populations. The phylogenetic tree and PCA revealed that a vast number of samples were clustered together and that cultivated individuals were distributed more centrally than wild individuals. The findings contribute to the further assessment of germplasm resources and assist to provide valuable SSRs for marker-assisted breeding of P. trituberculatus in the future.

Comparative Transcriptome Analysis Reveals the Process of Ovarian Development and Nutrition Metabolism in Chinese Mitten Crab, Eriocheir Sinensis

Ovarian development is a key physiological process that holds great significance in the reproduction of the Chinese mitten crab (Eriocheir sinensis), which is an economically important crab species for aquaculture. However, there is limited knowledge for the regulatory mechanisms of ovarian development. To study the molecular mechanisms of its ovarian development, transcriptome analysis was performed in the ovary and hepatopancreas of E. sinensis during ovarian stages I (oogonium proliferation), II (endogenous vitellogenesis), and III (exogenous vitellogenesis). The results showed that 5,520 and 226 genes were differentially expressed in the ovary and hepatopancreas, respectively. For KEGG enrichment analysis, the differentially expressed genes in the ovary were significantly clustered in phototransduction-fly, phagosome, and ECM-receptor interaction. Significantly enriched pathways in the hepatopancreas included fatty acid biosynthesis, fatty acid metabolism, and riboflavin metabolism. Further analysis showed that 25 genes and several pathways were mainly involved in oogenesis, including the ubiquitin-proteasome pathway, cyclic AMP-protein kinase A signaling pathway, and mitogen-activated protein kinase signaling pathway. Twenty-five candidate genes involved in vitellogenesis and endocrine regulation were identified, such as vitellogenin, vitellogenin receptor, estrogen sulfotransferase, ecdysone receptor, prostaglandin reductase 1, hematopoietic prostaglandin D synthase and juvenile hormone acid O-methyltransferase. Fifty-six genes related to nutritional metabolism were identified, such as fatty acid synthase, long-chain-fatty-acid-CoA ligase 4, 1-acyl-sn-glycerol-3-phosphate acyltransferase 4, fatty acid-binding protein, and glycerol-3-phosphate acyltransferase 1. These results highlight the genes involved in ovarian development and nutrition deposition, which enhance our understanding of the regulatory pathways and physiological processes of crustacean ovarian development.

De novo gonad transcriptome analysis of elongate loach (Leptobotia elongata) provides novel insights into sex-related genes

Elongate loach (Leptobotia elongata) is endemic in middle and upper reaches of the Yangtze River in China. Because of many anthropogenic factors such as overfishing and dam construction, the loach has become a highly endangered species. So far, the genomic resources which benefit for species conservation and utilization are still lacking in elongate loach. Therefore, the first gonad transcriptome analysis of the loach was conducted in this study, providing novel insights into sex-related genes. A total of 286,800,660 clean reads with a total length of 42.02 Gb were obtained. 18,975 differentially expressed genes (DEGs) were identified, where 12,976 DEGs, especially Sox9a, Sox9b, Igf2 and Fgfr2, were upregulated in the testis, and 5999 DEGs, especially Zp3, Eg2, Plk1, Ccnb1, Cdc20 and Mos, were upregulated in the ovary. Meanwhile, some testis-specific genes (i.e. Cald1, Atp1a, Muc2 and Ca2) and ovary-specific genes (i.e. Ca4, Tuba, Acp5, Ccna, Larp6 and Nop4) were identified and verified. According to Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of DEGs, we found a series of enrichment pathways related to reproduction in elongate loach, such as the MAPK signaling pathway, oxytocin signaling pathway and oocyte meiosis pathway. Twelve DEGs were randomly selected to verify RNA-seq results by qPCR. In conclusion, this study provides a data source to study the molecular characteristics and regulatory mechanisms of sex-related genes in elongate loach, which has a potential to improve the resource protection and aquaculture production of the loach.

Comparative Transcriptomics of Gonads Reveals the Molecular Mechanisms Underlying Gonadal Development in Giant Freshwater Prawns (Macrobrachium rosenbergii)

The giant freshwater prawn, Macrobrachium rosenbergii, is a prawn that has economic significance throughout the world. It exhibits sex-related growth dimorphism, whereby the males grow significantly more rapidly than the females. Therefore, a study on the molecular regulatory mechanism, which underlies the sexual differentiation of M. rosenbergii, is of both scientific and commercial importance. However, a scarcity of genomic and transcriptomic resources severely limits our knowledge of the sexual differentiation mechanisms in M. rosenbergii. Here, transcriptome sequencing of several gonadic samples of males and females in M. rosenbergii was performed to investigate the molecular basis underlying gonadal development. Our results showed that 2149 unigenes presented as differentially expressed genes (DEGs) in the ovaries of females compared to the testes of males, which contained 484 down-regulated and 1665 up-regulated genes. Enrichment analysis of DEGs revealed many of these genes to be related to sexual differentiation and gonadal development. From our transcriptome analyses, and as confirmed by quantitative real-time PCR, male-related genes (Mrr, MRPINK, IR, IAGBP, TESK1, and dsx) in the testes were significantly up-regulated, and female-related genes (ERR, Sxl3, cyclinB, Dmrt99B, PPP2A, and ADCY9) in the ovaries were also significantly up-regulated. This indicates the potential role these genes play in the gonadal development of M. rosenbergii. Furthermore, multiple signal transduction pathways relating to gonadal maturation and spermatogenesis, including MAPK, were identified herein. Our data also supports previous ideas that IAG and IAGBP-IR signaling schemes could help in the regulation of testis’ development in M. rosenbergii and the ERR gene could regulate ovarian development by affecting the expression of cyclinB, PPP2A, and ADCY9. The data from this study provides incredibly usefully genomic resources for future research on the sexual differentiation and practical aquaculture of M. rosenbergii.

Full-Length Transcriptome Analysis Provides New Insights Into the Diversity of Immune-Related Genes in Portunus trituberculatus

Generally, invertebrates were thought to solely rely on their non-specific innate immune system to fight against invading microorganisms. However, increasing studies have implied that the innate immune response of invertebrates displayed diversity and specificity owing to the hyper-variable immune molecules in organisms. In order to get an insight into the diversity of immune-related genes in Portunus trituberculatus, a full-length transcriptome analysis of several immune-related tissues (hemocytes, hepatopancreas and gills) in P. trituberculatus was performed and the diversity of several immune-related genes was analyzed. The full-length transcriptome analysis of P. trituberculatus was conducted using a combination of SMRT long-read sequencing and Illumina short-read sequencing. A total of 17,433 nonredundant full-length transcripts with average length of 2,271 bp and N50 length of 2,841 bp were obtained, among which 13,978 (80.18%) transcripts were annotated. Moreover, numerous transcript variants of various immune-related genes were identified, including pattern recognition receptors, antimicrobial peptides, heat shock proteins (HSPs), antioxidant enzymes and vital molecules in prophenoloxidase (proPO)-activating system. Based on the full-length transcriptome analysis, open reading frames (ORFs) of four C-type lectins (CTLs) were cloned, and tissue distributions showed that the four CTLs were ubiquitously expressed in all the tested tissues, and mainly expressed in hepatopancreas and gills. The transcription of the four CTLs significantly increased in several immune-related tissues (hemocytes, hepatopancreas and gills) of P. trituberculatus challenged with Vibrio alginolyticus and displayed different profiles. Moreover, the four CTLs displayed distinct bacterial binding and antibacterial activities. The recombinant protein PtCTL-1 (rPtCTL-1) and rPtCTL-3 displayed bacterial binding and antibacterial activities against all tested bacteria. rPtCTL-2 only showed bacterial binding and antibacterial activities against V. alginolyticus. No obvious bacterial binding or antibacterial activities for PtCTL-4 was observed against the tested bacteria. This study enriches the transcriptomic information on P. trituberculatus and provides new insights into the innate immune system of crustaceans. Additionally, our study provided candidates of antibiotic agents for the prevention and treatment of bacteriosis.

Comparative transcriptome profiling of horseshoe crab Tachypleus gigas hemocytes in response to lipopolysaccharides

Horseshoe crabs (HSCs) are living fossil species of marine arthropods with a long evolutionary history spanning approximately 500 million years. Their survival is helped by their innate immune system that comprises cellular and humoral immune components to protect them against invading pathogens. To help understand the genetic mechanisms involved, the present study utilised the Illumina HiSeq platform to perform transcriptomic analysis of hemocytes from the HSC, Tachypleus gigas, that were challenged with lipopolysaccharides (LPS). The high-throughput sequencing resulted in 352,077,208 and 386,749,136 raw reads corresponding to 282,490,910 and 305,709,830 high-quality mappable reads for the control and LPS-treated hemocyte samples, respectively. Based on the log-fold change of > 0.3 or < -0.3, 1338 genes were significantly upregulated and 215 genes were significantly downregulated following LPS stimulation. The differentially expressed genes (DEGs) were further identified to be associated with multiple pathways such as those related to immune defence, stress response, cytoskeleton function and signal transduction. This study provides insights into the underlying molecular and regulatory mechanisms in hemocytes exposed to LPS, which has relevance for the study of the immune response of HSCs to infection.

Comparative proteomics elucidates the dynamics of ovarian development in the Chinese mitten crab Eriocheir sinensis

Ovarian development is a complex physiological process for crustacean reproduction that is divided into the oogonium proliferation stage, endogenous vitellogenic stage, exogenous vitellogenic stage, and oocyte maturation stage. Proteomics analysis offers a feasible approach to reveal the proteins involved in the complex physiological processes of any organism. Therefore, this study performed a comparative proteomics analysis of the ovary and hepatopancreas at three key ovarian stages, including stages I (oogonium proliferation), II (endogenous vitellogenesis) and IV (exogenous vitellogenesis), of the Chinese mitten crab Eriocheir sinensis using a label-free quantitative approach. The results showed that a total of 2,224 proteins were identified, and some key proteins related to ovarian development and nutrition metabolism were differentially expressed. The 26 key proteins were mainly involved in the ubiquitin/proteasome pathway (UPP), cyclic AMP-protein kinase A (cAMP-PKA) signaling pathway, and mitogen-activated protein kinase (MAPK) signaling pathway during oogenesis. Fifteen differentially abundant proteins (DAPs) were found to participate in vitellogenesis and oocyte development, such as vitelline membrane outer layer protein 1 homolog, vitellogenin, vitellogenin receptor, heat shock 70 kDa protein cognate 3 and farnesyl pyrophosphate synthase. Forty-seven DAPs related to nutrition metabolism were identified, including the protein digestion, fatty acid metabolism, prostaglandin metabolism, lipid digestion and transportation, i.e. short-chain specific acyl-CoA dehydrogenase, acyl-CoA desaturase, fatty acid-binding protein, long-chain fatty acid CoA ligase 4, and hematopoietic prostaglandin D synthase. These results not only indicate proteins involved in ovarian development and nutrient deposition but also enhance the understanding of the regulatory pathways and physiological processes of crustacean ovarian development.

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.

Transcriptomic analyses provide insights into the adaptive responses to heat stress in the ark shells, Scapharca subcrenata

The ark shell, Scapharca subcrenata, is susceptible to high temperature which may lead to mass mortality in hot summers. Herein, we conducted the transcriptomic analyses of haemocytes in ark shells under thermal stress, to reveal the underlying molecular mechanisms of heat resistance in these animals. The results showed that a total of 7773, 11,500 and 13,046 unigenes were expressed differentially at 12, 24 and 48 h post thermal stress, respectively. The expression levels of key DEGs as revealed by RNA-seq were confirmed by quantitative real-time PCR. GO and KEGG enrichment analyses showed that the DEGs were mainly associated with apoptosis, NF-kappa B signaling pathway, TNF signaling pathway and RIG-I-like receptor signaling pathway. Among the DEGs, 40 were candidate heat stress response-related genes and 169 were identified to be involved in antioxidant defense, cell detoxification, protein metabolism and endoplasmic reticulum stress responses. It seemed that ark shells may adapt to short term thermal stress through regulation of protein metabolism, DNA replication and anti-apoptotic system. However, if the stress sustains, it may cause irreparable injury gradually in the animals due to oxygen limitation and metabolic dysregulation. Noteworthily, the expression of DEGs involved in protein biosynthesis and proteolysis was significantly elevated in ark shells under heat stress. These findings may provide preliminary insights into the molecular response of ark shells to acute thermal stress and lay the groundwork for marker-assisted selection of heat-resistant strains in S. subcrenata.

Effects of dietary icariin supplementation on the ovary development-related transcriptome of Chinese mitten crab (Eriocheir sinensis)

The Chinese mitten crab (Eriocheir sinensis) is an economically important aquaculture species in China, with distinct differences in ovarian maturation status between crabs fed with natural diets and artificial diets during the listing period, thus, leading to selling price differentiation. Our previous study showed that dietary supplementation with 100 mg/kg icariin can effectively promote ovarian development of E. sinensis. However, the internal molecular mechanism has not yet been elucidated because of a lack of comprehensive genome sequence information. We compared the ovary transcriptomes of E. sinensis fed with two diets containing 0 and 100 mg/kg ICA using the BGISEQ-500 platform. This yielded 12.54 Gb clean bases and 54,794 unigenes, 13,832 of which were found to be differentially expressed after icariin exposure. Twenty pathways closely related to gonadal development were selected through KEGG analysis. Seven differentially expressed genes relevant to vitellogenesis and oocyte maturation (serine/threonine-protein kinase mos-like, Eg2, cytoplasmic polyadenylation element-binding protein, cyclin B, vitellogenin 1, cathepsin D, and juvenile hormone esterase-like carboxylesterase 1) were validated by qRT-PCR, and four proteins (MEK1/2, ERK1/2, Cyclin B and Cdc2) associated with the progesterone mediated oocyte maturation pathway (i.e., MAPK/MPF pathway) were analyzed by western-blot. The results showed that icariin could promote the synthesis, processing and deposition of vitellogenin in oocytes, and that it also has the potential to promote oocyte maturation (resumption of Meiosis I) by altering the expression of the relevant genes and proteins.

Transcriptomic analysis and expression of C-type lectins in response to Vibrio parahaemolyticus challenge in Scapharca subcrenata

Little information is available on innate immune defense mechanisms of Scapharca subcrenata. C-type lectins (CTLs) are not only pattern recognition proteins that can bind pathogen-associated molecular patterns, but also crucial maternally-derived immune factors in mollusc egg. In this study, the comparative transcriptome analysis of Vibrio parahaemolyticus-infected and untreated hepatopancreas were performed to identify the key genes involved in maternal transfer of immunity. A total of 3514 and 9327 differentially expressed genes (DEGs) were identified at 6 and 48 h post challenge compared to control groups. Gene Ontology and Cluster of Orthologous Groups analysis showed that most DEGs were classified under regulation of signal transduction, regulation of the metabolic process of carbohydrates and secondary metabolites, while the processes of posttranscriptional modification and protein translation were inhibited manifestly. The DEGs were most enriched in pathways related to lysosome, phagosome and EMC-receptor interaction. Among the DEGs, 191 maternal immune-related genes that could provide developing embryos a better protection against pathogen infection were identified according to previous studies. Additionally, five CTLs (designated as SsCTL1-5) identified from the DEGs were cloned, and their expression patterns in different tissues and post immune stimulation were analyzed. These findings would be beneficial for understanding the innate immune defense mechanisms of S. subcrenata.

Uncloaking lncRNA-meditated gene expression as a potential regulator of CMS in cotton (Gossypium hirsutum L.)

Cytoplasmic male sterility is a well-proven mechanism for cotton hybrid production. Long non-coding RNAs belong to a class of transcriptional regulators that function in multiple biological processes. The cDNA libraries from the flower buds of the cotton CGMS, it's restorer (Rf) and maintainer lines were sequenced using high throughput NGS technique. A total of 1531 lncRNAs showed significant differential expression patterns between these three lines. Functional analysis of the co-expression network of lncRNA-mRNA using gene ontology vouchsafes that, lncRNAs play a crucial role in cytoplasmic male sterility and fertility restoration through pollen development, INO80 complex, development of anther wall tapetum, chromatin remodeling, and histone modification. Additionally, 94 lncRNAs were identified as putative precursors of 49 miRNAs. qRT-PCR affirms the concordance of expression pattern to RNA-seq data. These findings divulge the lncRNA driven miRNA-mediated regulation of gene expression profiling superintended for a better understanding of the CMS mechanisms of cotton.

Evolution of digestive enzyme genes associated with dietary diversity of crabs

Crabs feed on a wide range of items and display diverse feeding strategies. The primary objective of this study was to investigate 10 digestive enzyme genes in representative crabs to provide insights into the genetic basis of feeding habits among crab functional groups. Crabs were classified into three groups based on their feeding habits: herbivores (HV), omnivores (OV), and carnivores (CV). To test whether crabs' feeding adaptations matched adaptive evolution of digestive enzyme genes, we examined the 10 digestive enzyme genes of 12 crab species based on hepatopancreas transcriptome data. Each of the digestive enzyme genes was compared to orthologous sequences using both nucleotide- (i.e., PAML and Datamonkey) and protein-level (i.e., TreeSAAP) approaches. Positive selection genes were detected in HV crabs (AMYA, APN, and MGAM) and CV crabs (APN, CPB, PNLIP, RISC, TRY, and XPD). Additionally, a series of positive selection sites were localized in important functional regions of these digestive enzyme genes. This is the first study to characterize the molecular basis of crabs' digestive enzyme genes based on functional feeding group. Our data suggest that HV crabs have evolved an enhanced digestion capacity for carbohydrates, and CV crabs have acquired digestion capacity for proteins and lipids.

Comparative transcriptome analysis of the gills of Procambarus clarkii provides novel insights into the immune-related mechanism of copper stress tolerance

The red-swamp crayfish (Procambarus clarkii) is the most important economic shrimp species in China, and is an important model crustacean organism in many fields of research. In crustaceans, gills interface directly with the ambient environment and thus play a vital role in the toxicology. In the context of increasing environmental heavy metal pollution, the relationship between copper (Cu²⁺) stress and the immune response of P. clarkii has recently received considerable attention. However, impact of Cu²⁺ on the crayfish immune system is still not fully understood. In this study, we used Illumina sequencing technology to perform a transcriptome analysis of the gills of P. clarkii after 24 h of Cu²⁺ treatment. A total of 37,226,812 unigenes were assembled, and 1943 unigenes were significantly differentially expressed between the control and Cu²⁺ treatment groups. Functional categorization of differentially expressed genes (DEGs) revealed that genes related to antioxidant activity, detoxication, metabolic processes, biosynthetic processes, and immune system processes were differentially regulated during Cu²⁺ stress. In addition, DEGs in the immune system were classified as being related to the MAPK signaling pathway, purine metabolism, Toll and Imd signaling pathway, PI3K-Akt signaling pathway and Hippo signaling pathway. Five genes (CuZnSOD, CAT, IDH1, PHYH and DECR2) were significantly up-regulated in the peroxisome pathway, which plays an important role in reacting to oxidative stress. Importantly, qRT-PCR validation of the results for seven genes chosen at random (NDK, ATP6L, ATP5C1, RPS14, RPL22e, CTSF and HSP90A) confirmed the Illumina sequencing results. This study provides a valuable starting point for further studies to elucidate the molecular basis of the immune system's response to Cu²⁺ stress in crayfish.

Chromosome-level genome assembly reveals the unique genome evolution of the swimming crab (Portunus trituberculatus)

BACKGROUND

The swimming crab, Portunus trituberculatus, is an important commercial species in China and is widely distributed in the coastal waters of Asia-Pacific countries. Despite increasing interest in swimming crab research, a high-quality chromosome-level genome is still lacking.

FINDINGS

Here, we assembled the first chromosome-level reference genome of P. trituberculatus by combining the short reads, Nanopore long reads, and Hi-C data. The genome assembly size was 1.00 Gb with a contig N50 length of 4.12 Mb. In addition, BUSCO assessment indicated that 94.7% of core eukaryotic genes were present in the genome assembly. Approximately 54.52% of the genome was identified as repetitive sequences, with a total of 16,796 annotated protein-coding genes. In addition, we anchored contigs into chromosomes and identified 50 chromosomes with an N50 length of 21.80 Mb by Hi-C technology.

CONCLUSIONS

We anticipate that this chromosome-level assembly of the P. trituberculatus genome will not only promote study of basic development and evolution but also provide important resources for swimming crab reproduction.

Comparative transcriptome reveals the potential modulation mechanisms of estradiol affecting ovarian development of female Portunus trituberculatus

Estradiol is an important sex steroid hormone that is involved in the regulation of crustacean ovarian development. However, the molecular regulatory mechanisms of estradiol on ovarian development are largely unknown. This study performed transcriptome sequencing of ovary, hepatopancreas, brain ganglion, eyestalk, and mandibular organ of crabs after estradiol treatment (0.1μg g^-1 crab weight). A total of 23, 806 genes were annotated, and 316, 1300, 669, 142, 383 genes were expressed differently in ovary, hepatopancreas, brain ganglion, eyestalk, and mandibular organ respectively. Differentially expressed gene enrichment analysis revealed several crucial pathways including protein digestion and absorption, pancreatic secretion, insect hormone biosynthesis, drug metabolism-cytochrome P450 and signal transduction pathway. Through this study, some key genes in correlation with the ovarian development and nutrition metabolism were significantly affected by estradiol, such as vitelline membrane outer layer 1-like protein, heat shock protein 70, Wnt5, JHE-like carboxylesterase 1, cytochrome P302a1, crustacean hyperglycemic hormone, neuropeptide F2, trypsin, carboxypeptidase B, pancreatic triacylglycerol lipase-like, and lipid storage droplet protein. Moreover, RT-qPCR validation demonstrated that expression of transcripts related to ovarian development (vitelline membrane outer layer 1-like protein and cytochrome P302a1) and nutrition metabolism (trypsin, glucose dehydrogenase and lipid storage droplet protein) were significantly affected by estradiol treatment. This study not only has identified relevant genes and several pathways that are involved in estradiol regulation on ovarian development of P. trituberculatus, but also provided new insight into the understanding of the molecular function mechanisms of estradiol in crustacean.

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.

Integrated mRNA and miRNA expression profile analyses reveal the potential roles of sex-biased miRNA-mRNA pairs in gonad tissues of the Chinese concave-eared torrent frog (Odorrana tormota)

The Chinese concave-eared torrent frog (Odorrana tormota) is typically sexually dimorphic. Females are significantly less common than males in the wild. Until now, the molecular mechanisms of reproduction and sex differentiation of frogs remain unclear. Here, we integrated mRNA and microRNA (miRNA) expression profiles to reveal the molecular mechanisms of reproduction and sex differentiation in O. tormota. We identified 234 differentially expressed miRNAs (DEMs) and 18,551 differentially expressed transcripts. Of these, 12,053 mRNAs and 64 miRNAs were upregulated in testes, and 6,498 mRNAs and 170 miRNAs were upregulated in ovaries. Integrated analysis of the miRNA and mRNA expression profiles predicted 75,602 potential miRNA-mRNA interaction sites, with 42,065 negative miRNA-mRNA interactions. We found 36 differentially expressed genes (DEGs) related to reproduction and sex differentiation, of which 15 DEGs formed 92 negative miRNA-mRNA interactions with 34 known DEMs. Thus, miRNAs may play other important roles in O. tormota. Furthermore, Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analyses showed reproductive-related processes, such as the gonadotropinreleasing hormone signaling pathway and ovarian steroidogenesis. Based on functional annotation and the literature, the retinoic acid signaling pathway, the SOX9-AMH pathway, and the process of spermatogenesis may be involved in the molecular mechanisms of reproduction and sex differentiation in O. tormota, and may be regulated by miRNAs. The miRNA-mRNA pairs described may provide further understanding of the regulatory mechanisms associated with reproduction and sex differentiation, and the molecular mechanism of reproduction in O. tormota.

The crustacean ecdysone cassette: A gatekeeper for molt and metamorphosis

Arthropods have long been utilized as models to explore molecular function, and the findings derived from them can be applied throughout metazoa, including as a basis for medical research. This has led to the adoption of many representative insect models beyond Drosophila, as each lends its own unique perspective to questions in endocrinology and genetics. However, non-insect arthropods are yet to be realised for the potential insight they may provide in such studies. The Crustacea are among the most ancient arthropods from which insects descended, comprising a huge variety of life histories and ecological roles. Of the events in a typical crustacean development, metamorphosis is perhaps the most ubiquitous, challenging and highly studied. Despite this, our knowledge of the endocrinology which underpins metamorphosis is rudimentary at best; although several key molecules have been identified and studied in depth, the link between them is quite nebulous and leans heavily on well-explored insect models, which diverged from the Pancrustacea over 450 million years ago. As omics technologies become increasingly accessible, they bring the prospect of explorative molecular research which will allow us to uncover components and pathways unique to crustaceans. This review reconciles known components of crustacean metamorphosis and reflects on our findings in insects to outline a future search space, with focus given to the ecdysone cascade. To expand our knowledge of this ubiquitous endocrine system not only aids in our understanding of crustacean metamorphosis, but also provides a deeper insight into the adaptive capacity of arthropods throughout evolution.

Transcriptome sequencing and molecular markers discovery in the gonads of Portunus sanguinolentus

Crab culture has gained prominence in the last decade due to the large global market demand for live crabs and crab products. Portunus sanguinolentus is one of the economically important crab species in the Indo-Pacific region, with distinct differences in growth and size between male and female crabs, thus, leading to huge difference in their market values. The culture of P. sanguinolentus is still in its infancy, with crab supplies heavily dependent on wild catch. In order to unravel the molecular differences between male and female crabs, we generated a comprehensive transcriptomic dataset for P. sanguinolentus by sequencing the gonads of both sexes using the Illumina Hiseq 2500 system. Transcriptomes were assembled using Trinity de novo assembly followed by annotation. This transcriptomic data set for P. sanguinolentus would serve as an important reference data for genomic and genetic studies in this crab and related species.