A regulatory variant in the C1Q gene cluster is associated with tuberculosis susceptibility and C1qA plasma levels in a South African population.

Several genetic studies have implicated genes that encode for components of the innate immune response in tuberculosis (TB) susceptibility. The complement system is an early player in the innate immune response and provides the host with initial protection by promoting phagocytosis of apoptotic or necrotic cells. The C1q molecule is the first component of the classical pathway that leads to the activation of complement by binding to immune complexes and is encoded by the C1Q gene cluster. We investigated variants in this region to determine its association with TB susceptibility. Five single nucleotide polymorphisms (SNPs) (rs12033074, rs631090, rs172378, rs587585, and rs665691) were genotyped using TaqMan® SNP assays in 456 TB cases and 448 healthy controls and analysed by logistic regression models. The rs587585 variant showed a significant additive allelic association where the minor G allele was found more frequently in TB cases than in controls in both the discovery (p = 0.023; OR = 1.30; 95% CI, 1.04-1.64) and validation cohort (p = 0.038; OR = 1.31; 95% CI, 1.22-1.40). In addition, we detected increased C1qA expression when comparing cases and controls (p = 0.037) and linked this to a dosage effect of the G allele, which increased C1qA expression in TB cases. This is the first study to report the association of C1Q gene polymorphisms with progression to tuberculosis.

Proteogenomic Investigation of Strain Variation in Clinical Mycobacterium tuberculosis Isolates.

Mycobacterium tuberculosis consists of a large number of different strains that display unique virulence characteristics. Whole-genome sequencing has revealed substantial genetic diversity among clinical M. tuberculosis isolates, and elucidating the phenotypic variation encoded by this genetic diversity will be of the utmost importance to fully understand M. tuberculosis biology and pathogenicity. In this study, we integrated whole-genome sequencing and mass spectrometry (GeLC-MS/MS) to reveal strain-specific characteristics in the proteomes of two clinical M. tuberculosis Latin American-Mediterranean isolates. Using this approach, we identified 59 peptides containing single amino acid variants, which covered ∼9% of all coding nonsynonymous single nucleotide variants detected by whole-genome sequencing. Furthermore, we identified 29 distinct peptides that mapped to a hypothetical protein not present in the M. tuberculosis H37Rv reference proteome. Here, we provide evidence for the expression of this protein in the clinical M. tuberculosis SAWC3651 isolate. The strain-specific databases enabled confirmation of genomic differences (i.e., large genomic regions of difference and nonsynonymous single nucleotide variants) in these two clinical M. tuberculosis isolates and allowed strain differentiation at the proteome level. Our results contribute to the growing field of clinical microbial proteogenomics and can improve our understanding of phenotypic variation in clinical M. tuberculosis isolates.

Recombination in pe/ppe genes contributes to genetic variation in Mycobacterium tuberculosis lineages.

BACKGROUND: Approximately 10% of the Mycobacterium tuberculosis genome is made up of two families of genes that are poorly characterized due to their high GC content and highly repetitive nature. The PE and PPE families are typified by their highly conserved N-terminal domains that incorporate proline-glutamate (PE) and proline-proline-glutamate (PPE) signature motifs. They are hypothesised to be important virulence factors involved with host-pathogen interactions, but their high genetic variability and complexity of analysis means they are typically disregarded in genome studies. RESULTS: To elucidate the structure of these genes, 518 genomes from a diverse international collection of clinical isolates were de novo assembled. A further 21 reference M. tuberculosis complex genomes and long read sequence data were used to validate the approach. SNP analysis revealed that variation in the majority of the 168 pe/ppe genes studied was consistent with lineage. Several recombination hotspots were identified, notably pe_pgrs3 and pe_pgrs17. Evidence of positive selection was revealed in 65 pe/ppe genes, including epitopes potentially binding to major histocompatibility complex molecules. CONCLUSIONS: This, the first comprehensive study of the pe and ppe genes, provides important insight into M. tuberculosis diversity and has significant implications for vaccine development.

A metabolomics investigation of the function of the ESX-1 gene cluster in mycobacteria.

The ESX-1 gene cluster, encoding the Type-VII secretion (T7S) system and its virulence associated proteins, ESAT-6 and CFP-10, is thought to be responsible for the transport of extracellular proteins across the hydrophobic and highly impermeable, cell envelope of Mycobacterium, and is involved in virulence in Mycobacterium tuberculosis, the causative agent of tuberculosis. Using a GCxGC-TOFMS metabolomics approach, a M. smegmatis ESX-1 knock-out strain (ΔESX-1ms) was compared to that of the M. smegmatis wild-type parent strain, and the metabolite markers due to the presence or absence of the ESX-1 gene cluster were identified. A general increase in specific metabolites in the ΔESX-1ms, confirmed the roles previously described for ESX-1 in mycolic acid biosynthesis and cell wall integrity. However, a number of other metabolite markers identified indicates ESX-1 has an additional role the in cell envelope structure, altering the levels of antioxidants and energy metabolism. Furthermore, the metabolome profiles correlated with the metabolomic variation observed when comparing a hyper- and hypo-virulent Beijing strain of M. tuberculosis, suggesting that the pathways which modulate virulence in M. tuberculosis are also influenced by ESX-1, reaffirming the previously described association of ESX-1 with virulence and cell envelope biogenesis.

Novel cause of tuberculosis in meerkats, South Africa.

The organism that causes tuberculosis in meerkats (Suricata suricatta) has been poorly characterized. Our genetic analysis showed it to be a novel member of the Mycobacterium tuberculosis complex and closely related to the dassie bacillus. We have named this epidemiologically and genetically unique strain M. suricattae.

Rifampicin reduces susceptibility to ofloxacin in rifampicin-resistant Mycobacterium tuberculosis through efflux.

RATIONALE: Central dogma suggests that rifampicin resistance in Mycobacterium tuberculosis develops solely through rpoB gene mutations. OBJECTIVE: To determine whether rifampicin induces efflux pumps activation in rifampicin resistant M. tuberculosis strains thereby defining rifampicin resistance levels and reducing ofloxacin susceptibility. METHODS: Rifampicin and/or ofloxacin minimum inhibitory concentrations (MICs) were determined in rifampicin resistant strains by culture in BACTEC 12B medium. Verapamil and reserpine were included to determine their effect on rifampicin and ofloxacin susceptibility. RT-qPCR was applied to assess expression of efflux pump/transporter genes after rifampicin exposure. To determine whether verapamil could restore susceptibility to first-line drugs, BALB/c mice were infected with a MDR-TB strain and treated with first-line drugs with/without verapamil. MEASUREMENTS AND MAIN FINDINGS: Rifampicin MICs varied independently of rpoB mutation and genetic background. Addition reserpine and verapamil significantly restored rifampicin susceptibility (p = 0.0000). RT-qPCR demonstrated that rifampicin induced differential expression of efflux/transporter genes in MDR-TB isolates. Incubation of rifampicin mono-resistant strains in rifampicin (2 µg/ml) for 7 days induced ofloxacin resistance (MIC > 2 µg/ml) in strains with an rpoB531 mutation. Ofloxacin susceptibility was restored by exposure to efflux pump inhibitors. Studies in BALB/c mice showed that verapamil in combination with first-line drugs significantly reduced pulmonary CFUs after 1 and 2 months treatment (p < 0.05). CONCLUSION: Exposure of rifampicin resistant M. tuberculosis strains to rifampicin can potentially compromise the efficacy of the second-line treatment regimens containing ofloxacin, thereby emphasising the need for rapid diagnostics to guide treatment. Efflux pump inhibitors have the potential to improve the efficacy of anti-tuberculosis drug treatment.

The non-clonality of drug resistance in Beijing-genotype isolates of Mycobacterium tuberculosis from the Western Cape of South Africa.

BACKGROUND: The Beijing genotype of M. tuberculosis is a virulent strain that is disseminating worldwide and has a strong association with drug resistance. In the Western Cape of South Africa, epidemiological studies have identified the R220 cluster of the Beijing genotype as a major contributor to a recent outbreak of drug-resistant tuberculosis. Although the outbreak is considered to be due to clonal transmission, the relationship among drug resistant isolates has not yet been established. RESULTS: To better understand the evolution of drug resistance among these strains, 14 drug-resistant clinical isolates of the Beijing genotype were sequenced by whole-genome sequencing, including eight from R220 and six from a more ancestral Beijing cluster, R86, for comparison. While each cluster shares a distinct resistance mutation for isoniazid, mapping of other drug-resistance mutations onto a phylogenetic tree constructed from single nucleotide polymorphisms shows that resistance mutations to many drugs have arisen multiple times independently within each cluster of isolates. Thus, drug resistance among these isolates appears to be acquired, not clonally derived. This observation suggests that, although the Beijing genotype as a whole might have selective advantages enabling its rapid dissemination, the XDR isolates are relatively less fit and do not propagate well. Although it has been hypothesized that the increased frequency of drug resistance in some Beijing lineages might be caused by a mutator phenotype, no significant shift in synonymous substitution patterns is observed in the genomes. CONCLUSION: While MDR-TB is spreading by transmission in the Western Cape, our data suggests that further drug resistance (i.e. XDR-TB) at this stage is acquired.

Novel Mycobacterium tuberculosis complex pathogen, M. mungi.

Seven outbreaks involving increasing numbers of banded mongoose troops and high death rates have been documented. We identified a Mycobacterium tuberculosis complex pathogen, M. mungi sp. nov., as the causative agent among banded mongooses that live near humans in Chobe District, Botswana. Host spectrum and transmission dynamics remain unknown.

PPE and PE_PGRS proteins of Mycobacterium marinum are transported via the type VII secretion system ESX-5.

ESX-5 is one of the five type VII secretion systems found in mycobacteria. These secretion systems are also known as ESAT-6-like secretion systems. Here, we have determined the secretome of ESX-5 by a proteomic approach in two different strains of Mycobacterium marinum. Comparison of the secretion profile of wild-type strains and their ESX-5 mutants showed that a number of PE_PGRS and PPE-MPTR proteins are dependent on ESX-5 for transport. The PE and PPE protein families are unique to mycobacteria, are highly expanded in several pathogenic species, such as Mycobacterium tuberculosis and M. marinum, and certain family members are cell surface antigens associated with virulence. Using a monoclonal antibody directed against the PGRS domain we showed that nearly all PE_PGRS proteins that are recognized by this antibody are missing in the supernatant of ESX-5 mutants. In addition to PE_PGRS and PPE proteins, the ESX-5 secretion system is responsible for the secretion of a ESAT-6-like proteins. Together, these data show that ESX-5 is probably a major secretion pathway for mycobacteria and that this system is responsible for the secretion of recently evolved PE_PGRS and PPE proteins.

Evaluation of the Capilia TB assay for culture confirmation of Mycobacterium tuberculosis infections in Zambia and South Africa.

The performance and cost of the Capilia TB assay were evaluated for use in a resource-limited setting. The sensitivity and specificity were 99.6% and 99.5%, respectively. The incremental costs of the Capilia test were estimated to be $1.46 and $1.84 when the test was added to liquid and solid culture processes, respectively. These findings suggest that the Capilia TB assay represents a rapid, simple, and inexpensive Mycobacterium tuberculosis identification test that can be used in resource-limited settings.

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.

Frequent homologous recombination events in Mycobacterium tuberculosis PE/PPE multigene families: potential role in antigenic variability.

The PE and PPE (PE/PPE) multigene families of Mycobacterium tuberculosis are particularly GC-rich and share extensive homologous repetitive sequences. We hypothesized that they may undergo homologous recombination events, a mechanism rarely described in the natural evolution of mycobacteria. To test our hypothesis, we developed a specific oligonucleotide-based microarray targeting nearly all of the PE/PPE genes, aimed at detecting signals for homologous recombination. Such a microarray has never before been reported due to the multiplicity and highly repetitive and homologous nature of these sequences. Application of the microarray to a collection of M. tuberculosis clinical isolates (n = 33) representing prevalent spoligotype strain families in Tunisia allowed successful detection of six deleted genomic regions involving a total of two PE and seven PPE genes. Some of these deleted genes are known to be immunodominant or involved in virulence. The four precisely determined deletions were flanked by 400- to 500-bp stretches of nearly identical sequences lying mainly at the conserved N-terminal region of the PE/PPE genes. These highly homologous sequences thus serve as substrates to mediate both intergenic and intragenic homologous recombination events, indicating an important function in generating strain variation. Importantly, all recombination events yielded a new in-frame fusion chimeric gene. Hence, homologous recombination within and between PE/PPE genes likely increased their antigenic variability, which may have profound implications in pathogenicity and/or host adaptation. The finding of high prevalence (approximately 45% and approximately 58%) for at least two of the genomic deletions suggests that they likely confer advantageous biological attributes.

Application of sensitive and specific molecular methods to uncover global dissemination of the major RDRio Sublineage of the Latin American-Mediterranean Mycobacterium tuberculosis spoligotype family.

The Latin American-Mediterranean (LAM) family of Mycobacterium tuberculosis is believed to be the cause of approximately 15% of tuberculosis cases worldwide. Previously, we defined a prevalent sublineage of the LAM family in Brazil by a single characteristic genomic deletion designated RD(Rio). Using the Brazilian strains, we pinpoint an Ag85C(103) single nucleotide polymorphism (SNP) (screened by restriction fragment length polymorphism [RFLP] analysis) that correctly identified all LAM family strains. Importantly, all RD(Rio) strains concomitantly possessed the RD174 deletion. These genetic signatures, along with a newly developed multiplex PCR for rapid differentiation between "wild-type" and RD(Rio) strains, were then used to analyze an international collection of M. tuberculosis strains. RD(Rio) M. tuberculosis was identified from four continents involving 11 countries. Phylogenetic analysis of the IS6110-RFLP patterns from representative RD(Rio) and LAM strains from Brazil, along with all representative clusters from a South African database, confirmed their genetic relatedness and transcontinental transmission. The Ag85C(103) SNP RFLP, as compared to results obtained using a PCR method targeting a LAM-restricted IS6110 element, correctly identified 99.8% of LAM spoligotype strains. Together, these tests were more accurate than spoligotyping at categorizing strains with indefinable spoligotypes and segregated true LAM strains from those with convergent spoligotypes. The fact that RD(Rio) strains were identified worldwide highlights the importance of this LAM family sublineage and suggests that this strain is a global threat that should be specifically targeted by public health resources. Our provision of simple and robust molecular methods will assist the evaluation of the LAM family and the RD(Rio) sublineage.

Pulmonary infection due to the dassie bacillus (Mycobacterium tuberculosis complex sp.) in a free-living dassie (rock hyrax-Procavia capensis) from South Africa.

We report a case of extensive necrogranulomatous pneumonia due to infection with the dassie bacillus (Mycobacterium tuberculosis complex sp.) in a free-living pregnant adult female dassie (rock hyrax-Procavia capensis). A juvenile female dassie from the same colony also showed a focal lesion in the lungs suggestive of mycobacterial pneumonia. Our findings indicate the widespread occurrence of the dassie bacillus in free-living dassies and suggest very high infection rates in some populations. The introduction of South African dassies into novel environments should be considered in this light.

Geospatial distribution of Mycobacterium tuberculosis genotypes in Africa.

OBJECTIVE: To investigate the distribution of Mycobacterium tuberculosis genotypes across Africa. METHODS: The SITVIT2 global repository and PUBMED were searched for spoligotype and published genotype data respectively, of M. tuberculosis from Africa. M. tuberculosis lineages in Africa were described and compared across regions and with those from 7 European and 6 South-Asian countries. Further analysis of the major lineages and sub-lineages using Principal Component analysis (PCA) and hierarchical cluster analysis were done to describe clustering by geographical regions. Evolutionary relationships were assessed using phylogenetic tree analysis. RESULTS: A total of 14727 isolates from 35 African countries were included in the analysis and of these 13607 were assigned to one of 10 major lineages, whilst 1120 were unknown. There were differences in geographical distribution of major lineages and their sub-lineages with regional clustering. Southern African countries were grouped based on high prevalence of LAM11-ZWE strains; strains which have an origin in Portugal. The grouping of North African countries was due to the high percentage of LAM9 strains, which have an origin in the Eastern Mediterranean region. East African countries were grouped based on Central Asian (CAS) and East-African Indian (EAI) strain lineage possibly reflecting historic sea trade with Asia, while West African Countries were grouped based on Cameroon lineage of unknown origin. A high percentage of the Haarlem lineage isolates were observed in the Central African Republic, Guinea, Gambia and Tunisia, however, a mixed distribution prevented close clustering. CONCLUSIONS: This study highlighted that the TB epidemic in Africa is driven by regional epidemics characterized by genetically distinct lineages of M. tuberculosis. M. tuberculosis in these regions may have been introduced from either Europe or Asia and has spread through pastoralism, mining and war. The vast array of genotypes and their associated phenotypes should be considered when designing future vaccines, diagnostics and anti-TB drugs.

The role of IS6110 in the evolution of Mycobacterium tuberculosis.

Members of the Mycobacterium tuberculosis complex contain the transposable element IS6110 which, due to its high numerical and positional polymorphism, has become a widely used marker in epidemiological studies. Here, we review the evidence that IS6110 is not simply a passive or 'junk' DNA sequence, but that, through its transposable activity, it is able to generate genotypic variation that translates into strain-specific phenotypic variation. We also speculate on the role that this variation has played in the evolution of M. tuberculosis and conclude that the presence of a moderate IS6110 copy number within the genome may provide the pathogen with a selective advantage that has aided its virulence.

Evolution and expansion of the Mycobacterium tuberculosis PE and PPE multigene families and their association with the duplication of the ESAT-6 (esx) gene cluster regions.

BACKGROUND: The PE and PPE multigene families of Mycobacterium tuberculosis comprise about 10% of the coding potential of the genome. The function of the proteins encoded by these large gene families remains unknown, although they have been proposed to be involved in antigenic variation and disease pathogenesis. Interestingly, some members of the PE and PPE families are associated with the ESAT-6 (esx) gene cluster regions, which are regions of immunopathogenic importance, and encode a system dedicated to the secretion of members of the potent T-cell antigen ESAT-6 family. This study investigates the duplication characteristics of the PE and PPE gene families and their association with the ESAT-6 gene clusters, using a combination of phylogenetic analyses, DNA hybridization, and comparative genomics, in order to gain insight into their evolutionary history and distribution in the genus Mycobacterium. RESULTS: The results showed that the expansion of the PE and PPE gene families is linked to the duplications of the ESAT-6 gene clusters, and that members situated in and associated with the clusters represent the most ancestral copies of the two gene families. Furthermore, the emergence of the repeat protein PGRS and MPTR subfamilies is a recent evolutionary event, occurring at defined branching points in the evolution of the genus Mycobacterium. These gene subfamilies are thus present in multiple copies only in the members of the M. tuberculosis complex and close relatives. The study provides a complete analysis of all the PE and PPE genes found in the sequenced genomes of members of the genus Mycobacterium such as M. smegmatis, M. avium paratuberculosis, M. leprae, M. ulcerans, and M. tuberculosis. CONCLUSION: This work provides insight into the evolutionary history for the PE and PPE gene families of the mycobacteria, linking the expansion of these families to the duplications of the ESAT-6 (esx) gene cluster regions, and showing that they are composed of subgroups with distinct evolutionary (and possibly functional) differences.

Insights into the evolutionary history of tubercle bacilli as disclosed by genetic rearrangements within a PE_PGRS duplicated gene pair.

BACKGROUND: The highly homologous PE_PGRS (Proline-glutamic acid_polymorphic GC-rich repetitive sequence) genes are members of the PE multigene family which is found only in mycobacteria. PE genes are particularly abundant within the genomes of pathogenic mycobacteria where they seem to have expanded as a result of gene duplication events. PE_PGRS genes are characterized by their high GC content and extensive repetitive sequences, making them prone to recombination events and genetic variability. RESULTS: Comparative sequence analysis of Mycobacterium tuberculosis genes PE_PGRS17 (Rv0978c) and PE_PGRS18 (Rv0980c) revealed a striking genetic variation associated with this typical tandem duplicate. In comparison to the M. tuberculosis reference strain H37Rv, the variation (named the 12/40 polymorphism) consists of an in-frame 12-bp insertion invariably accompanied by a set of 40 single nucleotide polymorphisms (SNPs) that occurs either in PE_PGRS17 or in both genes. Sequence analysis of the paralogous genes in a representative set of worldwide distributed tubercle bacilli isolates revealed data which supported previously proposed evolutionary scenarios for the M. tuberculosis complex (MTBC) and confirmed the very ancient origin of "M. canettii" and other smooth tubercle bacilli. Strikingly, the identified polymorphism appears to be coincident with the emergence of the post-bottleneck successful clone from which the MTBC expanded. Furthermore, the findings provide direct and clear evidence for the natural occurrence of gene conversion in mycobacteria, which appears to be restricted to modern M. tuberculosis strains. CONCLUSION: This study provides a new perspective to explore the molecular events that accompanied the evolution, clonal expansion, and recent diversification of tubercle bacilli.