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.

Scots pine (Pinus sylvestris L.) wood is desired in woodworking industries due to its favorable timber characteristics and natural durability that is contributed by heartwood extractives. It has been discussed whether the Scots pine heartwood extractives (mainly stilbenes and resin acids) are synthesized in the cells of the transition zone between sapwood and heartwood, or if they are transported from the sapwood. Timing of heartwood formation during the yearly cycle has also not been unambiguously defined. We measured steady-state mRNA levels in Scots pine transition zone and sapwood using RNA sequencing. Year-round expression profiles of selected transcripts were further investigated by quantitative RT-PCR. Differentially accumulating transcripts suggest that, of the Scots pine heartwood extractives, stilbenes are synthesized in situ in the transition zone and gain their carbon-skeletons from Suc and triglycerides. Resin acids, on the other hand, are synthesized early in the spring mainly in the sapwood, meaning that they must be transported to the heartwood transition zone. Heartwood formation is marked by programmed cell death that occurs during the summer months in the transition zone.

Summary Tree bark is a highly specialized array of tissues that plays important roles in plant protection and development. Bark tissues develop from two lateral meristems; the phellogen (cork cambium) produces the outermost stem–environment barrier called the periderm, while the vascular cambium contributes with phloem tissues. Although bark is diverse in terms of tissues, functions and species, it remains understudied at higher resolution. We dissected the stem of silver birch (Betula pendula) into eight major tissue types, and characterized these by a combined transcriptomics and metabolomics approach. We further analyzed the varying bark types within the Betulaceae family. The two meristems had a distinct contribution to the stem transcriptomic landscape. Furthermore, inter‐ and intraspecies analyses illustrated the unique molecular profile of the phellem. We identified multiple tissue‐specific metabolic pathways, such as the mevalonate/betulin biosynthesis pathway, that displayed differential evolution within the Betulaceae. A detailed analysis of suberin and betulin biosynthesis pathways identified a set of underlying regulators and highlighted the important role of local, small‐scale gene duplication events in the evolution of metabolic pathways. This work reveals the transcriptome and metabolic diversity among bark tissues and provides insights to its development and evolution, as well as its biotechnological applications.

Apoplastic events such as monolignol oxidation and lignin polymerization are difficult to study in intact trees. To investigate the role of apoplastic hydrogen peroxide (H 2 O 2 ) in gymnosperm phenolic metabolism, an extracellular lignin-forming cell culture of Norway spruce (Picea abies) was used as a research model. Scavenging of apoplastic H 2 O 2 by potassium iodide repressed lignin formation, in line with peroxidases activating monolignols for lignin polymerization. Time-course analyses coupled to candidate substrate-product pair network propagation revealed differential accumulation of low-molecular-weight phenolics, including (glycosylated) oligolignols, (glycosylated) flavonoids, and proanthocyanidins, in lignin-forming and H 2 O 2 -scavenging cultures and supported that monolignols are oxidatively coupled not only in the cell wall but also in the cytoplasm, where they are coupled to other monolignols and proanthocyanidins. Dilignol glycoconjugates with reduced structures were found in the culture medium, suggesting that cells are able to transport glycosylated dilignols to the apoplast. Transcriptomic analyses revealed that scavenging of apoplastic H 2 O 2 resulted in remodulation of the transcriptome, with reduced carbon flux into the shikimate pathway propagating down to monolignol biosynthesis. Aggregated coexpression network analysis identified candidate enzymes and transcription factors for monolignol oxidation and apoplastic H 2 O 2 production in addition to potential H 2 O 2 receptors. The results presented indicate that the redox state of the apoplast has a profound influence on cellular metabolism.