Skip to main content
  • Research article
  • Open access
  • Published:

Physiological and transcriptomic analyses reveal the roles of secondary metabolism in the adaptive responses of Stylosanthes to manganese toxicity

Abstract

Background

As a heavy metal, manganese (Mn) can be toxic to plants. Stylo (Stylosanthes) is an important tropical legume that exhibits tolerance to high levels of Mn. However, little is known about the adaptive responses of stylo to Mn toxicity. Thus, this study integrated both physiological and transcriptomic analyses of stylo subjected to Mn toxicity.

Results

Results showed that excess Mn treatments increased malondialdehyde (MDA) levels in leaves of stylo, resulting in the reduction of leaf chlorophyll concentrations and plant dry weight. In contrast, the activities of enzymes, such as peroxidase (POD), phenylalanine ammonia-lyase (PAL) and polyphenol oxidase (PPO), were significantly increased in stylo leaves upon treatment with increasing Mn levels, particularly Mn levels greater than 400 μM. Transcriptome analysis revealed 2471 up-regulated and 1623 down-regulated genes in stylo leaves subjected to Mn toxicity. Among them, a set of excess Mn up-regulated genes, such as genes encoding PAL, cinnamyl-alcohol dehydrogenases (CADs), chalcone isomerase (CHI), chalcone synthase (CHS) and flavonol synthase (FLS), were enriched in secondary metabolic processes based on gene ontology (GO) analysis. Numerous genes associated with transcription factors (TFs), such as genes belonging to the C2H2 zinc finger transcription factor, WRKY and MYB families, were also regulated by Mn in stylo leaves. Furthermore, the C2H2 and MYB transcription factors were predicted to be involved in the transcriptional regulation of genes that participate in secondary metabolism in stylo during Mn exposure. Interestingly, the activation of secondary metabolism-related genes probably resulted in increased levels of secondary metabolites, including total phenols, flavonoids, tannins and anthocyanidins.

Conclusions

Taken together, this study reveals the roles of secondary metabolism in the adaptive responses of stylo to Mn toxicity, which is probably regulated by specific transcription factors.

Background

Of the mineral nutrients, manganese (Mn) is essential for plant growth and participates in a series of metabolic processes, such as photosynthesis, respiration, secondary metabolism and protein biosynthesis [1]. Mn acts as a cofactor of many enzymes, such as superoxide dismutase (SOD) and enzymes involved in the tricarboxylic acid cycle. Mn also plays roles in flavonoid and lignin biosynthesis [2]. As a trace element, Mn is only required in small amounts of 20–40 mg Mn per kilogram dry weight for most plants [3]. However, Mn is also considered to be a heavy metal that can cause phytotoxicity when it reaches the level of 150 mg per kilogram dry weight in plants [1].

In soils, available Mn levels fluctuate from 450 to 4000 mg per kilogram, and Mn solubility is mainly dependent on pH values and redox conditions [4, 5]. Hence, excess Mn toxicity generally occurs in acid soils due to the accumulation of bioactive divalent Mn (II) [5, 6]. Consequently, soil amelioration, such as lime application, is typically conducted to alleviate Mn toxicity by decreasing Mn availability, but this application is costly from both economic and environmental aspects [7]. For these reasons, breeding crop varieties with superior Mn tolerance represents a sustainable alternative agronomical strategy, but this strategy requires better understanding of how plants respond to Mn toxicity. Although morphological changes in plants grown under the condition of Mn toxicity vary among plant species, the appearances of Mn toxicity reported in most plants generally include leaf chlorosis, brown spots, crinkled leaves and brown roots, and ultimately plant growth inhibition [8, 9]. Adverse impacts caused by Mn toxicity have also been documented in plant cells at physiological levels, such as triggering oxidative stress, causing lipid peroxidation, inhibiting enzyme activity, impairing chlorophyll biosynthesis and photosynthesis and disturbing the uptake and translocation of other mineral elements [1, 9].

To counteract Mn toxicity, plants are equipped with sophisticated adaptive strategies to detoxify Mn, such as modified Mn translocation and distribution, sequestration of Mn into subcellular compartment, antioxidant system activation, and adjusted root organic acid exudation [9]. Cumulative results show that reactive oxygen species (ROS) scavenging systems involving antioxidant enzymes, including peroxidase (POD) and ascorbate peroxidase (APX), are regulated by the plant’s response to Mn toxicity, thereby alleviating excess Mn-induced oxidative stress [10, 11]. Furthermore, secondary metabolic processes and metabolites, such as phenolics, flavonoids and phenylalanine, are regulated by Mn stress in plants [2, 12, 13], suggesting the potential roles of secondary metabolism in the adaptation of plants to Mn toxicity. To date, some key genes have been characterized that participate in Mn uptake, translocation and distribution and help plants address environmental Mn stress [14, 15]. For example, as one of the natural resistance-associated macrophage protein (Nramp) members, OsNramp3 in rice (Oryza sativa) is a plasma membrane Mn transporter and responsible for Mn distribution from young leaves and panicles to old tissues, thereby protecting plants from Mn toxicity [16]. Metal tolerance proteins (MTPs) belonging to the cation diffusion facilitator (CDF) family, such as ShMTP1 from the Caribbean stylo (Stylosanthes hamata) [17], OsMTP8.1 and OsMTP8.2 from rice [18, 19], AtMTP8 from Arabidopsis [20], CsMTP8 from cucumber (Cucumis sativus) [21] and CasMTP8 from the tea plant (Camellia sinensis) [22], are involved in sequestering Mn into the vacuole for detoxification. Although the roles of the above genes have been functionally characterized, the transcriptome profiles of Mn-responsive genes in plants have not been fully elucidated. Studying the responses of plants to varying Mn concentrations is useful to determine how plants cope with Mn toxicity.

Stylo (Stylosanthes spp.) is a dominant tropical legume that is widely grown in tropical areas worldwide [23, 24]. Superior Mn tolerance is observed in the stylo compared to other legumes [25]. Recently, it has been documented that high Mn adaptability in stylo may be achieved by its fine regulation of proteins involved in specific pathways, such as defense response, photosynthesis and metabolism [11]. Furthermore, important roles of organic acids, such as malate, in stylo adaptation to Mn toxicity have been reported; a malate dehydrogenase (SgMDH1) enzyme that is up-regulated in response to Mn catalyzes malate synthesis and contributes to Mn detoxification [26]. Although stylo has considerable potential for Mn tolerance, the effects of excess Mn toxicity on profile alterations in the gene expression of stylo have not been reported, and the molecular responses of stylo to Mn stress remain largely unknown. Previous studies have paved the way for the current study dissecting the molecular responses of stylo to Mn toxicity. Accordingly, in this study, the effects of various Mn concentrations on the physiological changes in stylo were first investigated. Transcriptomic analysis of Mn-responsive genes in stylo leaves was further performed using an RNA-seq approach. The results of this study provide a platform to understand the adaptive responses of stylo to Mn toxicity and the genes involved.

Results

Effects of excess Mn stress on stylo growth

Thirty-day-old stylo plants were subjected to Mn treatments ranging from 5 to 800 μM MnSO4 for 10 d. Leaf chlorosis, a symptom of Mn toxicity, was observed in stylo leaves treated with greater than 200 μM Mn, especially 400 and 800 μM Mn (Fig. 1a). O2− levels were observed in leaves based on nitroblue tetrazolium chloride (NBT) staining (Fig. 1b). O2− accumulation was mainly observed in leaves treated with excess Mn compared with the control (5 μM), and the most intense blue color was observed in leaves treated with 800 μM Mn (Fig. 1b). Consistent with this finding, increases in malondialdehyde (MDA) levels were observed in stylo leaves subjected to high Mn stress. MDA concentrations in leaves exposed to 400 and 800 μM Mn treatments were 55.4 and 103.9% greater, respectively, compared with control conditions (Fig. 1c). Furthermore, the relative electrolyte leakage of stylo leaves was significantly increased under Mn treatments exceeding 200 μM (Additional file 1: Fig. S1). In contrast, leaf chlorophyll concentrations were decreased by 35.0–79.1% upon treatment with 200 to 800 μM Mn compared with the control (Fig. 1d). Similarly, the maximum quantum yield of photosystem II (Fv/Fm) significantly declined in stylo under Mn treatments from 200 to 800 μM Mn compared with the controls (Additional file 1: Fig. S1), suggesting that photosynthesis was inhibited by Mn toxicity.

Fig. 1
figure 1

Effects of different Mn treatments on stylo growth. a Stylo leaves with different Mn treatments. b NBT staining of stylo leaves treated with different Mn concentrations. c MDA concentrations in leaves. d Chlorophyll concentrations in leaves. Thirty-day-old stylo plants were treated with 5, 100, 200, 400 and 800 μM MnSO4 for 10 d. Values are the mean of three replicates with standard error bars. Different letters represent significant differences at P < 0.05. Bar = 1 cm

Stylo shoot and root growth were inhibited by Mn concentrations greater than 200 μM. The shoot dry weight was reduced by 29.4–50.0% with 200 to 800 μM Mn treatments, whereas root dry weight decreased by 18.3–40.2% with 200 to 800 μM Mn treatment compared with their respective controls (Fig. 2a, b). Additionally, stylo plant height decreased at Mn levels greater than 400 μM compared with the controls (Additional file 1: Fig. S1). Increases in Mn concentrations were found in both stylo shoots and roots under Mn stress. Mn concentrations in shoots and roots were increased by more than 2.2-fold and 3.5-fold with more than 100 μM Mn compared with their respective controls (Fig. 2c, d).

Fig. 2
figure 2

Plant dry weight and Mn concentrations in stylo under different Mn treatments. a Shoot dry weight. b Root dry weight. c Shoot Mn concentrations. d Root Mn concentrations. Thirty-day-old stylo plants were treated with 5 to 800 μM MnSO4 for 10 d. Values are the mean of three replicates with standard error bars. Different letters represent significant differences at P < 0.05

Enzyme activity response to excess Mn

Activities of peroxidase (POD), ascorbate peroxidase (APX), polyphenol oxidase (PPO) and phenylalanine ammonia-lyase (PAL) were analyzed in stylo leaves exposed to Mn treatments. The results showed that the tested enzymes were differentially regulated by Mn (Fig. 3). Compared to their respective controls, POD and APX activities increased as Mn treatments increased from 200 to 800 μM, peaking at 400 μM Mn (Fig. 3a, b). POD and APX activities increased by 2.9-fold and 0.41-fold with 400 μM Mn treatment compared to their respective controls (Fig. 3a, b).

Fig. 3
figure 3

Determination of enzyme activities. a POD activity. b APX activity. c PPO activity. d PAL activity. Thirty-day-old stylo plants were treated with 5 to 800 μM MnSO4 for 10 d. Values are the mean of three replicates with standard error bars. Different letters represent significant differences at P < 0.05

In addition, increases in PPO and PAL activities were also found in stylo leaves under Mn stress, especially 400 and 800 μM Mn treatments (Fig. 3c, d). PPO and PAL activities in stylo leaves treated with 400 μM Mn were 54.7 and 35.8% increased compared with their controls (Fig. 3c, d).

Transcriptome analysis of stylo leaves responding to Mn toxicity

In this study, comparative transcriptomic analysis in stylo leaves subjected to 5 and 400 μM Mn treatments was performed. Approximately 48.4 and 51.6 million clean reads were obtained from the libraries of the stylo leaves treated under control Mn (5 μM) and toxic Mn (400 μM) conditions (Additional file 2: Table S1), respectively. As the genome sequence for stylo is not available, a de novo assembly approach was employed. De novo assembly of the reads produced 234,557 transcripts corresponding to 102,872 unigenes from all samples (Additional file 2: Table S1). The mean lengths of transcripts and unigenes were 1132 and 869 bp, respectively (Additional file 2: Table S1).

Differentially expressed genes (DEGs) in stylo leaves exposed to Mn treatments were identified based on |log2(fold change)| ≥ 1 and q < 0.05. A total of 4094 DEGs were identified via comparison of the two Mn treatments (Additional file 3: Table S2). Among these genes, 2471 genes were up-regulated, while 1623 genes were down-regulated by Mn toxicity (Additional file 3: Table S2). Gene ontology (GO) analysis showed that the identified DEGs can be classified into 21 biological processes (BP), 17 cellular components (CC), and 11 molecular function (MF) terms (Fig. 4). The main categories in BP included cellular process, metabolic process, biological regulation and response to stimulus terms. The dominant categories in CC included cell part, organelle and membrane. Prominent MF categories included catalytic activity, transcription regulator activity and structural molecule activity (Fig. 4). Furthermore, DEGs were mainly involved in the following pathways based on Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis: ribosome, plant hormone signal transduction, plant-pathogen interaction, glutathione metabolism, phenylpropanoid biosynthesis, DNA replication, sesquiterpenoid and triterpenoid biosynthesis, and isoflavonoid biosynthesis (Fig. 5).

Fig. 4
figure 4

GO analysis of DEGs in stylo. The up or down-regulated genes were classified into biological process, cellular component and molecular function. X- and Y-axis indicate GO terms and the number of DEGs, respectively

Fig. 5
figure 5

KEGG enrichment analysis of the DEGs. R package ggplot2 was used to data visualization. Rich ratio is the number of significant genes divided by background genes of corresponding pathway term. The size of the dot represents the number of DEGs and the color indicates q-value (the corrected p-value) of the pathway term

Secondary metabolic pathways in response to Mn toxicity

According to GO analysis, a total of 1201 unigenes were predicted to be involved in metabolic processes in the BP categories, including 824 up-regulated unigenes and 377 down-regulated unigenes (Fig. 4). Among these DEGs, a set of 94 genes were enriched in seven secondary metabolic pathways (Additional file 4: Table S3). The top Mn-responsive genes belonged to phenylpropanoid metabolic process, flavonoid metabolic process, isoflavonoid metabolic process and lignin metabolic process (Fig. 6). Interestingly, a large number of genes in the above pathways were enhanced by Mn toxicity. For example, genes encoding PAL, β-glucosidases (GLUs), cinnamoyl-CoA reductase (CCR), cinnamyl-alcohol dehydrogenases (CADs), feruloyl-CoA 6-hydroxylases (F6Hs) and caffeic acid 3-O-methyltransferases (COMTs), were up-regulated in the phenylpropanoid biosynthesis pathway (Fig. 7). In addition, the transcript of genes encoding chalcone isomerases (CHIs), one homolog of chalcone synthase (CHS) and flavonol synthases (FLSs) associated with the flavonoid biosynthesis pathway, isoflavone 7-O-methyltransferase (I7OMT), 2-hydroxyisoflavanone dehydratase (HIDH), isoflavone 2-hydroxylase (I2H) and vestitone reductases (VRs) related to the isoflavonoid biosynthesis process was enhanced in stylo under Mn toxicity (Fig. 7). Furthermore, transcripts of PAL and a set of homologs of I7OMT and gene encoding anthranilate N-methyltransferase-like in the above pathways were increased by more than 8-fold under Mn stress (Additional file 4: Table S3).

Fig. 6
figure 6

DEGs related to secondary metabolism in stylo’s response to Mn toxicity

Fig. 7
figure 7

DEGs associated with phenylpropanoid, flavonoid and isoflavonoid biosynthesis pathways. DEGs were mapped to the reference pathways in KEGG. Gene transcripts are presented as a heatmap for log2(fold change) of gene transcription between control (5 μM) and Mn toxicity (400 μM) treatments. Gene expression levels range from red (up-regulated) to green (down-regulated). The copyright permission to use and modify these secondary metabolism pathways in the figure has been granted from KEGG

Transcription factors involved in stylo responses to Mn toxicity

A total of 123 DEGs were enriched in transcription factors (TFs) (Additional file 5: Tables S4). Among them, the largest group of TFs belonged to the AP2 family with 23 up-regulated and 1 down-regulated genes. Other DEGs encoding TFs included 22 MYBs, 16 ZFs, 15 HLHs, 11 NAMs, 9 WRKYs, 5 HSFs, 5 TCPs, 3 GRASs, 2 EINs, and one of each of the following genes: NF-X1, K-box, QLQ, PHD, SRF-TF, Homeobox, E2F-TDP, CCT, B3, HB and SBP (Fig. 8). Among them, 11 genes were significantly increased by greater than 8-fold under Mn toxicity: DREB group protein, protein PPLZ02, dehydration-responsive element-binding protein 1E-like, ethylene-responsive transcription factor 4 and TINY transcription factor belonging to the AP2 family; myb-related protein Zm1 belonging to the MYB family; zinc finger of C2H2 type belonging to the ZF family; transcription factor bHLH18-like belonging to the HLH family; hypothetical protein LR48 belonging to the WRKY family; heat shock factor protein HSF30-like belonging to the HSF family; hypothetical protein GLYMA belonging to the NF-X1 family, and MADS-box transcription factor 1-like isoform X1 belonging to the K-box family (Additional file 5: Table S4).

Fig. 8
figure 8

DEGs associated with transcription factors in stylo leaves. X- and Y-axis indicate the abbreviation and the number of TFs, respectively

Protein-protein interaction network analysis

To explore the candidate TFs involved in the transcriptional regulation of genes associated with secondary metabolism in stylo during Mn exposure, the protein-protein interaction networks of DEGs related to secondary metabolic pathways and TFs were further constructed. The interaction networks contained 305 edges with 72 nodes (Fig. 9). Among them, two genes CAD1/2 (TRINITY_DN21432_c0_g1 and TRINITY_DN45912_c0_g3) were regulated by six TFs, all of which belonged to C2H2 zinc finger transcription factors (TRINITY_DN12526_c0_g1, TRINITY_DN42276_c1_g1, RINITY_DN54332_c0_g1, TRINITY_DN13820_c0_g1, TRINITY_DN46140_c3_g3 and TRINITY_DN47158_c1_g1) (Fig. 9). Furthermore, HIDH (TRINITY_DN56723_c0_g1) can be regulated by the MYB gene (TRINITY_DN42931_c2_g2) and a NAC-like transcription factor (TRINITY_DN45594_c4_g1) (Fig. 9). In addition, genes encoding anthocyanidin reductase (TRINITY_DN52696_c0_g1) and hypothetical protein (TRINITY_DN45901_c1_g1) were the targets of the MYB gene (TRINITY_DN42931_c2_g2) (Fig. 9). These results suggest that the candidate TFs highlighted above may be involved in the regulation of secondary metabolism-related genes in the response of stylo to Mn toxicity.

Fig. 9
figure 9

Interaction networks of DEGs related to secondary metabolism pathways and transcription factors. Gene IDs are available in Additional file 4: Table S3 and Additional file 5: Table S4. The interaction networks were constructed using Cytoscape (version 3.8.0). Nodes and edges indicate genes and interactions, respectively

Validation of RNA-seq data using qRT-PCR

To confirm the RNA-seq data, 15 DEGs involved in secondary metabolism and TFs were further selected for qRT-PCR analysis (Fig. 10). These DEGs included 10 up-regulated and five down-regulated genes. The results showed a significant correlation (R2 = 0.887, p < 0.01) between the RNA-seq data and qRT-PCR results (Fig. 10), suggesting the reliability of the RNA-seq data.

Fig. 10
figure 10

Correlation analysis of gene expression between transcriptome data and qRT-PCR results. Ten up-regulated and five down-regulated DEGs were selected for qRT-PCR analysis. Transcriptome data were plotted against data from qRT-PCR. Data are presented on a log2 scale

Changes in secondary metabolite levels under Mn stress

Subsequently, the concentrations of four secondary metabolites, including total phenols, flavonoids, tannins and anthocyanidins, in stylo leaves treated with various Mn concentrations were further detected. Compared to their respective controls, the levels of total phenols, tannins and anthocyanidins in stylo leaves increased with 400 and 800 μM Mn treatments, while the concentrations of flavonoids increased with Mn treatments increasing from 200 to 800 μM, peaking at 800 μM (Fig. 11). Furthermore, the concentrations of total phenols, flavonoids, tannins and anthocyanidins in stylo leaves subjected to 400 μM Mn treatment were increased 1.1-, 1.6-, 2.1- and 7.4-fold, respectively, compared with their controls (Fig. 11). These results suggest that secondary metabolism may play a key role in stylo’s response to Mn toxicity.

Fig. 11
figure 11

Concentrations of secondary metabolites in stylo leaves subjected to different Mn treatments. a Total phenols. b Flavonoids. c Tannins. d Anthocyanins. Thirty-day-old stylo plants were treated with 5 to 800 μM MnSO4 for 10 d. Values are the mean of three replicates with standard error bars. Different letters represent significant differences at P < 0.05

Discussion

Although Mn is a trace element required for plant growth, Mn phytotoxicity generally occurs, especially in acid soil regions worldwide [8, 9]. Mn toxicity not only causes adverse impacts in plant cells but also threatens human health [27, 28]. In recent decades, genome-wide identification of metal responsive genes through transcriptomic techniques has been conducted to investigate the molecular mechanisms of plant adaptation to heavy metal stresses [29,30,31]. However, the transcriptome profiles of plants in response to Mn toxicity have not been well documented. In this study, the molecular responses of stylo to Mn toxicity were investigated through integration of physiological and transcriptomic analyses, providing potential Mn-tolerant strategies in plants for further study.

Response of stylo to Mn toxicity

Stylo growth was inhibited by high Mn stress, as reflected by the inhibition of leaf chlorophyll concentrations and plant dry weight, especially under 400 and 800 μM Mn treatments (Figs. 1 and 2). Similarly, decreases in plant growth are also observed in ryegrass (Lolium perenne), white clover (Trifolium repens) and polish wheat (Triticum polonicum) under Mn toxicity [32, 33]. Stylo growth inhibition might be due to the overproduction of ROS, as indicated by the accumulation of O2− and MDA (Fig. 1), which was supported by the results obtained from a previous study in stylo [11]. A considerable increase in oxidative stress is also found in polish wheat leaves under Mn toxicity, inhibiting shoot and root growth [33].

RNA-seq analysis was further conducted to identify DEGs in stylo’s response to Mn toxicity. The identified DEGs can be classified into numerous pathways (Fig. 4), suggesting various adjustments in stylo’s response to Mn toxicity. Among them, a large number of excess Mn-regulated genes were categorized as defense response genes, including homologous genes encoding antioxidant enzymes, such as POD and APX (Additional file 6: Table S5). This finding is consistent with a previous study demonstrating that homologous genes encoding antioxidant enzymes are up-regulated in stylo and citrus (Citrus grandis) in response to Mn toxicity [11, 34]. Furthermore, excess Mn up-regulation of POD and APX transcripts may result in increasing POD and APX activities in stylo leaves subjected to excess Mn toxicity (Fig. 3 and Additional file 6: Table S5), suggesting that these antioxidant enzymes may participate in ROS scavenging in stylo. In addition, 8 out of 10 genes encoding homologs to glutathione S-transferase (GST) were up-regulated by Mn toxicity (Additional file 6: Table S5). As a cellular detoxification enzyme, GST protects plants from oxidative damage [35, 36]. Homologous genes of GST play protective roles in plants against different stresses. For instance, the expression of OsGSTLs is enhanced by arsenic (As) stress in rice, and overexpression of OsGSTL2 increases tolerance to As, cadmium (Cd) and chromium (Cr) stresses in transgenic Arabidopsis [35].

Moreover, other DEGs belonging to defense response related proteins, such as pathogenesis-related proteins (PRs), disease-resistance proteins (DRPs), heat shock proteins (HSPs), chitinases, beta-1,3-glucanases (GLUs) and late embryogenesis abundant proteins (LEAs), were also identified in this study (Additional file 6: Table S5). Homologues of these genes are regulated by Mn or other metal stresses [11, 29, 30, 37, 38], suggesting that overlap between metal stress responses may be conserved in stylo and other plants.

DEGs involved in secondary metabolism

A group of DEGs involved in plant secondary metabolism, such as phenylpropanoid, flavonoid and isoflavonoid biosynthesis processes, were identified in this study (Figs. 6 and 7), suggesting that secondary metabolism is regulated by Mn toxicity. In the phenylpropanoid pathway, for example, the expression of PAL was up-regulated by Mn stress in the transcriptome data (Fig. 7). PAL is a key gene that catalyzes phenylalanine to form cinnamic acid in the phenylpropanoid pathway, which plays important roles in plant stress tolerance [39, 40]. Heavy metal treatments, such as lead (Pb), increase the transcript level of PAL or its enzyme activity in soybean (Glycine max), lupine (Lupinus luteus) and tomato (Lycopersicon esculentum) [41, 42]. Accordingly, increased PAL transcript levels may enhance its enzyme activity in stylo’s response to Mn toxicity (Figs. 3 and 7). Furthermore, some genes encoding proteins involved in different steps in the phenylpropanoid biosynthesis pathway, such as CCR, CADs, COMTs, F6Hs and PODs, were all up-regulated in stylo subjected to Mn toxicity (Fig. 7). Similar results are found in proteomic analysis of stylo, in which a group of proteins associated with the phenylpropanoid metabolic pathway are up-regulated under excess Mn stress [11]. Therefore, these results strongly suggest that the phenylpropanoid biosynthesis process is activated in stylo to counteract Mn stress.

The flavonoid biosynthesis process is performed by a set of key genes encoding CHSs, CHIs, FLSs and F3,5H, which were found to be regulated by Mn toxicity in this study (Fig. 7). CHS, CHI and FLS transcripts are regulated by abiotic stresses, such as UV light, wounding and heavy metal toxicity [42,43,44]. For example, CHS is a key gene involved in regulating flavonoid biosynthesis; CHS expression is enhanced in Pb-exposed tomato [42]. Similarly, CHI expression is increased in pea roots (Pisum sativum) under Cd stress [45]. In addition, ZmFLS1 in maize is induced by UV-B radiation, and high expression levels of ZmFLS1 protect plants from UV-B stress [43]. In addition, DEGs belonging to isoflavonoid metabolism related proteins, such as HIDH, I2H, I7OMT and VRs, were also regulated by Mn stress (Fig. 7). Isoflavonoids are a predominant subclass of flavonoid metabolites that play critical roles in plant defense [46, 47].

As important osmotic regulators and ROS scavengers in plant stress tolerance, secondary metabolites are regulated by Mn stress in many plants, such as phenylalanine in Populus cathayana [13], lignin and flavonoid in rice [2], and phenol and callose in cowpea (Vigna unguiculata) [12]. Similarly, in this study, significant increases in secondary metabolites, including total phenols, flavonoids, tannins and anthocyanins, were observed in stylo leaves exposed to 400 and 800 μM Mn (Fig. 11). The accumulation of secondary metabolites might be attributed to the regulation of various genes, such as PAL, CCR, CAD, CHS, CHI, FLS and HIDH (Fig. 7). Cumulative results suggest that increasing gene expression can cause the accumulation of secondary metabolites. For example, increases in CHS and PAL transcripts are closely related to the increased levels of phenols, flavonoids and anthocyanins in tomato under Pb stress [42]. Furthermore, exogenous expression of ZmFLS1 from maize increases anthocyanin accumulation in an Arabidopsis fls1 mutant [43]. Although the exact roles of the candidate DEGs described above warrant further research, the results presented in this study strongly suggest that adjustment of secondary metabolism and corresponding gene expression might be important for stylo adaptation to Mn toxicity.

DEGs related to transcription factors

In this study, the identification of 123 transcription factors suggested complex regulation in stylo’s response to Mn toxicity (Fig. 8 and Additional file 5: Table S4). Among those identified transcription factors, seven DEGs belonging to the C2H2 zinc finger transcription factor family were up-regulated by Mn toxicity in stylo (Additional file 5: Table S4). C2H2 transcription factors are implicated in low pH-stress and aluminum (Al) toxicity tolerance in many plants, such as AtSTOP1/2 in Arabidopsis, VuSTOP1 in rice bean (Vigna umbellata) and GmSTOP1s in soybean [48,49,50,51]. For example, GmSTOP1–1/1–3 expression is enhanced by Al toxicity in soybean roots; complementation analyses of GmSTOP1s in an Arabidopsis Atstop1 mutant suggest that GmSTOP1s play critical roles in Al tolerance [51]. Furthermore, it has been demonstrated that C2H2 transcription factors function in Al detoxification mainly through activation of Al-tolerant genes, such as gene encoding glutamate dehydrogenase 1 (GDH1), pectin methylesterase inhibitor (PMI), malic enzyme (ME), aluminum-activated malate transporter (ALMT), multidrug and toxic compound exclusion protein (MATE) and tonoplast dicarboxylate transporter (TDT) [48, 49, 51]; these homologous genes were also up-regulated by Mn toxicity in stylo (Additional file 3: Table S2), suggesting that the functions of C2H2 transcription factors may be conserved in Mn detoxification through regulation of downstream critical Mn-resistant genes. Interestingly, two DEGs, CAD1/2, which belong to the phenylpropanoid biosynthesis pathway, were predicted to be regulated by six C2H2 transcription factors (Fig. 9), suggesting that these C2H2 transcription factors may also play transcriptional regulatory roles in secondary metabolism in stylo during Mn exposure, which deserves further study.

Moreover, nine transcription factor genes belonging to the WRKY family were up-regulated by Mn toxicity in stylo (Additional file 5: Table S4). Although the function of WRKY transcription factor remains largely unknown in stylo’s response to Mn toxicity, WRKY46 in Arabidopsis a negative regulator of ALMT1; disruption of WRKY46 results in increased malate exudation, conferring Al tolerance [52]. However, WRKY46 expression in Arabidopsis is suppressed by Al toxicity [52], but WRKY gene transcripts in stylo were enhanced by Mn stress in this study (Additional file 5: Table S4), suggesting that there are other regulatory networks that are mediated by WRKY-mediated gene expression in stylo during Mn toxicity.

Additionally, 22 transcription factors belonging to the MYB family were identified in this study (Additional file 5: Table S4). Cumulative studies have demonstrated that MYB family transcription factors participate in the response to environmental stresses in plants [53]. For example, MYB49 expression in Arabidopsis is induced by Cd treatment; overexpression of MYB49 in Arabidopsis results in a significant increase in Cd accumulation, inhibiting primary root growth [54]. Furthermore, MYB49 interacts with the transcription factors bHLH38 and bHLH101, leading to the activation of IRON-REGULATED TRANSPORTER1 and mediating Cd uptake [54]. In addition, OsARM1, which encodes the R2R3 MYB transcription factor, is induced by As treatment in rice. OsARM1 knockout improves tolerance to As stress in rice through regulation of As transporters, thereby mediating As uptake and translocation [55]. Furthermore, in this study, the MYB transcription factor is probably involved in the regulation of secondary metabolic pathways by mediating the transcripts of genes belonging to HIDH, anthocyanidin reductase and hypothetical protein (Fig. 9). Thus, future efforts are required to investigate the functions of these transcription factors in stylo subjected to Mn toxicity.

DEGs related to transporters

Different kinds of DEGs belonging to transporters possibly involved in Mn translocation, distribution and sequestration were identified in this study (Additional file 7: Table S6). Among them, a set of genes encoding homologs to ATP-binding cassette (ABC) transporters were regulated by Mn toxicity. Plant ABC transporters constitute a large superfamily and are implicated in the transport of a variety of substrates, including hormones, secondary metabolites and heavy metal ions [56]. It has been demonstrated that the transcripts of OsABCG43, belonging to the rice ABC transporter, is up-regulated by Cd treatment; overexpression of OsABCG43 confers Cd tolerance in the Cd-sensitive yeast mutant ycf1acr3, suggesting that the role of OsABCG43 in sequestration of Cd into subcellular organelles (e.g., vacuoles) [57]. OsALS1, encoding a half-size ABC transporter in rice, is found to be localized at the tonoplast and is required for Al detoxification through sequestration of Al into the vacuoles [58]. A similar result has been reported in Arabidopsis that AtALS1 participates in Al sequestration [59]. Furthermore, studies on cyanobacterium (Synechocystis) also suggest the potential role of ABC transporter in Mn transport [60]. Therefore, genes encoding ABC transporters are probably involved in stylo leaves defense against Mn toxicity.

Furthermore, genes encoding cation/H+ exchangers were found to be up-regulated by Mn toxicity in this study (Additional file 7: Table S6). The cation/H+ exchanger plays important roles in the vacuolar accumulation of metals, such as Mn. For example, the cation/H+ exchanger in Arabidopsis, AtCAX2, is involved in conferring tolerance to Mn toxicity when it is heterologously expressed in a Mn-sensitive yeast mutant pmc1vcx1cnb or in tobacco (Nicotiana tabacum) via sequestration of Mn into vacuoles [61]. Similar functions have been assigned to AtCAX4 and AtCAX5, which localize to the vacuolar membrane [62, 63]. In addition, one gene encoding homologs to aquaporin protein, aquaporin TIP1–3-like, was also up-regulated by Mn toxicity in stylo leaves (Additional file 7: Table S6). Genes belonging to the aquaporin protein family are reported to be involved in plants response to metal stresses. For example, two genes encoding tonoplast-intrinsic protein (TIP) are regulated by Al toxicity in buckwheat (Fagopyrum esculentum) leaves [64]. It has been demonstrated that HmVALT, belonging to the tonoplast-intrinsic protein in Al hyperaccumulating hydrangea (Hydrangea macrophylla), is a tonoplast-localized Al transporter that involves in vacuolar Al transport; overexpression of HmVALT confers Al tolerance in Arabidopsis [65]. There results suggest the potential roles of this gene in vacuolar Mn transport.

In addition, a group of genes belonging to transporters for various substrates, including phosphate, potassium and copper, were found to be differentially regulated by Mn in stylo leaves (Additional file 7: Table S6). The balance between the accumulation of Mn and other abundant mineral elements is thought to be important for regulating Mn homeostasis in plants [33]. Thus, these DEGs related to transporters responsible for Mn translocation, distribution and sequestration are regarded to be critical for intracellular Mn detoxification in stylo, and their functions warrant further investigation.

Conclusion

In conclusion, this study showed that stylo growth was inhibited by high Mn due to oxidative stress. Consistently, transcriptome analysis showed that a set of DEGs involved in regulating ROS scavenging and defense responses might contribute to stylo’s response to Mn toxicity. Furthermore, the identified DEGs associated with secondary metabolic pathways suggested that fine-tuned regulation of secondary metabolism might represent an adaptive strategy of stylo to cope with Mn toxicity. In addition, activation of various transcription factors may help stylo to tolerate Mn toxicity through transcriptional regulation. Taken together, this study not only reveals new insights into the molecular mechanisms underlying stylo’s response to Mn toxicity but also provides candidate genes that may be useful in the development of crop varieties with Mn adaptability through genetic engineering.

Methods

Plant growth and treatments

The stylo (Stylosanthes guianensis) genotype ‘TF303’ used in this study was provided by the Tropical Pasture Research Center, Institute of Tropical Crop Genetic Resources (TCGRI), Chinese Academy of Tropical Agriculture Sciences (CATAS), Hainan, China (voucher specimen deposited in the National Tropical Forage Germplasm Genebank, TCGRI, CATAS). The plant material was authenticated by Guodao Liu, TCGRI, CATAS. All experimental procedures were performed at the TCGRI, CATAS, Hainan, China (19°30′N, 109°30′E). After seed germination, stylo seedlings were transferred to a modified Hoagland nutrient solution (pH 5.8) containing 2500 μM KNO3, 2500 μM Ca (NO3)2, 400 μM NH4NO3, 250 μΜ KH2PO4, 500 μM MgSO4, 250 μM K2SO4, 5 μM MnSO4, 0.5 μM ZnSO4, 1.5 μM CuSO4, 0.09 μM (NH4)6Mo7O24, 23 μM NaB4O7 and 80 μM Fe-Na-EDTA [26]. After 30 d of normal growth, stylo seedlings were transplanted into fresh nutrient solution supplied with 5, 100, 200, 400 and 800 μM MnSO4 (pH 5.0). The background Mn level in the nutrient solution is 5 μM and thus used as the control [26]. After 10 d of Mn treatments, shoots and roots were harvested for further study. An individual hydroponic box containing four stylo seedlings served as one biological replicate. All treatments included three biological replicates.

Analysis of chlorophyll concentrations

To analyze chlorophyll concentrations, approximately 0.05 g leaf sample was incubated in 10 mL of 80% (v/v) acetone for 6 h. The chlorophyll concentration was spectrophotometrically detected at 440, 645 and 663 nm using a UV-2010 spectrophotometer (Hitachi, Japan) and determined as previously described [26].

Determination of Mn concentrations

To measure Mn concentrations, approximately 0.07 g dry shoot and root samples were transferred to a muffle furnace and burned at 600 °C for 10 h. Then, the sample was adequately extracted in 7 mL of 100 mM HCl. Mn concentrations in each sample were detected using atomic absorption spectroscopy [11].

Detection of O2− and MDA levels

O2− levels in stylo leaves were detected using the histochemical method [66]. Leaves subjected to different Mn treatments were incubated in 0.25 mM NBT solution for 12 h. Samples were then decolorized in boiling ethanol for 6 h. O2− levels were assessed in leaves based on the intensity of the blue insoluble formazan.

For MDA analysis, approximately 0.1 g samples were extracted in 1 mL of 10% trichloracetic acid (TCA). Homogenates were centrifuged at 12,000×g for 10 min. The supernatants were collected, and 1 mL supernatant was added to 1 mL of 0.5% thiobarbituric acid (TBA). The mixture was then incubated at 100 °C for 20 min. Then, homogenates were centrifuged at 12,000×g for 10 min. The MDA concentration was detected as previously described [67].

Analysis of electrolyte leakage and Fv/Fm

Stylo leaves treated with various Mn concentrations were separately immersed in 10 mL of deionized water in a glass test tube overnight at room temperature. The initial conductivity (IC) of the solution was detected using a conductimeter (INESA Scientific Instrument Co., Ltd., Shanghai, China). Total conductivity (TC) of the solution was detected after boiling for 15 min. Relative electrolyte leakage was expressed as the percentage of the initial conductivity (IC) to the total conductivity (TC). Fv/Fm values of stylo leaves were measured using a pulse-modulated fluorometer Model FMS-2 (Hansatech Instruments Ltd., UK) according to the manufacturer’s instructions.

Analysis of enzyme activity

To determine POD activity, approximately 0.1 g leaf samples from different Mn treatments were extracted in 2 mL of 50 mM phosphate buffer (pH 7.8) at 4 °C. The extracts were centrifuged at 12,000×g for 15 min at 4 °C, and the supernatants were recovered for POD assay [66]. APX was extracted by grinding 0.1 g leaves in 2 mL of 50 mM phosphate buffer (pH 7.0) supplied with 1 mM AsA and 1 mM ethylenediaminetetraacetic acid (EDTA) at 4 °C. After centrifugation, the supernatants were collected for APX activity assays [66].

To detect PPO activity, approximately 0.1 g leaf samples were extracted in 1 mL of extraction buffer obtained from commercial assay kits (Solarbio Science and Technology Co., Ltd., Beijing, China) at 4 °C. The extracts were centrifuged at 12,000×g for 15 min at 4 °C, and supernatants were collected for PPO activity assessment at 410 nm. PPO activity was calculated according to the manufacturer’s protocol. To analyze PAL activity, 0.1 g leaf samples were extracted in 0.9 mL of extraction buffer obtained from commercial assay kits (Nanjing Jiancheng Bioengineering Institute, Jiangsu, China) at 4 °C. After centrifugation, the supernatants were used to determine PAL activity at 290 nm according to the protocol.

Transcriptome analysis of stylo leaves

RNA sequencing was performed by Annoroad Gene Technology Co., Ltd. (Beijing, China). Total RNA of stylo leaves subjected to control (5 μM) and Mn toxicity (400 μM) treatments was isolated using TRIzol reagent (Invitrogen, Carlsbad, CA, USA). RNA integrity and purity were analyzed using an Agilent 2100 RNA Nano 6000 Assay Kit (Agilent Technologies, Palo Alto, CA, USA). Sequencing libraries were generated using the NEBNext Ultra Directional RNA Library Prep kit (NEB, Beverly, MA, USA). The cDNA libraries of the three biological replicates were obtained using an Illumina Hiseq™ 4000 platform (Illumina, San Diego, CA, USA). PE150 was used as the sequencing program.

Raw image data were transformed to raw reads using the Consensus Assessment of Sequence and Variation (CASAVA, v.1.8.2) base recognition program. The raw reads were stored in FASTQ files. High-quality clean reads were obtained by filtering out adaptor sequences and low-quality reads. The obtained clean reads were further assembled into unigenes using Trinity software (version 2.4.0) [68]. The clean reads were mapped back to the assembled contigs through Bowtie2 v.2.2.3 software [69]. For annotation, all unigene sequences were aligned by BLAST search (E-value<1e-5) against the NCBI non-redundant protein sequences (Nr), NCBI non-redundant nucleotide sequences (Nt), Universal Protein (UniProt), Gene Ontology (GO), the Clusters of Orthologous Groups of protein database (COG), the Functional Annotation (eggNOG), the KEGG database and the Protein Family (Pfam).

The transcript abundance of each gene in each sample was calculated and represented by the reads per kilobase million mapped reads (RPKM) value using RNA-seq data with expectation maximization (RSEM) [70]. DESeq2 was used to calculate the differences in the expression abundance between 5 and 400 μM Mn treatments [71]. Genes with q < 0.05 and |log2(fold change)| ≥1 were considered DEGs. Subsequently, DEGs were extracted for GO and KEGG enrichment analysis [72, 73]. Protein-protein interaction network analysis for the identified DEGs was performed using STRING (version 11.0) (https://string-db.org/). The interaction networks were constructed using Cytoscape (version 3.8.0) (https://github.com/cytoscape/cytoscape) with default settings. Nodes and edges indicate genes and interactions, respectively. The sequencing data were deposited in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus under GEO series number GSE147431.

Quantitative real-time PCR (qRT-PCR) analysis

Total RNA from leaves was extracted using the TRIzol reagent (Invitrogen, USA) according to the manual. The first cDNA strand was synthesized from 2 μg of DNase I-treated RNA using oligo d(T), dNTPs, RNase inhibitor and reverse transcriptase in the HiScript III cDNA synthesis kit (Vazyme, China) according to the protocol. qRT-PCR analysis was performed in a 20 μL volume containing 2 μL 1:50 diluted cDNA, 0.2 μM gene primers, 6.4 μL ddH2O and 10 μL SYBR Green Master mix (Vazyme, China). The reaction was performed on a QuantStudio™ 6 Flex Real-Time System (Thermo Fisher Scientific, Waltham, MA, USA). qRT-PCR analysis conditions were as follows: 95 °C for 1 min, 40 cycles of 95 °C for 15 s, 58 °C for 15 s and 72 °C for 30 s. qRT-PCR primers of genes are listed in Additional file 8: Table S7. The housekeeping gene SgEF-1a (Accession No. JX164254) was used as an endogenous control to calculate the relative gene expression as the ratio of the transcript of the candidate gene to that of SgEF-1a [26]. Gene expression analysis included three biological replicates.

Determination of secondary metabolite concentrations

The concentrations of secondary metabolites, including total phenols, flavonoids, anthocyanins and tannins, in stylo leaves exposed to various Mn treatments were analyzed using commercial assay kits (Nanjing Jiancheng Bioengineering Co., Ltd., Jiangsu, China). To detect total phenols, flavonoids and anthocyanins, 0.1 g dry leaf samples were extracted and measured as previously described [66]. For tannin assays, approximately 0.1 g fresh leaf samples were extracted in 1 mL of ddH2O. The concentration of tannins was detected according to the manufacturer’s protocol.

Statistical analysis

Microsoft Excel 2003 (Microsoft Company, USA) was used in the data analysis. One-way analysis of variance (ANOVA) was performed using SPSS software, version 13.0 (SPSS Institute, Chicago, USA).

Availability of data and materials

The datasets for this study can be found in the NCBI Gene Expression Omnibus under GEO series number of GSE147431.

Abbreviations

APX:

Ascorbate peroxidase

CAD:

Cinnamyl-alcohol dehydrogenase

CALDH:

Coniferyl-aldehyde dehydrogenase

CCR:

Cinnamoyl-CoA reductase

CHI:

Chalcone isomerase

CHS:

Chalcone synthase

COMT:

Caffeoyl-CoA O-methyltransferase

CSE:

Caffeoylshikimate esterase

DEGs:

Differentially expressed genes

F3’5’H:

Flavonoid 3,5-hydroxylase

F6H:

Feruloyl-CoA 6-hydroxylase

FLS:

Flavonol synthase

GLU:

Beta-glucosidase

GO:

Gene Ontology

GST:

Glutathione S-transferase

HI4OMT:

2,7,4-Trihydroxyisoflavanone 4-O-methyltransferase/isoflavone 4-O-methyltransferase

HIDH:

2-Hydroxyisoflavanone dehydratase

I2H:

Isoflavone 2-hydroxylase

I7OMT:

sofilavone-7-O-methyltransferase

KEGG:

Kyoto Encyclopedia of Genes and Genomes

MDA:

Malondialdehyde

NBT:

Nitroblue tetrazolium chloride

PAL:

Phenylalanine ammonia-lyase

POD:

Peroxidase

PPO:

Polyphenol oxidase

qRT-PCR:

Quantitative real-time polymerase chain reaction

RNA-seq:

RNA sequencing

ROS:

Reactive oxygen species

SOD:

Superoxide dismutase

SPL19:

Serine carboxypeptidase-like 19

TFs:

Transcription factors

VR:

Vestitone reductase

References

  1. Millaleo R, Reyes-Díaz M, Ivanov AG, Mora ML, Alberdi M. Manganese as essential and toxic element for plants: transport, accumulation and resistance mechanism. J Soil Sci Plant Nutr. 2010;10:470–81.

    Article  Google Scholar 

  2. Lidon FC, Barreiro MG, Ramalho JC. Manganese accumulation in rice: implications for photosynthetic functioning. J Plant Physiol. 2004;161:1235–44.

    Article  CAS  PubMed  Google Scholar 

  3. Shao JF, Yamaji N, Shen RF, Ma JF. The key to Mn homeostasis in plants: regulation of Mn transporters. Trends Plant Sci. 2017;22:215–24.

    Article  CAS  PubMed  Google Scholar 

  4. Geszvain K, Butterfield C, Davis RE, Madison AS, Lee SW, Parker DL, et al. The molecular biogeochemistry of manganese (II) oxidation. Biochem Soc Trans. 2012;40:1244–8.

    Article  CAS  PubMed  Google Scholar 

  5. Sparrow LA, Uren NC. Manganese oxidation and reduction in soils: effects of temperature, water potential, pH and their interactions. Soil Res. 2014;52:483–94.

    Article  CAS  Google Scholar 

  6. Watmough SA, Eimers MC, Dillon PJ. Manganese cycling in Central Ontario forests: response to soil acidification. Appl Geochem. 2007;22:1241–7.

    Article  CAS  Google Scholar 

  7. Hue NV, Mai Y. Manganese toxicity in watermelon as affected by lime and compost amended to a Hawaiian acid Oxisol. HortScience. 2002;37:656–61.

    Article  CAS  Google Scholar 

  8. Fernando DR, Lynch JP. Manganese phytotoxicity: new light on an old problem. Ann Bot. 2015;116:313–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Li JF, Jia YD, Dong RS, Huang R, Liu PD, Li XY, et al. Advances in the mechanisms of plant tolerance to manganese toxicity. Int J Mol Sci. 2019;20:5096.

    Article  CAS  PubMed Central  Google Scholar 

  10. Ribera-Fonseca A, Inostroza-Blancheteau C, Cartes P, Rengel Z, Mora ML. Early induction of Fe-SOD gene expression is involved in tolerance to Mn toxicity in perennial ryegrass. Plant Physiol Biochem. 2013;73:77–82.

    Article  CAS  PubMed  Google Scholar 

  11. Liu P, Huang R, Hu X, Jia Y, Li J, Luo J, et al. Physiological responses and proteomic changes reveal insights into Stylosanthes response to manganese toxicity. BMC Plant Biol. 2019;19:212.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Fecht-Christoffers MM, Führs H, Braun HP, Horst WJ. The role of H2O2-producing and H2O2-consuming peroxidases in the leaf apoplast of Vigna unguiculata L. in manganese tolerance. Plant Physiol. 2006;140:1451–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Lei Y, Korpelainen H, Li C. Physiological and biochemical responses to high Mn concentrations in two contrasting Populus cathayana populations. Chemosphere. 2007;68:686–94.

    Article  CAS  PubMed  Google Scholar 

  14. Pittman JK. Managing the manganese: molecular mechanisms of manganese transport and homeostasis. New Phytol. 2005;167:733–42.

    Article  CAS  PubMed  Google Scholar 

  15. Socha AL, Guerinot ML. Mn-euvering manganese: the role of transporter gene family members in manganese uptake and mobilization in plants. Front Plant Sci. 2014;5:106.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Yamaji N, Sasaki A, Xia JX, Yokosho K, Ma JF. A node-based switch for preferential distribution of manganese in rice. Nat Commun. 2013;4:2442.

    Article  PubMed  CAS  Google Scholar 

  17. Delhaize E, Kataoka T, Hebb DM, White RG, Ryan PR. Genes encoding proteins of the cation diffusion facilitator family that confer manganese tolerance. Plant Cell. 2003;15:1131–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Chen Z, Fujii Y, Yamaji N, Masuda S, Takemoto Y, Kamiya T, et al. Mn tolerance in rice is mediated by MTP8.1, a member of the cation diffusion facilitator family. J Exp Bot. 2013;64:4375–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Takemoto Y, Tsunemitsu Y, Fujii KM, Mitani UN, Yamaji N, Ma JF, et al. The tonoplast-localized transporter MTP8.2 contributes to manganese detoxification in the shoots and roots of Oryza sativa L. Plant Cell Physiol. 2017;58:1573–82.

    Article  CAS  PubMed  Google Scholar 

  20. Eroglu S, Meier B, Wirén NV, Peiter E. The vacuolar manganese transporter MTP8 determines tolerance to iron deficiency-induced chlorosis in Arabidopsis. Plant Physiol. 2016;170:1030–45.

    Article  CAS  PubMed  Google Scholar 

  21. Migocka M, Papierniak A, Dziubińska EM, Poździk P, Posyniak E, Garbiec A, et al. Cucumber metal transport protein MTP8 confers increased tolerance to manganese when expressed in yeast and Arabidopsis thaliana. J Exp Bot. 2014;65:5367–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Li Q, Li Y, Wu X, Zhou L, Zhu X, Fang W. Metal transport protein 8 in Camellia sinensis confers superior manganese tolerance when expressed in yeast and Arabidopsis thaliana. Sci Rep. 2017;7:39915.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Schultze-Kraft R, Rao IM, Peters M, Clements RJ, Bai C, Liu GD. Tropical forage legumes for environmental benefits: an overview. Trop Grassl Forrajes. 2018;6:1–14.

    Article  Google Scholar 

  24. Guo PF, Liu PD, Lei J, Chen CH, Qiu H, Liu GD, et al. Improvement of plant regeneration and Agrobacterium-mediated genetic transformation of Stylosanthes guianensis. Trop Grassl Forrajes. 2019;7:480–92.

    Article  Google Scholar 

  25. Andrew CS, Hegarty MP. Comparative responses to manganese excess of eight tropical and four temperate pasture legume species. Aust J Agvic Res. 1969;20:687–96.

    Article  Google Scholar 

  26. Chen ZJ, Sun LL, Liu PD, Liu GD, Tian J, Liao H. Malate synthesis and secretion mediated by a manganese-enhanced malate dehydrogenase confers superior manganese tolerance in Stylosanthes guianensis. Plant Physiol. 2015;167:176–88.

    Article  CAS  PubMed  Google Scholar 

  27. Nagajyoti PC, Lee KD, Sreekanth TVM. Heavy metals, occurrence andtoxicity for plants: a review. Environ Chem Lett. 2010;8:199–216.

    Article  CAS  Google Scholar 

  28. Lucchini R, Placidi D, Cagna G, Fedrighi C, Oppini M, Peli M, et al. Manganese and developmental neurotoxicity. Adv Neurobiol. 2017;18:13–34.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Zhang P, Ding Z, Zhong Z, Tong H. Transcriptomic analysis for Indica and Japonica rice varieties under aluminum toxicity. Int J Mol Sci. 2019;20:997.

    Article  CAS  PubMed Central  Google Scholar 

  30. Jiang W, Wu Z, Wang T, Mantri N, Huang H, Li H, et al. Physiological and transcriptomic analyses of cadmium stress response in Dendrobium officinale seedling. Plant Physiol Biochem. 2020;148:152–65.

    Article  CAS  PubMed  Google Scholar 

  31. Wang YM, Yang Q, Xu H, Liu YJ, Yang HL. Physiological and transcriptomic analysis provide novel insight into cobalt stress responses in willow. Sci Rep. 2020;10:2308.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Rosas A, Rengel Z, Mora ML. Manganese supply and pH influence growth, carboxylate exudation and peroxidase activity of ryegrass and white clover. J Plant Nutr. 2007;30:253–70.

    Article  CAS  Google Scholar 

  33. Sheng H, Zeng J, Yan F, Wang X, Wang Y, Kang H. Effect of exogenous salicylic acid on manganese toxicity, mineral nutrients translocation and antioxidative system in polish wheat (Triticum polonicum L.). Acta Physiol Plant. 2015;37:32.

    Article  CAS  Google Scholar 

  34. Zhou CP, Qi YP, You X, Yang LT, Guo P, Ye X, et al. Leaf cDNA-AFLP analysis of two citrus species differing in manganese tolerance in response to long-term manganese-toxicity. BMC Genomics. 2013;14:621.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Kumar S, Asif MH, Chakrabarty D, Tripathi RD, Dubey RS, Trivedi PK. Expression of a rice lambda class of glutathione S-transferase, OsGSTL2, in Arabidopsis provides tolerance to heavy metal and other abiotic stresses. J Hazard Mater. 2013;248–249:228–37.

    Article  PubMed  CAS  Google Scholar 

  36. Dmitriev AA, Krasnov GS, Rozhmina TA, Kishlyan NV, Zyablitsin AV, Sadritdinova AF, et al. Glutathione S-transferases and UDP-glycosyltransferases are involved in response to aluminum stress in flax. Front Plant Sci. 2016;7:1920.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Mészáros P, Rybanský L, Spieß N, Socha P, Kuna R, Libantová J, et al. Plant chitinase responses to different metal-type stresses reveal specificity. Plant Cell Rep. 2014;33:1789–99.

    Article  PubMed  CAS  Google Scholar 

  38. Liu Y, Li D, Song Q, Zhang T, Li D, Yang X. The maize late embryogenesis abundant protein ZmDHN13 positively regulates copper tolerance in transgenic yeast and tobacco. Crop J. 2019;7:403–10.

    Article  Google Scholar 

  39. Cochrane FC, Davin LB, Lewis NG. The Arabidopsis phenylalanine ammonialyase gene family: kinetic characterization of the four PAL isoforms. Phytochemistry. 2004;65:1557–64.

    Article  CAS  PubMed  Google Scholar 

  40. Sharma A, Shahzad B, Rehman A, Bhardwaj R, Landi M, Zheng B. Response of phenylpropanoid pathway and the role of polyphenols in plants under abiotic stress. Molecules. 2019;24:2452.

    Article  CAS  PubMed Central  Google Scholar 

  41. Pawlak-Sprada S, Arasimowicz-Jelonek M, Podgórska M, Deckert J. Activation of phenylpropanoid pathway in legume plants exposed to heavy metals. Part I. effects of cadmium and lead on phenylalanine ammonia-lyase gene expression, enzyme activity and lignin content. Acta Biochim Pol. 2011;58:211–6.

    CAS  PubMed  Google Scholar 

  42. Bali S, Jamwal VL, Kohli SK, Kaur P, Tejpal R, Bhalla V, et al. Jasmonic acid application triggers detoxification of lead (Pb) toxicity in tomato through the modifications of secondary metabolites and gene expression. Chemosphere. 2019;235:734–48.

    Article  CAS  PubMed  Google Scholar 

  43. Falcone-Ferreyra ML, Rius SP, Emiliani J, Pourcel L, Feller A, Morohashi K, et al. Cloning and characterization of a UV-B-inducible maize flavonol synthase. Plant J. 2010;62:77–91.

    Article  PubMed  CAS  Google Scholar 

  44. Dao TTH, Linthorst HJM, Verpoorte R. Chalcone synthase and its functions in plant resistance. Phytochem Rev. 2011;10:397–412.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Rivera-Becerril F, Metwally A, Martin-Laurent F, Tuinen DV, Dietz KJ, Gianinazzi S, et al. Molecular responses to cadmium in roots of Pisum sativum L. Water Air Soil Poll. 2005;168:171–86.

    Article  CAS  Google Scholar 

  46. Dixon RA, Steele CL. Flavonoids and isoflavonoids-a gold mine for metabolic engineering. Trends Plant Sci. 1999;4:394–400.

    Article  CAS  PubMed  Google Scholar 

  47. Pawlak-Sprada S, Stobiecki M, Deckert J. Activation of phenylpropanoid pathway in legume plants exposed to heavy metals. Part II. Profiling of iso-flavonoids and their glycoconjugates induced in roots of lupine (Lupinus luteus) seedlings treated with cadmium and lead. Acta Biochim Pol. 2011;58:217–23.

    CAS  PubMed  Google Scholar 

  48. Sawaki Y, Iuchi S, Kobayashi Y, Kobayashi Y, Ikka T, Sakurai N, et al. STOP1 regulates multiple genes that protect Arabidopsis from proton and aluminum toxicities. Plant Physiol. 2009;150:281–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Kobayashi Y, Ohyama Y, Kobayashi Y, Ito H, Iuchi S, Fujita M, et al. STOP2 activates transcription of several genes for Al and low pH tolerance that are regulated by STOP1 in Arabidopsis. Mol Plant. 2014;7:311–22.

    Article  CAS  PubMed  Google Scholar 

  50. Fan W, Lou HQ, Gong YL, Liu MY, Cao MJ, Liu Y, et al. Characterization of an inducible C2H2-type zinc finger transcription factor VuSTOP1 in rice bean (Vigna umbellata) reveals differential regulation between low pH and aluminum tolerance mechanisms. New Phytol. 2015;208:456–68.

    Article  CAS  PubMed  Google Scholar 

  51. Wu W, Lin Y, Chen Q, Peng W, Peng J, Tian J, et al. Functional conservation and divergence of soybean GmSTOP1 members in proton and aluminum tolerance. Front Plant Sci. 2018;9:570.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Ding ZJ, Yan JY, Xu XY, Li GX, Zheng SJ. WRKY46 functions as a transcriptional repressor of ALMT1, regulating aluminum induced malate secretion in Arabidopsis. Plant J. 2013;76:825–35.

    Article  CAS  PubMed  Google Scholar 

  53. Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C, Lepiniec L. MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010;15:573–81.

    Article  CAS  PubMed  Google Scholar 

  54. Zhang P, Wang R, Ju Q, Li W, Tran LS, Xu J. The R2R3-MYB transcription factor MYB49 regulates cadmium accumulation. Plant Physiol. 2019;180:529–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Wang FZ, Chen MX, Yu LJ, Xie LJ, Yuan LB, Qi H, et al. OsARM1, an R2R3 MYB transcription factor, is involved in regulation of the response to arsenic stress in rice. Front Plant Sci. 2017;8:1868.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Jasinski M. The ATP-binding cassette transporters: structure, function, and gene family comparison between rice and Arabidopsis. Plant Physiol. 2003;131:1169–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Oda K, Otani M, Uraguchi S, Akihiro T, Fujiwara T. Rice ABCG43 is cd inducible and confers cd tolerance on yeast. Biosci Biotechnol Biochem. 2011;75:1211–3.

    Article  CAS  PubMed  Google Scholar 

  58. Huang CF, Yamaji N, Chen Z, Ma JF. A tonoplast-localized half-size ABC transporter is required for internal detoxification of aluminum in rice. Plant J. 2012;69:857–67.

    Article  CAS  PubMed  Google Scholar 

  59. Larsen PB, Cancel J, Rounds M, Ochoa V. Arabidopsis ALS1 encodes a root tip and stele localized half type ABC transporter required for root growth in an aluminum toxic environment. Planta. 2007;225:1447–58.

    Article  CAS  PubMed  Google Scholar 

  60. Bartsevich V, Pakrasi HB. Manganese transport in the cyanobacterium Synechocystis sp PCC 6803. J Biol Chem. 1996;271:26057–61.

    Article  CAS  PubMed  Google Scholar 

  61. Hirschi KD, Korenkov VD, Wilganowski NL, Wagner GJ. Expression of Arabidopsis CAX2 in tobacco. Altered metal accumulation and increased manganese tolerance. Plant Physiol. 2000;124:125–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Cheng NH, Pittman JK, Shigaki T, Hirschi KD. Characterization of CAX4, an Arabidopsis H+/cation antiporter. Plant Physiol. 2002;128:1245–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Edmond C, Shigaki T, Ewert S, Nelson MD, Connorton JM, Chalova V, et al. Comparative analysis of CAX2-like cation transporters indicates functional and regulatory diversity. Biochem J. 2009;418:145–54.

    Article  CAS  PubMed  Google Scholar 

  64. Chen WW, Xu JM, Jin JF, Lou HQ, Fan W, Yang JL. Genome-wide transcriptome analysis reveals conserved and distinct molecular mechanisms of Al resistance in buckwheat (Fagopyrum esculentum Moench) leaves. Int J Mol Sci. 2017;18:1859.

    Article  PubMed Central  CAS  Google Scholar 

  65. Negishi T, Oshima K, Hattori M, Kana M, Mano S, Nishimura M, et al. Tonoplast-and plasma membrane-localized aquaporin-family transporters in blue hydrangea sepals of aluminum hyperaccumulating plant. PLoS One. 2012;7:e43189.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Wang H, Chen ZJ, Liu GD, Bai CJ, Qiu H, Jia YX, et al. Alterations of growth, antioxidant system and gene expression in Stylosanthes guianensis during Colletotrichum gloeosporioides infection. Plant Physiol Biochem. 2017;118:256–66.

    Article  CAS  PubMed  Google Scholar 

  67. Guo Z, Huang M, Lu S, Zhao Y, Zhong Q. Differential response to paraquat induced oxidative stress in two rice cultivars on antioxidants and chlorophyll a fluorescence. Acta Physiol Plant. 2007;29:39–46.

    Article  CAS  Google Scholar 

  68. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:25–34.

    Article  CAS  Google Scholar 

  70. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011;12:323.

    Article  CAS  Google Scholar 

  71. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  72. Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Zhao N, Li W, Bai S, Guo W, Yuan G, Wang F, et al. Transcriptome profiling to identify genes involved in mesosulfuron-methyl resistance in Alopecurus aequalis. Front Plant Sci. 2017;8:1391.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank Mr. Chun Liu for providing critical comments about this work. We sincerely thank American Journal Experts (https://www.aje.com/) for critical review on this manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (31801951, 31672483, 31861143013), the Central Public-interest Scientific Institution Basal Research Fund for CATAS (1630032020004, 1630032019002, 1630032018004), the Modern Agro-industry Technology Research System (CARS-34), and the Integrated Demonstration of Key Techniques for the Industrial Development of Featured Crops in Rocky Desertification Areas of Yunnan-Guangxi-Guizhou Provinces (SMH2019–2021). The funding bodies had no roles in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

ZC and LL conceived and designed the research. YJ, XL, JL and RD performed the physiological experiments. YJ, XL and QL analyzed the RNA-Seq data. XH and PL carried out genes expression analysis. GL provided and prepared the plant material as well as analyzed and interpreted the results for this work. ZC, YJ and XL wrote the manuscript. GL, LL and ZC substantially revised the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Lijuan Luo or Zhijian Chen.

Ethics declarations

Ethics approval and consent to participate

This study did not involve field experiment. The stylo plant used in this study was provided by the Tropical Pasture Research Center, Institute of Tropical Crop Genetic Resources (TCGRI), Chinese Academy of Tropical Agriculture Sciences (CATAS), Hainan, China. All experimental procedures concerning on the plant were conducted in accordance with the institutional and national guidelines of China. This research did not include experiments concerning on human participants or animals.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Figure S1.

Effects of different Mn treatments on relative electrolyte leakage, Fv/Fm, plant height and total protein content.

Additional file 2: Table S1.

Summary of stylo leave transcriptomes in control (5 μM) and Mn toxicity (400 μM) libraries.

Additional file 3: Table S2.

Information of DEGs identified in stylo leaves subjected to control (5 μM) and Mn toxicity (400 μM) treatments.

Additional file 4: Table S3.

DEGs involved in secondary metabolism.

Additional file 5: Table S4.

Transcription factors responding to Mn toxicity.

Additional file 6: Table S5.

DEGs involved in defense responses.

Additional file 7: Table S6.

DEGs related to transporters.

Additional file 8: Table S7.

Primers used for qRT-PCR analysis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Jia, Y., Li, X., Liu, Q. et al. Physiological and transcriptomic analyses reveal the roles of secondary metabolism in the adaptive responses of Stylosanthes to manganese toxicity. BMC Genomics 21, 861 (2020). https://doi.org/10.1186/s12864-020-07279-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-020-07279-2

Keywords