Bacterial droplet-based single-cell RNA-seq reveals antibiotic-associated heterogeneous cellular states.

We introduce BacDrop, a highly scalable technology for bacterial single-cell RNA sequencing that has overcome many challenges hindering the development of scRNA-seq in bacteria. BacDrop can be applied to thousands to millions of cells from both gram-negative and gram-positive species. It features universal ribosomal RNA depletion and combinatorial barcodes that enable multiplexing and massively parallel sequencing. We applied BacDrop to study Klebsiella pneumoniae clinical isolates and to elucidate their heterogeneous responses to antibiotic stress. In an unperturbed population presumed to be homogeneous, we found within-population heterogeneity largely driven by the expression of mobile genetic elements that promote the evolution of antibiotic resistance. Under antibiotic perturbation, BacDrop revealed transcriptionally distinct subpopulations associated with different phenotypic outcomes including antibiotic persistence. BacDrop thus can capture cellular states that cannot be detected by bulk RNA-seq, which will unlock new microbiological insights into bacterial responses to perturbations and larger bacterial communities such as the microbiome.

A transcription factor atlas of directed differentiation.

Transcription factors (TFs) regulate gene programs, thereby controlling diverse cellular processes and cell states. To comprehensively understand TFs and the programs they control, we created a barcoded library of all annotated human TF splice isoforms (>3,500) and applied it to build a TF Atlas charting expression profiles of human embryonic stem cells (hESCs) overexpressing each TF at single-cell resolution. We mapped TF-induced expression profiles to reference cell types and validated candidate TFs for generation of diverse cell types, spanning all three germ layers and trophoblasts. Targeted screens with subsets of the library allowed us to create a tailored cellular disease model and integrate mRNA expression and chromatin accessibility data to identify downstream regulators. Finally, we characterized the effects of combinatorial TF overexpression by developing and validating a strategy for predicting combinations of TFs that produce target expression profiles matching reference cell types to accelerate cellular engineering efforts.

A human fetal lung cell atlas uncovers proximal-distal gradients of differentiation and key regulators of epithelial fates.

We present a multiomic cell atlas of human lung development that combines single-cell RNA and ATAC sequencing, high-throughput spatial transcriptomics, and single-cell imaging. Coupling single-cell methods with spatial analysis has allowed a comprehensive cellular survey of the epithelial, mesenchymal, endothelial, and erythrocyte/leukocyte compartments from 5-22 post-conception weeks. We identify previously uncharacterized cell states in all compartments. These include developmental-specific secretory progenitors and a subtype of neuroendocrine cell related to human small cell lung cancer. Our datasets are available through our web interface (https://lungcellatlas.org). To illustrate its general utility, we use our cell atlas to generate predictions about cell-cell signaling and transcription factor hierarchies which we rigorously test using organoid models.

Mapping transcriptomic vector fields of single cells.

Single-cell (sc)RNA-seq, together with RNA velocity and metabolic labeling, reveals cellular states and transitions at unprecedented resolution. Fully exploiting these data, however, requires kinetic models capable of unveiling governing regulatory functions. Here, we introduce an analytical framework dynamo (https://github.com/aristoteleo/dynamo-release), which infers absolute RNA velocity, reconstructs continuous vector fields that predict cell fates, employs differential geometry to extract underlying regulations, and ultimately predicts optimal reprogramming paths and perturbation outcomes. We highlight dynamo's power to overcome fundamental limitations of conventional splicing-based RNA velocity analyses to enable accurate velocity estimations on a metabolically labeled human hematopoiesis scRNA-seq dataset. Furthermore, differential geometry analyses reveal mechanisms driving early megakaryocyte appearance and elucidate asymmetrical regulation within the PU.1-GATA1 circuit. Leveraging the least-action-path method, dynamo accurately predicts drivers of numerous hematopoietic transitions. Finally, in silico perturbations predict cell-fate diversions induced by gene perturbations. Dynamo, thus, represents an important step in advancing quantitative and predictive theories of cell-state transitions.

Transient Osmotic Perturbation Causes Long-Term Alteration to the Gut Microbiota.

Osmotic diarrhea is a prevalent condition in humans caused by food intolerance, malabsorption, and widespread laxative use. Here, we assess the resilience of the gut ecosystem to osmotic perturbation at multiple length and timescales using mice as model hosts. Osmotic stress caused reproducible extinction of highly abundant taxa and expansion of less prevalent members in human and mouse microbiotas. Quantitative imaging revealed decimation of the mucus barrier during osmotic perturbation, followed by recovery. The immune system exhibited temporary changes in cytokine levels and a lasting IgG response against commensal bacteria. Increased osmolality prevented growth of commensal strains in vitro, revealing one mechanism contributing to extinction. Environmental availability of microbiota members mitigated extinction events, demonstrating how species reintroduction can affect community resilience. Our findings (1) demonstrate that even mild osmotic diarrhea can cause lasting changes to the microbiota and host and (2) lay the foundation for interventions that increase system-wide resilience.

Genome-scale exon perturbation screens uncover exons critical for cell fitness.

CRISPR-Cas technology has transformed functional genomics, yet understanding of how individual exons differentially shape cellular phenotypes remains limited. Here, we optimized and conducted massively parallel exon deletion and splice-site mutation screens in human cell lines to identify exons that regulate cellular fitness. Fitness-promoting exons are prevalent in essential and highly expressed genes and commonly overlap with protein domains and interaction interfaces. Conversely, fitness-suppressing exons are enriched in nonessential genes, exhibiting lower inclusion levels, and overlap with intrinsically disordered regions and disease-associated mutations. In-depth mechanistic investigation of the screen-hit TAF5 alternative exon-8 revealed that its inclusion is required for assembly of the TFIID general transcription initiation complex, thereby regulating global gene expression output. Collectively, our orthogonal exon perturbation screens established a comprehensive repository of phenotypically important exons and uncovered regulatory mechanisms governing cellular fitness and gene expression.

Memory effects of prior subculture may impact the quality of multiomic perturbation profiles.

Mass spectrometry-based omics technologies are increasingly used in perturbation studies to map drug effects to biological pathways by identifying significant molecular events. Significance is influenced by fold change and variation of each molecular parameter, but also by multiple testing corrections. While the fold change is largely determined by the biological system, the variation is determined by experimental workflows. Here, it is shown that memory effects of prior subculture can influence the variation of perturbation profiles using the two colon carcinoma cell lines SW480 and HCT116. These memory effects are largely driven by differences in growth states that persist into the perturbation experiment. In SW480 cells, memory effects combined with moderate treatment effects amplify the variation in multiple omics levels, including eicosadomics, proteomics, and phosphoproteomics. With stronger treatment effects, the memory effect was less pronounced, as demonstrated in HCT116 cells. Subculture homogeneity was controlled by real-time monitoring of cell growth. Controlled homogeneous subculture resulted in a perturbation network of 321 causal conjectures based on combined proteomic and phosphoproteomic data, compared to only 58 causal conjectures without controlling subculture homogeneity in SW480 cells. Some cellular responses and regulatory events were identified that extend the mode of action of arsenic trioxide (ATO) only when accounting for these memory effects. Controlled prior subculture led to the finding of a synergistic combination treatment of ATO with the thioredoxin reductase 1 inhibitor auranofin, which may prove useful in the management of NRF2-mediated resistance mechanisms.

Absolute energy levels of liquid water from many-body perturbation theory with effective vertex corrections.

We demonstrate the importance of addressing the Γ vertex and thus going beyond the GW approximation for achieving the energy levels of liquid water in many-body perturbation theory. In particular, we consider an effective vertex function in both the polarizability and the self-energy, which does not produce any computational overhead compared with the GW approximation. We yield the band gap, the ionization potential, and the electron affinity in good agreement with experiment and with a hybrid functional description. The achieved electronic structure and dielectric screening further lead to a good description of the optical absorption spectrum, as obtained through the solution of the Bethe-Salpeter equation. In particular, the experimental peak position of the exciton is accurately reproduced.

Applications of CRISPR technology in cellular immunotherapy.

CRISPR technology has transformed multiple fields, including cancer and immunology. CRISPR-based gene editing and screening empowers direct genomic manipulation of immune cells, opening doors to unbiased functional genetic screens. These screens aid in the discovery of novel factors that regulate and reprogram immune responses, offering novel drug targets. The engineering of immune cells using CRISPR has sparked a transformation in the cellular immunotherapy field, resulting in a multitude of ongoing clinical trials. In this review, we discuss the development and applications of CRISPR and related gene editing technologies in immune cells, focusing on functional genomics screening, gene editing-based cell therapies, as well as future directions in this rapidly advancing field.

CPADS: a web tool for comprehensive pancancer analysis of drug sensitivity.

Drug therapy is vital in cancer treatment. Accurate analysis of drug sensitivity for specific cancers can guide healthcare professionals in prescribing drugs, leading to improved patient survival and quality of life. However, there is a lack of web-based tools that offer comprehensive visualization and analysis of pancancer drug sensitivity. We gathered cancer drug sensitivity data from publicly available databases (GEO, TCGA and GDSC) and developed a web tool called Comprehensive Pancancer Analysis of Drug Sensitivity (CPADS) using Shiny. CPADS currently includes transcriptomic data from over 29 000 samples, encompassing 44 types of cancer, 288 drugs and more than 9000 gene perturbations. It allows easy execution of various analyses related to cancer drug sensitivity. With its large sample size and diverse drug range, CPADS offers a range of analysis methods, such as differential gene expression, gene correlation, pathway analysis, drug analysis and gene perturbation analysis. Additionally, it provides several visualization approaches. CPADS significantly aids physicians and researchers in exploring primary and secondary drug resistance at both gene and pathway levels. The integration of drug resistance and gene perturbation data also presents novel perspectives for identifying pivotal genes influencing drug resistance. Access CPADS at https://smuonco.shinyapps.io/CPADS/ or https://robinl-lab.com/CPADS.

scDecouple: decoupling cellular response from infected proportion bias in scCRISPR-seq.

Single-cell clustered regularly interspaced short palindromic repeats-sequencing (scCRISPR-seq) is an emerging high-throughput CRISPR screening technology where the true cellular response to perturbation is coupled with infected proportion bias of guide RNAs (gRNAs) across different cell clusters. The mixing of these effects introduces noise into scCRISPR-seq data analysis and thus obstacles to relevant studies. We developed scDecouple to decouple true cellular response of perturbation from the influence of infected proportion bias. scDecouple first models the distribution of gene expression profiles in perturbed cells and then iteratively finds the maximum likelihood of cell cluster proportions as well as the cellular response for each gRNA. We demonstrated its performance in a series of simulation experiments. By applying scDecouple to real scCRISPR-seq data, we found that scDecouple enhances the identification of biologically perturbation-related genes. scDecouple can benefit scCRISPR-seq data analysis, especially in the case of heterogeneous samples or complex gRNA libraries.

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.

Herb-CMap: a multimodal fusion framework for deciphering the mechanisms of action in traditional Chinese medicine using Suhuang antitussive capsule as a case study.

Herbal medicines, particularly traditional Chinese medicines (TCMs), are a rich source of natural products with significant therapeutic potential. However, understanding their mechanisms of action is challenging due to the complexity of their multi-ingredient compositions. We introduced Herb-CMap, a multimodal fusion framework leveraging protein-protein interactions and herb-perturbed gene expression signatures. Utilizing a network-based heat diffusion algorithm, Herb-CMap creates a connectivity map linking herb perturbations to their therapeutic targets, thereby facilitating the prioritization of active ingredients. As a case study, we applied Herb-CMap to Suhuang antitussive capsule (Suhuang), a TCM formula used for treating cough variant asthma (CVA). Using in vivo rat models, our analysis established the transcriptomic signatures of Suhuang and identified its key compounds, such as quercetin and luteolin, and their target genes, including IL17A, PIK3CB, PIK3CD, AKT1, and TNF. These drug-target interactions inhibit the IL-17 signaling pathway and deactivate PI3K, AKT, and NF-κB, effectively reducing lung inflammation and alleviating CVA. The study demonstrates the efficacy of Herb-CMap in elucidating the molecular mechanisms of herbal medicines, offering valuable insights for advancing drug discovery in TCM.

Dynamic extracellular vestibule of human SERT: Unveiling druggable potential with high-affinity allosteric inhibitors.

The serotonin transporter (SERT) tightly regulates synaptic serotonin levels and has been the primary target of antidepressants. Binding of inhibitors to the allosteric site of human SERT (hSERT) impedes the dissociation of antidepressants bound at the central site and may enhance the efficacy of such antidepressants to potentially reduce their dosage and side effects. Here, we report the identification of a series of high-affinity allosteric inhibitors of hSERT in a unique scaffold, with the lead compound, Lu AF88273 (3-(1-(2-(1H-indol-3-yl)ethyl)piperidin-4-yl)-6-chloro-1H-indole), having 2.1 nM allosteric potency in inhibiting imipramine dissociation. In addition, we find that Lu AF88273 also inhibits serotonin transport in a noncompetitive manner. The binding pose of Lu AF88273 in the allosteric site of hSERT is determined with extensive molecular dynamics simulations and rigorous absolute binding free energy perturbation (FEP) calculations, which show that a part of the compound occupies a dynamically formed small cavity. The predicted binding location and pose are validated by site-directed mutagenesis and can explain much of the structure-activity relationship of these inhibitors using the relative binding FEP calculations. Together, our findings provide a promising lead compound and the structural basis for the development of allosteric drugs targeting hSERT. Further, they demonstrate that the divergent allosteric sites of neurotransmitter transporters can be selectively targeted.

A mutagenesis study of autoantigen optimization for potential T1D vaccine design.

A previously reported autoreactive antigen, termed the X-idiotype, isolated from a unique cell population in Type 1 diabetes (T1D) patients, was found to stimulate their CD4+ T cells. This antigen was previously determined to bind more favorably than insulin and its mimic (insulin superagonist) to HLA-DQ8, supporting its strong role in CD4+ T cell activation. In this work, we probed HLA-X-idiotype-TCR binding and designed enhanced-reactive pHLA-TCR antigens using an in silico mutagenesis approach which we functionally validated by cell proliferation assays and flow cytometry. From a combination of single, double, and swap mutations, we identified antigen-binding sites p4 and p6 as potential mutation sites for HLA binding affinity enhancement. Site p6 is revealed to favor smaller but more hydrophobic residues than the native tyrosine, such as valine (Y6V) and isoleucine (Y6I), indicating a steric mechanism in binding affinity improvement. Meanwhile, site p4 methionine mutation to hydrophobic residues isoleucine (M4I) or leucine (M4L) modestly increases HLA binding affinity. Select p6 mutations to cysteine (Y6C) or isoleucine (Y6I) exhibit favorable TCR binding affinities, while a swap p5-p6 tyrosine-valine double mutant (V5Y_Y6V) and a p6-p7 glutamine-glutamine double mutant (Y6Q_Y7Q) exhibit enhanced HLA binding affinity but weakened TCR affinity. This work holds relevance to potential T1D antigen-based vaccine design and optimization.

Recurrent Neural Circuits Overcome Partial Inactivation by Compensation and Re-learning.

Technical advances in artificial manipulation of neural activity have precipitated a surge in studying the causal contribution of brain circuits to cognition and behavior. However, complexities of neural circuits challenge interpretation of experimental results, necessitating new theoretical frameworks for reasoning about causal effects. Here, we take a step in this direction, through the lens of recurrent neural networks trained to perform perceptual decisions. We show that understanding the dynamical system structure that underlies network solutions provides a precise account for the magnitude of behavioral effects due to perturbations. Our framework explains past empirical observations by clarifying the most sensitive features of behavior, and how complex circuits compensate and adapt to perturbations. In the process, we also identify strategies that can improve the interpretability of inactivation experiments.

Specific protonation of acidic residues confers K+ selectivity to the gastric proton pump.

The gastric proton pump (H+,K+-ATPase) transports a proton into the stomach lumen for every K+ ion exchanged in the opposite direction. In the lumen-facing state of the pump (E2), the pump selectively binds K+ despite the presence of a 10-fold higher concentration of Na+. The molecular basis for the ion selectivity of the pump is unknown. Using molecular dynamics simulations, free energy calculations, and Na+ and K+-dependent ATPase activity assays, we demonstrate that the K+ selectivity of the pump depends upon the simultaneous protonation of the acidic residues E343 and E795 in the ion-binding site. We also show that when E936 is protonated, the pump becomes Na+ sensitive. The protonation-mimetic mutant E936Q exhibits weak Na+-activated ATPase activity. A 2.5-Å resolution cryo-EM structure of the E936Q mutant in the K+-occluded E2-Pi form shows, however, no significant structural difference compared with wildtype except less-than-ideal coordination of K+ in the mutant. The selectivity toward a specific ion correlates with a more rigid and less fluctuating ion-binding site. Despite being exposed to a pH of 1, the fundamental principle driving the K+ ion selectivity of H+,K+-ATPase is similar to that of Na+,K+-ATPase: the ionization states of the acidic residues in the ion-binding sites determine ion selectivity. Unlike the Na+,K+-ATPase, however, protonation of an ion-binding glutamate residue (E936) confers Na+ sensitivity.

SCREE: a comprehensive pipeline for single-cell multi-modal CRISPR screen data processing and analysis.

Single-cell CRISPR screens have been widely used to investigate gene regulatory circuits in diverse biological systems. The recent development of single-cell CRISPR screens has enabled multimodal profiling of perturbed cells with both gene expression, chromatin accessibility and protein levels. However, current methods cannot meet the analysis requirements of different types of data and have limited functions. Here, we introduce Single-cell CRISPR screens data analysEs and perturbation modEling (SCREE) as a comprehensive and flexible pipeline to facilitate the analyses of various types of single-cell CRISPR screens data. SCREE performs read alignment, sgRNA assignment, quality control, clustering and visualization, perturbation enrichment evaluation, perturbation efficiency modeling, gene regulatory score calculation and functional analyses of perturbations for single-cell CRISPR screens with both RNA, ATAC and multimodal readout. SCREE is available at https://github.com/wanglabtongji/SCREE.

Comprehensive exploration of multi-modal and multi-branch imaging markers for autism diagnosis and interpretation: insights from an advanced deep learning model.

Autism spectrum disorder is a complex neurodevelopmental condition with diverse genetic and brain involvement. Despite magnetic resonance imaging advances, autism spectrum disorder diagnosis and understanding its neurogenetic factors remain challenging. We propose a dual-branch graph neural network that effectively extracts and fuses features from bimodalities, achieving 73.9% diagnostic accuracy. To explain the mechanism distinguishing autism spectrum disorder from healthy controls, we establish a perturbation model for brain imaging markers and perform a neuro-transcriptomic joint analysis using partial least squares regression and enrichment to identify potential genetic biomarkers. The perturbation model identifies brain imaging markers related to structural magnetic resonance imaging in the frontal, temporal, parietal, and occipital lobes, while functional magnetic resonance imaging markers primarily reside in the frontal, temporal, occipital lobes, and cerebellum. The neuro-transcriptomic joint analysis highlights genes associated with biological processes, such as "presynapse," "behavior," and "modulation of chemical synaptic transmission" in autism spectrum disorder's brain development. Different magnetic resonance imaging modalities offer complementary information for autism spectrum disorder diagnosis. Our dual-branch graph neural network achieves high accuracy and identifies abnormal brain regions and the neuro-transcriptomic analysis uncovers important genetic biomarkers. Overall, our study presents an effective approach for assisting in autism spectrum disorder diagnosis and identifying genetic biomarkers, showing potential for enhancing the diagnosis and treatment of this condition.

Proteomic Dynamics of Breast Cancer Cell Lines Identifies Potential Therapeutic Protein Targets.

Treatment and relevant targets for breast cancer (BC) remain limited, especially for triple-negative BC (TNBC). We identified 6091 proteins of 76 human BC cell lines using data-independent acquisition (DIA). Integrating our proteomic findings with prior multi-omics datasets, we found that including proteomics data improved drug sensitivity predictions and provided insights into the mechanisms of action. We subsequently profiled the proteomic changes in nine cell lines (five TNBC and four non-TNBC) treated with EGFR/AKT/mTOR inhibitors. In TNBC, metabolism pathways were dysregulated after EGFR/mTOR inhibitor treatment, while RNA modification and cell cycle pathways were affected by AKT inhibitor. This systematic multi-omics and in-depth analysis of the proteome of BC cells can help prioritize potential therapeutic targets and provide insights into adaptive resistance in TNBC.