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.

PBAC: A pathway-based attention convolution neural network for predicting clinical drug treatment responses.

Precise and personalized drug application is crucial in the clinical treatment of complex diseases. Although neural networks offer a new approach to improving drug strategies, their internal structure is difficult to interpret. Here, we propose PBAC (Pathway-Based Attention Convolution neural network), which integrates a deep learning framework and attention mechanism to address the complex biological pathway information, thereby provide a biology function-based robust drug responsiveness prediction model. PBAC has four layers: gene-pathway layer, attention layer, convolution layer and fully connected layer. PBAC improves the performance of predicting drug responsiveness by focusing on important pathways, helping us understand the mechanism of drug action in diseases. We validated the PBAC model using data from four chemotherapy drugs (Bortezomib, Cisplatin, Docetaxel and Paclitaxel) and 11 immunotherapy datasets. In the majority of datasets, PBAC exhibits superior performance compared to traditional machine learning methods and other research approaches (area under curve = 0.81, the area under the precision-recall curve = 0.73). Using PBAC attention layer output, we identified some pathways as potential core cancer regulators, providing good interpretability for drug treatment prediction. In summary, we presented PBAC, a powerful tool to predict drug responsiveness based on the biology pathway information and explore the potential cancer-driving pathways.

Heterogeneity analysis of the immune microenvironment in laryngeal carcinoma revealed potential prognostic biomarkers.

Laryngeal squamous cell cancer (LSCC) is the second most prevalent malignancy occurring in the head and neck with a high incidence and mortality rate. Immunotherapy has recently become an emerging treatment for cancer. It is therefore essential to explore the role of tumour immunity in laryngeal cancer. Our study first delineated and evaluated the comprehensive immune infiltration landscapes of the tumour microenvironment in LSCC. A hierarchical clustering method was applied to classify the LSCC samples into two groups (high- and low-infiltration groups). We found that individuals with low immune infiltration characteristics had significantly better survival than those in the high-infiltration group, possibly because of the elevated infiltration of immune suppressive cells, such as regulatory T cells and myeloid-derived suppressor cells, in the high-infiltration group. Differentially expressed genes between two groups were involved in some immune-related terms, such as antigen processing and presentation. A univariate Cox analysis and least absolute shrinkage and selection operator analysis were performed to identify an immune gene-set-based prognostic signature (IBPS) to assess the risk of LSCC. The prognostic model comprising six IBPSs was successfully verified to be robust in different cohorts. The expression of the six IBPSs was detected by immunohistochemistry in 110 cases of LSCC. In addition, different inflammatory profiles and immune checkpoint landscape of LSCC were found between two groups. Hence, our model could serve as a candidate immunotherapeutic biomarker and potential therapeutic target for laryngeal cancer.

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.

TF-Marker: a comprehensive manually curated database for transcription factors and related markers in specific cell and tissue types in human.

Transcription factors (TFs) play key roles in biological processes and are usually used as cell markers. The emerging importance of TFs and related markers in identifying specific cell types in human diseases increases the need for a comprehensive collection of human TFs and related markers sets. Here, we developed the TF-Marker database (TF-Marker, http://bio.liclab.net/TF-Marker/), aiming to provide cell/tissue-specific TFs and related markers for human. By manually curating thousands of published literature, 5905 entries including information about TFs and related markers were classified into five types according to their functions: (i) TF: TFs which regulate expression of the markers; (ii) T Marker: markers which are regulated by the TF; (iii) I Marker: markers which influence the activity of TFs; (iv) TFMarker: TFs which play roles as markers and (v) TF Pmarker: TFs which play roles as potential markers. The 5905 entries of TF-Marker include 1316 TFs, 1092 T Markers, 473 I Markers, 1600 TFMarkers and 1424 TF Pmarkers, involving 383 cell types and 95 tissue types in human. TF-Marker further provides a user-friendly interface to browse, query and visualize the detailed information about TFs and related markers. We believe TF-Marker will become a valuable resource to understand the regulation patterns of different tissues and cells.

TcoFBase: a comprehensive database for decoding the regulatory transcription co-factors in human and mouse.

Transcription co-factors (TcoFs) play crucial roles in gene expression regulation by communicating regulatory cues from enhancers to promoters. With the rapid accumulation of TcoF associated chromatin immunoprecipitation sequencing (ChIP-seq) data, the comprehensive collection and integrative analyses of these data are urgently required. Here, we developed the TcoFBase database (http://tcof.liclab.net/TcoFbase), which aimed to document a large number of available resources for mammalian TcoFs and provided annotations and enrichment analyses of TcoFs. TcoFBase curated 2322 TcoFs and 6759 TcoFs associated ChIP-seq data from over 500 tissues/cell types in human and mouse. Importantly, TcoFBase provided detailed and abundant (epi) genetic annotations of ChIP-seq based TcoF binding regions. Furthermore, TcoFBase supported regulatory annotation information and various functional annotations for TcoFs. Meanwhile, TcoFBase embedded five types of TcoF regulatory analyses for users, including TcoF gene set enrichment, TcoF binding genomic region annotation, TcoF regulatory network analysis, TcoF-TF co-occupancy analysis and TcoF regulatory axis analysis. TcoFBase was designed to be a useful resource that will help reveal the potential biological effects of TcoFs and elucidate TcoF-related regulatory mechanisms.

RNAactDrug: a comprehensive database of RNAs associated with drug sensitivity from multi-omics data.

Drug sensitivity has always been at the core of individualized cancer chemotherapy. However, we have been overwhelmed by large-scale pharmacogenomic data in the era of next-generation sequencing technology, which makes it increasingly challenging for researchers, especially those without bioinformatic experience, to perform data integration, exploration and analysis. To bridge this gap, we developed RNAactDrug, a comprehensive database of RNAs associated with drug sensitivity from multi-omics data, which allows users to explore drug sensitivity and RNA molecule associations directly. It provides association data between drug sensitivity and RNA molecules including mRNAs, long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) at four molecular levels (expression, copy number variation, mutation and methylation) from integrated analysis of three large-scale pharmacogenomic databases (GDSC, CellMiner and CCLE). RNAactDrug currently stores more than 4 924 200 associations of RNA molecules and drug sensitivity at four molecular levels covering more than 19 770 mRNAs, 11 119 lncRNAs, 438 miRNAs and 4155 drugs. A user-friendly interface enriched with various browsing sections augmented with advance search facility for querying the database is offered for users retrieving. RNAactDrug provides a comprehensive resource for RNA molecules acting in drug sensitivity, and it could be used to prioritize drug sensitivity-related RNA molecules, further promoting the identification of clinically actionable biomarkers in drug sensitivity and drug development more cost-efficiently by making this knowledge accessible to both basic researchers and clinical practitioners. Database URL: http://bio-bigdata.hrbmu.edu.cn/RNAactDrug.

Identification and comprehensive characterization of lncRNAs with copy number variations and their driving transcriptional perturbed subpathways reveal functional significance for cancer.

Numerous studies have shown that copy number variation (CNV) in lncRNA regions play critical roles in the initiation and progression of cancer. However, our knowledge about their functionalities is still limited. Here, we firstly provided a computational method to identify lncRNAs with copy number variation (lncRNAs-CNV) and their driving transcriptional perturbed subpathways by integrating multidimensional omics data of cancer. The high reliability and accuracy of our method have been demonstrated. Then, the method was applied to 14 cancer types, and a comprehensive characterization and analysis was performed. LncRNAs-CNV had high specificity in cancers, and those with high CNV level may perturb broad biological functions. Some core subpathways and cancer hallmarks widely perturbed by lncRNAs-CNV were revealed. Moreover, subpathways highlighted the functional diversity of lncRNAs-CNV in various cancers. Survival analysis indicated that functional lncRNAs-CNV could be candidate prognostic biomarkers for clinical applications, such as ST7-AS1, CDKN2B-AS1 and EGFR-AS1. In addition, cascade responses and a functional crosstalk model among lncRNAs-CNV, impacted genes, driving subpathways and cancer hallmarks were proposed for understanding the driving mechanism of lncRNAs-CNV. Finally, we developed a user-friendly web interface-LncCASE (http://bio-bigdata.hrbmu.edu.cn/LncCASE/) for exploring lncRNAs-CNV and their driving subpathways in various cancer types. Our study identified and systematically characterized lncRNAs-CNV and their driving subpathways and presented valuable resources for investigating the functionalities of non-coding variations and the mechanisms of tumorigenesis.

Characterizing heterogeneity of non-small cell lung tumour microenvironment to identify signature prognostic genes.

Growing evidence has highlighted the immune response as an important feature of carcinogenesis and therapeutic efficacy in non-small cell lung cancer (NSCLC). This study focused on the characterization of immune infiltration profiling in patients with NSCLC and its correlation with survival outcome. All TCGA samples were divided into three heterogeneous clusters based on immune cell profiles: cluster 1 ('low infiltration' cluster), cluster 2 ('heterogeneous infiltration' cluster) and cluster 3 ('high infiltration' cluster). The immune cells were responsible for a significantly favourable prognosis for the 'high infiltration' community. Cluster 1 had the lowest cytotoxic activity, tumour-infiltrating lymphocytes and interferon-gamma (IFN-γ), as well as immune checkpoint molecules expressions. In addition, MHC-I and immune co-stimulator were also found to have lower cluster 1 expressions, indicating a possible immune escape mechanism. A total of 43 differentially expressed genes (DEGs) that overlapped among the groups were determined based on three clusters. Finally, based on a univariate Cox regression model, prognostic immune-related genes were identified and combined to construct a risk score model able to predict overall survival (OS) rates in the validation datasets.

Comprehensive analysis of the whole coding and non-coding RNA transcriptome expression profiles and construction of the circRNA-lncRNA co-regulated ceRNA network in laryngeal squamous cell carcinoma.

Recently, accumulating evidence has demonstrated that non-coding RNAs (ncRNAs) play a vital role in oncogenicity. Nevertheless, the regulatory mechanisms and functions remain poorly understood, especially for lncRNAs and circRNAs. In this study, we simultaneously detected, for the first time, the expression profiles of the whole transcriptome, including miRNA, circRNA and lncRNA + mRNA, in five pairs of laryngeal squamous cell carcinoma (LSCC) and matched non-carcinoma tissues by microarrays. Five miRNAs, four circRNAs, three lncRNAs and five mRNAs that were dysregulated were selected to confirm the verification of the microarray data by quantitative real-time PCR (qRT-PCR) in 20 pairs of LSCC samples. We constructed LSCC-related competing endogenous RNA (ceRNA) networks of lncRNAs and circRNAs (circRNA or lncRNA-miRNA-mRNA) respectively. Functional annotation revealed the lncRNA-mediated ceRNA network were enriched for genes involved in the tumor-associated pathways. Hsa_circ_0033988 with the highest degree in the circRNA-mediated ceRNA network was associated with fatty acid degradation, which was responsible for the depletion of fat in tumor-associated cachexia. Finally, to clarify the ncRNA co-regulation mechanism, we constructed a circRNA-lncRNA co-regulated network by integrating the above two networks and identified 9 modules for further study. A subnetwork of module 2 with the most dysregulated microRNAs was extracted to establish the ncRNA-involved TGF-ß-associated pathway. In conclusion, our findings provide a high-throughput microarray data of the coding and non-coding RNAs and establish the foundation for further functional research on the ceRNA regulatory mechanism of non-coding RNAs in LSCC.

Identifying subpathway signatures for individualized anticancer drug response by integrating multi-omics data.

BACKGROUND: Individualized drug response prediction is vital for achieving personalized treatment of cancer and moving precision medicine forward. Large-scale multi-omics profiles provide unprecedented opportunities for precision cancer therapy. METHODS: In this study, we propose a pipeline to identify subpathway signatures for anticancer drug response of individuals by integrating the comprehensive contributions of multiple genetic and epigenetic (gene expression, copy number variation and DNA methylation) alterations. RESULTS: Totally, 46 subpathway signatures associated with individual responses to different anticancer drugs were identified based on five cancer-drug response datasets. We have validated the reliability of subpathway signatures in two independent datasets. Furthermore, we also demonstrated these multi-omics subpathway signatures could significantly improve the performance of anticancer drug response prediction. In-depth analysis of these 46 subpathway signatures uncovered the essential roles of three omics types and the functional associations underlying different anticancer drug responses. Patient stratification based on subpathway signatures involved in anticancer drug response identified subtypes with different clinical outcomes, implying their potential roles as prognostic biomarkers. In addition, a landscape of subpathways associated with cellular responses to 191 anticancer drugs from CellMiner was provided and the mechanism similarity of drug action was accurately unclosed based on these subpathways. Finally, we constructed a user-friendly web interface-CancerDAP ( http://bio-bigdata.hrbmu.edu.cn/CancerDAP/ ) available to explore 2751 subpathways relevant with 191 anticancer drugs response. CONCLUSIONS: Taken together, our study identified and systematically characterized subpathway signatures for individualized anticancer drug response prediction, which may promote the precise treatment of cancer and the study for molecular mechanisms of drug actions.

Characterizing and optimizing human anticancer drug targets based on topological properties in the context of biological pathways.

One of the challenging problems in drug discovery is to identify the novel targets for drugs. Most of the traditional methods for drug targets optimization focused on identifying the particular families of "druggable targets", but ignored their topological properties based on the biological pathways. In this study, we characterized the topological properties of human anticancer drug targets (ADTs) in the context of biological pathways. We found that the ADTs tended to present the following seven topological properties: influence the number of the pathways related to cancer, be localized at the start or end of the pathways, interact with cancer related genes, exhibit higher connectivity, vulnerability, betweenness, and closeness than other genes. We first ranked ADTs based on their topological property values respectively, then fused them into one global-rank using the joint cumulative distribution of an N-dimensional order statistic to optimize human ADTs. We applied the optimization method to 13 anticancer drugs, respectively. Results demonstrated that over 70% of known ADTs were ranked in the top 20%. Furthermore, the performance for mercaptopurine was significant: 6 known targets (ADSL, GMPR2, GMPR, HPRT1, AMPD3, AMPD2) were ranked in the top 15 and other four out of the top 15 (MAT2A, CDKN1A, AREG, JUN) have the potentialities to become new targets for cancer therapy.

Subpathway-GM: identification of metabolic subpathways via joint power of interesting genes and metabolites and their topologies within pathways.

Various 'omics' technologies, including microarrays and gas chromatography mass spectrometry, can be used to identify hundreds of interesting genes, proteins and metabolites, such as differential genes, proteins and metabolites associated with diseases. Identifying metabolic pathways has become an invaluable aid to understanding the genes and metabolites associated with studying conditions. However, the classical methods used to identify pathways fail to accurately consider joint power of interesting gene/metabolite and the key regions impacted by them within metabolic pathways. In this study, we propose a powerful analytical method referred to as Subpathway-GM for the identification of metabolic subpathways. This provides a more accurate level of pathway analysis by integrating information from genes and metabolites, and their positions and cascade regions within the given pathway. We analyzed two colorectal cancer and one metastatic prostate cancer data sets and demonstrated that Subpathway-GM was able to identify disease-relevant subpathways whose corresponding entire pathways might be ignored using classical entire pathway identification methods. Further analysis indicated that the power of a joint genes/metabolites and subpathway strategy based on their topologies may play a key role in reliably recalling disease-relevant subpathways and finding novel subpathways.

Topologically inferring risk-active pathways toward precise cancer classification by directed random walk.

MOTIVATION: The accurate prediction of disease status is a central challenge in clinical cancer research. Microarray-based gene biomarkers have been identified to predict outcome and outperform traditional clinical parameters. However, the robustness of the individual gene biomarkers is questioned because of their little reproducibility between different cohorts of patients. Substantial progress in treatment requires advances in methods to identify robust biomarkers. Several methods incorporating pathway information have been proposed to identify robust pathway markers and build classifiers at the level of functional categories rather than of individual genes. However, current methods consider the pathways as simple gene sets but ignore the pathway topological information, which is essential to infer a more robust pathway activity. RESULTS: Here, we propose a directed random walk (DRW)-based method to infer the pathway activity. DRW evaluates the topological importance of each gene by capturing the structure information embedded in the directed pathway network. The strategy of weighting genes by their topological importance greatly improved the reproducibility of pathway activities. Experiments on 18 cancer datasets showed that the proposed method yielded a more accurate and robust overall performance compared with several existing gene-based and pathway-based classification methods. The resulting risk-active pathways are more reliable in guiding therapeutic selection and the development of pathway-specific therapeutic strategies. AVAILABILITY: DRW is freely available at http://210.46.85.180:8080/DRWPClass/

The detection of risk pathways, regulated by miRNAs, via the integration of sample-matched miRNA-mRNA profiles and pathway structure.

The use of genome-wide, sample-matched miRNA (miRNAs)-mRNA expression data provides a powerful tool for the investigation of miRNAs and genes involved in diseases. The identification of miRNA-regulated pathways has been crucial for analysis of the role of miRNAs. However, the classical identification method fails to consider the structural information of pathways and the regulation of miRNAs simultaneously. We proposed a method that simultaneously integrated the change in gene expression and structural information in order to identify pathways. Our method used fold changes in miRNAs and gene products, along with the quantification of the regulatory effect on target genes, to measure the change in gene expression. Topological characteristics were investigated to measure the influence of gene products on entire pathways. Through the analysis of multiple myeloma and prostate cancer expression data, our method was proven to be effective and reliable in identifying disease risk pathways that are regulated by miRNAs. Further analysis showed that the structure of a pathway plays a crucial role in the recognition of the pathway as a factor in disease risk.

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.

Characterization of the metabolic alteration-modulated tumor microenvironment mediated by TP53 mutation and hypoxia.

BACKGROUND: TP53 mutation and hypoxia play an essential role in cancer progression. However, the metabolic reprogramming and tumor microenvironment (TME) heterogeneity mediated by them are still not fully understood. METHODS: The multi-omics data of 32 cancer types and immunotherapy cohorts were acquired to comprehensively characterize the metabolic reprogramming pattern and the TME across cancer types and explore immunotherapy candidates. An assessment model for metabolic reprogramming was established by integration of multiple machine learning methods, including lasso regression, neural network, elastic network, and survival support vector machine (SVM). Pharmacogenomics analysis and in vitro assay were conducted to identify potential therapeutic drugs. RESULTS: First, we identified metabolic subtype A (hypoxia-TP53 mutation subtype) and metabolic subtype B (non-hypoxia-TP53 wildtype subtype) in hepatocellular carcinoma (HCC) and showed that metabolic subtype A had an "immune inflamed" microenvironment. Next, we established an assessment model for metabolic reprogramming, which was more effective compared to the traditional prognostic indicators. Then, we identified a potential targeting drug, teniposide. Finally, we performed the pan-cancer analysis to illustrate the role of metabolic reprogramming in cancer and found that the metabolic alteration (MA) score was positively correlated with tumor mutational burden (TMB), neoantigen load, and homologous recombination deficiency (HRD) across cancer types. Meanwhile, we demonstrated that metabolic reprogramming mediated a potential immunotherapy-sensitive microenvironment in bladder cancer and validated it in an immunotherapy cohort. CONCLUSION: Metabolic alteration mediated by hypoxia and TP53 mutation is associated with TME modulation and tumor progression across cancer types. In this study, we analyzed the role of metabolic alteration in cancer and propose a predictive model for cancer prognosis and immunotherapy responsiveness. We also explored a potential therapeutic drug, teniposide.

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.

A novel cuproptosis-related diagnostic gene signature and differential expression validation in atherosclerosis.

Atherosclerosis (AS) is a major contributor to morbidity and mortality worldwide. However, the molecular mechanisms and mediator molecules involved remain largely unknown. Copper, which plays an essential role in cardiovascular disease, has been suggested as a potential risk factor. Copper homeostasis is closely related to the occurrence and development of AS. Recently, a new cell death pathway called cuproptosis has been discovered, which is driven by intracellular copper excess. However, no previous studies have reported a relationship between cuproptosis and AS. In this study, we integrated bulk and single-cell sequencing data to screen and identify key cuproptosis-related genes in AS. We used correlation analysis, enrichment analysis, random forest, and other bioinformatics methods to reveal their relationships. Our findings report, for the first time, the involvement of cuproptosis-related genes FDX1, SLC31A1, and GLS in atherogenesis. FDX1 and SLC31A1 were upregulated, while GLS was downregulated in atherosclerotic plaque. Receiver operating characteristic curves demonstrate their potential diagnostic value for AS. Additionally, we confirm that GLS is mainly expressed in vascular smooth muscle cells, and SLC31A1 is mainly localized in macrophages of atherosclerotic lesions in experiments. These findings shed light on the cuproptosis landscape and potential diagnostic biomarkers for AS, providing further evidence about the vital role of cuproptosis in atherosclerosis progression.

Landscape and significance of human super enhancer-driven core transcription regulatory circuitry.

A core transcription regulatory circuitry (CRC) is an interconnected self-regulatory circuitry that is formed by a group of core transcription factors (TFs). These core TFs collectively regulate gene expression by binding not only to their own super enhancers (SEs) but also to the SEs of one another. For most human tissue/cell types, a global view of CRCs and core TFs has not been generated. Here, we identified numerous CRCs using two identification methods and detailed the landscape of the CRCs driven by SEs in large cell/tissue samples. The comprehensive biological analyses, including sequence conservation, CRC activity and genome binding affinity were conducted for common TFs, moderate TFs, and specific TFs, which exhibit different biological features. The local module located from the common CRC network highlighted the essential functions and prognostic performance. The tissue-specific CRC network was highly related to cell identity. Core TFs in tissue-specific CRC networks exhibited disease markers, and had regulatory potential for cancer immunotherapy. Moreover, a user-friendly resource named CRCdb (http://www.licpathway.net/crcdb/index.html) was developed, which contained the detailed information of CRCs and core TFs used in this study, as well as other interesting results, such as the most representative CRC, frequency of TFs, and indegree/outdegree of TFs.