Tracking SARS-CoV-2 introductions in Mozambique using pandemic-scale phylogenies: a retrospective observational study.

BACKGROUND: From the start of the SARS-CoV-2 outbreak, global sequencing efforts have generated an unprecedented amount of genomic data. Nonetheless, unequal sampling between high-income and low-income countries hinders the implementation of genomic surveillance systems at the global and local level. Filling the knowledge gaps of genomic information and understanding pandemic dynamics in low-income countries is essential for public health decision making and to prepare for future pandemics. In this context, we aimed to discover the timing and origin of SARS-CoV-2 variant introductions in Mozambique, taking advantage of pandemic-scale phylogenies. METHODS: We did a retrospective, observational study in southern Mozambique. Patients from Manhiça presenting with respiratory symptoms were recruited, and those enrolled in clinical trials were excluded. Data were included from three sources: (1) a prospective hospital-based surveillance study (MozCOVID), recruiting patients living in Manhiça, attending the Manhiça district hospital, and fulfilling the criteria of suspected COVID-19 case according to WHO; (2) symptomatic and asymptomatic individuals with SARS-CoV-2 infection recruited by the National Surveillance system; and (3) sequences from SARS-CoV-2-infected Mozambican cases deposited on the Global Initiative on Sharing Avian Influenza Data database. Positive samples amenable for sequencing were analysed. We used Ultrafast Sample placement on Existing tRees to understand the dynamics of beta and delta waves, using available genomic data. This tool can reconstruct a phylogeny with millions of sequences by efficient sample placement in a tree. We reconstructed a phylogeny (~7·6 million sequences) adding new and publicly available beta and delta sequences. FINDINGS: A total of 5793 patients were recruited between Nov 1, 2020, and Aug 31, 2021. During this time, 133 328 COVID-19 cases were reported in Mozambique. 280 good quality new SARS-CoV-2 sequences were obtained after the inclusion criteria were applied and an additional 652 beta (B.1.351) and delta (B.1.617.2) public sequences were included from Mozambique. We evaluated 373 beta and 559 delta sequences. We identified 187 beta introductions (including 295 sequences), divided in 42 transmission groups and 145 unique introductions, mostly from South Africa, between August, 2020 and July, 2021. For delta, we identified 220 introductions (including 494 sequences), with 49 transmission groups and 171 unique introductions, mostly from the UK, India, and South Africa, between April and November, 2021. INTERPRETATION: The timing and origin of introductions suggests that movement restrictions effectively avoided introductions from non-African countries, but not from surrounding countries. Our results raise questions about the imbalance between the consequences of restrictions and health benefits. This new understanding of pandemic dynamics in Mozambique can be used to inform public health interventions to control the spread of new variants. FUNDING: European and Developing Countries Clinical Trials, European Research Council, Bill & Melinda Gates Foundation, and Agència de Gestió d'Ajuts Universitaris i de Recerca.

The mineralocorticoid receptor modulates timing and location of genomic binding by glucocorticoid receptor in response to synthetic glucocorticoids in keratinocytes.

Glucocorticoids (GCs) exert potent antiproliferative and anti-inflammatory properties, explaining their therapeutic efficacy for skin diseases. GCs act by binding to the GC receptor (GR) and the mineralocorticoid receptor (MR), co-expressed in classical and non-classical targets including keratinocytes. Using knockout mice, we previously demonstrated that GR and MR exert essential nonoverlapping functions in skin homeostasis. These closely related receptors may homo- or heterodimerize to regulate transcription, and theoretically bind identical GC-response elements (GRE). We assessed the contribution of MR to GR genomic binding and the transcriptional response to the synthetic GC dexamethasone (Dex) using control (CO) and MR knockout (MREKO ) keratinocytes. GR chromatin immunoprecipitation (ChIP)-seq identified peaks common and unique to both genotypes upon Dex treatment (1 h). GREs, AP-1, TEAD, and p53 motifs were enriched in CO and MREKO peaks. However, GR genomic binding was 35% reduced in MREKO , with significantly decreased GRE enrichment, and reduced nuclear GR. Surface plasmon resonance determined steady state affinity constants, suggesting preferred dimer formation as MR-MR > GR-MR ~ GR-GR; however, kinetic studies demonstrated that GR-containing dimers had the longest lifetimes. Despite GR-binding differences, RNA-seq identified largely similar subsets of differentially expressed genes in both genotypes upon Dex treatment (3 h). However, time-course experiments showed gene-dependent differences in the magnitude of expression, which correlated with earlier and more pronounced GR binding to GRE sites unique to CO including near Nr3c1. Our data show that endogenous MR has an impact on the kinetics and differential genomic binding of GR, affecting the time-course, specificity, and magnitude of GC transcriptional responses in keratinocytes.

Epidemiological and Genomic Analysis of a Large SARS-CoV-2 Outbreak in a Long-Term Care Facility in Catalonia, Spain.

Limiting outbreaks in long-term care facilities (LTCFs) is a cornerstone strategy to avoid an excess of COVID-19-related morbidity and mortality and to reduce its burden on the health system. We studied a large outbreak that occurred at an LTCF, combining methods of classical and genomic epidemiology analysis. The outbreak lasted for 31 days among residents, with an attack rate of 98% and 57% among residents and staff, respectively. The case fatality rate among residents was 16% (n = 15). Phylogenetic analysis of 59 SARS-CoV-2 isolates revealed the presence of two closely related viral variants in all cases (B.1.177 lineage), revealing a far more complex outbreak than initially thought and suggesting an initial spread driven by staff members. In turn, our results suggest that resident relocations to mitigate viral spread might have increased the risk of infection for staff members, creating secondary chains of transmission that were responsible for prolonging the outbreak. Our results highlight the importance of considering unnoticed chains of transmission early during an outbreak and making an adequate use and interpretation of diagnostic tests. Outbreak containment measures should be carefully tailored to each LTCF. IMPORTANCE The impact of COVID-19 on long-term care facilities (LTCFs) has been disproportionately large due to the high frailty of the residents. Here, we report epidemiological and genomic findings of a large outbreak that occurred at an LTCF, which ultimately affected almost all residents and nearly half of staff members. We found that the outbreak was initially driven by staff members; however, later resident relocation to limit the outbreak resulted in transmission from residents to staff members, evidencing the complexity and different phases of the outbreak. The phylogenetic analysis of SARS-CoV-2 isolates indicated that two closely related variants were responsible for the large outbreak. Our study highlights the importance of combining methods of classical and genomic epidemiology to take appropriate outbreak containment measures in LTCFs.

Fine-grain population structure and transmission patterns of Mycobacterium tuberculosis in southern Mozambique, a high TB/HIV burden area.

Genomic studies of the Mycobacterium tuberculosis complex (MTBC) might shed light on the dynamics of its transmission, especially in high-burden settings, where recent outbreaks are embedded in the complex natural history of the disease. To this end, we conducted a 1 year prospective surveillance-based study in Mozambique. We applied whole-genome sequencing (WGS) to 295 positive cultures. We fully characterized MTBC isolates by phylogenetics and dating evaluation, and carried out a molecular epidemiology analysis to investigate further associations with pre-defined transmission risk factors. The majority of strains (49.5%, 136/275) belonged to lineage (L) 4; 57.8 % of them (159/275) were in genomic transmission clusters (cut-off 5 SNPs), and a strikingly high proportion (45.5%) shared an identical genotype (0 SNP pairwise distance). We found two 'likely endemic' clades, comprising 67 strains, belonging to L1.2, which dated back to the late 19th century and were associated with recent spread among people living with human immunodeficiency virus (PLHIV). We describe for the first time the population structure of MTBC in our region, a high tuberculosis (TB)/HIV burden area. Clustering analysis revealed an unforeseen pattern of spread and high rates of progression to active TB, suggesting weaknesses in TB control activities. The long-term presence of local strains in Mozambique, which were responsible for large transmission among HIV/TB-coinfected patients, calls into question the role of HIV in TB transmission.

Gene evolutionary trajectories in Mycobacterium tuberculosis reveal temporal signs of selection.

Genetic differences between different Mycobacterium tuberculosis complex (MTBC) strains determine their ability to transmit within different host populations, their latency times, and their drug resistance profiles. Said differences usually emerge through de novo mutations and are maintained or discarded by the balance of evolutionary forces. Using a dataset of ∼5,000 strains representing global MTBC diversity, we determined the past and present selective forces that have shaped the current variability observed in the pathogen population. We identified regions that have evolved under changing types of selection since the time of the MTBC common ancestor. Our approach highlighted striking differences in the genome regions relevant for host­pathogen interaction and, in particular, suggested an adaptive role for the sensor protein of two-component systems. In addition, we applied our approach to successfully identify potential determinants of resistance to drugs administered as second-line tuberculosis treatments.

Spatial and temporal distribution of SARS-CoV-2 diversity circulating in wastewater.

Wastewater-based epidemiology (WBE) has proven to be an effective tool for epidemiological surveillance of SARS-CoV-2 during the current COVID-19 pandemic. Furthermore, combining WBE together with high-throughput sequencing techniques can be useful for the analysis of SARS-CoV-2 viral diversity present in a given sample. The present study focuses on the genomic analysis of SARS-CoV-2 in 76 sewage samples collected during the three epidemiological waves that occurred in Spain from 14 wastewater treatment plants distributed throughout the country. The results obtained demonstrate that the metagenomic analysis of SARS-CoV-2 in wastewater allows the detection of mutations that define the B.1.1.7 lineage and the ability of the technique to anticipate the detection of certain mutations before they are detected in clinical samples. The study proves the usefulness of sewage sequencing to track Variants of Concern that can complement clinical testing to help in decision-making and in the analysis of the evolution of the pandemic.

Recurrences of multidrug-resistant tuberculosis: Strains involved, within-host diversity, and fine-tuned allocation of reinfections.

Recurrent tuberculosis occurs due to exogenous reinfection or reactivation/persistence. We analysed 90 sequential MDR Mtb isolates obtained in Argentina from 27 patients with previously diagnosed MDR-TB that recurred in 2018 (1-10 years, 2-10 isolates per patient). Three long-term predominant strains were responsible for 63% of all MDR-TB recurrences. Most of the remaining patients were infected by strains different from each other. Reactivation/persistence of the same strain caused all but one recurrence, which was due to a reinfection with a predominant strain. One of the prevalent strains showed marked stability in the recurrences, while in another strain higher SNP-based diversity was observed. Comparisons of intra- versus inter-patient SNP distances identified two possible reinfections with closely related variants circulating in the community. Our results show a complex scenario of MDR-TB infections in settings with predominant MDR Mtb strains.

Overlapping prison/community tuberculosis outbreaks in Costa Rica revealed by alternative analysis of suboptimal material.

Costa Rica has a low incidence of tuberculosis. Thus, identifying transmission hotspots is key to implement interventions. A tuberculosis outbreak was suspected in a prison in Costa Rica. Given the suboptimal quality of the samples received in our laboratory in Madrid, we applied alternative schemes for their analysis. In the first scheme, we bypassed the standard approach of applying systematic mycobacterial interspersed repetitive units-variable number of tandem repeats (MIRU-VNTR) and used a strain-specific polymerase chain reaction (PCR) that allowed identifying a cluster involving six cases (C1). The second scheme followed the canonical MIRU-VNTR path coupled with a whole-genomic amplification step, by which a second unsuspected overlapping cluster (C2), was detected in the same prison. These findings justified the implementation of a surveillance programme adapted to local resources based on a tailored multiplex allele-specific oligonucleotide (ASO)-PCR targeting C1 and C2. Presence of the C2 strain at a different prison was determined. ASO-PCR was applied extensively and alerted to the active circulation of one of the strains within and beyond prisons. Our study shows that alternative methodological strategies may provide useful data in settings with lack of resources for performing systematic standard molecular epidemiology programmes and/or with suboptimal material for analysis.

Genomic analysis of Mycobacterium brumae sustains its nonpathogenic and immunogenic phenotype.

Mycobacterium brumae is a rapid-growing, non-pathogenic Mycobacterium species, originally isolated from environmental and human samples in Barcelona, Spain. Mycobacterium brumae is not pathogenic and it's in vitro phenotype and immunogenic properties have been well characterized. However, the knowledge of its underlying genetic composition is still incomplete. In this study, we first describe the 4 Mb genome of the M. brumae type strain ATCC 51384T assembling PacBio reads, and second, we assess the low intraspecies variability by comparing the type strain with Illumina reads from three additional strains. Mycobacterium brumae genome is composed of a circular chromosome with a high GC content of 69.2% and containing 3,791 CDSs, 97 pseudogenes, one prophage and no CRISPR loci. Mycobacterium brumae has shown no pathogenic potential in in vivo experiments, and our genomic analysis confirms its phylogenetic position with other non-pathogenic and rapid growing mycobacteria. Accordingly, we determined the absence of virulence-related genes, such as ESX-1 locus and most PE/PPE genes, among others. Although the immunogenic potential of M. brumae was proved to be as high as Mycobacterium bovis BCG, the only mycobacteria licensed to treat cancer, the genomic content of M. tuberculosis T cell and B cell antigens in M. brumae genome is considerably lower than those antigens present in M. bovis BCG genome. Overall, this work provides relevant genomic data on one of the species of the mycobacterial genus with high therapeutic potential.

Accuracy of an amplicon-sequencing nanopore approach to identify variants in tuberculosis drug-resistance-associated genes.

A rapid and accurate diagnostic assay represents an important means to detect Mycobacterium tuberculosis, identify drug-resistant strains and ensure treatment success. Currently employed techniques to diagnose drug-resistant tuberculosis include slow phenotypic tests or more rapid molecular assays that evaluate a limited range of drugs. Whole-genome-sequencing-based approaches can detect known drug-resistance-conferring mutations and novel variations; however, the dependence on growing samples in culture, and the associated delays in achieving results, represents a significant limitation. As an alternative, targeted sequencing strategies can be directly performed on clinical samples at high throughput. This study proposes a targeted sequencing assay to rapidly detect drug-resistant strains of M. tuberculosis using the Nanopore MinION sequencing platform. We designed a single-tube assay that targets nine genes associated with drug resistance to seven drugs and two phylogenetic-determining regions to determine strain lineage and tested it in nine clinical isolates and six sputa. The study's main aim is to calibrate MinNION variant calling to detect drug-resistance-associated mutations with different frequencies to match the accuracy of Illumina (the current gold-standard sequencing technology) from both culture and sputum samples. After calibrating Nanopore MinION variant calling, we demonstrated 100% agreement between Illumina WGS and our MinION set up to detect known drug resistance and phylogenetic variants in our dataset. Importantly, other variants in the amplicons are also detected, decreasing the recall. We identify minority variants and insertions/deletions as crucial bioinformatics challenges to fully reproduce Illumina WGS results.

Evolutionary and Phenotypic Characterization of Two Spike Mutations in European Lineage 20E of SARS-CoV-2.

We have detected two mutations in the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) at amino acid positions 1163 and 1167 that appeared independently in multiple transmission clusters and different genetic backgrounds. Furthermore, both mutations appeared together in a cluster of 1,627 sequences belonging to clade 20E. This cluster is characterized by 12 additional single nucleotide polymorphisms but no deletions. The available structural information on the S protein in the pre- and postfusion conformations predicts that both mutations confer rigidity, which could potentially decrease viral fitness. Accordingly, we observed reduced infectivity of this spike genotype relative to the ancestral 20E sequence in vitro, and the levels of viral RNA in nasopharyngeal swabs were not significantly higher. Furthermore, the mutations did not impact thermal stability or antibody neutralization by sera from vaccinated individuals but moderately reduce neutralization by convalescent-phase sera from the early stages of the pandemic. Despite multiple successful appearances of the two spike mutations during the first year of SARS-CoV-2 evolution, the genotype with both mutations was displaced upon the expansion of the 20I (Alpha) variant. The midterm fate of the genotype investigated was consistent with the lack of advantage observed in the clinical and experimental data. IMPORTANCE We observed repeated, independent emergence of mutations in the SARS-CoV-2 spike involving amino acids 1163 and 1167, within the HR2 functional motif. Conclusions derived from evolutionary and genomic diversity analysis suggest that the co-occurrence of both mutations might pose an advantage for the virus and therefore a threat to effective control of the epidemic. However, biological characterization, including in vitro experiments and analysis of clinical data, indicated no clear benefit in terms of stability or infectivity. In agreement with this, continuous epidemiological surveillance conducted months after the first observations revealed that both mutations did not successfully outcompete other variants and stopped circulating 9 months after their initial detection. Additionally, we evaluated the potential of both mutations to escape neutralizing antibodies, finding that the presence of these two mutations on their own is not likely to confer antibody escape. Our results provide an example of how newly emerged spike mutations can be assessed to better understand the risk posed by new variants and indicate that some spike mutations confer no clear advantage to the virus despite independently emerging multiple times and are eventually displaced by fitter variants.

An evolutionary functional genomics approach identifies novel candidate regions involved in isoniazid resistance in Mycobacterium tuberculosis.

Efforts to eradicate tuberculosis are hampered by the rise and spread of antibiotic resistance. Several large-scale projects have aimed to specifically link clinical mutations to resistance phenotypes, but they were limited in both their explanatory and predictive powers. Here, we combine functional genomics and phylogenetic associations using clinical strain genomes to decipher the architecture of isoniazid resistance and search for new resistance determinants. This approach has allowed us to confirm the main target route of the antibiotic, determine the clinical relevance of redox metabolism as an isoniazid resistance mechanism and identify novel candidate genes harboring resistance mutations in strains with previously unexplained isoniazid resistance. This approach can be useful for characterizing how the tuberculosis bacilli acquire resistance to new antibiotics and how to forestall them.

The first wave of the COVID-19 epidemic in Spain was associated with early introductions and fast spread of a dominating genetic variant.

The coronavirus disease 2019 (COVID-19) pandemic has affected the world radically since 2020. Spain was one of the European countries with the highest incidence during the first wave. As a part of a consortium to monitor and study the evolution of the epidemic, we sequenced 2,170 samples, diagnosed mostly before lockdown measures. Here, we identified at least 500 introductions from multiple international sources and documented the early rise of two dominant Spanish epidemic clades (SECs), probably amplified by superspreading events. Both SECs were related closely to the initial Asian variants of SARS-CoV-2 and spread widely across Spain. We inferred a substantial reduction in the effective reproductive number of both SECs due to public-health interventions (Re < 1), also reflected in the replacement of SECs by a new variant over the summer of 2020. In summary, we reveal a notable difference in the initial genetic makeup of SARS-CoV-2 in Spain compared with other European countries and show evidence to support the effectiveness of lockdown measures in controlling virus spread, even for the most successful genetic variants.

Complete Analysis of the Epidemiological Scenario around a SARS-CoV-2 Reinfection: Previous Infection Events and Subsequent Transmission.

The first descriptions of reinfection by SARS-CoV-2 have been recently reported. However, these studies focus exclusively on the reinfected case, without considering the epidemiological context of the event. Our objectives were to perform a complete analysis of the sequential infections and community transmission events around a SARS-CoV-2 reinfection, including the infection events preceding it, the exposure, and subsequent transmissions. Our analysis was supported by host genetics, viral whole-genome sequencing, phylogenomic viral population analysis, and refined epidemiological data obtained from interviews with the involved subjects. The reinfection involved a 53-year-old woman with asthma (Case A), with a first COVID-19 episode in April 2020 and a much more severe second episode 4-1/2 months later, with SARS-CoV-2 seroconversion in August, that required hospital admission. An extended genomic analysis allowed us to demonstrate that the strain involved in Case A's reinfection was circulating in the epidemiological context of Case A and was also transmitted subsequently from Case A to her family context. The reinfection was also supported by a phylogenetic analysis, including 348 strains from Madrid, which revealed that the strain involved in the reinfection was circulating by the time Case A suffered the second episode, August-September 2020, but absent at the time range corresponding to Case A's first episode. IMPORTANCE We present the first complete analysis of the epidemiological scenario around a reinfection by SARS-CoV-2, more severe than the first episode, including three cases preceding the reinfection, the reinfected case per se, and the subsequent transmission to another seven cases.

Integrative transnational analysis to dissect tuberculosis transmission events along the migratory route from Africa to Europe.

BACKGROUND: Growing international migration has increased the complexity of tuberculosis transmission patterns. Italy's decision to close its borders in 2018 made of Spain the new European porte entrée for migration from the Horn of Africa (HA). In one of the first rescues of migrants from this region at the end of 2018, tuberculosis was diagnosed in eight subjects, mainly unaccompanied minors. METHODS: Mycobacterium tuberculosis isolates from these recently arrived migrants were analysed by Mycobacterial Interspersed Repetitive-Unit/Variable-Number of Tandem Repeat (MIRU-VNTR) and subsequent whole genome sequencing (WGS) analysis. Data were compared with those from collections from other European countries receiving migrants from the HA and a strain-specific PCR was applied for a fast searching of common strains. Infections in a cellular model were performed to assess strain virulence. RESULTS: MIRU-VNTR analysis allowed identifying an epidemiological cluster involving three of the eight cases from Somalia (0 single-nucleotide polymorphisms between isolates, HA cluster). Following detailed interviews revealed that two of these cases had shared the same migratory route in most of the trip and had spent a long time at a detention camp in Libya. To confirm potential en route transmission for the three cases, we searched the same strain in collections from other European countries receiving migrants from the HA. MIRU-VNTR, WGS and a strain-specific PCR for the HA strain were applied. The same strain was identified in 12 cases from Eritrea diagnosed soon after their arrival in 2018 to the Netherlands, Belgium and Italy. Intracellular replication rate of the strain did not reveal abnormal virulence. CONCLUSIONS: Our study suggests a potential en route transmission of a pan-susceptible strain, which caused at least 15 tuberculosis cases in Somalian and Eritrean migrants diagnosed in four different European countries.

SARS-CoV-2 Outbreak on a Spanish Mink Farm: Epidemiological, Molecular, and Pathological Studies.

Farmed minks have been reported to be highly susceptible to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and may represent a risk to humans. In this study, we describe the first outbreak of SARS-CoV-2 occurred on a mink farm in Spain, between June and July 2020, involving 92,700 animals. The outbreak started shortly after some farm workers became seropositive for SARS-CoV-2. Minks showed no clinical signs compatible with SARS-CoV-2 infection throughout the outbreak. Samples from 98 minks were collected for histopathological, serological, and molecular studies. Twenty out of 98 (20.4%) minks were positive by RT-qPCR and 82 out 92 (89%) seroconverted. This finding may reflect a rapid spread of the virus at the farm with most of the animals overcoming the infection. Additionally, SARS-CoV-2 was detected by RT-qPCR in 30% of brain samples from positive minks. Sequencing analysis showed that the mink sequences were not closely related with the other mink SARS-CoV-2 sequences available, and that this mink outbreak has its probable origin in one of the genetic variants that were prevalent in Spain during the first COVID-19 epidemic wave. Histological studies revealed bronchointerstitial pneumonia in some animals. Immunostaining of viral nucleocapsid was also observed in nasal turbinate tissue. Farmed minks could therefore constitute an important SARS-CoV-2 reservoir, contributing to virus spread among minks and humans. Consequently, continuous surveillance of mink farms is needed.

Glucocorticoid-dependent transcription in skin requires epidermal expression of the glucocorticoid receptor and is modulated by the mineralocorticoid receptor.

Glucocorticoid (GC) actions are mediated through two closely related ligand-dependent transcription factors, the GC receptor (GR) and the mineralocorticoid receptor (MR). Given the wide and effective use of GCs to combat skin inflammatory diseases, it is important to understand the relative contribution of these receptors to the transcriptional response to topical GCs. We evaluated the gene expression profiles in the skin of mice with epidermal-specific loss of GR (GREKO), MR (MREKO), or both (double KO; DKO) in response to dexamethasone (Dex). The overall transcriptional response was abolished in GREKO and DKO skin suggesting dependence of the underlying dermis on the presence of epidermal GR. Indeed, the observed dermal GC resistance correlated with a constitutive decrease in GR activity and up-regulation of p38 activity in this skin compartment. Upon Dex treatment, more than 90% of differentially expressed genes (DEGs) in CO overlapped with MREKO. However, the number of DEGs was fourfold increased and the magnitude of response was higher in MREKO vs CO, affecting both gene induction and repression. Taken together our data reveal that, in the cutaneous transcriptional response to GCs mediated through endogenous receptors, epidermal GR is mandatory while epidermal MR acts as a chief modulator of gene expression.

Mycobacterium tuberculosis associated with severe tuberculosis evades cytosolic surveillance systems and modulates IL-1ß production.

Genetic diversity of Mycobacterium tuberculosis affects immune responses and clinical outcomes of tuberculosis (TB). However, how bacterial diversity orchestrates immune responses to direct distinct TB severities is unknown. Here we study 681 patients with pulmonary TB and show that M. tuberculosis isolates from cases with mild disease consistently induce robust cytokine responses in macrophages across multiple donors. By contrast, bacteria from patients with severe TB do not do so. Secretion of IL-1ß is a good surrogate of the differences observed, and thus to classify strains as probable drivers of different TB severities. Furthermore, we demonstrate that M. tuberculosis isolates that induce low levels of IL-1ß production can evade macrophage cytosolic surveillance systems, including cGAS and the inflammasome. Isolates exhibiting this evasion strategy carry candidate mutations, generating sigA recognition boxes or affecting components of the ESX-1 secretion system. Therefore, we provide evidence that M. tuberculosis strains manipulate host-pathogen interactions to drive variable TB severities.

Author Correction: Integrating standardized whole genome sequence analysis with a global Mycobacterium tuberculosis antibiotic resistance knowledgebase.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.