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.

Plants adapt to the environment by either long-term genome evolution or by acclimatization processes where the cellular processes and metabolism of the plant are adjusted within the existing potential in the genome. Here we studied the adaptation strategies in date palm, Phoenix dactylifera, under mild heat, drought and combined heat and drought by transcriptomic and metabolomic profiling. In transcriptomics data, combined heat and drought resembled heat response, whereas in metabolomics data it was more similar to drought. In both conditions, soluble carbohydrates, such as fucose, and glucose derivatives, were increased, suggesting a switch to carbohydrate metabolism and cell wall biogenesis. This result is consistent with the evidence from transcriptomics and cis-motif analysis. In addition, transcriptomics data showed transcriptional activation of genes related to reactive oxygen species in all three conditions (drought, heat, and combined heat and drought), suggesting increased activity of enzymatic antioxidant systems in cytosol, chloroplast and peroxisome. Finally, the genes that were differentially expressed in heat and combined heat and drought stresses were significantly enriched for circadian and diurnal rhythm motifs, suggesting new stress avoidance strategies.

The atmospheric pollutant ozone (O ) is a strong oxidant that causes extracellular reactive oxygen species (ROS) formation, has significant ecological relevance, and is used here as a non-invasive ROS inducer to study plant signalling. Previous genetic screens identified several mutants exhibiting enhanced O sensitivity, but few with enhanced tolerance. We found that loss-of-function mutants in Arabidopsis MLO2, a gene implicated in susceptibility to powdery mildew disease, exhibit enhanced dose-dependent tolerance to O and extracellular ROS, but a normal response to intracellular ROS. This phenotype is increased in a mlo2 mlo6 mlo12 triple mutant, reminiscent of the genetic redundancy of MLO genes in powdery mildew resistance. Stomatal assays revealed that enhanced O tolerance in mlo2 mutants is not caused by altered stomatal conductance. We explored modulation of the mlo2-associated O tolerance, powdery mildew resistance, and early senescence phenotypes by genetic epistasis analysis, involving mutants with known effects on ROS sensitivity or antifungal defence. Mining of publicly accessible microarray data suggests that these MLO proteins regulate accumulation of abiotic stress response transcripts, and transcript accumulation of MLO2 itself is O responsive. In summary, our data reveal MLO2 as a novel negative regulator in plant ROS responses, which links biotic and abiotic stress response pathways.

The photoreceptors UV RESISTANCE LOCUS 8 (UVR8) and CRYPTOCHROMES 1 and 2 (CRYs) play major roles in the perception of UV-B (280–315 nm) and UV-A/blue radiation (315–500 nm), respectively. However, it is poorly understood how they function in sunlight. The roles of UVR8 and CRYs were assessed in a factorial experiment with Arabidopsis thaliana wild-type and photoreceptor mutants exposed to sunlight for 6 h or 12 h under five types of filters with cut-offs in UV and blue-light regions. Transcriptome-wide responses triggered by UV-B and UV-A wavelengths shorter than 350 nm (UV-Asw) required UVR8 whereas those induced by blue and UV-A wavelengths longer than 350 nm (UV-Alw) required CRYs. UVR8 modulated gene expression in response to blue light while lack of CRYs drastically enhanced gene expression in response to UV-B and UV-Asw. These results agree with our estimates of photons absorbed by these photoreceptors in sunlight and with in vitro monomerization of UVR8 by wavelengths up to 335 nm. Motif enrichment analysis predicted complex signaling downstream of UVR8 and CRYs. Our results highlight that it is important to use UV waveband definitions specific to plants’ photomorphogenesis as is routinely done in the visible region.