De novo transcriptome assembly and analysis of differential gene expression in response to drought in European beech

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.

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.

Transcriptome profiling, physiological, and biochemical analyses provide new insights towards drought stress response in sugar maple (Acer saccharum Marshall) saplings

Sugar maple (Acer saccharum Marshall) is a temperate tree species in the northeastern parts of the United States and is economically important for its hardwood and syrup production. Sugar maple trees are highly vulnerable to changing climatic conditions, especially drought, so understanding the physiological, biochemical, and molecular responses is critical. The sugar maple saplings were subjected to drought stress for 7, 14, and 21 days and physiological data collected at 7, 14, and 21 days after stress (DAS) showed significantly reduced chlorophyll and Normalized Difference Vegetation Index with increasing drought stress time. The drought stress-induced biochemical changes revealed a higher accumulation of malondialdehyde, proline, and peroxidase activity in response to drought stress. Transcriptome analysis identified a total of 14,099 differentially expressed genes (DEGs); 328 were common among all stress periods. Among the DEGs, transcription factors (including NAC, HSF, ZFPs, GRFs, and ERF), chloroplast-related and stress-responsive genes such as peroxidases, membrane transporters, kinases, and protein detoxifiers were predominant. GO enrichment and KEGG pathway analysis revealed significantly enriched processes related to protein phosphorylation, transmembrane transport, nucleic acids, and metabolic, secondary metabolite biosynthesis pathways, circadian rhythm-plant, and carotenoid biosynthesis in response to drought stress. Time-series transcriptomic analysis revealed changes in gene regulation patterns in eight different clusters, and pathway analysis by individual clusters revealed a hub of stress-responsive pathways. In addition, qRT-PCR validation of selected DEGs revealed that the expression patterns were consistent with transcriptome analysis. The results from this study provide insights into the dynamics of physiological, biochemical, and gene responses to progressive drought stress and reveal the important stress-adaptive mechanisms of sugar maple saplings in response to drought stress.

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.

Response of Poplar Leaf Transcriptome to Changed Management and Environmental Conditions in Pure and Mixed with Black Locust Stands

Mixed cropping in short rotation coppice can be an alternative to monocultures. To design optimized mixtures, field trials are needed. Poplar, as an economically important and fast-growing species, and black locust, as a nitrogen-fixing species, are promising candidates for such studies. RNA sequencing (RNA-seq) was used to monitor effects of mixed and pure cultivations on the gene expression of poplar along with growth measurements during 2017 and 2018. Both biomass production and leaf transcriptomes revealed a strong competition pressure of black locust and the abiotic environment on poplar trees. Gene expression differed between the two study sites and pure and mixed stands. Shading effects from black locust caused the downregulation of photosynthesis and upregulation of shade avoidance genes in mixed stands in 2017. As a result of higher light availability after cutting black locust, plant organ development genes were upregulated in mixed stands in 2018. Drought conditions during the summer of 2018 and competition for water between the two species caused the upregulation of drought stress response genes in mixed stands and at the unfavorable growing site. Further investigations are required to discover the mechanisms of interspecific competition and to develop stand designs, which could increase the success and productivity of mixed plantations.

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

In the course of global climate change, Central Europe is experiencing more frequent and prolonged periods of drought. The drought years 2018 and 2019 affected European beeches (Fagus sylvatica L.) differently: even in the same stand, drought-damaged trees neighboured healthy trees, suggesting that the genotype rather than the environment was responsible for this conspicuous pattern. We used this natural experiment to study the genomic basis of drought resistance with Pool-GWAS. Contrasting the extreme phenotypes identified 106 significantly associated single-nucleotide polymorphisms (SNPs) throughout the genome. Most annotated genes with associated SNPs (>70%) were previously implicated in the drought reaction of plants. Non-synonymous substitutions led either to a functional amino acid exchange or premature termination. An SNP assay with 70 loci allowed predicting drought phenotype in 98.6% of a validation sample of 92 trees. Drought resistance in European beech is a moderately polygenic trait that should respond well to natural selection, selective management, and breeding.

Climate change is having a serious impact on many ecosystems. In the summer of 2018 and 2019, around two thirds of European beech trees were damaged or killed by extreme drought. It is critical to keep these beech woods healthy, as they are central to the survival of over 6,000 other species of animals and plants.

The level of damage caused by the drought varied between forests. However, not all the trees in each forest responded in the same way, with severely damaged trees often sitting next to fully healthy ones. This suggests that the genetic make-up of each tree determines how well it can adapt to drought rather than its local environment.

To investigate this further, Pfenninger et al. studied the genome of over 400 European beech trees from the Hesse region in Germany. The samples came from pairs of neighbouring trees that had responded differently to the droughts. The analysis found more than 80 parts of the genome that differed between healthy and damaged trees.

Pfenninger et al. then used this information to create a genetic test which can quickly and inexpensively predict how well an individual beech tree might survive in a drought. Applying this test to another 92 trees revealed that it can reliably detect which ones were healthy and which ones were damaged.

Beech forests are typically managed by private owners, agencies or breeders that could use this genetic test to select and reproduce trees that are better adapted to drought. The goal now is to develop the test so that it can be used more widely to manage European beech trees and potentially other species.

Time-course transcriptomic analysis of Petunia ×hybrida leaves under water deficit stress using RNA sequencing

Water deficit limits plant growth and development, resulting in quality loss of horticultural crops. However, there is limited information on gene regulation and signaling pathways related to water deficit stress response at multiple time points. The objective of this research was to investigate global gene expression patterns under water deficit stress to provide an insight into how petunia (Petunia ×hybrida 'Mitchell Diploid') responded in the process of stress. Nine-week-old petunias were irrigated daily or placed under water stress by withholding water. Stressed plants reduced stomatal conductance after five days of water deficit, indicating they perceived stress and initiated stress response mechanisms. To analyze transcriptomic changes at the early stage of water deficit, leaf tissue samples were collected 1, 3, and 5 days after water was withheld for RNA sequencing. Under water deficit stress, 154, 3611, and 980 genes were upregulated and 41, 2806, and 253 genes were downregulated on day 1, 3, and 5, respectively. Gene Ontology analysis revealed that redox homeostasis processes through sulfur and glutathione metabolism pathways, and hormone signal transduction, especially abscisic acid and ethylene, were enriched under water deficit stress. Thirty-four transcription factor families were identified, including members of AP2/ERF, NAC, MYB-related, C2H2, and bZIP families, and TFs in AP2/ERF family was the most abundant in petunia. Interestingly, only one member of GRFs was upregulated on day 1, while most of TFs were differentially expressed on day 3 and/or 5. The transcriptome data from this research will provide valuable molecular resources for understanding the early stages of water stress-responsive networks as well as engineering petunia with enhanced water stress tolerance.

TOA: A software package for automated functional annotation in non-model plant species

The increase of sequencing capacity provided by high-throughput platforms has made it possible to routinely obtain large sets of genomic and transcriptomic sequences from model and non-model organisms. Subsequent genomic analysis and gene discovery in next-generation sequencing experiments are, however, bottlenecked by functional annotation. One common way to perform functional annotation of sets of sequences obtained from next-generation sequencing experiments, is by searching for homologous sequences and accessing the related functional information deposited in genomic databases. Functional annotation is especially challenging for non-model organisms, like many plant species. In such cases, existing free and commercial general-purpose applications may not offer complete and accurate results. We present TOA (Taxonomy-oriented annotation), a Python-based user-friendly open source application designed to establish functional annotation pipelines geared towards non-model plant species that can run in Linux/Mac computers, HPCs and cloud servers. TOA performs homology searches against proteins stored in the PLAZA databases, NCBI RefSeq Plant, Nucleotide Database and Non-Redundant Protein Sequence Database, and outputs functional information from several ontology systems: Gene Ontology, InterPro, EC, KEGG, Mapman and MetaCyc. The software performance was validated by comparing the runtimes, total number of annotated sequences and accuracy of the functional information obtained for several plant benchmark data sets with TOA and other functional annotation solutions. TOA outperformed the other software in terms of number of annotated sequences and accuracy of the annotation and constitutes a good alternative to improve functional annotation in plants. TOA is especially recommended for gymnosperms or for low quality sequence data sets of non-model plants.

Transcriptome profiles of Quercus rubra responding to increased O3 stress

BACKGROUND

Climate plays an essential role in forest health, and climate change may increase forest productivity losses due to abiotic and biotic stress. Increased temperature leads to the increased formation of ozone (O3). Ozone is formed by the interaction of sunlight, molecular oxygen and by the reactions of chemicals commonly found in industrial and automobile emissions such as nitrogen oxides and volatile organic compounds. Although it is well known that productivity of Northern red oak (Quercus rubra) (NRO), an ecologically and economically important species in the forests of eastern North America, is reduced by exposure to O3, limited information is available on its responses to exogenous stimuli at the level of gene expression.

RESULTS

RNA sequencing yielded more than 323 million high-quality raw sequence reads. De novo assembly generated 52,662 unigenes, of which more than 42,000 sequences could be annotated through homology-based searches. A total of 4140 differential expressed genes (DEGs) were detected in response to O3 stress, as compared to their respective controls. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of the O3-response DEGs revealed perturbation of several biological pathways including energy, lipid, amino acid, carbohydrate and terpenoid metabolism as well as plant-pathogen interaction.

CONCLUSION

This study provides the first reference transcriptome for NRO and initial insights into the genomic responses of NRO to O3. Gene expression profiling reveals altered primary and secondary metabolism of NRO seedlings, including known defense responses such as terpenoid biosynthesis.

Abiotic genetic adaptation in the Fagaceae

Fagaceae can be found in tropical and temperate regions and contain species of major ecological and economic importance. In times of global climate change, tree populations need to adapt to rapidly changing environmental conditions. The predicted warmer and drier conditions will potentially result in locally maladapted populations. There is evidence that major genera of the Fagaceae are already negatively affected by climate change-related factors such as drought and associated biotic stressors. Therefore, knowledge of the mechanisms underlying adaptation is of great interest. In this review, we summarise current literature related to genetic adaptation to abiotic environmental conditions. We begin with an overview of genetic diversity in Fagaceae species and then summarise current knowledge related to drought stress tolerance, bud burst timing and frost tolerance in the Fagaceae. Finally, we discuss the role of hybridisation, epigenetics and phenotypic plasticity in adaptation.

Re-assembly, quality evaluation, and annotation of 678 microbial eukaryotic reference transcriptomes

BACKGROUND

De novo transcriptome assemblies are required prior to analyzing RNA sequencing data from a species without an existing reference genome or transcriptome. Despite the prevalence of transcriptomic studies, the effects of using different workflows, or "pipelines," on the resulting assemblies are poorly understood. Here, a pipeline was programmatically automated and used to assemble and annotate raw transcriptomic short-read data collected as part of the Marine Microbial Eukaryotic Transcriptome Sequencing Project. The resulting transcriptome assemblies were evaluated and compared against assemblies that were previously generated with a different pipeline developed by the National Center for Genome Research.

RESULTS

New transcriptome assemblies contained the majority of previous contigs as well as new content. On average, 7.8% of the annotated contigs in the new assemblies were novel gene names not found in the previous assemblies. Taxonomic trends were observed in the assembly metrics. Assemblies from the Dinoflagellata showed a higher number of contigs and unique k-mers than transcriptomes from other phyla, while assemblies from Ciliophora had a lower percentage of open reading frames compared to other phyla.

CONCLUSIONS

Given current bioinformatics approaches, there is no single "best" reference transcriptome for a particular set of raw data. As the optimum transcriptome is a moving target, improving (or not) with new tools and approaches, automated and programmable pipelines are invaluable for managing the computationally intensive tasks required for re-processing large sets of samples with revised pipelines and ensuring a common evaluation workflow is applied to all samples. Thus, re-assembling existing data with new tools using automated and programmable pipelines may yield more accurate identification of taxon-specific trends across samples in addition to novel and useful products for the community.

Physiological and Transcriptomic Analysis of Tree Peony (Paeonia section Moutan DC.) in Response to Drought Stress

Tree peony (Paeonia section Moutan DC.) is a famous ornamental plant, and P. ostii has been used for seed oil production in China because it is rich in α-linolenic acid. P. ostii has some resistance to drought, but lack of water can severely hinder its growth and development in arid areas. In order to clarify drought stress induced physiological and molecular changes of P. ostia, its physiological and transcriptomic analyses were performed under drought stress, and we found that P. ostii leaves drooped significantly 12 days after treatment and observed a significant increase in all detected physiological indices in response to drought treatment except leaf water content, chlorophyll, and carotenoid content. Meanwhile, the activity of three antioxidant enzymes basically increased under drought treatment. Moreover, drought treatment significantly reduced photosynthetic and chlorophyll fluorescence parameters except non-photochemical quenching (qN), and maintained more intact mesophyll cell structures. Additionally, many differentially expressed genes (DEGs) were found by transcriptome sequencing, which play an important role in P. ostia drought tolerance by controlling a variety of biological processes, including the reactive oxygen species (ROS) system, chlorophyll degradation and photosynthetic competency, fatty acid metabolism, proline metabolism, biosynthesis of secondary metabolism, and plant hormone metabolism. These results provide a better understanding of P. ostii responses to drought stress.

A reference genome of the European beech (Fagus sylvatica L.)

Background

The European beech is arguably the most important climax broad-leaved tree species in Central Europe, widely planted for its valuable wood. Here, we report the 542 Mb draft genome sequence of an up to 300-year-old individual (Bhaga) from an undisturbed stand in the Kellerwald-Edersee National Park in central Germany.

Findings

Using a hybrid assembly approach, Illumina reads with short- and long-insert libraries, coupled with long Pacific Biosciences reads, we obtained an assembled genome size of 542 Mb, in line with flow cytometric genome size estimation. The largest scaffold was of 1.15 Mb, the N50 length was 145 kb, and the L50 count was 983. The assembly contained 0.12% of Ns. A Benchmarking with Universal Single-Copy Orthologs (BUSCO) analysis retrieved 94% complete BUSCO genes, well in the range of other high-quality draft genomes of trees. A total of 62,012 protein-coding genes were predicted, assisted by transcriptome sequencing. In addition, we are reporting an efficient method for extracting high-molecular-weight DNA from dormant buds, by which contamination by environmental bacteria and fungi was kept at a minimum.

Conclusions

The assembled genome will be a valuable resource and reference for future population genomics studies on the evolution and past climate change adaptation of beech and will be helpful for identifying genes, e.g., involved in drought tolerance, in order to select and breed individuals to adapt forestry to climate change in Europe. A continuously updated genome browser and download page can be accessed from beechgenome.net, which will include future genome versions of the reference individual Bhaga, as new sequencing approaches develop.

Genome-Wide Identification and Functional Prediction of Novel Drought-Responsive lncRNAs in Pyrus betulifolia

Increasing evidence shows that long noncoding RNAs (lncRNAs) play important roles in developmental regulation and many other biological processes in plants. However, identification of lncRNAs in Pyrus betulifolia is limited compared with studies of functional gene expression. Using high-throughput sequencing technology, the transcriptome of P. betulifolia under drought stress was analyzed to identify lncRNAs. A total of 14,478 lncRNAs were identified, of which 251 were found to be drought-responsive. The putative target genes of these differentially expressed lncRNAs were significantly enriched in metabolic processes, organic substance metabolic processes, macromolecule metabolic processes, and heterocyclic compound binding. Real-time quantitative polymerase chain reaction validation suggested that the results of the RNA sequencing data analysis were reliable. This study will provide genetic resources for pear breeding and provide reference to other pomological studies.