Leaf Growth Response to Mild Drought: Natural Variation in Arabidopsis Sheds Light on Trait Architecture

Dual and spatially resolved drought responses in the Arabidopsis leaf mesophyll revealed by single-cell transcriptomics

Drought stress imposes severe challenges on agriculture by impacting crop performance. Understanding drought responses in plants at a cellular level is a crucial first step toward engineering improved drought resilience. However, the molecular responses to drought are complex as they depend on multiple factors, including the severity of drought, the profiled organ, its developmental stage or even the cell types therein. Thus, deciphering the transcriptional responses to drought is especially challenging. In this study, we investigated tissue-specific responses to mild drought (MD) in young Arabidopsis thaliana (Arabidopsis) leaves using single-cell RNA sequencing (scRNA-seq). To preserve transcriptional integrity during cell isolation, we inhibited RNA synthesis using the transcription inhibitor actinomycin D, and demonstrated the benefits of transcriptome fixation for studying mild stress responses at a single-cell level. We present a curated and validated single-cell atlas, comprising 50 797 high-quality cells from almost all known cell types present in the leaf. All cell type annotations were validated with a new library of reporter lines. The curated data are available to the broad community in an intuitive tool and a browsable single-cell atlas (http://www.single-cell.be/plant/leaf-drought). We show that the mesophyll contains two spatially separated cell populations with distinct responses to drought: one enriched in canonical abscisic acid-related drought-responsive genes, and another one enriched in genes involved in iron starvation responses. Our study thus reveals a dual adaptive mechanism of the leaf mesophyll in response to MD and provides a valuable resource for future research on stress responses.

The determination of leaf size on the basis of developmental traits

Mature leaf area (LA) is a showcase of diversity - varying enormously within and across species, and associated with the productivity and distribution of plants and ecosystems. Yet, it remains unclear how developmental processes determine variation in LA. We introduce a mathematical framework pinpointing the origin of variation in LA by quantifying six epidermal 'developmental traits': initial mean cell size and number (approximating values within the leaf primordium), and the maximum relative rates and durations of cell proliferation and expansion until leaf maturity. We analyzed a novel database of developmental trajectories of LA and epidermal anatomy, representing 12 eudicotyledonous species and 52 Arabidopsis experiments. Within and across species, mean primordium cell number and maximum relative cell proliferation rate were the strongest developmental determinants of LA. Trade-offs between developmental traits, consistent with evolutionary and metabolic scaling theory, strongly constrain LA variation. These include trade-offs between primordium cell number vs cell proliferation, primordium mean cell size vs cell expansion, and the durations vs maximum relative rates of cell proliferation and expansion. Mutant and wild-type comparisons showed these trade-offs have a genetic basis in Arabidopsis. Analyses of developmental traits underlying LA and its diversification highlight mechanisms for leaf evolution, and opportunities for breeding trait shifts.

Mitogen-Activated Protein Kinases 3/6 Reduce Auxin Signaling via Stabilizing Indoleacetic Acid-Induced Proteins 8/9 in Plant Abiotic Stress Adaptation

The balance between plant growth and stress response is a key issue in the field of biology. In this process, mitogen-activated protein kinase 3 (MPK3) and MPK6 contribute to the construction of plants' defense system during stress tolerance, while auxin, a growth-promoting hormone, is the key to maintaining plant growth. Nevertheless, the antagonistic or cooperative relationship between MPK3/6-mediated stress response and auxin-mediated plant growth remains unclear. Here, we demonstrate that stress-activated MPK3/6 interact with the auxin signaling repressors indoleacetic acid-induced protein 8 (IAA8) and IAA9, two key targets for regulating the auxin signaling output during stress responses. Protein phosphorylation mass spectrometry followed by a co-analysis with in vitro phosphorylation experiments revealed that MPK3/6 phosphorylated the S91, T94, and S152 residues of IAA8 and the S88 residue of IAA9. Phosphorylation significantly enhanced the protein stability of IAA8/9, thereby maintaining basal auxin signaling during plant stress adaptation. Collectively, MPK3/6-IAA8/9 are key modules that are turned on during plant stress adaptation to precisely reduce auxin signaling output, thereby preventing plants from improper vigorous growth under stress conditions.

Tensor decomposition reveals trans-regulated gene modules in maize drought response

When plants respond to drought stress, dynamic cellular changes occur, accompanied by alterations in gene expression, which often act through trans-regulation. However, the detection of trans-acting genetic variants and networks of genes is challenged by the large number of genes and markers. Using a tensor decomposition method, we identify trans-acting expression quantitative trait loci (trans-eQTLs) linked to gene modules, rather than individual genes, which were associated with maize drought response. Module-to-trait association analysis demonstrates that half of the modules are relevant to drought-related traits. Genome-wide association studies of the expression patterns of each module identify 286 trans-eQTLs linked to drought-responsive modules, the majority of which cannot be detected based on individual gene expression. Notably, the trans-eQTLs located in the regions selected during maize improvement tend towards relatively strong selection. We further prioritize the genes that affect the transcriptional regulation of multiple genes in trans, as exemplified by two transcription factor genes. Our analyses highlight that multidimensional reduction could facilitate the identification of trans-acting variations in gene expression in response to dynamic environments and serve as a promising technique for high-order data processing in future crop breeding.

AtZAT10/STZ1 improves drought tolerance and increases fiber yield in cotton

Drought poses a significant challenge to global crop productivity, necessitating innovative approaches to bolster plant resilience. Leveraging transgenic technology to bolster drought tolerance in crops emerges as a promising strategy for addressing the demands of a rapidly growing global populace. AtZAT10/STZ1, a C2H2-type zinc finger protein transcription factor has shown to significantly improve Arabidopsis' tolerance to various abiotic stresses. In this study, we reports that AtSTZ1 confers notable drought resistance in upland cotton (Gossypium hirsutum), amplifying cotton fiber yield under varying conditions, including irrigated and water-limited environments, in field trials. Notably, AtSTZ1-overexpressing transgenic cotton showcases enhanced drought resilience across critical growth stages, including seed germination, seedling establishment, and reproductive phases. Morphological analysis reveals an expanded root system characterized by an elongated taproot system, increased lateral roots, augmented root biomass, and enlarged cell dimensions from transgenic cotton plants. Additionally, higher contents of proline, chlorophyll, soluble sugars, and enhanced ROS-scavenging enzyme activities are observed in leaves of transgenic plants subjected to drought, underscoring improved physiological adaptations. Furthermore, transgenic lines exhibit heightened photosynthetic rate, increased water use efficiency, and larger stomatal and epidermal cell sizes, coupled with a decline in leaf stomatal conductance and density, as well as diminished transpiration rates compared to the wild type counterparts. Transcriptome profiling unveils 106 differentially expressed genes in transgenic cotton leaves post-drought treatment, including protein kinases, transcription factors, aquaporins, and heat shock proteins, indicative of an orchestrated stress response. Collectively, these findings underscore the capacity of AtSTZ1 to augment the expression of abiotic stress-related genes in cotton following drought conditions, thus presenting a compelling candidate for genetic manipulation aimed at enhancing crop resilience.

Natural variation in response to combined water and nitrogen deficiencies in Arabidopsis

Understanding plant responses to individual stresses does not mean that we understand real-world situations, where stresses usually combine and interact. These interactions arise at different levels, from stress exposure to the molecular networks of the stress response. Here, we built an in-depth multiomic description of plant responses to mild water (W) and nitrogen (N) limitations, either individually or combined, among 5 genetically different Arabidopsis (Arabidopsis thaliana) accessions. We highlight the different dynamics in stress response through integrative traits such as rosette growth and the physiological status of the plants. We also used transcriptomic and metabolomic profiling during a stage when the plant response was stabilized to determine the wide diversity in stress-induced changes among accessions, highlighting the limited reality of a "universal" stress response. The main effect of the W × N interaction was an attenuation of the N-deficiency syndrome when combined with mild drought, but to a variable extent depending on the accession. Other traits subject to W × N interactions are often accession specific. Multiomic analyses identified a subset of transcript-metabolite clusters that are critical to stress responses but essentially variable according to the genotype factor. Including intraspecific diversity in our descriptions of plant stress response places our findings in perspective.

Alternative splicing as a driver of natural variation in abscisic acid response

Abscisic acid (ABA) is a crucial player in plant responses to the environment. It accumulates under stress, activating downstream signaling to implement molecular responses that restore homeostasis. Natural variance in ABA sensitivity remains barely understood, and the ABA pathway has been mainly studied at the transcriptional level, despite evidence that posttranscriptional regulation, namely, via alternative splicing, contributes to plant stress tolerance. Here, we identified the Arabidopsis accession Kn-0 as less sensitive to ABA than the reference Col-0, as shown by reduced effects of the hormone on seedling establishment, root branching, and stomatal closure, as well as by decreased induction of ABA marker genes. An in-depth comparative transcriptome analysis of the ABA response in the two variants revealed lower expression changes and fewer genes affected for the least ABA-sensitive ecotype. Notably, Kn-0 exhibited reduced levels of the ABA-signaling SnRK2 protein kinases and lower basal expression of ABA-reactivation genes, consistent with our finding that Kn-0 contains less endogenous ABA than Col-0. ABA also markedly affected alternative splicing, primarily intron retention, with Kn-0 being less responsive regarding both the number and magnitude of alternative splicing events, particularly exon skipping. We find that alternative splicing introduces a more ecotype-specific layer of ABA regulation and identify ABA-responsive splicing changes in key ABA pathway regulators that provide a functional and mechanistic link to the differential sensitivity of the two ecotypes. Our results offer new insight into the natural variation of ABA responses and corroborate a key role for alternative splicing in implementing ABA-mediated stress responses.

Genome-wide association analysis reveals genes controlling an antagonistic effect of biotic and osmotic stress on Arabidopsis thaliana growth

While the response of Arabidopsis thaliana to drought, herbivory or fungal infection has been well-examined, the consequences of exposure to a series of such (a)biotic stresses are not well studied. This work reports on the genetic mechanisms underlying the Arabidopsis response to single osmotic stress, and to combinatorial stress, either fungal infection using Botrytis cinerea or herbivory using Pieris rapae caterpillars followed by an osmotic stress treatment. Several small-effect genetic loci associated with rosette dry weight (DW), rosette water content (WC), and the projected rosette leaf area in response to combinatorial stress were mapped using univariate and multi-environment genome-wide association approaches. A single-nucleotide polymorphism (SNP) associated with DROUGHT-INDUCED 19 (DI19) was identified by both approaches, supporting its potential involvement in the response to combinatorial stress. Several SNPs were found to be in linkage disequilibrium with known stress-responsive genes such as PEROXIDASE 34 (PRX34), BASIC LEUCINE ZIPPER 25 (bZIP25), RESISTANCE METHYLATED GENE 1 (RMG1) and WHITE RUST RESISTANCE 4 (WRR4). An antagonistic effect between biotic and osmotic stress was found for prx34 and arf4 mutants, which suggests PRX34 and ARF4 play an important role in the response to the combinatorial stress.

Variability in drought gene expression datasets highlight the need for community standardization

Physiologically relevant drought stress is difficult to apply consistently, and the heterogeneity in experimental design, growth conditions, and sampling schemes make it challenging to compare water deficit studies in plants. Here, we re-analyzed hundreds of drought gene expression experiments across diverse model and crop species and quantified the variability across studies. We found that drought studies are surprisingly uncomparable, even when accounting for differences in genotype, environment, drought severity, and method of drying. Many studies, including most Arabidopsis work, lack high-quality phenotypic and physiological datasets to accompany gene expression, making it impossible to assess the severity or in some cases the occurrence of water deficit stress events. From these datasets, we developed supervised learning classifiers that can accurately predict if RNA-seq samples have experienced a physiologically relevant drought stress, and suggest this can be used as a quality control for future studies. Together, our analyses highlight the need for more community standardization, and the importance of paired physiology data to quantify stress severity for reproducibility and future data analyses.

Natural variation in the SVP contributes to the pleiotropic adaption of Arabidopsis thaliana across contrasted habitats

Coordinated plant adaptation involves the interplay of multiple traits driven by habitat-specific selection pressures. Pleiotropic effects, wherein genetic variants of a single gene control multiple traits, can expedite such adaptations. Until present, only a limited number of genes have been reported to exhibit pleiotropy. Here, we create a recombinant inbred line (RIL) population derived from two Arabidopsis thaliana (A. thaliana) ecotypes originating from divergent habitats. Using this RIL population, we identify an allelic variation in a MADS-box transcription factor, SHORT VEGETATIVE PHASE (SVP), which exerts a pleiotropic effect on leaf size and drought-versus-humidity tolerance. Further investigation reveals that a natural null variant of the SVP protein disrupts its normal regulatory interactions with target genes, including GRF3, CYP707A1/3, and AtBG1, leading to increased leaf size, enhanced tolerance to humid conditions, and changes in flowering time of humid conditions in A. thaliana. Remarkably, polymorphic variations in this gene have been traced back to early A. thaliana populations, providing a genetic foundation and plasticity for subsequent colonization of diverse habitats by influencing multiple traits. These findings advance our understanding of how plants rapidly adapt to changing environments by virtue of the pleiotropic effects of individual genes on multiple trait alterations.

Genome-wide mapping of regulatory variants for temperature- and salinity-adaptive genes reveals genetic basis of genotype-by-environment interaction in Crassostrea ariakensis

Regulatory variants in gene expression serve as bridges linking genetic variation and phenotypic plasticity. Environmental conditions typically influence the effects of regulatory variants on phenotypic plasticity; however, such genotype-by-environment interactions (G × E) are poorly understood. This study aimed to investigate the genetic basis of G × E in estuarine oyster (Crassostrea ariakensis), which is an important model animal for studying environmental adaption owing to its high plasticity and large intraspecific divergence. Genome-wide mapping of expression quantitative trait loci (eQTLs) for 23 environmental adaptive genes was performed for 256 estuarine oysters. We identified 1194 eQTL single nucleotide polymorphisms (eSNPs), including 433 cis-eSNPs in four genes and 722 trans-eSNPs in eight genes. The expression variation explanation of cis-eSNPs (9.95%) was significantly higher than that of trans-eSNPs (9.15%). We specifically showed cis- and trans-eSNPs with high linkage disequilibrium (LD) for Traf7, Slc6a5, Ggt, and Dap3. For example, we identified a cis-regulatory LD block containing 68 cis-eSNP and a trans-regulatory LD block, including 20 trans-eSNPs in Traf7. A high proportion (85%) of 40 vital eSNPs exhibited significant G × E effects. We identified crossing and nonparallel interactions of G × E, with the tag cis-eSNPs of Baat and Slc6a5 as representatives. Our results indicated that cis-eQTLs are highly conserved. This study provides insights into the understanding of adaptive evolutionary mechanisms and phenotypic response prediction to variable environments, as well as the genetic improvement for superior adaptive traits for genetic resource conservation and aquaculture.

Phosphate starvation response precedes abscisic acid response under progressive mild drought in plants

Drought severely damages crop production, even under conditions so mild that the leaves show no signs of wilting. However, it is unclear how field-grown plants respond to mild drought. Here, we show through six years of field trials that ridges are a useful experimental tool to mimic mild drought stress in the field. Mild drought reduces inorganic phosphate levels in the leaves to activate the phosphate starvation response (PSR) in soybean plants in the field. Using Arabidopsis thaliana and its mutant plants grown in pots under controlled environments, we demonstrate that PSR occurs before abscisic acid response under progressive mild drought and that PSR plays a crucial role in plant growth under mild drought. Our observations in the field and laboratory using model crop and experimental plants provide insight into the molecular response to mild drought in field-grown plants and the relationship between nutrition and drought stress response.

MdbZIP74 negatively regulates osmotic tolerance and adaptability to moderate drought conditions of apple plants

Drought is the most prominent threat to global agricultural production. The basic leucine zipper (bZIP) family is related to the response to a series of abiotic stress. In this case, apple calli and the seedlings of MdbZIP74-RNAi transgenic lines were obtained. Under osmotic stress and moderate drought conditions, the content of malondialdehyde, relative water content and other stress-related assays were measured. MdbZIP74 was found to negatively regulate the osmotic tolerance of apple callus. The growth of MdbZIP74-RNAi calli enhanced resistance without significant production loss. The silencing of MdbZIP74 contributes to redox balance and the adaptability of apple seedlings to moderate drought conditions. Four related differentially expressed genes in the biosynthesis of cytokinin and catabolic pathway were identified through a transcriptome analysis of MdbZIP74-RNAi seedlings under moderate drought conditions. MdLOG8 was further identified as the target of MdbZIP74 involved in the drought adaptability of apple plants using a dual experiment. Further confirmation showed MdLOG8 was maintained in the MdbZIP74-RNAi seedlings presumably acting as the growth regulator to enhance drought adaptability. It was concluded that the correct regulation of cytokinin level under moderate drought conditions maintains the redox balance and avoids the situation of plants surviving with the minimal resources.

Predicting Gene Regulatory Interactions Using Natural Genetic Variation

Genome-wide association studies (GWAS) are a powerful tool to elucidate the genotype-phenotype map. Although GWAS are usually used to assess simple univariate associations between genetic markers and traits of interest, it is also possible to infer the underlying genetic architecture and to predict gene regulatory interactions. In this chapter, we describe the latest methods and tools to perform GWAS by calculating permutation-based significance thresholds. For this purpose, we first provide guidelines on univariate GWAS analyses that are extended in the second part of this chapter to more complex models that enable the inference of gene regulatory networks and how these networks vary.

StATL2-like could affect growth and cold tolerance of plant by interacting with StCBFs

Our results confirmed that StATL2-like could interact with StCBFs and regulate plant growth. Meanwhile, StATL2-like acted as a negative regulator on low-temperature tolerance in plants. As important transcription factors for resisting many kinds of stresses, C-repeat-binding factors (CBF) play a key role in plant low-temperature tolerance by increasing COR genes expressions. Here, we report that StATL2-like, a RING-H2 E3 ubiquitin in Solanum tuberosum L., interacted with StCBF1 and StCBF4, respectively. AtATL2 is a highly homologous gene of StATL2-like in Arabidopsis thaliana. Under normal conditions, atl2 Arabidopsis mutant showed a growth inhibition phenotype while overexpressed StATL2-like in wild type Arabidopsis and atl2 mutant promoted plant growth. Besides, atl2 mutant had better low-temperature tolerance compared with wild type and StATL2-like transgenic lines which demonstrated that StATL2-like acted as a negatively regulator on low-temperature tolerance in plant. Moreover, atl2 mutant improved the scavenging capacity of reactive oxygen species (ROS) and alleviate the damage of photosynthetic system II (PSII) compared with StATL2-like transgenic lines under cold conditions. These results suggested a new component in CBF-dependent pathway to regulate plant growth and response to low-temperature stress in potato plants.

Coordination of plant growth and abiotic stress responses by tryptophan synthase β subunit 1 through modulation of tryptophan and ABA homeostasis in Arabidopsis

To adapt to changing environments, plants have evolved elaborate regulatory mechanisms balancing their growth with stress responses. It is currently unclear whether and how the tryptophan (Trp), the growth-related hormone auxin, and the stress hormone abscisic acid (ABA) are coordinated in this trade-off. Here, we show that tryptophan synthase β subunit 1 (TSB1) is involved in the coordination of Trp and ABA, thereby affecting plant growth and abiotic stress responses. Plants experiencing high salinity or drought display reduced TSB1 expression, resulting in decreased Trp and auxin accumulation and thus reduced growth. In comparison with the wild type, amiR-TSB1 lines and TSB1 mutants exhibited repressed growth under non-stress conditions but had enhanced ABA accumulation and stress tolerance when subjected to salt or drought stress. Furthermore, we found that TSB1 interacts with and inhibits β-glucosidase 1 (BG1), which hydrolyses glucose-conjugated ABA into active ABA. Mutation of BG1 in the amiR-TSB1 lines compromised their increased ABA accumulation and enhanced stress tolerance. Moreover, stress-induced H₂O₂ disrupted the interaction between TSB1 and BG1 by sulfenylating cysteine-308 of TSB1, relieving the TSB1-mediated inhibition of BG1 activity. Taken together, we revealed that TSB1 serves as a key coordinator of plant growth and stress responses by balancing Trp and ABA homeostasis.

Past accomplishments and future challenges of the multi-omics characterization of leaf growth

The advent of omics technologies has revolutionized biology and advanced our understanding of all biological processes, including major developmental transitions in plants and animals. Here, we review the vast knowledge accumulated concerning leaf growth in terms of transcriptional regulation before turning our attention to the historically less well-characterized alterations at the protein and metabolite level. We will then discuss how the advent of biochemical methods coupled with metabolomics and proteomics can provide insight into the protein-protein and protein-metabolite interactome of the growing leaves. We finally highlight the substantial challenges in detection, spatial resolution, integration, and functional validation of the omics results, focusing on metabolomics as a prerequisite for a comprehensive understanding of small-molecule regulation of plant growth.

Water Shortage Affects Vegetative and Reproductive Stages of Common Bean (Phaseolus vulgaris) Chilean Landraces, Differentially Impacting Grain Yield Components

Water availability for agricultural use is currently a global problem that worsens with climate change in several regions of the world. Among grain legumes, common bean (Phaseolus vulgaris) is the most cultivated in the worldwide. The Chilean germplasm of common bean is characterized by tolerance to water stress. Here, we analyzed a selection of nine ancient Chilean landraces in regard to their drought tolerance, simulating optimal (OW) and restricted watering (RW) in a Mediterranean environment. Phenological, growth, and yield traits were recorded, and correlation analysis was performed. Accordingly, leaf temperature and osmotic potential were higher under RW, while the leaf chlorophyll content decreased in all landraces. Physiological maturity days and seed-filling days were lower in RW than in OW. This similarly occurred with the grain yield. The % yield reduction was negatively correlated with the % pod reduction and the relative rate of leaf expansion (RLAE) reduction. However, the 100-seed weight value was not significantly modified by water treatment (p > 0.05). For instance, landraces that preferred to fill the grain with a lower rate of leaf expansion showed a lower loss in grain yield under drought conditions. These results suggest that the resource partitioning between growing leaves, flowers, and developing pods in Chilean landraces is variable, affecting the common bean drought tolerance.

Jasmonate-Dependent Response of the Flower Abscission Zone Cells to Drought in Yellow Lupine

Lipid membranes, as primary places of the perception of environmental stimuli, are a source of various oxygenated polyunsaturated fatty acids-oxylipins-functioning as modulators of many signal transduction pathways, e.g., phytohormonal. Among exogenous factors acting on plant cells, special attention is given to drought, especially in highly sensitive crop species, such as yellow lupine. Here, we used this species to analyze the contribution of lipid-related enzymes and lipid-derived plant hormones in drought-evoked events taking place in a specialized group of cells-the flower abscission zone (AZ)-which is responsible for organ detachment from the plant body. We revealed that water deficits in the soil causes lipid peroxidation in these cells and the upregulation of phospholipase D, lipoxygenase, and, concomitantly, jasmonic acid (JA) strongly accumulates in AZ tissue. Furthermore, we followed key steps in JA conjugation and signaling under stressful conditions by monitoring the level and tissue localization of enzyme providing JA derivatives (JASMONATE RESISTANT1) and the JA receptor (CORONATINE INSENSITIVE1). Collectively, drought-triggered AZ activation during the process of flower abscission is closely associated with the lipid modifications, leading to the formation of JA, its conjugation, and induction of signaling pathways.

Meta-Analysis of Common and Differential Transcriptomic Responses to Biotic and Abiotic Stresses in Arabidopsis thaliana

Environmental stresses adversely affect crop growth and yield, resulting in major losses to plants. These stresses occur simultaneously in nature, and we therefore conducted a meta-analysis in this study to identify differential and shared genes, pathways, and transcriptomic mechanisms involved in Arabidopsis response to biotic and abiotic stresses. The results showed a total of 436/21 significant up-/downregulated differentially expressed genes (DEGs) in response to biotic stresses, while 476 and 71 significant DEGs were respectively up- and downregulated in response to abiotic stresses in Arabidopsis thaliana. In addition, 21 DEGs (2.09%) were commonly regulated in response to biotic and abiotic stresses. Except for WRKY45 and ATXTH22, which were respectively up-/down- and down-/upregulated in response to biotic and abiotic stresses, other common DEGs were upregulated in response to all biotic and abiotic treatments. Moreover, the transcription factors (TFs) bHLH, MYB, and WRKY were the common TFs in response to biotic and abiotic stresses. In addition, ath-miR414 and ath-miR5658 were identified to be commonly expressed in response to both biotic and abiotic stresses. The identified common genes and pathways during biotic and abiotic stresses may provide potential candidate targets for the development of stress resistance breeding programs and for the genetic manipulation of crop plants.

Spatial differences in stoichiometry of EGR phosphatase and Microtubule-associated Stress Protein 1 control root meristem activity during drought stress

During moderate severity drought and low water potential (ψw) stress, poorly understood signaling mechanisms restrict both meristem cell division and subsequent cell expansion. We found that the Arabidopsis thaliana Clade E Growth-Regulating 2 (EGR2) protein phosphatase and Microtubule-Associated Stress Protein 1 (MASP1) differed in their stoichiometry of protein accumulation across the root meristem and had opposing effects on root meristem activity at low ψw. Ectopic MASP1 or EGR expression increased or decreased, respectively, root meristem size and root elongation during low ψw stress. This, along with the ability of phosphomimic MASP1 to overcome the EGR-mediated suppression of root meristem size and the observation that ectopic EGR expression had no effect on unstressed plants, indicated that during low ψw EGR activation and attenuation of MASP1 phosphorylation in their overlapping zone of expression determines root meristem size and activity. Ectopic EGR expression also decreased root cell size at low ψw. Conversely, both the egr1-1 egr2-1 and egr1-1 egr2-1 masp1-1 mutants had similarly increased root cell size but only egr1-1egr2-1 had increased cell division. These observations demonstrated that EGRs affect meristem activity via MASP1 but affect cell expansion via other mechanisms. Interestingly, EGR2 was highly expressed in the root cortex, a cell type important for growth regulation and environmental response.

Inter-tissue and inter-organ signaling in drought stress response and phenotyping of drought tolerance

Plant response to drought stress includes systems for intracellular regulation of gene expression and signaling, as well as inter-tissue and inter-organ signaling, which helps entire plants acquire stress resistance. Plants sense water-deficit conditions both via the stomata of leaves and roots, and transfer water-deficit signals from roots to shoots via inter-organ signaling. Abscisic acid is an important phytohormone involved in the drought stress response and adaptation, and is synthesized mainly in vascular tissues and guard cells of leaves. In leaves, stress-induced abscisic acid is distributed to various tissues by transporters, which activates stomatal closure and expression of stress-related genes to acquire drought stress resistance. Moreover, the stepwise stress response at the whole-plant level is important for proper understanding of the physiological response to drought conditions. Drought stress is sensed by multiple types of sensors as molecular patterns of abiotic stress signals, which are transmitted via separate parallel signaling networks to induce downstream responses, including stomatal closure and synthesis of stress-related proteins and metabolites. Peptide molecules play important roles in the inter-organ signaling of dehydration from roots to shoots, as well as signaling of osmotic changes and reactive oxygen species/Ca²⁺ . In this review, we have summarized recent advances in research on complex plant drought stress responses, focusing on inter-tissue signaling in leaves and inter-organ signaling from roots to shoots. We have discussed the mechanisms via which drought stress adaptations and resistance are acquired at the whole-plant level, and have proposed the importance of quantitative phenotyping for measuring plant growth under drought conditions.

Genetic Dissection of Mature Root Characteristics by Genome-Wide Association Studies in Rapeseed (Brassica napus L.)

Roots are complicated quantitative characteristics that play an essential role in absorbing water and nutrients. To uncover the genetic variations for root-related traits in rapeseed, twelve mature root traits of a Brassica napus association panel were investigated in the field within three environments. All traits showed significant phenotypic variation among genotypes, with heritabilities ranging from 55.18% to 79.68%. Genome-wide association studies (GWAS) using 20,131 SNPs discovered 172 marker-trait associations, including 103 significant SNPs (-log10 (p) > 4.30) that explained 5.24-20.31% of the phenotypic variance. With the linkage disequilibrium r² > 0.2, these significant associations were binned into 40 quantitative trait loci (QTL) clusters. Among them, 14 important QTL clusters were discovered in two environments and/or with phenotypic contributions greater than 10%. By analyzing the genomic regions within 100 kb upstream and downstream of the peak SNPs within the 14 loci, 334 annotated genes were found. Among these, 32 genes were potentially associated with root development according to their expression analysis. Furthermore, the protein interaction network using the 334 annotated genes gave nine genes involved in a substantial number of interactions, including a key gene associated with root development, BnaC09g36350D. This research provides the groundwork for deciphering B. napus' genetic variations and improving its root system architecture.

Global Genetic Heterogeneity in Adaptive Traits

Understanding the genetic architecture of complex traits is a major objective in biology. The standard approach for doing so is genome-wide association studies (GWAS), which aim to identify genetic polymorphisms responsible for variation in traits of interest. In human genetics, consistency across studies is commonly used as an indicator of reliability. However, if traits are involved in adaptation to the local environment, we do not necessarily expect reproducibility. On the contrary, results may depend on where you sample, and sampling across a wide range of environments may decrease the power of GWAS because of increased genetic heterogeneity. In this study, we examine how sampling affects GWAS in the model plant species Arabidopsis thaliana. We show that traits like flowering time are indeed influenced by distinct genetic effects in local populations. Furthermore, using gene expression as a molecular phenotype, we show that some genes are globally affected by shared variants, whereas others are affected by variants specific to subpopulations. Remarkably, the former are essentially all cis-regulated, whereas the latter are predominately affected by trans-acting variants. Our result illustrate that conclusions about genetic architecture can be extremely sensitive to sampling and population structure.

Genome-wide association studies: assessing trait characteristics in model and crop plants

GWAS involves testing genetic variants across the genomes of many individuals of a population to identify genotype–phenotype association. It was initially developed and has proven highly successful in human disease genetics. In plants genome-wide association studies (GWAS) initially focused on single feature polymorphism and recombination and linkage disequilibrium but has now been embraced by a plethora of different disciplines with several thousand studies being published in model and crop species within the last decade or so. Here we will provide a comprehensive review of these studies providing cases studies on biotic resistance, abiotic tolerance, yield associated traits, and metabolic composition. We also detail current strategies of candidate gene validation as well as the functional study of haplotypes. Furthermore, we provide a critical evaluation of the GWAS strategy and its alternatives as well as future perspectives that are emerging with the emergence of pan-genomic datasets.

Using Network-Based Machine Learning to Predict Transcription Factors Involved in Drought Resistance

Gene regulatory networks underpin stress response pathways in plants. However, parsing these networks to prioritize key genes underlying a particular trait is challenging. Here, we have built the Gene Regulation and Association Network (GRAiN) of rice (Oryza sativa). GRAiN is an interactive query-based web-platform that allows users to study functional relationships between transcription factors (TFs) and genetic modules underlying abiotic-stress responses. We built GRAiN by applying a combination of different network inference algorithms to publicly available gene expression data. We propose a supervised machine learning framework that complements GRAiN in prioritizing genes that regulate stress signal transduction and modulate gene expression under drought conditions. Our framework converts intricate network connectivity patterns of 2160 TFs into a single drought score. We observed that TFs with the highest drought scores define the functional, structural, and evolutionary characteristics of drought resistance in rice. Our approach accurately predicted the function of OsbHLH148 TF, which we validated using in vitro protein-DNA binding assays and mRNA sequencing loss-of-function mutants grown under control and drought stress conditions. Our network and the complementary machine learning strategy lends itself to predicting key regulatory genes underlying other agricultural traits and will assist in the genetic engineering of desirable rice varieties.

Distinct cellular strategies determine sensitivity to mild drought of Arabidopsis natural accessions

The worldwide distribution of Arabidopsis (Arabidopsis thaliana) accessions imposes different types of evolutionary pressures, which contributes to various responses of these accessions to environmental stresses. Responses to drought stress have mostly been studied in the Columbia accession, which is predominantly used in plant research. However, the reactions to drought stress are complex and our understanding of the responses that contribute to maintaining plant growth during mild drought (MD) is very limited. Here, we studied the mechanisms with which natural accessions react to MD at a physiological and molecular level during early leaf development. We documented variations in MD responses among natural accessions and used transcriptome sequencing of a drought-sensitive accession, ICE163, and a drought-insensitive accession, Yeg-1, to gain insights into the mechanisms underlying this discrepancy. This revealed that ICE163 preferentially induces jasmonate- and anthocyanin-related pathways, which are beneficial in biotic stress defense, whereas Yeg-1 has a more pronounced activation of abscisic acid signaling, the classical abiotic stress response. Related physiological traits, including the content of proline, anthocyanins, and reactive oxygen species, stomatal closure, and cellular leaf parameters, were investigated and linked to the transcriptional responses. We can conclude that most of these processes constitute general drought response mechanisms that are regulated similarly in drought-insensitive and -sensitive accessions. However, the capacity to close stomata and maintain cell expansion under MD appeared to be major factors that allow to better sustain leaf growth under MD.

Phloem exudate metabolic content reflects the response to water-deficit stress in pea plants (Pisum sativum L.)

Drought stress impacts the quality and yield of Pisum sativum. Here, we show how short periods of limited water availability during the vegetative stage of pea alters phloem sap content and how these changes are connected to strategies used by plants to cope with water deficit. We have investigated the metabolic content of phloem sap exudates and explored how this reflects P. sativum physiological and developmental responses to drought. Our data show that drought is accompanied by phloem-mediated redirection of the components that are necessary for cellular respiration and the proper maintenance of carbon/nitrogen balance during stress. The metabolic content of phloem sap reveals a shift from anabolic to catabolic processes as well as the developmental plasticity of P. sativum plants subjected to drought. Our study underlines the importance of phloem-mediated transport for plant adaptation to unfavourable environmental conditions. We also show that phloem exudate analysis can be used as a useful proxy to study stress responses in plants. We propose that the decrease in oleic acid content within phloem sap could be considered as a potential marker of early signalling events mediating drought response.

Proximal Hyperspectral Imaging Detects Diurnal and Drought-Induced Changes in Maize Physiology

Hyperspectral imaging is a promising tool for non-destructive phenotyping of plant physiological traits, which has been transferred from remote to proximal sensing applications, and from manual laboratory setups to automated plant phenotyping platforms. Due to the higher resolution in proximal sensing, illumination variation and plant geometry result in increased non-biological variation in plant spectra that may mask subtle biological differences. Here, a better understanding of spectral measurements for proximal sensing and their application to study drought, developmental and diurnal responses was acquired in a drought case study of maize grown in a greenhouse phenotyping platform with a hyperspectral imaging setup. The use of brightness classification to reduce the illumination-induced non-biological variation is demonstrated, and allowed the detection of diurnal, developmental and early drought-induced changes in maize reflectance and physiology. Diurnal changes in transpiration rate and vapor pressure deficit were significantly correlated with red and red-edge reflectance. Drought-induced changes in effective quantum yield and water potential were accurately predicted using partial least squares regression and the newly developed Water Potential Index 2, respectively. The prediction accuracy of hyperspectral indices and partial least squares regression were similar, as long as a strong relationship between the physiological trait and reflectance was present. This demonstrates that current hyperspectral processing approaches can be used in automated plant phenotyping platforms to monitor physiological traits with a high temporal resolution.

Induction of Acquired Tolerance Through Gradual Progression of Drought Is the Key for Maintenance of Spikelet Fertility and Yield in Rice Under Semi-irrigated Aerobic Conditions

Plants have evolved several adaptive mechanisms to cope with water-limited conditions. While most of them are through constitutive traits, certain "acquired tolerance" traits also provide significant improvement in drought adaptation. Most abiotic stresses, especially drought, show a gradual progression of stress and hence provide an opportunity to upregulate specific protective mechanisms collectively referred to as "acquired tolerance" traits. Here, we demonstrate a significant genetic variability in acquired tolerance traits among rice germplasm accessions after standardizing a novel gradual stress progress protocol. Two contrasting genotypes, BPT 5204 (drought susceptible) and AC 39000 (tolerant), were used to standardize methodology for capturing acquired tolerance traits at seedling phase. Seedlings exposed to gradual progression of stress showed higher recovery with low free radical accumulation in both the genotypes compared to rapid stress. Further, the gradual stress progression protocol was used to examine the role of acquired tolerance at flowering phase using a set of 17 diverse rice genotypes. Significant diversity in free radical production and scavenging was observed among these genotypes. Association of these parameters with yield attributes showed that genotypes that managed free radical levels in cells were able to maintain high spikelet fertility and hence yield under stress. This study, besides emphasizing the importance of acquired tolerance, explains a high throughput phenotyping approach that significantly overcomes methodological constraints in assessing genetic variability in this important drought adaptive mechanism.

Using Network-Based Machine Learning to Predict Transcription Factors Involved in Drought Resistance

Gene regulatory networks underpin stress response pathways in plants. However, parsing these networks to prioritize key genes underlying a particular trait is challenging. Here, we have built the Gene Regulation and Association Network (GRAiN) of rice (Oryza sativa). GRAiN is an interactive query-based web-platform that allows users to study functional relationships between transcription factors (TFs) and genetic modules underlying abiotic-stress responses. We built GRAiN by applying a combination of different network inference algorithms to publicly available gene expression data. We propose a supervised machine learning framework that complements GRAiN in prioritizing genes that regulate stress signal transduction and modulate gene expression under drought conditions. Our framework converts intricate network connectivity patterns of 2160 TFs into a single drought score. We observed that TFs with the highest drought scores define the functional, structural, and evolutionary characteristics of drought resistance in rice. Our approach accurately predicted the function of OsbHLH148 TF, which we validated using in vitro protein-DNA binding assays and mRNA sequencing loss-of-function mutants grown under control and drought stress conditions. Our network and the complementary machine learning strategy lends itself to predicting key regulatory genes underlying other agricultural traits and will assist in the genetic engineering of desirable rice varieties.

Leaf shape in Populus tremula is a complex, omnigenic trait

Leaf shape is a defining feature of how we recognize and classify plant species. Although there is extensive variation in leaf shape within many species, few studies have disentangled the underlying genetic architecture. We characterized the genetic architecture of leaf shape variation in Eurasian aspen (Populus tremula L.) by performing genome-wide association study (GWAS) for physiognomy traits. To ascertain the roles of identified GWAS candidate genes within the leaf development transcriptional program, we generated RNA-Seq data that we used to perform gene co-expression network analyses from a developmental series, which is publicly available within the PlantGenIE resource. We additionally used existing gene expression measurements across the population to analyze GWAS candidate genes in the context of a population-wide co-expression network and to identify genes that were differentially expressed between groups of individuals with contrasting leaf shapes. These data were integrated with expression GWAS (eQTL) results to define a set of candidate genes associated with leaf shape variation. Our results identified no clear adaptive link to leaf shape variation and indicate that leaf shape traits are genetically complex, likely determined by numerous small-effect variations in gene expression. Genes associated with shape variation were peripheral within the population-wide co-expression network, were not highly connected within the leaf development co-expression network, and exhibited signatures of relaxed selection. As such, our results are consistent with the omnigenic model.

Natural variation in Arabidopsis thaliana rosette area unveils new genes involved in plant development

Growth is a complex trait influenced by multiple genes that act at different moments during the development of an organism. This makes it difficult to spot its underlying genetic mechanisms. Since plant growth is intimately related to the effective leaf surface area (ELSA), identifying genes controlling this trait will shed light on our understanding of plant growth. To find new genes with a significant contribution to plant growth, here we used the natural variation in Arabidopsis thaliana to perform a genome-wide association study of ELSA. To do this, the projected rosette area of 710 worldwide distributed natural accessions was measured and analyzed using the genome-wide efficient mixed model association algorithm. From this analysis, ten genes were identified having SNPs with a significant association with ELSA. To validate the implication of these genes into A. thaliana growth, six of them were further studied by phenotyping knock-out mutant plants. It was observed that rem1.2, orc1a, ppd1, and mcm4 mutants showed different degrees of reduction in rosette size, thus confirming the role of these genes in plant growth. Our study identified genes already known to be involved in plant growth but also assigned this role, for the first time, to other genes.

The Evolutionary Forces Shaping Cis- and Trans-Regulation of Gene Expression within a Population of Outcrossing Plants

Understanding the persistence of genetic variation within populations has long been a goal of evolutionary biology. One promising route toward achieving this goal is using population genetic approaches to describe how selection acts on the loci associated with trait variation. Gene expression provides a model trait for addressing the challenge of the maintenance of variation because it can be measured genome-wide without information about how gene expression affects traits. Previous work has shown that loci affecting the expression of nearby genes (local or cis-eQTLs) are under negative selection, but we lack a clear understanding of the selective forces acting on variants that affect the expression of genes in trans. Here, we identify loci that affect gene expression in trans using genomic and transcriptomic data from one population of the obligately outcrossing plant, Capsella grandiflora. The allele frequencies of trans-eQTLs are consistent with stronger negative selection acting on trans-eQTLs than cis-eQTLs, and stronger negative selection acting on trans-eQTLs associated with the expression of multiple genes. However, despite this general pattern, we still observe the presence of a trans-eQTL at intermediate frequency that affects the expression of a large number of genes in the same coexpression module. Overall, our work highlights the different selective pressures shaping variation in cis- and trans-regulation.

Plant growth under suboptimal water conditions: early responses and methods to study them

Drought stress forms a major environmental constraint during the life cycle of plants, often decreasing plant yield and in extreme cases threatening survival. The molecular and physiological responses induced by drought have been the topic of extensive research during the past decades. Because soil-based approaches to studying drought responses are often challenging due to low throughput and insufficient control of the conditions, osmotic stress assays in plates were developed to mimic drought. Addition of compounds such as polyethylene glycol, mannitol, sorbitol, or NaCl to controlled growth media has become increasingly popular since it offers the advantage of accurate control of stress level and onset. These osmotic stress assays enabled the discovery of very early stress responses, occurring within seconds or minutes following osmotic stress exposure. In this review, we construct a detailed timeline of early responses to osmotic stress, with a focus on how they initiate plant growth arrest. We further discuss the specific responses triggered by different types and severities of osmotic stress. Finally, we compare short-term plant responses under osmotic stress versus in-soil drought and discuss the advantages, disadvantages, and future of these plate-based proxies for drought.

Steady-State Levels of Cytokinins and Their Derivatives May Serve as a Unique Classifier of Arabidopsis Ecotypes

We determined steady-state (basal) endogenous levels of three plant hormones (abscisic acid, cytokinins and indole-3-acetic acid) in a collection of thirty different ecotypes of Arabidopsis that represent a broad genetic variability within this species. Hormone contents were analysed separately in plant shoots and roots after 21 days of cultivation on agar plates in a climate-controlled chamber. Using advanced statistical and machine learning methods, we tested if basal hormonal levels can be considered a unique ecotype-specific classifier. We also explored possible relationships between hormone levels and the prevalent environmental conditions in the site of origin for each ecotype. We found significant variations in basal hormonal levels and their ratios in both root and shoot among the ecotypes. We showed the prominent position of cytokinins (CK) among the other hormones. We found the content of CK and CK metabolites to be a reliable ecotype-specific identifier. Correlation with the mean temperature at the site of origin and the large variation in basal hormonal levels suggest that the high variability may potentially be in response to environmental factors. This study provides a starting point for ecotype-specific genetic maps of the CK metabolic and signalling network to explore its contribution to the adaptation of plants to local environmental conditions.

AraPheno and the AraGWAS Catalog 2020: a major database update including RNA-Seq and knockout mutation data for Arabidopsis thaliana

Genome-wide association studies (GWAS) are integral for studying genotype-phenotype relationships and gaining a deeper understanding of the genetic architecture underlying trait variation. A plethora of genetic associations between distinct loci and various traits have been successfully discovered and published for the model plant Arabidopsis thaliana. This success and the free availability of full genomes and phenotypic data for more than 1,000 different natural inbred lines led to the development of several data repositories. AraPheno (https://arapheno.1001genomes.org) serves as a central repository of population-scale phenotypes in A. thaliana, while the AraGWAS Catalog (https://aragwas.1001genomes.org) provides a publicly available, manually curated and standardized collection of marker-trait associations for all available phenotypes from AraPheno. In this major update, we introduce the next generation of both platforms, including new data, features and tools. We included novel results on associations between knockout-mutations and all AraPheno traits. Furthermore, AraPheno has been extended to display RNA-Seq data for hundreds of accessions, providing expression information for over 28 000 genes for these accessions. All data, including the imputed genotype matrix used for GWAS, are easily downloadable via the respective databases.

Leaf Production and Expansion: A Generalized Response to Drought Stresses from Cells to Whole Leaf Biomass-A Case Study in the Tomato Compound Leaf

It is clearly established that there is not a unique response to soil water deficit but that there are as many responses as soil water deficit characteristics: Drought intensity, drought duration, and drought position during plant cycle. For a same soil water deficit, responses can also differ on plant genotype within a same species. In spite of this variability, at least for leaf production and expansion processes, robust tendencies can be extracted from the literature when similar watering regimes are compared. Here, we present response curves and multi-scale dynamics analyses established on tomato plants exposed to different soil water deficit treatments. Results reinforce the trends already observed for other species: Reduction in plant leaf biomass under water stress was due to reduction in individual leaf biomass and areas whereas leaf production and specific leaf area were not affected. The dynamics of leaf expansion was modified both at the leaf and cell scales. Cell division and expansion were reduced by drought treatments as well as the endoreduplication process. Combining response curves analyses together with dynamic analyses of tomato compound leaf growth at different scales not only corroborate results on simple leaf responses to drought but also increases our knowledge on the cellular mechanisms behind leaf growth plasticity.

Analysis of Drought Tolerance and Associated Traits in Upland Cotton at the Seedling Stage

(1) Background: Upland cotton (Gossypium hirsutum L.) is the most important natural fiber worldwide, and it is extensively planted and plentifully used in the textile industry. Major cotton planting regions are frequently affected by abiotic stress, especially drought stress. Drought resistance is a complex, quantitative trait. A genome-wide association study (GWAS) constitutes an efficient method for dissecting the genetic architecture of complex traits. In this study, the drought resistance of a population of 316 upland cotton accessions was studied via GWAS. (2) Methods: GWAS methodology was employed to identify relationships between molecular markers or candidate genes and phenotypes of interest. (3) Results: A total of 8, 3, and 6 SNPs were associated with the euphylla wilting score (EWS), cotyledon wilting score (CWS), and leaf temperature (LT), respectively, based on a general linear model and a factored spectrally transformed linear mixed model. For these traits, 7 QTLs were found, of which 2 each were located on chromosomes A05, A11, and D03, and of which 1 was located on chromosome A01. Importantly, in the candidate regions WRKY70, GhCIPK6, SnRK2.6, and NET1A, which are involved in the response to abscisic acid (ABA), the mitogen-activated protein kinase (MAPK) signaling pathway and the calcium transduction pathway were identified in upland cotton at the seedling stage under drought stress according to annotation information and linkage disequilibrium (LD) block analysis. Moreover, RNA sequencing analysis showed that WRKY70, GhCIPK6, SnRK2.6, and NET1A were induced by drought stress, and the expression of these genes was significantly different between normal and drought stress conditions. (4) Conclusions: The present study should provide some genomic resources for drought resistance in upland cotton. Moreover, the germplasm of the different phenotypes, the detected SNPs and, the potential candidate genes will be helpful for molecular marker-assisted breeding studies about increased drought resistance in upland cotton.

PropaNet: Time-Varying Condition-Specific Transcriptional Network Construction by Network Propagation

Transcription factor (TF) has a significant influence on the state of a cell by regulating multiple down-stream genes. Thus, experimental and computational biologists have made great efforts to construct TF gene networks for regulatory interactions between TFs and their target genes. Now, an important research question is how to utilize TF networks to investigate the response of a plant to stress at the transcription control level using time-series transcriptome data. In this article, we present a new computational network, PropaNet, to investigate dynamics of TF networks from time-series transcriptome data using two state-of-the-art network analysis techniques, influence maximization and network propagation. PropaNet uses the influence maximization technique to produce a ranked list of TFs, in the order of TF that explains differentially expressed genes (DEGs) better at each time point. Then, a network propagation technique is used to select a group of TFs that explains DEGs best as a whole. For the analysis of Arabidopsis time series datasets from AtGenExpress, we used PlantRegMap as a template TF network and performed PropaNet analysis to investigate transcriptional dynamics of Arabidopsis under cold and heat stress. The time varying TF networks showed that Arabidopsis responded to cold and heat stress quite differently. For cold stress, bHLH and bZIP type TFs were the first responding TFs and the cold signal influenced histone variants, various genes involved in cell architecture, osmosis and restructuring of cells. However, the consequences of plants under heat stress were up-regulation of genes related to accelerating differentiation and starting re-differentiation. In terms of energy metabolism, plants under heat stress show elevated metabolic process and resulting in an exhausted status. We believe that PropaNet will be useful for the construction of condition-specific time-varying TF network for time-series data analysis in response to stress. PropaNet is available at http://biohealth.snu.ac.kr/software/PropaNet.

Distinctive phytohormonal and metabolic profiles of Arabidopsis thaliana and Eutrema salsugineum under similar soil drying

Arabidopsis and Eutrema show similar stomatal sensitivity to drying soil. In Arabidopsis, larger metabolic adjustments than in Eutrema occurred, with considerable differences in the phytohormonal responses of the two species. Although plants respond to soil drying via a series of concurrent physiological and molecular events, drought tolerance differs greatly within the plant kingdom. While Eutrema salsugineum (formerly Thellungiella salsuginea) is regarded as more stress tolerant than its close relative Arabidopsis thaliana, their responses to soil water deficit have not previously been directly compared. To ensure a similar rate of soil drying for the two species, daily soil water depletion was controlled to 5-10% of the soil water content. While partial stomatal closure occurred earlier in Arabidopsis (Day 4) than Eutrema (from Day 6 onwards), thereafter both species showed similar stomatal sensitivity to drying soil. However, both targeted and untargeted metabolite analysis revealed greater response to drought in Arabidopsis than Eutrema. Early peaks in foliar phytohormone concentrations and different sugar profiles between species were accompanied by opposing patterns in the bioactive cytokinin profiles. Untargeted analysis showed greater metabolic adjustment in Arabidopsis with more statistically significant changes in both early and severe drought stress. The distinct metabolic responses of each species during early drought, which occurred prior to leaf water status declining, seemed independent of later stomatal closure in response to drought. The two species also showed distinct water usage, with earlier reduction in water consumption in Eutrema (Day 3) than Arabidopsis (Day 6), likely reflecting temporal differences in growth responses. We propose Arabidopsis as a promising model to evaluate the mechanisms responsible for stress-induced growth inhibition under the mild/moderate soil drying that crop plants are typically exposed to.

The complex genetic architecture of shoot growth natural variation in Arabidopsis thaliana

One of the main outcomes of quantitative genetics approaches to natural variation is to reveal the genetic architecture underlying the phenotypic space. Complex genetic architectures are described as including numerous loci (or alleles) with small-effect and/or low-frequency in the populations, interactions with the genetic background, environment or age. Linkage or association mapping strategies will be more or less sensitive to this complexity, so that we still have an unclear picture of its extent. By combining high-throughput phenotyping under two environmental conditions with classical QTL mapping approaches in multiple Arabidopsis thaliana segregating populations as well as advanced near isogenic lines construction and survey, we have attempted to improve our understanding of quantitative phenotypic variation. Integrative traits such as those related to vegetative growth used in this work (highlighting either cumulative growth, growth rate or morphology) all showed complex and dynamic genetic architecture with respect to the segregating population and condition. The more resolutive our mapping approach, the more complexity we uncover, with several instances of QTLs visible in near isogenic lines but not detected with the initial QTL mapping, indicating that our phenotyping accuracy was less limiting than the mapping resolution with respect to the underlying genetic architecture. In an ultimate approach to resolve this complexity, we intensified our phenotyping effort to target specifically a 3Mb-region known to segregate for a major quantitative trait gene, using a series of selected lines recombined every 100kb. We discovered that at least 3 other independent QTLs had remained hidden in this region, some with trait- or condition-specific effects, or opposite allelic effects. If we were to extrapolate the figures obtained on this specific region in this particular cross to the genome- and species-scale, we would predict hundreds of causative loci of detectable phenotypic effect controlling these growth-related phenotypes.

Low Water Potential and At14a-Like1 (AFL1) Effects on Endocytosis and Actin Filament Organization

At14a-Like1 (AFL1) is a stress-induced protein of unknown function that promotes growth during low water potential stress and drought. Previous analysis indicated that AFL1 may have functions related to endocytosis and regulation of actin filament organization, processes for which the effects of low water potential are little known. We found that low water potential led to a decrease in endocytosis, as measured by uptake of the membrane-impermeable dye FM4-64. Ectopic expression of AFL1 reversed the decrease in FM4-64 uptake seen in wild type, while reduced AFL1 expression led to further inhibition of FM4-64 uptake. Increased AFL1 also made FM4-64 uptake less sensitive to the actin filament disruptor Latrunculin B (LatB). LatB decreased AFL1-Clathrin Light Chain colocalization, further indicating that effects of AFL1 on endocytosis may be related to actin filament organization or stability. Consistent with this hypothesis, ectopic AFL1 expression made actin filaments less sensitive to disruption by LatB or Cytochalasin D and led to increased actin filament skewness and decreased occupancy, indicative of more bundled actin filaments. This latter effect could be partially mimicked by the actin filament stabilizer Jasplakinolide (JASP). However, AFL1 did not substantially inhibit actin filament dynamics, indicating that AFL1 acts via a different mechanism than JASP-induced stabilization. AFL1 partially colocalized with actin filaments but not with microtubules, further indicating actin-filament-related function of AFL1. These data provide insight into endocytosis and actin filament responses to low water potential stress and demonstrate an involvement of AFL1 in these key cellular processes.

A new regulator of seed size control in Arabidopsis identified by a genome-wide association study

Organ size in plants is controlled by the interaction between genotype and the environment. Seed size, an important agronomic trait, largely determines yield and is an important focus of research. However, the genetic components underpinning natural variation of seed size in undomesticated species remain largely unidentified. Here we report a genome-wide association study (GWAS) of seed size in Arabidopsis thaliana, which identified 38 significantly associated loci, including one locus associated with CYCB1;4. Natural variations in CYCB1;4, which encodes a cyclin protein involved in the cell cycle, significantly influence seed size in A. thaliana. Transgenic plants with enhanced CYCB1;4 expression show normal development, exhibit increased seed size as a result of an accelerated cell cycle progression, and tend to produce higher yields. By contrast, cycb1;4 mutants have smaller seeds, and the effect is especially pronounced in a large-seed accession. The temporal and spatial expression pattern of CYCB1;4 suggests that this gene may function in both maternal tissues and zygotic tissues to coordinate the final size of seeds. Taken together, our results provide genetic insights into natural variation in seed size in Arabidopsis. Moreover, CYCB1;4 homologs in other crops could have great potential as targets for efforts aimed at yield improvement.

Linking genes with ecological strategies in Arabidopsis thaliana

Arabidopsis thaliana is the most prominent model system in plant molecular biology and genetics. Although its ecology was initially neglected, collections of various genotypes revealed a complex population structure, with high levels of genetic diversity and substantial levels of phenotypic variation. This helped identify the genes and gene pathways mediating phenotypic change. Population genetics studies further demonstrated that this variation generally contributes to local adaptation. Here, we review evidence showing that traits affecting plant life history, growth rate, and stress reactions are not only locally adapted, they also often co-vary. Co-variation between these traits indicates that they evolve as trait syndromes, and reveals the ecological diversification that took place within A. thaliana. We argue that examining traits and the gene that control them within the context of global summary schemes that describe major ecological strategies will contribute to resolve important questions in both molecular biology and ecology.

Heliaphen, an Outdoor High-Throughput Phenotyping Platform for Genetic Studies and Crop Modeling

Heliaphen is an outdoor platform designed for high-throughput phenotyping. It allows the automated management of drought scenarios and monitoring of plants throughout their lifecycles. A robot moving between plants growing in 15-L pots monitors the plant water status and phenotypes the leaf or whole-plant morphology. From these measurements, we can compute more complex traits, such as leaf expansion (LE) or transpiration rate (TR) in response to water deficit. Here, we illustrate the capabilities of the platform with two practical cases in sunflower (Helianthus annuus): a genetic and genomic study of the response of yield-related traits to drought, and a modeling study using measured parameters as inputs for a crop simulation. For the genetic study, classical measurements of thousand-kernel weight (TKW) were performed on a biparental population under automatically managed drought stress and control conditions. These data were used for an association study, which identified five genetic markers of the TKW drought response. A complementary transcriptomic analysis identified candidate genes associated with these markers that were differentially expressed in the parental backgrounds in drought conditions. For the simulation study, we used a crop simulation model to predict the impact on crop yield of two traits measured on the platform (LE and TR) for a large number of environments. We conducted simulations in 42 contrasting locations across Europe using 21 years of climate data. We defined the pattern of abiotic stresses occurring at the continental scale and identified ideotypes (i.e., genotypes with specific trait values) that are more adapted to specific environment types. This study exemplifies how phenotyping platforms can assist the identification of the genetic architecture controlling complex response traits and facilitate the estimation of ecophysiological model parameters to define ideotypes adapted to different environmental conditions.

MeABL5, an ABA Insensitive 5-Like Basic Leucine Zipper Transcription Factor, Positively Regulates MeCWINV3 in Cassava (Manihot esculenta Crantz)

The basic leucine zipper (bZIP) transcription factor family plays crucial roles in multiple biological processes, especially stress responses. Cassava (Manihot esculenta Crantz) is an important tropical crop with a strong tolerance to environmental stresses such as drought, heat, and low-fertility environments. Currently, limited information is available regarding the functional identification of bZIP transcription factors in response to abiotic stress in cassava. Herein, a gene encoding an ABA Insensitive 5 (ABI5)-like transcription factor, designated as MeABL5, was identified in cassava. Sequence and phylogenetic analysis showed that MeABL5 is a cassava bZIP transcription factor that is not included in the previously identified cassava bZIP family members, belongs to subfamily A, and has high sequence similarity to ABI5-like proteins. Subcellular localization and transactivation assays revealed that MeABL5 was a nuclear-localized protein and possessed transactivation activity. Furthermore, MeABL5 was able to specifically interact with the ABRE cis-element in the promoter of the cassava major cell wall invertase gene, MeCWINV3, in vitro and in vivo. MeABL5 and MeCWINV3 exhibited similar expression patterns in various organs or tissues and under abiotic stress in cassava. The expressions of MeABL5 and MeCWINV3 within cassava plantlets were both induced by exogenous abscisic acid (ABA), gibberellic acid (GA3), methyl jasmonate (MeJA), and heat. Overexpression of MeABL5 increased the activity of the MeCWINV3 gene, and the up-regulated expressions of MeCWINV3 were significantly activated under ABA-, salicylic acid (SA)-, and MeJA-induced conditions. Overall, these results suggest that MeABL5 is a positive regulator of MeCWINV3 and might participate in the robust resistance of cassava in response to abiotic stress. This study also provides a foundation for further research on ABA-mediated and stress-related signaling pathways in cassava.

Sodium-Related Adaptations to Drought: New Insights From the Xerophyte Plant Zygophyllum xanthoxylum

Understanding the unusual physiological mechanisms that enable drought tolerance in xerophytes will be of considerable benefit because of the potential to identify novel and key genetic elements for future crop improvements. These plants are interesting because they are well-adapted for life in arid zones; Zygophyllum xanthoxylum, for example, is a typical xerophytic shrub that inhabits central Asian deserts, accumulating substantial levels of sodium (Na+) in its succulent leaves while growing in soils that contain very low levels of this ion. The physiological importance of this unusual trait to drought adaptations remains poorly understood, however. Thus, 2-week-old Z. xanthoxylum plants were treated with 50 mM NaCl (Na) for 7 days in this study in order to investigate their drought tolerance, leaf osmotic potential (Ψs) related parameters, anatomical characteristics, and transpiration traits. The results demonstrated that NaCl treatment significantly enhanced both the survivability and durability of Z. xanthoxylum plants under extreme drought conditions. The bulk of the Na+ ions encapsulated in plants was overwhelmingly allocated to leaves rather than roots or stems under drought conditions; thus, compared to the control, significantly more Na+ compared to other solutes such as K+, Ca2+, Cl-, sugars, and proline accumulated in the leaves of NaCl-treated plants and led to a marked decrease (31%) in leaf Ψs. In addition, the accumulation of Na+ ions also resulted in mesophyll cell enlargement and leaf succulence, enabling the additional storage of water; Na+ ions also reduced the rate of water loss by decreasing stomatal density and down-regulating stomatal aperture size. The results of this study demonstrate that Z. xanthoxylum has evolved a notable ability to utilize Na+ ions to lower Ψs, swell its leaves, and decrease stomatal aperture sizes, in order to enable the additional uptake and storage of water and mitigate losses. These distinctive drought adaption characteristics mean that the xerophytic plant Z. xanthoxylum presents a fascinating case study for the potential identification of important and novel genetic elements that could improve crops. This report provides insights on the eco-physiological role of sodium accumulation in xerophytes adapted to extremely arid habitats.

Allele-specific expression and genetic determinants of transcriptomic variations in response to mild water deficit in tomato

Characterizing the natural diversity of gene expression across environments is an important step in understanding how genotype-by-environment interactions shape phenotypes. Here, we analyzed the impact of water deficit onto gene expression levels in tomato at the genome-wide scale. We sequenced the transcriptome of growing leaves and fruit pericarps at cell expansion stage in a cherry and a large fruited accession and their F₁ hybrid grown under two watering regimes. Gene expression levels were steadily affected by the genotype and the watering regime. Whereas phenotypes showed mostly additive inheritance, ~80% of the genes displayed non-additive inheritance. By comparing allele-specific expression (ASE) in the F₁ hybrid to the allelic expression in both parental lines, respectively, 3005 genes in leaf and 2857 genes in fruit deviated from 1:1 ratio independently of the watering regime. Among these genes, ~55% were controlled by cis factors, ~25% by trans factors and ~20% by a combination of both types of factors. A total of 328 genes in leaf and 113 in fruit exhibited significant ASE-by-watering regime interaction, among which ~80% presented trans-by-watering regime interaction, suggesting a response to water deficit mediated through a majority of trans-acting loci in tomato. We cross-validated the expression levels of 274 transcripts in fruit and leaves of 124 recombinant inbred lines (RILs) and identified 163 expression quantitative trait loci (eQTLs) mostly confirming the divergences identified by ASE. Combining phenotypic and expression data, we observed a complex network of variation between genes encoding enzymes involved in the sugar metabolism.