Putative Compensatory Mutations Research Articles

Compensatory mutations are associated with increased in vitro growth in resistant clinical samples of Mycobacterium tuberculosis.

Mutations in Mycobacterium tuberculosis associated with resistance to antibiotics often come with a fitness cost for the bacteria. Resistance to the first-line drug rifampicin leads to lower competitive fitness of M. tuberculosis populations when compared to susceptible populations. This fitness cost, introduced by resistance mutations in the RNA polymerase, can be alleviated by compensatory mutations (CMs) in other regions of the affected protein. CMs are of particular interest clinically since they could lock in resistance mutations, encouraging the spread of resistant strains worldwide. Here, we report the statistical inference of a comprehensive set of CMs in the RNA polymerase of M. tuberculosis, using over 70 000 M. tuberculosis genomes that were collated as part of the CRyPTIC project. The unprecedented size of this data set gave the statistical tests more power to investigate the association of putative CMs with resistance-conferring mutations. Overall, we propose 51 high-confidence CMs by means of statistical association testing and suggest hypotheses for how they exert their compensatory mechanism by mapping them onto the protein structure. In addition, we were able to show an association of CMs with higher in vitro growth densities, and hence presumably with higher fitness, in resistant samples in the more virulent M. tuberculosis lineage 2. Our results suggest the association of CM presence with significantly higher in vitro growth than for wild-type samples, although this association is confounded with lineage and sub-lineage affiliation. Our findings emphasize the integral role of CMs and lineage affiliation in resistance spread and increases the urgency of antibiotic stewardship, which implies accurate, cheap and widely accessible diagnostics for M. tuberculosis infections to not only improve patient outcomes but also prevent the spread of resistant strains.

Novel and reported compensatory mutations in rpoABC genes found in drug resistant tuberculosis outbreaks.

Rifampicin (RIF) is a key first-line drug used to treat tuberculosis, a primarily pulmonary disease caused by Mycobacterium tuberculosis. RIF resistance is caused by mutations in rpoB, at the cost of slower growth and reduced transcription efficiency. Antibiotic resistance to RIF is prevalent despite this fitness cost. Compensatory mutations in rpoABC genes have been shown to alleviate the fitness cost of rpoB:S450L, explaining how RIF resistant strains harbor this mutation can spread so rapidly. Unfortunately, the full set of RIF compensatory mutations is still unknown, particularly those compensating for rarer RIF resistance mutations. We performed an association study on a globally representative set of 4,309 whole genome sequenced clinical M. tuberculosis isolates to identify novel putative compensatory mutations, determine the prevalence of known and previously reported putative compensatory mutations, and determine which RIF resistance markers associate with these compensatory mutations. Of the 1,079 RIF resistant isolates, 638 carried previously reported putative and high-probability compensatory mutations. Our strict criteria identified 46 additional mutations in rpoABC for which no strong prior evidence of their compensatory role exists. Of these, 35 have previously been reported. As such, our independent corroboration adds to the mounting evidence that these 35 also carry a compensatory role. The remaining 11 are novel putative compensatory markers, reported here for the first time. Six of these 11 novel putative compensatory mutations had two or more mutation events. Most compensatory mutations appear to be specifically compensating for the fitness loss due to rpoB:S450L. However, an outbreak of 22 closely related isolates each carried three rpoB mutations, the rare RIFR markers D435G and L452P and the putative compensatory mutation I1106T. This suggests compensation may require specific combinations of rpoABC mutations. Here, we report only mutations that met our very strict criteria. It is highly likely that many additional rpoABC mutations compensate for rare resistance-causing mutations and therefore did not carry the statistical power to be reported here. These findings aid in the identification of RIF resistant M. tuberculosis strains with restored fitness, which pose a greater risk of causing resistant outbreaks.

Genomic analysis of lineage-specific transmission of multidrug resistance tuberculosis in China

ABSTRACT Objectives We investigated the genetic diversities and lineage-specific transmission dynamics of multidrug-resistant tuberculosis (MDR-TB), with the goal of determining the potential factors driving the MDR epidemics in China. Methods We curated a large nationwide Mycobacterium tuberculosis (M. tuberculosis) whole genome sequence data set, including 1313 MDR strains. We reconstructed the phylogeny and mapped the transmission networks of MDR-TB across China using Bayesian inference. To identify drug-resistance variants linked to enhanced transmissibility, we employed ordinary least-squares (OLS) regression analysis. Result The majority of MDR-TB strains in China belong to lineage 2.2.1. Transmission chain analysis has indicated that the repeated and frequent transmission of L2.2.1 plays a central role in the establishment of MDR epidemic in China, but no occurrence of a large predominant MDR outbreak was detected. Using OLS regression, the most common single nucleotide polymorphisms (SNPs) associated with resistance to isoniazid (katG_p.Ser315Thr and katG_p.Ser315Asn) and rifampicin (rpoB_p.Ser450Leu, rpoB_p.His445Tyr, rpoB_p.His445Arg, rpoB_p.His445Asp, and rpoB_p.His445Asn) were more likely to be found in L2 clustered strains. Several putative compensatory mutations in rpoA, rpoC, and katG were significantly associated with clustering. The eastern, central, and southern regions of China had a high level of connectivity for the migration of L2 MDR strains throughout the country. The skyline plot showed distinct population size expansion dynamics for MDR-TB lineages in China. Conclusion MDR-TB epidemic in China is predominantly driven by the spread of highly transmissible Beijing strains. A range of drug-resistance mutations of L2 MDR-TB strains displayed minimal fitness costs and may facilitate their transmission.

Genomic Characterization of Drug-Resistant Mycobacterium tuberculosis L2/Beijing Isolates from Astana, Kazakhstan.

Kazakhstan ranks among the countries with the highest number of MDR-TB patients per 100,000 population worldwide. The successful transmission of local MDR strains of Mycobacterium tuberculosis (Mtb) poses a significant threat to disease control. In this study, we employed whole-genome sequencing to examine drug resistance, compensatory mutations, population structure, and transmission patterns in a sample of 24 clinical isolates of L2/Beijing Mtb collected in Astana, Kazakhstan between 2021 and 2022. The genotypic prediction of Mtb susceptibility to anti-TB agents was consistent with the phenotypic susceptibility, except for bedaquiline. An analysis of resistance-associated genes characterized most of the isolates as pre-extensively drug-resistant tuberculosis (pre-XDR-TB) (n = 15; 62.5%). The phylogenetic analysis grouped the isolates into four transmission clusters; the dominant cluster was assigned to the "aggressive" Central Asia outbreak (CAO) clade of L2/Beijing (n = 15; 62.5%). Thirteen mutations with putative compensatory effects were observed exclusively in Mtb isolates containing the rpoB S450L mutation. The putative compensatory mutations had a stabilizing effect on RpoABC protein stability and dynamics. The high prevalence of the CAO clade in the population structure of Mtb may explain the rapid spread of MDR-TB in Kazakhstan.

SEQUENCE-FREE PHYLOGENETICS WITH MASS SPECTROMETRY.

An alternative, more rapid, sequence-free approach to build phylogenetic trees has been conceived and implemented. Molecular phylogenetics has continued to mostly focus on improvement in tree construction based on gene sequence alignments. Here protein-based phylogenies are constructed using numerical data sets ("phylonumerics") representing the masses of peptide segments recorded in a mass mapping experiment. This truly sequence-free method requires no gene sequences, nor their alignment, to build the trees affording a considerable time and cost-saving to conventional phylogenetics methods. The approach also calculates single point amino acid mutations from a comparison of mass pairs from different maps in the data set and displays these at branch nodes across the tree together with their frequency. Studies of the consecutive, and near-consecutive, ancestral and descendant mutations across interconnected branches of a mass tree allow putative adaptive, epistatic, and compensatory mutations to be identified in order to investigate mechanisms associated with evolutionary processes and pathways. A side-by-side comparison of this sequence-free approach and conventional gene sequence phylogenetics is discussed.

Fitness Cost and Compensatory Evolution in Levofloxacin-Resistant Mycobacterium aurum.

We isolated spontaneous levofloxacin-resistant strains of Mycobacterium aurum to study the fitness cost and compensatory evolution of fluoroquinolone resistance in mycobacteria. Five of six mutant strains with substantial growth defects showed restored fitness after being serially passaged for 18 growth cycles, along with increased cellular ATP level. Whole-genome sequencing identified putative compensatory mutations in the glgC gene that restored the fitness of the resistant strains, presumably by altering the bacterial energy metabolism.

Determination of potentially novel compensatory mutations in rpoc associated with rifampin resistance and rpob mutations in Mycobacterium tuberculosis Clinical isolates from peru.

Rifampicin (RIF) resistance in Mycobacterium tuberculosis is frequently caused by mutations in the rpoB gene. These mutations are associated with a fitness cost, which can be overcome by compensatory mutations in other genes, among which rpoC may be the most important. We analyzed 469 Peruvian M. tuberculosis clinical isolates to identify compensatory mutations in rpoC/rpoA associated with RIF resistance. The M. tuberculosis isolates were collected and tested for RIF susceptibility and spoligotyping. Samples were sequenced and aligned to the reference genome to identify mutations. By analyzing the sequences and the metadata, we identified a list of rpoC mutations exclusively associated with RIF resistance and mutations in rpoB. We then evaluated the distribution of these mutations along the protein sequence and tridimensional structure. One hundred and twenty-five strains were RIF susceptible and 346 were resistant. We identified 35 potential new compensatory mutations, some of which were distributed on the interface surface between rpoB and rpoC, arising in clusters and suggesting the presence of hotspots for compensatory mutations. This study identifies 35 putative novel compensatory mutations in the β' subunit of M. tuberculosis RNApol. Six of these (S428T, L507V, A734V, I997V, and V1252LM) are considered most likely to have a compensatory role, as they fall in the interaction zone of the two subunits and the mutation did not lead to any change in the protein's physical-chemical properties.

Complex Transmission Patterns and Age-Related Dynamics of a Selfish mtDNA Deletion.

Despite wide-ranging implications of selfish mitochondrial DNA (mtDNA) elements for human disease and topics in evolutionary biology (e.g., speciation), the forces controlling their formation, age-related accumulation, and offspring transmission remain largely unknown. Selfish mtDNA poses a significant challenge to genome integrity, mitochondrial function, and organismal fitness. For instance, numerous human diseases are associated with mtDNA mutations; however, few genetic systems can simultaneously represent pathogenic mitochondrial genome evolution and inheritance. The nematode Caenorhabditis briggsae is one such system. Natural C. briggsae isolates harbor varying levels of a large-scale deletion affecting the mitochondrial nduo-5 gene, termed nad5Δ. A subset of these isolates contains putative compensatory mutations that may reduce the risk of deletion formation. We studied the dynamics of nad5Δ heteroplasmy levels during animal development and transmission from mothers to offspring in genetically diverse C. briggsae natural isolates. Results support previous work demonstrating that nad5Δ is a selfish element and that heteroplasmy levels of this deletion can be quite plastic, exhibiting high degrees of inter-family variability and divergence between generations. The latter is consistent with a mitochondrial bottleneck effect, and contrasts with previous findings from a laboratory-derived model uaDf5 mtDNA deletion in C. elegans. However, we also found evidence for among-isolate differences in the ability to limit nad5Δ accumulation, the pattern of which suggested that forces other than the compensatory mutations are important in protecting individuals and populations from rampant mtDNA deletion expansion over short time scales.

Comprehensive identification of compensatory mutations in rifampicin-resistant Mycobacterium tuberculosis strains

Objective: To comprehensively identify compensatory mutations in rpoA, rpoB and rpoC genes of rifampicin-resistant Mycobacterium tuberculosis (RIF-r MTB) and to evaluate the effect of rifampicin-resistant mutation type and lineage background on occurrence of compensatory mutations. Methods: Published MTB whole genome sequencing data were searched and downloaded. RIF-r MTB was identified through known rifampicin-resistant mutations. Based on parallel evolutionary patterns, we identified putative compensatory mutations in the phylogenetic tree and calculated proportions of accumulating compensatory mutations in each rifampicin-resistant mutations' type and lineage background of RIF-r MTB. Statistic significance was analyzed by chi-square test. Results: A total of 8 453 global MTB whole genome sequencing data were downloaded form ENA (covering 12 countries), including 1 749 RIF-r MTB. Based on phylogenetic analysis, we totally identified 60 putative compensatory mutations (6 in rpoA gene, 16 in rpoB gene and 38 in rpoC gene), 11 of which were newly reported. RIF-R strains carrying rpoB S450L (41.7%, 279/669) had a significant higher chance to accumulate compensatory mutations than strains with other rpoB mutations (8.0%, 31/388, χ(2)=378.5, P<0.000 1). In addition, RIF-R strains from lineage 2 (34.0%, 223/656) had a significant higher chance to accumulate compensatory mutations than strains from other lineages [lineage1: 4.7%(2/43), 2/43, lineage3: 12.5%(4/32), 4/32, lineage4: 15.1%(78/517), 78/517; χ(2)=238.5, P<0.000 1]. Conclusions: Our study comprehensively identified putative rifampicin-resistant compensatory mutations of rifampicin resistance. RIF-R strains carrying rpoB S450L mutation or from lineage 2 had a significantly higher chance to accumulate compensatory mutations than strains either with other rpoB mutations or from other lineages.

Whole-Genome Sequencing Reveals the Origin and Rapid Evolution of an Emerging Outbreak Strain of Streptococcus pneumoniae 12F.

Streptococcus pneumoniaeis a major cause of community-acquired pneumonia and septicemia in adults. The global drug-susceptible capsular serotype 12F, clonal complex 218 caused several outbreaks in the United States between 1989 and 2008, as well as a recent large outbreak in Manitoba, Canada, that resulted in 36 cases of septicemia and 3 deaths. The evolutionary origin of the Canadian outbreak strain and its relationship to the historical US outbreak strains are not known. Whole-genome deep sequencing was performed on isolates from the Canadian outbreak (n = 36), the US outbreaks (n = 9), and nonoutbreak surveys (n = 21). Phylogenomic analysis and comparative genomics were used to assess evolutionary relationships and to detect gene content differences between the isolates. The Canadian outbreak was closely related to sporadic cases that occurred preoutbreak in cross-border geographic regions in Manitoba, North Dakota, and Iowa. The emerging Canadian strain differed from US strains by acquisition of a cell-surface protein and macrolide resistance determinants via incorporation of a 5.3-kb mega cassette harboringmsrDandmefE Furthermore, during 11 months of transmission, this clone evolved rapidly and acquired fluoroquinolone resistance through precise stepwise mutations in bothparCandgyrA, and putative compensatory mutations inuraAorIMPDHunder drug selection. Alarmingly, this drug-resistant clone appears to have spread quickly to other regions of Canada and the United States, and replaced drug-susceptible strains. Whole-genome sequencing revealed an independent emergence and secondary adaptation of a new virulent and drug-resistant pneumococcal epidemic clone. Ongoing molecular surveillance is required, and measures to prevent its spread should be developed.

Primer ID ultra-deep sequencing reveals dynamics of drug resistance-associated variants in breakthrough hepatitis C viruses: relevance to treatment outcome and resistance screening.

Use of direct-acting antiviral drugs (DAAs) that target HCV may be hampered by the rapid selection of viral strains that harbour drug resistance-associated variants (RAVs). These RAVs are often associated with a fitness cost and tend to occur on low-frequency strains within treatment-naive subjects. To address the clinical relevance of low frequency RAVs in the setting of DAAs, this study utilized a Primer ID ultra-deep sequencing approach to mitigate PCR errors and bias to accurately quantify viral sequences in subjects that failed DAA treatment. Subjects were enrolled in the follow-up study P05063, and had previous treatment with boceprevir and all had detectable RAVs at virological failure (VF) based on Sanger-based population sequencing. Twelve subjects had three time points available: baseline, VF and follow-up (median 830.5 days). Viral RNA was amplified using unique primer identifiers (Primer IDs) and sequenced using 454 ultra-deep sequencing. The sequencing strategy used in this study improved the detection of clinically relevant low frequency strains bearing RAVs compared to population sequencing and showed that these strains can persist for up to 2 years post-treatment failure. Strains carrying multiple RAVs were common in breakthrough viruses. Putative compensatory mutations were identified. The Primer ID ultra-deep sequencing approach identifies RAVs that can reduce drug sensitivity at levels below the detection threshold for population sequencing. The approach also removes PCR errors and biases, suggesting this sequencing strategy should become the standard approach by which to perform temporal quasispecies studies and resistance screening. ClinicalTrials.gov NCT00689390.

Diffuse Large B-Cell Lymphoma Lacking Surface B-Cell Receptor Expression Are Germinal Center B-Cell (GCB)-Type and Lack Evidence of Known Downstream Activating Mutations

Background: Diffuse large B-cell lymphoma (DLBCL) typically express monotypic surface immunoglobulin (sIg) light chain as part of the B-cell receptor (BCR) complex. Occasional cases are, however, sIg-negative and the clinical significance of this remains unclear. Distinct subtypes of DLBCL, including germinal center B-cell-like (GCB) and activated B-cell (ABC) type, exhibit differential dependence upon BCR-mediated survival signals. The role of activating mutations downstream of the BCR signaling complex (i.e. CD79A/B or CARD11), or in parallel costimulatory pathways (i.e. MYD88 mutations), are now well recognized in promoting such survival mechanisms. Surprisingly, however, less is known about the importance of sIg expression, or lack thereof, on various subsets of DLBCL or other aggressive large B cell lymphomas.Design: We sought to further characterize the clinicopathologic and genetic features of primary human surface Ig-negative DLBCL-not otherwise specified (NOS). All cases of large B-cell lymphoma that had flow cytometry performed at the University of Rochester from 2010-2013 were screened for evidence of a surface Ig-negative B-cell population. Additional available immunophenotypic data and/or clinicopathologic features were recorded. PCR heteroduplex analysis and sequencing of fresh frozen tumor samples was performed to identify putative compensatory activating mutations in MYD88 and CD79A /B, respectively. In addition, using multiplexed PCR and sequencing analysis we attempted to identify the presence or absence of functional immunoglobulin heavy chain (IGH)gene rearrangements.Results: Among 209 diagnostic specimens analyzed, 13 cases (6.2%) were found to contain a discrete population of CD19/CD20-positive B-cells lacking sIg light chain expression. These included 8 DLBCL-NOS, nearly all of which were GCB-type, including 6 out of 7 cases, or 83%, which were CD10-positive. This proportion is significantly greater than the 42% of all DLBCL-NOS that express CD10 (based on data from a large publicly available on-line database; p < 0.025, based on binomial probability). Furthermore, of the other high grade sIg-negative B cell lymphomas identified, 4 out of 5 were also CD10-positive. These included 2 cases of B-cell lymphoma, unclassifiable (BCL-U), intermediate between DLBCL and Burkitt lymphoma, and 2 cases of transformed follicular lymphoma (FL) - therefore also histologically and immunophenotypically consistent with derivation from a germinal center B-cell. The only exception was a non-GCB primary CNS lymphoma. Among 5 cases of sIg- DLBCL available for analysis, all were negative for the MYD88 L265P activating mutation, while only a single case contained an activating mutation in CD79A/B (CD79B Y196 ITAM mutation). In the same set, only 2 (of 5) cases produced amplicons by IGH PCR, and sequencing revealed only one of these had an intact reading frame.Conclusion: DLBCL-NOS lacking definitive sIg light chain expression are infrequent (less than 10%), but these tumors are predominantly GCB subtype. This observation lends support to the model that these tumors rely less on BCR-mediated survival signals, which may have implications for choosing targeted therapy. The mechanism by which these lymphoma survive despite the apparent lack of functional BCR is unknown, but compensatory activating mutations known to be frequent among ABC-type DLBCL appear to be lacking. DisclosuresNo relevant conflicts of interest to declare.

Evolution of Extensively Drug-Resistant Tuberculosis over Four Decades: Whole Genome Sequencing and Dating Analysis of Mycobacterium tuberculosis Isolates from KwaZulu-Natal.

BackgroundThe continued advance of antibiotic resistance threatens the treatment and control of many infectious diseases. This is exemplified by the largest global outbreak of extensively drug-resistant (XDR) tuberculosis (TB) identified in Tugela Ferry, KwaZulu-Natal, South Africa, in 2005 that continues today. It is unclear whether the emergence of XDR-TB in KwaZulu-Natal was due to recent inadequacies in TB control in conjunction with HIV or other factors. Understanding the origins of drug resistance in this fatal outbreak of XDR will inform the control and prevention of drug-resistant TB in other settings. In this study, we used whole genome sequencing and dating analysis to determine if XDR-TB had emerged recently or had ancient antecedents.Methods and FindingsWe performed whole genome sequencing and drug susceptibility testing on 337 clinical isolates of Mycobacterium tuberculosis collected in KwaZulu-Natal from 2008 to 2013, in addition to three historical isolates, collected from patients in the same province and including an isolate from the 2005 Tugela Ferry XDR outbreak, a multidrug-resistant (MDR) isolate from 1994, and a pansusceptible isolate from 1995. We utilized an array of whole genome comparative techniques to assess the relatedness among strains, to establish the order of acquisition of drug resistance mutations, including the timing of acquisitions leading to XDR-TB in the LAM4 spoligotype, and to calculate the number of independent evolutionary emergences of MDR and XDR. Our sequencing and analysis revealed a 50-member clone of XDR M. tuberculosis that was highly related to the Tugela Ferry XDR outbreak strain. We estimated that mutations conferring isoniazid and streptomycin resistance in this clone were acquired 50 y prior to the Tugela Ferry outbreak (katG S315T [isoniazid]; gidB 130 bp deletion [streptomycin]; 1957 [95% highest posterior density (HPD): 1937–1971]), with the subsequent emergence of MDR and XDR occurring 20 y (rpoB L452P [rifampicin]; pncA 1 bp insertion [pyrazinamide]; 1984 [95% HPD: 1974–1992]) and 10 y (rpoB D435G [rifampicin]; rrs 1400 [kanamycin]; gyrA A90V [ofloxacin]; 1995 [95% HPD: 1988–1999]) prior to the outbreak, respectively. We observed frequent de novo evolution of MDR and XDR, with 56 and nine independent evolutionary events, respectively. Isoniazid resistance evolved before rifampicin resistance 46 times, whereas rifampicin resistance evolved prior to isoniazid only twice. We identified additional putative compensatory mutations to rifampicin in this dataset. One major limitation of this study is that the conclusions with respect to ordering and timing of acquisition of mutations may not represent universal patterns of drug resistance emergence in other areas of the globe.ConclusionsIn the first whole genome-based analysis of the emergence of drug resistance among clinical isolates of M. tuberculosis, we show that the ancestral precursor of the LAM4 XDR outbreak strain in Tugela Ferry gained mutations to first-line drugs at the beginning of the antibiotic era. Subsequent accumulation of stepwise resistance mutations, occurring over decades and prior to the explosion of HIV in this region, yielded MDR and XDR, permitting the emergence of compensatory mutations. Our results suggest that drug-resistant strains circulating today reflect not only vulnerabilities of current TB control efforts but also those that date back 50 y. In drug-resistant TB, isoniazid resistance was overwhelmingly the initial resistance mutation to be acquired, which would not be detected by current rapid molecular diagnostics employed in South Africa that assess only rifampicin resistance.

Genomic epidemiology of multidrug-resistant Mycobacterium tuberculosis during transcontinental spread.

The transcontinental spread of multidrug-resistant (MDR) tuberculosis is poorly characterized in molecular epidemiologic studies. We used genomic sequencing to understand the establishment and dispersion of MDR Mycobacterium tuberculosis within a group of immigrants to the United States. We used a genomic epidemiology approach to study a genotypically matched (by spoligotype, IS6110 restriction fragment length polymorphism, and mycobacterial interspersed repetitive units-variable number of tandem repeat signature) lineage 2/Beijing MDR strain implicated in an outbreak of tuberculosis among refugees in Thailand and consecutive cases within California. All 46 MDR M. tuberculosis genomes from both Thailand and California were highly related, with a median difference of 10 single-nucleotide polymorphisms (SNPs). The Wat Tham Krabok (WTK) strain is a new sequence type distinguished from all known Beijing strains by 55 SNPs and a genomic deletion (Rv1267c) associated with increased fitness. Sequence data revealed a highly prevalent MDR strain that included several closely related but distinct allelic variants within Thailand, rather than the occurrence of a single outbreak. In California, sequencing data supported multiple independent introductions of WTK with subsequent transmission and reactivation within the state, as well as a potential super spreader with a prolonged infectious period. Twenty-seven drug resistance-conferring mutations and 4 putative compensatory mutations were found within WTK strains. Genomic sequencing has substantial epidemiologic value in both low- and high-burden settings in understanding transmission chains of highly prevalent MDR strains.

Comprehensive phenotypic characterization of rifampicin resistance mutations in Salmonella provides insight into the evolution of resistance in Mycobacterium tuberculosis.

Mutations in the β-subunit of RNA polymerase (RNAP), encoded by rpoB, are responsible for rifampicin resistance (Rif(R)). Although many mutations in rpoB can reduce susceptibility, only a few are frequent amongst Rif(R) clinical Mycobacterium tuberculosis (MTB) isolates. It has been suggested that there is a negative correlation between the fitness costs of Rif(R) mutations and their respective clinical frequency, but so far comparable fitness cost measurements have only been conducted for a very limited number of Rif(R) mutations. We tested this hypothesis using Salmonella and Mycobacterium smegmatis as model organisms. We constructed 122 different Rif(R) mutations in Salmonella. MICs and relative fitness costs in the presence and absence of rifampicin were determined for each mutant, including for a smaller number of Rif(R) M. smegmatis strains. Results were compared with available mutation frequency data from clinical MTB isolates. (i) Rif(R) mutations frequently found in MTB isolates have a fitness cost in Salmonella Typhimurium and M. smegmatis. (ii) Clinically frequent Rif(R) mutations have a high rifampicin MIC. (iii) There is a strong correlation between the magnitude of the fitness cost of a Rif(R) mutation in Salmonella Typhimurium or M. smegmatis and the frequency with which that mutation is associated with secondary (putative compensatory) mutations in RNAP of clinical MTB isolates. This suggests that the success of Rif(R) mutations in clinical MTB isolates may be dependent not only on a low initial fitness cost, but rather the results of three factors: (i) a high rifampicin MIC; (ii) a relatively low initial fitness cost; and (iii) the ability to additionally acquire compensatory mutations selected to further reduce fitness cost.

Genetic characterization of compensatory evolution in strains carrying rpoB Ser531Leu, the rifampicin resistance mutation most frequently found in clinical isolates

The evolution of rifampicin resistance in Mycobacterium tuberculosis is a major threat to effective tuberculosis therapy. Much is known about the initial emergence of rifampicin resistance, but the further evolution of these resistant strains has only lately been subject to investigation. Although resistance can be caused by many different mutations in rpoB, among clinical M. tuberculosis isolates the mutation rpoB S531L is overwhelmingly the most frequently found. Clinical isolates with rpoB S531L frequently carry additional mutations in genes for RNA polymerase subunits, and it has been speculated that these are fitness-compensatory mutations, ameliorating the fitness cost of the primary resistance mutation. We tested this hypothesis using Salmonella as a model organism. We created the rpoB S531L mutation in Salmonella and then evolved independent lineages with selection for mutants with increased relative fitness. Relative fitness associated with putative compensatory mutations was measured after genetic reconstruction in isogenic strains. Compensatory mutations were identified in genes coding for different subunits of RNA polymerase: rpoA, rpoB and rpoC. Genetic reconstructions demonstrated that each of these secondary mutations reduced the fitness cost of the rpoB S531L resistance mutation. The compensatory mutations identified in Salmonella cluster in similar locations to the additional mutations found in M. tuberculosis isolates. These new data strongly support the idea that many of the previously identified rpoA, rpoB and rpoC mutations in rifampicin-resistant M. tuberculosis (rpoB S531L) are indeed fitness-compensatory mutations.

Putative Compensatory Mutations in the rpoC Gene of Rifampin-Resistant Mycobacterium tuberculosis Are Associated with Ongoing Transmission

Rifampin resistance in clinical isolates of Mycobacterium tuberculosis arises primarily through the selection of bacterial variants harboring mutations in the 81-bp rifampin resistance-determining region of the rpoB gene. While these mutations were shown to infer a fitness cost in the absence of antibiotic pressure, compensatory mutations in rpoA and rpoC were identified which restore the fitness of rifampin-resistant bacteria carrying mutations in rpoB. To investigate the epidemiological relevance of these compensatory mutations, we analyzed 286 drug-resistant and 54 drug-susceptible clinical M. tuberculosis isolates from the Western Cape, South Africa, a high-incidence setting of multidrug-resistant tuberculosis. Sequencing of a portion of the RpoA-RpoC interaction region of the rpoC gene revealed that 23.5% of all rifampin-resistant isolates tested carried a nonsynonymous mutation in this region. These putative compensatory mutations in rpoC were associated with transmission, as 30.8% of all rifampin-resistant isolates with an IS6110 restriction fragment length polymorphism (RFLP) pattern belonging to a recognized RFLP cluster harbored putative rpoC mutations. Such mutations were present in only 9.4% of rifampin-resistant isolates with unique RFLP patterns (P < 0.01). Moreover, these putative compensatory mutations were associated with specific strain genotypes and the rpoB S531L rifampin resistance mutation. Among isolates harboring this rpoB mutation, 44.1% also harbored rpoC mutations, while only 4.1% of the isolates with other rpoB mutations exhibited mutations in rpoC (P < 0.001). Our study supports a role for rpoC mutations in the transmission of multidrug-resistant tuberculosis and illustrates how epistatic interactions between drug resistance-conferring mutations, compensatory mutations, and different strain genetic backgrounds might influence compensatory evolution in drug-resistant M. tuberculosis.

Temporal Association of HLA-B*81:01- and HLA-B*39:10-Mediated HIV-1 p24 Sequence Evolution with Disease Progression

HLA-B*81:01 and HLA-B*39:10 alleles have been associated with viremic control in HIV-1 subtype C infection. Both alleles restrict the TL9 epitope in p24 Gag, and cytotoxic-T-lymphocyte (CTL)-mediated escape mutations in this epitope have been associated with an in vitro fitness cost to the virus. We investigated the timing and impact of mutations in the TL9 epitope on disease progression in five B*81:01- and two B*39:10-positive subtype C-infected individuals. Whereas both B*39:10 participants sampled at 2 months postinfection had viruses with mutations in the TL9 epitope, in three of the five (3/5) B*81:01 participants, TL9 escape mutations were only detected 10 months after infection, taking an additional 10 to 15 months to reach fixation. In the two remaining B*81:01 individuals, one carried a TL9 escape variant at 2 weeks postinfection, whereas no escape mutations were detected in the virus from the other participant for up to 33 months postinfection, despite CTL targeting of the epitope. In all participants, escape mutations in TL9 were linked to coevolving residues in the region of Gag known to be associated with host tropism. Late escape in TL9, together with coevolution of putative compensatory mutations, coincided with a spontaneous increase in viral loads in two individuals who were otherwise controlling the infection. These results provide in vivo evidence of the detrimental impact of B*81:01-mediated viral evolution, in a single Gag p24 epitope, on the control of viremia.

Impact of HLA-B*81-Associated Mutations in HIV-1 Gag on Viral Replication Capacity

HIV-1 attenuation resulting from immune escape mutations selected in Gag may contribute to slower disease progression in HIV-1-infected individuals expressing certain HLA class I alleles. We previously showed that the protective allele HLA-B*81 and the HLA-B*81-selected Gag T186S mutation are strongly associated with a lower viral replication capacity of recombinant viruses encoding Gag-protease derived from individuals chronically infected with HIV-1 subtype C. In the present study, we directly tested the effect of this mutation on viral replication capacity. In addition, we investigated potential compensatory effects of various polymorphisms, including other HLA-B*81-associated mutations that significantly covary with the T186S mutation. Mutations were introduced into a reference subtype B backbone and into patient-derived subtype C sequences in subtype B and C backbones by site-directed mutagenesis. The exponential-phase growth of mutant and wild-type viruses was assayed by flow cytometry of a green fluorescent protein reporter T cell line or by measurement of HIV-1 reverse transcriptase activity in culture supernatants. Engineering of the T186S mutation alone into all patient-derived subtype C sequences failed to yield replication-competent viruses, while in the subtype B sequence, the T186S mutation resulted in impaired replication capacity. Only the T186S mutation in combination with the T190I mutation yielded replication-competent viruses for all virus backbones tested; however, these constructs replicated slower than the wild type, suggesting that only partial compensation is mediated by the T190I mutation. Constructs encoding the T186S mutation in combination with other putative compensatory mutations were attenuated or defective. These results suggest that the T186S mutation is deleterious to HIV-1 subtype C replication and likely requires complex compensatory pathways, which may contribute to the clinical benefit associated with HLA-B*81.

Muller's Ratchet and compensatory mutation in Caenorhabditis briggsae mitochondrial genome evolution

BackgroundThe theory of Muller' Ratchet predicts that small asexual populations are doomed to accumulate ever-increasing deleterious mutation loads as a consequence of the magnified power of genetic drift and mutation that accompanies small population size. Evidence for Muller's Ratchet and knowledge on its underlying molecular mechanisms, however, are lacking for natural populations.ResultsWe characterized mitochondrial genome evolutionary processes in Caenorhabditis briggsae natural isolates to show that numerous lineages experience a high incidence of nonsynonymous substitutions in protein-coding genes and accumulate unusual deleterious noncoding DNA stretches with associated heteroplasmic function-disrupting genome deletions. Isolate-specific deletion proportions correlated negatively with nematode fecundity, suggesting that these deletions might negatively affect C. briggsae fitness. However, putative compensatory mutations were also observed that are predicted to reduce heteroplasmy levels of deleterious deletions. Paradoxically, compensatory mutations were observed in one major intraspecific C. briggsae clade where population sizes are estimated to be very small (and selection is predicted to be relatively weak), but not in a second major clade where population size estimates are much larger and selection is expected to be more efficient.ConclusionThis study provides evidence that the mitochondrial genomes of animals evolving in nature are susceptible to Muller's Ratchet, suggests that context-dependent compensatory mutations can accumulate in small populations, and predicts that Muller's Ratchet can affect fundamental evolutionary forces such as the rate of mutation.