Detection of isoforms and genomic alterations by high-throughput full-length single-cell RNA sequencing in ovarian cancer

Understanding the complex background of cancer requires genotype-phenotype information in single-cell resolution. Here, we perform long-read single-cell RNA sequencing (scRNA-seq) on clinical samples from three ovarian cancer patients presenting with omental metastasis and increase the PacBio sequencing depth to 12,000 reads per cell. Our approach captures 152,000 isoforms, of which over 52,000 were not previously reported. Isoform-level analysis accounting for non-coding isoforms reveals 20% overestimation of protein-coding gene expression on average. We also detect cell type-specific isoform and poly-adenylation site usage in tumor and mesothelial cells, and find that mesothelial cells transition into cancer-associated fibroblasts in the metastasis, partly through the TGF-β/miR-29/Collagen axis. Furthermore, we identify gene fusions, including an experimentally validated IGF2BP2::TESPA1 fusion, which is misclassified as high TESPA1 expression in matched short-read data, and call mutations confirmed by targeted NGS cancer gene panel results. With these findings, we envision long-read scRNA-seq to become increasingly relevant in oncology and personalized medicine.

Transcriptional linkage analysis with in vivo AAV-Perturb-seq

The ever-growing compendium of genetic variants associated with human pathologies demands new methods to study genotype-phenotype relationships in complex tissues in a high-throughput manner1,2. Here we introduce adeno-associated virus (AAV)-mediated direct in vivo single-cell CRISPR screening, termed AAV-Perturb-seq, a tuneable and broadly applicable method for transcriptional linkage analysis as well as high-throughput and high-resolution phenotyping of genetic perturbations in vivo. We applied AAV-Perturb-seq using gene editing and transcriptional inhibition to systematically dissect the phenotypic landscape underlying 22q11.2 deletion syndrome3,4 genes in the adult mouse brain prefrontal cortex. We identified three 22q11.2-linked genes involved in known and previously undescribed pathways orchestrating neuronal functions in vivo that explain approximately 40% of the transcriptional changes observed in a 22q11.2-deletion mouse model. Our findings suggest that the 22q11.2-deletion syndrome transcriptional phenotype found in mature neurons may in part be due to the broad dysregulation of a class of genes associated with disease susceptibility that are important for dysfunctional RNA processing and synaptic function. Our study establishes a flexible and scalable direct in vivo method to facilitate causal understanding of biological and disease mechanisms with potential applications to identify genetic interventions and therapeutic targets for treating disease.

Dynamic thresholding and tissue dissociation optimization for CITE-seq identifies differential surface protein abundance in metastatic melanoma

Multi-omics profiling by CITE-seq bridges the RNA-protein gap in single-cell analysis but has been largely applied to liquid biopsies. Applying CITE-seq to clinically relevant solid biopsies to characterize healthy tissue and the tumor microenvironment is an essential next step in single-cell translational studies. In this study, gating of cell populations based on their transcriptome signatures for use in cell type-specific ridge plots allowed identification of positive antibody signals and setting of manual thresholds. Next, we compare five skin dissociation protocols by taking into account dissociation efficiency, captured cell type heterogeneity and recovered surface proteome. To assess the effect of enzymatic digestion on transcriptome and epitope expression in immune cell populations, we analyze peripheral blood mononuclear cells (PBMCs) with and without dissociation. To further assess the RNA-protein gap, RNA-protein we perform codetection and correlation analyses on thresholded protein values. Finally, in a proof-of-concept study, using protein abundance analysis on selected surface markers in a cohort of healthy skin, primary, and metastatic melanoma we identify CD56 surface marker expression on metastatic melanoma cells, which was further confirmed by multiplex immunohistochemistry. This work provides practical guidelines for processing and analysis of clinically relevant solid tissue biopsies for biomarker discovery.

gExcite - A start-to-end framework for single-cell gene expression, hashing, and antibody analysis

SUMMARY

Recently, CITE-seq emerged as a multimodal single-cell technology capturing gene expression and surface protein information from the same single-cells, which allows unprecedented insights into disease mechanisms and heterogeneity, as well as immune cell profiling. Multiple single-cell profiling methods exist, but they are typically focussed on either gene expression or antibody analysis, not their combination. Moreover, existing software suites are not easily scalable to a multitude of samples. To this end, we designed gExcite, a start-to-end workflow that provides both gene and antibody expression analysis, as well as hashing deconvolution. Embedded in the Snakemake workflow manager, gExcite facilitates reproducible and scalable analyses. We showcase the output of gExcite on a study of different dissociation protocols on PBMC samples.

AVAILABILITY

gExcite is open source available on github at https://github.com/ETH-NEXUS/gExcite_pipeline. The software is distributed under the GNU General Public License 3 (GPL3).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

ROS Induction Targets Persister Cancer Cells with Low Metabolic Activity in NRAS-Mutated Melanoma

Clinical management of melanomas with NRAS mutations is challenging. Targeting MAPK signaling is only beneficial to a small subset of patients due to resistance that arises through genetic, transcriptional, and metabolic adaptation. Identification of targetable vulnerabilities in NRAS-mutated melanoma could help improve patient treatment. Here, we used multiomics analyses to reveal that NRAS-mutated melanoma cells adopt a mesenchymal phenotype with a quiescent metabolic program to resist cellular stress induced by MEK inhibition. The metabolic alterations elevated baseline reactive oxygen species (ROS) levels, leading these cells to become highly sensitive to ROS induction. In vivo xenograft experiments and single-cell RNA sequencing demonstrated that intratumor heterogeneity necessitates the combination of a ROS inducer and a MEK inhibitor to inhibit both tumor growth and metastasis. Ex vivo pharmacoscopy of 62 human metastatic melanomas confirmed that MEK inhibitor-resistant tumors significantly benefited from the combination therapy. Finally, oxidative stress response and translational suppression corresponded with ROS-inducer sensitivity in 486 cancer cell lines, independent of cancer type. These findings link transcriptional plasticity to a metabolic phenotype that can be inhibited by ROS inducers in melanoma and other cancers.

SIGNIFICANCE

Metabolic reprogramming in drug-resistant NRAS-mutated melanoma cells confers sensitivity to ROS induction, which suppresses tumor growth and metastasis in combination with MAPK pathway inhibitors.