Female common carp grow faster than male individuals, implying that rearing females could be more profitable in aquaculture. Non-coding RNAs (ncRNAs) serve as versatile regulators with multiple functions in diverse biological processes. However, the roles of ncRNAs in the sex differentiation of common carp are less studied. In this study, we investigated the expression profiles of ncRNAs, including miRNAs, lncRNAs, and circRNAs, in the gonads to comprehend the roles of ncRNAs in sex differentiation in common carp. A substantial number of differentially expressed (DE) ncRNAs in ovaries and testes were identified. Some miRNAs, notably miR-205, miR-214, and miR-460-5p, might modulate hormone synthesis and thus maintain sex. A novel miRNA, novel_158, was predicted to suppress the expression of foxl3. DE lncRNAs were associated with oocyte meiosis, GnRH signaling pathways, and steroid biosynthesis, while DE circRNA target genes were enriched in the WNT signaling pathway and MAPK signaling pathway. We also analyzed ncRNA-mRNA interactions to shed light on the crosstalk between competing endogenous RNAs (ceRNAs), which is the critical mechanism by which lncRNAs and circRNAs function. Some lncRNAs and circRNAs may be able to competitively bind novel_313, a new miRNA, and thus regulate hsd17β3. Our research will provide a valuable resource for understanding the genetic basis of gonadal differentiation and development in common carp.
Epigenetic modifications are critical in precisely regulating gene expression. The common carp (Cyprinus carpio) is an economically important fish species, and females exhibit faster growth rates than males. However, the studies related to epigenetic modifications in the common carp gonads are limited. In this study, we conducted the Assay for Transposase Accessible Chromatin sequencing (ATAC-seq) and Bisulfite sequencing (BS-seq) to explore the roles of epigenetic modifications in the common carp gonads. We identified 84,207 more accessible regions and 77,922 less accessible regions in ovaries compared to testes, and some sex-biased genes showed differential chromatin accessibility in their promoter regions, such as sox9a and zp3. Motif enrichment analysis showed that transcription factors (TFs) associated with embryonic development and cell proliferation were heavily enriched in ovaries, and the TFs Foxl2 and SF1 were only identified in ovaries. We also analyzed the possible regulations between chromatin accessibility and gene expression. By BS-seq, we identified 2087 promoter differentially methylated genes (promoter-DMGs) and 5264 gene body differentially methylated genes (genebody-DMGs) in CG contexts. These genebody-DMGs were significantly enriched in the Wnt signaling pathway, TGF-beta signaling pathway, and GnRH signaling pathway, indicating that methylation in gene body regions could play an essential role in sex maintenance, just like methylation in promoter regions. Combined with transcriptomes, we revealed that the expression of dmrtb1-like, spag6, and fels was negatively correlated with their methylation levels in promoter regions. Our study on the epigenetic modifications of gonads contributes to elucidating the molecular mechanism of sex differentiation and sex maintenance in the common carp.
The common carp (Cyprinus carpio), a major economic freshwater fish, is suffering from a variety of bacterial infectious diseases because of its high-density, factory and intensive farming patterns. Aeromonas veronii is the causative agent of high mortality in common carp, causing severe economic losses in aquaculture. However, the regulatory mechanisms involved in the response of common carp to this bacterial infection remain poorly understood. In this study, we compared mortality, blood serum LZM (Lysozyme) and IgM (Immunoglobulin M) levels and transcriptome patterns of head kidney tissues after infection with different concentrations of Aeromonas veronii. We observed that mortality increased progressively with an increasing pathogen concentration. The concentrations of blood serum LZM and IgM significantly increased after infection. A total of 13 and 925 differentially expressed genes (DEGs) were identified after infection with low (T4) and high (T9) concentrations of bacterial suspension, respectively. KEGG and GO analyses of the DEGs highlighted multiple immune-related signaling pathways. Weighted gene co-expression network analysis (WGCNA) revealed that 136 and 83 hub genes were related to blood serum LZM and IgM, respectively. Finally, the gene expression in the head kidney was validated via RT–qPCR to be consistent with the transcriptome. These results provide insights into the mechanisms of the immune response to infection with different concentrations of Aeromonas veronii and offer useful information for further studies on immune defense mechanisms in common carp.
The common carp (Cyprinus carpio) is a crucial freshwater species cultivated worldwide for food consumption. Female carp have better growth performance than males, which fascinates scholars to uncover the mechanism of gonadal differentiation and produce mono-sex populations. However, the mechanism of ovarian development at single-cell resolution is limited. Here, we conducted single-nucleus RNA sequencing in adolescent common carp ovaries. Our study obtained transcriptional profiles of 13,155 nuclei and revealed 13 distinct cell clusters in the ovaries, including three subtypes of germ cells and four subtypes of granulosa cells. We subsequently performed pseudotime trajectory analysis to delineate potential mechanisms underlying the development of germ cells and granulosa cells. We identified 1250 dynamic expression genes in germ cells and 1815 in granulosa cells (q-value < 0.01), including zp3, eif4a2 and aspm in germ cells and fshr and esr1 in granulosa cells. The functional annotation showed that dynamic expression genes in germ cells were involved in sperm–egg recognition and some terms related to meiosis, such as sister chromatid segregation and homologous recombination. Genes expressed dynamically in granulosa cells were related to the TGF-β signaling pathway, response to gonadotropin, and development of primary female sexual characteristics. In addition, the dynamic genes expressed in granulosa cells might relate to the complex communication between different cell types. In summary, our study provided a transcriptome profile of common carp ovaries at single-nucleus resolution, and we further revealed the potential cell type-specific mechanisms underlying oogenesis and the differentiation of granulosa cells, which will facilitate breeding all-female common carp populations.
DNA methylation and chromatin accessibility play important roles in gene expression, but their function in subgenome expression dominance remains largely unknown. We conducted comprehensive analyses of the transcriptome, DNA methylation, and chromatin accessibility in liver and muscle tissues of allotetraploid common carp, aiming to reveal the function of epigenetic modifications in subgenome expression dominance. A noteworthy overlap in differential expressed genes (DEGs) as well as their functions was observed across the two subgenomes. In the promoter and gene body, the DNA methylation level of the B subgenome was significantly different than that of the A subgenome. Nevertheless, differences in DNA methylation did not align with changes in homoeologous biased expression across liver and muscle tissues. Moreover, the B subgenome exhibited a higher prevalence of open chromatin regions and greater chromatin accessibility, in comparison to the A subgenome. The expression levels of genes located proximally to open chromatin regions were significantly higher than others. Genes with higher chromatin accessibility in the B subgenome exhibited significantly elevated expression levels compared to the A subgenome. Contrastingly, genes without accessibility exhibited similar expression levels in both subgenomes. This study contributes to understanding the regulation of subgenome expression dominance in allotetraploid common carp.
The uplift of the Qinghai–Tibet Plateau has shaped a unique extreme environment, fostering distinct endemic aquatic organisms. Schizopygopsis pylzovi, a vulnerable endemic fish in the upper Yellow River, is a key model for studying the biogeographic patterns of plateau fish. To assess its genetic characteristics and evolutionary dynamics, we comprehensively evaluated 11 geographic populations of S. pylzovi using two complementary mitochondrial markers, the conserved COI gene and the fast-evolving D-loop region. A total of 142 COI and 143 D-loop sequences were analyzed, and sequences alignment, haplotype network construction, AMOVA, and neutrality tests were performed. AMOVA revealed that genetic variation was mainly distributed within populations, indicating weak population differentiation. Neutrality tests and mismatch distribution analysis suggested historical and recent population expansion events. Our findings highlight the value of joint analysis using COI and D-loop markers in revealing the genetic structure of S. pylzovi, provide new insights into the impact of plateau uplift on fish evolution, and establish a scientific basis for the conservation of this vulnerable species.