An esophagus cell atlas reveals dynamic rewiring during active eosinophilic esophagitis and remission.

Coordinated cell interactions within the esophagus maintain homeostasis, and disruption can lead to eosinophilic esophagitis (EoE), a chronic inflammatory disease with poorly understood pathogenesis. We profile 421,312 individual cells from the esophageal mucosa of 7 healthy and 15 EoE participants, revealing 60 cell subsets and functional alterations in cell states, compositions, and interactions that highlight previously unclear features of EoE. Active disease displays enrichment of ALOX15+ macrophages, PRDM16+ dendritic cells expressing the EoE risk gene ATP10A, and cycling mast cells, with concomitant reduction of TH17 cells. Ligand-receptor expression uncovers eosinophil recruitment programs, increased fibroblast interactions in disease, and IL-9+IL-4+IL-13+ TH2 and endothelial cells as potential mast cell interactors. Resolution of inflammation-associated signatures includes mast and CD4+ TRM cell contraction and cell type-specific downregulation of eosinophil chemoattractant, growth, and survival factors. These cellular alterations in EoE and remission advance our understanding of eosinophilic inflammation and opportunities for therapeutic intervention.

The improvement of agronomic performances in the cold weather conditions for perennial wheatgrass by crossing Thinopyrum intermedium with wheat-Th. intermedium partial amphiploids.

Thinopyrum intermedium (2n=6x=42, StStJrJrJvsJvs) is resistant or tolerant to biotic and abiotic stresses, making it suitable for developing perennial crops and forage. Through five cycles of selection, we developed 24 perennial wheatgrass lines, designated 19HSC-Q and 20HSC-Z, by crossing wheat-Th. intermedium partial amphiploids with Th. intermedium. The cold resistance, morphological performance, chromosome composition, and yield components of these perennial lines were investigated from 2019 to 2022. Six lines of 19HSC-Q had higher 1,000-kernel weight, grains per spike, and tiller number than Th. intermedium, as well as surviving -30°C in winter. Lines 19HSC-Q14, 19HSC-Q18, and 19HSC-Q20 had the best performances for grain number per spike and 1,000-kernel weight. The 20HSC-Z lines, 20HSC-Z1, 20HSC-Z2, and 20HSC-Z3, were able to survive in the cold winter in Harbin and had been grown for two years. Sequential multicolor GISH analysis revealed that the Jvs subgenome of Th. intermedium were divided into two karyotypes, three pairs of type-I Jvs chromosomes and four pairs of type-II Jvs chromosomes. Both Th. intermedium and the 24 advanced perennial wheatgrass lines had similar chromosome compositions, but the translocations among subgenome chromosomes were detected in some lines with prominent agronomic traits, such as 19HSC-Q11, 19HSC-Q14, 19HSC-Q18, 19HSC-Q20, and the three 20HSC-Z lines. The chromosome aberrations were distinguished into two types: the large fragment translocation with St-Jr, Jvs-St, Jr-IIJvs, and Jvs-Jr and the small fragment introgression of Jr-St, St-IJvs, and Jvs-Jr. These chromosomal variations can be used to further analyze the relationship between the subgenomes and phenotypes of Th. intermedium. The results of this study provide valuable materials for the next selection cycle of cold-resistant perennial wheatgrass.

Temporal single-cell atlas of non-neuronal retinal cells reveals dynamic, coordinated multicellular responses to central nervous system injury.

Non-neuronal cells are key to the complex cellular interplay that follows central nervous system insult. To understand this interplay, we generated a single-cell atlas of immune, glial and retinal pigment epithelial cells from adult mouse retina before and at multiple time points after axonal transection. We identified rare subsets in naive retina, including interferon (IFN)-response glia and border-associated macrophages, and delineated injury-induced changes in cell composition, expression programs and interactions. Computational analysis charted a three-phase multicellular inflammatory cascade after injury. In the early phase, retinal macroglia and microglia were reactivated, providing chemotactic signals concurrent with infiltration of CCR2+ monocytes from the circulation. These cells differentiated into macrophages in the intermediate phase, while an IFN-response program, likely driven by microglia-derived type I IFN, was activated across resident glia. The late phase indicated inflammatory resolution. Our findings provide a framework to decipher cellular circuitry, spatial relationships and molecular interactions following tissue injury.

Disruption of the IL-33-ST2-AKT signaling axis impairs neurodevelopment by inhibiting microglial metabolic adaptation and phagocytic function.

To accommodate the changing needs of the developing brain, microglia must undergo substantial morphological, phenotypic, and functional reprogramming. Here, we examined whether cellular metabolism regulates microglial function during neurodevelopment. Microglial mitochondria bioenergetics correlated with and were functionally coupled to phagocytic activity in the developing brain. Transcriptional profiling of microglia with diverse metabolic profiles revealed an activation signature wherein the interleukin (IL)-33 signaling axis is associated with phagocytic activity. Genetic perturbation of IL-33 or its receptor ST2 led to microglial dystrophy, impaired synaptic function, and behavioral abnormalities. Conditional deletion of Il33 from astrocytes or Il1rl1, encoding ST2, in microglia increased susceptibility to seizures. Mechanistically, IL-33 promoted mitochondrial activity and phagocytosis in an AKT-dependent manner. Mitochondrial metabolism and AKT activity were temporally regulated in vivo. Thus, a microglia-astrocyte circuit mediated by the IL-33-ST2-AKT signaling axis supports microglial metabolic adaptation and phagocytic function during early development, with implications for neurodevelopmental and neuropsychiatric disorders.

Deep generative model embedding of single-cell RNA-Seq profiles on hyperspheres and hyperbolic spaces.

Single-cell RNA-Seq (scRNA-seq) is invaluable for studying biological systems. Dimensionality reduction is a crucial step in interpreting the relation between cells in scRNA-seq data. However, current dimensionality reduction methods are often confounded by multiple simultaneous technical and biological variability, result in "crowding" of cells in the center of the latent space, or inadequately capture temporal relationships. Here, we introduce scPhere, a scalable deep generative model to embed cells into low-dimensional hyperspherical or hyperbolic spaces to accurately represent scRNA-seq data. ScPhere addresses multi-level, complex batch factors, facilitates the interactive visualization of large datasets, resolves cell crowding, and uncovers temporal trajectories. We demonstrate scPhere on nine large datasets in complex tissue from human patients or animal development. Our results show how scPhere facilitates the interpretation of scRNA-seq data by generating batch-invariant embeddings to map data from new individuals, identifies cell types affected by biological variables, infers cells' spatial positions in pre-defined biological specimens, and highlights complex cellular relations.

Chromatin Potential Identified by Shared Single-Cell Profiling of RNA and Chromatin.

Cell differentiation and function are regulated across multiple layers of gene regulation, including modulation of gene expression by changes in chromatin accessibility. However, differentiation is an asynchronous process precluding a temporal understanding of regulatory events leading to cell fate commitment. Here we developed simultaneous high-throughput ATAC and RNA expression with sequencing (SHARE-seq), a highly scalable approach for measurement of chromatin accessibility and gene expression in the same single cell, applicable to different tissues. Using 34,774 joint profiles from mouse skin, we develop a computational strategy to identify cis-regulatory interactions and define domains of regulatory chromatin (DORCs) that significantly overlap with super-enhancers. During lineage commitment, chromatin accessibility at DORCs precedes gene expression, suggesting that changes in chromatin accessibility may prime cells for lineage commitment. We computationally infer chromatin potential as a quantitative measure of chromatin lineage-priming and use it to predict cell fate outcomes. SHARE-seq is an extensible platform to study regulatory circuitry across diverse cells in tissues.

Systematic comparison of single-cell and single-nucleus RNA-sequencing methods.

The scale and capabilities of single-cell RNA-sequencing methods have expanded rapidly in recent years, enabling major discoveries and large-scale cell mapping efforts. However, these methods have not been systematically and comprehensively benchmarked. Here, we directly compare seven methods for single-cell and/or single-nucleus profiling-selecting representative methods based on their usage and our expertise and resources to prepare libraries-including two low-throughput and five high-throughput methods. We tested the methods on three types of samples: cell lines, peripheral blood mononuclear cells and brain tissue, generating 36 libraries in six separate experiments in a single center. To directly compare the methods and avoid processing differences introduced by the existing pipelines, we developed scumi, a flexible computational pipeline that can be used with any single-cell RNA-sequencing method. We evaluated the methods for both basic performance, such as the structure and alignment of reads, sensitivity and extent of multiplets, and for their ability to recover known biological information in the samples.

Author Correction: Systematic comparison of single-cell and single-nucleus RNA-sequencing methods.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Single-Cell Transcriptome Analysis Reveals Disease-Defining T-cell Subsets in the Tumor Microenvironment of Classic Hodgkin Lymphoma.

Hodgkin lymphoma is characterized by an extensively dominant tumor microenvironment (TME) composed of different types of noncancerous immune cells with rare malignant cells. Characterization of the cellular components and their spatial relationship is crucial to understanding cross-talk and therapeutic targeting in the TME. We performed single-cell RNA sequencing of more than 127,000 cells from 22 Hodgkin lymphoma tissue specimens and 5 reactive lymph nodes, profiling for the first time the phenotype of the Hodgkin lymphoma-specific immune microenvironment at single-cell resolution. Single-cell expression profiling identified a novel Hodgkin lymphoma-associated subset of T cells with prominent expression of the inhibitory receptor LAG3, and functional analyses established this LAG3+ T-cell population as a mediator of immunosuppression. Multiplexed spatial assessment of immune cells in the microenvironment also revealed increased LAG3+ T cells in the direct vicinity of MHC class II-deficient tumor cells. Our findings provide novel insights into TME biology and suggest new approaches to immune-checkpoint targeting in Hodgkin lymphoma. SIGNIFICANCE: We provide detailed functional and spatial characteristics of immune cells in classic Hodgkin lymphoma at single-cell resolution. Specifically, we identified a regulatory T-cell-like immunosuppressive subset of LAG3+ T cells contributing to the immune-escape phenotype. Our insights aid in the development of novel biomarkers and combination treatment strategies targeting immune checkpoints.See related commentary by Fisher and Oh, p. 342.This article is highlighted in the In This Issue feature, p. 327.

Transcriptional Atlas of Intestinal Immune Cells Reveals that Neuropeptide α-CGRP Modulates Group 2 Innate Lymphoid Cell Responses.

Signaling abnormalities in immune responses in the small intestine can trigger chronic type 2 inflammation involving interaction of multiple immune cell types. To systematically characterize this response, we analyzed 58,067 immune cells from the mouse small intestine by single-cell RNA sequencing (scRNA-seq) at steady state and after induction of a type 2 inflammatory reaction to ovalbumin (OVA). Computational analysis revealed broad shifts in both cell-type composition and cell programs in response to the inflammation, especially in group 2 innate lymphoid cells (ILC2s). Inflammation induced the expression of exon 5 of Calca, which encodes the alpha-calcitonin gene-related peptide (α-CGRP), in intestinal KLRG1+ ILC2s. α-CGRP antagonized KLRG1+ ILC2s proliferation but promoted IL-5 expression. Genetic perturbation of α-CGRP increased the proportion of intestinal KLRG1+ ILC2s. Our work highlights a model where α-CGRP-mediated neuronal signaling is critical for suppressing ILC2 expansion and maintaining homeostasis of the type 2 immune machinery.

Interpretable dimensionality reduction of single cell transcriptome data with deep generative models.

Single-cell RNA-sequencing has great potential to discover cell types, identify cell states, trace development lineages, and reconstruct the spatial organization of cells. However, dimension reduction to interpret structure in single-cell sequencing data remains a challenge. Existing algorithms are either not able to uncover the clustering structures in the data or lose global information such as groups of clusters that are close to each other. We present a robust statistical model, scvis, to capture and visualize the low-dimensional structures in single-cell gene expression data. Simulation results demonstrate that low-dimensional representations learned by scvis preserve both the local and global neighbor structures in the data. In addition, scvis is robust to the number of data points and learns a probabilistic parametric mapping function to add new data points to an existing embedding. We then use scvis to analyze four single-cell RNA-sequencing datasets, exemplifying interpretable two-dimensional representations of the high-dimensional single-cell RNA-sequencing data.

densityCut: an efficient and versatile topological approach for automatic clustering of biological data.

MOTIVATION: Many biological data processing problems can be formalized as clustering problems to partition data points into sensible and biologically interpretable groups. RESULTS: This article introduces densityCut, a novel density-based clustering algorithm, which is both time- and space-efficient and proceeds as follows: densityCut first roughly estimates the densities of data points from a K-nearest neighbour graph and then refines the densities via a random walk. A cluster consists of points falling into the basin of attraction of an estimated mode of the underlining density function. A post-processing step merges clusters and generates a hierarchical cluster tree. The number of clusters is selected from the most stable clustering in the hierarchical cluster tree. Experimental results on ten synthetic benchmark datasets and two microarray gene expression datasets demonstrate that densityCut performs better than state-of-the-art algorithms for clustering biological datasets. For applications, we focus on the recent cancer mutation clustering and single cell data analyses, namely to cluster variant allele frequencies of somatic mutations to reveal clonal architectures of individual tumours, to cluster single-cell gene expression data to uncover cell population compositions, and to cluster single-cell mass cytometry data to detect communities of cells of the same functional states or types. densityCut performs better than competing algorithms and is scalable to large datasets. AVAILABILITY AND IMPLEMENTATION: Data and the densityCut R package is available from https://bitbucket.org/jerry00/densitycut_dev CONTACT: : condon@cs.ubc.ca or sshah@bccrc.ca or jiaruid@cs.ubc.ca SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Identification of the epigenetic reader CBX2 as a potential drug target in advanced prostate cancer.

BACKGROUND: While localized prostate cancer (PCa) can be effectively cured, metastatic disease inevitably progresses to a lethal state called castration-resistant prostate cancer (CRPC). Emerging evidence suggests that aberrant epigenetic repression by the polycomb group (PcG) complexes fuels PCa progression, providing novel therapeutic opportunities. RESULTS: In the search for potential epigenetic drivers of CRPC, we analyzed the molecular profile of PcG members in patient-derived xenografts and clinical samples. Overall, our results identify the PcG protein and methyl-lysine reader CBX2 as a potential therapeutic target in advanced PCa. We report that CBX2 was recurrently up-regulated in metastatic CRPC and that elevated CBX2 expression was correlated with poor clinical outcome in PCa cohorts. Furthermore, CBX2 depletion abrogated cell viability and induced caspase 3-mediated apoptosis in metastatic PCa cell lines. Mechanistically explaining this phenotype, microarray analysis in CBX2-depleted cells revealed that CBX2 controls the expression of many key regulators of cell proliferation and metastasis. CONCLUSIONS: Taken together, this study provides the first evidence that CBX2 inhibition induces cancer cell death, positioning CBX2 as an attractive drug target in lethal CRPC.

Systematic analysis of somatic mutations impacting gene expression in 12 tumour types.

We present a novel hierarchical Bayes statistical model, xseq, to systematically quantify the impact of somatic mutations on expression profiles. We establish the theoretical framework and robust inference characteristics of the method using computational benchmarking. We then use xseq to analyse thousands of tumour data sets available through The Cancer Genome Atlas, to systematically quantify somatic mutations impacting expression profiles. We identify 30 novel cis-effect tumour suppressor gene candidates, enriched in loss-of-function mutations and biallelic inactivation. Analysis of trans-effects of mutations and copy number alterations with xseq identifies mutations in 150 genes impacting expression networks, with 89 novel predictions. We reveal two important novel characteristics of mutation impact on expression: (1) patients harbouring known driver mutations exhibit different downstream gene expression consequences; (2) expression patterns for some mutations are stable across tumour types. These results have critical implications for identification and interpretation of mutations with consequent impact on transcription in cancer.

Cis-regulatory somatic mutations and gene-expression alteration in B-cell lymphomas.

BACKGROUND: With the rapid increase of whole-genome sequencing of human cancers, an important opportunity to analyze and characterize somatic mutations lying within cis-regulatory regions has emerged. A focus on protein-coding regions to identify nonsense or missense mutations disruptive to protein structure and/or function has led to important insights; however, the impact on gene expression of mutations lying within cis-regulatory regions remains under-explored. We analyzed somatic mutations from 84 matched tumor-normal whole genomes from B-cell lymphomas with accompanying gene expression measurements to elucidate the extent to which these cancers are disrupted by cis-regulatory mutations. RESULTS: We characterize mutations overlapping a high quality set of well-annotated transcription factor binding sites (TFBSs), covering a similar portion of the genome as protein-coding exons. Our results indicate that cis-regulatory mutations overlapping predicted TFBSs are enriched in promoter regions of genes involved in apoptosis or growth/proliferation. By integrating gene expression data with mutation data, our computational approach culminates with identification of cis-regulatory mutations most likely to participate in dysregulation of the gene expression program. The impact can be measured along with protein-coding mutations to highlight key mutations disrupting gene expression and pathways in cancer. CONCLUSIONS: Our study yields specific genes with disrupted expression triggered by genomic mutations in either the coding or the regulatory space. It implies that mutated regulatory components of the genome contribute substantially to cancer pathways. Our analyses demonstrate that identifying genomically altered cis-regulatory elements coupled with analysis of gene expression data will augment biological interpretation of mutational landscapes of cancers.

An RCOR1 loss-associated gene expression signature identifies a prognostically significant DLBCL subgroup.

Effective treatment of diffuse large B-cell lymphoma (DLBCL) is plagued by heterogeneous responses to standard therapy, and molecular mechanisms underlying unfavorable outcomes in lymphoma patients remain elusive. Here, we profiled 148 genomes with 91 matching transcriptomes in a DLBCL cohort treated with rituximab plus cyclophosphamide, doxorubicin, vincristine, and prednisolone (R-CHOP) to uncover molecular subgroups linked to treatment failure. Systematic integration of high-resolution genotyping arrays and RNA sequencing data revealed novel deletions in RCOR1 to be associated with unfavorable progression-free survival (P = .001). Integration of expression data from the clinical samples with data from RCOR1 knockdowns in the lymphoma cell lines KM-H2 and Raji yielded an RCOR1 loss-associated gene signature comprising 233 genes. This signature identified a subgroup of patients with unfavorable overall survival (P = .023). The prognostic significance of the 233-gene signature for overall survival was reproduced in an independent cohort comprising 195 R-CHOP-treated patients (P = .039). Additionally, we discovered that within the International Prognostic Index low-risk group, the gene signature provides additional prognostic value that was independent of the cell-of-origin phenotype. We present a novel and reproducible molecular subgroup of DLBCL that impacts risk-stratification of R-CHOP-treated DLBCL patients and reveals a possible new avenue for therapeutic intervention strategies.

Intratumoral heterogeneity in a minority of ovarian low-grade serous carcinomas.

BACKGROUND: Ovarian low-grade serous carcinoma (LGSC) has fewer mutations than ovarian high-grade serous carcinoma (HGSC) and a less aggressive clinical course. However, an overwhelming majority of LGSC patients do not respond to conventional chemotherapy resulting in a poor long-term prognosis comparable to women diagnosed with HGSC. KRAS and BRAF mutations are common in LGSC, leading to clinical trials targeting the MAPK pathway. We assessed the stability of targetable somatic mutations over space and/or time in LGSC, with a view to inform stratified treatment strategies and clinical trial design. METHODS: Eleven LGSC cases with primary and recurrent paired samples were identified (stage IIB-IV). Tumor DNA was isolated from 1-4 formalin-fixed paraffin-embedded tumor blocks from both the primary and recurrence (n = 37 tumor and n = 7 normal samples). Mutational analysis was performed using the Ion Torrent AmpliSeqTM Cancer Panel, with targeted validation using Fluidigm-MiSeq, Sanger sequencing and/or Raindance Raindrop digital PCR. RESULTS: KRAS (3/11), BRAF (2/11) and/or NRAS (1/11) mutations were identified in five unique cases. A novel, non-synonymous mutation in SMAD4 was observed in one case. No somatic mutations were detected in the remaining six cases. In two cases with a single matched primary and recurrent sample, two KRAS hotspot mutations (G12V, G12R) were both stable over time. In three cases with multiple samplings from both the primary and recurrent surgery some mutations (NRAS Q61R, BRAF V600E, SMAD4 R361G) were stable across all samples, while others (KRAS G12V, BRAF G469V) were unstable. CONCLUSIONS: Overall, the majority of cases with detectable somatic mutations showed mutational stability over space and time while one of five cases showed both temporal and spatial mutational instability in presumed drivers of disease. Investigation of additional cases is required to confirm whether mutational heterogeneity in a minority of LGSC is a general phenomenon that should be factored into the design of clinical trials and stratified treatment for this patient population.

TITAN: inference of copy number architectures in clonal cell populations from tumor whole-genome sequence data.

The evolution of cancer genomes within a single tumor creates mixed cell populations with divergent somatic mutational landscapes. Inference of tumor subpopulations has been disproportionately focused on the assessment of somatic point mutations, whereas computational methods targeting evolutionary dynamics of copy number alterations (CNA) and loss of heterozygosity (LOH) in whole-genome sequencing data remain underdeveloped. We present a novel probabilistic model, TITAN, to infer CNA and LOH events while accounting for mixtures of cell populations, thereby estimating the proportion of cells harboring each event. We evaluate TITAN on idealized mixtures, simulating clonal populations from whole-genome sequences taken from genomically heterogeneous ovarian tumor sites collected from the same patient. In addition, we show in 23 whole genomes of breast tumors that the inference of CNA and LOH using TITAN critically informs population structure and the nature of the evolving cancer genome. Finally, we experimentally validated subclonal predictions using fluorescence in situ hybridization (FISH) and single-cell sequencing from an ovarian cancer patient sample, thereby recapitulating the key modeling assumptions of TITAN.

Ovarian and endometrial endometrioid carcinomas have distinct CTNNB1 and PTEN mutation profiles.

Ovarian endometrioid carcinomas and endometrial endometrioid carcinomas share many histological and molecular alterations. These similarities are likely due to a common endometrial epithelial precursor cell of origin, with most ovarian endometrioid carcinomas arising from endometriosis. To directly compare the mutation profiles of two morphologically similar tumor types, endometrial endometrioid carcinomas (n=307) and ovarian endometrioid carcinomas (n=33), we performed select exon capture sequencing on a panel of genes: ARID1A, PTEN, PIK3CA, KRAS, CTNNB1, PPP2R1A, TP53. We found that PTEN mutations are more frequent in low-grade endometrial endometrioid carcinomas (67%) compared with low-grade ovarian endometrioid carcinomas (17%) (P<0.0001). By contrast, CTNNB1 mutations are significantly different in low-grade ovarian endometrioid carcinomas (53%) compared with low-grade endometrial endometrioid carcinomas (28%) (P<0.0057). This difference in CTNNB1 mutation frequency may be reflective of the distinct microenvironments; the epithelial cells lining an endometriotic cyst within the ovary are exposed to a highly oxidative environment that promotes tumorigenesis. Understanding the distinct mutation patterns found in the PI3K and Wnt pathways of ovarian and endometrial endometrioid carcinomas may provide future opportunities for stratifying patients for targeted therapeutics.

Distinct evolutionary trajectories of primary high-grade serous ovarian cancers revealed through spatial mutational profiling.

High-grade serous ovarian cancer (HGSC) is characterized by poor outcome, often attributed to the emergence of treatment-resistant subclones. We sought to measure the degree of genomic diversity within primary, untreated HGSCs to examine the natural state of tumour evolution prior to therapy. We performed exome sequencing, copy number analysis, targeted amplicon deep sequencing and gene expression profiling on 31 spatially and temporally separated HGSC tumour specimens (six patients), including ovarian masses, distant metastases and fallopian tube lesions. We found widespread intratumoural variation in mutation, copy number and gene expression profiles, with key driver alterations in genes present in only a subset of samples (eg PIK3CA, CTNNB1, NF1). On average, only 51.5% of mutations were present in every sample of a given case (range 10.2-91.4%), with TP53 as the only somatic mutation consistently present in all samples. Complex segmental aneuploidies, such as whole-genome doubling, were present in a subset of samples from the same individual, with divergent copy number changes segregating independently of point mutation acquisition. Reconstruction of evolutionary histories showed one patient with mixed HGSC and endometrioid histology, with common aetiologic origin in the fallopian tube and subsequent selection of different driver mutations in the histologically distinct samples. In this patient, we observed mixed cell populations in the early fallopian tube lesion, indicating that diversity arises at early stages of tumourigenesis. Our results revealed that HGSCs exhibit highly individual evolutionary trajectories and diverse genomic tapestries prior to therapy, exposing an essential biological characteristic to inform future design of personalized therapeutic solutions and investigation of drug-resistance mechanisms.