9 páginas, 3 figuras ; Background Direct whole-genome sequencing of Mycobacterium tuberculosis from clinical specimens will be a major breakthrough in tuberculosis diagnosis and control. To date, direct whole-genome sequencing has never been used in genomic epidemiology, and its accuracy in transmission inference remains unknown. We investigated the technical challenges imposed by direct whole-genome sequencing, and used it to infer transmission clusters and predict drug resistance. Methods Using an optimised workflow, we did direct whole-genome sequencing for 37 clinical specimens from 23 tuberculosis patients. Nine sputum samples from nine patients who were infected with different non-tuberculous mycobacteria and culture-negative for tuberculosis were used as controls in the qPCR assays and pre-sequencing runs. Additionally, 780 clinical isolates in the region of Comunidad Valenciana (Spain) were whole-genome sequenced between Jan 1, 2014, and Dec 31, 2016. We analysed the genomic variants to build a tuberculosis transmission network for the region, including the clinical specimens, and to predict drug susceptibility profiles. Findings After sequencing 37 clinical specimens, 28 specimens (22 [85%] of 26 smear-positive and six [55%] of 11 smear-negative) met the quality criteria for downstream analysis. All 28 clinical specimens clustered with their matching culture isolates, with a median distance of 0 single nucleotide polymorphisms. Of the 28 clinical specimens, 16 (57%) were accurately assigned to ten transmission clusters in the region, and 12 (43%) were unique cases. Transmission inferences and drug-susceptibility predictions from direct whole-genome sequencing data were concordant with sequences from corresponding cultures and phenotypic drug-susceptibility testing. Complete genomic analysis, within a week of specimen receipt, cost €217 per sample (excluding personnel costs). Interpretation Direct whole-genome sequencing could be used to accurately delineate transmission clusters of tuberculosis and conduct culture-independent surveillance. Compared with conventional approaches, direct wholegenome sequencing allows researchers to do real-time genomic epidemiology and drug resistance surveillance in settings where culture and drug susceptibility testing are not available. ; This work was supported by projects of the European Research Council (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 IC], and BES-2014-071066 [to GAG]) ; Peer reviewed
Ruiz-Rodriguez, P.; Frances-Gomez, C.; Chiner-Oms, A.; Lopez, M. G.; Jimenez-Serrano, S.; Cancino-Munoz, I.; Ruiz-Hueso, P.; Torres-Puente, M.; Bracho, M. A.; D'Auria, G.; Martinez-Priego, L.; Guerreiro, M.; Montero-Alonso, M.; Gomez, M. D.; Pinana, J. L.; Seq, Covid-Spain Consortium; Gonzalez-Candelas, F.; Comas, I.; Marina, A.; Geller, R.; Coscolla, M. Evolutionary and Phenotypic Characterization of Two Spike Mutations in European Lineage 20E of SARS-CoV-2 [Preprint] MedRxiv; https://dx.doi.org/10.1101/2021.03.08.21253075 ; We have detected two mutations in the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) at amino acid positions 1163 and 1167 that appeared independently in multiple transmission clusters and different genetic backgrounds. Furthermore, both mutations appeared together in a cluster of 1,627 sequences belonging to clade 20E. This cluster is characterized by 12 additional single nucleotide polymorphisms but no deletions. The available structural information on the S protein in the pre- and postfusion conformations predicts that both mutations confer rigidity, which could potentially decrease viral fitness. Accordingly, we observed reduced infectivity of this spike genotype relative to the ancestral 20E sequence in vitro, and the levels of viral RNA in nasopharyngeal swabs were not significantly higher. Furthermore, the mutations did not impact thermal stability or antibody neutralization by sera from vaccinated individuals but moderately reduce neutralization by convalescent-phase sera from the early stages of the pandemic. Despite multiple successful appearances of the two spike mutations during the first year of SARS-CoV-2 evolution, the genotype with both mutations was displaced upon the expansion of the 20I (Alpha) variant. The midterm fate of the genotype investigated was consistent with the lack of advantage observed in the clinical and experimental data. IMPORTANCE We observed repeated, independent emergence of mutations in the SARS-CoV-2 spike involving amino acids 1163 and 1167, within the HR2 functional motif. Conclusions derived from evolutionary and genomic diversity analysis suggest that the co-occurrence of both mutations might pose an advantage for the virus and therefore a threat to effective control of the epidemic. However, biological characterization, including in vitro experiments and analysis of clinical data, indicated no clear benefit in terms of stability or infectivity. In agreement with this, continuous epidemiological surveillance conducted months after the first observations revealed that both mutations did not successfully outcompete other variants and stopped circulating 9 months after their initial detection. Additionally, we evaluated the potential of both mutations to escape neutralizing antibodies, finding that the presence of these two mutations on their own is not likely to confer antibody escape. Our results provide an example of how newly emerged spike mutations can be assessed to better understand the risk posed by new variants and indicate that some spike mutations confer no clear advantage to the virus despite independently emerging multiple times and are eventually displaced by fitter variants. ; This work was funded by the Instituto de Salud Carlos III project COV20/00140 and COV20/00437, Spanish National Research Council project CSIC-COV19-021 and CSIC-COVID19-082, and the Generalitat Valenciana (SEJI/2019/011 and Covid_19-SCI). Action co-financed by the European Union through the Operational Program of the European Regional Development Fund (ERDF) of the Valencian Community 2014-2020. M.C. and R.G. are supported by Ramon y Cajal program from Ministerio de Ciencia. ; Peer reviewed
Molecular epidemiology of SARS-CoV-2 aims to monitor the appearance of new variants with the potential to change the virulence or transmissibility of the virus. During the first year of SARS-CoV-2 evolution, numerous variants with possible public health impact have emerged. We have detected two mutations in the Spike protein at amino acid positions 1163 and 1167 that have appeared independently multiple times in different genetic backgrounds, indicating they may increase viral fitness. Interestingly, the majority of these sequences appear in transmission clusters, with the genotype encoding mutations at both positions increasing in frequency more than single-site mutants. This genetic outcome that we denote as Lineage B.1.177.637, belongs to clade 20E and includes 12 additional single nucleotide polymorphisms but no deletions with respect to the reference genome (first sequence in Wuhan). B.1.177.637 appeared after the first wave of the epidemic in Spain, and subsequently spread to eight additional countries, increasing in frequency among sequences in public databases. Positions 1163 and 1167 of the Spike protein are situated in the HR2 domain, which is implicated in the fusion of the host and viral membranes. To better understand the effect of these mutations on the virus, we examined whether B.1.177.637 altered infectivity, thermal stability, or antibody sensitivity. Unexpectedly, we observed reduced infectivity of this variant relative to the ancestral 20E variant in vitro while the levels of viral RNA in nasopharyngeal swabs did not vary significantly. In addition, we found the mutations do not impact thermal stability or antibody susceptibility in vaccinated individuals but display a moderate reduction in sensitivity to neutralization by convalescent sera from early stages of the pandemic. Altogether, this lineage could be considered a Variant of Interest (VOI), we denote VOI1163.7. Finally, we detected a sub-cluster of sequences within VOI1163.7 that have acquired two additional changes previously associated with antibody escape and it could be identified as VOI1163.7.V2. Overall, we have detected the spread of a new Spike variant that may be advantageous to the virus and whose continuous transmission poses risks by the acquisition of additional mutations that could affect pre-existing immunity. ; This work was funded by the Instituto de Salud Carlos III project COV20/00140 and COV20/00437, Spanish National Research Council project CSIC-COV19-021 and CSIC-COVID19-082, and the Generalitat Valenciana (SEJI/2019/011 and Covid_19-SCI). Action co-financed by the European Union through the Operational Program of the European Regional Development Fund (ERDF) of the Valencian Community 2014-2020. M.C. and R.G. are supported by Ramon y Cajal program from Ministerio de Ciencia. ; No
13 Pages, 1 Figure, 4 tables. The authors' affiliations are listed in the Supplementary Appendix, available at NEJM.org. Supplementary Material, available at http://dx.doi.org/10.1056/NEJMoa1800474 ; BACKGROUND: The World Health Organization recommends drug-susceptibility testing of Mycobacterium tuberculosis complex for all patients with tuberculosis to guide treatment decisions and improve outcomes. Whether DNA sequencing can be used to accurately predict profiles of susceptibility to first-line antituberculosis drugs has not been clear. METHODS: We obtained whole-genome sequences and associated phenotypes of resistance or susceptibility to the first-line antituberculosis drugs isoniazid, rifampin, ethambutol, and pyrazinamide for isolates from 16 countries across six continents. For each isolate, mutations associated with drug resistance and drug susceptibility were identified across nine genes, and individual phenotypes were predicted unless mutations of unknown association were also present. To identify how whole-genome sequencing might direct first-line drug therapy, complete susceptibility profiles were predicted. These profiles were predicted to be susceptible to all four drugs (i.e., pansusceptible) if they were predicted to be susceptible to isoniazid and to the other drugs or if they contained mutations of unknown association in genes that affect susceptibility to the other drugs. We simulated the way in which the negative predictive value changed with the prevalence of drug resistance. RESULTS: A total of 10,209 isolates were analyzed. The largest proportion of phenotypes was predicted for rifampin (9660 [95.4%] of 10,130) and the smallest was predicted for ethambutol (8794 [89.8%] of 9794). Resistance to isoniazid, rifampin, ethambutol, and pyrazinamide was correctly predicted with 97.1%, 97.5%, 94.6%, and 91.3% sensitivity, respectively, and susceptibility to these drugs was correctly predicted with 99.0%, 98.8%, 93.6%, and 96.8% specificity. Of the 7516 isolates with complete phenotypic drug-susceptibility profiles, 5865 (78.0%) had complete genotypic predictions, among which 5250 profiles (89.5%) were correctly predicted. Among the 4037 phenotypic profiles that were predicted to be pansusceptible, 3952 (97.9%) were correctly predicted. CONCLUSIONS: Genotypic predictions of the susceptibility of M. tuberculosis to first-line drugs were found to be correlated with phenotypic susceptibility to these drugs. (Funded by the Bill and Melinda Gates Foundation and others.). ; Supported by grants from the Bill and Melinda Gates Foundation (OPP1133541, to CRyPTIC, plus separate support to Dr. Rodwell), a Wellcome Trust/Newton Fund–MRC Collaborative Award (200205/Z/15/Z, to CRyPTIC), the National Institute for Health Research (NIHR) Oxford Biomedical Research Centre (BRC) and NIHR Health Protection Research Unit in Healthcare Associated Infections and Antimicrobial Resistance at University of Oxford in partnership with Public Health England, the NIHR Biomedical Research Centre at Barts, the NIHR Biomedical Research Centre at Imperial, the NIHR and NHS England (to the 100,000 Genomes Project, which is managed by Genomics England, a wholly owned company of the U.K. Department of Health), the Wellcome Trust, the Medical Research Council, Public Health England, a grant from the National Science and Technology Key Program of China (2014ZX10003002), a grant from the National Basic Research program of China (2014CB744403), a grant from the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB29020000), a grant from the European Commission Seventh Framework Program (FP7/2007-2013, to Borstel under grant agreement 278864 in the framework of the Patho-NGen-Trace project), the German Center for Infection Research (to Borstel), Leibniz Science Campus Evolutionary Medicine of the Lung (EvoLUNG), the Belgian Ministry of Social Affairs (to the Belgian Reference Center for Tuberculosis and Mycobacteria from Bacterial Diseases Service through a fund within the Health Insurance System), the French governmental program "Investing for the Future" (to Genoscreen), a grant from the European Commission Seventh Framework Program (FP7/2007-2013, to Genoscreen under grant agreement 278864 in the framework of the Patho-NGen-Trace project), grants from the Drug Resistant Tuberculosis Fund (R015833003, to Dr. Chaiprasert), the Faculty of Medicine, Siriraj Hospital, Mahidol University (to Dr. Chaiprasert), a grant from the Ministry of Economy and Competitiveness (MINECO), Spain (SAF2016-77346-R, to Dr. Comas), a grant from the European Research Council (638553-TB-ACCELERATE, to Dr. Comas), a grant from the BC Centre for Disease Control Foundation for Population and Public Health (to Dr. Gardy), a grant from the British Colombia Lung Association (to Dr. Gardy), grants from the Wellcome Trust and the Royal Society (101237/Z/13/Z and 102541/A/13/Z, to Drs. Wilson and Iqbal [Sir Henry Dale Fellows]), a grant from the National University of Singapore Yong Loo Lin School of Medicine Aspiration Fund (NUHSRO/2014/069/AF-New Idea/04, to Drs. Ong and Teo), a European Commission Seventh Framework Program European Genetic Network (EUROGEN) grant (201483, to Dr. Drobniewski), and the National Institute of Allergy and Infectious Diseases, National Institutes of Health (to Dr. Rodwell). Dr. T. Walker is an NIHR Academic Clinical Lecturer, and Drs. Crook, Peto, and Caulfield are NIHR Senior Investigators. No potential conflict of interest relevant to this article was reported. Disclosure forms provided by the authors are available with the full text of this article at NEJM.org. We thank Stéphanie Duthoy, Carina Hahn, Alamdar Hussain, Yannick Laurent, Mathilde Mairey, Vanessa Mohr, and Mahmood Qadir for technical assistance and George F. Gao, Director of the Chinese Center for Disease Control and Prevention, for directing the Chinese grant and sequencing program ; Peer reviewed