Haplotype-aware analysis of somatic copy number variations from single-cell transcriptomes

Somatic mutations associate with clonal expansion of CD8+ T cells

Somatic mutations in T cells can cause cancer but also have implications for immunological diseases and cell therapies. The mutation spectrum in nonmalignant T cells is unclear. Here, we examined somatic mutations in CD4+ and CD8+ T cells from 90 patients with hematological and immunological disorders and used T cell receptor (TCR) and single-cell sequencing to link mutations with T cell expansions and phenotypes. CD8+ cells had a higher mutation burden than CD4+ cells. Notably, the biggest variant allele frequency (VAF) of non-synonymous variants was higher than synonymous variants in CD8+ T cells, indicating non-random occurrence. The non-synonymous VAF in CD8+ T cells strongly correlated with the TCR frequency, but not age. We identified mutations in pathways essential for T cell function and often affected lymphoid neoplasia. Single-cell sequencing revealed cytotoxic TEMRA phenotypes of mutated T cells. Our findings suggest that somatic mutations contribute to CD8+ T cell expansions without malignant transformation.

HATCHet2: clone- and haplotype-specific copy number inference from bulk tumor sequencing data

Bulk DNA sequencing of multiple samples from the same tumor is becoming common, yet most methods to infer copy-number aberrations (CNAs) from this data analyze individual samples independently. We introduce HATCHet2, an algorithm to identify haplotype- and clone-specific CNAs simultaneously from multiple bulk samples. HATCHet2 extends the earlier HATCHet method by improving identification of focal CNAs and introducing a novel statistic, the minor haplotype B-allele frequency (mhBAF), that enables identification of mirrored-subclonal CNAs. We demonstrate HATCHet2's improved accuracy using simulations and a single-cell sequencing dataset. HATCHet2 analysis of 10 prostate cancer patients reveals previously unreported mirrored-subclonal CNAs affecting cancer genes.

Copy number alterations: a catastrophic orchestration of the breast cancer genome

Breast cancer (BCa) is a prevalent malignancy that predominantly affects women around the world. Somatic copy number alterations (CNAs) are tumor-specific amplifications or deletions of DNA segments that often drive BCa development and therapy resistance. Hence, the complex patterns of CNAs complement BCa classification systems. In addition, understanding the precise contributions of CNAs is essential for tailoring personalized treatment approaches. This review highlights how tumor evolution drives the acquisition of CNAs, which in turn shape the genomic landscapes of BCas. It also discusses advanced methodologies for identifying recurrent CNAs, studying CNAs in BCa and their clinical impact.

Integrated electrophysiological and genomic profiles of single cells reveal spiking tumor cells in human glioma

Prior studies have described the complex interplay that exists between glioma cells and neurons, however, the electrophysiological properties endogenous to tumor cells remain obscure. To address this, we employed Patch-sequencing on human glioma specimens and found that one third of patched cells in IDH mutant (IDH mut ) tumors demonstrate properties of both neurons and glia by firing single, short action potentials. To define these hybrid cells (HCs) and discern if they are tumor in origin, we developed a computational tool, Single Cell Rule Association Mining (SCRAM), to annotate each cell individually. SCRAM revealed that HCs represent tumor and non-tumor cells that feature GABAergic neuron and oligodendrocyte precursor cell signatures. These studies are the first to characterize the combined electrophysiological and molecular properties of human glioma cells and describe a new cell type in human glioma with unique electrophysiological and transcriptomic properties that are likely also present in the non-tumor mammalian brain.

Phasing millions of samples achieves near perfect accuracy, enabling parent-of-origin classification of variants

Haplotype phasing, the process of determining which genetic variants are physically located on the same chromosome, is crucial for various genetic analyses. In this study, we first benchmark SHAPEIT and Beagle, two state-of-the-art phasing methods, on two large datasets: > 8 million diverse, research-consented 23andMe, Inc. customers and the UK Biobank (UKB). We find that both perform exceptionally well. Beagle's median switch error rate (SER) (after excluding single SNP switches) in white British trios from UKB is 0.026% compared to 0.00% for European ancestry 23andMe research participants; 55.6% of European ancestry 23andMe research participants have zero non-single SNP switches, compared to 42.4% of white British trios. South Asian ancestry 23andMe research participants have the highest median SER amongst the 23andMe populations, but it is still remarkably low at 0.46%. We also investigate the relationship between identity-by-descent (IBD) and SER, finding that switch errors tend to occur in regions of little or no IBD segment coverage. SHAPEIT and Beagle excel at 'intra-chromosomal' phasing, but lack the ability to phase across chromosomes, motivating us to develop an inter-chromosomal phasing method, called HAPTIC ( HAP lotype TI ling and C lustering), that assigns paternal and maternal variants discretely genome-wide. Our approach uses identity-by-descent (IBD) segments to phase blocks of variants on different chromosomes. HAPTIC represents the segments a focal individual shares with their relatives as nodes in a signed graph and performs bipartite clustering on the signed graph using spectral clustering. We test HAPTIC on 1022 UKB trios, yielding a median phase error of 0.08% in regions covered by IBD segments (33.5% of sites). We also ran HAPTIC in the 23andMe database and found a median phase error rate (the rate of mismatching alleles between the inferred and true phase) of 0.92% in Europeans (93.8% of sites) and 0.09% in admixed Africans (92.7% of sites). HAPTIC's precision depends heavily on data from relatives, so will increase as datasets grow larger and more diverse. HAPTIC enables analyses that require the parent-of-origin of variants, such as association studies and ancestry inference of untyped parents.

CopyVAE: a variational autoencoder-based approach for copy number variation inference using single-cell transcriptomics

MOTIVATION

Copy number variations (CNVs) are common genetic alterations in tumour cells. The delineation of CNVs holds promise for enhancing our comprehension of cancer progression. Moreover, accurate inference of CNVs from single-cell sequencing data is essential for unravelling intratumoral heterogeneity. However, existing inference methods face limitations in resolution and sensitivity.

RESULTS

To address these challenges, we present CopyVAE, a deep learning framework based on a variational autoencoder architecture. Through experiments, we demonstrated that CopyVAE can accurately and reliably detect CNVs from data obtained using single-cell RNA sequencing. CopyVAE surpasses existing methods in terms of sensitivity and specificity. We also discussed CopyVAE's potential to advance our understanding of genetic alterations and their impact on disease advancement.

AVAILABILITY AND IMPLEMENTATION

CopyVAE is implemented and freely available under MIT license at https://github.com/kurtsemih/copyVAE.

Early Immune Remodeling Steers Clinical Response to First-Line Chemoimmunotherapy in Advanced Gastric Cancer

Adding anti-programmed cell death protein 1 (anti-PD-1) to 5-fluorouracil (5-FU)/platinum improves survival in some advanced gastroesophageal adenocarcinomas (GEA). To understand the effects of chemotherapy and immunotherapy, we conducted a phase II first-line trial (n = 47) sequentially adding pembrolizumab to 5-FU/platinum in advanced GEA. Using serial biopsy of the primary tumor at baseline, after one cycle of 5-FU/platinum, and after the addition of pembrolizumab, we transcriptionally profiled 358,067 single cells to identify evolving multicellular tumor microenvironment (TME) networks. Chemotherapy induced early on-treatment multicellular hubs with tumor-reactive T-cell and M1-like macrophage interactions in slow progressors. Faster progression featured increased MUC5A and MSLN containing treatment resistance programs in tumor cells and M2-like macrophages with immunosuppressive stromal interactions. After pembrolizumab, we observed increased CD8 T-cell infiltration and development of an immunity hub involving tumor-reactive CXCL13 T-cell program and epithelial interferon-stimulated gene programs. Strategies to drive increases in antitumor immune hub formation could expand the portion of patients benefiting from anti-PD-1 approaches.

SIGNIFICANCE

The benefit of 5-FU/platinum with anti-PD-1 in first-line advanced gastric cancer is limited to patient subgroups. Using a trial with sequential anti-PD-1, we show coordinated induction of multicellular TME hubs informs the ability of anti-PD-1 to potentiate T cell-driven responses. Differential TME hub development highlights features that underlie clinical outcomes. This article is featured in Selected Articles from This Issue, p. 695.

De novo detection of somatic mutations in high-throughput single-cell profiling data sets

Characterization of somatic mutations at single-cell resolution is essential to study cancer evolution, clonal mosaicism and cell plasticity. Here, we describe SComatic, an algorithm designed for the detection of somatic mutations in single-cell transcriptomic and ATAC-seq (assay for transposase-accessible chromatin sequence) data sets directly without requiring matched bulk or single-cell DNA sequencing data. SComatic distinguishes somatic mutations from polymorphisms, RNA-editing events and artefacts using filters and statistical tests parameterized on non-neoplastic samples. Using >2.6 million single cells from 688 single-cell RNA-seq (scRNA-seq) and single-cell ATAC-seq (scATAC-seq) data sets spanning cancer and non-neoplastic samples, we show that SComatic detects mutations in single cells accurately, even in differentiated cells from polyclonal tissues that are not amenable to mutation detection using existing methods. Validated against matched genome sequencing and scRNA-seq data, SComatic achieves F1 scores between 0.6 and 0.7 across diverse data sets, in comparison to 0.2-0.4 for the second-best performing method. In summary, SComatic permits de novo mutational signature analysis, and the study of clonal heterogeneity and mutational burdens at single-cell resolution.

Long-Term Follow-Up Defines the Population That Benefits from Early Interception in a High-Risk Smoldering Multiple Myeloma Clinical Trial Using the Combination of Ixazomib, Lenalidomide, and Dexamethasone

Background

Early therapeutic intervention in high-risk SMM (HR-SMM) has demonstrated benefit in previous studies of lenalidomide with or without dexamethasone. Triplets and quadruplet studies have been examined in this same population. However, to date, none of these studies examined the impact of depth of response on long-term outcomes of participants treated with lenalidomide-based therapy, and whether the use of the 20/2/20 model or the addition of genomic alterations can further define the population that would benefit the most from early therapeutic intervention. Here, we present the results of the phase II study of the combination of ixazomib, lenalidomide, and dexamethasone in patients with HR-SMM with long-term follow-up and baseline single-cell tumor and immune sequencing that help refine the population to be treated for early intervention studies.

Methods

This is a phase II trial of ixazomib, lenalidomide, and dexamethasone (IRD) in HR-SMM. Patients received 9 cycles of induction therapy with ixazomib 4mg on days 1, 8, and 15; lenalidomide 25mg on days 1-21; and dexamethasone 40mg on days 1, 8, 15, and 22. The induction phase was followed by maintenance with ixazomib 4mg on days 1, 8, and 15; and lenalidomide 15mg d1-21 for 15 cycles for 24 months of treatment. The primary endpoint was progression-free survival after 2 years of therapy. Secondary endpoints included depth of response, biochemical progression, and correlative studies included single-cell RNA sequencing and/or whole-genome sequencing of the tumor and single-cell sequencing of immune cells at baseline.

Results

Fifty-five patients, with a median age of 64, were enrolled in the study. The overall response rate was 93%, with 31% of patients achieving a complete response and 45% achieving a very good partial response or better. The most common grade 3 or greater treatment-related hematologic toxicities were neutropenia (16 patients; 29%), leukopenia (10 patients; 18%), lymphocytopenia (8 patients; 15%), and thrombocytopenia (4 patients; 7%). Non-hematologic grade 3 or greater toxicities included hypophosphatemia (7 patients; 13%), rash (5 patients; 9%), and hypokalemia (4 patients; 7%). After a median follow-up of 50 months, the median progression-free survival (PFS) was 48.6 months (95% CI: 39.9 - not reached; NR) and median overall survival has not been reached. Patients achieving VGPR or better had a significantly better progression-free survival (p<0.001) compared to those who did not achieve VGPR (median PFS 58.2 months vs. 31.3 months). Biochemical progression preceded or was concurrent with the development of SLiM-CRAB criteria in eight patients during follow-up, indicating that biochemical progression is a meaningful endpoint that correlates with the development of end-organ damage. High-risk 20/2/20 participants had the worst PFS compared to low- and intermediate-risk participants. The use of whole genome or single-cell sequencing of tumor cells identified high-risk aberrations that were not identified by FISH alone and aided in the identification of participants at risk of progression. scRNA-seq analysis revealed a positive correlation between MHC class I expression and response to proteasome inhibition and at the same time a decreased proportion of GZMB+ T cells within the clonally expanded CD8+ T cell population correlated with suboptimal response.

Conclusions

Ixazomib, lenalidomide and dexamethasone in HR-SMM demonstrates significant clinical activity with an overall favorable safety profile. Achievement of VGPR or greater led to significant improvement in time to progression, suggesting that achieving deep response is beneficial in HR-SMM. Biochemical progression correlates with end-organ damage. Patients with high-risk FISH and lack of deep response had poor outcomes. ClinicalTrials.gov identifier: ( NCT02916771 ).

Investigation of cellular communication and signaling pathways in tumor microenvironment for high TP53-expressing osteosarcoma cells through single-cell RNA sequencing

Osteosarcoma (OS) stands as the most prevalent primary bone cancer in children and adolescents, and its limited treatment options often result in unsatisfactory outcomes, particularly for metastatic cases. The tumor microenvironment (TME) has been recognized as a crucial determinant in OS progression. However, the intercellular dynamics between high TP53-expressing OS cells and neighboring cell types within the TME are yet to be thoroughly understood. In our study, we harnessed the single-cell RNA sequencing (scRNA-seq) technology in combination with the computational tool-Cellchat, aiming to elucidate the intercellular communication networks present within OS. Through meticulous quantitative inference and subsequent analysis of these networks, we succeeded in identifying significant signaling pathways connecting high TP53-expressing OS cells with proximate cell types, namely Macrophages, Monocytes, Endothelial Cells, and PVLs. This research brings forth a nuanced understanding of the intricate patterns and coordination involved in the TME's intercellular communication signals. These findings not only provide profound insights into the molecular mechanisms underpinning OS but also indicate potential therapeutic targets that could revolutionize treatment strategies.

Allele-specific transcriptional effects of subclonal copy number alterations enable genotype-phenotype mapping in cancer cells

Subclonal copy number alterations are a prevalent feature in tumors with high chromosomal instability and result in heterogeneous cancer cell populations with distinct phenotypes. However, the extent to which subclonal copy number alterations contribute to clone-specific phenotypes remains poorly understood. We develop TreeAlign, which computationally integrates independently sampled single-cell DNA and RNA sequencing data from the same cell population. TreeAlign accurately encodes dosage effects from subclonal copy number alterations, the impact of allelic imbalance on allele-specific transcription, and obviates the need to define genotypic clones from a phylogeny a priori, leading to highly granular definitions of clones with distinct expression programs. These improvements enable clone-clone gene expression comparisons with higher resolution and identification of expression programs that are genomically independent. Our approach sets the stage for dissecting the relative contribution of fixed genomic alterations and dynamic epigenetic processes on gene expression programs in cancer.

De novo detection of somatic variants in long-read single-cell RNA sequencing data

In cancer, genetic and transcriptomic variations generate clonal heterogeneity, possibly leading to treatment resistance. Long-read single-cell RNA sequencing (LR scRNA-seq) has the potential to detect genetic and transcriptomic variations simultaneously. Here, we present LongSom, a computational workflow leveraging LR scRNA-seq data to call de novo somatic single-nucleotide variants (SNVs), copy-number alterations (CNAs), and gene fusions to reconstruct the tumor clonal heterogeneity. For SNV calling, LongSom distinguishes somatic SNVs from germline polymorphisms by reannotating marker gene expression-based cell types using called variants and applying strict filters. Applying LongSom to ovarian cancer samples, we detected clinically relevant somatic SNVs that were validated against single-cell and bulk panel DNA-seq data and could not be detected with short-read (SR) scRNA-seq. Leveraging somatic SNVs and fusions, LongSom found subclones with different predicted treatment outcomes. In summary, LongSom enables de novo SNVs, CNAs, and fusions detection, thus enabling the study of cancer evolution, clonal heterogeneity, and treatment resistance.

Coordinated Immune Cell Networks in the Bone Marrow Microenvironment Define the Graft versus Leukemia Response with Adoptive Cellular Therapy

Understanding how intra-tumoral immune populations coordinate to generate anti-tumor responses following therapy can guide precise treatment prioritization. We performed systematic dissection of an established adoptive cellular therapy, donor lymphocyte infusion (DLI), by analyzing 348,905 single-cell transcriptomes from 74 longitudinal bone-marrow samples of 25 patients with relapsed myeloid leukemia; a subset was evaluated by protein-based spatial analysis. In acute myelogenous leukemia (AML) responders, diverse immune cell types within the bone-marrow microenvironment (BME) were predicted to interact with a clonally expanded population of ZNF683 + GZMB + CD8+ cytotoxic T lymphocytes (CTLs) which demonstrated in vitro specificity for autologous leukemia. This population, originating predominantly from the DLI product, expanded concurrently with NK and B cells. AML nonresponder BME revealed a paucity of crosstalk and elevated TIGIT expression in CD8+ CTLs. Our study highlights recipient BME differences as a key determinant of effective anti-leukemia response and opens new opportunities to modulate cell-based leukemia-directed therapy.

Somatic cancer driver mutations are enriched and associated with inflammatory states in Alzheimer's disease microglia

Alzheimer's disease (AD) is an age-associated neurodegenerative disorder characterized by progressive neuronal loss and pathological accumulation of the misfolded proteins amyloid-β and tau1,2. Neuroinflammation mediated by microglia and brain-resident macrophages plays a crucial role in AD pathogenesis1-5, though the mechanisms by which age, genes, and other risk factors interact remain largely unknown. Somatic mutations accumulate with age and lead to clonal expansion of many cell types, contributing to cancer and many non-cancer diseases6,7. Here we studied somatic mutation in normal aged and AD brains by three orthogonal methods and in three independent AD cohorts. Analysis of bulk RNA sequencing data from 866 samples from different brain regions revealed significantly higher (~two-fold) overall burdens of somatic single-nucleotide variants (sSNVs) in AD brains compared to age-matched controls. Molecular-barcoded deep (>1000X) gene panel sequencing of 311 prefrontal cortex samples showed enrichment of sSNVs and somatic insertions and deletions (sIndels) in cancer driver genes in AD brain compared to control, with recurrent, and often multiple, mutations in genes implicated in clonal hematopoiesis (CH)8,9. Pathogenic sSNVs were enriched in CSF1R+ microglia of AD brains, and the high proportion of microglia (up to 40%) carrying some sSNVs in cancer driver genes suggests mutation-driven microglial clonal expansion (MiCE). Analysis of single-nucleus RNA sequencing (snRNAseq) from temporal neocortex of 62 additional AD cases and controls exhibited nominally increased mosaic chromosomal alterations (mCAs) associated with CH10,11. Microglia carrying mCA showed upregulated pro-inflammatory genes, resembling the transcriptomic features of disease-associated microglia (DAM) in AD. Our results suggest that somatic driver mutations in microglia are common with normal aging but further enriched in AD brain, driving MiCE with inflammatory and DAM signatures. Our findings provide the first insights into microglial clonal dynamics in AD and identify potential new approaches to AD diagnosis and therapy.

Chromothripsis orchestrates leukemic transformation in blast phase MPN through targetable amplification of DYRK1A

Chromothripsis, the process of catastrophic shattering and haphazard repair of chromosomes, is a common event in cancer. Whether chromothripsis might constitute an actionable molecular event amenable to therapeutic targeting remains an open question. We describe recurrent chromothripsis of chromosome 21 in a subset of patients in blast phase of a myeloproliferative neoplasm (BP-MPN), which alongside other structural variants leads to amplification of a region of chromosome 21 in ∼25% of patients ('chr21amp'). We report that chr21amp BP-MPN has a particularly aggressive and treatment-resistant phenotype. The chr21amp event is highly clonal and present throughout the hematopoietic hierarchy. DYRK1A , a serine threonine kinase and transcription factor, is the only gene in the 2.7Mb minimally amplified region which showed both increased expression and chromatin accessibility compared to non-chr21amp BP-MPN controls. We demonstrate that DYRK1A is a central node at the nexus of multiple cellular functions critical for BP-MPN development, including DNA repair, STAT signalling and BCL2 overexpression. DYRK1A is essential for BP-MPN cell proliferation in vitro and in vivo , and DYRK1A inhibition synergises with BCL2 targeting to induce BP-MPN cell apoptosis. Collectively, these findings define the chr21amp event as a prognostic biomarker in BP-MPN and link chromothripsis to a druggable target.

STmut: a framework for visualizing somatic alterations in spatial transcriptomics data of cancer

Spatial transcriptomic technologies, such as the Visium platform, measure gene expression in different regions of tissues. Here, we describe new software, STmut, to visualize somatic point mutations, allelic imbalance, and copy number alterations in Visium data. STmut is tested on fresh-frozen Visium data, formalin-fixed paraffin-embedded (FFPE) Visium data, and tumors with and without matching DNA sequencing data. Copy number is inferred on all conditions, but the chemistry of the FFPE platform does not permit analyses of single nucleotide variants. Taken together, we propose solutions to add the genetic dimension to spatial transcriptomic data and describe the limitations of different datatypes.

A pan-cancer landscape of pathogenic somatic copy number variations

OBJECTIVE

Copy number variations (CNVs) play crucial roles in physiological and pathological processes, including cancer. However, the functional implications of somatic CNVs in tumor progression and evolution remain unclear. This study focuses on identifying CNV alterations with high pathogenic potential that drive and sustain tumorigenesis, distinguishing them from passenger alterations that accumulate during tumor growth. Our goal is to explore the variability of CNVs across different tumor types and infer their impact on tumor cell functions.

METHODS

Starting from 7352 copy number profiles across 33 different cancer types, we infer the pathogenicity of each CNV and perform both intra- and inter-tumor analyses to predict the functional impact of different genomic patterns. We evaluate the actionability of genes belonging to altered regions and we correlate the presence of pathogenic regions with genome instability patterns and patients' survival.

RESULTS

Our analysis uncovered large heterogeneity among different tumors suggesting in many cases distinct genetic drivers of tumorigenesis. Recurrent genomic alterations frequently coincide with dysfunctional homologous recombination pathways and negative regulation of the immune system. In certain tumors, the number of pathogenic CNVs emerged as a prognostic biomarker, highlighting their significance in cancer progression.

CONCLUSION

This study contributes to elucidate the functional impact of pathogenic CNVs in tumor progression and sheds light on their potential as prognostic markers in specific cancer types.

Acquired resistance to a GPRC5D-directed T-cell engager in multiple myeloma is mediated by genetic or epigenetic target inactivation

Bispecific antibodies targeting GPRC5D demonstrated promising efficacy in multiple myeloma, but acquired resistance usually occurs within a few months. Using a single-nucleus multi-omic strategy in three patients from the MYRACLE cohort (ClinicalTrials.gov registration: NCT03807128 ), we identified two resistance mechanisms, by bi-allelic genetic inactivation of GPRC5D or by long-range epigenetic silencing of its promoter and enhancer regions. Molecular profiling of target genes may help to guide the choice of immunotherapy and early detection of resistance in multiple myeloma.

A pan-tissue survey of mosaic chromosomal alterations in 948 individuals

Genetic mutations accumulate in an organism's body throughout its lifetime. While somatic single-nucleotide variants have been well characterized in the human body, the patterns and consequences of large chromosomal alterations in normal tissues remain largely unknown. Here, we present a pan-tissue survey of mosaic chromosomal alterations (mCAs) in 948 healthy individuals from the Genotype-Tissue Expression project, augmenting RNA-based allelic imbalance estimation with haplotype phasing. We found that approximately a quarter of the individuals carry a clonally-expanded mCA in at least one tissue, with incidence strongly correlated with age. The prevalence and genome-wide patterns of mCAs vary considerably across tissue types, suggesting tissue-specific mutagenic exposure and selection pressures. The mCA landscapes in normal adrenal and pituitary glands resemble those in tumors arising from these tissues, whereas the same is not true for the esophagus and skin. Together, our findings show a widespread age-dependent emergence of mCAs across normal human tissues with intricate connections to tumorigenesis.

Deciphering the heterogeneity and immunosuppressive function of regulatory T cells in osteosarcoma using single-cell RNA transcriptome

Osteosarcoma (OS) is a highly invasive malignant neoplasm with poor prognosis. The tumor microenvironment (TME) plays an essential role in the occurrence and development of OS. Regulatory T cells (Tregs) are known to facilitate immunosuppression, tumor progression, invasion, and metastasis. However, the effect of Tregs in the TME of OS remains unclear. In this study, single-cell RNA sequencing (scRNA-seq) data was used to identify Tregs and various other cell clusters in the TME of OS. Gene set variation analysis (GSVA) was used to investigate the signaling pathways in Tregs from OS and adjacent tissues. The CellChat and iTALK packages were used to analyze cellular communication. In addition, a prognostic model was established based on the Tregs-specific genes using bulk RNA-seq from the TARGET database, and it was verified using a Gene Expression Omnibus dataset. The pRRophetic package was used to predict drug sensitivity. Immunohistochemistry was used to verify the expression of candidate genes in OS. Based on the above methods, we showed that the OS samples were highly infiltrated with Tregs. GSVA revealed that oxidative phosphorylation, angiogenesis and mammalian target of rapamycin complex 1 (mTORC1) were highly activated in Tregs from OS compared with those from adjacent tissues. Using cellular communication analysis, we found that Tregs interacted with osteoblastic, endothelial, and myeloid cells via C-X-C motif chemokine ligand (CXCL) signaling; particularly, they strongly affected the expression of C-X-C motif chemokine receptor 4 (CXCR4) and interacted with other cell clusters through CXCL12/transforming growth factor β1 (TGFB1) to collectively enable tumor growth and progression. Subsequently, two Tregs-specific genes-CD320 and MAF-were screened through univariate, least absolute shrinkage and selection operator regression (LASSO) and multivariate analysis to construct a prognostic model, which showed excellent prognostic accuracy in two independent cohorts. In addition, drug sensitivity analysis demonstrated that OS patients at high Tregs risk were sensitive to sunitinib, sorafenib, and axitinib. We also used immunohistochemistry to validate that CD320 and MAF were significantly upregulated in OS tissues compared with adjacent tissues. Overall, this study reveals the heterogeneity of Tregs in the OS TME, providing new insights into the invasion and treatment of this cancer.

Neuronal DNA double-strand breaks lead to genome structural variations and 3D genome disruption in neurodegeneration

Persistent DNA double-strand breaks (DSBs) in neurons are an early pathological hallmark of neurodegenerative diseases including Alzheimer's disease (AD), with the potential to disrupt genome integrity. We used single-nucleus RNA-seq in human postmortem prefrontal cortex samples and found that excitatory neurons in AD were enriched for somatic mosaic gene fusions. Gene fusions were particularly enriched in excitatory neurons with DNA damage repair and senescence gene signatures. In addition, somatic genome structural variations and gene fusions were enriched in neurons burdened with DSBs in the CK-p25 mouse model of neurodegeneration. Neurons enriched for DSBs also had elevated levels of cohesin along with progressive multiscale disruption of the 3D genome organization aligned with transcriptional changes in synaptic, neuronal development, and histone genes. Overall, this study demonstrates the disruption of genome stability and the 3D genome organization by DSBs in neurons as pathological steps in the progression of neurodegenerative diseases.

Toward the functional interpretation of somatic structural variations: bulk- and single-cell approaches

Structural variants (SVs) are genomic rearrangements that can take many different forms such as copy number alterations, inversions and translocations. During cell development and aging, somatic SVs accumulate in the genome with potentially neutral, deleterious or pathological effects. Generation of somatic SVs is a key mutational process in cancer development and progression. Despite their importance, the detection of somatic SVs is challenging, making them less studied than somatic single-nucleotide variants. In this review, we summarize recent advances in whole-genome sequencing (WGS)-based approaches for detecting somatic SVs at the tissue and single-cell levels and discuss their advantages and limitations. First, we describe the state-of-the-art computational algorithms for somatic SV calling using bulk WGS data and compare the performance of somatic SV detectors in the presence or absence of a matched-normal control. We then discuss the unique features of cutting-edge single-cell-based techniques for analyzing somatic SVs. The advantages and disadvantages of bulk and single-cell approaches are highlighted, along with a discussion of their sensitivity to copy-neutral SVs, usefulness for functional inferences and experimental and computational costs. Finally, computational approaches for linking somatic SVs to their functional readouts, such as those obtained from single-cell transcriptome and epigenome analyses, are illustrated, with a discussion of the promise of these approaches in health and diseases.

Modeling specific aneuploidies: from karyotype manipulations to biological insights

An abnormal chromosome number, or aneuploidy, underlies developmental disorders and is a common feature of cancer, with different cancer types exhibiting distinct patterns of chromosomal gains and losses. To understand how specific aneuploidies emerge in certain tissues and how they contribute to disease development, various methods have been developed to alter the karyotype of mammalian cells and mice. In this review, we provide an overview of both classic and novel strategies for inducing or selecting specific chromosomal gains and losses in human and murine cell systems. We highlight how these customized aneuploidy models helped expanding our knowledge of the consequences of specific aneuploidies to (cancer) cell physiology.

Spatial transcriptomics reveals distinct and conserved tumor core and edge architectures that predict survival and targeted therapy response

The spatial organization of the tumor microenvironment has a profound impact on biology and therapy response. Here, we perform an integrative single-cell and spatial transcriptomic analysis on HPV-negative oral squamous cell carcinoma (OSCC) to comprehensively characterize malignant cells in tumor core (TC) and leading edge (LE) transcriptional architectures. We show that the TC and LE are characterized by unique transcriptional profiles, neighboring cellular compositions, and ligand-receptor interactions. We demonstrate that the gene expression profile associated with the LE is conserved across different cancers while the TC is tissue specific, highlighting common mechanisms underlying tumor progression and invasion. Additionally, we find our LE gene signature is associated with worse clinical outcomes while TC gene signature is associated with improved prognosis across multiple cancer types. Finally, using an in silico modeling approach, we describe spatially-regulated patterns of cell development in OSCC that are predictably associated with drug response. Our work provides pan-cancer insights into TC and LE biology and interactive spatial atlases ( http://www.pboselab.ca/spatial_OSCC/ ; http://www.pboselab.ca/dynamo_OSCC/ ) that can be foundational for developing novel targeted therapies.

Systematic characterization of therapeutic vulnerabilities in Multiple Myeloma with Amp1q reveals increased sensitivity to the combination of MCL1 and PI3K inhibitors

The development of targeted therapy for patients with Multiple Myeloma (MM) is hampered by the low frequency of actionable genetic abnormalities. Gain or amplification of chr1q (Amp1q) is the most frequent arm-level copy number gain in patients with MM, and it is associated with higher risk of progression and death despite recent advances in therapeutics. Thus, developing targeted therapy for patients with MM and Amp1q stands to benefit a large portion of patients in need of more effective management. Here, we employed large-scale dependency screens and drug screens to systematically characterize the therapeutic vulnerabilities of MM with Amp1q and showed increased sensitivity to the combination of MCL1 and PI3K inhibitors. Using single-cell RNA sequencing, we compared subclones with and without Amp1q within the same patient tumors and showed that Amp1q is associated with higher levels of MCL1 and the PI3K pathway. Furthermore, by isolating isogenic clones with different copy number for part of the chr1q arm, we showed increased sensitivity to MCL1 and PI3K inhibitors with arm-level gain. Lastly, we demonstrated synergy between MCL1 and PI3K inhibitors and dissected their mechanism of action in MM with Amp1q.

Single-cell profiling in multiple myeloma: insights, problems, and promises

In a short time, single-cell platforms have become the norm in many fields of research, including multiple myeloma (MM). In fact, the large amount of cellular heterogeneity in MM makes single-cell platforms particularly attractive because bulk assessments can miss valuable information about cellular subpopulations and cell-to-cell interactions. The decreasing cost and increasing accessibility of single-cell platform, combined with breakthroughs in obtaining multiomics data for the same cell and innovative computational programs for analyzing data, have allowed single-cell studies to make important insights into MM pathogenesis; yet, there is still much to be done. In this review, we will first focus on the types of single-cell profiling and the considerations for designing a single-cell profiling experiment. Then, we will discuss what have learned from single-cell profiling about myeloma clonal evolution, transcriptional reprogramming, and drug resistance, and about the MM microenvironment during precursor and advanced disease.

Clonally resolved single-cell multi-omics identifies routes of cellular differentiation in acute myeloid leukemia

Inter-patient variability and the similarity of healthy and leukemic stem cells (LSCs) have impeded the characterization of LSCs in acute myeloid leukemia (AML) and their differentiation landscape. Here, we introduce CloneTracer, a novel method that adds clonal resolution to single-cell RNA-seq datasets. Applied to samples from 19 AML patients, CloneTracer revealed routes of leukemic differentiation. Although residual healthy and preleukemic cells dominated the dormant stem cell compartment, active LSCs resembled their healthy counterpart and retained erythroid capacity. By contrast, downstream myeloid progenitors constituted a highly aberrant, disease-defining compartment: their gene expression and differentiation state affected both the chemotherapy response and leukemia's ability to differentiate into transcriptomically normal monocytes. Finally, we demonstrated the potential of CloneTracer to identify surface markers misregulated specifically in leukemic cells. Taken together, CloneTracer reveals a differentiation landscape that mimics its healthy counterpart and may determine biology and therapy response in AML.

Exploiting allele-specific transcriptional effects of subclonal copy number alterations for genotype-phenotype mapping in cancer cell populations

Somatic copy number alterations drive aberrant gene expression in cancer cells. In tumors with high levels of chromosomal instability, subclonal copy number alterations (CNAs) are a prevalent feature which often result in heterogeneous cancer cell populations with distinct phenotypes1. However, the extent to which subclonal CNAs contribute to clone-specific phenotypes remains poorly understood, in part due to the lack of methods to quantify how CNAs influence gene expression at a subclone level. We developed TreeAlign, which computationally integrates independently sampled single-cell DNA and RNA sequencing data from the same cell population and explicitly models gene dosage effects from subclonal alterations. We show through quantitative benchmarking data and application to human cancer data with single cell DNA and RNA libraries that TreeAlign accurately encodes clone-specific transcriptional effects of subclonal CNAs, the impact of allelic imbalance on allele-specific transcription, and obviates the need to arbitrarily define genotypic clones from a phylogenetic tree a priori. Combined, these advances lead to highly granular definitions of clones with distinct copy-number driven expression programs with increased resolution and accuracy over competing methods. The resulting improvement in assignment of transcriptional phenotypes to genomic clones enables clone-clone gene expression comparisons and explicit inference of genes that are mechanistically altered through CNAs, and identification of expression programs that are genomically independent. Our approach sets the stage for dissecting the relative contribution of fixed genomic alterations and dynamic epigenetic processes on gene expression programs in cancer.