Transcriptome dynamics of Arabidopsis during sequential biotic and abiotic stresses.

In nature, plants have to cope with a wide range of stress conditions that often occur simultaneously or in sequence. To investigate how plants cope with multi-stress conditions, we analyzed the dynamics of whole-transcriptome profiles of Arabidopsis thaliana exposed to six sequential double stresses inflicted by combinations of: (i) infection by the necrotrophic fungus Botrytis cinerea, (ii) herbivory by chewing larvae of Pieris rapae, and (iii) drought stress. Each of these stresses induced specific expression profiles over time, in which one-third of all differentially expressed genes was shared by at least two single stresses. Of these, 394 genes were differentially expressed during all three stress conditions, albeit often in opposite directions. When two stresses were applied in sequence, plants displayed transcriptome profiles that were very similar to the second stress, irrespective of the nature of the first stress. Nevertheless, significant first-stress signatures could be identified in the sequential stress profiles. Bioinformatic analysis of the dynamics of co-expressed gene clusters highlighted specific clusters and biological processes of which the timing of activation or repression was altered by a prior stress. The first-stress signatures in second stress transcriptional profiles were remarkably often related to responses to phytohormones, strengthening the notion that hormones are global modulators of interactions between different types of stress. Because prior stresses can affect the level of tolerance against a subsequent stress (e.g. prior herbivory strongly affected resistance to B. cinerea), the first-stress signatures can provide important leads for the identification of molecular players that are decisive in the interactions between stress response pathways.

Genome-wide association analysis reveals distinct genetic architectures for single and combined stress responses in Arabidopsis thaliana.

Plants are commonly exposed to abiotic and biotic stresses. We used 350 Arabidopsis thaliana accessions grown under controlled conditions. We employed genome-wide association analysis to investigate the genetic architecture and underlying loci involved in genetic variation in resistance to: two specialist insect herbivores, Pieris rapae and Plutella xylostella; and combinations of stresses, i.e. drought followed by P. rapae and infection by the fungal pathogen Botrytis cinerea followed by infestation by P. rapae. We found that genetic variation in resistance to combined stresses by drought plus P. rapae was limited compared with B. cinerea plus P. rapae or P. rapae alone. Resistance to the two caterpillars is controlled by different genetic components. There is limited overlap in the quantitative trait loci (QTLs) underlying resistance to combined stresses by drought plus P. rapae or B. cinerea plus P. rapae and P. rapae alone. Finally, several candidate genes involved in the biosynthesis of aliphatic glucosinolates and proteinase inhibitors were identified to be involved in resistance to P. rapae and P. xylostella, respectively. This study underlines the importance of investigating plant responses to combinations of stresses. The value of this approach for breeding plants for resistance to combinatorial stresses is discussed.

Genetic architecture of plant stress resistance: multi-trait genome-wide association mapping.

Plants are exposed to combinations of various biotic and abiotic stresses, but stress responses are usually investigated for single stresses only. Here, we investigated the genetic architecture underlying plant responses to 11 single stresses and several of their combinations by phenotyping 350 Arabidopsis thaliana accessions. A set of 214 000 single nucleotide polymorphisms (SNPs) was screened for marker-trait associations in genome-wide association (GWA) analyses using tailored multi-trait mixed models. Stress responses that share phytohormonal signaling pathways also share genetic architecture underlying these responses. After removing the effects of general robustness, for the 30 most significant SNPs, average quantitative trait locus (QTL) effect sizes were larger for dual stresses than for single stresses. Plants appear to deploy broad-spectrum defensive mechanisms influencing multiple traits in response to combined stresses. Association analyses identified QTLs with contrasting and with similar responses to biotic vs abiotic stresses, and below-ground vs above-ground stresses. Our approach allowed for an unprecedented comprehensive genetic analysis of how plants deal with a wide spectrum of stress conditions.

Natural variation in life history strategy of Arabidopsis thaliana determines stress responses to drought and insects of different feeding guilds.

Plants are sessile organisms and, consequently, are exposed to a plethora of stresses in their local habitat. As a result, different populations of a species are subject to different selection pressures leading to adaptation to local conditions and intraspecific divergence. The annual brassicaceous plant Arabidopsis thaliana is an attractive model for ecologists and evolutionary biologists due to the availability of a large collection of resequenced natural accessions. Accessions of A. thaliana display one of two different life cycle strategies: summer and winter annuals. We exposed a collection of 308 European Arabidopsis accessions, that have been genotyped for 250K SNPs, to a range of stresses: one abiotic stress (drought), four biotic stresses (Pieris rapae caterpillars, Plutella xylostella caterpillars, Frankliniella occidentalis thrips and Myzus persicae aphids) and two combined stresses (drought plus P. rapae and Botrytis cinerea fungus plus P. rapae). We identified heritable genetic variation for responses to the different stresses, estimated by narrow-sense heritability. We found that accessions displaying different life cycle strategies differ in their response to stresses. Winter annuals are more resistant to drought, aphids and thrips and summer annuals are more resistant to P. rapae and P. xylostella caterpillars. Summer annuals are also more resistant to the combined stresses of drought plus P. rapae and infection by the fungus Botryris cinerea plus herbivory by P. rapae. Adaptation to drought displayed a longitudinal gradient. Finally, trade-offs were recorded between the response to drought and responses to herbivory by caterpillars of the specialist herbivore P. rapae.

Effect of prior drought and pathogen stress on Arabidopsis transcriptome changes to caterpillar herbivory.

In nature, plants are exposed to biotic and abiotic stresses that often occur simultaneously. Therefore, plant responses to combinations of stresses are most representative of how plants respond to stresses. We used RNAseq to assess temporal changes in the transcriptome of Arabidopsis thaliana to herbivory by Pieris rapae caterpillars, either alone or in combination with prior exposure to drought or infection with the necrotrophic fungus Botrytis cinerea. Pre-exposure to drought stress or Botrytis infection resulted in a significantly different timing of the caterpillar-induced transcriptional changes. Additionally, the combination of drought and P. rapae induced an extensive downregulation of A. thaliana genes involved in defence against pathogens. Despite a more substantial growth reduction observed for plants exposed to drought plus P. rapae feeding compared with P. rapae feeding alone, this did not affect weight increase of this specialist caterpillar. Plants respond to combined stresses with phenotypic and transcriptional changes that differ from the single stress situation. The effect of a previous exposure to drought or B. cinerea infection on transcriptional changes to caterpillars is largely overridden by the stress imposed by caterpillars, indicating that plants shift their response to the most recent stress applied.

Genome-wide association analysis reveals genes controlling an antagonistic effect of biotic and osmotic stress on Arabidopsis thaliana growth.

While the response of Arabidopsis thaliana to drought, herbivory or fungal infection has been well-examined, the consequences of exposure to a series of such (a)biotic stresses are not well studied. This work reports on the genetic mechanisms underlying the Arabidopsis response to single osmotic stress, and to combinatorial stress, either fungal infection using Botrytis cinerea or herbivory using Pieris rapae caterpillars followed by an osmotic stress treatment. Several small-effect genetic loci associated with rosette dry weight (DW), rosette water content (WC), and the projected rosette leaf area in response to combinatorial stress were mapped using univariate and multi-environment genome-wide association approaches. A single-nucleotide polymorphism (SNP) associated with DROUGHT-INDUCED 19 (DI19) was identified by both approaches, supporting its potential involvement in the response to combinatorial stress. Several SNPs were found to be in linkage disequilibrium with known stress-responsive genes such as PEROXIDASE 34 (PRX34), BASIC LEUCINE ZIPPER 25 (bZIP25), RESISTANCE METHYLATED GENE 1 (RMG1) and WHITE RUST RESISTANCE 4 (WRR4). An antagonistic effect between biotic and osmotic stress was found for prx34 and arf4 mutants, which suggests PRX34 and ARF4 play an important role in the response to the combinatorial stress.

High-quality chromosome-scale genomes facilitate effective identification of large structural variations in hot and sweet peppers.

Pepper (Capsicum annuum) is an important vegetable crop that has been subjected to intensive breeding, resulting in limited genetic diversity, especially for sweet peppers. Previous studies have reported pepper draft genome assemblies using short read sequencing, but their capture of the extent of large structural variants (SVs), such as presence-absence variants (PAVs), inversions, and copy-number variants (CNVs) in the complex pepper genome falls short. In this study, we sequenced the genomes of representative sweet and hot pepper accessions by long-read and/or linked-read methods and advanced scaffolding technologies. First, we developed a high-quality reference genome for the sweet pepper cultivar 'Dempsey' and then used the reference genome to identify SVs in 11 other pepper accessions and constructed a graph-based pan-genome for pepper. We annotated an average of 42 972 gene families in each pepper accession, defining a set of 19 662 core and 23 115 non-core gene families. The new pepper pan-genome includes informative variants, 222 159 PAVs, 12 322 CNVs, and 16 032 inversions. Pan-genome analysis revealed PAVs associated with important agricultural traits, including potyvirus resistance, fruit color, pungency, and pepper fruit orientation. Comparatively, a large number of genes are affected by PAVs, which is positively correlated with the high frequency of transposable elements (TEs), indicating TEs play a key role in shaping the genomic landscape of peppers. The datasets presented herein provide a powerful new genomic resource for genetic analysis and genome-assisted breeding for pepper improvement.