Evolutionary implications of three novel members of the human sarcomeric myosin heavy chain gene family

Sarcomeric myosin heavy chain (MyHC) is the major contractile protein of striated muscle. Six tandemly linked skeletal MyHC genes on chromosome 17 and two cardiac MyHC genes on chromosome 14 have been previously described in the human genome. We report the identification of three novel human sarcomeric MyHC genes on chromosomes 3, 7, and 20, which are notable for their atypical size and intron-exon structure. Two of the encoded proteins are structurally most like the slow-beta MyHC, whereas the third one is closest to the adult fast IIb isoform. Data from pairwise comparisons of aligned coding sequences imply the existence of ancestral genomes with four sarcomeric genes before the emergence of a dedicated smooth muscle MyHC gene. To further address the evolutionary relationships of the distinct sarcomeric and nonsarcomeric rod sequences, we have identified and further annotated human genomic DNA sequences corresponding to 14 class-II MyHCs. An extensive analysis provides a timeline for intron gain and loss, gene contraction and expansion, and gene conversion among genes encoding class-II myosins. One of the novel human genes is found to have introns at positions shared only with the molluscan catchin/MyHC gene, providing evidence for the structure of a pre-Cambrian ancestral gene.     

Nonmuscle myosin IIB, a sarcomeric component in the extraocular muscles

Extraocular muscles (EOMs) are categorized as skeletal muscles; however, emerging evidence indicates that their gene expression profile, metabolic characteristics and functional properties are significantly different from the prototypical members of this muscle class. Gene expression profiling of developing and adult EOM suggest that many myofilament and cytoskeletal proteins have unique expression patterns in EOMs, including the maintained expression of embryonic and fetal isoforms of myosin heavy chains (MyHC), the presence of a unique EOM specific MyHC and mixtures of both cardiac and skeletal muscle isoforms of thick and thin filament accessory proteins. We demonstrate that nonmuscle myosin IIB (nmMyH IIB) is a sarcomeric component in ∼ 20% of the global layer fibers in adult rat EOMs. Comparisons of the myofibrillar distribution of nmMyHC IIB with sarcomeric MyHCs indicate that nmMyH IIB co-exists with slow MyHC isoforms. In longitudinal sections of adult rat EOM, nmMyHC IIB appears to be restricted to the A-bands. Although nmMyHC IIB has been previously identified as a component of skeletal and cardiac sarcomeres at the level of the Z-line, the novel distribution of this protein within the A band in EOMs is further evidence of both the EOMs complexity and unconventional phenotype.     

The sarcomeric myosin heavy chain gene family in the dog: analysis of isoform diversity and comparison with other mammalian species

Sarcomeric myosin heavy chains (MyHC) are the major contractile proteins of cardiac and skeletal muscles and belong to class II MyHC. In this study the sequences of nine sarcomeric MyHC isoforms were obtained by combining assembled contigs of the dog genome draft available in the NCBI database. With this information available the dog becomes the second species, after human, for which the sequences of all members of the sarcomeric MyHC gene family are identified. The newly determined sequences of canine MyHC isoforms were aligned with their orthologs in mammals, forming a set of 38 isoforms, to search for the molecular features that determine the structural and functional specificity of each type of isoform. In this way the structural motifs that allow identification of each isoform and are likely determinants of functional properties were identified in six specific regions (surface loop 1, loop 2, loop 3, converter, MLC binding region, and S2 proximal segment).     

Digoxigenin-labeled RNA probes for untranslated regions enable the isoform-specific gene expression analysis of myosin heavy chains in whole-mount in situ hybridization

Myosin heavy chains (MyHCs), which are encoded by myosin heavy chain (Myh) genes, are the most abundant proteins in myofiber. Among the 11 sarcomeric Myh isoform genes in the mammalian genome, seven are mainly expressed in skeletal muscle. Myh genes/MyHC proteins share a common role as force producing units with highly conserved sequences, but have distinct spatio-temporal expression patterns. As such, the expression patterns of Myh genes/MyHC proteins are considered as molecular signatures of specific fiber types or the regenerative status of mammalian skeletal muscles. Immunohistochemistry is widely used for identifying MyHC expression patterns; however, this method is costly and is not ideal for whole-mount samples, such as embryos. In situ hybridization (ISH) is another versatile method for the analysis of gene expression, but is not commonly applied for Myh genes, partly because of the highly homologous sequences of Myh genes. Here we demonstrate that an ISH analysis with the untranslated region (UTR) sequence of Myh genes is cost-effective and specific method for analyzing the Myh gene expression in whole-mount samples. Digoxigenin (DIG)-labeled antisense probes for UTR sequences, but not for protein coding sequences, specifically detected the expression patterns of respective Myh isoform genes in both embryo and adult skeletal muscle tissues. UTR probes also revealed the isoform gene-specific polarized localization of Myh mRNAs in embryonic myofibers, which implied a novel mRNA distribution mechanism. Our data suggested that the DIG-labeled UTR probe is a cost-effective and versatile method to specifically detect skeletal muscle Myh genes in a whole-mount analysis.     

The genome of the diploid anuran <Emphasis Type="Italic">Xenopus tropicalis</Emphasis> contains a novel array of sarcoplasmic myosin heavy chain genes expressed in larval muscle and larynx

The sarcomeric myosin heavy chain (MyHC) proteins are a family of molecular motors responsible for the transduction of chemical energy into mechanical work in striated muscle. The vertebrate genome contains multiple copies of the MyHC gene, and expression of different isoforms correlates with differences in the physiological properties of muscle fibers. Most MyHC isoforms are found in two arrays, one containing the "fast-twitch" skeletal muscle isoforms and the other the "slow-twitch" or cardiac isoforms. To extend our understanding of MyHC evolution, we have examined the genome of the anuran Xenopus tropicalis. The X. tropicalis genome includes15 full-length MyHC genes organized in seven genomic locations. One unique array of MyHC genes is similar to the mammalian fast-skeletal array, but is not found in amniotes. The isoforms in this array are expressed during larval stages and in muscles of the adult larynx. Duplication of the fast-skeletal MyHC array appears to have led to expression divergence of muscle proteins in the larval and adult stages of the anuran life cycle. A striking similarity of gene order between regions flanking X. tropicalis MyHC arrays and human arrays was evident; genomic organization of MyHC isoforms may thus be highly conserved across tetrapods.     

Absence of the functional Myosin heavy chain 2b isoform in equine skeletal muscles

Nucleotide sequences which included the full coding region for three types of myosin heavy chain (MyHC) isoforms were determined from equine skeletal muscles. The deduced amino acid sequences were 1937, 1938, and 1935 residues for the MyHC-2a, -2x, and -slow, respectively. No MyHC-2b isoform was amplified from the equine muscle cDNA except for one pseudogene fragment. One nucleotide was inserted in the coding region of the equine pseudogene product, a minute amount of which was expressed in the skeletal muscle. The 596 bp sequence of the equine MyHC pseudogene was categorized into the MyHC-2b genes on the phylogenetic tree of the mammalian MyHC genes. These results suggest that an ancestral MyHC-2b gene had lost its function and changed to a pseudogene during the course of horse history. The MyHC genes in some ungulates were analyzed through the PCR amplifications using the MyHC isoform-specific primers to confirm the presence of the MyHC-2b and -2x genes. The exon coding the 3' untranslated region of the MyHC-2x was successfully amplified from the all ungulates examined; however, that of the MyHC-2b gene was amplified only from horses, pigs and lesser mouse deer. The PCR analyses from rhinoceros, sika deer, moose, giraffes, water buffalo, bovine, Japanese serow and sheep genes implied the absence of the MyHC-2b-specific sequence in their genomes. These results suggest that the MyHC-2b gene independently lost its function in some ungulate species.     

A distinct subclass of mammalian striated myosins: structure and molecular evolution of "superfast" or masticatory myosin heavy chain

"Superfast" or masticatory myosin is the molecular motor in the powerful and specialized jaw-closing muscles of carnivores, folivores, and frugivores. This myosin presumably underpins the unusual high force and moderate shortening velocity of muscle fibers expressing it. Here, we report the cloning and sequencing of the cDNA encoding the full-length masticatory myosin heavy chain (MyHC) from cat temporalis muscle. This was obtained by immunoscreening a cDNA expression library and RACE-PCR (rapid amplification of cDNA ends–PCR). Sequence comparisons at the DNA and amino acid levels show that masticatory MyHC has less than 70% homology to known striated MyHCs, compared with 87–96% between other mammalian fast isoforms themselves. Nucleotide substitution rates at the nonsynonymous sites between masticatory MyHC and other mammalian striated MyHCs are considerably higher than between these striated MyHCs themselves. Phylogenetic analysis revealed that masticatory MyHC diverged from invertebrate MyHC before the avian cardiac MyHC subclass and the mammalian fast/developmental and slow/cardiac MyHC subclasses. Masticatory MyHC is thus a distinct new subclass of vertebrate striated myosins. The early divergence from invertebrate MyHC, combined with immunochemical evidence of its expression in reptilian and shark jaw-closing muscles, suggests that masticatory MyHC evolved in early gnathostomes, driven by benefits derived from powerful jaw closure. During the mammalian radiation, some taxa continued to express it, while others adapted to new types of food and eating habits by replacing masticatory MyHC with more appropriate isoforms normally found in limb and cardiac muscles.     

Postnatal myosin heavy chain isoform expression in normal mice and mice null for IIb or IId myosin heavy chains

The patterns of myosin heavy chain (MyHC) isoform expression in the embryo and in the adult mouse are reasonably well characterized and quite distinct. However, little is known about the transition between these two states, which involves major decreases and increases in the expression of several MyHC genes. In the present study, the expression of seven sarcomeric MyHCs was analyzed in the hindlimb muscles of wild-type mice and in mice null for the MyHC IIb or IId/x genes at several time points from 1 day of postnatal life (dpn) to 20 dpn. In early postnatal life, the developmental isoforms (embryonic and perinatal) comprise >90% of the total MyHC expression, while three adult fast isoforms (IIa, IIb, and IId) comprise <1% of the total MyHC protein. However, between 5 and 20 dpn their expression increases to comprise >90% of the total MyHC. Expression of each of the three adult fast isoforms occurs in a spatially and temporally distinct manner. We also show that alpha MyHC, which is almost exclusively expressed in the heart, is expressed in scattered fibers in all hindlimb muscles during postnatal development. Surprisingly, the timing and localization of expression of the MyHC isoforms is unchanged in IIb and IId/x null mice, although the magnitude of expression is altered for some isoforms. Together these data provide a comprehensive overview of the postnatal expression pattern of the sarcomeric MyHC isoforms in the mouse hindlimb.     

The evolutionary relationship of avian and mammalian myosin heavy-chain genes

Summary Sequence comparisons of avian and mammalian skeletal and cardiac myosin heavy-chain isoforms are used to examine the evolutionary relationships of sarcomeric myosin multigene families. Mammalian fast-myosin heavy-chain isoforms forms from different species, with comparable developmental expression, are more similar to each other than they are to other fast isoforms within the same genome. In contrast, the developmentally regulated chicken fast isoforms are more similar to each other than they are to myosin heavy-chain isoforms in other species. Extensive regions of nucleotide identity among the chicken fast myosin heavy chains and in the mouse and rat α- and β-cardiac myosin heavy-chain sequences suggest that geneconversion-like mechanisms have played a major role in the concerted evolution of these gene families. We also conclude that the chicken fast myosin heavy-chain multigene family has undergone recent expansion subsequent to the divergence of birds and mammals and that both the developmental regulation and the specialization of myosin isoforms have likely developed independently in birds and mammals.     

Interspecific and Intragenic Differences in Codon Usage Bias Among Vertebrate Myosin Heavy-Chain Genes

Synonymous codon usage bias is a broadly observed phenomenon in bacteria, plants, and invertebrates and may result from selection. However, the role of selective pressures in shaping codon bias is still controversial in vertebrates, particularly for mammals. The myosin heavy-chain (MyHC) gene family comprises multiple isoforms of the major force-producing contractile protein in cardiac and skeletal muscles. Slow and fast genes are tandemly arrayed on separate chromosomes, and have distinct patterns of functionality and expression in muscle. We analyze both full-length MyHC genes (~5400?bp) and a larger collection of partial sequences at the 3' end (~500?bp). The MyHC isoforms are an interesting system in which to study codon usage bias because of their length, expression, and critical importance to organismal mobility. Codon bias and GC content differs among MyHC genes with regards to functional type, isoform, and position within the gene. Codon bias even varies by isoform within a species. We find evidence in favor of both chromosomal influences on nucleotide composition and selection against nonsense errors (SANE) acting on codon usage in MyHC genes. Intragenic variation in codon bias and elongation rate is significant, with a strong trend for increasing codon bias and elongation rate towards the 3' end of the gene, although the trend is dependent upon the degeneracy class of the codons. Therefore, patterns of codon usage in MyHC genes are consistent with models supporting SANE as a major force shaping codon usage.     

Lamprey contractile protein genes mark different populations of skeletal muscles during development

Agnathan lampreys retain ancestral characteristics of vertebrates in the morphology of skeletal muscles derived from two mesodermal regions: trunk myotomes and unsegmented head mesoderm. During lamprey development, some populations of myoblasts migrate via pathways that differ from those of gnathostomes. To investigate the evolution of skeletal muscle differentiation in vertebrates, we characterize multiple contractile protein genes expressed in the muscle cells of the Japanese lamprey, Lethenteron japonicum. Lamprey actin gene LjMA2, and myosin heavy chain (MyHC) genes LjMyHC1 and LjMyHC2 are all expressed in the developing skeletal muscle cells of early embryos. However, LjMyHC1 and LjMyHC2 are expressed only in cells originating from myotomes, while LjMA2 is expressed in both myotomal and head musculature. Thus, in lampreys, myotomes and head mesoderm differ in the use of genes encoding contractile protein isoforms. Phylogenetic tree analyses including lamprey MyHCs suggest that the variety of muscle MyHC isoforms in different skeletal muscles may correspond to the morphological complexity of skeletal muscles of different vertebrate species. Another lamprey actin gene LjMA1 is likely to be the first smooth muscle actin gene isolated from non-tetrapods. We conclude that, in vertebrate evolution, the different regulatory systems for striated and smooth muscle-specific genes may have been established before the agnathan/gnathostome divergence.     

Effects of hypothyroidism on myosin heavy chain composition and fibre types of fast skeletal muscles in a small marsupial, Antechinus flavipes

Effects of drug-induced hypothyroidism on myosin heavy chain (MyHC) content and fibre types of fast skeletal muscles were studied in a small marsupial, Antechinus flavipes. SDS-PAGE of MyHCs from the tibialis anterior and gastrocnemius revealed four isoforms, 2B, 2X, 2A and slow, in that order of decreasing abundance. After 5 weeks treatment with methimazole, the functionally fastest 2B MyHC significantly decreased, while 2X, 2A and slow MyHCs increased. Immunohistochemistry using monospecific antibodies to each of the four MyHCs revealed decreased 2b and 2x fibres, and increased 2a and hybrid fibres co-expressing two or three MyHCs. In the normally homogeneously fast superficial regions of these muscles, evenly distributed slow-staining fibres appeared, resembling the distribution of slow primary myotubes in fast muscles during development. Hybrid fibres containing 2A and slow MyHCs were virtually absent. These results are more detailed but broadly similar to the earlier studies on eutherians. We hypothesize that hypothyroidism essentially reverses the effects of thyroid hormone on MyHC gene expression of muscle fibres during myogenesis, which differ according to the developmental origin of the fibre: it induces slow MyHC expression in 2b fibres derived from fast primary myotubes, and shifts fast MyHC expression in fibres of secondary origin towards 2A, but not slow, MyHC.     

Modulation of myosin filament organization by C-protein family members.

We have analyzed the interactions between two types of sarcomeric proteins: myosin heavy chain (MyHC) and members of an abundant thick filament-associated protein family (myosin-binding protein; MyBP). Previous work has demonstrated that when MyHC is transiently transfected into mammalian nonmuscle COS cells, the expressed protein forms spindle-shaped structures consisting of bundles of myosin thick filaments. Co-expression of MyHC and MyBP-C or -H modulates the MyHC structures, resulting in dramatically longer cables consisting of myosin and MyBP encircling the nucleus. Immunoelectron microscopy indicates that these cable structures are more uniform in diameter than the spindle structures consisting solely of MyHC, and that the myosin filaments are compacted in the presence of MyBP. Deletion analysis of MyBP-H indicates that cable formation is dependent on the carboxy terminal 24 amino acids. Neither the MyHC spindles nor the MyHC/MyBP cables associate with the endogenous actin cytoskeleton of the COS cell. While there is no apparent co-localization between these structures and the microtubule network, colchicine treatment of the cells promotes the formation of longer assemblages, suggesting that cytoskeletal architecture may physically impede or regulate polymer formation/extension. The data presented here contribute to a greater understanding of the interactions between the MyBP family and MyHC, and provide additional evidence for functional homology between MyBP-C and MyBP-H.     

Evolution of sarcomeric myosin heavy chain genes: evidence from fish

Myosin heavy chain (MYH) is a major structural protein, integral to the function of sarcomeric muscles. We investigated both exon-intron organization and amino acid sequence of sarcomeric MYH genes to infer their evolutionary history in vertebrates. Our results were consistent with the hypothesis that a multigene family encoded MYH proteins in the ancestral chordate, one gene ancestral to human MYH16 and its homologues and another ancestral to all other vertebrate sarcomeric MYH genes. We identified teleost homologues of mammalian skeletal and cardiac MYH genes, indicating that the ancestors of those genes were present before the divergence of actinopterygians and sarcopterygians. Indeed, the ancestral skeletal genes probably duplicated at least once before the divergence of teleosts and tetrapods. Fish homologues of mammalian skeletal MYH are expressed in skeletal tissue and homologues of mammalian cardiac genes are expressed in the heart but, unlike mammals, there is overlap between these expression domains. Our analyses inferred two other ancestral vertebrate MYH genes, giving rise to human MYH14 and MYH15 and their homologues. Relative to the skeletal and cardiac genes, MYH14 and MYH15 homologues are characterized by evolution of intron position, differences in evolutionary rate between the functionally differentiated head and rod of the myosin protein, and possible evolution of function among vertebrate classes. Tandem duplication and gene conversion appear to have played major roles in the evolution of at least cardiac and skeletal MYH genes in fish. One outcome of this high level of concerted evolution is that different fish taxa have different suites of MYH genes, i.e., true orthologs do not exist.     

The complete sequence of the human beta-myosin heavy chain gene and a comparative analysis of its product

We have isolated and sequenced the gene and the cDNA coding for the human cardiac beta-myosin heavy chain (designated MYH7). The gene is 22,883 bp long. The 1935 amino acids of this protein (Mr223,111) are encoded by 38 exons. The 5' untranslated region (86 bp) is split by two introns. The 3' untranslated region is 114 bp long. Three Alu repeats were identified within the gene and a fourth one in the 3' flanking intergenic region. The molecular organization of this gene reflects the conservative pattern with respect to size, coding ratio, and number or position of introns characteristic of vertebrate sarcomeric myosin heavy chain genes. The protein sequence of the human beta-heavy chain was compared with corresponding (homologous) sequences of rabbit, rat, and hamster as well as with the (heterologous) embryonic heavy chain sequences of rat, chicken, and man. The results show that protein subregions responsible for basic functions of myosin heavy chains (nucleotide binding and actin binding) are very similar in homologous and heterologous heavy chains. Regions that differ in their primary sequences in heterologous heavy chains appear to be highly conserved within mammalian beta-myosin heavy chains. Constant and variable subregions of heavy chains are discussed in terms of functional significance and evolutionary relatedness.     

Fiber types in rat laryngeal muscles and their transformations after denervation and reinnervation.

The intrinsic laryngeal muscles cricothyroid (CT) and thyroarythenoid (TA) differ in myosin expression. CT expresses limb myosin heavy chains (MyHCs) and TA expresses an MyHC found in extraocular (EO) muscles, in addition to limb isoforms. We used immunohistochemical (IHC) analyses with highly specific monoclonal antibodies (MAbs) against various MyHCs to study muscle fiber types in rat CT and TA and to investigate whether nerves to laryngeal muscles control MyHC expression. CT was found to have the full complement of limb fiber types. TA had three major fiber types: 2b/eo, co-expressing 2B and EO MyHCs, 2x/2b, co-expressing 2X and 2B MyHCs, and 2x, expressing 2X MyHC. Type 2a and slow fibers were absent. TA consisted of two divisions: the external division (TA-X), which is homogeneously 2b/eo, and the vocalis division (TA-V), composed principally of 2x and 2b/eo fibers with a minority of 2x/2b fibers. TA-V had two compartments that differ in fiber type composition. At 4 weeks after cutting and re-uniting the recurrent laryngeal nerve (RLN), many 2b/eo fibers in the TA-X began to express 2X MyHC, while EO and 2B MyHC expression in these fibers progressively declined. By 12 weeks, up to 16.5% of fibers in the TA-X were of type 2x. These findings suggest that nerve fibers originally innervating 2x fibers in TA-V and other muscles have randomly cross-innervated 2b/eo fibers in the TA-X and converted them into 2x fibers. We conclude that CT and TA are distinct muscle allotypes and that laryngeal muscle fibers are subject to neural regulation.     

Characterization of diverse forms of myosin heavy chain expressed in adult human skeletal muscle.

In an attempt to define myosin heavy chain (MHC) gene organization and expression in adult human skeletal muscle, we have isolated and characterized genomic sequences corresponding to different human sarcomeric MHC genes (1). In this report, we present the complete DNA sequence of two different adult human skeletal muscle MHC cDNA clones, one of which encodes the entire light meromyosin (LMM) segment of MHC and represents the longest described MHC cDNA sequence. Additionally, both clones provide new sequence data from a 228 amino acid segment of the MHC tail for which no protein or DNA sequence has been previously available. One clone encodes a "fast" form of skeletal muscle MHC while the other clone most closely resembles a MHC form described in rat cardiac ventricles. We show that the 3' untranslated region of skeletal MHC cDNAs are homologous from widely separated species as are cardiac MHC cDNAs. However, there is no homology between the 3' untranslated region of cardiac and skeletal muscle MHCs. Isotype-specific preservation of MHC 3' untranslated sequences during evolution suggests a functional role for these regions.     

Myosin heavy chain protein and gene expression in the masseter muscle of adult patients with distal or mesial malocclusion

The aim of this study was to determine the amount of myosin heavy chain (MyHC) proteins and MyHC mRNA in muscles of patients with different positions of the mandible. Ten adult patients for orthognathic surgery were divided into two groups: distal and mesial malocclusion. The mRNA expression of two MyHC isoforms of the anterior and posterior part of the right and left side of the human masseter muscle was analysed with a competitive RT-PCR assay. An exogenous template that includes oligonucleotide sequences specific for sarcomeric MyHC isoforms (1 and 2x) was constructed and utilized as competitor. Different isoforms of the MyHC protein were identified by Western blot analysis. In the total mRNA pool of the masseter muscle, the MyHC 1 mRNA level was 25.5 +/- 7.6% and the MyHC 2x mRNA was 2.5 +/- 1.2%. The anterior part of the masseter muscle from patients with distal occlusion contained more type 1 and 2x MyHC mRNA, as compared to patients with mesial occlusion (P < 0.05). No difference in the protein distribution was observed. The differences in mRNA expression may be caused by the enforced stress of the masticatory muscle in distal occlusion because of the disadvantageous pivot.