Knowledge of the historical distribution and postglacial phylogeography and evolution of a species is important to better understand its current distribution and population structure and potential fate in the future, especially under climate change conditions, and conservation of its genetic resources. We have addressed this issue in a wide-ranging and heavily exploited keystone forest tree species of eastern North America, eastern white pine (Pinus strobus). We examined the range-wide population genetic structure, tested various hypothetical population history and evolutionary scenarios and inferred the location of glacial refugium and post-glacial recolonization routes. Our hypothesis was that eastern white pine survived in a single glacial refugium and expanded through multiple post-glacial recolonization routes.
We studied the range-wide genetic diversity and population structure of 33 eastern white pine populations using 12 nuclear and 3 chloroplast microsatellite DNA markers. We used Approximate Bayesian Computation approach to test various evolutionary scenarios. We observed high levels of genetic diversity, and significant genetic differentiation (FST = 0.104) and population structure among eastern white pine populations across its range. A south to north trend of declining genetic diversity existed, consistent with repeated founder effects during post-glaciation migration northwards. We observed broad consensus from nuclear and chloroplast genetic markers supporting the presence of two main post-glacial recolonization routes that originated from a single southern refugium in the mid-Atlantic plain. One route gave rise to populations at the western margin of the species’ range in Minnesota and western Ontario. The second route gave rise to central-eastern populations, which branched into two subgroups: central and eastern. We observed minimal sharing of chloroplast haplotypes between recolonization routes but there was evidence of admixture between the western and west-central populations.
Our study reveals a single southern refugium, two recolonization routes and three genetically distinguishable lineages in eastern white pine that we suggest to be treated as separate Evolutionarily Significant Units. Like many wide-ranging North American species, eastern white pine retains the genetic signatures of post-glacial recolonization and evolution, and its contemporary population genetic structure reflects not just the modern distribution and effects of heavy exploitation but also routes northward from its glacial refugium.
Electronic supplementary material
The online version of this article (doi:10.1186/s12862-016-0624-1) contains supplementary material, which is available to authorized users.
Phylogeography; Population genetic structure; post-glacial migration; Pinus strobus; Genetic signatures of heavy exploitation; Molecular evolution; Evolutionarily significant units
Understanding the history of forests and their species' demographic responses to past disturbances is important for predicting impacts of future environmental changes. Tropical rainforests of the Guineo-Congolian region in Central Africa are believed to have survived the Pleistocene glacial periods in a few major refugia, essentially centred on mountainous regions close to the Atlantic Ocean. We tested this hypothesis by investigating the phylogeographic structure of a widespread, ancient rainforest tree species, Symphonia globulifera L. f. (Clusiaceae), using plastid DNA sequences (chloroplast DNA [cpDNA], psbA-trnH intergenic spacer) and nuclear microsatellites (simple sequence repeats, SSRs). SSRs identified four gene pools located in Benin, West Cameroon, South Cameroon and Gabon, and São Tomé. This structure was also apparent at cpDNA. Approximate Bayesian Computation detected recent bottlenecks approximately dated to the last glacial maximum in Benin, West Cameroon and São Tomé, and an older bottleneck in South Cameroon and Gabon, suggesting a genetic effect of Pleistocene cycles of forest contraction. CpDNA haplotype distribution indicated wide-ranging long-term persistence of S. globulifera both inside and outside of postulated forest refugia. Pollen flow was four times greater than that of seed in South Cameroon and Gabon, which probably enabled rapid population recovery after bottlenecks. Furthermore, our study suggested ecotypic differentiation—coastal or swamp vs terra firme—in S. globulifera. Comparison with other tree phylogeographic studies in Central Africa highlighted the relevance of species-specific responses to environmental change in forest trees.
forest refuge theory; phylogeography; psbA-trnH; SSR; tropical Africa; approximate bayesian computation
Geological events in the latter Cenozoic have influenced the distribution, abundance and genetic structure of tree populations in temperate and tropical North America. The biogeographical history of temperate vegetation that spans large ranges of latitude is complex, involving multiple latitudinal shifts that might have occurred via different migration routes. We determined the regional structuring of genetic variation of sugar maple (Acer saccharum subsp. saccharum) and its only subspecies in tropical America (Acer saccharum subsp. skutchii) using nuclear and chloroplast data. The studied populations span a geographic range from Maine, USA (46°N), to El Progreso, Guatemala (15°N). We examined genetic subdivisions, explored the locations of ancestral haplotypes, analyzed genetic data to explore the presence of a single or multiple glacial refugia, and tested whether genetic lineages are temporally consistent with a Pleistocene or older divergence.
Nuclear and chloroplast data indicated that populations in midwestern USA and western Mexico were highly differentiated from populations in the rest of the sites. The time of the most recent common ancestor of the western Mexico haplotype lineage was dated to the Pliocene (5.9 Ma, 95 % HPD: 4.3–7.3 Ma). Splits during the Pleistocene separated the rest of the phylogroups. The most frequent and widespread haplotype occurred in half of the sites (Guatemala, eastern Mexico, southeastern USA, and Ohio). Our data also suggested that multiple Pleistocene refugia (tropics-southeastern USA, midwestern, and northeastern USA), but not western Mexico (Jalisco), contributed to post-glacial northward expansion of ranges. Current southern Mexican and Guatemalan populations have reduced population sizes, genetic bottlenecks and tend toward homozygosity, as indicated using nuclear and chloroplast markers.
The divergence of western Mexican populations from the rest of the sugar maples likely resulted from orographic and volcanic barriers to gene flow. Past connectivity among populations in the southeastern USA and eastern Mexico and Guatemala possible occurred through gene flow during the Pleistocene. The time to the most common ancestor values revealed that populations from the Midwest and Northeast USA represented different haplotype lineages, indicating major divergence of haplotypes lineages before the Last Glacial Maximum and suggesting the existence of multiple glacial refugia.
Electronic supplementary material
The online version of this article (doi:10.1186/s12862-015-0518-7) contains supplementary material, which is available to authorized users.
Acer saccharum; Cloud forest; Divergence time; Gene flow; Genetic diversity; Last glacial maximum; Microsatellite; Phylogeography; Pleistocene; Pliocene; Range edge
Background and Aims
Few phylogeographic studies have been undertaken of species confined to narrow, linear coastal systems where past sea level and geomorphological changes may have had a profound effect on species population sizes and distributions. In this study, a phylogeographic analysis was conducted of Eucalyptus gomphocephala (tuart), a tree species restricted to a 400 × 10 km band of coastal sand-plain in south west Australia. Here, there is little known about the response of coastal vegetation to glacial/interglacial climate change, and a test was made as to whether this species was likely to have persisted widely through the Last Glacial Maximum (LGM), or conforms to a post-LGM dispersal model of recovery from few refugia.
The genetic structure over the entire range of tuart was assessed using seven nuclear (21 populations; n = 595) and four chloroplast (24 populations; n = 238) microsatellite markers designed for eucalypt species. Correlative palaeodistribution modelling was also conducted based on five climatic variables, within two LGM models.
The chloroplast markers generated six haplotypes, which were strongly geographically structured (GST = 0·86 and RST = 0·75). Nuclear microsatellite diversity was high (overall mean HE 0·75) and uniformly distributed (FST = 0·05), with a strong pattern of isolation by distance (r2 = 0·362, P = 0·001). Distribution models of E. gomphocephala during the LGM showed a wide distribution that extended at least 30 km westward from the current distribution to the palaeo-coastline.
The chloroplast and nuclear data suggest wide persistence of E. gomphocephala during the LGM. Palaeodistribution modelling supports the conclusions drawn from genetic data and indicates a widespread westward shift of E. gomphocephala onto the exposed continental shelf during the LGM. This study highlights the importance of the inclusion of complementary, non-genetic data (information on geomorphology and palaeoclimate) to interpret phylogeographic patterns.
Australian biogeography; climate change; coastal geomorphology; Eucalyptus gomphocephala; founder effects; Last Glacial Maximum; LGM; microsatellites; Myrtaceae; palaeodistribution modelling; phylogeography; southern hemisphere; south-western Australia; tuart
Background and Aims
Dendrobium species show tremendous morphological diversity and have broad geographical distribution. As repetitive sequence analysis is a useful tool to investigate the evolution of chromosomes and genomes, the aim of the present study was the characterization of repetitive sequences from Dendrobium moschatum for comparative molecular and cytogenetic studies in the related species Dendrobium aphyllum, Dendrobium aggregatum and representatives from other orchid genera.
In order to isolate highly repetitive sequences, a c0t-1 DNA plasmid library was established. Repeats were sequenced and used as probes for Southern hybridization. Sequence divergence was analysed using bioinformatic tools. Repetitive sequences were localized along orchid chromosomes by fluorescence in situ hybridization (FISH).
Characterization of the c0t-1 library resulted in the detection of repetitive sequences including the (GA)n dinucleotide DmoO11, numerous Arabidopsis-like telomeric repeats and the highly amplified dispersed repeat DmoF14. The DmoF14 repeat is conserved in six Dendrobium species but diversified in representative species of three other orchid genera. FISH analyses showed the genome-wide distribution of DmoF14 in D. moschatum, D. aphyllum and D. aggregatum. Hybridization with the telomeric repeats demonstrated Arabidopsis-like telomeres at the chromosome ends of Dendrobium species. However, FISH using the telomeric probe revealed two pairs of chromosomes with strong intercalary signals in D. aphyllum. FISH showed the terminal position of 5S and 18S–5·8S–25S rRNA genes and a characteristic number of rDNA sites in the three Dendrobium species.
The repeated sequences isolated from D. moschatum c0t-1 DNA constitute major DNA families of the D. moschatum, D. aphyllum and D. aggregatum genomes with DmoF14 representing an ancient component of orchid genomes. Large intercalary telomere-like arrays suggest chromosomal rearrangements in D. aphyllum while the number and localization of rRNA genes as well as the species-specific distribution pattern of an abundant microsatellite reflect the genomic diversity of the three Dendrobium species.
Orchidaceae; Dendrobium moschatum; Dendrobium aphyllum; Dendrobium aggregatum; repetitive DNA; FISH; c0t-1 DNA
The strawberry tree diverged into two groups during the Quaternary, but before the LGM, surviving in refugia located in the western end of the Mediterranean region and with the Eastern Mediterranean being colonized more recently. This migration was possible because Europe and North Africa were occasionally connected through the straits of Gibraltar and Sicily. Likewise, our evidence supports arrival in Ireland from northern Iberia in post-glacial times. Altogether, our results reveal the considerable ability of the strawberry tree for dispersal, allowing it to migrate over thousands of kilometres and cross stretches of sea, which may be crucial for its future survival.
Phylogeographic studies give us the opportunity to reconstruct the historical migrations of species and link them with climatic and geographic variation. They are, therefore, a key tool to understanding the relationships among biology, geology and history. One of the most interesting biogeographical areas of the world is the Mediterranean region. However, in this area, the description of concordant phylogeographic patterns is quite scarce, which limits the understanding of evolutionary patterns related to climate. Species with one-dimensional distribution ranges, such as the strawberry tree (Arbutus unedo), are particularly useful to unravel these patterns. Here, we describe its phylogeographic structure and check for concordance with patterns seen in other Mediterranean plants: longitudinal/latitudinal clines of diversity, evidence for glacial refugia and the role of sea straits in dispersal. We also identify the most likely source for the disjunct Irish population. With this aim, we sequenced four chloroplast non-coding fragments of A. unedo from 23 populations covering its whole distribution. We determined the genetic diversity, population structure, haplotype genealogy and time to the most recent common ancestor. The genealogy revealed two clades that separated during the last 700 ky but before the last glacial maximum. One clade occupies Atlantic Iberia and North Africa, while the other occurs in the Western Mediterranean. The Eastern Mediterranean is inhabited by newer haplotypes derived from both clades, while the Irish population is closely related to Iberian demes. The straits of Sicily and Gibraltar partially restricted the gene flow. We concluded that a vicariance event during the Late Quaternary in the western end of the species' range followed by eastward migration seems a likely explanation for the observed phylogeographic pattern. The role of straits indicates an occasional communication between Europe and North Africa, suggesting that the latter was a novel refugia. The East–West genetic split in Iberia is consistent with the refugia-within-refugia model. Finally, the strawberry tree possibly reached Ireland from Iberia instead of throughout the maritime fringe of France as previously thought.
Arbutus unedo; clades; cpDNA; Lusitanian; Mediterranean; refugia; strait; vicariance
Background and Aims
Under predicted climate change scenarios, increased temperatures are likely to predispose trees to leaf and other tissue damage, resulting in plant death and contraction of already narrow distribution ranges in many relictual species. The effects of predicted upward temperatures may be further exacerbated by changes in rainfall patterns and damage caused by frosts on trees that have been insufficiently cold-hardened. The Araucariaceae is a relictual family and the seven species found in Australia have limited natural distributions characterized by low frost intensity and frequency, and warm summer temperatures. The temperature limits for these species were determined in order to help understand how such species will fare in a changing climate.
Experiments were conducted using samples from representative trees of the Araucariaceae species occurring in Australia, Agathis (A. atropurpurea, A. microstachya and A. robusta), Arauacaria (A. bidwilli, A. cunninghamii and A. heterophylla) and Wollemia nobilis. Samples were collected from plants grown in a common garden environment. Lower and higher temperature limits were determined by subjecting detached winter-hardened leaves to temperatures from 0 to –17 °C and summer-exposed leaves to 25 to 63 °C, then measuring the efficiency of photosystem II (Fv/Fm) and visually rating leaf damage. The exotherm, a sharp rise in temperature indicating the point of ice nucleation within the cells of the leaf, was measured on detached leaves of winter-hardened and summer temperature-exposed leaves.
Lower temperature limits (indicated by FT50, the temperature at which PSII efficiency is 50 %, and LT50 the temperature at which 50 % visual leaf damage occurred) were approx. –5·5 to –7·5 °C for A. atropurpurea, A. microstachya and A. heterophylla, approx. –7 to –9 °C for A. robusta, A. bidwillii and A. cunninghamii, and –10·5 to –11 °C for W. nobilis. High temperature damage began at 47·5 °C for W. nobilis, and occurred in the range 48·5–52 °C for A. bidwillii and A. cunninghamii, and in the range 50·5–53·5 °C for A. robusta, A. microstachya and A. heterophylla. Winter-hardened leaves had ice nucleation temperatures of –5·5 °C or lower, with W. nobilis the lowest at –6·8 °C. All species had significantly higher ice nucleation temperatures in summer, with A. atropurpurea and A. heterophylla forming ice in the leaf at temperatures >3 °C higher in summer than in winter. Wollemia nobilis had lower FT50 and LT50 values than its ice nucleation temperature, indicating that the species has a degree of ice tolerance.
While lower temperature limits in the Australian Araucariaceae are generally unlikely to affect their survival in wild populations during normal winters, unseasonal frosts may have devastating effects on tree survival. Extreme high temperatures are not common in the areas of natural occurrence, but upward temperature shifts, in combination with localized radiant heating, may increase the heat experienced within a canopy by at least 10 °C and impact on tree survival, and may contribute to range contraction. Heat stress may explain why many landscape plantings of W. nobilis have failed in hotter areas of Australia.
Hardiness; frost; high temperature; climate change; Araucariaceae; rain forest
At the last glacial maximum, Fennoscandia was covered by an ice sheet while the tundra occupied most of the rest of northern Eurasia. More or less disjunct refugial populations of plants were dispersed in southern Europe, often trapped between mountain ranges and seas. Genetic and paleobotanical evidences indicate that these populations have contributed much to Holocene recolonization of more northern latitudes. Less supportive evidence has been found for the existence of glacial populations located closer to the ice margin. Scots pine (Pinus sylvestris L.) is a nordic conifer with a wide natural range covering much of Eurasia. Fractures in its extant genetic structure might be indicative of glacial vicariance and how different refugia contributed to the current distribution at the continental level. The population structure of Scots pine was investigated on much of its Eurasian natural range using maternally inherited mitochondrial DNA polymorphisms.
A novel polymorphic region of the Scots pine mitochondrial genome has been identified, the intron 1 of nad7, with three variants caused by insertions-deletions. From 986 trees distributed among 54 populations, four distinct multi-locus mitochondrial haplotypes (mitotypes) were detected based on the three nad7 intron 1 haplotypes and two previously reported size variants for nad1 intron B/C. Population differentiation was high (GST = 0.657) and the distribution of the mitotypes was geographically highly structured, suggesting at least four genetically distinct ancestral lineages. A cosmopolitan lineage was widely distributed in much of Europe throughout eastern Asia. A previously reported lineage limited to the Iberian Peninsula was confirmed. A new geographically restricted lineage was found confined to Asia Minor. A new lineage was restricted to more northern latitudes in northeastern Europe and the Baltic region.
The contribution of the various ancestral lineages to the current distribution of Scots pine was asymmetric and extant endemism reflected the presence of large geographic barriers to migration. The results suggest a complex biogeographical history with glacial refugia shared with temperate plant species in southern European Peninsulas and Asia Minor, and a genetically distinct glacial population located more North. These results confirm recent observations for cold tolerant species about the possible existence of refugial populations at mid-northern latitudes contributing significantly to the recolonization of northern Europe. Thus, Eurasian populations of nordic plant species might not be as genetically homogenous as assumed by simply considering them as offsets of glacial populations located in southern peninsulas. As such, they might have evolved distinctive genetic adaptations during glacial vicariance, worth evaluating and considering for conservation.
Molecular data from three chloroplast markers resolve individuals attributable to Radula buccinifera in six lineages belonging to two subgenera, indicating the species is polyphyletic as currently circumscribed. All lineages are morphologically diagnosable, but one pair exhibits such morphological overlap that they can be considered cryptic. Molecular and morphological data justify the re-instatement of a broadly circumscribed ecologically variable R. strangulata, of R. mittenii, and the description of five new species. Two species Radula mittenii Steph. and R. notabilis
sp. nov. are endemic to the Wet Tropics Bioregion of north-east Queensland, suggesting high diversity and high endemism might characterise the bryoflora of this relatively isolated wet-tropical region. Radula demissa
sp. nov. is endemic to southern temperate Australasia, and like R. strangulata occurs on both sides of the Tasman Sea. Radula imposita
sp. nov. is a twig and leaf epiphyte found in association with waterways in New South Wales and Queensland. Another species, R. pugioniformis
sp. nov., has been confused with Radula buccinifera but was not included in the molecular phylogeny. Morphological data suggest it may belong to subg. Odontoradula. Radula buccinifera is endemic to Australia including Western Australia and Tasmania, and to date is known from south of the Clarence River on the north coast of New South Wales. Nested within R. buccinifera is a morphologically distinct plant from Norfolk Island described as R. anisotoma
Radula australiana is resolved as monophyletic, sister to a species occurring in east coast Australian rainforests, and nesting among the R. buccinifera lineages with strong support. The molecular phylogeny suggests several long-distance dispersal events may have occurred. These include two east-west dispersal events from New Zealand to Tasmania and south-east Australia in R. strangulata, one east-west dispersal event from Tasmania to Western Australia in R. buccinifera, and at least one west-east dispersal from Australia to New Zealand in R. australiana. Another west-east dispersal event from Australia to Norfolk Island may have led to the budding speciation of R. anisotoma. In contrast, Radula demissa is phylogeographically subdivided into strongly supported clades either side of the Tasman Sea, suggesting long distance dispersal is infrequent in this species.
Radulaceae; Radula subg. Metaradula; Radula anisotoma; Radula australiana; Radula buccinifera; Radula demissa; Radula imposita; Radula mittenii; Radula notabilis; Radula pugioniformis; Radula strangulata; sp. nov.; morphology; DNA sequence data; Australia; New Zealand; Flora; liverwort; dispersal; biogeography; cryptic species
The aim of this research was to identify the effects of Pleistocene climate change on the distribution of fauna in Tasmania, and contrast this with biotic responses in other temperate regions in the Northern and Southern Hemisphere that experienced glacial activity during this epoch. This was achieved by examining the phylogeographic patterns in a widely distributed Tasmanian endemic reptile, Niveoscincus ocellatus. 204 individuals from 29 populations across the distributional range of N. ocellatus were surveyed for variation at two mitochondrial genes (ND2, ND4), and two nuclear genes (β-globin, RPS8). Phylogenetic relationships were reconstructed using a range of methods (maximum parsimony, Bayesian inference and haplotype networks), and the demographic histories of populations were assessed (AMOVA, Tajima’s D, Fu’s Fs, mismatch distributions, extended Bayesian skyline plots, and relaxed random walk analyses).
There was a high degree of mitochondrial haplotype diversity (96 unique haplotypes) and phylogeographic structure, where spatially distinct groups were associated with Tasmania’s Northeast and a large area covering Southeast and Central Tasmania. Phylogeographic structure was also present within each major group, but the degree varied regionally, being highest in the Northeast. Only the Southeastern group had a signature of demographic expansion, occurring during the Pleistocene but post-dating the Last Glacial Maximum. In contrast, nuclear DNA had low levels of variation and a lack of phylogeographic structure, and further loci should be surveyed to corroborate the mitochondrial inferences.
The phylogeographic patterns of N. ocellatus indicate Pleistocene range and demographic expansion in N. ocellatus, particularly in the Southeast and Central areas of Tasmania. Expansion in Central and Southeastern areas appears to have been more recent in both demographic and spatial contexts, than in Northeast Tasmania, which is consistent with inferences for other taxa of greater stability and persistence in Northeast Tasmania during the Last Glacial Maximum. These phylogeographic patterns indicate contrasting demographic histories of populations in close proximity to areas directly affected by glaciers in the Southern Hemisphere during the LGM.
Electronic supplementary material
The online version of this article (doi:10.1186/s12862-015-0397-y) contains supplementary material, which is available to authorized users.
Glacial refugia; Pleistocene; Tasmania; Phylogeography; Reptile; Recolonisation
Identify patterns of change in species distributions, diversity, concentrations of evolutionary history, and assembly of Australian rainforests.
We used the distribution records of all known rainforest woody species in Australia across their full continental extent. These were analysed using measures of species richness, phylogenetic diversity (PD), phylogenetic endemism (PE) and phylogenetic structure (net relatedness index; NRI). Phylogenetic structure was assessed using both continental and regional species pools. To test the influence of growth-form, freestanding and climbing plants were analysed independently, and in combination.
Species richness decreased along two generally orthogonal continental axes, corresponding with wet to seasonally dry and tropical to temperate habitats. The PE analyses identified four main areas of substantially restricted phylogenetic diversity, including parts of Cape York, Wet Tropics, Border Ranges, and Tasmania. The continental pool NRI results showed evenness (species less related than expected by chance) in groups of grid cells in coastally aligned areas of species rich tropical and sub-tropical rainforest, and in low diversity moist forest areas in the south-east of the Great Dividing Range and in Tasmania. Monsoon and drier vine forests, and moist forests inland from upland refugia showed phylogenetic clustering, reflecting lower diversity and more relatedness. Signals for evenness in Tasmania and clustering in northern monsoon forests weakened in analyses using regional species pools. For climbing plants, values for NRI by grid cell showed strong spatial structuring, with high diversity and PE concentrated in moist tropical and subtropical regions.
Concentrations of rainforest evolutionary history (phylo-diversity) were patchily distributed within a continuum of species distributions. Contrasting with previous concepts of rainforest community distribution, our findings of continuous distributions and continental connectivity have significant implications for interpreting rainforest evolutionary history and current day ecological processes, and for managing rainforest diversity in changing circumstances.
The aim of the present study was to use a combined phylogeographic and species distribution modelling approach to compare the glacial histories of two plant species with overlapping distributions, Orthilia secunda (one-sided wintergreen) and Monotropa hypopitys (yellow bird's nest). Phylogeographic analysis was carried out to determine the distribution of genetic variation across the range of each species and to test whether both correspond to the "classic" model of high diversity in the south, with decreasing diversity at higher latitudes, or whether the cold-adapted O. secunda might retain more genetic variation in northern populations. In addition, projected species distributions based on a future climate scenario were modelled to assess how changes in the species ranges might impact on total intraspecific diversity in both cases.
Palaeodistribution modelling and phylogeographic analysis using multiple genetic markers (chloroplast trnS-trnG region, nuclear ITS and microsatellites for O. secunda; chloroplast rps2, nuclear ITS and microsatellites for M. hypopitys) indicated that both species persisted throughout the Last Glacial Maximum in southern refugia. For both species, the majority of the genetic diversity was concentrated in these southerly populations, whereas those in recolonized areas generally exhibited lower levels of diversity, particularly in M. hypopitys. Species distribution modelling based on projected future climate indicated substantial changes in the ranges of both species, with a loss of southern and central populations, and a potential northward expansion for the temperate M. hypopitys.
Both Orthilia secunda and Monotropa hypopitys appear to have persisted through the LGM in Europe in southern refugia. The boreal O. secunda, however, has retained a larger proportion of its genetic diversity in more northerly populations outside these refugial areas than the temperate M. hypopitys. Given that future species distribution modelling suggests northern range shifts and loss of suitable habitat in the southern parts of the species' current distributions, extinction of genetically diverse rear edge populations could have a significant effect in the rangewide intraspecific diversity of both species, but particularly in M. hypopitys.
The ranges of arctic–alpine species have shifted extensively with Pleistocene climate changes and glaciations. Using sequence data from the trnH-psbA and trnT-trnL chloroplast DNA spacer regions, we investigated the phylogeography of the widespread, ancient (>3 million years) arctic–alpine plant Oxyria digyna (Polygonaceae). We identified 45 haplotypes and six highly divergent major lineages; estimated ages of these lineages (time to most recent common ancestor, TMRCA) ranged from ∼0.5 to 2.5 million years. One lineage is widespread in the arctic, a second is restricted to the southern Rocky Mountains of the western United States, and a third was found only in the Himalayan and Altai regions of Asia. Three other lineages are widespread in western North America, where they overlap extensively. The high genetic diversity and the presence of divergent major cpDNA lineages within Oxyria digyna reflect its age and suggest that it was widespread during much of its history. The distributions of individual lineages indicate repeated spread of Oxyria digyna through North America over multiple glacial cycles. During the Last Glacial Maximum it persisted in multiple refugia in western North America, including Beringia, south of the continental ice, and within the northern limits of the Cordilleran ice sheet. Our data contribute to a growing body of evidence that arctic–alpine species have migrated from different source regions over multiple glacial cycles and that cryptic refugia contributed to persistence through the Last Glacial Maximum.
Arctic-alpine plants; cpDNA; phylogeography; Pleistocene glaciations; refugia
The evolution of the savanna biome has been deeply marked by repeated contraction/expansion phases due to climate perturbations during the Quaternary period. In this study, we investigated the impact of the last glacial maximum (LGM) on the present genetic pattern of Vitellaria paradoxa (shea tree), a major African savanna tree. A range-wide sampling of the species enabled us to sample 374 individuals from 71 populations distributed throughout sub-Sahelian Africa. Trees were genotyped using 3 chloroplasts and 12 nuclear microsatellites, and were sequenced for 2 polymorphic chloroplast intergenic spacers. Analyses of genetic diversity and structure were based on frequency-based and Bayesian methods. Potential distributions of V. paradoxa at present, during the LGM and the last interglacial period, were examined using DIVA-GIS ecological niche modelling (ENM). Haplotypic and allelic richness varied significantly across the range according to chloroplast and nuclear microsatellites, which pointed to higher diversity in West Africa. A high but contrasted level of differentiation was revealed among populations with a clear phylogeographic signal, with both nuclear (FST=0.21; RST=0.28; RST>RST (permuted)) and chloroplast simple sequence repeats (SSRs) (GST=0.81; NST=0.90; NST>NST (permuted)). We identified a strong geographically related structure separating western and eastern populations, and a substructure in the eastern part of the area consistent with subspecies distinction. Using ENM, we deduced that perturbations during the LGM fragmented the potential eastern distribution of shea tree, but not its distribution in West Africa. Our main results suggest that climate variations are the major factor explaining the genetic pattern of V. paradoxa.
intergenic spacers; microsatellite markers; last glacial maximum; phylogeography; sub-Sahelian Africa; Vitellaria paradoxa
Background and Aims
The widely accepted paradigm that the modern genetic structure of plant species in the northern hemisphere has been largely determined by recolonization from refugia after the last glacial maximum fails to explain the presence of cold-tolerant species at intermediate latitudes. Another generally accepted paradigm is that mountain ridges act as important barriers causing genetic isolation of species, but this too has been challenged in recent studies. The aims of the work reported here were to determine the genetic diversity and distribution patterns of extant natural populations of an endangered cool temperate species, Faxinus mandshurica, and to examine whether these two paradigms are appropriate when applied to this species over a wide geographical scale.
1435 adult individuals were sampled from 30 natural populations across the main and central range of the species, covering major mountain ranges across North-east China (NEC). Genetic variation was estimated based on nine polymorphic nuclear microsatellite loci. Phylogeographical analyses were employed using various approaches, including Bayesian clustering, spatial analysis of molecular variance, Monmonier's algorithm, neighbor-joining trees, principal co-ordinate analysis and isolation by distance.
Genetic diversity within populations was relatively high, and no significant recent bottlenecks were detected in any of the populations. A significant negative correlation between intra-population genetic diversity and latitude was identified. In contrast, genetic differentiation among all the populations examined was extremely low and no clear geographic genetic structure was identified, with the exception of one distinct population.
The modern genetic structure in this species can be explained by extensive gene flow, an absence of mountains acting as barriers, and the presence of a wide refuge across NEC rather than multiple small refugia. Intra-population genetic variation along latitudes is probably associated with the systematically northward shifts of forest biomes in eastern China during the mid-Holocene. To determine important genetic patterns and identify resources for conservation, however, it will be necessary to examine differentially inherited genetic markers exposed to selection pressures (e.g. chloroplast DNA) and to investigate different generations.
Fraxinus mandshurica; nuclear microsatellites; latitude variation; historical migration; fossil pollen; spatial genetic structure; genetic barriers
The biogeographical relationships between far-separated populations, in particular, those in the mainland and islands, remain unclear for widespread species in eastern Asia where the current distribution of plants was greatly influenced by the Quaternary climate. Deciduous Oriental oak (Quercus variabilis) is one of the most widely distributed species in eastern Asia. In this study, leaf material of 528 Q. variabilis trees from 50 populations across the whole distribution (Mainland China, Korea Peninsular as well as Japan, Zhoushan and Taiwan Islands) was collected, and three cpDNA intergenic spacer fragments were sequenced using universal primers. A total of 26 haplotypes were detected, and it showed a weak phylogeographical structure in eastern Asia populations at species level, however, in the central-eastern region of Mainland China, the populations had more haplotypes than those in other regions, with a significant phylogeographical structure (NST = 0.751> GST = 0.690, P<0.05). Q. variabilis displayed high interpopulation and low intrapopulation genetic diversity across the distribution range. Both unimodal mismatch distribution and significant negative Fu’s FS indicated a demographic expansion of Q. variabilis populations in East Asia. A fossil calibrated phylogenetic tree showed a rapid speciation during Pleistocene, with a population augment occurred in Middle Pleistocene. Both diversity patterns and ecological niche modelling indicated there could be multiple glacial refugia and possible bottleneck or founder effects occurred in the southern Japan. We dated major spatial expansion of Q. variabilis population in eastern Asia to the last glacial cycle(s), a period with sea-level fluctuations and land bridges in East China Sea as possible dispersal corridors. This study showed that geographical heterogeneity combined with climate and sea-level changes have shaped the genetic structure of this wide-ranging tree species in East Asia.
Interglacial-glacial cycles of the Quaternary are widely recognized in shaping phylogeographic structure. Patterns from cold adapted species can be especially informative - in particular, uncovering additional glacial refugia, identifying likely recolonization patterns, and increasing our understanding of species’ responses to climate change. We investigated phylogenetic structure of the wolverine, a wide-ranging cold adapted carnivore, using a 318 bp of the mitochondrial DNA control region for 983 wolverines (n = 209 this study, n = 774 from GenBank) from across their full Holarctic distribution. Bayesian phylogenetic tree reconstruction and the distribution of observed pairwise haplotype differences (mismatch distribution) provided evidence of a single rapid population expansion across the wolverine’s Holarctic range. Even though molecular evidence corroborated a single refugium, significant subdivisions of population genetic structure (0.01< ΦST <0.99, P<0.05) were detected. Pairwise ΦST estimates separated Scandinavia from Russia and Mongolia, and identified five main divisions within North America - the Central Arctic, a western region, an eastern region consisting of Ontario and Quebec/Labrador, Manitoba, and California. These data are in contrast to the nearly panmictic structure observed in northwestern North America using nuclear microsatellites, but largely support the nuclear DNA separation of contemporary Manitoba and Ontario wolverines from northern populations. Historic samples (c. 1900) from the functionally extirpated eastern population of Quebec/Labrador displayed genetic similarities to contemporary Ontario wolverines. To understand these divergence patterns, four hypotheses were tested using Approximate Bayesian Computation (ABC). The most supported hypothesis was a single Beringia incursion during the last glacial maximum that established the northwestern population, followed by a west-to-east colonization during the Holocene. This pattern is suggestive of colonization occurring in accordance with glacial retreat, and supports expansion from a single refugium. These data are significant relative to current discussions on the conservation status of this species across its range.
The mesic habitats of eastern Australia harbour a highly diverse fauna. We examined the impact of climatic oscillations and recognised biogeographic barriers on the evolutionary history of the delicate skink (Lampropholis delicata), a species that occurs in moist habitats throughout eastern Australia. The delicate skink is a common and widespread species whose distribution spans 26° of latitude and nine major biogeographic barriers in eastern Australia. Sequence data were obtained from four mitochondrial genes (ND2, ND4, 12SrRNA, 16SrRNA) for 238 individuals from 120 populations across the entire native distribution of the species. The evolutionary history and diversification of the delicate skink was investigated using a range of phylogenetic (Maximum Likelihood, Bayesian) and phylogeographic analyses (genetic diversity, ΦST, AMOVA, Tajima's D, Fu's F statistic).
Nine geographically structured, genetically divergent clades were identified within the delicate skink. The main clades diverged during the late Miocene-Pliocene, coinciding with the decline and fragmentation of rainforest and other wet forest habitats in eastern Australia. Most of the phylogeographic breaks within the delicate skink were concordant with dry habitat or high elevation barriers, including several recognised biogeographic barriers in eastern Australia (Burdekin Gap, St Lawrence Gap, McPherson Range, Hunter Valley, southern New South Wales). Genetically divergent populations were also located in high elevation topographic isolates inland from the main range of L. delicata (Kroombit Tops, Blackdown Tablelands, Coolah Tops). The species colonised South Australia from southern New South Wales via an inland route, possibly along the Murray River system. There is evidence for recent expansion of the species range across eastern Victoria and into Tasmania, via the Bassian Isthmus, during the late Pleistocene.
The delicate skink is a single widespread, but genetically variable, species. This study provides the first detailed phylogeographic investigation of a widespread species whose distribution spans virtually all of the major biogeographic barriers in eastern Australia.
Tasmania is a natural laboratory for investigating the evolutionary processes of the Quaternary. It is a large island lying 40-44 degrees S, which was repeatedly glaciated and linked to southeastern continental Australia during the Quaternary. Climate change promoted both the isolation of species in glacial refugia, and an exchange between Tasmanian and mainland floras. Eucalyptus is a complex and diverse genus, which has increased in abundance in Australia over the past 100 kyr, probably in response to higher fire frequency. Morphological evidence suggests that gene flow may have occurred between many eucalypt species after changes in their distribution during the Quaternary. This paper summarizes recent genetic evidence for migration and introgressive hybridization in Tasmanian Eucalyptus. Maternally inherited chloroplast DNA reveals a long-term persistence of eucalypts in southeastern Tasmanian refugia, coupled with introgressive hybridization involving many species. Detailed analysis of the widespread species Eucalyptus globulus suggests that migration from mainland Australia was followed by introgression involving a rare Tasmanian endemic. The data support the hypothesis that changes in distribution of interfertile species during the Quaternary have promoted reticulate evolution in Eucalyptus.
East Asia has the most diverse temperate flora in the world primarily due to the lack of Pleistocene glaciation and the geographic heterogeneity. Although increasing phylogeography studies in this region provided more proofs in this issue, discrepancies and uncertainty still exist, especially in northern temperate deciduous broad-leaved and coniferous mixed forest region (II). And a widespread plant species could reduce the complexity to infer the relationship between diversity and physiographical pattern. Hence, we studied the evolution history of a widespread temperate tree, Acer mono, populations in region II and the influence of physiographic patterns on intraspecific genetic diversity. Analyses of chloroplast sequences and nuclear microsatellites indicated high levels of genetic diversity. The diversity distribution was spatially heterogeneous and a latitudinal cline existed in both markers. The spatial distribution pattern between genetic diversity within A. mono and the diversity at species level was generally consistent. Western subtropical evergreen broad-leaved forest subregion (IVb) had a unique ancient chloroplast clade (CP3) and a nuclear gene pool (GP5) with dominance indicating the critical role of this area in species diversification. Genetic data and ecological niche model results both suggested that populations in region II disappeared during the last glacial maximum (LGM) and recovered from south of Changbai Mt. and the Korean Peninsula. Two distribution centers were likely during the LGM, one in the north edge of warm temperate deciduous broad-leaved forest region (III) and another in the south edge of region III. This was reflected by the genetic pattern with two spatially independent genetic groups. This study highlights the key role of region III in sustaining genetic diversity in the northern range and connecting diversity between southern and northern range. We elucidated the diversity relationship between vegetation regions which could facilitate the understanding of biodiversity origin and maintenance in East Asia.
Acer mono; community assemblage; phylogeography; temperate mixed forests
The Patagonian steppe is an immense, cold, arid region, yet phylogeographically understudied. Nassauvia subgen. Strongyloma is a characteristic element of the steppe, exhibiting a continuum of morphological variation. This taxon provides a relevant phylogeographical model not only to understand how past environmental changes shaped the genetic structure of its populations, but also to explore phylogeographical scenarios at the large geographical scale of the Patagonian steppe. Here, we (1) assess demographic processes and historical events that shaped current geographic patterns of haplotypic diversity; (2) analyze hypotheses of isolation in refugia, fragmentation of populations, and/or colonization of available areas during Pleistocene glaciations; and (3) model extant and palaeoclimatic distributions to support inferred phylogeographical patterns. Chloroplast intergenic spacers, rpl32–trnL and trnQ–5′rps16, were sequenced for 372 individuals from 63 populations. Nested clade analysis, analyses of molecular variance, and neutrality tests were performed to assess genetic structure and range expansion. The present potential distribution was modelled and projected onto a last glacial maximum (LGM) model. Of 41 haplotypes observed, ten were shared among populations associated with different morphological variants. Populations with highest haplotype diversity and private haplotypes were found in central-western and south-eastern Patagonia, consistent with long-term persistence in refugia during Pleistocene. Palaeomodelling suggested a shift toward the palaeoseashore during LGM; new available areas over the exposed Atlantic submarine platform were colonized during glaciations with postglacial retraction of populations. A scenario of fragmentation and posterior range expansion may explain the observed patterns in the center of the steppe, which is supported by palaeomodelling. Northern Patagonian populations were isolated from southern populations by the Chubut and the Deseado river basins during glaciations. Pleistocene glaciations indirectly impacted the distribution, demography, and diversification of subgen. Strongyloma through decreased winter temperatures and water availability in different areas of its range.
Colonization; fragmentation; glacial refugia; Patagonian steppe; plant phylogeography; Pleistocene glaciations; range expansion; species distribution modelling
In East Asia, an increasing number of studies on temperate forest tree species find evidence for migration and gene exchange across the East China Sea (ECS) land bridge up until the last glacial maximum (LGM). However, it is less clear when and how lineages diverged in this region, whether in full isolation or in the face of post-divergence gene flow. Here, we investigate the effects of Quaternary changes in climate and sea level on the evolutionary and demographic history of Platycrater arguta, a rare temperate understorey shrub with disjunct distributions in East China (var. sinensis) and South Japan (var. arguta). Molecular data were obtained from 14 P. arguta populations to infer current patterns of molecular structure and diversity in relation to past (Last Interglacial and Last Glacial Maximum) and present distributions based on ecological niche modelling (ENM). A coalescent-based isolation-with-migration (IM) model was used to estimate lineage divergence times and population demographic parameters.
Combining information from nuclear/chloroplast sequence data with nuclear microsatellites, our IM analyses identify the two varieties as genetically distinct units that evolved in strict allopatry since the mid-Pleistocene, c. 0.89 (0.51–1.2) Ma. Together with Bayesian Skyeline Plots, our data further suggest that both lineages experienced post-divergence demographic growth, followed by refugial isolation, divergence, and in the case of var. arguta post-glacial admixture. However, past species distribution modelling indicates that the species’ overall distribution has not greatly changed over the last glacial cycles.
Our findings highlight the important influence of ancient sea-level changes on the diversification of East Asia’s temperate flora. Implicitly, they challenge the notion of general temperate forest expansion across the ECS land bridge, demonstrating instead its ‘filter’ effect owing to an unsuitable environment for certain species and their biological (e.g., recruitment) properties.
Background and Aims
Four species of Pinus subsection Australes occur in the Caribbean Basin: P. caribaea, P. cubensis, P. maestrensis and P. occidentalis. This study analyses the phylogeography of these species to assess possible colonization events from Central America to the islands and subsequent population expansions during glacial periods driven by both drier climate and larger emerged land areas.
Allele size data were obtained for plastid microsatellites for 314 individuals from 24 populations, covering the distribution range of subsection Australes in the Caribbean Basin.
In total, 113 plastid haplotypes were identified. The highest genetic diversity was found in populations of P. caribaea. Overall, Caribbean Basin populations fit the isolation by distance model. Significant phylogeographical structure was found (RST = 0·671 > permuted RST = 0·101; P < 0·0001). The haplotype network and a Bayesian analysis of population structure (BAPS) indicated different Central American origins for P. caribaea var. bahamensis and P. caribaea var. caribaea plastids, with Central America populations in northern and south-eastern groups. Sudden expansion times for BAPS clusters were close to three glacial maxima.
Central America contains ancestral plastid haplotypes. Population expansion has played a major role in the distribution of genetic diversity in P. caribaea var. hondurensis. Two colonization events gave rise to the P. caribaea var. bahamensis and P. caribaea var. caribaea lineages. Plastid variation in the eastern species (P. cubensis, P. maestrensis and P. occidentalis) evolved independently from that in P. caribaea var. caribaea. Incomplete lineage sorting between P. cubensis and P. maestrensis is apparent. Inferred expansion times for P. caribaea var. bahamensis and for the eastern lineages correspond to glacial maxima, whereas those for P. caribaea var. hondurensis correspond to the beginning of the temperature decrease that led to Marine Isotope Stage 8.
Caribbean Basin; historical demography; phylogeography; Pinus; plastid microsatellites; population genetics; subsection Australes
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.
The Qinghai-Tibetan Plateau (QTP) has become one of the hotspots for phylogeographical studies due to its high species diversity. However, most previous studies have focused on the effects of the Quaternary glaciations on phylogeographical structures and the locations of glacial refugia, and little is known about the effects of the aridization of interior Asia on plant population structure and speciation. Here the chloroplast DNA (cpDNA) trnT-trnF and trnS-trnfM sequences were used to investigate the differentiation and phylogeographical history of 14 Ephedra species from the QTP and northern China, based on a sampling of 107 populations. The phylogeographical analysis, together with phylogenetic reconstruction based on combined four cpDNA fragments (rbcL, rpl16, rps4, and trnS-trnfM), supports three main lineages (eastern QTP, southern QTP, and northern China) of these Ephedra species. Divergence of each lineage could be dated to the Middle or Late Miocene, and was very likely linked to the uplift of the QTP and the Asian aridification, given the high drought and/or cold tolerance of Ephedra. Most of the Ephedra species had low intraspecific variation and lacked a strong phylogeographical structure, which could be partially attributed to clonal reproduction and a relatively recent origin. In addition, ten of the detected 25 cpDNA haplotypes are shared among species, suggesting that a wide sampling of species is helpful to investigate the origin of observed haplotypes and make reliable phylogeographical inference. Moreover, the systematic positions of some Ephedra species are discussed.