Histopathological observation and comparative transcriptome analysis reveal immune response mechanisms to Aeromonas dhakensis infection in Macrobrachium rosenbergii.

The Macrobrachium rosenbergii industry is threatened by various Aeromonas, resulting in high mortality of adult prawns. However, there are few studies on the immune response of M. rosenbergii infected with Aeromonas dhakensis. In this study, we observed the hepatopancreas and gills histopathologically, performed a comparative transcriptome analysis of the hepatopancreas, and analyzed the candidate gene expression of immune-related genes in the hemolymph, hepatopancreas, and gills of M. rosenbergii that had been infected with A. dhakensis. Histopathology revealed the hepatopancreas was successively inflamed, followed by cellular vacuolation, lumen deformation, and finally tissue erosion; partial and severe inflammation of the gills occurred successively, and eventually the gill tissue atrophy and the gill filaments detached from the gill arch. Transcriptome analysis showed that a total of 77,742 unigenes and 8664 differentially expressed genes (DEGs), and the immune-related DEGs were mainly enriched in lysosome and phagosome pathways. In addition, 4 immune-related candidate genes (RhoA, CASP9, PKC, and DSCIGN) based on KEGG and PPI analysis were monitored at 6, 12, and 24h post injection (hpi) in hepatopancreas, hemolymph and gills. Their spatio-temporal expression results indicated that A. dhakensis have activated the immune system of M. rosenbergii. The present study may provide new information on the complex immune mechanism of M. rosenbergii.

Genomic evidence for the population genetic differentiation of Misgurnus anguillicaudatus in the Yangtze River basin of China.

Misgurnus anguillicaudatus, an important aquatic species, is mainly distributed in the Yangtze River basin. To reveal the population genetic structure of M. anguillicaudatus distributed in the Yangtze River basin, genotyping by sequencing (GBS) technique was employed to detect the genome wide genetic variations of M. anguillicaudatus. A total of 30.03 Gb raw data were yielded from 70 samples collected from 15 geographic sites located in the Yangtze River basin. Subsequently, 2092 high quality SNPs were genotyped across these samples and used for a series of genetic analysis. The results of genetic analysis showed that high levels of genetic diversity were observed and the populations from upper reaches (UR) were significantly differentiated from the middle and lower reaches (MLR) of Yangtze River basin. Meanwhile, no significant isolation by distance was detected among the populations. Ecological factors (e.g. complicated topography and climatic environment) and anthropogenic factors (e.g. aquaculture and agriculture cultivation) might account for the genetic disconnectivity between UR and MLR populations. This study provided valuable genetic data for the future breeding program and also for the conversation and scientific utilization of those abundant genetic resources stored in the Yangtze River basin.

Evolution and phylogeography analysis of diploid and polyploid Misgurnus anguillicaudatus populations across China.

The origin and evolution of polyploid organisms have been extensively studied in plants, but this topic remains only partially understood in vertebrates, where polyploidy is relatively rare. In this study, we used Misgurnus anguillicaudatus, a fish that comprises five ploidy levels in nature, as a model animal to improve our understanding of biogeographic history and evolution of polyploid vertebrates. After collecting samples from different geographical populations in China, their ploidy levels were determined using flow cytometry. Two mitochondrial markers ( cytochrome b and control region) were then used for phylogeographic analyses to unravel the possible origins of diploids and tetraploids in China. The results showed that diploids have wider geographical distribution than tetraploids and triploids. There was no clear allopatric geographical range or boundary to divide diploid and polyploid populations. Rather, the analysis of mitochondrial DNA sequences indicated that tetraploids were autopolyploids, with lower genetic diversity than diploids. This suggests that tetraploids originated from sympatric diploids via multiple independent polyploidization events. Genetic structure patterns were similar between diploids and tetraploids, whereas complex genetic differentiation was found among different regions. The potential origin of M. anguillicaudatus was deduced to be in the Pearl River basin, which exhibited the highest nucleotide diversity and genetic differentiation. These findings provide insights into the evolution of polyploidy in vertebrates.

Gene expression patterns indicate that a high-fat-high-carbohydrate diet causes mitochondrial dysfunction in fish.

Expensive and unsustainable fishmeal is increasingly being replaced with cheaper lipids and carbohydrates as sources of energy in aquaculture. Although it is known that the excess of lipids and carbohydrates has negative effects on nutrient utilization, growth, metabolic homeostasis, and health of fish, our current understanding of mechanisms behind these effects is limited. To improve the understanding of diet-induced metabolic disorders (both in fish and other vertebrates), we conducted an eight-week high-fat-high-carbohydrate diet feeding trial on blunt snout bream (Megalobrama amblycephala), and studied gene expression changes (transcriptome and qPCR) in the liver. Disproportionately large numbers of differentially expressed genes were associated with mitochondrial metabolism, neurodegenerative diseases (Alzheimer's, Huntington's, and Parkinson's), and functional categories indicative of liver dysfunction. A high-fat-high-carbohydrate diet may have caused mitochondrial dysfunction, and possibly downregulated the mitochondrial biogenesis in the liver. While the relationship between diet and neurodegenerative disorders is well-established in mammals, this is the first report of this connection in fish. We propose that fishes should be further explored as a potentially promising model to study the mechanisms of diet-associated neurodegenerative disorders in humans.

Selection shapes the patterns of codon usage in three closely related species of genus Misgurnus.

Neutrality plots revealed that selection probably dominates codon bias, whereas mutation plays only a minor role, in shaping the codon bias in three loaches, Misgurnus anguillicaudatus, M. mohoity, and M. bipartitus. These three species also clearly showed similar tendencies in the preferential usage of codons. Nineteen, nine, and 14 preferred codon pairs and 179, 182, and 173 avoided codon pairs were also detected in M. anguillicaudatus, M. bipartitus, and M. mohoity, respectively, and the most frequently avoided type of cP3-cA1 dinucleotide in these species was nnUAnn. The expression-linked patterns of codon usage revealed that higher expression was associated with higher GC3, lower ENC, and a smaller proportion of amino acids with high size/complexity (S/C) scores in these three species. These results elucidate selectively driven codon bias in Misgurnus species, and reveal the potential importance of expression-mediated selection in shaping the genome evolution of fish.

Genome survey sequencing of red swamp crayfish Procambarus clarkii.

Red swamp crayfish, Procambarus clarkii, presently is an important aquatic commercial species in China. The crayfish is a hot area of research focus, and its genetic improvement is quite urgent for the crayfish aquaculture in China. However, the knowledge of its genomic landscape is limited. In this study, a survey of P. clarkii genome was investigated based on Illumina's Solexa sequencing platform. Meanwhile, its genome size was estimated using flow cytometry. Interestingly, the genome size estimated is about 8.50 Gb by flow cytometry and 1.86 Gb with genome survey sequencing. Based on the assembled genome sequences, total of 136,962 genes and 152,268 exons were predicted, and the predicted genes ranged from 150 to 12,807 bp in length. The survey sequences could help accelerate the progress of gene discovery involved in genetic diversity and evolutionary analysis, even though it could not successfully applied for estimation of P. clarkii genome size.

Comprehensive Transcriptome Analysis Provides Evidence of Local Thermal Adaptation in Three Loaches (Genus: Misgurnus).

The geographic distribution of three Misgurnus species, M. anguillicaudatus, M. bipartitus, and M. mohoity, displays a specific pattern in China, coincident with temperature zones. In this study, we sequenced the transcriptomes of these three species and used the sequences to investigate the lineage-specific adaptations within the genus Misgurnus. In total, 51 orphan genes (19 in M. anguillicaudatus, 18 in M. bipartitus, and 14 in M. mohoity) that may contribute to the species-specific adaptations were identified. An analysis of 1392 one-to-one orthologous genes revealed significantly higher ratios of nonsynonymous-to-synonymous substitutions in the M. mohoity lineage than in M. anguillicaudatus. The genes displaying signatures of positive selection and rapid evolution in Misgurnus were involved in four function categories, (1) energy metabolism; (2) signal transduction; (3) membrane; and (4) cell proliferation or apoptosis, implying that these candidate genes play critical roles in the thermal adaptation of the fish to their living environments. We also detected more than five positively selected sites in cldn15lb and isca1, which function as important factors in paracellular Na⁺ transport and Fe/S cluster assembly, respectively. Overall, our study provides valuable insights into the adaptive evolution of loaches from different temperature zones in China and is a foundation for future studies to clarify the genetic basis of temperature adaptation in fishes.

Transcriptome analysis and microsatellite discovery in the blunt snout bream (Megalobrama amblycephala) after challenge with Aeromonas hydrophila.

The blunt snout bream, Megalobrama amblycephala, is a herbivorous freshwater fish species native to China and a major aquaculture species in Chinese freshwater polyculture systems. In recent years, the bacterium Aeromonas hydrophila has been reported to be its pathogen causing great losses of farmed fish. To understand the immune response of the blunt snout bream to A. hydrophila infection, we used the Solexa/Illumina technology to analyze the transcriptomic profile after artificial bacterial infection. Two nonnormalized cDNA libraries were synthesized from tissues collected from control blunt snout bream or those injected with A. hydrophila. After assembly, 155,052 unigenes (average length 692.8 bp) were isolated. All unigenes were annotated using BLASTX relative to several public databases: the National Center for Biotechnology Information nonreduntant (Nr) database, SwissProt, Eukaryotic Orthologous Groups of proteins (KOG), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Ontology (GO). The sequence similarity (86%) of the assembled unigenes was to zebrafish based on the Nr database. A number of unigenes (n = 30,482) were assigned to three GO categories: biological processes (25,242 unigenes), molecular functions (26,096 unigenes), and cellular components (22,778 unigenes). 20,909 unigenes were classified into 25 KOG categories and 28,744 unigenes were assigned into 315 specific signaling pathways. In total, 238 significantly differentially expressed unigenes (mapped to 125 genes) were identified: 101 upregulated genes and 24 downregulated genes. Another 303 unigenes were mapped to unknown or novel genes. Among the known expressed genes identified, 53 were immune-related genes and were distributed in 71 signaling pathways. The expression patterns of selected up- and downregulated genes from the control and injected groups were determined with reverse transcription-quantitative PCR (RT-qPCR). Microsatellites (n = 10,877), including di-to pentanucleotide repeat motifs, were also identified in the blunt snout bream transcriptome profiles. This study extends our understanding of the immune defense mechanisms of the blunt snout bream against A. hydrophila and provides useful data for further studies of the immunogenetics of this species.

A Comprehensive Analysis of Codon Usage Patterns in Blunt Snout Bream (Megalobrama amblycephala) Based on RNA-Seq Data.

Blunt snout bream (Megalobrama amblycephala) is an important fish species for its delicacy and high economic value in China. Codon usage analysis could be helpful to understand its codon biology, mRNA translation and vertebrate evolution. Based on RNA-Seq data for M. amblycephala, high-frequency codons (CUG, AGA, GUG, CAG and GAG), as well as low-frequency ones (NUA and NCG codons) were identified. A total of 724 high-frequency codon pairs were observed. Meanwhile, 14 preferred and 199 avoided neighboring codon pairs were also identified, but bias was almost not shown with one or more intervening codons inserted between the same pairs. Codon usage bias in the regions close to start and stop codons indicated apparent heterogeneity, which even occurs in the flanking nucleotide sequence. Codon usage bias (RSCU and SCUO) was related to GC3 (GC content of 3rd nucleotide in codon) bias. Six GO (Gene ontology) categories and the number of methylation targets were influenced by GC3. Codon usage patterns comparison among 23 vertebrates showed species specificities by using GC contents, codon usage and codon context analysis. This work provided new insights into fish biology and new information for breeding projects.

Leptin Genes in Blunt Snout Bream: Cloning, Phylogeny and Expression Correlated to Gonads Development.

To investigate the leptin related genes expression patterns and their possible function during the gonadal development in fish, the cDNA and genomic sequences of leptin, leptin receptor (leptinR), and leptin receptor overlapping transcript like-1 (leprotl1) were cloned and their expression levels were quantified in the different gonadal development stages of Megalobrama amblycephala. The results showed that the full length cDNA sequences of leptin, leptinR and leprotl1 were 953, 3432 and 1676 bp, coding 168, 1082, and 131 amino acid polypeptides, and the genomic sequences were 1836, 28,528 and 5480 bp, which respectively had 3, 15 and 4 exons, respectively. The phylogenetic analysis revealed that three genes were relatively conserved in fish species. Quantitative real-time PCR results showed that the three genes were ubiquitously expressed in all examined tissues during the different gonadal development stages. The leptin and leptinR took part in the onset of puberty, especially in female M. amblycephala, by increasing the expression levels in brain during the stage I to III of ovary. The expression levels of leptin and leptinR had significant differences between male and female in hypothalamic-pituitary-gonadal (HPG) axis tissues (p < 0.05). The leptinR had the same variation tendency with leptin, but the opposite changes of expression levels were found in leprotl1, which may resist the expression of leptinR for inhibiting the function of leptin in target organ. These findings revealed details about the possible role of these genes in regulating gonadal maturation in fish species.

Identification of MicroRNA for Intermuscular Bone Development in Blunt Snout Bream (Megalobrama amblycephala).

Intermuscular bone (IB), which occurs only in the myosepta of the lower teleosts, is attracting more attention of researchers due to its particular development and lack of genetic information. MicroRNAs (miRNAs) are emerging as important regulators for biological processes. In the present study, miRNAs from IBs and connective tissue (CT; encircled IBs) from six-month-old Megalobrama amblycephala were characterized and compared. The results revealed the sequences and expression levels of 218 known miRNA genes (belonging to 97 families). Of these miRNAs, 44 known microRNA sequences exhibited significant expression differences between the two libraries, with 24 and 20 differentially-expressed miRNAs exhibiting higher expression in the CT and IBs libraries, respectively. The expressions of 11 miRNAs were selected to validate in nine tissues. Among the high-ranked predicted gene targets, differentiation, cell cycle, metabolism, signal transduction and transcriptional regulation were implicated. The pathway analysis of differentially-expressed miRNAs indicated that they were abundantly involved in regulating the development and differentiation of IBs and CT. This study characterized the miRNA for IBs of teleosts for the first time, which provides an opportunity for further understanding of miRNA function in the regulation of IB development.

A chromosome-level genome assembly of East Asia endemic minnow Zacco platypus.

Zacco platypus is an endemic colorful freshwater minnow that is intensively distributed in East Asia. In this study, two adult female individuals collected from Haihe River basin were used for karyotypic study and genome sequencing, respectively. The karyotype formula of Z. platypus is 2N = 48 = 18 M + 24SM/ST + 6 T. We used PacBio long-read sequencing and Hi-C technology to assemble a chromosome-level genome of Z. platypus. As a result, an 814.87 Mb genome was assembled with the PacBio long reads. Subsequently, 98.64% assembled sequences were anchored into 24 chromosomes based on the Hi-C data. The chromosome-level assembly contained 54 scaffolds with a N50 length of 32.32 Mb. Repeat elements accounted for 52.35% in genome, and 24,779 protein-coding genes were predicted, with 92.11% were functionally annotated with the public databases. BUSCO analysis yielded a completeness score of 96.5%. This high-quality genome assembly provides valuable resources for future functional genomic research, comparative genomics, and evolutionary studies of genus Zacco.

Integrated microbiome and metabolome analyses reveal the effects of low pH on intestinal health and homeostasis of crayfish (Procambarus clarkii).

Low pH (LpH) poses a significant challenge to the health, immune response, and growth of aquatic animals worldwide. Crayfish (Procambarus clarkii) is a globally farmed freshwater species with a remarkable adaptability to various environmental stressors. However, the effects of LpH stress on the microbiota and host metabolism in crayfish intestines remain poorly understood. In this study, integrated analyses of antioxidant enzyme activity, histopathological damage, 16S rRNA gene sequencing, and liquid chromatography-mass spectrometry (LC-MS) were performed to investigate the physiology, histopathology, microbiota, and metabolite changes in crayfish intestines exposed to LpH treatment. The results showed that LpH stress induced obvious changes in superoxide dismutase and catalase activities and histopathological alterations in crayfish intestines. Furthermore, 16S rRNA gene sequencing analysis revealed that exposure to LpH caused significant alterations in the diversity and composition of the crayfish intestinal microbiota at the phylum and genus levels. At the genus level, 14 genera including Bacilloplasma, Citrobacter, Shewanella, Vibrio, RsaHf231, Erysipelatoclostridium, Anaerorhabdus, Dysgonomonas, Flavobacterium, Tyzzerella, Brachymonas, Muribaculaceae, Propionivibrio, and Comamonas, exhibited significant differences in their relative abundances. The LC-MS analysis revealed 859 differentially expressed metabolites in crayfish intestines in response to LpH, including 363 and 496 upregulated and downregulated metabolites, respectively. These identified metabolites exhibited significant enrichment in 24 Kyoto Encyclopedia of Genes and Genomes pathways (p < 0.05), including seven and 17 upregulated and downregulated pathways, respectively. These pathways are mainly associated with energy and amino acid metabolism. Correlation analysis revealed a strong correlation between the metabolites and intestinal microbiota of crayfish during LpH treatment. These findings suggest that LpH may induce significant oxidative stress, intestinal tissue damage, disruption of intestinal microbiota homeostasis, and alterations in the metabolism in crayfish. These findings provide valuable insights into how the microbial and metabolic processes of crayfish intestines respond to LpH stress.

Identification and characterization of microRNAs involved in growth of blunt snout bream (Megalobrama amblycephala) by Solexa sequencing.

BACKGROUND: Blunt snout bream (Megalobrama amblycephala) is an economically important fish species in the Chinese freshwater polyculture system for its delicacy and high economic value. MicroRNAs (miRNAs) play important roles in regulation of almost all biological processes in eukaryotes. Although previous studies have identified thousands of miRNAs from many species, little information is known for miRNAs of M. amblycephala. To investigate functions of miRNAs associated with growth of M. amblycephala, we adopted the Solexa sequencing technology to sequence two small RNA libraries prepared from four growth related tissues (brain, pituitary, liver and muscle) of M. amblycephala using individuals with relatively high and low growth rates. RESULTS: In this study, we have identified 347 conserved miRNAs (belonging to 123 families) and 22 novel miRNAs in M. amblycephala. Moreover, we observed sequence variants and seed edits of the miRNAs. Of the 5,166 single nucleotide substitutions observed in two libraries, the most abundant were G-to-U (15.9%), followed by U-to-C (12.1%), G-to-A (11.2%), and A to G (11.2%). Subsequently, we compared the expression patterns of miRNAs in the two libraries (big-size group with high growth rate versus small-size group with low growth rate). Results indicated that 27 miRNAs displayed significant differential expressions between the two libraries (p < 0.05). Of these, 16 were significantly up-regulated and 11 were significantly down-regulated in the big-size group compared to the small-size group. Furthermore, stem-loop RT-PCR was applied to validate and profile the expression of the differentially expressed miRNAs in ten tissues, and the result revealed that the conserved miRNAs expressed at higher levels than the novel miRNAs, especially in brain, liver and muscle. Also, targets prediction of differentially expressed miRNAs and KEGG pathway analysis suggested that differentially expressed miRNAs are involved in growth and metabolism, signal transduction, cell cycle, neural development and functions. CONCLUSIONS: The present study provides the first large-scale characterization of miRNAs in M. amblycephala and miRNA profile related to different growth performances. The discovery of miRNA resource from this study is expected to contribute to a better understanding of the miRNAs roles playing in regulating the growth biological processes and the study of miRNA function and phenotype-associated miRNA identification in fish.

DNA methylation mediates gonadal development via regulating the expression levels of cyp19a1a in loach Misgurnus anguillicaudatus.

DNA methylation-mediated transcriptional regulation has been considered to significantly impact some steroidogenic enzyme genes expression. To uncover the roles of DNA methylation on the regulation of aromatase gene expression during gametogenesis in Misgurnus anguillicaudatus, the expression profiles and cellular localization of cyp19a1a and cyp19a1b were analyzed, and the landscape of DNA methylation dynamics was investigated. We found that cyp19a1a was predominantly expressed in granulosa cells of oocytes, while cyp19a1b expression was enriched in radial glial cells of the forebrain. In ovary, cyp19a1a was highly expressed until the vitellogenesis stage. The average methylation levels, especially for two CpG sites within the cAMP response element, were negatively correlated with cyp19a1a expression levels, indicating that methylation could regulate cyp19a1a transcriptional activity by modulating the binding efficiency of cAMP to its response elements. Compared with in ovary, cyp19a1a showed lower expression in testis but was hypermethylated. Cyp19a1b in female brain weakly expressed before the vitellogenesis stage, but significantly elevated at the maturation stage. In both sexes, it maintained high methylation levels in brain despite the obvious fluctuation of the cyp19a1b expression. This study revealed that DNA methylation plays a key role in establishing cyp19a1a spatiotemporal expression patterns and thus mediates gonadal development in teleosts.

Identification and Molecular Characterization of a Novel Carlavirus Infecting Chrysanthemum morifolium in China.

Chrysanthemum (Chrysanthemum morifolium) is an important ornamental and medicinal plant suffering from many viruses and viroids worldwide. In this study, a new carlavirus, tentatively named Chinese isolate of Carya illinoinensis carlavirus 1 (CiCV1-CN), was identified from chrysanthemum plants in Zhejiang Province, China. The genome sequence of CiCV1-CN was 8795 nucleotides (nt) in length, with a 68-nt 5'-untranslated region (UTR) and a 76-nt 3'-UTR, which contained six predicted open reading frames (ORFs) that encode six corresponding proteins of various sizes. Phylogenetic analyses based on full-length genome and coat protein sequences revealed that CiCV1-CN is in an evolutionary branch with chrysanthemum virus R (CVR) in the Carlavirus genus. Pairwise sequence identity analysis showed that, except for CiCV1, CiCV1-CN has the highest whole-genome sequence identity of 71.3% to CVR-X6. At the amino acid level, the highest identities of predicted proteins encoded by the ORF1, ORF2, ORF3, ORF4, ORF5, and ORF6 of CiCV1-CN were 77.1% in the CVR-X21 ORF1, 80.3% in the CVR-X13 ORF2, 74.8% in the CVR-X21 ORF3, 60.9% in the CVR-BJ ORF4, 90.2% in the CVR-X6 and CVR-TX ORF5s, and 79.4% in the CVR-X21 ORF6. Furthermore, we also found a transient expression of the cysteine-rich protein (CRP) encoded by the ORF6 of CiCV1-CN in Nicotiana benthamiana plants using a potato virus X-based vector, which can result in a downward leaf curl and hypersensitive cell death over the time course. These results demonstrated that CiCV1-CN is a pathogenic virus and C. morifolium is a natural host of CiCV1.

Alterations of the Gut Microbiota and Metabolomics Associated with the Different Growth Performances of Macrobrachium rosenbergii Families.

To investigate the key gut microbiota and metabolites associated with the growth performance of Macrobrachium rosenbergii families, 16S rRNA sequencing and LC-MS metabolomic methods were used. In this study, 90 M. rosenbergii families were bred to evaluate growth performance. After 92 days of culture, high (H), medium (M), and low (L) experimental groups representing three levels of growth performance, respectively, were collected according to the weight gain and specific growth rate of families. The composition of gut microbiota showed that the relative abundance of Firmicutes, Lachnospiraceae, Lactobacillus, and Blautia were much higher in Group H than those in M and L groups. Meanwhile, compared to the M and L groups, Group H had significantly higher levels of spermidine, adenosine, and creatinine, and lower levels of L-citrulline. Correlation analysis showed that the abundances of Lactobacillus and Blautia were positively correlated with the levels of alpha-ketoglutaric acid and L-arginine. The abundance of Blautia was also positively correlated with the levels of adenosine, taurine, and spermidine. Notably, lots of metabolites related to the metabolism and biosynthesis of arginine, taurine, hypotaurine, and fatty acid were upregulated in Group H. This study contributes to figuring out the landscape of the gut microbiota and metabolites associated with prawn growth performance and provides a basis for selective breeding.

Analysis of mRNA and MicroRNA Expression Profiles of Nervous Tissues and Reproductive Tissues in Male Procambarus clarkii After Silencing IAG.

The insulin-like androgenic gland hormone gene (IAG), primarily expressed in the androgenic gland (AG), plays a crucial role in controlling male sex differentiation and maintaining male secondary sexual characteristics in decapods. In this study, we investigated the mRNA and microRNA expression profiles of male Procambarus clarkii to understand the transcriptomic regulatory mechanism of IAG after the injection of an efficient siRNA (GsiRNA) designed based on IAG. The results revealed that several differentially expressed genes were enriched in reproduction-related pathways, such as the wnt signaling pathway, MAPK signaling pathway, and GnRH signaling pathway. In the testis (Te), the injection of GsiRNA led to the up-regulation of many ovary-related genes and down-regulation of testis-related genes. Moreover, the brain (Br) and abdominal nerve cord (AN) appeared to be involved in the regulation of IAG, with numerous differentially expressed genes found in Br and AN. Notably, the expression of five neuropeptide genes, Crustacean hyperglycemic hormone, pigment-dispersing hormone, red pigment concentrating hormone precursor, corazonin, and gonadotropin-releasing hormone II receptor isoform X1 in Br/AN, was significantly changed. Additionally, three ovary-related miRNAs (miR-263a, miR-263b, miR-133) highly expressed in Te/AG showed significant up-regulation after GsiRNA injection. Furthermore, the long-term interference of GsiRNA was found to inhibit the development of male external sexual characteristics during the juvenile stage and delay it during the adult stage. This research provides valuable insights into the molecular regulatory mechanism and function of IAG in P. clarkii.

Integrated Transcriptomic and Metabolomic Analyses Reveal Low-Temperature Tolerance Mechanism in Giant Freshwater Prawn Macrobrachium rosenbergii.

Water temperature, as an important environmental factor, affects the growth and metabolism of aquatic animals and even their survival. The giant freshwater prawn (GFP) Macrobrachium rosenbergii is a kind of warm-water species, and its survival temperature ranges from 18 °C to 34 °C. In this study, we performed transcriptomic and metabolomic analyses to clarify the potential molecular mechanism of responding to low-temperature stress in adult GFP. The treatments with low-temperature stress showed that the lowest lethal temperature of the GFP was 12.3 °C. KEGG enrichment analyses revealed that the differentially expressed genes and metabolites were both enriched in lipid and energy metabolism pathways. Some key genes, such as phosphoenolpyruvate carboxykinase and fatty acid synthase, as well as the content of the metabolites dodecanoic acid and alpha-linolenic acid, were altered under low-temperature stress. Importantly, the levels of unsaturated fatty acids were decreased in LS (low-temperature sensitive group) vs. Con (control group). In LT (low-temperature tolerant group) vs. Con, the genes related to fatty acid synthesis and degradation were upregulated to cope with low-temperature stress. It suggested that the genes and metabolites associated with lipid metabolism and energy metabolism play vital roles in responding to low-temperature stress. This study provided a molecular basis for the selection of a low-temperature tolerant strain.

The Dynamics of Gene Expression Unraveling the Immune Response of Macrobrachium rosenbergii Infected by Aeromonas veronii.

To further investigate the immune response of Macrobrachium rosenbergii against Aeromonas veronii, comparative transcriptomic analyses of the M. rosenbergii hepatopancreas were conducted on challenge and control groups at 6, 12, and 24 h post-infection (hpi), independently. A total of 51,707 high-quality unigenes were collected from the RNA-seq data, and 8060 differentially expressed genes (DEGs) were discovered through paired comparisons. Among the three comparison groups, a KEGG pathway enrichment analysis showed that 173 immune-related DEGs were considerably clustered into 28 immune-related pathways, including the lysosome, the phagosome, etc. Moreover, the expression levels of the four key immune-related genes (TOLL, PAK1, GSK3ß, and IKKα) were evaluated at various stages following post-infection in the hepatopancreas, hemolymph, and gills. Both PAK1 and GSK3ß genes were highly up-regulated in all three tissues at 6 hpi with A. veronii; TOLL was up-regulated in the hepatopancreas and hemolymph but down-regulated in the gill at 6 hpi, and IKKα was up-regulated in hemolymph and gill, but down-regulated in the hepatopancreas at 6 hpi. These findings lay the groundwork for understanding the immune mechanism of M. rosenbergii after contracting A. veronii.