Next Article in Journal
Electroencephalography (EEG)-Based Neural Emotional Response to the Vegetation Density and Integrated Sound Environment in a Green Space
Next Article in Special Issue
Effects of Early, Small-Scale Nitrogen Addition on Germination and Early Growth of Scots Pine (Pinus sylvestris) Seedlings and on the Recruitment of the Root-Associated Fungal Community
Previous Article in Journal
Perceptions by Smallholder Farmers of Forest Plantations in Malaysia
Previous Article in Special Issue
High Seedling Mortality of Scots Pine Caused by Heterobasidion annosum s.s.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Drought Stress Described by Transcriptional Responses of Picea abies (L.) H. Karst. under Pathogen Heterobasidion parviporum Attack

by
Xenia Hao-Yi Yeoh
1,2,
Blessing Durodola
1,
Kathrin Blumenstein
1 and
Eeva Terhonen
1,3,*
1
Forest Pathology Research Group, Department of Forest Botany and Tree Physiology, Büsgen-Institute, Faculty of Forest Sciences and Forest Ecology, University of Göttingen, Büsgenweg 2, 37077 Göttingen, Germany
2
Faculty of Medicine, Technical University of Munich, Klinikum Rechts der Isar, Ismaninger Straße 22, 81675 Munich, Germany
3
Natural Resources Institute Finland (Luke), Forest Health and Biodiversity, Latokartanonkaari 9, FI-00790 Helsinki, Finland
*
Author to whom correspondence should be addressed.
Forests 2021, 12(10), 1379; https://doi.org/10.3390/f12101379
Submission received: 10 September 2021 / Revised: 1 October 2021 / Accepted: 7 October 2021 / Published: 10 October 2021
(This article belongs to the Special Issue Tree Host – Microbial Interactions)

Abstract

:
The major threats to the sustainable supply of forest tree products are adverse climate, pests and diseases. Climate change, exemplified by increased drought, poses a unique threat to global forest health. This is attributed to the unpredictable behavior of forest pathosystems, which can favor fungal pathogens over the host under persistent drought stress conditions in the future. Currently, the effects of drought on tree resistance against pathogens are hypothetical, thus research is needed to identify these correlations. Norway spruce (Picea abies (L.) H. Karst.) is one of the most economically important tree species in Europe and is considered highly vulnerable to changes in climate. Dedicated experiments to investigate how disturbances will affect the Norway spruce—Heterobasidion sp. pathosystem are important, in order to develop different strategies to limit the spread of H. annosum s.l. under the predicted climate change. Here, we report a transcriptional study to compare Norway spruce gene expressions to evaluate the effects of water availability and the infection of Heterobasidion parviporum. We performed inoculation studies of three-year-old saplings in a greenhouse (purchased from a nursery). Norway spruce saplings were treated in either high (+) or low (−) water groups: high water group received double the water amount than the low water group. RNA was extracted and sequenced. Similarly, we quantified gene expression levels of candidate genes in biotic stress and jasmonic acid (JA) signaling pathways using qRT-PCR, through which we discovered a unique preferential defense response of H. parviporum-infected Norway spruce under drought stress at the molecular level. Disturbances related to water availability, especially low water conditions can have negative effects on the tree host and benefit the infection ability of the pathogens in the host. From our RNA-seq analysis, 114 differentially expressed gene regions were identified between high (+) and low (−) water groups under pathogen attack. None of these gene pathways were identified to be differentially expressed from both non-treated and mock-control treatments between high (+) and low (−) water groups. Finally, only four genes were found to be associated with drought in all treatments.

1. Introduction

Norway spruce (Picea abies (L.) H. Karst.) is one of the most economically important forestry tree species in Europe. Currently, massive inputs of anthropogenic greenhouse gases (i.e., CO2, N2O, CH4) into the atmosphere have resulted in increasing atmospheric temperatures in an effect known as ‘global warming’ [1]. This effect has several consequences; one of which, climate change, is a potential driver in influencing forest health. Under global warming, extreme climate events, e.g., drought, are expected to increase in frequency, duration and intensity [2]. In fact, this has started to take its toll on Norway spruce that is sensitive to abiotic disturbances. Under climate change, there is an increase in drought-associated stress in this particular tree species, rendering them more susceptible to threats, such as pests and pathogens, which can compromise overall tree health [3,4,5]. This phenomenon can be explained by the disease triangle model that shows the relationship between three factors: (1) a susceptible plant host, (2) a pathogen that causes the disease in the plant host, and (3) an environment that favors the pathogen [6,7]. An environment that weakens the defense capability of the plant host and supports the growth and spread of the pathogen can greatly influence the success rate for a pathogen to infect its targeted host. This means that changes in climate do not limit the new challenges of Norway spruce to adjust only to abiotic disturbances, as the projected change in climate will favor certain pathogens in forests [2,5,8,9,10]. The change in tree resistance under global warming against fungal pathogens still remains hypothetical, thus research is needed to identify these correlations. Discriminating responses at the genomic level between biotic and abiotic stressors and understanding the potential trade-offs in plant gene expression will allow for more mechanistic predictions about future scenarios driven by complex shifts in climate change [11,12].
Heterobasidion parviporum Niemelä and Korhonen is one of the native fungal pathogen species (the other two being H. annosum s. str. (Fr.:Fr.) Bref. and H. abietinum Niemelä and Korhonen) belonging to a root-rot fungus, H. annosum sensu lato (s.l.) complex, which is considered to be the most devastating Norway spruce pathogen in Europe [13]. Being necrotrophic in nature, it mainly derives its nutrients by killing its tree host and feeding off its contents. There are two modes in which H. parviporum can infect Norway spruce: (1) through deposition of basidiospores on freshly cut stumps, or (2) via root contacts [14]. Detection of H. parviporum in Norway spruce is often challenging, as the infected tree might not necessarily show signs of symptoms (stringy white rot, necrosis in sapwood, presence of fruiting bodies) over decades. This further promotes the spread of H. parviporum to neighboring trees via root contacts. In addition, H. parviporum mainly colonizes the heartwood of Norway spruce, causing decay and even tree mortality [15].
The induction of a disease by a pathogen in a tree host is often linked to the environmental conditions in the ecosystem as well. In this particular Norway spruce—H. parviporum forestry pathosystem, climate plays an important role in promoting the spread of H. parviporum in several ways: (1) short winters that provide a longer period of hyphae growth and spore transmission, (2) warmer summers that encourage fungal hyphae to grow faster and produce more spores, and (3) windy days that can promote the dispersion of basidiospores [2,16,17]. Increased frequency and intensity of drought-associated stress in Norway spruce is expected in the following years [18]. Considering how Norway spruce is a valuable tree species to the boreal regions, it is imperative to study the interactions between this tree species and its associated fungal pathogens. This way, appropriate and effective measures can be implemented to mitigate this inevitable climate change.
Rapid advancements in molecular research have provided insights into the Norway spruce—H. parviporum forestry pathosystem through various studies, such as pathogen characterization and detection [19], constructing gene expression libraries [20] and screening of resistant trees based on genetic components [21]. Norway spruce resistance against H. annosum s.l. is a quantitative trait [21,22,23] and is associated with several known genes (defense responses) with variation [21,24]. Signaling molecules are a group of plant secondary metabolites with a critical role in defense processes [25]. Plant hormones, especially salicylic acid (SA) and jasmonic acid (JA) help mediate information about the attack beyond the point of invasion [25]. By altering the hormonal balance in plants, pathogens use the defensive machinery of plants to their advantage and either induce or suppress the processes relevant for cell death and accumulation of antimicrobial compounds [25,26]. Jasmonic acid (JA)-mediated signaling is the prioritized module in the Norway spruce defense-signaling network against H. parviporum [27,28,29]. Similarly, besides JA, the ethylene (ET) signaling pathway probably also plays a central role in the defense response of Norway spruce against H. parviporum [28]. In addition, the salicylic acid (SA) signaling pathway can, likewise, have an important role in this pathosystem as the SA-mediated hypersensitive response (HR) can be facilitated by necrotrophic pathogen infection [30]. Understanding the similarities and differences in transcriptional response between variable environmental conditions can thus provide useful information over how Norway spruce—H. parviporum pathosystems adapt to climate change.
We could show previously that H. parviporum benefits from host stress (i.e., drought) [18]. In this study, we performed RNA-seq analysis on samples derived from the same trees as in Terhonen et al. [18]. We concentrated on Norway spruce, to find out which genetic pathways are involved in the specific responses to unfavorable environmental stress (drought) with and without the pathogen, H. parviporum. In addition, we also built a transcriptional study based on previous studies [27,28,31,32,33] to investigate the interplay between H. parviporum infection and defense-related gene expression mediated by drought stress.

2. Materials and Methods

2.1. Plant and Fungal Material

The setting used in this experiment is from Terhonen et al. [18]. Summarily, plant material consisted of 18 three-year-old, apparently healthy and vital Norway spruce saplings purchased from the nursery Schlegel and Co Gartenprodukte (provenience of the saplings: No. 840 11, Thüringer Wald and Frankenwald, montane zone 600 m). Seedlings were potted into 3-L plastic pots filled with fertilized peat (Flora gard, TKS®2 Instant Plus, Hermann Meyer KG, Rellingen, Germany). The potted saplings were then acclimatized to the greenhouse conditions for 16-days prior to the water experiment, during which they received tap water, as required, to maintain moist soil. No additional fertilization was given during the experiment.
A H. parviporum strain was received from the North-West German Forest Research Institute, collected by G. Langer and colleagues. H. parviporum (strain NW-FVA 0459) was isolated from P. abies in 2010 (Germany, Lower Saxony, forest department Oerrel, Bobenwald). DNA of H. parviporum was extracted from 150 mg of the homogenized mycelium sample using the method developed by Keriö et al. [34] and the fungal strain used in this study was tested for the specificity with species-specific primers developed for H. parviporum (KJ-F and KJ-R; [19]). In brief, DNA template (100 ng), buffer (KCl extra buffer, 1X), 1.5 mM MgCl2, primers KJ-F and KJ-R (each concentration of 0.5 µM), a dNTP-mix (each deoxynucleotide in a concentration of 200 µM) and 20 U/mL of DNA-polymerase (VWR) was adjusted to 25 µL reaction with autoclaved Milli-Q H2O. The cycling conditions used were: Initial denaturation 10 min at 95 °C, followed by a 3-step cycling: Denaturation 30 s at 95 °C; annealing 35 s at 67 °C; extension 1 min at 72 °C, for 40 cycles. A final extension step for 7 min at 72 °C was applied. Taq DNA polymerase (Qiagen) was used for PCR amplification of ITS regions with the primer pair, ITS1-F and ITS4. Briefly, the PCR protocol was as follows: 1X CoralLoad PCR Buffer, 200 µM dNTP, 0.5 µM primer 1, 0.5 µM primer 2, 100ng template DNA, 0.2 U/µL DNA polymerase; the reaction was adjusted to 25 µL with autoclaved Milli-Q H2O. The PCR conditions used were 94 °C for 3 min; 30 cycles of 94 °C for 30 s, 55 °C for 1 min, 72 °C for 1 min; a final step at 72 °C for 10 min was applied before storing the PCR products at 4 °C. Possible contaminations were determined with a negative control using sterile water as a template in both PCR protocols. StainIN™ RED Nucleic Acid Stain was used to confirm DNA amplicons on a 2% agarose gel and the visual detection was made by ultraviolet transillumination. The expected band (fragment size around 350 bp) with KJ-F and KJ-R primers was visualized and the ITS region PCR products were purified and sequenced using the ITS4 primer at Microsynth SEQLAB (Göttingen, Germany). The FASTA files thus obtained were checked with BioEdit to confirm that the pathogen was not contaminated with other fungi strains. ITS region sequencing indicated no contamination in the culture of H. parviporum (see [18]).

2.2. Experimental Design

The following experiment is described in detail in [18]. In summary, the study was conducted at the Forest Botany and Tree Physiology greenhouses, in Göttingen, Germany (51°33′28.4″ N 9°57′30.5″ E) from early April until early August 2018. From this experiment, we chose 18 seedlings (three per water treatment group) randomly block-assigned to waterproof tables with either high water (high water group = +group) or low water (low water group = −group) availability treatments, which were experimentally inoculated with H. parviporum, mock-inoculated controls or left entirely untreated. The water treatment experiment was running for 35 days before the inoculations were performed. H. parviporum isolate (NW-FVA 0459) was inoculated using a puncher (Ø5 mm) through the bark to reach the sapwood surface. Equal-sized plugs from a pure culture of H. parviporum or control (2% MEA) were placed onto the exposed surface and sealed with Parafilm® [18]. The inoculation experiment was run for 70 days [18]. At the end of the experiment, we collected 30 samples aiming to sequence 18 samples. The stems of Norway spruce saplings (5 cm in both directions: above and below the infection point was collected, around 11 cm) were sampled into liquid nitrogen and stored individually into 15-mL Falcon tubes at −80 °C for long-term storage.

2.3. RNA Extraction

A total of 30 Norway spruce stems were selected, comprising of five sample replicates for each of the six treatment groups for small-scale totRNA extractions; a total of 18 (with three sample replicates for each of the six treatment groups) were selected for upscale totRNA extractions. For the first round of extractions, the surface of the inoculation point (bark, phloem, sapwood) was scraped with a scalpel and the sample was ground into a fine powder (Mixer Mill MM 400 from Retsch GmbH with a set program of 25.0 Hz for 20 s). The samples were handled with liquid nitrogen throughout the entire milling process. The ground product was then stored in 15-mL Falcon tubes at −80 °C for later use in the extractions. The extraction protocols used were modified from Zeng et al. [35] and unfortunately, the first try with 0.1–0.5 g of sample transferred to 2-mL Eppendorf tubes was unsuccessful as the RNA quantity was too low for sequencing. We subsequently re-did the scraping from the infection point and increased the area to collect more material (the sample included bark, phloem and sapwood). The new samples were similarly ground into a fine powder and extracted for totRNA, which was successful. From this, we established the following extraction protocol, which we had modified from Zeng et al. [35] for our study.
Day 1: 2–3 g of sample was transferred to a 50-mL Falcon tube, where 15 mL of Extraction buffer (2% CTAB, 2% PVP, 100 mM Tris-HCl, 25 mM EDTA, 2 M NaCl, DEPC-treated H2O, pre-heated to 65 °C) was added; 100 μL of β-mercaptoethanol was added to each tube and vigorously vortexed before incubating the tubes at 65 °C for 5 min. An equal volume of chloroform:isoamyl alcohol (24:1) was added to each tube and vortexed vigorously before centrifuging the tubes at 10,000× g for 15 min. The upper aqueous phase was transferred to a new sterile 50-mL Falcon tube. Chloroform:isoamyl alcohol was added, then the tubes were centrifuged and subsequently the transfer of the upper aqueous phase was repeated; 1/4 volume of 10 M LiCl was added to each tube and inverted eight times before precipitating the RNA overnight at 4 °C.
Day 2: The tubes were centrifuged using the benchtop Eppendorf 5810R Centrifuge (Rotor A-4-81 with a bucket for 7 × 50 mL conical tubes) with a speed of approximately 6480 rcf at 4 °C for 30 min. The supernatant was collected in a new 50-mL Falcon tube and set aside for gDNA extraction. The remaining pellet (containing totRNA) was briefly dried in a laminar flow hood; 1 mL of SSTE (1 M NaCl, 0.5% SDS, 10 mM Tris-HCl, 1 mM EDTA, DEPC-treated H2O, preheated to 65 °C) was then added to dissolve the pellet. Thereafter, an equal volume of chloroform:isoamyl alcohol was added and vortexed vigorously before centrifuging at 10,000× g for 15 min. The upper aqueous phase (maximum 500 μL) was transferred to a new 2-mL Eppendorf tube and three volumes of cold 100% absolute EtOH were added to each tube. The resultant mixture is incubated at −20 °C for 1–2 h before transferring to −80 °C overnight for RNA precipitation.
Day 3: totRNA samples were centrifuged at maximum centrifugal speed at 4 °C for 30 min. The supernatant was pipetted out and the remaining totRNA pellet was dried in a sterile hood for approximately 3 min. The pellet was then resuspended in 100 μL of DEPC-treated H2O, where it could be subsequently used in molecular analyses or stored at −80 °C for long-term storage.

2.4. Quantitation and Purity Analyses of Extracted totRNA

A Qubit 3.0 Fluorometer from Life Technologies (Catalogue number: Q33216) was used to quantitate totRNA. The analysis was performed according to the manufacturer’s instructions. Three biological replicates from each treatment (2 µg of total RNA) were sent for RNA sequencing (RNAseq) to Novogene. RNA integrity and quantitation were assessed using the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, CA, USA). Purified RNA was used for library construction using rRNA depletion method (GlobinZero kit) for Illumina® (NEB, Santiago, CA, USA), followed by Illumina paired-end (2  ×  150 bp) RNAseq in the facilities of the Novogene (HK) company limited. The raw reads have been submitted to NCBI SRA under Bioproject PRJNA761217.

2.5. RNA-seq Analysis

Quality control was made by Novogene through in-house scripts. Briefly, clean data (clean reads) were obtained by removing reads containing adapter and poly-N sequences and reads with low quality from raw data. The filtering process was as follows: (1) Discard reads with adapter contamination; (2) Discard reads when uncertain nucleotides constitute more than 10% of either read (N > 10%); (3) Discard reads when low quality nucleotides (base quality less than 20) constitute more than 50% of the read. At the same time, Q20, Q30 and GC content of the clean data were calculated. All downstream analyses were based on clean data with high quality. The quality was assessed using FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 17 February 2020)). Reference genome and gene model annotation files were downloaded from genome website browser (the processed reads were mapped against the genome assembly of Norway spruce (v 1.0) [36] downloaded from ftp://plantgenie.org/Data/ConGenIE/Picea_abies/ (accessed on 8 March 2020) directly. Paired-end clean reads were mapped to the reference genome using HISAT2 software [37]. The raw read counts table was loaded into R studio version 3.5.1 [38] and differential expression analysis between two conditions (three biological replicates per condition) was performed using DESeq2 R package [39]. Conditions were normal water versus low water in each treatment (treatments: non-treated, mock-inoculated, or H. parviporum infection). To track down only drought-related gene pathways, we compared all treatments to the non-treated (+) water treatment. The resulting p values were adjusted using Benjamini and Hochberg’s approach for controlling the False Discovery Rate (FDR). Genes with an adjusted p-value < 0.05 and the cutoff value for FDR < 0.05 were assigned as differentially expressed. We searched for statistically expressed genes between different conditions/treatments with Venny 2.1. Results were compared on genes expressed in the study by Chaudhary et al. [40]. When we found interesting genes, we visualized the expression in each treatment. Similarly, we studied how H. parviporum resistance candidate gene laccase PaLAC5 (MA_97119g0010 and MA_97119g0020) were expressed in each treatment [24].

2.6. GO Enrichment Test

The GO functions and enrichment analysis were performed using ConGenIE.org (http://congenie.org/ (accessed on 22 July 2020)). Annotate GeneList workflow was created on the ConGenIE platform selecting P. abies as the species of interest. The gene IDs were pasted in the search field and annotated. The enrichment analysis was carried out using the functional enrichment tools. The GO slims were run and based on the annotations, the transcripts were categorized in the biological process and molecular function groups with their respective descriptions. The significantly enriched GO terms were then categorized into Biological Process, Molecular Function and Cellular Component.

2.7. Transcriptional Responses Validated with Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR)

For cDNA preparation, we followed the protocol described by Raffaello and Asiegbu [41]. In this experiment, four genes were tested: LOX [27], ERF1 [28], CHIIV [33] and p/DIR32 [32]. A housekeeping gene, Elongation factor 1α (ELF1α) [31], was used as a reference against the tested genes. To prepare the loading samples for qRT-PCR, gene mastermixes comprising of 0.5 µM each of forward and reverse primers for each selected gene, qPCR Supermix (SsoFast™ EvaGreen® Supermix (BIO-RAD)), and DEPC-H2O as a diluent were made. A total of 20 µL per qRT-PCR reaction (including the pooled cDNA) were then loaded onto a 96-well V-bottom plate and sealed tight with plate foil. The qRT-PCR (iQ5™ Multicolor Real-Time PCR Detection System (BIO-RAD)) cycling conditions ran were set for an initial 2 min at 95 °C, then 40 cycles of 20 s at 95 °C, 30 s at 60 °C, and 30 s at 72 °C. According to a primer efficiency test prior to running the qRT-PCR, we decided on a cDNA concentration of 100 ng per pooled cDNA. The relative expressions were subsequently calculated with the method 2−ΔΔCt [42].

3. Results

3.1. RNA-seq

Our analysis of RNA-seq data identified 114 transcripts with a significantly different expression between high (+) and low (−) water groups under pathogen attack (Supplementary File S1).
Comparing H. parviporum infection under drought stress, the biological processes (in upregulated) significantly enriched were metabolic process, transport, protein folding, and translation. Molecular function GO terms in pathogen-infected seedlings were significant in transferase activity, hydrolase activity, and binding. Cellular component GO terms were significantly related to plant cell organelles (apoplast, chloroplast, Golgi apparatus, mitochondrion, cell wall and plasma membrane) (Supplementary File S1). The roles played by these genes are expressed as GO terms.
Drought substantially increased gene expression under H. parviporum inoculation (Table 1). Exclusively 377 genes were downregulated and 129 genes upregulated in the H. parviporum-infected (−) group (Supplementary File S1). The upregulated Gene Ontology (GO) terms assigned to the “Biological process” in the (−) water group were, protein folding, response to stress, oxidation-reduction process, rRNA processing, primary metabolic process, phosphorylation, pentose-phosphate shunt and Golgi organization. The terms assigned to the “Molecular function” were, ATP binding, oxidoreductase activity, nucleotide binding, unfolded protein binding, transferase activity and protein kinase activity. The downregulated genes expressed resulted in different gene functions and the most abundant within this group were metabolic process, oxidation-reduction process, response to wounding, response to fungus, phosphorylation, defense response to fungus, response to water deprivation, and response to oxidative stress. Only the Heterobasidion-infected (−) groups were found to be functionally enriched (Supplementary File S1).
We found four genes that were expressed significantly as higher or lower amounts (compared to non-treated (+) group) only under drought stress MA_189802g0010, MA_10265000g0010, MA_10435878g0010 and MA_178006g0010 (Figure 1, Table 2). These genes were not observed to be significantly expressed (up or downregulated) in normal water treatments. Only gene DEG MA_10435878g0010 had a role in the metabolic process (GO:0008152) (Table 2).
After comparing to dataset of Chaudhary et al. [40], we identified four candidate genes related to Heterobasidion sp. resistance (MA_110169g0010, MA_14707g0010, MA_15852g0010, MA_10427673g0020) that were significantly expressed (as higher or lower amounts compared to non-treated (+) group) only in H. parviporum infected plants (only in -group: MA_14707g0010, MA_15852g0010, MA_10427673g0020). Their expression levels in each treatment are presented in Figure 2. Palac5 genes MA_97119g0010 (Figure 3a) and MA_97119g0020 (Figure 3b) were not significantly expressed (compared to the non-treated (+) group) in H. parviporum infected plants, but these genes accumulated more under H. parviporum infection (Figure 3a,b).

3.2. qRT-PCR

Normalization between technical triplicate samples was achieved by determining the mean fold change of each gene (i.e., LOX, ERF1, CHIIV, p/DIR32) against the housekeeping gene, ELF1, for each treatment scheme (Table 3). The difference in gene expression between technical triplicates was mostly with a mean fold change value ( 2 C t ) of 1, irrespective of water treatment groups and inoculation status, indicating little non-biological variation. This is further illustrated with Figure 4, which generally demonstrated consistency between triplicates except for two treatment groups in p/DIR32 (SD: H. parviporum, normal = 0.788; mock control, low = 0.822) and one group in ERF1 (SD: mock-control, low = 0.565).
To obtain a meaningful interpretation of our results, we also compared the fold change in gene expression mediated by H. parviporum and physical wounding as singular and combinatorial factors to examine how Norway spruce responds differently under drought conditions (Table 4).
LOX (1.60) and ERF1 (2.81) were upregulated while CHIIV (0.19) and p/DIR32 (0.25) were downregulated when mediated solely by H. parviporum infection under normal water conditions. With the additional stress induced by low water conditions, CHIIV (7.16) and p/DIR32 (2.64) had an increase in transcription levels and were both upregulated as a result. The sole influence of a physical wound led to an increase in gene expression of LOX (6.92), ERF1 (2.31) and p/DIR32 (17.03) when subjected to drought stress. This was the opposite observation in CHIIV (0.91) where it is downregulated instead. With sufficient water supply, CHIIV (10.85) was abundantly expressed. Combinatorically, LOX, CHIIV and p/DIR32 were upregulated when wounded H. parviporum-infected plants were subjected to drought stress (LOX = 8.82, CHIIV = 6.50, p/DIR32 = 44.94). On the other hand, with the exception of CHIIV, the tested genes were downregulated when the plants received normal water treatment (LOX = 0.27, ERF1 = 0.87, p/DIR32 = 0.08). No change was observed for ERF1 even at low water treatment conditions (ERF1 = 1.06). These observations are further represented in Figure 5.
In general, LOX and p/DIR32 showed an increase in gene expression levels under drought stress when combinatorically influenced by H. parviporum infection and physical wounding in comparison to their respective singular factors that are independent of each other. Under normal water conditions, ERF1 is dominantly expressed in H. parviporum-infected Norway spruce; this is the same for CHIIV in Norway spruce that are physically wounded. This observation is flipped under low water conditions (Figure 5).

4. Discussion

4.1. Heterobasidion Resistance Related Genes

Drought induced 114 DEGs between H. parviporum treatments (−/+), similarly exclusively 377 DEGs were observed in (−) group of H. parviporum. Chaudhary et al. [40] identified several candidate genes related to Heterobasidion sp. resistance from which we found only four (MA_110169g0010, MA_14707g0010, MA_15852g0010, MA_10427673g0020) in H. parviporum infected plants to be statistically expressed. The accumulation of these genes was induced in drought treatment in almost all cases. The same trend could be found for Palac5 genes (MA_97119g0010 and MA_97119g0020) [24] as they accumulated more in H. parviporum infected plants under drought, and in non-treated (−) group plants these genes were at the minimum level. These results highlight the need for combined factor studies, abiotic and biotic stress in plant gene expression studies to resolve this combination effect as some pathways can be induced by the environment.

4.2. Changing Environment: Drought as a New Factor in Forest Pathosystems

Climatic extremes, such as drought affect trees directly, by inducing water deficiency, and indirectly, by making trees more susceptible to fungal pathogens [10,18]. This leads to unknown behavior of host-fungi interactions, resulting in increased disease outbreaks caused by native fungal pathogens [43]. Prolonged water and nutrient deficiencies are causing severe damage to Norway spruce and the drought disturbances are expected to increase in following years [2]. Norway spruce provenances react differentially to drought [44]. Trujillo-Moya et al. [44] found drought response indicators (24 genes, 29 SNPs) of Norway spruce genotypes. One of these 24 genes, MA_588952g0010 (no similarity to any known gene), was found exclusively to be significantly expressed in H. parviporum (-) treatment. Our analysis of RNA-seq data identified only four genes (MA_189802g0010, MA_10265000g0010, MA_10435878g0010 and MA_178006g0010) significantly expressed in drought stressed trees. The gene (MA_178006g0010) is a member of the BURP domain family (RD22-like). Members of this family have been described to be associated with several plant species subjected to stress conditions [45]. In Arabidopsis thaliana (L.) Heynh., the AtRD22 and AtUSPL1 (both members of the Arabidopsis thaliana BURP) are upregulated as part of the abscisic acid (ABA)-mediated moisture stress response [45]. ABA is a phytohormone that adjusts the drought stress response in plants. Additional ABA application enhances the tolerance of wheat [46] and tomato [47] seedlings to drought and ABA is accumulating under drought in maize [48].
Gene MA_10265000g0010 is highly homologous to intracellular pathogenesis-related PR10 proteins. These genes produce proteins that are toxic to invading fungal pathogens and act in defense signaling. PR10 gene down-regulation was observed in maize early in drought-induced plants [49]. PR-10 protein was also noted to accumulate in the needles of drought-stressed Pinus pinaster Ait. seedlings [50]. The sequence of gene MA_10265000g0010 has a 100% query coverage and 98.97% identity to Picea glauca (Moench) Voss putative intracellular pathogenesis-related protein (picg4) [51]. Ekramoddoullah [51] showed that the expression of PR-10 (including picg4) genes is induced by both biotic (white pine blister rust, Cronartium ribicola J.C.Fischer) and abiotic (cold) stress. Therefore, we suggest that the gene MA_10265000g0010 could be a possible marker of abiotic and biotic factors (i.e., drought stress and pathogen attack) for transcriptional response studies in conifers [51] as the expression seems to accumulate in the combination of these stressors.
Gene MA_10435878g0010 plays a role in metabolism and it belongs to the 2OG-Fe(II) oxygenase superfamily. It has been found to be significantly expressed in the biosynthesis of secondary metabolites and in metabolic pathways under different light spectra in Norway spruce [52]. It appears that this gene reacts to environmental change and is accumulates more when we combine biotic stress (H. parviporum) with abiotic stress (drought). MA_189802g0010 belongs to the thaumatin family and it responds to fungal organisms and salinity; the association of this gene with P. abies was previously unknown.
In addition to the RNA-seq analysis, another aim of this study was to uncover preliminary insights into the defense responses of H. parviporum-infected Norway spruce under drought stress. LOX [31], ERF1 [28], CHIIV [33], and p/DIR32 [32] have previously been demonstrated to be activated and expressed when Norway spruce is infected with H. parviporum; thus far, our observations from LOX and ERF1 through qRT-PCR concur with these studies. On the other hand, although we had similar observations [32] from p/DIR32 showing greater expression in wounded Norway spruce than H. parviporum-infected ones, they were only induced under drought stress conditions. This similar observation was also the same in CHIIV [33]. Plausible reasons for these findings are unfortunately unknown to date. Nevertheless, we discover that this defense-related gene response is more complex when additional abiotic stress factors are present. Summarily, we deduce that water availability plays a critical role in the gene switch between up-/down-regulation strategies employed by H. parviporum-infected Norway spruce. This is evident when p/DIR32 regulation is heavily influenced by water supply. Although CHIIV is more responsive to the presence of physical wounding under normal water conditions, this response can be easily overtaken by H. parviporum infection when subjected to drought stress. This suggests a preferential defense response of Norway spruce to H. parviporum infection over physical wounding, depending on water availability conditions. Likewise, this preferential response can also be influenced by the type of damage and infection present in the plant. For example, when affected by physical wounding, CHIIV is expressed under normal water conditions while ERF1 is expressed under drought stress conditions; this is the other way around when Norway spruce is infected with H. parviporum. Through these observations, the complexity of defense-related gene response in H. parviporum-infected Norway spruce is demonstrated and piques new research questions: within the Norway spruce—H. parviporum forest pathosystem, how does the host make its choice in activating a specific defense-related gene, and how does it relate back to the various abiotic signaling pathways involved? Currently, due to the explorative nature of our study, the RNA-seq and qRT-PCR dataset we have gathered thus far is, unfortunately, not sufficiently comprehensive to make sensible inferences to the above-mentioned questions. Therefore, opportunities for further in-depth research in this direction may be prompted.

4.3. Norway Spruce Resistance in a Changing Climate

Norway spruce is not only highly economic but also possesses ecological value in Europe, particularly in Northern and Central European countries. The current breeding activities concentrate on developing more productive and climate resilient trees in Scandinavia. Similarly, several markers are suggested for Heterobasidion annosum s.l. resistance [24,40,53,54]. These traits should all be studied for suitability in breeding programs as due to long tree generations, Norway spruce might not keep up with the changes in the local environment, making them more susceptible to infection by Heterobasidion species. In addition, with the new discovery of a preferential defense gene response in Norway spruce mediated by water availability, understanding how specific defense-related genes are regulated in the midst of a changing climate will be pivotal to support the successful breeding of these trees. The expression of these gene regions would be good to be tested in multifactorial analysis. Necrosis caused by fungal pathogens is increasing under host stress [3,10,18], and revealing the molecular determinants of resistance/virulence in a changing environment of forest pathosystems would be needed to mitigate this detrimental response. The major threats to the sustainable supply of forest tree products are adverse climates that benefit the pests and diseases at the cost of the host. Changes in the environment claim for the inclusion of new traits in breeding programs. Increasing the availability of Norway spruce and H. parviporum transcriptome data thus creates new opportunities to discover genes that underlie this tree-pathogen interaction on the projected climate change continuum.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/f12101379/s1.

Author Contributions

Conceptualization, K.B. and E.T.; Data curation, X.H.-Y.Y., K.B. and E.T.; Formal analysis, X.H.-Y.Y., B.D. and E.T.; Methodology, X.H.-Y.Y. and E.T.; Project administration, E.T.; Supervision, K.B. and E.T.; Writing—original draft, X.H.-Y.Y., B.D. and E.T.; Writing—Review & editing, X.H.-Y.Y., K.B. and E.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Faculty of Forest Sciences and Forest Ecology, Georg-August-Universität Göttingen, Germany.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article, the Supplementary Material and at the NCBI SRA under Bioproject PRJNA761217.

Acknowledgments

The work for RNA-seq data analysis has been performed under the Project HPC-EUROPA3 (INFRAIA-2016-1-730897), with the support of the EC Research Innovation Action under the H2020 Programme; in particular, Terhonen gratefully acknowledges the support of Fred Asiegbu and Zhen Zeng, Department of Forest Sciences, University of Helsinki, and the computer resources and technical support provided by CSC—IT CENTER FOR SCIENCE LTD, Finland. We acknowledge support by the Open Access Publication Funds of the Göttingen University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. NASA. The Causes of Climate Change. Available online: https://climate.nasa.gov/causes (accessed on 23 June 2019).
  2. Seidl, R.; Thom, D.; Kautz, M.; Martin-Benito, D.; Peltoniemi, M.; Vacchiano, G.; Wild, J.; Ascoli, D.; Petr, M.; Honkaniemi, J.; et al. Forest Disturbances under Climate Change. Nat. Clim. Chang. 2017, 7, 395–402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Linnakoski, R.; Sugano, J.; Junttila, S.; Pulkkinen, P.; Asiegbu, F.O.; Forbes, K.M. Effects of Water Availability on a Forestry Pathosystem: Fungal Strain-Specific Variation in Disease Severity. Sci. Rep. 2017, 7, 1–7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Linnakoski, R.; Kasanen, R.; Dounavi, A.; Forbes, K.M. Editorial: Forest Health Under Climate Change: Effects on Tree Resilience, and Pest and Pathogen Dynamics. Front. Plant Sci. 2019, 10, 1157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Terhonen, E.; Blumenstein, K.; Kovalchuk, A.; Asiegbu, F.O. Forest Tree Microbiomes and Associated Fungal Endophytes: Functional Roles and Impact on Forest Health. Forests 2019, 10, 42. [Google Scholar] [CrossRef] [Green Version]
  6. Moore, D.; Robson, G.D.; Trinci, A.P. 21st Century Guidebook to Fungi; Cambridge University Press: Cambridge, UK, 2020; ISBN 1-108-80784-4. [Google Scholar]
  7. McNew, G.L. The Nature, Origin, and Evolution of Parasitism. Plant Pathol. Adv. Treatise 1960, 2, 19–69. [Google Scholar]
  8. La Porta, N.; Capretti, P.; Thomsen, I.M.; Kasanen, R.; Hietala, A.M.; Von Weissenberg, K. Forest Pathogens with Higher Damage Potential Due to Climate Change in Europe. Can. J. Plant Pathol. 2008, 30, 177–195. [Google Scholar] [CrossRef]
  9. Brodde, L.; Adamson, K.; Julio Camarero, J.; Castaño, C.; Drenkhan, R.; Lehtijärvi, A.; Luchi, N.; Migliorini, D.; Sánchez-Miranda, Á.; Stenlid, J.; et al. Diplodia Tip Blight on Its Way to the North: Drivers of Disease Emergence in Northern Europe. Front. Plant Sci. 2019, 9, 1818. [Google Scholar] [CrossRef] [Green Version]
  10. Blumenstein, K.; Bußkamp, J.; Langer, G.J.; Schlößer, R.; Parra Rojas, N.M.; Terhonen, E. Sphaeropsis Sapinea and Associated Endophytes in Scots Pine: Interactions and Effect on the Host under Variable Water Content. Front. For. Glob. Chang. 2021, 4, 55. [Google Scholar] [CrossRef]
  11. Garrett, K.A.; Dendy, S.P.; Frank, E.E.; Rouse, M.N.; Travers, S.E. Climate Change Effects on Plant Disease: Genomes to Ecosystems. Annu. Rev. Phytopathol. 2006, 44, 489–509. [Google Scholar] [CrossRef] [Green Version]
  12. Bostock, R.M.; Pye, M.F.; Roubtsova, T.V. Predisposition in Plant Disease: Exploiting the Nexus in Abiotic and Biotic Stress Perception and Response. Annu. Rev. Phytopathol. 2014, 52, 517–549. [Google Scholar] [CrossRef] [Green Version]
  13. Asiegbu, F.O.; Adomas, A.; Stenlid, J. Conifer Root and Butt Rot Caused by Heterobasidion Annosum (Fr.) Bref. Sl. Mol. Plant Pathol. 2005, 6, 395–409. [Google Scholar] [CrossRef]
  14. Redfern, D.B. Spore Dispersal and Infection. In Heterobasidion Annosum: Ecology, Impact and Control; CAB International: Wallingford, UK, 1998; pp. 105–124. [Google Scholar]
  15. Gonthier, P.; Garbelotto, M. Reducing the Threat of Emerging Infectious Diseases of Forest Trees—Mini Review. CAB Rev. 2013, 8, 1–2. [Google Scholar] [CrossRef] [Green Version]
  16. Žemaitis, P.; Stakenas, V. Ecological Factors Influencing Frequency of Norway Spruce Butt Rot in Mature Stands in Lithuania. Russ. J. Ecol. 2016, 47, 355–363. [Google Scholar] [CrossRef]
  17. Piri, T.; Selander, A.; Hantula, J.; Kuitunen, P. Juurikääpätuhojen Tunnistaminen Ja Torjunta. 2019. Available online: http://urn.fi/URN:NBN:fi-fe2019091828606 (accessed on 2 August 2021).
  18. Terhonen, E.; Langer, G.J.; Bußkamp, J.; Rӑscuţoi, D.R.; Blumenstein, K. Low Water Availability Increases Necrosis in Picea Abies after Artificial Inoculation with Fungal Root Rot Pathogens Heterobasidion Parviporum and Heterobasidion Annosum. Forests 2019, 10, 55. [Google Scholar] [CrossRef] [Green Version]
  19. Hantula, J.; Vainio, E. Specific Primers for the Differentiation of Heterobasidion Annosum (s. Str.) and H. Parviporum Infected Stumps in Northern Europe. Silva Fenn. 2003, 37, 181–187. [Google Scholar] [CrossRef] [Green Version]
  20. Kovalchuk, A.; Zeng, Z.; Ghimire, R.P.; Kivimäenpää, M.; Raffaello, T.; Liu, M.; Mukrimin, M.; Kasanen, R.; Sun, H.; Julkunen-Tiitto, R. Dual RNA-Seq Analysis Provides New Insights into Interactions between Norway Spruce and Necrotrophic Pathogen Heterobasidion Annosum Sl. BMC Plant Biol. 2019, 19, 1–17. [Google Scholar] [CrossRef]
  21. Capador-Barreto, H.D.; Bernhardsson, C.; Milesi, P.; Vos, I.; Lundén, K.; Wu, H.X.; Karlsson, B.; Ingvarsson, P.K.; Stenlid, J.; Elfstrand, M. Killing Two Enemies with One Stone? Genomics of Resistance to Two Sympatric Pathogens in Norway Spruce. Mol. Ecol. 2021, 30, 4433–4447. [Google Scholar] [CrossRef]
  22. Swedjemark, G.; Stenlid, J.; Karlsson, B. Genetic Variation among Clones of Picea Abies in Resistance to Growth of Heterobasidion Annosum. Silvae Genetica 1998, 46, 369–374. [Google Scholar]
  23. Arnerup, J.; Swedjemark, G.; Elfstrand, M.; Karlsson, B.; Stenlid, J. Variation in Growth of Heterobasidion Parviporum in a Full-Sib Family of Picea Abies. Scand. J. For. Res. 2010, 25, 106–110. [Google Scholar] [CrossRef]
  24. Elfstrand, M.; Baison, J.; Lundén, K.; Zhou, L.; Vos, I.; Capador, H.D.; Åslund, M.S.; Chen, Z.; Chaudhary, R.; Olson, Å.; et al. Association Genetics Identifies a Specifically Regulated Norway Spruce Laccase Gene,PaLAC5, Linked to Heterobasidion Parviporum Resistance. Plant Cell Environ. 2020, 43, 1779–1791. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Kovalchuk, A.; Keriö, S.; Oghenekaro, A.; Jaber, E.; Raffaello, T.; Asiegbu, F. Antimicrobial Defenses and Resistance in Forest Trees: Challenges and Perspectives in a Genomic Era. Annu. Rev. Phytopathol. 2013, 51, 221–244. [Google Scholar] [CrossRef] [PubMed]
  26. Robert-Seilaniantz, A.; Grant, M.; Jones, J.D.G. Hormone Crosstalk in Plant Disease and Defense: More than Just Jasmonate-Salicylate Antagonism. Annu. Rev. Phytopathol. 2011, 49, 317–343. [Google Scholar] [CrossRef] [PubMed]
  27. Arnerup, J.; Lind, M.; Olson, Å.; Stenlid, J.; Elfstrand, M. The Pathogenic White-Rot Fungus Heterobasidion Parviporum Triggers Non-Specific Defence Responses in the Bark of Norway Spruce. Tree Physiol. 2011, 31, 1262–1272. [Google Scholar] [CrossRef] [Green Version]
  28. Arnerup, J.; Nemesio-Gorriz, M.; Lundén, K.; Asiegbu, F.O.; Stenlid, J.; Elfstrand, M. The Primary Module in Norway Spruce Defence Signalling against H. Annosum s.l. Seems to Be Jasmonate-Mediated Signalling without Antagonism of Salicylate-Mediated Signalling. Planta 2013, 237, 1037–1045. [Google Scholar] [CrossRef] [PubMed]
  29. Lunden, K.; Danielsson, M.; Durling, M.B.; Ihrmark, K.; Gorriz, M.N.; Stenlid, J.; Asiegbu, F.O.; Elfstrand, M. Transcriptional Responses Associated with Virulence and Defence in the Interaction between Heterobasidion Annosum s. s. and Norway Spruce. PLoS ONE 2015, 10, e0131182. [Google Scholar] [CrossRef]
  30. Likar, M.; Regvar, M. Early Defence Reactions in Norway Spruce Seedlings Inoculated with the Mycorrhizal Fungus Pisolithus Tinctorius (Persoon) Coker & Couch and the Pathogen Heterobasidion Annosum (Fr.) Bref. Trees 2008, 22, 861–868. [Google Scholar]
  31. Arnerup, J. Induced Defence Responses in Picea Abies Triggered by Heterobasidion Annosum s.l.; Acta Universitatis Agriculturae Sueciae: Uppsala, Sweden, 2011; p. 56. [Google Scholar]
  32. Oliva, J.; Rommel, S.; Fossdal, C.G.; Hietala, A.M.; Nemesio-Gorriz, M.; Solheim, H.; Elfstrand, M. Transcriptional Responses of Norway Spruce (Picea Abies) Inner Sapwood against Heterobasidion Parviporum. Tree Physiol. 2015, 35, 1007–1015. [Google Scholar] [CrossRef] [Green Version]
  33. Yaqoob, N.; Yakovlev, I.A.; Krokene, P.; Kvaalen, H.; Solheim, H.; Fossdal, C.G. Defence-Related Gene Expression in Bark and Sapwood of Norway Spruce in Response to Heterobasidion Parviporum and Methyl Jasmonate. Physiol. Mol. Plant Pathol. 2012, 77, 10–16. [Google Scholar] [CrossRef]
  34. Keriö, S.; Terhonen, E.; LeBoldus, J. Safe DNA-Extraction Protocol Suitable for Studying Tree-Fungus Interactions. BIO-PROTOCOL 2020, 10, 3634. [Google Scholar] [CrossRef]
  35. Zeng, Z.; Raffaello, T.; Liu, M.-X.; Asiegbu, F.O. Co-Extraction of Genomic DNA & Total RNA from Recalcitrant Woody Tissues for next-Generation Sequencing Studies. Future Sci. OA 2018, 4, FSO309. [Google Scholar]
  36. Nystedt, B.; Street, N.R.; Wetterbom, A.; Zuccolo, A.; Lin, Y.-C.; Scofield, D.G.; Vezzi, F.; Delhomme, N.; Giacomello, S.; Alexeyenko, A. The Norway Spruce Genome Sequence and Conifer Genome Evolution. Nature 2013, 497, 579–584. [Google Scholar] [CrossRef] [Green Version]
  37. Zhang, Y.; Park, C.; Bennett, C.; Thornton, M.; Kim, D. Rapid and Accurate Alignment of Nucleotide Conversion Sequencing Reads with HISAT-3N. Genome Res. 2021, gr-275193. [Google Scholar]
  38. R Core Team R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. Available online: https://www.r-project.org/ (accessed on 23 May 2020).
  39. Love, M.I.; Huber, W.; Anders, S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [Green Version]
  40. Chaudhary, R.; Lundén, K.; Dalman, K.; Dubey, M.; Nemesio-Gorriz, M.; Karlsson, B.; Stenlid, J.; Elfstrand, M. Combining Transcriptomics and Genetic Linkage Based Information to Identify Candidate Genes Associated with Heterobasidion-Resistance in Norway Spruce. Sci. Rep. 2020, 10, 12711. [Google Scholar] [CrossRef] [PubMed]
  41. Raffaello, T.; Asiegbu, F.O. Small Secreted Proteins from the Necrotrophic Conifer Pathogen Heterobasidion Annosum Sl. (HaSSPs) Induce Cell Death in Nicotiana Benthamiana. Sci. Rep. 2017, 7, 1–11. [Google Scholar] [CrossRef] [PubMed]
  42. Livak, K.J.; Schmittgen, T.D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  43. Prospero, S.; Botella, L.; Santini, A.; Robin, C. Biological Control of Emerging Forest Diseases: How Can We Move from Dreams to Reality? For. Ecol. Manag. 2021, 496, 119377. [Google Scholar] [CrossRef]
  44. Trujillo-Moya, C.; George, J.P.; Fluch, S.; Geburek, T.; Grabner, M.; Karanitsch-Ackerl, S.; Konrad, H.; Mayer, K.; Sehr, E.M.; Wischnitzki, E.; et al. Drought Sensitivity of Norway Spruce at the Species’ Warmest Fringe: Quantitative and Molecular Analysis Reveals High Genetic Variation among and within Provenances. G3 Genes Genomes Genet. 2018, 8, 1225–1245. [Google Scholar] [CrossRef] [Green Version]
  45. Harshavardhan, V.T.; Van Son, L.; Seiler, C.; Junker, A.; Weigelt-Fischer, K.; Klukas, C.; Altmann, T.; Sreenivasulu, N.; Bäumlein, H.; Kuhlmann, M. AtRD22 and AtUSPL1, Members of the Plant-Specific BURP Domain Family Involved in Arabidopsis Thaliana Drought Tolerance. PLoS ONE 2014, 9, e110065. [Google Scholar] [CrossRef] [Green Version]
  46. Wei, L.; Wang, L.; Yang, Y.; Wang, P.; Guo, T.; Kang, G. Abscisic Acid Enhances Tolerance of Wheat Seedlings to Drought and Regulates Transcript Levels of Genes Encoding Ascorbate-Glutathione Biosynthesis. Front. Plant Sci. 2015, 6, 458. [Google Scholar] [CrossRef]
  47. Yadegari, L.Z.; Heidari, R.; Rahmani, F.; Khara, J. Drought Tolerance Induced by Foliar Application of Abscisic Acid and Sulfonamide Compounds in Tomato. J. Stress Physiol. Biochem. 2014, 10, 326–334. [Google Scholar]
  48. Wang, C.; Yang, A.; Zhang, H.Y. and J. Influence of Water Stress on Endogenous Hormone Contents and Cell Damage of Maize Seedlings. J. Integr. Plant Biol. 2008, 50, 427. [Google Scholar] [CrossRef]
  49. Fountain, J.C.; Chen, Z.-Y.; Scully, B.T.; Kemerait, R.C.; Lee, R.D.; Guo, B. Pathogenesis-Related Gene Expressions in Different Maize Genotypes under Drought Stressed Conditions. AJPS 2010, 4, 433–440. [Google Scholar] [CrossRef]
  50. Dubos, C.; Plomion, C. Drought Differentially Affects Expression of a PR-10 Protein, in Needles of Maritime Pine (Pinus Pinaster Ait.) Seedlings. J. Exp. Bot. 2001, 52, 1143–1144. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Ekramoddoullah, A.K.M. Physiology and Molecular Biology of a Family of Pathogenesis-Related PR-10 Proteins in Conifers. J. Crop. Improv. 2004, 10, 261–280. [Google Scholar] [CrossRef]
  52. OuYang, F.; Mao, J.-F.; Wang, J.; Zhang, S.; Li, Y. Transcriptome Analysis Reveals That Red and Blue Light Regulate Growth and Phytohormone Metabolism in Norway Spruce [Picea Abies (L.) Karst.]. PLoS ONE 2015, 10, e0127896. [Google Scholar] [CrossRef] [Green Version]
  53. Nemesio-Gorriz, M.; Hammerbacher, A.; Ihrmark, K.; Källman, T.; Olson, Å.; Lascoux, M.; Stenlid, J.; Gershenzon, J.; Elfstrand, M. Different Alleles of a Gene Encoding Leucoanthocyanidin Reductase (PaLAR3) Influence Resistance against the Fungus Heterobasidion Parviporum in Picea Abies. Plant Physiol. 2016, 171, 2671–2681. [Google Scholar] [CrossRef] [PubMed]
  54. Edesi, J.; Tikkinen, M.; Elfstrand, M.; Olson, Å.; Varis, S.; Egertsdotter, U.; Aronen, T. Root Rot Resistance Locus PaLAR3 Is Delivered by Somatic Embryogenesis (SE) Pipeline in Norway Spruce (Picea Abies (L.) Karst.). Forests 2021, 12, 193. [Google Scholar] [CrossRef]
Figure 1. The gene expression levels of the four genes ((a) MA_189802g0010, (b) MA_10265000g0010, (c) MA_10435878g0010 and (d) MA_178006g0010) found exclusively to be significantly expressed only under drought (compared to (+) group non-treated).
Figure 1. The gene expression levels of the four genes ((a) MA_189802g0010, (b) MA_10265000g0010, (c) MA_10435878g0010 and (d) MA_178006g0010) found exclusively to be significantly expressed only under drought (compared to (+) group non-treated).
Forests 12 01379 g001
Figure 2. Heterobasidion resistance related genes and their expressions in each treatment (compared to (+) group non-treated). (a) MA_110169g0010 and (b) MA_14707g0010 were downregulated while (c) MA_15852g0010 and (d) MA_10427673g0020 were upregulated.
Figure 2. Heterobasidion resistance related genes and their expressions in each treatment (compared to (+) group non-treated). (a) MA_110169g0010 and (b) MA_14707g0010 were downregulated while (c) MA_15852g0010 and (d) MA_10427673g0020 were upregulated.
Forests 12 01379 g002
Figure 3. Palac5 (a) MA_97119g0010 and (b) MA_97119g0020 gene expression levels in each treatment.
Figure 3. Palac5 (a) MA_97119g0010 and (b) MA_97119g0020 gene expression levels in each treatment.
Forests 12 01379 g003
Figure 4. Overview of qRT-PCR mean fold change for each tested gene against ELF1. Variation between triplicate samples of each treatment group was represented by error bars (+/− SD of each triplicate) respectively.
Figure 4. Overview of qRT-PCR mean fold change for each tested gene against ELF1. Variation between triplicate samples of each treatment group was represented by error bars (+/− SD of each triplicate) respectively.
Forests 12 01379 g004
Figure 5. Influence of H. parviporum infection and physical wounding as singular and combinatorial factors in Norway spruce defense responses under drought stress conditions, as reflected by gene regulation in (a) LOX, (b) ERF1, (c) CHIIV, and (d) p/DIR32.
Figure 5. Influence of H. parviporum infection and physical wounding as singular and combinatorial factors in Norway spruce defense responses under drought stress conditions, as reflected by gene regulation in (a) LOX, (b) ERF1, (c) CHIIV, and (d) p/DIR32.
Forests 12 01379 g005
Table 1. Statistically Differentially Expressed Genes of Each Treatment When Compared to Non-treated (+) Water Group.
Table 1. Statistically Differentially Expressed Genes of Each Treatment When Compared to Non-treated (+) Water Group.
InoculationWaterSignificantUp-RegulatedDown-RegulatedNot SignificantSUM
Heterobasidion parviporum52113338829,59230,113
Heterobasidion parviporum+1841414329,94430,128
Non-treated82275529,71529,797
Mock-inoculated78116729,95030,028
Mock-inoculated+2021485429,95130,153
Table 2. Four Drought-related Candidate Genes, their Description and GO ID.
Table 2. Four Drought-related Candidate Genes, their Description and GO ID.
GeneDescriptionChromosomeSourceConfidenceTrinityGOGO-Description
MA_10265000g0010PR10MA_10265000AUGUSTUSHighcomp87387_c0_seq1
MA_10435878g0010ACCMA_10435878AUGUSTUSHighcomp94559_c0_seq1GO:0008152GO:0008152-metabolic process
MA_178006g0010BURP domain RD22-likeMA_178006AUGUSTUSHighcomp75759_c0_seq1
MA_189802g0010ThaumatinMA_189802AUGUSTUSHighcomp87856_c1_seq1
Table 3. Mean Fold Change of Each Candidate Gene in qRT-PCR for Each Treatment Scheme.
Table 3. Mean Fold Change of Each Candidate Gene in qRT-PCR for Each Treatment Scheme.
GeneMean Fold Change (2−ΔΔCt)SDInoculation TypeWater Treatment
LOX1.0060.129H. parviporumnormal
1.0880.291non-treatednormal
0.9620.206mock-controlnormal
0.9900.024H. parviporumlow
1.0520.122non-treatedlow
1.0380.102mock-controllow
ERF11.0090.162H. parviporumnormal
1.0170.210non-treatednormal
1.0040.069mock-controlnormal
1.0090.145H. parviporumlow
1.0020.073non-treatedlow
1.1120.565mock-controllow
CHIIV0.9980.020H. parviporumnormal
1.0330.267non-treatednormal
1.0200.251mock-controlnormal
1.0010.098H. parviporumlow
1.0190.214non-treatedlow
1.0280.171mock-controllow
p/DIR321.1820.788H. parviporumnormal
1.0340.324non-treatednormal
1.0250.288mock-controlnormal
1.0900.311H. parviporumlow
1.0130.139non-treatedlow
1.1720.822mock-controllow
Table 4. Comparative Quantification Determining Mean Fold Change (2−ΔΔCt) in Gene Expression Mediated by H. parviporum Infection only, Physical Wounding only, and the Combination of both H. parviporum Infection and Physical Wounding Caused by Inoculation.
Table 4. Comparative Quantification Determining Mean Fold Change (2−ΔΔCt) in Gene Expression Mediated by H. parviporum Infection only, Physical Wounding only, and the Combination of both H. parviporum Infection and Physical Wounding Caused by Inoculation.
GeneWater Treatment2−ΔΔCt
H. parviporumPhysical WoundH. parviporum + Physical Wound
LOXnormal1.600.170.27
low1.276.928.82
ERF1normal2.810.310.87
low0.462.311.06
CHIIVnormal0.1910.852.03
low7.160.916.50
p/DIR32normal0.250.330.08
low2.6417.0344.94
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yeoh, X.H.-Y.; Durodola, B.; Blumenstein, K.; Terhonen, E. Drought Stress Described by Transcriptional Responses of Picea abies (L.) H. Karst. under Pathogen Heterobasidion parviporum Attack. Forests 2021, 12, 1379. https://doi.org/10.3390/f12101379

AMA Style

Yeoh XH-Y, Durodola B, Blumenstein K, Terhonen E. Drought Stress Described by Transcriptional Responses of Picea abies (L.) H. Karst. under Pathogen Heterobasidion parviporum Attack. Forests. 2021; 12(10):1379. https://doi.org/10.3390/f12101379

Chicago/Turabian Style

Yeoh, Xenia Hao-Yi, Blessing Durodola, Kathrin Blumenstein, and Eeva Terhonen. 2021. "Drought Stress Described by Transcriptional Responses of Picea abies (L.) H. Karst. under Pathogen Heterobasidion parviporum Attack" Forests 12, no. 10: 1379. https://doi.org/10.3390/f12101379

APA Style

Yeoh, X. H. -Y., Durodola, B., Blumenstein, K., & Terhonen, E. (2021). Drought Stress Described by Transcriptional Responses of Picea abies (L.) H. Karst. under Pathogen Heterobasidion parviporum Attack. Forests, 12(10), 1379. https://doi.org/10.3390/f12101379

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop