KnockTF: a comprehensive human gene expression profile database with knockdown/knockout of transcription factors

TFTG: A comprehensive database for human transcription factors and their targets

Transcription factors (TFs) are major contributors to gene transcription, especially in controlling cell-specific gene expression and disease occurrence and development. Uncovering the relationship between TFs and their target genes is critical to understanding the mechanism of action of TFs. With the development of high-throughput sequencing techniques, a large amount of TF-related data has accumulated, which can be used to identify their target genes. In this study, we developed TFTG (Transcription Factor and Target Genes) database (http://tf.liclab.net/TFTG), which aimed to provide a large number of available human TF-target gene resources by multiple strategies, besides performing a comprehensive functional and epigenetic annotations and regulatory analyses of TFs. We identified extensive available TF-target genes by collecting and processing TF-associated ChIP-seq datasets, perturbation RNA-seq datasets and motifs. We also obtained experimentally confirmed relationships between TF and target genes from available resources. Overall, the target genes of TFs were obtained through integrating the relevant data of various TFs as well as fourteen identification strategies. Meanwhile, TFTG was embedded with user-friendly search, analysis, browsing, downloading and visualization functions. TFTG is designed to be a convenient resource for exploring human TF-target gene regulations, which will be useful for most users in the TF and gene expression regulation research.

PerturbDB for unraveling gene functions and regulatory networks

Perturb-Seq combines CRISPR (clustered regularly interspaced short palindromic repeats)-based genetic screens with single-cell RNA sequencing readouts for high-content phenotypic screens. Despite the rapid accumulation of Perturb-Seq datasets, there remains a lack of a user-friendly platform for their efficient reuse. Here, we developed PerturbDB (http://research.gzsys.org.cn/perturbdb), a platform to help users unveil gene functions using Perturb-Seq datasets. PerturbDB hosts 66 Perturb-Seq datasets, which encompass 4 518 521 single-cell transcriptomes derived from the knockdown of 10 194 genes across 19 different cell lines. All datasets were uniformly processed using the Mixscape algorithm. Genes were clustered by their perturbed transcriptomic phenotypes derived from Perturb-Seq data, resulting in 421 gene clusters, 157 of which were stable across different cellular contexts. Through integrating chemically perturbed transcriptomes with Perturb-Seq data, we identified 552 potential inhibitors targeting 1409 genes, including an mammalian target of rapamycin (mTOR) signaling inhibitor, retinol, which was experimentally verified. Moreover, we developed a 'Cancer' module to facilitate the understanding of the regulatory role of genes in cancer using Perturb-Seq data. An interactive web interface has also been developed, enabling users to visualize, analyze and download all the comprehensive datasets available in PerturbDB. PerturbDB will greatly drive gene functional studies and enhance our understanding of the regulatory roles of genes in diseases such as cancer.

NFAT activation by FKBP52 promotes cancer cell proliferation by suppressing p53

FK506-binding protein 52 (FKBP52) is a member of the FKBP family of proline isomerases. FKBP52 is up-regulated in various cancers and functions as a positive regulator of steroid hormone receptors. Depletion of FKBP52 is known to inhibit cell proliferation; however, the detailed mechanism remains poorly understood. In this study, we found that FKBP52 depletion decreased MDM2 transcription, leading to stabilization of p53, and suppressed cell proliferation. We identified NFATc1 and NFATc3 as transcription factors that regulate MDM2 We also found that FKBP52 associated with NFATc3 and facilitated its nuclear translocation. In addition, calcineurin, a well-known Ca2+ phosphatase essential for activation of NFAT, plays a role in MDM2 transcription. Supporting this notion, MDM2 expression was found to be regulated by intracellular Ca2+ Taken together, these findings reveal a new role of FKBP52 in promoting cell proliferation via the NFAT-MDM2-p53 axis, and indicate that inhibition of FKBP52 could be a new therapeutic tool to activate p53 and inhibit cell proliferation.

HAS2 facilitates glioma cell malignancy and suppresses ferroptosis in an FZD7-dependent manner

Glioma is the most common malignant tumor in the central nervous system, and it is crucial to uncover the factors that influence prognosis. In this study, we utilized Mfuzz to identify a gene set that showed a negative correlation with overall survival in patients with glioma. Gene Ontology (GO) enrichment analyses were then undertaken to gain insights into the functional characteristics and pathways associated with these genes. The expression distribution of Hyaluronan Synthase 2 (HAS2) was explored across multiple datasets, revealing its expression patterns. In vitro and in vivo experiments were carried out through gene knockdown and overexpression to validate the functionality of HAS2. Potential upstream transcription factors of HAS2 were predicted using transcriptional regulatory databases, and these predictions were experimentally validated using ChIP-PCR and dual-luciferase reporter gene assays. The results showed that elevated expression of HAS2 in glioma indicates poor prognosis. HAS2 was found to play a role in activating an antiferroptosis pathway in glioma cells. Inhibiting HAS2 significantly increased cellular sensitivity to ferroptosis-inducing agents. Finally, we determined that the oncogenic effect of HAS2 is mediated by the key receptor of the WNT pathway, FZD7.

Super-enhancer omics in stem cell

The hallmarks of stem cells, such as proliferation, self-renewal, development, differentiation, and regeneration, are critical to maintain stem cell identity which is sustained by genetic and epigenetic factors. Super-enhancers (SEs), which consist of clusters of active enhancers, play a central role in maintaining stemness hallmarks by specifically transcriptional model. The SE-navigated transcriptional complex, including SEs, non-coding RNAs, master transcriptional factors, Mediators and other co-activators, forms phase-separated condensates, which offers a toggle for directing diverse stem cell fate. With the burgeoning technologies of multiple-omics applied to examine different aspects of SE, we firstly raise the concept of "super-enhancer omics", inextricably linking to Pan-omics. In the review, we discuss the spatiotemporal organization and concepts of SEs, and describe links between SE-navigated transcriptional complex and stem cell features, such as stem cell identity, self-renewal, pluripotency, differentiation and development. We also elucidate the mechanism of stemness and oncogenic SEs modulating cancer stem cells via genomic and epigenetic alterations hijack in cancer stem cell. Additionally, we discuss the potential of targeting components of the SE complex using small molecule compounds, genome editing, and antisense oligonucleotides to treat SE-associated organ dysfunction and diseases, including cancer. This review also provides insights into the future of stem cell research through the paradigm of SEs.

Assessing next-generation sequencing-based computational methods for predicting transcriptional regulators with query gene sets

This article provides an in-depth review of computational methods for predicting transcriptional regulators (TRs) with query gene sets. Identification of TRs is of utmost importance in many biological applications, including but not limited to elucidating biological development mechanisms, identifying key disease genes, and predicting therapeutic targets. Various computational methods based on next-generation sequencing (NGS) data have been developed in the past decade, yet no systematic evaluation of NGS-based methods has been offered. We classified these methods into two categories based on shared characteristics, namely library-based and region-based methods. We further conducted benchmark studies to evaluate the accuracy, sensitivity, coverage, and usability of NGS-based methods with molecular experimental datasets. Results show that BART, ChIP-Atlas, and Lisa have relatively better performance. Besides, we point out the limitations of NGS-based methods and explore potential directions for further improvement.

TFTF: An R-Based Integrative Tool for Decoding Human Transcription Factor-Target Interactions

Transcription factors (TFs) are crucial in modulating gene expression and sculpting cellular and organismal phenotypes. The identification of TF-target gene interactions is pivotal for comprehending molecular pathways and disease etiologies but has been hindered by the demanding nature of traditional experimental approaches. This paper introduces a novel web application and package utilizing the R program, which predicts TF-target gene relationships and vice versa. Our application integrates the predictive power of various bioinformatic tools, leveraging their combined strengths to provide robust predictions. It merges databases for enhanced precision, incorporates gene expression correlation for accuracy, and employs pan-tissue correlation analysis for context-specific insights. The application also enables the integration of user data with established resources to analyze TF-target gene networks. Despite its current limitation to human data, it provides a platform to explore gene regulatory mechanisms comprehensively. This integrated, systematic approach offers researchers an invaluable tool for dissecting the complexities of gene regulation, with the potential for future expansions to include a broader range of species.

Mechanistic analysis of enhancer sequences in the estrogen receptor transcriptional program

Estrogen Receptor α (ERα) is a major lineage determining transcription factor (TF) in mammary gland development. Dysregulation of ERα-mediated transcriptional program results in cancer. Transcriptomic and epigenomic profiling of breast cancer cell lines has revealed large numbers of enhancers involved in this regulatory program, but how these enhancers encode function in their sequence remains poorly understood. A subset of ERα-bound enhancers are transcribed into short bidirectional RNA (enhancer RNA or eRNA), and this property is believed to be a reliable marker of active enhancers. We therefore analyze thousands of ERα-bound enhancers and build quantitative, mechanism-aware models to discriminate eRNAs from non-transcribing enhancers based on their sequence. Our thermodynamics-based models provide insights into the roles of specific TFs in ERα-mediated transcriptional program, many of which are supported by the literature. We use in silico perturbations to predict TF-enhancer regulatory relationships and integrate these findings with experimentally determined enhancer-promoter interactions to construct a gene regulatory network. We also demonstrate that the model can prioritize breast cancer-related sequence variants while providing mechanistic explanations for their function. Finally, we experimentally validate the model-proposed mechanisms underlying three such variants.

Ferroptosis inhibitors reduce celastrol toxicity and preserve its insulin sensitizing effects in insulin resistant HepG2 cells

OBJECTIVE

Research has shown that celastrol can effectively treat a variety of diseases, yet when passing a certain dosage threshold, celastrol becomes toxic, causing complications such as liver and kidney damage and erythrocytopenia, among others. With this dichotomy in mind, it is extremely important to find ways to preserve celastrol's efficacy while reducing or preventing its toxicity.

METHODS

In this study, insulin-resistant HepG2 (IR-HepG2) cells were prepared using palmitic acid and used for in vitro experiments. IR-HepG2 cells were treated with celastrol alone or in combination with N-acetylcysteine (NAC) or ferrostatin-1 (Fer-1) for 12, 24 or 48 h, at a range of doses. Cell counting kit-8 assay, Western blotting, quantitative reverse transcription-polymerase chain reaction, glucose consumption assessment, and flow cytometry were performed to measure celastrol's cytotoxicity and whether the cell death was linked to ferroptosis.

RESULTS

Celastrol treatment increased lipid oxidation and decreased expression of anti-ferroptosis proteins in IR-HepG2 cells. Celastrol downregulated glutathione peroxidase 4 (GPX4) mRNA. Molecular docking models predicted that solute carrier family 7 member 11 (SLC7A11) and GPX4 were covalently bound by celastrol. Importantly, we found for the first time that the application of ferroptosis inhibitors (especially NAC) was able to reduce celastrol's toxicity while preserving its ability to improve insulin sensitivity in IR-HepG2 cells.

CONCLUSION

One potential mechanism of celastrol's cytotoxicity is the induction of ferroptosis, which can be alleviated by treatment with ferroptosis inhibitors. These findings provide a new strategy to block celastrol's toxicity while preserving its therapeutic effects. Please cite this article as: Liu JJ, Zhang X, Qi MM, Chi YB, Cai BL, Peng B, Zhang DH. Ferroptosis inhibitors reduce celastrol toxicity and preserve its insulin sensitizing effects in insulin resistant HepG2 cells. J Integr Med. 2024; 22(3): 286-294.

Inferring gene regulatory networks from single-cell multiome data using atlas-scale external data

Existing methods for gene regulatory network (GRN) inference rely on gene expression data alone or on lower resolution bulk data. Despite the recent integration of chromatin accessibility and RNA sequencing data, learning complex mechanisms from limited independent data points still presents a daunting challenge. Here we present LINGER (Lifelong neural network for gene regulation), a machine-learning method to infer GRNs from single-cell paired gene expression and chromatin accessibility data. LINGER incorporates atlas-scale external bulk data across diverse cellular contexts and prior knowledge of transcription factor motifs as a manifold regularization. LINGER achieves a fourfold to sevenfold relative increase in accuracy over existing methods and reveals a complex regulatory landscape of genome-wide association studies, enabling enhanced interpretation of disease-associated variants and genes. Following the GRN inference from reference single-cell multiome data, LINGER enables the estimation of transcription factor activity solely from bulk or single-cell gene expression data, leveraging the abundance of available gene expression data to identify driver regulators from case-control studies.

Circular RNA CircSLC22A23 Promotes Gastric Cancer Progression by Activating HNRNPU Expression

BACKGROUND

Circular RNAs (CircRNAs) play essential roles in cancer occurrence as regulatory RNAs. However, circRNA-mediated regulation of gastric cancer (GC) remains poorly understood.

AIM

The purpose of this study was to investigate the molecular mechanism of circSLC22A23 (hsa_circ_0075504) underlying GC occurrence.

METHODS

CircSLC22A23 levels were first quantified by quantitative real-time reverse transcription-polymerase chain reaction in GC cell lines, 80 paired GC tissues and adjacent normal tissues, and 27 pairs of plasma samples from preoperative and postoperative patients with GC. Then circSLC22A23 was knocked-down with short hairpin RNA to analyze its oncogenic effects on the proliferation, migration, and invasion of GC cells. Finally, circRNA-binding proteins and their downstream target genes were identified by RNA pulldown, mass spectrometry, RNA immunoprecipitation, quantitative real-time reverse transcription-polymerase chain reaction, and Western blot assays.

RESULTS

CircSLC22A23 was found to be highly expressed in GC cells, GC tissues, and plasma from GC patients. Knockdown of circSLC22A23 inhibited GC cell proliferation, migration and invasion. RNA pulldown and RNA immunoprecipitation assays verified the interaction between circSLC22A23 and heterogeneous nuclear ribonucleoprotein U (HNRNPU). Knockdown of circSLC22A23 decreased HNRNPU protein levels. Moreover, rescue assays showed that the tumor suppressive effect of circSLC22A23 knockdown was reversed by HNRNPU overexpression. Finally, epidermal growth factor receptor (EGFR) was found to be one of the downstream target genes of HNRNPU that was up regulated by circSLC22A23.

CONCLUSION

CircSLC22A23 regulated the transcription of EGFR through activation of HNRNPU in GC cells, suggesting that circSLC22A23 may serve as a potential therapeutic target for the treatment of GC.

Assessing NGS-based computational methods for predicting transcriptional regulators with query gene sets

This article provides an in-depth review of computational methods for predicting transcriptional regulators with query gene sets. Identification of transcriptional regulators is of utmost importance in many biological applications, including but not limited to elucidating biological development mechanisms, identifying key disease genes, and predicting therapeutic targets. Various computational methods based on next-generation sequencing (NGS) data have been developed in the past decade, yet no systematic evaluation of NGS-based methods has been offered. We classified these methods into two categories based on shared characteristics, namely library-based and region-based methods. We further conducted benchmark studies to evaluate the accuracy, sensitivity, coverage, and usability of NGS-based methods with molecular experimental datasets. Results show that BART, ChIP-Atlas, and Lisa have relatively better performance. Besides, we point out the limitations of NGS-based methods and explore potential directions for further improvement.

Key points

An introduction to available computational methods for predicting functional TRs from a query gene set.A detailed walk-through along with practical concerns and limitations.A systematic benchmark of NGS-based methods in terms of accuracy, sensitivity, coverage, and usability, using 570 TR perturbation-derived gene sets.NGS-based methods outperform motif-based methods. Among NGS methods, those utilizing larger databases and adopting region-centric approaches demonstrate favorable performance. BART, ChIP-Atlas, and Lisa are recommended as these methods have overall better performance in evaluated scenarios.

TSCRE: a comprehensive database for tumor-specific cis-regulatory elements

Cis-regulatory elements (CREs) and super cis-regulatory elements (SCREs) are non-coding DNA regions which influence the transcription of nearby genes and play critical roles in development. Dysregulated CRE and SCRE activities have been reported to alter the expression of oncogenes and tumor suppressors, thereby regulating cancer hallmarks. To address the strong need for a comprehensive catalogue of dysregulated CREs and SCREs in human cancers, we present TSCRE (http://tscre.zsqylab.com/), an open resource providing tumor-specific and cell type-specific CREs and SCREs derived from the re-analysis of publicly available histone modification profiles. Currently, TSCRE contains 1 864 941 dysregulated CREs and 68 253 dysregulated SCREs identified from 1366 human patient samples spanning 17 different cancer types and 9 histone marks. Over 95% of these elements have been validated in public resources. TSCRE offers comprehensive annotations for each element, including associated genes, expression patterns, clinical prognosis, somatic mutations, transcript factor binding sites, cancer-type specificity, and drug response. Additionally, TSCRE integrates pathway and transcript factor enrichment analyses for each study, enabling in-depth functional and mechanistic investigations. Furthermore, TSCRE provides an interactive interface for users to explore any CRE and SCRE of interest. We believe TSCRE will be a highly valuable platform for the community to discover candidate cancer biomarkers.

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.

KnockTF 2.0: a comprehensive gene expression profile database with knockdown/knockout of transcription (co-)factors in multiple species

Transcription factors (TFs), transcription co-factors (TcoFs) and their target genes perform essential functions in diseases and biological processes. KnockTF 2.0 (http://www.licpathway.net/KnockTF/index.html) aims to provide comprehensive gene expression profile datasets before/after T(co)F knockdown/knockout across multiple tissue/cell types of different species. Compared with KnockTF 1.0, KnockTF 2.0 has the following improvements: (i) Newly added T(co)F knockdown/knockout datasets in mice, Arabidopsis thaliana and Zea mays and also an expanded scale of datasets in humans. Currently, KnockTF 2.0 stores 1468 manually curated RNA-seq and microarray datasets associated with 612 TFs and 172 TcoFs disrupted by different knockdown/knockout techniques, which are 2.5 times larger than those of KnockTF 1.0. (ii) Newly added (epi)genetic annotations for T(co)F target genes in humans and mice, such as super-enhancers, common SNPs, methylation sites and chromatin interactions. (iii) Newly embedded and updated search and analysis tools, including T(co)F Enrichment (GSEA), Pathway Downstream Analysis and Search by Target Gene (BLAST). KnockTF 2.0 is a comprehensive update of KnockTF 1.0, which provides more T(co)F knockdown/knockout datasets and (epi)genetic annotations across multiple species than KnockTF 1.0. KnockTF 2.0 facilitates not only the identification of functional T(co)Fs and target genes but also the investigation of their roles in the physiological and pathological processes.

scGRN: a comprehensive single-cell gene regulatory network platform of human and mouse

Gene regulatory networks (GRNs) are interpretable graph models encompassing the regulatory interactions between transcription factors (TFs) and their downstream target genes. Making sense of the topology and dynamics of GRNs is fundamental to interpreting the mechanisms of disease etiology and translating corresponding findings into novel therapies. Recent advances in single-cell multi-omics techniques have prompted the computational inference of GRNs from single-cell transcriptomic and epigenomic data at an unprecedented resolution. Here, we present scGRN (https://bio.liclab.net/scGRN/), a comprehensive single-cell multi-omics gene regulatory network platform of human and mouse. The current version of scGRN catalogs 237 051 cell type-specific GRNs (62 999 692 TF-target gene pairs), covering 160 tissues/cell lines and 1324 single-cell samples. scGRN is the first resource documenting large-scale cell type-specific GRN information of diverse human and mouse conditions inferred from single-cell multi-omics data. We have implemented multiple online tools for effective GRN analysis, including differential TF-target network analysis, TF enrichment analysis, and pathway downstream analysis. We also provided details about TF binding to promoters, super-enhancers and typical enhancers of target genes in GRNs. Taken together, scGRN is an integrative and useful platform for searching, browsing, analyzing, visualizing and downloading GRNs of interest, enabling insight into the differences in regulatory mechanisms across diverse conditions.

eRNAbase: a comprehensive database for decoding the regulatory eRNAs in human and mouse

Enhancer RNAs (eRNAs) transcribed from distal active enhancers serve as key regulators in gene transcriptional regulation. The accumulation of eRNAs from multiple sequencing assays has led to an urgent need to comprehensively collect and process these data to illustrate the regulatory landscape of eRNAs. To address this need, we developed the eRNAbase (http://bio.liclab.net/eRNAbase/index.php) to store the massive available resources of human and mouse eRNAs and provide comprehensive annotation and analyses for eRNAs. The current version of eRNAbase cataloged 10 399 928 eRNAs from 1012 samples, including 858 human samples and 154 mouse samples. These eRNAs were first identified and uniformly processed from 14 eRNA-related experiment types manually collected from GEO/SRA and ENCODE. Importantly, the eRNAbase provides detailed and abundant (epi)genetic annotations in eRNA regions, such as super enhancers, enhancers, common single nucleotide polymorphisms, expression quantitative trait loci, transcription factor binding sites, CRISPR/Cas9 target sites, DNase I hypersensitivity sites, chromatin accessibility regions, methylation sites, chromatin interactions regions, topologically associating domains and RNA spatial interactions. Furthermore, the eRNAbase provides users with three novel analyses including eRNA-mediated pathway regulatory analysis, eRNA-based variation interpretation analysis and eRNA-mediated TF-target gene analysis. Hence, eRNAbase is a powerful platform to query, browse and visualize regulatory cues associated with eRNAs.

LncSEA 2.0: an updated platform for long non-coding RNA related sets and enrichment analysis

Long non-coding RNAs (lncRNAs) possess a wide range of biological functions, and research has demonstrated their significance in regulating major biological processes such as development, differentiation, and immune response. The accelerating accumulation of lncRNA research has greatly expanded our understanding of lncRNA functions. Here, we introduce LncSEA 2.0 (http://bio.liclab.net/LncSEA/index.php), aiming to provide a more comprehensive set of functional lncRNAs and enhanced enrichment analysis capabilities. Compared with LncSEA 1.0, we have made the following improvements: (i) We updated the lncRNA sets for 11 categories and extremely expanded the lncRNA scopes for each set. (ii) We newly introduced 15 functional lncRNA categories from multiple resources. This update not only included a significant amount of downstream regulatory data for lncRNAs, but also covered numerous epigenetic regulatory data sets, including lncRNA-related transcription co-factor binding, chromatin regulator binding, and chromatin interaction data. (iii) We incorporated two new lncRNA set enrichment analysis functions based on GSEA and GSVA. (iv) We adopted the snakemake analysis pipeline to track data processing and analysis. In summary, LncSEA 2.0 offers a more comprehensive collection of lncRNA sets and a greater variety of enrichment analysis modules, assisting researchers in a more comprehensive study of the functional mechanisms of lncRNAs.

The feedback loop of EFTUD2/c-MYC impedes chemotherapeutic efficacy by enhancing EFTUD2 transcription and stabilizing c-MYC protein in colorectal cancer

BACKGROUND

Chemoresistance presents a significant obstacle in the treatment of colorectal cancer (CRC), yet the molecular basis underlying CRC chemoresistance remains poorly understood, impeding the development of new therapeutic interventions. Elongation factor Tu GTP binding domain containing 2 (EFTUD2) has emerged as a potential oncogenic factor implicated in various cancer types, where it fosters tumor growth and survival. However, its specific role in modulating the sensitivity of CRC cells to chemotherapy is still unclear.

METHODS

Public dataset analysis and in-house sample validation were conducted to assess the expression of EFTUD2 in 5-fluorouracil (5-FU) chemotherapy-resistant CRC cells and the potential of EFTUD2 as a prognostic indicator for CRC. Experiments both in vitro, including MTT assay, EdU cell proliferation assay, TUNEL assay, and clone formation assay and in vivo, using cell-derived xenograft models, were performed to elucidate the function of EFTUD2 in sensitivity of CRC cells to 5-FU treatment. The molecular mechanism on the reciprocal regulation between EFTUD2 and the oncogenic transcription factor c-MYC was investigated through molecular docking, ubiquitination assay, chromatin immunoprecipitation (ChIP), dual luciferase reporter assay, and co-immunoprecipitation (Co-IP).

RESULTS

We found that EFTUD2 expression was positively correlated with 5-FU resistance, higher pathological grade, and poor prognosis in CRC patients. We also demonstrated both in vitro and in vivo that knockdown of EFTUD2 sensitized CRC cells to 5-FU treatment, whereas overexpression of EFTUD2 impaired such sensitivity. Mechanistically, we uncovered that EFTUD2 physically interacted with and stabilized c-MYC protein by preventing its ubiquitin-mediated proteasomal degradation. Intriguingly, we found that c-MYC directly bound to the promoter region of EFTUD2 gene, activating its transcription. Leveraging rescue experiments, we further confirmed that the effect of EFTUD2 on 5-FU resistance was dependent on c-MYC stabilization.

CONCLUSION

Our findings revealed a positive feedback loop involving an EFTUD2/c-MYC axis that hampers the efficacy of 5-FU chemotherapy in CRC cells by increasing EFTUD2 transcription and stabilizing c-MYC oncoprotein. This study highlights the potential of EFTUD2 as a promising therapeutic target to surmount chemotherapy resistance in CRC patients.

Specifying cellular context of transcription factor regulons for exploring context-specific gene regulation programs

Understanding the role of transcription and transcription factors in cellular identity and disease, such as cancer and autoimmunity, is essential. However, comprehensive data resources for cell line-specific transcription factor-to-target gene annotations are currently limited. To address this, we developed a straightforward method to define regulons that capture the cell-specific aspects of TF binding and transcript expression levels. By integrating cellular transcriptome and transcription factor binding data, we generated regulons for four common cell lines comprising both proximal and distal cell line-specific regulatory events. Through systematic benchmarking involving transcription factor knockout experiments, we demonstrated performance on par with state-of-the-art methods, with our method being easily applicable to other cell types of interest. We present case studies using three cancer single-cell datasets to showcase the utility of these cell-type-specific regulons in exploring transcriptional dysregulation. In summary, this study provides a valuable tool and a resource for systematically exploring cell line-specific transcriptional regulations, emphasizing the utility of network analysis in deciphering disease mechanisms.

Decoding LINC00052 role in breast cancer by bioinformatic and experimental analyses

LncRNA is a group of transcripts with a length exceeding 200 nucleotides that contribute to tumour development. Our research group found that LINC00052 expression was repressed during the formation of breast cancer (BC) multicellular spheroids. Intriguingly, LINC00052 precise role in BC remains uncertain. We explored LINC00052 expression in BC patients` RNA samples (TCGA) in silico, as well as in an in-house patient cohort, and inferred its cellular and molecular mechanisms. In vitro studies evaluated LINC00052 relevance in BC cells viability, cell cycle and DNA damage. Results. Bioinformatic RNAseq analysis of BC patients showed that LINC00052 is overexpressed in samples from all BC molecular subtypes. A similar LINC00052 expression pattern was observed in an in-house patient cohort. In addition, higher LINC00052 levels are related to better BC patient´s overall survival. Remarkably, MCF-7 and ZR-75-1 cells treated with estradiol showed increased LINC00052 expression compared to control, while these changes were not observed in MDA-MB-231 cells. In parallel, bioinformatic analyses indicated that LINC00052 influences DNA damage and cell cycle. MCF-7 cells with low LINC00052 levels exhibited increased cellular protection against DNA damage and diminished growth capacity. Furthermore, in cisplatin-resistant MCF-7 cells, LINC00052 expression was downregulated. Conclusion. This work shows that LINC00052 expression is associated with better BC patient survival. Remarkably, LINC00052 expression can be regulated by Estradiol. Additionally, assays suggest that LINC00052 could modulate MCF-7 cells growth and DNA damage repair. Overall, this study highlights the need for further research to unravel LINC00052 molecular mechanisms and potential clinical applications in BC.

Identifying CD1c as a potential biomarker by the comprehensive exploration of tumor mutational burden and immune infiltration in diffuse large B cell lymphoma

Background

Tumor mutational burden (TMB) is a valuable prognostic biomarker. This study explored the predictive value of TMB and the potential association between TMB and immune infiltration in diffuse large B-cell lymphoma (DLBCL).

Methods

We downloaded the gene expression profile, somatic mutation, and clinical data of DLBCL patients from The Cancer Genome Atlas (TCGA) database. We classified the samples into high-and low-TMB groups to identify differentially expressed genes (DEGs). Functional enrichment analyses were performed to determine the biological functions of the DEGs. We utilized the cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) algorithm to estimate the abundance of 22 immune cells, and the significant difference was determined by the Wilcoxon rank-sum test between the high- and low-TMB group. Hub gene had been screened as the prognostic TMB-related immune biomarker by the combination of the Immunology Database and Analysis Portal (ImmPort) database and the univariate Cox analysis from the Gene Expression Omnibus (GEO) database including six DLBCL datasets. Various database applications such as Tumor Immune Estimation Resource (TIMER), CellMiner, konckTF, and Genotype-Tissue Expression (GTEx) verified the functions of the target gene. Wet assay confirmed the target gene expression at RNA and protein levels in DLBCL tissue and cell samples.

Results

Single nucleotide polymorphism (SNP) occurred more frequently than insertion and deletion, and C > T was the most common single nucleotide variant (SNV) in DLBCL. Survival analysis showed that the high-TMB group conferred poor survival outcomes. A total of 62 DEGs were obtained, and 13 TMB-related immune genes were identified. Univariate Cox analysis results illustrated that CD1c mutation was associated with lower TMB and manifested a satisfactory clinical prognosis by analysis of large samples from the GEO database. In addition, infiltration levels of immune cells in the high-TMB group were lower. Using the TIMER database, we systematically analyzed that the expression of CD1c was positively correlated with B cells, neutrophils, and dendritic cells and negatively correlated with CD8+ T cells, CD4+ T cells, and macrophages. Drug sensitivity showed a significant positive correlation between CD1c expression level and clinical drug sensitivity from the CellMiner database. CREB1, AHR, and TOX were used to comprehensively explore the regulation of CD1c-related transcription factors and signaling pathways by the KnockTF database. We searched the GETx database to compare the mRNA expression levels of CD1c between DLBCL and normal tissues, and the results suggested a significant difference between them. Moreover, wet experiments were conducted to verify the high expression of CD1c in DLBCL at the RNA and protein levels.

Conclusions

Higher TMB correlated with poor survival outcomes and inhibited the immune infiltrates in DLBCL. Our results suggest that CD1c is a TMB-related prognostic biomarker.

Computational tools for inferring transcription factor activity

Transcription factors (TFs) are essential players in orchestrating the regulatory landscape in cells. Still, their exact modes of action and dependencies on other regulatory aspects remain elusive. Since TFs act cell type-specific and each TF has its own characteristics, untangling their regulatory interactions from an experimental point of view is laborious and convoluted. Thus, there is an ongoing development of computational tools that estimate transcription factor activity (TFA) from a variety of data modalities, either based on a mapping of TFs to their putative target genes or in a genome-wide, gene-unspecific fashion. These tools can help to gain insights into TF regulation and to prioritize candidates for experimental validation. We want to give an overview of available computational tools that estimate TFA, illustrate examples of their application, debate common result validation strategies, and discuss assumptions and concomitant limitations.

Expanding the coverage of regulons from high-confidence prior knowledge for accurate estimation of transcription factor activities

Gene regulation plays a critical role in the cellular processes that underlie human health and disease. The regulatory relationship between transcription factors (TFs), key regulators of gene expression, and their target genes, the so called TF regulons, can be coupled with computational algorithms to estimate the activity of TFs. However, to interpret these findings accurately, regulons of high reliability and coverage are needed. In this study, we present and evaluate a collection of regulons created using the CollecTRI meta-resource containing signed TF-gene interactions for 1186 TFs. In this context, we introduce a workflow to integrate information from multiple resources and assign the sign of regulation to TF-gene interactions that could be applied to other comprehensive knowledge bases. We find that the signed CollecTRI-derived regulons outperform other public collections of regulatory interactions in accurately inferring changes in TF activities in perturbation experiments. Furthermore, we showcase the value of the regulons by examining TF activity profiles in three different cancer types and exploring TF activities at the level of single-cells. Overall, the CollecTRI-derived TF regulons enable the accurate and comprehensive estimation of TF activities and thereby help to interpret transcriptomics data.

Cis-Cardio: A comprehensive analysis platform for cardiovascular-relavant cis-regulation in human and mouse

Cis-regulatory elements are important molecular switches in controlling gene expression and are regarded as determinant hubs in the transcriptional regulatory network. Collection and processing of large-scale cis-regulatory data are urgent to decipher the potential mechanisms of cardiovascular diseases from a cis-regulatory element aspect. Here, we developed a novel web server, Cis-Cardio, which aims to document a large number of available cardiovascular-related cis-regulatory data and to provide analysis for unveiling the comprehensive mechanisms at a cis-regulation level. The current version of Cis-Cardio catalogs a total of 45,382,361 genomic regions from 1,013 human and mouse epigenetic datasets, including ATAC-seq, DNase-seq, Histone ChIP-seq, TF/TcoF ChIP-seq, RNA polymerase ChIP-seq, and Cohesin ChIP-seq. Importantly, Cis-Cardio provides six analysis tools, including region overlap analysis, element upstream/downstream analysis, transcription regulator enrichment analysis, variant interpretation, and protein-protein interaction-based co-regulatory analysis. Additionally, Cis-Cardio provides detailed and abundant (epi-) genetic annotations in cis-regulatory regions, such as super-enhancers, enhancers, transcription factor binding sites (TFBSs), methylation sites, common SNPs, risk SNPs, expression quantitative trait loci (eQTLs), motifs, DNase I hypersensitive sites (DHSs), and 3D chromatin interactions. In summary, Cis-Cardio is a valuable resource for elucidating and analyzing regulatory cues of cardiovascular-specific cis-regulatory elements. The platform is freely available at http://www.licpathway.net/Cis-Cardio/index.html.

The Landscape of m6A Regulators in Multiple Brain Regions of Alzheimer's Disease

Alzheimer's disease research has been conducted for many years, yet no effective cure methods have been found. N6-methyladenosine (m6A) RNA methylation, an essential post-transcriptional regulation mechanism, has been discovered to affect essential neurobiological processes, such as brain cell development and aging, which are closely related to neurodegenerative diseases such as Alzheimer's disease. The relationship between Alzheimer's disease and the m6A mechanism still needs further investigation. Our work evaluated the alteration profile of m6A regulators and their influences on Alzheimer's disease in 4 brain regions: the postcentral gyrus, superior frontal gyrus, hippocampus, and entorhinal cortex. We found that the expression levels of the m6A regulators FTO, ELAVL1, and YTHDF2 were altered in Alzheimer's disease and were related to pathological development and cognitive levels. We also assessed AD-related biological processes influenced by m6A regulators via GSEA and GSVA method. Biological Processes Gene Ontology terms including memory, cognition, and synapse-signaling were found to potentially be affected by m6A regulators in AD. We also found different m6A modification patterns in AD samples among different brain regions, mainly due to differences in m6A readers. Finally, we further evaluated the importance of AD-related regulators based on the WGCNA method, assessed their potential targets based on correlation relationships, and constructed diagnostic models in 3 of all 4 regions using hub regulators, including FTO, YTHDC1, YTHDC2, etc., and their potential targets. This work aims to provide a reference for the follow-up study of m6A and Alzheimer's disease.

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

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

Systematic discovery of transcription factors that improve hPSC-derived cardiomyocyte maturation via temporal analysis of bioengineered cardiac tissues

Human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs) have the potential to become powerful tools for disease modeling, drug testing, and transplantation; however, their immaturity limits their applications. Transcription factor (TF) overexpression can improve hPSC-CM maturity, but identifying these TFs has been elusive. Toward this, we establish here an experimental framework for systematic identification of maturation enhancing factors. Specifically, we performed temporal transcriptome RNAseq analyses of progressively matured hPSC-derived cardiomyocytes across 2D and 3D differentiation systems and further compared these bioengineered tissues to native fetal and adult-derived tissues. These analyses revealed 22 TFs whose expression did not increase in 2D differentiation systems but progressively increased in 3D culture systems and adult mature cell types. Individually overexpressing each of these TFs in immature hPSC-CMs identified five TFs (KLF15, ZBTB20, ESRRA, HOPX, and CAMTA2) as regulators of calcium handling, metabolic function, and hypertrophy. Notably, the combinatorial overexpression of KLF15, ESRRA, and HOPX improved all three maturation parameters simultaneously. Taken together, we introduce a new TF cocktail that can be used in solo or in conjunction with other strategies to improve hPSC-CM maturation and anticipate that our generalizable methodology can also be implemented to identify maturation-associated TFs for other stem cell progenies.

Characterizing the targets of transcription regulators by aggregating ChIP-seq and perturbation expression data sets

Mapping the gene targets of chromatin-associated transcription regulators (TRs) is a major goal of genomics research. ChIP-seq of TRs and experiments that perturb a TR and measure the differential abundance of gene transcripts are a primary means by which direct relationships are tested on a genomic scale. It has been reported that there is a poor overlap in the evidence across gene regulation strategies, emphasizing the need for integrating results from multiple experiments. Although research consortia interested in gene regulation have produced a valuable trove of high-quality data, there is an even greater volume of TR-specific data throughout the literature. In this study, we show a workflow for the identification, uniform processing, and aggregation of ChIP-seq and TR perturbation experiments for the ultimate purpose of ranking human and mouse TR-target interactions. Focusing on an initial set of eight regulators (ASCL1, HES1, MECP2, MEF2C, NEUROD1, PAX6, RUNX1, and TCF4), we identified 497 experiments suitable for analysis. We used this corpus to examine data concordance, to identify systematic patterns of the two data types, and to identify putative orthologous interactions between human and mouse. We build upon commonly used strategies to forward a procedure for aggregating and combining these two genomic methodologies, assessing these rankings against independent literature-curated evidence. Beyond a framework extensible to other TRs, our work also provides empirically ranked TR-target listings, as well as transparent experiment-level gene summaries for community use.

Maximizing the utility of public data

The human genome project galvanized the scientific community around an ambitious goal. Upon completion, the project delivered several discoveries, and a new era of research commenced. More importantly, novel technologies and analysis methods materialized during the project period. The cost reduction allowed many more labs to generate high-throughput datasets. The project also served as a model for other extensive collaborations that generated large datasets. These datasets were made public and continue to accumulate in repositories. As a result, the scientific community should consider how these data can be utilized effectively for the purposes of research and the public good. A dataset can be re-analyzed, curated, or integrated with other forms of data to enhance its utility. We highlight three important areas to achieve this goal in this brief perspective. We also emphasize the critical requirements for these strategies to be successful. We draw on our own experience and others in using publicly available datasets to support, develop, and extend our research interest. Finally, we underline the beneficiaries and discuss some risks involved in data reuse.

Integrative bioinformatics analysis to identify novel biomarkers associated with non-obstructive azoospermia

Aim

This study aimed to identify autophagy-related genes (ARGs) associated with non-obstructive azoospermia and explore the underlying molecular mechanisms.

Methods

Two datasets associated with azoospermia were downloaded from the Gene Expression Omnibus database, and ARGs were obtained from the Human Autophagy-dedicated Database. Autophagy-related differentially expressed genes were identified in the azoospermia and control groups. These genes were subjected to Gene Ontology and Kyoto Encyclopedia of Genes and Genomes, protein–protein interaction (PPI) network, and functional similarity analyses. After identifying the hub genes, immune infiltration and hub gene–RNA-binding protein (RBP)–transcription factor (TF)–miRNA–drug interactions were analyzed.

Results

A total 46 differentially expressed ARGs were identified between the azoospermia and control groups. These genes were enriched in autophagy-associated functions and pathways. Eight hub genes were selected from the PPI network. Functional similarity analysis revealed that HSPA5 may play a key role in azoospermia. Immune cell infiltration analysis revealed that activated dendritic cells were significantly decreased in the azoospermia group compared to those in the control groups. Hub genes, especially ATG3, KIAA0652, MAPK1, and EGFR were strongly correlated with immune cell infiltration. Finally, a hub gene–miRNA–TF–RBP–drug network was constructed.

Conclusion

The eight hub genes, including EGFR, HSPA5, ATG3, KIAA0652, and MAPK1, may serve as biomarkers for the diagnosis and treatment of azoospermia. The study findings suggest potential targets and mechanisms for the occurrence and development of this disease.

Discovery of Targets for Immune-Metabolic Antitumor Drugs Identifies Estrogen-Related Receptor Alpha

Drugs that kill tumors through multiple mechanisms have the potential for broad clinical benefits. Here, we first developed an in silico multiomics approach (BipotentR) to find cancer cell-specific regulators that simultaneously modulate tumor immunity and another oncogenic pathway and then used it to identify 38 candidate immune-metabolic regulators. We show the tumor activities of these regulators stratify patients with melanoma by their response to anti-PD-1 using machine learning and deep neural approaches, which improve the predictive power of current biomarkers. The topmost identified regulator, ESRRA, is activated in immunotherapy-resistant tumors. Its inhibition killed tumors by suppressing energy metabolism and activating two immune mechanisms: (i) cytokine induction, causing proinflammatory macrophage polarization, and (ii) antigen-presentation stimulation, recruiting CD8+ T cells into tumors. We also demonstrate a wide utility of BipotentR by applying it to angiogenesis and growth suppressor evasion pathways. BipotentR (http://bipotentr.dfci.harvard.edu/) provides a resource for evaluating patient response and discovering drug targets that act simultaneously through multiple mechanisms.

SIGNIFICANCE

BipotentR presents resources for evaluating patient response and identifying targets for drugs that can kill tumors through multiple mechanisms concurrently. Inhibition of the topmost candidate target killed tumors by suppressing energy metabolism and effects on two immune mechanisms. This article is highlighted in the In This Issue feature, p. 517.

High dose isoleucine stabilizes nuclear PTEN to suppress the proliferation of lung cancer

PURPOSE

Cancer cells require a supply of amino acids, particularly essential amino acids such as branched-chain amino acids (BCAAs, i.e., valine, leucine, and isoleucine), to meet the increased nutrient demands of malignant tumors. The cell-autonomous and non-autonomous roles of altered BCAA supply have been implicated in cancer progression. The critical proteins involved in BCAA uptake, transport, metabolism, etc. serve as potential therapeutic biomarkers in human cancers. Here, we summarize the potential anti-tumor mechanism of BCAA by exploring the chain reaction triggered by increased BCAA supply in the tumor.

METHOD

A system-wide strategy was employed to provide a generic solution to establish the links between BCAA and cancer based on comprehensive omics, molecular experimentation, and data analysis.

RESULTS

BCAA over-supplementation (900 mg/kg) significantly inhibited tumor growth and reduced tumor burden, with isoleucine having the most pronounced effect. Surprisingly, isoleucine inhibited tumor growth independently of mTORC1 activation, a classical amino acid sensor. Exploratory transcriptome analysis revealed that Phosphatase and tensin homolog (PTEN) is the critical factor in the anti-tumor effect of isoleucine. By inhibiting PTEN ubiquitination, isoleucine can promote PTEN nuclear import and maintain PTEN nuclear stability. Interestingly, this process was regulated by isoleucine-tRNA ligase, cytoplasmic (IARS), a direct target of isoleucine. We demonstrated the enhanced interaction between IARS and PTEN in the presence of excess isoleucine. At the same time, IARS knockout leads to loss of isoleucine tumor suppressor ability.

CONCLUSION

Overall, our results provide insights into the regulation of the IARS-PTEN anti-tumor axis by isoleucine and reveal a unique therapeutic approach based on enhancing cellular isoleucine supply.

Superenhancers as master gene regulators and novel therapeutic targets in brain tumors

Transcriptional deregulation, a cancer cell hallmark, is driven by epigenetic abnormalities in the majority of brain tumors, including adult glioblastoma and pediatric brain tumors. Epigenetic abnormalities can activate epigenetic regulatory elements to regulate the expression of oncogenes. Superenhancers (SEs), identified as novel epigenetic regulatory elements, are clusters of enhancers with cell-type specificity that can drive the aberrant transcription of oncogenes and promote tumor initiation and progression. As gene regulators, SEs are involved in tumorigenesis in a variety of tumors, including brain tumors. SEs are susceptible to inhibition by their key components, such as bromodomain protein 4 and cyclin-dependent kinase 7, providing new opportunities for antitumor therapy. In this review, we summarized the characteristics and identification, unique organizational structures, and activation mechanisms of SEs in tumors, as well as the clinical applications related to SEs in tumor therapy and prognostication. Based on a review of the literature, we discussed the relationship between SEs and different brain tumors and potential therapeutic targets, focusing on glioblastoma.

Localization of Noise in Biochemical Networks

Noise, or uncertainty in biochemical networks, has become an important aspect of many biological problems. Noise can arise and propagate from external factors and probabilistic chemical reactions occurring in small cellular compartments. For species survival, it is important to regulate such uncertainties in executing vital cell functions. Regulated noise can improve adaptability, whereas uncontrolled noise can cause diseases. Simulation can provide a detailed analysis of uncertainties, but parameters such as rate constants and initial conditions are usually unknown. A general understanding of noise dynamics from the perspective of network structure is highly desirable. In this study, we extended the previously developed law of localization for characterizing noise in terms of (co)variances and developed noise localization theory. With linear noise approximation, we can expand a biochemical network into an extended set of differential equations representing a fictitious network for pseudo-components consisting of variances and covariances, together with chemical species. Through localization analysis, perturbation responses at the steady state of pseudo-components can be summarized into a sensitivity matrix that only requires knowledge of network topology. Our work allows identification of buffering structures at the level of species, variances, and covariances and can provide insights into noise flow under non-steady-state conditions in the form of a pseudo-chemical reaction. We tested noise localization in various systems, and here we discuss its implications and potential applications. Results show that this theory is potentially applicable in discriminating models, scanning network topologies with interesting noise behavior, and designing and perturbing networks with the desired response.

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.

SEdb 2.0: a comprehensive super-enhancer database of human and mouse

Super-enhancers (SEs) are cell-specific DNA cis-regulatory elements that can supervise the transcriptional regulation processes of downstream genes. SEdb 2.0 (http://www.licpathway.net/sedb) aims to provide a comprehensive SE resource and annotate their potential roles in gene transcriptions. Compared with SEdb 1.0, we have made the following improvements: (i) Newly added the mouse SEs and expanded the scale of human SEs. SEdb 2.0 contained 1 167 518 SEs from 1739 human H3K27ac chromatin immunoprecipitation sequencing (ChIP-seq) samples and 550 226 SEs from 931 mouse H3K27ac ChIP-seq samples, which was five times that of SEdb 1.0. (ii) Newly added transcription factor binding sites (TFBSs) in SEs identified by TF motifs and TF ChIP-seq data. (iii) Added comprehensive (epi)genetic annotations of SEs, including chromatin accessibility regions, methylation sites, chromatin interaction regions and topologically associating domains (TADs). (iv) Newly embedded and updated search and analysis tools, including 'Search SE by TF-based', 'Differential-Overlapping-SE analysis' and 'SE-based TF-Gene analysis'. (v) Newly provided quality control (QC) metrics for ChIP-seq processing. In summary, SEdb 2.0 is a comprehensive update of SEdb 1.0, which curates more SEs and annotation information than SEdb 1.0. SEdb 2.0 provides a friendly platform for researchers to more comprehensively clarify the important role of SEs in the biological process.

Identification of Autophagy-Related Targets of Berberine against Non-Small Cell Lung Cancer and Their Correlation with Immune Cell Infiltration By Combining Network Pharmacology, Molecular Docking, and Experimental Verification

OBJECTIVE

Non-small cell lung cancer (NSCLC) is the most common lung cancer type with high incidence. This study aimed to reveal the anti-NSCLC mechanisms of berberine and identify novel therapeutic targets.

METHODS

Berberine-related targets were acquired from SuperPred, SwissTargetPrediction, and GeneCards. NSCLC-re-lated targets were collected from GeneCards and DisGeNET. Differentially expressed genes (DEGs) were identified GEO database, UCSC Xena, and limma. GO and KEGG analyses were performed using clusterProfiler. Autophagy-related genes and transcriptional factors were collected from HADb and KnockTF, respectively. STRING and Cytoscape were used for PPI network analysis. Immune cell infiltration in NSCLC was assessed using CIBERSORT, and its correlation with autophagy-related targets was evaluated. Molecular docking was conducted using PyMOL and AutoDock. qRT-PCR and CCK-8 assay was used for in vitro verification.

RESULTS

Thirty intersecting targets of berberine-related targets, NSCLC-related targets, and DEGs were obtained. GO and KEGG analyses revealed that the intersecting targets were mainly implicated in oxidative stress, focal adhesion, and cell-substrate junction, as well as AGE-RAGE, relaxin, FoxO, and estrogen signaling pathways. Significantly, CAPN1, IKBKB, and SIRT2 were identified as the foremost autophagy-related targets, and 21 corresponding transcriptional factors were obtained. PPI network analysis showed that CAPN1, IKBKB, and SIRT2 interacted with 50 other genes. Fifty immune cell types, such as B cells naive, T cells CD8, T cells CD4 naive, T cells follicular helper, and monocytes, were implicated in NSCLC pathogenesis, and CAPN1, IKBKB, and SIRT2 were related to immune cells. Molecular docking revealed the favorable binding activity of berberine with CAPN1, IKBKB, and SIRT2. In vitro assays showed lower CAPN1, IKBKB, and SIRT2 expression in NSCLC cells than that in normal cells. Notably, berberine inhibited the viability and elevated CAPN1, IKBKB, and SIRT2 expression in NSCLC cells.

CONCLUSIONS

Berberine might treat NSCLC mainly by targeting CAPN1, IKBKB, and SIRT2.

Drug screening and biomarker gene investigation in cancer therapy through the human transcriptional regulatory network

A complex and vast biological network regulates all biological functions in the human body in a sophisticated manner, and abnormalities in this network can lead to disease and even cancer. The construction of a high-quality human molecular interaction network is possible with the development of experimental techniques that facilitate the interpretation of the mechanisms of drug treatment for cancer. We collected 11 molecular interaction databases based on experimental sources and constructed a human protein-protein interaction (PPI) network and a human transcriptional regulatory network (HTRN). A random walk-based graph embedding method was used to calculate the diffusion profiles of drugs and cancers, and a pipeline was constructed by using five similarity comparison metrics combined with a rank aggregation algorithm, which can be implemented for drug screening and biomarker gene prediction. Taking NSCLC as an example, curcumin was identified as a potentially promising anticancer drug from 5450 natural small molecules, and combined with differentially expressed genes, survival analysis, and topological ranking, we obtained BIRC5 (survivin), which is both a biomarker for NSCLC and a key target for curcumin. Finally, the binding mode of curcumin and survivin was explored using molecular docking. This work has a guiding significance for antitumor drug screening and the identification of tumor markers.

Broad misappropriation of developmental splicing profile by cancer in multiple organs

Oncogenesis mimics key aspects of embryonic development. However, the underlying mechanisms are incompletely understood. Here, we demonstrate that the splicing events specifically active during human organogenesis, are broadly reactivated in the organ-specific tumor. Such events are associated with key oncogenic processes and predict proliferation rates in cancer cell lines as well as patient survival. Such events preferentially target nitrosylation and transmembrane-region domains, whose coordinated splicing in multiple genes respectively affect intracellular transport and N-linked glycosylation. We infer critical splicing factors potentially regulating embryonic splicing events and show that such factors are potential oncogenic drivers and are upregulated specifically in malignant cells. Multiple complementary analyses point to MYC and FOXM1 as potential transcriptional regulators of critical splicing factors in brain and liver. Our study provides a comprehensive demonstration of a splicing-mediated link between development and cancer, and suggest anti-cancer targets including splicing events, and their upstream splicing and transcriptional regulators.

Data analysis guidelines for single-cell RNA-seq in biomedical studies and clinical applications

The application of single-cell RNA sequencing (scRNA-seq) in biomedical research has advanced our understanding of the pathogenesis of disease and provided valuable insights into new diagnostic and therapeutic strategies. With the expansion of capacity for high-throughput scRNA-seq, including clinical samples, the analysis of these huge volumes of data has become a daunting prospect for researchers entering this field. Here, we review the workflow for typical scRNA-seq data analysis, covering raw data processing and quality control, basic data analysis applicable for almost all scRNA-seq data sets, and advanced data analysis that should be tailored to specific scientific questions. While summarizing the current methods for each analysis step, we also provide an online repository of software and wrapped-up scripts to support the implementation. Recommendations and caveats are pointed out for some specific analysis tasks and approaches. We hope this resource will be helpful to researchers engaging with scRNA-seq, in particular for emerging clinical applications.

GPSAdb: a comprehensive web resource for interactive exploration of genetic perturbation RNA-seq datasets

Gene knock-out/down methods are commonly used to explore the functions of genes of interest, but a database that systematically collects perturbed data is not available currently. Manual curation of all the available human cell line perturbed RNA-seq datasets enabled us to develop a comprehensive human perturbation database (GPSAdb, https://www.gpsadb.com/). The current version of GPSAdb collected 3048 RNA-seq datasets associated with 1458 genes, which were knocked out/down by siRNA, shRNA, CRISPR/Cas9, or CRISPRi. The database provides full exploration of these datasets and generated 6096 new perturbed gene sets (up and down separately). GPSAdb integrated the gene sets and developed an online tool, genetic perturbation similarity analysis (GPSA), to identify candidate causal perturbations from differential gene expression data. In summary, GPSAdb is a powerful platform that aims to assist life science researchers to easily access and analyze public perturbed data and explore differential gene expression data in depth.

IReNA: Integrated regulatory network analysis of single-cell transcriptomes and chromatin accessibility profiles

Recently, single-cell RNA sequencing (scRNA-seq) and single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) have been developed to separately measure transcriptomes and chromatin accessibility profiles at the single-cell resolution. However, few methods can reliably integrate these data to perform regulatory network analysis. Here, we developed integrated regulatory network analysis (IReNA) for network inference through the integrated analysis of scRNA-seq and scATAC-seq data, network modularization, transcription factor enrichment, and construction of simplified intermodular regulatory networks. Using public datasets, we showed that integrated network analysis of scRNA-seq data with scATAC-seq data is more precise to identify known regulators than scRNA-seq data analysis alone. Moreover, IReNA outperformed currently available methods in identifying known regulators. IReNA facilitates the systems-level understanding of biological regulatory mechanisms and is available at https://github.com/jiang-junyao/IReNA.

PertOrg 1.0: a comprehensive resource of multilevel alterations induced in model organisms by in vivo genetic perturbation

Genetically modified organisms (GMOs) can be generated to model human genetic disease or plant disease resistance, and they have contributed to the exploration and understanding of gene function, physiology, disease onset and drug target discovery. Here, PertOrg (http://www.inbirg.com/pertorg/) was introduced to provide multilevel alterations in GMOs. Raw data of 58 707 transcriptome profiles and associated information, such as phenotypic alterations, were collected and curated from studies involving in vivo genetic perturbation (e.g. knockdown, knockout and overexpression) in eight model organisms, including mouse, rat and zebrafish. The transcriptome profiles from before and after perturbation were organized into 10 116 comparison datasets, including 122 single-cell RNA-seq datasets. The raw data were checked and analysed using widely accepted and standardized pipelines to identify differentially expressed genes (DEGs) in perturbed organisms. As a result, 8 644 148 DEGs were identified and deposited as signatures of gene perturbations. Downstream functional enrichment analysis, cell type analysis and phenotypic alterations were also provided when available. Multiple search methods and analytical tools were created and implemented. Furthermore, case studies were presented to demonstrate how users can utilize the database. PertOrg 1.0 will be a valuable resource aiding in the exploration of gene functions, biological processes and disease models.

TREAT: Therapeutic RNAs exploration inspired by artificial intelligence technology

Recent advances in RNA engineering have enabled the development of RNA-based therapeutics for a broad spectrum of applications. Developing RNA therapeutics start with targeted RNA screening and move to the drug design and optimization. However, existing target screening tools ignore noncoding RNAs and their disease-relevant regulatory relationships. And designing therapeutic RNAs encounters high computational complexity of multi-objective optimization to overcome the immunogenicity, instability and inefficient translational production. To unlock the therapeutic potential of noncoding RNAs and enable one-stop screening and design of therapeutic RNAs, we have built the platform TREAT. It incorporates 43,087,953 regulatory relationships between coding and noncoding genes from 81 biological networks under different physiological conditions. TREAT introduces graph representation learning with Random Walk Diffusions to perform disease-relevant target screening, in addition to the commonly utilized Topological Degree and PageRank algorithms. Design and optimization of large RNAs or interfering RNAs are both available. To reduce the computational complexity of multi-objective optimization for large RNA, we stratified the features into local and global features. The local features are evaluated on the fixed-length or dynamic-length local bins, whereas the latter are inspired by AI language models of protein sequence. Then the global assessment is performed on refined candidates, thus reducing the enormous search space. Overall, TREAT is a one-stop platform for the screening and designing of therapeutic RNAs, with particular attention to noncoding RNAs and cutting-edge AI technology embedded, leading the progress of innovative therapeutics for challenging diseases. TREAT is freely accessible at https://rna.org.cn/treat.

CohesinDB: a comprehensive database for decoding cohesin-related epigenomes, 3D genomes and transcriptomes in human cells

Cohesin is a multifunctional protein responsible for transcriptional regulation and chromatin organization. Cohesin binds to chromatin at tens of thousands of distinct sites in a conserved or tissue-specific manner, whereas the function of cohesin varies greatly depending on the epigenetic properties of specific chromatin loci. Cohesin also extensively mediates cis-regulatory modules (CRMs) and chromatin loops. Even though next-generation sequencing technologies have provided a wealth of information on different aspects of cohesin, the integration and exploration of the resultant massive cohesin datasets are not straightforward. Here, we present CohesinDB (https://cohesindb.iqb.u-tokyo.ac.jp), a comprehensive multiomics cohesin database in human cells. CohesinDB includes 2043 epigenomics, transcriptomics and 3D genomics datasets from 530 studies involving 176 cell types. By integrating these large-scale data, CohesinDB summarizes three types of 'cohesin objects': 751 590 cohesin binding sites, 957 868 cohesin-related chromatin loops and 2 229 500 cohesin-related CRMs. Each cohesin object is annotated with locus, cell type, classification, function, 3D genomics and cis-regulatory information. CohesinDB features a user-friendly interface for browsing, searching, analyzing, visualizing and downloading the desired information. CohesinDB contributes a valuable resource for all researchers studying cohesin, epigenomics, transcriptional regulation and chromatin organization.

Nonparametric single-cell multiomic characterization of trio relationships between transcription factors, target genes, and cis-regulatory regions

The epigenetic control of gene expression is highly cell-type and context specific. Yet, despite its complexity, gene regulatory logic can be broken down into modular components consisting of a transcription factor (TF) activating or repressing the target gene expression through its binding to a cis-regulatory region. We propose a nonparametric approach, TRIPOD, to detect and characterize the three-way relationships between a TF, its target gene, and the accessibility of the TF's binding site using single-cell RNA and ATAC multiomic data. We apply TRIPOD to interrogate the cell-type-specific regulatory logic in peripheral blood mononuclear cells and contrast our results to detections from enhancer databases, cis-eQTL studies, ChIP-seq experiments, and TF knockdown/knockout studies. We then apply TRIPOD to mouse embryonic brain data and identify regulatory relationships, validated by ChIP-seq and PLAC-seq. Finally, we demonstrate TRIPOD on the SHARE-seq data of differentiating mouse hair follicle cells and identify lineage-specific regulation supported by histone marks and super-enhancer annotations. A record of this paper's transparent peer review process is included in the supplemental information.

GREAP: a comprehensive enrichment analysis software for human genomic regions

The rapid development of genomic high-throughput sequencing has identified a large number of DNA regulatory elements with abundant epigenetics markers, which promotes the rapid accumulation of functional genomic region data. The comprehensively understanding and research of human functional genomic regions is still a relatively urgent work at present. However, the existing analysis tools lack extensive annotation and enrichment analytical abilities for these regions. Here, we designed a novel software, Genomic Region sets Enrichment Analysis Platform (GREAP), which provides comprehensive region annotation and enrichment analysis capabilities. Currently, GREAP supports 85 370 genomic region reference sets, which cover 634 681 107 regions across 11 different data types, including super enhancers, transcription factors, accessible chromatins, etc. GREAP provides widespread annotation and enrichment analysis of genomic regions. To reflect the significance of enrichment analysis, we used the hypergeometric test and also provided a Locus Overlap Analysis. In summary, GREAP is a powerful platform that provides many types of genomic region sets for users and supports genomic region annotations and enrichment analyses. In addition, we developed a customizable genome browser containing >400 000 000 customizable tracks for visualization. The platform is freely available at http://www.liclab.net/Greap/view/index.

Identification and Validation of Autophagy-Related Genes in Primary Ovarian Insufficiency by Gene Expression Profile and Bioinformatic Analysis

Background

To investigate the relationship between primary ovarian insufficiency and autophagy, we detected and got the expression profile of human granulosa cell line SVOG, which was with or without LPS induced. The expression profile was analyzed with the focus on the autophagy genes, among which hub genes were identified.

Results

Totally, 6 genes were selected as candidate hub genes which might correlate with the process of primary ovarian insufficiency. The expression of hub genes was then validated by quantitative real-time PCR and two of them had significant expression change. Bioinformatics analysis was performed to observe the features of hub genes, including hub gene-RBP/TF/miRNA/drug network construction, functional analysis, and protein-protein interaction network. Pearson's correlation analysis was also performed to identify the correlation between hub genes and autophagy genes, among which there were four autophagy genes significantly correlated with hub genes, including ATG4B, ATG3, ATG13, and ULK1.

Conclusion

The results indicated that autophagy might play an essential role in the process and underlying molecular mechanism of primary ovarian insufficiency, which was revealed for the first time and may help to provide a molecular foundation for the development of diagnostic and therapeutic approaches for primary ovarian insufficiency.