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.

Large protein families are a prominent feature of plant genomes and their size variation is a key element for adaptation. However, gene and genome duplications pose difficulties for functional characterization and translational research. Here we infer the evolutionary history of the DOMAIN OF UNKNOWN FUNCTION (DUF) 26-containing proteins. The DUF26 emerged in secreted proteins. Domain duplications and rearrangements led to the appearance of CYSTEINE-RICH RECEPTOR-LIKE PROTEIN KINASES (CRKs) and PLASMODESMATA-LOCALIZED PROTEINS (PDLPs). The DUF26 is land plant-specific but structural analyses of PDLP ectodomains revealed strong similarity to fungal lectins and thus may constitute a group of plant carbohydrate-binding proteins. CRKs expanded through tandem duplications and preferential retention of duplicates following whole genome duplications, whereas PDLPs evolved according to the dosage balance hypothesis. We propose that new gene families mainly expand through small-scale duplications, while fractionation and genetic drift after whole genome multiplications drive families towards dosage balance.

We describe the genome contents of six Protomyces spp. that are pathogenic within the typical host range of the genus and a novel Protomyces strain (SC29) that was previously isolated from the phylloplane of wild Arabidopsis thaliana (Arabidopsis), an atypical or possible alternate host. Genome-wide phylogenetic analysis defined SC29 as a distinct Protomyces sp. Analysis of gene family expansions, gene retention, and gene loss patterns among these Protomyces spp. lead us to hypothesize that SC29 may have undergone a host jump. The role of phyllosphere residency in the lifecycle of Protomyces spp. was previously unknown. Genomic changes in SC29 and all other Protomyces spp. were consistent with adaptations to the plant phylloplane. As predicted by our analysis of its mating locus, SC29 did not cause disease on Arabidopsis as a single strain, but could persist in its phylloplane, while the closely related P. inouyei does not. SC29 treated Arabidopsis exhibited enhanced immunity against Botrytis cinerea infection, associated with activation of MAPK3/6, camalexin, and SA-signalling pathways. We conclude that SC29 is a novel Protomyces sp. able to survive in the Arabidopsis phylloplane and that phylloplane residency is an important element in the lifecycle of Protomyces spp.

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.

The following Supporting Information is available for this article: Figure S1. Solar spectrum at different times of the day when plants were moved outdoors. Figure S2. Photon irradiance for different wavebands in solar radiation.Figure S3. Multidimensional scaling of RNA-seq data.Figure S4. Comparison between RNA-seq and qRT-PCR data.Figure S5. Venn diagrams showing the number differentially expressed genes in RNA-seq data.Figure S6. Enrichment of KEGG pathways in RNA-seq data.Figure S7. In vitro absorption spectra of Arabidopsis UVR8 protein.Figure S8. Position weight matrices of the enriched DNA-binding motifs.Figure S9. Transcript abundance of seven genes measured using qRT-PCR.Table S1. Information of primers used and genes assessed in qRT-PCR. Table S2. Summary of the ANOVA of the qRT-PCR data.Methods S1. Description of the filters and the waveband contrasts.Dataset S1. Outcome of differential gene expression analysis for the three genotypes and multiple waveband contrasts combination. The dataset is included as a separate file in .Rda format and can be read using R.