Rates of molecular evolution and the fraction of nucleotide positions free to vary

Evolution: Morphological saturation and release in mammals

Mammal evolution in the Mesozoic was thought to be heavily constrained by competition and predation by dinosaurs. A new study suggests that placental mammals remained constrained for several million years after non-avian dinosaurs perished, perhaps due to competition from archaic mammals.

Model selection may not be a mandatory step for phylogeny reconstruction

Determining the most suitable model for phylogeny reconstruction constitutes a fundamental step in numerous evolutionary studies. Over the years, various criteria for model selection have been proposed, leading to debate over which criterion is preferable. However, the necessity of this procedure has not been questioned to date. Here, we demonstrate that although incongruency regarding the selected model is frequent over empirical and simulated data, all criteria lead to very similar inferences. When topologies and ancestral sequence reconstruction are the desired output, choosing one criterion over another is not crucial. Moreover, skipping model selection and using instead the most parameter-rich model, GTR+I+G, leads to similar inferences, thus rendering this time-consuming step nonessential, at least under current strategies of model selection.

Perspectives from the Avian Phylogenomics Project: Questions that Can Be Answered with Sequencing All Genomes of a Vertebrate Class

The rapid pace of advances in genome technology, with concomitant reductions in cost, makes it feasible that one day in our lifetime we will have available extant genomes of entire classes of species, including vertebrates. I recently helped cocoordinate the large-scale Avian Phylogenomics Project, which collected and sequenced genomes of 48 bird species representing most currently classified orders to address a range of questions in phylogenomics and comparative genomics. The consortium was able to answer questions not previously possible with just a few genomes. This success spurred on the creation of a project to sequence the genomes of at least one individual of all extant ∼10,500 bird species. The initiation of this project has led us to consider what questions now impossible to answer could be answered with all genomes, and could drive new questions now unimaginable. These include the generation of a highly resolved family tree of extant species, genome-wide association studies across species to identify genetic substrates of many complex traits, redefinition of species and the species concept, reconstruction of the genomes of common ancestors, and generation of new computational tools to address these questions. Here I present visions for the future by posing and answering questions regarding what scientists could potentially do with available genomes of an entire vertebrate class.

The impact of modelling rate heterogeneity among sites on phylogenetic estimates of intraspecific evolutionary rates and timescales

Phylogenetic analyses of DNA sequence data can provide estimates of evolutionary rates and timescales. Nearly all phylogenetic methods rely on accurate models of nucleotide substitution. A key feature of molecular evolution is the heterogeneity of substitution rates among sites, which is often modelled using a discrete gamma distribution. A widely used derivative of this is the gamma-invariable mixture model, which assumes that a proportion of sites in the sequence are completely resistant to change, while substitution rates at the remaining sites are gamma-distributed. For data sampled at the intraspecific level, however, biological assumptions involved in the invariable-sites model are commonly violated. We examined the use of these models in analyses of five intraspecific data sets. We show that using 6-10 rate categories for the discrete gamma distribution of rates among sites is sufficient to provide a good approximation of the marginal likelihood. Increasing the number of gamma rate categories did not have a substantial effect on estimates of the substitution rate or coalescence time, unless rates varied strongly among sites in a non-gamma-distributed manner. The assumption of a proportion of invariable sites provided a better approximation of the asymptotic marginal likelihood when the number of gamma categories was small, but had minimal impact on estimates of rates and coalescence times. However, the estimated proportion of invariable sites was highly susceptible to changes in the number of gamma rate categories. The concurrent use of gamma and invariable-site models for intraspecific data is not biologically meaningful and has been challenged on statistical grounds; here we have found that the assumption of a proportion of invariable sites has no obvious impact on Bayesian estimates of rates and timescales from intraspecific data.

Among-site rate variation and its impact on phylogenetic analyses

Although several decades of study have revealed the ubiquity of variation of evolutionary rates among sites, reliable methods for studying rate variation were not developed until very recently. Early methods fit theoretical distributions to the numbers of changes at sites inferred by parsimony and substantially underestimate the rate variation. Recent analyses show that failure to account for rate variation can have drastic effects, leading to biased dating of speciation events, biased estimation of the transition:transversion rate ratio, and incorrect reconstruction of phylogenies.

Use of phylogenetics in the molecular epidemiology and evolutionary studies of viral infections

Since DNA sequencing techniques first became available almost 30 years ago, the amount of nucleic acid sequence data has increased enormously. Phylogenetics, which is widely applied to compare and analyze such data, is particularly useful for the analysis of genes from rapidly evolving viruses. It has been used extensively to describe the molecular epidemiology and transmission of the human immunodeficiency virus (HIV), the origins and subsequent evolution of the severe acute respiratory syndrome (SARS)-associated coronavirus (SCoV), and, more recently, the evolving epidemiology of avian influenza as well as seasonal and pandemic human influenza viruses. Recent advances in phylogenetic methods can infer more in-depth information about the patterns of virus emergence, adding to the conventional approaches in viral epidemiology. Examples of this information include estimations (with confidence limits) of the actual time of the origin of a new viral strain or its emergence in a new species, viral recombination and reassortment events, the rate of population size change in a viral epidemic, and how the virus spreads and evolves within a specific population and geographical region. Such sequence-derived information obtained from the phylogenetic tree can assist in the design and implementation of public health and therapeutic interventions. However, application of many of these advanced phylogenetic methods are currently limited to specialized phylogeneticists and statisticians, mainly because of their mathematical basis and their dependence on the use of a large number of computer programs. This review attempts to bridge this gap by presenting conceptual, technical, and practical aspects of applying phylogenetic methods in studies of influenza, HIV, and SCoV. It aims to provide, with minimal mathematics and statistics, a practical overview of how phylogenetic methods can be incorporated into virological studies by clinical and laboratory specialists.

The phylogeny of the social wasp subfamily Polistinae: evidence from microsatellite flanking sequences, mitochondrial COI sequence, and morphological characters

BACKGROUND

Social wasps in the subfamily Polistinae (Hymenoptera: Vespidae) have been important in studies of the evolution of sociality, kin selection, and within colony conflicts of interest. These studies have generally been conducted within species, because a resolved phylogeny among species is lacking. We used nuclear DNA microsatellite flanking sequences, mitochondrial COI sequence, and morphological characters to generate a phylogeny for the Polistinae (Hymenoptera) using 69 species.

RESULTS

Our phylogeny is largely concordant with previous phylogenies at higher levels, and is more resolved at the species level. Our results support the monophyly of the New World subgenera of Polistini, while the Old World subgenera are a paraphyletic group. All genera for which we had more than one exemplar were supported as monophyletic except Polybia which is not resolved, and may be paraphyletic.

CONCLUSION

The combination of DNA sequences from flanks of microsatellite repeats with mtCOI sequences and morphological characters proved to be useful characters establishing relationships among the different subgenera and species of the Polistini. This is the first detailed hypothesis for the species of this important group.

A maximum-likelihood analysis of eight phylogenetic markers in gallwasps (Hymenoptera: Cynipidae): implications for insect phylogenetic studies

We assessed the utility of eight DNA sequence markers (5.8S rDNA, 18S rDNA, 28S rDNA, ITS regions, long-wavelength opsin, elongation factor 1-alpha, cytochrome b, and cytochrome oxidase I) in reconstructing phylogenetic relationships at various levels of divergence in gallwasps (Hymenoptera: Cynipidae), using a set of eight exemplar taxa. We report sequence divergence values and saturation levels and compare phylogenetic results of these sequences analyzed both separately and combined to a well-corroborated morphological phylogeny. Likelihood ratio tests were used to find the best evolutionary model fitting each of the markers. The likelihood model best explaining the data is, for most loci, parameter rich, with strong A-T bias for mitochondrial loci and strong rate heterogeneity for the majority of loci. Our data suggest that 28S rDNA, elongation factor 1-alpha, and long-wavelength opsin may be potentially useful markers for the resolution of cynipid and other insect within-family-level divergences (circa 50-100 mya old), whereas mitochondrial loci and ITS regions are most useful for lower-level phylogenetics. In contrast, the 18S rDNA marker is likely to be useful for the resolution of above-family-level relationships.

Accounting for evolutionary rate variation among sequence sites consistently changes universal phylogenies deduced from rRNA and protein-coding genes

Phylogenetic analyses of gene and protein sequences have led to two major competing views of the universal phylogeny, the evolutionary tree relating the three kinds of living organisms, Bacteria, Archaea, and Eukarya. In the first scheme, called "the archaebacterial tree, " organisms of the same type are clustered together. In the second scenario, called "the eocyte tree," the archaeal phylum of Crenarchaeota is more closely related to eukaryotes than are other Archaea. A major property of the evolution of functional ribosomal and protein-encoding genes is that the rate of nucleotide and amino acid substitution varies across sequence sites. Here, using distance-based and maximum-likelihood methods, we show that universal phylogenies of ribosomal RNAs and RNA polymerases built by ignoring this variation are biased toward the archaebacterial tree because of attraction between long branches. In contrast, taking among-site rate variability into account gives support for the eocyte tree.

Microsporidia are related to Fungi: evidence from the largest subunit of RNA polymerase II and other proteins

We have determined complete gene sequences encoding the largest subunit of the RNA polymerase II (RBP1) from two Microsporidia, Vairimorpha necatrix and Nosema locustae. Phylogenetic analyses of these and other RPB1 sequences strongly support the notion that Microsporidia are not early-diverging eukaryotes but instead are specifically related to Fungi. Our reexamination of elongation factors EF-1alpha and EF-2 sequence data that had previously been taken as support for an early (Archezoan) divergence of these amitochondriate protists show such support to be weak and likely caused by artifacts in phylogenetic analyses. These EF data sets are, in fact, not inconsistent with a Microsporidia + Fungi relationship. In addition, we show that none of these proteins strongly support a deep divergence of Parabasalia and Metamonada, the other amitochondriate protist groups currently thought to compose early branches. Thus, the phylogenetic placement among eukaryotes for these protist taxa is in need of further critical examination.

Evolution on a volcanic conveyor belt: using phylogeographic reconstructions and K-Ar-based ages of the Hawaiian Islands to estimate molecular evolutionary rates

The Hawaiian Islands form as the Pacific Plate moves over a 'hot spot' in the earth's mantle where magma extrudes through the crust to build huge shield volcanos. The islands subside and erode as the plate carries them to the north-west, eventually to become coral atolls and seamounts. Thus islands are ordered linearly by age, with the oldest islands in the north-west (e.g. Kauai at 5.1 Ma) and the youngest in the south-east (e.g. Hawaii at 0.43 Ma). K-Ar estimates of the date of an island's formation provide a maximum age for the taxa inhabiting the island. These ages can be used to calibrate rates of molecular change under the following assumptions: (i) K-Ar dates are accurate; (ii) tree topologies show that derivation of taxa parallels the timing of island formation; (iii) populations do not colonize long after island emergence; (iv) the coalescent point for sister taxa does not greatly predate the formation of the colonized younger island; (v) saturation effects and (vi) among-lineage rate variation are minimal or correctable; and (vii) unbiased standard errors of distances and regressions can be estimated from multiple pairwise comparisons. We use the approach to obtain overall corrected rate calibrations for: (i) part of the mitochondrial cytochrome b gene in Hawaiian drepanidines (0.016 sequence divergence/Myr); (ii) the Yp1 gene in Hawaiian Drosophila (0.019/Myr Kambysellis et al. 1995); and (iii) parts of the mitochondrial 12S and 16S rRNA and tRNAval in Laupala crickets (0.024-0.102/Myr, Shaw 1996). We discuss the reliability of the estimates given the assumptions (i-vii) above and contrast the results with previous calibrations of Adh in Hawaiian Drosophila and chloroplast DNA in lobeliods.

Correlation of functional domains and rates of nucleotide substitution in cytochrome b

Distinguishing noise from signal presents a problem when DNA sequences are used for phylogeny reconstruction. Multiple substitutions at sites are a primary cause of noise and this is compounded by variation in substitution rates among sites. For protein-coding genes, one method used to determine if data are noisy is to assess levels of saturation of substitutions by codon position. However, this procedure may not be a fine enough filter for assessing noise. Variation in substitution rates may also be caused by constraints on change imposed by the function of the protein product. Using a structural model of the cytochrome b protein as a template, I divided cyt b sequence data for species within the avian family Falconidae (falcons and caracaras) into three functional domains. Saturation of substitutions of sequences within these regions was assessed graphically. This qualitative determination of saturation was then used to differentially weight phylogenetic analysis, resulting in an hypothesis congruent with existing cladistic analyses and traditional morphology. These results demonstrate that saturation of substitutions is correlated with functional regions of cytochrome b and that using this information improves phylogenetic inference.

Modeling nucleotide evolution: a heterogeneous rate analysis

A new model of molecular evolution is introduced that allows for heterogeneous rates across the sequence positions. The development of this model was motivated by two issues: first, a number of studies have shown that the positions in a DNA sequence evolve at different rates, and second, it has been shown that not accounting for this heterogeneity can lead to biased estimates of evolutionary parameters. The authors generalize the Markovian model of molecular evolution to allow for heterogeneous rates and explore some of the consequences of such a model. In particular, they quantify the biases incurred by incorrectly assuming an equal-rate model and consider what can be learned about evolutionary parameters under a heterogeneous model.

High divergence of reproductive tract proteins and their association with postzygotic reproductive isolation in Drosophila melanogaster and Drosophila virilis group species

The possible association between gonadal protein divergence and postzygotic reproductive isolation was investigated among species of the Drosophila melanogaster and D. virilis groups. Protein divergence was scored by high-resolution two-dimensional electrophoresis (2DE). Close to 500 protein spots from gonadal tissues (testis and ovary) and nongonadal tissues (malpighian tubules and brain) were analyzed and protein divergence was calculated based on presence vs absence. Both testis and ovary proteins showed higher divergence than nongonadal proteins, and also a highly significant positive correlation with postzygotic reproductive isolation but a weaker correlation with postzygotic reproductive isolation. Particularly, a positive and significant correlation was found between proteins expressed in the testis and postzygotic reproductive isolation among closely related species such as the within-phylad species in the D. virilis group. The high levels of male-reproductive-tract protein divergence between species might be associated with F1 hybrid male sterility among closely related species. If so, a lower level of ovary protein divergence should be expected on the basis that F1 female hybrids are fully fertile. However, this is not necessarily true if relatively few genes are responsible for the reproductive isolation observed between closely related species, as recent studies seem to suggest. We suggest that the faster rate of evolution of gonadal proteins in comparison to nongonadal proteins and the association of that rate with postzygotic reproductive isolation may be the result of episodic and/or sexual selection on male and female molecular traits.

Molecular evolution of interleukin-3

Chimpanzee, tamarin, and marmoset interleukin-3 (IL-3) genes were cloned, sequenced, and expressed. Western blot analysis demonstrated that functional genes were isolated. IL-3 sequences were compared with those of mouse, rat, rhesus monkey, gibbon, and man. Multiple alignment of the IL-3 coding regions showed that only a few regions had been conserved during mammalian evolution, which are likely associated with functional domains of the IL-3 protein. Substitution rates for the various lineages were calculated and the numbers of synonymous and nonsynonymous substitutions were estimated separately. Distance matrices of the IL-3 coding regions were used to construct phylogenetic trees which revealed large differences in IL-3 evolution rate as well as a more rapid substitution rate for rodents and a rate slowdown during hominoid evolution. Extremes were rhesus monkey IL-3, which accumulated few synonymous substitutions, and gibbon IL-3, which had almost exclusively synonymous substitutions. In rhesus monkey IL-3, nonsynonymous substitutions outnumbered synonymous substitutions, which could not be readily explained by a random process of substitutions. We assume that during evolution of IL-3, the majority of the amino acid replacements and the impaired interspecies functional cross-reactivity originate from selection mechanisms with the most likely selective force being the structure of the heterodimeric IL.3 cell-surface receptor. Insight into IL-3 architecture and structural analysis of the IL-3 receptor are needed to analyze the unusually fast evolution of IL-3 in more detail.

Maximum likelihood phylogenetic estimation from DNA sequences with variable rates over sites: approximate methods

Two approximate methods are proposed for maximum likelihood phylogenetic estimation, which allow variable rates of substitution across nucleotide sites. Three data sets with quite different characteristics were analyzed to examine empirically the performance of these methods. The first, called the "discrete gamma model," uses several categories of rates to approximate the gamma distribution, with equal probability for each category. The mean of each category is used to represent all the rates falling in the category. The performance of this method is found to be quite good, and four such categories appear to be sufficient to produce both an optimum, or near-optimum fit by the model to the data, and also an acceptable approximation to the continuous distribution. The second method, called "fixed-rates model", classifies sites into several classes according to their rates predicted assuming the star tree. Sites in different classes are then assumed to be evolving at these fixed rates when other tree topologies are evaluated. Analyses of the data sets suggest that this method can produce reasonable results, but it seems to share some properties of a least-squares pairwise comparison; for example, interior branch lengths in nonbest trees are often found to be zero. The computational requirements of the two methods are comparable to that of Felsenstein's (1981, J Mol Evol 17:368-376) model, which assumes a single rate for all the sites.

Substitution rate variation among sites in hypervariable region 1 of human mitochondrial DNA

More than an order of magnitude difference in substitution rate exists among sites within hypervariable region 1 of the control region of human mitochondrial DNA. A two-rate Poisson mixture and a negative binomial distribution are used to describe the distribution of the inferred number of changes per nucleotide site in this region. When three data sets are pooled, however, the two-rate model cannot explain the data. The negative binomial distribution always fits, suggesting that substitution rates are approximately gamma distributed among sites. Simulations presented here provide support for the use of a biased, yet commonly employed, method of examining rate variation. The use of parsimony in the method to infer the number of changes at each site introduces systematic errors into the analysis. These errors preclude an unbiased quantification of variation in substitution rate but make the method conservative overall. The method can be used to distinguish sites with highly elevated rates, and 29 such sites are identified in hypervariable region 1. Variation does not appear to be clustered within this region. Simulations show that biases in rates of substitution among nucleotides and non-uniform base composition can mimic the effects of variation in rate among sites. However, these factors contribute little to the levels of rate variation observed in hypervariable region 1.

Evolution of nest construction in swallows (Hirundinidae): a molecular phylogenetic perspective

Nest construction is more diverse in the Hirundinidae than in any other family of oscine birds. To explore the evolution of this diversity, we superimposed nest data on a DNA-hybridization phylogeny of 17 swallow species. Nest construction is tightly linked to the inferred evolutionary history. Burrowing appears to be the primitive nesting mode, and burrowing ancestors gave rise to cavity-adopting and mud-nesting clades. Obligate cavity adoption is mostly confined to a monophyletic clade in the New World, and the diversification of obligate nest adopters appears to be tied to the richness of forest habitats and recent active mountain building there. Construction of mud nests originated only once in the history of the group, and mud-nesters have diversified principally in Africa, where a drier climatic history has favored their mode of nesting. The use of pure mud to construct a hanging nest is unique among all birds, and we infer that mud nests have increased in complexity during evolution from simple mud cups to fully enclosed retort-shaped nests. This increased complexity appears to have been the critical precursor for the evolution of high-density colonial mud-nesters.

Nucleotide sequence of the COX1 gene in Kluyveromyces lactis mitochondrial DNA: evidence for recent horizontal transfer of a group II intron

The cytochrome oxidase subunit 1 gene (COX1) in K. lactis K8 mtDNA spans 8,826 bp and contains five exons (termed E1-E5) totalling 1,602 bp that show 88% nucleotide base matching and 91% amino acid homology to the equivalent gene in S. cerevisiae. The four introns (termed K1 cox1.1-1.4) contain open reading frames encoding proteins of 786, 333, 319 and 395 amino acids respectively that potentially encode maturase enzymes. The first intron belongs to group II whereas the remaining three are group I type B. Introns K1 cox1.1, 1.3, and 1.4 are found at identical locations to introns Sc cox1.2, 1.5 a, and 1.5 b respectively from S. cerevisiae. Horizontal transfer of an intron between recent progenitors of K. lactis and S. cerevisiae is suggested by the observation that K1 cox1.1 and Sc cox1.2 show 96% base matching. Sequence comparisons between K1 cox1.3/Sc cox1.5 a and K1 cox1.4/Sc cox1.5 b suggest that these introns are likely to have been present in the ancestral COX1 gene of these yeasts. Intron K1 cox1.2 is not found in S. cerevisiae and appears at an unique location in K. lactis. A feature of the DNA sequences of the group I introns K1 cox1.2, 1.3, and 1.4 is the presence of 11 GC-rich clusters inserted into both coding and noncoding regions. Immediately downstream of the COX1 gene is the ATPase subunit 8 gene (A8) that shows 82.6% base matching to its counterpart in S. cerevisiae mtDNA.