Highly prolific Booroola sheep have a mutation in the intracellular kinase domain of bone morphogenetic protein IB receptor (ALK-6) that is expressed in both oocytes and granulosa cells

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.

Are two beneficial mutations (p.Q249R and 90-bp Indel) within the ovine BMPRIB gene associated with growth traits?

Background

The problem of achieving economic efficiency in sheep breeding can be largely solved by increasing sheep productivity. Recently, the BMPRIB gene has been revealed by GWAS as a potential candidate gene for sheep body morphometric traits. Therefore, the present study aimed to investigate whether genetic polymorphisms (p.Q249R SNP and 90-bp deletion) in the BMPRIB gene are associated with sheep growth traits.

Methods

PCR-based genotyping was performed on 1,875 sheep, including 1,191 Guiqian semi-fine wool (GQSFW), 560 Luxi Blackhead (LXBH), 55 Lanzhou fat-tailed (LZFT), and 69 Weining (WN) sheep. Genotype-phenotype association was assessed using the independent samples t-test and ANOVA. The significance level was set at αoriginal < 0.05. The threshold p-value for significance was adjusted after correction for multiple comparisons using the Bonferroni correction.

Results

After the Bonferroni correction, it was found that individuals with FecB+/FecB+ genotypes of the p.Q249R had significantly better growth traits in LXBH ewe lambs, including the body length, chest width, paunch girth, cannon circumference, and hip width (P<0.0005). Meanwhile, associations were observed between 90-bp deletion polymorphism and several growth traits (body length, body height, chest depth, and canon circumference) in GQSFW ewe adults after the Bonferroni correction (P < 0.0002), and individuals with the "DD" genotypes had greater growth traits.

Conclusion

Our findings align with the experimental observations from GWAS, which identified the BMPRIB gene as a potential candidate gene for body measurement traits. These findings not only confirm the previous study's results but also expand on them. Therefore, further investigations regarding the impact of BMPRIB polymorphisms on growth traits are necessary in other sheep breeds.

Detection of Booroola Polymorphism of Bone Morphogenetic Protein Receptor 1b and Embrapa Polymorphism of Growth Differentiation Factor 9 in Sheep in Thailand

This study aimed to investigate the appearance and frequencies of the Booroola polymorphism of the bone morphogenetic protein receptor 1b (BMPR1B) gene (FecB) and the Embrapa polymorphism of the growth differentiation factor 9 (GDF9) gene (FecGE) in sheep in Thailand. A total of 454 crossbred sheep blood samples were collected from four provinces in Thailand during August 2022 to July 2023. The polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method was used to identify the FecB and FecGE genotypes. The history of ewe birth types was collected from the owners to analyze the association between fecundity (Fec) genotypes and the history of birth types. The genotypic frequencies of FecB for homozygous genotype (B/B), heterozygous genotype (+/B), and wildtype (+/+) were 0.22%, 1.54%, and 98.24%, respectively. Meanwhile, the genotypic frequencies of FecGE for homozygous genotype (E/E), heterozygous genotype (+/E), and wildtype (+/+) were 0.00%, 2.42%, and 97.58%, respectively. Furthermore, three ewes exhibited both FecB and FecGE genotypes. Fisher's exact test revealed that possession of the FecB genotype was associated with multiple births (p < 0.01). Both FecB and FecGE mutations were identified in crossbred sheep in Thailand. Sheep containing FecB allele could be alternative candidates to be selected to improve the prolificacy of crossbred sheep in Thailand.

New insights into the (epi)genetics of twinning

Spontaneous dizygotic (DZ) twins, i.e. twins conceived without the use of ARTs, run in families and their prevalence varies widely around the globe. In contrast, monozygotic (MZ) twins occur at a constant rate across time and geographical regions and, with some rare exceptions, do not cluster in families. The leading hypothesis for MZ twins, which arise when a zygote splits during preimplantation stages of development, is random occurrence. We have found the first series of genes underlying the liability of being the mother of DZ twins and have shown that being an MZ twin is strongly associated with a stable DNA methylation signature in child and adult somatic tissues. Because identical twins keep this molecular signature across the lifespan, this discovery opens up completely new possibilities for the retrospective diagnosis of whether a person is an MZ twin whose co-twin may have vanished in the early stages of pregnancy. Here, we summarize the gene finding results for mothers of DZ twins based on genetic association studies followed by meta-analysis, and further present the striking epigenetic results for MZ twins.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity

STUDY QUESTION

Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins?

SUMMARY ANSWER

We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3.

WHAT IS KNOWN ALREADY

The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures.

STUDY DESIGN, SIZE, DURATION

We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls).

PARTICIPANTS/MATERIALS, SETTING, METHODS

Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation.

MAIN RESULTS AND THE ROLE OF CHANCE

This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle.

LARGE SCALE DATA

The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/).

LIMITATIONS, REASONS FOR CAUTION

Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility.

WIDER IMPLICATIONS OF THE FINDINGS

About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility.

STUDY FUNDING/COMPETING INTEREST(S)

Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Skłodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests.

TRIAL REGISTRATION NUMBER

N/A.

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.

Comparison of the growth performance between pure Moghani sheep and crosses with Texel or Booroola sheep carrying major genes contributing to muscularity and prolificacy

This study aimed to compare growth performance between Moghani sheep and crossbred lambs resulting from crossbreeding between Moghani pure breed ewes and the lines of rams e.g., Texel Tamlet, Texel Dalzell, Booroola Merino, and Booroola Romney. The first visible phenotypic characteristic was the presence of lean tail in all F1 crossbred lambs, whereas Moghani pure sheep is a well-known large fat-tailed breed. Moreover, the first generation of backcross (BC1) lambs from mating four types of F1 crossbred rams with Moghani pure ewes revealed lean-tailed to short fat-tailed. Comparative results showed that the F1 crossbred lambs had significantly (p < 0.0001) greater birth weight (BW) than the Moghani pure breed lambs. Despite no significant differences observed between Moghani pure breed sheep and its F1 crossbred lambs for body weight at pre-weaning, but F1 crossbred lambs achieved significantly (p < 0.0001) greater body weight after weaning compared to Moghani sheep. The growth performance of BC1 lambs was outperformed than F1 crossbred and Moghani sheep. These results encourage the continuation of the Moghani sheep crossbreeding programs to improve overall lamb growth, particularly post-weaning and to benefit from a better reproductive efficiency by elimination or reduction of the fat tail.

Assessment of FecGE genotypes on reproductive traits in Brazilian Morada Nova and Santa Inês sheep

The aim of this study was to evaluate the effect of the polymorphic FecGE allele on reproductive traits in Santa Inês and Morada Nova ewes. The traits evaluated were as follows: total progeny weights at birth (PWB) and weaning (PWW) and progeny survival rates at birth (PSRB) and weaning (PSRW). A total of 389 animals, belonging to two Santa Inês herds and one Morada Nova herd, were genotyped. There was a difference between the averages for all the traits studied regarding type of parturition, herd/breed, genotype/herd, and genotype/type of parturition. For each additional progeny, if the female was FecGE/E, the PWB decreased by 1.02 kg and the PWW by 3.16 kg, also with a 0.04% reduction in PSRB and no change in PSRW. If the female was FecGE/+, the reduction in PWB was 0.24 kg, with an increase in PSRW by 0.11%, but no change in PWW and PSRB. In general, these results demonstrate that FecG+/+ females have a better ability to increase their number of progenies without reducing PWB and PWW (also similar to FecGE/+). Thus, it is suggested that further studies on the association between the traits of interest and candidate genes in sheep should be carried out so that the regions which have the greatest effect on the expression of these traits are actually identified. It was not possible to verify the effect of the FecGE allele on the PWB, PWW, PSRB, and PSRW in these Morada Nova and Santa Inês herds.

CRISPR/Cas mediated disruption of BMPR-1B gene and introduction of FecB mutation into the Caprine embryos using Easi-CRISPR strategy

Bone Morphogenetic Proteins play a significant role in ovarian physiology and contribute to the reproductive fitness of mammals. The BMPR-1B/FecB mutation, a loss of function mutation increases litter size by 1-2 with each number of mutated alleles in sheep. Considering demand-supply gap of the meat industry, and low replacement rate of indigenous caprine species, the conservative BMPR-1B locus can be explored, and FecB mutated goats can be produced. The experiment one produced CRISPR/Cas mediated KO transferable caprine embryos, and experiment two generated caprine embryos with desired FecB mutation using Easi-CRISPR strategy. In the KO experiment, Cas9 and BMPR-1B guide RNA (100:100ng/ul) were electroporated into single stage caprine zygotes at 750V, 10 ms and 1pulse using Neon transfection system. In the second experiment, phosphorothioate (PS) modified single-stranded oligodeoxynucleotide (ssODN) was used as an HDR template along with CRISPR components (100:100ng/ul, ssODN 100ng/ul). The precise time and method of electroporation, RNP format of CRISPR components and PS modified asymmetric ssODN were the factors that affected the production of mosaicism free BMPR-1B edited caprine embryos. The editing efficiency of KO and KI experiments was 68.52 and 63.16% respectively, and successful production of goats with higher mean ovulation rate can be realized with addition of embryo transfer technology to these experiments.

FecB mutation and litter size are associated with a 90-base pair deletion in BMPR1B in East Friesian and Hu crossbred sheep

Litter size is a critical economic trait in livestock, but only a few studies have focused on associated indel mutations in BMPR1B, a key regulator of ovulation and litter size in sheep. We evaluated the effects of BMPR1B mutations on the reproductive performance of sheep. We used Hu, East Friesian, and East Friesian/Hu crossbred sheep as experimental subjects and identified a novel 90 bp deletion in BMPR1B, which coincides with the c.746A > G (FecB mutation) genotype. The correlation between the two loci and litter size was then evaluated. We identified three genotypes for the Del-90bp locus, namely, II, ID, and DD, and three genotypes for the c.746A > G locus, namely ++, B+, and BB. Both Del-90bp and c.746A > G significantly affected the litter size of Hu and East Friesian/Hu crossbred sheep. Linkage disequilibrium analysis revealed a strong linkage disequilibrium between these loci in Hu sheep and the F1 population (r2 > 0.33), which suggests that detecting this 90 bp deletion might be a simple method to identify the likely carriers of c.746A > G. However, the function of this 90-bp deletion still needs further exploration. We provide genetic data that can be used as a reference for the breeding of improved prolific traits in sheep.

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.

A meta-analysis of mRNA expression profiling studies in sheep with different FecB genotypes

The FecB mutation in the sheep BMPRIB is strongly correlated with high ovulation traits but its mechanism remains unclear. This study explored differentially expressed genes (DEGs) and their associated molecular mechanisms that may be involved in FecB mutation-induced high ovulation from the perspective of the hypothalamic-pituitary-gonadal (HPG) axis by conducting a systematic review and meta-analysis. The PubMed, EMBASE, CNKI, WanFang, and CBM databases were searched for eligible articles published before August 2022, focusing on mRNA sequencing of different tissues in the HPG axis in sheep with different FecB genotypes. A total of 6555 DEGs were identified from the analysis of six published articles and experimental results from our laboratory. The DEGs were screened by vote-counting rank and robust rank aggregation. Among these, in the follicular phase, FKBP5, CDCA7 and CRABP1 were upregulated in the hypothalamus. INSM2 was upregulated, while LDB3 was downregulated in the pituitary. CLU, SERPINA14, PENK, INHA and STAR were upregulated, while FERMT2 and NPY1R were downregulated in the ovary. On the HPG axis, TAC1 was upregulated and NPNT was downregulated. Many DEGs were found in sheep with different FecB genotypes. The genes FKBP5, CDCA7, CRABP1, INSM2, LDB3, CLU, SERPINA14, PENK, INHA, STAR, FERMT2, NPY1R, TAC1 and NPNT, may be associated with FecB mutation-induced high ovulation in different tissues. These candidate genes will further improve the mechanism of multiple fertility traits induced by the FecB mutation from the perspective of the HPG axis.

Detection of carrier Booroola (FecB) allele in BMPR1B gene of MEGA (Merino × Garut) sheep and its association with growth traits

BACKGROUND

Bone morphogenetic protein receptor 1B (BMPR1B) gene is one of candidate genes for reproductive and growth traits in sheep. The present study was aimed to detect the Booroola (Fec^B) allele in BMPR1B gene and its association with growth traits in MEGA (Merino × Garut) sheep. A total of 82DNA samples collected from individual lamb (mixed-sex) blood were genotyped for allelic polymorphism using a PCR-RFLP method.

RESULTS

The PCR analysis in BMPR1B gene resulted the amplicons with size of140 bp. The RFLP analysis with AvaII restriction enzymeresultedtwo allelic types of wildtype (A/Fec⁺) and mutant or Booroola (G/Fec^B) with frequency of 0.89 and 0.11, respectively. However, the genetic diversity in BMPR1B/AvaII gene of animal studies was categorized tolow category (PIC = 0.18)and under in a genetic equilibrium (χ² = 1.25).

CONCLUSIONS

Itshowed us that carrying Fec^B allele in the heterozygous sheep were not associated with growth traits in MEGA sheep.

Impact of FecB Mutation on Ovarian DNA Methylome in Small-Tail Han Sheep

Booroola fecundity (FecB) gene, a mutant of bone morphogenetic protein 1B (BMPR-1B) that was discovered in Booroola Merino, was the first prolificacy gene identified in sheep related to increased ovulation rate and litter size. The mechanism of FecB impact on reproduction is unclear.

METHODS

In this study, adult Han ewes with homozygous FecB(B)/FecB(B) mutations (Han BB group) and ewes with FecB(+)/FecB(+) wildtype (Han ++ group) were selected. Methylated DNA immunoprecipitation and high-throughput sequencing (MeDIP-seq) was used to identify differences in methylated genes in ovary tissue.

RESULTS

We examined differences in DNA methylation patterns between HanBB and Han ++ sheep. In both sheep, methylated reads were mainly distributed at the gene body regions, CpG islands and introns. The differentially methylated genes were enriched in neurotrophy in signaling pathway, Gonadotropin Releasing Hormone (GnRH) signaling pathway, Wnt signaling pathway, oocyte meiosis, vascular endothelial growth factor (VEGF) signaling pathway, etc. Differentially-methylated genes were co-analyzed with differentially-expressed mRNAs. Several genes which could be associated with female reproduction were identified, such as FOXP3 (forkhead box P3), TMEFF2 (Transmembrane Protein with EGF Like and Two Follistatin Like Domains 2) and ADAT2 (Adenosine Deaminase TRNA Specific 2).

CONCLUSIONS

We constructed a MeDIP-seq based methylomic study to investigate the ovarian DNA methylation differences between Small-Tail Han sheep with homozygous FecB mutant and wildtype, and successfully identified FecB gene-associated differentially-methylated genes. This study has provided information with which to understand the mechanisms of FecB gene-induced hyperprolificacy in sheep.

Paracrine regulation of granulosa cell development in the antral follicles in mammals

Background

Development of ovarian follicles is regulated by a complex interaction of intra- and extra-follicular signals. Oocyte-derived paracrine factors (ODPFs) play a central role in this process in cooperation with other signals.

Methods

This review provides an overview of the recent advances in our understanding of the paracrine regulation of antral follicle development in mammals. It specifically focuses on the regulation of granulosa cell development by ODPFs, along with other intrafollicular signals.

Main Findings

Bi-directional communication between oocytes and surrounding cumulus cells is a fundamental mechanism that determines cumulus cell differentiation. Along with estrogen, ODPFs promote the expression of forkhead box L2, a critical transcription factor required for mural granulosa cells. Follicle-stimulating hormone (FSH) facilitates these processes by stimulating estrogen production in mural granulosa cells.

Conclusion

Cooperative interactions among ODPFs, FSH, and estrogen are critical in determining the fate of cumulus and mural granulosa cells, as well as the development of oocytes.

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.

Two Retrotransposon Elements in Intron of Porcine BMPR1B Is Associated with Phenotypic Variation

It has been established that through binding to bone morphogenetic proteins (BMPs), bone morphogenetic protein receptor I B (BMPR1B) can mediate transforming growth factor β (TGF-β) signal transduction, and is involved in the regulation of several biological processes, such as bone and muscle formation and homeostasis, as well as folliculogenesis. Also known as FecB, BMPR1B has been reported as the major gene for sheep prolificacy. A number of previous studies have analyzed the relationship between single nucleotide polymorphisms (SNPs) in this gene and its related performance. In recent years, with the illustration of the effect of retrotransposon insertion on the expression of the proximal genes or phenotypic variation, retrotransposon insertion polymorphisms (RIPs) have been used as a novel type of molecular marker in the evaluation of evolution, population structure and breeding of plant and domestic animals. In this study, the RIPs in porcine BMPR1B gene were excavated, and thereafter verified using a comparative genome and polymerase chain reaction (PCR). The potential effects of phenotype, gene expression and functions related to RIPs were also explored. The results showed that 13 distinct RIPs were identified in introns of porcine BMPR1B. Among these, only BMPR1B-SINE-RIP9 and BMPR1B-LINE-RIP13 displayed a close relationship with the growth traits of Large White pigs. Moreover, the total number of BMPR1B-SINE+/+-RIP9 individuals born was found to be significantly higher than that of SINE−/− (p < 0.05). These two RIPs showed an obvious distribution pattern among Chinese indigenous breeds and Western commercial breeds. The expression of BMPR1B in ovaries of adult BMPR1B-SINE+/+-RIP9 Sushan pigs was found to be significantly higher in comparison to those of BMPR1B-SINE−/−-RIP9 (p < 0.05). SINE insertion of BMPR1B-SINE-RIP9 and LINE insertion of BMPR1B-LINE-RIP13 were observed to significantly increase the activity of Octamer binding transcription factor 4 (OCT4) minipromoter in CHO and C2C12 cells (p < 0.01). Therefore, these two RIPs could serve as useful molecular markers for modulating the growth or reproductive traits in assisted selection of pig breeding, while the mechanisms of the insertion function should be studied further.

Are Copy Number Variations within the FecB Gene Significantly Associated with Morphometric Traits in Goats?

The Booroola fecundity (FecB) gene is a major fertility-related gene first identified in Booroola sheep. Numerous studies have investigated whether the FecB gene is a major fecundity gene in goats or whether there are other genes that play a critical role in goat fertility. Nevertheless, little attention has been paid to the role of the FecB gene in the body morphometric traits of goats, despite the positive relationship discerned between litter size and growth. We identified five copy number variations (CNVs) within the FecB gene in 641 goats, including 318 Shaanbei white cashmere (SBWC) goats, 203 Guizhou Heima (GZHM) goats, and 120 Nubian goats, which exhibited different distributions among these populations. Our results revealed that these five CNVs were significantly associated with goat morphometric traits (p < 0.05). The normal type of CNV3 was the dominant type and displayed superior phenotypes in both litter size and morphometric traits, making it an effective marker for goat breeding. Consequently, LD blocks in the region of 10 Mb upstream and downstream from FecB and potential transcription factors (TFs) that could bind with the CNVs were analyzed via bioinformatics. Although no significant LD block was detected, our results illustrated that these CNVs could bind to growth-related TFs and indirectly affect the growth development of the goats. We identified potential markers to promote litter size and growth, and we offer a theoretical foundation for further breeding work.

Sequencing Reveals Population Structure and Selection Signatures for Reproductive Traits in Yunnan Semi-Fine Wool Sheep (Ovis aries)

Yunnan semi-fine wool sheep are among the most important cultivated sheep breeds in China. However, their population structure, genetic characteristics and traits of interest are poorly studied. In this study, we systematically studied the population characteristics and selection signatures of 40 Yunnan semi-fine wool sheep using SNPs obtained from whole-genome resequencing data. A total of 1393 Gb of clean data were acquired. The mapping rate against the reference genome was 91.23% on average (86.01%–92.26%), and the average sequence depth was 9.51X. After filtering, 28,593,198 SNPs and 4,725,259 indels with high quality were obtained. The heterozygosity rate, inbreeding coefficient and effective population size of the sheep were calculated to preliminarily explore their genetic characteristics. The average heterozygosity rate was 0.264, the average inbreeding coefficient was 0.0099, and the effective population size estimated from the heterozygote excess (HE) was 242.9. Based on the Tajima’s D and integrated haplotype score (iHS) approaches, 562 windows and 11,356 core SNPs showed selection signatures in the Yunnan semi-fine wool sheep population. After genome annotation and gene enrichment analysis, we found traces of early domestication in sensory organs, behavioural activity and the nervous system as well as adaptive changes in reproductive and wool traits under selection in this population. Some selected genes related to litter size, including FSHR, BMPR1B and OXT, were identified as being under selection. Specific missense mutations of the FSHR gene that differed from the reference genome were also identified in the population, and we found some SNP variations that may affect litter size. Our findings provide a theoretical basis for the conservation and utilization of Yunnan semi-fine wool sheep. Furthermore, our results reveal some changes common to sheep after domestication and provide a new opportunity to investigate the genetic variation influencing fecundity within a population evolving under artificial selection.

Endothelial cell-derived fibroblast growth factor-18 regulates ovarian function in sheep

Increasing the efficiency of farm animal reproduction is necessary to reduce the environmental impact of food production systems. One approach is to increase the number of healthy eggs (oocytes) produced per female for fertilization, thus it is important to understand factors that decrease oocyte health. One paracrine factor that decreases ovarian follicle growth is fibroblast growth factor 18 (FGF18) secreted by cells in the theca layer of the ovarian follicle, however the factors that regulate FGF18 secretion are unknown. In this study we hypothesized that FGF18 secretion is controled by intrafollicular factors and is linked to fertility, which we tested by using cell culture and sheep genetic models in vivo. Separation of theca cell populations revealed that FGF18 messenger RNA (mRNA) is located mainly in thecal endothelial rather than endocrine cells, and immunohistochemistry localized FGF18 protein to microvessels in the theca layer in situ. Culture of ovine theca-derived endothelial cells was used to demonstrate stimulation of FGF18 mRNA and protein abundance by bone morphogenetic protein 4 (BMP4), a growth factor derived from theca endocrine cells. Taking advantage of a sheep genetic model, we demonstrate reduced ovarian and peripheral FGF18 concentrations in the hyperprolific Booroola ewe harboring the FecB^B mutation in BMPR1B. These data suggest a novel control of fertility by follicular endothelial cells, in which theca endocrine cells secrete BMP4 that stimulates the secretion of FGF18 from thecal endothelial cells, which in turn diffuses into the granulosa cell layer and promotes apoptosis.

Comparative Transcriptomic Analysis of Hu Sheep Pituitary Gland Prolificacy at the Follicular and Luteal Phases

The pituitary gland directly regulates the reproduction of domestic animals. Research has increasingly focused on the potential regulatory mechanism of non-coding RNA in pituitary development. Little is known about the differential expression pattern of lncRNAs in Hu sheep, a famous sheep breed with high fecundity, and its role in the pituitary gland between the follicular phase and luteal phase. Herein, to identify the transcriptomic differences of the sheep pituitary gland during the estrus cycle, RNA sequencing (RNA-Seq) was performed. The results showed that 3529 lncRNAs and 16,651 mRNAs were identified in the pituitary gland. Among of them, 144 differentially expressed (DE) lncRNA transcripts and 557 DE mRNA transcripts were screened in the follicular and luteal phases. Moreover, GO and KEGG analyses demonstrated that 39 downregulated and 22 upregulated genes interacted with pituitary functions and reproduction. Lastly, the interaction of the candidate lncRNA XR_001039544.4 and its targeted gene LHB were validated in sheep pituitary cells in vitro. LncRNA XR_001039544.4 and LHB showed high expression levels in the luteal phase in Hu sheep. LncRNA XR_001039544.4 is mainly located in the cytoplasm, as determined by FISH analysis, indicating that XR_001039544.4 might act as competing endogenous RNAs for miRNAs to regulate LHB. LncRNA XR_001039544.4 knockdown significantly inhibited LH secretion and cell proliferation. LncRNA XR_001039544.4 may regulate the secretion of LH in the luteal-phase pituitary gland via affecting cell proliferation. Taken together, these findings provided genome-wide lncRNA- and mRNA-expression profiles for the sheep pituitary gland between the follicular and luteal phases, thereby contributing to the elucidation of the molecular mechanisms of pituitary function.

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.

A frameshift deletion in the GDF9 gene in Icelandic Loa prolific sheep

Several genetic variants have been shown to affect the mean number of offspring in different sheep breeds. Here, we analyzed samples from Icelandic sheep with the aim of identifying the genetic cause of the Icelandic Loa phenotype using three previously identified prolificacy genes as candidates. We demonstrate that a 4-bp frameshift deletion positioned in the mature region of the GDF9 protein in the Loa animals is a likely causal mutation for the observed increase in prolificacy; however, sequencing showed that not all ewes with a high number of offspring carried the deletion, suggesting the presence of a second mutation segregating within this group of animals.

Meta-analysis of genome-wide association studies for litter size in sheep

Litter size and ovulation rate are important reproduction traits in sheep and have important impacts on the profitability of farm animals. To investigate the genetic architecture of litter size, we report the first meta-analysis of genome-wide association studies (GWAS) using 522 ewes and 564,377 SNPs from six sheep breeds. We identified 29 significant associations for litter size which 27 of which have not been reported in individual GWAS for each population. However, we could confirm the role of BMPR1B in prolificacy. Our gene set analysis discovered biological pathways related to cell signaling, communication, and adhesion. Functional clustering and enrichment using protein databases identified epidermal growth factor-like domain affecting litter size. Through analyzing protein-protein interaction data, we could identify hub genes like CASK, PLCB4, RPTOR, GRIA2, and PLCB1 that were enriched in most of the significant pathways. These genes have a role in cell proliferation, cell adhesion, cell growth and survival, and autophagy. Notably, identified SNPs were scattered on several different chromosomes implying different genetic mechanisms underlying variation of prolificacy in each breed. Given the different layers that make up the follicles and the need for communication and transfer of hormones and nutrients through these layers to the oocyte, the significance of pathways related to cell signaling and communication seems logical. Our results provide genetic insights into the litter size variation in different sheep breeds.

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.

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.

Trends towards revealing the genetic architecture of sheep tail patterning: Promising genes and investigatory pathways

Different sheep breeds have evolved after initial domestication, generating various tail phenotypic patterns. The phenotypic diversity of sheep tail patterns offers ideal materials for comparative analysis of its genetic basis. Evolutionary biologists, animal geneticists, breeders, and producers have been curious to clearly understand the underlying genetics behind phenotypic differences in sheep tails. Understanding the causal gene(s) and mutation(s) underlying these differences will help probe an evolutionary riddle, improve animal production performance, promote animal welfare, and provide lessons that help comprehend human diseases related to fat deposition (i.e., obesity). Historically, fat tails have served as an adaptive response to aridification and climate change. However, the fat tail is currently associated with compromised mating and animal locomotion, fat distribution in the animal body, increased raising costs, reduced consumer preference, and other animal welfare issues such as tail docking. The developing genomic approaches provide unprecedented opportunities to determine causal variants underlying phenotypic differences among populations. In the last decade, researchers have performed several genomic investigations to assess the genomic causality underlying phenotypic variations in sheep tails. Various genes have been suggested with the prominence of several potentially significant causatives, including the BMP2 and PDGFD genes associated with the fat tail phenotype and the TBXT gene linked with the caudal vertebrae number and tail length. Although the potential genes related to sheep tail characteristics have been revealed, the causal variant(s) and mutation(s) of these high-ranking candidate genes are still elusive and need further investigation. The review discusses the potential genes, sheds light on a knowledge gap, and provides possible investigative approaches that could help determine the specific genomic causatives of sheep tail patterns. Besides, characterizing and revealing the genetic determinism of sheep tails will help solve issues compromising sheep breeding and welfare in the future.

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.

Meta-Analysis of Heifer Traits Identified Reproductive Pathways in Bos indicus Cattle

Fertility traits measured early in life define the reproductive potential of heifers. Knowledge of genetics and biology can help devise genomic selection methods to improve heifer fertility. In this study, we used ~2400 Brahman cattle to perform GWAS and multi-trait meta-analysis to determine genomic regions associated with heifer fertility. Heifer traits measured were pregnancy at first mating opportunity (PREG1, a binary trait), first conception score (FCS, score 1 to 3) and rebreeding score (REB, score 1 to 3.5). The heritability estimates were 0.17 (0.03) for PREG1, 0.11 (0.05) for FCS and 0.28 (0.05) for REB. The three traits were highly genetically correlated (0.75-0.83) as expected. Meta-analysis was performed using SNP effects estimated for each of the three traits, adjusted for standard error. We identified 1359 significant SNPs (p-value < 9.9 × 10-6 at FDR < 0.0001) in the multi-trait meta-analysis. Genomic regions of 0.5 Mb around each significant SNP from the meta-analysis were annotated to create a list of 2560 positional candidate genes. The most significant SNP was in the vicinity of a genomic region on chromosome 8, encompassing the genes SLC44A1, FSD1L, FKTN, TAL2 and TMEM38B. The genomic region in humans that contains homologs of these genes is associated with age at puberty in girls. Top significant SNPs pointed to additional fertility-related genes, again within a 0.5 Mb region, including ESR2, ITPR1, GNG2, RGS9BP, ANKRD27, TDRD12, GRM1, MTHFD1, PTGDR and NTNG1. Functional pathway enrichment analysis resulted in many positional candidate genes relating to known fertility pathways, including GnRH signaling, estrogen signaling, progesterone mediated oocyte maturation, cAMP signaling, calcium signaling, glutamatergic signaling, focal adhesion, PI3K-AKT signaling and ovarian steroidogenesis pathway. The comparison of results from this study with previous transcriptomics and proteomics studies on puberty of the same cattle breed (Brahman) but in a different population identified 392 genes in common from which some genes-BRAF, GABRA2, GABR1B, GAD1, FSHR, CNGA3, PDE10A, SNAP25, ESR2, GRIA2, ORAI1, EGFR, CHRNA5, VDAC2, ACVR2B, ORAI3, CYP11A1, GRIN2A, ATP2B3, CAMK2A, PLA2G, CAMK2D and MAPK3-are also part of the above-mentioned pathways. The biological functions of the positional candidate genes and their annotation to known pathways allowed integrating the results into a bigger picture of molecular mechanisms related to puberty in the hypothalamus-pituitary-ovarian axis. A reasonable number of genes, common between previous puberty studies and this study on early reproductive traits, corroborates the proposed molecular mechanisms. This study identified the polymorphism associated with early reproductive traits, and candidate genes that provided a visualization of the proposed mechanisms, coordinating the hypothalamic, pituitary, and ovarian functions for reproductive performance in Brahman cattle.

Expression of steroidogenic enzymes and TGFβ superfamily members in follicular cells of prepubertal gilts with distinct endocrine profiles

Regulation of the transforming growth factor beta (TGFβ) superfamily by gonadotrophins in swine follicular cells is not fully understood. This study evaluated the expression of steroidogenic enzymes and members of the TGFβ superfamily in prepubertal gilts allocated to three treatments: 1200 IU eCG at D -3 (eCG); 1200 IU eCG at D -6 plus 500 IU hCG at D -3 (eCG + hCG); and the control, composed of untreated gilts. Blood samples and ovaries were collected at slaughter (D0) and follicular cells were recovered thereafter. Relative gene expression was determined by real-time PCR. Serum progesterone levels were greater in the eCG + hCG group compared with the other groups (P < 0.01). No differences were observed in the expression of BMP15, BMPR1A, BMPR2, FSHR, GDF9, LHCGR and TGFBR1 (P > 0.05). Gilts from the eCG group presented numerically greater mean expression of CYP11A1 mRNA than in the control group that approached statistical significance (P = 0.08) and greater expression of CYP19A1 than in both the eCG and the control groups (P < 0.05). Expression of BMPR1B was lower in the eCG + hCG treatment group compared with the control (P < 0.05). In conclusion, eCG treatment increased the relative expression of steroidogenic enzymes, whereas treatment with eCG + hCG increased serum progesterone levels. Although most of the evaluated TGFβ members were not regulated after gonadotrophin treatment, the downregulation of BMPR1B observed after treatment with eCG + hCG and suggests a role in luteinization regulation.

Effect of nutritional stress on physiological parameters and seminal attributes of native-crossbred ram in semi-arid tropics

A prolific three-breed (Malpura, Patanwadi, and Garole) cross Avishaan sheep has been developed in the semi-arid zone to improve farmer's income. Nutritional scarcity is a major limitation in animal husbandry during the dearth period of semi-arid tropics. Therefore, before the inaugural launch of the breed into the field, a study was designed to evaluate the effect of nutritional stress on physiological parameters and seminal attributes of native-crossbred rams in semi-arid tropics. Thus, 16 native adapted (Malpura) and 16 native-crossbred rams were equally distributed into four groups, namely, native control (G1), native nutritional stress (G2), native-crossbred control (G3), and native-crossbred nutritional stress (G4). Both the control groups (G1 and G3) were kept on their maintenance requirement as per their body weight, whereas the nutritional stress groups (G2 and G4) were provided 30% less than their maintenance requirement. The body weight of G4 decline (P<0.05) as compared to their initial weight. The plasma glucose level of G2 and G4 reduced (P<0.05) in comparison with G1 and G3, respectively. The total motile sperm percentage, rapid motile sperm percentage, and sperm viability decrease significantly (P<0.05) within the acceptable limit in native-crossbred rams (G4) under nutritional scarcity. However, the similar blood biochemical along with acceptable seminal attributes of all the rams reflected that native-crossbred rams can cope with the nutritional scarcity in semi-arid tropics and have the potential to contribute to the sustainable small ruminant production system for livelihood security in this region.

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.

Stochastic nonlinear model for somatic cell population dynamics during ovarian follicle activation

In mammals, female germ cells are sheltered within somatic structures called ovarian follicles, which remain in a quiescent state until they get activated, all along reproductive life. We investigate the sequence of somatic cell events occurring just after follicle activation, starting by the awakening of precursor somatic cells, and their transformation into proliferative cells. We introduce a nonlinear stochastic model accounting for the joint dynamics of the two cell types, and allowing us to investigate the potential impact of a feedback from proliferative cells onto precursor cells. To tackle the key issue of whether cell proliferation is concomitant or posterior to cell awakening, we assess both the time needed for all precursor cells to awake, and the corresponding increase in the total cell number with respect to the initial cell number. Using the probabilistic theory of first passage times, we design a numerical scheme based on a rigorous finite state projection and coupling techniques to compute the mean extinction time and the cell number at extinction time. We find that the feedback term clearly lowers the number of proliferative cells at the extinction time. We calibrate the model parameters using an exact likelihood approach. We carry out a comprehensive comparison between the initial model and a series of submodels, which helps to select the critical cell events taking place during activation, and suggests that awakening is prominent over proliferation.

Frequency of BMP15 and GDF9 mutations increasing litter size and their phenotypic effects in Olkuska sheep population

Two mutations affecting the ovulation rate and litter size are segregating in Olkuska sheep population, FecXO in the BMP15 gene, and the G7 site mutation in GDF9 gene. Homozygous carriers of both mutations are hyperprolific, contrary to the sterility observed in homozygous carriers of most other BMP15 and GDF9 mutations. The objective of this study was to assess frequency and phenotypic effects of both mutations. Blood samples were obtained from 740 individuals, 111 rams and 629 ewes, out of which 91 rams and 561 ewes were successfully genotyped for the BMP15 and GDF9 loci. The reproductive performance included a number of lambs born/born alive and a number of lambs reared until 60 days of age, and for a subset of ewe ovulation rates. The study proved a high frequency of the FecXO mutation in two flocks that have been selected for many years for increased litter size (0.7–1.0 in breeding ewes and rams respectively), and a moderate frequency in another 19 private flocks (0.4–0.5). The frequency of the GDF9 mutation was low, with only 50 sheep out of 312 genotyped being carriers of the GDF9/G7 mutation, including three homozygous carriers. The FecXO mutation in the BMP15 had a significant effect on both litter size and the ovulation rate. The single copy in heterozygous carriers increased litter size by 0.255 (0.063), while the effect of two copies in homozygous genotypes was +0.874 (0.081) lambs born. Due to the low frequency of the GDF9 mutation, it can only be preliminarily concluded that litter size has been increased in double carriers of both the BM15 and GDF9 mutation, which may suggest their additive interaction. The positive effect of both mutations supports their direct use in selection programmes.

A nonsense mutation of bone morphogenetic protein-15 (BMP15) causes both infertility and increased litter size in pigs

Background

Atypical external genitalia are often a sign of reproductive organ pathologies and infertility with both environmental or genetic causes, including karyotypic abnormalities. Genome-wide association studies (GWAS) provide a means for identifying chromosomal regions harboring deleterious DNA-variants causing such phenotypes. We performed a GWAS to unravel the causes of incidental cases of atypically small vulvae in German Landrace gilts.

Results

A case-control GWAS involving Illumina porcine SNP60 BeadChip-called genotypes of 17 gilts with atypically small vulvae and 1818 control animals (fertile German Landrace sows) identified a significantly associated region on the X-chromosome (P = 8.81 × 10^− 43). Inspection of whole-genome sequencing data in the critical area allowed us to pinpoint a likely causal variant in the form of a nonsense mutation of bone morphogenetic protein-15 (BMP15; Sscrofa11.1_X:g.44618787C>T, BMP15:p.R212X). The mutant allele occurs at a frequency of 6.2% in the German Landrace breeding population. Homozygous gilts exhibit underdeveloped, most likely not functional ovaries and are not fertile. Male carriers do not seem to manifest defects. Heterozygous sows produce 0.41±0.02 (P=4.5 × 10^-83) piglets more than wildtype animals. However, the mutant allele’s positive effect on litter size accompanies a negative impact on lean meat growth.

Conclusion

Our results provide an example for the power of GWAS in identifying the genetic causes of a fuzzy phenotype and add to the list of natural deleterious BMP15 mutations that affect fertility in a dosage-dependent manner, the first time in a poly-ovulatory species. We advise eradicating the mutant allele from the German Landrace breeding population since the adverse effects on the lean meat growth outweigh the larger litter size in heterozygous sows.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07343-x.

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.

Modulation of granulosa cell function via CRISPR-Cas fuelled editing of BMPR-IB gene in goats (Capra hircus)

BMPs are multifunctional growth factors implicated in regulating the ovarian function as key intra-ovarian factors. Biological effects of BMPs are mediated through binding with membrane bound receptors like BMPR-IB and initiating downstream Smad signaling pathway. FecB mutation, regarded as a loss of function mutation in the BMPR-IB gene was identified in certain sheep breeds having high fecundity. Similar type of fecundity genes in goats have not been discovered so far. Hence, the current study was designed to investigate the effects of BMPR-IB gene modulation on granulosa cell function in goats. The BMPR-IB gene was knocked out using CRISPR-Cas technology in granulosa cells and cultured in vitro with BMP-4 stimulation for three different durations In addition, the FecB mutation was introduced in the BMPR-IB gene applying Easi-CRISPR followed by BMP-4/7 stimulation for 72 h. Steroidogenesis and cell viability were studied to explore the granulosa cell function on BMPR-IB gene modulation. BMPRs were found to be expressed stage specifically in granulosa cells of goats. Higher transcriptional abundance of R-Smads, LHR and FSHR indicating sensitisation of Smad signaling and increased gonadotropin sensitivity along with a significant reduction in the cell proliferation and viability was observed in granulosa cells upon BMPR-IB modulation. The inhibitory action of BMP-4/7 on P4 secretion was abolished in both KO and KI cells. Altogether, the study has revealed an altered Smad signaling, steroidogenesis and cell viability upon modulation of BMPR-IB gene in granulosa cells similar to that are documented in sheep breeds carrying the FecB mutation.

Spatiotemporal expression of activin receptor-like kinase-5 and bone morphogenetic protein receptor type II in the ovary of shortfinned eel, Anguilla australis

In the eel ovary, the expression of growth differentiation factor-9 (Gdf9) appears to be largely confined to the germ cell in early stages of oogenesis. However, both the target tissue and the function of Gdf9 in fish remain unknown. This study aimed to describe the abundance and localization of activin receptor-like kinase-5 (Alk5) and bone morphogenetic protein receptor type II (Bmpr2), which together mediate the Gdf9 signal, in the ovary of a basal teleost, the shortfinned eel, Anguilla australis, during early folliculogenesis. The cDNA encoding eel alk5 and bmpr2 genes were cloned, characterized and the transcript abundances of these receptors quantified by quantitative real-time PCR. Ovarian transcript abundance for both receptors, along with that of gdf9 and of its paralogue bmp15, increased from the previtellogenic to early vitellogenic stage. Localization of receptor mRNAs by in situ hybridization revealed that these receptors are located in the somatic cells surrounding the oocyte. Furthermore, tissue distribution analysis showed that the expression of alk5 and bmpr2 were highest in ovary and thyroid, respectively. Unexpectedly, however, bmpr2 mRNA levels were lower in the ovary than in any of the other 17 tissues examined, and indeed, lower than ovarian gdf9 transcript abundance. These findings, together with the ovarian expression pattern of Gdf9, suggest that Gdf9, and conceivably, Bmp15, from the oocyte can signal through receptors that are located on the somatic cells surrounding the oocyte; this, in turn, facilitates elucidation of the function of these growth factors during oogenesis in teleost fish.

Novel indel variations of the sheep FecB gene and their effects on litter size

Various studies had shown that the FecB gene (also known as BMPRIB, BMPR1B or ALK-6) was the major gene influencing sheep litter size, for which its SNPs' variations were reportedly linked. Yet, surprisingly, there was no published information on the insertion/deletion (indel) variation of this gene. Herein, using a population of Chinese Australian White sheep (n = 932), we identified five novel indels in the different introns of the FecB gene. Among them, the 12-bp indel was distinguished as a splicing region variation that was completely linked to the 17-bp indel. Analysis of variance revealed that only the 10-bp indel was significantly associated with sheep litter size (P = 0.010), for which the deletion/deletion (DD) genotype was the harmful one for fecundity. Nevertheless, the combined genotypes of these five indels were significantly (P = 0.033) correlated with the litter size. These findings would provide fresh insight into developing a sounder basis to accelerate molecular breeding in sheep via DNA markers in a marker-assisted selection strategy.

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.

Intrafollicular growth differentiation factor 9: bone morphogenetic 15 ratio determines litter size in mammals†

We previously showed that rat, pig, sheep, and red deer oocytes express species-specific ratios of GDF9: BMP15 mRNA (3.7, 0.5, 1.26, and 0.1, respectively), and with the exception of the pig, they are directly correlated to litter size. The purpose of this study was to determine the alternative mechanism that enables pig oocytes to secrete low ratios whilst maintaining a large litter size. Herein, we performed same- and cross-species coincubations of oocytes with granulosa cells (GCs) of rat, pig, sheep, and red deer to compare the proliferation rate, mRNA expression levels of growth factor receptors, and downstream signalling pathways in GCs. A decreased proliferation rate, lower Bmpr1b and Bmpr2 mRNA expression levels, and higher SMAD1/5/8 protein levels were exhibited in rat GCs cocultured with red deer oocytes, compared to all other species. Pig GCs unequivocally expressed GDF9 mRNA, suggesting that, similar to rat GCs, the proliferation of pig GCs is regulated mainly by GDF9, despite lower intraoocyte expression of GDF9 mRNA. In support, a higher basal proliferation, and their ability to proliferate readily when coincubated with red deer oocytes, was observed in pig GCs. In contrast, red deer GC proliferation is likely to be mainly regulated by BMP15 in vivo with only red deer oocytes capable of altering SMAD1/5/8 and pSMAD2/3 levels, while both GDF9 and BMP15 appear important for sheep GC proliferation. In summary, this study strengthens our hypothesis that the ratio of GDF9: BMP15 in the intrafollicular milieu is directly correlated with litter size, and that the GCs of each species have evolved to respond to these unique ratios.

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.