The evolutionary genomics of adaptation to stress in wild rhizobium bacteria.

Microbiota comprise the bulk of life's diversity, yet we know little about how populations of microbes accumulate adaptive diversity across natural landscapes. Adaptation to stressful soil conditions in plants provides seminal examples of adaptation in response to natural selection via allelic substitution. For microbes symbiotic with plants however, horizontal gene transfer allows for adaptation via gene gain and loss, which could generate fundamentally different evolutionary dynamics. We use comparative genomics and genetics to elucidate the evolutionary mechanisms of adaptation to physiologically stressful serpentine soils in rhizobial bacteria in western North American grasslands. In vitro experiments demonstrate that the presence of a locus of major effect, the nre operon, is necessary and sufficient to confer adaptation to nickel, a heavy metal enriched to toxic levels in serpentine soil, and a major axis of environmental soil chemistry variation. We find discordance between inferred evolutionary histories of the core genome and nreAXY genes, which often reside in putative genomic islands. This suggests that the evolutionary history of this adaptive variant is marked by frequent losses, and/or gains via horizontal acquisition across divergent rhizobium clades. However, different nre alleles confer distinct levels of nickel resistance, suggesting allelic substitution could also play a role in rhizobium adaptation to serpentine soil. These results illustrate that the interplay between evolution via gene gain and loss and evolution via allelic substitution may underlie adaptation in wild soil microbiota. Both processes are important to consider for understanding adaptive diversity in microbes and improving stress-adapted microbial inocula for human use.

Negotiating mutualism: A locus for exploitation by rhizobia has a broad effect size distribution and context-dependent effects on legume hosts.

In mutualisms, variation at genes determining partner fitness provides the raw material upon which coevolutionary selection acts, setting the dynamics and pace of coevolution. However, we know little about variation in the effects of genes that underlie symbiotic fitness in natural mutualist populations. In some species of legumes that form root nodule symbioses with nitrogen-fixing rhizobial bacteria, hosts secrete nodule-specific cysteine-rich (NCR) peptides that cause rhizobia to differentiate in the nodule environment. However, rhizobia can cleave NCR peptides through the expression of genes like the plasmid-borne Host range restriction peptidase (hrrP), whose product degrades specific NCR peptides. Although hrrP activity can confer host exploitation by depressing host fitness and enhancing symbiont fitness, the effects of hrrP on symbiosis phenotypes depend strongly on the genotypes of the interacting partners. However, the effects of hrrP have yet to be characterised in a natural population context, so its contribution to variation in wild mutualist populations is unknown. To understand the distribution of effects of hrrP in wild rhizobia, we measured mutualism phenotypes conferred by hrrP in 12 wild Ensifer medicae strains. To evaluate context dependency of hrrP effects, we compared hrrP effects across two Medicago polymorpha host genotypes and across two experimental years for five E. medicae strains. We show for the first time in a natural population context that hrrP has a wide distribution of effect sizes for many mutualism traits, ranging from strongly positive to strongly negative. Furthermore, we show that hrrP effect size varies across host genotypes and experiment years, suggesting that researchers should be cautious about extrapolating the role of genes in natural populations from controlled laboratory studies of single genetic variants.

Multiple Mutualism Effects generate synergistic selection and strengthen fitness alignment in the interaction between legumes, rhizobia and mycorrhizal fungi.

Nearly all organisms participate in multiple mutualisms, and complementarity within these complex interactions can result in synergistic fitness effects. However, it remains largely untested how multiple mutualisms impact eco-evolutionary dynamics in interacting species. We tested how multiple microbial mutualists-N-fixing bacteria and mycorrrhizal fungi-affected selection and heritability of traits in their shared host plant (Medicago truncatula), as well as fitness alignment between partners. Our results demonstrate for the first time that multiple mutualisms synergistically affect the selection and heritability of host traits and enhance fitness alignment between mutualists. Specifically, we found interaction with multiple microbial symbionts doubled the strength of natural selection on a plant architectural trait, resulted in 2- to 3-fold higher heritability of plant reproductive success, and more than doubled fitness alignment between N-fixing bacteria and plants. These findings show synergism generated by multiple mutualisms extends to key components of microevolutionary change, emphasising the importance of multiple mutualism effects on evolutionary trajectories.

Whose trait is it anyways? Coevolution of joint phenotypes and genetic architecture in mutualisms.

Evolutionary biologists typically envision a trait's genetic basis and fitness effects occurring within a single species. However, traits can be determined by and have fitness consequences for interacting species, thus evolving in multiple genomes. This is especially likely in mutualisms, where species exchange fitness benefits and can associate over long periods of time. Partners may experience evolutionary conflict over the value of a multi-genomic trait, but such conflicts may be ameliorated by mutualism's positive fitness feedbacks. Here, we develop a simulation model of a host-microbe mutualism to explore the evolution of a multi-genomic trait. Coevolutionary outcomes depend on whether hosts and microbes have similar or different optimal trait values, strengths of selection and fitness feedbacks. We show that genome-wide association studies can map joint traits to loci in multiple genomes and describe how fitness conflict and fitness feedback generate different multi-genomic architectures with distinct signals around segregating loci. Partner fitnesses can be positively correlated even when partners are in conflict over the value of a multi-genomic trait, and conflict can generate strong mutualistic dependency. While fitness alignment facilitates rapid adaptation to a new optimum, conflict maintains genetic variation and evolvability, with implications for applied microbiome science.

Electrochemically Active Biofilms as an Indicator of Soil Health.

Soil health is a complex phenomenon that reflects the ability of soil to support both plant growth and other ecosystem functions. To our knowledge, research on extracellular electron transfer processes in soil environments is limited and could provide novel knowledge and new ways of monitoring soil health. Electrochemical activities in the soil can be studied by inserting inert electrodes. Once the electrode is polarized to a favorable potential, nearby microorganisms attach to the electrodes and grow as biofilms. Biofilms are a major part of the soil and play critical roles in microbial activity and community dynamics. Our work aims to investigate the electrochemical behavior of healthy and unhealthy soils using chronoamperometry and cyclic voltammetry. We developed a bioelectrochemical soil reactor for electrochemical measurements using healthy and unhealthy soils taken from the Cook Agronomy Farm Long-Term Agroecological Research site; the soils showed similar physical and chemical characteristics, but there was higher plant growth where the healthy soil was taken. Using carbon cloth electrodes installed in these soil reactors, we explored the electrochemical signals in these two soils. First, we measured redox variations by depth and found that reducing conditions were prevalent in healthy soils. Current measurements showed distinct differences between healthy and unhealthy soils. Scanning electron microscopy images showed the presence of microbes attached to the electrode for healthy soil but not for unhealthy soil. Glucose addition stimulated current in both soil types and caused differences in cyclic voltammograms between the two soil types to converge. Our work demonstrates that we can use current as a proxy for microbial metabolic activity to distinguish healthy and unhealthy soil.

Research priorities for harnessing plant microbiomes in sustainable agriculture.

Feeding a growing world population amidst climate change requires optimizing the reliability, resource use, and environmental impacts of food production. One way to assist in achieving these goals is to integrate beneficial plant microbiomes-i.e., those enhancing plant growth, nutrient use efficiency, abiotic stress tolerance, and disease resistance-into agricultural production. This integration will require a large-scale effort among academic researchers, industry researchers, and farmers to understand and manage plant-microbiome interactions in the context of modern agricultural systems. Here, we identify priorities for research in this area: (1) develop model host-microbiome systems for crop plants and non-crop plants with associated microbial culture collections and reference genomes, (2) define core microbiomes and metagenomes in these model systems, (3) elucidate the rules of synthetic, functionally programmable microbiome assembly, (4) determine functional mechanisms of plant-microbiome interactions, and (5) characterize and refine plant genotype-by-environment-by-microbiome-by-management interactions. Meeting these goals should accelerate our ability to design and implement effective agricultural microbiome manipulations and management strategies, which, in turn, will pay dividends for both the consumers and producers of the world food supply.

A Genome Resource for the Apple Powdery Mildew Pathogen Podosphaera leucotricha.

Powdery mildew, caused by Podosphaera leucotricha, is an economically important disease of apple and pear trees. A single monoconidial strain (PuE-3) of this biotrophic fungus was used to extract DNA for Illumina sequencing. Data were assembled to form a draft genome of 43.8 Mb consisting of 8,921 contigs, 9,372 predicted genes, and 96.1% of complete benchmarking universal single copy orthologs (BUSCOs). This is the first reported genome sequence of P. leucotricha that will enable studies of the population biology, epidemiology, and fungicide resistance of this pathogen. Furthermore, this resource will be fundamental to uncover the genetic and molecular mechanisms of the apple-powdery mildew interaction, and support future pome fruit breeding efforts.

Identifying similar transcripts in a related organism from de Bruijn graphs of RNA-Seq data, with applications to the study of salt and waterlogging tolerance in Melilotus.

BACKGROUND: A popular strategy to study alternative splicing in non-model organisms starts from sequencing the entire transcriptome, then assembling the reads by using de novo transcriptome assembly algorithms to obtain predicted transcripts. A similarity search algorithm is then applied to a related organism to infer possible function of these predicted transcripts. While some of these predictions may be inaccurate and transcripts with low coverage are often missed, we observe that it is possible to obtain a more complete set of transcripts to facilitate possible functional assignments by starting the search from the intermediate de Bruijn graph that contains all branching possibilities. RESULTS: We develop an algorithm to extract similar transcripts in a related organism by starting the search from the de Bruijn graph that represents the transcriptome instead of from predicted transcripts. We show that our algorithm is able to recover more similar transcripts than existing algorithms, with large improvements in obtaining longer transcripts and a finer resolution of isoforms. We apply our algorithm to study salt and waterlogging tolerance in two Melilotus species by constructing new RNA-Seq libraries. CONCLUSIONS: We have developed an algorithm to identify paths in the de Bruijn graph that correspond to similar transcripts in a related organism directly. Our strategy bypasses the transcript prediction step in RNA-Seq data and makes use of support from evolutionary information.

Mutualists Stabilize the Coexistence of Congeneric Legumes.

Coexistence requires that stabilizing niche differences, which cause species to limit themselves more than others, outweigh relative fitness differences, which cause competitive exclusion. Interactions with shared mutualists, which can differentially affect host fitness and change in magnitude with host frequency, can satisfy these conditions for coexistence, yet empirical tests of mutualist effects on relative fitness and stabilizing niche differences are largely lacking within the framework of coexistence theory. Here, we show that N-fixing rhizobial mutualists mediate coexistence in four naturally co-occurring congeneric legume (Trifolium) species. Using experimental greenhouse communities, we quantified relative fitness and stabilizing niche differences for each species in the presence of rhizobia originating from conspecific or congeneric hosts. Rhizobia stabilized coexistence by increasing the self-limitation of Trifolium species grown with rhizobia isolated from conspecifics, thus allowing congeners to increase when rare. Greenhouse-measured invasion growth rates predicted natural, unmanipulated coexistence dynamics of Trifolium species over 2 years at our field sites. Our results demonstrate that interactions with shared mutualists can stabilize the coexistence of closely related species.

To Fix or Not To Fix: Controls on Free-Living Nitrogen Fixation in the Rhizosphere.

Free-living nitrogen fixation (FLNF) in the rhizosphere, or N fixation by heterotrophic bacteria living on/near root surfaces, is ubiquitous and a significant source of N in some terrestrial systems. FLNF is also of interest in crop production as an alternative to chemical fertilizer, potentially reducing production costs and ameliorating negative environmental impacts of fertilizer N additions. Despite this interest, a mechanistic understanding of controls (e.g., carbon, oxygen, nitrogen, and nutrient availability) on FLNF in the rhizosphere is lacking but necessary. FLNF is distinct from and occurs under more diverse and dynamic conditions than symbiotic N fixation; therefore, predicting FLNF rates and understanding controls on FLNF has proven difficult. This has led to large gaps in our understanding of FLNF, and studies aimed at identifying controls on FLNF are needed. Here, we provide a mechanistic overview of FLNF, including how various controls may influence FLNF in the rhizosphere in comparison with symbiotic N fixation occurring in plant nodules where environmental conditions are moderated by the plant. We apply this knowledge to a real-world example, the bioenergy crop switchgrass (Panicum virgatum), to provide context of how FLNF may function in a managed system. We also highlight future challenges to assessing FLNF and understanding how FLNF functions in the environment and significantly contributes to plant N availability and productivity.

Unfair trade underground revealed by integrating data with Nash bargaining models.

Mutually beneficial resource exchange is fundamental to global biogeochemical cycles and plant and animal nutrition. However, there is inherent potential conflict in mutualisms, as each organism benefits more when the exchange ratio ('price') minimizes its own costs and maximizes its benefits. Understanding the bargaining power that each partner has in these interactions is key to our ability to predict the exchange ratio and therefore the functionality of the cell, organism, community and ecosystem. We tested whether partners have symmetrical ('fair') or asymmetrical ('unfair') bargaining power in a legume-rhizobia nitrogen-fixing symbiosis using measurements of carbon and nitrogen dynamics in a mathematical modeling framework derived from economic theory. A model of symmetric bargaining power was not consistent with our data. Instead, our data indicate that the growth benefit to the plant (Medicago truncatula) has greater weight in determining trade dynamics than the benefit to the bacteria. Quantitative estimates of the relative power of the plant revealed that the plant's influence rises as soil nitrogen availability decreases and trade benefits to both partners increase. Our finding that M. truncatula legumes have more bargaining power than their rhizobial partner at lower nitrogen availabilities highlights the importance of context-dependence for the evolution of mutualism with increasing nutrient deposition.

Rapid establishment of a flowering cline in Medicago polymorpha after invasion of North America.

To establish and spread in a new location, an invasive species must be able to carry out its life cycle in novel environmental conditions. A key trait underlying fitness is the shift from vegetative to reproductive growth through floral development. In this study, we used a common garden experiment and genotyping-by-sequencing to test whether the latitudinal flowering cline of the North American invasive plant Medicago polymorpha was translocated from its European native range through multiple introductions, or whether the cline rapidly established due to evolution following a genetic bottleneck. Analysis of flowering time in 736 common garden plants showed a latitudinal flowering time cline in both the native and invaded ranges where genotypes from lower latitudes flowered earlier. Genotyping-by-sequencing of 9,658 SNPs in 446 individuals revealed two major subpopulations of M. polymorpha in the native range, only one of which is present in the invaded range. Additionally, native range populations have higher genetic diversity than invaded range populations, suggesting that a genetic bottleneck occurred during invasion. All invaded range individuals are closely related to plants collected from native range populations in Portugal and southern Spain, and population assignment tests assigned invaded range individuals to this same narrow source region. Taken together, our results suggest that latitudinal clinal variation in flowering time has rapidly evolved across the invaded range despite a genetic bottleneck following introduction.

Modelling nutritional mutualisms: challenges and opportunities for data integration.

Nutritional mutualisms are ancient, widespread, and profoundly influential in biological communities and ecosystems. Although much is known about these interactions, comprehensive answers to fundamental questions, such as how resource availability and structured interactions influence mutualism persistence, are still lacking. Mathematical modelling of nutritional mutualisms has great potential to facilitate the search for comprehensive answers to these and other fundamental questions by connecting the physiological and genomic underpinnings of mutualisms with ecological and evolutionary processes. In particular, when integrated with empirical data, models enable understanding of underlying mechanisms and generalisation of principles beyond the particulars of a given system. Here, we demonstrate how mathematical models can be integrated with data to address questions of mutualism persistence at four biological scales: cell, individual, population, and community. We highlight select studies where data has been or could be integrated with models to either inform model structure or test model predictions. We also point out opportunities to increase model rigour through tighter integration with data, and describe areas in which data is urgently needed. We focus on plant-microbe systems, for which a wealth of empirical data is available, but the principles and approaches can be generally applied to any nutritional mutualism.

Rubrobacter spartanus sp. nov., a moderately thermophilic oligotrophic bacterium isolated from volcanic soil.

Bacterial strain HPK2-2T was isolated from soil adjacent to the caldera of Kilauea Volcano in Hawaii Volcanoes National Park. HPK2-2T is a chemoorganoheterotroph that shows optimal growth at 50 °C (range 45-55 °C) and pH 8.0 (range 5.0-10.0). Sequence analysis of the 16S subunit of the rRNA gene showed that HPK2-2T is most closely related to the type strain of Rubrobactertaiwanensis (ATCC BAA-406T), with which it shared 94.5 % sequence identity. The major fatty acids detected in HPK2-2T were C18 : 0 14-methyl and C16 : 0 12-methyl; internally branched fatty acids such as these are characteristic of the genus Rubrobacter. The only respiratory quinone detected was MK-8, which is the major respiratory quinone for all members of the family Rubrobacteraceae examined thus far. We propose that HPK2-2T represents a novel species of the genus Rubrobacter, for which we propose the name Rubrobacterspartanus (type strain HPK2-2T; DSM 102139T; LMG 29988T).

Genetic variation of transgenerational plasticity of offspring germination in response to salinity stress and the seed transcriptome of Medicago truncatula.

BACKGROUND: Transgenerational plasticity provides phenotypic variation that contributes to adaptation. For plants, the timing of seed germination is critical for offspring survival in stressful environments, as germination timing can alter the environmental conditions a seedling experiences. Stored seed transcripts are important determinants of seed germination, but have not previously been linked with transgenerational plasticity of germination behavior. In this study we used RNAseq and growth chamber experiments of the model legume M. trucantula to test whether parental exposure to salinity stress influences the expression of stored seed transcripts and early offspring traits and test for genetic variation. RESULTS: We detected genotype-dependent parental environmental effects (transgenerational plasticity) on the expression levels of stored seed transcripts, seed size, and germination behavior of four M. truncatula genotypes. More than 50% of the transcripts detected in the mature, ungerminated seed transcriptome were annotated as regulating seed germination, some of which are involved in abiotic stress response and post-embryonic development. Some genotypes showed increased seed size in response to parental exposure to salinity stress, but no parental environmental influence on germination timing. In contrast, other genotypes showed no seed size differences across contrasting parental conditions but displayed transgenerational plasticity for germimation timing, with significantly delayed germination in saline conditions when parental plants were exposed to salinity. In genotypes that show significant transgenerational plastic germination response, we found significant coexpression networks derived from salt responsive transcripts involved in post-transcriptional regulation of the germination pathway. Consistent with the delayed germination response to saline conditions in these genotypes, we found genes associated with dormancy and up-regulation of abscisic acid (ABA). CONCLUSIONS: Our results demonstrate genetic variation in transgenerational plasticity within M. truncatula and show that parental exposure to salinity stress influences the expression of stored seed transcripts, seed weight, and germination behavior. Furthermore, we show that the parental environment influences gene expression to modulate biological pathways that are likely responsible for offspring germination responses to salinity stress.

Cheaters must prosper: reconciling theoretical and empirical perspectives on cheating in mutualism.

Cheating is a focal concept in the study of mutualism, with the majority of researchers considering cheating to be both prevalent and highly damaging. However, current definitions of cheating do not reliably capture the evolutionary threat that has been a central motivation for the study of cheating. We describe the development of the cheating concept and distill a relative-fitness-based definition of cheating that encapsulates the evolutionary threat posed by cheating, i.e. that cheaters will spread and erode the benefits of mutualism. We then describe experiments required to conclude that cheating is occurring and to quantify fitness conflict more generally. Next, we discuss how our definition and methods can generate comparability and integration of theory and experiments, which are currently divided by their respective prioritisations of fitness consequences and traits. To evaluate the current empirical evidence for cheating, we review the literature on several of the best-studied mutualisms. We find that although there are numerous observations of low-quality partners, there is currently very little support from fitness data that any of these meet our criteria to be considered cheaters. Finally, we highlight future directions for research on conflict in mutualisms, including novel research avenues opened by a relative-fitness-based definition of cheating.

Genotype-phenotype mapping in a post-GWAS world.

Understanding how metabolic reactions, cell signaling, and developmental pathways translate the genome of an organism into its phenotype is a grand challenge in biology. Genome-wide association studies (GWAS) statistically connect genotypes to phenotypes, without any recourse to known molecular interactions, whereas a molecular biology approach directly ties gene function to phenotype through gene regulatory networks (GRNs). Using natural variation in allele-specific expression, GWAS and GRN approaches can be merged into a single framework via structural equation modeling (SEM). This approach leverages the myriad of polymorphisms in natural populations to elucidate and quantitate the molecular pathways that underlie phenotypic variation. The SEM framework can be used to quantitate a GRN, evaluate its consistency across environments or sexes, identify the differences in GRNs between species, and annotate GRNs de novo in non-model organisms.

The ecological genomic basis of salinity adaptation in Tunisian Medicago truncatula.

BACKGROUND: As our world becomes warmer, agriculture is increasingly impacted by rising soil salinity and understanding plant adaptation to salt stress can help enable effective crop breeding. Salt tolerance is a complex plant phenotype and we know little about the pathways utilized by naturally tolerant plants. Legumes are important species in agricultural and natural ecosystems, since they engage in symbiotic nitrogen-fixation, but are especially vulnerable to salinity stress. RESULTS: Our studies of the model legume Medicago truncatula in field and greenhouse settings demonstrate that Tunisian populations are locally adapted to saline soils at the metapopulation level and that saline origin genotypes are less impacted by salt than non-saline origin genotypes; these populations thus likely contain adaptively diverged alleles. Whole genome resequencing of 39 wild accessions reveals ongoing migration and candidate genomic regions that assort non-randomly with soil salinity. Consistent with natural selection acting at these sites, saline alleles are typically rare in the range-wide species' gene pool and are also typically derived relative to the sister species M. littoralis. Candidate regions for adaptation contain genes that regulate physiological acclimation to salt stress, such as abscisic acid and jasmonic acid signaling, including a novel salt-tolerance candidate orthologous to the uncharacterized gene AtCIPK21. Unexpectedly, these regions also contain biotic stress genes and flowering time pathway genes. We show that flowering time is differentiated between saline and non-saline populations and may allow salt stress escape. CONCLUSIONS: This work nominates multiple potential pathways of adaptation to naturally stressful environments in a model legume. These candidates point to the importance of both tolerance and avoidance in natural legume populations. We have uncovered several promising targets that could be used to breed for enhanced salt tolerance in crop legumes to enhance food security in an era of increasing soil salinization.

Back to the future: Using herbarium specimens to isolate nodule-associated bacteria.

Herbarium specimens are increasingly being used as sources of information to understand the ecology and evolution of plants and their associated microbes. Most studies have used specimens as a source of genetic material using culture-independent approaches. We demonstrate that herbarium specimens can also be used to culture nodule-associated bacteria, opening the possibility of using specimens to understand plant-microbe interactions at new spatiotemporal scales. We used historic and contemporary nodules of a common legume, Medicago lupulina, to create a culture collection. We were able to recover historic bacteria in 15 genera from three specimens (collected in 1950, 2004, and 2015). This work is the first of its kind to isolate historic bacteria from herbarium specimens. Future work should include inoculating plants with historic strains to see if they produce nodules and if they affect plant phenotype and fitness. Although we were unable to recover any Ensifer, the main symbiont of Medicago lupulina, we recovered some other potential nodulating species, as well as many putative growth-promoting bacteria.