Search

14 páginas, 6 figuras. Supplementary information is available for this paper at https://doi.org/10.1038/s41467-020-15832-6. Sequence data are deposited on NCBI GEO database with the accession code GSE138580, 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, This work was financed by FCT-Fundação para a Ciência e a Tecnologia/ Ministério da Ciência, Tecnologia e Inovação grant POCI-01-0145-FEDER-028955 (to M.S.) and by the Northern Portugal Regional Operational Program (NORTE 2020), under the Portugal 2020 Partnership Agreement, through the European Regional Development Fund (FEDER) (NORTE-01-0145-FEDER-000013, to M.I.V., F.R., A.G.C., and N.S.O.). I.C.acknowledges the support of Ministerio de Ciencia, Innovación y Universidades (SAF2016-77346-R) and the European Research Council (638553-TB-ACCELERATE). H.N.B. acknowledges the support of Bolsa D. Manuel de Mello and of the Portuguese Society for Pneumology; A.B. and M.S. were also recipients of an International Exchanges Grant from the Royal Society. J.S. is funded by a research fellow NORTE-01-0145-FEDER-000012; B.C. and K.L.F. are funded by FCT Ph.D scholarships SFRH/BD/114403/2016 and SFRH/BD/114405/2016, respectively; M.I.V. is funded by FCT through DL 57/2016 (CRP) and MS through Estimulo Individual ao Emprego Científico. We thank the excellent support from the i3S scientific platforms, namely Animal facility, Advanced Light Microscopy, and BioSciences Screening, member of the national infrastructure PPBI-Portuguese Platform of Bioimaging (PPBI-POCI-01-0145-FEDER-022122)and Genomics (GenCore) part of the GenomePT project (POCI-01-0145-FEDER022184), supported by COMPETE 2020 - Operational Programme for Competitiveness and Internationalisation (POCI), Lisboa Portugal Regional Operational Programme (Lisboa2020), Algarve Portugal Regional Operational Programme (CRESC Algarve2020), under the PORTUGAL 2020 Partnership Agreement, through the European Regional Development Fund (ERDF), and by Fundação para a Ciência e a Tecnologia (FCT).

12 páginas, 4 figuras, 1 tabla, 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., This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programs 638553 (TB-ACCELERATE), SAF2016-77346-R from Ministerio de Economía y Competitividad (Spanish Government) and AICO/2018/113 from Generalitat Valenciana (to I.C.). M.M.-M. is recipient of a FPI grant from the Ministerio de Economía y Competitividad (Spanish Government, code BES-2017-079656). V.F. was recipient of a post-doctoral research grant from the Ministerio de Economía y Competitividad (Spanish Government, code FPDI-2013-18757). Action cofinanced by the European Union through the Operational Program of European Regional Development Fund (ERDF) of Valencia Region (Spain) 2014-2020.

11 páginas, 6 figuras., The Mycobacterium tuberculosis complex (MTBC) members display different host-specificities and virulence phenotypes. Here, we have performed a comprehensive RNAseq and methylome analysis of the main clades of the MTBC and discovered unique transcriptional profiles. The majority of genes differentially expressed between the clades encode proteins involved in host interaction and metabolic functions. A significant fraction of changes in gene expression can be explained by positive selection on single mutations that either create or disrupt transcriptional start sites (TSS). Furthermore, we show that clinical strains have different methyltransferases inactivated and thus different methylation patterns. Under the tested conditions, differential methylation has a minor direct role on transcriptomic differences between strains. However, disruption of a methyltransferase in one clinical strain revealed important expression differences suggesting indirect mechanisms of expression regulation. Our study demonstrates that variation in transcriptional profiles are mainly due to TSS mutations and have likely evolved due to differences in host characteristics., This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programmes 638553 (TB-ACCELERATE to I.C.) and 637730 (MtbTransReg to T.C.). M.B. was funded by the National Institute of Health grant AI119573. W.R.J. was funded by the National Institutes of Health grant AI26170. In addition, this work was funded by projects SAF2016-77346-R from Ministerio de Economía y Competitividad (Spanish Government) and AICO/2018/113 from Generalitat Valenciana (to I.C.), BFU2014-58656-R, BFU2017-89594-R from Ministerio de Economía y Competitividad (Spanish Government), PROMETEO/2016/122 from Generalitat Valenciana (to FGC), Wellcome Trust grant 098051 to the Sanger Institute. A.C.O. is recipient of a FPU fellowship from Ministerio de Educación y Ciencia FPU13/00913 (Spanish Government).

12 páginas, 2 figuras, 1 tabla., 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., This research was funded by Instituto de Salud Carlos III (Ref COV20/00140: SeqCOVID—Consorcio para la epidemiología genómica de SARS-CoV-2 en España) and by Consejo Superior de Investigaciones Científicas (CSIC) (PTI Salud Global). Miguel Servet Contracts CP15/00075 and CPII20/00001) to L.P.-L.

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

12 páginas, 5 figuras., Widespread and frequent resistance to the second-line tuberculosis (TB) medicine streptomycin, suggests ongoing transmission of low fitness cost streptomycin resistance mutations. To investigate this hypothesis, we studied a cohort of 681 individuals from a TB epidemic in Portugal. Whole-genome sequencing (WGS) analyses were combined with phenotypic growth studies in culture media and in mouse bone marrow derived macrophages. Streptomycin resistance was the most frequent resistance in the cohort accounting for 82.7% (n = 67) of the resistant Mycobacterium tuberculosis isolates. WGS of 149 clinical isolates identified 13 transmission clusters, including three clusters containing only streptomycin resistant isolates. The biggest cluster was formed by eight streptomycin resistant isolates with a maximum of five pairwise single nucleotide polymorphisms of difference. Interestingly, despite their genetic similarity, these isolates displayed different resistance levels to streptomycin, as measured both in culture media and in infected mouse bone marrow derived macrophages. The genetic bases underlying this phenotype are a combination of mutations in gid and other genes. This study suggests that specific streptomycin resistance mutations were transmitted in the cohort, with the resistant isolates evolving at the cluster level to allow low-to-high streptomycin resistance levels without a significative fitness cost. This is relevant not only to better understand transmission of streptomycin resistance in a clinical setting dominated by Lineage 4 M. tuberculosis infections, but mainly because it opens new prospects for the investigation of selection and spread of drug resistance in general., This work has been funded by of the projects NORTE-01-0145- FEDER-000039, NORTE-01-0145-FEDER-000013, and NORTE01-0145-FEDER-00002, supported by Norte Portugal Regional Operational Program (NORTE 2020), under the PORTUGAL 2020 Partnership Agreement, through the European Regional Development Fund (ERDF). It was also funded through the Foundation for Science and Technology (FCT) – project FCT IF/00474/2014, UIDB/50026/2020 and UIDP/50026/2020. DR and CM were funded by FCT Ph.D. scholarships (grant numbers SFRH/BD/135422/2017 and SFRH/BD/132797/2017, respectively). DR was partially funded by the PGCD – Graduate Program Science for Development. MS is supported by FCT through the Estimulo ao Emprego Científico.

9 páginas, 1 figurra, 3 tablas. Te online version contains supplementary material available at https://doi.org/10.1038/s41598-020-80919-5., Whole genome sequencing (WGS) has been shown to be superior to traditional procedures of genotyping in tuberculosis (TB), nevertheless, reports of its use in drug resistant TB (DR-TB) isolates circulating in Mexico, are practically unknown. Considering the above the main of this work was to identify and characterize the lineages and genomic transmission clusters present in 67 DR-TB isolates circulating in southeastern Mexico. The results show the presence of three major lineages: L1 (3%), L2 (3%) and L4 (94%), the last one included 16 sublineages. Sublineage 4.1.1.3 (X3) was predominant in 18 (27%) of the isolates, including one genomic cluster, formed by eleven multidrug resistant isolates and sharing the SIT 3278, which seems to be restricted to Mexico. By the use of WGS, it was possible to identify the high prevalence of L4 and a high number of sublineages circulating in the region, also was recognized the presence of a novel X3 sublineage, formed exclusively by multidrug resistant isolates and with restrictive circulation in Mexico for at least the past 17 years., AC-JR and CFM-M were CONACYT fellows of the Health Sciences Master of the Health Sciences Institute of Universidad Veracruzana. RZ-C and IC were funded by i-COOP-20017-COOPB2032.

16 p´ginas, 12 figuras, 1 tabla. The authors confirm all supporting data, code and protocols have been provided within the article or through supplementary data files. Supplementary data files can be found at 10.6084/m9.figshare.14562321., The occurrence of multiple strains of a bacterial pathogen such as M. tuberculosis or C. difficile within a single human host, referred to as a mixed infection, has important implications for both healthcare and public health. However, methods for detecting it, and especially determining the proportion and identities of the underlying strains, from WGS (whole-genome sequencing) data, have been limited. In this paper we introduce SplitStrains, a novel method for addressing these challenges. Grounded in a rigorous statistical model, SplitStrains not only demonstrates superior performance in proportion estimation to other existing methods on both simulated as well as real M. tuberculosis data, but also successfully determines the identity of the underlying strains. We conclude that SplitStrains is a powerful addition to the existing toolkit of analytical methods for data coming from bacterial pathogens and holds the promise of enabling previously inaccessible conclusions to be drawn in the realm of public health microbiology, This work has been funded in part by a CANSSI Collaborative Research Team grant, 'Statistical methods for challenging problems in public health microbiology' and a Genome Canada grant, 'Machine Learning Methods to Predict Drug Resistance in Pathogenic Bacteria'. LC acknowledges funding from a Sloan Foundation fellowship (FG-2016–6392) and the MRC Centre for Global Infectious Disease Analysis (reference MR/R015600/1), jointly funded by the UK Medical Research Council (MRC) and the UK Foreign, Commonwealth and Development Office (FCDO), under the MRC/FCDO Concordat agreement, and is part of the EDCTP2 programme supported by the European Union. ML acknowledges funding from a NSERC Discovery grant.

Following its emergence in late 2019, the spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)1,2 has been tracked via phylogenetic analysis of viral genome sequences in unprecedented detail3–5. While the virus spread globally in early 2020 before borders closed, intercontinental travel has since been greatly reduced. However, within Europe travel resumed in the summer of 2020. Here we report on a novel SARS-CoV-2 variant, 20E (EU1), that emerged in Spain in early summer, and subsequently spread across Europe. We find no evidence of 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 20E (EU1) was introduced hundreds of times to European countries by summertime travelers, likely undermining local efforts to keep SARS-CoV-2 cases low. Our results demonstrate how a variant can rapidly become dominant even in absence of a substantial transmission advantage in favorable epidemiological settings. Genomic surveillance is critical to understanding how travel can impact SARS-CoV-2 transmission, and thus for informing future containment strategies as travel resumes., This work was supported by the Swiss National Science Foundation (SNSF) through grant numbers 31CA30 196046 (to RAN, EBH, CLA), 31CA30 196267 (to TS), European Union’s Horizon 2020 research and innovation programme - project EpiPose (No 101003688) (MLR, CLA), core funding by the University of Basel and ETH Zürich, the National Institute of General Medical Sciences (R01GM120553 to DV), the National Institute of Allergy and Infectious Diseases (DP1AI158186 and HHSN272201700059C to DV), a Pew Biomedical Scholars Award (DV), an Investigators in the Pathogenesis of Infectious Disease Awards from the Burroughs Wellcome Fund (DV and JDB), a Fast Grants (DV), and NIAID grants R01AI141707 (JDB) and F30AI149928 (KHDC). SeqCOVID-SPAIN is funded by the Instituto de Salud Carlos III project COV20/00140, Spanish National Research Council and ERC StG 638553 to IC and BFU2017-89594R from MICIN to FGC. JDB is an Investigator of the Howard Hughes Medical Institute.

11 páginas, 5 figuras, 1 tabla, Polyclonal infections occur when at least two unrelated strains of the same pathogen are detected in an individual. This has been linked to worse clinical outcomes in tuberculosis, as undetected strains with different antibiotic resistance profiles can lead to treatment failure. Here, we examine the amount of polyclonal infections in sputum and surgical resections from patients with tuberculosis in the country of Georgia. For this purpose, we sequence and analyse the genomes of Mycobacterium tuberculosis isolated from the samples, acquired through an observational clinical study (NCT02715271). Access to the lung enhanced the detection of multiple strains (40% of surgery cases) as opposed to just using a sputum sample (0-5% in the general population). We show that polyclonal infections often involve genetically distant strains and can be associated with reversion of the patient's drug susceptibility profile over time. In addition, we find different patterns of genetic diversity within lesions and across patients, including mutational signatures known to be associated with oxidative damage; this suggests that reactive oxygen species may be acting as a selective pressure in the granuloma environment. Our results support the idea that the magnitude of polyclonal infections in high-burden tuberculosis settings is underestimated when only testing sputum samples., The authors were supported by projects SAF2016-77346-R and PID2019-104477RB-I00 awarded to IC and the grant BES-2017-079656 awarded to MM by the Spanish Ministry of Economy and Competitiveness and Ministry of Science, the ERC project 638553-TBACCELERATE awarded to IC, Spanish Government-FEDER Funds through CV contract CPII18/00031 and grant PI16/01511, and Generalitat Valencia Grant to I.C. (code PROMETEO/2020/012). The grant providers played no part in study design, data collection, and analysis, or the preparation of the manuscript.

SeqCOVID-Spain consortium: Álvaro Chiner-Oms, Irving Cancino-Muñoz, Mariana G. López, Manuela Torres-Puente, Inmaculada Gómez-Navarro, Santiago Jiménez-Serrano, Jordi Pérez-Tur, Darío García de Viedma, Laura Pérez-Lago, Marta Herranz, Jon Sicilia, Pilar Catalán-Alonso, Julia Suárez González, Patricia Muñoz, Mireia Coscolla, Paula Ruiz-Rodríguez, Fernando González-Candelas, Iñaki Comas, Lidia Ruiz-Roldán, María Alma Bracho, Neris García-González, Llúcia Martínez Priego, Inmaculada Galán-Vendrell, Paula Ruiz-Hueso, Griselda De Marco, María Loreto Ferrús-Abad, Sandra Carbó-Ramírez, Giuseppe D’Auria, Galo Adrian Goig, Juan Alberola, Jose Miguel Nogueira, Juan José Camarena, David Navarro, Eliseo Albert, Ignacio Torres, Maitane Aranzamendi Zaldumbide, Óscar Martínez Expósito, Nerea Antona Urieta, María de Toro, María Pilar Bea-Escudero, Jose Antonio Boga, Cristian Castelló-Abietar, Susana Rojo-Alba, Marta Elena Álvarez-Argüelles, Santiago Melón, Elisa Martró, Antoni E. Bordoy, Anna Not, Adrián Antuori, Anabel Fernández-Navarro, Andrés Canut-Blasco, Silvia Hernáez Crespo, Maria Luz Cordón Rodríguez, Maria Concepción Lecaroz Agara, Carmen Gómez-González, Amaia Aguirre-Quiñonero, José Israel López-Mirones, Marina Fernández-Torres, Maria Rosario Almela-Ferrer, Ana Carvajal, Juan Miguel Fregeneda-Grandes, Héctor Argüello, Gustavo Cilla Eguiluz, Milagrosa Montes Ros, Luis Piñeiro Vázquez, Ane Sorarrain, José María Marimón, José J. Costa-Alcalde, Rocío Trastoy, Gema Barbeito Castiñeiras, Amparo Coira, María Luisa Pérez del Molino, Antonio Aguilera, Begoña Palop-Borrás, Inmaculada de Toro Peinado, Maria Concepción Mediavilla Gradolph, Mercedes Pérez-Ruiz, Mirian Fernández-Alonso, Jose Luis del Pozo, Oscar González-Recio, Mónica Gutiérrez-Rivas, Jovita Fernández-Pinero, Miguel Ángel Jiménez Clavero, Begoña Fuster Escrivá, Concepción Gimeno Cardona, María Dolores Ocete Mochón, Rafael Medina-Gonzalez, José Antonio Lepe, Verónica González Galán, Ángel Rodríguez-Villodres, Nieves Gonzalo Jiménez, Jordi Reina, Carla López-Causapé, Maria Dolores Gómez-Ruiz, Eva M. Gonzalez-Barbera, José Luis López-Hontangas, Vicente Martín, Antonio J. Molina, Tania Fernandez-Villa, Ana Milagro Beamonte, Nieves Felisa Martínez-Cameo, Yolanda Gracia-Grataloup, Rosario Moreno-Muñoz, Maria Dolores Tirado Balaguer, José María Navarro-Marí, Irene Pedrosa-Corral, Sara Sanbonmatsu-Gámez, Antonio Oliver, Mónica Parra Grande, Bárbara Gómez Alonso, Francisco José Arjona Zaragozí, Maria Carmen Pérez González, Francisco Javier Chamizo López, Ana Bordes-Benítez, Núria Rabella, Ferran Navarro, Elisenda Miró, Antonio Rezusta, Alexander Tristancho, Encarnación Simarro Córdoba, Julia Lozano-Serra, Lorena Robles Fonseca, Álex Soriano, Francisco Javier Roig Sena, Hermelinda Vanaclocha Luna, Isabel Sanmartín, Daniel García-Souto, Ana Pequeño-Valtierra, Jose M. C. Tubio, Javier Temes, Jorge Rodríguez-Castro, Martín Santamarina García, Manuel Rodríguez-Iglesias, Fátima Galán-Sanchez, Salud Rodríguez-Pallares, José Manuel Azcona-Gutiérrez, Miriam Blasco-Alberdi, Alfredo Mayor, Alberto L. García-Basteiro, Gemma Moncunill, Carlota Dobaño, Pau Cisteró, Oriol Mitjà, Camila González-Beiras, Martí Vall-Mayans, Marc Corbacho-Monné, Andrea Alemany, Cristina Muñoz-Cuevas, Guadalupe Rodríguez-Rodríguez, Rafael Benito, Sonia Algarate, Jessica Bueno, Andrea Vergara-Gómez, Miguel J. Martínez, Jordi Vila, Elisa Rubio, Aida Peiró-Mestres, Jessica Navero-Castillejos, David Posada, Diana Valverde, Nuria Estévez, Iria Fernández-Silva, Loretta de Chiara, Pilar Gallego-García, Nair Varela, Ulises Gómez-Pinedo, Mónica Gozalo-Margüello, Maria Eliecer Cano García, José Manuel Méndez-Legaza, Jesus Rodríguez-Lozano, María Siller, Daniel Pablo-Marcos, Maria Montserrat Ruiz-García, Antonio Galiana, Judith Sánchez-Almendro, Maria Isabel Gascón Ros, Cristina Juana Torregrosa-Hetland, Eva María Pastor Boix, Paloma Cascales Ramos, Pedro Luis Garcinuño Enríquez, Salvador Raga Borja, Julia González Cantó, Olalla Martínez Macias, Adolfo de Salazar, Laura Viñuela González, Natalia Chueca, Federico García, Cristina Gómez-Camarasa, Amparo Farga Martí, Rocío Falcón, Victoria Domínguez-Márquez, Anna M. Planas, Israel Fernández-Cádenas, Maria Ángeles Marcos, Carmen Ezpeleta, Ana Navascués, Ana Miqueleiz Zapatero, Manuel Segovia, Antonio Moreno-Docón, Esther Viedma, Raúl Recio Martínez, Irene Muñoz-Gallego, Sara Gonzalez-Bodi, Maria Dolores Folgueira, Jesús Mingorance, Elias Dahdouh, Fernando Lázaro-Perona, María Rodríguez-Tejedor, María Pilar Romero-Gómez, Julio García-Rodríguez, Juan Carlos Galán, Mario Rodríguez-Dominguez, Laura Martínez-García, Melanie Abreu Di Berardino, Manuel Ponce-Alonso, Jose Maria González-Alba, Ivan Sanz-Muñoz, Diana Pérez San José, Maria Gil Fortuño, Juan B. Bellido-Blasco, Alberto Yagüe Muñoz, Noelia Hernández Pérez, Helena Buj Jordá, Óscar Pérez Olaso, Alejandro González Praetorius, Nora Mariela Martínez Ramírez, Aida Ramírez Marinero, Eduardo Padilla León, Alba Vilas Basil, Mireia Canal Aranda, Albert Bernet Sánchez, Alba Bellés Bellés, Eric López González, Iván Prats Sánchez, Mercè García-González, Miguel José Martínez-Lirola, Manuel Ángel Rodríguez Maresca, Maria Teresa Cabezas Fernández, María Eugenia Carrillo Gil, Maria Paz Ventero Martín, Carmen Molina Pardines, Nieves Orta Mira, María Navarro Cots, Inmaculada Vidal Catalá, Isabel García Nava, Soledad Illescas Fernández-Bermejo, José Martínez-Alarcón, Marta Torres-Narbona, Cristina Colmenarejo, Lidia García-Agudo, Jorge A. Pérez García, Martín Yago López, María Ángeles Goberna Bravo, Victoria Simón García, Gonzalo Llop Furquet, Agustín Iranzo Tatay, Sandra Moreno-Marro, Noelia Lozano Rodríguez, Amparo Broseta Tamarit, Juan José Badiola Díez, Amparo Martínez-Ramírez, Ana Dopazo, Sergio Callejas, Alberto Benguría, Begoña Aguado, Antonio Alcamí, Marta Bermejo Bermejo, Ricardo Ramos-Ruíz, Víctor Manuel Fernández Soria, Fernando Simón Soria & Mercedes Roig Cardells, 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, This work was mainly funded by the Instituto de Salud Carlos III project COV20/00140, with additional funding by Spanish National Research Council project CSIC-COV19-021, Ministerio de Ciencia project PID2019-104477RB-100, ERC StG 638553 and ERC CoG 101001038 to I.C., and BFU2017-89594R to F.G.C. M.C. is supported by Ramón y Cajal program from Ministerio de Ciencia and grants RTI2018-094399-A-I00 and Generalitat Valenciana (Regional Government) project SEJI/2019/011. We gratefully acknowledge Hospital Universitari Vall d’Hebron, Instituto de Salud Carlos III, IrsiCaixa AIDS Research Lab and all the international researchers and institutions that submitted sequenced SARS-CoV-2 genomes to the GISAID’s EpiCov Database (Supplementary Table 1), as an important part of our analyses has been made possible by the sharing of their work. We also thank Unidad de Bioinformática y Estadística, Centro de Investigación Príncipe Felipe, for allowing us to use the Computer Cluster to perform some of the bioinformatic analysis.

The COVID-19 pandemic has shaken the world since the beginning of 2020. Spain is among the European countries with the highest incidence of the disease during the first pandemic wave. We established a multidisciplinar consortium to monitor and study the evolution of the epidemic, with the aim of contributing to decision making and stopping rapid spreading across the country. We present the results for 2170 sequences from the first wave of the SARS-Cov-2 epidemic in Spain and representing 12% of diagnosed cases until 14th March. This effort allows us to document at least 500 initial introductions, between early February-March from multiple international sources. Importantly, we document the early raise of two dominant genetic variants in Spain (Spanish Epidemic Clades), named SEC7 and SEC8, likely amplified by superspreading events. In sharp contrast to other non-Asian countries those two variants were closely related to the initial variants of SARS-CoV-2 described in Asia and represented 40% of the genome sequences analyzed. The two dominant SECs were widely spread across the country compared to other genetic variants with SEC8 reaching a 60% prevalence just before the lockdown. Employing Bayesian phylodynamic analysis, we inferred a reduction in the effective reproductive number of these two SECs from around 2.5 to below 0.5 after the implementation of strict public-health interventions in mid March. The effects of lockdown on the genetic variants of the virus are reflected in the general replacement of preexisting SECs by a new variant at the beginning of the summer season. Our results reveal a significant difference in the genetic makeup of the epidemic in Spain and support the effectiveness of lockdown measures in controlling virus spread even for the most successful genetic variants. Finally, earlier control of SEC7 and particularly SEC8 might have reduced the incidence and impact of COVID-19 in our country., This work was funded by the Instituto de Salud Carlos III project COV20/00140, Spanish 593 National Research Council project CSIC-COV19-021 and ERC StG 638553 to IC, and BFU2017- 594 89594R to FGC. MC is supported by Ramón y Cajal program from Ministerio de Ciencia and 595 grants RTI2018-094399-A-I00 and SEJI/2019/011. 596 We gratefully acknowledge Hospital Universitari Vall d'Hebron, Instituto de Salud Carlos 597 III, IrsiCaixa AIDS Research Lab and all the international researchers and institutions that 598 submitted sequenced SARS-CoV-2 genomes to the GISAID’s EpiCov™ Database, as an 599 important part of our analyses have been made possible by the share of their work.

A variant of SARS-CoV-2 emerged in early summer 2020, presumably in Spain, and has since spread to multiple European countries. The variant was first observed in Spain in June and has been at frequencies above 40% since July. Outside of Spain, the frequency of this variant has increased from very low values prior to 15th July to 40-70% in Switzerland, Ireland, and the United Kingdom in September. It is also prevalent in Norway, Latvia, the Netherlands, and France. Little can be said about other European countries because few recent sequences are available. Sequences in this cluster (20A.EU1) differ from ancestral sequences at 6 or more positions, including the mutation A222V in the spike protein and A220V in the nucleoprotein. We show that this variant was exported from Spain to other European countries multiple times and that much of the diversity of this cluster in Spain is observed across Europe. It is currently unclear whether this variant is spreading because of a transmission advantage of the virus or whether high incidence in Spain followed by dissemination through tourists is sufficient to explain the rapid rise in multiple countries., This work was supported by the SNF through grant numbers 31CA30 196046 (to RAN, EBH) and 31CA30 196267 and core funding by the University of Basel. SeqCOVID-SPAIN is funded by the Instituto de Salud Carlos III project COV20/00140, Spanish National Research Council and ERC StG 638553 to IC.

9 páginas, 4 figuras. All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors. Sequencing reads have been submitted to the European Nucleotide Archive (EMBL-EBI) under study accession PRJEB34582 and PRJEB34609. The analysis scripts used in this study are available online at GitHub (https://github.com/StopTB/Single_Colony_Project)., Mycobacterium tuberculosis (Mtb) infection results in a spectrum of clinical and histopathologic manifestations. It has been proposed that the environmental and immune pressures associated with different contexts of infection have different consequences for the associated bacterial populations, affecting drug susceptibility and the emergence of resistance. However, there is little concrete evidence for this model. We prospectively collected sputum samples from 18 newly diagnosed and treatment-naïve patients with tuberculosis and sequenced 795 colony-derived Mtb isolates. Mutant accumulation rates varied considerably between different bacilli isolated from the same individual, and where high rates of mutation were observed, the mutational spectrum was consistent with reactive oxygen species-induced mutagenesis. Elevated bacterial mutation rates were identified in isolates from HIV-negative but not HIV-positive individuals, suggesting that they were immune-driven. These results support the model that mutagenesis of Mtb in vivo is modulated by the host environment, which could drive the emergence of variants associated with drug resistance in a host-dependent manner., This work was supported by the National Natural Science Foundation of China (91631301 and 81661128043 to Q.G., 81701975 to Q.L., and 31771416 to X.L.), the National Science and Technology Major Project of China (2017ZX10201302 to Q.G. and 2018ZX10714002-001-005 to Z.Z.), the Sanming Project of Medicine in Shenzhen (SZSM201611030 to Q.G.), European Research Council 638553-TB-AcCELERATE (to I.C.), the Key Research Program of the Chinese Academy of Sciences (KFZD-SW-220-1 to X.L.), and the CAS Light of West China Program (to X.L.). Y.F. is supported in part by NIH R01HG009524. Support was also received from NIH awards P01 AI132130 and AI142793 to S.M.F.

11 páginas, 6 figuras, 2 tablas. Accesible en PubMedCentral: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6295914/, A small number of high-burden countries account for the majority of tuberculosis cases worldwide. Detailed data are lacking from these regions. To explore the evolutionary history of Mycobacterium tuberculosis in China-the country with the third highest tuberculosis burden-we analysed a countrywide collection of 4,578 isolates. Little genetic diversity was detected, with 99.4% of the bacterial population belonging to lineage 2 and three sublineages of lineage 4. The deeply rooted phylogenetic positions and geographic restriction of these four genotypes indicate that their populations expanded in situ following a small number of introductions to China. Coalescent analyses suggest that these bacterial subpopulations emerged in China around 1,000 years ago, and expanded in parallel from the twelfth century onwards, and that the whole population peaked in the late eighteenth century. More recently, sublineage L2.3, which is indigenous to China and exhibited relatively high transmissibility and extensive global dissemination, came to dominate the population dynamics of M. tuberculosis in China. Our results indicate that historical expansion of four M. tuberculosis strains shaped the current tuberculosis epidemic in China, and highlight the long-term genetic continuity of the indigenous M. tuberculosis population., This work was supported by the Natural Science Foundation of China (91631301 and 81661128043 to Q.G. and 81701975 to Q.L.). C.S.P. was supported by the National Institutes of Health (grant 1R01AI113287-01A1). S.G. was supported by the Swiss National Science Foundation (grants IZRJZ3_164171, 310030_166687, IZLSZ3_170834 and CRSII5_177163). This work was also supported by MINECO research grant SAF2016-77346-R (to I.C.), the European Research Council (638553-TB-ACCELERATE to I.C.), the National Science and Technology Major Project of China (2017ZX10201302 to Q.G., 2018ZX10103001 to Y.Z.), the Sanming Project of Medicine in Shenzhen (SZSM201611030 to Q.G.), JSGG20170413142559220 (to Q.G.), and National Basic Research programme of China (2014CB744403 to Y.Z.).

10 páginas, 3 figuras, 3 tablas. Material suplementario en: https://advances.sciencemag.org/content/suppl/2020/02/10/6.7.eaay4453.DC1, The postreplicative mismatch repair (MMR) is an almost ubiquitous DNA repair essential for maintaining genome stability. It has been suggested that Mycobacteria have an alternative MMR in which NucS, an endonuclease with no structural homology to the canonical MMR proteins (MutS/MutL), is the key factor. Here, we analyze the spontaneous mutations accumulated in a neutral manner over thousands of generations by Mycobacterium smegmatis and its MMR-deficient derivative (ΔnucS). The base pair substitution rates per genome per generation are 0.004 and 0.165 for wild type and ΔnucS, respectively. By comparing the activity of different bacterial MMR pathways, we demonstrate that both MutS/L- and NucS-based systems display similar specificity and mutagenesis bias, revealing a functional evolutionary convergence. However, NucS is not able to repair indels in vivo. Our results provide an unparalleled view of how this mycobacterial system works in vivo to maintain genome stability and how it may affect Mycobacterium evolution., This work was supported by Plan Nacional de I + D + i 2013–2016 and the Instituto de Salud Carlos III, Subdireccion General de Redes y Centros de Investigacion Cooperativa, Ministerio de Economia, Industria y Competitividad, Spanish Network for Research in Infectious Diseases (REIPI RD16/0016/0009)—cofinanced by the European Development Regional Fund “A Way to Achieve Europe”—by Operative Program Intelligent Growth 2014-2020 and grants FIS PI17/00159 from the Instituto de Salud Carlos III and SAF2015-72793-EXP from the Spanish Ministry of Science and Competitiveness (MINECO)-FEDER. E.C.-S. was the recipient of a PFIS predoctoral research fellowship (FI18/00036) cofinanced by the Instituto de Salud Carlos III and the European Social Fund. I.C., A.C.-O., and M.T.-P. are funded by projects of the 581 European Research Council (ERC) (638553-TB-ACCELERATE) and Ministerio de Economía y Competitividad research grant SAF2016-77346-R.

8 páginas, 4 figuras, 1 tabla. Información suplementaria accesible en: https://doi.org/10.1038/s41598-020-59373-w, It is relevant to evaluate MDR-tuberculosis in prisons and its impact on the global epidemiology of this disease. However, systematic molecular epidemiology programs in prisons are lacking. A health-screening program performed on arrival for inmates transferred from Peruvian prisons to Spain led to the diagnosis of five MDR-TB cases from one of the biggest prisons in Latin America. They grouped into two MIRU-VNTR-clusters (Callao-1 and Callao-2), suggesting a reservoir of two prevalent MDR strains. A high-rate of overexposure was deduced because one of the five cases was coinfected by a pansusceptible strain. Callao-1 strain was also identified in 2018 in a community case in Spain who had been in the same Peruvian prison in 2002-5. A strain-specific-PCR tailored from WGS data was implemented in Peru, allowing the confirmation that these strains were currently responsible for the majority of the MDR cases in that prison, including a new mixed infection., This work was supported by ERANet-LAC [ELAC2015/T08–0664, E035-ERANet-LAC/J110-2016/FONDECYT, PER-2012-ELAC2015/T08-0664] and Instituto de Salud Carlos III [AC16/00057, FIS15/01554, FIS13/01207] and cofunded by European Regional Development Funds from the European Commission: “A way of making Europe”. FA holds a grant from IFARHUSENACYT [270-2016-293].

14 páginas, 2 figuras, 2 tablas, The Dot/Icm type IVB secretion system of Legionella pneumophila is essential for its pathogenesis by delivering >300 effector proteins into the host cell. However, their precise secretion mechanism and which components interact with the host cell is only partly understood. Here, we undertook evolutionary analyses of the Dot/Icm system of 58 Legionella species to identify those components that interact with the host and/or the substrates. We show that high recombination rates are acting on DotA, DotG, and IcmX, supporting exposure of these proteins to the host. Specific amino acids under positive selection on the periplasmic region of DotF, and the cytoplasmic domain of DotM, support a role of these regions in substrate binding. Diversifying selection acting on the signal peptide of DotC suggests its interaction with the host after cleavage. Positive selection acts on IcmR, IcmQ, and DotL revealing that these components are probably participating in effector recognition and/or translocation. Furthermore, our results predict the participation in host/effector interaction of DotV and IcmF. In contrast, DotB, DotO, most of the core subcomplex elements, and the chaperones IcmS-W show a high degree of conservation and not signs of recombination or positive selection suggesting that these proteins are under strong structural constraints and have an important role in maintaining the architecture/function of the system. Thus, our analyses of recombination and positive selection acting on the Dot/Icm secretion system predicted specific Dot/Icm components and regions implicated in host interaction and/or substrate recognition and translocation, which will guide further functional analyses., This work was supported by the Institut Pasteur, the Agence National de la Recherche (grant nANR-10-LABX-62-IBEID) and the Fondation pour la Recherche Medicale (FRM) (grant NDEQ20120323697).

14 páginas, 6 figuras, Models on how bacterial lineages differentiate increase our understanding of early bacterial speciation events and the genetic loci involved. Here, we analyze the population genomics events leading to the emergence of the tuberculosis pathogen. The emergence is characterized by a combination of recombination events involving core pathogenesis functions and purifying selection on early diverging loci. We identify the phoR gene, the sensor kinase of a two-component system involved in virulence, as a key functional player subject to pervasive positive selection after the divergence of the Mycobacterium tuberculosis complex from its ancestor. Previous evidence showed that phoR mutations played a central role in the adaptation of the pathogen to different host species. Now, we show that phoR mutations have been under selection during the early spread of human tuberculosis, during later expansions, and in ongoing transmission events. Our results show that linking pathogen evolution across evolutionary and epidemiological time scales points to past and present virulence determinants., This work was funded by projects of the European Research Council (638553-TB-ACCELERATE) and Ministerio de Economía y Competitividad (Spanish government) research grant SAF2016-77346-R (to I.C.); BFU2014-58656-R and BFU2017-89594-R from Ministerio de Economía y Competitividad (Spanish government) and PROMETEO/2016/122 from Generalitat Valenciana (to F.G.-C); and the Swiss National Science Foundation (grants 310030_166687, IZRJZ3_164171, IZLSZ3_170834, and CRSII5_177163), the European Research Council (309540-EVODRTB), and SystemsX.ch (to S.G.). Á.C.-O. is the recipient of an FPU fellowship from Ministerio de Ciencia, Innovación y Universidades FPU13/00913 (Spanish government). L.S.-B. was funded by Wellcome grant number 098051.

5 páginas, 2 figuras. AVAILABILITY AND IMPLEMENTATION: The database of non-tuberculous mycobacteria assemblies can be accessed at: 10.5281/zenodo.3374377. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online: http://dx.doi.org/10.1093/bioinformatics/btz729, MOTIVATION: Tuberculosis (TB) remains one of the main causes of death worldwide. The long and cumbersome process of culturing Mycobacterium tuberculosis complex (MTBC) bacteria has encouraged the development of specific molecular tools for detecting the pathogen. Most of these tools aim to become novel TB diagnostics, and big efforts and resources are invested in their development, looking for the endorsement of the main public health agencies. Surprisingly, no study has been conducted where the vast amount of genomic data available is used to identify the best MTBC diagnostic markers. RESULTS: In this work, we used large-scale comparative genomics to identify 40 MTBC-specific loci. We assessed their genetic diversity and physiological features to select 30 that are good targets for diagnostic purposes. Some of these markers could be used to assess the physiological status of the bacilli. Remarkably, none of the most used MTBC markers is in our catalog. Illustrating the translational potential of our work, we develop a specific qPCR assay for quantification and identification of MTBC DNA. Our rational design of targeted molecular assays for TB could be used in many other fields of clinical and basic research., This work was supported by projects of the European Research Council (ERC) (638553-TB-ACCELERATE), Ministerio de Economía y Competitividad and Ministerio de Ciencia, Innovación y Universidades (Spanish Government), SAF2013-43521-R, SAF2016-77346-R and SAF2017–92345–EXP (to I.C.), BES-2014-071066 (to G.A.G.), FPU 13/00913 (to A.C.O.).

Review: 5 páginas, 1 figura, Two strains of Mycobacterium tuberculosis complex can be separated as much as 2500 single nucleotide differences (Coscolla and Gagneux, 2014). In that limited amount of diversity, we find an astonishing range of clinical, epidemiological and biological phenotypes. The most striking is the strong host preferences depending on the infecting strain while more subtle differences can be found looking at different human tuberculosis isolates. Those subtle differences are the most difficult to spot given that analysis methods for so little diversity are limited and phenotypes like virulence are difficult to define and measure. Recent genomics advances allow to study the pathogen diversity at a resolution not available before from comparative species level, to global diversity to transmission in local settings. Here, we will review some of these recent results to highlight how population genomics approaches can aid not only to understand how MTBC evolved but also to identify relevant biomedical targets., IC group is currently funded by projects of the European ResearchCouncil (ERC) research grant 638553-TB-ACCELERATE, Ministerio deEconomía y Competitividad (Spanish Government) research grant SAF2016-77346-R and Consellería de Educació, Investigació, Cultura iEsport (Valencian government) research grant AICO/2018/113. ACO isrecipient of a FPU fellowship from Ministerio de Ciencia, Innovación y Universidades FPU13/00913 (Spanish Government).

11 págs, 4 figuras y fórmulas matemáticas. Material suplementario en: http://dx.doi.org/10.1093/molbev/msaa069, The dN/dS ratio provides evidence of adaptation or functional constraint in protein-coding genes by quantifying the relative excess or deficit of amino acid-replacing versus silent nucleotide variation. Inexpensive sequencing promises a better understanding of parameters, such as dN/dS, but analyzing very large data sets poses a major statistical challenge. Here, I introduce genomegaMap for estimating within-species genome-wide variation in dN/dS, and I apply it to 3,979 genes across 10,209 tuberculosis genomes to characterize the selection pressures shaping this global pathogen. GenomegaMap is a phylogeny-free method that addresses two major problems with existing approaches: 1) It is fast no matter how large the sample size and 2) it is robust to recombination, which causes phylogenetic methods to report artefactual signals of adaptation. GenomegaMap uses population genetics theory to approximate the distribution of allele frequencies under general, parent-dependent mutation models. Coalescent simulations show that substitution parameters are well estimated even when genomegaMap's simplifying assumption of independence among sites is violated. I demonstrate the ability of genomegaMap to detect genuine signatures of selection at antimicrobial resistance-conferring substitutions in Mycobacterium tuberculosis and describe a novel signature of selection in the cold-shock DEAD-box protein A gene deaD/csdA. The genomegaMap approach helps accelerate the exploitation of big data for gaining new insights into evolution within species., D.J.W. is a Sir Henry Dale Fellow, jointly funded by the Wellcome Trust and the Royal Society (grant no. 101237/Z/13/B) and is a Big Data Institute Robertson Fellow. The CRyPTIC Consortium was supported by grants from the Bill and Melinda Gates Foundation (OPP1133541) and a Wellcome Trust/Newton Fund-MRC Collaborative Award (200205/Z/15/Z). F.A.D. was supported by the Imperial Biomedical Research Centre

10 Pages, 1 Figure, 3 Tables. Supplementary information: http://dx.doi.org/10.1038/s41598-018-33731-1, Drug-resistant tuberculosis poses a persistent public health threat. The ReSeqTB platform is a collaborative, curated knowledgebase, designed to standardize and aggregate global Mycobacterium tuberculosis complex (MTBC) variant data from whole genome sequencing (WGS) with phenotypic drug susceptibility testing (DST) and clinical data. We developed a unified analysis variant pipeline (UVP) ( https://github.com/CPTR-ReSeqTB/UVP ) to identify variants and assign lineage from MTBC sequence data. Stringent thresholds and quality control measures were incorporated in this open source tool. The pipeline was validated using a well-characterized dataset of 90 diverse MTBC isolates with conventional DST and DNA Sanger sequencing data. The UVP exhibited 98.9% agreement with the variants identified using Sanger sequencing and was 100% concordant with conventional methods of assigning lineage. We analyzed 4636 publicly available MTBC isolates in the ReSeqTB platform representing all seven major MTBC lineages. The variants detected have an above 94% accuracy of predicting drug based on the accompanying DST results in the platform. The aggregation of variants over time in the platform will establish confidence-graded mutations statistically associated with phenotypic drug resistance. These tools serve as critical reference standards for future molecular diagnostic assay developers, researchers, public health agencies and clinicians working towards the control of drug-resistant tuberculosis., This study was supported by the Bill & Melinda Gates Foundation under grant agreement OPP1115887 to C-Path for developing the ReSeqTB drug resistance data sharing platform and under grant agreement FIND OPP1115209 to address how to score mutations in the ReSeqTB data sharing platform initiative. The South African MRC and the EDCTP support K. Dheda. I. Comas is supported by the Ministerio de Economía y Competitividad (Spanish Government) research grant SAF2016-77346-R and the European Research Council (ERC) (638553-TB-ACCELERATE). L. Chindelevitch acknowledges support by NSERC, Genome Canada, and the Sloan Foundation. Use of trade names is for identification only and does not constitute endorsement by the US Department of Health and Human Services, the US Public Health Service, or the Centers for Disease Control and Prevention. The findings and conclusions in this report are those of the authors and do not necessarily represent the views of the funding agency.

18 páginas, 11 figuras, 3 tablas. Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2019.02464/full#supplementary-material, The administration of a high fat content diet is an accelerating factor for metabolic syndrome, impaired glucose tolerance, and early type 2 diabetes. The present study aims to assess the impact of a high fat diet on tuberculosis progression and microbiota composition in an experimental animal model using a C3HeB/FeJ mouse strain submitted to single or multiple consecutive aerosol infections. These models allowed us to study the protection induced by Bacillus Calmette-Guérin vaccination as well as by the natural immunity induced by chemotherapy after a low dose Mycobacterium tuberculosis infection. Our results show that a high fat diet is able to trigger a pro-inflammatory response, which results in a faster progression toward active tuberculosis and an impaired protective effect of BCG vaccination, which is not the case for natural immunity. This may be related to dysbiosis and a reduction in the Firmicutes/Bacteroidetes ratio in the gut microbiota caused by a decrease in the abundance of the Porphyromonadaceae family and, in particular, the Barnesiella genus. It should also be noted that a high fat diet is also related to an increase in the genera Alistipes, Parasuterella, Mucispirillum, and Akkermansia, which have previously been related to dysbiotic processes. As diabetes mellitus type 2 is a risk factor for developing tuberculosis, these findings may prove useful in the search for new prophylactic strategies for this population subset., This work was supported by the European Commission Horizon 2020 research and innovation program under grant agreement TBVAC2020 No. 643381, European Research Council research grant 638553-TB-ACCELERATE and Ministerio de Economía y Competitividad (Spanish government) research grant SAF2017-92345- EXP (to IC). CV has a project Servet II contract funded by the Instituto de Salud Carlos III (ISCIII, CPII 18/00031). LA received an Improvement and Mobility Program grant from the Centro de Investigación Biomédica en Red de Enfermedades Respiratorias (CIBERES). GG contract was funded by Ministerio de Economia y Competitividad (Spanish Government) (BES-2014-071066). PC contract was funded by the Instituto de Salud Carlos III (IFI14/00015). The Experimental Tuberculosis Unit is accredited by the Catalan Agency for Management of University and Research Grants with code 2017 SGR500 and the IGTP is a member of the CERCA network of institutes. The funders played no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

12 páginas, 6 figuras, 4 tablas, The Mycobacterium tuberculosis complex (MTBC) comprises the species that causes tuberculosis (TB) which affects 10 million people every year. A robust classification of species, lineages, and sub-lineages is important to explore associations with drug resistance, epidemiological patterns or clinical outcomes. We present a rapid and easy-to-follow methodology to classify clinical TB samples into the main MTBC clades. Approaches are based on the identification of lineage and sub-lineage diagnostic SNP using a real-time PCR high resolution melting assay and classic Sanger sequencing from low-concentrated, low quality DNA. Thus, suitable for implementation in middle and low-income countries. Once we validated our molecular procedures, we characterized a total of 491 biological samples from human and cattle hosts, representing countries with different TB burden. Overall, we managed to genotype ~95% of all samples despite coming from unpurified and low-concentrated DNA. Our approach also allowed us to detect zoonotic cases in eight human samples from Nigeria. To conclude, the molecular techniques we have developed, are accurate, discriminative and reproducible. Furthermore, it costs less than other classic typing methods, resulting in an affordable alternative method in TB laboratories., This work was supported by European Research Council (grant, 638553-TB-ACCELERATE) and the Ministerio de Economía y Competitividad grants (SAF2013-43521-R and SAF2016-77346-R) both to I.C.

5 páginas, 2 figuras, Understanding Multidrug-resistant tuberculosis cases not detected by rapid phenotypic and genotypic routine clinical tests is essential to improve diagnostics and advance towards personalized TB treatments. Here, we combine WGS with single-colony phenotyping testing to diagnose a MDR strain infecting a patient over nine years. Our investigation revealed the failure of rapid testing assays and genome-based prediction tools to identify the MDR strain. The false negative outcome was caused by uncommon RIF and INH resistance mutations. Despite WGS data helped to personalize treatment, the patient developed XDR-TB, highlighting the importance of coupling new diagnostic methods with appropriate treatment regimens., This work was supported by the European Research Council (638553-TB-ACCELERATE to I. C.) and the Ministerio de Economía y Competitividad (SAF2016-77346-R to I. C.).

Artículo con 20 páginas, 5 figuras, 1 tabla. All the sequence data are deposited in the European Nucleotide Archive under the Bioproject number PRJEB29604 (https://www.ebi.ac.uk/ena/data/view/PRJEB29604) and the accession numbers ERR2099780 (https://www.ebi.ac.uk/ena/data/search?query=ERR2099780) and ERR2099784 (https://www.ebi.ac.uk/ena/data/search?query=ERR2099784)., 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., IC was supported by European Research Council (638553-TB-ACCELERATE), the Ministerio de Economía y Competitividad (SAF2016-77346-R). CC and YX were supported by the Engineering and Physical Sciences Research Council of the UK (EPSRC EP/K026003/1 (CC) and EPSRC EP/N014529/1 (CC and YX).

14 páginas, 7 figuras, BackgroundThe analysis of transmission of tuberculosis (TB) is challenging in areas with a large migrant population. Standard genotyping may fail to differentiate transmission within the host country from new importations, which is key from an epidemiological perspective.AimTo propose a new strategy to simplify and optimise cross-border surveillance of tuberculosis and to distinguish between recent transmission in the host country and new importationsMethodsWe selected 10 clusters, defined by 24-locus mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR), from a population in Spain rich in migrants from eastern Europe, north Africa and west Africa and reanalysed 66 isolates by whole-genome sequencing (WGS). A multiplex-allele-specific PCR was designed to target strain-specific marker single nucleotide polymorphisms (SNPs), identified from WGS data, to optimise the surveillance of the most complex cluster.ResultsIn five of 10 clusters not all isolates showed the short genetic distances expected for recent transmission and revealed a higher number of SNPs, thus suggesting independent importations of prevalent strains in the country of origin. In the most complex cluster, rich in Moroccan cases, a multiplex allele-specific oligonucleotide-PCR (ASO-PCR) targeting the marker SNPs for the transmission subcluster enabled us to prospectively identify new secondary cases. The ASO-PCR-based strategy was transferred and applied in Morocco, demonstrating that the strain was prevalent in the country.ConclusionWe provide a new model for optimising the analysis of cross-border surveillance of TB transmission in the scenario of global migration., This project was funded by ISCIII: ERANET-LAC (TRANS-TB-TRANS REF AC16/00057; ELAC2015/T08-0664), FIS (13/01207; 15/01554) and cofunded by ERDF Funds from the European Commission: “A way of making Europe”. Miguel Servet grant (CP15/00075) for LPL. Ministerio de Economía y Competitividad (grant SAF2016-77346-R), ERC (638553-TB-ACCELERATE) to IC. FPU13/00913 (Ministerio de Educación y Ciencia) to ACO.

2 Páginas

17 páginas, 1 figura, 7 tablas, BACKGROUND: Whole genome sequencing (WGS) has been proposed as a tool for diagnosing drug resistance in tuberculosis. However, reports of its effectiveness in endemic countries with important numbers of drug resistance are scarce. The goal of this study was to evaluate the effectiveness of this procedure in isolates from a tuberculosis endemic region in Mexico. METHODS: WGS analysis was performed in 81 tuberculosis positive clinical isolates with a known phenotypic profile of resistance against first-line drugs (isoniazid, rifampin, ethambutol, pyrazinamide and streptomycin). Mutations related to drug resistance were identified for each isolate; drug resistant genotypes were predicted and compared with the phenotypic profile. Genotypes and transmission clusters based on genetic distances were also characterized. FINDINGS: Prediction by WGS analysis of resistance against isoniazid, rifampicin, ethambutol, pyrazinamide and streptomycin showed sensitivity values of 84%, 96%, 71%, 75% and 29%, while specificity values were 100%, 94%, 90%, 90% and 98%, respectively. Prediction of multidrug resistance showed a sensitivity of 89% and specificity of 97%. Moreover, WGS analysis revealed polymorphisms related to second-line drug resistance, enabling classification of eight and two clinical isolates as pre- and extreme drug-resistant cases, respectively. Lastly, four lineages were identified in the population (L1, L2, L3 and L4). The most frequent of these was L4, which included 90% (77) of the isolates. Six transmission clusters were identified; the most frequent was TC6, which included 13 isolates with a L4.1.1 and a predominantly multidrug-resistant condition. CONCLUSIONS: The results illustrate the utility of WGS for establishing the potential for prediction of resistance against first and second line drugs in isolates of tuberculosis from the region. They also demonstrate the feasibility of this procedure for use as a tool to support the epidemiological surveillance of drug- and multidrug-resistant tuberculosis., CFM-M was a CONACYT fellow of the Maestría en Ciencias de la Salud of the Instituto de Ciencias de la Salud Universidad Veracruzana. RZ-C and IC were partially funded by I-COOP-2017- COOPB2032. RMS was funded by CONACyT-Problemas nacionales grant no. 2015-01-147. XS was funded by GACD-FONCICYT Conacyt No.264693.

8 páginas, 3 figuras, We show the distribution of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) genetic clades over time and between countries and outline potential genomic surveillance objectives. We applied three genomic nomenclature systems to all sequence data from the World Health Organization European Region available until 10 July 2020. We highlight the importance of real-time sequencing and data dissemination in a pandemic situation, compare the nomenclatures and lay a foundation for future European genomic surveillance of SARS-CoV-2., We gratefully acknowledge the authors, originating and submitting laboratories of the sequences from GISAID’s EpiCoV Database used in the phylogenetic analysis. We gratefully acknowledge all the staff working with sample collection, sample preparation, sequencing, data analysis and data sharing in all laboratories in the WHO European Region for making this work possible.

15 páginas, 6 figuras, 7 tablas., 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., This project has been funded by the European Research Council (638553-TB-ACCELERATE, 101001038-TB-RECONNECT), the Ministerio de Economía, Indústria y Competividad (SAF2016-77346-R, PID2019-104477RB-I00) and the Conselleria d'Educació, Investigació, Cultura i Esport (PROMETEO 2020/012).

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 now 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., This research work was supported by the European Commission–NextGenerationEU through the CSIC Global Health Platform. Additionally, authors would like to acknowledge economic support from the Spanish Ministry of Science and Innovation through Grants: PID2019-104757RB-I00 funded by MCIN/AEI/ 10.13039/501100011033, RTI2018-094399-A-I00, and “ERDF A way of making Europe”, by the “European Union”, Grant SEV 2017-0712 funded by MCIN/AEI /10.13039/501100011033, the “Comunidad Autónoma de Madrid" through Grant: S2017/BMD3817, and the European Union (EU) and Horizon 2020 through grants: Marie-Curie Fellowship EnLaCES (MSCA IF 2020, Proposal: 101024130) (to JK), HighResCells (ERC - 2018 - SyG, Proposal: 810057), and iNEXT-Discovery (Proposal: 871037). AM, VR, JB and JLL are funded by CIBERER-ISCIII (proposal: COV20/00437), Fondo Supera COVID-19 (proposal: CSICCOVID19-082), Banco Santander (Proposal: BlockAce), and CSIC PTI Salud Global (Proposal: 202080E110). VR is funded by the Spanish Ministry of Science and Innovation through Grant PID2020-120322RB-C21. IC is funded by project PID2019-104477RB-100, Fondo COVID COV20/00140 and ERC CoG 101001038. MC is funded by the RyC program from the Spanish Ministry of Science and Innovation, the Generalitat Valenciana (SEJI/2019/011).

12 páginas, 5 figuras, 2 tablas. DNA sequences generated as part of this study have been deposited at the European Nucleotide Archive under study accession number PRJEB32589 (Table S1) (2019), Tuberculosis (TB) surveillance is scarce in most African countries, even though it is the continent with the greatest disease incidence according to the World Health Organization. Liberia is within the 30 countries with the highest TB burden, probably as a consequence of the long civil war and the recent Ebola outbreak, both crippling the health system and depreciating the TB prevention and control programmes. Due to difficulties working in the country, there is a lack of resistance surveys and bacillus characterization. Here, we use genome sequencing of Mycobacterium tuberculosis clinical isolates to fill this gap. Our results highlight that the bacillus population structure is dominated by lineage 4 strains that harbour an outstanding genetic diversity, higher than in the rest of Africa as a whole. Coalescent analyses demonstrate that strains currently circulating in Liberia were introduced several times beginning in the early year 600 CE until very recently coinciding with migratory movements associated with the civil war and Ebola epidemics. A higher multidrug-resistant (MDR)-TB frequency (23.5 %) than current estimates was obtained together with non-catalogued drug-resistance mutations. Additionally, 39 % of strains were in genomic clusters revealing that ongoing transmission is a major contribution to the TB burden in the country. Our report emphasizes the importance of TB surveillance and control in African countries where bacillus diversity, MDR-TB prevalence and transmission are coalescing to jeopardize TB control programmes., This work was supported by the European Research Council (638553-TB-ACCELERATE) and the Ministerio de Economía y Competitividad (SAF2013-43521-R and SAF2016-77346-R to I.C.). Funding support from the John D. and Catherine T. MacArthur Foundation, USA, under the Higher Education Initiative in Africa (grant no. 97944-0-800/406/99) for the establishment of the Center for Control and Prevention of Zoonoses, University of Ibadan, Nigeria, is appreciated.