Transcriptional changes revealed water acidification leads to the immune response and ovary maturation delay in the Chinese mitten crab Eriocheir sinensis.

Nowadays, due to increasing carbon dioxide released, water acidification poses a series of serious impacts on aquatic organisms. To evaluate the effects of water acidification on crustaceans, we focused on the Chinese mitten crab Eriocheir sinensis, which is a spawning migration and farmed species in China. Based on histological and oocyte transparent liquid observation, we found that the acidified environment significantly delayed the ovarian maturation of E. sinensis. Moreover, RNA-seq was applied to obtain gene expression profile from the crab's gills and ovaries in response to acidified environment. Compared with control groups, a total of 5471 differentially expressed genes (DEGs) were identified in acidified gills and 485 DEGs were identified in acidified ovaries. Enrichment analysis indicated that some pathways also responded to the acidified environment, such as PI3K-Akt signaling pathway, Chemokine signaling pathway, apoptosis, and toll-like receptor signaling pathway. Subsequently, some DEGs involved in immune response (ALF, Cathepsin A, HSP70, HSP90, and catalase) and ovarian maturation (Cyclin B, Fem-1a, Fem-1b, and Fem-1c) were selected to further validate the influence of water acidification on gene expression by qRT-PCR. The results showed that the expression level of immune-related genes was significantly increased to response to the water acidification, while the ovarian maturation-related genes were significantly decreased. Overall, our data suggested that E. sinensis was sensitive to the reduced pH. This comparative transcriptome also provides valuable molecular information on the mechanisms of the crustaceans responding to acidified environment.

Identification of a novel germ cell marker MnTdrd from the oriental river prawn Macrobrachium nipponense.

Germ cell-specific genes play an important role in establishing the reproductive system in sexual organisms and have been used as valuable markers for studying gametogenesis and sex differentiation. Previously, we isolated a vasa transcript as a germ cell marker to trace the origin and migration of germ cells in the oriental river prawn Macrobrachium nipponense. Here, we identified a new germ cell-specific marker MnTdrd RNA and assessed its temporal and spatial expression during oogenesis and embryogenesis. MnTdrd transcripts were expressed in high abundance in unfertilized eggs and embryos at cleavage stage and then dropped significantly during late embryogenesis, suggesting that MnTdrd mRNA is maternally inherited. In situ hybridization of ovarian tissue showed that MnTdrd mRNA was initially present in the cytoplasm of previtellogenic oocyte and localized to the perinuclear region as the accumulation of yolk in vitellogenic oocyte. Whole-mount in situ hybridization of embryos showed that MnTdrd-positive signals were only localized in one blastomere until 16-cell stage. In the blastula, there were approximately 16 MnTdrd-positive blastomeres. During embryonized-zoea stage, the MnTdrd-positive cells aggregated as a cluster and migrated to the genital rudiment which would develop into primordial germ cells (PGCs). The localized expression pattern of MnTdrd transcripts resembled that of the previously identified germ cell marker vasa, supporting the preformation mode of germ cell specification. Therefore, we concluded that MnTdrd, together with vasa, is a component of the germ plasm and might have critical roles in germ cell formation and differentiation in the prawn. Thus, MnTdrd can be used as a novel germ cell marker to trace the origin and migration of germ cells.

Molecular characterization of ovary-specific gene Mrfem-1 and siRNA-mediated regulation on targeting Mrfem-1 in the giant freshwater prawn, Macrobrachium rosenbergii.

Characterized by ankyrin repeat motifs, the feminization-1 (fem-1) gene plays an essential role in sex determination/differentiation in Caenorhabditis elegans. However, there are only a few reports on fem-1 in crustaceans. In this study, a fem-1 gene (Mrfem-1) was first isolated from the giant freshwater prawn Macrobrachium rosenbergii. The full-length cDNA of Mrfem-1 was 2607 bp long, containing an open reading frame encoding 615 amino acids, and presenting eight ankyrin repeats. The full-length cDNA has been submitted to GenBank with the accession no. MT160093. According to the RT-PCR results, Mrfem-1 was exclusively expressed in the ovary. The expression level of Mrfem-1 had increased with ovarian maturation and reached the highest peak at vitellogenic stage. In situ hybridization results showed that positive signals were concentrated in the cytoplasm of previtellogenic stage, and scattered in the cytoplasm and follicular cells at vitellogenic stage, suggesting that Mrfem-1 might be associated with ovarian maturation. Moreover, two effective siRNAs targeting Mrfem-1 were found and their effectiveness verified in vitro. These results on Mrfem-1 will help us to better understand the fem family and provide a new resource for subsequent investigations of siRNA-mediated regulation on ovarian development in M. rosenbergii.

Alternative splice variants and differential relative abundance patterns of vasa mRNAs during gonadal development in the Chinese mitten crab Eriocheir sinensis.

Gonadal development usually involves alternative splicing of sex-related genes. Vasa, a highly conserved ATP-dependent RNA helicase present mainly in germ cells, has an important function in gonadal development. As an important sex-related gene, recent evidence indicates that different splice variants of vasa exist in many species. In this study, there was identification of two types of vasa splice variants in the Chinese mitten crab Eriocheir sinensis, termed Esvasa-l and Esvasa-s, respectively. Furthermore, splice variants of Esvasa-s were sub-divided into Esvasa-s1, Esvasa-s2, Esvasa-s3, Esvasa-s4, and Esvasa-s5, based on differing numbers of TGG repeats. Results from genomic structure analyses indicated that these forms are alternatively spliced transcripts from a single vasa gene. Results from tissue distribution assessments indicate the vasa splice variants were exclusively expressed in the gonads of male and female adult crabs. In situ hybridization results indicate Esvasa mRNA was mainly present in the cytoplasm of previtellogenic oocytes. As oocyte size increased, relative abundance of Esvasa mRNA decreased and became distributed near the cellular membrane. The Esvasa mRNA was not detectable in mature oocytes. In testis, Esvasa mRNA was detected in spermatids and spermatozoa, but not in spermatogonia and spermatocytes. Notably, results from qPCR analysis of Esvasa-l and Esvasa-s indicate there are different relative proportions during gametogenesis, implying that splice variants of the Esvasa gene may have different biological functions during crab gonadal development.

Identification and profiling of microRNAs during gonadal development in the giant freshwater prawn Macrobrachium rosenbergii.

As post-transcriptional regulators, microRNAs (miRNAs) play an important role in growth and reproductive processes. So far, there is limited information regarding crustacean miRNAs. To explore the potential role of miRNAs in the gonadal development of the prawn Macrobrachium rosenbergii, we constructed seven small RNA libraries from ovarian and testicular tissues at various stages using somatic tissue as the control. A total of 1,954 known and 129 novel miRNAs were retrieved. By comparing differentially expressed miRNAs (DEMs) between testes and ovaries, forty-one miRNAs were identified with sex-biased expression patterns, including 17 ovary-biased and 24 testis-biased patterns. Furthermore, the putative target genes of the sex-biased miRNAs, such as cyclin L1, mitogen-activated protein kinase 7 (MAPK 7), heat shock protein (HSP), and zinc finger protein, were significantly enriched in many reproduction-related pathways including the Gonadotropin-releasing hormone (GnRH) pathway, glycolysis, gluconeogenesis pathway, ovarian steroidogenesis, estrogen signaling pathway, MAPK pathway, Wnt pathway, and insulin signaling pathway, implicating potential regulatory roles of miRNAs in reproduction. These data aid in the further investigation of the mechanism of gonadal development and reproductive regulation mediated by miRNA in M. rosenbergii.

Construction of a Genomic Bacterial Artificial Chromosome (BAC) Library for the Prawn Macrobrachium rosenbergii and Initial Analysis of ZW Chromosome-Derived BAC Inserts.

Knowledge on sex determination has proven valuable for commercial production of the prawn Macrobrachium rosenbergii due to sex dimorphism of the male and female individuals. Previous studies indicated that prawn sex is determined by a ZW-ZZ chromosomal system, but no genomic information is available for the sex chromosome. Herein, we constructed a genomic bacterial artificial chromosome (BAC) library and identified the ZW-derived BAC clones for initial analysis of the sex chromosomal DNA sequence. The arrayed BAC library contains 200,448 clones with average insert size of 115.4 kb, corresponding to ∼ 4× coverage of the estimated 5.38 Gb genome. Based on a short female-specific marker, a Z- and a W-fragment were retrieved with the genomic walking method. Screening the BAC library using a ZW-specific marker as probe resulted in 12 positive clones. From these, a Z-derived (P331M17) and a W-derived (P122G2) BAC clones were randomly selected and sequenced by PacBio method. We report the construction of a large insert, deep-coverage, and high-quality BAC library for M. rosenbergii that provides a useful resource for positional cloning of target genes, genomic organization, and comparative genomics analysis. Our study not only confirmed the ZW/ZZ system but also discovered sex-linked genes on ZW chromosomes for the first time, contributing to a comprehensive understanding of the genomic structure of sex chromosomes in M. rosenbergii.

Molecular cloning and characterization of a gonadotropin-releasing hormone receptor homolog in the Chinese mitten crab, Eriocheir sinensis.

As an essential mediator in the Gonadotropin-releasing hormone (GnRH) signaling pathway, GnRH receptor (GnRHR) coupled to GnRH, plays an important role in activating the downstream pathway after stimulating a series of cascades to regulate reproduction. To detect the existence of GnRHR and potential GnRH signaling pathway, we cloned and characterized GnRHR in the Chinese mitten crab, Eriocheir sinensis (named EsGnRHR). The full-length EsGnRHR cDNA is 2038 bp in length, including an open reading frame (ORF) of 1566 bp, a 57 bp 5'-untranslated region (5'-UTR) and a 415 bp 3'-UTR. Prediction of transmembrane domains in protein sequence revealed that the EsGnRHR protein contained seven hydrophobic transmembrane regions (TMs). Reverse transcription PCR revealed that EsGnRHR was mainly expressed in the thoracic nerve group and ovary, and weakly distributed in the testis and brain. In situ hybridization further demonstrated that EsGnRHR mRNA was localized at the protocerebrum and deutocerebrum. In the ovary and testis, the hybridization signal was dominantly at the earlier developmental stages. The signal was mainly localized in the cytoplasm cell in the ovary, and in the epithelium cell in the testis. During the different stages of gonadal development, EsGnRHR displayed increasing trends in both female and male when analyzed by quantitative real-time PCR, suggesting that EsGnRHR was involved in controlling gonadal development. Our study provides important information for further research on the molecular mechanisms underlying crab development.

Identification of putative regulatory region of insulin-like androgenic gland hormone gene (IAG) in the prawn Macrobrachium nipponense and proteins that interact with IAG by using yeast two-hybrid system.

Insulin-like androgenic gland hormone gene (IAG) is a sex regulator specifically expressed in male crustaceans, controlling the male sexual differentiation, spermatogenesis and reproductive strategy. Our previous study reported the cloning and characterization of the prawn Macrobrachium nipponense IAG (MnIAG). In this study, we further identified a 2214-bp MnIAG 5'-flanking region, and analyzed its transcription factor binding sites and transcriptional activity. The results showed that there were two potential promoter core sequences, three TATA boxes and one CAAT box existing in the MnIAG 5'-flanking region as well as many potential transcription factor binding sites, such as SRY, Sox-5, GATA-1, etc. Notably, the transcriptional activity was weak in this region, and a negative regulatory region was found in -604 to -231bp. In addition, we constructed M. nipponense yeast libraries and identified proteins interacting with the MnIAG protein by yeast two hybridization assay. The yeast two-hybrid screening yielded ten positive clones, of which five were annotated by NCBI database, namely heat shock protein 21, NADH dehydrogenase, zinc finger protein, beta-N-acetylglucosaminidase and a hypothetical protein. The identification of MnIAG putative regulatory region and proteins that interact with IAG will facilitate our understanding of the regulatory role of MnIAG and provide a foundation for deep insight into the prawn sex differentiation mechanism and signaling transduction pathways.

Molecular characterization of a novel ovary-specific gene fem-1 homolog from the oriental river prawn, Macrobrachium nipponense.

The feminization-1 (fem-1) gene is characterized by one of the most common protein-protein interaction motifs, ankyrin repeat motifs, displays many expression patterns in vertebrates and invertebrates, and plays an essential role in the sex-determination/differentiation pathway in Caenorhabditis elegans. In this study, a fem-1 homolog, designated as Mnfem-1, was first cloned from the oriental river prawn Macrobrachium nipponense. The prawn Mnfem-1 gene consists of six exons and five introns. The full-length cDNA (2603bp) of Mnfem-1 contains an open reading frame (ORF) encoding a protein of 622 amino acids. The Mnfem-1 RNA and protein are exclusively expressed in the ovary in adult prawns as revealed by RT-PCR and immunofluorescence analysis, respectively. In situ hybridization results showed that strong positive signals were concentrated at the edge of the previtellogenic and vitellogenic oocyte. During embryogenesis, Mnfem-1 is highly expressed in both unfertilized eggs and embryos at cleavage stage and thereafter dropped to a low level from blastula to zoea, indicating that the Mnfem-1 in early embryos is maternal. After hatching, the Mnfem-1 expression significantly increased in the larvae at length of 2cm, an important stage of sex differentiation. Yeast two hybridization results showed that the Mnfem-1 protein can be potentially interactive with cathepsin L and proteins containing the domains of insulinase, ankyrin or ubiquitin. Our results suggested that Mnfem-1 could have roles in prawn ovarian development and sex determination/differentiation.

Molecular cloning and characterization of perlucin from the freshwater pearl mussel, Hyriopsis cumingii.

Perlucin is an important functional protein that regulates shell and pearl formation. In this study, we cloned the perlucin gene from the freshwater pearl mussel Hyriopsis cumingii, designated as Hcperlucin. The full-length cDNA transcribed from the Hcperlucin gene was 1460 bp long, encoding a putative signal peptide of 20 amino acids and a mature protein of 141 amino acids. The mature Hcperlucin peptide contained six conserved cysteine residues and a carbohydrate recognition domain, similar to other members of the C-type lectin families. In addition, a "QPS" and an invariant "WND" motif near the C-terminal region were also found, which are extremely important for polysaccharide recognition and calcium binding of lectins. The mRNA of Hcperlucin was constitutively expressed in all tested H. cumingii tissues, with the highest expression levels observed in the mantle, adductor, gill and hemocytes. In situ hybridization was used to detect the presence of Hcperlucin mRNA in the mantle, and the result showed that the mRNA was specifically expressed in the epithelial cells of the dorsal mantle pallial, an area known to express genes involved in the biosynthesis of the nacreous layer of the shell. The significant Hcperlucin mRNA expression was detected on day 14 post shell damage and implantation, suggesting that the Hcperlucin might be an important gene in shell nacreous layer and pearl formation. The change of perlucin expression in pearl sac also confirmed that the mantle transplantation results in a new expression pattern of perlucin genes in pearl sac cells that are required for pearl biomineralization. These findings could help better understanding the function of perlucin in the shell and pearl formation.

Molecular characterization and expression analysis of an insulin-like gene from the androgenic gland of the oriental river prawn, Macrobrachium nipponense.

The androgenic gland (AG), a male-specific endocrine organ in crustacean, is responsible for the maintenance of male characteristics and gender differentiation. In this study, an AG-specific gene, the Macrobrachium nipponesne insulin-like androgenic gland factor (MnIAG) was isolated from a transcriptome library of M. nipponesne and its full-length cDNA sequences were obtained by RACE method. The cDNA was 1,547 bp in length and encoded a precursor protein of 175 amino acids. The deduced precursor protein consisted of a signal peptide, B chain, C peptide and an A chain, which exhibited the same structural organization as that of previously identified insulin-like androgenic gland in crustacean. The mature peptide of the MnIAG owned two additional conserved cysteine residues, which were also found in the Palaemonidae species reported. Results of the tissue distribution and in situ hybridization showed the MnIAG expressed exclusively in androgenic gland. The quantitative RT-PCR results demonstrated that the MnIAG transcript was present at blastula stage and later developmental stages with low levels, which suggested that the primordial cells of the AG might form at these stages.

The complete mitochondrial genome of Channa argus, Channa maculata and hybrid snakehead fish [Channa maculata (♀) × Channa argus (♂)].

We sequenced and characterized the complete mitochondrial genome of Channa argus, Channa maculata and their hybrid [C. maculata (♀) and C. argus (♂)]. All the three mitochondrial genomes contained the typical complement of 13 protein-coding genes, 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs) and 1 control region. The entire mitochondrial DNA (mtDNA) molecule of C. maculata was 16,559 bp long while the complete mtDNA molecule of C. argus and hybrid snakehead fish was 16,558 bp long. This is the first report on the complete mitogenome sequence of C. maculata and hybrid snakehead fish.

[Genetic structure of wild Macrobrachium nipponense populations in Taihu Lake based on microsatellite analysis].

By using eight highly polymorphic microsatellite DNA loci, this paper analyzed the genetic structure of wild Macrobrachium nipponense populations in Taihu Lake. For the 15 M. nipponense populations in the Lake, there were at least three of the loci presenting heterozygosity deficiency and obvious deviation from Hardy-Weinberg equilibrium after Bonferroni correction. The observed heterozygosity values of the 15 populations were all above 0. 683, displaying a high genetic diversity, but the diversity varied obviously with site. For example, the genetic diversity of the eastern and southern populations at Dukou and Luxiang was higher than that of the western and northern populations at Huazhuang and Yangzhu. For the 15 populations, parts of the loci showed heterozygote excess and departure from mutation-drift equilibrium, suggesting that the population structure had experienced bottleneck effect and the population amount had declined. The AMOVA analysis across all the populations and loci showed that the genetic divergence among the 15 populations was at a lower level (F(ST) = 0.011 ). 98.9% of the genetic variation came from intra-population, and 1.1% came from inter-population, suggesting that all the M. nipponense populations in the Lake could be protected and managed as a single unit in genetic resource. However, the genetic distance between Huazhuang and Wutangmen populations reached 0.206, being close to the delimitation of species identification. Further studies would be needed for the sustainable utilization of the genetic resource of M. nipponense in Taihu Lake.

[Microsatellite analysis of genetic variation of the oriental river prawn Macrobrachium nipponense in Qiantang River].

Genetic diversity and genetic structure of 7 wild stocks of oriental river prawn Macrobrachium nipponense in Qiantang River, i.e. Wen-yan, Fu-yang, Chang-kou, Tong-lu, Xin-an-jiang, She-xian and Xiu-ning, were investigated using 10 microsatellite DNA markers. The result showed that all the 10 loci were highly polymorphic. There was a trend that the level of genetic diversity of wild stocks in downstream and midstream were higher than the upstream ones'. Sign test and Wilcoxon sign rank test results showed that the stocks in Qiantang River had no bottleneck effect, and the number of stocks had not declined recently. F(ST) among stocks ranged from 0.0201 to 0.1069. Analysis of molecular variance (AMOVA) revealed that a higher portion (93.48%) of variations existed within individuals, while lower portion (6.52%) existed among stocks. F(ST) and AMOVA analysis across all stocks and loci indicated the medium level of divergence among the stocks. The NJ clustering tree based on D(A) genetic distance demonstrated that the stocks of adjacent geographical position clustered together. 413 individuals obtained from six wild stocks could be divided into two potential populations based on the genetic structure. This study demonstrated that genetic diversity and genetic structure of M. nipponense stocks were relevant to geographical position where they survived.