Bioinformatic analysis based on the complete coding region of the MSTN gene within and among different species

Myostatin, encoded by the MSTN gene (previously GDF8), is a member of the transforming growth factor-β superfamily, which normally acts to limit skeletal muscle mass by regulating the number and growth of muscle fibers. In this study, a total of 84 myostatin gene sequences with known complete coding regions (CDS) and corresponding amino acid sequences were analyzed from 17 species, and differentiation within and among species was studied using comparative genomics and bioinformatics. Characteristics of the nucleotide and amino acid sequences were also predicted. The results indicated that a total of 569 polymorphic sites, including 53 singleton variable sites and 516 parsimony informative sites, which could be sorted into 44 haplotypes, were detected from 17 species. Observed genetic diversity was higher among species than within species, and Vulpes lagopus was more polymorphic than other species. There was clear differentiation of the myostatin gene among species and the reconstructed phylogenetic tree was consistent with the NCBI taxonomy. The myostatin gene was 375-aa long in most species, except for Mus musculus (376 aa) and Danio rerio (373 aa). The amino acid sequences of myostatin were deemed hydrophilic, and had theoretical pI values of <7.0, mostly due to the acidic polypeptide. The instability index of the myostatin protein was 40.48-51.63, indicating that the polypeptide is not stable. The G+C content of the CDS nucleotide sequence in different species was 40.60-51.69%. The predicted promoter region of the Ovis aries myostatin gene was 150-220 bp upstream of the start codon.