Plasticity of the Arabidopsis Root System under Nutrient Deficiencies

Plant roots show a particularly high variation in their morphological response to different nutrient deficiencies. Although such changes often determine the nutrient efficiency or stress tolerance of plants, it is surprising that a comprehensive and comparative analysis of root morphological responses to different nutrient deficiencies has not yet been conducted. Since one reason for this is an inherent difficulty in obtaining nutrient-deficient conditions in agar culture, we first identified conditions appropriate for producing nutrient-deficient plants on agar plates. Based on a careful selection of agar specifically for each nutrient being considered, we grew Arabidopsis (Arabidopsis thaliana) plants at four levels of deficiency for 12 nutrients and quantified seven root traits. In combination with measurements of biomass and elemental concentrations, we observed that the nutritional status and type of nutrient determined the extent and type of changes in root system architecture (RSA). The independent regulation of individual root traits further pointed to a differential sensitivity of root tissues to nutrient limitations. To capture the variation in RSA under different nutrient supplies, we used principal component analysis and developed a root plasticity chart representing the overall modulations in RSA under a given treatment. This systematic comparison of RSA responses to nutrient deficiencies provides a comprehensive view of the overall changes in root plasticity induced by the deficiency of single nutrients and provides a solid basis for the identification of nutrient-sensitive steps in the root developmental program.Plant survival and performance are highly dependent on the plant’s ability to efficiently explore the soil in the search for water and minerals. Thus, root growth and architecture are extremely relevant for the plant’s adaptation to the growth medium, as they determine the soil volume that a plant is able to explore at a given time. Root system architecture (RSA) represents the spatial arrangement of roots of different ages and orders (Lynch, 1995; Osmont et al., 2007) and is determined by genetic factors and the integration of environmental cues (Malamy, 2005). The genetic component determines the fundamental morphology and blueprint of a plant’s root system, whereas environmental cues shape root architecture by modifying the intrinsic genetic program. The existence of this additional level of regulation allows plants to display a high level of root plasticity, which reflects the shape, three-dimensional distribution, branching pattern, and age of the primary and postembryonically generated roots (Pacheco-Villalobos and Hardtke, 2012). The dynamic modulation of RSA is based on the intrinsic developmental nature of the different components of the root system. In fact, the primary root (PR) is established during embryogenesis, while the lateral roots (LRs) that originate from the PR develop postembryonically (Osmont et al., 2007; Péret et al., 2009). These highly dynamic changes in the overall RSA throughout time finally determine root plasticity and allow plants to efficiently adapt to environmental constraints.Nutrient availability can exert a profound impact on RSA by altering the number, length, angle, and diameter of roots and root hairs (for review, see Forde and Lorenzo, 2001; López-Bucio et al., 2003; Malamy, 2005; Osmont et al., 2007). In fact, plants can respond to the heterogenous availability of resources by allocating roots where the most favorable conditions are found (Zhang and Forde, 1998; Linkohr et al., 2002; Remans et al., 2006; Lima et al., 2010; Giehl et al., 2012). When grown under limited phosphorus (P) availability, roots exhibit a shallower architecture that results from the inhibition of PR elongation and the concomitant increase in LR formation (Williamson et al., 2001; López-Bucio et al., 2002; Sanchez-Calderon et al., 2005). Such an architectural rearrangement of the root is thought to improve the plant’s ability to forage P from the usually P-enriched topsoil horizon (Lynch and Brown, 2001; Rubio et al., 2003; Zhu et al., 2005). In contrast to low P, reduced nitrogen (N) availability stimulates PR and particularly LR elongation but not LR initiation (Linkohr et al., 2002; López-Bucio et al., 2003). However, it is noteworthy that under severe N shortage, LR formation is almost completely absent (Krouk et al., 2010), suggesting that plants require a certain level of N to sustain an active foraging strategy. These examples indicate that the availability of different nutrients can evoke distinct effects on RSA that depend upon which nutrient is supplied and the concentration of the supplied nutrient.Unfortunately, for the majority of the nutrients, a more detailed analysis of the architectural modifications under deficient conditions is still missing. In fact, most studies describe the effect of nutrient deficiencies on root growth and development only in terms of root biomass or total root length (Hermans and Verbruggen, 2005; Hermans et al., 2006; Richard-Molard et al., 2008; Jung et al., 2009; Cailliatte et al., 2010). Thus, important features of the root system are not comprehensible from these rather basic measurements. The characterization of RSA in more detail appears justified due to the positive correlations found between single root characteristics and plant yield, especially when the supply of water or mineral resources was limited (Landi et al., 2002; Tuberosa et al., 2002; Manschadi et al., 2006; Kirkegaard et al., 2007; Steele et al., 2007). Although a large number of studies have been conducted on the root development of grasses (Hochholdinger and Tuberosa, 2009; Iyer-Pascuzzi et al., 2010; Pacheco-Villalobos and Hardtke, 2012), our understanding of the molecular players involved in the regulation of root growth and development has benefited most from studies of the reference plant Arabidopsis (Arabidopsis thaliana) grown under controlled conditions to minimize variability. However, imposing consistent nutrient deficiencies presents an experimental challenge as long as plants are grown on agar medium, which is the method of choice to preserve the spatial arrangement of the root system and access a larger number of root traits.A major drawback of agar and agarose media is their inherent nutrient load, such that traces of nutrient contamination must often be made unavailable to plants, for example by adding chelating agents to lower the free activities of micronutrients (Bell et al., 1991; Yang et al., 1994; Rengel, 1999). Additionally, in many cases, symptoms of deficiency are only observed in mutants impaired in the uptake of the nutrient in question (Tomatsu et al., 2007; Mills et al., 2008; Assunção et al., 2010). In general, gelling agents may contribute considerable amounts of nutrients (Debergh, 1983; Scholten and Pierik, 1998), hampering the occurrence of deficiency for specific nutrients (Jain et al., 2009). Thus, it becomes crucial to select the most suitable gelling agent when particular nutrient deficiencies are to be obtained. This is particularly relevant as strategies depending upon the use of gelling media are being developed to overcome the bottleneck that often limits RSA traits from being characterized in high-throughput phenotyping studies (Iyer-Pascuzzi et al., 2010; Clark et al., 2011).In our approach to compare RSA under different nutrient deficiencies in Arabidopsis plants grown on solid medium, we first identified the most appropriate conditions for producing nutrient-deficient plants on agar plates. Once identified, these conditions allowed us to characterize the effects of 12 deficiencies at four intensity levels on the RSA by measuring seven root traits. These measurements, in combination with biomass and elemental concentrations, allowed us to determine the nutrient-specific effects on particular parameters of the RSA and thus to describe the root plasticity of Arabidopsis and analyze the underlying traits under different nutrient deficiencies.     

Natural Variation of Arabidopsis Root Architecture Reveals Complementing Adaptive Strategies to Potassium Starvation

Root architecture is a highly plastic and environmentally responsive trait that enables plants to counteract nutrient scarcities with different foraging strategies. In potassium (K) deficiency (low K), seedlings of the Arabidopsis (Arabidopsis thaliana) reference accession Columbia (Col-0) show a strong reduction of lateral root elongation. To date, it is not clear whether this is a direct consequence of the lack of K as an osmoticum or a triggered response to maintain the growth of other organs under limiting conditions. In this study, we made use of natural variation within Arabidopsis to look for novel root architectural responses to low K. A comprehensive set of 14 differentially responding root parameters were quantified in K-starved and K-replete plants. We identified a phenotypic gradient that links two extreme strategies of morphological adaptation to low K arising from a major tradeoff between main root (MR) and lateral root elongation. Accessions adopting strategy I (e.g. Col-0) maintained MR growth but compromised lateral root elongation, whereas strategy II genotypes (e.g. Catania-1) arrested MR elongation in favor of lateral branching. K resupply and histochemical staining resolved the temporal and spatial patterns of these responses. Quantitative trait locus analysis of K-dependent root architectures within a Col-0 × Catania-1 recombinant inbred line population identified several loci each of which determined a particular subset of root architectural parameters. Our results indicate the existence of genomic hubs in the coordinated control of root growth in stress conditions and provide resources to facilitate the identification of the underlying genes.The ability of plants to actively respond to nutrient scarcity with changes in root architecture is a fascinating phenomenon. Advances in root research and breeding efforts that focus on the enhancement of root traits have been recognized as principal goals to ensure those high yields necessary to feed an ever-growing human population (Hammer et al., 2009; Den Herder et al., 2010). Indeed, understanding the adaptations of root systems to environmental factors has been pointed out as a key issue in modern agriculture (Den Herder et al., 2010).Potassium (K) is the quantitatively most important cation for plant growth, as it serves as the major osmoticum for cell expansion (Leigh and Wyn Jones, 1984; Amtmann et al., 2006). Moreover, K is essential for many cellular and tissue processes, such as enzymatic activity, transport of minerals and metabolites, and regulation of stomatal aperture (Amtmann et al., 2006). Even in fertilized fields, rapid K uptake by plants can lead to K shortage in the root environment, especially early in the growth season. Root adaptations to K deficiency (low K) take place at the physiological (Armengaud et al., 2004; Shin and Schachtman, 2004; Alemán et al., 2011), metabolic (Armengaud et al., 2009a), and morphological levels. In a classic study, Drew (1975) showed an increase in overall lateral root (LR) growth of barley seedlings, even when K was supplied only to parts of the root system. Conversely, a typical response of Arabidopsis (Arabidopsis thaliana) Columbia (Col-0) seedlings to low K is the drastic reduction of LR elongation (Armengaud et al., 2004; Shin and Schachtman, 2004). Conflicting data have been published on the effect of low K on main root (MR) growth in the same species, ranging from no effect (Shin and Schachtman, 2004) to impaired MR elongation (Jung et al., 2009; Kim et al., 2010). Some components involved in K starvation responses have been identified, such as jasmonates (Armengaud et al., 2004, 2010), reactive oxygen species (Shin and Schachtman, 2004), and ethylene (Jung et al., 2009). However, the molecular identity of a root K sensor acting at the base of the signaling cascade is so far unknown.Genetic variation within species is a useful resource to dissect the genetic components determining phenotypes (Koornneef et al., 2004; Trontin et al., 2011; Weigel, 2012). Natural variation within Arabidopsis has been the basis for many studies on plant morphology, physiology, and development as well as stress response (Alonso-Blanco et al., 2009; Weigel, 2012). Natural variation of root traits such as primary root length (Mouchel et al., 2004; Loudet et al., 2005; Sergeeva et al., 2006), LR length (Loudet et al., 2005), and total root size (Fitz Gerald et al., 2006) have pinpointed genomic regions underlying the phenotypic variation via mapping of quantitative trait loci (QTLs) as a first step toward the identification of novel regulatory genes (Mouchel et al., 2004). This strategy has also been applied to environmental responses, such as growth responses to phosphate starvation (Reymond et al., 2006; Svistoonoff et al., 2007). However, despite their importance for plant growth and their strong effect on overall root architecture, root responses to K deficiency have not been genetically dissected.Here, we show that Arabidopsis accessions follow different strategies to adapt to K starvation. We present the quantification of a comprehensive set of root architectural parameters of Arabidopsis grown in K-sufficient and K-deficient media and the identification of genetic loci, each of which determines the response of a distinct subset of root architectural parameters to K starvation.     

Salt Stress Reduces Root Meristem Size by Nitric Oxide-Mediated Modulation of Auxin Accumulation and Signaling in Arabidopsis

The development of the plant root system is highly plastic, which allows the plant to adapt to various environmental stresses. Salt stress inhibits root elongation by reducing the size of the root meristem. However, the mechanism underlying this process remains unclear. In this study, we explored whether and how auxin and nitric oxide (NO) are involved in salt-mediated inhibition of root meristem growth in Arabidopsis (Arabidopsis thaliana) using physiological, pharmacological, and genetic approaches. We found that salt stress significantly reduced root meristem size by down-regulating the expression of PINFORMED (PIN) genes, thereby reducing auxin levels. In addition, salt stress promoted AUXIN RESISTANT3 (AXR3)/INDOLE-3-ACETIC ACID17 (IAA17) stabilization, which repressed auxin signaling during this process. Furthermore, salt stress stimulated NO accumulation, whereas blocking NO production with the inhibitor N^ω-nitro-l-arginine-methylester compromised the salt-mediated reduction of root meristem size, PIN down-regulation, and stabilization of AXR3/IAA17, indicating that NO is involved in salt-mediated inhibition of root meristem growth. Taken together, these findings suggest that salt stress inhibits root meristem growth by repressing PIN expression (thereby reducing auxin levels) and stabilizing IAA17 (thereby repressing auxin signaling) via increasing NO levels.Due to agricultural practices and climate change, soil salinity has become a serious factor limiting the productivity and quality of agricultural crops (Zhu, 2007). Worldwide, high salinity in the soil damages approximately 20% of total irrigated lands and takes 1.5 million ha out of production each year (Munns and Tester, 2008). In general, high salinity affects plant growth and development by reducing plant water potential, altering nutrient uptake, and increasing the accumulation of toxic ions (Hasegawa et al., 2000; Munns, 2002; Zhang and Shi, 2013). Together, these effects severely reduce plant growth and survival.Because the root is the first organ to sense high salinity, salt stress plays a direct, important role in modulating root system architecture (Wang et al., 2009). For instance, salt stress negatively regulates root hair formation and gravitropism (Sun et al., 2008; Wang et al., 2008). The role of salt in lateral root formation depends on the NaCl concentration. While high NaCl levels inhibit lateral root formation, lower NaCl levels stimulate lateral root formation in an auxin-dependent manner (Zolla et al., 2010; Ji et al., 2013). The root meristem plays an essential role in sustaining root growth (Perilli et al., 2012). Salt stress inhibits primary root elongation by suppressing root meristem activity (West et al., 2004). However, how this inhibition occurs remains largely unclear.Plant hormones are important intermediary signaling compounds that function downstream of environmental stimuli. Among plant hormones, indole-3-acetic acid (IAA) is thought to play a fundamental role in root system architecture by regulating cell division, expansion, and differentiation. In Arabidopsis (Arabidopsis thaliana) root tips, a distal auxin maximum is formed and maintained by polar auxin transport (PAT), which determines the orientation and extent of cell division in the root meristem as well as root pattern formation (Sabatini et al., 1999). PINFORMED (PIN) proteins, which are components of the auxin efflux machinery, regulate primary root elongation and root meristem size (Blilou et al., 2005; Dello Ioio et al., 2008; Yuan et al., 2013, 2014). The auxin signal transduction pathway is activated by direct binding of auxin to its receptor protein, TRANSPORT INHIBITOR RESPONSE1 (TIR1)/AUXIN SIGNALING F-BOX (AFB), promoting the degradation of Aux/IAA proteins, releasing auxin response factors (ARFs), and activating the expression of auxin-responsive genes (Gray et al., 2001; Dharmasiri et al., 2005a; Kepinski and Leyser, 2005). Aux/IAA proteins are short-lived, nuclear-localized proteins that play key roles in auxin signal activation and root growth modulation (Rouse et al., 1998). Other hormones and stresses often regulate auxin signaling by affecting Aux/IAA protein stability (Lim and Kunkel, 2004; Nemhauser et al., 2004; Wang et al., 2007; Kushwah and Laxmi, 2014).Nitric oxide (NO) is a signaling molecule with diverse biological functions in plants (He et al., 2004; Fernández-Marcos et al., 2011; Shi et al., 2012), including important roles in the regulation of root growth and development. NO functions downstream of auxin during the adventitious rooting process in cucumber (Cucumis sativus; Pagnussat et al., 2002). Exogenous auxin-induced NO biosynthesis is associated with nitrate reductase activity during lateral root formation, and NO is necessary for auxin-induced lateral root and root hair development (Pagnussat et al., 2002; Lombardo et al., 2006). Pharmacological and genetic analyses in Arabidopsis indicate that NO suppresses primary root growth and root meristem activity (Fernández-Marcos et al., 2011). Additionally, both exogenous application of the NO donor sodium nitroprusside (SNP) and overaccumulation of NO in the mutant chlorophyll a/b binding protein underexpressed1 (cue1)/nitric oxide overproducer1 (nox1) result in reduced PIN1 expression and auxin accumulation in root tips. The auxin receptors protein TIR1 is S-nitrosylated by NO, suggesting that this protein is a direct target of NO in the regulation of root development (Terrile et al., 2012).Because NO is a free radical, NO levels are dynamically regulated by endogenous and environmental cues. Many phytohormones, including abscisic acid, auxin, cytokinin, salicylic acid, jasmonic acid, and ethylene, induce NO biosynthesis (Zottini et al., 2007; Kolbert et al., 2008; Tun et al., 2008; García et al., 2011). In addition, many abiotic and biotic stresses or stimuli, such as cold, heat, salt, drought, heavy metals, and pathogens/elicitors, also stimulate NO biosynthesis (Zhao et al., 2009; Mandal et al., 2012). Salt stress stimulates NO and ONOO^– accumulation in roots (Corpas et al., 2009), but the contribution of NO to root meristem growth under salinity stress has yet to be examined in detail.In this study, we found that salt stress significantly down-regulated the expression of PIN genes and promoted AUXIN RESISTANT3 (AXR3)/IAA17 stabilization. Furthermore, salt stress stimulated NO accumulation, and pharmacological inhibition of NO biosynthesis compromised the salt-mediated reduction in root meristem size. Our results support a model in which salt stress reduces root meristem size by increasing NO accumulation, which represses PIN expression and stabilizes IAA17, thereby reducing auxin levels and repressing auxin signaling.     

Spore Density Determines Infection Strategy by the Plant Pathogenic Fungus Plectosphaerella cucumerina

Necrotrophic and biotrophic pathogens are resisted by different plant defenses. While necrotrophic pathogens are sensitive to jasmonic acid (JA)-dependent resistance, biotrophic pathogens are resisted by salicylic acid (SA)- and reactive oxygen species (ROS)-dependent resistance. Although many pathogens switch from biotrophy to necrotrophy during infection, little is known about the signals triggering this transition. This study is based on the observation that the early colonization pattern and symptom development by the ascomycete pathogen Plectosphaerella cucumerina (P. cucumerina) vary between inoculation methods. Using the Arabidopsis (Arabidopsis thaliana) defense response as a proxy for infection strategy, we examined whether P. cucumerina alternates between hemibiotrophic and necrotrophic lifestyles, depending on initial spore density and distribution on the leaf surface. Untargeted metabolome analysis revealed profound differences in metabolic defense signatures upon different inoculation methods. Quantification of JA and SA, marker gene expression, and cell death confirmed that infection from high spore densities activates JA-dependent defenses with excessive cell death, while infection from low spore densities induces SA-dependent defenses with lower levels of cell death. Phenotyping of Arabidopsis mutants in JA, SA, and ROS signaling confirmed that P. cucumerina is differentially resisted by JA- and SA/ROS-dependent defenses, depending on initial spore density and distribution on the leaf. Furthermore, in situ staining for early callose deposition at the infection sites revealed that necrotrophy by P. cucumerina is associated with elevated host defense. We conclude that P. cucumerina adapts to early-acting plant defenses by switching from a hemibiotrophic to a necrotrophic infection program, thereby gaining an advantage of immunity-related cell death in the host.Plant pathogens are often classified as necrotrophic or biotrophic, depending on their infection strategy (Glazebrook, 2005; Nishimura and Dangl, 2010). Necrotrophic pathogens kill living host cells and use the decayed plant tissue as a substrate to colonize the plant, whereas biotrophic pathogens parasitize living plant cells by employing effector molecules that suppress the host immune system (Pel and Pieterse, 2013). Despite this binary classification, the majority of pathogenic microbes employ a hemibiotrophic infection strategy, which is characterized by an initial biotrophic phase followed by a necrotrophic infection strategy at later stages of infection (Perfect and Green, 2001). The pathogenic fungi Magnaporthe grisea, Sclerotinia sclerotiorum, and Mycosphaerella graminicola, the oomycete Phytophthora infestans, and the bacterial pathogen Pseudomonas syringae are examples of hemibiotrophic plant pathogens (Perfect and Green, 2001; Koeck et al., 2011; van Kan et al., 2014; Kabbage et al., 2015).Despite considerable progress in our understanding of plant resistance to necrotrophic and biotrophic pathogens (Glazebrook, 2005; Mengiste, 2012; Lai and Mengiste, 2013), recent debate highlights the dynamic and complex interplay between plant-pathogenic microbes and their hosts, which is raising concerns about the use of infection strategies as a static tool to classify plant pathogens. For instance, the fungal genus Botrytis is often labeled as an archetypal necrotroph, even though there is evidence that it can behave as an endophytic fungus with a biotrophic lifestyle (van Kan et al., 2014). The rice blast fungus Magnaporthe oryzae, which is often classified as a hemibiotrophic leaf pathogen (Perfect and Green, 2001; Koeck et al., 2011), can adopt a purely biotrophic lifestyle when infecting root tissues (Marcel et al., 2010). It remains unclear which signals are responsible for the switch from biotrophy to necrotrophy and whether these signals rely solely on the physiological state of the pathogen, or whether host-derived signals play a role as well (Kabbage et al., 2015).The plant hormones salicylic acid (SA) and jasmonic acid (JA) play a central role in the activation of plant defenses (Glazebrook, 2005; Pieterse et al., 2009, 2012). The first evidence that biotrophic and necrotrophic pathogens are resisted by different immune responses came from Thomma et al. (1998), who demonstrated that Arabidopsis (Arabidopsis thaliana) genotypes impaired in SA signaling show enhanced susceptibility to the biotrophic pathogen Hyaloperonospora arabidopsidis (formerly known as Peronospora parastitica), while JA-insensitive genotypes were more susceptible to the necrotrophic fungus Alternaria brassicicola. In subsequent years, the differential effectiveness of SA- and JA-dependent defense mechanisms has been confirmed in different plant-pathogen interactions, while additional plant hormones, such as ethylene, abscisic acid (ABA), auxins, and cytokinins, have emerged as regulators of SA- and JA-dependent defenses (Bari and Jones, 2009; Cao et al., 2011; Pieterse et al., 2012). Moreover, SA- and JA-dependent defense pathways have been shown to act antagonistically on each other, which allows plants to prioritize an appropriate defense response to attack by biotrophic pathogens, necrotrophic pathogens, or herbivores (Koornneef and Pieterse, 2008; Pieterse et al., 2009; Verhage et al., 2010).In addition to plant hormones, reactive oxygen species (ROS) play an important regulatory role in plant defenses (Torres et al., 2006; Lehmann et al., 2015). Within minutes after the perception of pathogen-associated molecular patterns, NADPH oxidases and apoplastic peroxidases generate early ROS bursts (Torres et al., 2002; Daudi et al., 2012; O’Brien et al., 2012), which activate downstream defense signaling cascades (Apel and Hirt, 2004; Torres et al., 2006; Miller et al., 2009; Mittler et al., 2011; Lehmann et al., 2015). ROS play an important regulatory role in the deposition of callose (Luna et al., 2011; Pastor et al., 2013) and can also stimulate SA-dependent defenses (Chaouch et al., 2010; Yun and Chen, 2011; Wang et al., 2014; Mammarella et al., 2015). However, the spread of SA-induced apoptosis during hyperstimulation of the plant immune system is contained by the ROS-generating NADPH oxidase RBOHD (Torres et al., 2005), presumably to allow for the sufficient generation of SA-dependent defense signals from living cells that are adjacent to apoptotic cells. Nitric oxide (NO) plays an additional role in the regulation of SA/ROS-dependent defense (Trapet et al., 2015). This gaseous molecule can stimulate ROS production and cell death in the absence of SA while preventing excessive ROS production at high cellular SA levels via S-nitrosylation of RBOHD (Yun et al., 2011). Recently, it was shown that pathogen-induced accumulation of NO and ROS promotes the production of azelaic acid, a lipid derivative that primes distal plants for SA-dependent defenses (Wang et al., 2014). Hence, NO, ROS, and SA are intertwined in a complex regulatory network to mount local and systemic resistance against biotrophic pathogens. Interestingly, pathogens with a necrotrophic lifestyle can benefit from ROS/SA-dependent defenses and associated cell death (Govrin and Levine, 2000). For instance, Kabbage et al. (2013) demonstrated that S. sclerotiorum utilizes oxalic acid to repress oxidative defense signaling during initial biotrophic colonization, but it stimulates apoptosis at later stages to advance necrotrophic colonization. Moreover, SA-induced repression of JA-dependent resistance not only benefits necrotrophic pathogens but also hemibiotrophic pathogens after having switched from biotrophy to necrotrophy (Glazebrook, 2005; Pieterse et al., 2009, 2012).Plectosphaerella cucumerina ((P. cucumerina, anamorph Plectosporum tabacinum) anamorph Plectosporum tabacinum) is a filamentous ascomycete fungus that can survive saprophytically in soil by decomposing plant material (Palm et al., 1995). The fungus can cause sudden death and blight disease in a variety of crops (Chen et al., 1999; Harrington et al., 2000). Because P. cucumerina can infect Arabidopsis leaves, the P. cucumerina-Arabidopsis interaction has emerged as a popular model system in which to study plant defense reactions to necrotrophic fungi (Berrocal-Lobo et al., 2002; Ton and Mauch-Mani, 2004; Carlucci et al., 2012; Ramos et al., 2013). Various studies have shown that Arabidopsis deploys a wide range of inducible defense strategies against P. cucumerina, including JA-, SA-, ABA-, and auxin-dependent defenses, glucosinolates (Tierens et al., 2001; Sánchez-Vallet et al., 2010; Gamir et al., 2014; Pastor et al., 2014), callose deposition (García-Andrade et al., 2011; Gamir et al., 2012, 2014; Sánchez-Vallet et al., 2012), and ROS (Tierens et al., 2002; Sánchez-Vallet et al., 2010; Barna et al., 2012; Gamir et al., 2012, 2014; Pastor et al., 2014). Recent metabolomics studies have revealed large-scale metabolic changes in P. cucumerina-infected Arabidopsis, presumably to mobilize chemical defenses (Sánchez-Vallet et al., 2010; Gamir et al., 2014; Pastor et al., 2014). Furthermore, various chemical agents have been reported to induce resistance against P. cucumerina. These chemicals include β-amino-butyric acid, which primes callose deposition and SA-dependent defenses, benzothiadiazole (BTH or Bion; Görlach et al., 1996; Ton and Mauch-Mani, 2004), which activates SA-related defenses (Lawton et al., 1996; Ton and Mauch-Mani, 2004; Gamir et al., 2014; Luna et al., 2014), JA (Ton and Mauch-Mani, 2004), and ABA, which primes ROS and callose deposition (Ton and Mauch-Mani, 2004; Pastor et al., 2013). However, among all these studies, there is increasing controversy about the exact signaling pathways and defense responses contributing to plant resistance against P. cucumerina. While it is clear that JA and ethylene contribute to basal resistance against the fungus, the exact roles of SA, ABA, and ROS in P. cucumerina resistance vary between studies (Thomma et al., 1998; Ton and Mauch-Mani, 2004; Sánchez-Vallet et al., 2012; Gamir et al., 2014).This study is based on the observation that the disease phenotype during P. cucumerina infection differs according to the inoculation method used. We provide evidence that the fungus follows a hemibiotrophic infection strategy when infecting from relatively low spore densities on the leaf surface. By contrast, when challenged by localized host defense to relatively high spore densities, the fungus switches to a necrotrophic infection program. Our study has uncovered a novel strategy by which plant-pathogenic fungi can take advantage of the early immune response in the host plant.     

Focus Issue on Roots: Duplicate and Conquer: Multiple Homologs of PHOSPHORUS-STARVATION TOLERANCE1 Enhance Phosphorus Acquisition and Sorghum Performance on Low-Phosphorus Soils

Low soil phosphorus (P) availability is a major constraint for crop production in tropical regions. The rice (Oryza sativa) protein kinase, PHOSPHORUS-STARVATION TOLERANCE1 (OsPSTOL1), was previously shown to enhance P acquisition and grain yield in rice under P deficiency. We investigated the role of homologs of OsPSTOL1 in sorghum (Sorghum bicolor) performance under low P. Association mapping was undertaken in two sorghum association panels phenotyped for P uptake, root system morphology and architecture in hydroponics and grain yield and biomass accumulation under low-P conditions, in Brazil and/or in Mali. Root length and root surface area were positively correlated with grain yield under low P in the soil, emphasizing the importance of P acquisition efficiency in sorghum adaptation to low-P availability. SbPSTOL1 alleles reducing root diameter were associated with enhanced P uptake under low P in hydroponics, whereas Sb03g006765 and Sb03g0031680 alleles increasing root surface area also increased grain yield in a low-P soil. SbPSTOL1 genes colocalized with quantitative trait loci for traits underlying root morphology and dry weight accumulation under low P via linkage mapping. Consistent allelic effects for enhanced sorghum performance under low P between association panels, including enhanced grain yield under low P in the soil in Brazil, point toward a relatively stable role for Sb03g006765 across genetic backgrounds and environmental conditions. This study indicates that multiple SbPSTOL1 genes have a more general role in the root system, not only enhancing root morphology traits but also changing root system architecture, which leads to grain yield gain under low-P availability in the soil.Increasing food production is one of the major global challenges in dealing with continuously growing population and food consumption (Godfray et al., 2010). One of the major obstacles to improve crop production in tropical regions is phosphorus (P) deficiency caused by P fixation in the soil clays. P is one of the most important plant nutrients, contributing approximately 0.2% of a plant’s dry weight, and is a component of key organic molecules such as nucleic acids, phospholipids, and ATP (Schachtman et al., 1998). On tropical soils, even when the total amount of soil P is high, its bioavailability is low due to P fixation by aluminum and iron oxides in clay minerals (Marschner, 1995) and immobilization into organic forms (Schachtman et al., 1998). Approximately half of the world’s agricultural lands occurs on low-P soils (Lynch, 2011); hence, crop adaptation to P insufficiency should be a major breeding target to enable sustainable agricultural production worldwide. In addition, because phosphate rock fertilizer is a nonrenewable resource that is being depleted by agricultural demands, increasing fertilizer prices negatively impact agriculture, particularly for small-holder farmers in developing countries in the tropics and subtropics (Cordell et al., 2009; Sattari et al., 2012). In sorghum (Sorghum bicolor), breeding strategies for low-P adaptation have been developed based on multienvironment trials in West Africa, indicating the importance of undertaking selection in low-P soil conditions (Leiser et al., 2012a, 2012b). Therefore, developing crops with greater ability to grow and maintain satisfactory yields on soils with reduced P availability is expected to substantially improve food security worldwide.Tolerance to P deficiency in plants can be achieved by mechanisms underlying both P acquisition and P internal utilization efficiency (Parentoni and Souza Junior, 2008). One of the major mechanisms that plants have evolved to overcome low-P availability is to maximize the ability of the roots to acquire and absorb P from the soil. Plants can mobilize P through the exudation of organic acids, acid phosphatases, and ribonucleases, resulting in enhanced P availability and uptake (Hinsinger, 2001; Ryan et al., 2001; Dakora and Phillips, 2002; Hammond and White, 2008; Ma et al., 2009; Pang et al., 2009). Another strategy to cope with low-P availability is to increase the soil volume accessed by root systems by forming mycorrhizal symbioses (Li et al., 2012; Smith and Smith, 2012; Rai et al., 2013). Due to low-P mobility on tropical soils, changes in root architecture and morphology enhance P uptake by facilitating soil exploration (Williamson et al., 2001; Ho et al., 2005; Walk et al., 2006; Svistoonoff et al., 2007; Lynch, 2011; Ingram et al., 2012; Niu et al., 2013). Root structural changes leading to higher P uptake include increased root hair growth (Yan et al., 2004; Haling et al., 2013; Lan et al., 2013) and length and enhancing lateral root over primary root growth (Williamson et al., 2001; Wang et al., 2013). In addition, increased root surface area is achieved by a combination of reduced root diameter and enhanced elongation of relatively thinner roots (Fitter et al., 2002). There is both intraspecific and interspecific genetic variation for P deficiency tolerance in crop species (Lynch and Brown, 2001, 2012; Mudge et al., 2002; Paszkowski et al., 2002; Rausch and Bucher, 2002; Huang et al., 2011; Zhang et al., 2011; Leiser et al., 2014a) that can be explored to develop P-efficient cultivars.In rice (Oryza sativa), Phosphorus uptake1 (Pup1), a major quantitative trait locus (QTL) for P deficiency tolerance donated by an aus-type Indian variety, Kasalath, was mapped to the long arm of chromosome 12 (Ni et al., 1998; Wissuwa et al., 1998, 2002; Heuer et al., 2009). Near-isogenic lines bearing the Kasalath allele at Pup1 showed 3-fold higher P uptake and grain yield in low-P trials compared with the recurrent parent, cv Nipponbare, which is intolerant to P starvation (Wissuwa and Ae, 2001). Following high-resolution mapping of Pup1, comparative sequence analyses of homologous bacterial artificial chromosomes showed that a Kasalath genomic fragment contained several genes not present in cv Nipponbare, highlighting an approximately 90-kb deletion in the cv Nipponbare reference genome that encompassed the Pup1 locus (Heuer et al., 2009). Within this insertion/deletion, OsPupK46-2, a gene encoding a Ser/Thr kinase of the Receptor-like Protein Kinase LRK10L-2 subfamily, was found to enhance grain yield and P uptake in rice lines overexpressing this gene, indicating that this protein kinase underlies the Pup1 locus (Gamuyao et al., 2012). OsPupK46-2, which is now designated PHOSPHORUS-STARVATION TOLERANCE1 (OsPSTOL1), was found to be up-regulated in the root tissues of tolerant near-isogenic lines under P-deficient conditions and was shown to increase P uptake by a physiological mechanism based on the enhancement of early root growth and development. Furthermore, lines overexpressing OsPupK46-2 showed an approximately 30% grain yield increase in comparison with the null lines, suggesting that PSTOL1 has potential for molecular breeding applications to improve crop performance under low-P conditions. Consistent with the proposed physiological mechanism underlying OsPSTOL1, the superior performance of the transgenic lines was related to enhanced root dry weight, root length, and root surface area (Gamuyao et al., 2012).Sorghum is the world’s fifth most important cereal crop and is a staple food for more than half a billion people. It is widely adapted to harsh environmental conditions, and more specifically, to arid and semiarid regions of the world (Doumbia et al., 1993, 1998). It has been estimated that rice diverged from its most recent common ancestor with sorghum and maize (Zea mays) approximately 50 million years ago (Kellogg, 1998; Paterson et al., 2000, 2004; Paterson, 2008). About 60% of the genes in the sorghum genome are located in syntenic regions to rice (Paterson et al., 2009), emphasizing the potential for using comparative genomics for cross-species identification of genes underlying abiotic stress tolerance in the grass family. Here, we applied association analysis to specifically study the role of sorghum homologs of rice OsPSTOL1 in tolerance to P starvation in sorghum. Single-nucleotide polymorphisms (SNPs) within PSTOL1 homologs in sorghum, collectively designated SbPSTOL1, were significantly associated with grain yield under low-P conditions and also root morphology and root system architecture (RSA) traits phenotyped from hydroponically grown plants. Under low P, SbPSTOL1 genes increased biomass accumulation and P content in the African landrace panel and grain yield in the sorghum association panel phenotyped in a low-P Brazilian soil. This suggests a stable effect across environments and sorghum genotypes that potentially can be used for molecular breeding applications. QTL mapping with a large sorghum recombinant inbred line population was used to validate the association results, indicating that SbPSTOL1 homologs colocalize with QTLs related to root morphology and performance under low P. Our results indicate that SbPSTOL1 homologs have the ability to enhance P uptake and sorghum performance in low-P soils by a mechanism related not only to early root growth enhancement, as was the case for rice OsPSTOL1, but also by modulating RSA.     

Focus on Ethylene: Ethylene Biosynthesis Is Promoted by Very-Long-Chain Fatty Acids during Lysigenous Aerenchyma Formation in Rice Roots

In rice (Oryza sativa) roots, lysigenous aerenchyma, which is created by programmed cell death and lysis of cortical cells, is constitutively formed under aerobic conditions, and its formation is further induced under oxygen-deficient conditions. Ethylene is involved in the induction of aerenchyma formation. reduced culm number1 (rcn1) is a rice mutant in which the gene encoding the ATP-binding cassette transporter RCN1/OsABCG5 is defective. Here, we report that the induction of aerenchyma formation was reduced in roots of rcn1 grown in stagnant deoxygenated nutrient solution (i.e. under stagnant conditions, which mimic oxygen-deficient conditions in waterlogged soils). 1-Aminocyclopropane-1-carboxylic acid synthase (ACS) is a key enzyme in ethylene biosynthesis. Stagnant conditions hardly induced the expression of ACS1 in rcn1 roots, resulting in low ethylene production in the roots. Accumulation of saturated very-long-chain fatty acids (VLCFAs) of 24, 26, and 28 carbons was reduced in rcn1 roots. Exogenously supplied VLCFA (26 carbons) increased the expression level of ACS1 and induced aerenchyma formation in rcn1 roots. Moreover, in rice lines in which the gene encoding a fatty acid elongase, CUT1-LIKE (CUT1L; a homolog of the gene encoding Arabidopsis CUT1, which is required for cuticular wax production), was silenced, both ACS1 expression and aerenchyma formation were reduced. Interestingly, the expression of ACS1, CUT1L, and RCN1/OsABCG5 was induced predominantly in the outer part of roots under stagnant conditions. These results suggest that, in rice under oxygen-deficient conditions, VLCFAs increase ethylene production by promoting 1-aminocyclopropane-1-carboxylic acid biosynthesis in the outer part of roots, which, in turn, induces aerenchyma formation in the root cortex.Aerenchyma formation is a morphological adaptation of plants to complete submergence and waterlogging of the soil, and facilitates internal gas diffusion (Armstrong, 1979; Jackson and Armstrong, 1999; Colmer, 2003; Voesenek et al., 2006; Bailey-Serres and Voesenek, 2008; Licausi and Perata, 2009; Sauter, 2013; Voesenek and Bailey-Serres, 2015). To adapt to waterlogging in soil, rice (Oryza sativa) develops lysigenous aerenchyma in shoots (Matsukura et al., 2000; Colmer and Pedersen, 2008; Steffens et al., 2011) and roots (Jackson et al., 1985b; Justin and Armstrong, 1991; Kawai et al., 1998), which is formed by programmed cell death and subsequent lysis of some cortical cells (Jackson and Armstrong, 1999; Evans, 2004; Yamauchi et al., 2013). In rice roots, lysigenous aerenchyma is constitutively formed under aerobic conditions (Jackson et al., 1985b), and its formation is further induced under oxygen-deficient conditions (Colmer et al., 2006; Shiono et al., 2011). The former and latter are designated constitutive and inducible lysigenous aerenchyma formation, respectively (Colmer and Voesenek, 2009). The gaseous plant hormone ethylene regulates adaptive growth responses of plants to submergence (Voesenek and Blom, 1989; Voesenek et al., 1993; Visser et al., 1996a,b; Lorbiecke and Sauter, 1999; Hattori et al., 2009; Steffens and Sauter, 2009; van Veen et al., 2013). Ethylene also induces lysigenous aerenchyma formation in roots of some gramineous plants (Drew et al., 2000; Shiono et al., 2008). The treatment of roots with ethylene or its precursor (1-aminocyclopropane-1-carboxylic acid [ACC]) stimulates aerenchyma formation in rice (Justin and Armstrong, 1991; Colmer et al., 2006; Yukiyoshi and Karahara, 2014), maize (Zea mays; Drew et al., 1981; Jackson et al., 1985a; Takahashi et al., 2015), and wheat (Triticum aestivum; Yamauchi et al., 2014a,b). Moreover, treatment of roots with inhibitors of ethylene action or ethylene biosynthesis effectively blocks aerenchyma formation under hypoxic conditions in maize (Drew et al., 1981; Konings, 1982; Jackson et al., 1985a; Rajhi et al., 2011).Ethylene biosynthesis is accomplished by two main successive enzymatic reactions: conversion of S-adenosyl-Met to ACC by 1-aminocyclopropane-1-carboxylic acid synthase (ACS), and conversion of ACC to ethylene by 1-aminocyclopropane-1-carboxylic acid oxidase (ACO; Yang and Hoffman, 1984). The activities of both enzymes are enhanced during aerenchyma formation under hypoxic conditions in maize root (He et al., 1996). Since the ACC content in roots of maize is increased by oxygen deficiency and is strongly correlated with ethylene production (Atwell et al., 1988), ACC biosynthesis is essential for ethylene production during aerenchyma formation in roots. In fact, exogenously supplied ACC induced ethylene production in roots of maize (Drew et al., 1979; Konings, 1982; Atwell et al., 1988) and wheat (Yamauchi et al., 2014b), even under aerobic conditions. Ethylene production in plants is inversely related to oxygen concentration (Yang and Hoffman, 1984). Under anoxic conditions, the oxidation of ACC to ethylene by ACO, which requires oxygen, is almost completely repressed (Yip et al., 1988; Tonutti and Ramina, 1991). Indeed, anoxic conditions stimulate neither ethylene production nor aerenchyma formation in maize adventitious roots (Drew et al., 1979). Therefore, it is unlikely that the root tissues forming inducible aerenchyma are anoxic, and that the ACO-mediated step is repressed. Moreover, aerenchyma is constitutively formed in rice roots even under aerobic conditions (Jackson et al., 1985b), and thus, after the onset of waterlogging, oxygen can be immediately supplied to the apical regions of roots through the constitutively formed aerenchyma.Very-long-chain fatty acids (VLCFAs; ≥20 carbons) are major constituents of sphingolipids, cuticular waxes, and suberin in plants (Franke and Schreiber, 2007; Kunst and Samuels, 2009). In addition to their structural functions, VLCFAs directly or indirectly participate in several physiological processes (Zheng et al., 2005; Reina-Pinto et al., 2009; Roudier et al., 2010; Ito et al., 2011; Nobusawa et al., 2013; Tsuda et al., 2013), including the regulation of ethylene biosynthesis (Qin et al., 2007). During fiber cell elongation in cotton ovules, ethylene biosynthesis is enhanced by treatment with saturated VLCFAs, especially 24-carbon fatty acids, and is suppressed by an inhibitor of VLCFA biosynthesis (Qin et al., 2007). The first rate-limiting step in VLCFA biosynthesis is condensation of acyl-CoA with malonyl-CoA by β-ketoacyl-CoA synthase (KCS; Joubès et al., 2008). KCS enzymes are thought to determine the substrate and tissue specificities of fatty acid elongation (Joubès et al., 2008). The Arabidopsis (Arabidopsis thaliana) genome has 21 KCS genes (Joubès et al., 2008). In the Arabidopsis cut1 mutant, which has a defect in the gene encoding CUT1 that is required for cuticular wax production (i.e. one of the KCS genes), the expression of AtACO genes and growth of root cells were reduced when compared with the wild type (Qin et al., 2007). Furthermore, expression of the AtACO genes was rescued by exogenously supplied saturated VLCFAs (Qin et al., 2007). These observations imply that VLCFAs or their derivatives work as regulatory factors for gene expression during some physiological processes in plants.reduced culm number1 (rcn1) was first identified as a rice mutant with a low tillering rate in a paddy field (Takamure and Kinoshita, 1985; Yasuno et al., 2007). The rcn1 (rcn1-2) mutant has a single nucleotide substitution in the gene encoding a member of the ATP-binding cassette (ABC) transporter subfamily G, RCN1/OsABCG5, causing an Ala-684Pro substitution (Yasuno et al., 2009). The mutation results in several mutant phenotypes, although the substrates of RCN1/OsABCG5 have not been determined (Ureshi et al., 2012; Funabiki et al., 2013; Matsuda et al., 2014). We previously found that the rcn1 mutant has abnormal root morphology, such as shorter root length and brownish appearance of roots, under stagnant (deoxygenated) conditions (which mimics oxygen-deficient conditions in waterlogged soils). We also found that the rcn1 mutant accumulates less of the major suberin monomers originating from VLCFAs in the outer part of adventitious roots, and this results in a reduction of a functional apoplastic barrier in the root hypodermis (Shiono et al., 2014a).The objective of this study was to elucidate the molecular basis of inducible aerenchyma formation. To this end, we examined lysigenous aerenchyma formation and ACC, ethylene, and VLCFA accumulation and their biosyntheses in rcn1 roots. Based on the results of these studies, we propose that VLCFAs are involved in inducible aerenchyma formation through the enhancement of ethylene biosynthesis in rice roots.     

The Role of a Potassium Transporter OsHAK5 in Potassium Acquisition and Transport from Roots to Shoots in Rice at Low Potassium Supply Levels

In plants, K transporter (KT)/high affinity K transporter (HAK)/K uptake permease (KUP) is the largest potassium (K) transporter family; however, few of the members have had their physiological functions characterized in planta. Here, we studied OsHAK5 of the KT/HAK/KUP family in rice (Oryza sativa). We determined its cellular and tissue localization and analyzed its functions in rice using both OsHAK5 knockout mutants and overexpression lines in three genetic backgrounds. A β-glucuronidase reporter driven by the OsHAK5 native promoter indicated OsHAK5 expression in various tissue organs from root to seed, abundantly in root epidermis and stele, the vascular tissues, and mesophyll cells. Net K influx rate in roots and K transport from roots to aerial parts were severely impaired by OsHAK5 knockout but increased by OsHAK5 overexpression in 0.1 and 0.3 mm K external solution. The contribution of OsHAK5 to K mobilization within the rice plant was confirmed further by the change of K concentration in the xylem sap and K distribution in the transgenic lines when K was removed completely from the external solution. Overexpression of OsHAK5 increased the K-sodium concentration ratio in the shoots and salt stress tolerance (shoot growth), while knockout of OsHAK5 decreased the K-sodium concentration ratio in the shoots, resulting in sensitivity to salt stress. Taken together, these results demonstrate that OsHAK5 plays a major role in K acquisition by roots faced with low external K and in K upward transport from roots to shoots in K-deficient rice plants.Potassium (K) is one of the three most important macronutrients and the most abundant cation in plants. As a major osmoticum in the vacuole, K drives the generation of turgor pressure, enabling cell expansion. In the vascular tissue, K is an important participant in the generation of root pressure (for review, see Wegner, 2014 [including his new hypothesis]). In the phloem, K is critical for the transport of photoassimilates from source to sink (Marschner, 1996; Deeken et al., 2002; Gajdanowicz et al., 2011). In addition, enhancing K absorption and decreasing sodium (Na) accumulation is a major strategy of glycophytes in salt stress tolerance (Maathuis and Amtmann, 1999; Munns and Tester, 2008; Shabala and Cuin, 2008).Plants acquire K through K-permeable proteins at the root surface. Since available K concentration in the soil may vary by 100-fold, plants have developed multiple K uptake systems for adapting to this variability (Epstein et al., 1963; Grabov, 2007; Maathuis, 2009). In a classic K uptake experiment in barley (Hordeum vulgare), root K absorption has been described as a high-affinity and low-affinity biphasic transport process (Epstein et al., 1963). It is generally assumed that the low-affinity transport system (LATS) in the roots mediates K uptake in the millimolar range and that the activity of this system is insensitive to external K concentration (Maathuis and Sanders, 1997; Chérel et al., 2014). In contrast, the high-affinity transport system (HATS) was rapidly up-regulated when the supply of exogenous K was halted (Glass, 1976; Glass and Dunlop, 1978).The membrane transporters for K flux identified in plants are generally classified into three channels and three transporter families based on phylogenetic analysis (Mäser et al., 2001; Véry and Sentenac, 2003; Lebaudy et al., 2007; Alemán et al., 2011). For K uptake, it was predicted that, under most circumstances, K transporters function as HATS, while K-permeable channels mediate LATS (Maathuis and Sanders, 1997). However, a root-expressed K channel in Arabidopsis (Arabidopsis thaliana), Arabidopsis K Transporter1 (AKT1), mediates K absorption over a wide range of external K concentrations (Sentenac et al., 1992; Lagarde et al., 1996; Hirsch et al., 1998; Spalding et al., 1999), while evidence is accumulating that many K transporters, including members of the K transporter (KT)/high affinity K transporter (HAK)/K uptake permease (KUP) family, are low-affinity K transporters (Quintero and Blatt, 1997; Senn et al., 2001), implying that functions of plant K channels and transporters overlap at different K concentration ranges.Out of the three families of K transporters, cation proton antiporter (CPA), high affinity K/Na transporter (HKT), and KT/HAK/KUP, CPA was characterized as a K⁺(Na⁺)/H⁺ antiporter, HKT may cotransport Na and K or transport Na only (Rubio et al., 1995; Uozumi et al., 2000), while KT/HAK/KUP were predicted to be H⁺-coupled K⁺ symporters (Mäser et al., 2001; Lebaudy et al., 2007). KT/HAK/KUP were named by different researchers who first identified and cloned them (Quintero and Blatt, 1997; Santa-María et al., 1997). In plants, the KT/HAK/KUP family is the largest K transporter family, including 13 members in Arabidopsis and 27 members in the rice (Oryza sativa) genome (Rubio et al., 2000; Mäser et al., 2001; Bañuelos et al., 2002; Gupta et al., 2008). Sequence alignments show that genes of this family share relatively low homology to each other. The KT/HAK/KUP family was divided into four major clusters (Rubio et al., 2000; Gupta et al., 2008), and in cluster I and II, they were further separated into A and B groups. Genes of cluster I or II likely exist in all plants, cluster III is composed of genes from both Arabidopsis and rice, while cluster IV includes only four rice genes (Grabov, 2007; Gupta et al., 2008).The functions of KT/HAK/KUP were studied mostly in heterologous expression systems. Transporters of cluster I, such as AtHAK5, HvHAK1, OsHAK1, and OsHAK5, are localized in the plasma membrane (Kim et al., 1998; Bañuelos et al., 2002; Gierth et al., 2005) and exhibit high-affinity K uptake in the yeast Saccharomyces cerevisiae (Santa-María et al., 1997; Fu and Luan, 1998; Rubio et al., 2000) and in Escherichia coli (Horie et al., 2011). Transporters of cluster II, like AtKUP4 (TINY ROOT HAIRS1, TRH1), HvHAK2, OsHAK2, OsHAK7, and OsHAK10, could not complement the K uptake-deficient yeast (Saccharomyces cerevisiae) but were able to mediate K fluxes in a bacterial mutant; they might be tonoplast transporters (Senn et al., 2001; Bañuelos et al., 2002; Rodríguez-Navarro and Rubio, 2006). The function of transporters in clusters III and IV is even less known (Grabov, 2007).Existing data suggest that some KT/HAK/KUP transporters also may respond to salinity stress (Maathuis, 2009). The cluster I transporters of HvHAK1 mediate Na influx (Santa-María et al., 1997), while AtHAK5 expression is inhibited by Na (Rubio et al., 2000; Nieves-Cordones et al., 2010). Expression of OsHAK5 in tobacco (Nicotiana tabacum) BY2 cells enhanced the salt tolerance of these cells by accumulating more K without affecting their Na content (Horie et al., 2011).There are only scarce reports on the physiological function of KT/HAK/KUP in planta. In Arabidopsis, mutation of AtKUP2 (SHORT HYPOCOTYL3) resulted in a short hypocotyl, small leaves, and a short flowering stem (Elumalai et al., 2002), while a loss-of-function mutation of AtKUP4 (TRH1) resulted in short root hairs and a loss of gravity response in the root (Rigas et al., 2001; Desbrosses et al., 2003; Ahn et al., 2004). AtHAK5 is the only system currently known to mediate K uptake at concentrations below 0.01 mm (Rubio et al., 2010) and provides a cesium uptake pathway (Qi et al., 2008). AtHAK5 and AtAKT1 are the two major physiologically relevant molecular entities mediating K uptake into roots in the range between 0.01 and 0.05 mm (Pyo et al., 2010; Rubio et al., 2010). AtAKT1 may contribute to K uptake within the K concentrations that belong to the high-affinity system described by Epstein et al. (1963).Among all 27 members of the KT/HAK/KUP family in rice, OsHAK1, OsHAK5, OsHAK19, and OsHAK20 were grouped in cluster IB (Gupta et al., 2008). These four rice HAK members share 50.9% to 53.4% amino acid identity with AtHAK5. OsHAK1 was expressed in the whole plant, with maximum expression in roots, and was up-regulated by K deficiency; it mediated high-affinity K uptake in yeast (Bañuelos et al., 2002). In this study, we examined the tissue-specific localization and the physiological functions of OsHAK5 in response to variation in K supply and to salt stress in rice. By comparing K uptake and translocation in OsHAK5 knockout (KO) mutants and in OsHAK5-overexpressing lines with those in their respective wild-type lines supplied with different K concentrations, we found that OsHAK5 not only mediates high-affinity K acquisition but also participates in root-to-shoot K transport as well as in K-regulated salt tolerance.     

Multiscale Metabolic Modeling: Dynamic Flux Balance Analysis on a Whole-Plant Scale

Plant metabolism is characterized by a unique complexity on the cellular, tissue, and organ levels. On a whole-plant scale, changing source and sink relations accompanying plant development add another level of complexity to metabolism. With the aim of achieving a spatiotemporal resolution of source-sink interactions in crop plant metabolism, a multiscale metabolic modeling (MMM) approach was applied that integrates static organ-specific models with a whole-plant dynamic model. Allowing for a dynamic flux balance analysis on a whole-plant scale, the MMM approach was used to decipher the metabolic behavior of source and sink organs during the generative phase of the barley (Hordeum vulgare) plant. It reveals a sink-to-source shift of the barley stem caused by the senescence-related decrease in leaf source capacity, which is not sufficient to meet the nutrient requirements of sink organs such as the growing seed. The MMM platform represents a novel approach for the in silico analysis of metabolism on a whole-plant level, allowing for a systemic, spatiotemporally resolved understanding of metabolic processes involved in carbon partitioning, thus providing a novel tool for studying yield stability and crop improvement.Plants are of vital significance as a source of food (Grusak and DellaPenna, 1999; Rogalski and Carrer, 2011), feed (Lu et al., 2011), energy (Tilman et al., 2006; Parmar et al., 2011), and feedstocks for the chemical industry (Metzger and Bornscheuer, 2006; Kinghorn et al., 2011). Given the close connection between plant metabolism and the usability of plant products, there is a growing interest in understanding and predicting the behavior and regulation of plant metabolic processes. In order to increase crop quality and yield, there is a need for methods guiding the rational redesign of the plant metabolic network (Schwender, 2009).Mathematical modeling of plant metabolism offers new approaches to understand, predict, and modify complex plant metabolic processes. In plant research, the issue of metabolic modeling is constantly gaining attention, and different modeling approaches applied to plant metabolism exist, ranging from highly detailed quantitative to less complex qualitative approaches (for review, see Giersch, 2000; Morgan and Rhodes, 2002; Poolman et al., 2004; Rios-Estepa and Lange, 2007).A widely used modeling approach is flux balance analysis (FBA), which allows the prediction of metabolic capabilities and steady-state fluxes under different environmental and genetic backgrounds using (non)linear optimization (Orth et al., 2010). Assuming steady-state conditions, FBA has the advantage of not requiring the knowledge of kinetic parameters and, therefore, can be applied to model detailed, large-scale systems. In recent years, the FBA approach has been applied to several different plant species, such as maize (Zea mays; Dal’Molin et al., 2010; Saha et al., 2011), barley (Hordeum vulgare; Grafahrend-Belau et al., 2009b; Melkus et al., 2011; Rolletschek et al., 2011), rice (Oryza sativa; Lakshmanan et al., 2013), Arabidopsis (Arabidopsis thaliana; Poolman et al., 2009; de Oliveira Dal’Molin et al., 2010; Radrich et al., 2010; Williams et al., 2010; Mintz-Oron et al., 2012; Cheung et al., 2013), and rapeseed (Brassica napus; Hay and Schwender, 2011a, 2011b; Pilalis et al., 2011), as well as algae (Boyle and Morgan, 2009; Cogne et al., 2011; Dal’Molin et al., 2011) and photoautotrophic bacteria (Knoop et al., 2010; Montagud et al., 2010; Boyle and Morgan, 2011). These models have been used to study different aspects of metabolism, including the prediction of optimal metabolic yields and energy efficiencies (Dal’Molin et al., 2010; Boyle and Morgan, 2011), changes in flux under different environmental and genetic backgrounds (Grafahrend-Belau et al., 2009b; Dal’Molin et al., 2010; Melkus et al., 2011), and nonintuitive metabolic pathways that merit subsequent experimental investigations (Poolman et al., 2009; Knoop et al., 2010; Rolletschek et al., 2011). Although FBA of plant metabolic models was shown to be capable of reproducing experimentally determined flux distributions (Williams et al., 2010; Hay and Schwender, 2011b) and generating new insights into metabolic behavior, capacities, and efficiencies (Sweetlove and Ratcliffe, 2011), challenges remain to advance the utility and predictive power of the models.Given that many plant metabolic functions are based on interactions between different subcellular compartments, cell types, tissues, and organs, the reconstruction of organ-specific models and the integration of these models into interacting multiorgan and/or whole-plant models is a prerequisite to get insight into complex plant metabolic processes organized on a whole-plant scale (e.g. source-sink interactions). Almost all FBA models of plant metabolism are restricted to one cell type (Boyle and Morgan, 2009; Knoop et al., 2010; Montagud et al., 2010; Cogne et al., 2011; Dal’Molin et al., 2011), one tissue or one organ (Grafahrend-Belau et al., 2009b; Hay and Schwender, 2011a, 2011b; Pilalis et al., 2011; Mintz-Oron et al., 2012), and only one model exists taking into account the interaction between two cell types by specifying the interaction between mesophyll and bundle sheath cells in C4 photosynthesis (Dal’Molin et al., 2010). So far, no model representing metabolism at the whole-plant scale exists.Considering whole-plant metabolism raises the problem of taking into account temporal and environmental changes in metabolism during plant development and growth. Although classical static FBA is unable to predict the dynamics of metabolic processes, as the network analysis is based on steady-state solutions, time-dependent processes can be taken into account by extending the classical static FBA to a dynamic flux balance analysis (dFBA), as proposed by Mahadevan et al. (2002). The static (SOA) and dynamic optimization approaches introduced in this work provide a framework for analyzing the transience of metabolism by integrating kinetic expressions to dynamically constrain exchange fluxes. Due to the requirement of knowing or estimating a large number of kinetic parameters, so far dFBA has only been applied to a plant metabolic model once, to study the photosynthetic metabolism in the chloroplasts of C3 plants by a simplified model of five biochemical reactions (Luo et al., 2009). Integrating a dynamic model into a static FBA model is an alternative approach to perform dFBA.In this study, a multiscale metabolic modeling (MMM) approach was applied with the aim of achieving a spatiotemporal resolution of cereal crop plant metabolism. To provide a framework for the in silico analysis of the metabolic dynamics of barley on a whole-plant scale, the MMM approach integrates a static multiorgan FBA model and a dynamic whole-plant multiscale functional plant model (FPM) to perform dFBA. The performance of the novel whole-plant MMM approach was tested by studying source-sink interactions during the seed developmental phase of barley plants.     

Capturing Arabidopsis Root Architecture Dynamics with root-fit Reveals Diversity in Responses to Salinity

The plant root is the first organ to encounter salinity stress, but the effect of salinity on root system architecture (RSA) remains elusive. Both the reduction in main root (MR) elongation and the redistribution of the root mass between MRs and lateral roots (LRs) are likely to play crucial roles in water extraction efficiency and ion exclusion. To establish which RSA parameters are responsive to salt stress, we performed a detailed time course experiment in which Arabidopsis (Arabidopsis thaliana) seedlings were grown on agar plates under different salt stress conditions. We captured RSA dynamics with quadratic growth functions (root-fit) and summarized the salt-induced differences in RSA dynamics in three growth parameters: MR elongation, average LR elongation, and increase in number of LRs. In the ecotype Columbia-0 accession of Arabidopsis, salt stress affected MR elongation more severely than LR elongation and an increase in LRs, leading to a significantly altered RSA. By quantifying RSA dynamics of 31 different Arabidopsis accessions in control and mild salt stress conditions, different strategies for regulation of MR and LR meristems and root branching were revealed. Different RSA strategies partially correlated with natural variation in abscisic acid sensitivity and different Na⁺/K⁺ ratios in shoots of seedlings grown under mild salt stress. Applying root-fit to describe the dynamics of RSA allowed us to uncover the natural diversity in root morphology and cluster it into four response types that otherwise would have been overlooked.Salt stress is known to affect plant growth and productivity as a result of its osmotic and ionic stress components. Osmotic stress imposed by salinity is thought to act in the early stages of the response, by reducing cell expansion in growing tissues and causing stomatal closure to minimize water loss. The build-up of ions in photosynthetic tissues leads to toxicity in the later stages of salinity stress and can be reduced by limiting sodium transport into the shoot tissue and compartmentalization of sodium ions into the root stele and vacuoles (Munns and Tester, 2008). The effect of salt stress on plant development was studied in terms of ion accumulation, plant survival, and signaling (Munns et al., 2012; Hasegawa, 2013; Pierik and Testerink, 2014). Most studies focus on traits in the aboveground tissues, because minimizing salt accumulation in leaf tissue is crucial for plant survival and its productivity. This approach has led to the discovery of many genes underlying salinity tolerance (Munns and Tester, 2008; Munns et al., 2012; Hasegawa, 2013; Maathuis, 2014). Another way to estimate salinity stress tolerance is by studying the rate of main root (MR) elongation of seedlings transferred to medium supplemented with high salt concentration. This is how Salt Overly Sensitive mutants were identified, being a classical example of genes involved in salt stress signaling and tolerance (Hasegawa, 2013; Maathuis, 2014). The success of this approach is to be explained by the important role that the root plays in salinity tolerance. Roots not only provide anchorage and ensure water and nutrient uptake, but also act as a sensory system, integrating changes in nutrient availability, water content, and salinity to adjust root morphology to exploit available resources to the maximum capacity (Galvan-Ampudia et al., 2013; Gruber et al., 2013). Understanding the significance of environmental modifications of root system architecture (RSA) for plant productivity is one of the major challenges of modern agriculture (de Dorlodot et al., 2007; Den Herder et al., 2010; Pierik and Testerink, 2014).The RSA of dicotyledonous plants consists of an embryonically derived MR and lateral roots (LRs) that originate from xylem pole pericycle cells of the MR, or from LRs in the case of higher-order LRs. Root growth and branching is mainly guided through the antagonistic action of two plant hormones: auxin and cytokinins (Petricka et al., 2012). Under environmental stress conditions, the synthesis of abscisic acid (ABA), ethylene, and brassinosteroids is known to be induced and to modulate the growth of MRs and LRs (Achard et al., 2006; Osmont et al., 2007; Achard and Genschik, 2009; Duan et al., 2013; Geng et al., 2013). In general, lower concentrations of salt were observed to slightly induce MR and LR elongation, whereas higher concentrations resulted in decreased growth of both MRs and LRs (Wang et al., 2009; Zolla et al., 2010). The reduction of growth is a result of the inhibition of cell cycle progression and a reduction in root apical meristem size (West et al., 2004). However, conflicting results were presented for the effect of salinity on lateral root density (LRD; Wang et al., 2009; Zolla et al., 2010; Galvan-Ampudia and Testerink, 2011). Some studies suggest that mild salinity enhances LR initiation or emergence events, thereby affecting patterning, whereas other studies imply that salinity arrests LR development. The origin of those contradictory observations could be attributable to studying LR initiation and density at single time points, rather than observing the dynamics of LR development, because LR formation changes as a function of root growth rate (De Smet et al., 2012). The dynamics of LR growth and development were characterized previously for the MR region formed before the salt stress exposure, identifying the importance of ABA in early growth arrest of postemerged LRs in response to salt stress (Duan et al., 2013). The effect of salt on LR emergence and initiation was found to differ for MR regions formed prior and subsequent to salinity exposure (Duan et al., 2013), consistent with LR patterning being determined at the root tip (Moreno-Risueno et al., 2010). Yet the effect of salt stress on the reprogramming of the entire RSA on a longer timescale remains elusive.Natural variation in Arabidopsis (Arabidopsis thaliana) is a great source for dissecting the genetic components underlying phenotypic diversity (Trontin et al., 2011; Weigel, 2012). Genes underlying phenotypic plasticity of RSA to environmental stimuli were also found to have high allelic variation leading to differences in root development between different Arabidopsis accessions (Rosas et al., 2013). Supposedly, genes responsible for phenotypic plasticity of the root morphology to different environmental conditions are under strong selection for adaptation to local environments. Various populations of Arabidopsis accessions were used to study natural variation in ion accumulation and salinity tolerance (Rus et al., 2006; Jha et al., 2010; Katori et al., 2010; Roy et al., 2013). In addition, a number of studies focusing on the natural variation in RSA have been published, identifying quantitative trait loci and allelic variation for genes involved in RSA development under control conditions (Mouchel et al., 2004; Meijón et al., 2014) and nutrient-deficient conditions (Chevalier et al., 2003; Gujas et al., 2012; Gifford et al., 2013; Kellermeier et al., 2013; Rosas et al., 2013). Exploring natural variation not only expands the knowledge of genes and molecular mechanisms underlying biological processes, but also provides insight on how plants adapt to challenging environmental conditions (Weigel, 2012) and whether the mechanisms are evolutionarily conserved. The early growth arrest of newly emerged LRs upon exposure to salt stress was observed to be conserved among the most commonly used Arabidopsis accessions Columbia-0 (Col-0), Landsberg erecta, and Wassilewskija (Ws; Duan et al., 2013). By studying salt stress responses of the entire RSA and a wider natural variation in root responses to stress, one could identify new morphological traits that are under environmental selection and possibly contribute to stress tolerance.In this work, we not only identify the RSA components that are responsive to salt stress, but we also describe the natural variation in dynamics of salt-induced changes leading to redistribution of root mass and different root morphology. The growth dynamics of MRs and LRs under different salt stress conditions were described by fitting a set of quadratic growth functions (root-fit) to individual RSA components. Studying salt-induced changes in RSA dynamics of 31 Arabidopsis accessions revealed four major strategies conserved among the accessions. Those four strategies were due to differences in salt stress sensitivity of individual RSA components (i.e. growth rates of MRs and LRs, and increases in the number of emerged LRs). This diversity in root morphology responses caused by salt stress was observed to be partially associated with differences in ABA, but not ethylene sensitivity. In addition, we observed that a number of accessions exhibiting a relatively strong inhibition of LR elongation showed a smaller increase in the Na⁺/K⁺ ratio in shoot tissue after exposure to salt stress. Our results imply that different RSA strategies identified in this study reflect diverse adaptations to different soil conditions and thus might contribute to efficient water extraction and ion compartmentalization in their native environments.     

Quantitative 3D Analysis of Plant Roots Growing in Soil Using Magnetic Resonance Imaging

Precise measurements of root system architecture traits are an important requirement for plant phenotyping. Most of the current methods for analyzing root growth require either artificial growing conditions (e.g. hydroponics), are severely restricted in the fraction of roots detectable (e.g. rhizotrons), or are destructive (e.g. soil coring). On the other hand, modalities such as magnetic resonance imaging (MRI) are noninvasive and allow high-quality three-dimensional imaging of roots in soil. Here, we present a plant root imaging and analysis pipeline using MRI together with an advanced image visualization and analysis software toolbox named NMRooting. Pots up to 117 mm in diameter and 800 mm in height can be measured with the 4.7 T MRI instrument used here. For 1.5 l pots (81 mm diameter, 300 mm high), a fully automated system was developed enabling measurement of up to 18 pots per day. The most important root traits that can be nondestructively monitored over time are root mass, length, diameter, tip number, and growth angles (in two-dimensional polar coordinates) and spatial distribution. Various validation measurements for these traits were performed, showing that roots down to a diameter range between 200 μm and 300 μm can be quantitatively measured. Root fresh weight correlates linearly with root mass determined by MRI. We demonstrate the capabilities of MRI and the dedicated imaging pipeline in experimental series performed on soil-grown maize (Zea mays) and barley (Hordeum vulgare) plants.The root system is of critical importance for the survival, development, and performance of higher plants, because it is the major organ for anchorage, acquisition of water and nutrients, and carbon storage. Therefore, there is a long-standing interest in the scientific community regarding the structure, function, and development of root systems. Rising concerns about the environmental impact and increasing cost of fertilizers have initiated breeding programs for more resource-efficient cultivars and the development of methods for phenotyping root systems. The opaque nature of soils generally demands destructive methods such as root excavation for subsequent optical assessment (Lynch, 2007; Trachsel et al., 2011). Although efficient for screening large numbers of plants for a limited set of clearly discernible traits, this approach does not allow detailed monitoring of root development over time. Other approaches, such as rhizotrons or mini-rhizotron tubes, where root growth is observed along transparent windows (Nagel et al., 2009), monitor only a fraction of the roots. Methods in which the whole root system is visible are typically based on artificial media such as paper pouches (Chen et al., 2011; Le Marié et al., 2014), three-dimensional (3D) gels (Iyer-Pascuzzi et al., 2010), and hydro- or aeroponics (Herdel et al., 2001). Results may thus not be directly transferable to plants grown in natural 3D soil environments (Gregory et al., 2003). For example, roots are known to grow faster and thinner when the penetration resistance is low (Bengough et al., 2011; Chimungu et al., 2015). Computed tomography (CT; both x-ray and neutron) has been proposed to overcome the mentioned difficulties with studying roots in natural soil. CT has been successfully used to obtain high-resolution images of roots (Moradi et al., 2009; Flavel et al., 2012; Mooney et al., 2012). High resolution is necessary for segmenting roots due to a poor contrast between roots and soil (Jassogne et al., 2009; Mairhofer et al., 2012; Mairhofer et al., 2013). A first direct comparison (to our knowledge) of magnetic resonance imaging (MRI) and x-ray CT for 3D root imaging has recently been published (Metzner et al., 2015), showing that the two modalities pose different opportunities and limitations for root imaging.MRI is based on the magnetic moment of atomic nuclei like ¹H (protons), which are highly abundant in living tissues, mainly in water molecules. The magnetic moment can be manipulated using strong magnetic and radio frequency fields that have no known impact on plant development to produce 3D datasets of samples. MRI offers several contrast parameters that can be manipulated for discriminating different structures such as roots from soil background (Rogers and Bottomley, 1987; Jahnke et al., 2009). The basic principles of MRI are described in detail in several textbooks (Callaghan, 1993; Haacke et al., 1999) or review articles (Köckenberger et al., 2004; Blümler et al., 2009; van As et al., 2009; Borisjuk et al., 2012). Research applications to plant roots range from phytopathology (Hillnhütter et al., 2012), across storage root internal structures (Metzner et al., 2014) and water uptake modeling (Stingaciu et al., 2013), to coregistration with positron emission tomography for investigating structure-function relations (Jahnke et al., 2009). Water mobility in roots and soil has also been shown to be detectable with MRI (MacFall and Johnson, 2012; Gruwel, 2014). In particular for imaging roots with MRI, these studies generally explored the applicability of MRI but largely lacked validation of the data against conventional techniques of root visualization after harvest. Our goal was to develop MRI protocols to image roots of plants growing in soil to obtain global root parameters such as root length, mass, or root diameters; gather root growth angles and number of root tips; get spatial information on the distribution of root system architecture (RSA) parameters such as root length densities; and, wherever possible, verify these parameters against harvest data.