Bud phenology and growth are subject to divergent selection across a latitudinal gradient in Populus angustifolia and impact adaptation across the distributional range and associated arthropods

Temperate forest tree species that span large geographical areas and climatic gradients often have high levels of genetic variation. Such species are ideal for testing how neutral demographic factors and climate‐driven selection structure genetic variation within species, and how this genetic variation can affect ecological communities. Here, we quantified genetic variation in vegetative phenology and growth traits in narrowleaf cottonwood, Populus angustifolia, using three common gardens planted with genotypes originating from source populations spanning the species' range along the Rocky Mountains of North America (ca. 1700 km). We present three main findings. First, we found strong evidence of divergent selection (Q_ST > F_ST) on fall phenology (bud set) with adaptive consequences for frost avoidance. We also found evidence for selection on bud flush duration, tree height, and basal diameter, resulting in population differentiation. Second, we found strong associations with climate variables that were strongly correlated with latitude of origin. More strongly differentiated traits also showed stronger climate correlations, which emphasizes the role that climate has played in divergent selection throughout the range. We found population × garden interaction effects; for some traits, this accounted for more of the variance than either factor alone. Tree height was influenced by the difference in climate of the source and garden locations and declined with increasing transfer distance. Third, growth traits were correlated with dependent arthropod community diversity metrics. Synthesis. Overall, we conclude that climate has influenced genetic variation and structure in phenology and growth traits and leads to local adaptation in P. angustifolia, which can then impact dependent arthropod species. Importantly, relocation of genotypes far northward or southward often resulted in poor growth, likely due to a phenological mismatch with photoperiod, the proximate cue for fall growth cessation. Genotypes moved too far southward suffer from early growth cessation, whereas those moved too far northward are prone to fall frost and winter dieback. In the face of current and forecasted climate change, habitat restoration, forestry, and tree breeding efforts should utilize these findings to better match latitudinal and climatic source environments with management locations for optimal future outcomes.     

Population at the edge: increased divergence but not inbreeding towards northern range limit in Acer campestre

Ecological conditions shape natural distribution of plants. Populations are denser in optimal habitats but become more fragmented in the areas of suboptimal environmental conditions. Usually, fragmentation increases towards the limits of species distribution. Fragmented populations are often characterised by decreased genetic variation, and this effect is frequent in peripheral populations, mostly due to the reduced effective population size. Interestingly, the genetic consequences of fragmentation seem to be relatively weak in forest trees. Using microsatellite markers, we assessed the impact of population fragmentation on the genetic structure of a European tree species Acer campestre. Within the study area, this medium-size wind-dispersed and insect-pollinated tree reveals a gradual decrease in population density towards the northern range limit. Over the distance of 150 km, we detected the significant decrease in allelic richness, heterozygosity as well as an increase in the rate of population divergence along with latitude. On the other hand, we failed to show that the observed patterns of genetic structure result from the variation in population densities. Moreover, inbreeding levels revealed no association with both density and geographic location, suggesting that pollen limitation does not occur, even at the range margin. As we showed that there is no difference in a dispersal scale between low- and high-density populations in the study species, we argue that the genetic structure is a result of postglacial recolonization. However, unlike many other forest trees, A. campestre showed the sharp latitudinal genetic pattern at a very restricted spatial scale. Limited dispersal and high fragmentation are likely the reasons.     

Influence of late Quaternary climate change on present patterns of genetic variation in valley oak,Quercus lobata Née

Phylogeography and ecological niche models (ENMs) suggest that late Quaternary glacial cycles have played a prominent role in shaping present population genetic structure and diversity, but have not applied quantitative methods to dissect the relative contribution of past and present climate vs. other forces. We integrate multilocus phylogeography, climate‐based ENMs and multivariate statistical approaches to infer the effects of late Quaternary climate change on contemporary genetic variation of valley oak (Quercus lobata Née). ENMs indicated that valley oak maintained a stable distribution with local migration from the last interglacial period (~120 ka) to the Last Glacial Maximum (~21 ka, LGM) to the present compared with large‐scale range shifts for an eastern North American white oak (Quercus alba L.). Coast Range and Sierra Nevada foothill populations diverged in the late Pleistocene before the LGM [104 ka (28–1622)] and have occupied somewhat distinct climate niches, according to ENMs and coalescent analyses of divergence time. In accordance with neutral expectations for stable populations, nuclear microsatellite diversity positively correlated with niche stability from the LGM to present. Most strikingly, nuclear and chloroplast microsatellite variation significantly correlated with LGM climate, even after controlling for associations with geographic location and present climate using partial redundancy analyses. Variance partitioning showed that LGM climate uniquely explains a similar proportion of genetic variance as present climate (16% vs. 11–18%), and together, past and present climate explains more than geography (19%). Climate can influence local expansion–contraction dynamics, flowering phenology and thus gene flow, and/or impose selective pressures. These results highlight the lingering effect of past climate on genetic variation in species with stable distributions.     

Combining ecological niche modelling and morphology to assess the range‐wide population genetic structure of bobcats (Lynx rufus)

Despite a broad distribution, general habitat requirements, and a large dispersal potential, bobcats (Lynx rufus) exhibit a genetic division that longitudinally transects central North America. We investigated (1) whether the climate of the Last Glacial Maximum (LGM; 21 kya) isolated bobcats into refugia and also whether the current climate influences gene flow between the segregate populations and (2) whether the geographical patterns in cranial morphology reflect population identity. We created ecological niche models (ENMs) to evaluate climatic suitability and to estimate distributions of the disparate populations under both historical (LGM) and contemporary conditions. We used two‐dimensional geometric morphometric methods to evaluate variations in the cranium and mandible. These variations were then regressed across geographical variables to assess morphological differences throughout the range of the bobcat. ENMs projected onto LGM climate provided evidence of refugia during the LGM via increased suitability in the north‐west and south‐east portions of this species' range. Contemporarily, our models suggest that the Great Plains may be restricting bobcat migration and gene flow, effectively maintaining disparate populations. Morphological analyses identified a significant linear trend in shape variation across latitudinal and longitudinal gradients rather than distinct morphological divergence between lineages. Similar shape variations, however, did converge in approximate locations of assumed refugia. The findings of the present study provide a robust assessment of the biogeographical considerations for the population genetic structure of bobcats.     

Signatures of diversifying selection at EST‐SSR loci and association with climate in natural Eucalyptus populations

Understanding the environmental parameters that drive adaptation among populations is important in predicting how species may respond to global climatic changes and how gene pools might be managed to conserve adaptive genetic diversity. Here, we used Bayesian F_ST outlier tests and allele–climate association analyses to reveal two Eucalyptus EST‐SSR loci as strong candidates for diversifying selection in natural populations of a southwestern Australian forest tree, Eucalyptus gomphocephala (Myrtaceae). The Eucalyptus homolog of a CONSTANS‐like gene was an F_ST outlier, and allelic variation showed significant latitudinal clinal associations with annual and winter solar radiation, potential evaporation, summer precipitation and aridity. A second F_ST outlier locus, homologous to quinone oxidoreductase, was significantly associated with measures of temperature range, high summer temperature and summer solar radiation, with important implications for predicting the effect of temperature on natural populations in the context of climate change. We complemented these data with investigations into neutral population genetic structure and diversity throughout the species range. This study provides an investigation into selection signatures at gene‐homologous EST‐SSRs in natural Eucalyptus populations, and contributes to our understanding of the relationship between climate and adaptive genetic variation, informing the conservation of both putatively neutral and adaptive components of genetic diversity.     

The population genomic signature of environmental selection in the widespread insect-pollinated tree species Frangula alnus at different geographical scales

The evaluation of the molecular signatures of selection in species lacking an available closely related reference genome remains challenging, yet it may provide valuable fundamental insights into the capacity of populations to respond to environmental cues. We screened 25 native populations of the tree species Frangula alnus subsp. alnus (Rhamnaceae), covering three different geographical scales, for 183 annotated single-nucleotide polymorphisms (SNPs). Standard population genomic outlier screens were combined with individual-based and multivariate landscape genomic approaches to examine the strength of selection relative to neutral processes in shaping genomic variation, and to identify the main environmental agents driving selection. Our results demonstrate a more distinct signature of selection with increasing geographical distance, as indicated by the proportion of SNPs (i) showing exceptional patterns of genetic diversity and differentiation (outliers) and (ii) associated with climate. Both temperature and precipitation have an important role as selective agents in shaping adaptive genomic differentiation in F. alnus subsp. alnus, although their relative importance differed among spatial scales. At the ‘intermediate'' and ‘regional'' scales, where limited genetic clustering and high population diversity were observed, some indications of natural selection may suggest a major role for gene flow in safeguarding adaptability. High genetic diversity at loci under selection in particular, indicated considerable adaptive potential, which may nevertheless be compromised by the combined effects of climate change and habitat fragmentation.     

Testing the role of ancient and contemporary landscapes on structuring genetic variation in a specialist grasshopper

Understanding the processes underlying spatial patterns of genetic diversity and structure of natural populations is a central topic in evolutionary biogeography. In this study, we combine data on ancient and contemporary landscape composition to get a comprehensive view of the factors shaping genetic variation across the populations of the scrub‐legume grasshopper (Chorthippus binotatus binotatus) from the biogeographically complex region of southeast Iberia. First, we examined geographical patterns of genetic structure and employed an approximate Bayesian computation (ABC) approach to compare different plausible scenarios of population divergence. Second, we used a landscape genetic framework to test for the effects of (1) Late Miocene paleogeography, (2) Pleistocene climate fluctuations, and (3) contemporary topographic complexity on the spatial patterns of population genetic differentiation. Genetic structure and ABC analyses supported the presence of three genetic clusters and a sequential west‐to‐east splitting model that predated the last glacial maximum (LGM, c. 21 Kya). Landscape genetic analyses revealed that population genetic differentiation was primarily shaped by contemporary topographic complexity, but was not explained by any paleogeographic scenario or resistance distances based on climate suitability in the present or during the LGM. Overall, this study emphasizes the need of integrating information on ancient and contemporary landscape composition to get a comprehensive view of their relative importance to explain spatial patterns of genetic variation in organisms inhabiting regions with complex biogeographical histories.     

Using genetic diversity and mating system parameters estimated from genetic markers to determine strategies for the conservation of <Emphasis Type="Italic">Araucaria angustifolia</Emphasis> (Bert.) O. Kuntze (Araucariaceae)

In order to understand the impacts of forest fragmentation on Araucaria angustifolia populations, we evaluated the genetic diversity and mating system using SSR markers and open-pollinated seeds from four populations of varying sizes and spatial isolation, in and around one of the best-conserved Araucaria Forest remnants in Southern Brazil. The four population types of A. angustifolia include: (1) a continuous forest; (2) a physically isolated cluster located 2 km from the continuous forest; (3) an open population in a field located between the cluster and continuous forest; and (4) a fragment on a private property located 5 km from the cluster. Approximately 28 seeds were collected from ten reproductive trees in each population. We found higher amounts of alleles (113) and exclusive alleles (25) in the continuous forest than in the other populations. The multilocus paternity correlation was significantly higher and effective number of pollen donors was significantly lower in the private population, decreasing the diversity and consequently the variance effective size of families sampled from that population. However, despite its isolation from the other studied fragments, the private population had the second highest number of alleles as well as unique alleles from the other populations. Therefore, strategies for A. angustifolia conservation should focus not only on larger populations, such as those found in protected areas, but also include smaller and isolated fragments on private properties as these populations are able to maintain high levels of genetic diversity and functional connectivity between isolated stands across a landscape.     

Climatic signals in tree‐rings of Araucaria angustifolia in the southern Brazilian highlands

Araucaria angustifolia (Bertol.) O. Kuntze (Araucariaceae) is a Neotropical tree, widely distributed in subtropical mountain rain forests and nearby natural grasslands of Southern Brazil. This species produces annual growth rings, but its dendroclimatic potential is barely known. In the present paper, the long‐term growth patterns of A. angustifolia were investigated using annual growth ring time series and association to climate over the last century. Wood cores of A. angustifolia trees growing in forest and grassland habitats were obtained with an increment borer. The cores were surfaced, measured and cross‐dated. The dated ring‐width time series were standardized and submitted to correlation and principal component analysis to verify growth trends among sites and trees. Growth‐climate relationships were investigated using correlation and regression analyses, comparing the ordination axes scores to regional time series of precipitation and temperature. Due to anatomical irregularities, mainly partial rings, only 35 out of 60 trees were cross‐dated. The correlation and ordination analyses showed common tree‐growth trends within and between sites, indicative of a regional environmental force determining inter‐annual cambial activity variation. Despite growing in distinct habitats and disturbance regimes, A. angustifolia trees share a common long‐term growth pattern, which is significantly related to thermal conditions during the current and previous growing seasons. Moreover, site‐specific characteristics may have influenced opposite growth responses and association to climate conditions between forest and grassland trees.     

Adaptive population differentiation in phenology across a latitudinal gradient in European aspen (Populus tremula, L.): a comparison of neutral markers, candidate genes and phenotypic traits

A correct timing of growth cessation and dormancy induction represents a critical ecological and evolutionary trade-off between survival and growth in most forest trees (Rehfeldt et al. 1999; Horvath et al. 2003; Howe et al. 2003). We have studied the deciduous tree European Aspen (Populus tremula) across a latitudinal gradient and compared genetic differentiation in phenology traits with molecular markers. Trees from 12 different areas covering 10 latitudinal degrees were cloned and planted in two common gardens. Several phenology traits showed strong genetic differentiation and clinal variation across the latitudinal gradient, with Q(ST) values generally exceeding 0.5. This is in stark contrast to genetic differentiation at several classes of genetic markers (18 neutral SSRs, 7 SSRs located close to phenology candidate genes and 50 SNPs from five phenology candidate genes) that all showed F(ST) values around 0.015. We thus find strong evidence for adaptive divergence in phenology traits across the latitudinal gradient. However, the strong population structure seen at the quantitative traits is not reflected in underlying candidate genes. This result fit theoretical expectations that suggest that genetic differentiation at candidate loci is better described by F(ST) at neutral loci rather than by Q(ST) at the quantitative traits themselves.     

Climate‐driven reduction of genetic variation in plant phenology alters soil communities and nutrient pools

We examined the hypothesis that climate‐driven evolution of plant traits will influence associated soil microbiomes and ecosystem function across the landscape. Using a foundation tree species, Populus angustifolia, observational and common garden approaches, and a base population genetic collection that spans 17 river systems in the western United States, from AZ to MT, we show that (a) as mean annual temperature (MAT) increases, genetic and phenotypic variation for bud break phenology decline; (b) soil microbiomes, soil nitrogen (N), and soil carbon (C) vary in response to MAT and conditioning by trees; and (c) with losses of genetic variation due to warming, population‐level regulation of community and ecosystem functions strengthen. These results demonstrate a relationship between the potential evolutionary response of populations and subsequent shifts in ecosystem function along a large temperature gradient.     

Ecological interpretation of multiple population size structures in trees: The case of Araucaria angustifolia in South America

Abstract In this paper I examined the usefulness of tree population size distributions in evaluating the conservation status of populations of an endangered tree species. I set expectations derived from two complimentary views of the ecology of rainforest trees and examined whether they were met by size distributions of populations of the South American dominant conifer Araucaria angustifolia and its relationship with forest structural characteristics. Specifically, I evaluated the expectations that (i) A. angustifolia trees have larger diameter than average angiosperm trees and form a higher monospecific canopy layer above the shorter angiosperm canopies; (ii) A. angustifolia populations are characterized by size distributions with many large individuals and a long tail of relatively rare, small individuals (have symmetry coefficient <0); (iii) the symmetry of the size distribution of A. angustifolia populations is negatively related to the abundance of large (d.b.h. ≥ 10.0 cm) individuals in the population; and (iv) the abundance of A. angustifolia trees is negatively related to the abundance of angiosperm trees, as the successional accumulation of angiosperm stems would not be accompanied by the recruitment of new A. angustifolia. These expectations were evaluated using data on 25 populations of A. angustifolia sampled in the Rio Grande do Sul State in southern Brazil. The first, third and fourth expectations were met, while the second one was only partially met. This was because populations showed a bimodal distribution regarding symmetry in size distribution, with most populations showing normal or negative symmetry, and those logged, adult‐depleted populations showing positive skewness.     

Phylogeography of the endangered rosewood Dalbergia nigra (Fabaceae): insights into the evolutionary history and conservation of the Brazilian Atlantic Forest

The Brazilian rosewood (Dalbergia nigra) is an endangered tree endemic to the central Brazilian Atlantic Forest, one of the world''s most threatened biomes. The population diversity, phylogeographic structure and demographic history of this species were investigated using the variation in the chloroplast DNA (cpDNA) sequences of 185 individuals from 19 populations along the geographical range of the species. Fifteen haplotypes were detected in the analysis of 1297 bp from two non-coding sequences, trnV-trnM and trnL. We identified a strong genetic structure (F_ST=0.62, P<0.0001), with a latitudinal separation into three phylogeographic groups. The two northernmost groups showed evidence of having maintained historically larger populations than the southernmost group. Estimates of divergence times between these groups pointed to vicariance events in the Middle Pleistocene (ca. 350 000–780 000 years ago). The recurrence of past climatic changes in the central part of the Atlantic forest, with cycles of forest expansion and contraction, may have led to repeated vicariance events, resulting in the genetic differentiation of these groups. Based on comparisons among the populations of large reserves and small, disturbed fragments of the same phylogeographic group, we also found evidence of recent anthropogenic effects on genetic diversity. The results were also analysed with the aim of contributing to the conservation of D. nigra. We suggest that the three phylogeographic groups could be considered as three distinct management units. Based on the genetic diversity and uniqueness of the populations, we also indicate priority areas for conservation.     

Genetic analysis of admixture and patterns of introgression in foundation cottonwood trees (Salicaceae) in southwestern Colorado,USA

Cottonwoods are well known as foundation riparian trees that support diverse communities and drive ecosystem processes. Although hybridization naturally occurs when the distributions of two or more cottonwood species overlap, few cottonwood hybrid zones have been genetically characterized. We use genetic and genomic analyses to characterize patterns of admixture and introgression for a newly described hybrid zone at the intersection of three species (Populus L. Salicaceae—Populus deltoides, Populus fremontii, and Populus angustifolia) in southwestern Colorado, USA. Analysis of nuclear and chloroplast microsatellite marker data detected substantial genetic variation among individuals, revealing that (1) hybridization is occurring between two, not three, species (P. deltoides and P. angustifolia); (2) gene flow is bidirectional; (3) hybrids are not abundant (admixture detected in only 34 of 270 trees), with most being early-generation F₁ hybrids; (4) cytonuclear disequilibria exists and F₁ hybrids tend to retain P. deltoides—like chloroplasts; and (5) roughly 30 % of the nuclear markers deviated from a neutral pattern of introgression, suggesting that selection may play a role in shaping the genetic structure of the hybrid zone in this region. Overall, our results show that despite strong selection maintaining species divergence, transfer of allelic variants across species boundaries can occur. Our study assesses the fine-scale genetic structure of hybridization between P. angustifolia and P. deltoides and lays the foundation for examining how geographic differences in hybrid zone dynamics arise and may influence subsequent ecological and evolutionary processes.     

Pollen,wind and fire: how to investigate genetic effects of disturbance‐induced change in forest trees

Understanding the consequences of habitat disturbance on mating patterns although pollen and seed dispersal in forest trees has been a long‐standing theme of forest and conservation genetics. Forest ecosystems face global environmental pressures from timber exploitation to genetic pollution and climate change, and it is therefore essential to comprehend how disturbances may alter the dispersal of genes and their establishment in tree populations in order to formulate relevant recommendations for sustainable resource management practices and realistic predictions of potential adaptation to climate change by means of range shift or expansion (Kremer et al. 2012 ). However, obtaining reliable evidence of disturbance‐induced effects on gene dispersal processes from empirical evaluation of forest tree populations is difficult. Indeed, tree species share characteristics such as high longevity, long generation time and large reproductive population size, which may impede the experimenter's ability to assess parameters at the spatial and time scales at which any change may occur (Petit and Hampe 2006 ). It has been suggested that appropriate study designs should encompass comparison of populations before and after disturbance as well as account for demonstrated variation in conspecific density, that is, the spatial distribution of mates, and forest density, including all species and relating to alteration in landscape openness (Bacles & Jump 2011 ). However, more often than not, empirical studies aiming to assess the consequences of habitat disturbance on genetic processes in tree populations assume rather than quantify a change in tree densities in forests under disturbance and generally fail to account for population history, which may lead to inappropriate interpretation of a causal relationship between population genetic structure and habitat disturbance due to effects of unmonitored confounding variables (Gauzere et al. 2013). In this issue, Shohami and Nathan ( 2014 ) take advantage of the distinctive features of the fire‐adapted wind‐pollinated Aleppo pine Pinus halepensis (Fig. 1) to provide an elegant example of best practice. Thanks to long‐term monitoring of the study site, a natural stand in Israel, Shohami and Nathan witnessed the direct impact of habitat disturbance, here taking the shape of fire, on conspecific and forest densities and compared pre‐ and postdisturbance mating patterns estimated from cones of different ages sampled on the same surviving maternal individuals (Fig. 2). This excellent study design is all the more strong that Shohami and Nathan took further analytical steps to account for confounding variables, such as historical population genetic structure and possible interannual variation in wind conditions, thus giving high credibility to their findings of unequivocal fire‐induced alteration of mating patterns in P. halepensis. Most notably, the authors found, at the pollen pool level, a disruption of local genetic structure which, furthermore, they were able to attribute explicitly to enhanced pollen‐mediated gene immigration into the low‐density fire‐disturbed stand. This cleverly designed research provides a model approach to be followed if we are to advance our understanding of disturbance‐induced dispersal and genetic change in forest trees.     

Local adaptation shapes pattern of mitochondrial population structure in <Emphasis Type="Italic">Sebastiscus marmoratus</Emphasis>

Considering the homogeneity of marine environment, temperature might be a main environmental factor responsible for local adaptation across a wide latitudinal range. We reconstructed population structure and performed selection tests on a fragment of 1821 bp of the mitochondrial ATPase 6 (ATP6) and cytochrome b (CYTB) genes of Sebastiscus marmoratus sampled from five populations throughout the northwest Pacific Ocean across a large latitudinal range. In line with the results of F _st statistics, significant population structure was detected by applying PCoA and STRUCTURE analysis, which showed discrepant genetic patterns compared to previous studies using neutral markers. Selection tests showed that purifying selection was the predominant evolutionary force acting on the analyzed sequences. Valine substitution by Alanine in a single codon 372 in CYTB gene was detected predominantly (93%) in FA population and temperature was a potentially relevant environmental feature. These findings suggest that local adaptation might be shaping the population structure of marbled rockfish. Further, discrepant patterns of population structure suggested that discerning contributions of different mechanisms in shaping and maintaining patterns of population structure is essential to better understand complex population structure and dynamics of marine organisms.     

Population Genetics of Indigenous Quercus robur L. Populations and of Derived Half‐Sib Families has Implications for the Reproductive Management of the Species

Abstract: In the Netherlands indigenous Quercus robur L. populations are rare and have been maintained as patches in ancient woodland. For adequate conservation of these populations, information about genetic variation and population structure is necessary. In order to assess the genetic variation and structure of these populations, microsatellite polymorphisms were studied in two autochthonous populations. These two populations differed slightly for their gene diversity, which was as high as was observed for Q. robur populations in France and Germany. For reforestation purposes there is an interest in the genetic variation of a half‐sib family harvested from one tree. The gene diversity of the two studied half‐sib families ‐ obtained from a forest and an urban area ‐ was similar, but relatively low. This indicates that, for reforestation purposes, seeds should be harvested from many different trees in order to obtain a population with a genetic variation as high as was observed for an autochthonous population.     

Latitudinal Environmental Niches and Riverine Barriers Shaped the Phylogeography of the Central Chilean Endemic Dioscorea humilis (Dioscoreaceae)

The effects of Pleistocene glaciations and geographical barriers on the phylogeographic patterns of lowland plant species in Mediterranean-climate areas of Central Chile are poorly understood. We used Dioscorea humilis (Dioscoreaceae), a dioecious geophyte extending 530 km from the Valparaíso to the Bío-Bío Regions, as a case study to disentangle the spatio-temporal evolution of populations in conjunction with latitudinal environmental changes since the Last Inter-Glacial (LIG) to the present. We used nuclear microsatellite loci, chloroplast (cpDNA) sequences and environmental niche modelling (ENM) to construct current and past scenarios from bioclimatic and geographical variables and to infer the evolutionary history of the taxa. We found strong genetic differentiation at nuclear microsatellite loci between the two subspecies of D. humilis, probably predating the LIG. Bayesian analyses of population structure revealed strong genetic differentiation of the widespread D. humilis subsp. humilis into northern and southern population groups, separated by the Maipo river. ENM revealed that the ecological niche differentiation of both groups have been maintained up to present times although their respective geographical distributions apparently fluctuated in concert with the climatic oscillations of the Last Glacial Maximum (LGM) and the Holocene. Genetic data revealed signatures of eastern and western postglacial expansion of the northern populations from the central Chilean depression, whereas the southern ones experienced a rapid southward expansion after the LGM. This study describes the complex evolutionary histories of lowland Mediterranean Chilean plants mediated by the summed effects of spatial isolation caused by riverine geographical barriers and the climatic changes of the Quaternary.     

Contrasting genetic patterns between two coexisting Eleutherococcus species in northern China

Climate oscillations are the key factors to understand the patterns in modern biodiversity. East Asia harbors the most diverse temperate flora, largely because an extensive terrestrial ice cap was absent during repeated Pleistocene glaciation–interglacial cycles. Comparing the demographic histories of species that are codistributed and are close relatives may provide insight into how the process of climate change influences species ranges. In this study, we compared the spatial genetic structure and demographic histories of two coexisting Eleutherococcus species, Eleutherococcus senticosus and E. sessiliflorus. Both species are distributed in northern China, regions that are generally considered to be sensitive to climatic fluctuations. These regions once hosted temperate forest, but this temperate forest was replaced by tundra and taiga forest during the Last Glacial Maximum (LGM), according to pollen records. Using three chloroplast DNA fragments, we assessed the genetic structure of 20 and 9 natural populations of E. senticosus and E. sessiliflorus, respectively. Extremely contrasting genetic patterns were found between the two species; E. sessiliflorus had little genetic variation, whereas E. senticosus had considerably higher levels of genetic variation (15 haplotypes). We speculated that a recent severe bottleneck may have resulted in the extremely low genetic diversity in E. sessiliflorus. In E. senticosus, populations in Northeast China (NEC) harbored all of the haplotypes found in this species and included private haplotypes. The populations in NEC had higher levels of genetic diversity than did those from North China (NC). Therefore, we suggest that both the NC and NEC regions can sustain LGM refugia and that lineage admixture from multiple refugia took place after the LGM elevated the local genetic diversity in NEC. In NEC, multiple genetic hot spots were found in the Changbai Mountains and the Xiaoxing'an Range, which implied that multiple locations in NEC may sustain LGM refugia, even in the Xiaoxing'an Range.