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.     

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).     

Characterization of sarcomeric myosin heavy chain genes

Myosin heavy chain is encoded by a large multigene family. Using pMHC-25, a recombinant cDNA clone isolated from the rat myogenic cell line L6E9, four members of this family in the rat have been isolated and shown to be tissue-specific and developmentally regulated. The coding regions of these genes share regions of homology interspaced with regions of non-homology. Detailed analysis of one embryonic and one adult myosin heavy chain gene shows that the coding sequences are interrupted by numerous intervening sequences whose number, size, and distribution do not appear to be conserved in the same organism or between species.     

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 human beta-myosin heavy chain gene: sequence diversity and functional characteristics of the protein

The beta-myosin heavy chain gene (MYH7) encodes the motor protein that drives myocardial contraction. It has been proven to be a disease gene for hypertrophic cardiomyopathy (HCM). We analyzed the DNA sequence variation of MYH7 (about 16 kb) of eight individuals: six patients with HCM and two healthy controls. The overall DNA sequence identity was up to 97.2% compared to Jaenicke and coworkers (Jaenicke et al. [1990] Genomics 8:194-206), while the corresponding amino acid sequences revealed 100% identity. In HCM patients, eleven nucleotide substitutions were identified but no causative disease mutation was found: six were detected in coding, four in intronic, and one in 5' regulatory regions. The average nucleotide diversity across this locus was 0.015% with an average of 0.02% in the coding and 0.012% in the noncoding sequence. Analysis of the kinetic behaviour of beta-MHC in the intact contractile structure of normal individuals and HCM patients revealed apparent rate constants of tension development ranging between 1.58 s(-1) and 1.48 s(-1).     

Expression and DNA sequence analysis of a human embryonic skeletal muscle myosin heavy chain gene.

Vertebrate myosin heavy chains (MHC) are represented by multiple genes that are expressed in a spatially and temporally distinct pattern during development. In order to obtain molecular probes for developmentally regulated human MHC isoforms, we used monoclonal antibodies to screen an expression cDNA library constructed from primary human myotube cultures. A 3.4 kb cDNA was isolated that encodes one of the first MHCs to be transcribed in human skeletal muscle development. A portion of the corresponding gene encoding this isoform has also been isolated. Expression of this embryonic MHC is a hallmark of muscle regeneration after birth and is a characteristic marker of human muscular dystrophies. During normal human development, expression is restricted to the embryonic period of development prior to birth. In primary human muscle cell cultures, devoid of other cell types, mRNA accumulation begins as myotubes form, reaches a peak 2 days later and declines to undetectable levels within 10 days. The expression of the protein encoded by the embryonic skeletal MHC gene follows a similar time course, lagging behind the mRNA by approximately two days. Thus, expression of the human embryonic gene is efficiently induced and then repressed in cultured muscle cells, as it is in muscle tissue. The study of the regulation of a human MHC isoform with a central role in muscle development and in muscle regeneration in disease states is therefore amendable to analysis at a molecular level.     

The primary structure of skeletal muscle myosin heavy chain: IV. Sequence of the rod, and the complete 1,938-residue sequence of the heavy chain

In the preceding paper [Maita, T., Miyanishi, T., Matsuzono, K., Tanioka, Y., & Matsuda, G. (1991) J. Biochem. 110, 68-74], we reported the amino-terminal 837-residue sequence of the heavy chain of adult chicken pectoralis muscle myosin. This paper describes the carboxyl terminal 1,097-residue sequence and the linkage of the two sequences. Rod obtained by digesting myosin filaments with alpha-chymotrypsin was redigested with the protease at high KCl concentration, and two fragments, subfragment-2 and light meromyosin, were isolated and sequenced by conventional methods. The linkage of the two fragments was deduced from the sequence of an overlapping peptide obtained by cleaving the rod with cyanogen bromide. The rod contained 1,039 amino acid residues, but lacked the carboxyl-terminal 58 residues of the heavy chain. A carboxyl-terminal 63-residue peptide obtained by cleaving the whole heavy chain with cyanogen bromide was sequenced. Thus, the carboxyl terminal 1,097-residue sequence of the heavy chain was completed. The linkage of subfragment-1 and the rod was deduced from the sequence of an overlapping peptide between the two which was obtained by cleaving heavy meromyosin with cyanogen bromide. Comparing the sequence of the adult myosin thus determined with that of chicken embryonic myosin reported by Molina et al. [Molina, M.I., Kropp, K.E., Gulick, J., & Robbins, J. (1987) J. Biol. Chem. 262, 6478-6488], we found that the sequence homology is 94%.     

The complete sequence and expression patterns of the atrial myosin heavy chain in the developing chick

We have earlier reported partial cloning of a cDNA of a chick atrial myosin heavy chain (MHC) gene, CCSV2 and its expression pattern in embryonic chick hearts (Oana et al (1995) Eur J Cell Biol 67, 42-49). In this study, five overlapping cDNA clones (including CCSV2) which together encode the entire open reading frame of the chick atrial MHC gene were characterized, and both the entire nucleotide sequence consisting of 5825 bases and the deduced amino acid sequence consisting of 1931 amino acids determined. Reinvestigation of the nucleotide sequence of the previously reported and presumably different chick atrial specific MHC cDNA clone, AMHC1 (Yutzey et al (1994) Development 120, 871-883), revealed that our clone and AMHC1 encoded the same MHC. The chick atrial MHC gene was strongly expressed in developing chick atria from a very early stage (Hamburger and Hamilton stage 9, 29-33 h) to the adult stage. This gene was also expressed, although weakly, in the ventricle, somite (the precursor to skeletal muscle) and skeletal muscle during embryonic stages but not in adults.     

Characterization of the carp myosin heavy chain multigene family

We isolated partial coding sequences for 29 carp myosin heavy chain genes (MyoHCs) and determined the nucleotide sequences around the region encoding the loop 2 of the myosin molecule. The predicted amino acid sequences from the isolated genes all showed very high similarity to those of skeletal and cardiac muscles from higher vertebrates, but not to those of smooth and non-muscle counterparts. Among all clones isolated, carp MyoHC10, MyoHCI-1-3 and MyoHC30 showed exon-nucleotide sequences identical to those of cDNAs encoding the loop 2 region of the 10 degrees C-, intermediate- and 30 degrees C-type fast skeletal isoforms [Hirayama and Watabe, Euro. J. Biochem. 246 (1997) 380-387]. The loop 2 of 28 types of carp MyoHCs was encoded by two exons separated by an intron corresponding to that of the 16th in higher vertebrate MyoHCs, whilst this intron was not found in carp MyoHC30. Although carp MyoHC30 had a gene organization different from those of higher vertebrates and other carp MyoHCs, its predicted amino acid sequence for loop 2 showed the highest homology to those of higher vertebrates among carp MyoHCs. In the 28 carp MyoHCs containing the intron, a combination of different nucleotide sequences for the two resulted in 14 distinct series for the combined coding sequence. These different nucleotide sequences encoded nine distinct amino acid sequences. Phylogenetic analysis for the present loop 2 and light meromyosin previously reported for carp MyoHCs [Imai et al., J. Exp. Biol. 200 (1997) 27-34] revealed that carp MyoHCs have recently diverged and are more closely related to each other than to MyoHCs from other species.     

Fractal landscapes and molecular evolution: modeling the myosin heavy chain gene family.

Mapping nucleotide sequences onto a "DNA walk" produces a novel representation of DNA that can then be studied quantitatively using techniques derived from fractal landscape analysis. We used this method to analyze 11 complete genomic and cDNA myosin heavy chain (MHC) sequences belonging to 8 different species. Our analysis suggests an increase in fractal complexity for MHC genes with evolution with vertebrate > invertebrate > yeast. The increase in complexity is measured by the presence of long-range power-law correlations, which are quantified by the scaling exponent alpha. We develop a simple iterative model, based on known properties of polymeric sequences, that generates long-range nucleotide correlations from an initially noncorrelated coding region. This new model-as well as the DNA walk analysis-both support the intron-late theory of gene evolution.     

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.     

Cloning and sequence analysis of channel catfish heavy chain cDNA indicate phylogenetic diversity within the IgM immunoglobulin family

Catfish cDNA libraries were constructed using the poly(A+) RNA obtained from in vitro stimulated catfish leukocytes. Antigenic analysis with different antisera to catfish Ig resulted in the definition of cDNA clones encoding the catfish H chain. Sequence analysis confirmed that the catfish H chain was definitively identified, based on its similarities with chicken and mouse mu chains. Two clones were each shown to encode part of the CH2 domain, the complete CH3 and CH4 domains, the C-terminus, and a 184-bp 3' untranslated region before the poly(A+) tail. The conservation of domain size and structure is clearly evident. The two cysteines forming the intradomain disulfide bridge, as well as the tryptophans located within each domain, are absolutely conserved. There are four carbohydrate acceptor sites in the catfish H chain, only one of which is phylogenetically conserved. Of the six sequenced H chain clones, one was found to differ in a single base in the CH3, which results in the loss of a carbohydrate acceptor site. Whether this difference indicates isotypic variation between closely related genes or somatic mutation is unresolved. Amino acid sequence comparisons indicate that there is a approximately 24% similarity when the catfish H chain is aligned with mouse mu chains. This is considerably less than the approximately 40% amino acid conservation found between the chicken and mouse mu chain. The amino acid sequence of the catfish H chain is most conserved in the C-terminus (approximately 30%) and the CH4 (approximately 26%); there is less conservation in the CH3 (approximately 20%) when comparisons are made with mouse mu chain. The CH3 domain of the catfish H chain also has different hydropathy properties, when compared with the CH3 domain of the higher vertebrate mu chains. Finally, the sequence of the catfish H chain indicates an unusual arrangement of the cysteines that likely participate in intersubunit and inter-H chain disulfide linkages. The disulfide linkage of these cysteines during Ig polymerization may account for the unusual covalent architecture associated with the catfish tetramer.     

Erratum: Identification of the gene for fly non-muscle myosin heavy chain: Drosophila myosin heavy chains are encoded by a gene family

[This corrects the article on p. 913 in vol. 8, PMID: 2498088.].     

Regulation of human cardiac myosin heavy chain genes: the effect of catecholamine.

The 5'-flanking regions of the alpha- and beta-cardiac myosin heavy chain (MyHC) genes were excised from the cosmid human genomic clones using Hind III and Xbal for the alpha-MyHC gene, and the Hind III and Hind III sites for the beta-MyHC gene. These fragments were linked to chloramphenicol acetyl transferase (CAT) vector to generate a chimeric fusion gene. These fusion genes were subsequently transfected to neonatal rat cardiac cultured cells to analyze the CAT activity. The alpha-MyHC gene is preferentially expressed as compared to the beta-MyHC. In the presence of norepinephrine (NE) the beta-MyHC gene is remarkably induced (within 24 hours following the addition of norepinephrine to the cardiocyte culture). However, the alpha-MyHC is also induced. Specific alpha andrenergic antagonists such as terazosin (Tz) partially suppressed both the alpha- and beta-MyHC genes as revealed by the CAT activity. These findings suggest that catecholamine does activate the human cardiac MyHC genes but does not differentiate the specific expression of either the alpha- or beta-MyHC genes.     

Divergent evolution of the myosin heavy chain gene family in fish and tetrapods: evidence from comparative genomic analysis.

Myosin heavy chain genes (MYHs) are the most important functional domains of myosins, which are highly conserved throughout evolution. The human genome contains 15 MYHs, whereas the corresponding number in teleost appears to be much higher. Although teleosts comprise more than one-half of all vertebrate species, our knowledge of MYHs in teleosts is rather limited. A comprehensive analysis of the torafugu (Takifugu rubripes) genome database enabled us to detect at least 28 MYHs, almost twice as many as in humans. RT-PCR revealed that at least 16 torafugu MYH representatives (5 fast skeletal, 3 cardiac, 2 slow skeletal, 1 superfast, 2 smooth, and 3 nonmuscle types) are actually transcribed. Among these, MYH(M743-2) and MYH(M5) of fast and slow skeletal types, respectively, are expressed during development of torafugu embryos. Syntenic analysis reveals that torafugu fast skeletal MYHs are distributed across five genomic regions, three of which form clusters. Interestingly, while human fast skeletal MYHs form one cluster, its syntenic region in torafugu is duplicated, although each locus contains just a single MYH in torafugu. The results of the syntenic analysis were further confirmed by corresponding analysis of MYHs based on databases from Tetraodon, zebrafish, and medaka genomes. Phylogenetic analysis suggests that fast skeletal MYHs evolved independently in teleosts and tetrapods after fast skeletal MYHs had diverged from four ancestral MYHs.