Phylogeography of the world's tallest angiosperm, Eucalyptus regnans: evidence for multiple isolated Quaternary refugia

Aim There is a need for more Southern Hemisphere phylogeography studies, particularly in Australia, where, unlike much of Europe and North America, ice sheet cover was not extensive during the Last Glacial Maximum (LGM). This study examines the phylogeography of the south‐east Australian montane tree species Eucalyptus regnans. The work aimed to identify any major evolutionary divergences or disjunctions across the species’ range and to examine genetic signatures of past range contraction and expansion events. Location South‐eastern mainland Australia and the large island of Tasmania. Methods We determined the chloroplast DNA haplotypes of 410 E. regnans individuals (41 locations) based on five chloroplast microsatellites. Genetic structure was examined using analysis of molecular variance (AMOVA), and a statistical parsimony tree was constructed showing the number of nucleotide differences between haplotypes. Geographic structure in population genetic diversity was examined with the calculation of diversity parameters for the mainland and Tasmania, and for 10 regions. Regional analysis was conducted to test hypotheses that some areas within the species’ current distribution were refugia during the LGM and that other areas have been recolonized by E. regnans since the LGM. Results Among the 410 E. regnans individuals analysed, 31 haplotypes were identified. The statistical parsimony tree shows that haplotypes divided into two distinct groups corresponding to mainland Australia and Tasmania. The distribution of haplotypes across the range of E. regnans shows strong geographic patterns, with many populations and even certain regions in which a particular haplotype is fixed. Many locations had unique haplotypes, particularly those in East Gippsland in south‐eastern mainland Australia, north‐eastern Tasmania and south‐eastern Tasmania. Higher haplotype diversity was found in putative refugia, and lower haplotype diversity in areas likely to have been recolonized since the LGM. Main conclusions The data are consistent with the long‐term persistence of E. regnans in many regions and the recent recolonization of other regions, such as the Central Highlands of south‐eastern mainland Australia. This suggests that, in spite of the narrow ecological tolerances of the species and the harsh environmental conditions during the LGM, E. regnans was able to persist locally or contracted to many near‐coastal refugia, maintaining a diverse genetic structure.     

Multiple glacial refugia for cool‐temperate deciduous trees in northern East Asia: the Mongolian oak as a case study

In East Asia, temperate forests are predicted to have retracted southward to c. 30° N during the last glacial maximum (LGM) based on fossil pollen data, whereas phylogeographic studies have often suggested glacial in situ survival of cool‐temperate deciduous trees in their modern northern ranges. Here we report a study of the genetic diversity and structure of 29 natural Mongolian oak (Quercus mongolica) populations using 19 nuclear simple sequence repeat (nSSR) loci and four chloroplast DNA fragments. Bayesian clustering analysis with nSSRs revealed five groups, which were inferred by approximate Bayesian computation (ABC) to have diverged in multiple refugia through multiple glacial–interglacial cycles. Analysis of chloroplast DNA variation revealed four lineages that were largely but incompletely geographically disjunct. Ecological niche modelling (ENMs) indicated a southward range shift of the oak's distribution at the LGM, although high suitability scores were also evident in the Changbai Mts. (Northeast China), the Korean Peninsula, areas surrounding the Bohai Sea, and along the coast of the Russian Far East. In addition, endemic chloroplast DNA haplotypes and nuclear lineages occurred in high‐latitude northern areas where the ENM predicted no suitable habitat. The combined evidence from nuclear and chloroplast DNA, and the results of the ENM clearly demonstrate that multiple northern refugia, including cryptic ones, were maintained across the current distributional range of the Mongolian oak during the LGM or earlier glacial periods. Though spatially limited, postglacial expansions from these refugia have led to a pattern of decreased genetic diversity with increasing latitude.     

Phylogeography and ecological niche modelling of the New Zealand stick insect Clitarchus hookeri (White) support survival in multiple coastal refugia

Aim Increasing our understanding of the effects of the Last Glacial Maximum (LGM) and determining the location of refugia requires studies on widely distributed species with dense sampling of populations. We have reconstructed the biogeographic history of Clitarchus hookeri (White), a widespread species of New Zealand stick insect that exhibits geographic parthenogenesis, using phylogeographic analysis and ecological niche modelling. Location New Zealand. Methods We used DNA sequence data from the mitochondrial cytochrome c oxidase subunit I gene to reconstruct phylogenetic relationships among haplotypes from C. hookeri and two undescribed Clitarchus species. We also used distribution data from our own field surveys and museum records to reconstruct the geographic distribution of C. hookeri during the present and the LGM, using ecological niche modelling. Results The ecological niche models showed that the geographic distribution of C. hookeri has expanded dramatically since the LGM. Our model predicted large areas of suitable LGM habitat in upper North Island, and small patches along the east coast of South Island. The phylogeographic analysis shows that populations in the northern half of North Island contain much higher levels of genetic variation than those from southern North Island and South Island, and is congruent with the ecological niche model. The distribution of bisexual populations is also non-random, with males completely absent from South Island and very rare in southern North Island. Main conclusions During the LGM C. hookeri was most likely restricted to several refugia in upper North Island and one or more smaller refugia along the east coast of South Island. The unisexual populations predominate in post-glacial landscapes and are clearly favoured in the recolonization of such areas. Our study exemplifies the utility of integrating ecological niche modelling and phylogeographic analysis.     

Low but structured chloroplast diversity in Atherosperma moschatum (Atherospermataceae) suggests bottlenecks in response to the Pleistocene glacials

Background and Aims

The cool temperate rainforests of Australia were much reduced in range during the cold and dry glacial periods, although genetic evidence indicates that two key rainforest species, Nothofagus cunninghamii and Tasmannia lanceolata, survived within multiple locations and underwent only local range expansions at the end of the Last Glacial. To better understand the glacial response of a co-occurring but wind-dispersed and less cold-tolerant rainforest tree species, Atherosperma moschatum, a chloroplast phylogeographic study was undertaken.

Methods

A total of 3294 bp of chloroplast DNA sequence was obtained for 155 samples collected from across the species'' range.

Key Results

The distribution of six haplotypes observed in A. moschatum was geographically structured with an inferred ancestral haplotype restricted to Tasmania, while three non-overlapping and endemic haplotypes were found on the mainland of south-eastern Australia. Last glacial refugia for A. moschatum are likely to have occurred in at least one location in western Tasmania and in Victoria and within at least two locations in the Great Dividing Range of New South Wales. Nucleotide diversity of A. moschatum was lower (π = 0·00021) than either N. cunninghamii (0·00101) or T. lanceolata (0·00073), and was amongst the lowest recorded for any tree species.

Conclusions

This study provides evidence for past bottlenecks having impacted the chloroplast diversity of A. moschatum as a result of the species narrower climatic niche during glacials. This hypothesis is supported by the star-like haplotype network and similar estimated rates of chloroplast DNA substitution for A. moschatum and the two more cold tolerant and co-occurring species that have higher chloroplast diversity, N. cunninghamii and T. lanceolata.     

Shared phylogeographic patterns and widespread chloroplast haplotype sharing in Eucalyptus species with different ecological tolerances

We examined the phylogeography of three south-east Australian trees (Eucalyptus delegatensis, Eucalyptus obliqua, and Eucalyptus regnans) with different tolerances, in terms of cold, drought, fire and soil to explore whether species with different ecologies share major phylogeographic patterns. A second aim of this study was to examine geographic patterns of chloroplast DNA (cpDNA) haplotype sharing among the three study species. Trees of E. delegatensis (n?=?120), E. obliqua (n?=?265) and E. regnans (n?=?270) were genotyped with five cpDNA microsatellite markers. The species shared major phylogeographic disjunctions, and common patterns at proposed glacial refugia (generally high haplotype diversity) and areas thought to have been treeless during the Last Glacial Maximum (LGM) (low diversity). Inter-specific sharing of haplotypes was extensive, and fixation of shared, regional haplotypes was more frequent in areas postulated as having been treeless at the LGM. Despite ecological differences, chloroplast microsatellite data suggest the three species have responded to past climatic changes in a similar way, by persisting in multiple, generally common refugia. We suggest that the natural ability of eucalypt species to hybridise with others with quite different or broader ecological tolerances may provide an “insurance policy” for response to rapidly changing abiotic conditions.     

Phylogeographic evidence for two mesic refugia in a biodiversity hotspot

Phylogeographic studies of flora in species-rich south-western Australia point to complex evolutionary histories, reflecting patterns of persistence and resilience to climatic changes during the Pleistocene. We asked whether coastal areas of the mid-west and south, as well as granite outcrops and inland ranges, have acted as major refugia within this region during Pleistocene climatic fluctuations by analysing phylogeographic patterns in the shrub Calothamnus quadrifidus R.Br. (Myrtaceae). We determined variation in chloroplast DNA data for 41 populations across the geographic range. Relationships and major clades were resolved using parsimony and Bayesian analyses. We tested for demographic and spatial expansion of the major clades and estimated clade divergence dates using an uncorrelated, lognormal relaxed clock based on two conservative chloroplast mutation rates. Two distinct phylogeographic clades were identified showing divergence during the Pleistocene, consistent with other phylogeographic studies of south-west Australian flora, emphasising the impact of climatic oscillations in driving divergence in this landscape. The southern clade was more diverse, having higher haplotype diversity and greater genetic structure, while the northern clade showed evidence of fluctuation in population size. Regions of high haplotype diversity with adjacent areas of low diversity observed in each clade indicated the locations of two coastal refugia: one on the south coast and another along the mid-west coast. This is the first evidence for major Pleistocene refugia using chloroplast genetic data in a common, widespread species from this region.     

Identification of glacial refugia in south-eastern North America by phylogeographical analyses of a forest understorey plant, Trillium cuneatum

Aim We examine several hypotheses emerging from biogeographical and fossil records regarding glacial refugia of a southern thermophilic plant species. Specifically, we investigated the glacial history and post‐glacial colonization of a forest understorey species, Trillium cuneatum. We focused on the following questions: (1) Did T. cuneatum survive the Last Glacial Maximum (LGM) in multiple refugia, and (if so) where were they located, and is the modern genetic structure congruent with the fossil record‐based reconstruction of refugia for mesic deciduous forests? (2) What are the post‐glacial colonization patterns in the present geographical range? Location South‐eastern North America. Methods We sampled 45 populations of T. cuneatum throughout its current range. We conducted phylogeographical analyses based on maternally inherited chloroplast DNA (cpDNA haplotypes) and used TCS software to reconstruct intraspecific phylogeny. Results We detected six cpDNA haplotypes, geographically highly structured into non‐overlapping areas. With one exception, none of the populations had mixed haplotype composition. TCS analysis resulted in two intraspecific cpDNA lineages, with one clade subdivided further by shallower diversification. Main conclusions Our investigation revealed that T. cuneatum survived the LGM in multiple refugia, belonging to two (western, eastern) genealogical lineages geographically structured across south‐eastern North America. The western clade is confined to the south‐western corner of T. cuneatum’s modern range along the Lower Mississippi Valley, where fossil records document a major refugium of mesic deciduous forest. For the eastern clade, modern patterns of cpDNA haplotype distribution suggest cryptic vicariance, in the form of forest contractions and subsequent expansions associated with Pleistocene glacial cycles, rather than simple southern survival and subsequent northward colonization. The north–south partitioning of cpDNA haplotypes was unexpected, suggesting that populations of this rather southern thermophilic species may have survived in more northern locations than initially expected based on LGM climate reconstruction, and that the Appalachian Mountains functioned as a barrier to the dispersal of propagules originating in more southern refugia. Furthermore, our results reveal south‐west to north‐east directionality in historical migration through the Valley and Ridge region of north‐west Georgia.     

Northern glacial refugia for the pygmy shrew Sorex minutus in Europe revealed by phylogeographic analyses and species distribution modelling

The southern European peninsulas (Iberian, Italian and Balkan) are traditionally recognized as glacial refugia from where many species colonized central and northern Europe after the Last Glacial Maximum (LGM). However, evidence that some species had more northerly refugia is accumulating from phylogeographic, palaeontological and palynological studies, and more recently from species distribution modelling (SDM), but further studies are needed to test the idea of northern refugia in Europe. Here, we take a rarely implemented multidisciplinary approach to assess if the pygmy shrew Sorex minutus, a widespread Eurasian mammal species, had northern refugia during the LGM, and if these influenced its postglacial geographic distribution. First, we evaluated the phylogeographic and population expansion patterns using mtDNA sequence data from 123 pygmy shrews. Then, we used SDM to predict present and past (LGM) potential distributions using two different training data sets, two different algorithms (Maxent and GARP) and climate reconstructions for the LGM with two different general circulation models. An LGM distribution in the southern peninsulas was predicted by the SDM approaches, in line with the occurrence of lineages of S. minutus in these areas. The phylogeographic analyses also indicated a widespread and strictly northern‐central European lineage, not derived from southern peninsulas, and with a postglacial population expansion signature. This was consistent with the SDM predictions of suitable LGM conditions for S. minutus occurring across central and eastern Europe, from unglaciated parts of the British Isles to much of the eastern European Plain. Hence, S. minutus likely persisted in parts of central and eastern Europe during the LGM, from where it colonized other northern areas during the late‐glacial and postglacial periods. Our results provide new insights into the glacial and postglacial colonization history of the European mammal fauna, notably supporting glacial refugia further north than traditionally recognized.     

Extensive intraspecific genetic diversity of a freshwater crayfish in a biodiversity hotspot

1. The freshwater crayfish Cherax dispar (Decapoda: Parastacidae) inhabits coastal regions and islands of South East Queensland, Australia. We hypothesised that populations of C. dispar on different islands would be more genetically divergent from each other than populations from different drainages within the same island or on the mainland. 2. Phylogenetic and phylogeographic analyses were conducted on two mitochondrial genes (cytochrome oxidase subunit I & 16S ribosomal DNA) and one nuclear gene (Internal Transcribed Spacer region 2). Phylogeographic patterns were compared with those for other freshwater organisms in the area. 3. Deep genetic divergences were found within C. dispar, including four highly divergent (up to 20%) clades. The geographic distribution of each of the clades revealed strong latitudinal structuring along the coast rather than structuring among the islands. The high genetic divergence observed among the C. dispar clades was estimated to have pre‐dated island formation and may represent ancient river drainage patterns. 4. A restricted distribution was observed for the most divergent clade, which was discovered only on two of the sand islands (North Stradbroke Island and Moreton Island). Furthermore, strong phylogeographic structuring was observed within this clade on North Stradbroke Island, where no haplotypes were shared between samples from opposite sides of the island. This low connectivity within the island supports the idea that C. dispar rarely disperse terrestrially (i.e. across watersheds).     

Phylogeographic analysis of North American populations of the parasitic herbaceous plant Monotropa hypopitys L. reveals a complex history of range expansion from multiple late glacial refugia

Aim We carried out a phylogeographic study across the range of the herbaceous plant species Monotropa hypopitys L. in North America to determine whether its current disjunct distribution is due to recolonization from separate eastern and western refugia after the Last Glacial Maximum (LGM). Location North America: Pacific Northwest and north‐eastern USA/south‐eastern Canada. Methods Palaeodistribution modelling was carried out to determine suitable climatic regions for M. hypopitys at the LGM. We analysed between 155 and 176 individuals from 39 locations spanning the species’ entire range in North America. Sequence data were obtained for the chloroplast rps2 gene (n = 168) and for the nuclear ITS region (n = 158). Individuals were also genotyped for eight microsatellite loci (n = 176). Interpolation of diversity values was used to visualize the range‐wide distribution of genetic diversity for each of the three marker classes. Minimum spanning networks were constructed showing the relationships between the rps2 and ITS haplotypes, and the geographical distributions of these haplotypes were plotted. The numbers of genetic clusters based on the microsatellite data were estimated using Bayesian clustering approaches. Results The palaeodistribution modelling indicated suitable climate envelopes for M. hypopitys at the LGM in both the Pacific Northwest and south‐eastern USA. High levels of genetic diversity and endemic haplotypes were found in Oregon, the Alexander Archipelago, Wisconsin, and in the south‐eastern part of the species’ distribution range. Main conclusions Our results suggest a complex recolonization history for M. hypopitys in North America, involving persistence in separate eastern and western refugia. A generally high degree of congruence between the different marker classes analysed indicated the presence of multiple refugia, with at least two refugia in each area. In the west, putative refugia were identified in Oregon and the Alexander Archipelago, whereas eastern refugia may have been located in the southern part of the species’ current distribution, as well as in the ‘Driftless Area’. These findings are in contrast to a previous study on the related species Orthilia secunda, which has a similar disjunct distribution to M. hypopitys, but which appears to have recolonized solely from western refugia.     

Species distribution modelling and cpSSR reveal population history of the Neotropical annual herb Tithonia rotundifolia (Asteraceae)

• The impacts of the historical geologic and climatic events on the diversity and genetic structure of Neotropical taxa have recently become a subject of study. However, annual plants associated with tropical dry forests remain under‐studied. The exploration of additional taxa in contrasting environments will improve the current understanding of responses of the Neotropical biota to these events. Here, we explore the species distribution and geographic structure of the annual herb Tithonia rotundifolia.
• We sampled 175 individuals from 19 populations of T. rotundifolia. Species distribution modelling and six microsatellite chloroplast loci were used to infer its population history. We identified areas of historical climate suitability and then tested if there is genetic structuring among these areas.
• Haplotypes showed strong phylogeographic structure. Historical climatic suitability areas were found along the Pacific coast; however, a gap was found at the Isthmus of Tehuantepec (IT). Although Bayesian analysis showed population structuring, amova revealed that the IT is not its main driver. Instead, a subdivision into a higher number of regions had higher F_CT values. Also, populations to the east of the IT showed evidence of recent population expansion and migration in a south–north direction.
• Pleistocene climate fluctuations partially explain the geographic structure of T. rotundifolia. However, life‐history characteristics such as limited seed dispersal and the patchy distribution of suitable habitats explain the high haplotype diversity and population sub‐structuring and diversity. Lastly, the absence of geographic structure of some haplotypes may indicate long‐distance dispersal, or hybridisation with the closely related T. tubaeformis.

     

Significant Asia‐Europe divergence in the middle spotted woodpecker (Aves,Picidae)

Population divergence could be strongly affected by species’ ecology and might not be a direct response to climate‐driven environmental change. We tested this in the middle spotted woodpecker (Dendrocoptes medius), a non‐migratory, low‐dispersal habitat specialist associated with old deciduous forests of the Western Palearctic. We present the first phylogeographic study of this species integrating genetic data (three mitochondrial loci, one autosomal and one Z‐linked intron) with species distribution modelling. Based on this species’ ecology, we predicted that the middle spotted woodpecker could have colonized its current range from multiple Last Glacial Maximum (LGM) refugia and that strongly structured populations could be expected. Indeed, we discovered a strong genetic divergence between Asian and European populations, with a split estimated at around one million of years ago. This was surprising given only slight intraspecific variation in plumage and morphology. Although there was no significant phylogeographic structure within the Asian and European groups, we cannot exclude the possibility of multiple refugia within either group during the LGM. This has to be further investigated with more extensive geographic sampling and larger number of variable independently evolving markers. Future studies should also investigate potential differences in vocalizations and ecology between the two groups. Lineages showing similar level of genetic differentiation including woodpeckers are often treated as species‐level taxa. Comparison of our results with the phylogeographic history of other woodpeckers, suggests that sympatric species with similar life‐histories might have idiosyncratic phylogeographic patterns probably resulting from different ecological requirements or historic stochasticity.     

Quaternary range dynamics of ecologically divergent species (Edraianthus serpyllifolius and E. tenuifolius,Campanulaceae) within the Balkan refugium

Aim We investigated Quaternary range dynamics of two closely related but ecologically divergent species (cold‐tolerant Edraianthus serpyllifolius and thermophilic Edraianthus tenuifolius) with overlapping distribution ranges endemic to the western Balkan Peninsula, an important yet understudied Pleistocene refugium. Our aims were: to test predictions of the ‘refugia‐within‐refugia’ model of strong genetic subdivisions due to population isolation in separate refugia; to explore whether two ecologically divergent species reacted differently to Pleistocene climatic fluctuations; and to test predictions of the displacement refugia model of stronger differentiation among populations in the thermophilic E. tenuifolius compared with the cold‐tolerant E. serpyllifolius. Location The western Balkan Peninsula. Methods We gathered amplified fragment‐length polymorphism (AFLP) data and plastid DNA sequences from two to five individuals from 10 populations of E. serpyllifolius and 22 populations of E. tenuifolius, spanning their entire respective distribution areas. AFLP data were analysed using a Bayesian clustering approach and a distance‐based network approach. Plastid sequences were used to depict relationships among haplotypes in a statistical parsimony network, and to obtain age estimates in a Bayesian framework. Results In E. serpyllifolius, both AFLP and plastid sequence data showed clear geographic structure. Western populations showed high AFLP diversity and a high number of rare fragments. In E. tenuifolius, both markers congruently identified a major phylogeographic split along the lower Neretva valley in central Dalmatia. The most distinct and earliest diverging chloroplast DNA (cpDNA) haplotypes were found further south in the south‐easternmost populations. North‐western populations, identified as a separate cluster by Bayesian clustering, were characterized by low genetic diversity and a low number of rare AFLP markers. Main conclusions Clear evidence for multiple Pleistocene refugia is found not only in the high‐elevation E. serpyllifolius, but also in the lowland E. tenuifolius, despite the lack of obvious dispersal barriers, in line with the refugia‐within‐refugia model. Genealogical relationships and genetic diversity patterns support the hypothesis that cold‐adapted E. serpyllifolius responded to climatic oscillations mostly by elevational range shifts, whereas thermophilic E. tenuifolius did so mainly by latitudinal range shifts, with different phases (and probably extents) of range expansion. In contrast to the displacement refugia hypothesis, the two elevationally differentiated species do not differ in their genetic diversity.     

Genetic evidence for the origins of range disjunctions in the Australian dry sclerophyll plant Hardenbergia violacea

Aim The distribution of genetic variation in the Australian dry sclerophyll plant Hardenbergia violacea (Fabaceae) is examined in the context of Pleistocene climate change in order to identify likely refugia. Particular consideration is given to the origin of range disjunctions in South Australia and Tasmania, and to determining whether the Tasmanian population is indigenous or recently introduced from mainland Australia. Location Southeastern Australian mainland and Tasmania. Methods A combination of chloroplast polymerase chain reaction–restriction fragment length polymorphism and genomic amplified fragment length polymorphism (AFLP) marker systems was used to examine the genetic structure of 292 individuals from 13 populations across the range of H. violacea in southeastern Australia. Results Hardenbergia violacea populations in Tasmania and southern Victoria were characterized by low, almost monotypic chloroplast diversity. New South Wales showed higher haplotype diversity and haplotype sharing among widely distributed populations. Principal coordinates analysis (PCoA) of the AFLP data found a strong latitudinal cline in AFLP variation from northern New South Wales south to Tasmania. The Tasmanian population formed an isolated and somewhat disjunct genetic cluster at one end of this cline. However, the South Australian population was an exception to the clinal variation shown by all other populations, forming a highly disjunct cluster in the PCoA. Within‐population genetic diversity was low in both disjunct populations. Main conclusions The genetic evidence indicates that the Tasmanian population is likely to be indigenous and probably the product of vicariance, which was followed by range contraction at the Last Glacial Maximum or an earlier glacial event. The deep phylogenetic disjunction in South Australia is evidence of a much earlier separation on mainland Australia. The chloroplast structure indicates that, during the Pleistocene, H. violacea underwent broad‐scale recolonization in southern Victoria and Tasmania, possibly from a large continental refugium in eastern New South Wales. We conclude that H. violacea, and presumably the sclerophyll communities in which it occurs, have undergone multiple range contractions to large continental refugia during different Pleistocene glaciations in southeastern Australia.     

History and evolution of alpine plants endemic to the Qinghai-Tibetan Plateau: Aconitum gymnandrum (Ranunculaceae)

How Quaternary climatic oscillations affected range distributions and intraspecific divergence of alpine plants on the Qinghai‐Tibetan Plateau (QTP) remains largely unknown. Here, we report a survey of chloroplast DNA (cpDNA) and nuclear ribosomal internal transcribed spacer (ITS) DNA variation aimed at exploring the phylogeographical history of the QTP alpine endemic Aconitum gymnandrum. We sequenced three cpDNA fragments (rpl20–rps12 intergenic spacer, the trnV intron and psbA‐trnH spacer) and also the nuclear (ITS) region in 245 individuals from 23 populations sampled throughout the species’ range. Two distinct lineages, with eastern and western geographical distributions respectively, were identified from a phylogenetic analysis of ITS sequence variation. Based on a fast substitution rate, these were estimated to have diverged from each other in the early Pleistocene approximately 1.45 Ma. The analysis of cpDNA variation identified nine chlorotypes that clustered into two major clades that were broadly congruent in geographical distribution with the two ITS lineages. The east–west split of cpDNA divergence was supported by an amova which partitioned approximately half of the total variance between these two groups of populations. Analysis of the spatial distribution of chlorotypes showed that each clade was subdivided into two groups of populations such that a total of four population groups existed in the species. It is suggested that these different groups derive from four independent glacial refugia that existed during the Last Glacial Maximum (LGM), and that three of these refugia were located at high altitude on the QTP platform itself at that time. Coalescent simulation of chlorotype genealogies supported both an early Pleistocene origin of the two main cpDNA clades and also the ‘four‐refugia’ hypothesis during the LGM. Two previous phylogeographical studies of QTP alpine plants indicated that such plants retreated to refugia at the eastern/south‐eastern plateau edge during the LGM and/or previous glacial maxima. However, the results for A. gymnandrum suggest that at least some of these cold‐tolerant species may have also survived centrally on the QTP platform throughout the Quaternary.     

Late Pleistocene speciation of three closely related tree peonies endemic to the Qinling–Daba Mountains,a major glacial refugium in Central China

Determining the factors promoting speciation is a major task in ecological and evolutionary research and can be aided by phylogeographic analysis. The Qinling–Daba Mountains (QDM) located in central China form an important geographic barrier between southern subtropical and northern temperate regions, and exhibit complex topography, climatic, and ecological diversity. Surprisingly, few phylogeographic analyses and studies of plant speciation in this region have been conducted. To address this issue, we investigated the genetic divergence and evolutionary histories of three closely related tree peony species (Paeonia qiui, P. jishanensis, and P. rockii) endemic to the QDM. Forty populations of the three tree peony species were genotyped using 22 nuclear simple sequence repeat markers (nSSRs) and three chloroplast DNA sequences to assess genetic structure and phylogenetic relationships, supplemented by morphological characterization and ecological niche modeling (ENM). Morphological and molecular genetic analyses showed the three species to be clearly differentiated from each other. In addition, coalescent analyses using DIYABC conducted on nSSR variation indicated that the species diverged from each other in the late Pleistocene, while ecological niche modeling (ENM) suggested they occupied a larger area during the Last Glacial Maximum (LGM) than at present. The combined genetic evidence from nuclear and chloroplast DNA and the results of ENM indicate that each species persisted through the late Pleistocene in multiple refugia in the Qinling, Daba, and Taihang Mountains with divergence favored by restricted gene flow caused by geographic isolation, ecological divergence, and limited pollen and seed dispersal. Our study contributes to a growing understanding of the origin and population structure of tree peonies and provides insights into the high level of plant endemism present in the Qinling–Daba Mountains of Central China.     

Comparative phylogeography of red maple (Acer rubrum L.) and silver maple (Acer saccharinum L.): impacts of habitat specialization,hybridization and glacial history

Aim We analysed variation in chloroplast DNA (cpDNA) in red maple (Acer rubrum L.) and silver maple (Acer saccharinum L.) across a large part of their geographic ranges. Acer rubrum is one of the most common and morphologically variable deciduous trees of eastern North America, while its sister species A. saccharinum has a more restricted habitat distribution and displays markedly less morphological variation. Our objective was to infer the impact of biogeographic history on cpDNA diversity and phylogeographic structure in both species. Location Deciduous forests of eastern North America. Methods We sequenced 1289 to 1645 bp of non‐coding cpDNA from A. rubrum (n = 258) and A. saccharinum (n = 83). Maximum parsimony networks and spatial analysis of molecular variance (SAMOVA) were used to analyse phylogeographic structure. Rarefaction analyses were used to compare genetic diversity. Results A total of 40 cpDNA haplotypes were recovered from A. rubrum (38 haplotypes) and A. saccharinum (7 haplotypes). Five of the seven A. saccharinum haplotypes were shared with nearby samples of A. rubrum. SAMOVA recovered four phylogeographic groups for A. rubrum in: (1) south‐eastern USA, (2) the Gulf and south‐eastern Coastal Plain, (3) the lower Mississippi River Valley, and (4) the central and northern regions of eastern North America. Acer saccharinum had significantly lower haplotype diversity than A. rubrum, and novel haplotypes in post‐glaciated northern limits of its range were shared with A. rubrum. Main conclusions This is the first study of A. rubrum to report a distinct phylogeographic group centred on the lower Mississippi River, and the first to examine data comparatively with A. saccharinum. We hypothesized that A. rubrum would display stronger phylogeographic structure and greater haplotype diversity than A. saccharinum because of its greater geographic range, and ecological and morphological variation. This hypothesis was supported by the cpDNA analysis. The sharing of cpDNA and chloroplast simple sequence repeat (cpSSR) haplotypes in areas of geographic overlap provides evidence of introgression, which led to an increase in haplotype diversity in both species, and to novel phylogeographic structure in A. rubrum. We recommend that introgression be considered, along with other potential causes, as an explanation for the phylogeographic structure of cpDNA in plants.     

Genetic turnovers and northern survival during the last glacial maximum in European brown bears

The current phylogeographic pattern of European brown bears (Ursus arctos) has commonly been explained by postglacial recolonization out of geographically distinct refugia in southern Europe, a pattern well in accordance with the expansion/contraction model. Studies of ancient DNA from brown bear remains have questioned this pattern, but have failed to explain the glacial distribution of mitochondrial brown bear clades and their subsequent expansion across the European continent. We here present 136 new mitochondrial sequences generated from 346 remains from Europe, ranging in age between the Late Pleistocene and historical times. The genetic data show a high Late Pleistocene diversity across the continent and challenge the strict confinement of bears to traditional southern refugia during the last glacial maximum (LGM). The mitochondrial data further suggest a genetic turnover just before this time, as well as a steep demographic decline starting in the mid‐Holocene. Levels of stable nitrogen isotopes from the remains confirm a previously proposed shift toward increasing herbivory around the LGM in Europe. Overall, these results suggest that in addition to climate, anthropogenic impact and inter‐specific competition may have had more important effects on the brown bear's ecology, demography, and genetic structure than previously thought.     

Patterns of intraspecific diversification in the Piriqueta caroliniana complex in southeastern North America and the Bahamas

Abstract.— We use chloroplast DNA (cpDNA) variation and nested clade phylogeographic analyses to infer the historical processes that have contributed to the high level of morphological and ecological diversification present in a group of herbaceous perennials (the Piriqueta caroliniana complex) in North America and the Bahamas. The presence of morphologically distinct and intercompatible varieties (morphotypes) that can be distinguished based on suites of taxonomic characters (e.g., leaf shape, pubescence type, stature) and contrasting habitat affinities (from marshes to dry pinelands and sand scrub) makes this group particularly appropriate for studies of intraspecific diversification. To examine the distribution of haplotypes among populations, we sampled 467 individuals from 55 locations in Florida, Georgia, and the northern Bahamas (Grand Bahama and Abaco) and screened each individual for cpDNA variation using restriction fragment length polymorphism (RFLP) and heteroduplex analyses. We develop a one-step haplotype phylogeny for this group and use the geographic distributions of haplotypes and clades to test specific phylogeographic hypotheses using the methods developed by Templeton and his colleagues (Templeton 1998). In general, the distribution of haplotypes was strongly influenced by limited dispersal distances, with the more recently derived haplotypes having much lower levels of dispersion and lower frequencies in populations than the ancestral haplotypes. The patterns of clade and haplotype dispersion and displacement and the distribution of morphotypes imply at least three cases of long-distance dispersal and one case of historical fragmentation. The historical patterns inferred for populations of Piriqueta are consistent with known biogeographical events, historical vegetation change, and the concordant patterns of multiple Pleistocene refugia that have been observed for a number of other taxa in southeastern North America.     

Comparative phylogeography of four component species of deciduous broad-leaved forests in Japan based on chloroplast DNA variation

A phylogeographic study of four tree species (Padus grayana, Euonymus oxyphyllus, Magnolia hypoleuca, and Carpinus laxiflora) growing in Japanese deciduous broad-leaved forests was conducted based on chloroplast DNA (cpDNA) variations. Using nucleotide sequences of 702–1,059 bp of intergenic spacers of cpDNA, 20, 27, eight, and eight haplotypes were detected among 251, 251, 226, and 262 individuals sampled from 67, 79, 75, and 71 populations of the above species, respectively. The geographical pattern of the cpDNA variations was highly structured in each species, and the following three regional populations were genetically highly differentiated among all four species: (1) the Sea of Japan-side area, (2) the Kanto region, and (3) southwestern Japan. Based on some interspecific similarities among the phylogeographic patterns, the following migration scenario of Japanese deciduous broad-leaved forests was postulated. During the last glacial maximum (LGM), the forests were separately distributed in six regions. After LGM, as the climate warmed, the forests in eastern Japan separately expanded from each of the refugia along the Sea of Japan-side or along the Pacific Ocean-side. In contrast, those in southwestern Japan retreated and moved to high altitudes from each of the continuous forests.