Engineering topics

Ke, Ruian

Publications and source records attributed to Ke, Ruian.

At least 19 records

DOE OSTI2025

Identifying impacts of contact tracing on HIV epidemiological inference from phylogenetic data

Abstract Robust sampling methods are foundational to inferences using phylogenies. Yet the impact of using contact tracing, a type of non-uniform sampling used in public health applications such as infectious disease outbreak investigations, has not been investigated in the molecular epidemiology field. To understand how contact tracing influences a recovered phylogeny, we developed a new simulation tool called SEEPS (Sequence Evolution and Epidemiological Process Simulator) that allows for the simulation of contact tracing and the resulting transmission tree, pathogen phylogeny, and corresponding virus genetic sequences. Importantly, SEEPS takes within-host evolution into account when generating pathogen phylogenies and sequences from transmission histories. Using SEEPS, we demonstrate that contact tracing can significantly impact the structure of the resulting tree, as described by popular tree statistics. Contact tracing generates phylogenies that are less balanced than the underlying transmission process, less representative of the larger epidemiological process, and affects the internal/external branch length ratios that characterize specific epidemiological scenarios. We also examined real data from a 2007–2008 Swedish HIV-1 outbreak and the broader 1998–2010 European HIV-1 epidemic to highlight the differences in contact tracing and expected phylogenies. Aided by SEEPS, we show that the data collection of the Swedish outbreak was strongly influenced by contact tracing even after downsampling, while the broader European Union epidemic showed little evidence of universal contact tracing, agreeing with the known epidemiological information about sampling and spread. Overall, our results highlight the importance of including possible non-uniform sampling schemes when examining phylogenetic trees. For that, SEEPS serves as a useful tool to evaluate such impacts, thereby facilitating better phylogenetic inferences of the characteristics of a disease outbreak. SEEPS is available at https://github.com/MolEvolEpid/SEEPS.

Virology↗

DOE OSTI2024

CovTransformer: A transformer model for SARS-CoV-2 lineage frequency forecasting

With hundreds of SARS-CoV-2 lineages circulating in the global population, there is an ongoing need for predicting and forecasting lineage frequencies and thus identifying rapidly expanding lineages. Accurate prediction would allow for more focused experimental efforts to understand pathogenicity of future dominating lineages and characterize the extent of their immune escape. Here, we first show that the inherent noise and biases in lineage frequency data make a commonly-used regression-based approach unreliable. To address this weakness, we constructed a machine learning model for SARS-CoV-2 lineage frequency forecasting, called CovTransformer, based on the transformer architecture. We designed our model to navigate challenges such as a limited amount of data with high levels of noise and bias. We first trained and tested the model using data from the UK and the USA, and then tested the generalization ability of the model to many other countries and US states. Remarkably, the trained model makes accurate predictions two months into the future with high levels of accuracy both globally (in 31 countries with high levels of sequencing effort) and at the US-state level. Our model performed substantially better than a widely used forecasting tool, the multinomial regression model implemented in Nextstrain, demonstrating its utility in SARS-CoV-2 monitoring. Assuming a newly emerged lineage is identified and assigned, our test using retrospective data shows that our model is able to identify the dominating lineages 7 weeks in advance on average before they became dominant. Overall, our work demonstrates that transformer models represent a promising approach for SARS-CoV-2 forecasting and pandemic monitoring.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

The kinetics of SARS-CoV-2 infection based on a human challenge study

Studying the early events that occur after viral infection in humans is difficult unless one intentionally infects volunteers in a human challenge study. Here, we use data about severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in such a study in combination with mathematical modeling to gain insights into the relationship between the amount of virus in the upper respiratory tract and the immune response it generates. We propose a set of dynamic models of increasing complexity to dissect the roles of target cell limitation, innate immunity, and adaptive immunity in determining the observed viral kinetics. We introduce an approach for modeling the effect of humoral immunity that describes a decline in infectious virus after immune activation. We fit our models to viral load and infectious titer data from all the untreated infected participants in the study simultaneously. We found that a power-law with a power h < 1 describes the relationship between infectious virus and viral load. Viral replication at the early stage of infection is rapid, with a doubling time of ~2 h for viral RNA and ~3 h for infectious virus. We estimate that adaptive immunity is initiated ~7 to 10 d postinfection and appears to contribute to a multiphasic viral decline experienced by some participants; the viral rebound experienced by other participants is consistent with a decline in the interferon response. Altogether, we quantified the kinetics of SARS-CoV-2 infection, shedding light on the early dynamics of the virus and the potential role of innate and adaptive immunity in promoting viral decline during infection.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

CoVTransformer

The code is a transformer model to forecast SARS-CoV-2 lineage frequencies in the future.

Feng, Yinan↗

DOE OSTI2024

Understanding early HIV-1 rebound dynamics following antiretroviral therapy interruption: The importance of effector cell expansion

Most people living with HIV-1 experience rapid viral rebound once antiretroviral therapy is interrupted; however, a small fraction remain in viral remission for an extended duration. Understanding the factors that determine whether viral rebound is likely after treatment interruption can enable the development of optimal treatment regimens and therapeutic interventions to potentially achieve a functional cure for HIV-1. We built upon the theoretical framework proposed by Conway and Perelson to construct dynamic models of virus-immune interactions to study factors that influence viral rebound dynamics. We evaluated these models using viral load data from 24 individuals following antiretroviral therapy interruption. The best-performing model accurately captures the heterogeneity of viral dynamics and highlights the importance of the effector cell expansion rate. Our results show that post-treatment controllers and non-controllers can be distinguished based on the effector cell expansion rate in our models. Furthermore, these results demonstrate the potential of using dynamic models incorporating an effector cell response to understand early viral rebound dynamics post-antiretroviral therapy interruption.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

Influenza A viral burst size from thousands of infected single cells using droplet quantitative PCR (dqPCR)

An important aspect of how viruses spread and infect is the viral burst size, or the number of new viruses produced by each infected cell. Surprisingly, this value remains poorly characterized for influenza A virus (IAV), commonly known as the flu. In this study, we screened tens of thousands of cells using a microfluidic method called droplet quantitative PCR (dqPCR). The high-throughput capability of dqPCR enabled the measurement of a large population of infected cells producing progeny virus. By measuring the fully assembled and successfully released viruses from these infected cells, we discover that the viral burst sizes for both the seasonal H3N2 and the 2009 pandemic H1N1 strains vary significantly, with H3N2 ranging from 10 1 to 10 4 viruses per cell, and H1N1 ranging from 10 1 to 10 3 viruses per cell. Some infected cells produce average numbers of new viruses, while others generate extensive number of viruses. In fact, we find that only 10% of the single-cell infections are responsible for creating a significant portion of all the viruses. This small fraction produced approximately 60% of new viruses for H3N2 and 40% for H1N1. On average, each infected cell of the H3N2 flu strain produced 709 new viruses, whereas for H1N1, each infected cell produced 358 viruses. This novel method reveals insights into the flu virus and can lead to improved strategies for managing and preventing the spread of viruses.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

Modeling the emergence of viral resistance for SARS-CoV-2 during treatment with an anti-spike monoclonal antibody

To mitigate the loss of lives during the COVID-19 pandemic, emergency use authorization was given to several anti-SARS-CoV-2 monoclonal antibody (mAb) therapies for the treatment of mild-to-moderate COVID-19 in patients with a high risk of progressing to severe disease. Monoclonal antibodies used to treat SARS-CoV-2 target the spike protein of the virus and block its ability to enter and infect target cells. Monoclonal antibody therapy can thus accelerate the decline in viral load and lower hospitalization rates among high-risk patients with variants susceptible to mAb therapy. However, viral resistance has been observed, in some cases leading to a transient viral rebound that can be as large as 3–4 orders of magnitude. As mAbs represent a proven treatment choice for SARS-CoV-2 and other viral infections, evaluation of treatment-emergent mAb resistance can help uncover underlying pathobiology of SARS-CoV-2 infection and may also help in the development of the next generation of mAb therapies. Although resistance can be expected, the large rebounds observed are much more difficult to explain. We hypothesize replenishment of target cells is necessary to generate the high transient viral rebound. Thus, we formulated two models with different mechanisms for target cell replenishment (homeostatic proliferation and return from an innate immune response antiviral state) and fit them to data from persons with SARS-CoV-2 treated with a mAb. We showed that both models can explain the emergence of resistant virus associated with high transient viral rebounds. We found that variations in the target cell supply rate and adaptive immunity parameters have a strong impact on the magnitude or observability of the viral rebound associated with the emergence of resistant virus. Both variations in target cell supply rate and adaptive immunity parameters may explain why only some individuals develop observable transient resistant viral rebound. Our study highlights the conditions that can lead to resistance and subsequent viral rebound in mAb treatments during acute infection.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

How robust are estimates of key parameters in standard viral dynamic models?

Mathematical models of viral infection have been developed, fitted to data, and provide insight into disease pathogenesis for multiple agents that cause chronic infection, including HIV, hepatitis C, and B virus. However, for agents that cause acute infections or during the acute stage of agents that cause chronic infections, viral load data are often collected after symptoms develop, usually around or after the peak viral load. Consequently, we frequently lack data in the initial phase of viral growth, i.e., when pre-symptomatic transmission events occur. Missing data may make estimating the time of infection, the infectious period, and parameters in viral dynamic models, such as the cell infection rate, difficult. However, having extra information, such as the average time to peak viral load, may improve the robustness of the estimation. Here, we evaluated the robustness of estimates of key model parameters when viral load data prior to the viral load peak is missing, when we know the values of some parameters and/or the time from infection to peak viral load. Although estimates of the time of infection are sensitive to the quality and amount of available data, particularly pre-peak, other parameters important in understanding disease pathogenesis, such as the loss rate of infected cells, are less sensitive. Viral infectivity and the viral production rate are key parameters affecting the robustness of data fits. Fixing their values to literature values can help estimate the remaining model parameters when pre-peak data is missing or limited. We find a lack of data in the pre-peak growth phase underestimates the time to peak viral load by several days, leading to a shorter predicted growth phase. On the other hand, knowing the time of infection (e.g., from epidemiological data) and fixing it results in good estimates of dynamical parameters even in the absence of early data. While we provide ways to approximate model parameters in the absence of early viral load data, our results also suggest that these data, when available, are needed to estimate model parameters more precisely.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

AZD5582 plus SIV-specific antibodies reduce lymph node viral reservoirs in antiretroviral therapy-suppressed macaques

The main barrier to HIV cure is a persistent reservoir of latently infected CD4 + T cells harboring replication-competent provirus that fuels rebound viremia upon antiretroviral therapy (ART) interruption. A leading approach to target this reservoir involves agents that reactivate latent HIV proviruses followed by direct clearance of cells expressing induced viral antigens by immune effector cells and immunotherapeutics. We previously showed that AZD5582, an antagonist of inhibitor of apoptosis proteins and mimetic of the second mitochondrial-derived activator of caspases (IAPi/SMACm), induces systemic reversal of HIV/SIV latency but with no reduction in size of the viral reservoir. In this study, we investigated the effects of AZD5582 in combination with four SIV Env-specific Rhesus monoclonal antibodies (RhmAbs) ± N-803 (an IL-15 superagonist) in SIV-infected, ART-suppressed rhesus macaques. Here we confirm the efficacy of AZD5582 in inducing SIV reactivation, demonstrate enhancement of latency reversal when AZD5582 is used in combination with N-803 and show a reduction in total and replication-competent SIV-DNA in lymph-node-derived CD4 + T cells in macaques treated with AZD5582 + RhmAbs. Further exploration of this therapeutic approach may contribute to the goal of achieving an HIV cure.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2023

Swift and extensive Omicron outbreak in China after sudden exit from ‘zero-COVID’ policy

In late 2022, China transitioned from a strict ‘zero-COVID’ policy to rapidly abandoning nearly all interventions and data reporting. This raised great concern about the presumably-rapid but unreported spread of the SARS-CoV-2 Omicron variant in a very large population of very low pre-existing immunity. By modeling a combination of case count and survey data, we show that Omicron spread extremely rapidly, at a rate of 0.42/day (95% credibility interval: [0.35, 0.51]/day), translating to an epidemic doubling time of 1.6 days ([1.6, 2.0] days) after the full exit from zero-COVID on Dec. 7, 2022. Consequently, we estimate that the vast majority of the population (97% [95%, 99%], sensitivity analysis lower limit of 90%) was infected during December, with the nation-wide epidemic peaking on Dec. 23. Overall, our results highlight the extremely high transmissibility of the variant and the importance of proper design of intervention exit strategies to avoid large infection waves.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

A deep learning approach to real-time HIV outbreak detection using genetic data

Pathogen genomic sequence data are increasingly made available for epidemiological monitoring. A main interest is to identify and assess the potential of infectious disease outbreaks. While popular methods to analyze sequence data often involve phylogenetic tree inference, they are vulnerable to errors from recombination and impose a high computational cost, making it difficult to obtain real-time results when the number of sequences is in or above the thousands. Here, we propose an alternative strategy to outbreak detection using genomic data based on deep learning methods developed for image classification. The key idea is to use a pairwise genetic distance matrix calculated from viral sequences as an image, and develop convolutional neutral network (CNN) models to classify areas of the images that show signatures of active outbreak, leading to identification of subsets of sequences taken from an active outbreak. We showed that our method is efficient in finding HIV-1 outbreaks with R0 ≥ 2.5, and overall a specificity exceeding 98% and sensitivity better than 92%. We validated our approach using data from HIV-1 CRF01 in Europe, containing both endemic sequences and a well-known dual outbreak in intravenous drug users. Our model accurately identified known outbreak sequences in the background of slower spreading HIV. Importantly, we detected both outbreaks early on, before they were over, implying that had this method been applied in real-time as data became available, one would have been able to intervene and possibly prevent the extent of these outbreaks. This approach is scalable to processing hundreds of thousands of sequences, making it useful for current and future real-time epidemiological investigations, including public health monitoring using large databases and especially for rapid outbreak identification.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

A simple model of COVID-19 explains disease severity and the effect of treatments

Abstract Considerable effort has been made to better understand why some people suffer from severe COVID-19 while others remain asymptomatic. This has led to important clinical findings; people with severe COVID-19 generally experience persistently high levels of inflammation, slower viral load decay, display a dysregulated type-I interferon response, have less active natural killer cells and increased levels of neutrophil extracellular traps. How these findings are connected to the pathogenesis of COVID-19 remains unclear. We propose a mathematical model that sheds light on this issue by focusing on cells that trigger inflammation through molecular patterns: infected cells carrying pathogen-associated molecular patterns (PAMPs) and damaged cells producing damage-associated molecular patterns (DAMPs). The former signals the presence of pathogens while the latter signals danger such as hypoxia or lack of nutrients. Analyses show that SARS-CoV-2 infections can lead to a self-perpetuating feedback loop between DAMP expressing cells and inflammation, identifying the inability to quickly clear PAMPs and DAMPs as the main contributor to hyperinflammation. The model explains clinical findings and reveal conditions that can increase the likelihood of desired clinical outcome from treatment administration. In particular, the analysis suggest that antivirals need to be administered early during infection to have an impact on disease severity. The simplicity of the model and its high level of consistency with clinical findings motivate its use for the formulation of new treatment strategies.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Multi-species outbreak of SARS-CoV-2 Delta variant in a zoological institution, with the detection in two new families of carnivores

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has a worldwide distribution in humans and many other mammalian species. In late September 2021, 12 animals maintained by the Chicago Zoological Society's Brookfield Zoo were observed with variable clinical signs. The Delta variant of SARS-CoV-2 was detected in faeces and nasal swabs by qRT-PCR, including the first detection in animals from the families Procyonidae and Viverridae. Test positivity rate was 12.5% for 35 animals tested. All animals had been vaccinated with at least one dose of a recombinant vaccine designed for animals and all recovered with variable supportive treatment. Sequence analysis showed that six zoo animal strains were closely correlated with 18 human SARS-CoV-2 strains, suggestive of potential human-to-animal transmission events. This report documents the expanding host range of COVID-19 during the ongoing pandemic.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Disrupting autorepression circuitry generates “open-loop lethality” to yield escape-resistant antiviral agents

Across biological scales, gene-regulatory networks employ autorepression (negative feedback) to maintain homeostasis and minimize failure from aberrant expression. In this study, we present a proof of concept that disrupting transcriptional negative feedback dysregulates viral gene expression to therapeutically inhibit replication and confers a high evolutionary barrier to resistance. We find that nucleic-acid decoys mimicking cis-regulatory sites act as “feedback disruptors,” break homeostasis, and increase viral transcription factors to cytotoxic levels (termed “open-loop lethality”). Feedback disruptors against herpesviruses reduced viral replication >2-logs without activating innate immunity, showed sub-nM IC 50 , synergized with standard-of-care antivirals, and inhibited virus replication in mice. In contrast to approved antivirals where resistance rapidly emerged, no feedback-disruptor escape mutants evolved in long-term cultures. For SARS-CoV-2, disruption of a putative feedback circuit also generated open-loop lethality, reducing viral titers by >1-log. These results demonstrate that generating open-loop lethality, via negative-feedback disruption, may yield a class of antimicrobials with a high genetic barrier to resistance.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Longitudinal Analysis of SARS-CoV-2 Vaccine Breakthrough Infections Reveals Limited Infectious Virus Shedding and Restricted Tissue Distribution

Background. The global effort to vaccinate people against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) during an ongoing pandemic has raised questions about how vaccine breakthrough infections compare with infections in immunologically naive individuals and the potential for vaccinated individuals to transmit the virus. Methods. We examined viral dynamics and infectious virus shedding through daily longitudinal sampling in 23 adults infected with SARS-CoV-2 at varying stages of vaccination, including 6 fully vaccinated individuals. Results. The durations of both infectious virus shedding and symptoms were significantly reduced in vaccinated individuals compared with unvaccinated individuals. We also observed that breakthrough infections are associated with strong tissue compartmentalization and are only detectable in saliva in some cases. Conclusions. Vaccination shortens the duration of time of high transmission potential, minimizes symptom duration, and may restrict tissue dissemination.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Assessing the impact of human mobility to predict regional excess death in Ecuador

Abstract COVID-19 outbreaks have had high mortality in low- and middle-income countries such as Ecuador. Human mobility is an important factor influencing the spread of diseases possibly leading to a high burden of disease at the country level. Drastic control measures, such as complete lockdown, are effective epidemic controls, yet in practice one hopes that a partial shutdown would suffice. It is an open problem to determine how much mobility can be allowed while controlling an outbreak. In this paper, we use statistical models to relate human mobility to the excess death in Ecuador while controlling for demographic factors. The mobility index provided by GRANDATA, based on mobile phone users, represents the change of number of out-of-home events with respect to a benchmark date (March 2nd, 2020). The study confirms the global trend that more men are dying than expected compared to women, and that people under 30 show less deaths than expected, particularly individuals younger than 20 with a death rate reduction between 22 and 27%. The weekly median mobility time series shows a sharp decrease in human mobility immediately after a national lockdown was declared on March 17, 2020 and a progressive increase towards the pre-lockdown level within two months. Relating median mobility to excess deaths shows a lag in its effect: first, a decrease in mobility in the previous two to three weeks decreases excess death and, more novel, we found an increase of mobility variability four weeks prior increases the number of excess deaths.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Estimating the strength of selection for new SARS-CoV-2 variants

Abstract Controlling the SARS-CoV-2 pandemic becomes increasingly challenging as the virus adapts to human hosts through the continual emergence of more transmissible variants. Simply observing that a variant is increasing in frequency is relatively straightforward, but more sophisticated methodology is needed to determine whether a new variant is a global threat and the magnitude of its selective advantage. We present two models for quantifying the strength of selection for new and emerging variants of SARS-CoV-2 relative to the background of contemporaneous variants. These methods range from a detailed model of dynamics within one country to a broad analysis across all countries, and they include alternative explanations such as migration and drift. We find evidence for strong selection favoring the D614G spike mutation and B.1.1.7 (Alpha), weaker selection favoring B.1.351 (Beta), and no advantage of R.1 after it spreads beyond Japan. Cutting back data to earlier time horizons reveals that uncertainty is large very soon after emergence, but that estimates of selection stabilize after several weeks. Our results also show substantial heterogeneity among countries, demonstrating the need for a truly global perspective on the molecular epidemiology of SARS-CoV-2.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

In vivo kinetics of SARS-CoV-2 infection and its relationship with a person’s infectiousness

Significance Quantifying the kinetics of SARS-CoV-2 infection and individual infectiousness is important for understanding SARS-CoV-2 transmission and evaluating intervention strategies. Here, we developed within-host models of SARS-CoV-2 infection, and by fitting them to clinical data, we estimated key within-host viral dynamic parameters. We also developed a mechanistic model for viral transmission and show that the logarithm of the viral load in the upper respiratory tract serves as an appropriate surrogate for a person’s infectiousness. Using data on how viral load changes during infection, we further evaluated the effectiveness of PCR and antigen-based testing strategies for averting transmission and identifying infected individuals.

60 APPLIED LIFE SCIENCES↗