Filtered reproductive long non-coding RNAs by genome-wide analyses of goat ovary at different estrus periods

Multi-tissue transcriptomic characterization of endogenous retrovirus-derived transcripts in Capra hircus

Background

Transposable elements (TEs, or transposons) are repetitive genomic sequences, accounting for half of a mammal genome. Most TEs are transcriptionally silenced, whereas some TEs, especially endogenous retroviruses (ERVs, long terminal repeat retrotransposons), are physiologically expressed in certain conditions. However, the expression pattern of TEs in those less studied species, like goat (Capra hircus), remains unclear. To obtain an overview of the genomic and transcriptomic features of TEs and ERVs in goat, an important farm species, we herein analyzed transcriptomes of ten C. hircus tissues and cells under various physiological and pathological conditions.

Method

Distribution of classes, families, and subfamilies of TEs in the C. hircus genome were systematically annotated. The expression patterns of TE-derived transcripts in multiple tissues were investigated at subfamily and location levels. Differential expression of ERV-derived reads was measured under various physiological and pathological conditions, such as embryo development and virus infection challenges. Co-expression between ERV-reads and their proximal genes was also explored to decipher the expression regulation of ERV-derived transcripts.

Results

There are around 800 TE subfamilies in the goat genome, accounting for 49.1% of the goat genome sequence. TE-derived reads account for 10% of the transcriptome and their abundance are comparable in various goat tissues, while expression of ERVs are variable among tissues. We further characterized expression pattern of ERV reads in various tissues. Differential expression analysis showed that ERVs are highly active in 16-cell embryos, when the genome of the zygote begins to transcribe its own genes. We also recognized numerous activated ERV reads in response to RNA virus infection in lung, spleen, caecum, and immune cells. CapAeg_1.233:ERVK in chromosome 1 and 17 are dysregulated under endometrium development and infection conditions. They showed strong co-expression with their proximal gene OAS1 and TMPRSS2, indicating the impact of activated proximal gene expression on nearby ERVs.

Conclusion

We generated ERV transcriptomes across goat tissues, and identified ERVs activated in response to different physiological and pathological conditions.

Characterization of MicroRNA expression profiles in the ovarian tissue of goats during the sexual maturity period

BACKGROUND

The ovary is an important reproductive organ in mammals, and its development directly affects the sexual maturity and reproductive capacity of individuals. MicroRNAs (miRNAs) are recognized as regulators of reproductive physiological processes in various animals and have been shown to regulate ovarian development through typical targeting and translational repression. However, little is known about the regulatory role of miRNAs in ovarian tissue development during sexual maturity in goats. To comprehensively profile the different physiological stages of sexual maturation in goats, we performed small-RNA sequencing of ovarian tissue samples collected at four specific time points (1 day after birth (D1), 2 months old (M2), 4 months old (M4), and 6 months old (M6)). In addition, we used ELISAs to measure serum levels of reproductive hormones to study their temporal changes.

RESULTS

The results showed that serum levels of gonadotropin-releasing hormone, follicle-stimulating hormone, luteinizing hormone, oestradiol, progesterone, oxytocin, and prolactin were lower in goats at the D1 stage than at the other three developmental stages (P < 0.05). The secretion patterns of these seven hormones show a similar trend, with hormone levels reaching their peaks at 4 months of age. A total of 667 miRNAs were detected in 20 libraries, and 254 differentially expressed miRNAs and 3 groups of miRNA clusters that had unique expression patterns were identified (|log2-fold change|> 1, FDR < 0.05) in the 6 comparison groups. RT‒qPCR was employed to confirm that the expression pattern of the 15 selected miRNAs was consistent with the Illumina sequencing results. Gene ontology analyses revealed significant enrichment of GO terms such as cell proliferation regulation, epithelial cell development, and amino acid transport, as well as important signaling pathways including the MAPK signaling pathway, the PI3K-Akt signaling pathway, and the oestrogen signaling pathway. Further miRNA‒mRNA regulation network analysis revealed that 8 differentially expressed miRNAs (chi-miR-1343, chi-miR-328-3p, chi-miR-877-3p, chi-miR-296-3p, chi-miR-128-5p, chi-miR-331-3p, chi-miR-342-5p and chi-miR-34a) have important regulatory roles in ovarian cell proliferation, hormone secretion and metabolism-related biological processes.

CONCLUSIONS

Overall, our study investigated the changes in serum hormone and miRNA levels in the ovaries. These data provide a valuable resource for understanding the molecular regulatory mechanisms of miRNAs in ovarian tissue during the sexual maturity period in goats.

Analysis of serum reproductive hormones and ovarian genes in pubertal female goats

BACKGROUND

Age at puberty is an important factor affecting goat fertility, with endocrine and genetic factors playing a crucial role in the onset of puberty. To better understand the relationship between endocrine and genetic factors and mechanisms underlying puberty onset in goats, reproductive hormone levels were analyzed by ELISA and ultraperformance liquid chromatography-multiple reaction monitoring-multistage/mass spectrometry and RNA sequencing was performed to analyze ovarian genes.

RESULTS

Serum follicle stimulating hormone, luteinizing hormone, estradiol, 11-deoxycortisol, 11-deoxycorticosterone, corticosterone, cortisone, and cortisol levels were found to be higher but progesterone were lower in pubertal goats as compared to those in prepubertal goats (P < 0.05). A total of 18,139 genes were identified in cDNA libraries, and 75 differentially expressed genes (DEGs) were identified (|log₂ fold change|≥ 1, P ≤ 0.05), of which 32 were significantly up- and 43 were down-regulated in pubertal goats. Gene ontology enrichment analyses indicated that DEGs were mainly involved in "metabolic process," "signaling," "reproduction," and "growth." Further, DEGs were significantly enriched in 91 Kyoto Encyclopedia of Genes and Genomes pathways, including estrogen signaling pathway, steroid hormone biosynthesis, and cAMP signaling pathway. Bioinformatics analysis showed that PRLR and THBS1 were highly expressed in pubertal ovaries, and ZP3, ZP4, and ASTL showed low expression, suggesting their involvement in follicular development and lutealization.

CONCLUSIONS

To summarize, serum hormone changes and ovarian DEGs expression were investigated in our study. Further studies are warranted to comprehensively explore the functions of DEGs in goat puberty.

Oviduct Transcriptomic Reveals the Regulation of mRNAs and lncRNAs Related to Goat Prolificacy in the Luteal Phase

Simple Summary

The kidding number is an important reproductive trait in domestic goats. The oviduct, as one of the most major organs, is directly involved in the reproductive process, providing nutrition and a location for early embryonic development. The current study provides genome-wide expression profiles of mRNA and long noncoding RNAs (lncRNAs) expression in Yunshang black goat, a new breed of meat goat bred in China with a high kidding number. During the luteal phases, oviduct mRNAs and lncRNAs associated with high- and low-fecundity Yunshang black goats were identified, and their potential biological functions were predicted using GO, KEGG, and GSEA enrichment analysis. These findings shed light on the oviduct-based prolificacy mechanism in goats.

Abstract

The oviduct is associated with embryo development and transportation and regulates the pregnancy success of mammals. Previous studies have indicated a molecular mechanism of lncRNAs in gene regulation and reproduction. However, little is known about the function of lncRNAs in the oviduct in modulating goat kidding numbers. Therefore, we combined RNA sequencing (RNA-seq) to map the expression profiles of the oviduct at the luteal phase from high- and low-fecundity goats. The results showed that 2023 differentially expressed mRNAs (DEGs) and 377 differentially expressed lncRNAs (DELs) transcripts were screened, and 2109 regulated lncRNA-mRNA pairs were identified. Subsequently, the genes related to reproduction (IGF1, FGFRL1, and CREB1) and those associated with embryonic development and maturation (DHX34, LHX6) were identified. KEGG analysis of the DEGs revealed that the GnRH- and prolactin-signaling pathways, progesterone-mediated oocyte maturation, and oocyte meiosis were related to reproduction. GSEA and KEGG analyses of the target genes of DELs demonstrated that several biological processes and pathways might interact with oviduct functions and the prolificacy of goats. Furthermore, the co-expression network analysis showed that XLOC_029185, XLOC_040647, and XLOC_090025 were the cis-regulatory elements of the DEGs MUC1, PPP1R9A, and ALDOB, respectively; these factors might be associated with the success of pregnancy and glucolipid metabolism. In addition, the GATA4, LAMA2, SLC39A5, and S100G were trans-regulated by lncRNAs, predominantly mediating oviductal transport to the embryo and energy metabolism. Our findings could pave the way for a better understanding of the roles of mRNAs and lncRNAs in fecundity-related oviduct function in goats.

Identification and Characterization of lncRNAs Expression Profile Related to Goat Skeletal Muscle at Different Development Stages

LncRNAs are essential for regulating skeletal muscle. However, the expression profile and function of lncRNAs in goat muscle remains unclear. Here, an average of ~14.58 Gb high-quality reads were obtained from longissimus dorsi tissues of 1-month-old (n = 3) and 9-month-old (n = 3) Wu'an black goats using RNA sequencing. Of a total of 3441 lncRNAs, 1281 were lincRNAs, 805 were antisense lncRNAs, and 1355 were sense_overlapping lncRNAs. These lncRNAs shared some properties with goats, such as fewer exons, shorter transcript, and open reading frames (ORFs) length. Among them, 36 differentially expressed lncRNAs (DE lncRNA) were identified, and then 10 random lncRNAs were validated by RT-qPCR. Furthermore, 30 DE lncRNAs were neighboring 71 mRNAs and several genes were functionally enriched in muscle development-related pathways, such as APC, IFRD1, NKX2-5, and others. Additionally, 36 DE lncRNAs and 2684 mRNAs were included in co-expression interactions. A lncRNA-miRNA-mRNA network containing 4 lncRNAs, 3 miRNAs, and 8 mRNAs was finally constructed, of which XR_001296113.2 might regulate PDLIM7 expression by interaction with chi-miR-1296 to affect skeletal muscle development. This study revealed the expression profile of goat lncRNAs for further investigative studies and provides a fuller understanding of skeletal muscle development.

Hand2os1 Regulates the Secretion of Progesterone in Mice Corpus Luteum

The corpus luteum plays a key role in pregnancy maintenance and estrous cycle regulation by secreting progesterone. Hand2os1 is an lncRNA located upstream of Hand2, with which a bidirectional promoter is shared and is involved in the regulation of cardiac development and embryo implantation in mice. The aim of this study was to investigate the expression and regulation of Hand2os1 in the ovaries. Here, we used RNAscope to detect differential expression of Hand2os1 in the ovaries of cycling and pregnant mice. Hand2os1 was specifically detected in luteal cells during the proestrus and estrus phases, showing its highest expression in the corpus luteum at estrus. Additionally, Hand2os1 was strongly expressed in the corpus luteum on day 4 of pregnancy, but the positive signal progressively disappeared after day 8, was detected again on day 18, and gradually decreased after delivery. Hand2os1 significantly promoted the synthesis of progesterone and the expression of StAR and Cyp11a1. The decreased progesterone levels caused by Hand2os1 interference were rescued by the overexpression of StAR. Our findings suggest that Hand2os1 may regulate the secretion of progesterone in the mouse corpus luteum by affecting the key rate-limiting enzyme StAR, which may have an impact on the maintenance of pregnancy.

Quantitative proteomics analysis to assess protein expression levels in the ovaries of pubescent goats

Background

Changes in the abundance of ovarian proteins play a key role in the regulation of reproduction. However, to date, no studies have investigated such changes in pubescent goats. Herein we applied isobaric tags for relative and absolute quantitation (iTRAQ) and liquid chromatography–tandem mass spectrometry to analyze the expression levels of ovarian proteins in pre-pubertal (n = 3) and pubertal (n = 3) goats.

Results

Overall, 7,550 proteins were recognized; 301 (176 up- and 125 downregulated) were identified as differentially abundant proteins (DAPs). Five DAPs were randomly selected for expression level validation by Western blotting; the results of Western blotting and iTRAQ analysis were consistent. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis indicated that DAPs were enriched in olfactory transduction, glutathione metabolism, and calcium signaling pathways. Besides, gene ontology functional enrichment analysis revealed that several DAPs enriched in biological processes were associated with cellular process, biological regulation, metabolic process, and response to stimulus. Protein–protein interaction network showed that proteins interacting with CDK1, HSPA1A, and UCK2 were the most abundant.

Conclusions

We identified 301 DAPs, which were enriched in olfactory transduction, glutathione metabolism, and calcium signaling pathways, suggesting the involvement of these processes in the onset of puberty. Further studies are warranted to more comprehensively explore the function of the identified DAPs and aforementioned signaling pathways to gain novel, deeper insights into the mechanisms underlying the onset of puberty.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08699-y.

Identification, characterization and functional analysis of gonadal long noncoding RNAs in a protogynous hermaphroditic teleost fish, the ricefield eel (Monopterus albus)

Background

An increasing number of long noncoding RNAs (lncRNAs) have been found to play important roles in sex differentiation and gonad development by regulating gene expression at the epigenetic, transcriptional and posttranscriptional levels. The ricefield eel, Monopterus albus, is a protogynous hermaphroditic fish that undergoes a sequential sex change from female to male. However, the roles of lncRNA in the sex change is unclear.

Results

Herein, we performed RNA sequencing to analyse lncRNA expression patterns in five different stages of M. albus development to investigate the roles of lncRNAs in the sex change process. A total of 12,746 lncRNAs (1503 known lncRNAs and 11,243 new lncRNAs) and 2901 differentially expressed lncRNAs (DE-lncRNAs) were identified in the gonads. The target genes of the DE-lncRNAs included foxo1, foxm1, smad3, foxr1, camk4, ar and tgfb3, which were mainly enriched in signalling pathways related to gonadal development, such as the insulin signalling pathway, MAPK signalling pathway, and calcium signalling pathway. We selected 5 highly expressed DE-lncRNAs (LOC109952131, LOC109953466, LOC109954337, LOC109954360 and LOC109958454) for full length amplification and expression pattern verification. They were all expressed at higher levels in ovaries and intersex gonads than in testes, and exhibited specific time-dependent expression in ovarian tissue incubated with follicle-stimulating hormone (FSH) and human chorionic gonadotropin (hCG). The results of quantitative real-time PCR (qRT-PCR) analysis and a dual-luciferase assay showed that znf207, as the gene targeted by LOC109958454, was expressed in multiple tissues and gonadal developmental stages of M. albus, and its expression was also inhibited by the hormones FSH and hCG.

Conclusions

These results provide new insights into the role of lncRNAs in gonad development, especially regarding natural sex changes in fish, which will be useful for enhancing our understanding of sequential hermaphroditism and sex changes in the ricefield eel (M. albus) and other teleosts.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08679-2.

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.

Expression Profile Analysis to Identify Circular RNA Expression Signatures in the Prolificacy Trait of Yunshang Black Goat Pituitary in the Estrus Cycle

The pituitary gland is an important organ. It is a complex area of the brain involved in endocrine function and reproductive regulation. However, the function of the pituitary in goat reproduction is still unclear. Herein, RNA sequencing was used to explore the expression patterns of circle RNAs (circRNAs) in the pituitary of Yunshang black goats during the various estrus phases. Then the host genes of the circRNAs were predicted, and a competing endogenous RNA (ceRNA) network was constructed. The results showed a total of 6,705 circRNAs in the pituitary of Yunshang black goats, among which 388 differentially expressed (DE) circRNAs (214 were upregulated, while 174 were downregulated) were identified between high- and low-yield Yunshang black goats in the follicular phase (HF vs. LF); moreover, 361 DE circRNAs (136 were upregulated, while 225 were downregulated) were identified between high- and low-yield Yunshang black goats in the luteal phase (HL vs. LL). There were 65 DE circRNAs targeting 40 miRNAs in the HF vs. LF comparison and 46 DE circRNAs targeting 31 miRNAs in the HL vs. LL comparison. We identified chi_circ_0030920, chi_circ_0043017, chi_circ_0008353, chi_circ_0041580, and chi_circ_0016478 as the key circRNAs through functional enrichment analysis. The ceRNA network analysis showed that chi_circ_0031209 and chi_circ_0019448 might play an important role in reproduction by influencing the expression of prolactin receptor (PRLR) in high- and low-yielding goats during the luteal phase, whereas chi_circ_0014542 regulates the expression of WNT5A during the follicular phase. Our study provided the overall expression profiles of circRNAs in the goat pituitary during the estrus phase, which provides new insight into the mechanism of high-yield goats, which can be helpful to guide goat breeding.

NORHA, a novel follicular atresia-related lncRNA, promotes porcine granulosa cell apoptosis via the miR-183-96-182 cluster and FoxO1 axis

BACKGROUND

Follicular atresia has been shown to be strongly associated with a low follicle utilization rate and female infertility, which are regulated by many factors such as microRNAs (miRNAs), which constitute a class of noncoding RNAs (ncRNAs). However, little is known about long noncoding RNAs (lncRNAs), which constitute another ncRNA family that regulate follicular atresia.

RESULTS

A total of 77 differentially expressed lncRNAs, including 67 upregulated and 10 downregulated lncRNAs, were identified in early atretic follicles compared to healthy follicles by RNA-Sequencing. We characterized a noncoding RNA that was highly expressed in atretic follicles (NORHA). As an intergenic lncRNA, NORHA was one of the upregulated lncRNAs identified in the atretic follicles. To determine NORHA function, RT-PCR, flow cytometry and western blotting were performed, and the results showed that NORHA was involved in follicular atresia by influencing GC apoptosis with or without oxidative stress. To determine the mechanism of action, bioinformatics analysis, luciferase reporter assay and RNA immunoprecipitation assay were performed, and the results showed that NORHA acted as a 'sponge', that directly bound to the miR-183-96-182 cluster, and thus prevented its targeted inhibition of FoxO1, a major sensor and effector of oxidative stress.

CONCLUSIONS

We provide a comprehensive perspective of lncRNA regulation of follicular atresia, and demonstrate that NORHA, a novel lncRNA related to follicular atresia, induces GC apoptosis by influencing the activities of the miR-183-96-182 cluster and FoxO1 axis.

Analysis of Transcripts of Uncertain Coding Potential Using RNA Sequencing During the Preattachment Phase in Goat Endometrium

Transcripts of uncertain coding potential (TUCP) are part of long noncoding RNAs, which include short open reading frames and could be translated into small peptides. In recent years, a growing number of TUCPs has been implicated in multiple biological activities, such as embryogenesis and transcriptional regulation. However, the abundance of TUCPs and their roles in goat endometrium during pregnancy recognition (day 16) remain undocumented. In this study, bioinformatics analyses were conducted to identify the differentially expressed (DE) TUCPs between pregnant animals and corresponding nonpregnant controls. A total of 5551 TUCPs were identified; 114 TUCPs were DE in goat endometrium, of which 74 TUCPs were upregulated in pregnant endometrium, whereas 40 TUCPs were downregulated. The related genes of TUCP were predicted by using coexpression and colocalization methods. In summary, 419 genes were predicted by colocalization, and 9464 genes were predicted by coexpression. The kyoto encyclopedia of genes and genomes (KEGG) and gene ontology (GO) analysis showed that TUCPs, which are highly expressed in pregnant endometrium, were mainly associated with endometrial remodeling, nutrient synthesis, and transportation. However, TUCPs that were lowly expressed in pregnant endometrium were mainly associated with immune tolerance, which is necessary for the protection and development of the embryo in the uterus. These findings may be used for the comparative analysis of TUCP transcripts in endometrium and assist in the selection of applicable candidate genes associated with embryo implantation for further functional analyses.

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.

Transcriptome Analysis of Long Noncoding RNAs and mRNAs in Granulosa Cells of Jinghai Yellow Chickens Illuminated With Red Light

Jinghai Yellow chickens are a new indigenous breed with a dual purpose in China, but their egg laying performance is limited. Compared with white light (WL), exposure to red light (RL) can improve the egg laying performance of hens. Herein, to elucidate the molecular mechanism by which RL affects the egg laying performance, RNA sequencing was used to analyze long noncoding RNAs (lncRNAs) and mRNAs from granulosa cells of small yellow follicles from Jinghai Yellow chickens in RL and WL groups. A total of 12,466 lncRNAs were identified among the assembled transcripts, of which 168 lncRNAs were significantly different between the RL and WL groups (101 downregulated and 67 upregulated). Additionally, 1182 differentially expressed mRNAs were identified (958 downregulated and 224 upregulated). Integrated network analysis demonstrated that numerous differential mRNAs were involved in follicular development through steroid hormone synthesis, oocyte meiosis, and the PI3K-Akt signaling pathway. The impact of lncRNAs on cis and trans target mRNAs indicates that some lncRNAs play important roles in follicular development of small yellow follicles. The results provide a starting point for studies aimed at understanding the molecular mechanisms by which monochromatic light affects follicular development and egg production in hens.

Differential expression and prediction of function of lncRNAs in the ovaries of low and high fecundity Hanper sheep

Litter size is an important trait that determines the production efficiency of sheep bred for meat. Its detailed investigation can reveal the molecular mechanisms that control the fecundity of sheep and possibly accelerate the breeding process of new varieties of sheep that have high prolificacy. Long non-coding RNAs (lncRNAs) have proven to be an important factor in the regulation of follicular development. However, the mechanisms by which lncRNAs regulate litter size in sheep remain unclear. In the present study, ovarian tissues from the follicular (F) or luteal phase (L) of Hanper sheep that were either monotocous (M) or polytocous (P; FM, FP, LM and LP groups) were collected and sequenced to identify differentially expressed lncRNAs and predict their function. The results indicate that the number of up- and down-regulated lncRNAs in the follicular phase (FM vs. FP) was 95 and 111 and 109 and 49, respectively, in the luteal phase (LM vs. LP). The functional enrichment of the different lncRNAs coexpressed with mRNA was analysed. The results demonstrated that the KISS1-GnRH-LH/FSH-E2 and EGF-EGFR-RAS-PI3K signalling pathways promoted the initiation of the primordial period, follicular development and ovulation in the follicular phase (FM vs. FP). During the luteal phase (LM vs. LP), the production and development of the corpus luteum in ewes was influenced by the KITLG-KIT/FGF-FGFR/HGF-MET-RAS-ERK signalling pathway. STEM clustering functional enrichment analysis of the differentially expressed lncRNAs indicated that profile11 was principally enriched in the Cytokine-Jak-STAT, PDGF-PDGFR-PI3K and KITLG-KIT-RAS-ERK signalling pathways. By analysis of the differential expression of the lncRNAs and their expression in each group, lncRNAs Xist (loc101112291) and Gtl2 (loc101123329) were found to be highly expressed, suggesting that regulation of follicular development was mediated through methylation processes.

Transcriptome analysis of long noncoding RNAs ribonucleic acids from the livers of Hu sheep with different residual feed intake

Long noncoding RNAs (LncRNAs), as key regulators, have vital functions in various biological activities. However, in sheep, little has been reported concerning the genetic mechanism of LncRNA regulation of feed efficiency. In the present study, we explored the genome-wide expression of LncRNAs and transcripts of uncertain coding potential (TUCPs) in the livers of sheep with extreme residual feed intake (RFI) using RNA sequencing. We identified 1 523 TUCPs and 1 996 LncRNAs, among which 10 LncRNAs and 16 TUCPs were identified as being differentially expressed between the High-RFI and Low-RFI groups. Co-expression and co-localization methods were used to search for LncRNA and TUCP target genes, which identified 970/1 538 and 23/27 genes, respectively. Ontology and pathways analysis revealed that the LncRNAs/TUCPs that were highly expressed in the Low-RFI group are mostly concentrated in energy metabolism pathways. For example, LNC_000890 and TUCP_000582 might regulate liver tissue metabolic efficiency. The LncRNAs/TUCPs that were highly expressed in the High-RFI group are mostly enriched in immune function pathways. For example, TUCP_000832 might regulate animal health, thereby affecting feed efficiency. Subsequently, a co-expression network was established by applying the expression information of both the differentially expressed LncRNAs and TUCPs and their target mRNAs. The network indicated that differentially expressed genes targeted by the upregulated LncRNAs and TUCPs were mainly related to energy metabolism, while those genes targeted by the downregulated LncRNAs and TUCPs were mainly related to immune response. These results provide the basis for further study of LncRNA/TUCP-mediated regulation of feed efficiency.

Identification and Expression Analysis of Long Non-coding RNA in Large Yellow Croaker (Larimichthys crocea) in Response to Cryptocaryon irritans Infection

Large-scale transcription studies have revealed numerous lncRNAs (long non-coding RNAs). lncRNAs have been proposed to participate in the regulation of a diverse range of biological processes, including transcriptional regulation. Although lncRNAs have attracted increasing attention, the studies in large yellow croaker (Larimichthys crocea) are still rare, and they lack systematic analysis. In this study, 101 RNA-seq datasets varied in ages, sexes, and tissues were retrieved from the NCBI database to generate a comprehensive catalog of large yellow croaker transcriptome database. A set of 14,599 high-confidence lncRNAs from 13,673 loci were identified and characterized. Furthermore, RNA-seq datasets obtained from the infection of C. irritans were employed to investigate the differential expression pattern of lncRNAs and analyze potential biological functions. A total of 77 differentially expressed lncRNAs targeting to 567 protein-coding genes were identified by using expression analysis. Several immune genes, including TLR5, CD2AP, and MMP9, were highlighted. With GO enrichment and KEGG pathway analysis, the immune-related terms or pathways were enriched. This study created a comprehensive dataset of lncRNAs for large yellow croaker, which would be helpful for the researches of functional roles of lncRNAs in large yellow croaker.

Analysis and Screening of Reproductive Long Non-coding RNAs Through Genome-Wide Analyses of Goat Endometrium During the Pre-attachment Phase

Reproduction in goat is highly impeded by implantation failure. Of concern, the underlying mechanism leading to embryo implantation remains unclear. In this study, deep sequencing was employed through strand-specific Ribo-Zero RNA-Seq to characterize transcriptome changes in the endometrium during the maternal recognition of pregnancy. A total of 996 differential transcripts (115 lncRNAs and 881 mRNAs) existing between the pregnant and non-pregnant endometrium were revealed through bioinformatics analysis. The screening was performed on lncRNAs (XR_001918173.1, LNC_002760, and LNC_000599) and LNC_009053, to determine their potential role in regulating the synthesis of retinol and endometrium remolding through the proteasome pathway, respectively. The hypothesis of whether certain lncRNAs, namely, LNC_007223, LNC_005256, and LNC_010092 could play important roles in embryo implantation was tested. These novel findings are of paramount relevance to further elucidate the molecular mechanisms of embryo implantation and uncover new targets to improve goat reproduction.

Genome-Wide Network of lncRNA-mRNA During Ovine Oocyte Development From Germinal Vesicle to Metaphase II in vitro

Long non-coding RNA (lncRNA) is involved in many biological processes, and it has been closely investigated. However, research into the role of lncRNA in ovine ovarian development is scant and poorly understood, particularly in relation to the molecular mechanisms of ovine oocyte maturation. In the current study, RNA sequencing was performed with germinal vesicle (GV) and in vitro matured metaphase II (MII) stage oocytes, isolated from ewes. Through the use of bioinformatic analysis, abundant candidate lncRNAs in stage-specific ovine oocytes were identified, and their trans- and cis-regulatory effects were deeply dissected using computational prediction software. Functional enrichment analysis of these lncRNAs revealed that they were involved in the regulation of many key signaling pathways during ovine oocyte development, which was reflected by their targeted genes. From this study, multiple lncRNA-mRNA networks were presumed to be involved in key signaling pathways regarding ovine oocyte maturation and meiotic resumption. In particular, one novel lncRNA (MSTRG.17927) appeared to mediate the regulation of phosphatidylinositol 3-kinase signaling (PI3K) signaling during ovine oocyte maturation. Therefore, this research offers novel insights into the molecular mechanisms underlying ovine oocyte meiotic maturation regulated by lncRNA-mRNA networks from a genome-wide perspective.

Comprehensive analysis of mRNAs and miRNAs in the ovarian follicles of uniparous and multiple goats at estrus phase

Background

Fertility is an important economic trait in the production of meat goat, and follicular development plays an important role in fertility. Although many mRNAs and microRNAs (miRNAs) have been found to play critical roles in ovarian biological processes, the interaction between mRNAs and miRNAs in follicular development is not yet completely understood. In addition, less attention has been given to the study of single follicle (dominant or atretic follicle) in goats. This study aimed to identify mRNAs, miRNAs, and signaling pathways as well as their interaction networks in the ovarian follicles (large follicles and small follicles) of uniparous and multiple Chuanzhong black goats at estrus phase using RNA-sequencing (RNA-seq) technique.

Results

The results showed that there was a significant difference in the number of large follicles between uniparous and multiple goats (P < 0.05), but no difference in the number of small follicles was observed (P > 0.05). For the small follicles of uniparous and multiple goats at estrus phase, 289 differentially expressed mRNAs (DEmRNAs) and 16 DEmiRNAs were identified; and for the large follicles, 195 DEmRNAs and 7 DEmiRNAs were identified. The functional enrichment analysis showed that DE genes in small follicles were significantly enriched in ovarian steroidogenesis and steroid hormone biosynthesis, while in large follicles were significantly enriched in ABC transporters and steroid hormone biosynthesis. The results of quantitative real-time polymerase chain reaction were consistent with those of RNA-seq. Analysis of the mRNA-miRNA interaction network suggested that CD36 (miR-122, miR-200a, miR-141), TNFAIP6 (miR-141, miR-200a, miR-182), CYP11A1 (miR-122), SERPINA5 (miR-1, miR-206, miR-133a-3p, miR-133b), and PTGFR (miR-182, miR-122) might be related to fertility, but requires further research on follicular somatic cells.

Conclusions

This study was used for the first time to reveal the DEmRNAs and DEmiRNAs as well as their interaction in the follicles of uniparous and multiple goats at estrus phase using RNA-seq technology. Our findings provide new clues to uncover the molecular mechanisms and signaling networks of goat reproduction that could be potentially used to increase ovulation rate and kidding rate in goat.