- © 2017 American Society of Plant Biologists. All Rights Reserved.
Abstract
Arbuscular mycorrhizal (AM) associations enhance the phosphorous and nitrogen nutrition of host plants, but little is known about their role in potassium (K+) nutrition. Medicago truncatula plants were cocultured with the AM fungus Rhizophagus irregularis under high and low K+ regimes for 6 weeks. We determined how K+ deprivation affects plant development and mineral acquisition and how these negative effects are tempered by the AM colonization. The transcriptional response of AM roots under K+ deficiency was analyzed by whole-genome RNA sequencing. K+ deprivation decreased root biomass and external K+ uptake and modulated oxidative stress gene expression in M. truncatula roots. AM colonization induced specific transcriptional responses to K+ deprivation that seem to temper these negative effects. A gene network analysis revealed putative key regulators of these responses. This study confirmed that AM associations provide some tolerance to K+ deprivation to host plants, revealed that AM symbiosis modulates the expression of specific root genes to cope with this nutrient stress, and identified putative regulators participating in these tolerance mechanisms.
Potassium (K+) is a macronutrient required by all organisms for growth and development. In plants, K+ participates in many diverse processes, including plasma membrane polarization, growth, stomatal aperture regulation, and the acquisition of other nutrients (Wang and Wu, 2013). K+ represents 2% to 10% of the plant dry biomass, and its optimal cytoplasmic concentration for enzyme activities is around 100 to 200 mm (Leigh and Wyn Jones, 1984). Thus, a high cellular K+ concentration is vital for plants to maintain a wide range of physiological processes (Adams and Shin, 2014; Benito et al., 2014). Depending on soil type, the concentration of available K+ varies from 0.1 to 1 mm (Asher and Ozanne, 1967). Only two fractions of K+ are immediately accessible to plants: K+ in water solution and K+ in exchangeable form. These two fractions represent about 2.2% of the total amount of soil K+ (Zörb et al., 2014). This very low availability combined with the constitutive demand of plants leads to the formation of depletion zones near the root surface (Drew and Nye, 1969). Plants have developed efficient strategies to deal with limited K+ availability, including the establishment of symbiotic associations with microbes.
Arbuscular mycorrhizal (AM) fungi colonize the root cortex intercellularly and intracellularly (Brundrett, 2004; Harrison, 2005). These AM fungi belong to the fungal subphylum of Glomeromycotina and probably played a critical role in the emergence and evolution of land plants (Delaux et al., 2015; Spatafora et al., 2016). This AM symbiosis is the oldest mycorrhizal association and is found in about 65% of extant land plants, including many major crop species (Wang and Qiu, 2006). AM associations are characterized by the formation of arbuscules in the cortical cells of host plants, which are the site of nutrient exchange between the plant and the fungus (Smith and Read, 2008). The improvement of plant nutrient acquisition through AM associations becomes particularly important when hydromineral resources are scarce. In return for the nutrients provided by AM fungi, host plants provide up to 20% to 25% of their photosynthesis-derived carbohydrates to their symbionts (López et al., 2008).
The improvement of plant nutrition through AM symbiosis and the molecular basis of nutrient transfer are well characterized for phosphorus, nitrogen, and sulfur (Govindarajulu et al., 2005; Javot et al., 2007; Jin et al., 2012; Casieri et al., 2013; Courty et al., 2016; Garcia et al., 2016). However, the role of mycorrhizal symbioses in plant K+ nutrition is still understudied and misunderstood (Garcia and Zimmermann, 2014). A small number of studies reported AM fungal mediation of plant K+ nutrition. The assessment of K+ distribution in mycorrhizal plants using particle-induced x-ray emission experiments revealed a higher K+ concentration in root sections of Aster tripolium colonized by Rhizophagus irregularis than in noninoculated plants (Scheloske et al., 2004). The same trend of K+ enrichment for plants mycorrhized by various AM fungi also was observed in maize (Zea mays) root steles (Kaldorf et al., 1999), Pelargonium peltatum shoots (Perner et al., 2007), Lactuca sativa leaves (Baslam et al., 2013), and wheat (Triticum aestivum) stems (Oliveira et al., 2016). The expression of three K+ transport systems likely involved in phloem loading/unloading (ZmAKT2), xylem release (ZmSKOR), and sodium/K+ homeostasis (ZmSOS1) was differentially regulated in AM maize roots (Estrada et al., 2013). These findings suggested that AM symbioses could play a significant role in plant K+ nutrition by the reorganization of molecular responses.
If K+ amendments are not used frequently, long-term K+ deprivation can result in irreversible damage to the crop, including poor pest tolerance and significant yield losses. Soils from many parts of the world are increasingly K+ deficient, reinforcing the need for understanding how plants can cope with K+ deprivation (Moody and Bell, 2006; Römheld and Kirkby, 2010; Zörb et al., 2014). However, most studies have only investigated the short-term molecular plant responses to K+ deprivation, and not in the context of AM symbiosis. The transcriptional responses of plant roots to K+ deprivation were studied in various species, including Arabidopsis (Arabidopsis thaliana; Armengaud et al., 2004), rice (Oryza sativa; Ma et al., 2012; Shankar et al., 2013), wheat (Ruan et al., 2015), and sugarcane (Saccharum officinarum; Zeng et al., 2015). However, in all these studies, K+ was completely absent from the medium and the plants were analyzed after a relatively short period of time (1–2 weeks). Following these drastic and short-term treatments, the up-regulation of high-affinity K+ transporters often was reported. This response is a rapid and transient strategy but probably not a sustainable one if the stress is prolonged. The transcriptional response of plant roots colonized by an AM fungus during K+ deprivation is still unknown.
To address this question, we examined, to our knowledge for the first time, the effect of long-term K+ deprivation on the development, nutrition, gene expression, and AM formation using the model legume Medicago truncatula. To evaluate the impact of AM symbiosis on plant adaptation to long-term K+ deprivation, we inoculated M. truncatula plants with the AM fungus R. irregularis and watered them with high- or low-K+ solution. Plant biomass and the tissue content of various ions were determined after 6 weeks. The transcriptional profile of AM roots under K+ deprivation was investigated using RNA sequencing (RNA-seq) and revealed mycorrhiza-specific responses under K+ deprivation. Finally, a network-based prioritization analysis revealed putative regulators controlling AM symbiosis and K+ deprivation in M. truncatula.
RESULTS
AM Symbiosis Improves the K+ Acquisition of M. truncatula
The impact of AM associations on K+ nutrition is understudied and, in particular, was never investigated in legumes. M. truncatula plants were inoculated with spores of R. irregularis (AM) or were kept nonmycorrhized (NM) for 6 weeks. All plants were watered with either high- or low-K+ solution as described previously, and various physiological parameters were measured (Fig. 1). The quantification of fungal colonization revealed that around 30% of plant roots were colonized under both K+-sufficient and K+-deficient conditions, indicating that K+ availability did not affect the establishment of AM symbiosis (Fig. 1A). The expression of four AM-specific marker genes identified by RNA-seq was determined in plants growing at high or low K+. Transcripts coding for PHOSPHATE TRANSPORTER4 (PT4), REDUCED ARBUSCULAR MYCORRHIZA1 (RAM1) and RAM2, and STUNTED ARBUSCULE (STR) were up-regulated only in AM roots at high and low K+ levels, indicating that the symbiosis is functional (Fig. 1B). Interestingly, the up-regulation of these transcripts was more pronounced in M. truncatula roots at high K+ than at low K+, suggesting a slight alteration of the AM symbiosis under K+ deprivation.
Impact of AM symbiosis on plant fitness under K+ deprivation. A and B, The rate of fungal colonization (A) and the fragments per kilobase of transcript per million mapped reads (FPKM) of PT4, RAM1, RAM2, and STR (B) were measured on M. truncatula roots grown under high-K+ (+K) or low-K+ (−K) conditions after 6 weeks of coculture using the gridline intersection method. C, Shoot and root dry weight (DW) were measured for AM and NM plants grown under +K and −K conditions. D and E, K+ (D) and sodium (E) contents were determined by ICP-OES analysis in shoots of AM and NM plants under +K and −K conditions. Significant differences were obtained using Student’s t test between AM and NM plants (C–E) and between −K and +K conditions (C; *, P < 0.05 and **, P < 0.01). n = 6. Experiments were replicated three times independently.
Although no significant difference was observed in shoot dry biomass between AM and NM plants at high and low K+, AM plants displayed a significantly higher root dry biomass than NM plants under K+ deprivation (Fig. 1C). Sodium and K+ contents in plant roots and shoots were determined by inductively coupled plasma optical emission spectrometry (ICP-OES). The impossibility of separating the fungus from colonized tissues resulted in noncomparable ion content determination between AM and NM plant roots. As a consequence, only ion contents in the shoot are presented (Fig. 1, D and E). Both AM and NM plants watered with low-K+ solution displayed a significant lower shoot K+ content than those grown under K+-sufficient conditions (Fig. 1D). A significant reduction of K+ content was observed in NM plants in comparison with AM plants exclusively under K+ limitation (Fig. 1D).
Together, these results indicate that the mycorrhizal status helped M. truncatula deal efficiently with long-term K+ deprivation. Interestingly, the sodium content was twice higher in the shoot of NM plants than AM plants and whatever the external K+ availability, suggesting a buffering role of the fungus in preventing the accumulation of sodium cations into the plant (Fig. 1E).
AM Association Reverses Most Transcriptional Responses to K+ Deprivation and Reduces the Production of Reactive Oxygen Species in M. truncatula Roots
RNA-seq was performed on roots from AM and NM plants grown under a low- or high-K+ regime. RNAs from three to four biological replicates composed of four plants each were sequenced, and the reliability of our data set was analyzed using principal component analyses and multidimensional scaling (Supplemental Figs. S1 and S2). K-means clustering also was performed to identify transcripts following the same expression profile along each condition: AM or NM plants at high (AM+K or NM+K) or low (AM−K or NM−K) K+ levels. This analysis revealed groups of genes with condition-specific responses: for example, two clusters of genes that were mainly up-regulated (clusters 1 and 21) and two that were mainly down-regulated (clusters 8 and 20) in the AM−K condition (Supplemental Fig. S3; Supplemental Table S4).
The first step was to identify transcripts specifically regulated at a low-K+ level in NM plants. Taking the NM+K condition as a control, DESeq analysis allowed the identification of 353 and 317 up- and down-regulated genes, respectively, in M. truncatula roots under K+ deprivation (Fig. 2A; Supplemental Table S1). Gene Ontology (GO) enrichment analysis identified a very significant overrepresentation of metal-binding transcripts, particularly among up-regulated transcripts (Fig. 2B). Of these proteins annotated as metal binding, one-third encoded peroxidases (six and two were up- and down-regulated, respectively) or cytochrome P450 enzymes (two and 15 were up- and down-regulated, respectively). These proteins are known to play an important role in oxidative stress responses in plants, in particular under K+ deprivation (Apel and Hirt, 2004; Hernandez et al., 2012; Wang et al., 2012). Overall, 62 transcripts encoding proteins involved in oxidative stress were differentially regulated under K+ deprivation (Supplemental Table S2), and both up- and down-regulated transcript groups were significantly enriched for oxidoreductases (Fig. 2B). Also, two-thirds of the differentially expressed hydrolases were down-regulated (Fig. 2B; Supplemental Table S3). Overall, our data revealed that many enzymes were differentially regulated in M. truncatula roots under K+ deprivation, suggesting that posttranslational regulations play an important role in the long-term adaptation to K+ deprivation.
Transcriptional profiling of M. truncatula roots under K+ deprivation. A, Venn diagram showing the number of M. truncatula root genes up-regulated (yellow) and down-regulated (green) in response to K+ deprivation, based on q < 0.05 and a 1.5-fold change cutoff threshold. B, Significantly enriched GO molecular function terms for up-regulated (yellow) and down-regulated (green) genes analyzed with the AgriGO Web-based tool (http://bioinfo.cau.edu.cn/agriGO/analysis.php).
To evaluate how AM symbiosis affects the transcriptional response of M. truncatula to K+ deprivation, the expression of transcripts regulated in NM−K plants was investigated in AM−K plants. Investigating the expression of genes in AM plants at high and low K+ levels revealed that the expression of two-thirds of them (387 transcripts) was altered in AM plants under K+ deprivation (Fig. 3). As mentioned above, many of these transcripts code for proteins likely involved in oxidative stress responses. Thus, the production of reactive oxygen species (ROS) in M. truncatula roots colonized or not by R. irregularis at high or low K+ was determined using carboxy-2,7-difluorodihydrofluorescein diacetate (carboxy-H2DFFDA; Fig. 4). In NM plants, K+ deprivation significantly enhanced the production of ROS in roots. However, this effect was reduced significantly and visible in only a few cells when plants were mycorrhized (Fig. 4).
Clustering and heat map of the expression values for differentially expressed genes (DEGs) in NM plants under K+ deprivation. The heat map displays DEGs identified previously in NM plants at low K+ (Fig. 2; Supplemental Table S1). An important part of the transcriptional responses of M. truncatula roots under K+ deprivation was altered during mycorrhizal association. White and red areas indicate higher and lower expression values, respectively.
Impact of AM symbiosis and K+ deprivation on the production of ROS in M. truncatula. A, ROS fluorescence images are shown in AM and NM plants at high (+K) and low (−K) K+ levels. Green pseudocolor indicates ROS production. Bars = 250 μm. B, Quantified data from three representative images taken from six plants per condition. ROS were visualized by staining the roots of plants inoculated and watered with high- or low-K+ solutions for 6 weeks with 20 μm carboxy-H2DFFDA (DFFDA). **, P < 0.01.
Altogether, these observations revealed that AM symbiosis compensates the transcriptional response of M. truncatula roots and prevents the production of ROS in roots under K+-limiting conditions.
Specific Transcriptional Responses Are Triggered by AM Association under K+ Deprivation in M. truncatula Roots
Using NM+K plants as a control, DEGs regulated in AM+K, AM−K, and NM−K plants were identified to highlight specific and shared expression patterns in each condition (Supplemental Fig. S4; Supplemental Table S5). Although 230 transcripts were similarly and significantly regulated in both AM+K and AM−K plants, revealing molecular players that were involved in AM symbiosis independently of extracellular K+ availability, including various transport systems, 511 transcripts were regulated specifically in plant roots under AM−K conditions (Supplemental Fig. S5; Supplemental Table S5). A GO enrichment analysis revealed an overrepresentation of genes involved in redox homeostasis, plant cell wall formation, and responses to oxidative stress (Fig. 5A). Interestingly, five transcripts encoding for transport systems were up-regulated in AM plants under K+ deprivation: an ortholog of the plasma membrane K+/H+ exchanger AtCHX20 from Arabidopsis (Medtr7g099800.1), a putative major facilitator superfamily (MFS) transporter (Medtr2g081930.1), a mitochondrial phosphate transporter (Medtr6g033280.1), and vacuolar malate (Medtr4g133230.1) and ion (Medtr4g094332.1) transporters (Fig. 5B; Supplemental Fig. S6). Among the 30 most DEGs in AM−K plants (Table I), a clade A type 2C protein phosphatase (Medtr5g009370.1) was expressed and could be involved in the posttranslational regulation of K+ transport systems. These results suggest that AM symbiosis activates specific mechanisms to tolerate long-term K+ deprivation, including transport mechanisms to likely facilitate the acquisition of K+.
GO enrichment and heat map of the expression values for DEGs in transport activity in AM plants under K+ deprivation. A, Significantly enriched GO terms in AM−K plants specifically in comparison with NM+K plants. Dashed arrows indicate that all the GO levels are not represented. B, The heat map displays DEGs in transport activity specifically up-regulated in AM−K plants. White and red areas indicate higher and lower expression values, respectively.
Identification of Putative Regulators for the Long-Term Adaptation of M. truncatula to K+ Deprivation in AM and NM Roots
We sought to identify regulatory networks associated with the AM and K+ DEG sets described previously (Supplemental Fig. S4) by comparing them with our global M. truncatula regulatory network (Marx et al., 2016). This network was inferred using the MERLIN algorithm (Roy et al., 2013). We defined five condition-specific gene sets, each corresponding to a bar in Supplemental Figure S4 and Supplemental Table S5. To assess the connectivity patterns, we overlaid the gene sets on the regulatory network (Supplemental Methods S1). Figure 6 shows regulatory network connections involving the DEGs and their prioritized network neighbors, summarized into three major categories: genes differentially expressed in AM roots regardless of K+ level (Fig. 6Ai), genes specifically differentially expressed in the AM−K condition (Fig. 6Bi), and genes differentially expressed under K+ deprivation regardless of AM status (Fig. 6Ci). Genes within each gene set were statistically enriched for direct interactions compared with randomly chosen genes, controlling for the degree (Z ≥ 7.09; Supplemental Table S6). Within each condition, genes with the same direction of expression change (up- or down-regulated) represented distinct network components (Fig. 6, i images) and had significantly shorter network distances than pairs of genes with opposite expression patterns (Mann-Whitney rank-sum test P < 0.001; Supplemental Fig. S7). This result indicates that genes that were induced represent distinct components of a pathway compared with genes that were repressed. Distinct regulatory systems are likely involved in controlling the up- and down-regulated components of each DEG set.
Network analysis of DEGs. Shown are inferred regulatory interactions spanning DEGs and genes prioritized by the diffusion kernel. In each image, the same network is shown (i.e. same edges and nodes), but they are colored differently to highlight DEG sets or top prioritized genes. A, Network for AM plants, regardless of K+ level (bars I and II in Supplemental Fig. S4). i, Up- and down-regulated DEGs as dark and light green nodes. ii, Genes prioritized by the diffusion kernel are colored from white to red; size and red intensity are proportional to the prioritization score. B, Network for the AM−K condition (bar III in Supplemental Fig. S4). i, Up- and down-regulated DEGs as dark and light orange. ii, Prioritized genes. C, Network for K+ deprivation (−K), regardless of mycorrhizal status (bars IV and V in Supplemental Fig. S4). i, Up- and down-regulated DEGs as dark and light purple. ii, Prioritized genes. Regulators (transcription factors and signaling proteins) are indicated with triangles, whereas all other genes are represented as circles. Small white nodes are genes (DEG or prioritized) in other conditions. Black edges are between genes relevant to the condition of interest, and gray edges are between genes from other conditions. For details, see Supplemental Methods S1, Supplemental Figure S8, and Supplemental Tables S1 and S2.
Two complementary network-based approaches were used to identify important regulatory genes associated with responses to AM and low-K+ conditions: diffusion kernel-based prioritization and statistical enrichment of targets of regulators in DEG sets. While the diffusion kernel approach can prioritize any gene in the network, the gene set enrichment analysis focuses on regulatory proteins such as transcription factors and signaling proteins. The diffusion kernel-based approach identified several candidate genes from AM roots putatively involved in the response to K+ deprivation (Fig. 6Bii; Table II). For example, among the top 25 predicted genes for the AM−K DEG set were nine kinases, including a calcineurin B-like (CBL)-interacting kinase (Medtr2g105010.1) and an ethylene response factor (Medtr4g078710.1; Table II). Full prioritization results for all gene sets are available in Supplemental Table S7.
Shown are the top 25 prioritized genes for the AM−K DEG set. Prioritization (Rank and Score) columns give the rank and score from the diffusion kernel analysis. Target Enrichment indicates whether a prioritized gene in this list also has target enrichment in a DEG set; nonregulators are designated NA in this column.
Statistical enrichment of regulatory targets in the DEG sets identified 56 putative regulators (Supplemental Table S7; Supplemental Fig. S8; Supplemental Methods S1) that included both new and known proteins such as MYB, BHLH, and AP2/ERF transcription factors as well as chromatin-modifying enzymes. For example, targets from a CBL-interacting kinase (Medtr2g105010.1; nine targets), a Cys desulfurase (Medtr8g093560.1; seven targets), and a GATA transcription factor (Medtr3g109760.1; 24 targets) were enriched for up-regulated genes in response to K+ deprivation in both AM−K and NM−K conditions, suggesting their involvement in the tolerance to low-K+ stress. Similarly, targets from a dehydration-responsive element-binding protein (Medtr1g019110.1; six targets), a C2H2-type zinc finger protein (Medtr3g102980.1; three targets), and a leucine-rich repeat (LRR) receptor-like kinase (Medtr8g014970.1; 10 targets) were specifically enriched for up-regulated genes in the AM−K condition. Finally, 29 and 17 targets from a cyclin-dependent kinase (Medtr8g092290.1) and a Tyr kinase (Medtr7g116650.1), respectively, were enriched for up-regulated genes in AM roots regardless of extracellular K+ level, suggesting a major role of these proteins in the regulation of AM symbiosis.
Interestingly, the regulators identified by both methods were consistent with each other. The 56 regulators identified by the enrichment analysis received diffusion kernel percentile ranks of 95 or greater compared with all ranked transcription factors and signaling proteins in the regulatory network. Each method also was able to identify some unique prioritized genes that the other would miss or rank lower. Overall, prioritization and regulator enrichment approaches revealed promising candidates for further functional analyses of the regulatory network controlling responses to changes in AM fungi and K+ availability in M. truncatula.
DISCUSSION
K+ Deprivation Affects Growth, K+ Acquisition, and Gene Expression in M. truncatula Roots
K+ is an essential macronutrient for plants. Its limitation contributes to slowing down many vital processes. Plants elicit rapid responses to strong K+ deprivation by inducing plasma membrane hyperpolarization and the up-regulation of high-affinity transporters (Schachtman and Shin, 2007). Prolonged K+ deficiency has been reported previously to limit root development (Armengaud et al., 2004; Shin and Schachtman, 2004; Kellermeier et al., 2013), which was observed in our experiments (Fig. 1). Transcriptional analysis of plants facing short-term and extreme K+ deprivation revealed an increase in the expression of K+ transporters and channels from the HAK, HKT, and AKT families (Armengaud et al., 2004; Ma et al., 2012; Shankar et al., 2013; Ruan et al., 2015; Zeng et al., 2015). In our study, we analyzed the impact of long-term K+ deprivation in the model legume M. truncatula. In our case, no K+ transport system was differentially regulated in response to long-term K+ deprivation. Because the goal of our study was also to evaluate how a prolonged low-K+ treatment impacts AM symbiosis, and reciprocally how the AM status alters the plant responses to K+ deprivation, the stress condition was 0.05 mm external K+ over 6 weeks. This treatment was obviously very different from that in previous studies in which short and extreme K+ stresses were applied. Therefore, our results suggest that the up-regulation of HAK, HKT, and AKT transport systems observed in other plant species is probably a rapid and transient response to an intense stress.
It is well established that K+ deficiency contributes to the accumulation of ROS in plant roots, leading to a general oxidative stress (Shin and Schachtman, 2004; Kim et al., 2010; Hernandez et al., 2012). Long-term K+ deprivation in M. truncatula resulted in the differential expression of many transcripts related to oxidative stress, including peroxidases, cytochrome P450 enzymes, and reductases. Most of the peroxidases were up-regulated under K+ deficiency in M. truncatula roots. Interestingly, Kim et al. (2010) hypothesized that the up-regulation of peroxidases could be the origin of ROS production observed in Arabidopsis roots. The accumulation of ROS production in M. truncatula roots detected under long-term K+ deprivation supports this view and validates the in silico prediction we made with our RNA-seq approach. Overall, many enzymes were differentially regulated in stressed plants, suggesting that posttranscriptional mechanisms may be involved mainly in the tolerance of K+ deprivation (Schachtman and Shin, 2007).
AM Symbiosis Modulates the Plant Responses to K+ Deprivation
The role of endomycorrhizal and ectomycorrhizal associations in plant K+ nutrition has been debated for decades. Although some studies described an improvement of K+ acquisition in plants colonized by mycorrhizal fungi, others presented apparently contradictory results (for review, see Garcia and Zimmermann, 2014). The transport of nutrients from the external medium to the host plant through mycorrhizal structures has been investigated using radioactive isotopes (He et al., 2009; Smith et al., 2011). Because K+ isotopes are not stable enough to track K+ movement from the soil to the mycorrhizal plants, some studies have used rubidium as an analog tracer (Rygiewicz and Bledsoe, 1984; Hawkes and Casper, 2002). Further experiments will be needed in M. truncatula using rubidium to evaluate the actual transport of K+ from the soil to the mycorrhizal plants. However, in a previous publication on ectomycorrhizae, we demonstrated an improvement of maritime pine (Pinus pinaster) K+ nutrition when colonized by an ectomycorrhizal fungus only under K+ deprivation (Garcia et al., 2014). In this study, the shoot K+ content of M. truncatula also was significantly higher in AM plants compared with NM plants but only at low K+ levels and after 6 weeks of coculture. This effect should be investigated in other plant species. We can assume that under K+-sufficient conditions, host plants can acquire external K+ by themselves. However, when K+ ions in solution have been limiting for a relatively long period of time, AM fungi may help their hosts to acquire K+ from the soil. Our analysis revealed that M. truncatula plants colonized by R. irregularis displayed increases of biomass and K+ ion content as well as specific transcriptional responses to the low-K+ regime. These results confirmed that AM symbiosis was able to provide specific adaptation mechanisms to the host plant to tolerate long-term K+ deprivation.
ROS have been detected in legume roots colonized by both AM fungi and nitrogen-fixing bacteria, particularly in colonized cells (Salzer et al., 1999; Fester and Hause, 2005; Puppo et al., 2013). Moreover, some NADPH oxidase-encoding genes were recently found to be up-regulated in colonized cortical cells of M. truncatula roots, suggesting their role in arbuscule development (Belmondo et al., 2016a, 2016b). Our results revealed that the two-thirds of M. truncatula root transcriptional responses to low K+ were absent in mycorrhizal plants, including many enzymes putatively involved in ROS production. The ROS accumulation observed in NM roots was correlated with the external K+ availability but not with the AM symbiosis, suggesting specific ROS production in response to K+ deprivation. This accumulation was reduced significantly in AM plants at low K+, indicating that AM symbiosis helped M. truncatula to cope with long-term K+ deprivation.
Also, NM plants accumulated much more sodium in shoots than AM plants. Sodium can be toxic to the cells at a high level. By preventing its accumulation, AM symbiosis may protect M. truncatula from salt stress by buffering the uptake of sodium. Although the role of AM associations on salt tolerance was suggested in other plants (e.g. basil [Ocimum basilicum; Zuccarini and Okurowska, 2008], olive [Olea europaea; Porras-Soriano et al., 2009], and maize [Estrada et al., 2013]), further experiments will be needed to unravel the molecular mechanisms of this adaptation in M. truncatula.
The uptake of nutrients from inner symbiotic structures requires the specific expression and regulation of plasma membrane transport systems in colonized cortical cells (Casieri et al., 2013; Garcia et al., 2016). Some plant transporters were described to transport phosphorus or nitrogen from AM fungi to plant cells in arbuscules (Harrison et al., 2002; Breuillin-Sessoms et al., 2015). One study in the legume Lotus japonicus reported a 44-fold up-regulation of a K+ transporter in AM roots compared with NM roots (Guether et al., 2009). However, in transcriptome studies of mycorrhizal roots from M. truncatula, no K+ transporter was ever found up-regulated in AM roots under the standard K+ regime (Gomez et al., 2009; Gaude et al., 2012; this study). In our study, we found that some genes encoding putative transporters were up-regulated in mycorrhizal plants under K+ deprivation, particularly a putative K+/H+ exchanger (CHX; Medtr7g099800.1). A close homolog of this CHX protein in Arabidopsis was found expressed in guard cells (Padmanaban et al., 2007). Other members of the CHX family were regulated in starved K+ roots, such as AtCHX17 (Cellier et al., 2004). As a consequence, further analyses will be required in M. truncatula to determine if and how this CHX can transport K+ during AM symbiosis under long-term K+ deprivation.
Network Analysis Identified Putative Regulators Controlling the Tolerance of Mycorrhizal Roots to K+ Deprivation
The MERLIN algorithm (Roy et al., 2013) was used to infer an M. truncatula regulatory network from publicly available gene expression data representing a wide array of conditions and tissues. Such a network analysis allows us to identify regulatory connections between DEGs but also to predict which specific regulators (signaling proteins and transcription factors) may control the expression of genes of interest. In our study, this approach identified differentially regulated genes and prioritized putative regulators that could be involved in controlling the tolerance of AM plants to K+ deprivation. Among the top 30 up-regulated genes in AM−K roots, a clade A type 2C protein phosphatase (Medtr5g009370.1) was identified. This gene encodes an ortholog of AtPP2CA in Arabidopsis involved in the regulation of the weak-rectifying K+ channel AKT2 (Chérel et al., 2002) as well as of the uptake and efflux Shaker K+ channels AKT1 and GORK, respectively (Lee et al., 2007; Lan et al., 2011; Lefoulon et al., 2016). Both prioritization and target enrichment analyses suggested the importance of a CBL-interacting kinase (Medtr2g105010.1) in AM−K plant roots. CBL-interacting kinases are major regulators of nutrient transporters, abscisic acid responses, and K+ homeostasis (D’Angelo et al., 2006; Ho et al., 2009; Hashimoto et al., 2012; Liu et al., 2013; Ragel et al., 2015). This observation suggests a role of Medtr2g105010.1 in the posttranslational regulation of transport systems in M. truncatula roots colonized by AM fungi to promote the symbiotic acquisition of K+. The targets of regulators involved in plant adaptation to oxidative stress also were predicted in AM−K roots and included a Cys desulfurase (Medtr8g093560.1) and a dehydration-responsive element-binding protein (Medtr1g019110). Dehydration-responsive element-binding proteins are involved in salt and hydric stresses, which are consequences of K+ deprivation in plants, and the overexpression of one of them in Arabidopsis (AtDREBP2C) improved its tolerance to global oxidative stress (Hwang et al., 2012).
Altogether, our study revealed that AM symbiosis allows host plants to cope with long-term K+ deprivation. This tolerance mechanism probably involves the regulation of a gene network to protect host plants against oxidative stress but also mechanisms to facilitate K+ uptake from the soil through symbiotic structures. Further studies will be needed to validate these predictions and to fully decipher the role of AM associations in plant K+ nutrition.
MATERIALS AND METHODS
Plant Growth Conditions, Fungal Inoculation, and Long-Term K+ Deficiency Treatment
Medicago truncatula ‘Jemalong A17’ seeds were acid scarified and surface sterilized, plated on 1% (w/v) agar supplemented with 1 µg mL−1 GA3, vernalized at 4°C for 4 d, and allowed to germinate overnight at room temperature. The germinated seedlings were placed for 2 weeks on modified Fahräeus medium as described previously (Catoira et al., 2000). The plants were transferred to pots filled with Turface and inoculated with 400 spores of Rhizophagus irregularis. Plants were watered regularly with standard (+K; 3.75 mm K+) Long Ashton solution [3.75 mm KNO3, 2 mm Ca(NO3)2·4H2O, 7.5 μm NaH2PO4·H2O, 1 mm MgSO4·7H2O, 0.05 mm NaCl, 5 μm MnSO4, 0.5 μm CuSO4, 1 μm ZnSO4, 16.5 μm H3BO3, and 0.1 μm Na2MoO4] or low-K+ (−K; 0.05 mm K+) solution [0.05 mm KNO3, 3.85 mm Ca(NO3)2·4H2O, 7.5 μm NaH2PO4·H2O, 1 mm MgSO4·7H2O, 0.05 mm NaCl, 5 μm MnSO4, 0.5 μm CuSO4, 1 μm ZnSO4, 16.5 μm H3BO3, and 0.1 μm Na2MoO4]. The plants were harvested after 6 weeks of treatment.
Mycorrhizal Quantification, Dry Weight Determination, and Ion Content in Plants
Shoots and roots of 6-week-old M. truncatula plants growing in K+-sufficient or K+-deficient conditions were harvested separately. Mycorrhizal colonization rates were assessed using the gridline intersection method on roots stained previously with Sheaffer ink (McGonigle et al., 1990). Fresh weight was determined, and roots and shoots were dried at 60°C for 1 week for dry weight determination. The dried samples were sent to the Soil & Forage Analysis Laboratory (University of Wisconsin-Madison) to determine the K+, sodium, phosphorus, and calcium contents in plant tissues by ICP-OES.
RNA Isolation and Sequencing
Frozen root tissues were ground in liquid nitrogen, and total RNA was extracted from the resulting powder using the PureLink RNA Mini Kit (Thermo Fisher) and treated with TURBO DNase (Thermo Fisher). RNA quality was assessed using 2100 BioAnalyzer technology (Agilent Technologies). Six micrograms of total RNA was used to construct poly(A) selection libraries using the Illumina TruSeq RNA Sample Preparation kit; 100-nucleotide single-end reads were obtained on an Illumina HiSeq2000 platform using TruSeq version 3 (University of Wisconsin Biotechnology Center DNA Sequencing Facility).
RNA-seq Analysis, Annotations, Cutoff, and Analysis
RNA-seq data were processed to obtain DEGs using the TopHat2-SE, Cufflinks2, Cuffmerge2, and Cuffdiff2 applications from the iPlant Collaborative platform (http://www.iplantcollaborative.org). The DESeq Bioconductor package version 3.3 for R (Anders and Huber, 2010) was used for the analysis of DEGs. Cuffdiff2 sequencing data were analyzed using the cummeRbund package to generate principal component and multidimensional scaling analyses on biological and technical replicates. K-means clustering also was performed using the cummeRbund package for R with K = 25. For further analyses, DEGs were called by DESeq using a 1.5-fold change cutoff threshold and false discovery rate-adjusted q < 0.05. DEGs were annotated with the Mt4.0v1 version of the M. truncatula ‘Jemalong A17’ genome (http://jcvi.org/medicago/index.php). GO enrichment analysis was conducted with the AgriGO Web-based tool (http://bioinfo.cau.edu.cn/agriGO/analysis.php) using the Mt4.0v1 version of the M. truncatula genome. The Venn diagram was obtained using the Bioinformatics & Evolutionary Genomics Web tool (http://bioinformatics.psb.ugent.be/webtools/Venn/). Finally, clustering and heat maps were generated on R with scripts based on the heatmap.2 function as available in the gplots Bioconductor package. The data discussed in this article have been deposited in the National Center for Biotechnology Information’s Gene Expression Omnibus (Edgar et al., 2002) and are accessible through the Gene Expression Omnibus Series accession number GSE94266 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE94266).
ROS Detection and Measurement
To detect ROS, the root samples of AM and NM 6-week-old plants grown under K+-sufficient or K+-deficient conditions were incubated with 20 μm carboxy-H2DFFDA (Thermo Fisher) for 20 min in high- or low-K+ solution, respectively. The roots were washed twice with high- or low-K+ solution before microscopy visualization. All fluorescence images were captured using a Leica DMi8 microscope and a Leica DFC365 FX camera. ROS fluorescence was quantified and converted into pseudocolor images using the ImageJ software program. Background noise was subtracted from the fluorescence intensity value for quantification.
Regulatory Network Analysis and Visualization
The base network visualized in all parts of Figure 6 was selected from a global inferred M. truncatula regulatory network (Marx et al., 2016) by extracting any interaction between a DEG and top n computationally prioritized genes for each K+ and AM condition, where n is the size of the experimentally identified DEG set. For visual simplicity, connected components involving only one edge were omitted. Network images were created using Cytoscape (Shannon et al., 2003). For full details of the regulatory network analyses, including prioritization and regulator gene set enrichment, see Supplemental Methods S1.
Accession Numbers
Sequence data from this article can be found in the GenBank/EMBL data libraries under the accession numbers GSE94266.
Supplemental Data
The following supplemental materials are available.
Supplemental Figure S1. Principal component analysis of the changes in transcript abundance in AM and NM M. truncatula roots under high- and low-K+ regimes.
Supplemental Figure S2. Multidimensional scaling plot showing the two-dimensional distribution of the samples.
Supplemental Figure S3. K-means cluster analysis of DEGs of M. truncatula roots under high- and low-K+ regimes.
Supplemental Figure S4. Impact of AM symbiosis on the transcriptional profiling of M. truncatula roots under K+ deprivation.
Supplemental Figure S5. Representative clustering and heat map of expression values for DEG in transport activity specifically regulated by AM symbiosis independently of K+ availability.
Supplemental Figure S6. Phylogenetic analysis of the cation/H+ exchanger Medtr7g099800.1.
Supplemental Figure S7. Network connectivity statistics suggest differential regulation for induced and repressed gene sets.
Supplemental Figure S8. Predicted regulators of mycorrhizal and low-K+ gene sets represent relevant pathways and propose candidates for future study.
Supplemental Table S1. List of up- and down-regulated transcripts in NM M. truncatula roots under K+ deprivation.
Supplemental Table S2. List of differentially expressed oxidative stress-related transcripts in NM M. truncatula roots under K+ deprivation.
Supplemental Table S3. List of differentially expressed hydrolase transcripts in NM M. truncatula roots under K+ deprivation.
Supplemental Table S4. List of transcripts categorized in the 25 clusters of the K-means cluster analysis (Supplemental Fig. S4).
Supplemental Table S5. List of transcripts differentially regulated in M. truncatula roots in comparison with the NM−K condition.
Supplemental Table S6. Edge density analysis results for input gene sets.
Supplemental Table S7. Genes ranked by diffusion kernel and enrichment-based prioritization analyses.
Supplemental Methods S1. Regulatory network-based interpretation and prioritization of DEGs.
Acknowledgments
We thank Dr. Sabine Zimmermann for constructive comments on the article and the University of Wisconsin Biotechnology Center DNA Sequencing Facility for providing RNA-seq facilities and services.
Footnotes
The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantphysiol.org) is: Jean-Michel Ané (jeanmichel.ane{at}wisc.edu).
K.G. and J.-M.A. conceived the research plans; J.-M.A. supervised the experiments; K.G. designed and performed most of the experiments; K.G. and D.C. analyzed the data; K.G. and D.C. wrote the article with contributions of all the authors; J.-M.A. and S.R. supervised and complemented the writing.
↵1 This work was supported by the National Science Foundation (grant no. NSF-IOS 1331098 to J.-M.A. and CAREER award grant no. NSF-DBI 1350677), by a Sloan Foundation research fellowship, by University of Wisconsin-Madison startup funds to S.R., and by the Environmental Protection Agency (grant no. 83573701).
↵[OPEN] Articles can be viewed without a subscription.
Glossary
- AM
- arbuscular mycorrhizal
- RNA-seq
- RNA sequencing
- NM
- nonmycorrhized
- ICP-OES
- inductively coupled plasma optical emission spectrometry
- ROS
- reactive oxygen species
- GO
- Gene Ontology
- DEG
- differentially expressed gene
- Received December 22, 2016.
- Accepted January 31, 2017.
- Published February 3, 2017.