Search

Background: This study aimed to investigate clinicopathological and molecular tumour features associated with intratumoral pks+ Escherichia coli (pks+E.coli+), pks+E.coli- (non-E.coli bacteria harbouring the pks island), Enterotoxigenic Bacteroides fragilis (ETBF) and Fusobacterium nucleatum (F. nucleatum). Methods: We screened 1697 tumour-derived DNA samples from the Australasian Colorectal Cancer Family Registry, Melbourne Collaborative Cohort Study and the ANGELS study using targeted PCR. Results: Pks+E.coli+ was associated with male sex (P < 0.01) and APC:c.835-8 A > G somatic mutation (P = 0.03). The association between pks+E.coli+ and APC:c.835-8 A > G was specific to early-onset CRCs (diagnosed<45years, P = 0.02). The APC:c.835-A > G was not associated with pks+E.coli- (P = 0.36). F. nucleatum was associated with DNA mismatch repair deficiency (MMRd), BRAF:c.1799T>A p.V600E mutation, CpG island methylator phenotype, proximal tumour location, and high levels of tumour infiltrating lymphocytes (Ps < 0.01). In the stratified analysis by MMRd subgroups, F. nucleatum was associated with Lynch syndrome, MLH1 methylated and double MMR somatic mutated MMRd subgroups (Ps < 0.01). Conclusion: Intratumoral pks+E.coli+ but not pks+E.coli- are associated with CRCs harbouring the APC:c.835-8 A > G somatic mutation, suggesting that this mutation is specifically related to DNA damage from colibactin-producing E.coli exposures. F. nucleatum was associated with both hereditary and sporadic MMRd subtypes, suggesting the MMRd tumour microenvironment is important for F. nucleatum colonisation irrespective of its cause.

BACKGROUND: Obesity and internalising disorders, including depression and anxiety, often co-occur. There is evidence that familial confounding contributes to the co-occurrence of internalising disorders and obesity in adults. However, its impact on this association among young people is unclear. Our study investigated the extent to which familial factors confound the association between internalising disorders and obesity in adolescents and young adults. SUBJECTS/METHODS: We used a matched co-twin design to investigate the impact of confounding by familial factors on associations between internalising symptoms and obesity in a sample of 4018 twins aged 16 to 27 years. RESULTS: High levels of internalising symptoms compared to low levels increased the odds of obesity for the whole cohort (adjusted odds ratio [AOR] = 3.1, 95% confidence interval [CI]: 1.5, 6.8), and in females (AOR = 4.1, 95% CI 1.5, 11.1), but not in males (AOR = 2.8 95% CI 0.8, 10.0). We found evidence that internalising symptoms were associated with an increased between-pair odds of obesity (AOR 6.2, 95% CI 1.7, 22.8), using the paired analysis but not using a within-pair association, which controls for familial confounding. Sex-stratified analyses indicated high internalising symptoms were associated with increased between-pair odds of obesity for females (AOR 12.9, 95% CI 2.2, 76.8), but this attenuated to the null using within-pair analysis. We found no evidence of between or within-pair associations for males and weak evidence that sex modified the association between internalising symptoms and obesity (likelihood ratio test p = 0.051). CONCLUSIONS: Some familial factors shared by twins confound the association between internalising symptoms and obesity in adolescent and young adult females. Internalising symptoms and obesity were not associated for adolescent and young adult males. Therefore, prevention and treatment efforts should especially address familial shared determinants of obesity, pa

A polygenic risk score (PRS) combines the associations of multiple genetic variants that could be due to direct causal effects, indirect genetic effects, or other sources of familial confounding. We have developed new approaches to assess evidence for and against causation by using family data for pairs of relatives (Inference about Causation from Examination of FAmiliaL CONfounding [ICE FALCON]) or measures of family history (Inference about Causation from Examining Changes in Regression coefficients and Innovative STatistical AnaLyses [ICE CRISTAL]). Inference is made from the changes in regression coefficients of relatives' PRSs or PRS and family history before and after adjusting for each other. We applied these approaches to two breast cancer PRSs and multiple studies and found that (a) for breast cancer diagnosed at a young age, for example, <50 years, there was no evidence that the PRSs were causal, while (b) for breast cancer diagnosed at later ages, there was consistent evidence for causation explaining increasing amounts of the PRS-disease association. The genetic variants in the PRS might be in linkage disequilibrium with truly causal variants and not causal themselves. These PRSs cause minimal heritability of breast cancer at younger ages. There is also evidence for nongenetic factors shared by first-degree relatives that explain breast cancer familial aggregation. Familial associations are not necessarily due to genes, and genetic associations are not necessarily causal.

Young breast and bowel cancers (e.g., those diagnosed before age 40 or 50 years) have far greater morbidity and mortality in terms of years of life lost, and are increasing in incidence, but have been less studied. For breast and bowel cancers, the familial relative risks, and therefore the familial variances in age-specific log(incidence), are much greater at younger ages, but little of these familial variances has been explained. Studies of families and twins can address questions not easily answered by studies of unrelated individuals alone. We describe existing and emerging family and twin data that can provide special opportunities for discovery. We present designs and statistical analyses, including novel ideas such as the VALID (Variance in Age-specific Log Incidence Decomposition) model for causes of variation in risk, the DEPTH (DEPendency of association on the number of Top Hits) and other approaches to analyse genome-wide association study data, and the within-pair, ICE FALCON (Inference about Causation from Examining FAmiliaL CONfounding) and ICE CRISTAL (Inference about Causation from Examining Changes in Regression coefficients and Innovative STatistical AnaLysis) approaches to causation and familial confounding. Example applications to breast and colorectal cancer are presented. Motivated by the availability of the resources of the Breast and Colon Cancer Family Registries, we also present some ideas for future studies that could be applied to, and compared with, cancers diagnosed at older ages and address the challenges posed by young breast and bowel cancers.

BACKGROUND: We examined associations between adherence to adaptations of the 2018 World Cancer Research Fund/American Institute for Cancer Research (WCRF/AICR) cancer prevention recommendations and total, exposure-related and site-specific cancer risk. METHODS: A total of 20,001 participants ages 40 to 69 years at enrollment into the Melbourne Collaborative Cohort Study in 1990 to 1994, who had diet, body size, and lifestyle reassessed in 2003 to 2007 ("baseline"), were followed-up through June 2021. We constructed diet and standardized lifestyle scores based on core WCRF/AICR recommendations on diet, alcohol intake, body size and physical activity, and additional scores incorporating weight change, sedentary behavior, and smoking. Associations with cancer risk were estimated using Cox regression, adjusting for confounders. RESULTS: During follow-up (mean = 16 years), 4,710 incident cancers were diagnosed. For highest quintile ("most adherent") of the standardized lifestyle score, compared with lowest ("least adherent"), a HR of 0.82 [95% confidence interval (CI): 0.74-0.92] was observed for total cancer. This association was stronger with smoking included in the score (HR = 0.74; 95% CI: 0.67-0.81). A higher score was associated with lower breast and prostate cancer risk for the standardized score, and with lung, stomach, rectal, and pancreatic cancer risk when the score included smoking. Our analyses identified alcohol use, waist circumference and smoking as key drivers of associations with total cancer risk. CONCLUSIONS: Adherence to WCRF/AICR cancer prevention recommendations is associated with lower cancer risk. IMPACT: With <0.2% of our sample fully adherent to the recommendations, the study emphasizes the vast potential for preventing cancer through modulation of lifestyle habits.

BACKGROUND: Mental disorders and perceived discrimination share common risk factors. The association between having a mental disorder and experiencing discrimination is well-known, but the extent to which familial factors, such as genetic and shared environmental factors, might confound this association, including sex differences in familial confounding, remains unexplored. AIMS: We investigated potential unmeasured familial confounding in the association between mental disorders and perceived discrimination using a matched twin study design. METHOD: We examined data from 2044 same-sex twin pairs (n = 4088) aged 16-25 years from the German population-based study 'TwinLife'. We applied random-effects logistic regression to within-individual and within-and-between pair models of the association between mental disorder and perceived discrimination, and used likelihood ratio tests (LRTs) to compare these models. Multivariable models were adjusted for body mass index, educational attainment, and life satisfaction. RESULTS: There were 322 (8.1%) participants with a diagnosed mental disorder, and 15% (n = 604) of the cohort reported having experienced discrimination in the previous 12 months. Mental disorder and discrimination were associated in the adjusted within-individual model (adjusted odds ratio = 2.19, 95% confidence interval: 1.42-3.39, P<0.001). However, the within-and-between pair model showed that this association was explained by the within-pair mean (aOR = 4.24, 95% CI: 2.17-8.29, P<0.001) and not the within-pair difference (aOR = 1.26, 95% CI: 0.70-2.28, P = 0.4) of mental disorder. Therefore, this association was mostly explained by familial confounding, which is also supported by the LRTs for the unadjusted and adjusted models (P<0.001 and P = 0.03, respectively). This familial confounding was more prominent for males than females. CONCLUSIONS: Our findings show that the association between mental disorder and discrimination is at least partially explained

Pubertal timing varies considerably and is associated with later health outcomes. We performed multi-ancestry genetic analyses on ~800,000 women, identifying 1,080 signals for age at menarche. Collectively, these explained 11% of trait variance in an independent sample. Women at the top and bottom 1% of polygenic risk exhibited ~11 and ~14-fold higher risks of delayed and precocious puberty, respectively. We identified several genes harboring rare loss-of-function variants in ~200,000 women, including variants in ZNF483, which abolished the impact of polygenic risk. Variant-to-gene mapping approaches and mouse gonadotropin-releasing hormone neuron RNA sequencing implicated 665 genes, including an uncharacterized G-protein-coupled receptor, GPR83, which amplified the signaling of MC3R, a key nutritional sensor. Shared signals with menopause timing at genes involved in DNA damage response suggest that the ovarian reserve might signal centrally to trigger puberty. We also highlight body size-dependent and independent mechanisms that potentially link reproductive timing to later life disease.

BACKGROUND: Cirrus is an automated risk predictor for breast cancer that comprises texture-based mammographic features and is mostly independent of mammographic density. We investigated genetic and environmental variance of variation in Cirrus. METHODS: We measured Cirrus for 3,195 breast cancer-free participants, including 527 pairs of monozygotic (MZ) twins, 271 pairs of dizygotic (DZ) twins, and 1,599 siblings of twins. Multivariate normal models were used to estimate the variance and familial correlations of age-adjusted Cirrus as a function of age. The classic twin model was expanded to allow the shared environment effects to differ by zygosity. The SNP-based heritability was estimated for a subset of 2,356 participants. RESULTS: There was no evidence that the variance or familial correlations depended on age. The familial correlations were 0.52 (SE, 0.03) for MZ pairs and 0.16(SE, 0.03) for DZ and non-twin sister pairs combined. Shared environmental factors specific to MZ pairs accounted for 20% of the variance. Additive genetic factors accounted for 32% (SE = 5%) of the variance, consistent with the SNP-based heritability of 36% (SE = 16%). CONCLUSION: Cirrus is substantially familial due to genetic factors and an influence of shared environmental factors that was evident for MZ twin pairs only. The latter could be due to nongenetic factors operating in utero or in early life that are shared by MZ twins. IMPACT: Early-life factors, shared more by MZ pairs than DZ/non-twin sister pairs, could play a role in the variation in Cirrus, consistent with early life being recognized as a critical window of vulnerability to breast carcinogens.

BACKGROUND: 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. AIM: To assessed the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS. METHODS: Analyses were based on 82,701 women diagnosed with a first primary invasive BC including 963 CHEK2 c.1100delC carriers; median follow-up was 9.1 years. Differential associations with 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. RESULTS: 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)]. 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.

BACKGROUND: Mammogram risk scores based on texture and density defined by different brightness thresholds are associated with breast cancer risk differently and could reveal distinct information about breast cancer risk. We aimed to investigate causal relationships between these intercorrelated mammogram risk scores to determine their relevance to breast cancer aetiology. METHODS: We used digitised mammograms for 371 monozygotic twin pairs, aged 40-70 years without a prior diagnosis of breast cancer at the time of mammography, from the Australian Mammographic Density Twins and Sisters Study. We generated normalised, age-adjusted, and standardised risk scores based on textures using the Cirrus algorithm and on three spatially independent dense areas defined by increasing brightness threshold: light areas, bright areas, and brightest areas. Causal inference was made using the Inference about Causation from Examination of FAmilial CONfounding (ICE FALCON) method. RESULTS: The mammogram risk scores were correlated within twin pairs and with each other (r = 0.22-0.81; all P < 0.005). We estimated that 28-92% of the associations between the risk scores could be attributed to causal relationships between the scores, with the rest attributed to familial confounders shared by the scores. There was consistent evidence for positive causal effects: of Cirrus, light areas, and bright areas on the brightest areas (accounting for 34%, 55%, and 85% of the associations, respectively); and of light areas and bright areas on Cirrus (accounting for 37% and 28%, respectively). CONCLUSIONS: In a mammogram, the lighter (less dense) areas have a causal effect on the brightest (highly dense) areas, including through a causal pathway via textural features. These causal relationships help us gain insight into the relative aetiological importance of different mammographic features in breast cancer. For example our findings are consistent with the brightest areas being more aetiologically impor

BACKGROUND: Low back pain (LBP) is more likely to occur in people with a family history of this condition, highlighting the importance of accounting for familial factors when studying the individual risk of LBP. We conducted a study of opposite-sex twin pairs investigating sex differences in LBP while accounting for (genetic and shared environmental) familial factors. METHODS: We applied a matched co-twin control design to study 795 adult opposite-sex pairs from Australia, Spain, and the United States (US). We used mixed-effects logistic regression to assess the within-pair association between female sex and lifetime prevalence of LBP in unadjusted and adjusted models for body-mass-index, and depression, as well as interactions between female sex and age (

Enhanced survival following a diagnosis of cancer has led to a steep rise in the number of individuals diagnosed with a second primary cancer. We examined the association between pre-cancer cigarette smoking and risk of second cancer in 9785 participants diagnosed with first invasive cancer after enrolment in the Melbourne Collaborative Cohort Study. Follow-up was from date of first invasive cancer until diagnosis of second primary invasive cancer, death, or 31 July 2019, whichever came first. Data on cigarette smoking was collected at enrolment (1990-94) along with information on other lifestyle factors including body size, alcohol intake and diet. We estimated hazard ratios (HR) and 95 % confidence intervals (CI) for incident second cancer with several smoking measures, adjusted for potential confounders. After a mean follow-up of 7.3 years, 1658 second cancers were identified. All measures of smoking were associated with increased risk of second cancer. We observed a 44 % higher risk of second cancer for smokers of ≥ 20 cigarettes/day (HR=1.44, 95 % CI: 1.18-1.76), compared with never smokers. We also observed dose-dependent associations with number of cigarettes smoked (HR=1.05 per 10 cigarettes/day, 95 % CI: 1.01-1.09) and duration of smoking (HR=1.07 per 10 years, 95 % CI: 1.03-1.10). The risk of second cancer increased by 4 % per 10 pack-years of smoking (HR=1.04, 95 % CI: 1.02-1.06; p < 0.001). There was suggestive evidence of stronger associations with number of cigarettes smoked and pack-years of smoking for women (pinteraction<0.05), particularly for the highest risk categories of both variables. These associations with pre-diagnostic smoking were markedly stronger for second cancers known to be smoking-related than for others (phomogeneity<0.001). Our findings for pre-diagnostic cigarette smoking indicated increased risk of second primary cancer for cancer sites considered smoking-related, highlighting the importance of assessing smoking habits in cancer

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 analyses 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 to classical case-control methods. As an initiative from the ENIGMA Analytical Working Group, we provide user-friendly scripts and pre-formatted excel calculators for implementation of the method for rare variants in BRCA1, BRCA2 and other high-risk genes with known penetrance.

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.

Routine screening of tumors for DNA mismatch repair (MMR) deficiency (dMMR) in colorectal (CRC), endometrial (EC) and sebaceous skin (SST) tumors leads to a significant proportion of unresolved cases classified as suspected Lynch syndrome (SLS). SLS cases (n=135) were recruited from Family Cancer Clinics across Australia and New Zealand. Targeted panel sequencing was performed on tumor (n=137; 80xCRCs, 33xECs and 24xSSTs) and matched blood-derived DNA to assess for microsatellite instability status, tumor mutation burden, COSMIC tumor mutational signatures and to identify germline and somatic MMR gene variants. MMR immunohistochemistry (IHC) and MLH1 promoter methylation were repeated. In total, 86.9% of the 137 SLS tumors could be resolved into established subtypes. For 22.6% of these resolved SLS cases, primary MLH1 epimutations (2.2%) as well as previously undetected germline MMR pathogenic variants (1.5%), tumor MLH1 methylation (13.1%) or false positive dMMR IHC (5.8%) results were identified. Double somatic MMR gene mutations were the major cause of dMMR identified across each tumor type (73.9% of resolved cases, 64.2% overall, 70% of CRC, 45.5% of ECs and 70.8% of SSTs). The unresolved SLS tumors (13.1%) comprised tumors with only a single somatic (7.3%) or no somatic (5.8%) MMR gene mutations. A tumor-focused testing approach reclassified 86.9% of SLS into Lynch syndrome, sporadic dMMR or MMR-proficient cases. These findings support the incorporation of tumor sequencing and alternate MLH1 methylation assays into clinical diagnostics to reduce the number of SLS patients and provide more appropriate surveillance and screening recommendations.

Light-at-night triggers the decline of pineal gland melatonin biosynthesis and secretion and is an IARC-classified probable breast-cancer risk factor. We applied a large-scale molecular epidemiology approach to shed light on the putative role of melatonin in breast cancer. We investigated associations between breast-cancer risk and polymorphisms at genes of melatonin biosynthesis/signaling using a study population of 44,405 women from the Breast Cancer Association Consortium (22,992 cases, 21,413 population-based controls). Genotype data of 97 candidate single nucleotide polymorphisms (SNPs) at 18 defined gene regions were investigated for breast-cancer risk effects. We calculated adjusted odds ratios (ORs) and 95% confidence intervals (CI) by logistic regression for the main-effect analysis as well as stratified analyses by estrogen- and progesterone-receptor (ER, PR) status. SNP-SNP interactions were analyzed via a two-step procedure based on logic regression. The Bayesian false-discovery probability (BFDP) was used for all analyses to account for multiple testing. Noteworthy associations (BFDP < 0.8) included 10 linked SNPs in tryptophan hydroxylase 2 (TPH2) (e.g. rs1386492: OR = 1.07, 95% CI 1.02-1.12), and a SNP in the mitogen-activated protein kinase 8 (MAPK8) (rs10857561: OR = 1.11, 95% CI 1.04-1.18). The SNP-SNP interaction analysis revealed noteworthy interaction terms with TPH2- and MAPK-related SNPs (e.g. rs1386483R ∧ rs1473473D ∧ rs3729931D: OR = 1.20, 95% CI 1.09-1.32). In line with the light-at-night hypothesis that links shift work with elevated breast-cancer risks our results point to SNPs in TPH2 and MAPK-genes that may impact the intricate network of circadian regulation.

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.

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.

BACKGROUND: MLH1 epimutation is characterised by constitutional monoallelic MLH1 promoter hypermethylation, which can cause colorectal cancer (CRC). Tumour molecular profiles of MLH1 epimutation CRCs were used to classify germline MLH1 promoter variants of uncertain significance and MLH1 methylated early-onset CRCs (EOCRCs). Genome-wide DNA methylation and somatic mutational profiles of tumours from two germline MLH1: c.-11C > T and one MLH1: c.-[28A > G; 7C > T] carriers and three MLH1 methylated EOCRCs (< 45 years) were compared with 38 reference CRCs. Methylation-sensitive droplet digital PCR (ddPCR) was used to detect mosaic MLH1 methylation in blood, normal mucosa and buccal DNA. RESULTS: Genome-wide methylation-based Consensus Clustering identified four clusters where the tumour methylation profiles of germline MLH1: c.-11C > T carriers and MLH1 methylated EOCRCs clustered with the constitutional MLH1 epimutation CRCs but not with the sporadic MLH1 methylated CRCs. Furthermore, monoallelic MLH1 methylation and APC promoter hypermethylation in tumour were observed in both MLH1 epimutation and germline MLH1: c.-11C > T carriers and MLH1 methylated EOCRCs. Mosaic constitutional MLH1 methylation in MLH1: c.-11C > T carriers and 1 of 3 MLH1 methylated EOCRCs was identified by methylation-sensitive ddPCR. CONCLUSIONS: Mosaic MLH1 epimutation underlies the CRC aetiology in MLH1: c.-11C > T germline carriers and a subset of MLH1 methylated EOCRCs. Tumour profiling and ultra-sensitive ddPCR methylation testing can be used to identify mosaic MLH1 epimutation carriers.

Identifying tumor DNA mismatch repair deficiency (dMMR) is important for precision medicine. Tumor features, individually and in combination, derived from whole-exome sequenced (WES) colorectal cancers (CRCs) and panel-sequenced CRCs, endometrial cancers (ECs), and sebaceous skin tumors (SSTs) were assessed for their accuracy in detecting dMMR. CRCs (n = 300) with WES, where mismatch repair status was determined by immunohistochemistry, were assessed for microsatellite instability (MSMuTect, MANTIS, MSIseq, and MSISensor), Catalogue of Somatic Mutations in Cancer tumor mutational signatures, and somatic mutation counts. A 10-fold cross-validation approach (100 repeats) evaluated the dMMR prediction accuracy for i) individual features, ii) Lasso statistical model, and iii) an additive feature combination approach. Panel-sequenced tumors (29 CRCs, 22 ECs, and 20 SSTs) were assessed for the top performing dMMR predicting features/models using these three approaches. For WES CRCs, 10 features provided >80% dMMR prediction accuracy, with MSMuTect, MSIseq, and MANTIS achieving ≥99% accuracy. The Lasso model achieved 98.3% accuracy. The additive feature approach, with three or more of six of MSMuTect, MANTIS, MSIseq, MSISensor, insertion-deletion count, or tumor mutational signature small insertion/deletion 2 + small insertion/deletion 7 achieved 99.7% accuracy. For the panel-sequenced tumors, the additive feature combination approach of three or more of six achieved accuracies of 100%, 95.5%, and 100% for CRCs, ECs, and SSTs, respectively. The microsatellite instability calling tools performed well in WES CRCs; however, an approach combining tumor features may improve dMMR prediction in both WES and panel-sequenced data across tissue types.

BACKGROUND: Wildfire-related fine particulate matter (PM2.5) has many adverse health impacts, but its impacts on human epigenome are unknown. We aimed to evaluate the associations between long-term exposure to wildfire-related PM2.5 and blood DNA methylation, and whether the associations differ from those with non-wildfire-related PM2.5. METHODS: We studied 479 Australian women comprising 132 twin pairs and 215 of their sisters. Blood-derived DNA methylation was measured using the HumanMethylation450 BeadChip array. Data on 3-year (year of blood collection and previous two years) average wildfire-related and non-wildfire-related PM2.5 at 0.01°×0.01° spatial resolution were created by combining information from satellite observations, chemical transport models, and ground-based observations. Exposure data were linked to each participant's home address, assuming the address did not change during the exposure window. For DNA methylation of each cytosine-guanine dinucleotide (CpG), and for global DNA methylation represented by the average of all measured CpGs or CpGs in repetitive elements, we evaluated their associations with wildfire- or non-wildfire-related PM2.5 using a within-sibship analysis controlling for factors shared between siblings and other important covariates. Differentially methylated regions (DMRs) were defined by comb-p and DMRcate. RESULTS: The 3-year average wildfire-related PM2.5 (range: 0.3 to 7.6 µg/m3, mean: 1.6 µg/m3) was negatively, but not significantly (p-values greater than 0.05) associated with all seven global DNA methylation measures. There were 26 CpGs and 33 DMRs associated with wildfire-related PM2.5 (Bonferroni adjusted p-value < 0.05) mapped to 47 genes enriched for pathways related to inflammatory regulation and platelet activation. These genes have been related to many human diseases or phenotypes e.g., cancer, mental disorders, diabetes, obesity, asthma, blood pressure. These CpGs, DMRs and enriched pathways did not overlap with

BACKGROUND: Temperature fluctuations can affect human health independent of the effect of mean temperature. However, no study has evaluated whether short-term temperature fluctuations could affect DNA methylation. METHODS: Peripheral blood DNA methylation for 479 female siblings of 130 families were analysed. Gridded daily temperatures data were obtained, linked to each participant's home address, and used to calculate nine different metrics of short-term temperature fluctuations: temperature variabilities (TVs) within the day of blood draw and preceding one to seven days (TV 0-1 to TV 0-7), diurnal temperature range (DTR), and temperature change between neighbouring days (TCN). Within-sibship design was used to perform epigenome-wide association analyses, adjusting for daily mean temperatures, and other important covariates (e.g., smoking, alcohol use, cell-type proportions). Differentially methylated regions (DMRs) were further identified. Multiple-testing comparisons with a significant threshold of 0.01 for cytosine-guanine dinucleotides (CpGs) and 0.05 for DMRs were applied. RESULTS: Among 479 participants (mean age ± SD, 56.4 ± 7.9 years), we identified significant changes in methylation levels in 14 CpGs and 70 DMRs associated with temperature fluctuations. Almost all identified CpGs were associated with exposure to temperature fluctuations within three days. Differentially methylated signals were mapped to 68 genes that were linked to human diseases such as cancer (e.g., colorectal carcinoma, breast carcinoma, and metastatic neoplasms) and mental disorder (e.g., schizophrenia, mental depression, and bipolar disorder). The top three most significantly enriched gene ontology terms were Response to bacterium (TV 0-3), followed by Hydrolase activity, acting on ester bonds (TCN), and Oxidoreductase activity (TV 0-3). CONCLUSIONS: Short-term temperature fluctuations were associated with differentially methylated signals across the human genome, which provides evide

Background: Previous findings for the genetic and environmental contributions to DNA methylation variation were for limited age ranges only. We investigated the lifespan contributions and their implications for human health for the first time. Methods: 1,720 monozygotic twin (MZ) pairs and 1,107 dizygotic twin (DZ) pairs aged 0-92 years were included. Familial correlations (i.e., correlations between twins) for 353,681 methylation sites were estimated and modelled as a function of twin pair cohabitation history. Findings: The methylome average familial correlation was around zero at birth (MZ pair: -0.01; DZ pair: -0.04), increased with the time of twins living together during childhood at rates of 0.16 (95%CI: 0.12-0.20) for MZ pairs and 0.13 (95%CI: 0.07-0.20) for DZ pairs per decade, and decreased with the time of living apart during adulthood at rates of 0.026 (95%CI: 0.019-0.033) for MZ pairs and 0.027 (95%CI: 0.011-0.043) for DZ pairs per decade. Neither the increasing nor decreasing rate differed by zygosity (both P>0.1), consistent with cohabitation environment shared by twins, rather than genetic factors, influencing the methylation familial correlation changes. Familial correlations for 6.6% (23,386/353,681) sites changed with twin pair cohabitation history. These sites were enriched for high heritability, proximal promoters, and epigenetic/genetic associations with various early-life factors and late-life health conditions. Interpretation: Early life strongly influences DNA methylation variation across the lifespan, and the effects are stronger for heritable sites and sites biologically relevant to the regulation of gene expression. Early life could affect late-life health through influencing DNA methylation. Funding: Victorian Cancer Agency, Cancer Australia, Cure Cancer Foundation.

BACKGROUND: Genome-wide association studies (GWAS) have identified more than 200 susceptibility loci for breast cancer, but these variants explain less than a fifth of the disease risk. Although gene-environment interactions have been proposed to account for some of the remaining heritability, few studies have empirically assessed this. METHODS: We obtained genotype and risk factor data from 46,060 cases and 47,929 controls of European ancestry from population-based studies within the Breast Cancer Association Consortium (BCAC). We built gene expression prediction models for 4,864 genes with a significant (P<0.01) heritable component using the transcriptome and genotype data from the Genotype-Tissue Expression (GTEx) project. We leveraged predicted gene expression information to investigate the interactions between gene-centric genetic variation and 14 established risk factors in association with breast cancer risk, using a mixed-effects score test. RESULTS: After adjusting for number of tests using Bonferroni correction, no interaction remained statistically significant. The strongest interaction observed was between the predicted expression of the C13orf45 gene and age at first full-term pregnancy (PGXE=4.44×10-6). CONCLUSION: In this transcriptome-informed genome-wide gene-environment interaction study of breast cancer, we found no strong support for the role of gene expression in modifying the associations between established risk factors and breast cancer risk. IMPACT: Our study suggests a limited role of gene-environment interactions in breast cancer risk.

BACKGROUND: Genome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER) status, but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear. METHODS: Among 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes. RESULTS: Eighty-five of 173 variants were associated with at least one tumor feature (false discovery rate < 5%), most commonly ER and grade, followed by PR and HER2. Models for intrinsic-like subtypes found nearly all of these variants (83 of 85) associated at p < 0.05 with risk for at least one luminal-like subtype, and approximately half (41 of 85) of the variants were associated with risk of at least one non-luminal subtype, including 32 variants associated with triple-negative (TN) disease. Ten variants were associated with risk of all subtypes in different magnitude. Five variants were associated with risk of luminal A-like and TN subtypes in opposite directions. CONCLUSION: This report demonstrates a high level of complexity in the etiology heterogeneity of breast cancer susceptibility variants and can inform investigations of subtype-specific risk prediction.

BACKGROUND: Mammographic density (MD) phenotypes, including percent density (PMD), area of dense tissue (DA), and area of non-dense tissue (NDA), are associated with breast cancer risk. Twin studies suggest that MD phenotypes are highly heritable. However, only a small proportion of their variance is explained by identified genetic variants. METHODS: We conducted a genome-wide association study, as well as a transcriptome-wide association study (TWAS), of age- and BMI-adjusted DA, NDA, and PMD in up to 27,900 European-ancestry women from the MODE/BCAC consortia. RESULTS: We identified 28 genome-wide significant loci for MD phenotypes, including nine novel signals (5q11.2, 5q14.1, 5q31.1, 5q33.3, 5q35.1, 7p11.2, 8q24.13, 12p11.2, 16q12.2). Further, 45% of all known breast cancer SNPs were associated with at least one MD phenotype at p < 0.05. TWAS further identified two novel genes (SHOX2 and CRISPLD2) whose genetically predicted expression was significantly associated with MD phenotypes. CONCLUSIONS: Our findings provided novel insight into the genetic background of MD phenotypes, and further demonstrated their shared genetic basis with breast cancer.

Estimates from genome-wide association studies (GWAS) of unrelated individuals capture effects of inherited variation (direct effects), demography (population stratification, assortative mating) and relatives (indirect genetic effects). Family-based GWAS designs can control for demographic and indirect genetic effects, but large-scale family datasets have been lacking. We combined data from 178,086 siblings from 19 cohorts to generate population (between-family) and within-sibship (within-family) GWAS estimates for 25 phenotypes. Within-sibship GWAS estimates were smaller than population estimates for height, educational attainment, age at first birth, number of children, cognitive ability, depressive symptoms and smoking. Some differences were observed in downstream SNP heritability, genetic correlations and Mendelian randomization analyses. For example, the within-sibship genetic correlation between educational attainment and body mass index attenuated towards zero. In contrast, analyses of most molecular phenotypes (for example, low-density lipoprotein-cholesterol) were generally consistent. We also found within-sibship evidence of polygenic adaptation on taller height. Here, we illustrate the importance of family-based GWAS data for phenotypes influenced by demographic and indirect genetic effects.

BACKGROUND: Multigene panel tests for breast cancer predisposition routinely include ATM as it is now a well-established breast cancer predisposition gene. METHODS: We included ATM in a multigene panel test applied to the Australian Breast Cancer Family Registry (ABCFR), a population-based case-control-family study of breast cancer, with the purpose of estimating the prevalence and penetrance of heterozygous ATM pathogenic variants from the family data, using segregation analysis. RESULTS: The estimated breast cancer hazard ratio for carriers of pathogenic ATM variants in the ABCFR was 1.32 (95% confidence interval 0.45-3.87; P = 0.6). The estimated cumulative risk of breast cancer to age 80 years for heterozygous ATM pathogenic variant carriers was estimated to be 13% (95% CI 4.6-30). CONCLUSIONS: Although ATM has been definitively identified as a breast cancer predisposition gene, further evidence, such as variant-specific penetrance estimates, are needed to inform risk management strategies for carriers of pathogenic variants to increase the clinical utility of population testing of this gene.

In this issue, Kresovich and colleagues have published a hallmark paper in Molecular, Environmental, Genetic and Analytic Epidemiology. By applying artificial intelligence to the Sister Study they created a new methylation-based breast cancer risk score (mBCRS) based on blood DNA methylation. Using a prospective design and after accounting for age and questionnaire-based breast cancer risk factors, the Odds PER Adjusted standard deviation (OPERA) for mBCRS and polygenic risk score (PRS) was 1.58 (95% CI: 1.38, 1.81) and 1.58 (95% CI: 1.36, 1.83), respectively, and the corresponding area under the receiver operating curve was 0.63 for both. Therefore, mBCRS could be as powerful as the current best PRS in differentiating women of the same age in terms of their breast cancer risk. These risk scores are among the strongest known breast cancer risk-stratifiers, shaded only by new mammogram risk scores based on measures other than conventional mammographic density, such as Cirrocumulus and Cirrus, which when combined have an OPERA as high as 2.3. The combination of PRS and mBCRS with the other measured risk factors gave an OPERA of 2.2. OPERA has many advantages over changes in areas under the receiver operator curve because the latter depend on the order in which risk factors are considered. Although more replication is needed using prospective data to protect against reverse causation, there are many novel molecular and analytic aspects to this paper which uncovers a potential mechanism for how genetic and environmental factors combine to cause breast cancer.

Multiple myeloma (MM) is the second most common hematological cancer and causes significant mortality and morbidity. Knowledge regarding modifiable risk factors for MM remains limited. This analysis of an Australian population‐based case–control family study investigates whether smoking or alcohol consumption is associated with risk of MM and related diseases. Incident cases (n = 789) of MM were recruited via cancer registries in Victoria and New South Wales. Controls (n = 1,113) were either family members of cases (n = 696) or controls recruited for a similarly designed study of renal cancers (n = 417). Adjusted odds ratios (OR) and 95% confidence intervals (CI) were estimated using unconditional multivariable logistic regression. Heavy intake (>20 g ethanol/day) of alcohol had a lower risk of MM compared with nondrinkers (OR = 0.68, 95% CI: 0.50–0.93), and there was an inverse dose–response relationship for average daily alcohol intake (OR per 10 g ethanol per day = 0.92, 95% CI: 0.86–0.99); there was no evidence of an interaction with sex. There was no evidence of an association with MM risk for smoking‐related exposures (p > 0.18). The associations between smoking and alcohol with MM are similar to those with non‐Hodgkin lymphoma. Further research into potential underlying mechanisms is warranted.