Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Letter
  • Published:

Root-secreted bitter triterpene modulates the rhizosphere microbiota to improve plant fitness

Abstract

Underground microbial ecosystems have profound impacts on plant health1,2,3,4,5. Recently, essential roles have been shown for plant specialized metabolites in shaping the rhizosphere microbiome6,7,8,9. However, the potential mechanisms underlying the root-to-soil delivery of these metabolites remain to be elucidated10. Cucurbitacins, the characteristic bitter triterpenoids in cucurbit plants (such as melon and watermelon), are synthesized by operon-like gene clusters11. Here we report two Multidrug and Toxic Compound Extrusion (MATE) proteins involved in the transport of their respective cucurbitacins, a process co-regulated with cucurbitacin biosynthesis. We further show that the transport of cucurbitacin B from the roots of melon into the soil modulates the rhizosphere microbiome by selectively enriching for two bacterial genera, Enterobacter and Bacillus, and we demonstrate that this, in turn, leads to robust resistance against the soil-borne wilt fungal pathogen, Fusarium oxysporum. Our study offers insights into how transporters for specialized metabolites manipulate the rhizosphere microbiota and thereby affect crop fitness.

This is a preview of subscription content, access via your institution

Access options

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: Genome mining for bitterness transporters in cucurbits.
Fig. 2: Characterization of in vitro and in vivo functions of the two cucurbit MATEs.
Fig. 3: Roles of CuB in shaping the rhizosphere microbiome of melon.
Fig. 4: Characteristics of CuB functional activity against wilt disease of melon.

Similar content being viewed by others

Data availability

All data supporting the findings of this study are available in the paper, the extended data and the Supplementary Information. The amplicon sequencing data and shotgun sequencing data were deposited at NCBI under Bioproject No. PRJNA756119 and are publicly accessible at https://www.ncbi.nlm.nih.gov/search/all/?term=PRJNA756119. Source data are provided with this paper.

References

  1. Pieterse, C. M. et al. Induced systemic resistance by beneficial microbes. Annu. Rev. Phytopathol. 52, 347–375 (2014).

    Article  CAS  Google Scholar 

  2. Hacquard, S., Spaepen, S., Garrido-Oter, R. & Schulze-Lefert, P. Interplay between innate immunity and the plant microbiota. Annu. Rev. Phytopathol. 55, 565–589 (2017).

    Article  CAS  Google Scholar 

  3. Berendsen, R. L., Pieterse, C. M. & Bakker, P. A. The rhizosphere microbiome and plant health. Trends Plant Sci. 17, 478–486 (2012).

    Article  CAS  Google Scholar 

  4. Vannier, N., Agler, M. & Hacquard, S. Microbiota-mediated disease resistance in plants. PLoS Pathog. 15, e1007740 (2019).

    Article  CAS  Google Scholar 

  5. de la Fuente Cantó, C. et al. An extended root phenotype: the rhizosphere, its formation and impacts on plant fitness. Plant J. 103, 951–964 (2020).

    Article  Google Scholar 

  6. Stringlis, I. A. et al. MYB72-dependent coumarin exudation shapes root microbiome assembly to promote plant health. Proc. Natl Acad. Sci. USA 115, E5213–E5222 (2018).

    Article  CAS  Google Scholar 

  7. Hu, L. et al. Root exudate metabolites drive plant–soil feedbacks on growth and defense by shaping the rhizosphere microbiota. Nat. Commun. 9, 2738 (2018).

    Article  Google Scholar 

  8. Huang, A. C. et al. A specialized metabolic network selectively modulates Arabidopsis root microbiota. Science 364, eaau6389 (2019).

    Article  CAS  Google Scholar 

  9. Koprivova, A. et al. Root-specific camalexin biosynthesis controls the plant growth-promoting effects of multiple bacterial strains. Proc. Natl Acad. Sci. USA 116, 15735–15744 (2019).

    Article  CAS  Google Scholar 

  10. Sasse, J., Martinoia, E. & Northen, T. Feed your friends: do plant exudates shape the root microbiome? Trends Plant Sci. 23, 25–41 (2018).

    Article  CAS  Google Scholar 

  11. Zhou, Y. et al. Convergence and divergence of bitterness biosynthesis and regulation in Cucurbitaceae. Nat. Plants 2, 16183 (2016).

    Article  CAS  Google Scholar 

  12. Augustin, J. M., Kuzina, V., Andersen, S. B. & Bak, S. Molecular activities, biosynthesis and evolution of triterpenoid saponins. Phytochemistry 72, 435–457 (2011).

    Article  CAS  Google Scholar 

  13. Osbourn, A., Goss, R. J. & Field, R. A. The saponins—polar isoprenoids with important and diverse biological activities. Nat. Prod. Rep. 28, 1261–1268 (2011).

    Article  CAS  Google Scholar 

  14. Thimmappa, R., Geisler, K., Louveau, T., O’Maille, P. & Osbourn, A. Triterpene biosynthesis in plants. Annu. Rev. Plant Biol. 65, 225–257 (2014).

    Article  CAS  Google Scholar 

  15. Field, B. & Osbourn, A. E. Metabolic diversification—independent assembly of operon-like gene clusters in different plants. Science 320, 543–547 (2008).

    Article  CAS  Google Scholar 

  16. Kazachkova, Y. et al. The GORKY glycoalkaloid transporter is indispensable for preventing tomato bitterness. Nat. Plants 7, 468–480 (2021).

    Article  CAS  Google Scholar 

  17. Qi, X. et al. A different function for a member of an ancient and highly conserved cytochrome P450 family: from essential sterols to plant defense. Proc. Natl Acad. Sci. USA 103, 18848–18853 (2006).

    Article  CAS  Google Scholar 

  18. Mugford, S. T. et al. A serine carboxypeptidase-like acyltransferase is required for synthesis of antimicrobial compounds and disease resistance in oats. Plant Cell 21, 2473–2484 (2009).

    Article  CAS  Google Scholar 

  19. Boisson-Dernier, A. et al. Agrobacterium rhizogenes-transformed roots of Medicago truncatula for the study of nitrogen-fixing and endomycorrhizal symbiotic associations. Mol. Plant Microbe Interact. 14, 695–700 (2001).

    Article  CAS  Google Scholar 

  20. Macias, F. A. et al. Structure–activity relationships (SAR) studies of benzoxazinones, their degradation products and analogues: phytotoxicity on standard target species (STS). J. Agric. Food Chem. 53, 538–548 (2005).

    Article  CAS  Google Scholar 

  21. Ongena, M. & Jacques, P. Bacillus lipopeptides: versatile weapons for plant disease biocontrol. Trends Microbiol. 16, 115–125 (2008).

    Article  CAS  Google Scholar 

  22. Niu, B., Paulson, J. N., Zheng, X. & Kolter, R. Simplified and representative bacterial community of maize roots. Proc. Natl Acad. Sci. USA 114, E2450–E2459 (2017).

    Article  CAS  Google Scholar 

  23. Baetz, U. & Martinoia, E. Root exudates: the hidden part of plant defense. Trends Plant Sci. 19, 90–98 (2014).

    Article  CAS  Google Scholar 

  24. Yoo, S.-D., Cho, Y.-H. & Sheen, J. Arabidopsis mesophyll protoplasts: a versatile cell system for transient gene expression analysis. Nat. Protoc. 2, 1565–1572 (2007).

    Article  CAS  Google Scholar 

  25. Li, X. et al. Single-molecule analysis of PIP2;1 dynamics and partitioning reveals multiple modes of Arabidopsis plasma membrane aquaporin regulation. Plant Cell 23, 3780–3797 (2011).

    Article  CAS  Google Scholar 

  26. Chen, H. et al. An ACC oxidase gene essential for cucumber carpel development. Mol. Plant 9, 1315–1327 (2016).

    Article  CAS  Google Scholar 

  27. Lincoln, C., Long, J., Yamaguchi, J., Serikawa, K. & Hake, S. A knotted1-like homeobox gene in Arabidopsis is expressed in the vegetative meristem and dramatically alters leaf morphology when overexpressed in transgenic plants. Plant Cell 6, 1859–1876 (1994).

    CAS  PubMed  PubMed Central  Google Scholar 

  28. Tian, S. et al. Efficient CRISPR/Cas9-based gene knockout in watermelon. Plant Cell Rep. 36, 399–406 (2017).

    Article  CAS  Google Scholar 

  29. Pak, H.-K. et al. Hairy root induction in Oriental melon (Cucumis melo) by Agrobacterium rhizogenes and reproduction of the root-knot nematode (Meloidogyne incognita). Plant Cell Tiss. Org. Cult. 98, 219–228 (2009).

    Article  CAS  Google Scholar 

  30. Kajikawa, M., Morikawa, K., Abe, Y., Yokota, A. & Akashi, K. Establishment of a transgenic hairy root system in wild and domesticated watermelon (Citrullus lanatus) for studying root vigor under drought. Plant Cell Rep. 29, 771–778 (2010).

    Article  CAS  Google Scholar 

  31. Zhou, X. & Wu, F. p-Coumaric acid influenced cucumber rhizosphere soil microbial communities and the growth of Fusarium oxysporum f. sp. cucumerinum Owen. PLoS ONE 7, e48288 (2012).

    Article  CAS  Google Scholar 

  32. Wickham, H. Getting Started with ggplot2. In: ggplot2. Use R!. (Springer, Cham, 2016).

  33. Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).

    Article  CAS  Google Scholar 

  34. Garcia-Mas, J. et al. The genome of melon (Cucumis melo L.). Proc. Natl Acad. Sci. USA 109, 11872–11877 (2012).

    Article  CAS  Google Scholar 

  35. Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359 (2012).

    Article  CAS  Google Scholar 

  36. Menzel, P., Ng, K. L. & Krogh, A. Fast and sensitive taxonomic classification for metagenomics with Kaiju. Nat. Commun. 7, 11257 (2016).

    Article  CAS  Google Scholar 

  37. Klein, M. et al. The ATP-binding cassette (ABC) transporter Bpt1p mediates vacuolar sequestration of glutathione conjugates in yeast. FEBS Lett. 520, 63–67 (2002).

    Article  CAS  Google Scholar 

  38. Yamanaka, N., Marqués, G. & O’Connor, M. B. Vesicle-mediated steroid hormone secretion in Drosophila melanogaster. Cell 163, 907–919 (2015).

    Article  CAS  Google Scholar 

  39. Zhao, J. & Dixon, R. A. MATE transporters facilitate vacuolar uptake of epicatechin 3′-O-glucoside for proanthocyanidin biosynthesis in Medicago truncatula and Arabidopsis. Plant Cell 21, 2323–2340 (2009).

    Article  CAS  Google Scholar 

  40. Zhang, J. et al. NRT1.1B is associated with root microbiota composition and nitrogen use in field-grown rice. Nat. Biotechnol. 37, 676–684 (2019).

    Article  CAS  Google Scholar 

  41. Xu, Z. et al. Contribution of bacillomycin D in Bacillus amyloliquefaciens SQR9 to antifungal activity and biofilm formation. Appl. Environ. Microbiol. 79, 808–815 (2013).

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We thank Q. R. Shen (Nanjing Agricultural University), H. S. Wang and S. L. Wang (Institute of Vegetables and Flowers, CAAS) for providing experimental assistance. This work was supported by the Yunnan provincial, Shenzhen municipal and Dapeng district governments. This work was funded by the National Key Research and Development Program of China (grant no. 2018YFA0901800 to Y.S.), the Yunnan Science Fund (grant no. 202105AF150028, 202005AE160015 and 2019FJ004 to Y.S.), the Shenzhen Science and Technology Program (grant no. KQTD2016113010482651 to S.H.), the Special Funds for Science Technology Innovation and Industrial Development of Shenzhen Dapeng New District (grant no. RC201901-05 to S.H.), the BBSRC Institute Strategic Programme Grant ‘Molecules from Nature–Products and Pathways’ (grant no. BBS/E/J/000PR9790 to A.O.) and the John Innes Foundation.

Author information

Authors and Affiliations

Authors

Contributions

Y.S. and Y.M. conceived and designed the project. Y.M., Y. Zhong, X.W., S.T., D.L., Y. Zhou, Y.Q., B.Z., G.Z., X.C. and H.C. performed the experiments. W.X. and Y. Zhang performed the bioinformatics analysis on the microbiota sequencing results. Y.M., W.X., X.W., Y. Zhong and Y.S. analysed the data. Y.M., W.X., Y.L., L.L., A.O., W.J.L., Y.S. and S.H. wrote the paper.

Corresponding authors

Correspondence to Yongshuo Ma or Yi Shang.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature Plants thanks Reuben Peters and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data

Extended Data Fig. 1 Distribution of cucurbitacin and expression profiles of the candidate genes in the cultivar and wild relative.

For melon, gene expression is presented by FPKM value of the candidates; for watermelon, gene expression was measured by qRT-PCR assay. Data are represented as mean ± SD (n = 3, biological independent samples). Relative gene expression level is shown with the same scale. Potential bitterness transporters are indicated in bold font.

Source data

Extended Data Fig. 2 Interaction between the promoter of two candidate cucurbitacin transporter genes and the bitterness regulators of melon and watermelon.

a and d, Yeast one-hybrid analysis of the interaction between the melon bitterness regulator, CmBl or CmBt, and the promoter of CmMATE1 (a), and the interaction between watermelon bitterness regulator, ClBl or ClBt, and the promoter of ClMATE1 (d). P, promoter; AD, activation domain; 3-AT, 3-amino-1,2,4-triazole; SD-L-T, synthetic dropout media W/O leucine and tryptophan; SD-L-T-H, synthetic dropout media W/O leucine, tryptophan, and histidine. b and e, CmBl or CmBt, binds to the promoter region of CmMATE1 (b) and activates a higher luciferase activity relative to the control (empty vector) in N. benthamiana leaves. Similar result was also observed for ClBl or ClBt, and ClMATE1 (e). LUC, luciferase. Renilla LUC was used as reference. Data are represented as mean ± SD (n = 3, biological independent experiments). Statistical differences were analysed using two-sided Student’s t-test. ** P < 0.01. *** P < 0.001. c and f, EMSA shows that CmBl or CmBt was recruited to the promoter of CmMATE1 (c). Similar result was also observed for ClBl and ClMATE1, or ClBt and ClMATE1 (f). Promoter region from -2,000 upstream from the start codon is schematically presented. Black vertical lines indicate the location of an E-box and red horizontal lines indicate the region amplified for EMSA. Comp, competitor (unlabeled probe); His, His tag; +/-, presence/absence of protein or competitor; closed triangle, increasing amount of protein or competitor.

Source data

Extended Data Fig. 3 Subcellular and tissue localization of the transport and biosynthetic genes.

a, Subcellular localization of the two MATEs in cucumber protoplasts. CmMATE1-GFP or ClMATE1-GFP fluorescence is co-localized with a known Arabidopsis plasma membrane marker PIP2A-RFP fluorescence. From left to right: GFP signal (green), RFP signal (red), Differential Interference Contrast (DIC), and overlay of three signals (green, red and DIC) of the same protoplast. b and c, In situ localization of CmBi (b) or ClBi (c) transcripts in the outer region of young root tips in melon and watermelon, respectively. Both transverse (above the broken line) and longitudinal (below the broken line) sections of the root tip are presented. d, Subcellular localization of CmMATE1-GFP or ClMATE1-GFP fusion protein in yeast. All MATE-GFP fusion proteins were localized to the yeast tonoplast. From left to right: GFP fluorescence (green), FM4-64 staining of the vacuolar membrane (red), Differential Interference Contrast (DIC), merged images with GFP, FM4-64 and DIC from the same sample. For a to d, Three independent experiments were repeated with similar results.

Extended Data Fig. 4 Phenotypic analysis of melon and watermelon seedlings.

a, CuB level accumulated in the rhizosphere of soils planted with Cm104 or Cm70 treated with different levels of CuB solution. Data are represented as mean ± SD (n = 3). DW, dry weight of soil; ND, not detected. b and c, The growth of Cm70 and Cm104 was normal when planted in pots containing the steriled conducive soil after 8 days (b). However, in the non-sterile conducive soil, most of the Cm70 seedlings exhibited wilting symptoms and eventually died, whereas most Cm104 seedlings grew normally (c). Cm104, the melon with a bitter root; Cm70, the melon line with a non-bitter root. d, Progress of wilt disease on the Cm70 plants treated with CuC or CuE. Data are represented as mean ± SD for three independent experiments (n = 48 in total for each treatment). e, The growth of ClMATE1-edited and its wild-type (WT) watermelon plants was normal when planted on the same non-sterile conducive soil, after 8 days. For b, c, e, Three independent experiments were repeated with similar results.

Source data

Extended Data Fig. 5 Shannon diversity of the root-associated microbiome.

a and b, The bacterial (a) and fungal (b) diversity index of the soil, rhizosphere and endosphere communities of healthy and diseased seedlings (means ± SD, n = 3). The bacterial Shannon diversity was significantly lower in the healthy lines than the diseased lines, in the rhizosphere, but not in the endosphere communities, whereas the fungal Shannon diversity showed no significant changes between the two groups, regardless of rhizosphere or endosphere. Statistical differences were analysed using two-sided Student’s t-test. ***, significant difference (P < 0.001). n.s., Not significant.

Source data

Extended Data Fig. 6 Principal coordinates analysis of bacterial communities from the rhizosphere of healthy and diseased seedlings.

The three genera were selected as indicators correlating to plant disease susceptibility, and their significance was detected by the Mantel test.

Extended Data Fig. 7 The relationship between relative abundance of Fusarium and Bacillus (a) as well as Fusarium and Enterobacter (b).

The grey area along the solid lines shows the 95% confidence bands of the linear regression.

Extended Data Fig. 8 The relative abundance of detected Bacillus (a) and Enterobacter (b) species and F. oxysporum (c) in the rhizosphere of the Cm104, healthy and diseased Cm70 seedlings as determined by shotgun sequencing (n = 3).

Unpaired two-tail Student’s T-test was used to detect the difference of the relative abundance of species in the rhizosphere between group levels. *P < 0.05, **P < 0.01, *** P < 0.001. ns: not significant. The box plots show the first and third quartiles (boxes) and the median (middle line); whiskers show inner fences (that is 1.5 times the interquartile range below the first quartile and above the third quartile); points in the boxes show the mean values.

Source data

Extended Data Fig. 9 Growth profiles of the Bacillus (a) or Enterobacter (b) strains cultivated in TSB medium supplemented with CuB.

Note that growth of the Bacillus strains could be significantly inhibited in TSB medium supplemented with 150 ppm CuB, whereas CuB had little effect on the growth of the Enterobacter strains E-126 and E-128. Data are represented as mean ± SD (n = 6 for a and n = 8 for b).

Source data

Extended Data Fig. 10 CuB potentially provide a carbon source for promotion of Enterobacter growth.

a, The profiles of CuB content in TSB medium after 3 days of culture with the indicated Enterobacter strains. A significant decrease in CuB level was observed in the medium cultured with Enterobacter strains, compared to that of TSB medium without Enterobacter strain (CK). Data are represented as mean ± SD (n = 3). b and c, HPLC-QTOF-MS analysis of the E-159 medium supplemented with CuB. b, Metabolic profiling analysis of extracts prepared from E-159 medium. The medium supplemented with only CuB (CuB + medium) or without CuB but containing E-159 strain (E-159 + medium) was considered as negative control. CuB was dissolved in methanol, and equal-volume methanol was added for each treatment. Two specific peaks (indicated by CuB derivate), at retention time of 2.64 and 2.72 min, respectively, were detected in the E-159 medium supplemented with CuB. TIC, total ion chromatograms; EIC 558.3421, extracted a characteristic fragmental ion of the CuB derivates at m/z of 558.3421. c, Comparison of the mass spectrum between CuB derivates and CuB standard. Although the mass spectra of CuB derivates and CuB shared the same quasi-molecular ion peak at m/z 576, the ion signal in CuB derivates were deciphered as a hydron plus CuB derivate ([M + H]+), while the ion signal in CuB it was suggested to be an ammonium cation plus CuB ([M + NH4]+).

Source data

Supplementary information

Supplementary Information

Supplementary Figs. 1–6.

Reporting Summary

Supplementary Tables

Supplementary Tables 1–3.

Source data

Source Data Fig. 2

Statistical source data.

Source Data Fig. 3

Statistical source data.

Source Data Fig. 4

Statistical source data.

Source Data Extended Data Fig. 1

Statistical source data.

Source Data Extended Data Fig. 2

Statistical source data.

Source Data Extended Data Fig. 2

Unprocessed EMSA.

Source Data Extended Data Fig. 4

Statistical source data.

Source Data Extended Data Fig. 5

Statistical source data.

Source Data Extended Data Fig. 8

Statistical source data.

Source Data Extended Data Fig. 9

Statistical source data.

Source Data Extended Data Fig. 10

Statistical source data.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhong, Y., Xun, W., Wang, X. et al. Root-secreted bitter triterpene modulates the rhizosphere microbiota to improve plant fitness. Nat. Plants 8, 887–896 (2022). https://doi.org/10.1038/s41477-022-01201-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41477-022-01201-2

This article is cited by

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing