Multiple Gene Trees Research Articles

A roadmap of phylogenomic methods for studying polyploid plant genera.

Phylogenetic inference of polyploid species is the first step towards understanding their patterns of diversification. In this paper, we review the challenges and limitations of inferring species relationships of polyploid plants using traditional phylogenetic sequencing approaches, as well as the mischaracterization of the species tree from single or multiple gene trees. We provide a roadmap to infer interspecific relationships among polyploid lineages by comparing and evaluating the application of current phylogenetic, phylogenomic, transcriptomic, and whole-genome approaches using different sequencing platforms. For polyploid species tree reconstruction, we assess the following criteria: (1) the amount of prior information or tools required to capture the genetic region(s) of interest; (2) the probability of recovering homeologs for polyploid species; and (3) the time efficiency of downstream data analysis. Moreover, we discuss bioinformatic pipelines that can reconstruct networks of polyploid species relationships. In summary, although current phylogenomic approaches have improved our understanding of reticulate species relationships in polyploid-rich genera, the difficulties of recovering reliable orthologous genes and sorting all homeologous copies for allopolyploids remain a challenge. In the future, assembled long-read sequencing data will assist the recovery and identification of multiple gene copies, which can be particularly useful for reconstructing the multiple independent origins of polyploids.

Visualizing incompatibilities in phylogenetic trees using consensus outlines.

Phylogenetic analysis frequently leads to the creation of many phylogenetic trees, either from using multiple genes or methods, or through bootstrapping or Bayesian analysis. A consensus tree is often used to summarize what the trees have in common. Consensus networks were introduced to also allow the visualization of the main incompatibilities among the trees. However, in practice, such networks often contain a large number of nodes and edges, and can be non-planar, making them difficult to interpret. Here, we introduce the new concept of a phylogenetic consensus outline, which provides a planar visualization of incompatibilities in the input trees, without the complexities of a consensus network. Furthermore, we present an effective algorithm for its computation. We demonstrate its usage and explore how it compares to other methods on a Bayesian phylogenetic analysis of languages using data from a published database and on multiple gene trees from a published study on water lilies.

Comparative transcriptomic analysis unveils the deep phylogeny and secondary metabolite evolution of 116 Camellia plants.

Camellia plants include more than 200 species of great diversity and immense economic, ornamental, and cultural values. We sequenced the transcriptomes of 116 Camellia plants from almost all sections of the genus Camellia. We constructed a pan-transcriptome of Camellia plants with 89 394 gene families and then resolved the phylogeny of genus Camellia based on 405 high-quality low-copy core genes. Most of the inferred relationships are well supported by multiple nuclear gene trees and morphological traits. We provide strong evidence that Camellia plants shared a recent whole genome duplication event, followed by large expansions of transcription factor families associated with stress resistance and secondary metabolism. Secondary metabolites, particularly those associated with tea quality such as catechins and caffeine, were preferentially heavily accumulated in the Camellia plants from section Thea. We thoroughly examined the expression patterns of hundreds of genes associated with tea quality, and found that some of them exhibited significantly high expression and correlations with secondary metabolite accumulations in Thea species. We also released a web-accessible database for efficient retrieval of Camellia transcriptomes. The reported transcriptome sequences and obtained novel findings will facilitate the efficient conservation and utilization of Camellia germplasm towards a breeding program for cultivated tea, camellia, and oil-tea plants.

The British chalk specialist Lecidea lichenicola auct. revealed as a new genus of Lichinomycetes

The lichen, to which the name Lecidea lichenicola is found to have been misapplied, was first described from England and is an extreme specialist of chalk pebbles. It has long been known that it is not closely related to Lecidea in the strict sense, but its true evolutionary relationships have been unknown. Here we use metagenome-assembled genome data to place this fungus in a six-locus phylogeny of Ascomycota, and find strong support for its placement in the class Lichinomycetes. Multiple gene trees using existing data from Lichinomycetes support its further placement within the family Lichinaceae. Based on a revision of types and original descriptions, we conclude that the earliest name for this species is Lecidea obsoleta (syn. Thrombium cretaceum). We neotypify that name by a modern collection and accommodate it in the new genus Watsoniomyces. Type and other original material of L. lichenicola (syn. Discocera lichenicola) was re-examined and found not to be on chalk and to represent a different lichen, Trapelia glebulosa. Watsoniomyces is the first described member of Lichinomycetes with an endolithic thallus.

UglyTrees: a browser-based multispecies coalescent tree visualizer.

SummaryVisualization is a vital task in phylogenetics and yet there is a deficit in programs which visualize the multispecies coalescent (MSC) model. UglyTrees (UT) is an easy-to-use program for visualizing multiple gene trees embedded within a single species trees. The mapping between gene and species nodes is automatically detected allowing for ready access to the program. UT can scrape the contents of a website for MSC analyses, enabling the sharing of interactive MSC figures through optional parameters in the URL. If a posterior distribution is uploaded, the transitions between MSC states are animated allowing the visual tracking of trees throughout the sequence.Availability and implementationUT runs in all major web browsers including mobile devices, and is hosted at www.uglytrees.nz. The MIT-licensed code is available at https://github.com/UglyTrees/uglytrees.github.io.

The published complete mitochondrial genome of Eptesicus serotinus is a chimera of Vespertilio sinensis and Hypsugo alaschanicus (Mammalia: Chiroptera)

The mitogenome of Eptesicus serotinus (Serotine bat) was published in 2013 with GenBank accession number KF111725 and NCBI Reference Sequence number NC_022474. This sequence was placed with Vespertilio sinensis (Asian parti-colored bat) in a COI gene tree but with Hypsugo alashanicus (Alashanian pipistrelle) in a cytochrome b gene tree. Direct comparison of mitogenomes showed that 92.4% of this mitogenome is similar to Vespertilio sinensis, 5.9% to Hypsugo alaschanicus, and that 1.6% of the mitogenome could not be attributed to either species, or any other species. This mitogenome has been re-used in at least 17 phylogenies. Our findings suggest that mitogenomes are best verified with multiple gene trees, followed by direct comparison of sequences. We conclude that greater vigilance is warranted to ensure that problematic sequences do not enter the scientific record, and are not re-used in subsequent studies.

Modern Phylogenomics: Building Phylogenetic Trees Using the Multispecies Coalescent Model.

The multispecies coalescent (MSC) model provides a compelling framework for building phylogenetic trees from multilocus DNA sequence data. The pure MSC is best thought of as a special case of so-called "multispecies network coalescent" models, in which gene flow is allowed among branches of the tree, whereas MSC methods assume there is no gene flow between diverging species. Early implementations of the MSC, such as "parsimony" or "democratic vote" approaches to combining information from multiple gene trees, as well as concatenation, in which DNA sequences from multiple gene trees are combined into a single "supergene," were quickly shown to be inconsistent in some regions of tree space, in so far as they converged on the incorrect species tree as more gene trees and sequence data were accumulated. The anomaly zone, a region of tree space in which the most frequent gene tree is different from the species tree, is one such region where many so-called "coalescent" methods are inconsistent. Second-generation implementations of the MSC employed Bayesian or likelihood models; these are consistent in all regions of gene tree space, but Bayesian methods in particular are incapable of handling the large phylogenomic data sets currently available. Two-step methods, such as MP-EST and ASTRAL, in which gene trees are first estimated and then combined to estimate an overarching species tree, are currently popular in part because they can handle large phylogenomic data sets. These methods are consistent in the anomaly zone but can sometimes provide inappropriate measures of tree support or apportion error and signal in the data inappropriately. MP-EST in particular employs a likelihood model which can be conveniently manipulated to perform statistical tests of competing species trees, incorporating the likelihood of the collected gene trees on each species tree in a likelihood ratio test. Such tests provide a useful alternative to the multilocus bootstrap, which only indirectly tests the appropriateness of competing species trees. We illustrate these tests and implementations of the MSC with examples and suggest that MSC methods are a useful class of models effectively using information from multiple loci to build phylogenetic trees.

Gene tree parsimony for incomplete gene trees: addressing true biological loss

Motivation Species tree estimation from gene trees can be complicated by gene duplication and loss, and “gene tree parsimony” (GTP) is one approach for estimating species trees from multiple gene trees. In its standard formulation, the objective is to find a species tree that minimizes the total number of gene duplications and losses with respect to the input set of gene trees. Although much is known about GTP, little is known about how to treat inputs containing some incomplete gene trees (i.e., gene trees lacking one or more of the species).ResultsWe present new theory for GTP considering whether the incompleteness is due to gene birth and death (i.e., true biological loss) or taxon sampling, and present dynamic programming algorithms that can be used for an exact but exponential time solution for small numbers of taxa, or as a heuristic for larger numbers of taxa. We also prove that the “standard” calculations for duplications and losses exactly solve GTP when incompleteness results from taxon sampling, although they can be incorrect when incompleteness results from true biological loss. The software for the DP algorithm is freely available as open source code at https://github.com/smirarab/DynaDup.

The impact of phenotypic and molecular data on the inference of Colletotrichum diversity associated with Musa

ABSTRACTDeveloping a comprehensive and reliable taxonomy for the Colletotrichum gloeosporioides species complex will require adopting data standards on the basis of an understanding of how methodological choices impact morphological evaluations and phylogenetic inference. We explored the impact of methodological choices in a morphological and molecular evaluation of Colletotrichum species associated with banana in Brazil. The choice of alignment filtering algorithm has a significant impact on topological inference and the retention of phylogenetically informative sites. Similarly, the choice of phylogenetic marker affects the delimitation of species boundaries, particularly if low phylogenetic signal is confounded with strong discordance, and inference of the species tree from multiple-gene trees. According to both phylogenetic informativeness profiling and Bayesian concordance analyses, the most informative loci are DNA lyase (APN2), intergenic spacer (IGS) between DNA lyase and the mating-type locus MAT1-2-1 (APN2/MAT-IGS), calmodulin (CAL), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), glutamine synthetase (GS), β-tubulin (TUB2), and a new marker, the intergenic spacer between GAPDH and an hypothetical protein (GAP2-IGS). Cornmeal agar minimizes the variance in conidial dimensions compared with potato dextrose agar and synthetic nutrient-poor agar, such that species are more readily distinguishable based on phenotypic differences. We apply these insights to investigate the diversity of Colletotrichum species associated with banana anthracnose in Brazil and report C. musae, C. tropicale, C. theobromicola, and C. siamense in association with banana anthracnose. One lineage did not cluster with any previously described species and is described here as C. chrysophilum.

Fast Construction of Near Parsimonious Hybridization Networks for Multiple Phylogenetic Trees.

Hybridization networks represent plausible evolutionary histories of species that are affected by reticulate evolutionary processes. An established computational problem on hybridization networks is constructing the most parsimonious hybridization network such that each of the given phylogenetic trees (called gene trees) is "displayed" in the network. There have been several previous approaches, including an exact method and several heuristics, for this NP-hard problem. However, the exact method is only applicable to a limited range of data, and heuristic methods can be less accurate and also slow sometimes. In this paper, we develop a new algorithm for constructing near parsimonious networks for multiple binary gene trees. This method is more efficient for large numbers of gene trees than previous heuristics. This new method also produces more parsimonious results on many simulated datasets as well as a real biological dataset than a previous method. We also show that our method produces topologically more accurate networks for many datasets.

Speciation of two gobioid species, Pterogobius elapoides and Pterogobius zonoleucus revealed by multi-locus nuclear and mitochondrial DNA analyses

To understand how geographical differentiation of gobioid fish species led to speciation, two populations of the Pacific Ocean and the Sea of Japan for each of the two gobioid species, Pterogobius elapoides and Pterogobius zonoleucus, were studied in both morphological and molecular features. Analyzing mitochondrial genes, Akihito et al. (2008) suggested that P. zonoleucus does not form a monophyletic clade relative to P. elapoides, indicating that “Sea of Japan P. zonoleucus” and P. elapoides form a clade excluding “Pacific P. zonoleucus” as an outgroup. Because morphological classification clearly distinguish these two species and a gene tree may differ from a population tree, we examined three nuclear genes, S7RP, RAG1, and TBR1, in this work, in order to determine whether nuclear and mitochondrial trees are concordant, thus shedding light on the evolutionary history of this group of fishes. Importantly, nuclear trees were based on exactly the same individuals that were used for the previously published mtDNA trees. The tree based on RAG1 exon sequences suggested a closer relationship of P. elapoides with “Sea of Japan P. zonoleucus”, which was in agreement with the mitochondrial tree. In contrast, S7RP and TBR1 introns recovered a monophyletic P. zonoleucus. If the mitochondrial tree represents the population tree in which P. elapoides evolved from “Sea of Japan P. zonoleucus”, the population size of P. elapoides is expected to be smaller than that of “Sea of Japan P. zonoleucus”. This is because a smaller population of the new species is usually differentiated from a larger population of the ancestral species when the speciation occurred. However, we found no evidence of such a small population size during the evolution of P. elapoides. Therefore, we conclude that the monophyletic P. zonoleucus as suggested by S7RP and TBR1 most likely represents the population tree, which is consistent with the morphological classification. In this case, it is possible that the incongruent mitochondrial and RAG1 trees are either due to incomplete lineage sorting of ancestral polymorphisms or to introgression by hybridization. Because of a smaller effective population size of mitochondria compared with nuclear genes, the introgression might be a more likely scenario in explaining the incongruent mitochondrial tree than the incomplete lineage sorting. Because of smaller effective population size of “Sea of Japan P. zonoleucus” than that of P. elapoides, the direction of the introgression was likely to be from the latter to the former. This evolutionary work of the two gobioid species highlights the need of analyzing multiple gene trees for both nuclear and mitochondrial genes as well as scrutinization of morphological characteristics to obtain a population tree representing the organismal evolutionary history.

An algorithm for constructing parsimonious hybridization networks with multiple phylogenetic trees.

A phylogenetic network is a model for reticulate evolution. A hybridization network is one type of phylogenetic network for a set of discordant gene trees and "displays" each gene tree. A central computational problem on hybridization networks is: given a set of gene trees, reconstruct the minimum (i.e., most parsimonious) hybridization network that displays each given gene tree. This problem is known to be NP-hard, and existing approaches for this problem are either heuristics or making simplifying assumptions (e.g., work with only two input trees or assume some topological properties). In this article, we develop an exact algorithm (called PIRNC) for inferring the minimum hybridization networks from multiple gene trees. The PIRNC algorithm does not rely on structural assumptions (e.g., the so-called galled networks). To the best of our knowledge, PIRNC is the first exact algorithm implemented for this formulation. When the number of reticulation events is relatively small (say, four or fewer), PIRNC runs reasonably efficient even for moderately large datasets. For building more complex networks, we also develop a heuristic version of PIRNC called PIRNCH. Simulation shows that PIRNCH usually produces networks with fewer reticulation events than those by an existing method. PIRNC and PIRNCH have been implemented as part of the software package called PIRN and is available online.

EHive: An Artificial Intelligence workflow system for genomic analysis

BackgroundThe Ensembl project produces updates to its comparative genomics resources with each of its several releases per year. During each release cycle approximately two weeks are allocated to generate all the genomic alignments and the protein homology predictions. The number of calculations required for this task grows approximately quadratically with the number of species. We currently support 50 species in Ensembl and we expect the number to continue to grow in the future.ResultsWe present eHive, a new fault tolerant distributed processing system initially designed to support comparative genomic analysis, based on blackboard systems, network distributed autonomous agents, dataflow graphs and block-branch diagrams. In the eHive system a MySQL database serves as the central blackboard and the autonomous agent, a Perl script, queries the system and runs jobs as required. The system allows us to define dataflow and branching rules to suit all our production pipelines. We describe the implementation of three pipelines: (1) pairwise whole genome alignments, (2) multiple whole genome alignments and (3) gene trees with protein homology inference. Finally, we show the efficiency of the system in real case scenarios.ConclusionseHive allows us to produce computationally demanding results in a reliable and efficient way with minimal supervision and high throughput. Further documentation is available at: http://www.ensembl.org/info/docs/eHive/.

Incomplete Lineage Sorting: Consistent Phylogeny Estimation from Multiple Loci

We introduce a simple computationally efficient algorithm for reconstructing phylogenies from multiple gene trees in the presence of incomplete lineage sorting, that is, when the topology of the gene trees may differ from that of the species tree. We show that our technique is statistically consistent under standard stochastic assumptions, that is, it returns the correct tree given sufficiently many unlinked loci. We also show that it can tolerate moderate estimation errors.

Bayesian Inference of Species Trees from Multilocus Data

Until recently, it has been common practice for a phylogenetic analysis to use a single gene sequence from a single individual organism as a proxy for an entire species. With technological advances, it is now becoming more common to collect data sets containing multiple gene loci and multiple individuals per species. These data sets often reveal the need to directly model intraspecies polymorphism and incomplete lineage sorting in phylogenetic estimation procedures.For a single species, coalescent theory is widely used in contemporary population genetics to model intraspecific gene trees. Here, we present a Bayesian Markov chain Monte Carlo method for the multispecies coalescent. Our method coestimates multiple gene trees embedded in a shared species tree along with the effective population size of both extant and ancestral species. The inference is made possible by multilocus data from multiple individuals per species.Using a multiindividual data set and a series of simulations of rapid species radiations, we demonstrate the efficacy of our new method. These simulations give some insight into the behavior of the method as a function of sampled individuals, sampled loci, and sequence length. Finally, we compare our new method to both an existing method (BEST 2.2) with similar goals and the supermatrix (concatenation) method. We demonstrate that both BEST and our method have much better estimation accuracy for species tree topology than concatenation, and our method outperforms BEST in divergence time and population size estimation.

Reticulation, Data Combination, and Inferring Evolutionary History: An Example from Danthonioideae (Poaceae)

We explore the potential impact of conflicting gene trees on inferences of evolutionary history above the species level. When conflict between gene trees is discovered, it is common practice either to analyze the data separately or to combine the data having excluded the conflicting taxa or data partitions for those taxa (which are then recoded as missing). We demonstrate an alternative approach, which involves duplicating conflicting taxa in the matrix, such that each duplicate is represented by one partition only. This allows the combination of all available data in standard phylogenetic analyses, despite reticulations. We show how interpretation of contradictory gene trees can lead to conflicting inferences of both morphological evolution and biogeographic history, using the example of the pampas grasses, Cortaderia. The characteristic morphological syndrome of Cortaderia can be inferred as having arisen multiple times (chloroplast DNA [cpDNA]) or just once (nuclear ribosomal DNA [nrDNA]). The distributions of species of Cortaderia and related genera in Australia/New Guinea, New Zealand, and South America can be explained by few (nrDNA) or several (cpDNA) dispersals between the southern continents. These contradictions can be explained by past hybridization events, which have linked gains of complex morphologies with unrelated chloroplast lineages and have erased evidence of dispersals from the nuclear genome. Given the discrepancies between inferences based on the gene trees individually, we urge the use of approaches such as ours that take multiple gene trees into account.

The utility of nuclear conserved ortholog set II (COSII) genomic regions for species-level phylogenetic inference in Lycium (Solanaceae)

The identification of genomic regions with sufficient variation to elucidate fine-scale relationships among closely related species is a major goal of phylogenetic systematics. However, the accumulation of such multi-locus data sets brings its own challenges, given that gene trees do not necessarily represent the true species tree. Using genomic tools developed for Solanum (Solanaceae), we have evaluated the utility of nuclear conserved ortholog set II (COSII) regions for phylogenetic inference in tribe Lycieae (Solanaceae). Five COSII regions, with intronic contents ranging from 68% to 91%, were sequenced in 10 species. Their phylogenetic utility was assessed and compared with data from more commonly used nuclear (GBSSI, nrITS) and cpDNA spacer data. We compared the effectiveness of a traditional total evidence concatenation approach versus the recently developed Bayesian estimation of species trees (BEST) method to infer species trees given multiple independent gene trees. All of the sampled COSII regions had high numbers of parsimony-informative (PI) characters, and two of the COSII regions had more PI characters than the GBSSI, ITS, and cpDNA spacer data sets combined. COSII markers are a promising new tool for phylogenetic inference in Solanaceae, and should be explored in related groups. Both the concatenation and BEST approaches yielded similar topologies; however, when multiple individuals with polyphyletic alleles were included, BEST was clearly the more robust approach for inferring species trees in the presence of gene tree incongruence.

Maximum tree: a consistent estimator of the species tree

We propose a model based approach to use multiple gene trees to estimate the species tree. The coalescent process requires that gene divergences occur earlier than species divergences when there is any polymorphism in the ancestral species. Under this scenario, speciation times are restricted to be smaller than the corresponding gene split times. The maximum tree (MT) is the tree with the largest possible speciation times in the space of species trees restricted by available gene trees. If all populations have the same population size, the MT is the maximum likelihood estimate of the species tree. It can be shown the MT is a consistent estimator of the species tree even when the MT is built upon the estimates of the true gene trees if the gene tree estimates are statistically consistent. The MT converges in probability to the true species tree at an exponential rate.

Does organismal pedigree impact the magnitude of topological congruence among gene trees for unlinked loci?

One of the fundamental assumptions in the multi-locus approach to phylogeographic studies is that unlinked loci have independent genealogies. For this reason, congruence among gene trees from unlinked loci is normally interpreted as support for the existence of external forces that may have concordantly shaped the topology of multiple gene trees. However, it is also important to address and quantify the possibility that gene trees within a given species are all inherently constrained to some degree by their shared organismal pedigree, and thus in this strict sense are not entirely independent. Here we demonstrate by computer simulations that gene trees from a shared pedigree tend to display higher topological concordance than do gene trees from independent pedigrees with the same demographic parameters, but we also show that these constraining effects are normally minor in comparison to the much higher degree of topological concordance that can routinely emerge from external phylogeographic shaping forces. The topology-constraining effect of a shared pedigree decreases as effective population size increases, and becomes almost negligible in a random mating population of more than 1,000 individuals. Moreover, statistical detection of the pedigree effect requires a relatively large number of unlinked loci that far exceed what is typically used in current phylogeographic studies. Thus, with the possible exception of extremely small populations, multiple unlinked genes within a pedigree can indeed be assumed, for most practical purposes, to have independent genealogical histories.

Estimating Species Phylogeny from Gene-Tree Probabilities Despite Incomplete Lineage Sorting: An Example from Melanoplus Grasshoppers

Estimating phylogenetic relationships among closely related species can be extremely difficult when there is incongruence among gene trees and between the gene trees and the species tree. Here we show that incorporating a model of the stochastic loss of gene lineages by genetic drift into the phylogenetic estimation procedure can provide a robust estimate of species relationships, despite widespread incomplete sorting of ancestral polymorphism. This approach is applied to a group of montane Melanoplus grasshoppers for which genealogical discordance among loci and incomplete lineage sorting obscures any obvious phylogenetic relationships among species. Unlike traditional treatments where gene trees estimated using standard phylogenetic methods are implicitly equated with the species tree, with the coalescent-based approach the species tree is modeled probabilistically from the estimated gene trees. The estimated species phylogeny (the ESP) is calculated for the grasshoppers from multiple gene trees reconstructed for nuclear loci and a mitochondrial gene. This empirical application is coupled with a simulation study to explore the performance of the coalescent-based approach. Specifically, we test the accuracy of the ESP given the data based on analyses of simulated data matching the multilocus data collected in Melanoplus (i.e., data were simulated for each locus with the same number of base pairs and locus-specific mutational models). The results of the study show that ESPs can be computed using the coalescent-based approach long before reciprocal monophyly has been achieved, and that these statistical estimates are accurate. This contrasts with analyses of the empirical data collected in Melanoplus and simulated data based on concatenation of multiple loci, for which the incomplete lineage sorting of recently diverged species posed significant problems. The strengths and potential challenges associated with incorporating an explicit model of gene-lineage coalescence into the phylogenetic procedure to obtain an ESP, as illustrated by application to Melanoplus, versus concatenation and consensus approaches are discussed. This study represents a fundamental shift in how species relationships are estimated - the relationship between the gene trees and the species phylogeny is modeled probabilistically rather than equating gene trees with a species tree.