SCReadCounts: estimation of cell-level SNVs expression from scRNA-seq data

Biomedical Utility of Non-Enzymatic DNA Amplification Reaction: From Material Design to Diagnosis and Treatment

Nucleic acid nanotechnology has become a promising strategy for disease diagnosis and treatment, owing to remarkable programmability, precision, and biocompatibility. However, current biosensing and biotherapy approaches by nucleic acids exhibit limitations in sensitivity, specificity, versatility, and real-time monitoring. DNA amplification reactions present an advantageous strategy to enhance the performance of biosensing and biotherapy platforms. Non-enzymatic DNA amplification reaction (NEDAR), such as hybridization chain reaction and catalytic hairpin assembly, operate via strand displacement. NEDAR presents distinct advantages over traditional enzymatic DNA amplification reactions, including simplified procedures, milder reaction conditions, higher specificity, enhanced controllability, and excellent versatility. Consequently, research focusing on NEDAR-based biosensing and biotherapy has garnered significant attention. NEDAR demonstrates high efficacy in detecting multiple types of biomarkers, including nucleic acids, small molecules, and proteins, with high sensitivity and specificity, enabling the parallel detection of multiple targets. Besides, NEDAR can strengthen drug therapy, cellular behavior control, and cell encapsulation. Moreover, NEDAR holds promise for constructing assembled diagnosis-treatment nanoplatforms in the forms of pure DNA nanostructures and hybrid nanomaterials, which offer utility in disease monitoring and precise treatment. Thus, this paper aims to comprehensively elucidate the reaction mechanism of NEDAR and review the substantial advancements in NEDAR-based diagnosis and treatment over the past five years, encompassing NEDAR-based design strategies, applications, and prospects.

Computational methods for allele-specific expression in single cells

Allele-specific expression (ASE) is a powerful signal that can be used to investigate multiple molecular mechanisms, such as cis-regulatory effects and imprinting. Single-cell RNA-sequencing (scRNA-seq) enables ASE characterization at the resolution of individual cells. In this review, we highlight the computational methods for processing and analyzing single-cell ASE data. We first describe a bioinformatics pipeline to obtain ASE counts from raw reads synthesized from previous literature. We then discuss statistical methods for detecting allelic imbalance and its variability across conditions using scRNA-seq data. In addition, we describe other methods that use single-cell ASE to address specific biological questions. Finally, we discuss future directions and emphasize the need for an integrated, optimized bioinformatics pipeline, and further development of statistical methods for different technologies.

CanCellVar: A database for single-cell variants map in human cancer

Numerous variants, including both single-nucleotide variants (SNVs) in DNA and A>G RNA edits in mRNA as essential drivers of cellular proliferation and tumorigenesis, are commonly associated with cancer progression and growth. Thus, mining and summarizing single-cell variants will provide a refined and higher-resolution view of cancer and further contribute to precision medicine. Here, we established a database, CanCellVar, which aims to provide and visualize the comprehensive atlas of single-cell variants in tumor microenvironment. The current CanCellVar identified ∼3 million variants (∼1.4 million SNVs and ∼1.4 million A>G RNA edits) involved in 2,754,531 cells of 5 major cell types across 37 cancer types. CanCellVar provides the basic annotation information as well as cellular and molecular function properties of variants. In addition, the clinical relevance of variants can be obtained including tumor grade, treatment, metastasis, and others. Several flexible tools were also developed to aid retrieval and to analyze cell-cell interactions, gene expression, cell-development trajectories, regulation, and molecular structure affected by variants. Collectively, CanCellVar will serve as a valuable resource for investigating the functions and characteristics of single-cell variations and their roles in human tumor evolution and treatment.

ScSNViz: a user-friendly toolset for visualization and analysis of Cell-Specific Expressed SNVs

Motivation

Understanding genetic variation at the single-cell level is crucial for insights into cellular heterogeneity, clonal evolution, and gene expression regulation, but there is a scarcity of tools for visualizing and analyzing cell-level genetic variants.

Results

We introduce scSNViz, a comprehensive R-based toolset for visualization and analysis of cell-specific expressed Single Nucleotide Variants (sceSNVs) within cell-barcoded single-cell RNA-sequencing (scRNA-seq) data. ScSNViz offers 3D sceSNV visualization capabilities for dimensionally reduced scRNA-seq gene expression data, compatibility with popular scRNA-seq processing tools like Seurat, cell-type classification tools such as SingleR and scType, and trajectory inference computation using Slingshot. Furthermore, scSNViz conducts estimation, summary, and graphical representation of statistical metrics pertaining to sceSNVs distribution and expression across individual cells. It also provides support for the analysis of individual sceSNVs as well as sets comprising multiple expressed sceSNVs of interest.

Availability

ScSNViz is implemented as user-friendly R-scripts, freely available on https://horvathlab.github.io/NGS/scSNViz , supported by help utilities, and requiring no specialized bioinformatics skills for use.

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.

Single-cell multiomics reveals the interplay of clonal evolution and cellular plasticity in hepatoblastoma

Hepatoblastomas (HB) display heterogeneous cellular phenotypes that influence the clinical outcome, but the underlying mechanisms are poorly understood. Here, we use a single-cell multiomic strategy to unravel the molecular determinants of this plasticity. We identify a continuum of HB cell states between hepatocytic (scH), liver progenitor (scLP) and mesenchymal (scM) differentiation poles, with an intermediate scH/LP population bordering scLP and scH areas in spatial transcriptomics. Chromatin accessibility landscapes reveal the gene regulatory networks of each differentiation pole, and the sequence of transcription factor activations underlying cell state transitions. Single-cell mapping of somatic alterations reveals the clonal architecture of each tumor, showing that each genetic subclone displays its own range of cellular plasticity across differentiation states. The most scLP subclones, overexpressing stem cell and DNA repair genes, proliferate faster after neo-adjuvant chemotherapy. These results highlight how the interplay of clonal evolution and epigenetic plasticity shapes the potential of HB subclones to respond to chemotherapy.

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.

Preneoplastic liver colonization by 11p15.5 altered mosaic cells in young children with hepatoblastoma

Pediatric liver tumors are very rare tumors with the most common diagnosis being hepatoblastoma. While hepatoblastomas are predominantly sporadic, around 15% of cases develop as part of predisposition syndromes such as Beckwith-Wiedemann (11p15.5 locus altered). Here, we identify mosaic genetic alterations of 11p15.5 locus in the liver of hepatoblastoma patients without a clinical diagnosis of Beckwith-Wiedemann syndrome. We do not retrieve these alterations in children with other types of pediatric liver tumors. We show that mosaic 11p15.5 alterations in liver FFPE sections of hepatoblastoma patients display IGF2 overexpression and H19 downregulation together with an alteration of the liver zonation. Moreover, mosaic livers' microenvironment is enriched in extracellular matrix and angiogenesis. Spatial transcriptomics and single-nucleus RNAseq analyses identify a 60-gene signature in 11p15.5 altered hepatocytes. These data provide insights for 11p15.5 mosaicism detection and its functional consequences during the early steps of carcinogenesis.

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 targeted sequencing extension for transcript genotyping in single-cell transcriptomics

As no existing methods within the single-cell RNA sequencing repertoire combine genotyping of specific genomic loci with high throughput, we evaluated a straightforward, targeted sequencing approach as an extension to high-throughput droplet-based single-cell RNA sequencing. Overlaying standard gene expression data with transcript level genotype information provides a strategy to study the impact of genetic variants. Here, we describe this targeted sequencing extension, explain how to process the data and evaluate how technical parameters such as amount of input cDNA, number of amplification rounds, and sequencing depth influence the number of transcripts detected. Finally, we demonstrate how targeted sequencing can be used in two contexts: (1) simultaneous investigation of the presence of a somatic variant and its potential impact on the transcriptome of affected cells and (2) evaluation of allele-specific expression of a germline variant in ad hoc cell subsets. Through these and other comparable applications, our targeted sequencing extension has the potential to improve our understanding of functional effects caused by genetic variation.

Single-nuclei transcriptomics enable detection of somatic variants in patient brain tissue

Somatic variants are a major cause of human disease, including neurological disorders like focal epilepsies, but can be challenging to study due to their mosaicism in bulk tissue biopsies. Coupling single-cell genotype and transcriptomic data has potential to provide insight into the role somatic variants play in disease etiology, such as by determining what cell types are affected or how the mutations affect gene expression. Here, we asked whether commonly used single-nucleus 3'- or 5'-RNA-sequencing assays can be used to derive single-nucleus genotype data for a priori known variants that are located near to either end of a transcript. To that end, we compared performance of commercially available single-nuclei 3'- and 5'- gene expression kits using resected brain samples from three pediatric patients with focal epilepsy. We quantified the ability to detect genetic variants in single-nucleus datasets depending on distance from the transcript end. Finally, we demonstrated the ability to identify affected cell types in a patient with a RHEB somatic variant causing an epilepsy-associated cortical malformation. Our results demonstrate that single-nuclei 3' or 5'-RNA-sequencing data can be used to identify known somatic variants in single-nuclei when they are expressed within proximity to a transcript end.

SCExecute: custom cell barcode-stratified analyses of scRNA-seq data

MOTIVATION

In single-cell RNA-sequencing (scRNA-seq) data, stratification of sequencing reads by cellular barcode is necessary to study cell-specific features. However, apart from gene expression, the analyses of cell-specific features are not sufficiently supported by available tools designed for high-throughput sequencing data.

RESULTS

We introduce SCExecute, which executes a user-provided command on barcode-stratified, extracted on-the-fly, single-cell binary alignment map (scBAM) files. SCExecute extracts the alignments with each cell barcode from aligned, pooled single-cell sequencing data. Simple commands, monolithic programs, multi-command shell scripts or complex shell-based pipelines are then executed on each scBAM file. scBAM files can be restricted to specific barcodes and/or genomic regions of interest. We demonstrate SCExecute with two popular variant callers-GATK and Strelka2-executed in shell-scripts together with commands for BAM file manipulation and variant filtering, to detect single-cell-specific expressed single nucleotide variants from droplet scRNA-seq data (10X Genomics Chromium System).In conclusion, SCExecute facilitates custom cell-level analyses on barcoded scRNA-seq data using currently available tools and provides an effective solution for studying low (cellular) frequency transcriptome features.

AVAILABILITY AND IMPLEMENTATION

SCExecute is implemented in Python3 using the Pysam package and distributed for Linux, MacOS and Python environments from https://horvathlab.github.io/NGS/SCExecute.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

RDAClone: Deciphering Tumor Heterozygosity through Single-Cell Genomics Data Analysis with Robust Deep Autoencoder

Rapid advances in single-cell genomics sequencing (SCGS) have allowed researchers to characterize tumor heterozygosity with unprecedented resolution and reveal the phylogenetic relationships between tumor cells or clones. However, high sequencing error rates of current SCGS data, i.e., false positives, false negatives, and missing bases, severely limit its application. Here, we present a deep learning framework, RDAClone, to recover genotype matrices from noisy data with an extended robust deep autoencoder, cluster cells into subclones by the Louvain-Jaccard method, and further infer evolutionary relationships between subclones by the minimum spanning tree. Studies on both simulated and real datasets demonstrate its robustness and superiority in data denoising, cell clustering, and evolutionary tree reconstruction, particularly for large datasets.

Improved SNV Discovery in Barcode-Stratified scRNA-seq Alignments

Currently, the detection of single nucleotide variants (SNVs) from 10 x Genomics single-cell RNA sequencing data (scRNA-seq) is typically performed on the pooled sequencing reads across all cells in a sample. Here, we assess the gaining of information regarding SNV assessments from individual cell scRNA-seq data, wherein the alignments are split by cellular barcode prior to the variant call. We also reanalyze publicly available data on the MCF7 cell line during anticancer treatment. We assessed SNV calls by three variant callers—GATK, Strelka2, and Mutect2, in combination with a method for the cell-level tabulation of the sequencing read counts bearing variant alleles–SCReadCounts (single-cell read counts). Our analysis shows that variant calls on individual cell alignments identify at least a two-fold higher number of SNVs as compared to the pooled scRNA-seq; these SNVs are enriched in novel variants and in stop-codon and missense substitutions. Our study indicates an immense potential of SNV calls from individual cell scRNA-seq data and emphasizes the need for cell-level variant detection approaches and tools, which can contribute to the understanding of the cellular heterogeneity and the relationships to phenotypes, and help elucidate somatic mutation evolution and functionality.