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.

Digitizing omics profiles by divergence from a baseline.

Data collected from omics technologies have revealed pervasive heterogeneity and stochasticity of molecular states within and between phenotypes. A prominent example of such heterogeneity occurs between genome-wide mRNA, microRNA, and methylation profiles from one individual tumor to another, even within a cancer subtype. However, current methods in bioinformatics, such as detecting differentially expressed genes or CpG sites, are population-based and therefore do not effectively model intersample diversity. Here we introduce a unified theory to quantify sample-level heterogeneity that is applicable to a single omics profile. Specifically, we simplify an omics profile to a digital representation based on the omics profiles from a set of samples from a reference or baseline population (e.g., normal tissues). The state of any subprofile (e.g., expression vector for a subset of genes) is said to be "divergent" if it lies outside the estimated support of the baseline distribution and is consequently interpreted as "dysregulated" relative to that baseline. We focus on two cases: single features (e.g., individual genes) and distinguished subsets (e.g., regulatory pathways). Notably, since the divergence analysis is at the individual sample level, dysregulation can be analyzed probabilistically; for example, one can estimate the probability that a gene or pathway is divergent in some population. Finally, the reduction in complexity facilitates a more "personalized" and biologically interpretable analysis of variation, as illustrated by experiments involving tissue characterization, disease detection and progression, and disease-pathway associations.

Determining mesoscale chromatin structure parameters from spatially correlated cleavage data using a coarse-grained oligonucleosome model.

The three-dimensional structure of chromatin has emerged as an important feature of eukaryotic gene regulation. Recent technological advances in DNA sequencing-based assays have revealed locus- and chromatin state-specific structural patterns at the length scale of a few nucleosomes (~1 kb). However, interpreting these data sets remains challenging. Radiation-induced correlated cleavage of chromatin (RICC-seq) is one such chromatin structure assay that maps DNA-DNA-contacts at base pair resolution by sequencing single-stranded DNA fragments released from irradiated cells. Here, we develop a flexible modeling and simulation framework to enable the interpretation of RICC-seq data in terms of oligonucleosome structure ensembles. Nucleosomes are modeled as rigid bodies with excluded volume and adjustable DNA wrapping, connected by linker DNA modeled as a worm-like chain. We validate the model's parameters against cryo-electron microscopy and sedimentation data. Our results show that RICC-seq is sensitive to nucleosome spacing, nucleosomal DNA wrapping, and the strength of inter-nucleosome interactions. We show that nucleosome repeat lengths consistent with orthogonal assays can be extracted from experimental RICC-seq data using a 1D convolutional neural net trained on RICC-seq signal predicted from simulated ensembles. We thus provide a suite of analysis tools that add quantitative structural interpretability to RICC-seq experiments.

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.

DRUL for school: Opening Pre-K with safe, simple, sensitive saliva testing for SARS-CoV-2.

To address the need for simple, safe, sensitive, and scalable SARS-CoV-2 tests, we validated and implemented a PCR test that uses a saliva collection kit use at home. Individuals self-collected 300 µl saliva in vials containing Darnell Rockefeller University Laboratory (DRUL) buffer and extracted RNA was assayed by RT-PCR (the DRUL saliva assay). The limit of detection was confirmed to be 1 viral copy/µl in 20 of 20 replicate extractions. Viral RNA was stable in DRUL buffer at room temperature up to seven days after sample collection, and safety studies demonstrated that DRUL buffer immediately inactivated virus at concentrations up to 2.75x106 PFU/ml. Results from SARS-CoV-2 positive nasopharyngeal (NP) swab samples collected in viral transport media and assayed with a standard FDA Emergency Use Authorization (EUA) test were highly correlated with samples placed in DRUL buffer. Direct comparison of results from 162 individuals tested by FDA EUA oropharyngeal (OP) or NP swabs with co-collected saliva samples identified four otherwise unidentified positive cases in DRUL buffer. Over six months, we collected 3,724 samples from individuals ranging from 3 months to 92 years of age. This included collecting weekly samples over 10 weeks from teachers, children, and parents from a pre-school program, which allowed its safe reopening while at-risk pods were quarantined. In sum, we validated a simple, sensitive, stable, and safe PCR-based test using a self-collected saliva sample as a valuable tool for clinical diagnosis and screening at workplaces and schools.