Nonsynonymous Polymorphism Counts in Bacterial Genomes: a Comparative Examination

Reversions mask the contribution of adaptive evolution in microbiomes

When examining bacterial genomes for evidence of past selection, the results depend heavily on the mutational distance between chosen genomes. Even within a bacterial species, genomes separated by larger mutational distances exhibit stronger evidence of purifying selection as assessed by dN/dS, the normalized ratio of nonsynonymous to synonymous mutations. Here, we show that the classical interpretation of this scale dependence, weak purifying selection, leads to problematic mutation accumulation when applied to available gut microbiome data. We propose an alternative, adaptive reversion model with opposite implications for dynamical intuition and applications of dN/dS. Reversions that occur and sweep within-host populations are nearly guaranteed in microbiomes due to large population sizes, short generation times, and variable environments. Using analytical and simulation approaches, we show that adaptive reversion can explain the dN/dS decay given only dozens of locally fluctuating selective pressures, which is realistic in the context of Bacteroides genomes. The success of the adaptive reversion model argues for interpreting low values of dN/dS obtained from long timescales with caution as they may emerge even when adaptive sweeps are frequent. Our work thus inverts the interpretation of an old observation in bacterial evolution, illustrates the potential of mutational reversions to shape genomic landscapes over time, and highlights the importance of studying bacterial genomic evolution on short timescales.

Patterns of Change in Nucleotide Diversity Over Gene Length

Nucleotide diversity at a site is influenced by the relative strengths of neutral and selective population genetic processes. Therefore, attempts to estimate Effective population size based on the diversity of synonymous sites demand a better understanding of their selective constraints. The nucleotide diversity of a gene was previously found to correlate with its length. In this work, I measure nucleotide diversity at synonymous sites and uncover a pattern of low diversity towards the translation initiation site of a gene. The degree of reduction in diversity at the translation initiation site and the length of this region of reduced diversity can be quantified as "Effect Size" and "Effect Length" respectively, using parameters of an asymptotic regression model. Estimates of Effect Length across bacteria covaried with recombination rates as well as with a multitude of translation-associated traits such as the avoidance of mRNA secondary structure around translation initiation site, the number of rRNAs, and relative codon usage of ribosomal genes. Evolutionary simulations under purifying selection reproduce the observed patterns and diversity-length correlation and highlight that selective constraints on the 5'-region of a gene may be more extensive than previously believed. These results have implications for the estimation of effective population size, and relative mutation rates, and for genome scans of genes under positive selection based on "silent-site" diversity.

Patterns of change in nucleotide diversity over gene length

Nucleotide diversity at a site is influenced by the relative strengths of neutral and selective population genetic processes. Therefore, attempts to identify sites under positive selection require an understanding of the expected diversity in its absence. The nucleotide diversity of a gene was previously found to correlate with its length. In this work, I measure nucleotide diversity at synonymous sites and uncover a pattern of low diversity towards the translation initiation site (TIS) of a gene. The degree of reduction in diversity at the TIS and the length of this region of reduced diversity can be quantified as "Effect Size" and "Effect Length" respectively, using parameters of an asymptotic regression model. Estimates of Effect Length across bacteria covaried with recombination rates as well as with a multitude of fast-growth adaptations such as the avoidance of mRNA secondary structure around TIS, the number of rRNAs, and relative codon usage of ribosomal genes. Thus, the dependence of nucleotide diversity on gene length is governed by a combination of selective and non-selective processes. These results have implications for the estimation of effective population size and relative mutation rates based on "silent-site" diversity, and for pN/pS-based prediction of genes under selection.

Molecular Insight into Mycobacterium tuberculosis Resistance to Nitrofuranyl Amides Gained through Metagenomics-like Analysis of Spontaneous Mutants

We performed synthesis of new nitrofuranyl amides and investigated their anti-TB activity and primary genetic response of mycobacteria through whole-genome sequencing (WGS) of spontaneous resistant mutants. The in vitro activity was assessed on reference strain Mycobacterium tuberculosis H37Rv. The most active compound 11 was used for in vitro selection of spontaneous resistant mutants. The same mutations in six genes were detected in bacterial cultures grown under increased concentrations of 11 (2×, 4×, 8× MIC). The mutant positions were presented as mixed wild type and mutant alleles while increasing the concentration of the compound led to the semi-proportional and significant increase in mutant alleles. The identified genes belong to different categories and pathways. Some of them were previously reported as mediating drug resistance or drug tolerance, and counteracting oxidative and nitrosative stress, in particular: Rv0224c, fbiC, iniA, and Rv1592c. Gene-set interaction analysis revealed a certain weak interaction for gene pairs Rv1592–Rv1639c and Rv1592–Rv0224c. To conclude, this study experimentally demonstrated a multifaceted primary genetic response of M. tuberculosis to the action of nitrofurans. All three 11-treated subcultures independently presented the same six SNPs, which suggests their non-random occurrence and likely causative relationship between compound action and possible resistance mechanism.

Dissemination of a single ST11 clone of OXA-48-producing Klebsiella pneumoniae within a large polyclonal hospital outbreak determined by genomic sequencing

The population structure of a set of OXA-48-producing Klebsiella pneumoniae isolates belonging to sequence type 11 (ST11 Kp-OXA) and obtained from two hospitals in Madrid in the period from 2012 to 2015 was studied by genome sequencing. Overall, 97 ST11 Kp-OXA isolates were sequenced and their population structure and demography were studied by Bayesian phylodynamic analysis using core-genome SNVs. In total, 92 isolates were from Hospital La Paz, 57 of them from two selected units. The remaining five isolates were from different units of Hospital Doce de Octubre. Altogether, 96 out of the 97 ST11 Kp-OXA isolates could be ascribed to a single lineage that evolved into three sublineages. Demographic inference showed an expansion of the ST11 Kp-OXA in the first half of 2013 in agreement with the registered incidences. Dated phylogeny showed transmission clusters within hospital wards, between wards and between hospitals. The ST11 Kp-OXA outbreak in Hospital La Paz was largely due to the expansion of a single clone that was transmitted between different units and to Hospital Doce de Octubre. This clone diverged into three sub-lineages and spread out following a mixed mode of neutral core-genome evolution with some features of antibiotic selection, frequent large deletions and plasmid loss and gain events.