The Nature and Evolution of Genomic Diversity in the Mycobacterium tuberculosis Complex

Prevalence, drug resistance, and genotypic diversity of the RDRio subfamily of Mycobacterium tuberculosis in Ecuador: a retrospective analysis for years 2012-2016

Introduction

A major sublineage within the Mycobacterium tuberculosis (MTB) LAM family characterized by a new in-frame fusion gene Rv3346c/55c was discovered in Rio de Janeiro (Brazil) in 2007, called RDRio, associated to drug resistance. The few studies about prevalence of MTB RDRio strains in Latin America reported values ranging from 3% in Chile to 69.8% in Venezuela, although no information is available for countries like Ecuador.

Methods

A total of 814 MTB isolates from years 2012 to 2016 were screened by multiplex PCR for RDRio identification, followed by 24-loci MIRU-VNTR and spoligotyping.

Results

A total number of 17 MTB RDRio strains were identified, representing an overall prevalence of 2.09% among MTB strains in Ecuador. While 10.9% of the MTB isolates included in the study were multidrug resistance (MDR), 29.4% (5/17) of the RDRio strains were MDR.

Discussion

This is the first report of the prevalence of MTB RDRio in Ecuador, where a strong association with MDR was found, but also a very low prevalence compared to other countries in Latin America. It is important to improve molecular epidemiology tools as a part of MTB surveillance programs in Latin America to track the transmission of potentially dangerous MTB stains associated to MDR TB like MTB RDRio.

Strain structure analysis of Mycobacterium tuberculosis circulating among HIV negative, positive and drug resistant TB patients attending chest clinics in Western Kenya

BACKGROUND

Despite global tuberculosis (TB) interventions, the disease remains one of the major public health concerns. Kenya is ranked 15th among 22 high burden TB countries globally.

METHODS

A cross-sectional study was conducted in Western Kenya, which comprises 10 counties. A multistage sampling method was used where a single sub-county was randomly selected followed by sampling two high volume health facility from each sub-county. Identification of spoligotype profiles and their family distribution and lineage level were achieved by comparison with SITVIT database.

RESULTS

Lineage distribution pattern revealed that the most predominant lineage was CAS 220 (39.8%) followed by Beijing 128 (23.1%). The other lineages identified were T, LAM, H, X, S and MANU which were quantified as 87 (15.7%), 67 (12.1%), 16 (2.8%), 10 (1.8%), 8 (1.4%) and 5 (0.9%) respectively. CAS and Beijing strains were the most predominant lineage in both HIV negative and positive TB patients. The Beijing lineage was also the most predominant in resistant M. tuberculosis strains as compared to wild type. A total of 12 (2.0%) were orphaned M. tuberculosis strains which were spread across all the 10 counties of the study site. In multivariate logistic regression adjusting for potential cofounders three potential risk factors were significant. HIV status (OR = 1.52, CI = 0.29-3.68 and P value of 0.001), Alcohol use (OR = 0.59, CI = 0.43-3.12 and P-value =0.001) and cross border travel (OR = 0.61, CI = 0.49-3.87 and P value = 0.026). Most M. tuberculosis clinical isolates showed genetic clustering with multivariate logistic regression indicating three potential risk factors to clustering. HIV status (OR = 1.52, CI = 0.29-3.68 and P value of 0.001), Alcohol use (OR = 0.59, CI = 0.43-3.12 and P-value =0.001) and cross border travel (OR = 0.61, CI = 0.49-3.87 and P value = 0.026).

CONCLUSION

There exist diverse strains of M. tuberculosis across the 10 counties of Western Kenya. Predominant distribution of clustered genotype points to the fact that most TB cases in this region are as a result of resent transmission other than activation of latent TB.

Pathogenicity and virulence of Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis, an infectious disease with one of the highest morbidity and mortality rates worldwide. Leveraging its highly evolved repertoire of non-protein and protein virulence factors, Mtb invades through the airway, subverts host immunity, establishes its survival niche, and ultimately escapes in the setting of active disease to initiate another round of infection in a naive host. In this review, we will provide a concise synopsis of the infectious life cycle of Mtb and its clinical and epidemiologic significance. We will also take stock of its virulence factors and pathogenic mechanisms that modulate host immunity and facilitate its spread. Developing a greater understanding of the interface between Mtb virulence factors and host defences will enable progress toward improved vaccines and therapeutics to prevent and treat tuberculosis.

Genomic changes underpinning the emergence of a successful Mycobacterium tuberculosis Latin American and Mediterranean clonal complex

Introduction

The Latin American and Mediterranean sublineage (L4.3/LAM) is the most common generalist sublineage of Mycobacterium tuberculosis lineage 4 (L4), yet certain L4.3/LAM genotypes appear to be confined to particular geographic regions. This is typically the case of a L4.3/LAM clonal complex (CC), TUN4.3_CC1, which is the most preponderant in Tunisia (61.5% of L4.3/LAM).

Methods

Here, we used whole-genome sequencing data of 346 globally distributed L4 clinical strains, including 278 L4.3/LAM isolates, to reconstruct the evolutionary history of TUN4.3_CC1 and delineate critical genomic changes underpinning its success.

Results and Discussion

Phylogenomic coupled to phylogeographic analyses indicated that TUN4.3_CC1 has evolved locally, being confined mainly to North Africa. Maximum likelihood analyses using the site and branch-site models of the PAML package disclosed strong evidence of positive selection in the gene category "cell wall and cell processes" of TUN4.3_CC1. Collectively, the data indicate that TUN4.3_CC1 has inherited several mutations, which could have potentially contributed to its evolutionary success. Of particular interest are amino acid replacements at the esxK and eccC2 genes of the ESX/Type VII secretion system, which were found to be specific to TUN4.3_CC1, being common to almost all isolates. Because of its homoplastic nature, the esxK mutation could potentially have endowed TUN4.3_CC1 with a selective advantage. Moreover, we noticed the occurrence of additional, previously described homoplasic nonsense mutations in ponA1 and Rv0197. The mutation in the latter gene, a putative oxido-reductase, has previously been shown to be correlated with enhanced transmissibility in vivo. In sum, our findings unveiled several features underpinning the success of a locally evolved L4.3/LAM clonal complex, lending further support to the critical role of genes encoded by the ESX/type VII secretion system.

Drug Resistance in Tuberculous Lymphadenitis: Molecular Characterization

Background

Drug-resistant tuberculosis (TB) epidemic in high-TB-incidence countries, particularly Ethiopia, remains a significant challenge. As a result, we investigated the drug resistance, common gene mutation, and molecular characterization of mycobacterial isolates from patients with suspected tuberculous lymphadenitis (TBLN). Methodology. A cross-sectional study of 218 FNA samples from TBLN patients inoculated on Lowenstein-Jensen media was carried out. The culture isolates were identified as MTB by polymerase chain reaction (PCR) and the difference-9 (RD9) test region. In addition, the GenoType MTBDRplus assay tested the first and second-line MTB drugs, and the spoligotyping strain-dependent polymorphism test was determined.

Results

Among the 50 culture-positive isolates, 14% (7/50) had drug resistance caused by a gene mutation. Out of these, 4 (8%) isolates were mono-resistant to isoniazid drug, which is caused by a gene mutation in katG in the region of interrogated at codon 315 in the amino acid sequence of S315T1, and 3 (6%) isolates were resistant to both rifampicin and isoniazid drugs. The mutation was observed for katG (at codon 315 with a change in the sequence of amino acid S315T) and rpoB (at codon 530-533 with a change in the sequence of amino acid S531L (S450L)) genes. The most prevalent spoligotypes were orphan and SIT53 strains.

Conclusion

The predominance of INH mono-resistance poses a critical risk for the potential development of MDR-TB, as INH mono-resistance is a typical pathway to the occurrence of MDR-TB. The orphan and SIT53 (T) strains were the most common in the study area, and a drug-resistant strain caused by a common gene mutation could indicate the transmission of clonal-resistant strains in the community.

Performance of lipid fingerprint by routine matrix-assisted laser desorption/ionization time of flight for the diagnosis of Mycobacterium tuberculosis complex species

OBJECTIVES

Rapid detection of bacterial pathogens at species and sub-species levels is crucial for appropriate treatment, infection control, and public health management. Currently, one of the challenges in clinical microbiology is the discrimination of mycobacterial sub-species within the M. tuberculosis complex (MTBC). Our objective was to evaluate the ability of a biosafe mycobacterial lipid-based approach to identify MTBC cultures and sub-species.

METHODS

A blinded study was conducted using 90 mycobacterial clinical isolate strains comprising MTBC strains sub-cultured in Middlebrook 7H11 medium supplemented with 10% oleic-acid, dextrose, catalase growth supplement and incubated for up to 6 weeks at 37°C and using the following seven reference strains (M. tuberculosis H37Rv, M canettii, M. africanum, M. pinnipedii, M. caprae, M. bovis, and M. bovis BCG) grown under the same conditions, to set the reference lipid database and test it against the 90 MTBC clinical isolates. Cultured mycobacteria were heat-inactivated and loaded onto the matrix-assisted laser desorption/ionization target followed by the addition of the matrix. Acquisition of the data was performed using the positive ion mode.

RESULTS

Based on the identification of clear and defined lipid signatures from the seven reference strains, the method that we developed was fast (<10 minutes) and produced interpretable profiles for all but four isolates, caused by poor ionization giving an n = 86 with interpretable spectra. The sensitivity and specificity of the matrix-assisted laser desorption/ionization time of flight were 94.4 (95% CI, 86.4-98.5) and 94.4 (95% CI, 72.7-99.9), respectively.

CONCLUSIONS

Mycobacterial lipid profiling provides a means of rapid, safe, and accurate discrimination of species within the MTBC.

Evolutionary history and spread of the Mycobacterium tuberculosis Latin American and Mediterranean (L4.3/LAM) sublineage, Tunisia

BACKGROUND

To infer the origin and spread of the Mycobacterium tuberculosis Latin American and Mediterranean (L4.3/LAM) sublineage in a Mediterranean country, Tunisia, where it predominates.

METHODS

We combined Bayesian (STRUCTURE) and maximum likelihood (MIGRAINE) estimation approaches based on a global 24-loci mycobacterial interspersed repetitive units-variable numbers of tandem repeats (MIRU-VNTR24) genotyping dataset consisting of 1573 L4.3/LAM clinical strains from four continents, including 252 isolates originating from Tunisia.

RESULTS

Phylogenetic analyses coupled with Bayesian estimations suggested that the most predominant L4.3/LAM subpopulation in Tunisia (65.07%), which is dominated by a single clonal complex, TUN4.3_CC1 (94.51%), has evolved from an ancestral pool that is restricted to Europe and Africa, contrasting with the remaining L4.3/LAM subpopulations whose ancestry was traced all over the word. Maximum likelihood analysis revealed that TUN4.3_CC1 has been undergoing a demographic expansion since 131 years ago (CI95%: 90.7-205), thus explaining its preponderance relative to the second most predominant CC, TUN4.3_CC2, whose population was found under contraction.

CONCLUSIONS

The preponderance of L4.3/LAM in Tunisia stems from a 130-year expansion process of a locally evolved clone.

Human Pulmonary Tuberculosis: Understanding the Immune Response in the Bronchoalveolar System

Mycobacterium tuberculosis, the causal agent of one of the most devastating infectious diseases worldwide, can evade or modulate the host immune response and remain dormant for many years. In this review, we focus on identifying the local immune response induced in vivo by M. tuberculosis in the lungs of patients with active tuberculosis by analyzing data from untouched cells from bronchoalveolar lavage fluid (BALF) or exhaled breath condensate (EBC) samples. The most abundant resident cells in patients with active tuberculosis are macrophages and lymphocytes, which facilitate the recruitment of neutrophils. The cellular response is characterized by an inflammatory state and oxidative stress produced mainly by macrophages and T lymphocytes. In the alveolar microenvironment, the levels of cytokines such as interleukins (IL), chemokines, and matrix metalloproteinases (MMP) are increased compared with healthy patients. The production of cytokines such as interferon (IFN)-γ and IL-17 and specific immunoglobulin (Ig) A and G against M. tuberculosis indicate that the adaptive immune response is induced despite the presence of a chronic infection. The role of epithelial cells, the processing and presentation of antigens by macrophages and dendritic cells, as well as the role of tissue-resident memory T cells (Trm) for in situ vaccination remains to be understood.

Bacterial Strain-Dependent Dissociation of Cell Recruitment and Cell-to-Cell Spread in Early M. tuberculosis Infection

In the initial stage of respiratory infection, Mycobacterium tuberculosis traverses from alveolar macrophages to phenotypically diverse monocyte-derived phagocytes and neutrophils in the lung parenchyma. Here, we compare the in vivo kinetics of early bacterial growth and cell-to-cell spread of two strains of M. tuberculosis: a lineage 2 strain, 4334, and the widely studied lineage 4 strain H37Rv. Using flow cytometry, live cell sorting of phenotypic subsets, and quantitation of bacteria in cells of the distinct subsets, we found that 4334 induces less leukocyte influx into the lungs but demonstrates earlier population expansion and cell-to-cell spread. The earlier spread of 4334 to recruited cells, including monocyte-derived dendritic cells, is accompanied by earlier and greater magnitude of CD4+ T cell activation. The results provide evidence that strain-specific differences in interactions with lung leukocytes can shape adaptive immune responses in vivo. IMPORTANCE Tuberculosis is a leading infectious disease killer worldwide and is caused by Mycobacterium tuberculosis. After exposure to M. tuberculosis, outcomes range from apparent elimination to active disease. Early innate immune responses may contribute to differences in outcomes, yet it is not known how bacterial strains alter the early dynamics of innate immune and T cell responses. We infected mice with distinct strains of M. tuberculosis and discovered striking differences in innate cellular recruitment, cell-to-cell spread of bacteria in the lungs, and kinetics of initiation of antigen-specific CD4 T cell responses. We also found that M. tuberculosis can spread beyond alveolar macrophages even before a large influx of inflammatory cells. These results provide evidence that distinct strains of M. tuberculosis can exhibit differential kinetics in cell-to-cell spread which is not directly linked to early recruitment of phagocytes but is subsequently linked to adaptive immune responses.

Sub-Lineage Specific Phenolic Glycolipid Patterns in the Mycobacterium tuberculosis Complex Lineage 1

“Ancestral” Mycobacterium tuberculosis complex (MTBC) strains of Lineage 1 (L1, East African Indian) are a prominent tuberculosis (TB) cause in countries around the Indian Ocean. However, the pathobiology of L1 strains is insufficiently characterized. Here, we used whole genome sequencing (WGS) of 312 L1 strains from 43 countries to perform a characterization of the global L1 population structure and correlate this to the analysis of the synthesis of phenolic glycolipids (PGL) – known MTBC polyketide-derived virulence factors. Our results reveal the presence of eight major L1 sub-lineages, whose members have specific mutation signatures in PGL biosynthesis genes, e.g., pks15/1 or glycosyltransferases Rv2962c and/or Rv2958c. Sub-lineage specific PGL production was studied by NMR-based lipid profiling and strains with a completely abolished phenolphthiocerol dimycoserosate biosynthesis showed in average a more prominent growth in human macrophages. In conclusion, our results show a diverse population structure of L1 strains that is associated with the presence of specific PGL types. This includes the occurrence of mycoside B in one sub-lineage, representing the first description of a PGL in an M. tuberculosis lineage other than L2. Such differences may be important for the evolution of L1 strains, e.g., allowing adaption to different human populations.

Extracellular Matrix-Induced GM-CSF and Hypoxia Promote Immune Control of Mycobacterium tuberculosis in Human In Vitro Granulomas

Several in vitro cellular models have been developed with the aim to reproduce and dissect human granulomatous responses, the hallmark of tuberculosis (TB) immunopathogenesis. In that context, we compared two- (2D) versus three-dimensional (3D) granuloma models resulting from infection of human peripheral blood mononuclear cells with M. tuberculosis (Mtb) in the absence or presence of a collagen-based extracellular matrix (ECM). Granuloma formation was found to be significantly enhanced in the 2D model. This feature was associated with an earlier chemokine production and lymphocyte activation, but also a significantly increased bacterial burden. Remarkably, the reduction in Mtb burden in the 3D model correlated with an increase in GM-CSF production. GM-CSF, which is known to promote macrophage survival, was found to be inherently induced by the ECM. We observed that only 3D in vitro granulomas led to the accumulation of lipid inclusions within Mtb. Our data suggest that a hypoxic environment within the ECM could be responsible for this dormant-like Mtb phenotype. Furthermore, exposure to a TNF-α antagonist reverted Mtb dormancy, thereby mimicking the reactivation of TB observed in rheumatic patients receiving this therapy. To conclude, we showed that only in vitro granulomas generated in the presence of an ECM could recapitulate some clinically relevant features of granulomatous responses in TB. As such, this model constitutes a highly valuable tool to study the interplay between immunity and Mtb stress responses as well as to evaluate novel treatment strategies.

Genome-wide estimation of recombination, mutation and positive selection enlightens diversification drivers of Mycobacterium bovis

Genome sequencing has reinvigorated the infectious disease research field, shedding light on disease epidemiology, pathogenesis, host-pathogen interactions and also evolutionary processes exerted upon pathogens. Mycobacterium tuberculosis complex (MTBC), enclosing M. bovis as one of its animal-adapted members causing tuberculosis (TB) in terrestrial mammals, is a paradigmatic model of bacterial evolution. As other MTBC members, M. bovis is postulated as a strictly clonal, slowly evolving pathogen, with apparently no signs of recombination or horizontal gene transfer. In this work, we applied comparative genomics to a whole genome sequence (WGS) dataset composed by 70 M. bovis from different lineages (European and African) to gain insights into the evolutionary forces that shape genetic diversification in M. bovis. Three distinct approaches were used to estimate signs of recombination. Globally, a small number of recombinant events was identified and confirmed by two independent methods with solid support. Still, recombination reveals a weaker effect on M. bovis diversity compared with mutation (overall r/m = 0.037). The differential r/m average values obtained across the clonal complexes of M. bovis in our dataset are consistent with the general notion that the extent of recombination may vary widely among lineages assigned to the same taxonomical species. Based on this work, recombination in M. bovis cannot be excluded and should thus be a topic of further effort in future comparative genomics studies for which WGS of large datasets from different epidemiological scenarios across the world is crucial. A smaller M. bovis dataset (n = 42) from a multi-host TB endemic scenario was then subjected to additional analyses, with the identification of more than 1,800 sites wherein at least one strain showed a single nucleotide polymorphism (SNP). The majority (87.1%) was located in coding regions, with the global ratio of non-synonymous upon synonymous alterations (dN/dS) exceeding 1.5, suggesting that positive selection is an important evolutionary force exerted upon M. bovis. A higher percentage of SNPs was detected in genes enriched into "lipid metabolism", "cell wall and cell processes" and "intermediary metabolism and respiration" functional categories, revealing their underlying importance in M. bovis biology and evolution. A closer look on genes prone to horizontal gene transfer in the MTBC ancestor and included in the 3R (DNA repair, replication and recombination) system revealed a global average negative value for Taijima's D neutrality test, suggesting that past selective sweeps and population expansion after a recent bottleneck remain as major evolutionary drivers of the obligatory pathogen M. bovis in its struggle with the host.

Scientific considerations for global drug development

Requiring regional or in-country confirmatory clinical trials before approval of drugs already approved elsewhere delays access to medicines in low- and middle-income countries and raises drug costs. Here, we discuss the scientific and technological advances that may reduce the need for in-country or in-region clinical trials for drugs approved in other countries and limitations of these advances that could necessitate in-region clinical studies.

Host Genetic Diversity and Infectious Diseases. Focus on Wild Boar, Red Deer and Tuberculosis

Host genetic diversity tends to limit disease spread in nature and buffers populations against epidemics. Genetic diversity in wildlife is expected to receive increasing attention in contexts related to disease transmission and human health. Ungulates such as wild boar (Sus scrofa) and red deer (Cervus elaphus) are important zoonotic hosts that can be precursors to disease emergence and spread in humans. Tuberculosis is a zoonotic disease with relevant consequences and can present high prevalence in wild boar and red deer populations. Here, we review studies on the genetic diversity of ungulates and determine to what extent these studies consider its importance on the spread of disease. This assessment also focused on wild boar, red deer, and tuberculosis. We found a disconnection between studies treating genetic diversity and those dealing with infectious diseases. Contrarily, genetic diversity studies in ungulates are mainly concerned with conservation. Despite the existing disconnection between studies on genetic diversity and studies on disease emergence and spread, the knowledge gathered in each discipline can be applied to the other. The bidirectional applications are illustrated in wild boar and red deer populations from Spain, where TB is an important threat for wildlife, livestock, and humans.

The puzzle of the evolutionary natural history of tuberculosis

Several pieces of the puzzle of the natural history of tuberculosis are assembled in this review to illustrate the potential reservoirs and sources of the Mycobacterium tuberculosis complex (MTBC) mycobacteria, their transmission to animals and humans, and their fate in populations, in a co-evolutionary perspective. Millennia-old companions of mammalian and human populations, MTBC are detected in the soil, in which they infect and survive within vegetative amoebae and cysts, except for Mycobacterium canettii. Never detected in the sphere of plants, they are transmissible by transcutaneous, digestive and respiratory routes and cause an infection of the lymphatic system with secondary dissemination in most tissues, in which they determine a specific and non-pathognomonic granulomatous inflammatory reaction; in which MTBC survives in dormant form irrespective of MTBC species and mammalian species; indicating that the current epidemiology in mammalian populations is essentially governed by the probabilities of contact between mammalian species and MTBC species. Individual variabilities in clinical expression of tuberculosis are related to MTBC species, strain and inoculum; host genetic factors; acquired modulations of the inflammatory response; and probably human microbiota. This review of the literature suggests an evolutionary natural history of telluric environmental mycobacteria, satellites of unicellular eukaryotes, transmissible to mammals via the digestive and then respiratory tracts, in which they determine a fatal contagious infection that is primarily lymphatic and a quiescence-mimicking encysted form. This review opens perspectives for microbiological and translational medical research.

One Size Fits All? Not in In Vivo Modeling of Tuberculosis Chemotherapeutics

Tuberculosis (TB) remains a global health problem despite almost universal efforts to provide patients with highly effective chemotherapy, in part, because many infected individuals are not diagnosed and treated, others do not complete treatment, and a small proportion harbor Mycobacterium tuberculosis (Mtb) strains that have become resistant to drugs in the standard regimen. Development and approval of new drugs for TB have accelerated in the last 10 years, but more drugs are needed due to both Mtb's development of resistance and the desire to shorten therapy to 4 months or less. The drug development process needs predictive animal models that recapitulate the complex pathology and bacterial burden distribution of human disease. The human host response to pulmonary infection with Mtb is granulomatous inflammation usually resulting in contained lesions and limited bacterial replication. In those who develop progressive or active disease, regions of necrosis and cavitation can develop leading to lasting lung damage and possible death. This review describes the major vertebrate animal models used in evaluating compound activity against Mtb and the disease presentation that develops. Each of the models, including the zebrafish, various mice, guinea pigs, rabbits, and non-human primates provides data on number of Mtb bacteria and pathology resolution. The models where individual lesions can be dissected from the tissue or sampled can also provide data on lesion-specific bacterial loads and lesion-specific drug concentrations. With the inclusion of medical imaging, a compound's effect on resolution of pathology within individual lesions and animals can also be determined over time. Incorporation of measurement of drug exposure and drug distribution within animals and their tissues is important for choosing the best compounds to push toward the clinic and to the development of better regimens. We review the practical aspects of each model and the advantages and limitations of each in order to promote choosing a rational combination of them for a compound's development.

Drug Resistance Conferring Mutation and Genetic Diversity of Mycobacterium tuberculosis Isolates in Tuberculosis Lymphadenitis Patients; Ethiopia

Background

Tuberculosis lymphadenitis (TBLN) is a growing public health concern in Ethiopia. However, there is limited information available on gene mutations conferring drug resistance and genetic diversity of M. tuberculosis isolates from TBLN patients.

Methods

Drug resistance and genetic diversity analysis were done on 91 M. tuberculosis isolates from culture positive TBLN patients collected between 2016 and 2017. Detection of mutations conferring resistance was carried out using GenoType MTBDRplus VER 2.0. Thereafter, isolates were typed using spoligotyping.

Results

Out of the 91 strains, mutations conferring resistance to rifampicin (RIF) and isoniazid (INH) were observed in two (2.2%) and six (6.6%) isolates, respectively. The two RIF resistant isolates displayed a mutation at codon 531 in the rpoB gene with amino acid change of S531L. Among the six INH resistant strains, four isolates had shown mutation at the KatG gene at codon 315 with amino acid change of S315T, one isolate had a mutation at the inhA gene at codon 15 with amino acid change of C15T and one isolate had a mutation at the inhA gene with unknown amino acid change. All drug resistant isolates were from treatment naive TBLN patients. The dominantly identified Spoligo International Types (SITs) were SIT25, SIT149, and SIT53, respectively; these accounted for 43% of the total number of strains. The isolates were grouped into four main lineages; Lineage 1 (2, 2.2%), Lineage 3 (38, 41.7%), Lineage 4 (49, 53.8%) and Lineage 7 (2, 2.2%). Four out of six (66.7%) isolates with drug resistance conferring mutations belonged to clustered strains (strains with shared SIT).

Conclusion

The detection of drug resistant conferring mutation in treatment naïve TBLN patients together with detection of drug resistant isolates among clustered strains might suggest resistant strains' transmission in the community. This needs to be carefully considered to prevent the spread of drug resistant clones in the country.

Variation in Mycobacterium tuberculosis population structure in Iran: a systemic review and meta-analysis

Background

Acquiring comprehensive insight into the dynamics of Mycobacterium tuberculosis (Mtb) population structure is an essential step to adopt effective tuberculosis (TB) control strategies and improve therapeutic methods and vaccines. Accordingly, we performed this systematic review and meta-analysis to determine the overall prevalence of Mtb genotypes/ sublineages in Iran.

Methods

We carried out a comprehensive literature search using the international databases of MEDLINE and Scopus as well as Iranian databases. Articles published until April 2020 were selected based on the PRISMA flow diagram. The overall prevalence of the Mtb genotypes/sublineage in Iran was determined using the random effects or fixed effect model. The metafor R package and MedCalc software were employed for performing this meta-analysis.

Results

We identified 34 studies for inclusion in this study, containing 8329 clinical samples. Based on the pooled prevalence of the Mtb genotypes, NEW1 (21.94, 95% CI: 16.41–28.05%), CAS (19.21, 95% CI: 14.95–23.86%), EAI (12.95, 95% CI: 7.58–19.47%), and T (12.16, 95% CI: 9.18–15.50%) were characterized as the dominant circulating genotypes in Iran. West African (L 5/6), Cameroon, TUR and H37Rv were identified as genotypes with the lowest prevalence in Iran (< 2%). The highest pooled prevalence rates of multidrug-resistant strains were related to Beijing (2.52, 95% CI) and CAS (1.21, 95% CI).

Conclusions

This systematic review showed that Mtb populations are genetically diverse in Iran, and further studies are needed to gain a better insight into the national diversity of Mtb populations and their drug resistance pattern.

Independent genomic polymorphisms in the PknH serine threonine kinase locus during evolution of the Mycobacterium tuberculosis Complex affect virulence and host preference

Species belonging to the Mycobacterium tuberculosis Complex (MTBC) show more than 99% genetic identity but exhibit distinct host preference and virulence. The molecular genetic changes that underly host specificity and infection phenotype within MTBC members have not been fully elucidated. Here, we analysed RD900 genomic region across MTBC members using whole genome sequences from 60 different MTBC strains so as to determine its role in the context of MTBC evolutionary history. The RD900 region comprises two homologous genes, pknH1 and pknH2, encoding a serine/threonine protein kinase PknH flanking the tbd2 gene. Our analysis revealed that RD900 has been independently lost in different MTBC lineages and different strains, resulting in the generation of a single pknH gene. Importantly, all the analysed M. bovis and M. caprae strains carry a conserved deletion within a proline rich-region of pknH, independent of the presence or absence of RD900. We hypothesized that deletion of pknH proline rich-region in M. bovis may affect PknH function, having a potential role in its virulence and evolutionary adaptation. To explore this hypothesis, we constructed two M. bovis ‘knock-in’ strains containing the M. tuberculosis pknH gene. Evaluation of their virulence phenotype in mice revealed a reduced virulence of both M. bovis knock-in strains compared to the wild type, suggesting that PknH plays an important role in the differential virulence phenotype of M. bovis vs M. tuberculosis.

Tuberculosis is caused in humans and animals by organisms from the Mycobacterium tuberculosis Complex (MTBC), that share more than 99% genetic identity but exhibit distinct host preference and virulence. While Mycobacterium tuberculosis is the main causative agent of human TB, Mycobacterium bovis is responsible for bovine TB disease, that exacts a tremendous economic burden worldwide, as well as being a zoonotic threat. Unlike the human restriction of M. tuberculosis, M. bovis has a broader host range and it has been found to be more virulent than M. tuberculosis in different animal models. However, the molecular basis for host preference and virulence divergence between M. tuberculosis and M. bovis is not fully elucidated. Here we study the genetic variations of the genomic region RD900 in the context of MTBC phylogeny. RD900 contains two genes encoding orthologues of the serine/threonine kinase PknH, which is linked to the regulation of several bacterial processes including virulence. We found that M. bovis pknH genes show a conserved deletion that is not present in M. tuberculosis strains, and we evaluated the potential impact of these variations in the regulation of M. bovis vs M. tuberculosis virulence through the construction and in vivo characterization of M. bovis pknH mutant strains.

Successful expansion of Mycobacterium tuberculosis Latin American and Mediterranean sublineage (L4.3/LAM) in Tunisia mainly driven by a single, long-established clonal complex

OBJECTIVES

To explore the evolutionary history of Mycobacterium tuberculosis Latin American and Mediterranean (L4.3/LAM) sublineage in Tunisia, where it predominates.

METHODS

High-resolution genotyping of 252 L4.3/LAM clinical strains was undertaken, and whole-genome sequencing was performed on 31 representative isolates.

RESULTS

Genotyping data coupled with Bayesian analyses split the Tunisian L4.3/LAM strain collection into two divergent entities (65.07% vs 34.92%): a major subpopulation, dominated by a single clonal complex (CC), TUN4.3_CC1 (94.51%); and a minor subpopulation, dominated by TUN4.3_CC2 (42.04%). TUN4.3_CC1 is clearly thriving in Tunisia, accounting for 61.5% of the L4.3/LAM sublineage. TUN4.3_CC1 displayed higher mean allelic richness compared with TUN4.3_CC2 and predominated throughout the entire region, indicating a long-established history. The very low proportion of drug resistance among TUN4.3_CC1 isolates is indicative of their intrinsic ability to spread successfully in the host population. Genomic analyses further confirmed the clear genetic separation between the two main CCs (pairwise fixation index 0.56), and suggested the relatively ancient origin of TUN4.3_CC1. Consistent with its successful expansion, TUN4.3_CC1 showed reduced mean pairwise genetic distance between genomes.

CONCLUSIONS

These findings link the successful expansion of L4.3/LAM in Tunisia to a single long-established clone.

Mycobacterium tuberculosis complex and human coadaptation: a two-way street complicating host susceptibility to TB

For centuries, the Mycobacterium tuberculosis complex (MTBC) has infected numerous populations, both human and non-human, causing symptomatic tuberculosis (TB) in some hosts. Research investigating the MTBC and how it has evolved with its host over time is sparse and has not resulted in many significant findings. There are even fewer studies investigating adaptation of the human host susceptibility to TB and these have largely focused on genome-wide association and candidate gene association studies. However, results emanating from these association studies are rarely replicated and appear to be population specific. It is, therefore, necessary to relook at the approach taken to investigate the relationship between the MTBC and the human host. Understanding that the evolution of the pathogen is coupled to the evolution of the host might be the missing link needed to effectively investigate their relationship. We hypothesize that this knowledge will bolster future efforts in combating the disease.

The whole genome sequence data analyses of a Mycobacterium tuberculosis strain SBH321 isolated in Sabah, Malaysia, belongs to Ural family of Lineage 4

In 2019, 10 million new cases of tuberculosis have been reported worldwide. Our data reports genetic analyses of a Mycobacterium tuberculosis strain SBH321 isolated from a 31-year-old female with pulmonary tuberculosis. The genomic DNA of the strain was extracted from pure culture and subjected to sequencing using Illumina platform. M. tuberculosis strain SBH321 consists of 4,374,895 bp with G+C content of 65.59%. The comparative analysis by SNP-based phylogenetic analysis using maximum-likelihood method showed that our strain belonging to sublineage of the Ural family of Europe-America-Africa lineage (Lineage 4) and clustered with M. tuberculosis strain OFXR-4 from Taiwan. The whole genome sequence is deposited at DDBJ/ENA/GenBank under the accession WCJH00000000 (SRR10230353).

Multicenter analysis of sputum microbiota in tuberculosis patients

The impact of tuberculosis and of anti-tuberculosis therapy on composition and modification of human lung microbiota has been the object of several investigations. However, no clear outcome has been presented so far and the relationship between M. tuberculosis pulmonary infection and the resident lung microbiota remains vague. In this work we describe the results obtained from a multicenter study of the microbiota of sputum samples from patients with tuberculosis or unrelated lung diseases and healthy donors recruited in Switzerland, Italy and Bangladesh, with the ultimate goal of discovering a microbiota-based biomarker associated with tuberculosis. Bacterial 16S rDNA amplification, high-throughput sequencing and extensive bioinformatic analyses revealed patient-specific flora and high variability in taxon abundance. No common signature could be identified among the individuals enrolled except for minor differences which were not consistent among the different geographical settings. Moreover, anti-tuberculosis therapy did not cause any important variation in microbiota diversity, thus precluding its exploitation as a biomarker for the follow up of tuberculosis patients undergoing treatment.

Culturomics Discloses Anti-Tubercular Enterococci Exclusive of Pulmonary Tuberculosis: A Preliminary Report

Mycobacterium tuberculosis causes pulmonary tuberculosis, a deadly infection of which the clinical expression and prognosis are not fully understood at the individual level, apart from genetic susceptibility traits. We investigated whether individual gut microbiota may correlate with pulmonary tuberculosis status. Culturomics investigations of gut microbiota in two pulmonary tuberculosis patients and two controls in Burkina Faso found 60 different bacterial species in patients and 97 in controls, including 45 in common. Further analysis of the results at the individual level indicated seven bacteria, including Enterococcus mundtii and Enterococcus casseliflavus, which were exclusively cultured in controls. Blind quantitative PCR-based exploration of faeces samples in two cohorts in Burkina Faso and in France confirmed a nonsignificant association of E. mundtii and E. casseliflavus with controls. Further in vitro explorations found four E. mundtii and E. casseliflavus strains inhibiting the growth of M. tuberculosis strains representative of four different lineages as well as Mycobacterium africanum, Mycobacterium canettii, and Mycobacterium bovis, in an inoculum-dependent manner. Heat-killed E. mundtii or E. casseliflavus were ineffective. These unprecedented observations of direct interactions between gut E. mundtii and E. casseliflavus with M. tuberculosis complex mycobacteria suggest that gut microbiota may modulate the expression of pulmonary tuberculosis.

Prolonged infection triggered by dormant Mycobacterium tuberculosis: Immune and inflammatory responses in lungs of genetically susceptible and resistant mice

We developed an approach for substantial attenuation of Mycobacterium tuberculosis by prolonged culturing under gradually acidifying conditions. Bacteria subjected to acidification lost the capacity to form colonies on solid media, but readily resuscitated their growth in the murine host, providing a useful model to study in vivo development of infection mimicking latent and reactivation tuberculosis (TB) in humans. Here we characterize biomarkers of lung pathology and immune responses triggered by such attenuated bacteria in genetically TB-susceptible and resistant mice. In susceptible I/St mice, CFU counts in lungs and spleens were ~1.5-log higher than in resistant B6 mice, accompanied by diffuse pneumonia and excessive lung infiltration with highly activated CD44⁺CD62L^- T-lymphocytes resulting in death between months 7–9 post challenge. B6 mice were characterized by development of local inflammatory foci, higher production of pro-inflammatory IL-6 and IL-11 cytokines and a more balanced T-cell activation in their lungs. CFU counts remained stable in B6 mice during the whole 18-mo observation period, and all mice survived. Thus, we established a mouse model of fatal reactivation TB vs. indefinite mycobacterial possession after identical challenge and characterized the features of immune responses in the lung tissue underlining these polar phenotypes.

Global Distribution and Evolution of Mycobacterium bovis Lineages

Mycobacterium bovis is the main causative agent of zoonotic tuberculosis in humans and frequently devastates livestock and wildlife worldwide. Previous studies suggested the existence of genetic groups of M. bovis strains based on limited DNA markers (a.k.a. clonal complexes), and the evolution and ecology of this pathogen has been only marginally explored at the global level. We have screened over 2,600 publicly available M. bovis genomes and newly sequenced four wildlife M. bovis strains, gathering 1,969 genomes from 23 countries and at least 24 host species, including humans, to complete a phylogenomic analyses. We propose the existence of four distinct global lineages of M. bovis (Lb1, Lb2, Lb3, and Lb4) underlying the current disease distribution. These lineages are not fully represented by clonal complexes and are dispersed based on geographic location rather than host species. Our data divergence analysis agreed with previous studies reporting independent archeological data of ancient M. bovis (South Siberian infected skeletons at ∼2,000 years before present) and indicates that extant M. bovis originated between 715 and 3,556 years BP, with later emergence in the New World and Oceania, likely influenced by trades among countries.

Live attenuated TB vaccines representing the three modern Mycobacterium tuberculosis lineages reveal that the Euro-American genetic background confers optimal vaccine potential

Background

Human tuberculosis (TB) is caused by a plethora of Mycobacterium tuberculosis complex (MTBC) strains belonging to seven phylogenetic branches. Lineages 2, 3 and 4 are considered “modern” branches of the MTBC responsible for the majority of worldwide TB. Since the current BCG vaccine confers variable protection against pulmonary TB, new candidates are investigated. MTBVAC is the unique live attenuated vaccine based on M. tuberculosis in human clinical trials.

Methods

MTBVAC was originally constructed by unmarked phoP and fadD26 deletions in a clinical isolate belonging to L4. Here we construct new vaccines based on isogenic gene deletions in clinical isolates of the L2 and L3 modern lineages. These three vaccine candidates were characterized at molecular level and also in animal experiments of protection and safety.

Findings

Safety studies in immunocompromised mice showed that MTBVAC-L2 was less attenuated than BCG Pasteur, while the original MTBVAC was found even more attenuated than BCG and MTBVAC-L3 showed an intermediate phenotype. The three MTBVAC candidates showed similar or superior protection compared to BCG in immunocompetent mice vaccinated with each MTBVAC candidate and challenged with three representative strains of the modern lineages.

Interpretation

MTBVAC vaccines, based on double phoP and fadD26 deletions, protect against TB independently of the phylogenetic linage used as template strain for their construction. Nevertheless, lineage L4 confers the best safety profile.

Funding

European Commission (TBVAC2020, H2020-PHC-643381), Spanish Ministry of Science (RTI2018-097625-B-I00), Instituto de Salud Carlos III (PI18/0336), Gobierno de Aragón/Fondo Social Europeo and the French National Research Council (ANR-10-LABX-62-IBEID, ANR-16-CE35-0009, ANR-16-CE15-0003).

The Interplay of Human and Mycobacterium Tuberculosis Genomic Variability

Tuberculosis (TB), caused by the human pathogens Mycobacterium tuberculosis (Mtb) and Mycobacterium africanum, has plagued humanity for millennia and remains the deadliest infectious disease in the modern world. Mycobacterium tuberculosis and M. africanum can be subdivided phylogenetically into seven lineages exhibiting a low but significant degree of genomic diversity and preferential geographic distributions. Human genetic variability impacts all stages of TB pathogenesis ranging from susceptibility to infection with Mtb, progression of infection to disease, and the development of distinct clinical subtypes. The genetic study of severe childhood TB identified strong inborn single-gene errors revealing crucial pathways of vulnerability to TB. However, the identification of major TB-susceptibility genes on the population level has remained elusive. In particular, the replication of findings from candidate and genome-wide association studies across distinct human populations has proven difficult, thus hampering the characterization of reliable host molecular markers of susceptibility. Among the possible confounding factors of genetic association studies is Mtb genomic variability, which generally was not taken into account by human genetic studies. In support of this possibility, Mtb lineage was found to be a contributing factor to clinical presentation of TB and epidemiological spread of Mtb in exposed populations. The confluence of pathogen and human host genetic variability to TB pathogenesis led to the consideration of a possible coadaptation of Mtb strains and their human hosts, which should reveal itself in significant interaction effects between Mtb strain and TB-susceptibility/resistance alleles. Here, we present some of the most consistent findings of genetic susceptibility factors in human TB and review studies that point to genome-to-genome interaction between humans and Mtb lineages. The limited results available so far suggest that analyses considering joint human–Mtb genomic variability may provide improved power for the discovery of pathogenic drivers of the ongoing TB epidemic.

Bacterial Genotyping of Central Nervous System Tuberculosis in South Africa: Heterogenic Mycobacterium tuberculosis Infection and Predominance of Lineage 4

Tuberculous meningitis (TBM), the most severe extrapulmonary manifestation of tuberculosis, is caused by the pathogen Mycobacterium tuberculosis The M. tuberculosis complex includes seven lineages, all described to harbor a unique geographical dissemination pattern and clinical presentation. In this study, we set out to determine whether a certain M. tuberculosis lineage demonstrated tropism to cause TBM in patients from Cape Town, South Africa. DNA was extracted from formalin-fixed paraffin-embedded central nervous system (CNS) tissue from a unique neuropathological cohort of 83 TBM patients, collected between 1975 and 2012. M. tuberculosis lineages 1, 2, 3, and 4 were determined using an allele-specific PCR and Sanger sequencing. Of the 83 patient specimens tested, bacterial characterization could be performed on 46 specimens (55%). M. tuberculosis lineage 4 was present in 26 patient specimens (56%), and non-lineage 4 was identified in 10 cases (22%). Moreover, genomic heterogeneity was detected in the CNS specimens of 7 adults and 3 children. We could show that infection of the CNS is not restricted to a single M. tuberculosis lineage and that even young children with rapid progression of disease can harbor more than one M. tuberculosis lineage in the CNS.

Reference set of Mycobacterium tuberculosis clinical strains: A tool for research and product development

The Mycobacterium tuberculosis complex (MTBC) causes tuberculosis (TB) in humans and various other mammals. The human-adapted members of the MTBC comprise seven phylogenetic lineages that differ in their geographical distribution. There is growing evidence that this phylogeographic diversity modulates the outcome of TB infection and disease. For decades, TB research and development has focused on the two canonical MTBC laboratory strains H37Rv and Erdman, both of which belong to Lineage 4. Relying on only a few laboratory-adapted strains can be misleading as study results might not be directly transferrable to clinical settings where patients are infected with a diverse array of strains, including drug-resistant variants. Here, we argue for the need to expand TB research and development by incorporating the phylogenetic diversity of the MTBC. To facilitate such work, we have assembled a group of 20 genetically well-characterized clinical strains representing the seven known human-adapted MTBC lineages. With the "MTBC clinical strains reference set" we aim to provide a standardized resource for the TB community. We hope it will enable more direct comparisons between studies that explore the physiology of MTBC beyond the laboratory strains used thus far. We anticipate that detailed phenotypic analyses of this reference strain set will increase our understanding of TB biology and assist in the development of new control tools that are broadly effective.

A mycobacterial disease is associated with the silent mass mortality of the pen shell Pinna nobilis along the Tyrrhenian coastline of Italy

Disease is an increasing threat for marine bivalves worldwide. Recently, a mass mortality event (MME) impacting the bivalve Pinna nobilis was detected across a wide geographical area of the Spanish Mediterranean Sea and linked to a haplosporidian parasite. In 2017–2018, mass mortality events affecting the pen shell Pinna nobilis were recorded in two different regions of Italy, Campania and Sicily, in the Tyrrhenian Sea (Mediterranean Sea). Histopathological and molecular examinations of specimens showed the presence of Haplosporidium sp. in only one specimen in one area. Conversely, in all of the surveyed moribund animals, strong inflammatory lesions at the level of connective tissue surrounding the digestive system and gonads and linked to the presence of intracellular Zhiel-Neelsen-positive bacteria were observed. Molecular analysis of all of the diseased specimens (13) confirmed the presence of a Mycobacterium. Blast analysis of the sequences from all of the areas revealed that they were grouped together with the human mycobacterium M. sherrisii close to the group including M. shigaense, M. lentiflavum and M. simiae. Based on pathological and molecular findings, it is proposed that a mycobacterial disease is associated with the mortality episodes of Pinna nobilis, indicating that, at this time, Haplosporidium sp. is not responsible for these events in Campanian and Sicilian waters.

A New Phylogenetic Framework for the Animal-Adapted Mycobacterium tuberculosis Complex

Tuberculosis (TB) affects humans and other animals and is caused by bacteria from the Mycobacterium tuberculosis complex (MTBC). Previous studies have shown that there are at least nine members of the MTBC infecting animals other than humans; these have also been referred to as ecotypes. However, the ecology and the evolution of these animal-adapted MTBC ecotypes are poorly understood. Here we screened 12,886 publicly available MTBC genomes and newly sequenced 17 animal-adapted MTBC strains, gathering a total of 529 genomes of animal-adapted MTBC strains. Phylogenomic and comparative analyses confirm that the animal-adapted MTBC members are paraphyletic with some members more closely related to the human-adapted Mycobacterium africanum Lineage 6 than to other animal-adapted strains. Furthermore, we identified four main animal-adapted MTBC clades that might correspond to four main host shifts; two of these clades are hypothesized to reflect independent cattle domestication events. Contrary to what would be expected from an obligate pathogen, MTBC nucleotide diversity was not positively correlated with host phylogenetic distances, suggesting that host tropism in the animal-adapted MTBC seems to be driven by contact rates and demographic aspects of the host population rather by than host relatedness. By combining phylogenomics with ecological data, we propose an evolutionary scenario in which the ancestor of Lineage 6 and all animal-adapted MTBC ecotypes was a generalist pathogen that subsequently adapted to different host species. This study provides a new phylogenetic framework to better understand the evolution of the different ecotypes of the MTBC and guide future work aimed at elucidating the molecular mechanisms underlying host range.

Correlation of genetic diversity between hosts and parasites in Entamoeba nuttalli isolates from Tibetan and rhesus macaques in China

Entamoeba nuttalli infection is prevalent in captive and wild macaques. Recent studies have suggested that genotypes of E. nuttalli isolates are correlated with the geographical distribution of host macaques. Correlation of amoebic genotypes with genetic diversity of host macaques was analyzed in present study. Sixty fresh stool samples were obtained from wild Tibetan macaques living in Mount Huang (HS) of the An-hui Province in China. PCR analysis revealed that the most prevalent Entamoeba species was E. chattoni (E. polecki ST2) (86.7%) followed by E. nuttalli (58.3%) and E. coli (25%). Six E. nuttalli HS isolates were successfully cultured. The tRNA-linked short tandem repeat (STR) loci and serine-rich protein gene of E. nuttalli isolates from four different regions of China (Mount Long-hu, Gui-yang, Mount E-mei, and HS, the former three isolates were obtained in previous studies) were studied and high numbers of polymorphisms were detected. When genetic diversity of different populations of E. nuttalli isolates was compared with geographical distance, an r² value of 0.919 was assigned by a Mantel test based on the tRNA-STR loci. In host macaques, the mtDNA HVS-I gene was also highly polymorphic in each of the genomes. Multiple regression analysis using E. nuttalli tRNA-STR loci genetic, macaque mtDNA HVS-I gene, and geographic distances showed an r² value of 0.943, indicating that a higher relevance was demonstrated when geographic and host gene factors were considered. Analysis of genetic factor of host would benefit for better understanding of the evolution of E. nuttalli.