Mutations in BMPR-IB and BMP-15 genes are associated with litter size in Small Tailed Han sheep (Ovis aries)

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.

Multiomics Analyses Provide New Insight into Genetic Variation of Reproductive Adaptability in Tibetan Sheep

Domestication and artificial selection during production-oriented breeding have greatly shaped the level of genomic variability in sheep. However, the genetic variation associated with increased reproduction remains elusive. Here, two groups of samples from consecutively monotocous and polytocous sheep were collected for genome-wide association, transcriptomic, proteomic, and metabolomic analyses to explore the genetic variation in fecundity in Tibetan sheep. Genome-wide association study revealed strong associations between BMPR1B (p.Q249R) and litter size, as well as between PAPPA and lambing interval; these findings were validated in 1,130 individuals. Furthermore, we constructed the first single-cell atlas of Tibetan sheep ovary tissues and identified a specific mural granulosa cell subtype with PAPPA-specific expression and differential expression of BMPR1B between the two groups. Bulk RNA-seq indicated that BMPR1B and PAPPA expressions were similar between the two groups of sheep. 3D protein structure prediction and coimmunoprecipitation analysis indicated that mutation and mutually exclusive exons of BMPR1B are the main mechanisms for prolific Tibetan sheep. We propose that PAPPA is a key gene for stimulating ovarian follicular growth and development, and steroidogenesis. Our work reveals the genetic variation in reproductive performance in Tibetan sheep, providing insights and valuable genetic resources for the discovery of genes and regulatory mechanisms that improve reproductive success.

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.

LncRNA12097.1 contributes to endometrial cell growth by enhancing YES1 activating β-catenin via sponging miR-145-5p

Despite previous investigations elucidating the regulatory mechanisms of long non-coding RNAs (lncRNAs) in endometrial function and reproductive disorders, the precise pathways through which lncRNAs impact endometrial functions and fertility remain unclear. In this study, we performed an expression profile analysis of lncRNAs in the endometrial tissue of Hu sheep with different prolificacy, identifying 13,707 lncRNAs. We discovered a bidirectional lncRNA, designated lncRNA12097.1, exhibiting significant up-regulation exclusively in the endometrium of Hu sheep with high fecundity. Functional analyses revealed lncRNA12097.1 significantly enhanced proliferation and cell cycle progression in both endometrial epithelial cell (EEC) and stromal cells (ESC), while inhibiting apoptosis in these cell types. Mechanistically, we demonstrated a directly interaction between lncRNA12097.1 and miR-145-5p, with YES proto-oncogene 1 (YES1) being identified as a validated target of miR-145-5p. Interference with lncRNA12097.1 resulted in suppressed cell growth through down-regulation of YES1 expression, which could be rescued by miR-145-5p. Furthermore, lncRNA12097.1 functions as a competitive endogenous RNA (ceRNA) for miR-145-5p in ESCs, sequestering miR-145-5p and preventing its binding to the 3'UTR of YES1 mRNA. This interaction led to increased expression of YES1 and subsequent activation of downstream β-catenin signaling, thereby promoting ESC growth in Hu sheep. These findings provide novel molecular insights into the mechanisms underlying prolificacy in sheep.

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.

The Significance of a Multilocus Analysis for Assessing the Biodiversity of the Romanov Sheep Breed in a Comparative Aspect

The Romanov breed was evaluated using immunological and genetic markers. The seven blood group systems were characterized with a greater accuracy than in previous works on sheep in the Russian Federation, and were compared to eight ruminant species. Unlike other breeds, Romanov sheep shows a higher frequency of HB^A than HB^B alleles. There are 3-4 genotypes at the transferrin locus whereas in other breeds 6-11 genotypes have been found. At the albumin locus, the majority of the identified genotypes were heterozygotes, unlike in the other breeds studied. In the prealbumin locus, the Romanov breed was the only one where all the genotypes were heterozygous. We speculate that polymorphism at two loci (BMP-15 and BMPR1B) could effect on the high ovulation rates of Romanov sheep. Based on different genetic markers, the prevalence of heterozygotes in the Romanov sheep could determine their higher viability. A cluster analysis showed the close proximity of 12 populations of the Romanov breed, as the breeding stock come from the Yaroslavl region.

Effects of New Mutations in BMPRIB, GDF9, BMP15, LEPR, and B4GALNT2 Genes on Litter Size in Sheep

Prolificacy is a crucial characteristic of livestock, particularly for species such as sheep that have many births. The objectives of this study were as follows: (1) to investigate the genetic diversity of the 13 new and 7 known variants in the BMPRIB, GDF9, BMP15, LEPR, and B4GALNT2 genes in Ujimqin (UM), the F1 population of Dorper × Ujimqin crossbred (DPU), the F1 population of Suffolk × Ujimqin crossbred (SFKU), Sonid sheep (SN), Tan sheep (Tan), Hu sheep (Hu), and Small-tailed Han sheep (STH) sheep breeds/populations; (2) to perform an association analysis of the above 20 variants with litter size in 325 UM, 304 DPU, and 66 SFKU sheep populations; (3) to compare the frequencies of the litter-size-related alleles of these 20 variants among 8 sheep breeds/populations (the above seven sheep breeds + Mongolia sheep breed). With the use of the Sequenom MassARRAY®SNP assay technology, these 20 mutations were genotyped. The association analysis results showed that the c.746A>G (FecB) mutation in BMPR1B was significantly associated with the litter size of UM and DPU, the c.994A>G (FecGA) in GDF9 was significantly associated with the litter size of SFKU, and the c.31_33CTTinsdel (B1) in BMP15 was significantly associated with the litter size of UM. Our findings might provide valuable genetic markers for expanding sheep litter sizes.

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.

LncRNA-412.25 activates the LIF/STAT3 signaling pathway in ovarian granulosa cells of Hu sheep by sponging miR-346

Although long non-coding RNAs (lncRNAs) are reported to regulate follicular development and reproductive disease pathogenesis, the underlying mechanisms have not been elucidated. In this study, lncRNA expression profiling of different-sized healthy follicles from Hu sheep with different prolificacy revealed 50 613 lncRNAs. Numerous lncRNAs were differentially expressed among different comparison groups. This study characterized one novel transcript, lncRNA-412.25 (from healthy follicles with a diameter of >5 mm), which was predominantly expressed in the high prolificacy group and localized to the cytoplasm of granulosa cells (GCs). LncRNA-412.25 knockdown promoted and inhibited Hu sheep GC apoptosis and proliferation, respectively. Interestingly, lncRNA-412.25 could directly bind to miR-346, which can target the gene of leukemia inhibitory factor (LIF). Knockdown of lncRNA-412.25 promoted GC apoptosis by downregulating LIF expression, where this effect was attenuated by miR-346. Moreover, the miR-346 inhibitor mitigated the lncRNA-412.25 knockdown-induced downregulation of phosphorylated protein of signal transducer and activator of transcription 3 (STAT3), which was validated using immunofluorescence analysis. Our results demonstrated that lncRNA-412.25 regulates GC proliferation and apoptosis in Hu sheep by binding to miR-346 and then activating the LIF/STAT3 pathway. These findings provide novel insights into the mechanisms underlying prolificacy in sheep.

Molecular Characterization, Expression Profiles of SMAD4, SMAD5 and SMAD7 Genes and Lack of Association with Litter Size in Tibetan Sheep

SMAD4, SMAD5 and SMAD7 belonging to the transforming growth factor β (TGF-β) superfamily are indispensable for oocyte formation and development, ovarian organogenesis and folliculogenesis. However, only a few studies have investigated the characteristics of SMAD4, SMAD5 and SMAD7 in Tibetan sheep and the effect of their polymorphism on litter size. In this study, we examined the expression of SMAD4, SMAD5 and SMAD7 in 13 tissues of Tibetan sheep by reverse transcription-quantitative polymerase chain reaction. Further, cDNA of these genes was cloned, sequenced and subjected to bioinformatics analysis. DNA sequencing was also used to detect single nucleotide polymorphisms (SNPs). However, iM-LDRTM technology was used for SNP genotyping. Associations between polymorphisms and litter size were analyzed using data from genotyping of 433 Tibetan sheep. The results showed that the expression of SMAD4, SMAD5 and SMAD7 genes was ubiquitous in the tissues of Tibetan sheep, such as the ovary, uterus and oviduct, hypothalamus, hypophysis, heart, liver, spleen, lung, kidney, rumen, duodenum and longissimus dorsi. However, the expression was unbalanced and upregulated in the spleen, lung, ovary and uterus and downregulated in the longissimus dorsi. The bioinformatics analysis showed that SMAD4, SMAD5 and SMAD7 in Tibetan sheep encoded proteins of 533, 465 and 427 amino acids, respectively. Sequence homology analysis of the three proteins among other animals showed that the sequences of SMAD4, SMAD5 and SMAD7 of Tibetan sheep were similar to those in sheep, yak, cattle, dog, human, pig, chimpanzee, rhesus monkey and house mouse. Two synonymous mutations, g.51537A>G and g.319C>T, were detected in SMAD5 and SMAD7, respectively. The associations of these SNPs and litter size were determined, and it was found that both g.51537A>G and g.319C>T have no significant effect on the litter size of Tibetan sheep. The results provided novel insights into the molecular characterization, expression profiles and polymorphisms of SMAD4, SMAD5 and SMAD7 in Tibetan sheep, but our results do not support associations between these genes and the litter size of Tibetan sheep.

Expression Profiles of the Ovine IL18 Gene and Association of Its Polymorphism With Hematologic Parameters in Hu Lambs

Hematological traits are important indexes to evaluate health status and immunological conditions in human and livestock. In this study, we measured the hematologic indexes of 819 male Hu lambs and performed the descriptive statistical analysis. The results showed the coefficients of variation of partial indexes were >10%, and the heritability for mean erythrocyte volume (MCV), white blood cell count (WBC), hemoglobin concentration (HGB), hematocrit (HCT), and red blood cell (RBC) distribution-standard deviation (RDW_SD) were medium to high, ranging from 0.17 to 0.43. In addition, Interleukin 18 (IL18), as an important regulator of both innate and acquired immune responses, was selected as candidate gene and subjected to the expression profile analysis, single nucleotide polymorphism (SNP) scanning and association analysis by using quantitative real-time PCR (qRT-PCR), PCR amplification, Sanger sequencing, and KASP genotyping. The results of qRT-PCR indicated that IL18 is predominantly expressed in lymph and lung compared with that in the other tested tissues. In addition, three novel polymorphisms (g. 24991544 A > G, g. 24991651 A > G, and g. 24991749 C > T) were identified in IL18, and the three SNPs were in a strong linkage state. Therefore, only a SNP was genotyped and performed association analysis in the enlarged experimental population, the result of association analysis demonstrated that the polymorphism g. 24991651 A > G was significantly associated with RBC, MCV, MCHC, and RDW_CV. These results will provide the reference values and the novel genetic markers of hematological parameters in sheep.

Integrative Genome-Wide DNA Methylome and Transcriptome Analysis of Ovaries from Hu Sheep with High and Low Prolific

DNA methylation plays an important role in biological processes by affecting gene expression. However, how DNA methylation regulates phenotypic variation in Hu sheep remains unclear. Therefore, we generated genome-wide DNA methylation and transcriptomic profiles in the ovaries of Hu sheep with different prolificacies and genotypes (FecBB and FecB+). Results showed that ovary DNA methylome and transcriptome were significantly different between high prolificacy and low prolificacy Hu sheep. Comparative methylome analyses identified 10,644, 9,594, and 12,214 differentially methylated regions and 87, 1,121, and 2,375 genes, respectively, showing differential expression levels in three different comparison groups. Female reproduction-associated differentially methylated regions-related genes and differentially expressed genes were enriched, thereby the respective interaction networks were constructed. Furthermore, systematical integrative analyses revealed a negative correlation between DNA methylation around the transcriptional start site and gene expression levels, which was confirmed by testing the expression of integrin β2 subunit (ITGB2) and lysosome-associated protein transmembrane-4 beta (LAPTM4B) in vivo and in vitro. These findings demonstrated that DNA methylation influences the propensity for prolificacy by affecting gene expression in the ovaries, which may contribute to a greater understanding of the epigenome and transcriptome that will be useful for animal breeding.

An ancient positively selected BMPRIB missense variant increases litter size of Mongolian sheep populations following latitudinal gradient

New gene mutation origination is a driving force for the evolution of organisms. The effect of FecB mutation in BMPRIB gene on the litter size of sheep has been well known for a long time, each copy of the mutant allele increases litter size by 0.4-0.5. However, the origin and adaptive evolution mechanism of FecB mutation are still unclear. Here we carried on the thorough analysis on evolutionary features of BMPRIB gene and found that 150 species as a whole is under purifying selection while sheep lineage shows evidence of positive selection. The results of allele age estimation revealed that the FecB mutation in Mongolian sheep of China originated in Mongolian Plateau at about 5000 years ago. Due the two shape drops in temperature subsequently, Mongolian sheep migrated from north to south following the northern nomadic people. Accordingly, the FecB mutant allele frequency increased, with the lowest in sheep locating at Mongolian plateau (0.01) and the highest in sheep locating at Yangtze River valley (0.96). In conclusion, the FecB mutation in Mongolian sheep of China originated in Mongolian Plateau at about 5000 years ago, and the differentiated litter size of Mongolian sheep might be the result of adaptation to various environments during the migration following latitudinal gradient. This study may well exemplify selection on an ancient variation triggered by drastic ecological shifts, and is also helpful to analyze the adaptive evolution mechanism of economic traits of domestic animals and identify major genes and molecular markers.

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.

Novel InDel variations of the Cry2 gene are associated with litter size in Australian White sheep

Cryptochrome 2 (Cry2) gene regulates circadian rhythm and affects reproduction and pregnancy. Therefore, this study aimed to explore polymorphisms of the Cry2 gene and their associations with litter size at different parity in Australian White (AuW) ewes. Five putative insertion or deletion mutations within the Cry2 gene were selected to study their association with litter size. Two novel deletion mutations were identified in intronic region of Cry2 gene and were genotyped by agarose gel electrophoresis and DNA sequencing. The polymorphism information content (PIC) indicated that both mutations were low polymorphism in tested groups. Statistical analysis revealed that the P1-Del-6-bp was significantly correlated with litter size at third parity (P = 0.010), in which individuals with insertion/deletion (ID) genotype had larger litter size than insertion/insertion (II) genotype (P < 0.05). Whereas, the P2-Del-6-bp was significantly correlated with litter size at first parity (P = 0.036), in which individuals with insertion/insertion (II) genotype had larger litter size than insertion/deletion (ID) genotype (P < 0.05). Collectively, these findings may provide new insights to expedite molecular breeding in sheep through marker-assisted selection strategies (MAS).

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.

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.

Detection of mRNA Expression and Copy Number Variations Within the Goat Fec B Gene Associated With Litter Size

The Booroola fecundity (Fec ^B ) gene, as the first major fecundity gene identified in Booroola sheep, has attracted careful attention. So far, previous research have uncovered the Fec^B mutation (Q249R) as the main mutation by virtue of which sheep exhibits multiple lambing phenomena. This mutation is now being intensively studied and widely used. However, such effect of the Fec^B mutation has not been applied to goats, and similar types of the Fec ^B gene in goats still need to be studied. Thus, the current study attempted to verify potential mutations in the goat Fec ^B gene as well as investigate their functions related to fecundity. First, Fec ^B expression was investigated in six different goat tissues, and we found that Fec ^B expression was highest in the mammary gland, followed by the ovary. Next, the influence of the Fec ^B gene was analyzed from a new perspective, where five potential copy number variations (CNVs) (CNV1-5) within the Fec ^B gene were identified for the first time, and then their effects on litter size were measured. Our results point out that CNV3 (P = 3.44E-4) and CNV5 (P = 0.034) could significantly influence the litter size of goats. Identically, the combination genotype of CNV3 and CNV5 which consisted of their dominant genotypes was also significantly associated with goat litter size (P = 7.80E-5). Hence, CNV3 and CNV5 could serve as potential DNA molecular markers applied to DNA editing and DNA microarray. Additionally, the abovementioned study has laid a theoretical foundation for the detection of potential fertility-related quantitative trait loci within the goat Fec ^B gene.

Survey of the relationship between polymorphisms within the BMPR1B gene and sheep reproductive traits

The BMPRIB gene is one of the main genes that can be used as a molecular genetic marker for the early selection of highly productive ewes. It is well-documented that the p.Q249R (g.746A > G) is the first mutation in the kinase domain of the BMPR1B gene that is highly related to increased ovulation rate and litter size. It is likely that the presence of the p.Q249R mutation in the sheep population is one of the factors contributing to the outstanding productivity of the sheep. Moreover, in recent years, researchers have been explored other polymorphisms in the BMPR1B gene with respect to reproductive traits in sheep. Therefore, we carried out the current study to evaluate the association between polymorphisms in this gene and sheep litter size from all appropriate studies. As a result, among 41 polymorphisms in the ovine BMPRIB gene, eight variants, including p.Q249R (g.746A > G), g.29362047T > C, g.29427689G > A, BMPR1B-2 (ss:1960972599), g.29382337G > A, g.29382340G > A, rs1092293287 (10 bp insertion/deletion) and g.29380965A > G were found to be associated with litter size in sheep. This systematic analysis presents the most current data evidence for BMPRIB polymorphisms, highlighting the need for further large-scale studies to determine more important variants.

Polymorphisms in BMPRIB gene affect litter size in Chinese indigenous sheep breed

The BMPRIB gene belongs to the TGF-β superfamily and is considered to be a regulator of sheep reproductive performance. Single nucleotide polymorphisms (SNPs) of BMPRIB gene in the Small Tail Han, Hu, Mongolian, Oula, Gansu Alpine Fine-wool, Dorper and Australian White sheep were detected by Sanger sequencing. Five SNPs (rs427897187 G > A, rs418841713 A > G, rs159952533 T > C, rs429416173 C > A and rs403555643 A > G) of BMPRIB gene were identified. For rs427897187 G > A, further analysis revealed that genotype GG and GA had 0.26 (p < 0.05) and 0.33 (p < 0.05) litter size less than those with genotype AA in Oula sheep. For rs403555643 A > G, further analysis revealed that genotype GG and AG had 0.65 (p < 0.05) and 0.38 (p < 0.05) litter size more than those with genotype AA in Oula sheep, and genotype GG had 0.56 (p < 0.05) litter size more than those with genotype AA in Mongolian sheep. The results showed that rs427897187 G > A and rs403555643 A > G are potential molecular markers wich could improve litter size of Chinese indigenous sheep and be used in Chinese indigenous sheep breeding.

Association between novel variants in BMPR1B gene and litter size in Mongolia and Ujimqin sheep breeds

Prolificacy is an important trait of animals, specifically for sheep. The Bone morphogenetic protein receptor 1B (BMPR1B) is a major gene affecting the litter size of many sheep breeds. The well-known FecB mutation (Q249R) was associated fully with the hyper prolific phenotype of Booroola Merino. However, the identification of variation in all exonic regions of BMPR1B was rare. In this study, we sequenced all exonic regions of BMPR1B gene of Mongolia sheep breed, and ten novel variants were detected by direct sequencing. Among them, the litter size of the Mongolia ewes with the CC genotype was significantly higher (0.34 additional lambs, p < .05) than those with the TT genotype of the g.29346567C>T single nucleotide polymorphism (SNP). The litter size of the Mongolia ewes with the TT genotype was significantly higher (0.19 additional lambs, p < .05 and .31 additional lambs, p < .01, respectively) than those with the GT and GG genotypes of the c.1470G>T SNP. The silent c.1470G>T mutation is predicted to increase the stability of the mRNA secondary structure through reducing minimum free energy and is predicted to change the mRNA secondary structure of BMPR1B. Our findings may give potentially useful genetic markers for increasing litter size in sheep.

Impact of Booroola fecundity gene introgression on sheep production: Indian perspective

Surge of demand for meat from diverse sources in a time of decling land and other natural resources shifted the trend of producing more sheep per sheep. Prolificacy aspect in sheep can be used as an important arsenal to fulfil this objective. FecB gene mutation perhaps spread from Garole microsheep of India, has been widely used for introduction of prolificacy in monotocus sheep worldwide. The ovulation rate and litter size are two most important reproductive parameters influenced by FecB gene besides small influences on other sheep production parameters. While implementing the FecB introgression programmes precaution needs to be observed to avoid lamb mortality. The results of the experiments conducted worldwide on FecB introgression programmes suggest that high prolificacy is economically beneficial under intensive and semi-intensive management system of sheep rearing. There were variable effects of FecB genotype documented from different parts of the world. With goal oriented research efforts, the prolificacy in FecB carrier Avishaan sheep, a triple breed cross developed from Garole, Malpura and Patanwadi at ICAR-CSWRI are showing promising sheep for future. There are several challenges for using FecB gene in the future using assisted reproductive and genome editing technology. The effect of FecB gene on different sheep production parameters as reported in the literature are summarized in present work.

The expression and mutation of BMPR1B and its association with litter size in small-tail Han sheep (Ovis aries)

Previous studies have shown that BMPR1B promotes follicular development and ovarian granulosa cell proliferation, thereby affecting ovulation in mammals. In this study, the expression and polymorphism of the BMPR1B gene associated with litter size in small-tail Han (STH) sheep were determined. The expression of BMPR1B was detected in 14 tissues of STH sheep during the follicular phase as well as in the hypothalamic–pituitary–gonadal (HPG) axis of monotocous and polytocous STH sheep during the follicular and luteal phases using quantitative polymerase chain reaction (qPCR). Sequenom MassARRAY^® single nucleotide polymorphism (SNP) technology was also used to detect the polymorphism of SNPs in seven sheep breeds. Here, BMPR1B was highly expressed in hypothalamus, ovary, uterus, and oviduct tissue during the follicular phase, and BMPR1B was expressed significantly more in the hypothalamus of polytocous ewes than in monotocous ewes during both the follicular and luteal phases (P<0.05). For genotyping, we found that genotype and allele frequencies of three loci of the BMPR1B gene were extremely significantly different (P<0.01) between the monotocous and polytocous groups. Association analysis results showed that the g.29380965A>G locus had significant negative effects on the litter size of STH sheep, and the combination of g.29380965A>G and FecB (Fec – fecundity and B – Booroola; A746G) at the BMPR1B gene showed that the litter size of AG–GG, AA–GG, and GG–GG genotypes was significantly higher compared with other genotypes (P<0.05). This is the first study to find a new molecular marker affecting litter size and to systematically analyze the expression of BMPR1B in different fecundity and physiological periods of STH 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.

Polymorphism of fecundity genes (BMP15 and GDF9) and their association with litter size in Bangladeshi prolific Black Bengal goat

Goat farming in Bangladesh is primarily centred on indigenous Black Bengal goat, a highly prolific breed. Searching for genetic markers associated with prolificacy in this breed is vital for the country's goat breeding industry. However, there are no reports on polymorphisms associated with the fertility of Bangladeshi Black Bengal goats. This study investigated two major fecundity genes-bone morphogenetic protein 15 (BMP15) and growth differentiation factor 9 (GDF9) to detect any possible mutations in these two genes associated with litter size in Black Bengal goats. Blood samples were collected from 40 raised goats in Hathazari Government Goat Farm, Bangladesh. Genomic DNA was extracted; PCR amplification was performed; and sequencing of PCR products was performed to detect polymorphism loci in the target genes. Five SNPs viz. C735A, C743A, G754T, C781A and C808G were detected in exon 2 of BMP15 gene. A SNP (T1173A) was detected in GDF9 exon 2. Association results show that SNPs at the 735, 754 and 781 nucleotide positions of BMP15 exon 2 had a significant association with litter size in Black Bengal goat. The effect of parity was also highly significant (P < 0.001) on litter size. For the first time, this study explored SNP loci in fecundity genes in Bangladeshi prolific Black Bengal goats. Further studies with many genetically unrelated animals for assessing the association of these loci and others in the fecundity genes with litter size may be useful.

lncRNA FDNCR promotes apoptosis of granulosa cells by targeting the miR-543-3p/DCN/TGF-β signaling pathway in Hu sheep

Long non-coding RNAs (lncRNAs) regulate the development of follicles and reproductive diseases, but the mechanisms by which lncRNAs regulate ovarian functions and fertility remain elusive. We profiled the expression of lncRNAs in ovarian tissues of Hu sheep with different prolificacy and identified 21,327 lncRNAs. Many of the lncRNAs were differentially expressed in different groups. We further characterized an lncRNA that was predominantly expressed in the ovaries of the low prolificacy Fec^B+ (LPB+) group and mainly present in granulosa cells (GCs), and the expression of this lncRNA decreased during follicular development, which we named follicular development-associated lncRNA (FDNCR). Next, we found that FDNCR directly binds miR-543-3p, and decorin (DCN) was identified as a target of miR-543-3p. FDNCR overexpression promoted GC apoptosis through increased expression of DCN, which could be attenuated by miR-543-3p. Furthermore, miR-543-3p increased and FDNCR reduced the expression of transforming growth factor-β (TGF-β) pathway-related genes, including TGF-β1 and inhibin beta A (INHBA), which were upregulated upon DCN silencing. Our results demonstrated that FDNCR sponges miR-543-3p in GCs and prevents miR-543-3p from binding to the DCN 3′ UTR, resulting in DCN transactivation and TGF-β pathway inhibition and promotion of GC apoptosis in Hu sheep. These findings provide insights into the mechanisms underlying prolificacy in sheep.

Polymorphism Detection of GDF9 Gene and Its Association with Litter Size in Luzhong Mutton Sheep (Ovis aries)

Simple Summary

GDF9 and BMPR1B are two important reproduction genes. In this study, the whole coding region of GDF9 was sequenced, of which the mutations were detected in Luzhong mutton sheep. The results suggested that two single nucleotide polymorphisms (SNPs), g.41768501A > G and g.41768485 G > A in GDF9 gene were associated with litter size. The g.41768485 G > A is a missense mutation which is predicted to affect the tertiary structure of the protein. Thus, these two mutations may be potential effective genetic markers to improve the litter size in sheep.

Abstract

Litter size is one of the most important economic traits in sheep. GDF9 and BMPR1B are major genes affecting the litter size of sheep. In this study, the whole coding region of GDF9 was sequenced and all the SNPs (single nucleotide polymorphisms) were determined in Luzhong mutton ewes. The FecB mutation was genotyped using the Sequenom MassARRAY^®SNP assay technology. Then, the association analyses between polymorphic loci of GDF9 gene, FecB, and litter size were performed using a general linear model procedure. The results showed that eight SNPs were detected in GDF9 of Luzhong mutton sheep, including one novel mutation (g.41769606 T > G). The g.41768501A > G, g.41768485 G > A in GDF9 and FecB were significantly associated with litter size in Luzhong mutton ewes. The g.41768485 G > A is a missense mutation in the mature GDF9 protein region and is predicted to affect the tertiary structure of the protein. The results preliminarily demonstrated that GDF9 was a major gene affecting the fecundity of Luzhong mutton sheep and the two loci g.41768501A > G and g.41768485 G > A may be potential genetic markers for improving litter size.

Pineal gland transcriptomic profiling reveals the differential regulation of lncRNA and mRNA related to prolificacy in STH sheep with two FecB genotypes

BACKGROUND

Long noncoding RNA (lncRNA) has been identified as important regulator in hypothalamic-pituitary-ovarian axis associated with sheep prolificacy. However, little is known of their expression pattern and potential roles in the pineal gland of sheep. Herein, RNA-Seq was used to detect transcriptome expression pattern in pineal gland between follicular phase (FP) and luteal phase (LP) in FecB^BB (MM) and FecB⁺⁺ (ww) STH sheep, respectively, and differentially expressed (DE) lncRNAs and mRNAs associated with reproduction were identified.

RESULTS

Overall, 135 DE lncRNAs and 1360 DE mRNAs in pineal gland between MM and ww sheep were screened. Wherein, 39 DE lncRNAs and 764 DE mRNAs were identified (FP vs LP) in MM sheep, 96 DE lncRNAs and 596 DE mRNAs were identified (FP vs LP) in ww sheep. Moreover, GO and KEGG enrichment analysis indicated that the targets of DE lncRNAs and DE mRNAs were annotated to multiple biological processes such as phototransduction, circadian rhythm, melanogenesis, GSH metabolism and steroid biosynthesis, which directly or indirectly participate in hormone activities to affect sheep reproductive performance. Additionally, co-expression of lncRNAs-mRNAs and the network construction were performed based on correlation analysis, DE lncRNAs can modulate target genes involved in related pathways to affect sheep fecundity. Specifically, XLOC_466330, XLOC_532771, XLOC_028449 targeting RRM2B and GSTK1, XLOC_391199 targeting STMN1, XLOC_503926 targeting RAG2, XLOC_187711 targeting DLG4 were included.

CONCLUSION

All of these differential lncRNAs and mRNAs expression profiles in pineal gland provide a novel resource for elucidating regulatory mechanism underlying STH sheep prolificacy.

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.

Study on the correlation between BMPR1B protein in sheep blood and reproductive performance

Sheep reproductive performance is one of the important economic traits in sheep farming. The bone morphogenetic protein receptor 1B (BMPR1B) gene and protein may play an important role in sheep fertility. This study was to investigate the association of blood BMPR1B protein expression with reproductive performance in sheep. Mongolian sheep with single and twin births and polytocous Small Tail Han sheep were selected due to differences in birth numbers. The BMPR1B mRNA in sheep blood was measured by a reverse transcription-polymerase chain reaction as well as the BMPR1B protein was measured by enzyme-linked immunosorbent assay in blood samples of Mongolian and Small Tail Han sheep. The results demonstrated that blood BMPR1B concentration in Mongolian sheep with twin birth was higher (P < 0.05) than Small Tail Han sheep and Mongolian sheep with single birth. The protein concentration in the anestrus season was higher (P < 0.045) than those in the estrus season for both Mongolian and Small Tail Han sheep. Moreover, BMPR1B concentration in Mongolian sheep increased (P < 0.05) at the age of 6 to 12 mo and that in Small Tail Han sheep increased (P < 0.05) at the age of 3 to 6 mo. The result indicates that the increase in BMPR1B protein concentrations in the blood of Mongolian ewes and Small Tail Han ewes may be beneficial to follicular development, but too high or too low of this blood protein concentration in Mongolian and Small Tail Han sheep is not conducive to ovulation.

Investigating the Polymorphism of Bone Morphogenetic Protein Receptor-1B (BMPR1B) Gene in Markhoz Goat Breed

Reproductive traits in livestock species are genetically controlled by the action of single genes with a major effect, commonly known as fecundity genes. One of the genes involved in controlling prolificacy is BMPR1B (FecB), a dominant autosomal gene located in chromosome 6 responsible for the fecundity and twinning rate in sheep and goat species. Markhoz goat is a valuable Iranian genetic resource endangered by extinction. Increasing the genetic variability and reproductive performances of Markhoz goat could preserve and enhance its economic value. This study was carried out to detect possible polymorphisms in BMPR1B gene in a sample of 100 Markhoz goats from Iran. DNA samples were screened by PCR-RFLP to assess the presence of the previously reported FecB mutation. Finally, the amplicons from seven goats out of the 100 samples were sequenced. The results showed that all the analyzed individuals did not carry the previously reported FecB mutant allele. However, our findings revealed two novel possible mutations in exon 8 of BMPR1B gene (775A > G and 777G > A) that need further investigations.

miR-130a/TGF-β1 axis is involved in sow fertility by controlling granulosa cell apoptosis

TGF-β1 is a ligand of the TGF-β superfamily and an important cytokine that regulates ovarian functions including follicular development, steroid production, ovulation, luteinization, and female fertility. However, little is known about the regulation of TGF-β1 expression in ovary. Here, we identified that TGF-β1 is a functional target of miR-130a in porcine ovarian granulosa cells (GCs). The 3'-UTR sequence of TGF-β1 gene (1137 bp in length) in Large White (LW) pig was isolated, and multiple RNA regulatory elements (RREs), including several binding motifs of different miRNAs, were identified in this region. Luciferase activity assay showed that miR-130a dramatically suppresses the 3'-UTR luciferase activity of TGF-β1 gene, and further inhibits the expression of TGF-β1 in porcine GCs. FACS revealed that miR-130a acts as a pro-apoptotic factor and promotes GC apoptosis by inhibiting TGF-β1. Two novel linked mutations (-573G > A and -540T  >  C) were identified in the promoter region of ssc-miR-130a, but their polymorphisms are not associated with sow reproductive traits. Importantly, combined genotype analysis with a known mutation (c.1583 A  >  G) in the 3'-UTR of porcine TGF-β1 gene showed a significant association with reproductive performance in LW sow population. Overall, our findings defined a novel regulatory axis, miR-130a/TGF-β1 axis, which is involved in regulating sow fertility.

Polymorphisms in BMPR-IB gene and their association with litter size trait in Chinese Hu sheep

Identification and utilization of sheep major fecundity genes offer opportunities for the increase in litter size, as well as the improvement of production efficiency in livestock industry. BMPR-IB gene belongs to the TGF-β superfamily, and is also considered as a regulator for sheep reproductive performance due to its involvement in the mammalian gametogenesis pathway. This study aimed to detect the variations of BMPR-IB gene in Hu sheep (N = 934) and to evaluate their effects on the litter size trait. qRT-PCR results showed that the mRNA expression level of BMPR-IB in kidney was the highest. And in the tissues of ovary and pituitary, the expression levels of prolific group were significantly higher than that of non-prolific group (p < 0.05). Through DNA sequencing and PCR-RFLP, three SNPs were identified in the genomic region of BMPR-IB gene; the individuals with CC in g.29362047T > C, AA in g.29427689G > A and GG in FecB had better fecundity characterization. Additionally, association analysis indicated that two diplotypes of Hap2/2 and Hap2/4 showed larger litter size. Overall, our results verified several useful markers which would contribute to further development of sheep breeding strategies.

Highly efficient generation of sheep with a defined FecBB mutation via adenine base editing

Base editing has the potential to improve important economic traits in agriculture and can precisely convert single nucleotides in DNA or RNA sequences into minimal double-strand DNA breaks (DSB). Adenine base editors (ABE) have recently emerged as a base editing tool for the conversion of targeted A:T to G:C, but have not yet been used in sheep. ABEmax is one of the latest versions of ABE, which consists of a catalytically-impaired nuclease and a laboratory-evolved DNA-adenosine deaminase. The Booroola fecundity (FecB^B) mutation (g.A746G, p.Q249R) in the bone morphogenetic protein receptor 1B (BMPR1B) gene influences fecundity in many sheep breeds. In this study, by using ABEmax we successfully obtained lambs with defined point mutations that result in an amino acid substitution (p.Gln249Arg). The efficiency of the defined point mutations was 75% in newborn lambs, since six lambs were heterozygous at the FecB^B mutation site (g.A746G, p.Q249R), and two lambs were wild-type. We did not detect off-target mutations in the eight edited lambs. Here, we report the validation of the first gene-edited sheep generated by ABE and highlight its potential to improve economically important traits in livestock.

BMPR-1B, BMP-15 and GDF-9 genes structure and their relationship with litter size in six sheep breeds reared in Egypt

Objective

The aim of this work was to investigate three different mutations; Fec-B, FecX^G, Fec-G^H at three candidate genes; Bone morphogenetic protein receptor IB, Bone morphogenetic protein 15 and Growth Differentiation Factor 9, respectively, in six sheep breeds reared in Egypt namely; Rahmani, Barki, Rahmani X Barki cross, Awassi, Awassi X Suffolk cross, and Ossimi and their association with litter size.

Results

Genomic DNA of 132 sheep was investigated for the Fec-B, FecX^G, and Fec-G^H mutations by Restriction Fragment Length Polymorphism, Single-Stranded Conformation Polymorphism and DNA sequencing. The results revealed that all breeds did not carry Fec-B mutation. On the other side, the mutations of FecX^G, and Fec-G^H were detected in Rahmani, and Rahmani X Barki cross which is associated with the high twinning rate/litter size of Rahmani (1.28) and Rahmani X Barki cross (1.22). While, the average litter size for other breeds had almost a constant values rate over six parities, ranging between 1.00 and 1.04.

Novel Variants in GDF9 Gene Affect Promoter Activity and Litter Size in Mongolia Sheep

Litter size is an economically important trait in sheep breeding. The objectives of this study were as follows: (1) to ascertain if any of the 19 known variants in the BMPRIB, BMP15, and GDF9 genes are present and associated with the litter size of Mongolia sheep; (2) to identify novel variants in GDF9 and perform association analysis; and (3) to validate the effects of these GDF9 promoter variants on the activity of the gene. The results of the 19 known variants showed that the FecB^B affected the litter size of Mongolia sheep (p < 0.001). The association analysis results of novel variants showed that the g.46544883A>G (GenBank accession: NC_040256, the same below) in the 3' untranslated region (3' UTR), the c.1040T>C (Phe347Ser) in the exon 2, and the g.46547859C>T SNP in the promotor of GDF9 were significantly associated with litter size of Mongolia ewes (p < 0.01, p < 0.05, and p < 0.001, respectively). In addition, the GDF9 promoter activity analysis showed that the C allele at the -332 position (g.46547859C>T) could decrease luciferase activity compared with the T allele (p < 0.01). Our findings may facilitate effective marker-assisted selection to increase litter size in Mongolia sheep populations, as well as bring new insights into GDF9 expression.

Effects of Long-Chain Fatty Acyl-CoA Synthetase 1 on Diglyceride Synthesis and Arachidonic Acid Metabolism in Sheep Adipocytes

Long-chain fatty acyl-CoA synthetase (ACSLs) is an essential enzyme for the synthesis of fatty acyl-CoA. ACSL1 plays a key role in the synthesis of triglycerides, phospholipids, and cholesterol esters.

BACKGROUND

In the current study, triglyceride content did not increase after overexpression of the ACSL1 gene.

METHODS

RNA-seq and lipid metabolome profiling were performed to determine why triglyceride levels did not change with ACSL1 overexpression.

RESULTS

Fatty acyl-CoA produced by ACSL1 was determined to be involved in the diglyceride synthesis pathway, and diglyceride content significantly increased when ACSL1 was overexpressed. Moreover, the arachidonic acid (AA) content in sheep adipocytes significantly increased, and the level of cyclooxygenase 2 (COX2) expression, the downstream metabolic gene, was significantly downregulated. Knocking down the ACSL1 gene was associated with an increase in COX2 mRNA expression, as well as an increase in prostaglandin content, which is the downstream metabolite of AA.

CONCLUSIONS

The overexpression of the ACSL1 gene promotes the production of AA via downregulation of COX2 gene expression.

Integrated Hypothalamic Transcriptome Profiling Reveals the Reproductive Roles of mRNAs and miRNAs in Sheep

Early studies have provided a wealth of information on the functions of microRNAs (miRNAs). However, less is known regarding their functions in the hypothalamus involved in sheep reproduction. To explore the potential roles of hypothalamic messenger RNAs (mRNAs) and miRNAs in sheep without FecB mutation, in total, 172 and 235 differentially expressed genes (DEGs) and 42 and 79 differentially expressed miRNAs (DE miRNAs) were identified in polytocous sheep in the follicular phase versus monotocous sheep in the follicular phase (PF vs. MF) and polytocous sheep in the luteal phase versus monotocous sheep in the luteal phase (PL vs. ML), respectively, using RNA sequencing. We also identified several key mRNAs (e.g., POMC, GNRH1, PRL, GH, TRH, and TTR) and mRNA–miRNAs pairs (e.g., TRH co-regulated by oar-miR-379-5p, oar-miR-30b, oar-miR-152, oar-miR-495-3p, oar-miR-143, oar-miR-106b, oar-miR-218a, oar-miR-148a, and PRL regulated by oar-miR-432) through functional enrichment analysis, and the identified mRNAs and miRNAs may function, conceivably, by influencing gonadotropin-releasing hormone (GnRH) activities and nerve cell survival associated with reproductive hormone release via direct and indirect ways. This study represents an integral analysis between mRNAs and miRNAs in sheep hypothalamus and provides a valuable resource for elucidating sheep prolificacy.

Review on Genomic Regions and Candidate Genes Associated with Economically Important Production and Reproduction Traits in Sheep (Ovies aries)

Sheep (Ovis aries) is one of the most economically, culturally, and socially important domestic animals. They are reared primarily for meat, milk, wool, and fur production. Sheep were reared using natural selection for a long period of time to offer these traits. In fact, this production system has been slowing the productivity and production potential of the sheep. To improve production efficiency and productivity of this animal through genetic improvement technologies, understanding the genetic background of traits such as body growth, weight, carcass quality, fat percent, fertility, milk yield, wool quality, horn type, and coat color is essential. With the development and utilization of animal genotyping technologies and gene identification methods, many functional genes and genetic variants associated with economically important phenotypic traits have been identified and annotated. This is useful and presented an opportunity to increase the pace of animal genetic gain. Quantitative trait loci and genome wide association study have been playing an important role in identifying candidate genes and animal characterization. This review provides comprehensive information on the identified genomic regions and candidate genes associated with production and reproduction traits, and gene function in sheep.

Effect of the Booroola fecundity (FecB) gene on the reproductive performance of ewes under assisted reproduction

Reproductive traits are important factors in sheep production. The Booroola fecundity (FecB) gene-the first major gene for prolificacy identified in sheep-has a positive effect on ovulation rates and litter size under natural reproductive conditions. However, the effect of the FecB gene on reproductive performance under assisted reproduction, which uses many artificial hormones, remains unclear. In the present study, we evaluated the effect of FecB (BMPR-1B mutation) on reproductive performance under assisted reproduction, and examined offspring body weight at birth and weaning and survival rate at weaning. There were no differences among three genotype groups (homozygous carrier, BB; heterozygous carrier, B+; non-carrier, ++) in terms of estrus detection rate, time to estrus onset, or estrus duration following estrus synchronization (P > 0.05). The pregnancy rates at 60 d were similar among three genotype groups after artificial insemination (P > 0.05). However, the B allele had an additive effect on litter size (one copy resulted in an increase of 0.88 lambs and two copies produced an additional 0.41 lambs; P < 0.01), and increased lambing and fecundity rates (P < 0.01). After multiple ovulation, the average numbers of recovered embryos per ewe were 9.16 ± 0.79, 8.20 ± 0.77, and 8.44 ± 0.61 in the BB, B+, and ++ ewes, respectively (P > 0.05). There were no differences in the fertilization rate or numbers of grade 1-2 embryos among different groups (P > 0.05). The birth and weaning weights of lambs from BB and B+ ewes were lower than those of lambs born from ++ ewes (P < 0.01) owing to the high fecundity. The survival rate of lambs at weaning did not differ among groups (P > 0.05). Our results indicated that the presence of the B allele had an additive effect on litter size after artificial insemination, but it did not influence the parameters of estrus synchronization and multiple ovulation. Furthermore, the higher prolificacy in ewes carrying the B allele was associated with a reduction in offspring body weight at birth and weaning.

Genetic Effects of Single Nucleotide Polymorphisms in the Goat GDF9 Gene on Prolificacy: True or False Positive?

Simple Summary

As an important regulator factor, which was secreted by female oocytes, the growth differentiation factor 9 (GDF9) plays an essential role during the growth and differentiation of ovarian follicles. Single nucleotide polymorphisms (SNPs) within the GDF9 gene have been found to be involved in reproductive traits in livestock, and some of these mutations have been used as the effective makers in animal molecular breeding. However, it is remarkable that the SNPs of the goat GDF9 gene have not been systematically sorted and analyzed from the reported studies, which leads to an inability to find effective loci that could be applied in improving the prolificacy of goats via the molecular breeding method. In this study, we gathered and sorted 45 SNPs of the goat GDF9 gene from all relevant studies and the National Center for Biotechnology Information Search database (NCBI), and especially analyzed and discussed the relationship between part controversial and potentially effective SNPs and the reproductive traits. The results indicated that non-synonymous SNPs A240V, Q320P, and V397I and synonymous SNPs L61L, N121N, and L141L were six “true” positive SNPs in improving goat fertility. Nevertheless, the regulation pathways and the specific mechanism of these six SNPs on goat fecundity are not clear, which still need further study in more goat breeds and a large sample size. These results provided an effective tool for follow-up research studies on the molecular genetic breeding of goats’ reproductive traits.

Abstract

Goat reproductive traits are complex quantitative traits controlled by polygenes and multipoint. To date, some high-fertility candidate genes in livestock have been unearthed and the growth differentiation factor 9 (GDF9) gene is one of them, which plays a crucial role in early folliculogenesis. According to the relevant previous studies and the National Center for Biotechnology Information Search database (NCBI), a total of 45 single nucleotide polymorphisms (SNPs) have been detected in the goat GDF9 gene, but which one or which ones have important effects on goat fecundity is still uncertain. Hence, in order to find effective molecular markers for goat genetic breeding and accelerate the goat improvement, this study summarized and classified the above 45 SNPs into four kinds, as well as compared and analyzed the same SNP effects and the different SNPs linkage effects on the reproductive traits in different goat breeds. Since there were many SNPs in the goat GDF9 gene, only 15 SNPs have been identified in more than 30 goat breeds worldwide and they showed different effects on the litter size. Therefore, this study mainly chose these 15 SNPs and discussed their relationship with goat productivity. Results showed that three non-synonymous SNPs A240V, Q320P, and V397I and three synonymous ones L61L, N121N, and L141L played a “true” role in the litter size trait in many goat breeds around the world. However, the regulatory mechanisms still need further research. These results provide an effective tool for follow-up research developing the goat molecular breeding strategies and improving the goat reproductive traits.

Identification of mutations in BMP15 and GDF9 genes associated with prolificacy of Markhoz goats

The Markhoz is a local goat breed in the Kurdistan region of Iran. The mohair obtained from these animals plays an important cultural role and is used for making local clothes in the Kurdistan region. This breed is a low-fecundity local goat, and searching for genes associated with fertility of these goats is important for their breeding. Moreover, this research is complementary to prior studies of candidate genes associated with fertility. The growth differentiation factor 9 (GDF9) and bone morphogenetic protein 15 (BMP15) are attractive candidates expressed by the oocyte and are associated with increased ovulation rate in sheep. However, there are no reports on single nucleotide polymorphisms associated with fertility of Markhoz goats. Hence, we studied these candidate genes and found two novel mutations (g.233C>A and g.755T>G) in GDF9 exon I and in BMP15 exon II, respectively. Furthermore, we investigated their association with prolificacy. These nucleotide changes are detectable with the PCR-RFLP method and can be used in the screening for highly fecund goats. Both of the mutations significantly increased litter size in heterozygote form for BMP15 and homozygote form for GDF9 in this goat breed. Homozygote females for the BMP15 mutation were not identified in the Markhoz breed, which is similar to the situation found in Belclare sheep, small-tailed Han sheep, and Jining Grey goats.

Effect of BMPRIB gene on litter size of sheep in China: A meta-analysis

Genes or genetic markers related to litter size have been studied for many years in gene selection or marker-assisted selection experiments. The bone morphogenetic protein receptor IB (BMPRIB) gene is one of the candidate genes for increasing litter size of sheep. Results of studies with different sheep breeds in China have been inconsistent with some confirming significant associations between the BMPRIB gene polymorphism and litter size and with other results being inconsistent with there not being an association. In the present study, a meta-analysis was conducted evaluating 21 studies in which 5089 samples were utilized to analyze the genetic effects of BMPRIB genes on litter size in different sheep breeds in China. Results for weighted mean differences (WMD) among BMPRIB genotypes indicated there was an association between the BPRIB gene polymorphism and litter size. Effects of the BMPRIB gene on litter size are remarkably consistent in many sheep populations of China, with each gene copy being associated with an increase in litter size of 0.4 to 0.5 lambs. There are, however, some populations in which there is no effect of a second copy of the B allele. An example is the Zeller black sheep. Furthermore, in the Tan sheep and its crosses, there was no effect of the BMPRIB gene on litter size. In conclusion, with this study, there is a summarizing of magnitude of BMPRIB gene effects on litter size for sheep breeds in China, and these results provide reference information for consideration in indigenous sheep breeding programs.

Quantitative expression of mRNA encoding BMP/SMAD signalling genes in the ovaries of Booroola carrier and non-carrier GMM sheep

The GMM sheep is a carrier of Booroola fecundity (FecB) gene, which produces the twins and triplets in one lambing. The homozygous carrier GMM (FecB^BB ), non-carrier GMM and Malpura (FecB⁺⁺ ) ewes were synchronized by progesterone sponges, and the plasma progesterone concentration was measured by RIA. The results showed that the progesterone concentration did not differ significantly (p > .05) in homozygous carrier GMM (5.74 ± 1.2 ng/ml), non-carrier GMM (5.42 ± 1.4 ng/ml) and non-carrier Malpura ewes (5.67 ± 1.5 ng/ml). Further, quantitative expression of BMP factors/receptors and SMAD signalling genes were analysed in the ovaries of sheep by qRT-PCR. The study showed that the expression of BMP2 was slightly higher (p > .05) in carrier GMM than that of non-carrier GMM, but it was almost similar to Malpura ewes. Expression of BMP4 and BMP7 was significantly higher (p < .001; p < .05) in carrier GMM than that of non-carrier GMM and Malpura ewes. Although BMP6 expression was higher (p > .05) in carrier GMM than that of non-carrier GMM, but lower (p > .05) than the Malpura ewes. Expression of BMP15 (p < .05), GDF5 (p < .01) and GDF9 (p < .05) was significantly higher in carrier GMM than non-carrier GMM ewes. Surprisingly, BMPR1B expression was significantly higher (p < .001) in non-carrier GMM and Malpura than the carrier GMM ewes, while TGFβRI did not differ significantly (p > .05) among both GMM genotypes. On the other hand, expression of BMPR1A (p > .05) and BMPRII (p < .05) was higher in carrier GMM than the non-carrier GMM, but significantly lower (p < .001) than the Malpura ewes. It was interesting to note that the expression of SMAD1 (p > .05), SMAD2 (p < .001), SMAD3 (p < .05), SMAD4 (p < .001), SMAD5 (p < .001) and SMAD8 (p < .001) was lower in the carrier GMM than that of non-carrier GMM ewes. It is concluded that the FecB mutation alters the expression of BMPR1B and SMAD signalling genes in the ovaries of homozygous carrier GMM ewes.