Multiplexed single-cell characterization of alternative polyadenylation regulators.

Most mammalian genes have multiple polyA sites, representing a substantial source of transcript diversity regulated by the cleavage and polyadenylation (CPA) machinery. To better understand how these proteins govern polyA site choice, we introduce CPA-Perturb-seq, a multiplexed perturbation screen dataset of 42 CPA regulators with a 3' scRNA-seq readout that enables transcriptome-wide inference of polyA site usage. We develop a framework to detect perturbation-dependent changes in polyadenylation and characterize modules of co-regulated polyA sites. We find groups of intronic polyA sites regulated by distinct components of the nuclear RNA life cycle, including elongation, splicing, termination, and surveillance. We train and validate a deep neural network (APARENT-Perturb) for tandem polyA site usage, delineating a cis-regulatory code that predicts perturbation response and reveals interactions between regulatory complexes. Our work highlights the potential for multiplexed single-cell perturbation screens to further our understanding of post-transcriptional regulation.

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.

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.

Niche-Selective Inhibition of Pathogenic Th17 Cells by Targeting Metabolic Redundancy.

Targeting glycolysis has been considered therapeutically intractable owing to its essential housekeeping role. However, the context-dependent requirement for individual glycolytic steps has not been fully explored. We show that CRISPR-mediated targeting of glycolysis in T cells in mice results in global loss of Th17 cells, whereas deficiency of the glycolytic enzyme glucose phosphate isomerase (Gpi1) selectively eliminates inflammatory encephalitogenic and colitogenic Th17 cells, without substantially affecting homeostatic microbiota-specific Th17 cells. In homeostatic Th17 cells, partial blockade of glycolysis upon Gpi1 inactivation was compensated by pentose phosphate pathway flux and increased mitochondrial respiration. In contrast, inflammatory Th17 cells experience a hypoxic microenvironment known to limit mitochondrial respiration, which is incompatible with loss of Gpi1. Our study suggests that inhibiting glycolysis by targeting Gpi1 could be an effective therapeutic strategy with minimum toxicity for Th17-mediated autoimmune diseases, and, more generally, that metabolic redundancies can be exploited for selective targeting of disease processes.

Comprehensive Integration of Single-Cell Data.

Single-cell transcriptomics has transformed our ability to characterize cell states, but deep biological understanding requires more than a taxonomic listing of clusters. As new methods arise to measure distinct cellular modalities, a key analytical challenge is to integrate these datasets to better understand cellular identity and function. Here, we develop a strategy to "anchor" diverse datasets together, enabling us to integrate single-cell measurements not only across scRNA-seq technologies, but also across different modalities. After demonstrating improvement over existing methods for integrating scRNA-seq data, we anchor scRNA-seq experiments with scATAC-seq to explore chromatin differences in closely related interneuron subsets and project protein expression measurements onto a bone marrow atlas to characterize lymphocyte populations. Lastly, we harmonize in situ gene expression and scRNA-seq datasets, allowing transcriptome-wide imputation of spatial gene expression patterns. Our work presents a strategy for the assembly of harmonized references and transfer of information across datasets.

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.

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.

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.

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.

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.

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.

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.

Periodontal pathogens are a risk factor of oral cavity squamous cell carcinoma, independent of tobacco and alcohol and human papillomavirus.

Over the past decade, there has been a change in the epidemiology of oral cavity squamous cell cancer (OC-SCC). Many new cases of OC-SCC lack the recognized risk factors of smoking, alcohol and human papilloma virus. The aim of this study was to determine if the oral microbiome may be associated with OC-SCC in nonsmoking HPV negative patients. We compared the oral microbiome of HPV-negative nonsmoker OC-SCC(n = 18), premalignant lesions(PML) (n = 8) and normal control patients (n = 12). Their oral microbiome was sampled by oral wash and defined by 16S rRNA gene sequencing. We report that the periodontal pathogens Fusobacterium, Prevotella, Alloprevotella were enriched while commensal Streptococcus depleted in OC-SCC. Based on the four genera plus a marker genus Veillonella for PML, we classified the oral microbiome into two types. Gene/pathway analysis revealed a progressive increase of genes encoding HSP90 and ligands for TLRs 1, 2 and 4 along the controls→PML → OC-SCC progression sequence. Our findings suggest an association between periodontal pathogens and OC-SCC in non smoking HPV negative patients.

HPViewer: sensitive and specific genotyping of human papillomavirus in metagenomic DNA.

Motivation: Shotgun DNA sequencing provides sensitive detection of all 182 HPV types in tissue and body fluid. However, existing computational methods either produce false positives misidentifying HPV types due to shared sequences among HPV, human and prokaryotes, or produce false negative since they identify HPV by assembled contigs requiring large abundant of HPV reads. Results: We designed HPViewer with two custom HPV reference databases masking simple repeats and homology sequences respectively and one homology distance matrix to hybridize these two databases. It directly identified HPV from short DNA reads rather than assembled contigs. Using 100 100 simulated samples, we revealed that HPViewer was robust for samples containing either high or low number of HPV reads. Using 12 shotgun sequencing samples from respiratory papillomatosis, HPViewer was equal to VirusTAP, and Vipie and better than HPVDetector with the respect to specificity and was the most sensitive method in the detection of HPV types 6 and 11. We demonstrated that contigs-based approaches had disadvantages of detection of HPV. In 1573 sets of metagenomic data from 18 human body sites, HPViewer identified 104 types of HPV in a body-site associated pattern and 89 types of HPV co-occurring in one sample with other types of HPV. We demonstrated HPViewer was sensitive and specific for HPV detection in metagenomic data. Availability and implementation: HPViewer can be accessed at https://github.com/yuhanH/HPViewer/. Supplementary information: Supplementary data are available at Bioinformatics online.

Histone deacetylase HDA6 enhances brassinosteroid signaling by inhibiting the BIN2 kinase.

Glycogen synthase kinase 3 (GSK3)-like kinases play important roles in brassinosteroid (BR), abscisic acid, and auxin signaling to regulate many aspects of plant development and stress responses. The Arabidopsis thaliana GSK3-like kinase BR-INSENSITIVE 2 (BIN2) acts as a key negative regulator in the BR signaling pathway, but the mechanisms regulating BIN2 function remain unclear. Here we report that the histone deacetylase HDA6 can interact with and deacetylate BIN2 to repress its kinase activity. The hda6 mutant showed a BR-repressed phenotype in the dark and was less sensitive to BR biosynthesis inhibitors. Genetic analysis indicated that HDA6 regulates BR signaling through BIN2. Furthermore, we identified K189 of BIN2 as an acetylated site, which can be deacetylated by HDA6 to influence BIN2 activity. Glucose can affect the acetylation level of BIN2 in plants, indicating a connection to cellular energy status. These findings provide significant insights into the regulation of GSK3-like kinases in plant growth and development.

A subclass of HSP70s regulate development and abiotic stress responses in Arabidopsis thaliana.

Members of the HSP70 family function as molecular chaperones to maintain cellular homeostasis and help plants cope with environmental stimuli. However, due to functional redundancy and lack of effective chemical inhibitors, our knowledge of functions of individual HSP70s has remained limited. Here, we confirmed a subclass of HSP70s, including HSP70-1, -2, -3, -4, and -5, localized to the cytosol and nucleus in Arabidopsis thaliana. Histochemical analyses of promoter:GUS reporter lines showed that HSP70-1, -2, -3, and -4 genes were widely expressed, but HSP70-5 was not. In addition, individual HSP70 showed not only similar but also distinct transcriptions when treated by different abiotic stresses and phytohormones. No apparent phenotype was observed when individual HSP70 genes were overexpressed or knocked-out/down, but the double mutant hsp70-1 hsp70-4 and triple mutant hsp70-2 hsp70-4 hsp70-5 plants exhibited developmental phenotypes with shortened specific growth periods, curly and round leaves, twisted petioles, thin stems, and short siliques. Moreover, both mutants were hypersensitive to heat, cold, high glucose, salt and osmotic stress, but hyposensitive to abscisic acid. Genes related to flowering, and the cytokinin, brassinosteroid, and abscisic acid signaling pathways were differentially expressed in both mutants. Our studies suggest that, the individual HSP70 possibly performs both redundant and specific functions with the other members in the cytosolic/nuclear HSP70 subclass, and apart from enabling plants to cope with abiotic stresses, this subclass of cytosolic/nuclear HSP70 proteins also participates in diverse developmental processes and signaling pathways.

Proteomic analysis of Cucumis sativus cotyledons after glucohexaose treatment as a part of ROS accumulation related resistance mechanism.

BACKGROUND: Glucohexaose is a safe farm chemical used for pathogen prevention, which can induce systemic acquired resistance in cucumber. RESULTS: We found that glucohexaose treatment of cucumber plant induced an accumulation of the reactive oxidative species (ROS). Histochemistry showed sharp increases in O(2-) and H2O2 5 h after glucohexaose treatment. After 5 h, the O(2-) content decreased to a normal level, but the H2O2 content remained at a high level 10 h after glucohexaose treatment. And antioxidant enzymes were also changed after glucohexaose treatment. We also investigated the relationship between ROS accumulation and glucohexaose-induced proteome alteration using 2D electrophoresis coupled with MS/MS. 54 protein spots, which enhanced expression under glucohexaose treatment but suppressed the expression by application of DPI and DMTU, have been identified. CONCLUSION: Our study showed the accumulation of ROS is a part of mechanism of glucohexaose induced resistance in cucumber cotyledons. The up-regulated proteins identified by MS such as PP2C and antioxidation proteins are important in ROS signaling. It will be interesting to find out the regulatory mechanism underlying the induction of these proteins via ROS, and provide some clues to the mechanism of glucohexaose-induced resistance.

Regeneration of axons in injured spinal cord by activation of bone morphogenetic protein/Smad1 signaling pathway in adult neurons.

Axon growth potential is highest in young neurons but diminishes with age, thus becoming a significant obstacle to axonal regeneration after injury in maturity. The mechanism for the decline is incompletely understood, and no effective clinical treatment is available to rekindle innate growth capability. Here, we show that Smad1-dependent bone morphogenetic protein (BMP) signaling is developmentally regulated and governs axonal growth in dorsal root ganglion (DRG) neurons. Down-regulation of the pathway contributes to the age-related decline of the axon growth potential. Reactivating Smad1 selectively in adult DRG neurons results in sensory axon regeneration in a mouse model of spinal cord injury (SCI). Smad1 signaling can be effectively manipulated by an adeno-associated virus (AAV) vector encoding BMP4 delivered by a clinically applicable and minimally invasive technique, an approach devoid of unwanted abnormalities in mechanosensation or pain perception. Importantly, transected axons are able to regenerate even when the AAV treatment is delivered after SCI, thus mimicking a clinically relevant scenario. Together, our results identify a therapeutic target to promote axonal regeneration after SCI.

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.