Myostatin rapid sequence evolution in ruminants predates domestication

Beyond the Big Five: Investigating Myostatin Structure, Polymorphism and Expression in Camelus dromedarius

Myostatin, a negative regulator of skeletal muscle mass in animals, has been shown to play a role in determining muscular hypertrophy in several livestock species, and a high degree of polymorphism has been previously reported for this gene in humans and cattle. In this study, we provide a characterization of the myostatin gene in the dromedary (Camelus dromedarius) at the genomic, transcript and protein level. The gene was found to share high structural and sequence similarity with other mammals, notably Old World camelids. 3D modeling highlighted several non-conservative SNP variants compared to the bovine, as well as putative functional variants involved in the stability of the myostatin dimer. NGS data for nine dromedaries from various countries revealed 66 novel SNPs, all of them falling either upstream or downstream the coding region. The analysis also confirmed the presence of three previously described SNPs in intron 1, predicted here to alter both splicing and transcription factor binding sites (TFBS), thus possibly impacting myostatin processing and/or regulation. Several putative TFBS were identified in the myostatin upstream region, some of them belonging to the myogenic regulatory factor family. Patterns of SNP distribution across countries, as suggested by Bayesian clustering of the nine dromedaries using the 69 SNPs, pointed to weak geographic differentiation, in line with known recurrent gene flow at ancient trading centers along caravan routes. Myostatin expression was investigated in a set of 8 skeletal muscles, both at transcript and protein level, via Digital Droplet PCR and Western Blotting, respectively. No significant differences were observed at the transcript level, while, at the protein level, the only significant differences concerned the promyostatin dimer (75 kDa), in four pair-wise comparisons, all involving the tensor fasciae latae muscle. Beside the mentioned band at 75 kDa, additional bands were observed at around 40 and 25 kDa, corresponding to the promyostatin monomer and the active C-terminal myostatin dimer, respectively. Their weaker intensity suggests that the unprocessed myostatin dimers could act as important reservoirs of slowly available myostatin forms. Under this assumption, the sequential cleavage steps may contribute additional layers of control within an already complex regulatory framework.

The Adaptive Evolution Database (TAED): A New Release of a Database of Phylogenetically Indexed Gene Families from Chordates

With the large collections of gene and genome sequences, there is a need to generate curated comparative genomic databases that enable interpretation of results in an evolutionary context. Such resources can facilitate an understanding of the co-evolution of genes in the context of a genome mapped onto a phylogeny, of a protein structure, and of interactions within a pathway. A phylogenetically indexed gene family database, the adaptive evolution database (TAED), is presented that organizes gene families and their evolutionary histories in a species tree context. Gene families include alignments, phylogenetic trees, lineage-specific dN/dS ratios, reconciliation with the species tree to enable both the mapping and the identification of duplication events, mapping of gene families onto pathways, and mapping of amino acid substitutions onto protein structures. In addition to organization of the data, new phylogenetic visualization tools have been developed to aid in interpreting the data that are also available, including TreeThrasher and TAED Tree Viewer. A new resource of gene families organized by species and taxonomic lineage promises to be a valuable comparative genomics database for molecular biologists, evolutionary biologists, and ecologists. The new visualization tools and database framework will be of interest to both evolutionary biologists and bioinformaticians.

The salmonid myostatin gene family: a novel model for investigating mechanisms that influence duplicate gene fate

BACKGROUND

Most fishes possess two paralogs for myostatin, a muscle growth inhibitor, while salmonids are presumed to have four: mstn1a, mstn1b, mstn2a and mstn2b, a pseudogene. The mechanisms responsible for preserving these duplicates as well as the depth of mstn2b nonfunctionalization within the family remain unknown. We therefore characterized several genomic clones in order to better define species and gene phylogenies.

RESULTS

Gene organization and sequence conservation was particularly evident among paralog groupings and within salmonid subfamilies. All mstn2b sequences included in-frame stop codons, confirming its nonfunctionalization across taxa, although the indels and polymorphisms responsible often differed. For example, the specific indels within the Onchorhynchus tshawytscha and O. nerka genes were remarkably similar and differed equally from other mstn2b orthologs. A phylogenetic analysis weakly established a mstn2b clade including only these species, which coupled with a shared 51 base pair deletion might suggest a history involving hybridization or a shared phylogenetic history. Furthermore, mstn2 introns all lacked conserved splice site motifs, suggesting that the tissue-specific processing of mstn2a transcripts, but not those of mstn2b, is due to alternative cis regulation and is likely a common feature in salmonids. It also suggests that limited transcript processing may have contributed to mstn2b nonfunctionalization.

CONCLUSIONS

Previous studies revealed divergence within gene promoters while the current studies provide evidence for relaxed or positive selection in some coding sequence lineages. These results together suggest that the salmonid myostatin gene family is a novel resource for investigating mechanisms that regulate duplicate gene fate as paralog specific differences in gene expression, transcript processing and protein structure are all suggestive of active divergence.

Genes contributing to genetic variation of muscling in sheep

Selective breeding programs aiming to increase the productivity and profitability of the sheep meat industry use elite, progeny tested sires. The broad genetic traits of primary interest in the progeny of these sires include skeletal muscle yield, fat content, eating quality, and reproductive efficiency. Natural mutations in sheep that enhance muscling have been identified, while a number of genome scans have identified and confirmed quantitative trait loci (QTL) for skeletal muscle traits. The detailed phenotypic characteristics of sheep carrying these mutations or QTL affecting skeletal muscle show a number of common biological themes, particularly changes in developmental growth trajectories, alterations of whole animal morphology, and a shift toward fast twitch glycolytic fibers. The genetic, developmental, and biochemical mechanisms underpinning the actions of some of these genetic variants are described. This review critically assesses this research area, identifies gaps in knowledge, and highlights mechanistic linkages between genetic polymorphisms and skeletal muscle phenotypic changes. This knowledge may aid the discovery of new causal genetic variants and in some cases lead to the development of biochemical and immunological strategies aimed at enhancing skeletal muscle.

Age-related changes in craniofacial morphology in GDF-8 (myostatin)-deficient mice

It is well recognized that masticatory muscle function helps determine morphology, although the extent of function on final form is still debated. GDF-8 (myostatin), a transcription factor is a negative regulator of skeletal muscle growth. A recent study has shown that mice homozygous for the myostatin mutation had increased muscle mass and craniofacial dysmorphology in adulthood. However, it is unclear whether such dysmorphology is present at birth. This study examines the onset and relationship between hypermuscularity and craniofacial morphology in neonatal and adult mice with GDF-8 deficiency. Fifteen (8 wild-type and 7 GDF-8 -/-), 1-day-old and 16 (9 wt and 7 GDF-8 -/-), 180-day-old male CD-1 mice were used. Standardized radiographs were taken of each head, scanned, traced, and cephalometric landmarks identified. Significant mean differences were assessed using a group x age, two-way ANOVA. Myostatin-deficient mice had significantly (P < 0.01) smaller body and masseter muscle weights and craniofacial skeletons at 1 day of age and significantly greater body and masseter muscle weights at 180 days of age compared to controls. Myostatin-deficient mice showed significantly (P < 0.001) longer and "rocker-shaped" mandibles and shorter and wider crania compared to controls at 180 days. Significant correlations were noted between masseter muscle weight and all cephalometric measurements in 180-day-old Myostatin-deficient mice. Results suggest that in this mouse model, there may be both early systemic skeletal growth deficiencies and later compensatory changes from hypermuscularity. These findings reiterate the role that masticatory muscle function plays on the ontogeny of the cranial vault, base, and most notably the mandible.

Assessment of selection mapping near the myostatin gene (GDF-8) in cattle

Domestic species provide a unique opportunity to examine the effects of selection on the genome. The myostatin gene (GDF-8) has been under strong selection in a number of cattle breeds because of its influence on muscle conformation and association with the 'double-muscling' phenotype. This study examined genetic diversity near this gene in a set of breeds including some nearly fixed for the allele associated with double-muscling (MH), some where the allele is segregating at intermediate frequency and some where the allele is absent. A set of microsatellites and SNPs were used to examine patterns of diversity at the centromeric end of bovine chromosome 2, the region where GDF-8 is located, using various statistical methods. The putative position of a selected gene was moved across the genomic region to determine, by regression, a best position of reduced heterozygosity. Additional analyses examined extended homozygous regions and linkage disequilibrium patterns. While the SNP data was not found to be very informative for selection mapping in this dataset, analyses of the microsatellite data provided evidence of selection on GDF-8 in several breeds. These results suggested that, of the breeds examined, the allele was most recently introduced into the South Devon. Limitations to the selection-mapping approach were highlighted from the analysis of the SNP data and the situation where the MH allele was at intermediate frequency.

Clinical, agricultural, and evolutionary biology of myostatin: a comparative review

The discovery of myostatin and our introduction to the "Mighty Mouse" over a decade ago spurred both basic and applied research and impacted popular culture as well. The myostatin-null genotype produces "double muscling" in mice and livestock and was recently described in a child. The field's rapid growth is by no means surprising considering the potential benefits of enhancing muscle growth in clinical and agricultural settings. Indeed, several recent studies suggest that blocking myostatin's inhibitory effects could improve the clinical treatment of several muscle growth disorders, whereas comparative studies suggest that these actions are at least partly conserved. Thus, neutralizing myostatin's effects could also have agricultural significance. Extrapolating between studies that use different vertebrate models, particularly fish and mammals, is somewhat confusing because whole genome duplication events have resulted in the production and retention of up to four unique myostatin genes in some fish species. Such comparisons, however, suggest that myostatin's actions may not be limited to skeletal muscle per se, but may additionally influence other tissues including cardiac muscle, adipocytes, and the brain. Thus, therapeutic intervention in the clinic or on the farm must consider the potential of alternative side effects that could impact these or other tissues. In addition, the presence of multiple and actively diversifying myostatin genes in most fish species provides a unique opportunity to study adaptive molecular evolution. It may also provide insight into myostatin's nonmuscle actions as results from these and other comparative studies gain visibility in biomedical fields.

Complex microsatellite dynamics in the myostatin gene within ruminants

A microsatellite has previously been identified in myostatin in cattle. Sequencing of this region from other artiodactyls coupled with phylogenetic analysis has been used to uncover the potential origins of the microsatellite event, which appears either to have been born twice or to have been gained and lost within ruminants. While caprids and ovids share the ancestral state with pigs and other mammals, microsatellite activity (length polymorphism) is uncovered in both deer and bovids. The dynamic process of microsatellite evolution, including birth, is discussed here in light of several models. Finally, these models are evaluated in the context of patterns of microsatellite conservation between closely related mammalian genomes.

Análise filogenética do gene da miogenina

Estudou-se a filogenia do gene da miogenina, um membro da família MyoD, reguladora da miogênese, que ocorre durante o desenvolvimento embrionário, e sua história evolutiva em espécies domésticas que apresentem seqüências de DNA depositadas no Genbank, comparando-se o índice de substituição de nucleotídeos não-sinônimos pelo índice de substituição sinônima. Valores maiores do que um (1) indicaram que o gene sofreu mudanças que tornaram o organismo mais adaptado ao ambiente. As árvores filogenéticas foram obtidas por máxima verossimilhança, e os índices de substituição sinônima e não-sinônima foram analisadas pelo método de parcimônia. Os resultados indicaram que, provavelmente, o gene sofreu evolução adaptativa no grupo Ruminantia, Bos taurus e Ovis aries, depois que essas espécies divergiram do ancestral comum. Para as outras espécies analisadas, o gene parece ter evoluído de modo conservativo.

Functional analysis of the Myostatin gene promoter in sheep

Compared with the understanding for the functional mechanism of the myostatin gene, little is known about the regulatory mechanism of the myostatin gene transcription and expression. To better understand the function of the myostatin gene promoter (MSTNpro) in the transcriptional regulation of the myostatin gene and to further investigate the transcriptional regulation mechanism of the myostatin gene, the promoter region of the myostatin gene in sheep has been cloned in our recent study (AY918121). In this study, the wild (W) type MSTNPro(W)-EGFP vectors and E-box (E) (CANNTG) mutant (M) type MSTNPro(E(3+5+7)M)-EGFP vectors were constructed and the transcriptional regulation activities were compared by detecting the fluorescent strength of EGFP (enhanced green fluorescent protein) in C2C12 myoblasts (or myotubes) and sheep fibroblasts transfected with the vectors. Results showed that the 0.3-1.2 kb sheep myostatin promoter could activate the transcription and expression of EGFP gene in C2C12 myoblasts to different extent and the 1.2 kb promoter was the strongest. However, fluorescence was not observed in the sheep fibroblasts transfected with the 1.2 kb sheep myostatin promoter. These results suggested that the specific nature of the myostatin gene expression in skeletal muscle was attributed to the specific nature of the myostatin promoter activity. The increasing growth density of C2C12 myoblasts inhibited the transcriptional regulation activity of the wild type sheep myostatin promoter by a mechanism of feedback. The transcriptional regulation activity of the 1.2 kb wild type sheep myostatin promoter increased significantly after C2C12 myoblasts were differentiated, while the activity of 1.2 kb E(3+5+7)-mutant type myostatin promoter had no obvious change. This result suggested that MyoD may be responsible for the difference of the myostatin gene transcription and expression between growing and differentiating conditions by binding to E-box of the myostatin promoter.

The quest for natural selection in the age of comparative genomics

Continued genome sequencing has fueled progress in statistical methods for understanding the action of natural selection at the molecular level. This article reviews various statistical techniques (and their applicability) for detecting adaptation events and the functional divergence of proteins. As large-scale automated studies become more frequent, they provide a useful resource for generating biological null hypotheses for further experimental and statistical testing. Furthermore, they shed light on typical patterns of lineage-specific evolution of organisms, on the functional and structural evolution of protein families and on the interplay between the two. More complex models are being developed to better reflect the underlying biological and chemical processes and to complement simpler statistical models. Linking molecular processes to their statistical signatures in genomes can be demanding, and the proper application of statistical models is discussed.

Myostatin (MSTN) gene duplications in Atlantic salmon (Salmo salar): evidence for different selective pressure on teleost MSTN-1 and -2

Whereas the negative muscle regulator myostatin (MSTN) in mammals is almost exclusively expressed in the muscle by a single encoding gene, teleost fish possess at least two MSTN genes which are differentially expressed in both muscular and non-muscular tissues. Duplicated MSTN-1 genes have previously been identified in the tetraploid salmonid genome. From Atlantic salmon we succeeded in isolating the paralogous genes of MSTN-2, which shared about 70% identity with MSTN-1a and -1b. The salmon MSTN-2a cDNA encoded a predicted protein of 363 residues and included the conserved C-terminal bioactive domain. MSTN-2a seemed to be primarily expressed in the brain, and a functional role of teleost MSTN-2 in the neurogenesis similar to the inhibitory action of the closely related GDF-11 in the mammalian brain was proposed. In contrast, a frame-shift mutation in exon 1 of salmon MSTN-2b would lead to the synthesis of a putatively non-functional truncated protein. The absence of processed MSTN-2b mRNA in the examined tissues indicated that this gene has become a non-functional pseudogene. The differential, but partially overlapping, expression patterns of salmon MSTN-2a, -1a and -1b in muscular and non-muscular tissues are probably due to the different arrangement of the potential cis-acting regulatory elements identified in their putative promoter regions. Single and paired E-boxes in the MSTN-1b promoter were shown to bind both homo-and hetero-dimers of the myogenic regulatory factor MyoD and E47 in vitro of importance for initiating the myogenic program. Analyses of nucleotide substitution patterns indicated that the teleost MSTNs essentially have evolved under purifying selection, but a subset of amino acid sites under positive selective pressure were identified within the MSTN1 branch. The results may reflect the evolutionary forces related to adoption of the different functional roles proposed for the teleost MSTN isoforms. The phylogenetic analysis of multiple vertebrate MSTNs suggested at least two separate gene duplication events in the fish lineage. Linkage analysis of polymorphic microsatellites within intron 2 of salmon MSTN-1a and -1b mapped the two genes to different linkage groups in agreement with the tetraploid origin of the duplicated salmonid MSTN-1 and MSTN-2 genes.

Characterization of rainbow trout myostatin-2 genes (rtMSTN-2a and -2b): genomic organization, differential expression, and pseudogenization

Myostatin is an extremely potent negative regulator of vertebrate skeletal muscle development. A phylogenetic analysis suggests that salmonids should possess four distinct genes, although only MSTN-1 orthologs have been characterized. Described herein are the rainbow trout (rt) MSTN-2a and -2b genes and subsequence analysis of their promoters and their quantitative expression profiles. Both genes are similarly organized, contain several putative myogenic response elements, and are legitimate MSTN-2 orthologs based on Bayesian analyses. However, rtMSTN-2b contains two in-frame stop codons within the first exon and unspliced variants of both transcripts were expressed in a tissue-specific manner. Complete splicing of rtMSTN-2a occurred only in brain, where expression is highest, whereas rtMSTN-2b transcripts were mostly present in unspliced forms. The presence of stop codons in the rtMSTN-2b open reading frame and the expression of mostly unspliced transcripts indicate that this particular homolog is a pseudogene. These results confirm our previous phylogenetic analysis and suggest that all salmonids likely possess four distinct myostatin genes. The tissue-specific expression and differential processing of both rtMSTN-2 transcripts as well the pseudogenization of rtMSTN-2b may reflect compensatory and adaptive responses to tetraploidization and may help limit rtMSTN-2a's influences primarily to neural tissue.

Two novel missense mutations in the myostatin gene identified in Japanese patients with Duchenne muscular dystrophy

Background

Myostatin is a negative regulator of skeletal muscle growth. Truncating mutations in the myostatin gene have been reported to result in gross muscle hypertrophy. Duchenne muscular dystrophy (DMD), the most common lethal muscle wasting disease, is a result of an absence of muscle dystrophin. Although this disorder causes a rather uniform pattern of muscle wasting, afflicted patients display phenotypic variability. We hypothesized that genetic variation in myostatin is a modifier of the DMD phenotype.

Methods

We analyzed 102 Japanese DMD patients for mutations in the myostatin gene.

Results

Two polymorphisms that are commonly observed in Western countries, p.55A>T and p.153K>R, were not observed in these Japanese patients. An uncommon polymorphism of p.164E>K was uncovered in four cases; each patient was found to be heterozygous for this polymorphism, which had the highest frequency of the polymorphism observed in the Japanese patients. Remarkably, two patients were found to be heterozygous for one of two novel missense mutations (p.95D>H and p.156L>I). One DMD patient carrying a novel missense mutation of p.95D>H was not phenotypically different from the non-carriers. The other DMD patient was found to carry both a novel mutation (p.156L>I) and a known polymorphism (p.164E>K) in one allele, although his phenotype was not significantly modified. Any nucleotide change creating a target site for micro RNAs was not disclosed in the 3' untranslated region.

Conclusion

Our results indicate that heterozygous missense mutations including two novel mutations did not produce an apparent increase in muscle strength in Japanese DMD cases, even in a patient carrying two missense mutations.

Human adaptive evolution at Myostatin (GDF8), a regulator of muscle growth

Myostatin (GDF8) is a negative regulator of muscle growth in mammals, and loss-of-function mutations are associated with increased skeletal-muscle mass in mice, cattle, and humans. Here, we show that positive natural selection has acted on human nucleotide variation at GDF8, since the observed ratio of nonsynonymous:synonymous changes among humans is significantly greater than expected under the neutral model and is strikingly different from patterns observed across mammalian orders. Furthermore, extended haplotypes around GDF8 suggest that two amino acid variants have been subject to recent positive selection. Both mutations are rare among non-Africans yet are at frequencies of up to 31% in sub-Saharan Africans. These signatures of selection at the molecular level suggest that human variation at GDF8 is associated with functional differences.

Sequence diversity and rates of molecular evolution between sheep and cattle genes

Experiments that aim to identify genes of importance in sheep are currently inhibited by a paucity of genomic resources. One approach, therefore, is to exploit the wealth of data and associated capabilities becoming available for the bovine genome. Cross-species application of microarrays and comparative sequencing to identify single nucleotide polymorphisms are two possibilities; however, both are dependant on the level of nucleotide sequence similarity between the two species. This study used 120 gene orthologues consisting of over 60 kb of aligned sequence to estimate the gene diversity between cattle and sheep. Less than 3% of protein-coding nucleotide positions were found to be different, indicating that the prospect for successfully using cross-species strategies is high. Substitution at synonymous sites ranged between 6.9 and 7.7% (+/- 0.3%), and was higher than at non-synonymous sites (1.4-1.7 +/- 0.1%). The relative rate test was used to determine whether the observed mutation rates were constant between the two lineages. While the rate at synonymous sites appeared constant, the rate at non-synonymous sites was significantly higher within the caprinae lineage (sheep) when compared with bovinae (cattle; chi2 = 10.03; d.f. = 1, P < 0.01). This is the first demonstration that variable rates of molecular evolution may be present within the family Bovidae.

Microsatellite polymorphism in intron 1 of the bovine myostatin gene

Microsatellites within genes have become important because of their association with genetic diseases in humans. A novel microsatellite was identified in the first intron of the bovine myostatin gene. It is characterized by a mononucleotide core motif that exhibits polymorphic sequence variants (from 12 to 21 repeats) within and between some bovine breeds. Structural analysis of the microsatellite region in bovines as well as in closely related species permitted to trace the possible mechanisms for its structural evolution across the order Artiodactyla.

Myostatinand its implications on animal breeding: a review

Myostatin, or growth and differentiation factor 8 (GDF8), has been identified as the factor causing a phenotype known as double muscling, in which a series of mutations render the gene inactive, and therefore, unable to regulate muscle fibre deposition. This phenotype occurs at a high frequency in some breeds of cattle such as Belgian Blue and Peidmontese. Phylogenetic analysis has shown that there has been positive selection pressure for non-synonymous mutations within the myostatin gene family, around the time of the divergence of cattle, sheep and goats, and these positive selective pressures on non-ancestral myostatin are relatively recent. To date, there have been reports of nine mutations in coding regions of myostatin that cause non-synonymous changes, of which three cause missense mutations, including two in exon 1 and one in exon 2. The remaining six mutations, located in exons 2 and 3, result in premature stop codons, which are the mutations responsible for the double-muscling phenotype. Unfortunately, breed management problems exist for double-muscled cattle, such as birthing difficulties, which can be overcome through genetically controlled breeding programmes, as shown in this review.