Airpart: interpretable statistical models for analyzing allelic imbalance in single-cell datasets

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.

Detection of allele-specific expression in spatial transcriptomics with spASE

Spatial transcriptomics technologies permit the study of the spatial distribution of RNA at near-single-cell resolution genome-wide. However, the feasibility of studying spatial allele-specific expression (ASE) from these data remains uncharacterized. Here, we introduce spASE, a computational framework for detecting and estimating spatial ASE. To tackle the challenges presented by cell type mixtures and a low signal to noise ratio, we implement a hierarchical model involving additive mixtures of spatial smoothing splines. We apply our method to allele-resolved Visium and Slide-seq from the mouse cerebellum and hippocampus and report new insight into the landscape of spatial and cell type-specific ASE therein.

Single-cell genomics and regulatory networks for 388 human brains

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multiomics datasets into a resource comprising >2.8 million nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550,000 cell type-specific regulatory elements and >1.4 million single-cell expression quantitative trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Disentangling sex-dependent effects of APOE on diverse trajectories of cognitive decline in Alzheimer's disease

Current diagnostic systems for Alzheimer's disease (AD) rely upon clinical signs and symptoms, despite the fact that the multiplicity of clinical symptoms renders various neuropsychological assessments inadequate to reflect the underlying pathophysiological mechanisms. Since putative neuroimaging biomarkers play a crucial role in understanding the etiology of AD, we sought to stratify the diverse relationships between AD biomarkers and cognitive decline in the aging population and uncover risk factors contributing to the diversities in AD. To do so, we capitalized on a large amount of neuroimaging data from the ADNI study to examine the inflection points along the dynamic relationship between cognitive decline trajectories and whole-brain neuroimaging biomarkers, using a state-of-the-art statistical model of change point detection. Our findings indicated that the temporal relationship between AD biomarkers and cognitive decline may differ depending on the synergistic effect of genetic risk and biological sex. Specifically, tauopathy-PET biomarkers exhibit a more dynamic and age-dependent association with Mini-Mental State Examination scores (p<0.05), with inflection points at 72, 78, and 83 years old, compared with amyloid-PET and neurodegeneration (cortical thickness from MRI) biomarkers. In the landscape of health disparities in AD, our analysis indicated that biological sex moderates the rate of cognitive decline associated with APOE4 genotype. Meanwhile, we found that higher education levels may moderate the effect of APOE4, acting as a marker of cognitive reserve.

Single-cell genomics and regulatory networks for 388 human brains

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet, little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multi-omics datasets into a resource comprising >2.8M nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550K cell-type-specific regulatory elements and >1.4M single-cell expression-quantitative-trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Opportunities and tradeoffs in single-cell transcriptomic technologies

Recent technological and algorithmic advances enable single-cell transcriptomic analysis with remarkable depth and breadth. Nonetheless, a persistent challenge is the compromise between the ability to profile high numbers of cells and the achievement of full-length transcript coverage. Currently, the field is progressing and developing new and creative solutions that improve cellular throughput, gene detection sensitivity and full-length transcript capture. Furthermore, long-read sequencing approaches for single-cell transcripts are breaking frontiers that have previously blocked full transcriptome characterization. We here present a comprehensive overview of available options for single-cell transcriptome profiling, highlighting the key advantages and disadvantages of each approach.

Single-cell allele-specific expression analysis reveals dynamic and cell-type-specific regulatory effects

Differential allele-specific expression (ASE) is a powerful tool to study context-specific cis-regulation of gene expression. Such effects can reflect the interaction between genetic or epigenetic factors and a measured context or condition. Single-cell RNA sequencing (scRNA-seq) allows the measurement of ASE at individual-cell resolution, but there is a lack of statistical methods to analyze such data. We present Differential Allelic Expression using Single-Cell data (DAESC), a powerful method for differential ASE analysis using scRNA-seq from multiple individuals, with statistical behavior confirmed through simulation. DAESC accounts for non-independence between cells from the same individual and incorporates implicit haplotype phasing. Application to data from 105 induced pluripotent stem cell (iPSC) lines identifies 657 genes dynamically regulated during endoderm differentiation, with enrichment for changes in chromatin state. Application to a type-2 diabetes dataset identifies several differentially regulated genes between patients and controls in pancreatic endocrine cells. DAESC is a powerful method for single-cell ASE analysis and can uncover novel insights on gene regulation.

DNA methylation-associated allelic inactivation regulates Keratin 19 gene expression during pancreatic development and carcinogenesis

Within the pancreas, Keratin 19 (KRT19) labels the ductal lineage and is a determinant of pancreatic ductal adenocarcinoma (PDAC). To investigate KRT19 expression dynamics, we developed a human pluripotent stem cell (PSC)-based KRT19-mCherry reporter system in different genetic backgrounds to monitor KRT19 expression from its endogenous gene locus. A differentiation protocol to generate mature pancreatic duct-like organoids was applied. While KRT19/mCherry expression became evident at the early endoderm stage, mCherry signal was present in nearly all cells at the pancreatic endoderm (PE) and pancreatic progenitor (PP) stages. Interestingly, despite homogenous KRT19 expression, mCherry positivity dropped to 50% after ductal maturation, indicating a permanent switch from biallelic to monoallelic expression. DNA methylation profiling separated the distinct differentiation intermediates, with site-specific DNA methylation patterns occurring at the KRT19 locus during ductal maturation. Accordingly, the monoallelic switch was partially reverted upon treatment with a DNA-methyltransferase inhibitor. In human PDAC cohorts, high KRT19 levels correlate with low locus methylation and decreased survival. At the same time, activation of oncogenic KRASG12D signalling in our reporter system reversed monoallelic back to biallelic KRT19 expression in pancreatic duct-like organoids. Allelic reactivation was also detected in single-cell transcriptomes of human PDACs, which further revealed a positive correlation between KRT19 and KRAS expression. Accordingly, KRAS mutant PDACs had higher KRT19 mRNA but lower KRT19 gene locus DNA methylation than wildtype counterparts. KRT19 protein was additionally detected in plasma of PDAC patients, with higher concentrations correlating with shorter progression-free survival in gemcitabine/nabPaclitaxel-treated and opposing trends in FOLFIRINOX-treated patients. Apart from being an important pancreatic ductal lineage marker, KRT19 appears tightly controlled via a switch from biallelic to monoallelic expression during ductal lineage entry and is aberrantly expressed after oncogenic KRASG12D expression, indicating a role in PDAC development and malignancy. Soluble KRT19 might serve as a relevant biomarker to stratify treatment. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Modelling capture efficiency of single-cell RNA-sequencing data improves inference of transcriptome-wide burst kinetics

MOTIVATION

Gene expression is characterized by stochastic bursts of transcription that occur at brief and random periods of promoter activity. The kinetics of gene expression burstiness differs across the genome and is dependent on the promoter sequence, among other factors. Single-cell RNA sequencing (scRNA-seq) has made it possible to quantify the cell-to-cell variability in transcription at a global genome-wide level. However, scRNA-seq data are prone to technical variability, including low and variable capture efficiency of transcripts from individual cells.

RESULTS

Here, we propose a novel mathematical theory for the observed variability in scRNA-seq data. Our method captures burst kinetics and variability in both the cell size and capture efficiency, which allows us to propose several likelihood-based and simulation-based methods for the inference of burst kinetics from scRNA-seq data. Using both synthetic and real data, we show that the simulation-based methods provide an accurate, robust and flexible tool for inferring burst kinetics from scRNA-seq data. In particular, in a supervised manner, a simulation-based inference method based on neural networks proves to be accurate and useful when applied to both allele and nonallele-specific scRNA-seq data.

AVAILABILITY AND IMPLEMENTATION

The code for Neural Network and Approximate Bayesian Computation inference is available at https://github.com/WT215/nnRNA and https://github.com/WT215/Julia_ABC, respectively.