Parallel Evolution in Predatory Bdellovibrio sp. NC01 during Long-Term Coculture with a Single Prey Strain.

Experimental evolution provides a powerful tool for examining how Bdellovibrio evolves in response to unique selective pressures associated with its predatory lifestyle. We tested how Bdellovibrio sp. NC01 adapts to long-term coculture with Pseudomonas sp. NC02, which is less susceptible to predation compared to other Gram-negative bacteria. Analyzing six replicate Bdellovibrio populations across six time points spanning 40 passages and 2,880 h of coculture, we detected 30 to 40 new mutations in each population that exceeded a frequency of 5%. Nonsynonymous substitutions were the most abundant type of new mutation, followed by small indels and synonymous substitutions. After completing the final passage, we detected 20 high-frequency (>75%) mutations across all six evolved Bdellovibrio populations. Eighteen of these alter protein sequences, and most increased in frequency rapidly. Four genes acquired a high-frequency mutation in two or more evolved Bdellovibrio populations, reflecting parallel evolution and positive selection. The genes encode a sodium/phosphate cotransporter family protein (Bd2221), a metallophosphoesterase (Bd0054), a TonB family protein (Bd0396), and a hypothetical protein (Bd1601). Tested prey range and predation efficiency phenotypes did not differ significantly between evolved Bdellovibrio populations and the ancestor; however, all six evolved Bdellovibrio populations demonstrated enhanced starvation survival compared to the ancestor. These results suggest that, instead of evolving improved killing of Pseudomonas sp. NC02, Bdellovibrio evolved to better withstand nutrient limitation in the presence of this prey strain. The mutations identified here point to genes and functions that may be important for Bdellovibrio adaptation to the different selective pressures of long-term coculture with Pseudomonas. IMPORTANCE Bdellovibrio attack and kill Gram-negative bacteria, including drug-resistant pathogens of animals and plants. This lifestyle is unusual among bacteria, and it imposes unique selective pressures on Bdellovibrio. Determining how Bdellovibrio evolve in response to these pressures is valuable for understanding the mechanisms that govern predation. We applied experimental evolution to test how Bdellovibrio sp. NC01 evolved in response to long-term coculture with a single Pseudomonas strain, which NC01 can kill, but with low efficiency. Our experimental design imposed different selective pressures on the predatory bacteria and tracked the evolutionary trajectories of replicate Bdellovibrio populations. Using genome sequencing, we identified Bdellovibrio genes that acquired high-frequency mutations in two or more populations. Using phenotype assays, we determined that evolved Bdellovibrio populations did not improve their ability to kill Pseudomonas, but rather are better able to survive starvation. Overall, our results point to functions that may be important for Bdellovibrio adaptation.

Variation in genome content and predatory phenotypes between Bdellovibrio sp. NC01 isolated from soil and B. bacteriovorus type strain HD100.

Defining phenotypic and associated genotypic variation among Bdellovibrio may further our understanding of how this genus attacks and kills different Gram-negative bacteria. We isolated Bdellovibrio sp. NC01 from soil. Analysis of 16S rRNA gene sequences and average amino acid identity showed that NC01 belongs to a different species than the type species bacteriovorus. By clustering amino acid sequences from completely sequenced Bdellovibrio and comparing the resulting orthologue groups to a previously published analysis, we defined a 'core genome' of 778 protein-coding genes and identified four protein-coding genes that appeared to be missing only in NC01. To determine how horizontal gene transfer (HGT) may have impacted NC01 genome evolution, we performed genome-wide comparisons of Bdellovibrio nucleotide sequences, which indicated that eight NC01 genomic regions were likely acquired by HGT. To investigate how genome variation may impact predation, we compared protein-coding gene content between NC01 and the B. bacteriovorus type strain HD100, focusing on genes implicated as important in successful killing of prey. Of these, NC01 is missing ten genes that may play roles in lytic activity during predation. Compared to HD100, NC01 kills fewer tested prey strains and kills Escherichia coli ML35 less efficiently. NC01 causes a smaller log reduction in ML35, after which the prey population recovers and the NC01 population decreases. In addition, NC01 forms turbid plaques on lawns of E. coli ML35, in contrast to clear plaques formed by HD100. Linking phenotypic variation in interactions between Bdellovibrio and Gram-negative bacteria with underlying Bdellovibrio genome variation is valuable for understanding the ecological significance of predatory bacteria and evaluating their effectiveness in clinical applications.

CARD 2017: expansion and model-centric curation of the comprehensive antibiotic resistance database.

The Comprehensive Antibiotic Resistance Database (CARD; http://arpcard.mcmaster.ca) is a manually curated resource containing high quality reference data on the molecular basis of antimicrobial resistance (AMR), with an emphasis on the genes, proteins and mutations involved in AMR. CARD is ontologically structured, model centric, and spans the breadth of AMR drug classes and resistance mechanisms, including intrinsic, mutation-driven and acquired resistance. It is built upon the Antibiotic Resistance Ontology (ARO), a custom built, interconnected and hierarchical controlled vocabulary allowing advanced data sharing and organization. Its design allows the development of novel genome analysis tools, such as the Resistance Gene Identifier (RGI) for resistome prediction from raw genome sequence. Recent improvements include extensive curation of additional reference sequences and mutations, development of a unique Model Ontology and accompanying AMR detection models to power sequence analysis, new visualization tools, and expansion of the RGI for detection of emergent AMR threats. CARD curation is updated monthly based on an interplay of manual literature curation, computational text mining, and genome analysis.

Ras GTPase-like protein MglA, a controller of bacterial social-motility in Myxobacteria, has evolved to control bacterial predation by Bdellovibrio.

Bdellovibrio bacteriovorus invade Gram-negative bacteria in a predatory process requiring Type IV pili (T4P) at a single invasive pole, and also glide on surfaces to locate prey. Ras-like G-protein MglA, working with MglB and RomR in the deltaproteobacterium Myxococcus xanthus, regulates adventurous gliding and T4P-mediated social motility at both M. xanthus cell poles. Our bioinformatic analyses suggested that the GTPase activating protein (GAP)-encoding gene mglB was lost in Bdellovibrio, but critical residues for MglA(Bd) GTP-binding are conserved. Deletion of mglA(Bd) abolished prey-invasion, but not gliding, and reduced T4P formation. MglA(Bd) interacted with a previously uncharacterised tetratricopeptide repeat (TPR) domain protein Bd2492, which we show localises at the single invasive pole and is required for predation. Bd2492 and RomR also interacted with cyclic-di-GMP-binding receptor CdgA, required for rapid prey-invasion. Bd2492, RomR(Bd) and CdgA localize to the invasive pole and may facilitate MglA-docking. Bd2492 was encoded from an operon encoding a TamAB-like secretion system. The TamA protein and RomR were found, by gene deletion tests, to be essential for viability in both predatory and non-predatory modes. Control proteins, which regulate bipolar T4P-mediated social motility in swarming groups of deltaproteobacteria, have adapted in evolution to regulate the anti-social process of unipolar prey-invasion in the "lone-hunter" Bdellovibrio. Thus GTP-binding proteins and cyclic-di-GMP inputs combine at a regulatory hub, turning on prey-invasion and allowing invasion and killing of bacterial pathogens and consequent predatory growth of Bdellovibrio.

Large plasmids of Escherichia coli and Salmonella encode highly diverse arrays of accessory genes on common replicon families.

Plasmids are important in evolution and adaptation of host bacteria, yet we lack a comprehensive picture of their own natural variation. We used replicon typing and RFLP analysis to assess diversity and distribution of plasmids in the ECOR, SARA, SARB and SARC reference collections of Escherichia coli and Salmonella. Plasmids, especially large (≥30 kb) plasmids, are abundant in these collections. Host species and genotype clearly impact plasmid prevalence; plasmids are more abundant in ECOR than SAR, but, within ECOR, subgroup B2 strains have the fewest large plasmids. The majority of large plasmids have unique RFLP patterns, suggesting high variation, even within dominant replicon families IncF and IncI1. We found only four conserved plasmid types within ECOR, none of which are widely distributed. Within SAR, conserved plasmid types are primarily serovar-specific, including a pSLT-like plasmid in 13 Typhimurium strains. Conservation of pSLT contrasts with variability of other plasmids, suggesting evolution of serovar-specific virulence plasmids is distinct from that of most enterobacterial plasmids. We sequenced a conserved serovar Heidelberg plasmid but did not detect virulence or antibiotic resistance genes. Our data illustrate the high degree of natural variation in large plasmids of E. coli and Salmonella, even among plasmids sharing backbone genes.

Genome analysis of a simultaneously predatory and prey-independent, novel Bdellovibrio bacteriovorus from the River Tiber, supports in silico predictions of both ancient and recent lateral gene transfer from diverse bacteria.

BACKGROUND: Evolution equipped Bdellovibrio bacteriovorus predatory bacteria to invade other bacteria, digesting and replicating, sealed within them thus preventing nutrient-sharing with organisms in the surrounding environment. Bdellovibrio were previously described as "obligate predators" because only by mutations, often in gene bd0108, are 1 in ~1x10(7) of predatory lab strains of Bdellovibrio converted to prey-independent growth. A previous genomic analysis of B. bacteriovorus strain HD100 suggested that predatory consumption of prey DNA by lytic enzymes made Bdellovibrio less likely than other bacteria to acquire DNA by lateral gene transfer (LGT). However the Doolittle and Pan groups predicted, in silico, both ancient and recent lateral gene transfer into the B. bacteriovorus HD100 genome. RESULTS: To test these predictions, we isolated a predatory bacterium from the River Tiber- a good potential source of LGT as it is rich in diverse bacteria and organic pollutants- by enrichment culturing with E. coli prey cells. The isolate was identified as B. bacteriovorus and named as strain Tiberius. Unusually, this Tiberius strain showed simultaneous prey-independent growth on organic nutrients and predatory growth on live prey. Despite the prey-independent growth, the homolog of bd0108 did not have typical prey-independent-type mutations. The dual growth mode may reflect the high carbon content of the river, and gives B. bacteriovorus Tiberius extended non-predatory contact with the other bacteria present. The HD100 and Tiberius genomes were extensively syntenic despite their different cultured-terrestrial/freshly-isolated aquatic histories; but there were significant differences in gene content indicative of genomic flux and LGT. Gene content comparisons support previously published in silico predictions for LGT in strain HD100 with substantial conservation of genes predicted to have ancient LGT origins but little conservation of AT-rich genes predicted to be recently acquired. CONCLUSIONS: The natural niche and dual predatory, and prey-independent growth of the B. bacteriovorus Tiberius strain afforded it extensive non-predatory contact with other marine and freshwater bacteria from which LGT is evident in its genome. Thus despite their arsenal of DNA-lytic enzymes; Bdellovibrio are not always predatory in natural niches and their genomes are shaped by acquiring whole genes from other bacteria.

Complete Genome of Aeromonas encheleia Strain SOD01 Isolated from an Urban Freshwater Stream.

We isolated Aeromonas encheleia strain SOD01 from an urban freshwater stream in Providence, RI. De novo assembly of PacBio RSII data followed by polishing with Illumina MiSeq data generated a complete 4,450,115 bp genome with 61.8% GC content. PGAP annotation predicted 3,877 protein-coding genes, 127 tRNA, and 31 rRNA.

Chromosome and Plasmid Sequences of Pantoea sp. Strain SOD02 Isolated from an Urban Freshwater Stream.

After isolating Pantoea sp. strain SOD02 from an urban freshwater stream in Providence, RI, we used PacBio RSII data for de novo assembly and Illumina MiSeq data for polishing. This yielded complete circular sequences for a 4,227,027-bp chromosome with 54.7% GC and a 926,844-bp plasmid with 54.0% GC.

Unprecedented loss of ammonia assimilation capability in a urease-encoding bacterial mutualist.

BACKGROUND: Blochmannia are obligately intracellular bacterial mutualists of ants of the tribe Camponotini. Blochmannia perform key nutritional functions for the host, including synthesis of several essential amino acids. We used Illumina technology to sequence the genome of Blochmannia associated with Camponotus vafer. RESULTS: Although Blochmannia vafer retains many nutritional functions, it is missing glutamine synthetase (glnA), a component of the nitrogen recycling pathway encoded by the previously sequenced B. floridanus and B. pennsylvanicus. With the exception of Ureaplasma, B. vafer is the only sequenced bacterium to date that encodes urease but lacks the ability to assimilate ammonia into glutamine or glutamate. Loss of glnA occurred in a deletion hotspot near the putative replication origin. Overall, compared to the likely gene set of their common ancestor, 31 genes are missing or eroded in B. vafer, compared to 28 in B. floridanus and four in B. pennsylvanicus. Three genes (queA, visC and yggS) show convergent loss or erosion, suggesting relaxed selection for their functions. Eight B. vafer genes contain frameshifts in homopolymeric tracts that may be corrected by transcriptional slippage. Two of these encode DNA replication proteins: dnaX, which we infer is also frameshifted in B. floridanus, and dnaG. CONCLUSIONS: Comparing the B. vafer genome with B. pennsylvanicus and B. floridanus refines the core genes shared within the mutualist group, thereby clarifying functions required across ant host species. This third genome also allows us to track gene loss and erosion in a phylogenetic context to more fully understand processes of genome reduction.

Complete Genome Sequence of Bdellovibrio sp. Strain KM01, a Predatory Bacterium Isolated from Soil.

We report the complete genome sequence of the predatory bacterium Bdellovibrio sp. strain KM01, isolated from soil collected near a pond. The genome is 3,961,288 bp long with 45.5% GC content. Comparative genomics among Bdellovibrio strains will help us understand how genotypic differences affect differences in predatory phenotypes.

Over half of breakpoints in gene pairs involved in cancer-specific recurrent translocations are mapped to human chromosomal fragile sites.

BACKGROUND: Gene rearrangements such as chromosomal translocations have been shown to contribute to cancer development. Human chromosomal fragile sites are regions of the genome especially prone to breakage, and have been implicated in various chromosome abnormalities found in cancer. However, there has been no comprehensive and quantitative examination of the location of fragile sites in relation to all chromosomal aberrations. RESULTS: Using up-to-date databases containing all cancer-specific recurrent translocations, we have examined 444 unique pairs of genes involved in these translocations to determine the correlation of translocation breakpoints and fragile sites in the gene pairs. We found that over half (52%) of translocation breakpoints in at least one gene of these gene pairs are mapped to fragile sites. Among these, we examined the DNA sequences within and flanking three randomly selected pairs of translocation-prone genes, and found that they exhibit characteristic features of fragile DNA, with frequent AT-rich flexibility islands and the potential of forming highly stable secondary structures. CONCLUSION: Our study is the first to examine gene pairs involved in all recurrent chromosomal translocations observed in tumor cells, and to correlate the location of more than half of breakpoints to positions of known fragile sites. These results provide strong evidence to support a causative role for fragile sites in the generation of cancer-specific chromosomal rearrangements.

Antimicrobial Resistance Genes, Cassettes, and Plasmids Present in Salmonella enterica Associated With United States Food Animals.

The ability of antimicrobial resistance (AR) to transfer, on mobile genetic elements (MGEs) between bacteria, can cause the rapid establishment of multidrug resistance (MDR) in bacteria from animals, thus creating a foodborne risk to human health. To investigate MDR and its association with plasmids in Salmonella enterica, whole genome sequence (WGS) analysis was performed on 193 S. enterica isolated from sources associated with United States food animals between 1998 and 2011; 119 were resistant to at least one antibiotic tested. Isolates represented 86 serotypes and variants, as well as diverse phenotypic resistance profiles. A total of 923 AR genes and 212 plasmids were identified among the 193 strains. Every isolate contained at least one AR gene. At least one plasmid was detected in 157 isolates. Genes were identified for resistance to aminoglycosides (n = 472), ß-lactams (n = 84), tetracyclines (n = 171), sulfonamides (n = 91), phenicols (n = 42), trimethoprim (n = 8), macrolides (n = 5), fosfomycin (n = 48), and rifampicin (n = 2). Plasmid replicon types detected in the isolates were A/C (n = 32), ColE (n = 76), F (n = 43), HI1 (n = 4), HI2 (n = 20), I1 (n = 62), N (n = 4), Q (n = 7), and X (n = 35). Phenotypic resistance correlated with the AR genes identified in 95.4% of cases. Most AR genes were located on plasmids, with many plasmids harboring multiple AR genes. Six antibiotic resistance cassette structures (ARCs) and one pseudo-cassette were identified. ARCs contained between one and five resistance genes (ARC1: sul2, strAB, tetAR; ARC2: aac3-iid; ARC3: aph, sph; ARC4: cmy-2; ARC5: floR; ARC6: tetB; pseudo-ARC: aadA, aac3-VIa, sul1). These ARCs were present in multiple isolates and on plasmids of multiple replicon types. To determine the current distribution and frequency of these ARCs, the public NCBI database was analyzed, including WGS data on isolates collected by the USDA Food Safety and Inspection Service (FSIS) from 2014 to 2018. ARC1, ARC4, and ARC5 were significantly associated with cattle isolates, while ARC6 was significantly associated with chicken isolates. This study revealed that a diverse group of plasmids, carrying AR genes, are responsible for the phenotypic resistance seen in Salmonella isolated from United States food animals. It was also determined that many plasmids carry similar ARCs.

Bacterial communities in the digester bed and liquid effluent of a microflush composting toilet system.

Lack of access to clean water and sanitation is a major factor impacting public health in communities worldwide. To address this, the S-Lab at Providence College and the Global Sustainable Aid Project developed a microflush composting toilet system to isolate and treat human waste. Solid waste is composted within a filter-digester bed via an aerobic process involving microbes and invertebrates. Liquid waste may be sanitized by solar disinfection (SODIS) or slow sand filtration (SSF). Here, we used 16S rRNA amplicon sequencing of samples from a scaled-down test version of the system to better understand the bacterial component of the toilet system. Immediately after fecal matter was deposited in the test system, the bacterial community of the filter-digester bed at the site of deposition resembled that of the human gut at both the phylum and genus level, which was expected. Genus-level analysis of filter-digester bed samples collected over the next 30 days from the site of deposition showed reduced or undetectable levels of fecal-associated taxa, with the exception of Clostridium XI, which persisted at low abundance throughout the sampling period. Starting with the sample collected on day 4, the bacterial community of the filter-digester bed at the site of deposition was dominated by bacterial taxa commonly associated with environmental sources, reflecting a major shift in bacterial community composition. These data support the toilet system's capacity for processing solid human waste. We also analyzed how SODIS and SSF sanitization methods affected the bacterial community composition of liquid effluent collected on day 15 from the test system. Untreated and treated liquid effluent samples were dominated by Proteobacteria. At the genus level, the bacterial community of the untreated effluent included taxa commonly associated with environmental sources. In the SODIS-treated effluent, these genera increased in abundance, whereas in the SSF-treated effluent, they were greatly reduced or undetectable. By analyzing operational taxonomic units that were unclassified at the genus level, we observed that SSF appears to introduce new taxa into the treated effluent, likely from the biological film of microbes and small animals that constitutes the key element of SSF. These data will inform continued development of liquid waste handling strategies for the toilet system. Using the test system as an indicator of the performance of the full-scale version, we have shown the effectiveness of the microflush composting toilet system for containing and eliminating gut-associated bacteria, thereby improving sanitation and contributing to better public health in rural and peri-urban communities.

Prey Range and Genome Evolution of Halobacteriovorax marinus Predatory Bacteria from an Estuary.

Halobacteriovorax strains are saltwater-adapted predatory bacteria that attack Gram-negative bacteria and may play an important role in shaping microbial communities. To understand how Halobacteriovorax strains impact ecosystems and develop them as biocontrol agents, it is important to characterize variation in predation phenotypes and investigate Halobacteriovorax genome evolution. We isolated Halobacteriovorax marinus BE01 from an estuary in Rhode Island using Vibrio from the same site as prey. Small, fast-moving, attack-phase BE01 cells attach to and invade prey cells, consistent with the intraperiplasmic predation strategy of the H. marinus type strain, SJ. BE01 is a prey generalist, forming plaques on Vibrio strains from the estuary, Pseudomonas from soil, and Escherichia coli. Genome analysis revealed extremely high conservation of gene order and amino acid sequences between BE01 and SJ, suggesting strong selective pressure to maintain the genome in this H. marinus lineage. Despite this, we identified two regions of gene content difference that likely resulted from horizontal gene transfer. Analysis of modal codon usage frequencies supports the hypothesis that these regions were acquired from bacteria with different codon usage biases than H. marinus. In one of these regions, BE01 and SJ carry different genes associated with mobile genetic elements. Acquired functions in BE01 include the dnd operon, which encodes a pathway for DNA modification, and a suite of genes involved in membrane synthesis and regulation of gene expression that was likely acquired from another Halobacteriovorax lineage. This analysis provides further evidence that horizontal gene transfer plays an important role in genome evolution in predatory bacteria. IMPORTANCE Predatory bacteria attack and digest other bacteria and therefore may play a role in shaping microbial communities. To investigate phenotypic and genotypic variation in saltwater-adapted predatory bacteria, we isolated Halobacteriovorax marinus BE01 from an estuary in Rhode Island, assayed whether it could attack different prey bacteria, and sequenced and analyzed its genome. We found that BE01 is a prey generalist, attacking bacteria from different phylogenetic groups and environments. Gene order and amino acid sequences are highly conserved between BE01 and the H. marinus type strain, SJ. By comparative genomics, we detected two regions of gene content difference that likely occurred via horizontal gene transfer events. Acquired genes encode functions such as modification of DNA, membrane synthesis and regulation of gene expression. Understanding genome evolution and variation in predation phenotypes among predatory bacteria will inform their development as biocontrol agents and clarify how they impact microbial communities.

Complete Genome Sequence of Pseudomonas sp. Strain NC02, Isolated from Soil.

We report here the complete genome sequence of Pseudomonas sp. strain NC02, isolated from soil in eastern Massachusetts. We assembled PacBio reads into a single closed contig with 132× mean coverage and then polished this contig using Illumina MiSeq reads, yielding a 6,890,566-bp sequence with 61.1% GC content.

Complete Genome Sequence of Escherichia coli ML35.

We report here the complete genome sequence of Escherichia coli strain ML35. We assembled PacBio reads into a single closed contig with 169× mean coverage and then polished this contig using Illumina MiSeq reads, yielding a 4,918,774-bp sequence with 50.8% GC content.

Complete Genome Sequence of the Predatory Bacterium Ensifer adhaerens Casida A.

We report here the complete genome sequence of the facultative predatory bacterium Ensifer adhaerens strain Casida A. The genome was assembled into three circular contigs, with a main chromosome as well as two large secondary replicons, that totaled 7,267,502 bp with 6,641 predicted open reading frames.

Genome evolution in an ancient bacteria-ant symbiosis: parallel gene loss among Blochmannia spanning the origin of the ant tribe Camponotini.

Stable associations between bacterial endosymbionts and insect hosts provide opportunities to explore genome evolution in the context of established mutualisms and assess the roles of selection and genetic drift across host lineages and habitats. Blochmannia, obligate endosymbionts of ants of the tribe Camponotini, have coevolved with their ant hosts for ∼40 MY. To investigate early events in Blochmannia genome evolution across this ant host tribe, we sequenced Blochmannia from two divergent host lineages, Colobopsis obliquus and Polyrhachis turneri, and compared them with four published genomes from Blochmannia of Camponotus sensu stricto. Reconstructed gene content of the last common ancestor (LCA) of these six Blochmannia genomes is reduced (690 protein coding genes), consistent with rapid gene loss soon after establishment of the symbiosis. Differential gene loss among Blochmannia lineages has affected cellular functions and metabolic pathways, including DNA replication and repair, vitamin biosynthesis and membrane proteins. Blochmannia of P. turneri (i.e., B. turneri) encodes an intact DnaA chromosomal replication initiation protein, demonstrating that loss of dnaA was not essential for establishment of the symbiosis. Based on gene content, B. obliquus and B. turneri are unable to provision hosts with riboflavin. Of the six sequenced Blochmannia, B. obliquus is the earliest diverging lineage (i.e., the sister group of other Blochmannia sampled) and encodes the fewest protein-coding genes and the most pseudogenes. We identified 55 genes involved in parallel gene loss, including glutamine synthetase, which may participate in nitrogen recycling. Pathways for biosynthesis of coenzyme A, terpenoids and riboflavin were lost in multiple lineages, suggesting relaxed selection on the pathway after inactivation of one component. Analysis of Illumina read datasets did not detect evidence of plasmids encoding missing functions, nor the presence of coresident symbionts other than Wolbachia. Although gene order is strictly conserved in four Blochmannia of Camponotus sensu stricto, comparisons with deeply divergent lineages revealed inversions in eight genomic regions, indicating ongoing recombination despite ancestral loss of recA. In sum, the addition of two Blochmannia genomes of divergent host lineages enables reconstruction of early events in evolution of this symbiosis and suggests that Blochmannia lineages may experience distinct, host-associated selective pressures. Understanding how evolutionary forces shape genome reduction in this system may help to clarify forces driving gene loss in other bacteria, including intracellular pathogens.

Sequence context of indel mutations and their effect on protein evolution in a bacterial endosymbiont.

Indel mutations play key roles in genome and protein evolution, yet we lack a comprehensive understanding of how indels impact evolutionary processes. Genome-wide analyses enabled by next-generation sequencing can clarify the context and effect of indels, thereby integrating a more detailed consideration of indels with our knowledge of nucleotide substitutions. To this end, we sequenced Blochmannia chromaiodes, an obligate bacterial endosymbiont of carpenter ants, and compared it with the close relative, B. pennsylvanicus. The genetic distance between these species is small enough for accurate whole genome alignment but large enough to provide a meaningful spectrum of indel mutations. We found that indels are subjected to purifying selection in coding regions and even intergenic regions, which show a reduced rate of indel base pairs per kilobase compared with nonfunctional pseudogenes. Indels occur almost exclusively in repeat regions composed of homopolymers and multimeric simple sequence repeats, demonstrating the importance of sequence context for indel mutations. Despite purifying selection, some indels occur in protein-coding genes. Most are multiples of three, indicating selective pressure to maintain the reading frame. The deleterious effect of frameshift-inducing indels is minimized by either compensation from a nearby indel to restore reading frame or the indel's location near the 3'-end of the gene. We observed amino acid divergence exceeding nucleotide divergence in regions affected by frameshift-inducing indels, suggesting that these indels may either drive adaptive protein evolution or initiate gene degradation. Our results shed light on how indel mutations impact processes of molecular evolution underlying endosymbiont genome evolution.

Purifying selection, sequence composition, and context-specific indel mutations shape intraspecific variation in a bacterial endosymbiont.

Comparative genomics of closely related bacterial strains can clarify mutational processes and selective forces that impact genetic variation. Among primary bacterial endosymbionts of insects, such analyses have revealed ongoing genome reduction, raising questions about the ultimate evolutionary fate of these partnerships. Here, we explored genomic variation within Blochmannia vafer, an obligate mutualist of the ant Camponotus vafer. Polymorphism analysis of the Illumina data set used previously for de novo assembly revealed a second Bl. vafer genotype. To determine why a single ant colony contained two symbiont genotypes, we examined polymorphisms in 12 C. vafer mitochondrial sequences assembled from the Illumina data; the spectrum of variants suggests that the colony contained two maternal lineages, each harboring a distinct Bl. vafer genotype. Comparing the two Bl. vafer genotypes revealed that purifying selection purged most indels and nonsynonymous differences from protein-coding genes. We also discovered that indels occur frequently in multimeric simple sequence repeats, which are relatively abundant in Bl. vafer and may play a more substantial role in generating variation in this ant mutualist than in the aphid endosymbiont Buchnera. Finally, we explored how an apparent relocation of the origin of replication in Bl. vafer and the resulting shift in strand-associated mutational pressures may have caused accelerated gene loss and an elevated rate of indel polymorphisms in the region spanning the origin relocation. Combined, these results point to significant impacts of purifying selection on genomic polymorphisms as well as distinct patterns of indels associated with unusual genomic features of Blochmannia.