A Bayesian mixture model for the analysis of allelic expression in single cells

Single-cell RNA Sequencing Reveals the Diversity of the Immunological Landscape Response to Genital Herpes

Genital herpes (GH) is a common sexually transmitted disease, which is primarily caused by herpes simplex virus type 2 (HSV-2), and continues to be a global health concern. Although our understanding of the alterations in immune cell populations and immunomodulation in GH patients is still limited, it is evident that systemic intrinsic immunity, innate immunity, and adaptive immunity play crucial roles during HSV-2 infection and GH reactivation. To investigate the mechanisms underlying HSV-2 infection and recurrence, single-cell RNA sequencing (scRNA-seq) was performed on immune cells isolated from the peripheral blood of both healthy individuals and patients with recurrent GH. Furthermore, the systemic immune response in patients with recurrent GH showed activation of classical monocytes, CD4+ T cells, natural killer cells (NK cells), and plasmacytoid dendritic cells (pDCs), especially of genes associated with the Toll-like receptor signaling pathway and T cell activation. Circulating immune cells in GH patients show higher expression of genes associated with inflammation and antiviral responses both in the scRNA-Seq data set and in independent quantitative real-time polymerase chain reaction (qRT-PCR) analysis and ELISA experiments. This study demonstrated that localized genital herpes, resulting from HSV reactivation, may influence the functionality of circulating immune cells, suggesting a potential avenue for future research into the role of systemic immunity during HSV infection and recurrence.

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.

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.

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.

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.

scAllele: A versatile tool for the detection and analysis of variants in scRNA-seq

Single-cell RNA sequencing (scRNA-seq) data contain rich information at the gene, transcript, and nucleotide levels. Most analyses of scRNA-seq have focused on gene expression profiles, and it remains challenging to extract nucleotide variants and isoform-specific information. Here, we present scAllele, an integrative approach that detects single-nucleotide variants, insertions, deletions, and their allelic linkage with splicing patterns in scRNA-seq. We demonstrate that scAllele achieves better performance in identifying nucleotide variants than other commonly used tools. In addition, the read-specific variant calls by scAllele enables allele-specific splicing analysis, a unique feature not afforded by other methods. Applied to a lung cancer scRNA-seq dataset, scAllele identified variants with strong allelic linkage to alternative splicing, some of which are cancer specific and enriched in cancer-relevant pathways. scAllele represents a versatile tool to uncover multilayer information and previously unidentified biological insights from scRNA-seq data.

A Novel Prognostic Signature Revealed the Interaction of Immune Cells in Tumor Microenvironment Based on Single-Cell RNA Sequencing for Lung Adenocarcinoma

Background

The tumor immune microenvironment (TIME) played an important role in immunotherapy prognosis and treatment response. Immune cells constitute a large part of the tumor microenvironment and regulate tumor progression. Our research is dedicated to studying the infiltrating immune cell in lung adenocarcinoma (LUAD) and seeking potential targets.

Methods

The scRNA-seq data were collected from our FDZSH and two public datasets. The code for cell-type mapping algorithms was downloaded from the CIBERSORTx portal. The bioinformatics data of LUAD patients could be approached from The Cancer Genome Atlas (TCGA) portal. Weighted gene coexpression network analysis (WGCNA) and least absolute shrinkage and selection operator (LASSO) analyses were performed to construct a risk model. TIMER2 and TIDE helped with the immune infiltration estimation, while PROGENy helped the cancer-related pathways' enrichment analysis. GSE31210 dataset and IMVigor ICB therapy cohort validated our findings as the external validation datasets.

Results

We clustered the scRNA-seq dataset (integrating our FDZSH datasets and other public datasets) into 23 subpopulations. After curated cell annotation, we implemented Cibersort and WGCNA analysis to anchor the brown module and natural killer cell cluster1 due to the most relationship with tumor trait. The overlap of the brown module gene, natural killer cell signature, and DEGs of tumor and adjacent normal samples was screened by LASSO Cox regression. The obtained 5-gene risk model showed an excellent prognostic performance in the validation dataset. Furthermore, there was a correlation between risk score and tumor-infiltrating immune cells and tumor genomics abnormity. Patients with higher risk scores had a significantly lower immunotherapy response rate.

Conclusion

Our observations implied that immune cells played a pivotal role in TIME and established a 5-gene signature (including IDH2, ADRB2, SFTPC, CCDC69, and CCND2) on the basement of nature killer markers targeted by WGCNA analysis. The significance of clinical outcome and immunotherapy response prediction was validated robustly.

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

MOTIVATION

Allelic expression analysis aids in detection of cis-regulatory mechanisms of genetic variation, which produce allelic imbalance (AI) in heterozygotes. Measuring AI in bulk data lacking time or spatial resolution has the limitation that cell-type-specific (CTS), spatial- or time-dependent AI signals may be dampened or not detected.

RESULTS

We introduce a statistical method airpart for identifying differential CTS AI from single-cell RNA-sequencing data, or dynamics AI from other spatially or time-resolved datasets. airpart outputs discrete partitions of data, pointing to groups of genes and cells under common mechanisms of cis-genetic regulation. In order to account for low counts in single-cell data, our method uses a Generalized Fused Lasso with Binomial likelihood for partitioning groups of cells by AI signal, and a hierarchical Bayesian model for AI statistical inference. In simulation, airpart accurately detected partitions of cell types by their AI and had lower Root Mean Square Error (RMSE) of allelic ratio estimates than existing methods. In real data, airpart identified differential allelic imbalance patterns across cell states and could be used to define trends of AI signal over spatial or time axes.

AVAILABILITY AND IMPLEMENTATION

The airpart package is available as an R/Bioconductor package at https://bioconductor.org/packages/airpart.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Power calculator for detecting allelic imbalance using hierarchical Bayesian model

OBJECTIVE

Allelic imbalance (AI) is the differential expression of the two alleles in a diploid. AI can vary between tissues, treatments, and environments. Methods for testing AI exist, but methods are needed to estimate type I error and power for detecting AI and difference of AI between conditions. As the costs of the technology plummet, what is more important: reads or replicates?

RESULTS

We find that a minimum of 2400, 480, and 240 allele specific reads divided equally among 12, 5, and 3 replicates is needed to detect a 10, 20, and 30%, respectively, deviation from allelic balance in a condition with power > 80%. A minimum of 960 and 240 allele specific reads divided equally among 8 replicates is needed to detect a 20 or 30% difference in AI between conditions with comparable power. Higher numbers of replicates increase power more than adding coverage without affecting type I error. We provide a Python package that enables simulation of AI scenarios and enables individuals to estimate type I error and power in detecting AI and differences in AI between conditions.

A Comprehensive Characterization of Monoallelic Expression During Hematopoiesis and Leukemogenesis via Single-Cell RNA-Sequencing

Single-cell RNA-sequencing (scRNA-seq) is becoming a powerful tool to investigate monoallelic expression (MAE) in various developmental and pathological processes. However, our knowledge of MAE during hematopoiesis and leukemogenesis is limited. In this study, we conducted a systematic interrogation of MAEs in bone marrow mononuclear cells (BMMCs) at single-cell resolution to construct a MAE atlas of BMMCs. We identified 1,020 constitutive MAEs in BMMCs, which included imprinted genes such as MEG8, NAP1L5, and IRAIN. We classified the BMMCs into six cell types and identified 74 cell type specific MAEs including MTSS1, MOB1A, and TCF12. We further identified 114 random MAEs (rMAEs) at single-cell level, with 78.1% single-allele rMAE and 21.9% biallelic mosaic rMAE. Many MAEs identified in BMMCs have not been reported and are potentially hematopoietic specific, supporting MAEs are functional relevance. Comparison of BMMC samples from a leukemia patient with multiple clinical stages showed the fractions of constitutive MAE were correlated with fractions of leukemia cells in BMMCs. Further separation of the BMMCs into leukemia cells and normal cells showed that leukemia cells have much higher constitutive MAE and rMAEs than normal cells. We identified the leukemia cell-specific MAEs and relapsed leukemia cell-specific MAEs, which were enriched in immune-related functions. These results indicate MAE is prevalent and is an important gene regulation mechanism during hematopoiesis and leukemogenesis. As the first systematical interrogation of constitutive MAEs, cell type specific MAEs, and rMAEs during hematopoiesis and leukemogenesis, the study significantly increased our knowledge about the features and functions of MAEs.

Perspectives on Allele-Specific Expression

Diploidy has profound implications for population genetics and susceptibility to genetic diseases. Although two copies are present for most genes in the human genome, they are not necessarily both active or active at the same level in a given individual. Genomic imprinting, resulting in exclusive or biased expression in favor of the allele of paternal or maternal origin, is now believed to affect hundreds of human genes. A far greater number of genes display unequal expression of gene copies due to cis-acting genetic variants that perturb gene expression. The availability of data generated by RNA sequencing applied to large numbers of individuals and tissue types has generated unprecedented opportunities to assess the contribution of genetic variation to allelic imbalance in gene expression. Here we review the insights gained through the analysis of these data about the extent of the genetic contribution to allelic expression imbalance, the tools and statistical models for gene expression imbalance, and what the results obtained reveal about the contribution of genetic variants that alter gene expression to complex human diseases and phenotypes.

Transcriptional and imprinting complexity in Arabidopsis seeds at single-nucleus resolution

Seeds are a key life cycle stage for many plants. Seeds are also the basis of agriculture and the primary source of calories consumed by humans¹. Here, we employ single-nucleus RNA-sequencing to generate a transcriptional atlas of developing Arabidopsis thaliana seeds, with a focus on endosperm. Endosperm, the primary site of gene imprinting in flowering plants, mediates the relationship between the maternal parent and the embryo². We identify transcriptionally uncharacterized nuclei types in the chalazal endosperm, which interfaces with maternal tissue for nutrient unloading^3,4. We demonstrate that the extent of parental bias of maternally expressed imprinted genes varies with cell-cycle phase, and that imprinting of paternally expressed imprinted genes is strongest in chalazal endosperm. Thus, imprinting is spatially and temporally heterogeneous. Increased paternal expression in the chalazal region suggests that parental conflict, which is proposed to drive imprinting evolution, is fiercest at the boundary between filial and maternal tissues.

A New Hypothesis for Type 1 Diabetes Risk: The At-Risk Allele at rs3842753 Associates With Increased Beta-cell INS Messenger RNA in a Meta-Analysis of Single-Cell RNA-Sequencing Data

OBJECTIVES

Type 1 diabetes is characterized by the autoimmune destruction of insulin-secreting beta cells. Genetic variants upstream at the insulin (INS) locus contribute to ∼10% of type 1 diabetes heritable risk. Previous studies showed an association between rs3842753 C/C genotype and type 1 diabetes susceptibility, but the molecular mechanisms remain unclear. To date, no large-scale studies have looked at the effect of genetic variation at rs3842753 on INS mRNA at the single-cell level.

METHODS

We aligned all human islet single-cell RNA sequencing data sets available to us in year 2020 to the reference genome GRCh38.98 and genotyped rs3842753, integrating 2,315 β cells and 1,223 β-like cells from 13 A/A protected donors, 23 A/C heterozygous donors and 35 C/C at-risk donors, including adults without diabetes and with type 2 diabetes.

RESULTS

INS expression mean and variance were significantly higher in single β cells from females compared with males. On comparing across β cells and β-like cells, we found that rs3842753 C‒containing cells (either homozygous or heterozygous) had the highest INS expression. We also found that β cells with the rs3842753 C allele had significantly higher endoplasmic reticulum stress marker gene expression compared with the A/A homozygous genotype.

CONCLUSIONS

These findings support the emerging concept that inherited risk of type 1 diabetes may be associated with inborn, persistent elevated insulin production, which may lead to β-cell endoplasmic reticulum stress and fragility.

Natural genetic variation determines microglia heterogeneity in wild-derived mouse models of Alzheimer's disease

Genetic and genome-wide association studies suggest a central role for microglia in Alzheimer’s disease (AD). However, single-cell RNA sequencing (scRNA-seq) of microglia in mice, a key preclinical model, has shown mixed results regarding translatability to human studies. To address this, scRNA-seq of microglia from C57BL/6J (B6) and wild-derived strains (WSB/EiJ, CAST/EiJ, and PWK/PhJ) with and without APP/PS1 demonstrates that genetic diversity significantly alters features and dynamics of microglia in baseline neuroimmune functions and in response to amyloidosis. Results show significant variation in the abundance of microglial subtypes or states, including numbers of previously identified disease-associated and interferon-responding microglia, across the strains. For each subtype, significant differences in the expression of many genes are observed in wild-derived strains relative to B6, including 19 genes previously associated with human AD including Apoe, Trem2, and Sorl1. This resource is critical in the development of appropriately targeted therapeutics for AD and other neurological diseases.

Neuroinflammation is a key component of Alzheimer’s disease. Yang et al. perform single-cell sequencing of microglia in wild-derived mouse strains that carry amyloid and show that these strains differ from the commonly used strains, exhibiting significant variation in abundance of microglial subtypes, including numbers of disease-associated and interferon-responding microglia.

Allele-specific expression: applications in cancer and technical considerations

Allele-specific gene expression can influence disease traits. Non-coding germline genetic variants that alter regulatory elements can cause allele-specific gene expression and contribute to cancer susceptibility. In tumors, both somatic copy number alterations and somatic single nucleotide variants have been shown to lead to allele-specific expression of genes, many of which are considered drivers of tumor growth. Here, we review recent studies revealing the pervasive presence of this phenomenon in cancer susceptibility and progression. Furthermore, we underscore the importance of careful experimental design and computational analysis for accurate allelic expression quantification and avoidance of false positives. Finally, we discuss additional methodological challenges encountered in cancer studies and in the burgeoning field of single-cell transcriptomics.