Hopper, John L.

Abstract Studies have identified genetic and epidemiologic factors associated with mammographic density (MD) phenotypes. However, MD-associated genetic variants only account for a small proportion of the total estimated heritability. Interrogating interactions between genetic and epidemiologic factors could potentially identify additional MD-associated loci, expand our understanding of the genetic basis of MD phenotypes, and clarify how epidemiologic factors modulate relationships between genetic variants and MD. We conducted six separate genome-wide, gene–environment (GxE) interaction analyses, applying 2 degrees of freedom (df) and 1df interaction tests, for each of three MD phenotypes (percent density, dense area (DA), and nondense area (NDA)). The six epidemiologic factors considered were height, ever parous, parity, ever menopausal hormone therapy, ever breastfeeding, and months of breastfeeding. We included European ancestry participants from multiple studies within the Markers of Density consortium and the Breast Cancer Association Consortium (n = 4895-16 218 depending on specific analyses). We identified 11 loci with genome-wide significant (P < 5 × 10−8) interaction tests including two novel common genetic signals interacting with parity (8p21.2) and ever breastfeeding (19p13.2) for NDA. Our results suggest that epidemiologic risk factors might influence relationships between common genetic variants and MD phenotypes at particular genomic loci.

Abstract Menopausal users of hormone replacement therapy (HRT) are at increased breast cancer risk and decreased colorectal cancer (CRC) risk compared with individuals who have never used HRT, but these opposing associations may differ by familial risk of breast cancer and CRC. We harmonized data from 3 cohorts and generated separate breast cancer and CRC familial risk scores based on cancer family history. We defined moderate or strong family history as a risk score of 0.4 or higher, where 0.4 was equivalent to a 50-year-old woman with 1 parent diagnosed with either breast cancer or CRC at 55 years of age. Of 24 486 women assessed, 1243 and 405 were diagnosed with incident breast cancer and CRC, respectively. For breast cancer, menopausal HRT ever use versus never use hazard ratios were 1.27 (95% CI = 1.11 to 1.45) for a breast cancer familial risk score below 0.4 and 1.01 (95% CI = 0.82 to 1.25) for a breast cancer familial risk score of 0.4 or higher (Pdifference = .08). For CRC, menopausal HRT hazard ratios were 0.63 (95% CI = 0.50 to 0.78) for a CRC familial risk score below 0.4 and 1.21 (95% CI = 0.73 to 2.00) for a CRC familial risk score of 0.4 or higher (Pdifference = .03). Associations with menopausal HRT use that apply to the general population may not hold for women at moderate or strong familial risk of these cancers.

A large number of variants identified through clinical genetic testing in disease susceptibility genes are of uncertain significance (VUS). Following the recommendations of the American College of Medical Genetics and Genomics (ACMG) and Association for Molecular Pathology (AMP), the frequency in case-control datasets (PS4 criterion) can inform their interpretation. We present a novel case-control likelihood ratio-based method that incorporates gene-specific age-related penetrance. We demonstrate the utility of this method in the analysis of simulated and real datasets. In the analysis of simulated data, the likelihood ratio method was more powerful compared to other methods. Likelihood ratios were calculated for a case-control dataset of BRCA1 and BRCA2 variants from the Breast Cancer Association Consortium (BCAC) and compared with logistic regression results. A larger number of variants reached evidence in favor of pathogenicity, and a substantial number of variants had evidence against pathogenicity—findings that would not have been reached using other case-control analysis methods. Our novel method provides greater power to classify rare variants compared with classical case-control methods. As an initiative from the ENIGMA Analytical Working Group, we provide user-friendly scripts and preformatted Excel calculators for implementation of the method for rare variants in BRCA1, BRCA2, and other high-risk genes with known penetrance.

Abstract Background: DEPendency of association on the number of Top Hits (DEPTH) is an approach to identify candidate susceptibility regions by considering the risk signals from overlapping groups of sequential variants across the genome. Methods: We applied a DEPTH analysis using a sliding window of 200 SNPs to colorectal cancer data from the Colon Cancer Family Registry (CCFR; 5,735 cases and 3,688 controls), and Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO; 8,865 cases and 10,285 controls) studies. A DEPTH score > 1 was used to identify candidate susceptibility regions common to both analyses. We compared DEPTH results against those from conventional genome-wide association study (GWAS) analyses of these two studies as well as against 132 published susceptibility regions. Results: Initial DEPTH analysis revealed 2,622 (CCFR) and 3,686 (GECCO) candidate susceptibility regions, of which 569 were common to both studies. Bootstrapping revealed 40 and 49 candidate susceptibility regions in the CCFR and GECCO data sets, respectively. Notably, DEPTH identified at least 82 regions that would not be detected using conventional GWAS methods, nor had they been identified by previous colorectal cancer GWASs. We found four reproducible candidate susceptibility regions (2q22.2, 2q33.1, 6p21.32, 13q14.3). The highest DEPTH scores were in the human leukocyte antigen locus at 6p21 where the strongest associated SNPs were rs762216297, rs149490268, rs114741460, and rs199707618 for the CCFR data, and rs9270761 for the GECCO data. Conclusions: DEPTH can identify candidate susceptibility regions for colorectal cancer not identified using conventional analyses of larger datasets. Impact: DEPTH has potential as a powerful complementary tool to conventional GWAS analyses for discovering susceptibility regions within the genome.

Abstract Breast cancer (BC) patients with a germline CHEK2 c.1100delC variant have an increased risk of contralateral BC (CBC) and worse BC-specific survival (BCSS) compared to non-carriers. We aimed to assess the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS. Analyses were based on 82,701 women diagnosed with invasive BC including 963 CHEK2 c.1100delC carriers; median follow-up was 9.1 years. Differential associations of treatment by CHEK2 c.1100delC status were tested by including interaction terms in a multivariable Cox regression model. A multi-state model was used for further insight into the relation between CHEK2 c.1100delC status, treatment, CBC risk and death. There was no evidence for differential associations of therapy with CBC risk by CHEK2 c.1100delC status The strongest association with reduced CBC risk was observed for the combination of chemotherapy and endocrine therapy [HR(95%CI): 0.66 (0.55–0.78)]. No association was observed with radiotherapy. Results from the multi-state model showed shorter BCSS for CHEK2 c.1100delC carriers versus non-carriers also after accounting for CBC occurrence [HR(95%CI) :1.30 (1.09–1.56)]. In conclusion, systemic therapy was associated with reduced CBC risk irrespective of CHEK2 c.1100delC status. Moreover, CHEK2 c.1100delC carriers had shorter BCSS, which appears not to be fully explained by their CBC risk. (Main MS: 3201 words)

AbstractBackgroundReproductive factors have been shown to be differentially associated with risk of estrogen receptor (ER)-positive and ER-negative breast cancer. However, their associations with intrinsic-like subtypes are less clear.MethodsAnalyses included up to 23 353 cases and 71 072 controls pooled from 31 population-based case-control or cohort studies in the Breast Cancer Association Consortium across 16 countries on 4 continents. Polytomous logistic regression was used to estimate the association between reproductive factors and risk of breast cancer by intrinsic-like subtypes (luminal A-like, luminal B-like, luminal B-HER2–like, HER2-enriched–like, and triple-negative breast cancer) and by invasiveness. All statistical tests were 2-sided.ResultsCompared with nulliparous women, parous women had a lower risk of luminal A-like, luminal B-like, luminal B-HER2–like, and HER2-enriched–like disease. This association was apparent only after approximately 10 years since last birth and became stronger with increasing time (odds ratio [OR] = 0.59, 95% confidence interval [CI] = 0.49 to 0.71; and OR = 0.36, 95% CI = 0.28 to 0.46 for multiparous women with luminal A-like tumors 20 to less than 25 years after last birth and 45 to less than 50 years after last birth, respectively). In contrast, parous women had a higher risk of triple-negative breast cancer right after their last birth (for multiparous women: OR = 3.12, 95% CI = 2.02 to 4.83) that was attenuated with time but persisted for decades (OR = 1.03, 95% CI = 0.79 to 1.34, for multiparous women 25 to less than 30 years after last birth). Older age at first birth (Pheterogeneity < .001 for triple-negative compared with luminal A-like breast cancer) and breastfeeding (Pheterogeneity < .001 for triple-negative compared with luminal A-like breast cancer) were associated with lower risk of triple-negative breast cancer but not with other disease subtypes. Younger age at menarche was associated with higher risk of all subtypes; older age at menopause was associated with higher risk of luminal A-like but not triple-negative breast cancer. Associations for in situ tumors were similar to luminal A-like.ConclusionsThis large and comprehensive study demonstrates a distinct reproductive risk factor profile for triple-negative breast cancer compared with other subtypes, with implications for the understanding of disease etiology and risk prediction.

Salicylic acid (SA) has observationally been shown to decrease colorectal cancer (CRC) risk. Aspirin (acetylsalicylic acid, that rapidly deacetylates to SA) is an effective primary and secondary chemopreventive agent. Through a Mendelian randomization (MR) approach, we aimed to address whether levels of SA affected CRC risk, stratifying by aspirin use. A two-sample MR analysis was performed using GWAS summary statistics of SA (INTERVAL and EPIC-Norfolk, N = 14,149) and CRC (CCFR, CORECT, GECCO and UK Biobank, 55,168 cases and 65,160 controls). The DACHS study (4410 cases and 3441 controls) was used for replication and stratification of aspirin-use. SNPs proxying SA were selected via three methods: (1) functional SNPs that influence the activity of aspirin-metabolising enzymes; (2) pathway SNPs present in enzymes’ coding regions; and (3) genome-wide significant SNPs. We found no association between functional SNPs and SA levels. The pathway and genome-wide SNPs showed no association between SA and CRC risk (OR: 1.03, 95% CI: 0.84–1.27 and OR: 1.08, 95% CI: 0.86–1.34, respectively). Results remained unchanged upon aspirin use stratification. We found little evidence to suggest that an SD increase in genetically predicted SA protects against CRC risk in the general population and upon stratification by aspirin use.

Abstract Background: Fusobacterium nucleatum (F. nucleatum) activates oncogenic signaling pathways and induces inflammation to promote colorectal carcinogenesis. Methods: We characterized F. nucleatum and its subspecies in colorectal tumors and examined associations with tumor characteristics and colorectal cancer–specific survival. We conducted deep sequencing of nusA, nusG, and bacterial 16s rRNA genes in tumors from 1,994 patients with colorectal cancer and assessed associations between F. nucleatum presence and clinical characteristics, colorectal cancer–specific mortality, and somatic mutations. Results: F. nucleatum, which was present in 10.3% of tumors, was detected in a higher proportion of right-sided and advanced-stage tumors, particularly subspecies animalis. Presence of F. nucleatum was associated with higher colorectal cancer–specific mortality (HR, 1.97; P = 0.0004). This association was restricted to nonhypermutated, microsatellite-stable tumors (HR, 2.13; P = 0.0002) and those who received chemotherapy [HR, 1.92; confidence interval (CI), 1.07–3.45; P = 0.029). Only F. nucleatum subspecies animalis, the main subspecies detected (65.8%), was associated with colorectal cancer–specific mortality (HR, 2.16; P = 0.0016), subspecies vincentii and nucleatum were not (HR, 1.07; P = 0.86). Additional adjustment for tumor stage suggests that the effect of F. nucleatum on mortality is partly driven by a stage shift. Presence of F. nucleatum was associated with microsatellite instable tumors, tumors with POLE exonuclease domain mutations, and ERBB3 mutations, and suggestively associated with TP53 mutations. Conclusions: F. nucleatum, and particularly subspecies animalis, was associated with a higher colorectal cancer–specific mortality and specific somatic mutated genes. Impact: Our findings identify the F. nucleatum subspecies animalis as negatively impacting colorectal cancer mortality, which may occur through a stage shift and its effect on chemoresistance.

Abstract Background Individuals who carry pathogenic mutations in DNA mismatch repair (MMR) genes have high risks of cancer, and small studies have suggested that these risks depend on the sex of the parent from whom the mutation was inherited. We have conducted the first large study of such a parent-of-origin effect (POE). Methods Our study was based on all MMR gene mutation carriers and their relatives in the Colon Cancer Family Registry, comprising 18,226 people. The POE was estimated as a hazard ratio (HR) using a segregation analysis approach that adjusted for ascertainment. HR = 1 corresponds to no POE and HR > 1 corresponds to higher risks for maternal mutations. Results For all MMR genes combined, the estimated POE HRs were 1.02 (95% confidence interval (CI) 0.75-1.39, p = 0.9) for male colorectal cancer, 1.12 (95% CI 0.81-1.54, p = 0.5) for female colorectal cancer and 0.84 (95% CI 0.52-1.36, p = 0.5) for endometrial cancer. Separate results for each MMR gene were similar. Conclusions Despite being well-powered, our study did not find any evidence that cancer risks for MMR gene mutation carriers depend on the parent-of-origin of the mutation. Based on current evidence, we don’t recommend that POEs be incorporated into the clinical guidelines or advice for such carriers. Key messages MMR gene mutations inherited from the maternal and paternal side confer similar risks of developing colorectal and endometrial cancer.

Background. Lynch syndrome is the most common genetic predisposition for hereditary cancer. Carriers of pathogenic changes in mismatch repair (MMR) genes have an increased risk of developing colorectal (CRC), endometrial, ovarian, urinary tract, prostate, and other cancers, depending on which gene is malfunctioning. In Lynch syndrome, differences in cancer incidence (penetrance) according to the gene involved have led to the stratification of cancer surveillance. By contrast, any differences in penetrance determined by the type of pathogenic variant remain unknown. Objective. To determine cumulative incidences of cancer in carriers of truncating and missense or aberrant splicing pathogenic variants of the MLH1 and MSH2 genes. Methods. Carriers of pathogenic variants of MLH1 (path_MLH1) and MSH2 (path_MSH2) genes filed in the Prospective Lynch Syndrome Database (PLSD) were categorized as truncating or missense/aberrant splicing according to the InSiGHT criteria for pathogenicity. Results. Among 5199 carriers, 1045 had missense or aberrant splicing variants, and 3930 had truncating variants. Prospective observation years for the two groups were 8205 and 34,141 years, respectively, after which there were no significant differences in incidences for cancer overall or for colorectal cancer or endometrial cancers separately. Conclusion. Truncating and missense or aberrant splicing pathogenic variants were associated with similar average cumulative incidences of cancer in carriers of path MLH1 and path_MSH2.