Key mRNAs and lncRNAs of pituitary that affect the reproduction of FecB + + small tail han sheep.

BACKGROUND: The pituitary directly regulates the reproductive process through follicle-stimulating hormone (FSH) and luteinizing hormone (LH). Transcriptomic research on the pituitaries of ewes with different FecB (fecundity Booroola) genotypes has shown that some key genes and lncRNAs play an important role in pituitary function and sheep fecundity. Our previous study found that ewes with FecB + + genotypes (without FecB mutation) still had individuals with more than one offspring per birth. It is hoped to analyze this phenomenon from the perspective of the pituitary transcriptome. RESULTS: The 12 Small Tail Han Sheep were equally divided into polytocous sheep in the follicular phase (PF), polytocous sheep in the luteal phase (PL), monotocous sheep in the follicular phase (MF), and monotocous sheep in the luteal phase (ML). Pituitary tissues were collected after estrus synchronous treatment for transcriptomic analysis. A total of 384 differentially expressed genes (DEGs) (182 in PF vs. MF and 202 in PL vs. ML) and 844 differentially expressed lncRNAs (DELs) (427 in PF vs. MF and 417 in PL vs. ML) were obtained from the polytocous-monotocous comparison groups in the two phases. Functional enrichment analysis showed that the DEGs in the two phases were enriched in signaling pathways known to play an important role in sheep fecundity, such as calcium ion binding and cAMP signaling pathways. A total of 1322 target relationship pairs (551 pairs in PF vs. MF and 771 pairs in PL vs. ML) were obtained for the target genes prediction of DELs, of which 29 DEL-DEG target relationship pairs (nine pairs in PF vs. MF and twenty pairs in PL vs. ML). In addition, the competing endogenous RNA (ceRNA) networks were constructed to explore the regulatory relationships of DEGs, and some important regulatory relationship pairs were obtained. CONCLUSION: According to the analysis results, we hypothesized that the pituitary first receives steroid hormone signals from the ovary and uterus and that VAV3 (Vav Guanine Nucleotide Exchange Factor 3), GABRG1 (Gamma-Aminobutyric Acid A Receptor, Gamma 1), and FNDC1 (Fibronectin Type III Domain Containing 1) played an important role in this process. Subsequently, the reproductive process was regulated by gonadotropins, and IGFBP1 (Insulin-like Growth Factor Binding Protein 1) was directly involved in this process, ultimately affecting litter size. In addition, TGIF1 (Transforming Growth Factor-Beta-Induced Factor 1) and TMEFF2 (Transmembrane Protein With EGF Like And Two Follistatin Like Domains 2) compensated for the effect of the FecB mutation and function by acting on TGF-ß/SMAD signaling pathway, an important pathway for sheep reproduction. These results provided a reference for understanding the mechanism of multiple births in Small Tail Han Sheep without FecB mutation.

Detection of polymorphisms in six genes and their association analysis with litter size in sheep.

Litter size in sheep is a complex trait controlled by micro-effective polygenes. APAF1, CLSTN2, CTH, PLCB1, PLCB4, and CHST11 are all involved in mammalian reproduction. However, the effects of these genes on litter size in sheep are still unclear. Therefore, in this study, we used Sequenom MassARRAY® SNP assay technology to type the single nucleotide polymorphisms (SNPs) loci of six genes in five sheep breeds. The results showed that most sheep breeds contain three genotypes at each locus. Then, we conducted population genetic analysis on the SNPs of six genes and found that the polymorphic information content in all sheep breeds ranged from 0 to 0.37, and most sheep breeds were in Hardy-Weinberg equilibrium (p > 0.05). In addition, association analysis in Small Tail Han sheep indicated that the rs399534524 locus in CLSTN2 was highly associated with first parity litter size, and litter size in ewes with CT genotype was higher than that in ewes with CC genotype or TT genotype. Furthermore, the rs407142552 locus in CTH was highly associated with second parity litter size in Small Tail Han sheep, and litter size in ewes with CT genotype was higher than that in ewes with TT genotype. Finally, we predicted the CTH and CLSTN2 protein interaction network and found that HTR1E, NOM1, CCDC174 and ALPK3 interact with CLSTN2 and have been reported as candidate genes related to litter size in sheep. These results suggest that they may be useful genetic markers for increasing litter size in sheep.

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.

Thyroid transcriptomic profiling reveals the differential regulation of lncRNA and mRNA related to prolificacy in small tail han sheep with FecB BB genotype.

BACKGROUND: The thyroid gland is an important endocrine gland in animals that secretes thyroid hormones and acts on various organs throughout the body. lncRNAs are long non-coding RNAs that play an important role in animal reproduction; however, there is a lack of understanding of their expression patterns and potential roles in the thyroid gland of the Small Tail Han (STH) sheep. In this study, we used RNA-Seq technology to examine the transcriptome expression pattern of the thyroid from the luteal phase (LP) and follicular phase (FP) of FecB BB (MM) STH sheep. RESULTS: We identified a total of 122 and 1287 differential expression lncRNAs (DELs) and differential expression mRNAs (DEGs), respectively, which were significantly differentially expressed. These DELs target genes and DEGs can be enriched in several signalling pathways related to the animal reproduction process. CONCLUSIONS: The expression profiles of DELs and DEGs in thyroid glands provide a more comprehensive resource for elucidating the reproductive regulatory mechanisms of STH sheep.

Effect of melatonin on ATG2B-mediated autophagy regulation in sheep granulosa cells with different FecB genotypes.

Melatonin (MLT) protects cells by reducing reactive oxygen species (ROS) levels, which are key for inducing cellular autophagy. The aim of this study was to investigate the molecular mechanisms underlying MLT regulation of autophagy in granulosa cells (GCs) with BMPR-1B homozygous (FecB BB) and wild type (FecB ++) mutations. GCs collected from small-tailed Han sheep with different FecB genotypes were typed using a TaqMan probe assay, and autophagy levels were found to be significantly higher in GCs with FecB BB than the levels in those with FecB ++. Autophagy-related 2 homolog B (ATG2B) was associated with cell autophagy and was highly expressed in GCs with the FecB BB genotype in small-tailed Han sheep. Overexpression of ATG2B in the GCs of sheep with both FecB genotypes promoted GC autophagy, and the contrary was observed after the inhibition of ATG2B expression. Subsequently, treatment of GCs with different genotypes of FecB and MLT revealed a significant decrease in cellular autophagy and an increase in ATG2B expression. Addition of MLT to GCs with inhibited ATG2B expression revealed that MLT could protect GCs by decreasing ROS levels, especially in GCs with FecB ++ genotype. In conclusion, this study determined that autophagy levels were significantly higher in sheep GCs with FecB BB genotype than the levels in those with FecB ++ genotype, which may have contributed to the difference in lambing numbers between the two FecB genotypes. Autophagy was regulated by ATG2B and was able to protect GCs by reducing the high levels of ROS produced following inhibition of ATG2B through the addition of MLT in vitro.

ALAS1 associated with goat kidding number trait was regulated by the transcription factor ASCL2 to affect granulosa cell proliferation.

ALAS1 is a member of the α-oxoamine synthase family, which is the first rate-limiting enzyme for heme synthesis and is important for maintaining intracellular heme levels. In the ovary, ALAS1 is associated with the regulation of ovulation-related mitochondrial P450 cytochromes, steroid metabolism, and steroid hormone production. However, there are few studies on the relationship between ALAS1 and reproductive traits in goats. In this study, a mutation located in the promoter region of ALAS1 (g.48791372C>A) was found to be significantly (p < 0.05) associated with the kidding number of Yunshang black goats. Specifically, the mean kidding number in the first three litters and the kidding numbers of all three litters were significantly (p < 0.05) higher in individuals with the CA genotype or AA genotype than in those with the CC genotype. To further investigate the regulatory mechanism of ALAS1, the expression of ALAS1 in goat ovarian tissues with different genotypes was verified by real-time quantitative PCR. The results showed that the expression of ALAS1 was significantly higher in the ovaries of individuals with AA genotype than those with AC and CC genotypes (p < 0.01), and the expression trend of transcription factor ASCL2 was consistent with ALAS1. Additionally, the ALAS1 g.48791372C>A mutation created a new binding site for the transcription factor ASCL2. The luciferase activity assay indicated that the mutation increased the promoter activity of ALAS1. Overexpression of the transcription factor ASCL2 induced increased expression of ALAS1 in goat granulosa cells (p < 0.05). The opposite trend was shown for the inhibition of ASCL2 expression. The results of real-time quantitative PCR, EdU and Cell Counting Kit-8 assays indicated that the transcription factor ASCL2 increased the proliferation of goat granulosa cells by mediating the expression of ALAS1. In conclusion, the transcription factor ASCL2 positively regulated the transcriptional activity and expression levels of ALAS1, altering granulosa cell proliferation and the kidding number in goats.

Identification of differentially expressed genes associated with precocious puberty by suppression subtractive hybridization in goat pituitary tissues.

The aim of this study was to identify genes related to precocious puberty expressed in the pituitary of goats at different growth stages by suppression subtractive hybridization (SSH). The pituitary glands from Jining Gray (JG) goats (early puberty) and Liaoning Cashmere (LC) goats (late puberty) at 30, 90, and 180 days were used in this study. To identify differentially expressed genes (DEGs) in the pituitary glands, mRNA was extracted from these tissues, and SSH libraries were constructed and divided into the following groups: juvenile group (30-JG vs. 30-LC, API), puberty group (90-JG vs. 180-LC, BPI), and control group (90-JG vs. 90-LC, EPI). A total of 60, 49, and 58 DEGs were annotated by 222 Gene Ontology (GO) terms and 75 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Most of the DEGs were significantly enriched in GO terms related to 'structural constituent of ribosome', 'translation' and 'GTP binding', and numerous DEGs were also significantly enriched in KEGG terms related to the Jak-STAT signaling and oocyte meiosis pathways. Candidate genes associated with precocious puberty and sexual development were screened from the SSH libraries. These genes were analyzed to determine if they were expressed in the pituitary tissues of the goats at different growth stages and to identify genes that may influence the hypothalamic-pituitary-gonadal (HPG) axis. In this study, we found precocious puberty-related genes (such as PRLP0, EIF5A, and YWHAH) that may be interesting from an evolutionary perspective and that could be investigated for use in future goat breeding programs. Our results provide a valuable dataset that will facilitate further research into the reproductive biology of goats.

Photoperiod Induces the Epigenetic Change of the GNAQ Gene in OVX+E2 Ewes.

GNAQ, a member of the alpha subunit encoding the q-like G protein, is a critical gene in cell signaling, and multiple studies have shown that upregulation of GNAQ gene expression ultimately inhibits the proliferation of gonadotropin-releasing hormone (GnRH) neurons and GnRH secretion, and ultimately affects mammalian reproduction. Photoperiod is a key inducer which plays an important role in gene expression regulation by affecting epigenetic modification. However, fewer studies have confirmed how photoperiod induces epigenetic modifications of the GNAQ gene. In this study, we examined the expression and epigenetic changes of GNAQ in the hypothalamus in ovariectomized and estradiol-treated (OVX+E2) sheep under three photoperiod treatments (short photoperiod treatment for 42 days, SP42; long photoperiod treatment for 42 days, LP42; 42 days of short photoperiod followed by 42 days of long photoperiod, SP-LP42). The results showed that the expression of GNAQ was significantly higher in SP-LP42 than in SP42 and LP42 (p < 0.05). Whole genome methylation sequencing (WGBS) results showed that there are multiple differentially methylated regions (DMRs) and loci between different groups of GNAQ. Among them, the DNA methylation level of DMRs at the CpG1 locus in SP42 was significantly higher than that of SP-LP42 (p < 0.01). Subsequently, we confirmed that the core promoter region of the GNAQ gene was located with 1100 to 1500 bp upstream, and the DNA methylation level of all eight CpG sites in SP42 was significantly higher than those in LP42 (p < 0.01), and significantly higher than those in SP-LP42 (p < 0.01), except site 2 and site 4 in the first sequencing fragment (p < 0.05) in the core promoter region. The expression of acetylated GNAQ histone H3 was significantly higher than that of the control group under three different photoperiods (p < 0.01); the acetylation level of sheep hypothalamic GNAQ genomic protein H3 was significantly lower under SP42 than under SP-LP42 (p < 0.05). This suggests that acetylated histone H3 binds to the core promoter region of the GNAQ gene, implying that GNAQ is epigenetically regulated by photoperiod through histone acetylation. In summary, the results suggest that photoperiod can induce DNA methylation in the core promoter region and histone acetylation in the promoter region of the GNAQ gene, and hypothesize that the two may be key factors in regulating the differential expression of GNAQ under different photoperiods, thus regulating the hypothalamus-pituitary-gonadal axis (HPGA) through the seasonal estrus in sheep. The results of this study will provide some new information to understand the function of epigenetic modifications in reproduction in sheep.

The mechanism of circadian clock and its influence on animal circannual rhythm.

The circadian clock exists in almost all life forms, and is an internal activity generated by organisms adapting to the daily periodic changes of the external environment. The circadian clock is regulated by the transcription-translation-negative feedback loop in the body, which can regulate the activities of tissues and organs. Its normal maintenance is important for the health, growth, and reproduction of organisms. In contrast, due to the season changes of the environment, organisms have also formed annual cycle physiological changes in their bodies, such as seasonal estrus, etc. The annual rhythm of living things is mainly affected by environmental factors such as photoperiod, and is related to gene expression, hormone content, morphological changes of cell and tissues in vivo. Melatonin is an important signal to recognize the changes of photoperiod, and the circadian clock plays an important role in the pituitary to interpret the signal of melatonin and regulate the changes of downstream signals, which plays an important guiding role in the recognition of annual changes in the environment and the generation of the body's annual rhythm. In this review, we summarize the progress of research on the mechanism of action of circadian clocks in influencing annual rhythms, by introducing the mechanisms of circadian and annual rhythms generation in insects and mammals, and in the context of annual rhythms in birds, with the aim of providing a broader range of ideas for future research on the mechanism of annual rhythms influence.

Progress on the effect of FecB mutation on BMPR1B activity and BMP/SMAD pathway in sheep.

BMPR1B is the first major gene of litter size identified in sheep. However, the molecular mechanism of the FecB mutation that increases the ovulation rate in sheep is still unclear. In recent years, it has been demonstrated that BMPR1B activity is regulated by the small molecule repressor protein FKBP1A, which acts as a key activity switch of the BMPR1B in the BMP/SMAD pathway. The FecB mutation is located close to the binding site of FKBP1A and BMPR1B. In this review, we summarize the structure of BMPR1B and FKBP1A proteins, and clarify the spatial interactive domains of the two proteins with respect to the location of the FecB mutation. Then the relationship between the FecB mutation and the degree of affinity of the two proteins are predicted. Finally, the hypothesis that FecB mutation causes change of activity in BMP/SMAD pathway by affecting the intensity of the interactions between BMPR1B and FKBP1A is proposed. This hypothesis provides a new clue to investigate the molecular mechanism of FecB mutation affecting ovulation rate and litter size in sheep.

Transcriptome study digs out BMP2 involved in adipogenesis in sheep tails.

BACKGROUND: Hu sheep and Tibetan sheep in China are characterized by fat tails and thin tails, respectively. Several transcriptomes have been conducted in different sheep breeds to identify the differentially expressed genes (DEGs) underlying this trait. However, these studies identified different DEGs in different sheep breeds. RESULTS: Hence, RNA sequencing was performed on Hu sheep and Tibetan sheep. We obtained a total of 45.57 and 43.82 million sequencing reads, respectively. Two libraries mapped reads from 36.93 and 38.55 million reads after alignment to the reference sequences. 2108 DEGs were identified, including 1247 downregulated and 861 upregulated DEGs. GO and KEGG analyses of all DEGs demonstrated that pathways were enriched in the regulation of lipolysis in adipocytes and terms related to the chemokine signalling pathway, lysosomes, and glycosaminoglycan degradation. Eight genes were selected for validation by RT-qPCR. In addition, the transfection of BMP2 overexpression into preadipocytes resulted in increased PPAR-γ expression and expression. BMP2 potentially induces adipogenesis through LOX in preadipocytes. The number of lipid drops in BMP2 overexpression detected by oil red O staining was also greater than that in the negative control. CONCLUSION: In summary, these results showed that significant genes (BMP2, HOXA11, PPP1CC and LPIN1) are involved in the regulation of adipogenesis metabolism and suggested novel insights into metabolic molecules in sheep fat tails.

Association between single-nucleotide polymorphism in PIK3CD gene and litter size in Small Tail Han sheep.

The p110δ isoform of the PI3K catalytic subunit (encoded by the PIK3CD gene) is a key component of the PI3K pathway for follicle growth in mammalian ovarian granulosa cells. Nevertheless, little is known about the association of its polymorphisms with ovine litter size. In this study, the distribution of different genotypes of two SNPs in the PIK3CD gene was calculated in more than forty sheep breeds, and the associations between SNPs and litter size in Small Tail Han (STH) sheep were also analyzed. Besides, the mRNA expression of the PIK3CD gene was also detected in some reproduction-related tissues. The results showed that the "A" allele frequency was higher in rs412889931 (g.41926327G > A) in a typical polytocous sheep breed (p < 0.01). The association's analysis showed rs412889931 was correlated with ovine fecundity as assessed by three parity litter sizes (p < 0.05). Finally, we found the expression of PIK3CD in the ovary had significant differences in different fecundity sheep breeds, indicating that SNP may regulate the litter size by influencing the PIK3CD gene expression. The present results demonstrated that rs412889931 could be used in the marker-assisted selection of the litter size in sheep breeding.

Genome-wide identification and functional prediction of circular RNAs in response to heat stress in Chinese Holstein cows.

Heat stress (HS) leads to substantial economic loss of dairy industry each year. The negative effect of HS in dairy cows is becoming one of the more urgent issue due to accelerating side-effects of global warming. Various genes are involved in HS response but the information about the role of noncoding RNAs, especially circular RNAs (circRNAs) is largely unknown. In our study, we aimed to investigate the different expression profile of circRNAs between HS and Non-heat-stressed condition (NC) of Chinese Holstein cow's mammary gland. CircRNAs were identified using RNA sequencing and bioinformatics analysis. In total, 37405 circRNAs were detected and 95 were differentially expressed (DE), including 15 downregulated and 80 upregulated circRNAs in HS group compared to NC. Eight circRNAs were randomly selected to verify the RNA sequencing result. Further, Sanger sequencing validated the backsplicing site of the eight circRNAs. Moreover, results obtained from the Quantitative real time PCR (qRT-PCR) showed consistent expression trend with that of RNA sequencing. GO annotation and KEGG analysis suggested that these DE circRNAs probably involved in the energy metabolic regulation. Furthermore, we constructed ceRNA network and the result indicated that these DE circRNAs could regulate lactation through IGF1 and PRL signaling pathway.

m6A mRNA methylation analysis provides novel insights into heat stress responses in the liver tissue of sheep.

N6-methyladenosine (m6A) mRNA methylation varies in response to stress. However, no map of m6A mRNA methylation has been obtained for sheep, nor is it known what effect this has on regulating heat stress in sheep. Here, we obtained m6A methylation maps of sheep liver tissues with and without heat stress by MeRIP-seq. In total, 8306 m6A peaks associated with 2697 genes were detected in the heat stress group, and 12,958 m6A peaks associated with 5494 genes were detected in the control group. Peaks were mainly enriched in coding regions and near stop codons with classical RRACH motifs. Methylation levels of heat stress and control sheep were higher near stop codons, although methylation was significantly lower in heat stress sheep. GO and KEGG revealed that differential m6A-containing genes were significantly enriched in the stress response and fat metabolism. Our results showed that m6A mRNA methylation modifications regulate heat stress in sheep.

Integrative Proteomics and Transcriptomics Profiles of the Oviduct Reveal the Prolificacy-Related Candidate Biomarkers of Goats (Capra hircus) in Estrous Periods.

The oviduct is a dynamic reproductive organ for mammalian reproduction and is required for gamete storage, maturation, fertilization, and early embryonic development, and it directly affects fecundity. However, the molecular regulation of prolificacy occurring in estrous periods remain poorly understood. This study aims to gain a better understanding of the genes involved in regulating goat fecundity in the proteome and transcriptome levels of the oviducts. Twenty female Yunshang black goats (between 2 and 3 years old, weight 52.22 ± 0.43 kg) were divided into high- and low-fecundity groups in the follicular (FH and FL, five individuals per group) and luteal (LH and LL, five individuals per group) phases, respectively. The DIA-based high-resolution mass spectrometry (MS) method was used to quantify proteins in twenty oviducts. A total of 5409 proteins were quantified, and Weighted gene co-expression network analysis (WGCNA) determined that the tan module was highly associated with the high-fecundity trait in the luteal phase, and identified NUP107, ANXA11, COX2, AKP13, and ITF140 as hub proteins. Subsequently, 98 and 167 differentially abundant proteins (DAPs) were identified in the FH vs. FL and LH vs. LL comparison groups, respectively. Parallel reaction monitoring (PRM) was used to validate the results of the proteomics data, and the hub proteins were analyzed with Western blot (WB). In addition, biological adhesion and transporter activity processes were associated with oviductal function, and several proteins that play roles in oviductal communication with gametes or embryos were identified, including CAMSAP3, ITGAM, SYVN1, EMG1, ND5, RING1, CBS, PES1, ELP3, SEC24C, SPP1, and HSPA8. Correlation analysis of proteomics and transcriptomic revealed that the DAPs and differentially expressed genes (DEGs) are commonly involved in the metabolic processes at the follicular phase; they may prepare the oviductal microenvironment for gamete reception; and the MAP kinase activity, estrogen receptor binding, and angiotensin receptor binding terms were enriched in the luteal phase, which may be actively involved in reproductive processes. By generating the proteome data of the oviduct at two critical phases and integrating transcriptome analysis, we uncovered novel aspects of oviductal gene regulation of fecundity and provided a reference for other mammals.

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

BACKGROUND: Litter size is an important index of mammalian prolificacy and is determined by the ovulation rate. The ovary is a crucial organ for mammalian reproduction and is associated with follicular development, maturation and ovulation. However, prolificacy is influenced by multiple factors, and its molecular regulation in the follicular phase remains unclear. METHODS: Ten female goats with no significant differences in age and weight were randomly selected and divided into either the high-yielding group (n = 5, HF) or the low-yielding group (n = 5, LF). Ovarian tissues were collected from goats in the follicular phase and used to construct mRNA and miRNA sequencing libraries to analyze transcriptomic variation between high- and low-yield Yunshang black goats. Furthermore, integrated analysis of the differentially expressed (DE) miRNA-mRNA pairs was performed based on their correlation. The STRING database was used to construct a PPI network of the DEGs. RT-qPCR was used to validate the results of the predicted miRNA-mRNA pairs. Luciferase analysis and CCK-8 assay were used to detect the function of the miRNA-mRNA pairs and the proliferation of goat granulosa cells (GCs). RESULTS: A total of 43,779 known transcripts, 23,067 novel transcripts, 424 known miRNAs and 656 novel miRNAs were identified by RNA-seq in the ovaries from both groups. Through correlation analysis of the miRNA and mRNA expression profiles, 263 negatively correlated miRNA-mRNA pairs were identified in the LF vs. HF comparison. Annotation analysis of the DE miRNA-mRNA pairs identified targets related to biological processes such as "estrogen receptor binding (GO:0030331)", "oogenesis (GO:0048477)", "ovulation cycle process (GO:0022602)" and "ovarian follicle development (GO:0001541)". Subsequently, five KEGG pathways (oocyte meiosis, progesterone-mediated oocyte maturation, GnRH signaling pathway, Notch signaling pathway and TGF-ß signaling pathway) were identified in the interaction network related to follicular development, and a PPI network was also constructed. In the network, we found that CDK12, FAM91A1, PGS1, SERTM1, SPAG5, SYNE1, TMEM14A, WNT4, and CAMK2G were the key nodes, all of which were targets of the DE miRNAs. The PPI analysis showed that there was a clear interaction among the CAMK2G, SERTM1, TMEM14A, CDK12, SYNE1 and WNT4 genes. In addition, dual luciferase reporter and CCK-8 assays confirmed that miR-1271-3p suppressed the proliferation of GCs by inhibiting the expression of TXLNA. CONCLUSIONS: These results increase the understanding of the molecular mechanisms underlying goat prolificacy. These results also provide a basis for studying interactions between genes and miRNAs, as well as the functions of the pathways in ovarian tissues involved in goat prolificacy in the follicular phase.

Expression characteristics of pineal miRNAs at ovine different reproductive stages and the identification of miRNAs targeting the AANAT gene.

BACKGROUND: Many recent studies have shown that miRNAs play important roles in the regulation of animal reproduction, including seasonal reproduction. The pineal gland is a crucial hub in the regulation of seasonal reproduction. However, little is known about the expression characteristics of pineal miRNAs in different reproductive seasons (anestrus and breeding season). Therefore, the expression profiles and regulatory roles of ovine pineal miRNAs were investigated during different reproductive stages using Solexa sequencing technology and dual luciferase reporter assays. RESULTS: A total of 427 miRNAs were identified in the sheep pineal gland. Significant differences in miRNA expression were demonstrated between anestrus and the breeding season in terms of the frequency distributions of miRNA lengths, number of expressed miRNAs, and specifically and highly expressed miRNAs in each reproductive stage. KEGG analysis of the differentially expressed (DE) miRNAs between anestrus and the breeding season indicated that they are significantly enriched in pathways related to protein synthesis, secretion and uptake. Furthermore, transcriptome analysis revealed that many target genes of DE miRNAs in the ribosome pathway showed relatively low expression in the breeding season. On the other hand, analyses combining miRNA-gene expression data with target relationship validation in vitro implied that miR-89 may participate in the negative regulation of aralkylamine N-acetyltransferase (AANAT) mRNA expression by targeting its 3'UTR at a unique binding site. CONCLUSIONS: Our results provide new insights into the expression characteristics of sheep pineal miRNAs at different reproductive stages and into the negative regulatory effects of pineal miRNAs on AANAT mRNA expression.

Combined approaches identify known and novel genes associated with sheep litter size and non-seasonal breeding.

Improvement of ewe reproduction is considerable by appropriately increasing litter size and sustaining non-seasonal breeding. However, their genetic makeups have not been entirely elucidated. Genome-wide analyses of 821 individuals were performed by combining three genomic approaches (genome-wide association study, XP-nSL, and runs of homozygosity). Consequently, 35 candidate genes including three domestication genes (TSHR, GTF2A1, and KITLG) were identified. Other than the FecB mutation at BMPR1B, we described a significant association of a missense mutation rs406686139 at seasonal lambing-associated TSHR gene with litter size. Some promising novel genes may be relevant for sheep reproduction by multitude biological processes, such as FETUB functioning in fertilization, HNRNPA1 in oogenesis, DCUN1D1 in spermatogenesis, and HRG in fertility outcome. The present study suggests that improvement of ewe reproduction is attributed to selective breeding, and casts light on the genetic basis and improvement of sheep reproduction.

Expression and functional analysis of the Follistatin-like 3 (FSTL3) gene in the sheep ovary during the oestrous cycle.

Follistatin-like 3 (FSTL3) is a regulator of cellular apoptosis and was previously identified via RNA-Seq to be associated with follicular development in mammalian ovaries. However, the mechanism underlying the FSTL3 regulation of oestrus in sheep remained poorly understood. In this study, the oestrogen (E2) and progesterone (P4) concentrations in blood were detected, and the expression level and functional analysis of FSTL3 in the ovary were studied during the different reproductive stage in Aohan fine wool sheep (seasonal breeding breed in China). The concentrations of E2 and P4 at the anestrus were significantly lower compared to dioestrus, proestrus and oestrus stages. Higher expression levels of FSTL3 were observed in the sheep ovary, hypothalamus, and thyroid. During different reproductive stages, higher expression levels were found during the stages of dioestrus and proestrus, while lower levels were found during the oestrus and anestrus stages. Functional analysis of FSTL3 was performed in primary granulosa cells (GCs) of sheep. The concentration of E2 increased significantly after RNAi interference of FSTL3, while the P4 level decreased. FSTL3 can decrease P4 levels, which might be involved in mediating oestrous cycle in sheep.

Trend analysis of the role of circular RNA in goat skeletal muscle development.

BACKGROUND: Circular RNA (circRNA) is produced during the splicing of mRNA (in addition to linear splicing) and is part of the gene regulatory network. The temporal expression patterns the different developmental stages were inseparable from these molecules' function. RESULTS: Skeletal muscles of Anhui white goat (AWG) across seven fetal to postnatal development stages were sequenced and 21 RNA sequencing libraries were constructed. We thereby identified 9090 circRNAs and analyzed their molecular properties, temporal expression patterns, and potential functions at the different stages. CircRNAs showed complexities and diversity of formation as the same host gene produces multiple isoforms of these nucleic acids with different expression profiles. The differential expression of 2881 circRNAs (DECs, P < 0.05) was identified and four were randomly selected and validated by qPCR. Moreover, 1118 DECs under strict selected (SDECs, |log2FC| > 2 and P-adj value < 0.01) showed 4 expression trends (Clusters 0, 19, 16 and 18). Cluster 0 molecules had increasing expression at all stages with effects on muscle through metabolism, regulation of enzyme activity, and biosynthesis. Cluster 16 circRNAs had high expression in the early and late stages and are involved in "Wnt signaling pathway", "AMPK signaling pathway" and others. Cluster 18 molecules were mainly expressed at F120 and participate in "cytoskeletal protein binding", "Notch signaling pathway" and so on. Cluster 19 circRNAs were down-regulated at all stages and related to muscle structure and development. Lastly, the SDECs divided the period of skeletal muscle development into three transitional stages: stage 1 (F45 to F90), which related to muscle satellite cell proliferation and muscle fiber structure; stage 2 (F90 to B1), in which the attachment of the cytoplasmic surface to the actin cytoskeleton initiates; and stage 3, which involved the "cGMP-PKG signaling pathway". Moreover, the paraffin sections messages also validated that there are three transitional stages of skeletal muscle development. CONCLUSION: Our current study provides a catalog of goat muscle-related circRNAs that can stratify skeletal muscle development fetus 45 days to newborn 90 days into three developmental stages. These findings better our understanding of functional transitions during mammalian muscle development.