Estimating gene regulatory networks with pandaR

Evaluation of altered cell-cell communication between glia and neurons in the hippocampus of 3xTg-AD mice at two time points

Alzheimer's disease (AD) is the most common form of dementia and is characterized by progressive memory loss and cognitive decline, affecting behavior, speech, and motor abilities. The neuropathology of AD includes the formation of extracellular amyloid-β plaques and intracellular neurofibrillary tangles of phosphorylated tau, along with neuronal loss. Although neuronal loss is an AD hallmark, cell-cell communication between neuronal and non-neuronal cell populations maintains neuronal health and brain homeostasis. To study changes in cell-cell communication during disease progression, we performed snRNA-sequencing of the hippocampus from female 3xTg-AD and wild-type littermates at 6 and 12 months. We inferred differential cell-cell communication between 3xTg-AD and wild-type mice across time points and between senders (astrocytes, microglia, oligodendrocytes, and OPCs) and receivers (excitatory and inhibitory neurons) of interest. We also assessed the downstream effects of altered glia-neuron communication using pseudobulk differential gene expression, functional enrichment, and gene regulatory analyses. We found that glia-neuron communication is increasingly dysregulated in 12-month 3xTg-AD mice. We also identified 23 AD-associated ligand-receptor pairs that are upregulated in the 12-month-old 3xTg-AD hippocampus. Our results suggest increased AD association of interactions originating from microglia. Signaling mediators were not significantly differentially expressed but showed altered gene regulation and transcription factor activity. Our findings indicate that altered glia-neuron communication is increasingly dysregulated and affects the gene regulatory mechanisms in neurons of 12-month-old 3xTg-AD mice.

Emerging Signatures of Hematological Malignancies from Gene Expression and Transcription Factor-Gene Regulations

Hematological malignancies are a diverse group of cancers developing in the peripheral blood, the bone marrow or the lymphatic system. Due to their heterogeneity, the identification of novel and advanced molecular signatures is essential for enhancing their characterization and facilitate its translation to new pharmaceutical solutions and eventually to clinical applications. In this study, we collected publicly available microarray data for more than five thousand subjects, across thirteen hematological malignancies. Using PANDA to estimate gene regulatory networks (GRNs), we performed hierarchical clustering and network analysis to explore transcription factor (TF) interactions and their implications on biological pathways. Our findings reveal distinct clustering patterns among leukemias and lymphomas, with notable differences in gene and TF expression profiles. Gene Set Enrichment Analysis (GSEA) identified 57 significantly enriched KEGG pathways, highlighting both common and unique biological processes across HMs. We also identified potential drug targets within these pathways, emphasizing the role of TFs such as CEBPB and NFE2L1 in disease pathophysiology. Our comprehensive analysis enhances the understanding of the molecular landscape of HMs and suggests new avenues for targeted therapeutic strategies. These findings also motivate the adoption of regulatory networks, combined with modern biotechnological possibilities, for insightful pan-cancer exploratory studies.

Evaluation of altered cell-cell communication between glia and neurons in the hippocampus of 3xTg-AD mice at two time points

Alzheimer's disease (AD) is the most common form of dementia and is characterized by progressive memory loss and cognitive decline, affecting behavior, speech, and motor abilities. The neuropathology of AD includes the formation of extracellular amyloid-β plaque and intracellular neurofibrillary tangles of phosphorylated tau, along with neuronal loss. While neuronal loss is an AD hallmark, cell-cell communication between neuronal and non-neuronal cell populations maintains neuronal health and brain homeostasis. To study changes in cellcell communication during disease progression, we performed snRNA-sequencing of the hippocampus from female 3xTg-AD and wild-type littermates at 6 and 12 months. We inferred differential cell-cell communication between 3xTg-AD and wild-type mice across time points and between senders (astrocytes, microglia, oligodendrocytes, and OPCs) and receivers (excitatory and inhibitory neurons) of interest. We also assessed the downstream effects of altered glia-neuron communication using pseudobulk differential gene expression, functional enrichment, and gene regulatory analyses. We found that glia-neuron communication is increasingly dysregulated in 12-month 3xTg-AD mice. We also identified 23 AD-associated ligand-receptor pairs that are upregulated in the 12-month-old 3xTg-AD hippocampus. Our results suggest increased AD association of interactions originating from microglia. Signaling mediators were not significantly differentially expressed but showed altered gene regulation and TF activity. Our findings indicate that altered glia-neuron communication is increasingly dysregulated and affects the gene regulatory mechanisms in neurons of 12-month-old 3xTg-AD mice.

Cell-type-specific gene expression and regulation in the cerebral cortex and kidney of atypical Setbp1S858R Schinzel Giedion Syndrome mice

Schinzel Giedion Syndrome (SGS) is an ultra-rare autosomal dominant Mendelian disease presenting with abnormalities spanning multiple organ systems. The most notable phenotypes involve severe developmental delay, progressive brain atrophy, and drug-resistant seizures. SGS is caused by spontaneous variants in SETBP1, which encodes for the epigenetic hub SETBP1 transcription factor (TF). SETBP1 variants causing classical SGS cluster at the degron, disrupting SETBP1 protein degradation and resulting in toxic accumulation, while those located outside cause milder atypical SGS. Due to the multisystem phenotype, we evaluated gene expression and regulatory programs altered in atypical SGS by snRNA-seq of the cerebral cortex and kidney of Setbp1S858R heterozygous mice (corresponds to the human likely pathogenic SETBP1S867R variant) compared to matched wild-type mice by constructing cell-type-specific regulatory networks. Setbp1 was differentially expressed in excitatory neurons, but known SETBP1 targets were differentially expressed and regulated in many cell types. Our findings suggest molecular drivers underlying neurodevelopmental phenotypes in classical SGS also drive atypical SGS, persist after birth, and are present in the kidney. Our results indicate SETBP1's role as an epigenetic hub leads to cell-type-specific differences in TF activity, gene targeting, and regulatory rewiring. This research provides a framework for investigating cell-type-specific variant impact on gene expression and regulation.

Architecture and topologies of gene regulatory networks associated with breast cancer, adjacent normal, and normal tissues

Most cancer studies employ adjacent normal tissues to tumors (ANTs) as controls, which are not completely normal and represent a pre-cancerous state. However, the regulatory landscape of ANTs compared to tumor and non-tumor-bearing normal tissues is largely unexplored. Among cancers, breast cancer is the most commonly diagnosed cancer and a leading cause of death in women worldwide, with a lack of sufficient treatment regimens for various reasons. Hence, we aimed to gain deeper insights into normal, pre-cancerous, and cancerous regulatory systems of breast tissues towards identifying ANT and subtype-specific candidate genes. For this, we constructed and analyzed eight gene regulatory networks (GRNs), including five subtypes (viz., Basal, Her2, Luminal A, Luminal B, and Normal-Like), one ANT, and two normal tissue networks. Whereas several topological properties of these GRNs enabled us to identify tumor-related features of ANT, escape velocity centrality (EVC+) identified 24 functionally significant common genes, including well-known genes such as E2F1, FOXA1, JUN, BRCA1, GATA3, ERBB2, and ERBB3 across all six tissues including subtypes and ANT. Similarly, the EVC+ also helped us to identify tissue-specific key genes (Basal: 18, Her2: 6, Luminal A: 5, Luminal B: 5, Normal-Like: 2, and ANT: 7). Additionally, differentially correlated switching gene pairs along with functional, pathway, and disease annotations highlighted the cancer-associated role of these genes. In a nutshell, the present study revealed ANT and subtype-specific regulatory features and key candidate genes, which can be explored further using in vitro and in vivo experiments for better and effective disease management at an early stage.

The Network Zoo: a multilingual package for the inference and analysis of gene regulatory networks

Inference and analysis of gene regulatory networks (GRNs) require software that integrates multi-omic data from various sources. The Network Zoo (netZoo; netzoo.github.io) is a collection of open-source methods to infer GRNs, conduct differential network analyses, estimate community structure, and explore the transitions between biological states. The netZoo builds on our ongoing development of network methods, harmonizing the implementations in various computing languages and between methods to allow better integration of these tools into analytical pipelines. We demonstrate the utility using multi-omic data from the Cancer Cell Line Encyclopedia. We will continue to expand the netZoo to incorporate additional methods.

Network-based approaches for modeling disease regulation and progression

Molecular interaction networks lay the foundation for studying how biological functions are controlled by the complex interplay of genes and proteins. Investigating perturbed processes using biological networks has been instrumental in uncovering mechanisms that underlie complex disease phenotypes. Rapid advances in omics technologies have prompted the generation of high-throughput datasets, enabling large-scale, network-based analyses. Consequently, various modeling techniques, including network enrichment, differential network extraction, and network inference, have proven to be useful for gaining new mechanistic insights. We provide an overview of recent network-based methods and their core ideas to facilitate the discovery of disease modules or candidate mechanisms. Knowledge generated from these computational efforts will benefit biomedical research, especially drug development and precision medicine. We further discuss current challenges and provide perspectives in the field, highlighting the need for more integrative and dynamic network approaches to model disease development and progression.

Identifying predictive signalling networks for Vedolizumab response in ulcerative colitis

BACKGROUND

In ulcerative colitis (UC), the molecular mechanisms that drive disease development and patient response to therapy are not well understood. A significant proportion of patients with UC fail to respond adequately to biologic therapy. Therefore, there is an unmet need for biomarkers that can predict patients' responsiveness to the available UC therapies as well as ascertain the most effective individualised therapy. Our study focused on identifying predictive signalling pathways that predict anti-integrin therapy response in patients with UC.

METHODS

We retrieved and pre-processed two publicly accessible gene expression datasets (GSE73661 and GSE72819) of UC patients treated with anti-integrin therapies: (1) 12 non-IBD controls and 41 UC patients treated with Vedolizumab therapy, and (2) 70 samples with 58 non-responder and 12 responder UC patient samples treated with Etrolizumab therapy without non-IBD controls. We used a diffusion-based signalling model which is mainly focused on the T-cell receptor signalling network. The diffusion model uses network connectivity between receptors and transcription factors.

RESULTS

The network diffusion scores were able to separate VDZ responder and non-responder patients before treatment better than the original gene expression. On both anti-integrin treatment datasets, the diffusion model demonstrated high predictive performance for discriminating responders from non-responders in comparison with 'nnet'. We have found 48 receptor-TF pairs identified as the best predictors for VDZ therapy response with AUC ≥ 0.76. Among these receptor-TF predictors pairs, FFAR2-NRF1, FFAR2-RELB, FFAR2-EGR1, and FFAR2-NFKB1 are the top best predictors. For Etrolizumab, we have identified 40 best receptor-TF pairs and CD40-NFKB2 as the best predictor receptor-TF pair (AUC = 0.72). We also identified subnetworks that highlight the network interactions, connecting receptors and transcription factors involved in cytokine and fatty acid signalling. The findings suggest that anti-integrin therapy responses in cytokine and fatty acid signalling can stratify UC patient subgroups.

CONCLUSIONS

We identified signalling pathways that may predict the efficacy of anti-integrin therapy in UC patients and personalised therapy alternatives. Our results may lead to the advancement of a promising clinical decision-making tool for the stratification of UC patients.

Identifying anti-TNF response biomarkers in ulcerative colitis using a diffusion-based signalling model

Motivation

Resistance to anti-TNF therapy in subgroups of ulcerative colitis (UC) patients is a major challenge and incurs significant treatment costs. Identification of patients at risk of nonresponse to anti-TNF is of major clinical importance. To date, no quantitative computational framework exists to develop a complex biomarker for the prognosis of UC treatment. Modelling patient-wise receptor to transcription factor (TF) network connectivity may enable personalized treatment.

Results

We present an approach for quantitative diffusion analysis between receptors and TFs using gene expression data. Key TFs were identified using pandaR. Network connectivities between immune-specific receptor-TF pairs were quantified using network diffusion in UC patients and controls. The patient-specific network could be considered a complex biomarker that separates anti-TNF treatment-resistant and responder patients both in the gene expression dataset used for model development and separate independent test datasets. The model was further validated in rheumatoid arthritis where it successfully discriminated resistant and responder patients to tocilizumab treatment. Our model may contribute to prognostic biomarkers that may identify treatment-resistant and responder subpopulations of UC patients.

Availability and implementation

Software is available at https://github.com/Amy3100/receptor2tfDiffusion.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Different Biological Pathways Between Good and Poor Inhaled Corticosteroid Responses in Asthma

Gene regulatory networks address how transcription factors (TFs) and their regulatory roles in gene expression determine the responsiveness to anti-asthma therapy. The purpose of this study was to assess gene regulatory networks of adult patients with asthma who showed good or poor lung function improvements in response to inhaled corticosteroids (ICSs). A total of 47 patients with asthma were recruited and classified as good responders (GRs) and poor responders (PRs) based on their responses to ICSs. Genome-wide gene expression was measured using peripheral blood mononuclear cells obtained in a stable state. We used Passing Attributes between Networks for Data Assimilations to construct the gene regulatory networks associated with GRs and PRs to ICSs. We identified the top-10 TFs that showed large differences in high-confidence edges between the GR and PR aggregate networks. These top-10 TFs and their differentially-connected genes in the PR and GR aggregate networks were significantly enriched in distinct biological pathways, such as TGF-β signaling, cell cycle, and IL-4 and IL-13 signaling pathways. We identified multiple TFs and related biological pathways influencing ICS responses in asthma. Our results provide potential targets to overcome insensitivity to corticosteroids in patients with asthma.

Understanding the Molecular Mechanisms of Asthma through Transcriptomics

The transcriptome represents the complete set of RNA transcripts that are produced by the genome under a specific circumstance or in a specific cell. High-throughput methods, including microarray and bulk RNA sequencing, as well as recent advances in biostatistics based on machine learning approaches provides a quick and effective way of identifying novel genes and pathways related to asthma, which is a heterogeneous disease with diverse pathophysiological mechanisms. In this manuscript, we briefly review how to analyze transcriptome data and then provide a summary of recent transcriptome studies focusing on asthma pathogenesis and asthma drug responses. Studies reviewed here are classified into 2 classes based on the tissues utilized: blood and airway cells.

Sufu- and Spop-mediated downregulation of Hedgehog signaling promotes beta cell differentiation through organ-specific niche signals

Human embryonic stem cell-derived beta cells offer a promising cell-based therapy for diabetes. However, efficient stem cell to beta cell differentiation has proven difficult, possibly due to the lack of cross-talk with the appropriate mesenchymal niche. To define organ-specific niche signals, we isolated pancreatic and gastrointestinal stromal cells, and analyzed their gene expression during development. Our genetic studies reveal the importance of tightly regulated Hedgehog signaling in the pancreatic mesenchyme: inactivation of mesenchymal signaling leads to annular pancreas, whereas stroma-specific activation of signaling via loss of Hedgehog regulators, Sufu and Spop, impairs pancreatic growth and beta cell genesis. Genetic rescue and transcriptome analyses show that these Sufu and Spop knockout defects occur through Gli2-mediated activation of gastrointestinal stromal signals such as Wnt ligands. Importantly, inhibition of Wnt signaling in organoid and human stem cell cultures significantly promotes insulin-producing cell generation, altogether revealing the requirement for organ-specific regulation of stromal niche signals.

Tissue-aware RNA-Seq processing and normalization for heterogeneous and sparse data

BACKGROUND

Although ultrahigh-throughput RNA-Sequencing has become the dominant technology for genome-wide transcriptional profiling, the vast majority of RNA-Seq studies typically profile only tens of samples, and most analytical pipelines are optimized for these smaller studies. However, projects are generating ever-larger data sets comprising RNA-Seq data from hundreds or thousands of samples, often collected at multiple centers and from diverse tissues. These complex data sets present significant analytical challenges due to batch and tissue effects, but provide the opportunity to revisit the assumptions and methods that we use to preprocess, normalize, and filter RNA-Seq data - critical first steps for any subsequent analysis.

RESULTS

We find that analysis of large RNA-Seq data sets requires both careful quality control and the need to account for sparsity due to the heterogeneity intrinsic in multi-group studies. We developed Yet Another RNA Normalization software pipeline (YARN), that includes quality control and preprocessing, gene filtering, and normalization steps designed to facilitate downstream analysis of large, heterogeneous RNA-Seq data sets and we demonstrate its use with data from the Genotype-Tissue Expression (GTEx) project.

CONCLUSIONS

An R package instantiating YARN is available at http://bioconductor.org/packages/yarn .