Selection on soil microbiomes reveals reproducible impacts on plant function

Soil microorganisms found in the root zone impact plant growth and development, but the potential to harness these benefits is hampered by the sheer abundance and diversity of the players influencing desirable plant traits. Here, we report a high level of reproducibility of soil microbiomes in altering plant flowering time and soil functions when partnered within and between plant hosts. We used a multi-generation experimental system using Arabidopsis thaliana Col to select for soil microbiomes inducing earlier or later flowering times of their hosts. We then inoculated the selected microbiomes from the tenth generation of plantings into the soils of three additional A. thaliana genotypes (Ler, Be, RLD) and a related crucifer (Brassica rapa). With the exception of Ler, all other plant hosts showed a shift in flowering time corresponding with the inoculation of early- or late-flowering microbiomes. Analysis of the soil microbial community using 16 S rRNA gene sequencing showed distinct microbiota profiles assembling by flowering time treatment. Plant hosts grown with the late-flowering-associated microbiomes showed consequent increases in inflorescence biomass for three A. thaliana genotypes and an increase in total biomass for B. rapa. The increase in biomass was correlated with two- to five-fold enhancement of microbial extracellular enzyme activities associated with nitrogen mineralization in soils. The reproducibility of the flowering phenotype across plant hosts suggests that microbiomes can be selected to modify plant traits and coordinate changes in soil resource pools.     

Global Genetic Heterogeneity in Adaptive Traits

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

Genome-Wide Association Analysis of Adaptation Using Environmentally Predicted Traits

Current methods for studying the genetic basis of adaptation evaluate genetic associations with ecologically relevant traits or single environmental variables, under the implicit assumption that natural selection imposes correlations between phenotypes, environments and genotypes. In practice, observed trait and environmental data are manifestations of unknown selective forces and are only indirectly associated with adaptive genetic variation. In theory, improved estimation of these forces could enable more powerful detection of loci under selection. Here we present an approach in which we approximate adaptive variation by modeling phenotypes as a function of the environment and using the predicted trait in multivariate and univariate genome-wide association analysis (GWAS). Based on computer simulations and published flowering time data from the model plant Arabidopsis thaliana, we find that environmentally predicted traits lead to higher recovery of functional loci in multivariate GWAS and are more strongly correlated to allele frequencies at adaptive loci than individual environmental variables. Our results provide an example of the use of environmental data to obtain independent and meaningful information on adaptive genetic variation.     

Genome-wide association mapping in Arabidopsis identifies previously known flowering time and pathogen resistance genes

There is currently tremendous interest in the possibility of using genome-wide association mapping to identify genes responsible for natural variation, particularly for human disease susceptibility. The model plant Arabidopsis thaliana is in many ways an ideal candidate for such studies, because it is a highly selfing hermaphrodite. As a result, the species largely exists as a collection of naturally occurring inbred lines, or accessions, which can be genotyped once and phenotyped repeatedly. Furthermore, linkage disequilibrium in such a species will be much more extensive than in a comparable outcrossing species. We tested the feasibility of genome-wide association mapping in A. thaliana by searching for associations with flowering time and pathogen resistance in a sample of 95 accessions for which genome-wide polymorphism data were available. In spite of an extremely high rate of false positives due to population structure, we were able to identify known major genes for all phenotypes tested, thus demonstrating the potential of genome-wide association mapping in A. thaliana and other species with similar patterns of variation. The rate of false positives differed strongly between traits, with more clinal traits showing the highest rate. However, the false positive rates were always substantial regardless of the trait, highlighting the necessity of an appropriate genomic control in association studies.     

Direct and indirect selection on flowering time,water‐use efficiency (WUE, δ13C), and WUE plasticity to drought in Arabidopsis thaliana

Flowering time and water-use efficiency (WUE) are two ecological traits that are important for plant drought response. To understand the evolutionary significance of natural genetic variation in flowering time, WUE, and WUE plasticity to drought in Arabidopsis thaliana, we addressed the following questions: (1) How are ecophysiological traits genetically correlated within and between different soil moisture environments? (2) Does terminal drought select for early flowering and drought escape? (3) Is WUE plasticity to drought adaptive and/or costly? We measured a suite of ecophysiological and reproductive traits on 234 spring flowering accessions of A. thaliana grown in well-watered and season-ending soil drying treatments, and quantified patterns of genetic variation, correlation, and selection within each treatment. WUE and flowering time were consistently positively genetically correlated. WUE was correlated with WUE plasticity, but the direction changed between treatments. Selection generally favored early flowering and low WUE, with drought favoring earlier flowering significantly more than well-watered conditions. Selection for lower WUE was marginally stronger under drought. There were no net fitness costs of WUE plasticity. WUE plasticity (per se) was globally neutral, but locally favored under drought. Strong genetic correlation between WUE and flowering time may facilitate the evolution of drought escape, or constrain independent evolution of these traits. Terminal drought favored drought escape in these spring flowering accessions of A. thaliana. WUE plasticity may be favored over completely fixed development in environments with periodic drought.     

QTL mapping of naturally-occurring variation in flowering time of Arabidopsis thaliana

A segregating F₂ population of Arabidopsis thaliana derived from a cross between the late-flowering ecotype Hannover/Münden (HM) and the early-flowering ecotype Wassilewskija (WS) was analyzed for flowering time and other morphological traits. Two unlinked quantitative trait loci (QTLs) affecting days to first flower (DFF-a and DFF-b) mapped to chromosome 5. QTLs which affect node number (NN), leaf length at flowering (LLF), and leaf length at 35 days (LL35) also mapped to chromosome 5; LLF-a, LL35-a, NN-a map to the same region of chromosome 5 as DFF-a; LLF-b and LL35-bmap to the same region of chromosome 5 as DFF-b. Another QTL affecting leaf length at flowering (LLF-c) maps to chromosome 3. The proximity of DFF-a, LLF-a, LL35-a and NN-a, as well as the similarity in gene action among these QTLs (additivity), suggest that they may be pleiotropic consequences of a single gene at this locus. Similarly, LL35-b and LLF-b map near each other and both display recessive gene action, again suggesting the possibility of pleiotropy. DFF-b, which also maps near LL35-b and LLF-b, displays largely additive gene action (although recessive gene action could not be ruled out). This suggests that DFF-b may represent a different gene from LL35-b and/or LLF-b. DFF-a maps near two previously identified mutants: co (which also affects flowering time and displays gene action consistent with additivity) and flc. Similar map locations and gene actions of QTLs affecting the correlated traits DFF, LLF, LL35 and NN suggest that these genomic regions harbor naturally occurring allelic variants involved in the general transition of the plant from vegetative to reproductive growth.     

Demographic and Genetic Patterns of Variation among Populations of Arabidopsis thaliana from Contrasting Native Environments

Background

Understanding the relationship between environment and genetics requires the integration of knowledge on the demographic behavior of natural populations. However, the demographic performance and genetic composition of Arabidopsis thaliana populations in the species'' native environments remain largely uncharacterized. This information, in combination with the advances on the study of gene function, will improve our understanding on the genetic mechanisms underlying adaptive evolution in A. thaliana.

Methodology/Principal Findings

We report the extent of environmental, demographic, and genetic variation among 10 A. thaliana populations from Mediterranean (coastal) and Pyrenean (montane) native environments in northeast Spain. Geographic, climatic, landscape, and soil data were compared. Demographic traits, including the dynamics of the soil seed bank and the attributes of aboveground individuals followed over a complete season, were also analyzed. Genetic data based on genome-wide SNP markers were used to describe genetic diversity, differentiation, and structure. Coastal and montane populations significantly differed in terms of environmental, demographic, and genetic characteristics. Montane populations, at higher altitude and farther from the sea, are exposed to colder winters and prolonged spring moisture compared to coastal populations. Montane populations showed stronger secondary seed dormancy, higher seedling/juvenile mortality in winter, and initiated flowering later than coastal populations. Montane and coastal regions were genetically differentiated, montane populations bearing lower genetic diversity than coastal ones. No significant isolation-by-distance pattern and no shared multilocus genotypes among populations were detected.

Conclusions/Significance

Between-region variation in climatic patterns can account for differences in demographic traits, such as secondary seed dormancy, plant mortality, and recruitment, between coastal and montane A. thaliana populations. In addition, differences in plant mortality can partly account for differences in the genetic composition of coastal and montane populations. This study shows how the interplay between variation in environmental, demographic, and genetic parameters may operate in natural A. thaliana populations.     

Pleiotropy of FRIGIDA enhances the potential for multivariate adaptation

An evolutionary response to selection requires genetic variation; however, even if it exists, then the genetic details of the variation can constrain adaptation. In the simplest case, unlinked loci and uncorrelated phenotypes respond directly to multivariate selection and permit unrestricted paths to adaptive peaks. By contrast, ‘antagonistic’ pleiotropic loci may constrain adaptation by affecting variation of many traits and limiting the direction of trait correlations to vectors that are not favoured by selection. However, certain pleiotropic configurations may improve the conditions for adaptive evolution. Here, we present evidence that the Arabidopsis thaliana gene FRI (FRIGIDA) exhibits ‘adaptive’ pleiotropy, producing trait correlations along an axis that results in two adaptive strategies. Derived, low expression FRI alleles confer a ‘drought escape’ strategy owing to fast growth, low water use efficiency and early flowering. By contrast, a dehydration avoidance strategy is conferred by the ancestral phenotype of late flowering, slow growth and efficient water use during photosynthesis. The dehydration avoidant phenotype was recovered when genotypes with null FRI alleles were transformed with functional alleles. Our findings indicate that the well-documented effects of FRI on phenology result from differences in physiology, not only a simple developmental switch.     

Foliar fungi alter reproductive timing and allocation in Arabidopsis under normal and water-stressed conditions

Microbes influence plant phenotypes but most known examples of this are from the study of below-ground microbes and plant disease modification. To examine the potential importance of phyllosphere microbes on non-disease related plant traits, we used sterile Arabidopsis clones to test the effects of foliar fungi on flowering phenology and reproductive allocation under conditions of varying water stress. We inoculated the sterile plants with fully-factorial combinations of four fungal isolates, then measured flowering time and reproductive allocation for each treatment group under normal and water-stressed conditions. All plants inoculated with foliar fungi had significantly later flowering and greater seed mass than the sterile control groups. The magnitude of this effect depended on the specific fungi present, but individual fungal effects diminished as inoculum richness increased. Above-ground microbes likely influence other plant traits as well and should be considered in any study measuring plant phenotypes.     

Quantitative Trait Loci (QTL) and Mendelian Trait Loci (MTL) Analysis in Prunus: a Breeding Perspective and Beyond

Trait loci analysis, a classic procedure in quantitative (quantitative trait loci, QTL) and qualitative (Mendelian trait loci, MTL) genetics, continues to be the most important approach in studies of gene labeling in Prunus species from the Rosaceae family. Since 2011, the number of published Prunus QTLs and MTLs has doubled. With increased genomic resources, such as whole genome sequences and high-density genotyping platforms, trait loci analysis can be more readily converted to markers that can be directly utilized in marker-assisted breeding. To provide this important resource to the community and to integrate it with other genomic, genetic, and breeding data, a global review of the QTLs and MTLs linked to agronomic traits in Prunus has been performed and the data made available in the Genome Database for Rosaceae. We describe detailed information on 760 main QTLs and MTLs linked to a total of 110 agronomic traits related to tree development, pest and disease resistance, flowering, ripening, and fruit and seed quality. Access to these trait loci enables the application of this information in the post-genomic era, characterized by the availability of a high-quality peach reference genome and new high-throughput DNA and RNA analysis technologies.     

QTL analysis of flowering time and ripening traits suggests an impact of a genomic region on linkage group 1 in Vitis

In the recent past, genetic analyses of grapevine focused mainly on the identification of resistance loci for major diseases such as powdery and downy mildew. Currently, breeding programs make intensive use of these results by applying molecular markers linked to the resistance traits. However, modern genetics also allows to address additional agronomic traits that have considerable impact on the selection of grapevine cultivars. In this study, we have used linkage mapping for the identification and characterization of flowering time and ripening traits in a mapping population from a cross of V3125 (‘Schiava Grossa’ × ‘Riesling’) and the interspecific rootstock cultivar ‘Börner’ (Vitis riparia × Vitis cinerea). Comparison of the flowering time QTL mapping with data derived from a second independent segregating population identified several common QTLs. Especially a large region on linkage group 1 proved to be of special interest given the genetic divergence of the parents of the two populations. The proximity of the QTL region contains two CONSTANS-like genes. In accordance with data from other plants such as Arabidopsis thaliana and Oryza sativa, we hypothesize that these genes are major contributors to control the time of flowering in Vitis.     

Mass spectrometry analysis of the variants of histone H3 and H4 of soybean and their post-translational modifications

Background

Functional genomics tools provide researchers with the ability to apply high-throughput techniques to determine the function and interaction of a diverse range of genes. Mutagenised plant populations are one such resource that facilitate gene characterisation. They allow complex physiological responses to be correlated with the expression of single genes in planta, through either reverse genetics where target genes are mutagenised to assay the affect, or through forward genetics where populations of mutant lines are screened to identify those whose phenotype diverges from wild type for a particular trait. One limitation of these types of populations is the prevalence of gene redundancy within plant genomes, which can mask the affect of individual genes. Activation or enhancer populations, which not only provide knock-out but also dominant activation mutations, can facilitate the study of such genes.

Results

We have developed a population of almost 50,000 activation tagged A. thaliana lines that have been archived as individual lines to the T₃ generation. The population is an excellent tool for both reverse and forward genetic screens and has been used successfully to identify a number of novel mutants. Insertion site sequences have been generated and mapped for 15,507 lines to enable further application of the population, while providing a clear distribution of T-DNA insertions across the genome. The population is being screened for a number of biochemical and developmental phenotypes, provisional data identifying novel alleles and genes controlling steps in proanthocyanidin biosynthesis and trichome development is presented.

Conclusion

This publicly available population provides an additional tool for plant researcher's to assist with determining gene function for the many as yet uncharacterised genes annotated within the Arabidopsis genome sequence http://aafc-aac.usask.ca/FST. The presence of enhancer elements on the inserted T-DNA molecule allows both knock-out and dominant activation phenotypes to be identified for traits of interest.     

Next-generation sequencing as a tool to quickly identify causative EMS-generated mutations

The advent of next generation sequencing has influenced every aspect of biological research. Many labs are now using whole genome sequencing in Arabidopsis thaliana as a means to quickly identify EMS-generated mutations present in isolated mutants. Following identification of these mutations, examination of T-DNA insertional alleles defective in candidate genes or complementation of the mutant phenotype with a wild type copy of candidate genes can be used to verify which mutation is causative for the phenotype of interest. Here, we discuss the benefits and pitfalls of using this method to identify mutations underlying phenotypes.     

Identification of two phenotypes of <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis> under <Emphasis Type="Italic">in vitro</Emphasis> salt stress conditions

This study describes two phenotypes of Arabidopsis thaliana (ecotype Columbia) developed in vitro under salt stress (75 mM NaCl). The phenotypes 01 and 02 appeared visibly distinguishable by rosette morphology and competence to produce flowers. Phenotype 01, sensible to salt stress, accumulated high quantities of Na⁺, showed a slight reduction in dry mass, and high protein and chlorophyll contents. Moreover, its anatomy exhibited some xeromorphic traits. Phenotype 02, clearly salt tolerant, showed a morphology similar to control plants, displaying typical phyllotactic rosette and flowering stalk production. Accumulation of Na⁺, protein and chlorophyll contents were close to control plants. Reversion experiments on NaCl free MS medium, showed a partially recovered phenotype 01. A threshold salt stress concentration that permits the simultaneous development of two phenotypes, was found.     

Disruption of the Homogentisate Solanesyltransferase Gene Results in Albino and Dwarf Phenotypes and Root,Trichome and Stomata Defects in Arabidopsis thaliana

Homogentisate solanesyltransferase (HST) plays an important role in plastoquinone (PQ) biosynthesis and acts as the electron acceptor in the carotenoids and abscisic acid (ABA) biosynthesis pathways. We isolated and identified a T-DNA insertion mutant of the HST gene that displayed the albino and dwarf phenotypes. PCR analyses and functional complementation also confirmed that the mutant phenotypes were caused by disruption of the HST gene. The mutants also had some developmental defects, including trichome development and stomata closure defects. Chloroplast development was also arrested and chlorophyll (Chl) was almost absent. Developmental defects in the chloroplasts were consistent with the SDS-PAGE result and the RNAi transgenic phenotype. Exogenous gibberellin (GA) could partially rescue the dwarf phenotype and the root development defects and exogenous ABA could rescue the stomata closure defects. Further analysis showed that ABA and GA levels were both very low in the pds2-1 mutants, which suggested that biosynthesis inhibition by GAs and ABA contributed to the pds2-1 mutants'' phenotypes. An early flowering phenotype was found in pds2-1 mutants, which showed that disruption of the HST gene promoted flowering by partially regulating plant hormones. RNA-sequencing showed that disruption of the HST gene resulted in expression changes to many of the genes involved in flowering time regulation and in the biosynthesis of PQ, Chl, GAs, ABA and carotenoids. These results suggest that HST is essential for chloroplast development, hormone biosynthesis, pigment accumulation and plant development.     

A stress-inducible protein regulates drought tolerance and flowering time in Brachypodium and Arabidopsis

To cope with environmental stresses and ensure maximal reproductive success, plants have developed strategies to adjust the timing of their transition to reproductive growth. This has a substantial impact on the stress resilience of crops and ultimately on agricultural productivity. Here, we report a previously uncharacterized, plant-specific gene family designated as Regulator of Flowering and Stress (RFS). Overexpression of the BdRFS gene in Brachypodium distachyon delayed flowering, increased biomass accumulation, and promoted drought tolerance, whereas clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9)-mediated knockout mutants exhibited opposite phenotypes. A double T-DNA insertional mutant in the two Arabidopsis (Arabidopsis thaliana) homologs replicated the effects on flowering and water deprivation seen in the B. distachyon CRISPR knockout lines, highlighting the functional conservation of the family between monocots and dicots. Lipid analysis of B. distachyon and Arabidopsis revealed that digalactosyldiacylglycerol (DGDG) and phosphatidylcholine (PC) contents were significantly, and reciprocally, altered in overexpressor and knockout mutants. Importantly, alteration of C16:0-containing PC, a Flowering Locus T-interacting lipid, associated with flowering phenotype, with elevated levels corresponding to earlier flowering. Co-immunoprecipitation analysis suggested that BdRFS interacts with phospholipase Dα1 as well as several other abscisic acid-related proteins. Furthermore, reduction of C18:3 fatty acids in DGDG corresponded with reduced jasmonic acid metabolites in CRISPR mutants. Collectively, we suggest that stress-inducible RFS proteins represent a regulatory component of lipid metabolism that impacts several agronomic traits of biotechnological importance.

A novel protein family modulates phospholipid metabolism and impacts multiple traits of biotechnological significance, likely through interaction with phospholipase and abscisic acid-related proteins.     

Raman imaging of cell wall polymers in Arabidopsis thaliana

We present chemical images of Arabidopsis thaliana stem cross-sections acquired by confocal Raman microscopy. Using green light (532 nm) from a continuous wave laser, the spatial distributions of cell wall polymers in Arabidopsis are visualized for the first time with lateral resolution that is sub-μm. Our results facilitate the label-free in situ characterization and screening of cell wall composition in this plant biology and genetics model organism, contributing ultimately towards an understanding of the molecular biology of many plant traits.     

Mutational effects on constraints on character evolution and phenotypic plasticity inArabidopsis thaliana

Although the concept of genetic constraints plays an important role in our understanding of the evolution of natural populations, there are still few empirical investigations probing the nature and limits of constraints in plant and animal species, aside from some studies inDrosophila. In the work reported here, we use an induced mutation - artificial selection protocol to analyse constraints on character means and phenotypic plasticity to nutrients inArabidopsis thaliana, an annual crucifer. We induced point mutations in a highly inbred line characterized by an extreme phenotype (very fast life cycle, early flowering, reduced leaf production) and little plasticity. We then selected individuals with increased leaf numbers. The goals were to determine if: (i) it is possible to increase leaf production; (ii) this has an effect on reproductive fitness; (iii) a mutation-selection process simultaneously alters the environmental insensitivity of the plant, thereby allowing phenotypic plasticity; and (iv) changes in the target trait affect other characters or their plasticities. The results demonstrate that: (a) mutations do increase leaf number; (b) this yields a much higher reproductive fitness, owing to the extension of the very short life cycle of the base inbred line; (c) there are no changes in plasticity of leaf number or of any other trait, possibly because few loci are involved in the control of plasticity; (d) changes in leaf number are related to alterations in three other traits comprising a strong set of covarying characters inA. thaliana. Two uncorrelated traits are capable of independent evolution from the constrained set. We therefore suggest that environmentally insensitive ecotypes of A.thaliana can quickly evolve to form ecologically specialized, relatively environmentally invariant genotypes.     

Combining genome-wide association mapping and transcriptional networks to identify novel genes controlling glucosinolates in Arabidopsis thaliana

Background

Genome-wide association (GWA) is gaining popularity as a means to study the architecture of complex quantitative traits, partially due to the improvement of high-throughput low-cost genotyping and phenotyping technologies. Glucosinolate (GSL) secondary metabolites within Arabidopsis spp. can serve as a model system to understand the genomic architecture of adaptive quantitative traits. GSL are key anti-herbivory defenses that impart adaptive advantages within field trials. While little is known about how variation in the external or internal environment of an organism may influence the efficiency of GWA, GSL variation is known to be highly dependent upon the external stresses and developmental processes of the plant lending it to be an excellent model for studying conditional GWA.

Methodology/Principal Findings

To understand how development and environment can influence GWA, we conducted a study using 96 Arabidopsis thaliana accessions, >40 GSL phenotypes across three conditions (one developmental comparison and one environmental comparison) and ∼230,000 SNPs. Developmental stage had dramatic effects on the outcome of GWA, with each stage identifying different loci associated with GSL traits. Further, while the molecular bases of numerous quantitative trait loci (QTL) controlling GSL traits have been identified, there is currently no estimate of how many additional genes may control natural variation in these traits. We developed a novel co-expression network approach to prioritize the thousands of GWA candidates and successfully validated a large number of these genes as influencing GSL accumulation within A. thaliana using single gene isogenic lines.

Conclusions/Significance

Together, these results suggest that complex traits imparting environmentally contingent adaptive advantages are likely influenced by up to thousands of loci that are sensitive to fluctuations in the environment or developmental state of the organism. Additionally, while GWA is highly conditional upon genetics, the use of additional genomic information can rapidly identify causal loci en masse.