Multiple gene organization of pufferfish Fugu rubripes tropomyosin isoforms and tissue distribution of their transcripts

Qualitative and quantitative evaluation of TPM transcripts and proteins in developing striated chicken muscles indicate TPM4α is the major sarcomeric cardiac tropomyosin from early embryonic life to adulthood

The chicken has been used since the 1980s as an animal model for developmental studies regarding tropomyosin isoform diversity in striated muscles, however, the pattern of expression of transcripts as well as the corresponding TPM proteins of various tropomyosin isoforms in avian hearts are not well documented. In this study, using conventional and qRT-PCR, we report the expression of transcripts for various sarcomeric TPM isoforms in striated muscles through development. Transcripts of both TPM1α and TPM1κ, the two sarcomeric isoforms of the TPM1 gene, are expressed in embryonic chicken hearts but disappear in post hatch stages. TPM1α transcripts are expressed in embryonic and adult skeletal muscle. The sarcomeric isoform of the TPM2 gene is expressed mostly in embryonic skeletal muscles. As reported earlier, TPM3α is expressed in embryonic heart and skeletal muscle but significantly lower in adult striated muscle. TPM4α transcripts are expressed from embryonic to adult chicken hearts but not in skeletal muscle. Our 2D Western blot analyses using CH1 monoclonal antibody followed by mass spectra evaluations found TPM4α protein is the major sarcomeric tropomysin expressed in embryonic chicken hearts. However, in 7-day-old embryonic hearts, a minute quantity of TPM1α or TPM1κ is also expressed. This finding suggests that sarcomeric TPM1 protein may play some important role in cardiac contractility and/or cardiac morphogenesis during embryogenesis. Since only the transcripts of TPM4α are expressed in adult chicken hearts, it is logical to presume that TPM4α is the only sarcomeric TPM protein produced in adult cardiac tissues.

Identification, characterization, and expression of sarcomeric tropomyosin isoforms in zebrafish

Tropomyosin is a component of thin filaments that constitute myofibrils, the contractile apparatus of striated muscles. In vertebrates, except for fish, four TPM genes TPM1, TPM2, TPM3, and TPM4 are known. In zebrafish, there are six TPM genes that include the paralogs of the TPM1 (TPM1-1 and TPM1-2), the paralogs of the TPM4 gene (TPM4-1 and TPM4-2), and the two single copy genes TPM2 and TPM3. In this study, we have identified, cloned, and sequenced the TPM1-1κ isoform of the TPM1-1 gene and also discovered a new isoform TPM1-2ν of the TPM1-2. Further, we have cloned and sequenced the sarcomeric isoform of the TPM4-2 gene designated as TPM4-2α. Using conventional RT-PCR, we have shown the expression of the sarcomeric isoforms of TPM1-1, TPM1-2, TPM2, TPM3, TPM4-1, and TPM4-2 in heart and skeletal muscles. By qRT-PCR using both relative expression as well as the absolute copy number, we have shown that TPM1-1α, TPM1-2α, and TPM1-2ν are expressed mostly in skeletal muscle; the level of expression of TPM1-1κ is significantly lower compared to TPM1-1α in skeletal muscle. In addition, both TPM4-1α and TPM4-2α are predominantly expressed in heart. 2D Western blot analyses using anti-TPM antibody followed by Mass Spectrometry of the proteins from the antibody-stained spots show that TPM1-1α and TPM3α are expressed in skeletal muscle whereas TPM4-1α and TPM3α are expressed in zebrafish heart. To the best of our knowledge, this is by far the most comprehensive analysis of tropomyosin expression in zebrafish, one of the most popular animal models for gene expression study.

Population genomics of local adaptation versus speciation in coral reef fishes (Hypoplectrus spp, Serranidae)

Are the population genomic patterns underlying local adaptation and the early stages of speciation similar? Addressing this question requires a system in which (i) local adaptation and the early stages of speciation can be clearly identified and distinguished, (ii) the amount of genetic divergence driven by the two processes is similar, and (iii) comparisons can be repeated both taxonomically (for local adaptation) and geographically (for speciation). Here, we report just such a situation in the hamlets (Hypoplectrus spp), brightly colored reef fishes from the wider Caribbean. Close to 100,000 SNPs genotyped in 126 individuals from three sympatric species sampled in three repeated populations provide genome-wide levels of divergence that are comparable among allopatric populations (F st estimate = 0.0042) and sympatric species (F st estimate = 0.0038). Population genetic, clustering, and phylogenetic analyses reveal very similar patterns for local adaptation and speciation, with a large fraction of the genome undifferentiated (F st estimate ≈ 0), a very small proportion of F st outlier loci (0.05-0.07%), and remarkably few repeated outliers (1-3). Nevertheless, different loci appear to be involved in the two processes in Hypoplectrus, with only 7% of the most differentiated SNPs and outliers shared between populations and species comparisons. In particular, a tropomyosin (Tpm4) and a previously identified hox (HoxCa) locus emerge as candidate loci (repeated outliers) for local adaptation and speciation, respectively. We conclude that marine populations may be locally adapted notwithstanding shallow levels of genetic divergence, and that from a population genomic perspective, this process does not appear to differ fundamentally from the early stages of speciation.

Tissue expression of PPAR-α isoforms in Scophthalmus maximus and transcriptional response of target genes in the heart after exposure to WY-14643

Peroxisome proliferator-activated receptors (PPARs) are involved in the regulation of lipid and carbohydrate metabolism and can be activated either by natural ligands as fatty acids or by synthetic ligands including several environmental chemicals. In this study, two PPARα isoforms (α1 and α2) were analyzed in turbot (Scophthalmus maximus) for a different tissue distribution. PPARα1 was ubiquitously expressed, while the PPARα2 was predominantly expressed in the heart. Following this result, turbot juveniles were exposed by injection to a synthetic selective PPARα agonist, WY-14643, for 14 days. Suppression subtractive hybridization (SSH) was performed with pools of heart samples of control and exposed fish to get insights into PPARα-regulated genes in the heart of juvenile turbot. Four genes were positively identified in the forward-subtracted and 12 genes in the reverse-subtracted cDNA SSH library, corresponding to the down-regulated and up-regulated genes in response to the WY-14643 treatment, respectively. The confirmation of these results in individual samples of juvenile turbot exposed to WY-14643 revealed a statistically significant mRNA induction of two cardiac muscle proteins (myosin light chain 2 and tropomyosin 4), which were shown to be involved in heart contraction and heartbeat regulation in other teleost species. Herewith, we showed for the first time that PPARα2 is predominantly expressed in the heart and that a PPARα agonist can induce the mRNA expression of cardiac muscle proteins in teleosts.

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

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

Fast skeletal muscle myosin heavy chain gene cluster of medaka Oryzias latipes enrolled in temperature adaptation

To disclose mechanisms involved in temperature acclimation of fish muscle, we subjected eurythermal fish of medaka Oryzias latipes to cloning of myosin heavy chain genes (MYHs). We cloned cDNAs encoding fast skeletal muscle myosin heavy chain (MYH) isoforms from cDNA libraries of medaka acclimated to 10 and 30 degrees C and observed that different MYH cDNA clones are expressed in the two temperature-acclimated fish. Subsequently, we isolated several overlapping MYH contigs by shotgun cloning strategy from a medaka genomic library. Contig assembly of the complete medaka MYH (mMYH) locus of 219 kbp revealed a cluster of tandemly arrayed 11 mMYHs, in which eight genes are actually transcribed, with the remaining three being pseudogenes. Expression analysis of the transcribed genes revealed that two genes were each highly expressed in medaka acclimated to 10 and 30 degrees C, whereas comparatively lower expression levels of the three genes were exclusively observed in medaka acclimated to 30 degrees C. cDNAs of the remaining genes were too underrepresented in the libraries to determine the expression levels, and the transcripts could only be obtained by reverse transcription-polymerase chain reaction. Deduced amino acid sequences in the loop 1 and loop 2 regions of mMYHs were highly variable, suggesting that these isoforms were functionally different. The present findings consolidate our knowledge on teleost MYH multigene family and would provide further insight into the mechanisms by which expressions of individual MYH molecules are fine-tuned with environmental temperature fluctuations with further functional analysis of the genes concerned.

Physiological changes and differential gene expression in mummichogs (Fundulus heteroclitus) exposed to arsenic

Arsenic has been detected as a contaminant in water bodies around the world. Although a number of studies have shown toxicity to adult fish, little is known about its effects on the offspring. However, human epidemiological studies have shown that arsenic increases the number of stillbirths and prematurely born infants. We examined changes in the morphology and gene expression in juvenile mummichogs (Fundulus heteroclitus) whose parents were exposed to 230 ppb arsenic for 10 days immediately prior to spawning. The hatchlings of exposed fish had a 2.8-fold increased incidence of curved or stunted tails. Total RNA from 6-week-old hatchlings, reared in clean water, was used to construct a cDNA subtractive hybridization library. Using this library, we found 13 genes whose expression was altered in the hatchlings as a result of arsenic exposure. We confirmed differential expression by real-time PCR and found significant up-regulation of myosin light chain 2 (4.2-fold), type II keratin (1.5-fold), tropomyosin (3.1-fold) and parvalbumin (3.5-fold) in the hatchlings whose parents were exposed to arsenic. These genes are important during embryogenesis and their differential expression may be linked to the morphological changes observed in the hatchlings.

Diversity of the pufferfish Takifugu rubripes fast skeletal myosin heavy chain genes

Myosin is a highly conserved, ubiquitous actin-based molecular motor that is distributed as diverse as from prokaryotes to mammalian tissues. Among various types in the myosin family proteins, class II, also called sarcomeric, myosin is a classical, conventional molecule that has been extensively studies on its functional and structural properties. It consists of two heavy chains (MYH) of about 200 kDa and four light chains of about 20 kDa. The exon-intron organization was determined for the major subunit of MYH, which contains ATP-hydrolysis and actin-binding sites, from torafugu (tiger pufferfish) Takifugu rubripes fast skeletal muscles. Comprehensive investigation for fast skeletal MYHs based on the fugu (torafugu) genome database and subsequent construction of their physical map revealed that torafugu contains at least 8 putative skeletal MYHs. Furthermore, genomic structural analysis revealed that skeletal MYHs are not clustered in a single locus, but rather spread to at least four loci, with two of them locating at the mammalian syntenic regions. Such arrangement of torafugu MYHs are in a marked contrast to mammalian fast skeletal MYHs that are clustered in a single locus. These data suggest that an ancient segmental duplication or whole-genome duplication occurred in fish lineage as in many other reported torafugu genes.

Fish fast skeletal muscle tropomyosins show species-specific thermal stability

Tropomyosin (TM) was isolated from the fast skeletal muscle of six fish species, whose amino acid sequences of this protein have already been revealed. The thermal stability of these TMs was measured by differential scanning calorimetry (DSC) and circular dichroism (CD), while the molecular weights were measured by mass spectrometry. The results showed clear differences in thermostability among these fish TMs, though the identity of amino acid sequences was more than 93.3%. Therefore, only a few amino acid substitutions could affect the overall stability of the TM molecule. Especially, several residues located on the molecular surface were considered to be responsible for such stability difference. In contrast, the molecular weights of these TMs as measured by mass spectrometry were higher than those calculated from amino acid composition, suggesting the presence of post-translational modification(s) which could also affect their thermal stability.

The complete genomic sequence of the carp fast skeletal myosin heavy chain gene

We have determined the complete DNA nucleotide sequence of the carp Cyprinus carpio fast skeletal myosin heavy chain (MYH) gene. Introns and exons were predicted by comparison with the corresponding carp MYH cDNAs previously reported. The gene encoded the entire mRNA transcript and contained 5958 nucleotides (nt) including 77 nt 5'-untranslated region, 5796 nt coding region for 1931 amino acids, and 85 nt 3'-untranslated region. The coding region was split by 38 introns and the complete gene contained 11,385 nt. This integration of the carp fast skeletal MYH gene was comparable to those of the rat and chicken embryonic MYH genes, which have 41 and 40 exons, respectively. However, the entire gene size of carp MYH was about half those of rat and chicken due to much smaller size of carp introns. We have also demonstrated that this carp MYH gene belonged to so-called intermediate type in a multigene family of carp fast skeletal muscle MYH in comparison of its nucleotide and deduced amino acid sequences to those of carp MYH cDNAs reported previously.

Characterization of the pufferfish Takifugu rubripes apolipoprotein multigene family

We have characterized the apolipoprotein multigene family of the pufferfish Takifugu rubripes. The pufferfish mainly contains 28-kDa, 27-kDa, and 14-kDa apolipoproteins in its plasma and was designated apo-28 kDa, apo-27 kDa, and apo-14 kDa, respectively. N-terminal amino acid sequencing revealed that pufferfish apo-28 kDa and apo-27 kDa have an identical amino acid sequence except an additional propeptide in the former; and both are homologues of apoA-I from other animals. The sequence of pufferfish apo-14 kDa is homologous to that of eel apo-14 kDa previously reported, both being apparently specific to fish. In silico screening, using the publicly available Fugu genome database confirmed the pufferfish apoA-I and apo-14 kDa genes. The database further contained the genes encoding four types of apoA-IV, one apoC-II and two types of apoE. Thus, pufferfish contains nine genes encoding apolipoprotein multigene family. Two apoA-IV and one apoE genes were tandemly arrayed and located on one scaffold. Thus two sets of these genes formed two gene clusters. The apoC-II and apo-14 kDa genes are also located on a single scaffold. apoA-I and apo-14 kDa gene transcripts were mainly expressed in liver and less abundantly in brain. The transcripts of the former gene were also observed in intestine. In contrast, the transcripts encoding four apoA-IVs, one apoC-II, and two apoEs were mainly expressed in intestine. These structural details of pufferfish apolipoproteins and tissue distribution of their gene transcripts provide a novel evidence for better understanding of evolutionary relationships of apolipoprotein multigene family.

Identification of three duplicated Spin genes in medaka (Oryzias latipes)

Gene and genomic duplications are very important and frequent events in fish evolution, and the divergence of duplicated genes in sequences and functions is a focus of research on gene evolution. Here, we report the identification and characterization of three duplicated Spindlin (Spin) genes from medaka (Oryzias latipes): OlSpinA, OlSpinB, and OlSpinC. Molecular cloning, genomic DNA Blast analysis and phylogenetic relationship analysis demonstrated that the three duplicated OlSpin genes should belong to gene duplication. Furthermore, Western blot analysis revealed significant expression differences of the three OlSpins among different tissues and during embryogenesis in medaka, and suggested that sequence and functional divergence might have occurred in evolution among them.

Subfunction partitioning, the teleost radiation and the annotation of the human genome

Half of all vertebrate species are teleost fish. What accounts for this explosion of biodiversity? Recent evidence and advances in evolutionary theory suggest that genomic features could have played a significant role in the teleost radiation. This review examines evidence for an ancient whole-genome duplication (tetraploidization) event that probably occurred just before the teleost radiation. The partitioning of ancestral subfunctions between gene copies arising from this duplication could have contributed to the genetic isolation of populations, to lineage-specific diversification of developmental programs, and ultimately to phenotypic variation among teleost fish. Beyond its importance for understanding mechanisms that generate biodiversity, the partitioning of subfunctions between teleost co-orthologs of human genes can facilitate the identification of tissue-specific conserved noncoding regions and can simplify the analysis of ancestral gene functions obscured by pleiotropy or haploinsufficiency. Applying these principles on a genomic scale can accelerate the functional annotation of the human genome and understanding of the roles of human genes in health and disease.

Microarray gene expression profiling during the segmentation phase of zebrafish development

We analyzed 15,512 unique transcripts from wild-type Danio rerio using a long oligonucleotide microarray containing >16,000 65-mers probes. Total RNA was isolated from staged embryos at 2 h intervals over a 24-h period. On average, at any given time point, 27% of the probe set detected corresponding transcripts in embryonic RNA. There were two predominant patterns in the nearly 4000 genes that changed expression in at least one time point during the first 24 hpf. At 12 hpf, we detected 420 up-regulated and 386 down-regulated genes. By 24 hpf, the number of up- and down-regulated genes had increased to 954 and 766, respectively. While the majority of these genes maintained their new level of expression for the duration of the time course, we identified five genes with phasic regulation over the 24-h time course. Two of these genes, germ cell nuclear factor and mesogenin, have been identified as being expressed during gastrulation (5 1/4 to 10 h postfertilization) and subsequently repressed. A cluster containing 36 distinct ribosomal proteins was up-regulated at 12 h, indicating a capability for de novo protein synthesis during and after this stage. Twenty-three muscle-specific genes were up-regulated late during the initial 24 hpf, corresponding to the development and differentiation of the somites.