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

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

Evolutionary analyses of genes in Echinodermata offer insights towards the origin of metazoan phyla

Despite recent studies discussing the evolutionary impacts of gene duplications and losses among metazoans, the genomic basis for the evolution of phyla remains enigmatic. Here, we employ phylogenomic approaches to search for orthologous genes without known functions among echinoderms, and subsequently use them to guide the identification of their homologs across other metazoans. Our final set of 14 genes was obtained via a suite of homology prediction tools, gene expression data, gene ontology, and generating the Strongylocentrotus purpuratus phylome. The gene set was subjected to selection pressure analyses, which indicated that they are highly conserved and under negative selection. Their presence across broad taxonomic depths suggests that genes required to form a phylum are ancestral to that phylum. Therefore, rather than de novo gene genesis, we posit that evolutionary forces such as selection on existing genomic elements over large timescales may drive divergence and contribute to the emergence of phyla.     

The organization of cytoplasmic ribosomal protein genes in the Arabidopsis genome

Eukaryotic ribosomes are made of two components, four ribosomal RNAs, and approximately 80 ribosomal proteins (r-proteins). The exact number of r-proteins and r-protein genes in higher plants is not known. The strong conservation in eukaryotic r-protein primary sequence allowed us to use the well-characterized rat (Rattus norvegicus) r-protein set to identify orthologues on the five haploid chromosomes of Arabidopsis. By use of the numerous expressed sequence tag (EST) accessions and the complete genomic sequence of this species, we identified 249 genes (including some pseudogenes) corresponding to 80 (32 small subunit and 48 large subunit) cytoplasmic r-protein types. None of the r-protein genes are single copy and most are encoded by three or four expressed genes, indicative of the internal duplication of the Arabidopsis genome. The r-proteins are distributed throughout the genome. Inspection of genes in the vicinity of r-protein gene family members confirms extensive duplications of large chromosome fragments and sheds light on the evolutionary history of the Arabidopsis genome. Examination of large duplicated regions indicated that a significant fraction of the r-protein genes have been either lost from one of the duplicated fragments or inserted after the initial duplication event. Only 52 r-protein genes lack a matching EST accession, and 19 of these contain incomplete open reading frames, confirming that most genes are expressed. Assessment of cognate EST numbers suggests that r-protein gene family members are differentially expressed.     

Phylogenetic and expression analysis of the basic helix-loop-helix transcription factor gene family: genomic approach to cellular differentiation

A phylogenetic analysis of seven different species (human, mouse, rat, worm, fly, yeast, and plant) utilizing all (541) basic helix-loop-helix (bHLH) genes identified, including expressed sequence tags (EST), was performed. A super-tree involving six clades and a structural categorization involving the entire coding sequence was established. A nomenclature was developed based on clade distribution to discuss the functional and ancestral relationships of all the genes. The position/location of specific genes on the phylogenetic tree in relation to known bHLH factors allows for predictions of the potential functions of uncharacterized bHLH factors, including EST's. A genomic analysis using microarrays for four different mouse cell types (i.e. Sertoli, Schwann, thymic, and muscle) was performed and considered all known bHLH family members on the microarray for comparison. Cell-specific groups of bHLH genes helped clarify those bHLH genes potentially involved in cell specific differentiation. This phylogenetic and genomic analysis of the bHLH gene family has revealed unique aspects of the evolution and functional relationships of the different genes in the bHLH gene family.     

Analysis of 142 genes resolves the rapid diversification of the rice genus

Background

The completion of rice genome sequencing has made rice and its wild relatives an attractive system for biological studies. Despite great efforts, phylogenetic relationships among genome types and species in the rice genus have not been fully resolved. To take full advantage of rice genome resources for biological research and rice breeding, we will benefit from the availability of a robust phylogeny of the rice genus.

Results

Through screening rice genome sequences, we sampled and sequenced 142 single-copy genes to clarify the relationships among all diploid genome types of the rice genus. The analysis identified two short internal branches around which most previous phylogenetic inconsistency emerged. These represent two episodes of rapid speciation that occurred approximately 5 and 10 million years ago (Mya) and gave rise to almost the entire diversity of the genus. The known chromosomal distribution of the sampled genes allowed the documentation of whole-genome sorting of ancestral alleles during the rapid speciation, which was responsible primarily for extensive incongruence between gene phylogenies and persisting phylogenetic ambiguity in the genus. Random sample analysis showed that 120 genes with an average length of 874 bp were needed to resolve both short branches with 95% confidence.

Conclusion

Our phylogenomic analysis successfully resolved the phylogeny of rice genome types, which lays a solid foundation for comparative and functional genomic studies of rice and its relatives. This study also highlights that organismal genomes might be mosaics of conflicting genealogies because of rapid speciation and demonstrates the power of phylogenomics in the reconstruction of rapid diversification.     

Identification and phylogenetic sorting of bacterial lineages with universally conserved genes and proteins

Molecular characterizations of bacteria often employ ribosomal DNA (rDNA) to establish the identity and relationships among organisms, but the use of rRNA sequences can be problematic as the result of alignment ambiguities caused by indels, the lack of informative characters, and varying functional constraints over the molecule. Although protein-coding regions have been used as an alternative to rRNA, there is neither consensus among the genes examined nor ways to rapidly obtain sequence information for such genes from uncharacterized bacterial species. To standardize the set of protein-coding loci assayed in bacterial genomes, we examined over 100 widely distributed genes to identify sets of universal primers for use in the PCR amplification of protein coding regions that are common to virtually all bacteria. From this set, we developed primer sets that each target of 10 genes spanning an array of genomic locations and functional categories. Although many of the primers contain sequence degeneracies that aid in targeting genes across diverse taxa, most are adequate for direct sequencing of amplification products, thereby eliminating intermediate cloning before sequence determination. We foresee the analysis of these protein-coding regions as being complementary to ribosomal DNA for answering questions pertaining to bacterial identification, classification, phylogenetics and evolution.     

Phylogenetic distribution of catalase-peroxidases: are there patches of order in chaos?

Hydrogen peroxide features in many biological oxidative processes and must be continuously degraded enzymatically either via a catalatic or a peroxidatic mechanism. For this purpose ancestral bacteria evolved a battery of different heme and non-heme enzymes, among which heme-containing catalase-peroxidases (CP) are one of the most widespread representatives. They are unique since they can follow both H(2)O(2)-degrading mechanisms, the catalase activity being clearly dominant. With the fast increasing amount of genomic data available, we were able to perform an extensive search for CP and found almost 300 sequences covering a large range of microorganisms. Most of them were encoded by bacterial genomes, but we could also find some in eukaryotic organisms other than fungi, which has never been shown until now. Our screen also reveals that approximately 60% of the bacteria do not possess CP genes. Chaotic distribution among species and incongruous phylogenetic reconstruction indicated existence of numerous lateral gene transfers in addition to duplication events and regular speciation. The results obtained show an impressively complex gene transmission pattern, and give some new insights about the role of CP and the origin of life on earth. Finally, we propose for the first time bacterial candidates that may have participated in the transfer of CP from bacteria to eukaryotes.     

Genomics and homeostasis

The Cannon lecture this year illustrates how knowledge of DNA sequences of complex living organisms is beginning to shape the landscape of physiology in the 21st century. Enormous challenges and opportunities now exist for physiologists to relate the galaxy of genes to normal and pathological functions. The first extensive genomic systems biology map for cardiovascular and renal function was completed last year as well as a new hypothesis-generating tool ("physiological profiling") that enables us to hypothesize relationships between specific genes responsible for the regulation of regulatory pathways. Techniques of chromosomal substitution (consomic and congenic rats) are beginning to confirm statistical results from linkage analysis studies, narrow the regions of genetic interest for positional cloning, and provide genetically well-defined control strains for physiological studies. Patterns of gene expression identified by microarray and mapping of expressed genes to chromosomal sites are adding to the understanding of systems physiology. The previously unimaginable goal of connecting approximately 36,000 genes to the complex functions of mammalian systems is indeed well underway.     

Extracting three-way gene interactions from microarray data

MOTIVATION: It is an important and difficult task to extract gene network information from high-throughput genomic data. A common approach is to cluster genes using pairwise correlation as a distance metric. However, pairwise correlation is clearly too simplistic to describe the complex relationships among real genes since co-expression relationships are often restricted to a specific set of biological conditions/processes. In this study, we described a three-way gene interaction model that captures the dynamic nature of co-expression relationship between a gene pair through the introduction of a controller gene. RESULTS: We surveyed 0.4 billion possible three-way interactions among 1000 genes in a microarray dataset containing 678 human cancer samples. To test the reproducibility and statistical significance of our results, we randomly split the samples into a training set and a testing set. We found that the gene triplets with the strongest interactions (i.e. with the smallest P-values from appropriate statistical tests) in the training set also had the strongest interactions in the testing set. A distinctive pattern of three-way interaction emerged from these gene triplets: depending on the third gene being expressed or not, the remaining two genes can be either co-expressed or mutually exclusive (i.e. expression of either one of them would repress the other). Such three-way interactions can exist without apparent pairwise correlations. The identified three-way interactions may constitute candidates for further experimentation using techniques such as RNA interference, so that novel gene network or pathways could be identified.     

Lifestyle evolution in symbiotic bacteria: insights from genomics

Bacteria that live only in eukaryotic cells and tissues, including chronic pathogens and mutualistic bacteriocyte associates, often possess a distinctive set of genomic traits, including reduced genome size, biased nucleotide base composition and fast polypeptide evolution. These phylogenetically diverse bacteria have lost certain functional categories of genes, including DNA repair genes, which affect mutational patterns. However, pathogens and mutualistic symbionts retain loci that underlie their unique interaction types, such as genes enabling nutrient provisioning by mutualistic bacteria-inhabiting animals. Recent genomic studies suggest that many of these bacteria are irreversibly specialized, precluding shifts between pathogenesis and mutualism.     

Incorporating 16S Gene Copy Number Information Improves Estimates of Microbial Diversity and Abundance

The abundance of different SSU rRNA (“16S”) gene sequences in environmental samples is widely used in studies of microbial ecology as a measure of microbial community structure and diversity. However, the genomic copy number of the 16S gene varies greatly – from one in many species to up to 15 in some bacteria and to hundreds in some microbial eukaryotes. As a result of this variation the relative abundance of 16S genes in environmental samples can be attributed both to variation in the relative abundance of different organisms, and to variation in genomic 16S copy number among those organisms. Despite this fact, many studies assume that the abundance of 16S gene sequences is a surrogate measure of the relative abundance of the organisms containing those sequences. Here we present a method that uses data on sequences and genomic copy number of 16S genes along with phylogenetic placement and ancestral state estimation to estimate organismal abundances from environmental DNA sequence data. We use theory and simulations to demonstrate that 16S genomic copy number can be accurately estimated from the short reads typically obtained from high-throughput environmental sequencing of the 16S gene, and that organismal abundances in microbial communities are more strongly correlated with estimated abundances obtained from our method than with gene abundances. We re-analyze several published empirical data sets and demonstrate that the use of gene abundance versus estimated organismal abundance can lead to different inferences about community diversity and structure and the identity of the dominant taxa in microbial communities. Our approach will allow microbial ecologists to make more accurate inferences about microbial diversity and abundance based on 16S sequence data.     

Novel Relationships Among Ten Fish Model Species Revealed Based on a Phylogenomic Analysis Using ESTs

The power of comparative phylogenomic analyses also depends on the amount of data that are included in such studies. We used expressed sequence tags (ESTs) from fish model species as a proof of principle approach in order to test the reliability of using ESTs for phylogenetic inference. As expected, the robustness increases with the amount of sequences. Although some progress has been made in the elucidation of the phylogeny of teleosts, relationships among the main lineages of the derived fish (Euteleostei) remain poorly defined and are still debated. We performed a phylogenomic analysis of a set of 42 of orthologous genes from 10 available fish model systems from seven different orders (Salmoniformes, Siluriformes, Cypriniformes, Tetraodontiformes, Cyprinodontiformes, Beloniformes, and Perciformes) of euteleostean fish to estimate divergence times and evolutionary relationships among those lineages. All 10 fish species serve as models for developmental, aquaculture, genomic, and comparative genetic studies. The phylogenetic signal and the strength of the contribution of each of the 42 orthologous genes were estimated with randomly chosen data subsets. Our study revealed a molecular phylogeny of higher-level relationships of derived teleosts, which indicates that the use of multiple genes produces robust phylogenies, a finding that is expected to apply to other phylogenetic issues among distantly related taxa. Our phylogenomic analyses confirm that the euteleostean superorders Ostariophysi and Acanthopterygii are monophyletic and the Protacanthopterygii and Ostariophysi are sister clades. In addition, and contrary to the traditional phylogenetic hypothesis, our analyses determine that killifish (Cyprinodontiformes), medaka (Beloniformes), and cichlids (Perciformes) appear to be more closely related to each other than either of them is to pufferfish (Tetraodontiformes). All 10 lineages split before or during the fragmentation of the supercontinent Pangea in the Jurassic. [Reviewing Editor: Dr. Rafael Zardoya]     

Heavy chain variable region gene families evolved early in phylogeny. Ig complexity in fish

The V regions of channel catfish H chain cDNA clones have been analyzed. Based upon sequence relationships and hybridization analyses, five different groups of VH genes are identified whose definition is consistent with that of five different VH families. Genomic Southern blots indicate that as many as 100 different germ-line VH genes are likely represented by these families. The sequence diversity between identified members of these different families is similar in magnitude to the divergence represented between members of different human or mouse VH families. The FR regions are the most conserved regions when members of different catfish VH families are compared; specific amino acid positions appear to be highly conserved in phylogeny. Equally important is that diversity is represented in complementarity-determining regions CDR1 and CDR2 in members of the different families as well as in members of the same VH family. These results suggest that an extensive repertoire of VH genes can contribute to antibody diversity in this lower vertebrate. Sequence comparisons indicate that one of the catfish VH families shares considerable structural similarity to several higher vertebrate VH gene families--a relationship which suggests that this VH family may be ancestral to some VH gene families of higher vertebrates. Characteristic of the genomic organization of higher vertebrate H chains, catfish appear to have different VH families wherein a VH gene likely undergoes functional recombination with putative DH gene segments and one of apparently several different JH segments. The recombined V region is expressed with the same C region gene. These combined results suggest that bony fishes are the earliest known phylogenetic representatives to have evolved extensive V region gene families.     

Evolution of the globin gene family in deuterostomes: lineage-specific patterns of diversification and attrition

In the Metazoa, globin proteins display an underlying unity in tertiary structure that belies an extraordinary diversity in primary structures, biochemical properties, and physiological functions. Phylogenetic reconstructions can reveal which of these functions represent novel, lineage-specific innovations, and which represent ancestral functions that are shared with homologous globin proteins in other eukaryotes and even prokaryotes. To date, our understanding of globin diversity in deuterostomes has been hindered by a dearth of genomic sequence data from the Ambulacraria (echinoderms + hemichordates), the sister group of chordates, and the phylum Xenacoelomorpha, which includes xenoturbellids, acoelomorphs, and nemertodermatids. Here, we report the results of a phylogenetic and comparative genomic analysis of the globin gene repertoire of deuterostomes. We first characterized the globin genes of the acorn worm, Saccoglossus kowalevskii, a representative of the phylum Hemichordata. We then integrated genomic sequence data from the acorn worm into a comprehensive analysis of conserved synteny and phylogenetic relationships among globin genes from representatives of the eight lineages that comprise the superphylum Deuterostomia. The primary aims were 1) to unravel the evolutionary history of the globin gene superfamily in deuterostomes and 2) to use the estimated phylogeny to gain insights into the functional evolution of deuterostome globins. Results of our analyses indicate that the deuterostome common ancestor possessed a repertoire of at least four distinct globin paralogs and that different subsets of these ancestral genes have been retained in each of the descendant organismal lineages. In each major deuterostome group, a different subset of ancestral precursor genes underwent lineage-specific expansions of functional diversity through repeated rounds of gene duplication and divergence. By integrating results of the phylogenetic analysis with available functional data, we discovered that circulating oxygen-transport hemoglobins evolved independently in several deuterostome lineages and that intracellular nerve globins evolved independently in chordates and acoelomorph worms.     

Genomic structure and phylogeny of the plant pathogen Ralstonia solanacearum inferred from gene distribution analysis

In the present study, we investigated the gene distribution among strains of the highly polymorphic plant pathogenic beta-proteobacterium Ralstonia solanacearum, paying particular attention to the status of known or candidate pathogenicity genes. Based on the use of comparative genomic hybridization on a pangenomic microarray for the GMI1000 reference strain, we have defined the conditions that allowed comparison of the repertoires of genes among a collection of 18 strains that are representative of the biodiversity of the R. solanacearum species. This identified a list of 2,690 core genes present in all tested strains. As a corollary, a list of 2,338 variable genes within the R. solanacearum species has been defined. The hierarchical clustering based on the distribution of variable genes is fully consistent with the phylotype classification that was previously defined from the nucleotide sequence analysis of four genes. The presence of numerous pathogenicity-related genes in the core genome indicates that R. solanacearum is an ancestral pathogen. The results establish the long coevolution of the two replicons that constitute the bacterial genome. We also demonstrate the clustering of variable genes in genomic islands. Most genomic islands are included in regions with an alternative codon usage, suggesting that they originate from acquisition of foreign genes through lateral gene transfers. Other genomic islands correspond to genes that have the same base composition as core genes, suggesting that they either might be ancestral genes lost by deletion in certain strains or might originate from horizontal gene transfers.     

Contiguous genomic DNA sequence comprising the 19-kD zein gene family from maize

A new approach has been undertaken to analyze the sequences and linear organization of the 19-kD zein genes in maize (Zea mays). A high-coverage, large-insert genomic library of the inbred line B73 based on bacterial artificial chromosomes was used to isolate a redundant set of clones containing members of the 19-kD zein gene family, which previously had been estimated to consist of 50 members. The redundant set of clones was used to create bins of overlapping clones that represented five distinct genomic regions. Representative clones containing the entire set of 19-kD zein genes were chosen from each region and sequenced. Seven bacterial artificial chromosome clones yielded 1,160 kb of genomic DNA. Three of them formed a contiguous sequence of 478 kb, the longest contiguous sequenced region of the maize genome. Altogether, these DNA sequences provide the linear organization of 25 19-kD zein genes, one-half the number previously estimated. It is suggested that the difference is because of haplotypes exhibiting different degrees of gene amplification in the zein multigene family. About one-half the genes present in B73 appear to be expressed. Because some active genes have only been duplicated recently, they are so conserved in their sequence that previous cDNA sequence analysis resulted in "unigenes" that were actually derived from different gene copies. This analysis also shows that the 22- and 19-kD zein gene families shared a common ancestor. Although both ancestral genes had the same incremental gene amplification, the 19-kD zein branch exhibited a greater degree of far-distance gene translocations than the 22-kD zein gene family.     

Phylogenetic relationships among early-diverging eudicots based on four genes: were the eudicots ancestrally woody?

Based on analyses of combined data sets of three genes (18S rDNA, rbcL, and atpB), phylogenetic relationships among the early-diverging eudicot lineages (Ranunculales, Proteales, Trochodendraceae, Sabiaceae, and Buxaceae) remain unclear, as are relationships within Ranunculales, especially the placement of Eupteleaceae. To clarify relationships among these early-diverging eudicot lineages, we added entire sequences of 26S rDNA to the existing three-gene data set. In the combined analyses of four genes based on parsimony, ML, and Bayesian analysis, Ranunculales are strongly supported as a clade and are sister to other eudicots. Proteales appear as sister to the remaining eudicots, which are weakly (59%) supported as a clade. Relationships among Trochodendraceae, Buxaceae (including Didymeles), Sabiaceae, and Proteales remain unclear. Within Ranunculales, Eupteleaceae are sister to all other Ranunculales, with bootstrap support of 70% in parsimony analysis and with posterior probability of 1.00 in Bayesian analysis. Our character reconstructions indicate that the woody habit is ancestral, not only for the basal angiosperms, but also for the eudicots. Furthermore, Ranunculales may not be ancestrally herbaceous, as long maintained. The woody habit appears to have been ancestral for several major clades of eudicots, including Caryophyllales, and asterids.     

A phylogenetic hypothesis for species of the genus Taenia (Eucestoda : Taeniidae)

Cladistic analysis of a numerical data matrix describing 27 characters for species of Taenia resulted in 4 most parsimonious phylogenetic trees (174 steps; consistency index = 0.28; homoplasy index = 0.72; retention index = 0.48). Monophyly for Taenia is diagnosed by the metacestode that is either a cysticercus or a form derived from a bladder-like larva; no other unequivocal synapomorphies are evident. Tree structure provides no support for recognition of a diversity of tribes or genera within the Taeniinae: Fimbriotaeniini and Taeniini have no phylogenetic basis. Hydatigera, Fimbriotaenia, Fossor, Monordotaenia, Multiceps, Taeniarhynchus, Tetratirotaenia must be subsumed within Taenia as synonyms. Taenia saginata and Taenia asiatica are sister species and distantly related to Taenia solium. Cospeciation with respect to carnivorous definitive hosts and Taenia appears to be limited. Although felids are putative ancestral hosts, contemporary associations appear to have resulted from extensive host-switching among felids, canids, hyaenids, and others. In contrast, relationships with herbivorous intermediate hosts are indicative of more pervasive coevolution; rodents as intermediate hosts are postulated as ancestral for the Taeniidae, Taenia + Echinococcus. Patterns appear consistent with rapid shifts between phylogenetically unrelated carnivores but among those that historically exploited a common prey resource within communities in specific biogeographic regions.     

Horizontal gene transfer in the molecular evolution of mannose PTS transporters

The phosphoenolpyruvate:carbohydrate phosphotransferase system (PTS) represents hitherto the only example of group translocation transport systems. PTS transporters are exclusively found in bacteria and can be grouped on the basis of sequence and structure into six classes. We have analyzed the evolution of mannose-class PTS transporters. These transporters have a limited distribution among bacteria being mostly harbored by species associated to animals. The results obtained indicate that these genes have undergone a complex evolutionary history, including extensive horizontal gene transfer events, duplications, and nonorthologous displacements. The phylogenetic analysis revealed an early diversification to specialize in different transport capabilities, but these events have also occurred relatively recently. In addition, these transporters can be further divided into seven groups and this division correlates with their transport capabilities. Finally, the consideration of the genomic context allowed us to propose putative functional roles for some uncharacterized PTS transporters. The functional role and distribution of mannose-class PTS transporters suggest that their expansion may have played a significant role in the establishment of symbiotic relationships between animals and some bacteria.