ceRNA crosstalk mediated by ncRNAs is a novel regulatory mechanism in fish sex determination and differentiation

The potential regulatory role of the non-coding RNAs in regulating the exogenous estrogen-induced feminization in Takifugu rubripes gonad

Estrogen plays a pivotal role in the early stage of sex differentiation in teleost. However, the underlying mechanisms of estrogen-induced feminization process are still needed for further clarification. Here, the comparative analysis of whole-transcriptome RNA sequencing was conducted between 17beta-Estradiol induced feminized XY (E-XY) gonads and control gonads (C) in Takifugu rubripes. A total of 57 miRNAs, 65 lncRNAs, and 4 circRNAs were found to be expressed at lower levels in control-XY (C-XY) than that in control-XX (C-XX), and were up-regulated in XY during E2-induced feminization process. The expression levels of 24 miRNAs, and 55 lncRNAs were higher in C-XY than that in C-XX, and were down-regulated in E2-treated XY. Furthermore, a correlation analysis was performed between miRNA-seq and mRNA-seq data. In C-XX/C-XY, 114 differential expression (DE) miRNAs were predicted to target to 904 differential expression genes (DEGs), while in C-XY/E-XY, 226 DEmiRNAs were predicted to target to 2,048 DEGs. In C-XX/C-XY, and C-XY/E-XY, KEGG pathway enrichment analysis showed that those targeted genes were mainly enriched in MAPK signaling, calcium signaling, steroid hormone biosynthesis and ovarian steroidogenesis pathway. Additionally, the competitive endogenous RNA (ceRNA) regulatory network was constructed by 24 miRNAs, 21 lncRNAs, 4 circRNAs and 5 key sex-related genes. These findings suggested that the expression of critical genes in sex differentiation were altered in E2-treated XY T. rubripes may via the lncRNA-miRNA-mRNA regulation network to facilitate the differentiation and maintenance of ovaries. Our results provide a new insight into the comprehensive understanding of the effects of estrogen signaling pathways on sex differentiation in teleost gonads.

Molecular functional characterization of the setdb1 and its potential target gene sox5 illuminate the histone modification-mediated orchestration of gonadal development in Chinese tongue sole (Cynoglossus semilaevis)

SET (SuVar3-9, Enhancer of Zeste, Trithorax) domain bifurcated histone lysine methyltransferase 1, setdb1, is the predominant histone lysine methyltransferase catalyzing H3K9me3. Prior studies have illustrated that setdb1 and H3K9me3 critically regulate sex differentiation and gametogenesis. However, the molecular details by which setdb1 is involved in these processes in fish have been poorly reported. Here, we cloned and characterized the setdb1 ORF (open reading frame) sequence from Chinese tongue sole (Cynoglossus semilaevis). The setdb1 ORF sequence was 3,669 bp, encoding a 1,222-amino-acid protein. Phylogenetic analysis showed that setdb1 was structurally conserved. qRT-PCR revealed that setdb1 had a high expression level in the testes at 12 mpf (months post fertilization). Single-cell RNA-seq data at 24 mpf indicated that setdb1 was generally expressed in spermatogenic cells at each stage except for sperm and was centrally expressed in oogonia. H3K9me3 modification was observed in gonads with the immunofluorescence technique. Furthermore, the overexpression experiment suggested that sox5 was a candidate target of setdb1. sox5 was abundantly expressed in male and pseudomale gonads at 24 mpf. Single-cell RNA-seq data showed that sox5 was mainly expressed in spermatogonia and its expression gradually declined with differentiation. Taken together, our findings imply that setdb1 regulates sox5 transcription in gonads, which provides molecular clues into histone modification-mediated orchestration of sex differentiation and gametogenesis.

Integrated analysis of lncRNA, miRNA and mRNA expression profiles reveals regulatory pathways associated with pig testis function

Long noncoding RNA (lncRNA) and microRNA (miRNA) are known to play pivotal roles in mammalian testicular function and spermatogenesis. However, their impact on porcine male reproduction has yet to be well unraveled. Here, we sequenced and identified lncRNA and miRNA expressed in the testes of Chinese indigenous Banna mini-pig inbred line (BMI) and introduced Western Duroc (DU) and Large White (LW) pigs. By pairwise comparison (BMI vs DU, BMI vs LW, and DU vs LW), we found the gene expression differences in the testes between Chinese local pigs and introduced Western commercial breeds were more striking than those between introduced commercial breeds. Furthermore, we found 1622 co-differentially expressed genes (co-DEGs), 122 co-differentially expressed lncRNAs (co-DELs), 39 co-differentially expressed miRNAs (co-DEMs) in BMI vs introduced commercial breeds (DU and LW). Functional analysis revealed that these co-DEGs and co-DELs/co-DEMs target genes were enriched in male sexual function pathways, including MAPK, AMPK, TGF-β/Smad, Hippo, NF-kappa B, and PI3K/Akt signaling pathways. Additionally, we established 10,536 lncRNA-mRNA, 11,248 miRNA-mRNA pairs, and 62 ceRNA (lncRNA-miRNA-mRNA) networks. The ssc-miR-1343 had the most interactive factors in the ceRNA network, including 20 mRNAs and 3 lncRNAs, consisting of 56 ceRNA pairs. These factors played extremely important roles in the regulation of testis function as key nodes in the interactive regulatory network. Our results provide insight into the functional roles of lncRNAs and miRNAs in porcine testis and offer a valuable resource for understanding the differences between Chinese indigenous and introduced Western pigs.

Identification of lncRNA-based regulatory mechanisms of Takifugu rubripes growth traits in fast and slow-growing family lines

Family selection is an important method in fish aquaculture because growth is the most important economic trait. Fast-and slow-growing families of tiger puffer fish (Takifugu rubripes) have been established through family selection. The development of teleost fish is primarily controlled by the growth hormone (GH)-insulin-like growth factor 1 (IGF-1) axis that includes the hypothalamus-pituitary-liver. In this study, the molecular mechanisms underlying T. rubripes growth were analyzed by comparing transcriptomes from fast- and slow-growing families. The expressions of 214 lncRNAs were upregulated, and those of 226 were downregulated in the brain tissues of the fast-growing T. rubripes family compared to those of the slow-growing family. Differentially expressed lncRNAs centrally regulate mitogen-activated protein kinase (MAPK) and forkhead box O (FoxO) signaling pathways. Based on the results of lncRNA-gene network construction, we found that lncRNA3133.13, lncRNA23169.1, lncRNA23145.1, and lncRNA23141.3 regulated all four genes (igf1, mdm2, flt3, and cwf19l1). In addition, lncRNA7184.10 may be a negative regulator of rasgrp2 and a positive regulator of gadd45ga, foxo3b, and dusp5. These target genes are associated with the growth and development of organisms through the PI3K/AKT and MAPK/ERK pathways. Overall, transcriptomic analyses of fast- and slow-growing families of T. rubripes provided insights into the molecular mechanisms of teleost fish growth rates. Further, these analyses provide evidence for key genes related to growth regulation and the lncRNA expression regulatory network that will provide a framework for improving puffer fish germplasm resources.

The ncRNA-TGF-β axis: Unveiling new frontiers in colorectal cancer research

Colorectal cancer (CRC) poses a substantial global challenge, necessitating a deeper understanding of the molecular underpinnings governing its onset and progression. The transforming growth factor beta (TGF-β) network has been a well-recognized cornerstone in advancing CRC. Nevertheless, a recent study has highlighted the growing importance of non-coding RNAs (ncRNAs) in this context. This comprehensive review aims to present an extensive examination of the interaction between ncRNAs and TGF-signaling. Noncoding RNAs (ncRNAs), encompassing circular RNAs (circRNAs), long-ncRNAs (lncRNAs), and microRNAs (miRNAs), have surfaced as pivotal modulators governing various aspects of TGF-β signaling. MiRNAs have been discovered to target elements within the TGF-β signaling, either enhancing or inhibiting signaling, depending on the context. LncRNAs have been associated with CRC progression, functioning as miRNA sponges or directly influencing TGF-β pathway elements. Even circRNAs, a relatively recent addition to the ncRNA family, have impacted CRC, affecting TGF-β signaling through diverse mechanisms. This review encompasses recent progress in comprehending specific ncRNAs involved in TGF-β signaling, their functional roles, and their clinical relevance in CRC. We investigate the possibility of ncRNAs as targets for detection, prognosis, and therapy. Additionally, we explore the interaction of TGF-β and other pathways in CRC and the role of ncRNAs within this intricate network. As we unveil the intricate regulatory function of ncRNAs in the TGF-β signaling in CRC, we gain valuable insights into the disease's pathogenesis. Incorporating these discoveries into clinical settings holds promise for more precise diagnosis, prognosis, and targeted therapeutic approaches, ultimately enhancing the care of CRC patients. This comprehensive review underscores the ever-evolving landscape of ncRNA research in CRC and the potential for novel interventions in the battle against this formidable disease.

Histomorphic analysis and expression of mRNA and miRNA in embryonic gonadal differentiation in Chinese soft-shelled turtle (Pelodiscus sinensis)

The Chinese soft-shelled turtle (Pelodiscus sinensis) is extensively cultured in Asia for its nutritional and medical value. Gonadal differentiation is fantastic in turtles, whereas morphologic, mRNA, and miRNA expressions were insufficient in the turtle. In this study, ovaries and testes histomorphology analysis of 14-23 stage embryos were performed, and mRNA and miRNA expression profiles were analyzed. Histomorphology analysis revealed that gonads were undifferentiated at embryonic stage 14. Ovarian morphological differentiation became evident from stage 15, which was characterized by the development of the cortical region and degeneration of the medullary region. Concurrently, testicular morphological differentiation was apparent from stage 15, marked by the development of the medullary region and degeneration of the cortical region. qRT-PCR results showed that Cyp19a1 and Foxl2 exhibited female-specific expression at stage 15 and the expression increased throughout most of the embryonic development. Dmrt1, Amh, and Sox9 displayed male-specific expression at stage 15 and tended to increase substantially at later developmental stages. The expression of miR-8356 and miR-3299 in ZZ gonads were significantly higher than that in ZW gonads at stage 15, 17 and 19, and they had the highest expression at stage 15. While the expression of miR-8085 and miR-7982 had the highest expression at stage 19. Furthermore, chromatin remodeler genes showed differential expression in female and male P. sinensis gonads. These results of master sex-differentiation genes and morphological characteristics would provide a reference for the research of sex differentiation and sex reversal in turtles. Additionally, the expression of chromatin remodeler genes indicated they might be involved in gonadal differentiation of P. sinensis.

Comprehensive analysis of differentially expressed mRNA, lncRNA and miRNA, and their ceRNA networks in the regulation of shell color in the Manila clam (Ruditapes philippinarum)

The regulatory mechanism of ceRNA network plays an important role in molecular function and biological processes, however, the molecular mechanism in the shell color of Ruditapes philippinarum has not yet been reported. In this study, we performed transcriptome sequencing on the mantle of R. philippinarum with different shell colors, and screened for mRNA, miRNA, and lncRNA. A total of 61 mRNAs, 3725 lncRNAs and 90 miRNAs were obtained from all the shell color comparison groups (all mRNAs, lncRNAs and miRNAs P < 0.05), and 7 mRNAs, 8 lncRNAs, and 4 miRNAs of the porphyrin pathway and melanin pathway were screened for competitive endogenous RNA (ceRNA) network construction. The results indicate that the ceRNA network composed of mRNA and lncRNA, centered around efu-miR-101, mle-bantam-3p, egr-miR-9-5p, and sma-miR-75p, may play a crucial regulatory role in shell color formation. This study reveals for the first time the mechanism of ceRNA regulatory networks in the shell color of R. philippinarum and providing important reference data for molecular breeding of shell color in R. philippinarum.

Differential Expression of miRNAs, lncRNAs, and circRNAs between Ovaries and Testes in Common Carp (Cyprinus carpio)

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.

Exosome delivery to the testes for dmrt1 suppression: A powerful tool for sex-determining gene studies

Exosomes are endosome-derived extracellular vesicles about 100 nm in diameter. They are emerging as promising delivery platforms due to their advantages in biocompatibility and engineerability. However, research into and applications for engineered exosomes are still limited to a few areas of medicine in mammals. Here, we expanded the scope of their applications to sex-determining gene studies in early vertebrates. An integrated strategy for constructing the exosome-based delivery system was developed for efficient regulation of dmrt1, which is one of the most widely used sex-determining genes in metazoans. By combining classical methods in molecular biology and the latest technology in bioinformatics, isomiR-124a was identified as a dmrt1 inhibitor and was loaded into exosomes and a testis-targeting peptide was used to modify exosomal surface for efficient delivery. Results showed that isomiR-124a was efficiently delivered to the testes by engineered exosomes and revealed that dmrt1 played important roles in maintaining the regular structure and function of testis in juvenile fish. This is the first de novo development of an exosome-based delivery system applied in the study of sex-determining gene, which indicates an attractive prospect for the future applications of engineered exosomes in exploring more extensive biological conundrums.

Regulatory mechanism of LncRNAs in gonadal differentiation of hermaphroditic fish, Monopterus albus

BACKGROUND

Monopterus albus is a hermaphroditic fish with sex reversal from ovaries to testes via the ovotestes in the process of gonadal development, but the molecular mechanism of the sex reversal was unknown.

METHODS

We produced transcriptomes containing mRNAs and lncRNAs in the crucial stages of the gonad, including the ovary, ovotestis and testis. The expression of the crucial lncRNAs and their target genes was detected using qRT‒PCR and in situ hybridization. The methylation level and activity of the lncRNA promoter were analysed by applying bisulfite sequencing PCR and dual-luciferase reporter assays, respectively.

RESULTS

This effort revealed that gonadal development was a dynamic expression change. Regulatory networks of lncRNAs and their target genes were constructed through integrated analysis of lncRNA and mRNA data. The expression and DNA methylation of the lncRNAs MSTRG.38036 and MSTRG.12998 and their target genes Psmβ8 and Ptk2β were detected in developing gonads and sex reversal gonads. The results showed that lncRNAs and their target genes exhibited consistent expression profiles and that the DNA methylation levels were negatively regulated lncRNA expression. Furthermore, we found that Ptk2β probably regulates cyp19a1 expression via the Ptk2β/EGFR/STAT3 pathway to reprogram sex differentiation.

CONCLUSIONS

This study provides novel insight from lncRNA to explore the potential molecular mechanism by which DNA methylation regulates lncRNA expression to facilitate target gene transcription to reprogram sex differentiation in M. albus, which will also enrich the sex differentiation mechanism of teleosts.

Comprehensive analysis on the regulation of differentially expressed of mRNA and ncRNA in different ovarian stages of ark shell Scapharca broughtonii

BACKGROUND

Ovarian development is an important prerequisite and basis for animal reproduction. In many vertebrates, it is regulated by multiple genes and influenced by sex steroid hormones and environmental factors. However, relative information is limited in shellfish. To explore the biological functions and molecular mechanisms of mRNA and non-coding RNA that regulate ovarian development in Scapharca broughtonii, we performed whole transcriptome sequencing analysis on ovaries at three developmental stages. Furthermore, the biological processes involved in the differential expression of mRNA and ncRNA were analyzed.

RESULTS

A total of 11,342 mRNAs, 6897 lncRNAs, 135 circRNAs, and 275 miRNAs were differentially expressed. By mapping the differentially expressed RNAs from the three developmental stages of Venn diagram, multiple groups of shared mRNAs and lncRNAs were found to be associated with ovarian development, with some mRNA and ncRNA functions associated with steroid hormone. In addition, we constructed and visualized the lncRNA/circRNA-miRNA-mRNA network based on ceRNA targeting relationships.

CONCLUSIONS

These findings may facilitate our further understanding the mRNA and ncRNAs roles in the regulation of shellfish reproduction.

High-density lipoprotein regulates angiogenesis by long non-coding RNA HDRACA

Normal high-density lipoprotein (nHDL) can induce angiogenesis in healthy individuals. However, HDL from patients with coronary artery disease undergoes various modifications, becomes dysfunctional (dHDL), and loses its ability to promote angiogenesis. Here, we identified a long non-coding RNA, HDRACA, that is involved in the regulation of angiogenesis by HDL. In this study, we showed that nHDL downregulates the expression of HDRACA in endothelial cells by activating WW domain-containing E3 ubiquitin protein ligase 2, which catalyzes the ubiquitination and subsequent degradation of its transcription factor, Kruppel-like factor 5, via sphingosine 1-phosphate (S1P) receptor 1. In contrast, dHDL with lower levels of S1P than nHDL were much less effective in decreasing the expression of HDRACA. HDRACA was able to bind to Ras-interacting protein 1 (RAIN) to hinder the interaction between RAIN and vigilin, which led to an increase in the binding between the vigilin protein and proliferating cell nuclear antigen (PCNA) mRNA, resulting in a decrease in the expression of PCNA and inhibition of angiogenesis. The expression of human HDRACA in a hindlimb ischemia mouse model inhibited the recovery of angiogenesis. Taken together, these findings suggest that HDRACA is involved in the HDL regulation of angiogenesis, which nHDL inhibits the expression of HDRACA to induce angiogenesis, and that dHDL is much less effective in inhibiting HDRACA expression, which provides an explanation for the decreased ability of dHDL to stimulate angiogenesis.

Phosphoproteomics Reveal New Candidates in Abnormal Spermatogenesis of Pseudomales in Cynoglossus semilaevis

Phosphorylation is a post-translational modification that contributes to versatile protein functions in spermatogenesis, and the variations they generate usually results in abnormal spermatogenesis or sperm dysfunction. The sex-reversal phenomenon exists in Chinese tongue sole under certain conditions such that individuals with a ZW genotype can acquire a male phenotype and are thus called pseudomales. Pseudomale tongue sole can reach sexual maturity but produce only Z-type sperm, and the Z sperm carries paternal epigenetic information. Whether phosphorylation plays a role in the sperm abnormality of pseudomales is unknown. In this study, a phosphoproteomic analysis was performed to compare protein phosphorylation profiles between pseudomale and male testes. Altogether, we identified 14,253 phosphopeptides matching with 4843 proteins, with 1329 differentially phosphorylated peptides corresponding to 1045 differentially phosphorylated proteins (DPPs). Phosphorylation at 781 sites was upregulated and at 548 sites was downregulated. Four motifs were identified among differentially phosphorylated peptides, which were "SP", "SD", "RxxS", and "TP". Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses suggested that the cell cycle and DNA/RNA processing were significantly enriched with the genes encoding DPPs. To analyze DPP function in depth, a protein-protein interaction network was constructed, and Ran-binding protein 2 was found to play a central role in spermatogenesis by regulating several processes such as the cell cycle, eukaryotic translation, ubiquitination, and minichromosome maintenance. In kinase-associated network analyses, two "mitogen-activated protein kinase (Mapk)-centered" clusters were identified that may account for abnormal spermatogenesis in pseudomales. One cluster was centered on Mapk6, which predominantly regulated the cell cycle by interacting with several cyclin-dependent kinases, and the other was centered on the "testis-expressed kinase 1-like (Tesk1l)/Pim1l-Mapk4l- testis-expressed 14 (Tex14)" kinase cascade, which might contribute to spermatogenesis by regulating β-catenin. Taken together, these data suggested the new candidates involved in pseudomale sperm abnormalities and provided clues to discover the phosphorylated regulatory mechanism underlying tongue sole spermatogenesis.

BmNPV circular RNA-encoded peptide VSP39 promotes viral replication

CircRNAs are covalently closed single-stranded circular RNA molecules, which are not easily degraded by endonucleases and play vital roles in many biological processes. Currently, most studies on circRNAs focus on endogenous circRNAs in cells, and there are few studies on virus-encoded circRNAs. In this study, a viral circRNA (circRNA-000010) derived from the region (-/bp: 114514-115,319) of the complementary strand of Bombyx mori Nucleopolyhedrovirus (BmNPV) genome was identified with the circRNA-sequencing. The authenticity of viral circRNA-000010 was further confirmed by reverse transcription PCR, reverse transcription-rolling circle amplification (TCA), in situ hybridization, immunofluorescent staining, and Northern blotting. The results of overexpression and knockdown experiments showed that circRNA-000010 promoted viral replication. Furthermore, a viral small peptide VSP39 with 39 amino acid residues translated by circRNA-000010 but not its linear molecule was confirmed. Finally, VSP39 was found to promote viral replication. Our findings indicated that a viral circRNA encoded by BmNPV promoted viral replication. These findings will provide new clues for further understanding coding information of the BmNPV genome and open a new insight for investigating host-virus interactions.

CircRNA Identification and CircRNA-miRNA-mRNA Network in Cynoglossus semilaevis Sexual Size Dimorphism

Sexual size dimorphism (SSD), which is the sexual differences in body size, has been widely reported in various species including fishes. For Chinese tongue sole (Cynoglossus semilaevis), a flatfish exhibiting typically female-biased SSD, little is known for its epigenetic regulation mechanism, especially the role of circRNAs. Here, we identified the differently expressed abundances of circRNAs in females, males, and pseudo-males to explore the potential functions of circRNAs in Chinese tongue sole SSD. In total, 14,745 novel circRNAs were screened, among which 1461 DE circRNAs were identified from the brain, gonad, liver, and muscle in female, male, and pseudo-male individuals. The ceRNA network was subsequently constructed, including 10 circRNAs, 26 mRNAs, and 11 miRNAs. These DE mRNAs were mainly related to the mRNA surveillance pathway, metabolic pathways, and cellular senescence. Importantly, the ceRNA network has revealed that several circRNAs such as novel_circ_004374 and novel_circ_014597 may regulate homeodomain interacting protein kinase 2 (hipk2) expression by sponging miR-130-x. It is also worth exploring whether or how novel_circ_008696 regulates SET Domain Containing 2, histone lysine methyltransferase (setd2), which in turn affects the epigenetic patterns of different sexual individuals. The present study not only enriches the knowledge on the potential roles of circRNA in the physiological process, but also provides new clues for the explanation of fish SSD. In future studies, the precise function and involvement of circRNAs in female-biased SSD will require more efforts.