An Amygdala Circuit Mediates Experience-Dependent Momentary Arrests during Exploration.

Exploration of novel environments ensures survival and evolutionary fitness. It is expressed through exploratory bouts and arrests that change dynamically based on experience. Neural circuits mediating exploratory behavior should therefore integrate experience and use it to select the proper behavioral output. Using a spatial exploration assay, we uncovered an experience-dependent increase in momentary arrests in locations where animals arrested previously. Calcium imaging in freely exploring mice revealed a genetically and projection-defined neuronal ensemble in the basolateral amygdala that is active during self-paced behavioral arrests. This ensemble was recruited in an experience-dependent manner, and closed-loop optogenetic manipulation of these neurons revealed that they are sufficient and necessary to drive experience-dependent arrests during exploration. Projection-specific imaging and optogenetic experiments revealed that these arrests are effected by basolateral amygdala neurons projecting to the central amygdala, uncovering an amygdala circuit that mediates momentary arrests in familiar places but not avoidance or anxiety/fear-like behaviors.

Single-cell epigenetic, transcriptional, and protein profiling of latent and active HIV-1 reservoir revealed that IKZF3 promotes HIV-1 persistence.

Understanding how HIV-1-infected cells proliferate and persist is key to HIV-1 eradication, but the heterogeneity and rarity of HIV-1-infected cells hamper mechanistic interrogations. Here, we used single-cell DOGMA-seq to simultaneously capture transcription factor accessibility, transcriptome, surface proteins, HIV-1 DNA, and HIV-1 RNA in memory CD4+ T cells from six people living with HIV-1 during viremia and after suppressive antiretroviral therapy. We identified increased transcription factor accessibility in latent HIV-1-infected cells (RORC) and transcriptionally active HIV-1-infected cells (interferon regulatory transcription factor [IRF] and activator protein 1 [AP-1]). A proliferation program (IKZF3, IL21, BIRC5, and MKI67 co-expression) promoted the survival of transcriptionally active HIV-1-infected cells. Both latent and transcriptionally active HIV-1-infected cells had increased IKZF3 (Aiolos) expression. Distinct epigenetic programs drove the heterogeneous cellular states of HIV-1-infected cells: IRF:activation, Eomes:cytotoxic effector differentiation, AP-1:migration, and cell death. Our study revealed the single-cell epigenetic, transcriptional, and protein states of latent and transcriptionally active HIV-1-infected cells and cellular programs promoting HIV-1 persistence.

Single-cell multiomics reveals persistence of HIV-1 in expanded cytotoxic T cell clones.

Understanding the drivers and markers of clonally expanding HIV-1-infected CD4+ T cells is essential for HIV-1 eradication. We used single-cell ECCITE-seq, which captures surface protein expression, cellular transcriptome, HIV-1 RNA, and TCR sequences within the same single cell to track clonal expansion dynamics in longitudinally archived samples from six HIV-1-infected individuals (during viremia and after suppressive antiretroviral therapy) and two uninfected individuals, in unstimulated conditions and after CMV and HIV-1 antigen stimulation. Despite antiretroviral therapy, persistent antigen and TNF responses shaped T cell clonal expansion. HIV-1 resided in Th1-polarized, antigen-responding T cells expressing BCL2 and SERPINB9 that may resist cell death. HIV-1 RNA+ T cell clones were larger in clone size, established during viremia, persistent after viral suppression, and enriched in GZMB+ cytotoxic effector memory Th1 cells. Targeting HIV-1-infected cytotoxic CD4+ T cells and drivers of clonal expansion provides another direction for HIV-1 eradication.

Brain Repair by Cell Replacement via In Situ Neuronal Reprogramming.

Neurodegenerative diseases, characterized by progressive neural loss, have been some of the most challenging medical problems in aging societies. Treatment strategies such as symptom management have little impact on disease progression, while intervention with specific disease mechanisms may only slow down disease progression. One therapeutic strategy that has the potential to reverse the disease phenotype is to replenish neurons and rebuild the pathway lost to degeneration. Although it is generally believed that the central nervous system has lost the capability to regenerate, increasing evidence indicates that the brain is more plastic than previously thought, containing perhaps the biggest repertoire of cells with latent neurogenic programs in the body. This review focuses on key advances in generating new neurons through in situ neuronal reprogramming, which is tied to fundamental questions regarding adult neurogenesis, cell source, and mechanisms for neuronal reprogramming, as well as the ability of new neurons to integrate into the existing circuitry.

Unraveling the mesoscale organization induced by network-driven processes.

Complex systems are characterized by emergent patterns created by the nontrivial interplay between dynamical processes and the networks of interactions on which these processes unfold. Topological or dynamical descriptors alone are not enough to fully embrace this interplay in all its complexity, and many times one has to resort to dynamics-specific approaches that limit a comprehension of general principles. To address this challenge, we employ a metric-that we name Jacobian distance-which captures the spatiotemporal spreading of perturbations, enabling us to uncover the latent geometry inherent in network-driven processes. We compute the Jacobian distance for a broad set of nonlinear dynamical models on synthetic and real-world networks of high interest for applications from biological to ecological and social contexts. We show, analytically and computationally, that the process-driven latent geometry of a complex network is sensitive to both the specific features of the dynamics and the topological properties of the network. This translates into potential mismatches between the functional and the topological mesoscale organization, which we explain by means of the spectrum of the Jacobian matrix. Finally, we demonstrate that the Jacobian distance offers a clear advantage with respect to traditional methods when studying human brain networks. In particular, we show that it outperforms classical network communication models in explaining functional communities from structural data, therefore highlighting its potential in linking structure and function in the brain.

Stochastic machine learning via sigma profiles to build a digital chemical space.

This work establishes a different paradigm on digital molecular spaces and their efficient navigation by exploiting sigma profiles. To do so, the remarkable capability of Gaussian processes (GPs), a type of stochastic machine learning model, to correlate and predict physicochemical properties from sigma profiles is demonstrated, outperforming state-of-the-art neural networks previously published. The amount of chemical information encoded in sigma profiles eases the learning burden of machine learning models, permitting the training of GPs on small datasets which, due to their negligible computational cost and ease of implementation, are ideal models to be combined with optimization tools such as gradient search or Bayesian optimization (BO). Gradient search is used to efficiently navigate the sigma profile digital space, quickly converging to local extrema of target physicochemical properties. While this requires the availability of pretrained GP models on existing datasets, such limitations are eliminated with the implementation of BO, which can find global extrema with a limited number of iterations. A remarkable example of this is that of BO toward boiling temperature optimization. Holding no knowledge of chemistry except for the sigma profile and boiling temperature of carbon monoxide (the worst possible initial guess), BO finds the global maximum of the available boiling temperature dataset (over 1,000 molecules encompassing more than 40 families of organic and inorganic compounds) in just 15 iterations (i.e., 15 property measurements), cementing sigma profiles as a powerful digital chemical space for molecular optimization and discovery, particularly when little to no experimental data is initially available.

Bayesian functional analysis for untargeted metabolomics data with matching uncertainty and small sample sizes.

Untargeted metabolomics based on liquid chromatography-mass spectrometry technology is quickly gaining widespread application, given its ability to depict the global metabolic pattern in biological samples. However, the data are noisy and plagued by the lack of clear identity of data features measured from samples. Multiple potential matchings exist between data features and known metabolites, while the truth can only be one-to-one matches. Some existing methods attempt to reduce the matching uncertainty, but are far from being able to remove the uncertainty for most features. The existence of the uncertainty causes major difficulty in downstream functional analysis. To address these issues, we develop a novel approach for Bayesian Analysis of Untargeted Metabolomics data (BAUM) to integrate previously separate tasks into a single framework, including matching uncertainty inference, metabolite selection and functional analysis. By incorporating the knowledge graph between variables and using relatively simple assumptions, BAUM can analyze datasets with small sample sizes. By allowing different confidence levels of feature-metabolite matching, the method is applicable to datasets in which feature identities are partially known. Simulation studies demonstrate that, compared with other existing methods, BAUM achieves better accuracy in selecting important metabolites that tend to be functionally consistent and assigning confidence scores to feature-metabolite matches. We analyze a COVID-19 metabolomics dataset and a mouse brain metabolomics dataset using BAUM. Even with a very small sample size of 16 mice per group, BAUM is robust and stable. It finds pathways that conform to existing knowledge, as well as novel pathways that are biologically plausible.

GLDADec: marker-gene guided LDA modeling for bulk gene expression deconvolution.

Inferring cell type proportions from bulk transcriptome data is crucial in immunology and oncology. Here, we introduce guided LDA deconvolution (GLDADec), a bulk deconvolution method that guides topics using cell type-specific marker gene names to estimate topic distributions for each sample. Through benchmarking using blood-derived datasets, we demonstrate its high estimation performance and robustness. Moreover, we apply GLDADec to heterogeneous tissue bulk data and perform comprehensive cell type analysis in a data-driven manner. We show that GLDADec outperforms existing methods in estimation performance and evaluate its biological interpretability by examining enrichment of biological processes for topics. Finally, we apply GLDADec to The Cancer Genome Atlas tumor samples, enabling subtype stratification and survival analysis based on estimated cell type proportions, thus proving its practical utility in clinical settings. This approach, utilizing marker gene names as partial prior information, can be applied to various scenarios for bulk data deconvolution. GLDADec is available as an open-source Python package at https://github.com/mizuno-group/GLDADec.

CTLA-4+PD-1- Memory CD4+ T Cells Critically Contribute to Viral Persistence in Antiretroviral Therapy-Suppressed, SIV-Infected Rhesus Macaques.

Antiretroviral therapy (ART) suppresses viral replication in HIV-infected individuals but does not eliminate the reservoir of latently infected cells. Recent work identified PD-1+ follicular helper T (Tfh) cells as an important cellular compartment for viral persistence. Here, using ART-treated, SIV-infected rhesus macaques, we show that CTLA-4+PD-1- memory CD4+ T cells, which share phenotypic markers with regulatory T cells, were enriched in SIV DNA in blood, lymph nodes (LN), spleen, and gut, and contained replication-competent and infectious virus. In contrast to PD-1+ Tfh cells, SIV-enriched CTLA-4+PD-1- CD4+ T cells were found outside the B cell follicle of the LN, predicted the size of the persistent viral reservoir during ART, and significantly increased their contribution to the SIV reservoir with prolonged ART-mediated viral suppression. We have shown that CTLA-4+PD-1- memory CD4+ T cells are a previously unrecognized component of the SIV and HIV reservoir that should be therapeutically targeted for a functional HIV-1 cure.

Biphasic decay of intact SHIV genomes following initiation of antiretroviral therapy complicates analysis of interventions targeting the reservoir.

The latent reservoir for HIV-1 in resting CD4+ T cells persists despite antiretroviral therapy (ART) and precludes cure. Reservoir-targeting interventions are evaluated in ART-treated macaques infected with simian immunodeficiency virus (SIV) or simian-human immunodeficiency virus (SHIV). Efficacy is determined by reservoir measurements before and after the intervention. However, most proviruses persisting in the setting of ART are defective. In addition, intact HIV-1 and SIV genomes undergo complex, multiphasic decay observable when new infection events are blocked by ART. Intervention-induced elimination of latently infected cells must be distinguished from natural decay. Here, we address these issues for SHIV. We describe an intact proviral DNA assay that allows digital counting of SHIV genomes lacking common fatal defects. We show that intact SHIV genomes in circulating CD4+ T cells undergo biphasic decay during the first year of ART, with a rapid first phase (t1/2 = 30.1 d) and a slower second phase (t1/2 = 8.1 mo) that is still more rapid that the slow decay observed in people with HIV-1 on long-term ART (t1/2 = 3.7 y). In SHIV models, most interventions are tested during 2nd phase decay. Natural 2nd phase decay must be considered in evaluating interventions as most infected cells present at this time do not become part of the stable reservoir. In addition, for interventions tested during 2nd phase decay, a caveat is that the intervention may not be equally effective in people with HIV on long-term ART whose reservoirs are dominated by latently infected cells with a slower decay rate.

Ants combine object affordance with latent learning to make efficient foraging decisions.

The affordance of an object refers to its functional properties. For example, a bowl has the affordance of holding water, but a sieve does not. Here, we report that ants learn the affordance of a novel object without this attribute being rewarded, and use the memory of this affordance to avoid predicted, but never experienced, crowding. Ants were trained to feeders, which could support either only one ant or many. Two feeders were encountered, each of identical design but differently scented. After training, on the outward journey, half the ants encounter nestmates, which had fed on food matching one of the training feeders. Encountering returning nestmates reduced preference for the feeder matching the scent of the encountered nestmates, but only for ants trained on a limited-access feeder; ants trained on an unlimited feeder were unaffected. In other words, only if ants know that the food access is limited, and receive information that this feeder is heavily visited, do they reduce their preference for this feeder. To achieve this, the ants had to combine memories of the feeders' affordance with the presence of nestmates. Then they had to use semantic knowledge that restricted food access combined with nestmate presence predicts a likelihood of crowding, or a rule such as "if the food is restricted and there are nestmates on the path, go to another food source." Regardless of the mechanism, these results demonstrate that ants latently learn the affordance of their surroundings, an unexpected cognitive ability for an invertebrate.

Identifying Distinct Neural Features between the Initial and Corrective Phases of Precise Reaching Using AutoLFADS.

Many initial movements require subsequent corrective movements, but how the motor cortex transitions to make corrections and how similar the encoding is to initial movements is unclear. In our study, we explored how the brain's motor cortex signals both initial and corrective movements during a precision reaching task. We recorded a large population of neurons from two male rhesus macaques across multiple sessions to examine the neural firing rates during not only initial movements but also subsequent corrective movements. AutoLFADS, an autoencoder-based deep-learning model, was applied to provide a clearer picture of neurons' activity on individual corrective movements across sessions. Decoding of reach velocity generalized poorly from initial to corrective submovements. Unlike initial movements, it was challenging to predict the velocity of corrective movements using traditional linear methods in a single, global neural space. We identified several locations in the neural space where corrective submovements originated after the initial reaches, signifying firing rates different than the baseline before initial movements. To improve corrective movement decoding, we demonstrate that a state-dependent decoder incorporating the population firing rates at the initiation of correction improved performance, highlighting the diverse neural features of corrective movements. In summary, we show neural differences between initial and corrective submovements and how the neural activity encodes specific combinations of velocity and position. These findings are inconsistent with assumptions that neural correlations with kinematic features are global and independent, emphasizing that traditional methods often fall short in describing these diverse neural processes for online corrective movements.

Making Connections: Neurodevelopmental Changes in Brain Connectivity After Adverse Experiences in Early Adolescence.

Adverse childhood experiences have been linked to detrimental mental health outcomes in adulthood. This study investigates a potential neurodevelopmental pathway between adversity and mental health outcomes: brain connectivity. We used data from the prospective, longitudinal Adolescent Brain Cognitive Development (ABCD) study (N ≍ 12.000, participants aged 9-13 years, male and female) and assessed structural brain connectivity using fractional anisotropy (FA) of white matter tracts. The adverse experiences modeled included family conflict and traumatic experiences. K-means clustering and latent basis growth models were used to determine subgroups based on total levels and trajectories of brain connectivity. Multinomial regression was used to determine associations between cluster membership and adverse experiences. The results showed that higher family conflict was associated with higher FA levels across brain tracts (e.g., t (3) = -3.81, ß = -0.09, p bonf = 0.003) and within the corpus callosum (CC), fornix, and anterior thalamic radiations (ATR). A decreasing FA trajectory across two brain imaging timepoints was linked to lower socioeconomic status and neighborhood safety. Socioeconomic status was related to FA across brain tracts (e.g., t (3) = 3.44, ß = 0.10, p bonf = 0.01), the CC and the ATR. Neighborhood safety was associated with FA in the Fornix and ATR (e.g., t (1) = 3.48, ß = 0.09, p bonf = 0.01). There is a complex and multifaceted relationship between adverse experiences and brain development, where adverse experiences during early adolescence are related to brain connectivity. These findings underscore the importance of studying adverse experiences beyond early childhood to understand lifespan developmental outcomes.

The ß8 integrin cytoplasmic domain activates extracellular matrix adhesion to promote brain neurovascular development.

In the developing mammalian brain, neuroepithelial cells interact with blood vessels to regulate angiogenesis, blood-brain barrier maturation and other key neurovascular functions. Genetic studies in mice have shown that neurovascular development is controlled, in part, by Itgb8, which encodes the neuroepithelial cell-expressed integrin ß8 subunit. However, these studies have involved complete loss-of-function Itgb8 mutations, and have not discerned the relative roles for the ß8 integrin extracellular matrix (ECM) binding region versus the intracellular signaling tail. Here, Cre/lox strategies have been employed to selectively delete the cytoplasmic tail of murine Itgb8 without perturbing its transmembrane and extracellular domains. We report that the ß8 integrin cytoplasmic domain is essential for inside-out modulation of adhesion, including activation of latent-TGFßs in the ECM. Quantitative sequencing of the brain endothelial cell transcriptome identifies TGFß-regulated genes with putative links to blood vessel morphogenesis, including several genes linked to Wnt/ß-catenin signaling. These results reveal that the ß8 integrin cytoplasmic domain is essential for the regulation of TGFß-dependent gene expression in endothelial cells and suggest that cross-talk between TGFßs and Wnt pathways is crucial for neurovascular development.

Identification of complier and noncomplier average causal effects in the presence of latent missing-at-random (LMAR) outcomes: a unifying view and choices of assumptions.

The study of treatment effects is often complicated by noncompliance and missing data. In the one-sided noncompliance setting where of interest are the complier and noncomplier average causal effects, we address outcome missingness of the latent missing at random type (LMAR, also known as latent ignorability). That is, conditional on covariates and treatment assigned, the missingness may depend on compliance type. Within the instrumental variable (IV) approach to noncompliance, methods have been proposed for handling LMAR outcome that additionally invoke an exclusion restriction-type assumption on missingness, but no solution has been proposed for when a non-IV approach is used. This article focuses on effect identification in the presence of LMAR outcomes, with a view to flexibly accommodate different principal identification approaches. We show that under treatment assignment ignorability and LMAR only, effect nonidentifiability boils down to a set of two connected mixture equations involving unidentified stratum-specific response probabilities and outcome means. This clarifies that (except for a special case) effect identification generally requires two additional assumptions: a specific missingness mechanism assumption and a principal identification assumption. This provides a template for identifying effects based on separate choices of these assumptions. We consider a range of specific missingness assumptions, including those that have appeared in the literature and some new ones. Incidentally, we find an issue in the existing assumptions, and propose a modification of the assumptions to avoid the issue. Results under different assumptions are illustrated using data from the Baltimore Experience Corps Trial.

Direct estimation and inference of higher-level correlations from lower-level measurements with applications in gene-pathway and proteomics studies.

This paper tackles the challenge of estimating correlations between higher-level biological variables (e.g. proteins and gene pathways) when only lower-level measurements are directly observed (e.g. peptides and individual genes). Existing methods typically aggregate lower-level data into higher-level variables and then estimate correlations based on the aggregated data. However, different data aggregation methods can yield varying correlation estimates as they target different higher-level quantities. Our solution is a latent factor model that directly estimates these higher-level correlations from lower-level data without the need for data aggregation. We further introduce a shrinkage estimator to ensure the positive definiteness and improve the accuracy of the estimated correlation matrix. Furthermore, we establish the asymptotic normality of our estimator, enabling efficient computation of P-values for the identification of significant correlations. The effectiveness of our approach is demonstrated through comprehensive simulations and the analysis of proteomics and gene expression datasets. We develop the R package highcor for implementing our method.

Tree-informed Bayesian multi-source domain adaptation: cross-population probabilistic cause-of-death assignment using verbal autopsy.

Determining causes of deaths (CODs) occurred outside of civil registration and vital statistics systems is challenging. A technique called verbal autopsy (VA) is widely adopted to gather information on deaths in practice. A VA consists of interviewing relatives of a deceased person about symptoms of the deceased in the period leading to the death, often resulting in multivariate binary responses. While statistical methods have been devised for estimating the cause-specific mortality fractions (CSMFs) for a study population, continued expansion of VA to new populations (or "domains") necessitates approaches that recognize between-domain differences while capitalizing on potential similarities. In this article, we propose such a domain-adaptive method that integrates external between-domain similarity information encoded by a prespecified rooted weighted tree. Given a cause, we use latent class models to characterize the conditional distributions of the responses that may vary by domain. We specify a logistic stick-breaking Gaussian diffusion process prior along the tree for class mixing weights with node-specific spike-and-slab priors to pool information between the domains in a data-driven way. The posterior inference is conducted via a scalable variational Bayes algorithm. Simulation studies show that the domain adaptation enabled by the proposed method improves CSMF estimation and individual COD assignment. We also illustrate and evaluate the method using a validation dataset. The article concludes with a discussion of limitations and future directions.

Mendelian randomization analysis using multiple biomarkers of an underlying common exposure.

Mendelian randomization (MR) analysis is increasingly popular for testing the causal effect of exposures on disease outcomes using data from genome-wide association studies. In some settings, the underlying exposure, such as systematic inflammation, may not be directly observable, but measurements can be available on multiple biomarkers or other types of traits that are co-regulated by the exposure. We propose a method for MR analysis on latent exposures (MRLE), which tests the significance for, and the direction of, the effect of a latent exposure by leveraging information from multiple related traits. The method is developed by constructing a set of estimating functions based on the second-order moments of GWAS summary association statistics for the observable traits, under a structural equation model where genetic variants are assumed to have indirect effects through the latent exposure and potentially direct effects on the traits. Simulation studies show that MRLE has well-controlled type I error rates and enhanced power compared to single-trait MR tests under various types of pleiotropy. Applications of MRLE using genetic association statistics across five inflammatory biomarkers (CRP, IL-6, IL-8, TNF-α, and MCP-1) provide evidence for potential causal effects of inflammation on increasing the risk of coronary artery disease, colorectal cancer, and rheumatoid arthritis, while standard MR analysis for individual biomarkers fails to detect consistent evidence for such effects.

DifferentialRegulation: a Bayesian hierarchical approach to identify differentially regulated genes.

Although transcriptomics data is typically used to analyze mature spliced mRNA, recent attention has focused on jointly investigating spliced and unspliced (or precursor-) mRNA, which can be used to study gene regulation and changes in gene expression production. Nonetheless, most methods for spliced/unspliced inference (such as RNA velocity tools) focus on individual samples, and rarely allow comparisons between groups of samples (e.g. healthy vs. diseased). Furthermore, this kind of inference is challenging, because spliced and unspliced mRNA abundance is characterized by a high degree of quantification uncertainty, due to the prevalence of multi-mapping reads, ie reads compatible with multiple transcripts (or genes), and/or with both their spliced and unspliced versions. Here, we present DifferentialRegulation, a Bayesian hierarchical method to discover changes between experimental conditions with respect to the relative abundance of unspliced mRNA (over the total mRNA). We model the quantification uncertainty via a latent variable approach, where reads are allocated to their gene/transcript of origin, and to the respective splice version. We designed several benchmarks where our approach shows good performance, in terms of sensitivity and error control, vs. state-of-the-art competitors. Importantly, our tool is flexible, and works with both bulk and single-cell RNA-sequencing data. DifferentialRegulation is distributed as a Bioconductor R package.

Multi-view spectral clustering with latent representation learning for applications on multi-omics cancer subtyping.

Driven by multi-omics data, some multi-view clustering algorithms have been successfully applied to cancer subtypes prediction, aiming to identify subtypes with biometric differences in the same cancer, thereby improving the clinical prognosis of patients and designing personalized treatment plan. Due to the fact that the number of patients in omics data is much smaller than the number of genes, multi-view spectral clustering based on similarity learning has been widely developed. However, these algorithms still suffer some problems, such as over-reliance on the quality of pre-defined similarity matrices for clustering results, inability to reasonably handle noise and redundant information in high-dimensional omics data, ignoring complementary information between omics data, etc. This paper proposes multi-view spectral clustering with latent representation learning (MSCLRL) method to alleviate the above problems. First, MSCLRL generates a corresponding low-dimensional latent representation for each omics data, which can effectively retain the unique information of each omics and improve the robustness and accuracy of the similarity matrix. Second, the obtained latent representations are assigned appropriate weights by MSCLRL, and global similarity learning is performed to generate an integrated similarity matrix. Third, the integrated similarity matrix is used to feed back and update the low-dimensional representation of each omics. Finally, the final integrated similarity matrix is used for clustering. In 10 benchmark multi-omics datasets and 2 separate cancer case studies, the experiments confirmed that the proposed method obtained statistically and biologically meaningful cancer subtypes.