A long-term genetic legacy of refugial isolation has been postulated and was demonstrated for maternal refugial lineages for numerous plant and animal species. The lineages were assumed to have remained separated from each other for several glacial periods. The conifer Abies alba Miller, silver fir, is an excellent model to test whether pollen-mediated gene flow may eliminate the genetic imprints of Pleistocene refugial isolation. Two DNA markers with contrasting modes of inheritance were applied to 100 populations covering the entire range of silver fir in Europe. The markers exhibited each two highly conserved alleles based on an insertion͞deletion of 80 bp in the fourth intron of the mitochondrial nad5 gene and on a synonymous substitution in the chloroplast psbC gene. The geographical distribution of the maternally inherited mitochondrial variation supported the existence of at least two refugia with two recolonizing maternal lineages remaining largely separated throughout the range. The cline of the nad5 allele frequencies was much steeper than the one of the two psbC alleles. The psbC cline was as wide as the whole range of the species. Our results provide striking evidence that even a species with very long generation times and heavy pollen grains was able to establish a highly efficient pollen-mediated gene flow between refugia. Therefore we postulate that an exchange of genetic information between refugia by range-wide paternal introgression is possible in wind-pollinated plant species.

Because of heterogeneous topographies, high-mountain areas could harbor a significant pool of cryptic forest refugia (glacial microrefugia unrecognized by palaeodata), which, as a result of poor accessibility, have been largely overlooked. The juniper forests of the southern Tibetan Plateau, with one of the highest tree lines worldwide, are ideal for assessing the potential of high-mountain areas to harbor glacial refugia. Genetic evidence for Last Glacial Maximum (LGM) endurance of these microrefugia is presented using paternally inherited chloroplast markers. Five-hundred and ninety individuals from 102 populations of the Juniperus tibetica complex were sequenced at three polymorphic chloroplast regions. Significant interpopulation differentiation and phylogeographic structure were detected (G(ST) = 0.49, N(ST) = 0.72, N(ST) > G(ST), P < 0.01), indicating limited among-population gene flow. Of 62 haplotypes recovered, 40 were restricted to single populations. These private haplotypes and overall degrees of diversity were evenly spread among plateau and edge populations, strongly supporting the existence of LGM microrefugia throughout the present distribution range, partly well above 3500 m. These results mark the highest LGM tree lines known, illustrating the potential significance of high-mountain areas for glacial refugia. Furthermore, as the close vicinity of orographic rear-edge and leading-edge populations potentially allows gene flow, surviving populations could preserve the complete spectrum of rear-edge and leading-edge adaptations.

No abstract

Understanding local adaptation in forest trees is currently a key research and societal priority. Geographically and ecologically marginal populations provide ideal case studies, because environmental stress along with reduced gene flow can facilitate the establishment of locally adapted populations. We sampled European silver fir (Abies alba Mill.) trees in the French Mediterranean Alps, along the margin of its distribution range, from pairs of high- and low-elevation plots on four different mountains situated along a 170-km east-west transect. The analysis of 267 SNP loci from 175 candidate genes suggested a neutral pattern of east-west isolation by distance among mountain sites. F(ST) outlier tests revealed 16 SNPs that showed patterns of divergent selection. Plot climate was characterized using both in situ measurements and gridded data that revealed marked differences between and within mountains with different trends depending on the season. Association between allelic frequencies and bioclimatic variables revealed eight genes that contained candidate SNPs, of which two were also detected using F(ST) outlier methods. All SNPs were associated with winter drought, and one of them showed strong evidence of selection with respect to elevation. Q(ST)-F(ST) tests for fitness-related traits measured in a common garden suggested adaptive divergence for the date of bud flush and for growth rate. Overall, our results suggest a complex adaptive picture for A. alba in the southern French Alps where, during the east-to-west Holocene recolonization, locally advantageous genetic variants established at both the landscape and local scales.

We present a novel measurement setup for monitoring changes in leaf water status using nondestructive terahertz time-domain spectroscopy (THz-TDS). Previous studies on a variety of plants showed the principal applicability of THz-TDS. In such setups, decreasing leaf water content directly correlates with increasing THz transmission. Our new system allows for continuous, nondestructive monitoring of the water status of multiple individual plants each at the same constant leaf position. It overcomes previous drawbacks, which were mainly due to the necessity of relocating the plants. Using needles of silver fir (Abies alba) seedlings as test subjects, we show that the transmission varies along the main axis of a single needle due to a variation in thickness. Therefore, the relocation of plants during the measuring period, which was necessary in the previous THz-TDS setups, should be avoided. Furthermore, we show a highly significant correlation between gravimetric water content and respective THz transmission. By monitoring the relative change in transmission, we were able to narrow down the permanent wilting point of the seedlings. Thus, we established groups of plants with well-defined levels of water stress that could not be detected visually. This opens up the possibility for a broad range of genetic and physiological experiments.

Genetic association studies in forest trees would greatly benefit from information on the response of trees to environmental stressors over time, which can be provided by dendroecological analysis. Here, we jointly analysed dendroecological and genetic data of surviving silver fir trees to explore the genetic basis of their response to the iconic stress episode of the 1970s and 1980s that led to large-scale forest dieback in Central Europe and has been attributed to air pollution. Specifically, we derived dendrophenotypic measures from 190 trees in the Bavarian Forest that characterize the resistance, resilience and recovery during this growth depression, and in the drought year in 1976. By focusing on relative growth changes of trees and by standardizing the dendrophenotypes within stands, we accounted for variation introduced by micro- and macroscale environmental differences. We associated the dendrophenotypes with single nucleotide polymorphisms (SNPs) in candidate genes using general linear models (GLMs) and the machine learning algorithm random forest with subsequent feature selection. Most trees at our study sites experienced a severe growth decline from 1974 until the mid-1980s with minimum values during the drought year. Fifteen genes were associated with the dendrophenotypes, including genes linked to photosynthesis and drought stress. With our study, we show that dendrophenotypes can be a powerful resource for genetic association studies that permit to account for micro- and macroenvironmental variation when data are derived from natural populations. We call for a wider collaboration of dendroecologists and forest geneticists to integrate individual tree-level dendrophenotypes in genetic association studies.

No abstract

Silver fir ( Abies alba Mill.) is a keystone conifer of European montane forest ecosystems that has experienced large fluctuations in population size during during the Quaternary and, more recently, due to land-use change. To forecast the species’ future distribution and survival, it is important to investigate the genetic basis of adaptation to environmental change, notably to extreme events. For this purpose, we here provide a first draft genome assembly and annotation of the silver fir genome, established through a community-based initiative. DNA obtained from haploid megagametophyte and diploid needle tissue was used to construct and sequence Illumina paired-end and mate-pair libraries, respectively, to high depth. The assembled A. alba genome sequence accounted for over 37 million scaffolds corresponding to 18.16 Gb, with a scaffold N50 of 14,051 bp. Despite the fragmented nature of the assembly, a total of 50,757 full-length genes were functionally annotated in the nuclear genome. The chloroplast genome was also assembled into a single scaffold (120,908 bp) that shows a high collinearity with both the A. koreana and A. sibirica complete chloroplast genomes. This first genome assembly of silver fir is an important genomic resource that is now publicly available in support of a new generation of research. By genome-enabling this important conifer, this resource will open the gate for new research and more precise genetic monitoring of European silver fir forests.