Evolution of pathogen response genes associated with increased disease susceptibility during adaptation to an extreme drought in a Brassica rapa plant population

doi:10.21203/rs.2.16267/v3

Download PDF

Research article

Evolution of pathogen response genes associated with increased disease susceptibility during adaptation to an extreme drought in a Brassica rapa plant population

https://doi.org/10.21203/rs.2.16267/v3

This work is licensed under a CC BY 4.0 License

Journal Publication

published 21 Apr, 2021

Read the published version in BMC Ecology and Evolution →

You are reading this latest preprint version

BACKGROUND: 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.

RESULTS: 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, 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 F_ST 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.

CONCLUSIONS: 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.

Evolutionary Biology

Evolutionary Developmental Biology

Climate change

Brassica rapa

Alternaria brassicae

drought

resurrection approach

necrotrophic fungal pathogen response

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-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 [8], stimulate a response pathway generally mediated by jasmonic acid (JA) and ethylene (ET) [9]. 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 [10]. ABA is thus an example of a regulatory system that may potentially be involved in adaptation to drought [11], 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 [12], 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) [15]. 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-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 [17]. 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 [22]. We previously found that this population experienced a genetically-based evolutionary increase in Alternaria blackspot disease susceptibility [4], together with an adaptive shift to escape drought through earlier flowering [17]. 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-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 [26], 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 F_ST 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 (~25x 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 [26]. 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 [4]. 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 1,182 candidate NFPR genes were significantly differentiated between ancestors and descendants (Figure 1), with an F_STabove background compared to all the genes, as determined by outlier F_ST 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 [26]. 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 F_ST for all NFPR genes was not significantly elevated above background (ANOVA: F _1,35199 = 0.008, P = 0.929). In a Manhattan plot comparing F_ST of NFPR to all other genes across the genome, the NFPR and all genes’ trend lines tracked one another (Figure 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 F_ST 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 F_ST values) in a Venn diagram (Figure 3). 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).

HISTORICAL SELECTION

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 (F_ST) NFPR genes and the rest of the NFPR genes (ANOVA: ancestors F _1,1180 = 2.398, P = 0.122; descendants F _1,1180 = 0.551, P = 0.458). This was consistent with our genome-wide result, where we found little overlap between genes with high F_ST and negative Tajima’s D [26]. 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 (F_ST) NFPR genes were in regions of reduced Tajima’s D (red bars in Figure 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 F_ST 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 (Figure 2). 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 1,182 of COI1/JA genes which are differentially expressed in response to pathogens in van Wees et al. (2003) [27] and Mukherjee et al. (2009) [28] 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 (F_ST). We found that pathogen response categories were not significantly enriched and were not highly ranked out of the 3,136 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 [17] and increased fungal disease susceptibility [4] 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 [31]. 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-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 [26].

EVOLUTION OF NFPR GENES

Using an outlier F_ST 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 [37]. 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 F_ST 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 F_ST 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 [26].

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 [4]. 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-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 [41]. Furthermore, it is possible that antagonistic pleiotropy was playing a role in the evolution observed [42]. 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 [4]. Indeed, the increased susceptibility we observed in this population was accompanied by an evolutionary shift to thinner leaves (increased specific leave area) [4], which could indicate a trade-off between drought escape and disease susceptibility in agreement with the growth and defense trade-off [43, 44].

HISTORICAL SELECTION

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 F_ST and Tajima’s D, and there were no genes that were outliers for both high F_ST (>0.2) and low Tajima’s D (<-2). This finding concurs with other prior work in this system [26]. 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) (Figure 2A; F_ST = 0.21; q = 1.75 x 10^-5) which has been found to be differentially expressed in Arabidopsis in response to infection with A. brassicicola [27] and cabbage leaf curl virus [45]. 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 [46]. 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)(F_ST = 0.20; q = 7.30 x 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 F_ST, including Bra012746 and Bra033968, and Bra011282, have pleiotropic roles in both abiotic and biotic stress response [49-51]. Studies have shown that ERF4 is induced by ABA, and ERF4 suppresses expression of PDF1.2 [52], a defense effector that is elicited by JA signaling [6]. This gene showed a shift in allele frequency at three non-synonymous sites, and evolution at these sites could have a functional effect (Figure 2B). 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 [27] and also had a significant F_ST 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 [53], which works primarily through leaf structure [54]. 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 [55]. 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 [4], 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 [17]. 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 [26].

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.

PLANT SAMPLES

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 [17].

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. (2016). 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. (2016) 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 [4] due to the greater evolutionary shifts in multiple traits experienced by that population [17].

BIOINFORMATIC ANALYSIS AND VALIDATION

We previously used computational methods described in Franks et al. (2016) 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/) [56], calling sequence variants, and calculating population genetic estimates including divergence (F_ST) [57, 58] and Tajima's D [59] 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 8X 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.5x-1x) population-level (>20 individuals) approach is used [60]. 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 [26].

In our previous work, statistical significance of genes, results which are used here, were determined for the F_ST values of genes by comparing the experimental F_ST 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 [61]. We also report the top 50 most differentiated NFPR genes, including those with a q value of >0.05.

PLANT-PATHOGEN SYSTEM

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 [4]. The molecular mechanisms of pathogen response in the Brassica rapa – Alternaria brassicae system is not well understood [62], 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 [63] 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 (~1,100) 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 2005 for further details). Additionally, a recent B. rapa expression study was used to identify candidate drought response genes [64].

FUNCTIONAL ANALYSIS

In the current study, for our a priori analysis approach, a database of 1,182 well-characterized candidate NFPR genes (Additional file 2) was assembled from two previous A. thaliana-A. brassicicola studies by Mukherjee et al. (2009) and van Wees et al. (2003); 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 [36], including genes shown to be involved in both compatible interactions [28] and incompatible interactions [27]. Briefly, to identify candidate NFPR genes, van Wees et al. (2003) used mutants with non-functional copies of genes involved in JA and SA signaling followed by inoculation and expression profiling, and Mukherjee et al. (2009) 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 7,076 genes candidate drought response genes (Additional file 3) which were differentially expressed in Guo et al. (2017) study in response to a polyethylene glycol (PEG) treatment, as compared to controls (no PEG), and that had a >2 fold change at 12 hours 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 [65], with annotations from Arabidopsis thaliana. To analyze our F_ST 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.

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

Not applicable.

AVAILABILITY OF DATA AND MATERIAL

The dataset supporting the conclusions of this article are included within the article (and its additional files).

COMPETING INTERESTS

The authors have no conflict of interest to declare.

FUNDING

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.

AUTHOR CONTRIBUTIONS

NO, SF, NK, ST, and JR designed the experiments; NO managed the project, and lead the analysis. NO lead the writing of the manuscript and all authors contributed to writing and editing; NO, SF, NK, ST, and JR contributed to data analyses; and all authors have read, edited, and approved of the final manuscript.

ACKNOWLEDGEMENTS

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.

COI1: coronatine-insensitive 1

JA: Jasmonic Acid

ET: ethylene

ABA: Abscisic Acid

NFPR: necrotrophic fungal pathogen response

SNP: Single nucleotide polymorphism

SSH: suppression subtractive hybridization

PEG treatment: polyethylene glycol treatment

GO Analysis: gene ontology analysis

FDR: False discover rate

Desprez-Loustau ML, Marcais B, Nageleisen L, Piou D, Vannini A: Interactive effects of drought and pathogens in forest trees. Annals of Forest Science 2006, 63:597-612.
Garrett K, Thomas-Sharma S, Forbes G, Nopsa J: Climate change and plant pathogen invasions. In: Invasive Species and Global Climate Change. Edited by Ziska L, Dukes J: CABI; 2014.
Pautasso M, Doring T, Garbelotto M, Lorenzo P, Jeger M: Impacts of climate change on plant disease - opinions and trends. European Journal of Plant Pathology 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-248.
Garrett KA, Dendy SP, Frank EE, Rouse MN, Travers SE: Climate change effects on plant disease: genomes to ecosystems. Annual Review of Phytopathology 2006, 44:489-509.
Glazebrook J: Contrasting mechanisms of defense against biotrophic and necrotrophic pathogens. Annual Review of Phytopathology 2005, 43:205 - 227.
Mengiste T: Plant immunity to necrotrophs. Annual Review of Phytopathology 2012, 50:267-294.
Agrios GN: Plant Pathology, 5th edn. Burlington, MA: Academic Press; 2005.
McDowell J, Dangl J: Signal transduction in the plant immune response. Trends in Biochemical Sciences 2000, 25:79 - 82.
Mauch-Mani B, Mauch F: The role of abscisic acid in plant-pathogen interactions. Current Opinion in Plant Biology 2005, 8:409 - 414.
Shinozaki K, Yamaguchi-Shinozaki K, Seki M: Regulatory network of gene expression in the drought and cold stress responses. Current Opinion in Plant Biology 2003, 6:410-417.
Thrall PH, Burdon JJ: Evolution of virulence in a plant host-pathogen metapopulation. Science 2003, 299:1735 - 1737.
Stam R, Scheikl D, Tellier A: Pooled Enrichment Sequencing Identifies Diversity and Evolutionary Pressures at NLR Resistance Genes within a Wild Tomato Population. Genome Biology and Evolution 2016, 8(5):1501-1515.
Stam R, Silva-Arias G, Tellier A: Subsets of NLR genes show differential signatures of adaptation during colonization of new habitats. New Phytologist 2019, 224:367-379.
Jorgensen T, Emerson B: Functional variation in disease resistance gene in populations of Arabidopsis thaliana. Molecular Ecology 2008, 17(22):4912 - 4923.
Laine A, Burdon J, Dodds P, Thrall P: Spatial variation in disease resistance: from molecules to metapopulations. Journal of Ecology 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. Proceedings of the National Academy of Sciences 2007, 104(4):1278-1282.
Franks SJ, Weis AE: A change in climate causes rapid evolution of multiple life-history traits and their interactions in an annual plant. Journal of Evolutionary Biology 2008, 21:1321-1334.
Franks SJ: Plasticity and evolution in drought avoidance and escape in the annual plant Brassica rapa. New Phytologist 2011, 190(1):249-257.
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 - 873.
Franks S, Hamann E, Weis AE: Using the resurrection approach to understand contemporary evolution in changing environments. Evolutionary Applications 2018, 11(1):17-28.
Koike S, Gladders P, Paulus A: Vegetable Diseases: a Color Handbook. USA: Academic Press; 2006.
Penuelas J, Filella I: Responses to a warming world. Science 2001, 294:793 - 795.
Hoffmann A, Sgro C: Climate change and evolutionary adaptation. Nature 2011, 470:479 - 485.
Franks SJ, Weber JJ, Aitken SN: Evolutionary and plastic responses to climate change in terrestrial plant populations. Evolutionary Applications 2014, 7(1):123 -139.
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. Molecular Ecology 2016, 25(15):3622-3631.
van Wees S, Chang H, Zhu T, Glazebrook J: Characterization of the early response of Arabidopsis to Alternaria brassicicola infection using expression profiling. Plant Physiology 2003, 132:606 - 617.
Mukherjee A, Lev S, Gepstein S, Horwitz B: A compatible interaction of Alternaria brassicicola with Arabidopsis thaliana ecotype DiG: evidence for a specific transcriptional signature. BMC Plant Biology 2009, 9.
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. Journal of 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-257.
Mitchell-Olds T, Schmitt J: Genetic mechanisms and evolutionary significance of natural variation in Arabidopsis. Nature 2006, 441:947-952.
Metcalf C, Mitchell-Olds T: Life history in a model system: opening the black box with Arabidopsis thaliana. Ecology Letters 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-223.
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. Nature Communications 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. Vegetable Crops Research Bulletin 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 -1094.
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-1204.
Wu K, Zhang L, Zhou C, Yu C, Chaikam V: HDA6 is required for jasmonate response, senescence and flowering in Arabidopsis. Journal of Experimental Botany 2008, 59(2):225-234.
Garg R, Tyagi A, Jain M: Microarray analysis reveals overlapping and specific transcriptional responses to different plant hormones in rice Plant Signaling & Behavior 2012, 7(8).
McKay J, Rchards J, Mitchell-Olds T: Genetics of drought adaptation in Arabidopsis thaliana: I. Pleiotropy contributes to genetic correlations among ecological traits. Molecular Ecology 2003, 12:1137-1151.
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. The American Naturalist 1987, 130(4):570 - 581.
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 Physiology 2008, 148:436 - 454.
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 Research 2011.
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. Journal of Cell Science 2002, 115(24):4891 - 4900.
Nemhauser J, Hong F, Chory J: Different plant hormones regulate similar processes through largely nonoverlapping transcriptional responses. Cell 2006, 126:467 - 475.
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. The Plant Cell 2006, 18:3289 - 3302.
Huang D, Wu W, Abrams S, Cutler A: The relationship of drought-related gene expression in Arabidopsis thaliana to hormonal and environmental factors. Journal of Experimental Botany 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. Molecular Biology Reports 2011, 38:417 - 427.
McGrath K, Dombrecht B, Manners J, PM S, 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 Physiology 2005, 139:949 - 959.
Vance C, Kirk T, Sherwood R: Lignification as a mechanism of disease resistance. Annual Review of Phytopathology 1980, 18:259 - 288.
Taiz L, Zeiger E: Plant Physiology, 4th edn. Sunderland, Mass.: 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. Functional Integrative Genomics 2007, 7:181 - 191.
Wang X: The genome of the mesopolyploid crop species Brassica rapa. Nature Genetics 2011.
Hartl DL, Clark A: Principles of Population Genetics: 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. Nature Genetics 2007, 39:1321-1328.
Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 1989, 123:585-595.
Buerkle A, Gompert Z: Population genomics based on low coverage sequencing: how low should we go? Molecular Ecology 2013, 22:3028 - 3035.
Storey J, Tibshirani R: Statistical significance for genomewide studies. Proceedings of the National Academy of Sciences 2003, 100:9440.
Lawrence C, Mitchell T, Craven K, Cho Y, Cramer R, Kim K: At death's door: Alternaria pathogenecity mechanisms. Plant Pathology 2008, 24(2):101-111.
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). Comparative Biochemistry and Physiology Part D: Genomics and Proteomics 2012, 7(4):328 - 339.
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 Molecular Biology Reporter 2017, 35:416-430.
Gillis J, Mistry M, Pavlidis P: Gene function analysis in complex data sets using ErmineJ Nature Protocols 2010, 5(6):1148-1159.

Table 1. The 20 most highly differentiated (high F_ST) necrotrophic fungal pathogen response genes. Significantly differentiated genes (based on q value) are bolded. Asterisks show significantly differentiated genes that are involved in both drought and pathogen stress response (van Wees et al., 2003; Mukherjee et al., 2009), as determined by expression analysis (1), or where genes have been shown to be COI1 dependent (2). Descriptions were retrieved from http://brassicadb.org/brad/.

*A. thaliana* ID	*B. rapa* ID	F_ST	q value	Gene Name	Notes	Description
AT2G21520	Bra030295	0.21	1.75 x 10^-5		1	Sec14p-like phosphatidylinositol transfer family
AT3G15210	Bra038107*	0.20	7.30 x 10^-5	ERF4, RAP2.5	1,2	ethylene responsive element binding factor 4
AT5G54160	Bra026320*	0.20	1.45 x 10-⁴	OMT1	1	O-methyltransferase 1
AT2G02390	Bra017417	0.19	4.16 x 10^-4	GST18,GSTZ1	1	glutathione S-transferase zeta 1
AT4G01960	Bra036302	0.17	4.48 x 10^-3		1
AT4G39270	Bra033611*	0.17	5.63 x 10^-3		1	Leucine-rich repeat protein kinase family
AT4G31550	Bra011282*	0.16	0.02	WRKY11	1	WRKY DNA-binding protein 11
AT1G23040	Bra012356*	0.16	0.02		1	hydroxyproline-rich glycoprotein family protein
AT2G33120	Bra021828	0.16	0.02	SAR1,VAMP722	1	synaptobrevin-related protein 1
AT1G67980	Bra033968	0.16	0.03	CCOAMT	1,2	caffeoyl-CoA 3-O-methyltransferase
AT2G17720	Bra024495	0.15	0.05		1	2-oxoglutarate (2OG) and Fe(II)-dependent oxygenase superfamily protein
AT1G45145	Bra036335	0.15	0.06	ATH5,LIV1,TRX5	1,2	thioredoxin H-type 5
AT1G74020	Bra003820	0.15	0.07	SS2	1,2	strictosidine synthase 2
AT1G23730	Bra012384	0.15	0.07	BCA3	1	beta carbonic anhydrase 3
AT4G12720	Bra012746	0.14	0.11	GFG1,NUDT7	1	MutT/nudix family protein
AT2G02930	Bra018543	0.14	0.14	GST16,GSTF3	1	glutathione S-transferase F3
AT1G10140	Bra031702	0.14	0.16		1,2	Uncharacterised conserved protein UCP031279
AT3G02520	Bra040592	0.14	0.17	GF14 NU,GRF7	1	general regulatory factor 7
AT5G65750	Bra024417	0.14	0.19		1	2-oxoglutarate dehydrogenase, E1 component
AT1G16520	Bra026043	0.14	0.22		1

Download PDF

Journal Publication

published 21 Apr, 2021

Read the published version in BMC Ecology and Evolution →

Editorial decision: Accept
06 Apr, 2021
Reviewers invited by journal
22 Nov, 2020
Editor assigned by journal
29 Oct, 2020
Editor invited by journal
29 Oct, 2020
Submission checks completed at journal
28 Oct, 2020

You are reading this latest preprint version

Evolution of pathogen response genes associated with increased disease susceptibility during adaptation to an extreme drought in a Brassica rapa plant population

Status:

Journal Publication

Version 3

Abstract

Figures

Background

Results

Discussion

Conclusions

Methods

Declarations

Abbreviations

References

Tables

Supplementary Files

Status:

Journal Publication

Version 3