Drought response in Arabidopsis displays synergistic coordination between stems and leaves.

The synergy between drought-responsive traits across different organs is crucial in the whole-plant mechanism influencing drought resilience. These organ interactions, however, are poorly understood, limiting our understanding of drought response strategies at the whole-plant level. Therefore, we need more integrative studies, especially on herbaceous species that represent many important food crops but remain underexplored in their drought response. We investigated inflorescence stems and rosette leaves of six Arabidopsis thaliana genotypes with contrasting drought tolerance, and combined anatomical observations with hydraulic measurements and gene expression studies to assess differences in drought response. The soc1ful double mutant was the most drought-tolerant genotype based on its synergistic combination of low stomatal conductance, largest stomatal safety margin, more stable leaf water potential during non-watering, reduced transcript levels of drought stress marker genes, and reduced loss of chlorophyll content in leaves, in combination with stems showing the highest embolism resistance, most pronounced lignification, and thickest intervessel pit membranes. In contrast, the most sensitive Cvi ecotype shows the opposite extreme of the same set of traits. The remaining four genotypes show variations in this drought syndrome. Our results reveal that anatomical, ecophysiological, and molecular adaptations across organs are intertwined, and multiple (differentially combined) strategies can be applied to acquire a certain level of drought tolerance.

Epigenetic mapping of the Arabidopsis metabolome reveals mediators of the epigenotype-phenotype map.

Identifying the sources of natural variation underlying metabolic differences between plants will enable a better understanding of plant metabolism and provide insights into the regulatory networks that govern plant growth and morphology. So far, however, the contribution of epigenetic variation to metabolic diversity has been largely ignored. In the present study, we utilized a panel of Arabidopsis thaliana epigenetic recombinant inbred lines (epiRILs) to assess the impact of epigenetic variation on the metabolic composition. Thirty epigenetic QTL (QTLepi) were detected, which partly overlap with QTLepi linked to growth and morphology. In an effort to identify causal candidate genes in the QTLepi regions and their putative trans-targets, we performed in silico small RNA and qPCR analyses. Differentially expressed genes were further studied by phenotypic and metabolic analyses of knockout mutants. Three genes were detected that recapitulated the detected QTLepi effects, providing evidence for epigenetic regulation in cis and in trans These results indicate that epigenetic mechanisms impact metabolic diversity, possibly via small RNAs, and thus aid in further disentangling the complex epigenotype-phenotype map.

Publisher Correction: A cheminformatics approach to characterize metabolomes in stable-isotope-labeled organisms.

In the originally published Supplementary Information for this paper, the files presented as Supplementary Tables 3, 4, and 7 were duplicates of Supplementary Tables 5, 6, and 9, respectively. All Supplementary Table files are now correct online.

A cheminformatics approach to characterize metabolomes in stable-isotope-labeled organisms.

We report a computational approach (implemented in MS-DIAL 3.0; http://prime.psc.riken.jp/) for metabolite structure characterization using fully 13C-labeled and non-labeled plants and LC-MS/MS. Our approach facilitates carbon number determination and metabolite classification for unknown molecules. Applying our method to 31 tissues from 12 plant species, we assigned 1,092 structures and 344 formulae to 3,604 carbon-determined metabolite ions, 69 of which were found to represent structures currently not listed in metabolome databases.

Meiotic crossover reduction by virus-induced gene silencing enables the efficient generation of chromosome substitution lines and reverse breeding in Arabidopsis thaliana.

Plant breeding applications exploiting meiotic mutant phenotypes (like the increase or decrease of crossover (CO) recombination) have been proposed over the last years. As recessive meiotic mutations in breeding lines may affect fertility or have other pleiotropic effects, transient silencing techniques may be preferred. Reverse breeding is a breeding technique that would benefit from the transient downregulation of CO formation. The technique is essentially the opposite of plant hybridization: a method to extract parental lines from a hybrid. The method can also be used to efficiently generate chromosome substitution lines (CSLs). For successful reverse breeding, the two homologous chromosome sets of a heterozygous plant must be divided over two haploid complements, which can be achieved by the suppression of meiotic CO recombination and the subsequent production of doubled haploid plants. Here we show the feasibility of transiently reducing CO formation using virus-induced gene silencing (VIGS) by targeting the meiotic gene MSH5 in a wild-type heterozygote of Arabidopsis thaliana. The application of VIGS (rather than using lengthy stable transformation) generates transgene-free offspring with the desired genetic composition: we obtained parental lines from a wild-type heterozygous F1 in two generations. In addition, we obtained 20 (of the 32 possible) CSLs in one experiment. Our results demonstrate that meiosis can be modulated at will in A. thaliana to generate CSLs and parental lines rapidly for hybrid breeding. Furthermore, we illustrate how the modification of meiosis using VIGS can open routes to develop efficient plant breeding strategies.

Intervessel pit membrane thickness best explains variation in embolism resistance amongst stems of Arabidopsis thaliana accessions.

BACKGROUND AND AIMS: The ability to avoid drought-induced embolisms in the xylem is one of the essential traits for plants to survive periods of water shortage. Over the past three decades, hydraulic studies have been focusing on trees, which limits our ability to understand how herbs tolerate drought. Here we investigate the embolism resistance in inflorescence stems of four Arabidopsis thaliana accessions that differ in growth form and drought response. We assess functional traits underlying the variation in embolism resistance amongst the accessions studied using detailed anatomical observations. METHODS: Vulnerability to xylem embolism was evaluated via vulnerability curves using the centrifuge technique and linked with detailed anatomical observations in stems using light microscopy and transmission electron microscopy. KEY RESULTS: The data show significant differences in stem P50, varying 2-fold from -1.58 MPa in the Cape Verde Island accession to -3.07 MPa in the woody soc1 ful double mutant. Out of all the anatomical traits measured, intervessel pit membrane thickness (TPM) best explains the differences in P50, as well as P12 and P88. The association between embolism resistance and TPM can be functionally explained by the air-seeding hypothesis. There is no evidence that the correlation between increased woodiness and increased embolism resistance is directly related to functional aspects. However, we found that increased woodiness is strongly linked to other lignification characters, explaining why mechanical stem reinforcement is indirectly related to increased embolism resistance. CONCLUSIONS: The woodier or more lignified accessions are more resistant to embolism than the herbaceous accessions, confirming the link between increased stem lignification and increased embolism resistance, as also observed in other lineages. Intervessel pit membrane thickness and, to a lesser extent, theoretical vessel implosion resistance and vessel wall thickness are the missing functional links between stem lignification and embolism resistance.

Converging phenomics and genomics to study natural variation in plant photosynthetic efficiency.

In recent years developments in plant phenomic approaches and facilities have gradually caught up with genomic approaches. An opportunity lies ahead to dissect complex, quantitative traits when both genotype and phenotype can be assessed at a high level of detail. This is especially true for the study of natural variation in photosynthetic efficiency, for which forward genetics studies have yielded only a little progress in our understanding of the genetic layout of the trait. High-throughput phenotyping, primarily from chlorophyll fluorescence imaging, should help to dissect the genetics of photosynthesis at the different levels of both plant physiology and development. Specific emphasis should be directed towards understanding the acclimation of the photosynthetic machinery in fluctuating environments, which may be crucial for the identification of genetic variation for relevant traits in food crops. Facilities should preferably be designed to accommodate phenotyping of photosynthesis-related traits in such environments. The use of forward genetics to study the genetic architecture of photosynthesis is likely to lead to the discovery of novel traits and/or genes that may be targeted in breeding or bio-engineering approaches to improve crop photosynthetic efficiency. In the near future, big data approaches will play a pivotal role in data processing and streamlining the phenotype-to-gene identification pipeline.

Genome-Wide Association Mapping Reveals That Specific and Pleiotropic Regulatory Mechanisms Fine-Tune Central Metabolism and Growth in Arabidopsis.

Central metabolism is a coordinated network that is regulated at multiple levels by resource availability and by environmental and developmental cues. Its genetic architecture has been investigated by mapping metabolite quantitative trait loci (QTL). A more direct approach is to identify enzyme activity QTL, which distinguishes between cis-QTL in structural genes encoding enzymes and regulatory trans-QTL. Using genome-wide association studies, we mapped QTL for 24 enzyme activities, nine metabolites, three structural components, and biomass in Arabidopsis thaliana We detected strong cis-QTL for five enzyme activities. A cis-QTL for UDP-glucose pyrophosphorylase activity in the UGP1 promoter is maintained through balancing selection. Variation in acid invertase activity reflects multiple evolutionary events in the promoter and coding region of VAC-INVcis-QTL were also detected for ADP-glucose pyrophosphorylase, fumarase, and phosphoglucose isomerase activity. We detected many trans-QTL, including transcription factors, E3 ligases, protein targeting components, and protein kinases, and validated some by knockout analysis. trans-QTL are more frequent but tend to have smaller individual effects than cis-QTL. We detected many colocalized QTL, including a multitrait QTL on chromosome 4 that affects six enzyme activities, three metabolites, protein, and biomass. These traits are coordinately modified by different ACCELERATED CELL DEATH6 alleles, revealing a trade-off between metabolism and defense against biotic stress.

Genetic Components of Root Architecture Remodeling in Response to Salt Stress.

Salinity of the soil is highly detrimental to plant growth. Plants respond by a redistribution of root mass between main and lateral roots, yet the genetic machinery underlying this process is still largely unknown. Here, we describe the natural variation among 347 Arabidopsis thaliana accessions in root system architecture (RSA) and identify the traits with highest natural variation in their response to salt. Salt-induced changes in RSA were associated with 100 genetic loci using genome-wide association studies. Two candidate loci associated with lateral root development were validated and further investigated. Changes in CYP79B2 expression in salt stress positively correlated with lateral root development in accessions, and cyp79b2 cyp79b3 double mutants developed fewer and shorter lateral roots under salt stress, but not in control conditions. By contrast, high HKT1 expression in the root repressed lateral root development, which could be partially rescued by addition of potassium. The collected data and multivariate analysis of multiple RSA traits, available through the Salt_NV_Root App, capture root responses to salinity. Together, our results provide a better understanding of effective RSA remodeling responses, and the genetic components involved, for plant performance in stress conditions.

Parental DNA Methylation States Are Associated with Heterosis in Epigenetic Hybrids.

Despite the importance and wide exploitation of heterosis in commercial crop breeding, the molecular mechanisms behind this phenomenon are not completely understood. Recent studies have implicated changes in DNA methylation and small RNAs in hybrid performance; however, it remains unclear whether epigenetic changes are a cause or a consequence of heterosis. Here, we analyze a large panel of over 500 Arabidopsis (Arabidopsis thaliana) epigenetic hybrid plants (epiHybrids), which we derived from near-isogenic but epigenetically divergent parents. This proof-of-principle experimental system allowed us to quantify the contribution of parental methylation differences to heterosis. We measured traits such as leaf area, growth rate, flowering time, main stem branching, rosette branching, and final plant height and observed several strong positive and negative heterotic phenotypes among the epiHybrids. Using an epigenetic quantitative trait locus mapping approach, we were able to identify specific differentially methylated regions in the parental genomes that are associated with hybrid performance. Sequencing of methylomes, transcriptomes, and genomes of selected parent-epiHybrid combinations further showed that these parental differentially methylated regions most likely mediate the remodeling of methylation and transcriptional states at specific loci in the hybrids. Taken together, our data suggest that locus-specific epigenetic divergence between the parental lines can directly or indirectly trigger heterosis in Arabidopsis hybrids independent of genetic changes. These results add to a growing body of evidence that points to epigenetic factors as one of the key determinants of hybrid performance.

AraQTL - workbench and archive for systems genetics in Arabidopsis thaliana.

Genetical genomics studies uncover genome-wide genetic interactions between genes and their transcriptional regulators. High-throughput measurement of gene expression in recombinant inbred line populations has enabled investigation of the genetic architecture of variation in gene expression. This has the potential to enrich our understanding of the molecular mechanisms affected by and underlying natural variation. Moreover, it contributes to the systems biology of natural variation, as a substantial number of experiments have resulted in a valuable amount of interconnectable phenotypic, molecular and genotypic data. A number of genetical genomics studies have been published for Arabidopsis thaliana, uncovering many expression quantitative trait loci (eQTLs). However, these complex data are not easily accessible to the plant research community, leaving most of the valuable genetic interactions unexplored as cross-analysis of these studies is a major effort. We address this problem with AraQTL (http://www.bioinformatics.nl/Ara QTL/), an easily accessible workbench and database for comparative analysis and meta-analysis of all published Arabidopsis eQTL datasets. AraQTL provides a workbench for comparing, re-using and extending upon the results of these experiments. For example, one can easily screen a physical region for specific local eQTLs that could harbour candidate genes for phenotypic QTLs, or detect gene-by-environment interactions by comparing eQTLs under different conditions.

Epigenetic basis of morphological variation and phenotypic plasticity in Arabidopsis thaliana.

Epigenetics is receiving growing attention in the plant science community. Epigenetic modifications are thought to play a particularly important role in fluctuating environments. It is hypothesized that epigenetics contributes to plant phenotypic plasticity because epigenetic modifications, in contrast to DNA sequence variation, are more likely to be reversible. The population of decrease in DNA methylation 1-2 (ddm1-2)-derived epigenetic recombinant inbred lines (epiRILs) in Arabidopsis thaliana is well suited for studying this hypothesis, as DNA methylation differences are maximized and DNA sequence variation is minimized. Here, we report on the extensive heritable epigenetic variation in plant growth and morphology in neutral and saline conditions detected among the epiRILs. Plant performance, in terms of branching and leaf area, was both reduced and enhanced by different quantitative trait loci (QTLs) in the ddm1-2 inherited epigenotypes. The variation in plasticity associated significantly with certain genomic regions in which the ddm1-2 inherited epigenotypes caused an increased sensitivity to environmental changes, probably due to impaired genetic regulation in the epiRILs. Many of the QTLs for morphology and plasticity overlapped, suggesting major pleiotropic effects. These findings indicate that epigenetics contributes substantially to variation in plant growth, morphology, and plasticity, especially under stress conditions.

Genome-Wide Association Mapping of Fertility Reduction upon Heat Stress Reveals Developmental Stage-Specific QTLs in Arabidopsis thaliana.

For crops that are grown for their fruits or seeds, elevated temperatures that occur during flowering and seed or fruit set have a stronger effect on yield than high temperatures during the vegetative stage. Even short-term exposure to heat can have a large impact on yield. In this study, we used Arabidopsis thaliana to study the effect of short-term heat exposure on flower and seed development. The impact of a single hot day (35°C) was determined in more than 250 natural accessions by measuring the lengths of the siliques along the main inflorescence. Two sensitive developmental stages were identified, one before anthesis, during male and female meiosis, and one after anthesis, during fertilization and early embryo development. In addition, we observed a correlation between flowering time and heat tolerance. Genome-wide association mapping revealed four quantitative trait loci (QTLs) strongly associated with the heat response. These QTLs were developmental stage specific, as different QTLs were detected before and after anthesis. For a number of QTLs, T-DNA insertion knockout lines could validate assigned candidate genes. Our findings show that the regulation of complex traits can be highly dependent on the developmental timing.

Genetic architecture of plant stress resistance: multi-trait genome-wide association mapping.

Plants are exposed to combinations of various biotic and abiotic stresses, but stress responses are usually investigated for single stresses only. Here, we investigated the genetic architecture underlying plant responses to 11 single stresses and several of their combinations by phenotyping 350 Arabidopsis thaliana accessions. A set of 214 000 single nucleotide polymorphisms (SNPs) was screened for marker-trait associations in genome-wide association (GWA) analyses using tailored multi-trait mixed models. Stress responses that share phytohormonal signaling pathways also share genetic architecture underlying these responses. After removing the effects of general robustness, for the 30 most significant SNPs, average quantitative trait locus (QTL) effect sizes were larger for dual stresses than for single stresses. Plants appear to deploy broad-spectrum defensive mechanisms influencing multiple traits in response to combined stresses. Association analyses identified QTLs with contrasting and with similar responses to biotic vs abiotic stresses, and below-ground vs above-ground stresses. Our approach allowed for an unprecedented comprehensive genetic analysis of how plants deal with a wide spectrum of stress conditions.

Natural variation in life history strategy of Arabidopsis thaliana determines stress responses to drought and insects of different feeding guilds.

Plants are sessile organisms and, consequently, are exposed to a plethora of stresses in their local habitat. As a result, different populations of a species are subject to different selection pressures leading to adaptation to local conditions and intraspecific divergence. The annual brassicaceous plant Arabidopsis thaliana is an attractive model for ecologists and evolutionary biologists due to the availability of a large collection of resequenced natural accessions. Accessions of A. thaliana display one of two different life cycle strategies: summer and winter annuals. We exposed a collection of 308 European Arabidopsis accessions, that have been genotyped for 250K SNPs, to a range of stresses: one abiotic stress (drought), four biotic stresses (Pieris rapae caterpillars, Plutella xylostella caterpillars, Frankliniella occidentalis thrips and Myzus persicae aphids) and two combined stresses (drought plus P. rapae and Botrytis cinerea fungus plus P. rapae). We identified heritable genetic variation for responses to the different stresses, estimated by narrow-sense heritability. We found that accessions displaying different life cycle strategies differ in their response to stresses. Winter annuals are more resistant to drought, aphids and thrips and summer annuals are more resistant to P. rapae and P. xylostella caterpillars. Summer annuals are also more resistant to the combined stresses of drought plus P. rapae and infection by the fungus Botryris cinerea plus herbivory by P. rapae. Adaptation to drought displayed a longitudinal gradient. Finally, trade-offs were recorded between the response to drought and responses to herbivory by caterpillars of the specialist herbivore P. rapae.

Genome-Wide Association Mapping and Genomic Prediction Elucidate the Genetic Architecture of Morphological Traits in Arabidopsis.

Quantitative traits in plants are controlled by a large number of genes and their interaction with the environment. To disentangle the genetic architecture of such traits, natural variation within species can be explored by studying genotype-phenotype relationships. Genome-wide association studies that link phenotypes to thousands of single nucleotide polymorphism markers are nowadays common practice for such analyses. In many cases, however, the identified individual loci cannot fully explain the heritability estimates, suggesting missing heritability. We analyzed 349 Arabidopsis accessions and found extensive variation and high heritabilities for different morphological traits. The number of significant genome-wide associations was, however, very low. The application of genomic prediction models that take into account the effects of all individual loci may greatly enhance the elucidation of the genetic architecture of quantitative traits in plants. Here, genomic prediction models revealed different genetic architectures for the morphological traits. Integrating genomic prediction and association mapping enabled the assignment of many plausible candidate genes explaining the observed variation. These genes were analyzed for functional and sequence diversity, and good indications that natural allelic variation in many of these genes contributes to phenotypic variation were obtained. For ACS11, an ethylene biosynthesis gene, haplotype differences explaining variation in the ratio of petiole and leaf length could be identified.

Genetic analysis of metabolome-phenotype interactions: from model to crop species.

The past decade has seen increased interest from the scientific community, and particularly plant biologists, in integrating metabolic approaches into research aimed at unraveling phenotypic diversity and its underlying genetic variation. Advances in plant metabolomics have enabled large-scale analyses that have identified qualitative and quantitative variation in the metabolic content of various species, and this variation has been linked to genetic factors through genetic-mapping approaches, providing a glimpse of the genetic architecture of the plant metabolome. Parallel analyses of morphological phenotypes and physiological performance characteristics have further enhanced our understanding of the complex molecular mechanisms regulating these quantitative traits. This review aims to illustrate the advantages of including assessments of phenotypic and metabolic diversity in investigations of the genetic basis of complex traits, and the value of this approach in studying agriculturally important crops. We highlight the ground-breaking work on model species and discuss recent achievements in important crop species.

Genome-wide association mapping of time-dependent growth responses to moderate drought stress in Arabidopsis.

Large areas of arable land are often confronted with irregular rainfall resulting in limited water availability for part(s) of the growing seasons, which demands research for drought tolerance of plants. Natural variation was observed for biomass accumulation upon controlled moderate drought stress in 324 natural accessions of Arabidopsis. Improved performance under drought stress was correlated with early flowering and lack of vernalization requirement, indicating overlap in the regulatory networks of flowering time and drought response or correlated responses of these traits to natural selection. In addition, plant size was negatively correlated with relative water content (RWC) independent of the absolute water content (WC), indicating a prominent role for soluble compounds. Growth in control and drought conditions was determined over time and was modelled by an exponential function. Genome-wide association (GWA) mapping of temporal plant size data and of model parameters resulted in the detection of six time-dependent quantitative trait loci (QTLs) strongly associated with drought. Most QTLs would not have been identified if plant size was determined at a single time point. Analysis of earlier reported gene expression changes upon drought enabled us to identify for each QTL the most likely candidates.

Arabidopsis semidwarfs evolved from independent mutations in GA20ox1, ortholog to green revolution dwarf alleles in rice and barley.

Understanding the genetic bases of natural variation for developmental and stress-related traits is a major goal of current plant biology. Variation in plant hormone levels and signaling might underlie such phenotypic variation occurring even within the same species. Here we report the genetic and molecular basis of semidwarf individuals found in natural Arabidopsis thaliana populations. Allelism tests demonstrate that independent loss-of-function mutations at GA locus 5 (GA5), which encodes gibberellin 20-oxidase 1 (GA20ox1) involved in the last steps of gibberellin biosynthesis, are found in different populations from southern, western, and northern Europe; central Asia; and Japan. Sequencing of GA5 identified 21 different loss-of-function alleles causing semidwarfness without any obvious general tradeoff affecting plant performance traits. GA5 shows signatures of purifying selection, whereas GA5 loss-of-function alleles can also exhibit patterns of positive selection in specific populations as shown by Fay and Wu's H statistics. These results suggest that antagonistic pleiotropy might underlie the occurrence of GA5 loss-of-function mutations in nature. Furthermore, because GA5 is the ortholog of rice SD1 and barley Sdw1/Denso green revolution genes, this study illustrates the occurrence of conserved adaptive evolution between wild A.thaliana and domesticated plants.

Genome-wide association mapping of growth dynamics detects time-specific and general quantitative trait loci.

Growth is a complex trait determined by the interplay between many genes, some of which play a role at a specific moment during development whereas others play a more general role. To identify the genetic basis of growth, natural variation in Arabidopsis rosette growth was followed in 324 accessions by a combination of top-view imaging, high-throughput image analysis, modelling of growth dynamics, and end-point fresh weight determination. Genome-wide association (GWA) mapping of the temporal growth data resulted in the detection of time-specific quantitative trait loci (QTLs), whereas mapping of model parameters resulted in another set of QTLs related to the whole growth curve. The positive correlation between projected leaf area (PLA) at different time points during the course of the experiment suggested the existence of general growth factors with a function in multiple developmental stages or with prolonged downstream effects. Many QTLs could not be identified when growth was evaluated only at a single time point. Eleven candidate genes were identified, which were annotated to be involved in the determination of cell number and size, seed germination, embryo development, developmental phase transition, or senescence. For eight of these, a mutant or overexpression phenotype related to growth has been reported, supporting the identification of true positives. In addition, the detection of QTLs without obvious candidate genes implies the annotation of novel functions for underlying genes.