RNA-seq Transcriptome Analysis in Ovarian Tissue of Pelibuey Breed to Explore the Regulation of Prolificacy

Thyroid transcriptomic profiling reveals the differential regulation of lncRNA and mRNA related to prolificacy in Small Tail Han sheep with FecB++ genotype

The thyroid gland is an important endocrine gland in animals, which mainly secretes thyroid hormones and acts on various organs of the body. Long-chain non-coding RNA (lncRNA) plays an important role in animal reproduction. However, there is still a lack of understanding of their expression patterns and potential roles in the thyroid of Small Tail Han (STH) sheep. In this study, RNA-seq was used to examine the transcriptome expression patterns of lncRNAs and mRNAs in the follicular phase (ww_FT) and luteal phase (ww_LT) in FecB++ genotype STH Sheep. A total of 17,217 lncRNAs and 39,112 mRNAs were identified including 96 differentially expressed lncRNAs (DELs) and 1054 differentially expressed mRNAs (DEGs). Functional analysis of genes with significant differences in expression level showed that these genes could be enriched in Ras signalling pathway, hedgehog (HH) signalling pathway, ATP-binding cassette (ABC) transporters and other signalling pathways related to animal reproduction. In addition, through correlation analysis for lncRNA-mRNA co-expression and network construction, we found that LNC_009115 and LNC_005796 trans target NIK-related kinase (NRK) and poly(A)-specific ribonuclease (PARN). LNC_007189 and LNC_002045 trans target progesterone-induced blocking factor 1 (PIBF1), LNC_009013 trans targets small mothers against decapentaplegic (SMAD1) are related to animal reproduction. These genes add new resources for elucidating the regulatory mechanisms of reproduction in sheep with different reproductive cycles of the FecB++ genotype STH sheep.

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.

CREB1 Is Involved in miR-134-5p-Mediated Endometrial Stromal Cell Proliferation, Apoptosis, and Autophagy

The successful establishment of endometrial receptivity is a key factor in ensuring the fertility of ewes and their economic benefits. Hu sheep have attracted attention due to their high fecundity and year-round estrus. In this study, we found that in the luteal phase, the uterine gland density, uterine coefficient, and number of uterine caruncles of high-fertility Hu sheep were higher than those of low-fertility Hu sheep. Thousands of differentially expressed genes were identified in the endometrium of Hu sheep with different fertility potential using RNA sequencing (RNA-Seq). Several genes involved in endometrial receptivity were screened using bioinformatics analysis. The qRT-PCR analysis further revealed the differential expression of cAMP reactive element binding protein-1 (CREB1) in the Hu sheep endometrium during the estrous cycle. Functionally, our results suggested that CREB1 significantly affected the expression level of endometrial receptivity marker genes, promoted cell proliferation by facilitating the transition from the G1 phase to the S phase, and inhibited cell apoptosis and autophagy. Moreover, we observed a negative linear correlation between miR-134-5p and CREB1 in the endometrium. In addition, CREB1 overexpression prevented the negative effect of miR-134-5p on endometrial stromal cell (ESC) growth. Taken together, these data indicated that CREB1 was regulated by miR-134-5p and may promote the establishment of uterine receptivity by regulating the function of ESCs. Moreover, this study provides new theoretical references for identifying candidate genes associated with fertility.

Genome wide identification and characterization of fertility associated novel CircRNAs as ceRNA reveal their regulatory roles in sheep fecundity

Reproductive traits play a vital role in determining the production efficiency of sheep. Maximizing the production is of paramount importance for breeders worldwide due to the growing population. Circular RNAs (circRNAs) act as miRNA sponges by absorbing miRNA activity through miRNA response elements (MREs) and participate in ceRNA regulatory networks (ceRNETs) to regulate mRNA expression. Despite of extensive research on role of circRNAs as miRNA sponges in various species, their specific regulatory roles and mechanism in sheep ovarian tissue are still not well understood. In this study, we performed whole genome sequencing of circRNAs, miRNA and mRNA employing bioinformatic techniques on ovine tissues of two contrasting sheep breeds "Small tail Han (X_LC) and Dolang sheep (D_LC)", which results into identification of 9,878 circRNAs with a total length of 23,522,667 nt and an average length of 2,381.32 nt. Among them, 44 differentially expressed circRNAs (DECs) were identified. Moreover, correlation between miRNA-mRNA and lncRNA-miRNA provided us with to prediction of miRNA binding sites on nine differentially expressed circRNAs and 165 differentially expressed mRNAs using miRanda. miRNA-mRNA and lncRNA-miRNA pairs with negative correlation were selected to determine the ceRNA score along with positively correlated pairs from lncRNA and mRNA network. Integration of ceRNA score and positively correlated pairs exhibit a significant ternary relationship among circRNAs-miRNA-mRNA demonestrated by ceRNA, comprising of 50 regulatory pairs sharring common nodes and predicted potential differentially expressed circRNAs-miRNAs-mRNAs regulatory axis. Based on functional enrichment analysis shortlisted key ceRNA regulatory pairs associated with reproduction including circRNA_3257-novel579_mature-EPHA3, circRNA_8396-novel130_mature-LOC101102473, circRNA_4140- novel34_mature > novel661_mature-KCNK9, and circRNA_8312-novel339_mature-LOC101110545. Furthermore, expression profiling, functional enrichments and qRT-PCR analysis of key target genes infer their implication in reproduction and metabolism. ceRNA target mRNAs evolutionary trajectories, expression profiling, functional enrichments, subcellular localizations following genomic organizations will provide new insights underlying molecular mechanisms of reproduction, and establish a solid foundation for future research. Graphical abstract summarizing the scheme of study.

ADAMTS Proteases: Importance in Animal Reproduction

Many reproductive physiological processes, such as folliculogenesis, ovulation, implantation, and fertilization, require the synthesis, remodeling, and degradation of the extracellular matrix (ECM). The ADAMTS (A Disintegrin and Metalloproteinase with Thrombospondin Motifs) family genes code for key metalloproteinases in the remodeling process of different ECM. Several genes of this family encode for proteins with important functions in reproductive processes; in particular, ADAMTS1, 4, 5 and 9 are genes that are differentially expressed in cell types and the physiological stages of reproductive tissues. ADAMTS enzymes degrade proteoglycans in the ECM of the follicles so that the oocytes can be released and regulate follicle development during folliculogenesis, favoring the action of essential growth factors, such as FGF-2, FGF-7 and GDF-9. The transcriptional regulation of ADAMTS1 and 9 in preovulatory follicles occurs because of the gonadotropin surge in preovulatory follicles, via the progesterone/progesterone receptor complex. In addition, in the case of ADAMTS1, pathways involving protein kinase A (PKA), extracellular signal regulated protein kinase (ERK1/2) and the epidermal growth factor receptor (EGFR) might contribute to ECM regulation. Different Omic studies indicate the importance of genes of the ADAMTS family from a reproductive aspect. ADAMTS genes could serve as biomarkers for genetic improvement and contribute to enhance fertility and animal reproduction; however, more research related to these genes, the synthesis of proteins encoded by these genes, and regulation in farm animals is needed.

Comprehensive Analysis of Differentially Expressed CircRNAs in the Ovaries of Low- and High-Fertility Sheep

CircRNAs are essential in regulating follicle growth and development and the female reproductive system at multiple levels. However, the molecular mechanism by which circRNAs regulate reproduction in sheep is unclear and requires further exploration. In this study, RNA sequencing was performed to reveal the circRNA expression profiles in the ovaries of Cele black sheep and Hetian sheep during estrus. Analysis of the number of circRNAs in their host genes revealed that 5031 genes could produce 20,835 circRNAs. Among the differentially expressed circRNAs (DEcircRNA), 75 were upregulated, and 105 were downregulated. Functional enrichment analysis showed that the host genes of DEcircRNA were involved in several pathways, including the MAPK and Hippo signaling pathway. In addition, we constructed a subnetwork of competitive endogenous RNA (ceRNA) containing 4 mRNAs, 4 microRNAs (miRNAs), and 10 circRNAs, potentially related to follicle development. Functional circRNAs (e.g., novel_circ_0003851, novel_circ_0015526, novel_circ_0008117) were found to act as ceRNAs for follicle growth and development-related mRNAs (CUEDC1, KPNB1, ZFPM2) by sponging functional miRNAs (miR-29a, miR-29b, miR-17-5p). Finally, through an RNA pull-down assay, oar-miR-125b was selected and confirmed as the target miRNA of novel-circ-0041512. We analyzed the overall expression of circRNAs in sheep ovaries. Further, we explored the potential mechanisms underlying the circRNA functions, providing a theoretical basis for the genetic progress of reproductive traits in sheep.

Temperature Incubation Influences Gonadal Gene Expression during Leopard Gecko Development

During development, sexual differentiation results in physiological, anatomical and metabolic differences that implicate not only the gonads but also other body structures. Sex in Leopard geckos is determined by egg incubation temperature. Based on the premise that the developmental decision of gender does not depend on a single gene, we performed an analysis on E. macularius to gain insights into the genes that may be involved in gonads' sexual differentiation during the thermosensitive period. All the genes were identified as differentially expressed at stage 30 during the labile phase of sex differentiation. In this way, the expression of genes known to be involved in gonadal sexual differentiation, such as WNT4, SOX9, DMRT1, Erα, Erβ, GnRH, P450 aromatase, PRLand PRL-R, was investigated. Other genes putatively involved in sex differentiation were sought by differential display. Our findings indicate that embryo exposure to a sex-determining temperature induces differential expression of several genes that are involved not only in gonadal differentiation, but also in several biological pathways (ALDOC, FREM1, BBIP1, CA5A, NADH5, L1 non-LTR retrotransposons, PKM). Our data perfectly fit within the new studies conducted in developmental biology, which indicate that in the developing embryo, in addition to gonadal differentiation, sex-specific tissue and metabolic polarization take place in all organisms.

Unrevealing functional candidate genes for bovine fertility through RNA sequencing meta-analysis and regulatory elements networks of co-expressed genes and lncRNAs

The high genetic heterogeneity and environmental effects of subfertility in livestock species make the elucidation of the genetic mechanisms associated with reproductive efficiency a difficult task. Network and co-expression network meta-analyses were applied alongside genetic variant calling and long non-coding RNA (lncRNA) characterization to identify functionally relevant target genes and regulatory subnetworks associated with fertility in dairy cattle. In total, 505 lncRNAs (441 previously annotated in the bovine reference genome ARS-UCD 1.2 and 64 novel lncRNAs) were identified. Seven differentially expressed genes between high-fertile (HF) and sub-fertile (SF) Holstein cows were identified in the network meta-analysis (CA5A, ENSBTAG00000051149, ENSBTAG00000003272, DEFB7, DIO2, TRPV3, and COL4A4). Additionally, seven functional candidate differentially co-expressed (DcoExp) modules with a differential regulatory pattern (|z-score|>2) were identified between HF and SF cows. The functional candidate genes and DcoExp modules identified were associated with fertility relevant processes such as the regulation of embryonic implantation and proliferation, interaction and molecule transfer between the fetus and the cow, and the immune system. These results help to better understand the genetic mechanisms associated with reproductive efficiency in dairy cattle through the identification of potential biomarkers and genetic variants associated with differentially expressed regulatory gene and lncRNAs regulatory element networks.

lncRNA-miRNA-mRNA ceRNA Network Involved in Sheep Prolificacy: An Integrated Approach

Understanding the molecular pattern of fertility is considered as an important step in breeding of different species, and despite the high importance of the fertility, little success has been achieved in dissecting the interactome basis of sheep fertility. However, the complex mechanisms associated with prolificacy in sheep have not been fully understood. Therefore, this study aimed to use competitive endogenous RNA (ceRNA) networks to evaluate this trait to better understand the molecular mechanisms responsible for fertility. A competitive endogenous RNA (ceRNA) network of the corpus luteum was constructed between Romanov and Baluchi sheep breeds with either good or poor genetic merit for prolificacy using whole-transcriptome analysis. First, the main list of lncRNAs, miRNAs, and mRNA related to the corpus luteum that alter with the breed were extracted, then miRNA−mRNA and lncRNA−mRNA interactions were predicted, and the ceRNA network was constructed by integrating these interactions with the other gene regulatory networks and the protein−protein interaction (PPI). A total of 264 mRNAs, 14 lncRNAs, and 34 miRNAs were identified by combining the GO and KEGG enrichment analyses. In total, 44, 7, 7, and 6 mRNAs, lncRNAs, miRNAs, and crucial modules, respectively, were disclosed through clustering for the corpus luteum ceRNA network. All these RNAs involved in biological processes, namely proteolysis, actin cytoskeleton organization, immune system process, cell adhesion, cell differentiation, and lipid metabolic process, have an overexpression pattern (Padj < 0.01). This study increases our understanding of the contribution of different breed transcriptomes to phenotypic fertility differences and constructed a ceRNA network in sheep (Ovis aries) to provide insights into further research on the molecular mechanism and identify new biomarkers for genetic improvement.

MiR-124 and Small Molecules Synergistically Regulate the Generation of Neuronal Cells from Rat Cortical Reactive Astrocytes

Irreversible neuron loss caused by central nervous system injuries usually leads to persistent neurological dysfunction. Reactive astrocytes, because of their high proliferative capacity, proximity to neuronal lineage, and significant involvement in glial scarring, are ideal starting cells for neuronal regeneration. Having previously identified several small molecules as important regulators of astrocyte-to-neuron reprogramming, we established herein that miR-124, ruxolitinib, SB203580, and forskolin could co-regulate rat cortical reactive astrocyte-to-neuron conversion. The induced cells had reduced astroglial properties, displayed typical neuronal morphologies, and expressed neuronal markers, reflecting 25.9% of cholinergic neurons and 22.3% of glutamatergic neurons. Gene analysis revealed that induced neuron gene expression patterns were more similar to that of primary neurons than of initial reactive astrocytes. On the molecular level, miR-124-driven neuronal differentiation of reactive astrocytes was via targeting of the SOX9-NFIA-HES1 axis to inhibit HES1 expression. In conclusion, we present a novel approach to inducing endogenous rat cortical reactive astrocytes into neurons through co-regulation involving miR-124 and three small molecules. Thus, our research has potential implications for inhibiting glial scar formation and promoting neuronal regeneration after central nervous system injury or disease.

Genome-wide transcriptome profiling uncovers differential miRNAs and lncRNAs in ovaries of Hu sheep at different developmental stages

Ovary development is an important determinant of the procreative capacity of female animals. Here, we performed genome-wide sequencing of long non-coding RNAs (lncRNAs) and mRNAs on ovaries of 1, 3 and 8 months old Hu sheep to assess their expression profiles and roles in ovarian development. We identified 37,309 lncRNAs, 45,404 messenger RNAs (mRNAs) and 330 novel micro RNAs (miRNAs) from the transcriptomic analysis. Six thousand, seven hundred and sixteen (6716) mRNAs and 1972 lncRNAs were significantly and differentially expressed in ovaries of 1 month and 3 months old Hu sheep (H1 vs H3). These mRNAs and target genes of lncRNAs were primarily enriched in the TGF-β and PI3K-Akt signalling pathways which are closely associated with ovarian follicular development and steroid hormone biosynthesis regulation. We identified MSTRG.162061.1, MSTRG.222844.7, MSTRG.335777.1, MSTRG.334059.16, MSTRG.188947.6 and MSTRG.24344.3 as vital genes in ovary development by regulating CTNNB1, CCNA2, CDK2, CDC20, CDK1 and EGFR expressions. A total of 2903 mRNAs and 636 lncRNAs were differentially expressed in 3 and 8 months old ovaries of Hu sheep (H3 vs H8); and were predominantly enriched in PI3K-Akt, progesterone-mediated oocyte maturation, estrogen metabolism, ovulation from the ovarian follicle and oogenesis pathways. These lncRNAs were also found to regulate FGF7, PRLR, PTK2, AMH and INHBA expressions during follicular development. Our result indicates the identified genes participate in the development of the final stages of follicles and ovary development in Hu sheep.

Genome-wide association study and pathway analysis identify NTRK2 as a novel candidate gene for litter size in sheep

Litter size is one of the most important economic traits in sheep. Identification of gene variants that are associated with the prolificacy rate is an important step in breeding program success and profitability of the farm. So, to identify genetic mechanisms underlying the variation in litter size in Iranian Baluchi sheep, a two-step genome-wide association study (GWAS) was performed. GWAS was conducted using genotype data from 91 Baluchi sheep. Estimated breeding values (EBVs) for litter size calculated for 3848 ewes and then used as the response variable. Besides, a pathway analysis using GO and KEGG databases were applied as a complementary approach. A total of three single nucleotide polymorphisms (SNPs) associated with litter size were identified, one each on OAR2, OAR10, and OAR25. The SNP on OAR2 is located within a novel putative candidate gene, Neurotrophic receptor tyrosine kinase 2. This gene product works as a receptor which is essential for follicular assembly, early follicular growth, and oocyte survival. The SNP on OAR25 is located within RAB4A which is involved in blood vessel formation and proliferation through angiogenesis. The SNP on OAR10 was not associated with any gene in the 1Mb span. Moreover, gene-set analysis using the KEGG database identified several pathways, such as Ovarian steroidogenesis, Steroid hormone biosynthesis, Calcium signaling pathway, and Chemokine signaling. Also, pathway analysis using the GO database revealed several functional terms, such as cellular carbohydrate metabolic, biological adhesion, cell adhesion, cell junction, and cell-cell adherens junction, among others. This is the first study that reports the NTRK2 gene affecting litter size in sheep and our study of this gene functions showed that this gene could be a good candidate for further analysis.

Ovarian transcriptomic analysis of black Muscovy duck at the early, peak and late egg-laying stages

Ovarian development is a complex process involving many genes and pathways. A well-developed ovary is essential for poultry to keep high egg production and egg fertility. In order to better understand the mechanism of egg production performance, a comparative transcriptomic analysis was performed on ovaries of black Muscovy ducks at the early (BE), peak (BP) and late laying (BL) stages. 1683 DEGs were identified from BL-vs-BE, BL-vs-BP and BP-vs-BE, and the up-regulated genes were 41, 835, 260, the down-regulated genes were 60, 255, 730, respectively. Besides, there were 32, 20 and 424 DEGs co-expressed in the two comparison groups, and 11 DEGs were co-expressed in the three comparison groups. HOXA10, HtrA3, StAR, ZP2 and TAT were found to be involved in the regulation of ovarian development were significantly differentially expressed at different laying stages, which helped to regulate ovarian maturation and egg production. Moreover, we discovered several important functional pathways, such as steroid hormone biosynthesis and ovarian steroidogenesis, that appear to be much more active in the BP ovary compared to those of the BE and BL. Furthermore, 17 coding and 244 non-coding new transcripts were detected in the three comparison groups, the gene structures were optimized and the gene annotation informations were improved. These findings will provide a solid foundation on ovarian development in black Muscovy ducks and other poultry animals at different laying stages, and help to understand the complex molecular and cellular mechanisms of ovary.

Comparative transcriptome of reproductive axis in Chinese indigenous sheep with different FecB genotypes and prolificacies

In the present study, there were 128, 263, and 139 differentially expressed genes (DEGs) from the ovary, pituitary, and hypothalamus, respectively, with the FecB B+/FecB ++ genotype. Results indicated that there was a synchronized, and differential increase in mRNA transcript abundances for genes encoding for ovarian steroidogenesis-signaling proteins in the ovary, TGF-β-signaling in hypothalamus, cAMP-signaling-transduction, and dopaminergic-synapse pathways in the pituitary. The values for these variables were associated with mean prolificacy in indigenous sheep and results provided indications for candidate regulator proteins that coordinate processes related to prolificacy in ewes. Furthermore, the results when there was evaluation of the FecB-genotype-specific co-expression modules implied that CYP17 was the hub gene connecting the pathways and for modulation of folliculogenesis and ovulation in the ewe with the FecB B+ genotype.

Gene Expression Profiling of Corpus luteum Reveals Important Insights about Early Pregnancy in Domestic Sheep

The majority of pregnancy loss in ruminants occurs during the preimplantation stage, which is thus the most critical period determining reproductive success. Here, we performed a comparative transcriptome study by sequencing total mRNA from corpus luteum (CL) collected during the preimplantation stage of pregnancy in Finnsheep, Texel and F1 crosses. A total of 21,287 genes were expressed in our data. Highly expressed autosomal genes in the CL were associated with biological processes such as progesterone formation (STAR, CYP11A1, and HSD3B1) and embryo implantation (e.g., TIMP1, TIMP2 and TCTP). Among the list of differentially expressed genes, sialic acid-binding immunoglobulin (Ig)-like lectins (SIGLEC3, SIGLEC14, SIGLEC8), ribosomal proteins (RPL17, RPL34, RPS3A, MRPS33) and chemokines (CCL5, CCL24, CXCL13, CXCL9) were upregulated in Finnsheep, while four multidrug resistance-associated proteins (MRPs) were upregulated in Texel ewes. A total of 17 known genes and two uncharacterized non-coding RNAs (ncRNAs) were differentially expressed in breed-wise comparisons owing to the flushing diet effect. The significantly upregulated TXNL1 gene indicated potential for embryonic diapause in Finnsheep and F1. Moreover, we report, for the first time in any species, several genes that are active in the CL during early pregnancy (including TXNL1, SIGLEC14, SIGLEC8, MRP4, and CA5A).