Accuracy of the No-Biopsy Approach for the Diagnosis of Celiac Disease in Adults: A Systematic Review and Meta-Analysis

BACKGROUND & AIMS

Current international guidelines recommend duodenal biopsies to confirm the diagnosis of celiac disease in adult patients. However, growing evidence suggests that immunoglobulin A (IgA) anti-tissue transglutaminase (tTg) antibody levels ≥10 times the upper limit of normal (ULN) can accurately predict celiac disease, eliminating the need for biopsy. We performed a systematic review and meta-analysis to evaluate the accuracy of the no-biopsy approach to confirm the diagnosis of celiac disease in adults.

METHODS

We systematically searched MEDLINE, EMBASE, Cochrane Library, and Web of Science from January 1998 to October 2023 for studies reporting the sensitivity and specificity of IgA-tTG ≥10×ULN against duodenal biopsies (Marsh grade ≥2) in adults with suspected celiac disease. We used a bivariate random effects model to calculate the summary estimates of sensitivity, specificity, and positive and negative likelihood ratios. The positive and negative likelihood ratios were used to calculate the positive predictive value of the no-biopsy approach across different pretest probabilities of celiac disease. The methodological quality of the included studies was evaluated using the QUADAS-2 tool. This study was registered with PROSPERO, number CRD42023398812.

RESULTS

A total of 18 studies comprising 12,103 participants from 15 countries were included. The pooled prevalence of biopsy-proven celiac disease in the included studies was 62% (95% confidence interval [CI], 40%-83%). The proportion of patients with IgA-tTG ≥10×ULN was 32% (95% CI, 24%-40%). The summary sensitivity of IgA-tTG ≥10×ULN was 51% (95% CI, 42%-60%), and the summary specificity was 100% (95% CI, 98%-100%). The area under the summary receiver operating characteristic curve was 0.83 (95% CI, 0.77 - 0.89). The positive predictive value of the no-biopsy approach to identify patients with celiac disease was 65%, 88%, 95%, and 99% if celiac disease prevalence was 1%, 4%, 10%, and 40%, respectively. Between-study heterogeneity was moderate (I2 =30.3%), and additional sensitivity analyses did not significantly alter our findings. Only 1 study had a low risk of bias across all domains.

CONCLUSION

The results of this meta-analysis suggest that selected adult patients with IgA-tTG ≥10×ULN and a moderate to high pretest probability of celiac disease could be diagnosed without undergoing invasive endoscopy and duodenal biopsy.

New entity of adult ultra-short coeliac disease: the first international cohort and case-control study

BACKGROUND

Ultra-short coeliac disease (USCD) is defined as villous atrophy only present in the duodenal bulb (D1) with concurrent positive coeliac serology. We present the first, multicentre, international study of patients with USCD.

METHODS

Patients with USCD were identified from 10 tertiary hospitals (6 from Europe, 2 from Asia, 1 from North America and 1 from Australasia) and compared with age-matched and sex-matched patients with conventional coeliac disease.

FINDINGS

Patients with USCD (n=137, median age 27 years, IQR 21-43 years; 73% female) were younger than those with conventional coeliac disease (27 vs 38 years, respectively, p<0.001). Immunoglobulin A-tissue transglutaminase (IgA-tTG) titres at index gastroscopy were lower in patients with USCD versus conventional coeliac disease (1.8×upper limit of normal (ULN) (IQR 1.1-5.9) vs 12.6×ULN (IQR 3.3-18.3), p<0.001).Patients with USCD had the same number of symptoms overall (median 3 (IQR 2-4) vs 3 (IQR 1-4), p=0.875). Patients with USCD experienced less iron deficiency (41.8% vs 22.4%, p=0.006).Both USCD and conventional coeliac disease had the same intraepithelial lymphocytes immunophenotype staining pattern; positive for CD3 and CD8, but not CD4.At follow-up having commenced a gluten-free diet (GFD) (median of 1181 days IQR: 440-2160 days) both USCD and the age-matched and sex-matched controls experienced a similar reduction in IgA-tTG titres (0.5 ULN (IQR 0.2-1.4) vs 0.7 ULN (IQR 0.2-2.6), p=0.312). 95.7% of patients with USCD reported a clinical improvement in their symptoms.

INTERPRETATION

Patients with USCD are younger, have a similar symptomatic burden and benefit from a GFD. This study endorses the recommendation of D1 sampling as part of the endoscopic coeliac disease diagnostic workup.

Differences in polygenic score distributions in European ancestry populations: implications for breast cancer risk prediction

The 313-variant polygenic risk score (PRS313) provides a promising tool for breast cancer risk prediction. However, evaluation of the PRS313 across different European populations which could influence risk estimation has not been performed. Here, we explored the distribution of PRS313 across European populations using genotype data from 94,072 females without breast cancer, of European-ancestry from 21 countries participating in the Breast Cancer Association Consortium (BCAC) and 225,105 female participants from the UK Biobank. The mean PRS313 differed markedly across European countries, being highest in south-eastern Europe and lowest in north-western Europe. Using the overall European PRS313 distribution to categorise individuals leads to overestimation and underestimation of risk in some individuals from south-eastern and north-western countries, respectively. Adjustment for principal components explained most of the observed heterogeneity in mean PRS. Country-specific PRS distributions may be used to calibrate risk categories in individuals from different countries.

Evaluation of European-based polygenic risk score for breast cancer in Ashkenazi Jewish women in Israel

BACKGROUND

Polygenic risk score (PRS), calculated based on genome-wide association studies (GWASs), can improve breast cancer (BC) risk assessment. To date, most BC GWASs have been performed in individuals of European (EUR) ancestry, and the generalisation of EUR-based PRS to other populations is a major challenge. In this study, we examined the performance of EUR-based BC PRS models in Ashkenazi Jewish (AJ) women.

METHODS

We generated PRSs based on data on EUR women from the Breast Cancer Association Consortium (BCAC). We tested the performance of the PRSs in a cohort of 2161 AJ women from Israel (1437 cases and 724 controls) from BCAC (BCAC cohort from Israel (BCAC-IL)). In addition, we tested the performance of these EUR-based BC PRSs, as well as the established 313-SNP EUR BC PRS, in an independent cohort of 181 AJ women from Hadassah Medical Center (HMC) in Israel.

RESULTS

In the BCAC-IL cohort, the highest OR per 1 SD was 1.56 (±0.09). The OR for AJ women at the top 10% of the PRS distribution compared with the middle quintile was 2.10 (±0.24). In the HMC cohort, the OR per 1 SD of the EUR-based PRS that performed best in the BCAC-IL cohort was 1.58±0.27. The OR per 1 SD of the commonly used 313-SNP BC PRS was 1.64 (±0.28).

CONCLUSIONS

Extant EUR GWAS data can be used for generating PRSs that identify AJ women with markedly elevated risk of BC and therefore hold promise for improving BC risk assessment in AJ women.

Associations between AI-Assisted Tumor Amphiregulin and Epiregulin IHC and Outcomes from Anti-EGFR Therapy in the Routine Management of Metastatic Colorectal Cancer

PURPOSE

High tumor production of the EGFR ligands, amphiregulin (AREG) and epiregulin (EREG), predicted benefit from anti-EGFR therapy for metastatic colorectal cancer (mCRC) in a retrospective analysis of clinical trial data. Here, AREG/EREG IHC was analyzed in a cohort of patients who received anti-EGFR therapy as part of routine care, including key clinical contexts not investigated in the previous analysis.

EXPERIMENTAL DESIGN

Patients who received panitumumab or cetuximab ± chemotherapy for treatment of RAS wild-type mCRC at eight UK cancer centers were eligible. Archival formalin-fixed paraffin-embedded tumor tissue was analyzed for AREG and EREG IHC in six regional laboratories using previously developed artificial intelligence technologies. Primary endpoints were progression-free survival (PFS) and overall survival (OS).

RESULTS

A total of 494 of 541 patients (91.3%) had adequate tissue for analysis. A total of 45 were excluded after central extended RAS testing, leaving 449 patients in the primary analysis population. After adjustment for additional prognostic factors, high AREG/EREG expression (n = 360; 80.2%) was associated with significantly prolonged PFS [median: 8.5 vs. 4.4 months; HR, 0.73; 95% confidence interval (CI), 0.56-0.95; P = 0.02] and OS [median: 16.4 vs. 8.9 months; HR, 0.66 95% CI, 0.50-0.86; P = 0.002]. The significant OS benefit was maintained among patients with right primary tumor location (PTL), those receiving cetuximab or panitumumab, those with an oxaliplatin- or irinotecan-based chemotherapy backbone, and those with tumor tissue obtained by biopsy or surgical resection.

CONCLUSIONS

High tumor AREG/EREG expression was associated with superior survival outcomes from anti-EGFR therapy in mCRC, including in right PTL disease. AREG/EREG IHC assessment could aid therapeutic decisions in routine practice. See related commentary by Randon and Pietrantonio, p. 4021.

Polymorphisms in genes of melatonin biosynthesis and signaling support the light-at-night hypothesis for breast cancer

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.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Association of the CHEK2 c.1100delC variant, radiotherapy, and systemic treatment with contralateral breast cancer risk and breast cancer-specific survival

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.

Association of the CHEK2 c.1100delC variant, radiotherapy, and systemic treatment with contralateral breast cancer risk and breast cancer-specific survival

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

Physical activity, sedentary time and breast cancer risk: a Mendelian randomisation study

OBJECTIVES

Physical inactivity and sedentary behaviour are associated with higher breast cancer risk in observational studies, but ascribing causality is difficult. Mendelian randomisation (MR) assesses causality by simulating randomised trial groups using genotype. We assessed whether lifelong physical activity or sedentary time, assessed using genotype, may be causally associated with breast cancer risk overall, pre/post-menopause, and by case-groups defined by tumour characteristics.

METHODS

We performed two-sample inverse-variance-weighted MR using individual-level Breast Cancer Association Consortium case-control data from 130 957 European-ancestry women (69 838 invasive cases), and published UK Biobank data (n=91 105-377 234). Genetic instruments were single nucleotide polymorphisms (SNPs) associated in UK Biobank with wrist-worn accelerometer-measured overall physical activity (nsnps=5) or sedentary time (nsnps=6), or accelerometer-measured (nsnps=1) or self-reported (nsnps=5) vigorous physical activity.

RESULTS

Greater genetically-predicted overall activity was associated with lower breast cancer overall risk (OR=0.59; 95% confidence interval (CI) 0.42 to 0.83 per-standard deviation (SD;~8 milligravities acceleration)) and for most case-groups. Genetically-predicted vigorous activity was associated with lower risk of pre/perimenopausal breast cancer (OR=0.62; 95% CI 0.45 to 0.87,≥3 vs. 0 self-reported days/week), with consistent estimates for most case-groups. Greater genetically-predicted sedentary time was associated with higher hormone-receptor-negative tumour risk (OR=1.77; 95% CI 1.07 to 2.92 per-SD (~7% time spent sedentary)), with elevated estimates for most case-groups. Results were robust to sensitivity analyses examining pleiotropy (including weighted-median-MR, MR-Egger).

CONCLUSION

Our study provides strong evidence that greater overall physical activity, greater vigorous activity, and lower sedentary time are likely to reduce breast cancer risk. More widespread adoption of active lifestyles may reduce the burden from the most common cancer in women.

Incorporating progesterone receptor expression into the PREDICT breast prognostic model

BACKGROUND

Predict Breast (www.predict.nhs.uk) is an online prognostication and treatment benefit tool for early invasive breast cancer. The aim of this study was to incorporate the prognostic effect of progesterone receptor (PR) status into a new version of PREDICT and to compare its performance to the current version (2.2).

METHOD

The prognostic effect of PR status was based on the analysis of data from 45,088 European patients with breast cancer from 49 studies in the Breast Cancer Association Consortium. Cox proportional hazard models were used to estimate the hazard ratio for PR status. Data from a New Zealand study of 11,365 patients with early invasive breast cancer were used for external validation. Model calibration and discrimination were used to test the model performance.

RESULTS

Having a PR-positive tumour was associated with a 23% and 28% lower risk of dying from breast cancer for women with oestrogen receptor (ER)-negative and ER-positive breast cancer, respectively. The area under the ROC curve increased with the addition of PR status from 0.807 to 0.809 for patients with ER-negative tumours (p = 0.023) and from 0.898 to 0.902 for patients with ER-positive tumours (p = 2.3 × 10-6) in the New Zealand cohort. Model calibration was modest with 940 observed deaths compared to 1151 predicted.

CONCLUSION

The inclusion of the prognostic effect of PR status to PREDICT Breast has led to an improvement of model performance and more accurate absolute treatment benefit predictions for individual patients. Further studies should determine whether the baseline hazard function requires recalibration.

Non-Responsive and Refractory Coeliac Disease: Experience from the NHS England National Centre

We characterised the aetiology of non-responsive coeliac disease (NRCD) and provided contemporary mortality data in refractory coeliac disease (RCD) from our centre. We also measured urine gluten immunogenic peptides (GIPs) in patients with established RCD1 to evaluate gluten exposure in these individuals. Methods: This was a longitudinal cohort study conducted in Sheffield, UK. Between 1998 and 2019, we evaluated 285 adult (≥16 years) patients with NRCD or RCD. Patients with established RCD1 and persisting mucosal inflammation and/or ongoing symptoms provided three urine samples for GIP analysis. Results: The most common cause of NRCD across the cohort was gluten exposure (72/285; 25.3%). RCD accounted for 65/285 patients (22.8%), 54/65 patients (83.1%) had RCD1 and 11/65 patients (16.9%) had RCD2. The estimated 5-year survival was 90% for RCD1 and 58% for RCD2 (p = 0.016). A total of 36/54 (66.7%) patients with RCD1 underwent urinary GIP testing and 17/36 (47.2%) had at least one positive urinary GIP test. Conclusion: The contemporary mortality data in RCD2 remains poor; patients with suspected RCD2 should be referred to a recognised national centre for consideration of novel therapies. The high frequency of urinary GIP positivity suggests that gluten exposure may be common in RCD1; further studies with matched controls are warranted to assess this further.

Rare germline copy number variants (CNVs) and breast cancer risk

Germline copy number variants (CNVs) are pervasive in the human genome but potential disease associations with rare CNVs have not been comprehensively assessed in large datasets. We analysed rare CNVs in genes and non-coding regions for 86,788 breast cancer cases and 76,122 controls of European ancestry with genome-wide array data. Gene burden tests detected the strongest association for deletions in BRCA1 (P = 3.7E-18). Nine other genes were associated with a p-value < 0.01 including known susceptibility genes CHEK2 (P = 0.0008), ATM (P = 0.002) and BRCA2 (P = 0.008). Outside the known genes we detected associations with p-values < 0.001 for either overall or subtype-specific breast cancer at nine deletion regions and four duplication regions. Three of the deletion regions were in established common susceptibility loci. To the best of our knowledge, this is the first genome-wide analysis of rare CNVs in a large breast cancer case-control dataset. We detected associations with exonic deletions in established breast cancer susceptibility genes. We also detected suggestive associations with non-coding CNVs in known and novel loci with large effects sizes. Larger sample sizes will be required to reach robust levels of statistical significance.

Common variants in breast cancer risk loci predispose to distinct tumor subtypes

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.

Clinical classification and long-term outcomes of seronegative coeliac disease: a 20-year multicentre follow-up study

BACKGROUND

Seronegative coeliac disease is poorly defined.

AIMS

To study clinical phenotypes and long-term outcomes of seronegative coeliac disease in a multicentre cohort over 20 years.

METHODS

Seronegative coeliac disease was diagnosed in HLA-DQ2/DQ8-positive patients with villous atrophy (VA), negative IgA endomysial (EmA), tissue transglutaminase (tTG) and deamidated-gliadin antibodies (DGP), clinical and histological response to a gluten-free diet (GFD), and no alternative causes for VA. In patients with IgA deficiency, coeliac disease was diagnosed through VA, positive IgG EmA/tTG/DGP and clinical/histological response to a GFD (coeliac disease+IgAd). Patients with seropositive coeliac disease served as controls.

RESULTS

Of 227 patients previously diagnosed with seronegative coeliac disease, true seronegative coeliac disease was confirmed in 84, coeliac disease+IgAd in 48, and excluded in 55. Lack of follow-up duodenal biopsy precluded diagnosing seronegative coeliac disease in 40 patients. 2084 patients with seropositive coeliac disease served as controls. True seronegative coeliac disease had more severe symptoms at diagnosis and a higher risk of complications (HR 10.87, 95% CI 6.11-19.33, P < 0.001) and mortality (HR 2.18, 95% CI 1.12-4.26, P < 0.01) than seropositive coeliac disease. There were no differences between true seronegative coeliac disease and coeliac disease+IgAd. On multivariate analysis, age at diagnosis, lack of clinical response to a GFD, true seronegative coeliac disease, coeliac disease+IgAd, and classical presentation predicted complications. Age at diagnosis, complications and absence of clinical response to a GFD predicted mortality.

CONCLUSIONS

Seronegative coeliac disease has a more aggressive disease phenotype than seropositive coeliac disease. These data argue against over-reliance on serology for the diagnosis of coeliac disease and support a strict clinical and histologic follow-up in seronegative coeliac disease.

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.

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.

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.

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).

What is the Optimal Method Assessing for Persistent Villous Atrophy in Adult Coeliac Disease?

BACKGROUND AND AIMS

Methods of assessing gluten-free diet (GFD) adherence in adults with coeliac disease (CD) include serological testing, dietitian evaluation, questionnaires and repeat duodenal biopsies. Persisting villous atrophy (VA) is associated with CD complications, however gastroscopy with biopsies is expensive and invasive. This study aimed to assess the abilities of a duodenal bulb (D1) biopsy and the Celiac Dietary Adherence Test (CDAT) to detect persisting VA in adults with CD.

METHODS

A prospective observational study of adult CD patients referred for follow-up duodenal biopsies was performed. Quadrantic biopsies were taken from the second part of the duodenum (D2), in addition to a D1 biopsy. Patients underwent follow-up serological testing, and completed the CDAT and Biagi Score. These non-invasive adherence markers were compared against duodenal histology.

RESULTS

368 patients (mean age 51.0 years, 70.1% female) had D1 and D2 biopsies taken at follow-up gastroscopy. Compared to D2 biopsies alone, additional D1 biopsies increased detection of VA by 10.4% (p<0.0001). 201 patients (mean age 50.3 years, 67.7% female) completed adherence questionnaires and serology. When detecting VA, sensitivities and specificities of these markers were 39.7% and 94.2% for IgA- tTG, 38.1% and 96.4% for IgA-EMA, 55.6% and 52.2% for CDAT and 20.6% and 96.4% for the Biagi score.

CONCLUSIONS

Bulbar biopsies increase detection of persisting VA by 10.4%. Serology, CDAT and Biagi performed poorly when predicting VA. The gold standard for predicting persisting VA remains repeat biopsy.

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.

Accuracy of a no-biopsy approach for the diagnosis of coeliac disease across different adult cohorts

OBJECTIVE

We aimed to determine the predictive capacity and diagnostic yield of a 10-fold increase in serum IgA antitissue transglutaminase (tTG) antibody levels for detecting small intestinal injury diagnostic of coeliac disease (CD) in adult patients.

DESIGN

The study comprised three adult cohorts. Cohort 1: 740 patients assessed in the specialist CD clinic at a UK centre; cohort 2: 532 patients with low suspicion for CD referred for upper GI endoscopy at a UK centre; cohort 3: 145 patients with raised tTG titres from multiple international sites. Marsh 3 histology was used as a reference standard against which we determined the performance characteristics of an IgA tTG titre of ≥10×ULN for a diagnosis of CD.

RESULTS

Cohort 1: the sensitivity, specificity, positive predictive value (PPV) and negative predictive value (NPV) for IgA tTG levels of ≥10×ULN at identifying individuals with Marsh 3 lesions were 54.0%, 90.0%, 98.7% and 12.5%, respectively. Cohort 2: the sensitivity, specificity, PPV and NPV for IgA tTG levels of ≥10×ULN at identifying individuals with Marsh 3 lesions were 50.0%, 100.0%, 100.0% and 98.3%, respectively. Cohort 3: the sensitivity, specificity, PPV and NPV for IgA tTG levels of ≥10×ULN at identifying individuals with Marsh 3 lesions were 30.0%, 83.0%, 95.2% and 9.5%, respectively.

CONCLUSION

Our results show that IgA tTG titres of ≥10×ULN have a strong predictive value at identifying adults with intestinal changes diagnostic of CD. This study supports the use of a no-biopsy approach for the diagnosis of adult CD.

Breast Cancer Risk Factors and Survival by Tumor Subtype: Pooled Analyses from the Breast Cancer Association Consortium

BACKGROUND

It is not known whether modifiable lifestyle factors that predict survival after invasive breast cancer differ by subtype.

METHODS

We analyzed data for 121,435 women diagnosed with breast cancer from 67 studies in the Breast Cancer Association Consortium with 16,890 deaths (8,554 breast cancer specific) over 10 years. Cox regression was used to estimate associations between risk factors and 10-year all-cause mortality and breast cancer-specific mortality overall, by estrogen receptor (ER) status, and by intrinsic-like subtype.

RESULTS

There was no evidence of heterogeneous associations between risk factors and mortality by subtype (P adj > 0.30). The strongest associations were between all-cause mortality and BMI ≥30 versus 18.5-25 kg/m2 [HR (95% confidence interval (CI), 1.19 (1.06-1.34)]; current versus never smoking [1.37 (1.27-1.47)], high versus low physical activity [0.43 (0.21-0.86)], age ≥30 years versus <20 years at first pregnancy [0.79 (0.72-0.86)]; >0-<5 years versus ≥10 years since last full-term birth [1.31 (1.11-1.55)]; ever versus never use of oral contraceptives [0.91 (0.87-0.96)]; ever versus never use of menopausal hormone therapy, including current estrogen-progestin therapy [0.61 (0.54-0.69)]. Similar associations with breast cancer mortality were weaker; for example, 1.11 (1.02-1.21) for current versus never smoking.

CONCLUSIONS

We confirm associations between modifiable lifestyle factors and 10-year all-cause mortality. There was no strong evidence that associations differed by ER status or intrinsic-like subtype.

IMPACT

Given the large dataset and lack of evidence that associations between modifiable risk factors and 10-year mortality differed by subtype, these associations could be cautiously used in prognostication models to inform patient-centered care.

Combined Associations of a Polygenic Risk Score and Classical Risk Factors With Breast Cancer Risk

We evaluated the joint associations between a new 313-variant PRS (PRS313) and questionnaire-based breast cancer risk factors for women of European ancestry, using 72 284 cases and 80 354 controls from the Breast Cancer Association Consortium. Interactions were evaluated using standard logistic regression and a newly developed case-only method for breast cancer risk overall and by estrogen receptor status. After accounting for multiple testing, we did not find evidence that per-standard deviation PRS313 odds ratio differed across strata defined by individual risk factors. Goodness-of-fit tests did not reject the assumption of a multiplicative model between PRS313 and each risk factor. Variation in projected absolute lifetime risk of breast cancer associated with classical risk factors was greater for women with higher genetic risk (PRS313 and family history) and, on average, 17.5% higher in the highest vs lowest deciles of genetic risk. These findings have implications for risk prevention for women at increased risk of breast cancer.

CYP3A7*1C allele: linking premenopausal oestrone and progesterone levels with risk of hormone receptor-positive breast cancers

BACKGROUND

Epidemiological studies provide strong evidence for a role of endogenous sex hormones in the aetiology of breast cancer. The aim of this analysis was to identify genetic variants that are associated with urinary sex-hormone levels and breast cancer risk.

METHODS

We carried out a genome-wide association study of urinary oestrone-3-glucuronide and pregnanediol-3-glucuronide levels in 560 premenopausal women, with additional analysis of progesterone levels in 298 premenopausal women. To test for the association with breast cancer risk, we carried out follow-up genotyping in 90,916 cases and 89,893 controls from the Breast Cancer Association Consortium. All women were of European ancestry.

RESULTS

For pregnanediol-3-glucuronide, there were no genome-wide significant associations; for oestrone-3-glucuronide, we identified a single peak mapping to the CYP3A locus, annotated by rs45446698. The minor rs45446698-C allele was associated with lower oestrone-3-glucuronide (-49.2%, 95% CI -56.1% to -41.1%, P = 3.1 × 10-18); in follow-up analyses, rs45446698-C was also associated with lower progesterone (-26.7%, 95% CI -39.4% to -11.6%, P = 0.001) and reduced risk of oestrogen and progesterone receptor-positive breast cancer (OR = 0.86, 95% CI 0.82-0.91, P = 6.9 × 10-8).

CONCLUSIONS

The CYP3A7*1C allele is associated with reduced risk of hormone receptor-positive breast cancer possibly mediated via an effect on the metabolism of endogenous sex hormones in premenopausal women.

Breast Cancer Polygenic Risk Score and Contralateral Breast Cancer Risk

Previous research has shown that polygenic risk scores (PRSs) can be used to stratify women according to their risk of developing primary invasive breast cancer. This study aimed to evaluate the association between a recently validated PRS of 313 germline variants (PRS313) and contralateral breast cancer (CBC) risk. We included 56,068 women of European ancestry diagnosed with first invasive breast cancer from 1990 onward with follow-up from the Breast Cancer Association Consortium. Metachronous CBC risk (N = 1,027) according to the distribution of PRS313 was quantified using Cox regression analyses. We assessed PRS313 interaction with age at first diagnosis, family history, morphology, ER status, PR status, and HER2 status, and (neo)adjuvant therapy. In studies of Asian women, with limited follow-up, CBC risk associated with PRS313 was assessed using logistic regression for 340 women with CBC compared with 12,133 women with unilateral breast cancer. Higher PRS313 was associated with increased CBC risk: hazard ratio per standard deviation (SD) = 1.25 (95%CI = 1.18-1.33) for Europeans, and an OR per SD = 1.15 (95%CI = 1.02-1.29) for Asians. The absolute lifetime risks of CBC, accounting for death as competing risk, were 12.4% for European women at the 10th percentile and 20.5% at the 90th percentile of PRS313. We found no evidence of confounding by or interaction with individual characteristics, characteristics of the primary tumor, or treatment. The C-index for the PRS313 alone was 0.563 (95%CI = 0.547-0.586). In conclusion, PRS313 is an independent factor associated with CBC risk and can be incorporated into CBC risk prediction models to help improve stratification and optimize surveillance and treatment strategies.

Transcriptome-wide association study of breast cancer risk by estrogen-receptor status

Previous transcriptome-wide association studies (TWAS) have identified breast cancer risk genes by integrating data from expression quantitative loci and genome-wide association studies (GWAS), but analyses of breast cancer subtype-specific associations have been limited. In this study, we conducted a TWAS using gene expression data from GTEx and summary statistics from the hitherto largest GWAS meta-analysis conducted for breast cancer overall, and by estrogen receptor subtypes (ER+ and ER-). We further compared associations with ER+ and ER- subtypes, using a case-only TWAS approach. We also conducted multigene conditional analyses in regions with multiple TWAS associations. Two genes, STXBP4 and HIST2H2BA, were specifically associated with ER+ but not with ER- breast cancer. We further identified 30 TWAS-significant genes associated with overall breast cancer risk, including four that were not identified in previous studies. Conditional analyses identified single independent breast-cancer gene in three of six regions harboring multiple TWAS-significant genes. Our study provides new information on breast cancer genetics and biology, particularly about genomic differences between ER+ and ER- breast cancer.

Genome-wide association study identifies 32 novel breast cancer susceptibility loci from overall and subtype-specific analyses

Breast cancer susceptibility variants frequently show heterogeneity in associations by tumor subtype1-3. To identify novel loci, we performed a genome-wide association study including 133,384 breast cancer cases and 113,789 controls, plus 18,908 BRCA1 mutation carriers (9,414 with breast cancer) of European ancestry, using both standard and novel methodologies that account for underlying tumor heterogeneity by estrogen receptor, progesterone receptor and human epidermal growth factor receptor 2 status and tumor grade. We identified 32 novel susceptibility loci (P < 5.0 × 10-8), 15 of which showed evidence for associations with at least one tumor feature (false discovery rate < 0.05). Five loci showed associations (P < 0.05) in opposite directions between luminal and non-luminal subtypes. In silico analyses showed that these five loci contained cell-specific enhancers that differed between normal luminal and basal mammary cells. The genetic correlations between five intrinsic-like subtypes ranged from 0.35 to 0.80. The proportion of genome-wide chip heritability explained by all known susceptibility loci was 54.2% for luminal A-like disease and 37.6% for triple-negative disease. The odds ratios of polygenic risk scores, which included 330 variants, for the highest 1% of quantiles compared with middle quantiles were 5.63 and 3.02 for luminal A-like and triple-negative disease, respectively. These findings provide an improved understanding of genetic predisposition to breast cancer subtypes and will inform the development of subtype-specific polygenic risk scores.

A network analysis to identify mediators of germline-driven differences in breast cancer prognosis

Identifying the underlying genetic drivers of the heritability of breast cancer prognosis remains elusive. We adapt a network-based approach to handle underpowered complex datasets to provide new insights into the potential function of germline variants in breast cancer prognosis. This network-based analysis studies ~7.3 million variants in 84,457 breast cancer patients in relation to breast cancer survival and confirms the results on 12,381 independent patients. Aggregating the prognostic effects of genetic variants across multiple genes, we identify four gene modules associated with survival in estrogen receptor (ER)-negative and one in ER-positive disease. The modules show biological enrichment for cancer-related processes such as G-alpha signaling, circadian clock, angiogenesis, and Rho-GTPases in apoptosis.

Fine-mapping of 150 breast cancer risk regions identifies 191 likely target genes

Genome-wide association studies have identified breast cancer risk variants in over 150 genomic regions, but the mechanisms underlying risk remain largely unknown. These regions were explored by combining association analysis with in silico genomic feature annotations. We defined 205 independent risk-associated signals with the set of credible causal variants in each one. In parallel, we used a Bayesian approach (PAINTOR) that combines genetic association, linkage disequilibrium and enriched genomic features to determine variants with high posterior probabilities of being causal. Potentially causal variants were significantly over-represented in active gene regulatory regions and transcription factor binding sites. We applied our INQUSIT pipeline for prioritizing genes as targets of those potentially causal variants, using gene expression (expression quantitative trait loci), chromatin interaction and functional annotations. Known cancer drivers, transcription factors and genes in the developmental, apoptosis, immune system and DNA integrity checkpoint gene ontology pathways were over-represented among the highest-confidence target genes.

Capsule endoscopy in the management of refractory coeliac disease

BACKGROUND AND AIMS

There is no literature on the role of repeat small bowel capsule endoscopy (SBCE) in patients with refractory coeliac disease (RCD) following treatment with steroids +/- immunosuppressants.

METHODS

The findings on SBCEs from a group of patients with histologically proven RCD (n=23) were compared to the findings from patients with uncomplicated coeliac disease (n=48). All patients had concurrent duodenal histology and serology taken at the time of SBCE.

RESULTS

Patients with RCD had a greater extent of mucosal involvement on SBCE than patients with uncomplicated CD (42.4+/-34.1% vs 9.7+/-21.7%, p=0.0001). Following treatment with steroids and / or immunosuppressants, patients with RCD had an improvement in the extent of affected small bowel mucosa (42.4+/-34.1% vs 26.4+/-28.9% p=0.012). There was no statistical difference in histology and serology taken at the time of the first and second SBCE in patients with RCD.

CONCLUSIONS

Our study suggests that SBCE is valuable in documenting the extent of mucosal involvement in patients with RCD. This is the first study that delineates the value of a second look SBCE to assess improvement in the extent of disease in the small bowel following treatment.

Genome-wide association and transcriptome studies identify target genes and risk loci for breast cancer

Genome-wide association studies (GWAS) have identified more than 170 breast cancer susceptibility loci. Here we hypothesize that some risk-associated variants might act in non-breast tissues, specifically adipose tissue and immune cells from blood and spleen. Using expression quantitative trait loci (eQTL) reported in these tissues, we identify 26 previously unreported, likely target genes of overall breast cancer risk variants, and 17 for estrogen receptor (ER)-negative breast cancer, several with a known immune function. We determine the directional effect of gene expression on disease risk measured based on single and multiple eQTL. In addition, using a gene-based test of association that considers eQTL from multiple tissues, we identify seven (and four) regions with variants associated with overall (and ER-negative) breast cancer risk, which were not reported in previous GWAS. Further investigation of the function of the implicated genes in breast and immune cells may provide insights into the etiology of breast cancer.

Genome-wide association study of germline variants and breast cancer-specific mortality

BACKGROUND

We examined the associations between germline variants and breast cancer mortality using a large meta-analysis of women of European ancestry.

METHODS

Meta-analyses included summary estimates based on Cox models of twelve datasets using ~10.4 million variants for 96,661 women with breast cancer and 7697 events (breast cancer-specific deaths). Oestrogen receptor (ER)-specific analyses were based on 64,171 ER-positive (4116) and 16,172 ER-negative (2125) patients. We evaluated the probability of a signal to be a true positive using the Bayesian false discovery probability (BFDP).

RESULTS

We did not find any variant associated with breast cancer-specific mortality at P < 5 × 10-8. For ER-positive disease, the most significantly associated variant was chr7:rs4717568 (BFDP = 7%, P = 1.28 × 10-7, hazard ratio [HR] = 0.88, 95% confidence interval [CI] = 0.84-0.92); the closest gene is AUTS2. For ER-negative disease, the most significant variant was chr7:rs67918676 (BFDP = 11%, P = 1.38 × 10-7, HR = 1.27, 95% CI = 1.16-1.39); located within a long intergenic non-coding RNA gene (AC004009.3), close to the HOXA gene cluster.

CONCLUSIONS

We uncovered germline variants on chromosome 7 at BFDP < 15% close to genes for which there is biological evidence related to breast cancer outcome. However, the paucity of variants associated with mortality at genome-wide significance underpins the challenge in providing genetic-based individualised prognostic information for breast cancer patients.

Polygenic Risk Scores for Prediction of Breast Cancer and Breast Cancer Subtypes

Stratification of women according to their risk of breast cancer based on polygenic risk scores (PRSs) could improve screening and prevention strategies. Our aim was to develop PRSs, optimized for prediction of estrogen receptor (ER)-specific disease, from the largest available genome-wide association dataset and to empirically validate the PRSs in prospective studies. The development dataset comprised 94,075 case subjects and 75,017 control subjects of European ancestry from 69 studies, divided into training and validation sets. Samples were genotyped using genome-wide arrays, and single-nucleotide polymorphisms (SNPs) were selected by stepwise regression or lasso penalized regression. The best performing PRSs were validated in an independent test set comprising 11,428 case subjects and 18,323 control subjects from 10 prospective studies and 190,040 women from UK Biobank (3,215 incident breast cancers). For the best PRSs (313 SNPs), the odds ratio for overall disease per 1 standard deviation in ten prospective studies was 1.61 (95%CI: 1.57-1.65) with area under receiver-operator curve (AUC) = 0.630 (95%CI: 0.628-0.651). The lifetime risk of overall breast cancer in the top centile of the PRSs was 32.6%. Compared with women in the middle quintile, those in the highest 1% of risk had 4.37- and 2.78-fold risks, and those in the lowest 1% of risk had 0.16- and 0.27-fold risks, of developing ER-positive and ER-negative disease, respectively. Goodness-of-fit tests indicated that this PRS was well calibrated and predicts disease risk accurately in the tails of the distribution. This PRS is a powerful and reliable predictor of breast cancer risk that may improve breast cancer prevention programs.

A transcriptome-wide association study of 229,000 women identifies new candidate susceptibility genes for breast cancer

The breast cancer risk variants identified in genome-wide association studies explain only a small fraction of the familial relative risk, and the genes responsible for these associations remain largely unknown. To identify novel risk loci and likely causal genes, we performed a transcriptome-wide association study evaluating associations of genetically predicted gene expression with breast cancer risk in 122,977 cases and 105,974 controls of European ancestry. We used data from the Genotype-Tissue Expression Project to establish genetic models to predict gene expression in breast tissue and evaluated model performance using data from The Cancer Genome Atlas. Of the 8,597 genes evaluated, significant associations were identified for 48 at a Bonferroni-corrected threshold of P < 5.82 × 10-6, including 14 genes at loci not yet reported for breast cancer. We silenced 13 genes and showed an effect for 11 on cell proliferation and/or colony-forming efficiency. Our study provides new insights into breast cancer genetics and biology.

The BRCA2 c.68-7T > A variant is not pathogenic: A model for clinical calibration of spliceogenicity

Although the spliceogenic nature of the BRCA2 c.68-7T > A variant has been demonstrated, its association with cancer risk remains controversial. In this study, we accurately quantified by real-time PCR and digital PCR (dPCR), the BRCA2 isoforms retaining or missing exon 3. In addition, the combined odds ratio for causality of the variant was estimated using genetic and clinical data, and its associated cancer risk was estimated by case-control analysis in 83,636 individuals. Co-occurrence in trans with pathogenic BRCA2 variants was assessed in 5,382 families. Exon 3 exclusion rate was 4.5-fold higher in variant carriers (13%) than controls (3%), indicating an exclusion rate for the c.68-7T > A allele of approximately 20%. The posterior probability of pathogenicity was 7.44 × 10-115 . There was neither evidence for increased risk of breast cancer (OR 1.03; 95% CI 0.86-1.24) nor for a deleterious effect of the variant when co-occurring with pathogenic variants. Our data provide for the first time robust evidence of the nonpathogenicity of the BRCA2 c.68-7T > A. Genetic and quantitative transcript analyses together inform the threshold for the ratio between functional and altered BRCA2 isoforms compatible with normal cell function. These findings might be exploited to assess the relevance for cancer risk of other BRCA2 spliceogenic variants.

Publisher Correction: Evidence that breast cancer risk at the 2q35 locus is mediated through IGFBP5 regulation

This corrects the article DOI: 10.1038/ncomms5999.

Joint associations of a polygenic risk score and environmental risk factors for breast cancer in the Breast Cancer Association Consortium

Background

Polygenic risk scores (PRS) for breast cancer can be used to stratify the population into groups at substantially different levels of risk. Combining PRS and environmental risk factors will improve risk prediction; however, integrating PRS into risk prediction models requires evaluation of their joint association with known environmental risk factors.

Methods

Analyses were based on data from 20 studies; datasets analysed ranged from 3453 to 23 104 invasive breast cancer cases and similar numbers of controls, depending on the analysed environmental risk factor. We evaluated joint associations of a 77-single nucleotide polymorphism (SNP) PRS with reproductive history, alcohol consumption, menopausal hormone therapy (MHT), height and body mass index (BMI). We tested the null hypothesis of multiplicative joint associations for PRS and each of the environmental factors, and performed global and tail-based goodness-of-fit tests in logistic regression models. The outcomes were breast cancer overall and by estrogen receptor (ER) status.

Results

The strongest evidence for a non-multiplicative joint associations with the 77-SNP PRS was for alcohol consumption (P-interaction = 0.009), adult height (P-interaction = 0.025) and current use of combined MHT (P-interaction = 0.038) in ER-positive disease. Risk associations for these factors by percentiles of PRS did not follow a clear dose-response. In addition, global and tail-based goodness of fit tests showed little evidence for departures from a multiplicative risk model, with alcohol consumption showing the strongest evidence for ER-positive disease (P = 0.013 for global and 0.18 for tail-based tests).

Conclusions

The combined effects of the 77-SNP PRS and environmental risk factors for breast cancer are generally well described by a multiplicative model. Larger studies are required to confirm possible departures from the multiplicative model for individual risk factors, and assess models specific for ER-negative disease.

Outcome of elective withdrawal of anti-tumour necrosis factor-α therapy in patients with Crohn's disease in established remission

BACKGROUND AND AIMS

Outcomes of cessation of anti-TNF therapy for Crohn's disease (CD) in clinical and/or endoscopic remission in routine clinical practice is uncertain. This study aimed to evaluate clinical outcomes and factors associated with relapse in CD patients following formal disease assessment and elective anti-TNF withdrawal.

METHODS

Prospective observational study of CD patients in whom anti-TNF therapy was stopped electively after ≥12months and follow-up of ≥6months. Investigations at assessment prior to cessation included ≥1 of clinical assessment, endoscopic and/or imaging. Relapse was defined as recurrent symptoms of CD requiring medical or surgical therapy.

RESULTS

Eighty-six patients received anti-TNF for a median duration of 23 (12-80) months for severe active luminal (70%), fistulating perianal (25.5%) and other fistulating disease (4.5%). Relapse rates at 90,180 and 365days were 4.7%, 18.6% and 36%, respectively. If anti-TNF dose escalation occurred 6months prior to withdrawal, 88% (7/8) relapsed. Based on multivariate analysis, risk factors for relapse include ileocolonic disease at diagnosis and previous anti-TNF therapy. An elevated faecal calprotectin (FC) is likely to predict relapse (p=0.02), with a PPV of 66.7% at >50μg/g. Of 36 patients who relapsed, 31 were retreated with anti-TNF, with an overall recapture rate of 93%.

CONCLUSION

Relapse rates at 1year following elective withdrawal of anti-TNF are 36%, with high retreatment response rate. Predictors of relapse include ileocolonic involvement, previous anti-TNF therapy and raised FC. Endoscopic/radiologic assessment prior to cessation of therapy does not appear to predict those at lower risk of relapse.

Body mass index and breast cancer survival: a Mendelian randomization analysis

Background

There is increasing evidence that elevated body mass index (BMI) is associated with reduced survival for women with breast cancer. However, the underlying reasons remain unclear. We conducted a Mendelian randomization analysis to investigate a possible causal role of BMI in survival from breast cancer.

Methods

We used individual-level data from six large breast cancer case-cohorts including a total of 36 210 individuals (2475 events) of European ancestry. We created a BMI genetic risk score (GRS) based on genotypes at 94 known BMI-associated genetic variants. Association between the BMI genetic score and breast cancer survival was analysed by Cox regression for each study separately. Study-specific hazard ratios were pooled using fixed-effect meta-analysis.

Results

BMI genetic score was found to be associated with reduced breast cancer-specific survival for estrogen receptor (ER)-positive cases [hazard ratio (HR) = 1.11, per one-unit increment of GRS, 95% confidence interval (CI) 1.01-1.22, P = 0.03). We observed no association for ER-negative cases (HR = 1.00, per one-unit increment of GRS, 95% CI 0.89-1.13, P = 0.95).

Conclusions

Our findings suggest a causal effect of increased BMI on reduced breast cancer survival for ER-positive breast cancer. There is no evidence of a causal effect of higher BMI on survival for ER-negative breast cancer cases.

Identification of ten variants associated with risk of estrogen-receptor-negative breast cancer

Most common breast cancer susceptibility variants have been identified through genome-wide association studies (GWAS) of predominantly estrogen receptor (ER)-positive disease. We conducted a GWAS using 21,468 ER-negative cases and 100,594 controls combined with 18,908 BRCA1 mutation carriers (9,414 with breast cancer), all of European origin. We identified independent associations at P < 5 × 10-8 with ten variants at nine new loci. At P < 0.05, we replicated associations with 10 of 11 variants previously reported in ER-negative disease or BRCA1 mutation carrier GWAS and observed consistent associations with ER-negative disease for 105 susceptibility variants identified by other studies. These 125 variants explain approximately 16% of the familial risk of this breast cancer subtype. There was high genetic correlation (0.72) between risk of ER-negative breast cancer and breast cancer risk for BRCA1 mutation carriers. These findings may lead to improved risk prediction and inform further fine-mapping and functional work to better understand the biological basis of ER-negative breast cancer.

PHIP - a novel candidate breast cancer susceptibility locus on 6q14.1

Most non-BRCA1/2 breast cancer families have no identified genetic cause. We used linkage and haplotype analyses in familial and sporadic breast cancer cases to identify a susceptibility locus on chromosome 6q. Two independent genome-wide linkage analysis studies suggested a 3 Mb locus on chromosome 6q and two unrelated Swedish families with a LOD >2 together seemed to share a haplotype in 6q14.1. We hypothesized that this region harbored a rare high-risk founder allele contributing to breast cancer in these two families. Sequencing of DNA and RNA from the two families did not detect any pathogenic mutations. Finally, 29 SNPs in the region were analyzed in 44,214 cases and 43,532 controls from BCAC, and the original haplotypes in the two families were suggested as low-risk alleles for European and Swedish women specifically. There was also some support for one additional independent moderate-risk allele in Swedish familial samples. The results were consistent with our previous findings in familial breast cancer and supported a breast cancer susceptibility locus at 6q14.1 around the PHIP gene.

Reproductive profiles and risk of breast cancer subtypes: a multi-center case-only study

BACKGROUND

Previous studies have shown that reproductive factors are differentially associated with breast cancer (BC) risk by subtypes. The aim of this study was to investigate associations between reproductive factors and BC subtypes, and whether these vary by age at diagnosis.

METHODS

We used pooled data on tumor markers (estrogen and progesterone receptor, human epidermal growth factor receptor-2 (HER2)) and reproductive risk factors (parity, age at first full-time pregnancy (FFTP) and age at menarche) from 28,095 patients with invasive BC from 34 studies participating in the Breast Cancer Association Consortium (BCAC). In a case-only analysis, we used logistic regression to assess associations between reproductive factors and BC subtype compared to luminal A tumors as a reference. The interaction between age and parity in BC subtype risk was also tested, across all ages and, because age was modeled non-linearly, specifically at ages 35, 55 and 75 years.

RESULTS

Parous women were more likely to be diagnosed with triple negative BC (TNBC) than with luminal A BC, irrespective of age (OR for parity = 1.38, 95% CI 1.16-1.65, p = 0.0004; p for interaction with age = 0.076). Parous women were also more likely to be diagnosed with luminal and non-luminal HER2-like BCs and this effect was slightly more pronounced at an early age (p for interaction with age = 0.037 and 0.030, respectively). For instance, women diagnosed at age 35 were 1.48 (CI 1.01-2.16) more likely to have luminal HER2-like BC than luminal A BC, while this association was not significant at age 75 (OR = 0.72, CI 0.45-1.14). While age at menarche was not significantly associated with BC subtype, increasing age at FFTP was non-linearly associated with TNBC relative to luminal A BC. An age at FFTP of 25 versus 20 years lowered the risk for TNBC (OR = 0.78, CI 0.70-0.88, p < 0.0001), but this effect was not apparent at a later FFTP.

CONCLUSIONS

Our main findings suggest that parity is associated with TNBC across all ages at BC diagnosis, whereas the association with luminal HER2-like BC was present only for early onset BC.

Association analysis identifies 65 new breast cancer risk loci

Breast cancer risk is influenced by rare coding variants in susceptibility genes, such as BRCA1, and many common, mostly non-coding variants. However, much of the genetic contribution to breast cancer risk remains unknown. Here we report the results of a genome-wide association study of breast cancer in 122,977 cases and 105,974 controls of European ancestry and 14,068 cases and 13,104 controls of East Asian ancestry. We identified 65 new loci that are associated with overall breast cancer risk at P < 5 × 10-8. The majority of credible risk single-nucleotide polymorphisms in these loci fall in distal regulatory elements, and by integrating in silico data to predict target genes in breast cells at each locus, we demonstrate a strong overlap between candidate target genes and somatic driver genes in breast tumours. We also find that heritability of breast cancer due to all single-nucleotide polymorphisms in regulatory features was 2-5-fold enriched relative to the genome-wide average, with strong enrichment for particular transcription factor binding sites. These results provide further insight into genetic susceptibility to breast cancer and will improve the use of genetic risk scores for individualized screening and prevention.

The clinical and phenotypical assessment of seronegative villous atrophy; a prospective UK centre experience evaluating 200 adult cases over a 15-year period (2000-2015)

BACKGROUND

Seronegative villous atrophy (SNVA) is commonly attributed to coeliac disease (CD). However, there are other causes of SNVA. More recently angiotensin-2-receptor-blockers (A2RBs) have been reported as an association but data on SNVA have been limited to centres evaluating complex case referrals and not SNVA in general.

OBJECTIVES

To provide clinical outcomes and associations in a large prospective study overseeing all newcomers with SNVA.

DESIGN

Over a 15-year period (2000-2015) we evaluated 200 adult patients with SNVA at a UK centre. A diagnosis of either seronegative CD (SNCD) or seronegative non-CD (SN-non-CD) was reached. Baseline comparisons were made between the groups, with 343 seropositive CD subjects serving as controls.

RESULTS

Of the 200 SNVA cases, SNCD represented 31% (n=62) and SN-non-CD 69% (n=138). The human leucocyte antigen (HLA)-DQ2 and/or DQ8 genotype was present in 61%, with a 51% positive predictive value for SNCD. The breakdown of identifiable causes in the SN-non-CD group comprised infections (27%, n=54), inflammatory/immune-mediated disorders (17.5%, n=35) and drugs (6.5%, n=13; two cases related to A2RBs). However, no cause was found in 18% (n=36) and of these 72% (n=26/36) spontaneously normalised duodenal histology while consuming a gluten-enriched diet. Following multivariable logistic regression analysis an independent factor associated with SN-non-CD was non-white ethnicity (OR 10.8, 95% CI 2.2 to 52.8); in fact, 66% of non-whites had GI infections. On immunohistochemistry all groups stained positive for CD8-T-cytotoxic intraepithelial lymphocytes. However, additional CD4-T helper intraepithelial lymphocytes were occasionally seen in SN-non-CD mimicking the changes associated with refractory CD.

CONCLUSIONS

Most patients with SNVA do not have CD, in particular those who are not white. Furthermore, a subgroup with no obvious aetiology will show spontaneous histological resolution while consuming gluten. These findings suggest caution in empirically prescribing a gluten-free diet without investigation.

Unrestricted faecal calprotectin testing performs poorly in the diagnosis of inflammatory bowel disease in patients in primary care

Background

Faecal calprotectin (FC) measurement distinguishes patients with inflammatory bowel disease (IBD) from those with irritable bowel syndrome but evidence of its performance in primary care is limited.

Aims

To assess the yield of IBD from FC testing in primary care.

Methods

Retrospective review of hospital records to assess the outcome following FC testing in primary care. Investigations for all patients undergoing FC testing in a single laboratory for 6 months from 1 October 2013 to 28 February 2014 were reviewed.

Results

410 patients (162 male; median age 42; range 16–91) were included. FC>50 µg/g was considered positive (FC+). 148/410 (36.1%; median age 44 (17–91)) were FC+ (median FC 116.5 µg/g (51–1770)). 122/148 FC-positive patients (82.4%) underwent further investigation. 97 (65.5%) underwent lower gastrointestinal endoscopy (LGIE), of which 7 (7.2%) had IBD. 49/262 (18.7%) FC-negative (FC−) patients (FC ≤50 µg/g) (median age 47 (19–76)) also underwent LGIE, of whom 3 (6.1%) had IBD.

IBD was diagnosed in 11/410 (2.7%; 4 ulcerative colitis, 3 Crohn’s disease, 4 microscopic colitis). 8/11 were FC+ (range 67–1170) and 3 FC−. At a 50 µg/g threshold, sensitivity for detecting IBD was 72.7%, specificity 64.9%, positive predictive value (PPV) 5.41% and negative predictive value 98.9%. Increasing the threshold to 100 µg/g reduced the sensitivity of the test for detecting IBD to 54.6%.

Conclusions

FC testing in primary care has low sensitivity and specificity with poor PPV for diagnosing IBD. Its use needs to be directed to those with a higher pretest probability of disease. Local services and laboratories should advise general practitioners accordingly.

BRCA2 Hypomorphic Missense Variants Confer Moderate Risks of Breast Cancer

Breast cancer risks conferred by many germline missense variants in the BRCA1 and BRCA2 genes, often referred to as variants of uncertain significance (VUS), have not been established. In this study, associations between 19 BRCA1 and 33 BRCA2 missense substitution variants and breast cancer risk were investigated through a breast cancer case-control study using genotyping data from 38 studies of predominantly European ancestry (41,890 cases and 41,607 controls) and nine studies of Asian ancestry (6,269 cases and 6,624 controls). The BRCA2 c.9104A>C, p.Tyr3035Ser (OR = 2.52; P = 0.04), and BRCA1 c.5096G>A, p.Arg1699Gln (OR = 4.29; P = 0.009) variant were associated with moderately increased risks of breast cancer among Europeans, whereas BRCA2 c.7522G>A, p.Gly2508Ser (OR = 2.68; P = 0.004), and c.8187G>T, p.Lys2729Asn (OR = 1.4; P = 0.004) were associated with moderate and low risks of breast cancer among Asians. Functional characterization of the BRCA2 variants using four quantitative assays showed reduced BRCA2 activity for p.Tyr3035Ser compared with wild-type. Overall, our results show how BRCA2 missense variants that influence protein function can confer clinically relevant, moderately increased risks of breast cancer, with potential implications for risk management guidelines in women with these specific variants. Cancer Res; 77(11); 2789-99. ©2017 AACR.

Genetic modifiers of CHEK2*1100delC-associated breast cancer risk

PURPOSE

CHEK2*1100delC is a founder variant in European populations that confers a two- to threefold increased risk of breast cancer (BC). Epidemiologic and family studies have suggested that the risk associated with CHEK2*1100delC is modified by other genetic factors in a multiplicative fashion. We have investigated this empirically using data from the Breast Cancer Association Consortium (BCAC).

METHODS

Using genotype data from 39,139 (624 1100delC carriers) BC patients and 40,063 (224) healthy controls from 32 BCAC studies, we analyzed the combined risk effects of CHEK2*1100delC and 77 common variants in terms of a polygenic risk score (PRS) and pairwise interaction.

RESULTS

The PRS conferred odds ratios (OR) of 1.59 (95% CI: 1.21-2.09) per standard deviation for BC for CHEK2*1100delC carriers and 1.58 (1.55-1.62) for noncarriers. No evidence of deviation from the multiplicative model was found. The OR for the highest quintile of the PRS was 2.03 (0.86-4.78) for CHEK2*1100delC carriers, placing them in the high risk category according to UK NICE guidelines. The OR for the lowest quintile was 0.52 (0.16-1.74), indicating a lifetime risk close to the population average.

CONCLUSION

Our results confirm the multiplicative nature of risk effects conferred by CHEK2*1100delC and the common susceptibility variants. Furthermore, the PRS could identify carriers at a high lifetime risk for clinical actions.Genet Med advance online publication 06 October 2016.

Association of breast cancer risk with genetic variants showing differential allelic expression: Identification of a novel breast cancer susceptibility locus at 4q21

There are significant inter-individual differences in the levels of gene expression. Through modulation of gene expression, cis-acting variants represent an important source of phenotypic variation. Consequently, cis-regulatory SNPs associated with differential allelic expression are functional candidates for further investigation as disease-causing variants. To investigate whether common variants associated with differential allelic expression were involved in breast cancer susceptibility, a list of genes was established on the basis of their involvement in cancer related pathways and/or mechanisms. Thereafter, using data from a genome-wide map of allelic expression associated SNPs, 313 genetic variants were selected and their association with breast cancer risk was then evaluated in 46,451 breast cancer cases and 42,599 controls of European ancestry ascertained from 41 studies participating in the Breast Cancer Association Consortium. The associations were evaluated with overall breast cancer risk and with estrogen receptor negative and positive disease. One novel breast cancer susceptibility locus on 4q21 (rs11099601) was identified (OR = 1.05, P = 5.6x10-6). rs11099601 lies in a 135 kb linkage disequilibrium block containing several genes, including, HELQ, encoding the protein HEL308 a DNA dependant ATPase and DNA Helicase involved in DNA repair, MRPS18C encoding the Mitochondrial Ribosomal Protein S18C and FAM175A (ABRAXAS), encoding a BRCA1 BRCT domain-interacting protein involved in DNA damage response and double-strand break (DSB) repair. Expression QTL analysis in breast cancer tissue showed rs11099601 to be associated with HELQ (P = 8.28x10-14), MRPS18C (P = 1.94x10-27) and FAM175A (P = 3.83x10-3), explaining about 20%, 14% and 1%, respectively of the variance inexpression of these genes in breast carcinomas.

rs2735383, located at a microRNA binding site in the 3'UTR of NBS1, is not associated with breast cancer risk

NBS1, also known as NBN, plays an important role in maintaining genomic stability. Interestingly, rs2735383 G > C, located in a microRNA binding site in the 3'-untranslated region (UTR) of NBS1, was shown to be associated with increased susceptibility to lung and colorectal cancer. However, the relation between rs2735383 and susceptibility to breast cancer is not yet clear. Therefore, we genotyped rs2735383 in 1,170 familial non-BRCA1/2 breast cancer cases and 1,077 controls using PCR-based restriction fragment length polymorphism (RFLP-PCR) analysis, but found no association between rs2735383CC and breast cancer risk (OR = 1.214, 95% CI = 0.936-1.574, P = 0.144). Because we could not exclude a small effect size due to a limited sample size, we further analyzed imputed rs2735383 genotypes (r2 > 0.999) of 47,640 breast cancer cases and 46,656 controls from the Breast Cancer Association Consortium (BCAC). However, rs2735383CC was not associated with overall breast cancer risk in European (OR = 1.014, 95% CI = 0.969-1.060, P = 0.556) nor in Asian women (OR = 0.998, 95% CI = 0.905-1.100, P = 0.961). Subgroup analyses by age, age at menarche, age at menopause, menopausal status, number of pregnancies, breast feeding, family history and receptor status also did not reveal a significant association. This study therefore does not support the involvement of the genotype at NBS1 rs2735383 in breast cancer susceptibility.