Spread of a SARS-CoV-2 variant through Europe in the summer of 2020.

Following its emergence in late 2019, the spread of SARS-CoV-21,2 has been tracked by phylogenetic analysis of viral genome sequences in unprecedented detail3-5. Although the virus spread globally in early 2020 before borders closed, intercontinental travel has since been greatly reduced. However, travel within Europe resumed in the summer of 2020. Here we report on a SARS-CoV-2 variant, 20E (EU1), that was identified in Spain in early summer 2020 and subsequently spread across Europe. We find no evidence that this variant has increased transmissibility, but instead demonstrate how rising incidence in Spain, resumption of travel, and lack of effective screening and containment may explain the variant's success. Despite travel restrictions, we estimate that 20E (EU1) was introduced hundreds of times to European countries by summertime travellers, which is likely to have undermined local efforts to minimize infection with SARS-CoV-2. Our results illustrate how a variant can rapidly become dominant even in the absence of a substantial transmission advantage in favourable epidemiological settings. Genomic surveillance is critical for understanding how travel can affect transmission of SARS-CoV-2, and thus for informing future containment strategies as travel resumes.

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.

Monitoring of First-line Drug Resistance Mutations Outside the Scope of Xpert MTB/RIF Ultra is Needed for Successful Control of DR-TB in Southern Mozambique.

Multidrug-resistant(MDR) tuberculosis in Southern Africa is of great concern, exacerbated by the spread of a clone harboring a mutation missed by Xpert Ultra. In Southern Mozambique, the presence of such mutation and rising cases of non-MDR isoniazid resistance highlights the need to ensure accurate detection of antimicrobial-resistance in the country.

Correction: The structural role of SARS-CoV-2 genetic background in the emergence and success of spike mutations: The case of the spike A222V mutation.

[This corrects the article DOI: 10.1371/journal.ppat.1010631.].

The structural role of SARS-CoV-2 genetic background in the emergence and success of spike mutations: The case of the spike A222V mutation.

The S:A222V point mutation, within the G clade, was characteristic of the 20E (EU1) SARS-CoV-2 variant identified in Spain in early summer 2020. This mutation has since reappeared in the Delta subvariant AY.4.2, raising questions about its specific effect on viral infection. We report combined serological, functional, structural and computational studies characterizing the impact of this mutation. Our results reveal that S:A222V promotes an increased RBD opening and slightly increases ACE2 binding as compared to the parent S:D614G clade. Finally, S:A222V does not reduce sera neutralization capacity, suggesting it does not affect vaccine effectiveness.

Dispersal history of SARS-CoV-2 in Galicia, Spain.

The dynamics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission are influenced by a variety of factors, including social restrictions and the emergence of distinct variants. In this study, we delve into the origins and dissemination of the Alpha, Delta, and Omicron-BA.1 variants of concern in Galicia, northwest Spain. For this, we leveraged genomic data collected by the EPICOVIGAL Consortium and from the GISAID database, along with mobility information from other Spanish regions and foreign countries. Our analysis indicates that initial introductions during the Alpha phase were predominantly from other Spanish regions and France. However, as the pandemic progressed, introductions from Portugal and the United States became increasingly significant. The number of detected introductions varied from 96 and 101 for Alpha and Delta to 39 for Omicron-BA.1. Most of these introductions left a low number of descendants (<10), suggesting a limited impact on the evolution of the pandemic in Galicia. Notably, Galicia's major coastal cities emerged as critical hubs for viral transmission, highlighting their role in sustaining and spreading the virus. This research emphasizes the critical role of regional connectivity in the spread of SARS-CoV-2 and offers essential insights for enhancing public health strategies and surveillance measures.

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.

Cumulative incidence of SARS-CoV-2 infection in the general population of the Valencian Community (Spain) after the surge of the Omicron BA.1 variant.

Studies investigating the cumulative incidence of and immune status against SARS-CoV-2 infection provide valuable information for shaping public health decision-making. A cross-sectional study on 935 participants, conducted in the Valencian Community (VC), measuring anti-SARS-CoV-2-receptor binding domain-RBD-total antibodies and anti-Nucleocapsid (N)-IgGs via electrochemiluminescence assays. Quantitation of neutralizing antibodies (NtAb) against ancestral and Omicron BA.1 and BA.2 variants and enumeration of SARS-CoV-2-S specific-IFNγ-producing CD4+ and CD8+ T cells was performed in 100 and 137 participants, respectively. The weighted cumulative incidence was 51.9% (95% confidence interval [CI]: 48.7-55.1) and was inversely related to age. Anti-RBD total antibodies were detected in 97% of participants; vaccinated and SARS-CoV-2-experienced (VAC-ex; n = 442) presented higher levels (p < 0.001) than vaccinated/naïve (VAC-n; n = 472) and nonvaccinated/experienced (UNVAC-ex; n = 63) subjects. Antibody levels correlated inversely with time elapsed since last vaccine dose in VAC-n (Rho, -0.52; p < 0.001) but not in VAC-ex (rho -0.02; p = 0.57). Heterologous booster shots resulted in increased anti-RBD antibody levels compared with homologous schedules in VAC-n, but not in VAC-ex. NtAbs against Omicron BA.1 were detected in 94%, 75%, and 50% of VAC-ex, VAC-n and UNVAC-ex groups, respectively. For Omicron BA.2, the figures were 97%, 84%, and 40%, respectively. SARS-CoV-2-S-reactive IFN-γ T cells were detected in 73%, 75%, and 64% of VAC-ex, VAC-n and UNVAC-ex, respectively. Median frequencies for both T-cell subsets were comparable across groups. In summary, by April 2022, around half of the VC population had been infected with SARS-CoV-2 and, due to extensive vaccination, displayed hybrid immunity.

Severe Acute Respiratory Syndrome Coronavirus 2 Adaptive Immunity in Nursing Home Residents Following a Third Dose of the Comirnaty Coronavirus Disease 2019 Vaccine.

A third Comirnaty vaccine dose increased severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) receptor-binding domain antibody levels (median, 93-fold) and neutralizing antibody titers against Wuhan-Hu-1 (median, 57-fold), Beta (me 22-fold), Delta, (median, 43-fold), and Omicron (median, 8-fold) variants, but had less impact on S-reactive T-cell immunity in nursing home residents.

Whole genome sequencing analysis to evaluate the influence of T2DM on polymorphisms associated with drug resistance in M. tuberculosis.

BACKGROUND: Type 2 diabetes mellitus (T2DM) has been associated with treatment failure, and the development of drug resistance in tuberculosis (TB). Also, whole-genome sequencing has provided a better understanding and allowed the growth of knowledge about polymorphisms in genes associated with drug resistance. Considering the above, this study analyzes genome sequences to evaluate the influence of type 2 diabetes mellitus in the development of mutations related to tuberculosis drug resistance. M. tuberculosis isolates from individuals with (n = 74), and without (n = 74) type 2 diabetes mellitus was recovered from online repositories, and further analyzed. RESULTS: The results showed the presence of 431 SNPs with similar proportions between diabetics, and non-diabetics individuals (48% vs. 52%), but with no significant relationship. A greater number of mutations associated with rifampicin resistance was observed in the T2DM-TB individuals (23.2% vs. 16%), and the exclusive presence of rpoBQ432L, rpoBQ432P, rpoBS441L, and rpoBH445L variants. While these variants are not private to T2DM-TB cases they are globally rare highlighting a potential role of T2DM. The phylogenetic analysis showed 12 sublineages, being 4.1.1.3, and 4.1.2.1 the most prevalent in T2DM-TB individuals but not differing from those most prevalent in their geographic location. Four clonal complexes were found, however, no significant relationship with T2DM was observed. Samples size and potential sampling biases prevented us to look for significant associations. CONCLUSIONS: The occurrence of globally rare rifampicin variants identified only in isolates from individuals with T2DM could be due to the hyperglycemic environment within the host. Therefore, further studies about the dynamics of SNPs' generation associated with antibiotic resistance in patients with diabetes mellitus are necessary.

Systematic Genomic and Clinical Analysis of Severe Acute Respiratory Syndrome Coronavirus 2 Reinfections and Recurrences Involving the Same Strain.

Estimates of the burden of severe acute respiratory syndrome coronavirus 2 reinfections are limited by the scarcity of population-level studies incorporating genomic support. We conducted a systematic study of reinfections in Madrid, Spain, supported by genomic viral analysis and host genetic analysis, to cleanse laboratory errors and to discriminate between reinfections and recurrences involving the same strain. Among the 41,195 cases diagnosed (March 2020-March 2021), 93 (0.23%) had 2 positive reverse transcription PCR tests (55-346 days apart). After eliminating cases with specimens not stored, of suboptimal sequence quality, or belonging to different persons, we obtained valid data from 22 cases. Of those, 4 (0.01%) cases were recurrences involving the same strain; case-patients were 39-93 years of age, and 3 were immunosuppressed. Eighteen (0.04%) cases were reinfections; patients were 19-84 years of age, and most had no relevant clinical history. The second episode was more severe in 8 cases.

SARS-CoV-2 Delta-variant breakthrough infections in nursing home residents at midterm after Comirnaty® COVID-19 vaccination.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Delta variant breakthrough infections in nursing home residents following vaccination with Comirnaty® COVID-19 vaccine were characterized. In total, 201 participants (median age, 87 years; range, 64-100; 133 female) from two nursing homes in the Valencian community (Spain) were included. SARS-CoV-2-Spike (S) antibody responses were determined by a lateral flow immunocromatography (LFIC) assay and by quantitative electrochemiluminescent assay in LFIC-negative participants. SARS-CoV-2-S-IFNγ T cells were enumerated by flow cytometry in 10 participants. Nasopharyngeal SARS-CoV-2 RNA loads were quantified by real-time polymerase chain reaction assays. Vaccine breakthrough COVID-19 due to the Delta variant occurred in 39 residents (median age, 87 years; range, 69-96; 31 female) at a median of 6.5 months after vaccination (nine requiring hospitalization). Breakthrough infections occurred at a higher rate (p < 0.0001) in residents who had not been previously infected with SARS-CoV-2 (naïve) (33/108; 18%) than in those with prior diagnosis of SARS-CoV-2 infection (experienced) (6/93; 6.4%), and were more likely (p < 0.0001) to develop in residents who tested negative by LFIC (20/49) at 3 months after vaccination as compared to their LFIC-positive counterparts (19/142). Among LFIC-negative residents, a trend towards lower plasma anti-RBD antibody levels was noticed in those developing breakthrough infection (p = 0.16). SARS-CoV-2 RNA loads in nasopharyngeal specimens were lower in SARS-CoV-2-experienced residents (p < 0.001) and in those testing positive by LFIC (p = 0.13). The frequency of SARS-CoV-2-S-reactive T cells at 3 months was similar in LFIC-negative residents with (n = 7) or without (n = 3) breakthrough infection. Prior history of SARS-CoV-2 infection and detection of S-reactive antibodies by LFIC at 3 months is associated with a lower risk of Delta-variant breakthrough infection in nursing home residents at midterm after Comirnaty® COVID-19 vaccination.

SARS-CoV-2 Omicron BA.1 variant breakthrough infections in nursing home residents after an homologous third dose of the Comirnaty® COVID-19 vaccine: Looking for correlates of protection.

We investigated whether peripheral blood levels of SARS-CoV-2 Spike (S) receptor binding domain antibodies (anti-RBD), neutralizing antibodies (NtAb) targeting Omicron S, and S-reactive-interferon (IFN)-γ-producing CD4+ and CD8+ T cells measured after a homologous booster dose (3D) with the Comirnaty® vaccine was associated with the likelihood of subsequent breakthrough infections due to the Omicron variant. An observational study including 146 nursing home residents (median age, 80 years; range, 66-99; 109 female) evaluated for an immunological response after 3D (at a median of 16 days). Anti-RBD total antibodies were measured by chemiluminescent immunoassay. NtAb were quantified by an Omicron S pseudotyped virus neutralization assay. SARS-CoV-2-S specific-IFNγ-producing CD4+ and CD8+ T cells were enumerated by whole-blood flow cytometry for intracellular cytokine staining. In total, 33/146 participants contracted breakthrough Omicron infection (symptomatic in 30/33) within 4 months after 3D. Anti-RBD antibody levels were comparable in infected and uninfected participants (21 123 vs. 24 723 BAU/ml; p = 0.34). Likewise, NtAb titers (reciprocal IC50 titer, 157 vs. 95; p = 0.32) and frequency of virus-reactive CD4+ (p = 0.82) and CD8+ (p = 0.91) T cells were similar across participants in both groups. anti-RBD antibody levels and NtAb titers estimated at around the time of infection were also comparable (3445 vs. 4345 BAU/ml; p = 0.59 and 188.5 vs. 88.9; p = 0.70, respectively). Having detectable NtAb against Omicron or SARS-CoV-2-S-reactive-IFNγ-producing CD4+ or CD8+ T cells after 3D was not correlated with increased protection from breakthrough infection (OR, 1.50; p = 0.54; OR, 0.0; p = 0.99 and OR 3.70; p = 0.23, respectively). None of the immune parameters evaluated herein, including NtAb titers against the Omicron variant, may reliably predict at the individual level the risk of contracting COVID-19 due to the Omicron variant in nursing home residents.

Proper Assignation of Reactivation in a COVID-19 Recurrence Initially Interpreted as a Reinfection.

A 77-year-old man (case R) with previous diagnosis of a mild COVID-19 episode was hospitalized 35 days later. On day 23 postadmission, he developed a second COVID-19 episode, now severe, and finally died. Initially, case R's COVID-19 recurrence was interpreted as a reinfection due to the exposure to a SARS-CoV-2 RT-PCR-positive roommate. However, whole-genome sequencing indicated that case R's recurrence corresponded to a reactivation of the strain involved in his first episode. Case R's reactivation had major consequences, leading to a more severe episode, and causing subsequent transmission to another 2 hospitalized patients, 1 of them with fatal outcome.

Prediction of Susceptibility to First-Line Tuberculosis Drugs by DNA Sequencing.

BACKGROUND: The World Health Organization recommends drug-susceptibility testing of Mycobacterium tuberculosis complex for all patients with tuberculosis to guide treatment decisions and improve outcomes. Whether DNA sequencing can be used to accurately predict profiles of susceptibility to first-line antituberculosis drugs has not been clear. METHODS: We obtained whole-genome sequences and associated phenotypes of resistance or susceptibility to the first-line antituberculosis drugs isoniazid, rifampin, ethambutol, and pyrazinamide for isolates from 16 countries across six continents. For each isolate, mutations associated with drug resistance and drug susceptibility were identified across nine genes, and individual phenotypes were predicted unless mutations of unknown association were also present. To identify how whole-genome sequencing might direct first-line drug therapy, complete susceptibility profiles were predicted. These profiles were predicted to be susceptible to all four drugs (i.e., pansusceptible) if they were predicted to be susceptible to isoniazid and to the other drugs or if they contained mutations of unknown association in genes that affect susceptibility to the other drugs. We simulated the way in which the negative predictive value changed with the prevalence of drug resistance. RESULTS: A total of 10,209 isolates were analyzed. The largest proportion of phenotypes was predicted for rifampin (9660 [95.4%] of 10,130) and the smallest was predicted for ethambutol (8794 [89.8%] of 9794). Resistance to isoniazid, rifampin, ethambutol, and pyrazinamide was correctly predicted with 97.1%, 97.5%, 94.6%, and 91.3% sensitivity, respectively, and susceptibility to these drugs was correctly predicted with 99.0%, 98.8%, 93.6%, and 96.8% specificity. Of the 7516 isolates with complete phenotypic drug-susceptibility profiles, 5865 (78.0%) had complete genotypic predictions, among which 5250 profiles (89.5%) were correctly predicted. Among the 4037 phenotypic profiles that were predicted to be pansusceptible, 3952 (97.9%) were correctly predicted. CONCLUSIONS: Genotypic predictions of the susceptibility of M. tuberculosis to first-line drugs were found to be correlated with phenotypic susceptibility to these drugs. (Funded by the Bill and Melinda Gates Foundation and others.).

Contaminant DNA in bacterial sequencing experiments is a major source of false genetic variability.

BACKGROUND: Contaminant DNA is a well-known confounding factor in molecular biology and in genomic repositories. Strikingly, analysis workflows for whole-genome sequencing (WGS) data commonly do not account for errors potentially introduced by contamination, which could lead to the wrong assessment of allele frequency both in basic and clinical research. RESULTS: We used a taxonomic filter to remove contaminant reads from more than 4000 bacterial samples from 20 different studies and performed a comprehensive evaluation of the extent and impact of contaminant DNA in WGS. We found that contamination is pervasive and can introduce large biases in variant analysis. We showed that these biases can result in hundreds of false positive and negative SNPs, even for samples with slight contamination. Studies investigating complex biological traits from sequencing data can be completely biased if contamination is neglected during the bioinformatic analysis, and we demonstrate that removing contaminant reads with a taxonomic classifier permits more accurate variant calling. We used both real and simulated data to evaluate and implement reliable, contamination-aware analysis pipelines. CONCLUSION: As sequencing technologies consolidate as precision tools that are increasingly adopted in the research and clinical context, our results urge for the implementation of contamination-aware analysis pipelines. Taxonomic classifiers are a powerful tool to implement such pipelines.

Pre-Columbian mycobacterial genomes reveal seals as a source of New World human tuberculosis.

Modern strains of Mycobacterium tuberculosis from the Americas are closely related to those from Europe, supporting the assumption that human tuberculosis was introduced post-contact. This notion, however, is incompatible with archaeological evidence of pre-contact tuberculosis in the New World. Comparative genomics of modern isolates suggests that M. tuberculosis attained its worldwide distribution following human dispersals out of Africa during the Pleistocene epoch, although this has yet to be confirmed with ancient calibration points. Here we present three 1,000-year-old mycobacterial genomes from Peruvian human skeletons, revealing that a member of the M. tuberculosis complex caused human disease before contact. The ancient strains are distinct from known human-adapted forms and are most closely related to those adapted to seals and sea lions. Two independent dating approaches suggest a most recent common ancestor for the M. tuberculosis complex less than 6,000 years ago, which supports a Holocene dispersal of the disease. Our results implicate sea mammals as having played a role in transmitting the disease to humans across the ocean.

Cryptic Resistance Mutations Associated With Misdiagnoses of Multidrug-Resistant Tuberculosis.

Understanding why some multidrug-resistant tuberculosis cases are not detected by rapid phenotypic and genotypic routine clinical tests is essential to improve diagnostic assays and advance toward personalized tuberculosis treatment. Here, we combine whole-genome sequencing with single-colony phenotyping to identify a multidrug-resistant strain that had infected a patient for 9 years. Our investigation revealed the failure of rapid testing and genome-based prediction tools to identify the multidrug-resistant strain. The false-negative findings were caused by uncommon rifampicin and isoniazid resistance mutations. Although whole-genome sequencing data helped to personalize treatment, the patient developed extensively drug-resistant tuberculosis, highlighting the importance of coupling new diagnostic methods with appropriate treatment regimens.

Simplified Model to Survey Tuberculosis Transmission in Countries Without Systematic Molecular Epidemiology Programs.

Systematic molecular/genomic epidemiology studies for tuberculosis surveillance cannot be implemented in many countries. We selected Panama as a model for an alternative strategy. Mycobacterial interspersed repetitive unit-variable-number tandem-repeat (MIRU-VNTR) analysis revealed a high proportion (50%) of Mycobacterium tuberculosis isolates included in 6 clusters (A-F) in 2 provinces (Panama and Colon). Cluster A corresponded to the Beijing sublineage. Whole-genome sequencing (WGS) differentiated clusters due to active recent transmission, with low single-nucleotide polymorphism-based diversity (cluster C), from clusters involving long-term prevalent strains with higher diversity (clusters A, B). Prospective application in Panama of 3 tailored strain-specific PCRs targeting marker single-nucleotide polymorphisms identified from WGS data revealed that 31.4% of incident cases involved strains A-C and that the Beijing strain was highly represented and restricted mainly to Colon. Rational integration of MIRU-VNTR, WGS, and tailored strain-specific PCRs could be a new model for tuberculosis surveillance in countries without molecular/genomic epidemiology programs.

High-resolution mapping of tuberculosis transmission: Whole genome sequencing and phylogenetic modelling of a cohort from Valencia Region, Spain.

BACKGROUND: Whole genome sequencing provides better delineation of transmission clusters in Mycobacterium tuberculosis than traditional methods. However, its ability to reveal individual transmission links within clusters is limited. Here, we used a 2-step approach based on Bayesian transmission reconstruction to (1) identify likely index and missing cases, (2) determine risk factors associated with transmitters, and (3) estimate when transmission happened. METHODS AND FINDINGS: We developed our transmission reconstruction method using genomic and epidemiological data from a population-based study from Valencia Region, Spain. Tuberculosis (TB) incidence during the study period was 8.4 cases per 100,000 people. While the study is ongoing, the sampling frame for this work includes notified TB cases between 1 January 2014 and 31 December 2016. We identified a total of 21 transmission clusters that fulfilled the criteria for analysis. These contained a total of 117 individuals diagnosed with active TB (109 with epidemiological data). Demographic characteristics of the study population were as follows: 80/109 (73%) individuals were Spanish-born, 76/109 (70%) individuals were men, and the mean age was 42.51 years (SD 18.46). We found that 66/109 (61%) TB patients were sputum positive at diagnosis, and 10/109 (9%) were HIV positive. We used the data to reveal individual transmission links, and to identify index cases, missing cases, likely transmitters, and associated transmission risk factors. Our Bayesian inference approach suggests that at least 60% of index cases are likely misidentified by local public health. Our data also suggest that factors associated with likely transmitters are different to those of simply being in a transmission cluster, highlighting the importance of differentiating between these 2 phenomena. Our data suggest that type 2 diabetes mellitus is a risk factor associated with being a transmitter (odds ratio 0.19 [95% CI 0.02-1.10], p < 0.003). Finally, we used the most likely timing for transmission events to study when TB transmission occurred; we identified that 5/14 (35.7%) cases likely transmitted TB well before symptom onset, and these were largely sputum negative at diagnosis. Limited within-cluster diversity does not allow us to extrapolate our findings to the whole TB population in Valencia Region. CONCLUSIONS: In this study, we found that index cases are often misidentified, with downstream consequences for epidemiological investigations because likely transmitters can be missed. Our findings regarding inferred transmission timing suggest that TB transmission can occur before patient symptom onset, suggesting also that TB transmits during sub-clinical disease. This result has direct implications for diagnosing TB and reducing transmission. Overall, we show that a transition to individual-based genomic epidemiology will likely close some of the knowledge gaps in TB transmission and may redirect efforts towards cost-effective contact investigations for improved TB control.