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

LncRNA STAT3-AS regulates endometrial receptivity via the STAT3 signaling pathway

Endometrial receptivity is critical for the successful establishment of pregnancy in ruminants. Interferon tau (IFNT) plays a key role in promoting embryo attachment by activating the Janus kinase/signal transducer and activator of transcription pathway, which induces the expression of a series of interferon-stimulated genes (ISGs). In our previous study, sequencing analysis of goat endometrial epithelial cells (gEECs) treated with 20 ng/mL IFNT revealed a differentially expressed long non-coding RNA located on the STAT3 antisense chain, which we designated STAT3-AS. The aim of this study was to investigate the role and mechanism of STAT3-AS in establishing endometrial receptivity in goats. The results showed that STAT3-AS was expressed in both the nucleus and cytoplasm of gEECs, and its expression increased significantly in the uterus on day 15 of pregnancy. STAT3-AS expression was upregulated in gEECs treated with IFNT alone or in combination with progesterone and estradiol. Knockdown of STAT3-AS using specific short interfering RNA significantly inhibited the expression of classical ISGs such as interferon-stimulated gene 15 and 2',5'-oligodenylate synthetase 2, as well as uterine endometrial receptivity-related genes including homeobox gene A11, integrin beta 3, and vascular endothelial growth factor. Moreover, gEEC proliferation and the STAT3 pathway were suppressed in the absence of STAT3-AS. However, pretreatment with the STAT3 activator RO8191 restored the effect of silencing STAT3-AS on endometrial receptivity. Overall, these results suggest that STAT3-AS is an important regulator of endometrial receptivity in goats and that it regulates endometrial receptivity through the STAT3 pathway.

Transcriptomic Analysis of the Developing Testis and Spermatogenesis in Qianbei Ma Goats

Reproductive competence in male mammals depends on testicular function. Testicular development and spermatogenesis in goats involve highly complex physiological processes. In this study, six testes were, respectively, obtained from each age group, immature (1 month), sexually mature (6 months) and physically mature (12 months old) Qianbei Ma goats. RNA-Seq was performed to assess testicular mRNA expression in Qianbei Ma goats at different developmental stages. Totally, 18 libraries were constructed to screen genes and pathways involved in testis development and spermatogenesis. Totally, 9724 upregulated and 4153 downregulated DEGs were found between immature (I) and sexually mature (S) samples; 7 upregulated and 3 downregulated DEGs were found between sexually mature (S) and physically mature (P) samples, and about 4% of the DEGs underwent alternative splicing events between I and S. Select genes were assessed by qRT-PCR, corroborating RNA-Seq findings. The detected genes have key roles in multiple developmental stages of goat testicular development and spermatogenesis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to determine differentially expressed genes (DEGs). GO analysis revealed DEGs between S and P contributed to "reproduction process", "channel activity" and "cell periphery part" between I and S, and in "ion transport process", "channel activity" and "transporter complex part". KEGG analysis suggested the involvement of "glycerolipid metabolism", "steroid hormone biosynthesis" and "MAPK signaling pathway" in testis development and spermatogenesis. Genes including IGF1, TGFB1, TGFBR1 and EGFR may control the development of the testis from immature to sexually mature, which might be important candidate genes for the development of goat testis. The current study provides novel insights into goat testicular development and spermatogenesis.

Uterus proliferative period ceRNA network of Yunshang black goat reveals candidate genes on different kidding number trait

Pregnancy loss that occurs in the uterus is an important and widespread problem in humans and farm animals and is also a key factor affecting the fecundity of livestock. Understanding the differences in the fecundity of goats may be helpful in guiding the breeding of goats with high fecundity. In this study, we performed RNA sequencing and bioinformatics analysis to study the uterus of Yunshang black goats with high and low fecundity in the proliferative period. We identified mRNAs, long non-coding RNAs (lncRNAs), and microRNAs (miRNAs) by analyzing the uterine transcriptomes. The target genes of the identified miRNAs and lncRNAs were predicted, and miRNA-mRNA interaction and competitive endogenous RNA (ceRNA) networks were constructed. By comparisons between low- and high-fecundity groups, we identified 1,674 differentially expressed mRNAs (914 were upregulated, and 760 were downregulated), 288 differentially expressed lncRNAs (149 were upregulated, and 139 were downregulated), and 17 differentially expressed miRNAs (4 were upregulated, and 13 were downregulated). In addition, 49 miRNA-mRNA pairs and 45 miRNA-lncRNA pairs were predicted in the interaction networks. We successfully constructed a ceRNA interaction network with 108 edges that contained 19 miRNAs, 11 mRNAs, and 73 lncRNAs. Five candidate genes (PLEKHA7, FAT2, FN1, SYK, and ITPR2) that were annotated as cell adhesion or calcium membrane channel protein were identified. Our results provide the overall expression profiles of mRNAs, lncRNAs, and miRNAs in the goat uterus during the proliferative period and are a valuable reference for studies into the mechanisms associated with the high fecundity, which may be helpful to guide goat to reduce pregnancy loss.

Interferon-τ -induced ISG15-AS regulates endometrial receptivity during early goat pregnancy

Endometrial receptivity is a critical process for the successful establishment of pregnancy in ruminants. However, the biological role of long non-coding RNAs (lncRNAs) in the development of endometrial receptivity is poorly understood. In this study, we performed RNA-seq analysis of immortalised goat endometrial epithelial cells (gEECs) treated with interferon-τ (IFNT). Transcriptome profiles showed that 8069 high-confidence putative lncRNAs, including 6498 intronic lncRNA transcripts, 1078 lincRNAs and 493 antisense lncRNAs were identified in gEECs with or without IFNT treatment. Functional clustering analysis was performed by using cis and trans lncRNAs prediction. GO and KEGG analyses revealed that differentially expressed lncRNAs may regulate tissue remodelling and immune responses. Subsequently, six of the 21 differentially expressed antisense lncRNAs were validated using qRT-PCR. Through functional screening and co-expression analysis of lncRNAs in gEECs, we identified that ISG15-AS was mainly expressed in the luminal and glandular epithelium on days 5 and 15 and was strongly upregulated on day 18 of pregnancy in vivo. Similarly, ISG15-AS was abundant in the nucleus and cytoplasm, and was significantly upregulated after treatment with IFNT in gEECs. In addition, ISG15 is an IFNT-responsive gene, that displayed an evident increase in vivo and in vitro. Moreover, sense ISG15 was significantly upregulated following ISG15-AS silencing. The key genes related to ISGylation and endometrial receptivity in gEECs dramatically increased after ISG15-AS inhibition. Collectively, our results indicate that a novel antisense lncRNA, ISG15-AS, may be important in regulating endometrial receptivity through ISGylation.

Analyses of lncRNAs, circRNAs, and the Interactions between ncRNAs and mRNAs in Goat Submandibular Glands Reveal Their Potential Function in Immune Regulation

As part of one of the main ruminants, goat salivary glands hardly secrete digestive enzymes, but play an important role in immunity. The immune function of goat salivary glands significantly changes with age, while the expression profile and specific function of non-coding RNA during this process are unknown. In this study, transcriptome sequencing was performed on submandibular gland (SMG) tissues of 1-month-old, 12-month-old, and 24-month-old goats, revealing the expression patterns of lncRNA and circRNA at different ages. A total of 369 lncRNAs and 1699 circRNAs were found to be differentially expressed. Functional enrichment analyses showed that the lncRNA regulated target mRNAs and circRNA host genes were significantly enriched in immune-related GO terms and pathways. CeRNA network analysis showed that the key differentially expressed circRNAs and lncRNAs mainly regulate the key immune-related genes ITGB2, LCP2, PTPRC, SYK, and ZAP70 through competitive binding with miR-141-x, miR-29-y, and chi-miR-29b-3p, thereby affecting the natural killer cell-mediated cytotoxicity pathway, the T cell receptor signaling pathway, and other immune-related pathways. It should be noted that the expression of key circRNAs, lncRNAs, and key immune-related genes in goat SMGs decreased significantly with the growth of the goat. This is the first reporting of lncRNAs, circRNAs, and ceRNA network regulation in goat SMGs. Our study contributes to the knowledge of changes in the expression of non-coding RNAs during SMG development in goats and provides new insights into the relationship between non-coding RNAs and salivary gland immune function in goats.

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.

New Roles for EVs, miRNA and lncRNA in Bovine Embryo Implantation

The sine qua non of new life is fertilization. However, approximately 50% of fertilized eggs/blastocysts in cattle and up to 75% of those from human assisted reproductive procedures fail during the first 3 to 4 weeks of pregnancy, including peri-implantation periods. In these periods, blastocyst hatching and implantation to the maternal endometrium proceeds, during which physiological events such as epithelial-mesenchymal transition (EMT) and trophoblast cell fusion occur. Quite recently, extracellular vesicles (EVs) with micro RNAs (miRNAs) and long non-coding RNAs (lncRNAs) have been found to play a pivotal role for the establishment of the proper uterine environment required for peri-implantation processes to proceed. New findings of EVs, miRNA, and lncRNAs will be described and discussed to elucidate their connections with conceptus implantation to the maternal endometrium.

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.

Comprehensive Analysis of Long Noncoding RNA Modified by m6A Methylation in Oxidative and Glycolytic Skeletal Muscles

N⁶-methyladenosine (m⁶A) is the most common modification in eukaryotic RNAs. Accumulating evidence shows m⁶A methylation plays vital roles in various biological processes, including muscle and fat differentiation. However, there is a lack of research on lncRNAs’ m⁶A modification in regulating pig muscle-fiber-type conversion. In this study, we identified novel and differentially expressed lncRNAs in oxidative and glycolytic skeletal muscles through RNA-seq, and further reported the m⁶A-methylation patterns of lncRNAs via MeRIP-seq. We found that most lncRNAs have one m⁶A peak, and the m⁶A peaks were preferentially enriched in the last exon of the lncRNAs. Interestingly, we found that lncRNAs’ m⁶A levels were positively correlated with their expression homeostasis and levels. Furthermore, we performed conjoint analysis of MeRIP-seq and RNA-seq data and obtained 305 differentially expressed and differentially m⁶A-modified lncRNAs (dme-lncRNAs). Through QTL enrichment analysis of dme-lncRNAs and PPI analysis for their cis-genes, we finally identified seven key m⁶A-modified lncRNAs that may play a potential role in muscle-fiber-type conversion. Notably, inhibition of one of the key lncRNAs, MSTRG.14200.1, delayed satellite cell differentiation and stimulated fast-to-slow muscle-fiber conversion. Our study comprehensively analyzed m⁶A modifications on lncRNAs in oxidative and glycolytic skeletal muscles and provided new targets for the study of pig muscle-fiber-type conversion.

Characterization of lncRNA functioning in ovine conceptuses and endometria during the peri-implantation period

In ruminants, RNA-sequence analyses have revealed many characteristics of transcripts expressed in conceptuses (embryo and extraembryonic membrane) during peri-implantation periods; however, lncRNA profiles are yet characterized. In this study, we aimed to characterize the lncRNA expression profile in conceptuses during peri-implantation periods in sheep. We analyzed the RNA-sequence data of ovine conceptuses and endometria obtained from pregnant animals on days 15, 17, 19 and 21 (day 0 = day of estrus, n = 3 or 4/day). We predicted the protein coding ability of the assembled transcripts to identify the lncRNA candidates. This analysis identified 8808 lncRNAs, 3423 of which were novel lncRNAs. Gene ontology analysis revealed that lncRNA target genes were enriched for biological processes involved in the respiratory electron transport chain (RETC). qPCR analysis demonstrated that the expression levels on transcripts encoding RETC such as mitochondrially encoded cytochrome c oxidase II (MTCO2) and mitochondria DNA copy number in conceptuses were not increased on P21, although western blotting analysis and immunohistochemistry demonstrated that MTCO2 protein in conceptuses was increased on P21. NAD/NADH assay revealed that NADH level in conceptuses was increased on P21. These results indicate that lncRNAs could regulate the RETC through post-transcriptional levels in the conceptuses. Therefore, lncRNA is a potential new regulator in ovine conceptus development during peri-implantation periods.

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.

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.

The Complicated Effects of Extracellular Vesicles and Their Cargos on Embryo Implantation

As a rate-limiting step in pregnancy, embryo implantation is highly dependent on intercellular communication. Extracellular vesicles (EVs) are newly identified to be important in the course of intercellular communication. EVs have been isolated from a wide variety of biofluids and tissues, including plasma, liver, uterine, semen, embryo, etc. The present and future use of EVs not only as biomarkers, but also as targeting drug delivery system, is promisingly pave the way for advanced comprehension of implantation failure in reproductive diseases. However, as the precise mechanisms of EVs in embryo implantation has not been elucidated yet. Herein, we summarize the current knowledge on the diverse effects of EVs from various sources and their cargos such as microRNA, long non-coding RNA, protein, etc. on embryo implantation, and the potential mechanisms of EVs in reproductive diseases such as recurrent implantation failure, polycystic ovary syndrome and endometriosis. It is essential to note that many of the biologically plausible functions of EVs in embryo implantation discussed in present literatures still need further research in vivo.