1
|
Hong Z, Shimagaki KS, Barton JP. popDMS infers mutation effects from deep mutational scanning data. Bioinformatics 2024; 40:btae499. [PMID: 39115383 PMCID: PMC11335369 DOI: 10.1093/bioinformatics/btae499] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2024] [Revised: 07/10/2024] [Accepted: 08/06/2024] [Indexed: 08/22/2024] Open
Abstract
SUMMARY Deep mutational scanning (DMS) experiments provide a powerful method to measure the functional effects of genetic mutations at massive scales. However, the data generated from these experiments can be difficult to analyze, with significant variation between experimental replicates. To overcome this challenge, we developed popDMS, a computational method based on population genetics theory, to infer the functional effects of mutations from DMS data. Through extensive tests, we found that the functional effects of single mutations and epistasis inferred by popDMS are highly consistent across replicates, comparing favorably with existing methods. Our approach is flexible and can be widely applied to DMS data that includes multiple time points, multiple replicates, and different experimental conditions. AVAILABILITY AND IMPLEMENTATION popDMS is implemented in Python and Julia, and is freely available on GitHub at https://github.com/bartonlab/popDMS.
Collapse
Affiliation(s)
- Zhenchen Hong
- Department of Physics and Astronomy, University of California, Riverside, CA 92521, United States
| | - Kai S Shimagaki
- Department of Computational and Systems Biology, University of Pittsburgh School of Medicine, PA 15260, United States
| | - John P Barton
- Department of Physics and Astronomy, University of California, Riverside, CA 92521, United States
- Department of Computational and Systems Biology, University of Pittsburgh School of Medicine, PA 15260, United States
- Department of Physics and Astronomy, University of Pittsburgh, PA 15260, United States
| |
Collapse
|
2
|
Bose D, Deb Adhikary N, Xiao P, Rogers KA, Ferrell DE, Cheng-Mayer C, Chang TL, Villinger F. SHIV-C109p5 NHP induces rapid disease progression in elderly macaques with extensive GI viral replication. J Virol 2024; 98:e0165223. [PMID: 38299866 PMCID: PMC10878093 DOI: 10.1128/jvi.01652-23] [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: 10/23/2023] [Accepted: 01/02/2024] [Indexed: 02/02/2024] Open
Abstract
CCR5-tropic simian/human immunodeficiency viruses (SHIV) with clade C transmitted/founder envelopes represent a critical tool for the investigation of HIV experimental vaccines and microbicides in nonhuman primates, although many such isolates lead to spontaneous viral control post infection. Here, we generated a high-titer stock of pathogenic SHIV-C109p5 by serial passage in two rhesus macaques (RM) and tested its virulence in aged monkeys. The co-receptor usage was confirmed before infecting five geriatric rhesus macaques (four female and one male). Plasma viral loads were monitored by reverse transcriptase-quantitative PCR (RT-qPCR), cytokines by multiplex analysis, and biomarkers of gastrointestinal damage by enzyme-linked immunosorbent assay. Antibodies and cell-mediated responses were also measured. Viral dissemination into tissues was determined by RNAscope. Intravenous SHIV-C109p5 infection of aged RMs leads to high plasma viremia and rapid disease progression; rapid decrease in CD4+ T cells, CD4+CD8+ T cells, and plasmacytoid dendritic cells; and wasting necessitating euthanasia between 3 and 12 weeks post infection. Virus-specific cellular immune responses were detected only in the two monkeys that survived 4 weeks post infection. These were Gag-specific TNFα+CD8+, MIP1β+CD4+, Env-specific IFN-γ+CD4+, and CD107a+ T cell responses. Four out of five monkeys had elevated intestinal fatty acid binding protein levels at the viral peak, while regenerating islet-derived protein 3α showed marked increases at later time points in the three animals surviving the longest, suggesting gut antimicrobial peptide production in response to microbial translocation post infection. Plasma levels of monocyte chemoattractant protein-1, interleukin-15, and interleukin-12/23 were also elevated. Viral replication in gut and secondary lymphoid tissues was extensive.IMPORTANCESimian/human immunodeficiency viruses (SHIV) are important reagents to study prevention of virus acquisition in nonhuman primate models of HIV infection, especially those representing transmitted/founder (T/F) viruses. However, many R5-tropic SHIV have limited fitness in vivo leading to many monkeys spontaneously controlling the virus post acute infection. Here, we report the generation of a pathogenic SHIV clade C T/F stock by in vivo passage leading to sustained viral load set points, a necessity to study pathogenicity. Unexpectedly, administration of this SHIV to elderly rhesus macaques led to extensive viral replication and fast disease progression, despite maintenance of a strict R5 tropism. Such age-dependent rapid disease progression had previously been reported for simian immunodeficiency virus but not for R5-tropic SHIV infections.
Collapse
Affiliation(s)
- Deepanwita Bose
- New Iberia Research Center, University of Louisiana at Lafayette, New Iberia, Louisiana, USA
| | - Nihar Deb Adhikary
- New Iberia Research Center, University of Louisiana at Lafayette, New Iberia, Louisiana, USA
| | - Peng Xiao
- New Iberia Research Center, University of Louisiana at Lafayette, New Iberia, Louisiana, USA
| | - Kenneth A. Rogers
- New Iberia Research Center, University of Louisiana at Lafayette, New Iberia, Louisiana, USA
| | - Douglas E. Ferrell
- New Iberia Research Center, University of Louisiana at Lafayette, New Iberia, Louisiana, USA
| | | | - Theresa L. Chang
- The Public Health Research Institute, Rutgers New Jersey Medical School, Newark, New Jersey, USA
| | - Francois Villinger
- New Iberia Research Center, University of Louisiana at Lafayette, New Iberia, Louisiana, USA
| |
Collapse
|
3
|
Hong Z, Barton JP. popDMS infers mutation effects from deep mutational scanning data. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.01.29.577759. [PMID: 38352383 PMCID: PMC10862717 DOI: 10.1101/2024.01.29.577759] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 02/21/2024]
Abstract
Deep mutational scanning (DMS) experiments provide a powerful method to measure the functional effects of genetic mutations at massive scales. However, the data generated from these experiments can be difficult to analyze, with significant variation between experimental replicates. To overcome this challenge, we developed popDMS, a computational method based on population genetics theory, to infer the functional effects of mutations from DMS data. Through extensive tests, we found that the functional effects of single mutations and epistasis inferred by popDMS are highly consistent across replicates, comparing favorably with existing methods. Our approach is flexible and can be widely applied to DMS data that includes multiple time points, multiple replicates, and different experimental conditions.
Collapse
Affiliation(s)
- Zhenchen Hong
- Department of Physics and Astronomy, University of California, Riverside, USA
| | - John P. Barton
- Department of Physics and Astronomy, University of California, Riverside, USA
- Department of Computational and Systems Biology, University of Pittsburgh School of Medicine, USA
- Department of Physics and Astronomy, University of Pittsburgh, USA
| |
Collapse
|
4
|
Thadani NN, Gurev S, Notin P, Youssef N, Rollins NJ, Ritter D, Sander C, Gal Y, Marks DS. Learning from prepandemic data to forecast viral escape. Nature 2023; 622:818-825. [PMID: 37821700 PMCID: PMC10599991 DOI: 10.1038/s41586-023-06617-0] [Citation(s) in RCA: 54] [Impact Index Per Article: 27.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2022] [Accepted: 09/06/2023] [Indexed: 10/13/2023]
Abstract
Effective pandemic preparedness relies on anticipating viral mutations that are able to evade host immune responses to facilitate vaccine and therapeutic design. However, current strategies for viral evolution prediction are not available early in a pandemic-experimental approaches require host polyclonal antibodies to test against1-16, and existing computational methods draw heavily from current strain prevalence to make reliable predictions of variants of concern17-19. To address this, we developed EVEscape, a generalizable modular framework that combines fitness predictions from a deep learning model of historical sequences with biophysical and structural information. EVEscape quantifies the viral escape potential of mutations at scale and has the advantage of being applicable before surveillance sequencing, experimental scans or three-dimensional structures of antibody complexes are available. We demonstrate that EVEscape, trained on sequences available before 2020, is as accurate as high-throughput experimental scans at anticipating pandemic variation for SARS-CoV-2 and is generalizable to other viruses including influenza, HIV and understudied viruses with pandemic potential such as Lassa and Nipah. We provide continually revised escape scores for all current strains of SARS-CoV-2 and predict probable further mutations to forecast emerging strains as a tool for continuing vaccine development ( evescape.org ).
Collapse
Affiliation(s)
- Nicole N Thadani
- Marks Group, Department of Systems Biology, Harvard Medical School, Boston, MA, USA
| | - Sarah Gurev
- Marks Group, Department of Systems Biology, Harvard Medical School, Boston, MA, USA
- Department of Electrical Engineering and Computer Science, MIT, Cambridge, MA, USA
| | - Pascal Notin
- OATML Group, Department of Computer Science, University of Oxford, Oxford, UK
| | - Noor Youssef
- Marks Group, Department of Systems Biology, Harvard Medical School, Boston, MA, USA
| | - Nathan J Rollins
- Marks Group, Department of Systems Biology, Harvard Medical School, Boston, MA, USA
- Seismic Therapeutic, Watertown, MA, USA
| | - Daniel Ritter
- Marks Group, Department of Systems Biology, Harvard Medical School, Boston, MA, USA
| | - Chris Sander
- Marks Group, Department of Systems Biology, Harvard Medical School, Boston, MA, USA
- Broad Institute of Harvard and MIT, Cambridge, MA, USA
| | - Yarin Gal
- OATML Group, Department of Computer Science, University of Oxford, Oxford, UK
| | - Debora S Marks
- Marks Group, Department of Systems Biology, Harvard Medical School, Boston, MA, USA.
- Broad Institute of Harvard and MIT, Cambridge, MA, USA.
| |
Collapse
|
5
|
Haddox HK, Galloway JG, Dadonaite B, Bloom JD, Matsen IV FA, DeWitt WS. Jointly modeling deep mutational scans identifies shifted mutational effects among SARS-CoV-2 spike homologs. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.07.31.551037. [PMID: 37577604 PMCID: PMC10418112 DOI: 10.1101/2023.07.31.551037] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/15/2023]
Abstract
Deep mutational scanning (DMS) is a high-throughput experimental technique that measures the effects of thousands of mutations to a protein. These experiments can be performed on multiple homologs of a protein or on the same protein selected under multiple conditions. It is often of biological interest to identify mutations with shifted effects across homologs or conditions. However, it is challenging to determine if observed shifts arise from biological signal or experimental noise. Here, we describe a method for jointly inferring mutational effects across multiple DMS experiments while also identifying mutations that have shifted in their effects among experiments. A key aspect of our method is to regularize the inferred shifts, so that they are nonzero only when strongly supported by the data. We apply this method to DMS experiments that measure how mutations to spike proteins from SARS-CoV-2 variants (Delta, Omicron BA.1, and Omicron BA.2) affect cell entry. Most mutational effects are conserved between these spike homologs, but a fraction have markedly shifted. We experimentally validate a subset of the mutations inferred to have shifted effects, and confirm differences of > 1,000-fold in the impact of the same mutation on spike-mediated viral infection across spikes from different SARS-CoV-2 variants. Overall, our work establishes a general approach for comparing sets of DMS experiments to identify biologically important shifts in mutational effects.
Collapse
Affiliation(s)
- Hugh K. Haddox
- Computational Biology Program, Fred Hutchinson Cancer Research Center, Seattle, WA 98102, USA
| | - Jared G. Galloway
- Computational Biology Program, Fred Hutchinson Cancer Research Center, Seattle, WA 98102, USA
| | - Bernadeta Dadonaite
- Basic Sciences Division, Fred Hutchinson Cancer Research Center, Seattle, WA 98109, USA
| | - Jesse D. Bloom
- Computational Biology Program, Fred Hutchinson Cancer Research Center, Seattle, WA 98102, USA
- Basic Sciences Division, Fred Hutchinson Cancer Research Center, Seattle, WA 98109, USA
- Department of Genome Sciences, University of Washington, Seattle, WA 98195, USA
- Howard Hughes Medical Institute, Seattle, WA 98109, USA
| | - Frederick A. Matsen IV
- Computational Biology Program, Fred Hutchinson Cancer Research Center, Seattle, WA 98102, USA
- Basic Sciences Division, Fred Hutchinson Cancer Research Center, Seattle, WA 98109, USA
- Howard Hughes Medical Institute, Seattle, WA 98109, USA
- Department of Statistics, University of Washington, Seattle, WA 98195, USA
| | - William S. DeWitt
- Department of Electrical Engineering and Computer Sciences, University of California, Berkeley, CA 94720, USA
| |
Collapse
|