Evidence for time dependency of molecular rate estimates

A method to correct for the effects of purifying selection on genealogical inference

Accurate reconstruction of the divergence times among individuals is an essential step toward inferring population parameters from genetic data. However, our ability to reconstruct accurate genealogies is often thwarted by the evolutionary forces we hope to detect, most prominently natural selection. Here, I demonstrate that purifying selection acting at many linked sites can systematically bias current methods of genealogical reconstruction, and I present a new method that corrects for this bias by allowing a class of sites to have a time-dependent rate. The parameters influencing the time dependency can be estimated from the data, allowing for a general method to detect the presence of selected sites and correcting for their distortion of the apparent mutation rate. The method works well under a variety of scenarios, including gamma-distributed selection coefficients as well as entirely neutral evolution. I also compare the performance of the new method to relaxed clock models, and I demonstrate the method on a data set from the mitochondrion of the North Atlantic whale-"louse" Cyamus ovalis.

Molecular evidence of the survival of subterranean amphipods (Arthropoda) during Ice Age underneath glaciers in Iceland

A Two endemic groundwater arthropod crustacean species, Crangonyx islandicus and Crymostygius thingvallensis, were recently discovered on the mid-Atlantic volcanic island of Iceland. The extent of morphological differences from closest relatives, endemism, along with the geographic isolation of Iceland and its complete coverage by glaciers 21,000 years ago, suggests that these two species have survived glaciation periods in sub-glacial refugia. Here we provide strong support for this hypothesis by an analysis of mitochondrial genetic variation within Crangonyx islandicus. Our results show that the species is divided into several distinct monophyletic groups that are found along the volcanic zone in Iceland, which have been separated by 0.5 to around 5 million years. The genetic divergence between groups reflects geographic distances between sampling sites, indicating that divergence occurred after the colonization of Iceland. The genetic patterns, as well as the dependency of genetic variation on distances from the tectonic plate boundary and altitude, points to recent expansion from several refugia within Iceland. This presents the first genetic evidence of multicellular organisms as complex as crustacean amphipods which have survived glaciations beneath an ice sheet. This survival may be explained by geothermal heat linked to volcanic activities, which may have maintained favourable habitats in fissures along the tectonic plate boundary in Iceland during glaciations.

Phylogeographic history of grey wolves in Europe

Background

While it is generally accepted that patterns of intra-specific genetic differentiation are substantially affected by glacial history, population genetic processes occurring during Pleistocene glaciations are still poorly understood. In this study, we address the question of the genetic consequences of Pleistocene glaciations for European grey wolves. Combining our data with data from published studies, we analysed phylogenetic relationships and geographic distribution of mitochondrial DNA haplotypes for 947 contemporary European wolves. We also compared the contemporary wolf sequences with published sequences of 24 ancient European wolves.

Results

We found that haplotypes representing two haplogroups, 1 and 2, overlap geographically, but substantially differ in frequency between populations from south-western and eastern Europe. A comparison between haplotypes from Europe and other continents showed that both haplogroups are spread throughout Eurasia, while only haplogroup 1 occurs in contemporary North American wolves. All ancient wolf samples from western Europe that dated from between 44,000 and 1,200 years B.P. belonged to haplogroup 2, suggesting the long-term predominance of this haplogroup in this region. Moreover, a comparison of current and past frequencies and distributions of the two haplogroups in Europe suggested that haplogroup 2 became outnumbered by haplogroup 1 during the last several thousand years.

Conclusions

Parallel haplogroup replacement, with haplogroup 2 being totally replaced by haplogroup 1, has been reported for North American grey wolves. Taking into account the similarity of diets reported for the late Pleistocene wolves from Europe and North America, the correspondence between these haplogroup frequency changes may suggest that they were associated with ecological changes occurring after the Last Glacial Maximum.

How accurate and robust are the phylogenetic estimates of Austronesian language relationships?

We recently used computational phylogenetic methods on lexical data to test between two scenarios for the peopling of the Pacific. Our analyses of lexical data supported a pulse-pause scenario of Pacific settlement in which the Austronesian speakers originated in Taiwan around 5,200 years ago and rapidly spread through the Pacific in a series of expansion pulses and settlement pauses. We claimed that there was high congruence between traditional language subgroups and those observed in the language phylogenies, and that the estimated age of the Austronesian expansion at 5,200 years ago was consistent with the archaeological evidence. However, the congruence between the language phylogenies and the evidence from historical linguistics was not quantitatively assessed using tree comparison metrics. The robustness of the divergence time estimates to different calibration points was also not investigated exhaustively. Here we address these limitations by using a systematic tree comparison metric to calculate the similarity between the Bayesian phylogenetic trees and the subgroups proposed by historical linguistics, and by re-estimating the age of the Austronesian expansion using only the most robust calibrations. The results show that the Austronesian language phylogenies are highly congruent with the traditional subgroupings, and the date estimates are robust even when calculated using a restricted set of historical calibrations.

Revisiting the insect mitochondrial molecular clock: the mid-Aegean trench calibration

Phylogenetic trees in insects are frequently dated by applying a "standard" mitochondrial DNA (mtDNA) clock estimated at 2.3% My(-1), but despite its wide use reliable calibration points have been lacking. Here, we used a well-established biogeographic barrier, the mid-Aegean trench separating the western and eastern Aegean archipelago, to estimate substitution rates in tenebrionid beetles. Cytochrome oxidase I (cox1) for six codistributed genera across 28 islands (444 individuals) on both sides of the mid-Aegean trench revealed 60 independently coalescing entities delimited with a mixed Yule-coalescent model. One representative per entity was used for phylogenetic analysis of mitochondrial (cox1, 16S rRNA) and nuclear (Mp20, 28S rRNA) genes. Six nodes marked geographically congruent east-west splits whose separation was largely contemporaneous and likely to reflect the formation of the mid-Aegean trench at 9-12 Mya. Based on these "known" dates, a divergence rate of 3.54% My(-1) for the cox1 gene (2.69% when combined with the 16S rRNA gene) was obtained under the preferred partitioning scheme and substitution model selected using Bayes factors. An extensive survey suggests that discrepancies in mtDNA substitution rates in the entomological literature can be attributed to the use of different substitution models, the use of different mitochondrial gene regions, mixing of intraspecific with interspecific data, and not accounting for variance in coalescent times or postseparation gene flow. Different treatments of these factors in the literature confound estimates of mtDNA substitution rates in opposing directions and obscure lineage-specific differences in rates when comparing data from various sources.

Do complex population histories drive higher estimates of substitution rate in phylogenetic reconstructions?

Our curiosity about biodiversity compels us to reconstruct the evolutionary past of species. Molecular evolutionary theory now allows parameterization of mathematically sophisticated and detailed models of DNA evolution, which have resulted in a wealth of phylogenetic histories. But reconstructing how species and population histories have played out is critically dependent on the assumptions we make, such as the clock-like accumulation of genetic differences over time and the rate of accumulation of such differences. An important stumbling block in the reconstruction of evolutionary history has been the discordance in estimates of substitution rate between phylogenetic and pedigree-based studies. Ancient genetic data recovered directly from the past are intermediate in time scale between phylogenetics-based and pedigree-based calibrations of substitution rate. Recent analyses of such ancient genetic data suggest that substitution rates are closer to the higher, pedigree-based estimates. In this issue, Navascués & Emerson (2009) model genetic data from contemporary and ancient populations that deviate from a simple demographic history (including changes in population size and structure) using serial coalescent simulations. Furthermore, they show that when these data are used for calibration, we are likely to arrive at upwardly biased estimates of mutation rate.

Insights into the evolutionary history of an emerging livestock pathogen: porcine circovirus 2

Porcine circovirus 2 (PCV2) is the primary etiological agent of postweaning multisystemic wasting syndrome (PMWS), one of the most economically important emerging swine diseases worldwide. Virulent PCV2 was first identified following nearly simultaneous outbreaks of PMWS in North America and Europe in the 1990s and has since achieved global distribution. However, the processes responsible for the emergence and spread of PCV2 remain poorly understood. Here, phylogenetic and cophylogenetic inferences were utilized to address key questions on the time scale, processes, and geographic diffusion of emerging PCV2. The results of these analyses suggest that the two genotypes of PCV2 (PCV2a and PCV2b) are likely to have emerged from a common ancestor approximately 100 years ago and have been on independent evolutionary trajectories since that time, despite cocirculating in the same host species and geographic regions. The patterns of geographic movement of PCV2 that we recovered appear to mimic those of the global pig trade and suggest that the movement of asymptomatic animals is likely to have facilitated the rapid spread of virulent PCV2 around the globe. We further estimated the rate of nucleotide substitution for PCV2 to be on the order of 1.2 x 10(-3) substitutions/site/year, the highest yet recorded for a single-stranded DNA virus. This high rate of evolution may allow PCV2 to maintain evolutionary dynamics closer to those of single-stranded RNA viruses than to those of double-stranded DNA viruses, further facilitating the rapid emergence of PCV2 worldwide.

Molecular insights into the colonization and chromosomal diversification of Madeiran house mice

The colonization history of Madeiran house mice was investigated by analysing the complete mitochondrial (mt) D-loop sequences of 156 mice from the island of Madeira and mainland Portugal, extending on previous studies. The numbers of mtDNA haplotypes from Madeira and mainland Portugal were substantially increased (17 and 14 new haplotypes respectively), and phylogenetic analysis confirmed the previously reported link between the Madeiran archipelago and northern Europe. Sequence analysis revealed the presence of four mtDNA lineages in mainland Portugal, of which one was particularly common and widespread (termed the 'Portugal Main Clade'). There was no support for population bottlenecks during the formation of the six Robertsonian chromosome races on the island of Madeira, and D-loop sequence variation was not found to be structured according to karyotype. The colonization time of the Madeiran archipelago by Mus musculus domesticus was approached using two molecular dating methods (mismatch distribution and Bayesian skyline plot). Time estimates based on D-loop sequence variation at mainland sites (including previously published data from France and Turkey) were evaluated in the context of the zooarchaeological record of M. m. domesticus. A range of values for mutation rate (mu) and number of mouse generations per year was considered in these analyses because of the uncertainty surrounding these two parameters. The colonization of Portugal and Madeira by house mice is discussed in the context of the best-supported parameter values. In keeping with recent studies, our results suggest that mutation rate estimates based on interspecific divergence lead to gross overestimates concerning the timing of recent within-species events.

Quantitative prediction of molecular clock and ka/ks at short timescales

Recent empirical studies of taxa including humans, fish, and birds have shown elevated rates of molecular evolution between species that diverged recently. Using the Moran model, we calculate expected divergence as a function of time. Our findings suggest that the observed phenomenon of elevated rates at short timescales is consistent with standard population genetics theory. The apparent acceleration of the molecular clock at short timescales can be explained by segregating polymorphisms present at the time of the ancestral population, both neutral and slightly deleterious, and not newly arising slightly deleterious mutations as has been previously hypothesized. Our work also suggests that the duration of the rate elevation depends on the effective population size, providing a method to correct time estimates of recent divergence events. Our model concords with estimates of divergence obtained from African cichlid fish and humans. As an additional application of our model, we calculate that K(a)/K(s) is elevated within a population before decaying slowly to its long-term value. Similar to the molecular clock, the duration and magnitude of K(a)/K(s) elevation depend on the effective population size. Unlike the molecular clock, however, K(a)/K(s) elevation is caused by newly arising slightly deleterious mutations. This elevation, although not as severe in magnitude as had been previously predicted in models neglecting ancestral polymorphism, persists slightly longer.

Bayesian hypothesis testing supports long-distance Pleistocene migrations in a European high mountain plant (Androsace vitaliana, Primulaceae)

Colonization of the south-western European mountain ranges is suggested to have predominantly progressed from the Iberian Peninsula eastwards, but this hypothesis has never been tested in a statistical framework. Here, we test this hypothesis using Androsace vitaliana, a high elevation species with eight mostly allopatric subspecies, which is widely but disjunctly distributed across all major south-western European mountain ranges. To this end, we use plastid and nuclear sequence data as well as fingerprint (amplified fragment length polymorphisms) data and employ Bayesian methods, which allow co-estimation of genealogy and divergence times using explicit demographic models, as well as hypothesis testing via Bayes factors. Irrespective of the ambiguity concerning where A. vitaliana started to diversify -- both the Alps and the mountain ranges of the Iberian Peninsula outside the Pyrenees were possible -- colonization routes were not simply unidirectional, but involved Pleistocene connections between the Alps and mountain ranges of the Iberian Peninsula bypassing the interjacent Pyrenees via long-distance dispersal. In contrast, the species' post-glacial history is shaped by regional gene pool homogenization resulting in the genetic pattern showing good congruence with geographical proximity in agreement with a vicariance model, but only partly supporting current taxonomy.

The perils of plenty: what are we going to do with all these genes?

This new century's biology promises more of everything--more genes, more organisms, more species and, in short, more data. The flood of data challenges us to find better and quicker ways to summarize and analyse. Here, we present preliminary results and proofs of concept from three of our research projects that are motivated by our search for solutions to the perils of plenty. First, we discuss how models of evolution can accommodate change to better reflect the dynamics of sequence diversity, particularly when it is becoming a lot easier to obtain sequences at different times and across intervals where the probability of new mutations contributing to this diversity is high. Second, we describe our work on the use of a single locus for species delimitation; this research targets the new DNA-barcoding approach that aims to catalogue the entirety of life. We have developed a single-locus test based on the coalescent that tests the null hypothesis of panmixis. Finally, we discuss new sequencing technologies, the types of data available and the efficacy of alignment-free methods to estimate pairwise distances for phylogenetic analyses.

Correcting for purifying selection: an improved human mitochondrial molecular clock

There is currently no calibration available for the whole human mtDNA genome, incorporating both coding and control regions. Furthermore, as several authors have pointed out recently, linear molecular clocks that incorporate selectable characters are in any case problematic. We here confirm a modest effect of purifying selection on the mtDNA coding region and propose an improved molecular clock for dating human mtDNA, based on a worldwide phylogeny of > 2000 complete mtDNA genomes and calibrating against recent evidence for the divergence time of humans and chimpanzees. We focus on a time-dependent mutation rate based on the entire mtDNA genome and supported by a neutral clock based on synonymous mutations alone. We show that the corrected rate is further corroborated by archaeological dating for the settlement of the Canary Islands and Remote Oceania and also, given certain phylogeographic assumptions, by the timing of the first modern human settlement of Europe and resettlement after the Last Glacial Maximum. The corrected rate yields an age of modern human expansion in the Americas at approximately 15 kya that-unlike the uncorrected clock-matches the archaeological evidence, but continues to indicate an out-of-Africa dispersal at around 55-70 kya, 5-20 ky before any clear archaeological record, suggesting the need for archaeological research efforts focusing on this time window. We also present improved rates for the mtDNA control region, and the first comprehensive estimates of positional mutation rates for human mtDNA, which are essential for defining mutation models in phylogenetic analyses.

The influence of climatic oscillations during the Quaternary Era on the genetic structure of Asian black bears in Japan

The Asian black bear (Ursus thibetanus) inhabits two of the main islands, Honshu and Shikoku, in Japan. To determine how climatic oscillations during the Quaternary Era affected the genetic structure of the black bear populations in Japan, we examined their phylogeographic relationships and compared their genetic structure. We analysed an approximately 700-bp sequence in the D-loop region of mitochondrial DNA collected from 589 bears in this study with 108 bears from a previous study. We observed a total of 57 haplotypes and categorized them into three clusters (Eastern, Western and Southern) based on the spatial distribution of the haplotypes. All but 2 of the 41 haplotypes in the Eastern cluster were distributed locally. Genetic diversity was generally low in northern Japan and high in central Japan. Demographic tests rejected the expansion model in northern populations. Haplotypes of the Western and Southern clusters were unique to local populations. We conclude that the extant genetic structure of the Asian black bear populations arose as follows: first, populations became small and genetic drift decreased genetic diversity in the northern area during the last glacial period, whereas large continuous populations existed in the southern part of central Japan. These patterns were essentially maintained until the present time. In western and southern Japan, the effects of climatic oscillations were smaller, and thus, local structure was maintained.

Language phylogenies reveal expansion pulses and pauses in Pacific settlement

Debates about human prehistory often center on the role that population expansions play in shaping biological and cultural diversity. Hypotheses on the origin of the Austronesian settlers of the Pacific are divided between a recent "pulse-pause" expansion from Taiwan and an older "slow-boat" diffusion from Wallacea. We used lexical data and Bayesian phylogenetic methods to construct a phylogeny of 400 languages. In agreement with the pulse-pause scenario, the language trees place the Austronesian origin in Taiwan approximately 5230 years ago and reveal a series of settlement pauses and expansion pulses linked to technological and social innovations. These results are robust to assumptions about the rooting and calibration of the trees and demonstrate the combined power of linguistic scholarship, database technologies, and computational phylogenetic methods for resolving questions about human prehistory.

Island biogeography of Galápagos lava lizards (Tropiduridae: Microlophus): species diversity and colonization of the archipelago

The "lava lizards" (Microlophus) are distributed throughout the Galápagos Archipelago, and consist of radiations derived from two independent colonizations. The "Eastern Radiation" includes M. bivittatus and M. habeli endemic to San Cristobal and Marchena Islands. The "Western Radiation" includes five to seven historically recognized species distributed across almost the entire Archipelago. We combine dense geographic sampling and multilocus sequence data to estimate a phylogenetic hypothesis for the Western Radiation, to delimit species boundaries in this radiation, and to estimate a time frame for colonization events. Our phylogenetic hypothesis rejects two earlier topologies for the Western Radiation and paraphyly of M. albemarlensis, while providing strong support for single colonizations on each island. The colonization history implied by our phylogeny is consistent with general expectations of an east-to-west route predicted by the putative age of island groups, and prevailing ocean currents in the Archipelago. Additionally, combined evidence suggests that M. indefatigabilis from Santa Fe should be recognized as a full species. Finally, molecular divergence estimates suggest that the two colonization events likely occurred on the oldest existing islands, and the Western Radiation represents a recent radiation that, in most cases, has produced species that are considerably younger than the islands they inhabit.

A recalibrated molecular clock and independent origins for the cholera pandemic clones

Cholera, caused by Vibrio cholerae, erupted globally from South Asia in 7 pandemics, but there were also local outbreaks between the 6(th) (1899-1923) and 7(th) (1961-present) pandemics. All the above are serotype O1, whereas environmental or invertebrate isolates are antigenically diverse. The pre 7th pandemic isolates mentioned above, and other minor pathogenic clones, are related to the 7(th) pandemic clone, while the 6(th) pandemic clone is in the same lineage but more distantly related, and non-pathogenic isolates show no clonal structure. To understand the origins and relationships of the pandemic clones, we sequenced the genomes of a 1937 prepandemic strain and a 6(th) pandemic isolate, and compared them with the published 7(th) pandemic genome. We distinguished mutational and recombinational events, and allocated these and other events, to specific branches in the evolutionary tree. There were more mutational than recombinational events, but more genes, and 44 times more base pairs, changed by recombination. We used the mutational single-nucleotide polymorphisms and known isolation dates of the prepandemic and 7(th) pandemic isolates to estimate the mutation rate, and found it to be 100 fold higher than usually assumed. We then used this to estimate the divergence date of the 6(th) and 7(th) pandemic clones to be about 1880. While there is a large margin of error, this is far more realistic than the 10,000-50,000 years ago estimated using the usual assumptions. We conclude that the 2 pandemic clones gained pandemic potential independently, and overall there were 29 insertions or deletions of one or more genes. There were also substantial changes in the major integron, attributed to gain of individual cassettes including copying from within, or loss of blocks of cassettes. The approaches used open up new avenues for analysing the origin and history of other important pathogens.

Mutation and evolutionary rates in adélie penguins from the antarctic

Precise estimations of molecular rates are fundamental to our understanding of the processes of evolution. In principle, mutation and evolutionary rates for neutral regions of the same species are expected to be equal. However, a number of recent studies have shown that mutation rates estimated from pedigree material are much faster than evolutionary rates measured over longer time periods. To resolve this apparent contradiction, we have examined the hypervariable region (HVR I) of the mitochondrial genome using families of Adélie penguins (Pygoscelis adeliae) from the Antarctic. We sequenced 344 bps of the HVR I from penguins comprising 508 families with 915 chicks, together with both their parents. All of the 62 germline heteroplasmies that we detected in mothers were also detected in their offspring, consistent with maternal inheritance. These data give an estimated mutation rate (micro) of 0.55 mutations/site/Myrs (HPD 95% confidence interval of 0.29-0.88 mutations/site/Myrs) after accounting for the persistence of these heteroplasmies and the sensitivity of current detection methods. In comparison, the rate of evolution (k) of the same HVR I region, determined using DNA sequences from 162 known age sub-fossil bones spanning a 37,000-year period, was 0.86 substitutions/site/Myrs (HPD 95% confidence interval of 0.53 and 1.17). Importantly, the latter rate is not statistically different from our estimate of the mutation rate. These results are in contrast to the view that molecular rates are time dependent.

A reevaluation of the Native American mtDNA genome diversity and its bearing on the models of early colonization of Beringia

The Americas were the last continents to be populated by humans, and their colonization represents a very interesting chapter in our species' evolution in which important issues are still contentious or largely unknown. One difficult topic concerns the details of the early peopling of Beringia, such as for how long it was colonized before people moved into the Americas and the demography of this occupation. A recent work using mitochondrial genome (mtDNA) data presented evidence for a so called "three-stage model" consisting of a very early expansion into Beringia followed by approximately 20,000 years of population stability before the final entry into the Americas. However, these results are in disagreement with other recent studies using similar data and methods. Here, we reanalyze their data to check the robustness of this model and test the ability of Native American mtDNA to discriminate details of the early colonization of Beringia. We apply the Bayesian Skyline Plot approach to recover the past demographic dynamic underpinning these events using different mtDNA data sets. Our results refute the specific details of the "three-stage model", since the early stage of expansion into Beringia followed by a long period of stasis could not be reproduced in any mtDNA data set cleaned from non-Native American haplotypes. Nevertheless, they are consistent with a moderate population bottleneck in Beringia associated with the Last Glacial Maximum followed by a strong population growth around 18,000 years ago as suggested by other recent studies. We suggest that this bottleneck erased the signals of ancient demographic history from recent Native American mtDNA pool, and conclude that the proposed early expansion and occupation of Beringia is an artifact caused by the misincorporation of non-Native American haplotypes.

The effect of ancient DNA damage on inferences of demographic histories

The field of ancient DNA (aDNA) is casting new light on many evolutionary questions. However, problems associated with the postmortem instability of DNA may complicate the interpretation of aDNA data. For example, in population genetic studies, the inclusion of damaged DNA may inflate estimates of diversity. In this paper, we examine the effect of DNA damage on population genetic estimates of ancestral population size. We simulate data using standard coalescent simulations that include postmortem damage and show that estimates of effective population sizes are inflated around, or right after, the sampling time of the ancestral DNA sequences. This bias leads to estimates of increasing, and then decreasing, population sizes, as observed in several recently published studies. We reanalyze a recently published data set of DNA sequences from the Bison (Bison bison/Bison priscus) and show that the signal for a change in effective population size in this data set vanishes once the effects of putative damage are removed. Our results suggest that population genetic analyses of aDNA sequences, which do not accurately account for damage, should be interpreted with great caution.

Mutation patterns of mtDNA: empirical inferences for the coding region

Background

Human mitochondrial DNA (mtDNA) has been extensively used in population and evolutionary genetics studies. Thus, a valid estimate of human mtDNA evolutionary rate is important in many research fields. The small number of estimations performed for the coding region of the molecule, showed important differences between phylogenetic and empirical approaches. We analyzed a portion of the coding region of mtDNA (tRNA^Leu, ND1 and tRNA^Ile genes), using individuals belonging to extended families from the Azores Islands (Portugal) with the main aim of providing empirical estimations of the mutation rate of the coding region of mtDNA under different assumptions, and hence to better understand the mtDNA evolutionary process.

Results

Heteroplasmy was detected in 6.5% (3/46) of the families analyzed. In all of the families the presence of mtDNA heteroplasmy resulted from three new point mutations, and no cases of insertions or deletions were identified. Major differences were found in the proportion and type of heteroplasmy found in the genes studied when compared to those obtained in a previous report for the D-loop. Our empirical estimation of mtDNA coding region mutation rate, calculated taking into account the sex of individuals carrying new mutations, the probability of intra-individual fixation of mutations present in heteroplasmy and, to the possible extent, the effect of selection, is similar to that obtained using phylogenetic approaches.

Conclusion

Based on our results, the discrepancy previously reported between the human mtDNA coding region mutation rates observed along evolutionary timescales and estimations obtained using family pedigrees can be resolved when correcting for the previously cited factors.

Pleistocene evolution of closely related sand martins Riparia riparia and R. diluta

Climatic fluctuations during the Quaternary resulted in a dynamic history of species' range shifts, fragmentations and expansions. Some of these events left traces in the genetic structures of plants and animals. Recent avian phylogeographic studies demonstrated that Holarctic birds responded idiosyncratically to Pleistocene climate fluctuations. We present phylogeographic analyses of the Holarctic collared sand martin (Riparia riparia) and the Asian pale sand martin (Riparia diluta), which were considered conspecific until recently. Mitochondrial and nuclear sequences confirm species status of the pale sand martin; the two species diverged sometime between late Pliocene and middle Pleistocene, but precise dates could not be provided without calibration of the substitution rate. Within the pale sand martin, we found two mitochondrial clades that are likely to have diverged in the Pleistocene, one from Central Siberia, and the other restricted to Mongolia. The two clades were sympatric with the collared sand martin in Buryatiya and Mongolia, respectively. The mitochondrial gene genealogy and phi(st) analysis of the collared sand martin haplotypes indicate recent, but not ongoing, gene exchange between North America and Eurasia, and restricted gene flow between western and eastern Siberia that likely resulted from historic fragmentation of the species' range during the last glacial maximum.

A Bayesian evaluation of human mitochondrial substitution rates

Accurate estimates of mitochondrial substitution rates are central to molecular studies of human evolution, but meaningful comparisons of published studies are problematic because of the wide range of methodologies and data sets employed. These differences are nowhere more pronounced than among rates estimated from phylogenies, genealogies, and pedigrees. By using a data set comprising mitochondrial genomes from 177 humans, we estimate substitution rates for various data partitions by using Bayesian phylogenetic analysis with a relaxed molecular clock. We compare the effect of multiple internal calibrations with the customary human-chimpanzee split. The analyses reveal wide variation among estimated substitution rates and divergence times made with different partitions and calibrations, with evidence of substitutional saturation, natural selection, and significant rate heterogeneity among lineages and among sites. Collectively, the results support dates for migration out of Africa and the common mitochondrial ancestor of humans that are considerably more recent than most previous estimates. Our results also demonstrate that human mitochondrial genomes exhibit a number of molecular evolutionary complexities that necessitate the use of sophisticated analytical models for genetic analyses.

Shrimps down under: evolutionary relationships of subterranean crustaceans from Western Australia (Decapoda: Atyidae: Stygiocaris)

BACKGROUND

We investigated the large and small scale evolutionary relationships of the endemic Western Australian subterranean shrimp genus Stygiocaris (Atyidae) using nuclear and mitochondrial genes. Stygiocaris is part of the unique cave biota of the coastal, anchialine, limestones of the Cape Range and Barrow Island, most of whose nearest evolutionary relations are found in coastal caves of the distant North Atlantic. The dominance of atyids in tropical waters and their food resources suggest they are pivotal in understanding these groundwater ecosystems.

METHODOLOGY/PRINCIPLE FINDINGS

Our nuclear and mitochondrial analyses all recovered the Mexican cave genus Typhlatya as the sister taxon of Stygiocaris, rather than any of the numerous surface and cave atyids from Australia or the Indo-Pacific region. The two described Stygiocaris species were recovered as monophyletic, and a third, cryptic, species was discovered at a single site, which has very different physiochemical properties from the sites hosting the two described species.

CONCLUSIONS/SIGNIFICANCE

Our findings suggest that Stygiocaris and Typhlatya may descend from a common ancestor that lived in the coastal marine habitat of the ancient Tethys Sea, and were subsequently separated by plate tectonic movements. This vicariant process is commonly thought to explain the many disjunct anchialine faunas, but has rarely been demonstrated using phylogenetic techniques. The Cape Range's geological dynamism, which is probably responsible for the speciation of the various Stygiocaris species, has also led to geographic population structure within species. In particular, Stygiocaris lancifera is split into northern and southern groups, which correspond to population splits within other sympatric subterranean taxa.

The effect of inappropriate calibration: three case studies in molecular ecology

Time-scales estimated from sequence data play an important role in molecular ecology. They can be used to draw correlations between evolutionary and palaeoclimatic events, to measure the tempo of speciation, and to study the demographic history of an endangered species. In all of these studies, it is paramount to have accurate estimates of time-scales and substitution rates. Molecular ecological studies typically focus on intraspecific data that have evolved on genealogical scales, but often these studies inappropriately employ deep fossil calibrations or canonical substitution rates (e.g., 1% per million years for birds and mammals) for calibrating estimates of divergence times. These approaches can yield misleading estimates of molecular time-scales, with significant impacts on subsequent evolutionary and ecological inferences. We illustrate this calibration problem using three case studies: avian speciation in the late Pleistocene, the demographic history of bowhead whales, and the Pleistocene biogeography of brown bears. For each data set, we compare the date estimates that are obtained using internal and external calibration points. In all three cases, the conclusions are significantly altered by the application of revised, internally-calibrated substitution rates. Collectively, the results emphasise the importance of judicious selection of calibrations for analyses of recent evolutionary events.

Under the volcano: phylogeography and evolution of the cave-dwelling Palmorchestia hypogaea (Amphipoda, Crustacea) at La Palma (Canary Islands)

BACKGROUND

The amphipod crustacean Palmorchestia hypogaea occurs only in La Palma (Canary Islands) and is one of the few terrestrial amphipods in the world that have adapted to a strictly troglobitic life in volcanic cave habitats. A surface-dwelling closely related species (Palmorchestia epigaea) lives in the humid laurel forest on the same island. Previous studies have suggested that an ancestral littoral Orchestia species colonized the humid forests of La Palma and that subsequent drought episodes in the Canaries reduced the distribution of P. epigaea favouring the colonization of lava tubes through an adaptive shift. This was followed by dispersal via the hypogean crevicular system.

RESULTS

P. hypogaea and P. epigaea did not form reciprocally monophyletic mitochondrial DNA clades. They showed geographically highly structured and genetically divergent populations with current gene flow limited to geographically close surface locations. Coalescence times using Bayesian estimations assuming a non-correlated relaxed clock with a normal prior distribution of the age of La Palma, together with the lack of association of habitat type with ancestral and recent haplotypes, suggest that their adaptation to cave life is relatively ancient.

CONCLUSION

The data gathered here provide evidence for multiple invasions of the volcanic cave systems that have acted as refuges. A re-evaluation of the taxonomic status of the extant species of Palmorchestia is needed, as the division of the two species by habitat and ecology is unnatural. The information obtained here, and that from previous studies on hypogean fauna, shows the importance of factors such as the uncoupling of morphological and genetic evolution, the role of climatic change and regressive evolution as key processes in leading to subterranean biodiversity.

Elevated substitution rates estimated from ancient DNA sequences

Ancient DNA sequences are able to offer valuable insights into molecular evolutionary processes, which are not directly accessible via modern DNA. They are particularly suitable for the estimation of substitution rates because their ages provide calibrating information in phylogenetic analyses, circumventing the difficult task of choosing independent calibration points. The substitution rates obtained from such datasets have typically been high, falling between the rates estimated from pedigrees and species phylogenies. Many of these estimates have been made using a Bayesian phylogenetic method that explicitly accommodates heterochronous data. Stimulated by recent criticism of this method, we present a comprehensive simulation study that validates its performance. For datasets of moderate size, it produces accurate estimates of rates, while appearing robust to assumptions about demographic history. We then analyse a large collection of 749 ancient and 727 modern DNA sequences from 19 species of animals, plants and bacteria. Our new estimates confirm that the substitution rates estimated from ancient DNA sequences are elevated above long-term phylogenetic levels.