RTNduals: an R/Bioconductor package for analysis of co-regulation and inference of dual regulons

NetREm: Network Regression Embeddings reveal cell-type transcription factor coordination for gene regulation

Transcription factor (TF) coordination plays a key role in target gene (TG) regulation via protein-protein interactions (PPIs) and DNA co-binding to regulatory elements. Single-cell technologies facilitate gene expression measurement for individual cells and cell-type identification, yet the connection between TF coordination and TG regulation of various cell types remains unclear. To address this, we have developed a novel computational approach, Network Regression Embeddings (NetREm), to reveal cell-type TF-TF coordination activities for TG regulation. NetREm leverages network-constrained regularization using prior knowledge of direct and/or indirect PPIs among TFs to analyze single-cell gene expression data. We test NetREm by simulation data and benchmark its performance in 4 real-world applications that have gold standard TF-TG networks available: mouse (mESCs) and simulated human (hESCs) embryonic stem (ESCs), human hematopoietic stem (HSCs), and mouse dendritic (mDCs) cells. Further, we use NetREm to prioritize valid novel TF-TF coordination links in human Peripheral Blood Mononuclear cell (PBMC) sub-types. We apply NetREm to analyze various cell types in both central (CNS) and peripheral (PNS) nerve system (NS) (e.g. neuronal, glial, Schwann cells (SCs)) as well as in Alzheimers disease (AD). Our findings uncover cell-type coordinating TFs and identify new TF-TG candidate links. We validate our top predictions using Cut&Run and knockout loss-of-function expression data in rat/mouse models and compare results with additional functional genomic data, including expression quantitative trait loci (eQTL) and Genome-Wide Association Studies (GWAS) to link genetic variants (single nucleotide polymorphisms (SNPs)) to TF coordination.

Development of gene panel for predicting recurrence in early-stage cervical cancer patients

Cervical cancer (CC) is a common malignancy affecting women worldwide. Our objective was to develop a consensus-based gene panel using multi-omics data that could effectively predict recurrence in early-stage cervical cancer patients. We utilized the "Multi-Omics Consensus Integration Analysis (MOVICS)" package for consensus clustering design to integrate multiple omics datasets and improve the molecular classification landscape of early-stage CC. We identified the "resting and naive" tumor microenvironment (TME) as cancer subtype (CS) 2. Leveraging the feature genes from the CS classifier, we employed machine learning algorithms to identify a gene panel, including ALDH1A1, CLDN10, MUC13, and C10orf99, which could generate a consensus machine learning-driven score (CMLS) for each patient. Stratifying patients into high and low CMLS groups resulted in Kaplan-Meier curves demonstrating a significant difference in recurrence rates between the two groups. This difference remained significant even after adjusting for clinical features in multivariate Cox regression analysis, with the risk ratio of CMLS surpassing that of clinical characteristics. Furthermore, the TME exhibited notable differences between the different CMLS groups, suggesting that patients with low CMLS may exhibit a better response to immunotherapy. This study highlights the potential of the CMLS approach in predicting recurrence in early-stage cervical cancer patients and provides a screening model for selecting patients suitable for immunotherapy.

Transcriptional Profile of Human Pancreatic Acinar Ductal Metaplasia

BACKGROUND AND AIMS

Aberrant acinar to ductal metaplasia (ADM), one of the earliest events involved in exocrine pancreatic cancer development, is typically studied using pancreata from genetically engineered mouse models.

METHODS

We used primary, human pancreatic acinar cells from organ donors to evaluate the transcriptional and pathway profiles during the course of ADM.

RESULTS

Following 6 days of three-dimensional culture on Matrigel, acinar cells underwent morphological and molecular changes indicative of ADM. mRNA from 14 donors' paired cells (day 0, acinar phenotype and day 6, ductal phenotype) was subjected to whole transcriptome sequencing. Acinar cell specific genes were significantly downregulated in the samples from the day 6 cultures while ductal cell-specific genes were upregulated. Several regulons of ADM were identified including transcription factors with reduced activity (PTF1A, RBPJL, and BHLHA15) and those ductal and progenitor transcription factors with increased activity (HNF1B, SOX11, and SOX4). Cells with the ductal phenotype contained higher expression of genes increased in pancreatic cancer while cells with an acinar phenotype had lower expression of cancer-associated genes.

CONCLUSION

Our findings support the relevancy of human in vitro models to study pancreas cancer pathogenesis and exocrine cell plasticity.

Comprehensive Characterization of the Regulatory Landscape of Adrenocortical Carcinoma: Novel Transcription Factors and Targets Associated with Prognosis

We reconstructed a transcriptional regulatory network for adrenocortical carcinoma (ACC) using transcriptomic and clinical data from The Cancer Genome Atlas (TCGA)-ACC cohort. We investigated the association of transcriptional regulatory units (regulons) with overall survival, molecular phenotypes, and immune signatures. We annotated the ACC regulons with cancer hallmarks and assessed single sample regulon activities in the European Network for the Study of Adrenal Tumors (ENSAT) cohort. We found 369 regulons associated with overall survival and subdivided them into four clusters: RC1 and RC2, associated with good prognosis, and RC3 and RC4, associated with worse outcomes. The RC1 and RC3 regulons were highly correlated with the 'Steroid Phenotype,' while the RC2 and RC4 regulons were highly correlated with a molecular proliferation signature. We selected two regulons, NR5A1 (steroidogenic factor 1, SF-1) and CENPA (Centromeric Protein A), that were consistently associated with overall survival for further downstream analyses. The CENPA regulon was the primary regulator of MKI-67 (a marker of proliferation KI-67), while the NR5A1 regulon is a well-described transcription factor (TF) in ACC tumorigenesis. We also found that the ZBTB4 (Zinc finger and BTB domain-containing protein 4) regulon, which is negatively associated with CENPA in our transcriptional regulatory network, is also a druggable anti-tumorigenic TF. We anticipate that the ACC regulons may be used as a reference for further investigations concerning the complex molecular interactions in ACC tumors.

A transcriptional network of cell cycle dysregulation in noninvasive papillary urothelial carcinoma

Human cancers display a restricted set of expression profiles, despite diverse mutational drivers. This has led to the hypothesis that select sets of transcription factors act on similar target genes as an integrated network, buffering a tumor’s transcriptional state. Noninvasive papillary urothelial carcinoma (NIPUC) with higher cell cycle activity has higher risk of recurrence and progression. In this paper, we describe a transcriptional network of cell cycle dysregulation in NIPUC, which was delineated using the ARACNe algorithm applied to expression data from a new cohort (n = 81, RNA sequencing), and two previously published cohorts. The transcriptional network comprised 121 transcription factors, including the pluripotency factors SOX2 and SALL4, the sex hormone binding receptors ESR1 and PGR, and multiple homeobox factors. Of these 121 transcription factors, 65 and 56 were more active in tumors with greater and less cell cycle activity, respectively. When clustered by activity of these transcription factors, tumors divided into High Cell Cycle versus Low Cell Cycle groups. Tumors in the High Cell Cycle group demonstrated greater mutational burden and copy number instability. A putative mutational driver of cell cycle dysregulation, such as homozygous loss of CDKN2A, was found in only 50% of High Cell Cycle NIPUC, suggesting a prominent role of transcription factor activity in driving cell cycle dysregulation. Activity of the 121 transcription factors strongly associated with expression of EZH2 and other members of the PRC2 complex, suggesting regulation by this complex influences expression of the transcription factors in this network. Activity of transcription factors in this network also associated with signatures of pluripotency and epithelial-to-mesenchymal transition (EMT), suggesting they play a role in driving evolution to invasive carcinoma. Consistent with this, these transcription factors differed in activity between NIPUC and invasive urothelial carcinoma.

A Transcription Factor Signature Can Identify the CMS4 Subtype and Stratify the Prognostic Risk of Colorectal Cancer

Background

Colorectal cancer (CRC) is a heterogeneous disease, and current classification systems are insufficient for stratifying patients with different risks. This study aims to develop a generalized, individualized prognostic consensus molecular subtype (CMS)-transcription factors (TFs)-based signature that can predict the prognosis of CRC.

Methods

We obtained differentially expressed TF signature and target genes between the CMS4 and other CMS subtypes of CRC from The Cancer Genome Atlas (TCGA) database. A multi-dimensional network inference integrative analysis was conducted to identify the master genes and establish a CMS4-TFs-based signature. For validation, an in-house clinical cohort (n = 351) and another independent public CRC cohort (n = 565) were applied. Gene set enrichment analysis (GSEA) and prediction of immune cell infiltration were performed to interpret the biological significance of the model.

Results

A CMS4-TFs-based signature termed TF-9 that includes nine TF master genes was developed. Patients in the TF-9 high-risk group have significantly worse survival, regardless of clinical characteristics. The TF-9 achieved the highest mean C-index (0.65) compared to all other signatures reported (0.51 to 0.57). Immune infiltration revealed that the microenvironment in the high-risk group was highly immune suppressed, as evidenced by the overexpression of TIM3, CD39, and CD40, suggesting that high-risk patients may not directly benefit from the immune checkpoint inhibitors.

Conclusions

The TF-9 signature allows a more precise categorization of patients with relevant clinical and biological implications, which may be a valuable tool for improving the tailoring of therapeutic interventions in CRC patients.

Identification and Verification of m7G Modification Patterns and Characterization of Tumor Microenvironment Infiltration via Multi-Omics Analysis in Clear Cell Renal Cell Carcinoma

The epigenetic modification of tumorigenesis and progression in neoplasm has been demonstrated in recent studies. Nevertheless, the underlying association of N7-methylguanosine (m⁷G) regulation with molecular heterogeneity and tumor microenvironment (TME) in clear cell renal cell carcinoma (ccRCC) remains unknown. We explored the expression profiles and genetic variation features of m⁷G regulators and identified their correlations with patient outcomes in pan-cancer. Three distinct m⁷G modification patterns, including MGCS1, MGCS2, and MGCS3, were further determined and systematically characterized via multi-omics data in ccRCC. Compared with the other two subtypes, patients in MGCS3 exhibited a lower clinical stage/grade and better prognosis. MGCS1 showed the lowest enrichment of metabolic activities. MGCS2 was characterized by the suppression of immunity. We then established and validated a scoring tool named m7Sig, which could predict the prognosis of ccRCC patients. This study revealed that m⁷G modification played a vital role in the formation of the tumor microenvironment in ccRCC. Evaluating the m⁷G modification landscape helps us to raise awareness and strengthen the understanding of ccRCC’s characterization and, furthermore, to guide future clinical decision making.

Stage-stratified molecular profiling of non-muscle-invasive bladder cancer enhances biological, clinical, and therapeutic insight

Understanding the molecular determinants that underpin the clinical heterogeneity of non-muscle-invasive bladder cancer (NMIBC) is essential for prognostication and therapy development. Stage T1 disease in particular presents a high risk of progression and requires improved understanding. We present a detailed multi-omics study containing gene expression, copy number, and mutational profiles that show relationships to immune infiltration, disease recurrence, and progression to muscle invasion. We compare expression and genomic subtypes derived from all NMIBCs with those derived from the individual disease stages Ta and T1. We show that sufficient molecular heterogeneity exists within the separate stages to allow subclassification and that this is more clinically meaningful for stage T1 disease than that derived from all NMIBCs. This provides improved biological understanding and identifies subtypes of T1 tumors that may benefit from chemo- or immunotherapy.