Sequence and polymorphism analysis of the camel
(Camelus dromedarius) myostatin gene

Stefania Muzzachi, Ahmad Oulmouden, Youcef Cherifi, Mohamed Habib

Yahyaoui, Mohamed Zayed, Pamela Burger, Giovanni Michele Lacalandra,

Bernard Faye, Elena Ciani

► To cite this version:

Stefania Muzzachi, Ahmad Oulmouden, Youcef Cherifi, Mohamed Habib Yahyaoui, Mohamed Zayed, et al.. Sequence and polymorphism analysis of the camel (Camelus dromedarius) myostatin gene. Emirates Journal of Food and Agriculture, 2015, 27 (4), pp.367-373. 10.9755/ejfa.v27i4.19910. hal-01211990

HAL Id: hal-01211990

https://hal.science/hal-01211990v1

Submitted on 27 May 2020

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
REGULAR ARTICLE

Sequence and polymorphism analysis of the camel (*Camelus dromedarius*)
myostatin gene

S. Muzzachi<fnref n="1" />, A. Oulmouden<fnref n="2" />, Y. Cherifi<fnref n="3" />, H. Yahyaoui<fnref n="4" />, M. A. Zayed<fnref n="5" />, P. Burger<fnref n="6" />, G. M. Lacalandra<fnref n="7" />, B. Faye<fnref n="8" />,<fnref n="9" />
and E. Ciani<fnref n="1" />*

Abstract

Myostatin (MSTN), a negative regulator of skeletal muscle development in mammals, represents a key target for genetic investigations in meat-producing animals, with mutations responsible for increased skeletal-muscle mass currently described in several livestock species. Dromedary camels play a major economic role as suppliers of meat for human consumption across several countries. Notwithstanding, a comprehensive characterization of the sequence variability at the *Camelus dromedarius* MSTN locus was still lacking. Here we present the first extensive sequence and polymorphism analysis of the MSTN gene in the *C. dromedarius* species. Out of more than 3.6 kb of nucleotide sequence screened on 22 animals from 3 different Northern African regions, only 3 variant sites in the first intron were detected. The low observed diversity may reflect the evolutionary history of the species, likely developed as domesticates from a low variable wild ancestor population. Sequence identity among *C. dromedarius* and other *Cetartiodactyla* highlighted a tree topology consistent with previous reports of a closer relationship between *Tylopoda* and *Suiformes*. A close similarity between *C. ferus* and *C. dromedarius* was observed within *Tylopoda*. A markedly higher sequence identity between *C. dromedarius* and the other vertebrate species was observed at the MSTN locus compared to other genes, thus confirming it as a highly conserved target across mammals.

Key words: *Camelus dromedarius*, Myostatin, Sequence, Polymorphism

Introduction

Myostatin (MSTN), a TGF-β family member,
is a negative regulator of skeletal muscle
development and growth in mammals. The gene

sequence has been shown to be highly conserved
across vertebrate species (McPherron and Lee,
1997). Loss-of-function mutations associated with
increased skeletal-muscle mass and other
pleiotropic effects have been detected in several
species, including cattle (Grobet et al., 1997), mice
(Szabò et al., 1998), humans (Schuelke et al.,
2004), sheep (Clop et al., 2006), dogs (Mosher et
al., 2007), chickens (Ye et al., 2007; McFarland et
al., 2007; Yang et al., 2003) and putatively pigs
(Stinckens et al., 2008). These evidences support
the hypothesis that also MSTN function could be
highly conserved among vertebrates. A total of
eight replacement changes at the myostatin gene

Received 05 May 2014; Revised 21 November 2014; Accepted
10 December 2014; Published Online 01 April 2015

Elena Ciani
Department of Biosciences, Biotechnologies, Biopharmaceutics,
Viale Amendola 165/ 70126 Bari, Italy

Email: elena.ciani@uniba.it

<footnote n="1">¹Department of Biosciences, Biotechnologies, Biopharmaceutics, Viale Amendola 165/ 70126 Bari, Italy</footnote>

<footnote n="2">²Unité de Génétique Moléculaire Animale, Université de Limoges Faculté des Sciences et Techniques, INRA UMR 1061,
123 Avenue Albert Thomas, 87060 Limoges, France</footnote>

<footnote n="3">³Laboratoire de Génétique Moléculaire et Cellulaire, Faculté des Sciences de la Nature et de la Vie, Département de Génétique Moléculaire Appliquée. Université des Sciences et de la Technologie d'Oran Mohamed Boudiaf BP 1505 El Mnaouer, Oran, Algerie</footnote>

<footnote n="4">⁴Laboratoire Elevage et Faune Sauvage, Institut des Régions Arides, km 22.5 Route Djorf 4119 Medenine, Tunisie</footnote>

<footnote n="5">⁵Meat Technology Unit, Animal Breeding Department, Desert Research Center, 1 Mathaf El-Matarya St., El-Matarya, Cairo, Egypt</footnote>

<footnote n="6">⁶Institut für Populationsgenetik, Vetmeduni Vienna, Veterinärplatz 1, 1210 Wien, Austria</footnote>

<footnote n="7">⁷Department of Emergency and Organ Transplantation (D.E.T.O.), Veterinary Clinics and Animal Production Section,
University of Bari, Valenzano, Bari, Italy</footnote>

<footnote n="8">⁸CIRAD-ES, Campus International de Baillarguet, 34398 Montpellier Cedex 1, France</footnote>

<footnote n="9">⁹ FAO/MOA Camel project, P.O.Box 761Al-Kharj, 11942 Saudi Arabia</footnote>

<footnote>*Corresponding Author</footnote>
have been reported in humans (Saunders et al., 2006) while in cattle, one of the most extensively surveyed species, 20 different mutations have been found, 14 of them being located in the coding sequence and six in the intronic regions (for a review, see Stinckens et al., 2011; Rodgers and Garikipati, 2008). In humans, the myostatin gene has been reported to be under the effect of strong positive diversifying selection, as supported by the excess of non-synonymous polymorphisms, to an extent previously observed only for a few other genes, like G6PD and the major histocompatibility complex. Moreover, the frequency patterns observed at two amino acid variants suggested the detected selection signature as a recent phenomenon (Saunders et al., 2006). A recent phenomenon has been postulated also for cattle, where homozygous segments associated to the different MSTN alleles resulted to be longer than those observed in the wild-type counterpart; indeed, the estimate that these alleles had a recent common ancestor less than 400 years ago, together with the evidence of their increasing frequency or fixation have been interpreted as a proof of recent anthropic selection in cattle. A selection pressure on MSTN has been postulated also in bovid where, instead, the observed accelerated evolution of the gene sequence has been explained as the consequence of remote naturally occurring selection (Tellgren et al., 2004). Finally, putative naturally occurring selection has been evoked in salmonids (Ostbye et al., 2007). Based on the above, we decided to partially sequence the MSTN gene in an economical relevant livestock species, *Camelus dromedarius*, for which only a 256 bp region in the first exon of the MSTN locus was available up to now. Moreover, a preliminary sequence polymorphism analysis was carried out in a population sample of 22 animals from three different Northern African geographic regions.

## Materials and Methods

### Sample collection

A total of 22 *Camelus dromedarius* blood samples from Tunisia (5) Egypt (7) and Algeria (10) were collected from different farms in different areas in order to increase the probability to capture a higher genetic variability at the MSTN locus.

### DNA extraction from blood and muscle

Genomic DNA was extracted from blood with FlexiGene DNA Kit (Qiagen), according to the manufacturer's instructions. Briefly 2 ml of blood were mixed with buffer FG1 and centrifuged for 5 min at 2000g. After supernatant was discarded and buffer FG2 /protease mixture was added. Samples

were incubated at 65°C for 10 min. After incubation, isopropanol was added for DNA precipitation and samples were centrifuged for 3 min at 2000g. Then ethanol was added and samples were centrifuged for 3 min at 2000g. Supernatant was eliminated and DNA was dissolved with buffer FG3 and incubated for 1 h at 65°C in a water bath. Genomic DNA concentration was evaluated with the Nanodrop 2000C spectrophotometer (Thermo Scientific). Samples were stored a -20°C until use.

## PCR amplification

Heterologous oligonucleotide primers to amplify the *Camelus dromedarius* MSTN locus were designed as follows. First, BLAST analysis using the *Vicugna pacos* MSTN sequence, available at the Ensemble database (www.ensembl.org) was carried out, which returned, for the *Camelus ferus* species, a contig (GenBank Accession No AGVR01040332.1) and a scaffold (GenBank Accession No NW_006211358.1), both including the whole-gene MSTN sequence. Alignment of the above sequences was carried out using ClustalW (www.genome.jp/tools/clustalw/) in order to *in silico* detect possible variant sites. Primers were then designed on the *Camelus ferus* consensus sequence using Primer3 (www.primer3.ut.ee/, version 4.0.0). Primers and PCR conditions for each target region are provided below.

**Exon 1** – The sequence of the forward primer was 5'CCTTGGCATTACTCAAAAGCAA3'; the sequence of the reverse primer was 5'CTTACATACACACGAACAGCC3'. The PCR thermal profile consisted of an initial step at 95°C for 15 min, followed by 35 cycles of 30 s at 94°C, 90 s at 57°C, 90 s at 72°C; the final extension step was carried out at 72°C for 10 min. The sequence also covers putatively 58 bp of the 5'UTR region and 91 bp of the first intron.

**Intron 1** (partial) - The sequence of the forward primer was 5'TTTGGGCATTTGCTGACAACCTAGAATGA C3'; the sequence of the reverse primer was 5'CTTTGGGAATCTGAGTAGTTACACTTTACT G3'. The PCR thermal profile consisted of an initial step at 95°C for 15 min, followed by 35 cycles of 30 s at 94°C, 90 s at 61°C, 90 s at 72°C; the final extension step was carried out at 72°C for 10 min.

**Exon 2** - The sequence of the forward primer was 5'ATCAGATAGTCCTGGAGTAGATCTGCCTT A3'; the sequence of the reverse primer was 5'GTCTGGCTTATGAGCATAGGAAACTGGTA G3'. The PCR thermal profile consisted of an initial step at 95°C for 15 min, followed by 35 cycles of
30 s at 94°C, 90 s at 61°C, 90 s at 72°C; the final extension step was carried out at 72°C for 10 min. The sequence also covers putatively 717 bp of the first intron and 111 bp of the second intron.

*Intron 2* (partial) - The sequence of the forward primer was 5'TTACTCTTCTGGCTTACCTGTGCTTGTGT3'; the sequence of the reverse primer was 5'AAATTTGTCACATCACCAACTTTAATGAT T3'. The PCR thermal profile consisted of an initial step at 95°C for 15 min, followed by 35 cycles of 30 s at 94°C, 90 s at 61°C, 90 s at 72°C; the final extension step was carried out at 72°C for 10 min.

*Exon 3* - The sequence of the forward primer was 5'AGCCATAAGGGAAGAATCAAGCC3'; the sequence of the reverse primer was 5'GATCTCCTAGCTCGAAAACTTAGG3'. The PCR thermal profile consisted of an initial step at 95°C for 15 min, followed by 35 cycles of 30 s at 94°C, 90 s at 58°C, 90 s at 72°C; the final extension step was carried out at 72°C for 10 min. The sequence also covers putatively 154 bp of the second intron and 212 bp of the 3'-UTR region.

For all the samples, DNA amplification was carried out in a reaction mixture at a final volume of 25 µl, consisting of Multiplex PCR Master Mix 1X (Qiagen), oligonucleotide primers 0.2 µM (MWG), DNA template 2 ng/µl, and RNAse free-water.

## Amplicon purification and sequence analysis

PCR products were separated on a 2% agarose gel electrophoresis at 90V for 1 hour. The fragments were then excised from the agarose gel and purified by using the QIAquick Gel Extraction Kit (Qiagen) according to the manufacturer's instructions. Briefly, DNA was excised from the agarose gel and 3 volumes of buffer QG were added to 1 volume of gel. Gel slices were incubated at 50°C for 10 min. After incubation, isopropanol was added to the sample and the mix was transferred to the column and centrifuged for 1 min at 10000g. For each column 0.5 ml of buffer QG was added and centrifuged for 1 min at 10000g. A washing step was carried out adding 0.75 ml of buffer PE and centrifuging for 1 min at 10000g. The flow-through was discarded and the column was centrifuged for an additional 1 min at 10000g. DNA elution was carried out adding 30 µl of water to the membrane and centrifuging the column 1 min at 10000g. The concentration of the purified PCR products was evaluated with the Nanodrop 2000C spectrophotometer (Thermo Scientific). Sanger sequencing of the purified fragments was performed by the Eurofins MWG Operon service.

Detection of sequence polymorphisms and differentially fixed nucleotide sites was carried out through alignment of *Camelus dromedarius* sequences against the *Camelus ferus* contig sequence (GenBank Accession No AGVR01040332.1). Both sequence alignment and pair-wise sequence identity percentages among *Camelus dromedarius* and other species belonging to the *Cetartiodactyla* clade were obtained, for the MSTN locus, using ClustalW2 (Larkin et al., 2007). The software Splitstree (Huson and Bryant, 2006) was used to construct a NJ tree from the pair-wise sequence identity matrix.

## Results and Discussion

In the present study we report for the first time the nucleotide sequence of the three exons of the MSTN locus for the *Camelus dromedarius* species obtained after *in vitro* amplification of genomic DNA using heterologous oligonucleotide primers designed from the publicly available sequence of the *Vicugna pacos* myostatin gene. Previously, only a 256 bp region in the first exon of the *Camelus dromedarius* MSTN locus had been published (Shah et al., 2006). Here we sequenced more than 3.6 kb of nucleotide sequence, including the three exons, part of intron 1 and intron 2 and part of the 3' and 5' ends of the *Camelus dromedarius* myostatin gene. The NJ tree obtained from the pair-wise sequence identity among *Camelus dromedarius* and other species belonging to the *Cetartiodactyla* clade is provided in Figure 1. Since coding and non-coding regions of a genome evolve at different rates and alignments of non-coding regions are often difficult to interpret, we used here pair-wise sequence identity values estimated after removal of non-exonic sequences. The tree topology was found to be sound and in line with previous reports such that of a closer relationship between *Tylopoda* and *Suiformes* previously described, among others, by Kacskovics et al. (2006). As what concerns *Tylopoda*, a close similarity between *C. ferus* and *C. dromedarius* was observed, with a sequence identity of 99.73%, only lower to that observed among *Orcinus orca* and *Tursiops truncatu* (99.91%). As expected, a slightly lower identity than that observed between *C. dromedarius* and *C. ferus* (Old World camelids) was observed between *C. dromedarius* and the New World camelid, *Vicugna pacos* (Table 2). The lowest value was observed between *C. dromedarius* vs. *Bubalus bubalis* (93.59%), followed by Bos taurus (93.68%), *Capra hircus* and *Ovis aries* (94.04% both), thus supporting the evidence that, despite being generally assimilated to ruminants,
camelids have peculiar evolutionary (and physiological) features. Interestingly, when compared to the available literature (Dahiya et al., 2014; Nagarajan et al., 2012, 2011; Du et al., 2009; Kacskovics et al., 2006), the sequence identity values observed for the MSTN gene between *C. dromedarius* and the other vertebrate species looks markedly higher, thus confirming previous evidences of the myostatin gene as a highly conserved gene across mammals. Further, to characterize the *Camelus dromedarius* MSTN gene in terms of sequence polymorphism, we sequenced a total of 22 animals from three different geographic regions (Algeria, Tunisia and Egypt). The results of the polymorphism analysis are represented in Table 1. As can be observed, we only detected three variant nucleotide sites, notably two transitions (798_G/A and 799_C/T) and one transversion (486_G/C). All the polymorphic sites were located in the first intron. Shah et al. (2006) screened a 256 bp region in the first exon of the *Camelus dromedarius* MSTN locus in 12 samples from six different Pakistani breeds without observing any sequence polymorphism. On the contrary, they found a C→T transition responsible for an Alanine to Valine amino acidic substitution when a 422 bp region of the first exon of the MYF5 gene was screened in the same animals. Vaccarelli et al. (2012) detected, at the T-cell receptor (TCR) gamma genes, SNPs with a frequency of one every 114 to 485 nucleotides, depending on the considered rearrangements, and somatic hypermutation was suggested on that locus as a mechanism to diversify the *Camelus dromedarius*

TCR gene repertoire face to the wide spectrum of antigenic determinants presented by a variety of diverse and evolving pathogens. High nucleotide sequence diversity as a consequence of somatic hypermutation had been previously reported also for the T-cell receptor delta chains (Antonacci et al., 2011). We did not observe any polymorphism in the second exon of the myostatin, which is known from other species to harbor functional mutations (Baron et al., 2012). The low diversity observed at the MSTN locus in *Camelus dromedarius* may reflect the evolutionary history of this species, which likely developed as domesticates from a low variable wild ancestor population. Indeed, a recent world-wide survey of genetic variability in the *Camelus dromedarius* species carried out adopting both maternally inherited mtDNA and bi-parentally inherited STR markers (Pamela Burger, unpublished data) revealed, consistently with previous published data from STR markers, low levels of genetic variation compared to other livestock species, which has been supposed to derive from the limited geographical distribution of the wild ancestor on the Arabian Peninsula and to the brief co-existence of wild and early domesticated individuals. The extensive dispersal of dromedary camels out of the Arabian Peninsula and the recurrent gene flow phenomena at ancient trading centers along the routes connecting the Arabian Peninsula with the Sahara would also explain the substantial homogeneity observed among our Northern African samples.

<figure><img src="image_2.png" /><figcaption>Figure 1. NJ tree obtained from pair-wise sequence identity at the three myostatin exons among *Camelus dromedarius* and other species belonging to the *Cetartiodactyla* clade.</figcaption></figure>
<table><caption>Table 1. Results from the polymorphism analysis at the MSTN locus carried out on 22 *Camelus dromedarius* samples from three different geographic regions (Algeria, Tunisia and Egypt).</caption><thead><tr><th>Target region</th><th>Size (bp)</th><th>N</th><th>Polymorphism<sup>a</sup></th></tr></thead><tbody><tr><td>5' end</td><td>80</td><td>17</td><td>-</td></tr><tr><td>Exon 1</td><td>373</td><td>17</td><td>-</td></tr><tr><td>Intron 1</td><td>1466</td><td>13</td><td>486_G/C<br>798_G/A<br>799_C/T</td></tr><tr><td>Exon 2</td><td>378</td><td>18</td><td>-</td></tr><tr><td>Intron 2</td><td>722</td><td>18</td><td>-</td></tr><tr><td>Exon 3</td><td>381</td><td>17</td><td>-</td></tr><tr><td>3' end</td><td>236</td><td>17</td><td>-</td></tr></tbody></table>

<table><caption>Table 2. Matrix of pair-wise percentage sequence identity estimated for the three *Camelus dromedarius* MSTN exons.</caption><thead><tr><th>Species</th><th>C. dromedarius</th><th>C. ferus</th><th>V. pacos</th><th>S. scrofa</th><th>T. truncates</th><th>O. orca</th><th>B. taurus</th><th>B. bubalis</th><th>C. hircus</th><th>O. aries</th></tr></thead><tbody><tr><td>C. dromedarius</td><td>100.00</td><td>99.73</td><td>99.29</td><td>97.42</td><td>97.06</td><td>96.97</td><td>93.68</td><td>93.59</td><td>94.04</td><td>94.04</td></tr><tr><td>C. ferus</td><td></td><td>100.00</td><td>99.55</td><td>97.69</td><td>97.33</td><td>97.24</td><td>93.95</td><td>93.86</td><td>94.31</td><td>94.30</td></tr><tr><td>V. pacos</td><td></td><td></td><td>100.00</td><td>97.78</td><td>97.59</td><td>97.50</td><td>94.22</td><td>94.13</td><td>94.57</td><td>94.57</td></tr><tr><td>S. scrofa</td><td></td><td></td><td></td><td>100.00</td><td>98.05</td><td>97.96</td><td>94.49</td><td>94.41</td><td>95.03</td><td>94.85</td></tr><tr><td>T. truncatus</td><td></td><td></td><td></td><td></td><td>100.00</td><td>99.91</td><td>95.91</td><td>95.56</td><td>96.00</td><td>96.00</td></tr><tr><td>O. orca</td><td></td><td></td><td></td><td></td><td></td><td>100.00</td><td>95.83</td><td>95.47</td><td>95.91</td><td>95.91</td></tr><tr><td>B. taurus</td><td></td><td></td><td></td><td></td><td></td><td></td><td>100.00</td><td>99.02</td><td>96.53</td><td>96.52</td></tr><tr><td>B. bubalis</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>100.00</td><td>96.44</td><td>96.44</td></tr><tr><td>C. hircus</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>100.00</td><td>99.65</td></tr><tr><td>O. aries</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>100.00</td></tr></tbody></table>

## Conclusions

This work represented the first extensive characterization of the myostatin gene in the dromedary camel. Interestingly, a very low diversity was observed at the MSTN locus in our Northern African population sample, which may reflect the peculiar evolutionary history of this species. Further ongoing studies from recent available whole genome sequences on the three *Camelus* species (*C. bactrianus*, *C. ferus* and *C. dromedarius*) will help to better clarify the MSTN patterns of evolution within the interesting *Tylopoda* genus.

## Acknowledgements

This study has been achieved within the PROCAMED project funded by European Union within the ENPI-CBC-MED, reference number I.B/1.1/493. The content of the present document is

under the responsibility of the PROCAMED partners and could not be considered as the position of European Union.

## Author Contributions

Study conception and design: A. O., E. C. Acquisition of data: Y. C., H. Y., A. Z., B. F. Analysis and interpretation of data: S. M., A. O. Drafting of manuscript: S. M., E. C. Critical revision: H. Y., P. B., G. M. L., B. F., E. C.

## References

* Antonacci R., M. Mineccia, M. P. Lefranc, H. M. Ashmaoui, C. Lanave, B. Piccinni, G. Pesole, M. S. Hassanane, S. Massari and S. Ciccarese. 2011. Expression and genomic analyses of *Camelus dromedarius* T cell receptor delta (TRD) genes reveal a variable domain repertoire enlargement due to CDR3

<footnote><sup>a</sup>Positions are defined in relation to the first nucleotide in the forward primer.</footnote>
diversification and somatic mutation. Mol. Immunol. 48(12-13):1384-1396.

Baron E. E., M. S. Lopes, D. Mendonça and A. da Câmara Machado. 2012. SNP identification and polymorphism analysis in exon 2 of the horse myostatin gene. Anim. Genet. 43(2):229-232.

Clop, A., F. Marcq, H. Takeda, D. Pirottin, X. Tordoir, B. Bibé, J. Bouix, F. Caiment, J. M. Elsen, F. Eychenne, C. Larzul, E. Laville, F. Meish, D. Milenkovic, J. Tobin, C. Charlier and M. Georges. 2006. A mutation creating a potential illegitimate microRNA target site in the myostatin gene affects muscularity in sheep. Nat. Genet. 38(7):813-818.

Dahiya, S. S., G. Nagarajan, V. K. Bharti, S. K. Swami, S. C. Mehta, F. C. Tuteja, S. D. Narnaware and N. V. Patil. 2014. Sequence analysis of the Toll-like receptor 2 gene of old world camels. J. Adv. Res. 5(6):695-704.

Du, J., S. Gao, H. Chang, G. Cong, T. Lin, J. Shao, Z. Liu, X. Liu and X. Cai. 2009. Bactrian camel (*Camelus bactrianus*) integrins avb3 and avb6 as FMDV receptors: Molecular cloning, sequence analysis and comparison with other species. Veteri. Immunol. Immunopathol. 131:190-199.

Grobet, L., L. J. Martin, D. Poncelet, D. Pirottin, B. Brouwers, J. Riquet, A. Schoeberlein, S. Dunner, F. Menissier, J. Massabanda, R. Fries, R. Hanset and M. Georges. 1997. A deletion in the bovine myostatin gene causes the double-muscled phenotype in cattle. Nat. Genet. 17:71-74.

Huson, D. H. and D. Bryant. 2006. Application of phylogenetic networks in evolutionary studies. Mol. Biol. Evol. 23(2):254-267.

Kacskovics, I., B. Mayer, Z. Kis, L. V. Frenyó, Y. Zhao, S. Muyldermans and L. Hammarström. 2006. Cloning and characterization of the dromedary (*Camelus dromedarius*) neonatal Fc receptor (drFcRn). 2006. Dev. Comp. Immunol. 30(12):1203-1215.

Larkin, M. A., G. Blackshields, N. P. Brown, R. Chenna, P. A. McGettigan, H. McWilliam, F. Valentin, I. M. Wallace, A. Wilm, R. Lopez, J. D. Thompson, T. J. Gibson and D. G. Higgins. 2007. ClustalW and ClustalX version 2. Bioinformat. 23(21):2947-2948.

McFarland, D. C., S. G. Velleman, J. E. Pessall, C.

Liu. 2007. The role of myostatin in chicken (*Gallus domesticus*) myogenic satellite cell proliferation and differentiation. Gen. Comp. Endocrinol. 151(3):351-357.

McPherron, A. C. and S. J. Lee. 1997. Double muscling in cattle due to mutations in the myostatin gene. Proc. Natl. Acad. Sci. USA. 94:12457-12461.

Mosher, D. S., P. Quignon, C. D. Bustamante, N. B. Sutter, C. S. Mellersh, H. G. Parker and E. A. Ostrander. 2007. A mutation in the myostatin gene increases muscle mass and enhances racing performance in heterozygote dogs. PLoS Genet. 25:3(5):e79.

Nagarajan, G., S. K. Swami, S. K. Ghorui, K. M. Pathak, R. K. Singh and N. V. Patil. 2011. Cloning and phylogenetic analysis of Interleukin-6 (IL-6) and tumor necrosis factor-α (TNF-α) from Indian dromedaries (*Camelus dromedarius*). Comp. Immunol. Microbiol. Dis. 34(3):291-298.

Nagarajan, G., S. K. Swami, S. K. Ghorui, K. M. Pathak, R. K. Singh and N. V. Patil. 2012. Cloning and sequence analysis of IL-2, IL-4 and IFN-γ from Indian Dromedary camels (*Camelus dromedarius*). Res. Vet. Sci. 92(3):420-426.

Ostbye, T. K., O. F. Wetten, A. Tooming-Klunderud, K. S. Jakobsen, A. Yafe, S. Etzioni, T. Moen and O. Andersen. 2007. Myostatin (MSTN) gene duplications in Atlantic salmon (*Salmo salar*): evidence for different selective pressure on teleost MSTN-1 and -2. Gene 403:159-169.

Rodgers, B. D., and D. K. Garikipati. 2008. Clinical, agricultural, and evolutionary biology of myostatin: A comparative review. Endocr. Rev. 29(5):513-534.

Saunders, M. A., J. M. Good, E. C. Lawrence, R. E. Ferrel, W. H. Li and M. W. Nachman. 2006. Human adaptive evolution at myostatin (GDF8), a regulator of muscle growth. Am. J. Hum. Genet. 79(6):1089-1097.

Schuelke, M., K. R. Wagner, L. E. Stolz, C. Hubner, T. Riebel, W. Komen, T. Braun, J. F. Tobin and S. J. Lee. 2004. Myostatin mutation associated with gross muscle hypertrophy in a child. New Engl. J. Med. 350:2682-2688.

Shah, M. G., A. S. Qureshi, M. Reissmann and H. J. Schwartz. 2006. Sequencing and sequence
analysis of myostatin gene in the exon 1 of the Camel (*Camelus dromedarius*). Pak. Vet. J. 26(4):176-178.

Stinckens, A., T. Luyten, J. Bijttebier, K. Van den Maagdenberg, D. Dieltiens, S. Janssens, S. De Smet, M. Georges and N. Buys. 2008. Characterization of the complete porcine MSTN gene and expression levels in pig breeds differing in muscularity. Anim. Genet. 39:586-596.

Stinckens, A., M. Georges and N. Buys. 2011. Mutations in the myostatin gene leading to hypermuscularity in mammals: indications for a similar mechanism in fish? Anim. Genet. 42(3):229-234.

Szabo, G., G. Dallmann, G. Muller, L. Patthy, M. Soller and L. Varga. 1998. A deletion in the myostatin gene causes the compact (Cmpt) hypermuscular mutation in mice. Mamm. Gen. 9:671-672.

Tellgren, A., A. C. Berglund, P. Savolainen, C. M. Janis and D. A. Liberles. 2004. Myostatin rapid sequence evolution in ruminants

predates domestication. Mol. Phylogenet. Evol. 33(3):782-790.

Vaccarelli, G., R. Antonacci, G. Tasco, F. Yang, L. Giordano, H. M. El Ashmaoui, M. S. Hassanane, S. Massari, R. Casadio and S. Ciccarese. 2012. Generation of diversity by somatic mutation in the *Camelus dromedarius* T-cell receptor gamma variable domains. Eur. J. Immunol. 42(12):3416-3428.

Yang, X. Y., J. Hou, X. R. An, H. Guan, K. M. Gou, S. H. Yang, L. D. Chen and Y. F. Chen. 2003. Molecular cloning, expression mutation of myostatin and study on biochemical activity of its C-terminal peptide. Sheng wu gong cheng xue bao. Chin. J. Biotechnol. 19(4):480-483.

Ye, X., S. R. Brown, K. Nones, L. L. Coutinho, J. C. Dekkers and S. J. Lamont. 2007. Associations of myostatin gene polymorphisms with performance and mortality traits in broiler chickens. Genet. Sel. Evol. 39:73-89.