Evolution of pathogen response genes associated with increased disease susceptibility during adaptation to an extreme drought in a Brassica rapa plant population
BMC Ecology and Evolution volume 21, Article number: 61 (2021)
Pathogens are key components in natural and agricultural plant systems. There is evidence of evolutionary changes in disease susceptibility as a consequence of climate change, but we know little about the underlying genetic basis of this evolution. To address this, we took advantage of a historical seed collection of a Brassica rapa population, which we previously demonstrated evolved an increase in disease susceptibility to a necrotrophic fungal pathogen following a drought.
Previously, we combined a resurrection experiment with genome-wide sequencing of 124 pooled ancestral and descendant plants. Here, using these previously generated sequence data (Franks et al. in Mol Ecol 25(15):3622–3631, 2016), we show that well-characterized necrotrophic fungal pathogen response (NFPR) genes have evolved, as indicated by changes in allele frequency, between ancestors and descendants, with several of them identified as extreme FST outliers. The jasmonic acid (JA) signaling pathway in particular seems to underlie the evolution of disease susceptibility, in addition to its well characterized role in plastic disease response. We identify a list of 260 genes that are both NFPR genes and are differentially expressed in response to drought, based on publicly available data. We present evidence that five of these genes evolved between ancestors and descendants, suggesting that the drought acted as the evolutionary driver, and that the accompanying increase in disease susceptibility may have been a consequence of genetic pleiotropy.
Our study provides evidence that for this population, standing variation in NFPR genes is affected by natural selection related to climate change. Our results reveal potentially important candidates that may underlie trait evolution in both crops and natural systems. Additionally, this trade-off between adaptation to biotic and abiotic stresses is an example of how climate change can have diverse and unexpected consequences.
Pathogens exert strong selection on plants, often causing them to evolve increased resistance. Simultaneously, plant systems are exposed to a myriad of other selective pressures. Persistent climate change, such as long-term drought, may result in evolution of the hosts, pathogens, or both, leading to changes in virulence or susceptibility [1,2,3,4,5]. More broadly, exploration of how multiple, and at times opposing, selective pressures shape the genomes of plants is an important area of ongoing research.
The genetics of plant pathogen response is exceedingly complex and specific to whether the pathogen is necrotrophic or biotrophic [6, 7]. While response to biotrophs has been extensively studied and found to be generally mediated by a salicylic acid-dependent pathway, response to necrotrophic pathogens (the focus here) is less well understood and tends to have a more complex genetic architecture. Necrotrophic pathogens, which first kill tissue and then extract nutrients , stimulate a response pathway generally mediated by jasmonic acid (JA) and ethylene (ET) . Such biotic stress response pathways have been shown to interact with abiotic stress response pathways. In particular, abscisic acid (ABA), which regulates growth and development in response to water availability, has been shown to suppress biotic stress response pathways, leaving plants more susceptible to pathogens . ABA is thus an example of a regulatory system that may potentially be involved in adaptation to drought , and which may also result in pleiotropic evolution of pathogen susceptibility.
The most common approach used to characterize the genes responsible for pathogen response are expression studies in model plant-pathogen systems. However, low genetic variation, and distant phylogenetic relationships often limit the translation of such studies to crop breeding or characterization of evolving natural systems. Therefore, there is great promise in studying plant-pathogen genetics of natural systems. Some progress has been made in plant-pathogen systems with a simple genetic architecture including exploring genetic variation in resistance in natural populations for the Linum-Melampsora system , as well as exploration of resistance gene diversity within and between wild tomato populations [13, 14], and variation in specific resistance genes such as the RPW8 (Recognition of Powdery Mildew) . But mostly, thus far, variation in resistance has been mapped at a spatial scale looking at local adaptation, rather than being followed over time [3, 12, 16]. The lack of focus on evolutionary trajectories limits our ability to predict evolutionary responses over short time scales. We addressed this limitation by taking advantage of evolutionary events, measuring changes in allele frequency, which allowed us to compare the genomes of ancestors and descendants to identify the genes involved in evolutionary changes in disease susceptibility in a natural system.
The system we used is the well-studied example of the adaptive evolution of earlier flowering in the annual herbaceous plant, Brassica rapa L. (syn. campestris) (Brassicaceae, field mustard) in response to a natural, extended drought [17,18,19]. Previous work used a resurrection approach [20, 21] in which B. rapa seed was collected in southern California in 1997, before an extended drought, and then again in 2004, post-drought. During the drought period, precipitation was substantially below-average, particularly later in spring, resulting in an abbreviated growing season . Pre- and post-drought seeds were then grown under common conditions in the greenhouse and challenged with Alternaria brassicae. A. brassicae, which causes Alternaria blackspot disease, is a necrotrophic pathogenic sac fungus that causes damping off, leaf spots, defoliation and reduced seed yield . We previously found that this population experienced a genetically-based evolutionary increase in Alternaria blackspot disease susceptibility , together with an adaptive shift to escape drought through earlier flowering . While the shift to earlier flowering and drought escape was adaptive, the drivers for an increase in disease susceptibility remain unclear. Pathogen susceptibility may have evolved either due to direct selection, or indirectly as a result of selection on common genes [4, 23,24,25].
To explore the shift in disease susceptibility following a drought, we hypothesized that: (1) known necrotrophic fungal pathogen response genes (including on JA and ET pathways) evolved, and (2) pleiotropic genes, underlying pathogen and drought response evolved. To test these hypotheses we used genome sequences, generated in our previous work , of an ancestral (pre-drought) and descendant (post-drought) B. rapa natural population, and explored whether necrotrophic fungal pathogen response (NFPR) genes evolved. We assessed signatures of evolution in two complementary ways. First, we conducted genome-wide outlier FST analysis to determine if allele frequencies of previously characterized NFPR genes were differentiated between ancestors and descendants, indicating rapid contemporary evolution at these loci. Second, we used site frequency spectrum analysis (Tajima’s D) to look for historical signatures of selective sweeps or other strong or recurrent selective patterns at NFPR genes. Finally, we identified a list of genes that were both NFPR genes and were differentially expressed in response to drought, based on publicly available data. We used this list to consider the potential role of pleiotropic evolution of both NFPR and drought response.
Previously we identified 5,812,602 SNPs (out of 235,128,010 sites), including SNPs in 35,202 genes using a low-coverage approach (~ 25× per pooled library). The general results of our sequencing and our observation of global and flowering time related allele frequency shifts have been published elsewhere . In that study, we found an expected heterozygosity of 0.5–0.6 across the pools, with no major reductions in genetic diversity or heterozygosity, indicating that the drought likely did not cause a genetic bottleneck. Also, we observed an evolutionary shift in disease susceptibility for this population following the drought event . In the present study, we re-analyzed our genomics data for candidate NFPR genes in order to explore the genetic basis of this evolutionary shift in disease susceptibility. We found evidence of evolution in NFPR genes, demonstrated by a change in allele frequencies at these pathogen response loci between ancestor and descendant populations, and that the JA signaling pathway seemed to underlie the evolution of disease susceptibility. We found that a number of these NFPR genes were also responsive to drought, suggesting the potential role of pleiotropic evolution.
Evolution of NFPR and NFPR/drought response genes
To explore our first hypothesis, that pathogen response genes evolved, we used an a priori approach, and found that 11 genes from of our database of 1182 candidate NFPR genes were significantly differentiated between ancestors and descendants (Fig. 1), with an FST above background compared to all the genes, as determined by outlier FST analysis (q < 0.05). These genes were found throughout the genome on the majority of B. rapa’s 10 chromosomes, including chromosomes 1, 2, 4, 5, 6, 7, and 9. In contrast, in the entire genome, we saw 434 out of 35,202 genes evolved . The NFPR genes that evolved were not over-represented among all genes that showed significant evolution (Chi-square = 0.846, p = 0.358). The mean FST for all NFPR genes was not significantly elevated above background (ANOVA: F1,35199 = 0.008, p = 0.929). In a Manhattan plot comparing FST of NFPR to all other genes across the genome, the NFPR and all genes’ trend lines tracked one another (Fig. 1), in fact 95% confidence intervals (not shown) overlapped across the whole genome. We report the 20 most differentiated NFPR genes (Table 1; see Additional file 1 for top 50). This includes the 11 genes that were significant, plus additional genes that could potentially have been involved in the shift in susceptibility despite not being statistically significant after false discovery correction.
To explore our second hypothesis, that genes involved in both NFPR and drought response evolved, we also used an a priori approach and found that 5 (2.0%) of the genes from our joint list of 260 genes that responded to both conditions were significantly differentiated between ancestors and descendants, with an FST above background compared to all genes (q < 0.05) (Table 1). These five genes did not constitute an over-representation among all genes that showed significant evolution (Chi-square = 0.974, p = 0.324). We show the overlap in NFPR, drought response, and genes that evolved (demonstrated by significant FST values) in a Venn diagram (Fig. 2). This list of genes that evolved, and that are associated with both NFPR and drought response, are candidates genes for pleiotropic evolution (i.e. where genetic change at a single locus affected both phenotypes).
To determine if the genes that were highly differentiated between 1997 and 2004 were also under selection historically, we calculated the Tajima’s D value for each gene. We found that none of the highly differentiated genes also had a significantly reduced Tajima's D value. There was no difference in Tajima’s D between significantly differentiated (FST) NFPR genes and the rest of the NFPR genes (ANOVA: ancestors F1,1180 = 2.398, p = 0.122; descendants F1,1180 = 0.551, p = 0.458). This was consistent with our genome-wide result, where we found little overlap between genes with high FST and negative Tajima’s D . These findings were not surprising because these estimates detect selection acting at different time scales. Visual inspection of Tajima's D values over the genome suggested that three of our significant (FST) NFPR genes were in regions of reduced Tajima’s D (red bars in Fig. 1), which could indicate they were linked to regions that may have experienced historical selective sweeps, although the values of Tajima’s D for the genes of interest were themselves not significantly reduced.
We next explored the nature of segregating sites in two NFPR candidate genes with the highest FST values, Bra030295 and Bra038107. We explored whether polymorphism in these genes would result in changes in the amino acid sequence, and therefore potentially be a direct subject of selection. We made alignments of the ancestral, descendant, and reference sequences for these genes (Fig. 3). We determined whether sites were in exons or introns and, if in exons, whether they were synonymous or nonsynonymous. These genes were also in regions that had reduced Tajima’s D, by visual inspection, or have important functions (see “Discussion”). We found that the two genes together had 10 SNPs in exons, and 5 of these exonic SNPs were nonsynonymous. Bra030295 had two and Bra038107 had three nonsynonymous sites that experienced an allele shift between ancestors to descendants.
Jasmonic acid and the evolution of disease susceptibility
Since the JA pathway is such an important mediator of necrotrophic fungal pathogen response, we specifically explored whether JA (and downstream gene COI1) may be involved in the evolution of disease susceptibility. We found that 17 (34%) out of our top 50 most NFPR differentiated genes were COI1 and/or JA dependent (Additional file 1). This is similar to the ratio 441 (37%) out of 1182 of COI1/JA genes which are differentially expressed in response to pathogens in van Wees et al.  and Mukherjee et al.  pathogen studies; we found that these ratios are statistically equivalent (Chi-square = 0.225, p = 0.635), despite the very different approach of analysis using resurrection/sequencing instead of gene expression.
Gene ontology analysis
We used gene ontology analysis to determine if pathogen response functional categories are overrepresented globally among all genes when ranked by their level of evolutionary differentiation between ancestors and descendants (FST). We found that pathogen response categories were not significantly enriched and were not highly ranked out of the 3136 categories. The top pathogen response category was “plant-type hypersensitive response” (GO:0009626), ranked 54th (p = 0.014, Benjamini–Hochberg FDR corrected p = 0.835), followed closely by “host programmed cell death induced by a symbiont” (GO:0034050), ranked 56th (p = 0.015, Benjamini–Hochberg FDR corrected p = 0.822).
In previous work, we found that a natural B. rapa plant population in southern California evolved earlier flowering  and increased fungal disease susceptibility  following a 7-year natural drought. A shift in susceptibility for this plant pathogen system could have widespread implications because this fungus is an important agricultural problem for B. rapa crops [29, 30]. Also, we have found widespread infection (22.2% of plants observed to be infected) in all 6 of our B. rapa field populations on the west coast of the USA . In the current study, we compared the genomics of ancestral and descendant plants from one population to elucidate the underlying genetic and molecular basis of this shift in disease susceptibility. This is the first study, to our knowledge, to examine the genetic basis of evolution of susceptibility to pathogens in a natural plant population following drought [32,33,34,35]. There are many challenges to studying evolution in natural systems. Using a resurrection genomics approach we were able to both directly measure the disease susceptibility of ancestors and descendants in common conditions, and directly assess evolution by measuring changes in allele frequencies between these ancestral and descendant populations .
Evolution of NFPR genes
Using an outlier FST analysis approach, we found evidence for our first hypothesis, that a number of NFPR genes evolved over the course of just 7 years of drought. This is consistent with selection acting on a trait that is controlled by multiple loci, such as defense response to a necrotrophic fungus. This is expected to be more complex than response to a biotrophic pathogen, which in many cases is controlled by a single gene [6, 36]. We did not see significant evolution across the wider NFPR system compared to the rest of the genome. It is possible, given the genetically complex nature of NFPR, that specific causative genes could evolve, underlying a shift in susceptibility, without evolution across the wider system, especially given the short time frame between ancestors and descendants studied. In addition, we found that a substantial proportion of these NFPR genes are JA/COI1 dependent. Studies in the closely related A. thaliana–A. brassicicola model plant pathogen system have identified COI1-dependent genes, part of the JA signalling pathway, as playing a central role in defense against necrotrophic fungal pathogens, suggesting that JA signalling mediates the defense response. JA activity is dependent on the downstream COI1 gene, which encodes an F-box protein involved in proteolysis . JA signalling might play an important role in the evolution of susceptibility in addition to its known role in a plastic response to necrotrophic fungal pathogens.
There are limitations to using an outlier FST analysis to characterize the ultimate molecular or developmental mechanisms by which a population phenotypically evolves (e.g. driving mechanism of evolution). Additionally, our pooling sample size may limit us to identifying alleles with larger changes in allele frequency. However, our outlier FST approach, together with a q-value threshold correction, is effective at characterizing loci that showed significant shifts, and thereby can identify putative candidate genes underlying genetically-based phenotypic shifts observed. This approach has also been successfully utilized in other studies in these populations .
Potential signatures of pleiotropic evolution
In previous work, we found adaptive evolution in response to drought (through earlier flowering) also resulted in the evolution of pathogen susceptibility in this population . In exploration of this associated response, in the current study, we also found evidence for our second hypothesis: that a number of genes involved in both drought and pathogen response evolved, suggesting a signature of pleiotropic evolution. In the current study, we considered genes to be candidates for pleiotropic evolution for drought and necrotrophic fungal pathogens if they had been found to be differentially expressed in response to independent drought stress and inoculation with a necrotrophic fungal pathogen, and if in our study they evolved between ancestors and descendants. A subset of five of these genes evolved in our study. Allelic variation in any of these five genes may have simultaneously affected both drought response and NFPR response, i.e. caused pleiotropic phenotypic changes. For this population we found clear evidence that JA pathway genes, which affect both pathogen response, flowering time, and stress resistance [38,39,40], evolved, suggesting evidence for evolution in pleiotropic genes in this system. These findings agree with other studies which have found that pleiotropy plays a role in shaping multiple ecological traits for plant populations displaying adaptation to local water availability . Furthermore, it is possible that antagonistic pleiotropy was playing a role in the evolution observed . Under this scenario, pleiotropic genes underlying both drought and pathogen response could have evolved under positive selection due to drought, despite negative selection by pathogens, if the negative selection was not strong enough to overcome the positive selection, or if it was not acting in the contemporary environment . Indeed, the increased susceptibility we observed in this population was accompanied by an evolutionary shift to thinner leaves (increased specific leave area) , which could indicate a trade-off between drought escape and disease susceptibility in agreement with the growth and defense trade-off [43, 44].
In addition to examining evolution following the recent drought, we tested whether NFPR genes were under historical directional selection by looking for genes with a low Tajima's D (< − 2). Loci under balancing selection would not have a low Tajima's D. Genome-wide, we found no association between FST and Tajima’s D, and there were no genes that were outliers for both high FST (> 0.2) and low Tajima’s D (< − 2). This finding concurs with other prior work in this system . This indicates that evolution in response to the drought affected different genes in this population than have historically been the target of directional selection.
Top candidate genes underlying evolutionary shift in susceptibility
The most highly differentiated gene in this study was Sec14p-like phosphatidylinositol transfer family protein (Bra030295) (Fig. 3a; FST = 0.21; q = 1.75 × 10–5) which has been found to be differentially expressed in Arabidopsis in response to infection with A. brassicicola  and cabbage leaf curl virus . Little is known about the specific function of this gene, although it is thought to play a role in the movement of substances across cell membranes and has been found to be expressed widely in 23 plant structures . The shift in allele frequency at two non-synonymous segregating sites in this gene indicate a potential functional effect, with one of these polymorphisms being shared between ancestors and descendants (potentially indicating a previous shift due to the recurrent droughts in this region), and the other polymorphism found only in the descendant population, indicating a better functional candidate in response to this drought event.
Several of the NFPR genes with the highest level of differentiation between ancestors and descendants demonstrate the pleiotropic nature of the genes subject to evolution in response to the drought in this population. We found that the second most highly differentiated NFPR gene was the COI1 dependent gene ERF4 (Bra038107) (FST = 0.20; q = 7.30 × 10–5). In addition to a demonstrated role in the antagonistic regulation of JA and ET, this gene is also induced by ABA [47, 48]. This is interesting because this gene, along with other genes mediated by ABA that we found to have high FST, including Bra012746 and Bra033968, and Bra011282, have pleiotropic roles in both abiotic and biotic stress response [49,50,51]. Studies have shown that ERF4 is induced by ABA, and ERF4 suppresses expression of PDF1.2 , a defense effector that is elicited by JA signaling . This gene showed a shift in allele frequency at three non-synonymous sites, and evolution at these sites could have a functional effect (Fig. 3b). Two of these polymorphisms were shared between ancestors and descendants (potentially indicating a previous shift due to the recurrent droughts in this region), and the third polymorphism was found only in the descendant population, indicating a better functional candidate in response to this drought event.
Caffeoyl-CoA 3-O-methyltransferase (Bra033968) is an NFPR gene that is COI1 dependent  and also had a significant FST value. It is a lignin biosynthesis gene, which is pleiotropic, playing a role in both drought and pathogen response. Lignification is a primary defense mechanism against pathogens , which works primarily through leaf structure . It has also been shown to be suppressed in Arabidopsis by ABA application, along with other pathogen response genes, resulting in a decrease in lignin accumulation and an increase in disease susceptibility to Pseudomonas syringae pathovar (pv.) tomato . This gene is also interesting because lignin plays a key role in leaf structure, and we found evidence that leaf structure played a role in the increase in disease susceptibility seen , with earlier flowering plants having thinner leaves that displayed greater disease severity.
Limitations of this study
This study was limited due to it being an observational study focused on one natural population that underwent an evolutionary shift in response to an extended drought. Alternative explanations to pleiotropy include that drought response and pathogen susceptibility were acted on separately, or that if there is a large cost of disease resistance and these pathogens become less common during a drought, then selection could favor reduced resistance due to these costs. The current study design does not allow us to differentiate between these alternative explanations. Experimental evolution studies with additional populations, including a non-drought control population will be important to address some of these findings further. Additionally, epigenetics, especially longer-term stable epigenetics could have played a role in the shift observed, and was not tested in the present study. Such epigenetic effects may act in addition to or in concert with the observed genetic changes, and would not replace or negate the contribution of genetic change. Although seeds from prior generations emerging from the seedbank could have influenced our findings, this effect is likely minimal and would only have caused our analyses to be more conservative because less evolution would be observed than actually occurred . Finally, in addition to pleiotropy, correlated change in these two phenotypes could have been caused by drift, linkage, or by simultaneous change in selective pressure on both flowering time and pathogen load. While we cannot eliminate these from consideration, they likely played a minor role in the observed patterns due to data presented in a previous study showing an adaptive evolution in response to drought, lack of evidence of a bottleneck, and the short 7-year time frame over which the evolution occurred, too short for drift alone to explain the extensive allelic differentiation observed .
The results from our resurrection genomics study demonstrate that many well-characterized candidate NFPR genes evolved, with a change in allele frequency, in just 7 years of drought. We conclude that COI1 and the JA pathway are likely to be involved in the evolution of disease susceptibility in this population. These genes are thus candidates for future investigation of variation in disease susceptibility in other crop and wild populations of this and related species. Several implicated genes are potentially involved in both defense and drought resistance, suggesting that antagonistic pleiotropy potentially caused the observed evolutionary increase in disease susceptibility following a natural drought. That this previously unexplained increase in disease susceptibility is a result of adaptation to drought represents a surprising response of a plant population facing changes in local climate, and demonstrates that tradeoffs among fitness components can explain seemingly maladaptive responses.
In prior work, a large number of seeds (> 10,000) were bulk collected from a natural feral B. rapa population (Arboretum) in Southern California near the University of California Irvine campus. Seeds were collected in 1997, before an extended natural drought, and then again in 2004, post-drought. Before generating data, plants were previously grown from seed and cross-pollinated to control for maternal and storage effects .
Sample preparation and shotgun sequencing
The pooled sequence data used in this study was generated in our previous work and described in more detail there in Franks et al. . To summarize, previously, 500 ancestral and descendant plants (total N = 1000) were grown up from F2 seeds, tissue was collected, genomic DNA was extracted, quality-checked, and quantified by real-time PCR. Samples were pooled in equimolar amounts, in duplicate, within year (ancestral N = 50, descendant N = 74), and sequenced on a HiSeq 2000 (Illumina). Note that Franks et al.  also included a second population (Back Bay) not studied here, and for which we do not have disease susceptibility data. We have focused on the Arboretum population for disease susceptibility studies  due to the greater evolutionary shifts in multiple traits experienced by that population .
Bioinformatic analysis and validation
We previously used computational methods described in Franks et al.  to trim, clean and prepare our reads for analysis including aligning sequence reads to the B. rapa genome version 1.18 (http://brassicadb.org/brad/) , calling sequence variants, and calculating population genetic estimates including divergence (FST) [57, 58] and Tajima's D  across 100-kb windows. In this previous work, steps were taken to improve confidence in calling SNPs and population genetic estimates, including using a minimum of 8× coverage across each window, and removal of singletons SNPs. Pooling design, including choosing number of samples to pool to generate sufficient coverage, was informed by modelling and empirical studies that show precision and accuracy of population parameters (e.g. allele frequency) to be maximized when a low coverage (i.e. 0.5×–1×) population-level (> 20 individuals) approach is used . Population genetic estimates were compared between pooling replicates and found to be highly correlated, so these technical replicates were merged for all reported analyses. SNPs for these pooled whole-genome sequences were also validated in our previous study using Kompetitive Allele-Specific PCR genotyping (KASP) on SNPs in 10 genes for 116 individual plants, with allele frequencies obtained by KAPS-genotyping on individuals showing high correlation with pooled Illumina sequencing .
In our previous work, statistical significance of genes, results which are used here, were determined for the FST values of genes by comparing the experimental FST distribution to a null distribution created by resampling and regressing the experimental distribution on a null distribution to calculate p values for residuals. After multiple test correction, genes were considered significant if they had a q value of less than 0.05 . We also report the top 50 most differentiated NFPR genes, including those with a q value of > 0.05.
We considered previously published observations from the compatible plant-pathogen system Brassica rapa–Alternara brassicae in this study, as described in detail in our previous work . The molecular mechanisms of pathogen response in the Brassica rapa–Alternaria brassicae system is not well understood , however, extensive research has been conducted on the closely related plant pathogen system Arabidopsis thaliana–Alternaria brassicicola, which is a model system used for studying diseases caused by necrotrophs [6, 27, 28]. In the current study, we took advantage of this extensive model system work and used an a priori approach  to characterize the evolution of a set of well-researched candidate disease-related genes in our resurrected ancestral and descendent natural populations, as described in more detail in the next section. Studies in the A. thaliana–A. brassicicola system have identified a large number (~ 1100) of candidate NFPR genes [6, 27, 28]. Also, these studies found that JA is required for resistance, indicated by susceptibility in mutants with non-functional copies of genes involved in the JA pathway, including a mutation in the downstream COI1 gene (see Glazebrook  for further details). Additionally, a recent B. rapa expression study was used to identify candidate drought response genes .
In the current study, for our a priori analysis approach, a database of 1182 well-characterized candidate NFPR genes (Additional file 2) was assembled from two previous A. thaliana–A. brassicicola studies by Mukherjee et al.  and van Wees et al. ; results from these studies were used for functional annotations in our study. These studies were chosen in order to span the range of compatibility that exists between B. rapa–A. brassicae , including genes shown to be involved in both compatible interactions  and incompatible interactions . Briefly, to identify candidate NFPR genes, van Wees et al.  used mutants with non-functional copies of genes involved in JA and SA signaling followed by inoculation and expression profiling, and Mukherjee et al.  used suppression subtractive hybridization (SSH) to identify differential expression between compatible and incompatible plant lines challenged with the pathogen. For both studies, which measured differential expression across the whole transcriptome, NFPR candidate genes were determined as genes with differential expression between plants that were inoculated with the fungal pathogen and those that received a mock inoculation. We used this database for functional annotations in our study.
We next considered 7076 genes candidate drought response genes (Additional file 3) which were differentially expressed in Guo et al.  study in response to a polyethylene glycol (PEG) treatment, as compared to controls (no PEG), and that had a > twofold change at 12 h post-treatment. We then identified genes differentially expressed in common among both this drought study and the NFPR studies (above). In total, we found 260 genes were significantly differentially expressed in both types of conditions.
Lastly, we conducted a functional analysis by analyzing gene ontology (GO) using ErmineJ , with annotations from Arabidopsis thaliana. To analyze our FST distribution, which was non-normal, we used ranked gene scores in a receiver operator characteristic (ROC) analysis and default ErmineJ settings. p values were corrected using a Bejamini–Hochberg false discover rate (FDR) correction.
Availability of data and materials
The dataset supporting the conclusions of this article are included within the article (and its additional files).
Necrotrophic fungal pathogen response
Single nucleotide polymorphism
Suppression subtractive hybridization
- PEG treatment:
Polyethylene glycol treatment
- GO Analysis:
Gene ontology analysis
False discover rate
Desprez-Loustau ML, Marcais B, Nageleisen L, Piou D, Vannini A. Interactive effects of drought and pathogens in forest trees. Ann For Sci. 2006;63:597–612.
Garrett K, Thomas-Sharma S, Forbes G, Nopsa J. Climate change and plant pathogen invasions. In: Ziska L, Dukes J, editors. Invasive species and global climate change. Wallingford: CABI; 2014.
Pautasso M, Doring T, Garbelotto M, Lorenzo P, Jeger M. Impacts of climate change on plant disease—opinions and trends. Eur J Plant Pathol. 2012;133:295–313.
O’Hara N, Rest J, Franks S. Increased susceptibility to fungal disease accompanies adaptation to drought in Brassica rapa. Evolution. 2016;70(1):241–8.
Garrett KA, Dendy SP, Frank EE, Rouse MN, Travers SE. Climate change effects on plant disease: genomes to ecosystems. Annu Rev Phytopathol. 2006;44:489–509.
Glazebrook J. Contrasting mechanisms of defense against biotrophic and necrotrophic pathogens. Annu Rev Phytopathol. 2005;43:205–27.
Mengiste T. Plant immunity to necrotrophs. Annu Rev Phytopathol. 2012;50:267–94.
Agrios GN. Plant pathology. 5th ed. Burlington: Academic Press; 2005.
McDowell J, Dangl J. Signal transduction in the plant immune response. Trends Biochem Sci. 2000;25:79–82.
Mauch-Mani B, Mauch F. The role of abscisic acid in plant-pathogen interactions. Curr Opin Plant Biol. 2005;8:409–14.
Shinozaki K, Yamaguchi-Shinozaki K, Seki M. Regulatory network of gene expression in the drought and cold stress responses. Curr Opin Plant Biol. 2003;6:410–7.
Thrall PH, Burdon JJ. Evolution of virulence in a plant host-pathogen metapopulation. Science. 2003;299:1735–7.
Stam R, Scheikl D, Tellier A. Pooled enrichment sequencing identifies diversity and evolutionary pressures at NLR resistance genes within a wild tomato population. Genome Biol Evol. 2016;8(5):1501–15.
Stam R, Silva-Arias G, Tellier A. Subsets of NLR genes show differential signatures of adaptation during colonization of new habitats. New Phytol. 2019;224:367–79.
Jorgensen T, Emerson B. Functional variation in disease resistance gene in populations of Arabidopsis thaliana. Mol Ecol. 2008;17(22):4912–23.
Laine A, Burdon J, Dodds P, Thrall P. Spatial variation in disease resistance: from molecules to metapopulations. J Ecol. 2011;99(1):96–112.
Franks SJ, Sheina S, Weis AE. Rapid evolution of flowering time by an annual plant in response to a climate fluctuation. Proc Natl Acad Sci. 2007;104(4):1278–82.
Franks SJ, Weis AE. A change in climate causes rapid evolution of multiple life-history traits and their interactions in an annual plant. J Evol Biol. 2008;21:1321–34.
Franks SJ. Plasticity and evolution in drought avoidance and escape in the annual plant Brassica rapa. New Phytol. 2011;190(1):249–57.
Franks SJ, Avise JC, Bradshaw WE, Conner JK, Etterson JR, Mazer SJ, Shaw RG, Weis AE. The resurrection initiative: storing ancestral genotypes to capture evolution in action. Bioscience. 2008;58(9):870–3.
Franks S, Hamann E, Weis AE. Using the resurrection approach to understand contemporary evolution in changing environments. Evol Appl. 2018;11(1):17–28.
Koike S, Gladders P, Paulus A. Vegetable diseases: a color handbook. New York: Academic Press; 2006.
Penuelas J, Filella I. Responses to a warming world. Science. 2001;294:793–5.
Hoffmann A, Sgro C. Climate change and evolutionary adaptation. Nature. 2011;470:479–85.
Franks SJ, Weber JJ, Aitken SN. Evolutionary and plastic responses to climate change in terrestrial plant populations. Evol Appl. 2014;7(1):123–39.
Franks S, Kane N, O’Hara N, Tittes S, Rest J. Rapid genome-wide evolution in Brassica rapa populations following drought revealed by sequencing of ancestral and descendant gene pools. Mol Ecol. 2016;25(15):3622–31.
van Wees S, Chang H, Zhu T, Glazebrook J. Characterization of the early response of Arabidopsis to Alternaria brassicicola infection using expression profiling. Plant Physiol. 2003;132:606–17.
Mukherjee A, Lev S, Gepstein S, Horwitz B. A compatible interaction of Alternaria brassicicola with Arabidopsisthaliana ecotype DiG: evidence for a specific transcriptional signature. BMC Plant Biol. 2009;9:31.
Rotem J. The genus alternaria: biology, epidemiology, and pathogenicity. St. Paul: APS Press; 1994.
Meena PD, Awasthi RP, Chattopadhyay C, Kolte SJ, Kumar A. Alternaria blight: a chronic disease in rapeseed-mustard. J Oilseed Brassica. 2010;1(1):1–11.
O’Hara N, Rest J, Franks S. Factors affecting the disease severity of Alternaria blackspot in natural Brassica rapa populations on the California and Oregon coasts. Madrono. 2016;63:249–57.
Mitchell-Olds T, Schmitt J. Genetic mechanisms and evolutionary significance of natural variation in Arabidopsis. Nature. 2006;441:947–52.
Metcalf C, Mitchell-Olds T. Life history in a model system: opening the black box with Arabidopsisthaliana. Ecol Lett. 2009;12:593–600.
Walsh B. Quantitative genetics, version 3.0: where have we gone since 1987 and where are we headed? Genetica. 2009;136:213–23.
Rosas U, Mei Y, Xie Q, Banta J, Zhou R, Seufferheld G, et al. Variation in Arabidopsis flowering time associated with cis-regulatory variation in CONSTANS. Nat Commun. 2014;5:3651.
Nowicki M, Nowakowska M, Niezgoda A, Kozik E. Alternaria black spot of crucifers: symptoms, importance of disease, and perspectives of resistance breeding. Veg Crops Res Bull. 2012;76:5–19.
Xie D, Feys B, James S, Nieto-Rostro M, Turner J. COI1: an Arabidopsis gene required for jasmonate-regulated defense and fertility. Science. 1998;280:1091–4.
Zhou C, Zhang L, Duan J, Miki B, Wu K. HISTONE DEACETYLASE19 is involved in jasmonic acid and ethylene signaling of pathogen response in Arabidopsis. Plant Cell. 2005;17:1196–204.
Wu K, Zhang L, Zhou C, Yu C, Chaikam V. HDA6 is required for jasmonate response, senescence and flowering in Arabidopsis. J Exp Bot. 2008;59(2):225–34.
Garg R, Tyagi A, Jain M. Microarray analysis reveals overlapping and specific transcriptional responses to different plant hormones in rice. Plant Signal Behav. 2012;7(8):951.
McKay J, Rchards J, Mitchell-Olds T. Genetics of drought adaptation in Arabidopsis thaliana: I. Pleiotropy contributes to genetic correlations among ecological traits. Mol Ecol. 2003;12:1137–51.
Williams G. Pleiotropy, natural selection, and the evolution of senescence. Evolution. 1957;11(4):398–411.
Simms EL, Rausher MD. Costs and benefits of plant resistance to herbivory. Am Nat. 1987;130(4):570–81.
Mole S. Trade-offs and constraints in plant-herbivore defense theory: a life-history perspective. Oikos. 1994;71:3–12.
Ascencio-Ibanez J, Sozzani R, Lee T, Chu T, Wolfinger R, Cella R, Hanley-Bowdoin L. Global analysis of Arabidopsis gene expression uncovers a complex array of changes impacting pathogen response and cell cycle during geminivirus infection. Plant Physiol. 2008;148:436–54.
Lamesch P, Berardini T, Li D, Swarbreck D, Wilks C, Sasidharan R, Muller R, Dreher K, Alexander D, Garcia-Hernandez M, et al. The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools. Nucleic Acids Res. 2011;40:D1202–10.
Hoth S, Morgante M, Sanchez J, Hanafey M, Tingey S, Chua N. Genome-wide gene expression profiling in Arabidopsis thaliana reveals new targets of abscisic acid and largely impaired gene regulation in the abi1-1 mutant. J Cell Sci. 2002;115(24):4891–900.
Nemhauser J, Hong F, Chory J. Different plant hormones regulate similar processes through largely nonoverlapping transcriptional responses. Cell. 2006;126:467–75.
Journot-Catalino N, Somssich I, Roby D, Kroj T. The transcription factors WRKY11 and WRKY17 act as negative regulators of basal resistance in Arabidopsis thaliana. Plant Cell. 2006;18:3289–302.
Huang D, Wu W, Abrams S, Cutler A. The relationship of drought-related gene expression in Arabidopsis thaliana to hormonal and environmental factors. J Exp Bot. 2008;59(11):2991–3007.
Liu H, Yang W, Dongcheng L, Han Y, Zhang A, Shaohua L. Ectopic expression of a grapevine transcription factor VvWRKY11 contributes to osmotic stress tolerance in Arabidopsis. Mol Biol Rep. 2011;38:417–27.
McGrath K, Dombrecht B, Manners J, Schenk PM, Edgar C, Maclean D, Scheible W, Udvardi M, Kazan K. Repressor- and activator-type ethylene response factors functioning in jasmonate signaling and disease resistance identified via a genome-wide screen of Arabidopsis transcription factor gene expression. Plant Physiol. 2005;139:949–59.
Vance C, Kirk T, Sherwood R. Lignification as a mechanism of disease resistance. Annu Rev Phytopathol. 1980;18:259–88.
Taiz L, Zeiger E. Plant physiology. 4th ed. Sunderland: Sinauer; 2006.
Mohr P, Cahill D. Suppression by ABA of salicylic acid and lignin accumulation and the expression of multiple genes, in Arabidopsis infected with Pseudomonas syringae pv. tomato. Funct Integr Genom. 2007;7:181–91.
Wang X. The genome of the mesopolyploid crop species Brassica rapa. Nat Genet. 2011;43:1035–9.
Hartl DL, Clark A. Principles of population genetics. Oxford: Oxford University Press; 2007.
Karlsson E, Baranowska I, Wade C, Hillbertz N, Zody M, Anderson N, et al. Efficient mapping of mendelian traits in dogs through genome-wide association. Nat Genet. 2007;39:1321–8.
Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.
Buerkle A, Gompert Z. Population genomics based on low coverage sequencing: how low should we go? Mol Ecol. 2013;22:3028–35.
Storey J, Tibshirani R. Statistical significance for genomewide studies. Proc Natl Acad Sci. 2003;100:9440.
Lawrence C, Mitchell T, Craven K, Cho Y, Cramer R, Kim K. At death’s door: Alternaria pathogenecity mechanisms. Plant Pathol. 2008;24(2):101–11.
Marra J, Eo S, Hale M, Waser P, DeWoody J. A priori and a posteriori approaches for finding genes of evolutionary interest in non-model species: osmoregulatory genes in the kidney transcriptome of the desert rodent Dipodomys spectabilis (banner-tailed kangaroo rat). Comp Biochem Physiol D Genom Proteom. 2012;7(4):328–39.
Guo Y, Samans B, Chen S, Kibret K, Hatzig S, Turner N, et al. Drought-tolerant brassica rapa shows rapid expression of gene networks for general stress responses and programmed cell death under simulated drought stress. Plant Mol Biol Rep. 2017;35:416–30.
Gillis J, Mistry M, Pavlidis P. Gene function analysis in complex data sets using ErmineJ. Nat Protoc. 2010;5(6):1148–59.
We thank A.E. Weis and E. Austen for providing seeds, R. Welt for growing plants and DNA extraction, and G. Amidi-Abraham for laboratory assistance.
Funding was provided by a fellowship from AAUW and the Elisha M. Friedman Postdoctoral Fellowship to N. O’Hara, a Grant (DEB-1142784) from the National Science Foundation to S. Franks, and a Grant (R01GM108904) from NIGMS-NIH to J.S.R. Aside from providing funds, funders played no role in this study.
Ethics approval and consent to participate
All fungal work referenced in O'Hara et al. 2016 was conducted in sterile conditions, and was permitted under APHIS license #P526P-11-00130.
Consent for publication
The authors have no conflict of interest to declare.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table of top 50 most differentiated necrotrophic fungal pathogen response genes including FST value, gene description, and notes on pathogen response pathway from previous studies.
About this article
Cite this article
O’Hara, N.B., Franks, S.J., Kane, N.C. et al. Evolution of pathogen response genes associated with increased disease susceptibility during adaptation to an extreme drought in a Brassica rapa plant population. BMC Ecol Evo 21, 61 (2021). https://doi.org/10.1186/s12862-021-01789-7