A Reference Genome Sequence for the European Silver Fir (Abies alba Mill.): A Community-Generated Genomic Resource.

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.

Detection of somatic epigenetic variation in Norway spruce via targeted bisulfite sequencing.

Epigenetic mechanisms represent a possible mechanism for achieving a rapid response of long-lived trees to changing environmental conditions. However, our knowledge on plant epigenetics is largely limited to a few model species. With increasing availability of genomic resources for many tree species, it is now possible to adopt approaches from model species that permit to obtain single-base pair resolution data on methylation at a reasonable cost. Here, we used targeted bisulfite sequencing (TBS) to study methylation patterns in the conifer species Norway spruce (Picea abies). To circumvent the challenge of disentangling epigenetic and genetic differences, we focused on four clone pairs, where clone members were growing in different climatic conditions for 24 years. We targeted >26.000 genes using TBS and determined the performance and reproducibility of this approach. We characterized gene body methylation and compared methylation patterns between environments. We found highly comparable capture efficiency and coverage across libraries. Methylation levels were relatively constant across gene bodies, with 21.3 ± 0.3%, 11.0 ± 0.4% and 1.3 ± 0.2% in the CG, CHG, and CHH context, respectively. The variance in methylation profiles did not reveal consistent changes between environments, yet we could identify 334 differentially methylated positions (DMPs) between environments. This supports that changes in methylation patterns are a possible pathway for a plant to respond to environmental change. After this successful application of TBS in Norway spruce, we are confident that this approach can contribute to broaden our knowledge of methylation patterns in natural tree populations.

Local Adaptation in European Firs Assessed through Extensive Sampling across Altitudinal Gradients in Southern Europe.

BACKGROUND: Local adaptation is a key driver of phenotypic and genetic divergence at loci responsible for adaptive traits variations in forest tree populations. Its experimental assessment requires rigorous sampling strategies such as those involving population pairs replicated across broad spatial scales. METHODS: A hierarchical Bayesian model of selection (HBM) that explicitly considers both the replication of the environmental contrast and the hierarchical genetic structure among replicated study sites is introduced. Its power was assessed through simulations and compared to classical 'within-site' approaches (FDIST, BAYESCAN) and a simplified, within-site, version of the model introduced here (SBM). RESULTS: HBM demonstrates that hierarchical approaches are very powerful to detect replicated patterns of adaptive divergence with low false-discovery (FDR) and false-non-discovery (FNR) rates compared to the analysis of different sites separately through within-site approaches. The hypothesis of local adaptation to altitude was further addressed by analyzing replicated Abies alba population pairs (low and high elevations) across the species' southern distribution range, where the effects of climatic selection are expected to be the strongest. For comparison, a single population pair from the closely related species A. cephalonica was also analyzed. The hierarchical model did not detect any pattern of adaptive divergence to altitude replicated in the different study sites. Instead, idiosyncratic patterns of local adaptation among sites were detected by within-site approaches. CONCLUSION: Hierarchical approaches may miss idiosyncratic patterns of adaptation among sites, and we strongly recommend the use of both hierarchical (multi-site) and classical (within-site) approaches when addressing the question of adaptation across broad spatial scales.

Evidence of divergent selection for drought and cold tolerance at landscape and local scales in Abies alba Mill. in the French Mediterranean Alps.

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.

Within-range translocations and their consequences in European larch.

In contrast to biological invasions, translocations of individuals within a species range are understudied, due to difficulties in systematically detecting them. This results in limited knowledge about the corresponding processes and uncertainties regarding the status of extant populations. European larch, a forest tree whose fragmented native distribution is restricted to the Alps and to other Central European mountains, has been massively planted for at least 300 years. Here we focus on the genetic characterization of translocations having taken place within its native range. Microsatellite variation at 13 nuclear loci and sequence data of two mitochondrial DNA fragments were analyzed on the basis of a comprehensive range-wide population sample. Two complementary methods (Geneclass and Structure) were used to infer translocation events based on nuclear data whereas mitochondrial data were used for validation of these inferences. Using Geneclass, we found translocation events in a majority of populations. Additional cases of translocation and many instances of admixture were identified using Structure, thanks to the clear-cut ancestral genetic structure detected in this species. In particular, a strong divide between Alpine and Central European populations, also apparent at mitochondrial markers, helped uncover details on translocation events and related processes. Translocations and associated admixture events were found to be heterogeneously distributed across the species range, with a particularly high frequency in Central Europe. Furthermore, translocations frequently involved multiple geographic sources, some of which were over-represented. Our study illustrates the importance of range-wide investigations for tracing translocations back to their origins and for revealing some of their consequences. It provides some first clues for developing suitable conservation and management strategies.

Differential Gene Expression Reveals Candidate Genes for Drought Stress Response in Abies alba (Pinaceae).

Increasing drought periods as a result of global climate change pose a threat to many tree species by possibly outpacing their adaptive capabilities. Revealing the genetic basis of drought stress response is therefore implemental for future conservation strategies and risk assessment. Access to informative genomic regions is however challenging, especially for conifers, partially due to their large genomes, which puts constraints on the feasibility of whole genome scans. Candidate genes offer a valuable tool to reduce the complexity of the analysis and the amount of sequencing work and costs. For this study we combined an improved drought stress phenotyping of needles via a novel terahertz water monitoring technique with Massive Analysis of cDNA Ends to identify candidate genes for drought stress response in European silver fir (Abies alba Mill.). A pooled cDNA library was constructed from the cotyledons of six drought stressed and six well-watered silver fir seedlings, respectively. Differential expression analyses of these libraries revealed 296 candidate genes for drought stress response in silver fir (247 up- and 49 down-regulated) of which a subset was validated by RT-qPCR of the twelve individual cotyledons. A majority of these genes code for currently uncharacterized proteins and hint on new genomic resources to be explored in conifers. Furthermore, we could show that some traditional reference genes from model plant species (GAPDH and eIF4A2) are not suitable for differential analysis and we propose a new reference gene, TPC1, for drought stress expression profiling in needles of conifer seedlings.

Characterization of 23 polymorphic SSR markers in Salix humboldtiana (Salicaceae) using next-generation sequencing and cross-amplification from related species.

PREMISE OF THE STUDY: We present a set of 23 polymorphic nuclear microsatellite loci, 18 of which are identified for the first time within the riparian species Salix humboldtiana (Salicaceae) using next-generation sequencing. METHODS AND RESULTS: To characterize the 23 loci, up to 60 individuals were sampled and genotyped at each locus. The number of alleles ranged from two to eight, with an average of 4.43 alleles per locus. The effective number of alleles ranged from 1.15 to 3.09 per locus, and allelic richness ranged from 2.00 to 7.73 alleles per locus. CONCLUSIONS: The new marker set will be used for future studies of genetic diversity and differentiation as well as for unraveling spatial genetic structures in S. humboldtiana populations in northern Patagonia, Argentina.

Monitoring plant drought stress response using terahertz time-domain spectroscopy.

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.

Annotation and re-sequencing of genes from de novo transcriptome assembly of Abies alba (Pinaceae).

PREMISE OF THE STUDY: We present a protocol for the annotation of transcriptome sequence data and the identification of candidate genes therein using the example of the nonmodel conifer Abies alba. • METHODS AND RESULTS: A normalized cDNA library was built from an A. alba seedling. The sequencing on a 454 platform yielded more than 1.5 million reads that were de novo assembled into 25149 contigs. Two complementary approaches were applied to annotate gene fragments that code for (1) well-known proteins and (2) proteins that are potentially adaptively relevant. Primer development and testing yielded 88 amplicons that could successfully be resequenced from genomic DNA. • CONCLUSIONS: The annotation workflow offers an efficient way to identify potential adaptively relevant genes from the large quantity of transcriptome sequence data. The primer set presented should be prioritized for single-nucleotide polymorphism detection in adaptively relevant genes in A. alba.

Identification and characterization of 12 microsatellite loci in Cnidium dubium (Apiaceae) using next-generation sequencing.

PREMISE OF THE STUDY: We present a set of 12 polymorphic nuclear microsatellite loci, identified for the first time within the herbaceous plant species Cnidium dubium (Schkuhr) Thell. using next-generation sequencing. METHODS AND RESULTS: To characterize these loci, 40 plants were analyzed genetically. The number of alleles ranged from 2 to 23, with an average of 8.3 alleles per locus. Mean observed and expected heterozygosities ranged from 0.25 to 0.95 and 0.30 to 0.94, respectively. CONCLUSIONS: All 12 nuclear microsatellite loci showed scorable and polymorphic fragments after polymerase chain reaction. The new marker set can be used for future studies of genetic diversity and differentiation as well as estimation of gene flow and spatial genetic structures.

Tree endurance on the Tibetan Plateau marks the world's highest known tree line of the Last Glacial Maximum.

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.

Authenticated DNA from ancient wood remains.

BACKGROUND: The reconstruction of biological processes and human activities during the last glacial cycle relies mainly on data from biological remains. Highly abundant tissues, such as wood, are candidates for a genetic analysis of past populations. While well-authenticated DNA has now been recovered from various fossil remains, the final 'proof' is still missing for wood, despite some promising studies. SCOPE: The goal of this study was to determine if ancient wood can be analysed routinely in studies of archaeology and palaeogenetics. An experiment was designed which included blind testing, independent replicates, extensive contamination controls and rigorous statistical tests. Ten samples of ancient wood from major European forest tree genera were analysed with plastid DNA markers. CONCLUSIONS: Authentic DNA was retrieved from wood samples up to 1,000 years of age. A new tool for real-time vegetation history and archaeology is ready to use.

Wind-dispersed pollen mediates postglacial gene flow among refugia.

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.