Association between telomere length and Plasmodium falciparum malaria endemicity in sub-Saharan Africans

Leukocyte telomere length (LTL) varies significantly across human populations, with individuals of African ancestry having longer LTL than non-Africans. However, the genetic and environmental drivers of LTL variation in Africans remain largely unknown. We report here on the relationship between LTL, genetics, and a variety of environmental and climatic factors in ethnically diverse African adults (n = 1,818) originating from Botswana, Tanzania, Ethiopia, and Cameroon. We observe significant variation in LTL among populations, finding that the San hunter-gatherers from Botswana have the longest leukocyte telomeres and that the Fulani pastoralists from Cameroon have the shortest telomeres. Genetic factors explain ∼50% of LTL variation among individuals. Moreover, we observe a significant negative association between Plasmodium falciparum malaria endemicity and LTL while adjusting for age, sex, and genetics. Within Africa, adults from populations indigenous to areas with high malaria exposure have shorter LTL than those in populations indigenous to areas with low malaria exposure. Finally, we explore to what degree the genetic architecture underlying LTL in Africa covaries with malaria exposure.

Geographic variation of mutagenic exposures in kidney cancer genomes

International differences in the incidence of many cancer types indicate the existence of carcinogen exposures that have not yet been identified by conventional epidemiology make a substantial contribution to cancer burden1. In clear cell renal cell carcinoma, obesity, hypertension and tobacco smoking are risk factors, but they do not explain the geographical variation in its incidence2. Underlying causes can be inferred by sequencing the genomes of cancers from populations with different incidence rates and detecting differences in patterns of somatic mutations. Here we sequenced 962 clear cell renal cell carcinomas from 11 countries with varying incidence. The somatic mutation profiles differed between countries. In Romania, Serbia and Thailand, mutational signatures characteristic of aristolochic acid compounds were present in most cases and but these were rare elsewhere. In Japan, a mutational signature of unknown cause was found in more than 70% of cases but in less than 2% elsewhere. A further mutational signature of unknown cause was ubiquitous but exhibited higher mutation loads in countries with higher incidence rates of kidney cancer. Known signatures of tobacco smoking correlated with tobacco consumption, but no signature was associated with obesity or hypertension, suggesting that non-mutagenic mechanisms of action underlie these risk factors. The results of this study indicate the existence of multiple, geographically variable, mutagenic exposures that potentially affect tens of millions of people and illustrate the opportunities for new insights into cancer causation through large-scale global cancer genomics.

Multi-ancestry genome-wide association study of kidney cancer identifies 63 susceptibility regions

Here, in a multi-ancestry genome-wide association study meta-analysis of kidney cancer (29,020 cases and 835,670 controls), we identified 63 susceptibility regions (50 novel) containing 108 independent risk loci. In analyses stratified by subtype, 52 regions (78 loci) were associated with clear cell renal cell carcinoma (RCC) and 6 regions (7 loci) with papillary RCC. Notably, we report a variant common in African ancestry individuals ( rs7629500 ) in the 3' untranslated region of VHL, nearly tripling clear cell RCC risk (odds ratio 2.72, 95% confidence interval 2.23-3.30). In cis-expression quantitative trait locus analyses, 48 variants from 34 regions point toward 83 candidate genes. Enrichment of hypoxia-inducible factor-binding sites underscores the importance of hypoxia-related mechanisms in kidney cancer. Our results advance understanding of the genetic architecture of kidney cancer, provide clues for functional investigation and enable generation of a validated polygenic risk score with an estimated area under the curve of 0.65 (0.74 including risk factors) among European ancestry individuals.

Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

Novel breast cancer susceptibility loci under linkage peaks identified in African ancestry consortia

BACKGROUND

Expansion of genome-wide association studies across population groups is needed to improve our understanding of shared and unique genetic contributions to breast cancer. We performed association and replication studies guided by a priori linkage findings from African ancestry (AA) relative pairs.

METHODS

We performed fixed-effect inverse-variance weighted meta-analysis under three significant AA breast cancer linkage peaks (3q26-27, 12q22-23, and 16q21-22) in 9241 AA cases and 10 193 AA controls. We examined associations with overall breast cancer as well as estrogen receptor (ER)-positive and negative subtypes (193,132 SNPs). We replicated associations in the African-ancestry Breast Cancer Genetic Consortium (AABCG).

RESULTS

In AA women, we identified two associations on chr12q for overall breast cancer (rs1420647, OR = 1.15, p = 2.50×10-6; rs12322371, OR = 1.14, p = 3.15×10-6), and one for ER-negative breast cancer (rs77006600, OR = 1.67, p = 3.51×10-6). On chr3, we identified two associations with ER-negative disease (rs184090918, OR = 3.70, p = 1.23×10-5; rs76959804, OR = 3.57, p = 1.77×10-5) and on chr16q we identified an association with ER-negative disease (rs34147411, OR = 1.62, p = 8.82×10-6). In the replication study, the chr3 associations were significant and effect sizes were larger (rs184090918, OR: 6.66, 95% CI: 1.43, 31.01; rs76959804, OR: 5.24, 95% CI: 1.70, 16.16).

CONCLUSION

The two chr3 SNPs are upstream to open chromatin ENSR00000710716, a regulatory feature that is actively regulated in mammary tissues, providing evidence that variants in this chr3 region may have a regulatory role in our target organ. Our study provides support for breast cancer variant discovery using prioritization based on linkage evidence.

APOBEC shapes tumor evolution and age at onset of lung cancer in smokers

APOBEC enzymes are part of the innate immunity and are responsible for restricting viruses and retroelements by deaminating cytosine residues1,2. Most solid tumors harbor different levels of somatic mutations attributed to the off-target activities of APOBEC3A (A3A) and/or APOBEC3B (A3B)3-6. However, how APOBEC3A/B enzymes shape the tumor evolution in the presence of exogenous mutagenic processes is largely unknown. Here, by combining deep whole-genome sequencing with multi-omics profiling of 309 lung cancers from smokers with detailed tobacco smoking information, we identify two subtypes defined by low (LAS) and high (HAS) APOBEC mutagenesis. LAS are enriched for A3B-like mutagenesis and KRAS mutations, whereas HAS for A3A-like mutagenesis and TP53 mutations. Unlike APOBEC3A, APOBEC3B expression is strongly associated with an upregulation of the base excision repair pathway. Hypermutation by unrepaired A3A and tobacco smoking mutagenesis combined with TP53-induced genomic instability can trigger senescence7, apoptosis8, and cell regeneration9, as indicated by high expression of pulmonary healing signaling pathway, stemness markers and distal cell-of-origin in HAS. The expected association of tobacco smoking variables (e.g., time to first cigarette) with genomic/epigenomic changes are not observed in HAS, a plausible consequence of frequent cell senescence or apoptosis. HAS have more neoantigens, slower clonal expansion, and older age at onset compared to LAS, particularly in heavy smokers, consistent with high proportions of newly generated, unmutated cells and frequent immuno-editing. These findings show how heterogeneity in mutational burden across co-occurring mutational processes and cell types contributes to tumor development, with important clinical implications.

Pan-cancer mutational signature analysis of 111,711 targeted sequenced tumors using SATS

Tumor mutational signatures have the potential to inform cancer diagnosis and treatment. However, their detection in targeted sequenced tumors is hampered by sparse mutations and variability in targeted gene panels. Here we present SATS, a scalable mutational signature analyzer addressing these challenges by leveraging tumor mutational burdens from targeted gene panels. Through analyzing simulated data, pseudo-targeted sequencing data generated by down-sampling whole exome and genome data, and samples with matched whole genome sequencing and targeted sequencing, we showed that SATS can accurately detect common mutational signatures and estimate signature burdens. Applying SATS to 111,711 targeted sequenced tumors from the AACR Project GENIE, we generated a pan-cancer catalogue of mutational signatures tailored to targeted sequencing, enabling estimation of signature burdens within individual tumors. Integrating signatures with clinical data, we demonstrated SATS's clinical utility, including identifying signatures enriched in early-onset hypermutated colorectal cancers and signatures associated with cancer prognosis and immunotherapy response.

Polygenic risk scores, radiation treatment exposures and subsequent cancer risk in childhood cancer survivors

Survivors of childhood cancer are at increased risk for subsequent cancers attributable to the late effects of radiotherapy and other treatment exposures; thus, further understanding of the impact of genetic predisposition on risk is needed. Combining genotype data for 11,220 5-year survivors from the Childhood Cancer Survivor Study and the St Jude Lifetime Cohort, we found that cancer-specific polygenic risk scores (PRSs) derived from general population, genome-wide association study, cancer loci identified survivors of European ancestry at increased risk of subsequent basal cell carcinoma (odds ratio per s.d. of the PRS: OR = 1.37, 95% confidence interval (CI) = 1.29-1.46), female breast cancer (OR = 1.42, 95% CI = 1.27-1.58), thyroid cancer (OR = 1.48, 95% CI = 1.31-1.67), squamous cell carcinoma (OR = 1.20, 95% CI = 1.00-1.44) and melanoma (OR = 1.60, 95% CI = 1.31-1.96); however, the association for colorectal cancer was not significant (OR = 1.19, 95% CI = 0.94-1.52). An investigation of joint associations between PRSs and radiotherapy found more than additive increased risks of basal cell carcinoma, and breast and thyroid cancers. For survivors with radiotherapy exposure, the cumulative incidence of subsequent cancer by age 50 years was increased for those with high versus low PRS. These findings suggest a degree of shared genetic etiology for these malignancy types in the general population and survivors, which remains evident in the context of strong radiotherapy-related risk.

State of the science and future directions for research on HIV and cancer: Summary of a joint workshop sponsored by IARC and NCI

An estimated 38 million people live with human immunodeficiency virus (HIV) worldwide and are at excess risk for multiple cancer types. Elevated cancer risks in people living with HIV (PLWH) are driven primarily by increased exposure to carcinogens, most notably oncogenic viruses acquired through shared transmission routes, plus acceleration of viral carcinogenesis by HIV-related immunosuppression. In the era of widespread antiretroviral therapy (ART), life expectancy of PLWH has increased, with cancer now a leading cause of co-morbidity and death. Furthermore, the types of cancers occurring among PLWH are shifting over time and vary in their relative burden in different parts of the world. In this context, the International Agency for Research on Cancer (IARC) and the US National Cancer Institute (NCI) convened a meeting in September 2022 of multinational and multidisciplinary experts to focus on cancer in PLWH. This report summarizes the proceedings, including a review of the state of the science of cancer descriptive epidemiology, etiology, molecular tumor characterization, primary and secondary prevention, treatment disparities and survival in PLWH around the world. A consensus of key research priorities and recommendations in these domains is also presented.

Association between circulating inflammatory markers and adult cancer risk: a Mendelian randomization analysis

BACKGROUND

Tumour-promoting inflammation is a "hallmark" of cancer and conventional epidemiological studies have reported links between various inflammatory markers and cancer risk. The causal nature of these relationships and, thus, the suitability of these markers as intervention targets for cancer prevention is unclear.

METHODS

We meta-analysed 6 genome-wide association studies of circulating inflammatory markers comprising 59,969 participants of European ancestry. We then used combined cis-Mendelian randomization and colocalisation analysis to evaluate the causal role of 66 circulating inflammatory markers in risk of 30 adult cancers in 338,294 cancer cases and up to 1,238,345 controls. Genetic instruments for inflammatory markers were constructed using genome-wide significant (P < 5.0 × 10-8) cis-acting SNPs (i.e., in or ±250 kb from the gene encoding the relevant protein) in weak linkage disequilibrium (LD, r2 < 0.10). Effect estimates were generated using inverse-variance weighted random-effects models and standard errors were inflated to account for weak LD between variants with reference to the 1000 Genomes Phase 3 CEU panel. A false discovery rate (FDR)-corrected P-value ("q-value") <0.05 was used as a threshold to define "strong evidence" to support associations and 0.05 ≤ q-value < 0.20 to define "suggestive evidence". A colocalisation posterior probability (PPH4) >70% was employed to indicate support for shared causal variants across inflammatory markers and cancer outcomes. Findings were replicated in the FinnGen study and then pooled using meta-analysis.

FINDINGS

We found strong evidence to support an association of genetically-proxied circulating pro-adrenomedullin concentrations with increased breast cancer risk (OR: 1.19, 95% CI: 1.10-1.29, q-value = 0.033, PPH4 = 84.3%) and suggestive evidence to support associations of interleukin-23 receptor concentrations with increased pancreatic cancer risk (OR: 1.42, 95% CI: 1.20-1.69, q-value = 0.055, PPH4 = 73.9%), prothrombin concentrations with decreased basal cell carcinoma risk (OR: 0.66, 95% CI: 0.53-0.81, q-value = 0.067, PPH4 = 81.8%), and interleukin-1 receptor-like 1 concentrations with decreased triple-negative breast cancer risk (OR: 0.92, 95% CI: 0.88-0.97, q-value = 0.15, PPH4 = 85.6%). These findings were replicated in pooled analyses with the FinnGen study. Though suggestive evidence was found to support an association of macrophage migration inhibitory factor concentrations with increased bladder cancer risk (OR: 2.46, 95% CI: 1.48-4.10, q-value = 0.072, PPH4 = 76.1%), this finding was not replicated when pooled with the FinnGen study. For 22 of 30 cancer outcomes examined, there was little evidence (q-value ≥0.20) that any of the 66 circulating inflammatory markers examined were associated with cancer risk.

INTERPRETATION

Our comprehensive joint Mendelian randomization and colocalisation analysis of the role of circulating inflammatory markers in cancer risk identified potential roles for 4 circulating inflammatory markers in risk of 4 site-specific cancers. Contrary to reports from some prior conventional epidemiological studies, we found little evidence of association of circulating inflammatory markers with the majority of site-specific cancers evaluated.

FUNDING

Cancer Research UK (C68933/A28534, C18281/A29019, PPRCPJT∖100005), World Cancer Research Fund (IIG_FULL_2020_022), National Institute for Health Research (NIHR202411, BRC-1215-20011), Medical Research Council (MC_UU_00011/1, MC_UU_00011/3, MC_UU_00011/6, and MC_UU_00011/4), Academy of Finland Project 326291, European Union's Horizon 2020 grant agreement no. 848158 (EarlyCause), French National Cancer Institute (INCa SHSESP20, 2020-076), Versus Arthritis (21173, 21754, 21755), National Institutes of Health (U19 CA203654), National Cancer Institute (U19CA203654).

Human leukocyte antigen-DQA1*04:01 and rs2040406 variants are associated with elevated risk of childhood Burkitt lymphoma

Burkitt lymphoma (BL) is responsible for many childhood cancers in sub-Saharan Africa, where it is linked to recurrent or chronic infection by Epstein-Barr virus or Plasmodium falciparum. However, whether human leukocyte antigen (HLA) polymorphisms, which regulate immune response, are associated with BL has not been well investigated, which limits our understanding of BL etiology. Here we investigate this association among 4,645 children aged 0-15 years, 800 with BL, enrolled in Uganda, Tanzania, Kenya, and Malawi. HLA alleles are imputed with accuracy >90% for HLA class I and 85-89% for class II alleles. BL risk is elevated with HLA-DQA1*04:01 (adjusted odds ratio [OR] = 1.61, 95% confidence interval [CI] = 1.32-1.97, P = 3.71 × 10-6), with rs2040406(G) in HLA-DQA1 region (OR = 1.43, 95% CI = 1.26-1.63, P = 4.62 × 10-8), and with amino acid Gln at position 53 versus other variants in HLA-DQA1 (OR = 1.36, P = 2.06 × 10-6). The associations with HLA-DQA1*04:01 (OR = 1.29, P = 0.03) and rs2040406(G) (OR = 1.68, P = 0.019) persist in mutually adjusted models. The higher risk rs2040406(G) variant for BL is associated with decreased HLA-DQB1 expression in eQTLs in EBV transformed lymphocytes. Our results support the role of HLA variation in the etiology of BL and suggest that a promising area of research might be understanding the link between HLA variation and EBV control.

Sickle cell allele HBB-rs334(T) is associated with decreased risk of childhood Burkitt lymphoma in East Africa

Burkitt lymphoma (BL) is an aggressive B-cell lymphoma that significantly contributes to childhood cancer burden in sub-Saharan Africa. Plasmodium falciparum, which causes malaria, is geographically associated with BL, but the evidence remains insufficient for causal inference. Inference could be strengthened by demonstrating that mendelian genes known to protect against malaria-such as the sickle cell trait variant, HBB-rs334(T)-also protect against BL. We investigated this hypothesis among 800 BL cases and 3845 controls in four East African countries using genome-scan data to detect polymorphisms in 22 genes known to affect malaria risk. We fit generalized linear mixed models to estimate odds ratios (OR) and 95% confidence intervals (95% CI), controlling for age, sex, country, and ancestry. The ORs of the loci with BL and P. falciparum infection among controls were correlated (Spearman's ρ = 0.37, p = .039). HBB-rs334(T) was associated with lower P. falciparum infection risk among controls (OR = 0.752, 95% CI 0.628-0.9; p = .00189) and BL risk (OR = 0.687, 95% CI 0.533-0.885; p = .0037). ABO-rs8176703(T) was associated with decreased risk of BL (OR = 0.591, 95% CI 0.379-0.992; p = .00271), but not of P. falciparum infection. Our results increase support for the etiological correlation between P. falciparum and BL risk.

Observational and genetic associations between cardiorespiratory fitness and cancer: a UK Biobank and international consortia study

BACKGROUND

The association of fitness with cancer risk is not clear.

METHODS

We used Cox proportional hazards models to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for risk of lung, colorectal, endometrial, breast, and prostate cancer in a subset of UK Biobank participants who completed a submaximal fitness test in 2009-12 (N = 72,572). We also investigated relationships using two-sample Mendelian randomisation (MR), odds ratios (ORs) were estimated using the inverse-variance weighted method.

RESULTS

After a median of 11 years of follow-up, 4290 cancers of interest were diagnosed. A 3.5 ml O2⋅min-1⋅kg-1 total-body mass increase in fitness (equivalent to 1 metabolic equivalent of task (MET), approximately 0.5 standard deviation (SD)) was associated with lower risks of endometrial (HR = 0.81, 95% CI: 0.73-0.89), colorectal (0.94, 0.90-0.99), and breast cancer (0.96, 0.92-0.99). In MR analyses, a 0.5 SD increase in genetically predicted O2⋅min-1⋅kg-1 fat-free mass was associated with a lower risk of breast cancer (OR = 0.92, 95% CI: 0.86-0.98). After adjusting for adiposity, both the observational and genetic associations were attenuated.

DISCUSSION

Higher fitness levels may reduce risks of endometrial, colorectal, and breast cancer, though relationships with adiposity are complex and may mediate these relationships. Increasing fitness, including via changes in body composition, may be an effective strategy for cancer prevention.

JAK2 V617F mutation and associated chromosomal alterations in primary and secondary myelofibrosis and post-HCT outcomes

JAK2 V617F is the most common driver mutation in primary or secondary myelofibrosis for which allogeneic hematopoietic cell transplantation (HCT) is the only curative treatment. Knowledge of the prognostic utility of JAK2 alterations in the HCT setting is limited. We identified all patients with MF who received HCT between 2000 and 2016 and had a pre-HCT blood sample (N = 973) available at the Center of International Blood and Marrow Transplant Research biorepository. PacBio sequencing and single nucleotide polymorphism-array genotyping were used to identify JAK2V617F mutation and associated mosaic chromosomal alterations (mCAs), respectively. Cox proportional hazard models were used for HCT outcome analyses. Genomic testing was complete for 924 patients with MF (634 primary MF [PMF], 135 postpolycythemia vera [PPV-MF], and 155 postessential thrombocytopenia [PET-MF]). JAK2V617F affected 562 patients (57.6% of PMF, 97% of PPV-MF, and 42.6% of PET-MF). Almost all patients with mCAs involving the JAK2 region (97.9%) were JAK2V617-positive. In PMF, JAK2V617F mutation status, allele burden, or identified mCAs were not associated with disease progression/relapse, nonrelapse mortality (NRM), or overall survival. Almost all PPV-MF were JAK2V617F-positive (97%), with no association between HCT outcomes and mutation allele burden or mCAs. In PET-MF, JAK2V617F high mutation allele burden (≥60%) was associated with excess risk of NRM, restricted to transplants received in the era of JAK inhibitors (2013-2016; hazard ratio = 7.65; 95% confidence interval = 2.10-27.82; P = .002). However, allele burden was not associated with post-HCT disease progression/relapse or survival. Our findings support the concept that HCT can mitigate the known negative effect of JAK2V617F in patients with MF, particularly for PMF and PPV-MF.

Mosaic chromosomal alterations in peripheral blood leukocytes of children in sub-Saharan Africa

In high-income countries, mosaic chromosomal alterations in peripheral blood leukocytes are associated with an elevated risk of adverse health outcomes, including hematologic malignancies. We investigate mosaic chromosomal alterations in sub-Saharan Africa among 931 children with Burkitt lymphoma, an aggressive lymphoma commonly characterized by immunoglobulin-MYC chromosomal rearrangements, 3822 Burkitt lymphoma-free children, and 674 cancer-free men from Ghana. We find autosomal and X chromosome mosaic chromosomal alterations in 3.4% and 1.7% of Burkitt lymphoma-free children, and 8.4% and 3.7% of children with Burkitt lymphoma (P-values = 5.7×10-11 and 3.74×10-2, respectively). Autosomal mosaic chromosomal alterations are detected in 14.0% of Ghanaian men and increase with age. Mosaic chromosomal alterations in Burkitt lymphoma cases include gains on chromosomes 1q and 8, the latter spanning MYC, while mosaic chromosomal alterations in Burkitt lymphoma-free children include copy-neutral loss of heterozygosity on chromosomes 10, 14, and 16. Our results highlight mosaic chromosomal alterations in sub-Saharan African populations as a promising area of research.

Characterizing prostate cancer risk through multi-ancestry genome-wide discovery of 187 novel risk variants

The transferability and clinical value of genetic risk scores (GRSs) across populations remain limited due to an imbalance in genetic studies across ancestrally diverse populations. Here we conducted a multi-ancestry genome-wide association study of 156,319 prostate cancer cases and 788,443 controls of European, African, Asian and Hispanic men, reflecting a 57% increase in the number of non-European cases over previous prostate cancer genome-wide association studies. We identified 187 novel risk variants for prostate cancer, increasing the total number of risk variants to 451. An externally replicated multi-ancestry GRS was associated with risk that ranged from 1.8 (per standard deviation) in African ancestry men to 2.2 in European ancestry men. The GRS was associated with a greater risk of aggressive versus non-aggressive disease in men of African ancestry (P = 0.03). Our study presents novel prostate cancer susceptibility loci and a GRS with effective risk stratification across ancestry groups.

Health diplomacy in action: The cancer legacy of the Good Friday Agreement

2023 marks the 25th anniversary of the Good Friday Agreement, which led peace in Northern Ireland. As well as its impact on peace and reconciliation, the Good Friday Agreement has also had a lasting positive impact on cancer research and cancer care across the island of Ireland. Pursuant to the Good Friday Agreement, a Memorandum of Understanding (MOU) was signed between the respective Departments of Health in Ireland, Northern Ireland and the US National Cancer Institute (NCI), giving rise to the Ireland - Northern Ireland - National Cancer Institute Cancer Consortium, an unparalleled tripartite agreement designed to nurture and develop linkages between cancer researchers, physicians and allied healthcare professionals across Ireland, Northern Ireland and the US, delivering world class research and better care for cancer patients on the island of Ireland and driving research and innovation in the US.

Serum biomarkers are altered in UK Biobank participants with mosaic chromosomal alterations

Age-related clonal expansion of cells harbouring mosaic chromosomal alterations (mCAs) is one manifestation of clonal haematopoiesis. Identifying factors that influence the generation and promotion of clonal expansion of mCAs are key to investigate the role of mCAs in health and disease. Herein, we report on widely measured serum biomarkers and their possible association with mCAs, which could provide new insights into molecular alterations that promote acquisition and clonal expansion. We performed a cross-sectional investigation of the association of 32 widely measured serum biomarkers with autosomal mCAs, mosaic loss of the Y chromosome, and mosaic loss of the X chromosome in 436 784 cancer-free participants from the UK Biobank. mCAs were associated with a range of commonly measured serum biomarkers such as lipid levels, circulating sex hormones, blood sugar homeostasis, inflammation and immune function, vitamins and minerals, kidney function, and liver function. Biomarker levels in participants with mCAs were estimated to differ by up to 5% relative to mCA-free participants, and individuals with higher cell fraction mCAs had greater deviation in mean biomarker values. Polygenic scores associated with sex hormone binding globulin, vitamin D, and total cholesterol were also associated with mCAs. Overall, we observed commonly used clinical serum biomarkers related to disease risk are associated with mCAs, suggesting mechanisms involved in these diseases could be related to mCA proliferation and clonal expansion.

Genome-wide association study of prostate-specific antigen levels in 392,522 men identifies new loci and improves cross-ancestry prediction

We conducted a multi-ancestry genome-wide association study of prostate-specific antigen (PSA) levels in 296,754 men (211,342 European ancestry; 58,236 African ancestry; 23,546 Hispanic/Latino; 3,630 Asian ancestry; 96.5% of participants were from the Million Veteran Program). We identified 318 independent genome-wide significant (p≤5e-8) variants, 184 of which were novel. Most demonstrated evidence of replication in an independent cohort (n=95,768). Meta-analyzing discovery and replication (n=392,522) identified 447 variants, of which a further 111 were novel. Out-of-sample variance in PSA explained by our new polygenic risk score reached 16.9% (95% CI=16.1%-17.8%) in European ancestry, 9.5% (95% CI=7.0%-12.2%) in African ancestry, 18.6% (95% CI=15.8%-21.4%) in Hispanic/Latino, and 15.3% (95% CI=12.7%-18.1%) in Asian ancestry, and lower for higher age. Our study highlights how including proportionally more participants from underrepresented populations improves genetic prediction of PSA levels, with potential to personalize prostate cancer screening.

Heterogeneous genetic architectures and evolutionary genomics of prostate cancer in Sub-Saharan Africa

Men of African descent have the highest prostate cancer (CaP) incidence and mortality rates, yet the genetic basis of CaP in African men has been understudied. We used genomic data from 3,963 CaP cases and 3,509 controls recruited in Ghana, Nigeria, Senegal, South Africa, and Uganda, to infer ancestry-specific genetic architectures and fine-mapped disease associations. Fifteen independent associations at 8q24.21, 6q22.1, and 11q13.3 reached genome-wide significance, including four novel associations. Intriguingly, multiple lead SNPs are private alleles, a pattern arising from recent mutations and the out-of-Africa bottleneck. These African-specific alleles contribute to haplotypes with odds ratios above 2.4. We found that the genetic architecture of CaP differs across Africa, with effect size differences contributing more to this heterogeneity than allele frequency differences. Population genetic analyses reveal that African CaP associations are largely governed by neutral evolution. Collectively, our findings emphasize the utility of conducting genetic studies that use diverse populations.

PRScalc, a privacy-preserving calculation of raw polygenic risk scores from direct-to-consumer genomics data

Motivation

Currently, the Polygenic Score (PGS) Catalog curates over 400 publications on over 500 traits corresponding to over 3000 polygenic risk scores (PRSs). To assess the feasibility of privately calculating the underlying multivariate relative risk for individuals with consumer genomics data, we developed an in-browserPRS calculator for genomic data that does not circulate any data or engage in any computation outside of the user's personal device.

Results

A prototype personal risk score calculator, created for research purposes, was developed to demonstrate how the PGS Catalog can be privately and readily applied to readily available direct-to-consumer genetic testing services, such as 23andMe. No software download, installation, or configuration is needed. The PRS web calculator matches individual PGS catalog entries with an individual's 23andMe genome data composed of 600k to 1.4 M single-nucleotide polymorphisms (SNPs). Beta coefficients provide researchers with a convenient assessment of risk associated with matched SNPs. This in-browser application was tested in a variety of personal devices, including smartphones, establishing the feasibility of privately calculating personal risk scores with up to a few thousand reference genetic variations and from the full 23andMe SNP data file (compressed or not).

Availability and implementation

The PRScalc web application is developed in JavaScript, HTML, and CSS and is available at GitHub repository (https://episphere.github.io/prs) under an MIT license. The datasets were derived from sources in the public domain: [PGS Catalog, Personal Genome Project].

Combining Asian and European genome-wide association studies of colorectal cancer improves risk prediction across racial and ethnic populations

Polygenic risk scores (PRS) have great potential to guide precision colorectal cancer (CRC) prevention by identifying those at higher risk to undertake targeted screening. However, current PRS using European ancestry data have sub-optimal performance in non-European ancestry populations, limiting their utility among these populations. Towards addressing this deficiency, we expand PRS development for CRC by incorporating Asian ancestry data (21,731 cases; 47,444 controls) into European ancestry training datasets (78,473 cases; 107,143 controls). The AUC estimates (95% CI) of PRS are 0.63(0.62-0.64), 0.59(0.57-0.61), 0.62(0.60-0.63), and 0.65(0.63-0.66) in independent datasets including 1681-3651 cases and 8696-115,105 controls of Asian, Black/African American, Latinx/Hispanic, and non-Hispanic White, respectively. They are significantly better than the European-centric PRS in all four major US racial and ethnic groups (p-values < 0.05). Further inclusion of non-European ancestry populations, especially Black/African American and Latinx/Hispanic, is needed to improve the risk prediction and enhance equity in applying PRS in clinical practice.

Shared and distinct genetic etiologies for different types of clonal hematopoiesis

Clonal hematopoiesis (CH)-age-related expansion of mutated hematopoietic clones-can differ in frequency and cellular fitness by CH type (e.g., mutations in driver genes (CHIP), gains/losses and copy-neutral loss of chromosomal segments (mCAs), and loss of sex chromosomes). Co-occurring CH raises questions as to their origin, selection, and impact. We integrate sequence and genotype array data in up to 482,378 UK Biobank participants to demonstrate shared genetic architecture across CH types. Our analysis suggests a cellular evolutionary trade-off between different types of CH, with LOY occurring at lower rates in individuals carrying mutations in established CHIP genes. We observed co-occurrence of CHIP and mCAs with overlap at TET2, DNMT3A, and JAK2, in which CHIP precedes mCA acquisition. Furthermore, individuals carrying overlapping CH had high risk of future lymphoid and myeloid malignancies. Finally, we leverage shared genetic architecture of CH traits to identify 15 novel loci associated with leukemia risk.

Genetic Susceptibility to Nonalcoholic Fatty Liver Disease and Risk for Pancreatic Cancer: Mendelian Randomization

BACKGROUND

There are conflicting data on whether nonalcoholic fatty liver disease (NAFLD) is associated with susceptibility to pancreatic cancer. Using Mendelian randomization (MR), we investigated the relationship between genetic predisposition to NAFLD and risk for pancreatic cancer.

METHODS

Data from genome-wide association studies (GWAS) within the Pancreatic Cancer Cohort Consortium (PanScan; cases n = 5,090, controls n = 8,733) and the Pancreatic Cancer Case Control Consortium (PanC4; cases n = 4,163, controls n = 3,792) were analyzed. We used data on 68 genetic variants with four different MR methods [inverse variance weighting (IVW), MR-Egger, simple median, and penalized weighted median] separately to predict genetic heritability of NAFLD. We then assessed the relationship between each of the four MR methods and pancreatic cancer risk, using logistic regression to calculate ORs and 95% confidence intervals (CI), adjusting for PC risk factors, including obesity and diabetes.

RESULTS

No association was found between genetically predicted NAFLD and pancreatic cancer risk in the PanScan or PanC4 samples [e.g., PanScan, IVW OR, 1.04; 95% confidence interval (CI), 0.88-1.22; MR-Egger OR, 0.89; 95% CI, 0.65-1.21; PanC4, IVW OR, 1.07; 95% CI, 0.90-1.27; MR-Egger OR, 0.93; 95% CI, 0.67-1.28]. None of the four MR methods indicated an association between genetically predicted NAFLD and pancreatic cancer risk in either sample.

CONCLUSIONS

Genetic predisposition to NAFLD is not associated with pancreatic cancer risk.

IMPACT

Given the close relationship between NAFLD and metabolic conditions, it is plausible that any association between NAFLD and pancreatic cancer might reflect host metabolic perturbations (e.g., obesity, diabetes, or metabolic syndrome) and does not necessarily reflect a causal relationship between NAFLD and pancreatic cancer.

The Childhood Cancer Data Initiative: Using the Power of Data to Learn From and Improve Outcomes for Every Child and Young Adult With Pediatric Cancer

Data-driven basic, translational, and clinical research has resulted in improved outcomes for children, adolescents, and young adults (AYAs) with pediatric cancers. However, challenges in sharing data between institutions, particularly in research, prevent addressing substantial unmet needs in children and AYA patients diagnosed with certain pediatric cancers. Systematically collecting and sharing data from every child and AYA can enable greater understanding of pediatric cancers, improve survivorship, and accelerate development of new and more effective therapies. To accomplish this goal, the Childhood Cancer Data Initiative (CCDI) was launched in 2019 at the National Cancer Institute. CCDI is a collaborative community endeavor supported by a 10-year, $50-million (in US dollars) annual federal investment. CCDI aims to learn from every patient diagnosed with a pediatric cancer by designing and building a data ecosystem that facilitates data collection, sharing, and analysis for researchers, clinicians, and patients across the cancer community. For example, CCDI's Molecular Characterization Initiative provides comprehensive clinical molecular characterization for children and AYAs with newly diagnosed cancers. Through these efforts, the CCDI strives to provide clinical benefit to patients and improvements in diagnosis and care through data-focused research support and to build expandable, sustainable data resources and workflows to advance research well past the planned 10 years of the initiative. Importantly, if CCDI demonstrates the success of this model for pediatric cancers, similar approaches can be applied to adults, transforming both clinical research and treatment to improve outcomes for all patients with cancer.

Genetically inferred birthweight, height, and puberty timing and risk of osteosarcoma

INTRODUCTION

Several studies have linked increased risk of osteosarcoma with tall stature, high birthweight, and early puberty, although evidence is inconsistent. We used genetic risk scores (GRS) based on established genetic loci for these traits and evaluated associations between genetically inferred birthweight, height, and puberty timing with osteosarcoma.

METHODS

Using genotype data from two genome-wide association studies, totaling 1039 cases and 2923 controls of European ancestry, association analyses were conducted using logistic regression for each study and meta-analyzed to estimate pooled odds ratios (ORs) and 95% confidence intervals (CIs). Subgroup analyses were conducted by case diagnosis age, metastasis status, tumor location, tumor histology, and presence of a known pathogenic variant in a cancer susceptibility gene.

RESULTS

Genetically inferred higher birthweight was associated with an increased risk of osteosarcoma (OR =1.59, 95% CI 1.07-2.38, P = 0.02). This association was strongest in cases without metastatic disease (OR =2.46, 95% CI 1.44-4.19, P = 9.5 ×10-04). Although there was no overall association between osteosarcoma and genetically inferred taller stature (OR=1.06, 95% CI 0.96-1.17, P = 0.28), the GRS for taller stature was associated with an increased risk of osteosarcoma in 154 cases with a known pathogenic cancer susceptibility gene variant (OR=1.29, 95% CI 1.03-1.63, P = 0.03). There were no significant associations between the GRS for puberty timing and osteosarcoma.

CONCLUSION

A genetic propensity to higher birthweight was associated with increased osteosarcoma risk, suggesting that shared genetic factors or biological pathways that affect birthweight may contribute to osteosarcoma pathogenesis.

A genome-wide gene-environment interaction study of breast cancer risk for women of European ancestry

BACKGROUND

Genome-wide studies of gene-environment interactions (G×E) may identify variants associated with disease risk in conjunction with lifestyle/environmental exposures. We conducted a genome-wide G×E analysis of ~ 7.6 million common variants and seven lifestyle/environmental risk factors for breast cancer risk overall and for estrogen receptor positive (ER +) breast cancer.

METHODS

Analyses were conducted using 72,285 breast cancer cases and 80,354 controls of European ancestry from the Breast Cancer Association Consortium. Gene-environment interactions were evaluated using standard unconditional logistic regression models and likelihood ratio tests for breast cancer risk overall and for ER + breast cancer. Bayesian False Discovery Probability was employed to assess the noteworthiness of each SNP-risk factor pairs.

RESULTS

Assuming a 1 × 10-5 prior probability of a true association for each SNP-risk factor pairs and a Bayesian False Discovery Probability < 15%, we identified two independent SNP-risk factor pairs: rs80018847(9p13)-LINGO2 and adult height in association with overall breast cancer risk (ORint = 0.94, 95% CI 0.92-0.96), and rs4770552(13q12)-SPATA13 and age at menarche for ER + breast cancer risk (ORint = 0.91, 95% CI 0.88-0.94).

CONCLUSIONS

Overall, the contribution of G×E interactions to the heritability of breast cancer is very small. At the population level, multiplicative G×E interactions do not make an important contribution to risk prediction in breast cancer.

Probing the diabetes and colorectal cancer relationship using gene - environment interaction analyses

BACKGROUND

Diabetes is an established risk factor for colorectal cancer. However, the mechanisms underlying this relationship still require investigation and it is not known if the association is modified by genetic variants. To address these questions, we undertook a genome-wide gene-environment interaction analysis.

METHODS

We used data from 3 genetic consortia (CCFR, CORECT, GECCO; 31,318 colorectal cancer cases/41,499 controls) and undertook genome-wide gene-environment interaction analyses with colorectal cancer risk, including interaction tests of genetics(G)xdiabetes (1-degree of freedom; d.f.) and joint testing of Gxdiabetes, G-colorectal cancer association (2-d.f. joint test) and G-diabetes correlation (3-d.f. joint test).

RESULTS

Based on the joint tests, we found that the association of diabetes with colorectal cancer risk is modified by loci on chromosomes 8q24.11 (rs3802177, SLC30A8 - ORAA: 1.62, 95% CI: 1.34-1.96; ORAG: 1.41, 95% CI: 1.30-1.54; ORGG: 1.22, 95% CI: 1.13-1.31; p-value3-d.f.: 5.46 × 10-11) and 13q14.13 (rs9526201, LRCH1 - ORGG: 2.11, 95% CI: 1.56-2.83; ORGA: 1.52, 95% CI: 1.38-1.68; ORAA: 1.13, 95% CI: 1.06-1.21; p-value2-d.f.: 7.84 × 10-09).

DISCUSSION

These results suggest that variation in genes related to insulin signaling (SLC30A8) and immune function (LRCH1) may modify the association of diabetes with colorectal cancer risk and provide novel insights into the biology underlying the diabetes and colorectal cancer relationship.

Evaluating approaches for constructing polygenic risk scores for prostate cancer in men of African and European ancestry

Genome-wide polygenic risk scores (GW-PRSs) have been reported to have better predictive ability than PRSs based on genome-wide significance thresholds across numerous traits. We compared the predictive ability of several GW-PRS approaches to a recently developed PRS of 269 established prostate cancer-risk variants from multi-ancestry GWASs and fine-mapping studies (PRS269). GW-PRS models were trained with a large and diverse prostate cancer GWAS of 107,247 cases and 127,006 controls that we previously used to develop the multi-ancestry PRS269. Resulting models were independently tested in 1,586 cases and 1,047 controls of African ancestry from the California Uganda Study and 8,046 cases and 191,825 controls of European ancestry from the UK Biobank and further validated in 13,643 cases and 210,214 controls of European ancestry and 6,353 cases and 53,362 controls of African ancestry from the Million Veteran Program. In the testing data, the best performing GW-PRS approach had AUCs of 0.656 (95% CI = 0.635-0.677) in African and 0.844 (95% CI = 0.840-0.848) in European ancestry men and corresponding prostate cancer ORs of 1.83 (95% CI = 1.67-2.00) and 2.19 (95% CI = 2.14-2.25), respectively, for each SD unit increase in the GW-PRS. Compared to the GW-PRS, in African and European ancestry men, the PRS269 had larger or similar AUCs (AUC = 0.679, 95% CI = 0.659-0.700 and AUC = 0.845, 95% CI = 0.841-0.849, respectively) and comparable prostate cancer ORs (OR = 2.05, 95% CI = 1.87-2.26 and OR = 2.21, 95% CI = 2.16-2.26, respectively). Findings were similar in the validation studies. This investigation suggests that current GW-PRS approaches may not improve the ability to predict prostate cancer risk compared to the PRS269 developed from multi-ancestry GWASs and fine-mapping.

Genome-wide Association Study of Bladder Cancer Reveals New Biological and Translational Insights

BACKGROUND

Genomic regions identified by genome-wide association studies (GWAS) for bladder cancer risk provide new insights into etiology.

OBJECTIVE

To identify new susceptibility variants for bladder cancer in a meta-analysis of new and existing genome-wide genotype data.

DESIGN, SETTING, AND PARTICIPANTS

Data from 32 studies that includes 13,790 bladder cancer cases and 343,502 controls of European ancestry were used for meta-analysis.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSES

Log-additive associations of genetic variants were assessed using logistic regression models. A fixed-effects model was used for meta-analysis of the results. Stratified analyses were conducted to evaluate effect modification by sex and smoking status. A polygenic risk score (PRS) was generated on the basis of known and novel susceptibility variants and tested for interaction with smoking.

RESULTS AND LIMITATIONS

Multiple novel bladder cancer susceptibility loci (6p.22.3, 7q36.3, 8q21.13, 9p21.3, 10q22.1, 19q13.33) as well as improved signals in three known regions (4p16.3, 5p15.33, 11p15.5) were identified, bringing the number of independent markers at genome-wide significance (p < 5 × 10-8) to 24. The 4p16.3 (FGFR3/TACC3) locus was associated with a stronger risk for women than for men (p-interaction = 0.002). Bladder cancer risk was increased by interactions between smoking status and genetic variants at 8p22 (NAT2; multiplicative p value for interaction [pM-I] = 0.004), 8q21.13 (PAG1; pM-I = 0.01), and 9p21.3 (LOC107987026/MTAP/CDKN2A; pM-I = 0.02). The PRS based on the 24 independent GWAS markers (odds ratio per standard deviation increase 1.49, 95% confidence interval 1.44-1.53), which also showed comparable results in two prospective cohorts (UK Biobank, PLCO trial), revealed an approximately fourfold difference in the lifetime risk of bladder cancer according to the PRS (e.g., 1st vs 10th decile) for both smokers and nonsmokers.

CONCLUSIONS

We report novel loci associated with risk of bladder cancer that provide clues to its biological underpinnings. Using 24 independent markers, we constructed a PRS to stratify lifetime risk. The PRS combined with smoking history, and other established risk factors, has the potential to inform future screening efforts for bladder cancer.

PATIENT SUMMARY

We identified new genetic markers that provide biological insights into the genetic causes of bladder cancer. These genetic risk factors combined with lifestyle risk factors, such as smoking, may inform future preventive and screening strategies for bladder cancer.

Evidence of Novel Susceptibility Variants for Prostate Cancer and a Multiancestry Polygenic Risk Score Associated with Aggressive Disease in Men of African Ancestry

BACKGROUND

Genetic factors play an important role in prostate cancer (PCa) susceptibility.

OBJECTIVE

To discover common genetic variants contributing to the risk of PCa in men of African ancestry.

DESIGN, SETTING, AND PARTICIPANTS

We conducted a meta-analysis of ten genome-wide association studies consisting of 19378 cases and 61620 controls of African ancestry.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Common genotyped and imputed variants were tested for their association with PCa risk. Novel susceptibility loci were identified and incorporated into a multiancestry polygenic risk score (PRS). The PRS was evaluated for associations with PCa risk and disease aggressiveness.

RESULTS AND LIMITATIONS

Nine novel susceptibility loci for PCa were identified, of which seven were only found or substantially more common in men of African ancestry, including an African-specific stop-gain variant in the prostate-specific gene anoctamin 7 (ANO7). A multiancestry PRS of 278 risk variants conferred strong associations with PCa risk in African ancestry studies (odds ratios [ORs] >3 and >5 for men in the top PRS decile and percentile, respectively). More importantly, compared with men in the 40-60% PRS category, men in the top PRS decile had a significantly higher risk of aggressive PCa (OR = 1.23, 95% confidence interval = 1.10-1.38, p = 4.4 × 10-4).

CONCLUSIONS

This study demonstrates the importance of large-scale genetic studies in men of African ancestry for a better understanding of PCa susceptibility in this high-risk population and suggests a potential clinical utility of PRS in differentiating between the risks of developing aggressive and nonaggressive disease in men of African ancestry.

PATIENT SUMMARY

In this large genetic study in men of African ancestry, we discovered nine novel prostate cancer (PCa) risk variants. We also showed that a multiancestry polygenic risk score was effective in stratifying PCa risk, and was able to differentiate risk of aggressive and nonaggressive disease.

Germline genetic variants and pediatric rhabdomyosarcoma outcomes: a report from the Children's Oncology Group

BACKGROUND

Relative to other pediatric cancers, survival for rhabdomyosarcoma (RMS) has not improved in recent decades, suggesting the need to enhance risk stratification. Therefore, we conducted a genome-wide association study for event-free survival (EFS) and overall survival (OS) to identify genetic variants associated with outcomes in individuals with RMS.

METHODS

The study included 920 individuals with newly diagnosed RMS who were enrolled in Children's Oncology Group protocols. To assess the association of each single nucleotide polymorphism (SNP) with EFS and OS, we estimated hazard ratios (HRs) and 95% confidence intervals (CIs) using multivariable Cox proportional hazards models, adjusted for clinical covariates. All statistical tests were two sided. We also performed stratified analyses by histological subtype (alveolar and embryonal RMS) and carried out sensitivity analyses of statistically significant SNPs by PAX3/7-FOXO1 fusion status and genetic ancestry group.

RESULTS

We identified that rs17321084 was associated with worse EFS (HR = 2.01, 95% CI = 1.59 to 2.53, P = 5.39 × 10-9) and rs10094840 was associated with worse OS (HR = 1.84, 95% CI = 1.48 to 2.27, P = 2.13 × 10-8). Using publicly available data, we found that rs17321084 lies in a binding region for transcription factors GATA2 and GATA3, and rs10094840 is associated with SPAG1 and RNF19A expression. We also identified that CTNNA3 rs2135732 (HR = 3.75, 95% CI = 2.34 to 5.99, P = 3.54 × 10-8) and MED31 rs74504320 (HR = 3.21, 95% CI = 2.12 to 4.86, P = 3.60 × 10-8) were associated with worse OS among individuals with alveolar RMS.

CONCLUSIONS

We demonstrated that common germline variants are associated with EFS and OS among individuals with RMS. Additional replication and investigation of these SNP effects may further support their consideration in risk stratification protocols.

Genome-wide analyses characterize shared heritability among cancers and identify novel cancer susceptibility regions

BACKGROUND

The shared inherited genetic contribution to risk of different cancers is not fully known. In this study, we leverage results from 12 cancer genome-wide association studies (GWAS) to quantify pairwise genome-wide genetic correlations across cancers and identify novel cancer susceptibility loci.

METHODS

We collected GWAS summary statistics for 12 solid cancers based on 376 759 participants with cancer and 532 864 participants without cancer of European ancestry. The included cancer types were breast, colorectal, endometrial, esophageal, glioma, head and neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancers. We conducted cross-cancer GWAS and transcriptome-wide association studies to discover novel cancer susceptibility loci. Finally, we assessed the extent of variant-specific pleiotropy among cancers at known and newly identified cancer susceptibility loci.

RESULTS

We observed widespread but modest genome-wide genetic correlations across cancers. In cross-cancer GWAS and transcriptome-wide association studies, we identified 15 novel cancer susceptibility loci. Additionally, we identified multiple variants at 77 distinct loci with strong evidence of being associated with at least 2 cancer types by testing for pleiotropy at known cancer susceptibility loci.

CONCLUSIONS

Overall, these results suggest that some genetic risk variants are shared among cancers, though much of cancer heritability is cancer-specific and thus tissue-specific. The increase in statistical power associated with larger sample sizes in cross-disease analysis allows for the identification of novel susceptibility regions. Future studies incorporating data on multiple cancer types are likely to identify additional regions associated with the risk of multiple cancer types.

Moving Toward Findable, Accessible, Interoperable, Reusable Practices in Epidemiologic Research

Data sharing is essential for reproducibility of epidemiologic research, replication of findings, pooled analyses in consortia efforts, and maximizing study value to address multiple research questions. However, barriers related to confidentiality, costs, and incentives often limit the extent and speed of data sharing. Epidemiological practices that follow Findable, Accessible, Interoperable, Reusable (FAIR) principles can address these barriers by making data resources findable with the necessary metadata, accessible to authorized users, and interoperable with other data, to optimize the reuse of resources with appropriate credit to its creators. We provide an overview of these principles and describe approaches for implementation in epidemiology. Increasing degrees of FAIRness can be achieved by moving data and code from on-site locations to remote, accessible ("Cloud") data servers, using machine-readable and nonproprietary files, and developing open-source code. Adoption of these practices will improve daily work and collaborative analyses and facilitate compliance with data sharing policies from funders and scientific journals. Achieving a high degree of FAIRness will require funding, training, organizational support, recognition, and incentives for sharing research resources, both data and code. However, these costs are outweighed by the benefits of making research more reproducible, impactful, and equitable by facilitating the reuse of precious research resources by the scientific community.

Genetically adjusted PSA levels for prostate cancer screening

Prostate-specific antigen (PSA) screening for prostate cancer remains controversial because it increases overdiagnosis and overtreatment of clinically insignificant tumors. Accounting for genetic determinants of constitutive, non-cancer-related PSA variation has potential to improve screening utility. In this study, we discovered 128 genome-wide significant associations (P < 5 × 10-8) in a multi-ancestry meta-analysis of 95,768 men and developed a PSA polygenic score (PGSPSA) that explains 9.61% of constitutive PSA variation. We found that, in men of European ancestry, using PGS-adjusted PSA would avoid up to 31% of negative prostate biopsies but also result in 12% fewer biopsies in patients with prostate cancer, mostly with Gleason score <7 tumors. Genetically adjusted PSA was more predictive of aggressive prostate cancer (odds ratio (OR) = 3.44, P = 6.2 × 10-14, area under the curve (AUC) = 0.755) than unadjusted PSA (OR = 3.31, P = 1.1 × 10-12, AUC = 0.738) in 106 cases and 23,667 controls. Compared to a prostate cancer PGS alone (AUC = 0.712), including genetically adjusted PSA improved detection of aggressive disease (AUC = 0.786, P = 7.2 × 10-4). Our findings highlight the potential utility of incorporating PGS for personalized biomarkers in prostate cancer screening.

Genome-wide association study of lung adenocarcinoma in East Asia and comparison with a European population

Lung adenocarcinoma is the most common type of lung cancer. Known risk variants explain only a small fraction of lung adenocarcinoma heritability. Here, we conducted a two-stage genome-wide association study of lung adenocarcinoma of East Asian ancestry (21,658 cases and 150,676 controls; 54.5% never-smokers) and identified 12 novel susceptibility variants, bringing the total number to 28 at 25 independent loci. Transcriptome-wide association analyses together with colocalization studies using a Taiwanese lung expression quantitative trait loci dataset (n = 115) identified novel candidate genes, including FADS1 at 11q12 and ELF5 at 11p13. In a multi-ancestry meta-analysis of East Asian and European studies, four loci were identified at 2p11, 4q32, 16q23, and 18q12. At the same time, most of our findings in East Asian populations showed no evidence of association in European populations. In our studies drawn from East Asian populations, a polygenic risk score based on the 25 loci had a stronger association in never-smokers vs. individuals with a history of smoking (Pinteraction = 0.0058). These findings provide new insights into the etiology of lung adenocarcinoma in individuals from East Asian populations, which could be important in developing translational applications.

Evaluating Approaches for Constructing Polygenic Risk Scores for Prostate Cancer in Men of African and European Ancestry

Genome-wide polygenic risk scores (GW-PRS) have been reported to have better predictive ability than PRS based on genome-wide significance thresholds across numerous traits. We compared the predictive ability of several GW-PRS approaches to a recently developed PRS of 269 established prostate cancer risk variants from multi-ancestry GWAS and fine-mapping studies (PRS 269 ). GW-PRS models were trained using a large and diverse prostate cancer GWAS of 107,247 cases and 127,006 controls used to develop the multi-ancestry PRS 269 . Resulting models were independently tested in 1,586 cases and 1,047 controls of African ancestry from the California/Uganda Study and 8,046 cases and 191,825 controls of European ancestry from the UK Biobank and further validated in 13,643 cases and 210,214 controls of European ancestry and 6,353 cases and 53,362 controls of African ancestry from the Million Veteran Program. In the testing data, the best performing GW-PRS approach had AUCs of 0.656 (95% CI=0.635-0.677) in African and 0.844 (95% CI=0.840-0.848) in European ancestry men and corresponding prostate cancer OR of 1.83 (95% CI=1.67-2.00) and 2.19 (95% CI=2.14-2.25), respectively, for each SD unit increase in the GW-PRS. However, compared to the GW-PRS, in African and European ancestry men, the PRS 269 had larger or similar AUCs (AUC=0.679, 95% CI=0.659-0.700 and AUC=0.845, 95% CI=0.841-0.849, respectively) and comparable prostate cancer OR (OR=2.05, 95% CI=1.87-2.26 and OR=2.21, 95% CI=2.16-2.26, respectively). Findings were similar in the validation data. This investigation suggests that current GW-PRS approaches may not improve the ability to predict prostate cancer risk compared to the multi-ancestry PRS 269 constructed with fine-mapping.

Polymorphisms within Autophagy-Related Genes as Susceptibility Biomarkers for Multiple Myeloma: A Meta-Analysis of Three Large Cohorts and Functional Characterization

Multiple myeloma (MM) arises following malignant proliferation of plasma cells in the bone marrow, that secrete high amounts of specific monoclonal immunoglobulins or light chains, resulting in the massive production of unfolded or misfolded proteins. Autophagy can have a dual role in tumorigenesis, by eliminating these abnormal proteins to avoid cancer development, but also ensuring MM cell survival and promoting resistance to treatments. To date no studies have determined the impact of genetic variation in autophagy-related genes on MM risk. We performed meta-analysis of germline genetic data on 234 autophagy-related genes from three independent study populations including 13,387 subjects of European ancestry (6863 MM patients and 6524 controls) and examined correlations of statistically significant single nucleotide polymorphisms (SNPs; p < 1 × 10^-9) with immune responses in whole blood, peripheral blood mononuclear cells (PBMCs), and monocyte-derived macrophages (MDM) from a large population of healthy donors from the Human Functional Genomic Project (HFGP). We identified SNPs in six loci, CD46, IKBKE, PARK2, ULK4, ATG5, and CDKN2A associated with MM risk (p = 4.47 × 10^-4-5.79 × 10^-14). Mechanistically, we found that the ULK4_rs6599175 SNP correlated with circulating concentrations of vitamin D3 (p = 4.0 × 10^-4), whereas the IKBKE_rs17433804 SNP correlated with the number of transitional CD24⁺CD38⁺ B cells (p = 4.8 × 10^-4) and circulating serum concentrations of Monocyte Chemoattractant Protein (MCP)-2 (p = 3.6 × 10^-4). We also found that the CD46_rs1142469 SNP correlated with numbers of CD19⁺ B cells, CD19⁺CD3^- B cells, CD5⁺IgD^- cells, IgM^- cells, IgD^-IgM^- cells, and CD4^-CD8^- PBMCs (p = 4.9 × 10^-4-8.6 × 10^-4) and circulating concentrations of interleukin (IL)-20 (p = 0.00082). Finally, we observed that the CDKN2A_rs2811710 SNP correlated with levels of CD4⁺EMCD45RO⁺CD27^- cells (p = 9.3 × 10^-4). These results suggest that genetic variants within these six loci influence MM risk through the modulation of specific subsets of immune cells, as well as vitamin D3^-, MCP-2^-, and IL20-dependent pathways.

Association between circulating inflammatory markers and adult cancer risk: a Mendelian randomization analysis

Background

Tumour-promoting inflammation is a "hallmark" of cancer and conventional epidemiological studies have reported links between various inflammatory markers and cancer risk. The causal nature of these relationships and, thus, the suitability of these markers as intervention targets for cancer prevention is unclear.

Methods

We meta-analysed 6 genome-wide association studies of circulating inflammatory markers comprising 59,969 participants of European ancestry. We then used combined cis-Mendelian randomization and colocalisation analysis to evaluate the causal role of 66 circulating inflammatory markers in risk of 30 adult cancers in 338,162 cancer cases and up to 824,556 controls. Genetic instruments for inflammatory markers were constructed using genome-wide significant (P < 5.0 x 10-8) cis-acting SNPs (i.e. in or ±250 kb from the gene encoding the relevant protein) in weak linkage disequilibrium (LD, r2 < 0.10). Effect estimates were generated using inverse-variance weighted random-effects models and standard errors were inflated to account for weak LD between variants with reference to the 1000 Genomes Phase 3 CEU panel. A false discovery rate (FDR)-corrected P-value ("q-value") < 0.05 was used as a threshold to define "strong evidence" to support associations and 0.05 ≤ q-value < 0.20 to define "suggestive evidence". A colocalisation posterior probability (PPH4) > 70% was employed to indicate support for shared causal variants across inflammatory markers and cancer outcomes.

Results

We found strong evidence to support an association of genetically-proxied circulating pro-adrenomedullin concentrations with increased breast cancer risk (OR 1.19, 95% CI 1.10-1.29, q-value=0.033, PPH4=84.3%) and suggestive evidence to support associations of interleukin-23 receptor concentrations with increased pancreatic cancer risk (OR 1.42, 95% CI 1.20-1.69, q-value=0.055, PPH4=73.9%), prothrombin concentrations with decreased basal cell carcinoma risk (OR 0.66, 95% CI 0.53-0.81, q-value=0.067, PPH4=81.8%), macrophage migration inhibitory factor concentrations with increased bladder cancer risk (OR 1.14, 95% CI 1.05-1.23, q-value=0.072, PPH4=76.1%), and interleukin-1 receptor-like 1 concentrations with decreased triple-negative breast cancer risk (OR 0.92, 95% CI 0.88-0.97, q-value=0.15), PPH4=85.6%). For 22 of 30 cancer outcomes examined, there was little evidence (q-value ≥ 0.20) that any of the 66 circulating inflammatory markers examined were associated with cancer risk.

Conclusion

Our comprehensive joint Mendelian randomization and colocalisation analysis of the role of circulating inflammatory markers in cancer risk identified potential roles for 5 circulating inflammatory markers in risk of 5 site-specific cancers. Contrary to reports from some prior conventional epidemiological studies, we found little evidence of association of circulating inflammatory markers with the majority of site-specific cancers evaluated.

Oncogenic CDK13 mutations impede nuclear RNA surveillance

RNA surveillance pathways detect and degrade defective transcripts to ensure RNA fidelity. We found that disrupted nuclear RNA surveillance is oncogenic. Cyclin-dependent kinase 13 (CDK13) is mutated in melanoma, and patient-mutated CDK13 accelerates zebrafish melanoma. CDK13 mutation causes aberrant RNA stabilization. CDK13 is required for ZC3H14 phosphorylation, which is necessary and sufficient to promote nuclear RNA degradation. Mutant CDK13 fails to activate nuclear RNA surveillance, causing aberrant protein-coding transcripts to be stabilized and translated. Forced aberrant RNA expression accelerates melanoma in zebrafish. We found recurrent mutations in genes encoding nuclear RNA surveillance components in many malignancies, establishing nuclear RNA surveillance as a tumor-suppressive pathway. Activating nuclear RNA surveillance is crucial to avoid accumulation of aberrant RNAs and their ensuing consequences in development and disease.

Abstract 3508: Association between clonal hematopoiesis and risk of prostate cancer in a large sample of African ancestry men

Clonal hematopoiesis of indeterminate potential (CHIP) has been associated with inflammation, which is a risk factor for cancer, including prostate cancer. We previously reported weak evidence of an association between CHIP and prostate cancer risk in men of European ancestry. However, little is known for African ancestry populations. We investigated the association of age-related CHIP with overall and aggressive prostate cancer risk in a large whole-exome sequencing study of 12,049 African ancestry men, including 7,176 prostate cancer cases (of which 3,283 had aggressive disease and 1,074 had metastatic disease) and 4,873 controls. Somatic variant calling was carried out using GATK Mutect2, and only variants with minor allele frequencies (MAF) &lt;0.1% and a variant allelic fraction (VAF) &gt;5% were included. Variants with a MAF ≥0.1% in gnomAD were excluded. CHIP variants were identified from a list of pre-specified mutations in 74 genes. Associations were tested using regression models adjusting for age, sub-study, and top 10 principal components, with statistical significance tested by the likelihood ratio test and applying a Bonferroni correction to account for multiple testing. In total, 998 variants in 57 CHIP genes were identified. Consistent with previous results, we observed a strong association between CHIP and age at blood draw. CHIP genes in aggregate were not statistically significantly associated with risks of total (OR=1.12, 95% CI=0.97-1.28), aggressive (OR=1.14, 95% CI=0.92-1.43) or metastatic (OR=1.17, 95% CI=0.91-1.49) prostate cancer. We observed that carriers of variants in DNMT3A, which is the gene that harbors the most CHIP driver mutations, had a nominally elevated risk of prostate cancer compared to non-carriers (OR=1.35, 95% CI=1.08-1.68, p=0.007). Additionally, carriers of variants in EZH2, which is implicated in cancer progression, showed a suggestive association with aggressive prostate cancer (OR=7.33, 95% CI=1.01-53.21, p=0.029). After adjusting for age at blood draw, CHIP genes in aggregate were not associated with age at prostate cancer diagnosis. However, we found that EZH2 variants carriers were diagnosed 12.9 years earlier on average than non-carriers (95% CI=6.1-19.7, adjusted p=0.01). A prostate cancer polygenic risk score (PRS) constructed using 269 risk variants was not associated with CHIP carrier status in aggregate (OR=0.99, 95% CI=0.92-1.06, p=0.70) or with any individual gene (all adjusted p&gt;0.05). In summary, overall CHIP is not likely to be a risk factor of prostate cancer or aggressive disease in men of African ancestry. However, our results do confirm the association of CHIP in DNMT3A with prostate cancer risk as reported in previous studies in men of European ancestry. Future work will be needed to evaluate the biological causality of DNMT3A- and EZH2- related CHIP on prostate cancer.

Citation Format: Anqi Wang, Yili Xu, Xin Sheng, Raymond Hughley, Ben Adusei, Mohamed Jalloh, Serigne Magueye Gueye, Andrew A Adjei, James Mensah, Pedro W. Fernandez, Akin Olupelumi Adebiyi, Oseremen Inokhoife Aisuodionoe-Shadrach, Lindsay Petersen, Maureen Joffe, Jeannette T. Bensen, James L. Mohler, Jack A. Taylor, Eboneé N. Butler, Sue A. Ingles, Benjamin A. Rybicki, Janet L. Stanford, Wei Zheng, Sonja I. Berndt, Chad D. Huff, Joseph Lachance, Luc Multigner, Caroline Andrews, Timothy R. Rebbeck, Laurent Brureau, Stephen J. Chanock, Adam de Smith, Fei Chen, Burcu F. Darst, David V. Conti, Christopher A. Haiman. Association between clonal hematopoiesis and risk of prostate cancer in a large sample of African ancestry men [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3508.

Abstract 5722: The mutational signatures of 100,477 targeted sequenced tumors

Mutational signatures, the footprints of somatic mutations, are associated with the causes of cancer and have been well-studied for tumors with whole exome or genome sequencing. However, due to the low number of detected mutations, mutational signatures were insufficiently explored for many tumors sequenced by targeted panels in clinics. It impedes the clinical application of mutational signatures. Here, we present a new method, SALMON (Signature Analyzer for Low Mutation cOuNts), to identify mutational signatures in targeted sequenced tumors based on tumor mutational burden. SALMON adjusts for panel size differences and uses a large number of targeted sequenced tumors rather than a large number of mutations per tumor (as with whole exome or genome sequencing) to overcome the challenges of mutational signature analysis for targeted sequencing. Extensive simulations and pseudo-targeted sequenced data show that SALMON can accurately detect spiky or common signatures. We applied SALMON to investigate the pan-cancer patterns of mutational signatures for 100,477 targeted sequenced tumors in AACR Project GENIE, including 14,428 lung and 11,389 breast tumors. We detected well-established signatures in tumor types that have not previously been associated with these signatures, such as the smoking signature in ovarian tumors. Interestingly, analysis of thousands of tumors per cancer type from diverse populations revealed gender discrepancies, self-described race differences, subtype heterogeneity, and metastatic enrichment of mutational signatures. For instance, most sex-biased signatures are more frequently present in males for non-gender-specific cancers. Thiopurine treatment-induced signature in glioma is enriched in Black patients. And endometrioid ovarian or uterine cancers have a higher prevalence of polymerase epsilon (POLE) deficiency-related signatures than non-endometrioid ovarian or uterine cancers, respectively. Our study demonstrates the feasibility and utility of mutational signature analysis for targeted sequenced tumors, enabling precision applications of mutational signatures in the clinical setting.

Citation Format: Donghyuk Lee, Min Hua, Difei Wang, Lei Song, Tongwu Zhang, Kai Yu, Xiaohong R. Yang, Jianxin Shi, Stephen J. Chanock, Maria Teresa Landi, Bin Zhu. The mutational signatures of 100,477 targeted sequenced tumors. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 5722.

Targeted long-read sequencing of the Ewing sarcoma 6p25.1 susceptibility locus identifies germline-somatic interactions with EWSR1-FLI1 binding

Ewing sarcoma (EwS) is a rare bone and soft tissue malignancy driven by chromosomal translocations encoding chimeric transcription factors, such as EWSR1-FLI1, that bind GGAA motifs forming novel enhancers that alter nearby expression. We propose that germline microsatellite variation at the 6p25.1 EwS susceptibility locus could impact downstream gene expression and EwS biology. We performed targeted long-read sequencing of EwS blood DNA to characterize variation and genomic features important for EWSR1-FLI1 binding. We identified 50 microsatellite alleles at 6p25.1 and observed that EwS-affected individuals had longer alleles (>135 bp) with more GGAA repeats. The 6p25.1 GGAA microsatellite showed chromatin features of an EWSR1-FLI1 enhancer and regulated expression of RREB1, a transcription factor associated with RAS/MAPK signaling. RREB1 knockdown reduced proliferation and clonogenic potential and reduced expression of cell cycle and DNA replication genes. Our integrative analysis at 6p25.1 details increased binding of longer GGAA microsatellite alleles with acquired EWSR-FLI1 to promote Ewing sarcomagenesis by RREB1-mediated proliferation.

A Multimedia Strategy to Integrate Introductory Broad-Based Radiation Science Education in US Medical Schools

US physicians in multiple specialties who order or conduct radiological procedures lack formal radiation science education and thus sometimes order procedures of limited benefit or fail to order what is necessary. To this end, a multidisciplinary expert group proposed an introductory broad-based radiation science educational program for US medical schools. Suggested preclinical elements of the curriculum include foundational education on ionizing and nonionizing radiation (eg, definitions, dose metrics, and risk measures) and short- and long-term radiation-related health effects as well as introduction to radiology, radiation therapy, and radiation protection concepts. Recommended clinical elements of the curriculum would impart knowledge and practical experience in radiology, fluoroscopically guided procedures, nuclear medicine, radiation oncology, and identification of patient subgroups requiring special considerations when selecting specific ionizing or nonionizing diagnostic or therapeutic radiation procedures. Critical components of the clinical program would also include educational material and direct experience with patient-centered communication on benefits of, risks of, and shared decision making about ionizing and nonionizing radiation procedures and on health effects and safety requirements for environmental and occupational exposure to ionizing and nonionizing radiation. Overarching is the introduction to evidence-based guidelines for procedures that maximize clinical benefit while limiting unnecessary risk. The content would be further developed, directed, and integrated within the curriculum by local faculties and would address multiple standard elements of the Liaison Committee on Medical Education and Core Entrustable Professional Activities for Entering Residency of the Association of American Medical Colleges.

Population analyses of mosaic X chromosome loss identify genetic drivers and widespread signatures of cellular selection

Mosaic loss of the X chromosome (mLOX) is the most commonly occurring clonal somatic alteration detected in the leukocytes of women, yet little is known about its genetic determinants or phenotypic consequences. To address this, we estimated mLOX in >900,000 women across eight biobanks, identifying 10% of women with detectable X loss in approximately 2% of their leukocytes. Out of 1,253 diseases examined, women with mLOX had an elevated risk of myeloid and lymphoid leukemias and pneumonia. Genetic analyses identified 49 common variants influencing mLOX, implicating genes with established roles in chromosomal missegregation, cancer predisposition, and autoimmune diseases. Complementary exome-sequence analyses identified rare missense variants in FBXO10 which confer a two-fold increased risk of mLOX. A small fraction of these associations were shared with mosaic Y chromosome loss in men, suggesting different biological processes drive the formation and clonal expansion of sex chromosome missegregation events. Allelic shift analyses identified alleles on the X chromosome which are preferentially retained, demonstrating that variation at many loci across the X chromosome is under cellular selection. A novel polygenic score including 44 independent X chromosome allelic shift loci correctly inferred the retained X chromosomes in 80.7% of mLOX cases in the top decile. Collectively our results support a model where germline variants predispose women to acquiring mLOX, with the allelic content of the X chromosome possibly shaping the magnitude of subsequent clonal expansion.

FANCM missense variants and breast cancer risk: a case-control association study of 75,156 European women

Evidence from literature, including the BRIDGES study, indicates that germline protein truncating variants (PTVs) in FANCM confer moderately increased risk of ER-negative and triple-negative breast cancer (TNBC), especially for women with a family history of the disease. Association between FANCM missense variants (MVs) and breast cancer risk has been postulated. In this study, we further used the BRIDGES study to test 689 FANCM MVs for association with breast cancer risk, overall and in ER-negative and TNBC subtypes, in 39,885 cases (7566 selected for family history) and 35,271 controls of European ancestry. Sixteen common MVs were tested individually; the remaining rare 673 MVs were tested by burden analyses considering their position and pathogenicity score. We also conducted a meta-analysis of our results and those from published studies. We did not find evidence for association for any of the 16 variants individually tested. The rare MVs were significantly associated with increased risk of ER-negative breast cancer by burden analysis comparing familial cases to controls (OR = 1.48; 95% CI 1.07-2.04; P = 0.017). Higher ORs were found for the subgroup of MVs located in functional domains or predicted to be pathogenic. The meta-analysis indicated that FANCM MVs overall are associated with breast cancer risk (OR = 1.22; 95% CI 1.08-1.38; P = 0.002). Our results support the definition from previous analyses of FANCM as a moderate-risk breast cancer gene and provide evidence that FANCM MVs could be low/moderate risk factors for ER-negative and TNBC subtypes. Further genetic and functional analyses are necessary to clarify better the increased risks due to FANCM MVs.

Aggregation tests identify new gene associations with breast cancer in populations with diverse ancestry

BACKGROUND

Low-frequency variants play an important role in breast cancer (BC) susceptibility. Gene-based methods can increase power by combining multiple variants in the same gene and help identify target genes.

METHODS

We evaluated the potential of gene-based aggregation in the Breast Cancer Association Consortium cohorts including 83,471 cases and 59,199 controls. Low-frequency variants were aggregated for individual genes' coding and regulatory regions. Association results in European ancestry samples were compared to single-marker association results in the same cohort. Gene-based associations were also combined in meta-analysis across individuals with European, Asian, African, and Latin American and Hispanic ancestry.

RESULTS

In European ancestry samples, 14 genes were significantly associated (q < 0.05) with BC. Of those, two genes, FMNL3 (P = 6.11 × 10-6) and AC058822.1 (P = 1.47 × 10-4), represent new associations. High FMNL3 expression has previously been linked to poor prognosis in several other cancers. Meta-analysis of samples with diverse ancestry discovered further associations including established candidate genes ESR1 and CBLB. Furthermore, literature review and database query found further support for a biologically plausible link with cancer for genes CBLB, FMNL3, FGFR2, LSP1, MAP3K1, and SRGAP2C.

CONCLUSIONS

Using extended gene-based aggregation tests including coding and regulatory variation, we report identification of plausible target genes for previously identified single-marker associations with BC as well as the discovery of novel genes implicated in BC development. Including multi ancestral cohorts in this study enabled the identification of otherwise missed disease associations as ESR1 (P = 1.31 × 10-5), demonstrating the importance of diversifying study cohorts.