Search PubMedSearch

Biomedical subjects

Tongwu Zhang

Publications and source records attributed to Tongwu Zhang.

10 recordsLinked to original sources

Functional characterization of the 9q34.13 locus identifies RAPGEF1 as a candidate gene modulating risk for melanoma and nevi via RAS activation.

Genome-wide association studies identified a melanoma- and nevus count-associated locus on chromosome band 9q34.13. Fine-mapping and melanocyte expression data collectively suggest two potential risk genes with opposite associations with risk: higher levels of Rap guanine nucleotide exchange factor 1 (RAPGEF1) and lower levels of uridine-cytidine kinase 1 (UCK1). Colocalization analyses and conditional transcriptome-wide association studies (TWASs) suggest multiple causal cis-regulatory sequence variants in partial linkage disequilibrium (LD) to each other. Melanocyte capture-HiC and CRISPR inhibition demonstrated regulatory interactions between fine-mapped variants and the RAPGEF1 and UCK1 promoters. Focusing on RAPGEF1, we demonstrate that RAPGEF1 expression promotes melanocyte growth and drives colony formation of human immortalized melanocytes. Following treatment with human epidermal growth factor (EGF), RAPGEF1 overexpression activated both RAP1 and RAS. Further, we show that RAPGEF1 expression is significantly enriched in melanomas that lack strongly activating RAS-MAPK pathway mutations, which suggests that RAPGEF1 may promote oncogenic RAS-MAPK pathway signaling in melanomas. Furthermore, in these tumors, we provide preliminary evidence to support the prognostic relevance of RAPGEF1 expression in individuals whose melanomas lack RAS or BRAF mutations. Together with other recent studies, these data suggest that germline variation influencing RAS activation may play a key role in nevus development and melanoma risk.

GWAS

Identification of immune cell type-specific susceptibility genes in multiple cancers using transcriptome-wide association studies.

BACKGROUND: Transcriptome-wide association studies (TWAS) integrate gene expression and genome-wide association studies (GWAS) to identify disease susceptibility genes. Because gene expression varies substantially across cell types within tissues, cell type-specific prediction models may enhance the power of TWAS. METHODS: We conducted cell type-specific TWAS leveraging single-cell RNA sequencing data from the OneK1K cohort (14 immune cell types, 1.27 million cells) and GWAS summary statistics for 7 cancers (>290 000 cases in total). To improve prediction accuracy, we developed a modeling framework that incorporates shared gene expression effects across cell types. RESULTS: At a false discovery rate of 5%, we identified 106 (Bonferroni 5%: 13) previously unreported loci for breast cancer, 51 (4) loci for prostate cancer, 11 (4) loci for lung cancer, 39 (5) loci for melanoma, 9 (1) loci for ovarian cancer, and 2 (1) loci for diffuse large B-cell lymphoma, with most genes exhibiting cell type specificity. Gene set analyses confirmed joint associations of unreported genes with breast and prostate cancer risk in UK Biobank data. Additional lung tissue single-cell RNA sequencing data with 113 individuals validated 18 of 32 (56.3%) statistically significant genes for lung cancer. Across cancers, 139 statistically significant genes were shared by at least 2 cancer types and were primarily enriched in specific immune cell types. CONCLUSION: Cell type-specific TWAS improve the identification of novel cancer susceptibility loci and provide insights into the immune landscape of cancer etiology.

Humans

Detection of House Dust Mite-derived DNA in Human Lung Tumors by Whole-Genome Sequencing.

Lung cancer in never-smokers (LCINS) accounts for an increasing proportion of lung cancer cases, yet its risk factors remain poorly understood. House dust mites (HDM) are common aeroallergens that induce airway inflammation, but their potential contribution to lung cancer is unknown. We analyzed unmapped whole-genome sequencing reads from 783 lung cancers from the Sherlock-Lung (n = 621 never-smokers) and EAGLE (n = 162 smokers) cohorts, including 328 matched adjacent normal lung tissues. After removal of human sequences, reads were aligned to reference genomes from the two major HDM species and confirmed by BLAST. Samples with top BLAST matches were classified as HDM-detected. Associations between HDM detection and genomic, microbiome, and bulk RNA-seq-derived immune features were evaluated. HDM-derived DNA was detected at low abundance in a subset of tumors and adjacent normal tissues, with higher detection frequencies in tumors than matched normal tissues and in smokers than never-smokers. In LCINS tumors, HDM detection was not associated with tumor mutational burden or recurrent driver alterations but was associated with modest differences in immune cell composition and a limited but reproducible bacterial co-detection pattern. These findings provide a foundation for investigating aeroallergen-derived DNA signatures and their potential relationship to the lung tumor microenvironment.

Environmental exposure

Methylation profiling of normal tissue adjacent to breast tumors reveals two distinct groups with divergent tumor microenvironment features.

We previously identified diverse genetic evolutionary patterns in whole-genome sequencing of paired normal tissue adjacent to tumor (NAT) and tumor tissues from Hong Kong breast cancer (HKBC) patients. Here, we investigated whether DNA methylation (DNAm) contributes to NAT heterogeneity and shapes the tumor microenvironment (TME). Genome-wide DNAm profiling was performed on paired NAT and tumor tissues from 188 HKBC patients using the Infinium 850 K array. RNA-seq data were available for 76 NATs and 177 tumors. Cellular composition was inferred using MethylCIBERSORT, CIBERSORTx, and EpiDISH, and histopathologic features were assessed on 115 H&E-stained sections. Unsupervised clustering identified two distinct NAT subtypes with divergent TME characteristics. Cluster 1 (N = 139) showed higher epithelial and fibroblast content and enrichment of estrogen response pathways. Cluster 2 (N = 49) exhibited an immune-metabolic phenotype characterized by increased fat and immune cells, stromal disruption, inflammatory pathway activation, and greater macrophage infiltration. Cluster 2 patients also demonstrated significantly younger epigenetic age estimated using multiple epigenetic clocks. These DNAm-defined NAT subtypes and associated TME features were validated in 97 NAT samples from TCGA breast cancer patients. Overall, our findings identify DNAm-driven NAT heterogeneity with distinct TME landscapes, providing new insights into field cancerization and tumor evolution in breast cancer.

Journal Article

Disruption of CTCF binding by germline non-coding variants in CDKN2B suppress CDKN2A expression and predispose to melanoma.

Some melanoma-prone families linked to the 9p21 locus, harboring the established susceptibility gene CDKN2A, lack pathogenic protein-coding variants. Using whole-exome and targeted sequencing, we identified three rare single-nucleotide variants in two melanoma-prone families and one sporadic melanoma case. Variants map to a conserved CTCF-bound region within the first intron of CDKN2B that physically interacts with CDKN2A. Analysis of UK Biobank showed significant enrichment of variants in this region in melanoma cases. Variants result in diminished CTCF binding in vitro. CTCF ChIP-seq in fibroblasts from the carriers of the largest family demonstrated loss of CTCF binding, accompanied by weakened promoter interactions and allele-specific reduction of CDKN2A p16 transcript expression from the variant haplotype. CRISPR-based perturbation of this region and editing of the large family variant into melanocytes resulted in reduced expression of p14 and p16 CDKN2A transcripts. These findings suggest that non-coding regulatory variants function as high-penetrance susceptibility alleles in melanoma families by altering CDKN2A function.

Journal Article

Functional characterization of the 9q34.13 locus identifies RAPGEF1 as modulating risk for melanoma and nevi via RAS activation.

Genome-wide association studies identified a melanoma- and nevus count-associated locus on chromosome band 9q34.13. Fine-mapping and melanocyte expression data collectively suggest two potential causal genes with opposite association with risk: higher levels of Rap guanine nucleotide exchange factor 1 (RAPGEF1) and lower levels of uridine-cytidine kinase 1 (UCK1). Colocalization analyses and conditional TWAS suggest multiple causal cis-regulatory sequence variants in partial linkage disequilibrium (LD) to each other. Melanocyte capture-HiC and CRISPR-inhibition demonstrated regulatory interactions between fine-mapped variants and the RAPGEF1 and UCK1 promoters. Focusing on RAPGEF1, we demonstrate RAPGEF1 expression promotes melanocyte growth and drives malignant transformation of human immortalized melanocytes. Following treatment with human EGF, RAPGEF1 overexpression activated both RAP1 and RAS. Further, we show RAPGEF1 expression is significantly enriched in melanomas lacking strongly activating RAS-MAPK mutations, suggesting that RAPGEF1 may promote oncogenic RAS-MAPK signaling in melanomas. Furthermore, in these tumors, we provide preliminary evidence to support the prognostic relevance of RAPGEF1 expression in patients lacking RAS or BRAF mutations. Together with other recent studies, these data suggest that germline variation influencing RAS activation may play a key role in nevus development and melanoma risk.

Journal Article

STK11 Mutations and Deletions Define an Aggressive Molecular Subgroup of Cervical Adenocarcinoma.

Cervical adenocarcinoma accounts for 15%-20% of cervical cancers and is associated with poorer survival and reduced response to screening and immunotherapy compared with squamous cell carcinoma (SCC). The genomic drivers underlying this molecular subgroup remain incompletely characterized. Whole-exome sequencing was performed on 302 invasive cervical cancers from Guatemala and Venezuela. Structural variation analysis was conducted using SNP-array and whole-genome sequencing data. Findings were replicated in more than 4600 additional cervical cancer samples from TCGA, AACR Project GENIE, MSKCC, and Caris datasets. TP53 mutations were more frequent in adenocarcinoma than SCC, particularly in HPV-negative tumors. STK11 alterations, including mutations and focal deletions, were significantly enriched in HPV-positive adenocarcinomas compared with SCC and affected 23% of adenocarcinomas overall. Whole-genome analyses identified recurrent focal deletions, inversions, chromosomal rearrangements, and breakage-fusion-bridge events involving chromosome 19p and STK11 that were not detected by exome sequencing alone. STK11 alterations were associated with younger age at diagnosis, poorer overall survival, and inferior outcomes following immune checkpoint inhibitor (ICI) therapy. STK11 alterations significantly co-occurred with YAP1 amplification but were largely mutually exclusive with PIK3CA mutation. Cervical adenocarcinomas also demonstrated significantly lower CD274 (PD-L1) expression than SCC. STK11 alterations define a distinct molecular subgroup of cervical adenocarcinoma characterized by structural disruption of chromosome 19p, younger age at onset, and poorer clinical outcomes. These findings have implications for molecular classification and future targeted therapeutic approaches in cervical cancer.

Humans

Panorama of Chromosomal Instability in Lung Cancer.

Lung cancer is a highly heterogeneous disease primarily driven by tobacco smoking. About 20% of lung cancers occur among patients who have never smoked (LCINS) with differences in patient ancestry, sex, tumor histology, and clinical features. Our understanding of chromosomal instability in lung cancer, especially LCINS, is still limited. Here, we perform a comprehensive study of 182,429 somatic structural variations (SVs) detected in 1,209 whole-genome sequenced lung cancers, of which 864 LCINS. SVs are more abundant in tumors from patients who have smoked (LCSS); however, they are more complex and play more important roles in tumorigenesis in LCINS. EGFR mutations and KRAS mutations profoundly and independently shape the SV landscape. EGFR-mutant tumors have higher SV burden and more cancer-driving SVs. In contrast, KRAS mutations are associated with lower SV burden and less driver SVs. We decompose 16 SV signatures for both complex and simple SVs that likely represent divergent molecular mechanisms. The SV breakpoints have distinct distributions across the genome depending on the signatures due to mutagenic mechanisms and positive selection. Many established cancer-driving genes are recurrently rearranged by multiple SV signatures suggesting functional convergence of these genome instability mechanisms.

Journal Article

Stratifying Lung Adenocarcinoma Risk with Multi-ancestry Polygenic Risk Scores in East Asian Never-Smokers.

BACKGROUND: Lung adenocarcinoma (LUAD) 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 LUAD risk prediction. METHODS: PRSs were developed using genome-wide association study summary statistics from East Asian (8,002 cases; 20,782 controls) and European (2,058 cases; 5,575 controls) populations. Single-ancestry models included PRS-25, PRS-CT, and LDpred2; multi-ancestry models included LDpred2+PRS-EUR128, PRS-CSx, and CT-SLEB. Performance was evaluated in independent East Asian data from the Female Lung Cancer Consortium (FLCCA) and externally validated in the Nanjing Lung Cancer Cohort (NJLCC). We assessed predictive accuracy via AUC, with 10-year and (age 30-80) absolute risks estimates. RESULTS: The best multi-ancestry PRS, using East Asian and European data via CT-SLEB (clumping and thresholding, super learning, empirical Bayes), outperformed the best East Asian-only PRS (LDpred2; AUC=0.629, 95% CI:0.618,0.641), achieving an AUC of 0.640 (95% CI:0.629,0.653) and odds ratio of 1.71 (95% CI:1.61,1.82) per SD increase. NJLCC Validation confirmed robust performance (AUC =0.649, 95% CI: 0.623, 0.676). The top 20% PRS group had a 3.92-fold higher LUAD 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 LUAD risk at age 50 (0.42%) by age 41, nine years earlier. CONCLUSIONS: Multi-ancestry PRS approaches enhance LUAD risk stratification in East Asian never-smokers, with consistent external validation, suggesting future clinical utility.

East Asian never smokers

A prognostic signature for lung adenocarcinoma in people who have never smoked.

Knowledge of tumor cell dynamics can inform prognosis and treatment yet is largely lacking for lung adenocarcinoma in people who have never smoked (NS-LUAD). With RNA-seq data from 684 NS-LUAD and validation in an independent dataset, we identified three subtypes with distinct phenotypic traits and cell compositions. Additional genomic and histological data further characterized the subtypes. 'Steady', marked by low proliferation, high alveolar cell fraction, moderate-to-well differentiation, and fewer driver genes' alterations, is linked to prolonged survival and low immune evasion. 'Proliferative' shows high proliferation markers, TP53 mutations, and gene fusions. 'Chaotic', with high epithelial-to-mesenchymal transition markers, has the worst prognosis even within stage I tumors. Lacking known molecular or histological characteristics, this aggressive subtype is solely identified by transcriptomic data. A 60-gene signature recapitulates the overall classification and strongly predicts survival even within subgroups based on tumor stage or known genomic features, emphasizing its potential for improving NS-LUAD prognostication in clinical settings.

Journal Article