Rapid genotyping of MSTN gene polymorphism using high-resolution melting for association study in rabbits.
The myostatin (MSTN) gene is known as [GDF.sub.8] (growth differentiation factor 8) belonging to TGF-[beta] (transforming growth factor-beta) family, which is also considered as a negative regulator of skeletal muscle growth. It has been proved that the MSTN knockout mice have significantly larger body weight and more widespread increase in skeletal muscle mass than that of the wild-type (McPherron et al., 1997). The polymorphisms of this gene have been reported to be associated with production traits which make it a candidate gene to enhance the productivity of livestock and fowl (Lu et al., 2011; Han et al., 2012; Zhang et al., 2012). Rabbits have traits of economic importance like high rate of reproduction, early maturity, high growth rate and efficient feed utilization (Bindu et al., 2012). In previous studies, the rabbit MSTN was considered as one of the candidate genes of meat production traits and its single nucleotide polymorphisms (SNP) were only focused on the exon, intron and 3'-untranslated region (Fontanesi et al., 2008; Fontanesi et al., 2011; Kurkute et al., 2011). The mutation in 5'regulatory region of MSTN could lead to changes of the gene expression and thereby influence growth and carcass traits. However, the research of 5'regulatory region of rabbit MSTN is limited.
Direct sequencing technique is considered as the gold-standard approach for genotyping analysis and expected to have almost 100% sensitivity (Laurie et al., 2009). Many techniques have been developed to reduce the use of direct sequencing and control costs while maintaining sensitivity and specificity. The sensitivity of early methods, such as single-stranded conformation polymorphism analysis (SSCP) (Kunhareang et al., 2009) and denaturing gradient gel electrophoresis (DGGE) (Kirk et al., 2002), are low (70% to 80%) and require labor intensive and considerable technical skills. Denaturing high performance liquid chromatography (DHPLC) has high sensitivity (>96%) if amplicons are well-designed and analysis temperatures are optimized (Bunn et al., 2002).
High-resolution melting (HRM) analysis is a well established, closed tube, rapid, and cost-effective SNP genotyping technique based on the analysis of the melting profile of PCR products, using intercalating fluorescent dyes to monitor the transition during heating (Martino et al., 2010). As the temperature increases, the fluorescence decreases, producing a characteristic melting profile (Montgomery et al., 2007). The presence of a single nucleotide mutation in the amplicon will change the melting temperature (Tm) and the shape of the melting curve.
MSTN gene and its genetic polymorphism have been widely proposed to significantly associate with economic traits in farm animals. Therefore, developing a simple and economical genotyping method for MSTN gene is necessary for marker-assisted selection and breeding. HRM, characterized by high genotyping sensitivity and specificity, is being employed for MSTN genotyping in rabbit. In this study, direct sequencing was used to identify single nucleotide polymorphisms in the 5'regulatory region of rabbit MSTN. Then, genotyping of 222 randomly sampled individuals from the Ira rabbits, Champagne rabbits and Tianfu black rabbits was performed by HRM. We first compared the genotyping results of 47 samples with direct sequencing. Meanwhile, the association between MSTN polymorphism and rabbit production traits was subsequently analyzed by least square method.
MATERIAL AND METHODS
Animals, growth traits and carcass traits collection
This study was conducted in a total of 222 commercial meat rabbits including 88 Ira rabbits, 86 Champagne rabbits and 48 Tianfu black rabbits. The nutritional levels and feeding management were described previously (Zhang et al., 2011), with CP = 16%, DE = 10.8 MJ/kg. The carcasses were dissected according to the norms of the World Rabbit Science Association, and the data of growth traits and carcass traits were collected. The growth traits included 35-d-weight, 70-d-weight, and 84-d-weight; the carcass traits included eviscerated weight, semi-eviscerated weight, eviscerated slaughter rate, and semi-eviscerated slaughter rate (Blasco et al., 1993).
Genomic DNA isolation and mutation screening
Ear tissues were collected for genomic DNA extraction using AxyPrep Genomic DNA Miniprep Kit (Axygen, USA) and the DNA was stored at -20[degrees]C. A polymerase chain reaction (PCR) primer pair (Forward 5'-GAC TTGTGACAGACAGGGTTT-3' and Reverse 5'-AATCT TTATGGCTTCTATTCTTG-3') was newly designed in this study to amplify the 5'regulatory region according to the reference sequence of rabbit MSTN (GenBank accession number GU938461.1) by Primer 5.0, the PCR product size was 437 bp. The PCR reaction was performed with the following condition: one denaturation cycle at 94[degrees]C for 5 min; 35 cycles of 94[degrees]C for 30 s, 56[degrees]C for 30 s and 72[degrees]C for 40 s; and ended with an extension cycle at 72[degrees]C for 10 min. The 30 [micro]l reaction volume included 15 [micro]l 2xTaq PCR MasterMix (TIANGEN, Beijing, China), 3 [micro]l DNA template (20 ng/[micro]l), 9.6 [micro]l dd[H.sub.2]O, 1.2 [micro]l of each primer (10 pmol/[micro]l). PCR products from 16 rabbits were randomly selected to scan the genetic polymorphism in the 5'regulatory region using direct sequencing method by Shanghai Invitrogen Biotechnology Co., Ltd (Shanghai, P. R. China).
Genotyping using HRM
A primer pair (Forward 5'-TCATAGACCCTGACGAC ACT-3' and Reverse 5'-CTTTTCAGTAATGGCAAGC-3') was newly designed to amplify the small fragment containing the SNP revealed by direct sequencing, and the product size was 141 bp. The 10 [micro]l reaction volume included 5 [micro]l 1 x SsoFast EvaGreen supermix (Bio-Rad, USA), 1 [micro]l DNA template (20 ng/[micro]l), 3.2 [micro]l dd[H.sub.2]O, 0.4 [micro]l of each primer (10 pmol/[micro]l), and was pipetted into each well on the 96-well PCR plate. The cycling condition was: one denaturation cycle at 98[degrees]C for 3 min; 40 cycles of 5 s at 98[degrees]C and 25 s at 57[degrees]C; following the step of 95[degrees]C for 10 s; a melting step ramped from 65[degrees]C to 95[degrees]C rising by 0.2[degrees]C each step with date acquisition for 10 s. All samples were amplified in duplicate, each plate included DNA from sequencing-verified wild-type, heterozygote, homozygote and no template control (NTC). The PCR reaction was performed on Bio-Rad CFX96 real-time PCR detection system (Bio-Rad, Inc., Hercules, CA, USA). The amplicon was subsequently subjected to HRM analysis (Martino et al., 2010).
The DNA sequences were assembled and analyzed using the DNAstar program (DNAS Inc, Madison, WI, USA). HRM curve data were analyzed using the Precision Melt Analysis[TM] software (Martino et al., 2010). Allele frequencies, genotype frequencies, heterozygosity (He), effective number of alleles (Ne) and Polymorphic Information Content (PIC) were estimated based on the study of previously described (Bindu et al., 2012). The effects of genotype on the traits were analyzed by the general linear model (GLM) procedure of SAS 9.2 program according to the following statistical model:
[Y.sub.ijk] = [mu] + [G.sub.i] + [S.sub.j] + [B.sub.k] + [e.sub.ijk]
Where [Y.sub.ijk] is studied traits, | is the overall mean of observations, [G.sub.i] is the fixed genotype effect of the MSTN mutation site, [S.sub.j] is gender effect, [B.sub.k] is breed effect, and [e.sub.ijk] is the residual error.
RESULTS AND DISCUSSION
Identification of SNP in rabbit MSTN
In previous studies, SNPs in rabbit MSTN was focused on the exon, intron and 3'-untranslated region. Fontanesi et al. (2008) sequenced the entire MSTN in four different breeds of rabbit and only detected one SNP in intron 2 (c.747+34C>T). Then Fontanesi et al. (2011) re-sequenced the MSTN including the three coding exons in 14 rabbits with different conformations and muscle mass, and identified another three SNPs of MSTN in exon 1 (c.108C>T), exon 2 (c.713T>A) and 3'-untranslated region (c.*194A>G). Kurkute et al. (2011) detected all three exons of MSTN in three rabbit breeds and found two mutations in exon 1 (c.1A>G) and exon 3 (c.799G>A).
In the present study, a T/C variant at -125 bp (relative to ATG start codon) of rabbit MSTN was identified by direct sequencing method (Figure 1). The T/C variant is a newly discovered variant in 5'regulatory region of MSTN. In silico analysis revealed that the allele C increased a transcription factor binding site of AP-2 and decreased a transcription factor binding site of CAAT (http://tfbind.hgc.jp/). The CAAT box is a motif present in the conserved sequence upstream of the start points in eukaryotic transcription units, which is recognized by a large group of transcription factors. AP-2 transcription factors are localized predominantly in the nucleus where they bind to target sequences and regulate transcription of target genes (Dawid et al., 2005). The change of these transcription factor binding sites might affect RNA transcription, which should be further investigated.
Genotyping of SNP in rabbit MSTN by HRM
HRM was performed on sequencing-verified samples in order to determine the Tm and characterize of the melting curve profiles of different genotypes. In this study, each genotype produced unique melt curves (amplitude and/or shape of the curves) and could be discriminated directly (Figure 2). The genotyping of SNP in MSTN by HRM was done as following: first, the curves were auto grouped by the software based on the melt transition, and then the observer confirmed the difference between Tm and melt curves. There was no difference in the classification by the observer compared to that grouped by the software.
The homozygous CC (79.4[degrees]C) melted at a higher temperature than the homozygous TT (79.0[degrees]C) due to the higher stability of the G [equivalent to] C pairing compared to A = T. The difference in temperature could unambiguously differentiate the two homozygous sequences. The heterozygotes were identified based on the changes in the shapes of the melt curves as compared to homozygotes. The presence of both forms of homozygotes together with heteroduplexes results in a change in the melting pattern characterized by a flatter curve shape/melting transition (Raghavendra et al., 2011), as observed in the normalized melting curves (Figure 2).
Validation of sensitivity and specificity of genotyping by HRM
We selected 47 samples which were on the same plate to be sequenced for verifying the sensitivity and specificity of genotyping by HRM. Sequencing was performed by Shanghai Invitrogen Biotechnology Co., Ltd (Shanghai, P. R. China). The sensitivity and specificity of HRM were 0. 96 (45 true positive/(45 true positive+2 false negative)) and 0.98 (92 true negative/(92 true negative+2 false positive)), respectively. Recent studies have also reported the sensitivity of HRM close to 100% except in the situation where amplicons have a high GC content (Garritano et al., 2009; Laurie et al., 2009).
Direct sequencing is the gold-standard approach for detecting SNPs, but involves additional steps leading to a higher turnaround time. HRM is a closed tube, fast method following PCR amplification and eliminates the possibility of template contamination, where samples can be analyzed in a single run (Raghavendra et al., 2011). HRM does not have 100% sensitivity; however, mutation scanning is a trade-off between the speed, cost, and sensitivity. It is worth to emphasizing how rapid this technique is: a plate of 96 samples can be melted in 10 min (Laurie et al., 2009). SsoFast EvaGreen supermix was used in this study. EvaGreen dye is a fluorescent nucleic acid dye with spectral properties similar to SYBR Green I and fluorescein. Unlike SYBR Green I, EvaGreen dye exhibits very low PCR inhibition. Thus, it can be used at high concentrations to generate greater fluorescent signals and provide increased sensitivity. Consequently, HRM is a better candidate to detect sequence mutations on a large scale, as it is much quicker than direct sequencing.
Genotype and allele frequencies in 5'regulatory region of rabbit MSTN
The genotype and allele frequencies in the 5'regulatory region of MSTN in different rabbit breeds is shown in Table 1. Allele C was the predominant allele with frequencies of 0.56 in Iraq, 0.58 in Champagne and 0.76 in Tianfu rabbits, respectively. The frequencies of genotype TC for Iraq and Champagne rabbits were the highest (0.64 and 0.59, respectively), but the frequency of genotype CC was the highest in Tianfu rabbits (0.56). The polymorphic site showed high heterozygosity (He = 0.48) and a high effective number of alleles (Ne = 1.91), and the genetic diversity was reasonably informative (0.25 < PIC < 0.50) in the three rabbits breeds. The genetic diversity revealed higher genetic variation and a selective potential that could be expected to get more genetic progress. In previous studies, a C-T transition at position 34 of MSTN intron 2 has been reported in New Zealand White, Soviet Chinchilla, Belgian Hare and several other rabbit breeds (Fontanesi et al., 2008; Bindu et al., 2012). In the present study, we randomly selected three rabbit breeds that were different from previous studies. The results showed similar frequencies and genetic characteristics.
Association between MSTN genotypes and production traits in rabbits
Some studies have showed that mutations in MSTN can inactivate its expression and further affect muscle fiber development and cause dramatic muscularity and a "double muscling" phenomenon in cattle (Wiener et al., 2002; Esmailizadeh et al., 2008; Grisolia et al., 2009). In the present study, the association analysis (Table 2) showed that different genotypes had a significant effect on the 84-d-weight of rabbits, individuals with the genotype TC indicated a significant association with higher 84-d-weight than genotype CC (p = 0.047). But MSTN genotypes had no significant effect on other production traits. Similar findings have been reported previously. Fontanesi et al. (2011) indicated that MSTN is not an important source of variability for performance traits, at least in the rabbit population analyzed. Bindu et al. (2012) also showed that a single nucleotide polymorphism (C-T transition) at position 34 of intron 2 of rabbit MSTN was associated with higher body weight but the difference in body weight among different genetic groups was not statistically significant.
As previously mentioned, the T/C variant at -125 bp caused the change of transcription factor binding sites which might affect RNA transcription. Sadkowski et al. (2008) also showed that polymorphism in the 5'-flanking region of MSTN was associated with the level of mRNA in skeletal muscle of bulls. Lin et al. (2002) showed that MSTN expression levels of Asian pigs increased markedly at 3 moths and maintained a constant level thereafter. Sun et al. (2012) also demonstrated that MSTN expression levels in muscle of Hu sheep were lowest at two days old, and tended to increase with the increasing age before 60 days. For these reasons, we could assume that MSTN expression levels might be relatively low in the initial growth process. Therefore, the T/C variant at -125 bp of rabbit MSTN had no significant effect on 35-d-weight and 70-d-weight. As the rabbit grows, the expression levels should gradually increase and the inhibitory effect on muscle growth should also gradually become obvious. At 84-d, the rabbit growth reached a relatively stable period and the inhibitory effect on muscle growth of MSTN might be most obvious. That caused the individuals with the genotype TC to have a significantly higher 84-d-weight than genotype CC (p = 0.047).
These results showed that HRM could be effectively used for genotyping analysis of the MSTN gene. The T/C variant in 5'regulatory region of MSTN might be one of the candidate SNP locus affecting the trait of 84-d-weight.
Bindu, K. A., A. Raveendran, S. Antony and K. V. Raghunandanan. 2012. Association of myostatin gene (MSTN) polymorphism with economic traits in rabbits. In: Fibre Production in South American Camelids and Other Fibre Animals (Ed. M. A. Perez-Cabal, J. P. Gutierrez, I. Cervantes and M. J. Alcalde). Wageningen Academic Publishers, Wageningen, Netherlands. pp. 131-133.
Blasco, A., J. Ouhayoun and G. Masoero. 1993. Harmonization of criteria and terminology in rabbit meat research. World Rabbit Sci. 1:03-10.
Bunn, C. F., C. J. Lintott, R. S. Scott and P. M. George. 2002. Comparison of SSCP and DHPLC for the detection of LDLR mutations in a New Zealand cohort. Hum. Mutat. 19:311.
Dawid, E., S. Buhl, S. Weber, R. Jager and H. Schorle. 2005. The AP-2 family of transcription factors. Genome Biol. 6:246.
Esmailizadeh, A. K., C. D. K. Bottema, G. S. Sellick, A. P. Verbyla, C. A. Morris, N. G. Cullen and W. S. Pitchford. 2008. Effects of the myostatin F94L substitution on beef traits. J. Anim. Sci. 86:1038-1046.
Fontanesi, L., M. Tazzoli, E. Scotti and V. Russo. 2008. Analysis of candidate genes for meat production traits in domestic rabbit breeds. In: 9th World Rabbit Congress, Verona, Italy. pp. 79-84.
Fontanesi, L., E. Scotti, A. Frabetti, D. Fornasini, A. Picconi and V. Russo. 2011. Identification of polymorphisms in the rabbit (Oryctolagus cuniculus) myostatin (MSTN) gene and association analysis with finishing weight in a commercial rabbit population. Anim. Genet. 42:339.
Garritano, S., F. Gemignani, C. Voegele, T. Nguyen-Dumont, F. L. Calvez-Kelm, D. D. Silva, F. Lesueur, S. Landi and S. V. Tavtigian. 2009. Determining the effectiveness of High Resolution Melting analysis for SNP genotyping and mutation scanning at the TP53 locus. BMC. Genet. 10:5.
Grisolia, A. B., G. T. D'Angelo, L. R. P. Neto, F. Siqueira and J. F. Garcia. 2009. Myostatin (GDF8) single nucleotide polymorphisms in Nellore cattle. Genet. Mol. Res. 8:822-830.
Han, S. H., I. C. Cho, M. S. Ko, E. Y. Kim, S. P. Park, S. S. Lee and H. S. Oh. 2012. A promoter polymorphism of MSTN g.-371T>A and its associations with carcass traits in Korean cattle. Mol. Biol. Rep. 39:3767-3772.
Kirk, B. W., M. Feinsod, R. Favis, R. M. Kliman and F. Barany. 2002. Single nucleotide polymorphism seeking long term association with complex disease. Nucleic Acids Res. 30: 3295-3311.
Kunhareang, S., H. Zhou and J. G. H. Hickford. 2009. Allelic variation in the porcine MYF5 gene detected by PCR-SSCP. Mol. Biotechnol. 41:208-212.
Kurkute, A. S., A. K. Tripathi, N. Shabir, C. V. Jawale, U. V. Ramani, A. M. Pande, D. N. Rank and C. G. Joshi. 2011. Molecular cloning and characterization of rabbit myostatin gene. IIOAB J. 2:1-7.
Laurie, A. D. and P. M. George. 2009. Evaluation of high-resolution melting analysis for screening the LDL receptor gene. Clin. Biochem. 42:528-535.
Lin, C. S., Y. C. Wu, Y. L. Sun and M. C. Huang. 2002. Postnatal expression of growth/differentiation factor-8 (GDF-8) gene in European and Asian pigs. Asian-Aust. J. Anim. 15:1244-1249.
Lu, J., S. Hou, W. Huang, J. Yu and W. Wang. 2011. Polymorphisms in the myostatin gene and their association with growth and carcass traits in duck. Afr. J. Biotechnol. 10: 11309-11312.
Martino, A., T. Mancuso and A. M. Rossi. 2010. Application of high-resolution melting to large-scale, high-throughput SNP genotyping: a comparison with the TaqMan[R] method. J. Biomol. Screen 15:623-629.
McPherron, A. C., A. M. Lawler and S. J. Lee. 1997. Regulation of skeletal muscle mass in mice by a new TGF-[beta] superfamily member. Nature 387:83-90.
Montgomery, J., C. T. Wittwer, R. Palais and L. Zhou. 2007. Simultaneous mutation scanning and genotyping by high-resolution DNA melting analysis. Nat. Protoc. 2:59-66.
Raghavendra, A., A. Siji, T. S. Sridhar, K. Phadke and A. Vasudevan. 2011. Evaluation of High Resolution Melting analysis as an alternate tool to screen for risk alleles associated with small kidneys in Indian newborns. BMC. Nephrol. 12:60.
Sadkowski, T., M. Jank, L. Zwierzchowski, E. Siadkowska, J. Oprzadek and T. Motyl. 2008. Gene expression profiling in skeletal muscle of Holstein-Friesian bulls with single-nucleotide polymorphism in the myostatin gene 5'-flanking region. J. Appl. Genet. 49:237-250.
Sun, W., D. Li, P. Wang, H. H. Musa, J. T. Ding, B. C. Li, Y. H. Ma, W. J. Guan, M. X. Chu, L. Chen, Y. F. Zhang, W. Z. Wu and H. Zhou. 2012. Postnatal expression of myostain (MSTN) and myogenin (MYoG) genes in Hu sheep of China. Afr. J. Biotechnol. 11:12246-12251.
Wiener, P., J. A. Smith, A. M. Lewis, J. A. Woolliams and J. L. Williams. 2002. Muscle-related traits in cattle: the role of the myostatin gene in the South Devon breed. Genet. Sel. Evol. 34: 221-232.
Zhang, C., Y. Liu, D. Xu, Q. Wen, X. Li, W. Zhang and L. Yang. 2012. Polymorphisms of myostatin gene (MSTN) in four goat breeds and their effects on Boer goat growth performance. Mol. Biol. Rep. 39:3081-3087.
Zhang, G. W., H. Z. Wang, S. Y. Chen, Z. C. Li, W. X. Zhang and S. J. Lai. 2011. A reduced incidence of digestive disorders in rabbits is associated with allelic diversity at the TLR4 locus. Vet. Immunol. Immunopathol. 144:482-486.
Jin Peng, Gong-Wei Zhang, Wen-Xiu Zhang, Yun-Fu Liu, Yu Yang and Song-Jia Lai *
Institute of Animal Genetics and Breeding, Sichuan Agricultural University, Chengdu Campus, Chengdu 611130, China
* Corresponding Author: Song-Jia Lai. Tel: +86-28-86290987, Fax: +86-28-86290987, E-mail: email@example.com
Submitted Jul. 9, 2012; Accepted Sept. 10, 2012; Revised Sept. 24, 2012
Table 1. Genotype and allele frequencies in 5'regulatory region of rabbit MSTN Breeds (n) Genotype frequency (n) TT TC CC Iraq rabbit (88) 0.12 (11) 0.64 (56) 0.24 (21) Champagne rabbit (86) 0.13 (11) 0.59 (51) 0.28 (24) Tianfu rabbit (48) 0.04 (2) 0.40 (19) 0.56 (27) Total (222) 0.11 (24) 0.57 (126) 0.32 (72) Breeds (n) Allele Genetic characteristic frequency T C He Ne PIC Iraq rabbit (88) 0.44 0.56 0.49 1.97 0.37 Champagne rabbit (86) 0.42 0.58 0.49 1.96 0.37 Tianfu rabbit (48) 0.24 0.76 0.36 1.57 0.30 Total (222) 0.39 0.61 0.48 1.91 0.36 He = Heterozygosity, Ne = Effective number of alleles, PIC = Polymorphism information content, n = Number. Table 2. The association analysis between MSTN genotypes and production traits in rabbits Traits Genotypes TT 35-d-weight (g) 789.62[+ or -]23.67 (a) 70-d-weight (g) 2,166.59[+ or -]35.15 (a) 84-d-weight (g) 2,615.55[+ or -]48.27 (ab) Eviscerated weight (g) 1,314.82[+ or -]28.47 (a) Eviscerated slaughter 0.53[+ or -]0.0047 (a) percentage (%) Semi-eviscerated 1,420.15[+ or -]30.64 (a) weight (g) Semi-eviscerated 0.57[+ or -]0.0048 (a) slaughter percentage (%) Traits Genotypes TC 35-d-weight (g) 796.57[+ or -]10.73 (a) 70-d-weight (g) 2,171.64[+ or -]15.92 (a) 84-d-weight (g) 2,630.34[+ or -]21.87 (a) Eviscerated weight (g) 1,316.75[+ or -]12.90 (a) Eviscerated slaughter 0.52[+ or -]0.0021 (a) percentage (%) Semi-eviscerated 1,426.00[+ or -]13.88 (a) weight (g) Semi-eviscerated 0.57[+ or -]0.0022 (a) slaughter percentage (%) Traits Genotypes CC 35-d-weight (g) 801.64[+ or -]13.35 (a) 70-d-weight (g) 2,144.00[+ or -]19.82 (a) 84-d-weight (g) 2,560.22[+ or -]27.21 (b) Eviscerated weight (g) 1,321.60[+ or -]16.05 (a) Eviscerated slaughter 0.53[+ or -]0.0026 (a) percentage (%) Semi-eviscerated 1,428.92[+ or -]17.27 (a) weight (g) Semi-eviscerated 0.57[+ or -]0.0027 (a) slaughter percentage (%) Values are presented by the least squares means [+ or -]standard error. In the same row, the same letter means not significantly, different letters mean significant difference at 0.05 levels.
|Printer friendly Cite/link Email Feedback|
|Author:||Peng, Jin; Zhang, Gong-Wei; Zhang, Wen-Xiu; Liu, Yun-Fu; Yang, Yu; Lai, Song-Jia|
|Publication:||Asian - Australasian Journal of Animal Sciences|
|Date:||Jan 1, 2013|
|Previous Article:||Gene co-expression analysis to characterize genes related to marbling trait in Hanwoo (Korean) cattle.|
|Next Article:||A study on BMPR-IB genes of Bayanbulak sheep.|