Corrigendum: CARB-ES-19 multicenter study of carbapenemase-producing Klebsiella pneumoniae and Escherichia coli from all Spanish provinces reveals interregional spread of high-risk clones such as ST307/OXA-48 and ST512/KPC-3.

[This corrects the article DOI: 10.3389/fmicb.2022.918362.].

Widespread Detection of Yersiniabactin Gene Cluster and Its Encoding Integrative Conjugative Elements (ICEKp) among Nonoutbreak OXA-48-Producing Klebsiella pneumoniae Clinical Isolates from Spain and the Netherlands.

In this study, we determined the presence of virulence factors in nonoutbreak, high-risk clones and other isolates belonging to less common sequence types associated with the spread of OXA-48-producing Klebsiella pneumoniae clinical isolates from The Netherlands (n = 61) and Spain (n = 53). Most isolates shared a chromosomally encoded core of virulence factors, including the enterobactin gene cluster, fimbrial fim and mrk gene clusters, and urea metabolism genes (ureAD). We observed a high diversity of K-Locus and K/O loci combinations, KL17 and KL24 (both 16%), and the O1/O2v1 locus (51%) being the most prevalent in our study. The most prevalent accessory virulence factor was the yersiniabactin gene cluster (66.7%). We found seven yersiniabactin lineages-ybt 9, ybt 10, ybt 13, ybt 14, ybt 16, ybt 17, and ybt 27-which were chromosomally embedded in seven integrative conjugative elements (ICEKp): ICEKp3, ICEKp4, ICEKp2, ICEKp5, ICEKp12, ICEKp10, and ICEKp22, respectively. Multidrug-resistant lineages-ST11, ST101, and ST405-were associated with ybt 10/ICEKp4, ybt 9/ICEKp3, and ybt 27/ICEKp22, respectively. The fimbrial adhesin kpi operon (kpiABCDEFG) was predominant among ST14, ST15, and ST405 isolates, as well as the ferric uptake system kfuABC, which was also predominant among ST101 isolates. No convergence of hypervirulence and resistance was observed in this collection of OXA-48-producing K. pneumoniae clinical isolates. Nevertheless, two isolates, ST133 and ST792, were positive for the genotoxin colibactin gene cluster (ICEKp10). In this study, the integrative conjugative element, ICEKp, was the major vehicle for yersiniabactin and colibactin gene clusters spreading. IMPORTANCE Convergence of multidrug resistance and hypervirulence in Klebsiella pneumoniae isolates has been reported mostly related to sporadic cases or small outbreaks. Nevertheless, little is known about the real prevalence of carbapenem-resistant hypervirulent K. pneumoniae since these two phenomena are often separately studied. In this study, we gathered information on the virulent content of nonoutbreak, high-risk clones (i.e., ST11, ST15, and ST405) and other less common STs associated with the spread of OXA-48-producing K. pneumoniae clinical isolates. The study of virulence content in nonoutbreak isolates can help us to expand information on the genomic landscape of virulence factors in K. pneumoniae population by identifying virulence markers and their mechanisms of spread. Surveillance should focus not only on antimicrobial resistance but also on virulence characteristics to avoid the spread of multidrug and (hyper)virulent K. pneumoniae that may cause untreatable and more severe infections.

First confirmation of importation and transmission in Spain of the newly identified SARS-CoV-2 B.1.1.7 variant.

INTRODUCTION: A newly identified SARS-CoV-2 variant, VOC202012/01 originating lineage B.1.1.7, recently emerged in the United Kingdom. The rapid spread in the UK of this new variant has caused other countries to be vigilant. MATERIAL AND METHODS: We based our initial screening of B.1.1.7 on the dropout of the S gene signal in the TaqPath assay, caused by the 69/70 deletion. Subsequently, we confirmed the B.1.1.7 candidates by whole genome sequencing. RESULTS: We describe the first three imported cases of this variant from London to Madrid, subsequent post-arrival household transmission to three relatives, and the two first cases without epidemiological links to UK. One case required hospitalization. In all cases, drop-out of gene S was correctly associated to the B.1.1.7 variant, as all the corresponding sequences carried the 17 lineage-marker mutations. CONCLUSION: The first identifications of the SARS-CoV-2 B.1.1.7 variant in Spain indicate the role of independent introductions from the UK coexisting with post-arrival transmission in the community, since the early steps of this new variant in our country.

CARB-ES-19 Multicenter Study of Carbapenemase-Producing Klebsiella pneumoniae and Escherichia coli From All Spanish Provinces Reveals Interregional Spread of High-Risk Clones Such as ST307/OXA-48 and ST512/KPC-3.

Objectives: CARB-ES-19 is a comprehensive, multicenter, nationwide study integrating whole-genome sequencing (WGS) in the surveillance of carbapenemase-producing K. pneumoniae (CP-Kpn) and E. coli (CP-Eco) to determine their incidence, geographical distribution, phylogeny, and resistance mechanisms in Spain. Methods: In total, 71 hospitals, representing all 50 Spanish provinces, collected the first 10 isolates per hospital (February to May 2019); CPE isolates were first identified according to EUCAST (meropenem MIC > 0.12 mg/L with immunochromatography, colorimetric tests, carbapenem inactivation, or carbapenem hydrolysis with MALDI-TOF). Prevalence and incidence were calculated according to population denominators. Antibiotic susceptibility testing was performed using the microdilution method (EUCAST). All 403 isolates collected were sequenced for high-resolution single-nucleotide polymorphism (SNP) typing, core genome multilocus sequence typing (cgMLST), and resistome analysis. Results: In total, 377 (93.5%) CP-Kpn and 26 (6.5%) CP-Eco isolates were collected from 62 (87.3%) hospitals in 46 (92%) provinces. CP-Kpn was more prevalent in the blood (5.8%, 50/853) than in the urine (1.4%, 201/14,464). The cumulative incidence for both CP-Kpn and CP-Eco was 0.05 per 100 admitted patients. The main carbapenemase genes identified in CP-Kpn were bla OXA-48 (263/377), bla KPC-3 (62/377), bla VIM-1 (28/377), and bla NDM-1 (12/377). All isolates were susceptible to at least two antibiotics. Interregional dissemination of eight high-risk CP-Kpn clones was detected, mainly ST307/OXA-48 (16.4%), ST11/OXA-48 (16.4%), and ST512-ST258/KPC (13.8%). ST512/KPC and ST15/OXA-48 were the most frequent bacteremia-causative clones. The average number of acquired resistance genes was higher in CP-Kpn (7.9) than in CP-Eco (5.5). Conclusion: This study serves as a first step toward WGS integration in the surveillance of carbapenemase-producing Enterobacterales in Spain. We detected important epidemiological changes, including increased CP-Kpn and CP-Eco prevalence and incidence compared to previous studies, wide interregional dissemination, and increased dissemination of high-risk clones, such as ST307/OXA-48 and ST512/KPC-3.

First importations of SARS-CoV-2 P.1 and P.2 variants from Brazil to Spain and early community transmission.

INTRODUCTION: SARS-CoV-2variants of concern (VOC) have been described in the UK (B.1.1.7), South Africa (B.1.351) and Brazil (P.1). Among them, the most scarce information has been obtained from the P.1 variant and more data on its global presence and about its spreading dynamics are needed. METHODS: Whole genome sequencing was performed prospectively on travellers arriving from Brazil and on a random selection of SARS-CoV-2 positive cases from our population. RESULTS: In this study we report the first SARS-CoV-2 P.1 and P.2 variants exported from Brazil to Spain. The case infected with the P.1 variant, who had only stayed in Rio de Janeiro, required hospitalisation. The two P.2 cases remained asymptomatic. A wider distribution for P.1 variant beyond the Brazilian Amazonia should be considered. The exportation of the P.2 variant, carrying the E484K mutation, deserves attention. One month after the first description of P.1 and P.2 importations from Brazil to Madrid, these variants were identified circulating in the community, in cases without a travel history, and involved in household transmissions CONCLUSION: Whole genome sequencing of SARS-CoV-2 positive travellers arriving from Brazil allowed us to identify the first importations of P.1 and P.2 variants to Spain and their early community transmission.

Insights into the Complexity of a Dormant Mycobacterium tuberculosis Cluster Once Transmission Is Resumed.

Genotyping tools help identify the complexity in Mycobacterium tuberculosis transmission clusters. We carried out a thorough analysis of the epidemiological and bacteriological complexity of a cluster in Almería, Spain. The cluster, initially associated with Moroccan migrants and with no secondary cases identified in 4 years, then reappeared in Spanish-born individuals. In one case, two Mycobacterium tuberculosis clonal variants were identified. We reanalyzed the cluster, supported by the characterization of multiple cultured isolates and respiratory specimens, whole-genome sequencing, and epidemiological case interviews. Our findings showed that the cluster, which was initially thought to have restarted activity with just a single case harboring a small degree of within-host diversity, was in fact currently growing due to coincidental reactivation of past exposures, with clonal diversity transmitted throughout the cluster. In one case, within-host diversity was amplified, probably due to prolonged diagnostic delay. IMPORTANCE The precise study of the dynamics of tuberculosis transmission in socio-epidemiologically complex scenarios may require more thorough analysis than the standard molecular epidemiology strategies. Our study illustrates the epidemiological and bacteriological complexity present in a transmission cluster in a challenging epidemiological setting with a high proportion of migrant cases. The combination of whole-genome sequencing, refined and refocused epidemiological interviews, and in-depth analysis of the bacterial composition of sputa and cultured isolates was crucial in order to correctly reinterpret the true nature of this cluster. Our global approach allowed us to reinterpret correctly the unnoticed epidemiological and bacteriological complexity involved in the Mycobacterium tuberculosis transmission event under study, which had been overlooked by the usual molecular epidemiology approaches.

Nosocomial Transmission of SARS-CoV-2 Involving Vaccinated Health Care Workers.

COVID-19 vaccination has proven to be effective at preventing symptomatic disease but there are scarce data to fully understand whether vaccinated individuals can still behave as SARS-CoV-2 transmission vectors. Based on viral genome sequencing and detailed epidemiological interviews, we report a nosocomial transmission event involving two vaccinated health care-workers (HCWs) and four patients, one of them with fatal outcome. Strict transmission control measures, as during the prevaccination period, must be kept between HCWs and HCWs-patients in nosocomial settings. IMPORTANCE COVID-19 vaccination has proven to be effective at preventing symptomatic disease. Although some transmission events involving vaccinated cases have also been reported, scarce information is still available to fully understand whether vaccinated individuals may still behave as vectors in SARS-CoV-2 transmission events. Here, we report a SARS-CoV-2 nosocomial transmission event, supported on whole genome sequencing, in early March 2021 involving two vaccinated HCWs and four patients in our institution. Strict transmission control measures between HCWs and HCWs - patients in nosocomial settings must not be relaxed, and should be kept as strictly as during the prevaccination period.

First importations of SARS-CoV-2 P.1 and P.2 variants from Brazil to Spain and early community transmission.

Introduction: SARS-CoV-2variants of concern (VOC) have been described in the UK (B.1.1.7), South Africa (B.1.351) and Brazil (P.1). Among them, the most scarce information has been obtained from the P.1 variant and more data on its global presence and about its spreading dynamics are needed. Methods: Whole genome sequencing was performed prospectively on travellers arriving from Brazil and on a random selection of SARS-CoV-2 positive cases from our population. Results: In this study we report the first SARS-CoV-2 P.1 and P.2 variants exported from Brazil to Spain. The case infected with the P.1 variant, who had only stayed in Rio de Janeiro, required hospitalisation. The two P.2 cases remained asymptomatic. A wider distribution for P.1 variant beyond the Brazilian Amazonia should be considered. The exportation of the P.2 variant, carrying the E484K mutation, deserves attention. One month after the first description of P.1 and P.2 importations from Brazil to Madrid, these variants were identified circulating in the community, in cases without a travel history, and involved in household transmissions. Conclusion: Whole genome sequencing of SARS-CoV-2 positive travellers arriving from Brazil allowed us to identify the first importations of P.1 and P.2 variants to Spain and their early community transmission.

Introducción: Se han descrito «variantes de preocupación¼ (VOC) de SARS-CoV-2 en el Reino Unido (B.1.1.7), Sudáfrica (B.1.351) y Brasil (P.1). Entre ellas, se dispone de información más escasa para la variante P.1 y se necesitan más datos sobre su presencia global y sobre su dinámica de expansión. Métodos: Se realizó secuenciación del genoma completo de forma prospectiva de SARS-CoV-2 en viajeros procedentes de Brasil y en una selección aleatoria de casos positivos de SARS-CoV-2 de nuestra población. Resultados: En este estudio reportamos las primeras variantes de SARS-CoV-2 P.1 y P.2 exportadas desde Brasil a España. El caso infectado por la variante P.1, que solo había permanecido en Río de Janeiro, requirió hospitalización. Los 2 casos de la variante P.2 permanecieron asintomáticos. Se debe considerar una distribución más amplia para la variante P.1 más allá de la Amazonía brasileña. La exportación de la variante P.2, que porta la mutación E484K, merece asimismo atención adicional. Un mes después de la primera descripción de las importaciones de P.1 y P.2 de Brasil a Madrid, se identificaron estas variantes circulando en la comunidad, en casos sin antecedentes de viaje, e implicadas en transmisiones domiciliarias. Conclusión: La secuenciación de genoma completo de viajeros positivos para SARS-CoV-2 procedentes de Brasil nos permitió identificar las primeras importaciones de variantes P.1 y P.2 a España y su transmisión comunitaria precoz.

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.

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 superinfection and reinfection with three different strains.

We report a corona virus disease (COVID-19) case with unprecedented viral complexity. In the first severe episode, two different severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) strains (superinfection) were identified within a week. Three months after discharge, the patient was readmitted and was infected in a nosocomial outbreak with a different strain, suffering a second milder COVID-19 episode.

Detection of Minority Variants and Mixed Infections in Mycobacterium tuberculosis by Direct Whole-Genome Sequencing on Noncultured Specimens Using a Specific-DNA Capture Strategy.

Detection of mixed Mycobacterium tuberculosis (MTB) infections is essential, particularly when resistance mutations are present in minority bacterial populations that may affect patients' disease evolution and treatment. Whole-genome sequencing (WGS) has extended the amount of key information available for the diagnosis of MTB infection, including the identification of mixed infections. Having genomic information at diagnosis for early intervention requires carrying out WGS directly on the clinical samples. However, few studies have been successful with this approach due to the low representation of MTB DNA in sputa. In this study, we evaluated the ability of a strategy based on specific MTB DNA enrichment by using a newly designed capture platform (MycoCap) to detect minority variants and mixed infections by WGS on controlled mixtures of MTB DNAs in a simulated sputum genetic background. A pilot study was carried out with 12 samples containing 98% of a DNA pool from sputa of patients without MTB infection and 2% of MTB DNA mixtures at different proportions. Our strategy allowed us to generate sequences with a quality equivalent to those obtained from culture: 62.5× depth coverage and 95% breadth coverage (for at least 20× reads). Assessment of minority variant detection was carried out by manual analysis and allowed us to identify heterozygous positions up to a 95:5 ratio. The strategy also automatically distinguished mixed infections up to a 90:10 proportion. Our strategy efficiently captures MTB DNA in a nonspecific genetic background, allows detection of minority variants and mixed infections, and is a promising tool for performing WGS directly on clinical samples. IMPORTANCE We present a new strategy to identify mixed infections and minority variants in Mycobacterium tuberculosis by whole-genome sequencing. The objective of the strategy is the direct detection in patient sputum; in this way, minority populations of resistant strains can be identified at the time of diagnosis, facilitating identification of the most appropriate treatment for the patient from the first moment. For this, a platform for capturing M. tuberculosis-specific DNA was designed to enrich the clinical sample and obtain quality sequences.

SARS-CoV-2 B.1.1.7 Decline Is Not Driven by the Introduction of a More Successful Variant.

The SARS-CoV-2 variant of concern (VOC) Delta (B.617.2 lineage) displaced the predominant VOC Alpha (B.1.1.7 lineage) in the United Kingdom. In Madrid, recent start of the decline of predominant VOC Alpha suggested an equivalent phenomenon. However, 11 different variants, none overrepresented in frequency, occupied progressively over a period of 7 weeks the niche previously dominated by VOC Alpha. Only after these 7 weeks, VOC Delta started to emerge. Viral competition due to the entry of VOC Delta is not the major force driving the start of VOC Alpha decline in Madrid. IMPORTANCE Our data indicate that the dynamics of SARS-CoV-2 VOCs turnover in our setting differ from those proposed for other countries. A systematic genomic analysis, updated on a weekly basis, of representative randomly selected samples of SARS-CoV-2 circulating variants allowed us to define a lapse of 7 weeks between the start of VOC Alpha decline and the final emergence of VOC Delta. During this period, VOC Alpha showed a sustained decline, while 11 VOCs, variants of interest (VOIs), and other identified variants, none overrepresented, occupied the niche left by VOC Alpha. Only after these 7 weeks, emergence of VOC Delta occurred, indicating that viral competition involving VOC Delta was not the exclusive direct driving force behind the starting of VOC Alpha decline.

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.

Different Within-Host Viral Evolution Dynamics in Severely Immunosuppressed Cases with Persistent SARS-CoV-2.

A successful Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variant, B.1.1.7, has recently been reported in the UK, causing global alarm. Most likely, the new variant emerged in a persistently infected patient, justifying a special focus on these cases. Our aim in this study was to explore certain clinical profiles involving severe immunosuppression that may help explain the prolonged persistence of viable viruses. We present three severely immunosuppressed cases (A, B, and C) with a history of lymphoma and prolonged SARS-CoV-2 shedding (2, 4, and 6 months), two of whom finally died. Whole-genome sequencing of 9 and 10 specimens from Cases A and B revealed extensive within-patient acquisition of diversity, 12 and 28 new single nucleotide polymorphisms, respectively, which suggests ongoing SARS-CoV-2 replication. This diversity was not observed for Case C after analysing 5 sequential nasopharyngeal specimens and one plasma specimen, and was only observed in one bronchoaspirate specimen, although viral viability was still considered based on constant low Ct values throughout the disease and recovery of the virus in cell cultures. The acquired viral diversity in Cases A and B followed different dynamics. For Case A, new single nucleotide polymorphisms were quickly fixed (13-15 days) after emerging as minority variants, while for Case B, higher diversity was observed at a slower emergence: fixation pace (1-2 months). Slower SARS-CoV-2 evolutionary pace was observed for Case A following the administration of hyperimmune plasma. This work adds knowledge on SARS-CoV-2 prolonged shedding in severely immunocompromised patients and demonstrates viral viability, noteworthy acquired intra-patient diversity, and different SARS-CoV-2 evolutionary dynamics in persistent cases.

Overlapping of Independent SARS-CoV-2 Nosocomial Transmissions in a Complex Outbreak.

SARS-CoV-2 nosocomial outbreaks in the first COVID-19 wave were likely associated with a shortage of personal protective equipment and scarce indications on control measures. Having covered these limitations, updates on current SARS-CoV-2 nosocomial outbreaks are required. We carried out an in-depth analysis of a 27-day nosocomial outbreak in a gastroenterology ward in our hospital, potentially involving 15 patients and 3 health care workers. Patients had stayed in one of three neighboring rooms in the ward. The severity of the infections in six of the cases and a high fatality rate made the clinicians suspect the possible involvement of a single virulent strain persisting in those rooms. Whole-genome sequencing (WGS) of the strains from 12 patients and 1 health care worker revealed an unexpected complexity. Five different SARS-CoV-2 strains were identified, two infecting a single patient each, ruling out their relationship with the outbreak; the remaining three strains were involved in three independent, overlapping, limited transmission clusters with three, three, and five cases. Whole-genome sequencing was key to understand the complexity of this outbreak. IMPORTANCE We report a complex epidemiological scenario of a nosocomial COVID-19 outbreak in the second wave, based on WGS analysis. Initially, standard epidemiological findings led to the assumption of a homogeneous outbreak caused by a single SARS-CoV-2 strain. The discriminatory power of WGS offered a strikingly different perspective consisting of five introductions of different strains, with only half of them causing secondary cases in three independent overlapping clusters. Our study exemplifies how complex the SARS-CoV-2 transmission in the nosocomial setting during the second COVID-19 wave occurred and leads to extending the analysis of outbreaks beyond the initial epidemiological assumptions.

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.

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.