Root System Depth in Arabidopsis Is Shaped by EXOCYST70A3 via the Dynamic Modulation of Auxin Transport.

Root system architecture (RSA), the distribution of roots in soil, plays a major role in plant survival. RSA is shaped by multiple developmental processes that are largely governed by the phytohormone auxin, suggesting that auxin regulates responses of roots that are important for local adaptation. However, auxin has a central role in numerous processes, and it is unclear which molecular mechanisms contribute to the variation in RSA for environmental adaptation. Using natural variation in Arabidopsis, we identify EXOCYST70A3 as a modulator of the auxin system that causes variation in RSA by acting on PIN4 protein distribution. Allelic variation and genetic perturbation of EXOCYST70A3 lead to alteration of root gravitropic responses, resulting in a different RSA depth profile and drought resistance. Overall our findings suggest that the local modulation of the pleiotropic auxin pathway can gives rise to distinct RSAs that can be adaptive in specific environments.

Plant genetic effects on microbial hubs impact host fitness in repeated field trials.

Although complex interactions between hosts and microbial associates are increasingly well documented, we still know little about how and why hosts shape microbial communities in nature. In addition, host genetic effects on microbial communities vary widely depending on the environment, obscuring conclusions about which microbes are impacted and which plant functions are important. We characterized the leaf microbiota of 200 Arabidopsis thaliana genotypes in eight field experiments and detected consistent host effects on specific, broadly distributed microbial species (operational taxonomic unit [OTUs]). Host genetic effects disproportionately influenced central ecological hubs, with heritability of particular OTUs declining with their distance from the nearest hub within the microbial network. These host effects could reflect either OTUs preferentially associating with specific genotypes or differential microbial success within them. Host genetics associated with microbial hubs explained over 10% of the variation in lifetime seed production among host genotypes across sites and years. We successfully cultured one of these microbial hubs and demonstrated its growth-promoting effects on plants in sterile conditions. Finally, genome-wide association mapping identified many putatively causal genes with small effects on the relative abundance of microbial hubs across sites and years, and these genes were enriched for those involved in the synthesis of specialized metabolites, auxins, and the immune system. Using untargeted metabolomics, we corroborate the consistent association between variation in specialized metabolites and microbial hubs across field sites. Together, our results reveal that host genetic variation impacts the microbial communities in consistent ways across environments and that these effects contribute to fitness variation among host genotypes.

Multiple FLC haplotypes defined by independent cis-regulatory variation underpin life history diversity in Arabidopsis thaliana.

Relating molecular variation to phenotypic diversity is a central goal in evolutionary biology. In Arabidopsis thaliana, FLOWERING LOCUS C (FLC) is a major determinant of variation in vernalization--the acceleration of flowering by prolonged cold. Here, through analysis of 1307 A. thaliana accessions, we identify five predominant FLC haplotypes defined by noncoding sequence variation. Genetic and transgenic experiments show that they are functionally distinct, varying in FLC expression level and rate of epigenetic silencing. Allelic heterogeneity at this single locus accounts for a large proportion of natural variation in vernalization that contributes to adaptation of A. thaliana.

Revisiting a GWAS peak in Arabidopsis thaliana reveals possible confounding by genetic heterogeneity.

Genome-wide association studies (GWAS) have become a standard approach for exploring the genetic basis of phenotypic variation. However, correlation is not causation, and only a tiny fraction of all associations have been experimentally confirmed. One practical problem is that a peak of association does not always pinpoint a causal gene, but may instead be tagging multiple causal variants. In this study, we reanalyze a previously reported peak associated with flowering time traits in Swedish Arabidopsis thaliana population. The peak appeared to pinpoint the AOP2/AOP3 cluster of glucosinolate biosynthesis genes, which is known to be responsible for natural variation in herbivore resistance. Here we propose an alternative hypothesis, by demonstrating that the AOP2/AOP3 flowering association can be wholly accounted for by allelic variation in two flanking genes with clear roles in regulating flowering: NDX1, a regulator of the main flowering time controller FLC, and GA1, which plays a central role in gibberellin synthesis and is required for flowering under some conditions. In other words, we propose that the AOP2/AOP3 flowering-time association may be yet another example of a spurious, "synthetic" association, arising from trying to fit a single-locus model in the presence of two statistically associated causative loci. We conclude that caution is needed when using GWAS for fine-mapping.

A genome-wide association study identifies variants underlying the Arabidopsis thaliana shade avoidance response.

Shade avoidance is an ecologically and molecularly well-understood set of plant developmental responses that occur when the ratio of red to far-red light (R:FR) is reduced as a result of foliar shade. Here, a genome-wide association study (GWAS) in Arabidopsis thaliana was used to identify variants underlying one of these responses: increased hypocotyl elongation. Four hypocotyl phenotypes were included in the study, including height in high R:FR conditions (simulated sun), height in low R:FR conditions (simulated shade), and two different indices of the response of height to low R:FR. GWAS results showed that variation in these traits is controlled by many loci of small to moderate effect. A known PHYC variant contributing to hypocotyl height variation was identified and lists of significantly associated genes were enriched in a priori candidates, suggesting that this GWAS was capable of generating meaningful results. Using metadata such as expression data, GO terms, and other annotation, we were also able to identify variants in candidate de novo genes. Patterns of significance among our four phenotypes allowed us to categorize associations into three groups: those that affected hypocotyl height without influencing shade avoidance, those that affected shade avoidance in a height-dependent fashion, and those that exerted specific control over shade avoidance. This grouping allowed for the development of explicit hypotheses about the genetics underlying shade avoidance variation. Additionally, the response to shade did not exhibit any marked geographic distribution, suggesting that variation in low R:FR-induced hypocotyl elongation may represent a response to local conditions.

Rapid creation of Arabidopsis doubled haploid lines for quantitative trait locus mapping.

Quantitative trait loci (QTL) mapping is a powerful tool for investigating the genetic basis of natural variation. QTL can be mapped using a number of different population designs, but recombinant inbred lines (RILs) are among the most effective. Unfortunately, homozygous RIL populations are time consuming to construct, typically requiring at least six generations of selfing starting from a heterozygous F(1). Haploid plants produced from an F(1) combine the two parental genomes and have only one allele at every locus. Converting these sterile haploids into fertile diploids (termed "doubled haploids," DHs) produces immortal homozygous lines in only two steps. Here we describe a unique technique for rapidly creating recombinant doubled haploid populations in Arabidopsis thaliana: centromere-mediated genome elimination. We generated a population of 238 doubled haploid lines that combine two parental genomes and genotyped them by reduced representation Illumina sequencing. The recombination rate and parental allele frequencies in our population are similar to those found in existing RIL sets. We phenotyped this population for traits related to flowering time and for petiole length and successfully mapped QTL controlling each trait. Our work demonstrates that doubled haploid populations offer a rapid, easy alternative to RILs for Arabidopsis genetic analysis.

Amino acid polymorphisms in Arabidopsis phytochrome B cause differential responses to light.

Plants have a sophisticated system for sensing and responding to their light environment. The light responses of populations and species native to different habitats show adaptive variation; understanding the mechanisms underlying photomorphogenic variation is therefore of significant interest. In Arabidopsis thaliana, phytochrome B (PHYB) is the dominant photoreceptor for red light and plays a major role in white light. Because PHYB has been proposed as a candidate gene for several quantitative trait loci (QTLs) affecting light response, we have investigated sequence and functional variation in Arabidopsis PHYB. We examined PHYB sequences in 33 A. thaliana individuals and in the close relative Arabidopsis lyrata. From 14 nonsynonymous polymorphisms, we chose 5 for further study based on previous QTL studies. In a larger collection of A. thaliana accessions, one of these five polymorphisms, I143L, was associated with variation in red light response. We used transgenic analysis to test this association and confirmed experimentally that natural PHYB polymorphisms cause differential plant responses to light. Furthermore, our results show that allelic variation of PHYB activity is due to amino acid rather than regulatory changes. Together with earlier studies linking variation in light sensitivity to photoreceptor genes, our work suggests that photoreceptors may be a common target of natural selection.

Ultra-High-Density QTL Marker Mapping for Seedling Photomorphogenesis Mediating Arabidopsis Establishment in Southern Patagonia.

Arabidopsis thaliana shows a wide range of genetic and trait variation among wild accessions. Because of its unparalleled biological and genomic resources, Arabidopsis has a high potential for the identification of genes underlying ecologically important complex traits, thus providing new insights on genome evolution. Previous research suggested that distinct light responses were crucial for Arabidopsis establishment in a peculiar ecological niche of southern Patagonia. The aim of this study was to explore the genetic basis of contrasting light-associated physiological traits that may have mediated the rapid adaptation to this new environment. From a biparental cross between the photomorphogenic contrasting accessions Patagonia (Pat) and Columbia (Col-0), we generated a novel recombinant inbred line (RIL) population, which was entirely next-generation sequenced to achieve ultra-high-density saturating molecular markers resulting in supreme mapping sensitivity. We validated the quality of the RIL population by quantitative trait loci (QTL) mapping for seedling de-etiolation, finding seven QTLs for hypocotyl length in the dark and continuous blue light (Bc), continuous red light (Rc), and continuous far-red light (FRc). The most relevant QTLs, Rc1 and Bc1, were mapped close together to chromosome V; the former for Rc and Rc/dark, and the latter for Bc, FRc, and dark treatments. The additive effects of both QTLs were confirmed by independent heterogeneous inbred families (HIFs), and we explored TZP and ABA1 as potential candidate genes for Rc1 and Bc1QTLs, respectively. We conclude that the Pat × Col-0 RIL population is a valuable novel genetic resource to explore other adaptive traits in Arabidopsis.

Molecular variation in a functionally divergent homolog of FCA regulates flowering time in Arabidopsis thaliana.

The identification and functional characterization of natural variants in plants are essential for understanding phenotypic adaptation. Here we identify a molecular variation in At2g47310 that contributes to the natural variation in flowering time in Arabidopsis thaliana accessions. This gene, which we term SISTER of FCA (SSF), functions in an antagonistic manner to its close homolog FCA. Genome-wide association analysis screens two major haplotypes of SSF associated with the natural variation in FLC expression, and a single polymorphism, SSF-N414D, is identified as a main contributor. The SSF414N protein variant interacts more strongly with CUL1, a component of the E3 ubiquitination complex, than the SSF414D form, mediating differences in SSF protein degradation and FLC expression. FCA and SSF appear to have arisen through gene duplication after dicot-monocot divergence, with the SSF-N414D polymorphism emerging relatively recently within A. thaliana. This work provides a good example for deciphering the functional importance of natural polymorphisms in different organisms.

The Aquilegia genome provides insight into adaptive radiation and reveals an extraordinarily polymorphic chromosome with a unique history.

The columbine genus Aquilegia is a classic example of an adaptive radiation, involving a wide variety of pollinators and habitats. Here we present the genome assembly of A. coerulea 'Goldsmith', complemented by high-coverage sequencing data from 10 wild species covering the world-wide distribution. Our analyses reveal extensive allele sharing among species and demonstrate that introgression and selection played a role in the Aquilegia radiation. We also present the remarkable discovery that the evolutionary history of an entire chromosome differs from that of the rest of the genome - a phenomenon that we do not fully understand, but which highlights the need to consider chromosomes in an evolutionary context.

The Generation of Doubled Haploid Lines for QTL Mapping.

Recombinant inbred lines (RILs) are an essential tool for quantitative trait locus (QTL) mapping in Arabidopsis thaliana. Conventionally, the development of these lines is a time-consuming and tedious process requiring six to eight generations of selfing. Here, we describe an alternative approach: the rapid generation of RILs in A. thaliana via the creation of doubled haploids. In this method, F1 plants are crossed to an engineered haploid inducer to produce haploid plants. The chromosomes of these haploids then spontaneously double, generating immortalized homozygous F2 lines called doubled haploid RILs (DH RILs). Finally, DH RILs are genotyped using low-coverage whole-genome sequencing and are ready to be used for QTL mapping.

Verification of Arabidopsis stock collections using SNPmatch, a tool for genotyping high-plexed samples.

Large-scale studies such as the Arabidopsis thaliana '1,001 Genomes' Project require routine genotyping of stocks to avoid sample contamination. To genotype samples efficiently and economically, sequencing must be inexpensive and data processing simple. Here we present SNPmatch, a tool that identifies strains (or inbred lines, or accessions) by matching them to a SNP database. We tested the tool by performing low-coverage resequencing of over 2,000 strains from our lab seed stock collection. SNPmatch correctly genotyped samples from 1-fold coverage sequencing data, and could also identify the parents of F1 or F2 individuals. SNPmatch can be run either on the command line or through AraGeno (https://arageno.gmi.oeaw.ac.at), a web interface that permits sample genotyping from a user-uploaded VCF or BED file.

Multiple alleles at a single locus control seed dormancy in Swedish Arabidopsis.

Seed dormancy is a complex life history trait that determines the timing of germination and is crucial for local adaptation. Genetic studies of dormancy are challenging, because the trait is highly plastic and strongly influenced by the maternal environment. Using a combination of statistical and experimental approaches, we show that multiple alleles at the previously identified dormancy locus DELAY OF GERMINATION1 jointly explain as much as 57% of the variation observed in Swedish Arabidopsis thaliana, but give rise to spurious associations that seriously mislead genome-wide association studies unless modeled correctly. Field experiments confirm that the major alleles affect germination as well as survival under natural conditions, and demonstrate that locally adaptive traits can sometimes be dissected genetically.

Anti-leukemia effect of perillyl alcohol in Bcr/Abl-transformed cells indirectly inhibits signaling through Mek in a Ras- and Raf-independent fashion.

PURPOSE: Perillyl alcohol (POH) displays preventive and therapeutic activity against a wide variety of tumor models, and it has been suggested that this might be associated with the ability of POH to interfere with Ras prenylation. POH also selectively induces G(1) arrest and apoptosis in Bcr/Abl-transformed hematopoietic cells. Because signaling through Ras is necessary for Bcr/Abl transformation, we examined whether POH induces its anti-leukemia effect by inhibiting Ras signaling. EXPERIMENTAL DESIGN: The ability of POH to inhibit posttranslational farnesylation and signaling from Ras as well as signaling through the Raf-Mek-Erk cascade was examined in Bcr/Abl-transformed and mock-transformed cells and related to the anti-leukemia effect of POH. RESULTS: POH does not affect Ras prenylation or Ras activity, but it blocks signaling downstream of Ras by reversing the state of activation of the Erk kinase, Mek. POH affects Mek activity only when it is added to intact cells. Treatment of either cell lysates or of purified Mek with POH has no effect on Mek activity. Inhibition of the Mek-Erk pathway seems to be related to the POH anti-leukemia effect for the following reasons: (a) the concentration of POH needed to block the Erk pathway, as well the kinetics with which POH inhibits this signaling cascade, both correlate with the anti-leukemia effect of POH; (b) both U0126 (a specific Mek inhibitor) and POH induce similar anti-leukemia effects; and (c) mock-transformed hematopoietic cells are simultaneously resistant to POH anti-leukemia effects and inhibition of the Mek-Erk pathway. CONCLUSION: Blocking Mek is sufficient to induce growth arrest and apoptosis in Bcr/Abl-transformed cells; therefore, POH represents a novel small molecule inhibitor of Mek that might be effective for treating Bcr/Abl leukemias.

Native environment modulates leaf size and response to simulated foliar shade across wild tomato species.

The laminae of leaves optimize photosynthetic rates by serving as a platform for both light capture and gas exchange, while minimizing water losses associated with thermoregulation and transpiration. Many have speculated that plants maximize photosynthetic output and minimize associated costs through leaf size, complexity, and shape, but a unifying theory linking the plethora of observed leaf forms with the environment remains elusive. Additionally, the leaf itself is a plastic structure, responsive to its surroundings, further complicating the relationship. Despite extensive knowledge of the genetic mechanisms underlying angiosperm leaf development, little is known about how phenotypic plasticity and selective pressures converge to create the diversity of leaf shapes and sizes across lineages. Here, we use wild tomato accessions, collected from locales with diverse levels of foliar shade, temperature, and precipitation, as a model to assay the extent of shade avoidance in leaf traits and the degree to which these leaf traits correlate with environmental factors. We find that leaf size is correlated with measures of foliar shade across the wild tomato species sampled and that leaf size and serration correlate in a species-dependent fashion with temperature and precipitation. We use far-red induced changes in leaf length as a proxy measure of the shade avoidance response, and find that shade avoidance in leaves negatively correlates with the level of foliar shade recorded at the point of origin of an accession. The direction and magnitude of these correlations varies across the leaf series, suggesting that heterochronic and/or ontogenic programs are a mechanism by which selective pressures can alter leaf size and form. This study highlights the value of wild tomato accessions for studies of both morphological and light-regulated development of compound leaves, and promises to be useful in the future identification of genes regulating potentially adaptive plastic leaf traits.