What Type of Mutation Occurred in the Belgian Blue Myostatin
J Biomed Biotechnol. 2010; 2010: 542945.
Analysis of Horse Myostatin Gene and Identification of Single Nucleotide Polymorphisms in Breeds of Different Morphological Types
Stefania Dall'Olio
1DIPROVAL, Sezione di Allevamenti Zootecnici, Faculty of Agriculture, University of Bologna, 42123 Reggio Emilia, Italy
Luca Fontanesi
1DIPROVAL, Sezione di Allevamenti Zootecnici, Faculty of Agriculture, University of Bologna, 42123 Reggio Emilia, Italy
Leonardo Nanni Costa
1DIPROVAL, Sezione di Allevamenti Zootecnici, Faculty of Agriculture, University of Bologna, 42123 Reggio Emilia, Italy
Marco Tassinari
2DIMORFIPA, Sezione di Zootecnia, Faculty of Veterinary Medicine, University of Bologna, 40064 Ozzano Emilia, Italy
Laura Minieri
3Department of Human and General Physiology, Faculty of Medicine, University of Bologna, 40127 Bologna, Italy
Adalberto Falaschini
2DIMORFIPA, Sezione di Zootecnia, Faculty of Veterinary Medicine, University of Bologna, 40064 Ozzano Emilia, Italy
Received 2009 Aug 3; Revised 2010 Feb 20; Accepted 2010 Apr 27.
Abstract
Myostatin (MSTN) is a negative modulator of muscle mass. We characterized the horse (Equus caballus) MSTN gene and identified and analysed single nucleotide polymorphisms (SNPs) in breeds of different morphological types. Sequencing of coding, untranslated, intronic, and regulatory regions of MSTN gene in 12 horses from 10 breeds revealed seven SNPs: two in the promoter, four in intron 1, and one in intron 2. The SNPs of the promoter ({"type":"entrez-nucleotide","attrs":{"text":"GQ183900","term_id":"260408325","term_text":"GQ183900"}}GQ183900:g.26T>C and {"type":"entrez-nucleotide","attrs":{"text":"GQ183900","term_id":"260408325","term_text":"GQ183900"}}GQ183900:g.156T>C, the latter located within a conserved TATA-box like motif) were screened in 396 horses from 16 breeds. The g.26C and the g.156C alleles presented higher frequency in heavy (brachymorphic type) than in light breeds (dolichomorphic type such as Italian Trotter breed). The significant difference of allele frequencies for the SNPs at the promoter and analysis of molecular variance (AMOVA) on haplotypes indicates that these polymorphisms could be associated with variability of morphology traits in horse breeds.
1. Introduction
Myostatin, encoded by the MSTN gene (previously referred to as GDF8), is a member of the transforming growth factor β superfamily that normally acts to limit skeletal muscle mass by regulating both the number and growth of muscle fibres [1, 2]. MSTN is synthesized as precursor and upon proteolytic processing gives an N-terminal latency-associated peptide, termed myostatin propeptide or LAP-fragment, and a smaller mature peptide at the C-terminus [3]. The MSTN gene, composed of three exons and two introns, has been characterized in rodents [1], humans [4], and several livestock species [3, 5–8]. Natural mutations that decrease the amounts of myostatin and/or inhibit its function have been identified in a human subject [9] and in several cattle [2, 3, 10–13], sheep [14–16], and dog [17] breeds. In Belgian Blue, Piedmontese, Marchigiana, and other cattle breeds, loss-of-function mutations within the coding sequence of the MSTN gene determine increased skeletal muscle mass, relevant in shoulders and thighs, and the produced phenotype is known as "double-muscling" [2, 3, 10–13, 18]. These polymorphisms have, in several cases, effects on growth, reproductive, performances, and carcass quality traits [3, 18, 19]. In the Whippet dog breed a mutation in the third exon determining a premature stop codon causes an increased muscle mass phenotype in homozygous state and enhanced racing performance in heterozygous dogs [17]. In two Norwegian sheep breeds, two different mutations in the MSTN coding region are associated with carcass conformation and fatness [15, 16]. In addition, in other sheep and in pigs, mutations identified in non coding regulatory regions affect the level of MSTN gene expression and/or are associated with growth, muscle mass, and other carcass traits [8, 14, 20, 21].
In horse (Equus caballus), only few studies examined the MSTN gene so far. Hosoyama et al. [22] isolated and sequenced an MSTN cDNA from a Thoroughbred horse and Caetano et al. [23] mapped this gene to equine chromosome 18. Mutations in the equine MSTN gene have been identified only recently in Thoroughbred breed [24].
Different horse breeds present a variety of morphological phenotypes that have been used to group breeds into a few classes. However, no system provides a robust classification in which each breed could have an unequivocal assignment. Based on size traits and build, horse breeds are categorized in draught (or heavy), light, and pony (or animals that mature at less than 148 cm high, usually used as riding school and children's mounts) [25]. Considering skeletal structure, proportions, zoometrical indices, length, and volume of muscling, that, in turn, reflect the selective goals and uses of the horse breeds, they are categorized in brachymorphic, mesomorphic, dolichomorphic, and intermediate types (such as meso-dolichomorphic) [26]. Brachymorphic horses (corresponding to draught horses), traditionally referred to as cold-blooded horses in relation to their quiet and calm temperament, are tall in stature, heavy boned, and extremely muscular with short and thick muscles and slow twitch oxidative fibers for slow contraction. They most likely develop strength and power, and their conformation is well suited for pulling carriage, draught power and meat production. Dolichomorphic horses (corresponding to light horses) are characterized by longer bodies and long and thin muscles mainly constituted by fast twitch glycolytic fibers. They are selected for sport purposes, fast running, and high speeds. Examples of Italian breeds representative of these two extreme phenotypes are reported in Figure 1. Mesomorphic type is characterized by a lighter physical structure than brachymorphic but still powerful and compact with massive muscling. This group also includes some breeds with draft-type qualities and classified as ponies based on their withers height (such as Bardigiano and Haflinger breeds). The mesomorphic horses are usually used for pleasure and riding. In addition, several breeds (like local breeds with influence of Oriental, Thoroughbred, and Iberian halfbreed and descendents) have characteristics of both mesomorphic and dolichomorphic types (referred to as meso-dolichomorphic) [26].
Horses of extreme and opposite morphological types:(a) Rapid Heavy Draft (brachymorphic type or heavy) and (b) Italian Trotter (dolichomorphic type or light).
For the important pleiotropic effects of the MSNT gene, including its role on muscle mass development, polymorphisms in this gene could contribute to explain the morphological variability among horse breeds. Here we sequenced the MSTN gene, including regulatory regions, in several horse breeds and identified a few polymorphisms that were used to evaluate their potential association with different morphological types.
2. Materials and Methods
2.1. Animals and Horse Breeds Classification Based on Different Morphological Types
A total of 396 minimal related horses belonging to 16 breeds were sampled in different farms or stables. Details of the horse breeds involved in the analysis are given in Table 1.
Table 1
Numbers, types, origin, uses, and morphological information of the analysed horse breeds*.
| International breed name (Home page of Breed National Associations) | Types§ | Number of horses | Origin of the samples | Uses (adult animals, males/females: withers height, cm; live weight, kg) |
|---|---|---|---|---|
| Rapid Heavy Draft (http://www.anacaitpr.it/) | B | 81 | North East of Italy | Meat, draught power (160/155 cm; 700/570 kg) |
| Noric or Pinzgauer (http://www.aia.it/tecnico/equini/a_norico.htm) | B | 26 | Alps, Northern Italy | Draught power, driving, carriage, dressage, meat (155/153 cm; 700/650 kg) |
| Bardigiano (http://www.bardigiano.it/) | M | 34 | Emilia Romagna region, Northern Italy | Meat, agricultural work, hobby (143/142 cm; 530/480 kg) |
| Haflinger or Avelignese (http://www.haflinger.it/) | M | 31 | Northern Italy | Hobby, meat, draught power (142/140 cm; 450/450 kg) |
| Lipizzan or Lipizzaner (http://www.aia.it/tecnico/equini/a_lipizzano.htm) | M | 12 | Rome, Central Italy | Dressage, riding, carriage (164/150 cm; 550/480 kg) |
| Murgese (http://www.aia.it/tecnico/equini/download/DM%20MURGESE.pdf) | M | 12 | Apulia region, Southern Italy | Riding, equestrian tourism, draught power (164/162 cm; 550/480 kg) |
| Tolfetano (http://www.aia.it/tecnico/equini/a_tolfetan.htm) | M | 7 | Latium region, Central Italy | Riding (160/150 cm) |
| Uruguayan Creole or Criollo (http://www.caballoscriollos.com.uy/homepage_en.php) | M | 35 | Uruguay | Draught power, government of livestock, endurance (144/144 cm; 420/420 kg) |
| Italian Saddle (http://www.aia.it/ANACSI/Frame.htm) | M-D | 30 | Country wide, Italia | Riding, equestrian tourism, drift competitions (jumping, dressage, complete, endurance) (min 156 cm). |
| Maremmano (http://www.anamcavallomaremmano.com) | M-D | 13 | Tuscany region, Central Italy | Riding, equestrian tourism, draught power, government of livestock (165/162 cm; 500/450 kg) |
| Quarter Horse (http://www.aiqh.it) | M-D | 15 | Emilia Romagna region, Italy | Riding, equestrian tourism (160 cm; 550/350 kg) |
| Salernitano or Salernitano-Persano (http://www.agraria.org/equini/salernitano.htm) | M-D | 4 | Campania region, Southern Italy | Riding, equestrian tourism (158/150 cm; 500/450 kg) |
| Spanish Purebred or Andalusian (http://www.aipre.com) | M-D | 10 | Emilia Romagna region, Italy | Sport, hobby (160/155 cm; 570 kg) |
| Ventasso (http://www.aia.it/tecnico/equini/a_ventass.htm) | M-D | 8 | Emilia-Romagna region, Northern Italy | Equestrian tourism (164/152 cm) |
| Italian Trotter (http://www.anact.it/) | D | 67 | Country wide, Italy | Selected for competitive merit as trotter at short and medium distances (160–145 cm) |
| Throroughbred (http://www.anacpurosangue.com/) | D | 11 | Country wide, Italy | Racehorses (gallop) and used for genetic improvement of Italian local breeds. |
These horse breeds were classified as brachymorphic (B), mesomorphic, (M), meso-dolichomorphic (M-D), and dolichomorphic (D) (Table 1) as indicated in the homepage of their own Breed National Associations based on linear measures (height at withers, chest girth, and cannon circumference), structure and anamorphosis index (AI= (chest girth)2∗100/height at wither), and based on bibliographic data [27–30].
All horses were registered in the Stud Books or in the Italian Anagraphic Register constituted for local ethnic groups (Noric, Salernitano, Tolfetano, and Ventasso). The Lipizzan samples (Lipizzan Italian Stud, Monterotondo, Italy) included all six classical stallion lines: Conversano, Favory, Maestoso, Neapolitano, Pluto, and Siglavy.
2.2. PCR and Sequencing
Genomic DNA was extracted from hair roots following standard procedures. Ten primer pairs (Table 2) that amplify different MSTN regions were designed using Primer 3 (http://frodo.wi.mit.edu/primer3/input.htm) software. PCR reactions were performed in a final volume of 20μL containing 10–80 ng of equine genomic DNA, 250 mM of each dNTP, 10 pmol of each primer, 1 U of EuroTaq DNA polymerase (EuroClone Ltd., Paington, Devon, UK), or 1 U of TaKaRa Ex Taq DNA Polymerase (TaKaRa Bio Inc., Shiga, Japan) and 1× PCR buffer with MgCl2 concentration specific for each primer pair (Table 2). PCR conditions were: an initial step at 95°C for 5 minutes, 35 cycles of 95°C for 30 s, specific annealing temperature for each primer pair for 30 s, 72°C for specific reaction times for different primer pair (Table 2), and a final step at 72°C for 9 minutes. Genomic DNA obtained from 12 horses of 10 breeds (1 Bardigiano, 1 Haflinger, 1 Italian Saddle, 2 Italian Trotter, 1 Noric, 2 Rapid Heavy Draft, 1 Salernitano, 2 Throroughbred, and 1 Ventasso) constituted the sequencing panel. PCR fragments obtained from the sequencing panel with primer pairs 1–10 were purified using the QIAquick PCR Purification Kit (Qiagen, Düsseldorf, Germany) and sequenced on both strands using the BigDye Cycle Sequencing kit v.3.1 (Applied Biosystems, Foster City, CA, USA). Sequencing reactions were electrophoresed in a capillary sequencer (Applied Biosystems).
Table 2
Primers, PCR conditions, amplified regions, and use of the PCR products.
| Primer pair | F: Forward sequence (5′-3′) R: Reverse sequence (5′-3′) | Primer annealing regions | Amplified region (bp)* | PCR conditions§ | Use of the PCR products |
|---|---|---|---|---|---|
| 1# | 1F: TCAGGGAAACAAGTTTCTCAAAT | Promoter | 1–484 | 58/1.5/45/ | Sequencing, |
| 1R: TGCTCCACAATGAATCTCG | Promoter | (484) | E | PCR-RFLP | |
| 2# | 2F: TGAATCAGCTCACCCTTGAC | Promoter | 369–727 | 62/0.8/45/ | Sequencing |
| 2R: CCAGCAACAATCAGCATAAA | Exon 1 | (360) | E | ||
| 3 ¶ | 3F TGTGCTGATTCTTGCTGGTC | Exon 1 | 710–947 | 58/1.5/45/ | Sequencing |
| 3R: ATCAATCAGTTCCCGGAGTG | Exon 1 | (238) | E | ||
| 4 ¶ | 4F: GACCCGTCAAGACTCCTACA | Exon 2 | 2982–3226 | 58/1.0/45/ | Sequencing |
| 4F: TGGGAAGGTTACAGCAAGA | Exon 2 | (245) | E | ||
| 5 ¶ | 5F: AGGCCAATTACTGCTCTGGA | Exon 3 | 5430–5744 | 58/1.5/45/ | Sequencing |
| 5R: ATACTCTAGGCTTATAGCCT | 3' UTR | (315) | E | ||
| 6 ¶ | 6F: CACTCCGGGAACTGATTGAT | Exon 1 | 928–3098 | 58/0.8/130 | Sequencing |
| 6R: CGCCTGGGTTCATGTCAAGT | Exon 2 | (2171) | /T | ||
| 7 ¶ | 7F: AGGCAGGCACATTGCTTAAT | Intron 1 | 1808–2287 | 58/0.8/45/ | Sequencing |
| 7R: GAATGTTATATTCAGGCTATCTCAA | Intron 1 | (480) | E | ||
| 8 ¶ | 8F: AAATGTGACATAAGCAAATGATTAG | Intron 2 | 5152–5534 | 62/0.8/45/ | Sequencing |
| 8R: AGCAGGGGCCTGCTGAACCTCTGGG | Exon 3 | (383) | E | ||
| 9 ¶ | 9F: TGCAAAATTGGCTCAAACAG | Exon 2 | 3135–5361 | 55/0.8/130 | Sequencing |
| 9R: CAGCATCGAGATTCTGTGGA | Exon 3 | (2227) | /T | ||
| 10 ¶ | 10F: CCCCCAGAAGAGTGTCAAAT | Intron 2 | 3627–4822 | 60/0.8/130 | Sequencing |
| 10R: TCTTTACTTGGGGAAACTTGGA | Intron 2 | (1196) | /T | ||
| 11 | 1F: TCAGGGAAACAAGTTTCTCAAAT | Promoter | 1–204 | 53/1.2/35/ | PCR-RFLP |
| 11R: ACTTCCTCAGAAATTAAGATTTAAT | Promoter | (204) | E |
2.3. Sequence Analysis, Polymorphism Identification, and Genotyping
Sequences were aligned and processed with the help of the BioEdit software v.7.0.5.2. Polymorphisms were identified by visual inspection of the electropherograms and sequences were aligned with ClustalW2 program (http://www.ebi.ac.uk/Tools/clustalw2/index.html) and using BLASTN (http://blast.ncbi.nlm.nih.gov/Blast.cgi). The horse MSTN promoter sequence was analysed in silico for the presence of putative transcription factor binding sites using MatInspector (http://www.genomatix.de/) bioinformatics tool. This region of horse MSTN gene was aligned with that of cattle ({"type":"entrez-nucleotide","attrs":{"text":"AJ310751","term_id":"15072712","term_text":"AJ310751"}}AJ310751), goat ({"type":"entrez-nucleotide","attrs":{"text":"AY827576","term_id":"56181514","term_text":"AY827576"}}AY827576), human ({"type":"entrez-nucleotide","attrs":{"text":"AX058992","term_id":"12311262","term_text":"AX058992"}}AX058992), mouse ({"type":"entrez-nucleotide","attrs":{"text":"AY204900","term_id":"28797966","term_text":"AY204900"}}AY204900), pig (AY8641281), and sheep ({"type":"entrez-nucleotide","attrs":{"text":"DQ530260","term_id":"105921189","term_text":"DQ530260"}}DQ530260) to identify evolutionary conserved motifs.
PCR-RFLP protocols were designed to genotype two identified SNPs (g.26T>C and g.156T>C) in the sampled horses. To genotype the g.26T>C SNP, the amplified products of 484 bp obtained with primer pair 1 (Table 2) were digested with RsaI (recognition sequence: GT↓AC). Briefly, 5–10μ L of PCR reaction was restricted with 2.5 U of RsaI (Fermentas, Vilnius, Lithuania) at 37°C overnight and the resulting fragments (g.26T allele=484 bp; g.26C allele=437 bp + 47 bp) were resolved on 2.0 % agarose gels stained with ethidium bromide. The g.156T>C SNP was genotyped amplifying a fragment of 204 bp with primer pair 11 (Table 2) that inserted an artificial restriction site (with a mismatched reverse primer) for SspI (recognition sequence: AAT↓ATT) when allele g.156T occurred. The obtained fragments (g.156T allele=179 bp + 25 bp; g.156C allele=204 bp) resulting from digestion of 5–10μL of PCR reaction with 2.5 U of SspI (Fermentas) at 37°C overnight were electrophoresed in 3.5% agarose gels and visualized with ethidium bromide.
2.4. Statistical Analysis
Allele and genotype frequencies, observed and expected heterozygosity, and F st were calculated using PopGene software v. 1.32 [31]. F st is a measure of population differentiation based on genotypic data. Allele frequencies among the four groups (B, M, M-D, and D) were compared using Fisher's exact test. The haplotypes of the two promoter SNPs were reconstructed using PHASE program v. 2.0 [32]. ARLEQUIN software v. 3.1 (http://cmpg.unibe.ch/software/arlequin3) was used for the analysis of molecular variance (AMOVA) testing the effect of the morphological types in population differentiation with a model including types (four levels: B, M, M-D and D; and two levels: B + M and M-D + D), types/breeds, individuals/breeds, and individuals.
3. Results and Discussion
3.1. Horse MSTN Genomic Structure and Sequence Analysis
Sequenced fragments of the horse MSTN gene were assembled into one sequence of 5724 bp (submitted to GenBank under accession number {"type":"entrez-nucleotide","attrs":{"text":"GQ183900","term_id":"260408325","term_text":"GQ183900"}}GQ183900) that resulted 100% identical with that that was, in the meantime, annotated in the EquCab2 horse genome assembly derived from a Thoroughbred horse (http://www.ensembl.org/Equus_caballus/Search/, Ensembl release 52-Dec 2009). Our sequence contained 671 bp upstream from the ATG start codon, 538 bp of the promoter, and the entire 5′-untranslated region (UTR) of 133 bp, the three exons (except 33 bp of exon 1), the two intervening introns, and 80 bp of the 3′-UTR (Figure 2). The transcription start site of the first exon was deduced from human and bovine MSTN exon 1 sequences [4, 6]. The coding regions of exons 1, 2, and 3 of the horse MSTN gene contained 373, 374, and 381 bp, respectively. Introns 1 and 2 included 1829 bp and 2016 bp, respectively, almost the same length reported in cattle (1840 bp and 2033 bp, respectively) and pig (1809 bp and 1980 bp, respectively) [6, 8]. Intron 1 is a type 1 intron as it interrupts a codon between the first and second exon whereas intron 2 is a type 0 intron as it divides the coding sequence between two codons as in other species [6, 8]. The analysed proximal promoter region and the 5′-UTR of the horse MSTN exhibited a degree of identity with the corresponding regions of other species ranging from 77% (mouse) to 90% (pig). Putative consensus DNA sequences known as transcription factor binding sites, DNA-binding motifs, or cis-regulatory elements were identified in the positive strand of horse promoter (Figure 3). Considering the general transcription factors, three different putative TATA boxes (TATA-1, TATA-2, and TATA-3) and one CCAAT box were detected. Among muscle-specific transcription factors, four E-boxes (named E1, E2, E3, and E4 boxes, Figure 3), one putative site for myocyte specific enhancer factor 2 (MEF2 or MEB1) and consensus sequences for FoxO and SMAD binding sites (CAAAATA and CAGACA, respectively) family sequences were identified. The alignment of MSTN promoter sequences across different species (horse, cattle, goat, human, mouse, pig, and sheep) revealed that these DNA-binding motives, particularly close to the TATA-1 surrounding sequence, were highly conserved across species. In particular, TATA-1 was conserved in all examined species except mouse, the second TATA sequence (TATA-2) was conserved across all seven species, and TATA-3 was conserved in all species except pig and mouse. The MEF2 and E-boxes were conserved in all the considered mammals except in human and mouse for E-box4. The E-boxes can be activated by the myogenic regulatory factors (MRFs: MyoD, Myf5, myogenin, and MRF4). MyoD upregulates MSTN transcription [33] and at the same time MSTN inhibits MyoD expression and activity regulating the differentiation of myoblasts into myotubes [34]. MyoD and MRF4 play competitive roles in myogenesis and might act as molecular switches to determine myogenic differentiation and cell proliferation, respectively [35]. Additional E-boxes were identified in the analysed region (such as an E-box located near the TATA-2 in pig and an additional E-box in all mammalian but cattle) and in the distal region of the promoter of the other mammals (not included in Figure 3). In cattle, Spiller et al. [33] showed the importance of three functional E-boxes (E3, E4, and E6) of which the E6, occupied by MyoD in vitro and in vivo, resulted crucial for the MSTN promoter activity. The close position of functional E-boxes suggests that they might function as a cluster to better sustain the stability of DNA-protein. Across the MSTN promoter sequences of all considered livestock species we identified the conserved position of sites matching the consensus for FoxO binding and the adjacent SMAD box whose presence was not evidenced in previous works [8, 33, 36]. Recent data demonstrated that these factors appear to act through independent pathways but additively to regulate the expression of MSTN and contribute to control muscle cell growth and differentiation [37, 38]. In addition, FoxO transcription factors plays a critical role in development of muscle atrophy by stimulating proteolysis and by increasing myostatin expression. Putative E-boxes were identified both in intron 1 (six boxes) and in intron 2 (six boxes) and one putative E-box was located in the 3′-UTR at seven nucleotides downstream of the TGA stop codon (data not shown). The presence of E-boxes in the introns and 3′-UTR of equine MSTN gene has not been described yet even if their occurrence has been highlighted recently in introns of porcine MSTN gene [8].
Sequence analysis of the proximal promoter region and 5′-UTR of MSTN gene. Approximately 670 bp of the MSTN 5′-flanking region and 5′-UTR from horse, cattle, goat, human, mouse, pig, and sheep is shown. The consensus for transcription factor binding motives is shown in bold text and underlined. Numbering of the horse sequence is relative to the ATG start codon. Numbering of E-boxes is relative to the horse sequence. The two identified SNPs of horse MSTN gene are indicated with Y (T or C according to the IUPAC nomenclature) and in red.
3.2. Identification of Polymorphisms in the Horse MSTN Gene
Sequencing of the panel of horses of different morphological types revealed a total of seven single nucleotide polymorphisms (SNPs) (Figure 2). Two transitions were located in the promoter region at -646 ({"type":"entrez-nucleotide","attrs":{"text":"GQ183900","term_id":"260408325","term_text":"GQ183900"}}GQ183900:g.26T>C) and -516 ({"type":"entrez-nucleotide","attrs":{"text":"GQ183900","term_id":"260408325","term_text":"GQ183900"}}GQ183900:g.156T>C) bp upstream from the start codon. The g.26T>C SNP was within a conserved position (except in mouse) but not within an identified known functional motif while the g.156T>C polymorphism was within a TATA box-like (TATA-3; YATAAA, Figure 3). Sequence alignments of the MSTN promoter regions of different species indicate that the g.26T and g.156T alleles derive from an ancestral MSTN sequence as most close species present the indicated nucleotides (Figure 3 and data not shown). The other five SNPs were in intronic regions: four were localized in intron 1 and one in intron 2 (Figure 2). Three of the SNPs of intron 1 (g.1634T>G, g.2115A>G, and g.2327A>C) were also recently identified in Thoroughbred breeds [24]. One of which (g.2115A>G; indicated by [24] as g.66493737C>T) has been associated with sprinting ability and racing stamina in Thoroughbred horses [24]. The remaining SNPs were not reported by others and represent new polymorphisms of the horse MSTN gene. None of these intronic SNPs resided within splice sites or within particularly conserved sequence elements. No indels and synonymous or nonsynonymous substitutions were identified.
3.3. Analysis of Polymorphisms in Breeds with Different Morphological Types and Genetic Diversity Parameters
Allele frequencies for the two SNPs located in the promoter region (g.26T>C and g.156T>C) are shown in Table 3. The g.26T>C SNP was polymorphic in 6 out of 16 breeds with higher observed frequency of the g.26C allele in the Lipizzan breed (0.21). For the g.156T>C polymorphism, the mutant g.156C allele, which changes the predicted TATA box3-like, was detected in 11 out 16 breeds and was identified in homozygous condition in a few Bardigiano, Haflinger, Noric, Rapid Heavy Draft, and Uruguayan Creole horses.
Table 3
Minor allele frequency (MAF) and haplotype (Hap*) frequencies of the two promoter SNPs of the MSTN gene in different horse breeds.
| Horse breeds | Types§ | Number of horses | MAF g.26C | MAF g.156C | Hap [T:T] | Hap [T:C] | Hap [C:T] |
|---|---|---|---|---|---|---|---|
| Rapid Heavy Draft | B | 81 | 0.12 | 0.15 | 0.73 | 0.15 | 0.12 |
| Noric | B | 26 | 0.06 | 0.40 | 0.54 | 0.40 | 0.06 |
| Bardigiano | M | 34 | — | 0.15 | 0.85 | 0.15 | — |
| Haflinger | M | 31 | — | 0.37 | 0.63 | 0.37 | — |
| Lipizzan | M | 12 | 0.21 | — | 0.79 | — | 0.21 |
| Murgese | M | 12 | — | 0.17 | 0.83 | 0.17 | — |
| Tolfetano | M | 7 | — | 0.29 | 0.71 | 0.29 | — |
| Uruguayan Creole | M | 35 | — | 0.24 | 0.76 | 0.24 | — |
| Italian Saddle | M-D | 30 | 0.05 | 0.10 | 0.85 | 0.10 | 0.05 |
| Maremmano | M-D | 13 | — | 0.12 | 0.88 | 0.12 | — |
| Quarter Horse | M-D | 15 | — | — | 1.00 | — | — |
| Salernitano | M-D | 4 | — | 0.25 | — | — | — |
| Spanish Purebred | M-D | 10 | — | — | 1.00 | — | — |
| Ventasso | M-D | 8 | 0.12 | — | 0.88 | — | 0.12 |
| Italian Trotter | D | 67 | 0.01 | — | 0.99 | — | 0.01 |
| Thoroughbred | D | 11 | — | 0.05 | 0.95 | 0.05 | — |
Haplotype analyses of the two mutations showed the presence of three haplotypes: [g.26T:g.156T], [g.26T:g.156C], and [g.26C:g.156T] (Table 3). The [T:T] haplotype could be the wild type according to its presence in all breeds and higher frequency (from 0.54 to 1.00). The [T:C] haplotype was observed in 10 breeds (frequency from 0.05 to 0.40), whereas the [C:T] haplotype was identified only in 6 breeds (frequency from 0.01 to 0.21) (Table 3).
In order to evaluate if the two SNPs in the promoter region could account for a quote of variability related to morphological types, we classified the analysed horse breeds in four groups (brachymorphic, B; mesomorphic, M; meso-dolichomorphic, M-D; and dolichomorphic, D) (see Materials and Methods). Several descriptive statistics summarizing the genetic diversity of these groups are reported in Table 4. The B group showed the highest observed and expected heterozygosity (0.24 ± 0.09 and 0.26 ± 0.10, respectively), whereas the D group had the lowest values (0.02 ± 0.01 for both measures). For the g.26T>C SNP, differences in allele frequencies were significant between B and the other three groups (P < .0001 for B versus M and B versus D; P < .05 for B versus M-D). For the g.156T>C polymorphism, only the comparison between B and M groups was not significant. In particular, differences in allele frequencies were highly significant between the B and D groups and between the M and D groups (P = 3.78E − 11 and P = 3.95E − 9, respectively. For the remaining comparisons: P < .0001 for B versus M-D and for M versus M-D, P < .01 for M-D versus D). The overall F st value showed that the genetic differences among the groups accounted for 6.1% (3.6% for the g.26T>C SNP and 7.0% for the g.156T>C SNP) of the genetic variation. The AMOVA on haplotypes confirmed that a proportion of the total molecular variance was associated to morphological types of the horses. Using the four morphological types the molecular variance explained was 6.40% (P < .05). Grouping these four types into two groups (B + M and M-D + D) according to their similarities on morphological types the quote of explained molecular variance was 10.6% (P < .01). It could be possible that differences of allele and haplotype frequencies among types are influenced by phylogenetic closeness rather than any association with morphological types. This issue should be further investigated as, to our knowledge, there are no studies analyzing this question that include most of the breeds we investigated. However, Di Stasio et al. [39] analysed genetic relationships among only three breeds included in our study and evidenced significant genetic differentiation among Bardigiano, Haflinger, and Maremmano, suggesting that the results we obtained might not be biased by a putative common origin of the breeds. The association of the two promoter SNPs with morphological types could be due to linkage disequilibrium with alleles in other chromosome 18 loci that affect the variability of morphological traits in horses. However, based on our results it cannot be excluded that MSTN SNPs could influence morphological traits, that are indirectly related to muscle mass. A few SNPs in the promoter region of the swine MSTN gene were associated with muscularity, growth, and meat quality traits [8, 20, 21]. One of them, with high frequency in the muscled Belgian Pietrain breed, was associated with MSTN expression level, suggesting that promoter polymorphisms could contribute to muscle mass in this pig breed [8]. To demonstrate the putative functional role of the identified horse MSTN promoter SNPs, expression studies in skeletal muscle of animals with different genotypes should be performed. However, it is worth to point out that in vivo RNA expression studies in horses are very complicated as it is quite difficult to standardize temporary and permanent environmental factors (i.e. age, sex, management, feeding, etc.) that are major sources of variability in such experiments. For these reasons in vitro assays might be needed to clarify if the identified SNPs could alter MSTN gene expression. In addition association analysis in breeds segregating for the two promoter SNPs and for which estimated breeding values for several conformational and performance traits are available could be useful to further evaluate the association of these polymorphic sites with phenotypic traits.
Table 4
Minor allele frequency (MAF), observed and expected heterozygosity (Het*), and haplotype (Hap§) frequencies of the two promoter SNPs of the MSTN gene in horses with different morphological types.
| Types# | Number of horses | MAF g.26C | MAF g.156C | Observed Het (s.d.)§ | Expected Het (s.d.)§ | Hap [T:T] | Hap [T:C] | Hap [C:T] |
|---|---|---|---|---|---|---|---|---|
| B | 107 | 0.11 | 0.21 | 0.24 (0.09) | 0.26 (0.10) | 0.69 | 0.22 | 0.09 |
| M | 131 | 0.02 | 0.22 | 0.18 (0.20) | 0.19 (0.22) | 0.76 | 0.22 | 0.02 |
| M-D | 80 | 0.04 | 0.07 | 0.11 (0.03) | 0.11 (0.03) | 0.89 | 0.08 | 0.03 |
| D | 78 | 0.01 | 0.01 | 0.02 (0.01) | 0.02 (0.01) | 0.98 | 0.01 | 0.01 |
| Total | 396 | 0.05 | 0.15 | 0.15 (0.09) | 0.17 (0.11) | 0.81 | 0.15 | 0.04 |
Acknowledgments
The authors thank horse breeders, Italian Horse National Breeders Associations, Associazione Italiana Allevatori (A.I.A.), Istituto Sperimentale per la Zootecnia of Rome (Italy), and Association Rural de Uruguay for providing horse samples and genealogical information. They also thanks three anonymous reviewers for their comments that made it possible to improve the manuscript. This study was funded by the University of Bologna RFO.
References
1. McPherron AC, Lawler AM, Lee S-J. Regulation of skeletal muscle mass in mice by a new TGF-β superfamily member. Nature. 1997;387(6628):83–90. [PubMed] [Google Scholar]
2. McPherron AC, Lee S-J. Double muscling in cattle due to mutations in the myostatin gene. Proceedings of the National Academy of Sciences of the United States of America. 1997;94(23):12457–12461. [PMC free article] [PubMed] [Google Scholar]
3. Joulia-Ekaza D, Cabello G. Myostatin regulation of muscle development: molecular basis, natural mutations, physiopathological aspects. Experimental Cell Research. 2006;312(13):2401–2414. [PubMed] [Google Scholar]
4. Gonzalez-Cadavid NF, Taylor WE, Yarasheski K, et al. Organization of the human myostatin gene and expression in healthy men and HIV-infected men with muscle wasting. Proceedings of the National Academy of Sciences of the United States of America. 1998;95(25):14938–14943. [PMC free article] [PubMed] [Google Scholar]
5. Stratil A, Kopecny M. Genomic organization, sequence and polymorphism of the porcine myostatin (GDF8; MSTN) gene. Animal Genetics. 1999;30(6):468–469. [PubMed] [Google Scholar]
6. Jeanplong F, Sharma M, Somers WG, Bass JJ, Kambadur R. Genomic organization and neonatal expression of the bovine myostatin gene. Molecular and Cellular Biochemistry. 2001;220(1-2):31–37. [PubMed] [Google Scholar]
7. Du R, Chen Y-F, An X-R, et al. Cloning and sequence analysis of myostatin promoter in sheep. DNA Sequence. 2005;16(6):412–417. [PubMed] [Google Scholar]
8. Stinckens A, Luyten T, Bijttebier J, et al. Characterization of the complete porcine MSTN gene and expression levels in pig breeds differing in muscularity. Animal Genetics. 2008;39(6):586–596. [PubMed] [Google Scholar]
9. Schuelke M, Wagner KR, Stolz LE, et al. Myostatin mutation associated with gross muscle hypertrophy in a child. New England Journal of Medicine. 2004;350(26):2682–2688. [PubMed] [Google Scholar]
10. Kambadur R, Sharma M, Smith TPL, Bass JJ. Mutations in myostatin (GDF8) in double-muscled Belgian Blue and Piedmontese cattle. Genome Research. 1997;7(9):910–916. [PubMed] [Google Scholar]
11. Grobet L, Martin LJR, Poncelet D, et al. A deletion in the bovine myostatin gene causes the double-muscled phenotype in cattle. Nature Genetics. 1997;17(1):71–74. [PubMed] [Google Scholar]
12. Grobet L, Poncelet D, Royo LJ, et al. Molecular definition of an allelic series of mutations disrupting the myostatin function and causing double-muscling in cattle. Mammalian Genome. 1998;9(3):210–213. [PubMed] [Google Scholar]
13. Marchitelli C, Savarese MC, Crisà A, Nardone A, Marsan PA, Valentini A. Double muscling in Marchigiana beef breed is caused by a stop codon in the third exon of myostatin gene. Mammalian Genome. 2003;14(6):392–395. [PubMed] [Google Scholar]
14. Clop A, Marcq F, Takeda H, et al. A mutation creating a potential illegitimate microRNA target site in the myostatin gene affects muscularity in sheep. Nature Genetics. 2006;38(7):813–818. [PubMed] [Google Scholar]
15. Boman IA, Klemetsdal G, Blichfeldt T, Nafstad O, Våge DI. A frameshift mutation in the coding region of the myostatin gene (MSTN) affects carcass conformation and fatness in Norwegian White Sheep (Ovis aries) Animal Genetics. 2009;40(4):418–422. [PubMed] [Google Scholar]
16. Boman IA, Våge DI. An insertion in the coding region of the myostatin (MSTN) gene affects carcass conformation and fatness in the Norwegian Spælsau (Ovis aries) BMC Research Notes. 2009;2, article 98 [PMC free article] [PubMed] [Google Scholar]
17. Mosher DS, Quignon P, Bustamante CD, et al. A mutation in the myostatin gene increases muscle mass and enhances racing performance in heterozygote dogs. PLoS Genetics. 2007;3(5):779–786. [PMC free article] [PubMed] [Google Scholar]
18. Bellinge RHS, Liberles DA, Iaschi SPA, O'Brien PA, Tay GK. Myostatin and its implications on animal breeding: a review. Animal Genetics. 2005;36(1):1–6. [PubMed] [Google Scholar]
19. Gill JL, Bishop SC, McCorquodale C, Williams JL, Wiener P. Associations between the 11-bp deletion in the myostatin gene and carcass quality in Angus-sired cattle. Animal Genetics. 2009;40(1):97–100. [PubMed] [Google Scholar]
20. Guimaraes SEF, Stahl CH, Lonergan SM, Geiger B, Rothschild MF. Myostatin promoter analysis and expression pattern in pigs. Livestock Science. 2007;112(1-2):143–150. [Google Scholar]
21. Yu L, Tang H, Wang J, et al. Polymorphisms in the 5′ regulatory region of myostatin gene are associated with early growth traits in Yorkshire pigs. Science in China C. 2007;50(5):642–647. [PubMed] [Google Scholar]
22. Hosoyama T, Kawada S, Oshiumi R, et al. Molecular cloning of equine (thoroughbred) myostatin cDNA and detection of myostatin precursor proteins in the serum. Journal of Reproduction and Development. 2002;48(4):335–342. [Google Scholar]
23. Caetano AR, Pomp D, Murray JD, Bowling AT. Comparative mapping of 18 equine type I genes assigned by somatic cell hybrid analysis. Mammalian Genome. 1999;10(3):271–276. [PubMed] [Google Scholar]
24. Hill EW, Gu J, Eivers SS, et al. A sequence polymorphism in MSTN predicts sprinting ability and racing stamina in thoroughbred horses. PLoS ONE. 2010;5(1, article e8645) [PMC free article] [PubMed] [Google Scholar]
25. Bowling AT, Ruvinsky A. Genetic aspect of domestication, breeds and their origins. In: Bowling AT, Ruvinsky A, editors. The Genetics of the Horse. Oxfordshire, UK: CAB International; 2000. [Google Scholar]
26. Aparicio G. Zootecnia especial. Etnologia compendiada. Deposito Legal. 1960;16:25–27. [Google Scholar]
27. Falaschini A, Rizzi S, Pasquini M. Morphological evolution of the Haflinger horse. Italian Journal of Animal Science. 2003;2(1):595–597. [Google Scholar]
28. Sabbioni A, Beretti V, Zanon A, et al. Morphological evolution of Bardigiano horse. Italian Journal of Animal Science. 2005;4(2):412–414. [Google Scholar]
29. Dario C, Carnicella D, Dario M, Bufano G. Morphological evolution and heritability estimates for some biometric traits in the Murgese horse breed. Genetics and Molecular Research. 2006;5(2):309–314. [PubMed] [Google Scholar]
30. Falaschini A, Gallinelli C. Morphological evolution of the Throroughbred reared in Italy. Annual Meeting of the Italian Society for Veterinary Sciences. 2006;60:497–498. [Google Scholar]
31. Yeh FC, Yang R, Boyle T. POPGENE Version 1.31. Microsoft Windows-based freeware for population genetic analysis. 1999, http://www.ualberta.ca/~fyeh/
32. Stephens M, Smith NJ, Donnelly P. A new statistical method for haplotype reconstruction from population data. American Journal of Human Genetics. 2001;68(4):978–989. [PMC free article] [PubMed] [Google Scholar]
33. Spiller MP, Kambadur R, Jeanplong F, et al. The myostatin gene is a downstream target gene of basic helix-loop-helix transcription factor MyoD. Molecular and Cellular Biology. 2002;22(20):7066–7082. [PMC free article] [PubMed] [Google Scholar]
34. Langley B, Thomas M, Bishop A, Sharma M, Gilmour S, Kambadur R. Myostatin inhibits myoblast differentiation by down-regulating MyoD expression. Journal of Biological Chemistry. 2002;277(51):49831–49840. [PubMed] [Google Scholar]
35. Jin X, Kim J-G, Oh M-J, et al. Opposite roles of MRF4 and MyoD in cell proliferation and myogenic differentiation. Biochemical and Biophysical Research Communications. 2007;364(3):476–482. [PubMed] [Google Scholar]
36. Ma K, Mallidis C, Artaza J, Taylor W, Gonzalez-Cadavid N, Bhasin S. Characterization of 5′-regulatory region of human myostatin gene: regulation by dexamethasone in vitro. American Journal of Physiology. 2001;281(6):E1128–E1136. [PubMed] [Google Scholar]
37. Allen DL, Unterman TG. Regulation of myostatin expression and myoblast differentiation by FoxO and SMAD transcription factors. American Journal of Physiology. 2007;292(1):C188–C199. [PubMed] [Google Scholar]
38. Allen DL, Du M. Comparative functional analysis of the cow and mouse myostatin genes reveals novel regulatory elements in their upstream promoter regions. Comparative Biochemistry and Physiology. 2008;150(4):432–439. [PubMed] [Google Scholar]
39. Di Stasio L, Perrotta G, Blasi M, Lisa C. Genetic characterization of the Bardigiano horse using microsatellite markers. Italian Journal of Animal Science. 2008;7(2):243–250. [Google Scholar]
Articles from Journal of Biomedicine and Biotechnology are provided here courtesy of Hindawi Limited
Source: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2913906/
Postar um comentário for "What Type of Mutation Occurred in the Belgian Blue Myostatin"