Mendelian randomisation study of smoking exposure in relation to breast cancer risk

BACKGROUND

Despite a modest association between tobacco smoking and breast cancer risk reported by recent epidemiological studies, it is still equivocal whether smoking is causally related to breast cancer risk.

METHODS

We applied Mendelian randomisation (MR) to evaluate a potential causal effect of cigarette smoking on breast cancer risk. Both individual-level data as well as summary statistics for 164 single-nucleotide polymorphisms (SNPs) reported in genome-wide association studies of lifetime smoking index (LSI) or cigarette per day (CPD) were used to obtain MR effect estimates. Data from 108,420 invasive breast cancer cases and 87,681 controls were used for the LSI analysis and for the CPD analysis conducted among ever-smokers from 26,147 cancer cases and 26,072 controls. Sensitivity analyses were conducted to address pleiotropy.

RESULTS

Genetically predicted LSI was associated with increased breast cancer risk (OR 1.18 per SD, 95% CI: 1.07-1.30, P = 0.11 × 10-2), but there was no evidence of association for genetically predicted CPD (OR 1.02, 95% CI: 0.78-1.19, P = 0.85). The sensitivity analyses yielded similar results and showed no strong evidence of pleiotropic effect.

CONCLUSION

Our MR study provides supportive evidence for a potential causal association with breast cancer risk for lifetime smoking exposure but not cigarettes per day among smokers.

Inflammation-Related Marker Profiling of Dietary Patterns and All-cause Mortality in the Melbourne Collaborative Cohort Study

BACKGROUND

Nutritional epidemiology research using self-reported dietary intake is prone to measurement error. Objective methods are being explored to overcome this limitation.

OBJECTIVES

We aimed to examine 1) the association between plasma markers related to inflammation and derive marker scores for dietary patterns [Mediterranean dietary score (MDS), energy-adjusted Dietary Inflammatory Index (E-DIITM), Alternative Healthy Eating Index 2010 (AHEI)] and 2) the associations of these marker scores with mortality.

METHODS

Weighted marker scores were derived from the cross-sectional association between 30 plasma markers and each dietary score (assessed using food-frequency questionnaires) using linear regression for 770 participants in the Melbourne Collaborative Cohort Study (aged 50-82 y). Prospective associations between marker scores and mortality (n = 249 deaths) were assessed using Cox regression (median follow-up: 14.4 y).

RESULTS

The MDS, E-DII, and AHEI were associated (P < 0.05) with 9, 14, and 11 plasma markers, respectively. Healthier diets (higher MDS and AHEI, and lower anti-inflammatory, E-DII) were associated with lower concentrations of kynurenines, neopterin, IFN-γ, cytokines, and C-reactive protein. Five of 6 markers common to the 3 dietary scores were components of the kynurenine pathway. The 3 dietary-based marker scores were highly correlated (Spearman ρ: -0.74, -0.82, and 0.93). Inverse associations (for 1-SD increment) were observed with all-cause mortality for the MDS marker score (HR: 0.84; 95% CI: 0.72-0.98) and the AHEI marker score (HR: 0.76; 95% CI: 0.66-0.89), whereas a positive association was observed with the E-DII marker score (HR: 1.18; 95% CI: 1.01-1.39). The same magnitude of effect was not observed for the respective dietary patterns.

CONCLUSIONS

Markers involved in inflammation-related processes are associated with dietary quality, including a substantial overlap between markers associated with the MDS, the E-DII, and the AHEI, especially kynurenines. Unfavorable marker scores, reflecting poorer-quality diets, were associated with increased mortality.

The EUS molecular evaluation of pancreatic cancer: A prospective multicenter cohort trial

Background and Objectives:

Patients with locally advanced or metastatic pancreatic ductal adenocarcinoma (A-PDAC) are not candidates for surgical resection and are often offered palliative chemotherapy. The ready availability of a safe and effective tumor sampling technique to provide material for both diagnosis and comprehensive genetic profiling is critical for informing precision medicine in A-PDAC, thus potentially increasing survival. The aim of this study is to examine the feasibility and benefits of routine comprehensive genomic profiling (CGP) of A-PDAC using EUS-FNA material.

Methods:

This is a prospective cohort study to test the clinical utility of fresh frozen or archival EUS-FNA samples in providing genetic material for CGP. The results of the CGP will be reviewed at a molecular tumor board. The proportion of participants that have a change in their treatment recommendations based on their individual genomic profiling will be assessed. Correlations between CGP and stage, prognosis, response to treatment and overall survival will also be investigated. This study will open to recruitment in 2020, with a target accrual of 150 A-PDAC patients within 36 months, with a 2-year follow-up. It is expected that the majority of participants will be those who have already consented for their tissue to be biobanked in the Victorian Pancreatic Cancer Biobank at the time of diagnostic EUS-FNA. Patients without archival or biobanked material that is suitable for CGP may be offered a EUS-FNA procedure for the purposes of obtaining fresh frozen material.

Discussion:

This trial is expected to provide crucial data regarding the feasibility of routine CGP of A-PDAC using EUS-FNA material. It will also provide important information about the impact of this methodology on patients’ survival.

Ambient temperature and genome-wide DNA methylation: A twin and family study in Australia

Little is known about the association between ambient temperature and DNA methylation, which is a potential biological process through which ambient temperature affects health. This study aimed to evaluate the association between ambient temperature and DNA methylation across human genome. We included 479 Australian women, including 132 twin pairs and 215 sisters of these twins. Blood-derived DNA methylation was measured using the HumanMethylation450 BeadChip array. Data on average ambient temperature during eight different exposure windows [lag0d (the blood draw day), lag0-7d (the current day and previous seven days prior to blood draw), lag0-14d, lag0-21d, lag0-28d, lag0-90d, lag0-180d, and lag0-365d)] was linked to each participant's home address. For each cytosine-guanine dinucleotide (CpG), we evaluated the association between its methylation level and temperature using generalized estimating equations (GEE), adjusting for important covariates. We used comb-p and DMRcate to identify differentially methylated regions (DMRs). We identified 31 CpGs at which blood DNA methylation were significantly associated with ambient temperature with false discovery rate [FDR] < 0.05. There were 82 significant DMRs identified by both comb-p (Sidak p-value < 0.01) and DMRcate (FDR < 0.01). Most of these CpGs and DMRs only showed association with temperature during one specific exposure window. These CpGs and DMRs were mapped to 85 genes. These related genes have been related to many human chronic diseases or phenotypes (e.g., diabetes, arthritis, breast cancer, depression, asthma, body height) in previous studies. The signals of short-term windows (lag0d and lag0-21d) showed enrichment in biological processes related to cell adhesion. In conclusion, short-, medium-, and long-term exposures to ambient temperature were all associated with blood DNA methylation, but the target genomic loci varied by exposure window. These differential methylation signals may serve as potential biomarkers to understand the health impacts of temperature.

Smoking Methylation Marks for Prediction of Urothelial Cancer Risk

BACKGROUND

Self-reported information may not accurately capture smoking exposure. We aimed to evaluate whether smoking-associated DNA methylation markers improve urothelial cell carcinoma (UCC) risk prediction.

METHODS

Conditional logistic regression was used to assess associations between blood-based methylation and UCC risk using two matched case-control samples: 404 pairs from the Melbourne Collaborative Cohort Study (MCCS) and 440 pairs from the Women's Health Initiative (WHI) cohort. Results were pooled using fixed-effects meta-analysis. We developed methylation-based predictors of UCC and evaluated their prediction accuracy on two replication data sets using the area under the curve (AUC).

RESULTS

The meta-analysis identified associations (P < 4.7 × 10^-5) for 29 of 1,061 smoking-associated methylation sites, but these were substantially attenuated after adjustment for self-reported smoking. Nominally significant associations (P < 0.05) were found for 387 (36%) and 86 (8%) of smoking-associated markers without/with adjustment for self-reported smoking, respectively, with same direction of association as with smoking for 387 (100%) and 79 (92%) markers. A Lasso-based predictor was associated with UCC risk in one replication data set in MCCS [N = 134; odds ratio per SD (OR) = 1.37; 95% CI, 1.00-1.90] after confounder adjustment; AUC = 0.66, compared with AUC = 0.64 without methylation information. Limited evidence of replication was found in the second testing data set in WHI (N = 440; OR = 1.09; 95% CI, 0.91-1.30).

CONCLUSIONS

Combination of smoking-associated methylation marks may provide some improvement to UCC risk prediction. Our findings need further evaluation using larger data sets.

IMPACT

DNA methylation may be associated with UCC risk beyond traditional smoking assessment and could contribute to some improvements in stratification of UCC risk in the general population.

Residential surrounding greenness and DNA methylation: An epigenome-wide association study

BACKGROUND

DNA methylation is a potential biological mechanism through which residential greenness affects health, but little is known about its association with greenness and whether the association could be modified by genetic background. We aimed to evaluate the association between surrounding greenness and genome-wide DNA methylation and potential gene-greenness interaction effects on DNA methylation.

METHODS

We measured blood-derived DNA methylation using the HumanMethylation450 BeadChip array (Illumina) for 479 Australian women, including 66 monozygotic, 66 dizygotic twin pairs, and 215 sisters of these twins. Surrounding greenness was represented by Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) within 300, 500, 1000 or 2000 m surrounding participants' home addresses. For each cytosine-guanine dinucleotide (CpG), the associations between its methylation level and NDVI or EVI were evaluated by generalized estimating equations, after adjusting for age, education, marital status, area-level socioeconomic status, smoking behavior, cell-type proportions, and familial clustering. We used comb-p and DMRcate to identify significant differentially methylated regions (DMRs). For each significant CpG, we evaluated the interaction effects of greenness and single-nucleotide polymorphisms (SNPs) within ±1 Mb window on its methylation level.

RESULTS

We found associations between surrounding greenness and blood DNA methylation for one CpG (cg04720477, mapped to the promoter region of CNP gene) with false discovery rate [FDR] < 0.05, and for another 9 CpGs with 0.05 ≤ FDR < 0.10. For two of these CpGs, we found 33 SNPs significantly (FDR < 0.05) modified the greenness-methylation association. There were 35 significant DMRs related to surrounding greenness that were identified by both comb-p (Sidak p-value < 0.01) and DMRcate (FDR < 0.01). Those CpGs and DMRs were mapped to genes related to many human diseases, such as mental health disorders and neoplasms as well as nutritional and metabolic diseases.

CONCLUSIONS

Surrounding greenness was associated with blood DNA methylation of many loci across human genome, and this association could be modified by genetic variations.

Independent prognostic impact of plasma NCOA2 alterations in metastatic castration-resistant prostate cancer

BACKGROUND

The androgen receptor (AR) pathway-associated gene nuclear receptor coactivator 2 (NCOA2) has an established oncogenic role in early prostate cancer and likewise is a driver of metastatic disease and castration-resistant prostate cancer. However, its significance as a biomarker in metastatic castration-resistant prostate cancer (mCRPC), both alone and in conjunction with co-occurring AR alterations using a liquid biopsy approach has not been investigated.

METHODS

Ninety-one patients were included in this study, (n = 68 receiving an androgen receptor pathway inhibitor and n = 23 receiving taxane chemotherapy). Up to 30 ml of peripheral blood was collected before commencing treatment from each patient. Plasma cell-free DNA, along with a matched germline sample, underwent targeted next-generation sequencing using a validated, highly sensitive in-house prostate cancer panel. Variants in AR and NCOA2 were identified and correlated with clinical outcomes.

RESULTS

Plasma AR and NCOA2 aberrations were identified in 35% and 13% of the cohort, respectively, whilst 8% had concurrent AR and NCOA2 alterations. NCOA2 copy number gain and any NCOA2 aberration predicted for lower prostate-specific antigen (PSA) response rates. Likewise, median overall survival was shorter for NCOA2 gain (10.1 vs. 18.3 months; p = .004), remaining significant after adjusting for covariates including circulating tumor DNA fraction and tumor suppressor gene alterations. Importantly, dual AR and NCOA2 aberrations were also associated with inferior outcomes, including no PSA responses in patients treated with AR pathway inhibitors (0% vs. 64%; p = .02).

CONCLUSIONS

These data highlight the importance of identifying multiple markers of AR pathway modulation in mCRPC and represent the first instance of the assessment of plasma NCOA2 status as a prognostic biomarker for standard-of-care therapies. Further assessment is warranted to determine if NCOA2 aberrations are a marker of primary resistance to AR pathway inhibitors.

Association of germline genetic variants with breast cancer-specific survival in patient subgroups defined by clinic-pathological variables related to tumor biology and type of systemic treatment

BACKGROUND

Given the high heterogeneity among breast tumors, associations between common germline genetic variants and survival that may exist within specific subgroups could go undetected in an unstratified set of breast cancer patients.

METHODS

We performed genome-wide association analyses within 15 subgroups of breast cancer patients based on prognostic factors, including hormone receptors, tumor grade, age, and type of systemic treatment. Analyses were based on 91,686 female patients of European ancestry from the Breast Cancer Association Consortium, including 7531 breast cancer-specific deaths over a median follow-up of 8.1 years. Cox regression was used to assess associations of common germline variants with 15-year and 5-year breast cancer-specific survival. We assessed the probability of these associations being true positives via the Bayesian false discovery probability (BFDP < 0.15).

RESULTS

Evidence of associations with breast cancer-specific survival was observed in three patient subgroups, with variant rs5934618 in patients with grade 3 tumors (15-year-hazard ratio (HR) [95% confidence interval (CI)] 1.32 [1.20, 1.45], P = 1.4E-08, BFDP = 0.01, per G allele); variant rs4679741 in patients with ER-positive tumors treated with endocrine therapy (15-year-HR [95% CI] 1.18 [1.11, 1.26], P = 1.6E-07, BFDP = 0.09, per G allele); variants rs1106333 (15-year-HR [95% CI] 1.68 [1.39,2.03], P = 5.6E-08, BFDP = 0.12, per A allele) and rs78754389 (5-year-HR [95% CI] 1.79 [1.46,2.20], P = 1.7E-08, BFDP = 0.07, per A allele), in patients with ER-negative tumors treated with chemotherapy.

CONCLUSIONS

We found evidence of four loci associated with breast cancer-specific survival within three patient subgroups. There was limited evidence for the existence of associations in other patient subgroups. However, the power for many subgroups is limited due to the low number of events. Even so, our results suggest that the impact of common germline genetic variants on breast cancer-specific survival might be limited.

Androgens alter the heterogeneity of small extracellular vesicles and the small RNA cargo in prostate cancer

Proliferation and survival of prostate cancer cells are driven by the androgen receptor (AR) upon binding to androgen steroid hormones. Manipulating the AR signalling axis is the focus for prostate cancer therapy; thus, it is crucial to understand the role of androgens and AR on extracellular vesicle (EV) secretion and cargo. In this study, we report that plasma-derived circulating vesicles consisting of CD9 and double-positive for CD9 and Prostate Specific Membrane Antigen (PSMA) are increased in patients with advanced metastatic prostate cancer, whereas double positives for CD9 and CD63 small extracellular vesicles (S-EVs) are significantly higher in patients with localised prostate cancer. Androgen manipulation by dihydrotestosterone (DHT) and the clinical antagonist enzalutamide (ENZ) altered the heterogeneity and size of CD9 positive S-EVs in AR expressing prostate cancer cells, while assessment of the total number and protein cargo of total S-EVs was unaltered across different treatment groups. Furthermore, hormone stimulation caused strong and specific effects on the small RNA cargo of S-EVs. A total of 543 small RNAs were found to be regulated by androgens including miR-19-3p and miR-361-5p. Analysis of S-EVs heterogeneity and small RNA cargo may provide clinical utility for prostate cancer and be informative to understand further the mechanism of resistance to androgen targeted therapy in castration-resistant prostate cancer.

Surrounding Greenness and Biological Aging Based on DNA Methylation: A Twin and Family Study in Australia

BACKGROUND

High surrounding greenness has many health benefits and might contribute to slower biological aging. However, very few studies have evaluated this from the perspective of epigenetics.

OBJECTIVES

We aimed to evaluate the association between surrounding greenness and biological aging based on DNA methylation.

METHODS

We derived Horvath's DNA methylation age (DNAmAge), Hannum's DNAmAge, PhenoAge, and GrimAge based on DNA methylation measured in peripheral blood samples from 479 Australian women in 130 families. Measures of DNAmAge acceleration (DNAmAgeAC) were derived from the residuals after regressing each DNAmAge metric on chronological age. Greenness was represented by satellite-derived Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) metrics within 300-, 500-, 1,000-, and 2,000-m buffers surrounding participant addresses. Greenness-DNAmAgeAC associations were estimated using a within-sibship design fitted by linear mixed effect models, adjusting for familial clustering and important covariates.

RESULTS

Greenness metrics were associated with significantly lower DNAmAgeAC based on GrimAge acceleration, suggesting slower biological aging with higher greenness based on both NDVI and EVI in 300-2,000m buffer areas. For example, each interquartile range increase in NDVI within 1,000m was associated with a 0.59 (95% CI: 0.18, 1.01)-year decrease in GrimAge acceleration. Greenness was also inversely associated with three of the eight components of GrimAge, specifically, DNA methylation-based surrogates of serum cystatin-C, serum growth differentiation factor 15, and smoking pack years. Associations between greenness and biological aging measured by Horvath's and Hannum's DNAmAgeAC were less consistent, and depended on neighborhood socioeconomic status. No significant associations were estimated for PhenoAge acceleration.

DISCUSSION

Higher surrounding greenness was associated with slower biological aging, as indicated by GrimAge age acceleration, in Australian women. Associations were also evident for three individual components of GrimAge, but were inconsistent for other measures of biological aging. Additional studies are needed to confirm our results. https://doi.org/10.1289/EHP8793.

Genetic insights into biological mechanisms governing human ovarian ageing

Reproductive longevity is essential for fertility and influences healthy ageing in women1,2, but insights into its underlying biological mechanisms and treatments to preserve it are limited. Here we identify 290 genetic determinants of ovarian ageing, assessed using normal variation in age at natural menopause (ANM) in about 200,000 women of European ancestry. These common alleles were associated with clinical extremes of ANM; women in the top 1% of genetic susceptibility have an equivalent risk of premature ovarian insufficiency to those carrying monogenic FMR1 premutations3. The identified loci implicate a broad range of DNA damage response (DDR) processes and include loss-of-function variants in key DDR-associated genes. Integration with experimental models demonstrates that these DDR processes act across the life-course to shape the ovarian reserve and its rate of depletion. Furthermore, we demonstrate that experimental manipulation of DDR pathways highlighted by human genetics increases fertility and extends reproductive life in mice. Causal inference analyses using the identified genetic variants indicate that extending reproductive life in women improves bone health and reduces risk of type 2 diabetes, but increases the risk of hormone-sensitive cancers. These findings provide insight into the mechanisms that govern ovarian ageing, when they act, and how they might be targeted by therapeutic approaches to extend fertility and prevent disease.

Transcriptomic changes in peripheral blood mononuclear cells with weight loss: systematic literature review and primary data synthesis

Background and objectives

Peripheral blood mononuclear cells (PBMCs) have shown promise as a tissue sensitive to subtle and possibly systemic transcriptomic changes, and as such may be useful in identifying responses to weight loss interventions. The primary aim was to comprehensively evaluate the transcriptomic changes that may occur during weight loss and to determine if there is a consistent response across intervention types in human populations of all ages.

Methods

Included studies were randomised control trials or cohort studies that administered an intervention primarily designed to decrease weight in any overweight or obese human population. A systematic search of the literature was conducted to obtain studies and gene expression databases were interrogated to locate corresponding transcriptomic datasets. Datasets were normalised using the ArrayAnalysis online tool and differential gene expression was determined using the limma package in R. Over-represented pathways were explored using the PathVisio software. Heatmaps and hierarchical clustering were utilised to visualise gene expression.

Results

Seven papers met the inclusion criteria, five of which had raw gene expression data available. Of these, three could be grouped into high responders (HR, ≥ 5% body weight loss) and low responders (LR). No genes were consistently differentially expressed between high and low responders across studies. Adolescents had the largest transcriptomic response to weight loss followed by adults who underwent bariatric surgery. Seven pathways were altered in two out of four studies following the intervention and the pathway ‘cytoplasmic ribosomal proteins’ (WikiPathways: WP477) was altered between HR and LR at baseline in the two datasets with both groups. Pathways related to ‘toll-like receptor signalling’ were altered in HR response to the weight loss intervention in two out of three datasets.

Conclusions

Transcriptomic changes in PBMCs do occur in response to weight change. Transparent and standardised data reporting is needed to realise the potential of transcriptomics for investigating phenotypic features.

Registration number

PROSPERO: CRD42019106582

Supplementary Information

The online version contains supplementary material available at 10.1186/s12263-021-00692-6.

Functional annotation of the 2q35 breast cancer risk locus implicates a structural variant in influencing activity of a long-range enhancer element

A combination of genetic and functional approaches has identified three independent breast cancer risk loci at 2q35. A recent fine-scale mapping analysis to refine these associations resulted in 1 (signal 1), 5 (signal 2), and 42 (signal 3) credible causal variants at these loci. We used publicly available in silico DNase I and ChIP-seq data with in vitro reporter gene and CRISPR assays to annotate signals 2 and 3. We identified putative regulatory elements that enhanced cell-type-specific transcription from the IGFBP5 promoter at both signals (30- to 40-fold increased expression by the putative regulatory element at signal 2, 2- to 3-fold by the putative regulatory element at signal 3). We further identified one of the five credible causal variants at signal 2, a 1.4 kb deletion (esv3594306), as the likely causal variant; the deletion allele of this variant was associated with an average additional increase in IGFBP5 expression of 1.3-fold (MCF-7) and 2.2-fold (T-47D). We propose a model in which the deletion allele of esv3594306 juxtaposes two transcription factor binding regions (annotated by estrogen receptor alpha ChIP-seq peaks) to generate a single extended regulatory element. This regulatory element increases cell-type-specific expression of the tumor suppressor gene IGFBP5 and, thereby, reduces risk of estrogen receptor-positive breast cancer (odds ratio = 0.77, 95% CI 0.74-0.81, p = 3.1 × 10-31).

Genome-wide association studies identify 137 genetic loci for DNA methylation biomarkers of aging

BACKGROUND

Biological aging estimators derived from DNA methylation data are heritable and correlate with morbidity and mortality. Consequently, identification of genetic and environmental contributors to the variation in these measures in populations has become a major goal in the field.

RESULTS

Leveraging DNA methylation and SNP data from more than 40,000 individuals, we identify 137 genome-wide significant loci, of which 113 are novel, from genome-wide association study (GWAS) meta-analyses of four epigenetic clocks and epigenetic surrogate markers for granulocyte proportions and plasminogen activator inhibitor 1 levels, respectively. We find evidence for shared genetic loci associated with the Horvath clock and expression of transcripts encoding genes linked to lipid metabolism and immune function. Notably, these loci are independent of those reported to regulate DNA methylation levels at constituent clock CpGs. A polygenic score for GrimAge acceleration showed strong associations with adiposity-related traits, educational attainment, parental longevity, and C-reactive protein levels.

CONCLUSION

This study illuminates the genetic architecture underlying epigenetic aging and its shared genetic contributions with lifestyle factors and longevity.

Genome-wide homozygosity and risk of four non-Hodgkin lymphoma subtypes

AIM

Recessive genetic variation is thought to play a role in non-Hodgkin lymphoma (NHL) etiology. Runs of homozygosity (ROH), defined based on long, continuous segments of homozygous SNPs, can be used to estimate both measured and unmeasured recessive genetic variation. We sought to examine genome-wide homozygosity and NHL risk.

METHODS

We used data from eight genome-wide association studies of four common NHL subtypes: 3061 chronic lymphocytic leukemia (CLL), 3814 diffuse large B-cell lymphoma (DLBCL), 2784 follicular lymphoma (FL), and 808 marginal zone lymphoma (MZL) cases, as well as 9374 controls. We examined the effect of homozygous variation on risk by: (1) estimating the fraction of the autosome containing runs of homozygosity (FROH); (2) calculating an inbreeding coefficient derived from the correlation among uniting gametes (F3); and (3) examining specific autosomal regions containing ROH. For each, we calculated beta coefficients and standard errors using logistic regression and combined estimates across studies using random-effects meta-analysis.

RESULTS

We discovered positive associations between FROH and CLL (β = 21.1, SE = 4.41, P = 1.6 × 10-6) and FL (β = 11.4, SE = 5.82, P = 0.02) but not DLBCL (P = 1.0) or MZL (P = 0.91). For F3, we observed an association with CLL (β = 27.5, SE = 6.51, P = 2.4 × 10-5). We did not find evidence of associations with specific ROH, suggesting that the associations observed with FROH and F3 for CLL and FL risk were not driven by a single region of homozygosity.

CONCLUSION

Our findings support the role of recessive genetic variation in the etiology of CLL and FL; additional research is needed to identify the specific loci associated with NHL risk.

Methylation marks of prenatal exposure to maternal smoking and risk of cancer in adulthood

BACKGROUND

Prenatal exposure to maternal smoking is detrimental to child health but its association with risk of cancer has seldom been investigated. Maternal smoking induces widespread and long-lasting DNA methylation changes, which we study here for association with risk of cancer in adulthood.

METHODS

Eight prospective case-control studies nested within the Melbourne Collaborative Cohort Study were used to assess associations between maternal-smoking-associated methylation marks in blood and risk of several cancers: breast (n = 406 cases), colorectal (n = 814), gastric (n = 166), kidney (n = 139), lung (n = 327), prostate (n = 847) and urothelial (n = 404) cancer and B-cell lymphoma (n = 426). We used conditional logistic regression models to estimate odds ratios (ORs) and 95% confidence intervals (CIs) for associations between cancer and five methylation scores calculated as weighted averages for 568, 19, 15, 28 and 17 CpG sites. Models were adjusted for confounders, including personal smoking history (smoking status, pack-years, age at starting and quitting) and methylation scores for personal smoking.

RESULTS

All methylation scores for maternal smoking were strongly positively associated with risk of urothelial cancer. Risk estimates were only slightly attenuated after adjustment for smoking history, other potential confounders and methylation scores for personal smoking. Potential negative associations were observed with risk of lung cancer and B-cell lymphoma. No associations were observed for other cancers.

CONCLUSIONS

We found that methylation marks of prenatal exposure to maternal smoking are associated with increased risk of urothelial cancer. Our study demonstrates the potential for using DNA methylation to investigate the impact of early-life, unmeasured exposures on later-life cancer risk.

Prospective Evaluation of the Addition of Polygenic Risk Scores to Breast Cancer Risk Models

Background

The Breast and Ovarian Analysis of Disease Incidence and Carrier Estimation Algorithm and the International Breast Cancer Intervention Study breast cancer risk models are used to provide advice on screening intervals and chemoprevention. We evaluated the performance of these models, which now incorporate polygenic risk scores (PRSs), using a prospective cohort study.

Methods

We used a case-cohort design, involving women in the Melbourne Collaborative Cohort Study aged 50-75 years when surveyed in 2003-2007, of whom 408 had a first primary breast cancer diagnosed within 10 years (cases), and 2783 were from the subcohort. Ten-year risks were calculated based on lifestyle factors, family history data, and a 313-variant PRS. Discrimination was assessed using a C-statistic compared with 0.50 and calibration using the ratio of expected to observed number of cases (E/O).

Results

When the PRS was added to models with lifestyle factors and family history, the C-statistic (95% confidence interval [CI]) increased from 0.57 (0.54 to 0.60) to 0.62 (0.60 to 0.65) using IBIS and from 0.56 (0.53 to 0.59) to 0.62 (0.59 to 0.64) using BOADICEA. IBIS underpredicted risk (E/O = 0.62, 95% CI = 0.48 to 0.80) for women in the lowest risk category (<1.7%) and overpredicted risk (E/O = 1.40, 95% CI = 1.18 to 1.67) in the highest risk category (≥5%), using the Hosmer-Lemeshow test for calibration in quantiles of risk and a 2-sided P value less than  .001. BOADICEA underpredicted risk (E/O = 0.82, 95% CI = 0.67 to 0.99) in the second highest risk category (3.4%-5%); the Hosmer-Lemeshow test and a 2-sided P value was equal to .02.

Conclusions

Although the inclusion of a 313 genetic variant PRS doubles discriminatory accuracy (relative to reference 0.50), models with and without this PRS have relatively modest discrimination and might require recalibration before their clinical and wider use are promoted.

DNA Methylation Signatures and the Contribution of Age-Associated Methylomic Drift to Carcinogenesis in Early-Onset Colorectal Cancer

Simple Summary

The role of DNA methylation in the carcinogenesis of colorectal cancer (CRC) diagnosed <50 years of age (early-onset CRC or EOCRC) is currently unknown. In this study, we investigated the genome-wide DNA methylation of 97 tumour and 54 normal colonic mucosa samples from people with EOCRC with the aim of identifying unique DNA methylation signatures and determining the role of ageing-related DNA methylation drift and age-acceleration in EOCRC aetiology. We found extensive DNA methylation alterations associated with EOCRC carcinogenesis, including a unique signature comprising 234 loci compared with CRCs from people >50 years of age. CpGs that undergo ageing-related methylation drift were significantly altered in EOCRC, and accelerated ageing was also evident in normal mucosa from people with EOCRC. Our study is the first study to identify unique DNA methylation changes in EOCRC, contributing novel information that may aid future efforts towards EOCRC prevention.

Abstract

We investigated aberrant DNA methylation (DNAm) changes and the contribution of ageing-associated methylomic drift and age acceleration to early-onset colorectal cancer (EOCRC) carcinogenesis. Genome-wide DNAm profiling using the Infinium HM450K on 97 EOCRC tumour and 54 normal colonic mucosa samples was compared with: (1) intermediate-onset CRC (IOCRC; diagnosed between 50–70 years; 343 tumour and 35 normal); and (2) late-onset CRC (LOCRC; >70 years; 318 tumour and 40 normal). CpGs associated with age-related methylation drift were identified using a public dataset of 231 normal mucosa samples from people without CRC. DNAm-age was estimated using epiTOC2. Common to all three age-of-onset groups, 88,385 (20% of all CpGs) CpGs were differentially methylated between tumour and normal mucosa. We identified 234 differentially methylated genes that were unique to the EOCRC group; 13 of these DMRs/genes were replicated in EOCRC compared with LOCRCs from TCGA. In normal mucosa from people without CRC, we identified 28,154 CpGs that undergo ageing-related DNAm drift, and of those, 65% were aberrantly methylated in EOCRC tumours. Based on the mitotic-based DNAm clock epiTOC2, we identified age acceleration in normal mucosa of people with EOCRC compared with normal mucosa from the IOCRC, LOCRC groups (p = 3.7 × 10⁻¹⁶) and young people without CRC (p = 5.8 × 10⁻⁶). EOCRC acquires unique DNAm alterations at 234 loci. CpGs associated with ageing-associated drift were widely affected in EOCRC without needing the decades-long accrual of DNAm drift as commonly seen in intermediate- and late-onset CRCs. Accelerated ageing in normal mucosa from people with EOCRC potentially underlies the earlier age of diagnosis in CRC carcinogenesis.

Mismatch repair gene pathogenic germline variants in a population-based cohort of breast cancer

The advent of gene panel testing is challenging the previous practice of using clinically defined cancer family syndromes to inform single-gene genetic screening. Individual and family cancer histories that would have previously indicated testing of a single gene or a small number of related genes are now, increasingly, leading to screening across gene panels that contain larger numbers of genes. We have applied a gene panel test that included four DNA mismatch repair (MMR) genes (MLH1, MSH2, MSH6 and PMS2) to an Australian population-based case-control-family study of breast cancer. Altogether, eight pathogenic variants in MMR genes were identified: six in 1421 case-families (0.4%, 4 MSH6 and 2 PMS2) and two in 833 control-families (0.2%, one each of MLH1 and MSH2). This testing highlights the current and future challenges for clinical genetics in the context of anticipated gene panel-based population-based screening that includes the MMR genes. This testing is likely to provide additional opportunities for cancer prevention via cascade testing for Lynch syndrome and precision medicine for breast cancer treatment.

Population-based estimates of the age-specific cumulative risk of breast cancer for pathogenic variants in CHEK2: Findings from the Australian Breast Cancer Family Registry

551

Background: Case-control studies of breast cancer have consistently shown that pathogenic variants in CHEK2 are associated with about a 3-fold increased risk of breast cancer. Information about the recurrent protein truncating variant CHEK2c.1100delC dominates this estimate. There have been no formal estimates of age-specific cumulative risk of breast cancer for all CHEK2 pathogenic (including likely pathogenic) variants combined. Methods: We conducted a genetic screen of CHEK2 in an Australian population-based case-control-family study of breast cancer. This study is focused on disease at an early age and participants were unselected for family history. The age-specific cumulative risk (penetrance) of breast cancer was estimated using segregation analysis. Results: The estimated hazard ratio for carriers of pathogenic CHEK2 variants (combined) was 4.9 (95% CI 2.5-9.5; p < 0.0001) relative to non-carriers. The HR for carriers of the CHEK2 c.1100delC variant was estimated to be 3.5 (95% CI 1.02-11.6) and the HR for carriers of all other CHEK2 variants combined was estimated to be 5.7 (95% CI 2.5-12.9). The age-specific cumulative risk of breast cancer was estimated to be 18% (95% CI 11-30%) and 33% (95% CI 21-48%) to age 60 and 80 years, respectively. Conclusions: These findings provide important information for the clinical management of breast cancer risk for women carrying pathogenic variants in CHEK2.

Author Correction: A case-only study to identify genetic modifiers of breast cancer risk for BRCA1/BRCA2 mutation carriers

A Correction to this paper has been published: https://doi.org/10.1038/s41467-021-23162-4

Gene-Environment Interactions Relevant to Estrogen and Risk of Breast Cancer: Can Gene-Environment Interactions Be Detected Only among Candidate SNPs from Genome-Wide Association Studies?

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.

Evaluation of the association of heterozygous germline variants in NTHL1 with breast cancer predisposition: an international multi-center study of 47,180 subjects

Bi-allelic loss-of-function (LoF) variants in the base excision repair (BER) gene NTHL1 cause a high-risk hereditary multi-tumor syndrome that includes breast cancer, but the contribution of heterozygous variants to hereditary breast cancer is unknown. An analysis of 4985 women with breast cancer, enriched for familial features, and 4786 cancer-free women revealed significant enrichment for NTHL1 LoF variants. Immunohistochemistry confirmed reduced NTHL1 expression in tumors from heterozygous carriers but the NTHL1 bi-allelic loss characteristic mutational signature (SBS 30) was not present. The analysis was extended to 27,421 breast cancer cases and 19,759 controls from 10 international studies revealing 138 cases and 93 controls with a heterozygous LoF variant (OR 1.06, 95% CI: 0.82-1.39) and 316 cases and 179 controls with a missense variant (OR 1.31, 95% CI: 1.09-1.57). Missense variants selected for deleterious features by a number of in silico bioinformatic prediction tools or located within the endonuclease III functional domain showed a stronger association with breast cancer. Somatic sequencing of breast cancers from carriers indicated that the risk associated with NTHL1 appears to operate through haploinsufficiency, consistent with other described low-penetrance breast cancer genes. Data from this very large international multicenter study suggests that heterozygous pathogenic germline coding variants in NTHL1 may be associated with low- to moderate- increased risk of breast cancer.

Rare genetic variants: making the connection with breast cancer susceptibility

The practice of clinical genetics in the context of breast cancer predisposition has reached another critical point in its evolution. For the past two decades, genetic testing offered to women attending clinics has been limited to BRCA1 and BRCA2 unless other syndromic indicators have been evident (e.g. PTEN and TP53 for Cowden and Li-Fraumeni syndrome, respectively). Women (and their families) who are concerned about their personal and/or family history of breast and ovarian cancer have enthusiastically engaged with clinical genetics services, anticipating a genetic cause for their cancer predisposition will be identified and to receive clinical guidance for their risk management and treatment options. Genetic testing laboratories have demonstrated similar enthusiasm for transitioning from single gene to gene panel testing that now provide opportunities for the large number of women found not to carry mutations in BRCA1 and BRCA2, enabling them to undergo additional genetic testing. However, these panel tests have limited clinical utility until more is understood about the cancer risks (if any) associated with the genetic variation observed in the genes included on these panels. New data is urgently needed to improve the interpretation of the genetic variation data that is already reported from these panels and to inform the selection of genes included in gene panel tests in the future. To address this issue, large internationally coordinated research studies are required to provide the evidence-base from which clinical genetics for breast cancer susceptibility can be practiced in the era of gene panel testing and oncogenetic practice.

Two significant steps associated with this process include i) validating the genes on these panels (and those likely to be added in the future) as bona fide1 breast cancer predisposition genes and ii) interpreting the variation, on a variant-by-variant basis in terms of their likely “pathogenicity”—a process commonly referred to as “variant classification” that will enable this new genetic information to be used at an individual level in clinical genetics services. Neither of these fundamental steps have been achieved for the majority of genes included on the panels.

We are thus at a critical point for translational research in breast cancer clinical genetics—how can rare genetic variants be interpreted such that they can be used in clinical genetics services and oncogenetic practice to identify and to inform the management of families that carry these variants?

Epigenetic Drift Association with Cancer Risk and Survival, and Modification by Sex

Simple Summary

Ageing is the strongest cancer risk factor, and men and women exhibit different risk profiles in terms of incidence and survival. DNA methylation is known to strongly vary by age and sex. Epigenetic drift refers to age-related DNA methylation changes and the tendency for increasing discordance between epigenomes over time, but it remains unknown to what extent the epigenetic drift contributes to cancer risk and survival. The aims of this study were to identify age-associated, sex-associated and sexually dimorphic age-associated (age-by-sex-associated) DNA methylation markers and investigate whether age- and age-by-sex-associated markers are associated with cancer risk and survival. Our study, which used a total of 2754 matched case–control pairs with DNA methylation in pre-diagnostic blood, is the first large study to examine the association between sex-specific epigenetic drift and cancer development and progression. The results may be useful for cancer early diagnosis and prediction of prognosis.

Abstract

To investigate age- and sex-specific DNA methylation alterations related to cancer risk and survival, we used matched case–control studies of colorectal (n = 835), gastric (n = 170), kidney (n = 143), lung (n = 332), prostate (n = 869) and urothelial (n = 428) cancers, and mature B-cell lymphoma (n = 438). Linear mixed-effects models were conducted to identify age-, sex- and age-by-sex-associated methylation markers using a discovery (controls)-replication (cases) strategy. Replication was further examined using summary statistics from Generation Scotland (GS). Associations between replicated markers and risk of and survival from cancer were assessed using conditional logistic regression and Cox models (hazard ratios (HR)), respectively. We found 32,659, 23,141 and 48 CpGs with replicated associations for age, sex and age-by-sex, respectively. The replication rates for these CpGs using GS summary data were 94%, 86% and 91%, respectively. Significant associations for cancer risk and survival were identified at some individual age-related CpGs. Opposite to previous findings using epigenetic clocks, there was a strong negative trend in the association between epigenetic drift and risk of colorectal cancer. Methylation at two CpGs overlapping TMEM49 and ARX genes was associated with survival of overall (HR = 0.91, p = 7.7 × 10⁻⁴) and colorectal (HR = 1.52, p = 1.8 × 10⁻⁴) cancer, respectively, with significant age-by-sex interaction. Our results may provide markers for cancer early detection and prognosis prediction.