Guillaume de Lafontaine

Climate refugia: joint inference from fossil records, species distribution models and phylogeography

Climate refugia, locations where taxa survive periods of regionally adverse climate, are thought to be critical for maintaining biodiversity through the glacial-interglacial climate changes of the Quaternary. A critical research need is to better integrate and reconcile the three major lines of evidence used to infer the existence of past refugia - fossil records, species distribution models and phylogeographic surveys - in order to characterize the complex spatiotemporal trajectories of species and populations in and out of refugia. Here we review the complementary strengths, limitations and new advances for these three approaches. We provide case studies to illustrate their combined application, and point the way towards new opportunities for synthesizing these disparate lines of evidence. Case studies with European beech, Qinghai spruce and Douglas-fir illustrate how the combination of these three approaches successfully resolves complex species histories not attainable from any one approach. Promising new statistical techniques can capitalize on the strengths of each method and provide a robust quantitative reconstruction of species history. Studying past refugia can help identify contemporary refugia and clarify their conservation significance, in particular by elucidating the fine-scale processes and the particular geographic locations that buffer species against rapidly changing climate.

Stronger spatial genetic structure in recolonized areas than in refugia in the European beech

Extant rear‐edge populations located in former glacial refugia remain understudied despite their high conservation value. These populations should have experienced strong genetic drift due to their small size and long isolation. Moreover, the prolonged action of isolation by distance in refugial areas should result in stronger regional spatial genetic structure (SGS) than in recolonized areas, but empirical tests of this prediction are scarce. To fill this gap, we first used a set of 16 microsatellite markers to investigate the genetic structure of European beech in France in 65 populations from three refugial areas and one control recolonized (nonrefugial) area. Then, using the same approach, we reanalysed published isozyme data from 375 populations distributed across the entire species range. We found stronger genetic differentiation among populations in refugia than in recolonized areas. However, contrary to expectations, regional SGS was lower within refugia than within recolonized areas. Published studies presenting similar analyses suggest that our results could have generality across different biogeographical settings and types of organisms. Strong and prolonged genetic drift in refugial areas could have erased the signature of range expansions that is still visible in recolonized areas. Our results therefore suggest that Pleistocene population isolation has played a key role in increasing the genetic complexity of extant rear‐edge populations.

Cryptic no more: soil macrofossils uncover Pleistocene forest microrefugia within a periglacial desert.

Despite their critical importance for understanding the local effects of global climate change on biodiversity, glacial microrefugia are not well studied because they are difficult to detect by using classical palaeoecological or population genetics approaches. We used soil macrofossil charcoal analysis to uncover the presence of cryptic glacial refugia for European beech (Fagus sylvatica) and other tree species in the Landes de Gascogne (southwestern France). Using botanical identification and direct radiocarbon dating (140 (14) C-dates) of macrofossil charcoal extracted from mineral soils, we reconstructed the glacial and postglacial history of all extant beech stands in the region (n = 11). Soil charcoal macrofossils were found in all sites, allowing the identification of up to at least 14 distinct fire events per site. There was direct evidence of the presence of beech during the last glacial period at three sites. Beech was detected during Heinrich stadial-1, one of the coldest and driest intervals of the last glacial period in Western Europe. Together with previous results on the genetic structure of the species in the region, these findings suggest that beech persisted in situ in several microrefugia through full glacial and interglacial periods up to the present day.

Long-term fire and forest history of subalpine balsam fir (Abies balsamea) and white spruce (Picea glauca) stands in eastern Canada inferred from soil charcoal analysis

The northernmost balsam fir forest in eastern Canada forms disjunct stands far beyond the extensive balsam fir forest zone of southern Canada. The northern balsam fir stands are distributed in the subalpine belt of high plateaus and coexist locally with white spruce stands. These subalpine stands contrast greatly with black spruce forest stands located in lowlands. Given that subalpine stands are remnants of an earlier northern expansion of the balsam fir forest, the main objective of this study is to assess whether white spruce stands are distinct communities having diverged from the balsam fir forest community earlier in the Holocene or if they rather correspond to a different stage of the chronosequence within the subalpine belt. Macrofossil analysis of charcoal in mineral soils was used to compare the stand-scale fire histories and taxonomic fossil composition of subalpine, old-growth balsam fir stands and white spruce stands. No significant differences of mean number of observed fires (mean = 6.35 fires per site), Holocene fire recurrence at the landscape scale and mean fire-return interval (mean = 580 years) were found between white spruce stands and balsam fir stands. The botanical composition of charcoal fragments from mineral soils showed that Abies, Betula and Picea were present throughout the fire period from 5600 cal. BP to present, and no difference was found in the fossil composition of the balsam fir and white spruce stands. No historical change in the botanical composition of charcoal from soils of both stand types was observed indicating that the initial floristic composition remained through the period of recurrent fires. Charcoal data suggest that white spruce stands are not divergent community types. Rather, the two community types are arranged along a chronosequence of different successional stages within the subalpine relict flora.

The Origin and Dynamics of Subalpine White Spruce and Balsam Fir Stands in Boreal Eastern North America

Associations among the few tree species in the North American boreal landscape are the result of complex interactions between climate, biota, and historical disturbances during the Holocene. The closed-crown boreal forest of eastern North America is subdivided into two ecological regions having distinct tree species associations; the balsam fir zone and the black spruce zone, south and north of 49°N, respectively. Subalpine old-growth stands dominated by trees species typical of the balsam fir forest flora (either balsam fir or white spruce) are found on high plateaus, some of which are isolated within the black spruce zone. Here we identified the ecological processes responsible for the distinct forest associations in the subalpine belt across the eastern boreal landscape. Extensive radiocarbon dating, species composition, and size structure analyses indicated contrasted origin and dynamics of the subalpine forests between the two ecological regions. In the black spruce zone, the subalpine belt is a mosaic of post-fire white spruce or balsam fir stands coexisting at similar elevation on the high plateaus. With increasing time without wildfire, the subalpine forests become structurally similar to the balsam fir forest of the fir zone. These results concur with the hypothesis that the subalpine forests of this area are protected remnants of an historical northern expansion of the fir zone. Its replacement by the fire-prone black spruce forest flora was caused by recurrent fires. In the subalpine belt of the fir zone, no fire was recorded for several millennia. Harsh climate at high altitude is the primary factor explaining white spruce dominance over balsam fir forming a distinct subalpine white spruce belt above the balsam fir dominated forest.

Integrating phylogeography and paleoecology to investigate the origin and dynamics of hybrid zones: insights from two widespread North American firs

Secondary contact between closely related taxa routinely occurs during postglacial migrations. After initial contact, the location of hybrid zones may shift geographically or remain spatially stable over time in response to various selective pressures or neutral processes. Studying the extent and direction of introgression using markers having contrasted levels of gene flow can help unravel the historical dynamics of hybrid zones. Thanks to their contrasted maternal and paternal inheritance, resulting in different levels of gene flow for mitochondrial and chloroplast DNA (mtDNA and cpDNA), the Pinaceae stand out as a relevant biological model for this purpose. The objective of the study was to assess whether the hybrid zone between Abies balsamea and Abies lasiocarpa (two largely distributed Pinaceae) has moved or remained stable over time by analysing the distribution of cytoplasmic DNA variation as well as published palaeobotanical data. Interspecific gene flow was higher for cpDNA than mtDNA markers; hence, the geographic distribution of mitotypes was more congruent with species distributions than chlorotypes. This genetic signature was contrary to expectations under a moving hybrid zone scenario, as well as empirical observations in other conifers. Genetic evidence for this rare instance of stable hybrid zone was corroborated by the colonization chronology derived from published fossil data, indicating that the two fir species initially came into contact in the area corresponding to the current sympatric zone 11 kyr ago. While an explanatory analysis suggested the putative influence of various environmental factors on the relative abundance of cytoplasmic genome combinations, further research appears necessary to assess the role of both demographic history and selective factors in driving the dynamics of hybrid zones.

Tracking the progression of speciation: variable patterns of introgression across the genome provide insights on the species delimitation between progenitor–derivative spruces (Picea mariana × P. rubens)

The genic species concept implies that while most of the genome can be exchanged somewhat freely between species through introgression, some genomic regions remain impermeable to interspecific gene flow. Hence, interspecific differences can be maintained despite ongoing gene exchange within contact zones. This study assessed the heterogeneous patterns of introgression at gene loci across the hybrid zone of an incipient progenitor–derivative species pair, Picea mariana (black spruce) and Picea rubens (red spruce). The spruce taxa likely diverged in geographic isolation during the Pleistocene and came into secondary contact during late Holocene. A total of 300 SNPs distributed across the 12 linkage groups (LG) of black spruce were genotyped for 385 individual trees from 33 populations distributed across the allopatric zone of each species and within the zone of sympatry. An integrative framework combining three population genomic approaches was used to scan the genomes, revealing heterogeneous patterns of introgression. A total of 23 SNPs scattered over 10 LG were considered impermeable to introgression and putatively under diverging selection. These loci revealed the existence of impermeable genomic regions forming the species boundary and are thus indicative of ongoing speciation between these two genetic lineages. Another 238 SNPs reflected selectively neutral diffusion across the porous species barrier. Finally, 39 highly permeable SNPs suggested ancestral polymorphism along with balancing selection. The heterogeneous patterns of introgression across the genome indicated that the speciation process between black spruce and red spruce is young and incomplete, albeit some interspecific differences are maintained, allowing ongoing species divergence even in sympatry. The approach developed in this study can be used to track the progression of ongoing speciation processes.

Species richness along a production gradient: a multivariate approach.

The specific shape of the relationship between plant diversity and productivity and the causal mechanism(s) behind the observed pattern(s) are still highly debated. Recent advances suggest that the relationship depends on several environmental variables and may change with the observational scale. In this study, a multivariate, multiscale approach was used to identify the variables that determine the relationship between species richness and annual production along a forest/old field edge in southern Québec (Canada). Various relationships between richness and production were found at different distances to the edge. In the forest, most relationships were positive and linear, while in the old field the relationship shifted from positive linear to non-significant with increasing distance from the edge. In the forest or in the old field, the shape of the relationship (all distances from the edge combined) was unimodal. Path analyses showed that species richness was determined mostly by production, which was influenced by different limiting resources, depending on the community (forest or old field). An increasing range in production created by pooling across community types can confound the resources and/or conditions determining the diversity-productivity relationship.

The critical role of local refugia in postglacial colonization of Chinese pine: joint inferences from DNA analyses, pollen records, and species distribution modeling

The importance of long-distance migration from low to high latitudes relative to local spread from northern refugia after the Last Glacial Maximum (LGM) remains a focus of debate for many temperate tree species. We assessed the dynamics of Chinese pine Pinus tabulaeformis, a widespread species endemic to northern China, since the LGM by integrating cytoplasmic DNA data, mapped pollen records and ecological niche modeling. Genetic variation among 544 individuals from 50 populations spanning the entire natural species range revealed eight genetic clusters with distinct geographic distribution, indicating glacial lineages likely originating from multiple local microrefugia. Palynological evidence suggested that the northernmost part of the natural distribution originated from local postglacial spread. Niche modeling indicated high probability of the species being present in the area of the Loess Plateau and coastal areas north of the Yangtze River during the LGM. The three lines of evidence jointly suggest that the species persisted through the last glaciation in the mountains surrounding the Loess Plateau of northern China and that the current distribution of the species originated primarily from the spread of local refugial populations, instead of long-distance migration. These results cast doubt on the notion that Chinese pine migrated from areas south of the Yangtze River and underscore the importance of northern refugia.

Asymmetry matters: A genomic assessment of directional biases in gene flow between hybridizing spruces

Abstract Assessing directional bias in interspecific gene flow might be important in determining the evolutionary trajectory of closely related species pairs. Using a set of 300 single nucleotide polymorphisms (SNPs) having variable propensity to cross species boundary, we evaluated the genomic extent and direction of interspecific gene flow in a progenitor‐derivative spruce species pair (black spruce and red spruce). A higher rate of gene flow was found from black spruce toward red spruce purebreds than vice versa. This asymmetry could reflect the historical gene flow between the two taxa at the time of species inception and during postglacial colonization. A clear asymmetry in introgression was depicted by a greater gene flow between red spruce and hybrids than between black spruce and hybrids. While backcrossing toward red spruce was invariably high across the genome, the actual species boundary is between hybrids and black spruce where gene flow is impeded at those genomic regions impermeable to introgression. Associations between hybrid index and climatic variables (total annual precipitation and mean annual temperature) were tested, as these might indicate a role for exogenous selection in maintaining the species boundary. While an apparent association was found between the hybrid index and precipitation, it collapsed when considered in light of the directional bias in interspecific gene flow. Hence, considering asymmetrical patterns of introgression allowed us to falsify an apparent role for exogenous selection. Although this was not formerly tested here, we suggest that this pattern could result from asymmetrical endogenous selection, a contention that deserves further investigations.