Search

Rhizobium strains that nodulate Egyptian winter clover belong to two symbiotic species, Rhizobium aegyptiacum sv. trifolii and Rhizobium bangladeshense sv. trifolii. Rhizobia from Egyptian clover are often tolerant of extreme environmental conditions and are generally poorly characterized, despite the importance of this crop species. Here, we report on the comparative genomic analysis of four stress-tolerant, highly effective nitrogen-fixing, Egyptian-clover-nodulating rhizobia, R. aegyptiacum sv. trifolii strain Rhiz950 and Rhizobium bangladeshense sv. trifolii strains Rhiz1002, Rhiz1017 and Rhiz1024, as well as other diverse Rhizobium strains. While Rhizobium bangladeshense sv. trifolii strains Rhiz1002, Rhiz1017, and Rhiz1024 had similar genome sizes of around 6.6 Mb, the salt tolerant strain Rhiz950 had a slightly larger genome size of 7.4 Mb. All four strains lacked genes responsible for methylation and sulfation of nod factors (nodS and nodHPQ). They also lacked nodZ, and instead had the functionally-equivalent nod factor fucosylation gene nolK. No strains contained homocitrate synthase (encoded by nifV) that is essential for free-living nitrogen fixing by diazotrophs. A number of genes unique to the Egyptian clover strains and to other rhizobia are potential candidates for variation in stress tolerance and host-range.

Leguminous plants form symbiotic relationships with rhizobia. These nitrogen-fixing bacteria live in specialized root organs called nodules. Although rhizobia form the most notable host relationship within root nodules, other bacterial endophytes also inhabit nodules and can influence host–rhizobia interactions as well as exert effects of their own. In this study, we elucidate the effect of intraspecific host genetic variation and between generational legacies on root nodule endophyte communities in Medicago truncatula. Although the diversity of endophytes in nodules was similar across hosts, both nodule endophyte composition and gene functional groups differed. In contrast, neither the presence nor identity of a host in the previous generation (either A17 or R108) had a significant effect on the nodule endophyte diversity or composition. However, whether or not a host was present significantly altered gene functional groups. We conclude that genetic variation within a legume host species can play both direct and indirect effects on establishment of nodule microbiomes. Further studies, including genome-wide association studies and functional assays, can open the door for engineering and optimizing nodule endophyte communities that promote growth or have other beneficial qualities.

Changing climate and land‐use practices have the potential to bring previously isolated populations of pest insects into new sympatry. This heightens the need to better understand how differing patterns of host–plant association, and unique endosymbionts, serve to promote genetic isolation or integration. We addressed these factors in populations of potato psyllid, Bactericera cockerelli (Šulc), a generalist herbivore that vectors a bacterial pathogen (Candidatus Liberibacter solanacearum, causal pathogen of zebra chip disease) of potato (Solanum tuberosum L.). Genome‐wide SNP data revealed two major genetic clusters—psyllids collected from potato crops were genetically similar to psyllids found on a common weed, Lycium spp., but dissimilar from those found on another common non‐crop host, Solanum dulcamara L. Most psyllids found on Lycium spp. and potato represented a single mitochondrial cytochrome oxidase I (COI) haplotype that has been suggested to not be native to the region, and whose arrival may have been concurrent with zebra chip disease first emerging. The putatively introduced COI haplotype usually co‐occurred with endosymbiotic Wolbachia, while the putatively resident COI haplotype generally did not. Genetic intermediates between the two genetic populations of insects were rare, consistent with recent sympatry or reproductive isolation, although admixture patterns of apparent hybrids were consistent with introgression of genes from introduced into resident populations. Our results suggest that both host–plant associations and endosymbionts are shaping the population genetic structure of sympatric psyllid populations associated with different non‐crop hosts. It is of future interest to explicitly examine vectorial capacity of the two populations and their potential hybrids, as population structure and hybridization might alter regional vector capacity and disease outbreaks.

Coevolution is predicted to depend on how the genetic diversity of interacting species is geographically structured. Plant-microbe symbioses such as the legume-rhizobium mutualism are ecologically and economically important, but distinct life history and dispersal mechanisms for these host and microbial partners, plus dynamic genome composition in bacteria, present challenges for understanding spatial genetic processes in these systems. Here we study the model rhizobium Ensifer meliloti using a hierarchically-structured sample of 191 strains from 21 sites in the native range and compare its population structure to that of its host plant Medicago truncatula. We find high local genomic variation and minimal isolation by distance across the rhizobium genome, particularly at the two symbiosis elements pSymA and pSymB, which have evolutionary histories and population structures that are similar to each other but distinct from both the chromosome and the host. While the chromosome displays weak isolation by distance, it is uncorrelated with hosts. Patterns of discordant population structure among elements with the bacterial genome has implications for bacterial adaptation to life in the soil versus symbiosis, while discordant population genetic structure of hosts and microbes might restrict local adaptation of species to each other and give rise to phenotypic mismatches in coevolutionary traits.

Genetic studies of legume symbiosis with nitrogen-fixing rhizobial bacteria have traditionally focused on nodule and nitrogen-fixation phenotypes when hosts are inoculated with a single rhizobial strain. These approaches overlook the potential effect of host genes on rhizobial fitness (i.e. how many rhizobia are released from host nodules) and strain-specific effects of host genes (i.e. genome × genome interactions). Using Medicago truncatula mutants in the recently described nodule-specific PLAT domain (NPD) gene family, we show how inoculating plants with a mixed inoculum of 68 rhizobial strains (Ensifer meliloti) via a select-and-resequence approach can be used to efficiently assay host mutants for strain-specific effects of late-acting host genes on interacting bacteria. The deletion of a single NPD gene (npd2) or all five members of the NPD gene family (npd1-5) differentially altered the frequency of rhizobial strains in nodules even though npd2 mutants had no visible nodule morphology or N-fixation phenotype. Also, npd1-5 nodules were less diverse and had larger populations of colony-forming rhizobia despite their smaller size. Lastly, NPD mutations disrupt a positive correlation between strain fitness and wild-type host biomass. These changes indicate that the effects of NPD proteins are strain dependent and that NPD family members are not redundant with regard to their effects on rhizobial strains. Association analyses of the rhizobial strains in the mixed inoculation indicate that rhizobial genes involved in chromosome segregation, cell division, GABA metabolism, efflux systems, and stress tolerance play an important role in the strain-specific effects of NPD genes.

Understanding the genomic basis of adaptation is critical for understanding evolutionary processes and predicting how species will respond to environmental change. Spinner dolphins in the eastern tropical Pacific (ETP) present a unique system for studying adaptation. Within this large geographical region are four spinner dolphin ecotypes with weak neutral genetic divergence and no obvious barriers to gene flow, but strong spatial variation in morphology, behaviour and habitat. These ecotypes have large population sizes, which could reduce the effects of drift and facilitate selection. To identify genomic regions putatively under divergent selective pressures between ecotypes, we used genome scans with 8994 RADseq single nucleotide polymorphisms (SNPs) to identify population differentiation outliers and genotype–environment association outliers. Gene ontology enrichment analyses indicated that outlier SNPs from both types of analyses were associated with multiple genes involved in social behaviour and hippocampus development, including 15 genes associated with the human social disorder autism. Evidence for divergent selection on social behaviour is supported by previous evidence that these spinner dolphin ecotypes differ in mating systems and associated social behaviours. In particular, three of the ETP ecotypes probably have a polygynous mating system characterized by strong premating competition among males, whereas the fourth ecotype probably has a polygynandrous mating system characterized by strong postmating competition such as sperm competition. Our results provide evidence that selection for social behaviour may be an evolutionary force driving diversification of spinner dolphins in the ETP, potentially as a result of divergent sexual selection associated with different mating systems. Future studies should further investigate the potential adaptive role of the candidate genes identified here, and could probably find further signatures of selection using whole genome sequence data.

Horizontal transfer (HT) alters the repertoire of symbiosis genes in rhizobial genomes and may play an important role in the on-going evolution of the rhizobia–legume symbiosis. To gain insight into the extent of HT of symbiosis genes with different functional roles (nodulation, N-fixation, host benefit and rhizobial fitness), we conducted comparative genomic and selection analyses of the full-genome sequences from 27 rhizobial genomes. We find that symbiosis genes experience high rates of HT among rhizobial lineages but also bear signatures of purifying selection (low Ka : Ks). HT and purifying selection appear to be particularly strong in genes involved in initiating the symbiosis (e.g. nodulation) and in genome-wide association candidates for mediating benefits provided to the host. These patterns are consistent with rhizobia adapting to the host environment through the loss and gain of symbiosis genes, but not with host-imposed positive selection driving divergence of symbiosis genes through recurring bouts of positive selection.

The environmental context of the nitrogen-fixing mutualism between leguminous plants and rhizobial bacteria varies over space and time. Variation in resource availability, population density, and composition likely affect the ecology and evolution of rhizobia and their symbiotic interactions with hosts. We examined how host genotype, nitrogen addition, rhizobial density, and community complexity affected selection on 68 rhizobia strains in the Ensifer meliloti - Medicago truncatula mutualism. As expected, the host genotype had the most substantial effect on the size, number, and strain composition of root nodules (the symbiotic organ). The understudied environmental variable of rhizobial density had a more significant effect on strain frequency in nodules than the addition of low nitrogen levels. Higher inoculum density resulted in a nodule community that was less diverse and more beneficial but only in the context of the more selective host genotype. Higher density resulted in more diverse and less beneficial nodule communities with the less selective host. Density effects on strain composition deserve additional scrutiny as they can create eco-evolutionary feedback. Lastly, we found that relative strain rankings were stable across increasing community complexity (community complexity (2, 3, 8, or 68 strains). This unexpected result suggests that higher-order interactions between strains are rare in the context of host nodule formation and development. Taken together, our empirical work highlights the importance of developing new theoretical predictions that incorporate density dependence. Further, it has translational relevance for overcoming establishment barriers in bio-inoculants and motivating host breeding programs that maintain beneficial plant-microbe interactions across diverse agro-ecological contexts.IMPORTANCELegume cash, forage, and cover crops establish beneficial associations with rhizobial bacteria who perform biological nitrogen fixation (BNF)—providing Nitrogen (N) fertilizer to plants without the economic and greenhouse gas emission costs of chemical N inputs. Here, for the first time, we examine the relative influence of three environmental factors that vary in agricultural fields on strain relative fitness in nodules when scores rhizobial strains compete. In addition to manipulating Nitrogen, we also use two biotic variables that have rarely been examined: the rhizobial community’s density and complexity. Taken together, our results suggest 1) breeding legume varieties that select beneficial strains despite environmental variation are possible, 2) changes in rhizobial population densities that occur routinely in agricultural fields could drive evolutionary changes in rhizobia populations, and 3) the lack of higher-order interactions between strains will allow the high-throughput assessments of rhizobia winners and losers during plant interactions.

Tasmanian devils (Sarcophilus harrisii) are evolving in response to a unique transmissible cancer, devil facial tumour disease (DFTD), first described in 1996. Persistence of wild populations and the recent emergence of a second independently evolved transmissible cancer suggest that transmissible cancers may be a recurrent feature in devils. Here we compared signatures of selection across temporal scales to determine whether genes or gene pathways under contemporary selection (6-8 generations) have also been subject to historical selection (65-85 million years), and test for recurrent selection in devils. First, we used a targeted sequencing approach, RAD-capture, to identify genomic regions subject to rapid evolution in approximately 2,500 devils in six populations as DFTD spread across the species range. We documented genome-wide contemporary evolution, including 186 candidate genes related to cell cycling and immune response. Then we used a molecular evolution approach to identify historical positive selection in devils compared to other marsupials and found evidence of selection in 1,773 genes. However, we found limited overlap across time scales, with historical selection detected in only 16 contemporary candidate genes, and no overlap in enriched functional gene sets. Our results are consistent with a novel, multi-locus evolutionary response of devils to DFTD. Our results can inform management actions to conserve adaptive potential of devils by identifying high priority targets for genetic monitoring and maintenance of functional diversity in managed populations.

Natural enemies often move among habitats to track prey and resources. Indeed, biocontrol often depends on natural enemies dispersing into crops after disturbances such as tillage and pesticide applications. However, the small size of many natural enemies makes it difficult to observe such movements. Here we used genetic relatedness among entomopathogenic nematodes, tiny, soil-dwelling, and thus cryptic natural enemies, to infer their movement across an agricultural landscape. We collected strains of two nematode species, Heterorhabditis bacteriophora and Steinernema feltiae, by placing sentinel hosts into eight irrigated Solanum tuberosum fields across arid central Washington State. We then used restriction site associated DNA sequencing to generate single nucleotide polymorphisms and infer relatedness among strains across our study sites. We identified 3,367 and 138,286 polymorphic loci for H. bacteriophora and S. feltiae, respectively. Genetic differentiation for both species increased with greater distance between sites, although there was considerable variation. While strains collected from the same field were generally more closely related than those from different sites, for both species, strains from different fields were sometimes quite closely related. Altogether, our results suggest a surprising amount of genetic similarity among nematodes from distant sites, despite the presumably limited distances that can be traversed by individual worms. This is consistent with the nematodes moving, at least occasionally, between far-apart locations. More generally, we suggest that recent advances in population genomics are providing powerful new tools for mapping natural enemy movement across broad landscapes.

Emerging infectious diseases are increasingly recognized in species' declines and extinctions. Landscape genetics can be used as a tool to predict disease emergence and spread. The Tasmanian devil is threatened with extinction by a nearly 100% fatal transmissible cancer, which has spread across 95% of the species' geographic range in 20 years. Here, we present a landscape genetic analysis in the last remaining uninfected parts of the Tasmanian devil's geographic range to: describe population genetic structure, characterize genetic diversity, and test the influence of landscape variables on Tasmanian devil gene flow to assess the potential for disease spread. In contrast to previous genetic studies on Tasmanian devils showing evidence for two genetic populations island-wide, our genetic based assignment tests and spatial principal components analyses suggest at least two, and possibly three, populations in a study area that is approximately 15% of the size of the overall species' geographic range. Positive spatial autocorrelation declined at about 40 km, in contrast to 80 km in eastern populations, highlighting the need for range-wide genetic studies. Strong genetic structure was found between devils in the northern part of the study area and those found south of Macquarie Harbour, with weaker structure found between the northeastern and northwestern portion of our study area. Consistent with previous work, we found low overall genetic diversity, likely owing to a combination of founder effects and extreme weather events thousands of years ago that likely caused large-scale population declines. We also found possible signs of recent bottlenecks, perhaps resulting from forest clearing for dairy farming in the central part of the study area. This human disturbance also may have contributed to weak genetic structuring detected between the northeastern and northwestern part of the study area. Individual-based least cost path modeling showed limited influence of landscape variables on gene flow, with weak effects of variation in elevation in the northeast. In the northwest, however, landscape genetic models did not perform better than the null isolation-by-distance model. At the larger spatial scale of the northern part of the study area, elevation and temperatures were negatively correlated with gene flow, consistent with low dispersal suitability of higher elevation habitats that have lower temperatures and dense, wet vegetation. Overall, Tasmanian devils are a highly vagile species for which dispersal and gene flow appear to be influenced little by landscape features, and spread of devil facial tumor disease to the remaining portion of the devil's geographic range seems imminent. Nonetheless, strong genetic structure found between the northern and southern portions of our study area, combined with low densities and limited possible colonization of DFTD from the east suggest there is some time for implementation of management strategies.

Tasmanian devils face a combination of threats to persistence, including Devil Facial Tumor Disease (DFTD), an epidemic transmissible cancer. We used RAD sequencing to investigate genome-wide patterns of genetic diversity and geographic population structure. Consistent with previous results, we found very low genetic diversity in the species as a whole, and we detected two broad genetic clusters occupying the northwestern portion of the range, and the central and eastern portions. However, these two groups overlap across a broad geographic area, and differentiation between them is modest (FST = 0.1081). Our results refine the geographic extent of the zone of mixed ancestry and substructure within it, potentially informing management of genetic variation that existed in pre-diseased populations of the species. DFTD has spread across both genetic clusters, but recent evidence points to a genomic response to selection imposed by DFTD. Any allelic variation for resistance to DFTD may be able to spread across the devil population under selection by DFTD, and/or be present as standing variation in both genetic regions.

Landscape genomics studies focus on identifying candidate genes under selection via spatial variation in abiotic environmental variables, but rarely by biotic factors such as disease. The Tasmanian devil (Sarcophilus harrisii) is found only on the environmentally heterogeneous island of Tasmania and is threatened with extinction by a nearly 100% fatal, transmissible cancer, devil facial tumor disease (DFTD). Devils persist in regions of long-term infection despite epidemiological model predictions of species’ extinction, suggesting possible adaptation to DFTD. Here, we test the extent to which spatial variation and genetic diversity are associated with the abiotic environment and/or by DFTD. We employ genetic-environment association (GEAs) analyses using a RAD-capture panel consisting of 6,886 SNPs from 3,286 individuals sampled pre- and post-disease arrival. Pre-disease, we find significant correlations of allele frequencies with environmental variables, including 349 unique loci linked to 64 genes, suggesting local adaptation to abiotic environment. Following DFTD arrival, however, we detected few of the pre-disease candidate loci, but instead frequencies of candidate loci linked to 14 genes correlated with disease prevalence. Loss of apparent signal of abiotic local adaptation following disease arrival suggests swamping by the strong selection imposed by DFTD. Further support for this result comes from the fact that post-disease candidate loci are in linkage disequilibrium with genes putatively involved in immune response, tumor suppression and apoptosis. This suggests the strength GEA associations of loci with the abiotic environment are swamped resulting from the rapid onset of a biotic selective pressure.

Whether hybridization generates or erodes species diversity has long been debated, but to date most studies have been conducted at small taxonomic scales. Salamanders (order Caudata) represent a taxonomic order in which hybridization plays a prevalent ecological and evolutionary role. We employed a recently developed model of trait-dependent diversification to test the hypothesis that hybridization impacts the diversification dynamics of species that are currently hybridizing. We find strong evidence supporting this hypothesis, showing that hybridizing salamander lineages have significantly greater net-diversification rates than non-hybridizing lineages. This pattern is driven by concurrently increased speciation rates and decreased extinction rates in hybridizing lineages. Our results support the hypothesis that hybridization can act as a generative force in macroevolutionary diversification.

Understanding factors that cause species' geographic range limits is a major focus in ecology and evolution. The central marginal hypothesis (CMH) predicts that species cannot adapt to conditions beyond current geographic range edges because genetic diversity decreases from core to edge due to smaller, more isolated edge populations. We employed a population genomics framework using 24 235-33 112 SNP loci to test major predictions of the CMH in the ongoing invasion of the cane toad (Rhinella marina) in Australia. Cane toad tissue samples were collected along broad-scale, core-to-edge transects across their invasive range. Geographic and ecological core areas were identified using GIS and habitat suitability indices from ecological niche modelling. Bayesian clustering analyses revealed three genetic clusters, in the northwest invasion-front region, northeast precipitation-limited region and southeast cold temperature-limited region. Core-to-edge patterns of genetic diversity and differentiation were consistent with the CMH in the southeast, but were not supported in the northeast and showed mixed support in the northwest. Results suggest cold temperatures are a likely contributor to southeastern range limits, consistent with CMH predictions. In the northeast and northwest, ecological processes consisting of a steep physiological barrier and ongoing invasion dynamics, respectively, are more likely explanations for population genomic patterns than the CMH.

Genome-wide association analyses are a powerful approach for identifying gene function. These analyses are becoming commonplace in studies of humans, domesticated animals, and crop plants but have rarely been conducted in bacteria. We applied association analyses to 20 traits measured in Ensifer meliloti, an agriculturally and ecologically important bacterium because it fixes nitrogen when in symbiosis with leguminous plants. We identified candidate alleles and gene presence-absence variants underlying variation in symbiosis traits, antibiotic resistance, and use of various carbon sources; some of these candidates are in genes previously known to affect these traits whereas others were in genes that have not been well characterized. Our results point to the potential power of association analyses in bacteria, but also to the need to carefully evaluate the potential for false associations., Genome-wide association studies (GWAS) can identify genetic variants responsible for naturally occurring and quantitative phenotypic variation. Association studies therefore provide a powerful complement to approaches that rely on de novo mutations for characterizing gene function. Although bacteria should be amenable to GWAS, few GWAS have been conducted on bacteria, and the extent to which nonindependence among genomic variants (e.g., linkage disequilibrium [LD]) and the genetic architecture of phenotypic traits will affect GWAS performance is unclear. We apply association analyses to identify candidate genes underlying variation in 20 biochemical, growth, and symbiotic phenotypes among 153 strains of Ensifer meliloti. For 11 traits, we find genotype-phenotype associations that are stronger than expected by chance, with the candidates in relatively small linkage groups, indicating that LD does not preclude resolving association candidates to relatively small genomic regions. The significant candidates show an enrichment for nucleotide polymorphisms (SNPs) over gene presence-absence variation (PAV), and for five traits, candidates are enriched in large linkage groups, a possible signature of epistasis. Many of the variants most strongly associated with symbiosis phenotypes were in genes previously identified as being involved in nitrogen fixation or nodulation. For other traits, apparently strong associations were not stronger than the range of associations detected in permuted data. In sum, our data show that GWAS in bacteria may be a powerful tool for characterizing genetic architecture and identifying genes responsible for phenotypic variation. However, careful evaluation of candidates is necessary to avoid false signals of association. IMPORTANCE Genome-wide association analyses are a powerful approach for identifying gene function. These analyses are becoming commonplace in studies of humans, domesticated animals, and crop plants but have rarely been conducted in bacteria. We applied association analyses to 20 traits measured in Ensifer meliloti, an agriculturally and ecologically important bacterium because it fixes nitrogen when in symbiosis with leguminous plants. We identified candidate alleles and gene presence-absence variants underlying variation in symbiosis traits, antibiotic resistance, and use of various carbon sources; some of these candidates are in genes previously known to affect these traits whereas others were in genes that have not been well characterized. Our results point to the potential power of association analyses in bacteria, but also to the need to carefully evaluate the potential for false associations.

Ensifer meliloti (formerly Rhizobium meliloti and Sinorhizobium meliloti) is a model bacterium for understanding legume–rhizobial symbioses. The tripartite genome of E. meliloti consists of a chromosome, pSymA and pSymB, and in some instances strain-specific accessory plasmids. The majority of previous sequencing studies have relied on the use of assemblies generated from short read sequencing, which leads to gaps and assembly errors. Here we used PacBio-based, long-read assemblies and were able to assemble, de novo, complete circular replicons. In this study, we sequenced, de novo-assembled and analysed 10 E. meliloti strains. Sequence comparisons were also done with data from six previously published genomes. We identified genome differences between the replicons, including mol% G+C and gene content, nucleotide repeats, and transposon-associated loci. Additionally, genomic rearrangements both within and between replicons were identified, providing insight into evolutionary processes at the structural level. There were few cases of inter-replicon gene transfer of core genes between the main replicons. Accessory plasmids were more similar to pSymA than to either pSymB or the chromosome, with respect to gene content, transposon content and G+C content. In our population, the accessory plasmids appeared to share an open genome with pSymA, which contains many nodulation- and nitrogen fixation-related genes. This may explain previous observations that horizontal gene transfer has a greater effect on the content of pSymA than pSymB, or the chromosome, and why some rhizobia show unstable nodulation phenotypes on legume hosts.

Identifying the genetic architecture of complex phenotypes is a central goal of modern biology, particularly for disease-related traits. Genome-wide association methods are a classical approach for identifying the genomic basis of variation in disease phenotypes, but such analyses are particularly challenging in natural populations due to sample size difficulties. Extensive mark-recapture data, strong linkage disequilibrium, and a lethal transmissible cancer make the Tasmanian devil (Sarcophilus harrisii) an ideal model for such an association study. We used a RAD-capture approach to genotype 624 devils at ~16,000 loci and then used association analyses to assess the heritability of three cancer-related phenotypes: infection case-control (where cases were infected devils and controls were devils that were never infected), age of first infection, and survival following infection. The SNP array explained much of the phenotypic variance for female survival (>80%) and female case-control (>61%). We found that a few large-effect SNPs explained much of the variance for female survival (~5 SNPs explained >61% of the total variance) whereas more SNPs (~56) of smaller effect explained less of the variance for female case-control (~23% of the total variance). By contrast, these same SNPs did not account for a significant proportion of phenotypic variance in males, suggesting that the genetic bases of these traits and/or selection differ across sexes. Loci involved with cell adhesion and cell-cycle regulation underlay trait variation, suggesting that the devil immune system is rapidly evolving to recognize and potentially suppress cancer growth through these pathways. Overall, our study provided necessary data for genomics-based conservation and management in Tasmanian devils.

Herbivores often move among spatially interspersed host plants, tracking high-quality resources through space and time. This dispersal is of particular interest for vectors of plant pathogens. Existing molecular tools to track such movement have yielded important insights, but often provide insufficient genetic resolution to infer spread at finer spatiotemporal scales. Here, we explore the use of Nextera-tagmented reductively-amplified DNA (NextRAD) sequencing to infer movement of a highly-mobile winged insect, the potato psyllid (Bactericera cockerelli), among host plants. The psyllid vectors the pathogen that causes zebra chip disease in potato (Solanum tuberosum), but understanding and managing the spread of this pathogen is limited by uncertainty about the insect's host plant(s) outside of the growing season. We identified 1,978 polymorphic loci among psyllids separated spatiotemporally on potato or in patches of bittersweet nightshade (S. dulcumara), a weedy plant proposed to be the source of potato-colonizing psyllids. A subset of the psyllids on potato exhibited genetic similarity to insects on nightshade, consistent with regular movement between these two host plants. However, a second subset of potato-collected psyllids was genetically distinct from those collected on bittersweet nightshade; this suggests that a currently unrecognized source, i.e., other nightshade patches or a third host-plant species, could be contributing to psyllid populations in potato. Oftentimes, dispersal of vectors of pathogens must be tracked at a fine scale in order to understand, predict, and manage disease spread. We demonstrate that emerging sequencing technologies that detect genome-wide SNPs of a vector can be used to infer such localized movement.

Herbivores often move among spatially interspersed host plants, tracking high-quality resources through space and time. This dispersal is of particular interest for vectors of plant pathogens. Existing molecular tools to track such movement have yielded important insights, but often provide insufficient genetic resolution to infer spread at finer spatiotemporal scales. Here, we explore the use of Nextera-tagmented reductively-amplified DNA (NextRAD) sequencing to infer movement of a highly-mobile winged insect, the potato psyllid (Bactericera cockerelli), among host plants. The psyllid vectors the pathogen that causes zebra chip disease in potato (Solanum tuberosum), but understanding and managing the spread of this pathogen is limited by uncertainty about the insect’s host plant(s) outside of the growing season. We identified 8,443 polymorphic loci among psyllids separated spatiotemporally on potato or in patches of bittersweet nightshade (S. dulcumara), a weedy plant proposed to be the source of potato-colonizing psyllids. A subset of the psyllids on potato exhibited close genetic similarity to insects on nightshade, consistent with regular movement between these two host plants. However, a second subset of potato-collected psyllids was genetically distinct from those collected on bittersweet nightshade; this suggests that a currently unrecognized host-plant species could be contributing to psyllid populations in potato. Oftentimes, dispersal of vectors of plant or animal pathogens must be tracked at a relatively fine scale in order to understand, predict, and manage disease spread. We demonstrate that emerging sequencing technologies that detect SNPs across a vector’s entire genome can be used to infer such localized movement.

Ranaviruses, a genus of the Iridoviridae, are large double-stranded DNA viruses that infect cold-blooded vertebrates worldwide. Ranaviruses have caused severe epizootics in commercial frog and fish populations, and are currently classified as notifiable pathogens in international trade. Previous work shows that a ranavirus that infects tiger salamanders throughout Western North America (Ambystoma tigrinum virus, or ATV) is in high prevalence among salamanders in the fishing bait trade. Bait ATV strains have elevated virulence and are transported long distances by humans, providing widespread opportunities for pathogen pollution. We sequenced the genomes of 15 strains of ATV collected from tiger salamanders across western North America and performed phylogenetic and population genomic analyses and tests for recombination. We find that ATV forms a monophyletic clade within the rest of the Ranaviruses and that it likely emerged within the last several thousand years, before human activities influenced its spread. We also identify several genes under strong positive selection, some of which appear to be involved in viral virulence and/or host immune evasion. In addition, we provide support for the pathogen pollution hypothesis with evidence of recombination among ATV strains, and potential bait-endemic strain recombination.

Structural variation, including variation in gene copy number and presence or absence of genes, is a widespread and important source of genomic variation. We used whole-genome DNA sequences from 48 strains of Sinorhizobium (recently renamed Ensifer), including 20 strains of Sinorhizobium meliloti and 12 strains of S. medicae that were the focus of the analyses, to study the fitness effects of new structural variants created by duplication and horizontal gene transfer. We find that derived duplicated and horizontally transferred (HT) genes segregate at lower frequency than synonymous and nonsynonymous nucleotide variants in S. meliloti and S. medicae. Furthermore, the relative frequencies of different types of variants are more similar in S. medicae than in S. meliloti, the species with the larger effective population size. These results are consistent with the hypothesis that most duplications and HT genes have deleterious effects. Diversity of duplications, as measured by segregating duplicated genes per gene, is greater than nucleotide diversity, consistent with a high rate of duplication. Our results suggest that the vast majority of structural variants found among closely related bacterial strains are short-lived and unlikely to be involved in species-wide adaptation.

Background The sinorhizobia are amongst the most well studied members of nitrogen-fixing root nodule bacteria and contribute substantial amounts of fixed nitrogen to the biosphere. While the alfalfa symbiont Sinorhizobium meliloti RM 1021 was one of the first rhizobial strains to be completely sequenced, little information is available about the genomes of this large and diverse species group. Results Here we report the draft assembly and annotation of 48 strains of Sinorhizobium comprising five genospecies. While S. meliloti and S. medicae are taxonomically related, they displayed different nodulation patterns on diverse Medicago host plants, and have differences in gene content, including those involved in conjugation and organic sulfur utilization. Genes involved in Nod factor and polysaccharide biosynthesis, denitrification and type III, IV, and VI secretion systems also vary within and between species. Symbiotic phenotyping and mutational analyses indicated that some type IV secretion genes are symbiosis-related and involved in nitrogen fixation efficiency. Moreover, there is a correlation between the presence of type IV secretion systems, heme biosynthesis and microaerobic denitrification genes, and symbiotic efficiency. Conclusions Our results suggest that each Sinorhizobium strain uses a slightly different strategy to obtain maximum compatibility with a host plant. This large genome data set provides useful information to better understand the functional features of five Sinorhizobium species, especially compatibility in legume-Sinorhizobium interactions. The diversity of genes present in the accessory genomes of members of this genus indicates that each bacterium has adopted slightly different strategies to interact with diverse plant genera and soil environments.

International audience; The concept of a trade-off has long played a prominent role in understanding the evolution of organismal interactions such as mutualism, parasitism, and competition. Given the complexity inherent to interactions between different evolutionary entities, ecological factors may especially limit the power of trade-off models to predict evolutionary change. Here, we use four case studies to examine the importance of ecological context for the study of trade-offs in organismal interactions: (1) resource-based mutualisms, (2) parasite transmission and virulence, (3) plant biological invasions, and (4) host range evolution in parasites and parasitoids. In the first two case studies, mechanistic trade-off models have long provided a strong theoretical framework but face the challenge of testing assumptions under ecologically realistic conditions. Work under the second two case studies often has a strong ecological grounding, but faces challenges in identifying or quantifying the underlying genetic mechanism of the trade-off. Attention is given to recent studies that have bridged the gap between evolutionary mechanism and ecological realism. Finally, we explore the distinction between ecological factors that mask the underlying evolutionary trade-offs, and factors that actually change the trade-off relationship between fitness-related traits important to organismal interactions.

The symbiosis between rhizobial bacteria and legume plants has served as a model for investigating the genetics of nitrogen fixation and the evolution of facultative mutualism. We used deep sequence coverage (>100×) to characterize genomic diversity at the nucleotide level among 12 Sinorhizobium medicae and 32 S. meliloti strains. Although these species are closely related and share host plants, based on the ratio of shared polymorphisms to fixed differences we found that horizontal gene transfer (HGT) between these species was confined almost exclusively to plasmid genes. Three multi-genic regions that show the strongest evidence of HGT harbor genes directly involved in establishing or maintaining the mutualism with host plants. In both species, nucleotide diversity is 1.5–2.5 times greater on the plasmids than chromosomes. Interestingly, nucleotide diversity in S. meliloti but not S. medicae is highly structured along the chromosome – with mean diversity (θπ) on one half of the chromosome five times greater than mean diversity on the other half. Based on the ratio of plasmid to chromosome diversity, this appears to be due to severely reduced diversity on the chromosome half with less diversity, which is consistent with extensive hitchhiking along with a selective sweep. Frequency-spectrum based tests identified 82 genes with a signature of adaptive evolution in one species or another but none of the genes were identified in both species. Based upon available functional information, several genes identified as targets of selection are likely to alter the symbiosis with the host plant, making them attractive targets for further functional characterization., Author Summary Facultative mutualisms are relationships between two species that can live independently, but derive benefits when living together with their mutualistic partners. The facultative mutualism between rhizobial bacteria and legume plants contributes approximately half of all biologically fixed nitrogen, an essential plant nutrient, and is an important source of nitrogen to both natural and agricultural ecosystems. We resequenced the genomes of 44 strains of two closely related species of the genus Sinorhizobium that form facultative mutualisms with the model legme Medicago truncatula. These data provide one of the most complete examinations of genomic diversity segregating within microbial species that are not causative agents of human illness. Our analyses reveal that horizontal gene transfer, a common source of new genes in microbial species, disproportionately affects genes with direct roles in the rhizobia-plant symbiosis. Analyses of nucleotide diversity segregating within each species suggests that strong selection, along with genetic hitchhiking has sharply reduced diversity along an entire chromosome half in S. meliloti. Despite the two species' ecological similarity, we did not find evidence for selection acting on the same genetic targets. In addition to providing insight into the evolutionary history of rhizobial, this study shows the feasibility and potential power of applying population genomic analyses to microbial species.