Articles | Volume 66, issue 4
https://doi.org/10.5194/aab-66-401-2023
https://doi.org/10.5194/aab-66-401-2023
Original study
 | 
07 Dec 2023
Original study |  | 07 Dec 2023

Detection of genetic variations in the GDF9 and BMP15 genes in Kazakh meat–wool sheep

Makpal Amandykova, Zarina Orazymbetova, Tilek Kapassuly, Altynay Kozhakhmet, Saltanat Khamzina, Kairat Iskakov, and Kairat Dossybayev
Abstract

Kazakh meat–wool sheep are of great interest because of the intrabreed multifetal type's high productivity of 140 %–160 %. Genes encoding growth differentiation factor-9 (GDF9) and bone morphogenetic protein 15 (BMP15) are promising candidates for studying sheep productivity, as they affect fertility in mammals, including sheep. Thus, the purpose of this study was to assess the fertility of the Kazakh meat–wool sheep breed based on GDF9 and BMP15 candidate genes of fecundity for the selection of animals with valuable genotypes. We selected 300 heads of the Kazakh meat–wool sheep breed from two populations for PCR-RFLP (polymerase chain reaction–restriction fragment length polymorphism) analysis, 15 of which were subsequently used for sequencing of exon regions of the GDF9 and BMP15 genes. The sheep populations were tested for G1 and G8 mutations of the GDF9 gene and B2 and B4 mutations of the BMP15 gene. The PCR-RFLP analysis revealed that 59 (19.7 %) of the 300 Kazakh meat–wool breed sheep were heterozygous carriers of the G1 mutation (genotype AG) of the GDF9 gene, and sequencing analysis supported these results. The comparative phylogenetic analysis showed a clear separation of Kazakh meat–wool sheep wild types and carriers of the G1 mutation. This mutation was reported to have a relationship with the animals' litter size in other sheep breeds. For this reason, similar relationships should be investigated in Kazakh meat–wool sheep. However, G8, B2, and B4 mutations were not detected among the studied animal populations, showing that these mutations are not characteristic of the Kazakh meat–wool sheep breed.

1 Introduction

Sheep breeding is one of the leading livestock sectors in the Republic of Kazakhstan. In 2011, after long-term selection and breeding work, scientists and practitioners, sheep breeders in Kazakhstan, tested a multiparous intrabreed type of the Kazakh meat–wool sheep breed. This type was created by crossing Kazakh fine-wool ewes with Finnish Landrace rams. The fertility of ewes is 140 %–160 % (Makhatov et al., 2021).

Sheep fertility is a principal breeding trait, as it largely determines the economic effect of animal husbandry. High precocity, combined with high fertility, can provide a quick return on investment in the industry. Today, four major fecundity genes affecting follicle growth and ovulation frequency are known. Three loci belonging to the TGF-β (transforming growth factor β) family include the Booroola BMPR1B gene (bone morphogenetic protein receptor 1B), BMP15 (bone morphogenetic protein 15), and GDF9 (growth differentiation factor 9) (Monsivais et al., 2017), and the fourth is the B4GALNT2 gene encoding a glycosylation enzyme that does not belong to the BMP (bone morphogenetic protein) family (Ben et al., 2019). These genes regulate the expression and secretion of hormones that affect follicle growth and ovulation frequency (McNatty et al., 2003). The GDF9 and BMP15 genes belong to the TGF-β family, a transforming growth factor, and are principal factors in folliculogenesis regulation and ovulation in sheep (Liu et al., 2019; Chu et al., 2005; Dong et al., 1996; Nilsson and Skinner, 2002). The polymorphism of these genes may significantly impact the productivity of sheep, including their reproductive function and livestock output.

GDF9 gene mutations (FecGH, FecGT, FecGE, FecGF, and FecGV) led to hyperproliferation in heterozygous sheep and sterility in homozygous sheep (Melo et al., 2008; El Fiky et al., 2017; Mullen et al., 2013). Mutations in the FecGE and FecGF of the GDF9 gene affect fecundity traits like ovulation rate and litter size, while mutations in FecGH, FecGT, and FecGV cause increased ovulation rate and litter size in heterozygote ewes and infertility in homozygote carriers (Muhaghegh Dolatabady and Habibizad, 2019). Many studies of the GDF9 gene mutation aim at revealing its association with an increase in litter size. For example, studies in short-tailed sheep have shown that individuals with simultaneous mutations in the GDF9 and BMP15 genes have a higher ovulation frequency than animals carrying only one mutation each (Hanrahan et al., 2004).

BMP15 was the first gene to be associated with fertility. All BMP15 mutations (FecXI, FecXH, FecXB, and FecXG) result in the same phenotype. Heterozygous sheep are highly fertile, while homozygous females are infertile because of inhibition of follicle development at the primary stage (Hanrahan et al., 2004; Galloway et al., 2000; Otsuka et al., 2011). Demars et al. (2013) identified two new mutations in the BMP15 gene associated with a rise in litter size (LS) and ovulation rate (OR). Interestingly, homozygous FecXGr/FecXGr Grivette sheep and homozygous FecXO/FecXO Olkuska sheep exhibited hyperproliferation, which is the opposite of the sterility expressed by all other known BMP15 homozygous mutations (Demars et al., 2013).

Eight point mutations of G1, G2, G3, G4, G5, G6, G7, and G8 have been detected in the GDF9 gene of Belclare and Cambridge sheep breeds (Hanrahan et al., 2004), and only five mutations will lead to deduced amino acid exchanges: G1 – R87H, G4 – E241K, G6 – V332I, G7– V371M, and G8 – S395F (Polley et al., 2010). The first mutation (G1) occurs in exon 1, and others (G4, G6, G7, and G8) occur in exon 2 of the GDF9 gene. As for the BMP15 gene there were identified four polymorphisms across the coding region (B1, B2, B3, and B4). But only two of them led to changes in the amino acid sequence, potentially affecting the fertility and sterility of sheep. Hanrahan et al. (2004) showed that only the G8 mutation in GDF9 and the B2 and B4 mutations in BMP15 were associated with the sterility phenotype. All sheep homozygous for G8 were sterile, all sheep homozygous for B2 or homozygous for B4 were sterile, and sheep heterozygous for both B2 and B4 simultaneously (i.e., one allele with each mutation, B2/B4) were sterile (Hanrahan et al., 2004).

In this paper, we investigate the polymorphisms of the GDF9 and BMP15 genes in the Kazakh meat–wool sheep breed and identify their relationship with the reproductive function and productivity of animals. To the best of our knowledge, studies on mutations in the GDF9 and BMP15 genes have not been determined in Kazakh native sheep breeds. This is especially crucial for the Kazakh meat–wool sheep breed, considering it is multifetal. Thus, the purpose of this work was to assess the genetic polymorphism of the Kazakh meat–wool sheep breed based on the GDF9 and BMP15 candidate genes responsible for fertility, as well as to determine the frequency of various alleles and genotypes for these genes in the studied populations.

2 Materials and methods

2.1 Sample collection

The objects of the study were 2- to 3-year-old sheep of Kazakh meat–wool breed bred in the Kuatzhan (n=75 females) and Dukeev (n=225 (of them 155 females and 70 males)) collective farms in the Zhambyl district of the Almaty region of the Republic of Kazakhstan. The sheep's peripheral blood was used as a biological material. Blood samples were collected by a veterinarian in EDTA vacuum tubes. All animal care and experiments were approved by the Local Ethics Committee of the Institute of Genetics and Physiology SC MSHE RK (19 October 2021, Almaty, Kazakhstan). Further, the blood samples were delivered to the Laboratory of Animal Genetics and Cytogenetics of the RSE Institute of Genetics and Physiology in containers with refrigerant and stored in a freezer (at −25C) until they were used for DNA extraction.

2.2 DNA extraction and quality assessment

DNA isolation was performed using the DNA-sorb-B kit for DNA extraction from clinical materials (AmpliSens, Moscow, Russia) according to the manufacturer's protocol. The determination of quantitative indicators was carried out on the NanoDrop One (ThermoScientific, USA), while the concentration of isolated DNA was 60–100 ng µL−1 on average. DNA qualitative characteristics were checked through agarose gel electrophoresis. Purified DNA was stored in a freezer at −25C for a week before being used for PCR-RFLP (polymerase chain reaction–restriction fragment length polymorphism) analysis.

2.3 PCR-RFLP analysis

PCR-RFLP analysis was performed to check for G1 and G8 mutations for the GDF9 gene and B2 and B4 mutations for the BMP15 gene among the two studied populations of Kazakh meat–wool sheep breed. The PCR-RFLP analysis conditions for the four mutations are shown in Table 1.

Table 1Primer sequences and PCR-RFLP analysis conditions for GDF9 and BMP15 genotyping.

Download Print Version | Download XLSX

The digested fragments were checked using a 4 % agarose gel and visualized with ethidium bromide staining. The gels were scored for the presence or absence of the mutations by iBright Imaging Systems (ThermoScientific, USA). The statistical analysis of the obtained data was conducted using POPGENE 1.32 (https://mybiosoftware.com/popgene-1-32-population-genetic-analysis-2.html, last access: 15 November 2023).

2.4 Sequencing of GDF9 and BMP15 genes and data analysis

Sequencing of exon regions of the GDF9 and BMP15 genes was carried out to build a genetic tree in comparison with sequence data of these genes from other sheep populations and to confirm the existence or absence of the desired mutations in the studied populations. For this purpose, 15 samples selected according to agarose gel visuals were subjected to sequence analysis. First, the GDF9 and BMP15 sheep genes were amplified using PCR with primers designed on the IDT website (https://eu.idtdna.com/pages/tools/primerquest, Integrated DNA technologies, 2023) from published sheep sequences (sheep genomic BMP15 exon 1, AF236078; sheep genomic BMP15 exon 2, AF236079; sheep genomic GDF9 exons 1 and 2, AF078545). The resulting PCR products were sequenced on an ABI 3500XL sequencer (Applied Biosystems) using BigDye™ Terminator v3.1 Cycle Sequencing Kit and the BigDye XTerminator™ Purification Kit. Detailed information on the primers is shown in Table 2.

Table 2Primers used for the amplification of the GDF9 and BMP15 genes.

Download Print Version | Download XLSX

Obtained sequence results belonging to the mutant and wild-type alleles were aligned and compared with the GDF9 sequences of Ovis aries from Sudan, USA, Brazil, Türkiye, Mexico, Norway, China, Iraq, and Egypt (GenBank accession numbers listed in the phylogenetic tree), using BioEdit 7.7 (Alzohairy, 2011). The phylogenetic tree was constructed using the neighbor-joining method on MEGA software (https://www.megasoftware.net/, Molecular Evolutionary Genetics Analisys, 2023).

3 Results and discussion

3.1 PCR-RFLP analysis

The G1 mutation causes an increase in ovulation rate or multiple births in heterozygous ewes. However, for homozygous ewes, it causes sterility because of the arrested follicular development (Hanrahan et al., 2004). Some authors reported that heterozygous ewes with the G1 mutation were fertile and had a large litter size, while wild-type ewes had a smaller litter size (Abdoli et al., 2013). Another study revealed that ewes of the AG genotype have a higher number of lambs (1.88) compared to AA genotypes (1.22) (Gorlov et al., 2018). Paz et al. (2015) reported that the G1 mutation increased the ovulation rate in heterozygous ewes in Chilota sheep. Apart from its effects on ovulation rates, it was reported that lamb weight at birth in heterozygous ewes was higher than in homozygous ewes (Getmantseva et al., 2019). Liandris et al. (2012) revealed that the G1 and G8 mutations of the GDF9 gene were only significantly over-presented in the highly prolific Greek sheep breed Chios.

Regarding the sterility of sheep, Hanrahan et al. (2004) revealed that only the G8 change in GDF9 and the B2 and B4 changes in BMP15 were associated with the sterility phenotype. Thus, homozygous individuals of sheep that do not have any of the three mutations sought are fertile (+/+), sheep heterozygous for one of the three mutations have an increased ovulation rate (G8/+; B2/+ or B4/+), and homozygous carriers of one of the three mutations of the individual are sterile (G8/G8; B2/B2 or B4/B4). However, it is also crucial to consider the relationship of these mutations to each other: individuals that are heterozygous at the same time for mutations B2 (B2/+) and B4 (B4/+) are sterile, and heterozygous sheep for G8 (G8/+) and B2 (B2/+) or G8 (G8/+) and B4 (B4/+) simultaneously have an even higher ovulation rate. Also, studies report that the effect of the GDF9 mutation is significantly higher than that of the BMP15 mutations. However, the GDF9 and BMP15 mutations have an additive effect while appearing in one animal simultaneously, implying that GDF9 and BMP15 are likely to work independently (Hanrahan et al., 2004).

Considering this information, we selected four mutations (G1, G8, B2, and B4) for PCR-RFLP analysis. We analyzed 300 individuals of the Kazakh meat–wool sheep breed from two distinct populations for these mutations. As for the G1 mutation, the nucleotide substitution G to A in the GDF9 exon1 region disrupts the cleavage site of the HhaI restriction enzyme (GCGC to GCAC) at nucleotide 260 of the 462 bp PCR product. Thus, after the treatment of a 462 bp PCR product with the Hha1 enzyme, three fragments of 254, 156, and 52 bp were observed for the wild homozygous GG genotype; four fragments of 410, 254, 156, and 52 bp for the mutant heterozygous GA genotype; and two fragments of 410 and 52 bp for the mutant homozygous AA genotype (Fig. 1).

https://aab.copernicus.org/articles/66/401/2023/aab-66-401-2023-f01

Figure 1Determination of genotypes for the GDF9 gene (mutation G1). M – GeneRuler DNA Ladders (ThermoScientific, USA). 1–13 – digested fragments.

Download

As for the G8 mutation, the 139 bp PCR fragment was observed with two DdeI restriction sites in the AA genotype at 31 and 108 bp (wild type). The results indicated no polymorphism in this exon among the studied sheep of the Kazakh meat–wool breed. The G8 mutation was also not found in other populations of sheep considered to be prolific: Kermani (Khodabakhshzadeh et al., 2016) and Shal sheep in Iran (Ghaffari et al., 2009), Merino sheep in Germany (Chu et al., 2005) and China (Guan et al., 2005), and others (Bai et al., 2007; Sun et al., 2009). Thus, it can be assumed that the G8 mutation is not associated with litter size in Kazakh meat–wool sheep, as in the abovementioned sheep breeds.

Table 3Distribution of genotypes for the GDF9 gene (mutation G1) in two populations of Kazakh meat–wool sheep breed.

Download Print Version | Download XLSX

The following data were used to detect the mutations in the BMP15 gene: processing of a 141 bp fragment of the BMP15 gene with the Hinf1 restriction enzyme results in two fragments of 106 and 35 bp (in the presence of mutation B2), and processing of a 153 bp fragment with Dde1 restriction enzyme forms two fragments of 118 and 35 bp (in the presence of mutation B4). Thus, in the mutation's presence, three fragments should be detected on the electrophoregram: a fragment with the original length of the PCR product and two fragments subjected to restriction enzyme processing. However, the third 25–35 bp fragment cannot be seen visually. The wild type, which does not have the desired mutation, is not subjected to restriction enzyme treatment. Therefore, only one fragment is detected on the electrophoregram, which kept the original size of the PCR product. The study results showed that B2 and B4 single-nucleotide polymorphisms (SNPs) of the BMP15 gene were monomorphic in the Kazakh meat–wool sheep breed populations, which is consistent with the results of the Iran-Black sheep population (Rezaei et al., 2020). Liandris et al. (2012) reported that the B4 mutation of the BMP15 gene was only significantly over-presented in the less prolific Greek sheep breed Karagouniki. Thus, we may suggest that this mutation does not determine the level of fertility in some sheep breeds, including the Kazakh meat–wool breed. This assumption should be tested in further studies aimed at studying the effect of this mutation on the fertility of the Kazakh meat–wool sheep breed by comparing the genotypic data with the phenotypic ones. It is also crucial to understand that the presence or absence of these mutations at the same time is not an indicator of increased fertility or sterility in animals, as in the study of Iranian fat-tailed sheep, where prolificacy genes were polymorphic in the studied population (Abdoli et al., 2013).

Regarding the G1 mutation of the GDF9 gene, 78.7 % of the studied sheep were homozygous for the wild type (GG genotype). Thus, the studied mutation occurred only in 1.6 % of sheep in the homozygous form (G1/G1, genotype AA) and 19.7 % of animals in the heterozygous state (G1/+, genotype GA). There were no differences in these parameters among the two populations (Table 3).

These results are consistent with Kirikçi et al. (2021), where 21 % of the studied Turkish Karayaka sheep were heterozygous GA genotype carriers. However, in our study, the mutant homozygous AA genotype was also found, albeit in a small amount. The homozygous mutant genotype AA was also found among Sudanese desert sheep breeds, and its frequency was equal to 0.06. It was reported that the ewes with heterozygous (GA) and homozygous wild-type (GG) genotypes had more lambs than the homozygous (AA) genotypes (Abdelgadir et al., 2021) compared with the Kazakh meat–wool breed, where this indicator was equal to 0.01. A few number of homozygous individuals for the mutant genotype can be explained by the low frequency (0.12) of the A allele in the Kazakh meat–wool sheep breed, and mutant genotypes could suffer from embryonic death and reproductive defects. The obtained result might also be due to differences in the biological effects of mutations by species (Dinçel et al., 2018; Javanmard et al., 2011).

The data obtained were analyzed for the determination of several statistical indicators: the identified and effective number of alleles, Nei's gene diversity, Shannon's information index, and original measures of Nei's genetic identity and genetic distance (Table 4).

Table 4Summary of the genetic variation statistics for all loci of the G1 mutation of the GDF9 gene.

Na: observed number of alleles. Ne: effective number of alleles. h: Nei's gene diversity. I: Shannon's information index.

Download Print Version | Download XLSX

The average value of the observed number of alleles (Na), representing the actual number of alleles found in the studied populations, was the same (2.0000) for all identified genotypes in both populations. The average effective number of alleles (Ne), describing the number of alleles with the same frequency required to achieve the same expected heterozygosity as in the population under study, varied from 1.0089 to 1.9998, while in the second population it was higher, possibly because of the difference in the studied animals' number among the populations. In terms of gene diversity (h), it can be noted that a significant difference in this indicator was observed for the AA locus, 0.0396 and 0.0089, in the first and second populations, respectively. Shannon's index (I), which measures gene diversity, showed that the genetic diversity of the Dukeev collective farm population was higher than that of the Kuatzhan collective farm population, at 0.3324 and 0.2969, respectively. Both populations were found to be in Hardy–Weinberg disequilibrium (χ2=2.553) because of the heterozygosity deficit. Possible causes of heterozygosity deficits are inbreeding, selection, sub-population structure, and genetic drift.

https://aab.copernicus.org/articles/66/401/2023/aab-66-401-2023-f02

Figure 2(a, b) Identified SNP in the exon 1 of GDF9 gene in Kazakh meat–wool sheep breed.

Download

3.2 Sequencing of GDF9 and BMP15 genes

Sequencing of exon regions of GDF9 in 15 Kazakh meat–wool sheep was performed to confirm the G1 mutation of the GDF9 gene detected by PCR-RFLP analysis. Sequencing results for both alleles were compared to the GDF9 sequencing results for Sudanese sheep (GenBank accession: KY310682.1) to show the nucleotide change (Fig. 2).

Abdelgadir et al. (2021) revealed that the presence of one copy of the G1 mutation of GDF9 gene increased litter size in the studied Sudanese Desert sheep, and this locus may be used as a biomarker for litter size improvement through genotypic selection and allele or gene introgression.

Sequences of exon 1 of the GDF9 gene of different sheep breeds with high identities were downloaded from NCBI in the FASTA format and used for phylogenetic tree construction. Thus, wild individuals (genotype GG) and carriers of the G1 mutation (genotype AG) settled down on separate branches of the phylogenetic tree (Fig. 3).

https://aab.copernicus.org/articles/66/401/2023/aab-66-401-2023-f03

Figure 3Phylogenetic tree based on GDF9 gene of different sheep breeds.

Download

Remarkably, the carriers of the G allele (wild type) were in the same branch of the phylogenetic tree with Sudanese Dubasi, Shugor, and Watish sheep breeds, which had the lowest frequencies of G1 mutant genotypes (Abdelgadir et al., 2021). Also, there were other sheep breeds in which the G1 mutation was not identified, such as Brazilian Santa Inês sheep (Silva et al., 2011) and Norwegian white sheep (Våge et al., 2013) (Muñoz-García et al., 2021; Bodensteiner et al., 1999), whereas the carriers of the A allele (mutant type) are on a close branch with the carriers of this mutation in Turkish sheep breeds (Kirikçi et al., 2021). Egyptian sheep breeds, in which the G1 mutation was also not found (Saleh, 2020), are on a branch of the phylogenetic tree separate from all other sheep breeds. The sequences of the exon2 of the GDF9 gene and the BMP15 (exon1 and exon2) were also used for phylogenetic tree construction to determine the relatedness of the Kazakh meat–wool breed with other world sheep populations according to the studied genes (Figs. S1 and S2 in the Supplement). Thus, it was revealed that the Kazakh meat–wool sheep breed has a close phylogenetic relationship with sheep of Mexican origin in the GDF9 gene (exon 2) and with sheep of Iranian and Brazilian origin in the BMP15 gene.

The present results have important practical implications for the Kazakh meat–wool sheep breed. The application of a marker-assisted selection scheme could be pursued, with introgression of the favored allele(s) showing positive dominance effects on litter size. It could lead to a significant improvement in fecundity, with apparent economic implications for the breeders.

4 Conclusion

Considering the findings from the present study, fertility traits might be improved in the Kazakh meat–wool sheep breed, as evidenced in other sheep breeds. However, it is necessary to perform further studies to provide exact evidence of the effects of the G1 mutation on fertility traits in the Kazakh meat–wool sheep breed. In conclusion, we determined that exon 1 of the GDF9 gene was polymorphic in the Kazakh meat–wool sheep breed and found the G1 mutation (c.260G>A) on this gene. This mutation was reported to have the relationships with the animals' litter size in other sheep breeds. For this reason, similar relationships should be investigated in Kazakh meat–wool sheep.

Data availability

GDF9 and BMP15 gene sequencing data of Kazakh meat–wool sheep used in this study are available from GenBank accession no. PRJNA1001622 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA1001622/, Institute of Genetics and Physiology, 2023).

Supplement

The supplement related to this article is available online at: https://doi.org/10.5194/aab-66-401-2023-supplement.

Author contributions

MA: conceptualization, methodology, data curation, formal analysis, software, project administration, visualization, writing (original draft preparation, review, and editing). ZO: conceptualization, writing (original draft preparation). TK: conceptualization, methodology, software, formal analysis, investigation, data curation. AK: methodology, formal analysis, investigation. SK: methodology. KI: conceptualization, project administration, supervision, validation. KD: conceptualization, methodology, project administration, software, data curation, supervision, validation, writing (review and editing).

Competing interests

The contact author has declared that none of the authors has any competing interests.

Ethical statement

For experimental animals, all protocols were performed in accordance with the “Guide to Animal Experimentation” and approved by the Local Ethics Committee of the Institute of Genetics and Physiology SC MSHE RK (19 October 2021, Almaty, Kazakhstan).

Disclaimer

Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

Acknowledgements

The authors express sincere thanks to the heads of the Kuatzhan and Dukeev collective farms for providing biomaterial for this research work and Amina Omarova for assistance in uploading sequencing data to the GenBank database.

Financial support

This research was supported by the Scientific Committee of Ministry of Science and Higher Education of the Republic of Kazakhstan, grant no. AP09058265, “Evaluation of productivity of the Kazakh meat–wool sheep breed by fecundity genes”.

Review statement

This paper was edited by Joachim Weitzel and reviewed by Bin Tong and two anonymous referees.

References

Abdelgadir, A. Z., Musa, L. M. A., Jawasreh, K. I., Saleem, A. O., El-Hag, F., and Ahmed, M. A.: G1 point mutation in growth differentiation factor 9 gene affects litter size in Sudanese desert sheep, Vet. World, 14, 104–112, https://doi.org/10.14202/vetworld.2021.104-112, 2021. 

Abdoli, R., Zamani, P., Deljou, A., and Rezvan, H.: Association of BMPR-1B and GDF9 genes polymorphisms and secondary protein structure changes with reproduction traits in Mehraban ewes, Gene, 524, 296–303, https://doi.org/10.1016/j.gene.2013.03.133, 2013. 

Alzohairy, A. M.: BioEdit: An important software for molecular biology, Gerf. Bull. Biosci., 2, 60–61, 2011. 

Bai, J., Shi, H. C., Liu, M. J., and Li, F.: Study on GDF9 and BMP15 as candidate genes for fecundity in Cele Black sheep and Duolang sheep, Grass-Feed Livestock, 2, 1–4, 2007. 

Ben, J. S., Ruesche, J., Sarry, J., Woloszyn, F., Lassoued, N., and Fabre, S.: The high prolificacy of D'man sheep is associated with the segregation of the FecLL mutation in the B4GALNT2 gene, Reprod. Domest. Anim., 54, 531–537, https://doi.org/10.1111/rda.13391, 2019. 

Bodensteiner, K. J., Clay, C. M., Moeller, C. L., and Sawyer, H. R.: Molecular cloning of the ovine Growth/Differentiation factor-9 gene and expression of growth/differentiation factor-9 in ovine and bovine ovaries, Biol. Reprod., 60, 381–386, https://doi.org/10.1095/biolreprod60.2.381, 1999. 

Chu, M. X., Sang, L. H., Wang, J. Y., Fang, L., and Ye, S. C.: Study on BMP15 and GDF9 as candidate genes for prolificacy of Small Tail Han sheep, Yi Chuan Xue Bao, 32, 38–45, 2005. 

Demars, J., Fabre, S., Sarry, J., Rossetti, R., Gilbert, H., Persani, L., Tosser-Klopp, G., Mulsant, P., Nowak, Z., Drobik, W., Martyniuk, E., and Bodin, L.: Genome-Wide Association Studies Identify Two Novel BMP15 Mutations Responsible for an Atypical Hyperprolificacy Phenotype in Sheep, PLoS Genet., 9, e1003482, https://doi.org/10.1371/journal.pgen.1003482, 2013. 

Dinçel, D., Ardıçlı, S., Şamlı, H., and Balcı, F.: Genotype frequency of FecXB (Belclare) mutation of BMP15 gene in Chios (Sakiz) sheep, Uludağ University Journal of the Faculty of Veterinary Medicine, 37, 87–91, https://doi.org/10.30782/uluvfd.413857, 2018. 

Dong, J., Albertini, D. F., Nishimori, K., Kumar, T. R., Lu, N., and Matzuk, M. M.: Growth differentiation factor-9 is required during early ovarian folliculogenesis, Nature, 383, 531–535, https://doi.org/10.1038/383531a0, 1996. 

El Fiky, Z. A., Hassan, G. M., and Nassar, M. I.: Genetic polymorphism of growth differentiation factor 9 (GDF9) gene related to fecundity in two Egyptian sheep breeds, J. Assist. Reprod. Gen., 34, 1683–1690, https://doi.org/10.1007/s10815-017-1007-2, 2017. 

Galloway, S. M., McNatty, K. P., Cambridge, L. M., Laitinen, M. P., Juengel, J. L., Jokiranta, T. S., McLaren, R. J., Luiro, K., Dodds, K. G., Montgomery, G. W., Beattie, A. E., Davis, G. H., and Ritvos, O.: Mutations in an oocyte-derived growth factor gene (BMP15) cause increased ovulation rate and infertility in a dosage-sensitive manner, Nat. Genet., 25, 279–283, https://doi.org/10.1038/77033, 2000. 

Getmantseva, L., Bakoev, N., Shirokova, N., Kolosova, M., Bakoev, S., Kolosov, A., Usatov, A. V., Shevtsova, V., and Kolosov, Y. A.: Effect of the GDF9 gene on the weight of lambs at birth, Bulg. J. Agric. Sci., 25, 153–157, 2019. 

Ghaffari, M., Nejati-Javaremi, A., and Rahimi-Mianji, G.: Lack of polymorphism in the oocyte derived growth factor (GDF9) gene in the Shal breed of sheep, S. Afr. J. Anim. Sci., 39, 355–360, https://doi.org/10.4314/sajas.v39i4.51127, 2009. 

Gorlov, I. F., Kolosov, Y. A., Shirokova, N. V., Getmantseva, L. V., Slozhenkina, M. I., Mosolova, N. I., Bakoev, N. F., Leonova, M. A., Kolosov, A. Y., and Zlobina, E. Y.: GDF9 gene polymorphism and its association with litter size in two Russian sheep breeds, Rendiconti Lincei Fisiche e Naturali, 29, 61–66, https://doi.org/10.1007/s12210-017-0659-2, 2018. 

Guan, F., Ai, J. T., Pang, X. S., Liu, S. R., Shi, G. Q., and Yang, L. G.: Detection of the polymorphisms of GDF9 and BMP15 genes in sheep, Life Sci. Res., 9, 184–188, 2005. 

Hanrahan, J. P., Gregan, S. M., Mulsant, P., Mullen, M., Davis, G. H., Powell, R., and Galloway, S. M.: 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., 70, 900–909, https://doi.org/10.1095/biolreprod.103.023093, 2004. 

Institute of Genetics and Physiology: Evaluation of productivity of the Kazakh meat-wool sheep breed by fecundity genes, NCBI, https://www.ncbi.nlm.nih.gov/bioproject/PRJNA1001622/, last access: 6 December 2023. 

Integrated DNA technologies: PrimerQuest™ Tool, https://eu.idtdna.com/pages/tools/primerquest, last access: 25 July 2023. 

Javanmard, A., Azadzadeh, N., and Esmailizadeh, A. K.: Mutations in bone morphogenetic protein 15 and growth differentiation factor 9 genes are associated with increased litter size in fat-tailed sheep breeds, Vet. Res. Commun., 35, 157–167, https://doi.org/10.1007/s11259-011-9467-9, 2011. 

Khodabakhshzadeh, R., Mohammadabadi, M. R., Esmailizadeh, A. K., Shahrebabak, H. M., Bordbar, F., and Namin S. A.: Identification of point mutations in exon 2 of GDF9 gene in Kermani sheep, Pol. J. Vet. Sci., 19, 281–289, https://doi.org/10.1515/pjvs-2016-0035, 2016. 

Kirikçi, K., Çam, M. A., and Mercan, L.: Investigation of G1 (c.260G>A) polymorphism in exon 1 of GDF9 gene in Turkish sheep breed Karayaka, Turk. J. Vet. Anim. Sci., 45, 21, https://doi.org/10.3906/vet-2009-40, 2021. 

Liandris, E., Kominakis, A., Andreadou, M., Kapeoldassi, K., Chadio, S., Tsiligianni, T., Gazouli, M., and Ikonomopoulos, I.: Associations between single nucleotide polymorphisms of GDF9 and BMP15 genes and litter size in two dairy sheep breeds of Greece, Small Ruminant Res., 107, 16–21, https://doi.org/10.1016/j.smallrumres.2012.04.004, 2012. 

Liu, M. N., Zhang, K., and Xu, T. M.: The role of BMP15 and GDF9 in the pathogenesis of primary ovarian insufficiency, Hum. Fertil. (Camb.), 5, 325–332, https://doi.org/10.1080/14647273.2019.1672107, 2019. 

Makhatov, B. M., Meldebekov, A. M., and Kulmanova, G. A.: Status and development prospects of sheep breeding in Kazakhstan [Sostoyaniye i perspektivy razvitiya ovcevodstva v Kazakhstane], Bulletin of Osh State University, 1–2, 382–389, https://doi.org/10.52754/16947452_2021_1_2_382, 2021. 

McNatty, K. P., Juengel, J. L., Wilson, T., Galloway, S. M., Davis, G. H., Hudson, N. L., Moeller, C. L., Cranfield, M., Reader, K. L., Laitinen, M. P., Groome, N. P., Sawyer, H. R., and Ritvos, O.: Oocyte-derived growth factors and ovulation rate in sheep, Reprod. Suppl., 61, 339–351, 2003. 

Melo, E., Silva, B., Castro, E., Silva, T., Paiva, S., Sartori, R., Franco, M., Souza, C. J., and Neves, J.: A novel mutation in the growth and differentiation factor 9 (GDF9) gene is associated, in homozygosis, with increased ovulation rate in Santa Ines sheep, Biol. Reprod., 78, 141, https://doi.org/10.1093/biolreprod/78.s1.141b, 2008. 

Molecular Evolutionary Genetics Analisys: MEGA software, https://www.megasoftware.net/, last access: 25 July 2023. 

Monsivais, D., Matzuk, M. M., and Pangas, S. A.: The TGF-β Family in the Reproductive Tract, CSH Perspect. Biol., 9, a022251, https://doi.org/10.1101/cshperspect.a022251, 2017. 

Muhaghegh Dolatabady, M. and Habibizad, J.: Single nucleotide polymorphisms (SNPs) of GDF9 gene in Bahmaei and Lak Ghashghaei sheep breeds and its association with litter size, Iranian Journal of Applied Animal Science, 9, 427–432, 2019. 

Mullen, M. P., Hanrahan, J. P., Howard, D. J., and Powell, R.: 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, PLoS One, 8, e53172, https://doi.org/10.1371/journal.pone.0053172, 2013.  

Muñoz-García, C., Vaquera-Huerta, H., Gallegos-Sánchez, J., Becerril-Pérez, C. M., Tarango-Arámbula, L. A., Bravo-Vinaja, Á., and Cortez-Romero, C.: Influence of FecGE mutation on the reproductive variables of Pelibuey ewes in the anestrus period, Trop. Anim. Health Pro., 53, 328, https://doi.org/10.1007/s11250-021-02755-7, 2021. 

Nilsson, E. E. and Skinner, M. K.: Growth and differentiation factor 9 stimulates progression of early primary but not primordial rat ovarian follicle development, Biol. Reprod., 67, 1018–1024, https://doi.org/10.1095/biolreprod.101.002527, 2002. 

Otsuka, F., McTavish, K. J., and Shimasaki, S.: Integral role of GDF-9 and BMP-15 in ovarian function, Mol. Reprod. Dev., 78, 9–21, https://doi.org/10.1002/mrd.21265, 2011. 

Paz, E., Quiñones, J., Bravo, S., Montaldo, H., and Sepúlveda, N.: Genotyping of BMPR1B, BMP15 and GDF9 genes in Chilean sheep breeds and association with prolificacy, Anim. Genet., 46, 98–99, https://doi.org/10.1111/age.12254, 2015. 

Polley, S., De, S., Brahma, B., Mukherjee, A., Vinesh, P. V., Batabyal, S., Arora, J. S., Pan, S., Samanta, A. K., Datta, T. K., and Goswami, S. L.: Polymorphism of BMPR1B, BMP15 and GDF9 fecundity genes in prolific Garole, Trop. Anim. Health Pro., 42, 985–993, https://doi.org/10.1007/s11250-009-9518-1, 2010. 

Rezaei, V., Esmailizadeh, A., Mehrgardi, A. A., and Qanatqestani, M. D.: Allelic polymorphism in exon 1 of GDF9 and exon 2 of BMP15 genes and its impact on litter size at lambing in Iran-Black sheep, Journal of Livestock Science and Technologies, 8, 57–65, https://doi.org/10.22103/jlst.2020.15007.1299, 2020. 

Saleh, A.: Utilization of molecular markers to detect some genes and mutations affecting economic traits in Egiptian sheep breeds, Figshare [thesis], https://doi.org/10.6084/m9.figshare.13026575.v2, 2020. 

Silva, B. D., Castro, E. A., Souza, C. J., Paiva, S. R., Sartori, R., Franco, M. M., Azevedo, H. C., Silva, T. A., Vieira, A. M., Neves, J. P., and Melo, E. O.: 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., 42, 89–92, https://doi.org/10.1111/j.1365-2052.2010.02078.x, 2011. 

Sun, H., Tian, X., and Wang, Y.: BMPR-IB, BMP15 and GDF9 as candidate genes for prolific trait in Tan sheep, Acta. Agricul. Boreali-Occident Sin., 18, 17–21, 2009. 

Våge, D. I., Husdal, M., Kent, M. P., Klemetsdal, G., and Boman, I. A.: A missense mutation in growth differentiation factor 9 (GDF9) is strongly associated with litter size in sheep, BMC Genet., 14, 1, https://doi.org/10.1186/1471-2156-14-1, 2013. 

Download
Short summary
Sheep breeding is one of the leading livestock sectors in agricultural countries, including Kazakhstan. A comprehensive study of the genetic aspects of these animals will solve many issues of the effectiveness of sheep breeding. This is important for such little-studied breeds as the Kazakh meat–wool sheep breed, which are multifetal. In this study, for the first time, we characterized this breed by looking at fecundity genes and found mutations that positively affect sheep fertility.