Identification of a line of sheep carrying a putative autosomal gene increasing ovulation rate in sheep that does not appear to interact with mutations in the transforming growth factor beta superfamily

Genetic regulation of ovulation rate and multiple births

Ovulation rate in many mammalian species is controlled to regulate the numbers of offspring and maximise reproductive success. Pathways that regulate ovulation rate still respond to genetic and environmental factors and show considerable variation within and between species. Genetic segregation, positional cloning, and association studies have discovered numerous mutations and genetic risk factors that contribute to this variation. Notable among the discoveries has been the role of mutations in bone morphogenetic protein 15 (BMP15 ), growth differentiation factor 9 (GDF9 ) and bone morphogenetic protein receptor type 1B (BMPR1B ) from the intra-ovarian signalling pathway contributing to the evidence that signalling from the oocyte is the key driver in follicle regulation rather than circulating gonadotrophin concentrations. Multiple variants in different domains of BMP15 and GDF9 result in partial or complete loss of function of the proteins providing insights into their functional roles and differential regulation contributing to species differences in ovulation rate. Early success encouraged many more studies in prolific strains of sheep, cattle and goats providing a valuable catalogue of genetic variants of large effect increasing ovulation rate and litter size. More recently, genetic association studies are beginning to identify genetic risk factors with smaller effects. Most genes implicated are from pathways with defined roles in regulation of the ovarian function. However, some genomic regions suggest regulation by novel genes. Continuing genetic and related functional studies will add further to our understanding of the detailed regulation of ovulation rate and litter size with implications for health and animal production systems.

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.

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

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

Trio, a novel high fecundity allele: I. Transcriptome analysis of granulosa cells from carriers and noncarriers of a major gene for bovine ovulation rate

A major gene for bovine ovulation rate has been mapped to a 1.2 Mb region of chromosome 10. Screening of coding regions of positional candidate genes within this region failed to reveal a causative polymorphism, leading to the hypothesis that the phenotype results from differences in candidate gene expression rather than alteration of gene structure. This study tested differences in expression of positional candidate genes in granulosa cells between carriers and noncarriers of the high fecundity allele, as well as characterizing differences in the transcriptomic profile between genotypes. Five carriers and five noncarriers, female descendants of "Trio," a carrier of the high fecundity allele were initially used in an RNA-seq analysis of gene expression. Four of ten samples were contaminated with theca cells, so that six samples were used in the final analysis (three of each genotype). Of 14 973 genes expressed, 143 were differentially expressed (false discovery rate P < 0.05) in carriers versus noncarriers. Among the positional candidate genes, SMAD6 was 6.6-fold overexpressed in the carriers compared to noncarriers (P < 5 × 10-5). This result was replicated in an independent group of 12 females (7 carriers and 5 noncarriers) using quantitative real-time PCR; SMAD6 was 9.3-fold overexpressed in carriers versus noncarriers (P = 1.17 × 10-6). Association of overexpression of SMAD6, an inhibitor of the BMP/SMAD signaling pathway, with high ovulation rate corresponds well with disabling mutations in ligands (BMP15 and GDF9) and a receptor (BMPR1B) of this pathway that cause increased ovulation rate in sheep.

Variant <i>GDF9</i> mRNA is likely not the main cause of larger litter size in Iranian Lori-Bakhtyari, Shal, Ghezel, and Afshari sheep breeds

Abstract. This study was carried out to screen the GDF9 gene and evaluate the polymorphism effect on litter size of four Iranian sheep breeds using the PCR-RFLP and PCR-SSCP methods. First, sequencing of the GDF9 gene in 16 twin-birth, 4 triplet-birth, and 2 infertile ewes showed that, in addition to G2, G3, G4, G5, and G6 mutations that have been previously reported in other breeds, a new G0 mutation, called C25T, exists in the GDF9 sequence of 1 out of 22 ewes and causes L9F substitution in the signal peptide region. None of the triplet-birth or infertile ewes carried G1, G4, G7, FecGE, G8, or FecGT mutations. In the second experiment, a large dataset was used: 605 individuals including 496 ewes (145 Afshari, 54 Shal, 126 Ghezel, and 171 Lori-Bakhtyari sheep), and 109 rams (26 Afshari, 23 Shal, 10 Ghezel, and 50 Lori-Bakhtyari sheep. There were no sheep carrying the G7, G8, or Thoka mutations. Among all 109 rams that were used in this study, none of them were homozygous for the G1 mutation. Moreover, abundance of heterozygote rams (G1/G+) varied from 0.0 (Afshari) to 28.6 % (Lori-Bakhtyari and Ghezel). The highest and the lowest frequencies of the G4 mutation were 30.6 and 3.0 % in Shal and Afshari breeds, respectively. Moreover, G4 abundance varied from 0.0 to 42.3 %, from 3.0 to 26.9, and from 3.0 to 30.6 % in rams, ewes, and overall, respectively. There was a significant difference in the abundance of G1 and G4 mutations between breeds. However, neither the G1 nor the G4 mutation was associated with litter size in Afshari, Ghezel, Lori-Bakhtyari, or Shal sheep breeds. In conclusion, the results of this study showed that GDF9 G1 and G4 mutations are not the reason for higher litter size in Iranian sheep. Moreover, the GDF9 G0 and G6 mutations do not cause triplet births or infertility in Iranian ewes. Therefore, it is unlikely that variant GDF9 mRNA induces larger litter size or infertility in Iranian ewes.

A Major Gene for Bovine Ovulation Rate

Half-sib daughters sired by a bull believed to be a carrier of a major gene for high ovulation rate were evaluated for ovulation rate and genotyped in an effort to both test the hypothesis of segregation of a major gene and to map the gene’s location. A total of 131 daughters were produced over four consecutive years at a University of Wisconsin-Madison research farm. All were evaluated for ovulation rate over an average of four estrous cycles using transrectal ultrasonography. The sire and all daughters were genotyped using a 3K SNP chip and the genotype and phenotype data were used in a linkage analysis. Subsequently, daughters recombinant within the QTL region and the sire were genotyped successively with 50K and 777K SNP chips to refine the location of the causative polymorphism. Positional candidate genes within the fine-mapped region were examined for polymorphism by Sanger sequencing of PCR amplicons encompassing coding and 5’ and 3’ flanking regions of the genes. Sire DNA was used as template in the PCR reactions. Strong evidence of a major gene for ovulation rate was observed (p<1x10^-28) with the gene localized to bovine chromosome 10. Fine-mapping subsequently reduced the location to a 1.2 Mb region between 13.6 and 14.8 Mb on chromosome 10. The location identified does not correspond to that for any previously identified major gene for ovulation rate. This region contains three candidate genes, SMAD3, SMAD6 and IQCH. While candidate gene screening failed to identify the causative polymorphism, three polymorphisms were identified that can be used as a haplotype to track inheritance of the high ovulation rate allele in descendants of the carrier sire.

Mutations in the leptin receptor gene associated with delayed onset of puberty are also associated with decreased ovulation and lambing rates in prolific Davisdale sheep

The aim of this study was to determine if single nucleotide polymorphisms (SNPs) in the leptin receptor (LEPR) gene associated with delayed onset of puberty are associated with changes in other reproductive traits in adult ewes. The ovulation rate of ewes homozygous for the SNPs was ~15% lower (PPLEPR SNPs than their wild-type or heterozygous contemporaries. Partial failure of multiple ovulations was also increased (PLEPR had on average 0.2 fewer lambs at mid-pregnancy and at birth compared with the wild-type or heterozygous ewes (PLEPR were strongly associated with poorer reproductive performance in Davisdale ewes, which is likely to be linked to both a reduced number of ova available for fertilisation and an increased number of ewes failing to become pregnant. Increased partial failure of multiple ovulations in ewes with high ovulation rates (i.e. 3 or greater) may also contribute to the poor reproductive performance.

Using sheep lines with mutations in single genes to better understand ovarian function

Livestock populations have been subjected to strong selection pressure to improve reproductive success, and this has led to the identification of lines of animals with increased fecundity. These animals provide a rich biological resource for discovery of genes and regulatory mechanisms that underpin improved reproductive success. To date, three genes, all related to the transforming growth factor β pathway, have been identified as having mutations that lead to alterations in ovulation in sheep. In addition, several other sheep lines have been identified with putative mutations in single genes with major effects on ovulation rate. This review is focused on the identification of the mutations affecting ovulation rate and how these discoveries have provided new insights into control of ovarian function.

Investigation of prolific sheep from UK and Ireland for evidence on origin of the mutations in BMP15 (FecX(G), FecX(B)) and GDF9 (FecG(H)) in Belclare and Cambridge sheep

This paper concerns the likely origin of three mutations with large effects on ovulation rate identified in the Belclare and Cambridge sheep breeds; two in the BMP15 gene (FecX^G and FecX^B) and the third (FecG^H) in GDF9. All three mutations segregate in Belclare sheep while one, FecX^B, has not been found in the Cambridge. Both Belclare and Cambridge breeds are relatively recently developed composites that have common ancestry through the use of genetic material from the Finnish Landrace and Lleyn breeds. The development of both composites also involved major contributions from exceptionally prolific ewes screened from flocks in Ireland (Belclare) and Britain (Cambridge) during the 1960s. The objective of the current study was to establish the likely origin of the mutations (FecX^G, FecX^B and FecG^H) through analysis of DNA from Finnish Landrace and Lleyn sheep, and Galway and Texel breeds which contributed to the development of the Belclare breed. Ewes with exceptionally high prolificacy (hyper-prolific ewes) in current flocks on Irish farms were identified to simulate the screening of ewes from Irish flocks in the 1960s. DNA was obtained from: prolific ewes in extant flocks of Lleyn sheep (n = 44) on the Lleyn peninsula in Wales; hyper-prolific ewes (n = 41); prolific Galway (n = 41) ewes; Finnish Landrace (n = 124) and Texel (n = 19) ewes. The FecX^G mutation was identified in Lleyn but not in Finnish Landrace, Galway or Texel sheep; FecX^B was only found among the hyper-prolific ewes. The FecG^H mutation was identified in the sample of Lleyn sheep. It was concluded from these findings that the Lleyn breed was the most likely source of the FecX^G and FecG^H mutations in Belclare and Cambridge sheep and that the FecX^B mutation came from the High Fertility line that was developed using prolific ewes selected from commercial flocks in Ireland in the 1960′s and subsequently used in the genesis of the Belclare.

Genetic control of multiple births in low ovulating mammalian species

In mammals, litter size is a highly variable trait. Some species such as humans or cattle are monotocous, with one or sometimes two newborns per birth, whereas others, the polytocous species such as mice or pigs, are highly prolific and often produce a dozen newborns at each farrowing. In monotocous species, however, two or three newborns per birth may sometime be unwanted. In more polytocous species such as sheep or pigs, litter size is studied in order to increase livestock prolificacy. By contrast, twinning rates in humans or cattle may increase birth difficulties and health problems in the newborns. In this context, the aim of our review was to provide a clearer understanding of the genetic and physiological factors that control multiple births in low-ovulating mammalian species, with particular focus on three species: sheep, cattle, and humans, where knowledge of the ovulation rate in one may enlighten findings in the others. This article therefore reviews the phenotypic and genetic variability observed with respect to ovulation and twinning rates. It then presents the QTL and major genes that have been identified in each species. Finally, we draw a picture of the diversity of the physiological mechanisms underlying multiple ovulation. Although several major genes have been discovered in sheep, QTL detection methods in humans or cattle have suggested that the determinism of litter size is complex and probably involves several genes in order to explain variations in the number of ovulations.