Phylogenetic relationships of ascomycetes and basidiomycetes based on comparative genomics analysis

The relationship between ascomycetes and basidiomycetes, the two main phyla of non-flagellated fungi, has rarely been investigated. In this study, we performed a comparative genomics analysis of genome sequences of 55 ascomycetes and 26 basidiomycetes species and detected 81 universal markers, 875 homologous genes and a conserved contig in the glucose-regulated protein gene. In dendrograms based on simple sequence repeat markers and homologous genes, ascomycetes and basidiomycetes formed distinct clusters, with each set of taxa having a high coefficient of relatedness. Ascomycetes and basidiomycetes also constituted distinct groups in a phylogenetic tree based on a conserved contig in the glucose-regulated protein gene. These results provide evidence that basidiomycetes may be derived from ascomycetes but are definitely genetically differentiated at the genomic level. The phylogenetic relationships of ascomycetes and basidiomycetes uncovered in this study provide new insights for future research related to fungal classification and evolution.     

Estimation of bacterial species phylogeny through oligonucleotide frequency distances

Classification of bacteria is mainly based on sequence comparisons of certain homologous genes such as 16S rRNA. Recently there are challenges to classify bacteria using oligonucleotide frequency pattern of nonhomologous sequences. However, the evolutionary significance of oligonucleotides longer than tetra-nucleotide is not studied well. We performed phylogenetic analysis by using the Euclidean distances calculated from the di to deca-nucleotide frequencies in bacterial genomes, and compared these oligonucleotide frequency-based tree topologies with those for 16S rRNA gene and concatenated seven genes. When oligonucleotide frequency-based trees were constructed for bacterial species with similar GC content, their topologies at genus and family level were congruent with those based on homologous genes. Our results suggest that oligonucleotide frequency is useful not only for classification of bacteria, but also for estimation of their phylogenetic relationships for closely related species.     

Free from mitochondrial DNA: Nuclear genes and the inference of species trees among closely related darter lineages (Teleostei: Percidae: Etheostomatinae)

Investigations into the phylogenetics of closely related animal species are dominated by the use of mitochondrial DNA (mtDNA) sequence data. However, the near-ubiquitous use of mtDNA to infer phylogeny among closely related animal lineages is tempered by an increasing number of studies that document high rates of transfer of mtDNA genomes among closely related species through hybridization, leading to substantial discordance between phylogenies inferred from mtDNA and nuclear gene sequences. In addition, the recent development of methods that simultaneously infer a species phylogeny and estimate divergence times, while accounting for incongruence among individual gene trees, has ushered in a new era in the investigation of phylogeny among closely related species. In this study we assess if DNA sequence data sampled from a modest number of nuclear genes can resolve relationships of a species-rich clade of North American freshwater teleost fishes, the darters. We articulate and expand on a recently introduced method to infer a time-calibrated multi-species coalescent phylogeny using the computer program ^*BEAST. Our analyses result in well-resolved and strongly supported time-calibrated darter species tree. Contrary to the expectation that mtDNA will provide greater phylogenetic resolution than nuclear gene data; the darter species tree inferred exclusively from nuclear genes exhibits a higher frequency of strongly supported nodes than the mtDNA time-calibrated gene tree.     

Mitochondrial genomes and divergence times of crocodile newts: inter-islands distribution of Echinotriton andersoni and the origin of a unique repetitive sequence found in Tylototriton mt genomes

Crocodile newts, which constitute the genera Echinotriton and Tylototriton, are known as living fossils, and these genera comprise many endangered species. To identify mitochondrial (mt) genes suitable for future population genetic analyses for endangered taxa, we determined the complete nucleotide sequences of the mt genomes of the Japanese crocodile newt Echinotriton andersoni and Himalayan crocodile newt Tylototriton verrucosus. Although the control region (CR) is known as the most variable mtDNA region in many animal taxa, the CRs of crocodile newts are highly conservative. Rather, the genes of NADH dehydrogenase subunits and ATPase subunit 6 were found to have high sequence divergences and to be usable for population genetics studies. To estimate the inter-population divergence ages of E. andersoni endemic to the Ryukyu Islands, we performed molecular dating analysis using whole and partial mt genomic data. The estimated divergence ages of the inter-island individuals are older than the paleogeographic segmentation ages of the islands, suggesting that the lineage splits of E. andersoni populations were not caused by vicariant events. Our phylogenetic analysis with partial mt sequence data also suggests the existence of at least two more undescribed species in the genus Tylototriton. We also found unusual repeat sequences containing the 3' region of cytochrome apoenzyme b gene, whole tRNA-Thr gene, and a noncoding region (the T-P noncoding region characteristic in caudate mtDNAs) from T. verrucosus mtDNA. Similar repeat sequences were found in two other Tylototriton species. The Tylototriton taxa with the repeats become a monophyletic group, indicating a single origin of the repeat sequences. The intra-and inter-specific comparisons of the repeat sequences suggest the occurrences of homologous recombination-based concerted evolution among the repeat sequences.     

Trees of trees: an approach to comparing multiple alternative phylogenies

Phylogenetic analysis very commonly produces several alternative trees for a given fixed set of taxa. For example, different sets of orthologous genes may be analyzed, or the analysis may sample from a distribution of probable trees. This article describes an approach to comparing and visualizing multiple alternative phylogenies via the idea of a "tree of trees" or "meta-tree." A meta-tree clusters phylogenies with similar topologies together in the same way that a phylogeny clusters species with similar DNA sequences. Leaf nodes on a meta-tree correspond to the original set of phylogenies given by some analysis, whereas interior nodes correspond to certain consensus topologies. The construction of meta-trees is motivated by analogy with construction of a most parsimonious tree for DNA data, but instead of using DNA letters, in a meta-tree the characters are partitions or splits of the set of taxa. An efficient algorithm for meta-tree construction is described that makes use of a known relationship between the majority consensus and parsimony in terms of gain and loss of splits. To illustrate these ideas meta-trees are constructed for two datasets: a set of gene trees for species of yeast and trees from a bootstrap analysis of a set of gene trees in ray-finned fish. A software tool for constructing meta-trees and comparing alternative phylogenies is available online, and the source code can be obtained from the author.     

Molecular and evolutionary relationships among enteric bacteria

Classification of bacterial species into genera has traditionally relied upon variation in phenotypic characteristics. However, these phenotypes often have a multifactorial genetic basis, making unambiguous taxonomic placement of new species difficult. By designing evolutionarily conserved oligonucleotide primers, it is possible to amplify homologous regions of genes in diverse taxa using the polymerase chain reaction and determine their nucleotide sequences. We have constructed a phylogeny of some enteric bacteria, including five species classified as members of the genus Escherichia, based on nucleotide sequence variation at the loci encoding glyceraldehyde-3-phosphate dehydrogenase and outer membrane protein 3A, and compared this genealogy with the relationships inferred by biotyping. The DNA sequences of these genes defined congruent and robust phylogenetic trees indicating that they are an accurate reflection of the evolutionary history of the bacterial species. The five species of Escherichia were found to be distantly related and, contrary to their placement in the same genus, do not form a monophyletic group. These data provide a framework which allows the relationships of additional species of enteric bacteria to be inferred. These procedures have general applicability for analysis of the classification, evolution, and epidemiology of bacterial taxa.     

Using Homolog Groups to Create a Whole-Genomic Tree of Free-Living Organisms: An Update

Genomic trees have been constructed based on the presence and absence of families of protein-encoding genes observed in 27 complete genomes, including genomes of 15 free-living organisms. This method does not rely on the identification of suspected orthologs in each genome, nor the specific alignment used to compare gene sequences because the protein-encoding gene families are formed by grouping any protein with a pairwise similarity score greater than a preset value. Because of this all inclusive grouping, this method is resilient to some effects of lateral gene transfer because transfers of genes are masked when the recipient genome already has a homolog (not necessarily an ortholog) of the incoming gene. Of 71 genes suspected to have been laterally transferred to the genome of Aeropyrum pernix, only approximately 7 to 15 represent genes where a lateral gene transfer appears to have generated homoplasy in our character dataset. The genomic tree of the 15 free-living taxa includes six different bacterial orders, six different archaeal orders, and two different eukaryotic kingdoms. The results are remarkably similar to results obtained by analysis of rRNA. Inclusion of the other 12 genomes resulted in a tree only broadly similar to that suggested by rRNA with at least some of the differences due to artifacts caused by the small genome size of many of these species. Very small genomes, such as those of the two Mycoplasma genomes included, fall to the base of the Bacterial domain, a result expected due to the substantial gene loss inherent to these lineages. Finally, artificial ``partial genomes' were generated by randomly selecting ORFs from the complete genomes in order to test our ability to recover the tree generated by the whole genome sequences when only partial data are available. The results indicated that partial genomic data, when sampled randomly, could robustly recover the tree generated by the whole genome sequences. Received: 30 May 2001 / Accepted: 10 October 2001     

Genome-scale phylogenetics: inferring the plant tree of life from 18,896 gene trees

Phylogenetic analyses using genome-scale data sets must confront incongruence among gene trees, which in plants is exacerbated by frequent gene duplications and losses. Gene tree parsimony (GTP) is a phylogenetic optimization criterion in which a species tree that minimizes the number of gene duplications induced among a set of gene trees is selected. The run time performance of previous implementations has limited its use on large-scale data sets. We used new software that incorporates recent algorithmic advances to examine the performance of GTP on a plant data set consisting of 18,896 gene trees containing 510,922 protein sequences from 136 plant taxa (giving a combined alignment length of >2.9 million characters). The relationships inferred from the GTP analysis were largely consistent with previous large-scale studies of backbone plant phylogeny and resolved some controversial nodes. The placement of taxa that were present in few gene trees generally varied the most among GTP bootstrap replicates. Excluding these taxa either before or after the GTP analysis revealed high levels of phylogenetic support across plants. The analyses supported magnoliids sister to a eudicot + monocot clade and did not support the eurosid I and II clades. This study presents a nuclear genomic perspective on the broad-scale phylogenic relationships among plants, and it demonstrates that nuclear genes with a history of duplication and loss can be phylogenetically informative for resolving the plant tree of life.     

Reciprocal illumination in the gene content tree of life

Phylogenies based on gene content rely on statements of primary homology to characterize gene presence or absence. These statements (hypotheses) are usually determined by techniques based on threshold similarity or distance measurements between genes. This fundamental but problematic step can be examined by evaluating each homology hypothesis by the extent to which it is corroborated by the rest of the data. Here we test the effects of varying the stringency for making primary homology statements using a range of similarity (e-value) cutoffs in 166 fully sequenced and annotated genomes spanning the tree of life. By evaluating each resulting data set with tree-based measurements of character consistency and information content, we find a set of homology statements that optimizes overall corroboration. The resulting data set produces well-resolved and well-supported trees of life and greatly ameliorates previously noted inconsistencies such as the misclassification of small genomes. The method presented here, which can be used to test any technique for recognizing primary homology, provides an objective framework for evaluating phylogenetic hypotheses and data sets for the tree of life. It also can serve as a technique for identifying well-corroborated sets of homologous genes for functional genomic applications.     

Inferring angiosperm phylogeny from EST data with widespread gene duplication

Background

Most studies inferring species phylogenies use sequences from single copy genes or sets of orthologs culled from gene families. For taxa such as plants, with very high levels of gene duplication in their nuclear genomes, this has limited the exploitation of nuclear sequences for phylogenetic studies, such as those available in large EST libraries. One rarely used method of inference, gene tree parsimony, can infer species trees from gene families undergoing duplication and loss, but its performance has not been evaluated at a phylogenomic scale for EST data in plants.

Results

A gene tree parsimony analysis based on EST data was undertaken for six angiosperm model species and Pinus, an outgroup. Although a large fraction of the tentative consensus sequences obtained from the TIGR database of ESTs was assembled into homologous clusters too small to be phylogenetically informative, some 557 clusters contained promising levels of information. Based on maximum likelihood estimates of the gene trees obtained from these clusters, gene tree parsimony correctly inferred the accepted species tree with strong statistical support. A slight variant of this species tree was obtained when maximum parsimony was used to infer the individual gene trees instead.

Conclusion

Despite the complexity of the EST data and the relatively small fraction eventually used in inferring a species tree, the gene tree parsimony method performed well in the face of very high apparent rates of duplication.

     

Use of whole genome sequence data to infer baculovirus phylogeny

Several phylogenetic methods based on whole genome sequence data were evaluated using data from nine complete baculovirus genomes. The utility of three independent character sets was assessed. The first data set comprised the sequences of the 63 genes common to these viruses. The second set of characters was based on gene order, and phylogenies were inferred using both breakpoint distance analysis and a novel method developed here, termed neighbor pair analysis. The third set recorded gene content by scoring gene presence or absence in each genome. All three data sets yielded phylogenies supporting the separation of the Nucleopolyhedrovirus (NPV) and Granulovirus (GV) genera, the division of the NPVs into groups I and II, and species relationships within group I NPVs. Generation of phylogenies based on the combined sequences of all 63 shared genes proved to be the most effective approach to resolving the relationships among the group II NPVs and the GVs. The history of gene acquisitions and losses that have accompanied baculovirus diversification was visualized by mapping the gene content data onto the phylogenetic tree. This analysis highlighted the fluid nature of baculovirus genomes, with evidence of frequent genome rearrangements and multiple gene content changes during their evolution. Of more than 416 genes identified in the genomes analyzed, only 63 are present in all nine genomes, and 200 genes are found only in a single genome. Despite this fluidity, the whole genome-based methods we describe are sufficiently powerful to recover the underlying phylogeny of the viruses.     

Conditionally Neutral Phylogenetic Markers of Major Taxa: A New Aspect of the Evolution of Macromolecules

The current phase of molecular phylogenetics can be named the 18S rRNA gene era, which is now approaching the end. To date, almost all phyla of metazoans and many taxa of protists are represented in databases of 18S rRNA gene sequences. The elements of the phylogenetic tree of Metazoa inferred from 18S rRNA genes are characterized by unequal validity: some of them seem to be well grounded; others are not adequately supported, and probably will be revised later. The validity of phylogenetic reconstruction is influenced by two main factors: (1) erroneous grouping of long branches that occur because of abnormally high evolution rate; (2) deficit of phylogenetically informative characters. A method for overcoming these difficulties is suggested in addition to known tools: using phylogenetic markers that are stable within individual taxa and evolve by punctuated equilibrium. These markers are least influenced by the convergence caused by a high evolution rate of the entire gene. The nature of these markers of ancient taxa, paradoxical from the perspective of neutral evolution, is discussed, as well as their importance for establishing monophyly of both new large-scale taxonomic groups of invertebrates (Bilateria + Rhombozoa + Orthonectida + Myxozoa + Cnidaria + Placozoa and Echinodermata + Hemichordata) and some major taxa of Nematoda.     

Automated ortholog inference from phylogenetic trees and calculation of orthology reliability

MOTIVATION: Orthologous proteins in different species are likely to have similar biochemical function and biological role. When annotating a newly sequenced genome by sequence homology, the most precise and reliable functional information can thus be derived from orthologs in other species. A standard method of finding orthologs is to compare the sequence tree with the species tree. However, since the topology of phylogenetic tree is not always reliable one might get incorrect assignments. RESULTS: Here we present a novel method that resolves this problem by analyzing a set of bootstrap trees instead of the optimal tree. The frequency of orthology assignments in the bootstrap trees can be interpreted as a support value for the possible orthology of the sequences. Our method is efficient enough to analyze data in the scale of whole genomes. It is implemented in Java and calculates orthology support levels for all pairwise combinations of homologous sequences of two species. The method was tested on simulated datasets and on real data of homologous proteins.     

Conditionally neutral phylogenetic markers of major taxa: a new aspect of the evolution of macromolecules

The current phase of molecular phylogenetics can be named the 18S rRNA gene era, which is now approaching the end. To date, almost all phyla of metazoans and many taxa of protists are represented in databases of 18S rRNA gene sequences. The elements of the phylogenetic tree of Metazoa inferred from 18S rRNA genes are characterized by unequal validity: some of them seem to be well grounded; others are not adequately supported, and probably will be revised later. The validity of phylogenetic reconstruction is influenced by two main factors: (1) erroneous grouping of long branches that occur because of abnormally high evolution rate; (2) deficit of phylogenetically informative characters. A method for overcoming these difficulties is suggested in addition to known tools: using phylogenetic markers that are stable within individual taxa and evolve by punctuated equilibrium. These markers are least influenced by the convergence caused by a high evolution rate of the entire gene. The nature of these markers of ancient taxa, paradoxical from the perspective of neutral evolution, is discussed, as well as their importance for establishing monophyly of both new large-scale taxonomic groups of invertebrates (Bilateria + Rhombozoa + Orthonectida + Myxozoa + Cnidaria + Placozoa and Echinodermata + Hemichordata) and some major taxa of Nematoda.     

Genotypic differentiation of twelve Clostridium species by polymorphism analysis of the triosephosphate isomerase (tpi) gene

Housekeeping genes encoding metabolic enzymes may provide alternative markers to 16S ribosomal DNA (rDNA) for genotypic and phylogenetic characterization of bacterial species. We have developed a PCR-restriction fragment length polymorphism (PCR-RFLP) assay, targeting the triosephosphate isomerase (tpi) gene, which allows the differentiation of twelve pathogenic Clostridium species. Degenerate primers constructed from alignments of tpi sequences of various gram-positive bacteria allowed the amplification of a 501 bp target region in the twelve Clostridium type strains. A phylogenetic tree constructed from the nucleotidic sequences of these tpi amplicons was well correlated with that inferred from analysis of 16S rDNA gene sequences. The analysis of tpi sequences revealed restriction sites of enzyme AluI that could be species-specific. Indeed, AluI digestion of amplicons from the twelve type strains provided distinct restriction patterns. A total of 127 strains (three to sixteen strains for each species) was further analyzed by PCR-RFLP of the tpi gene, and confirmed that each species could be characterized by one to three restriction types (RTs). The differences between RTs within species could be explained by point mutations in AluI restriction sites of the tpi sequences. PCR-restriction analysis of the tpi gene offers an accurate tool for species identification within the genus Clostridium, and provides an alternative marker to 16S rDNA for phylogenetic analyses.     

N. I. Vavilov's law of genetic homology in modern genetics

The paper considers the place in modern genetics and the significance of the main N.I. Vavilov's generalization--his law of homologous series in variation. Recent analysis of amino acid sequences of gene products (proteins) and especially of nucleotide sequences of genes shows the high degree of molecular homology between genes of closely related species and the retention of the homology through the course of evolution. The study of homologous genes disposition in chromosomes shows conservation of the similar orders of genes in many organisms, especially in mammals. Thus, the law of genetic homology has been confirmed by modern genetic researchers. It is a foundation-stone of comparative genetics--new and rapid developing branch of genetics which involves studies on similarity and differences in heredity and variation in organisms of different taxa.     

Signature genes as a phylogenomic tool

Gene content has been shown to contain a strong phylogenetic signal, yet its usage for phylogenetic questions is hampered by horizontal gene transfer and parallel gene loss and until now required completely sequenced genomes. Here, we introduce an approach that allows the phylogenetic signal in gene content to be applied to any set of sequences, using signature genes for phylogenetic classification. The hundreds of publicly available genomes allow us to identify signature genes at various taxonomic depths, and we show how the presence of signature genes in an unspecified sample can be used to characterize its taxonomic composition. We identify 8,362 signature genes specific for 112 prokaryotic taxa. We show that these signature genes can be used to address phylogenetic questions on the basis of gene content in cases where classic gene content or sequence analyses provide an ambiguous answer, such as for Nanoarchaeum equitans, and even in cases where complete genomes are not available, such as for metagenomics data. Cross-validation experiments leaving out up to 30% of the species show that approximately 92% of the signature genes correctly place the species in a related clade. Analyses of metagenomics data sets with the signature gene approach are in good agreement with the previously reported species distributions based on phylogenetic analysis of marker genes. Summarizing, signature genes can complement traditional sequence-based methods in addressing taxonomic questions.     

Delimiting species in recent radiations

Despite considerable effort from the systematics community, delimiting species boundaries in recent radiations remains a daunting challenge. We argue that genealogical approaches, although sometimes useful, may not solve this important problem, because recently derived species often have not had sufficient time to achieve monophyly. Instead, we suggest that population genetic approaches that rely on large sets of informative markers like single nucleotide polymorphisms (SNPs) provide an alternative framework for delimiting very recently derived species. We address two major challenges in applying such markers to species delimitation: discovering markers in nonmodel systems and using them to delimit recently derived species. Using turtles as a test case, we explore the utility of a single, relatively low-coverage genomic resource as an aid in gene and marker discovery. We exploit an end-sequenced bacterial artificial chromosome (BAC) library from an individual painted turtle (Chrysemys picta) and outline a novel protocol that efficiently identifies primer pairs that amplify homologous sequences across the tree of living turtles. Preliminary data using this library to discover SNPs in Emydura macquarii, a species that diverged from C. picta approximately 210 million years ago, indicate that sequences identified from the Chrysemys BAC library provide useful SNPs even in this very distantly related taxon. Several recent methods in wide use in the population genetics literature allow one to discover potential species, or test existing species hypotheses, with SNP data and may be particularly informative for very recently derived species. As BAC and other genomic resources become increasingly available for scattered taxa across the tree of life, we are optimistic that these resources will provide abundant, inexpensive markers that will help delimit boundaries in problematic, recent species radiations.     

Reconciling gene and genome duplication events: using multiple nuclear gene families to infer the phylogeny of the aquatic plant family Pontederiaceae

Most plant phylogenetic inference has used DNA sequence data from the plastid genome. This genome represents a single genealogical sample with no recombination among genes, potentially limiting the resolution of evolutionary relationships in some contexts. In contrast, nuclear DNA is inherently more difficult to employ for phylogeny reconstruction because major mutational events in the genome, including polyploidization, gene duplication, and gene extinction can result in homologous gene copies that are difficult to identify as orthologs or paralogs. Gene tree parsimony (GTP) can be used to infer the rooted species tree by fitting gene genealogies to species trees while simultaneously minimizing the estimated number of duplications needed to reconcile conflicts among them. Here, we use GTP for five nuclear gene families and a previously published plastid data set to reconstruct the phylogenetic backbone of the aquatic plant family Pontederiaceae. Plastid-based phylogenetic studies strongly supported extensive paraphyly of Eichhornia (one of the four major genera) but also depicted considerable ambiguity concerning the true root placement for the family. Our results indicate that species trees inferred from the nuclear genes (alone and in combination with the plastid data) are highly congruent with gene trees inferred from plastid data alone. Consideration of optimal and suboptimal gene tree reconciliations place the root of the family at (or near) a branch leading to the rare and locally restricted E. meyeri. We also explore methods to incorporate uncertainty in individual gene trees during reconciliation by considering their individual bootstrap profiles and relate inferred excesses of gene duplication events on individual branches to whole-genome duplication events inferred for the same branches. Our study improves understanding of the phylogenetic history of Pontederiaceae and also demonstrates the utility of GTP for phylogenetic analysis.     

From Arabidopsis thaliana to Brassica napus: development of amplified consensus genetic markers (ACGM) for construction of a gene map

The evolution of genomes can be studied by comparing maps of homologous genes which show changes in nucleic acid sequences and chromosome rearrangements. In this study, we developed a set of 32 amplified consensus gene markers (ACGMs) that amplified gene sequences from Arabidopsis thaliana and Brassica napus. Our methodology, based on PCR, facilitated the rapid sequencing of homologous genes from various species of the same phylogenetic family and the detection of intragenic polymorphism. We found that such polymorphism principally concerned intron sequences and we used it to attribute a Brassica oleracea or Brassica rapa origin to the B. napus sequences and to map 43 rapeseed genes. We confirm that the genetic position of homologous genes varied between B. napus and A. thaliana. ACGMs are a useful tool for genome evolution studies and for the further development of single nucleotide polymorphism suitable for use in genetic mapping and genetic diversity analyses.