Genetic variation in the nuclear and organellar genomes modulates stochastic variation in the metabolome, growth, and defense

Highly expressed cell wall genes contribute to robustness of sepal size

Reproducibility in organ size and shape is a fascinating trait of living organisms. The mechanisms underlying such robustness remain, however, to be elucidated. Taking the sepal of Arabidopsis as a model, we investigated whether variability of gene expression plays a role in variation of organ size and shape. Previous work from our team identified cell-wall related genes as being enriched among the genes whose expression is highly variable. We then hypothesized that the variation of measured morphological parameters in cell-wall related single knockout mutants could be correlated with the variation in gene expression of the corresponding gene (the knocked-out gene) in wild-type plants. We analyzed sepal size and shape from 16 cell-wall mutants and found that sepal size variability correlates positively, not with gene expression variation, but with mean gene expression of the corresponding gene in wild type. These findings support a contribution of cell-wall related genes to the robustness of sepal size.

Embracing plant plasticity or robustness as a means of ensuring food security

The dual challenges of global population explosion and environmental deterioration represent major hurdles for 21st Century agriculture culminating in an unprecedented demand for food security. In this Review, we revisit historical concepts of plasticity and canalization before integrating them with contemporary studies of genotype-environment interactions (G×E) that are currently being carried out at the genome-wide level. In doing so we address both fundamental questions regarding G×E and potential strategies to best secure yields in both current and future climate scenarios.

Imprinting but not cytonuclear interactions determines seed size heterosis in Arabidopsis hybrids

The parent-of-origin effect on seeds can result from imprinting (unequal expression of paternal and maternal alleles) or combinational effects between cytoplasmic and nuclear genomes, but their relative contributions remain unknown. To discern these confounding factors, we produced cytoplasmic-nuclear substitution (CNS) lines using recurrent backcrossing in Arabidopsis (Arabidopsis thaliana) ecotypes Col-0 and C24. These CNS lines differed only in the nuclear genome (imprinting) or cytoplasm. The CNS reciprocal hybrids with the same cytoplasm displayed ∼20% seed size difference, whereas the seed size was similar between the reciprocal hybrids with fixed imprinting. Transcriptome analyses in the endosperm of CNS hybrids using laser-capture microdissection identified 104 maternally expressed genes (MEGs) and 90 paternally expressed genes (PEGs). These imprinted genes were involved in pectin catabolism and cell wall modification in the endosperm. Homeodomain Glabrous9 (HDG9), an epiallele and one of 11 cross-specific imprinted genes, affected seed size. In the embryo, there were a handful of imprinted genes in the CNS hybrids but only 1 was expressed at higher levels than in the endosperm. AT4G13495 was found to encode a long-noncoding RNA (lncRNA), but no obvious seed phenotype was observed in lncRNA knockout lines. Nuclear RNA Polymerase D1 (NRPD1), encoding the largest subunit of RNA Pol IV, was involved in the biogenesis of small interfering RNAs. Seed size and embryos were larger in the cross using nrpd1 as the maternal parent than in the reciprocal cross, supporting a role of the maternal NRPD1 allele in seed development. Although limited ecotypes were tested, these results suggest that imprinting and the maternal NRPD1-mediated small RNA pathway play roles in seed size heterosis in plant hybrids.

LTP2 hypomorphs show genotype-by-environment interaction in early seedling traits in Arabidopsis thaliana

Isogenic individuals can display seemingly stochastic phenotypic differences, limiting the accuracy of genotype-to-phenotype predictions. The extent of this phenotypic variation depends in part on genetic background, raising questions about the genes involved in controlling stochastic phenotypic variation. Focusing on early seedling traits in Arabidopsis thaliana, we found that hypomorphs of the cuticle-related gene LIPID TRANSFER PROTEIN 2 (LTP2) greatly increased variation in seedling phenotypes, including hypocotyl length, gravitropism and cuticle permeability. Many ltp2 hypocotyls were significantly shorter than wild-type hypocotyls while others resembled the wild-type. Differences in epidermal properties and gene expression between ltp2 seedlings with long and short hypocotyls suggest a loss of cuticle integrity as the primary determinant of the observed phenotypic variation. We identified environmental conditions that reveal or mask the increased variation in ltp2 hypomorphs and found that increased expression of its closest paralog LTP1 is necessary for ltp2 phenotypes. Our results illustrate how decreased expression of a single gene can generate starkly increased phenotypic variation in isogenic individuals in response to an environmental challenge.

LTP2 hypomorphs show genotype-by-environment interaction in early seedling traits in Arabidopsis thaliana

Isogenic individuals can display seemingly stochastic phenotypic differences, limiting the accuracy of genotype-to-phenotype predictions. The extent of this phenotypic variation depends in part on genetic background, raising questions about the genes involved in controlling stochastic phenotypic variation. Focusing on early seedling traits in Arabidopsis thaliana, we found that hypomorphs of the cuticle-related gene LTP2 greatly increased variation in seedling phenotypes, including hypocotyl length, gravitropism and cuticle permeability. Many ltp2 hypocotyls were significantly shorter than wild-type hypocotyls while others resembled the wild type. Differences in epidermal properties and gene expression between ltp2 seedlings with long and short hypocotyls suggest a loss of cuticle integrity as the primary determinant of the observed phenotypic variation. We identified environmental conditions that reveal or mask the increased variation in ltp2 hypomorphs, and found that increased expression of its closest paralog LTP1 is necessary for ltp2 phenotypes. Our results illustrate how decreased expression of a single gene can generate starkly increased phenotypic variation in isogenic individuals in response to an environmental challenge.

Imprinting but not cytonuclear interactions affects parent-of-origin effect on seed size in Arabidopsis hybrids

The parent-of-origin effect on seed size can result from imprinting or a combinational effect between cytoplasmic and nuclear genomes, but their relative contributions remain unknown. To discern these confounding effects, we generated cytoplasmic-nuclear substitution (CNS) lines using recurrent backcrossing in the Arabidopsis thaliana ecotypes Col-0 and C24. These CNS lines differ only in the nuclear genome (imprinting) or in the cytoplasm. The CNS reciprocal hybrids with the same cytoplasm display a ~20% seed size difference as observed in the conventional hybrids. However, seed size is similar between the reciprocal cybrids with fixed imprinting. Transcriptome analyses in the endosperm of CNS hybrids using laser-capture microdissection have identified 104 maternally expressed genes (MEGs) and 90 paternally-expressed genes (PEGs). These imprinted genes are involved in pectin catabolism and cell wall modification in the endosperm. HDG9, an epiallele and one of 11 cross-specific imprinted genes, controls seed size. In the embryo, a handful of imprinted genes is found in the CNS hybrids but only one is expressed higher in the embryo than endosperm. AT4G13495 encodes a long-noncoding RNA (lncRNA), but no obvious seed phenotype is observed in the lncRNA knockout lines. NRPD1, encoding the largest subunit of RNA Pol IV, is involved in the biogenesis of small interfering RNAs. Seed size and embryo is larger in the cross using nrpd1 as the maternal parent than in the reciprocal cross. In spite of limited ecotypes tested, these results suggest potential roles of imprinting and NRPD1-mediated small RNA pathway in seed size variation in hybrids.

Plant biochemical genetics in the multiomics era

Our understanding of plant biology has been revolutionized by modern genetics and biochemistry. However, biochemical genetics can be traced back to the foundation of Mendelian genetics; indeed, one of Mendel's milestone discoveries of seven characteristics of pea plants later came to be ascribed to a mutation in a starch branching enzyme. Here, we review both current and historical strategies for the elucidation of plant metabolic pathways and the genes that encode their component enzymes and regulators. We use this historical review to discuss a range of classical genetic phenomena including epistasis, canalization, and heterosis as viewed through the lens of contemporary high-throughput data obtained via the array of approaches currently adopted in multiomics studies.

PANTOTHENATE KINASE4, LOSS OF GDU2, and TRANSPOSON PROTEIN1 affect the canalization of tomato fruit metabolism

Most studies investigating quantitative traits focus on mean levels per genotype rather than the variation between different individuals of one genotype or the variation elicited by different environments. Consequently, the genes that govern this effect are not well understood. The concept, named canalization, which describes a lack of variation, is well known in the context of developmental processes but is poorly studied for quantitative traits such as metabolism. In this study, we selected 8 putative candidate genes from previously identified canalized metabolic quantitative trait loci and created genome-edited tomato (Solanum lycopersicum) mutants of these genes for experimental validation. Most lines showed wild-type morphology, except for an ARF-like protein mutant showing aberrant phenotypes in the form of scarred fruit cuticles. In greenhouse trials with different irrigation conditions, whole-plant traits showed a general increase of their level toward the more optimal irrigation conditions, whereas most metabolic traits showed an increase toward the other end of the gradient. Mutants of a PANTOTHENATE KINASE 4, the AIRP ubiquitin gene LOSS OF GDU2, and TRANSPOSON PROTEIN 1 grown under these conditions showed an overall improved plant performance. Additional effects, on both target and other metabolites in tomato fruits, regarding the mean level at specific conditions and, ergo, the cross-environment coefficient of variation, were observed. However, variation between individuals remained unaffected. In conclusion, this study supports the idea of distinct sets of genes regulating different types of variation.

Genetic and epigenetic control of the plant metabolome

Plant metabolites are mainly produced through chemical reactions catalysed by enzymes encoded in the genome. Mutations in enzyme-encoding or transcription factor-encoding genes can alter the metabolome by changing the enzyme's catalytic activity or abundance, respectively. Insertion of transposable elements into non-coding regions has also been reported to affect transcription and ultimately metabolite content. In addition to genetic mutations, transgenerational epigenetic variations have also been found to affect metabolic content by controlling the transcription of metabolism-related genes. However, the majority of cases reported so far, in which epigenetic mechanisms are associated with metabolism, are non-transgenerational, and are triggered by developmental signals or environmental stress. Although, accumulating research has provided evidence of strong genetic control of the metabolome, epigenetic control has been largely untouched. Here, we provide a review of the genetic and epigenetic control of metabolism with a focus on epigenetics. We discuss both transgenerational and non-transgenerational epigenetic marks regulating metabolism as well as prospects of the field of metabolic control where intricate interactions between genetics and epigenetics are involved.

Methylome decoding of RdDM-mediated reprogramming effects in the Arabidopsis MSH1 system

BACKGROUND

Plants undergo programmed chromatin changes in response to environment, influencing heritable phenotypic plasticity. The RNA-directed DNA methylation (RdDM) pathway is an essential component of this reprogramming process. The relationship of epigenomic changes to gene networks on a genome-wide basis has been elusive, particularly for intragenic DNA methylation repatterning.

RESULTS

Epigenomic reprogramming is tractable to detailed study and cross-species modeling in the MSH1 system, where perturbation of the plant-specific gene MSH1 triggers at least four distinct nongenetic states to impact plant stress response and growth vigor. Within this system, we have defined RdDM target loci toward decoding phenotype-relevant methylome data. We analyze intragenic methylome repatterning associated with phenotype transitions, identifying state-specific cytosine methylation changes in pivotal growth-versus-stress, chromatin remodeling, and RNA spliceosome gene networks that encompass 871 genes. Over 77% of these genes, and 81% of their central network hubs, are functionally confirmed as RdDM targets based on analysis of mutant datasets and sRNA cluster associations. These dcl2/dcl3/dcl4-sensitive gene methylation sites, many present as singular cytosines, reside within identifiable sequence motifs. These data reflect intragenic methylation repatterning that is targeted and amenable to prediction.

CONCLUSIONS

A prevailing assumption that biologically relevant DNA methylation variation occurs predominantly in density-defined differentially methylated regions overlooks behavioral features of intragenic, single-site cytosine methylation variation. RdDM-dependent methylation changes within identifiable sequence motifs reveal gene hubs within networks discriminating stress response and growth vigor epigenetic phenotypes. This study uncovers components of a methylome "code" for de novo intragenic methylation repatterning during plant phenotype transitions.

Adaptive significance of flowering time variation across natural seasonal environments in Arabidopsis thaliana

The relevance of flowering time variation and plasticity to climate adaptation requires a comprehensive empirical assessment. We investigated natural selection and the genetic architecture of flowering time in Arabidopsis through field experiments in Europe across multiple sites and seasons. We estimated selection for flowering time, plasticity and canalization. Loci associated with flowering time, plasticity and canalization by genome-wide association studies were tested for a geographic signature of climate adaptation. Selection favored early flowering and increased canalization, except at the northernmost site, but was rarely detected for plasticity. Genome-wide association studies revealed significant associations with flowering traits and supported a substantial polygenic inheritance. Alleles associated with late flowering, including functional FRIGIDA variants, were more common in regions experiencing high annual temperature variation. Flowering time plasticity to fall vs spring and summer environments was associated with GIGANTEA SUPPRESSOR 5, which promotes early flowering under decreasing day length and temperature. The finding that late flowering genotypes and alleles are associated with climate is evidence for past adaptation. Real-time phenotypic selection analysis, however, reveals pervasive contemporary selection for rapid flowering in agricultural settings across most of the species range. The response to this selection may involve genetic shifts in environmental cuing compared to the ancestral state.

Plasticity of rosette size in response to nitrogen availability is controlled by an RCC1-family protein

Nitrogen (N) is fundamental to plant growth, development and yield. Genes underlying N utilization and assimilation are well-characterized, but mechanisms underpinning plasticity of different phenotypes in response to N remain elusive. Here, using Arabidopsis thaliana accessions, we dissected the genetic architecture of plasticity in early and late rosette diameter, flowering time and yield, in response to three levels of N in the soil. Furthermore, we found that the plasticity in levels of primary metabolites were related with the plasticities of the studied traits. Genome-wide association analysis identified three significant associations for phenotypic plasticity, one for early rosette diameter and two for flowering time. We confirmed that the gene At1g19880, hereafter named as PLASTICITY OF ROSETTE TO NITROGEN 1 (PROTON1), encoding for a regulator of chromatin condensation 1 (RCC1) family protein, conferred plasticity of rosette diameter in response to N. Treatment of PROTON1 T-DNA line with salt implied that the reduced plasticity of early rosette diameter was not a general growth response to stress. We further showed that plasticities of growth and flowering-related traits differed between environmental cues, indicating decoupled genetic programs regulating these traits. Our findings provide a prospective to identify genes that stabilize performance under fluctuating environments.

An ABA-GA bistable switch can account for natural variation in the variability of Arabidopsis seed germination time

Genetically identical plants growing in the same conditions can display heterogeneous phenotypes. Here we use Arabidopsis seed germination time as a model system to examine phenotypic variability and its underlying mechanisms. We show extensive variation in seed germination time variability between Arabidopsis accessions and use a multiparent recombinant inbred population to identify two genetic loci involved in this trait. Both loci include genes implicated in modulating abscisic acid (ABA) sensitivity. Mutually antagonistic regulation between ABA, which represses germination, and gibberellic acid (GA), which promotes germination, underlies the decision to germinate and can act as a bistable switch. A simple stochastic model of the ABA-GA network shows that modulating ABA sensitivity can generate the range of germination time distributions we observe experimentally. We validate the model by testing its predictions on the effects of exogenous hormone addition. Our work provides a foundation for understanding the mechanism and functional role of phenotypic variability in germination time.

Metabolomic Analysis of Natural Variation in Arabidopsis

Methodological advances in coupled-mass spectrometry (gas chromatography and liquid chromatography; GC-MS and LC-MS) have rendered the profiling of highly complex plant extracts relatively facile and allowed that their high-throughput use aids the investigation of a range of biological questions. Among these is the elucidation of the genetic factors underlying metabolite abundance. For this purpose genome-wide association studies (GWAS) are being widely adopted in Arabidopsis with the resultant quantitative trait loci being subjected to cross-validation by the use of recombinant inbred lines, introgression lines, and T-DNA insertional knockout lines.

Does Gene Expression Noise Play a Functional Role in Plants?

Gene expression in individual cells can be surprisingly noisy. In unicellular organisms this noise can be functional; for example, by allowing a subfraction of the population to prepare for environmental stress. The role of gene expression noise in multicellular organisms has, however, remained unclear. In this review, we discuss how new techniques are revealing an unexpected level of variability in gene expression between and within genetically identical plants. We describe recent progress as well as speculate on the function of transcriptional noise as a mechanism for generating functional phenotypic diversity in plants.

Outstanding questions in flower metabolism

The great diversity of flowers, their color, odor, taste, and shape, is mostly a result of the metabolic processes that occur in this reproductive organ when the flower and its tissues develop, grow, and finally die. Some of these metabolites serve to advertise flowers to animal pollinators, other confer protection towards abiotic stresses, and a large proportion of the molecules of the central metabolic pathways have bioenergetic and signaling functions that support growth and the transition to fruits and seeds. Although recent studies have advanced our general understanding of flower metabolism, several questions still await an answer. Here, we have compiled a list of open questions on flower metabolism encompassing molecular aspects, as well as topics of relevance for agriculture and the ecosystem. These questions include the study of flower metabolism through development, the biochemistry of nectar and its relevance to promoting plant-pollinator interaction, recycling of metabolic resources after flowers whiter and die, as well as the manipulation of flower metabolism by pathogens. We hope with this review to stimulate discussion on the topic of flower metabolism and set a reference point to return to in the future when assessing progress in the field.

Segregation of an MSH1 RNAi transgene produces heritable non-genetic memory in association with methylome reprogramming

MSH1 is a plant-specific protein. RNAi suppression of MSH1 results in phenotype variability for developmental and stress response pathways. Segregation of the RNAi transgene produces non-genetic msh1 ‘memory’ with multi-generational inheritance. First-generation memory versus non-memory comparison, and six-generation inheritance studies, identifies gene-associated, heritable methylation repatterning. Genome-wide methylome analysis integrated with RNAseq and network-based enrichment studies identifies altered circadian clock networks, and phytohormone and stress response pathways that intersect with circadian control. A total of 373 differentially methylated loci comprising these networks are sufficient to discriminate memory from nonmemory full sibs. Methylation inhibitor 5-azacytidine diminishes the differences between memory and wild type for growth, gene expression and methylation patterning. The msh1 reprogramming is dependent on functional HISTONE DEACETYLASE 6 and methyltransferase MET1, and transition to memory requires the RNA-directed DNA methylation pathway. This system of phenotypic plasticity may serve as a potent model for defining accelerated plant adaptation during environmental change.

Segregation of an MSH1 RNAi transgene produces non-genetic memory that displays transgenerational inheritance in Arabidopsis. Here, the authors compare memory and non-memory full-sib progenies to show the involvement of DNA methylation reprogramming, involving the RdDM pathway, in transition to a heritable memory state.

The Role of Transcriptional Regulation in Hybrid Vigor

The genetic basis of hybrid vigor in plants remains largely unsolved but strong evidence suggests that variation in transcriptional regulation can explain many aspects of this phenomenon. Natural variation in transcriptional regulation is highly abundant in virtually all species and thus a potential source of heterotic variability. Allele Specific Expression (ASE), which is tightly linked to parent of origin effects and modulated by complex interactions in cis and in trans, is generally considered to play a key role in explaining the differences between hybrids and parental lines. Here we discuss the recent developments in elucidating the role of transcriptional variation in a number of aspects of hybrid vigor, thereby bridging old paradigms and hypotheses with contemporary research in various species.

Life History Variation as a Model for Understanding Trade-Offs in Plant-Environment Interactions

All plants must allocate limited resources to survival, growth, and reproduction. In natural species, allocation strategies reflect trade-offs between survivorship risk and subsequent fitness benefits and are therefore central to a species' ecology. Artificial selection on allocation has generated high-yielding crops that often invest the bare minimum in defense or longevity. Ecological, genetic, and evolutionary analyses of plant life history - particularly with respect to longevity and resource allocation along an axis from annual to perennial species - provides a framework to evaluate trade-offs in plant-environment interactions in natural and managed systems. Recent efforts to develop new model plant systems for research and to increase agricultural resilience and efficiency by developing herbaceous perennial crops motivates our critical assessment of traditional assumptions regarding differences between annual and perennial plant species. Here, we review our present understanding of the genetic basis of physiological, developmental, and anatomical differences in wild and crop species and reach two broad conclusions. First, that perenniality and annuality should be considered syndromes comprised of many interacting traits, and that elucidating the genetic basis of these traits is required to assess models of evolution and to develop successful breeding strategies. Modern phenomic and biotechnology tools will facilitate these enquiries. Second, many classic assumptions about the difference between the two syndromes are supported by limited evidence. Throughout this Review, we highlight key knowledge gaps in the proximate and ultimate mechanisms driving life history variation, and suggest empirical approaches to parameterize trade-offs and to make progress in this critical area of direct relevance to ecology and plant performance in a changing world.

Organellar protein multi-functionality and phenotypic plasticity in plants

With the increasing impact of climate instability on agricultural and ecological systems has come a heightened sense of urgency to understand plant adaptation mechanisms in more detail. Plant species have a remarkable ability to disperse their progeny to a wide range of environments, demonstrating extraordinary resiliency mechanisms that incorporate epigenetics and transgenerational stability. Surprisingly, some of the underlying versatility of plants to adapt to abiotic and biotic stress emerges from the neofunctionalization of organelles and organellar proteins. We describe evidence of possible plastid specialization and multi-functional organellar protein features that serve to enhance plant phenotypic plasticity. These features appear to rely on, for example, spatio-temporal regulation of plastid composition, and unusual interorganellar protein targeting and retrograde signalling features that facilitate multi-functionalization. Although we report in detail on three such specializations, involving MSH1, WHIRLY1 and CUE1 proteins in Arabidopsis, there is ample reason to believe that these represent only a fraction of what is yet to be discovered as we begin to elaborate cross-species diversity. Recent observations suggest that plant proteins previously defined in one context may soon be rediscovered in new roles and that much more detailed investigation of proteins that show subcellular multi-targeting may be warranted.

This article is part of the theme issue ‘Linking the mitochondrial genotype to phenotype: a complex endeavour’.

Organellar carbon metabolism is coordinated with distinct developmental phases of secondary xylem

Subcellular compartmentation of plant biosynthetic pathways in the mitochondria and plastids requires coordinated regulation of nuclear encoded genes, and the role of these genes has been largely ignored by wood researchers. In this study, we constructed a targeted systems genetics coexpression network of xylogenesis in Eucalyptus using plastid and mitochondrial carbon metabolic genes and compared the resulting clusters to the aspen xylem developmental series. The constructed network clusters reveal the organization of transcriptional modules regulating subcellular metabolic functions in plastids and mitochondria. Overlapping genes between the plastid and mitochondrial networks implicate the common transcriptional regulation of carbon metabolism during xylem secondary growth. We show that the central processes of organellar carbon metabolism are distinctly coordinated across the developmental stages of wood formation and are specifically associated with primary growth and secondary cell wall deposition. We also demonstrate that, during xylogenesis, plastid-targeted carbon metabolism is partially regulated by the central clock for carbon allocation towards primary and secondary xylem growth, and we discuss these networks in the context of previously established associations with wood-related complex traits. This study provides a new resolution into the integration and transcriptional regulation of plastid- and mitochondrial-localized carbon metabolism during xylogenesis.

Guidelines for Sample Normalization to Minimize Batch Variation for Large-Scale Metabolic Profiling of Plant Natural Genetic Variance

Recent methodological advances in both liquid chromatography-mass spectrometry (LC-MS) and gas chromatography-mass spectrometry (GC-MS) have facilitated the profiling highly complex mixtures of primary and secondary metabolites in order to investigate a diverse range of biological questions. These techniques usually face a large number of potential sources of technical and biological variation. In this chapter we describe guidelines and normalization procedures to reduce the analytical variation, which are essential for the high-throughput evaluation of metabolic variance used in broad genetic populations which commonly entail the evaluation of hundreds or thousands of samples. This chapter specifically deals with handling of large-scale plant samples for metabolomics analysis of quantitative trait loci (mQTL) in order to reduce analytical error as well as batch-to-batch variation.

Genetic basis of plasticity in plants

The ability of an organism to change its phenotype in response to different environments, termed plasticity, is a particularly important characteristic to enable sessile plants to adapt to rapid changes in their surroundings. Plasticity is a quantitative trait that can provide a fitness advantage and mitigate negative effects due to environmental perturbations. Yet, its genetic basis is not fully understood. Alongside technological limitations, the main challenge in studying plasticity has been the selection of suitable approaches for quantification of phenotypic plasticity. Here, we propose a categorization of the existing quantitative measures of phenotypic plasticity into nominal and relative approaches. Moreover, we highlight the recent advances in the understanding of the genetic architecture underlying phenotypic plasticity in plants. We identify four pillars for future research to uncover the genetic basis of phenotypic plasticity, with emphasis on development of computational approaches and theories. These developments will allow us to perform specific experiments to validate the causal genes for plasticity and to discover their role in plant fitness and evolution.

Metabolomics 20 years on: what have we learned and what hurdles remain?

The term metabolome was coined in 1998, by analogy to genome, transcriptome and proteome. The first research papers using the terms metabolomics, metabonomics, metabolic profiling or metabolite profiling were published shortly thereafter. In this short review we reflect on the major achievements brought about by the use of these approaches, and document the knowledge and technology gaps that are currently constraining its further development. Finally, we detail why we think that the time is ripe to refocus our efforts on the understanding of metabolic function.

Comparison of the Relative Potential for Epigenetic and Genetic Variation To Contribute to Trait Stability

The theoretical ability of epigenetic variation to influence the heritable variation of complex traits is gaining traction in the study of adaptation. This theory posits that epigenetic marks can control adaptive phenotypes but the relative potential of epigenetic variation in comparison to genetic variation in these traits is not presently understood. To compare the potential of epigenetic and genetic variation in adaptive traits, we analyzed the influence of DNA methylation variation on the accumulation of chemical defense compounds glucosinolates from the order Brassicales. Several decades of work on glucosinolates has generated extensive knowledge about their synthesis, regulation, genetic variation and contribution to fitness establishing this pathway as a model pathway for complex adaptive traits. Using high-throughput phenotyping with a randomized block design of ddm1 derived Arabidopsis thaliana epigenetic Recombinant Inbred Lines, we measured the correlation between DNA methylation variation and mean glucosinolate variation and within line stochastic variation. Using this information, we identified epigenetic Quantitative Trait Loci that contained specific Differentially Methylated Regions associated with glucosinolate traits. This showed that variation in DNA methylation correlates both with levels and variance of glucosinolates and flowering time with trait-specific loci. By conducting a meta-analysis comparing the results to different genetically variable populations, we conclude that the influence of DNA methylation variation on these adaptive traits is much lower than the corresponding impact of standing genetic variation. As such, selective pressure on these traits should mainly affect standing genetic variation to lead to adaptation.

Unsung and understudied: plastids involved in secondary growth

Plastids represent the only subcellular compartment where aromatic amino acid precursors for lignin can be synthesized during secondary growth in vascular plants. Despite this, aside from a general shared understanding that plastid-localized metabolism occurs during secondary growth, virtually no research has been performed on understanding their biology. Of particular importance will be insight into their ontogeny, morphology and ultrastructure, and (given the complex cytonuclear communication required) their nuclear-encoded and organellar-encoded regulation. Updating and integrating this knowledge will contribute to our fundamental understanding of a ubiquitous developmental process in vascular plants, and a major terrestrial carbon sink, as well as carbon-related plant biotechnology. Given available evidence, we propose a new name for a distinct plastid derivative-the 'xyloplast', is required.

Canalization of Tomato Fruit Metabolism

To explore the genetic robustness (canalization) of metabolism, we examined the levels of fruit metabolites in multiple harvests of a tomato introgression line (IL) population. The IL partitions the whole genome of the wild species Solanum pennellii in the background of the cultivated tomato (Solanum lycopersicum). We identified several metabolite quantitative trait loci that reduce variability for both primary and secondary metabolites, which we named canalization metabolite quantitative trait loci (cmQTL). We validated nine cmQTL using an independent population of backcross inbred lines, derived from the same parents, which allows increased resolution in mapping the QTL previously identified in the ILs. These cmQTL showed little overlap with QTL for the metabolite levels themselves. Moreover, the intervals they mapped to harbored few metabolism-associated genes, suggesting that the canalization of metabolism is largely controlled by regulatory genes.

The Genetics of Plant Metabolism

Plant metabolic studies have traditionally focused on the role and regulation of the enzymes catalyzing key reactions within specific pathways. Within the past 20 years, reverse genetic approaches have allowed direct determination of the effects of the deficiency, or surplus, of a given protein on the biochemistry of a plant. In parallel, top-down approaches have also been taken, which rely on screening broad, natural genetic diversity for metabolic diversity. Here, we compare and contrast the various strategies that have been adopted to enhance our understanding of the natural diversity of metabolism. We also detail how these approaches have enhanced our understanding of both specific and global aspects of the genetic regulation of metabolism. Finally, we discuss how such approaches are providing important insights into the evolution of plant secondary metabolism.

An Integrative Genetic Study of Rice Metabolism, Growth and Stochastic Variation Reveals Potential C/N Partitioning Loci

Studying the genetic basis of variation in plant metabolism has been greatly facilitated by genomic and metabolic profiling advances. In this study, we use metabolomics and growth measurements to map QTL in rice, a major staple crop. Previous rice metabolism studies have largely focused on identifying genes controlling major effect loci. To complement these studies, we conducted a replicated metabolomics analysis on a japonica (Lemont) by indica (Teqing) rice recombinant inbred line population and focused on the genetic variation for primary metabolism. Using independent replicated studies, we show that in contrast to other rice studies, the heritability of primary metabolism is similar to Arabidopsis. The vast majority of metabolic QTLs had small to moderate effects with significant polygenic epistasis. Two metabolomics QTL hotspots had opposing effects on carbon and nitrogen rich metabolites suggesting that they may influence carbon and nitrogen partitioning, with one locus co-localizing with SUSIBA2 (WRKY78). Comparing QTLs for metabolomic and a variety of growth related traits identified few overlaps. Interestingly, the rice population displayed fewer loci controlling stochastic variation for metabolism than was found in Arabidopsis. Thus, it is possible that domestication has differentially impacted stochastic metabolite variation more than average metabolite variation.

Developmental mechanisms underlying variable, invariant and plastic phenotypes

BACKGROUND

Discussions of phenotypic robustness often consider scenarios where invariant phenotypes are optimal and assume that developmental mechanisms have evolved to buffer the phenotypes of specific traits against stochastic and environmental perturbations. However, plastic plant phenotypes that vary between environments or variable phenotypes that vary stochastically within an environment may also be advantageous in some scenarios.

SCOPE

Here the conditions under which invariant, plastic and variable phenotypes of specific traits may confer a selective advantage in plants are examined. Drawing on work from microbes and multicellular organisms, the mechanisms that may give rise to each type of phenotype are discussed.

CONCLUSION

In contrast to the view of robustness as being the ability of a genotype to produce a single, invariant phenotype, changes in a phenotype in response to the environment, or phenotypic variability within an environment, may also be delivered consistently (i.e. robustly). Thus, for some plant traits, mechanisms have probably evolved to produce plasticity or variability in a reliable manner.

Cytonuclear interactions affect adaptive traits of the annual plant Arabidopsis thaliana in the field

Although the contribution of cytonuclear interactions to plant fitness variation is relatively well documented at the interspecific level, the prevalence of cytonuclear interactions at the intraspecific level remains poorly investigated. In this study, we set up a field experiment to explore the range of effects that cytonuclear interactions have on fitness-related traits in Arabidopsis thaliana To do so, we created a unique series of 56 cytolines resulting from cytoplasmic substitutions among eight natural accessions reflecting within-species genetic diversity. An assessment of these cytolines and their parental lines scored for 28 adaptive whole-organism phenotypes showed that a large proportion of phenotypic traits (23 of 28) were affected by cytonuclear interactions. The effects of these interactions varied from slight but frequent across cytolines to strong in some specific parental pairs. Two parental pairs accounted for half of the significant pairwise interactions. In one parental pair, Ct-1/Sha, we observed symmetrical phenotypic responses between the two nuclear backgrounds when combined with specific cytoplasms, suggesting nuclear differentiation at loci involved in cytonuclear epistasis. In contrast, asymmetrical phenotypic responses were observed in another parental pair, Cvi-0/Sha. In the Cvi-0 nuclear background, fecundity and phenology-related traits were strongly affected by the Sha cytoplasm, leading to a modified reproductive strategy without penalizing total seed production. These results indicate that natural variation in cytoplasmic and nuclear genomes interact to shape integrative traits that contribute to adaptation, thereby suggesting that cytonuclear interactions can play a major role in the evolutionary dynamics ofA. thaliana.

Comprehensive Transcriptome Analysis of Six Catfish Species from an Altitude Gradient Reveals Adaptive Evolution in Tibetan Fishes

Glyptosternoid fishes (Siluriformes), one of the three broad fish lineages (the two other are schizothoracines and Triplophysa), have a limited distribution in the rivers in the Tibetan Plateau and peripheral regions. To investigate the genetic mechanisms underlying adaptation to the Tibetan Plateau in several fish species from gradient altitudes, a total of 20,659,183–37,166,756 sequence reads from six species of catfish were generated by Illumina sequencing, resulting in six assemblies. Analysis of the 1,656 orthologs among the six assembled catfish unigene sets provided consistent evidence for genome-wide accelerated evolution in the three glyptosternoid lineages living at high altitudes. A large number of genes refer to functional categories related to hypoxia and energy metabolism exhibited rapid evolution in the glyptosternoid lineages relative to yellowhead catfish living in plains areas. Genes showing signatures of rapid evolution and positive selection in the glyptosternoid lineages were also enriched in functions associated with energy metabolism and hypoxia. Our analyses provide novel insights into highland adaptation in fishes and can serve as a foundation for future studies aiming to identify candidate genes underlying the genetic basis of adaptation in Tibetan fishes.

Natural Variation of Plant Metabolism: Genetic Mechanisms, Interpretive Caveats, and Evolutionary and Mechanistic Insights

Combining quantitative genetics studies with metabolomics/metabolic profiling platforms, genomics, and transcriptomics is creating significant progress in identifying the causal genes controlling natural variation in metabolite accumulations and profiles. In this review, we discuss key mechanistic and evolutionary insights that are arising from these studies. This includes the potential role of transport and other processes in leading to a separation of the site of mechanistic causation and metabolic consequence. A reilluminated observation is the potential for genomic variation in the organelle to alter phenotypic variation alone and in epistatic interaction with the nuclear genetic variation. These studies are also highlighting new aspects of metabolic pleiotropy both in terms of the breadth of loci altering metabolic variation as well as the potential for broader effects on plant defense regulation of the metabolic variation than has previously been predicted. We also illustrate caveats that can be overlooked when translating quantitative genetics descriptors such as heritability and per-locus r(2) to mechanistic or evolutionary interpretations.

Integrative Approaches to Enhance Understanding of Plant Metabolic Pathway Structure and Regulation

Huge insight into molecular mechanisms and biological network coordination have been achieved following the application of various profiling technologies. Our knowledge of how the different molecular entities of the cell interact with one another suggests that, nevertheless, integration of data from different techniques could drive a more comprehensive understanding of the data emanating from different techniques. Here, we provide an overview of how such data integration is being used to aid the understanding of metabolic pathway structure and regulation. We choose to focus on the pairwise integration of large-scale metabolite data with that of the transcriptomic, proteomics, whole-genome sequence, growth- and yield-associated phenotypes, and archival functional genomic data sets. In doing so, we attempt to provide an update on approaches that integrate data obtained at different levels to reach a better understanding of either single gene function or metabolic pathway structure and regulation within the context of a broader biological process.

Whole genome resequencing of Botrytis cinerea isolates identifies high levels of standing diversity

How standing genetic variation within a pathogen contributes to diversity in host/pathogen interactions is poorly understood, partly because most studied pathogens are host-specific, clonally reproducing organisms which complicates genetic analysis. In contrast, Botrytis cinerea is a sexually reproducing, true haploid ascomycete that can infect a wide range of diverse plant hosts. While previous work had shown significant genomic variation between two isolates, we proceeded to assess the level and frequency of standing variation in a population of B. cinerea. To begin measuring standing genetic variation in B. cinerea, we re-sequenced the genomes of 13 different isolates and aligned them to the previously sequenced T4 reference genome. In addition one of these isolates was resequenced from four independently repeated cultures. A high level of genetic diversity was found within the 13 isolates. Within this variation, we could identify clusters of genes with major effect polymorphisms, i.e., polymorphisms that lead to a predicted functional knockout, that surrounded genes involved in controlling vegetative incompatibility. The genotype at these loci was able to partially predict the interaction of these isolates in vegetative fusion assays showing that these loci control vegetative incompatibility. This suggests that the vegetative incompatibility loci within B. cinerea are associated with regions of increased genetic diversity. The genome re-sequencing of four clones from the one isolate (Grape) that had been independently propagated over 10 years showed no detectable spontaneous mutation. This suggests that B. cinerea does not display an elevated spontaneous mutation rate. Future work will allow us to test if, and how, this diversity may be contributing to the pathogen's broad host range.