Evidence of Site-Specific and Male-Biased Germline Mutation Rate in a Wild Songbird.

Germline mutations are the ultimate source of genetic variation and the raw material for organismal evolution. Despite their significance, the frequency and genomic locations of mutations, as well as potential sex bias, are yet to be widely investigated in most species. To address these gaps, we conducted whole-genome sequencing of 12 great reed warblers (Acrocephalus arundinaceus) in a pedigree spanning 3 generations to identify single-nucleotide de novo mutations (DNMs) and estimate the germline mutation rate. We detected 82 DNMs within the pedigree, primarily enriched at CpG sites but otherwise randomly located along the chromosomes. Furthermore, we observed a pronounced sex bias in DNM occurrence, with male warblers exhibiting three times more mutations than females. After correction for false negatives and adjusting for callable sites, we obtained a mutation rate of 7.16 × 10-9 mutations per site per generation (m/s/g) for the autosomes and 5.10 × 10-9 m/s/g for the Z chromosome. To demonstrate the utility of species-specific mutation rates, we applied our autosomal mutation rate in models reconstructing the demographic history of the great reed warbler. We uncovered signs of drastic population size reductions predating the last glacial period (LGP) and reduced gene flow between western and eastern populations during the LGP. In conclusion, our results provide one of the few direct estimates of the mutation rate in wild songbirds and evidence for male-driven mutations in accordance with theoretical expectations.

Solar heating may explain extreme diel flight altitude changes in migrating birds.

Great reed warblers, Acrocephalus arundinaceus,1 and great snipes, Gallinago media,2 exhibit a diel cycle in flight altitudes-flying much higher during the day than the night-when performing migratory flights covering both night and day. One hypothesis proposed to explain this behavior is that the birds face additional heating by solar radiation during daytime and hence must climb to very high, and thus also very cold, altitudes to avoid overheating during daytime flights.1,2 Yet, solar heat gain in birds has been shown to drastically decrease with wind speed,3,4 and the quantitative heating effect by solar radiation on a bird flying with an airspeed of 10 m/s or more is unknown. We analyzed temperature data from multisensor data loggers (MDLs)5,6 placed without direct exposure to solar radiation on great reed warblers (the logger covered by feathers on the back) and great snipes (the logger on the leg, covered from the sun by the tail). We found that logger temperatures were significantly higher (5.9°C-8.8°C in great reed warblers and 4.8°C-5.4°C in great snipes) during the day than during the night in birds flying at the same altitudes (and thus also the same expected ambient air temperatures). These results strongly indicate that the heat balance of the flying birds is indeed affected by solar radiation, which is in accordance with the hypothesis that solar radiation is a key factor causing the remarkable diel cycles in flight altitude observed in these two long-distance migrant bird species.1,2.

Heritability and Parental Effects in Telomere Length in a Color Polymorphic Long-Lived Bird.

AbstractRelative telomere length (RTL), an indicator of senescence, has been shown to be heritable but can also be affected by environmental factors, such as parental effects. Investigating heritability as well as parental effects and rearing environment can help us to understand the factors affecting offspring telomeres. Moreover, how phenotypic parental traits linked with fitness can impact offspring RTL is still unclear. A phenotypic marker closely associated with physiological traits and fitness is melanin-based color polymorphism, which in tawny owl (Strix aluco) is highly heritable and strongly associated with adult telomere shortening and survival. We studied narrow-sense heritability (h2) of RTL, as well as the impact of parental age and color morph and their interaction on offspring RTL. Offspring RTL at fledging was strongly positively correlated with both mother RTL and father RTL at breeding. Offspring RTL was also negatively associated with father age, suggesting that older fathers sired offspring with shorter telomeres. Parental color morph did not explain offspring RTL, and there were no interactive effects of parental morph and age, despite previously documented morph-specific senescence patterns. Our results suggest that RTL is highly heritable and affected by paternal age but not related to color polymorphism. This suggests that either morph-specific telomere shortening as an adult does not result in significantly shorter telomeres in their gametes, or that parents compensate morph-specific senescence via parental care. Morph-specific patterns of telomere dynamics in polymorphic species may thus emerge from different life history strategies adopted in adulthood.

Telomeres in ecology and evolution: A review and classification of hypotheses.

Research on telomeres in the fields of ecology and evolution has been rapidly expanding over the last two decades. This has resulted in the formulation of a multitude of, often name-given, hypotheses related to the associations between telomeres and life-history traits or fitness-facilitating processes (and the mechanisms underlying them). However, the differences (or similarities) between the various hypotheses, which can originate from different research fields, are often not obvious. Our aim here is therefore to give an overview of the hypotheses that are of interest in ecology and evolution and to provide two frameworks that help discriminate among them. We group the hypotheses (i) based on their association with different research questions, and (ii) using a hierarchical approach that builds on the assumptions they make, such as about causality of telomere length/shortening and/or the proposed functional consequences of telomere shortening on organism performance. Both our frameworks show that there exist parallel lines of thoughts in different research fields. Moreover, they also clearly illustrate that there are in many cases competing hypotheses within clusters, and that some of these even have contradictory assumptions and/or predictions. We also touch upon two topics in telomere research that would benefit from further conceptualization. This review should help researchers, both those familiar with and those new to the subject, to identify future avenues of research.

Avian Neo-Sex Chromosomes Reveal Dynamics of Recombination Suppression and W Degeneration.

How the avian sex chromosomes first evolved from autosomes remains elusive as 100 million years (My) of divergence and degeneration obscure their evolutionary history. The Sylvioidea group of songbirds is interesting for understanding avian sex chromosome evolution because a chromosome fusion event ∼24 Ma formed "neo-sex chromosomes" consisting of an added (new) and an ancestral (old) part. Here, we report the complete female genome (ZW) of one Sylvioidea species, the great reed warbler (Acrocephalus arundinaceus). Our long-read assembly shows that the added region has been translocated to both Z and W, and whereas the added-Z has retained its gene order the added-W part has been heavily rearranged. Phylogenetic analyses show that recombination between the homologous added-Z and -W regions continued after the fusion event, and that recombination suppression across this region took several million years to be completed. Moreover, recombination suppression was initiated across multiple positions over the added-Z, which is not consistent with a simple linear progression starting from the fusion point. As expected following recombination suppression, the added-W show signs of degeneration including repeat accumulation and gene loss. Finally, we present evidence for nonrandom maintenance of slowly evolving and dosage-sensitive genes on both ancestral- and added-W, a process causing correlated evolution among orthologous genes across broad taxonomic groups, regardless of sex linkage.

Extreme altitudes during diurnal flights in a nocturnal songbird migrant.

Billions of nocturnally migrating songbirds fly across oceans and deserts on their annual journeys. Using multisensor data loggers, we show that great reed warblers (Acrocephalus arundinaceus) regularly prolong their otherwise strictly nocturnal flights into daytime when crossing the Mediterranean Sea and the Sahara Desert. Unexpectedly, when prolonging their flights, they climbed steeply at dawn, from a mean of 2394 meters above sea level to reach extreme cruising altitudes (mean 5367 and maximum 6267 meters above sea level) during daytime flights. This previously unknown behavior of using exceedingly high flight altitudes when migrating during daytime could be caused by diel variation in ambient temperature, winds, predation, vision range, and solar radiation. Our finding of this notable behavior provides new perspectives on constraints in bird flight and might help to explain the evolution of nocturnal migration.

A genome-wide linkage map for the house sparrow (Passer domesticus) provides insights into the evolutionary history of the avian genome.

The house sparrow is an important model species for studying physiological, ecological and evolutionary processes in wild populations. Here, we present a medium density, genome wide linkage map for house sparrow (Passer domesticus) that has aided the assembly of the house sparrow reference genome, and that will provide an important resource for ongoing mapping of genes controlling important traits in the ecology and evolution of this species. Using a custom house sparrow 10 K iSelect Illumina SNP chip we have assigned 6,498 SNPs to 29 autosomal linkage groups, based on a mean of 430 informative meioses per SNP. The map was constructed by combining the information from linkage with that of the physical position of SNPs within scaffold sequences in an iterative process. Averaged between the sexes; the linkage map had a total length of 2,004 cM, with a longer map for females (2,240 cM) than males (1,801 cM). Additionally, recombination rates also varied along the chromosomes. Comparison of the linkage map to the reference genomes of zebra finch, collared flycatcher and chicken, showed a chromosome fusion of the two avian chromosomes 8 and 4A in house sparrow. Lastly, information from the linkage map was utilized to conduct analysis of linkage disequilibrium (LD) in eight populations with different effective population sizes (Ne ) in order to quantify the background level LD. Together, these results aid the design of future association studies, facilitate the development of new genomic tools and support the body of research that describes the evolution of the avian genome.

Adaptive responses of animals to climate change are most likely insufficient.

Biological responses to climate change have been widely documented across taxa and regions, but it remains unclear whether species are maintaining a good match between phenotype and environment, i.e. whether observed trait changes are adaptive. Here we reviewed 10,090 abstracts and extracted data from 71 studies reported in 58 relevant publications, to assess quantitatively whether phenotypic trait changes associated with climate change are adaptive in animals. A meta-analysis focussing on birds, the taxon best represented in our dataset, suggests that global warming has not systematically affected morphological traits, but has advanced phenological traits. We demonstrate that these advances are adaptive for some species, but imperfect as evidenced by the observed consistent selection for earlier timing. Application of a theoretical model indicates that the evolutionary load imposed by incomplete adaptive responses to ongoing climate change may already be threatening the persistence of species.

Evidence for sexual conflict over major histocompatibility complex diversity in a wild songbird.

Sex differences in parasite load and immune responses are found across a wide range of animals, with females generally having lower parasite loads and stronger immune responses than males. Intrigued by these general patterns, we investigated if there was any sign of sex-specific selection on an essential component of adaptive immunity that is known to affect fitness, the major histocompatibility complex class I (MHC-I) genes, in a 20-year study of great reed warblers. Our analyses on fitness related to MHC-I diversity showed a highly significant interaction between MHC-I diversity and sex, where males with higher, and females with lower, MHC-I diversity were more successful in recruiting offspring. Importantly, mean MHC-I diversity did not differ between males and females, and consequently neither sex reached its MHC-I fitness optimum. Thus, there is an unresolved genetic sexual conflict over MHC-I diversity in great reed warblers. Selection from pathogens is known to maintain MHC diversity, but previous theory ignores that the immune environments are considerably different in males and females. Our results suggest that sexually antagonistic selection is an important, previously neglected, force in the evolution of vertebrate adaptive immunity, and have implications for evolutionary understanding of costs of immune responses and autoimmune diseases.

Sex differences in life history, behavior, and physiology along a slow-fast continuum: a meta-analysis.

The pace-of-life syndrome (POLS) hypothesis predicts that behavior and physiology covary with life history. Evidence for such covariation is contradictory, possibly because systematic sources of variation (e.g. sex) have been neglected. Sexes often experience different selection pressures leading to sex-specific allocation between reproduction and self-maintenance, facilitating divergence in life-history. Sex-specific differences in means and possibly variances may therefore play a key role in the POLS framework. We investigate whether sexes differ in means and variances along the fast-slow pace-of-life continuum for life history and physiological and behavioral traits. In addition, we test whether social and environmental characteristics such as breeding strategy, mating system, and study environment explain heterogeneity between the sexes. Using meta-analytic methods, we found that populations with a polygynous mating system or for studies conducted on wild populations, males had a faster pace-of-life for developmental life-history traits (e.g., growth rate), behavior, and physiology. In contrast, adult life-history traits (e.g., lifespan) were shifted towards faster pace-of-life in females, deviating from the other trait categories. Phenotypic variances were similar between the sexes across trait categories and were not affected by mating system or study environment. Breeding strategy did not influence sex differences in variances or means. We discuss our results in the light of sex-specific selection that might drive sex-specific differences in pace-of-life and ultimately POLS.

Contrasting results from GWAS and QTL mapping on wing length in great reed warblers.

A major goal in evolutionary biology is to understand the genetic basis of adaptive traits. In migratory birds, wing morphology is such a trait. Our previous work on the great reed warbler (Acrocephalus arundinaceus) shows that wing length is highly heritable and under sexually antagonistic selection. Moreover, a quantitative trait locus (QTL) mapping analysis detected a pronounced QTL for wing length on chromosome 2, suggesting that wing morphology is partly controlled by genes with large effects. Here, we re-evaluate the genetic basis of wing length in great reed warblers using a genomewide association study (GWAS) approach based on restriction site-associated DNA sequencing (RADseq) data. We use GWAS models that account for relatedness between individuals and include covariates (sex, age and tarsus length). The resulting association landscape was flat with no peaks on chromosome 2 or elsewhere, which is in line with expectations for polygenic traits. Analysis of the distribution of p-values did not reveal biases, and the inflation factor was low. Effect sizes were however not uniformly distributed on some chromosomes, and the Z chromosome had weaker associations than autosomes. The level of linkage disequilibrium (LD) in the population decayed to background levels within c. 1 kbp. There could be several reasons to why our QTL study and GWAS gave contrasting results including differences in how associations are modelled (cosegregation in pedigree vs. LD associations), how covariates are accounted for in the models, type of marker used (multi- vs. biallelic), difference in power or a combination of these. Our study highlights that the genetic architecture even of highly heritable traits is difficult to characterize in wild populations.

Evolution of sex-specific pace-of-life syndromes: genetic architecture and physiological mechanisms.

Sex differences in life history, physiology, and behavior are nearly ubiquitous across taxa, owing to sex-specific selection that arises from different reproductive strategies of the sexes. The pace-of-life syndrome (POLS) hypothesis predicts that most variation in such traits among individuals, populations, and species falls along a slow-fast pace-of-life continuum. As a result of their different reproductive roles and environment, the sexes also commonly differ in pace-of-life, with important consequences for the evolution of POLS. Here, we outline mechanisms for how males and females can evolve differences in POLS traits and in how such traits can covary differently despite constraints resulting from a shared genome. We review the current knowledge of the genetic basis of POLS traits and suggest candidate genes and pathways for future studies. Pleiotropic effects may govern many of the genetic correlations, but little is still known about the mechanisms involved in trade-offs between current and future reproduction and their integration with behavioral variation. We highlight the importance of metabolic and hormonal pathways in mediating sex differences in POLS traits; however, there is still a shortage of studies that test for sex specificity in molecular effects and their evolutionary causes. Considering whether and how sexual dimorphism evolves in POLS traits provides a more holistic framework to understand how behavioral variation is integrated with life histories and physiology, and we call for studies that focus on examining the sex-specific genetic architecture of this integration.

Did natural selection make the Dutch taller? A cautionary note on the importance of quantification in understanding evolution.

One of the main achievements of the modern synthesis is a rigorous mathematical theory for evolution by natural selection. Combining this theory with statistical models makes it possible to estimate the relevant parameters so as to quantify selection and evolution in nature. Although quantification is a sign of a mature science, statistical models are unfortunately often interpreted independently of the motivating mathematical theory. Without a link to theory, numerical results do not represent proper quantifications, because they lack the connections that designate their biological meaning. Here, we want to raise awareness and exemplify this problem by examining a recent study on natural selection in a contemporary human population. Stulp et al. (2015) concluded that natural selection may partly explain the increasing stature of the Dutch population. This conclusion was based on a qualitative assessment of the presence of selection on height. Here, we provide a quantitative interpretation of these results using standard evolutionary theory to show that natural selection has had a minuscule effect.

Maternal and genetic factors determine early life telomere length.

In a broad range of species--including humans--it has been demonstrated that telomere length declines throughout life and that it may be involved in cell and organismal senescence. This potential link to ageing and thus to fitness has triggered recent interest in understanding how variation in telomere length is inherited and maintained. However, previous studies suffer from two main drawbacks that limit the possibility of understanding the relative importance of genetic, parental and environmental influences on telomere length variation. These studies have been based on (i) telomere lengths measured at different time points in different individuals, despite the fact that telomere length changes over life, and (ii) parent-offspring regression techniques, which do not enable differentiation between genetic and parental components of inheritance. To overcome these drawbacks, in our study of a songbird, the great reed warbler, we have analysed telomere length measured early in life in both parents and offspring and applied statistical models (so-called 'animal models') that are based on long-term pedigree data. Our results showed a significant heritability of telomere length on the maternal but not on the paternal side, and that the mother's age was positively correlated with their offspring's telomere length. Furthermore, the pedigree-based analyses revealed a significant heritability and an equally large maternal effect. Our study demonstrates strong maternal influence on telomere length and future studies now need to elucidate possible underlying factors, including which types of maternal effects are involved.

Assessing multivariate constraints to evolution across ten long-term avian studies.

BACKGROUND: In a rapidly changing world, it is of fundamental importance to understand processes constraining or facilitating adaptation through microevolution. As different traits of an organism covary, genetic correlations are expected to affect evolutionary trajectories. However, only limited empirical data are available. METHODOLOGY/PRINCIPAL FINDINGS: We investigate the extent to which multivariate constraints affect the rate of adaptation, focusing on four morphological traits often shown to harbour large amounts of genetic variance and considered to be subject to limited evolutionary constraints. Our data set includes unique long-term data for seven bird species and a total of 10 populations. We estimate population-specific matrices of genetic correlations and multivariate selection coefficients to predict evolutionary responses to selection. Using Bayesian methods that facilitate the propagation of errors in estimates, we compare (1) the rate of adaptation based on predicted response to selection when including genetic correlations with predictions from models where these genetic correlations were set to zero and (2) the multivariate evolvability in the direction of current selection to the average evolvability in random directions of the phenotypic space. We show that genetic correlations on average decrease the predicted rate of adaptation by 28%. Multivariate evolvability in the direction of current selection was systematically lower than average evolvability in random directions of space. These significant reductions in the rate of adaptation and reduced evolvability were due to a general nonalignment of selection and genetic variance, notably orthogonality of directional selection with the size axis along which most (60%) of the genetic variance is found. CONCLUSIONS: These results suggest that genetic correlations can impose significant constraints on the evolution of avian morphology in wild populations. This could have important impacts on evolutionary dynamics and hence population persistence in the face of rapid environmental change.

Intralocus sexual conflict over wing length in a wild migratory bird.

Intralocus sexual conflict (ISC) occurs when males and females have different adaptive peaks but are constrained from evolving sexual dimorphism because of shared genes. Implications of this conflict on evolutionary dynamics in wild populations have not been investigated in detail. In comprehensive analyses of selection, heritability, and genetic correlations, we found evidence for an ISC over wing length, a key trait for flight performance and migration, in a long-term study of wild great reed warblers (Acrocephalus arundinaceus). We found moderate sexual dimorphism, high heritability, moderate sexually antagonistic selection, and strong positive cross-sex genetic correlation in wing length, together supporting the presence of ISC. A negative genetic correlation between male wing length and female fitness indicated that females inheriting alleles for longer wings from their male relatives also inherited lower fitness. Moreover, cross-sex genetic correlations imposed constraint on the predicted microevolutionary trajectory of wing length (based on selection gradients), especially in females where the predicted response was reversed. The degree of sexual dimorphism in wing length did not change over time, suggesting no sign of conflict resolution. Our study provides novel insight into how an ISC over a fitness trait can affect microevolution in a wild population under natural selection.

Annual cycle and migration strategies of a trans-Saharan migratory songbird: a geolocator study in the great reed warbler.

Recent technological advancements now allow us to obtain geographical position data for a wide range of animal movements. Here we used light-level geolocators to study the annual migration cycle in great reed warblers (Acrocephalus arundinaceus), a passerine bird breeding in Eurasia and wintering in sub-Saharan Africa. We were specifically interested in seasonal strategies in routes and schedules of migration. We found that the great reed warblers (all males, no females were included) migrated from the Swedish breeding site in early August. After spending up to three weeks at scattered stopover sites in central to south-eastern Europe, they resumed migration and crossed the Mediterranean Sea and Sahara Desert without lengthy stopovers. They then spread out over a large overwintering area and each bird utilised two (or even three) main wintering sites that were spatially separated by a distinct mid-winter movement. Spring migration initiation date differed widely between individuals (1-27 April). Several males took a more westerly route over the Sahara in spring than in autumn, and in general there were fewer long-distance travels and more frequent shorter stopovers, including one in northern Africa, in spring. The shorter stopovers made spring migration on average faster than autumn migration. There was a strong correlation between the spring departure dates from wintering sites and the arrival dates at the breeding ground. All males had a high migration speed in spring despite large variation in departure dates, indicating a time-minimization strategy to achieve an early arrival at the breeding site; the latter being decisive for high reproductive success in great reed warblers. Our results have important implications for the understanding of long-distance migrants' ability to predict conditions at distant breeding sites and adapt to rapid environmental change.

Hybridization but no evidence for backcrossing and introgression in a sympatric population of great reed warblers and clamorous reed warblers.

Hybridization is observed frequently in birds, but often it is not known whether the hybrids are fertile and if backcrossing occurs. The breeding ranges of the great reed warbler (Acrocephalus arundinaceus) and the clamorous reed warbler (A. stentoreus) overlap in southern Kazakhstan and a previous study has documented hybridization in a sympatric population. In the present study, we first present a large set of novel microsatellite loci isolated and characterised in great reed warblers. Secondly, we evaluate whether hybridization in the sympatric breeding population has been followed by backcrossing and introgression.We isolated 181 unique microsatellite loci in great reed warblers. Of 41 loci evaluated, 40 amplified and 30 were polymorphic. Bayesian clustering analyses based on genotype data from 23 autosomal loci recognised two well-defined genetic clusters corresponding to the two species. Individuals clustered to a very high extent to either of these clusters (admixture proportions ≥ 0.984) with the exception of four previously suggested arundinaceus-stentoreus hybrid birds that showed mixed ancestry (admixture proportions 0.495-0.619). Analyses of simulated hybrids and backcrossed individuals showed that the sampled birds do not correspond to first-fourth-generation backcrosses, and that fifth or higher generation backcrosses to a high extent resemble 'pure' birds at this set of markers.We conclude that these novel microsatellite loci provide a useful molecular resource for Acrocephalus warblers. The time to reach reproductive isolation is believed to be very long in birds, approximately 5 Myrs, and with an estimated divergence time of 2 Myrs between these warblers, some backcrossing and introgression could have been expected. However, there was no evidence for backcrossing and introgression suggesting that hybrids are either infertile or their progeny inviable. Very low levels of introgression cannot be excluded, which still may be an important factor as a source of new genetic variation.

A strong quantitative trait locus for wing length on chromosome 2 in a wild population of great reed warblers.

Wing length is a key character for essential behaviours related to bird flight such as migration and foraging. In the present study, we initiate the search for the genes underlying wing length in birds by studying a long-distance migrant, the great reed warbler (Acrocephalus arundinaceus). In this species wing length is an evolutionary interesting trait with pronounced latitudinal gradient and sex-specific selection regimes in local populations. We performed a quantitative trait locus (QTL) scan for wing length in great reed warblers using phenotypic, genotypic, pedigree and linkage map data from our long-term study population in Sweden. We applied the linkage analysis mapping method implemented in GridQTL (a new web-based software) and detected a genome-wide significant QTL for wing length on chromosome 2, to our knowledge, the first detected QTL in wild birds. The QTL extended over 25 cM and accounted for a substantial part (37%) of the phenotypic variance of the trait. A genome scan for tarsus length (a body-size-related trait) did not show any signal, implying that the wing-length QTL on chromosome 2 was not associated with body size. Our results provide a first important step into understanding the genetic architecture of avian wing length, and give opportunities to study the evolutionary dynamics of wing length at the locus level.

Postglacial colonisation patterns and the role of isolation and expansion in driving diversification in a passerine bird.

Pleistocene glacial cycles play a major role in diversification and speciation, although the relative importance of isolation and expansion in driving diversification remains debated. We analysed mitochondrial DNA sequence data from 15 great reed warbler (Acrocephalus arundinaceus) populations distributed over the vast Eurasian breeding range of the species, and revealed unexpected postglacial expansion patterns from two glacial refugia. There were 58 different haplotypes forming two major clades, A and B. Clade A dominated in Western Europe with declining frequencies towards Eastern Europe and the Middle East, but showed a surprising increase in frequency in Western and Central Asia. Clade B dominated in the Middle East, with declining frequencies towards north in Central and Eastern Europe and was absent from Western Europe and Central Asia. A parsimonious explanation for these patterns is independent postglacial expansions from two isolated refugia, and mismatch distribution analyses confirmed this suggestion. Gene flow analyses showed that clade A colonised both Europe and Asia from a refugium in Europe, and that clade B expanded much later and colonised parts of Europe from a refugium in the Middle East. Great reed warblers in the eastern parts of the range have slightly paler plumage than western birds (sometimes treated as separate subspecies; A. a. zarudnyi and A. a. arundinaceus, respectively) and our results suggest that the plumage diversification took place during the easterly expansion of clade A. This supports the postglacial expansion hypothesis proposing that postglacial expansions drive diversification in comparatively short time periods. However, there is no indication of any (strong) reproductive isolation between clades and our data show that the refugia populations became separated during the last glaciation. This is in line with the Pleistocene speciation hypothesis invoking that much longer periods of time in isolation are needed for speciation to occur.