Phylogenetic trees based on gene content

Comparing gene content between species can be a useful approach for reconstructing phylogenetic trees. In this paper, we derive a maximum-likelihood estimation of evolutionary distance between species under a simple model of gene genesis and gene loss. Using simulated data on a biological tree with 107 taxa (and on a number of randomly generated trees), we compare the accuracy of tree reconstruction using this ML distance measure to an earlier ad hoc distance. We then compare these distance-based approaches to a character-based tree reconstruction method (Dollo parsimony) which seems well suited to the analysis of gene content data. To simplify simulations, we give a formal proof of the well-known 'fact' that the Dollo parsimony score is independent of the choice of root. Our results show a consistent trend, with the character-based method and ML distance measure outperforming the earlier ad hoc distance method. AVAILABILITY: http://www.ab.informatik.uni-tuebingen.de/software/genecontent/welcome_en.html     

Testing hybridization hypotheses and evaluating the evolutionary potential of hybrids in mangrove plant species

Natural hybridization is of marked importance from global to local biological diversity. In mangroves, species ranges overlap extensively with one another and species share a long overlap of flowering time. Although hybridization has been suggested, patterns of hybridization and the evolutionary potential of hybrids are not yet fully understood. This study provides molecular evidence for the parental origins and status of hybrids in the dominant mangrove genus Rhizophora based on comparisons of chloroplast and nuclear phylogenies and estimations of genetic relatedness and structure from inter‐simple sequence repeat (ISSR) markers. Phylogenetic analyses indicate that almost all species can act as maternal parents to hybrids and that hybridization can be bidirectional. Bayesian analyses indicate that hybrids are simple F₁s, and no trace of backcrossing was detected within populations. Hybridization, for the most part, occurs almost only locally and dispersal of hybrid individuals is limited beyond the hybrid sites.     

Testing fundamental evolutionary hypotheses

Sober and Steel (J. Theor. Biol. 218, 395-408) give important limits on the use of current models with sequence data for studying ancient aspects of evolution; but they go too far in suggesting that several fundamental aspects of evolutionary theory cannot be tested in a normal scientific manner. To the contrary, we show examples of how some alternatives to the theory of descent can be formulated in such a way that they lead to predictions that can be evaluated (and rejected). The critical factor is a logical formulation of the alternatives, even though not all possible alternatives can be tested simultaneously. Similarly, some of the limits using DNA sequence data can be overcome by other types of sequence derived characters. The uniqueness (or not) of the origin of life, though still difficult, is similarly amenable to the testing of alternative hypotheses.     

Testing hypotheses in order

In certain circumstances, one wishes to test one hypothesisonly if certain other hypotheses have been rejected. This orderingof hypotheses simplifies the task of controlling the probabilityof rejecting any true hypothesis. In an example from an observationalstudy, a treated group is shown to be further from both of twocontrol groups than the two control groups are from each other.     

Molecular trees of trypanosomes incongruent with fossil records of hosts

Molecular trees of trypanosomes have confirmed conventionally accepted genera, but often produce topologies that are incongruent with knowledge of the evolution, systematics, and biogeography of hosts and vectors. These distorted topologies result largely from incorrect assumptions about molecular clocks. A host-based phylogenetic tree could serve as a broad outline against which the reasonability of molecular phylogenies could be evaluated. The host-based tree of trypanosomes presented here supports the " invertebrate first " hypothesis of trypanosome evolution, supports the monophyly of Trypanosomatidae, and indicates the digenetic lifestyle arose three times. An area cladogram of Leishmania supports origination in the Palaearctic during the Palaeocene.     

Testing some hypotheses on human evolution and sexual dimorphism

Research on human evolution and sexual dimorphism motivates an interesting test problem. In studying hominid phylogeny it is of interest to test whether parallel evolution plays a role. With regard to sexual dimorphism it is of interest to known whether the directions of sexual dimorphism in the populations being compared are the same. We show that testing these two problems gives rise to the same type of hypothesis testing, viz. the problem of testing the hypothesis that the means of independent, normally distributed random vectors with unit covariance matrices are situated on a straight line through the origin. A test is proposed and applied to study the sexual dimorphism of 20 recent skull populations. In this example the hypothesis of equal directions of sexual dimorphism is rejected. The classical theory of constructing multiple discriminant functions (canonical variates) is adapted to the problem of comparing sexual dimorphisms.     

Testing some hypotheses on Human Evolution and sexual dimorphism

Research on human evolution and sexual dimorphism motivates an interesting test problem. In studying hominid phylogeny it is of interest to test whether parallel evolution plays a role. With regard to sexual dimorphism it is of interest to know whether the directions of sexual dimorphism in the populations being compared are the same. We show that testing these two problems gives rise to the same type of hypothesis testing, viz. the problem of testing the hypothesis that the means of independent, normally distributed random vectors with unit covariance matrices are situated on a straight line through the origin. A test is proposed and applied to study the sexual dimorphism of 20 recent skull populations. In this example the hypothesis of equal directions of sexual dimorphism is rejected. The classical theory of constructing multiple discriminant functions (canonical variates) is adapted to the problem of comparing sexual dimorphisms.     

Testing hypotheses regarding the genetics of adaptation

Many of the hypotheses regarding the genetics of adaptation require that one know specific details about the genetic basis of complex traits, such as the number and effects of the loci involved. Developments in molecular biology have made it possible to create relatively dense maps of markers that can potentially be used to map genes underlying specific traits. However, there are a number of reasons to doubt that such mapping will provide the level of resolution necessary to specifically address many evolutionary questions. Moreover, evolutionary change is built upon the substitution of individual mutations, many of which may now be cosegregating in the same allele. In order for this developing area not to become a mirage that traps the efforts of an entire field, the genetic dissection of adaptive traits should be conducted within a strict hypothesis-testing framework and within systems that promise a reasonable chance of identifying the specific genetic changes of interest. Continuing advances in molecular technology may lead the way here, but some form of genetic testing is likely to be forever required.     

JML: testing hybridization from species trees

I introduce the software JML that tests for the presence of hybridization in multispecies sequence data sets by posterior predictive checking following Joly, McLenachan and Lockhart (2009, American Naturalist 174, e54). Although their method could potentially be applied on any data set, the lack of appropriate software made its application difficult. The software JML thus fills a need for an easy application of the method but also includes improvements such as the possibility to incorporate uncertainty in the species tree topology. The JML software uses a posterior distribution of species trees, population sizes and branch lengths to simulate replicate sequence data sets using the coalescent with no migration. A test quantity, defined as the minimum pairwise sequence distance between sequences of two species, is then evaluated on the simulated data sets and compared to the one estimated from the original data. Because the test quantity is a good predictor of hybridization events, departure from the bifurcating species tree model could be interpreted as evidence of hybridization. Software performance in terms of computing time is evaluated for several parameters. I also show an application example of the software for detecting hybridization among native diploid North American roses.     

Testing hypotheses about interclass correlations from familial data

Testing problems concerning interclass correlations from familial data are considered in the case where the number of siblings varies among families. Under the assumption of multivariate normality, two test procedures are proposed for testing the hypothesis that an interclass correlation is equal to a specified value. To compare the properties of the tests, including a likelihood ratio test, Monte Carlo experiments are performed. Several test statistics are derived for testing whether two variables about a parent and child are uncorrelated. The proposed tests are compared with previous test procedures, using Monte Carlo simulation. A general procedure for finding confidence intervals for interclass correlations is also derived.     

Estimating species trees from unrooted gene trees

In this study, we develop a distance method for inferring unrooted species trees from a collection of unrooted gene trees. The species tree is estimated by the neighbor joining (NJ) tree built from a distance matrix in which the distance between two species is defined as the average number of internodes between two species across gene trees, that is, average gene-tree internode distance. The distance method is named NJ(st) to distinguish it from the original NJ method. Under the coalescent model, we show that if gene trees are known or estimated correctly, the NJ(st) method is statistically consistent in estimating unrooted species trees. The simulation results suggest that NJ(st) and STAR (another coalescence-based method for inferring species trees) perform almost equally well in estimating topologies of species trees, whereas the Bayesian coalescence-based method, BEST, outperforms both NJ(st) and STAR. Unlike BEST and STAR, the NJ(st) method can take unrooted gene trees to infer species trees without using an outgroup. In addition, the NJ(st) method can handle missing data and is thus useful in phylogenomic studies in which data sets often contain missing loci for some individuals.     

Using patterns of genetic structure based on microsatellite loci to test hypotheses of current hybridization,ancient hybridization and incomplete lineage sorting in Conradina (Lamiaceae)

Hybridization and/or incomplete sorting of ancestral polymorphism are commonly implicated to explain discordant phylogenetic analyses of closely related species complexes. One genus in which these phenomena have been suggested to have played major roles based on phylogenetic data is Conradina, a genus of mints (Lamiaceae) endemic to the southeastern USA containing several endangered species. The goals of this study were to use microsatellite data to better understand patterns of genetic structure in Conradina, to test hypotheses of recent or ancient hybridization and incomplete lineage sorting, and to clarify species boundaries. Individuals from 55 populations representing all Conradina species were genotyped using 10 microsatellite loci. Analyses of the patterns of genetic structure in Conradina revealed a clear differentiation of populations following recognized species boundaries, indicating that species have diverged from one another genetically and interspecific hybridization has not occurred recently. Neither ancient hybridization nor incomplete lineage sorting is supported as the sole cause of species nonmonophyly, suggesting that both may have contributed to patterns found in phylogenetic trees; however, analyses of other types of data may be more appropriate to distinguish between these two hypotheses. Because all described species appear to be valid entities, the current listing status of most endangered species of Conradina is appropriate; however, populations of Conradina canescens are genetically differentiated into three groups, each of which may merit species status, and several recently discovered populations of Conradina in Dunn's Creek State Park in Florida are highly differentiated genetically and also appear to represent a new species.     

Relationships between gene trees and species trees

It is well known that a phylogenetic tree (gene tree) constructed from DNA sequences for a genetic locus does not necessarily agree with the tree that represents the actual evolutionary pathway of the species involved (species tree). One of the important factors that cause this difference is genetic polymorphism in the ancestral species. Under the assumption of neutral mutations, this problem can be studied by evaluating the probability (P) that a gene tree has the same topology as that of the species tree. When one gene (allele) is used from each of the species involved, the probability can be expressed as a simple function of Ti = ti/(2N), where ti is the evolutionary time measured in generations for the ith internodal branch of the species tree and N is the effective population size. When any of the Ti's is less than 1, the probability P becomes considerably less than 1.0. This probability cannot be substantially increased by increasing the number of alleles sampled from a locus. To increase the probability, one has to use DNA sequences from many different loci that have evolved independently of each other.      

Testing hypotheses of adaptive variation in cricket ovipositor lengths

We experimentally tested a series of hypotheses proposed by Masaki (1979, 1986) for the evolution of ovipositor length in crickets. Female crickets use the ovipositor to bury eggs in the soil, where it was hypothesized to protect their eggs from desiccation, cold and other disturbance. However, we found no effect of depth on the overwinter survival of eggs of three species of Nemobiinae. The probability of hatchlings reaching the soil surface was negatively correlated with depth documenting a significant cost to females laying eggs deep in the soil. Hatchling survival may be an important agent of selection on ovipositor length in habitats of different soil moistures. Hatchling survival in the soil was also correlated with body size, which may impose a constraint on egg-size fecundity trade-offs. Females of a bivoltine population of Allonemobius socius lay eggs at shallower depths when reared under summer compared to fall conditions and, therefore, may be able to respond to selection through behavioral plasticity when morphological adaptation is constrained by allometry.     

Testing evolutionary hypotheses about the phylotypic period of zebrafish

Vertebrate embryos pass through a period of morphological similarity, the phylotypic period. Since Haeckel's biogenetic law of recapitulation, proximate and ultimate evolutionary causes of such similarity of embryos were discussed. We test predictions about changes in phenotypic and genetic variances that were derived from three hypotheses about the evolutionary origin of the phylotypic stage, i.e. random, epigenetic effects, and stabilizing selection. The random hypothesis predicts increasing values for phenotypic variances and stable or increasing values for genetic variances; the epigenetic effects hypothesis predicts declining values for phenotypic variances but stable or increasing values of genetic variances, and the stabilizing selection predicts stable phenotypic variances but decreasing genetic variances. We studied zebrafish as a model species, because it can be bred in large numbers as necessary for a quantitative genetics breeding design. A half-sib breeding scheme provided estimates of additive genetic variances from 11 embryonic characters from 12 through to 24 hr after fertilization, i.e. before, during (15-19 hr), and after the phylotypic period. Because additive genetic variances are size dependent, we calculated narrow-sense heritabilities as a size independent gauge of genetic contributions to the phenotype. The results show declining phenotypic variances and stable heritabilities. In conclusion, we reject the random and the stabilizing selection hypotheses and favor ideas about epigenetic effects that constrain the early embryonic development. Additive genetic variance during the phylotypic stage makes it accessible for evolution, thus explaining in a simple and straightforward way why the phylotypic period differs among vertebrates in timing, duration, and morphologies.     

Testing allele homogeneity: the problem of nested hypotheses

Background

Coronary artery disease (CAD), and one of its intermediate risk factors, dyslipidemia, possess a demonstrable genetic component, although the genetic architecture is incompletely defined. We previously reported a linkage peak on chromosome 5q31-33 for early-onset CAD where the strength of evidence for linkage was increased in families with higher mean low density lipoprotein-cholesterol (LDL-C). Therefore, we sought to fine-map the peak using association mapping of LDL-C as an intermediate disease-related trait to further define the etiology of this linkage peak. The study populations consisted of 1908 individuals from the CATHGEN biorepository of patients undergoing cardiac catheterization; 254 families (N = 827 individuals) from the GENECARD familial study of early-onset CAD; and 162 aorta samples harvested from deceased donors. Linkage disequilibrium-tagged SNPs were selected with an average of one SNP per 20 kb for 126.6-160.2 MB (region of highest linkage) and less dense spacing (one SNP per 50 kb) for the flanking regions (117.7-126.6 and 160.2-167.5 MB) and genotyped on all samples using a custom Illumina array. Association analysis of each SNP with LDL-C was performed using multivariable linear regression (CATHGEN) and the quantitative trait transmission disequilibrium test (QTDT; GENECARD). SNPs associated with the intermediate quantitative trait, LDL-C, were then assessed for association with CAD (i.e., a qualitative phenotype) using linkage and association in the presence of linkage (APL; GENECARD) and logistic regression (CATHGEN and aortas).

Results

We identified four genes with SNPs that showed the strongest and most consistent associations with LDL-C and CAD: EBF1, PPP2R2B, SPOCK1, and PRELID2. The most significant results for association of SNPs with LDL-C were: EBF1, rs6865969, p = 0.01; PPP2R2B, rs2125443, p = 0.005; SPOCK1, rs17600115, p = 0.003; and PRELID2, rs10074645, p = 0.0002). The most significant results for CAD were EBF1, rs6865969, p = 0.007; PPP2R2B, rs7736604, p = 0.0003; SPOCK1, rs17170899, p = 0.004; and PRELID2, rs7713855, p = 0.003.

Conclusion

Using an intermediate disease-related quantitative trait of LDL-C we have identified four novel CAD genes, EBF1, PRELID2, SPOCK1, and PPP2R2B. These four genes should be further examined in future functional studies as candidate susceptibility loci for cardiovascular disease mediated through LDL-cholesterol pathways.     

Testing evolutionary hypotheses for phenotypic divergence using landscape genetics

Understanding the evolutionary causes of phenotypic variation among populations has long been a central theme in evolutionary biology. Several factors can influence phenotypic divergence, including geographic isolation, genetic drift, divergent natural or sexual selection, and phenotypic plasticity. But the relative importance of these factors in generating phenotypic divergence in nature is still a tantalizing and unresolved problem in evolutionary biology. The origin and maintenance of phenotypic divergence is also at the root of many ongoing debates in evolutionary biology, such as the extent to which gene flow constrains adaptive divergence ( Garant et al. 2007 ) and the relative importance of genetic drift, natural selection, and sexual selection in initiating reproductive isolation and speciation ( Coyne & Orr 2004 ). In this issue, Wang & Summers (2010) test the causes of one of the most fantastic examples of phenotypic divergence in nature: colour pattern divergence among populations of the strawberry poison frog (Dendrobates pumilio) in Panama and Costa Rica ( Fig. 1 ). This study provides a beautiful example of the use of the emerging field of landscape genetics to differentiate among hypotheses for phenotypic divergence. Using landscape genetic analyses, Wang & Summers were able to reject the hypotheses that colour pattern divergence is due to isolation‐by‐distance (IBD) or landscape resistance. Instead, the hypothesis left standing is that colour divergence is due to divergent selection, in turn driving reproductive isolation among populations with different colour morphs. More generally, this study provides a wonderful example of how the emerging field of landscape genetics, which has primarily been applied to questions in conservation and ecology, now plays an essential role in evolutionary research.

Figure 1 Open in figure viewer PowerPoint Divergent colour morphs observed among populations of the strawberry poison frog, Dendrobates pumilio. Frogs are from San Cristobal (upper left), Cerro Brujo (upper right), Bastimentos (lower right), and Agua (lower left).