Integrated analyses of miRNA-mRNA expression profiles of ovaries reveal the crucial interaction networks that regulate the prolificacy of goats in the follicular phase

Characterization and Comparative Analysis of Whole-Transcriptome Sequencing in High- and Low-Fecundity Chongming White Goat Ovaries during the Estrus Phase

Reproductive performance is one of the most important economic traits in the goat industry. Increasing the number of goats is an effective measure to improve production efficiency and reduce production costs. Ovaries are important reproductive organs in female mammals that directly affect the estrous cycle and reproductive abilities. Understanding the complex transcription network of non-coding RNAs (lncRNAs, circRNAs, and miRNAs) and messenger RNA (mRNA) could lead to significant insights into the ovarian regulation of the reproductive processes of animals. However, the whole-transcriptome analysis of the non-coding RNAs and mRNA of the ovaries in Chongming white goats between high-fecundity (HP) and low-fecundity (LP) groups is limited. In this study, a whole-transcriptome sequencing approach was used to identify lncRNA, circRNA, miRNA, and mRNA expression in the ovaries of Chongming white goats during the estrus phase using RNA-Seq technology. More than 20,000 messenger RNAs (mRNAs), 10,000 long non-coding RNAs (lncRNAs), 3500 circular RNAs (circRNAs), and 1000 micro RNAs (miRNAs) were identified. A total of 1024 differential transcripts (724 mRNAs, 112 lncRNAs, 178 circRNAs, and 10 miRNAs) existing between the HP and the LP groups were revealed through a bioinformatics analysis. They were enriched in the prolactin signaling pathway, the Jak-STAT signaling pathway, and the GnRH signaling pathway, as well as various metabolic pathways. Differentially expressed mRNAs (such as LYPD6, VEGFA, NOS3, TNXB, and EPHA2) and miRNAs (such as miR-10a-5p) play key roles in the regulation of goat ovaries during the estrus phase. The enrichment of pathways related to reproduction, such as the Hippo, Hedgehog, PI3K-AKT, and MAPK signaling pathways, suggests that they might be involved in the prolificacy of goat ovaries. Overall, we identified several gene modules associated with goat fecundity and provided a basis for a molecular mechanism in the ovaries of Chongming white goats.

Transcriptome Analysis Reveals Differentially Expressed circRNAs Associated with Fecundity in Small-Tail Han Sheep Thyroid with Different FecB Genotypes

Litter size is an economically important trait in sheep, and it is a complex trait controlled by multiple genes in multiple organs. Among them, the regulation of lamb number trait by the thyroid gland is a very important part. However, the molecular mechanisms of the thyroid gland in sheep reproduction remain unclear. Here, RNA-seq was used to detect transcriptome expression patterns in the thyroid gland between follicular phase (FP) and luteal phase (LP) in FecB BB (MM) and FecB ++ (ww) STH sheep, respectively, and to identify differentially expressed circRNAs (DECs) associated with reproduction. Bioinformatic analysis of the source genes of these DECs revealed that they can be enriched in multiple signaling pathways involved in the reproductive process of animals. We found that the source genes of these DECs, such as GNAQ, VEGFC, MAPK1, STAT1, and HSD17B7, may play important roles in the reproductive process of animals. To better understand the function of these DECs, we constructed circRNA-miRNA co-expression networks. Dual luciferase reporter assays suggested that a ceRNA regulatory mechanism between circ_0003259-oar-miR-133-TXLNA and circ_0012128-oar-miR-370-3p-FGFR1 may hold. All of these DEC expression profiles in the thyroid gland provide a novel resource for elucidating the regulatory mechanisms underlying STH sheep prolificacy.

Landscape genomics reveals adaptive divergence of indigenous sheep in different ecological environments of Xinjiang, China

Sheep are important livestock animals that have evolved under various ecological pressures. Xinjiang is a region with diverse and harsh environments that have shaped many local sheep breeds with unique characteristics and environmental adaptability. However, these breeds are losing ecological flexibility due to the promotion of intensive farming practices. Here we sequenced 14 local sheep breeds from Xinjiang and analyzed their genetic structure and gene flow with other sheep breeds from neighboring regions. The Tibetan Plateau was the geographic origin of Xinjiang native sheep evolution. We performed genome-environment association analysis and identified Bio9: Mean Temperature of Driest Quarter and Bio15: Precipitation Seasonality as the key environmental factors affecting Xinjiang local sheep and the key genes involved in their survival and adaptation. We classified Xinjiang native sheep breeds into six groups based on their differential genes by pairwise selective sweep analysis and Community Network Analysis. We analyzed transcriptome expression data of 832 sheep tissues and detected tissue-specific enrichment of six group-specific genes in different biological systems. Our results revealed the genetic basis of year-round estrus, drought tolerance, hypoxia resistance, and cold tolerance traits of Xinjiang sheep breeds. Moreover, we proposed conservation strategies for Xinjiang local sheep breeds and provided theoretical guidance for breeding new sheep breeds under global extreme environments.

Oogenesis in Women: From Molecular Regulatory Pathways and Maternal Age to Stem Cells

It is a well-known fact that the reproductive organs in women, especially oocytes, are exposed to numerous regulatory pathways and environmental stimuli. The maternal age is one cornerstone that influences the process of oocyte fertilization. More precisely, the longer a given oocyte is in the waiting-line to be ovulated from menarche to menopause, the longer the duration from oogenesis to fertilization, and therefore, the lower the chances of success to form a viable embryo. The age of menarche in girls ranges from 10 to 16 years, and the age of menopause in women ranges from approximately 45 to 55 years. Researchers are paying attention to the regulatory pathways that are impacting the oocyte at the very beginning during oogenesis in fetal life to discover genes and proteins that could be crucial for the oocyte’s lifespan. Due to the general trend in industrialized countries in the last three decades, women are giving birth to their first child in their thirties. Therefore, maternal age has become an important factor impacting oocytes developmental competence, since the higher a woman’s age, the higher the chances of miscarriage due to several causes, such as aneuploidy. Meiotic failures during oogenesis, such as, for instance, chromosome segregation failures or chromosomal non-disjunction, are influencing the latter-mentioned aging-related phenomenon too. These errors early in life of women can lead to sub- or infertility. It cannot be neglected that oogenesis is a precisely orchestrated process, during which the oogonia and primary oocytes are formed, and RNA synthesis takes place. These RNAs are crucial for oocyte growth and maturation. In this review, we intend to describe the relevance of regulatory pathways during the oogenesis in women. Furthermore, we focus on molecular pathways of oocyte developmental competence with regard to maternal effects during embryogenesis. On the background of transcriptional mechanisms that enable the transition from a silenced oocyte to a transcriptionally active embryo, we will briefly discuss the potential of induced pluripotent stem cells.

Temporal transcriptomic changes in microRNAs involved in the host immune response and metabolism during Neospora caninum infection

BACKGROUND

Neospora caninum infection is a major cause of abortion in cattle, which results in serious economic losses to the cattle industry. However, there are no effective drugs or vaccines for the control of N. caninum infections. There is increasing evidence that microRNAs (miRNAs) are involved in many physiological and pathological processes, and dysregulated expression of host miRNAs and the biological implications of this have been reported for infections by various protozoan parasites. However, to our knowledge, there is presently no published information on host miRNA expression during N. caninum infection.

METHODS

The expression profiles of miRNAs were investigated by RNA sequencing (RNA-seq) in caprine endometrial epithelial cells (EECs) infected with N. caninum at 24 h post infection (pi) and 48 hpi, and the functions of differentially expressed (DE) miRNAs were predicted by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. The transcriptome data were validated by using quantitative real-time polymerase chain reaction. One of the upregulated DEmiRNAs, namely chi-miR-146a, was selected to study the effect of DEmiRNAs on the propagation of N. caninum tachyzoites in caprine EECs.

RESULTS

RNA-seq showed 18 (17 up- and one downregulated) and 79 (54 up- and 25 downregulated) DEmiRNAs at 24 hpi and 48 hpi, respectively. Quantitative real-time polymerase chain reaction analysis of 13 randomly selected DEmiRNAs (10 up- and three downregulated miRNAs) confirmed the validity of the RNA-seq data. A total of 7835 messenger RNAs were predicted to be potential targets for 66 DEmiRNAs, and GO and KEGG enrichment analysis of these predicted targets revealed that DEmiRNAs altered by N. caninum infection may be involved in host immune responses (e.g. Fc gamma R-mediated phagocytosis, Toll-like receptor signaling pathway, tumor necrosis factor signaling pathway, transforming growth factor-β signaling pathway, mitogen-activated protein kinase signaling pathway) and metabolic pathways (e.g. lysine degradation, insulin signaling pathway, AMP-activated protein kinase signaling pathway, Rap1 signaling pathway, calcium signaling pathway). Upregulated chi-miR-146a was found to promote N. caninum propagation in caprine EECs.

CONCLUSIONS

This is, to our knowledge, the first report on the expression profiles of host miRNAs during infection with N. caninum, and shows that chi-miR-146a may promote N. caninum propagation in host cells. The novel findings of the present study should help to elucidate the interactions between host cells and N. caninum.

Identification of the genetic basis of sow pelvic organ prolapse

Introduction: Pelvic organ prolapse (POP) is one contributor to recent increases in sow mortality that have been observed in some populations and environments, leading to financial losses and welfare concerns. Methods: With inconsistent previous reports, the objective here was to investigate the role of genetics on susceptibility to POP, using data on 30,429 purebred sows, of which 14,186 were genotyped (25K), collected from 2012 to 2022 in two US multiplier farms with a high POP incidence of 7.1% among culled and dead sows and ranging from 2% to 4% of all sows present by parity. Given the low incidence of POP for parities 1 and >6, only data from parities 2 to 6 were retained for analyses. Genetic analyses were conducted both across parities, using cull data (culled for POP versus another reason), and by parity, using farrowing data. (culled for POP versus culled for another reason or not culled). Results and Discussion: Estimates of heritability from univariate logit models on the underlying scale were 0.35 ± 0.02 for the across-parity analysis and ranged from 0.41 ± 0.03 in parity 2 to 0.15 ± 0.07 in parity 6 for the by-parity analyses. Estimates of genetic correlations of POP between parities based on bivariate linear models indicated a similar genetic basis of POP across parities but less similar with increasing distance between parities. Genome wide association analyses revealed six 1 Mb windows that explained more than 1% of the genetic variance in the across-parity data. Most regions were confirmed in several by-parity analyses. Functional analyses of the identified genomic regions showed a potential role of several genes on chromosomes 1, 3, 7, 10, 12, and 14 in susceptibility to POP, including the Estrogen Receptor gene. Gene set enrichment analyses showed that genomic regions that explained more variation for POP were enriched for several terms from custom transcriptome and gene ontology libraries. Conclusion: The influence of genetics on susceptibility to POP in this population and environment was confirmed and several candidate genes and biological processes were identified that can be targeted to better understand and mitigate the incidence of POP.

Integrated Transcriptome Analysis Reveals the Crucial mRNAs and miRNAs Related to Fecundity in the Hypothalamus of Yunshang Black Goats during the Luteal Phase

A normal estrus cycle is essential for the breeding of goats, and the luteal phase accounts for most of the estrus cycle. The corpus luteum (CL) formed during the luteal phase is a transient endocrine gland that is crucial for the reproductive cycle and pregnancy maintenance, and is controlled by many regulatory factors. However, the molecular mechanism of the hypothalamus effect on the reproductive performance of different litter sizes during the luteal phase of goats has not been elucidated. In this study, RNA-sequencing was used to analyze the mRNA and miRNA expression profiles of the hypothalamic tissues with the high-fecundity goats during the luteal phase (LP-HF) and low-fecundity goats during the luteal phase (LP-LF). The RNA-seq results found that there were 1963 differentially expressed genes (DEGs) (890 up-regulated and 1073 down-regulated). The miRNA-seq identified 57 differentially expressed miRNAs (DEMs), including 11 up-regulated and 46 down-regulated, of which 199 DEGs were predicted to be potential target genes of DEMs. Meanwhile, the functional enrichment analysis identified several mRNA-miRNA pairs involved in the regulation of the hypothalamic activity, such as the common target gene MEA1 of novel-miR-972, novel-miR-125 and novel-miR-403, which can play a certain role as a related gene of the reproductive development in the hypothalamic-pituitary-gonadal (HPG) axis and its regulated network, by regulating the androgen secretion. While another target gene ADIPOR2 of the novel-miR-403, is distributed in the hypothalamus and affects the reproductive system through a central role on the HPG axis and a peripheral role in the gonadal tissue. An annotation analysis of the DE miRNA-mRNA pairs identified targets related to biological processes, such as anion binding (GO:0043168) and small molecule binding (GO: 0036094). Subsequently, the KEGG(Kyoto Encyclopedia of Genes and Genomes) pathways were performed to analyze the miRNA-mRNA pairs with negatively correlated miRNAs. We found that the GnRH signaling pathway (ko04912), the estrogen signaling pathway (ko04915), the Fc gamma R-mediated phagocytosis (ko04666), and the IL-17 signaling pathway (ko04657), etc., were directly and indirectly associated with the reproductive process. These targeting interactions may be closely related to the reproductive performance of goats. The results of this study provide a reference for further research on the molecular regulation mechanism for the high fertility in goats.

Why Is It So Difficult To Have Competent Oocytes from In vitro Cultured Preantral Follicles?

The developmental competence of oocytes is acquired gradually during follicular development, mainly through oocyte accumulation of RNA molecules and proteins that will be used during fertilization and early embryonic development. Several attempts to develop in vitro culture systems to support preantral follicle development up to maturation are reported in the literature, but oocyte competence has not yet been achieved in human and domestic animals. The difficulties to have fertilizable oocytes are related to thousands of mRNAs and proteins that need to be synthesized, long-term duration of follicular development, size of preovulatory follicles, composition of in vitro culture medium, and the need of multi-step culture systems. The development of a culture system that maintains bidirectional communication between the oocyte and granulosa cells and that meets the metabolic demands of each stage of follicle growth is the key to sustain an extended culture period. This review discusses the physiological and molecular mechanisms that determine acquisition of oocyte competence in vitro, like oocyte transcriptional activity, follicle and oocyte sizes, and length and regulation of follicular development in murine, human, and domestic animal species. The state of art of in vitro follicular development and the challenges to have complete follicular development in vitro are also highlighted.

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.

Differentially Expressed Circular RNA Profile Signatures Identified in Prolificacy Trait of Yunshang Black Goat Ovary at Estrus Cycle

CircRNAs acting as miRNA sponges play important roles in the growth process of animal individuals. The prolificacy trait of goats is involved in many pathways, however, the variation of circRNA expression profiles in the different phases of the estrus cycle at high and low fecundity groups is still unknown. Here, we analyzed the circRNA profiles of ovarian tissues among high and low fecundity groups in the follicular phase (HF vs LF), high and low fecundity groups in the luteal phase (HL vs LL), and high and low fecundity in the whole estrus cycle (HF vs HL and LF vs LL) using RNA-seq. A total of 283 (114 upregulated and 169 downregulated), 559 (299 upregulated and 260 downregulated), 449 (254 upregulated and 195 downregulated), and 314 (210 upregulated and 104 downregulated) differentially expressed (DE) circRNAs were screened in HF vs LF, HF vs HL, HL vs LL, and LF vs LL groups, respectively. Enrichment analysis suggested that the targeting of DE circRNAs was mainly enriched in oocyte meiosis, the GnRH signaling pathway, and estrogen signaling pathway. After integrating our previous study of miRNA-seq, there were 56 miRNAs that could target to 192 DE circRNAs, including the miR-133 family (including miR-133a-3p and miR-133b), miR-129-3p, and miR-21, which also had important influence on the prolificacy trait of goats. Then, 18 circRNAs with coding potential were obtained by four software predictions, and 9 circRNAs were validated by RT-qPCR. Together, circRNAs play a key role in the prolificacy trait and the transformation of the follicular phase to the luteal phase in the estrus cycle of goats.