Molecular Pathways of Colon Inflammation Induced by Cancer Immunotherapy.

Checkpoint blockade with antibodies specific for the PD-1 and CTLA-4 inhibitory receptors can induce durable responses in a wide range of human cancers. However, the immunological mechanisms responsible for severe inflammatory side effects remain poorly understood. Here we report a comprehensive single-cell analysis of immune cell populations in colitis, a common and severe side effect of checkpoint blockade. We observed a striking accumulation of CD8 T cells with highly cytotoxic and proliferative states and no evidence of regulatory T cell depletion. T cell receptor (TCR) sequence analysis demonstrated that a substantial fraction of colitis-associated CD8 T cells originated from tissue-resident populations, explaining the frequently early onset of colitis symptoms following treatment initiation. Our analysis also identified cytokines, chemokines, and surface receptors that could serve as therapeutic targets for colitis and potentially other inflammatory side effects of checkpoint blockade.

Direct Promoter Repression by BCL11A Controls the Fetal to Adult Hemoglobin Switch.

Fetal hemoglobin (HbF, α2γ2) level is genetically controlled and modifies severity of adult hemoglobin (HbA, α2ß2) disorders, sickle cell disease, and ß-thalassemia. Common genetic variation affects expression of BCL11A, a regulator of HbF silencing. To uncover how BCL11A supports the developmental switch from γ- to ß- globin, we use a functional assay and protein binding microarray to establish a requirement for a zinc-finger cluster in BCL11A in repression and identify a preferred DNA recognition sequence. This motif appears in embryonic and fetal-expressed globin promoters and is duplicated in γ-globin promoters. The more distal of the duplicated motifs is mutated in individuals with hereditary persistence of HbF. Using the CUT&RUN approach to map protein binding sites in erythroid cells, we demonstrate BCL11A occupancy preferentially at the distal motif, which can be disrupted by editing the promoter. Our findings reveal that direct γ-globin gene promoter repression by BCL11A underlies hemoglobin switching.

Mapping the Mouse Cell Atlas by Microwell-Seq.

Single-cell RNA sequencing (scRNA-seq) technologies are poised to reshape the current cell-type classification system. However, a transcriptome-based single-cell atlas has not been achieved for complex mammalian systems. Here, we developed Microwell-seq, a high-throughput and low-cost scRNA-seq platform using simple, inexpensive devices. Using Microwell-seq, we analyzed more than 400,000 single cells covering all of the major mouse organs and constructed a basic scheme for a mouse cell atlas (MCA). We reveal a single-cell hierarchy for many tissues that have not been well characterized previously. We built a web-based "single-cell MCA analysis" pipeline that accurately defines cell types based on single-cell digital expression. Our study demonstrates the wide applicability of the Microwell-seq technology and MCA resource.

Altered nucleocytoplasmic export of adenosine-rich circRNAs by PABPC1 contributes to neuronal function.

Circular RNAs (circRNAs) are upregulated during neurogenesis. Where and how circRNAs are localized and what roles they play during this process have remained elusive. Comparing the nuclear and cytoplasmic circRNAs between H9 cells and H9-derived forebrain (FB) neurons, we identify that a subset of adenosine (A)-rich circRNAs are restricted in H9 nuclei but exported to cytosols upon differentiation. Such a subcellular relocation of circRNAs is modulated by the poly(A)-binding protein PABPC1. In the H9 nucleus, newly produced (A)-rich circRNAs are bound by PABPC1 and trapped by the nuclear basket protein TPR to prevent their export. Modulating (A)-rich motifs in circRNAs alters their subcellular localization, and introducing (A)-rich circRNAs in H9 cytosols results in mRNA translation suppression. Moreover, decreased nuclear PABPC1 upon neuronal differentiation enables the export of (A)-rich circRNAs, including circRTN4(2,3), which is required for neurite outgrowth. These findings uncover subcellular localization features of circRNAs, linking their processing and function during neurogenesis.

Acquired Tissue-Specific Promoter Bivalency Is a Basis for PRC2 Necessity in Adult Cells.

Bivalent promoters in embryonic stem cells (ESCs) carry methylation marks on two lysine residues, K4 and K27, in histone3 (H3). K4me2/3 is generally considered to promote transcription, and Polycomb Repressive Complex 2 (PRC2) places K27me3, which is erased at lineage-restricted genes when ESCs differentiate in culture. Molecular defects in various PRC2 null adult tissues lack a unifying explanation. We found that epigenomes in adult mouse intestine and other self-renewing tissues show fewer and distinct bivalent promoters compared to ESCs. Groups of tissue-specific genes that carry bivalent marks are repressed, despite the presence of promoter H3K4me2/3. These are the predominant genes de-repressed in PRC2-deficient adult cells, where aberrant expression is proportional to the H3K4me2/3 levels observed at their promoters in wild-type cells. Thus, in adult animals, PRC2 specifically represses genes with acquired, tissue-restricted promoter bivalency. These findings provide new insights into specificity in chromatin-based gene regulation.

MLL::AF9 degradation induces rapid changes in transcriptional elongation and subsequent loss of an active chromatin landscape.

MLL rearrangements produce fusion oncoproteins that drive leukemia development, but the direct effects of MLL-fusion inactivation remain poorly defined. We designed models with degradable MLL::AF9 where treatment with small molecules induces rapid degradation. We leveraged the kinetics of this system to identify a core subset of MLL::AF9 target genes where MLL::AF9 degradation induces changes in transcriptional elongation within 15 minutes. MLL::AF9 degradation subsequently causes loss of a transcriptionally active chromatin landscape. We used this insight to assess the effectiveness of small molecules that target members of the MLL::AF9 multiprotein complex, specifically DOT1L and MENIN. Combined DOT1L/MENIN inhibition resembles MLL::AF9 degradation, whereas single-agent treatment has more modest effects on MLL::AF9 occupancy and gene expression. Our data show that MLL::AF9 degradation leads to decreases in transcriptional elongation prior to changes in chromatin landscape at select loci and that combined inhibition of chromatin complexes releases the MLL::AF9 oncoprotein from chromatin globally.

The histone demethylase UTX regulates the lineage-specific epigenetic program of invariant natural killer T cells.

Invariant natural killer T cells (iNKT cells) are innate-like lymphocytes that protect against infection, autoimmune disease and cancer. However, little is known about the epigenetic regulation of iNKT cell development. Here we found that the H3K27me3 histone demethylase UTX was an essential cell-intrinsic factor that controlled an iNKT-cell lineage-specific gene-expression program and epigenetic landscape in a demethylase-activity-dependent manner. UTX-deficient iNKT cells exhibited impaired expression of iNKT cell signature genes due to a decrease in activation-associated H3K4me3 marks and an increase in repressive H3K27me3 marks within the promoters occupied by UTX. We found that JunB regulated iNKT cell development and that the expression of genes that were targets of both JunB and the iNKT cell master transcription factor PLZF was UTX dependent. We identified iNKT cell super-enhancers and demonstrated that UTX-mediated regulation of super-enhancer accessibility was a key mechanism for commitment to the iNKT cell lineage. Our findings reveal how UTX regulates the development of iNKT cells through multiple epigenetic mechanisms.

Reprogramming committed murine blood cells to induced hematopoietic stem cells with defined factors.

Hematopoietic stem cells (HSCs) sustain blood formation throughout life and are the functional units of bone marrow transplantation. We show that transient expression of six transcription factors Run1t1, Hlf, Lmo2, Prdm5, Pbx1, and Zfp37 imparts multilineage transplantation potential onto otherwise committed lymphoid and myeloid progenitors and myeloid effector cells. Inclusion of Mycn and Meis1 and use of polycistronic viruses increase reprogramming efficacy. The reprogrammed cells, designated induced-HSCs (iHSCs), possess clonal multilineage differentiation potential, reconstitute stem/progenitor compartments, and are serially transplantable. Single-cell analysis revealed that iHSCs derived under optimal conditions exhibit a gene expression profile that is highly similar to endogenous HSCs. These findings demonstrate that expression of a set of defined factors is sufficient to activate the gene networks governing HSC functional identity in committed blood cells. Our results raise the prospect that blood cell reprogramming may be a strategy for derivation of transplantable stem cells for clinical application.

CDK7 inhibition suppresses super-enhancer-linked oncogenic transcription in MYCN-driven cancer.

The MYC oncoproteins are thought to stimulate tumor cell growth and proliferation through amplification of gene transcription, a mechanism that has thwarted most efforts to inhibit MYC function as potential cancer therapy. Using a covalent inhibitor of cyclin-dependent kinase 7 (CDK7) to disrupt the transcription of amplified MYCN in neuroblastoma cells, we demonstrate downregulation of the oncoprotein with consequent massive suppression of MYCN-driven global transcriptional amplification. This response translated to significant tumor regression in a mouse model of high-risk neuroblastoma, without the introduction of systemic toxicity. The striking treatment selectivity of MYCN-overexpressing cells correlated with preferential downregulation of super-enhancer-associated genes, including MYCN and other known oncogenic drivers in neuroblastoma. These results indicate that CDK7 inhibition, by selectively targeting the mechanisms that promote global transcriptional amplification in tumor cells, may be useful therapy for cancers that are driven by MYC family oncoproteins.

A vaccine targeting resistant tumours by dual T cell plus NK cell attack.

Most cancer vaccines target peptide antigens, necessitating personalization owing to the vast inter-individual diversity in major histocompatibility complex (MHC) molecules that present peptides to T cells. Furthermore, tumours frequently escape T cell-mediated immunity through mechanisms that interfere with peptide presentation1. Here we report a cancer vaccine that induces a coordinated attack by diverse T cell and natural killer (NK) cell populations. The vaccine targets the MICA and MICB (MICA/B) stress proteins expressed by many human cancers as a result of DNA damage2. MICA/B serve as ligands for the activating NKG2D receptor on T cells and NK cells, but tumours evade immune recognition by proteolytic MICA/B cleavage3,4. Vaccine-induced antibodies increase the density of MICA/B proteins on the surface of tumour cells by inhibiting proteolytic shedding, enhance presentation of tumour antigens by dendritic cells to T cells and augment the cytotoxic function of NK cells. Notably, this vaccine maintains efficacy against MHC class I-deficient tumours resistant to cytotoxic T cells through the coordinated action of NK cells and CD4+ T cells. The vaccine is also efficacious in a clinically important setting: immunization following surgical removal of primary, highly metastatic tumours inhibits the later outgrowth of metastases. This vaccine design enables protective immunity even against tumours with common escape mutations.

Integrated spatial genomics reveals global architecture of single nuclei.

Identifying the relationships between chromosome structures, nuclear bodies, chromatin states and gene expression is an overarching goal of nuclear-organization studies1-4. Because individual cells appear to be highly variable at all these levels5, it is essential to map different modalities in the same cells. Here we report the imaging of 3,660 chromosomal loci in single mouse embryonic stem (ES) cells using DNA seqFISH+, along with 17 chromatin marks and subnuclear structures by sequential immunofluorescence and the expression profile of 70 RNAs. Many loci were invariably associated with immunofluorescence marks in single mouse ES cells. These loci form 'fixed points' in the nuclear organizations of single cells and often appear on the surfaces of nuclear bodies and zones defined by combinatorial chromatin marks. Furthermore, highly expressed genes appear to be pre-positioned to active nuclear zones, independent of bursting dynamics in single cells. Our analysis also uncovered several distinct mouse ES cell subpopulations with characteristic combinatorial chromatin states. Using clonal analysis, we show that the global levels of some chromatin marks, such as H3 trimethylation at lysine 27 (H3K27me3) and macroH2A1 (mH2A1), are heritable over at least 3-4 generations, whereas other marks fluctuate on a faster time scale. This seqFISH+-based spatial multimodal approach can be used to explore nuclear organization and cell states in diverse biological systems.

Large-scale benchmarking of circRNA detection tools reveals large differences in sensitivity but not in precision.

The detection of circular RNA molecules (circRNAs) is typically based on short-read RNA sequencing data processed using computational tools. Numerous such tools have been developed, but a systematic comparison with orthogonal validation is missing. Here, we set up a circRNA detection tool benchmarking study, in which 16 tools detected more than 315,000 unique circRNAs in three deeply sequenced human cell types. Next, 1,516 predicted circRNAs were validated using three orthogonal methods. Generally, tool-specific precision is high and similar (median of 98.8%, 96.3% and 95.5% for qPCR, RNase R and amplicon sequencing, respectively) whereas the sensitivity and number of predicted circRNAs (ranging from 1,372 to 58,032) are the most significant differentiators. Of note, precision values are lower when evaluating low-abundance circRNAs. We also show that the tools can be used complementarily to increase detection sensitivity. Finally, we offer recommendations for future circRNA detection and validation.

Targeting of the CD161 inhibitory receptor enhances T-cell-mediated immunity against hematological malignancies.

ABSTRACT: The CD161 inhibitory receptor is highly upregulated by tumor-infiltrating T cells in multiple human solid tumor types, and its ligand, CLEC2D, is expressed by both tumor cells and infiltrating myeloid cells. Here, we assessed the role of the CD161 receptor in hematological malignancies. Systematic analysis of CLEC2D expression using the Cancer Cell Line Encyclopedia revealed that CLEC2D messenger RNA was most abundant in hematological malignancies, including B-cell and T-cell lymphomas as well as lymphocytic and myelogenous leukemias. CLEC2D protein was detected by flow cytometry on a panel of cell lines representing a diverse set of hematological malignancies. We, therefore, used yeast display to generate a panel of high-affinity, fully human CD161 monoclonal antibodies (mAbs) that blocked CLEC2D binding. These mAbs were specific for CD161 and had a similar affinity for human and nonhuman primate CD161, a property relevant for clinical translation. A high-affinity CD161 mAb enhanced key aspects of T-cell function, including cytotoxicity, cytokine production, and proliferation, against B-cell lines originating from patients with acute lymphoblastic leukemia, diffuse large B-cell lymphoma, and Burkitt lymphoma. In humanized mouse models, this CD161 mAb enhanced T-cell-mediated immunity, resulting in a significant survival benefit. Single cell RNA-seq data demonstrated that CD161 mAb treatment enhanced expression of cytotoxicity genes by CD4 T cells as well as a tissue-residency program by CD4 and CD8 T cells that is associated with favorable survival outcomes in multiple human cancer types. These fully human mAbs, thus, represent potential immunotherapy agents for hematological malignancies.

RNAlight: a machine learning model to identify nucleotide features determining RNA subcellular localization.

Different RNAs have distinct subcellular localizations. However, nucleotide features that determine these distinct distributions of lncRNAs and mRNAs have yet to be fully addressed. Here, we develop RNAlight, a machine learning model based on LightGBM, to identify nucleotide k-mers contributing to the subcellular localizations of mRNAs and lncRNAs. With the Tree SHAP algorithm, RNAlight extracts nucleotide features for cytoplasmic or nuclear localization of RNAs, indicating the sequence basis for distinct RNA subcellular localizations. By assembling k-mers to sequence features and subsequently mapping to known RBP-associated motifs, different types of sequence features and their associated RBPs were additionally uncovered for lncRNAs and mRNAs with distinct subcellular localizations. Finally, we extended RNAlight to precisely predict the subcellular localizations of other types of RNAs, including snRNAs, snoRNAs and different circular RNA transcripts, suggesting the generality of using RNAlight for RNA subcellular localization prediction.

Transcriptome-scale super-resolved imaging in tissues by RNA seqFISH.

Imaging the transcriptome in situ with high accuracy has been a major challenge in single-cell biology, which is particularly hindered by the limits of optical resolution and the density of transcripts in single cells1-5. Here we demonstrate an evolution of sequential fluorescence in situ hybridization (seqFISH+). We show that seqFISH+ can image mRNAs for 10,000 genes in single cells-with high accuracy and sub-diffraction-limit resolution-in the cortex, subventricular zone and olfactory bulb of mouse brain, using a standard confocal microscope. The transcriptome-level profiling of seqFISH+ allows unbiased identification of cell classes and their spatial organization in tissues. In addition, seqFISH+ reveals subcellular mRNA localization patterns in cells and ligand-receptor pairs across neighbouring cells. This technology demonstrates the ability to generate spatial cell atlases and to perform discovery-driven studies of biological processes in situ.

BORIS promotes chromatin regulatory interactions in treatment-resistant cancer cells.

The CCCTC-binding factor (CTCF), which anchors DNA loops that organize the genome into structural domains, has a central role in gene control by facilitating or constraining interactions between genes and their regulatory elements1,2. In cancer cells, the disruption of CTCF binding at specific loci by somatic mutation3,4 or DNA hypermethylation5 results in the loss of loop anchors and consequent activation of oncogenes. By contrast, the germ-cell-specific paralogue of CTCF, BORIS (brother of the regulator of imprinted sites, also known as CTCFL)6, is overexpressed in several cancers7-9, but its contributions to the malignant phenotype remain unclear. Here we show that aberrant upregulation of BORIS promotes chromatin interactions in ALK-mutated, MYCN-amplified neuroblastoma10 cells that develop resistance to ALK inhibition. These cells are reprogrammed to a distinct phenotypic state during the acquisition of resistance, a process defined by the initial loss of MYCN expression followed by subsequent overexpression of BORIS and a concomitant switch in cellular dependence from MYCN to BORIS. The resultant BORIS-regulated alterations in chromatin looping lead to the formation of super-enhancers that drive the ectopic expression of a subset of proneural transcription factors that ultimately define the resistance phenotype. These results identify a previously unrecognized role of BORIS-to promote regulatory chromatin interactions that support specific cancer phenotypes.

Construction of a cross-species cell landscape at single-cell level.

Individual cells are basic units of life. Despite extensive efforts to characterize the cellular heterogeneity of different organisms, cross-species comparisons of landscape dynamics have not been achieved. Here, we applied single-cell RNA sequencing (scRNA-seq) to map organism-level cell landscapes at multiple life stages for mice, zebrafish and Drosophila. By integrating the comprehensive dataset of > 2.6 million single cells, we constructed a cross-species cell landscape and identified signatures and common pathways that changed throughout the life span. We identified structural inflammation and mitochondrial dysfunction as the most common hallmarks of organism aging, and found that pharmacological activation of mitochondrial metabolism alleviated aging phenotypes in mice. The cross-species cell landscape with other published datasets were stored in an integrated online portal-Cell Landscape. Our work provides a valuable resource for studying lineage development, maturation and aging.

How many cell types are there in nature? How do they change during the life cycle? These are two fundamental questions that researchers have been trying to understand in the area of biology. In this study, single-cell mRNA sequencing data were used to profile over 2.6 million individual cells from mice, zebrafish and Drosophila at different life stages, 1.3 million of which were newly collected. The comprehensive datasets allow investigators to construct a cross-species cell landscape that helps to reveal the conservation and diversity of cell taxonomies at genetic and regulatory levels. The resources in this study are assembled into a publicly available website at http://bis.zju.edu.cn/cellatlas/.

miR-309a is a regulator of ovarian development in the oriental fruit fly Bactrocera dorsalis.

Fecundity is arguably one of the most important life history traits, as it is closely tied to fitness. Most arthropods are recognized for their extreme reproductive capacity. For example, a single female of the oriental fruit fly Bactrocera dorsalis, a highly invasive species that is one of the most destructive agricultural pests worldwide, can lay more than 3000 eggs during its life span. The ovary is crucial for insect reproduction and its development requires further investigation at the molecular level. We report here that miR-309a is a regulator of ovarian development in B. dorsalis. Our bioinformatics and molecular studies have revealed that miR-309a binds the transcription factor pannier (GATA-binding factor A/pnr), and this activates yolk vitellogenin 2 (Vg 2) and vitellogenin receptor (VgR) advancing ovarian development. We further show that miR-309a is under the control of juvenile hormone (JH) and independent from 20-hydroxyecdysone. Thus, we identified a JH-controlled miR-309a/pnr axis that regulates Vg2 and VgR to control the ovarian development. This study has further enhanced our understanding of molecular mechanisms governing ovarian development and insect reproduction. It provides a background for identifying targets for controlling important Dipteran pests.

Orthogonal multimodality integration and clustering in single-cell data.

Multimodal integration combines information from different sources or modalities to gain a more comprehensive understanding of a phenomenon. The challenges in multi-omics data analysis lie in the complexity, high dimensionality, and heterogeneity of the data, which demands sophisticated computational tools and visualization methods for proper interpretation and visualization of multi-omics data. In this paper, we propose a novel method, termed Orthogonal Multimodality Integration and Clustering (OMIC), for analyzing CITE-seq. Our approach enables researchers to integrate multiple sources of information while accounting for the dependence among them. We demonstrate the effectiveness of our approach using CITE-seq data sets for cell clustering. Our results show that our approach outperforms existing methods in terms of accuracy, computational efficiency, and interpretability. We conclude that our proposed OMIC method provides a powerful tool for multimodal data analysis that greatly improves the feasibility and reliability of integrated data.