Genetic determinants of switchgrass-root-associated microbiota in field sites spanning its natural range.

A fundamental goal in plant microbiome research is to determine the relative impacts of host and environmental effects on root microbiota composition, particularly how host genotype impacts bacterial community composition. Most studies characterizing the effect of plant genotype on root microbiota undersample host genetic diversity and grow plants outside of their native ranges, making the associations between host and microbes difficult to interpret. Here, we characterized the root microbiota of a large diversity panel of switchgrass, a North American native C4 bioenergy crop, in three field locations spanning its native range. Our data, composed of 1,961 samples, suggest that field location is the primary determinant of microbiome composition; however, substantial heritable variation is widespread across bacterial taxa, especially those in the Sphingomonadaceae family. Despite diverse compositions, relatively few highly prevalent taxa make up the majority of the switchgrass root microbiota, a large fraction of which is shared across sites. Local genotypes preferentially recruit/filter for local microbes, supporting the idea of affinity between local plants and their microbiota. Using genome-wide association, we identified loci impacting the abundance of >400 microbial strains and found an enrichment of genes involved in immune responses, signaling pathways, and secondary metabolism. We found loci associated with over half of the core microbiota (i.e., microbes in >80% of samples), regardless of field location. Finally, we show a genetic relationship between a basal plant immunity pathway and relative abundances of root microbiota. This study brings us closer to harnessing and manipulating beneficial microbial associations via host genetics.

The genetic basis for panicle trait variation in switchgrass (Panicum virgatum).

KEY MESSAGE: We investigate the genetic basis of panicle architecture in switchgrass in two mapping populations across a latitudinal gradient, and find many stable, repeatable genetic effects and limited genetic interactions with the environment. Grass species exhibit large diversity in panicle architecture influenced by genes, the environment, and their interaction. The genetic study of panicle architecture in perennial grasses is limited. In this study, we evaluate the genetic basis of panicle architecture including panicle length, primary branching number, and secondary branching number in an outcrossed switchgrass QTL population grown across ten field sites in the central USA through multi-environment mixed QTL analysis. We also evaluate genetic effects in a diversity panel of switchgrass grown at three of the ten field sites using genome-wide association (GWAS) and multivariate adaptive shrinkage. Furthermore, we search for candidate genes underlying panicle traits in both of these independent mapping populations. Overall, 18 QTL were detected in the QTL mapping population for the three panicle traits, and 146 unlinked genomic regions in the diversity panel affected one or more panicle trait. Twelve of the QTL exhibited consistent effects (i.e., no QTL by environment interactions or no QTL × E), and most (four of six) of the effects with QTL × E exhibited site-specific effects. Most (59.3%) significant partially linked diversity panel SNPs had significant effects in all panicle traits and all field sites and showed pervasive pleiotropy and limited environment interactions. Panicle QTL co-localized with significant SNPs found using GWAS, providing additional power to distinguish between true and false associations in the diversity panel.

QTL × environment interactions underlie ionome divergence in switchgrass.

Ionomics measures elemental concentrations in biological organisms and provides a snapshot of physiology under different conditions. In this study, we evaluate genetic variation of the ionome in outbred, perennial switchgrass in three environments across the species' native range, and explore patterns of genotype-by-environment interactions. We grew 725 clonally replicated genotypes of a large full sib family from a four-way linkage mapping population, created from deeply diverged upland and lowland switchgrass ecotypes, at three common gardens. Concentrations of 18 mineral elements were determined in whole post-anthesis tillers using ion coupled plasma mass spectrometry (ICP-MS). These measurements were used to identify quantitative trait loci (QTL) with and without QTL-by-environment interactions (QTLxE) using a multi-environment QTL mapping approach. We found that element concentrations varied significantly both within and between switchgrass ecotypes, and GxE was present at both the trait and QTL level. Concentrations of 14 of the 18 elements were under some genetic control, and 77 QTL were detected for these elements. Seventy-four percent of QTL colocalized multiple elements, half of QTL exhibited significant QTLxE, and roughly equal numbers of QTL had significant differences in magnitude and sign of their effects across environments. The switchgrass ionome is under moderate genetic control and by loci with highly variable effects across environments.

Genomic mechanisms of climate adaptation in polyploid bioenergy switchgrass.

Long-term climate change and periodic environmental extremes threaten food and fuel security1 and global crop productivity2-4. Although molecular and adaptive breeding strategies can buffer the effects of climatic stress and improve crop resilience5, these approaches require sufficient knowledge of the genes that underlie productivity and adaptation6-knowledge that has been limited to a small number of well-studied model systems. Here we present the assembly and annotation of the large and complex genome of the polyploid bioenergy crop switchgrass (Panicum virgatum). Analysis of biomass and survival among 732 resequenced genotypes, which were grown across 10 common gardens that span 1,800 km of latitude, jointly revealed extensive genomic evidence of climate adaptation. Climate-gene-biomass associations were abundant but varied considerably among deeply diverged gene pools. Furthermore, we found that gene flow accelerated climate adaptation during the postglacial colonization of northern habitats through introgression of alleles from a pre-adapted northern gene pool. The polyploid nature of switchgrass also enhanced adaptive potential through the fractionation of gene function, as there was an increased level of heritable genetic diversity on the nondominant subgenome. In addition to investigating patterns of climate adaptation, the genome resources and gene-trait associations developed here provide breeders with the necessary tools to increase switchgrass yield for the sustainable production of bioenergy.

Genetic Associations in Four Decades of Multienvironment Trials Reveal Agronomic Trait Evolution in Common Bean.

Multienvironment trials (METs) are widely used to assess the performance of promising crop germplasm. Though seldom designed to elucidate genetic mechanisms, MET data sets are often much larger than could be duplicated for genetic research and, given proper interpretation, may offer valuable insights into the genetics of adaptation across time and space. The Cooperative Dry Bean Nursery (CDBN) is a MET for common bean (Phaseolus vulgaris) grown for > 70 years in the United States and Canada, consisting of 20-50 entries each year at 10-20 locations. The CDBN provides a rich source of phenotypic data across entries, years, and locations that is amenable to genetic analysis. To study stable genetic effects segregating in this MET, we conducted genome-wide association studies (GWAS) using best linear unbiased predictions derived across years and locations for 21 CDBN phenotypes and genotypic data (1.2 million SNPs) for 327 CDBN genotypes. The value of this approach was confirmed by the discovery of three candidate genes and genomic regions previously identified in balanced GWAS. Multivariate adaptive shrinkage (mash) analysis, which increased our power to detect significant correlated effects, found significant effects for all phenotypes. Mash found two large genomic regions with effects on multiple phenotypes, supporting a hypothesis of pleiotropic or linked effects that were likely selected on in pursuit of a crop ideotype. Overall, our results demonstrate that statistical genomics approaches can be used on MET phenotypic data to discover significant genetic effects and to define genomic regions associated with crop improvement.

Population Genetics of the Highly Polymorphic RPP8 Gene Family.

Plant nucleotide-binding domain and leucine-rich repeat containing (NLR) genes provide some of the most extreme examples of polymorphism in eukaryotic genomes, rivalling even the vertebrate major histocompatibility complex. Surprisingly, this is also true in Arabidopsis thaliana, a predominantly selfing species with low heterozygosity. Here, we investigate how gene duplication and intergenic exchange contribute to this extraordinary variation. RPP8 is a three-locus system that is configured chromosomally as either a direct-repeat tandem duplication or as a single copy locus, plus a locus 2 Mb distant. We sequenced 48 RPP8 alleles from 37 accessions of A. thaliana and 12 RPP8 alleles from Arabidopsis lyrata to investigate the patterns of interlocus shared variation. The tandem duplicates display fixed differences and share less variation with each other than either shares with the distant paralog. A high level of shared polymorphism among alleles at one of the tandem duplicates, the single-copy locus and the distal locus, must involve both classical crossing over and intergenic gene conversion. Despite these polymorphism-enhancing mechanisms, the observed nucleotide diversity could not be replicated under neutral forward-in-time simulations. Only by adding balancing selection to the simulations do they approach the level of polymorphism observed at RPP8. In this NLR gene triad, genetic architecture, gene function and selection all combine to generate diversity.

The genomic landscape of molecular responses to natural drought stress in Panicum hallii.

Environmental stress is a major driver of ecological community dynamics and agricultural productivity. This is especially true for soil water availability, because drought is the greatest abiotic inhibitor of worldwide crop yields. Here, we test the genetic basis of drought responses in the genetic model for C4 perennial grasses, Panicum hallii, through population genomics, field-scale gene-expression (eQTL) analysis, and comparison of two complete genomes. While gene expression networks are dominated by local cis-regulatory elements, we observe three genomic hotspots of unlinked trans-regulatory loci. These regulatory hubs are four times more drought responsive than the genome-wide average. Additionally, cis- and trans-regulatory networks are more likely to have opposing effects than expected under neutral evolution, supporting a strong influence of compensatory evolution and stabilizing selection. These results implicate trans-regulatory evolution as a driver of drought responses and demonstrate the potential for crop improvement in drought-prone regions through modification of gene regulatory networks.

Genetic architecture and pleiotropy shape costs of Rps2-mediated resistance in Arabidopsis thaliana.

The mounting evidence that R genes incur large fitness costs raises a question: how can there be a 5-10% fitness reduction for all 149 R genes in the Arabidopsis thaliana genome? The R genes tested to date segregate for insertion-deletion (indel) polymorphisms where susceptible alleles are complete deletions. Since costs of resistance are measured as the differential fitness of isolines carrying resistant and susceptible alleles, indels reveal costs that may be masked when susceptible alleles are expressed. Rps2 segregates for two expressed clades of alleles, one resistant and one susceptible. Plants with resistant Rps2 are not less fit than those with a susceptible Rps2 allele in the absence of disease. Instead, all alleles provide a fitness benefit relative to an artificial deletion because of the role of RPS2 as a negative regulator of defence. Our results highlight the interplay between genomic architecture and the magnitude of costs of resistance.

Modulation of R-gene expression across environments.

Some environments are more conducive to pathogen growth than others, and, as a consequence, plants might be expected to invest more in resistance when pathogen growth is favored. Resistance (R-) genes in Arabidopsis thaliana have unusually extensive variation in basal expression when comparing the same R-gene among accessions collected from different environments. R-gene expression variation was characterized to explore whether R-gene expression is up-regulated in environments favoring pathogen proliferation and down-regulated when risks of infection are low; down-regulation would follow if costs of R-gene expression negatively impact plant fitness in the absence of disease. Quantitative reverse transcription-PCR was used to quantify the expression of 13 R-gene loci in plants grown in eight environmental conditions for each of 12 A. thaliana accessions, and large effects of the environment on R-gene expression were found. Surprisingly, almost every change in the environment--be it a change in biotic or abiotic conditions--led to an increase in R-gene expression, a response that was distinct from the average transcriptome response and from that of other stress response genes. These changes in expression are functional in that environmental change prior to infection affected levels of specific disease resistance to isolates of Pseudomonas syringae. In addition, there are strong latitudinal clines in basal R-gene expression and clines in R-gene expression plasticity correlated with drought and high temperatures. These results suggest that variation in R-gene expression across environments may be shaped by natural selection to reduce fitness costs of R-gene expression in permissive or predictable environments.

Unique features of the m6A methylome in Arabidopsis thaliana.

Recent discoveries of reversible N(6)-methyladenosine (m(6)A) methylation on messenger RNA (mRNA) and mapping of m(6)A methylomes in mammals and yeast have revealed potential regulatory functions of this RNA modification. In plants, defects in m(6)A methyltransferase cause an embryo-lethal phenotype, suggesting a critical role of m(6)A in plant development. Here, we profile m(6)A transcriptome-wide in two accessions of Arabidopsis thaliana and reveal that m(6)A is a highly conserved modification of mRNA in plants. Distinct from mammals, m(6)A in A. thaliana is enriched not only around the stop codon and within 3'-untranslated regions, but also around the start codon. Gene ontology analysis indicates that the unique distribution pattern of m(6)A in A. thaliana is associated with plant-specific pathways involving the chloroplast. We also discover a positive correlation between m(6)A deposition and mRNA abundance, suggesting a regulatory role of m(6)A in plant gene expression.

Extensive loss of RNA editing sites in rapidly evolving Silene mitochondrial genomes: selection vs. retroprocessing as the driving force.

Theoretical arguments suggest that mutation rates influence the proliferation and maintenance of RNA editing. We identified RNA editing sites in five species within the angiosperm genus Silene that exhibit highly divergent mitochondrial mutation rates. We found that mutational acceleration has been associated with rapid loss of mitochondrial editing sites. In contrast, we did not find a significant difference in the frequency of editing in chloroplast genes, which lack the mutation rate variation observed in the mitochondrial genome. As found in other angiosperms, the rate of substitution at RNA editing sites in Silene greatly exceeds the rate at synonymous sites, a pattern that has previously been interpreted as evidence for selection against RNA editing. Alternatively, we suggest that editing sites may experience higher rates of C-to-T mutation than other portions of the genome. Such a pattern could be caused by gene conversion with reverse-transcribed mRNA (i.e., retroprocessing). If so, the genomic distribution of RNA editing site losses in Silene suggests that such conversions must be occurring at a local scale such that only one or two editing sites are affected at a time. Because preferential substitution at editing sites appears to occur in angiosperms regardless of the mutation rate, we conclude that mitochondrial rate accelerations within Silene have "fast-forwarded" a preexisting pattern but have not fundamentally changed the evolutionary forces acting on RNA editing sites.