Shi, Jianxin

Abstract Background: Epidemiologic and other studies have identified environmental, occupational, and genetic risk factors for lymphoid and myeloid malignancies, but most studies have been conducted in Western populations. Investigations in populations with differing exposure patterns, distributions of disease subtypes, and genetic architecture are needed to fully understand hematopoietic tumor etiology. Methods: We conducted a large, multicenter, hospital-based case-control study of lymphoid and myeloid neoplasms in East Asia (AsiaLymph), including 5,671 lymphoid cases, 1,879 myeloid cases, and 3,858 controls. Participants completed a computer-assisted personal interview and provided biospecimens. Cases underwent central pathology review and were coded into the WHO Classification. Herein, we describe the study methods in detail and report association results for education, body mass index, and family history. Results: Greater BMI at age 20 but not at age 40 was associated with odds of total lymphoid neoplasm (per 5 kg/m2 increase: OR [95% CI]: 1.18 [1.09-1.27]), total myeloid neoplasm (OR [95% CI]: 1.17 [1.05-1.29]), and several subtypes. Greater educational attainment was associated with increased odds of total lymphoid neoplasm (for college vs. less than primary education: OR [95% CI]: 1.41 [1.21-1.65]) but not myeloid neoplasm (OR [95% CI]: 1.17 [0.95-1.45]; p-heterogeneity=0.02). There were also positive associations between family history of hematologic cancer in first degree relatives and odds of lymphoid neoplasm (OR [95% CI]: 1.53 [1.15-2.03]) and myeloid neoplasm (OR [95% CI]: 1.65 [1.11-2.44]) and specific subtypes. None of the evaluated risk factors were associated with NK/T-cell lymphoma, which supports a distinct etiology for this subtype. Conclusion: The AsiaLymph Study is one of the largest molecular epidemiology studies of both lymphoid and myeloid neoplasms with standardized World Health Organization classification of histopathologic subtypes and will serve as a valuable resource for etiologic investigations into risk factors for these malignancies. Citation Format: Qing Lan, Lauren M. Hurwitz, John K. Chan, Tai Hing Lam, Kexin Chen, Yok Lam Kwong, Xu Caigang, Brian CH Chiu, Raymond Liang, Ip Dennis, Wei Hu, Bryan Bassig, Mark Purdue, Jun Xu, Sarah Locke, Sophia S. Wang, James R. Cerhan, Sonja Berndt, Jonathan N. Hofmann, Jianxin Shi, Kai Yu, Shahinaz Gadalla, Lisa J. McReynolds, Rena Jones, Hongji Dai, Zhangyan Lyu, Lugui Qiu, Wei Liu, Huilai Zhang, Xianhuo Wang, Lindsay M. Morton, Stephen Chanock, Martha Linet, Melissa C. Friesen, Roel Vermeulen, Nathaniel Rothman. A multi-center hospital-based case-control study of lymphoid and myeloid neoplasms (AsiaLymph): Study design and initial findings for education, body mass index, and family history [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2026; Part 1 (Regular Abstracts); 2026 Apr 17-22; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2026;86(7 Suppl):Abstract nr 5051.

Abstract Chromosome 5p15.33 harbors several independent association signals which demonstrate antagonistic pleiotropy across cancer types, with causal mechanisms largely unresolved. To identify functional variants and enhancer elements at this locus, we performed statistical fine-mapping followed by massively parallel reporter assays (MPRA) and proliferation based CRISPRi screens. This approach identified eight multi-cancer functional variants (MCFVs) across three GWAS signals. Targeting rs421629 (part of the CLPTM1L signal marked by rs465498) with CRISPRi revealed opposing effects on TERT expression in pancreatic versus lung cancer cells, consistent with the antagonistic pleiotropy observed for this signal. Furthermore, CRISPRi nominated an intronic CLPTM1L variable number tandem repeat (VNTR) as a potent enhancer. Long-read sequencing established VNTR polymorphisms as potential causal variants for the rs465498 signal. We showed that Hippo-pathway transcription factors mediate VNTR enhancer activity in lung and pancreatic cancer cells. Together, these findings indicate that cancer susceptibility at 5p15.33 may be mediated by both SNPs and VNTRs and provide an integrated framework for resolving complex pleiotropic loci.

Abstract Understanding tumor cell dynamics can improve prognosis and treatment but remains limited for lung adenocarcinoma in people who have never smoked (NS-LUAD). With RNA sequencing data from 684 NS-LUAD cases and validation in an independent dataset, we identified three subtypes with distinct phenotypic traits and cell compositions. Additional genomic and histologic data further characterized the subtypes. “Steady,” marked by low proliferation, high alveolar cell fraction, moderate-to-well differentiation, and fewer driver gene 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 histologic characteristics, this aggressive subtype is solely identified by transcriptomic data. A 60-gene signature recapitulates the classification and predicts survival even within subgroups based on tumor stage or known genomic features, emphasizing its potential for improving early-stage NS-LUAD prognostication in clinical settings. Significance: The transcriptome of 684 NS-LUAD identifies three subtypes with different cellular dynamics and genomic and morphologic features. A 60-gene signature accurately stratifies subjects for mortality risk, even in stage I, offering a potential clinically applicable tool for treatment decision-making in patients with NS-LUAD. See related commentary by Azizi et al., p. 423

ABSTRACT Understanding lung cancer evolution can identify tools for intercepting its growth. In a landscape analysis of 1024 lung adenocarcinomas (LUAD) with deep whole-genome sequencing integrated with multiomic data, we identified 542 LUAD that displayed diverse clonal architecture. In this group, we observed an interplay between mobile elements, endogenous and exogenous mutational processes, distinct driver genes, and epidemiological features. Our results revealed divergent evolutionary trajectories based on tobacco smoking exposure, ancestry, and sex. LUAD from smokers showed an abundance of tobacco-related C:G>A:T driver mutations in KRAS plus short subclonal diversification. LUAD in never smokers showed early occurrence of copy number alterations and EGFR mutations associated with SBS5 and SBS40a mutational signatures. Tumors harboring EGFR mutations exhibited long latency, particularly in females of European-ancestry (EU_N). In EU_N, EGFR mutations preceded the occurrence of other driver genes, including TP53 and RBM10. Tumors from Asian never smokers showed a short clonal evolution and presented with heterogeneous repetitive patterns for the inferred mutational order. Importantly, we found that the mutational signature ID2 is a marker of a previously unrecognized mechanism for LUAD evolution. Tumors with ID2 showed short latency and high L1 retrotransposon activity linked to L1 promoter demethylation. These tumors exhibited an aggressive phenotype, characterized by increased genomic instability, elevated hypoxia scores, low burden of neoantigens, propensity to develop metastasis, and poor overall survival. Re-activated L1 retrotransposition-induced mutagenesis can contribute to the origin of the mutational signature ID2, including through the regulation of the transcriptional factor ZNF695 , a member of the KZFP family. The complex nature of LUAD evolution creates both challenges and opportunities for screening and treatment plans.

Abstract Background People with cancer experience high rates of venous thromboembolism (VTE). Risk of subsequent cancer is also increased in people experiencing their first VTE. The causal mechanisms underlying this association are not completely understood, and it is unknown whether VTE is itself a risk factor for cancer. Methods We used data from large genome-wide association study meta-analyses to perform bidirectional Mendelian randomization analyses to estimate causal associations between genetic liability to VTE and risk of 18 different cancers. Results We found no conclusive evidence that genetic liability to VTE was causally associated with an increased incidence of cancer, or vice versa. We observed an association between liability to VTE and pancreatic cancer risk [odds ratio for pancreatic cancer: 1.23 (95% confidence interval: 1.08–1.40) per log-odds increase in VTE risk, P = 0.002]. However, sensitivity analyses revealed this association was predominantly driven by a variant proxying non-O blood group, with inadequate evidence to suggest a causal relationship. Conclusions These findings do not support the hypothesis that genetic liability to VTE is a cause of cancer. Existing observational epidemiological associations between VTE and cancer are therefore more likely to be driven by pathophysiological changes which occur in the setting of active cancer and anti-cancer treatments. Further work is required to explore and synthesize evidence for these mechanisms.

Abstract Genome wide association studies (GWAS) in 9,013 pancreatic cancer patients and 12,452 controls of European ancestry have discovered over 20 risk loci in the human genome. Here, we fine mapped one such locus on chromosome 16q23.1 in the vicinity of two Chymotrypsinogen precursor genes, CTRB1 and CTRB2. We fine-mapped this locus to rs72802365 (P=2.51x10-17, OR=1.36, 95% CI=1.31-1.40) and identified colocalization (PP=87%) with aberrant exon 5→7 CTRB2 splicing in pancreatic tissues. Imputation of a 584 bp insertion/deletion variant overlapping exon 6 of CTRB2 into the GWAS datasets resulted in a highly significant association with pancreatic cancer risk (P=2.83x10-16, OR=1.36, 95% CI=1.31-1.42), indicating that it may underlie this signal. Exon skipping attributable to the deletion (risk) allele introduces a premature stop codon in exon 7 of CTRB2, yielding a truncated chymotrypsin B2 protein that lacks chymotrypsin activity, is poorly secreted, and accumulates intracellularly in the endoplasmic reticulum (ER). The intracellular accumulation of a nonfunctional chymotrypsinogen B2 protein may lead to ER stress and pancreatic inflammation, which in turn, may explain the increased pancreatic cancer risk in carriers of CTRB2 exon 6 deletion alleles. Citation Format: Ashley Jermusyk, Jun Zhong, Naomi Gordon, Katelyn Connelly, Sumeth Perera, Ehssan Abdolalizadeh, Tongwu Zhang, Aidan O'Brien, Jason Hoskins, Irene Collins, Daina Eiser, Chen Yuan, Harvey Risch, Eric J. Jacobs, Donghui Li, Mengmeng Du, Rachael Stolzenberg-Solomon, Alison P. Klein, Jill P. Smith, Brian M. Wolpin, Stephen J. Chanock, Jianxin Shi, Gloria M. Petersen, Christopher Westlake, Laufey T. Amundadottir. An exception to the rule: A coding functional variant at a pancreatic cancer GWAS locus [abstract]. In: Proceedings of the AACR Virtual Special Conference on Pancreatic Cancer; 2020 Sep 29-30. Philadelphia (PA): AACR; Cancer Res 2020;80(22 Suppl):Abstract nr IA-11.

Abstract Registry-based epidemiologic studies suggest associations between chronic inflammatory intestinal diseases and pancreatic ductal adenocarcinoma (PDAC). As genetic susceptibility contributes to a large proportion of chronic inflammatory intestinal diseases, we hypothesize that the genomic regions surrounding established genome-wide associated variants for these chronic inflammatory diseases are associated with PDAC. We examined the association between PDAC and genomic regions (±500 kb) surrounding established common susceptibility variants for ulcerative colitis, Crohn's disease, inflammatory bowel disease, celiac disease, chronic pancreatitis, and primary sclerosing cholangitis. We analyzed summary statistics from genome-wide association studies data for 8,384 cases and 11,955 controls of European descent from two large consortium studies using the summary data-based adaptive rank truncated product method to examine the overall association of combined genomic regions for each inflammatory disease group. Combined genomic susceptibility regions for ulcerative colitis, Crohn disease, inflammatory bowel disease, and chronic pancreatitis were associated with PDAC at P values < 0.05 (0.0040, 0.0057, 0.011, and 3.4 × 10−6, respectively). After excluding the 20 PDAC susceptibility regions (±500 kb) previously identified by GWAS, the genomic regions for ulcerative colitis, Crohn disease, and inflammatory bowel disease remained associated with PDAC (P = 0.0029, 0.0057, and 0.0098, respectively). Genomic regions for celiac disease (P = 0.22) and primary sclerosing cholangitis (P = 0.078) were not associated with PDAC. Our results support the hypothesis that genomic regions surrounding variants associated with inflammatory intestinal diseases, particularly, ulcerative colitis, Crohn disease, inflammatory bowel disease, and chronic pancreatitis are associated with PDAC. Significance: The joint effects of common variants in genomic regions containing susceptibility loci for inflammatory bowel disease and chronic pancreatitis are associated with PDAC and may provide insights to understanding pancreatic cancer etiology.

A major direction for future microbiome research is implementation of fecal sample collections in large-scale, prospective epidemiologic studies. Studying microbiome-disease associations likely requires microbial data to be pooled from multiple studies. Our findings suggest collection methods that are most optimal to be used standardly across future WGSS microbiome studies.

Abstract Background Although 20 pancreatic cancer susceptibility loci have been identified through genome-wide association studies in individuals of European ancestry, much of its heritability remains unexplained and the genes responsible largely unknown. Methods To discover novel pancreatic cancer risk loci and possible causal genes, we performed a pancreatic cancer transcriptome-wide association study in Europeans using three approaches: FUSION, MetaXcan, and Summary-MulTiXcan. We integrated genome-wide association studies summary statistics from 9040 pancreatic cancer cases and 12 496 controls, with gene expression prediction models built using transcriptome data from histologically normal pancreatic tissue samples (NCI Laboratory of Translational Genomics [n = 95] and Genotype-Tissue Expression v7 [n = 174] datasets) and data from 48 different tissues (Genotype-Tissue Expression v7, n = 74–421 samples). Results We identified 25 genes whose genetically predicted expression was statistically significantly associated with pancreatic cancer risk (false discovery rate < .05), including 14 candidate genes at 11 novel loci (1p36.12: CELA3B; 9q31.1: SMC2, SMC2-AS1; 10q23.31: RP11-80H5.9; 12q13.13: SMUG1; 14q32.33: BTBD6; 15q23: HEXA; 15q26.1: RCCD1; 17q12: PNMT, CDK12, PGAP3; 17q22: SUPT4H1; 18q11.22: RP11-888D10.3; and 19p13.11: PGPEP1) and 11 at six known risk loci (5p15.33: TERT, CLPTM1L, ZDHHC11B; 7p14.1: INHBA; 9q34.2: ABO; 13q12.2: PDX1; 13q22.1: KLF5; and 16q23.1: WDR59, CFDP1, BCAR1, TMEM170A). The association for 12 of these genes (CELA3B, SMC2, and PNMT at novel risk loci and TERT, CLPTM1L, INHBA, ABO, PDX1, KLF5, WDR59, CFDP1, and BCAR1 at known loci) remained statistically significant after Bonferroni correction. Conclusions By integrating gene expression and genotype data, we identified novel pancreatic cancer risk loci and candidate functional genes that warrant further investigation.

Abstract Pancreatic Ductal Adenocarcinoma (PDAC) is the third leading cause of cancer-related deaths in United States with a 5-year survival rate of only 8%. Inherited predisposition plays an important role in PDAC risk. Rare, moderately to highly penetrant mutations in hereditary cancer and pancreatitis genes, identified in families with a high incidence of disease, account for a small fraction of cases. At the other end of the spectrum, common risk variants of low penetrance have been discovered through genome-wide association studies (GWAS). Currently, over 20 common PDAC risk signals have been identified through GWAS in populations of European ancestry. However, these loci explain only a small fraction of genetic heritability for PDAC, and the causal genes for the large majority of these are unknown. Most susceptibility alleles discovered through GWAS are located in noncoding regions of the genome and are likely to function through allele specific regulation of gene expression. Transcriptome-wide association studies integrate expression quantitative trait loci (eQTL) with GWAS data to discover genes whose cis-regulated expression is associated with complex traits. To identify novel susceptibility loci and likely causal genes, we applied this approach to PDAC. Firstly, we used two reference datasets generated from histologically normal pancreatic tissue samples from individuals of European ancestry, the DCEG/LTG (n=95) and Genotype-Tissue Expression (GTEx v7, n=174), to build robust predictors for gene expression (n=2,872-5,906; r2>0.01), using 5 different models (Top1, Lasso, Enet, Blup and Bslmm) in FUSION and PrediXcan. Gene expression predictors (n=2,043-21,422) were also generated by Enet in 48 different tissues (n=74-421 samples) from 608 individuals of European ancestry (GTEx). We then applied these predictor models to PDAC GWAS summary statistics from 9,040 PC cases and 12,496 controls. We discovered 24 genes whose genetically predicted expression is significantly associated with PDAC risk (FDR < 0.05), including 13 candidate genes at 10 novel loci (CELA3B, SMC2, SMC2-AS1, RP11-80H5.9, SMUG1, BTBD6, RCCD1, PNMT, CDK12, PGAP3, SUPT4H1, RP11-888D10.3 and PGPEP1) and 11 at 5 known risk loci (TERT, CLPTM1L, ZDHHC11B, INHBA, ABO, PDX1, KLF5, WDR59, CFDP1, BCAR1 and TMEM170A). Pathway enrichment analysis for genes differentially expressed in the top vs. bottom quartile of CELA3B expression highlighted higher expression of inflammatory and immune response genes in both the DCEG/LTG and GTEx pancreatic transcriptome datasets, indicating that low expression of this gene may be associated with an inflammatory state in the pancreas. Our findings showcase the utility of integrating multiple datasets to generate hypotheses and provide novel clues for the genetic mechanisms underlying pancreatic cancer risk. Citation Format: Jun Zhong, Ashley Jermusyk, Lang Wu, Jason W. Hoskins, Irene Collins, Mingfeng Zhang, Song Lei, Charles C. Chung, Tongwu Zhang, Wenming Xiao, Rachael Stolzenberg-Solomon, Alison P. Klein, Brian M. Wolpin, Xiao-Ou Shu, Stephen J. Chanock, Sara Olson, Nilanjan Chatterjee, Jill Smith, Jianxin Shi, Peter Kraft, Gloria M. Petersen, Wei Zheng, Laufey T. Amundadottir. Large-scale transcriptome-wide association study (TWAS) identifies novel candidate susceptibility genes for pancreatic cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 1591.