The Impact of Genome Region of Difference 4 (RD4) on Mycobacterial Virulence and BCG Efficacy

Comprehensive analysis of Mycobacterium tuberculosis genomes reveals genetic variations in bacterial virulence

Tuberculosis, a disease caused by Mycobacterium tuberculosis (Mtb), is a significant health problem worldwide. Here, we developed a method to detect large insertions and deletions (indels), which have been generally understudied. Leveraging this framework, we performed a comprehensive analysis of single nucleotide variants and small and large indels across 1,960 Mtb clinical isolates. Comparing the distribution of variants demonstrated that gene disruptive variants are underrepresented in genes essential for bacterial survival. An evolutionary analysis revealed that Mtb genomes are enriched in partially deleterious mutations. Genome-wide association studies identified small and large deletions in eccB2 significantly associated with patient prognosis. Additionally, we unveil significant associations with antibiotic resistance in 23 non-canonical genes. Among these, large indels are primarily found in genetic regions of Rv1216c, Rv1217c, fadD11, and ctpD. This study provides a comprehensive catalog of genetic variations and highlights their potential impact for the future treatment and risk prediction of tuberculosis.

Accuracy of molecular diagnostic assays for detection of Mycobacterium bovis: A systematic review and meta-analysis

Bovine tuberculosis (bovine TB) is a chronic wasting disease of cattle caused primarily by Mycobacterium bovis. Controlling bovine TB requires highly sensitive, specific, quick, and reliable diagnostic methods. This systematic review and meta-analysis evaluated molecular diagnostic tests for M. bovis detection to inform the selection of the most viable assay. On a per-test basis, loop-mediated isothermal amplification (LAMP) showed the highest overall sensitivity of 99.0% [95% CI: 86.2%-99.9%] and specificity of 99.8% [95% CI: 96.2%-100.00%]. Quantitative real-time polymerase chain reaction (qPCR) outperformed conventional PCR and nested PCR (nPCR) with a diagnostic specificity of 96.6% [95% CI: 88.9%-99.0%], while the diagnostic sensitivity of 70.8% [95% CI: 58.6-80.5%] was comparable to that of nPCR at 71.4% [95% CI: 60.7-80.2%]. Test sensitivity was higher with the input of milk samples (90.9% [95% CI: 56.0%-98.7%]), while specificity improved with tests based on major M. bovis antigens (97.8% [95% CI: 92.3%-99.4%]), the IS6110 insertion sequence (95.4% [95% CI: 87.6%-98.4%]), and the RD4 gene (90.7% [95% CI: 52.2%-98.9%]). The design of the currently available molecular diagnostic assays, while mostly based on nonspecific gene targets, prevents them from being accurate enough to diagnose M. bovis infections in cattle, despite their promise. Future assay development should focus on the RD4 region since it is the only target identified by genome sequence data as being distinctive for detecting M. bovis. The availability of a sufficiently accurate diagnostic test combined with the routine screening of milk samples can decrease the risk of zoonotic transmissions of M. bovis.

Recent progress in the genotyping of bovine tuberculosis and its rapid diagnosis via nanoparticle-based electrochemical biosensors

Bovine tuberculosis (bTB) is considered a worldwide infectious zoonotic disease. Mycobacterium bovis causes bTB disease. It is one of the Mycobacterium tuberculosis complex (MTBC) members. MTBC is a clonal complex of close relatives with approximately 99.95% similarity. M. bovis is a spillover pathogen that can transmit from animals to humans and rarely from humans to animals with contact. Genotyping techniques are important to discriminate and differentiate between MTBC species. Spoligotyping and mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) are widely used but they have some limitations. As an alternative, whole genome sequencing approaches have been utilized due to their high-resolution power. They are employed in typing M. bovis and explain the evolutionary and phylogenetic relationships between isolates. The control of bTB disease has attracted a large amount of attention. Rapid and proper diagnosis is necessary for monitoring the disease as an initial step for its control and treatment. Nanotechnology has a potential impact on the rapid diagnosis and treatment of bTB through the use of nanocarrier and metal nanoparticles (NPs). Special attention has been paid to voltammetric and impedimetric electrochemical strategies as facile, sensitive, and selective methods for the efficient detection of tuberculosis. The efficacy of these sensors is enhanced in the presence of NPs, which act as recognition and/or redox probes. Gold, silver, copper, cobalt, graphene, and magnetic NPs, as well as polypyrrole nanowires and multiwalled carbon nanotubes have been employed for detecting tuberculosis. Overall, NP-based electrochemical sensors represent a promising tool for the diagnosis of bTB.

Mapping the phylogeny and lineage history of geographically distinct BCG vaccine strains

The bacillus Calmette-Guérin (BCG) vaccine has been in use for prevention of tuberculosis for over a century. It remains the only widely available tuberculosis vaccine and its protective efficacy has varied across geographical regions. Since it was developed, the BCG vaccine strain has been shared across different laboratories around the world, where use of differing culture methods has resulted in genetically distinct strains over time. Whilst differing BCG vaccine efficacy around the world is well documented, and the reasons for this may be multifactorial, it has been hypothesized that genetic differences in BCG vaccine strains contribute to this variation. Isolates from an historic archive of lyophilized BCG strains were regrown, DNA was extracted and then whole-genome sequenced using Oxford Nanopore Technologies. The resulting whole-genome data were plotted on a phylogenetic tree and analysed to identify the presence or absence of regions of difference (RDs) and single-nucleotide polymorphisms (SNPs) relating to virulence, growth and cell wall structure. Of 50 strains available, 36 were revived in culture and 39 were sequenced. Morphology differed between the strains distributed before and after 1934. There was phylogenetic association amongst certain geographically classified strains, most notably BCG-Russia, BCG-Japan and BCG-Danish. RD2, RD171 and RD713 deletions were associated with late strains (seeded after 1927). When mapped to BCG-Pasteur 1172, the SNPs in sigK, plaA, mmaA3 and eccC5 were associated with early strains. Whilst BCG-Russia, BCG-Japan and BCG-Danish showed strong geographical isolate clustering, the late strains, including BCG-Pasteur, showed more variation. A wide range of SNPs were seen within geographically classified strains, and as much intra-strain variation as between-strain variation was seen. The date of distribution from the original Pasteur laboratory (early pre-1927 or late post-1927) gave the strongest association with genetic differences in regions of difference and virulence-related SNPs, which agrees with the previous literature.

Cryptococcus neoformans Infection Induces IL-17 Production by Promoting STAT3 Phosphorylation in CD4+ T Cells

Cryptococcus neoformans infection in the central nervous system is a severe infectious disease with poor outcomes and high mortality. It has been estimated that there are 220,000 new cases each year. Over 90% of C. neoformans meningitis cases were diagnosed in AIDS patients with CD4⁺ T cell count <100 cells/μl; however, the mechanism of cryptococcal meningitis in patients with normal immune functions remains unclear. IL-17 is a pro-inflammatory cytokine and plays an important role in anti-fungal immunity. Here we report that significantly high levels of IL-17 were predominantly detected in the cerebrospinal fluid of patients with either AIDS- or non-AIDS-associated C. neoformans meningitis but not in patients with tuberculous meningitis or non-neurosyphilis. Antifungal therapy minimized the IL-17 level in the cerebrospinal fluid. An in vitro mechanistic study showed that C. neoformans stimulation of healthy peripheral blood mononuclear cells prompted IL-17 production, and CD4⁺ T cells were the predominant IL-17-producing cells. IL-17 production by C. neoformans stimulation was STAT3 signaling dependent. Inhibition of STAT3 phosphorylation attenuated the C. neoformans-mediated IL-17 expression. Our data highlighted the significance of CD4⁺ T cells in antifungal immunity and suggested IL-17 as a diagnostic biomarker of C. neoformans infection and STAT3 as a checkpoint for antifungal targeted therapies.

BCG-induced trained immunity enhances acellular pertussis vaccination responses in an explorative randomized clinical trial

Acellular pertussis (aP) booster vaccines are central to pertussis immunization programs, although their effectiveness varies. The Bacille Calmette-Guérin (BCG) vaccine is a prototype inducer of trained immunity, which enhances immune responses to subsequent infections or vaccinations. While previous clinical studies have demonstrated that trained immunity can protect against heterologous infections, its effect on aP vaccines in humans is unknown. We conducted a clinical study in order to determine the immunological effects of trained immunity on pertussis vaccination. Healthy female volunteers were randomly assigned to either receive BCG followed by a booster dose of tetanus-diphteria-pertussis inactivated polio vaccine (Tdap-IPV) 3 months later (BCG-trained), BCG + Tdap-IPV concurrently, or Tdap-IPV followed by BCG 3 months later. Primary outcomes were pertussis-specific humoral, T- and B-cell responses and were quantified at baseline of Tdap-IPV vaccination and 2 weeks thereafter. As a secondary outcome in the BCG-trained cohort, ex vivo leukocyte responses were measured in response to unrelated stimuli before and after BCG vaccination. BCG vaccination 3 months prior to, but not concurrent with, Tdap-IPV improves pertussis-specific Th1-cell and humoral responses, and also increases total memory B cell responses. These responses were correlated with enhanced IL-6 and IL-1β production at the baseline of Tdap-IPV vaccination in the BCG-trained cohort. Our study demonstrates that prior BCG vaccination potentiates immune responses to pertussis vaccines and that biomarkers of trained immunity are the most reliable correlates of those responses.

Interaction analysis of Mycobacterium tuberculosis between the host environment and highly mutated genes from population genetic structure comparison

We aimed to investigate the genetic and demographic differences and interactions between areas where observed genomic variations in Mycobacterium tuberculosis (M. tb) were distributed uniformly in cold and hot spots.The cold and hot spot areas were identified using the reported incidence of TB over the previous 5 years. Whole genome sequencing was performed on 291 M. tb isolates between January and June 2018. Analysis of molecular variance and a multifactor dimensionality reduction (MDR) model was applied to test gene-gene-environment interactions. Adjusted odds ratios (OR) and 95% confidence intervals (CI) were computed to test the extent to which genetic mutation affects the TB epidemic using a multivariate logistic regression model.The percentage of the Beijing family strain in hot spots was significantly higher than that in cold spots (64.63% vs 50.69%, P = .022), among the elderly, people with a low BMI, and those having a history of contact with a TB patient (all P < .05). Individuals from cold spot areas had a higher frequency of out-of-town traveling (P < .05). The mutation of Rv1186c, Rv3900c, Rv1508c, Rv0210, and an Intergenic Region (SNP site: 3847237) showed a significant difference between cold and hot spots. (P < .001). The MDR model displayed a clear negative interaction effect of age groups with BMI (interaction entropy: -3.55%) and mutation of Rv0210 (interaction entropy: -2.39%). Through the mutations of Rv0210 and BMI had a low independent effect (interaction entropy: -1.46%).Our data suggests a statistically significant role of age, BMI and the polymorphisms of Rv0210 genes in the transmission and development of M. tb. The results provide clues for the study of susceptibility genes of M. tb in different populations. The characteristic strains showed a local epidemic. Strengthening genotype monitoring of strains in various regions can be used as an early warning signal of epidemic spillover.

Insights into the ancestry evolution of the Mycobacterium tuberculosis complex from analysis of Mycobacterium riyadhense

Current evolutionary scenarios posit the emergence of Mycobacterium tuberculosis from an environmental saprophyte through a cumulative process of genome adaptation. Mycobacterium riyadhense, a related bacillus, is being increasingly isolated from human clinical cases with tuberculosis-like symptoms in various parts of the world. To elucidate the evolutionary relationship between M. riyadhense and other mycobacterial species, including members of the M. tuberculosis complex (MTBC), eight clinical isolates of M. riyadhense were sequenced and analyzed. We show, among other features, that M. riyadhense shares a large number of conserved orthologs with M. tuberculosis and shows the expansion of toxin/antitoxin pairs, PE/PPE family proteins compared with other non-tuberculous mycobacteria. We observed M. riyadhense lacks wecE gene which may result in the absence of lipooligosaccharides (LOS) IV. Comparative transcriptomic analysis of infected macrophages reveals genes encoding inducers of Type I IFN responses, such as cytosolic DNA sensors, were relatively less expressed by macrophages infected with M. riyadhense or M. kansasii, compared to BCG or M. tuberculosis. Overall, our work sheds new light on the evolution of M. riyadhense, its relationship to the MTBC, and its potential as a system for the study of mycobacterial virulence and pathogenesis.

A Comprehensive Map of Mycobacterium tuberculosis Complex Regions of Difference

Mycobacterium tuberculosis complex (MTBC) species are classic examples of genetically monomorphic microorganisms due to their low genetic variability. Whole-genome sequencing made it possible to describe both the main species within the complex and M. tuberculosis lineages and sublineages. This differentiation is based on single nucleotide polymorphisms (SNPs) and large sequence polymorphisms in the so-called regions of difference (RDs). Although a number of studies have been performed to elucidate RD localizations, their distribution among MTBC species, and their role in the bacterial life cycle, there are some inconsistencies and ambiguities in the localization of RDs in different members of the complex. To address this issue, we conducted a thorough search for all possible deletions in the WGS data collection comprising 721 samples representing the full MTBC diversity. Discovered deletions were compared with a list of all previously described RDs. As with the SNP-based analysis, we confirmed the specificities of 79 regions at the species, lineage, or sublineage level, 17 of which are described for the first time. We also present RDscan (https://github.com/dbespiatykh/RDscan), an open-source workflow, which detects deletions from short-read sequencing data and correlates the results with high-specificity RDs, curated in this study. Testing of the workflow on a collection comprising ∼7,000 samples showed a high specificity of the found RDs. This study provides novel details that can contribute to a better understanding of the species differentiation within the MTBC and can help to determine how individual clusters evolve within various MTBC species.

IMPORTANCE Reductive genome evolution is one of the most important and intriguing adaptation strategies of different living organisms to their environment. Mycobacterium offers several notorious examples of either naturally reduced (Mycobacterium leprae) or laboratory-reduced (Mycobacterium bovis BCG) genomes. Mycobacterium tuberculosis complex has its phylogeny unambiguously framed by large sequence polymorphisms that present unidirectional unique event changes. In the present study, we curated all known regions of difference and analyzed both Mycobacterium tuberculosis and animal-adapted MTBC species. For 79 loci, we have shown a relationship with phylogenetic units, which can serve as a marker for diagnosing or studying biological effects. Moreover, intersections were found for some loci, which may indicate the nonrandomness of these processes and the involvement of these regions in the adaptation of bacteria to external conditions.

Pan and Core Genome Analysis of 183 Mycobacterium tuberculosis Strains Revealed a High Inter-Species Diversity among the Human Adapted Strains

Tuberculosis (TB) is an airborne communicable disease with high morbidity and mortality rates, especially in developing countries. The causal agents of TB belong to the complex Mycobacterium tuberculosis (MTBc), which is composed of different human and animal TB associated species. Some animal associated species have zoonotic potential and add to the burden of TB management. The BCG ("Bacillus Calmette-Guérin") vaccine is widely used for the prevention against TB, but its use is limited in immunocompromised patients and animals due to the adverse effects and disseminated life-threatening complications. In this study, we aimed to carry out a comparative genome analysis between the human adapted species including BCG vaccine strains to identify and pinpoint the conserved genes related to the virulence across all the species, which could add a new value for vaccine development. For this purpose, the sequences of 183 Mycobacterium tuberculosis (MTB) strains were retrieved from the freely available WGS dataset at NCBI. The species included: 168 sensu stricto MTB species with other human MTB complex associated strains: M. tuberculosis var. africanum (3), M. tuberculosis var. bovis (2 draft genomes) and 10 BCG species, which enabled the analysis of core genome which contains the conserved genes and some virulence factor determinants. Further, a phylogenetic tree was constructed including the genomes of human (183); animals MTB adapted strains (6) and the environmental Mycobacterium strain "M. canettii". Our results showed that the core genome consists of 1166 conserved genes among these species, which represents a small portion of the pangenome (7036 genes). The remaining genes in the pangenome (5870) are accessory genes, adding a high inter-species diversity. Further, the core genome includes several virulence-associated genes and this could explain the rare infectiousness potential of some attenuated vaccine strains in some patients. This study reveals that low number of conserved genes in human adapted MTBc species and high inter-species diversity of the pan-genome could be considered for vaccine candidate development.

Integrating fish models in tuberculosis vaccine development

Tuberculosis is a chronic infection by Mycobacterium tuberculosis that results in over 1.5 million deaths worldwide each year. Currently, there is only one vaccine against tuberculosis, the Bacillus Calmette–Guérin (BCG) vaccine. Despite widespread vaccination programmes, over 10 million new M. tuberculosis infections are diagnosed yearly, with almost half a million cases caused by antibiotic-resistant strains. Novel vaccination strategies concentrate mainly on replacing BCG or boosting its efficacy and depend on animal models that accurately recapitulate the human disease. However, efforts to produce new vaccines against an M. tuberculosis infection have encountered several challenges, including the complexity of M. tuberculosis pathogenesis and limited knowledge of the protective immune responses. The preclinical evaluation of novel tuberculosis vaccine candidates is also hampered by the lack of an appropriate animal model that could accurately predict the protective effect of vaccines in humans. Here, we review the role of zebrafish (Danio rerio) and other fish models in the development of novel vaccines against tuberculosis and discuss how these models complement the more traditional mammalian models of tuberculosis.

Summary: In this Review, we discuss how zebrafish (Danio rerio) and other fish models can complement the more traditional mammalian models in the development of novel vaccines against tuberculosis.

Immunodominant Mycobacterium tuberculosis Protein Rv1507A Elicits Th1 Response and Modulates Host Macrophage Effector Functions

Mycobacterium tuberculosis (M. tb) persists as latent infection in nearly a quarter of the global population and remains the leading cause of death among infectious diseases. While BCG is the only vaccine for TB, its inability to provide complete protection makes it imperative to engineer BCG such that it expresses immunodominant antigens that can enhance its protective potential. In-silico comparative genomic analysis of Mycobacterium species identified M. tb Rv1507A as a “signature protein” found exclusively in M. tb. In-vitro (cell lines) and in-vivo experiments carried out in mice, using purified recombinant Rv1507A revealed it to be a pro-inflammatory molecule, eliciting significantly high levels of IL-6, TNF-α, and IL-12. There was increased expression of activation markers CD69, CD80, CD86, antigen presentation molecules (MHC I/MHCII), and associated Th1 type of immune response. Rv1507A knocked-in M. smegmatis also induced significantly higher pro-inflammatory Th1 response and higher survivability under stress conditions, both in-vitro (macrophage RAW264.7 cells) and in-vivo (mice). Sera derived from human TB patients showed significantly enhanced B-cell response against M. tb Rv1507A. The ability of M. tb Rv1507A to induce immuno-modulatory effect, B cell response, and significant memory response, renders it a putative vaccine candidate that demands further exploration.

Molecular epidemiological study of multidrug-resistant tuberculosis isolated from sputum samples in Eastern Cape, South Africa

Drug-resistant tuberculosis prevalence is still a global challenge. Making it imperative to examine the molecular epidemiology of drug resistant tuberculosis. Molecular epidemiology methods can evaluate transmission patterns and risk factors, ascertain transmission cases of multidrug-resistant tuberculosis (MDR-TB) and furthermore determine transmission patterns in a human populace. This work focuses on MDR-TB isolates in distinguishing them into several species and genotyping the MDR-TB isolates, mainly for epidemiological studies using the genomic regions of difference and the spoligotyping techniques. A total of 184 deoxyribonucleic acid isolated from sputum samples that showed resistance against the two major first-line anti-tuberculosis drugs (Rifampicin and Isoniazid) were examined. The deoxyribonucleic acid samples were amplified with primers specific for each flanking region of the genomic regions of difference for the identification of different MTBC species. Isolates were further characterized into different lineages using the spoligotyping commercial kit. The M. tuberculosis species was detected in 83.7% (154/184) of the deoxyribonucleic acid isolates, followed by the M. caprae in 8.7% (16/184) and the least detected species was the M. africanum in 2.2% (4/184). Nineteen spoligotype international types (SITs) were identified in this study. The pre-existing shared types were from 94.6% (174/184) isolates with 1.1% (2/184) isolates recognized as orphans and 4.3% (8/184) isolates were not found in the SITVIT database. The predominant family (spoligotype) was the Beijing with 67.4% (124/184) strains. This study gives a general overview of drug resistant strains and the circulating strains in the Eastern Cape, South Africa and it shows that the common Mycobacteria in the province is the Beijing strain.

Matrix Assisted Laser Desorption Ionization-Time-of-Flight mass spectrometry identification of Mycobacterium bovis in Bovinae

In this study, Matrix Assisted Laser Desorption Ionization-Time-of-Flight (MALDI-TOF) mass spectrometry was used to identify Mycobacterium bovis from cattle and buffalo tissue isolates from the North and South regions of Brazil, grown in solid medium and previously identified by Polymerase Chain Reaction (PCR) based on Region of Difference 4 (RD4), sequencing and spoligotyping. For this purpose, the protein extraction protocol and the mass spectra reference database were optimized for the identification of 80 clinical isolates of mycobacteria. As a result of this optimization, it was possible to identify and differentiate M. bovis from other members of the Mycobacterium tuberculosis complex with 100% specificity, 90.91% sensitivity and 91.25% reliability. MALDI-TOF MS methodology described herein provides successful identification of M. bovis within bovine/bubaline clinical samples, demonstrating its usefulness for bovine tuberculosis diagnosis in the future.

Pan-Genomic Study of Mycobacterium tuberculosis Reflecting the Primary/Secondary Genes, Generality/Individuality, and the Interconversion Through Copy Number Variations

Tuberculosis (TB) has surpassed HIV as the leading infectious disease killer worldwide since 2014. The main pathogen, Mycobacterium tuberculosis (Mtb), contains ~4,000 genes that account for ~90% of the genome. However, it is still unclear which of these genes are primary/secondary, which are responsible for generality/individuality, and which interconvert during evolution. Here we utilized a pan-genomic analysis of 36 Mtb genomes to address these questions. We identified 3,679 Mtb core (i.e., primary) genes, determining their phenotypic generality (e.g., virulence, slow growth, dormancy). We also observed 1,122 dispensable and 964 strain-specific secondary genes, reflecting partially shared and lineage-/strain-specific individualities. Among which, five L2 lineage-specific genes might be related to the increased virulence of the L2 lineage. Notably, we discovered 28 Mtb “Super Core Genes” (SCGs: more than a copy in at least 90% strains), which might be of increased importance, and reflected the “super phenotype generality.” Most SCGs encode PE/PPE, virulence factors, antigens, and transposases, and have been verified as playing crucial roles in Mtb pathogenicity. Further investigation of the 28 SCGs demonstrated the interconversion among SCGs, single-copy core, dispensable, and strain-specific genes through copy number variations (CNVs) during evolution; different mutations on different copies highlight the delicate adaptive-evolution regulation amongst Mtb lineages. This reflects that the importance of genes varied through CNVs, which might be driven by selective pressure from environment/host-adaptation. In addition, compared with Mycobacterium bovis (Mbo), Mtb possesses 48 specific single core genes that partially reflect the differences between Mtb and Mbo individuality.