Mapping the genetic basis of symbiotic variation in legume-rhizobium interactions in Medicago truncatula

Genetic diversity, population structure, and taxonomic confirmation in annual medic (Medicago spp.) collections from Crimea, Ukraine

Annual medic (Medicago spp.) germplasm was collected from the Crimean Peninsula of Ukraine in 2008 to fill gaps in geographic coverage in the United States department of Agriculture, Agricultural Research Service, National Plant Germplasm System (NPGS) temperate-adapted forage legume collection. A total of 102 accessions across 10 Medicago species were collected. To assess genetic diversity, population structure, and to confirm taxonomic identities, the collections were phenotypically and genetically characterized. Phenotyping included the use of 24 descriptor traits while genetic characterization was accomplished using a 3K Diversity Array Technologies (DArTag) panel developed for alfalfa (Medicago sativa L.). For both field and molecular characterizations, a reference set of 92 geographically diverse and species-representative accessions were obtained from the NPGS collection. Phenotypic descriptors showed consistency among replicated plants within accessions, some variation across accessions within species, and evident distinctions between species. Because the DArTag panel was developed for cultivated alfalfa, the transferability of markers to the species being evaluated was limited, resulting in an average of ~1,500 marker loci detected per species. From these loci, 448 markers were present in 95% of the samples. Principal component and phylogenetic analysis based on a larger set of 2,396 selected markers clustered accessions by species and predicted evolutionary relationships among species. Additionally, the markers aided in the taxonomic identity of a few accessions that were likely mislabeled. The genotyping results also showed that sampling individual plants for these mostly self-pollinating species is sufficient due to high reproducibility between single (n=3) and pooled (n=7) biological replicate leaf samples. The phenotyping and the 2,396 Single Nucleotide Polymorphism (SNP) marker set were useful in estimating population structure in the Crimean and reference accessions, highlighting novel and unique genetic diversity captured in the Crimean accessions. This research not only demonstrated the utility of the DArTag marker panel in evaluating the Crimean germplasm but also highlighted its broader application in assessing genetic resources within the Medicago genus. Furthermore, we anticipate that our findings will underscore the importance of leveraging genetic resources and advanced genotyping tools for sustainable crop improvement and biodiversity conservation in annual medic species.

Wolbachia strains wMel and wAlbB differentially affect Aedes aegypti traits related to fecundity

Two Wolbachia strains, wMel and wAlbB, have been transinfected into Aedes aegypti mosquitoes for population replacement with the aim of reducing dengue transmission. Epidemiological data from various endemic sites suggest a pronounced decrease in dengue transmission after implementing this strategy. In this study, we investigated the impact of the Wolbachia strains wMel and wAlbB on Ae. aegypti fitness in a common genetic background. We found that Ae. aegypti females infected with the wMel strain exhibited several significant differences compared with those infected with the wAlbB strain. Specifically, wMel-infected females laid significantly fewer eggs, ingested a lower amount of blood, had a reduced egg production rate, and exhibited a decreased Wolbachia density at a later age compared with mosquitoes infected with the wAlbB strain. Conversely, the wAlbB strain showed only mild negative effects when compared with Wolbachia-uninfected specimens. These differential effects on Ae. aegypti fitness following infection with either wMel or wAlbB may have important implications for the success of population replacement strategies in invading native Ae. aegypti populations in endemic settings. Further research is needed to better understand the underlying mechanisms responsible for these differences in fitness effects and their potential impact on the long-term efficacy of Wolbachia-based dengue control programs.IMPORTANCEThe transmission of arboviruses such as dengue, Zika, and chikungunya is on the rise globally. Among the most promising strategies to reduce arbovirus burden is the release of one out of two strains of Wolbachia-infected Aedes aegypti: wMel and wAlbB. One critical aspect of whether this approach will succeed involves the fitness cost of either Wolbachia strains on mosquito life history traits. For instance, we found that wMel-infected Ae. aegypti females laid significantly fewer eggs, ingested a lower amount of blood, had a reduced egg production rate, and exhibited a decreased Wolbachia density at a later age compared with mosquitoes infected with the wAlbB strain. Conversely, the wAlbB strain showed only mild negative effects when compared with Wolbachia-uninfected specimens. These differential effects on mosquito fitness following infection with either wMel or wAlbB may have important implications for the success of population replacement strategies in invading native Ae. aegypti populations.

Genome-Wide Association Studies across Environmental and Genetic Contexts Reveal Complex Genetic Architecture of Symbiotic Extended Phenotypes

A goal of modern biology is to develop the genotype-phenotype (G→P) map, a predictive understanding of how genomic information generates trait variation that forms the basis of both natural and managed communities. As microbiome research advances, however, it has become clear that many of these traits are symbiotic extended phenotypes, being governed by genetic variation encoded not only by the host's own genome, but also by the genomes of myriad cryptic symbionts. Building a reliable G→P map therefore requires accounting for the multitude of interacting genes and even genomes involved in symbiosis. Here, we use naturally occurring genetic variation in 191 strains of the model microbial symbiont Sinorhizobium meliloti paired with two genotypes of the host Medicago truncatula in four genome-wide association studies (GWAS) to determine the genomic architecture of a key symbiotic extended phenotype-partner quality, or the fitness benefit conferred to a host by a particular symbiont genotype, within and across environmental contexts and host genotypes. We define three novel categories of loci in rhizobium genomes that must be accounted for if we want to build a reliable G→P map of partner quality; namely, (i) loci whose identities depend on the environment, (ii) those that depend on the host genotype with which rhizobia interact, and (iii) universal loci that are likely important in all or most environments. IMPORTANCE Given the rapid rise of research on how microbiomes can be harnessed to improve host health, understanding the contribution of microbial genetic variation to host phenotypic variation is pressing, and will better enable us to predict the evolution of (and select more precisely for) symbiotic extended phenotypes that impact host health. We uncover extensive context-dependency in both the identity and functions of symbiont loci that control host growth, which makes predicting the genes and pathways important for determining symbiotic outcomes under different conditions more challenging. Despite this context-dependency, we also resolve a core set of universal loci that are likely important in all or most environments, and thus, serve as excellent targets both for genetic engineering and future coevolutionary studies of symbiosis.

Effects of Domestication on Plant-Microbiome Interactions

Through the process of domestication, selection is targeted on a limited number of plant traits that are typically associated with yield. As an unintended consequence, domesticated plants often perform poorly compared to their wild progenitors for a multitude of traits that were not under selection during domestication, including abiotic and biotic stress tolerance. Over the past decade, advances in sequencing technology have allowed for the rigorous characterization of host-associated microbial communities, termed the microbiome. It is now clear that nearly every conceivable plant interaction with the environment is mediated by interactions with the microbiome. For this reason, plant-microbiome interactions are an area of great promise for plant breeding and crop improvement. Here, we review the literature to assess the potential impact that domestication has had on plant-microbiome interactions and the current understanding of the genetic basis of microbiome variation to inform plant breeding efforts. Overall, we find limited evidence that domestication impacts the diversity of microbiomes, but domestication is often associated with shifts in the abundance and composition of microbial communities, including taxa of known functional significance. Moreover, genome-wide association studies and mutant analysis have not revealed a consistent set of core candidate genes or genetic pathways that confer variation in microbiomes across systems. However, such studies do implicate a consistent role for plant immunity, root traits, root and leaf exudates and cell wall integrity as key traits that control microbiome colonization and assembly. Therefore, selection on these key traits may pose the most immediate promise for enhancing plant-microbiome interactions through breeding.

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.

The Potential for Genotype-by-Environment Interactions to Maintain Genetic Variation in a Model Legume-Rhizobia Mutualism

The maintenance of genetic variation in mutualism-related traits is key for understanding mutualism evolution, yet the mechanisms maintaining variation remain unclear. We asked whether genotype-by-environment (G×E) interaction is a potential mechanism maintaining variation in the model legume-rhizobia system, Medicago truncatula-Ensifer meliloti. We planted 50 legume genotypes in a greenhouse under ambient light and shade to reflect reduced carbon availability for plants. We found an expected reduction under shaded conditions for plant performance traits, such as leaf number, aboveground and belowground biomass, and a mutualism-related trait, nodule number. We also found G×E for nodule number, with ∼83% of this interaction due to shifts in genotype fitness rank order across light environments, coupled with strong positive directional selection on nodule number regardless of light environment. Our results suggest that G×E can maintain genetic variation in a mutualism-related trait that is under consistent positive directional selection across light environments.

Networks Underpinning Symbiosis Revealed Through Cross-Species eQTL Mapping

Organisms engage in extensive cross-species molecular dialog, yet the underlying molecular actors are known for only a few interactions. Many techniques have been designed to uncover genes involved in signaling between organisms. Typically, these focus on only one of the partners. We developed an expression quantitative trait locus (eQTL) mapping-based approach to identify cause-and-effect relationships between genes from two partners engaged in an interspecific interaction. We demonstrated the approach by assaying expression of 98 isogenic plants (Medicago truncatula), each inoculated with a genetically distinct line of the diploid parasitic nematode Meloidogyne hapla With this design, systematic differences in gene expression across host plants could be mapped to genetic polymorphisms of their infecting parasites. The effects of parasite genotypes on plant gene expression were often substantial, with up to 90-fold (P = 3.2 × 10-52) changes in expression levels caused by individual parasite loci. Mapped loci included a number of pleiotropic sites, including one 87-kb parasite locus that modulated expression of >60 host genes. The 213 host genes identified were substantially enriched for transcription factors. We distilled higher-order connections between polymorphisms and genes from both species via network inference. To replicate our results and test whether effects were conserved across a broader host range, we performed a confirmatory experiment using M. hapla-infected tomato. This revealed that homologous genes were similarly affected. Finally, to validate the broader utility of cross-species eQTL mapping, we applied the strategy to data from a Salmonella infection study, successfully identifying polymorphisms in the human genome affecting bacterial expression.

No evidence for adaptation to local rhizobial mutualists in the legume Medicago lupulina

Local adaptation is a common but not ubiquitous feature of species interactions, and understanding the circumstances under which it evolves illuminates the factors that influence adaptive population divergence. Antagonistic species interactions dominate the local adaptation literature relative to mutualistic ones, preventing an overall assessment of adaptation within interspecific interactions. Here, we tested whether the legume Medicago lupulina is adapted to the locally abundant species of mutualistic nitrogen-fixing rhizobial bacteria that vary in frequency across its eastern North American range. We reciprocally inoculated northern and southern M. lupulina genotypes with the northern (Ensifer medicae) or southern bacterium (E. meliloti) in a greenhouse experiment. Despite producing different numbers of root nodules (the structures in which the plants house the bacteria), neither northern nor southern plants produced more seeds, flowered earlier, or were more likely to flower when inoculated with their local rhizobia. We then used a pre-existing dataset to perform a genome scan for loci that showed elevated differentiation between field-collected plants that hosted different bacteria. None of the loci we identified belonged to the well-characterized suite of legume-rhizobia symbiosis genes, suggesting that the rhizobia do not drive genetic divergence between M. lupulina populations. Our results demonstrate that symbiont local adaptation has not evolved in this mutualism despite large-scale geographic variation in the identity of the interacting species.

Testing GxG interactions between coinfecting microbial parasite genotypes within hosts

Host-parasite interactions represent one of the strongest selection pressures in nature. They are often governed by genotype-specific (GxG) interactions resulting in host genotypes that differ in resistance and parasite genotypes that differ in virulence depending on the antagonist's genotype. Another type of GxG interactions, which is often neglected but which certainly influences host-parasite interactions, are those between coinfecting parasite genotypes. Mechanistically, within-host parasite interactions may range from competition for limited host resources to cooperation for more efficient host exploitation. The exact type of interaction, i.e., whether competitive or cooperative, is known to affect life-history traits such as virulence. However, the latter has been shown for chosen genotype combinations only, not considering whether the specific genotype combination per se may influence the interaction (i.e., GxG interactions). Here, we want to test for the presence of GxG interactions between coinfections of the bacterium Bacillus thuringiensis infecting the nematode Caenorhabditis elegans by combining two non-pathogenic and five pathogenic strains in all possible ways. Furthermore, we evaluate whether the type of interaction, reflected by the direction of virulence change of multiple compared to single infections, is genotype-specific. Generally, we found no indication for GxG interactions between non-pathogenic and pathogenic bacterial strains, indicating that virulence of pathogenic strains is equally affected by both non-pathogenic strains. Specific genotype combinations, however, differ in the strength of virulence change, indicating that the interaction type between coinfecting parasite strains and thus the virulence mechanism is specific for different genotype combinations. Such interactions are expected to influence host-parasite interactions and to have strong implications for coevolution.

Connecting functional and statistical definitions of genotype by genotype interactions in coevolutionary studies

Predicting how species interactions evolve requires that we understand the mechanistic basis of coevolution, and thus the functional genotype-by-genotype interactions (G × G) that drive reciprocal natural selection. Theory on host-parasite coevolution provides testable hypotheses for empiricists, but depends upon models of functional G × G that remain loosely tethered to the molecular details of any particular system. In practice, reciprocal cross-infection studies are often used to partition the variation in infection or fitness in a population that is attributable to G × G (statistical G × G). Here we use simulations to demonstrate that within-population statistical G × G likely tells us little about the existence of coevolution, its strength, or the genetic basis of functional G × G. Combined with studies of multiple populations or points in time, mapping and molecular techniques can bridge the gap between natural variation and mechanistic models of coevolution, while model-based statistics can formally confront coevolutionary models with cross-infection data. Together these approaches provide a robust framework for inferring the infection genetics underlying statistical G × G, helping unravel the genetic basis of coevolution.

Genetic variability and QTL mapping of freezing tolerance and related traits in Medicago truncatula

Freezing is a major environmental limitation to crop productivity for a number of species including legumes. We investigated the genetic determinism of freezing tolerance in the model legume Medicago truncatula Gaertn (M. truncatula). After having observed a large variation for freezing tolerance among 15 M. truncatula accessions, the progeny of a F6 recombinant inbred line population, derived from a cross between two accessions, was acclimated to low above-freezing temperatures and assessed for: (a) number of leaves (NOL), leaf area (LA), chlorophyll content index (CCI), shoot and root dry weights (SDW and RDW) at the end of the acclimation period and (b) visual freezing damage (FD) during the freezing treatment and 2 weeks after regrowth and foliar electrolyte leakage (EL) 2 weeks after regrowth. Consistent QTL positions with additive effects for FD were found on LG1, LG4 and LG6, the latter being the most explanatory (R (2) ≈ 40 %). QTL for NOL, QTL for EL, NOL and RDW, and QTL for EL and CCI colocalized with FD QTL on LG1, LG4 and LG6, respectively. Favorable alleles for these additive effects were brought by the same parent suggesting that this accession contributes to superior freezing tolerance by affecting plants' capacity to maintain growth at low above-freezing temperatures. No epistatic effects were found between FD QTL, but for each of the studied traits, 3-6 epistatic effects were detected between loci not detected directly as QTL. These results open the way to the assessment of syntenic relationships between QTL for frost tolerance in M. truncatula and cultivated legume species.

The genetic basis of speciation in the Giliopsis lineage of Ipomopsis (Polemoniaceae)

One of the most powerful drivers of speciation in plants is pollinator-mediated disruptive selection, which leads to the divergence of floral traits adapted to the morphology and behavior of different pollinators. Despite the widespread importance of this speciation mechanism, its genetic basis has been explored in only a few groups. Here, we characterize the genetic basis of pollinator-mediated divergence of two species in genus Ipomopsis, I. guttata and I. tenuifolia, using quantitative trait locus (QTL) analyses of floral traits and other variable phenotypes. We detected one to six QTLs per trait, with each QTL generally explaining small to modest amounts of the phenotypic variance of a backcross hybrid population. In contrast, flowering time and anthocyanin abundance (a metric of color variation) were controlled by a few QTLs of relatively large effect. QTLs were strongly clustered within linkage groups, with 26 of 37 QTLs localized to six marker-interval 'hotspots,' all of which harbored pleiotropic QTLs. In contrast to other studies that have examined the genetic basis of pollinator shifts, our results indicate that, in general, mutations of small to modest effect on phenotype were involved. Thus, the evolutionary transition between the distinct pollination modes of I. guttata and I. tenuifolia likely proceeded incrementally, rather than saltationally.