Chromosome-level genome assemblies of two parasitoid biocontrol wasps reveal the parthenogenesis mechanism and an associated novel virus.

BACKGROUND: Biocontrol is a key technology for the control of pest species. Microctonus parasitoid wasps (Hymenoptera: Braconidae) have been released in Aotearoa New Zealand as biocontrol agents, targeting three different pest weevil species. Despite their value as biocontrol agents, no genome assemblies are currently available for these Microctonus wasps, limiting investigations into key biological differences between the different species and strains. METHODS AND FINDINGS: Here we present high-quality genomes for Microctonus hyperodae and Microctonus aethiopoides, assembled with short read sequencing and Hi-C scaffolding. These assemblies have total lengths of 106.7 Mb for M. hyperodae and 129.2 Mb for M. aethiopoides, with scaffold N50 values of 9 Mb and 23 Mb respectively. With these assemblies we investigated differences in reproductive mechanisms, and association with viruses between Microctonus wasps. Meiosis-specific genes are conserved in asexual Microctonus, with in-situ hybridisation validating expression of one of these genes in the ovaries of asexual Microctonus aethiopoides. This implies asexual reproduction in these Microctonus wasps involves meiosis, with the potential for sexual reproduction maintained. Investigation of viral gene content revealed candidate genes that may be involved in virus-like particle production in M. aethiopoides, as well as a novel virus infecting M. hyperodae, for which a complete genome was assembled. CONCLUSION AND SIGNIFICANCE: These are the first published genomes for Microctonus wasps which have been deployed as biocontrol agents, in Aotearoa New Zealand. These assemblies will be valuable resources for continued investigation and monitoring of these biocontrol systems. Understanding the biology underpinning Microctonus biocontrol is crucial if we are to maintain its efficacy, or in the case of M. hyperodae to understand what may have influenced the significant decline of biocontrol efficacy. The potential for sexual reproduction in asexual Microctonus is significant given that empirical modelling suggests this asexual reproduction is likely to have contributed to biocontrol decline. Furthermore the identification of a novel virus in M. hyperodae highlights a previously unknown aspect of this biocontrol system, which may contribute to premature mortality of the host pest. These findings have potential to be exploited in future in attempt to increase the effectiveness of M. hyperodae biocontrol.

Galba: genome annotation with miniprot and AUGUSTUS.

BACKGROUND: The Earth Biogenome Project has rapidly increased the number of available eukaryotic genomes, but most released genomes continue to lack annotation of protein-coding genes. In addition, no transcriptome data is available for some genomes. RESULTS: Various gene annotation tools have been developed but each has its limitations. Here, we introduce GALBA, a fully automated pipeline that utilizes miniprot, a rapid protein-to-genome aligner, in combination with AUGUSTUS to predict genes with high accuracy. Accuracy results indicate that GALBA is particularly strong in the annotation of large vertebrate genomes. We also present use cases in insects, vertebrates, and a land plant. GALBA is fully open source and available as a docker image for easy execution with Singularity in high-performance computing environments. CONCLUSIONS: Our pipeline addresses the critical need for accurate gene annotation in newly sequenced genomes, and we believe that GALBA will greatly facilitate genome annotation for diverse organisms.

Species-wide genomics of kakapo provides tools to accelerate recovery.

The kakapo is a critically endangered, intensively managed, long-lived nocturnal parrot endemic to Aotearoa New Zealand. We generated and analysed whole-genome sequence data for nearly all individuals living in early 2018 (169 individuals) to generate a high-quality species-wide genetic variant callset. We leverage extensive long-term metadata to quantify genome-wide diversity of the species over time and present new approaches using probabilistic programming, combined with a phenotype dataset spanning five decades, to disentangle phenotypic variance into environmental and genetic effects while quantifying uncertainty in small populations. We find associations for growth, disease susceptibility, clutch size and egg fertility within genic regions previously shown to influence these traits in other species. Finally, we generate breeding values to predict phenotype and illustrate that active management over the past 45 years has maintained both genome-wide diversity and diversity in breeding values and, hence, evolutionary potential. We provide new pathways for informing future conservation management decisions for kakapo, including prioritizing individuals for translocation and monitoring individuals with poor growth or high disease risk. Overall, by explicitly addressing the challenge of the small sample size, we provide a template for the inclusion of genomic data that will be transformational for species recovery efforts around the globe.

Germline progenitors and oocyte production in the honeybee queen ovary.

Understanding the reproduction of honeybee queens is crucial to support populations of this economically important insect. Here we examine the structure of the honeybee ovary to determine the nature of the germline progenitors in the ovary. Using a panel of marker genes that mark somatic or germline tissue in other insects we determine which cells in the honeybee ovary are somatic and which germline. We examine patterns of cell division and demonstrate that, unlike Drosophila, there is no evidence of single germline stem cells that provide the germline in honeybees. Germline progenitors are clustered in groups of 8 cells, joined by a polyfusome, and collections of these, in each ovariole, appear to maintain the germline during reproduction. We also show that these 8-cell clusters can divide and that their division occurs such that the numbers of germline progenitors are relatively constant over the reproductive life of queen honeybees. This information helps us to understand the diversity of structures in insect reproduction, and provide information to better support honeybee reproduction.

Assembly of female and male hihi genomes (stitchbird; Notiomystis cincta) enables characterization of the W chromosome and resources for conservation genomics.

A high-quality reference genome can be a valuable resource for threatened species by providing a foundation to assess their evolutionary potential to adapt to future pressures such as environmental change. We assembled the genome of a female hihi (Notiomysits cincta), a threatened passerine bird endemic to Aotearoa New Zealand. The assembled genome is 1.06 Gb, and is of high quality and highly contiguous, with a contig N50 of 7.0 Mb, estimated QV of 44 and a BUSCO completeness of 96.8%. A male assembly of comparable quality was generated in parallel. A population linkage map was used to scaffold the autosomal contigs into chromosomes. Female and male sequence coverage and comparative genomics analyses were used to identify Z-, and W-linked contigs. In total, 94.6% of the assembly length was assigned to putative nuclear chromosome scaffolds. Native DNA methylation was highly correlated between sexes, with the W chromosome contigs more highly methylated than autosomal chromosomes and Z contigs. 43 differentially methylated regions were identified, and these may represent interesting candidates for the establishment or maintenance of sex differences. By generating a high-quality reference assembly of the heterogametic sex, we have created a resource that enables characterization of genome-wide diversity and facilitates the investigation of female-specific evolutionary processes. The reference genomes will form the basis for fine-scale assessment of the impacts of low genetic diversity and inbreeding on the adaptive potential of the species and will therefore enable tailored and informed conservation management of this threatened taonga (treasured) species.

GALBA: Genome Annotation with Miniprot and AUGUSTUS.

The Earth Biogenome Project has rapidly increased the number of available eukaryotic genomes, but most released genomes continue to lack annotation of protein-coding genes. In addition, no transcriptome data is available for some genomes. Various gene annotation tools have been developed but each has its limitations. Here, we introduce GALBA, a fully automated pipeline that utilizes miniprot, a rapid protein- to-genome aligner, in combination with AUGUSTUS to predict genes with high accuracy. Accuracy results indicate that GALBA is particularly strong in the annotation of large vertebrate genomes. We also present use cases in insects, vertebrates, and a previously unannotated land plant. GALBA is fully open source and available as a docker image for easy execution with Singularity in high-performance computing environments. Our pipeline addresses the critical need for accurate gene annotation in newly sequenced genomes, and we believe that GALBA will greatly facilitate genome annotation for diverse organisms.

The promise and challenges of characterizing genome-wide structural variants: A case study in a critically endangered parrot.

There is growing interest in the role of structural variants (SVs) as drivers of local adaptation and speciation. From a biodiversity genomics perspective, the characterization of genome-wide SVs provides an exciting opportunity to complement single nucleotide polymorphisms (SNPs). However, little is known about the impacts of SV discovery and genotyping strategies on the characterization of genome-wide SV diversity within and among populations. Here, we explore a near whole-species resequence data set, and long-read sequence data for a subset of highly represented individuals in the critically endangered kakapo (Strigops habroptilus). We demonstrate that even when using a highly contiguous reference genome, different discovery and genotyping strategies can significantly impact the type, size and location of SVs characterized genome-wide. Further, we found that the mean number of SVs in each of two kakapo lineages differed both within and across generations. These combined results suggest that genome-wide characterization of SVs remains challenging at the population-scale. We are optimistic that increased accessibility to long-read sequencing and advancements in bioinformatic approaches including multireference approaches like genome graphs will alleviate at least some of the challenges associated with resolving SV characteristics below the species level. In the meantime, we address caveats, highlight considerations, and provide recommendations for the characterization of genome-wide SVs in biodiversity genomic research.

Hidden impacts of conservation management on fertility of the critically endangered kakapo.

Background: Animal conservation often requires intensive management actions to improve reproductive output, yet any adverse effects of these may not be immediately apparent, particularly in threatened species with small populations and long lifespans. Hand-rearing is an example of a conservation management strategy which, while boosting populations, can cause long-term demographic and behavioural problems. It is used in the recovery of the critically endangered kakapo (Strigops habroptilus), a flightless parrot endemic to New Zealand, to improve the slow population growth that is due to infrequent breeding, low fertility and low hatching success. Methods: We applied Bayesian mixed models to examine whether hand-rearing and other factors were associated with clutch fertility in kakapo. We used projection predictive variable selection to compare the relative contributions to fertility from the parents' rearing environment, their age and previous copulation experience, the parental kinship, and the number of mates and copulations for each clutch. We also explored how the incidence of repeated copulations and multiple mates varied with kakapo density. Results: The rearing status of the clutch father and the number of mates and copulations of the clutch mother were the dominant factors in predicting fertility. Clutches were less likely to be fertile if the father was hand-reared compared to wild-reared, but there was no similar effect for mothers. Clutches produced by females copulating with different males were more likely to be fertile than those from repeated copulations with one male, which in turn had a higher probability of fertility than those from a single copulation. The likelihood of multiple copulations and mates increased with female:male adult sex ratio, perhaps as a result of mate guarding by females. Parental kinship, copulation experience and age all had negligible associations with clutch fertility. Conclusions: These results provide a rare assessment of factors affecting fertility in a wild threatened bird species, with implications for conservation management. The increased fertility due to multiple mates and copulations, combined with the evidence for mate guarding and previous results of kakapo sperm morphology, suggests that an evolutionary mechanism exists to optimise fertility through sperm competition in kakapo. The high frequency of clutches produced from single copulations in the contemporary population may therefore represent an unnatural state, perhaps due to too few females. This suggests that opportunity for sperm competition should be maximised by increasing population densities, optimising sex ratios, and using artificial insemination. The lower fertility of hand-reared males may result from behavioural defects due to lack of exposure to conspecifics at critical development stages, as seen in other taxa. This potential negative impact of hand-rearing must be balanced against the short-term benefits it provides.

Genomic signatures of parallel alpine adaptation in recently evolved flightless insects.

Natural selection along elevational gradients has potential to drive predictable adaptations across distinct lineages, but the extent of such repeated evolution remains poorly studied for many widespread alpine taxa. We present parallel genomic analyses of two recently evolved flightless alpine insect lineages to test for molecular signatures of repeated alpine adaptation. Specifically, we compare low-elevation vs. alpine stonefly ecotypes from parallel stream populations in which flightless upland ecotypes have been independently derived. We map 67,922 polymorphic genetic markers, generated across 176 Zelandoperla fenestrata specimens from two independent alpine stream populations in New Zealand's Rock and Pillar Range, to a newly developed plecopteran reference genome. Genome-wide scans revealed 31 regions with outlier single nucleotide polymorphisms (SNPs) differentiating lowland vs. alpine ecotypes in Lug Creek, and 37 regions with outliers differentiating ecotypes in Six Mile Creek. Of these regions, 13% (8/60) yielded outlier SNPs across both within-stream ecotype comparisons, implying comparable genomic shifts contribute to this repeated alpine adaptation. Candidate genes closely linked to repeated outlier regions include several with documented roles in insect wing-development (e.g., dishevelled), suggesting that they may contribute to repeated alpine wing reduction. Additional candidate genes have been shown to influence insect fecundity (e.g., ovo) and lifespan (e.g., Mrp4), implying that they might contribute to life history differentiation between upland and lowland ecotypes. Additional outlier genes have potential roles in the evolution of reproductive isolation among ecotypes (hedgehog and Desaturase 1). These results demonstrate how replicated outlier tests across independent lineages can potentially contribute to the discovery of genes underpinning repeated adaptation.

Expanding the conservation genomics toolbox: Incorporating structural variants to enhance genomic studies for species of conservation concern.

Structural variants (SVs) are large rearrangements (>50 bp) within the genome that impact gene function and the content and structure of chromosomes. As a result, SVs are a significant source of functional genomic variation, that is, variation at genomic regions underpinning phenotype differences, that can have large effects on individual and population fitness. While there are increasing opportunities to investigate functional genomic variation in threatened species via single nucleotide polymorphism (SNP) data sets, SVs remain understudied despite their potential influence on fitness traits of conservation interest. In this future-focused Opinion, we contend that characterizing SVs offers the conservation genomics community an exciting opportunity to complement SNP-based approaches to enhance species recovery. We also leverage the existing literature-predominantly in human health, agriculture and ecoevolutionary biology-to identify approaches for readily characterizing SVs and consider how integrating these into the conservation genomics toolbox may transform the way we manage some of the world's most threatened species.

An exploration of assembly strategies and quality metrics on the accuracy of the rewarewa (Knightia excelsa) genome.

We used long read sequencing data generated from Knightia excelsa, a nectar-producing Proteaceae tree endemic to Aotearoa (New Zealand), to explore how sequencing data type, volume and workflows can impact final assembly accuracy and chromosome reconstruction. Establishing a high-quality genome for this species has specific cultural importance to Maori and commercial importance to honey producers in Aotearoa. Assemblies were produced by five long read assemblers using data subsampled based on read lengths, two polishing strategies and two Hi-C mapping methods. Our results from subsampling the data by read length showed that each assembler tested performed differently depending on the coverage and the read length of the data. Subsampling highlighted that input data with longer read lengths but perhaps lower coverage constructed more contiguous, kmers and gene-complete assemblies than short read length input data with higher coverage. The final genome assembly was constructed into 14 pseudochromosomes using an initial flye long read assembly, a racon/medaka/pilon combined polishing strategy, salsa2 and allhic scaffolding, juicebox curation, and Macadamia linkage map validation. We highlighted the importance of developing assembly workflows based on the volume and read length of sequencing data and established a robust set of quality metrics for generating high-quality assemblies. Scaffolding analyses highlighted that problems found in the initial assemblies could not be resolved accurately by Hi-C data and that assembly scaffolding was more successful when the underlying contig assembly was of higher accuracy. These findings provide insight into how quality assessment tools can be implemented throughout genome assembly pipelines to inform the de novo reconstruction of a high-quality genome assembly for nonmodel organisms.

Population genomics of the critically endangered kakapo.

The kakapo is a flightless parrot endemic to New Zealand. Once common in the archipelago, only 201 individuals remain today, most of them descending from an isolated island population. We report the first genome-wide analyses of the species, including a high-quality genome assembly for kakapo, one of the first chromosome-level reference genomes sequenced by the Vertebrate Genomes Project (VGP). We also sequenced and analyzed 35 modern genomes from the sole surviving island population and 14 genomes from the extinct mainland population. While theory suggests that such a small population is likely to have accumulated deleterious mutations through genetic drift, our analyses on the impact of the long-term small population size in kakapo indicate that present-day island kakapo have a reduced number of harmful mutations compared to mainland individuals. We hypothesize that this reduced mutational load is due to the island population having been subjected to a combination of genetic drift and purging of deleterious mutations, through increased inbreeding and purifying selection, since its isolation from the mainland ∼10,000 years ago. Our results provide evidence that small populations can survive even when isolated for hundreds of generations. This work provides key insights into kakapo breeding and recovery and more generally into the application of genetic tools in conservation efforts for endangered species.

Genomics Reveals Widespread Ecological Speciation in Flightless Insects.

Recent genomic analyses have highlighted parallel divergence in response to ecological gradients, but the extent to which altitude can underpin such repeated speciation remains unclear. Wing reduction and flight loss have apparently evolved repeatedly in montane insect assemblages and have been suggested as important drivers of hexapod diversification. We test this hypothesis using genomic analyses of a widespread wing-polymorphic stonefly species complex in New Zealand. We identified over 50,000 polymorphic genetic markers generated across almost 200 Zelandoperla fenestrata stonefly specimens using a newly generated plecopteran reference genome, to reveal widespread parallel speciation between sympatric full-winged and wing-reduced ecotypes. Rather than the existence of a single, widespread, flightless taxon (Zelandoperla pennulata), evolutionary genomic data reveal that wing-reduced upland lineages have speciated repeatedly and independently from full-winged Z. fenestrata. This repeated evolution of reproductive isolation between local ecotype pairs that lack mitochondrial DNA differentiation suggests that ecological speciation has evolved recently. A cluster of outlier single-nucleotide polymorphisms detected in independently wing-reduced lineages, tightly linked in an approximately 85 kb genomic region that includes the developmental "supergene" doublesex, suggests that this "island of divergence" may play a key role in rapid ecological speciation. [Ecological speciation; genome assembly; genomic island of differentiation; genotyping-by-sequencing; incipient species; plecoptera; wing reduction.].

High-Quality Assemblies for Three Invasive Social Wasps from the Vespula Genus.

Social wasps of the genus Vespula have spread to nearly all landmasses worldwide and have become significant pests in their introduced ranges, affecting economies and biodiversity. Comprehensive genome assemblies and annotations for these species are required to develop the next generation of control strategies and monitor existing chemical control. We sequenced and annotated the genomes of the common wasp (Vespula vulgaris), German wasp (Vespula germanica), and the western yellowjacket (Vespula pensylvanica). Our chromosome-level Vespula assemblies each contain 176-179 Mb of total sequence assembled into 25 scaffolds, with 10-200 unanchored scaffolds, and 16,566-18,948 genes. We annotated gene sets relevant to the applied management of invasive wasp populations, including genes associated with spermatogenesis and development, pesticide resistance, olfactory receptors, immunity and venom. These genomes provide evidence for active DNA methylation in Vespidae and tandem duplications of venom genes. Our genomic resources will contribute to the development of next-generation control strategies, and monitoring potential resistance to chemical control.

The potential for a CRISPR gene drive to eradicate or suppress globally invasive social wasps.

CRISPR gene drives have potential for widespread and cost-efficient pest control, but are highly controversial. We examined a potential gene drive targeting spermatogenesis to control the invasive common wasp (Vespula vulgaris) in New Zealand. Vespula wasps are haplodiploid. Their life cycle makes gene drive production challenging, as nests are initiated by single fertilized queens in spring followed by several cohorts of sterile female workers and the production of reproductives in autumn. We show that different spermatogenesis genes have different levels of variation between introduced and native ranges, enabling a potential 'precision drive' that could target the reduced genetic diversity and genotypes within the invaded range. In vitro testing showed guide-RNA target specificity and efficacy that was dependent on the gene target within Vespula, but no cross-reactivity in other Hymenoptera. Mathematical modelling incorporating the genetic and life history traits of Vespula wasps identified characteristics for a male sterility drive to achieve population control. There was a trade-off between drive infiltration and impact: a drive causing complete male sterility would not spread, while partial sterility could be effective in limiting population size if the homing rate is high. Our results indicate that gene drives may offer viable suppression for wasps and other haplodiploid pests.

The Pacific Biosciences de novo assembled genome dataset from a parthenogenetic New Zealand wild population of the longhorned tick, Haemaphysalis longicornis Neumann, 1901.

The longhorned tick, Haemaphysalis longicornis, feeds upon a wide range of bird and mammalian hosts. Mammalian hosts include cattle, deer, sheep, goats, humans, and horses. This tick is known to transmit a number of pathogens causing tick-borne diseases, and was the vector of a recent serious outbreak of oriental theileriosis in New Zealand. A New Zealand-USA consortium was established to sequence, assemble, and annotate the genome of this tick, using ticks obtained from New Zealand's North Island. In New Zealand, the tick is considered exclusively parthenogenetic and this trait was deemed useful for genome assembly. Very high molecular weight genomic DNA was sequenced on the Illumina HiSeq4000 and the long-read Pac Bio Sequel platforms. Twenty-eight SMRT cells produced a total of 21.3 million reads which were assembled with Canu on a reserved supercomputer node with access to 12 TB of RAM, running continuously for over 24 days. The final assembly dataset consisted of 34,211 contigs with an average contig length of 215,205 bp. The quality of the annotated genome was assessed by BUSCO analysis, an approach that provides quantitative measures for the quality of an assembled genome. Over 95% of the BUSCO gene set was found in the assembled genome. Only 48 of the 1066 BUSCO genes were missing and only 9 were present in a fragmented condition. The raw sequencing reads and the assembled contigs/scaffolds are archived at the National Center for Biotechnology Information.

Genome-Wide Association Analyses in the Model Rhizobium Ensifer meliloti.

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.

The complete replicons of 16 Ensifer meliloti strains offer insights into intra- and inter-replicon gene transfer, transposon-associated loci, and repeat elements.

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.

Select and resequence reveals relative fitness of bacteria in symbiotic and free-living environments.

Assays to accurately estimate relative fitness of bacteria growing in multistrain communities can advance our understanding of how selection shapes diversity within a lineage. Here, we present a variant of the "evolve and resequence" approach both to estimate relative fitness and to identify genetic variants responsible for fitness variation of symbiotic bacteria in free-living and host environments. We demonstrate the utility of this approach by characterizing selection by two plant hosts and in two free-living environments (sterilized soil and liquid media) acting on synthetic communities of the facultatively symbiotic bacterium Ensifer meliloti We find (i) selection that hosts exert on rhizobial communities depends on competition among strains, (ii) selection is stronger inside hosts than in either free-living environment, and (iii) a positive host-dependent relationship between relative strain fitness in multistrain communities and host benefits provided by strains in single-strain experiments. The greatest changes in allele frequencies in response to plant hosts are in genes associated with motility, regulation of nitrogen fixation, and host/rhizobia signaling. The approach we present provides a powerful complement to experimental evolution and forward genetic screens for characterizing selection in bacterial populations, identifying gene function, and surveying the functional importance of naturally occurring genomic variation.