Wang, Junwen

Abstract Expression quantitative trait locus (eQTL) mapping is used to identify the functional link between a genomic variant and a gene’s expression. A significant eQTL association does not mean a causal relationship or mechanism, and further investigation is needed to understand how a single-nucleotide polymorphism (SNP) impacts gene expression. One of the most plausible explanations for eQTL is that a genomic variant affects transcription factor (TF) binding and thus impacts its regulation on target genes (TGs). However, the current eQTL model does not provide information on the TF and how its regulation is mediated by the SNP’s genotypes. Here, we propose a new method called differential co-expression QTL (ceQTL) among different alleles using Chow statistics to specifically detect eQTLs that are bound by a particular TF. We start with building a trio of TF, its TG, and related SNP, and then test the significant coefficient difference among different genotypes of the SNP. We applied this ceQTL model to simulated data and the lung tissue datasets from the genotype-tissue expression project. The simulated data results showed that the model was robust to detect true ceQTLs at variable sample sizes and different minor allele frequencies as measured by area under the curve. Our tool also performed a TF binding affinity analysis to add another layer of evidence for functional interpretation. In summary, ceQTL analysis provides a more interpretable and biological insight into the mechanism of eQTL and transcriptomic regulation, which would help us better understand how genomic variants affect phenotypes and diseases.

Abstract Glioblastoma (GBM) invasion is shaped by dynamic clonal evolution within the spatial and cellular architecture of the brain. To investigate this process, we integrated multiparametric MRI, whole-exome sequencing (WES), bulk RNA-seq and single-nucleus RNA-seq (snRNA-seq) from 69 IDH-wildtype GBM patients, using spatially resolved biopsies from contrast-enhancing (CE) core tumor and non-enhancing (NE) invasive regions. Phylogenetic analysis revealed distinct clonal architectures. CE-enriched clones carried truncal mutations in conanical drivers of cell cycle and RTK signaling (TP53, PTEN, EGFR, PDGFRA/B), and private alterations in DNA repair genes (CDH11, PAX7, RAD51B). In contrast, NE-enriched clones exhibited private mutations in MAPK and Hippo pathways components (STAG2, BRAF, MAP2K3, CTNNA1), suggesting region-specific evolutionary trajectories. Furthermore, snRNAseq cell-cell communication analysis uncovered enriched ligand-receptor signaling in NE regions between malignant cells and brain-resident populations, including PTN—PTPRZ1 (myeloid), NRG-ERBB4 and SLIT-ROBO (neuronal, astrocytic, and oligodendrocytic). Bulk RNAseq deconvolution confirmed increased representation of neurons and glial cells in biopsies harboring NE-specific clones. These findings define a unique population of GBM clones localized at the infiltrative margin, marked by private MAPK/Hippo mutations and enhanced intracellular communication with the surrounding brain microenvironment. The recurrent engagement of neural-glial signaling pathways suggests that these NE-associated clones co-opt normal cellular networks to support invasion, immune evasion, and persistence following surgical resection. This study establishes a mechanistic link between spatial clonal evolution and tumor-host communication at the GBM infiltrative region and provide a framework for targeting cell-cell signaling circuits that sustain invasive tumor ecosystems and escape therapeutic treatment.

Abstract Background Lung adenocarcinoma in never-smokers is a major public health burden, especially among East Asian women. Polygenic risk scores (PRSs) are promising for risk stratification but are primarily developed in European-ancestry populations. We aimed to develop and validate single- and multi-ancestry PRSs for East Asian never-smokers to improve lung adenocarcinoma risk prediction. Methods PRSs were developed using genome-wide association study summary statistics from East Asian (8002 cases; 20 782 controls) and European (2058 cases; 5575 controls) populations. Single-ancestry models included a PRS with 25 single nucleotide variants (PRS-25), PRS clumping and thresholding (CT), and LDpred2; multi-ancestry models included LDpred2 plus European PRS-128, PRS-CSx, and clumping and thresholding, super learning, empirical Bayes (CT-SLEB). Performance was evaluated in independent East Asian data from the Female Lung Cancer Consortium and externally validated in the Nanjing Lung Cancer Cohort. We assessed predictive accuracy via area under the curve operating procedure (AUC), with 10-year and (ages 30-80 years) absolute risks estimates. Results The best multi-ancestry PRS, using East Asian and European data via CT-SLEB, outperformed the best East Asian-only (LDpred2; AUC = 0.629, 95% confidence interval [CI] = 0.618 to 0.641), achieving an AUC of 0.640 (95% CI = 0.629 to 0.653) and odds ratio of 1.71 (95% CI = 1.61 to 1.82) per standard deviation increase. Nanjing Lung Cancer Cohort validation confirmed robust performance (AUC = 0.649, 95% CI = 0.623 to 0.676). The top 20% PRS group had a 3.92-fold higher lung adenocarcinoma risk than the bottom 20%. Further, the top 5% PRS group reached a 6.69% lifetime absolute risk. Notably, this group reached the average population 10-year lung adenocarcinoma risk at age 50 years (0.42%) by age 41 years, 9 years earlier. Conclusions Multi-ancestry PRS approaches enhance lung adenocarcinoma risk stratification in East Asian never-smokers, with consistent external validation, suggesting future clinical utility.

Abstract Single-cell sequencing technologies provide us with information at the level of individual cells. Combining single-cell RNA-seq and single-cell TCR-seq profiling enables the exploration of cell heterogeneity and T-cell receptor repertoires simultaneously. Integrating both types of data can play a crucial role in enhancing our understanding of T-cell-mediated immunity and, in turn, facilitate the advancement of immunotherapy. Here, we present immunopipe, a comprehensive and flexible pipeline to perform integrated analysis of scRNA-seq and scTCR-seq data. In addition to the command line tool, we provide a user-friendly web interface for pipeline configuration and execution monitoring, benefiting researchers without extensive programming experience. With its comprehensive functionality and ease of use, immunopipe empowers researchers to uncover valuable insights from scRNA-seq and scTCR-seq data, ultimately advancing the understanding of immune responses and immunotherapy development.

Abstract Approximately two-thirds of hypopharyngeal squamous cell carcinoma (HPSCC) are diagnosed in advanced stage, with the worst prognosis among head and neck squamous cell carcinomas (HNSCCs). Identifying biomarkers for high-risk patients who need aggressive treatment is crucial. We present mutational, transcriptomic and proteomic studies of 103 Chinese HPSCC patients, observing a higher prevalence and poorer prognosis in males comparing to females. Genes in estrogen response pathways are up-regulated and proteins phosphorylated by protein kinase C (PKC) and Cyclin-dependent kinases (CDKs) are aberrantly regulated. We identify aberrant copy number regions including SOX2 (3q26.33), FGFR (8p11.23), CCND1 (11q13.3), CDKN2A/2B (9p21.3), and MYC (8q24.21). HPV status combined with highly mutated genes as predictive markers was assessed, including SYNE1 in HPV(-), MUC4 in HPV(+) patients. We confirm the presence of four known gene expression (GE) subtypes of HNSCC in our HPSCC cohort and present the genomic alterations, patient characteristics and survival status of each subtype. A predictive model involving clinical factors and transcriptomic expression of six genes was established and cross-site validated. The findings open new opportunities for stratifying high-risk patients and molecular targets for novel therapeutic strategies. Citation Format: Yanxin Ren, Wei Xiong, Chun Feng, Dan Yu, Xiaoyan Wang, Qing Yang, Siting Yu, Hongjiang Zhang, Bangyun Huo, Honglu Jiang, Zuli Li, Junlin Wang, Yu-Xiong Su, Ping Yang, Yong Liao, Qi Zhong, Junwen Wang. Multi-omics insights into the molecular signature and prognosis of hypopharyngeal squamous cell carcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 4999.

Abstract The interactions between transcription factors (TFs) and the target genes could provide a basis for constructing gene regulatory networks (GRNs) for mechanistic understanding of various biological complex processes. From gene expression data, particularly single-cell transcriptomic data containing rich cell-to-cell variations, it is highly desirable to infer TF–gene interactions (TGIs) using deep learning technologies. Numerous models or software including deep learning–based algorithms have been designed to identify transcriptional regulatory relationships between TFs and the downstream genes. However, these methods do not significantly improve predictions of TGIs due to some limitations regarding constructing underlying interactive structures linking regulatory components. In this study, we introduce a deep learning framework, DeepTGI, that encodes gene expression profiles from single-cell and/or bulk transcriptomic data and predicts TGIs with high accuracy. Our approach could fuse the features extracted from Auto-encoder with self-attention mechanism and other networks and could transform multihead attention modules to define representative features. By comparing it with other models or methods, DeepTGI exhibits its superiority to identify more potential TGIs and to reconstruct the GRNs and, therefore, could provide broader perspectives for discovery of more biological meaningful TGIs and for understanding transcriptional gene regulatory mechanisms.

Abstract Single-cell technologies enable researchers to investigate cell functions at an individual cell level and study cellular processes with higher resolution. Several multi-omics single-cell sequencing techniques have been developed to explore various aspects of cellular behavior. Using NEAT-seq as an example, this method simultaneously obtains three kinds of omics data for each cell: gene expression, chromatin accessibility, and protein expression of transcription factors (TFs). Consequently, NEAT-seq offers a more comprehensive understanding of cellular activities in multiple modalities. However, there is a lack of tools available for effectively integrating the three types of omics data. To address this gap, we propose a novel pipeline called MultiSC for the analysis of MULTIomic Single-Cell data. Our pipeline leverages a multimodal constraint autoencoder (single-cell hierarchical constraint autoencoder) to integrate the multi-omics data during the clustering process and a matrix factorization–based model (scMF) to predict target genes regulated by a TF. Moreover, we utilize multivariate linear regression models to predict gene regulatory networks from the multi-omics data. Additional functionalities, including differential expression, mediation analysis, and causal inference, are also incorporated into the MultiSC pipeline. Extensive experiments were conducted to evaluate the performance of MultiSC. The results demonstrate that our pipeline enables researchers to gain a comprehensive view of cell activities and gene regulatory networks by fully leveraging the potential of multiomics single-cell data. By employing MultiSC, researchers can effectively integrate and analyze diverse omics data types, enhancing their understanding of cellular processes.

Little is known about the expression patterns and functions of circular RNAs (circRNAs) in the heart of large mammals. In this study, we examined the expression profiles of circRNAs, microRNAs (miRNAs), and messenger RNAs (mRNAs) in neonatal pig hearts. Pig heart samples collected on postnatal days 1 (P1), 3 (P3), 7 (P7), and 28 (P28) were sent for total RNA sequencing. Our data revealed a total of 7,000 circRNAs in the 24 pig hearts. Pathway enrichment analysis of hallmark gene sets demonstrated that differentially expressed circRNAs were engaged in different pathways. The most significant difference was observed between P1 and the other 3 groups (P3, P7, and P28) in pathways related to cell cycle and muscle development. Out of the 10 circRNAs that were validated through real-time quantitative PCR to verify their expression, 6 exhibited significant effects on cell cycle activity in human induced pluripotent stem cell–derived cardiomyocytes following small interfering RNA–mediated knockdown. circRNA-miRNA-mRNA networks were constructed to understand the potential mechanisms of circRNAs in the heart. In conclusion, our study provided a data set for exploring the roles of circRNAs in pig hearts. In addition, we identified several circRNAs that regulate cardiomyocyte cell cycle.