ARACNe-AP: gene network reverse engineering through adaptive partitioning inference of mutual information

Predicting transcription factor activity using prior biological information

Dysregulation of normal transcription factor activity is a common driver of disease. Therefore, the detection of aberrant transcription factor activity is important to understand disease pathogenesis. We have developed Priori, a method to predict transcription factor activity from RNA sequencing data. Priori has two key advantages over existing methods. First, Priori utilizes literature-supported regulatory information to identify transcription factor-target gene relationships. It then applies linear models to determine the impact of transcription factor regulation on the expression of its target genes. Second, results from a third-party benchmarking pipeline reveals that Priori detects aberrant activity from 124 single-gene perturbation experiments with higher sensitivity and specificity than 11 other methods. We applied Priori and other top-performing methods to predict transcription factor activity from two large primary patient datasets. Our work demonstrates that Priori uniquely discovered significant determinants of survival in breast cancer and identified mediators of drug response in leukemia.

Leucine rich repeat LGI family member 3: Integrative analyses support its prognostic association with pancreatic adenocarcinoma

Leucine rich repeat LGI family member 3 (LGI3) is a member of the LGI protein family. Previous studies of our group have reported that LGI3 is expressed in adipose tissue, skin and brain, and serves as a multifunctional cytokine. LGI3 may also be involved in cytokine networks in various cancers. This study aimed to analyze differentially expressed genes in pancreatic adenocarcinoma (PAC) tissues and PAC cohort data in order to evaluate the prognostic role of LGI3. The expression microarray and the PAC cohort data were analyzed by bioinformatic methods for differential expression, protein-protein interactions, functional enrichment and pathway analyses, gene co-expression network analysis, and prognostic association analysis. Results showed that LGI3 expression was significantly reduced in PAC tissues. Nineteen upregulated genes and 31 downregulated genes in PAC tissues were identified as LGI3-regulated genes. Protein-protein interaction network analysis demonstrated that 92% (46/50) of the LGI3-regulated genes that were altered in PACs belonged to a protein-protein interaction network cluster. Functional enrichment and gene co-expression network analyses demonstrated that these genes in the network cluster were associated with various processes including inflammatory and immune responses, metabolic processes, cell differentiation, and angiogenesis. PAC cohort analyses revealed that low expression levels of LGI3 were significantly associated with poor PAC prognosis. Analysis of favorable or unfavorable prognostic gene products in PAC showed that 93 LGI3-regulated genes were differentially associated with PAC prognosis. LGI3 expression was correlated with the tumor-infiltration levels of various immune cells. Taken together, these results suggested that LGI3 may be a potential prognostic marker of PAC.

The cycling and aging mouse female reproductive tract at single-cell resolution

The female reproductive tract (FRT) undergoes extensive remodeling during reproductive cycling. This recurrent remodeling and how it shapes organ-specific aging remains poorly explored. Using single-cell and spatial transcriptomics, we systematically characterized morphological and gene expression changes occurring in ovary, oviduct, uterus, cervix, and vagina at each phase of the mouse estrous cycle, during decidualization, and into aging. These analyses reveal that fibroblasts play central-and highly organ-specific-roles in FRT remodeling by orchestrating extracellular matrix (ECM) reorganization and inflammation. Our results suggest a model wherein recurrent FRT remodeling over reproductive lifespan drives the gradual, age-related development of fibrosis and chronic inflammation. This hypothesis was directly tested using chemical ablation of cycling, which reduced fibrotic accumulation during aging. Our atlas provides extensive detail into how estrus, pregnancy, and aging shape the organs of the female reproductive tract and reveals the unexpected cost of the recurrent remodeling required for reproduction.

SMARCAL1 is a dual regulator of innate immune signaling and PD-L1 expression that promotes tumor immune evasion

Genomic instability can trigger cancer-intrinsic innate immune responses that promote tumor rejection. However, cancer cells often evade these responses by overexpressing immune checkpoint regulators, such as PD-L1. Here, we identify the SNF2-family DNA translocase SMARCAL1 as a factor that favors tumor immune evasion by a dual mechanism involving both the suppression of innate immune signaling and the induction of PD-L1-mediated immune checkpoint responses. Mechanistically, SMARCAL1 limits endogenous DNA damage, thereby suppressing cGAS-STING-dependent signaling during cancer cell growth. Simultaneously, it cooperates with the AP-1 family member JUN to maintain chromatin accessibility at a PD-L1 transcriptional regulatory element, thereby promoting PD-L1 expression in cancer cells. SMARCAL1 loss hinders the ability of tumor cells to induce PD-L1 in response to genomic instability, enhances anti-tumor immune responses and sensitizes tumors to immune checkpoint blockade in a mouse melanoma model. Collectively, these studies uncover SMARCAL1 as a promising target for cancer immunotherapy.

Amphiregulin from regulatory T cells promotes liver fibrosis and insulin resistance in non-alcoholic steatohepatitis

Production of amphiregulin (Areg) by regulatory T (Treg) cells promotes repair after acute tissue injury. Here, we examined the function of Treg cells in non-alcoholic steatohepatitis (NASH), a setting of chronic liver injury. Areg-producing Treg cells were enriched in the livers of mice and humans with NASH. Deletion of Areg in Treg cells, but not in myeloid cells, reduced NASH-induced liver fibrosis. Chronic liver damage induced transcriptional changes associated with Treg cell activation. Mechanistically, Treg cell-derived Areg activated pro-fibrotic transcriptional programs in hepatic stellate cells via epidermal growth factor receptor (EGFR) signaling. Deletion of Areg in Treg cells protected mice from NASH-dependent glucose intolerance, which also was dependent on EGFR signaling on hepatic stellate cells. Areg from Treg cells promoted hepatocyte gluconeogenesis through hepatocyte detection of hepatic stellate cell-derived interleukin-6. Our findings reveal a maladaptive role for Treg cell-mediated tissue repair functions in chronic liver disease and link liver damage to NASH-dependent glucose intolerance.

Identification of a gene network driving the attenuated response to lipopolysaccharide of monocytes from hypertensive coronary artery disease patients

Introduction

The impact of cardiovascular disease (CVD) risk factors, encompassing various biological determinants and unhealthy lifestyles, on the functional dynamics of circulating monocytes-a pivotal cell type in CVD pathophysiology remains elusive. In this study, we aimed to elucidate the influence of CVD risk factors on monocyte transcriptional responses to an infectious stimulus.

Methods

We conducted a comparative analysis of monocyte gene expression profiles from the CTMM - CIRCULATING CELLS Cohort of coronary artery disease (CAD) patients, at baseline and after lipopolysaccharide (LPS) stimulation. Gene co-expression analysis was used to identify gene modules and their correlations with CVD risk factors, while pivotal transcription factors controlling the hub genes in these modules were identified by regulatory network analyses. The identified gene module was subjected to a drug repurposing screen, utilizing the LINCS L1000 database.

Results

Monocyte responsiveness to LPS showed a highly significant, negative correlation with blood pressure levels (ρ< -0.4; P<10-80). We identified a ZNF12/ZBTB43-driven gene module closely linked to diastolic blood pressure, suggesting that monocyte responses to infectious stimuli, such as LPS, are attenuated in CAD patients with elevated diastolic blood pressure. This attenuation appears associated with a dampening of the LPS-induced suppression of oxidative phosphorylation. Finally, we identified the serine-threonine inhibitor MW-STK33-97 as a drug candidate capable of reversing this aberrant LPS response.

Conclusions

Monocyte responses to infectious stimuli may be hampered in CAD patients with high diastolic blood pressure and this attenuated inflammatory response may be reversed by the serine-threonine inhibitor MW-STK33-97. Whether the identified gene module is a mere indicator of, or causal factor in diastolic blood pressure and the associated dampened LPS responses remains to be determined.

Identifying interactions among air pollutant emissions on diabetes prevalence in Northeast China using a complex network

BACKGROUND

Low air quality related to ambient air pollution is the largest environmental risk to health worldwide. Interactions between air pollution emissions may affect associations between air pollution exposure and chronic diseases. Therefore, this study aimed to quantify interactions among air pollution emissions and assess their effects on the association between air pollution and diabetes.

METHODS

After constructing long-term emission networks for six air pollutants based on data collected from routine monitoring stations in Northeast China, a mutual information network was used to quantify interactions among air pollution emissions. Multiple linear regression analysis was then used to explore the influence of emission interactions on the association between air pollution exposure and the prevalence of diabetes based on data reported from the Northeast Natural Cohort Study in China.

RESULTS

Complex network analysis detected three major emission sources in Northeast China located in Shenyang and Changchun. The effects of particulate matter (PM2.5 and PM10) and ground-level ozone (O3) emissions were limited to certain communities but could spread to other communities through emissions in Inner Mongolia. Emissions of sulfur dioxide (SO2), nitrogen dioxide (NO2), and carbon monoxide (CO) significantly influenced other communities. These results indicated that air pollutants in different geographic areas can interact directly or indirectly. Adjusting for interactions between emissions changed associations between air pollution emissions and diabetes prevalence, especially for PM2.5, NO2, and CO.

CONCLUSIONS

Complex network analysis is suitable for quantifying interactions among air pollution emissions and suggests that the effects of PM2.5 and NO2 emissions on health outcomes may have been overestimated in previous population studies while those of CO may have been underestimated. Further studies examining associations between air pollution and chronic diseases should consider controlling for the effects of interactions among pollution emissions.

From genomics to metabolomics: Deciphering sanguinarine biosynthesis in Dicranostigma leptopodum

Dicranostigma leptopodum (Maxim) Fedde (DLF) is a renowned medicinal plant in China, known to be rich in alkaloids. However, the unavailability of a reference genome has impeded investigation into its plant metabolism and genetic breeding potential. Here we present a high-quality chromosomal-level genome assembly for DLF, derived using a combination of Nanopore long-read sequencing, Illumina short-read sequencing and Hi-C technologies. Our assembly genome spans a size of 621.81 Mb with an impressive contig N50 of 93.04 Mb. We show that the species-specific whole-genome duplication (WGD) of DLF and Papaver somniferum corresponded to two rounds of WGDs of Papaver setigerum. Furthermore, we integrated comprehensive homology searching, gene family analyses and construction of a gene-to-metabolite network. These efforts led to the discovery of co-expressed transcription factors, including NAC and bZIP, alongside sanguinarine (SAN) pathway genes CYP719 (CFS and SPS). Notably, we identified P6H as a promising gene for enhancing SAN production. By providing the first reference genome for Dicranostigma, our study confirms the genomic underpinning of SAN biosynthesis and establishes a foundation for advancing functional genomic research on Papaveraceae species. Our findings underscore the pivotal role of high-quality genome assemblies in elucidating genetic variations underlying the evolutionary origin of secondary metabolites.

ADGRE5-centered Tsurv model in T cells recognizes responders to neoadjuvant cancer immunotherapy

Background

Neoadjuvant immunotherapy with anti-programmed death-1 (neo-antiPD1) has revolutionized perioperative methods for improvement of overall survival (OS), while approaches for major pathologic response patients' (MPR) recognition along with methods for overcoming non-MPR resistance are still in urgent need.

Methods

We utilized and integrated publicly-available immune checkpoint inhibitors regimens (ICIs) single-cell (sc) data as the discovery datasets, and innovatively developed a cell-communication analysis pipeline, along with a VIPER-based-SCENIC process, to thoroughly dissect MPR-responding subsets. Besides, we further employed our own non-small cell lung cancer (NSCLC) ICIs cohort's sc data for validation in-silico. Afterward, we resorted to ICIs-resistant murine models developed by us with multimodal investigation, including bulk-RNA-sequencing, Chip-sequencing and high-dimensional cytometry by time of flight (CYTOF) to consolidate our findings in-vivo. To comprehensively explore mechanisms, we adopted 3D ex-vivo hydrogel models for analysis. Furthermore, we constructed an ADGRE5-centered Tsurv model from our discovery dataset by machine learning (ML) algorithms for a wide range of tumor types (NSCLC, melanoma, urothelial cancer, etc.) and verified it in peripheral blood mononuclear cells (PBMCs) sc datasets.

Results

Through a meta-analysis of multimodal sequential sc sequencing data from pre-ICIs and post-ICIs, we identified an MPR-expanding T cells meta-cluster (MPR-E) in the tumor microenvironment (TME), characterized by a stem-like CD8+ T cluster (survT) with STAT5-ADGRE5 axis enhancement compared to non-MPR or pre-ICIs TME. Through multi-omics analysis of murine TME, we further confirmed the existence of survT with silenced function and immune checkpoints (ICs) in MPR-E. After verification of the STAT5-ADGRE5 axis of survT in independent ICIs cohorts, an ADGRE5-centered Tsurv model was then developed through ML for identification of MPR patients pre-ICIs and post-ICIs, both in TME and PBMCs, which was further verified in pan-cancer immunotherapy cohorts. Mechanistically, we unveiled ICIs stimulated ADGRE5 upregulation in a STAT5-IL32 dependent manner in a 3D ex-vivo system (3D-HYGTIC) developed by us previously, which marked Tsurv with better survival flexibility, enhanced stemness and potential cytotoxicity within TME.

Conclusion

Our research provides insights into mechanisms underlying MPR in neo-antiPD1 and a well-performed model for the identification of non-MPR.

BIC-LP: A Hybrid Higher-Order Dynamic Bayesian Network Score Function for Gene Regulatory Network Reconstruction

Reconstructing gene regulatory networks(GRNs) is an increasingly hot topic in bioinformatics. Dynamic Bayesian network(DBN) is a stochastic graph model commonly used as a vital model for GRN reconstruction. But probabilistic characteristics of biological networks and the existence of data noise bring great challenges to GRN reconstruction and always lead to many false positive/negative edges. ScoreLasso is a hybrid DBN score function combining DBN and linear regression with good performance. Its performance is, however, limited by first-order assumption and ignorance of the initial network of DBN. In this article, an integrated model based on higher-order DBN model, higher-order Lasso linear regression model and Pearson correlation model is proposed. Based on this, a hybrid higher-order DBN score function for GRN reconstruction is proposed, namely BIC-LP. BIC-LP score function is constructed by adding terms based on Lasso linear regression coefficients and Pearson correlation coefficients on classical BIC score function. Therefore, it could capture more information from dataset and curb information loss, compared with both many existing Bayesian family score functions and many state-of-the-art methods for GRN reconstruction. Experimental results show that BIC-LP can reasonably eliminate some false positive edges while retaining most true positive edges, so as to achieve better GRN reconstruction performance.

Systematic single cell RNA sequencing analysis reveals unique transcriptional regulatory networks of Atoh1-mediated hair cell conversion in adult mouse cochleae

Regeneration of mammalian cochlear hair cells (HCs) by modulating molecular pathways or transcription factors is a promising approach to hearing restoration; however, immaturity of the regenerated HCs in vivo remains a major challenge. Here, we analyzed a single cell RNA sequencing (scRNA-seq) dataset during Atoh1-induced supporting cell (SC) to hair cell (HC) conversion in adult mouse cochleae (Yamashita et al. (2018)) using multiple high-throughput sequencing analytical tools (WGCNA, SCENIC, ARACNE, and VIPER). Instead of focusing on differentially expressed genes, we established independent expression modules and confirmed the existence of multiple conversion stages. Gene regulatory network (GRN) analysis uncovered previously unidentified key regulators, including Nhlh1, Lhx3, Barhl1 and Nfia, that guide converted HC differentiation. Comparison of the late-stage converted HCs with the scRNA-seq data from neonatal mouse cochleae (Kolla et al. (2020)) revealed that they closely resemble postnatal day 1 wild-type OHCs, in contrast to other developmental stages. Using ARACNE and VIPER, we discovered multiple key regulators likely to promote conversion to a more mature OHC-like state, including Zbtb20, Nfia, Zmiz1, Gm14418, Bhlhe40, Six2, Fosb and Klf9. Our findings provide insights into the regulation of HC regeneration in adult mammalian cochleae in vivo and demonstrate an approach for analyzing GRNs in large scRNA-seq datasets.

Temporal evolution reveals bifurcated lineages in aggressive neuroendocrine small cell prostate cancer trans-differentiation

Trans-differentiation from an adenocarcinoma to a small cell neuroendocrine state is associated with therapy resistance in multiple cancer types. To gain insight into the underlying molecular events of the trans-differentiation, we perform a multi-omics time course analysis of a pan-small cell neuroendocrine cancer model (termed PARCB), a forward genetic transformation using human prostate basal cells and identify a shared developmental, arc-like, and entropy-high trajectory among all transformation model replicates. Further mapping with single cell resolution reveals two distinct lineages defined by mutually exclusive expression of ASCL1 or ASCL2. Temporal regulation by groups of transcription factors across developmental stages reveals that cellular reprogramming precedes the induction of neuronal programs. TFAP4 and ASCL1/2 feedback are identified as potential regulators of ASCL1 and ASCL2 expression. Our study provides temporal transcriptional patterns and uncovers pan-tissue parallels between prostate and lung cancers, as well as connections to normal neuroendocrine cell states.

The non-specific lethal complex regulates genes and pathways genetically linked to Parkinson's disease

Genetic variants conferring risks for Parkinson's disease have been highlighted through genome-wide association studies, yet exploration of their specific disease mechanisms is lacking. Two Parkinson's disease candidate genes, KAT8 and KANSL1, identified through genome-wide studies and a PINK1-mitophagy screen, encode part of the histone acetylating non-specific lethal complex. This complex localizes to the nucleus, where it plays a role in transcriptional activation, and to mitochondria, where it has been suggested to have a role in mitochondrial transcription. In this study, we sought to identify whether the non-specific lethal complex has potential regulatory relationships with other genes associated with Parkinson's disease in human brain. Correlation in the expression of non-specific lethal genes and Parkinson's disease-associated genes was investigated in primary gene co-expression networks using publicly-available transcriptomic data from multiple brain regions (provided by the Genotype-Tissue Expression Consortium and UK Brain Expression Consortium), whilst secondary networks were used to examine cell type specificity. Reverse engineering of gene regulatory networks generated regulons of the complex, which were tested for heritability using stratified linkage disequilibrium score regression. Prioritized gene targets were then validated in vitro using a QuantiGene multiplex assay and publicly-available chromatin immunoprecipitation-sequencing data. Significant clustering of non-specific lethal genes was revealed alongside Parkinson's disease-associated genes in frontal cortex primary co-expression modules, amongst other brain regions. Both primary and secondary co-expression modules containing these genes were enriched for mainly neuronal cell types. Regulons of the complex contained Parkinson's disease-associated genes and were enriched for biological pathways genetically linked to disease. When examined in a neuroblastoma cell line, 41% of prioritized gene targets showed significant changes in mRNA expression following KANSL1 or KAT8 perturbation. KANSL1 and H4K8 chromatin immunoprecipitation-sequencing data demonstrated non-specific lethal complex activity at many of these genes. In conclusion, genes encoding the non-specific lethal complex are highly correlated with and regulate genes associated with Parkinson's disease. Overall, these findings reveal a potentially wider role for this protein complex in regulating genes and pathways implicated in Parkinson's disease.

Gene essentiality in cancer is better predicted by mRNA abundance than by gene regulatory network-inferred activity

Gene regulatory networks (GRNs) are often deregulated in tumor cells, resulting in altered transcriptional programs that facilitate tumor growth. These altered networks may make tumor cells vulnerable to the inhibition of specific regulatory proteins. Consequently, the reconstruction of GRNs in tumors is often proposed as a means to identify therapeutic targets. While there are examples of individual targets identified using GRNs, the extent to which GRNs can be used to predict sensitivity to targeted intervention in general remains unknown. Here we use the results of genome-wide CRISPR screens to systematically assess the ability of GRNs to predict sensitivity to gene inhibition in cancer cell lines. Using GRNs derived from multiple sources, including GRNs reconstructed from tumor transcriptomes and from curated databases, we infer regulatory gene activity in cancer cell lines from ten cancer types. We then ask, in each cancer type, if the inferred regulatory activity of each gene is predictive of sensitivity to CRISPR perturbation of that gene. We observe slight variation in the correlation between gene regulatory activity and gene sensitivity depending on the source of the GRN and the activity estimation method used. However, we find that there is consistently a stronger relationship between mRNA abundance and gene sensitivity than there is between regulatory gene activity and gene sensitivity. This is true both when gene sensitivity is treated as a binary and a quantitative property. Overall, our results suggest that gene sensitivity is better predicted by measured expression than by GRN-inferred activity.

A Bayesian noisy logic model for inference of transcription factor activity from single cell and bulk transcriptomic data

The advent of high-throughput sequencing has made it possible to measure the expression of genes at relatively low cost. However, direct measurement of regulatory mechanisms, such as transcription factor (TF) activity is still not readily feasible in a high-throughput manner. Consequently, there is a need for computational approaches that can reliably estimate regulator activity from observable gene expression data. In this work, we present a noisy Boolean logic Bayesian model for TF activity inference from differential gene expression data and causal graphs. Our approach provides a flexible framework to incorporate biologically motivated TF-gene regulation logic models. Using simulations and controlled over-expression experiments in cell cultures, we demonstrate that our method can accurately identify TF activity. Moreover, we apply our method to bulk and single cell transcriptomics measurements to investigate transcriptional regulation of fibroblast phenotypic plasticity. Finally, to facilitate usage, we provide user-friendly software packages and a web-interface to query TF activity from user input differential gene expression data: https://umbibio.math.umb.edu/nlbayes/.

Anti-PD-1 immunotherapy with androgen deprivation therapy induces robust immune infiltration in metastatic castration-sensitive prostate cancer

When compared to other malignancies, the tumor microenvironment (TME) of primary and castration-resistant prostate cancer (CRPC) is relatively devoid of immune infiltrates. While androgen deprivation therapy (ADT) induces a complex immune infiltrate in localized prostate cancer, the composition of the TME in metastatic castration-sensitive prostate cancer (mCSPC), and the effects of ADT and other treatments in this context are poorly understood. Here, we perform a comprehensive single-cell RNA sequencing (scRNA-seq) profiling of metastatic sites from patients participating in a phase 2 clinical trial (NCT03951831) that evaluated standard-of-care chemo-hormonal therapy combined with anti-PD-1 immunotherapy. We perform a longitudinal, protein activity-based analysis of TME subpopulations, revealing immune subpopulations conserved across multiple metastatic sites. We also observe dynamic changes in these immune subpopulations in response to treatment and a correlation with clinical outcomes. Our study uncovers a therapy-resistant, transcriptionally distinct tumor subpopulation that expands in cell number in treatment-refractory patients.

Computational single cell oncology: state of the art

Single cell computational analysis has emerged as a powerful tool in the field of oncology, enabling researchers to decipher the complex cellular heterogeneity that characterizes cancer. By leveraging computational algorithms and bioinformatics approaches, this methodology provides insights into the underlying genetic, epigenetic and transcriptomic variations among individual cancer cells. In this paper, we present a comprehensive overview of single cell computational analysis in oncology, discussing the key computational techniques employed for data processing, analysis, and interpretation. We explore the challenges associated with single cell data, including data quality control, normalization, dimensionality reduction, clustering, and trajectory inference. Furthermore, we highlight the applications of single cell computational analysis, including the identification of novel cell states, the characterization of tumor subtypes, the discovery of biomarkers, and the prediction of therapy response. Finally, we address the future directions and potential advancements in the field, including the development of machine learning and deep learning approaches for single cell analysis. Overall, this paper aims to provide a roadmap for researchers interested in leveraging computational methods to unlock the full potential of single cell analysis in understanding cancer biology with the goal of advancing precision oncology. For this purpose, we also include a notebook that instructs on how to apply the recommended tools in the Preprocessing and Quality Control section.

Melatonin improves the vitrification of sheep morulae by modulating transcriptome

Embryo vitrification technology is widely used in livestock production, but freezing injury has been a key factor hindering the efficiency of embryo production. There is an urgent need to further analyze the molecular mechanism of embryo damage by the vitrification process. In the study, morulae were collected from Hu sheep uterine horns after superovulation and sperm transfusion. Morulae were Cryotop vitrified and warmed. Nine morulae were in the vitrified control group (frozen), and seven morulae were vitrified and warmed with 10-5 M melatonin (melatonin). Eleven non-frozen morulae were used as controls (fresh). After warming, each embryo was sequenced separately for library construction and gene expression analysis. p < 0.05 was used to differentiate differentially expressed genes (DEG). The results showed that differentiated differentially expressed genes (DEG) in vitrified morulae were mainly enriched in protein kinase activity, adhesion processes, calcium signaling pathways and Wnt, PI3K/AKT, Ras, ErbB, and MAPK signaling pathways compared to controls. Importantly, melatonin treatment upregulated the expression of key pathways that increase the resistance of morulae against vitrification induced damage. These pathways include kinase activity pathway, ErbB, and PI3K/Akt signaling pathway. It is worth mentioning that melatonin upregulates the expression of XPA, which is a key transcription factor for DNA repair. In conclusion, vitrification affected the transcriptome of in vivo-derived Hu sheep morulae, and melatonin had a protective effect on the vitrification process. For the first time, the transcriptome profiles caused by vitrification and melatonin in sheep morulae were analyzed in single embryo level. These data obtained from the single embryo level provide an important molecular mechanism for further optimizing the cryopreservation of embryos or other cells.

Adjustment of spurious correlations in co-expression measurements from RNA-Sequencing data

MOTIVATION

Gene co-expression measurements are widely used in computational biology to identify coordinated expression patterns across a group of samples. Coordinated expression of genes may indicate that they are controlled by the same transcriptional regulatory program, or involved in common biological processes. Gene co-expression is generally estimated from RNA-Sequencing data, which are commonly normalized to remove technical variability. Here, we demonstrate that certain normalization methods, in particular quantile-based methods, can introduce false-positive associations between genes. These false-positive associations can consequently hamper downstream co-expression network analysis. Quantile-based normalization can, however, be extremely powerful. In particular, when preprocessing large-scale heterogeneous data, quantile-based normalization methods such as smooth quantile normalization can be applied to remove technical variability while maintaining global differences in expression for samples with different biological attributes.

RESULTS

We developed SNAIL (Smooth-quantile Normalization Adaptation for the Inference of co-expression Links), a normalization method based on smooth quantile normalization specifically designed for modeling of co-expression measurements. We show that SNAIL avoids formation of false-positive associations in co-expression as well as in downstream network analyses. Using SNAIL, one can avoid arbitrary gene filtering and retain associations to genes that only express in small subgroups of samples. This highlights the method's potential future impact on network modeling and other association-based approaches in large-scale heterogeneous data.

AVAILABILITY AND IMPLEMENTATION

The implementation of the SNAIL algorithm and code to reproduce the analyses described in this work can be found in the GitHub repository https://github.com/kuijjerlab/PySNAIL.

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.

Multi-omics regulatory network inference in the presence of missing data

A key problem in systems biology is the discovery of regulatory mechanisms that drive phenotypic behaviour of complex biological systems in the form of multi-level networks. Modern multi-omics profiling techniques probe these fundamental regulatory networks but are often hampered by experimental restrictions leading to missing data or partially measured omics types for subsets of individuals due to cost restrictions. In such scenarios, in which missing data is present, classical computational approaches to infer regulatory networks are limited. In recent years, approaches have been proposed to infer sparse regression models in the presence of missing information. Nevertheless, these methods have not been adopted for regulatory network inference yet. In this study, we integrated regression-based methods that can handle missingness into KiMONo, a Knowledge guided Multi-Omics Network inference approach, and benchmarked their performance on commonly encountered missing data scenarios in single- and multi-omics studies. Overall, two-step approaches that explicitly handle missingness performed best for a wide range of random- and block-missingness scenarios on imbalanced omics-layers dimensions, while methods implicitly handling missingness performed best on balanced omics-layers dimensions. Our results show that robust multi-omics network inference in the presence of missing data with KiMONo is feasible and thus allows users to leverage available multi-omics data to its full extent.

Master regulator activity QTL protocol to implicate regulatory pathways potentially mediating GWAS signals using eQTL data

Here, we present a protocol to identify transcriptional regulators potentially mediating downstream biological effects of germline variants associated with complex traits of interest, which enables functional hypothesis generation independent of colocalizing expression quantitative trait loci (eQTLs). We describe steps for tissue-/cell-type-specific co-expression network modeling, expression regulator activity inference, and identification of representative phenotypic master regulators. Finally, we detail activity QTL and eQTL analyses. This protocol requires genotype, expression, and relevant covariables and phenotype data from existing eQTL datasets. For complete details on the use and execution of this protocol, please refer to Hoskins et al.1.

Proteomic Atlas of Atherosclerosis: The Contribution of Proteoglycans to Sex Differences, Plaque Phenotypes, and Outcomes

BACKGROUND

Using proteomics, we aimed to reveal molecular types of human atherosclerotic lesions and study their associations with histology, imaging, and cardiovascular outcomes.

METHODS

Two hundred nineteen carotid endarterectomy samples were procured from 120 patients. A sequential protein extraction protocol was employed in conjunction with multiplexed, discovery proteomics. To focus on extracellular proteins, parallel reaction monitoring was employed for targeted proteomics. Proteomic signatures were integrated with bulk, single-cell, and spatial RNA-sequencing data, and validated in 200 patients from the Athero-Express Biobank study.

RESULTS

This extensive proteomics analysis identified plaque inflammation and calcification signatures, which were inversely correlated and validated using targeted proteomics. The inflammation signature was characterized by the presence of neutrophil-derived proteins, such as S100A8/9 (calprotectin) and myeloperoxidase, whereas the calcification signature included fetuin-A, osteopontin, and gamma-carboxylated proteins. The proteomics data also revealed sex differences in atherosclerosis, with large-aggregating proteoglycans versican and aggrecan being more abundant in females and exhibiting an inverse correlation with estradiol levels. The integration of RNA-sequencing data attributed the inflammation signature predominantly to neutrophils and macrophages, and the calcification and sex signatures to smooth muscle cells, except for certain plasma proteins that were not expressed but retained in plaques, such as fetuin-A. Dimensionality reduction and machine learning techniques were applied to identify 4 distinct plaque phenotypes based on proteomics data. A protein signature of 4 key proteins (calponin, protein C, serpin H1, and versican) predicted future cardiovascular mortality with an area under the curve of 75% and 67.5% in the discovery and validation cohort, respectively, surpassing the prognostic performance of imaging and histology.

CONCLUSIONS

Plaque proteomics redefined clinically relevant patient groups with distinct outcomes, identifying subgroups of male and female patients with elevated risk of future cardiovascular events.

Identifying genetic regulatory variants that affect transcription factor activity

Genetic variants affecting gene expression levels in humans have been mapped in the Genotype-Tissue Expression (GTEx) project. Trans-acting variants impacting many genes simultaneously through a shared transcription factor (TF) are of particular interest. Here, we developed a generalized linear model (GLM) to estimate protein-level TF activity levels in an individual-specific manner from GTEx RNA sequencing (RNA-seq) profiles. It uses observed differential gene expression after TF perturbation as a predictor and, by analyzing differential expression within pairs of neighboring genes, controls for the confounding effect of variation in chromatin state along the genome. We inferred genotype-specific activities for 55 TFs across 49 tissues. Subsequently performing genome-wide association analysis on this virtual trait revealed TF activity quantitative trait loci (aQTLs) that, as a set, are enriched for functional features. Altogether, the set of tools we introduce here highlights the potential of genetic association studies for cellular endophenotypes based on a network-based multi-omics approach. The transparent peer review record is available.

Pancancer network analysis reveals key master regulators for cancer invasiveness

BACKGROUND

Tumor invasiveness reflects numerous biological changes, including tumorigenesis, progression, and metastasis. To decipher the role of transcriptional regulators (TR) involved in tumor invasiveness, we performed a systematic network-based pan-cancer assessment of master regulators of cancer invasiveness.

MATERIALS AND METHODS

We stratified patients in The Cancer Genome Atlas (TCGA) into invasiveness high (INV-H) and low (INV-L) groups using consensus clustering based on an established robust 24-gene signature to determine the prognostic association of invasiveness with overall survival (OS) across 32 different cancers. We devise a network-based protocol to identify TRs as master regulators (MRs) unique to INV-H and INV-L phenotypes. We validated the activity of MRs coherently associated with INV-H phenotype and worse OS across cancers in TCGA on a series of additional datasets in the Prediction of Clinical Outcomes from the Genomic Profiles (PRECOG) repository.

RESULTS

Based on the 24-gene signature, we defined the invasiveness score for each patient sample and stratified patients into INV-H and INV-L clusters. We observed that invasiveness was associated with worse survival outcomes in almost all cancers and had a significant association with OS in ten out of 32 cancers. Our network-based framework identified common invasiveness-associated MRs specific to INV-H and INV-L groups across the ten prognostic cancers, including COL1A1, which is also part of the 24-gene signature, thus acting as a positive control. Downstream pathway analysis of MRs specific to INV-H phenotype resulted in the identification of several enriched pathways, including Epithelial into Mesenchymal Transition, TGF-β signaling pathway, regulation of Toll-like receptors, cytokines, and inflammatory response, and selective expression of chemokine receptors during T-cell polarization. Most of these pathways have connotations of inflammatory immune response and feasibility for metastasis.

CONCLUSION

Our pan-cancer study provides a comprehensive master regulator analysis of tumor invasiveness and can suggest more precise therapeutic strategies by targeting the identified MRs and downstream enriched pathways for patients across multiple cancers.

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.

Uncovering Protein Networks in Cardiovascular Proteomics

Biological networks have been widely used in many different diseases to identify potential biomarkers and design drug targets. In the present review, we describe the main computational techniques for reconstructing and analyzing different types of protein networks and summarize the previous applications of such techniques in cardiovascular diseases. Existing tools are critically compared, discussing when each method is preferred such as the use of co-expression networks for functional annotation of protein clusters and the use of directed networks for inferring regulatory associations. Finally, we are presenting examples of reconstructing protein networks of different types (regulatory, co-expression, and protein-protein interaction networks). We demonstrate the necessity to reconstruct networks separately for each cardiovascular tissue type and disease entity and provide illustrative examples of the importance of taking into consideration relevant post-translational modifications. Finally, we demonstrate and discuss how the findings of protein networks could be interpreted using single-cell RNA-sequencing data.

Molecular mechanisms regulating natural menopause in the female ovary: a study based on transcriptomic data

Introduction

Natural menopause is an inevitable biological process with significant implications for women's health. However, the molecular mechanisms underlying menopause are not well understood. This study aimed to investigate the molecular and cellular changes occurring in the ovary before and after perimenopause.

Methods

Single-cell sequencing data from the GTEx V8 cohort (30-39: 14 individuals; 40-49: 37 individuals; 50-59: 61 individuals) and transcriptome sequencing data from ovarian tissue were analyzed. Seurat was used for single-cell sequencing data analysis, while harmony was employed for data integration. Cell differentiation trajectories were inferred using CytoTrace. CIBERSORTX assessed cell infiltration scores in ovarian tissue. WGCNA evaluated co-expression network characteristics in pre- and post-perimenopausal ovarian tissue. Functional enrichment analysis of co-expression modules was conducted using ClusterprofileR and Metascape. DESeq2 performed differential expression analysis. Master regulator analysis and signaling pathway activity analysis were carried out using MsViper and Progeny, respectively. Machine learning models were constructed using Orange3.

Results

We identified the differentiation trajectory of follicular cells in the ovary as ARID5B+ Granulosa -> JUN+ Granulosa -> KRT18+ Granulosa -> MT-CO2+ Granulosa -> GSTA1+ Granulosa -> HMGB1+ Granulosa. Genes driving Granulosa differentiation, including RBP1, TMSB10, SERPINE2, and TMSB4X, were enriched in ATP-dependent activity regulation pathways. Genes involved in maintaining the Granulosa state, such as DCN, ARID5B, EIF1, and HSP90AB1, were enriched in the response to unfolded protein and chaperone-mediated protein complex assembly pathways. Increased contents of terminally differentiated HMGB1+ Granulosa and GSTA1+ Granulosa were observed in the ovaries of individuals aged 50-69. Signaling pathway activity analysis indicated a gradual decrease in TGFb and MAPK pathway activity with menopause progression, while p53 pathway activity increased. Master regulator analysis revealed significant activation of transcription factors FOXR1, OTX2, MYBL2, HNF1A, and FOXN4 in the 30-39 age group, and GLI1, SMAD1, SMAD7, APP, and EGR1 in the 40-49 age group. Additionally, a diagnostic model based on 16 transcription factors (Logistic Regression L2) achieved reliable performance in determining ovarian status before and after perimenopause.

Conclusion

This study provides insights into the molecular and cellular mechanisms underlying natural menopause in the ovary. The findings contribute to our understanding of perimenopausal changes and offer a foundation for health management strategies for women during this transition.

Identification of transcriptional programs using dense vector representations defined by mutual information with GeneVector

Deciphering individual cell phenotypes from cell-specific transcriptional processes requires high dimensional single cell RNA sequencing. However, current dimensionality reduction methods aggregate sparse gene information across cells, without directly measuring the relationships that exist between genes. By performing dimensionality reduction with respect to gene co-expression, low-dimensional features can model these gene-specific relationships and leverage shared signal to overcome sparsity. We describe GeneVector, a scalable framework for dimensionality reduction implemented as a vector space model using mutual information between gene expression. Unlike other methods, including principal component analysis and variational autoencoders, GeneVector uses latent space arithmetic in a lower dimensional gene embedding to identify transcriptional programs and classify cell types. In this work, we show in four single cell RNA-seq datasets that GeneVector was able to capture phenotype-specific pathways, perform batch effect correction, interactively annotate cell types, and identify pathway variation with treatment over time.

NSD2 maintains lineage plasticity and castration-resistance in neuroendocrine prostate cancer

The clinical use of potent androgen receptor (AR) inhibitors has promoted the emergence of novel subtypes of metastatic castration-resistant prostate cancer (mCRPC), including neuroendocrine prostate cancer (CRPC-NE), which is highly aggressive and lethal 1 . These mCRPC subtypes display increased lineage plasticity and often lack AR expression 2-5 . Here we show that neuroendocrine differentiation and castration-resistance in CRPC-NE are maintained by the activity of Nuclear Receptor Binding SET Domain Protein 2 (NSD2) 6 , which catalyzes histone H3 lysine 36 dimethylation (H3K36me2). We find that organoid lines established from genetically-engineered mice 7 recapitulate key features of human CRPC-NE, and can display transdifferentiation to neuroendocrine states in culture. CRPC-NE organoids express elevated levels of NSD2 and H3K36me2 marks, but relatively low levels of H3K27me3, consistent with antagonism of EZH2 activity by H3K36me2. Human CRPC-NE but not primary NEPC tumors expresses high levels of NSD2, consistent with a key role for NSD2 in lineage plasticity, and high NSD2 expression in mCRPC correlates with poor survival outcomes. Notably, CRISPR/Cas9 targeting of NSD2 or expression of a dominant-negative oncohistone H3.3K36M mutant results in loss of neuroendocrine phenotypes and restores responsiveness to the AR inhibitor enzalutamide in mouse and human CRPC-NE organoids and grafts. Our findings indicate that NSD2 inhibition can reverse lineage plasticity and castration-resistance, and provide a potential new therapeutic target for CRPC-NE.

Discovering cellular programs of intrinsic and extrinsic drivers of metabolic traits using LipocyteProfiler

A primary obstacle in translating genetic associations with disease into therapeutic strategies is elucidating the cellular programs affected by genetic risk variants and effector genes. Here, we introduce LipocyteProfiler, a cardiometabolic-disease-oriented high-content image-based profiling tool that enables evaluation of thousands of morphological and cellular profiles that can be systematically linked to genes and genetic variants relevant to cardiometabolic disease. We show that LipocyteProfiler allows surveillance of diverse cellular programs by generating rich context- and process-specific cellular profiles across hepatocyte and adipocyte cell-state transitions. We use LipocyteProfiler to identify known and novel cellular mechanisms altered by polygenic risk of metabolic disease, including insulin resistance, fat distribution, and the polygenic contribution to lipodystrophy. LipocyteProfiler paves the way for large-scale forward and reverse deep phenotypic profiling in lipocytes and provides a framework for the unbiased identification of causal relationships between genetic variants and cellular programs relevant to human disease.

Mechanistic insights into the interactions between cancer drivers and the tumour immune microenvironment

BACKGROUND

The crosstalk between cancer and the tumour immune microenvironment (TIME) has attracted significant interest in the latest years because of its impact on cancer evolution and response to treatment. Despite this, cancer-specific tumour-TIME interactions and their mechanistic insights are still poorly understood.

METHODS

Here, we compute the significant interactions occurring between cancer-specific genetic drivers and five anti- and pro-tumour TIME features in 32 cancer types using Lasso regularised ordinal regression. Focusing on head and neck squamous cancer (HNSC), we rebuild the functional networks linking specific TIME driver alterations to the TIME state they associate with.

RESULTS

The 477 TIME drivers that we identify are multifunctional genes whose alterations are selected early in cancer evolution and recur across and within cancer types. Tumour suppressors and oncogenes have an opposite effect on the TIME and the overall anti-tumour TIME driver burden is predictive of response to immunotherapy. TIME driver alterations predict the immune profiles of HNSC molecular subtypes, and perturbations in keratinization, apoptosis and interferon signalling underpin specific driver-TIME interactions.

CONCLUSIONS

Overall, our study delivers a comprehensive resource of TIME drivers, gives mechanistic insights into their immune-regulatory role, and provides an additional framework for patient prioritisation to immunotherapy. The full list of TIME drivers and associated properties are available at http://www.network-cancer-genes.org .

Therapeutic blood-brain barrier modulation and stroke treatment by a bioengineered FZD4-selective WNT surrogate in mice

Derangements of the blood-brain barrier (BBB) or blood-retinal barrier (BRB) occur in disorders ranging from stroke, cancer, diabetic retinopathy, and Alzheimer's disease. The Norrin/FZD4/TSPAN12 pathway activates WNT/β-catenin signaling, which is essential for BBB and BRB function. However, systemic pharmacologic FZD4 stimulation is hindered by obligate palmitoylation and insolubility of native WNTs and suboptimal properties of the FZD4-selective ligand Norrin. Here, we develop L6-F4-2, a non-lipidated, FZD4-specific surrogate which significantly improves subpicomolar affinity versus native Norrin. In Norrin knockout (NdpKO) mice, L6-F4-2 not only potently reverses neonatal retinal angiogenesis deficits, but also restores BRB and BBB function. In adult C57Bl/6J mice, post-stroke systemic delivery of L6-F4-2 strongly reduces BBB permeability, infarction, and edema, while improving neurologic score and capillary pericyte coverage. Our findings reveal systemic efficacy of a bioengineered FZD4-selective WNT surrogate during ischemic BBB dysfunction, with potential applicability to adult CNS disorders characterized by an aberrant blood-brain barrier.

MCPNet: a parallel maximum capacity-based genome-scale gene network construction framework

MOTIVATION

Gene network reconstruction from gene expression profiles is a compute- and data-intensive problem. Numerous methods based on diverse approaches including mutual information, random forests, Bayesian networks, correlation measures, as well as their transforms and filters such as data processing inequality, have been proposed. However, an effective gene network reconstruction method that performs well in all three aspects of computational efficiency, data size scalability, and output quality remains elusive. Simple techniques such as Pearson correlation are fast to compute but ignore indirect interactions, while more robust methods such as Bayesian networks are prohibitively time consuming to apply to tens of thousands of genes.

RESULTS

We developed maximum capacity path (MCP) score, a novel maximum-capacity-path-based metric to quantify the relative strengths of direct and indirect gene-gene interactions. We further present MCPNet, an efficient, parallelized gene network reconstruction software based on MCP score, to reverse engineer networks in unsupervised and ensemble manners. Using synthetic and real Saccharomyces cervisiae datasets as well as real Arabidopsis thaliana datasets, we demonstrate that MCPNet produces better quality networks as measured by AUPRC, is significantly faster than all other gene network reconstruction software, and also scales well to tens of thousands of genes and hundreds of CPU cores. Thus, MCPNet represents a new gene network reconstruction tool that simultaneously achieves quality, performance, and scalability requirements.

AVAILABILITY AND IMPLEMENTATION

Source code freely available for download at https://doi.org/10.5281/zenodo.6499747 and https://github.com/AluruLab/MCPNet, implemented in C++ and supported on Linux.

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.

A Bayesian Noisy Logic Model for Inference of Transcription Factor Activity from Single Cell and Bulk Transcriptomic Data

The advent of high-throughput sequencing has made it possible to measure the expression of genes at relatively low cost. However, direct measurement of regulatory mechanisms, such as Transcription Factor (TF) activity is still not readily feasible in a high-throughput manner. Consequently, there is a need for computational approaches that can reliably estimate regulator activity from observable gene expression data. In this work, we present a noisy Boolean logic Bayesian model for TF activity inference from differential gene expression data and causal graphs. Our approach provides a flexible framework to incorporate biologically motivated TF-gene regulation logic models. Using simulations and controlled over-expression experiments in cell cultures, we demonstrate that our method can accurately identify TF activity. Moreover, we apply our method to bulk and single cell transcriptomics measurements to investigate transcriptional regulation of fibroblast phenotypic plasticity. Finally, to facilitate usage, we provide user-friendly software packages and a web-interface to query TF activity from user input differential gene expression data: https://umbibio.math.umb.edu/nlbayes/.

A gene regulatory network inference model based on pseudo-siamese network

MOTIVATION

Gene regulatory networks (GRNs) arise from the intricate interactions between transcription factors (TFs) and their target genes during the growth and development of organisms. The inference of GRNs can unveil the underlying gene interactions in living systems and facilitate the investigation of the relationship between gene expression patterns and phenotypic traits. Although several machine-learning models have been proposed for inferring GRNs from single-cell RNA sequencing (scRNA-seq) data, some of these models, such as Boolean and tree-based networks, suffer from sensitivity to noise and may encounter difficulties in handling the high noise and dimensionality of actual scRNA-seq data, as well as the sparse nature of gene regulation relationships. Thus, inferring large-scale information from GRNs remains a formidable challenge.

RESULTS

This study proposes a multilevel, multi-structure framework called a pseudo-Siamese GRN (PSGRN) for inferring large-scale GRNs from time-series expression datasets. Based on the pseudo-Siamese network, we applied a gated recurrent unit to capture the time features of each TF and target matrix and learn the spatial features of the matrices after merging by applying the DenseNet framework. Finally, we applied a sigmoid function to evaluate interactions. We constructed two maize sub-datasets, including gene expression levels and GRNs, using existing open-source maize multi-omics data and compared them to other GRN inference methods, including GENIE3, GRNBoost2, nonlinear ordinary differential equations, CNNC, and DGRNS. Our results show that PSGRN outperforms state-of-the-art methods. This study proposed a new framework: a PSGRN that allows GRNs to be inferred from scRNA-seq data, elucidating the temporal and spatial features of TFs and their target genes. The results show the model's robustness and generalization, laying a theoretical foundation for maize genotype-phenotype associations with implications for breeding work.

Bayesian Machine Learning Enables Identification of Transcriptional Network Disruptions Associated with Drug-Resistant Prostate Cancer

Survival rates of patients with metastatic castration-resistant prostate cancer (mCRPC) are low due to lack of response or acquired resistance to available therapies, such as abiraterone (Abi). A better understanding of the underlying molecular mechanisms is needed to identify effective targets to overcome resistance. Given the complexity of the transcriptional dynamics in cells, differential gene expression analysis of bulk transcriptomics data cannot provide sufficient detailed insights into resistance mechanisms. Incorporating network structures could overcome this limitation to provide a global and functional perspective of Abi resistance in mCRPC. Here, we developed TraRe, a computational method using sparse Bayesian models to examine phenotypically driven transcriptional mechanistic differences at three distinct levels: transcriptional networks, specific regulons, and individual transcription factors (TF). TraRe was applied to transcriptomic data from 46 patients with mCRPC with Abi-response clinical data and uncovered abrogated immune response transcriptional modules that showed strong differential regulation in Abi-responsive compared with Abi-resistant patients. These modules were replicated in an independent mCRPC study. Furthermore, key rewiring predictions and their associated TFs were experimentally validated in two prostate cancer cell lines with different Abi-resistance features. Among them, ELK3, MXD1, and MYB played a differential role in cell survival in Abi-sensitive and Abi-resistant cells. Moreover, ELK3 regulated cell migration capacity, which could have a direct impact on mCRPC. Collectively, these findings shed light on the underlying transcriptional mechanisms driving Abi response, demonstrating that TraRe is a promising tool for generating novel hypotheses based on identified transcriptional network disruptions.

SIGNIFICANCE

The computational method TraRe built on Bayesian machine learning models for investigating transcriptional network structures shows that disruption of ELK3, MXD1, and MYB signaling cascades impacts abiraterone resistance in prostate cancer.

Epigenomic mapping identifies an enhancer repertoire that regulates cell identity in bladder cancer through distinct transcription factor networks

Muscle-invasive bladder cancer (BLCA) is an aggressive disease. Consensus BLCA transcriptomic subtypes have been proposed, with two major Luminal and Basal subgroups, presenting distinct molecular and clinical characteristics. However, how these distinct subtypes are regulated remains unclear. We hypothesized that epigenetic activation of distinct super-enhancers could drive the transcriptional programs of BLCA subtypes. Through integrated RNA-sequencing and epigenomic profiling of histone marks in primary tumours, cancer cell lines, and normal human urothelia, we established the first integrated epigenetic map of BLCA and demonstrated the link between subtype and epigenetic control. We identified the repertoire of activated super-enhancers and highlighted Basal, Luminal and Normal-associated SEs. We revealed super-enhancer-regulated networks of candidate master transcription factors for Luminal and Basal subgroups including FOXA1 and ZBED2, respectively. FOXA1 CRISPR-Cas9 mutation triggered a shift from Luminal to Basal phenotype, confirming its role in Luminal identity regulation and induced ZBED2 overexpression. In parallel, we showed that both FOXA1 and ZBED2 play concordant roles in preventing inflammatory response in cancer cells through STAT2 inhibition. Our study furthers the understanding of epigenetic regulation of muscle-invasive BLCA and identifies a co-regulated network of super-enhancers and associated transcription factors providing potential targets for the treatment of this aggressive disease.

NaRnEA: An Information Theoretic Framework for Gene Set Analysis

Gene sets are being increasingly leveraged to make high-level biological inferences from transcriptomic data; however, existing gene set analysis methods rely on overly conservative, heuristic approaches for quantifying the statistical significance of gene set enrichment. We created Nonparametric analytical-Rank-based Enrichment Analysis (NaRnEA) to facilitate accurate and robust gene set analysis with an optimal null model derived using the information theoretic Principle of Maximum Entropy. By measuring the differential activity of ~2500 transcriptional regulatory proteins based on the differential expression of each protein's transcriptional targets between primary tumors and normal tissue samples in three cohorts from The Cancer Genome Atlas (TCGA), we demonstrate that NaRnEA critically improves in two widely used gene set analysis methods: Gene Set Enrichment Analysis (GSEA) and analytical-Rank-based Enrichment Analysis (aREA). We show that the NaRnEA-inferred differential protein activity is significantly correlated with differential protein abundance inferred from independent, phenotype-matched mass spectrometry data in the Clinical Proteomic Tumor Analysis Consortium (CPTAC), confirming the statistical and biological accuracy of our approach. Additionally, our analysis crucially demonstrates that the sample-shuffling empirical null models leveraged by GSEA and aREA for gene set analysis are overly conservative, a shortcoming that is avoided by the newly developed Maximum Entropy analytical null model employed by NaRnEA.

Gene Functional Networks from Time Expression Profiles: A Constructive Approach Demonstrated in Chili Pepper (Capsicum annuum L.)

Gene co-expression networks are powerful tools to understand functional interactions between genes. However, large co-expression networks are difficult to interpret and do not guarantee that the relations found will be true for different genotypes. Statistically verified time expression profiles give information about significant changes in expressions through time, and genes with highly correlated time expression profiles, which are annotated in the same biological process, are likely to be functionally connected. A method to obtain robust networks of functionally related genes will be useful to understand the complexity of the transcriptome, leading to biologically relevant insights. We present an algorithm to construct gene functional networks for genes annotated in a given biological process or other aspects of interest. We assume that there are genome-wide time expression profiles for a set of representative genotypes of the species of interest. The method is based on the correlation of time expression profiles, bound by a set of thresholds that assure both, a given false discovery rate, and the discard of correlation outliers. The novelty of the method consists in that a gene expression relation must be repeatedly found in a given set of independent genotypes to be considered valid. This automatically discards relations particular to specific genotypes, assuring a network robustness, which can be set a priori. Additionally, we present an algorithm to find transcription factors candidates for regulating hub genes within a network. The algorithms are demonstrated with data from a large experiment studying gene expression during the development of the fruit in a diverse set of chili pepper genotypes. The algorithm is implemented and demonstrated in a new version of the publicly available R package "Salsa" (version 1.0).

PrismEXP: gene annotation prediction from stratified gene-gene co-expression matrices

Background

Gene-gene co-expression correlations measured by mRNA-sequencing (RNA-seq) can be used to predict gene annotations based on the co-variance structure within these data. In our prior work, we showed that uniformly aligned RNA-seq co-expression data from thousands of diverse studies is highly predictive of both gene annotations and protein-protein interactions. However, the performance of the predictions varies depending on whether the gene annotations and interactions are cell type and tissue specific or agnostic. Tissue and cell type-specific gene-gene co-expression data can be useful for making more accurate predictions because many genes perform their functions in unique ways in different cellular contexts. However, identifying the optimal tissues and cell types to partition the global gene-gene co-expression matrix is challenging.

Results

Here we introduce and validate an approach called PRediction of gene Insights from Stratified Mammalian gene co-EXPression (PrismEXP) for improved gene annotation predictions based on RNA-seq gene-gene co-expression data. Using uniformly aligned data from ARCHS4, we apply PrismEXP to predict a wide variety of gene annotations including pathway membership, Gene Ontology terms, as well as human and mouse phenotypes. Predictions made with PrismEXP outperform predictions made with the global cross-tissue co-expression correlation matrix approach on all tested domains, and training using one annotation domain can be used to predict annotations in other domains.

Conclusions

By demonstrating the utility of PrismEXP predictions in multiple use cases we show how PrismEXP can be used to enhance unsupervised machine learning methods to better understand the roles of understudied genes and proteins. To make PrismEXP accessible, it is provided via a user-friendly web interface, a Python package, and an Appyter. AVAILABILITY. The PrismEXP web-based application, with pre-computed PrismEXP predictions, is available from: https://maayanlab.cloud/prismexp; PrismEXP is also available as an Appyter: https://appyters.maayanlab.cloud/PrismEXP/; and as Python package: https://github.com/maayanlab/prismexp.

Network-based elucidation of colon cancer drug resistance by phosphoproteomic time-series analysis

Aberrant signaling pathway activity is a hallmark of tumorigenesis and progression, which has guided targeted inhibitor design for over 30 years. Yet, adaptive resistance mechanisms, induced by rapid, context-specific signaling network rewiring, continue to challenge therapeutic efficacy. By leveraging progress in proteomic technologies and network-based methodologies, over the past decade, we developed VESPA-an algorithm designed to elucidate mechanisms of cell response and adaptation to drug perturbations-and used it to analyze 7-point phosphoproteomic time series from colorectal cancer cells treated with clinically-relevant inhibitors and control media. Interrogation of tumor-specific enzyme/substrate interactions accurately inferred kinase and phosphatase activity, based on their inferred substrate phosphorylation state, effectively accounting for signal cross-talk and sparse phosphoproteome coverage. The analysis elucidated time-dependent signaling pathway response to each drug perturbation and, more importantly, cell adaptive response and rewiring that was experimentally confirmed by CRISPRko assays, suggesting broad applicability to cancer and other diseases.

The role of non-additive gene action on gene expression variation in plant domestication

BACKGROUND

Plant domestication is a remarkable example of rapid phenotypic transformation of polygenic traits, such as organ size. Evidence from a handful of study cases suggests this transformation is due to gene regulatory changes that result in non-additive phenotypes. Employing data from published genetic crosses, we estimated the role of non-additive gene action in the modulation of transcriptional landscapes in three domesticated plants: maize, sunflower, and chili pepper. Using A. thaliana, we assessed the correlation between gene regulatory network (GRN) connectivity properties, transcript abundance variation, and gene action. Finally, we investigated the propagation of non-additive gene action in GRNs.

RESULTS

We compared crosses between domesticated plants and their wild relatives to a set of control crosses that included a pair of subspecies evolving under natural selection and a set of inbred lines evolving under domestication. We found abundance differences on a higher portion of transcripts in crosses between domesticated-wild plants relative to the control crosses. These transcripts showed non-additive gene action more often in crosses of domesticated-wild plants than in our control crosses. This pattern was strong for genes associated with cell cycle and cell fate determination, which control organ size. We found weak but significant negative correlations between the number of targets of trans-acting genes (Out-degree) and both the magnitude of transcript abundance difference a well as the absolute degree of dominance. Likewise, we found that the number of regulators that control a gene's expression (In-degree) is weakly but negatively correlated with the magnitude of transcript abundance differences. We observed that dominant-recessive gene action is highly propagable through GRNs. Finally, we found that transgressive gene action is driven by trans-acting regulators showing additive gene action.

CONCLUSIONS

Our study highlights the role of non-additive gene action on modulating domestication-related traits, such as organ size via regulatory divergence. We propose that GRNs are shaped by regulatory changes at genes with modest connectivity, which reduces the effects of antagonistic pleiotropy. Finally, we provide empirical evidence of the propagation of non-additive gene action in GRNs, which suggests a transcriptional epistatic model for the control of polygenic traits, such as organ size.

Machine learning bridges omics sciences and plant breeding

Some of the biological knowledge obtained from fundamental research will be implemented in applied plant breeding. To bridge basic research and breeding practice, machine learning (ML) holds great promise to translate biological knowledge and omics data into precision-designed plant breeding. Here, we review ML for multi-omics analysis in plants, including data dimensionality reduction, inference of gene-regulation networks, and gene discovery and prioritization. These applications will facilitate understanding trait regulation mechanisms and identifying target genes potentially applicable to knowledge-driven molecular design breeding. We also highlight applications of deep learning in plant phenomics and ML in genomic selection-assisted breeding, such as various ML algorithms that model the correlations among genotypes (genes), phenotypes (traits), and environments, to ultimately achieve data-driven genomic design breeding.

A generic parallel framework for inferring large-scale gene regulatory networks from expression profiles: application to Alzheimer's disease network

The inference of large-scale gene regulatory networks is essential for understanding comprehensive interactions among genes. Most existing methods are limited to reconstructing networks with a few hundred nodes. Therefore, parallel computing paradigms must be leveraged to construct large networks. We propose a generic parallel framework that enables any existing method, without re-engineering, to infer large networks in parallel, guaranteeing quality output. The framework is tested on 15 inference methods (not limited to) employing in silico benchmarks and real-world large expression matrices, followed by qualitative and speedup assessment. The framework does not compromise the quality of the base serial inference method. We rank the candidate methods and use the top-performing method to infer an Alzheimer's Disease (AD) affected network from large expression profiles of a triple transgenic mouse model consisting of 45,101 genes. The resultant network is further explored to obtain hub genes that emerge functionally related to the disease. We partition the network into 41 modules and conduct pathway enrichment analysis, revealing that a good number of participating genes are collectively responsible for several brain disorders, including AD. Finally, we extract the interactions of a few known AD genes and observe that they are periphery genes connected to the network's hub genes. Availability: The R implementation of the framework is downloadable from https://github.com/Netralab/GenericParallelFramework.

CRdb: a comprehensive resource for deciphering chromatin regulators in human

Chromatin regulators (CRs) regulate epigenetic patterns on a partial or global scale, playing a critical role in affecting multi-target gene expression. As chromatin immunoprecipitation sequencing (ChIP-seq) data associated with CRs are rapidly accumulating, a comprehensive resource of CRs needs to be built urgently for collecting, integrating, and processing these data, which can provide abundant annotated information on CR upstream and downstream regulatory analyses as well as CR-related analysis functions. This study established an integrative CR resource, named CRdb (http://cr.liclab.net/crdb/), with the aim of curating a large number of available resources for CRs and providing extensive annotations and analyses of CRs to help biological researchers clarify the regulation mechanism and function of CRs. The CRdb database comprised a total of 647 CRs and 2,591 ChIP-seq samples from more than 300 human tissues and cell types. These samples have been manually curated from NCBI GEO/SRA and ENCODE. Importantly, CRdb provided the abundant and detailed genetic annotations in CR-binding regions based on ChIP-seq. Furthermore, CRdb supported various functional annotations and upstream regulatory information on CRs. In particular, it embedded four types of CR regulatory analyses: CR gene set enrichment, CR-binding genomic region annotation, CR-TF co-occupancy analysis, and CR regulatory axis analysis. CRdb is a useful and powerful resource that can help in exploring the potential functions of CRs and their regulatory mechanism in diseases and biological processes.

TEAD4 is a master regulator of high-risk nasopharyngeal carcinoma

The molecular basis underlying nasopharyngeal carcinoma (NPC) remains unclear. Recent progress in transcriptional regulatory network analysis helps identify the master regulator (MR) proteins that transcriptionally define malignant tumor phenotypes. Here, we investigated transcription factor-target interactions and identified TEA domain transcription factor 4 (TEAD4) as an MR of high-risk NPC. Precisely, TEAD4 promoted NPC migration, invasion and cisplatin resistance, depending on its autopalmitoylation. Mechanistically, YTHDF2 (YTH domain family 2) recognized WTAP (Wilms tumor 1-associating protein)-mediated TEAD4 m6A methylation to facilitate its stability and led to aberrant up-regulation of TEAD4. Up-regulated TEAD4 further drove NPC progression by transcriptionally activating BZW2 (basic leucine zipper and W2 domains 2) to induce the oncogenic AKT pathway. Moreover, the transcriptional activity of TEAD4 was independent of its canonical coactivators YAP/TAZ. Clinically, TEAD4 serves as an independent predictor of unfavorable prognosis and cisplatin response in NPC. Our data revealed the crucial role of TEAD4 in driving tumor malignancy, thus, may provide therapeutic vulnerability in NPC.

Multi-omics inference of differential breast cancer-related transcriptional regulatory network gene hubs between young Black and White patients

OBJECTIVE

Breast cancers (BrCA) are a leading cause of illness and mortality worldwide. Black women have a higher incidence rate relative to white women prior to age 40 years, and a lower incidence rate after 50 years. The objective of this study is to identify -omics differences between the two breast cancer cohorts to better understand the disparities observed in patient outcomes.

MATERIALS AND METHODS

Using Standard SQL, we queried ISB-CGC hosted Google BigQuery tables storing TCGA BrCA gene expression, methylation, and somatic mutation data and analyzed the combined multi-omics results using a variety of methods.

RESULTS

Among Stage II patients 50 years or younger, genes PIK3CA and CDH1 are more frequently mutated in White (W50) than in Black or African American patients (BAA50), while HUWE1, HYDIN, and FBXW7 mutations are more frequent in BAA50. Over-representation analysis (ORA) and Gene Set Enrichment Analysis (GSEA) results indicate that, among others, the Reactome Signaling by ROBO Receptors gene set is enriched in BAA50. Using the Virtual Inference of Protein-activity by Enriched Regulon analysis (VIPER) algorithm, putative top 20 master regulators identified include NUPR1, NFKBIL1, ZBTB17, TEAD1, EP300, TRAF6, CACTIN, and MID2. CACTIN and MID2 are of prognostic value. We identified driver genes, such as OTUB1, with suppressed expression whose DNA methylation status were inversely correlated with gene expression. Networks capturing microRNA and gene expression correlations identified notable microRNA hubs, such as miR-93 and miR-92a-2, expressed at higher levels in BAA50 than in W50.

DISCUSSION/CONCLUSION

The results point to several driver genes as being involved in the observed differences between the cohorts. The findings here form the basis for further mechanistic exploration.

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.