Assessing the performance of single-copy genes for recovering robust phylogenies

Phylogenies involving nonmodel species are based on a few genes, mostly chosen following historical or practical criteria. Because gene trees are sometimes incongruent with species trees, the resulting phylogenies may not accurately reflect the evolutionary relationships among species. The increase in availability of genome sequences now provides large numbers of genes that could be used for building phylogenies. However, for practical reasons only a few genes can be sequenced for a wide range of species. Here we asked whether we can identify a few genes, among the single-copy genes common to most fungal genomes, that are sufficient for recovering accurate and well-supported phylogenies. Fungi represent a model group for phylogenomics because many complete fungal genomes are available. An automated procedure was developed to extract single-copy orthologous genes from complete fungal genomes using a Markov Clustering Algorithm (Tribe-MCL). Using 21 complete, publicly available fungal genomes with reliable protein predictions, 246 single-copy orthologous gene clusters were identified. We inferred the maximum likelihood trees using the individual orthologous sequences and constructed a reference tree from concatenated protein alignments. The topologies of the individual gene trees were compared to that of the reference tree using three different methods. The performance of individual genes in recovering the reference tree was highly variable. Gene size and the number of variable sites were highly correlated and significantly affected the performance of the genes, but the average substitution rate did not. Two genes recovered exactly the same topology as the reference tree, and when concatenated provided high bootstrap values. The genes typically used for fungal phylogenies did not perform well, which suggests that current fungal phylogenies based on these genes may not accurately reflect the evolutionary relationships among species. Analyses on subsets of species showed that the phylogenetic performance did not seem to depend strongly on the sample. We expect that the best-performing genes identified here will be very useful for phylogenetic studies of fungi, at least at a large taxonomic scale. Furthermore, we compare the method developed here for finding genes for building robust phylogenies with previous ones and we advocate that our method could be applied to other groups of organisms when more complete genomes are available.     

A phylogenomics approach for selecting robust sets of phylogenetic markers

Reconstructing the evolutionary relationships of species is a major goal in biology. Despite the increasing number of completely sequenced genomes, a large number of phylogenetic projects rely on targeted sequencing and analysis of a relatively small sample of marker genes. The selection of these phylogenetic markers should ideally be based on accurate predictions of their combined, rather than individual, potential to accurately resolve the phylogeny of interest. Here we present and validate a new phylogenomics strategy to efficiently select a minimal set of stable markers able to reconstruct the underlying species phylogeny. In contrast to previous approaches, our methodology does not only rely on the ability of individual genes to reconstruct a known phylogeny, but it also explores the combined power of sets of concatenated genes to accurately infer phylogenetic relationships of species not previously analyzed. We applied our approach to two broad sets of cyanobacterial and ascomycetous fungal species, and provide two minimal sets of six and four genes, respectively, necessary to fully resolve the target phylogenies. This approach paves the way for the informed selection of phylogenetic markers in the effort of reconstructing the tree of life.     

A fungal phylogeny based on 42 complete genomes derived from supertree and combined gene analysis

Background  

To date, most fungal phylogenies have been derived from single gene comparisons, or from concatenated alignments of a small number of genes. The increase in fungal genome sequencing presents an opportunity to reconstruct evolutionary events using entire genomes. As a tool for future comparative, phylogenomic and phylogenetic studies, we used both supertrees and concatenated alignments to infer relationships between 42 species of fungi for which complete genome sequences are available.     

New molecular markers for fungal phylogenetics: two genes for species-level systematics in the Sordariomycetes (Ascomycota)

Although significant progress has been made resolving deep branches of the fungal tree of life, many fungal systematists are interested in species-level questions to both define species and assess fungal biodiversity. Fungal genome sequences are a useful resource to systematic biologists for developing new phylogenetic markers that better represent the whole genome. Here we report primers for two newly identified single-copy protein-coding genes, FG1093 and MS204, for use with ascomycetes. Although fungi were the focus of this study, this methodological approach could be easily applied to marker development for studies of other organisms. The tests used here to assess phylogenetic informativeness are computationally rapid, require only rudimentary datasets to evaluate existing or newly developed markers, and can be applied to other non-model organisms to assist in experimental design of phylogenetic studies. Phylogenetic utility of the markers was tested in two genera, Gnomoniopsis and Ophiognomonia (Gnomoniaceae, Diaporthales). The phylogenetic performance of β-tubulin, ITS, and tef-1α was compared with FG1093 and MS204. Phylogenies inferred from FG1093 and MS204 were largely in agreement with β-tubulin, ITS, and tef-1α although some topological conflict was observed. Resolution and support for branches differed based on the combination of markers used for each genus. Based on two independent tests of phylogenetic performance, FG1093 and MS204 were determined to be equal to or better than β-tubulin, ITS, and tef-1α in resolving species relationships. Differences were found in site-specific rate of evolution in all five markers. In addition, isolates from 15 orders and 22 families of Ascomycota were screened using primers for FG1093 and MS204 to demonstrate primer utility across a wide diversity of ascomycetes. The primer sets for the newly identified genes FG1093 and MS204 and methods used to develop them are useful additions to the ascomycete systematists' toolbox.     

Phylogeny of the glomeromycota (arbuscular mycorrhizal fungi): recent developments and new gene markers

The fungal symbionts of arbuscular mycorrhiza form a monophyletic group in the true Fungi, the phylum Glomeromycota. Fewer than 200 described species currently are included in this group. The only member of this clade known to form a different type of symbiosis is Geosiphon pyriformis, which associates with cyanobacteria. Because none of these fungi has been cultivated without their plant hosts or cyanobacterial partners, progress in obtaining multigene phylogenies has been slow and the nuclear-encoded ribosomal RNA genes have remained the only widely accessible molecular markers. rDNA phylogenies have revealed considerable polyphyly of some glomeromycotan genera that has been used to reassess taxonomic concepts. Environmental studies using phylogenetic methods for molecular identification have recovered an amazing diversity of unknown phylotypes, suggesting considerable cryptic species diversity. Protein gene sequences that have become available recently have challenged the rDNA-supported sister group relationship of the Glomeromycota with Asco/Basidiomycota. However the number of taxa analyzed with these new markers is still too small to provide a comprehensive picture of intraphylum relationships. We use nuclear-encoded rDNA and rpb1 protein gene sequences to reassess the phylogeny of the Glomeromycota and discuss possible implications.     

A new promising phylogenetic marker to study the diversity of fungal communities: The Glycoside Hydrolase 63 gene

In molecular ecology, the development of efficient molecular markers for fungi remains an important research domain. Nuclear ribosomal internal transcribed spacer (ITS) region was proposed as universal DNA barcode marker for fungi, but this marker was criticized for Indel‐induced alignment problems and its potential lack of phylogenetic resolution. Our main aim was to develop a new phylogenetic gene and a putative functional marker, from single‐copy gene, to describe fungal diversity. Thus, we developed a series of primers to amplify a polymorphic region of the Glycoside Hydrolase GH63 gene, encoding exo‐acting α‐glucosidases, in basidiomycetes. These primers were validated on 125 different fungal genomic DNAs, and GH63 amplification yield was compared with that of already published functional markers targeting genes coding for laccases, N‐acetylhexosaminidases, cellobiohydrolases and class II peroxidases. Specific amplicons were recovered for 95% of the fungal species tested, and GH63 amplification success was strikingly higher than rates obtained with other functional genes. We downloaded the GH63 sequences from 483 fungal genomes publicly available at the JGI mycocosm database. GH63 was present in 461 fungal genomes belonging to all phyla, except Microsporidia and Neocallimastigomycota divisions. Moreover, the phylogenetic trees built with both GH63 and Rpb1 protein sequences revealed that GH63 is also a promising phylogenetic marker. Finally, a very high proportion of GH63 proteins was predicted to be secreted. This molecular tool could be a new phylogenetic marker of fungal species as well as potential indicator of functional diversity of basidiomycetes fungal communities in term of secretory capacities.     

Complete Sequencing of Five Araliaceae Chloroplast Genomes and the Phylogenetic Implications

Background

The ginseng family (Araliaceae) includes a number of economically important plant species. Previously phylogenetic studies circumscribed three major clades within the core ginseng plant family, yet the internal relationships of each major group have been poorly resolved perhaps due to rapid radiation of these lineages. Recent studies have shown that phyogenomics based on chloroplast genomes provides a viable way to resolve complex relationships.

Methodology/Principal Findings

We report the complete nucleotide sequences of five Araliaceae chloroplast genomes using next-generation sequencing technology. The five chloroplast genomes are 156,333–156,459 bp in length including a pair of inverted repeats (25,551–26,108 bp) separated by the large single-copy (86,028–86,566 bp) and small single-copy (18,021–19,117 bp) regions. Each chloroplast genome contains the same 114 unique genes consisting of 30 transfer RNA genes, four ribosomal RNA genes, and 80 protein coding genes. Gene size, content, and order, AT content, and IR/SC boundary structure are similar among all Araliaceae chloroplast genomes. A total of 140 repeats were identified in the five chloroplast genomes with palindromic repeat as the most common type. Phylogenomic analyses using parsimony, likelihood, and Bayesian inference based on the complete chloroplast genomes strongly supported the monophyly of the Asian Palmate group and the Aralia-Panax group. Furthermore, the relationships among the sampled taxa within the Asian Palmate group were well resolved. Twenty-six DNA markers with the percentage of variable sites higher than 5% were identified, which may be useful for phylogenetic studies of Araliaceae.

Conclusion

The chloroplast genomes of Araliaceae are highly conserved in all aspects of genome features. The large-scale phylogenomic data based on the complete chloroplast DNA sequences is shown to be effective for the phylogenetic reconstruction of Araliaceae.     

Strong phylogenetic congruence between Tulasnella fungi and their associated Drakaeinae orchids

The study of congruency between phylogenies of interacting species can provide a powerful approach for understanding the evolutionary history of symbiotic associations. Orchid mycorrhizal fungi can survive independently of orchids making cospeciation unlikely, leading us to predict that any congruence would arise from host-switches to closely related fungal species. The Australasian orchid subtribe Drakaeinae is an iconic group of sexually deceptive orchids that consists of approximately 66 species. In this study, we investigated the evolutionary relationships between representatives of all six Drakaeinae orchid genera (39 species) and their mycorrhizal fungi. We used an exome capture dataset to generate the first well-resolved phylogeny of the Drakaeinae genera. A total of 10 closely related Tulasnella Operational Taxonomic Units (OTUs) and previously described species were associated with the Drakaeinae orchids. Three of them were shared among orchid genera, with each genus associating with 1–6 Tulasnella lineages. Cophylogenetic analyses show Drakaeinae orchids and their Tulasnella associates exhibit significant congruence (p < 0.001) in the topology of their phylogenetic trees. An event-based method also revealed significant congruence in Drakaeinae–Tulasnella relationships, with duplications (35), losses (25), and failure to diverge (9) the most frequent events, with minimal evidence for cospeciation (1) and host-switches (2). The high number of duplications suggests that the orchids speciate independently from the fungi, and the fungal species association of the ancestral orchid species is typically maintained in the daughter species. For the Drakaeinae–Tulasnella interaction, a pattern of phylogenetic niche conservatism rather than coevolution likely explains the observed phylogenetic congruency in orchid and fungal phylogenies. Given that many orchid genera are characterized by sharing of fungal species between closely related orchid species, we predict that these findings may apply to a wide range of orchid lineages.     

An empirical demonstration of using pentatricopeptide repeat (PPR) genes as plant phylogenetic tools: Phylogeny of Verbenaceae and the Verbena complex

The pentatricopeptide repeat (PPR) gene family, with hundreds of members in land plant genomes, has been recognized as a tremendous resource for plant phylogenetic studies based on publicly available genomic data from model organisms. However, whether this appealing nuclear gene marker system can be readily applied to non-model organisms remains questionable, particularly given the potential uncertainties in designing specific primers to only amplify the locus of interest from the sea of PPR genes. Here we demonstrate empirically the use of PPR genes in the family Verbenaceae and the Verbena complex. We also lay out a general scheme to design locus-specific primers to amplify and sequence PPR genes in non-model organisms. Intergeneric relationships within the family Verbenaceae were fully resolved with strong support. Relationships among the closely related genera within the Verbena complex and among some species groups within each genus were also well resolved, but resolution among very closely related species was limited. Our results suggest that PPR genes can be readily employed in non-model organisms. They may be best used to resolve relationships in a spectrum from among distantly related genera to among not-so-closely related congeneric species, but may have limited use among very closely related species.     

Evidence from beta-tubulin phylogeny that microsporidia evolved from within the fungi

Microsporidia are obligate intracellular parasites that were thought to be an ancient eukaryotic lineage based on molecular phylogenies using ribosomal RNA and translation elongation factors. However, this ancient origin of microsporidia has been contested recently, as several other molecular phylogenies suggest that microsporidia are closely related to fungi. Most of the protein trees that place microsporidia with fungi are not well sampled, however, and it is impossible to resolve whether microsporidia evolved from a fungus or from a protistan relative of fungi. We have sequenced beta-tubulins from 3 microsporidia, 4 chytrid fungi, and 12 zygomycete fungi, expanding the representation of beta-tubulin to include all four fungal divisions and a wide diversity of microsporidia. In phylogenetic trees including these new sequences, the overall topology of the fungal beta-tubulins generally matched the expected relationships among the four fungal divisions, although the zygomycetes were polyphyletic in some analyses. The microsporidia consistently fell within this fungal diversification, and not as a sister group to fungi. Overall, beta-tubulin phylogeny suggests that microsporidia evolved from a fungus sometime after the divergence of chytrids. We also found that chytrid alpha- and beta-tubulins are much less divergent than are tubulins from other fungi or microsporidia. In trees in which the only fungal representatives were the chytrids, microsporidia still branched with fungi (i.e., with chytrids), suggesting that the affiliation between microsporidian and fungal tubulins is not an artifact of long-branch attraction.     

A consistent phylogenetic backbone for the fungi

The kingdom of fungi provides model organisms for biotechnology, cell biology, genetics, and life sciences in general. Only when their phylogenetic relationships are stably resolved, can individual results from fungal research be integrated into a holistic picture of biology. However, and despite recent progress, many deep relationships within the fungi remain unclear. Here, we present the first phylogenomic study of an entire eukaryotic kingdom that uses a consistency criterion to strengthen phylogenetic conclusions. We reason that branches (splits) recovered with independent data and different tree reconstruction methods are likely to reflect true evolutionary relationships. Two complementary phylogenomic data sets based on 99 fungal genomes and 109 fungal expressed sequence tag (EST) sets analyzed with four different tree reconstruction methods shed light from different angles on the fungal tree of life. Eleven additional data sets address specifically the phylogenetic position of Blastocladiomycota, Ustilaginomycotina, and Dothideomycetes, respectively. The combined evidence from the resulting trees supports the deep-level stability of the fungal groups toward a comprehensive natural system of the fungi. In addition, our analysis reveals methodologically interesting aspects. Enrichment for EST encoded data-a common practice in phylogenomic analyses-introduces a strong bias toward slowly evolving and functionally correlated genes. Consequently, the generalization of phylogenomic data sets as collections of randomly selected genes cannot be taken for granted. A thorough characterization of the data to assess possible influences on the tree reconstruction should therefore become a standard in phylogenomic analyses.     

Assessment of phylogenomic and orthology approaches for phylogenetic inference

MOTIVATION: Phylogenomics integrates the vast amount of phylogenetic information contained in complete genome sequences, and is rapidly becoming the standard for reliably inferring species phylogenies. There are, however, fundamental differences between the ways in which phylogenomic approaches like gene content, superalignment, superdistance and supertree integrate the phylogenetic information from separate orthologous groups. Furthermore, they all depend on the method by which the orthologous groups are initially determined. Here, we systematically compare these four phylogenomic approaches, in parallel with three approaches for large-scale orthology determination: pairwise orthology, cluster orthology and tree-based orthology. RESULTS: Including various phylogenetic methods, we apply a total of 54 fully automated phylogenomic procedures to the fungi, the eukaryotic clade with the largest number of sequenced genomes, for which we retrieved a golden standard phylogeny from the literature. Phylogenomic trees based on gene content show, relative to the other methods, a bias in the tree topology that parallels convergence in lifestyle among the species compared, indicating convergence in gene content. CONCLUSIONS: Complete genomes are no guarantee for good or even consistent phylogenies. However, the large amounts of data in genomes enable us to carefully select the data most suitable for phylogenomic inference. In terms of performance, the superalignment approach, combined with restrictive orthology, is the most successful in recovering a fungal phylogeny that agrees with current taxonomic views, and allows us to obtain a high-resolution phylogeny. We provide solid support for what has grown to be a common practice in phylogenomics during its advance in recent years. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Amplification of mitochondrial small subunit ribosomal DNA of polypores and its potential for phylogenetic analysis

There has been a systematic need to seek adequate phylogenetic markers that can be applied in phylogenetic analyses of fungal taxa at various levels. The mitochondrial small subunit ribosomal DNA (mt SSU rDNA) is generally considered to be one of the molecules that are appropriate for phylogenetic analyses at a family level. In order to obtain universal primers for polypores of Hymenomycetes, mt SSU rRNA genes were cloned from Bjerkandera adusta, Ganoderma lucidum, Phlebiopsis gigantea, and Phellinus laevigatus and their sequences were determined. Based on the conserved sequences of cloned genes from polypores and Agrocybe aegerita, PCR primers were designed for amplification and sequencing of mt SSU rDNAs. New primers allowed effective amplification and sequencing of almost full-sized genes from representative species of polypores and related species. Phylogenetic relationships were resolved quite efficiently by mt SSU rDNA sequences, and they proved to be more useful in phylogenetic reconstruction of Ganoderma than nuclear internal transcribed spacer (ITS) rDNA sequences.     

Making the most of mitochondrial genomes--markers for phylogeny, molecular ecology and barcodes in Schistosoma (Platyhelminthes: Digenea)

An increasing number of complete sequences of mitochondrial (mt) genomes provides the opportunity to optimise the choice of molecular markers for phylogenetic and ecological studies. This is particularly the case where mt genomes from closely related taxa have been sequenced; e.g., within Schistosoma. These blood flukes include species that are the causative agents of schistosomiasis, where there has been a need to optimise markers for species and strain recognition. For many phylogenetic and population genetic studies, the choice of nucleotide sequences depends primarily on suitable PCR primers. Complete mt genomes allow individual gene or other mt markers to be assessed relative to one another for potential information content, prior to broad-scale sampling. We assess the phylogenetic utility of individual genes and identify regions that contain the greatest interspecific variation for molecular ecological and diagnostic markers. We show that variable characters are not randomly distributed along the genome and there is a positive correlation between polymorphism and divergence. The mt genomes of African and Asian schistosomes were compared with the available intraspecific dataset of Schistosoma mansoni through sliding window analyses, in order to assess whether the observed polymorphism was at a level predicted from interspecific comparisons. We found a positive correlation except for the two genes (cox1 and nad1) adjoining the putative control region in S. mansoni. The genes nad1, nad4, nad5, cox1 and cox3 resolved phylogenies that were consistent with a benchmark phylogeny and in general, longer genes performed better in phylogenetic reconstruction. Considering the information content of entire mt genome sequences, partial cox1 would not be the ideal marker for either species identification (barcoding) or population studies with Schistosoma species. Instead, we suggest the use of cox3 and nad5 for both phylogenetic and population studies. Five primer pairs designed against Schistosoma mekongi and Schistosoma malayensis were tested successfully against Schistosoma japonicum. In combination, these fragments encompass 20-27% of the variation amongst the genomes (average total length approximately 14,000bp), thus providing an efficient means of encapsulating the greatest amount of variation within the shortest sequence. Comparative mitogenomics provides the basis of a rational approach to molecular marker selection and optimisation.     

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.     

Chloroplast genomic diversity in Bulbophyllum section Macrocaulia (Orchidaceae,Epidendroideae, Malaxideae): Insights into species divergence and adaptive evolution

Bulbophyllum is the largest genus in Orchidaceae with a pantropical distribution. Due to highly significant diversifications, it is considered to be one of the most taxonomically and phylogenetically complex taxa. The diversification pattern and evolutionary adaptation of chloroplast genomes are poorly understood in this species-rich genus, and suitable molecular markers are necessary for species determination and phylogenetic analysis. A natural Asian section Macrocaulia was selected to estimate the interspecific divergence of chloroplast genomes in this study. Here, we sequenced the complete chloroplast genome of four Bulbophyllum species, including three species from section Macrocaulia. The four chloroplast genomes had a typical quadripartite structure with a genome size ranged from 156,182 to 158,524 bp. The chloroplast genomes included 113 unique genes encoding 79 proteins, 30 tRNAs and 4 rRNAs. Comparison of the four chloroplast genomes showed that the three species from section Macrocaulia had similar structure and gene contents, and shared a number of indels, which mainly contribute to its monophyly. In addition, interspecific divergence level was also great. Several exclusive indels and polymorphism SSR loci might be used for taxonomical identification and determining interspecific polymorphisms. A total of 20 intergenic regions and three coding genes of the most variable hotspot regions were proposed as candidate effective molecular markers for future phylogenetic relationships at different taxonomical levels and species divergence in Bulbophyllum. All of chloroplast genes in four Bulbophyllum species were under purifying selection, while 13 sites within six genes exhibited site-specific selection. A whole chloroplast genome phylogenetic analysis based on Maximum Likelihood, Bayesian and Parsimony methods all supported the monophyly of section Macrocaulia and the genus of Bulbophyllum. Our findings provide valuable molecular markers to use in accurately identifying species, clarifying taxonomy, and resolving the phylogeny and evolution of the genus Bulbophyllum. The molecular markers developed in this study will also contribute to further research of conservation of Bulbophyllum species.     

Visualization of the phylogenetic content of five genomes using dekapentagonal maps

The methods presented here summarize phylogenetic relationships of genomes in visually appealing and informative figures. Dekapentagonal maps depict phylogenetic information for orthologous genes present in five genomes, and provide a pre-screen for putatively horizontally transferred genes. If the majority of individual gene phylogenies are unresolved, bipartition histograms provide a means of uncovering and analyzing the plurality consensus. Analyses of genomes representing five photosynthetic bacterial phyla and of the prokaryotic contributions to the eukaryotic cell illustrate the utility of the methods.