Reverse engineering of regulatory networks in human B cells

Approaches for Benchmarking Single-Cell Gene Regulatory Network Methods

Gene regulatory networks are powerful tools for modeling genetic interactions that control the expression of genes driving cell differentiation, and single-cell sequencing offers a unique opportunity to build these networks with high-resolution genomic data. There are many proposed computational methods to build these networks using single-cell data, and different approaches are used to benchmark these methods. However, a comprehensive discussion specifically focusing on benchmarking approaches is missing. In this article, we lay the GRN terminology, present an overview of common gold-standard studies and data sets, and define the performance metrics for benchmarking network construction methodologies. We also point out the advantages and limitations of different benchmarking approaches, suggest alternative ground truth data sets that can be used for benchmarking, and specify additional considerations in this context.

PPM1D/Wip1 is amplified, overexpressed, and mutated in human non-Hodgkin's lymphomas

BACKGROUND

Wip1, is a p53-dependent Ser/Thr phosphatase involved in the timely termination of DDR. The PPM1D gene encoding Wip1 is deregulated and thus gained an oncogene character in common human solid tumors and cell lines. This study assessed the oncogenic potential of the PPM1D gene in human non- Hodgkin's lymphomas (NHL), the most common hematological malignancy worldwide.

METHODS AND RESULTS

FFPE human lymphoid hyperplasia (LH) (n = 17) and NHL tumor lymph node samples (n = 65) and human NHL cell lines were used to assess the oncogenic potential of the PPM1D gene in the present study. Copy number gain and mRNA expression analysis of the PPM1D/Wip1 gene were assessed by qRT-PCR analysis. Mutational analysis of Exon 6 of the PPM1D gene was performed by PCR amplification and Sanger sequencing. Expressions of Wip1 and p53 proteins were assessed by immunohistochemistry and Western blot analysis.

CONCLUSIONS

We found that PPM1D gained gene copy number in NHL tumors by 0.7-8 times compared to the control (p < 0.01). Increased PPM1D/Wip1 gene copy number was associated with higher mRNA and protein expression in human NHL samples (p < 0.01). Overexpression of Wip1 in NHL tumors and NHL cell lines was associated with amplification level and was unaffected by p53 status. Furthermore, a heterozygous type insertion mutation was detected in exon 6 (c.1553_1554insA) of the PPM1D gene particularly in DLBCL samples. Wip1 may have oncogenic potential, perhaps playing a role in the onset and progression of human NHL. The possible significance of Wip1 overexpression to chemotherapy response in NHL remains an intriguing question that requires more exploration.

WENDY: Covariance dynamics based gene regulatory network inference

Determining gene regulatory network (GRN) structure is a central problem in biology, with a variety of inference methods available for different types of data. For a widely prevalent and challenging use case, namely single-cell gene expression data measured after intervention at multiple time points with unknown joint distributions, there is only one known specifically developed method, which does not fully utilize the rich information contained in this data type. We develop an inference method for the GRN in this case, netWork infErence by covariaNce DYnamics, dubbed WENDY. The core idea of WENDY is to model the dynamics of the covariance matrix, and solve this dynamics as an optimization problem to determine the regulatory relationships. To evaluate its effectiveness, we compare WENDY with other inference methods using synthetic data and experimental data. Our results demonstrate that WENDY performs well across different data sets.

Low cytoplasmic NUB1 protein exerts hypoxic cell death with poorer prognosis in oestrogen receptor negative breast cancer patients

Current prognostic biomarkers fall short in stratifying Oestrogen receptor (ER)-negative breast cancer patients regarding tumour progression risk at diagnosis. The role of AIPL1 in activating its tumour suppressor client protein, NEDD8 Ultimate Buster-1 (NUB1) remains unknown in cancer. Our study demonstrated how downregulated AIPL1 results in the deactivated NUB1 protein under hypoxic conditions. We examined the AIPL1-NUB1 pathwayin vitro using cell lines i.e. MCF-7, MDA-MB-231, RCC4 etc. NUB1 expression was assessed using Oncomine, and cBioPortal was performed to assess NUB1's prognostic significance in human cancers. In the John Radcliffe Hospital cohort (n = 122), immunohistochemistry analysis revealed downregulated AIPL1 (Log2 fold change=-0.28; p < 0.001) and upregulated NUB1 transcripts (Log2 fold change=0.59; p < 0.001) compared to adjacent normal tissues. In severe chronic hypoxia, multimerised AIPL1 localisedin the cytoplasm while NUB1 protein migrated to the nucleus, where the absence of NUB1 nuclear localisation led to cell cycle arrest. Biopsies showed that patients with lower cytoplasmic NUB1 expression (n = 57) had poorer overall survival compared to those with higher cytoplasmic expression (n = 57), HR=1.78; 95 % CI=1.01-3.35, p = 0.048. Low NUB1 protein levels in both normoxic and hypoxic conditions were associated with cell cycle arrest and upregulation ofp21 and p27 in breast cancer cell lines, correlating significantly withpoorer survival outcomes in all breast cancer and ER-negative breast cancer patients.

Screening of hub genes for sepsis-induced myopathy by weighted gene co-expression network analysis and protein-protein interaction network construction

Sepsis-induced myopathy is one of the serious complications of sepsis, which severely affects the respiratory and peripheral motor systems of patients, reduces their quality of life, and jeopardizes their lives, as evidenced by muscle atrophy, loss of strength, and impaired regeneration after injury. The pathogenesis of sepsis-induced myopathy is complex, mainly including cytokine action, enhances free radical production in muscle, increases muscle protein hydrolysis, and decreases skeletal muscle protein synthesis, etc. The above mechanisms have been demonstrated in existing studies. However, it is still unclear how the overall pattern of gene co-expression affects the pathological process of sepsis-induced myopathy. Therefore, we intend to identify hub genes and signaling pathways. Weighted gene co-expression network analysis was our main approach to study gene expression profiles: skeletal muscle transcriptome in ICU patients with sepsis-induced multi-organ failure (GSE13205). After data pre-processing, about 15,181 genes were used to identify 13 co-expression modules. Then, 16 genes (FEM1B, KLHDC3, GPX3, NIFK, GNL2, EBNA1BP2, PES1, FBP2, PFKP, BYSL, HEATR1, WDR75, TBL3, and WDR43) were selected as the hub genes including 3 up-regulated genes and 13 down-regulated genes. Then, Gene Set Enrichment Analysis was performed to show that the hub genes were closely associated with skeletal muscle dysfunction, necrotic and apoptotic skeletal myoblasts, and apoptosis in sepsis-induced myopathy. Overall, 16 candidate biomarkers were certified as reliable features for more in-depth exploration of sepsis-induced myopathy in basic and clinical studies.

Constructing the dynamic transcriptional regulatory networks to identify phenotype-specific transcription regulators

The transcriptional regulatory network (TRN) is a graph framework that helps understand the complex transcriptional regulation mechanisms in the transcription process. Identifying the phenotype-specific transcription regulators is vital to reveal the functional roles of transcription elements in associating the specific phenotypes. Although many methods have been developed towards detecting the phenotype-specific transcription elements based on the static TRN in the past decade, most of them are not satisfactory for elucidating the phenotype-related functional roles of transcription regulators in multiple levels, as the dynamic characteristics of transcription regulators are usually ignored in static models. In this study, we introduce a novel framework called DTGN to identify the phenotype-specific transcription factors (TFs) and pathways by constructing dynamic TRNs. We first design a graph autoencoder model to integrate the phenotype-oriented time-series gene expression data and static TRN to learn the temporal representations of genes. Then, based on the learned temporal representations of genes, we develop a statistical method to construct a series of dynamic TRNs associated with the development of specific phenotypes. Finally, we identify the phenotype-specific TFs and pathways from the constructed dynamic TRNs. Results from multiple phenotypic datasets show that the proposed DTGN framework outperforms most existing methods in identifying phenotype-specific TFs and pathways. Our framework offers a new approach to exploring the functional roles of transcription regulators that associate with specific phenotypes in a dynamic model.

GRNMOPT: Inference of gene regulatory networks based on a multi-objective optimization approach

BACKGROUND AND OBJECTIVE

The reconstruction of gene regulatory networks (GRNs) stands as a vital approach in deciphering complex biological processes. The application of nonlinear ordinary differential equations (ODEs) models has demonstrated considerable efficacy in predicting GRNs. Notably, the decay rate and time delay are pivotal in authentic gene regulation, yet their systematic determination in ODEs models remains underexplored. The development of a comprehensive optimization framework for the effective estimation of these key parameters is essential for accurate GRN inference.

METHOD

This study introduces GRNMOPT, an innovative methodology for inferring GRNs from time-series and steady-state data. GRNMOPT employs a combined use of decay rate and time delay in constructing ODEs models to authentically represent gene regulatory processes. It incorporates a multi-objective optimization approach, optimizing decay rate and time delay concurrently to derive Pareto optimal sets for these factors, thereby maximizing accuracy metrics such as AUROC (Area Under the Receiver Operating Characteristic curve) and AUPR (Area Under the Precision-Recall curve). Additionally, the use of XGBoost for calculating feature importance aids in identifying potential regulatory gene links.

RESULTS

Comprehensive experimental evaluations on two simulated datasets from DREAM4 and three real gene expression datasets (Yeast, In vivo Reverse-engineering and Modeling Assessment [IRMA], and Escherichia coli [E. coli]) reveal that GRNMOPT performs commendably across varying network scales. Furthermore, cross-validation experiments substantiate the robustness of GRNMOPT.

CONCLUSION

We propose a novel approach called GRNMOPT to infer GRNs based on a multi-objective optimization framework, which effectively improves inference accuracy and provides a powerful tool for GRNs inference.

A forward genetic screen identifies Sirtuin1 as a driver of neuroendocrine prostate cancer

Although localized prostate cancer is relatively indolent, advanced prostate cancer manifests with aggressive and often lethal variants, including neuroendocrine prostate cancer (NEPC). To identify drivers of aggressive prostate cancer, we leveraged Sleeping Beauty (SB) transposon mutagenesis in a mouse model based on prostate-specific loss-of-function of Pten and Tp53 . Compared with control mice, SB mice developed more aggressive prostate tumors, with increased incidence of metastasis. Notably, a significant percentage of the SB prostate tumors display NEPC phenotypes, and the transcriptomic features of these SB mouse tumors recapitulated those of human NEPC. We identified common SB transposon insertion sites (CIS) and prioritized associated CIS-genes differentially expressed in NEPC versus non-NEPC SB tumors. Integrated analysis of CIS-genes encoding for proteins representing upstream, post-translational modulators of master regulators controlling the transcriptional state of SB -mouse and human NEPC tumors identified sirtuin 1 ( Sirt1 ) as a candidate mechanistic determinant of NEPC. Gain-of-function studies in human prostate cancer cell lines confirmed that SIRT1 promotes NEPC, while its loss-of-function or pharmacological inhibition abrogates NEPC. This integrative analysis is generalizable and can be used to identify novel cancer drivers for other malignancies.

Summary

Using an unbiased forward mutagenesis screen in an autochthonous mouse model, we have investigated mechanistic determinants of aggressive prostate cancer. SIRT1 emerged as a key regulator of neuroendocrine prostate cancer differentiation and a potential target for therapeutic intervention.

Voice from both sides: a molecular dialogue between transcriptional activators and repressors in seed-to-seedling transition and crop adaptation

High-quality seeds provide valuable nutrients to human society and ensure successful seedling establishment. During maturation, seeds accumulate storage compounds that are required to sustain seedling growth during germination. This review focuses on the epigenetic repression of the embryonic and seed maturation programs in seedlings. We begin with an extensive overview of mutants affecting these processes, illustrating the roles of core proteins and accessory components in the epigenetic machinery by comparing mutants at both phenotypic and molecular levels. We highlight how omics assays help uncover target-specific functional specialization and coordination among various epigenetic mechanisms. Furthermore, we provide an in-depth discussion on the Seed dormancy 4 (Sdr4) transcriptional corepressor family, comparing and contrasting their regulation of seed germination in the dicotyledonous species Arabidopsis and two monocotyledonous crops, rice and wheat. Finally, we compare the similarities in the activation and repression of the embryonic and seed maturation programs through a shared set of cis-regulatory elements and discuss the challenges in applying knowledge largely gained in model species to crops.

Transcriptome Profiling of Mouse Embryonic Fibroblast Spontaneous Immortalization: A Comparative Analysis

Cell immortalization, a hallmark of cancer development, is a process that cells can undergo on their path to carcinogenesis. Spontaneously immortalized mouse embryonic fibroblasts (MEFs) have been used for decades; however, changes in the global transcriptome during this process have been poorly described. In our research, we characterized the poly-A RNA transcriptome changes after spontaneous immortalization. To this end, differentially expressed genes (DEGs) were screened using DESeq2 and characterized by gene ontology enrichment analysis and protein-protein interaction (PPI) network analysis to identify the potential hub genes. In our study, we identified changes in the expression of genes involved in proliferation regulation, cell adhesion, immune response and transcriptional regulation in immortalized MEFs. In addition, we performed a comparative analysis with previously reported MEF immortalization data, where we propose a predicted gene regulatory network model in immortalized MEFs based on the altered expression of Mapk11, Cdh1, Chl1, Zic1, Hoxd10 and the novel hub genes Il6 and Itgb2.

Elucidating Compound Mechanism of Action and Polypharmacology with a Large-scale Perturbational Profile Compendium

The Mechanism of Action (MoA) of a drug is generally represented as a small, non-tissue-specific repertoire of high-affinity binding targets. Yet, drug activity and polypharmacology are increasingly associated with a broad range of off-target and tissue-specific effector proteins. To address this challenge, we have leveraged a microfluidics-based Plate-Seq technology to survey drug perturbational profiles representing >700 FDA-approved and experimental oncology drugs, in cell lines selected as high-fidelity models of 23 aggressive tumor subtypes. Built on this dataset, we implemented an efficient computational framework to define a tissue-specific protein activity landscape of these drugs and reported almost 50 million differential protein activities derived from drug perturbations vs. vehicle controls. These analyses revealed thousands of highly reproducible and novel, drug-mediated modulation of tissue-specific targets, leading to generation of a proteome-wide drug functional network, characterization of MoA-related drug clusters and off-target effects, dramatical expansion of druggable human proteome, and identification and experimental validation of novel, tissue-specific inhibitors of undruggable oncoproteins, most never reported before. The drug perturbation profile resource described here represents the first, large-scale, whole-genome-wide, RNA-Seq based dataset assembled to date, with the proposed framework, which is easily extended to elucidating the MoA of novel small-molecule libraries, facilitates mechanistic exploration of drug functions, supports systematic and quantitative approaches to precision oncology, and serves as a rich data foundation for drug discovery.

Computational single-cell methods for predicting cancer risk

Despite recent biotechnological breakthroughs, cancer risk prediction remains a formidable computational and experimental challenge. Addressing it is critical in order to improve prevention, early detection and survival rates. Here, I briefly summarize some key emerging theoretical and computational challenges as well as recent computational advances that promise to help realize the goals of cancer-risk prediction. The focus is on computational strategies based on single-cell data, in particular on bottom-up network modeling approaches that aim to estimate cancer stemness and dedifferentiation at single-cell resolution from a systems-biological perspective. I will describe two promising methods, a tissue and cell-lineage independent one based on the concept of diffusion network entropy, and a tissue and cell-lineage specific one that uses transcription factor regulons. Application of these tools to single-cell and single-nucleus RNA-seq data from stages prior to invasive cancer reveal that they can successfully delineate the heterogeneous inter-cellular cancer-risk landscape, identifying those cells that are more likely to turn cancerous. Bottom-up systems biological modeling of single-cell omic data is a novel computational analysis paradigm that promises to facilitate the development of preventive, early detection and cancer-risk prediction strategies.

Ras-dependent activation of BMAL2 regulates hypoxic metabolism in pancreatic cancer

To identify drivers of malignancy in human pancreatic ductal adenocarcinoma (PDAC), we performed regulatory network analysis on a large collection of expression profiles from laser capture microdissected samples of PDAC and benign precursors. We discovered that BMAL2 plays a role in the initiation, progression, post resection survival, and KRAS activity in PDAC. Functional analysis of BMAL2 target genes led us to hypothesize that it plays a role in regulating the response to hypoxia, a critical but poorly understood feature of PDAC physiology. Knockout of BMAL2 in multiple human PDAC cell lines revealed effects on viability and invasion, particularly under hypoxic conditions. Loss of BMAL2 also affected glycolysis and other metabolic processes. We found that BMAL2 directly regulates hypoxia-responsive target genes. We also found that BMAL2 is necessary for the stabilization of HIF1A upon exposure to hypoxia, but destabilizes HIF2A under hypoxia. These data demonstrate that BMAL2 is a master transcriptional regulator of hypoxia responses in PDAC and may serve as a long-sought molecular switch that distinguishes HIF1A- and HIF2A-dependent modes of hypoxic metabolism.

Isthmus progenitor cells contribute to homeostatic cellular turnover and support regeneration following intestinal injury

The currently accepted intestinal epithelial cell organization model proposes that Lgr5+ crypt-base columnar (CBC) cells represent the sole intestinal stem cell (ISC) compartment. However, previous studies have indicated that Lgr5+ cells are dispensable for intestinal regeneration, leading to two major hypotheses: one favoring the presence of a quiescent reserve ISC and the other calling for differentiated cell plasticity. To investigate these possibilities, we studied crypt epithelial cells in an unbiased fashion via high-resolution single-cell profiling. These studies, combined with in vivo lineage tracing, show that Lgr5 is not a specific ISC marker and that stemness potential exists beyond the crypt base and resides in the isthmus region, where undifferentiated cells participate in intestinal homeostasis and regeneration following irradiation (IR) injury. Our results provide an alternative model of intestinal epithelial cell organization, suggesting that stemness potential is not restricted to CBC cells, and neither de-differentiation nor reserve ISC are drivers of intestinal regeneration.

Comprehensive pan-cancer analysis of mitochondrial outer membrane permeabilisation activity reveals positive immunomodulation and assists in identifying potential therapeutic targets for immunotherapy resistance

BACKGROUND

Mitochondrial outer membrane permeabilisation (MOMP) plays a pivotal role in cellular death and immune activation. A deeper understanding of the impact of tumour MOMP on immunity will aid in guiding more effective immunotherapeutic strategies.

METHODS

A comprehensive pan-cancer dataset comprising 30 cancer-type transcriptomic cohorts, 20 immunotherapy transcriptomic cohorts and three immunotherapy scRNA-seq datasets was collected and analysed to determine the influence of tumour MOMP activity on clinical prognosis, immune infiltration and immunotherapy effectiveness. Leveraging 65 scRNA-Seq datasets, the MOMP signature (MOMP.Sig) was developed to accurately reflect tumour MOMP activity. The clinical predictive value of MOMP.Sig was explored through machine learning models. Integration of the MOMP.Sig model and a pan-cancer immunotherapy CRISPR screen further investigated potential targets to overcome immunotherapy resistance, which subsequently underwent clinical validation.

RESULTS

Our research revealed that elevated MOMP activity reduces mortality risk in cancer patients, drives the formation of an anti-tumour immune environment and enhances the response to immunotherapy. This finding emphasises the potential clinical application value of MOMP activity in immunotherapy. MOMP.Sig, offering a more precise indicator of tumour cell MOMP activity, demonstrated outstanding predictive efficacy in machine-learning models. Moreover, with the assistance of the MOMP.Sig model, FOXO1 was identified as a core modulator that promotes immune resistance. Finally, these findings were successfully validated in clinical immunotherapy cohorts of skin cutaneous melanoma and triple-negative breast cancer patients.

CONCLUSIONS

This study enhances our understanding of MOMP activity in immune modulation, providing valuable insights for more effective immunotherapeutic strategies across diverse tumours.

Inference of gene regulatory networks based on directed graph convolutional networks

Inferring gene regulatory network (GRN) is one of the important challenges in systems biology, and many outstanding computational methods have been proposed; however there remains some challenges especially in real datasets. In this study, we propose Directed Graph Convolutional neural network-based method for GRN inference (DGCGRN). To better understand and process the directed graph structure data of GRN, a directed graph convolutional neural network is conducted which retains the structural information of the directed graph while also making full use of neighbor node features. The local augmentation strategy is adopted in graph neural network to solve the problem of poor prediction accuracy caused by a large number of low-degree nodes in GRN. In addition, for real data such as E.coli, sequence features are obtained by extracting hidden features using Bi-GRU and calculating the statistical physicochemical characteristics of gene sequence. At the training stage, a dynamic update strategy is used to convert the obtained edge prediction scores into edge weights to guide the subsequent training process of the model. The results on synthetic benchmark datasets and real datasets show that the prediction performance of DGCGRN is significantly better than existing models. Furthermore, the case studies on bladder uroepithelial carcinoma and lung cancer cells also illustrate the performance of the proposed model.

Identification and Pharmacological Targeting of Treatment-Resistant, Stem-like Breast Cancer Cells for Combination Therapy

Tumors frequently harbor isogenic yet epigenetically distinct subpopulations of multi-potent cells with high tumor-initiating potential-often called Cancer Stem-Like Cells (CSLCs). These can display preferential resistance to standard-of-care chemotherapy. Single-cell analyses can help elucidate Master Regulator (MR) proteins responsible for governing the transcriptional state of these cells, thus revealing complementary dependencies that may be leveraged via combination therapy. Interrogation of single-cell RNA sequencing profiles from seven metastatic breast cancer patients, using perturbational profiles of clinically relevant drugs, identified drugs predicted to invert the activity of MR proteins governing the transcriptional state of chemoresistant CSLCs, which were then validated by CROP-seq assays. The top drug, the anthelmintic albendazole, depleted this subpopulation in vivo without noticeable cytotoxicity. Moreover, sequential cycles of albendazole and paclitaxel-a commonly used chemotherapeutic -displayed significant synergy in a patient-derived xenograft (PDX) from a TNBC patient, suggesting that network-based approaches can help develop mechanism-based combinatorial therapies targeting complementary subpopulations.

The Wnt-dependent master regulator NKX1-2 controls mouse pre-implantation development

Embryo size, specification, and homeostasis are regulated by a complex gene regulatory and signaling network. Here we used gene expression signatures of Wnt-activated mouse embryonic stem cell (mESC) clones to reverse engineer an mESC regulatory network. We identify NKX1-2 as a novel master regulator of preimplantation embryo development. We find that Nkx1-2 inhibition reduces nascent RNA synthesis, downregulates genes controlling ribosome biogenesis, RNA translation, and transport, and induces severe alteration of nucleolus structure, resulting in the exclusion of RNA polymerase I from nucleoli. In turn, NKX1-2 loss of function leads to chromosome missegregation in the 2- to 4-cell embryo stages, severe decrease in blastomere numbers, alterations of tight junctions (TJs), and impairment of microlumen coarsening. Overall, these changes impair the blastocoel expansion-collapse cycle and embryo cavitation, leading to altered lineage specification and developmental arrest.

Subtype-selective prenylated isoflavonoids disrupt regulatory drivers of MYCN-amplified cancers

Transcription factors have proven difficult to target with small molecules because they lack pockets necessary for potent binding. Disruption of protein expression can suppress targets and enable therapeutic intervention. To this end, we developed a drug discovery workflow that incorporates cell-line-selective screening and high-throughput expression profiling followed by regulatory network analysis to identify compounds that suppress regulatory drivers of disease. Applying this approach to neuroblastoma (NBL), we screened bioactive molecules in cell lines representing its MYC-dependent (MYCNA) and mesenchymal (MES) subtypes to identify selective compounds, followed by PLATESeq profiling of treated cells. This revealed compounds that disrupt a sub-network of MYCNA-specific regulatory proteins, resulting in MYCN degradation in vivo. The top hit was isopomiferin, a prenylated isoflavonoid that inhibited casein kinase 2 (CK2) in cells. Isopomiferin and its structural analogs inhibited MYC and MYCN in NBL and lung cancer cells, highlighting the general MYC-inhibiting potential of this unique scaffold.

Data-driven modeling of core gene regulatory network underlying leukemogenesis in IDH mutant AML

Acute myeloid leukemia (AML) is characterized by uncontrolled proliferation of poorly differentiated myeloid cells, with a heterogenous mutational landscape. Mutations in IDH1 and IDH2 are found in 20% of the AML cases. Although much effort has been made to identify genes associated with leukemogenesis, the regulatory mechanism of AML state transition is still not fully understood. To alleviate this issue, here we develop a new computational approach that integrates genomic data from diverse sources, including gene expression and ATAC-seq datasets, curated gene regulatory interaction databases, and mathematical modeling to establish models of context-specific core gene regulatory networks (GRNs) for a mechanistic understanding of tumorigenesis of AML with IDH mutations. The approach adopts a new optimization procedure to identify the top network according to its accuracy in capturing gene expression states and its flexibility to allow sufficient control of state transitions. From GRN modeling, we identify key regulators associated with the function of IDH mutations, such as DNA methyltransferase DNMT1, and network destabilizers, such as E2F1. The constructed core regulatory network and outcomes of in-silico network perturbations are supported by survival data from AML patients. We expect that the combined bioinformatics and systems-biology modeling approach will be generally applicable to elucidate the gene regulation of disease progression.

β-cell Jagged1 is sufficient but not necessary for islet Notch activity and insulin secretory defects in obese mice

OBJECTIVE

Notch signaling, re-activated in β cells from obese mice and causal to β cell dysfunction, is determined in part by transmembrane ligand availability in a neighboring cell. We hypothesized that β cell expression of Jagged1 determines the maladaptive Notch response and resultant insulin secretory defects in obese mice.

METHODS

We assessed expression of Notch pathway components in high-fat diet-fed (HFD) or leptin receptor-deficient (db/db) mice, and performed single-cell RNA sequencing (scRNA-Seq) in islets from patients with and without type 2 diabetes (T2D). We generated and performed glucose tolerance testing in inducible, β cell-specific Jagged1 gain-of- and loss-of-function mice. We also tested effects of monoclonal neutralizing antibodies to Jagged1 in glucose-stimulated insulin secretion (GSIS) assays in isolated islets.

RESULTS

Jag1 was the only Notch ligand that tracked with increased Notch activity in HFD-fed and db/db mice, as well as in metabolically-inflexible β cells enriched in patients with T2D. Neutralizing antibodies to block Jagged1 in islets isolated from HFD-fed and db/db mice potentiated GSIS ex vivo. To demonstrate if β cell Jagged1 is sufficient to cause glucose tolerance in vivo, we generated inducible β cell-specific Jag1 transgenic (β-Jag1TG) and loss-of-function (iβ-Jag1KO) mice. While forced Jagged1 impaired glucose intolerance due to reduced GSIS, loss of β cell Jagged1 did not protect against HFD-induced insulin secretory defects.

CONCLUSIONS

Jagged1 is increased in islets from obese mice and in patients with T2D, and neutralizing Jagged1 antibodies lead to improved GSIS, suggesting that inhibition of Jagged1-Notch signaling may have therapeutic benefit. However, genetic loss-of-function experiments suggest that β cells are not a likely source of the Jagged1 signal.

Tumor Explants Elucidate a Cascade of Paracrine SHH, WNT, and VEGF Signals Driving Pancreatic Cancer Angiosuppression

The sparse vascularity of pancreatic ductal adenocarcinoma (PDAC) presents a mystery: What prevents this aggressive malignancy from undergoing neoangiogenesis to counteract hypoxia and better support growth? An incidental finding from prior work on paracrine communication between malignant PDAC cells and fibroblasts revealed that inhibition of the Hedgehog (HH) pathway partially relieved angiosuppression, increasing tumor vascularity through unknown mechanisms. Initial efforts to study this phenotype were hindered by difficulties replicating the complex interactions of multiple cell types in vitro. Here we identify a cascade of paracrine signals between multiple cell types that act sequentially to suppress angiogenesis in PDAC. Malignant epithelial cells promote HH signaling in fibroblasts, leading to inhibition of noncanonical WNT signaling in fibroblasts and epithelial cells, thereby limiting VEGFR2-dependent activation of endothelial hypersprouting. This cascade was elucidated using human and murine PDAC explant models, which effectively retain the complex cellular interactions of native tumor tissues.

SIGNIFICANCE

We present a key mechanism of tumor angiosuppression, a process that sculpts the physiologic, cellular, and metabolic environment of PDAC. We further present a computational and experimental framework for the dissection of complex signaling cascades that propagate among multiple cell types in the tissue environment. This article is featured in Selected Articles from This Issue, p. 201.

The Prognostic Hub Gene POLE2 Promotes BLCA Cell Growth via the PI3K/AKT Signaling Pathway

BACKGROUND

BLCA is a common urothelial malignancy characterized by a high recurrence rate. Despite its prevalence, the molecular mechanisms underlying its development remain unclear.

AIMS

This study aimed to explore new prognostic biomarkers and investigate the underlying mechanism of bladder cancer (BLCA).

OBJECTIVE

The objective of this study is to identify key prognostic biomarkers for BLCA and to elucidate their roles in the disease.

METHODS

We first collected the overlapping DEGs from GSE42089 and TCGA-BLCA samples for the subsequent weighted gene co-expression network analysis (WGCNA) to find a key module. Then, key module genes were analyzed by the MCODE algorithm, prognostic risk model, expression and immunohistochemical staining to identify the prognostic hub gene. Finally, the hub gene was subjected to clinical feature analysis, as well as cellular function assays.

RESULTS

In WGCNA on 1037 overlapping genes, the blue module was the key module. After a series of bioinformatics analyses, POLE2 was identified as a prognostic hub gene in BLCA from potential genes (TROAP, POLE2, ANLN, and E2F8). POLE2 level was increased in BLCA and related to different clinical features of BLCA patients. Cellular assays showed that si-POLE2 inhibited BLCA proliferation, and si-POLE2+ 740Y-P in BLCA cells up-regulated the PI3K and AKT protein levels.

CONCLUSION

In conclusion, POLE2 was identified to be a promising prognostic biomarker as an oncogene in BLCA. It was also found that POLE2 exerts a promoting function by the PI3K/AKT signaling pathway in BLCA.

Regulus infers signed regulatory relations from few samples' information using discretization and likelihood constraints

MOTIVATION

Transcriptional regulation is performed by transcription factors (TF) binding to DNA in context-dependent regulatory regions and determines the activation or inhibition of gene expression. Current methods of transcriptional regulatory circuits inference, based on one or all of TF, regions and genes activity measurements require a large number of samples for ranking the candidate TF-gene regulation relations and rarely predict whether they are activations or inhibitions. We hypothesize that transcriptional regulatory circuits can be inferred from fewer samples by (1) fully integrating information on TF binding, gene expression and regulatory regions accessibility, (2) reducing data complexity and (3) using biology-based likelihood constraints to determine the global consistency between a candidate TF-gene relation and patterns of genes expressions and region activations, as well as qualify regulations as activations or inhibitions.

RESULTS

We introduce Regulus, a method which computes TF-gene relations from gene expressions, regulatory region activities and TF binding sites data, together with the genomic locations of all entities. After aggregating gene expressions and region activities into patterns, data are integrated into a RDF (Resource Description Framework) endpoint. A dedicated SPARQL (SPARQL Protocol and RDF Query Language) query retrieves all potential relations between expressed TF and genes involving active regulatory regions. These TF-region-gene relations are then filtered using biological likelihood constraints allowing to qualify them as activation or inhibition. Regulus provides signed relations consistent with public databases and, when applied to biological data, identifies both known and potential new regulators. Regulus is devoted to context-specific transcriptional circuits inference in human settings where samples are scarce and cell populations are closely related, using discretization into patterns and likelihood reasoning to decipher the most robust regulatory relations.

Plasticity of intragraft alloreactive T cell clones in human gut correlates with transplant outcomes

The site of transition between tissue-resident memory (TRM) and circulating phenotypes of T cells is unknown. We integrated clonotype, alloreactivity, and gene expression profiles of graft-repopulating recipient T cells in the intestinal mucosa at the single-cell level after human intestinal transplantation. Host-versus-graft (HvG)-reactive T cells were mainly distributed to TRM, effector T (Teff)/TRM, and T follicular helper compartments. RNA velocity analysis demonstrated a trajectory from TRM to Teff/TRM clusters in association with rejection. By integrating pre- and post-transplantation (Tx) mixed lymphocyte reaction-determined alloreactive repertoires, we observed that pre-existing HvG-reactive T cells that demonstrated tolerance in the circulation were dominated by TRM profiles in quiescent allografts. Putative de novo HvG-reactive clones showed a transcriptional profile skewed to cytotoxic effectors in rejecting grafts. Inferred protein regulon network analysis revealed upstream regulators that accounted for the effector and tolerant T cell states. We demonstrate Teff/TRM interchangeability for individual T cell clones with known (allo)recognition in the human gut, providing novel insight into TRM biology.

A review on gene regulatory network reconstruction algorithms based on single cell RNA sequencing

BACKGROUND

Understanding gene regulatory networks (GRNs) is essential for unraveling the molecular mechanisms governing cellular behavior. With the advent of high-throughput transcriptome measurement technology, researchers have aimed to reverse engineer the biological systems, extracting gene regulatory rules from their outputs, which represented by gene expression data. Bulk RNA sequencing, a widely used method for measuring gene expression, has been employed for GRN reconstruction. However, it falls short in capturing dynamic changes in gene expression at the level of individual cells since it averages gene expression across mixed cell populations.

OBJECTIVE

In this review, we provide an overview of 15 GRN reconstruction tools and discuss their respective strengths and limitations, particularly in the context of single cell RNA sequencing (scRNA-seq).

METHODS

Recent advancements in scRNA-seq break new ground of GRN reconstruction. They offer snapshots of the individual cell transcriptomes and capturing dynamic changes. We emphasize how these technological breakthroughs have enhanced GRN reconstruction.

CONCLUSION

GRN reconstructors can be classified based on their requirement for cellular trajectory, which represents a dynamical cellular process including differentiation, aging, or disease progression. Benchmarking studies support the superiority of GRN reconstructors that do not require trajectory analysis in identifying regulator-target relationships. However, methods equipped with trajectory analysis demonstrate better performance in identifying key regulatory factors. In conclusion, researchers should select a suitable GRN reconstructor based on their specific research objectives.

The Role of Transcription Factors in the Loss of Inter-Chromosomal Co-Expression for Breast Cancer Subtypes

Breast cancer encompasses a diverse array of subtypes, each exhibiting distinct clinical characteristics and treatment responses. Unraveling the underlying regulatory mechanisms that govern gene expression patterns in these subtypes is essential for advancing our understanding of breast cancer biology. Gene co-expression networks (GCNs) help us identify groups of genes that work in coordination. Previous research has revealed a marked reduction in the interaction of genes located on different chromosomes within GCNs for breast cancer, as well as for lung, kidney, and hematopoietic cancers. However, the reasons behind why genes on the same chromosome often co-express remain unclear. In this study, we investigate the role of transcription factors in shaping gene co-expression networks within the four main breast cancer subtypes: Luminal A, Luminal B, HER2+, and Basal, along with normal breast tissue. We identify communities within each GCN and calculate the transcription factors that may regulate these communities, comparing the results across different phenotypes. Our findings indicate that, in general, regulatory behavior is to a large extent similar among breast cancer molecular subtypes and even in healthy networks. This suggests that transcription factor motif usage does not fully determine long-range co-expression patterns. Specific transcription factor motifs, such as CCGGAAG, appear frequently across all phenotypes, even involving multiple highly connected transcription factors. Additionally, certain transcription factors exhibit unique actions in specific subtypes but with limited influence. Our research demonstrates that the loss of inter-chromosomal co-expression is not solely attributable to transcription factor regulation. Although the exact mechanism responsible for this phenomenon remains elusive, this work contributes to a better understanding of gene expression regulatory programs in breast cancer.

A multi-layered network model identifies Akt1 as a common modulator of neurodegeneration

The accumulation of misfolded and aggregated proteins is a hallmark of neurodegenerative proteinopathies. Although multiple genetic loci have been associated with specific neurodegenerative diseases (NDs), molecular mechanisms that may have a broader relevance for most or all proteinopathies remain poorly resolved. In this study, we developed a multi-layered network expansion (MLnet) model to predict protein modifiers that are common to a group of diseases and, therefore, may have broader pathophysiological relevance for that group. When applied to the four NDs Alzheimer's disease (AD), Huntington's disease, and spinocerebellar ataxia types 1 and 3, we predicted multiple members of the insulin pathway, including PDK1, Akt1, InR, and sgg (GSK-3β), as common modifiers. We validated these modifiers with the help of four Drosophila ND models. Further evaluation of Akt1 in human cell-based ND models revealed that activation of Akt1 signaling by the small molecule SC79 increased cell viability in all models. Moreover, treatment of AD model mice with SC79 enhanced their long-term memory and ameliorated dysregulated anxiety levels, which are commonly affected in AD patients. These findings validate MLnet as a valuable tool to uncover molecular pathways and proteins involved in the pathophysiology of entire disease groups and identify potential therapeutic targets that have relevance across disease boundaries. MLnet can be used for any group of diseases and is available as a web tool at http://ssbio.cau.ac.kr/software/mlnet.

DeepGAMI: deep biologically guided auxiliary learning for multimodal integration and imputation to improve genotype-phenotype prediction

BACKGROUND

Genotypes are strongly associated with disease phenotypes, particularly in brain disorders. However, the molecular and cellular mechanisms behind this association remain elusive. With emerging multimodal data for these mechanisms, machine learning methods can be applied for phenotype prediction at different scales, but due to the black-box nature of machine learning, integrating these modalities and interpreting biological mechanisms can be challenging. Additionally, the partial availability of these multimodal data presents a challenge in developing these predictive models.

METHOD

To address these challenges, we developed DeepGAMI, an interpretable neural network model to improve genotype-phenotype prediction from multimodal data. DeepGAMI leverages functional genomic information, such as eQTLs and gene regulation, to guide neural network connections. Additionally, it includes an auxiliary learning layer for cross-modal imputation allowing the imputation of latent features of missing modalities and thus predicting phenotypes from a single modality. Finally, DeepGAMI uses integrated gradient to prioritize multimodal features for various phenotypes.

RESULTS

We applied DeepGAMI to several multimodal datasets including genotype and bulk and cell-type gene expression data in brain diseases, and gene expression and electrophysiology data of mouse neuronal cells. Using cross-validation and independent validation, DeepGAMI outperformed existing methods for classifying disease types, and cellular and clinical phenotypes, even using single modalities (e.g., AUC score of 0.79 for Schizophrenia and 0.73 for cognitive impairment in Alzheimer's disease).

CONCLUSION

We demonstrated that DeepGAMI improves phenotype prediction and prioritizes phenotypic features and networks in multiple multimodal datasets in complex brains and brain diseases. Also, it prioritized disease-associated variants, genes, and regulatory networks linked to different phenotypes, providing novel insights into the interpretation of gene regulatory mechanisms. DeepGAMI is open-source and available for general use.

Demographic confounders distort inference of gene regulatory and gene co-expression networks in cancer

Gene regulatory networks (GRNs) and gene co-expression networks (GCNs) allow genome-wide exploration of molecular regulation patterns in health and disease. The standard approach for obtaining GRNs and GCNs is to infer them from gene expression data, using computational network inference methods. However, since network inference methods are usually applied on aggregate data, distortion of the networks by demographic confounders might remain undetected, especially because gene expression patterns are known to vary between different demographic groups. In this paper, we present a computational framework to systematically evaluate the influence of demographic confounders on network inference from gene expression data. Our framework compares similarities between networks inferred for different demographic groups with similarity distributions obtained for random splits of the expression data. Moreover, it allows to quantify to which extent demographic groups are represented by networks inferred from the aggregate data in a confounder-agnostic way. We apply our framework to test four widely used GRN and GCN inference methods as to their robustness w. r. t. confounding by age, ethnicity and sex in cancer. Our findings based on more than $ {44000}$ inferred networks indicate that age and sex confounders play an important role in network inference for certain cancer types, emphasizing the importance of incorporating an assessment of the effect of demographic confounders into network inference workflows. Our framework is available as a Python package on GitHub: https://github.com/bionetslab/grn-confounders.

Actionable loss of SLF2 drives B-cell lymphomagenesis and impairs the DNA damage response

The DNA damage response (DDR) acts as a barrier to malignant transformation and is often impaired during tumorigenesis. Exploiting the impaired DDR can be a promising therapeutic strategy; however, the mechanisms of inactivation and corresponding biomarkers are incompletely understood. Starting from an unbiased screening approach, we identified the SMC5-SMC6 Complex Localization Factor 2 (SLF2) as a regulator of the DDR and biomarker for a B-cell lymphoma (BCL) patient subgroup with an adverse prognosis. SLF2-deficiency leads to loss of DDR factors including Claspin (CLSPN) and consequently impairs CHK1 activation. In line with this mechanism, genetic deletion of Slf2 drives lymphomagenesis in vivo. Tumor cells lacking SLF2 are characterized by a high level of DNA damage, which leads to alterations of the post-translational SUMOylation pathway as a safeguard. The resulting co-dependency confers synthetic lethality to a clinically applicable SUMOylation inhibitor (SUMOi), and inhibitors of the DDR pathway act highly synergistic with SUMOi. Together, our results identify SLF2 as a DDR regulator and reveal co-targeting of the DDR and SUMOylation as a promising strategy for treating aggressive lymphoma.

Dictys: dynamic gene regulatory network dissects developmental continuum with single-cell multiomics

Gene regulatory networks (GRNs) are key determinants of cell function and identity and are dynamically rewired during development and disease. Despite decades of advancement, challenges remain in GRN inference, including dynamic rewiring, causal inference, feedback loop modeling and context specificity. To address these challenges, we develop Dictys, a dynamic GRN inference and analysis method that leverages multiomic single-cell assays of chromatin accessibility and gene expression, context-specific transcription factor footprinting, stochastic process network and efficient probabilistic modeling of single-cell RNA-sequencing read counts. Dictys improves GRN reconstruction accuracy and reproducibility and enables the inference and comparative analysis of context-specific and dynamic GRNs across developmental contexts. Dictys' network analyses recover unique insights in human blood and mouse skin development with cell-type-specific and dynamic GRNs. Its dynamic network visualizations enable time-resolved discovery and investigation of developmental driver transcription factors and their regulated targets. Dictys is available as a free, open-source and user-friendly Python package.

Blockade of CCR1 induces a phenotypic shift in macrophages and triggers a favorable antilymphoma activity

Tumor-associated macrophages (TAMs) within the tumor microenvironment (TME) play an important role in tumor growth and progression. TAMs have been involved in producing immunosuppressive TME via various factors; however, the underlying mechanisms remain unclear in B-cell lymphoma, including mantle cell lymphoma (MCL). We identified that chemokine receptor-1 (CCR1) is highly expressed on monocytes (Mo) and macrophages (MΦ), and CCR1 pharmacological inhibition or CCR1 siRNA abolished lymphoma-mediated Mo/MΦ migration in a chemotaxis assay. The deficiency of host CCR1 (CCR1 KO) was associated with decreased infiltration of peritoneal-MΦ compared with WT-CCR1. Functional studies indicated that the genetic depletion of CCR1 or treatment inhibited protumor MΦ (M2-like) phenotype by decreasing CD206 and IL-10 expression. Moreover, CCR1 depletion reprogrammed MΦ toward an MHCII+/TNFα+ immunogenic phenotype. Mechanistically, protumor MΦ driven-IL-10 provides a positive feedback loop to tumor-CCL3 by regulating the CCL3 promoter via STAT1 signaling. Therapeutic in vivo targeting of CCR1 with CCR1 antagonist BX-471 significantly reduced FC-muMCL1 mouse tumors in the syngeneic MCL model by the depletion of M2-TAMs and increased infiltration of cytotoxic CD8+ T cells. Our study established that CCR1 exerts a pivotal role in macrophage programming, thus shaping protumor TME and lymphoma progression. CCR1 inhibition through CCR1 antagonists may be a promising therapeutic strategy to reprogram macrophages in lymphoma-TME and achieve better clinical outcomes in patients.

Tff2 defines transit-amplifying pancreatic acinar progenitors that lack regenerative potential and are protective against Kras-driven carcinogenesis

While adult pancreatic stem cells are thought not to exist, it is now appreciated that the acinar compartment harbors progenitors, including tissue-repairing facultative progenitors (FPs). Here, we study a pancreatic acinar population marked by trefoil factor 2 (Tff2) expression. Long-term lineage tracing and single-cell RNA sequencing (scRNA-seq) analysis of Tff2-DTR-CreERT2-targeted cells defines a transit-amplifying progenitor (TAP) population that contributes to normal homeostasis. Following acute and chronic injury, Tff2+ cells, distinct from FPs, undergo depopulation but are eventually replenished. At baseline, oncogenic KrasG12D-targeted Tff2+ cells are resistant to PDAC initiation. However, KrasG12D activation in Tff2+ cells leads to survival and clonal expansion following pancreatitis and a cancer stem/progenitor cell-like state. Selective ablation of Tff2+ cells prior to KrasG12D activation in Mist1+ acinar or Dclk1+ FP cells results in enhanced tumorigenesis, which can be partially rescued by adenoviral Tff2 treatment. Together, Tff2 defines a pancreatic TAP population that protects against Kras-driven carcinogenesis.

Reversing pancreatic β-cell dedifferentiation in the treatment of type 2 diabetes

The maintenance of glucose homeostasis is fundamental for survival and health. Diabetes develops when glucose homeostasis fails. Type 2 diabetes (T2D) is characterized by insulin resistance and pancreatic β-cell failure. The failure of β-cells to compensate for insulin resistance results in hyperglycemia, which in turn drives altered lipid metabolism and β-cell failure. Thus, insulin secretion by pancreatic β-cells is a primary component of glucose homeostasis. Impaired β-cell function and reduced β-cell mass are found in diabetes. Both features stem from a failure to maintain β-cell identity, which causes β-cells to dedifferentiate into nonfunctional endocrine progenitor-like cells or to trans-differentiate into other endocrine cell types. In this regard, one of the key issues in achieving disease modification is how to reestablish β-cell identity. In this review, we focus on the causes and implications of β-cell failure, as well as its potential reversibility as a T2D treatment.

Emerging Insights into Liver X Receptor α in the Tumorigenesis and Therapeutics of Human Cancers

Liver X receptor α (LXRα), a member of the nuclear receptor superfamily, is identified as a protein activated by ligands that interacts with the promoters of specific genes. It regulates cholesterol, bile acid, and lipid metabolism in normal physiological processes, and it participates in the development of some related diseases. However, many studies have demonstrated that LXRα is also involved in regulating numerous human malignancies. Aberrant LXRα expression is emerging as a fundamental and pivotal factor in cancer cell proliferation, invasion, apoptosis, and metastasis. Herein, we outline the expression levels of LXRα between tumor tissues and normal tissues via the Oncomine and Tumor Immune Estimation Resource (TIMER) 2.0 databases; summarize emerging insights into the roles of LXRα in the development, progression, and treatment of different human cancers and their diversified mechanisms; and highlight that LXRα can be a biomarker and therapeutic target in diverse cancers.

An integrative bioinformatics investigation and experimental validation of chromobox family in diffuse large B-cell lymphoma

BACKGROUND

Diffuse large B-cell lymphoma (DLBCL) is one of the most aggressive malignant tumors. Chromobox (CBX) family plays the role of oncogenes in various malignancies.

METHODS

The transcriptional and protein levels of CBX family were confirmed by GEPIA, Oncomine, CCLE, and HPA database. Screening of co-expressed genes and gene function enrichment analysis were performed by GeneMANIA and DAVID 6.8. The prognostic value, immune cell infiltration and drug sensitivity analysis of CBX family in DLBCL were performed by Genomicscape, TIMER2.0, and GSCALite database. Confirmatory Tests of CBX family protein expression in DLBCL were performed by immunohistochemistry.

RESULTS

The mRNA and protein expressions of CBX1/2/3/5/6 were higher in DLBCL tissues than control groups. Enrichment analysis showed that the functions of CBX family were mainly related to chromatin remodeling, methylation-dependent protein binding, and VEGF signaling pathway. The high mRNA expressions of CBX2/3/5/6 were identified to be associated with short overall survival (OS) in DLBCL patients. Multivariate COX regression indicated that CBX3 was independent prognostic marker. Immune infiltration analysis revealed that the mRNA expressions of CBX family (especially CBX1, CBX5, and CBX6) in DLBCL were significantly correlated with the infiltration of most immune cells (including B cells, CD8 + T cells, CD4 + T cells, neutrophils, monocytes, macrophages, and Treg cells). Meanwhile, there was a strong correlation between the expression levels of CBX1/5/6 and surface markers of immune cells, such as the widely studied PVR-like protein receptor/ligand and PDL-1 immune checkpoint. Notably, our study found that DLBCL cells with CBX1 over-expression were resistant to the common anti-tumor drugs, but CBX2/5 had two polarities. Finally, we confirmed the higher expressions of CBX1/2/3/5/6 in DLBCL tissues compared with control groups by immunohistochemistry.

CONCLUSION

We provided a detailed analysis of the relationship between the CBX family and the prognosis of DLBCL. Distinguished from other studies, We found that high mRNA expressions of CBX2/3/5/6 were associated with poor prognosis in DLBCL patients, and Multivariate COX regression indicated that CBX3 was independent prognostic marker. Besides, our study also found an association between the CBX family and anti-tumour drug resistance, and provided a relationship between CBX family expression and immune cell infiltration.

A Transcriptome-Based Precision Oncology Platform for Patient-Therapy Alignment in a Diverse Set of Treatment-Resistant Malignancies

Predicting in vivo response to antineoplastics remains an elusive challenge. We performed a first-of-kind evaluation of two transcriptome-based precision cancer medicine methodologies to predict tumor sensitivity to a comprehensive repertoire of clinically relevant oncology drugs, whose mechanism of action we experimentally assessed in cognate cell lines. We enrolled patients with histologically distinct, poor-prognosis malignancies who had progressed on multiple therapies, and developed low-passage, patient-derived xenograft models that were used to validate 35 patient-specific drug predictions. Both OncoTarget, which identifies high-affinity inhibitors of individual master regulator (MR) proteins, and OncoTreat, which identifies drugs that invert the transcriptional activity of hyperconnected MR modules, produced highly significant 30-day disease control rates (68% and 91%, respectively). Moreover, of 18 OncoTreat-predicted drugs, 15 induced the predicted MR-module activity inversion in vivo. Predicted drugs significantly outperformed antineoplastic drugs selected as unpredicted controls, suggesting these methods may substantively complement existing precision cancer medicine approaches, as also illustrated by a case study.

SIGNIFICANCE

Complementary precision cancer medicine paradigms are needed to broaden the clinical benefit realized through genetic profiling and immunotherapy. In this first-in-class application, we introduce two transcriptome-based tumor-agnostic systems biology tools to predict drug response in vivo. OncoTarget and OncoTreat are scalable for the design of basket and umbrella clinical trials. This article is highlighted in the In This Issue feature, p. 1275.

Systematic elucidation and pharmacological targeting of tumor-infiltrating regulatory T cell master regulators

Due to their immunosuppressive role, tumor-infiltrating regulatory T cells (TI-Tregs) represent attractive immuno-oncology targets. Analysis of TI vs. peripheral Tregs (P-Tregs) from 36 patients, across four malignancies, identified 17 candidate master regulators (MRs) as mechanistic determinants of TI-Treg transcriptional state. Pooled CRISPR-Cas9 screening in vivo, using a chimeric hematopoietic stem cell transplant model, confirmed the essentiality of eight MRs in TI-Treg recruitment and/or retention without affecting other T cell subtypes, and targeting one of the most significant MRs (Trps1) by CRISPR KO significantly reduced ectopic tumor growth. Analysis of drugs capable of inverting TI-Treg MR activity identified low-dose gemcitabine as the top prediction. Indeed, gemcitabine treatment inhibited tumor growth in immunocompetent but not immunocompromised allografts, increased anti-PD-1 efficacy, and depleted MR-expressing TI-Tregs in vivo. This study provides key insight into Treg signaling, specifically in the context of cancer, and a generalizable strategy to systematically elucidate and target MR proteins in immunosuppressive subpopulations.

OX40 shapes an inflamed tumor immune microenvironment and predicts response to immunochemotherapy in diffuse large B-cell lymphoma

OX40 enhances the T-cell activation via costimulatory signaling. However, its molecular characteristics and value in predicting response to immunochemotherapy in DLBCL remain largely unexplored. Here, we performed an integrative analysis of sequencing and multiplex immunofluorescence staining, and discovered abnormally higher expression of OX40 in DLBCL patients. Elevated OX40 could activate T cells leading to a higher immune score for tumor immune microenvironment (TiME). OX40 upregulation simultaneously happened with immune-related genes including PD-1, CTLA4 and TIGIT et,al. Patients with high OX40 expression exhibited a lower Ann Arbor stage and IPI score and more easily achieved a complete response/partial response. The analysis of infiltrated T-cell subset revealed that patients with a greater number of CD4+/OX40+ or CD8+/OX40+ T cells had a longer OS. Our findings indicated that OX40 shapes an inflamed tumor immune microenvironment and predicts response to immunochemotherapy, providing insights for the application of OX40 agonist in DLBCL patients.

Systematic Pan-cancer Functional Inference and Validation of Hyper, Hypo and Neomorphic Mutations

While the functional effects of many recurrent cancer mutations have been characterized, the TCGA repository comprises more than 10M non-recurrent events, whose function is unknown. We propose that the context specific activity of transcription factor (TF) proteins-as measured by expression of their transcriptional targets-provides a sensitive and accurate reporter assay to assess the functional role of oncoprotein mutations. Analysis of differentially active TFs in samples harboring mutations of unknown significance-compared to established gain (GOF/hypermorph) or loss (LOF/hypomorph) of function-helped functionally characterize 577,866 individual mutational events across TCGA cohorts, including identification of mutations that are either neomorphic (gain of novel function) or phenocopy other mutations ( mutational mimicry ). Validation using mutation knock-in assays confirmed 15 out of 15 predicted gain and loss of function mutations and 15 of 20 predicted neomorphic mutations. This could help determine targeted therapy in patients with mutations of unknown significance in established oncoproteins.

WMDS.net: a network control framework for identifying key players in transcriptome programs

MOTIVATION

Mammalian cells can be transcriptionally reprogramed to other cellular phenotypes. Controllability of such complex transitions in transcriptional networks underlying cellular phenotypes is an inherent biological characteristic. This network controllability can be interpreted by operating a few key regulators to guide the transcriptional program from one state to another. Finding the key regulators in the transcriptional program can provide key insights into the network state transition underlying cellular phenotypes.

RESULTS

To address this challenge, here, we proposed to identify the key regulators in the transcriptional co-expression network as a minimum dominating set (MDS) of driver nodes that can fully control the network state transition. Based on the theory of structural controllability, we developed a weighted MDS network model (WMDS.net) to find the driver nodes of differential gene co-expression networks. The weight of WMDS.net integrates the degree of nodes in the network and the significance of gene co-expression difference between two physiological states into the measurement of node controllability of the transcriptional network. To confirm its validity, we applied WMDS.net to the discovery of cancer driver genes in RNA-seq datasets from The Cancer Genome Atlas. WMDS.net is powerful among various cancer datasets and outperformed the other top-tier tools with a better balance between precision and recall.

AVAILABILITY AND IMPLEMENTATION

https://github.com/chaofen123/WMDS.net.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

OncoLoop: A Network-Based Precision Cancer Medicine Framework

Prioritizing treatments for individual patients with cancer remains challenging, and performing coclinical studies using patient-derived models in real time is often unfeasible. To circumvent these challenges, we introduce OncoLoop, a precision medicine framework that predicts drug sensitivity in human tumors and their preexisting high-fidelity (cognate) model(s) by leveraging drug perturbation profiles. As a proof of concept, we applied OncoLoop to prostate cancer using genetically engineered mouse models (GEMM) that recapitulate a broad spectrum of disease states, including castration-resistant, metastatic, and neuroendocrine prostate cancer. Interrogation of human prostate cancer cohorts by Master Regulator (MR) conservation analysis revealed that most patients with advanced prostate cancer were represented by at least one cognate GEMM-derived tumor (GEMM-DT). Drugs predicted to invert MR activity in patients and their cognate GEMM-DTs were successfully validated in allograft, syngeneic, and patient-derived xenograft (PDX) models of tumors and metastasis. Furthermore, OncoLoop-predicted drugs enhanced the efficacy of clinically relevant drugs, namely, the PD-1 inhibitor nivolumab and the AR inhibitor enzalutamide.

SIGNIFICANCE

OncoLoop is a transcriptomic-based experimental and computational framework that can support rapid-turnaround coclinical studies to identify and validate drugs for individual patients, which can then be readily adapted to clinical practice. This framework should be applicable in many cancer contexts for which appropriate models and drug perturbation data are available. This article is highlighted in the In This Issue feature, p. 247.

Identifying key multifunctional components shared by critical cancer and normal liver pathways via SparseGMM

Despite the abundance of multimodal data, suitable statistical models that can improve our understanding of diseases with genetic underpinnings are challenging to develop. Here, we present SparseGMM, a statistical approach for gene regulatory network discovery. SparseGMM uses latent variable modeling with sparsity constraints to learn Gaussian mixtures from multiomic data. By combining coexpression patterns with a Bayesian framework, SparseGMM quantitatively measures confidence in regulators and uncertainty in target gene assignment by computing gene entropy. We apply SparseGMM to liver cancer and normal liver tissue data and evaluate discovered gene modules in an independent single-cell RNA sequencing (scRNA-seq) dataset. SparseGMM identifies PROCR as a regulator of angiogenesis and PDCD1LG2 and HNF4A as regulators of immune response and blood coagulation in cancer. Furthermore, we show that more genes have significantly higher entropy in cancer compared with normal liver. Among high-entropy genes are key multifunctional components shared by critical pathways, including p53 and estrogen signaling.

14-3-3 and Smad2/3 are crucial mediators of atypical-PKCs: Implications for neuroblastoma progression

Neuroblastoma (NB) is a cancer that develops in the neuroblasts. It is the most common cancer in children under the age of 1 year, accounting for approximately 6% of all cancers. The prognosis of NB is linked to both age and degree of cell differentiation. This results in a range of survival rates for patients, with outcomes ranging from recurrence and mortality to high survival rates and tumor regression. Our previous work indicated that PKC-ι promotes cell proliferation in NB cells through the PKC-ι/Cdk7/Cdk2 cascade. We report on two atypical protein kinase inhibitors as potential therapeutic candidates against BE(2)-C and BE(2)-M17 cells: a PKC-ι-specific 5-amino-1-2,3-dihydroxy-4-(methylcyclopentyl)-1H-imidazole-4-carboxamide and a PKC-ζ specific 8-hydroxy-1,3,6-naphthalenetrisulfonic acid. Both compounds induced apoptosis and retarded the epithelial-mesenchymal transition (EMT) of NB cells. Proteins 14-3-3 and Smad2/3 acted as central regulators of aPKC-driven progression in BE(2)-C and BE(2)-M17 cells in relation to the Akt1/NF-κB and TGF-β pathways. Data indicates that aPKCs upregulate Akt1/NF-κB and TGF-β pathways in NB cells through an association with 14-3-3 and Smad2/3 that can be diminished by aPKC inhibitors. In summary, both inhibitors appear to be promising potential neuroblastoma therapeutics and merit further research.

First Report of FARSA in the Regulation of Cell Cycle and Survival in Mantle Cell Lymphoma Cells via PI3K-AKT and FOXO1-RAG1 Axes

Cancer-associated factors have been largely identified in the understanding of tumorigenesis and progression. However, aminoacyl-transfer RNA (tRNA) synthetases (aaRSs) have so far been neglected in cancer research due to their canonical activities in protein translation and synthesis. FARSA, the alpha subunit of the phenylalanyl-tRNA synthetase is elevated across many cancer types, but its function in mantle cell lymphoma (MCL) remains undetermined. Herein, we found the lowest levels of FARSA in patients with MCL compared with other subtypes of lymphomas, and the same lower levels of FARSA were observed in chemoresistant MCL cell lines. Unexpectedly, despite the essential catalytic roles of FARSA, knockdown of FARSA in MCL cells did not lead to cell death but resulted in accelerated cell proliferation and cell cycle, whereas overexpression of FARSA induced remarkable cell-cycle arrest and overwhelming apoptosis. Further RNA sequencing (RNA-seq) analysis and validation experiments confirmed a strong connection between FARSA and cell cycle in MCL cells. Importantly, FARSA leads to the alteration of cell cycle and survival via both PI3K-AKT and FOXO1-RAG1 axes, highlighting a FARSA-mediated regulatory network in MCL cells. Our findings, for the first time, reveal the noncanonical roles of FARSA in MCL cells, and provide novel insights into understanding the pathogenesis and progression of B-cell malignancies.

A pseudo-value regression approach for differential network analysis of co-expression data

BACKGROUND

The differential network (DN) analysis identifies changes in measures of association among genes under two or more experimental conditions. In this article, we introduce a pseudo-value regression approach for network analysis (PRANA). This is a novel method of differential network analysis that also adjusts for additional clinical covariates. We start from mutual information criteria, followed by pseudo-value calculations, which are then entered into a robust regression model.

RESULTS

This article assesses the model performances of PRANA in a multivariable setting, followed by a comparison to dnapath and DINGO in both univariable and multivariable settings through variety of simulations. Performance in terms of precision, recall, and F1 score of differentially connected (DC) genes is assessed. By and large, PRANA outperformed dnapath and DINGO, neither of which is equipped to adjust for available covariates such as patient-age. Lastly, we employ PRANA in a real data application from the Gene Expression Omnibus database to identify DC genes that are associated with chronic obstructive pulmonary disease to demonstrate its utility.

CONCLUSION

To the best of our knowledge, this is the first attempt of utilizing a regression modeling for DN analysis by collective gene expression levels between two or more groups with the inclusion of additional clinical covariates. By and large, adjusting for available covariates improves accuracy of a DN analysis.

A polyamine-centric, blood-based metabolite panel predictive of poor response to CAR-T cell therapy in large B cell lymphoma

Anti-CD19 chimeric antigen receptor (CAR) T cell therapy for relapsed or refractory (r/r) large B cell lymphoma (LBCL) results in durable response in only a subset of patients. MYC overexpression in LBCL tumors is associated with poor response to treatment. We tested whether an MYC-driven polyamine signature, as a liquid biopsy, is predictive of response to anti-CD19 CAR-T therapy in patients with r/r LBCL. Elevated plasma acetylated polyamines were associated with non-durable response. Concordantly, increased expression of spermidine synthase, a key enzyme that regulates levels of acetylated spermidine, was prognostic for survival in r/r LBCL. A broad metabolite screen identified additional markers that resulted in a 6-marker panel (6MetP) consisting of acetylspermidine, diacetylspermidine, and lysophospholipids, which was validated in an independent set from another institution as predictive of non-durable response to CAR-T therapy. A polyamine centric metabolomics liquid biopsy panel has predictive value for response to CAR-T therapy in r/r LBCL.

An integrated in silico-in vitro approach for identifying therapeutic targets against osteoarthritis

BACKGROUND

Without the availability of disease-modifying drugs, there is an unmet therapeutic need for osteoarthritic patients. During osteoarthritis, the homeostasis of articular chondrocytes is dysregulated and a phenotypical transition called hypertrophy occurs, leading to cartilage degeneration. Targeting this phenotypic transition has emerged as a potential therapeutic strategy. Chondrocyte phenotype maintenance and switch are controlled by an intricate network of intracellular factors, each influenced by a myriad of feedback mechanisms, making it challenging to intuitively predict treatment outcomes, while in silico modeling can help unravel that complexity. In this study, we aim to develop a virtual articular chondrocyte to guide experiments in order to rationalize the identification of potential drug targets via screening of combination therapies through computational modeling and simulations.

RESULTS

We developed a signal transduction network model using knowledge-based and data-driven (machine learning) modeling technologies. The in silico high-throughput screening of (pairwise) perturbations operated with that network model highlighted conditions potentially affecting the hypertrophic switch. A selection of promising combinations was further tested in a murine cell line and primary human chondrocytes, which notably highlighted a previously unreported synergistic effect between the protein kinase A and the fibroblast growth factor receptor 1.

CONCLUSIONS

Here, we provide a virtual articular chondrocyte in the form of a signal transduction interactive knowledge base and of an executable computational model. Our in silico-in vitro strategy opens new routes for developing osteoarthritis targeting therapies by refining the early stages of drug target discovery.

Correlation scan: identifying genomic regions that affect genetic correlations applied to fertility traits

Although the genetic correlations between complex traits have been estimated for more than a century, only recently we have started to map and understand the precise localization of the genomic region(s) that underpin these correlations. Reproductive traits are often genetically correlated. Yet, we don't fully understand the complexities, synergism, or trade-offs between male and female fertility. In this study, we used reproductive traits in two cattle populations (Brahman; BB, Tropical Composite; TC) to develop a novel framework termed correlation scan (CS). This framework was used to identify local regions associated with the genetic correlations between male and female fertility traits. Animals were genotyped with bovine high-density single nucleotide polymorphisms (SNPs) chip assay. The data used consisted of ~1000 individual records measured through frequent ovarian scanning for age at first corpus luteum (AGECL) and a laboratory assay for serum levels of insulin growth hormone (IGF1 measured in bulls, IGF1b, or cows, IGF1c). The methodology developed herein used correlations of 500-SNP effects in a 100-SNPs sliding window in each chromosome to identify local genomic regions that either drive or antagonize the genetic correlations between traits. We used Fisher's Z-statistics through a permutation method to confirm which regions of the genome harboured significant correlations. About 30% of the total genomic regions were identified as driving and antagonizing genetic correlations between male and female fertility traits in the two populations. These regions confirmed the polygenic nature of the traits being studied and pointed to genes of interest. For BB, the most important chromosome in terms of local regions is often located on bovine chromosome (BTA) 14. However, the important regions are spread across few different BTA's in TC. Quantitative trait loci (QTLs) and functional enrichment analysis revealed many significant windows co-localized with known QTLs related to milk production and fertility traits, especially puberty. In general, the enriched reproductive QTLs driving the genetic correlations between male and female fertility are the same for both cattle populations, while the antagonizing regions were population specific. Moreover, most of the antagonizing regions were mapped to chromosome X. These results suggest regions of chromosome X for further investigation into the trade-offs between male and female fertility. We compared the CS with two other recently proposed methods that map local genomic correlations. Some genomic regions were significant across methods. Yet, many significant regions identified with the CS were overlooked by other methods.