Matsuo, Keitaro

Background: Breast cancer polygenic risk scores (PRS) and traditional risk models (e.g., the Gail model [Gail]) are known to contribute largely independent information, but it is unclear how the overlap varies by ancestry, age, disease type (invasive breast cancer, DCIS), and risk threshold. Methods: In a retrospective case–control study, we evaluated risk prediction performance in 180,398 women (161,849 of European ancestry; 18,549 of Asian ancestry). Odds ratios (ORs) from logistic regression models and the area under the receiver operating characteristic curve (AUC) were estimated. Results: PRS for invasive disease showed a stronger association in younger (<50 years) women (OR = 2.51, AUC = 0.622) than in women ≥ 50 years (OR = 2.06, AUC = 0.653) of European ancestry. PRS performance in Asians was lower (OR range = 1.62–1.64, AUC = 0.551–0.600). Gail performance was modest across groups and poor in younger Asian women (OR = 0.94–0.99, AUC = 0.523–0.533). Age interactions were observed for both PRS (p < 0.001) and Gail (p < 0.001) in Europeans, whereas in Asians, age interaction was observed only for Gail (invasive: p < 0.001; DCIS: p = 0.002). PRS identified more high-risk individuals than Gail in Asian populations, especially ≥50 years, while Gail identified more in Europeans. Overlap between PRS, Gail, and family history was limited at higher thresholds. Calibration analysis, comparing empirical and model-based ROC curves, showed divergence for both PRS and Gail (p < 0.001), which indicates miscalibration. In Europeans, family history and prior biopsies drove Gail discrimination. In younger Asians, age at first live birth was influential. Conclusions: PRS adds value to risk stratification beyond traditional tools, especially in younger women and Asian ancestry populations.

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

Abstract Background The role of ovulation in epithelial ovarian cancer (EOC) is supported by the consistent protective effects of parity and oral contraceptive use. Whether these factors protect through anovulation alone remains unclear. We explored the association between lifetime ovulatory years (LOY) and EOC. Methods LOY was calculated using 12 algorithms. Odds ratios (ORs) and 95% confidence intervals (CIs) estimated the association between LOY or LOY components and EOC among 26 204 control participants and 21 267 case patients from 25 studies. To assess whether LOY components act through ovulation suppression alone, we compared beta coefficients obtained from regression models with expected estimates assuming 1 year of ovulation suppression has the same effect regardless of source. Results LOY was associated with increased EOC risk (OR per year increase = 1.014, 95% CI = 1.009 to 1.020 to OR per year increase = 1.044, 95% CI = 1.041 to 1.048). Individual LOY components, except age at menarche, also associated with EOC. The estimated model coefficient for oral contraceptive use and pregnancies were 4.45 times and 12- to 15-fold greater than expected, respectively. LOY was associated with high-grade serous, low-grade serous, endometrioid, and clear cell histotypes (ORs per year increase = 1.054, 1.040, 1.065, and 1.098, respectively) but not mucinous tumors. Estimated coefficients of LOY components were close to expected estimates for high-grade serous but larger than expected for low-grade serous, endometrioid, and clear cell histotypes. Conclusions LOY is positively associated with nonmucinous EOC. Differences between estimated and expected model coefficients for LOY components suggest factors beyond ovulation underlie the associations between LOY components and EOC in general and for non-HGSOC.

Abstract Background: Ovarian cancer is the deadliest gynecologic cancer. Standard treatments for advanced stage high-grade serous ovarian cancer, the most common type, include (1) primary cytoreductive surgery (PCS) followed by adjuvant chemotherapy or (2) neoadjuvant chemotherapy (NACT) with interval debulking surgery. Patients in whom PCS is unlikely to yield optimal cytoreduction to no visible residual disease (R0) or who have medical contraindications to PCS are recommended to undergo NACT. Previous studies on associations between environmental factors and risk of any macroscopic residual disease have yielded inconsistent results. We aimed to (1) comprehensively examine the associations between demographic, lifestyle and reproductive factors and risk of residual disease after PCS; and (2) develop and internally validate a risk prediction model based on these factors. Methods: We used pooled data on 3,492 women who had PCS following a diagnosis with advanced stage high-grade serous ovarian, fallopian tube or primary peritoneal cancers from the Ovarian Cancer Association Consortium. Fifteen exposures of interest included age at diagnosis, menopausal status, race/ethnicity, education level, first-degree family history of ovarian cancer, endometriosis, smoking, body mass index, parity, incomplete pregnancy, tubal ligation, use of combined oral contraceptives, depot-medroxyprogesterone acetate, estrogen (ET), and combined estrogen-progestin therapy. We fit logistic regression models to examine each exposure-residual disease association in 80% of the data (n=2,794; random selection of participants). We developed a risk prediction model including the factors with p-values ≤0.2. Area under the receiver operating characteristic curve (AUC) was used to assess the model’s discrimination in the remaining 20% of the data (n=698). Results: Of the 3,492 participants, 2,003 (57%) had residual disease following PCS. Older age at diagnosis was associated with an increased risk of residual disease (OR=1.07 per five years, 95% CI 1.03-1.11). In contrast, a family history of ovarian cancer (OR=0.63, 95% CI 0.42-0.95) or ET use for 5+ years (OR=0.61, 95% CI 0.39-0.96) were associated with achieving R0. These above factors were included in the risk prediction model as were race/ethnicity, personal history of endometriosis, and smoking, whose p-values ≤0.2. The model showed modest performance in the validation set (AUC=0.64). Conclusions: Younger age at diagnosis, family history of ovarian cancer, and long-term ET use were associated with achieving R0 following PCS. Future studies incorporating genetic and clinical factors to improve the risk prediction for residual disease following PCS are warranted. Citation Format: Minh Tung Phung, Andrew Berchuck, Ellen L. Goode, Marc T. Goodman, Gillian E. Hanley, Jean Richardson, Bronwyn Grout, Anne Chase, Cindy McKinnon Deurloo, Beth Y. Karlan, Toon Van Gorp, Keitaro Matsuo, Karen McLean, Malcolm C. Pike, Joellen M. Schildkraut, Kathryn L. Terry, Anna DeFazio, Penelope M. Webb, Paul D. P. Pharoah, Susan J. Ramus, Celeste Leigh Pearce. Environmental factors associated with residual disease after ovarian cancer primary cytoreduction surgery [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 2248.

Objective: Studies on low-grade serous ovarian cancer (LGSC) are limited by a low number of cases. The aim of this study was to define the prognostic significance of age, stage, and CA-125 levels on survival in a multi-institutional cohort of women with pathologically confirmed LGSC. Methods: Women with LGSC were identified from the collaborative Ovarian Cancer Association Consortium (OCAC). Cases of newly diagnosed primary LGSC were included if peri-operative CA-125 levels were available. Age at diagnosis, FIGO stage, pre- and post-treatment CA-125 levels, residual disease, adjuvant chemotherapy, disease recurrence, and vital status were collected by the participating institutions. Progression-free (PFS) and overall survival (OS) were calculated. Multivariable (MVA) Cox proportional hazard models were used and hazard ratios (HR) calculated. Results: A total of 176 women with LGSC were included in this study; 82% had stage III/IV disease. The median PFS was 2.3 years and the median OS was 6.4 years. Age at diagnosis was not significantly associated with worse PFS (p = 0.23) or OS (p = 0.3) (HR per year: 0.99; 95%CI, 0.96–1.01 and 0.98; 95%CI 0.95–1.01). FIGO stage III/IV was independently associated with PFS (HR 4.26, 95%CI 1.43–12.73) and OS (HR 1.69, 95%CI 0.56–5.05). Elevated CA-125 (≥35 U/mL) at diagnosis was not significantly associated with worse PFS (p = 0.87) or OS (p = 0.78) in MVA. Elevated CA-125 (≥35 U/mL) after completion of primary treatment was independently associated with worse PFS (HR 2.81, 95%CI 1.36–5.81) and OS (HR 6.62, 95%CI 2.45–17.92). In the MVA, residual disease was independently associated with PFS (0.022), but not OS (0.85). Conclusion: Advanced LGSC was associated with poor long-term prognosis. FIGO stage and abnormal post-treatment CA-125 level are key prognostic factors inversely associated with PFS and OS. Highlights: 1. Through a multi-center collaborative effort, data from 176 women with low-grade serous ovarian cancer were analyzed. 2. Although low-grade serous ovarian cancer is often considered indolent, the progression-free and overall survival are poor. 3. Elevated post-treatment CA-125 levels are independently associated with poor survival.

In this study we aim to examine gene–environment interactions (GxEs) between genes involved with estrogen metabolism and environmental factors related to estrogen exposure. GxE analyses were conducted with 1970 Korean breast cancer cases and 2052 controls in the case-control study, the Seoul Breast Cancer Study (SEBCS). A total of 11,555 SNPs from the 137 candidate genes were included in the GxE analyses with eight established environmental factors. A replication test was conducted by using an independent population from the Breast Cancer Association Consortium (BCAC), with 62,485 Europeans and 9047 Asians. The GxE tests were performed by using two-step methods in GxEScan software. Two interactions were found in the SEBCS. The first interaction was shown between rs13035764 of NCOA1 and age at menarche in the GE|2df model (p-2df = 1.2 × 10−3). The age at menarche before 14 years old was associated with the high risk of breast cancer, and the risk was higher when subjects had homozygous minor allele G. The second GxE was shown between rs851998 near ESR1 and height in the GE|2df model (p-2df = 1.1 × 10−4). Height taller than 160 cm was associated with a high risk of breast cancer, and the risk increased when the minor allele was added. The findings were not replicated in the BCAC. These results would suggest specificity in Koreans for breast cancer risk.

Abstract Background: It is not known whether modifiable lifestyle factors that predict survival after invasive breast cancer differ by subtype. Methods: We analyzed data for 121,435 women diagnosed with breast cancer from 67 studies in the Breast Cancer Association Consortium with 16,890 deaths (8,554 breast cancer specific) over 10 years. Cox regression was used to estimate associations between risk factors and 10-year all-cause mortality and breast cancer–specific mortality overall, by estrogen receptor (ER) status, and by intrinsic-like subtype. Results: There was no evidence of heterogeneous associations between risk factors and mortality by subtype (Padj &gt; 0.30). The strongest associations were between all-cause mortality and BMI ≥30 versus 18.5–25 kg/m2 [HR (95% confidence interval (CI), 1.19 (1.06–1.34)]; current versus never smoking [1.37 (1.27–1.47)], high versus low physical activity [0.43 (0.21–0.86)], age ≥30 years versus &lt;20 years at first pregnancy [0.79 (0.72–0.86)]; &gt;0–&lt;5 years versus ≥10 years since last full-term birth [1.31 (1.11–1.55)]; ever versus never use of oral contraceptives [0.91 (0.87–0.96)]; ever versus never use of menopausal hormone therapy, including current estrogen–progestin therapy [0.61 (0.54–0.69)]. Similar associations with breast cancer mortality were weaker; for example, 1.11 (1.02–1.21) for current versus never smoking. Conclusions: We confirm associations between modifiable lifestyle factors and 10-year all-cause mortality. There was no strong evidence that associations differed by ER status or intrinsic-like subtype. Impact: Given the large dataset and lack of evidence that associations between modifiable risk factors and 10-year mortality differed by subtype, these associations could be cautiously used in prognostication models to inform patient-centered care.

Abstract Background: Lung cancer remains the leading cause of cancer death worldwide, with 15% to 20% occurring in never smokers. To assess genetic determinants for prognosis among never smokers, we conducted a genome-wide investigation in the International Lung Cancer Consortium (ILCCO). Methods: Genomic and clinical data from 1,569 never-smoking patients with lung cancer of European ancestry from 10 ILCCO studies were included. HRs and 95% confidence intervals of overall survival were estimated. We assessed whether the associations were mediated through mRNA expression–based 1,553 normal lung tissues from the lung expression quantitative trait loci (eQTL) dataset and Genotype-Tissue Expression (GTEx). For cross-ethnicity generalization, we assessed the associations in a Japanese study (N = 887). Results: One locus at 13q22.2 was associated with lung adenocarcinoma survival at genome-wide level, with carriers of rs12875562-T allele exhibiting poor prognosis [HR = 1.71 (1.41–2.07), P = 3.60 × 10−8], and altered mRNA expression of LMO7DN in lung tissue (GTEx, P = 9.40 × 10−7; Lung eQTL dataset, P = 0.003). Furthermore, 2 of 11 independent loci that reached the suggestive significance level (P &lt; 10−6) were significant eQTL affecting mRNA expression of nearby genes in lung tissues, including CAPZB at 1p36.13 and UBAC1 at 9q34.3. One locus encoding NWD2/KIAA1239 at 4p14 showed associations in both European [HR = 0.50 (0.38–0.66), P = 6.92 × 10−7] and Japanese populations [HR = 0.79 (0.67–0.94), P = 0.007]. Conclusions: Based on the largest genomic investigation on the lung cancer prognosis of never smokers to date, we observed that lung cancer prognosis is affected by inherited genetic variants. Impact: We identified one locus near LMO7DN at genome-wide level and several potential prognostic genes with cis-effect on mRNA expression. Further functional genomics work is required to understand their role in tumor progression.

Abstract DNA methylation is instrumental for gene regulation. Global changes in the epigenetic landscape have been recognized as a hallmark of cancer. However, the role of DNA methylation in epithelial ovarian cancer (EOC) remains unclear. In this study, high-density genetic and DNA methylation data in white blood cells from the Framingham Heart Study (N = 1,595) were used to build genetic models to predict DNA methylation levels. These prediction models were then applied to the summary statistics of a genome-wide association study (GWAS) of ovarian cancer including 22,406 EOC cases and 40,941 controls to investigate genetically predicted DNA methylation levels in association with EOC risk. Among 62,938 CpG sites investigated, genetically predicted methylation levels at 89 CpG were significantly associated with EOC risk at a Bonferroni-corrected threshold of P &lt; 7.94 × 10−7. Of them, 87 were located at GWAS-identified EOC susceptibility regions and two resided in a genomic region not previously reported to be associated with EOC risk. Integrative analyses of genetic, methylation, and gene expression data identified consistent directions of associations across 12 CpG, five genes, and EOC risk, suggesting that methylation at these 12 CpG may influence EOC risk by regulating expression of these five genes, namely MAPT, HOXB3, ABHD8, ARHGAP27, and SKAP1. We identified novel DNA methylation markers associated with EOC risk and propose that methylation at multiple CpG may affect EOC risk via regulation of gene expression. Significance: Identification of novel DNA methylation markers associated with EOC risk suggests that methylation at multiple CpG may affect EOC risk through regulation of gene expression.

Abstract Genome-wide association studies have identified 40 ovarian cancer risk loci. However, the mechanisms underlying these associations remain elusive. In this study, we conducted a two-pronged approach to identify candidate causal SNPs and assess underlying biological mechanisms at chromosome 9p22.2, the first and most statistically significant associated locus for ovarian cancer susceptibility. Three transcriptional regulatory elements with allele-specific effects and a scaffold/matrix attachment region were characterized and, through physical DNA interactions, BNC2 was established as the most likely target gene. We determined the consensus binding sequence for BNC2 in vitro, verified its enrichment in BNC2 ChIP-seq regions, and validated a set of its downstream target genes. Fine-mapping by dense regional genotyping in over 15,000 ovarian cancer cases and 30,000 controls identified SNPs in the scaffold/matrix attachment region as among the most likely causal variants. This study reveals a comprehensive regulatory landscape at 9p22.2 and proposes a likely mechanism of susceptibility to ovarian cancer. Significance: Mapping the 9p22.2 ovarian cancer risk locus identifies BNC2 as an ovarian cancer risk gene. See related commentary by Choi and Brown, p. 439