Search

Single-cell sequencing technologies have revolutionized biomedical research by enabling deconvolution of cell type-specific properties in highly heterogeneous tissue. While robust tools have been developed to handle bioinformatic challenges posed by single-cell RNA and ATAC data, options for emergent modalities such as methylation are much more limited, impeding the utility of results. Here we present Amethyst, a comprehensive R package for atlas-scale single-cell methylation sequencing data analysis. Amethyst begins with base-level methylation calls and expedites batch integration, doublet detection, dimensionality reduction, clustering, cell type annotation, differentially methylated region calling, and interpretation of results, facilitating rapid data interaction in a local environment. We introduce the workflow using published single-cell methylation human peripheral blood mononuclear cell (PBMC) and human cortex data. We further leverage Amethyst on an atlas-scale brain dataset to describe a noncanonical methylation pattern in human astrocytes and oligodendrocytes, challenging the notion that this form of methylation is principally relevant to neurons in the brain. Tools such as Amethyst will increase accessibility to single-cell methylation data analysis, catalyzing research progress across diverse contexts.

Mapping the spatial interactions of cancer, immune, and stromal cell states presents novel opportunities for patient stratification and for advancing immunotherapy. While single-cell studies revealed significant molecular heterogeneity in prostate cancer cells, the impact of spatial stromal cell heterogeneity remains poorly understood. Here, we used cyclic immunofluorescent imaging on whole-tissue sections to uncover novel spatial associations between cancer and stromal cells in low- and high-grade prostate tumors and tumor-adjacent normal tissues. Our results provide a spatial map of single cells and recurrent cellular neighborhoods in the prostate tumor microenvironment of treatment-naive patients. We report unique populations of mast cells that show distinct spatial associations with M2 macrophages and regulatory T cells. Our results show disease-specific neighborhoods that are primarily driven by androgen receptor-positive (AR+) stromal cells and identify inflammatory gene networks active in AR+ prostate stroma., Competing Interests: The authors declare no competing interests., (© 2024 The Authors.)

DNA methylation is a key component of the mammalian epigenome, playing a regulatory role in development, disease, and other processes. Robust, high-throughput single-cell DNA methylation assays are now possible (sciMET); however, the genome-wide nature of DNA methylation results in a high sequencing burden per cell. Here, we leverage target enrichment with sciMET to capture sufficient information per cell for cell type assignment using substantially fewer sequence reads (sciMET-cap). Accumulated off-target coverage enables genome-wide differentially methylated region (DMR) calling for clusters with as few as 115 cells. We characterize sciMET-cap on human PBMCs and brain (middle frontal gyrus)., (© 2024. The Author(s).)

Abstract: Identifying and targeting microenvironment-driven pathways that are active across acute myeloid leukemia (AML) genetic subtypes should allow the development of more broadly effective therapies. The proinflammatory cytokine interleukin-1β (IL-1β) is abundant in the AML microenvironment and promotes leukemic growth. Through RNA-sequencing analysis, we identify that IL-1β-upregulated ASF1B (antisilencing function-1B), a histone chaperone, in AML progenitors compared with healthy progenitors. ASF1B, along with its paralogous protein ASF1A, recruits H3-H4 histones onto the replication fork during S-phase, a process regulated by Tousled-like kinase 1 and 2 (TLKs). Although ASF1s and TLKs are known to be overexpressed in multiple solid tumors and associated with poor prognosis, their functional roles in hematopoiesis and inflammation-driven leukemia remain unexplored. In this study, we identify that ASF1s and TLKs are overexpressed in multiple genetic subtypes of AML. We demonstrate that depletion of ASF1s significantly reduces leukemic cell growth in both in vitro and in vivo models using human cells. Using a murine model, we show that overexpression of ASF1B accelerates leukemia progression. Moreover, Asf1b or Tlk2 deletion delayed leukemia progression, whereas these proteins are dispensable for normal hematopoiesis. Through proteomics and phosphoproteomics analyses, we uncover that the TLK-ASF1 pathway promotes leukemogenesis by affecting the cell cycle and DNA damage pathways. Collectively, our findings identify the TLK1-ASF1 pathway as a novel mediator of inflammatory signaling and a promising therapeutic target for AML treatment across diverse genetic subtypes. Selective inhibition of this pathway offers potential opportunities to intervene effectively, address intratumoral heterogeneity, and ultimately improve clinical outcomes in AML.

Single-cell assay for transposase-accessible chromatin by sequencing (scATAC-seq) has emerged as a powerful tool for dissecting regulatory landscapes and cellular heterogeneity. However, an exploration of systemic biases among scATAC-seq technologies has remained absent. In this study, we benchmark the performance of eight scATAC-seq methods across 47 experiments using human peripheral blood mononuclear cells (PBMCs) as a reference sample and develop PUMATAC, a universal preprocessing pipeline, to handle the various sequencing data formats. Our analyses reveal significant differences in sequencing library complexity and tagmentation specificity, which impact cell-type annotation, genotype demultiplexing, peak calling, differential region accessibility and transcription factor motif enrichment. Our findings underscore the importance of sample extraction, method selection, data processing and total cost of experiments, offering valuable guidance for future research. Finally, our data and analysis pipeline encompasses 169,000 PBMC scATAC-seq profiles and a best practices code repository for scATAC-seq data analysis, which are freely available to extend this benchmarking effort to future protocols., (© 2023. The Author(s).)

We develop a large-scale single-cell ATAC-seq method by combining Tn5-based pre-indexing with 10× Genomics barcoding, enabling the indexing of up to 200,000 nuclei across multiple samples in a single reaction. We profile 449,953 nuclei across diverse tissues, including the human cortex, mouse brain, human lung, mouse lung, mouse liver, and lung tissue from a club cell secretory protein knockout (CC16 -/- ) model. Our study of CC16 -/- nuclei uncovers previously underappreciated technical artifacts derived from remnant 129 mouse strain genetic material, which cause profound cell-type-specific changes in regulatory elements near many genes, thereby confounding the interpretation of this commonly referenced mouse model., (© 2024. The Author(s).)

Background: Cancer health research relies on large-scale cohorts to derive generalizable results for different populations. While traditional epidemiological cohorts often use costly random sampling or self-motivated, preselected groups, a shift toward health system-based cohorts has emerged. However, such cohorts depend on participants remaining within a single system. Recent consumer engagement models using smartphone-based communication, driving projects, and social media have begun to upend these paradigms., Methods: We initiated the Healthy Oregon Project (HOP) to support basic and clinical cancer research. HOP study employs a novel, cost-effective remote recruitment approach to effectively establish a large-scale cohort for population-based studies. The recruitment leverages the unique email account, the HOP website, and social media platforms to direct smartphone users to the study app, which facilitates saliva sample collection and survey administration. Monthly newsletters further facilitate engagement and outreach to broader communities., Results: By the end of 2022, the HOP has enrolled approximately 35,000 participants aged 18-100 years (median = 44.2 years), comprising more than 1% of the Oregon adult population. Among those who have app access, ∼87% provided consent to genetic screening. The HOP monthly email newsletters have an average open rate of 38%. Efforts continue to be made to improve survey response rates., Conclusion: This study underscores the efficacy of remote recruitment approaches in establishing large-scale cohorts for population-based cancer studies. The implementation of the study facilitates the collection of extensive survey and biological data into a repository that can be broadly shared and supports collaborative clinical and translational research., Competing Interests: None., (© The Author(s) 2024.)

The Healthy Oregon Project (HOP) is a statewide effort that aims to build a large research repository and influence the health of Oregonians through providing no-cost genetic screening to participants for a next-generation sequencing 32-gene panel comprising genes related to inherited cancers and familial hypercholesterolemia. This type of unbiased population screening can detect at-risk individuals who may otherwise be missed by conventional medical approaches. However, challenges exist for this type of high-throughput testing in an academic setting, including developing a low-cost high-efficiency test and scaling up the clinical laboratory for processing large numbers of samples. Modifications to our academic clinical laboratory including efficient test design, robotics, and a streamlined analysis approach increased our ability to test more than 1,000 samples per month for HOP using only one dedicated HOP laboratory technologist. Additionally, enrollment using a HIPAA-compliant smartphone app and sample collection using mouthwash increased efficiency and reduced cost. Here, we present our experience three years into HOP and discuss the lessons learned, including our successes, challenges, opportunities, and future directions, as well as the genetic screening results for the first 13,670 participants tested. Overall, we have identified 730 pathogenic/likely pathogenic variants in 710 participants in 24 of the 32 genes on the panel. The carrier rate for pathogenic/likely pathogenic variants in the inherited cancer genes on the panel for an unselected population was 5.0% and for familial hypercholesterolemia was 0.3%. Our laboratory experience described here may provide a useful model for population screening projects in other states., Competing Interests: Declaration of interests C.S.R. is a Director on the American Board of Genetics and Genomics. P.T.S. has the following conflicts: employee of OHSU, consultant to Natera, Twinstrand, Foresight Diagnostics, and Lab Corp, and equity position in Convergent Genomics. T.D.O., A.B.P., C.C.D., G.G., J.H.L., S. McCabe, A.K., R.W., J.T., and C.S.R. are employees of the Knight Diagnostic Laboratories (KDL)., (Copyright © 2023 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Many acute myeloid leukemia (AML) patients exhibit hallmarks of immune exhaustion, such as increased myeloid-derived suppressor cells, suppressive regulatory T cells and dysfunctional T cells. Similarly, we have identified the same immune-related features, including exhausted CD8 + T cells (TEx) in a mouse model of AML. Here we show that inhibitors that target bromodomain and extra-terminal domain (BET) proteins affect tumor-intrinsic factors but also rescue T cell exhaustion and ICB resistance. Ex vivo treatment of cells from AML mice and AML patients with BET inhibitors (BETi) reversed CD8 + T cell exhaustion by restoring proliferative capacity and expansion of the more functional precursor-exhausted T cells. This reversal was enhanced by combined BETi and anti-PD1 treatment. BETi synergized with anti-PD1 in vivo, resulting in the reduction of circulating leukemia cells, enrichment of CD8 + T cells in the bone marrow, and increase in expression of Tcf7, Slamf6, and Cxcr5 in CD8 + T cells. Finally, we profiled the epigenomes of in vivo JQ1-treated AML-derived CD8 + T cells by single-cell ATAC-seq and found that JQ1 increases Tcf7 accessibility specifically in Tex cells, suggesting that BETi likely acts mechanistically by relieving repression of progenitor programs in Tex CD8 + T cells and maintaining a pool of anti-PD1 responsive CD8 + T cells., (© 2023. The Author(s).)

DNA methylation is a key epigenetic property that drives gene regulatory programs in development and disease. Current single-cell methods that produce high quality methylomes are expensive and low throughput without the aid of extensive automation. We previously described a proof-of-principle technique that enabled high cell throughput; however, it produced only low-coverage profiles and was a difficult protocol that required custom sequencing primers and recipes and frequently produced libraries with excessive adapter contamination. Here, we describe a greatly improved version that generates high-coverage profiles (~15-fold increase) using a robust protocol that does not require custom sequencing capabilities, includes multiple stopping points, and exhibits minimal adapter contamination. We demonstrate two versions of sciMETv2 on primary human cortex, a high coverage and rapid version, identifying distinct cell types using CH methylation patterns. These datasets are able to be directly integrated with one another as well as with existing snmC-seq2 datasets with little discernible bias. Finally, we demonstrate the ability to determine cell types using CG methylation alone, which is the dominant context for DNA methylation in most cell types other than neurons and the most applicable analysis outside of brain tissue., (© 2022. The Author(s).)

Single-cell RNA sequencing (scRNA-seq) distinguishes cell types, states and lineages within the context of heterogeneous tissues. However, current single-cell data cannot directly link cell clusters with specific phenotypes. Here we present Scissor, a method that identifies cell subpopulations from single-cell data that are associated with a given phenotype. Scissor integrates phenotype-associated bulk expression data and single-cell data by first quantifying the similarity between each single cell and each bulk sample. It then optimizes a regression model on the correlation matrix with the sample phenotype to identify relevant subpopulations. Applied to a lung cancer scRNA-seq dataset, Scissor identified subsets of cells associated with worse survival and with TP53 mutations. In melanoma, Scissor discerned a T cell subpopulation with low PDCD1/CTLA4 and high TCF7 expression associated with an immunotherapy response. Beyond cancer, Scissor was effective in interpreting facioscapulohumeral muscular dystrophy and Alzheimer's disease datasets. Scissor identifies biologically and clinically relevant cell subpopulations from single-cell assays by leveraging phenotype and bulk-omics datasets., (© 2021. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Identifying precise molecular subtypes attributable to specific stages of localized prostate cancer has proven difficult due to high levels of heterogeneity. Bulk assays represent a population-average, which mask the heterogeneity that exists at the single-cell level. In this work, we sequence the accessible chromatin regions of 14,424 single-cells from 18 flash-frozen prostate tumours. We observe shared chromatin features among low-grade prostate cancer cells are lost in high-grade tumours. Despite this loss, high-grade tumours exhibit an enrichment for FOXA1, HOXB13 and CDX2 transcription factor binding sites, indicating a shared trans-regulatory programme. We identify two unique genes encoding neuronal adhesion molecules that are highly accessible in high-grade prostate tumours. We show NRXN1 and NLGN1 expression in epithelial, endothelial, immune and neuronal cells in prostate cancer using cyclic immunofluorescence. Our results provide a deeper understanding of the active gene regulatory networks in primary prostate tumours, critical for molecular stratification of the disease., (© 2021. The Author(s).)

Single-cell combinatorial indexing (sci) with transposase-based library construction increases the throughput of single-cell genomics assays but produces sparse coverage in terms of usable reads per cell. We develop symmetrical strand sci ('s3'), a uracil-based adapter switching approach that improves the rate of conversion of source DNA into viable sequencing library fragments following tagmentation. We apply this chemistry to assay chromatin accessibility (s3-assay for transposase-accessible chromatin, s3-ATAC) in human cortical and mouse whole-brain tissues, with mouse datasets demonstrating a six- to 13-fold improvement in usable reads per cell compared with other available methods. Application of s3 to single-cell whole-genome sequencing (s3-WGS) and to whole-genome plus chromatin conformation (s3-GCC) yields 148- and 14.8-fold improvements, respectively, in usable reads per cell compared with sci-DNA-sequencing and sci-HiC. We show that s3-WGS and s3-GCC resolve subclonal genomic alterations in patient-derived pancreatic cancer cell lines. We expect that the s3 platform will be compatible with other transposase-based techniques, including sci-MET or CUT&Tag., (© 2021. The Author(s), under exclusive licence to Springer Nature America, Inc.)

SARS-CoV-2 and its variants continue to infect hundreds of thousands every day despite the rollout of effective vaccines. Therefore, it is essential to understand the levels of protection that these vaccines provide in the face of emerging variants. Here, we report two demographically balanced cohorts of BNT162b2 vaccine recipients and COVID-19 patients, from which we evaluate neutralizing antibody titers against SARS-CoV-2 as well as the B.1.1.7 (alpha) and B.1.351 (beta) variants. We show that both B.1.1.7 and B.1.351 are less well neutralized by serum from vaccinated individuals, and that B.1.351, but not B.1.1.7, is less well neutralized by convalescent serum. We also find that the levels of variant-specific anti-spike antibodies are proportional to neutralizing activities. Together, our results demonstrate the escape of the emerging SARS-CoV-2 variants from neutralization by serum antibodies, which may lead to reduced protection from re-infection or increased risk of vaccine breakthrough., (© 2021. The Author(s).)

High-throughput single-cell epigenomic assays can resolve cell type heterogeneity in complex tissues, however, spatial orientation is lost. Here, we present single-cell combinatorial indexing on Microbiopsies Assigned to Positions for the Assay for Transposase Accessible Chromatin, or sciMAP-ATAC, as a method for highly scalable, spatially resolved, single-cell profiling of chromatin states. sciMAP-ATAC produces data of equivalent quality to non-spatial sci-ATAC and retains the positional information of each cell within a 214 micron cubic region, with up to hundreds of tracked positions in a single experiment. We apply sciMAP-ATAC to assess cortical lamination in the adult mouse primary somatosensory cortex and in the human primary visual cortex, where we produce spatial trajectories and integrate our data with non-spatial single-nucleus RNA and other chromatin accessibility single-cell datasets. Finally, we characterize the spatially progressive nature of cerebral ischemic infarction in the mouse brain using a model of transient middle cerebral artery occlusion.

Intratumoral heterogeneity in cancers arises from genomic instability and epigenomic plasticity and is associated with resistance to cytotoxic and targeted therapies. We show here that cell-state heterogeneity, defined by differentiation-state marker expression, is high in triple-negative and basal-like breast cancer subtypes, and that drug tolerant persister (DTP) cell populations with altered marker expression emerge during treatment with a wide range of pathway-targeted therapeutic compounds. We show that MEK and PI3K/mTOR inhibitor-driven DTP states arise through distinct cell-state transitions rather than by Darwinian selection of preexisting subpopulations, and that these transitions involve dynamic remodeling of open chromatin architecture. Increased activity of many chromatin modifier enzymes, including BRD4, is observed in DTP cells. Co-treatment with the PI3K/mTOR inhibitor BEZ235 and the BET inhibitor JQ1 prevents changes to the open chromatin architecture, inhibits the acquisition of a DTP state, and results in robust cell death in vitro and xenograft regression in vivo.

Aberrations of the DNA methylome contribute to onset and progression of diseases. Whole genome bisulfite sequencing (WGBS) is the only analytical method covering the complete methylome. Alternative methods requiring less DNA than WGBS analyze only a minor portion of the methylome and do not cover important regulatory features like enhancers and noncoding RNAs. In tagmentation-based WGBS (TWGBS), several DNA and time-consuming steps of the conventional WGBS library preparation are circumvented by the use of a hyperactive transposase, which simultaneously fragments DNA and appends sequencing adapters. TWGBS requires only nanogram amounts of DNA and, thus, is well suited to study precious biological specimens such as sorted cells or micro-dissected tissue samples.