Genome-wide association studies (GWAS) have identified numerous genetic variants associated with diseases and traits. However, the functional interpretation of these variants remains challenging. Expression quantitative trait loci (eQTLs) have been widely used to identify mutations linked to disease, yet they explain only 20–50% of disease-related variants. Single-cell eQTLs (sc-eQTLs) studies provide an immense opportunity to identify new disease risk genes with expanded eQTL scales and transcriptional regulation at a much finer resolution. However, there is no comprehensive database dedicated to single-cell eQTLs that users can use to search, analyse and visualize them. Therefore, we developed the scQTLbase (http://bioinfo.szbl.ac.cn/scQTLbase), the first integrated human sc-eQTLs portal, featuring 304 datasets spanning 57 cell types and 95 cell states. It contains ∼16 million SNPs significantly associated with cell-type/state gene expression and ∼0.69 million disease-associated sc-eQTLs from 3 333 traits/diseases. In addition, scQTLbase offers sc-eQTL search, gene expression visualization in UMAP plots, a genome browser, and colocalization visualization based on the GWAS dataset of interest. scQTLbase provides a one-stop portal for sc-eQTLs that will significantly advance the discovery of disease susceptibility genes.

Background. Ovarian cancer (OC) is a malignancy exhibiting high mortality in female tumors. Glycosylation is a posttranslational modification of proteins but research has failed to demonstrate a systematic link between glycosylation‐related signatures and tumor environment of OC. Purpose. This study is aimed at developing a novel model with glycosylation‐related messenger RNAs (GRmRNAs) to predict the prognosis and immune function in OC patients. Methods. The transcriptional profiles and clinical phenotypes of OC patients were collected from the Gene Expression Omnibus and The Cancer Genome Atlas databases. A weighted gene coexpression network analysis and machine learning were performed to find the optimal survival‐related GRmRNAs. Least absolute shrinkage and selection operator regression (LASSO) and Cox regression were carried out to calculate the coefficients of each GRmRNA and compute the risk score of each patient as well as develop a prognostic model. A nomogram model was constructed, and several algorithms were used to investigate the relationship between risk subtypes and immune‐infiltrating levels. Results. A total of four signatures (ALG8, DCTN4, DCTN6, and UBB) were determined to calculate the risk scores, classifying patients into the high‐and low‐risk groups. High‐risk patients exhibited significantly poorer survival outcomes, and the established nomogram model had a promising prediction for OC patients’ prognosis. Tumor purity and tumor mutation burden were negatively correlated with risk scores. In addition, risk scores held statistical associations with pathway signatures such as Wnt, Hippo, and reactive oxygen species, and nonsynonymous mutation counts. Conclusion. The currently established risk scores based on GRmRNAs can accurately predict the prognosis, the immune microenvironment, and the immunotherapeutic efficacy of OC patients.

Purpose. This study is aimed at systematically analyzing the expression, function, and prognostic value of six transmembrane epithelial antigen of the prostate 1 (STEAP1) in various cancers. Methods. The expressions of STEAP1 between normal and tumor tissues were analyzed using TCGA and GTEx. Clinicopathologic data was collected from GEPIA and TCGA. Prognostic analysis was conducted by Cox proportional hazard regression and Kaplan‐Meier survival. DNA methylation, mutation features, and molecular subtypes of cancers were also investigated. The top‐100 coexpressed genes with STEAP1 were involved in functional enrichment analysis. ESTIMATE algorithm was used to analyze the correlation between STEAP1 and immunity value. The relationships of STEAP1 and biomarkers including tumor mutational burden (TMB), microsatellite instability (MSI), and stemness score as well as chemosensitivity were also illustrated. Results. Among 33 cancers, STEAP1 was overexpressed in 19 cancers such as cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), colon adenocarcinoma, and lymphoid neoplasm diffuse large B cell lymphoma while was downregulated in 5 cancers such as adrenocortical carcinoma, breast invasive carcinoma (BRCA), and kidney chromophobe renal cell carcinoma. STEAP1 has significant prognostic relationships in multiple cancers. 15 cancers exhibited differences of DNA methylation including bladder urothelial carcinoma, BRCA, and CESC. STEAP1 expression was positively correlated to immune molecules especially in thyroid carcinoma and negatively especially in uveal melanoma. STEAP1 was associated with TMB and MSI in certain cancers. In addition, STEAP1 was connected with increased chemosensitivity of drugs such as trametinib and pimasertib. Conclusions. STEAP1 was an underlying target for prognostic prediction in different cancer types and a potential biomarker of TMB, MSI, tumor microenvironment, and chemosensitivity.

Glioma is a malignancy with the highest mortality in central nervous system disorders. Here, we implemented the computational tools based on CRISPR/Cas9 to predict the clinical outcomes and biological characteristics of low-grade glioma (LGG). The transcriptional expression profiles and clinical phenotypes of LGG patients were retrieved from The Cancer Genome Atlas and Chinese Glioma Genome Atlas. The CERES algorithm was used to screen for LGG-lethal genes. Cox regression and random survival forest were adopted for survival-related gene selection. Nonnegative matrix factorization distinguished patients into different clusters. Single-sample gene set enrichment analysis was employed to create a novel CRISPR/Cas9 screening potential index (CCSPI), and patients were stratified into low- and high-CCSPI groups. Survival analysis, area under the curve values (AUCs), nomogram, and tumor microenvironment exploration were included for the model validation. A total of 20 essential genes in LGG were used to classify patients into two clusters and construct the CCSPI system. High-CCSPI patients were associated with a worse prognosis of both training and validation set (p < 0.0001) and higher immune fractions than low-CCSPI individuals. The CCSPI system had a promising performance with 1-, 3-, and 5-year AUCs of 0.816, 0.779, 0.724, respectively, and the C-index of the nomogram model reached 0.743 (95% CI = 0.725–0.760). Immune-infiltrating cells and immune checkpoints such as PD-1/PD-L1 and POLD3 were positively associated with CCSPI. In conclusion, the CCSPI had prognostic value in LGG, and the model will deepen our cognition of the interaction between the CNS and immune system in different LGG subtypes.

Accurately deciphering tumor cellular heterogeneity is crucial for developing effective treatments. While single-cell epigenomics assays offer powerful insights into studying tumor heterogeneity beyond the transcriptional level, their application remains a significant challenge. To address this, we introduce Multimodal-based Analysis of scATAC-Seq data (MAAS), a method designed to identify functional tumor cell subpopulations and reconstruct their lineages from single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq). MAAS uniquely integrates multimodal information, encompassing chromatin accessibility, copy number variations (CNVs), and single-nucleotide variants, through self-expressive multimodal matrix factorization. MAAS outperforms existing methods in accuracy and robustness on both simulated and real datasets. When applied to multiple glioma scATAC-seq datasets, MAAS uncovered a previously hidden tumor cell subpopulation resistant to temozolomide treatment. Furthermore, MAAS effectively detected clinically relevant subpopulations with low CNVs, such as those found in pediatric ependymoma and B-cell lymphoma. Additionally, we developed a novel MAAS-derived clinical signature that provides superior prognostic prediction than traditional clinicopathologic features across multiple cancer types. In summary, MAAS significantly advances the identification of clinically relevant tumor cell subpopulations, thereby accelerating the discovery of potential therapeutic targets.