A missense mutation in growth differentiation factor 9 (GDF9) is strongly associated with litter size in sheep
© Våge et al.; licensee BioMed Central Ltd. 2013
Received: 29 June 2012
Accepted: 28 December 2012
Published: 2 January 2013
A genome wide association study for litter size in Norwegian White Sheep (NWS) was conducted using the recently developed ovine 50K SNP chip from Illumina. After genotyping 378 progeny tested artificial insemination (AI) rams, a GWAS analysis was performed on estimated breeding values (EBVs) for litter size.
A QTL-region was identified on sheep chromosome 5, close to the growth differentiation factor 9 (GDF9), which is known to be a strong candidate gene for increased ovulation rate/litter size. Sequencing of the GDF9 coding region in the most extreme sires (high and low BLUP values) revealed a single nucleotide polymorphism (c.1111G>A), responsible for a Val→Met substitution at position 371 (V371M). This polymorphism has previously been identified in Belclare and Cambridge sheep, but was not found to be associated with fertility. In our NWS-population the c.1111G>A SNP showed stronger association with litter size than any other single SNP on the Illumina 50K ovine SNP chip. Based on the estimated breeding values, daughters of AI rams homozygous for c.1111A will produce minimum 0.46 - 0.57 additional lambs compared to daughters of wild-type rams.
We have identified a missense mutation in the bioactive part of the GDF9 protein that shows strong association with litter size in NWS. Based on the NWS breeding history and the marked increase in the c.1111A allele frequency in the AI ram population since 1983, we hypothesize that c.1111A allele originate from Finnish landrace imported to Norway around 1970. Because of the widespread use of Finnish landrace and the fact that the ewes homozygous for the c.1111A allele are reported to be fertile, we expect the commercial impact of this mutation to be high.
Litter size is an economically important trait in sheep breeding. To date, polymorphisms in three different genes have been associated with increased ovulation rate/litter size in sheep. These are the growth differentiation factor 9 (GDF9), bone morphogenetic protein 15 (BMP15) and bone morphogenetic protein receptor, type IB (BMPR1B).
The BMP15 gene was first reported to be associated with increased ovulation rate in Hanna (FecXH) and Inverdale (FecXI) sheep . In Hanna sheep a glutamine in position 291 was replaced by a premature stop codon (Q291X), while in the Inverdale sheep a valine residue in position 299 was replaced by aspartic acid (V299D). In the Merino Booroola sheep (FecBB), a glutamine to arginine mutation in position 249 (Q249R) of the BMPR1B gene was found to increase ovulation rate and litter size [2–4]. In Cambridge and Belclare sheep two additional mutations affecting ovulation rate were identified in BMP15, a substitution of glutamine with a premature stop codon in position 239 (Q239X) (FecXG) and a change of serine to isoleucine at position 367 (S367I) (FecXB) . In the same study, a substitution of serine with phenylalanine in position 395 (S395F) (FecGH) of the GDF9 gene was found to be associated with increased ovulation rate.
A common feature of the GDF9 and BMP15 mutations described above is that increased ovulation rate/litter size is observed for individuals being heterozygous for these mutations, while individuals being homozygous are sterile, indicating a dose dependent function of these proteins. Later, additional mutations with a similar effect have been identified in BMP15, a cysteine to tyrosine change in position 323 (C323Y) in Lacaune sheep (FecXL)  and a 17 basepair deletion in the Rasa Aragonesa sheep (FecXR) [7, 8]. Also in GDF9 an additional mutation was found, a substitution of serine with arginine in position 427 (S427R) (FecTT) .
In 2011 Silva et al.,  reported a GDF9 mutation, a substitution of phenylalanine with cysteine at position 345 (F345C) (FecGE). In contrast to previously reported mutations in both GDF9 and BMP15, the F345C mutation in GDF9 did not cause sterility in the homozygous individuals, but rather significantly higher prolificacy compared to heterozygous individuals. GDF9 and BMP15 are known to influence ovulation rate in a dose-responsive manner , indicating that the F345C variant has not completely lost its biological function.
In the present study we have used 378 rams selected for artificial insemination (AI rams) to search for polymorphisms affecting litter size in the Norwegian White Sheep (NWS) breed. All animals were included in a genome wide association study using the Illumina ovine 50K SNP chip. A QTL was detected on chromosome 5 which highlighted a previously detected candidate gene; GDF9. Sequencing of this gene in AI rams with extreme EBVs for litter size revealed a non-synonymous mutation in NWS.
378 AI rams were genotyped using Illumina's 50K ovine SNP array. Association testing was performed by a linear mixed model (GEMMA) , using estimated breeding values (EBVs) for daughter litter size as phenotypes. An initial test showed overlapping results when using EBV1, EBV2 and EBV3 as separate phenotypes (results not shown). We decided to use EBV1 phenotypes only for the association analysis, since this estimate is generated from the largest number of daughters compared to EBV2 and EBV3. Allowing for 5% missing data per SNP and 1% minor allele frequency, 47 986 SNPs were included in analysis.
Five Oar5 SNPs on the Illumina ovine 50K SNP array showing significant association with litter size, together with GDF9 c.1111G>A
Position at Oar v3.1
EBV contrasts between GDF9 c.1111 genotypes
GG vs GA
GA vs AA
GG vs AA
By combining a genome wide association study and a candidate gene approach in AI rams of the Norwegian White Sheep breed, we successfully identified a SNP causing an amino acid change in the mature GDF9 protein. This SNP showed a stronger association to EBVs for litter size than any of the SNPs present at the Illumina 50K SNP array (Figure 1B). This sequence variant has previously been reported to exist in Belclare and Cambride sheep (the G7 polymorphism) , but these authors focused on polymorphisms that caused female sterility in the homozygous state. Animals homozygous for c.1111A were found to be fertile in that study, and no additional attempt to correlate this polymorphism with ovulation rate/litter size was performed. Also, the number of animals carrying this mutation was most likely too low to detect any genotype dependent variation in fertility.
GDF9 is mediating its effect by binding to transforming growth factor, beta receptor 1 (TGFBR1) [14, 15] and bone morphogenetic protein receptor, type II (BMPR2) . As pointed out by Hanrahan et al. , the c.1111G>A polymorphism represents a relatively conservative change in the sense that one nonpolar amino acid (V) is substituted by another (M). However, the side chain of methionine is structurally different from that of valine, so depending on the location relative to the receptor binding region of GDF9, a reduced (but not lost) binding capacity can possibly be explained by this change.
Daughters of rams being homozygous for c.1111A gave birth to 0.46 - 0.57 additional lambs compared to daughters of c.1111G homozygous rams, while daughters of heterozygous rams gave 0.20 - 0.25 additional lambs. These figures can be considered as conservative, since they build on EBVs that are regressed towards the mean to compensate for a non-infinite number of progeny tested offspring. Based on the development of the c.1111A allele frequency in the AI rams over time (Figure 5), we can assume a similar but lagging frequency development in the ewe population, and thus low overall frequency. Basically, it will then be the heterozygous effect that is contained in the contrast of GG with AA rams, while the GA group with 50% heterozygote daughters will be logically intermediate. Also, when using the animal model for calculating EBVs without modeling the allele effect of c.1111A, the sire will mainly provide the allele, while the effect will be shared with the dam and therefore be underestimated when the frequency is low among the ewes.
The EBVs were estimated based on the daughters’ performance in terms of number of lambs born, but only ewes that gave births were included. Since GDF9 is known to influence oocyte maturation and ovulation rate, the ewe’s genotype will be determinative for her fertility. Therefore, this study should be followed up by genotyping a large number of ewes with known reproductive performance to obtain a better estimate of the separate genotypic effects, and to confirm that homozygous ewes are fertile as reported for Cambridge and Belclare sheep .
No single event in the breeding history of NWS has influenced litter size more than the introduction of Finnish landrace. Finnish landrace sheep are well known for their high fertility, and have been crossed to several breeds to investigate this trait [17, 18]. Also in Norway, Finnish landrace was imported in late 1960s and early 1970s to improve fertility [19–21]. The significant phenotypic effect of the GDF9 c.1111A - allele observed in this study could therefore indicate that this allele originate from Finnish landrace. The development in allele frequency in the AI rams population since 1990 also resemble the pattern observed for another imported allele, namely the Texel MSTN c.*1232A - allele . In both cases a new breed has been imported to improve key traits, litter size and meatiness respectively, and after some latency the organized breeding system catches the causal allele and rapidly increase its frequency, first in the AI ram population and subsequently in the whole ewe population.
We have identified a missense mutation in the bioactive part of the GDF9 - protein that is strongly associated with litter size in the NWS population. The observed amino acid change from valine to methionine might be the functional explanation of this trait, although the possibility of mutations outside the coding region in close linkage disequlibrium with c.1111G>A should not be excluded. According to the estimated breeding values, daughters of AI rams homozygous for c.1111A will produce minimum 0.46 - 0.57 additional lambs compared to daughters of rams not having this allele. Based on the breeding history of NWS and the recent increase in the allele frequency of c.1111A in this population, we hypothesize that this allele originates from imports of Finnish landrace. Because of the widespread use of Finnish landrace internationally, and the fact that this mutation is reported to also give fertile ewes in the homozygous state, we expect the commercial impact of this mutation to be high.
A total of 378 AI rams of the Norwegian White Sheep (NWS) breed, born between 1983 and 2008, were genotyped by the Illumina 50K SNP array and used for the GWAS - study. In addition to the 378 AI rams genotyped by the 50K SNP array, 21 AI rams born in 2009 and 4 rams born in 2007 and 2008 were genotyped for the GDF9 c.1111G>A mutation and used for calculating allele frequencies. These additional 25 individuals were not included in the association analysis, nor in the analysis of genotype effects.
Progeny testing and EBV calculations
All rams are progeny tested in so-called ram circles [19, 23]. A ram circle consists of several farms that exchange rams during the breeding season to ensure that their offspring for progeny testing are born in different production environments (flocks). From 1991 the best linear unbiased prediction (BLUP) method and the animal model replaced the selection index approach for calculating breeding values. The phenotypic trait “number of lambs born” is referred to as litter size in this paper, and includes both live born and stillborn lambs. Separate BLUP EBVs for litter size are estimated, depending on whether daughters are one, two or three years of age. These data are referred to as EBV1, EBV2 and EBV3, respectively, and mean number of daughters per ram were 235, 212 and 148 respectively. Only ewes giving births are included in the EBV-calculation, leaving out sterile individuals. The phenotypic data used for prediction of breeding values was retrieved from the National Sheep Recording System (SRS), comprising all sheep born after 1987 in the ram circles.
Genotyping by the ovineSNP50
DNA was isolated from semen using Qiagen MagAttract DNA extraction kit according to manufacturers protocols (Qiagen, Germany). DNA concentration was determined using PicoGreen reagent (Invitrogen, USA) and quality was assessed by gel electrophoresis using a 1% TAE agarose gel. Genotyping was performed using the OvineSNP50 array from Illumina (Illumina, USA) according to manufacturers recommendations. Raw data was converted to genotypes using the Illumina's Genotyping Module (version 1.9.4) within the Genome Studio Software (version 2011.1). Automatic clustering was performed with a call-threshold of 0.15, samples with call rates below 98% were excluded before performing manual re-clustering. SNP call frequency, minor allele frequency and pedigree error frequency were used as criteria to sort SNPs and facilitate manual adjustments. Of the 54 241 SNP assays present on the array, 79 were regarded as failing, with another 4 155 being monomorphic in our study.
Estimated breeding values (EBVs) for litter size (measured as number-of-lambs-born when daughters are one year of age, i.e. EBV1) were used as phenotypes of the 378 genotyped AI rams included in the study. A linear mixed-model algorithm (GEMMA), that calculates exact values of standard test statistics, was used for the association analysis . A standard relatedness matrix was estimated from SNP genotypes to account for population stratification and sample structure by this software. Only informative SNPs (> 95% data and > 1% MAF) were included in the study (n = 47 986). The significance of the associations was evaluated with likelihood-ratio test, and a conservative threshold for significance of p < 10-6 was applied using the Bonferroni correction for multiple testing, since (0.05/47986) = 1.04 × 10-6. The distribution of obtained versus expected genome-wide p-values were visualised by QQ-plots. The Manhattan plot and QQ-plots were made by a R-code provided at: http://gettinggeneticsdone.blogspot.no/2011/04/annotated-manhattan-plots-and-qq-plots.html.
PCR amplification, cloning and sequencing of GDF9
Primers used for amplification, sequencing and genotyping of the ovine GDF9 gene
1744 - 1763a
2187 - 2206a
3205 - 3224a
4283 - 4302a
3863 - 3882a
4443 - 4462a
3983 - 4002a
4039 - 4058a
3996 - 4014a
Genotyping of the c.1111G>A polymorphism in ovine GDF9
The Sequenon massARRAY platform (SEQUENOM, San Diego, USA) was used for genotyping the c.1111G>A SNP according to manufacturers recommendations. Amplification primers (GDF9F and GDF9R) and the extension primer (GDF9E) are shown in Table 3. All animals included in the GWAS - study were genotyped by this assay, including those used in sequencing.
To visualise the amino acid conservation across species at the sheep GDF9 371 position we used the COBALT program available at NCBI to align the following GDF9 protein sequences: Sheep (NP_001136360.1), Cattle (NP_777106.1), Pig (NP_001001909.1), Cat (NP_001159372.1), Human (NP_005251.1), Mouse (NP_032136.2), Chicken (NP_996871.2) and Zebra fish (NP_001012383.1).
where EBVi kl = the estimated breeding value for number of lambs born of ewes at one, two and three years of age, for a ram born in year k, with genotype l; μ = overall mean of the trait; BY k = the fixed effect of the kth birth year of ram (1994-2008); G l = the fixed effect of genotype l (GG, GA and AA, respectively) and e kl is a random error term ~ N(0, Iσ e 2), where I is the identity matrix and σ e 2 is the residual variance. The birth year effect was included to account for significant (p < 0.001) genetic change over time, caused by other genes than the SNP in question.
The project has received financial support from the Norwegian Research Council (190217), the Agricultural Agreement Research Fund, the Foundation for Research Levy on Agricultural Products, Animalia - Meat and Poultry Research Centre and the Norwegian Association of Sheep and Goat Breeders (NSG). We will thank Thor Blichfeldt from NSG for stimulating and useful discussions regarding the manuscript.
- Galloway SM, McNatty KP, Cambridge LM, Laitinen MP, Juengel JL, Jokiranta TS, McLaren RJ, Luiro K, Dodds KG, Montgomery GW, et al: Mutations in an oocyte-derived growth factor gene (BMP15) cause increased ovulation rate and infertility in a dosage-sensitive manner. Nat Genet. 2000, 25 (3): 279-283. 10.1038/77033.View ArticlePubMed
- Wilson T, Wu XY, Juengel JL, Ross IK, Lumsden JM, Lord EA, Dodds KG, Walling GA, McEwan JC, O'Connell AR, et al: 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. Biol Reprod. 2001, 64 (4): 1225-1235. 10.1095/biolreprod64.4.1225.View ArticlePubMed
- Mulsant P, Lecerf F, Fabre S, Schibler L, Monget P, Lanneluc I, Pisselet C, Riquet J, Monniaux D, Callebaut I, et al: Mutation in bone morphogenetic protein receptor-IB is associated with increased ovulation rate in Booroola Merino ewes. Proc Natl Acad Sci U S A. 2001, 98 (9): 5104-5109. 10.1073/pnas.091577598.PubMed CentralView ArticlePubMed
- Souza CJ, MacDougall C, Campbell BK, McNeilly AS, Baird DT: The Booroola (FecB) phenotype is associated with a mutation in the bone morphogenetic receptor type 1 B (BMPR1B) gene. J Endocrinol. 2001, 169 (2): R1-6. 10.1677/joe.0.169R001.View ArticlePubMed
- Hanrahan JP, Gregan SM, Mulsant P, Mullen M, Davis GH, Powell R, Galloway SM: Mutations in the genes for oocyte-derived growth factors GDF9 and BMP15 are associated with both increased ovulation rate and sterility in Cambridge and Belclare sheep (Ovis aries). Biol Reprod. 2004, 70 (4): 900-909.View ArticlePubMed
- Bodin L, Di Pasquale E, Fabre S, Bontoux M, Monget P, Persani L, Mulsant P: A novel mutation in the bone morphogenetic protein 15 gene causing defective protein secretion is associated with both increased ovulation rate and sterility in Lacaune sheep. Endocrinology. 2007, 148 (1): 393-400.View ArticlePubMed
- Martinez-Royo A, Jurado JJ, Smulders JP, Marti JI, Alabart JL, Roche A, Fantova E, Bodin L, Mulsant P, Serrano M, et al: A deletion in the bone morphogenetic protein 15 gene causes sterility and increased prolificacy in Rasa Aragonesa sheep. Anim Genet. 2008, 39 (3): 294-297. 10.1111/j.1365-2052.2008.01707.x.View ArticlePubMed
- Monteagudo LV, Ponz R, Tejedor MT, Lavina A, Sierra I: A 17 bp deletion in the Bone Morphogenetic Protein 15 (BMP15) gene is associated to increased prolificacy in the Rasa Aragonesa sheep breed. Anim Reprod Sci. 2009, 110 (1–2): 139-146.View ArticlePubMed
- Nicol L, Bishop SC, Pong-Wong R, Bendixen C, Holm LE, Rhind SM, McNeilly AS: Homozygosity for a single base-pair mutation in the oocyte-specific GDF9 gene results in sterility in Thoka sheep. Reproduction. 2009, 138 (6): 921-933. 10.1530/REP-09-0193.View ArticlePubMed
- Silva BD, Castro EA, Souza CJ, Paiva SR, Sartori R, Franco MM, Azevedo HC, Silva TA, Vieira AM, Neves JP, et al: A new polymorphism in the Growth and Differentiation Factor 9 (GDF9) gene is associated with increased ovulation rate and prolificacy in homozygous sheep. Anim Genet. 2011, 42 (1): 89-92. 10.1111/j.1365-2052.2010.02078.x.View ArticlePubMed
- Moore RK, Erickson GF, Shimasaki S: Are BMP-15 and GDF-9 primary determinants of ovulation quota in mammals?. Trends Endocrinol Metab. 2004, 15 (8): 356-361.View ArticlePubMed
- Zhou X, Stephens M: Genome-wide efficient mixed-model analysis for association studies. Nat Genet. 2012, 44 (7): 821-824. 10.1038/ng.2310.PubMed CentralView ArticlePubMed
- McPherron AC, Lee SJ: GDF-3 and GDF-9: two new members of the transforming growth factor-beta superfamily containing a novel pattern of cysteines. J Biol Chem. 1993, 268 (5): 3444-3449.PubMed
- Mazerbourg S, Klein C, Roh J, Kaivo-Oja N, Mottershead DG, Korchynskyi O, Ritvos O, Hsueh AJ: Growth differentiation factor-9 signaling is mediated by the type I receptor, activin receptor-like kinase 5. Mol Endocrinol. 2004, 18 (3): 653-665.View ArticlePubMed
- Kaivo-Oja N, Mottershead DG, Mazerbourg S, Myllymaa S, Duprat S, Gilchrist RB, Groome NP, Hsueh AJ, Ritvos O: Adenoviral gene transfer allows Smad-responsive gene promoter analyses and delineation of type I receptor usage of transforming growth factor-beta family ligands in cultured human granulosa luteal cells. J Clin Endocrinol Metab. 2005, 90 (1): 271-278.View ArticlePubMed
- Vitt UA, Mazerbourg S, Klein C, Hsueh AJ: Bone morphogenetic protein receptor type II is a receptor for growth differentiation factor-9. Biol Reprod. 2002, 67 (2): 473-480. 10.1095/biolreprod67.2.473.View ArticlePubMed
- Donald HP, Read JL: Performance of Finnish Landrace Sheep in Britain. Anim Prod. 1967, 9: 471-10.1017/S0003356100042033.View Article
- Ercanbrack SK, Knight AD: Lifetime (seven years) production of 1/4 and 1/2 Finnish Landrace ewes from Rambouillet, Targhee and Columbia dams under range conditions. J Anim Sci. 1985, 61 (1): 66-77.PubMed
- Eikje LS, Ådnøy T, Klemetsdal G: The Norwegian sheep breeding scheme: description, genetic and phenotypic change. Animal. 2008, 2 (02): 167-176.View ArticlePubMed
- Adnoy T: Selection for Prolificacy in Finnsheep and in Norwegian Sheep. J Agr Sci Finland. 1988, 60 (6): 518-522.
- Våbenø A: Finsk landrase i norsk saueavl. Historikk og betydning for ulike egenskaper hos norsk kvit sau. Bioforsk report. 2012, Bioforsk, Ås, Norway, 1-27. ISBN: 978-82-71-00901-6, 7
- Boman IA, Klemetsdal G, Nafstad O, Blichfeldt T, Vage DI: Selection based on progeny testing induces rapid changes in myostatin allele frequencies - a case study in sheep. J Anim Breed Genet = Zeitschrift fur Tierzuchtung und Zuchtungsbiologie. 2011, 128 (1): 52-55. 10.1111/j.1439-0388.2010.00879.x.View Article
- Olesen I, Svendsen M, Klemetsdal G, Steine TA: Application of a multiple-trait animal-model for genetic evaluation of maternal and lamb traits in Norwegian sheep. Animal Science. 1995, 60: 457-469. 10.1017/S1357729800013333.View Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.