CALIBRATING DIVERGENCE TIMES ON SPECIES TREES VERSUS GENE TREES: IMPLICATIONS FOR SPECIATION HISTORY OF APHELOCOMA JAYS

Phylogenetic codivergence supports coevolution of mimetic Heliconius butterflies

The unpalatable and warning-patterned butterflies Heliconius erato and Heliconius melpomene provide the best studied example of mutualistic Müllerian mimicry, thought-but rarely demonstrated-to promote coevolution. Some of the strongest available evidence for coevolution comes from phylogenetic codivergence, the parallel divergence of ecologically associated lineages. Early evolutionary reconstructions suggested codivergence between mimetic populations of H. erato and H. melpomene, and this was initially hailed as one of the most striking known cases of coevolution. However, subsequent molecular phylogenetic analyses found discrepancies in phylogenetic branching patterns and timing (topological and temporal incongruence) that argued against codivergence. We present the first explicit cophylogenetic test of codivergence between mimetic populations of H. erato and H. melpomene, and re-examine the timing of these radiations. We find statistically significant topological congruence between multilocus coalescent population phylogenies of H. erato and H. melpomene. Cophylogenetic historical reconstructions support repeated codivergence of mimetic populations, from the base of the sampled radiations. Pairwise distance correlation tests, based on our coalescent analyses plus recently published AFLP and wing colour pattern gene data, also suggest that the phylogenies of H. erato and H. melpomene show significant topological congruence. Divergence time estimates, based on a Bayesian coalescent model, suggest that the evolutionary radiations of H. erato and H. melpomene occurred over the same time period, and are compatible with a series of temporally congruent codivergence events. Our results suggest that differences in within-species genetic divergence are the result of a greater overall effective population size for H. erato relative to H. melpomene and do not imply incongruence in the timing of their phylogenetic radiations. Repeated codivergence between Müllerian co-mimics, predicted to exert mutual selection pressures, strongly suggests coevolution. Our results therefore support a history of reciprocal coevolution between Müllerian co-mimics characterised by phylogenetic codivergence and parallel phenotypic change.

Bayesian phylogenetics with BEAUti and the BEAST 1.7

Computational evolutionary biology, statistical phylogenetics and coalescent-based population genetics are becoming increasingly central to the analysis and understanding of molecular sequence data. We present the Bayesian Evolutionary Analysis by Sampling Trees (BEAST) software package version 1.7, which implements a family of Markov chain Monte Carlo (MCMC) algorithms for Bayesian phylogenetic inference, divergence time dating, coalescent analysis, phylogeography and related molecular evolutionary analyses. This package includes an enhanced graphical user interface program called Bayesian Evolutionary Analysis Utility (BEAUti) that enables access to advanced models for molecular sequence and phenotypic trait evolution that were previously available to developers only. The package also provides new tools for visualizing and summarizing multispecies coalescent and phylogeographic analyses. BEAUti and BEAST 1.7 are open source under the GNU lesser general public license and available at http://beast-mcmc.googlecode.com and http://beast.bio.ed.ac.uk

Miocene dispersal drives island radiations in the palm tribe Trachycarpeae (Arecaceae)

The study of three island groups of the palm tribe Trachycarpeae (Arecaceae/Palmae) permits both the analysis of each independent radiation and comparisons across the tribe to address general processes that drive island diversification. Phylogenetic relationships of Trachycarpeae were inferred from three plastid and three low-copy nuclear genes. The incongruent topological position of Brahea in CISP5 was hypothesized to be caused by a gene duplication event and was addressed using uninode coding. The resulting phylogenetic trees were well-resolved and the genera were all highly supported except for Johannesteijsmannia and Serenoa. Divergence time analysis estimated the stem of the tribe to be approximately 86 Ma and the crown to be 38 Ma, indicating that significant extinction may have occurred along this branch. Historical biogeographic analysis suggested that Trachycarpeae are of southern North American, Central American, or Caribbean origin and supports previous hypotheses of a Laurasian origin. The biogeography and disjunctions within the tribe were interpreted with respect to divergence times, the fossil record, and geological factors such as the formation of the Greater Antilles--Aves Ridge, the Bering and the North Atlantic land bridges, tectonic movement in Southeast Asia, climatic shifts between the Eocene and Pliocene, and volcanism in the Pacific basin. In considering the three major island radiations within Trachycarpeae, Miocene dispersal appears to have been the driving force in allopatric speciation and is highlighted here as an emerging pattern across the tree of life.

A multilocus phylogeny of Asian noodlefishes Salangidae (Teleostei: Osmeriformes) with a revised classification of the family

A group of small and transparent Asian noodlefishes (Osmeriformes: Salangidae) are commercially important fishery species, however, interrelationships among these fishes remain unresolved in previous studies using mitochondrial markers. We re-examine phylogenetic relationships of Salangidae by including complete taxon sampling, based on seven nuclear loci and one mitochondrial gene using a multilocus coalescence-based species-tree method. Our results show a well-resolved phylogeny of Salangidae that does not agree with previous hypotheses. The topology test suggests that our hypothesis represents the most likely phylogeny. Using the inferred species-tree as criterion, we recombine the rank of subfamilies and genera in the Salangidae, and erect a new genus Neosalangichthys. Our revised classification of Salangidae is well supported by reinterpreting previously proposed diagnostic characters. Finally, re-defined synapomorphic characters are used to erect a key to the genera of Salangidae.

Speciation in an avian complex endemic to the mountains of Middle America (Ergaticus, Aves: Parulidae)

The implementation of the phylogeographic approach for the study of biodiversity is critical in poorly sampled regions like the montane systems of Middle America, as complex evolutionary histories often result in the presence of independent lineages not properly considered by traditional taxonomy. Herein we sequenced 2370 bp of mtDNA (ND2, cyt b and ATPase) from 81 individuals of Ergaticus, a complex of birds endemic to the montane forests of Middle America. Although current taxonomy recognizes two species, the results reveal considerable genetic structure with the presence of four mtDNA lineages. Two of these lineages within Ergaticus ruber evidence the need of a revaluation of the species limits for this taxon. The general phylogeographic pattern can be explained as a consequence of relative isolation of the populations in different mountain ranges separated by low elevation barriers. Most population groups did not show signals of demographic expansion with the exception of the one corresponding to clade 1. The divergence time estimates point to the Pleistocene as an important time period for the diversification of this complex.

Recent speciation in three closely related sympatric specialists: inferences using multi-locus sequence, post-mating isolation and endosymbiont data

Shifting between unrelated host plants is relatively rare for phytophagous insects, and distinct host specificity may play crucial roles in reproductive isolation. However, the isolation status and the relationship between parental divergence and post-mating isolation among closely related sympatric specialists are still poorly understood. Here, multi-locus sequence were used to estimate the relationship among three host plant–specific closely related flea beetles, Altica cirsicola, A. fragariae and A. viridicyanea (abbreviated as AC, AF and AV respectively). The tree topologies were inconsistent using different gene or different combinations of gene fragments. The relationship of AF+(AC+AV) was supported, however, by both gene tree and species tree based on concatenated data. Post-mating reproductive data on the results of crossing these three species are best interpreted in the light of a well established phylogeny. Nuclear-induced but not Wolbachia-induced unidirectional cytoplasmic incompatibility, which was detected in AC-AF and AF-AV but not in AC-AV, may also suggest more close genetic affinity between AC and AV. Prevalence of Wolbachia in these three beetles, and the endosymbiont in most individuals of AV and AC sharing a same wsp haplotype may give another evidence of AF+(AC+AV). Our study also suggested that these three flea beetles diverged in a relative short time (0.94 My), which may be the result of shifting between unrelated host plants and distinct host specificity. Incomplete post-mating isolation while almost complete lineage sorting indicated that effective pre-mating isolation among these three species should have evolved.

Gone with the plate: the opening of the Western Mediterranean basin drove the diversification of ground-dweller spiders

BACKGROUND

The major islands of the Western Mediterranean--Corsica, Sardinia, and the Balearic Islands--are continental terrenes that drifted towards their present day location following a retreat from their original position on the eastern Iberian Peninsula about 30 million years ago. Several studies have taken advantage of this well-dated geological scenario to calibrate molecular rates in species for which distributions seemed to match this tectonic event. Nevertheless, the use of external calibration points has revealed that most of the present-day fauna on these islands post-dated the opening of the western Mediterranean basin. In this study, we use sequence information of the cox1, nad1, 16S, L1, and 12S mitochondrial genes and the 18S, 28S, and h3 nuclear genes, along with relaxed clock models and a combination of biogeographic and fossil external calibration points, to test alternative historical scenarios of the evolutionary history of the ground-dweller spider genus Parachtes (Dysderidae), which is endemic to the region.

RESULTS

We analyse 49 specimens representing populations of most Parachtes species and close relatives. Our results reveal that both the sequence of species formation in Parachtes and the estimated divergence times match the geochronological sequence of separation of the main islands, suggesting that the diversification of the group was driven by Tertiary plate tectonics. In addition, the confirmation that Parachtes diversification matches well-dated geological events provides a model framework to infer substitution rates of molecular markers. Divergence rates estimates ranged from 3.5% My(-1) (nad1) to 0.12% My(-1) (28S), and the average divergence rate for the mitochondrial genes was 2.25% My(-1), very close to the "standard" arthropod mitochondrial rate (2.3% My(-1)).

CONCLUSIONS

Our study provides the first unequivocal evidence of terrestrial endemic fauna of the major western Mediterranean islands, whose origin can be traced back to the Oligocene separation of these islands from the continent. Moreover, our study provides useful information on the divergence rate estimates of the most commonly used genes for phylogenetic inference in non-model arthropods.

Multilocus resolution of phylogeny and timescale in the extant adaptive radiation of Hawaiian honeycreepers

Evolutionary theory has gained tremendous insight from studies of adaptive radiations. High rates of speciation, morphological divergence, and hybridization, combined with low sequence variability, however, have prevented phylogenetic reconstruction for many radiations. The Hawaiian honeycreepers are an exceptional adaptive radiation, with high phenotypic diversity and speciation that occurred within the geologically constrained setting of the Hawaiian Islands. Here we analyze a new data set of 13 nuclear loci and pyrosequencing of mitochondrial genomes that resolves the Hawaiian honeycreeper phylogeny. We show that they are a sister taxon to Eurasian rosefinches (Carpodacus) and probably came to Hawaii from Asia. We use island ages to calibrate DNA substitution rates, which vary substantially among gene regions, and calculate divergence times, showing that the radiation began roughly when the oldest of the current large Hawaiian Islands (Kauai and Niihau) formed, ~5.7 million years ago (mya). We show that most of the lineages that gave rise to distinctive morphologies diverged after Oahu emerged (4.0-3.7 mya) but before the formation of Maui and adjacent islands (2.4-1.9 mya). Thus, the formation of Oahu, and subsequent cycles of colonization and speciation between Kauai and Oahu, played key roles in generating the morphological diversity of the extant honeycreepers.

Isolation-driven divergence: speciation in a widespread North American songbird (Aves: Certhiidae)

Lineage, or true 'species', trees may differ from gene trees because of stochastic processes in molecular evolution leading to gene-tree heterogeneity. Problems with inferring species trees because of excessive incomplete lineage sorting may be exacerbated in lineages with rapid diversification or recent divergences necessitating the use of multiple loci and individuals. Many recent multilocus studies that investigate divergence times identify lineage splitting to be more recent than single-locus studies, forcing the revision of biogeographic scenarios driving divergence. Here, we use 21 nuclear loci from regional populations to re-evaluate hypotheses identified in an mtDNA phylogeographic study of the Brown Creeper (Certhia americana), as well as identify processes driving divergence. Nuclear phylogeographic analyses identified hierarchical genetic structure, supporting a basal split at approximately 32°N latitude, splitting northern and southern populations, with mixed patterns of genealogical concordance and discordance between data sets within the major lineages. Coalescent-based analyses identify isolation, with little to no gene flow, as the primary driver of divergence between lineages. Recent isolation appears to have caused genetic bottlenecks in populations in the Sierra Madre Oriental and coastal mountain ranges of California, which may be targets for conservation concerns.

Rate variation and estimation of divergence times using strict and relaxed clocks

BACKGROUND

Understanding causes of biological diversity may be greatly enhanced by knowledge of divergence times. Strict and relaxed clock models are used in Bayesian estimation of divergence times. We examined whether: i) strict clock models are generally more appropriate in shallow phylogenies where rate variation is expected to be low, ii) the likelihood ratio test of the clock (LRT) reliably informs which model is appropriate for dating divergence times. Strict and relaxed models were used to analyse sequences simulated under different levels of rate variation. Published shallow phylogenies (Black bass, Primate-sucking lice, Podarcis lizards, Gallotiinae lizards, and Caprinae mammals) were also analysed to determine natural levels of rate variation relative to the performance of the different models.

RESULTS

Strict clock analyses performed well on data simulated under the independent rates model when the standard deviation of log rate on branches, σ, was low (≤ 0.1), but were inappropriate when σ>0.1 (95% of rates fall within 0.0082-0.0121 subs/site/Ma when σ = 0.1, for a mean rate of 0.01). The independent rates relaxed clock model performed well at all levels of rate variation, although posterior intervals on times were significantly wider than for the strict clock. The strict clock is therefore superior when rate variation is low. The performance of a correlated rates relaxed clock model was similar to the strict clock. Increased numbers of independent loci led to slightly narrower posteriors under the relaxed clock while older root ages provided proportionately narrower posteriors. The LRT had low power for σ = 0.01-0.1, but high power for σ = 0.5-2.0. Posterior means of σ2 were useful for assessing rate variation in published datasets. Estimates of natural levels of rate variation ranged from 0.05-3.38 for different partitions. Differences in divergence times between relaxed and strict clock analyses were greater in two datasets with higher σ2 for one or more partitions, supporting the simulation results.

CONCLUSIONS

The strict clock can be superior for trees with shallow roots because of low levels of rate variation between branches. The LRT allows robust assessment of suitability of the clock model as does examination of posteriors on σ2.

Refugial isolation and divergence in the Narrowheaded Gartersnake species complex (Thamnophis rufipunctatus) as revealed by multilocus DNA sequence data

Glacial-interglacial cycles of the Pleistocene are hypothesized as one of the foremost contributors to biological diversification. This is especially true for cold-adapted montane species, where range shifts have had a pronounced effect on population-level divergence. Gartersnakes of the Thamnophis rufipunctatus species complex are restricted to cold headwater streams in the highlands of the Sierra Madre Occidental and southwestern USA. We used coalescent and multilocus phylogenetic approaches to test whether genetic diversification of this montane-restricted species complex is consistent with two prevailing models of range fluctuation for species affected by Pleistocene climate changes. Our concatenated nuDNA and multilocus species analyses recovered evidence for the persistence of multiple lineages that are restricted geographically, despite a mtDNA signature consistent with either more recent connectivity (and introgression) or recent expansion (and incomplete lineage sorting). Divergence times estimated using a relaxed molecular clock and fossil calibrations fall within the Late Pleistocene, and zero gene flow scenarios among current geographically isolated lineages could not be rejected. These results suggest that increased climate shifts in the Late Pleistocene have driven diversification and current range retraction patterns and that the differences between markers reflect the stochasticity of gene lineages (i.e. ancestral polymorphism) rather than gene flow and introgression. These results have important implications for the conservation of T. rufipunctatus (sensu novo), which is restricted to two drainage systems in the southwestern US and has undergone a recent and dramatic decline.

Evolutionary patterns of diversification in the Andean hummingbird genus Adelomyia

The patterns of genetic diversity and morphological variation are of central importance in understanding the evolutionary process that drive diversification. We use molecular, morphological, and ecological data to explore the influence of geography and ecology in promoting speciation in the widespread Andean hummingbird genus Adelomyia. Six monophyletic clades were recovered which show distributional limits at well-defined geographic barriers. Percentage sequence divergence ranged between 5.8% and 8.2% between phylogroups separated by large (>4000 km) and small (<50 km) distances respectively, suggesting that geographic isolation may be influential at very different scales. We show that morphological traits in independent phylogroups are more related to environmental heterogeneity than to geographic barriers. We provide a molecular reconstruction of relationships within Adelomyia and recommend its use in future comparative studies of historical biogeography and diversification in the Andes.

Estimating phylogenetic relationships despite discordant gene trees across loci: the species tree of a diverse species group of feather mites (Acari: Proctophyllodidae)

With the increased availability of multilocus sequence data, the lack of concordance of gene trees estimated for independent loci has focused attention on both the biological processes producing the discord and the methodologies used to estimate phylogenetic relationships. What has emerged is a suite of new analytical tools for phylogenetic inference--species tree approaches. In contrast to traditional phylogenetic methods that are stymied by the idiosyncrasies of gene trees, approaches for estimating species trees explicitly take into account the cause of discord among loci and, in the process, provides a direct estimate of phylogenetic history (i.e. the history of species divergence, not divergence of specific loci). We illustrate the utility of species tree estimates with an analysis of a diverse group of feather mites, the pinnatus species group (genus Proctophyllodes). Discord among four sequenced nuclear loci is consistent with theoretical expectations, given the short time separating speciation events (as evident by short internodes relative to terminal branch lengths in the trees). Nevertheless, many of the relationships are well resolved in a Bayesian estimate of the species tree; the analysis also highlights ambiguous aspects of the phylogeny that require additional loci. The broad utility of species tree approaches is discussed, and specifically, their application to groups with high speciation rates--a history of diversification with particular prevalence in host/parasite systems where species interactions can drive rapid diversification.