Leveraging single-cell sequencing to classify and characterize tumor subgroups in bulk RNA-sequencing data.

PURPOSE: Accurate classification of cancer subgroups is essential for precision medicine, tailoring treatments to individual patients based on their cancer subtypes. In recent years, advances in high-throughput sequencing technologies have enabled the generation of large-scale transcriptomic data from cancer samples. These data have provided opportunities for developing computational methods that can improve cancer subtyping and enable better personalized treatment strategies. METHODS: Here in this study, we evaluated different feature selection schemes in the context of meningioma classification. To integrate interpretable features from the bulk (n = 77 samples) and single-cell profiling (∼ 10 K cells), we developed an algorithm named CLIPPR which combines the top-performing single-cell models, RNA-inferred copy number variation (CNV) signals, and the initial bulk model to create a meta-model. RESULTS: While the scheme relying solely on bulk transcriptomic data showed good classification accuracy, it exhibited confusion between malignant and benign molecular classes in approximately ∼ 8% of meningioma samples. In contrast, models trained on features learned from meningioma single-cell data accurately resolved the sub-groups confused by bulk-transcriptomic data but showed limited overall accuracy. CLIPPR showed superior overall accuracy and resolved benign-malignant confusion as validated on n = 789 bulk meningioma samples gathered from multiple institutions. Finally, we showed the generalizability of our algorithm using our in-house single-cell (∼ 200 K cells) and bulk TCGA glioma data (n = 711 samples). CONCLUSION: Overall, our algorithm CLIPPR synergizes the resolution of single-cell data with the depth of bulk sequencing and enables improved cancer sub-group diagnoses and insights into their biology.

Integrated clinical genomic analysis reveals xenobiotic metabolic genes are downregulated in meningiomas of current smokers.

INTRODUCTION: Meningiomas are the most common primary intracranial tumor. Recently, various genetic classification systems for meningioma have been described. We sought to identify clinical drivers of different molecular changes in meningioma. As such, clinical and genomic consequences of smoking in patients with meningiomas remain unexplored. METHODS: 88 tumor samples were analyzed in this study. Whole exome sequencing (WES) was used to assess somatic mutation burden. RNA sequencing data was used to identify differentially expressed genes (DEG) and genes sets (GSEA). RESULTS: Fifty-seven patients had no history of smoking, twenty-two were past smokers, and nine were current smokers. The clinical data showed no major differences in natural history across smoking status. WES revealed absence of AKT1 mutation rate in current or past smokers compared to non-smokers (p = 0.046). Current smokers had increased mutation rate in NOTCH2 compared to past and never smokers (p < 0.05). Mutational signature from current and past smokers showed disrupted DNA mismatch repair (cosine-similarity = 0.759 and 0.783). DEG analysis revealed the xenobiotic metabolic genes UGT2A1 and UGT2A2 were both significantly downregulated in current smokers compared to past (Log2FC = - 3.97, padj = 0.0347 and Log2FC = - 4.18, padj = 0.0304) and never smokers (Log2FC = - 3.86, padj = 0.0235 and Log2FC = - 4.20, padj = 0.0149). GSEA analysis of current smokers showed downregulation of xenobiotic metabolism and enrichment for G2M checkpoint, E2F targets, and mitotic spindle compared to past and never smokers (FDR < 25% each). CONCLUSION: In this study, we conducted a comparative analysis of meningioma patients based on their smoking history, examining both their clinical trajectories and molecular changes. Meningiomas from current smokers were more likely to harbor NOTCH2 mutations, and AKT1 mutations were absent in current or past smokers. Moreover, both current and past smokers exhibited a mutational signature associated with DNA mismatch repair. Meningiomas from current smokers demonstrate downregulation of xenobiotic metabolic enzymes UGT2A1 and UGT2A2, which are downregulated in other smoking related cancers. Furthermore, current smokers exhibited downregulation xenobiotic metabolic gene sets, as well as enrichment in gene sets related to mitotic spindle, E2F targets, and G2M checkpoint, which are hallmark pathways involved in cell division and DNA replication control. In aggregate, our results demonstrate novel alterations in meningioma molecular biology in response to systemic carcinogens.

Identification of novel fusion transcripts in meningioma.

INTRODUCTION: Meningiomas are the most common primary intracranial tumor. Recent next generation sequencing analyses have elaborated the molecular drivers of this disease. We aimed to identify and characterize novel fusion genes in meningiomas. METHODS: We performed a secondary analysis of our RNA sequencing data of 145 primary meningioma from 140 patients to detect fusion genes. Semi-quantitative rt-PCR was performed to confirm transcription of the fusion genes in the original tumors. Whole exome sequencing was performed to identify copy number variations within each tumor sample. Comparative RNA seq analysis was performed to assess the clonality of the fusion constructs within the tumor. RESULTS: We detected six fusion events (NOTCH3-SETBP1, NF2-SPATA13, SLC6A3-AGBL3, PHF19-FOXP2 in two patients, and ITPK1-FBP2) in five out of 145 tumor samples. All but one event (NF2-SPATA13) led to extremely short reading frames, making these events de facto null alleles. Three of the five patients had a history of childhood radiation. Four out of six fusion events were detected in expression type C tumors, which represent the most aggressive meningioma. We validated the presence of the RNA transcripts in the tumor tissue by semi-quantitative RT PCR. All but the two PHF19-FOXP2 fusions demonstrated high degrees of clonality. CONCLUSIONS: Fusion genes occur infrequently in meningiomas and are more likely to be found in tumors with greater degree of genomic instability (expression type C) or in patients with history of cranial irradiation.

Comparative analysis of the transcriptome across distant species.

The transcriptome is the readout of the genome. Identifying common features in it across distant species can reveal fundamental principles. To this end, the ENCODE and modENCODE consortia have generated large amounts of matched RNA-sequencing data for human, worm and fly. Uniform processing and comprehensive annotation of these data allow comparison across metazoan phyla, extending beyond earlier within-phylum transcriptome comparisons and revealing ancient, conserved features. Specifically, we discover co-expression modules shared across animals, many of which are enriched in developmental genes. Moreover, we use expression patterns to align the stages in worm and fly development and find a novel pairing between worm embryo and fly pupae, in addition to the embryo-to-embryo and larvae-to-larvae pairings. Furthermore, we find that the extent of non-canonical, non-coding transcription is similar in each organism, per base pair. Finally, we find in all three organisms that the gene-expression levels, both coding and non-coding, can be quantitatively predicted from chromatin features at the promoter using a 'universal model' based on a single set of organism-independent parameters.

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.

Enhancing Fairness in Disease Prediction by Optimizing Multiple Domain Adversarial Networks.

Predictive models in biomedicine need to ensure equitable and reliable outcomes for the populations they are applied to. Unfortunately, biases in medical predictions can lead to unfair treatment and widening disparities, underscoring the need for effective techniques to address these issues. To enhance fairness, we introduce a framework based on a Multiple Domain Adversarial Neural Network (MDANN), which incorporates multiple adversarial components. In an MDANN, an adversarial module is applied to learn a fair pattern by negative gradients back-propagating across multiple sensitive features (i.e., characteristics of individuals that should not be used to discriminate unfairly between individuals when making predictions or decisions.) We leverage loss functions based on the Area Under the Receiver Operating Characteristic Curve (AUC) to address the class imbalance, promoting equitable classification performance for minority groups (e.g., a subset of the population that is underrepresented or disadvantaged.) Moreover, we utilize pre-trained convolutional autoencoders (CAEs) to extract deep representations of data, aiming to enhance prediction accuracy and fairness. Combining these mechanisms, we alleviate biases and disparities to provide reliable and equitable disease prediction. We empirically demonstrate that the MDANN approach leads to better accuracy and fairness in predicting disease progression using brain imaging data for Alzheimer's Disease and Autism populations than state-of-the-art techniques.

Single-Cell Profiling of CD8+ T Cells in Acute Myeloid Leukemia Reveals a Continuous Spectrum of Differentiation and Clonal Hyperexpansion.

Comprehensive investigation of CD8+ T cells in acute myeloid leukemia (AML) is essential for developing immunotherapeutic strategies beyond immune checkpoint blockade. Herein, we performed single-cell RNA profiling of CD8+ T cells from 3 healthy bone marrow donors and 23 newly diagnosed (NewlyDx) and 8 relapsed/refractory (RelRef) patients with AML. Cells coexpressing canonical exhaustion markers formed a cluster constituting <1% of all CD8+ T cells. We identified two effector CD8+ T-cell subsets characterized by distinct cytokine and metabolic profiles that were differentially enriched in NewlyDx and RelRef patients. We refined a 25-gene CD8-derived signature correlating with therapy resistance, including genes associated with activation, chemoresistance, and terminal differentiation. Pseudotemporal trajectory analysis supported enrichment of a terminally differentiated state in CD8+ T cells with high CD8-derived signature expression at relapse or refractory disease. Higher expression of the 25-gene CD8 AML signature correlated with poorer outcomes in previously untreated patients with AML, suggesting that the bona fide state of CD8+ T cells and their degree of differentiation are clinically relevant. Immune clonotype tracking revealed more phenotypic transitions in CD8 clonotypes in NewlyDx than in RelRef patients. Furthermore, CD8+ T cells from RelRef patients had a higher degree of clonal hyperexpansion associated with terminal differentiation and higher CD8-derived signature expression. Clonotype-derived antigen prediction revealed that most previously unreported clonotypes were patient-specific, suggesting significant heterogeneity in AML immunogenicity. Thus, immunologic reconstitution in AML is likely to be most successful at earlier disease stages when CD8+ T cells are less differentiated and have greater capacity for clonotype transitions.

Single-Cell Profiling of CD8+ T Cells in Acute Myeloid Leukemia Reveals a Continuous Spectrum of Differentiation and Clonal Hyperexpansion.

Comprehensive investigation of CD8+ T cells in acute myeloid leukemia (AML) is essential for developing immunotherapeutic strategies beyond immune checkpoint blockade. Herein, we performed single-cell RNA profiling of CD8+ T cells from 3 healthy bone marrow donors and 23 newly diagnosed (NewlyDx) and 8 relapsed/refractory (RelRef) AML patients. Cells co-expressing canonical exhaustion markers formed a cluster constituting <1% of all CD8+ T cells. We identified two effector CD8+ T cell subsets characterized by distinct cytokine and metabolic profiles that were differentially enriched in NewlyDx and RelRef patients. We refined a 25-gene CD8-derived signature correlating with therapy resistance, including genes associated with activation, chemoresistance, and terminal differentiation. Pseudotemporal trajectory analysis supported enrichment of a terminally differentiated state in CD8+ T cells with high CD8-derived signature expression at relapse or refractory disease. Higher expression of the 25-gene CD8 AML signature correlated with poorer outcomes in previously untreated AML patients, suggesting that the bona fide state of CD8+ T cells and their degree of differentiation are clinically relevant. Immune clonotype tracking revealed more phenotypic transitions in CD8 clonotypes in NewlyDx than in RelRef patients. Furthermore, CD8+ T cells from RelRef patients had a higher degree of clonal hyperexpansion associated with terminal differentiation and higher CD8-derived signature expression. Clonotype-derived antigen prediction revealed that most previously unreported clonotypes were patient-specific, suggesting significant heterogeneity in AML immunogenicity. Thus, immunologic reconstitution in AML is likely to be most successful at earlier disease stages when CD8+ T cells are less differentiated and have greater capacity for clonotype transitions.

Targeted gene expression profiling predicts meningioma outcomes and radiotherapy responses.

Surgery is the mainstay of treatment for meningioma, the most common primary intracranial tumor, but improvements in meningioma risk stratification are needed and indications for postoperative radiotherapy are controversial. Here we develop a targeted gene expression biomarker that predicts meningioma outcomes and radiotherapy responses. Using a discovery cohort of 173 meningiomas, we developed a 34-gene expression risk score and performed clinical and analytical validation of this biomarker on independent meningiomas from 12 institutions across 3 continents (N = 1,856), including 103 meningiomas from a prospective clinical trial. The gene expression biomarker improved discrimination of outcomes compared with all other systems tested (N = 9) in the clinical validation cohort for local recurrence (5-year area under the curve (AUC) 0.81) and overall survival (5-year AUC 0.80). The increase in AUC compared with the standard of care, World Health Organization 2021 grade, was 0.11 for local recurrence (95% confidence interval 0.07 to 0.17, P < 0.001). The gene expression biomarker identified meningiomas benefiting from postoperative radiotherapy (hazard ratio 0.54, 95% confidence interval 0.37 to 0.78, P = 0.0001) and suggested postoperative management could be refined for 29.8% of patients. In sum, our results identify a targeted gene expression biomarker that improves discrimination of meningioma outcomes, including prediction of postoperative radiotherapy responses.

Multiple approaches converge on three biological subtypes of meningioma and extract new insights from published studies.

One-fifth of meningiomas classified as benign by World Health Organization (WHO) histopathological grading will behave malignantly. To better diagnose these tumors, several groups turned to DNA methylation, whereas we combined RNA-sequencing (RNA-seq) and cytogenetics. Both approaches were more accurate than histopathology in identifying aggressive tumors, but whether they revealed similar tumor types was unclear. We therefore performed unbiased DNA methylation, RNA-seq, and cytogenetic profiling on 110 primary meningiomas WHO grade I and II). Each technique distinguished the same three groups (two benign and one malignant) as our previous molecular classification; integrating these methods into one classifier further improved accuracy. Computational modeling revealed strong correlations between transcription and cytogenetic changes, particularly loss of chromosome 1p, in malignant tumors. Applying our classifier to data from previous studies also resolved certain anomalies entailed by grouping tumors by WHO grade. Accurate classification will therefore elucidate meningioma biology as well as improve diagnosis and prognosis.

Racial and Socioeconomic Disparities in Patients With Meningioma: A Retrospective Cohort Study.

BACKGROUND: Meningiomas are the most common intracranial neoplasms. Although genomic analysis has helped elucidate differences in survival, there is evidence that racial disparities may influence outcomes. African Americans have a higher incidence of meningiomas and poorer survival outcomes. The etiology of these disparities remains unclear, but may include a combination of pathophysiology and other factors. OBJECTIVE: To determine factors that contribute to different clinical outcomes in racial populations. METHODS: We retrospectively reviewed 305 patients who underwent resection for meningiomas at a single tertiary care facility. We used descriptive statistics and univariate, multivariable, and Kaplan-Meier analyses to study clinical, radiographical, and histopathological differences. RESULTS: Minority patients were more likely to present through the emergency department than an outpatient clinic (P < .0001). They were more likely to present with more advanced clinical symptoms with lower Karnofsky Performance scores, more frequently had peritumoral edema (P = .0031), and experienced longer postoperative stays in the hospital (P = .0053), and African-American patients had higher hospitalization costs (P = .046) and were more likely to be publicly insured. Extent of resection was an independent predictor of recurrence freedom (P = .039). Presentation in clinic setting trended toward an association with recurrence-free survival (P = .055). We observed no significant difference in gross total resection rates, postoperative recurrence, or recurrence-free survival. CONCLUSION: Minority patients are more likely to present with severe symptoms, require longer perioperative hospitalization, and generate higher hospitalization costs. This may be due to socioeconomic factors that affect access to health care. Targeting barriers to access, especially to subspecialty care, may facilitate more appropriate and timely diagnosis, thereby improving patient care and outcomes.

TurboFold: iterative probabilistic estimation of secondary structures for multiple RNA sequences.

BACKGROUND: The prediction of secondary structure, i.e. the set of canonical base pairs between nucleotides, is a first step in developing an understanding of the function of an RNA sequence. The most accurate computational methods predict conserved structures for a set of homologous RNA sequences. These methods usually suffer from high computational complexity. In this paper, TurboFold, a novel and efficient method for secondary structure prediction for multiple RNA sequences, is presented. RESULTS: TurboFold takes, as input, a set of homologous RNA sequences and outputs estimates of the base pairing probabilities for each sequence. The base pairing probabilities for a sequence are estimated by combining intrinsic information, derived from the sequence itself via the nearest neighbor thermodynamic model, with extrinsic information, derived from the other sequences in the input set. For a given sequence, the extrinsic information is computed by using pairwise-sequence-alignment-based probabilities for co-incidence with each of the other sequences, along with estimated base pairing probabilities, from the previous iteration, for the other sequences. The extrinsic information is introduced as free energy modifications for base pairing in a partition function computation based on the nearest neighbor thermodynamic model. This process yields updated estimates of base pairing probability. The updated base pairing probabilities in turn are used to recompute extrinsic information, resulting in the overall iterative estimation procedure that defines TurboFold.TurboFold is benchmarked on a number of ncRNA datasets and compared against alternative secondary structure prediction methods. The iterative procedure in TurboFold is shown to improve estimates of base pairing probability with each iteration, though only small gains are obtained beyond three iterations. Secondary structures composed of base pairs with estimated probabilities higher than a significance threshold are shown to be more accurate for TurboFold than for alternative methods that estimate base pairing probabilities. TurboFold-MEA, which uses base pairing probabilities from TurboFold in a maximum expected accuracy algorithm for secondary structure prediction, has accuracy comparable to the best performing secondary structure prediction methods. The computational and memory requirements for TurboFold are modest and, in terms of sequence length and number of sequences, scale much more favorably than joint alignment and folding algorithms. CONCLUSIONS: TurboFold is an iterative probabilistic method for predicting secondary structures for multiple RNA sequences that efficiently and accurately combines the information from the comparative analysis between sequences with the thermodynamic folding model. Unlike most other multi-sequence structure prediction methods, TurboFold does not enforce strict commonality of structures and is therefore useful for predicting structures for homologous sequences that have diverged significantly. TurboFold can be downloaded as part of the RNAstructure package at http://rna.urmc.rochester.edu.

CaSpER identifies and visualizes CNV events by integrative analysis of single-cell or bulk RNA-sequencing data.

RNA sequencing experiments generate large amounts of information about expression levels of genes. Although they are mainly used for quantifying expression levels, they contain much more biologically important information such as copy number variants (CNVs). Here, we present CaSpER, a signal processing approach for identification, visualization, and integrative analysis of focal and large-scale CNV events in multiscale resolution using either bulk or single-cell RNA sequencing data. CaSpER integrates the multiscale smoothing of expression signal and allelic shift signals for CNV calling. The allelic shift signal measures the loss-of-heterozygosity (LOH) which is valuable for CNV identification. CaSpER employs an efficient methodology for the generation of a genome-wide B-allele frequency (BAF) signal profile from the reads and utilizes it for correction of CNVs calls. CaSpER increases the utility of RNA-sequencing datasets and complements other tools for complete characterization and visualization of the genomic and transcriptomic landscape of single cell and bulk RNA sequencing data.

HDAC inhibitors elicit metabolic reprogramming by targeting super-enhancers in glioblastoma models.

The Warburg effect is a tumor-related phenomenon that could potentially be targeted therapeutically. Here, we showed that glioblastoma (GBM) cultures and patients' tumors harbored super-enhancers in several genes related to the Warburg effect. By conducting a transcriptome analysis followed by ChIP-Seq coupled with a comprehensive metabolite analysis in GBM models, we found that FDA-approved global (panobinostat, vorinostat) and selective (romidepsin) histone deacetylase (HDAC) inhibitors elicited metabolic reprogramming in concert with disruption of several Warburg effect-related super-enhancers. Extracellular flux and carbon-tracing analyses revealed that HDAC inhibitors blunted glycolysis in a c-Myc-dependent manner and lowered ATP levels. This resulted in the engagement of oxidative phosphorylation (OXPHOS) driven by elevated fatty acid oxidation (FAO), rendering GBM cells dependent on these pathways. Mechanistically, interference with HDAC1/-2 elicited a suppression of c-Myc protein levels and a concomitant increase in 2 transcriptional drivers of oxidative metabolism, PGC1α and PPARD, suggesting an inverse relationship. Rescue and ChIP experiments indicated that c-Myc bound to the promoter regions of PGC1α and PPARD to counteract their upregulation driven by HDAC1/-2 inhibition. Finally, we demonstrated that combination treatment with HDAC and FAO inhibitors extended animal survival in patient-derived xenograft model systems in vivo more potently than single treatments in the absence of toxicity.

Patient-Derived Orthotopic Xenograft (PDOX) Mouse Models of Primary and Recurrent Meningioma.

BACKGROUND: Meningiomas constitute one-third of all primary brain tumors. Although typically benign, about 20% of these tumors recur despite surgery and radiation, and may ultimately prove fatal. There are currently no effective chemotherapies for meningioma. We, therefore, set out to develop patient-derived orthotopic xenograft (PDOX) mouse models of human meningioma using tumor. METHOD: Of nine patients, four had World Health Organization (WHO) grade I tumors, five had WHO grade II tumors, and in this second group two patients also had recurrent (WHO grade III) meningioma. We also classified the tumors according to our recently developed molecular classification system (Types A, B, and C, with C being the most aggressive). We transplanted all 11 surgical samples into the skull base of immunodeficient (SCID) mice. Only the primary and recurrent tumor cells from one patient-both molecular Type C, despite being WHO grades II and III, respectively-led to the formation of meningioma in the resulting mouse models. We characterized the xenografts by histopathology and RNA-seq and compared them with the original tumors. We performed an in vitro drug screen using 60 anti-cancer drugs followed by in vivo validation. RESULTS: The PDOX models established from the primary and recurrent tumors from patient K29 (K29P-PDOX and K29R-PDOX, respectively) replicated the histopathology and key gene expression profiles of the original samples. Although these xenografts could not be subtransplanted, the cryopreserved primary tumor cells were able to reliably generate PDOX tumors. Drug screening in K29P and K29R tumor cell lines revealed eight compounds that were active on both tumors, including three histone deacetylase (HDAC) inhibitors. We tested the HDAC inhibitor Panobinostat in K29R-PDOX mice, and it significantly prolonged mouse survival (p < 0.05) by inducing histone H3 acetylation and apoptosis. CONCLUSION: Meningiomas are not very amenable to PDOX modeling, for reasons that remain unclear. Yet at least some of the most malignant tumors can be modeled, and cryopreserved primary tumor cells can create large panels of tumors that can be used for preclinical drug testing.