Transcriptional response of a target plant to benzoxazinoid and diterpene allelochemicals highlights commonalities in detoxification
BMC Plant Biology volume 22, Article number: 402 (2022)
Plants growing in proximity to other plants are exposed to a variety of metabolites that these neighbors release into the environment. Some species produce allelochemicals to inhibit growth of neighboring plants, which in turn have evolved ways to detoxify these compounds.
In order to understand how the allelochemical-receiving target plants respond to chemically diverse compounds, we performed whole-genome transcriptome analysis of Arabidopsis thaliana exposed to either the benzoxazinoid derivative 2-amino- 3H-phenoxazin-3-one (APO) or momilactone B. These two allelochemicals belong to two very different compound classes, benzoxazinoids and diterpenes, respectively, produced by different Poaceae crop species.
Despite their distinct chemical nature, we observed similar molecular responses of A. thaliana to these allelochemicals. In particular, many of the same or closely related genes belonging to the three-phase detoxification pathway were upregulated in both treatments. Further, we observed an overlap between genes upregulated by allelochemicals and those involved in herbicide detoxification. Our findings highlight the overlap in the transcriptional response of a target plant to natural and synthetic phytotoxic compounds and illustrate how herbicide resistance could arise via pathways involved in plant-plant interaction.
In the competition for nutrients and space, some plants secrete growth-inhibiting chemicals into the soil to inhibit the germination or growth of neighboring plants [1, 2]. These phytotoxic compounds, known as allelochemicals, can be released from the roots of the donor plant or leach into the soil from leaves or decaying plant material. Several species of the Poaceae, the economically most important order of plants, produce such allelochemicals, but their chemical nature can be quite diverse. For example, maize, wheat, and rye all produce benzoxazinoids (BX) [3,4,5], while many rice cultivars produce diterpenes such as momilactone B [6,7,8]. There is currently no example of a species that is capable of producing both momilactone B and BX compounds. This might suggest that—despite their very different chemical characteristics—these two types of allelochemicals exert redundant activities on target plants.
The two main forms of BX in Poaceae are DIBOA (2,4-dihydroxy-2H-1,4-benzoxazin-3(4H)-one) and its C-7-methoxy derivative DIMBOA, for which the biosynthetic pathways have been characterized in maize and partly in wheat (reviewed in ). BXs enter the soil either via exudation from the plant roots or via decomposing plant material that is incorporated into the soil. In an agricultural context, this can happen for example when grasses are used as cover crop or when straw is left on the field post-harvest as green manure . In the soil, microbes rapidly degrade DIBOA to the final, stable product 2-amino-3H-phenoxazin-3-one (APO) . APO and the analogous, DIMBOA-derived AMPO are potent phytotoxins  that inhibit histone deacetylases and slow down root growth .
As alluded to above, rice produces a very different class of allelochemicals than wheat, maize, and rye. Upon environmental cues, ranging from neighbor proximity to fungal attack to abiotic stress, rice plants produce the diterpenes momilactone A and B (reviewed in ). Momilactone A contributes to resistance to fungal pathogens [14, 15], while momilactone B has stronger allelochemical properties, inhibiting germination and root growth of a broad range of target species . The biosynthetic genes required for the production of momilactone B are known, and the pathway has recently been reconstituted in [17,18,19,20]. Like DIBOA and DIMBOA, momilactone B is exuded from the roots by a yet unknown mechanism. Although its growth-inhibitory properties are well known, its mode of action remains unresolved .
As some plants release allelochemicals to inhibit growth of neighboring plants, their neighbors in turn have evolved ways to detoxify these compounds . As for xenobiotics such as herbicides, allelochemicals taken up by the target plant may be detoxified through a three-phase detoxification system [22, 23], by which the compounds are metabolically activated (phase I) to allow conjugation with sugars or amino acids (phase II). Soluble conjugates may then be transported to and stored in the vacuole or in the apoplast (phase III). Different classes of enzymes and proteins contribute to the different stages of detoxification: cytochrome P450s (CYPs), hydroxylases, and peroxidases activate xenobiotics in phase I, while phase II-detoxifying enzymes include UDP-dependent glucosyltransferases (UGTs), glutathione-S-transferases (GSTs), and quinone oxidoreductases . Membrane transporters such as ATP-binding cassette (ABC), multi-antimicrobial extrusion (MATE) and major facilitator superfamily (MFS) transporters then remove the conjugates from the cytosol in phase III by sequestering them either in the vacuole or exporting them into the apoplast .
Here, we asked whether chemically distinct allelochemicals from different donor species but with similar growth arrest effects and potency towards the model plant Arabidopsis thaliana [12, 16, 19, 26] trigger similar molecular responses in the target plant. We explored the transcriptional response of A. thaliana to two different and agriculturally relevant allelochemicals, the benzoxazinoid APO and the diterpene momilactone B. We profiled the immediate, short-term transcriptional response of A. thaliana seedlings to half-maximal effect concentrations of either compound by time series mRNA sequencing (RNA-seq). We hypothesized that by determining the commonalities and differences in molecular response to these two compound classes, we might be able to discern their respective mode of action and shed light on coping mechanisms of the target plant. The respective transcriptional response to the two allelochemicals showed substantial functional overlap. In particular, we identified components of the three-phase detoxification pathway to be up-regulated in both treatments. This suggests that, despite their distinct chemical characteristics, both types of allelochemicals trigger similar xenobiotic detoxification responses in the target plant, and that the involved enzymes have promiscuous functionality in detoxifying xenobiotics and other phytotoxic compounds.
Different allelochemicals trigger similar transcriptional activation
To determine the transcriptional response of A. thaliana to the two allelochemicals APO and momilactone B, we exposed A. thaliana Columbia-0 (Col-0) seedlings to the respective half-maximal effect concentration (EC50) of either compound (3.5 µM for APO  and 4 µM for momilactone B; Supplemental Fig. 1). We sampled plant material for RNA extraction and RNA-seq library preparation after 1 h, 6 h, and 24 h of treatment (Supplemental Table 1). As control, we treated seedlings with the equivalent concentration of solvent (dimethyl sulfoxide, DMSO). Sequencing reads were mapped to the A. thaliana TAIR10 reference genome (arabidopsis.org) . In a principal component analysis (PCA) of read counts of the 1,000 genes with the highest variance after variance stabilizing transformation using the DESeq2 package , samples clustered by time point; treated groups separated from untreated ones, indicating transcriptional responses to either allelochemical (Fig. 1A). To investigate which genes were differentially expressed in response to the treatments, we performed differential expression analysis relative to the untreated 0 h baseline sample and combined this with clustering of genes using weighted gene correlation network analysis (WGCNA)  (Fig. 1B). For APO-treated samples, we observed three clusters (A3, A6, and A7) that contained genes that were more strongly up-regulated in APO than in control treatments. For the momilactone B-treated samples, the genes contained in cluster M1 were up-regulated upon momilactone B treatment.
We proceeded to investigate if the up-regulated genes were involved in specific metabolic or cellular pathways and performed over-representation analysis (ORA) of the genes contained in the clusters A3, A6, A7, and M1, respectively, using the publicly available GOslim annotation of A. thaliana (arabidopsis.org) (Fig. 2). We found that both treatments up-regulated genes involved in the activation of xenobiotics, such as, e.g., cytochromes, oxidoreductases, or in glutathione-mediated detoxification, glycosylation, and transport, which are pillars of the three-phase xenobiotic detoxification system.
General detoxification pathway genes are up-regulated upon allelochemical treatment
We next asked if the two allelochemicals simply activated the same general detoxification response, or if there was actual overlap among the differentially expressed genes. To this end, we performed an in-depth analysis of the genes involved in the different phases of the detoxification process: cytochrome P450 oxidases (CYPs) (phase I), glutathione-S-transferases (GSTs) and UPD-dependent glycosyltransferases (UGTs) (phase II), and transporters (phase III).
CYPs play a well-established role in the metabolic activation of xenobiotic compounds in phase I of the detoxification . In our data, five CYP genes (CYP71B15, CYP72A8, CYP81D8, CYP81D11, and CYP89A5) were up-regulated in both the APO and the momilactone B treatment, while eleven CYP genes were up-regulated only in APO and another eleven only in momilactone B treatment (Fig. 3 and Supplemental Table 2). The up-regulated genes were distributed among all CYP clans. Noticeably, both treatments up-regulated genes from the CYP81 family; CYP81D8 and CYP81D11 were up-regulated in both treatments, while other members were up-regulated upon either individual treatment (Fig. 3 and Supplemental Table 2).
To further understand the degree to which the observed up-regulation was compound-specific, we added public information on the A. thaliana transcriptome response to the herbicide safener fenclorim  and to 2-benzoxazolinone (BOA), an intermediate in the bioconversion of DIBOA to APO . Several of the genes up-regulated by APO and/or momilactone B were also up-regulated upon fenclorim or BOA treatment (Fig. 3).
Phase II detoxification involves the conjugation of sugar or glutathione to the activated compounds. Our initial analysis had already identified several UGTs and GSTs as being significantly up-regulated upon allelochemical treatment (Fig. 1C). Multiple UGTs were up-regulated by both treatments (Fig. 4A), with many belonging to the UGT73 family. While most members of this family were up-regulated upon momilactone B treatment and by fenclorim , APO and its precursor BOA up-regulated only members of the UGT73B family. In addition to UGT73s, a few other UGTs from different families were up-regulated upon momilactone B treatment, many of these overlapping with UGTs up-regulated by herbicide safener treatment (Fig. 4A) .
Similar to the UGTs, several of the up-regulated GSTs overlapped between the APO and momilactone B treatments (Fig. 4B). GSTUs 1, 4, 7, 8, 22, 24, and 25 were up-regulated by both APO and momilactone B, while tau class GSTUs 9, 11, and 12, phi class GSTs GSTF 6, and 7, and lambda class GSTL1 were up-regulated only in response to momilactone B. Up-regulated in response to APO were tau class GSTUs 2, 6, 17, and 19, GSTF8 from the phi class, and DHAR2 from the dehydroascorbate reductase clade. Most of the shared GSTs were also found to be up-regulated in response to fenclorim and BOA (Fig. 4B).
The third and final phase of the detoxification process involves transport of the conjugated compounds into the vacuole or the apoplast. We found that treatment with APO and momilactone B both up-regulated transcripts coding for ABC-type transporters (Table 1), namely AT2G47000, AT1G02520, AT3G47780, and AT3G59140, all of which were reported to be up-regulated in the general detoxification process , which might suggest the sequestration of conjugated forms of these compounds into the vacuole or the apoplast.
Genes encoding plant cell wall constituents are downregulated in response to allelochemicals
Downregulated genes in cluster A10 were enriched for genes involved in cell wall integrity (structural constituent of cell wall and xyloglucan glycosylases) or in the peroxisome (heme binding and purine transport). The downregulated genes not assigned to a cluster in the momilactone B samples contained hydrolases and glycosylases (Supplemental Fig. 2).
We treated A. thaliana with the two chemically very distinct, agriculturally relevant allelochemicals APO and momilactone B to understand the transcriptional response to these compounds. Both compounds are produced by Poaceae species; BX are produced by wheat, maize and rye, while momilactones are produced by rice. The genes involved in both pathways are localized in biosynthetic gene clusters [17, 35]. Compounds appear to be produced in a species-specific manner; to date, no species has been identified that is capable of producing both classes of allelochemicals. The genome of Echinochloa crus-galli, a notorious weed in rice fields, contains biosynthetic gene clusters for both BX and momilactones . However, neither has the presence of momilactone B in E. crus-galli been reported in the literature nor were we able to detect it by liquid chromatography–mass spectrometry (LC–MS) analysis (Supplemental Fig. 3). That these two biosynthetic pathways seem to be mutually exclusive in grasses may indicate that grasses produce either BX or momilactone B, but not both. One can therefore speculate that the two compound classes fulfill similar biological functions, and that there is no benefit in producing both. Although the two compounds have very different chemical structures, there was a significant overlap of genes up-regulated in response to both treatments, noticeably of genes that could be attributed the general detoxification response. The fact that up-regulated genes in our dataset overlapped with those from two published A. thaliana transcriptomic analyses in response to the herbicide safener fenclorim  and the BX BOA  (Figs. 3 and 4) further supported the notion that allelochemicals evoke the general detoxification response.
Among the genes up-regulated by either treatment or by both, UGTs, CYPs, and GSTs were overrepresented. One ongoing question about especially UGTs and CYPs involved in detoxification of herbicides and xenobiotics is whether their activity is specific for detoxification or if they usually have a different in planta role but are promiscuous and thus can act on both endogenous and exogenous substrates. With more genes from these families being characterized and their functions being identified, evidence points towards the latter [37, 38]. Our data further support this notion, even though further investigations are necessary to show that up-regulation of these genes in response to BX and momilactone B is accompanied by the according chemical modifications to these substrates.
CYPs involved in diterpene metabolism typically belong to the CYP71 and CYP85 clans; some can be found in the CYP72 clan , including those that were up-regulated by momilactone B. Momilactone B also up-regulated expression of CYP710A1 from the CYP710 clan. Only few CYPs involved in diterpenoid metabolism have been identified in Arabidopsis, and the vast majority of CYP functions is still unknown. CYP708A2 involved in thalianol to 7ß-hydroxythalianol biosynthesis  was up-regulated upon APO treatment, while several CYPs involved in camalexin and glucosinolate biosynthesis [41, 42] were enriched in the momilactone B up-regulated modules. This supports potentially significant enzyme promiscuity among CYPs that could explain their role in xenobiotic detoxification. The CYPs induced by both APO and momilactone B, as well as by fenclorim and BOA were all CYP81 family members. CYP81s have been linked to stress resistance and xenobiotic detoxification in A. thaliana and other plants [32,33,34, 43, 44].
UGT73Bs are consistently up-regulated upon allelochemical treatment (Fig. 4). UGT73B3 and UGT73B5 are also involved in A. thaliana redox response to pathogenic Pseudomonas pathogens [45, 46]. UGT73Bs up-regulated upon fenclorim treatment were shown to glycosylate a number of different xenobiotics in vitro, indicating that these UGTs play a direct role in protection against such compounds . Momilactone B and fenclorim treatments also up-regulated UGT73Cs; among these, UGT73C5 encodes a known brassinosteroid-O-glycosyltransferase  that can also glycosylate and detoxify the mycotoxin deoxynivalenol , while UGT73C6 encodes a UDP-glucose:flavonol-3-O-glycoside-7-O-glucosyltransferase .
Momilactone B is a diterpenoid, so we expected UGTs known to glycosylate diterpenoids to be up-regulated. These include UGTs from families 73, 74, 75, 76 and 85 . Apart from the UGT73s mentioned above, several of the genes up-regulated by momilactone B belong to this UGT family. However, many of these UGTs were also induced by fenclorim treatment  (Fig. 3A), which means that their transcriptional activation is most likely not a diterpenoid-specific response.
The majority of GSTs up-regulated in our study belong to the plant specific tau- and phi-classes that are known to be responsible for glutathione mediated detoxification of xenobiotics . Of the seven tau-class GSTs (GSTU1, 4, 7, 8, 22, 24 and 25) up-regulated by both APO and momilactone B treatment, all except GSTU4, 9 and 24 were also activated by fenclorim and BOA treatment (Fig. 4B), supporting a robust detoxification response of these GSTs. Arabidopsis tau- and phi-class GSTs heterologously expressed in E. coli [51, 52] and in yeast  showed widespread ability to conjugate different herbicides to GSH. GSTLs and DHAR contain a cysteine instead of the active serine residue usually found in the catalytic site of GSTs, and presumably have no GSH conjugating activity , but can be involved in recycling specialized metabolites . GSTL1 and DHAR2 were up-regulated upon momilactone B, and upon both momilactone B and APO treatment, respectively. Since they likely do not conjugate GSH, they may rather be involved in recycling GSH  in response to allelochemical treatment. GSTs have been known to be up-regulated by xenobiotics without playing a direct role in their conjugation , and part of the GST response to allelochemicals may be based on their general antioxidant function related to environmental stress . Whether the GSTs up-regulated in our data conjugate allelochemicals to GSH therefore remains to be determined.
The final step in the three-phase detoxification is the transport of conjugated compounds into the vacuole or apoplast. All of the four ABC transporters up-regulated by both APO and momilactone B treatment were reported to be up-regulated in the general detoxification process , and AT2g4700 and AT3G59140 were further found to be up-regulated in response to methanol toxicity  and BOA treatment . Since their discovery in plants , the MRP-class of ABC transporters has been associated with transport of GSH-conjugated xenobiotics , although MRP transporters are now known to transport a variety of substrates . In our data, MRP14 (AT3G59140) was up-regulated in response to both APO and momilactone B treatment; MRP2 (AT2g34660) was up-regulated by APO, and MRP12 (AT1G30420), MRP9 (AT3G13090), MRP7 (AT3G13100), and MRP3 (AT3g13080) by momilactone B. MRP3 transports GSH-conjugated xenobiotics and chlorophyll catabolites , but no specific transport activity has been reported for the other five MRP transporters. Nonetheless, the transcriptional activation of MRP transporters supports our notion that allelochemicals, like herbicides, are detoxified through the three-phase detoxification system.
In summary, our data show that chemically diverse, phytotoxic compounds that are employed in inhibitory plant-plant interaction can trigger similar detoxification responses in the target plant, and provide insights into possible mechanisms of allelochemical tolerance. Further studies are necessary to investigate the metabolic processes and the chemical nature of potential conjugates. The detoxification of possibly conjugated allelochemicals by plant transporters may also explain how herbicide-responsive detoxification systems are maintained in plant populations that are typically not exposed to synthetic herbicides: since plants have to cope with a wide array of compounds released by other organisms, they have evolved and retained an arsenal of promiscuous enzymes that are able to detoxify harmful molecules with a certain degree of agnosticism towards their origin.
Material and methods
Sample preparations and transcriptome sequencing
A. thaliana grown hydroponically in ½ Murashige and Skoog (MS) media (pH 5.8) without sugar for 3 weeks (APO) or 6 days (momilactone B) were treated with either 3.5 µM APO, 4 µM momilactone B, or the equivalent concentration of the solvent DMSO. Tissue from four biological replicates per treatment (minimum 20 pooled seedlings per replicate) was collected 1, 6, and 24 h after beginning of the treatment and flash-frozen in liquid nitrogen. The tissue was lysed using a Retsch Mill (Retsch) before extracting total RNA using the RNeasy Plant Mini kit (Qiagen). RNA yield was measured using Qubit (Invitrogen) and RNA integrity was confirmed by agarose gel electrophoresis. Barcoded mRNA libraries were generated using either TruSeq RNA Sample Prep kit v2 (Illumina) (APO) or NEBNext RNA Ultra II Directional Library Kit (New England Biolabs) (momilactone B) following the manufacturers’ instructions. Libraries were sequenced on an Illumina HiSeqV4 as 100 bp single-end reads (APO) or an Illumina NovaSeq as 150 bp paired-end reads (momilactone B).
Read mapping and quantification
All reads were mapped to the A. thaliana Col-0 TAIR10 reference genome (arabidopsis.org) (Supplemental Table 1). Mapping and feature counting were done using the nf-core pipeline RNAseq v3.5 [64, 65] with default parameters. In brief, reads were mapped using STAR and quantified using salmon [66, 67]. The counts per gene were further analyzed using R v4.1.2.
Data analysis was carried out using R v4.1.2 and is documented in (https://github.com/nschan/Knoch_et_al_transcriptomes). In brief, count tables were imported into R. Differential expression analysis was carried out using DESeq2 , weighted gene correlation network analysis was performed using WGCNA , and the beta parameter for WGCNA was picked automatically using CEMiTool . Over-representation analysis of GO terms was performed using the enricher function from the clusterProfiler package . Only genes that were contained in the clusters A3, A6, A7 or M1 and had a log2FoldChange > 0 and adjusted p-value below 0.01 at any time-point in the differential analysis were included in the overrepresentation analysis for upregulated genes, while the analysis of downregulated genes included genes contained in cluster A10, or not assigned to a cluster but included in the WGCNA analysis (momilactone B) with a log2FoldChange < 0 and an adjusted p-value below 0.01 in the differential expression analysis.
For phylogenetic analysis, protein sequences were aligned using the AlignSeqs function from DECIPHER , adjusted using the AdjustAlignment function, and maximum likelihood trees were fitted and bootstrapped 1000 times using the phangorn package  using the WAG substitution matrix.
Momilactone B dose response
A. thaliana Col-0 seeds were sterilized with chlorine gas for 1 h and stratified for 6 d in the dark at 4˚C. 20 seeds per plate were sown on ½ MS media supplemented with various concentrations of momilactone B (0 µM, 0.1 µM, 1 µM, 4 µM, 10 µM and 20 µM) dissolved in DMSO, to a final DMSO concentration of 0.1% in all plates. Seedlings were grown in a 16 h/8 h light/dark cycle chamber at 21˚C with a light intensity of 50 µM/m2/sec. After 5 d of growth, the seedlings were imaged using a fixed camera and a ruler for scale. The primary root length was traced using ImageJ and primary root length was calculated using the scale as a reference. Primary root length was plotted as relative percentage of growth compared to the control sample. The drc package in RStudio  was used to fit a dose–response model and to calculate the half-maximal-effect concentration; data was plotted using the ggplot2 package .
Momilactone content in Oryza sativa cv. Kitaake and Echinochloa crus-galli
Oryza sativa cv. Kitaake and Echinochloa crus-galli were grown in sand in the greenhouse for three weeks. Roots were harvested, weighed, snap-frozen in liquid nitrogen and homogenized. Metabolites were extracted in 10 times MeOH to sample weight, and momilactone A and B were measured on a Vanquish HPLC (Thermo Fisher Scientific) coupled via electrospray ionization to an TSQ Altis (Thermo Fisher Scientific) mass spectrometer, and quantified using authentic standards provided by Kazunori Okada, University of Tokyo .
Availability of data and materials
Sequencing reads have been deposited in the ENA Short Read Archive (www.ebi.ac.uk/ena/) under study accession number PRJEB51016.
Analysis code is available at https://github.com/nschan/Knoch_et_al_transcriptomes
Der MH. Einfluss einer Pflanze auf die Andere. Allelopathie Nature. 1938;141:493–493.
Rice E. Allelopathy. 2nd edition. Orlando, Florida: Academic Press; 1983.
Niculaes C, Abramov A, Hannemann L, Frey M. Plant Protection by Benzoxazinoids—Recent Insights into Biosynthesis and Function. Agronomy. 2018;8:143.
Niemeyer HM. Hydroxamic Acids Derived from 2-Hydroxy-2H-1,4-Benzoxazin-3(4H)-one: Key Defense Chemicals of Cereals. J Agric Food Chem. 2009;57:1677–96.
Sicker D, Frey M, Schulz M, Gierl A. Role of natural benzoxazinones in the survival strategy of plants. Int Rev Cytol. 2000;198:319–46.
Kato-Noguchi H, Peters RJ. The Role of Momilactones in Rice Allelopathy. J Chem Ecol. 2013;39:175–85.
Chung IM, Jung TK, Kim SH. Evaluation of allelopathic potential and quantification of momilactone A, B from rice hull extracts and assessment of inhibitory bioactivity on paddy field weeds. J Agric Food Chem. 2006;54:2527–36.
Kariya K, Ube N, Ueno M, Teraishi M, Okumoto Y, Mori N, et al. Natural variation of diterpenoid phytoalexins in cultivated and wild rice species. Phytochemistry. 2020;180:112518.
Krogh SS, Mensz SJM, Nielsen ST, Mortensen AG, Christophersen C, Fomsgaard IS. Fate of Benzoxazinone Allelochemicals in Soil after Incorporation of Wheat and Rye Sprouts. J Agric Food Chem. 2006;54:1064–74.
Macias FA, Marín D, Oliveros-Bastidas A, Castellano D, Simonet AM, Molinillo JMG. Structure-activity relationship (SAR) studies of benzoxazinones, their degradation products, and analogues. Phytotoxicity on problematic weeds Avena fatua L. and Lolium rigidum Gaud. J Agric Food Chem. 2006;54:1040–8.
Atwal AS, Teather RM, Liss SN, Collins FW. Antimicrobial activity of 2-aminophenoxazin-3-one under anaerobic conditions. Can J Microbiol. 1992;38:1084–8.
Venturelli S, Belz RG, Kämper A, Berger A, von Horn K, Wegner A, et al. Plants Release Precursors of Histone Deacetylase Inhibitors to Suppress Growth of Competitors. Plant Cell. 2015;27:3175–89.
Serra Serra N, Shanmuganathan R, Becker C. Allelopathy in rice: a story of momilactones, kin recognition, and weed management. J Exp Bot. 2021;72:4022–37.
Kato-Noguchi H, Hasegawa M, Ino T, Ota K, Kujime H. Contribution of momilactone A and B to rice allelopathy. J Plant Physiol. 2010;167:787–91.
Hasegawa M, Mitsuhara I, Seo S, Imai T, Koga J, Okada K, et al. Phytoalexin Accumulation in the Interaction Between Rice and the Blast Fungus. Mol Plant-Microbe Interactions. 2010;23:1000–11.
Kato-Noguchi H, Ota K, Kujime H. Absorption of momilactone A and B by Arabidopsis thaliana L. and the growth inhibitory effects. J Plant Physiol. 2012;169:1471–6.
Shimura K, Okada A, Okada K, Jikumaru Y, Ko K-W, Toyomasu T, et al. Identification of a Biosynthetic Gene Cluster in Rice for Momilactones. J Biol Chem. 2007;282:34013–8.
Wang Q, Hillwig ML, Okada K, Yamazaki K, Wu Y, Swaminathan S, et al. Characterization of CYP76M5–8 Indicates Metabolic Plasticity within a Plant Biosynthetic Gene Cluster. J Biol Chem. 2012;287:6159–68.
De La Peña R, Sattely ES. Rerouting plant terpene biosynthesis enables momilactone pathway elucidation. Nat Chem Biol. 2021;17:205–12.
Kitaoka N, Zhang J, Oyagbenro RK, Brown B, Wu Y, Yang B, et al. Interdependent evolution of biosynthetic gene clusters for momilactone production in rice. Plant Cell. 2021;33(2):290–305.
Bais HP, Weir TL, Perry LG, Gilroy S, Vivanco JM. The role of root exudates in rhizosphere interactions with plants and other organisms. Annu Rev Plant Biol. 2006;57:233–66.
Sandermann H. Plant metabolism of xenobiotics. Trends Biochem Sci. 1992;17:82–4.
Coleman J, Blake-Kalff M, Davies E. Detoxification of xenobiotics by plants: chemical modification and vacuolar compartmentation. Trends Plant Sci. 1997;2:144–51.
Nguyen T, Sherratt PJ, Pickett CB. Regulatory Mechanisms Controlling Gene Expression Mediated by the Antioxidant Response Element. Annu Rev Pharmacol Toxicol. 2003;43:233–60.
Kang J, Park J, Choi H, Burla B, Kretzschmar T, Lee Y, et al. Plant ABC Transporters. Arab Book Am Soc Plant Biol. 2011;9:e0153.
Venturelli S, Petersen S, Langenecker T, Weigel D, Lauer UM, Becker C. Allelochemicals of the phenoxazinone class act at physiologically relevant concentrations. Plant Signal Behav. 2016;11:e1176818.
Berardini TZ, Reiser L, Li D, Mezheritsky Y, Muller R, Strait E, et al. The arabidopsis information resource: Making and mining the “gold standard” annotated reference plant genome. Genesis. 2015;53:474–85.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.
Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J R Stat Soc Ser B Methodol. 1995;57:289–300.
Siminszky B. Plant cytochrome P450-mediated herbicide metabolism. Phytochem Rev. 2006;5:445–58.
Baerson SR, Sánchez-Moreiras A, Pedrol-Bonjoch N, Schulz M, Kagan IA, Agarwal AK, et al. Detoxification and transcriptome response in Arabidopsis seedlings exposed to the allelochemical benzoxazolin-2(3H)-one. J Biol Chem. 2005;280:21867–81.
Brazier-Hicks M, Gershater M, Dixon D, Edwards R. Substrate specificity and safener inducibility of the plant UDP-glucose-dependent family 1 glycosyltransferase super-family. Plant Biotechnol J. 2018;16:337–48.
Mueller S, Hilbert B, Dueckershoff K, Roitsch T, Krischke M, Mueller MJ, et al. General Detoxification and Stress Responses Are Mediated by Oxidized Lipids through TGA Transcription Factors in Arabidopsis. Plant Cell. 2008;20:768–85.
Frey M, Chomet P, Glawischnig E, Stettner C, Grün S, Winklmair A, et al. Analysis of a Chemical Plant Defense Mechanism in Grasses. Science. 1997;277:696–9.
Guo L, Qiu J, Ye C, Jin G, Mao L, Zhang H, et al. Echinochloa crus-galli genome analysis provides insight into its adaptation and invasiveness as a weed. Nat Commun. 2017;8:1031.
Meßner B, Thulke O, Schäffner AR. Arabidopsis glucosyltransferases with activities toward both endogenous and xenobiotic substrates. Planta. 2003;217:138–46.
Poppenberger B, Fujioka S, Soeno K, George GL, Vaistij FE, Hiranuma S, et al. The UGT73C5 of Arabidopsis thaliana glucosylates brassinosteroids. Proc Natl Acad Sci. 2005;102:15253–8.
Bathe U, Tissier A. Cytochrome P450 enzymes: A driving force of plant diterpene diversity. Phytochemistry. 2019;161:149–62.
Castillo DA, Kolesnikova MD, Matsuda SPT. An Effective Strategy for Exploring Unknown Metabolic Pathways by Genome Mining. J Am Chem Soc. 2013;135:5885–94.
Nafisi M, Sønderby IE, Hansen BG, Geu-Flores F, Nour-Eldin HH, Nørholm MHH, et al. Cytochromes P450 in the biosynthesis of glucosinolates and indole alkaloids. Phytochem Rev. 2006;5:331–46.
Pfalz M, Mikkelsen MD, Bednarek P, Olsen CE, Halkier BA, Kroymann J. Metabolic Engineering in Nicotiana benthamiana Reveals Key Enzyme Functions in Arabidopsis Indole Glucosinolate Modification[W]. Plant Cell. 2011;23:716–29.
Pan G, Zhang X, Liu K, Zhang J, Wu X, Zhu J, et al. Map-based cloning of a novel rice cytochrome P450 gene CYP81A6 that confers resistance to two different classes of herbicides. Plant Mol Biol. 2006;61:933–43.
Liu C-J, Huhman D, Sumner LW, Dixon RA. Regiospecific hydroxylation of isoflavones by cytochrome P450 81E enzymes from Medicago truncatula. Plant J. 2003;36:471–84.
Langlois-Meurinne M, Gachon CMM, Saindrenan P. Pathogen-Responsive Expression of Glycosyltransferase Genes UGT73B3 and UGT73B5 Is Necessary for Resistance to Pseudomonas syringae pv tomato in Arabidopsis. Plant Physiol. 2005;139:1890–901.
Simon C, Langlois-Meurinne M, Didierlaurent L, Chaouch S, Bellvert F, Massoud K, et al. The secondary metabolism glycosyltransferases UGT73B3 and UGT73B5 are components of redox status in resistance of Arabidopsis to P seudomonas syringae pv tomato: Glycosyltransferases and redox status. Plant Cell Environ. 2014;37:1114–29.
Poppenberger B, Berthiller F, Lucyshyn D, Sieberer T, Schuhmacher R, Krska R, et al. Detoxification of the Fusarium Mycotoxin Deoxynivalenol by a UDP-glucosyltransferase from Arabidopsis thaliana. J Biol Chem. 2003;278:47905–14.
Jones P, Messner B, Nakajima J-I, Schäffner AR, Saito K. UGT73C6 and UGT78D1, Glycosyltransferases Involved in Flavonol Glycoside Biosynthesis in Arabidopsis thaliana. J Biol Chem. 2003;278:43910–8.
Wilson AE, Tian L. Phylogenomic analysis of UDP-dependent glycosyltransferases provides insights into the evolutionary landscape of glycosylation in plant metabolism. Plant J. 2019;100:1273–88.
Cummins I, Dixon DP, Freitag-Pohl S, Skipsey M, Edwards R. Multiple roles for plant glutathione transferases in xenobiotic detoxification. Drug Metab Rev. 2011;43:266–80.
Dixon DP, Skipsey M, Edwards R. Roles for glutathione transferases in plant secondary metabolism. Phytochemistry. 2010;71:338–50.
Dixon DP, Hawkins T, Hussey PJ, Edwards R. Enzyme activities and subcellular localization of members of the Arabidopsis glutathione transferase superfamily. J Exp Bot. 2009;60:1207–18.
Krajewski MP, Kanawati B, Fekete A, Kowalski N, Schmitt-Kopplin P, Grill E. Analysis of Arabidopsis glutathione-transferases in yeast. Phytochemistry. 2013;91:198–207.
Dixon DP, Davis BG, Edwards R. Functional Divergence in the Glutathione Transferase Superfamily in Plants. J Biol Chem. 2002;277:30859–69.
Bjarnholt N, Neilson EHJ, Crocoll C, Jørgensen K, Motawia MS, Olsen CE, et al. Glutathione transferases catalyze recycling of auto-toxic cyanogenic glucosides in sorghum. Plant J. 2018;94:1109–25.
Rahantaniaina M-S, Li S, Chatel-Innocenti G, Tuzet A, Issakidis-Bourguet E, Mhamdi A, et al. Cytosolic and Chloroplastic DHARs Cooperate in Oxidative Stress-Driven Activation of the Salicylic Acid Pathway. Plant Physiol. 2017;174:956–71.
Rylott EL, Bruce NC. Plants disarm soil: engineering plants for the phytoremediation of explosives. Trends Biotechnol. 2009;27:73–81.
Ezaki B, Gardner RC, Ezaki Y, Matsumoto H. Expression of Aluminum-Induced Genes in Transgenic Arabidopsis Plants Can Ameliorate Aluminum Stress and/or Oxidative Stress. Plant Physiol. 2000;122:657–66.
Downie A, Miyazaki S, Bohnert H, John P, Coleman J, Parry M, et al. Expression profiling of the response of Arabidopsis thaliana to methanol stimulation. Phytochemistry. 2004;65:2305–16.
Martinoia E, Grill E, Tommasini R, Kreuz K, Amrhein N. ATP-dependent glutathione S-conjugate “export” pump in the vacuolar membrane of plants. Nature. 1993;364:247–9.
Liu G, Sánchez-Fernández R, Li Z-S, Rea PA. Enhanced Multispecificity of Arabidopsis Vacuolar Multidrug Resistance-associated Protein-type ATP-binding Cassette Transporter, AtMRP2. J Biol Chem. 2001;276:8648–56.
Klein M, Burla B, Martinoia E. The multidrug resistance-associated protein (MRP/ABCC) subfamily of ATP-binding cassette transporters in plants. FEBS Lett. 2006;580:1112–22.
Tommasini R, Vogt E, Fromenteau M, Hörtensteiner S, Matile P, Amrhein N, et al. An ABC-transporter of Arabidopsis thaliana has both glutathione-conjugate and chlorophyll catabolite transport activity. Plant J. 1998;13:773–80.
Ewels PA, Peltzer A, Fillinger S, Patel H, Alneberg J, Wilm A, et al. The nf-core framework for community-curated bioinformatics pipelines. Nat Biotechnol. 2020;38:276–8.
Ewels PA, Peltzer A, Fillinger S, Patel H, Alneberg J, Wilm A, et al. The nf-core framework for community-curated bioinformatics pipelines. Nat Biotechnol. 2020;38(3):276–8.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.
Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017;14:417–9.
Russo PST, Ferreira GR, Cardozo LE, Bürger MC, Arias-Carrasco R, Maruyama SR, et al. CEMiTool: a Bioconductor package for performing comprehensive modular co-expression analyses. BMC Bioinformatics. 2018;19:56.
Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R Package for Comparing Biological Themes Among Gene Clusters. OMICS J Integr Biol. 2012;16:284–7.
Wright ES. Using DECIPHER v2.0 to Analyze Big Biological Sequence Data in R. R J. 2016;8:352–9.
Schliep KP. phangorn: phylogenetic analysis in R. Bioinformatics. 2011;27:592–3.
Ritz C, Baty F, Streibig JC, Gerhard D. Dose-Response Analysis Using R. PLOS ONE. 2015;10:e0146021.
Pedersen HW Danielle Navarro, and Thomas Lin. ggplot2: Elegant Graphics for Data Analysis. Basel: Springer-Verlag New York; 2009.
Tomita K, Yashiroda Y, Matsuo Y, Piotrowski JS, Li SC, Okamoto R, et al. Genome-wide screening of genes associated with momilactone B sensitivity in the fission yeast Schizosaccharomyces pombe. G3 GenesGenomesGenetics. 2021;11:156.
We would like to thank the Next-Generation Sequencing unit at the Vienna Biocenter Core Facilities (VBCF) for their help in generating the sequencing data, and Katharina Jandrasits for technical lab support. LC-MS/MS analysis was performed by the Metabolomics Facility at Vienna BioCenter Core Facilities (VBCF), a member of the Vienna BioCenter (VBC), Austria, and funded by the City of Vienna through the Vienna Business Agency.
Open Access funding enabled and organized by Projekt DEAL. This research was supported by the Austrian Academy of Sciences (ÖAW) via funding to C.B. and via a DOC fellowship to N.S.S., grant ID 25652; by the Austrian Science Fund (FWF) via a Lise Meitner Fellowship grant to E.K., grant ID FWF 2484-B21; by the Deutsche Forschungsgemeinschaft (DFG) and the Japan Society for the Promotion of Science (JSPS) via a joint grant to C.B. and K.O., grant ID 457739273; and by the European Union’s Horizon 2020 research and innovation programme by the European Research Council (ERC), Grant Agreement No. 716823 ‘FEAR-SAP’.
The authors declare no competing interests.
Ethics approval and consent to participate
The research complied with national and international regulations.
The seeds for A. thaliana Col-0 were obtained from the Nottingham Arabidopsis Stock Centre (NASC) under the standard material transfer agreement (MTA).
Consent for publication
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Momilactone B dose response. (A) Growthphenotype of A. thaliana seedlings grown on ½ MS-agar with varyingconcentrations of momilactone B. Scale bar = 1 cm. (B) Primary root lengthplotted as relative percentage of growth compared to the control sample. Supplemental Figure 2.Overrepresentation analysis of down-regulated genes. Lollipop plot of geneontology (GO) terms of genes included in clusters A5 and A10 (APO) or not partof a cluster (“Not correlated”, momilactone B). Genes with a negative log foldchange and an adjusted p-value < 0.01 were included in an overrepresentationanalysis of GO terms. Orange bars indicate the number of genes belonging to aparticular GO-term relative to the total number of genes belonging to the term,blue bars indicate number of genes belonging to the GO-term compared to thetotal number of genes in the genome. Circle fill color indicates the p-value ofthe hypergeometric test, adjusted for multiple comparisons using the method ofBenjamini-Hochberg . SupplementalFigure 3. Momilactone content in Oryza sativa cv. Kitaake (Rice) andEchinochloa crus-galli (ECG). Momilactone A and B in MeOH extracts from rootsof three week old rice and E. crus-galli were measured by LCMS. Boxplotssummarizing 18 replicates are shown, boxes indicate 1st to 3rd quartile,horizontal line indicates the median, whiskers extend to quartile1-1,5*IQR andquartile3+1,5*IQR. Supplemental Table 1.RNA-seq mapping statistics. SupplementalTable 2. Differentially expressed CYP45s.
About this article
Cite this article
Knoch, E., Kovács, J., Deiber, S. et al. Transcriptional response of a target plant to benzoxazinoid and diterpene allelochemicals highlights commonalities in detoxification. BMC Plant Biol 22, 402 (2022). https://doi.org/10.1186/s12870-022-03780-w
- Momilactone B