-
PDF
- Split View
-
Views
-
Cite
Cite
Akifumi Sugiyama, Flavonoids and saponins in plant rhizospheres: roles, dynamics, and the potential for agriculture, Bioscience, Biotechnology, and Biochemistry, Volume 85, Issue 9, September 2021, Pages 1919–1931, https://doi.org/10.1093/bbb/zbab106
- Share Icon Share
ABSTRACT
Plants are in constant interaction with a myriad of soil microorganisms in the rhizosphere, an area of soil in close contact with plant roots. Recent research has highlighted the importance of plant-specialized metabolites (PSMs) in shaping and modulating the rhizosphere microbiota; however, the molecular mechanisms underlying the establishment and function of the microbiota mostly remain unaddressed. Flavonoids and saponins are a group of PSMs whose biosynthetic pathways have largely been revealed. Although these PSMs are abundantly secreted into the rhizosphere and exert various functions, the secretion mechanisms have not been clarified. This review summarizes the roles of flavonoids and saponins in the rhizosphere with a special focus on interactions between plants and the rhizosphere microbiota. Furthermore, this review introduces recent advancements in the dynamics of these metabolites in the rhizosphere and indicates potential applications of PSMs for crop production and discusses perspectives in this emerging research field.

Flavonoids and saponins are secreted into rhizospheres. These plant-specialized metabolites exert various functions in mediating interactions with microbiota.
Plants produce a diverse array of low molecular weight compounds, the number of which exceeds 200 000 (Fernie 2007; Afendi et al.2012). A large portion of these metabolites are not necessarily essential for growth and development but often have biological activities such as protection against biotic and abiotic stress and modulation of interactions with other organisms. These metabolites are called plant-specialized metabolites (PSMs), and it is thought that plants acquired the ability to biosynthesize these metabolites during evolution for adaptation to the ecosystem (Kim and Buell 2015; Lichman, Godden and Buell 2020).
In nature, plants are in constant interaction with a myriad of microorganisms ranging from pathogens, commensals, and beneficial microbes. It is believed that plants have evolved to benefit from beneficial microorganisms while protecting themselves from pathogens; however, how to distinguish these microbes at the molecular level remains unclear. How do plants engage with beneficial microorganisms while at the same time restricting pathogens is among the top unanswered questions about molecular plant–microbe interactions (Harris et al.2020), and answering this question is necessary for the utilization of these microbes for sustainable crop production (Finkel et al.2017; Zhang, Vivanco and Shen 2017). The rhizosphere, the region of soil in proximity to plant roots, is a hotspot for these interactions, and it harbors unique repositories of microbes and metabolites (Hartmann, Rothballer and Schmid 2008). Plant roots exude a substantial amount (up to 40%) of photosynthesis-derived carbons into the rhizosphere (Lynch and Whipps 1990; Badri and Vivanco 2009; Haichar et al.2016). The majority of them are primary metabolites including sugars, amino acids, and organic acids, which nourish rhizosphere microbes (Canarini et al.2019). In addition to the primary metabolites, plant roots secrete PSMs that have important roles in the interactions between plants and soil microbes.
Root-derived metabolites, along with the physicochemical properties of soil and environment factors, influence soil microbial communities and affect the formation of rhizosphere microbiota characterized by abundant and active microbial communities having reduced diversity, as compared with bulk soil (Vieira et al.2020; Wang et al.2020). Recent evidence supports the importance of the rhizosphere microbiota in plant growth, immunity, and fitness, emphasizing the potential for utilization in crop production (Compant et al.2019; Canto et al.2020). The involvement of PSMs in shaping the root/rhizosphere microbiota was demonstrated using plant mutants disrupted in a particular biosynthetic pathway and has been summarized in recent reviews (Pascale et al.2020; Jacoby, Koprivova and Kopriva 2021; Pang et al.2021). In addition to the loss-of-function approach, we have employed an artificial rhizosphere treatment method to reveal the functions of flavonoids and saponins in shaping the rhizosphere microbiota (Fujimatsu et al.2020; Okutani et al.2020; Nakayasu et al.2021) (Figure 1). In this review, we focus on 2 types of PSMs—flavonoids and saponins—that are abundantly secreted from roots to summarize their functions and dynamics in the rhizosphere and highlight future challenges in harnessing rhizosphere microbiota for plant robustness and sustainable crop production.

Soybean rhizosphere. Isoflavones and soyasaponins are secreted into the rhizosphere (①) and affect the bacterial communities (②).
Flavonoids
Synthesis and accumulation
Flavonoids are phenolic compounds comprising more than 8000 distinct molecules (Pietta 2000). Flavonoids confer protection from UV-B radiation, pathogen, and herbivores, attract pollinators, regulate auxin transport, and modulate reactive oxygen species and fertility in plants (Ferreyra, Rius and Casati 2012). Additionally, flavonoids are secreted from roots and exert diverse roles in the rhizosphere (see below). Flavonoids in plants are biosynthesized through the phenylpropanoid and the acetate–malonate pathways to form molecules composed of 15 carbon atoms arranged as C6–C3–C6. Chalcone synthase (CHS) catalyzes the first step in flavonoid biosynthesis using p-coumaroyl-CoA and malonyl-CoA as substrates (Yonekura-Sakakibara, Higashi and Nakabayashi. 2019). Based on the oxidation and substitution of the C (C3) ring, flavonoids are further categorized into the subgroups, including anthocyanins, flavonols, flavanones, flavones, flavanols, isoflavones, chalcones, catechins, and aurones (Panche, Diwan and Chandra 2016). Isoflavones are predominantly found in legume plants (Figure 2). Daidzein and genistein in soybean are biosynthesized from liquiritigenin and naringenin, respectively, by isoflavone synthase (IFS), a P450 protein, to form 2-hydroxyisoflavanone (Akashi, Aoki and Ayabe 1999; Jung et al.2000), followed by dehydration to produce daidzein and genistein either spontaneously or via 2-hydroxyisoflavanone dehydratase (HID) (Akashi, Aoki and Ayabe 2005). Glycitein is biosynthesized from liquiritigenin by flavonoid 6-hydroxylase to form 6-hydroxyliquiritigenin (Latunde-Dada et al.2001). 6-Hydroxyliquiritigenin is then converted to 6-hydroxydaidzein by IFS and HID and subsequently converted to glycitein by isoflavone O-methyltransferase (Uchida et al.2020).

Although direct biochemical evidence is lacking, most flavonoids in plant cells are assumed to be accumulated in vacuoles in the form of glucosides. Vacuolar sequestration of flavonoids involves transporters, glutathione S-transferase, and vesicle trafficking (Zhao 2015). Both ATP-binding cassette (ABC) transporters and multidrug and toxin extrusion (MATE) transporters have been identified as being involved in the vacuolar sequestration of flavonoids (Zhao 2015). Most of the research on vacuolar sequestration has analyzed the sequestration of anthocyanins and proanthocyanidins (Zhao 2015; Shitan and Yazaki 2020); however, biochemical transport analysis using yeast membrane vesicles revealed that AtABCC2 of Arabidopsis thaliana transports flavonoid glucosides such as luteolin 7-O-glucoside and apigenin 7-O-glucoside, in addition to anthocyanins (Behrens et al.2019), and MtMATE2 of Medicago truncatula transports apigenin 7-O-glucoside, apigenin 7-O-glucoside malonate, kaempferol 7-O-glucoside, and kaempferol 7-O-glucoside malonate, in addition to anthocyanins (Zhao et al.2011).
Secretion into the rhizosphere
Root exudates and root border cells are the primary sources of rhizosphere flavonoids (Hassan and Mathesius 2012; Sasse, Martinoia and Northen 2018). ABC transporters have been suggested to mediate the secretion of flavonoids into the rhizosphere. A biochemical transport assay using plasma membrane vesicles of soybean revealed the involvement of ABC transporters in the secretion of genistein, an isoflavone acting as a signal for nod gene expression in rhizobia (Sugiyama, Shitan and Yazaki 2007). Tobacco BY-2 cells expressing MtABCG10 showed the efflux of isoliquiritigenin from cells (Biala et al.2017).
Soybean is a suitable model plant to study flavonoid secretion because of its relatively large leaves and roots and the importance of secreted flavonoids in its interactions with rhizosphere microbes (Kosslak et al.1987; Sugiyama 2019; Okutani et al.2020). Three types of isoflavone aglycones—daidzein, genistein, and glycitein—are biosynthesized in the cytosol of soybean, and their glucosides and malonylglucosides are presumably accumulated in the vacuoles (Figure 2). These aglycones are also secreted from roots and found in root exudates (Pueppke et al.1998). In addition to the ATP-dependent active transport of isoflavone aglycones, secretion of isoflavone glucosides stored in vacuoles into the apoplast has also been proposed (Suzuki et al.2006; Sugiyama 2019). Secreted isoflavone glucosides are hydrolyzed to aglycones by isoflavone conjugate hydrolyzing β-glucosidase (ICHG) (Suzuki et al.2006).
The secretion of isoflavones from soybean roots has been analyzed both in hydroponic conditions and in the field. Both daidzein and genistein induce nod genes of the rhizobia, leading to a symbiosis for biological nitrogen fixation (Kosslak et al.1987; Pueppke et al.1998). Under nitrogen deficiency, the secretion of these isoflavones is increased approximately 10-fold in hydroponic conditions (Sugiyama et al.2016) (Figure 3). Daidzein is the predominant isoflavone in root exudates throughout the growth stages; a higher amount of daidzein is secreted during the vegetative stages than the reproductive stages, and the secretion of malonyldaidzin and daidzin is increased during the reproductive stages (Sugiyama et al.2016). The amount of daidzein and genistein in root exudates is stable during the day. However, genes for transcription factors and biosynthesis of isoflavone metabolism show diurnal regulation with increased expression during the daytime (Matsuda et al.2020). GmMYB176, a transcription factor of isoflavone biosynthesis highly expressed in roots, induces isoflavone biosynthetic genes from dawn to noon, followed by the induction of isoflavone biosynthetic genes at noon, and a slight increase of daidzein aglycone in roots in the afternoon (Matsuda et al.2020). Conceivably, the secretion pathway from vacuolar glucosides to apoplast would be induced during the nighttime when the expression of ICHG is increased to maintain the daidzein level in root exudates.
In field-grown soybean, the secretion of daidzein is also higher in the early vegetative stages than in the reproductive stages, but the amount is increased up to 10 000-fold, compared with hydroponic conditions (Toyofuku et al.2021) (Figure 3). The increase of secretion is possibly due to the presence of microbes in the field; however, soil particles in the field may influence the root morphology and the rate of isoflavone secretion because the particle size and chemistry of growth substrates such as sand, clay, and glass beads, affect root morphology and exudation in Brachypodium distachyon (Sasse et al.2020). It remains unclear whether both pathways of isoflavone secretion, i.e. active transport of isoflavone aglycones and secretion of vacuole-stored isoflavone glucosides, operate simultaneously or if one operates conditionally. Although higher daidzein secretion is suggested during the early vegetative stages, the rhizosphere isoflavone contents are slightly higher in the reproductive stage due to relative stability of the compounds in the soil (Sugiyama et al.2017). It is of particular importance to analyze the microbial communities and physicochemical properties of soils together with degradation kinetics since the stability of flavonoids depends on the soil (Sugiyama and Yazaki 2014).
Roles in the rhizosphere
Flavonoids display a broad range of biological activities not only in plants but in animals and microbes. The biological activities in humans include antioxidative activity, free radical scavenging capacity, anti-inflammatory capacity, and anticancer activities (Gorniak, Bartoszewski and Kroliczewski 2019), and flavonoids are widely used as phytomedicines. Flavonoids inhibit a range of root pathogens in the rhizosphere (Hassan and Mathesius 2012; Mierziak, Kostyn and Kulma 2014) because of their antipathogenic properties, such as the ability to disrupt membrane integrity (Weinstein and Albersheim 1983; Wu et al.2019) and inhibit DNA gyrase (Wu et al.2013). These defensive flavonoids can be exuded into the soil either constitutively or inductively; for example, an increased level of glyceollin I, a phytoalexin of soybean, is found in the root exudate upon exposure to pathogens and nonsymbiotic rhizobia (Schmidt, Parniske and Werner 1992; Lozovaya et al.2004), isoflavones in the root exudates of white lupin are increased upon treatments with various elicitors (Gagnon and Ibrahim 1997). Other roles of flavonoids in the rhizosphere include mediation of allelopathy, chelation, and reduction of metals in soil (Cesco et al.2012; Hassan and Mathesius 2012; Weston and Mathesius 2013). The following sections highlight the functions of flavonoids as regulators of the rhizosphere microbiome.
Roles as chemoattractants
In addition to functioning as a nutrient source for rhizosphere microorganisms, plant metabolites (particularly PSMs) exert both attractive and repellent effects on soil microbes and shape the rhizosphere microbiota (Pascale et al.2020; Jacoby, Koprivova and Kopriva 2021). Flavonoids have been thought of as an attractant for the rhizobia, as the chemical signals from plant roots to symbiotic rhizobia have been determined to be flavonoids. Several reports have shown a chemotactic activity of nod gene-inducing flavonoids, including luteolin, 4′,7-dihydroxyflavone, 4′,7-dihydroxyflavanone, and 4,4′-dihydroxy-2-methoxychalcone for Ensifer meliloti (Caetanoanolles, Cristestes and Bauer 1988; Dharmatilake and Bauer 1992) and also apigenin and luteolin for Rhizobium leguminosarum (Aguilar et al.1988); however, a recent comprehensive analysis revealed that E. meliloti did not show a chemotactic response to flavonoids such as hyperoside, luteolin, luteolin-7-O-glucoside, quercetin, and chrysoeriol in alfalfa seed exudates (Compton et al.2020). Amino acids and quaternary ammonium compounds, which are 10-fold more abundant than flavonoids in alfalfa root exudates, are the primary chemoattractants of E. meliloti (Compton et al.2020). Bacterial chemoreceptors for these metabolites have also been identified in E. meliloti (Webb et al.2014; Webb et al.2017). In the case of soybean-rhizobia interactions, Bradyrhizobium japonicum is most attracted to succinate, glutamate, and malonate and is not attracted to luteolin, daidzein, or genistein (Barbour, Hattermann and Stacey 1991). B. japonicum is also attracted to cinnamic acid and hydroxycinnamic acids, such as p-coumaric acid, caffeic acid, ferulic acid, and sinapinic acid (Kape, Parniske and Werner 1991). Together, these studies do not favor the contribution of flavonoids in the recruitment of rhizobia to the proximity of plant roots but indicate the contribution of amino acids, dicarboxylic acids, and quaternary ammonium compounds. Future research is needed to highlight the functions of these chemoattractants in the plant rhizosphere with consideration for the secretion, degradation, and distribution within the context of the soil community where multiple interactions occur.
Roles as nod gene inducers
Once attracted to the vicinity of the root surface, chemical signal exchange between legume plants and the rhizobia occurs (Janczarek et al.2015; Liu and Murray 2016). Luteolin and 7,4′-dihydroxyflavone were the first signaling compounds discovered in alfalfa and white clover, respectively, using a nod-lacZ expression system, which was also used to identify other signaling flavonoids in other legume species, such as daidzein and genistein in soybean (Peters, Frost and Long 1986; Redmond et al.1986; Kosslak et al.1987; Kape, Parniske and Werner 1991). Different classes of flavonoids, including flavones, flavonols, flavanones, isoflavones, and chalcones, have since been identified as plant signals for induction of nod genes in the rhizobia in addition to nonflavonoid metabolites (Janczarek et al.2015; Liu and Murray 2016). Most of these flavonoids induce nod genes at low micromolar concentrations within the range of the rhizospheric concentrations, at least for daidzein (Sugiyama et al.2017; Toyofuku et al.2021). The specificity of flavonoid profiles in each legume plant and the specific perception of flavonoids by NodD, a LysR-type transcription regulator, in the rhizobia are responsible for the first level of host specificity for legume–rhizobia symbiosis. Upon perception of the signal, the rhizobia synthesize and secrete lipo-chitooligosaccharide Nod factors (NFs) to be recognized by receptors on the root surface (Gourion et al.2015; Buhian and Bensmihen 2018). Genetic signaling pathways governing NF detection and nodule organogenesis are now quite well understood in the model legume plants, M. truncatula and Lotus japonicus (Roy et al.2020). In addition to NFs, type III secretion systems, together with exported proteins, are induced by genistein in Sinorhizobium fredii USDA257, a symbiont of soybean and other legume plants (Krishnan et al.2003), and B. elkanii SEMIA587 (de Campos et al.2011). Genistein also induces the expression of resistance-nodulation-division (RND) efflux pumps in B. japonicum (Takeshima et al.2013). A B. japonicum mutant deficient in this efflux pump is sensitive to genistein (but not to daidzein) and shows reduced nodulation and nitrogen fixation when this mutant is inoculated in soybean roots. The expression of this efflux pump is negatively regulated by a TerR-like regulator (BdtR) (Han et al.2020). Mutation of this regulator results in higher extracellular genistein levels and decreased susceptibility to genistein because of induction of the efflux pump. In contrast, the induction of nod genes is reduced in the mutant. These results suggest that the rhizobia maintain intracellular genistein homeostasis to induce nod genes for nodulation while the toxic effect of this isoflavone is alleviated.
Roles in modulating the microbiota
Flavonoids in the rhizosphere impact the microbiota as well as the rhizobia. The application of pure flavonoid compounds modifies both bacterial and fungal communities in soil. Daidzein and genistein applied to the soil resulted in different microbial community structures, as revealed by phospholipid fatty acid profiling (Guo et al.2011). The treatment of soils with 7,4′-dihydroxyflavone, a nod gene-inducing flavonoid of alfalfa roots, modified bacterial communities at a concentration found in root exudates, with an increase in Acidobacteria, suggesting multiple functions of this flavonoid in the rhizosphere beyond the establishment of symbiosis (Szoboszlay, White-Monsant and Moe 2016). Both bacterial and fungal communities of the peanut rhizosphere were modified when luteolin was continuously applied to the soil (Wang et al.2018). The growth and nodule formation of peanuts treated with luteolin were reduced, suggesting the inhibitory effects of luteolin in continuous monocropped peanut systems, leading to reduced productivity. Daidzein also modifies bacterial communities in soil (Okutani et al.2020). The relative abundance of Comamonadaceae was increased in a concentration-dependent manner, and the bacterial communities became more similar to the rhizosphere of soybean grown in fields rather than bulk soil. The rhizosphere bacterial communities of soybean hairy roots silenced with the IFS gene showed a slight change, depending on the gene silencing and hairy root transformation (White et al.2017). Comamonadaceae were not reduced in the IFS-silenced hairy root rhizosphere, suggesting the involvement of multiple metabolites in modulating the rhizosphere microbiota. In maize, flavones such as apigenin and luteolin are secreted from the roots and promote the enrichment of Oxalobacteraceae in the rhizosphere (Yu et al.2021). Oxalobacteraceae isolates belonging to the genus Massilia improve the growth of maize under nitrogen-deficient conditions via alteration of root development, suggesting a network of root architecture and the microbial taxa in the rhizosphere, resulting in improved plant growth under nutrient deficiency.
Saponins
Synthesis and accumulation
Saponins are a group of PSMs widely distributed in higher plants (Vincken et al.2007). They contain an aglycone hydrophobic backbone bound to hydrophilic saccharides such as glycosides, resulting in amphiphilicity and the formation of a soap-like foam when agitated in water. The name saponin is derived from the Latin word “sapo,” meaning soap, and they are typically subdivided into triterpenoid and steroid glycosides based on the carbon skeletons (Vincken et al.2007). Both types of compounds are biosynthesized from a common precursor, 2,3-oxidosqualene, via multiple reactions such as cyclization, oxidation, and glycosylation. In the plant kingdom, dicotyledonous plants mainly accumulate triterpenoid saponins, while monocotyledonous plants mainly synthesize steroidal saponins (with some exceptions) (Sparg, Light and van Staden 2004; Moses, Papadopoulou and Osbourn 2014).
The triterpenoid and steroidal aglycone backbones are synthesized from isopentenyl diphosphate units derived from the mevalonate pathway. Condensation of 2 farnesyl diphosphate by squalene synthase (SQS) generates squalene, which is then epoxidized to 2,3-oxidosqualene catalyzed by squalene epoxidase (SQE). Oxidosqualene cyclases (OSCs) catalyze the cyclization of 2,3-oxidosqualene, and they are positioned at a key metabolic branch point between primary metabolism for plant sterols (phytosterols) and brassinosteroid hormones and specialized metabolism for triterpenoids. 2,3-Oxidosqualene can be cyclized into a diverse range of compounds with triterpene backbones, including dammaranes, tirucallanes, lupanes, hopanes, oleananes, taraxasteranes, ursanes, lanostanes, and cucurbitanes (Vincken et al.2007). OSCs catalyzing the cyclization of 2,3-oxidosqualene are either specific or multifunctional, leading to either a single product or multiple products from a single reaction (Moses, Papadopoulou and Osbourn 2014). Following cyclization, triterpene aglycones are oxidized by cytochrome P450 and further modified by an array of transferases such as UDP-dependent glycosyltransferases (UGTs), acyltransferases, malonyltransferases, and methyltransferases (Thimmappa et al.2014; Seki, Tamura and Muranaka 2015).
Soyasaponins are triterpenoid saponins commonly found in legume plants. They are composed of aglycone and oligosaccharide moieties. Soyasaponins are subdivided into 4 groups based on the aglycone structures: soyasaponin group A, group B, and group E (derived from soyasapogenol A, soyasapogenol B, and soyasapogenol E, respectively), and DDMP saponin consisting of soyasapogenol B with a DDMP (2,3-dihydro-2,5-dihydroxy-6-methyl-4H-pyran-4-one) residue at the C-22 position of the aglycone. Soyasapogenol B is biosynthesized by the hydroxylation of β-amyrin at C-22 and C-24, soyasapogenol A has an additional hydroxylation at C-21, and soyasapogenol E has a carbonyl group at C-22 (Zhang and Popovich 2009; Seki, Tamura and Muranaka 2015; Krishnamurthy et al.2019). These soyasapogenols are further diversified by glycosylation catalyzed by UGTs to form the soyasaponin groups A, B, and E. DDMP saponin is biosynthesized from group B saponin, possibly via a UGT encoded by the Sg-9 locus (Sundaramoorthy et al.2019) (Figure 4).

Secretion of soyasaponin from soybean roots and tomatine from tomato roots: cycloartenol synthase (CAS), β-amyrin synthase (BAS), UDP-sugar-dependent glycosyltransferase (UGT), arabinose (Ara), galactose (Gal), glucose (Glc), glucuronic acid (GlcA), rhamnose (Rha), and xylose (Xyl). CAS and BAS are members of OSC family. Bold arrows represent multiple reactions.
Cycloartenol synthase (CAS) catalyzes the cyclization of 2,3-oxidosqualene to cycloartenol, which is a precursor for phytosterols including cholesterol, campesterol, and β-sitosterol. Cholesterol is then oxidized and glycosylated to form steroidal saponins (Figure 4). Steroidal glycoalkaloids, typically found in species of Solanum as toxic substances, are also biosynthesized from cholesterol as a precursor but incorporate an amine group at C-26 to generate aglycones such as tomatidine and solanidine (Harrison 1990; Cardenas et al.2015), which are subsequently glycosylated at the C3 hydroxy group and accumulated as α-tomatine in tomato (Solanum lycopersicum) and α-solanine in potato (Solanum tuberosum) (Friedman 2006). Genes involved in saponin biosynthesis have been identified from various plant species (Thimmappa et al.2014; Lee et al.2019; Chung et al.2020; Jozwiak et al.2020) and are often organized as a cluster in the genome (Itkin et al.2013; Nutzmann, Scazzocchio and Osbourn 2018; Akiyama et al.2021).
Secretion into the rhizosphere
Transporters of saponins for vacuolar accumulation or secretion to the rhizosphere have not been identified so far (Francisco and Martinoia 2018). The transporter responsible for the relocation of α-tomatine from the vacuole to the cytosol has recently been identified in tomato (Kazachkova et al.2021) (Figure 3). Secretion of soyasaponins from soybean roots and tomatine from tomato roots have been investigated both in hydroponic culture and in field conditions. The secretion of soyasaponins from soybean roots was discovered during a metabolomic analysis of soybean root exudates (Tsuno et al.2018) as a strong peak on the total ion chromatogram that did not show clear UV absorption typical to aromatic rings. A detailed assessment of m/z values and MS/MS spectra, together with authentic standard specimens, identified this peak as soyasaponin Bb.
Soyasaponins are secreted during growth in hydroponic culture but peak at the early growth stage, as is the case for isoflavones of soybean (Tsuno et al.2018) (Figure 4). Soybean roots secrete both group A and group B soyasaponins, but the secretion of DDMP saponins is limited despite the predominant accumulation in soybean roots. The differential composition of soyasaponins between roots and root exudates suggests the involvement of regulatory mechanisms such as transporters and apoplastic enzymes to be elucidated in future studies. Other legume plants including, but not limited to, L. japonicus, alfalfa (Medicago sativa), and pea (Pisum sativum), secrete soyasaponins at varying concentrations and compositions, but DDMP saponins are not detectable from root exudates in these legume species (Tsuno et al.2018). The contents of the major soyasaponins, soyasaponin Ab and soyasaponin Bb, in root exudates have no apparent diurnal pattern, although the gene expression levels of β‐amyrin synthase, cytochrome P450, and UDP‐glucuronosyltransferase involved in soyasaponin biosynthesis (Seki, Tamura and Muranaka 2015; Krishnamurthy et al.2019) are higher at night (Matsuda et al.2020). The amount of soyasaponins in the rhizosphere of soybean grown in field is slightly increased during the growth stages, but the composition is stable, with group B soyasaponins representing about 60% of the total soyasaponins, followed by group A, E, and DDMP soyasaponins. Soyasaponin Bb represents up to 70% of group B soyasaponins, and soyasapogenols are limited throughout the growth stages (Fujimatsu et al.2020).
Tomatine and its aglycone tomatidine were also found in root exudates of tomato (Kirwa et al.2018). In hydroponically grown tomato, the concentrations of tomatine and tomatidine in root exudates are higher during the early growth stages than the later growth stages, consistent with the secretion of isoflavones and soyasaponins during growth (Nakayasu et al.2021) (Figure 4). The regulations of tomatine secretion remain elusive, but it is known that it is regulated systemically by the addition of glycosylated azelaic and possibly by soil microbiota (Korenblum et al.2020). In field-grown tomatoes, rhizosphere tomatine contents are comparable between the flowering and green-fruit stages. Neither tomatine nor tomatidine is detectable in bulk soil, suggesting that tomatine is secreted from field-grown tomato plants and accumulates in the rhizosphere throughout the growth stages (Nakayasu et al.2021).
Roles in the rhizosphere
Saponins exert a diverse range of biological properties relevant to rhizosphere interactions including antibacterial, antifungal, and insecticidal activities, in addition to various pharmaceutical effects in humans (Vincken et al.2007; Augustin et al.2011; Cheok, Salman and Sulaiman 2014). The ecological significance of saponins has long been of particular interest for their relevance in crop production. Biological activities against pathogens have been reported; for example, minutosides extracted from Allium minutiflorum have antimicrobial activities against various fungal and oomycete pathogens such as Alternaria alternate, Botrytis cinerea, Fusarium oxysporum, Fusarium solani, Pythium ultimum, and Rhizoctonia solani (Barile et al.2007), alliospirosides extracted from Allium cepa have antifungal activities against a range of fungi including B. cinerea and Colletotrichum gloeosporioides (Teshima et al.2013); and aescin from Aesculus hippocastanum shows antifungal activities against Microdochium nivale, Pyrenophora teres, and Leptosphaeria maculans (Trda et al.2019). Additionally, steroidal glycoalkaloids such as α-solanine, α-chaconine, and α-tomatine show hatching stimulation activity in potato cyst nematode eggs, albeit weaker than solanoeclepin A, a triterpene secreted from potato roots (Shimizu et al.2020).
Roles as allelochemicals
Saponins have also long been recognized as allelochemicals (Oleszek and Górski 1992). Saponins in the roots of alfalfa (M. sativa) and red clover (Trifolium pratense) shows inhibitory effects on seed germination and seedling growth of wheat (Triticum aestivum) (Oleszek and Jurzysta 1987; Oleszek 1993), and medicagenic acid saponins exhibit autotoxic effects and significantly lower the germination rate of alfalfa seeds as well (Ghimire et al.2019). Ginsenosides from Panax notoginseng exhibit autotoxicity and these saponins accumulated in soil are suggested to cause the replant failure of this plant (Yang et al.2015).
Roles in modulating microbiota
The function of saponins in modulating the rhizosphere microbiota has recently been identified for saponins. An oat sad1 mutant deficient in avenacin, a triterpenoid saponin, harbors altered rhizosphere microbial communities, compared with the wildtype (Turner et al.2013). We used tomato as a model system to investigate the role of saponins in shaping the rhizosphere microbiota. Tomatine and tomatidine treatments of field soil in vitro resulted in the enrichment of 20 and 2 families, respectively, and the depletion of 35 and 78 families, respectively. The bacterial communities of both tomatine- and tomatidine-treated soils are more similar to those of the tomato rhizosphere than of bulk soil, as shown by a smaller weighted UniFrac distance between PSM-treated soil and the rhizosphere than between PSM-treated soils and bulk soil. Sphingomonadaceae was the only family enriched in both PSM-treated soil and tomato rhizosphere soil. Among Sphingomonadaceae, the genus Sphingobium was particularly enriched, and this increase was attributable to 1 amplicon sequence variant (ASV) among the 15 ASVs observed in these soil samples (Figure 4). Sphingomonadaceae was also enriched in soil treated with soyasaponin Bb (Fujimatsu et al.2020); operational taxonomic units annotated as Novosphingobium were particularly enriched in soyasaponin Bb-treated soil as well as the soybean rhizosphere (Figure 4). It is of particular interest to investigate whether differences in the carbon skeletons of saponins affect the influence of saponins in rhizosphere bacterial communities. While the effect on bacterial communities is not known, ginsenosides (triterpenoid saponins in Sanqi ginseng) alter fungal communities with the enrichment of potentially pathogenic taxa including Alternaria and Fusarium, and with the depletion of potentially beneficial taxa such as Acremonium, Mucor, and Ochroconis (Li et al.2020), which is a potential cause of replanting failure in Sanqi ginseng.
Dynamics of flavonoids and saponins in the rhizosphere
The dynamics and interactions of PSMs and microbes in the rhizosphere are of particular importance for gaining insight into the functions of PSMs. Traditionally, the amount of PSMs secreted from roots has been analyzed using hydroponic cultures because of the ease of applicability compared with field sampling (Oburger and Jones 2018). Although research in our laboratory has revealed that the secretion of daidzein from roots follows a similar pattern during the growth stages (i.e. higher secretion in early growth stages than later growth stages), the amount of daidzein secreted from roots was much higher in the field (Toyofuku et al.2021), pointing toward the importance of measurement of metabolites in field-grown plants. Along with the differences in the amount and rate of secretion, it is technically challenging to analyze the positional secretion rate using a hydroponic system. In hydroponic culture, the total amount of metabolites secreted from whole roots is usually measured; in reality, the secretion rate varies based on root area (Pearson and Parkinson 1960; McDougall and Rovira 1970; Weisskopf et al.2005). Recently, analysis using matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) showed a distribution of various metabolites including PSMs in the rhizosphere (Velickovic et al.2020). In this study, a polyvinylidene fluoride (PVDF) membrane was placed against the soil-root interface of plants grown in a rhizobox. The PVDF membrane was then mounted on a MALDI target plate for mass imaging. Biosensors are also useful for revealing the distribution of metabolites in the rhizosphere. Fux fusion bioreceptors have been developed in R. leguminosarum to detect specifically for sugars, polyols, amino acids, organic acids, and flavonoids (Pini et al.2017). Although a bacterial association with the plant root system needs to be validated, these biosensors are potentially applicable for a broader range of bacteria and metabolites to investigate the spatiotemporal distribution of PSMs in the rhizosphere.
Degradation
The amount of secretion and the stability of PSMs need to be analyzed to understand the dynamics of PSMs in the rhizosphere. Once secreted, PSMs are degraded by soil microbes. Despite the rather toxic properties of flavonoids and saponins, these PSMs are also a carbon source for microbes possessing catabolic enzymes. The stability of PSMs in the rhizosphere varies depending on the metabolites and soil microbial communities. In contrast to primary metabolites that are degraded within a few minutes (Gunina and Kuzyakov 2015), flavonoids are rather stable in the soil. The majority of naringenin and formononetin is degraded within 96 h (Shaw and Hooker 2008), while the half-life of apigenin, kaempferol, and their derivatives ranges from 4.4 to 14.7 days with a seasonal variation (Sosa et al.2010). The half-life of daidzein is about 7 days, measured in soil collected from a soybean farm (Okutani et al.2020) (Figure 3). The half-life of saponins in the soil is comparable with flavonoids. The allelopathic effects of saponins from alfalfa and red clover last for 1-2 weeks depending on the soil (Oleszek and Jurzysta 1987), and most soyasaponin Bb and tomatine in soil are degraded within a week (Fujimatsu et al.2020; Nakayasu et al.2021).
Aerobic flavonoid biodegradation has been reported for a wide range of bacterial species (Shaw, Morris and Hooker 2006). Nod gene-inducing daidzein and genistein are also to be degraded by rhizobial species via multiple C-ring fission, although genes involved in this pathway have not been identified (Rao and Cooper 1994, 1995). The metabolic pathways for flavonoid degradation have been characterized in intestinal microbes; daidzein, when consumed as a part of soy products, is converted to a reduction product, equol, and a C-ring cleavage product, O-demethylangolensin, by anaerobic intestinal bacteria (Feng et al.2018). These products have not been identified in rhizosphere soil to our knowledge. Screening of mutants defective in naringenin catabolism in Herbaspirillum seropedicae revealed the involvement of a monooxygenase, FdeE, in the cleavage of naringenin, possibly together with FdeD, a putative Rieske protein in fed operon (Marin et al.2013; Marin et al.2016). Other proteins encoded by the fde operon include FdeH, a cupin family protein that contains quercetin 2,3-dioxygenase for cleaving quercetin into 2-protocatechuoyl-phloroglucinol carboxylic acid and carbon monoxide, identified in Bacillus subtilis (Bowater et al.2004; Barney et al.2004) and Streptomyces sp. (Merkens et al.2007). Most of the genes involved in flavonoid catabolism, especially those for isoflavones, have not been identified.
Intestinal and fecal microorganisms capable of metabolizing saponins have also been reported (Hu et al.2004; Dong et al.2017; del Hierro et al.2018). Sapogenins such as oleanolic acid, hederagenin, serjanic acid, diosgenin, and soyasapogenol B accumulate when plant extracts rich in saponins (e.g. quinoa, lentil, and fenugreek) are fermented by gut microbiota (del Hierro et al.2020). Sapogenins are absorbed by the gastrointestinal tract, and no gut microorganisms capable of degrading sapogenin have been identified. In contrast, it is not well understood what kind of microorganisms decompose compounds in the soil. Microorganisms capable of degrading tomatine into tomatidine have been reported (Ford, McCance and Drysdale 1977; Okmen et al.2013), and Sphingobium spp. isolated from tomatine-treated soil degrade tomatine and tomatidine and use them as a carbon source (Nakayasu et al.2021). Saponins such as ginsenosides and soyasaponins, are also metabolized by intestinal and fecal microbes (Hu et al.2004; Dong et al.2017). Yet, the pathways and genes involved in the degradation of aglycones have not been characterized, particularly for rhizosphere microbes.
Adsorption
In addition to the degradation by soil microbes, adsorption by organic matter and clay minerals reduces the distribution of metabolites in the rhizosphere. We used daidzein as a model to investigate the adsorption to gray lowland soil collected from a soybean farm. Possible adsorption sites are humic substances and clay minerals. Humic substances contain hydroxyl and phenolic hydroxyl groups involved in the formation of complexes with organic substances (Pei Gan and Yau Li 2013). Decomposition of organic matter in gray lowland soil reduced daidzein adsorption, suggesting the involvement of humic substances in limiting the distribution of daidzein in the rhizosphere. In contrast, adsorption of daidzein to clay minerals such as kaolinite, a 1:1-type silicate mineral, is much lower than that to gray lowland soil, and adsorption to illite, a 2:1-type silicate material, is undetectable (Okutani et al.2020) (Figure 3).
Simulation
The distribution of mineral ions and water in the soil surrounding plant roots has been simulated (Duncan et al.2018; Zarebanadkouki et al.2014; Vereecken et al.2016). For PSMs in rhizosphere soil, the advection-diffusion (dispersion) equation has recently been used to simulate the dynamics of daidzein in the rhizosphere (Okutani et al.2020). A single root of diameter 2 mm and length 10 cm was set in the center of soil with a diameter of 20 cm and a depth of 20 cm. The root length and diameter were assumed to be constant, and daidzein secretion was assumed to be equal from all parts of the roots. Daidzein distribution was predicted to within a few millimeters of the root surface during the early growth stages (14 days). Although daidzein distribution was shown to be within 2 mm in a rhizobox experiment (Okutani et al.2020), this simulation is based on constant soil environmental conditions, the absence of root growth, and equal secretion from all parts of the roots. It is necessary to incorporate, at a minimum, the soil water contents (i.e. wet conditions on rainy days and dry conditions on sunny days) and root growth with differing secretion rates of metabolites to precisely simulate the distribution of metabolites in rhizosphere soil in field-grown conditions. Also, simulation based on the secreted amount of daidzein in hydroponic culture underestimates the amount of rhizosphere daidzein. When the daidzein secretion rate measured in field-grown soybean was used for the simulation, the estimated daidzein concentration in the soybean rhizosphere was within the range of that from the rhizosphere in the soybean field (Toyofuku et al.2021). So far, the adsorption coefficient has been analyzed for flavonoids and saponins in addition to organic xenobiotics and estrogens in soil (Shaw and Hooker 2008; Caron et al.2010; Fujimatsu et al.2020); thus, simulations based on a fluid model can be applied to predict the rhizosphere distribution of PSMs.
Conclusion and future perspectives
Research in the past few decades has identified a vast array of roles of flavonoids and saponins in the rhizosphere, especially for mediating interactions between plants and microbes. Recent advancement of multiomics analysis revealed a tight network between host plants and microbiota, which opens up a holistic approach toward a comprehensive understanding of plants and microbes. This concept considers the host plant and its microbiota as a unique biological entity called “holobiont,” in which the host and microbiota interact to affect morphology, development, and physiology among others (Rosenberg and Zilber-Rosenberg 2016; Hassani, Duran and Hacquard 2018). It remains largely unclear how the metabolic network in the holobiont is established and how it affects plant growth, fitness, and robustness to changing environments. For crop species, the domestication process has affected root microbiota, mediated at least partially by the alteration of root exudates (Iannucci et al.2017; Escudero-Martinez and Bulgarelli 2019). The metabolic network in the holobiont provides a valuable basis for designing an optimized microbiota to confer robustness against both biotic and abiotic stresses and, eventually, to improve crop yields. Flavonoids and saponins are key metabolites for such an approach since these bioactive compounds mediate the interaction between the plant and the microbiota, and recent evidence has revealed the possibility of fine-tuning the interactions (Fujimatsu et al.2020; Nakayasu et al.2021).
The direct application of metabolites is the first step toward this goal. Flavonoids applied to seeds improve the nodulation in several legume crops (Mabood, Zhou and Smith 2014), and saponins applied to seeds improve salinity stress tolerance in quinoa and soybean (Yang et al.2018; Soliman et al.2020). There remain possibilities to utilize these PSMs in agriculture directly, but an obstacle is the stability of these metabolites in the soil. Multiple applications are necessary to exert bioactive properties during a crop season, and the application of PSMs to the rhizosphere is technically unfeasible. Designing a holobiont–metabolic network utilizing both plant breeding and microbial inoculation would be promising to circumvent these problems. Designing rhizosphere microbiomes for crop productions has been proposed in reviews (Mueller and Sachs 2015; de Souza, Armanhi and Arruda 2020; Pascale et al.2020), and the holobiont–metabolic network could be a key component.
Metagenomic analysis revealed key bacterial taxa in the microbiome for beneficial traits. By comparing the rhizosphere metagenomes of resistant and susceptible tomato varieties to pathogenic Ralstonia solanacearum, it was found that Flavobacterium was abundant in the rhizosphere of resistant varieties, and isolated Flavobacterium suppressed R. solanacearum when inoculated in pots (Kwak et al.2018). Metagenomics on disease suppressive soil identified a consortium of Chitinophaga and Flavobacterium for the suppression of the fungal root pathogen R. solani, and gene clusters encoding the production of nonribosomal peptide synthetases and polyketide synthases in Flavobacterium were found to be essential for disease suppression (Carrion et al.2019). In contrast to metagenomics, the identities and functions of metabolites in the rhizosphere largely remained unknown because of the instability and limited amounts of metabolites in soil and difficulties with extraction from the rhizosphere. In addition to PSMs, microbial metabolites exert diverse ranges of influence on host plants and co-occurring microbes (Backer et al.2018; Weisskopf et al.2019). The rhizosphere metabolome is a prominent approach for uncovering the novel metabolites in the rhizosphere. We recently found okaramine A, B, and C in the rhizosphere of hairy vetch (Vicia villosa), a cover crop or a green manure crop (Sakurai et al.2020). Okaramines were first identified as insecticides from okara inoculated with P. simplicissimum AK-40 (Hayashi et al.1989), but they had not been identified in nature. Okaramine B was also detected in the rhizosphere of soybean grown after hairy vetch but not in soybean without previous hairy vetch cultivation, suggesting an interspecies soilborne legacy or an indirect defense of plants against pests (Sakurai et al.2020; Matsuda 2018).
Multiomics analyses integrating metagenomics and rhizosphere metabolomics are prerequisites for designing microbiomes based on holobiont–metabolic networks. The rhizosphere is vital for plant growth and crop production. Unlike the accumulation of microbial genomics, metabolomics has only revealed the tip of the iceberg. It will be of particular importance to design a holobiont–metabolic network by combining the identification of new key metabolites, metagenomic analysis, isolation, and functional characterization, and plant breeding to enable plants to grow more robustly (Figure 5).

Integration of multiomics to design holobiont–metabolic networks that promote plant growth and suppress pathogens.
Acknowledgments
All research in our laboratory was performed with graduate students, postdocs, technicians, and collaborators. I would like to thank all members who participated in this research. I would like to express my sincere gratitude to Dr. Kazufumi Yazaki (professor at Kyoto University), Dr. Jiro Sekiya (emeritus professor at Kyoto University), and Dr. Nobukazu Shitan (professor at Kobe Pharmaceutical University) for continuous guidance and encouragement. I would like to thank Dr. Masaru Nakayasu for the critical reading of the manuscript, and Science Graphics Co., Ltd. and Keiko Kanai for preparing the illustration.
Data availability
The data generated during this study are available from the corresponding author upon reasonable request.
Funding
This work was supported in part by JST CREST grant JPMJCR17O2, JSPS KAKENHI grants 18H02313 and 21H02329, and funds from the Research Institute for Sustainable Humanosphere, Kyoto University.
Disclosure statement
No potential conflict of interest was reported by the author.
References
Author notes
This review was written in response to the author's receipt of the JSBBA Award for Young Scientists in 2020