Zhu, Qianqian

Abstract Background Breast cancer genome-wide association studies (GWAS) have identified more than 200 susceptibility loci, but most studies are dominated by European and East Asian populations. Methods We analyzed breast cancer GWAS summary statistics from African (AFR), East Asian (EAS), European (EUR), and Hispanic/Latina (H/L) samples (159,297 cases and 212,102 controls). We estimated logit-scale SNP-based heritability, polygenicity, and cross-ancestry genetic correlation, partitioned heritability across functional annotations, and integrated GWAS results with the Tabula Sapiens single-cell atlas using scDRS+. Results The logit-scale heritability of breast cancer ranged from ℎ 2 =0.47 (SE = 0.07) in EAS to AFR ℎ 2 =0.61 (SE = 0.10), with no significant differences across ancestries (p=0.63). The estimated number of susceptibility markers in a sparse normal-mixture effects model also varied from 4,446 (SE = 3,100) in EAS to 8,308 (SE = 2,751) in AFR, but differences were not significant across ancestries (p=0.55). Cross-sample genetic correlations varied, with the strongest correlation between EUR and EAS (𝜌 = 0.79, SE = 0.08) and weakest between AFR and H/L (𝜌 = 0.26, SE = 0.24). Regulatory annotations were enriched for breast cancer heritability across samples. Integration with single-cell expression profiles implicated ancestry-shared associations with innate immune, secretory epithelial, and stromal cell types. Conclusion These results indicate substantial cross-ancestry sharing of breast cancer polygenic architecture, highlight a consistent contribution of regulatory variation, and identify convergent cellular contexts that motivate functional follow-up and inform expectations for the transferability and attainable performance of common-variant risk prediction across populations.

Abstract Introduction: Breast cancer genome-wide association studies (GWAS) have identified over 200 susceptibility loci, many replicated in diverse populations. However, cross-ancestry evaluation of breast cancer genetic architecture remains limited. We examined breast cancer genetic architecture using GWAS summary results from European (EUR; cases (ca) = 188,474, controls (co) = 96,201), East Asian (EAS; ca = 20,393, co = 86,329), African American (AA; ca = 9,235, co = 10,184), and US Hispanic/Latina and Latin American (H/L; ca = 2,396, co =7,468) studies. Methods: GWAS results were derived from the Breast Cancer Association Consortium, the African Ancestry Breast Cancer Genetics Consortium (AABCG), Biobank Japan, and a meta-analysis of five studies of H/L women with breast cancer. Linkage disequilibrium (LD) scores were generated for EUR, EAS, and AMR using the 1000 Genomes Project, while AA LD scores were derived from a subset of AABCG controls. Heritability was estimated for each ancestry group using LD score regression. Genetic correlations between populations were assessed via Popcorn. Polygenicity was analyzed using GENESIS; however, due to the limited sample size in the H/L studies, these estimates were restricted to AA, EAS, and EUR. Enrichment of heritability across regulatory elements was evaluated using stratified LD score regression across all populations. Results: The logit-scale heritability and corresponding standard error (SE) were 0.466 (0.066) for EAS, 0.501 (0.050) for EUR, 0.588 (0.360) for H/L, and 0.614 (0.095) for AA. The estimated number of independent susceptibility loci was 4,446 for EAS, 5,235 for EUR, and 8,308 for AA. Using a clumping and thresholding approach, an optimal set of common variants were projected to explain 38.6% (EUR), 39.4% (EAS), and 26.2% (AA) of genetic variance for samples of 100,000 cases and 100,000 controls with AUCs from corresponding estimated polygenic risk scores (PRSs) of 0.621 (EUR), 0.634 (EAS), and 0.611 (AA). Genetic correlations were strongest between EUR and EAS (ρ = 0.79, SE = 0.08) and EUR and H/L (ρ = 0.68, SE = 0.21), and weakest between AA and EAS (ρ = 0.42, SE = 0.14) and AA and H/L (ρ = 0.26, SE = 0.24). Among 73 genomic features, we found significant (P < 0.05/73) enrichment in heritability in EUR for ‘ancient promoter’ regions, transcription factor binding sites, H3K4me3, ‘super-enhancers’, H3K4me1, H3K27ac; in EAS for ‘super-enhancers’; and in AA for H3K27ac. No genomic features were significantly enriched in H/L, likely due to limited power. Conclusion: These findings suggest a shared breast cancer genetic architecture across diverse populations, as well as the potential for a similar level of breast cancer risk stratification by PRS in these populations. Expanding GWAS in underrepresented populations is essential to improve genetic risk predictions and foster equitable cancer prevention. Citation Format: James L. Li, Maria Zanti, Jacob Williams, Om Jahagirdar, Guochong Jia, Qiang Hu, Jean-Tristan Brandenburg, Li Yan, Weang-Kee Ho, Jingmei Li, José P. Miranda, Devika Godbole, Julie-Alexia Dias, Leila Dorling, Wenlong C. Chen, Nicholas Boddicker, Ying Wang, Alicia Martin, Martin J. Zhang, Yan Zhang, Joe Dennis, Esther M. John, Gabriela Torres-Mejia, Larry Kushi, Jeffrey Weitzel, Susan L. Neuhausen, Luis Carvajal-Carmona, Christopher Haiman, Elad Ziv, Laura Fejerman, Wei Zheng, Dezheng Huo, Douglas Easton, Nilanjan Chatterjee, Peter Kraft, Montserrat Garcia-Closas, Wendy Wong, Kyriaki Michailidou, Qianqian Zhu, Diptavo Dutta, Thomas U. Ahearn, Haoyu Zhang. Genetic architecture of breast cancer across diverse populations: Assessing heritability, genetic correlation, and polygenicity [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 2282.

Abstract The first two authors and last two authors contributed equally. Genome-wide association studies (GWAS) have identified risk loci for Acute Lymphoblastic Leukemia (ALL), Chronic Lymphoblastic Leukemia (CLL) and Non-Hodgkin Lymphoma, however an Acute Myeloid Leukemia (AML) GWAS has not been published to date. We performed a GWAS to identify AML and Myelodysplastic Syndrome (MDS) risk loci using a nested case-control study design in the DISCOVeRY-BMT cohorts which includes almost 2000 AML and MDS patients as cases and 2813 unrelated donors as controls. Genotyping was performed using the Illumina Human OmniExpress BeadChip and imputed using the Haplotype Reference Consortium, yielding > 8 million high-quality variants for statistical analysis. Logistic regression models with AML (de novo AML with normal cytogenetics, de novo AML with abnormal cytogenetics and therapy-related AML) and MDS cases and European American healthy donor controls adjusted for age and sex were used to test the association of each SNP with disease status. To identify the strongest association signal with disease we conducted a summary statistic SNP-based association analysis (ASSET) using non-overlapping AML and MDS cases implemented in R statistical software. ASSET uses an exhaustive search for SNPs with small but common pleiotropic effects across groups of traits while accounting for the multiple tests required by the subset search, as well as any shared controls between groups. This approach allowed us to further investigate the heterogeneity within AML subtypes and to gain increased power by pooling subtypes that show pleiotropic effects. ASSET genome wide (GW) significance is defined as P<5.0x10-8, however, we present results of each subset analysis. ASSET analyses identified an association of T allele at rs12203592, Interferon Regulatory Factor 4 (IRF4), with increased risk of de novo AML and MDS (Figure 1) (OR=1.37; 95% CI, 1.25-1.51, Pmeta=1.59x10-10). The variant has been reported significant in GWAS of skin pigmentation, hair color, non-melanoma skin cancer, squamous cell carcinoma, actinic keratosis, and childhood ALL; GW case-control studies of neuroblastoma and breast cancer show suggestive association signals (P<5 x 10-6). IRF4 belongs to the IRF family of transcription factors and is a key regulator of differentiation stages in hematopoiesis. The T allele at rs12203592 is associated with significantly increased expression of IRF4 in whole blood and lung tissue and in vitro mouse studies have shown the T allele to be associated with higher levels of IRF4 transcription. rs12203592 is <80bp from an IRF4 transcription start site and in an important position within NF-κB motifs in multiple blood and immune cell lines, supporting the hypothesis that this SNP modulates NF-κB repression of IRF4 expression. Another variant in the IRF4 regulatory region, rs62389423 (A allele), showed a putative association with subsets of de novo AML and MDS (OR=1.36; 95% CI,1.21-1.52, Pmeta=1.2x10-7). Although not in linkage disequilibrium with rs12203592, the A allele has also been previously associated with susceptibility to skin cancer and melanoma in multiple GWAS. Analysis of 596 UKbiobank CLL cases and >300,000 controls, shows the A allele correlates with a 50% increased risk of CLL (95% CI, 1.32-1.73, P= 2.4x10-9).Several other GWAS also show an association between additional IRF4 variants and CLL. Analyses in individual subtypes also revealed an intronic variant, rs10098598, in RAD21 to be associated with de novo AML (OR=2.18; 95% CI, 1.60-2.97; Pmeta=8.9x10-7). RAD21, on8q24, a region known to be associated with multiple cancer types, is involved in DNA double-stranded break repair and in chromatid cohesion in mitotic cells. Loss of function would theoretically lead to chromosomal instability and tumorigenesis. Components of the cohesion complex (including RAD21) are somatically mutated in ~12% of de novo AML and MDS patients, where it is frequently an early event. We provide the first GW evidence of association between a common variant and AML susceptibility. This SNP has been shown to be associated with multiple phenotypes, suggesting that there are pleiotropic effects at work. Our RAD21 finding is consistent with the role of cohesion in leukemogenesis and provides some evidence of its role in de novo AML susceptibility. Replication and further subset analyses with genome-wide data in >2000 AML cases are ongoing. Disclosures Griffiths: Alexion Inc.: Honoraria, Research Funding; Astex/Otsuka Pharmaceuticals: Honoraria, Research Funding; Celgene, Inc: Honoraria, Research Funding; Novartis, Inc.: Research Funding; Pfizer, Inc.: Research Funding. McCarthy:Bristol Myers Squibb: Consultancy, Honoraria; Celgene: Consultancy, Honoraria, Research Funding; Karyopharm: Consultancy, Honoraria; Takeda: Consultancy, Honoraria; Janssen: Consultancy, Honoraria.

Abstract Both genome wide association studies (GWAS) of common variation and exome wide association studies (EXWAS) of rare variation have successfully identified disease susceptibility variants for a variety of diseases. One GWAS of inherited susceptibility to Acute Myeloid Leukemia (AML) has been conducted, but no EXWAS have been performed to measure risk of AML attributable to low-frequency constitutional genetic variation. We performed the first EXWAS of risk of AML as a nested case-control study in the DISCOVeRY-BMT (Determining the Influence of Susceptibility Conveying Variants Related to one-Year mortality after BMT) cohorts. The DISCOVeRY-BMT parent study examined transplant-related mortality in leukemia patients undergoing unrelated donor allogeneic BMT. To identify low frequency variants and genes contributing to increased susceptibility to AML we used genotype data from the Illumina HumanExome BeadChip typed in the DISCOVeRY-BMT cohorts; the HumanExome BeadChip contains 242,901 variants, which are mainly protein-coding variants. The optimal sequence kernel association test (SKAT-O) was used to analyze gene-level associations with risk of AML. These gene-based tests evaluate the cumulative effects of multiple single gene variants on risk of AML. Analyses were performed in all European American AML cases and two subtypes: 1) de novo AML, 2) de novo AML with normal cytogenetics. Models were adjusted for age at transplant and principal components to control for population stratification. For gene-based tests at least 2 variants with minor allele frequency (MAF) ≤ 5%, were required to be present in the gene. This yielded a total of 13,687 genes tested, and a Bonferroni corrected significance level of P<3.65 x 10-6. Association tests were performed in 1,189 AML cases reported to CIBMTR 2000-08 (Cohort 1) and 327 AML cases reported to CIBMTR from 2009-11 (Cohort 2). Controls in Cohorts 1 (n=1,986) and 2 (n= 515) were 10/10 HLA-matched unrelated donors who passed a comprehensive medical exam and deemed healthy. We used metaSKAT to combine Cohorts 1 and 2 and obtain p-values of association with AML. We present the results of gene-level tests significant in both cohorts. The likely pathogenicity of these variants was determined in silico using SIFT, PolyPhen and MutationTaster. Patient characteristics are in Table 1. DNMT3A, on chromosome 2, was associated in the gene-based test with risk of AML (Pmeta=1.70x10-9, Table 2). Three missense variants at MAF <1% comprise both overall AML and de novo AML gene-based association: exm177559 (Asn->Ser), exm177507 (Arg->His), and exm177543 (Arg->Trp). Normal cytogenetics de novo AML gene-based assocations consisted of only 2 of these variants: exm177559 and exm177507 (Table 2). While prevalence of exm177507 is <1% for all AML cases, in de novo AML with normal cytogenetics the MAF was higher at 3%. The other 2 variants had a MAF<1% irrespective of subtype. Somatically, DNMT3A is most frequently mutated in hematologic malignancies, with >30% of de novo AML cases with a normal karyotype and >10% of MDS patients having DNMT3A mutations. Although these are germline gene associations all three of the variants found have been reported somatically in hematologic malignancies. In 200 AML cases from The Cancer Genome Atlas (TCGA) p.R882H (represented as exm177507 on the exome chip) was a frequent somatic mutation (25%). Exm177543 (p.R635W) and exm177559 (p.N501S) are reported in the Catalogue of Somatic Mutations in Cancer (COSMIC) as somatic mutations involved in hematopoietic and lymphoid tissue in both cell lines and humans. Exm177507 and exm177543 show evidence of pathogenicity in all three in silico tools, while exm177559 was reported as deleterious and disease causing by Sift and MutationTaster, respectively. Our results show that multiple potentially pathogenic missense germline variants in DNMT3A comprise the gene-based association with AML, specifically de novo AML with normal cytogenetics. Given the functional nature of these variants it is possible germline risk stratification could be informative in determining AML risk, and subsequently development of AML harboring DNMT3A mutations. Confirmation of these findings in additional cohorts could have implications for individualized risk screening, prediction and prognosis. Additional cytogenetic subgroup analyses, including treatment-related AML, are underway. Disclosures Hahn: Novartis: Equity Ownership; NIH: Research Funding. McCarthy:Sanofi: Honoraria, Membership on an entity's Board of Directors or advisory committees; Celgene: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; Janssen: Honoraria, Membership on an entity's Board of Directors or advisory committees; Onyx: Honoraria, Membership on an entity's Board of Directors or advisory committees; Bristol Myers Squibb: Honoraria, Membership on an entity's Board of Directors or advisory committees; The Binding Site: Honoraria, Membership on an entity's Board of Directors or advisory committees; Karyopharm: Honoraria, Membership on an entity's Board of Directors or advisory committees; Gamida Cell: Honoraria, Membership on an entity's Board of Directors or advisory committees. Sucheston-Campbell:NIH/NCI: Research Funding.

Abstract #the authors contributed equally to this work While survival outcomes after HLA-matched unrelated donor (URD) blood and marrow transplant (BMT) have significantly improved over the last two decades, about 40% of patients die before 1-year post URD BMT. Previously we performed a genome-wide association study (GWAS) named DISCOVeRY-BMT (Determining the Influence of Susceptibility COnveying Variants Related to one-Year mortality after BMT) to investigate the contributions of common genetic variants on survival outcomes. To address the specific contributions of low frequency (≤ 1%) exonic variants on survival outcomes, we used the DISCOVeRY-BMT cohorts typed with the Illumina HumanExome BeadChip containing 242,901 variants, which are mainly protein-coding variants. We used a gene-based test (the optimal sequence kernel association test (SKAT-O)) to evaluate the cumulative effects of multiple genetic variants within a gene on overall survival (OS), transplant-related mortality (TRM), and disease-related mortality (DRM) 1 year after URD BMT. SKAT-O measures the overall significance of the gene and allows the variants to have different directions and magnitude of effects. We used a 2-stage study design: SKAT-O analysis was carried out in Cohort 1 consisting of 1,972 AML, ALL or MDS patients reported to the Center for International Blood and Marrow Transplant Research from 2000 to 2008 and their 2,006 10/10 HLA-matched URD. Genes with a suggestive association with each outcome (P < 1×10-3) entered the final meta-analysis including Cohort 1 and an independent cohort (Cohort 2), consisting of 503 recipients and 520 10/10 HLA-matched URD from 2009 to 2011. The analysis included recipients and URD of European descent due to the low frequency of other populations. All association models included age at BMT, diagnosis (AML, ALL, MDS), disease status at BMT, cell source (peripheral blood, marrow), year of BMT, and principal components from EIGENSTRAT to control for population stratification. Only coding variants (missense and nonsense) with minor allele frequency ≤ 1% were included in our analysis, which resulted in > 12,000 genes tested in our gene-based analysis in recipients and donors. Exome-wide significance was set at Pmeta< 4.10×10-6 after Bonferroni correction for the total number of genes tested. The likely pathogenicity of these variants was determined in silico using SIFT and PolyPhen. In recipients, we identified OR51D1, a member of the olfactory receptor gene family, as significantly associated with both OS and TRM. The OS association with ORD51D1 is driven by TRM as the same six missense variants contribute to both TRM and OS, but not to DRM. In donors, four (ALPP, EMID1, SLC44A5, LRP1), one (HHAT), and two (LYZL4, NT5E) genes were significantly associated with OS, TRM, and DRM, respectively (Table). ALPP encodes an alkaline phosphatase. SLC44A5 is a member of the solute carrier gene family. LRP1 (low density lipoprotein receptor-related protein 1) encodes an endocytic receptor involved in several cellular processes, including intracellular signaling, lipid homeostasis, and clearance of apoptotic cells. HHAT encodes an enzyme that acts within the hedgehog secretory pathway, which is involved in hematologic malignancy. LYZL4 is a member of a family of lysozyme-like genes, which have a protective role in host defense. The encoded protein of NT5E, CD270, is used as a determinant of lymphocyte differentiation. The likely pathogenicity of several of the functional variants in all significant genes identified was predicted to be both (possibly or probably) damaging and deleterious (Table), thus indicating that the amino acid substitution of the variant may not be well tolerated. One damaging variant in LRP1 has been identified as a somatic mutation in hematopoietic and lymphoid tissue. Our study is the first analysis of the low-frequency coding variant contribution to URD BMT survival outcomes. We found that the missense variants contributing to many of these significant gene associations shows evidence of pathogenicity and thus it is biologically plausible these variants are contributing to survival outcomes. Further confirmation of these findings, and studies of the functional consequences of protein-coding changes in these genes, may provide more individualized risk prediction and prognosis as well as alternative donor selection strategies. Table. Table. Disclosures McCarthy: Gamida Cell: Honoraria, Membership on an entity's Board of Directors or advisory committees; Onyx: Honoraria, Membership on an entity's Board of Directors or advisory committees; Janssen: Honoraria, Membership on an entity's Board of Directors or advisory committees; Celgene: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; Bristol Myers Squibb: Honoraria, Membership on an entity's Board of Directors or advisory committees; Karyopharm: Honoraria, Membership on an entity's Board of Directors or advisory committees; The Binding Site: Honoraria, Membership on an entity's Board of Directors or advisory committees; Sanofi: Honoraria, Membership on an entity's Board of Directors or advisory committees. Hahn:Novartis: Equity Ownership; NIH: Research Funding. Sucheston-Campbell:NIH/NCI: Research Funding.

Abstract Over the past decade, there have been numerous published candidate gene association studies of non-HLA variants, genes and/or pathways in relation to BMT survival outcomes. Most of these studies investigated a small number of candidate single nucleotide polymorphisms (SNPs), chosen based on known gene function, in a limited number of cases without replication. Here we test the association of previously reported SNPs with survival outcomes (overall survival (OS), progression-free survival (PFS), death due to disease (DD) or transplant-related mortality (TRM)) after URD BMT using data from the DISCOVeRY-BMT (Determining the Influence of Susceptibility Conveying Variants Related to one-Year mortality after BMT) study which included 2,052 European descent recipient-unrelated donor pairs (Cohort 1) and 763 recipient-unrelated donor pairs (Cohort 1) reported to the Center for International Blood and Marrow Transplant Research (CIBMTR) between 2000-2009 and 2009-2011, respectively. DISCOVeRY-BMT is a Genome Wide Association Study (GWAS) that includes over 9 million typed or imputed SNPs. This represents the largest study to date and is well powered to test both SNP and gene-based association hypotheses. We conducted an extensive literature review of PubMed using MESH terms related to transplant and hematologic disease which yielded 76 candidate gene association studies that investigated a total of 46 genes with any survival outcome. Of the 46 genes, 11 genes (CTLA4, CYP2B6, CYP2C19, GSTM1, GSTT1, IL23R, IL6, MTHFR, NOD2/CARD15, TNFSRF1B, and VDR) were studied in at least 2 publications. A total of 16 SNPs in 8 of these genes reported a significant (P< 0.05) association with at least one survival outcome (OS, PFS, DD or TRM). We tested these 16 SNPs for an association with the survival outcome originally reported in each article (Table 1). All association models included age at BMT, diagnosis (AML, ALL, MDS), disease status at BMT, cell source (blood, bone marrow) and year of BMT. P-values from Cohorts 1 and 2 were combined using a fixed-effects model in METAL software. Since the candidate gene studies have a gene-based hypothesis, we also performed gene-level testing using the VErsatile Gene-Based Association Study (VEGAS2) software to measure the aggregate effect of all available typed and imputed SNPs from our GWAS data in these 8 genes. The candidate gene association studies included heterogeneous patient populations of autologous and/or allogeneic (related or unrelated donors) BMT for a variety of hematologic diseases. Most candidate gene studies performed genotyping on both donors and recipients, and SNPs were tested either independently for donor or recipient genomes, or considering the total number of variant alleles the recipient-donor pair had at a given locus (0-4) (Table 1). The 16 SNPs from the previous candidate gene association studies were not significantly associated at a nominal P<0.05 to either the survival outcome reported in the original published paper, or to any survival outcome (Table 1). Likewise, the VEGAS2 gene-based analysis did not reveal any statistically significant associations for any of the 8 genes at P<0.05. Our study also corroborated prior reports of SNPs and genes that were not associated with survival outcomes. Candidate gene association studies use a priori biological knowledge to select genes that are likely associated with outcome and subsequently SNPs within the gene thought likely to affect function based on knowledge at the time of the investigation. Our findings suggest that neither the genetic variation tested in these candidate genes, nor other individual SNPs in the gene independently or collectively are not associated with survival outcomes in recipients of URD BMT. Our failure to replicate the prior candidate gene studies may be due to these SNPs being dependent on the underlying disease (ie. AML, ALL, etc), the donor relation (related, unrelated, autologous), or that there is no true association of these SNPs with the outcomes studied. #These authors contributed equally Table 1. Table 1. Disclosures McCarthy: Gamida Cell: Honoraria, Membership on an entity's Board of Directors or advisory committees; Celgene: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; Karyopharm: Honoraria, Membership on an entity's Board of Directors or advisory committees; Onyx: Honoraria, Membership on an entity's Board of Directors or advisory committees; Bristol Myers Squibb: Honoraria, Membership on an entity's Board of Directors or advisory committees; Janssen: Honoraria, Membership on an entity's Board of Directors or advisory committees; The Binding Site: Honoraria, Membership on an entity's Board of Directors or advisory committees; Sanofi: Honoraria, Membership on an entity's Board of Directors or advisory committees. Hahn:Novartis: Equity Ownership; NIH: Research Funding. Sucheston-Campbell:NIH/NCI: Research Funding.

Key PointsIKZF1 associations with high-risk B-ALL may differ by age and sex. A novel variant on chromosome 14, rs189434316, is associated with over a 3.5-fold risk of normal cytogenetic B-ALL.

Key Points We conducted the first exome-wide association study between germ line variant genotype and survival outcomes after unrelated-donor BMT. A number of novel genes were found to significantly affect survival outcomes.

Key PointsCandidate SNP associations with survival outcomes after URD transplant are most likely false-positive findings. Over 85% of candidate SNPs are not linked to a biochemical function; of those that are, about half are not linked to the candidate gene.