Role of mRNAs and long non-coding RNAs in regulating the litter size trait in Chuanzhong black goats

Transcriptome analysis reveals pituitary lncRNA, circRNA and mRNA affecting fertility in high- and low-yielding goats

The pituitary gland serves as the central endocrine regulator of growth, reproduction, and metabolism and plays a crucial role in the reproductive process of female animals. Transcriptome analysis was conducted using pituitary gland samples from Leizhou goats with varying levels of fecundity to investigate the effects of long noncoding RNA (lncRNA), circular RNA (circRNA), and mRNA regulation on pituitary hormone secretion and its association with goat fecundity. The analysis aimed to identify lncRNAs, circRNAs, and mRNAs that influence the fertility of Leizhou goats. GO and KEGG enrichment analyses were performed on differentially expressed lncRNAs, circRNAs, and mRNAs and revealed considerable enrichment in pathways, such as regulation of hormone secretion, germ cell development, and gonadotropin-releasing hormone secretion. The pituitary lncRNAs (ENSCHIT00000010293, ENSCHIT00000010304, ENSCHIT00000010306, ENSCHIT00000010290, ENSCHIT00000010298, ENSCHIT00000006769, ENSCHIT00000006767, ENSCHIT00000006921, and ENSCHIT00000001330) and circRNAs (chicirc_029285, chicirc_026618, chicirc_129655, chicirc_018248, chicirc_122554, chicirc_087101, and chicirc_078945) identified as differentially expressed regulated hormone secretion in the pituitary through their respective host genes. Additionally, differential mRNAs (GABBR2, SYCP1, HNF4A, CBLN1, and CDKN1A) influenced goat fecundity by affecting hormone secretion in the pituitary gland. These findings contribute to the understanding of the molecular mechanisms underlying pituitary regulation of fecundity in Leizhou goats.

Transcriptome analysis reveals mRNAs and long non-coding RNAs associated with fecundity in the hypothalamus of high-and low-fecundity goat

As an important organ that coordinates the neuroendocrine system, the hypothalamus synthesizes and secretes reproductive hormones that act on the goat organism, thereby precisely regulating follicular development and reproductive processes in goats. However, it is still elusive to explore the mechanism of hypothalamic effects on goat fertility alone. Therefore, RNA-seq was used to analyze the gene expression in hypothalamic tissues of goats in high fertility group (HFG: litter size per litter ≥2) and low fertility group (LFG: litter size per litter = 1), and identified the differential lncRNAs and mRNAs and their associated pathways related to their fertility. The results showed that a total of 23 lncRNAs and 57 mRNAs were differentially expressed in the hypothalamic tissue of high and low fertility goats. GO terms and KEGG functional annotation suggest that DE lncRNAs and DE mRNAs were significantly enriched in hormone-related pathways regulating ovarian development, hormone synthesis and secretion, regulation of reproductive processes, Estrogen signaling pathway, Oxytocin signaling pathway and GnRH signaling pathway. And we constructed a co-expression network of lncRNAs and target genes, and identified reproduction-related genes such as NMUR2, FEZF1, and WT1. The sequencing results of the hypothalamic transcriptome have broadened our understanding of lncRNA and mRNA in goat hypothalamic tissue and provided some new insights into the molecular mechanisms of follicle development and regulation of its fertility in goats.

Population Structure and Selection Signatures Underlying Domestication Inferred from Genome-Wide Copy Number Variations in Chinese Indigenous Pigs

Single nucleotide polymorphism was widely used to perform genetic and evolution research in pigs. However, little is known about the effect of copy number variation (CNV) on characteristics in pigs. This study performed a genome-wide comparison of CNVs between Wannan black pigs (WBP) and Asian wild boars (AWB), using whole genome resequencing data. By using Manta, we detected in total 28,720 CNVs that covered approximately 1.98% of the pig genome length. We identified 288 selected CNVs (top 1%) by performing Fst statistics. Functional enrichment analyses for genes located in selected CNVs were found to be muscle related (NDN, TMOD4, SFRP1, and SMYD3), reproduction related (GJA1, CYP26B1, WNT5A, SRD5A2, PTPN11, SPEF2, and CCNB1), residual feed intake (RFI) related (MAP3K5), and ear size related (WIF1). This study provides essential information on selected CNVs in Wannan black pigs for further research on the genetic basis of the complex phenotypic and provides essential information for direction in the protection and utilization of Wannan black pig.

Integrated Analysis of mRNAs and Long Non-Coding RNAs Expression of Oviduct That Provides Novel Insights into the Prolificacy Mechanism of Goat (Capra hircus)

Artificial directional selection has replaced natural selection and resulted in trait differences across breeds in domestic animal breeding. However, the molecular mechanism by which the oviduct regulates litter size remains largely elusive in goats during the follicular phase. Accumulating data have linked lncRNAs to reproductive activities; however, little is known about the modulation mechanism in the oviduct. Herein, RNA-seq was used to measure mRNA and lncRNA expression levels in low- and high-fecundity goats. We observed distinctive differences in mRNA and lncRNA in terms of different kidding numbers and detected the differential expression of 1640 mRNA transcripts and 271 lncRNA transcripts. Enrichment analysis of differentially expressed mRNAs (DEGs) suggested that multiple pathways, such as the AMPK, PI3K–Akt, calcium signaling pathway, oocyte meiosis, ABC transporter, and ECM–receptor interaction pathways, directly or indirectly affected goat reproduction. Additionally, coexpression of differentially expressed lncRNAs (DEL)-genes analysis showed that XLOC_021615, XLOC_119780, and XLOC_076450 were trans-acting as the DEGs ATAD2, DEPDC5, and TRPM6, respectively, and could regulate embryo development. Moreover, XLOC_020079, XLOC_107361, XLOC_169844, XLOC_252348 were the trans-regulated elements of the DEGs ARHGEF2 and RAPGEF6, and the target DEGs CPEB3 of XLOC_089239, XLOC_090063, XLOC_107409, XLOC_153574, XLOC_211271, XLOC_251687 were associated with prolificacy. Collectively, our study has offered a thorough dissection of the oviduct lncRNA and mRNA landscapes in goats. These results could serve as potential targets of the oviduct affecting fertility in goats.

Transcription Analysis for Core Networks of lncRNAs–mRNAs: Implication for Potential Role in Sterility of Crassostrea gigas

Simple Summary

This study reveals the expression profiles of lncRNA in the gonads of the Pacific oyster Crassostrea gigas. The potential function of lncRNAs was predicted in the case of antisense and cis-regulatory mechanisms based on their physical positions and their coexpression relationships in the case of trans regulation. Sterility-related DEGs and DELs were chosen for subsequent analysis, demonstrating that trans-regulatory lncRNAs might play a vital role in the gametogenesis of C. gigas. We constructed core networks of lncRNAs–mRNAs for triploid sterile females and hermaphrodites based on pathway results, in which 28 lncRNAs and their 54 trans-regulatory genes were detected. Among 28 sterility-specific lncRNAs, MSTRG.79882.3 and MSTRG.79882.4 for triploid sterile females and MSTRG.33704.1, MSTRG.63844.1, and MSTRG.5675.1 for hermaphrodites play the most significant role.

Abstract

Long noncoding RNA (lncRNA), a type of non-protein-coding transcript, is emerging as a crucial regulator of gene expression. However, few roles of lncRNA in the reproductive process of the Pacific oyster (Crassostrea gigas) have been defined, especially in the regulatory mechanism of sterile triploids gametogenesis. To uncover the potential role of lncRNA, the gonads of diploids, sterile triploids, and partially sterile triploids underwent RNA sequencing. A total of 9618 reliable lncRNAs were identified. The target relationship between lncRNA and mRNA was predicted based on cis, trans, and antisense regulation with bioinformatic software. We chose differentially expressed lncRNAs and mRNAs when sterile triploids were compared to partially sterile triploids and diploids for subsequent functional enrichment analysis. Findings revealed that trans-regulatory lncRNAs might play a significant role in the gametogenesis of C. gigas. Combining pathway results, we constructed core networks of lncRNAs–mRNAs for triploid sterile females and hermaphrodites. Fifty-four genes related to cell division, germline-cell maintenance, and glycogen metabolism were found to be associated with sterility. A total of 28 candidate lncRNAs were predicted to trans-regulate these genes. We speculated that MSTRG.79882.3 and MSTRG.79882.4 for triploid sterile females and MSTRG.33704.1, MSTRG.63844.1, and MSTRG.5675.1 for hermaphrodites were highly important as they were predicted to regulate more sterility-specific genes than others. Our work collectively identified sterility-related lncRNAs and implicated the potential mechanism of lncRNA-mediated regulation in the gametogenesis of sterile triploid oysters.

Brain and Pituitary Transcriptome Analyses Reveal the Differential Regulation of Reproduction-Related LncRNAs and mRNAs in Cynoglossus semilaevis

Long noncoding RNAs (lncRNAs) have been identified to be involved in half-smooth tongue sole (Cynoglossus semilaevis) reproduction. However, studies of their roles in reproduction have focused mainly on the ovary, and their expression patterns and potential roles in the brain and pituitary are unclear. Thus, to explore the mRNAs and lncRNAs that are closely associated with reproduction in the brain and pituitary, we collected tongue sole brain and pituitary tissues at three stages for RNA sequencing (RNA-seq), the 5,135 and 5,630 differentially expressed (DE) mRNAs and 378 and 532 DE lncRNAs were identified in the brain and pituitary, respectively. The RNA-seq results were verified by RT-qPCR. Moreover, enrichment analyses were performed to analyze the functions of DE mRNAs and lncRNAs. Interestingly, their involvement in pathways related to metabolism, signal transduction and endocrine signaling was revealed. LncRNA-target gene interaction networks were constructed based on antisense, cis and trans regulatory mechanisms. Moreover, we constructed competing endogenous RNA (ceRNA) networks. In summary, this study provides mRNA and lncRNA expression profiles in the brain and pituitary to understand the molecular mechanisms regulating tongue sole reproduction.

LncRNAs in domesticated animals: from dog to livestock species

Animal genomes are pervasively transcribed into multiple RNA molecules, of which many will not be translated into proteins. One major component of this transcribed non-coding genome is the long non-coding RNAs (lncRNAs), which are defined as transcripts longer than 200 nucleotides with low coding-potential capabilities. Domestic animals constitute a unique resource for studying the genetic and epigenetic basis of phenotypic variations involving protein-coding and non-coding RNAs, such as lncRNAs. This review presents the current knowledge regarding transcriptome-based catalogues of lncRNAs in major domesticated animals (pets and livestock species), covering a broad phylogenetic scale (from dogs to chicken), and in comparison with human and mouse lncRNA catalogues. Furthermore, we describe different methods to extract known or discover novel lncRNAs and explore comparative genomics approaches to strengthen the annotation of lncRNAs. We then detail different strategies contributing to a better understanding of lncRNA functions, from genetic studies such as GWAS to molecular biology experiments and give some case examples in domestic animals. Finally, we discuss the limitations of current lncRNA annotations and suggest research directions to improve them and their functional characterisation.

Whole Transcriptome Analysis: Implication to Estrous Cycle Regulation

Simple Summary

The databases of mRNA and non-coding-RNAs (miRNA, circRNA, lncRNA) in the ovary of Xinong Sannen goat were reported in this study. The differential expression of mRNA and non-coding RNAs were analyzed, and the comprehensive analysis of the four databases provided RNA networks that regulate estrous cycle, which is essential to improve reproduction.

Abstract

Estrous cycle is one of the placental mammal characteristics after sexual maturity, including estrus stage (ES) and diestrus stage (DS). Estrous cycle is important in female physiology and its disorder may lead to diseases, such as polycystic ovary syndrome, ovarian carcinoma, anxiety, and epilepsy. In the latest years, effects of non-coding RNAs and messenger RNA (mRNA) on estrous cycle have started to arouse much concern, however, a whole transcriptome analysis among non-coding RNAs and mRNA has not been reported. Here, we report a whole transcriptome analysis of goat ovary in estrus and diestrus periods. Estrus synchronization was conducted to induce the estrus phase and on day 32, the goats shifted into the diestrus stage. The ovary RNA of estrus and diestrus stages was respectively collected to perform RNA-sequencing. Then, the circular RNA (circRNA), microRNA (miRNA), long non-coding RNA (lncRNA), and mRNA databases of goat ovary were acquired, and the differential expressions between estrus and diestrus stages were screened to construct circRNA-miRNA-mRNA/lncRNA and lncRNA-miRNA/mRNA networks, thus providing potential pathways that are involved in the regulation of estrous cycle. Differentially expressed mRNAs, such as MMP9, TIMP1, 3BHSD, and PTGIS, and differentially expressed miRNAs that play key roles in the regulation of estrous cycle, such as miR-21-3p, miR-202-3p, and miR-223-3p, were extracted from the network. Our data provided the miRNA, circRNA, lncRNA, and mRNA databases of goat ovary and each differentially expressed profile between ES and DS. Networks among differentially expressed miRNAs, circRNAs, lncRNAs, and mRNAs were constructed to provide valuable resources for the study of estrous cycle and related diseases.

Integrated lncRNA and mRNA Transcriptome Analyses in the Ovary of Cynoglossus semilaevis Reveal Genes and Pathways Potentially Involved in Reproduction

Long non-coding RNAs (lncRNAs) have been reported to be involved in multiple biological processes. However, the roles of lncRNAs in the reproduction of half-smooth tongue sole (Cynoglossus semilaevis) are unclear, especially in the molecular regulatory mechanism driving ovarian development and ovulation. Thus, to explore the mRNA and lncRNA mechanisms regulating reproduction, we collected tongue sole ovaries in three stages for RNA sequencing. In stage IV vs. V, we identified 312 differentially expressed (DE) mRNAs and 58 DE lncRNAs. In stage V vs. VI, we identified 1,059 DE mRNAs and 187 DE lncRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed that DE mRNAs were enriched in ECM-receptor interaction, oocyte meiosis and steroid hormone biosynthesis pathways. Furthermore, we carried out gene set enrichment analysis (GSEA) to identify potential reproduction related-pathways additionally, such as fatty metabolism and retinol metabolism. Based on enrichment analysis, DE mRNAs with a potential role in reproduction were selected and classified into six categories, including signal transduction, cell growth and death, immune response, metabolism, transport and catabolism, and cell junction. The interactions of DE lncRNAs and mRNAs were predicted according to antisense, cis-, and trans-regulatory mechanisms. We constructed a competing endogenous RNA (ceRNA) network. Several lncRNAs were predicted to regulate genes related to reproduction including cyp17a1, cyp19a1, mmp14, pgr, and hsd17b1. The functional enrichment analysis of these target genes of lncRNAs revealed that they were involved in several signaling pathways, such as the TGF-beta, Wnt signaling, and MAPK signaling pathways and reproduction related-pathways such as the progesterone-mediated oocyte maturation, oocyte meiosis, and GnRH signaling pathway. RT-qPCR analysis showed that two lncRNAs (XR_522278.2 and XR_522171.2) were mainly expressed in the ovary. Dual-fluorescence in situ hybridization experiments showed that both XR_522278.2 and XR_522171.2 colocalized with their target genes cyp17a1 and cyp19a1, respectively, in the follicular cell layer. The results further demonstrated that lncRNAs might be involved in the biological processes by modulating gene expression. Taken together, this study provides lncRNA profiles in the ovary of tongue sole and further insight into the role of lncRNA involvement in regulating reproduction in tongue sole.

A Preliminary Study on the Characteristics of microRNAs in Ovarian Stroma and Follicles of Chuanzhong Black Goat during Estrus

microRNAs (miRNAs) play a significant role in ovarian follicular maturity, but miRNA expression patterns in ovarian stroma (OS), large follicles (LF), and small follicles (SF) have been rarely explored. We herein aimed to identify miRNAs, their target genes and signaling pathways, as well as their interaction networks in OS, LF, and SF of Chuanzhong black goats at the estrus phase using small RNA-sequencing. We found that the miRNA expression profiles of LF and SF were more similar than those of OS—32, 16, and 29 differentially expressed miRNAs were identified in OS vs. LF, OS vs. SF, and LF vs. SF, respectively. Analyses of functional enrichment and the miRNA-targeted gene interaction network suggested that miR-182 (SMC3), miR-122 (SGO1), and miR-206 (AURKA) were involved in ovarian organogenesis and hormone secretion by oocyte meiosis. Furthermore, miR-202-5p (EREG) and miR-485-3p (FLT3) were involved in follicular maturation through the MAPK signaling pathway, and miR-2404 (BMP7 and CDKN1C) played a key role in follicular development through the TGF-β signaling pathway and cell cycle; nevertheless, further research is warranted. To our knowledge, this is the first study to investigate miRNA expression patterns in OS, LF, and SF of Chuanzhong black goats during estrus. Our findings provide a theoretical basis to elucidate the role of miRNAs in follicular maturation. These key miRNAs might provide candidate biomarkers for the diagnosis of follicular maturation and will assist in developing new therapeutic targets for female goat infertility.