Genomic approaches to analyzing natural variation in Arabidopsis thaliana

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.

Epistasis × environment interactions among Arabidopsis thaliana glucosinolate genes impact complex traits and fitness in the field

Despite the growing number of studies showing that genotype × environment and epistatic interactions control fitness, the influences of epistasis × environment interactions on adaptive trait evolution remain largely uncharacterized. Across three field trials, we quantified aliphatic glucosinolate (GSL) defense chemistry, leaf damage, and relative fitness using mutant lines of Arabidopsis thaliana varying at pairs of causal aliphatic GSL defense genes to test the impact of epistatic and epistasis × environment interactions on adaptive trait variation. We found that aliphatic GSL accumulation was primarily influenced by additive and epistatic genetic variation, leaf damage was primarily influenced by environmental variation and relative fitness was primarily influenced by epistasis and epistasis × environment interactions. Epistasis × environment interactions accounted for up to 48% of the relative fitness variation in the field. At a single field site, the impact of epistasis on relative fitness varied significantly over 2 yr, showing that epistasis × environment interactions within a location can be temporally dynamic. These results suggest that the environmental dependency of epistasis can profoundly influence the response to selection, shaping the adaptive trajectories of natural populations in complex ways, and deserves further consideration in future evolutionary studies.

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

Recent studies are starting to show that genetic control over stochastic variation is a key evolutionary solution of single celled organisms in the face of unpredictable environments. This has been expanded to show that genetic variation can alter stochastic variation in transcriptional processes within multi-cellular eukaryotes. However, little is known about how genetic diversity can control stochastic variation within more non-cell autonomous phenotypes. Using an Arabidopsis reciprocal RIL population, we showed that there is significant genetic diversity influencing stochastic variation in the plant metabolome, defense chemistry, and growth. This genetic diversity included loci specific for the stochastic variation of each phenotypic class that did not affect the other phenotypic classes or the average phenotype. This suggests that the organism's networks are established so that noise can exist in one phenotypic level like metabolism and not permeate up or down to different phenotypic levels. Further, the genomic variation within the plastid and mitochondria also had significant effects on the stochastic variation of all phenotypic classes. The genetic influence over stochastic variation within the metabolome was highly metabolite specific, with neighboring metabolites in the same metabolic pathway frequently showing different levels of noise. As expected from bet-hedging theory, there was more genetic diversity and a wider range of stochastic variation for defense chemistry than found for primary metabolism. Thus, it is possible to begin dissecting the stochastic variation of whole organismal phenotypes in multi-cellular organisms. Further, there are loci that modulate stochastic variation at different phenotypic levels. Finding the identity of these genes will be key to developing complete models linking genotype to phenotype.

Arabidopsis thaliana as a model organism in systems biology

Significant progress has been made in identification of genes and gene networks involved in key biological processes. Yet, how these genes and networks are coordinated over increasing levels of biological complexity, from cells to tissues to organs, remains unclear. To address complex biological questions, biologists are increasingly using high-throughput tools and systems biology approaches to examine complex biological systems at a global scale. A system is a network of interacting and interdependent components that shape the system's unique properties. Systems biology studies the organization of system components and their interactions, with the idea that unique properties of that system can be observed only through study of the system as a whole. The application of systems biology approaches to questions in plant biology has been informative. In this review, we give examples of how systems biology is currently being used in Arabidopsis to investigate the transcriptional networks regulating root development, the metabolic response to stress, and the genetic regulation of metabolic variability. From these studies, we are beginning obtain sufficient data to generate more accurate models for system function. Further investigation of plant systems will require data gathering from specific cells and tissues, continued improvement in metabolic technologies, and novel computational methods for data visualization and modeling.

The effect of growth conditions on the seed size/number trade-off

BACKGROUND

If the amount of resources allocated to reproduction (K) is fixed, then an increase in seed mass (S) can only be achieved by a decrease in seed number (n = K/S). Thus, log(n) = log(K)-log(S) producing a slope of -1 when seed mass and number are plotted on log-log axes. However, in comparative studies, empirical support for a slope of -1 is limited and contentious, leading some to question the utility of this concept.

METHODOLOGY/PRINCIPAL FINDINGS

First, we show that the expected slope depends on whether genotypes and species producing seeds of different mass are expected to reach the same adult size and that this in turn depends partly on the nature of growth. Second, we present experimental results using a population of recombinant inbred lines (RILs) of Arabidopsis thaliana. When these RILs are grown in large pots with plentiful nutrients, they exhibit a trade-off between seed size and number with a slope of -1.68 (+/-0.18) on log-log axes. This occurs because of genetic correlations between seed mass and adult size so that, under the right growth conditions, lines producing lighter seeds have the genetic potential to produce larger rosettes and hence a greater total mass of seeds. We re-grew lines in small pots (10 and 40 mm diameter) in a nutrient-poor substrate so that final adult size was heavily restricted by pot size.

CONCLUSIONS/SIGNIFICANCE

Under our growth conditions, small-seeded lines were unable to produce a greater total mass of seeds. Hence a trade-off emerged between seed mass and seed number with a slope of -1.166+/-0.319 on log-log axes in 40-mm diameter pots (close to the expected value of -1), although the slope was 0.132+/-0.263 in 10-mm diameter pots, demonstrating that the nature of the trade-off is sensitive to the growth conditions.

Phytochrome B and histone deacetylase 6 control light-induced chromatin compaction in Arabidopsis thaliana

Natural genetic variation in Arabidopsis thaliana exists for many traits and often reflects acclimation to local environments. Studying natural variation has proven valuable in the characterization of phenotypic traits and, in particular, in identifying genetic factors controlling these traits. It has been previously shown that chromatin compaction changes during development and biotic stress. To gain more insight into the genetic control of chromatin compaction, we investigated the nuclear phenotype of 21 selected Arabidopsis accessions from different geographic origins and habitats. We show natural variation in chromatin compaction and demonstrate a positive correlation with latitude of geographic origin. The level of compaction appeared to be dependent on light intensity. A novel approach, combining Quantitative Trait Locus (QTL) mapping and microscopic examination, pointed at PHYTOCHROME-B (PHYB) and HISTONE DEACETYLASE-6 (HDA6) as positive regulators of light-controlled chromatin compaction. Indeed, mutant analyses demonstrate that both factors affect global chromatin organization. HDA6, in addition, strongly promotes the light-mediated compaction of the Nucleolar Organizing Regions (NORs). The accession Cape Verde Islands-0 (Cvi-0), which shows sequence polymorphism in the PHYB gene and in the HDA6 promotor, resembles the hda6 mutant in having reduced chromatin compaction and decreased methylation levels of DNA and histone H3K9 at the NORs. We provide evidence that chromatin organization is controlled by light intensity. We propose that chromatin plasticity is associated with acclimation of Arabidopsis to its environment. The polymorphic alleles such as PHYB and HDA6 control this process.

The habitat of the plant model species Arabidopsis thaliana can be found throughout the Northern hemisphere. As a consequence, individual populations have acclimated to a great diversity of environmental conditions. This is reflected by a wealth of natural genetic variation in many phenotypic traits. We utilized this natural variation via a novel approach, combining microscopic examination, quantitative genetics, and analysis of environmental parameters, to understand the regulation of nuclear chromatin compaction in leaf mesophyll cells. We show that the level of chromatin compaction among natural Arabidopsis thaliana accessions correlates with latitude of origin and depends on local light intensity. Our study provides evidence that the photoreceptor PHYTOCHROME-B (PHYB) and the histone modifier HISTONE DEACETYLASE 6 (HDA6) are positive regulators of global chromatin organization in a light-dependent manner. In addition, HDA6 specifically controls light-mediated chromatin compaction of the Nucleolar Organizing Regions (NORs). We propose that the observed light-controlled plasticity of chromatin plays a role in acclimation and survival of plants in their natural environment.

Quantitative genetic analysis of thermal dissipation in Arabidopsis

Feedback deexcitation is a photosynthetic regulatory mechanism that can protect plants from high light stress by harmlessly dissipating excess absorbed light energy as heat. To understand the genetic basis for intraspecies differences in thermal dissipation capacity, we investigated natural variation in Arabidopsis (Arabidopsis thaliana). We determined the variation in the amount of thermal dissipation by measuring nonphotochemical quenching (NPQ) of chlorophyll fluorescence in Arabidopsis accessions of diverse origins. Ll-1 and Sf-2 were selected as high NPQ Arabidopsis accessions, and Columbia-0 (Col-0) and Wassilewskija-2 were selected as relatively low NPQ accessions. In spite of significant differences in NPQ, previously identified NPQ factors were indistinguishable between the high and the low NPQ accessions. Intermediate levels of NPQ in Ll-1 x Col-0 F1 and Sf-2 x Col-0 F1 compared to NPQ levels in their parental lines and continuous distribution of NPQ in F2 indicated that the variation in NPQ is under the control of multiple nuclear factors. To identify genetic factors responsible for the NPQ variation, we developed a polymorphic molecular marker set for Sf-2 x Col-0 at approximately 10-centimorgan intervals. From quantitative trait locus (QTL) mapping with undistorted genotype data and NPQ measurements in an F2 mapping population, we identified two high NPQ QTLs, HQE1 (high qE 1, for high energy-dependent quenching 1) and HQE2, on chromosomes 1 and 2, and the phenotype of HQE2 was validated by analysis of near isogenic lines. Neither QTL maps to a gene that had been identified previously in extensive forward genetics screens using induced mutants, suggesting that quantitative genetics can be used to find new genes affecting thermal dissipation.

Quantitative genomics: analyzing intraspecific variation using global gene expression polymorphisms or eQTLs

Scientific inquiries in fields ranging from ecology to plant breeding assess phenotypic variation within a plant species either to explain its presence or utilize its consequences. Frequently this natural genetic variation is studied via mapping quantitative trait loci (QTLs); however, elucidation of the underlying molecular mechanisms is a continuing bottleneck. The genomic analysis of transcripts as individual phenotypes has led to the emerging field of expression QTL analysis. This field has begun both to delve into the ecological/evolutionary significance of this transcript variation as well as to use specific eQTLs to speed up our analysis of the molecular basis of quantitative traits. This review introduces eQTL analysis and begins to illustrate how these data can be applied to multiple research fields.

Brachypodium distachyon: making hay with a wild grass

Brachypodium distachyon is a wild grass with a short life cycle. Although it is related to small grain cereals such as wheat, its genome is only a fraction of the size. A draft genome sequence is currently available, and molecular and genetic tools have been developed for transformation, mutagenesis and gene mapping. Accessions collected from across its ancestral range show a surprising degree of phenotypic variation in many traits, including those implicated in domestication of the cereals. Thus, given its rapid cycling time and ease of cultivation, Brachypodium will be a useful model for investigating problems in grass biology.

Identifying the molecular basis of QTLs: eQTLs add a new dimension

Natural genetic variation within plant species is at the core of plant science ranging from agriculture to evolution. Whereas much progress has been made in mapping quantitative trait loci (QTLs) controlling this natural variation, the elucidation of the underlying molecular mechanisms has remained a bottleneck. Recent systems biology tools have significantly shortened the time required to proceed from a mapped locus to testing of candidate genes. These tools enable research on natural variation to move from simple reductionistic studies focused on individual genes to integrative studies connecting molecular variation at multiple loci with physiological consequences. This review focuses on recent examples that demonstrate how expression QTL data can be used for gene discovery and exploited to untangle complex regulatory networks.

Genetic analysis of the physiological responses to low boron stress in Arabidopsis thaliana

Boron (B) is an essential micronutrient for higher plants. There is wide genetic variation in the response to B deficiency among plant species and cultivars. The objective of this study was to identify quantitative trait loci (QTL) that control B efficiency in natural Arabidopsis accessions. The B efficiency coefficient (BEC) and seed yield under low B conditions (SYLB) were investigated by solution culture in two separate experiments in an Arabidopsis recombinant inbred line (RIL) population. Both of the traits studied exhibited high transgressive variation in the RIL population, and, in total, five and three QTL were identified for BEC and SYLB, respectively. Three of the five QTL, including the QTL, AtBE1-2, that has a large effect on the BEC, were found at the interval of the corresponding QTL for SYLB in both experiments. The close genetic relationship between BEC and SYLB was further confirmed by conditional QTL mapping in the RIL population and unconditional QTL mapping in an AtBE1-2-segregated F(2) population. Epistatic interactions for the tested traits were analysed, and were found to be widespread in the detected QTL of Arabidopsis in the RIL population. Comparison of the QTL interval for B efficiency with reported B-related genes showed that 10 B-related genes, together with one BOR1 homolog (BOR5, At1g74810) were located in the QTL region of AtBE1-2. These results suggest that natural variation in B efficiency in Arabidopsis has a complex molecular basis. They also provide a basis for fine mapping and cloning of the B-efficiency genes, with the ultimate aim of discovering the physiological mechanism of action of the genes.

Whole-plant mineral partitioning throughout the life cycle in Arabidopsis thaliana ecotypes Columbia, Landsberg erecta, Cape Verde Islands, and the mutant line ysl1ysl3

Minimal information exists on whole-plant dynamics of mineral flow through Arabidopsis thaliana or on the source tissues responsible for mineral export to developing seeds. Understanding these phenomena in a model plant could help in the development of nutritionally enhanced crop cultivars. A whole-plant partitioning study, using sequential harvests, was conducted to characterize growth and mineral concentrations and contents of rosettes, cauline leaves, stems, immature fruit, mature fruit hulls, and seeds of three WT lines (Col-0, Ler, and Cvi) and one mutant line (Col-0::ysl1ysl3). Shoot mineral content increased throughout the life cycle for all minerals, although tissue-specific mineral partitioning differed between genotypes. In particular, iron (Fe), zinc (Zn), and copper (Cu) were aberrantly distributed in ysl1ysl3. Remobilization was observed for several minerals from various tissues, including cauline leaves and silique hulls, but the amounts were generally far below the total mineral accretion observed in seeds. When YSL1 and YSL3 are nonfunctional, Cu, Fe, and Zn are not effectively remobilized from, or do not effectively pass through, leaf and maternal fruit tissues. With respect to seed mineral accretion in Arabidopsis, continued uptake and translocation of minerals to source tissues during seed fill are as important, if not more important, than remobilization of previously stored minerals.

Development of a near-isogenic line population of Arabidopsis thaliana and comparison of mapping power with a recombinant inbred line population

In Arabidopsis recombinant inbred line (RIL) populations are widely used for quantitative trait locus (QTL) analyses. However, mapping analyses with this type of population can be limited because of the masking effects of major QTL and epistatic interactions of multiple QTL. An alternative type of immortal experimental population commonly used in plant species are sets of introgression lines. Here we introduce the development of a genomewide coverage near-isogenic line (NIL) population of Arabidopsis thaliana, by introgressing genomic regions from the Cape Verde Islands (Cvi) accession into the Landsberg erecta (Ler) genetic background. We have empirically compared the QTL mapping power of this new population with an already existing RIL population derived from the same parents. For that, we analyzed and mapped QTL affecting six developmental traits with different heritability. Overall, in the NIL population smaller-effect QTL than in the RIL population could be detected although the localization resolution was lower. Furthermore, we estimated the effect of population size and of the number of replicates on the detection power of QTL affecting the developmental traits. In general, population size is more important than the number of replicates to increase the mapping power of RILs, whereas for NILs several replicates are absolutely required. These analyses are expected to facilitate experimental design for QTL mapping using these two common types of segregating populations.

Natural variation for carbohydrate content in Arabidopsis. Interaction with complex traits dissected by quantitative genetics

Besides being a metabolic fuel, carbohydrates play important roles in plant growth and development, in stress responses, and as signal molecules. We exploited natural variation in Arabidopsis (Arabidopsis thaliana) to decipher the genetic architecture determining carbohydrate content. A quantitative trait locus (QTL) approach in the Bay-0 x Shahdara progeny grown in two contrasting nitrogen environments led to the identification of 39 QTLs for starch, glucose, fructose, and sucrose contents representing at least 14 distinct polymorphic loci. A major QTL for fructose content (FR3.4) and a QTL for starch content (ST3.4) were confirmed in heterogeneous inbred families. Several genes associated with carbon (C) metabolism colocalize with the identified QTL. QTLs for senescence-related traits, and for flowering time, water status, and nitrogen-related traits, previously detected with the same genetic material, colocalize with C-related QTLs. These colocalizations reflect the complex interactions of C metabolism with other physiological processes. QTL fine-mapping and cloning could thus lead soon to the identification of genes potentially involved in the control of different connected physiological processes.

Mapping quantitative trait loci associated with selenate tolerance in Arabidopsis thaliana

Selenium is essential for many organisms, but is toxic at higher levels. To investigate the genetic basis of selenate tolerance in Arabidopsis thaliana, quantitative trait loci (QTL) associated with selenate tolerance in accessions Landsberg erecta and Columbia were mapped using recombinant inbred lines (RILs). The selenate tolerance index (TI(D10) = root growth + 30 microm selenate/root growth control x 100%) was fourfold higher for parental line Col-4 (59%) than for parent Ler-0 (15%). Among the 96 F8 RILs, TI(D10) ranged from 11 to 75% (mean 37%). Using composite interval mapping, three QTL were found on chromosomes 1, 3 and 5, which together explained 24% of variation in TI(D10) and 32% of the phenotypic variation for the difference in root length +/- Se (RL(D10)). Highly significant epistatic interactions between the QTL and markers on chromosome 2 explained additional variation for both traits. Potential candidate genes for Se tolerance in each of the QTL regions are discussed. These results offer insight into the genetic basis of selenate tolerance, and may be useful for identification of selenate-tolerance genes.

Role of FRIGIDA and FLOWERING LOCUS C in determining variation in flowering time of Arabidopsis

Arabidopsis (Arabidopsis thaliana) accessions provide an excellent resource to dissect the molecular basis of adaptation. We have selected 192 Arabidopsis accessions collected to represent worldwide and local variation and analyzed two adaptively important traits, flowering time and vernalization response. There was huge variation in the flowering habit of the different accessions, with no simple relationship to latitude of collection site and considerable diversity occurring within local regions. We explored the contribution to this variation from the two genes FRIGIDA (FRI) and FLOWERING LOCUS C (FLC), previously shown to be important determinants in natural variation of flowering time. A correlation of FLC expression with flowering time and vernalization was observed, but it was not as strong as anticipated due to many late-flowering/vernalization-requiring accessions being associated with low FLC expression and early-flowering accessions with high FLC expression. Sequence analysis of FRI revealed which accessions were likely to carry functional alleles, and, from comparison of flowering time with allelic type, we estimate that approximately 70% of flowering time variation can be accounted for by allelic variation of FRI. The maintenance and propagation of 20 independent nonfunctional FRI haplotypes suggest that the loss-of-function mutations can confer a strong selective advantage. Accessions with a common FRI haplotype were, in some cases, associated with very different FLC levels and wide variation in flowering time, suggesting additional variation at FLC itself or other genes regulating FLC. These data reveal how useful these Arabidopsis accessions will be in dissecting the complex molecular variation that has led to the adaptive phenotypic variation in flowering time.

Natural variation in expression of self-incompatibility in Arabidopsis thaliana: implications for the evolution of selfing

The switch from an out-crossing to a self-fertilizing mating system is one of the most prevalent evolutionary trends in plant reproduction and is thought to have occurred repeatedly in flowering plants. However, little is known about the evolution of self-fertility and the genetic architecture of selfing. Here, we establish Arabidopsis thaliana as a model for genetic analysis of the switch to self-fertility in the crucifer family, where the ancestral out-crossing mode of mating is determined by self-incompatibility (SI), a genetic system controlled by the S locus. We show that A. thaliana ecotypes exhibit S-locus polymorphisms and differ in their ability to express the SI trait upon transformation with S-locus genes derived from the obligate out-crosser Arabidopsis lyrata. Remarkably, at least one ecotype was reverted to a stable, self-incompatible phenotype identical to that of naturally self-incompatible species. These ecotype differences are heritable and reflect the fixation in different A. thaliana populations of independent mutations that caused or enforced the switch to self-fertility. Their continued analysis promises to identify the loci that were the targets of natural selection for selfing and to contribute to a mechanistic understanding of the SI response.

Plant development: slowing root growth naturally

What determines organ size? A screen for plant strains showing natural variation in root length has identified a novel gene that modulates root growth rate by controlling cell division and elongation at the root tip.