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

Genome-wide identification and expressional profile of the Dmrt gene family in the swimming crab (Portunus trituberculatus)

The swimming crab, Portunus trituberculatus is one of crucial aquaculture crabs with significant differences in growth and economic performance between male and female swimming crabs. Consequently, the culture of female populations presents higher economic value. The doublesex and mab-3 related transcription factor (Dmrt) gene family are known to play crucial role in gonad differentiation and development. However, there is limited information about this gene family in Portunus trituberculatus. In this study, we identified seven members of the Dmrt gene family in P. trituberculatus based on the published transcriptome and genome data and designated as Ptdmrt-1, Ptdoublesex (Ptdsx), Ptidmrt-1, Ptdmrt-11E, Ptidmrt-2, Ptdmrt-99B, and Ptdmrt-3 based on the homology analysis results, respectively. These Ptdmrt genes distributed across 6 chromosomes and were predicted to encode 283 aa, 288 aa, 529 aa, 436 aa, 523 aa, 224 aa, and 435 aa protein precursors, respectively. The expression patterns of these dmrt genes were characterized by qRT-PCR and gonad transcriptome data. The results showed that five members (Ptdmrt-99B, Ptdsx, Ptdmrt-1, Ptdmrt-3, and Ptdmrt-11E) were differentially expressed between the testis and ovary. In addition, their expression patterns from zoea 2 to juvenile 1 were characterized by published transcriptome data and the results showed that they were lowly expressed and did not exhibit notable difference except for Ptdsx during early development. Overall, majority of Ptdmrt genes were involved in gonad differentiation in the swimming crab. Current findings provide a solid foundation for further exploration of the roles of these genes in gonad development and differentiation in P. trituberculatus.

Comparative transcriptomic characterization of the ovary in the spawning process of the mud crab Scylla paramamosain

Oviposition is induced upon mating in most insects. Spawning is a physiological process that is fundamental for the reproduction of Scylla paramamosain. However, the molecular mechanisms underlying the spawning process in this species are poorly understood. Herein, comprehensive ovary transcriptomic analysis was conducted at the germinal vesicle breakdown stage (GVBD), spawning stage, 0.5 h post-spawning stage, and 24 h post-spawning stage of S. paramamosain for gene discovery. A total of 67,230 unigenes were generated, and 27,975 (41.61%) unigenes were annotated. Meanwhile, the differentially expressed genes (DEGs) between the different groups were identified, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was subsequently conducted. These results suggested that octopamine (OA) and tyramine (TA) could induce oviposition, while dopamine (DA) and serotonin (5-hydroxytryptamine [5-HT]) inhibit oviposition. The 20-hydroxyecdysone (20E) and methyl farnesoate (MF) signal pathways might be positively associated with oviposition. Furthermore, numerous transcripts that encode neuropeptides and their G-protein-coupled receptors (GPCRs), such as CNMamide, RYamide, ecdysis-triggering hormone (ETH), GPA2/GPB5 receptor, and Moody receptor, appear to be differentially expressed during the spawning process. Eleven unigenes were selected for qRT-PCR and the pattern was found to be consistent with the transcriptome expression pattern. Our work is the first spawning-related investigation of S. paramamosain focusing on the ovary at the whole transcriptome level. These findings assist in improving our understanding of spawning regulation in S. paramamosain and provide information for oviposition studies in other crustaceans.

Effects of dietary melatonin on antioxidant and immune function of the Pacific white shrimp (Litopenaeus vannamei), as determined by transcriptomic analysis

Melatonin (MT) is regarded as an antioxidant and immunostimulant that can efficiently scavenge free radicals and activate antioxidant enzymes. The aim of this study was to investigate the effects of dietary MT on the growth performance and immune function of the Pacific white shrimp (Litopenaeus vannamei). Six groups of L. vannamei were supplemented with dietary MT at 0, 22.5, 41.2, 82.7, 165.1, and 329.2 mg/kg levels for 2 months. RNA-Seq analysis was performed to obtain transcriptome data of the control group and the group supplemented with dietary MT at 82.7 mg/kg BW. In total, 1220 DEGs (799 up-regulated and 421 down-regulated) were identified. Pathways and genes related to growth performance and immune function were verified by real-time quantitative polymerase chain reaction. The total hemocyte count, phagocytosis rate, and respiratory burst were significantly increased in the MT (82.7 mg/kg BW) group as compared to the control group. Analysis of antioxidant-related enzymes in the hepatopancreas showed that dietary MT (82.7 mg/kg BW) significantly increased activities of superoxide dismutase (SOD), glutathione peroxidase (GPx), and catalase, while dietary MT at 41.2 mg/kg BW significantly increased activities of glutathione S-transferase, lysozyme (LZM), and phenoloxidase (PO). At the transcriptional level, dietary MT up-regulated expression levels of genes associated with antioxidant immunity and growth, which included PO, SOD, LZM, GPx, chitin synthase, ecdysone receptor, calcium-calmodulin dependent protein kinase I, and retinoid X receptor. In conclusion, dietary MT may improve the growth performance and immune function of L. vannamei to some extent.

A Novel Rhipicephalus microplus Estrogen Related Receptor (RmERR), a Molecular and In Silico Characterization of a Potential Protein Binding Estrogen

The search for targets to control ticks and tick-borne diseases has been an ongoing problem, and so far, we still need efficient, non-chemical alternatives for this purpose. This search must consider new alternatives. For example genomics analysis is a widely applied tool in veterinary health studies to control pathogens. On the other hand, we propose that regulation of endocrine mechanisms represents a feasible alternative to biologically controlling tick infestations. Thus, we performed the molecular identification of an estrogen-related receptor gene of Rhipicephalus microplus called RmERR by RT-PCR in tick ovaries, embryonic cells, and hemolymph, which allowed us to analyze its expression and propose potential functions in endocrine mechanisms and developmental stages. In addition, we performed an in silico characterization to explore the molecular interactions of RmERR with different estrogens, estrogenic antagonists, and endocrine disruptor Bisphenol A (BPA), finding potential interactions predicted by docking analysis and supported by negative values of ΔG (which suggests the potential interaction of RmERR with the molecules evaluated). Additionally, phylogenetic reconstruction revealed that RmERR is grouped with other tick species but is phylogenetically distant from host vertebrates' ERRs. In summary, this study allowed for the identification of an ERR in cattle tick R. microplus for the first time and suggested its interaction with different estrogens, supporting the idea of a probable transregulation process in ticks. The elucidation of this interaction and its mechanisms unveiled its potential as a target to develop tick control strategies.

Transcriptomics and the origin of obligate parthenogenesis

Despite the presence of obligately parthenogenetic (OP) lineages derived from sexual ancestors in diverse phylogenetic groups, the genetic mechanisms giving rise to the OP lineages remain poorly understood. The freshwater microcrustacean Daphnia pulex typically reproduces via cyclical parthenogenesis. However, some populations of OP D. pulex have emerged due to ancestral hybridization and introgression events between two cyclically parthenogenetic (CP) species D. pulex and D. pulicaria. These OP hybrids produce both subitaneous and resting eggs parthenogenetically, deviating from CP isolates where resting eggs are produced via conventional meiosis and mating. This study examines the genome-wide expression and alternative splicing patterns of early subitaneous versus early resting egg production in OP D. pulex isolates to gain insight into the genes and mechanisms underlying this transition to obligate parthenogenesis. Our differential expression and functional enrichment analyses revealed a downregulation of meiosis and cell cycle genes during early resting egg production, as well as divergent expression patterns of metabolism, biosynthesis, and signaling pathways between the two reproductive modes. These results provide important gene candidates for future experimental verification, including the CDC20 gene that activates the anaphase-promoting complex in meiosis.

Histological, Physiological and Transcriptomic Analysis Reveal the Acute Alkalinity Stress of the Gill and Hepatopancreas of Litopenaeus vannamei

The pacific white shrimp (Litopenaeus vannamei) has gradually become a promising economic species in the development of saline-alkali water fishery. The study related to the stress reaction of pacific white shrimp under alkalinity stress is still limited, which is also a critical limiting factor for its saline-alkaline aquaculture. In this study, we aim to analyse the stress reaction of pacific white shrimp under acute alkalinity stress between control group (alkalinity:40 mg/L) and treatment group (alkalinity:350 mg/L) through histological observation, physiological determination and transcriptome. In the present study, during the process of acute alkalinity stress, the activities of Na+-K+-ATPase, carbonic anhydrase, sodium/hydrogen exchanger in gill related to homeostasis were significantly changed, the activities of superoxide dismutase and catalase related to antioxidant were decreased in both gill and hepatopancreas, and the activities of protease, lipase and amylase in hepatopancreas were decreased. At the same time, different degrees of histological damages were occured in the gill and hepatopancreas under acute alkalinity stress. There were 194 and 236 different expressed genes identified in gill and hepatopancreas respectively. Functional enrichment assessment indicated that the alkalinity stress-related genes in both gill and hepatopancreas were primarily involved in fatty acid metabolism, glycolysis/gluconeogenesis, glycerophospholipid metabolism. The results indicated that the functions of homeostasis regulation, antioxidation and digestion of pacific white shrimp were decreased under acute alkalinity stress, at the same time, the energy metabolism in gill and hepatopancreas were modified to cope with alkalinity stress. This work provides important clues for understanding the response mechanism of pacific white shrimp under acute alkalinity stress.

Evaluation of the Role of Soybean Lecithin, Egg Yolk Lecithin, and Krill Oil in Promoting Ovarian Development in the Female Redclaw Crayfish Cherax quadricarinatus

The optimal supplementation of lipid nutrients in the diet showed crucial physiological functions in gonadal development and maturation in adult female aquatic animals. Four isonitrogenous and isolipidic diets with no extra lecithin supplementation (control), 2% soybean lecithin (SL), egg yolk lecithin (EL), or krill oil (KO) supplementation were formulated for Cherax quadricarinatus (72.32 ± 3.58 g). Ovary development and physiological characteristics of crayfish were evaluated after a 10-week feeding trial. The results indicated that SL, EL, or KO supplementation all significantly increased the gonadosomatic index, especially in the KO group. Crayfish fed the diet with SL showed the highest hepatosomatic index compared with those fed the other experimental diets. KO was more efficient than SL and EL in promoting triacylglycerol and cholesterol deposition in the ovary and hepatopancreas but also showed the lowest concentration of low-density lipoprotein cholesterol in the serum. KO significantly increased yolk granule deposition and accelerated oocyte maturation than other experimental groups. Furthermore, dietary phospholipids significantly enhanced the gonad-stimulating hormone concentration in the ovary and reduced the secretion of gonad-inhibiting hormones in the eyestalk. KO supplementation also significantly improved organic antioxidant capacity. From the ovarian lipidomics results, phosphatidylcholine and phosphatidylethanolamine are two main glycerophospholipids that respond to different dietary phospholipids. Polyunsaturated fatty acids (especially C18:2n-6, C18:3n-3, C20:4n-6, C20:5n-3, and C22:6n-3) were pivotal participants during ovarian development of crayfish regardless of lipid type. Combined with the ovarian transcriptome, the best positive function of KO was due to activated steroid hormone biosynthesis, sphingolipid signaling, retinol metabolism, lipolysis, starch and sucrose metabolism, vitamin digestion and absorption, and pancreatic secretion. As a consequence, dietary supplementation with SL, EL, or KO all improved the ovarian development quality of C. quadricarinatus, especially KO, which was the optimum choice for promoting ovary development in adult female C. quadricarinatus.

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.

Reproductive toxicity induced by benzo[a]pyrene exposure: first exploration highlighting the multi-stage molecular mechanism in female scallop Chlamys farreri

Reproductive toxicity induced by benzo[a]pyrene (B[a]P) exposure has received great ecotoxicological concerns. However, huge gaps on the molecular mechanism still exist in bivalves. In this study, reproduction-related indicators were investigated in female scallops Chlamys farreri during life cycle of proliferative, growth, mature, and spawn stages, under gradient concentrations of B[a]P at 0, 0.04, 0.4, and 4 μg/L. Meanwhile, a multi-stage ovarian transcriptome analysis under 4 μg/L B[a]P exposure was also conducted to elucidate the potential molecular mechanisms. The results indicated that life-cycle exposure to 0.4 and 4 μg/L B[a]P significantly decreased GSI and sex steroid levels. Even 0.04 μg/L B[a]P could play the adverse role in DNA integrity at the mature and spawn stages. Ovarian histological sections showed that B[a]P inhibited the maturation and release of oocytes. Through the functional enrichment analysis of differentially expressed genes (DEGs) from transcriptome data, 18 genes involved in endocrine disruption effects, DNA damage and repair, and oogenesis were selected and further determined by qRT-PCR. The downregulation of genes involved in steroidogenic and estrogen signaling pathways indicated that B[a]P could cause endocrine disruption through both receptor-dependent and receptor-independent pathways. The variations of gene expressions involved in DNA single-strand break and repair implied the presence of toxic mechanisms similar with vertebrates. Additionally, the changes of gene expressions of cell cycle, apoptosis, and cell adhesion suggested that exposure to B[a]P possibly caused the reproductive toxicity effects by affecting oogenesis. Taken together, this study was a pioneer in combining genome-wide transcriptomic analysis with its corresponding reproductive indicators (GSI, sex steroid levels, DNA single-strand break, and histological sections) to explore the bivalves' toxic mechanisms under B[a]P exposure. Meanwhile, some genes involved in estrogen signaling pathway and DNA damage were firstly analyzed in bivalves, and the expression data might be useful in establishing new hypotheses and discovering new biomarkers for marine biomonitoring.

Comparative transcriptomic analysis of ovaries from high and low egg-laying Lingyun black-bone chickens

Egg‐laying rate is mainly determined by ovarian function and regulated by the hypothalamic‐pituitary‐gonadal axis; however, the mechanism by which the ovary regulates the egg‐laying rate is still poorly understood. The purpose of this study was to compare the differences in the transcriptomes of the ovary of Lingyun black‐bone chickens with relatively high and low egg‐laying rates and screen candidate genes related to the egg‐laying rate. RNA‐sequencing (RNA‐Seq) was conducted to explore the chicken transcriptome from the ovarian tissue of six Lingyun black‐bone chickens with high (group G, n = 3) and low (group D, n = 3) egg‐laying rates. The results showed that 235 differentially expressed genes (DEGs) were identified between the chickens with high and low egg‐laying rates; among them, 209 DEGs were up‐regulated and 26 DEGs were down‐regulated. Gene Ontology analysis showed that the up‐regulated 209 DEGs were enriched in 50 GO terms and the down‐regulated 26 DEGs were enriched in 40 GO terms. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that up‐regulated DEGs were significantly enriched in 25 pathways and down‐regulated DEGs were significantly enriched in three pathways. Among the pathways, we found the longevity regulating pathway‐multiple species pathway, Estrogen signalling pathway and PPAR signalling pathway may have an essential function in regulating the egg‐laying rate. The results highlighted DEGs in the ovarian tissues of relatively high and low laying Lingyun black‐bone chicken and identified essential candidate genes related to the egg‐laying rate, thereby providing a theoretical basis for improving the egg‐laying rate of Lingyun black‐bone chicken.