A computational method for direct imputation of cell type-specific expression profiles and cellular compositions from bulk-tissue RNA-Seq in brain disorders

The importance of cell type-specific gene expression in disease-relevant tissues is increasingly recognized in genetic studies of complex diseases. However, most gene expression studies are conducted on bulk tissues, without examining cell type-specific expression profiles. Several computational methods are available for cell type deconvolution (i.e. inference of cellular composition) from bulk RNA-Seq data, but few of them impute cell type-specific expression profiles. We hypothesize that with external prior information such as single cell RNA-seq and population-wide expression profiles, it can be computationally tractable to estimate both cellular composition and cell type-specific expression from bulk RNA-Seq data. Here we introduce CellR, which addresses cross-individual gene expression variations to adjust the weights of cell-specific gene markers. It then transforms the deconvolution problem into a linear programming model while taking into account inter/intra cellular correlations and uses a multi-variate stochastic search algorithm to estimate the cell type-specific expression profiles. Analyses on several complex diseases such as schizophrenia, Alzheimer’s disease, Huntington’s disease and type 2 diabetes validated the efficiency of CellR, while revealing how specific cell types contribute to different diseases. In summary, CellR compares favorably against competing approaches, enabling cell type-specific re-analysis of gene expression data on bulk tissues in complex diseases.

Integrated analysis on transcriptome and behaviors defines HTT repeat-dependent network modules in Huntington's disease

Huntington's disease (HD) is caused by a CAG repeat expansion in the huntingtin (HTT) gene. Knock-in mice carrying a CAG repeat-expanded Htt will develop HD phenotypes. Previous studies suggested dysregulated molecular networks in a CAG length genotype- and the age-dependent manner in brain tissues from knock-in mice carrying expanded Htt CAG repeats. Furthermore, a large-scale phenome analysis defined a behavioral signature for HD genotype in knock-in mice carrying expanded Htt CAG repeats. However, an integrated analysis correlating phenotype features with genotypes (CAG repeat expansions) was not conducted previously. In this study, we revealed the landscape of the behavioral features and gene expression correlations based on 445 mRNA samples and 445 microRNA samples, together with behavioral features (396 PhenoCube behaviors and 111 NeuroCube behaviors) in Htt CAG-knock-in mice. We identified 37 behavioral features that were significantly associated with CAG repeat length including the number of steps and hind limb stand duration. The behavioral features were associated with several gene coexpression groups involved in neuronal dysfunctions, which were also supported by the single-cell RNA sequencing data in the striatum and the spatial gene expression in the brain. We also identified 15 chemicals with significant responses for genes with enriched behavioral features, most of them are agonist or antagonist for dopamine receptors and serotonin receptors used for neurology/psychiatry. Our study provides further evidence that abnormal neuronal signal transduction in the striatum plays an important role in causing HD-related phenotypic behaviors and provided rich information for the further pharmacotherapeutic intervention possibility for HD.

Cell Type-Specific Annotation and Fine Mapping of Variants Associated With Brain Disorders

Common genetic variants confer susceptibility to a large number of complex brain disorders. Given that such variants predominantly localize in non-coding regions of the human genome, there is a significant challenge to predict and characterize their functional consequences. More importantly, most available computational methods, generally defined as context-free methods, output prediction scores regarding the functionality of genetic variants irrespective of the context, i.e., the tissue or cell-type affected by a disease, limiting the ability to predict the functional consequences of common variants on brain disorders. In this study, we introduce a comparative multi-step pipeline to investigate the relative effectiveness of context-specific and context-free approaches to prioritize disease causal variants. As an experimental case, we focused on schizophrenia (SCZ), a debilitating neuropsychiatric disease for which a large number of susceptibility variants is identified from genome-wide association studies. We tested over two dozen available methods and examined potential associations between the cell/tissue-specific mapping scores and open chromatin accessibility, and provided a prioritized map of SCZ risk loci for in vitro or in-vivo functional analysis. We found extensive differences between context-free and tissue-specific approaches and showed how they may play complementary roles. As a proof of concept, we found a few sets of genes, through a consensus mapping of both categories, including FURIN to be among the top hits. We showed that the genetic variants in this gene and related genes collectively dysregulate gene expression patterns in stem cell-derived neurons and characterize SCZ phenotypic manifestations, while genes which were not shared among highly prioritized candidates in both approaches did not demonstrate such characteristics. In conclusion, by combining context-free and tissue-specific predictions, our pipeline enables prioritization of the most likely disease-causal common variants in complex brain disorders.

Cell-Type-Specific Proteogenomic Signal Diffusion for Integrating Multi-Omics Data Predicts Novel Schizophrenia Risk Genes

Accumulation of diverse types of omics data on schizophrenia (SCZ) requires a systems approach to model the interplay between genome, transcriptome, and proteome. We introduce Markov affinity-based proteogenomic signal diffusion (MAPSD), a method to model intra-cellular protein trafficking paradigms and tissue-wise single-cell protein abundances. MAPSD integrates multi-omics data to amplify the signals at SCZ risk loci with small effect sizes, and reveal convergent disease-associated gene modules in the brain. We predicted a set of high-confidence SCZ risk loci followed by characterizing the subcellular localization of proteins encoded by candidate SCZ risk genes, and illustrated that most are enriched in neuronal cells in the cerebral cortex as well as Purkinje cells in the cerebellum. We demonstrated how the identified genes may be involved in neurodevelopment, how they may alter SCZ-related biological pathways, and how they facilitate drug repurposing. MAPSD is applicable in other polygenic diseases and can facilitate our understanding of disease mechanisms.

The proteome and its dynamics: A missing piece for integrative multi-omics in schizophrenia

The complex and heterogeneous pathophysiology of schizophrenia can be deconstructed by integration of large-scale datasets encompassing genes through behavioral phenotypes. Genome-wide datasets are now available for genetic, epigenetic and transcriptomic variations in schizophrenia, which are then analyzed by newly devised systems biology algorithms. A missing piece, however, is the inclusion of information on the proteome and its dynamics in schizophrenia. Proteomics has lagged behind omics of the genome, transcriptome and epigenome since analytic platforms were relatively less robust for proteins. There has been remarkable progress, however, in the instrumentation of liquid chromatography (LC) and mass spectrometry (MS) (LCMS), experimental paradigms and bioinformatics of the proteome. Here, we present a summary of methodological innovations of recent years in MS based proteomics and the power of new generation proteomics, review proteomics studies that have been conducted in schizophrenia to date, and propose how such data can be analyzed and integrated with other omics results. The function of a protein is determined by multiple molecular properties, i.e., subcellular localization, posttranslational modification (PTMs) and protein-protein interactions (PPIs). Incorporation of these properties poses additional challenges in proteomics and their integration with other omics; yet is a critical next step to close the loop of multi-omics integration. In sum, the recent advent of high-throughput proteome characterization technologies and novel mathematical approaches enable us to incorporate functional properties of the proteome to offer a comprehensive multi-omics based understanding of schizophrenia pathophysiology.

Deconvolution of transcriptional networks identifies TCF4 as a master regulator in schizophrenia

Applying tissue-specific deconvolution of transcriptional networks to identify their master regulators (MRs) in neuropsychiatric disorders has been largely unexplored. Here, using two schizophrenia (SCZ) case-control RNA-seq datasets, one on postmortem dorsolateral prefrontal cortex (DLPFC) and another on cultured olfactory neuroepithelium, we deconvolved the transcriptional networks and identified TCF4 as a top candidate MR that may be dysregulated in SCZ. We validated TCF4 as a MR through enrichment analysis of TCF4-binding sites in induced pluripotent stem cell (hiPSC)-derived neurons and in neuroblastoma cells. We further validated the predicted TCF4 targets by knocking down TCF4 in hiPSC-derived neural progenitor cells (NPCs) and glutamatergic neurons (Glut_Ns). The perturbed TCF4 gene network in NPCs was more enriched for pathways involved in neuronal activity and SCZ-associated risk genes, compared to Glut_Ns. Our results suggest that TCF4 may serve as a MR of a gene network dysregulated in SCZ at early stages of neurodevelopment.

Uncovering the Functional Link Between SHANK3 Deletions and Deficiency in Neurodevelopment Using iPSC-Derived Human Neurons

SHANK3 mutations, including de novo deletions, have been associated with autism spectrum disorders (ASD). However, the effects of SHANK3 loss of function on neurodevelopment remain poorly understood. Here we generated human induced pluripotent stem cells (iPSC) in vitro, followed by neuro-differentiation and lentivirus-mediated shRNA expression to evaluate how SHANK3 knockdown affects the in vitro neurodevelopmental process at multiple time points (up to 4 weeks). We found that SHANK3 knockdown impaired both early stage of neuronal development and mature neuronal function, as demonstrated by a reduction in neuronal soma size, growth cone area, neurite length and branch numbers. Notably, electrophysiology analyses showed defects in excitatory and inhibitory synaptic transmission. Furthermore, transcriptome analyses revealed that multiple biological pathways related to neuron projection, motility and regulation of neurogenesis were disrupted in cells with SHANK3 knockdown. In conclusion, utilizing a human iPSC-based neural induction model, this study presented combined morphological, electrophysiological and transcription evidence that support that SHANK3 as an intrinsic, cell autonomous factor that controls cellular function development in human neurons.

Graph-based semi-supervised learning with genomic data integration using condition-responsive genes applied to phenotype classification

Objective

Data integration methods that combine data from different molecular levels such as genome, epigenome, transcriptome, etc., have received a great deal of interest in the past few years. It has been demonstrated that the synergistic effects of different biological data types can boost learning capabilities and lead to a better understanding of the underlying interactions among molecular levels.

Methods

In this paper we present a graph-based semi-supervised classification algorithm that incorporates latent biological knowledge in the form of biological pathways with gene expression and DNA methylation data. The process of graph construction from biological pathways is based on detecting condition-responsive genes, where 3 sets of genes are finally extracted: all condition responsive genes, high-frequency condition-responsive genes, and P-value-filtered genes.

Results

The proposed approach is applied to ovarian cancer data downloaded from the Human Genome Atlas. Extensive numerical experiments demonstrate superior performance of the proposed approach compared to other state-of-the-art algorithms, including the latest graph-based classification techniques.

Conclusions

Simulation results demonstrate that integrating various data types enhances classification performance and leads to a better understanding of interrelations between diverse omics data types. The proposed approach outperforms many of the state-of-the-art data integration algorithms.

Transcriptional network analysis on brains reveals a potential regulatory role of PPP1R3F in autism spectrum disorders

OBJECTIVE

This study aims at identifying master regulators of transcriptional networks in autism spectrum disorders (ASDs).

RESULTS

With two sets of independent RNA-Seq data generated on cerebellum from patients with ASDs and control subjects (N = 39 and 45 for set 1, N = 24 and 38 for set 2, respectively), we carried out a network deconvolution of transcriptomic data, followed by virtual protein activity analysis. We identified PPP1R3F (Protein Phosphatase 1 Regulatory Subunit 3F) as a candidate master regulator affecting a large body of downstream genes that are associated with the disease phenotype. Pathway analysis on the identified targets of PPP1R3F in both datasets indicated alteration of endocytosis pathway. Despite a limited sample size, our study represents one of the first applications of network deconvolution approach to brain transcriptomic data to generate hypotheses that may be further validated by large-scale studies.

Next-generation sequencing in drug development: target identification and genetically stratified clinical trials

Next-generation sequencing (NGS) enabled high-throughput analysis of genotype-phenotype relationships on human populations, ushering in a new era of genetics-informed drug development. The year 2017 was remarkable, with the first FDA-approved gene therapy for cancer (Kymriah™) and for inherited diseases (LUXTURNA™), the first multiplex NGS panel for companion diagnostics (MSK-IMPACT™) and the first drug targeting a genetic signature rather than a disease (Keytruda^®). We envision that population-scale NGS with paired electronic health records (EHRs) will become a routine measure in the drug development process for the identification of novel drug targets, and that genetically stratified clinical trials could be widely adopted to improve power in precision-medicine-guided drug development.

Sparse Pathway-Induced Dynamic Network Biomarker Discovery for Early Warning Signal Detection in Complex Diseases

In many complex diseases, the transition process from the healthy stage to the catastrophic stage does not occur gradually. Recent studies indicate that the initiation and progression of such diseases are comprised of three steps including healthy stage, pre-disease stage, and disease stage. It has been demonstrated that a certain set of trajectories can be observed in the genetic signatures at the molecular level, which might be used to detect the pre-disease stage and to take necessary medical interventions. In this paper, we propose two optimization-based algorithms for extracting the dynamic network biomarkers responsible for catastrophic transition into the disease stage, and to open new horizons to reverse the disease progression at an early stage through pinpointing molecular signatures provided by high-throughput microarray data. The first algorithm relies on meta-heuristic intelligent search to characterize dynamic network biomarkers represented as a complete graph. The second algorithm induces sparsity on the adjacency matrix of the genes by taking into account the biological signaling and metabolic pathways, since not all the genes in the ineractome are biologically linked. Comprehensive numerical and meta-analytical experiments verify the effectiveness of the results of the proposed approaches in terms of network size, biological meaningfulness, and verifiability.

Alpha-plane based automatic general type-2 fuzzy clustering based on simulated annealing meta-heuristic algorithm for analyzing gene expression data

This paper considers microarray gene expression data clustering using a novel two stage meta-heuristic algorithm based on the concept of α-planes in general type-2 fuzzy sets. The main aim of this research is to present a powerful data clustering approach capable of dealing with highly uncertain environments. In this regard, first, a new objective function using α-planes for general type-2 fuzzy c-means clustering algorithm is represented. Then, based on the philosophy of the meta-heuristic optimization framework 'Simulated Annealing', a two stage optimization algorithm is proposed. The first stage of the proposed approach is devoted to the annealing process accompanied by its proposed perturbation mechanisms. After termination of the first stage, its output is inserted to the second stage where it is checked with other possible local optima through a heuristic algorithm. The output of this stage is then re-entered to the first stage until no better solution is obtained. The proposed approach has been evaluated using several synthesized datasets and three microarray gene expression datasets. Extensive experiments demonstrate the capabilities of the proposed approach compared with some of the state-of-the-art techniques in the literature.

A hybrid fuzzy-ontology based intelligent system to determine level of severity and treatment recommendation for Benign Prostatic Hyperplasia

This paper deals with application of fuzzy intelligent systems in diagnosing severity level and recommending appropriate therapies for patients having Benign Prostatic Hyperplasia. Such an intelligent system can have remarkable impacts on correct diagnosis of the disease and reducing risk of mortality. This system captures various factors from the patients using two modules. The first module determines severity level of the Benign Prostatic Hyperplasia and the second module, which is a decision making unit, obtains output of the first module accompanied by some external knowledge and makes an appropriate treatment decision based on its ontology model and a fuzzy type-1 system. In order to validate efficiency and accuracy of the developed system, a case study is conducted by 44 participants. Then the results are compared with the recommendations of a panel of experts on the experimental data. Then precision and accuracy of the results were investigated based on a statistical analysis.