Forest ecosystems maintain a large share of Northern Hemisphere biodiversity. Boreal forests comprise roughly 30% of global forest area 1 and contain the highest tree density among climate zones 2 . Moreover, boreal regions are undergoing extensive climate change. Annual temperature increases exceeding 1.5 °C are projected to result in a warming of 4-11 °C by the end of this century, with little concomitant increase in precipitation 1 . At this pace, climate zones will shift northward at a greater speed than trees can migrate 3 . To understand how future populations of forest trees may respond to climate change, it is essential to uncover past and present signatures of molecular adaptation in their genomes. Silver birch, B. pendula, is a pioneer species in boreal forests of Eurasia. Flowering of the species can be artificially accelerated 4 , giving it an advantage over other tree species with published genome sequences (such as poplar 5 , spruce 6 , and pine 7 ) for the optimization of fiber and biomass production.Here we sequenced 150 birch individuals and assembled a B. pendula reference genome from a fourth-generation inbred line, resulting in a high-quality assembly of 435 Mb that was linked to chromosomes using a dense genetic map. We analyzed SNPs in the genomes of 80 birch individuals spanning most of the geographic range of B. pendula, as well as seven other members of Betulaceae. Population genomic analyses of these data provide insights into the deep-time evolution of the birch family and on recent natural selection acting on silver birch.Genome sequencing and population genomic analyses provide insights into the adaptive landscape of silver birch Silver birch (Betula pendula) is a pioneer boreal tree that can be induced to flower within 1 year. Its rapid life cycle, small (440-Mb) genome, and advanced germplasm resources make birch an attractive model for forest biotechnology. We assembled and chromosomally anchored the nuclear genome of an inbred B. pendula individual. Gene duplicates from the paleohexaploid event were enriched for transcriptional regulation, whereas tandem duplicates were overrepresented by environmental responses. Population resequencing of 80 individuals showed effective population size crashes at major points of climatic upheaval. Selective sweeps were enriched among polyploid duplicates encoding key developmental and physiological triggering functions, suggesting that local adaptation has tuned the timing of and cross-talk between fundamental plant processes. Variation around the tightlylinked light response genes PHYC and FRS10 correlated with latitude and longitude and temperature, and with precipitation for PHYC. Similar associations characterized the growth-promoting cytokinin response regulator ARR1, and the wood development genes KAK and MED5A.A full list of affiliations appears at the end of the paper.

Populations of Noccaea caerulescens show tremendous differences in their capacity to hyperaccumulate and hypertolerate metals. To explore the differences that could contribute to these traits, we undertook SOLiD high-throughput sequencing of the root transcriptomes of three phenotypically well-characterized N. caerulescens accessions, i.e., Ganges, La Calamine, and Monte Prinzera. Genes with possible contribution to zinc, cadmium, and nickel hyperaccumulation and hypertolerance were predicted. The most significant differences between the accessions were related to metal ion (di-, trivalent inorganic cation) transmembrane transporter activity, iron and calcium ion binding, (inorganic) anion transmembrane transporter activity, and antioxidant activity. Analysis of correlation between the expression profile of each gene and the metal-related characteristics of the accessions disclosed both previously characterized (HMA4, HMA3) and new candidate genes (e.g., for nickel IRT1, ZIP10, and PDF2.3) as possible contributors to the hyperaccumulation/tolerance phenotype. A number of unknown Noccaea-specific transcripts also showed correlation with Zn(2+), Cd(2+), or Ni(2+) hyperaccumulation/tolerance. This study shows that N. caerulescens populations have evolved great diversity in the expression of metal-related genes, facilitating adaptation to various metalliferous soils. The information will be helpful in the development of improved plants for metal phytoremediation.

Several populations with diVerent metal tolerance, uptake and root-to-shoot transport are known for the metal hyperaccumulator plant Thlaspi caerulescens. In this study, genes diVerentially expressed under various Zn exposures were identiWed from the shoots of two T. caerulescens accessions (calaminous and non-calaminous) using Xuorescent diVerential display RT-PCR. cDNA fragments from 16 Zn-responsive genes, including those encoding metallothionein (MT) type 2 and type 3, MRP-like transporter, pectin methylesterase (PME) and Ole e 1-like gene as well as several unknown genes, were eventually isolated. The fulllength MT2 and MT3 sequences diVer from those previously isolated from other Thlaspi accessions, possibly representing new alleles or isoforms. Besides the diVerential expression in Zn exposures, the gene expression was dependent on the accession. Thlaspi homologues of ClpP protease and MRP transporter were induced at high Zn concentrations. MT2 and PME were expressed at higher levels in the calaminous accession. The MTs and MRP transporter expressed in transgenic yeasts were capable of conferring Cu and Cd tolerance, whereas the Ole e 1-like gene enhanced toxicity to these metals. The MTs increased yeast intracellular Cd content. As no signiWcant diVerences were found between Arabidopsis and Thlaspi MTs, they apparently do not diVer in their capacity to bind metals. However, the higher levels of MT2 in the calaminous accession may contribute to the Zn-adapted phenotype.

The metal hyperaccumulator Noccaea caerulescens is an established model to study the adaptation of plants to metalliferous soils. Various comparators have been used in these studies. The choice of suitable comparators is important and depends on the hypothesis to be tested and methods to be used. In high-throughput analyses such as microarray, N. caerulescens has been compared to non-tolerant, non-accumulator plants like Arabidopsis thaliana or Thlaspi arvense rather than to the related hypertolerant or hyperaccumulator plants. An underutilized source is N. caerulescens populations with considerable variation in their capacity to accumulate and tolerate metals. Whole transcriptome sequencing (RNA-Seq) is revealing interesting variation in their gene expression profiles. Combining physiological characteristics of N. caerulescens accessions with their RNA-Seq has a great potential to provide detailed insight into the underlying molecular mechanisms, including entirely new gene products. In this review we will critically consider comparative transcriptome analyses carried out to explore metal hyperaccumulation and hypertolerance of N. caerulescens, and demonstrate the potential of RNA-Seq analysis as a tool in evolutionary genomics.

These authors contributed equally to this work. SUMMARYCalamine accessions of the zinc/cadmium/nickel hyperaccumulator, Noccaea caerulescens, exhibit striking variation in foliar cadmium accumulation in nature. The Ganges accession (GA) from Southern France displays foliar cadmium hyperaccumulation (>1000 lg g À1 DW), whereas the accession La Calamine (LC) from Belgium, with similar local soil metal composition, does not (<100 lg g À1 DW). All calamine accessions are cadmium hypertolerant. To find out the differences between LC and GA in their basic adaptation mechanisms, we bypassed the cadmium excluding phenotype of LC by exposing the plants to 50 lM cadmium in hydroponics, achieving equal cadmium accumulation in the shoots. The iron content increased in the roots of both accessions. GA exhibited significant decreases in manganese and zinc contents in the roots and shoots, approaching those in LC. Altogether 702 genes responded differently to cadmium exposure between the accessions, 157 and 545 in the roots and shoots, respectively. Cadmium-exposed LC showed a stress response and had decreased levels of a wide range of photosynthesis-related transcripts. GA showed less changes, mainly exhibiting an iron deficiency-like response. This included increased expression of genes encoding five iron deficiency-regulated bHLH transcription factors, ferric reduction oxidase FRO2, iron transporters IRT1 and OPT3, and nicotianamine synthase NAS1, and decreased expression of genes encoding ferritins and NEET (a NEET family iron-sulfur protein), which is possibly involved in iron transfer, distribution and/or management. The function of the IRT1 gene in the accessions was compared. We conclude that the major difference between the two accessions is in the way they cope with iron under cadmium exposure.