Comparison of methods for rooting phylogenetic trees: a case study using Orcuttieae (Poaceae: Chloridoideae).

DNA sequence data (cpDNA trnL intron and nrDNA ITS1 and ITS2) were analyzed to identify relationships within Orcuttieae, a small tribe of endangered grasses endemic to vernal pools in California and Baja California. The tribe includes three genera: Orcuttia, Tuctoria, and Neostapfia. All three genera carry out C(4) photosynthesis but aquatic taxa of Orcuttia lack Kranz anatomy. The unusual habitat preference of the tribe is coupled with the atypical development of C(4) photosynthesis without Kranz anatomy. Furthermore, the tribe has no known close relatives and has been noted to be phylogenetically isolated within the subfamily Chloridoideae. In this study we examine the problem of inferring the root of the tribe in the absence of an identified outgroup, analyze the phylogenetic relationships of the constituent taxa, and evaluate the evolutionary development of C(4) photosynthesis. We compare four methods for inferring the root of the tree: (1) the outgroup method, (2) midpoint rooting, the imposition of a molecular clock for both (3) maximum likelihood (ML) and (4) Bayesian analysis. We examine the consequences of each method for the inferred phylogenetic relationships. Three of the methods (outgroup rooting and the ML and Bayesian molecular clock analyses) suggest that the root of Orcuttieae is between Neostapfia and the Tuctoria/Orcuttia lineage, while midpoint rooting gives a different root. The Bayesian method additionally provides information about probabilities associated with other possible root locations. Assuming that the true root of Orcuttieae is between Neostapfia and the Tuctoria/Orcuttia lineage, our data indicate Neostapfia and Orcuttia are both monophyletic, while Tuctoria is paraphyletic (with no synapomorphies in either dataset) and forming a grade between the other two genera and needs taxonomic revision. Our data support the hypothesis that Orcuttieae was derived from a terrestrial ancestor and evolved specializations to an aquatic environment, including C(4) photosynthesis without Kranz anatomy.

Identifying hybridization events in the presence of coalescence via model selection.

As DNA sequences have become more readily available, it has become increasingly desirable to infer species phylogenies from multigene data sets. Much recent work has centered around the recognition that substantial incongruence in single-gene phylogenies necessitates the development of statistical procedures to estimate species phylogenies that appropriately model the process of evolution at the level of the individual genes. One process that gives rise to variation in the histories of individual genes is incomplete lineage sorting, which is commonly modeled by the coalescent, and thus much current work is focused on proper estimation of species phylogenies under the coalescent model. A second common source of discord in single-gene phylogenies is hybridization, a process that is ubiquitous in many groups of plants and animals. Although methods to incorporate hybridization into phylogenetic estimation have also been developed, only a handful of methods that address both coalescence and hybridization have been proposed. Here, I propose an extension of an existing model that incorporates both of these processes simultaneously by utilizing gene trees for inference in a likelihood framework. The model allows examination of the evidence for hybridization in the presence of incomplete lineage sorting due to deep coalescence via model selection using standard information criteria (e.g., Akaike information criterion and Bayesian information criterion). The potential of the method is evaluated using simulated data.

Detecting hybrid speciation in the presence of incomplete lineage sorting using gene tree incongruence: a model.

The application of phylogenetic inference methods, to data for a set of independent genes sampled randomly throughout the genome, often results in substantial incongruence in the single-gene phylogenetic estimates. Among the processes known to produce discord between single-gene phylogenies, two of the best studied in a phylogenetic context are hybridization and incomplete lineage sorting. Much recent attention has focused on the development of methods for estimating species phylogenies in the presence of incomplete lineage sorting, but phylogenetic models that allow for hybridization have been more limited. Here we propose a model that allows incongruence in single-gene phylogenies to be due to both hybridization and incomplete lineage sorting, with the goal of determining the contribution of hybridization to observed gene tree incongruence in the presence of incomplete lineage sorting. Using our model, we propose methods for estimating the extent of the role of hybridization in both a likelihood and a Bayesian framework. The performance of our methods is examined using both simulated and empirical data.

Coalescent time distributions in trees of arbitrary size.

The relationship between speciation times and the corresponding times of gene divergence is of interest in phylogenetic inference as a means of understanding the past evolutionary dynamics of populations and of estimating the timing of speciation events. It has long been recognized that gene divergence times might substantially pre-date speciation events. Although the distribution of the difference between these has previously been studied for the case of two populations, this distribution has not been explicitly computed for larger species phylogenies. Here we derive a simple method for computing this distribution for trees of arbitrary size. A two-stage procedure is proposed which (i) considers the probability distribution of the time from the speciation event at the root of the species tree to the gene coalescent time conditionally on the number of gene lineages available at the root; and (ii) calculates the probability mass function for the number of gene lineages at the root. This two-stage approach dramatically simplifies numerical analysis, because in the first step the conditional distribution does not depend on an underlying species tree, while in the second step the pattern of gene coalescence prior to the species tree root is irrelevant. In addition, the algorithm provides intuition concerning the properties of the distribution with respect to the various features of the underlying species tree. The methodology is complemented by developing probabilistic formulae and software, written in R. The method and software are tested on five-taxon species trees with varying levels of symmetry. The examples demonstrate that more symmetric species trees tend to have larger mean coalescent times and are more likely to have a unimodal gamma-like distribution with a long right tail, while asymmetric trees tend to have smaller mean coalescent times with an exponential-like distribution. In addition, species trees with longer branches generally have shorter mean coalescent times, with branches closest to the root of the tree being most influential.

Vitamin D receptor gene polymorphisms and susceptibility M. tuberculosis in native Paraguayans.

Tuberculosis (TB) is a significant health problem for most of the world's populations, and prevalence among indigenous groups is typically higher than among their nonindigenous neighbors. Native South Americans experience high rates of TB, but while research in several other world populations indicates that susceptibility is multifactorial, polygenic, and population-specific, little work has been undertaken to investigate factors involved in Native American susceptibility. We conducted a family-based association study to examine immunologically relevant polymorphisms of a candidate gene, the vitamin D receptor, in conjunction with three measures of TB status in two Native Paraguayan populations, the Aché and the Avá. This is the first large-scale genetic analysis of Native South Americans to examine susceptibility to both infection and disease following exposure to M. tuberculosis. These two types of susceptibility reflect differences in innate and acquired immunity that have proven difficult to elucidate in other populations. Our results indicate that among the Aché, the FokI F allele protects individuals from infection, while the TaqI t allele protects against active disease but not infection. In particular, FF homozygotes are 17 times more likely to test positive for exposure to TB, but no more likely to have ever been diagnosed with active TB. TT individuals are 42 times less likely to mount a delayed-type hypersensitivity response, and the T allele was significantly more likely to have been transmitted to offspring who have been diagnosed with active TB. This ongoing research is of vital importance to indigenous groups of the Americas, because if there is a population-specific component to TB susceptibility, it will likely prove most effective to incorporate this into future treatment and prevention strategies.

Inconsistency of phylogenetic estimates from concatenated data under coalescence.

Although multiple gene sequences are becoming increasingly available for molecular phylogenetic inference, the analysis of such data has largely relied on inference methods designed for single genes. One of the common approaches to analyzing data from multiple genes is concatenation of the individual gene data to form a single supergene to which traditional phylogenetic inference procedures - e.g., maximum parsimony (MP) or maximum likelihood (ML) - are applied. Recent empirical studies have demonstrated that concatenation of sequences from multiple genes prior to phylogenetic analysis often results in inference of a single, well-supported phylogeny. Theoretical work, however, has shown that the coalescent can produce substantial variation in single-gene histories. Using simulation, we combine these ideas to examine the performance of the concatenation approach under conditions in which the coalescent produces a high level of discord among individual gene trees and show that it leads to statistically inconsistent estimation in this setting. Furthermore, use of the bootstrap to measure support for the inferred phylogeny can result in moderate to strong support for an incorrect tree under these conditions. These results highlight the importance of incorporating variation in gene histories into multilocus phylogenetics.