1
|
Iglhaut C, Pečerska J, Gil M, Anisimova M. Please Mind the Gap: Indel-Aware Parsimony for Fast and Accurate Ancestral Sequence Reconstruction and Multiple Sequence Alignment Including Long Indels. Mol Biol Evol 2024; 41:msae109. [PMID: 38842253 PMCID: PMC11221656 DOI: 10.1093/molbev/msae109] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2024] [Revised: 05/30/2024] [Accepted: 06/03/2024] [Indexed: 06/07/2024] Open
Abstract
Despite having important biological implications, insertion, and deletion (indel) events are often disregarded or mishandled during phylogenetic inference. In multiple sequence alignment, indels are represented as gaps and are estimated without considering the distinct evolutionary history of insertions and deletions. Consequently, indels are usually excluded from subsequent inference steps, such as ancestral sequence reconstruction and phylogenetic tree search. Here, we introduce indel-aware parsimony (indelMaP), a novel way to treat gaps under the parsimony criterion by considering insertions and deletions as separate evolutionary events and accounting for long indels. By identifying the precise location of an evolutionary event on the tree, we can separate overlapping indel events and use affine gap penalties for long indel modeling. Our indel-aware approach harnesses the phylogenetic signal from indels, including them into all inference stages. Validation and comparison to state-of-the-art inference tools on simulated data show that indelMaP is most suitable for densely sampled datasets with closely to moderately related sequences, where it can reach alignment quality comparable to probabilistic methods and accurately infer ancestral sequences, including indel patterns. Due to its remarkable speed, our method is well suited for epidemiological datasets, eliminating the need for downsampling and enabling the exploitation of the additional information provided by dense taxonomic sampling. Moreover, indelMaP offers new insights into the indel patterns of biologically significant sequences and advances our understanding of genetic variability by considering gaps as crucial evolutionary signals rather than mere artefacts.
Collapse
Affiliation(s)
- Clara Iglhaut
- Institute of Computational Life Science, Zurich University of Applied Science, Wädenswil, Switzerland
- Faculty of Mathematics and Science, University of Zurich, Zürich, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Jūlija Pečerska
- Institute of Computational Life Science, Zurich University of Applied Science, Wädenswil, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Manuel Gil
- Institute of Computational Life Science, Zurich University of Applied Science, Wädenswil, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Maria Anisimova
- Institute of Computational Life Science, Zurich University of Applied Science, Wädenswil, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| |
Collapse
|
2
|
Miton CM, Tokuriki N. Insertions and Deletions (Indels): A Missing Piece of the Protein Engineering Jigsaw. Biochemistry 2023; 62:148-157. [PMID: 35830609 DOI: 10.1021/acs.biochem.2c00188] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Over the years, protein engineers have studied nature and borrowed its tricks to accelerate protein evolution in the test tube. While there have been considerable advances, our ability to generate new proteins in the laboratory is seemingly limited. One explanation for these shortcomings may be that insertions and deletions (indels), which frequently arise in nature, are largely overlooked during protein engineering campaigns. The profound effect of indels on protein structures, by way of drastic backbone alterations, could be perceived as "saltation" events that bring about significant phenotypic changes in a single mutational step. Should we leverage these effects to accelerate protein engineering and gain access to unexplored regions of adaptive landscapes? In this Perspective, we describe the role played by indels in the functional diversification of proteins in nature and discuss their untapped potential for protein engineering, despite their often-destabilizing nature. We hope to spark a renewed interest in indels, emphasizing that their wider study and use may prove insightful and shape the future of protein engineering by unlocking unique functional changes that substitutions alone could never achieve.
Collapse
Affiliation(s)
- Charlotte M Miton
- Michael Smith Laboratories, University of British Columbia, Vancouver, V6T 1Z4 BC, Canada
| | - Nobuhiko Tokuriki
- Michael Smith Laboratories, University of British Columbia, Vancouver, V6T 1Z4 BC, Canada
| |
Collapse
|
3
|
Ndashimye E, Reyes PS, Arts EJ. New antiretroviral inhibitors and HIV-1 drug resistance: more focus on 90% HIV-1 isolates? FEMS Microbiol Rev 2023; 47:fuac040. [PMID: 36130204 PMCID: PMC9841967 DOI: 10.1093/femsre/fuac040] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2022] [Revised: 09/13/2022] [Accepted: 09/18/2022] [Indexed: 01/21/2023] Open
Abstract
Combined HIV antiretroviral therapy (cART) has been effective except if drug resistance emerges. As cART has been rolled out in low-income countries, drug resistance has emerged at higher rates than observed in high income countries due to factors including initial use of these less tolerated cART regimens, intermittent disruptions in drug supply, and insufficient treatment monitoring. These socioeconomic factors impacting drug resistance are compounded by viral mechanistic differences by divergent HIV-1 non-B subtypes compared to HIV-1 subtype B that largely infects the high-income countries (just 10% of 37 million infected). This review compares the inhibition and resistance of diverse HIV-1 subtypes and strains to the various approved drugs as well as novel inhibitors in clinical trials. Initial sequence variations and differences in replicative fitness between HIV-1 subtypes pushes strains through different fitness landscapes to escape from drug selective pressure. The discussions here provide insight to patient care givers and policy makers on how best to use currently approved ART options and reduce the emergence of drug resistance in ∼33 million individuals infected with HIV-1 subtype A, C, D, G, and recombinants forms. Unfortunately, over 98% of the literature on cART resistance relates to HIV-1 subtype B.
Collapse
Affiliation(s)
- Emmanuel Ndashimye
- Department of Microbiology and Immunology, Western University Schulich School of Medicine & Dentistry, Western University, N6A 3K7, London, Ontario, Canada
- Joint Clinical Research Centre, -Center for AIDS Research Laboratories, 256, Kampala, Uganda
| | - Paul S Reyes
- Department of Microbiology and Immunology, Western University Schulich School of Medicine & Dentistry, Western University, N6A 3K7, London, Ontario, Canada
| | - Eric J Arts
- Department of Microbiology and Immunology, Western University Schulich School of Medicine & Dentistry, Western University, N6A 3K7, London, Ontario, Canada
| |
Collapse
|
4
|
Alisoltani A, Jaroszewski L, Iyer M, Iranzadeh A, Godzik A. Increased Frequency of Indels in Hypervariable Regions of SARS-CoV-2 Proteins—A Possible Signature of Adaptive Selection. Front Genet 2022; 13:875406. [PMID: 35719386 PMCID: PMC9201826 DOI: 10.3389/fgene.2022.875406] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2022] [Accepted: 04/08/2022] [Indexed: 11/29/2022] Open
Abstract
Most attention in the surveillance of evolving SARS-CoV-2 genome has been centered on nucleotide substitutions in the spike glycoprotein. We show that, as the pandemic extends into its second year, the numbers and ratio of genomes with in-frame insertions and deletions (indels) increases significantly, especially among the variants of concern (VOCs). Monitoring of the SARS-CoV-2 genome evolution shows that co-occurrence (i.e., highly correlated presence) of indels, especially deletions on spike N-terminal domain and non-structural protein 6 (NSP6) is a shared feature in several VOCs such as Alpha, Beta, Delta, and Omicron. Indels distribution is correlated with spike mutations associated with immune escape and growth in the number of genomes with indels coincides with the increasing population resistance due to vaccination and previous infections. Indels occur most frequently in the spike, but also in other proteins, especially those involved in interactions with the host immune system. We also showed that indels concentrate in regions of individual SARS-CoV-2 proteins known as hypervariable regions (HVRs) that are mostly located in specific loop regions. Structural analysis suggests that indels remodel viral proteins’ surfaces at common epitopes and interaction interfaces, affecting the virus’ interactions with host proteins. We hypothesize that the increased frequency of indels, the non-random distribution of them and their independent co-occurrence in several VOCs is another mechanism of response to elevated global population immunity.
Collapse
Affiliation(s)
- Arghavan Alisoltani
- Biosciences Division, School of Medicine, University of California, Riverside, Riverside, CA, United States
| | - Lukasz Jaroszewski
- Biosciences Division, School of Medicine, University of California, Riverside, Riverside, CA, United States
| | - Mallika Iyer
- Graduate School of Biomedical Sciences, Sanford Burnham Prebys Medical Discovery Institute, La Jolla, CA, United States
| | - Arash Iranzadeh
- Computational Biology Division, Department of Integrative Biomedical Sciences, University of Cape Town, Cape Town, South Africa
| | - Adam Godzik
- Biosciences Division, School of Medicine, University of California, Riverside, Riverside, CA, United States
- *Correspondence: Adam Godzik,
| |
Collapse
|
5
|
Large Evolutionary Rate Heterogeneity among and within HIV-1 Subtypes and CRFs. Viruses 2021; 13:v13091689. [PMID: 34578270 PMCID: PMC8473000 DOI: 10.3390/v13091689] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2021] [Revised: 08/19/2021] [Accepted: 08/23/2021] [Indexed: 01/06/2023] Open
Abstract
HIV-1 is a fast-evolving, genetically diverse virus presently classified into several groups and subtypes. The virus evolves rapidly because of an error-prone polymerase, high rates of recombination, and selection in response to the host immune system and clinical management of the infection. The rate of evolution is also influenced by the rate of virus spread in a population and nature of the outbreak, among other factors. HIV-1 evolution is thus driven by a range of complex genetic, social, and epidemiological factors that complicates disease management and prevention. Here, we quantify the evolutionary (substitution) rate heterogeneity among major HIV-1 subtypes and recombinants by analyzing the largest collection of HIV-1 genetic data spanning the widest possible geographical (100 countries) and temporal (1981–2019) spread. We show that HIV-1 substitution rates vary substantially, sometimes by several folds, both across the virus genome and between major subtypes and recombinants, but also within a subtype. Across subtypes, rates ranged 3.5-fold from 1.34 × 10−3 to 4.72 × 10−3 in env and 2.3-fold from 0.95 × 10−3 to 2.18 × 10−3 substitutions site−1 year−1 in pol. Within the subtype, 3-fold rate variation was observed in env in different human populations. It is possible that HIV-1 lineages in different parts of the world are operating under different selection pressures leading to substantial rate heterogeneity within and between subtypes. We further highlight how such rate heterogeneity can complicate HIV-1 phylodynamic studies, specifically, inferences on epidemiological linkage of transmission clusters based on genetic distance or phylogenetic data, and can mislead estimates about the timing of HIV-1 lineages.
Collapse
|
6
|
Kurniasari MD, Huruta AD, Tsai HT, Lee CW. Forecasting future HIV infection cases: evidence from Indonesia. SOCIAL WORK IN PUBLIC HEALTH 2021; 36:12-25. [PMID: 33272146 DOI: 10.1080/19371918.2020.1851332] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Countries throughout the world, including Indonesia, are facing a complex problem with regards to HIV infection incidences and its prevalence. This is despite some local governments in some provinces of Indonesia working together with the Social Ministry of Indonesia to eradicate prostitution. There are high numbers of HIV sub-types in Indonesia such as HIV-1 (CRF01_AE and B). The forecast was conducted with the Autoregressive Integrated Moving Average model. The ARMA (1,1) was observed to be the best model for forecasting the number of HIV patients in Indonesia. This forecasting has been done since March 2019. Based on its dynamic forecasting with ARMA (1,1), this study proved the number of HIV-positive patients, from 2019 to 2030, had increased from 22,679 to 36,255, by almost 37% within 12 years. Indonesia is facing a growing trend in the number of new HIV cases, until 2030 which caused by stopped their follow-up treatments or they have ceased in consuming the Antiretroviral drugs even though the Indonesian government was provided national health insurance which covers the Antiretroviral drug and a limited number of health-care services providing the Antiretroviral therapy. Therefore, investigations focusing on estimate the number of HIV patients in Indonesia is an important finding. The information can be used as a resource for policy and decision making for plans and programs.
Collapse
Affiliation(s)
- Maria Dyah Kurniasari
- School of Nursing, College of Nursing, Taipei Medical University, Taipei City, Republic of Tiongkok
- Department of Nursing, Faculty of Medicine and Health Science, Universitas Kristen Satya Wacana, Salatiga City, Indonesia
| | - Andrian Dolfriandra Huruta
- Program in Business, College of Business, Chung Yuan Christian University, Taoyuan City, Taiwan, ROC
- Department of Economics, Faculty of Economic and Business, Universitas Kristen Satya Wacana, Salatiga City, Indonesia
| | - Hsiu Ting Tsai
- School of Nursing, College of Nursing, Taipei Medical University, Taipei City, Republic of Tiongkok
- School of Nursing, Post-Baccalaureate Program in Nursing, Taipei Medical University, Taipei City, Republic of China
| | - Cheng-Wen Lee
- Department of International Business, College of Business, Chung Yuan Christian University, Taoyuan City, Taiwan, ROC
| |
Collapse
|
7
|
Abstract
Viruses are widely used as vectors for heterologous gene expression in cultured cells or natural hosts, and therefore a large number of viruses with exogenous sequences inserted into their genomes have been engineered. Many of these engineered viruses are viable and express heterologous proteins at high levels, but the inserted sequences often prove to be unstable over time and are rapidly lost, limiting heterologous protein expression. Although virologists are aware that inserted sequences can be unstable, processes leading to insert instability are rarely considered from an evolutionary perspective. Here, we review experimental work on the stability of inserted sequences over a broad range of viruses, and we present some theoretical considerations concerning insert stability. Different virus genome organizations strongly impact insert stability, and factors such as the position of insertion can have a strong effect. In addition, we argue that insert stability not only depends on the characteristics of a particular genome, but that it will also depend on the host environment and the demography of a virus population. The interplay between all factors affecting stability is complex, which makes it challenging to develop a general model to predict the stability of genomic insertions. We highlight key questions and future directions, finding that insert stability is a surprisingly complex problem and that there is need for mechanism-based, predictive models. Combining theoretical models with experimental tests for stability under varying conditions can lead to improved engineering of viral modified genomes, which is a valuable tool for understanding genome evolution as well as for biotechnological applications, such as gene therapy.
Collapse
Affiliation(s)
- Anouk Willemsen
- Laboratory MIVEGEC (UMR CNRS IRD University of Montpellier), Centre National de la Recherche Scientifique (CNRS), 911 Avenue Agropolis, BP 64501, 34394 Montpellier cedex 5, France
| | - Mark P Zwart
- Netherlands Institute of Ecology (NIOO-KNAW), Postbus 50, 6700 AB, Wageningen, The Netherlands
| |
Collapse
|