Gene expression imputation and cell-type deconvolution in human brain with spatiotemporal precision and its implications for brain-related disorders

Differentially Altered Metabolic Pathways in the Amygdala of Subjects with Schizophrenia, Bipolar Disorder and Major Depressive Disorder

Background and hypothesis

A growing number of studies implicate a key role for metabolic processes in psychiatric disorders. Recent studies suggest that ketogenic diet may be therapeutically effective for subgroups of people with schizophrenia (SCZ), bipolar disorder (BPD) and possibly major depressive disorder (MDD). Despite this promise, there is currently limited information regarding brain energy metabolism pathways across these disorders, limiting our understanding of how brain metabolic pathways are altered and who may benefit from ketogenic diets. We conducted gene expression profiling on the amygdala, a key region involved in in the regulation of mood and appetitive behaviors, to test the hypothesis that amygdala metabolic pathways are differentially altered between these disorders.

Study Design

We used a cohort of subjects diagnosed with SCZ, BPD or MDD, and non-psychiatrically ill control subjects (n=15/group), together with our bioinformatic 3-pod analysis consisting of full transcriptome pathway analysis, targeted pathway analysis, leading-edge gene analysis and iLINCS perturbagen analysis.

Study Results

We identified differential expression of metabolic pathways in each disorder. Subjects with SCZ displayed downregulation of mitochondrial respiration and nucleotide metabolism pathways. In comparison, we observed upregulation of mitochondrial respiration pathways in subjects with MDD, while subjects with BPD displayed enrichment of pathways involved in carbohydrate metabolism. Several pathways associated with brain metabolism including immune system processes and calcium ion transport were also differentially altered between diagnosis groups.

Conclusion

Our findings suggest metabolic pathways are differentially altered in the amygdala in these disorders, which may impact approaches for therapeutic strategies.

lute: estimating the cell composition of heterogeneous tissue with varying cell sizes using gene expression

Relative cell type fraction estimates in bulk RNA-sequencing data are important to control for cell composition differences across heterogenous tissue samples. Current computational tools estimate relative RNA abundances rather than cell type proportions in tissues with varying cell sizes, leading to biased estimates. We present lute, a computational tool to accurately deconvolute cell types with varying sizes. Our software wraps existing deconvolution algorithms in a standardized framework. Using simulated and real datasets, we demonstrate how lute adjusts for differences in cell sizes to improve the accuracy of cell composition. Software is available from https://bioconductor.org/packages/lute.

Robust differences in cortical cell type proportions across healthy human aging inferred through cross-dataset transcriptome analyses

Age-related declines in cognitive function are driven by cell type-specific changes in the brain. However, it remains challenging to study cellular differences associated with healthy aging as traditional approaches scale poorly to the sample sizes needed to capture aging and cellular heterogeneity. Here, we employed cellular deconvolution to estimate relative cell type proportions using frontal cortex bulk gene expression from individuals without psychiatric conditions or brain pathologies. Our analyses comprised 8 datasets and 6 cohorts (1142 subjects and 1429 samples) with ages of death spanning 15-90 years. We found aging associated with profound differences in cellular proportions, with the largest changes reflecting fewer somatostatin- and vasoactive intestinal peptide-expressing interneurons, more astrocytes and other non-neuronal cells, and a suggestive "U-shaped" quadratic relationship for microglia. Cell type associations with age were markedly robust across bulk-and single nucleus datasets. Altogether, we present a comprehensive account of proportional differences in cortical cell types associated with healthy aging.

deCS: A Tool for Systematic Cell Type Annotations of Single-cell RNA Sequencing Data among Human Tissues

Single-cell RNA sequencing (scRNA-seq) is revolutionizing the study of complex and dynamic cellular mechanisms. However, cell type annotation remains a main challenge as it largely relies on a priori knowledge and manual curation, which is cumbersome and subjective. The increasing number of scRNA-seq datasets, as well as numerous published genetic studies, has motivated us to build a comprehensive human cell type reference atlas.Here, we present decoding Cell type Specificity (deCS), an automatic cell type annotation method augmented by a comprehensive collection of human cell type expression profiles and marker genes. We used deCS to annotate scRNA-seq data from various tissue types and systematically evaluated the annotation accuracy under different conditions, including reference panels, sequencing depth, and feature selection strategies. Our results demonstrate that expanding the references is critical for improving annotation accuracy. Compared to many existing state-of-the-art annotation tools, deCS significantly reduced computation time and increased accuracy. deCS can be integrated into the standard scRNA-seq analytical pipeline to enhance cell type annotation. Finally, we demonstrated the broad utility of deCS to identify trait-cell type associations in 51 human complex traits, providing deep insights into the cellular mechanisms underlying disease pathogenesis. All documents for deCS, including source code, user manual, demo data, and tutorials, are freely available at https://github.com/bsml320/deCS.

Evaluating performance and applications of sample-wise cell deconvolution methods on human brain transcriptomic data

Sample-wise deconvolution methods have been developed to estimate cell-type proportions and gene expressions in bulk-tissue samples. However, the performance of these methods and their biological applications has not been evaluated, particularly on human brain transcriptomic data. Here, nine deconvolution methods were evaluated with sample-matched data from bulk-tissue RNAseq, single-cell/nuclei (sc/sn) RNAseq, and immunohistochemistry. A total of 1,130,767 nuclei/cells from 149 adult postmortem brains and 72 organoid samples were used. The results showed the best performance of dtangle for estimating cell proportions and bMIND for estimating sample-wise cell-type gene expression. For eight brain cell types, 25,273 cell-type eQTLs were identified with deconvoluted expressions (decon-eQTLs). The results showed that decon-eQTLs explained more schizophrenia GWAS heritability than bulk-tissue or single-cell eQTLs alone. Differential gene expression associated with multiple phenotypes were also examined using the deconvoluted data. Our findings, which were replicated in bulk-tissue RNAseq and sc/snRNAseq data, provided new insights into the biological applications of deconvoluted data.

Cell-Type-Specific Profibrotic Scores across Multi-Organ Systems Predict Cancer Prognosis

Simple Summary

Fibrosis is a major player and contributor in the tumor microenvironment. Profibrotic changes precede the early development and establishment of a variety of human diseases, such as fibrosis and cancer. Being able to measure such early signals at the single cell level is critically useful for identifying new mechanisms and potential drug targets for a wide range of diseases. This study was designed to computationally identify profibrotic cell populations using single-cell transcriptomic data and to identify gene signatures that could predict cancer prognosis.

Abstract

Fibrosis is a major cause of mortality. Key profibrotic mechanisms are common pathways involved in tumorigenesis. Characterizing the profibrotic phenotype will help reveal the underlying mechanisms of early development and progression of a variety of human diseases, such as fibrosis and cancer. Fibroblasts have been center stage in response to various stimuli, such as viral infections. However, a comprehensive catalog of cell types involved in this process is currently lacking. Here, we deployed single-cell transcriptomic data across multi-organ systems (i.e., heart, kidney, liver, and lung) to identify novel profibrotic cell populations based on ECM pathway activity at single-cell resolution. In addition to fibroblasts, we also reported that epithelial, endothelial, myeloid, natural killer T, and secretory cells, as well as proximal convoluted tubule cells of the nephron, were significantly actively involved. Cell-type-specific gene signatures were enriched in viral infection pathways, enhanced glycolysis, and carcinogenesis, among others; they were validated using independent datasets in this study. By projecting the signatures into bulk TCGA tumor samples, we could predict prognosis in the patients using profibrotic scores. Our profibrotic cellular phenotype is useful for identifying new mechanisms and potential drug targets at the cell-type level for a wide range of diseases involved in ECM pathway activation.

DeepFun: a deep learning sequence-based model to decipher non-coding variant effect in a tissue- and cell type-specific manner

More than 90% of the genetic variants identified from genome-wide association studies (GWAS) are located in non-coding regions of the human genome. Here, we present a user-friendly web server, DeepFun (https://bioinfo.uth.edu/deepfun/), to assess the functional activity of non-coding genetic variants. This new server is built on a convolutional neural network (CNN) framework that has been extensively evaluated. Specifically, we collected chromatin profiles from ENCODE and Roadmap projects to construct the feature space, including 1548 DNase I accessibility, 1536 histone mark, and 4795 transcription factor binding profiles covering 225 tissues or cell types. With such comprehensive epigenomics annotations, DeepFun expands the functionality of existing non-coding variant prioritizing tools to provide a more specific functional assessment on non-coding variants in a tissue- and cell type-specific manner. By using the datasets from various GWAS studies, we conducted independent validations and demonstrated the functions of the DeepFun web server in predicting the effect of a non-coding variant in a specific tissue or cell type, as well as visualizing the potential motifs in the region around variants. We expect our server will be widely used in genetics, functional genomics, and disease studies.

Cell-type deconvolution analysis identifies cancer-associated myofibroblast component as a poor prognostic factor in multiple cancer types

Cancer-associated fibroblasts (CAFs) constitute a prominent component of the tumor microenvironment and play critical roles in cancer progression and drug resistance. Although recent studies indicate CAFs may consist of several CAF subtypes, the breadth of CAF heterogeneity and functional roles of CAF subtypes in cancer progression remain unclear. In this study, we implemented a cell-type deconvolutional approach to comprehensively characterize cell-type alternations across 18 cancer types from The Cancer Genome Atlas (TCGA). Pan-cancer survival analysis using deconvoluted CAF subtypes revealed myofibroblastic CAF (myCAF) composition as a poor prognostic factor in nine cancer types. Patients with higher myCAF compositions tend to have worse response to six antineoplastic drugs predicted by a lncRNA-based Elastic Net prediction model (LENP). In addition, integrative mutational analysis identified 14 and 413 genes associated with the differentiation degree of myCAF and inflammatory CAF (iCAF), respectively, with significant enrichment of genes involved in fibroblast and extracellular matrix (ECM)-related pathways. In summary, our findings systematically illustrated the complex roles of CAF subtypes in patient prognosis and drug response, and identified putative driver genes in CAF-subtype differentiation. These results provided novel therapeutic perspectives for targeting CAF subtypes in tumor microenvironment and arranging treatment scheme based on the CAF compositions in different cancer types.