In vitro analysis of rates and spectra of mutations in a polymorphic region of the Rv0746 PE_PGRS gene of Mycobacterium tuberculosis

PE_PGRS38 Interaction With HAUSP Downregulates Antimycobacterial Host Defense via TRAF6

Mycobacterium tuberculosis (Mtb) is the causative pathogen of tuberculosis (TB), which manipulates the host immunity to ensure survival and colonization in the host. Mtb possess a unique family of proteins, named PE_PGRS, associated with Mtb pathogenesis. Thus, elucidation of the functions of PE_PGRS proteins is necessary to understand TB pathogenesis. Here, we investigated the role of PE_PGRS38 binding to herpesvirus-associated ubiquitin-specific protease (HAUSP, USP7) in regulating the activity of various substrate proteins by modulating their state of ubiquitination. We constructed the recombinant PE_PGRS38 expressed in M. smegmatis (Ms_PE_PGRS38) to investigate the role of PE_PGRS38. We found that Ms_PE_PGRS38 regulated the cytokine levels in murine bone marrow-derived macrophages by inhibiting the deubiquitination of tumor necrosis factor receptor-associated factor (TRAF) 6 by HAUSP. Furthermore, the PE domain in PE_PGRS38 was identified as essential for mediating TRAF6 deubiquitination. Ms_PE_PGRS38 increased the intracellular burden of bacteria by manipulating cytokine levels in vitro and in vivo. Overall, we revealed that the interplay between HAUSP and PE_PGRS38 regulated the inflammatory response to increase the survival of mycobacteria.

DNA-Dependent Binding of Nargenicin to DnaE1 Inhibits Replication in Mycobacterium tuberculosis

Natural products provide a rich source of potential antimicrobials for treating infectious diseases for which drug resistance has emerged. Foremost among these diseases is tuberculosis. Assessment of the antimycobacterial activity of nargenicin, a natural product that targets the replicative DNA polymerase of Staphylococcus aureus, revealed that it is a bactericidal genotoxin that induces a DNA damage response in Mycobacterium tuberculosis (Mtb) and inhibits growth by blocking the replicative DNA polymerase, DnaE1. Cryo-electron microscopy revealed that binding of nargenicin to Mtb DnaE1 requires the DNA substrate such that nargenicin is wedged between the terminal base pair and the polymerase and occupies the position of both the incoming nucleotide and templating base. Comparative analysis across three bacterial species suggests that the activity of nargenicin is partly attributable to the DNA binding affinity of the replicative polymerase. This work has laid the foundation for target-led drug discovery efforts focused on Mtb DnaE1.

Nucleotide Metabolism and DNA Replication

The development and application of a highly versatile suite of tools for mycobacterial genetics, coupled with widespread use of "omics" approaches to elucidate the structure, function, and regulation of mycobacterial proteins, has led to spectacular advances in our understanding of the metabolism and physiology of mycobacteria. In this article, we provide an update on nucleotide metabolism and DNA replication in mycobacteria, highlighting key findings from the past 10 to 15 years. In the first section, we focus on nucleotide metabolism, ranging from the biosynthesis, salvage, and interconversion of purine and pyrimidine ribonucleotides to the formation of deoxyribonucleotides. The second part of the article is devoted to DNA replication, with a focus on replication initiation and elongation, as well as DNA unwinding. We provide an overview of replication fidelity and mutation rates in mycobacteria and summarize evidence suggesting that DNA replication occurs during states of low metabolic activity, and conclude by suggesting directions for future research to address key outstanding questions. Although this article focuses primarily on observations from Mycobacterium tuberculosis, it is interspersed, where appropriate, with insights from, and comparisons with, other mycobacterial species as well as better characterized bacterial models such as Escherichia coli. Finally, a common theme underlying almost all studies of mycobacterial metabolism is the potential to identify and validate functions or pathways that can be exploited for tuberculosis drug discovery. In this context, we have specifically highlighted those processes in mycobacterial DNA replication that might satisfy this critical requirement.

A combinatorial role for MutY and Fpg DNA glycosylases in mutation avoidance in Mycobacterium smegmatis

Hydroxyl radical (OH) among reactive oxygen species cause damage to nucleobases with thymine being the most susceptible, whilst in contrast, the singlet oxygen ((1)02) targets only guanine bases. The high GC content of mycobacterial genomes predisposes these organisms to oxidative damage of guanine. The exposure of cellular DNA to OH and one-electron oxidants results in the formation of two main degradation products, the pro-mutagenic 8-oxo-7,8-dihydroguanine (8-oxoGua) and the cytotoxic 2,6-diamino-4-hydroxy-5-formamidopyrimidine (FapyGua). These lesions are repaired through the base excision repair (BER) pathway and we previously, demonstrated a combinatorial role for the mycobacterial Endonuclease III (Nth) and the Nei family of DNA glycosylases in mutagenesis. In addition, the formamidopyrimidine (Fpg/MutM) and MutY DNA glycosylases have also been implicated in mutation avoidance and BER in mycobacteria. In this study, we further investigate the combined role of MutY and the Fpg/Nei DNA glycosylases in Mycobacterium smegmatis and demonstrate that deletion of mutY resulted in enhanced sensitivity to oxidative stress, an effect which was not exacerbated in Δfpg1 Δfpg2 or Δnei1 Δnei2 double mutant backgrounds. However, combinatorial loss of the mutY, fpg1 and fpg2 genes resulted in a significant increase in mutation rates suggesting interplay between these enzymes. Consistent with this, there was a significant increase in C → A mutations with a corresponding change in cell morphology of rifampicin resistant mutants in the Δfpg1 Δfpg2 ΔmutY deletion mutant. In contrast, deletion of mutY together with the nei homologues did not result in any growth/survival defects or changes in mutation rates. Taken together these data indicate that the mycobacterial mutY, in combination with the Fpg DNA N-glycosylases, plays an important role in controlling mutagenesis under oxidative stress.

An overview to understand the role of PE_PGRS family proteins in Mycobacterium tuberculosis H37 Rv and their potential as new drug targets

Tuberculosis has long been the scourge of humanity, claiming millions of lives. The family of PE_PGRS gene has been attributed to the Mycobacterium tuberculosis pathogenesis over the past few decades. The gene of PE_PGRS family proteins are most often clustered in a region of the genome often as overlapping genes and role in cell surface markers, adhesion and invasion of defense cells of the host (macrophage and dendritic cells). The proline-glutamic acid (PE) domain is responsible for the cellular localization of these proteins on bacterial cells. This gene family shows immense genetic variability in terms of multiple insertion-deletions and single-nucleotide polymorphisms as seen in PE_PGRS9, PE_PGRS17, PE_PGRS18, and PE_PGRS33. In spite of variability, there are indications of shared epitopes in these proteins. Few of these gene sequences that have been studied from evolutionary perspective show indication of positive selection and also landmarks of recent evolutionary events. Many of these proteins show calcium-binding motifs and consequently seen to be responsible in inhibition of phagolysosome formation via a calmodulin-kinase-dependent pathway. A number of PE_PGRS genes were tested for its expression with different growth conditions in vitro and in vivo, among which the contrast in expressivity was seen vividly in PE_PGRS16 (upregulated) and PE_PGRS26 (downregulated) in bacteria persisting in macrophages. Similarly, PE_PGRS33 has been indicated in macrophagial necrosis by a tumor necrosis factor-α-induced pathway. These PE_PGRS family genes may be an interesting subject for research and development. Their fibronectin-binding and calcium-binding property may be strongly implicated in immunopathogenesis of virulent M. tuberculosis strain. In this review, an attempt has been made to evaluate and present data for better understanding of in vivo pathogen functions, for understanding the physiological significance of PE_PGRS gene family, and their potential as new drug targets.

The role of glutamine oxoglutarate aminotransferase and glutamate dehydrogenase in nitrogen metabolism in Mycobacterium bovis BCG

Recent evidence suggests that the regulation of intracellular glutamate levels could play an important role in the ability of pathogenic slow-growing mycobacteria to grow in vivo. However, little is known about the in vitro requirement for the enzymes which catalyse glutamate production and degradation in the slow-growing mycobacteria, namely; glutamine oxoglutarate aminotransferase (GOGAT) and glutamate dehydrogenase (GDH), respectively. We report that allelic replacement of the Mycobacterium bovis BCG gltBD-operon encoding for the large (gltB) and small (gltD) subunits of GOGAT with a hygromycin resistance cassette resulted in glutamate auxotrophy and that deletion of the GDH encoding-gene (gdh) led to a marked growth deficiency in the presence of L-glutamate as a sole nitrogen source as well as reduction in growth when cultured in an excess of L-asparagine.

The contribution of Nth and Nei DNA glycosylases to mutagenesis in Mycobacterium smegmatis

The increased prevalence of drug resistant strains of Mycobacterium tuberculosis (Mtb) indicates that significant mutagenesis occurs during tuberculosis disease in humans. DNA damage by host-derived reactive oxygen/nitrogen species is hypothesized to be critical for the mutagenic process in Mtb thus, highlighting an important role for DNA repair enzymes in maintenance of genome fidelity. Formamidopyrimidine (Fpg/MutM/Fapy) and EndonucleaseVIII (Nei) constitute the Fpg/Nei family of DNA glycosylases and together with EndonucleaseIII (Nth) are central to the base excision repair pathway in bacteria. In this study we assess the contribution of Nei and Nth DNA repair enzymes in Mycobacterium smegmatis (Msm), which retains a single nth homologue and duplications of the Fpg (fpg1 and fpg2) and Nei (nei1 and nei2) homologues. Using an Escherichia coli nth deletion mutant, we confirm the functionality of the mycobacterial nth gene in the base excision repair pathway. Msm mutants lacking nei1, nei2 and nth individually or in combination did not display aberrant growth in broth culture. Deletion of nth individually results in increased UV-induced mutagenesis and combinatorial deletion with the nei homologues results in reduced survival under oxidative stress conditions and an increase in spontaneous mutagenesis to rifampicin. Deletion of nth together with the fpg homolgues did not result in any growth/survival defects or changes in mutation rate. Furthermore, no differential emergence of the common rifampicin resistance conferring genotypes were noted. Collectively, these data confirm a role for Nth in base excision repair in mycobacteria and further highlight a novel interplay between the Nth and Nei homologues in spontaneous mutagenesis.

Evolution of smooth tubercle Bacilli PE and PE_PGRS genes: evidence for a prominent role of recombination and imprint of positive selection

BACKGROUND

PE and PE_PGRS are two mycobateria-restricted multigene families encoding membrane associated and secreted proteins that have expanded mainly in the pathogenic species, notably the Mycobacterium tuberculosis complex (MTBC). Several lines of evidence attribute to PE and PE_PGRS genes critical roles in mycobacterial pathogenicity. To get more insight into the nature of these genes, we sought to address their evolutionary trajectories in the group of smooth tubercle bacilli (STB), the putative ancestor of the clonal MTBC.

METHODOLOGY/PRINCIPAL FINDINGS

By focussing on six polymorphic STB PE/PE_PGRS genes, we demonstrate significant incongruence among single gene genealogies and detect strong signals of recombination using various approaches. Coalescent-based estimation of population recombination and mutation rates (ρ and θ, respectively) indicates that the two mechanisms are of roughly equal importance in generating diversity (ρ/θ = 1.457), a finding in a marked contrast to house keeping genes (HKG) whose evolution is chiefly brought about by mutation (ρ/θ = 0.012). In comparison to HKG, we found 15 times higher mean rate of nonsynonymous substitutions, with strong evidence of positive selection acting on PE_PGRS62 (dN/dS = 1.42), a gene that has previously been shown to be essential for mycobacterial survival in macrophages and granulomas. Imprint of positive selection operating on specific amino acid residues or along branches of PE_PGRS62 phylogenetic tree was further demonstrated using maximum likelihood- and covarion-based approaches, respectively. Strikingly, PE_PGR62 proved highly conserved in present-day MTBC strains.

CONCLUSIONS/SIGNIFICANCE

Overall the data indicate that, in STB, PE/PE_PGRS genes have undergone a strong diversification process that is speeded up by recombination, with evidence of positive selection. The finding that positive selection involved an essential PE_PGRS gene whose sequence appears to be driven to fixation in present-day MTBC strains lends further support to the critical role of PE/PE_PGRS genes in the evolution of mycobacterial pathogenicity.

DNA metabolism in mycobacterial pathogenesis

Fundamental aspects of the lifestyle of Mycobacterium tuberculosis implicate DNA metabolism in bacillary survival and adaptive evolution. The environments encountered by M. tuberculosis during successive cycles of infection and transmission are genotoxic. Moreover, as an obligate pathogen, M. tuberculosis has the ability to persist for extended periods in a subclinical state, suggesting that active DNA repair is critical to maintain genome integrity and bacterial viability during prolonged infection. In this chapter, we provide an overview of the major DNA metabolic pathways identified in M. tuberculosis, and situate key recent findings within the context of mycobacterial pathogenesis. Unlike many other bacterial pathogens, M. tuberculosis is genetically secluded, and appears to rely solely on chromosomal mutagenesis to drive its microevolution within the human host. In turn, this implies that a balance between high versus relaxed fidelity mechanisms of DNA metabolism ensures the maintenance of genome integrity, while accommodating the evolutionary imperative to adapt to hostile and fluctuating environments. The inferred relationship between mycobacterial DNA repair and genome dynamics is considered in the light of emerging data from whole-genome sequencing studies of clinical M. tuberculosis isolates which have revealed the potential for considerable heterogeneity within and between different bacterial and host populations.

A TetR-like regulator broadly affects the expressions of diverse genes in Mycobacterium smegmatis

Transcriptional regulation plays a critical role in the life cycle of Mycobacterium smegmatis and its related species, M. tuberculosis, the causative microbe for tuberculosis. However, the key transcriptional factors involved in broad regulation of diverse genes remain to be characterized in mycobacteria. In the present study, a TetR-like family transcriptional factor, Ms6564, was characterized in M. smegmatis as a master regulator. A conserved 19 bp-palindromic motif was identified for Ms6564 binding using DNaseI footprinting and EMSA. A total of 339 potential target genes for Ms6564 were further characterized by searching the M. smegmatis genome based on the sequence motif. Notably, Ms6564 bound with the promoters of 37 cell cycle and DNA damage/repair genes and regulated positively their expressions. The Ms6564-overexpressed recombinant strain yielded 5-fold lower mutation rates and mutation frequencies, whereas deletion of Ms6564 resulted in ∼5-fold higher mutation rates for the mutant strain compared with the wild-type strain. These findings suggested that Ms6564 may function as a global regulator and might be a sensor necessary for activation of DNA damage/repair genes.

ClpR protein-like regulator specifically recognizes RecA protein-independent promoter motif and broadly regulates expression of DNA damage-inducible genes in mycobacteria

The RecA-dependent DNA damage response pathway (SOS response) appears to be the major DNA repair mechanism in most bacteria, but it has been suggested that a RecA-independent mechanism is responsible for controlling expression of most damage-inducible DNA repair genes in Mycobacterium tuberculosis. The specific reparative responses and molecular mediators involved in the DNA repair mechanism remain largely unclear in this pathogen and its related species. In this study, a mycobacterial ClpR-like regulator, corresponding to Rv2745c in M. tuberculosis and to Ms2694 in M. smegmatis mc(2)155, was found to interact with the promoter regions of multiple damage-inducible DNA repair genes. Specific binding of the ClpR-like factor to the conserved RecA-independent promoter RecA-NDp motif was then confirmed using in vitro electrophoretic mobility shift assays as well as in vivo chromatin immunoprecipitation experiments. The ClpR knock-out experiments, in combination with quantitative real time PCR assays, demonstrated that the expression of these RecA-independent genes were significantly down-regulated in the mutant strain of M. smegmatis in response to a DNA-damaging agent compared with the wild type strain. Furthermore, the ClpR-like factor was shown to contribute to mycobacterial genomic stability. These results enhance our understanding of the function of the ClpR regulator and the regulatory mechanism of RecA-independent DNA repair in mycobacteria.

Spontaneous emergence of multiple drug resistance in tuberculosis before and during therapy

The emergence of drug resistance in M. tuberculosis undermines the efficacy of tuberculosis (TB) treatment in individuals and of TB control programs in populations. Multiple drug resistance is often attributed to sequential functional monotherapy, and standard initial treatment regimens have therefore been designed to include simultaneous use of four different antibiotics. Despite the widespread use of combination therapy, highly resistant M. tb strains have emerged in many settings. Here we use a stochastic birth-death model to estimate the probability of the emergence of multidrug resistance during the growth of a population of initially drug sensitive TB bacilli within an infected host. We find that the probability of the emergence of resistance to the two principal anti-TB drugs before initiation of therapy ranges from 10⁻⁵ to 10⁻⁴; while rare, this is several orders of magnitude higher than previous estimates. This finding suggests that multidrug resistant M. tb may not be an entirely “man-made” phenomenon and may help explain how highly drug resistant forms of TB have independently emerged in many settings.

Role of the DinB homologs Rv1537 and Rv3056 in Mycobacterium tuberculosis

The environment encountered by Mycobacterium tuberculosis during infection is genotoxic. Most bacteria tolerate DNA damage by engaging specialized DNA polymerases that catalyze translesion synthesis (TLS) across sites of damage. M. tuberculosis possesses two putative members of the DinB class of Y-family DNA polymerases, DinB1 (Rv1537) and DinB2 (Rv3056); however, their role in damage tolerance, mutagenesis, and survival is unknown. Here, both dinB1 and dinB2 are shown to be expressed in vitro in a growth phase-dependent manner, with dinB2 levels 12- to 40-fold higher than those of dinB1. Yeast two-hybrid analyses revealed that DinB1, but not DinB2, interacts with the beta-clamp, consistent with its canonical C-terminal beta-binding motif. However, knockout of dinB1, dinB2, or both had no effect on the susceptibility of M. tuberculosis to compounds that form N(2)-dG adducts and alkylating agents. Similarly, deletion of these genes individually or in combination did not affect the rate of spontaneous mutation to rifampin resistance or the spectrum of resistance-conferring rpoB mutations and had no impact on growth or survival in human or mouse macrophages or in mice. Moreover, neither gene conferred a mutator phenotype when expressed ectopically in Mycobacterium smegmatis. The lack of the effect of altering the complements or expression levels of dinB1 and/or dinB2 under conditions predicted to be phenotypically revealing suggests that the DinB homologs from M. tuberculosis do not behave like their counterparts from other organisms.

Evidence for a rapid rate of molecular evolution at the hypervariable and immunogenic Mycobacterium tuberculosis PPE38 gene region

BACKGROUND

PPE38 (Rv2352c) is a member of the large PPE gene family of Mycobacterium tuberculosis and related mycobacteria. The function of PPE proteins is unknown but evidence suggests that many are cell-surface associated and recognised by the host immune system. Previous studies targeting other PPE gene members suggest that some display high levels of polymorphism and it is thought that this might represent a means of providing antigenic variation. We have analysed the genetic variability of the PPE38 genomic region on a cohort of M. tuberculosis clinical isolates representing all of the major phylogenetic lineages, along with the ancestral M. tuberculosis complex (MTBC) member M. canettii, and supplemented this with analysis of publicly available whole genome sequences representing additional M. tuberculosis clinical isolates, other MTBC members and non tuberculous mycobacteria (NTM). Where possible we have extended this analysis to include the adjacent plcABC and PPE39/40 genomic regions.

RESULTS

We show that the ancestral MTBC PPE38 region comprises 2 homologous PPE genes (PPE38 and PPE71), separated by 2 esat-6 (esx)-like genes and that this structure derives from an esx/esx/PPE duplication in the common ancestor of M. tuberculosis and M. marinum. We also demonstrate that this region of the genome is hypervariable due to frequent IS6110 integration, IS6110-associated recombination, and homologous recombination and gene conversion events between PPE38 and PPE71. These mutations result in combinations of gene deletion, gene truncation and gene disruption in the majority of clinical isolates. These mutations were generally found to be IS6110 strain lineage-specific, although examples of additional within-lineage and even within-cluster mutations were observed. Furthermore, we provide evidence that the published M. tuberculosis H37Rv whole genome sequence is inaccurate regarding this region.

CONCLUSION

Our results show that this antigen-encoding region of the M. tuberculosis genome is hypervariable. The observation that numerous different mutations have become fixed within specific lineages demonstrates that this genomic region is undergoing rapid molecular evolution and that further lineage-specific evolutionary expansion and diversification has occurred subsequent to the lineage-defining mutational events. We predict that functional loss of these genes could aid immune evasion. Finally, we also show that the PPE38 region of the published M. tuberculosis H37Rv whole genome sequence is not representative of the ATCC H37Rv reference strain.

Function and regulation of class I ribonucleotide reductase-encoding genes in mycobacteria

Ribonucleotide reductases (RNRs) are crucial to all living cells, since they provide deoxyribonucleotides (dNTPs) for DNA synthesis and repair. In Mycobacterium tuberculosis, a class Ib RNR comprising nrdE- and nrdF2-encoded subunits is essential for growth in vitro. Interestingly, the genome of this obligate human pathogen also contains the nrdF1 (Rv1981c) and nrdB (Rv0233) genes, encoding an alternate class Ib RNR small (R2) subunit and a putative class Ic RNR R2 subunit, respectively. However, the role(s) of these subunits in dNTP provision during M. tuberculosis pathogenesis is unknown. In this study, we demonstrate that nrdF1 and nrdB are dispensable for the growth and survival of M. tuberculosis after exposure to various stresses in vitro and, further, that neither gene is required for growth and survival in mice. These observations argue against a specialist role for the alternate R2 subunits under the conditions tested. Through the construction of nrdR-deficient mutants of M. tuberculosis and Mycobacterium smegmatis, we establish that the genes encoding the essential class Ib RNR subunits are specifically regulated by an NrdR-type repressor. Moreover, a strain of M. smegmatis mc(2)155 lacking the 56-kb chromosomal region, which includes duplicates of nrdHIE and nrdF2, and a mutant retaining only one copy of nrdF2 are shown to be hypersensitive to the class I RNR inhibitor hydroxyurea as a result of depleted levels of the target. Together, our observations identify a potential vulnerability in dNTP provision in mycobacteria and thereby offer a compelling rationale for pursuing the class Ib RNR as a target for drug discovery.