ESX-1-Independent Horizontal Gene Transfer by Mycobacterium tuberculosis Complex Strains

Pangenome-wide association study reveals the selective absence of CRISPR genes (Rv2816c-19c) in drug-resistant Mycobacterium tuberculosis

The presence of intermittently dispersed insertion sequences and transposases in the Mycobacterium tuberculosis (Mtb) genome makes intra-genome recombination events inevitable. Understanding their effect on the gene repertoires (GR), which may contribute to the development of drug-resistant Mtb, is critical. In this study, publicly available WGS data of clinical Mtb isolates (endemic region n = 2,601; non-endemic region n = 1,130) were de novo assembled, filtered, scaffolded into assemblies, and functionally annotated. Out of 2,601 Mtb WGS data sets from endemic regions, 2,184 (drug resistant/sensitive: 1,386/798) qualified as high quality. We identified 3,784 core genes, 123 softcore genes, 224 shell genes, and 762 cloud genes in the pangenome of Mtb clinical isolates from endemic regions. Sets of 33 and 39 genes showed positive and negative associations (P < 0.01) with drug resistance status, respectively. Gene ontology clustering showed compromised immunity to phages and impaired DNA repair in drug-resistant Mtb clinical isolates compared to the sensitive ones. Multidrug efflux pump repressor genes (Rv3830c and Rv3855c) and CRISPR genes (Rv2816c-19c) were absent in the drug-resistant Mtb. A separate WGS data analysis of drug-resistant Mtb clinical isolates from the Netherlands (n = 1130) also showed the absence of CRISPR genes (Rv2816c-17c). This study highlights the role of CRISPR genes in drug resistance development in Mtb clinical isolates and helps in understanding its evolutionary trajectory and as useful targets for diagnostics development.IMPORTANCEThe results from the present Pan-GWAS study comparing gene sets in drug-resistant and drug-sensitive Mtb clinical isolates revealed intricate presence-absence patterns of genes encoding DNA-binding proteins having gene regulatory as well as DNA modification and DNA repair roles. Apart from the genes with known functions, some uncharacterized and hypothetical genes that seem to have a potential role in drug resistance development in Mtb were identified. We have been able to extrapolate many findings of the present study with the existing literature on the molecular aspects of drug-resistant Mtb, further strengthening the relevance of the results presented in this study.

Evolution and emergence of Mycobacterium tuberculosis

Tuberculosis (TB) remains one of the deadliest infectious diseases in human history, prevailing even in the 21st century. The causative agents of TB are represented by a group of closely related bacteria belonging to the Mycobacterium tuberculosis complex (MTBC), which can be subdivided into several lineages of human- and animal-adapted strains, thought to have shared a last common ancestor emerged by clonal expansion from a pool of recombinogenic Mycobacterium canettii-like tubercle bacilli. A better understanding of how MTBC populations evolved from less virulent mycobacteria may allow for discovering improved TB control strategies and future epidemiologic trends. In this review, we highlight new insights into the evolution of mycobacteria at the genus level, describing different milestones in the evolution of mycobacteria, with a focus on the genomic events that have likely enabled the emergence and the dominance of the MTBC. We also review the recent literature describing the various MTBC lineages and highlight their particularities and differences with a focus on host preferences and geographic distribution. Finally, we discuss on putative mechanisms driving the evolution of tubercle bacilli and mycobacteria in general, by taking the mycobacteria-specific distributive conjugal transfer as an example.

A smooth tubercle bacillus from Ethiopia phylogenetically close to the Mycobacterium tuberculosis complex

The Mycobacterium tuberculosis complex (MTBC) includes several human- and animal-adapted pathogens. It is thought to have originated in East Africa from a recombinogenic Mycobacterium canettii-like ancestral pool. Here, we describe the discovery of a clinical tuberculosis strain isolated in Ethiopia that shares archetypal phenotypic and genomic features of M. canettii strains, but represents a phylogenetic branch much closer to the MTBC clade than to the M. canettii strains. Analysis of genomic traces of horizontal gene transfer in this isolate and previously identified M. canettii strains indicates a persistent albeit decreased recombinogenic lifestyle near the emergence of the MTBC. Our findings support that the MTBC emergence from its putative free-living M. canettii-like progenitor is evolutionarily very recent, and suggest the existence of a continuum of further extant derivatives from ancestral stages, close to the root of the MTBC, along the Great Rift Valley.

Horizontal Gene Transfer and Drug Resistance Involving Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) acquires drug resistance at a rate comparable to that of bacterial pathogens that replicate much faster and have a higher mutation rate. One explanation for this rapid acquisition of drug resistance in Mtb is that drug resistance may evolve in other fast-replicating mycobacteria and then be transferred to Mtb through horizontal gene transfer (HGT). This paper aims to address three questions. First, does HGT occur between Mtb and other mycobacterial species? Second, what genes after HGT tend to survive in the recipient genome? Third, does HGT contribute to antibiotic resistance in Mtb? I present a conceptual framework for detecting HGT and analyze 39 ribosomal protein genes, 23S and 16S ribosomal RNA genes, as well as several genes targeted by antibiotics against Mtb, from 43 genomes representing all major groups within Mycobacterium. I also included mgtC and the insertion sequence IS6110 that were previously reported to be involved in HGT. The insertion sequence IS6110 shows clearly that the Mtb complex participates in HGT. However, the horizontal transferability of genes depends on gene function, as was previously hypothesized. HGT is not observed in functionally important genes such as ribosomal protein genes, rRNA genes, and other genes chosen as drug targets. This pattern can be explained by differential selection against functionally important and unimportant genes after HGT. Functionally unimportant genes such as IS6110 are not strongly selected against, so HGT events involving such genes are visible. For functionally important genes, a horizontally transferred diverged homologue from a different species may not work as well as the native counterpart, so the HGT event involving such genes is strongly selected against and eliminated, rendering them invisible to us. In short, while HGT involving the Mtb complex occurs, antibiotic resistance in the Mtb complex arose from mutations in those drug-targeted genes within the Mtb complex and was not gained through HGT.

Cold Cas: reevaluating the occurrence of CRISPR/Cas systems in Mycobacteriaceae

Bacterial CRISPR/Cas systems target foreign genetic elements such as phages and regulate gene expression by some pathogens, even in the host. The system is a marker for evolutionary history and has been used for inferences in Mycobacterium tuberculosis for 30 years. However, knowledge about mycobacterial CRISPR/Cas systems remains limited. It is believed that Type III-A Cas systems are exclusive to Mycobacterium canettii and the M. tuberculosis complex (MTBC) of organisms and that very few of the >200 diverse species of non-tuberculous mycobacteria (NTM) possess any CRISPR/Cas system. This study sought unreported CRISPR/Cas loci across NTM to better understand mycobacterial evolution, particularly in species phylogenetically near the MTBC. An analysis of available mycobacterial genomes revealed that Cas systems are widespread across Mycobacteriaceae and that some species contain multiple types. The phylogeny of Cas loci shows scattered presence in many NTM, with variation even within species, suggesting gains/losses of these loci occur frequently. Cas Type III-A systems were identified in pathogenic Mycobacterium heckeshornense and the geological environmental isolate Mycobacterium SM1. In summary, mycobacterial CRISPR/Cas systems are numerous, Type III-A systems are unreliable as markers for MTBC evolution, and mycobacterial horizontal gene transfer appears to be a frequent source of genetic variation.

Comparative genome analysis reveals high-level drug resistance markers in a clinical isolate of Mycobacterium fortuitum subsp. fortuitum MF GZ001

Introduction

Infections caused by non-tuberculosis mycobacteria are significantly worsening across the globe. M. fortuitum complex is a rapidly growing pathogenic species that is of clinical relevance to both humans and animals. This pathogen has the potential to create adverse effects on human healthcare.

Methods

The MF GZ001 clinical strain was collected from the sputum of a 45-year-old male patient with a pulmonary infection. The morphological studies, comparative genomic analysis, and drug resistance profiles along with variants detection were performed in this study. In addition, comparative analysis of virulence genes led us to understand the pathogenicity of this organism.

Results

Bacterial growth kinetics and morphology confirmed that MF GZ001 is a rapidly growing species with a rough morphotype. The MF GZ001 contains 6413573 bp genome size with 66.18 % high G+C content. MF GZ001 possesses a larger genome than other related mycobacteria and included 6156 protein-coding genes. Molecular phylogenetic tree, collinearity, and comparative genomic analysis suggested that MF GZ001 is a novel member of the M. fortuitum complex. We carried out the drug resistance profile analysis and found single nucleotide polymorphism (SNP) mutations in key drug resistance genes such as rpoB, katG, AAC(2')-Ib, gyrA, gyrB, embB, pncA, blaF, thyA, embC, embR, and iniA. In addition, the MF GZ001strain contains mutations in iniA, iniC, pncA, and ribD which conferred resistance to isoniazid, ethambutol, pyrazinamide, and para-aminosalicylic acid respectively, which are not frequently observed in rapidly growing mycobacteria. A wide variety of predicted putative potential virulence genes were found in MF GZ001, most of which are shared with well-recognized mycobacterial species with high pathogenic profiles such as M. tuberculosis and M. abscessus.

Discussion

Our identified novel features of a pathogenic member of the M. fortuitum complex will provide the foundation for further investigation of mycobacterial pathogenicity and effective treatment.

Evolution of Tuberculosis Pathogenesis

Mycobacterium tuberculosis is a globally distributed, lethal pathogen of humans. The virulence armamentarium of M. tuberculosis appears to have been developed on a scaffold of antiphagocytic defenses found among diverse, mostly free-living species of Mycobacterium. Pathoadaptation was further aided by the modularity, flexibility, and interactivity characterizing mycobacterial effectors and their regulators. During emergence of M. tuberculosis, novel genetic material was acquired, created, and integrated with existing tools. The major mutational mechanisms underlying these adaptations are discussed in this review, with examples. During its evolution, M. tuberculosis lost the ability and/or opportunity to engage in lateral gene transfer, but despite this it has retained the adaptability that characterizes mycobacteria. M. tuberculosis exemplifies the evolutionary genomic mechanisms underlying adoption of the pathogenic niche, and studies of its evolution have uncovered a rich array of discoveries about how new pathogens are made. Expected final online publication date for the Annual Review of Microbiology, Volume 76 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

The MtZ Strain: Molecular Characteristics and Outbreak Investigation of the Most Successful Mycobacterium tuberculosis Strain in Aragon Using Whole-Genome Sequencing

Since 2004, a tuberculosis surveillance protocol has been carried out in Aragon, thereby managing to detect all tuberculosis outbreaks that take place in the community. The largest outbreak was caused by a strain named Mycobacterium tuberculosis Zaragoza (MtZ), causing 242 cases as of 2020. The main objective of this work was to analyze this outbreak and the molecular characteristics of this successful strain that could be related to its greater transmission. To do this, we first applied whole-genome sequencing to 57 of the isolates. This revealed two principal transmission clusters and six subclusters arising from them. The MtZ strain belongs to L4.8 and had eight specific single nucleotide polymorphisms (SNPs) in genes considered to be virulence factors [ptpA, mc3D, mc3F, VapB41, pks15 (two SNPs), virS, and VapC50]. Second, a transcriptomic study was carried out to better understand the multiple IS6110 copies present in its genome. This allowed us to observe three effects of IS6110: the disruption of the gene in which the IS6110 is inserted (desA3), the overexpression of a gene (ppe38), and the absence of transcription of genes (cut1:Rv1765c) due to the recombination of two IS6110 copies. Finally, because of the disruption of ppe38 and ppe71 genes by an IS6110, a study of PE_PGRS secretion was carried out, showing that MtZ secretes these factors in higher amounts than the reference strain, thereby differing from the hypervirulent phenotype described for the Beijing strains. In conclusion, MtZ consists of several SNPs in genes related to virulence, pathogenesis, and survival, as well as other genomic polymorphisms, which may be implicated in its success among our population.

The impact of frequently neglected model violations on bacterial recombination rate estimation: a case study in Mycobacterium canettii and Mycobacterium tuberculosis

Mycobacterium canettii is a causative agent of tuberculosis in humans, along with the members of the Mycobacterium tuberculosis complex. Frequently used as an outgroup to the M. tuberculosis complex in phylogenetic analyses, M. canettii is thought to offer the best proxy for the progenitor species that gave rise to the complex. Here, we leverage whole-genome sequencing data and biologically relevant population genomic models to compare the evolutionary dynamics driving variation in the recombining M. canettii with that in the nonrecombining M. tuberculosis complex, and discuss differences in observed genomic diversity in the light of expected levels of Hill-Robertson interference. In doing so, we highlight the methodological challenges of estimating recombination rates through traditional population genetic approaches using sequences called from populations of microorganisms and evaluate the likely mis-inference that arises owing to a neglect of common model violations including purifying selection, background selection, progeny skew, and population size change. In addition, we compare performance when full within-host polymorphism data are utilized, versus the more common approach of basing analyses on within-host consensus sequences.

A Polymorphic Gene within the Mycobacterium smegmatis esx1 Locus Determines Mycobacterial Self-Identity and Conjugal Compatibility

Mycobacteria mediate horizontal gene transfer (HGT) by a process called distributive conjugal transfer (DCT) that is mechanistically distinct from oriT-mediated plasmid transfer. The transfer of multiple, independent donor chromosome segments generates transconjugants with genomes that are mosaic blends of their parents. Previously, we had characterized contact-dependent conjugation between two independent isolates of Mycobacterium smegmatis. Here, we expand our analyses to include five independent isolates of M. smegmatis and establish that DCT is both active and prevalent among natural isolates of M. smegmatis. Two of these five strains were recipients but exhibited distinct conjugal compatibilities with donor strains, suggesting an ability to distinguish between potential donor partners. We determined that a single gene, Msmeg0070, was responsible for conferring mating compatibility using a combination of comparative DNA sequence analysis, bacterial genome-wide association studies (GWAS), and targeted mutagenesis. Msmeg0070 maps within the esx1 secretion locus, and we establish that it confers mycobacterial self-identity with parallels to kin recognition. Similar to other kin model systems, orthologs of Msmeg0070 are highly polymorphic. The identification of a kin recognition system in M. smegmatis reinforces the concept that communication between cells is an important checkpoint prior to DCT commitment and implies that there are likely to be other, unanticipated forms of social behaviors in mycobacteria. IMPORTANCE Conjugation, unlike other forms of HGT, requires direct interaction between two viable bacteria, which must be capable of distinguishing between mating types to allow successful DNA transfer from donor to recipient. We show that the conjugal compatibility of Mycobacterium smegmatis isolates is determined by a single, polymorphic gene located within the conserved esx1 secretion locus. This gene confers self-identity; the expression of identical Msmeg0070 proteins in both donor-recipient partners prevents DNA transfer. The presence of this polymorphic locus in many environmental mycobacteria suggests that kin identification is important in promoting beneficial gene flow between nonkin mycobacteria. Cell-cell communication, mediated by kin recognition and ESX secretion, is a key checkpoint in mycobacterial conjugation and likely plays a more global role in mycobacterial biology.

Crosstalk between the ancestral type VII secretion system ESX-4 and other T7SS in Mycobacterium marinum

The type VII secretion system (T7SS) of Mycobacterium tuberculosis secretes three substrate classes: Esx, Esp, and PE/PPE proteins, that play important roles in bacterial physiology and host interaction. Five subtypes of T7SS, namely ESX-1 to ESX-5, are present in M. tb. ESX-4 is the progenitor of T7SS but its function is not understood. We investigated the ESX-4 system in Mycobacterium marinum. We show that ESX-4 of M. marinum does not secrete its cognate substrates, EsxT and EsxU, under the conditions tested. Paradoxically, the deletion of eccC4, an essential component of ESX-4, resulted in elevated secretion of protein substrates of ESX-1 and ESX-5. Consequently, the ΔeccC4 mutant was more efficient in inducing actin cytoskeleton rearrangement, which led to enhanced phagocytosis by macrophages. Our results reveal an intimate crosstalk between the progenitor of T7SS and its more recent duplication and expansion, and provide new insight into the evolution of T7SS in mycobacteria.

PPE38-Secretion-Dependent Proteins of M. tuberculosis Alter NF-kB Signalling and Inflammatory Responses in Macrophages

It was previously shown that secretion of PE-PGRS and PPE-MPTR proteins is abolished in clinical M. tuberculosis isolates with a deletion in the ppe38-71 operon, which is associated with increased virulence. Here we investigate the proteins dependent on PPE38 for their secretion and their role in the innate immune response using temporal proteomics and protein turnover analysis in a macrophage infection model. A decreased pro-inflammatory response was observed in macrophages infected with PPE38-deficient M. tuberculosis CDC1551 as compared to wild type bacteria. We could show that dampening of the pro-inflammatory response is associated with activation of a RelB/p50 pathway, while the canonical inflammatory pathway is active during infection with wild type M. tuberculosis CDC1551. These results indicate a molecular mechanism by which M. tuberculosis PE/PPE proteins controlled by PPE38 have an effect on modulating macrophage responses through NF-kB signalling.