- Open Access
The ancient sarcomeric myosins found in specialized muscles
Skeletal Musclevolume 9, Article number: 7 (2019)
Striated muscles express an array of sarcomeric myosin motors that are tuned to accomplish specific tasks. Each myosin isoform found in muscle fibers confers unique contractile properties to the fiber in order to meet the demands of the muscle. The sarcomeric myosin heavy chain (MYH) genes expressed in the major cardiac and skeletal muscles have been studied for decades. However, three ancient myosins, MYH7b, MYH15, and MYH16, remained uncharacterized due to their unique expression patterns in common mammalian model organisms and due to their relatively recent discovery in these genomes. This article reviews the literature surrounding these three ancient sarcomeric myosins and the specialized muscles in which they are expressed. Further study of these ancient myosins and how they contribute to the functions of the specialized muscles may provide novel insight into the history of striated muscle evolution.
Striated muscles are comprised of many heterogeneous and highly specialized fibers. Skeletal muscle fibers are extremely adaptable and can meet the varying demands of muscles by responding to changes in environmental cues in several ways. These muscle fibers can remodel their structure and contractile properties, largely by reprogramming the gene expression profiles of sarcomeric components including myosin motor proteins. Myosin produces the force necessary for a variety of cellular movements by hydrolyzing ATP and interacting with actin filaments. The myosin superfamily encodes 18 distinct classes of myosin motors, which are found ubiquitously in eukaryotes and participate in a variety of cellular motile processes (see Hartman and Spudich  for review). A subset of class II myosins, those that power muscle contraction in striated muscles, will be the focus of this review. The class II myosin heavy chain (MyHC) protein is comprised of two functional units, a globular motor domain that contains the catalytic ATPase site and binds actin and an α-helical coiled-coil rod domain that dimerizes and assembles into bipolar thick filaments (Fig. 1a). Striated muscle has a highly organized ultrastructure, consisting of a repeating contractile unit called the sarcomere. Within the sarcomere, myosins assemble into bipolar thick filament structures where the myosin motor domains protrude off the surface to interact with neighboring actin filaments (Fig. 1b). Myosin thick filaments and actin thin filaments interdigitate and slide past one another to cause sarcomere shortening in a concerted motion resulting in muscle contraction. This review describes the most ancient sarcomeric myosins, which have largely escaped characterization due to their unique expression patterns, variable expression levels across species, and relatively recent annotation, and provides an overview of the specialized skeletal muscles in which they function.
Sarcomeric myosin genes
Vertebrates express 11 sarcomeric myosin heavy chain (MYH) genes in their striated muscles (Fig. 2); however, the expression level and pattern varies greatly across species, developmental timepoint, and muscle type [2, 3]. All 11 MYH genes are expressed in mammalian striated muscles, albeit some in extremely low abundance and some localized to just one or two specialized skeletal muscle fibers. These genes, their protein products, and muscle type of expression are listed in Fig. 2. Though humans encode MYH4 and mRNA has been detected in extraocular  and jaw muscles , as well as in cases of Duchenne muscular dystrophy , this isoform is not expressed at the protein level in humans. Furthermore, MYH16 is a pseudogene in humans and does not produce any functional protein . However, other mammals express MYH4 and MYH16 at the protein level in fast-twitch skeletal muscle fibers (known as Type II fibers), alongside the other skeletal muscle-specific myosin isoforms [8, 9]. MYH6 (α-MyHC) and MYH7 (β-MyHC) are known as the cardiac myosin isoforms and are expressed in the mammalian heart. However, MYH7 is also the dominant isoform expressed in slow-twitch skeletal muscle fibers (also known as Type I fibers), and both MYH6 and MYH7 are found in certain specialized skeletal muscles. Finally, recent genomic analysis identified MYH7b and MYH15, which are expressed in various specialized muscles of mammals such as the extraocular muscles and muscle spindles. Sequence identity is extremely high among the human sarcomeric myosin genes (Fig. 3), and yet, each myosin’s contractile properties, such as ATP turnover and force production, vary from isoform to isoform [10,11,12]. These variations in motor properties confer unique characteristics to different muscle types, and the varied myosin composition within skeletal muscle fibers allows for a wide range of contractile velocities and forces among different muscle types [11, 13]. For example, MyHC-embryonic and MyHC-perinatal are expressed in developing skeletal muscles when the demands on muscle are lower due to decreased load . In adult skeletal muscle, Type II fibers express varying ratios of the adult skeletal muscle myosin isoforms, which have relatively high ATP-turnover rates and confer fast contractile properties to these fibers . In contrast, β-MyHC, the major human cardiac isoform also found in slow Type I skeletal muscle fibers, is the slowest ATPase of all the human sarcomeric myosin isoforms characterized to date [11, 14]. Thus, sarcomeric myosin genes have adapted to have diverse functions, which meet the demands of various muscle types as they have evolved.
In addition to the high amino acid sequence conservation of sarcomeric myosins, the organization of the cardiac and skeletal muscle myosins in two tandemly linked genomic clusters is conserved across mammalian species . In the human genome, a 350 kilobase (kb) segment on chromosome 17 contains the six skeletal muscle myosin genes (MYH3, MYH2, MYH1, MYH4, MYH8, MYH13, in order of tandem linkage) . The two cardiac myosin isoforms, MYH6 and MYH7, are in a second cluster on chromosome 14, separated by 4.5 kb . Conversely, the three most recently discovered ancient myosins are physically unlinked to any other MYH gene in the genome; MYH7b is encoded on chromosome 20, MYH15 is found on chromosome 3, and MYH16 is on chromosome 7 . The broad conservation and clustering of the MYH genes suggests two things. First, this genomic organization is important for the genes’ regulation, though they are not organized temporally . Secondly, they suggest that this gene family resulted from gene duplication events of ancestral myosins [15, 17].
Discovery of the ancient myosins
Prior to the technological advances that define the genomics age, the major sarcomeric myosins expressed in skeletal and cardiac muscle were extensively studied and characterized due to their abundance in and accessibility of these tissues. However, as the human genome was being annotated, three additional class II myosin genes were discovered: MYH7b,Footnote 1 MYH15, and MYH16. In 2002, Desjardins et al. identified these three novel MYH sequences and determined their genomic locations by mapping these genes to cDNA databases. Deduced protein sequences of these genes revealed conserved sequence motifs and homology to known sarcomeric myosins, indicating that these genes are class II sarcomeric myosins . The same study predicted each myosin motor’s contractile speed by comparing the sequences to characterized sarcomeric myosin motor domains; MYH7b and MYH15 were suggested to be slow isoforms whereas the protein encoded by MYH16 (MyHC-masticatory) was predicted to be a fast isoform. However, activity assays have since demonstrated that MyHC-masticatory is more forceful rather than fast . Molecular evolutionary analysis of the three novel myosins indicated that these genes are ancient and predate the well-studied skeletal and cardiac isoforms and the divergence of a smooth muscle MYH gene . In fact, the ancient myosins exhibit a lower sequence identity in their full-length sequences and individual motor and rod domains to the sarcomeric class II myosins than any other sarcomeric myosin (Fig. 3), supporting the notion that they are more distantly related to the well-characterized skeletal and cardiac myosin isoforms. Since the initial discovery of the three ancient myosins, orthologs have been found in distant species including fish, chickens, snakes, and frogs, indicating that these ancient myosins were present in a common ancestor of vertebrates [19,20,21]. Interestingly, MYH7b and MYH15 play a prominent role in the heart and skeletal muscles of certain species like chicken and snakes, whereas in mammals, the only muscles that these myosin motors function in are highly specialized muscles (Leinwand unpublished ). The MYH16 gene is expressed exclusively in muscles that originate from the first pharyngeal arch; MYH16 is primarily expressed in the muscles of mastication but is also present in the tensor veli palatini and tensor tympani of certain species . While MYH16 expression has been observed in the jaw muscles of some vertebrates including cats  and crocodiles , a frameshift mutation led to the loss of MYH16 expression in humans . Though these ancient myosins are less studied, they are clearly set apart from the striated isoforms that have been extensively characterized and shown to play major roles in mammalian heart and skeletal muscle function. Little is known about why these myosins are absent from conventional striated tissues of mammals, unlike in reptiles and birds, and the evolutionary pressures that resulted in differential expression of these myosins across species. The three ancient myosins will be further discussed in the subsequent sections of this review.
Evolutionary perspective of the myosin genes
In order to understand the evolutionary relationships between the myosin genes within and across species, it is important to consider a broad perspective of muscle evolution itself. Smooth and striated muscle cells are unique to specific members of the animal kingdom [25, 26]. Originally, animal striated muscles were presumed to share a common origin based solely on the ultrastructural similarity of the highly ordered striated muscle tissue across phyla [27, 28]. However, molecular phylogenetic analysis encompassing animals, fungi, plants, and protists now strongly suggests that striated muscles are a result of convergent evolution. In one study, Steinmetz and colleagues discovered a core set of muscle proteins in organisms that predate the evolution of multicellular organisms . Previously, researchers hypothesized that a gene duplication event gave rise to two distinct MYH orthologs in bilaterians, accounting for the presence of both visceral smooth muscle and somatic striated muscle [30, 31]. More recent genome mining revealed that this duplication event occurred before the origin of muscle cells and that at least two distinct myosin isoforms were present in a shared common ancestor of all animals . These two ancient myosin isoforms are referred to as SM-MHC and ST-MHC and have since diversified from each other with respect to assembly, contractile function, and cell type . Smooth and nonmuscle myosin orthologs resulted from duplications of SM-MHC, while ST-MHC underwent duplication events to produce the current suite of sarcomeric myosin isoforms present today including MYH7b, MYH15, and MYH16, which were the first sarcomeric myosin genes to diverge from the ancestral myosin genes (Fig. 2) [17, 29].
These three ancient myosins are expressed in mammalian specialized muscles (Fig. 4), which is logical as the majority of these tissues are thought to predate conventional striated muscles. In this review, we define specialized muscles as skeletal muscles that have adopted a unique structural organization to perform a specific function, while cardiac muscle and trunk/limb skeletal muscles are referred to as conventional muscles. Specialized muscles include extraocular muscle, muscle spindles, and masticatory muscle, which will be further discussed in this review . It is worth noting that the intrinsic laryngeal muscles are also considered to be specialized muscles, but are not known to express any of the ancient myosins and are therefore outside the scope of this review. The three ancient myosins appear to have niche roles in specialized mammalian muscles, but given their diverse expression levels and variable abundance across species, several questions remain. Are these myosins required for the specific function of the specialized muscles in which they reside? Could lack of selective pressure lead to the replacement of these myosins by a more diversified set of sarcomeric myosins? Can these myosins aid our understanding of early striated muscle evolution? In this review, we will first introduce the specialized muscles and then turn our focus to MYH7b, MYH15, and MYH16 and their roles in these specialized muscle types.
The specialized muscles
Extraocular muscles (EOMs) are the muscles that surround the eyeball and control a variety of diverse and complex eye movements. These movements are necessary for eyesight precision and range from fast movements to slower pursuit and vergence movements to fixation . Six EOMs control eye movement across vertebrates including four rectus muscles (medial, lateral, inferior, and superior) and two oblique muscles (superior and inferior). Each EOM originates from the posterior of the eyeball near the optic nerve and attaches to the front of the eyeball with the exception of the inferior oblique muscle . The medial and lateral rectus muscles control horizontal eye movements. Vertical eye movement is achieved through the coordinated action of the superior and inferior rectus and oblique muscles, the latter of which also function in torsional movements [34,35,36]. Mammals possess an additional principal EOM called the levator palpebrae, which controls eyelid elevation . Each EOM compartmentalizes into two layers, the thin orbital layer that lies along the EOM surface facing the bony orbital wall and the thicker global layer that faces the eyeball globe . The orbital layer and global layer have distinct features as does the marginal zone, a third layer described in human EOM that lies on the orbital layer’s outer surface . Within these layers, there are six distinct fiber types that participate in all eye movements, and they are classified by fiber location (global or orbital layer), innervation (singly or multiply innervated fibers), and histochemical features [34, 37].
Mammalian EOMs express an unusually high number of striated muscle myosin genes including developmental and adult isoforms as well as an EOM-specific isoform. In total, mammalian EOMs express 10 of the 11 sarcomeric myosin isoforms: the adult skeletal myosins, MyHC-IIa, MyHC-IIx/d, and MyHC-IIb; the developmental isoforms, MyHC-embryonic and MyHC-perinatal; the cardiac isoforms, β-MyHC and α-MyHC; the EOM-specific isoform, MyHC-extraocular; and two of the ancient myosins, MYH7b (Fig. 4a) and MYH15 [21, 39,40,41,42,43,44]. Interestingly, the extraocular and laryngeal muscles are the only tissues that express MYH13, a gene first identified by Wieczorek et al. and later discovered to evolutionarily predate the fast skeletal isoforms [45,46,47]. The MYH13 protein product, MyHC-extraocular, was categorized as a fast myosin and observed in the innervation zone of both the EOM orbital and global layers, likely contributing to fast eye movements [47, 48]. Myosin isoform composition within the EOMs varies with respect to the orbital and global layers and within the muscle fiber longitudinally [41, 43, 49]. EOMs also exhibit heterogeneity in myosin composition across species. In rodents, specific MYH isoforms appear to be preferentially expressed in a certain layer , whereas in humans most MYH isoforms are observed in both layers . Furthermore, many EOM fibers co-express multiple myosins, resulting in extremely complex expression patterns across species [43, 44, 48,49,50]. Based on this complex expression profile, it seems that EOM function cannot be assigned to a specific layer or set of myosin isoforms. Thus, it is likely that the comprehensive myosin expression pattern in the EOMs is critical for maintaining the demanding and diverse eye movements.
Gene expression profiling has revealed distinct differences between EOMs and other striated muscles, indicating that EOMs are highly specialized muscles. DNA microarray was used to identify differential gene expression between EOM and representative cardiac, skeletal, and smooth muscles, revealing that EOM is substantially different than other striated muscles . Expression data also revealed that EOM has a distinct energy metabolism profile where genes related to glycogen metabolism are downregulated relative to skeletal muscles [51, 52]. Gene expression profiling has also identified a regulatory gene found preferentially in EOMs called the paired-like homeodomain transcription factor 2 (Pitx2) [52, 53]. Unlike in conventional skeletal muscles, Pitx2 is required for EOM differentiation and development in mice and appears to regulate gene expression related to MYH isoforms, contractility, and EOM fiber size [53,54,55,56]. Pitx2 is also highly expressed in a population of EOM myogenic precursor cells in healthy and muscular dystrophy mice, suggesting a role for Pitx2 in the EOM’s regenerative capabilities that may contribute to preferential EOM sparing in certain diseases .
One yet unresolved question in the field pertains to why EOM is preferentially involved in or spared in certain diseases. The EOMs are preferentially impacted in myasthenia gravis, an autoimmune disease that targets acetylcholine receptors of the neuromuscular junction. This is hypothesized to be attributable to the differences in the number and isoforms of acetylcholine receptors and innervation properties at the neuromuscular junction of EOM versus typical skeletal muscles [34, 58]. Likewise, differences in the immune response between EOMs and skeletal muscle likely underlie the selective involvement of the EOM in the autoimmune disease Graves’ ophthalmopathy [59, 60]. Conversely, despite the widespread muscle degeneration that characterizes muscular dystrophy, the EOMs are spared in this disease and patients retain normal EOM function . Expression profiling revealed that there are no clear changes in gene expression in EOMs of a muscular dystrophy mouse model, indicating that EOM remains unaffected in the disease state. Thus, the characteristics that distinguish EOM from conventional skeletal muscles may be protective in this disease [62, 63]. A current hypothesis in the field predicts that a population of myogenic precursor cells mediates the EOM’s regenerative capacity and contributes to EOM sparing in muscular dystrophy . In fact, EOMs of aged and muscular dystrophy mice retain an enriched population of muscle precursor cells with proliferative capacity compared to limb skeletal muscle [64, 65]. Though the precise mechanism remains undefined, the EOM’s susceptibility and protection from certain diseases likely stems from a combination of its unique features, including distinct myosin isoform expression patterns, fiber types, gene expression profile, and developmental requirements.
Muscle spindles are the innervated sensory structures buried within conventional striated muscles that relay information about muscle length and stretch to the central nervous system [66,67,68]. Muscle spindles provide essential sensory input about the spatial placement, extension, and contraction rate of muscles to the central nervous system, which is required for effective voluntary control of striated muscles [66, 69]. They provide a static and a dynamic stretch response by detecting both the initiation and the continuation of stretch, respectively [66, 69, 70]. Muscle spindles are differentially distributed throughout the body’s muscles and are typically focused along axial regions and in smaller muscles, though the reason for this distribution is unclear [71, 72]. Mice deficient in Egr3, a transcription factor essential for muscle spindle development, have severe phenotypes including loss of voluntary muscle control, scoliosis, resting tremors, and eyelid sagging, even though the extrafusal muscle fibers are normal . Thus, muscle spindles are necessary for normal locomotion.
The general muscle spindle morphology consists of a bundle of intrafusal muscle fibers that are separated from the extrafusal fibers of the surrounding skeletal muscle . There are three types of intrafusal muscle fibers: bag1, bag2, and chain fibers [66,67,68, 73]. Bag and chain type fibers are distinguished by their nuclear distribution within the muscle fiber; bag fibers have centrally located nuclei, whereas chain fiber nuclei are equally distributed along their length . In addition, electron microscopy studies have shown distinct sarcomeric, myofibril, and mitochondrial organization between these two muscle spindle fiber types. Chain fibers have clearly defined myofibril units and numerous, large mitochondria that span one to two sarcomeres each (Fig. 4b) [74, 75]. The organization in bag fibers is less clear, with less well-defined sarcomeres that are tightly packed and have fewer, smaller mitochondria than chain-type fibers (Fig. 4c) [74, 75]. Bag1 fibers are distinguished from bag2 fibers by a lower ATPase activity and contain a larger proportion of slow myosin isoforms such as β-MyHC [66, 76]. Functionally, bag1 fibers are solely responsible for the dynamic stretch response, whereas both bag2 and chain fibers contribute to the static stretch response [66, 67]. Stretch information, and therefore proprioception, is passed through the afferent fibers that innervate the muscle spindles to the central nervous system [66, 67].
The muscle spindle acts as the primary stretch receptor of the muscle to relay feedback information to the nervous system, which in turn regulates the length and tone of the extrafusal fibers. The myosin-actin crossbridge is essential to this muscle spindle function. The initial stretch response recorded by the spindle is thought to be the result of breaking static crossbridges by stretching the muscle [66, 67, 77]. This causes changes in spindle stiffness, triggering mechanosensitive ion channels to initiate a signal cascade [67, 77]. The crossbridge reformation rate also creates signals downstream of the original stretch, which influences whether or not the muscle remains extended . Myosin composition therefore plays an important role in muscle spindle response rate. There is a wide range of sarcomeric myosins expressed in muscle spindles across species, including β-MyHC, α-MyHC, MyHC-embryonic, MyHC-perinatal, and MyHC-IIa . Two of the ancient myosins, MYH7b and MYH15, have also been found in muscle spindles . Keeping in mind that bag1 and bag2 fibers are differentiated by their variable ATPase content, it is probable that the myosin expression profile of each intrafusal muscle fiber is an important determinant of its activity . The need to finely tune the crossbridge dynamics of different muscle spindles present in the body could explain why the two ancient myosins, MYH7b and MYH15, are expressed in muscle spindles. It is possible that the expression levels of these ancient myosins contribute to the precise contractile needs of muscle spindles in order to garner a robust stretch response in different muscles.
The masticatory muscles move the jaw to accomplish a wide range of activities including swallowing, chewing, speaking, and infantile suckling. The diverse functions of the jaw muscles are in part due to their unique structure, which allows the muscle to exert precise control over the mandible (i.e., detached lower jaw bone). In humans, four muscles comprise the masticatory apparatus: the masseter, temporalis, medial pterygoid, and lateral pterygoid. The first three muscles function to elevate the jaw, while the lateral pterygoid is a jaw-opening muscle . The masseter muscle is the most superficial masticatory muscle, originating from the zygomatic arch (i.e., cheekbone) and attaching to the mandible. The masseter is the most powerful of the four muscles allowing for forceful jaw closures. The temporalis muscle originates in a shallow depression on the side of the skull known as the temporal fossa and converges into a tendon that attaches to the mandible. The temporalis muscle can retrude the jaw as well as close the jaw. The last of the jaw-closing muscles is the medial pterygoid. Both the medial pterygoid and lateral pterygoid are attached to the inner surface of the jaw and can protrude the mandible. The lateral pterygoid is responsible for moving the jaw laterally and is also involved in protruding and depressing the mandible .
The masticatory unit in humans has evolved to be a highly specialized muscle group. The complex muscle fiber architecture, heterogeneity, and myosin composition of the masticatory muscles allow for a wide range of forces and contractile velocities. The mastication muscles are made up of a variable population of pure and hybrid muscle fibers. Pure fibers expressing a single myosin isoform are found for slow β-MyHC, MyHC-IIa, and MyHC-IIx/d. Hybrid fibers expressing a combination of these isoforms, as well as MyHC-embryonic and α-MyHC are also present . Immunohistochemistry studies on male and female human cadavers demonstrated that all four masticatory muscles express approximately 70% β-MyHC. The temporalis and masseter muscles show a relatively equal distribution of the remaining MyHCs present, while the pterygoids show a marked increase in MyHC-IIa and α-MyHC compared to MyHC-IIx/d and MyHC-embryonic . MYH4 transcripts are abundantly expressed without accompanying MyHC-IIb protein expression , similar to previous reports of MYH4 expression in other human skeletal muscles [6, 8]. The fiber-type population within the four jaw muscles varies greatly, as does fiber size itself [80, 81]. Interestingly, differences in fiber composition between males and females have been observed in the masseter muscle of rabbits  and mice . In humans, a significant increase in the cross-sectional area (CSA) of masseter Type I and Type II fibers has been observed in males compared to females [84,85,86]. Hence, the masseter muscle appears to be sexually dimorphic in certain species of mammals. Diverse expression of specific myosin motors in the masticatory muscles allows for energy consumption at the cellular and animal level to be fine-tuned via the non-redundant ATPase properties inherent to different isoforms . This complex sarcomeric myosin expression pattern suggests that the masticatory muscle is highly adaptable. Indeed, changes in fiber type, CSA, and the contractile mechanics of jaw muscles have been observed in response to diet [88,89,90,91,92,93,94,95,96,97], aging [98,99,100,101,102,103,104,105], and craniomandibular disorders .
Alterations in craniofacial morphology have been correlated with the composition and function of temporalis and masseter muscles, but it is unclear whether the jaw muscles are the driving force behind craniofacial morphological alterations or if the muscle function adapts to changes in external stimuli from the skeleton. There are three basic facial forms categorized for humans: long, average, and short . There is greater facial morphology variation in individuals with weak jaw muscles , while stronger and thicker jaw muscles produce more uniform facial morphological features suggesting that if the masticatory muscle is too weak, it exerts less influence over facial morphology . The jaw muscles are hypothesized to determine facial dimensions of humans . The jaw-closing muscles of short-faced individuals have thicker fibers [109, 110], while those of long-faced individuals have thinner fibers [111, 112], and produce lower molar bite forces [113,114,115]. Differences in fiber type are also evident in people with different facial forms. Rowlerson et al. observed that patients with an open bite presented with a higher proportion of slow Type I fibers, while patients affected by deep bite presented with a higher proportion of fast Type II fibers . Boyd et al. observed that a longer face correlated with more Type II fibers in masseter muscles . In addition, several animal studies indicate that the jaw muscles may play a role in skull growth and development (for review see Kiliaridis ). Surgical manipulation of the temporalis muscle in young rabbits directly resulted in changes to the local skull morphology and cranial development, although this was not correlated with any specific myosin isoform . A similar study in primates revealed that altering the temporalis and masseter muscle attachments could alter craniofacial growth patterns . Thus, the masticatory muscles are a complex unit that attach to the skull and jaw and can influence the phenotypic characteristics of the cranium.
The jaw made its first appearance in the evolutionary timeline over 400 million years ago when gnathostomes, or jawed vertebrates, first diverged . It is thought that the jaw evolved through adaption of ancient gill cartilages . The predatory success of early gnathostomes is largely attributed to the presence of the jaw , and today, gnathostomes make up 99% of vertebrate species. With the radiation of vertebrate and mammalian speciation, the jaw has had ample time and opportunity to evolve to the functional needs of a wide range of animals. The MYH16 gene is the oldest of the sarcomeric myosins present in vertebrates and is primarily expressed in jaw muscles. MyHC masticatory, which is associated with a high contraction force, has been observed in certain carnivorous species, such as sharks and cats (Fig. 4d, e) , where forceful jaw closures would be advantageous for predation. The jaw muscles have evolved to serve functions that are distinct from conventional skeletal muscle (locomotion and postural tonicity). By expressing a wide variety of myosins, masticatory muscle is able to regulate the force, speed, and energy efficiency of contraction in order to adapt to the ever-changing environments and demands of various animals.
The ancient myosins
Human MYH7b is a 27-kb gene found on chromosome 20. MYH7b shares the highest sequence identity to MYH7 (β-MyHC) and MYH6 (α-MyHC), reaching 69% sequence identity at the amino acid level (Fig. 3). In addition, each of these three myosin genes harbors an intronic microRNA; MYH6 encodes miR-208a in intron 27, MYH7 encodes miR-208b in intron 31, and MYH7b encodes miR-499 in intron 19 [124, 125]. miR-208a is cardiac myocyte-specific and has known roles in cardiac stress response regulation, whereas miR-208b and miR-499 play regulatory roles in the heart and have redundant roles in skeletal muscle fiber-type specification [124, 125]. Though MYH7b shares a high sequence homology and features with the two human cardiac myosin isoforms, this myosin’s expression pattern and regulation are unique in mammals.
MYH7b was first identified by Nagase et al. in an effort to categorize previously unknown human genes . The MYH7b sequence (designated KIAA1512) was cloned from a human fetal brain cDNA library and mapped to chromosome 20. Expression profiling by the authors revealed that MYH7b RNA is highly expressed in the heart, skeletal muscle, adult, and fetal brain and more lowly expressed in the ovary, kidney, lung, liver, pancreas, and spleen . As the human genome was being annotated, Desjardins et al. identified MYH7b as one of three ancient myosins belonging to the sarcomeric myosin family of genes . MYH7b transcripts were confirmed in several conventional mammalian muscles including the heart, soleus, tibialis anterior, quadriceps, and diaphragm [21, 127]. A 2012 study by Warkman et al. reported the presence of MYH7b protein in the mouse heart; however, this finding was later retracted as this result was due to non-specific antibody reactivity . To date, MYH7b has not been detected at the protein level in mammalian cardiac or skeletal muscle. This discrepancy between MYH7b RNA and protein expression in conventional muscles is due to a non-productive splicing event in which the transcript is produced, but undergoes some nonsense mediated decay, while miR-499 expression is retained in these tissues . Despite this expression pattern in the heart and skeletal muscle, MYH7b protein is present in specialized muscles of mammals including the EOMs and muscle spindles. In rats, MYH7b was detected in EOM orbital and global layers (Fig. 4a) and observed in varying abundance in muscle spindle bag1 and bag2 fibers . Both of these specialized muscles express many sarcomeric myosins and carry out highly specific tasks. Thus, MYH7b may have a role that uniquely fits a contractile requirement of these muscles.
MYH7b transcripts capable of encoding protein have also been observed in the mouse brain . Corroborating this finding, Rubio et al. identified MYH7b in a yeast two-hybrid screen designed to identify actin regulators in the rat brain. Subsequently, this group used siRNA to target MYH7b in cultured hippocampal neurons and observed alterations in dendritic spine morphology and excitatory synapse strength, indicating that MYH7b contributes to dendrite structure and function . MYH7b is also implicated in several diseases that involve both muscle and nonmuscle phenotypes. In 2014, Haraksingh et al. reported that compound heterozygous mutations in MYH7b are associated with sensorineural hearing loss. Exome sequencing on a family with sensorineural hearing loss in three children revealed two mutations in MYH7b: a maternally inherited mutation in the MYH7b motor domain and a paternally inherited mutation in the MYH7b rod domain . These are the only known associations between MYH7b and deafness, and thus is an unusual occurrence. Whole exome sequencing has also identified missense mutations in MYH7b associated with left ventricular non-compaction cardiomyopathy . Recently, a genome-wide association study (GWAS) identified associations between congenital heart defects involving left-sided lesions and chromosome 20q11, a region that includes MYH7b . Another GWAS aiming to identify new cutaneous melanoma risk loci identified significant associations with 20q11 and MYH7b . However, the variants associated with MYH7b are hypothesized to be positional markers for genes associated with melanoma-like ASIP, a pigmentation gene in the same chromosomal region as MYH7b, rather than a disease modifying variant [134, 135].
Though MYH7b does not appear to play a major role at the protein level in conventional muscles of mammals, MYH7b orthologs are found in fish, chickens, snakes, and frogs and appear to have predominant roles in the conventional muscles of these species. In chickens, the MYH7b ortholog is expressed in the Purkinje fibers of the heart, as well as slower skeletal muscles such as the anterior latissimus dorsi (Leinwand unpublished [21, 136]). The fish genome also contains MYH7b orthologs, and it appears that gene duplication has resulted in multiple MYH7b paralogs [19, 20]. For example, torafugu have two gene duplicates of the MYH7b ortholog; one is specifically expressed in slow skeletal muscle whereas the other has broader expression in adult and embryonic tissue . Syntenic organization of MYH7b is consistent in Xenopus, chickens, mice, and humans, further validating the concept that this myosin is derived from an ancient myosin found in an ancestor common to these species [21, 137].
MYH15 is encoded on human chromosome 3. Its genomic sequence is unusually large, around 142 kb, and genomic analysis led to the conclusion that it is one of three ancient myosins found outside of the canonical sarcomeric myosin clusters . Since the identification of MYH15 in humans, orthologs have been found in chickens, snakes, and various mammals, though no such ortholog has been found in teleost fish [20, 138]. The MYH15 ortholog found in chickens is the major myosin found in the adult ventricle and is also important for embryonic development (Leinwand unpublished ). These data suggest that MYH15 acts as a conventional sarcomeric myosin. However, the expression pattern of MYH15 in mammals does not follow that of a typical sarcomeric myosin. Instead, MYH15 transcripts have been observed in several brain regions, the testis, and the pituitary gland [139, 140]. MYH15 RNA has also been described as a long non-coding RNA found in the ovaries of sheep and is suggested to be important for fecundity in these animals [141, 142]. MYH15 protein has been identified in mammalian muscle tissue via immunostaining in rat EOMs and muscle spindles . However, in humans, MYH15 protein has only been observed in the pulmonary vascular epithelium and alveolar macrophages . This unusual expression pattern in mammals suggests a role for MYH15 outside of muscle.
The role of MYH15 in muscle spindles is unclear. Muscle spindle myosin expression profiles vary regionally throughout the spindle, and isoform variation influences the myosin-actin crossbridge dynamics known to be important for the response of the muscle spindle [66, 68, 76, 77]. Thus, the predicted slow ATPase rate of MYH15 may contribute in some way to the specificity of muscle spindle signaling. There are further clues for this myosin’s function in the pulmonary vascular system, as disease-associated single-nucleotide polymorphisms (SNPs) in MYH15 correlate with noncardioembolic stroke, coronary heart disease, and chronic obstructive pulmonary disease in humans [143,144,145,146,147]. In addition, a SNP in MYH15 has been linked to bovine pulmonary hypertension, supporting a role for MYH15 in the pulmonary system . The associated diseases in humans are all caused by obstructions in the vascular system, which suggests that MYH15 may have a role in clearing such blockages. Even though MYH15 shares the heptad repeat rod motif universal to sarcomeric myosins and is the major ventricular MyHC in chicken hearts, the identified expression pattern in mammals suggests the potential for a role outside of the sarcomere . Given that ancestral myosins were present before the evolution of a muscle cell, it is reasonable to contemplate that the first myosins to diverge (i.e., MYH15) could retain a contractile or motile role outside of a muscle cell.
Though the human MYH16 gene is a pseudogene, it has 42 predicted exons and encompasses 68 kb, spanning three times the physical size of a typical sarcomeric myosin . Phylogenetic analyses revealed that MYH16 predates all other vertebrate MYH isoforms and is the most distantly related MYH gene compared to all others in vertebrates [17, 149]. The human MYH16 sequence has five additional introns that are only shared with the striated muscle gene of mollusks, further supporting the notion that MYH16 is one of the most ancient myosins in the animal kingdom . Zhu et al. showed that MYH16 and its syntenic region are not present in the mouse genome, indicating MYH16 was lost in mice due to a genomic deletion . In fact, the MYH16 gene has been independently lost from several other mammalian lineages [22, 151, 152], indicative of the susceptibility of jaw-closing muscles to selective pressures . A trend of MYH16 loss is seen in mammals that typically favor faster contractions over forceful jaw contractions, which is often true for smaller animals. MYH16 produces protein in many other vertebrate species where it is expressed in masticatory muscle. More specifically, MYH16 orthologs have been found in the masseter and temporalis muscles of certain primates [22, 151], carnivores [9, 22,23,24], marsupials , rodents , bats , crocodiles , and sharks .
Human MYH16 became a pseudogene due to a two-nucleotide frameshift deletion in exon 18 that results in a premature stop codon . Stedman et al. predicted the MYH16 mutation occurred 2.4 million years ago (MYA), under the assumption that the neutral mutation rate remained constant since the human-chimp divergence (6–7 MYA) . According to this timeline, the genetic inactivating mutation is predicted to precede the reduction in mass of the mandible in the human ancestor, Homo erectus, which occurred around 2.0 MYA  simultaneously with a shift toward a larger cranium . This led many to hypothesize that the loss of MYH16 in early Homo resulted in small jaw muscles and less tension applied to the cranium during development, allowing for increased cranial capacity and encephalization. Mechanical measurements of pure MyHC-masticatory fibers isolated from the masseter and temporalis muscles of dogs and cats showed that MyHC-masticatory was the most forceful of all the myosin proteins [18, 157]. Contrary to the prediction of MyHC-masticatory constituting a “superfast” myosin , one study showed that MyHC-masticatory fibers had an intermediate contraction velocity, most similar to that of MyHC-IIa . The notion that loss of MYH16 corresponded to the loss of the most forceful myosin in the human jaw has lent itself to the hypothesis that encephalization resulted from reduced tension in the jaw. A more comprehensive analysis of MYH16 evolutionary history later called into question the proposed date of the MYH16 mutation, suggesting the gene may have been inactivated 5.3 MYA, which precedes the first appearance of Homo and the modern day cranium . In 2015, the same authors revisited this question by analyzing nuclear genomic sequences of multiple hominin lineages and concluded that the inactivating mutation occurred before the Human-Neanderthal divergence 0.6 MYA (as they share this mutation) , and after the human-chimpanzee divergence, consistent with the 2.4 MYA estimate by Stedman et al. .
Despite this timeline, some have argued that MYH16 inactivation would not contribute to an increased cranial capacity, mainly because the majority of brain growth precedes the formation of the masticatory apparatus during fetal development . Other researchers have instead proposed that smaller jaws were the cause of the mutation. A dietary shift to consuming softer foods could have resulted in smaller jaw muscles, leading to a lower dependency on the MYH16 isoform. The lack of selective pressure could have in turn resulted in the loss of MYH16 [9, 160]. Nevertheless, the loss of MYH16 in humans has been proposed to be associated with a marked reduction in masticatory muscle mass, which could have allowed for human encephalization [7, 162]. While loss of the MyHC-masticatory in humans might not have directly resulted in encephalization, it is a fascinating example of how genetic inactivation can lead to the acquisition of human-specific phenotypes .
The discovery of the three ancient myosins, MYH7b, MYH15, and MYH16, completed the inventory of the mammalian sarcomeric myosin genes. For decades, these myosins escaped detection due to their limited expression in the specialized tissues of mammals including the EOMs, muscle spindles, and the masticatory muscles. Even 17 years after their annotation in human genome, many questions about the ancient myosins and their role in specialized muscles remain unanswered. One such question is whether the ancient myosins have an essential role in specialized muscles or whether they are being eliminated through evolutionary processes in mammals. Take, for example, one of the characterized sarcomeric myosins, MYH13 (MyHC-extraocular), that is expressed in two highly specialized muscles (EOM and laryngeal muscle). Given that MyHC-extraocular has the fastest contractile rate of the sarcomeric myosins and is expressed exclusively in muscles that require superfast contractions (e.g., eye movements and phonation), it appears that this myosin has evolved to fulfill unique roles in these specialized muscles. Therefore, the argument can be made that MYH13 expression is required for proper muscle function in specialized muscles. Likewise, MYH16 may too be required for the proper function of the specialized masticatory muscles in certain species given its exclusivity in this tissue type and its high force-generating capability compared to the other sarcomeric myosins. A different argument could be made for MYH7b and MYH15 function in the EOM and muscle spindles. These specialized muscles are known to contain multiple myosin isoforms, which is hypothesized to allow the muscle to be functionally adaptable. In contrast to the idea that each myosin performs a specific task, perhaps this varied expression of myosins, including the ancient ones, is the basis for their plasticity and adaptability wherein there is no reliance on one isoform for a specific function. Lastly, certain species may selectively express the ancient myosins due to muscle type and demand, as is seen with the well-characterized myosins. Preferential expression of these ancient myosins at the protein level in conventional muscle of more distantly related species (e.g., snakes and birds that express MYH7b and MYH15 protein in their cardiac and skeletal muscle) contrasts with that of mammals in which these myosin motors are restricted to specialized muscles. The majority of studies encompassing the ancient myosins are done in mammals, but there may be more to be learned from studying these isoforms in diverse species. Doing so may help answer the question of whether these myosins require niche roles to remain evolutionarily relevant. Finally, a major question remains as to whether vertebrates will continue to evolve more functionally distinct myosins to satisfy the ever-changing needs of the muscle and whether the functional repertoire of myosins will grow or if certain myosin isoforms will become obsolete in the future.
The nomenclature at the time designated this myosin isoform as MYH14; however, current nomenclature now refers to this myosin as MYH7b to avoid confusion with the gene encoding nonmuscle myosin IIc.
Million years ago
Myosin heavy chain (gene)
Myosin heavy chain (protein)
Paired-like homeodomain transcription factor 2
Hartman MA, Spudich JA. The myosin superfamily at a glance. J Cell Sci. 2012;125:1627–32.
Mascarello F, Toniolo L, Cancellara P, Reggiani C, Maccatrozzo L. Expression and identification of 10 sarcomeric MyHC isoforms in human skeletal muscles of different embryological origin. Diversity and similarity in mammalian species. Ann Anat. 2016;207:9–20.
Schiaffino S, Rossi AC, Smerdu V, Leinwand LA, Reggiani C. Developmental myosins: expression patterns and functional significance. Skelet Muscle. 2015;5:1–14.
Horton MJ, Rosen C, Close JM, Sciote JJ. Quantification of myosin heavy chain RNA in human laryngeal muscles: differential expression in the vertical and horizontal posterior cricoarytenoid and thyroarytenoid. Laryngoscope. 2008;118:472–7.
Horton MJ, Brandon CA, Morris TJ, Braun TW, Yaw KM, Sciote JJ. Abundant expression of myosin heavy-chain IIB RNA in a subset of human masseter muscle fibres. Arch Oral Biol. 2001;46:1039–50.
Harrison BC, Allen DL, Leinwand LA. IIb or not IIb? Regulation of myosin heavy chain gene expression in mice and men. Skelet Muscle. 2011;1:1–9.
Stedman HH, Kozyak BW, Nelson A, Thesier DM, Su LT, Low DW, et al. Myosin gene mutation correlates with anatomical changes in the human lineage. Nature. 2004;428:415–8.
Schiaffino S, Reggiani C. Fiber types in mammalian skeletal muscles. Physiol Rev. 2011;91:1447–531.
Hoh JF. Superfast’ or masticatory myosin and the evolution of jaw-closing muscles of vertebrates. J Exp Biol. 2002;205:2203–10.
Resnicow D, Deacon JC, Warrick HM, Spudich JA, Leinwand LA. Functional diversity among a family of human skeletal muscle myosin motors. Proc Natl Acad Sci. 2010;107:1053–8.
Aksel T, Choe Yu E, Sutton S, Ruppel KM, Spudich JA. Ensemble force changes that result from human cardiac myosin mutations and a small-molecule effector. Cell Rep. 2015;11:910–20.
Harris DE, Work SS, Wright RK, Alpert NR, Warshaw DM. Smooth, cardiac and skeletal muscle myosin force and motion generation assessed by cross-bridge mechanical interactions in vitro. J Muscle Res Cell Motil. 1994;15:11–9.
Barany M. ATPase activity of myosin correlated with speed of muscle shortening. J Gen Physiol. 1967;50:197–218.
Walklate J, Geeves MA. Temperature manifold for a stopped-flow machine to allow measurements from −10 to +40°C. Anal Biochem. 2015;476:11–6.
Weiss A, McDonough D, Wertman B, Acakpo-Satchivi L, Montgomery K, Kucherlapati R, et al. Organization of human and mouse skeletal myosin heavy chain gene clusters is highly conserved. Proc Natl Acad Sci. 1999;96:2958–63.
Saez LJ, Gianola KM, McNally EM, Feghali R, Eddy R, Shows TB, et al. Human cardiac myosin heavy chain genes and their linkage in the genome. Nucleic Acids Res. 1987;15:5443–59.
Desjardins PR, Burkman JM, Shrager JB, Allmond LA, Stedman HH. Evolutionary implications of three novel members of the human sarcomeric myosin heavy chain gene family. Mol Biol Evol. 2002;19:375–93.
Toniolo L, Cancellara P, Maccatrozzo L, Patruno M, Mascarello F, Reggiani C. Masticatory myosin unveiled: first determination of contractile parameters of muscle fibers from carnivore jaw muscles. Am J Physiol Cell Physiol. 2008;295:C1535–42.
McGuigan K, Phillips PC, Postlethwait JH. Evolution of sarcomeric myosin heavy chain genes: evidence from fish. Mol Biol Evol. 2004;21:1042–56.
Ikeda D, Ono Y, Snell P, Edwards YJK, Elgar G, Watabe S. Divergent evolution of the myosin heavy chain gene family in fish and tetrapods: evidence from comparative genomic analysis. Physiol Genomics. 2007;32:1–15.
Rossi AC, Mammucari C, Argentini C, Reggiani C, Schiaffino S. Two novel/ancient myosins in mammalian skeletal muscles: MYH14/7b and MYH15 are expressed in extraocular muscles and muscle spindles. J Physiol. 2010;588:353–64.
Rowlerson A, Mascarello F, Veggetti A, Carpene E. The fibre-type composition of the first branchial arch muscles in Carnivora and Primates. J Muscle Res Cell Motil. 1983;4:443–72.
Rowlerson A, Pope B, Murray J, Whalen RB, Weeds AG. A novel myosin present in cat jaw-closing muscles. J Muscle Res Cell Motil. 1981;2:415–38.
Hoh JFY, Lim JHY, Kang LDH, Lucas CA. Expression of superfast myosin in the jaw-closing muscles of Crocodylus porosus. In: Grigg GC, Seebacher F, Franklin CE, editors. Crocodilian Biology and Evolution: Surrey Beatty & Sons; 2001. p. 156–64.
Schmidt-Rhaesa A. The evolution of organ systems. Oxford: Oxford University Press; 2007.
Mackie GO, Singla CL. Impulse propagation and contraction in the tunic of a compound ascidian. Biol Bull. 1987;173:188–204.
Schuchert P, Reber-Müller S, Schmid V. Life stage specific expression of a myosin heavy chain in the hydrozoan Podocoryne carnea. Differentiation. 1993;54:11–8.
Seipel K, Schmid V. Evolution of striated muscle: jellyfish and the origin of triploblasty. Dev Biol. 2005;282:14–26.
Steinmetz PRH, Kraus JEM, Larroux C, Hammel JU, Amon-Hassenzahl A, Houliston E, et al. Independent evolution of striated muscles in cnidarians and bilaterians. Nature. 2012;487:231–4.
Goodson H, Spudich J. Molecular evolution of the myosin family: relationships derived from comparisons of amino acid sequences. Proc Natl Acad Sci. 1993;90:659–63.
OOta S, Saitou N. Phylogenetic relationship of muscle tissues deduced from superimposition of gene trees. Mol Biol Evol. 1999;16:856–67.
Brunet T, Fischer AH, Steinmetz PR, Lauri A, Bertucci P, Arendt D. The evolutionary origin of bilaterian smooth and striated myocytes. Elife. 2016;5:1–24.
Sciote JJ, Horton MJ, Rowlerson AM, Link J. Specialized cranial muscles: how different are they from limb and abdominal muscles? Cells Tissues Organs. 2003;174:73–86.
Porter JD, Baker RS, Ragusa RJ, Brueckner JK. Extraocular muscles: basic and clinical aspects of structure and function. Surv Ophthalmol. 1995;39:451–84.
Spencer RF, Porter JD. Biological organization of the extraocular muscles. Prog Brain Res. 2006;151:43–80.
Purves D, Augustine G, Fitzpatrick D, Hall W, LaMantia A, White L. The actions and innervation of extraocular muscles. In: Neuroscience. 5th ed: Sinauer Associates, Inc; 2012. p. 436–40.
Porter JD. Extraocular muscle: cellular adaptations for a diverse functional repertoire. Ann N Y Acad Sci. 2002;956:7–16.
Wasicky R, Ziya-Ghazvini F, Blumer R, Lukas JR, Mayr R. Muscle fiber types of human extraocular muscles: a histochemical and immunohistochemical study. Investig Ophthalmol Vis Sci. 2000;41:980–90.
Jacoby J, Ko K, Weiss C, Rushbrook JI. Systematic variation in myosin expression along extraocular muscle fibres of the adult rat. J Muscle Res Cell Motil. 1990;11:24–40.
Briggs MM, Schachat F. The superfast extraocular myosin (MYH13) is localized to the innervation zone in both the global and orbital layers of rabbit extraocular muscle. J Exp Biol. 2002;205:3133–42.
Rubinstein NA, Porter JD, Hoh JFY. The development of longitudinal variation of myosin isoforms in the orbital fibers of extraocular muscles of rats. Investig Ophthalmol Vis Sci. 2004;45:3067–72.
Lim SJ, Jung HH, Cho YA. Postnatal development of myosin heavy chain isoforms in rat extraocular muscles. Mol Vis. 2006;12:243–50.
Zhou Y, Liu D, Kaminski HJ. Myosin heavy chain expression in mouse extraocular muscle: more complex than expected. Investig Ophthalmol Vis Sci. 2010;51:6355–63.
Park K, Lim J, Sohn S, Yeul S. Myosin heavy chain isoform expression in human extraocular muscles: longitudinal variation and patterns of expression in global and orbital layers. Muscle Nerve. 2012;45:713–20.
Wieczorek DF, Periasamy M, Butler-Browne GS, Whalen RG, Nadal-Ginard B. Co-expression of multiple myosin heavy chain genes, in addition to a tissue-specific one, in extraocular musculature. J Cell Biol. 1985;101:618–29.
Winters LM, Briggs MM, Schachat F. The human extraocular muscle myosin heavy chain gene (MYH13) maps to the cluster of fast and developmental myosin genes on chromosome 17. Genomics. 1998;54:188–9.
Briggs MM, Schachat F. Early specialization of the superfast myosin in extraocular and laryngeal muscles. J Exp Biol. 2000;203:2485–94.
Schachat F, Briggs MM. Phylogenetic implications of the superfast myosin in extraocular muscles. J Exp Biol. 2002;205:2189–201.
Kjellgren D, Thornell L, Andersen J, Pedrosa-Domellof F. Myosin heavy chain isoforms in human extraocular muscles. Investig Ophthalmol Vis Sci. 2003;44:1419–25.
Budak MT, Bogdanovich S, Wiesen MHJ, Lozynska O, Khurana TS, Rubinstein NA. Layer-specific differences of gene expression in extraocular muscles identified by laser-capture microscopy. Physiol Genomics. 2004;20:55–65.
Khanna S, Merriam AP, Gong B, Leahy P, Porter JD. Comprehensive expression profiling by muscle tissue class and identification of the molecular niche of extraocular muscle. FASEB J. 2003;17:1370–2.
Porter JD, Khanna S, Rao JS, Li J, Andrade FH. Extraocular muscle is defined by a fundamentally distinct gene expression profile. Proc Natl Acad Sci. 2001;98:12062–7.
Zhou Y, Cheng G, Dieter L, Hjalt TA, Andrade FH, Stahl JS, et al. An altered phenotype in a conditional knockout of Pitx2 in extraocular muscle. Investig Ophthalmol Vis Sci. 2009;50:4531–41.
Kitamura K, Miura H, Miyagawa-Tomita S, Yanazawa M, Katoh-Fukui Y, Suzuki R, et al. Mouse Pitx2 deficiency leads to anomalies of the ventral body wall, heart, extra- and periocular mesoderm and right pulmonary isomerism. Development. 1999;126:5749–58.
Diehl AG, Zareparsi S, Qian M, Khanna R, Angeles R, Gage PJ. Extraocular muscle morphogenesis and gene expression are regulated by Pitx2 gene dose. Investig Ophthalmol Vis Sci. 2006;47:1785–93.
Zhou Y, Liu D, Kaminski HJ. Pitx2 regulates myosin heavy chain isoform expression and multi-innervation in extraocular muscle. J Physiol. 2011;589:4601–14.
Hebert SL, Daniel ML, McLoon LK. The role of Pitx2 in maintaining the phenotype of myogenic precursor cells in the extraocular muscles. PLoS One. 2013;8
Conti-Fine BM, Milani M, Kaminski HJ. Myasthenia gravis: past, present, and future. J Clin Invest. 2006;116:2843–54.
Yu Wai Man CY, Chinnery PF, Griffiths PG. Extraocular muscles have fundamentally distinct properties that make them selectively vulnerable to certain disorders. Neuromuscul Disord. 2005;15:17–23.
Prabhakar BS, Bahn RS, Smith TJ. Current perspective on the pathogenesis of graves’ disease and ophthalmopathy. Endocr Rev. 2003;24:802–35.
Kaminski HJ, Mazen A, Leigh RJ, Katirji MB, Ruff RL. Extraocular muscles are spared in advanced Duchenne dystrophy. Ann Neurol. 1992;32:586–8.
Andrade FH, Porter JD, Kaminski HJ. Eye muscle sparing by the muscular dystrophies: lessons to be learned? Microsc Res Tech. 2000;48:192–203.
Porter JD, Merriam AP, Khanna S, Andrade FH, Richmonds CR, Leahy P, et al. Constitutive properties, not molecular adaptations, mediate extraocular muscle sparing in dystrophic mdx mice. FASEB J. 2003;17:893–5.
Kallestad KM, Hebert SL, Mcdonald AA, Daniel ML, Cu SR, Mcloon LK. Sparing of extraocular muscle in aging and muscular dystrophies: a myogenic precursor cell hypothesis. Exp Cell Res. 2011;317:873–85.
Verma M, Fitzpatrick KR, McLoon LK. Extraocular muscle repair and regeneration. Curr Ophthalmol Rep. 2017;5:207–15.
Hunt C. Mammalian muscle spindle: peripheral mechanisms. Physiol Rev. 1990;70:643–63.
Matthews PBC. Evolving views on the internal operation and functional role of the muscle spindle. J Physiol. 1981;320:1–30.
Hulliger M. The mammalian muscle spindle and its central control. Rev Physiol Biochem Pharmacol. 1984;101:1–110.
Tourtellotte WG, Milbrandt J. Sensory ataxia and muscle spindle agenesis in mice lacking the transcription factor Egr3. Nat Genet. 1998;20:87–91.
Blum KP, Lamotte D’Incamps B, Zytnicki D, Ting LH. Force encoding in muscle spindles during stretch of passive muscle. PLoS Comput Biol. 2017;13:1–24.
Banks R, Stacey M. Quantitative studies on mammalian muscle spindles. Boston: Springer; 1988.
Banks RW. An allometric analysis of the number of muscle spindles in mammalian skeletal muscles. J Anat. 2006;208:753–68.
Matthews PBC. Where anatomy led, physiology followed: a survey of our developing understanding of the muscle spindle, what it does and how it works. J Anat. 2015;227:104–14.
Ovalle WK. Fine structure of rat intrafusal muscle fibers. J Cell Biol. 1971;51:83–103.
Thornell LE, Carlsson L, Eriksson PO, Liu JX, Österlund C, Stål P, et al. Fibre typing of intrafusal fibres. J Anat. 2015;227:136–56.
Liu J, Eriksson P, Thornell L, Pedrosa-Domellof F. Fiber content and myosin heavy chain composition of muscle spindles in aged human biceps Brachii. J Histochem Cytochem. 2005;53:445–54.
Nichols T, Cope T. Cross-bridge mechanisms underlying the history-dependent properties of muscle spindles and stretch reflexes. Can J Phys Pharm. 2004;82:569–76.
Miller AJ. Craniomandibular muscles: their role in function and form: CRC Press; 1991.
Bredman JJ, Wessels A, Weijs WA, Korfage JA, Soffers CA, Moorman AF. Demonstration of “cardiac-specific” myosin heavy chain in masticatory muscles of human and rabbit. Histochem J. 1991;23:160–70.
Korfage JA, Van Eiiden TM. Myosin isoform composition of the human medial and lateral pterygoid muscles. J Dent Res. 2000;79:161–9.
Korfage JA, Koolstra JH, Langenbach GEJ, van Eijden TMGJ. Fiber-type composition of the human jaw muscles (part 2) role of hybrid fibers and factors responsible for inter-individual variation. J Dent Res. 2005;84:784–93.
English AW, Eason J, Pol M, Schwartz G, Shirley A. Different phenotypes among slow/beta myosin heavy chain-containing fibres of rabbit masseter muscle: a novel type of diversity in adult muscle. J Muscle Res Cell Motil. 1998;19:525–35.
Eason JM, Schwartz GA, Pavlath GK, English AW. Sexually dimorphic expression of myosin heavy chains in the adult mouse masseter. J Appl Physiol. 2000;89:251–8.
Boom HPW, van Spronsen PH, van Ginkel FC, van Schijndel RA, Castelijns JA, Tuinzing DB. A comparison of human jaw muscle cross-sectional area and volume in long- and short-face subjects, using MRI. Arch Oral Biol. 2008;53:273–81.
Kitai N, Kreiborg S, Bakke M, Paulsen HU, Møller E, Darvann TA, et al. Three-dimensional magnetic resonance image of the mandible and masticatory muscles in a case of juvenile chronic arthritis treated with the Herbst appliance. Angle Orthod. 2002;72:81–7.
Sciote JJ, Raoul G, Ferri J, Close J, Horton MJ, Rowlerson A. Masseter function and skeletal malocclusion. Rev Stomatol Chir Maxillofac Chir Orale. 2013;114:79–85.
Schiaffino S, Reggiani C. Molecular diversity of myofibrillar proteins: gene regulation and functional significance. Physiol Rev. 1996;76:371–423.
Ravosa MJ, Ning J, Costley DB, Daniel AN, Stock SR, Stack MS. Masticatory biomechanics and masseter fiber-type plasticity. J Musculoskelet Neuronal Interact. 2010;10:46–55.
Ciochon RL, Nisbett RA, Corruccini RS. Dietary consistency and craniofacial development related to masticatory function in minipigs. J Craniofac Genet Dev Biol. 1997;17:96–102.
Whiteley AT, Kendrick GS, Matthews JL. The effects of function on osseous and muscle tissues in the craniofacial area of the rat. Angle Orthod. 1966;36:13–7.
Kiliaridis S, Shyu BC. Isometric muscle tension generated by masseter stimulation after prolonged alteration of the consistency of the diet fed to growing rats. Arch Oral Biol. 1988;33:467–72.
Sato I, Konishi K. Effects of soft diet on rat masseter muscle mitochondrial development. Okajimas Folia Anat Jpn. 2004;81:67–74.
Langenbach G, van de Pavert S, Savalle W, Korfage H, van Eijden T. Influence of food consistency on the rabbit masseter muscle fibres. Eur J Oral Sci. 2003;111:81–4.
He T. Craniofacial morphology and growth in the ferret: effects from alteration of masticatory function. Swedish Dent J Suppl. 2004;165:1–72.
Kitagawa Y, Mitera K, Ogasawara T, Nojyo Y, Miyauchi K, Sano K. Alterations in enzyme histochemical characteristics of the masseter muscle caused by long-term soft diet in growing rabbits. Oral Dis. 2004;10:271–6.
Miehe B, Fanghänel J, Kubein-Meesenburg D, Nägerl H, Schwestka-Polly R. Masticatory musculature under altered occlusal relationships - a model study with experimental animals. Ann Anat. 1999;181:37–40.
Saito T, Ohnuki Y, Yamane A, Saeki Y. Effects of diet consistency on the myosin heavy chain mRNAs of rat masseter muscle during postnatal development. Arch Oral Biol. 2002;47:109–15.
Smith A, Weber CM, Newton J, Denny M. Developmental and age-related changes in reflexes of the human jaw-closing system. Electroencephalogr Clin Neurophysiol. 1991;81:118–28.
Kossioni AE, Karkazis HC. EMG study on the effect of ageing on the human masseteric jaw-jerk reflex. Gerodontology. 1994;11:30–8.
Newton JP, Yemm R, Abel RW, Menhinick S. Changes in human jaw muscles with age and dental state. Gerodontology. 1993;10:16–22.
Monemi M, Eriksson PO, Eriksson A, Thornell LE. Adverse changes in fibre type composition of the human masseter versus biceps brachii muscle during aging. J Neurol Sci. 1998;154:35–48.
Eriksson PO, Thornell LE. Histochemical and morphological muscle-fibre characteristics of the human masseter, the medial pterygoid and the temporal muscles. Arch Oral Biol. 1983;28:781–95.
Monemi M, Eriksson PO, Dubail I, Butler-Browne GS, Thornell LE. Fetal myosin heavy chain increases in human masseter muscle during aging. FEBS Lett. 1996;386:87–90.
Monemi M, Liu JX, Thornell L, Eriksson P. Myosin heavy chain composition of the human lateral pterygoid and digastric muscles in young adults and elderly. J Muscle Res Cell Motil. 2000;21:303–12.
Österlund C, Lindström M, Thornell LE, Eriksson PO. Remarkable heterogeneity in myosin heavy-chain composition of the human young masseter compared with young biceps brachii. Histochem Cell Biol. 2012;138:669–82.
Ricketts RM. A foundation for cephalometric communication. Am J Orthod. 1960;46:330–57.
Ingervall B, Helkimo E. Masticatory muscle force and facial morphology in man. Arch Oral Biol. 1978;23:203–6.
Kiliaridis S. Masticatory muscle influence on craniofacial growth. Acta Odontol Scand. 1995;53:196–202.
Weijs WA, Hillen B. Relationships between masticatory muscle cross-section and skull shape. J Dent Res. 1984;63:1154–7.
Kiliaridis S, Kälebo P. Masseter muscle thickness measured by ultrasonography and its relation to facial morphology. J Dent Res. 1991;70:1262–5.
Van Spronsen PH, Weijs WA, Valk J, Prahl-Andersen B, Van Ginkel FC. A comparison of jaw muscle cross-sections of long-face and normal adults. J Dent Res. 1992;71:1279–85.
Benington PC, Gardener JE, Hunt NP. Masseter muscle volume measured using ultrasonography and its relationship with facial morphology. Eur J Orthod. 1999;21:659–70.
Proffit WR, Fields HW, Nixon WL. Occlusal forces in normal- and long-face adults. J Dent Res. 1983;62:566–70.
Hunt NP, Cunningham SJ. The influence of orthognathic surgery on occlusal force in patients with vertical facial deformities. Int J Oral Maxillofac Surg. 1997;26:87–91.
Piancino MG, Isola G, Merlo A, Dalessandri D, Debernardi C, Bracco P. Chewing pattern and muscular activation in open bite patients. J Electromyogr Kinesiol. 2012;22:273–9.
Rowlerson A, Raoul G, Daniel Y, Close J, Maurage CA, Ferri J, et al. Fiber-type differences in masseter muscle associated with different facial morphologies. Am J Orthod Dentofac Orthop. 2005;127:37–46.
Boyd SB, Gonyea WJ, Finn RA, Woodard CE, Bell WH. Histochemical study of the masseter muscle in patients with vertical maxillary excess. J Oral Maxillofac Surg. 1984;42:75–83.
Brennan M, Antonyshyn O. The effects of temporalis muscle manipulation on skull growth: an experimental study. Plast Reconstr Surg. 1996;97:13–24.
Hohl TH. Masticatory muscle transposition in primates: effects on craniofacial growth. J Maxillofac Surg. 1983;11:149–56.
Zhu M, Zhao W, Jia L, Lu J, Qiao T, Qu Q. The oldest articulated osteichthyan reveals mosaic gnathostome characters. Nature. 2009;458:469–74.
Forey P, Janvier P. Agnathans and the origin of jawed vertebrates. Nature. 1993;361:129–34.
Brazeau MD, Friedman M. The origin and early phylogenetic history of jawed vertebrates. Nature. 2015;520:490–7.
Kang LHD, Rughani A, Walker ML, Bestak R, Hoh JFY. Expression of masticatory-specific isoforms of myosin heavy-chain, myosin-binding protein-C and tropomyosin in muscle fibers and satellite cell cultures of cat masticatory muscle. J Histochem Cytochem. 2010;58:623–34.
van Rooij E, Liu N, Olson EN. MicroRNAs flex their muscles. Trends Genet. 2008;24:159–66.
van Rooij E, Quiat D, Johnson BA, Sutherland LB, Qi X, Richardson JA, et al. A family of microRNAs encoded by myosin genes governs myosin expression and muscle performance. Dev Cell. 2009;17:662–73.
Nagase T, Kikuno R, Ishikawa K, Hirosawa M, Ohara O. Prediction of the coding sequences of unidentified human genes. IX. The complete sequences of 100 new cDNA clones from brain which can code for large proteins in vitro. DNA Res. 2000;7:143–50.
Bell ML, Buvoli M, Leinwand LA. Uncoupling of expression of an intronic MicroRNA and its myosin host gene by exon skipping. Mol Cell Biol. 2010;30:1937–45.
Warkman AS, Whitman SA, Miller MK, Garriock RJ, Schwach CM, Gregorio CC, et al. Developmental expression and cardiac transcriptional regulation of Myh7b, a third myosin heavy chain in the vertebrate heart. Cytoskeleton. 2012;69:324–35.
Rubio MD, Johnson R, Miller CA, Huganir RL, Rumbaugh G. Regulation of synapse structure and function by distinct myosin II motors. J Neurosci. 2011;31:1448–60.
Haraksingh RR, Jahanbani F, Rodriguez-Paris J, Gelernter J, Nadeau KC, Oghalai JS, et al. Exome sequencing and genome-wide copy number variant mapping reveal novel associations with sensorineural hereditary hearing loss. BMC Genomics. 2014;15:1–19.
Esposito T, Sampaolo S, Limongelli G, Varone A, Formicola D, Diodato D, et al. Digenic mutational inheritance of the integrin alpha 7 and the myosin heavy chain 7B genes causes congenital myopathy with left ventricular non-compact cardiomyopathy. Orphanet J Rare Dis. 2013;8:1–13.
Hanchard NA, Swaminathan S, Bucasas K, Furthner D, Fernbach S, Azamian MS, et al. A genome-wide association study of congenital cardiovascular left-sided lesions shows association with a locus on chromosome 20. Hum Mol Genet. 2016;25:2331–41.
Brown KM, MacGregor S, Montgomery GW, Craig DW, Zhao ZZ, Iyadurai K, et al. Common sequence variants on 20q11.22 confer melanoma susceptibility. Nat Genet. 2008;40:838–40.
Chatzinasiou F, Lill CM, Kypreou K, Stefanaki I, Nicolaou V, Spyrou G, et al. Comprehensive field synopsis and systematic meta-analyses of genetic association studies in cutaneous melanoma. J Natl Cancer Inst. 2011;103:1227–35.
Antonopoulou K, Stefanaki I, Lill CM, Chatzinasiou F, Kypreou KP, Karagianni F, et al. Updated field synopsis and systematic meta-analyses of genetic association studies in cutaneous melanoma: the MelGene database. J Invest Dermatol. 2015;135:1074–9.
Machida S, Noda S, Takao A, Nakazawa M, Matsuoka R. Expression of slow skeletal myosin heavy chain 2 gene in Purkinje fiber cells in chick heart. Biol Cell. 2002;94:389–99.
Nasipak BT, Kelley DB. The genome of the diploid anuran Xenopus tropicalis contains a novel array of sarcoplasmic myosin heavy chain genes expressed in larval muscle and larynx. Dev Genes Evol. 2008;218:389–97.
Garriock RJ, Meadows SM, Krieg PA. Developmental expression and comparative genomic analysis of Xenopus cardiac myosin heavy chain genes. Dev Dyn. 2005;233:1287–93.
Lonsdale J, Thomas J, Salvatore M, Phillips R, Lo E, Shad S, et al. The genotype-tissue expression (GTEx) project. Nat Genet. 2013;45:580–5.
Mele M, Ferreira P, Reverter F, DeLuca D, Monlong J, Sammeth M, et al. The human transcriptome across tissues and individuals. Science. 2015;348:660–5.
Miao X, Luo Q, Zhao H, Qin X. Co-expression analysis and identification of fecundity-related long non-coding RNAs in sheep ovaries. Sci Rep. 2016;6:1–10.
Miao X, Luo Q, Zhao H, Qin X. An integrated analysis of miRNAs and methylated genes encoding mRNAs and lncRNAs in sheep breeds with different fecundity. Front Physiol. 2017;8:1–14.
Hansel NN, Paré PD, Rafaels N, Sin DD, Sandford A, Daley D, et al. Genome-wide association study identification of novel loci associated with airway responsiveness in chronic obstructive pulmonary disease. Am J Respir Cell Mol Biol. 2015;53:226–34.
Bare LA, Morrison AC, Rowland CM, Shiffman D, Luke MM, Iakoubova OA, et al. Five common gene variants identify elevated genetic risk for coronary heart disease. Genet Med. 2007;9:682–9.
Luke MM, Lalouschek W, Rowland CM, Catanese JJ, Bolonick JI, Bui ND, et al. Polymorphisms associated with both noncardioembolic stroke and coronary heart disease: Vienna stroke registry. Cerebrovasc Dis. 2009;28:499–504.
Morrison AC, Bare LA, Chambless LE, Ellis SG, Malloy M, Kane JP, et al. Prediction of coronary heart disease risk using a genetic risk score: the atherosclerosis risk in communities study. Am J Epidemiol. 2007;166:28–35.
Yoshino S, Cilluffo R, Best P, Atkinson E, Aoki T, Cunningham J, et al. Single nucleotide polymorphisms associated with abnormal coronary microvascular function. Coron Artery Dis. 2014;25:281–9.
Neary MT, Neary JM, Lund GK, Holt TN, Garry FB, Mohun TJ, et al. Myosin heavy chain 15 is associated with bovine pulmonary arterial pressure. Pulm Circ. 2014;4:496–503.
Qin H, Hsu MKH, Morris BJ, Hoh JFY. A distinct subclass of mammalian striated myosins: structure and molecular evolution of “superfast” or masticatory myosin heavy chain. J Mol Evol. 2002;55:544–52.
Zhu J, Sanborn JZ, Diekhans M, Lowe CB, Pringle TH, Haussler D. Comparative genomics search for losses of long-established genes on the human lineage. PLoS Comput Biol. 2007;3:2498–509.
Kang L, Hughes S, Pettigrew J, Hoh J. Jaw-specific myosin heavy chain gene expression in sheep, dog, monkey, flying fox and microbat jaw-closing muscles. Basic Appl Myol. 1994;4:381–92.
Hoh JF, Kim Y, Sieber LG, Zhong WW, Lucas CA. Jaw-closing muscles of kangaroos express alpha-cardiac myosin heavy chain. J Muscle Res Cell Motil. 2000;21:673–80.
Reiser PJ, Bicer S, Chen Q, Zhu L, Quan N. Masticatory ('superfast’) myosin heavy chain and embryonic/atrial myosin light chain 1 in rodent jaw-closing muscles. J Exp Biol. 2009;212:2511–9.
Brunet M, Guy F, Pilbeam D, Mackaye HT, Likius A, Ahounta D, et al. A new hominid from the Upper Miocene of Chad, Central Africa. Nature. 2002;418:145–51.
Walker A, Leakey RE, editors. The Nariokotome Homo erectus skeleton. Cambridge: Harvard University Press; 1993.
Tobias PV, Gorge O. Volume 4: the skulls, endocasts, and teeth of homo habilis. Cambridge: University Press; 1991.
Kato C, Saeki Y, Yanagisawa K. Ca2+ sensitivities and transient tension responses to step-length stretches in feline mechanically-stripped single-fibre jaw-muscle preparations. Arch Oral Biol. 1985;30:429–32.
Perry GH, Verrelli BC, Stone AC. Comparative analyses reveal a complex history of molecular evolution for human MYH16. Mol Biol Evol. 2005;22:379–82.
Prüfer K, Racimo F, Patterson N, Jay F, Sankararaman S, Sawyer S, et al. The complete genome sequence of a Neanderthal from the Altai Mountains. Nature. 2014;505:43–9.
Perry GH, Kistler L, Kelaita MA, Sams AJ. Insights into hominin phenotypic and dietary evolution from ancient DNA sequence data. J Hum Evol. 2015;79:55–63.
McCollum MA, Sherwood CC, Vinyard CJ, Lovejoy CO, Schachat F. Of muscle-bound crania and human brain evolution: the story behind the MYH16 headlines. J Hum Evol. 2006;50:232–6.
Currie P. Human genetics: muscling in on hominid evolution. Nature. 2004;428:373–4.
Oh HJ, Choi D, Goh CJ, Hahn Y. Loss of gene function and evolution of human phenotypes. BMB Rep. 2015;48:373–9.
We would like to thank Dr. William K. Ovalle and Dr. Alberto C. Rossi for providing us with high resolution reprints and images for our re-use in compliance with the publishers permission. Thank you to Dr. Maureen Bjerke for the cartoon depiction of the myosin hetero-hexamer.
The authors are supported by the following funding sources: American Heart Association 18TPA34170194 (L.A. Leinwand and A. Karabina), National Institutes of Health R01GM029090 (L.A. Leinwand, L.A. Lee, and L.J. Broadwell), National Institutes of Health R01HL117138 (A. Karabina and L.A. Leinwand), National Institutes of Health T32GM008759 (L.A. Lee), and National Institutes of Health T32GN065103 (L.J. Broadwell).
Availability of data and materials
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.