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.

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.

Screening of Litter-Size-Associated SNPs in NOX4, PDE11A and GHR Genes of Sheep

In previous studies, NOX4, PDE11A and GHR genes have been screened as important candidate genes for litter size in sheep by using the GWAS method; however, neither their effects on litter size nor the loci associated with litter size have been identified. In this study, three candidate loci (c.1057-4C > T in NOX4, c.1983C > T in PDE11A and c.1618C > T in GHR) were first screened based on our previous resequencing data of 10 sheep breeds. After the three loci were genotyped using Sequenom MassARRAY technology, we carried out population genetics analysis on the three loci and performed association analysis between the polymorphism of the three loci and the litter size of sheep. The results of population genetics analysis suggested that c.1057-4C > T in NOX4 and c.1983C > T in PDE11A may be subject to natural or artificial selection. The results of association analysis indicated that litter size was significantly associated with c.1057-4C > T in NOX4 and c.1983C > T in PDE11A (p < 0.05) in Small Tail Han sheep, and there was no significant interaction effect between the two loci on the litter size. In summary, c.1057-4C > T in NOX4 and c.1983 C > T in PDE11A can be considered candidate molecular markers for improving litter size in sheep.

Transcriptomic Analysis Reveals Differentially Expressed Circular RNAs Associated with Fecundity in the Sheep Hypothalamus with Different FecB Genotypes

Circular RNAs (circRNAs) are a specific type of noncoding RNA, and some have defined roles in cellular and biological processes. However, little is known about the role of circRNAs in follicular development in sheep with FecB (fecundity Booroola) mutations. Here, the expression profiles of circRNAs were investigated using RNA sequencing (RNA-seq) in the follicular phase (F) and the luteal phase (L) of FecB mutant homozygous (BB) and wild-type (WW) Small Tail Han sheep. A total of 38,979 circRNAs were identified, and 314, 343, 336, and 296 of them were differentially expressed (DE) between BB_F and BB_L, WW_F and WW_L, BB_F and WW_F, and BB_L and WW_L, respectively. The length, type, and chromosome distribution of the circRNAs and the expression characteristic between the circRNAs and their host genes in the sheep hypothalamus were ascertained. Enrichment analysis showed that the host genes of DE circRNAs in the follicular and luteal phases were annotated to MAPK, gap junctions, progesterone-mediated oocyte maturation, oocyte meiosis, and other hormone-related signaling pathways, and the different FecB genotypes were annotated to the gap junctions, circadian entrainment, MAPK, and other hormone-related signaling pathways. The competing endogenous RNA network prediction revealed that the 129 target miRNAs might be bound to 336 DE circRNAs. oar_circ_0000523 and oar_circ_0028984, which were specifically expressed during the follicular phase in the BB genotype sheep, probably acted as miRNA sponges involved in the regulation of LH synthesis and secretion. This study reveals the expression profiles and characterization of circRNAs at two phases of follicular development considering different FecB genotypes, thereby providing an improved understanding of the roles of circRNAs in the sheep hypothalamus and their involvement in follicular development and ovulation.

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.

Proteomic Analysis Identifies Distinct Protein Patterns for High Ovulation in FecB Mutant Small Tail Han Sheep Granulosa Cells

The Booroola fecundity (FecB) mutation in the bone morphogenetic protein receptor type 1B (BMPR1B) gene increases ovulation in sheep. However, its effect on follicular maturation is not fully understood. Therefore, we collected granulosa cells (GCs) at a critical stage of follicle maturation from nine wild-type (WW), nine heterozygous FecB mutant (WB), and nine homozygous FecB mutant (BB) Small Tail Han sheep. The GCs of three ewes were selected at random from each genotype and consolidated into a single group, yielding a total of nine groups (three groups per genotype) for proteomic analysis. The tandem mass tag technique was utilized to ascertain the specific proteins linked to multiple ovulation in the various FecB genotypes. Using a general linear model, we identified 199 proteins significantly affected by the FecB mutation with the LIMMA package (p < 0.05). The differential abundance of proteins was enriched in pathways related to cholesterol metabolism, carbohydrate metabolism, amino acid biosynthesis, and glutathione metabolism. These pathways are involved in important processes for GC-regulated 'conservation' of oocyte maturation. Further, the sparse partial least-squares discriminant analysis and the Fuzzy-C-mean clustering method were combined to estimate weights and cluster differential abundance proteins according to ovulation to screen important ovulation-related proteins. Among them, ZP2 and ZP3 were found to be enriched in the cellular component catalog term "egg coat", as well as some apolipoproteins, such as APOA1, APOA2, and APOA4, enriched in several Gene Ontology terms related to cholesterol metabolism and lipoprotein transport. A higher abundance of these essential proteins for oocyte maturation was observed in BB and WB genotypes compared with WW ewes. These proteins had a high weight in the model for discriminating sheep with different FecB genotypes. These findings provide new insight that the FecB mutant in GCs improves nutrient metabolism, leading to better oocyte maturation by altering the abundance of important proteins (ZP2, ZP3, and APOA1) in favor of increased ovulation or better oocyte quality.

Polymorphism and expression of GLUD1 in relation to reproductive performance in Jining Grey goats

Understanding the molecular mechanism of mammalian reproduction (puberty and prolificacy) will play a part in improving animal reproductive performance. GLUD1 (glutamate dehydrogenase 1) is important for mammalian reproduction, as shown in previous studies; however, its roles in puberty and prolificacy have rarely been reported. In this study, we designed seven pairs of primers (P1 to P7) for cloning and sequencing genomic DNA of Jining Grey goats and Liaoning Cashmere goats. Primer 8 (P8) was designed to detect single nucleotide polymorphism (SNP) of the GLUD1 in both sexually precocious and high-fecundity breeds (Jining Grey, Nanjiang Brown and Matou goats) and sexually late-maturing and low-fecundity breeds (Liaoning Cashmere, Inner Mongolia Cashmere and Taihang goats) by PCR-RFLP (restriction fragment length polymorphism). The real-time quantitative polymerase chain reaction (RT-qPCR) technique was used to detect the expression of GLUD1 in a variety of tissues. The results showed that the A197C mutation was only found in the amplification product of P6. For this SNP locus, only two genotypes (AA and AC) were detected in Nanjiang Brown goats, while three genotypes (AA, AC and CC) were detected in the other five breeds. In Jining Grey goats, the frequency of genotypes AA, AC and CC was 0.69, 0.26 and 0.05, respectively. In Jining Grey goats, AA genotype had 0.54 (P < 0.05 ) and 0.3 (P < 0.05 ) more kids than the CC and AC genotype, respectively, and no significant difference (P > 0.05 ) was found in kidding number between the AC and CC genotype. GLUD1 was expressed in five tissues of different developmental stages. The expression level of GLUD1 in the hypothalamus was higher than that in the other four tissues except during puberty of Liaoning Cashmere goats. In puberty in goats, GLUD1 expression was significantly higher in ovaries than that in the juvenile period (P < 0.01 ). RT-qPCR results showed that the expression of GLUD1 in ovaries may relate to the puberty of goats. The present study preliminarily indicated that there might be an association between the 197 locus of GLUD1 and sexual precocity in goats, and allele A of GLUD1 was a potential DNA marker for improving kidding number in Jining Grey goats.

Thyroid Transcriptomics Revealed the Reproductive Regulation of miRNA in the Follicular and Luteal Phases in Small-Tail Han Sheep with Different FecB Genotypes

MicroRNA (miRNA) is a type of endogenous short-stranded ncRNA that influences many biological processes such as animal growth, development and metabolism. The thyroid gland is an important endocrine gland in sheep, and an increasing number of studies have shown that the thyroid gland plays an important role in animal reproduction, but the molecular mechanisms of the thyroid gland in sheep reproduction are poorly understood. In this study, RNA-seq was used to detect transcriptome expression patterns in the thyroid gland between the follicular phase (FP) and luteal phase (LP) in FecB BB (MM) and FecB ++ (ww) small-tail Han (STH) sheep, respectively, and to identify differentially expressed miRNAs (DEMs) associated with reproduction. Bioinformatic analysis of the target genes of these DEMs revealed that they can be enriched in multiple GO terms associated with the reproductive process in animals and in the KEGG signaling pathway. The miRNA-mRNA coexpression network revealed that oar-miR-133 and oar-miR-370-3p may play an important role in sheep reproduction. The results of the dual-luciferase reporter assay suggest a possible targeting relationship between novel-51 and TARBP2. These results provided a novel resource for elucidating regulatory mechanisms underlying STH sheep prolificacy.

Comparative Transcriptomics Identify Key Pituitary Circular RNAs That Participate in Sheep (Ovis aries) Reproduction

CircRNAs have been found to play key roles in many biological processes and have diverse biological functions. There have been studies on circRNAs in sheep pituitary, and some important circRNAs have been found. But there are still few studies on circRNAs in sheep pituitary with different fecundity. In this study, we obtained the circRNAs expression profiles in the pituitary of FecB ++ genotype Small Tail Han sheep with different fecundity and estrous phases. A total of 34,878 circRNAs were identified in 12 pituitary samples, 300 differentially expressed circRNAs (DE circRNAs) (down: 104; up: 196) were identified in polytocous sheep in the follicular phase (PF) and monotocous sheep in the follicular phase (MF) (PF vs. MF), and 347 DE circRNAs (down: 162; up: 185) were identified in polytocous sheep in the luteal phase (PL) and monotocous sheep in the luteal phase (ML) (PL vs. ML). Cortisol synthesis and secretion pathway (follicular phase) and estrogen signaling pathway (luteal phase) were obtained by functional enrichment analysis of circRNAs source genes. Competing endogenous RNA (ceRNA) network analysis of key DE circRNAs revealed that oar-circ-0022776 (source gene ITPR2, follicular phase) targeted oar-miR-432, oar-circ-0009003 (source gene ITPR1, luteal phase) and oar-circ-0003113 (source gene PLCB1, luteal phase) targeted oar-miR-370-3p. We also explored the coding ability of DE circRNAs. In conclusion, our study shows that changes in the pituitary circRNAs may be related to the response of the pituitary to steroid hormones and regulate the reproductive process of sheep by affecting the pituitary function. Results of this study provide some new information for understanding the functions of circRNAs and the fecundity of FecB ++ genotype sheep.

Comparative proteomics of ovaries elucidated the potential targets related to ovine prolificacy

Small Tail Han (STH) sheep, a unique Chinese breed, is recognized for its early maturity, year-round estrus, and prolificacy. However, the molecular mechanism of its high prolificacy has not been fully elucidated. The Proteomics approach is feasible and effective to reveal the proteins involved in the complex physiological processes of any organism. Given this, we performed the protein expression profiling of ovarian tissues during the luteal phase using polytocous STH sheep (litter size ≥2, three consecutive lambings) and monotocous STH sheep (litter size =1, three consecutive lambings) (PL vs. ML), and the follicular phase using polytocous STH sheep (litter size ≥2, three consecutive lambings) and monotocous STH sheep (litter size =1, three consecutive lambings) (PF vs. MF), respectively. Parallel Reaction Monitoring (PRM) was conducted to validate the differentially abundant proteins (DAPs). The tandem mass tag (TMT) quantitative proteomic results showed that a total of 5,237 proteins were identified, of which 49 and 44 showed differential abundance in the PL vs. ML and PF vs. MF groups, respectively. Enrichments analyses indicated that the DAPs including TIA1 cytotoxic granule-associated RNA-binding protein-like 1 (TIAL1), nicotinamide phosphoribosyltransferase (NAMPT), and cellular retinoic acid-binding protein 1 (CRABP1) were enriched at the luteal phase, while TIAL1, inhibin beta-a-subunit (A2ICA4), and W5PG55 were enriched at the follicular phase, potentially mediating reproductive processes in polytocous ewes. Furthermore, six DAPs were verified using PRM, confirming the accuracy of the TMT data acquired in this study. Together, our work expanded the database of indigenous sheep breeds and provided new ovarian candidate molecular targets, which will help in the study of the genetic mechanisms of ovine prolificacy.

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 (Fec^B BB) and wild type (Fec^B ++) mutations. GCs collected from small-tailed Han sheep with different Fec^B genotypes were typed using a TaqMan probe assay, and autophagy levels were found to be significantly higher in GCs with Fec^B BB than the levels in those with Fec^B ++. Autophagy-related 2 homolog B (ATG2B) was associated with cell autophagy and was highly expressed in GCs with the Fec^B BB genotype in small-tailed Han sheep. Overexpression of ATG2B in the GCs of sheep with both Fec^B genotypes promoted GC autophagy, and the contrary was observed after the inhibition of ATG2B expression. Subsequently, treatment of GCs with different genotypes of Fec^B 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 Fec^B ++ genotype. In conclusion, this study determined that autophagy levels were significantly higher in sheep GCs with Fec^B BB genotype than the levels in those with Fec^B ++ genotype, which may have contributed to the difference in lambing numbers between the two Fec^B 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. This article is protected by copyright. All rights reserved.

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.

The transcription factor RUNX1 affects the maturation of porcine oocytes via the BMP15/TGF-β signaling pathway

Bone morphogenetic protein 15 (BMP15) is specifically expressed in oocytes in pigs at all stages from early stages to ovulation and has an important role in oocyte maturation. However, there are few reports on the molecular mechanisms by which BMP15 affects oocyte maturation. In this study, we identified the core promoter region of BMP15 using a dual luciferase activity assay and successfully predicted the DNA binding motif of the transcription factor RUNX1. The effect of BMP15 and RUNX1 on oocyte maturation was examined using the first polar body extrusion rate, a reactive oxygen species (ROS) assay and total glutathione (GSH) content at three time points of 12, 24 and 48 h of in vitro culture of porcine isolated oocytes. Subsequently, the effect of the transcription factor RUNX1 on the TGF-β signaling pathway (BMPR1B and ALK5) was further verified using RT-qPCR and Western blotting. We found that the overexpression of BMP15 significantly increased the first polar body extrusion rate (P < 0.01) and total glutathione content of oocytes cultured in vitro for 24 h and decreased reactive oxygen levels (P < 0.01), whereas interference with BMP15 decreased the first polar body extrusion rate (P < 0.01), increased reactive oxygen levels in oocytes cultured in vitro for 24 h (P < 0.01), and decreased glutathione content (P < 0.01). The dual luciferase activity assay and online software prediction showed that RUNX1 is a potential transcription factor binding to the core promoter region (-1203/-1423 bp) of BMP15. Overexpression of RUNX1 significantly increased the expression of BMP15 and oocyte maturation rate, while inhibition of RUNX1 decreased the expression of BMP15 and the oocyte maturation rate. Moreover, the expression of BMPR1B and ALK5 in the TGF-β signaling pathway increased significantly after overexpression of RUNX1, whereas their expression decreased after inhibition of RUNX1. Overall, our results suggest that the transcription factor RUNX1 positively regulates the expression of BMP15 and influences oocyte maturation through the TGF-β signaling pathway. This study provides a theoretical basis for further complementing the BMP15/TGF-β signaling pathway to regulate mammalian oocyte maturation.

Integrated Proteotranscriptomics of the Hypothalamus Reveals Altered Regulation Associated with the FecB Mutation in the BMPR1B Gene That Affects Prolificacy in Small Tail Han Sheep

The litter size and ovulation rate are different among ewes of different FecB genotypes in Small Tail Han sheep. These variants in reproductive phenotypes may be regulated by hormones released by the hypothalamic-pituitary-ovarian axis. However, there have been few reports on the hypothalamus regarding regulating an increase in ovulation in sheep with FecB mutation at different estrous stages. Thus, we examined the abundance of hypothalamus tissue protein profiles of six FecB mutant homozygous (BB) and six wild-type (WW) ewes at the luteal and follicular phases. We determined this abundance by tandem mass tag-based quantitative analysis and parallel reaction monitoring methods. Furthermore, an integrated proteotranscriptomic analysis was performed by the Data Integration Analysis for Biomarker discovery using the latent variable approaches for Omics studies (DIABLO) framework to examine biological processes and pathway alterations by the FecB mutant. The abundance of 154 proteins was different between the two estrous stages. Growth hormone and prolactin were particularly enriched in the neuroactive ligand-receptor interactions, the prolactin signaling pathway, and the PI3K-Akt signaling pathway which are related to hypothalamic function and reproduction. We combined proteome and transcriptome data from different estrous stages and genotypes. There is a high correlation (Pearson correlation coefficient = 0.99) between the two datasets in the first two components. We applied the traditional single-omic multivariate approach to obtain differentially abundant proteins and differentially expressed genes. The major fertility related biomarkers were selected using the two approaches mentioned above. Several key pathways (GABAergic synapse, neuroactive ligand-receptor interaction, estrogen and MAPK signaling pathways) were enriched, which are central to gonadotrophin-releasing hormone (GnRH) secretion and reproduction. A higher level of gamma-aminobutyric acid type A receptor subunit alpha1 (GABRA1) and gamma-aminobutyric acid type A receptor subunit beta2 (GABRB2) expression was observed in BB ewes as compared to WW ewes. This finding suggested that a greater production of GnRH during follicular development in BB ewes may explain the higher mature follicle number in mutant ewes. FKBP prolyl isomerase 1A (FKBP1A), which was a major feature factor in the proteome selected by DIABLO, was an important switch for activating the transforming growth factor beta (TGFβ) pathway, and its expression was higher in the WW ewes than in the BB ewes. We suggest that BB sheep maintain TGFβ pathway activity by reducing FKBP1A protein abundance. This innovative data integration in the hypothalamus may provide fresh insight into the mechanisms by which the FecB mutation affects sheep fertility, while providing novel biomarkers related to reproductive endocrinology in sheep breeding.

Characteristics of piRNAs and their comparative profiling in testes of sheep with different fertility

As a novel class of small RNAs, piRNAs are highly expressed in the animal gonads and their main known role is to inhibit transposon activity for ensuring the correctness and integrity of genome. In order to explore the characteristics of piRNAs in sheep testis and their possible regulatory roles on male reproduction, deep sequencing technology was used to sequence small RNAs and identify piRNAs in testes of sheep. The length of piRNAs in sheep testes showed a unimodal distribution between 26 and 31 nt, with a peak at 29 nt. These piRNAs exhibited obvious ping-pong signature and strand specificity. In the genome, they were mainly aligned to CDS, intron, repetitive sequence regions and unannotated regions. Furthermore, in transposon analysis, piRNAs were aligned predominantly to LINE, SINE, and LTR types of retrotransposon in sheep testes, and the piRNAs derived from each type showed obvious ping-pong signature. The piRNA clusters identified in sheep testes were mainly distributed on chromosomes 3, 7, 15, 17, 18 and 20. The results combining semen determination with pathway enrichment analysis implied that differentially expressed piRNAs between the testes of rams with different fertility might participate in spermatogenesis by regulating multiple pathways closely related to stabilization of blood-testis barrier and renewal and differentiation of spermatogonial stem cell. Taken together, the study provided new insights into the characteristics, origin and expression patterns of piRNAs in sheep testes tissue, which would help us better understand the role of piRNAs in sheep reproduction.

Characterization of circular RNA profiles of oviduct reveal the potential mechanism in prolificacy trait of goat in the estrus cycle

The mammalian oviduct is functionally highly diverse during the estrus cycle. It provides a suitable milieu for oocyte maturation, sperm capacitation, fertilization, early embryo development and transportation. While there have been many studies of molecular mechanisms on the kidding number of goats, a systematic analysis by which the underlying circular RNAs (circRNAs) changes in the oviduct related to prolificacy traits is lacking. Herein, we present a comprehensive circRNA atlas of the oviduct among high- and low-fecundity goats in the follicular phase (FH vs. FL), luteal phase (LH vs. LL), and estrus cycle (FH vs. LH; FL vs. LL) to unravel their potential regulatory mechanisms in improving kidding number. We generated RNA sequencing data, and identified 4,078 circRNAs from twenty sampled Yunshang black goats. Many of these circRNAs are exon-derived and differentially expressed between each comparison group. Subsequently, eight differentially expressed (DE) circRNAs were validated by RT‒qPCR, which was consistent with the RNA-seq data. GO and KEGG enrichment analyses suggested that numerous host genes of DE circRNAs were involved in the hormone secretion, gamete production, fertilization, and embryo development processes. The competing endogenous RNA (ceRNA) interaction network analysis revealed that 2,673 circRNA–miRNA–mRNA axes (including 15 DE circRNAs, 14 miRNAs, and 1,699 mRNAs) were formed, and several target genes derived from the ceRNA network were associated with oviduct functions and reproduction, including SMAD1, BMPR1B, IGF1, REV1, and BMP2K. Furthermore, miR-15a-5p, miR-181b-5p, miR-23b-5p, miR-204-3p, and miR-145-5p might play important roles in reproduction. Finally, a novel circRNA, circIQCG, was identified as potentially involved in embryo development. Overall, our study provides a resource of circRNAs to understand the oviductal function and its connection to prolificacy trait of goats in the differentiation estrus cycle.

Expression characteristics of piRNAs in ovine luteal phase and follicular phase ovaries

PIWI-interacting RNAs (piRNAs), as a novel class of small non-coding RNAs that have been shown to be indispensable in germline integrity and stem cell development. However, the expressed characteristics and regulatory roles of piRNAs during different reproductive phases of animals remain unknown. In this study, we investigated the piRNAs expression profiles in ovaries of sheep during the luteal phase (LP) and follicular phase (FP) using the Solexa sequencing technique. A total of 85,219 and 1,27,156 piRNAs tags were identified in ovine ovaries across the two phases. Most expressed piRNAs start with uracil. piRNAs with a length of 24 nt or 27–29 nts accounted for the largest proportion. The obvious ping-pong signature appeared in the FP ovary. The piRNA clusters in the sheep ovary were unevenly distributed on the chromosomes, with high density on Chr 3 and 1. For genome distribution, piRNAs in sheep ovary were mainly derived from intron, CDS, and repeat sequence regions. Compared to the LP ovary, a greater number of expressed piRNA clusters were detected in the FP ovary. Simultaneously, we identified 271 differentially expressed (DE) piRNAs between LP and FP ovaries, with 96 piRNAs upregulated and 175 piRNAs downregulated, respectively. Functional enrichment analysis (GO and KEGG) indicated that their target genes were enriched in reproduction-related pathways including oocyte meiosis, PI3K-Akt, Wnt, and TGF-β signaling pathways. Together, our results highlighted the sequence and expression characteristics of the piRNAs in the sheep ovary, which will help us understand the roles of piRNAs in the ovine estrus cycle.

Identification of photoperiod-induced specific miRNAs in the adrenal glands of Sunite sheep (Ovis aries)

In seasonal estrus, it is well known that melatonin-regulated biorhythm plays a key role. Some studies indicate that the adrenal gland plays an important role in reproduction in mammals, but the molecular mechanism is not clear. This study used an artificially controlled light photoperiod model, combined with RNA-seq technology and bioinformatics analysis, to analyze the messenger RNA (mRNA) and microRNA (miRNA) of ewe (Sunite) adrenal glands under different photoperiod treatments. After identification, the key candidate genes GRHL2, CENPF, FGF16 and SLC25A30 that photoperiod affects reproduction were confirmed. The miRNAs (oar-miR-544-3p, oar-miR-411b-5p, oar-miR-376e-3p, oar-miR-376d, oar-miR-376b-3p, oar-miR-376a-3p) were specifically expressed in the adrenal gland. The candidate mRNA-miRNA pairs (e.g., SLC25A30 coagulated by novel miRNA554, novel miRNA555 and novel miRNA559) may affect seasonal estrus. In summary, we constructed relation network of the mRNAs and miRNAs of sheep adrenal glands using RNA sequencing and bioinformatics analysis, thereby, providing a valuable genetic variation resource for sheep genome research, which will contribute to the study of complex traits 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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08657-8.

Oar-miR-432 Regulates Fat Differentiation and Promotes the Expression of BMP2 in Ovine Preadipocytes

The fat tail is a unique characteristic of sheep that represents energy reserves and is a complex adaptative mechanism of fat-tailed sheep to environmental stress. MicroRNA plays a significant role as regulators at the posttranscriptional level, but no studies have explained the molecular mechanisms of miRNA which regulate fat deposition in sheep tails. In this study, mRNA and miRNA analysis examined tail fat tissue from three Hu fat-tailed and three Tibetan thin-tailed sheep. After aligning to the reference sequences, 2,108 differentially expressed genes and 105 differential expression miRNAs were identified, including 1,247 up- and 861 downregulated genes and 43 up- and 62 downregulated miRNAs. Among these differentially expressed miRNAs, oar-miR-432 was one of the most downregulated miRNAs between Hu sheep and Tibetan sheep, and 712 genes were predicted to be targeted by oar-miR-432, 80 of which overlapped with DEGs. The Gene Ontology analysis on these genes showed that BMP2, LEP, GRK5, BMP7, and RORC were enriched in fat cell differentiation terms. The genes for BMP2 targeted by oar-miR-432 were examined using dual-luciferase assay. The oar-miR-432 mimic transfected into preadipocytes resulted in increased expression of BMP2. The marker gene PPAR-γ of fat differentiation had a lower expression than the negative control on days 0, 2, and 4 after induced differentiation. The decrease in the number of lipids in the oar-miR-432 mimic group detected by oil red O stain was also less than that in the negative control. This is the first study to reveal the fat mechanisms by which oar-miR-432 inhibits fat differentiation and promotes the expression of BMP2 in sheep tails.

Thyroid Transcriptomic Profiling Reveals the Follicular Phase Differential Regulation of lncRNA and mRNA Related to Prolificacy in Small Tail Han Sheep with Two FecB Genotypes

Long non-coding RNA (lncRNA) accounts for a large proportion of RNA in animals. The thyroid gland has been established as an important gland involved in animal reproduction, however, little is known of its gene expression patterns and potential roles in the sheep. Herein, RNA-Seq was used to detect reproduction-related differentially expressed lncRNAs (DELs) and mRNAs (DEGs) in the follicular phase (FT) FecBBB (MM) and FecB++ (ww) genotypes of Small Tail Han (STH) sheep thyroids. Overall, 29 DELs and 448 DEGs in thyroid between MM and ww sheep were screened. Moreover, GO and KEGG enrichment analysis showed that targets of DELs and DEGs were annotated in biological transitions, such as cell cycle, oocyte meiosis and methylation, which in turn affect reproductive performance in sheep. In addition, we constructed co-expression and networks of lncRNAs-mRNAs. Specifically, XLOC_075176 targeted MYB, XLOC_014695 targeted VCAN, 106991527 targeted CASR, XLOC_075176 targeted KIFC1, XLOC_360232 targeted BRCA2. All these differential lncRNAs and mRNAs expression profiles in the thyroid provide a new resource for elucidating the regulatory mechanism underlying STH sheep prolificacy.

Transcriptomic Changes of Photoperiodic Response in the Hypothalamus Were Identified in Ovariectomized and Estradiol-Treated Sheep

Accurate timing of seasonal changes is an essential ability for an animal’s survival, and the change in the photoperiod is the key factor affecting reproductive seasonality in mammals. Emerging evidence has suggested that multiple hypothalamic genes participate in the photoperiod-induced regulation of reproductive activities in sheep, but the mechanism is still unclear. In this study, we initially examined the plasma level of two major reproductive hormones, namely, follicle-stimulating hormone (FSH) and prolactin (PRL), under different photoperiods in ovariectomized and estradiol-treated (OVX + E2) sheep using radioimmunoassay (RIA). Of the two hormones, the concentration of PRL significantly increased with the extension of the photoperiod, while FSH showed the opposite trend. Subsequently, an examination of the transcriptomic variation between the short photoperiod (SP) and long photoperiod (LP) was conducted. Differential expression analyses and functional annotation showed that several key genes in the insulin secretion (VAMP2, PRKACB, PRKCG, and PLCB1), GnRH (MAPK13, CGA, CDC42, ATF4, and LHB) pathways, and circadian entrainment (KCNJ5, PER1, GNB2, MTNR1A, and RASD1), as well as numerous lncRNAs, including XR_173257.3, XR_173415.3, XR_001435315.1, XR_001024596.2, and XR_001023464.2, were shown potentially vital for the hypothalamic photoperiodic response. Four of the differentially expressed mRNAs and lncRNAs were validated by qPCR. The constructed mRNA–mRNA interaction networks further revealed that transcripts potentially participated in hypothalamic thyroid hormone synthesis, endocrine resistance, and neuroactive ligand–receptor interactions. The interactome analysis of lncRNAs and their targets implied that XR_173257.3 and its target arylalkylamine N-acetyltransferase (AANAT) and XR_173415.3 and its target TH might participate in the regulation of seasonal reproduction. Together, the changes in reproductive hormones and transcriptome will help to determine the important photoperiod-induced lncRNAs and mRNAs and provide a valuable resource for further research on reproductive seasonality in sheep.

Comparative Transcriptomics Reveals the Key lncRNA and mRNA of Sunite Sheep Adrenal Gland Affecting Seasonal Reproduction

The hypothalamic–pituitary–adrenal (HPA) axis plays an important role in the growth and development of mammals. Recently, lncRNA transcripts have emerged as an area of importance in sheep photoperiod and seasonal estrus studies. This research aims to identify lncRNA and mRNA that are differentially expressed in the sheep adrenal gland in long (LP) or short (SP) photoperiods using transcriptome sequencing and bioinformatics analysis based on the OVX + E₂ (Bilateral ovariectomy and estradiol-implanted) model. We found significant differences in the expression of lncRNAs in LP42 (where LP is for 42 days) vs. SP-LP42 (where SP is for 42 days followed by LP for 42 days) (n = 304), SP42 (where SP is for 42 days) vs. SP-LP42 (n = 1,110) and SP42 vs. LP42 (n = 928). Cluster analysis and enrichment analysis identified SP42 vs. LP42 as a comparable group of interest and found the following candidate genes related to reproductive phenotype: FGF16, PLGF, CDKN1A, SEMA7A, EDG1, CACNA1C and ADCY5. FGF16 (Up-regulated lncRNA MSTRG.242136 and MSTRG.236582) is the only up-regulated gene that is closely related to oocyte maturation. However, EDG1 (Down-regulated lncRNA MSTRG.43609) and CACNA1C may be related to precocious puberty in sheep. PLGF (Down-regulated lncRNA MSTRG.146618 and MSTRG.247208) and CDKN1A (Up-regulated lncRNA MSTRG.203610 and MSTRG.129663) are involved in the growth and differentiation of placental and retinal vessels, and SEMA7A (Up-regulated lncRNA MSTRG.250579) is essential for the development of gonadotropin-releasing hormone (GnRH) neurons. These results identify novel candidate genes that may regulate sheep seasonality and may lead to new methods for the management of sheep reproduction. This study provides a basis for further explanation of the basic molecular mechanism of the adrenal gland, but also provides a new idea for a comprehensive understanding of seasonal estrus characteristics in Sunite sheep.

P2X7 receptor knockdown suppresses osteoclast differentiation by inhibiting autophagy and Ca2+/calcineurin signaling

Bone is continuously remodeled in a dynamic process maintained by osteoclasts and osteoblasts, and imbalances in the relative activities of these cell types can cause various pathological conditions, including rheumatoid arthritis and osteoporosis. Osteoclasts are multinucleated cells that serve an important role in regulating the development of osteoporosis. Furthermore, P2X7 receptor activation has a vital role in physiological and pathological reactions in bone, including bone disease. Therefore, the present study aimed to investigate the effect of P2X7 receptor on osteoclast differentiation and to explore the underlying molecular mechanism by western blotting and tartrate‑resistant acid phosphatase staining. The results indicated that the expression levels of P2X7 receptor and intracellular Ca2+ concentration levels were very high in mature osteoclasts. Furthermore, P2X7 receptor overexpression increased the number of multinucleated osteoclasts and the expression of osteoclastogenesis‑related proteins. P2X7 receptor overexpression was also associated with downstream activation of Ca²⁺/calcineurin/nuclear factor of activated T cells c1 (NFATc1) signaling and increased expression of autophagy‑related proteins during osteoclast differentiation. By contrast, knockdown of P2X7 receptor exerted the opposite effects. Notably, FK506 (a Ca²⁺/calcineurin/NFATc1 signaling inhibitor) abrogated P2X7 receptor overexpression‑induced osteoclast differentiation and activation of autophagy. Moreover, 3‑MA (an autophagy inhibitor) significantly suppressed P2X7 receptor overexpression‑induced osteoclast differentiation. In conclusion, P2X7 receptor knockdown may suppress osteoclast differentiation by modulating autophagy and the Ca²⁺/calcineurin/NFATc1 signaling pathway.

Photoperiods induced the circRNA differential expression in the thyroid gland of OVX+E2 ewes

Circular RNAs (circRNAs) are non-coding RNAs newly identified and play important roles in RNA regulation. However, little is known regarding photoperiods induced circRNAs in the thyroid gland. In this study, we performed a comprehensive analysis of circRNA profiles in the thyroid gland of OVX+E₂ ewes at different photoperiods by whole transcriptome sequencing. A total of 37,470 novel circRNAs were detected in different photoperiods (42 days of short photoperiod treatment, SP42; 42 days of long photoperiod treatment, LP42; SP42 transfer to LP42, SPLP42), with a total of 817 circRNAs for SP42-LP42 (down: 132; up: 114), LP42-SPLP42 (down: 136; up: 112) and SP42-SPLP42 (down: 182; up: 141) having differentially expressed. Functional enrichment annotation analysis of DE-circRNAs for GO and KEGG by R package, features that influence photoperiod response in Sunite ewes through the Inositol phosphate metabolism, cGMP-PKG signaling pathway, Calcium signaling pathway, MAPK signaling pathway, and Oocyte meiosis. In addition, competitive endogenous RNA (ceRNA) network analysis revealed target binding sites for identified miRNAs in DE-cirRNAs such as oar-miR-10b, oar-miR-200c, oar-miR-21, oar-miR-370-3p, oar-miR-377-3p, oar-miR-181a, oar-miR-432, and oar-miR-495-3p. These results of this study will provide some new information for understanding circRNA function as well as the changes in the sheep thyroid gland under different photoperiods.

Single-nucleotide polymorphisms in <i>FLT3</i>, <i>NLRP5</i>, and <i>TGIF1</i> are associated with litter size in Small-tailed Han sheep

Previous studies have indicated that FLT3, NLRP5, and TGIF1 play a pivotal role in sheep fecundity. Nevertheless, little is known about the association of the polymorphisms of these genes with litter size (LS). In this study, the selected single-nucleotide polymorphisms (SNPs) were genotyped using a Sequenom MassARRAY® platform, and the distribution of different genotypes of the SNPs in the seven sheep breeds (Small-tailed Han, Hu, Cele Black, Suffolk, Tan, Prairie Tibetan, and Sunite sheep) were analyzed. The reliability of the estimated allele frequency for all seven SNPs was at least 0.9545. Given the association of the TGIF1 g.37866222C > T polymorphism with LS in Small-tailed Han sheep (p<0.05), fecundity differences might be caused by the change in amino acid from proline (Pro) to serine (Ser), which has an impact on secondary, tertiary protein structures with concomitant TGIF1 functionality changes. The FLT3 rs421947730 locus has a great effect on the LS (p<0.05), indicating that the locus of FLT3 in synergy with KILTG is likely to facilitate ovarian follicle maturation and ovulation. Moreover, NLRP5 rs426897754 is associated with the LS of the second and third parities (p<0.05). We speculate that a synonymous variant of NLRP5 may be involved in folliculogenesis accompanied by BMP15, FSHR, BMPR1B, AMH, and GDF9, resulting in the different fecundity of Small-tailed Han sheep. Our studies provide valuable genetic markers for sheep breeding.

Detection of Novel Variations Related to Litter Size in BMP15 Gene of Luzhong Mutton Sheep (Ovis aries)

Simple Summary

BMP15 is a critical gene in sheep reproduction. Most of its variations have been reported in European sheep. In this study, the entire open reading frame (ORF) region of BMP15 was sequenced in 154 Luzhong mutton sheep. Among 13 identified variations, six were novel. Four SNPs (ENSOART00000010201.1:c.352+342C>A, c.352+1232T>C, c.352+1165A>G and c.353-2036T>A) were significantly associated with litter size, and could be used as candidate genetic markers for improving litter size. The results also suggested possible interaction between BMP15 and FecB/GDF9.

Abstract

Litter size is an important economic trait in the mutton sheep industry. BMP15 is one of the key candidate genes for litter size in sheep. In this study, the entire ORF region of BMP15 was sequenced in 154 Luzhong mutton ewes, and the novel variations were determined. The association between polymorphism in BMP15 and litter size was analyzed using a general linear model. Six out of a total of thirteen variations were identified to be novel. Association analysis indicated that four (SNPs ENSOART00000010201.1:c.352+342C>A, c.352+1232T>C, c.352+1165A>G and c.353-2036T>A) were significantly associated with litter size. The joint analysis among three major genes (BMP15, BMPR1B and GDF9) exhibited significant interaction effects in three combinations (FecB and c.352+1232T>C of BMP15; FecB and c.352+1165A>G of BMP15; c.352+342C>A of BMP15 and ENSOART00000014382.1:c.994G>A of GDF9). For the SNPs c.352+1232T>C and c.352+342C>A, the global distribution of allele frequencies showed that the highest variation frequency occurs in Western Europe. In conclusion, the results demonstrated that BMP15 is a major gene for litter size in Luzhong mutton sheep and candidate SNPs associated with litter size were identified.

Effects of FecB Mutation on Estrus, Ovulation, and Endocrine Characteristics in Small Tail Han Sheep

The Booroola fecundity gene (FecB) has a mutation that was found to increase the ovulation rate and litter size in Booroola Merino sheep. This mutation is also associated with the fecundity of small-tail han (STH) sheep, an important maternal breed used to produce hybrid offspring for mutton production in China. Previous research showed that the FecB gene affects reproduction in STH sheep, based on litter size records. However, the effects of this gene on estrus, ovulation, and endocrine characteristics in these sheep remain unclear. Here, we analyzed the traits mentioned earlier and compared them among the three FecB genotypes of STH ewes using estrus synchronization. Overall, 53 pluriparous ewes were selected from among 890 STH ewes and subjected to FecB genotyping for experiments to characterize estrous and ovulation rates. FecB heterozygous (+B) ewes presented an earlier onset of estrus (42.9 ± 2.2 h) and a shorter estrous cycle (17.2 ± 0.2 days) (P ≤ 0.05). The ovulation rates increased with the increasing copy number of the B allele (P ≤ 0.01). Ovulation time showed no significant differences among the three FecB genotypes. The serum concentrations of follicle-stimulating hormone (FSH), luteinizing hormone, estrogen (E2), and progesterone (P4) were measured in 19 of the ewes. Serum concentrations of E2 and FSH dramatically varied around the time of behavioral estrus. In FecB mutant homozygous (BB) ewes, E2 concentration had two peaks, which were higher (P ≤ 0.05) than those of ++ genotypes. FSH concentration of BB ewes was higher (P ≤ 0.05) than that of the ++ ewes just after estrus. The expression of the estrogen receptor 1 (ESR1) gene in the +B genotype was higher than in the other genotypes. Based on the data for the reproductive performance of STH ewes with the three FecB genotypes, our study suggests that the development of follicles in ewes with the B allele is dependent on the response to FSH regulated by E2 in the early stage. +B ewes, exhibiting moderate ovulation and litter size and a shorter estrous cycle, can be highly recommended in sheep crossbreeding systems for commercial mutton production. Moreover, this study provides useful information to conserve better and use the genetic resources of STH sheep in China.

Comparison of expression patterns of six canonical clock genes of follicular phase and luteal phase in Small-tailed Han sheep

The circadian rhythm is a biological rhythm that is closely related to the rhythmic expression of a series of clock genes. Results from several studies have indicated that clock genes are associated with the estrous cycle in female animals. Until now, the relationship between estrus cycle transition and clock gene expression in reproductive-axis-related tissues has remained unknown in Small-tailed Han (STH) sheep. This study was conducted to analyze the expression patterns of six canonical clock genes (Clock, BMAL1, Per1, Per2, Cry1, and Cry2) in the follicle phase and luteal phase of STH sheep. We found that all six genes were expressed in the brain, cerebellum, hypothalamus, pituitary, ovary, uterus, and oviduct in follicle and luteal phases. The results indicated that Clock expression was significantly higher in the cerebellum, hypothalamus, and uterus of the luteal phase than that of the follicle phase, whereas BMAL1 expression was significantly higher in the hypothalamus of the luteal phase than that of the follicle phase. Per1 expression was significantly higher in the brain, cerebellum, hypothalamus, and pituitary of the luteal phase than that of the follicle phase, and Per2 expression was significantly higher in the hypothalamus, pituitary, and uterus of the luteal phase than that of the follicle phase. Cry1 expression was significantly higher in the brain, cerebellum, and hypothalamus of the luteal phase than that of the follicle phase, whereas Cry2 expression was significantly higher in the pituitary of the luteal phase than that of the follicle phase. The clock gene expression in all tissues was different between follicle and luteal phases, but all clock gene mRNA levels were found to exhibit higher expression among seven tissues in the luteal phase. Our results suggest that estrous cycles may be associated with clock gene expression in the STH sheep. This is the first study to systematically analyze the expression patterns of clock genes of different estrous cycle in ewes, which could form a basis for further studies to develop the relationship between clock genes and the estrous cycle.

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.

Identification of Photoperiod-Induced LncRNAs and mRNAs in Pituitary Pars Tuberalis of Sheep

The pituitary pars tuberalis (PT) is the regulating center of seasonal reproduction, which can sense the melatonin signal and eventually cause downstream changes of GnRH secretion through TSHβ. Recently, lncRNAs have been identified in animal reproductive-related tissues, and they play important roles in reproductive regulation. Therefore, in this study, we expect to identify photoperiod-induced lncRNAs and genes in pituitary PT of sheep by comparison of expression profiles between short photoperiod (SP) and long photoperiod (LP). Through RNA-Seq, a total of 55,472 lncRNAs were identified in pituitary PT of Sunite ewes. The number of differentially expressed (DE) genes and lncRNAs between SP and LP increased gradually with the extension of LP (from LP7 to LP42). The notable LP-induced candidate genes included EYA3, TSHB, SIX1, DCT, VMO1, AREG, SUV39H2, and EZH2, and SP-induced genes involved ENSOARG00000012585, CHGA, FOS, SOCS3, and TH. In enriched pathways for DE genes and lncRNA target genes between SP and LP, the reproduction- and circadian-related pathways were highlighted. In addition, the interactome analysis of lncRNAs and their targets implied that MSTRG.209166 and its trans-target TSHB, MSTRG.288068 and its cis-target SIX1, and ENSOARG00000026131 and its cis-target TH might participate in regulation of seasonal reproduction. Together, these results will help to determine important photoperiod-induced lncRNAs and genes and give us some new insights into the epigenetic regulation of seasonal reproduction in sheep.

Genome-wide association study and inbreeding depression on body size traits in Qira black sheep (Ovis aries)

Qira black sheep is a famous indigenous sheep breed in China. The objectives of this study are to identify candidate genes related to body size, and to estimate the level of inbreeding depression on body size based on runs of homozygosity in Qira black sheep. Here, 188 adult Qira black sheep were genotyped with a high density (630 K) SNP chip and genome-wide association study for body weight and body size traits (including withers height, body slanting length, tail length, chest girth, chest width, and chest depth) were performed using an additive linear model. In consequence, 12 genome- and chromosome-wide significant SNPs and, accordingly, six candidate genes involved in muscle differentiation, metabolism and cell processes were identified. Of them, ZNF704 (zinc finger protein 704) was identified for body weight; AK2 (adenylate kinase 2) and PARK2 (parkin RBR E3 ubiquitin protein ligase) for tail length; MOCOS (molybdenum cofactor sulfurase) and ELP2 (elongator acetyltransferase complex subunit 2) for chest width; and MFAP1 (microfibril associated protein 1) for chest girth. Additionally, inbreeding depressions on body size were observed in the current herd. These results will provide insightful understandings into the genetic mechanisms of adult body size, and into the conservation and utilization of Qira black sheep.

Transcriptome Analysis Reveals Differentially Expressed Genes and Long Non-coding RNAs Associated With Fecundity in Sheep Hypothalamus With Different FecB Genotypes

Small-tailed Han sheep, with different FecB genotypes, manifest distinct ovulation rates and fecundities, which are due to differences in reproductive hormones secreted by the hypothalamic-pituitary-ovarian axis. Nevertheless, the function of the hypothalamus against a FecB mutant background on increasing ovulation rate is rarely reported. Therefore, we determined the expression profiles of hypothalamus tissue collected from six wild-type (WW) and six FecB mutant homozygous (BB) ewes at the follicular and luteal phases by whole-transcriptome sequencing. We identified 53 differentially expressed mRNAs (DEGs) and 40 differentially expressed long non-coding RNAs (DELs) between the two estrus states. Functional annotation analysis revealed that one of the DEGs, PRL, was particularly enriched in the hypothalamic function, hormone-related, and reproductive pathways. The lncRNA-target gene interaction networks and KEGG analysis in combination suggest that the lncRNAs LINC-676 and WNT3-AS cis-acting on DRD2 and WNT9B in different phases may induce gonadotropin-releasing hormone (GnRH) secretion. Furthermore, there were differences of regulatory elements and WNT gene family members involved in the follicular-luteal transition in the reproductive process between wild-type (WNT7A) and FecB mutant sheep (WNT9B). We combined the DEG and DEL data sets screened from different estrus states and genotypes. The overlap of these two sets was identified to select the mRNAs and lncRNAs that have major effects on ovulation. Among the overlapping molecules, seven DEGs and four DELs were involved in the follicular-luteal transition regulated by FecB mutation. Functional annotation analysis showed that two DEGs (FKBP5 and KITLG) were enriched in melanogenesis, oxytocin, and GnRH secretion. LINC-219386 and IGF2-AS were highly expressed in the BB ewes compared with WW ewes, modulating their target genes (DMXL2 and IGF2) to produce more GnRH during follicular development, which explains why mutated ewes produced more mature follicles. These results from expression profiling of the hypothalamus with the FecB mutation at different estrus states provide new insights into how the hypothalamus regulates ovulation under the effect of the FecB mutation.

Identification of genes associated with litter size combining genomic approaches in Luzhong mutton sheep

Litter size is one of the most important reproductive traits of sheep, which has pronounced effects on the profit of husbandry enterprises and enthusiasm of breeders. Despite the importance of litter size, the underlying genetic mechanisms have not been entirely elucidated. Therefore, based on a high-density SNP chip, genome-wide comparative analysis was performed between two groups with different fecundity to reveal candidate genes linked to litter size via detection of homozygosity and selection signatures in Luzhong mutton sheep. Consequently, nine promising genes were identified from six runs of homozygosity islands, and functionally linked to reproduction (ACTL7A, ACTL7B, and ELP1), embryonic development (KLF5 and PIBF1), and cell cycle (DACH1, BORA, DIS3, and MZT1). A total of 128 genes were observed under selection, of which HECW1 and HTR1E were related to total lambs born, GABRG3, LRP1B, and MACROD2 to teat number, and AGBL1 to reproductive seasonality. Additionally, the presence of inbreeding depression implies the urgency of reasonable mating system to increase litter size in the present herd. These findings provide a comprehensive insight to the genetic makeup of litter size, and also contribute to implementation of marker-assisted selection in sheep.

Expression analysis and single-nucleotide polymorphisms of SYNDIG1L and UNC13C genes associated with thoracic vertebral numbers in sheep (Ovis aries)

The objective of the current study was to analyze expression levels of synapse differentiation inducing 1-like (SYNDIG1L) and unc-13 homolog C (UNC13C) genes in different tissues, while single-nucleotide polymorphisms (SNPs) of two genes were associated with multiple thoracic vertebrae traits in both Small-tailed Han sheep (STH) and Sunite sheep (SNT). The expression levels of SYNDIG1L and UNC13C were analyzed in the brain, cerebellum, heart, liver, spleen, lung, kidney, adrenal gland, uterine horn, longissimus muscle, and abdominal adipose tissues of two sheep breeds with different thoracic vertebral number (TVN) sheep (T13 groups and T14 groups) by real-time quantitative polymerase chain reaction (RT-qPCR). Meanwhile, the polymorphisms of UNC13C gene g.52919279C>T and SYNDIG1L gene g.82573325C>A in T14 and T13 were genotyped by the Sequenom MassARRAY^® SNP assay, and association analysis was performed with the TVN. The results demonstrated that UNC13C gene was extensively expressed in 11 tissues. The expression of UNC13C gene in longissimus muscle of T14 groups of STH was significantly higher than that of T13 groups (P<0.05). SYNDIG1L gene was overexpressed in brain and cerebellum tissues, and the expression level of UNC13C gene in the brain and cerebellum of T13 groups in SNT was significantly higher than that of T14 groups (P<0.01). Association analysis showed that SNPs found in the UNC13C gene had no significant effects on TVN for both two genes. The polymorphism of SYNDIG1L g.82573325C>A was significantly correlated with the TVN in both STH (P<0.05) and SNT (P<0.01). Taken together, the SYNDIG1L gene was related to thoracic vertebral development, and this variation may be potentially used as a molecular marker to select the multiple thoracic vertebrae in sheep.

Photoperiod induced the pituitary differential regulation of lncRNAs and mRNAs related to reproduction in sheep

The pituitary is a vital endocrine organ that regulates animal seasonal reproduction by controlling the synthesis and secretion of the hormone. The change of photoperiod is the key factor affecting the function of the pituitary in animals, but the mechanism is unclear. Here, we studied the transcriptomic variation in pars distalis (PD) of the pituitary between short photoperiod (SP) and long photoperiod (LP) using RNA sequencing based on the OVX+E₂ sheep. 346 differentially expressed (DE) lncRNAs and 186 DE-mRNA were found in the PD. Moreover, function annotation analysis indicated that the reproductive hormones and photoperiod response-related pathways including aldosterone synthesis and secretion, insulin secretion, thyroid hormone synthesis, and circadian entrainment were enriched. The interaction analysis of mRNA-lncRNA suggested that MSTRG.240648, MSTRG.85500, MSTRG.32448, and MSTRG.304959 targeted CREB3L1 and DUSP6, which may be involved in the photoperiodic regulation of the PD. These findings provide resources for further study on the seasonal reproductive in ewes.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07536-y.

Litter Size of Sheep (Ovis aries): Inbreeding Depression and Homozygous Regions

Ovine litter size (LS) is an important trait showing variability within breeds. It remains largely unknown whether inbreeding depression on LS exists based on genomic homozygous regions, and whether the homozygous regions resulted from inbreeding are significantly associated with LS in sheep. We here reanalyze a set of single nucleotide polymorphism (SNP) chip of six breeds to characterize the patterns of runs of homozygosity (ROH), to evaluate inbreeding levels and inbreeding depressions on LS, and to identify candidate homozygous regions responsible for LS. Consequently, unique ROH patterns were observed among six sheep populations. Inbreeding depression on LS was only found in Hu sheep, where a significant reduction of 0.016, 0.02, and 0.02 per 1% elevated inbreeding FROH4-8, FROH>8 and the total inbreeding measure was observed, respectively. Nine significantly homozygous regions were found for LS in Hu sheep, where some promising genes for LS possibly via regulation of the development of oocytes (NGF, AKT1, and SYCP1), fertilization (SPAG17, MORC1, TDRD9, ZFYVE21, ADGRB3, and CKB), embryo implantation (PPP1R13B, INF2, and VANGL1) and development (DPPA2, DPPA4, CDCA4, CSDE1, and ADSSL1), and reproductive health (NRG3, BAG5, CKB, and XRCC3) were identified. These results from the present study would provide insights into the genetic management and complementary understandings of LS in sheep.

Genetic diversity and phylogenetic relationship of nine sheep populations based on microsatellite markers

The objective of this study was to assess the genetic diversity and phylogenetic relationship of nine sheep populations, including two famous high prolific populations and seven popular mutton populations raised in China. Overall, these sheep populations in this study exhibited a rich genetic diversity. Both the expected heterozygosity and Nei's unbiased gene diversity ranged from 0.64 to 0.75, with the lowest value found in Dorset sheep (DST) and the highest in Hu sheep (HUS) and Ba Han sheep (BAS). The polymorphic information content (PIC) varied between 0.59 in DST and 0.71 in HUS and BAS. Specifically, for individual breeds, the small-tail Han sheep (STH) and the four introduced populations did not display the expected diversity; therefore more attention should be paid to the maintenance of diversity during management of these populations. The results of un-weighted pair-group method (UPGMA) phylogenetic tree and structure analysis indicated that the nine investigated populations can be divided into two groups. Suffolk (SUF) and DST were clustered in one group, and the other group can be further divided into three clusters: German Mutton Merino (GMM)–BAS–Bamei Mutton sheep (BAM), HUS–STH and Du Han (DOS)–Dorper (DOP). This clustering result is consistent with sheep breeding history. TreeMix analysis also hinted at the possible gene flow from GMM to SUF. Together, an in-depth view of genetic diversity and genetic relationship will have important implications for breed-specific management.

Insight Into Pituitary lncRNA and mRNA at Two Estrous Stages in Small Tail Han Sheep With Different FecB Genotypes

The pituitary is a remarkably dynamic organ with roles in hormone (FSH and LH) synthesis and secretion. In animals with the FecB (fecundity Booroola) mutation, the pituitary experiences hormone fluctuations during the follicular-luteal transition, which is implicated in the expression and regulation of many genes and regulators. Long non-coding RNAs (lncRNAs) are a novel type of regulatory factors for the reproductive process. Nevertheless, the expression patterns of lncRNAs and their roles in FecB-mediated follicular development and ovulation remain obscure. Thus, we profiled the pituitary transcriptome during the follicular (F, 45 h after evacuation vaginal sponges) and luteal (L, 216 h after evacuation vaginal sponges) phases in FecB-mutant homozygous (BB) and wild-type (WW) Small Tail Han sheep. We identified 78 differentially expressed genes (DEGs) and 41 differentially expressed lncRNAs (DELs) between BB_F and BB_L, 32 DEGs and 26 DELs between BB_F and WW_F, 16 DEGs and 29 DELs between BB_L and WW_L, and 50 DEGs and 18 DELs between WW_F and WW_L. The results of real-time quantitative PCR (RT-qPCR) correlated well with the transcriptome data. In both the follicular and luteal phases, DEGs (GRID2, glutamate ionotropic receptor delta type subunit 2; ST14, ST14 transmembrane serine protease matriptase) were enriched in hormone synthesis, secretion, and action. MSTRG.47470 and MSTRG.101530 were the trans-regulated elements of ID1 (inhibitor of DNA binding 3, HLH protein) and the DEG ID3 (inhibitor of DNA binding 3, HLH protein), and EEF2 (eukaryotic translation elongation factor 2), respectively; these factors might be involved in melatonin and peptide hormone secretion. In the FecB-mediated follicular phase, MSTRG.125392 targeted seizure-related 6 homolog like (SEZ6L), and MSTRG.125394 and MSTRG.83276 targeted the DEG KCNQ3 (potassium voltage-gated channel subfamily Q member 3) in cis, while MSTRG.55861 targeted FKBP4 (FKBP prolyl isomerase 4) in trans. In the FecB-mediated luteal phase, LOC105613905, MSTRG.81536, and MSTRG.150434 modulated TGFB1, SMAD3, OXT, respectively, in trans. We postulated that the FecB mutation in pituitary tissue elevated the expression of certain genes associated with pituitary development and hormone secretion. Furthermore, this study provides new insights into how the pituitary regulates follicular development and ovulation, illustrated by the effect of the FecB mutation.

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.

Expression analysis of DIO2, EYA3, KISS1 and GPR54 genes in year-round estrous and seasonally estrous rams

The expression characteristics of the hypothalamic-pituitary-gonadal (HPG) axis-related candidate genes, DIO2, EYA3, KISS1 and GPR54, were analyzed in year-round estrous rams (small-tail Han sheep, STH) and seasonally estrous rams (Sunite sheep, SNT) using qPCR. The results were as follows: DIO2 was mainly expressed in pituitary, and KISS1 was specifically expressed in hypothalamus in the two groups. However, EYA3 and GPR54 were widely expressed in the cerebrum, cerebellum, hypothalamus, pituitary, testis, epididymis, vas deferens and adrenal gland tissues in both breeds, with significant differences in the cerebellum, hypothalamus, pituitary, testis and vas deferens tissues. We speculated that DIO2 and KISS1 may have positive roles in different regions in ram year-round estrus. Moreover, the expression patterns of EYA3 and GPR54 suggested that they may regulate the estrous mode of ram via testis and vas deferens. This is the first study to systematically analyze the expression patterns of HPG axis-related genes in rams.

Luzhong mutton sheep: inbreeding and selection signatures

Intense artificial selection has been imposed to Luzhong mutton sheep population in the past years. Improvements on growth and reproductive performance are two breeding goals in the present herd. Although some progresses were phenotypically observed possibly due to inbreeding induced by strong selection in terms of these traits, the genomic evaluation was poorly understood. Therefore, a high-density SNP array was used to characterize the pattern of runs of homozygosity (ROH), estimate inbreeding and inbreeding depressions on early growth performance and litter size based upon ROH, and scan positive selection signatures of recent population. Consequently, a low inbreeding level was observed which had negative effects on litter size, but not on early growth performance. And 160 genes were under selection, of which some were reported to be linked to several traits of sheep including body weight, litter size, carcass and meat quality, milk yield and composition, fiber quality and health, and the top genes were associated with growth (growth hormone [GH]- growth hormone receptor [GHR]- Insulin-like growth factor 1 [IGF1] axis) and litter size (bone morphogenic proteins [BMPs]-associated). The effectiveness of previous breeding measures was highlighted, but purging selection was proposed to alleviate the inbreeding depression on litter size, providing some genomic insights to breeding management of Luzhong mutton sheep.

Identification and Characterization of Hypothalamic Alternative Splicing Events and Variants in Ovine Fecundity-Related Genes

Simple Summary

Previous studies revealed that alternative splicing (AS) events and gene variants played key roles in reproduction. However, their location and distribution in hypothalamic fecundity-related genes in sheep without the FecB mutation remain largely unknown. In this study, we performed a correlation analysis of transcriptomics and proteomics, and the results suggested several differentially expressed genes (DEGs)/differentially expressed proteins (DEPs), including galectin 3 (LGALS3), aspartoacylase (ASPA) and transthyretin (TTR), could be candidate genes influencing ovine litter size. Further analysis suggested that AS events, single nucleotide polymorphisms (SNPs) and microRNA (miRNA)-binding sites existed in key DEGs/DEPs, such as ASPA and TTR. This study provides a new insight into ovine and even other mammalian reproduction.

Abstract

Previous studies revealed that alternative splicing (AS) events and gene variants played key roles in reproduction; however, their location and distribution in hypothalamic fecundity-related genes in sheep without the FecB mutation remain largely unknown. Therefore, in this study, we described the hypothalamic AS events and variants in differentially expressed genes (DEGs) in Small Tail Han sheep without the FecB mutation at polytocous sheep in the follicular phase vs. monotocous sheep in the follicular phase (PF vs. MF) and polytocous sheep in the luteal phase vs. monotocous sheep in the luteal phase (PL vs. ML) via an RNA-seq study for the first time. We found 39 DEGs with AS events (AS DEGs) in PF vs. MF, while 42 AS DEGs were identified in PL vs. ML. No DEGs with single nucleotide polymorphisms (SNPs) were observed in PF vs. MF, but five were identified in PL vs. ML. We also performed a correlation analysis of transcriptomics and proteomics, and the results suggested several key DEGs/differentially expressed proteins (DEPs), such as galectin 3 (LGALS3) in PF vs. MF and aspartoacylase (ASPA) and transthyretin (TTR) in PL vs. ML, could be candidate genes influencing ovine litter size. In addition, further analyses suggested that AS events, SNPs and miRNA-binding sites existed in key DEGs/DEPs, such as ASPA and TTR. All in all, this study provides a new insight into ovine and even other mammalian reproduction.