1
|
Viana R, Moyo S, Amoako DG, Tegally H, Scheepers C, Althaus CL, Anyaneji UJ, Bester PA, Boni MF, Chand M, Choga WT, Colquhoun R, Davids M, Deforche K, Doolabh D, du Plessis L, Engelbrecht S, Everatt J, Giandhari J, Giovanetti M, Hardie D, Hill V, Hsiao NY, Iranzadeh A, Ismail A, Joseph C, Joseph R, Koopile L, Kosakovsky Pond SL, Kraemer MUG, Kuate-Lere L, Laguda-Akingba O, Lesetedi-Mafoko O, Lessells RJ, Lockman S, Lucaci AG, Maharaj A, Mahlangu B, Maponga T, Mahlakwane K, Makatini Z, Marais G, Maruapula D, Masupu K, Matshaba M, Mayaphi S, Mbhele N, Mbulawa MB, Mendes A, Mlisana K, Mnguni A, Mohale T, Moir M, Moruisi K, Mosepele M, Motsatsi G, Motswaledi MS, Mphoyakgosi T, Msomi N, Mwangi PN, Naidoo Y, Ntuli N, Nyaga M, Olubayo L, Pillay S, Radibe B, Ramphal Y, Ramphal U, San JE, Scott L, Shapiro R, Singh L, Smith-Lawrence P, Stevens W, Strydom A, Subramoney K, Tebeila N, Tshiabuila D, Tsui J, van Wyk S, Weaver S, Wibmer CK, Wilkinson E, Wolter N, Zarebski AE, Zuze B, Goedhals D, Preiser W, Treurnicht F, Venter M, Williamson C, Pybus OG, Bhiman J, Glass A, Martin DP, Rambaut A, Gaseitsiwe S, von Gottberg A, de Oliveira T. Rapid epidemic expansion of the SARS-CoV-2 Omicron variant in southern Africa. Nature 2022; 603:679-686. [PMID: 35042229 PMCID: PMC8942855 DOI: 10.1038/s41586-022-04411-y] [Citation(s) in RCA: 918] [Impact Index Per Article: 459.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2021] [Accepted: 01/07/2022] [Indexed: 01/02/2023]
Abstract
The SARS-CoV-2 epidemic in southern Africa has been characterized by three distinct waves. The first was associated with a mix of SARS-CoV-2 lineages, while the second and third waves were driven by the Beta (B.1.351) and Delta (B.1.617.2) variants, respectively1-3. In November 2021, genomic surveillance teams in South Africa and Botswana detected a new SARS-CoV-2 variant associated with a rapid resurgence of infections in Gauteng province, South Africa. Within three days of the first genome being uploaded, it was designated a variant of concern (Omicron, B.1.1.529) by the World Health Organization and, within three weeks, had been identified in 87 countries. The Omicron variant is exceptional for carrying over 30 mutations in the spike glycoprotein, which are predicted to influence antibody neutralization and spike function4. Here we describe the genomic profile and early transmission dynamics of Omicron, highlighting the rapid spread in regions with high levels of population immunity.
Collapse
Affiliation(s)
| | - Sikhulile Moyo
- Botswana Harvard AIDS Institute Partnership, Botswana Harvard HIV Reference Laboratory, Gaborone, Botswana
- Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Botswana Presidential COVID-19 Taskforce, Gaborone, Botswana
| | - Daniel G Amoako
- National Institute for Communicable Diseases (NICD) of the National Health Laboratory Service (NHLS), Johannesburg, South Africa
| | - Houriiyah Tegally
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), Nelson R. Mandela School of Medicine, University of KwaZulu-Natal, Durban, South Africa
| | - Cathrine Scheepers
- National Institute for Communicable Diseases (NICD) of the National Health Laboratory Service (NHLS), Johannesburg, South Africa
- South African Medical Research Council Antibody Immunity Research Unit, School of Pathology, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
| | - Christian L Althaus
- Institute of Social and Preventive Medicine, University of Bern, Bern, Switzerland
| | - Ugochukwu J Anyaneji
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), Nelson R. Mandela School of Medicine, University of KwaZulu-Natal, Durban, South Africa
| | - Phillip A Bester
- Division of Virology, National Health Laboratory Service, Bloemfontein, South Africa
- Division of Virology, University of the Free State, Bloemfontein, South Africa
| | - Maciej F Boni
- Center for Infectious Disease Dynamics, Department of Biology, Pennsylvania State University, University Park, PA, USA
| | | | | | - Rachel Colquhoun
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK
| | - Michaela Davids
- Zoonotic Arbo and Respiratory Virus Program, Centre for Viral Zoonoses, Department of Medical Virology, University of Pretoria, Pretoria, South Africa
| | | | - Deelan Doolabh
- Division of Medical Virology, Faculty of Health Sciences, University of Cape Town, Cape Town, South Africa
| | - Louis du Plessis
- Department of Zoology, University of Oxford, Oxford, UK
- Department of Biosystems Science and Engineering, ETH Zurich, Zurich, Switzerland
| | - Susan Engelbrecht
- Division of Medical Virology, Faculty of Medicine and Health Sciences, Stellenbosch University, Tygerberg, Cape Town, South Africa
| | - Josie Everatt
- National Institute for Communicable Diseases (NICD) of the National Health Laboratory Service (NHLS), Johannesburg, South Africa
| | - Jennifer Giandhari
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), Nelson R. Mandela School of Medicine, University of KwaZulu-Natal, Durban, South Africa
| | - Marta Giovanetti
- Laboratorio de Flavivirus, Fundacao Oswaldo Cruz, Rio de Janeiro, Brazil
- Laboratório de Genética Celular e Molecular, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil
| | - Diana Hardie
- Division of Medical Virology, Faculty of Health Sciences, University of Cape Town, Cape Town, South Africa
- Division of Virology, NHLS Groote Schuur Laboratory, Cape Town, South Africa
| | - Verity Hill
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK
| | - Nei-Yuan Hsiao
- Division of Medical Virology, Faculty of Health Sciences, University of Cape Town, Cape Town, South Africa
- Division of Virology, NHLS Groote Schuur Laboratory, Cape Town, South Africa
- Wellcome Centre for Infectious Diseases Research in Africa (CIDRI-Africa), Cape Town, South Africa
| | - Arash Iranzadeh
- Division of Computational Biology, Faculty of Health Sciences, University of Cape Town, Cape Town, South Africa
| | - Arshad Ismail
- National Institute for Communicable Diseases (NICD) of the National Health Laboratory Service (NHLS), Johannesburg, South Africa
| | | | - Rageema Joseph
- Division of Medical Virology, Faculty of Health Sciences, University of Cape Town, Cape Town, South Africa
| | - Legodile Koopile
- Botswana Harvard AIDS Institute Partnership, Botswana Harvard HIV Reference Laboratory, Gaborone, Botswana
| | - Sergei L Kosakovsky Pond
- Institute for Genomics and Evolutionary Medicine, Department of Biology, Temple University, Philadelphia, PA, USA
| | | | - Lesego Kuate-Lere
- Health Services Management, Ministry of Health and Wellness, Gaborone, Botswana
| | - Oluwakemi Laguda-Akingba
- NHLS Port Elizabeth Laboratory, Port Elizabeth, South Africa
- Faculty of Health Sciences, Walter Sisulu University, Mthatha, South Africa
| | - Onalethatha Lesetedi-Mafoko
- Public Health Department, Integrated Disease Surveillance and Response, Ministry of Health and Wellness, Gaborone, Botswana
| | - Richard J Lessells
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), Nelson R. Mandela School of Medicine, University of KwaZulu-Natal, Durban, South Africa
| | - Shahin Lockman
- Botswana Harvard AIDS Institute Partnership, Botswana Harvard HIV Reference Laboratory, Gaborone, Botswana
- Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Alexander G Lucaci
- Institute for Genomics and Evolutionary Medicine, Department of Biology, Temple University, Philadelphia, PA, USA
| | - Arisha Maharaj
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), Nelson R. Mandela School of Medicine, University of KwaZulu-Natal, Durban, South Africa
| | - Boitshoko Mahlangu
- National Institute for Communicable Diseases (NICD) of the National Health Laboratory Service (NHLS), Johannesburg, South Africa
| | - Tongai Maponga
- Division of Medical Virology, Faculty of Medicine and Health Sciences, Stellenbosch University, Tygerberg, Cape Town, South Africa
| | - Kamela Mahlakwane
- Division of Medical Virology, Faculty of Medicine and Health Sciences, Stellenbosch University, Tygerberg, Cape Town, South Africa
- NHLS Tygerberg Laboratory, Tygerberg Hospital, Cape Town, South Africa
| | - Zinhle Makatini
- Department of Virology, Charlotte Maxeke Johannesburg Academic Hospital, Johannesburg, South Africa
| | - Gert Marais
- Division of Medical Virology, Faculty of Health Sciences, University of Cape Town, Cape Town, South Africa
- Division of Virology, NHLS Groote Schuur Laboratory, Cape Town, South Africa
| | - Dorcas Maruapula
- Botswana Harvard AIDS Institute Partnership, Botswana Harvard HIV Reference Laboratory, Gaborone, Botswana
| | - Kereng Masupu
- Botswana Presidential COVID-19 Taskforce, Gaborone, Botswana
| | - Mogomotsi Matshaba
- Botswana Presidential COVID-19 Taskforce, Gaborone, Botswana
- Botswana-Baylor Children's Clinical Centre of Excellence, Gaborone, Botswana
- Baylor College of Medicine, Houston, TX, USA
| | - Simnikiwe Mayaphi
- Department of Medical Virology, University of Pretoria, Pretoria, South Africa
| | - Nokuzola Mbhele
- Division of Medical Virology, Faculty of Health Sciences, University of Cape Town, Cape Town, South Africa
| | - Mpaphi B Mbulawa
- National Health Laboratory, Health Services Management, Ministry of Health and Wellness, Gaborone, Botswana
| | - Adriano Mendes
- Zoonotic Arbo and Respiratory Virus Program, Centre for Viral Zoonoses, Department of Medical Virology, University of Pretoria, Pretoria, South Africa
| | - Koleka Mlisana
- National Health Laboratory Service (NHLS), Johannesburg, South Africa
- Centre for the AIDS Programme of Research in South Africa (CAPRISA), Durban, South Africa
| | - Anele Mnguni
- National Institute for Communicable Diseases (NICD) of the National Health Laboratory Service (NHLS), Johannesburg, South Africa
| | - Thabo Mohale
- National Institute for Communicable Diseases (NICD) of the National Health Laboratory Service (NHLS), Johannesburg, South Africa
| | - Monika Moir
- Centre for Epidemic Response and Innovation (CERI), School of Data Science and Computational Thinking, Stellenbosch University, Stellenbosch, South Africa
| | - Kgomotso Moruisi
- Health Services Management, Ministry of Health and Wellness, Gaborone, Botswana
| | - Mosepele Mosepele
- Botswana Presidential COVID-19 Taskforce, Gaborone, Botswana
- Department of Medicine, Faculty of Medicine, University of Botswana, Gaborone, Botswana
| | - Gerald Motsatsi
- National Institute for Communicable Diseases (NICD) of the National Health Laboratory Service (NHLS), Johannesburg, South Africa
| | - Modisa S Motswaledi
- Botswana Presidential COVID-19 Taskforce, Gaborone, Botswana
- Department of Medical Laboratory Sciences, School of Allied Health Professions, Faculty of Health Sciences, University of Botswana, Gaborone, Botswana
| | - Thongbotho Mphoyakgosi
- National Health Laboratory, Health Services Management, Ministry of Health and Wellness, Gaborone, Botswana
| | - Nokukhanya Msomi
- Discipline of Virology, School of Laboratory Medicine and Medical Sciences and National Health Laboratory Service (NHLS), University of KwaZulu-Natal, Durban, South Africa
| | - Peter N Mwangi
- Division of Virology, University of the Free State, Bloemfontein, South Africa
- Next Generation Sequencing Unit, Division of Virology, Faculty of Health Sciences, University of the Free State, Bloemfontein, South Africa
| | - Yeshnee Naidoo
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), Nelson R. Mandela School of Medicine, University of KwaZulu-Natal, Durban, South Africa
| | - Noxolo Ntuli
- National Institute for Communicable Diseases (NICD) of the National Health Laboratory Service (NHLS), Johannesburg, South Africa
| | - Martin Nyaga
- Division of Virology, University of the Free State, Bloemfontein, South Africa
- Next Generation Sequencing Unit, Division of Virology, Faculty of Health Sciences, University of the Free State, Bloemfontein, South Africa
| | - Lucier Olubayo
- Wellcome Centre for Infectious Diseases Research in Africa (CIDRI-Africa), Cape Town, South Africa
- Division of Computational Biology, Faculty of Health Sciences, University of Cape Town, Cape Town, South Africa
| | - Sureshnee Pillay
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), Nelson R. Mandela School of Medicine, University of KwaZulu-Natal, Durban, South Africa
| | - Botshelo Radibe
- Botswana Harvard AIDS Institute Partnership, Botswana Harvard HIV Reference Laboratory, Gaborone, Botswana
| | - Yajna Ramphal
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), Nelson R. Mandela School of Medicine, University of KwaZulu-Natal, Durban, South Africa
| | - Upasana Ramphal
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), Nelson R. Mandela School of Medicine, University of KwaZulu-Natal, Durban, South Africa
| | - James E San
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), Nelson R. Mandela School of Medicine, University of KwaZulu-Natal, Durban, South Africa
| | - Lesley Scott
- Department of Molecular Medicine and Haematology, University of the Witwatersrand, Johannesburg, South Africa
| | - Roger Shapiro
- Botswana Harvard AIDS Institute Partnership, Botswana Harvard HIV Reference Laboratory, Gaborone, Botswana
- Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Lavanya Singh
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), Nelson R. Mandela School of Medicine, University of KwaZulu-Natal, Durban, South Africa
| | | | - Wendy Stevens
- Department of Molecular Medicine and Haematology, University of the Witwatersrand, Johannesburg, South Africa
| | - Amy Strydom
- Zoonotic Arbo and Respiratory Virus Program, Centre for Viral Zoonoses, Department of Medical Virology, University of Pretoria, Pretoria, South Africa
| | - Kathleen Subramoney
- Department of Virology, Charlotte Maxeke Johannesburg Academic Hospital, Johannesburg, South Africa
| | - Naume Tebeila
- National Institute for Communicable Diseases (NICD) of the National Health Laboratory Service (NHLS), Johannesburg, South Africa
| | - Derek Tshiabuila
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), Nelson R. Mandela School of Medicine, University of KwaZulu-Natal, Durban, South Africa
| | - Joseph Tsui
- Department of Zoology, University of Oxford, Oxford, UK
| | - Stephanie van Wyk
- Centre for Epidemic Response and Innovation (CERI), School of Data Science and Computational Thinking, Stellenbosch University, Stellenbosch, South Africa
| | - Steven Weaver
- Institute for Genomics and Evolutionary Medicine, Department of Biology, Temple University, Philadelphia, PA, USA
| | - Constantinos K Wibmer
- National Institute for Communicable Diseases (NICD) of the National Health Laboratory Service (NHLS), Johannesburg, South Africa
| | - Eduan Wilkinson
- Centre for Epidemic Response and Innovation (CERI), School of Data Science and Computational Thinking, Stellenbosch University, Stellenbosch, South Africa
| | - Nicole Wolter
- National Institute for Communicable Diseases (NICD) of the National Health Laboratory Service (NHLS), Johannesburg, South Africa
- School of Pathology, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
| | | | - Boitumelo Zuze
- Botswana Harvard AIDS Institute Partnership, Botswana Harvard HIV Reference Laboratory, Gaborone, Botswana
| | - Dominique Goedhals
- Division of Virology, University of the Free State, Bloemfontein, South Africa
- PathCare Vermaak, Pretoria, South Africa
| | - Wolfgang Preiser
- Division of Medical Virology, Faculty of Medicine and Health Sciences, Stellenbosch University, Tygerberg, Cape Town, South Africa
- NHLS Tygerberg Laboratory, Tygerberg Hospital, Cape Town, South Africa
| | - Florette Treurnicht
- Department of Virology, Charlotte Maxeke Johannesburg Academic Hospital, Johannesburg, South Africa
| | - Marietje Venter
- Zoonotic Arbo and Respiratory Virus Program, Centre for Viral Zoonoses, Department of Medical Virology, University of Pretoria, Pretoria, South Africa
| | - Carolyn Williamson
- Division of Medical Virology, Faculty of Health Sciences, University of Cape Town, Cape Town, South Africa
- Division of Virology, NHLS Groote Schuur Laboratory, Cape Town, South Africa
- Wellcome Centre for Infectious Diseases Research in Africa (CIDRI-Africa), Cape Town, South Africa
- Institute of Infectious Disease and Molecular Medicine, University of Cape Town, Cape Town, South Africa
| | | | - Jinal Bhiman
- National Institute for Communicable Diseases (NICD) of the National Health Laboratory Service (NHLS), Johannesburg, South Africa
- South African Medical Research Council Antibody Immunity Research Unit, School of Pathology, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
| | - Allison Glass
- Lancet Laboratories, Johannesburg, South Africa
- Department of Molecular Pathology, School of Pathology, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
| | - Darren P Martin
- Wellcome Centre for Infectious Diseases Research in Africa (CIDRI-Africa), Cape Town, South Africa
- Institute of Infectious Disease and Molecular Medicine, University of Cape Town, Cape Town, South Africa
| | - Andrew Rambaut
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK
| | - Simani Gaseitsiwe
- Botswana Harvard AIDS Institute Partnership, Botswana Harvard HIV Reference Laboratory, Gaborone, Botswana
- Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - Anne von Gottberg
- National Institute for Communicable Diseases (NICD) of the National Health Laboratory Service (NHLS), Johannesburg, South Africa
- School of Pathology, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
| | - Tulio de Oliveira
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), Nelson R. Mandela School of Medicine, University of KwaZulu-Natal, Durban, South Africa.
- Centre for Epidemic Response and Innovation (CERI), School of Data Science and Computational Thinking, Stellenbosch University, Stellenbosch, South Africa.
- Department of Global Health, University of Washington, Seattle, WA, USA.
| |
Collapse
|
2
|
Deforche K, Vercauteren J, Müller V, Vandamme AM. Behavioral changes before lockdown and decreased retail and recreation mobility during lockdown contributed most to controlling COVID-19 in Western countries. BMC Public Health 2021; 21:654. [PMID: 33823820 PMCID: PMC8022318 DOI: 10.1186/s12889-021-10676-1] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2020] [Accepted: 03/22/2021] [Indexed: 12/15/2022] Open
Abstract
BACKGROUND The COVID-19 pandemic has prompted a lockdown in many countries to control the exponential spread of the SARS-CoV-2 virus, hereby reducing the time-varying basic reproduction number (Rt) to below one. Governments are looking for evidence to balance the demand of their citizens to ease some of the restriction, against the fear of a new peak in infections. In this study, we wanted to quantify the relative contribution of mobility restrictions, and that of behavioral changes that occurred already before the lockdowns, on the reduction of transmission during lockdowns in Western countries in early 2020. METHODS Incidence data of cases and deaths from the first wave of infections for 35 Western countries (32 European, plus Israel, USA and Canada) were analyzed using epidemiological compartment models in a Bayesian framework. Mobility data was used to estimate the timing of changes associated with a lockdown, and was correlated with estimated reductions of Rt. RESULTS Across all countries, the initial median estimate for Rt was 3.6 (95% IQR 2.4-5.2), and it was reduced to 0.78 (95% IQR 0.58-1.01) during lockdown. 48% (18-65%) of the reduction occurred already in the week before lockdown, with lockdown itself causing the remaining drop in transmission. A lower Rt during lockdown was independently associated with an increased time spent at home (0.21 per 10% more time, p < 0.007), and decreased mobility related to retail and recreation (0.07 per 10% less mobility, p < 0.008). CONCLUSIONS In a Western population unaware of the risk, SARS-CoV-2 can be highly contagious with a reproduction number R0 > 5. Our results are consistent with evidence that recreational activities (including restaurant and bar visits) enable super-spreading events. Exiting from lockdown therefore requires continued physical distancing and tight control on this kind of activities.
Collapse
Affiliation(s)
| | - Jurgen Vercauteren
- Department of Microbiology, Immunology and Transplantation, Rega Institute for Medical Research, Clinical and Epidemiological Virology, Katholieke Universiteit Leuven, 3000, Leuven, Belgium.
| | - Viktor Müller
- Insitute of Biology, Eötvös Loránd University, Budapest, Hungary
| | - Anne-Mieke Vandamme
- Department of Microbiology, Immunology and Transplantation, Rega Institute for Medical Research, Clinical and Epidemiological Virology, Katholieke Universiteit Leuven, 3000, Leuven, Belgium
- Center for Global Health and Tropical Medicine, Unidade de Microbiologia, Instituto de Higiene e Medicina Tropical, Universidade Nova de Lisboa, Lisbon, Portugal
| |
Collapse
|
3
|
Cleemput S, Dumon W, Fonseca V, Abdool Karim W, Giovanetti M, Alcantara LC, Deforche K, de Oliveira T. Genome Detective Coronavirus Typing Tool for rapid identification and characterization of novel coronavirus genomes. Bioinformatics 2020; 36:3552-3555. [PMID: 32108862 PMCID: PMC7112083 DOI: 10.1093/bioinformatics/btaa145] [Citation(s) in RCA: 88] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2020] [Revised: 02/21/2020] [Accepted: 02/25/2020] [Indexed: 11/30/2022] Open
Abstract
Summary Genome detective is a web-based, user-friendly software application to quickly and accurately assemble all known virus genomes from next-generation sequencing datasets. This application allows the identification of phylogenetic clusters and genotypes from assembled genomes in FASTA format. Since its release in 2019, we have produced a number of typing tools for emergent viruses that have caused large outbreaks, such as Zika and Yellow Fever Virus in Brazil. Here, we present the Genome Detective Coronavirus Typing Tool that can accurately identify the novel severe acute respiratory syndrome (SARS)-related coronavirus (SARS-CoV-2) sequences isolated in China and around the world. The tool can accept up to 2000 sequences per submission and the analysis of a new whole-genome sequence will take approximately 1 min. The tool has been tested and validated with hundreds of whole genomes from 10 coronavirus species, and correctly classified all of the SARS-related coronavirus (SARSr-CoV) and all of the available public data for SARS-CoV-2. The tool also allows tracking of new viral mutations as the outbreak expands globally, which may help to accelerate the development of novel diagnostics, drugs and vaccines to stop the COVID-19 disease. Availability and implementation https://www.genomedetective.com/app/typingtool/cov Contact koen@emweb.be or deoliveira@ukzn.ac.za Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
| | | | - Vagner Fonseca
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), School of Laboratory Medicine and Medical Sciences, College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa.,Laboratório de Genética Celular e Molecular, Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil.,Coordenação Geral dos Laboratórios de Saúde Pública/Secretaria de Vigilância em Saúde, Ministério da Saúde, (CGLAB/SVS-MS), Brasília, Brazil
| | - Wasim Abdool Karim
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), School of Laboratory Medicine and Medical Sciences, College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Marta Giovanetti
- Laboratório de Flavivírus, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz, Rio de Janeiro, RJ, Brazil
| | - Luiz Carlos Alcantara
- Laboratório de Genética Celular e Molecular, Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil.,Laboratório de Flavivírus, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz, Rio de Janeiro, RJ, Brazil
| | | | - Tulio de Oliveira
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), School of Laboratory Medicine and Medical Sciences, College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa.,Centre for the AIDS Programme of Research in South Africa (CAPRISA), Durban, South Africa.,Department of Global Health, University of Washington, Seattle, WA, USA
| |
Collapse
|
4
|
Libin PJK, Deforche K, Abecasis AB, Theys K. VIRULIGN: fast codon-correct alignment and annotation of viral genomes. Bioinformatics 2020; 35:1763-1765. [PMID: 30295730 PMCID: PMC6513156 DOI: 10.1093/bioinformatics/bty851] [Citation(s) in RCA: 37] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2018] [Revised: 09/24/2018] [Accepted: 10/05/2018] [Indexed: 12/11/2022] Open
Abstract
Summary Virus sequence data are an essential resource for reconstructing spatiotemporal dynamics of viral spread as well as to inform treatment and prevention strategies. However, the potential benefit of these applications critically depends on accurate and correctly annotated alignments of genetically heterogeneous data. VIRULIGN was built for fast codon-correct alignments of large datasets, with standardized and formalized genome annotation and various alignment export formats. Availability and implementation VIRULIGN is freely available at https://github.com/rega-cev/virulign as an open source software project. Supplementary information Supplementary data is available at Bioinformatics online.
Collapse
Affiliation(s)
- Pieter J K Libin
- KU Leuven, Rega Institute for Medical, Laboratorium of Clinical and Evolutionary Virology, Leuven, Belgium.,Artificial Intelligence Lab, Department of Computer Science, Vrije Universiteit Brussel, Brussels, Belgium
| | | | - Ana B Abecasis
- Center for Global Health and Tropical Medicine, Institute for Hygiene and Tropical Medicine, Lisboa, Portugal
| | - Kristof Theys
- KU Leuven, Rega Institute for Medical, Laboratorium of Clinical and Evolutionary Virology, Leuven, Belgium
| |
Collapse
|
5
|
Cleemput S, Dumon W, Fonseca V, Karim WA, Giovanetti M, Alcantara LC, Deforche K, de Oliveira T. Genome Detective Coronavirus Typing Tool for rapid identification and characterization of novel coronavirus genomes. bioRxiv 2020:2020.01.31.928796. [PMID: 32511309 PMCID: PMC7217295 DOI: 10.1101/2020.01.31.928796] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Genome Detective is a web-based, user-friendly software application to quickly and accurately assemble all known virus genomes from next generation sequencing datasets. This application allows the identification of phylogenetic clusters and genotypes from assembled genomes in FASTA format. Since its release in 2019, we have produced a number of typing tools for emergent viruses that have caused large outbreaks, such as Zika and Yellow Fever Virus in Brazil. Here, we present The Genome Detective Coronavirus Typing Tool that can accurately identify novel coronavirus (2019-nCoV) sequences isolated in China and around the world. The tool can accept up to 2,000 sequences per submission and the analysis of a new whole genome sequence will take approximately one minute. The tool has been tested and validated with hundreds of whole genomes from ten coronavirus species, and correctly classified all of the SARS-related coronavirus (SARSr-CoV) and all of the available public data for 2019-nCoV. The tool also allows tracking of new viral mutations as the outbreak expands globally, which may help to accelerate the development of novel diagnostics, drugs and vaccines.
Collapse
Affiliation(s)
| | | | - Vagner Fonseca
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), School of Laboratory Medicine and Medical Sciences, College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
- Laboratório de Genética Celular e Molecular, Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil
- Coordenação Geral dos Laboratórios de Saúde Pública/Secretaria de Vigilância em Saúde, Ministério da Saúde, (CGLAB/SVS-MS) Brasília, Brazil
| | - Wasim Abdool Karim
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), School of Laboratory Medicine and Medical Sciences, College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Marta Giovanetti
- Laboratório de Flavivírus, Instituto Oswaldo CruZ, Fundação Oswaldo CruZ, Rio de Janeiro, RJ, Brazil
| | - Luiz Carlos Alcantara
- Laboratório de Genética Celular e Molecular, Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil
- Laboratório de Flavivírus, Instituto Oswaldo CruZ, Fundação Oswaldo CruZ, Rio de Janeiro, RJ, Brazil
| | | | - Tulio de Oliveira
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), School of Laboratory Medicine and Medical Sciences, College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
- Centre for the AIDS Programme of Research in South Africa (CAPRISA), Durban, South Africa
- Department of Global Health, University of Washington, Seattle, USA
| |
Collapse
|
6
|
Jariani A, Warth C, Deforche K, Libin P, Drummond AJ, Rambaut A, Matsen Iv FA, Theys K. SANTA-SIM: simulating viral sequence evolution dynamics under selection and recombination. Virus Evol 2019; 5:vez003. [PMID: 30863552 PMCID: PMC6407609 DOI: 10.1093/ve/vez003] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Abstract
Simulations are widely used to provide expectations and predictive distributions under known conditions against which to compare empirical data. Such simulations are also invaluable for testing and comparing the behaviour and power of inference methods. We describe SANTA-SIM, a software package to simulate the evolution of a population of gene sequences forwards through time. It models the underlying biological processes as discrete components: replication, recombination, point mutations, insertion-deletions, and selection under various fitness models and population size dynamics. The software is designed to be intuitive to work with for a wide range of users and executable in a cross-platform manner.
Collapse
Affiliation(s)
- Abbas Jariani
- Laboratory for Genetics and Genomics, Center of Microbial and Plant Genetics, KU Leuven, Leuven, Belgium.,Laboratory for Systems Biology, VIB, Leuven, Belgium
| | - Christopher Warth
- Program in Computational Biology, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | | | - Pieter Libin
- Laboratory Clinical and Evolutionary Virology, Rega Institute for Medical Research - KU Leuven, Leuven, Belgium.,Artifical Intelligence Lab, Department of Computer Science, Vrije Universiteit Brussel, Brussels, Belgium
| | - Alexei J Drummond
- Centre for Computational Evolution, University of Auckland, Auckland, New Zealand.,Department of Biosystems Science and Engineering, Eidgenössische Technische Hochschule Zurich, Basel, Switzerland
| | - Andrew Rambaut
- Ashworth Laboratories, Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK
| | - Frederick A Matsen Iv
- Program in Computational Biology, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | - Kristof Theys
- Laboratory Clinical and Evolutionary Virology, Rega Institute for Medical Research - KU Leuven, Leuven, Belgium
| |
Collapse
|
7
|
Vilsker M, Moosa Y, Nooij S, Fonseca V, Ghysens Y, Dumon K, Pauwels R, Alcantara LC, Vanden Eynden E, Vandamme AM, Deforche K, de Oliveira T. Genome Detective: an automated system for virus identification from high-throughput sequencing data. Bioinformatics 2019; 35:871-873. [PMID: 30124794 PMCID: PMC6524403 DOI: 10.1093/bioinformatics/bty695] [Citation(s) in RCA: 201] [Impact Index Per Article: 40.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2018] [Revised: 07/09/2018] [Accepted: 08/14/2018] [Indexed: 12/20/2022] Open
Abstract
SUMMARY Genome Detective is an easy to use web-based software application that assembles the genomes of viruses quickly and accurately. The application uses a novel alignment method that constructs genomes by reference-based linking of de novo contigs by combining amino-acids and nucleotide scores. The software was optimized using synthetic datasets to represent the great diversity of virus genomes. The application was then validated with next generation sequencing data of hundreds of viruses. User time is minimal and it is limited to the time required to upload the data. AVAILABILITY AND IMPLEMENTATION Available online: http://www.genomedetective.com/app/typingtool/virus/. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
| | - Yumna Moosa
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), School of Laboratory Medicine and Medical Sciences, Nelson R Mandela School of Medicine, College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Sam Nooij
- The Dutch National Institute for Public Health and the Environment (RIVM), Bilthoven, The Netherlands
| | - Vagner Fonseca
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), School of Laboratory Medicine and Medical Sciences, Nelson R Mandela School of Medicine, College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
- Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais, Belo Horizonte, MG, Brazil
- Laboratory of Hematology Genetic and computational Biology, Goncalo Moniz Research Center, Oswaldo Cruz Foundation (LHGB/CPqGM/FIOCRUZ), Bahia, Brazil
| | | | | | | | - Luiz Carlos Alcantara
- Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais, Belo Horizonte, MG, Brazil
- Laboratory of Hematology Genetic and computational Biology, Goncalo Moniz Research Center, Oswaldo Cruz Foundation (LHGB/CPqGM/FIOCRUZ), Bahia, Brazil
- Laboratório de Flavivírus, IOC, Fundação Oswaldo Cruz
| | - Ewout Vanden Eynden
- KU Leuven, Department of Microbiology and Immunology, Rega Institute for Medical Research, Clinical and Epidemiological Virology, Leuven, Belgium
| | - Anne-Mieke Vandamme
- KU Leuven, Department of Microbiology and Immunology, Rega Institute for Medical Research, Clinical and Epidemiological Virology, Leuven, Belgium
- Center for Global Health and Tropical Medicine, Unidade de Microbiologia, Instituto de Higiene e Medicina Tropical, Universidade Nova de Lisboa, Lisbon, Portugal
| | | | - Tulio de Oliveira
- KwaZulu-Natal Research Innovation and Sequencing Platform (KRISP), School of Laboratory Medicine and Medical Sciences, Nelson R Mandela School of Medicine, College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| |
Collapse
|
8
|
Moosa Y, Vilsker M, Vanden Eyden E, Fonseca V, Nooij S, Deforche K, de Oliveira T. A34 Automated profiling of the human virome from raw metagenomic data. Virus Evol 2017; 3:vew036.033. [PMID: 28845264 PMCID: PMC5565974 DOI: 10.1093/ve/vew036.033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yumna Moosa
- Nelson R Mandela School of Medicine College of Health Sciences, University of KwaZulu-Natal (UKZN), South Africa
| | | | | | | | - Sam Nooij
- Rijksinstituut voor Volksgezondheid en Milieu (RIVM), The Netherlands
| | | | - Tulio de Oliveira
- Nelson R Mandela School of Medicine College of Health Sciences, University of KwaZulu-Natal (UKZN), South Africa
| |
Collapse
|
9
|
Vilsker M, Deforche K. A16 A distributed pan-viral typing framework. Virus Evol 2017; 3:vew036.015. [PMID: 28845257 PMCID: PMC5565943 DOI: 10.1093/ve/vew036.015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
|
10
|
Armitage AE, Deforche K, Welch JJ, Van Laethem K, Camacho R, Rambaut A, Iversen AKN. Possible footprints of APOBEC3F and/or other APOBEC3 deaminases, but not APOBEC3G, on HIV-1 from patients with acute/early and chronic infections. J Virol 2014; 88:12882-94. [PMID: 25165112 PMCID: PMC4248940 DOI: 10.1128/jvi.01460-14] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2014] [Accepted: 08/21/2014] [Indexed: 12/19/2022] Open
Abstract
UNLABELLED Members of the apolipoprotein B mRNA-editing enzyme-catalytic polypeptide-like-3 (APOBEC3) innate cellular cytidine deaminase family, particularly APOBEC3F and APOBEC3G, can cause extensive and lethal G-to-A mutations in HIV-1 plus-strand DNA (termed hypermutation). It is unclear if APOBEC3-induced mutations in vivo are always lethal or can occur at sublethal levels that increase HIV-1 diversification and viral adaptation to the host. The viral accessory protein Vif counteracts APOBEC3 activity by binding to APOBEC3 and promoting proteasome degradation; however, the efficiency of this interaction varies, since a range of hypermutation frequencies are observed in HIV-1 patient DNA. Therefore, we examined "footprints" of APOBEC3G and APOBEC3F activity in longitudinal HIV-1 RNA pol sequences from approximately 3,000 chronically infected patients by determining whether G-to-A mutations occurred in motifs that were favored or disfavored by these deaminases. G-to-A mutations were more frequent in APOBEC3G-disfavored than in APOBEC3G-favored contexts. In contrast, mutations in APOBEC3F-disfavored contexts were relatively rare, whereas mutations in contexts favoring APOBEC3F (and possibly other deaminases) occurred 16% more often than average G-to-A mutations. These results were supported by analyses of >500 HIV-1 env sequences from acute/early infection. IMPORTANCE Collectively, our results suggest that APOBEC3G-induced mutagenesis is lethal to HIV-1, whereas mutagenesis caused by APOBEC3F and/or other deaminases may result in sublethal mutations that might facilitate viral diversification. Therefore, Vif-specific cytotoxic T lymphocyte (CTL) responses and drugs that manipulate the interplay between Vif and APOBEC3 may have beneficial or detrimental clinical effects depending on how they affect the binding of Vif to various members of the APOBEC3 family.
Collapse
Affiliation(s)
- Andrew E Armitage
- MRC Human Immunology Unit, Weatherall Institute of Molecular Medicine, Oxford University, Oxford, United Kingdom
| | - Koen Deforche
- KU Leuven-University of Leuven, Department of Microbiology and Immunology, Rega Institute for Medical Research, Leuven, Belgium
| | - John J Welch
- Department of Genetics, University of Cambridge, Cambridge, United Kingdom
| | - Kristel Van Laethem
- KU Leuven-University of Leuven, Department of Microbiology and Immunology, Rega Institute for Medical Research, Leuven, Belgium
| | - Ricardo Camacho
- KU Leuven-University of Leuven, Department of Microbiology and Immunology, Rega Institute for Medical Research, Leuven, Belgium Centro de Malária e outras Doenças Tropicais, Instituto de Higiene e Medicina Tropical, Universidade Nova de Lisboa, Lisbon, Portugal
| | - Andrew Rambaut
- Institute of Evolutionary Biology. University of Edinburgh, Edinburgh, United Kingdom
| | - Astrid K N Iversen
- MRC Human Immunology Unit, Weatherall Institute of Molecular Medicine, Oxford University, Oxford, United Kingdom Nuffield Department of Clinical Neurosciences, Division of Clinical Neurology, Weatherall Institute of Molecular Medicine, Oxford University, Oxford, United Kingdom
| |
Collapse
|
11
|
Pineda-Peña AC, Faria NR, Imbrechts S, Libin P, Abecasis AB, Deforche K, Gómez-López A, Camacho RJ, de Oliveira T, Vandamme AM. Automated subtyping of HIV-1 genetic sequences for clinical and surveillance purposes: performance evaluation of the new REGA version 3 and seven other tools. Infect Genet Evol 2013; 19:337-48. [PMID: 23660484 DOI: 10.1016/j.meegid.2013.04.032] [Citation(s) in RCA: 275] [Impact Index Per Article: 25.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/09/2012] [Revised: 04/10/2013] [Accepted: 04/28/2013] [Indexed: 12/26/2022]
Abstract
BACKGROUND To investigate differences in pathogenesis, diagnosis and resistance pathways between HIV-1 subtypes, an accurate subtyping tool for large datasets is needed. We aimed to evaluate the performance of automated subtyping tools to classify the different subtypes and circulating recombinant forms using pol, the most sequenced region in clinical practice. We also present the upgraded version 3 of the Rega HIV subtyping tool (REGAv3). METHODOLOGY HIV-1 pol sequences (PR+RT) for 4674 patients retrieved from the Portuguese HIV Drug Resistance Database, and 1872 pol sequences trimmed from full-length genomes retrieved from the Los Alamos database were classified with statistical-based tools such as COMET, jpHMM and STAR; similarity-based tools such as NCBI and Stanford; and phylogenetic-based tools such as REGA version 2 (REGAv2), REGAv3, and SCUEAL. The performance of these tools, for pol, and for PR and RT separately, was compared in terms of reproducibility, sensitivity and specificity with respect to the gold standard which was manual phylogenetic analysis of the pol region. RESULTS The sensitivity and specificity for subtypes B and C was more than 96% for seven tools, but was variable for other subtypes such as A, D, F and G. With regard to the most common circulating recombinant forms (CRFs), the sensitivity and specificity for CRF01_AE was ~99% with statistical-based tools, with phylogenetic-based tools and with Stanford, one of the similarity based tools. CRF02_AG was correctly identified for more than 96% by COMET, REGAv3, Stanford and STAR. All the tools reached a specificity of more than 97% for most of the subtypes and the two main CRFs (CRF01_AE and CRF02_AG). Other CRFs were identified only by COMET, REGAv2, REGAv3, and SCUEAL and with variable sensitivity. When analyzing sequences for PR and RT separately, the performance for PR was generally lower and variable between the tools. Similarity and statistical-based tools were 100% reproducible, but this was lower for phylogenetic-based tools such as REGA (~99%) and SCUEAL (~96%). CONCLUSIONS REGAv3 had an improved performance for subtype B and CRF02_AG compared to REGAv2 and is now able to also identify all epidemiologically relevant CRFs. In general the best performing tools, in alphabetical order, were COMET, jpHMM, REGAv3, and SCUEAL when analyzing pure subtypes in the pol region, and COMET and REGAv3 when analyzing most of the CRFs. Based on this study, we recommend to confirm subtyping with 2 well performing tools, and be cautious with the interpretation of short sequences.
Collapse
Affiliation(s)
- Andrea-Clemencia Pineda-Peña
- Laboratory for Clinical and Epidemiological Virology, Rega Institute for Medical Research, Department of Microbiology and Immunology, University of Leuven, Belgium; Clinical and Molecular Infectious Diseases Group, Faculty of Sciences and Mathematics, Universidad del Rosario, Bogotá, Colombia.
| | | | | | | | | | | | | | | | | | | |
Collapse
|
12
|
Libin P, Beheydt G, Deforche K, Imbrechts S, Ferreira F, Van Laethem K, Theys K, Carvalho AP, Cavaco-Silva J, Lapadula G, Torti C, Assel M, Wesner S, Snoeck J, Ruelle J, De Bel A, Lacor P, De Munter P, Van Wijngaerden E, Zazzi M, Kaiser R, Ayouba A, Peeters M, de Oliveira T, Alcantara LCJ, Grossman Z, Sloot P, Otelea D, Paraschiv S, Boucher C, Camacho RJ, Vandamme AM. RegaDB: community-driven data management and analysis for infectious diseases. Bioinformatics 2013; 29:1477-80. [PMID: 23645815 PMCID: PMC3661054 DOI: 10.1093/bioinformatics/btt162] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Summary: RegaDB is a free and open source data management and analysis environment for infectious diseases. RegaDB allows clinicians to store, manage and analyse patient data, including viral genetic sequences. Moreover, RegaDB provides researchers with a mechanism to collect data in a uniform format and offers them a canvas to make newly developed bioinformatics tools available to clinicians and virologists through a user friendly interface. Availability and implementation: Source code, binaries and documentation are available on http://rega.kuleuven.be/cev/regadb. RegaDB is written in the Java programming language, using a web-service-oriented architecture. Contact:pieter.libin@rega.kuleuven.be
Collapse
Affiliation(s)
- Pieter Libin
- Department of Microbiology and Immunology, Rega Institute for Medical Research, KU Leuven, Leuven, Belgium.
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
13
|
Sangeda RZ, Theys K, Beheydt G, Rhee SY, Deforche K, Vercauteren J, Libin P, Imbrechts S, Grossman Z, Camacho RJ, Van Laethem K, Pironti A, Zazzi M, Sönnerborg A, Incardona F, De Luca A, Torti C, Ruiz L, Van de Vijver DAMC, Shafer RW, Bruzzone B, Van Wijngaerden E, Vandamme AM. HIV-1 fitness landscape models for indinavir treatment pressure using observed evolution in longitudinal sequence data are predictive for treatment failure. Infect Genet Evol 2013; 19:349-60. [PMID: 23523594 DOI: 10.1016/j.meegid.2013.03.014] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/09/2012] [Revised: 02/25/2013] [Accepted: 03/04/2013] [Indexed: 11/29/2022]
Abstract
We previously modeled the in vivo evolution of human immunodeficiency virus-1 (HIV-1) under drug selective pressure from cross-sectional viral sequences. These fitness landscapes (FLs) were made by using first a Bayesian network (BN) to map epistatic substitutions, followed by scaling the fitness landscape based on an HIV evolution simulator trying to evolve the sequences from treatment naïve patients into sequences from patients failing treatment. In this study, we compared four FLs trained with different sequence populations. Epistatic interactions were learned from three different cross-sectional BNs, trained with sequence from patients experienced with indinavir (BNT), all protease inhibitors (PIs) (BNP) or all PI except indinavir (BND). Scaling the fitness landscape was done using cross-sectional data from drug naïve and indinavir experienced patients (Fcross using BNT) and using longitudinal sequences from patients failing indinavir (FlongT using BNT, FlongP using BNP, FlongD using BND). Evaluation to predict the failing sequence and therapy outcome was performed on independent sequences of patients on indinavir. Parameters included estimated fitness (LogF), the number of generations (GF) or mutations (MF) to reach the fitness threshold (average fitness when a major resistance mutation appeared), the number of generations (GR) or mutations (MR) to reach a major resistance mutation and compared to genotypic susceptibility score (GSS) from Rega and HIVdb algorithms. In pairwise FL comparisons we found significant correlation between fitness values for individual sequences, and this correlation improved after correcting for the subtype. Furthermore, FLs could predict the failing sequence under indinavir-containing combinations. At 12 and 48 weeks, all parameters from all FLs and indinavir GSS (both for Rega and HIVdb) were predictive of therapy outcome, except MR for FlongT and FlongP. The fitness landscapes have similar predictive power for treatment response under indinavir-containing regimen as standard rules-based algorithms, and additionally allow predicting genetic evolution under indinavir selective pressure.
Collapse
Affiliation(s)
- Raphael Z Sangeda
- Rega Institute for Medical Research, Department of Microbiology and Immunology, KU Leuven, 3000 Leuven, Belgium
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
14
|
Theys K, Deforche K, Vercauteren J, Libin P, van de Vijver DAMC, Albert J, Åsjö B, Balotta C, Bruckova M, Camacho RJ, Clotet B, Coughlan S, Grossman Z, Hamouda O, Horban A, Korn K, Kostrikis LG, Kücherer C, Nielsen C, Paraskevis D, Poljak M, Puchhammer-Stockl E, Riva C, Ruiz L, Liitsola K, Schmit JC, Schuurman R, Sönnerborg A, Stanekova D, Stanojevic M, Struck D, Van Laethem K, Wensing AMJ, Boucher CAB, Vandamme AM. Treatment-associated polymorphisms in protease are significantly associated with higher viral load and lower CD4 count in newly diagnosed drug-naive HIV-1 infected patients. Retrovirology 2012; 9:81. [PMID: 23031662 PMCID: PMC3487874 DOI: 10.1186/1742-4690-9-81] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2012] [Accepted: 08/23/2012] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND The effect of drug resistance transmission on disease progression in the newly infected patient is not well understood. Major drug resistance mutations severely impair viral fitness in a drug free environment, and therefore are expected to revert quickly. Compensatory mutations, often already polymorphic in wild-type viruses, do not tend to revert after transmission. While compensatory mutations increase fitness during treatment, their presence may also modulate viral fitness and virulence in absence of therapy and major resistance mutations. We previously designed a modeling technique that quantifies genotypic footprints of in vivo treatment selective pressure, including both drug resistance mutations and polymorphic compensatory mutations, through the quantitative description of a fitness landscape from virus genetic sequences. RESULTS Genotypic correlates of viral load and CD4 cell count were evaluated in subtype B sequences from recently diagnosed treatment-naive patients enrolled in the SPREAD programme. The association of surveillance drug resistance mutations, reported compensatory mutations and fitness estimated from drug selective pressure fitness landscapes with baseline viral load and CD4 cell count was evaluated using regression techniques. Protease genotypic variability estimated to increase fitness during treatment was associated with higher viral load and lower CD4 cell counts also in treatment-naive patients, which could primarily be attributed to well-known compensatory mutations at highly polymorphic positions. By contrast, treatment-related mutations in reverse transcriptase could not explain viral load or CD4 cell count variability. CONCLUSIONS These results suggest that polymorphic compensatory mutations in protease, reported to be selected during treatment, may improve the replicative capacity of HIV-1 even in absence of drug selective pressure or major resistance mutations. The presence of this polymorphic variation may either reflect a history of drug selective pressure, i.e. transmission from a treated patient, or merely be a result of diversity in wild-type virus. Our findings suggest that transmitted drug resistance has the potential to contribute to faster disease progression in the newly infected host and to shape the HIV-1 epidemic at a population level.
Collapse
Affiliation(s)
- Kristof Theys
- Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium
| | | | - Jurgen Vercauteren
- Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium
| | | | | | - Jan Albert
- Clinical Microbiology, Karolinska University Hospital and Department of Microbiology, Tumor and Cell Biology, Karolinska Institute, Stockholm, Sweden
| | - Birgitta Åsjö
- Section for Microbiology and Immunology, Gade institute, University of Bergen, Bergen, Norway
| | | | - Marie Bruckova
- National Institute of Public Health, Prague, Czech Republic
| | - Ricardo J Camacho
- Centro de Malária e outras Doenças Tropicais, Instituto de Higiene e Medicina Tropical, Universidade Nova de Lisboa, Lisbon, Portugal
- Centro Hospitalar de Lisboa Ocidental, Lisbon, Portugal
| | - Bonaventura Clotet
- irsiCaixa AIDS Research Institute & Lluita contra la SIDA Foundation, Hospital Universitari “Germans Trias i Pujol”, Badalona, Spain
| | | | - Zehava Grossman
- Sheba Medical Center, Tel-Hashomer, and School of Public Health, Tel-Aviv University, Tel-Aviv, Israel
| | | | - Andrzei Horban
- Warsaw Medical University and Hospital for Infectious Diseases, Warsaw, Poland
| | - Klaus Korn
- Institut für Klinische und Molekulare Virologie, University of Erlangen, Erlangen, Germany
| | | | | | | | - Dimitrios Paraskevis
- National Retrovirus Reference Center, Department of Hygiene Epidemiology of Medical Statistics, University of Athens, Medical School, Athens, Greece
| | | | | | | | - Lidia Ruiz
- irsiCaixa AIDS Research Institute & Lluita contra la SIDA Foundation, Hospital Universitari “Germans Trias i Pujol”, Badalona, Spain
| | - Kirsi Liitsola
- National Institute of Health and Welfare, Helsinki, Finland
| | - Jean-Claude Schmit
- Centre Hospitalier de Luxembourg and Centre de Recherche Public de la Santé, Luxembourg, Luxembourg
| | - Rob Schuurman
- Department of Medical Microbiology, University Medical Center Utrecht, Utrecht, the Netherland
| | - Anders Sönnerborg
- Divisions of Infectious Diseases and Clinical Virology, Karolinska Institutet, Stockholm, Sweden
| | | | - Maja Stanojevic
- School of Medicine, University of Belgrade, Belgrade, Serbia
| | - Daniel Struck
- Centre Hospitalier de Luxembourg and Centre de Recherche Public de la Santé, Luxembourg, Luxembourg
| | - Kristel Van Laethem
- Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium
| | - Annemarie MJ Wensing
- Department of Medical Microbiology, University Medical Center Utrecht, Utrecht, the Netherland
| | - Charles AB Boucher
- Department of Virology, Erasmus Medical Center, Rotterdam, the Netherlands
- Department of Medical Microbiology, University Medical Center Utrecht, Utrecht, the Netherland
| | - Anne-Mieke Vandamme
- Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium
- Centro de Malária e outras Doenças Tropicais, Instituto de Higiene e Medicina Tropical, Universidade Nova de Lisboa, Lisbon, Portugal
| |
Collapse
|
15
|
Armitage AE, Deforche K, Chang CH, Wee E, Kramer B, Welch JJ, Gerstoft J, Fugger L, McMichael A, Rambaut A, Iversen AKN. APOBEC3G-induced hypermutation of human immunodeficiency virus type-1 is typically a discrete "all or nothing" phenomenon. PLoS Genet 2012; 8:e1002550. [PMID: 22457633 PMCID: PMC3310730 DOI: 10.1371/journal.pgen.1002550] [Citation(s) in RCA: 56] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2011] [Accepted: 01/07/2012] [Indexed: 11/18/2022] Open
Abstract
The rapid evolution of Human Immunodeficiency Virus (HIV-1) allows studies of ongoing host-pathogen interactions. One key selective host factor is APOBEC3G (hA3G) that can cause extensive and inactivating Guanosine-to-Adenosine (G-to-A) mutation on HIV plus-strand DNA (termed hypermutation). HIV can inhibit this innate anti-viral defense through binding of the viral protein Vif to hA3G, but binding efficiency varies and hypermutation frequencies fluctuate in patients. A pivotal question is whether hA3G-induced G-to-A mutation is always lethal to the virus or if it may occur at sub-lethal frequencies that could increase viral diversification. We show in vitro that limiting-levels of hA3G-activity (i.e. when only a single hA3G-unit is likely to act on HIV) produce hypermutation frequencies similar to those in patients and demonstrate in silico that potentially non-lethal G-to-A mutation rates are ∼10-fold lower than the lowest observed hypermutation levels in vitro and in vivo. Our results suggest that even a single incorporated hA3G-unit is likely to cause extensive and inactivating levels of HIV hypermutation and that hypermutation therefore is typically a discrete "all or nothing" phenomenon. Thus, therapeutic measures that inhibit the interaction between Vif and hA3G will likely not increase virus diversification but expand the fraction of hypermutated proviruses within the infected host.
Collapse
Affiliation(s)
- Andrew E. Armitage
- MRC Human Immunology Unit, Weatherall Institute of Molecular Medicine, Oxford University, Oxford, United Kingdom
| | - Koen Deforche
- Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium
| | - Chih-hao Chang
- MRC Human Immunology Unit, Weatherall Institute of Molecular Medicine, Oxford University, Oxford, United Kingdom
| | - Edmund Wee
- MRC Human Immunology Unit, Weatherall Institute of Molecular Medicine, Oxford University, Oxford, United Kingdom
| | - Beatrice Kramer
- Department of Infectious Diseases, King's College London School of Medicine, London, United Kingdom
| | - John J. Welch
- Department of Genetics, University of Cambridge, Cambridge, United Kingdom
| | - Jan Gerstoft
- Department of Infectious Diseases, Rigshospitalet, The National University Hospital, Copenhagen, Denmark
| | - Lars Fugger
- MRC Human Immunology Unit, Weatherall Institute of Molecular Medicine, Oxford University, Oxford, United Kingdom
- Department of Clinical Neurology, John Radcliffe Hospital, Oxford University, Oxford, United Kingdom
| | - Andrew McMichael
- MRC Human Immunology Unit, Weatherall Institute of Molecular Medicine, Oxford University, Oxford, United Kingdom
| | - Andrew Rambaut
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, United Kingdom
- * E-mail: (AKNI); (AR)
| | - Astrid K. N. Iversen
- MRC Human Immunology Unit, Weatherall Institute of Molecular Medicine, Oxford University, Oxford, United Kingdom
- The Peter Medawar Building for Pathogen Research, Nuffield Department of Medicine, Oxford University, Oxford, United Kingdom
- * E-mail: (AKNI); (AR)
| |
Collapse
|
16
|
Araujo THA, Souza-Brito LI, Libin P, Vandamme AM, Deforche K, Galvao-Castro B, Alcantara LC. Development of HTLV-1 Database for data management and datamining HTLV-1 sequences. Retrovirology 2011. [PMCID: PMC3112795 DOI: 10.1186/1742-4690-8-s1-a79] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
|
17
|
Kroneman A, Vennema H, Deforche K, v d Avoort H, Peñaranda S, Oberste MS, Vinjé J, Koopmans M. An automated genotyping tool for enteroviruses and noroviruses. J Clin Virol 2011; 51:121-5. [PMID: 21514213 DOI: 10.1016/j.jcv.2011.03.006] [Citation(s) in RCA: 580] [Impact Index Per Article: 44.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2011] [Accepted: 03/21/2011] [Indexed: 10/18/2022]
Abstract
BACKGROUND Molecular techniques are established as routine in virological laboratories and virus typing through (partial) sequence analysis is increasingly common. Quality assurance for the use of typing data requires harmonization of genotype nomenclature, and agreement on target genes, depending on the level of resolution required, and robustness of methods. OBJECTIVE To develop and validate web-based open-access typing-tools for enteroviruses and noroviruses. STUDY DESIGN An automated web-based typing algorithm was developed, starting with BLAST analysis of the query sequence against a reference set of sequences from viruses in the family Picornaviridae or Caliciviridae. The second step is phylogenetic analysis of the query sequence and a sub-set of the reference sequences, to assign the enterovirus type or norovirus genotype and/or variant, with profile alignment, construction of phylogenetic trees and bootstrap validation. Typing is performed on VP1 sequences of Human enterovirus A to D, and ORF1 and ORF2 sequences of genogroup I and II noroviruses. For validation, we used the tools to automatically type sequences in the RIVM and CDC enterovirus databases and the FBVE norovirus database. RESULTS Using the typing-tools, 785(99%) of 795 Enterovirus VP1 sequences, and 8154(98.5%) of 8342 norovirus sequences were typed in accordance with previously used methods. Subtyping into variants was achieved for 4439(78.4%) of 5838 NoV GII.4 sequences. DISCUSSION AND CONCLUSIONS The online typing-tools reliably assign genotypes for enteroviruses and noroviruses. The use of phylogenetic methods makes these tools robust to ongoing evolution. This should facilitate standardized genotyping and nomenclature in clinical and public health laboratories, thus supporting inter-laboratory comparisons.
Collapse
Affiliation(s)
- A Kroneman
- Laboratory of infectious diseases, National Institute for Public Health and the Environment, Antonie van Leeuwenhoeklaan 9, 3720BA Bilthoven, The Netherlands.
| | | | | | | | | | | | | | | |
Collapse
|
18
|
Theys K, Deforche K, Beheydt G, Moreau Y, van Laethem K, Lemey P, Camacho RJ, Rhee SY, Shafer RW, Van Wijngaerden E, Vandamme AM. Estimating the individualized HIV-1 genetic barrier to resistance using a nelfinavir fitness landscape. BMC Bioinformatics 2010; 11:409. [PMID: 20682040 PMCID: PMC2921410 DOI: 10.1186/1471-2105-11-409] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2009] [Accepted: 08/03/2010] [Indexed: 11/10/2022] Open
Abstract
Background Failure on Highly Active Anti-Retroviral Treatment is often accompanied with development of antiviral resistance to one or more drugs included in the treatment. In general, the virus is more likely to develop resistance to drugs with a lower genetic barrier. Previously, we developed a method to reverse engineer, from clinical sequence data, a fitness landscape experienced by HIV-1 under nelfinavir (NFV) treatment. By simulation of evolution over this landscape, the individualized genetic barrier to NFV resistance may be estimated for an isolate. Results We investigated the association of estimated genetic barrier with risk of development of NFV resistance at virological failure, in 201 patients that were predicted fully susceptible to NFV at baseline, and found that a higher estimated genetic barrier was indeed associated with lower odds for development of resistance at failure (OR 0.62 (0.45 - 0.94), per additional mutation needed, p = .02). Conclusions Thus, variation in individualized genetic barrier to NFV resistance may impact effective treatment options available after treatment failure. If similar results apply for other drugs, then estimated genetic barrier may be a new clinical tool for choice of treatment regimen, which allows consideration of available treatment options after virological failure.
Collapse
Affiliation(s)
- Kristof Theys
- Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium.
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
19
|
Theys K, Deforche K, Libin P, Camacho RJ, Van Laethem K, Vandamme AM. Resistance pathways of human immunodeficiency virus type 1 against the combination of zidovudine and lamivudine. J Gen Virol 2010; 91:1898-1908. [PMID: 20410311 DOI: 10.1099/vir.0.022657-0] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
A better understanding of human immunodeficiency virus type 1 drug-resistance evolution under the selective pressure of combination treatment is important for the design of long-term effective treatment strategies. We applied Bayesian network learning to sequences from patients treated with the reverse transcriptase inhibitor combination of zidovudine (AZT) and lamivudine (3TC) to identify the role of many treatment-selected mutations in the development of resistance. Based on the Bayesian network structure, an in vivo fitness landscape was built, reflecting the necessary selective pressure under treatment, to evolve naive sequences to sequences obtained from patients treated with the combination. This landscape, combined with an evolutionary model, was used to predict resistance evolution in longitudinal sequence pairs. In our analysis, mutations 41L, 70R, 184V and 215F/Y were identified as major resistance mutations to the combination of AZT and 3TC, as they were associated directly with treatment experience. The network also suggested a possible role in resistance development for a number of novel mutations. Estimated fitness, using the landscape, correlated significantly with in vitro resistance phenotype in genotype-phenotype pairs (R(2)=0.70). Variation in predicted evolution under selective pressure correlated significantly with observed in vivo evolution during AZT plus 3CT treatment. In conclusion, we confirmed current knowledge on resistance development to the combination of AZT and 3CT, but additional novel mutations were identified. Moreover, a model to predict resistance evolution during AZT and 3CT treatment has been built and validated.
Collapse
Affiliation(s)
- K Theys
- Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium
| | | | - P Libin
- MyBioData, Rotselaar, Belgium
| | - R J Camacho
- Instituto de Higiene e Medicina Tropical, Universidade Nova de Lisboa, Lisbon, Portugal
| | - K Van Laethem
- Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium
| | - A-M Vandamme
- Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium
| |
Collapse
|
20
|
Alcantara LCJ, Cassol S, Libin P, Deforche K, Pybus OG, Van Ranst M, Galvão-Castro B, Vandamme AM, de Oliveira T. A standardized framework for accurate, high-throughput genotyping of recombinant and non-recombinant viral sequences. Nucleic Acids Res 2009; 37:W634-42. [PMID: 19483099 PMCID: PMC2703899 DOI: 10.1093/nar/gkp455] [Citation(s) in RCA: 133] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Human immunodeficiency virus type-1 (HIV-1), hepatitis B and C and other rapidly evolving viruses are characterized by extremely high levels of genetic diversity. To facilitate diagnosis and the development of prevention and treatment strategies that efficiently target the diversity of these viruses, and other pathogens such as human T-lymphotropic virus type-1 (HTLV-1), human herpes virus type-8 (HHV8) and human papillomavirus (HPV), we developed a rapid high-throughput-genotyping system. The method involves the alignment of a query sequence with a carefully selected set of pre-defined reference strains, followed by phylogenetic analysis of multiple overlapping segments of the alignment using a sliding window. Each segment of the query sequence is assigned the genotype and sub-genotype of the reference strain with the highest bootstrap (>70%) and bootscanning (>90%) scores. Results from all windows are combined and displayed graphically using color-coded genotypes. The new Virus-Genotyping Tools provide accurate classification of recombinant and non-recombinant viruses and are currently being assessed for their diagnostic utility. They have incorporated into several HIV drug resistance algorithms including the Stanford (http://hivdb.stanford.edu) and two European databases (http://www.umcutrecht.nl/subsite/spread-programme/ and http://www.hivrdb.org.uk/) and have been successfully used to genotype a large number of sequences in these and other databases. The tools are a PHP/JAVA web application and are freely accessible on a number of servers including: http://bioafrica.mrc.ac.za/rega-genotype/html/, http://lasp.cpqgm.fiocruz.br/virus-genotype/html/, http://jose.med.kuleuven.be/genotypetool/html/.
Collapse
|
21
|
Theys K, Vercauteren J, Abecasis AB, Libin P, Deforche K, Vandamme AM, Camacho R. The rise and fall of K65R in a Portuguese HIV-1 Drug Resistance database, despite continuously increasing use of tenofovir. Infect Genet Evol 2008; 9:683-8. [PMID: 19038365 DOI: 10.1016/j.meegid.2008.10.013] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/17/2008] [Revised: 09/09/2008] [Accepted: 10/26/2008] [Indexed: 11/16/2022]
Abstract
OBJECTIVE The overall prevalence of the K65R mutation in HIV-1 reverse transcriptase has increased in treatment-experienced patients, mostly attributed to the increasing use of tenofovir (TDF). A number of TDF-based regimens are associated with high rate of early virological failure. In this study, we evaluated the impact of these combinations on K65R selection over time. METHODS Treatment-experienced patients who had a genotypic resistance test at time of failure in the period 2002-2005 in a hospital in Lisbon were included. Incidence of K65R was calculated and compared with proportions of failing therapies in the respective year of sampling including TDF or didanosine plus stavudine and their respective frequency of K65R selection were analyzed using classical statistics and Bayesian Network Learning. RESULTS The overall rate of K65R in the database was consistent with earlier reports. The incidence of K65R increased in the period 2002-2004 but showed a sharp, significant decrease in 2005, despite a continuous increase of patients failing TDF-based regimens. The proportion of patients failing certain not-recommended regimens decreased sharply in 2005 and therefore correlated well with K65R incidence trend. CONCLUSIONS This study suggests that the use of certain TDF-based regimens caused the increase in K65R incidence over the last years.
Collapse
Affiliation(s)
- Kristof Theys
- Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium.
| | | | | | | | | | | | | |
Collapse
|
22
|
Deforche K, Cozzi-Lepri A, Theys K, Clotet B, Camacho RJ, Kjaer J, Van Laethem K, Phillips A, Moreau Y, Lundgren JD, Vandamme AM. Modelled in vivo HIV Fitness under drug Selective Pressure and Estimated Genetic Barrier Towards Resistance are Predictive for Virological Response. Antivir Ther 2008. [DOI: 10.1177/135965350801300316] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Background A method has been developed to estimate a fitness landscape experienced by HIV-1 under treatment selective pressure as a function of the genotypic sequence thereby also estimating the genetic barrier to resistance. Methods We evaluated the performance of two estimated fitness landscapes (nelfinavir [NFV] and zidovudine [AZT] plus lamivudine [3TC]) to predict week 12 viral load (VL) change for 176 treatment change episodes (TCEs) and probability of week 48 virological failure for 90 TCEs, in treatment experienced patients starting these drugs in combination. Results A higher genetic barrier for AZT plus 3TC, (quantified per additional mutation required to develop resistance against these drugs) was associated with a 0.54 (95% confidence interval [CI] 0.30–0.77) larger log10 VL reduction at 12 weeks ( P<0.0001) and a 0.39 (95% CI 0.23–0.66) lower odds of virological failure at 48 weeks ( P=0.0005), in analyses adjusting for the pre-TCE VL and the exact time-lag between the TCE and the date of determining response VL. The strength of these associations was comparable with those seen with expert interpretation systems (Rega, ANRS and HIVDB). A higher genetic barrier to NFV resistance was the only genotypic predictor that tended to be associated with a 0.19 (95% CI 0–0.39) higher log10 VL reduction at 12 weeks ( P=0.05) and a 0.63 (95% CI 0.36–1.09) lower odds of virological failure at 48 weeks ( P=0.10) per additional mutation. Conclusions These results suggest that an estimated genetic barrier derived from fitness landscapes may contribute to an improvement of predicted treatment outcome for NFV and this approach should be explored for other drugs.
Collapse
Affiliation(s)
| | - Koen Deforche
- Rega Institute, Katholieke Universiteit Leuven, Leuven, Belgium
| | | | - Kristof Theys
- Rega Institute, Katholieke Universiteit Leuven, Leuven, Belgium
| | - Bonaventura Clotet
- Irsicaixa Foundation, Hospital Universitari Germans Trias i Pujol, Badalona, Catalonia, Spain
| | | | - Jesper Kjaer
- Copenhagen HIV Programme, University of Copenhagen, Copenhagen, Denmark
| | | | - Andrew Phillips
- Royal Free and University College Medical School, University College London, London, UK
| | - Yves Moreau
- Rega Institute, Katholieke Universiteit Leuven, Leuven, Belgium
| | - Jens D Lundgren
- Copenhagen HIV Programme, University of Copenhagen, Copenhagen, Denmark
- Centre for Viral Diseases (CVD)/KMA, Rigshospitalet, Copenhagen, Denmark
| | | |
Collapse
|
23
|
Gonzalez LMF, Santos AF, Abecasis AB, Van Laethem K, Soares EA, Deforche K, Tanuri A, Camacho R, Vandamme AM, Soares MA. Impact of HIV-1 protease mutations A71V/T and T74S on M89I/V-mediated protease inhibitor resistance in subtype G isolates. J Antimicrob Chemother 2008; 61:1201-4. [DOI: 10.1093/jac/dkn099] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
|
24
|
Vercauteren J, Derdelinckx I, Sasse A, Bogaert M, Ceunen H, De Roo A, De Wit S, Deforche K, Echahidi F, Fransen K, Goffard JC, Goubau P, Goudeseune E, Yombi JC, Lacor P, Liesnard C, Moutschen M, Pierard D, Rens R, Schrooten Y, Vaira D, van den Heuvel A, van der Gucht B, van Ranst M, van Wijngaerden E, Vandercam B, Vekemans M, Verhofstede C, Clumeck N, Vandamme AM, van Laethem K. Prevalence and epidemiology of HIV type 1 drug resistance among newly diagnosed therapy-naive patients in Belgium from 2003 to 2006. AIDS Res Hum Retroviruses 2008; 24:355-62. [PMID: 18327983 DOI: 10.1089/aid.2007.0212] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
This study is the first prospective study to assess the prevalence, epidemiology, and risk factors of HIV-1 drug resistance in newly diagnosed HIV-infected patients in Belgium. In January 2003 it was initiated as part of the pan-European SPREAD program, and continued thereafter for four inclusion rounds until December 2006. Epidemiological, clinical, and behavioral data were collected using a standardized questionnaire and genotypic resistance testing was done on a sample taken within 6 months of diagnosis. Two hundred and eighty-five patients were included. The overall prevalence of transmitted HIV-1 drug resistance in Belgium was 9.5% (27/285, 95% CI: 6.6-13.4). Being infected in Belgium, which largely coincided with harboring a subtype B virus, was found to be significantly associated with transmission of drug resistance. The relatively high rate of baseline resistance might jeopardize the success of first line treatment as more than 1 out of 10 (30/285, 10.5%) viruses did not score as fully susceptible to one of the recommended first-line regimens, i.e., zidovudine, lamivudine, and efavirenz. Our results support the implementation of genotypic resistance testing as a standard of care in all treatment-naive patients in Belgium.
Collapse
Affiliation(s)
- Jurgen Vercauteren
- Rega Institute for Medical Research and University Hospitals, Katholieke Universiteit Leuven, Leuven, Belgium
| | - Inge Derdelinckx
- Rega Institute for Medical Research and University Hospitals, Katholieke Universiteit Leuven, Leuven, Belgium
| | - André Sasse
- Scientific Institute of Public Health, Brussels, Belgium
| | | | - Helga Ceunen
- Rega Institute for Medical Research and University Hospitals, Katholieke Universiteit Leuven, Leuven, Belgium
| | - Ann De Roo
- Prince Leopold Institute of Tropical Medicine, Antwerp, Belgium
| | | | - Koen Deforche
- Rega Institute for Medical Research and University Hospitals, Katholieke Universiteit Leuven, Leuven, Belgium
| | | | - Katrien Fransen
- Prince Leopold Institute of Tropical Medicine, Antwerp, Belgium
| | | | - Patrick Goubau
- Cliniques Universitaires Saint-Luc, Université Catholique de Louvain, Brussels, Belgium
| | | | - Jean-Cyr Yombi
- Cliniques Universitaires Saint-Luc, Université Catholique de Louvain, Brussels, Belgium
| | | | - Corinne Liesnard
- Cliniques Universitaires de Bruxelles, Hôpital Erasme, Brussels, Belgium
| | - Michel Moutschen
- Liège AIDS Reference Center, University of Liège, Liège, Belgium
| | | | - Roeland Rens
- Rega Institute for Medical Research and University Hospitals, Katholieke Universiteit Leuven, Leuven, Belgium
| | - Yoeri Schrooten
- Rega Institute for Medical Research and University Hospitals, Katholieke Universiteit Leuven, Leuven, Belgium
| | - Dolores Vaira
- Liège AIDS Reference Center, University of Liège, Liège, Belgium
| | | | | | - Marc van Ranst
- Rega Institute for Medical Research and University Hospitals, Katholieke Universiteit Leuven, Leuven, Belgium
| | - Eric van Wijngaerden
- Rega Institute for Medical Research and University Hospitals, Katholieke Universiteit Leuven, Leuven, Belgium
| | - Bernard Vandercam
- Cliniques Universitaires Saint-Luc, Université Catholique de Louvain, Brussels, Belgium
| | - Marc Vekemans
- Prince Leopold Institute of Tropical Medicine, Antwerp, Belgium
| | | | | | - Anne-Mieke Vandamme
- Rega Institute for Medical Research and University Hospitals, Katholieke Universiteit Leuven, Leuven, Belgium
| | - Kristel van Laethem
- Rega Institute for Medical Research and University Hospitals, Katholieke Universiteit Leuven, Leuven, Belgium
| |
Collapse
|
25
|
Vercauteren J, Deforche K, Theys K, Debruyne M, Duque LM, Peres S, Carvalho AP, Mansinho K, Vandamme AM, Camacho R. The incidence of multidrug and full class resistance in HIV-1 infected patients is decreasing over time (2001-2006) in Portugal. Retrovirology 2008; 5:12. [PMID: 18241328 PMCID: PMC2265747 DOI: 10.1186/1742-4690-5-12] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2007] [Accepted: 02/01/2008] [Indexed: 11/20/2022] Open
Abstract
Despite improvements in HIV treatment, the prevalence of multidrug resistance and full class resistance is still reported to be increasing. However, to investigate whether current treatment strategies are still selecting for multidrug and full class resistance, the incidence, instead of the prevalence, is more informative. Temporal trends in multidrug resistance (MDR defined as at most 1 drug fully susceptible) and full class resistance (FCR defined as no drug in this class fully susceptible) in Portugal based on 3394 viral isolates genotyped from 2000 to 2006 were examined using the Rega 6.4.1 interpretation system. From July 2001 to July 2006 there was a significant decreasing trend of MDR with 5.7%, 5.2%, 3.8%, 3.4% and 2.7% for the consecutive years (P = 0.003). Multivariate analysis showed that for every consecutive year the odds of having a new MDR case decreased with 20% (P = 0.003). Furthermore, a decline was observed for NRTI- and PI-FCR (both P < 0.001), whereas for NNRTI-FCR a parabolic trend over time was seen (P < 0.001), with a maximum incidence in 2003–'04. Similar trends were obtained when scoring resistance for only one drug within a class or by using another interpretation system. In conclusion, the incidence of multidrug and full class resistance is decreasing over time in Portugal, with the exception of NNRTI full class resistance which showed an initial rise, but subsequently also a decline. This is most probably reflecting the changing drug prescription, the increasing efficiency of HAART and the improved management of HIV drug resistance. This work was presented in part at the Eighth International Congress on Drug Therapy in HIV Infection, Glasgow (UK), 12-16 November 2006 (PL5.5); and at the Fifth European HIV Drug Resistance Workshop, Cascais (Portugal), 28-30 March 2007 (Abstract 1).
Collapse
|
26
|
Ho YS, Abecasis AB, Theys K, Deforche K, Dwyer DE, Charleston M, Vandamme AM, Saksena NK. HIV-1 gp120 N-linked glycosylation differs between plasma and leukocyte compartments. Virol J 2008; 5:14. [PMID: 18215327 PMCID: PMC2265691 DOI: 10.1186/1743-422x-5-14] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2007] [Accepted: 01/23/2008] [Indexed: 02/03/2023] Open
Abstract
Background N-linked glycosylation is a major mechanism for minimizing virus neutralizing antibody response and is present on the Human Immunodeficiency Virus (HIV) envelope glycoprotein. Although it is known that glycosylation changes can dramatically influence virus recognition by the host antibody, the actual contribution of compartmental differences in N-linked glycosylation patterns remains unclear. Methodology and Principal Findings We amplified the env gp120 C2-V5 region and analyzed 305 clones derived from plasma and other compartments from 15 HIV-1 patients. Bioinformatics and Bayesian network analyses were used to examine N-linked glycosylation differences between compartments. We found evidence for cellspecific single amino acid changes particular to monocytes, and significant variation was found in the total number of N-linked glycosylation sites between patients. Further, significant differences in the number of glycosylation sites were observed between plasma and cellular compartments. Bayesian network analyses showed an interdependency between N-linked glycosylation sites found in our study, which may have immense functional relevance. Conclusion Our analyses have identified single cell/compartment-specific amino acid changes and differences in N-linked glycosylation patterns between plasma and diverse blood leukocytes. Bayesian network analyses showed associations inferring alternative glycosylation pathways. We believe that these studies will provide crucial insights into the host immune response and its ability in controlling HIV replication in vivo. These findings could also have relevance in shielding and evasion of HIV-1 from neutralizing antibodies.
Collapse
Affiliation(s)
- Yung Shwen Ho
- Retroviral Genetics Laboratory, Centre for Virus Research, Westmead Millennium Institute, Westmead Hospital, University of Sydney, Westmead NSW, 2145 Sydney, Australia.
| | | | | | | | | | | | | | | |
Collapse
|
27
|
Deforche K, Cozzi-Lepri A, Theys K, Clotet B, Camacho RJ, Kjaer J, Van Laethem K, Phillips A, Moreau Y, Lundgren JD, Vandamme AM. Modelled in vivo HIV fitness under drug selective pressure and estimated genetic barrier towards resistance are predictive for virological response. Antivir Ther 2008; 13:399-407. [PMID: 18572753] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
BACKGROUND A method has been developed to estimate a fitness landscape experienced by HIV-1 under treatment selective pressure as a function of the genotypic sequence thereby also estimating the genetic barrier to resistance. METHODS We evaluated the performance of two estimated fitness landscapes (nelfinavir [NFV] and zidovudine [AZT] plus lamivudine [3TC]) to predict week 12 viral load (VL) change for 176 treatment change episodes (TCEs) and probability of week 48 virological failure for 90 TCEs, in treatment experienced patients starting these drugs in combination. RESULTS A higher genetic barrier for AZT plus 3TC, (quantified per additional mutation required to develop resistance against these drugs) was associated with a 0.54 (95% confidence interval [CI] 0.30-0.77) larger log10 VL reduction at 12 weeks (P < 0.0001) and a 0.39 (95%/ CI 0.23-0.66) lower odds of virological failure at 48 weeks (P = 0.0005), in analyses adjusting for the pre-TCE VL and the exact time-lag between the TCE and the date of determining response VL. The strength of these associations was comparable with those seen with expert interpretation systems (Rega, ANRS and HIVDB). A higher genetic barrier to NFV resistance was the only genotypic predictor that tended to be associated with a 0.19 (95% CI 0-0.39) higher log10 VL reduction at 12 weeks (P = 0.05) and a 0.63 (95% CI 0.36-1.09) lower odds of virological failure at 48 weeks ( P = 0.10) per additional mutation. CONCLUSIONS These results suggest that an estimated genetic barrier derived from fitness landscapes may contribute to an improvement of predicted treatment outcome for NFV and this approach should be explored for other drugs.
Collapse
Affiliation(s)
- Koen Deforche
- Rega Institute, Katholieke Universiteit Leuven, Leuven, Belgium
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
28
|
Deforche K, Camacho R, Laethem KV, Shapiro B, Moreau Y, Rambaut A, Vandamme AM, Lemey P. Estimating the relative contribution of dNTP pool imbalance and APOBEC3G/3F editing to HIV evolution in vivo. J Comput Biol 2007; 14:1105-14. [PMID: 17985990 DOI: 10.1089/cmb.2007.0073] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
The human immunodeficiency virus (HIV) has a genome that is rich in adenine, and its rapid evolution shows an observed bias of guanine (G) to adenine (A) mutations. Two mechanisms have been proposed to explain these properties: (1) an imbalance in dNTP pool concentrations which drives the misincorporation process during reverse transcription, and (2) cytidine deamination by the APOBEC3G/3F restriction factor, causing G to A mutations most notably in specific dinucleotide contexts. Although crucial to understanding HIV evolution, current estimates on misincorporation bias during the replication cycle are based on scarce in vitro measurements. In this work, HIV partial pol sequences obtained for drug resistance testing purposes are analyzed using likelihood methods to estimate various models of HIV misincorporation bias in vivo. The technique is robust to selection on the amino acid sequence and selection against CpG dinucleotides. A model where misincorporations are explained only by an imbalance in dNTP pool concentrations, together with a preference for transitions versus transversions, explained 98% (95% confidence interval [C.I.] 93-100) of the observed variation in freely estimated misincorporation rates. Although dinucleotide context was responsible for variation in misincorporation probabilities, this variation was not specific for G to A mutations implying that the footprint of APOBEC3G/3F editing could not be detected. These results indicate that an imbalance in dNTP pool concentrations explains most of the bias in HIV nucleotide misincorporations, while the effect of editing by APOBEC3G/3F on HIV evolution, based on its dinucleotide specificity, could not be observed in this study.
Collapse
Affiliation(s)
- Koen Deforche
- Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium.
| | | | | | | | | | | | | | | |
Collapse
|
29
|
Deforche K, Camacho R, Van Laethem K, Lemey P, Rambaut A, Moreau Y, Vandamme AM. Estimation of an in vivo fitness landscape experienced by HIV-1 under drug selective pressure useful for prediction of drug resistance evolution during treatment. ACTA ACUST UNITED AC 2007; 24:34-41. [PMID: 18024973 DOI: 10.1093/bioinformatics/btm540] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Abstract
MOTIVATION HIV-1 antiviral resistance is a major cause of antiviral treatment failure. The in vivo fitness landscape experienced by the virus in presence of treatment could in principle be used to determine both the susceptibility of the virus to the treatment and the genetic barrier to resistance. We propose a method to estimate this fitness landscape from cross-sectional clinical genetic sequence data of different subtypes, by reverse engineering the required selective pressure for HIV-1 sequences obtained from treatment naive patients, to evolve towards sequences obtained from treated patients. The method was evaluated for recovering 10 random fictive selective pressures in simulation experiments, and for modeling the selective pressure under treatment with the protease inhibitor nelfinavir. RESULTS The estimated fitness function under nelfinavir treatment considered fitness contributions of 114 mutations at 48 sites. Estimated fitness correlated significantly with the in vitro resistance phenotype in 519 matched genotype-phenotype pairs (R(2) = 0.47 (0.41 - 0.54)) and variation in predicted evolution under nelfinavir selective pressure correlated significantly with observed in vivo evolution during nelfinavir treatment for 39 mutations (with FDR = 0.05). AVAILABILITY The software is available on request from the authors, and data sets are available from http://jose.med.kuleuven.be/~kdforc0/nfv-fitness-data/.
Collapse
Affiliation(s)
- K Deforche
- Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium
| | | | | | | | | | | | | |
Collapse
|
30
|
Deforche K, Camacho R, Grossman Z, Silander T, Soares MA, Moreau Y, Shafer RW, Van Laethem K, Carvalho AP, Wynhoven B, Cane P, Snoeck J, Clarke J, Sirivichayakul S, Ariyoshi K, Holguin A, Rudich H, Rodrigues R, Bouzas MB, Cahn P, Brigido LF, Soriano V, Sugiura W, Phanuphak P, Morris L, Weber J, Pillay D, Tanuri A, Harrigan PR, Shapiro JM, Katzenstein DA, Kantor R, Vandamme AM. Bayesian network analysis of resistance pathways against HIV-1 protease inhibitors. Infect Genet Evol 2006; 7:382-90. [PMID: 17127103 DOI: 10.1016/j.meegid.2006.09.004] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/17/2006] [Revised: 09/08/2006] [Accepted: 09/11/2006] [Indexed: 11/23/2022]
Abstract
Interpretation of Human Immunodeficiency Virus 1 (HIV-1) genotypic drug resistance is still a major challenge in the follow-up of antiviral therapy in infected patients. Because of the high degree of HIV-1 natural variation, complex interactions and stochastic behaviour of evolution, the role of resistance mutations is in many cases not well understood. Using Bayesian network learning of HIV-1 sequence data from diverse subtypes (A, B, C, F and G), we could determine the specific role of many resistance mutations against the protease inhibitors (PIs) nelfinavir (NFV), indinavir (IDV), and saquinavir (SQV). Such networks visualize relationships between treatment, selection of resistance mutations and presence of polymorphisms in a graphical way. The analysis identified 30N, 88S, and 90M for nelfinavir, 90M for saquinavir, and 82A/T and 46I/L for indinavir as most probable major resistance mutations. Moreover we found striking similarities for the role of many mutations against all of these drugs. For example, for all three inhibitors, we found that the novel mutation 89I was minor and associated with mutations at positions 90 and 71. Bayesian network learning provides an autonomous method to gain insight in the role of resistance mutations and the influence of HIV-1 natural variation. We successfully applied the method to three protease inhibitors. The analysis shows differences with current knowledge especially concerning resistance development in several non-B subtypes.
Collapse
Affiliation(s)
- K Deforche
- Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium.
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
31
|
Deforche K, Silander T, Camacho R, Grossman Z, Soares MA, Van Laethem K, Kantor R, Moreau Y, Vandamme AM. Analysis of HIV-1 pol sequences using Bayesian Networks: implications for drug resistance. ACTA ACUST UNITED AC 2006; 22:2975-9. [PMID: 17021157 DOI: 10.1093/bioinformatics/btl508] [Citation(s) in RCA: 52] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
Abstract
Human Immunodeficiency Virus-1 (HIV-1) antiviral resistance is a major cause of antiviral therapy failure and compromises future treatment options. As a consequence, resistance testing is the standard of care. Because of the high degree of HIV-1 natural variation and complex interactions, the role of resistance mutations is in many cases insufficiently understood. We applied a probabilistic model, Bayesian networks, to analyze direct influences between protein residues and exposure to treatment in clinical HIV-1 protease sequences from diverse subtypes. We can determine the specific role of many resistance mutations against the protease inhibitor nelfinavir, and determine relationships between resistance mutations and polymorphisms. We can show for example that in addition to the well-known major mutations 90M and 30N for nelfinavir resistance, 88S should not be treated as 88D but instead considered as a major mutation and explain the subtype-dependent prevalence of the 30N resistance pathway.
Collapse
Affiliation(s)
- K Deforche
- Rega Institute for Medical Research, Katholieke Universiteit Leuven Leuven, Belgium.
| | | | | | | | | | | | | | | | | |
Collapse
|
32
|
Abecasis AB, Deforche K, Bacheler LT, McKenna P, Carvalho AP, Gomes P, Vandamme AM, Camacho RJ. Investigation of Baseline Susceptibility to Protease Inhibitors in HIV-1 Subtypes C, F, G and Crf02_Ag. Antivir Ther 2006. [DOI: 10.1177/135965350601100512] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
Objective To compare baseline susceptibility to protease inhibitors among HIV-1 isolates of subtypes C, F, G and CRF02_AG, and to identify polymorphisms that determine the differences in susceptibility. Methods A total of 42 samples of drug-naive patients infected with subtypes G ( n=19), CRF02_AG ( n=10), F ( n=6) and C ( n=7) were phenotyped and genotyped with the Antivirogram and the ViroSeq 2.0 genotyping system, respectively. A Bayesian network approach was used for a preliminary analysis of the collected data and the dependencies indicated by the network were statistically confirmed. Results CRF02_AG samples were found to be more susceptible to nelfinavir and ritonavir than other subtypes. Hypersusceptibility to these drugs was associated with the 70R polymorphism. 37D/S/T was associated with reduced susceptibility to indinavir and 89M with reduced susceptibility to lopinavir. Susceptibility to tipranavir was the lowest among the subtype F samples and the highest for subtype G samples, with samples carrying 57R being more susceptible than samples carrying 57K. Conclusions Our study suggests that there are baseline susceptibility differences between subtypes and these differences are due to naturally occurring polymorphisms in these subtypes. The predictive value for phenotype of these polymorphisms was even valid in subtypes where these polymorphisms are less prevalent. Taking into account such polymorphisms should improve current algorithms for interpretation of genotyping results in a subtype-independent way.
Collapse
Affiliation(s)
- Ana B Abecasis
- Virology Laboratory, Hospital Egas Moniz, Lisbon, Portugal
- Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium
| | - Koen Deforche
- Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium
| | | | | | | | - Perpétua Gomes
- Virology Laboratory, Hospital Egas Moniz, Lisbon, Portugal
| | - Anne-Mieke Vandamme
- Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium
| | | |
Collapse
|
33
|
Snoeck J, Kantor R, Shafer RW, Van Laethem K, Deforche K, Carvalho AP, Wynhoven B, Soares MA, Cane P, Clarke J, Pillay C, Sirivichayakul S, Ariyoshi K, Holguin A, Rudich H, Rodrigues R, Bouzas MB, Brun-Vézinet F, Reid C, Cahn P, Brigido LF, Grossman Z, Soriano V, Sugiura W, Phanuphak P, Morris L, Weber J, Pillay D, Tanuri A, Harrigan RP, Camacho R, Schapiro JM, Katzenstein D, Vandamme AM. Discordances between interpretation algorithms for genotypic resistance to protease and reverse transcriptase inhibitors of human immunodeficiency virus are subtype dependent. Antimicrob Agents Chemother 2006; 50:694-701. [PMID: 16436728 PMCID: PMC1366873 DOI: 10.1128/aac.50.2.694-701.2006] [Citation(s) in RCA: 66] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
The major limitation of drug resistance genotyping for human immunodeficiency virus remains the interpretation of the results. We evaluated the concordance in predicting therapy response between four different interpretation algorithms (Rega 6.3, HIVDB-08/04, ANRS [07/04], and VGI 8.0). Sequences were gathered through a worldwide effort to establish a database of non-B subtype sequences, and demographic and clinical information about the patients was gathered. The most concordant results were found for nonnucleoside reverse transcriptase (RT) inhibitors (93%), followed by protease inhibitors (84%) and nucleoside RT inhibitor (NRTIs) (76%). For therapy-naive patients, for nelfinavir, especially for subtypes C and G, the discordances were driven mainly by the protease (PRO) mutational pattern 82I/V + 63P + 36I/V for subtype C and 82I + 63P + 36I + 20I for subtype G. Subtype F displayed more discordances for ritonavir in untreated patients due to the combined presence of PRO 20R and 10I/V. In therapy-experienced patients, subtype G displayed a lot of discordances for saquinavir and indinavir due to mutational patterns involving PRO 90 M and 82I. Subtype F had more discordance for nelfinavir attributable to the presence of PRO 88S and 82A + 54V. For the NRTIs lamivudine and emtricitabine, CRF01_AE had more discordances than subtype B due to the presence of RT mutational patterns 65R + 115 M and 118I + 215Y, respectively. Overall, the different algorithms agreed well on the level of resistance scored, but some of the discordances could be attributed to specific (subtype-dependent) combinations of mutations. It is not yet known whether therapy response is subtype dependent, but the advice given to clinicians based on a genotypic interpretation algorithm differs according to the subtype.
Collapse
Affiliation(s)
- Joke Snoeck
- Rega Institute for Medical Research, Minderbroedersstraat 10, 3000 Leuven, Belgium
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
34
|
Abecasis AB, Deforche K, Bacheler LT, McKenna P, Carvalho AP, Gomes P, Vandamme AM, Camacho RJ. Investigation of baseline susceptibility to protease inhibitors in HIV-1 subtypes C, F, G and CRF02_AG. Antivir Ther 2006; 11:581-9. [PMID: 16964826] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
OBJECTIVE To compare baseline susceptibility to protease inhibitors among HIV-1 isolates of subtypes C, F, G and CRF02_AG, and to identify polymorphisms that determine the differences in susceptibility. METHODS A total of 42 samples of drug-naive patients infected with subtypes G (n=19), CRF02_AG (n = 10), F (n = 6) and C (n = 7) were phenotyped and genotyped with the Antivirogram and the ViroSeq 2.0 genotyping system, respectively. A Bayesian network approach was used for a preliminary analysis of the collected data and the dependencies indicated by the network were statistically confirmed. RESULTS CRF02_AG samples were found to be more susceptible to nelfinavir and ritonavir than other subtypes. Hypersusceptibility to these drugs was associated with the 70R polymorphism. 37D/S/T was associated with reduced susceptibility to indinavir and 89M with reduced susceptibility to lopinavir. Susceptibility to tipranavir was the lowest among the subtype F samples and the highest for subtype G samples, with samples carrying 57R being more susceptible than samples carrying 57K. CONCLUSIONS Our study suggests that there are baseline susceptibility differences between subtypes and these differences are due to naturally occurring polymorphisms in these subtypes. The predictive value for phenotype of these polymorphisms was even valid in subtypes where these polymorphisms are less prevalent. Taking into account such polymorphisms should improve current algorithms for interpretation of genotyping results in a subtype-independent way.
Collapse
Affiliation(s)
- Ana B Abecasis
- Virology Laboratory, Hospital Egas Moniz, Lisbon, Portugal.
| | | | | | | | | | | | | | | |
Collapse
|
35
|
Abecasis AB, Deforche K, Snoeck J, Bacheler LT, McKenna P, Carvalho AP, Gomes P, Camacho RJ, Vandamme AM. Protease mutation M89I/V is linked to therapy failure in patients infected with the HIV-1 non-B subtypes C, F or G. AIDS 2005; 19:1799-806. [PMID: 16227787 DOI: 10.1097/01.aids.0000188422.95162.b7] [Citation(s) in RCA: 57] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
OBJECTIVE To investigate whether and how mutations at position 89 of HIV-1 protease were associated with protease inhibitor (PI) failure, and what is the impact of the HIV-1 subtype. METHODS In a database containing pol nucleotide sequences and treatment history, the correlation between PI experience and mutations at codon 89 was determined separately for subtype B and several non-B subtypes. A Bayesian network model was used to map the resistance pathways in which M89I/V is involved for subtype G. The phenotypic effect of M89I/V for several PIs was also measured. RESULTS The analysis showed that for the subtypes C, F and G in which the wild-type codon at 89 was M compared to L for subtype B, M89I/V was significantly more frequently observed in PI-treated patients displaying major resistance mutations to PIs than in drug-naive patients. M89I/V was strongly associated with PI resistance mutations at codons 71, 74 and 90. Phenotypically, M89I/V alone did not confer a reduced susceptibility to PIs. However, when combined with L90M, a significantly reduced susceptibility to nelfinavir was observed (P < 0.05) in comparison with strains with L90M alone. CONCLUSIONS The results of the present study show that M89I/V is associated with PI experience in subtypes C, F and G but not in subtype B. M89I/V should be considered a secondary PI mutation with an important effect on nelfinavir susceptibility in the presence of L90M.
Collapse
Affiliation(s)
- Ana Barroso Abecasis
- Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium.
| | | | | | | | | | | | | | | | | |
Collapse
|
36
|
Van Laethem K, Schrooten Y, Dedecker S, Van Heeswijck L, Deforche K, Van Wijngaerden E, Van Ranst M, Vandamme AM. A genotypic assay for the amplification and sequencing of gag and protease from diverse human immunodeficiency virus type 1 group M subtypes. J Virol Methods 2005; 132:181-6. [PMID: 16271771 DOI: 10.1016/j.jviromet.2005.10.008] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2005] [Revised: 09/13/2005] [Accepted: 10/03/2005] [Indexed: 10/25/2022]
Abstract
In human immunodeficiency virus type 1 (HIV-1), an interaction exists between the in vivo evolution of Gag protein and protease to escape from antiretroviral drug selective pressure. Therefore, it was decided to develop a genotypic assay for the amplification and sequencing of HIV-1 gag and protease. As the HIV-1 pandemic is characterised by a large genetic diversity, the assay developed was evaluated on a panel of 28 genetically divergent samples belonging to the following subtypes A1, B, C, D, F1, F2, G, H, J, CRF01-AE, CRF02-AG and CRF13-cpx. The assay displayed a detection limit ranging between 500 RNA copies/ml and 5000 RNA copies/ml plasma. Full-length sequences could be obtained for 25 samples. The population sequences of the three other samples lacked a part of the sequence because of heterogeneous signal, probably due to the presence of quasi-species with insertions/deletions of a different length.
Collapse
Affiliation(s)
- Kristel Van Laethem
- Rega Institute for Medical Research, Katholieke Universiteit Leuven, Leuven, Belgium.
| | | | | | | | | | | | | | | |
Collapse
|
37
|
de Oliveira T, Deforche K, Cassol S, Salminen M, Paraskevis D, Seebregts C, Snoeck J, van Rensburg EJ, Wensing AMJ, van de Vijver DA, Boucher CA, Camacho R, Vandamme AM. An automated genotyping system for analysis of HIV-1 and other microbial sequences. Bioinformatics 2005; 21:3797-800. [PMID: 16076886 DOI: 10.1093/bioinformatics/bti607] [Citation(s) in RCA: 402] [Impact Index Per Article: 21.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
MOTIVATION Genetic analysis of HIV-1 is important not only for vaccine development, but also to guide treatment strategies, track the emergence of new viral variants and ensure that diagnostic assays are contemporary and fully optimized. However, most genotyping methods are laborious and complex, and involve the use of multiple software applications. Here, we describe the development of an automated genotyping system that can be easily applied to HIV-1 and other rapidly evolving viral pathogens. RESULTS The new REGA subtyping tool, developed using Java programming and PERL scripts, combines phylogenetic analyses with boot-scanning methods for the genetic subtyping of full-length and subgenomic fragments of HIV-1. When used to investigate the subtype of previously published reference datasets that were analysed using manual phylogenetic methods, the automated method correctly identified 97.5-100% of non-recombinant and circulating recombinant forms of HIV-1, including 108 full-length, 108 gag and 221 env sequences downloaded from the Los Alamos database.
Collapse
Affiliation(s)
- Tulio de Oliveira
- Evolution Group at the Zoology Department, University of Oxford, UK.
| | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
38
|
Paraskevis D, Deforche K, Abecasis A, Camacho R, Vandamme AM. Analysis of complex HIV-1 intersubtype recombinants using a Bayesian scanning method. Infect Genet Evol 2005; 5:219-24. [PMID: 15737912 DOI: 10.1016/j.meegid.2004.07.006] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/14/2004] [Accepted: 07/01/2004] [Indexed: 11/26/2022]
Abstract
The increased complexity of HIV-1 genetic heterogeneity raises the issue for reliable classification and analysis of these sequences. Until now, bootscanning analysis has been the main method used for the analysis of potential HIV-1 intersubtype recombinants. We show evidence that in some cases of complex recombinants, where three or more segments with discordant phylogenetic signal may exist in protease (PR) and partial reverse transcriptase (RT) region, Bayesian scanning provides a clearer picture than bootscanning plots about the boundaries of potential recombination. Thus, a recently developed Bayesian scanning tool can facilitate the analysis and classification of HIV-1 mosaic sequences.
Collapse
Affiliation(s)
- D Paraskevis
- Laboratory for Clinical and Epidemiological Virology, Rega Institute for Medical Research, Katholieke Universiteit Leuven, Minderbroedersstraat 10, B-3000 Leuven, Belgium.
| | | | | | | | | |
Collapse
|
39
|
Paraskevis D, Deforche K, Lemey P, Magiorkinis G, Hatzakis A, Vandamme AM. SlidingBayes: exploring recombination using a sliding window approach based on Bayesian phylogenetic inference. Bioinformatics 2004; 21:1274-5. [PMID: 15546940 DOI: 10.1093/bioinformatics/bti139] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
We developed a software tool (SlidingBayes) for recombination analysis based on Bayesian phylogenetic inference. Sliding-Bayes provides a powerful approach for detecting potential recombination, especially between highly divergent sequences and complex HIV-1 recombinants for which simpler methods like neighbor joining (NJ) may be less powerful. SlidingBayes guides Markov Chain Monte Carlo (MCMC) sampling performed by MrBayes in a sliding window across the alignment (Bayesian scanning). The tool can be used for nucleotide and amino acid sequences and combines all the modeling possibilities of MrBayes with the ability to plot the posterior probability support for clustering of various combinations of taxa.
Collapse
Affiliation(s)
- D Paraskevis
- Laboratory for Clinical and Epidemiological Virology, Rega Institute for Medical Research, Katholieke Universiteit Leuven, Minderbroedersstraat 10, B-3000 Leuven, Belgium.
| | | | | | | | | | | |
Collapse
|