Genomic basis for drought resistance in European beech forests threatened by climate change

Research on Plant Landscape Design of Urban Industrial Site Green Space Based on Green Infrastructure Concept

With the acceleration of the global urbanization process, more and more industrial plants are being abandoned, which puts great pressure on urban ecology and land resource management. These abandoned industrial spaces not only lead to persistent pollution problems, but also exacerbate the urban heat island effect, leading to a worsening microclimate. To address these issues, the concept of green infrastructure (GI) has emerged as a sustainable ecological restoration strategy, and it is an important tool for urban renewal and industrial land transformation. In this study, the landscape environment of the industrial site of Henrichshutte in Germany and the surrounding industrial plant was taken as an example, and ecological restoration and plant landscape design were carried out using the GI concept. Two climate simulation tools, ENVI-met and WindPerfect DX, were comprehensively adopted to simulate the environment of the site in detail. Based on an analysis of the potential temperature, PMV, wind speed, and UTCI data of the site, it was demonstrated that the plant landscape improved the microclimate of the industrial plant. The results show that the reasonable allocation of plants can effectively reduce surface temperature and building temperature, increase air humidity, alleviate the local heat island effect, and enhance the thermal comfort of the human body. The simulation results highlight the practical application value of the GI concept in improving the ecological benefit, social function, and landscape aesthetics of industrial land. This study provides a new idea for the ecological restoration and environmental optimization of urban industrial land through the combination of green infrastructure and plant landscape design, and emphasizes the important role of green infrastructure in alleviating the urban heat island effect and promoting the sustainable development of urban landscape spaces.

The reference genome of an endangered Asteraceae, Deinandra increscens subsp. villosa, endemic to the Central Coast of California

We present a reference genome for the federally endangered Gaviota tarplant, Deinandra increscens subsp. villosa (Madiinae, Asteraceae), an annual herb endemic to the Central California coast. Generating PacBio HiFi, Oxford Nanopore Technologies, and Dovetail Omni-C data, we assembled a haploid consensus genome of 1.67 Gb as 28.7 K scaffolds with a scaffold N50 of 74.9 Mb. We annotated repeat content in 74.8% of the genome. Long terminal repeats (LTRs) covered 44.0% of the genome with Copia families predominant at 22.9% followed by Gypsy at 14.2%. Both Gypsy and Copia elements were common in ancestral peaks of LTRs, and the most abundant element was a Gypsy element containing nested Copia/Angela sequence similarity, reflecting a complex evolutionary history of repeat activity. Gene annotation produced 33,257 genes and 68,942 transcripts, of which 99% were functionally annotated. BUSCO scores for the annotated proteins were 96.0% complete of which 77.6% was single copy and 18.4% duplicates. Whole genome duplication synonymous mutation rates of Gaviota tarplant and sunflower (Helianthus annuus) shared peaks that correspond to the last Asteraceae polyploidization event and subsequent divergence from a common ancestor at ∼27 MYA. Regions of high-density tandem genes were identified, pointing to potentially important loci of environmental adaptation in this species.

Rolling down that mountain: microgeographical adaptive divergence during a fast population expansion along a steep environmental gradient in European beech

Forest tree populations harbour high genetic diversity thanks to large effective population sizes and strong gene flow, allowing them to diversify through adaptation to local environmental pressures within dispersal distance. Many tree populations also experienced historical demographic fluctuations, including spatial population contraction or expansions at various temporal scales, which may constrain their ability to adapt to environmental variations. Our aim is to investigate how recent contraction and expansion events interfere with local adaptation, by studying patterns of adaptive divergence between closely related stands undergoing environmentally contrasted conditions, and having or not recently expanded. To investigate genome-wide signatures of local adaptation while accounting for demography, we analysed divergence in a European beech population by testing pairwise differentiation among four tree stands at ~35k Single Nucleotide Polymorphisms from ~9k genomic regions. We applied three divergence outlier search methods resting on different assumptions and targeting either single SNPs or contiguous genomic regions, while accounting for the effect of population size variations on genetic divergence. We found 27 signals of selective signatures in 19 target regions. Putatively adaptive divergence involved all stand pairs. We retrieved signals both when comparing old-growth stands and recently colonised areas and when comparing stands within the old-growth area. Therefore, adaptive divergence processes have taken place both over short time spans, under strong environmental contrasts, and over short ecological gradients, in populations that have been stable in the long term. This suggests that standing genetic variation supports local, microgeographic divergence processes, which can maintain genetic diversity at the landscape level.

Genetic response of a perennial grass to warm and wet environments interacts and is associated with trait means as well as plasticity

The potential for rapid evolution is an important mechanism allowing species to adapt to changing climatic conditions. Although such potential has been largely studied in various short-lived organisms, to what extent we can observe similar patterns in long-lived plant species, which often dominate natural systems, is largely unexplored. We explored the potential for rapid evolution in Festuca rubra, a long-lived grass with extensive clonal growth dominating in alpine grasslands. We used a field sowing experiment simulating expected climate change in our model region. Specifically, we exposed seeds from five independent seed sources to novel climatic conditions by shifting them along a natural climatic grid and explored the genetic profiles of established seedlings after 3 years. Data on genetic profiles of plants selected under different novel conditions indicate that different climate shifts select significantly different pools of genotypes from common seed pools. Increasing soil moisture was more important than increasing temperature or the interaction of the two climatic factors in selecting pressure. This can indicate negative genetic interaction in response to the combined effects or that the effects of different climates are interactive rather than additive. The selected alleles were found in genomic regions, likely affecting the function of specific genes or their expression. Many of these were also linked to morphological traits (mainly to trait plasticity), suggesting these changes may have a consequence on plant performance. Overall, these data indicate that even long-lived plant species may experience strong selection by climate, and their populations thus have the potential to rapidly adapt to these novel conditions.

A genome-wide genetic association study reveals SNPs significantly associated with environmental variables and specific leaf area in European beech

European beech is negatively affected by climate change and a further growth decline is predicted for large parts of its distribution range. Despite the importance of this species, little is known about its genetic adaptation and especially the genetic basis of its physiological traits. Here, we used genotyping by sequencing to identify SNPs in 43 German European beech populations growing under different environmental conditions. In total, 28 of these populations were located along a precipitation and temperature gradient in northern Germany, and single tree-based hydraulic and morphological traits were available. We obtained a set of 13,493 high-quality SNPs that were used for environmental and SNP-trait association analysis. In total, 22 SNPs were identified that were significantly associated with environmental variables or specific leaf area (SLA). Several SNPs were located in genes related to stress response. The majority of the significant SNPs were located in non-coding (intergenic and intronic) regions. These may be in linkage disequilibrium with the causative coding or regulatory regions. Our study gives insights into the genetic basis of abiotic adaptation in European beech, and provides genetic resources that can be used in future studies on this species. Besides clear patterns of local adaptation to environmental conditions of the investigated populations, the analyzed morphological and hydraulic traits explained most of the explainable genetic variation. Thus, they could successfully be altered in tree breeding programs, which may help to increase the adaptation of European beech to changing environmental conditions in the future.

No effect of snow on shrub xylem traits: Insights from a snow-manipulation experiment on Disko Island, Greenland

Widespread shrubification across the Arctic has been generally attributed to increasing air temperatures, but responses vary across species and sites. Wood structures related to the plant hydraulic architecture may respond to local environmental conditions and potentially impact shrub growth, but these relationships remain understudied. Using methods of dendroanatomy, we analysed shrub ring width (RW) and xylem anatomical traits of 80 individuals of Salix glauca L. and Betula nana L. at a snow manipulation experiment in Western Greenland. We assessed how their responses differed between treatments (increased versus ambient snow depth) and soil moisture regimes (wet and dry). Despite an increase in snow depth due to snow fences (28-39 %), neither RW nor anatomical traits in either species showed significant responses to this increase. In contrast, irrespective of the snow treatment, the xylem specific hydraulic conductivity (Ks) and earlywood vessel size (LA95) for the study period were larger in S. glauca (p < 0.1, p < 0.01) and B. nana (p < 0.01, p < 0.001) at the wet than the dry site, while both species had larger vessel groups at the dry than the wet site (p < 0.01). RW of B. nana was higher at the wet site (p < 0.01), but no differences were observed for S. glauca. Additionally, B. nana Ks and LA95 showed different trends over the study period, with decreases observed at the dry site (p < 0.001), while for other responses no difference was observed. Our results indicate that, taking into account ontogenetic and allometric trends, hydraulic related xylem traits of both species, along with B. nana growth, were influenced by soil moisture. These findings suggest that soil moisture regime, but not snow cover, may determine xylem responses to future climate change and thus add to the heterogeneity of Arctic shrub dynamics, though more long-term species- and site- specific studies are needed.

First whole-genome sequence and assembly of the Ecuadorian brown-headed spider monkey (Ateles fusciceps fusciceps), a critically endangered species, using Oxford Nanopore Technologies

The Ecuadorian brown-headed spider monkey (Ateles fusciceps fusciceps) is currently considered one of the most endangered primates in the world and is classified as critically endangered [International union for conservation of nature (IUCN)]. It faces multiple threats, the most significant one being habitat loss due to deforestation in western Ecuador. Genomic tools are keys for the management of endangered species, but this requires a reference genome, which until now was unavailable for A. f. fusciceps. The present study reports the first whole-genome sequence and assembly of A. f. fusciceps generated using Oxford Nanopore long reads. DNA was extracted from a subadult male, and libraries were prepared for sequencing following the Ligation Sequencing Kit SQK-LSK112 workflow. Sequencing was performed using a MinION Mk1C sequencer. The sequencing reads were processed to generate a genome assembly. Two different assemblers were used to obtain draft genomes using raw reads, of which the Flye assembly was found to be superior. The final assembly has a total length of 2.63 Gb and contains 3,861 contigs, with an N50 of 7,560,531 bp. The assembly was analyzed for annotation completeness based on primate ortholog prediction using a high-resolution database, and was found to be 84.3% complete, with a low number of duplicated genes indicating a precise assembly. The annotation of the assembly predicted 31,417 protein-coding genes, comparable with other mammal assemblies. A reference genome for this critically endangered species will allow researchers to gain insight into the genetics of its populations and thus aid conservation and management efforts of this vulnerable species.

Slower growth prior to the 2018 drought and a high growth sensitivity to previous year summer conditions predisposed European beech to crown dieback

The record-breaking drought in 2018 caused premature leaf discoloration and shedding (early browning) in many beech (Fagus sylvatica L.) dominated forests in Central Europe. However, a high degree of variability in drought response among individual beech trees was observed. While some trees were severely impacted by the prolonged water deficits and high temperatures, others remained vital with no or only minor signs of crown vitality loss. Why some beech trees were more susceptible to drought-induced crown damage than others and whether growth recovery is possible are poorly understood. Here, we aimed to identify growth characteristics associated with the variability in drought response between individual beech trees based on a sample of 470 trees in northern Switzerland. By combining tree growth measurements and crown condition assessments, we also investigated the possible link between crown dieback and growth recovery after drought. Beech trees with early browning exhibited an overall lower growth vigor before the 2018 drought than co-occurring vital beech trees. This lower vigor is mainly indicated by lower overall growth rates, stronger growth declines in the past decades, and higher growth-climate sensitivity. Particularly, warm previous year summer conditions negatively affected current growth of the early-browning trees. These findings suggest that the affected trees had less access to critical resources and were physiologically limited in their growth predisposing them to early browning. Following the 2018 drought, observed growth recovery potential corresponded to the amount of crown dieback and the local climatic water balance. Overall, our findings emphasize that beech-dominated forests in Central Europe are under increasing pressure from severe droughts, ultimately reducing the competitive ability of this species, especially on lowland sites with shallow soils and low water holding capacity.

A simulation study comparing advanced marker-assisted selection with genomic selection in tree breeding programs

Advances in DNA sequencing technologies allow the sequencing of whole genomes of thousands of individuals and provide several million single nucleotide polymorphisms (SNPs) per individual. These data combined with precise and high-throughput phenotyping enable genome-wide association studies (GWAS) and the identification of SNPs underlying traits with complex genetic architectures. The identified causal SNPs and estimated allelic effects could then be used for advanced marker-assisted selection (MAS) in breeding programs. But could such MAS compete with the broadly used genomic selection (GS)? This question is of particular interest for the lengthy tree breeding strategies. Here, with our new software "SNPscan breeder," we simulated a simple tree breeding program and compared the impact of different selection criteria on genetic gain and inbreeding. Further, we assessed different genetic architectures and different levels of kinship among individuals of the breeding population. Interestingly, apart from progeny testing, GS using gBLUP performed best under almost all simulated scenarios. MAS based on GWAS results outperformed GS only if the allelic effects were estimated in large populations (ca. 10,000 individuals) of unrelated individuals. Notably, GWAS using 3,000 extreme phenotypes performed as good as the use of 10,000 phenotypes. GS increased inbreeding and thus reduced genetic diversity more strongly compared to progeny testing and GWAS-based selection. We discuss the practical implications for tree breeding programs. In conclusion, our analyses further support the potential of GS for forest tree breeding and improvement, although MAS may gain relevance with decreasing sequencing costs in the future.

First genome-wide data from Italian European beech (Fagus sylvatica L.): Strong and ancient differentiation between Alps and Apennines

The European beech (Fagus sylvatica L.) is one of the most widespread forest trees in Europe whose distribution and intraspecific diversity has been largely shaped by repeated glacial cycles. Previous studies, mainly based on palaeobotanical evidence and a limited set of chloroplast and nuclear genetic markers, highlighted a complex phylogeographic scenario, with southern and western Europe characterized by a rather heterogeneous genetic structure, as a result of recolonization from different glacial refugia. Despite its ecological and economic importance, the genome of this broad-leaved tree has only recently been assembled, and its intra-species genomic diversity is still largely unexplored. Here, we performed whole-genome resequencing of nine Italian beech individuals sampled from two stands located in the Alpine and Apennine mountain ranges. We investigated patterns of genetic diversity at chloroplast, mitochondrial and nuclear genomes and we used chloroplast genomes to reconstruct a temporally-resolved phylogeny. Results allowed us to test European beech differentiation on a whole-genome level and to accurately date their divergence time. Our results showed comparable, relatively high levels of genomic diversity in the two populations and highlighted a clear differentiation at chloroplast, mitochondrial and nuclear genomes. The molecular clock analysis indicated an ancient split between the Alpine and Apennine populations, occurred between the Günz and the Riss glaciations (approximately 660 kyrs ago), suggesting a long history of separation for the two gene pools. This information has important conservation implications in the context of adaptation to ongoing climate changes.

How genomics can help biodiversity conservation

The availability of public genomic resources can greatly assist biodiversity assessment, conservation, and restoration efforts by providing evidence for scientifically informed management decisions. Here we survey the main approaches and applications in biodiversity and conservation genomics, considering practical factors, such as cost, time, prerequisite skills, and current shortcomings of applications. Most approaches perform best in combination with reference genomes from the target species or closely related species. We review case studies to illustrate how reference genomes can facilitate biodiversity research and conservation across the tree of life. We conclude that the time is ripe to view reference genomes as fundamental resources and to integrate their use as a best practice in conservation genomics.

Gene abundance linked to climate zone: Parallel evolution of gene content along elevation gradients in lichenized fungi

Introduction

Intraspecific genomic variability affects a species' adaptive potential toward climatic conditions. Variation in gene content across populations and environments may point at genomic adaptations to specific environments. The lichen symbiosis, a stable association of fungal and photobiont partners, offers an excellent system to study environmentally driven gene content variation. Many of these species have remarkable environmental tolerances, and often form populations across different climate zones. Here, we combine comparative and population genomics to assess the presence and absence of genes in high and low elevation genomes of two lichenized fungi of the genus Umbilicaria.

Methods

The two species have non-overlapping ranges, but occupy similar climatic niches in North America (U. phaea) and Europe (U. pustulata): high elevation populations are located in the cold temperate zone and low elevation populations in the Mediterranean zone. We assessed gene content variation along replicated elevation gradients in each of the two species, based on a total of 2050 individuals across 26 populations. Specifically, we assessed shared orthologs across species within the same climate zone, and tracked, which genes increase or decrease in abundance within populations along elevation.

Results

In total, we found 16 orthogroups with shared orthologous genes in genomes at low elevation and 13 at high elevation. Coverage analysis revealed one ortholog that is exclusive to genomes at low elevation. Conserved domain search revealed domains common to the protein kinase superfamily. We traced the discovered ortholog in populations along five replicated elevation gradients on both continents and found that the number of this protein kinase gene linearly declined in abundance with increasing elevation, and was absent in the highest populations.

Discussion

We consider the parallel loss of an ortholog in two species and in two geographic settings a rare find, and a step forward in understanding the genomic underpinnings of climatic tolerances in lichenized fungi. In addition, the tracking of gene content variation provides a widely applicable framework for retrieving biogeographical determinants of gene presence/absence patterns. Our work provides insights into gene content variation of lichenized fungi in relation to climatic gradients, suggesting a new research direction with implications for understanding evolutionary trajectories of complex symbioses in relation to climatic change.

Whole genome sequencing and RNA-seq evaluation allowed to detect Cd adaptation footprint in Chironomus riparius

Evolutionary adaptation and phenotypic plasticity are important processes on how organisms respond to pollutant exposure. We dissected here the contribution of both processes to increased tolerance in Chironomus riparius to cadmium (Cd) exposure in a multi-generation experiment and inferred the underlying genomic basis. We simulated environmentally realistic conditions by continuously increasing contaminant concentration in six replicates initiated with 1000 larvae each, three pre-exposed to Cd and three not exposed to Cd (no-Cd) over eight generations. We measured life-cycle traits, transcriptomic responses and genome-wide allele frequency changes from this evolve and resequencing (E&R) experiment. Overall, life cycle tests revealed little phenotypic adaptation to Cd exposure, but a slightly increase in survival in the first larval stage was observed. Population genomic analyses showed a strong genome-wide selective response in all replicates, highlighting two main biological functions involved in development and growth of the chironomids. Emphasizing that laboratory conditions continually exert selective pressure. However, the integration of the transcriptomic to the genomic data allowed to distinguish pathways specifically selected by the Cd exposure related to microtubules and organelles and cellular movement. Those pathways could be functionally related to an excretion of metals. Thus, our results indicate that genetic adaptation to Cd in C. riparius can happen within few generations under an environmentally relevant exposure scenario, but substantial phenotypic tolerance might take more time to arise. With our approach, we introduce an experimental setup to fill the existing gap in evolutionary ecotoxicology to investigate these early signs of genetic adaptation.

An RNA-sequencing transcriptome of the rodent Schwann cell response to peripheral nerve injury

BACKGROUND

The important contribution of glia to mechanisms of injury and repair of the nervous system is increasingly recognized. In stark contrast to the central nervous system (CNS), the peripheral nervous system (PNS) has a remarkable capacity for regeneration after injury. Schwann cells are recognized as key contributors to PNS regeneration, but the molecular underpinnings of the Schwann cell response to injury and how they interact with the inflammatory response remain incompletely understood.

METHODS

We completed bulk RNA-sequencing of Schwann cells purified acutely using immunopanning from the naïve and injured rodent sciatic nerve at 3, 5, and 7 days post-injury. We used qRT-PCR and in situ hybridization to assess cell purity and probe dataset integrity. Finally, we used bioinformatic analysis to probe Schwann cell-specific injury-induced modulation of cellular pathways.

RESULTS

Our data confirm Schwann cell purity and validate RNAseq dataset integrity. Bioinformatic analysis identifies discrete modules of genes that follow distinct patterns of regulation in the 1st days after injury and their corresponding molecular pathways. These findings enable improved differentiation of myeloid and glial components of neuroinflammation after peripheral nerve injury and highlight novel molecular aspects of the Schwann cell injury response such as acute downregulation of the AGE/RAGE pathway and of secreted molecules Sparcl1 and Sema5a.

CONCLUSIONS

We provide a helpful resource for further deciphering the Schwann cell injury response and a depth of transcriptional data that can complement the findings of recent single cell sequencing approaches. As more data become available on the response of CNS glia to injury, we anticipate that this dataset will provide a valuable platform for understanding key differences in the PNS and CNS glial responses to injury and for designing approaches to ameliorate CNS regeneration.

A Chromosome-Level Genome Assembly of the European Beech (Fagus sylvatica) Reveals Anomalies for Organelle DNA Integration, Repeat Content and Distribution of SNPs

The European Beech is the dominant climax tree in most regions of Central Europe and valued for its ecological versatility and hardwood timber. Even though a draft genome has been published recently, higher resolution is required for studying aspects of genome architecture and recombination. Here, we present a chromosome-level assembly of the more than 300 year-old reference individual, Bhaga, from the Kellerwald-Edersee National Park (Germany). Its nuclear genome of 541 Mb was resolved into 12 chromosomes varying in length between 28 and 73 Mb. Multiple nuclear insertions of parts of the chloroplast genome were observed, with one region on chromosome 11 spanning more than 2 Mb which fragments up to 54,784 bp long and covering the whole chloroplast genome were inserted randomly. Unlike in Arabidopsis thaliana, ribosomal cistrons are present in Fagus sylvatica only in four major regions, in line with FISH studies. On most assembled chromosomes, telomeric repeats were found at both ends, while centromeric repeats were found to be scattered throughout the genome apart from their main occurrence per chromosome. The genome-wide distribution of SNPs was evaluated using a second individual from Jamy Nature Reserve (Poland). SNPs, repeat elements and duplicated genes were unevenly distributed in the genomes, with one major anomaly on chromosome 4. The genome presented here adds to the available highly resolved plant genomes and we hope it will serve as a valuable basis for future research on genome architecture and for understanding the past and future of European Beech populations in a changing climate.

Growth and Wood Trait Relationships of Alnus glutinosa in Peatland Forest Stands With Contrasting Water Regimes

Human-driven peatland drainage has occurred in Europe for centuries, causing habitat degradation and leading to the emission of greenhouse gases. As such, in the last decades, there has been an increase in policies aiming at restoring these habitats through rewetting. Alder (Alnus glutinosa L.) is a widespread species in temperate forest peatlands with a seemingly high waterlogging tolerance. Yet, little is known about its specific response in growth and wood traits relevant for tree functioning when dealing with changing water table levels. In this study, we investigated the effects of rewetting and extreme flooding on alder growth and wood traits in a peatland forest in northern Germany. We took increment cores from several trees at a drained and a rewetted stand and analyzed changes in ring width, wood density, and xylem anatomical traits related to the hydraulic functioning, growth, and mechanical support for the period 1994-2018. This period included both the rewetting action and an extreme flooding event. We additionally used climate-growth and climate-density correlations to identify the stand-specific responses to climatic conditions. Our results showed that alder growth declined after an extreme flooding in the rewetted stand, whereas the opposite occurred in the drained stand. These changes were accompanied by changes in wood traits related to growth (i.e., number of vessels), but not in wood density and hydraulic-related traits. We found poor climate-growth and climate-density correlations, indicating that water table fluctuations have a stronger effect than climate on alder growth. Our results show detrimental effects on the growth of sudden water table changes leading to permanent waterlogging, but little implications for its wood density and hydraulic architecture. Rewetting actions should thus account for the loss of carbon allocation into wood and ensure suitable conditions for alder growth in temperate peatland forests.

Co-Inference of Data Mislabelings Reveals Improved Models in Genomics and Breast Cancer Diagnostics

Mislabeling of cases as well as controls in case–control studies is a frequent source of strong bias in prognostic and diagnostic tests and algorithms. Common data processing methods available to the researchers in the biomedical community do not allow for consistent and robust treatment of labeled data in the situations where both, the case and the control groups, contain a non-negligible proportion of mislabeled data instances. This is an especially prominent issue in studies regarding late-onset conditions, where individuals who may convert to cases may populate the control group, and for screening studies that often have high false-positive/-negative rates. To address this problem, we propose a method for a simultaneous robust inference of Lasso reduced discriminative models and of latent group-specific mislabeling risks, not requiring any exactly labeled data. We apply it to a standard breast cancer imaging dataset and infer the mislabeling probabilities (being rates of false-negative and false-positive core-needle biopsies) together with a small set of simple diagnostic rules, outperforming the state-of-the-art BI-RADS diagnostics on these data. The inferred mislabeling rates for breast cancer biopsies agree with the published purely empirical studies. Applying the method to human genomic data from a healthy-ageing cohort reveals a previously unreported compact combination of single-nucleotide polymorphisms that are strongly associated with a healthy-ageing phenotype for Caucasians. It determines that 7.5% of Caucasians in the 1000 Genomes dataset (selected as a control group) carry a pattern characteristic of healthy ageing.