Dictionary learning for integrative, multimodal and scalable single-cell analysis.

Mapping single-cell sequencing profiles to comprehensive reference datasets provides a powerful alternative to unsupervised analysis. However, most reference datasets are constructed from single-cell RNA-sequencing data and cannot be used to annotate datasets that do not measure gene expression. Here we introduce 'bridge integration', a method to integrate single-cell datasets across modalities using a multiomic dataset as a molecular bridge. Each cell in the multiomic dataset constitutes an element in a 'dictionary', which is used to reconstruct unimodal datasets and transform them into a shared space. Our procedure accurately integrates transcriptomic data with independent single-cell measurements of chromatin accessibility, histone modifications, DNA methylation and protein levels. Moreover, we demonstrate how dictionary learning can be combined with sketching techniques to improve computational scalability and harmonize 8.6 million human immune cell profiles from sequencing and mass cytometry experiments. Our approach, implemented in version 5 of our Seurat toolkit ( http://www.satijalab.org/seurat ), broadens the utility of single-cell reference datasets and facilitates comparisons across diverse molecular modalities.

Two-dimensional materials as sodium-ion battery anodes: The mass transfer and storage mechanisms of "fat" Na.

Sodium-ion batteries (SIBs) with abundant resource and high safety are attracting intensive interest from both research and industry communities in meeting the ever-increasing energy demands. Despite the rapid advance of SIBs, it is difficult yet necessary to enhance the cycling and rate performance at anode due to the sluggish kinetics of "fat" Na+. This review provides an overview of two-dimensional (2D) nanomaterials with a short ion diffusion pathway and a superior active sites exposure from the perspectives of synthesis, material chemistry, and structure engineering. We present the design principle of ideal carbon materials in SIBs. Moreover, we discuss the structure and chemistry regulations of different 2D materials to promote the efficient ion mass transfer and storage according to the different mechanisms of alloying, conversion, and insertion. Finally, we propose the remaining challenges and the possible solutions, in hope of guiding the future development of this booming field.

A Universal Method for Regulating Carbon Microcrystalline Structure for High-Capacity Sodium Storage: Binding Energy As Descriptor.

Sodium-ion batteries (SIBs) are attracting worldwide attention due to their multiple merits including abundant reserve and safety. However, industrialization is challenged by the scarcity of high-performance carbon anodes with high specific capacities. Here, we report the metal-assisted microcrystalline structure regulation of carbon materials to achieve high-capacity sodium storage. Systematic investigations of in situ thermal-treatment X-ray diffraction and multiple spectroscopies uncover the regulation mechanism of constructing steric hindrance (C-O-C bonds) to restrain the aromatic polycondensation reaction. The carbon precursor of polycyclic aromatic hydrocarbon-type pitch contributes to a high carbon yield rate (40%) compared with those of resin and biomass precursors. The as-synthesized carbon materials deliver high capacities of up to 390 mAh g-1, surpassing many reported carbon anodes for SIBs. Through correlating specific capacity with ID/IG values in Raman spectra and theoretical calculation of carbon materials regulated by different metal elements (Mn, Nb, Ce, Cr, and V), we identify and propose the binding energy as the descriptor for characterizing the capability of regulating the carbon microcrystalline structure to promote sodium storage. This work provides a universal method for regulating the carbon structure, which may lead to the controlled design and fabrication of carbon materials for energy storage and conversion and beyond.

Chemical and spatial dual-confinement engineering for stable Na-S batteries with approximately 100% capacity retention.

Sodium-sulfur (Na-S) batteries are attracting intensive attention due to the merits like high energy and low cost, while the poor stability of sulfur cathode limits the further development. Here, we report a chemical and spatial dual-confinement approach to improve the stability of Na-S batteries. It refers to covalently bond sulfur to carbon at forms of C-S/N-C=S bonds with high strength for locking sulfur. Meanwhile, sulfur is examined to be S1-S2 small species produced by thermally cutting S8 large molecules followed by sealing in the confined pores of carbon materials. Hence, the sulfur cathode achieves a good stability of maintaining a high-capacity retention of 97.64% after 1000 cycles. Experimental and theoretical results show that Na+ is hosted via a coordination structure (N···Na···S) without breaking the C-S bond, thus impeding the formation and dissolution of sodium polysulfide to ensure a good cycling stability. This work provides a promising method for addressing the S-triggered stability problem of Na-S batteries and other S-based batteries.

Multimodal single-cell datasets characterize antigen-specific CD8+ T cells across SARS-CoV-2 vaccination and infection.

The immune response to SARS-CoV-2 antigen after infection or vaccination is defined by the durable production of antibodies and T cells. Population-based monitoring typically focuses on antibody titer, but there is a need for improved characterization and quantification of T cell responses. Here, we used multimodal sequencing technologies to perform a longitudinal analysis of circulating human leukocytes collected before and after immunization with the mRNA vaccine BNT162b2. Our data indicated distinct subpopulations of CD8+ T cells, which reliably appeared 28 days after prime vaccination. Using a suite of cross-modality integration tools, we defined their transcriptome, accessible chromatin landscape and immunophenotype, and we identified unique biomarkers within each modality. We further showed that this vaccine-induced population was SARS-CoV-2 antigen-specific and capable of rapid clonal expansion. Moreover, we identified these CD8+ T cell populations in scRNA-seq datasets from COVID-19 patients and found that their relative frequency and differentiation outcomes were predictive of subsequent clinical outcomes.

CPA-Perturb-seq: Multiplexed single-cell characterization of alternative polyadenylation regulators.

Most mammalian genes have multiple polyA sites, representing a substantial source of transcript diversity that is governed by the cleavage and polyadenylation (CPA) regulatory machinery. To better understand how these proteins govern polyA site choice we introduce CPA-Perturb-seq, a multiplexed perturbation screen dataset of 42 known CPA regulators with a 3' scRNA-seq readout that enables transcriptome-wide inference of polyA site usage. We develop a statistical framework to specifically identify perturbation-dependent changes in intronic and tandem polyadenylation, and discover modules of co-regulated polyA sites exhibiting distinct functional properties. By training a multi-task deep neural network (APARENT-Perturb) on our dataset, we delineate a cis-regulatory code that predicts responsiveness to perturbation and reveals interactions between distinct regulatory complexes. Finally, we leverage our framework to re-analyze published scRNA-seq datasets, identifying new regulators that affect the relative abundance of alternatively polyadenylated transcripts, and characterizing extensive cellular heterogeneity in 3' UTR length amongst antibody-producing cells. Our work highlights the potential for multiplexed single-cell perturbation screens to further our understanding of post-transcriptional regulation in vitro and in vivo.

Green means go: Green light promotes hypocotyl elongation via brassinosteroid signaling.

Although many studies have elucidated the mechanisms by which different wavelengths of light (blue, red, far-red, or ultraviolet-B [UV-B]) regulate plant development, whether and how green light regulates plant development remains largely unknown. Previous studies reported that green light participates in regulating growth and development in land plants, but these studies have reported conflicting results, likely due to technical problems. For example, commercial green light-emitting diode light sources emit a little blue or red light. Here, using a pure green light source, we determined that unlike blue, red, far-red, or UV-B light, which inhibits hypocotyl elongation, green light promotes hypocotyl elongation in Arabidopsis thaliana and several other plants during the first 2-3 d after planting. Phytochromes, cryptochromes, and other known photoreceptors do not mediate green-light-promoted hypocotyl elongation, but the brassinosteroid (BR) signaling pathway is involved in this process. Green light promotes the DNA binding activity of BRI1-EMS-SUPPRESSOR 1 (BES1), a master transcription factor of the BR pathway, thus regulating gene transcription to promote hypocotyl elongation. Our results indicate that pure green light promotes elongation via BR signaling and acts as a shade signal to enable plants to adapt their development to a green-light-dominant environment under a canopy.

Multimodal characterization of antigen-specific CD8 + T cells across SARS-CoV-2 vaccination and infection.

The human immune response to SARS-CoV-2 antigen after infection or vaccination is defined by the durable production of antibodies and T cells. Population-based monitoring typically focuses on antibody titer, but there is a need for improved characterization and quantification of T cell responses. Here, we utilize multimodal sequencing technologies to perform a longitudinal analysis of circulating human leukocytes collected before and after BNT162b2 immunization. Our data reveal distinct subpopulations of CD8 + T cells which reliably appear 28 days after prime vaccination (7 days post boost). Using a suite of cross-modality integration tools, we define their transcriptome, accessible chromatin landscape, and immunophenotype, and identify unique biomarkers within each modality. By leveraging DNA-oligo-tagged peptide-MHC multimers and T cell receptor sequencing, we demonstrate that this vaccine-induced population is SARS-CoV-2 antigen-specific and capable of rapid clonal expansion. Moreover, we also identify these CD8 + populations in scRNA-seq datasets from COVID-19 patients and find that their relative frequency and differentiation outcomes are predictive of subsequent clinical outcomes. Our work contributes to our understanding of T cell immunity, and highlights the potential for integrative and multimodal analysis to characterize rare cell populations.

Efficient combinatorial targeting of RNA transcripts in single cells with Cas13 RNA Perturb-seq.

Pooled CRISPR screens coupled with single-cell RNA-sequencing have enabled systematic interrogation of gene function and regulatory networks. Here, we introduce Cas13 RNA Perturb-seq (CaRPool-seq), which leverages the RNA-targeting CRISPR-Cas13d system and enables efficient combinatorial perturbations alongside multimodal single-cell profiling. CaRPool-seq encodes multiple perturbations on a cleavable CRISPR array that is associated with a detectable barcode sequence, allowing for the simultaneous targeting of multiple genes. We compared CaRPool-seq to existing Cas9-based methods, highlighting its unique strength to efficiently profile combinatorially perturbed cells. Finally, we apply CaRPool-seq to perform multiplexed combinatorial perturbations of myeloid differentiation regulators in an acute myeloid leukemia (AML) model system and identify extensive interactions between different chromatin regulators that can enhance or suppress AML differentiation phenotypes.

Oral Microbiome in Nonsmoker Patients with Oral Cavity Squamous Cell Carcinoma, Defined by Metagenomic Shotgun Sequencing.

Objectives: Smoking is the commonest cause of oral cavity squamous cell carcinoma (OC-SCC), but the etiology of OC-SCC in nonsmokers is unknown. Our primary goal was to use metagenomic shotgun sequencing (MSS) to define the taxonomic composition and functional potential of oral metagenome in nonsmokers with OC-SCC. Methods: We conducted a case-control study with 42 OC-SCC case and 45 control nonsmokers. MSS was performed on DNA extracted from mouthwash samples. Taxonomic analysis and pathway analysis were done using MetaPhlAn2 and HUMAnN2, respectively. Statistical difference was determined using the Mann-Whitney test controlling false discovery rate. Results: There was no significant difference in age, sex, race, or alcohol consumption between OC-SCC and control patients. There was a significant difference in beta diversity between OC-SCC and controls. At the phylum level, Bacteroidetes and Synergistetes were overly represented in OC-SCC while Actinobacteria and Firmicutes were overly represented in controls. At the genus level, Fusobacterium was overly represented in OC-SCC compared with controls, while Corynebacterium, Streptococcus, Actinomyces, Cryptobacterium, and Selenomonas were overly represented in controls. Bacterial pathway analysis identified overrepresentation in OC-SCC of pathways related to metabolism of flavin, biotin, thiamin, heme, sugars, fatty acids, peptidoglycans, and tRNA and overrepresentation of nucleotides and essential amino acids in controls. Conclusions: The oral microbiome in nonsmoker patients with OC-SCC is significantly different from that of nonsmoker control patients in taxonomic compositions and functional potentials. Our study's MSS findings matched with previous 16S-based methods in taxonomic differentiation but varied greatly in functional differentiation of microbiomes in OC-SCC and controls.

A RORγt+ cell instructs gut microbiota-specific Treg cell differentiation.

The mutualistic relationship of gut-resident microbiota and the host immune system promotes homeostasis that ensures maintenance of the microbial community and of a largely non-aggressive immune cell compartment1,2. The consequences of disturbing this balance include proximal inflammatory conditions, such as Crohn's disease, and systemic illnesses. This equilibrium is achieved in part through the induction of both effector and suppressor arms of the adaptive immune system. Helicobacter species induce T regulatory (Treg) and T follicular helper (TFH) cells under homeostatic conditions, but induce inflammatory T helper 17 (TH17) cells when induced Treg (iTreg) cells are compromised3,4. How Helicobacter and other gut bacteria direct T cells to adopt distinct functions remains poorly understood. Here we investigated the cells and molecular components required for iTreg cell differentiation. We found that antigen presentation by cells expressing RORγt, rather than by classical dendritic cells, was required and sufficient for induction of Treg cells. These RORγt+ cells-probably type 3 innate lymphoid cells and/or Janus cells5-require the antigen-presentation machinery, the chemokine receptor CCR7 and the TGFß activator αv integrin. In the absence of any of these factors, there was expansion of pathogenic TH17 cells instead of iTreg cells, induced by CCR7-independent antigen-presenting cells. Thus, intestinal commensal microbes and their products target multiple antigen-presenting cells with pre-determined features suited to directing appropriate T cell differentiation programmes, rather than a common antigen-presenting cell that they endow with appropriate functions.

Progressive dysbiosis of human orodigestive microbiota along the sequence of gastroesophageal reflux, Barrett's esophagus and esophageal adenocarcinoma.

The incidence of esophageal adenocarcinoma (EA) has drastically increased in the United States since 1970s for unclear reasons. We hypothesized that the widespread usage of antibiotics has increased the procarcinogenic potential of the orodigestive microbiota along the sequence of gastroesophageal reflux (GR), Barrett's esophagus (BE) and EA phenotypes. This case control study included normal controls (NC) and three disease phenotypes GR, BE and EA. Microbiota in the mouth, esophagus, and stomach, and rectum were analyzed using 16S rRNA gene sequencing. Overall, we discovered 44 significant pairwise differences in abundance of microbial taxa between the four phenotypes, with 12 differences in the mouth, 21 in the esophagus, two in the stomach, and nine in the rectum. Along the GR→BE→EA sequence, oral and esophageal microbiota were more diversified, the dominant genus Streptococcus was progressively depleted while six other genera Atopobium, Actinomyces, Veillonella, Ralstonia, Burkholderia and Lautropia progressively enriched. In NC, Streptococcus appeared to control populations of other genera in the foregut via numerous negative and positive connections, while in disease states, the rich network was markedly simplified. Inferred gene functional content showed a progressive enrichment through the stages of EA development in genes encoding antibiotic resistance, ligands of Toll-like and NOD-like receptors, nitrate-nitrite-nitric oxide pathway and acetaldehyde metabolism. The orodigestive microbiota is in a progressive dysbiotic state along the GR-BE-EA sequence. The increasing dysbiosis and antibiotic and procarcinogenic genes in the disease states warrants further study to define their roles in EA pathogenesis.

Characterizing cellular heterogeneity in chromatin state with scCUT&Tag-pro.

Technologies that profile chromatin modifications at single-cell resolution offer enormous promise for functional genomic characterization, but the sparsity of the measurements and integrating multiple binding maps represent substantial challenges. Here we introduce single-cell (sc)CUT&Tag-pro, a multimodal assay for profiling protein-DNA interactions coupled with the abundance of surface proteins in single cells. In addition, we introduce single-cell ChromHMM, which integrates data from multiple experiments to infer and annotate chromatin states based on combinatorial histone modification patterns. We apply these tools to perform an integrated analysis across nine different molecular modalities in circulating human immune cells. We demonstrate how these two approaches can characterize dynamic changes in the function of individual genomic elements across both discrete cell states and continuous developmental trajectories, nominate associated motifs and regulators that establish chromatin states and identify extensive and cell-type-specific regulatory priming. Finally, we demonstrate how our integrated reference can serve as a scaffold to map and improve the interpretation of additional scCUT&Tag datasets.

Oral and gastric microbiome in relation to gastric intestinal metaplasia.

Evidence suggests that Helicobacter pylori plays a role in gastric cancer (GC) initiation. However, epidemiologic studies on the specific role of other bacteria in the development of GC are lacking. We conducted a case-control study of 89 cases with gastric intestinal metaplasia (IM) and 89 matched controls who underwent upper gastrointestinal endoscopy at three sites affiliated with NYU Langone Health. We performed shotgun metagenomic sequencing using oral wash samples from 89 case-control pairs and antral mucosal brushing samples from 55 case-control pairs. We examined the associations of relative abundances of bacterial taxa and functional pathways with IM using conditional logistic regression with and without elastic-net penalty. Compared with controls, oral species Peptostreptococcus stomatis, Johnsonella ignava, Neisseria elongata and Neisseria flavescens were enriched in cases (odds ratios [ORs] = 1.29-1.50, P = .004-.01) while Lactobacillus gasseri, Streptococcus mutans, S parasanguinis and S sanguinis were under-represented (ORs = 0.66-0.76, P = .006-.042) in cases. Species J ignava and Filifactor alocis in the gastric microbiota were enriched (ORs = 3.27 and 1.43, P = .005 and .035, respectively), while S mutans, S parasanguinis and S sanguinis were under-represented (ORs = 0.61-0.75, P = .024-.046), in cases compared with controls. The lipopolysaccharide and ubiquinol biosynthesis pathways were more abundant in IM, while the sugar degradation pathways were under-represented in IM. The findings suggest potential roles of certain oral and gastric microbiota, which are correlated with regulation of pathways associated with inflammation, in the development of gastric precancerous lesions.

Integrated analysis of multimodal single-cell data.

The simultaneous measurement of multiple modalities represents an exciting frontier for single-cell genomics and necessitates computational methods that can define cellular states based on multimodal data. Here, we introduce "weighted-nearest neighbor" analysis, an unsupervised framework to learn the relative utility of each data type in each cell, enabling an integrative analysis of multiple modalities. We apply our procedure to a CITE-seq dataset of 211,000 human peripheral blood mononuclear cells (PBMCs) with panels extending to 228 antibodies to construct a multimodal reference atlas of the circulating immune system. Multimodal analysis substantially improves our ability to resolve cell states, allowing us to identify and validate previously unreported lymphoid subpopulations. Moreover, we demonstrate how to leverage this reference to rapidly map new datasets and to interpret immune responses to vaccination and coronavirus disease 2019 (COVID-19). Our approach represents a broadly applicable strategy to analyze single-cell multimodal datasets and to look beyond the transcriptome toward a unified and multimodal definition of cellular identity.

Scalable, multimodal profiling of chromatin accessibility, gene expression and protein levels in single cells.

Recent technological advances have enabled massively parallel chromatin profiling with scATAC-seq (single-cell assay for transposase accessible chromatin by sequencing). Here we present ATAC with select antigen profiling by sequencing (ASAP-seq), a tool to simultaneously profile accessible chromatin and protein levels. Our approach pairs sparse scATAC-seq data with robust detection of hundreds of cell surface and intracellular protein markers and optional capture of mitochondrial DNA for clonal tracking, capturing three distinct modalities in single cells. ASAP-seq uses a bridging approach that repurposes antibody:oligonucleotide conjugates designed for existing technologies that pair protein measurements with single-cell RNA sequencing. Together with DOGMA-seq, an adaptation of CITE-seq (cellular indexing of transcriptomes and epitopes by sequencing) for measuring gene activity across the central dogma of gene regulation, we demonstrate the utility of systematic multi-omic profiling by revealing coordinated and distinct changes in chromatin, RNA and surface proteins during native hematopoietic differentiation and peripheral blood mononuclear cell stimulation and as a combinatorial decoder and reporter of multiplexed perturbations in primary T cells.

Alterations in neural oscillations related to working memory deficit in temporal lobe epilepsy.

OBJECTIVE: Working memory (WM) is critical for higher level cognition, but the underlying neural mechanisms are not fully understood. Impaired WM affects routine daily activities and is observed in patients with temporal lobe epilepsy (TLE). This study investigated neural oscillations associated with different WM phases, to determine the specific neural activity linked with the phases of WM impairment. METHODS: Patients with TLE (n = 52) and healthy volunteers (n = 35) completed a WM task, during which 34-channel electroencephalogram signals were recorded. Characteristic neural oscillation patterns during each WM phase were compared between the 2 groups. RESULTS: Patients with TLE showed decreased theta power during the encoding phase of WM, which was associated with reduced accuracy in the WM task. Altered theta power in the frontal region of the brain during the encoding phase was associated with a longer reaction time. CONCLUSIONS: Alterations in theta oscillation are related to WM impairment in patients with TLE and may serve as an early marker for evaluating WM deficits. SIGNIFICANCE: This study provides an early marker for evaluating WM deficits in TLE.

Characterizing the molecular regulation of inhibitory immune checkpoints with multimodal single-cell screens.

The expression of inhibitory immune checkpoint molecules, such as programmed death-ligand (PD-L)1, is frequently observed in human cancers and can lead to the suppression of T cell-mediated immune responses. Here, we apply expanded CRISPR-compatible (EC)CITE-seq, a technology that combines pooled CRISPR screens with single-cell mRNA and surface protein measurements, to explore the molecular networks that regulate PD-L1 expression. We also develop a computational framework, mixscape, that substantially improves the signal-to-noise ratio in single-cell perturbation screens by identifying and removing confounding sources of variation. Applying these tools, we identify and validate regulators of PD-L1 and leverage our multimodal data to identify both transcriptional and post-transcriptional modes of regulation. Specifically, we discover that the Kelch-like protein KEAP1 and the transcriptional activator NRF2 mediate the upregulation of PD-L1 after interferon (IFN)-γ stimulation. Our results identify a new mechanism for the regulation of immune checkpoints and present a powerful analytical framework for the analysis of multimodal single-cell perturbation screens.

Case control study comparing the HPV genome in patients with oral cavity squamous cell carcinoma to normal patients using metagenomic shotgun sequencing.

The aim of this study was to carry out a case control study comparing the HPV genome in patients with oral cavity squamous cell carcinoma (OC-SCC) to normal patients using metagenomic shotgun sequencing. We recruited 50 OC-SCC cases which were then matched with a control patient by age, gender, race, smoking status and alcohol status. DNA was extracted from oral wash samples from all patients and whole genome shotgun sequencing performed. The raw sequence data was cleaned, reads aligned with the human genome (GRCH38), nonhuman reads identified and then HPV genotypes identified using HPViewer. In the 50 patients with OC-SCC, the most common subsite was tongue in 26 (52%). All patients were treated with primary resection and neck dissection. All but 2 tumors were negative on p16 immunohistochemistry. There were no statistically significant differences between the cases and controls in terms of gender, age, race/ethnicity, alcohol drinking, and cigarette smoking. There was no statistically significant difference between the cancer samples and control samples in the nonhuman DNA reads (medians 4,228,072 vs. 5,719,715, P value = 0.324). HPV was detected in 5 cases (10%) of OC-SCC (genotypes 10, 16, 98) but only 1 tumor sample (genotype 16) yielded a high number of reads to suggest a role in the etiology of OC-SCC. HPV was detected in 4 control patients (genotypes 16, 22, 76, 200) but all had only 1-2 HPV reads per human genome. Genotypes of HPV are rarely found in patients with oral cancer.