- Open Access
Differential expression of myosin heavy chain isoforms in cardiac segments of gnathostome vertebrates and its evolutionary implications
Frontiers in Zoologyvolume 16, Article number: 18 (2019)
Immunohistochemical studies of hearts from the lesser spotted dogfish, Scyliorhinus canicula (Chondrichthyes) revealed that the pan-myosin heavy chain (pan-MyHC) antibody MF20 homogeneously labels all the myocardium, while the pan-MyHC antibody A4.1025 labels the myocardium of the inflow (sinus venosus and atrium) but not the outflow (ventricle and conus arteriosus) cardiac segments, as opposed to other vertebrates. We hypothesized that the conventional pattern of cardiac MyHC isoform distribution present in most vertebrates, i.e. MYH6 in the inflow and MYH7 in the outflow segments, has evolved from a primitive pattern that persists in Chondrichthyes. In order to test this hypothesis, we conducted protein detection techniques to identify the MyHC isoforms expressed in adult dogfish cardiac segments and to assess the pan-MyHC antibodies reactivity against the cardiac segments of representative species from different vertebrate groups.
Western and slot blot results confirmed the specificity of MF20 and A4.1025 for MyHC in dogfish and their differential reactivity against distinct myocardial segments. HPLC-ESI-MS/MS and ESI-Quadrupole-Orbitrap revealed abundance of MYH6 and MYH2 in the inflow and of MYH7 and MYH7B in the outflow segments. Immunoprecipitation showed higher affinity of A4.1025 for MYH2 and MYH6 than for MYH7 and almost no affinity for MYH7B. Immunohistochemistry showed that A4.1025 signals are restricted to the inflow myocardial segments of elasmobranchs, homogeneous in all myocardial segments of teleosts and acipenseriforms, and low in the ventricle of polypteriforms.
The cardiac inflow and outflow segments of the dogfish show predominance of fast- and slow-twitch MyHC isoforms respectively, what can be considered a synapomorphy of gnathostomes. The myocardium of the dogfish contains two isomyosins (MYH2 and MYH7B) not expressed in the adult heart of other vertebrates. We propose that these isomyosins lost their function in cardiac contraction during the evolution of gnathostomes, the later acquiring a regulatory role in myogenesis through its intronic miRNA. Loss of MYH2 and MYH7B expression in the heart possibly occurred before the origin of Osteichthyes, being the latter reacquired in polypteriforms. We raise the hypothesis that the slow tonic MYH7B facilitates the peristaltic contraction of the conus arteriosus of fish with a primitive cardiac anatomical design and of the vertebrate embryo.
Myosins form a large and diverse super-family of actin-dependent motor proteins responsible for contractile movement in cells and tissues . The myosin class II or conventional sarcomeric myosins constitute the group of myosins responsible for contraction in striated muscle, including the myocardium . They are hexameric proteins composed of two myosin heavy chains (MyHC) and four myosin light chains (MLC) . While MLC have regulatory functions, MyHC is responsible for the physical interaction with actin, producing the contraction of the sarcomere [20, 21].
Different types of muscle fibres contain distinct types of myosins, which are formed by specific MyHC isoforms or isomyosins . The sequences of cardiac MyHC isoforms, especially those corresponding to the motor domain including the ATP-binding site, appear to be very similar in all vertebrates . In the myocardium of adult vertebrates, two isomyosins are expressed across different cardiac segments. It has been shown that orthologs of human myh6 are predominantly expressed in the atrial myocardium of mammals (α-MyHC), birds (aMyHC), amphibians (α-MyHC) and teleosts (aMyHC), whereas orthologs of human myh7 are predominantly expressed in the ventricular myocardium of mammals (β-MyHC), birds (vMyHC) and teleosts (vMyHC) [4, 6, 15, 17, 34, 36, 38, 39, 49, 59, 60]. The identity and location of MyHC isoforms in adult hearts of chondrichthyans remain unknown.
Different pan-sarcomeric MyHC antibodies are routinely used to detect myosins in histological sections of myocardium, the most widely used being MF20 and A4.1025. The MF20 monoclonal antibody recognizes an epitope located in the tail domain of the MyHC . Immunohistochemistry with MF20 antibodies has been extensively used to highlight the myocardium of different groups of vertebrates, such as chondrichthyans , teleosts , Xenopus , birds  and mammals . The A4.1025 antibody is another pan-sarcomeric MyHC monoclonal antibody which reacts with an epitope near the ATP-binding site of the MyHC in a wide variety of species [7, 10]. Immunohistochemical studies with A4.1025 are less common in the literature than those with MF20, but myocardium-specific reactivity has been confirmed in a wide range of species, i.e. Xenopus [1, 48], mouse embryos , dogfish embryos  and zebrafish embryos .
In an ongoing research on the histomorphology of the chondrichthyan heart, using the lesser spotted dogfish (Scyliorhinus canicula) as an animal model, we have found that MF20 immunohistochemical signals show a homogeneous distribution throughout the heart, similar to that found in other vertebrates. Nevertheless, A4.1025 signals showed a heterogeneous pattern in the different cardiac segments. These preliminary results suggested that the dogfish may express distinct MyHC isoforms in a segment-specific manner. To test this hypothesis, we performed western blot analyses to demonstrate the specificity of MF20 and A4.1025 for the cardiac MyHC in the dogfish, slot blot combined with colorimetric analyses to verify segment-specific MyHC isoform expression and proteomic analyses using nano HPLC-ESI-MS/MS and ESI-Quadrupole-Orbitrap to identify and quantify the isomyosins differentially distributed in dogfish hearts. In addition, we calculated the relative affinity of the A4.1025 antibody for the cardiac MyHC isoforms by immunoprecipitation combined with protein quantification. Finally, we performed immunohistochemistry with MF20 and A4.1025 antibodies in the heart of different gnatosthome species (elasmobranchs, batoids, polypteriforms, acipenseriforms, teleosts and mammals) in order to assess the MyHC specific immunohistochemical signals in the cardiac segments of different taxa across the phylogeny.
Anatomical and immunohistochemical findings
The heart of the dogfish is composed of five segments contained in the pericardial cavity, namely, the myocardial sinus venosus, atrium, ventricle and conus arteriosus and the non-myocardial bulbus arteriosus (Fig. 1). We will refer to the inflow segments when describing the sinus venosus plus atrium, and to the outflow segments when describing the ventricle plus the conus arteriosus. Immunohistochemical staining with the antibody MF20 revealed a homogeneous signal in all myocardial segments (Fig. 1a). However, immunohistochemistry in consecutive sections with the antibody A4.1025 showed signals restricted to the myocardium of the inflow segments (Fig. 1b). A less intense mark was detected in the ventricular apex.
Both MF20 and A4.1025 signals were specifically restricted to myocardial tissue. This is evident when comparing the strong labelling of the atrial wall with the absence of labelling of the conal and atrioventricular valves and the wall of the bulbus arteriosus (Fig. 1). In addition, the inflow segment restricted A4.1025 staining was evident when observing the limit between the atrial and ventricular walls.
To find out whether the distinct A4.1025 immunohistochemical pattern detected in the dogfish occurs in other vertebrates, we examined hearts from elasmobranchs (gulper shark and sand ray), actinopterygians (gray bichir, Adriatic sturgeon, silver arowana, zebrafish) and mammals (Syrian hamster). In all these species, MF20 staining was strong and homogeneous in the wall of all myocardial segments (not shown). Immunohistochemistry with A4.1025 in the hearts from sturgeons, arowanas, zebrafish and hamsters also revealed a strong and homogeneous signal in the whole myocardium (Fig. 2). However, the hearts of the gulper shark (Fig. 2a), sand ray (Fig. 2b), dogfish (Fig. 2c) and gray bichir (Fig. 2d) showed heterogeneous immunoreactivity against A4.1025 antibodies. In general, the staining pattern in the heart of the gulper shark (Fig. 2a) and sand ray (Fig. 2b) was similar to that detected in the dogfish, with intense signals in the inflow region and low or absent signals in the ventricle and conus (Fig. 2c). The pattern of A4.1025 staining in the heart of the gray bichir was different to that found in the other species examined. In this species, A4.1025 signal was strong in the sinus venosus, atrium and conus arteriosus, but faint in the ventricle (Fig. 2d).
Slot and western blotting
The fig. 3 shows western blots of MyHC extractions from different cardiac regions of hamsters and dogfish. In both species, a clear single band with a molecular weight between 200 kDa and 250 kDa was detected with both antibodies for all the cardiac segments, except for the bulbus arteriosus, which displayed a faint signal. In extracts from hamster hearts, the four bands obtained using both antibodies showed a similar strong intensity (Fig. 3a, b). In extracts from the dogfish hearts, the bands revealed with the MF20 antibody (Fig. 3a) corresponding to the atrium, ventricle and conus arteriosus showed a similar strong intensity. By contrast, those revealed with the A4.1025 antibody (Fig. 3b) corresponding to ventricle and conus arteriosus were markedly less intense.
To quantify the intensity of the signals detected from the protein extractions of the dogfish, we performed colorimetric analyses in bands obtained by means of slot blot assays (Fig. 3c). For these blots, we used a constant amount of MyHC proteins from the same extractions as in western blots. Atrial A4.1025 signals were almost three-fold more intense than ventricular bands (p=0.00001), whereas the intensity of MF20 atrial and ventricular signals showed no statistical differences (p=0.74183).
Proteomic analysis by nano HPLC-ESI-MS/MS
The bands at 200-250 kDa corresponding to the myocardial inflow and outflow segments of dogfish samples were excised separately and subjected to peptide sequence analysis by nano HPLC-ESI-MS/MS. The resulting sequences were identified by MASCOT searching in a Callorrinchus milii (ghost or elephant shark) and in a chordate databases. A total of 275 sequence matches were found in the chordate database, 256 (93.03%) of which corresponded to isomyosins. In the C. milii database a total of 129 sequence matches were found, all of them corresponding to isomyosins.
From the inflow segments sample, the sequence matches detected in the C. milii database corresponded to four myosin isoforms: two cardiac isoforms (MYH6, MYH7), a skeletal isoform (MYH2) and a non-muscle isoform (MYH10) (Table 1). In the chordate database, the sequence matches corresponded to eight myosin isoforms: three cardiac isoforms (MYH6, MYH7, MYSC), four skeletal isoforms (MYH1B, MYH2, MYH8, MYH13) and a non-muscle isoform (MYH9) of a single species (Table 1).
From the outflow segments sample, the sequence matches detected in the C. milii database corresponded to six myosin isoforms: four cardiac isoforms (MYH6, MYH7, MYH7-like, MYH7B), a skeletal isoform (MYH2) and a non-muscle isoform (MYH10) (Table 1). In the chordate database, the sequence matches corresponded to six myosin isoforms: two cardiac isoforms (MYH6, MYH7B) of two species and four skeletal isoforms (MYH2, MYH4, MYH8; MYHemb) of four species (Table 1). All individual ions scores were >20 to indicate identity or extensive homology (p<0.01).
Myosin heavy chain isoforms quantification by ESI-Quadrupole-Orbitrap
The bands at 200-250 kDa corresponding to the myocardial inflow and outflow segments of two additional dogfish samples were excised separately and subjected to peptide quantification by ESI-Quadrupole-Orbitrap. The resulting sequences were identified by SEQUEST searching in the same C. milii database.
The sequence matches detected in the C. milii database corresponded to the same isomyosins detected by nano HPLC-ESI-MS/MS (Table 2; MYH2, MYH6, MYH7 for the inflow segments and MYH2, MYH6, MYH7, MYH7B for the outflow segments). Traces of MYH9, MYH10 and MYH11 were additionally found in both segments, what can be explained by the higher sensitivity of this technique. All individual ions scores were >20 to indicate identity or extensive homology (p<0.01).
The table 2 shows the average total and percental abundance of each MyHC isoform in the inflow and outflow segments. In the inflow segments MYH6 (42.7%) and MYH2 (34.3%) were the major isoforms, followed by MYH7 (22.0%), whereas MYH7B was almost absent. The non-sarcomeric isoforms MYH9, MYH10 and MYH11 represented less than 1% of the total myosins. In the outflow segments MYH7 was the major isoform (45.3%), followed by MYH7B (24.0%), MYH6 (20.4%) and MYH2 (8.3%). Non-sarcomeric isoforms (MY9, MYH10 and MYH11) represented around 2% of the total myosins.
A4.1025 immunoprecipitation and isoform antibody affinity analysis
One of the myosin extractions of the dogfish outflow used for ESI-Quadrupole-Orbitrap was immunoprecipitated using an affinity column coated with A4.1025. The bands at 200-250 kDa corresponding to the initial extraction and the same extraction after column elution were excised separately and subjected to peptide quantification by ESI-Quadrupole-Orbitrap. The resulting sequence were identified by SEQUEST searching in the C. milii database.
The cardiac MyHC sequence matches detected in the C. milii database from both bands corresponded to the same isomyosins detected by nano HPLC-ESI-MS/MS and ESI-Quadrupole-Orbitrap in previous experiments (compare Table 1 with Tables 2 and 3). The table 3 shows the total and percentual abundance of each cardiac MyHC isoform in both the initial extraction and the eluate. The isomyosin abundances of the initial extraction showed similar values to those obtained in the previous experiments (MYH7=49.1%; MYH7B=28.3%; MYH6=15.4%; MYH2=7.2%). However, the isomyosin abundances after A4.1025 column elution changed significantly. While MYH7 showed a relative abundance similar to that in the original extract (46.7%), MYH2 and MYH6 were overrepresented (15.3% and 31.4%, respectively) and MYH7B was reduced almost four folds (6.5%). The relative affinity of the A4.1025 for the different cardiac MYH isoforms in the dogfish outflow segments was estimated by dividing the relative amount of each isoform in the eluate by the relative amount in the initial extraction. MYH2 and MYH6 showed more than double affinity for A4.1025 than MYH7, and almost 10 times more affinity than MYH7B.
In most vertebrates studied until now, i.e. mammals, birds and teleosts, two MyHC isoforms, MYH6 and MYH7, are expressed in the adult heart. Orthologs of MYH7 are predominantly expressed in the ventricle, whereas orthologs of MYH6 are predominantly expressed in the atrium and, in fish species, also in the sinus venosus. This data was obtained by means of in situ hybridization [4, 6, 49, 59, 60] and electrophoretic separation [34, 39]. One exception is Xenopus, in which, as Garriock et al.  described after sequence comparison and in situ hybridization, an ortholog of MYH15 is the main ventricular MyHC isoform. MYH15 has become a pseudogene in mammals . However, the isomyosin distribution in the heart of chondrichthyans, the earliest phylogenetical group of living gnatosthomes, remains unknown. Knowledge of MYH isoform expression in the chondrichthyan heart would inform about (1) what was the primitive condition of the trait in gnatosthomes, and (2) when the cardiac MYH distribution shared by most gnatosthomes did appear.
In an ongoing immunohistochemical study of the chondrichthyan heart using the dogfish as an animal model, we detected a striking pattern of immunohistochemical signals with the two anti-pan-MyHC antibodies MF20 and A4.1025. While MF20 signals were homogeneously distributed in the four myocardial segments of the chondrichthyan heart, namely, the sinus venosus, atrium, ventricle, and conus arteriosus, A4.1025 signals were strong in the sinus venosus and atrium, but faint in the conus arteriosus and the ventricle (Fig. 1). This segment-specific immunohistochemical pattern differs from the homogeneous pattern detected with A4.1025 in all the vertebrate species examined, i.e. hamsters (present results), mice [29, 52], Xenopus [1, 44, 48] and zebrafish [22, 24, 44]. This data already suggests that the cardiac segments of the dogfish heart express distinct MyHC isoforms that show differential affinity for MF20 and A4.1025 antibodies.
In order to corroborate our immunohistochemical results, we performed western blot analyses using MF20 and A4.1025 antibodies. Single 220 kDa bands, the approximate molecular weight of MyHC , were obtained after applying a myosin-specific extraction method in dogfish and hamster hearts. This confirms the specificity of both antibodies for MyHC in both species. In addition, the different intensity between the A4.1025 bands corresponding to the inflow and outflow myocardial segments of the dogfish heart strengthens the immunohistochemical findings. Western blot results were corroborated by colorimetric quantification of bands obtained in slot blot analyses. In the dogfish, the intensity of A4.1025 signals in the outflow segments was almost three-fold significantly lower than that in the inflow segments, whereas the intensity of MF20 signals did not significantly differ between blots.
The wall of the bulbus arteriosus does not contain myocardium. It is mainly composed of smooth muscle cells, elastin and collagen [13, 30]. However, our western blot analyses of bulbus extracts showed a narrow band using MF20 or A4.1025 antibodies (Fig. 3). In this regard, we cannot exclude that myocardial tissue of the distal portion of the conus arteriosus was removed together with the bulbus during dissection for protein extraction.
MF20 and A4.1025 antibodies recognize different epitopes of the MyHC protein. MF20 is specific for the tail region , whereas the epitope recognized by A4.1025 is in the head region [7, 10]. Thus, our immunohistochemical, western blot and slot blot results strongly suggest that in contrast to the mammalian heart, in the dogfish heart there is at least one MyHC isoform with a different affinity for MF20 and A4.1025 antibodies. In addition, it can be deduced that this or these isoforms are differentially expressed in the inflow versus the outflow myocardial segments in this species.
To test our hypothesis and identify the MyHC isoform/s differentially expressed in the dogfish cardiac segments, we performed HPLC-ESI-MS/MS analyses in the samples of the inflow and outflow segments. The genomic DNA of the dogfish has not been fully sequenced, and there is no record of the peptide sequence of the MyHC isoforms in this species. Although the whole genome of three additional elasmobranch species (Chiloscyllium punctatum, Scyliorhinus torazame and Rhincodon typus) has been recently sequenced , these sequences are unassembled, and hence they are not accessible to the MS search engine. Therefore, we compared the trypsinised peptides of the inflow and the outflow tract myocardium of the dogfish with protein sequences obtained from the genomic database of the elephant shark (Callorhinchus milii) provided by Venkatesh et al. . C. milii, which belongs to the holocephalans, the sister group of the elasmobranchs, is the only chondrichthyan with the MYH protein sequences included in the UniProtKB/Swiss-Prot database. The MS search engine used in this study (MASCOT and SEQUEST) acquires its maximal efficiency when the highly reliable UniProtKB/Swiss-Prot database is used for comparisons. Samples from both the inflow and outflow segments of the dogfish contained peptides corresponding to MYH10, MYH2, MYH6 and MYH7, whereas peptides corresponding to MYH7B were only detected in the outflow segments (Table 1). In order to confirm these results, we used a chordate peptide database. As it might be expected taking into account the high number of entries in the chordate database and the high level of conservation of MyHC isoforms, several additional sequences matched with peptides in our samples (MYH1B, MYH4, MYH8, MYH9, MYH13, MYHemb and MYSC). Despite this, MYH7B was again found to match with peptides from the outflow segments but not with those from the inflow segments (Table 1).
In order to confirm the distribution of MyHC isoforms revealed by HPLC-ESI-MS/MS and to quantify the relative abundance of these isoforms, we performed ESI-Quadrupole-Orbitrap assays in additional dogfish samples. The results confirmed the data obtained by the previous experiments. The dogfish cardiac outflow segments contain four main sarcomeric isomyosins: MYH2, MYH6, MYH7 and MYH7B, whereas the inflow segments contain only MYH2, MYH6 and MYH7. In addition, the relative abundance of these isoforms is segment-specific. The fast-twitch MYH2 and MYH6 account for more than 75% of the total MyHC isoforms in the inflow segments, where MYH7 (22%) is the only slow-twitch isoform. By contrast, the two slow-twitch isoforms MYH7 and MYH7B are predominant in the outflow segments (70%), whereas MYH2 and MYH6 account for less than 30% of total isomyosins in these segments.
According to the contraction speed, there are two types of skeletal isomyosins: slow and fast-twitch isoforms. In all the vertebrates studied until now, the fast-twitch MyHC that characterizes the inflow segments is MYH6, whereas the slow-twitch MYH7 predominates in the ventricle . Our results show that the myocardium of chondrichthyans contains two additional MyHC isoforms, which are not expressed in the heart of other vertebrates, MYH2 (fast-twitch), abundant in the inflow segments, and MYH7B (slow-twitch), concentrated in the outflow segments. Nevertheless, the myocardium of chondrichthyans shows the same segment-specific composition of isomyosins as other vertebrates with respect to their contraction speed: the inflow segments contains predominantly fast-twitch isomyosins (MYH2 and MYH6; 77%), while the outflow segments contains predominantly slow-twitch isomyosins (MYH7 and MYH7B; 70%). Thus, this segment-specific distribution of fast vs. slow MyHC isoforms can be considered a shared trait of Gnathostomes. New studies in cyclostomes may help to elucidate whether this trait is a synapomorphy or a symplesiomorphy for Gnathostomes. However, the composition of specific cardiac isomyosins has changed during the evolution of the group, reducing its variability in derived taxa.
MYH2 is a fast-twitch skeletal muscle MyHC (), which has not been detected in the myocardium of any adult vertebrate studied until now. Our results point to two possible evolutionary scenarios. In early gnathostomes, MYH2 may have been a cardiac myosin isoform that lost its function in cardiac contraction throughout evolution. Given that MYH2 is absent in the heart of both mammals and teleosts, loss of cardiac expression of this isomyosin might have taken place before the origin of Osteichthyes. Alternatively, MYH2 cardiac expression may have been absent in early gnathostomes, being an independent acquisition in chondrichthyans. MYH7B (firstly known as MYH14) is a poorly understood slow tonic MyHC isoform which, together with MYH15 and MYH16 is considered an ancient MyHC isoform that could have given rise to the current main cardiac isomyosins, MYH6 and MYH7, after a duplication process [12, 53]. Although MYH7B transcripts have been detected in the heart of multiple vertebrate species, including human, mouse , rat , chick, Xenopus  and zebrafish , the cardiac expression of the protein has been the subject of debate. Warkman et al.  self-dismissed a previous positive result in the adult human heart due to a failure in their early experimental design. In addition, it was not clear if an ortholog of MYH7B (ssMHC) is really expressed in the cardiac conduction system of the chick , because later experiments were not able to reproduce this result . In mammals, MYH7B protein has only been found in the extraocular musculature of humans and rats , and in the soleus muscle of mice . Absence of MYH7B expression in the cardiac muscle relies on the presence of a premature termination codon in exon 7 of the gene. The activation of a splicing machinery in skeletal muscle that excludes this exon allows MYH7B translation in this tissue . Therefore, our results clearly show for the first time the expression of MYH2 and MYH7B protein in the adult myocardium of an extant vertebrate, namely, the lesser spotted dogfish.
Similar to myh6 and myh7, myh7b DNA sequence includes an intronic miRNA (MyomiR) named miR-499 , which is present in the myh7b gene of most vertebrate species . The importance of miR-499 relies on its function as a post-transcriptional regulation element of myh7b . It has been also proposed that miR-499 is able to interact with Sox6 to regulate the pathway responsible for the differentiation of the slow and the fast musculature . We propose that in teleosts, birds and mammals MYH6 and MYH7 became mainly restricted to the inflow and outflow segments respectively, whereas MYH2 and MYH7B disappeared from the heart. The mRNA of MYH7B probably acquired a regulatory role in cardiac and extracardiac myogenesis through its intronic miRNA, miR-499 [5, 28, 35, 55].
Using a combination of Orbitrap-Quadrupole and A4.1025 immunoprecipitation methodologies, we were able to calculate the relative affinity of the MyHC isoforms of the dogfish heart for the A4.1025 antibody. The slow-twitch MyHC isoforms showed a reduced affinity compared to the fast-twitch isoforms, being MYH7B particularly unable to bind the antibody (Table 3). These data allow interpretation of the A4.1025 immunohistochemical results in terms of MyHC isoform content. The strikingly low A4.1025 signals in the cardiac outflow segments of the dogfish (Fig. 1) is most probably due to three factors: 1) the high proportion of slow-twitch MyHC isoforms in these segments; 2) the low affinity of these isoforms (particularly MYH7B) for A4.1025; 3) the low concentration of total MyHC in the outflow compared with the inflow segments (see ).
If we assume that the A4.1025 antibody has similar affinities for MyHC isoforms of different gnathostomes, given the high level of conservation of the ATP-binding site in the MyHC protein sequences among vertebrates, our initial immunohistochemical study of different vertebrate species (Fig. 2) can be reinterpreted in terms of specific isomyosin abundance. The sand ray and the gulper shark showed the same A4.1025 immunohistochemical pattern as the dogfish, suggesting a distinctly shared MyHC distribution in the heart of chondrichthyans (MYH2 and MYH6 in the inflow and MYH7 and MYH7B in the outflow segments). The silver arowana, the sturgeon, the zebrafish and the hamster showed the expected A4.1025 homogeneous immunohistochemical pattern in all the myocardial segments, suggesting that the conventional MyHC distribution (MYH6 in the inflow and MYH7 in the outflow segments) appeared at the base of Osteichthyes. By contrast, the gray bichir, a representative of polypteriforms, the earliest diverged group of extant actinopterygians [27, 51], showed a different A4.1025 immunohistochemical pattern of staining compared to all the other species studied. A4.1025 signals were intense in all the myocardial segments except for the ventricle. This can be interpreted as a relatively high abundance of slow-twitch isomyosins, particularly MYH7B, only in the ventricle but not in the conus arteriosus. From the evolutionary viewpoint, this finding implies a reacquisition of MYH7B expression in the heart of polypteriforms (Fig. 4). In addition, it will be interesting to investigate possible differences in the type of myocardial contraction between bichirs and other actinopterygians. Further studies are required to test the hypothesis depicted in fig. 4, by analyzing the MyHC isoform distribution in the heart of additional vertebrate species, particularly in sarcopterygians and other early diverged actinopterygians (as lepisosteiforms) would be of high interest.
In gnathostomes there are two basic cardiac anatomical designs, primitive and modern with respect to the type of cardiac contraction [18, 30]. In mammals, birds and most teleosts, which show the modern design, the main contractile segments are the atrium and the ventricle, whereas in chondrichthyans and early actinopterygians, which show the primitive design, the conus arteriosus is interposed between the ventricle and the bulbus arteriosus. While atrium and ventricle work as synchronous, fast-contracting segments, the conus arteriosus is a slow-contracting segment with peristaltoid contraction. It has been hypothesized that the conus arteriosus serves as an accessory pumping chamber , as an elastic reservoir to minimize pressure fluctuations [9, 26], and to allow a correct function of the conal valvular apparatus  in elasmobranchs. It has been also proposed that the myocardium of the conus arteriosus of teleosts plays an active function in the performance of the outflow tract valves . Birds and mammals also show the primitive cardiac anatomical design at early embryonic stages [11, 16, 37]. In these embryos, MYH isoform coexpression in the ventricle and the conus arteriosus, including MYH7B at least in mammals , allows a slow conduction velocity and a peristaltoid contraction . After cardiac septation, a shift to segment-specific MYH single isoform expression allows development of the fast conduction velocity and synchronous contraction of atrium and ventricle, necessary for the adult type of one-way valves that start to develop at the end of septation (Moorman & Lammers, 1994). According to our results, MYH7B is expressed in the ventricle and/or the conus arteriosus of Chondrichthyes and Polypteriforms, which are fish with a primitive cardiac anatomical design, i.e. with a well-developed conus arteriosus. It can then be hypothesized that the slow tonic MYH7B is required for the peristaltic contraction of the conus arteriosus in fish with a primitive cardiac anatomical design. MYH7B peptides would have disappeared from the adult heart once the modern design, including a synchronous and fast contraction, has developed in tetrapods and teleosts. The MyHC isoform expression pattern in polypteriforms and acipenseriforms, both with a long conus but a different MYH7B cardiac content compared to chondrichthyans, may represent transitional conditions between the primitive and the modern cardiac anatomical design. In this context, it would be interesting to further investigate the anatomy, type of contraction and precise MYH isoform composition in these two actinopterygian groups.
Similar to the other gnathostomes, the cardiac inflow segments of the dogfish show predominance of fast-twitch MyHC isoforms whereas the outflow segments show predominance of slow-twitch isoforms. Therefore, the segment-specific distribution of fast- vs.slow-twitch isomyosins in the heart can be considered a synapomorphy or a symplesiomorphy for gnathostomes, pending the analysis of cyclostomes.
Besides the two conventional cardiac MyHC isoforms (MYH6 and MYH7), the myocardium of the dogfish contains two additional isomyosins (MYH2 and MYH7B), not expressed in the heart of other vertebrates. These two isoforms may have lost their function in cardiac contraction during the evolution of gnathostomes, the latter acquiring a regulatory role in cardiac and extracardiac myogenesis through its intronic miRNA. Our results suggest that this change in the distribution of cardiac MyHC isoforms occurred before the origin of Osteichthyes, and that polypteriforms reacquired ventricular MYH7B expression.
We raise the hypothesis that the slow tonic MYH7B may facilitate the peristaltic contraction of the conus arteriosus of fish with a primitive cardiac anatomical design and of the vertebrate embryo.
The sample examined consisted of 31 adult hearts of the following eight gnathostome species: Leucoraja circularis (sandy ray; n=2) and Centrophorus granulosus (gulper shark; n=2) both collected by scientific vessels of the Spanish Institute of Oceanography fleet; Scyliorhinus canicula (lesser spotted dogfish; n=15), Polypterus senegalus (gray bichir, n=2), Acipenser naccarii (Adriatic sturgeon, n=2), Osteoglossum bicirrhossum (silver arowana, n=2), and Danio rerio (zebrafish, n=2) from local suppliers, and Mesocricetus auratus (golden hamster, n=4) from an outbred strain commercialized by Janvier (RjHan: AURA; France).
The animals were handled in accordance with the European and Spanish guidelines for animal welfare. The animals processed in the laboratory were killed by anaesthetic overdose with 0.04% MS-222 (Sigma-Aldrich Chemical Co., Poole, UK) in freshwater for the actinopterygians, or by CO2 inhalation for hamsters. The ventral region of the pericardial cavity was exposed by means of a longitudinal incision along the anterior, midventral line of the animal, and the heart was removed and rinsed in elasmobranch buffer (16.38 g L-1 NaCl, 0.89 g L-1 KCl, 1.11 g L-1 CaCl2, 0.38 g L-1 NaHCO3, 0.06 g L-1 NaHPO4, 21.6 g L-1 urea, pH 7.2) for the chondrichthyans or 0.1M phosphate buffer saline (PBS; pH 7.3) for the actinopterygians and mammals. For the proteomic analyses of samples from S. canicula and M. auratus, inflow and outflow segments were dissected through the atrioventricular and conoventricular limits or through the atrioventricular and ventriculoarterial limits, respectively.
Five hearts from S. canicula and two from each of the other species used in this work were fixed in methanol/acetone/water (2:2:1) or in 4% paraformaldehyde diluted in elasmobranch buffer for chondrichthyans or in PBS for actinopterygians and mammals. Then, the hearts were dehydrated in a graded series of ethanol and embedded in Histosec paraffin (Merck KGaA; Darmstadt, Germany). Serial sagittal sections of hearts cut at 10 μm were dewaxed, hydrated and washed in Tris-PBS (TPBS, pH 7.8). For immunoperoxidase, endogenous peroxidase activity was quenched by incubation for 30 min with 3% hydrogen peroxide in TPBS. After washing with TPBS, non-specific binding sites were saturated for 30 min with 10% sheep serum, 1% bovine serum albumin (BSA) and 0.1% Triton X-100 in TPBS (SBT). The slides were then incubated overnight at 4°C in the primary antibody diluted in SBT. Negative control slides were incubated in SBT only. After incubation, the slides were washed in TPBS, incubated for 1.5 h at room temperature in biotin-conjugatedanti-mouse IgG (Sigma-Aldrich Chemical Co., Poole, UK) diluted 1:1000 in SBT, washed again and incubated for 1 h in ExtraAvidin®-peroxidase complex (Sigma-Aldrich Chemical Co., Poole, UK) diluted 1:250 in TPBS. Peroxidase activity was detected using Sigma Fast 3,3’-diaminobenzidine tablets (Sigma-Aldrich Chemical Co., Poole, UK) following the indications of the supplier. Some sections were counterstained with hematoxylin staining. The sections were observed with a Leica DMSL light microscope or with an Olympus VS120 virtual microscopy slide scanning system (Olympus, Tokyo, Japan) equipped with the VS-ASW software (Olympus, Tokyo, Japan) and viewed using the free of charge software OlyVIA (Informer Technologies, Inc., Walnut, CA, USA).
In selected sections, ExtraAvidin®-peroxidase and further incubations were replaced by ExtrAvidin-Cy3™ (Sigma-Aldrich Chemical Co., Poole, UK) incubation and counterstaining with 4’,6-diamidino-2-phenylindole dihydroclhoride (DAPI) (Sigma-Aldrich Chemical Co., Poole, UK). They were observed and photographed with an Olympus VS120 virtual microscopy slide scanning system (Olympus, Tokyo, Japan).
The primary antibodies were MF20 (DSHB Cat# MF-20, RRID: AB_2147781) and A4.1025 (DSHB Cat# A4.1025, RRID: AB_528356) both anti-MyHC monoclonal antibodies (Developmental Studies Hybridoma Bank, University of Iowa). They were used at dilutions of 1:20 and 1:200 of the respective supernatants.
Western blot, slot blot and colorimetric analysis
For western and slot blot analyses, MyHC was extracted and quantified from tissue lysates following a specific protocol for MyHC . Protein concentrations were evaluated with a NanoDrop ND-1000 spectrophotometer (Nanodrop technologies, Wilmintogton, DE, USA). A total of five adult hearts from S. canicula and two from M. auratus were used.
For western blot analysis, 1 μg of protein was diluted in 10 μL of sample buffer (Tris 0.312 M pH 6.8 with 5% SDS, 0.25% bromophenol blue, 50% glycerol and 2-mercaptoethanol 0.9 M), heated at 100°C for 5 min and resolved by SDS-PAGE in an 8% acrylamide gel. Electrophoresis was performed at 125 V for 2 h. Proteins were transferred onto nitrocellulose membranes using a Trans-Blot SD semi-dry blotting system (BioRad Laboratories, California, USA) run at 15V for 1.5 h. After washing with TPBS, non-specific binding sites were saturated for 30 min with 5% milk in TPBS. The membranes were then incubated overnight at 4°C in the primary antibody diluted in SBT. MF20 and A4.1025 were used at 1:1000 and 1:2000 dilution of the supernatants, respectively. After incubation, the membranes were washed in TPBS, incubated for 1.5 h at room temperature in biotin-conjugatedanti-mouse IgG (Sigma-Aldrich Chemical Co., Poole, UK) diluted 1:1000 in SBT, washed again and incubated for 1 h in avidin-peroxidase complex (Sigma-Aldrich Chemical Co., Poole, UK) diluted 1:250 in TPBS. Peroxidase activity was detected using the same protocol as for immunohistochemistry. The bands obtained were photographed with a Gel Doc™ XR+ System (BioRad Laboratories, California, USA) equipped with the Image Lab software (BioRad Laboratories, California, USA).
For slot blot analysis, nitrocellulose membranes were pre-soaked in PBS for five min. Fifteen μL of myosin extraction (3 μg/μL) were added to the blotter (Hoefer Scientific Instrument, San Francisco, USA) applying vacuum. The membranes were removed and washed in TPBS. Non-specific binding sites were saturated for 30 min with 5% milk in TPBS. The membranes were then incubated overnight at 4°C in the primary antibody diluted in SBT. MF20 and A4.1025 were used at 1:1000 and 1:2000 dilution of the supernatants, respectively. After incubation, the membranes were washed in TPBS, incubated for 1.5 h at room temperature in biotin-conjugatedanti-mouse IgG (Sigma-Aldrich Chemical Co., Poole, UK) diluted 1:1000 in BSA, washed again and incubated for 1 h in avidin-peroxidase complex (Sigma-Aldrich Chemical Co., Poole, UK) diluted 1:250 in TPBS. Peroxidase activity was detected using the same protocol as for immunohistochemistry. The membranes were photographed in black and white. The intensity of the signal in each band was measured using the software FIJI . Average data from five specimens were obtained and statistical differences were tested using t-student test.
Sample preparation for proteomic analysis
All proteomic analyses were performed at the Proteomics Core Facility of the University of Málaga. The bands identified on the western blot were isolated from the gel, destained, reduced with DTT, alkylated with iodoacetamide and digested in-gel with trypsin (Promega, Southampton, UK) automatically in a DigestPro MSI (INTAVIS Bioanalytical Instruments AG, Köln, Germany). Bands corresponding to extractions from atrial, sinusal, ventricular and conal myocardium were used.
Nano HPLC-ESI-MS/MS analysis
For the analysis of the trypsinized peptides, liquid chromatography was performed using a Nano-HPLC Agilent 1200 comprising a 1200 Series G1379B Degasser, a 1200 Series G2226A Nano pump, 1200 Series G1330B Thermostatted controller, 1200 Series G1377A Micro well plate Autosampler, a 1200 Series G1316B column compartment, a 1200 Series 4208A Gameboy controller and a 1200 Series Solvent Tray (Agilent Technologies, Santa Clara, CA, USA).
Samples were concentrated and desalted in the trapping column C18, ZORBAX 300SB-C18 (Agilent Technologies, Germany), 5x0.3 mm, 5 μm i.d. using a gradient of 98% H2O: 2% acetonitrile (ACN) containing 0.1% formic acid (FA), with a flow rate of 20 μL/min for 6 min. The trapping column was connected in line with the analytical column ZORBAX 300SB-C18 (Agilent Technologies, Germany), 150x0.075mm, 3.5 μm i.d. Elution of the trapping column samples was performed using a gradient of two solvents: A (0.1% FA in H2O) and B (0.1% FA in 80% ACN and H2O), with a flow rate of 300 nL/min. The elution gradient scheme was: isocratic conditions of 98% A: 2% B for 6 min, lineal increase to 50% B in 120 min, lineal increase to 90% B in 2 min, isocratic conditions B 90% for 3 min and back to the initial conditions in 3 min. The ejection volume was 20 μl. Nano-HPLC was interfaced directly with a 3D high capacity ion trap mass spectrometer AmaZon speed ETD (Bruker Daltonics; Billerica, Massachusetts, USA) using CaptiveSpray nanoBooster (Bruker Daltonics; Billerica, Massachusetts, USA) and nano-ESI mode. Ion trap was turned for a target mass of 300-1500 m/z and a smart ICC target of 200.000.
The raw HPLC-MS/MS data obtained were processed using ProteinScape 3 (Brunker Daltonics; Billerica, Massachusetts, USA). For protein identification, MS/MS spectra were compared with UniProtKB/Swiss-Prot database (Release 2017_02, 15-Feb-17, number of entries 553655). Searches were done using Mascot v.3.1 (http://www.matrixscience.com; Matrix Science, London, U.K.). Search parameters were set as follows: carbamidomethyl cysteine as fixed modification and oxidized methionine as variable modification. Peptide mass tolerance was set at ± 0.6 Da in the precursor ions and ± 0.5 Da to the fragmented mases in the MS/MS mode, 1 missed cleavage was allowed. Percolator  was used to assess peptide and protein identification confidence: Protein identifications contain at least two peptides matches that meet or exceed the threshold values for 99% confidence level. The ion score is -10*Log (P), where P is the probability that the observed match is a random event. Peptide identification was performed manually. After database searching, a set of peptide matches was ranked according to their corresponding Mascot scores.
For the quantification of the trypsinized peptides of inflow and outflow segments of three dogfish specimens we used an ESI-Orbitrap-Quadrupole system. The Quadrupole-Orbitrap (Q Exactive Plus Thermo Scientific, Waltham, MA, USA) was coupled online through a nanoESI emitter (10 μm tip; New Objective; Woburn, MA, USA) using a FlexIon nanospray apparatus (Thermo Scientific, Waltham, MA, USA). Data were acquired using XCalibur v3.0 in data-dependent acquisition mode, using a Top20 method. Survey mass range was 300−1750 m/z. The quadrupole isolation window was set to 1.5 m/z; the AGC target was set to 1 × 106, a maximum injection time of 120 ms, and a resolution of 70000. In MS2 scans, target AGC was set to 5 × 105 with a maximum injection time of 60 ms. Protein quantification was achieved by the Label Free Quantification method, using the results from three different samples.
A4.1025 immunoprecipitation and affinity analysis
To estimate the relative affinity of each MyHC isoform for the antibody A4.1025 we conducted an immunoprecipitation assay, using one of the extracts previously subjected to ESI-Quadrupole-Orbitrap quantification. For the immunoprecipitation, A4.1025-coated magnetic microbeads (Miltenyi Biotec, Auburn, CA) were used according to the manufacturer instructions. Quantification of the proteins from the eluate was performed by ESI-Orbitrap-Quadrupole as we previously described. In order to calculate the relative affinity of the four MyHC isoform identified in the previous proteomic analyses, we used extracts from the outflow segments, in which the four sarcomeric MyHC isoforms are significantly represented. The relative affinity of each isomyosin was calculated by dividing the concentration of each isoform in the eluate by the initial concentration of the isoform in the same extract.
Availability of data and materials
The datasets analysed during the current study are available in the UniProtKB/Swiss-Prot repository, https://www.uniprot.org/.
Bovine serum albumin
High-performance liquid chromatography-electrospray ionisation tandem mass spectrometry
Myosin Light Chain
Myosin Heavy Chain
Skeletal muscle microRNA
Phosphate buffer saline
1% bovine serum albumin (BSA) and 0.1% Triton X-100 in TPBS
Sodium dodecyl sulfate polyacrylamide gel electrophoresis
Tris- Phosphate buffer saline
Abu-Daya A, Santer AK, Wells DE, Mohun TJ, Zimmerman LB. Absence of heartbeat in the Xenopus tropicalis mutation muzak is caused by a nonsense mutation in cardiac myosin myh6. Dev Biol. 2009;336(1):20–9. https://doi.org/10.1016/j.ydbio.2009.09.019.
Bader D, Masaki T, Fischman DA. Immunochemical analysis of myosin heavy chain during avian myogenesis in vivo and in vitro. J Cell Biol. 1982;96:763–70. https://doi.org/10.1083/jcb.95.3.763.
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(8):1937–45. https://doi.org/10.1128/MCB.01370-09.
Berdougo E, Coleman H, Lee DH, Stainier DYR, Yelon D. Mutation of weak atrium/atrial myosin heavy chain disrupts atrial function and influences ventricular morphogenesis in zebrafish. Development. 2003;130(24):6121–9. https://doi.org/10.1242/dev.00838.
Bhuiyan SS, Kinoshita S, Wongwarangkana C, Asaduzzaman M, Asakawa S, Watabe S. Evolution of the myosin heavy chain gene MYH14 and its intronic microRNA miR-499: muscle-specific miR-499 expression persists in the absence of the ancestral host gene. BMC Evol Biol. 2013;13(1):142. https://doi.org/10.1186/1471-2148-13-142.
Bisaha JG, Bader D. Identification and characterization of a ventricular-specific avian myosin heavy chain, VMYHC1: Expression in differentiating cardiac and skeletal muscle. Dev Biol. 1991;148:355–64. https://doi.org/10.1016/0012-1606(91)90343-2.
Blagden CS, Currie PD, Ingham PW, Hughes SM. Notochord induction of zebrafish slow muscle mediated by Sonic hedgehog. Genes Dev. 1997;11:2163–75. https://doi.org/10.1101/gad.11.17.2163.
Broughton RE, Betancur RR, Li C, Arratia G, Ortí G. Multi-locus phylogenetic analysis reveals the pattern and tempo of bony fish evolution. PLoS Curr. 2013;5. https://doi.org/10.1371/currents.tol.2ca8041495ffafd0c92756e75247483e.
Daniel JF. The circulation of blood in ammocoetes. Oakland: University of California Press; 1934.
Dan-Goor M, Silberstein L, Kessel M, Muhlrad A. Localization of epitopes and functional effects of two novel monoclonal antibodies against skeletal muscle myosin. J Muscle Res Cell Motil. 1990;11(3):216–26.
De Jong F, Opthof T, Wilde AA, Janse MJ, Charles R, Lamers WH, Moorman AF. Persisting zones of slow impulse conduction in developing chicken hearts. Circ Res. 1992;71(2):240–50.
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(4):375–93. https://doi.org/10.1093/oxfordjournals.molbev.a004093.
Durán AC, Fernández B, Grimes AC, Rodríguez C, Arqué JM, Sans-Coma V. Chondrichtyans have a bulbus arteriosus at the arterial pole of the heart: morphological and evolutionary implications. J Anat. 2008;213:597–606. https://doi.org/10.1111/j.1469-7580.2008.00973.x.
Durán AC, López-Unzu MA, Rodríguez C, Fernández B, Lorenzale M, Linares A, Salmerón F, Sans-Coma V. Structure and vascularization of the ventricular myocardium in Holocephali: their evolutionary significance. J Anat. 2015;226(6):501–10. https://doi.org/10.1111/joa.12317.
England J, Loughna S. Heavy and light roles: myosin in the morphogenesis of the heart. Cell Mol Life Sci. 2013;70:1221–39. https://doi.org/10.1007/s00018-012-1131-1.
Franco D, Gallego A, Habets PE, Sans-Coma V, Moorman AF. Species-specific differences of myosin content in the developing cardiac chambers of fish, birds, and mammals. Anat Rec. 2002;268(1):27–37. https://doi.org/10.1002/ar.10126.
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. https://doi.org/10.1002/dvdy.20460.
Grimes AC. The Cardiac Outflow Tract of Primitive Fishes. In: Zaccone G, Dabrowski K, Hedrick MS, Fernandes JMO, Icardo JM, editors. Phylogeny, Anatomy and Physiology of Ancient Fishes. Boca Raton: CRC Press; 2015. p. 151–78.
Hara Y, Yamaguchi K, Onimaru K, Kadota M, Koyanagi M, Keeley SD, Tatsumi K, Tanaka K, Motone F, Kageyama Y, Nozu R, Adachi N, Nishimura O, Nakagawa R, Tanegashima C, Kiyatake I, Matsumoto R, Murakumo K, Nishida K, Terakita A, Kuratani S, Sato K, Hyodo S, Kuraku S. Shark genomes provide insights into elasmobranch evolution and the origin of vertebrates. Nat Ecol Evol. 2018;2:1761–71.
Harrington WF. On the origin of the contractile force in skeletal muscle. Proc Nat Acad Sci. 1979;76:5066–70.
Harrington WF, Rodgers ME. Myosin. Annu Rev Biochem. 1984;53:35–73. https://doi.org/10.1146/annurev.bi.53.070184.000343.
Hinits Y, Pan L, Walker C, Dowd J, Moens CB, Hughes SM. Zebrafish Mef2ca and Mef2cb are essential for both first and second heart field cardiomyocyte differentiation. Dev Biol. 2012;369:199–210. https://doi.org/10.1016/j.ydbio.2012.06.019.
Hodge TP, Cross R, Lendrick-Jones J. Role of the COOH-terminal nonhelical tailpiece in the assembly of a vertebrate nonmuscle myosin rod. J Cell Biol. 1992;118(5):1085–95. https://doi.org/10.1016/j.ydbio.2012.06.019.
Hu N, Sedmera D, Yost HJ, Clark EB. Structure and function of the developing zebrafish heart. Anat Rec. 2000;260:148–57. https://doi.org/10.1002/1097-0185(20001001)260:2<148::AID-AR50>3.0.CO;2-X.
Ikeda D, Ono Y, Snell P, Edwards YJ, Elgar G, Watabe S. Divergent evolution of the myosin heavy chain gene family in fish and tetrapods: evidence from comparative genomic analysis. Physiol Genom. 2007;32(1):1–5. https://doi.org/10.1152/physiolgenomics.00278.2006.
Keith A. On the fate of the bulbus cordis in the human heart. Lancet. 1924;204(5286):1267–73.
Kikugawa K, Katoh K, Kuraku S, Sakurai H, Ishida O, Iwabe N, Miyata T. Basal jawed vertebrate phylogeny inferred from multiple nuclear DNA-coded genes. BMC Biol. 2004;2(1):3. https://doi.org/10.1186/1741-7007-2-3.
Kinoshita S, Siddique BS, Ceyhun SB, Asaduzzaman M, Asakawa S, Watabe S. Species-specific expression variation of fish MYH14, an ancient vertebrate myosin heavy chain gene orthologue. Fish Sci. 2011;77:847–53. https://doi.org/10.1007/s12562-011-0375-2.
Liang S, Li HC, Wang YX, Wu SS, Cai YJ, Cui HL, Yang YP, Ya J. Pulmonary endoderm, second heart field and the morphogenesis of distal outflow tract in mouse embryonic heart. Dev Growth Differ. 2014;56:276–92. https://doi.org/10.1111/dgd.12129.
Lorenzale M, López-Unzu MA, Rodríguez C, Fernández B, Durán AC, Sans-Coma V. The anatomical components of the cardiac outflow tract of chondrichthyans and actinopterygians. Biol Rev. 2018;93:1604–19. https://doi.org/10.1111/brv.12411.
Pette D, Staron RS. Myosin isoforms, muscle fiber types and transitions. Microsc Res Tech. 2000;50:500–9. https://doi.org/10.1002/1097-0029(20000915)50:6<500::AID-JEMT7>3.0.CO;2-7.
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(6):389–99. https://doi.org/10.1016/S0248-4900(02)00010-2.
Margossian SS, Lowey S. Preparation of myosin and its subfragments from rabbit skeletal muscle. Methods Enzymol. 1982;85(B):55–71.
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. 2009;207:9–20. https://doi.org/10.1016/j.aanat.2016.02.007.
McCarthy JJ, Esser KA, Peterson CA, Dupont-Versteegden EE. Evidence of MyomiR network regulation of β-myosin heavy chain gene expression during skeletal muscle atrophy. Physiol genom. 2009;39(3):219–26. https://doi.org/10.1152/physiolgenomics.00042.2009.
Moore LA, Tidyman WE, Arrizubieta MJ, Bandman E. Gene conversions within the skeletal myosin multigene family. J Mol Biol. 1992;223:368–87. https://doi.org/10.1016/0022-2836(92)90741-2.
Moorman AF, Lamers WH. Molecular anatomy of the developing heart. Trends Cardiovasc Med. 1994;4(6):257–64. https://doi.org/10.1016/1050-1738(94)90029-9.
Morkin E. Control of cardiac myosin heavy chain gene expression. Microsc Res Tech. 2000;50:522–31. https://doi.org/10.1002/1097-0029(20000915)50:6<522::AID-JEMT9>3.0.CO;2-U.
Reiser PJ, Kline WO. Electrophoretic separation and quantitation of cardiac myosin heavy chain isoforms in eight mammalian species. Am J Physiol. 1998;274(3):H1048–53. https://doi.org/10.1152/ajpheart.1998.274.3.H1048.
Rossi A, 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(2):353–64. https://doi.org/10.1113/jphysiol.2009.181008.
Satchell GH, Jones X. The function in the conus arteriosus in the Port Jackson shark, Heterodontus portusjacksoni. J Exp Biol. 1967;46:373–82.
Schib JL, Icardo JM, Durán AC, Guerrero A, López D, Colvee E, De Andrés AV, Sans-Coma V. The conus arteriosus of the adult gilthead seabream (Sparus auratus). J Anat. 2002;201(5):395–404. https://doi.org/10.1046/j.0021-8782.2002.00110.x.
Schindelin J, Arganda-Carreras I, Frise E. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012;9(7):676–82. https://doi.org/10.1038/nmeth.2019.
Sedmera D, Reckova M, deAlmeida A, Sedmerova M, Biermann M, Vlejnik J, Sarre A, Raddatz E, McCarthy RA, Gourdie RG, Thompson RP. Functional and morphological for a ventricular conduction system in zebrafish and Xenopus hearts. Am J Physiol Heart Circ Physiol. 2002;284:H1152–60. https://doi.org/10.1152/ajpheart.00870.2002.
Sedmera D, Wessels A, Trusk TC, Thompson RP, Hewett KW, Gourdie RG. Changes in activation sequence of embryonic chick atria correlate with developing myocardial architecture. Am J Physiol-Heart Circul Physiol. 2006;291:H1646–52. https://doi.org/10.1152/ajpheart.01007.2005.
Sedmera D. Function and form in the developing cardiovascular system. Car Res. 2011;91:252–9. https://doi.org/10.1093/cvr/cvr062.
Sellers JR. Myosin: a diverse superfamily. Biochem Biophys Acta. 2000;1496(1):3–22. https://doi.org/10.1016/S0167-4889(00)00005-7.
Smith SJ, Towers N, Saldanha JW, Shang CA, Mahmood SR, Taylor WR, Mohun TJ. The cardiac-restricted protein ADP-ribosylhydrolase-like 1 isessential for heart chamber outgrowth and acts on muscle actin filament assembly. Dev Biol. 2016;416:373–88. https://doi.org/10.1016/j.ydbio.2016.05.006.
Somi S, Klein ATJ, Houweling AC, Ruijter JM, Buffing AAM, Moorman AFM, Van den Hoff MJB. Atrial and ventricular myosin heavy-chain expression in the developing chicken heart: strengths and limitations of non-radioactive in situ hybridization. J Histochem Cytochem. 2006;54(6):649–64. https://doi.org/10.1369/jhc.5A6846.2006.
Sudak FN. Intrapericardial and intracardiac pressures and the events of the cardiac cycle in Mustelus canis (Mitchill). Comp Biochem Physiol. 1965;14(4):689–705. https://doi.org/10.1016/0010-406X(65)90255-0.
Takeuchi M, Okabe M, Aizawa S. The genus Polypterus (Bichirs): a fish group diverged at the stem of ray-finned fishes (Actinopterygii). In: Crotty DA, Gann A, editors. Emerging Model Organisms: A Laboratory Manual, vol. 1. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press; 2009. p. 447–67. https://doi.org/10.1101/pdb.emo117.
Thomas PS, Sridurongrit S, Ruiz-Lozano P, Kaartinen V. Deficient signaling via Alk2 (Acvr1) leads to bicuspid aortic valve development. PLoS One. 2016;7(4):e35539. https://doi.org/10.1371/journal.pone.0035539.
Van Rooij E, Quiat D, Johnson AB, Sutherland LB, Qi X, Richardson JA, Kelm RJ, Olson EN. A family of micrornas encoded by myosin genes governs myosin expression and muscle performance. Dev Cell. 2009;17:662–73. https://doi.org/10.1016/j.devcel.2009.10.013.
Venkatesh B, Lee AP, Ravi V, Maurya AK, Lian MM, Swann JB, Ohta Y, Flajnik MF, Sutoh Y, Kasahara M, et al. Elephant shark genome provides unique insights into gnathostome evolution. Nature. 2014;505:174–9. https://doi.org/10.1038/nature12826.
Wang X, Ono Y, Tan SC, Chai RJ, Parkin C, Ingham PW. Prdm1a and miR-499 act sequentially to restrict Sox6 activity to the fast-twitch muscle lineage in the zebrafish embryo. Development. 2011;138(20):4399–404. https://doi.org/10.1242/dev.070516.
Warkman AS, Whitman SA, Miller MK, Garriock RJ, Schwach CM, Gregorio CC, Krieg PA. Developmental expression and cardiac transcriptional regulation of myh7b, a third myosin heavy chain in the vertebrate heart. Cytoskeleton. 2012a;69:324–35. https://doi.org/10.1002/cm.21029.
Warkman AS, Whitman SA, Miller MK, Garriock RJ, Schwach CM, Gregorio CC, Krieg PA. Correction: developmental expression and cardiac transcriptional regulation of myh7b, a third myosin heavy chain in the vertebrate heart. Cytoskeleton. 2012b;69:1086. https://doi.org/10.1002/cm.21089.
Wright JC, Collins MO, Yu L, Käll L, Brosch M, Choudhary JS. Enhanced peptide identification by electron transfer dissociation using an improved MASCOT percolator. Mol Cell Proteomics. 2012;11(8):478–91. https://doi.org/10.1074/mcp.O111.014522.
Yelon D, Horne SA, Stainier YR. Restricted expression of cardiac myosin genes reveals regulated aspects of heart tube assembly in zebrafish. Dev Biol. 1999;214:23–37. https://doi.org/10.1006/dbio.1999.9406.
Yutzey KE, Rhee JT, Bader D. Expression of the atrial-specific myosin heavy chain AMYHC1 and the establishment of anteroposterior polarity in the developing chicken heart. Development. 1994;120:871–83.
Ziermman JM, Freitas R, Diogo R. Muscle development in the shark Scyliorhinus canicula: implications for the evolution of the gnathostome head and paired appendage musculature. Front Zool. 2017;14:31–48. https://doi.org/10.1186/s12983-017-0216-y.
We thank Dr. Jorge Cubo for his valuable suggestions and discussions, and Dr. Javier Márquez for helping with the interpretation of the proteomics results. We also thank Mrs. Remedios Crespillo and Dr. David Navas from the “Servicios Centrales de Apoyo a la Investigación” (SCAI, University of Málaga, Spain) and Dr. Casimiro Cárdenas and Dr. Mercedes Martín from the Proteomics Core Facility (SCAI, University of Málaga, Spain) for their kind assistance in proteomics and microscopy techniques. Dr. Ana Ramos (Spanish Institute of Oceanography, head of the MAURIT campaign), Dr. Alberto Domezain (Sierra Nevada Fishery at Riofrío, Granada, Spain), Dr. Miguel Lorenzale, Dr. Cristina Rodríguez and Dr. M. Carmen Fernández (University of Málaga, Málaga, Spain), Mr. Sebastián Martín, Mr. Marcos Fernández (Caleta de Vélez, Málaga, Spain) and ‘Cofradía de Pescadores Caleta de Vélez’ (Málaga, Spain) kindly contributed in obtaining the material.
This work was supported by Subdirección General de Proyectos de Investigación, Ministerio de Economía y Competitividad (Grant number CGL2017-85090-P), Ministerio de Educación, Cultura y Deporte (Fellowship number FPU15/03209 to M.A.L.-U.), Sistema Nacional de Garantía Juvenil (Contract UMAJ175 to M.T.S.-N.), Universidad de Málaga and FEDER funds.
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.