Amine Namouchi

After the bottleneck: Genome-wide diversification of the Mycobacterium tuberculosis complex by mutation, recombination, and natural selection

Many of the most virulent bacterial pathogens show low genetic diversity and sexual isolation. Accordingly, Mycobacterium tuberculosis, the deadliest human pathogen, is thought to be clonal and evolve by genetic drift. Yet, its genome shows few of the concomitant signs of genome degradation. We analyzed 24 genomes and found an excess of genetic diversity in regions encoding key adaptive functions including the type VII secretion system and the ancient horizontally transferred virulence-related regions. Four different approaches showed evident signs of recombination in M. tuberculosis. Recombination tracts add a high density of polymorphisms, and many are thus predicted to arise from outside the clade. Some of these tracts match Mycobacterium canettii sequences. Recombination introduced an excess of non-synonymous diversity in general and even more in genes expected to be under positive or diversifying selection, e.g., cell wall component genes. Mutations leading to non-synonymous SNPs are effectively purged in MTBC, which shows dominance of purifying selection. MTBC mutation bias toward AT nucleotides is not compensated by biased gene conversion, suggesting the action of natural selection also on synonymous changes. Together, all of these observations point to a strong imprint of recombination and selection in the genome affecting both non-synonymous and synonymous positions. Hence, contrary to some other pathogens and previous proposals concerning M. tuberculosis, this lineage may have come out of its ancestral bottleneck as a very successful pathogen that is rapidly diversifying by the action of mutation, recombination, and natural selection.

Tuberculosis due to resistant Haarlem strain, Tunisia.

Multidrug-resistant tuberculosis was diagnosed in 21 HIV-negative, nonhospitalized male patients residing in northern Tunisia. A detailed investigation showed accelerated transmission of a Mycobacterium tuberculosis clone of the Haarlem type in 90% of all patients. This finding highlights the epidemic potential of this prevalent genotype.

A Very Virulent Genotype of Infectious Bursal Disease Virus Predominantly Associated with Recurrent Infectious Bursal Disease Outbreaks in Tunisian Vaccinated Flocks

Outbreaks of infectious bursal disease (IBD) still continue to afflict the Tunisian poultry industry even in those flocks where the vaccination program is strictly applied. To characterize the viruses that circumvent protection provided by vaccination, field isolates of infectious bursal disease virus (IBDV) obtained from vaccinated flocks that have repeatedly experienced IBDV outbreak episodes were analyzed from bursal samples by reverse transcription coupled with polymerase chain reaction and dideoxynucleotide sequencing of the VP2 hypervariable region. Although sequence data were obtained from samples collected from three distinct flocks over a period of 3 years, only limited sequence variation has been observed. The few nucleotide changes were silent and the deduced amino acid sequences were identical. Thus, the virus population that predominates in the field seems to represent a homogeneous antigenic pool. Compared with the VP2 sequences of several IBDV strains, this predominant pool was found to be closely related to the very virulent (vv) IBDV viruses described in Europe and Asia. Sequence and phylogenetic analysis of the precursor polyprotein coding sequence of a representative Tunisian isolate further confirmed its assignment to the vv genotype. The deduced amino acid sequence of the whole polyprotein of the Tunisian isolate was found to be identical to a South Korean IBDV strain. Alignment of the polyprotein amino acid sequence of 35 IBDV strains identified additional mutations outside the VP2 variable domain and which occur frequently in vv strains. Based on this comparative analysis, the set of amino acid residues that should represent a typical vv profile involves Ala222, Ile242, Ile256, Ile294, Leu451, Tyr680, N685, Ser715, Asp751, Val990, and Ala1005. Such a combination of amino acid changes was observed for the majority of vvIBDV strains that define a distinct phylogroup.

Mycobacterium tuberculosis Lineage 7 Strains Are Associated with Prolonged Patient Delay in Seeking Treatment for Pulmonary Tuberculosis in Amhara Region, Ethiopia

ABSTRACT Recent genotyping studies of Mycobacterium tuberculosis in Ethiopia have reported the identification of a new phylogenetically distinct M. tuberculosis lineage, lineage 7. We therefore investigated the genetic diversity and association of specific M. tuberculosis lineages with sociodemographic and clinical parameters among pulmonary TB patients in the Amhara Region, Ethiopia. DNA was isolated from M. tuberculosis-positive sputum specimens (n = 240) and analyzed by PCR and 24-locus mycobacterial interspersed repetitive unit–variable-number tandem-repeat (MIRU-VNTR) analysis and spoligotyping. Bioinformatic analysis assigned the M. tuberculosis genotypes to global lineages, and associations between patient characteristics and genotype were evaluated using logistic regression analysis. The study revealed a high diversity of modern and premodern M. tuberculosis lineages, among which approximately 25% were not previously reported. Among the M. tuberculosis strains (n = 138) assigned to seven subgroups, the largest cluster belonged to the lineage Central Asian (CAS) (n = 60; 26.0%), the second largest to lineage 7 (n = 36; 15.6%), and the third largest to the lineage Haarlem (n = 35; 15.2%). Four sublineages were new in the MIRU-VNTRplus database, designated NW-ETH3, NW-ETH1, NW-ETH2, and NW-ETH4, which included 24 (10.4%), 18 (7.8%), 8 (3.5%), and 5 (2.2%) isolates, respectively. Notably, patient delay in seeking treatment was significantly longer among patients infected with lineage 7 strains (Mann-Whitney test, P < 0.008) than in patients infected with CAS strains (adjusted odds ratio [AOR], 4.7; 95% confidence interval [CI], 1.6 to 13.5). Lineage 7 strains also grew more slowly than other M. tuberculosis strains. Cases of Haarlem (OR, 2.8; 95% CI, 1.2 to 6.6) and NW-ETH3 (OR, 2.8; 95% CI, 1.0 to 7.3) infection appeared in defined clusters. Intensified active case finding and contact tracing activities in the study region are needed to expedite diagnosis and treatment of TB.

Evidence for the critical role of a secondary site rpoB mutation in the compensatory evolution and successful transmission of an MDR tuberculosis outbreak strain

BACKGROUND MDR Mycobacterium tuberculosis clinical strains that cause large outbreaks, particularly among HIV-negative patients, are likely to have undergone the most successful compensatory evolution. Hence, mutations secondary to the acquisition of drug resistance are worthy of consideration in these highly transmissible strains. Here, we assessed the role of a mutation within rpoB, rpoB V615M, secondary to the rifampicin resistance-conferring mutation rpoB S531L, which is associated with a major MDR tuberculosis outbreak strain that evolved in an HIV-negative context in northern Tunisia. METHODS Using BCG as a model organism, we engineered strains harbouring either the rpoB S531L mutation alone or the double mutation rpoB S531L, V615M. Individual and competitive in vitro growth assays were performed in order to assess the relative fitness of each BCG mutant. RESULTS The rpoB V615M mutation was found to be invariably associated with rpoB S531L. Structural analysis mapped rpoB V615M to the same bridge helix region as rpoB compensatory mutations previously described in Salmonella. Compared with the rpoB single-mutant BCG, the double mutant displayed improved growth characteristics and fitness rates equivalent to WT BCG. Strikingly, the rpoB double mutation conferred high-level resistance to rifampicin. CONCLUSIONS Here, we demonstrated the fitness compensatory role of a mutation within rpoB, secondary to the rifampicin resistance mutation rpoB S531L, which is characteristic of an MDR M. tuberculosis major outbreak strain. The finding that this secondary mutation concomitantly increased the resistance level to rifampicin argues for its significant contribution to the successful transmission of the MDR-TB strain.

The Mycobacterium tuberculosis transcriptional landscape under genotoxic stress

BackgroundAs an intracellular human pathogen, Mycobacterium tuberculosis (Mtb) is facing multiple stressful stimuli inside the macrophage and the granuloma. Understanding Mtb responses to stress is essential to identify new virulence factors and pathways that play a role in the survival of the tubercle bacillus. The main goal of this study was to map the regulatory networks of differentially expressed (DE) transcripts in Mtb upon various forms of genotoxic stress. We exposed Mtb cells to oxidative (H2O2 or paraquat), nitrosative (DETA/NO), or alkylation (MNNG) stress or mitomycin C, inducing double-strand breaks in the DNA. Total RNA was isolated from treated and untreated cells and subjected to high-throughput deep sequencing. The data generated was analysed to identify DE genes encoding mRNAs, non-coding RNAs (ncRNAs), and the genes potentially targeted by ncRNAs.ResultsThe most significant transcriptomic alteration with more than 700 DE genes was seen under nitrosative stress. In addition to genes that belong to the replication, recombination and repair (3R) group, mainly found under mitomycin C stress, we identified DE genes important for bacterial virulence and survival, such as genes of the type VII secretion system (T7SS) and the proline-glutamic acid/proline-proline-glutamic acid (PE/PPE) family. By predicting the structures of hypothetical proteins (HPs) encoded by DE genes, we found that some of these HPs might be involved in mycobacterial genome maintenance. We also applied a state-of-the-art method to predict potential target genes of the identified ncRNAs and found that some of these could regulate several genes that might be directly involved in the response to genotoxic stress.ConclusionsOur study reflects the complexity of the response of Mtb in handling genotoxic stress. In addition to genes involved in genome maintenance, other potential key players, such as the members of the T7SS and PE/PPE gene family, were identified. This plethora of responses is detected not only at the level of DE genes encoding mRNAs but also at the level of ncRNAs and their potential targets.

Evolutionary Trends of the Transposase-Encoding Open Reading Frames A and B (orfA and orfB) of the Mycobacterial IS6110 Insertion Sequence

Background The IS6110 insertion sequence, a member of the IS3 family of insertion sequences, was found to be specific to the Mycobacterium tuberculosis complex (MTBC). Although IS6110 has been extensively characterized as a transposable genetic marker, the evolutionary history of its own transposase-encoding sequence has not, to the best of our knowledge, been investigated. Methodology/Principal Findings Here we explored the evolution of the IS6110 sequence by analysing the genetic variability and the selective forces acting on its transposase-encoding open reading frames (ORFs) A and B (orfA and orfB). For this purpose, we used a strain collection consisting of smooth tubercle bacilli (STB), an early branching lineage of the MTBC, and present-day M. tuberculosis strains representing the full breadth of genetic diversity in Tunisia. In each ORF, we found a major haplotype that dominated over a flat distribution of rare descendent haplotypes, consisting mainly of single- and double-nucleotide variant singletons. The predominant haplotypes consisted of both ancestral and present-day strains, suggesting that IS6110 acquisition predated the emergence of the MTBC. There was no evidence of recombination and both ORFs were subjected to strict purifying selection, as demonstrated by their dN/dS ratios (0.29 and 0.51, respectively), as well as their significantly negative Tajima’s D statistics. Strikingly, the purifying selection acting on orfA proved much more stringent, suggesting its critical role in regulating the transpositional process. Maximum likelihood analyses further excluded any possibility of positive selection acting on single amino acid residues. Conclusions/Significance Taken together our data fit with an evolutionary scenario according to which the observed variability pattern of the IS6110 transposase-encoding ORFs is generated mainly through random point mutations that accrued on a functionally optimal IS6110 copy, whose acquisition predated the emergence of the MTBC complex. Background selection acting against deleterious mutations led to an excess of low-frequency variants.

Characterization of the Neisseria meningitidis Helicase RecG

Neisseria meningitidis (Nm) is a Gram-negative oral commensal that opportunistically can cause septicaemia and/or meningitis. Here, we overexpressed, purified and characterized the Nm DNA repair/recombination helicase RecG (RecGNm) and examined its role during genotoxic stress. RecGNm possessed ATP-dependent DNA binding and unwinding activities in vitro on a variety of DNA model substrates including a Holliday junction (HJ). Database searching of the Nm genomes identified 49 single nucleotide polymorphisms (SNPs) in the recGNm including 37 non-synonymous SNPs (nsSNPs), and 7 of the nsSNPs were located in the codons for conserved active site residues of RecGNm. A transient reduction in transformation of DNA was observed in the Nm ΔrecG strain as compared to the wildtype. The gene encoding recGNm also contained an unusually high number of the DNA uptake sequence (DUS) that facilitate transformation in neisserial species. The differentially abundant protein profiles of the Nm wildtype and ΔrecG strains suggest that expression of RecGNm might be linked to expression of other proteins involved in DNA repair, recombination and replication, pilus biogenesis, glycan biosynthesis and ribosomal activity. This might explain the growth defect that was observed in the Nm ΔrecG null mutant.

The helicase DinG responds to stress due to DNA double strand breaks

Neisseria meningitidis (Nm) is a Gram-negative nasopharyngeal commensal that can cause septicaemia and meningitis. The neisserial DNA damage-inducible protein DinG is a helicase related to the mammalian helicases XPD and FANCJ. These helicases belong to superfamily 2, are ATP dependent and exert 5′ → 3′ directionality. To better understand the role of DinG in neisserial genome maintenance, the Nm DinG (DinGNm) enzymatic activities were assessed in vitro and phenotypical characterization of a dinG null mutant (NmΔdinG) was performed. Like its homologues, DinGNm possesses 5′ → 3′ directionality and prefers DNA substrates containing a 5′-overhang. ATPase activity of DinGNm is strictly DNA-dependent and DNA unwinding activity requires nucleoside triphosphate and divalent metal cations. DinGNm directly binds SSBNm with a Kd of 313 nM. Genotoxic stress analysis demonstrated that NmΔdinG was more sensitive to double-strand DNA breaks (DSB) induced by mitomycin C (MMC) than the Nm wildtype, defining the role of neisserial DinG in DSB repair. Notably, when NmΔdinG cells grown under MMC stress assessed by quantitative mass spectrometry, 134 proteins were shown to be differentially abundant (DA) compared to unstressed NmΔdinG cells. Among the DNA replication, repair and recombination proteins affected, polymerase III subunits and recombinational repair proteins RuvA, RuvB, RecB and RecD were significantly down regulated while TopA and SSB were upregulated under stress condition. Most of the other DA proteins detected are involved in metabolic functions. The present study shows that the helicase DinG is probably involved in regulating metabolic pathways as well as in genome maintenance.

Phenotypic and genomic comparison of Mycobacterium aurum and surrogate model species to Mycobacterium tuberculosis: implications for drug discovery

BackgroundTuberculosis (TB) is caused by Mycobacterium tuberculosis and represents one of the major challenges facing drug discovery initiatives worldwide. The considerable rise in bacterial drug resistance in recent years has led to the need of new drugs and drug regimens. Model systems are regularly used to speed-up the drug discovery process and circumvent biosafety issues associated with manipulating M. tuberculosis. These include the use of strains such as Mycobacterium smegmatis and Mycobacterium marinum that can be handled in biosafety level 2 facilities, making high-throughput screening feasible. However, each of these model species have their own limitations.ResultsWe report and describe the first complete genome sequence of Mycobacterium aurum ATCC23366, an environmental mycobacterium that can also grow in the gut of humans and animals as part of the microbiota. This species shows a comparable resistance profile to that of M. tuberculosis for several anti-TB drugs. The aims of this study were to (i) determine the drug resistance profile of a recently proposed model species, Mycobacterium aurum, strain ATCC23366, for anti-TB drug discovery as well as Mycobacterium smegmatis and Mycobacterium marinum (ii) sequence and annotate the complete genome sequence of this species obtained using Pacific Bioscience technology (iii) perform comparative genomics analyses of the various surrogate strains with M. tuberculosis (iv) discuss how the choice of the surrogate model used for drug screening can affect the drug discovery process.ConclusionsWe describe the complete genome sequence of M. aurum, a surrogate model for anti-tuberculosis drug discovery. Most of the genes already reported to be associated with drug resistance are shared between all the surrogate strains and M. tuberculosis. We consider that M. aurum might be used in high-throughput screening for tuberculosis drug discovery. We also highly recommend the use of different model species during the drug discovery screening process.