Longitudinal Change in Serum Neurofilament Light Chain in Type 2 Diabetes and Early Diabetic Polyneuropathy: ADDITION-Denmark

OBJECTIVE

To investigate the longitudinal development of neurofilament light chain (NfL) levels in type 2 diabetes with and without diabetic polyneuropathy (+/-DPN) and to explore the predictive potential of NfL as a biomarker for DPN.

RESEARCH DESIGN AND METHODS

We performed retrospective longitudinal case-control analysis of data from 178 participants of the Anglo-Danish-Dutch Study of Intensive Treatment in People with Screen-Detected Diabetes in Primary Care-Denmark (ADDITION-Denmark) cohort of people with screen-detected type 2 diabetes. Biobank samples acquired at the ADDITION-Denmark 5- and 10-year follow-ups were analyzed for serum NfL (s-NfL) using single-molecule array, and the results were compared with established reference material to obtain NfL z-scores. DPN was diagnosed according to Toronto criteria for confirmed DPN at the 10-year follow-up.

RESULTS

s-NfL increased over time in +DPN (N = 39) and -DPN participants (N = 139) at levels above normal age-induced s-NfL increase. Longitudinal s-NfL change was greater in +DPN than in -DPN participants (17.4% [95% CI 4.3; 32.2] or 0.31 SD [95% CI 0.03; 0.60] higher s-NfL or NfL z-score increase in +DPN compared with -DPN). s-NfL at the 5-year follow-up was positively associated with nerve conduction studies at the 10-year follow-up (P = 0.02 to <0.001), but not with DPN risk. Areas under the curve (AUCs) for s-NfL were not inferior to AUCs for the Michigan Neuropathy Screening Instrument questionnaire score or vibration detection thresholds. Higher yearly s-NfL increase was associated with higher DPN risk (odds ratio 1.36 [95% CI 1.08; 1.71] per 1 ng/L/year).

CONCLUSIONS

Our findings suggest that preceding s-NfL trajectories differ slightly between those with and without DPN and imply a possible biomarker value of s-NfL trajectories in DPN.

All exons are not created equal-exon vulnerability determines the effect of exonic mutations on splicing

It is now widely accepted that aberrant splicing of constitutive exons is often caused by mutations affecting cis-acting splicing regulatory elements (SREs), but there is a misconception that all exons have an equal dependency on SREs and thus a similar vulnerability to aberrant splicing. We demonstrate that some exons are more likely to be affected by exonic splicing mutations (ESMs) due to an inherent vulnerability, which is context dependent and influenced by the strength of exon definition. We have developed VulExMap, a tool which is based on empirical data that can designate whether a constitutive exon is vulnerable. Using VulExMap, we find that only 25% of all exons can be categorized as vulnerable, whereas two-thirds of 359 previously reported ESMs in 75 disease genes are located in vulnerable exons. Because VulExMap analysis is based on empirical data on splicing of exons in their endogenous context, it includes all features important in determining the vulnerability. We believe that VulExMap will be an important tool when assessing the effect of exonic mutations by pinpointing whether they are located in exons vulnerable to ESMs.

Molecular portraits of patients with intrahepatic cholangiocarcinoma who diverge as rapid progressors or long survivors on chemotherapy

OBJECTIVE

Cytotoxic agents are the cornerstone of treatment for patients with advanced intrahepatic cholangiocarcinoma (iCCA), despite heterogeneous benefit. We hypothesised that the pretreatment molecular profiles of diagnostic biopsies can predict patient benefit from chemotherapy and define molecular bases of innate chemoresistance.

DESIGN

We identified a cohort of advanced iCCA patients with comparable baseline characteristics who diverged as extreme outliers on chemotherapy (survival <6 m in rapid progressors, RP; survival >23 m in long survivors, LS). Diagnostic biopsies were characterised by digital pathology, then subjected to whole-transcriptome profiling of bulk and geospatially macrodissected tissue regions. Spatial transcriptomics of tumour-infiltrating myeloid cells was performed using targeted digital spatial profiling (GeoMx). Transcriptome signatures were evaluated in multiple cohorts of resected cancers. Signatures were also characterised using in vitro cell lines, in vivo mouse models and single cell RNA-sequencing data.

RESULTS

Pretreatment transcriptome profiles differentiated patients who would become RPs or LSs on chemotherapy. Biologically, this signature originated from altered tumour-myeloid dynamics, implicating tumour-induced immune tolerogenicity with poor response to chemotherapy. The central role of the liver microenviroment was confrmed by the association of the RPLS transcriptome signature with clinical outcome in iCCA but not extrahepatic CCA, and in liver metastasis from colorectal cancer, but not in the matched primary bowel tumours.

CONCLUSIONS

The RPLS signature could be a novel metric of chemotherapy outcome in iCCA. Further development and validation of this transcriptomic signature is warranted to develop precision chemotherapy strategies in these settings.

Use of sugar in coffee and tea and long-term risk of mortality in older adult Danish men: 32 years of follow-up from a prospective cohort study

BACKGROUND

Tea and coffee are the most consumed beverages worldwide and very often sweetened with sugar. However, the association between the use of sugar in tea or coffee and adverse events is currently unclear.

OBJECTIVES

To investigate the association between the addition of sugar to coffee or tea, and the risk of all-cause mortality, cardiovascular mortality, cancer mortality and incident diabetes mellitus.

METHODS

Participants from the prospective Copenhagen Male Study, included from 1985 to 1986, without cardiovascular disease, cancer or diabetes mellitus at inclusion, who reported regular coffee or tea consumption were included. Self-reported number of cups of coffee and tea and use of sugar were derived from the study questionnaires. Quantity of sugar use was not reported. Primary outcome was all-cause mortality and secondary endpoints were cardiovascular mortality, cancer mortality and incident diabetes mellitus, all assessed through the Danish national registries. The association between adding sugar and all-cause mortality was analyzed by Cox regression analysis. Age, smoking status, daily alcohol intake, systolic blood pressure, body mass index, number of cups of coffee and/or tea consumed per day and socioeconomic status were included as covariates. Vital status of patients up and until 22.03.2017 was assessed. Sugar could be added to either coffee, tea or both.

RESULTS

In total, 2923 men (mean age at inclusion: 63±5 years) were included, of which 1007 (34.5%) added sugar. In 32 years of follow-up, 2581 participants (88.3%) died, 1677 in the non-sugar group (87.5%) versus 904 in the sugar group (89.9%). Hazard ratio of the sugar group compared to the non-sugar group was 1.06 (95% CI 0.98;1.16) for all-cause mortality. An interaction term between number of cups of coffee and/or tea per day and adding sugar was 0.99 (0.96;1.01). A subgroup analysis of coffee-only drinkers showed a hazard ratio of 1.11 (0.99;1.26). The interaction term was 0.98 (0.94;1.02). Hazard ratios for the sugar group compared to the non-sugar group were 1.11 (95% CI 0.97;1.26) for cardiovascular disease mortality, 1.01 (95% CI 0.87;1.17) for cancer mortality and 1.04 (95% CI 0.79;1.36) for incident diabetes mellitus.

CONCLUSION

In the present population of Danish men, use of sugar in tea and/or coffee was not significantly associated with increased risk of mortality or incident diabetes.

Hormonal, immune, and oxidative stress responses to blood flow-restricted exercise

INTRODUCTION

Heavy-load free-flow resistance exercise (HL-FFRE) is a widely used training modality. Recently, low-load blood-flow restricted resistance exercise (LL-BFRRE) has gained attention in both athletic and clinical settings as an alternative when conventional HL-FFRE is contraindicated or not tolerated. LL-BFRRE has been shown to result in physiological adaptations in muscle and connective tissue that are comparable to those induced by HL-FFRE. The underlying mechanisms remain unclear; however, evidence suggests that LL-BFRRE involves elevated metabolic stress compared to conventional free-flow resistance exercise (FFRE).

AIM

The aim was to evaluate the initial (<10 min post-exercise), intermediate (10-20 min), and late (>30 min) hormonal, immune, and oxidative stress responses observed following acute sessions of LL-BFRRE compared to FFRE in healthy adults.

METHODS

A systematic literature search of randomized and non-randomized studies was conducted in PubMed, Embase, Cochrane Central, CINAHL, and SPORTDiscus. The Cochrane Risk of Bias (RoB2, ROBINS-1) and TESTEX were used to evaluate risk of bias and study quality. Data extractions were based on mean change within groups.

RESULTS

A total of 12525 hits were identified, of which 29 articles were included. LL-BFRRE demonstrated greater acute increases in growth hormone responses when compared to overall FFRE at intermediate (SMD 2.04; 95% CI 0.87, 3.22) and late (SMD 2.64; 95% CI 1.13, 4.16) post-exercise phases. LL-BFRRE also demonstrated greater increase in testosterone responses compared to late LL-FFRE.

CONCLUSION

These results indicate that LL-BFRRE can induce increased or similar hormone and immune responses compared to LL-FFRE and HL-FFRE along with attenuated oxidative stress responses compared to HL-FFRE.

The impact of AUD on death for men with different IQ-scores: a register-based cohort study of 645 955 men

AIM

To investigate whether the effect of alcohol use disorder (AUD) on death by natural and unnatural causes, respectively, differs according to intelligence quotient (IQ) scores.

METHODS

We followed 654 955 Danish men, including 75 267 brothers, born between 1939 and 1959 from their 25th birthday, 1 January 1970, or date of conscription (whichever came last) until 31 December 2018. The exposure of AUD was defined by first registered treatment (diagnosis since 1969, prescription medicine since 1994, or other treatment since 2006), and the outcomes of death by natural and unnatural causes, respectively, were obtained from nationwide registers since 1970. Information on IQ score was retrieved at conscription from the Danish Conscription Database.

RESULTS AND CONCLUSION

In total, 86 106 men were defined with an AUD. AUD combined with the highest, middle, and lowest IQ score tertiles, respectively, were associated with a 5.90 (95% confidence interval [CI] 5.75; 6.01), 6.88 (95% CI: 6.73; 7.04), and 7.53 (95% CI: 7.38; 7.68) times higher hazard of death by natural causes compared with no AUD and the highest IQ score tertile. The risk of death by unnatural causes was comparable for men with AUD regardless of IQ score tertile. A within-brother analysis showed that the impact of AUD on death by natural and unnatural causes, respectively, did not vary between men with different IQ score tertiles, but were hampered by statistical uncertainty. Our study indicates a need of special focus on men with lower levels of IQ score and AUD for prevention of death by natural causes.

Prenatal and postnatal exposure to acetaminophen in relation to autism spectrum and attention-deficit and hyperactivity symptoms in childhood: Meta-analysis in six European population-based cohorts

The potential etiological role of early acetaminophen exposure on Autism Spectrum Conditions (ASC) and Attention-Deficit/Hyperactivity Disorder (ADHD) is inconclusive. We aimed to study this association in a collaborative study of six European population-based birth/child cohorts. A total of 73,881 mother-child pairs were included in the study. Prenatal and postnatal (up to 18 months) acetaminophen exposure was assessed through maternal questionnaires or interviews. ASC and ADHD symptoms were assessed at 4-12 years of age using validated instruments. Children were classified as having borderline/clinical symptoms using recommended cutoffs for each instrument. Hospital diagnoses were also available in one cohort. Analyses were adjusted for child and maternal characteristics along with indications for acetaminophen use. Adjusted cohort-specific effect estimates were combined using random-effects meta-analysis. The proportion of children having borderline/clinical symptoms ranged between 0.9 and 12.9% for ASC and between 1.2 and 12.2% for ADHD. Results indicated that children prenatally exposed to acetaminophen were 19% and 21% more likely to subsequently have borderline or clinical ASC (OR = 1.19, 95% CI 1.07-1.33) and ADHD symptoms (OR = 1.21, 95% CI 1.07-1.36) compared to non-exposed children. Boys and girls showed higher odds for ASC and ADHD symptoms after prenatal exposure, though these associations were slightly stronger among boys. Postnatal exposure to acetaminophen was not associated with ASC or ADHD symptoms. These results replicate previous work and support providing clear information to pregnant women and their partners about potential long-term risks of acetaminophen use.

Genome-Wide Association Analysis of Pancreatic Beta-Cell Glucose Sensitivity

CONTEXT

Pancreatic beta-cell glucose sensitivity is the slope of the plasma glucose-insulin secretion relationship and is a key predictor of deteriorating glucose tolerance and development of type 2 diabetes. However, there are no large-scale studies looking at the genetic determinants of beta-cell glucose sensitivity.

OBJECTIVE

To understand the genetic determinants of pancreatic beta-cell glucose sensitivity using genome-wide meta-analysis and candidate gene studies.

DESIGN

We performed a genome-wide meta-analysis for beta-cell glucose sensitivity in subjects with type 2 diabetes and nondiabetic subjects from 6 independent cohorts (n = 5706). Beta-cell glucose sensitivity was calculated from mixed meal and oral glucose tolerance tests, and its associations between known glycemia-related single nucleotide polymorphisms (SNPs) and genome-wide association study (GWAS) SNPs were estimated using linear regression models.

RESULTS

Beta-cell glucose sensitivity was moderately heritable (h2 ranged from 34% to 55%) using SNP and family-based analyses. GWAS meta-analysis identified multiple correlated SNPs in the CDKAL1 gene and GIPR-QPCTL gene loci that reached genome-wide significance, with SNP rs2238691 in GIPR-QPCTL (P value = 2.64 × 10-9) and rs9368219 in the CDKAL1 (P value = 3.15 × 10-9) showing the strongest association with beta-cell glucose sensitivity. These loci surpassed genome-wide significance when the GWAS meta-analysis was repeated after exclusion of the diabetic subjects. After correction for multiple testing, glycemia-associated SNPs in or near the HHEX and IGF2B2 loci were also associated with beta-cell glucose sensitivity.

CONCLUSION

We show that, variation at the GIPR-QPCTL and CDKAL1 loci are key determinants of pancreatic beta-cell glucose sensitivity.

Genetic deletion of connexin 37 causes polyuria and polydipsia

The connexin 37 (Cx37) channel is clustered at gap junctions between cells in the renal vasculature or the renal tubule where it is abundant in basolateral cell interdigitations and infoldings of epithelial cells in the proximal tubule, thick ascending limb, distal convoluted tubule and collecting duct; however, physiological data regarding its role are limited. In this study, we investigated the role of Cx37 in fluid homeostasis using mice with a global deletion of Cx37 (Cx37-/- mice). Under baseline conditions, Cx37-/- had ~40% higher fluid intake associated with ~40% lower urine osmolality compared to wild-type (WT) mice. No differences were observed between genotypes in urinary adenosine triphosphate or prostaglandin E2, paracrine factors that alter renal water handling. After 18-hours of water deprivation, plasma aldosterone and urine osmolality increased significantly in Cx37-/- and WT mice; however, the latter remained ~375 mmol/kg lower in Cx37-/- mice, an effect associated with a more pronounced body weight loss despite higher urinary AVP/creatinine ratios compared to WT mice. Consistent with this, fluid intake in the first 3 hours after water deprivation was 37% greater in Cx37-/- vs WT mice. Cx37-/- mice showed significantly lower renal AQP2 abundance and AQP2 phosphorylation at serine 256 than WT mice in response to vehicle or dDAVP, suggesting a partial contribution of the kidney to the lower urine osmolality. The abundance and responses of the vasopressin V2 receptor, AQP3, NHE3, NKCC2, NCC, H+-ATPase, αENaC, γENaC or Na+/K+-ATPase were not significantly different between genotypes. In summary, these results demonstrate that Cx37 is important for body water handling.

Patient-reported outcomes after 10-year follow-up of intensive, multifactorial treatment in individuals with screen-detected type 2 diabetes: the ADDITION-Europe trial

AIMS

To present the longer-term impact of multifactorial treatment of type 2 diabetes on self-reported health status, diabetes-specific quality of life, and diabetes treatment satisfaction at 10-year follow up of the ADDITION-Europe trial.

METHODS

The ADDITION-Europe trial enrolled 3057 individuals with screen-detected type 2 diabetes from four centres [Denmark, the UK (Cambridge and Leicester) and the Netherlands], between 2001 and 2006. Participants were randomized at general practice level to intensive treatment or to routine care . The trial ended in 2009 and a 10-year follow-up was performed at the end of 2014. We measured self-reported health status (36-item Short-Form Health Survey and EQ-5D), diabetes-specific quality of life (Audit of Diabetes-Dependent Quality of Life questionnaire), and diabetes treatment satisfaction (Diabetes Treatment Satisfaction Questionnaire) at different time points during the study period. A mixed-effects model was applied to estimate the effect of intensive treatment (intention-to-treat analyses) on patient-reported outcome measures for each centre. Centre-specific estimates were pooled using a fixed effects meta-analysis.

RESULTS

There was no difference in patient-reported outcome measures between the routine care and intensive treatment arms in this 10-year follow-up study [EQ-5D: -0.01 (95% CI -0.03, 0.01); Physical Composite Score (36-item Short-Form Health Survey): -0.27 (95% CI -1.11, 0.57), Audit of Diabetes-Dependent Quality of Life questionnaire: -0.01 (95% CI -0.11, 0.10); and Diabetes Treatment Satisfaction Questionnaire: -0.20 (95% CI -0.70, 0.29)].

CONCLUSIONS

Intensive, multifactorial treatment of individuals with screen-detected type 2 diabetes did not affect self-reported health status, diabetes-specific quality of life, or diabetes treatment satisfaction at 10-year follow-up compared to routine care.

Pan-cancer analysis of whole genomes identifies driver rearrangements promoted by LINE-1 retrotransposition

About half of all cancers have somatic integrations of retrotransposons. Here, to characterize their role in oncogenesis, we analyzed the patterns and mechanisms of somatic retrotransposition in 2,954 cancer genomes from 38 histological cancer subtypes within the framework of the Pan-Cancer Analysis of Whole Genomes (PCAWG) project. We identified 19,166 somatically acquired retrotransposition events, which affected 35% of samples and spanned a range of event types. Long interspersed nuclear element (LINE-1; L1 hereafter) insertions emerged as the first most frequent type of somatic structural variation in esophageal adenocarcinoma, and the second most frequent in head-and-neck and colorectal cancers. Aberrant L1 integrations can delete megabase-scale regions of a chromosome, which sometimes leads to the removal of tumor-suppressor genes, and can induce complex translocations and large-scale duplications. Somatic retrotranspositions can also initiate breakage-fusion-bridge cycles, leading to high-level amplification of oncogenes. These observations illuminate a relevant role of L1 retrotransposition in remodeling the cancer genome, with potential implications for the development of human tumors.

TRAIP drives replisome disassembly and mitotic DNA repair synthesis at sites of incomplete DNA replication

The faithful segregation of eukaryotic chromosomes in mitosis requires that the genome be duplicated completely prior to anaphase. However, cells with large genomes sometimes fail to complete replication during interphase and instead enter mitosis with regions of incompletely replicated DNA. These regions are processed in early mitosis via a process known as mitotic DNA repair synthesis (MiDAS), but little is known about how cells switch from conventional DNA replication to MiDAS. Using the early embryo of the nematode Caenorhabditis elegans as a model system, we show that the TRAIP ubiquitin ligase drives replisome disassembly in response to incomplete DNA replication, thereby providing access to replication forks for other factors. Moreover, TRAIP is essential for MiDAS in human cells, and is important in both systems to prevent mitotic segregation errors. Our data indicate that TRAIP is a master regulator of the processing of incomplete DNA replication during mitosis in metazoa.

GWAS on longitudinal growth traits reveals different genetic factors influencing infant, child, and adult BMI

Early childhood growth patterns are associated with adult health, yet the genetic factors and the developmental stages involved are not fully understood. Here, we combine genome-wide association studies with modeling of longitudinal growth traits to study the genetics of infant and child growth, followed by functional, pathway, genetic correlation, risk score, and colocalization analyses to determine how developmental timings, molecular pathways, and genetic determinants of these traits overlap with those of adult health. We found a robust overlap between the genetics of child and adult body mass index (BMI), with variants associated with adult BMI acting as early as 4 to 6 years old. However, we demonstrated a completely distinct genetic makeup for peak BMI during infancy, influenced by variation at the LEPR/LEPROT locus. These findings suggest that different genetic factors control infant and child BMI. In light of the obesity epidemic, these findings are important to inform the timing and targets of prevention strategies.

Perinatal Whole Blood Zinc Status and Cytokines, Adipokines, and Other Immune Response Proteins

(1) Background: Zinc is an essential micronutrient and zinc deficiency is associated with immune dysfunction. The neonatal immune system is immature, and therefore an optimal neonatal zinc status may be important. The aim of this study was to investigate the possible association between neonatal whole blood (WB)-Zinc content and several immune markers. (2) Methods: In total, 398 healthy newborns (199 who later developed type 1 diabetes and 199 controls) from the Danish Newborn Screening Biobank had neonatal dried blood spots (NDBS) analyzed for WB-Zinc content and (i) cytokines: Interleukin (IL)-1β, IL-4, IL-6, IL-8, IL-10, IL-12 (p70), interferon gamma, tumor necrosis factor alpha, and transforming growth factor beta; (ii) adipokines: leptin and adiponectin; (iii) other immune response proteins: C-reactive protein (CRP), and mannose-binding lectin (MBL), and soluble triggering receptors expressed on myeloid cells1 (sTREM-1). WB-Zinc content was determined using laser ablation inductively coupled plasma mass spectrometry. For each analyte, the relative change in mean level was modelled by a robust log-normal model regression. (3) Results: No association was found between WB-Zinc content and all the immune response markers in either the unadjusted or adjusted models overall or when stratifying by case status. (4) Conclusions: In healthy Danish neonates, WB-Zinc content was not associated with cytokines, adipokines, CRP, MBL or sTREM, which does not indicate a strong immunological function of neonatal zinc status.

Patient survival and tumor characteristics associated with CHEK2:p.I157T - findings from the Breast Cancer Association Consortium

BACKGROUND

P.I157T is a CHEK2 missense mutation associated with a modest increase in breast cancer risk. Previously, another CHEK2 mutation, the protein truncating c.1100delC has been associated with poor prognosis of breast cancer patients. Here, we have investigated patient survival and characteristics of breast tumors of germ line p.I157T carriers.

METHODS

We included in the analyses 26,801 European female breast cancer patients from 15 studies participating in the Breast Cancer Association Consortium. We analyzed the association between p.I157T and the clinico-pathological breast cancer characteristics by comparing the p.I157T carrier tumors to non-carrier and c.1100delC carrier tumors. Similarly, we investigated the p.I157T associated risk of early death, breast cancer-associated death, distant metastasis, locoregional relapse and second breast cancer using Cox proportional hazards models. Additionally, we explored the p.I157T-associated genomic gene expression profile using data from breast tumors of 183 Finnish female breast cancer patients (ten p.I157T carriers) (GEO: GSE24450). Differential gene expression analysis was performed using a moderated t test. Functional enrichment was investigated using the DAVID functional annotation tool and gene set enrichment analysis (GSEA). The tumors were classified into molecular subtypes according to the St Gallen 2013 criteria and the PAM50 gene expression signature.

RESULTS

P.I157T was not associated with increased risk of early death, breast cancer-associated death or distant metastasis relapse, and there was a significant difference in prognosis associated with the two CHEK2 mutations, p.I157T and c.1100delC. Furthermore, p.I157T was associated with lobular histological type and clinico-pathological markers of good prognosis, such as ER and PR expression, low TP53 expression and low grade. Gene expression analysis suggested luminal A to be the most common subtype for p.I157T carriers and CDH1 (cadherin 1) target genes to be significantly enriched among genes, whose expression differed between p.I157T and non-carrier tumors.

CONCLUSIONS

Our analyses suggest that there are fundamental differences in breast tumors of CHEK2:p.I157T and c.1100delC carriers. The poor prognosis associated with c.1100delC cannot be generalized to other CHEK2 mutations.

Identification of independent association signals and putative functional variants for breast cancer risk through fine-scale mapping of the 12p11 locus

BACKGROUND

Multiple recent genome-wide association studies (GWAS) have identified a single nucleotide polymorphism (SNP), rs10771399, at 12p11 that is associated with breast cancer risk.

METHOD

We performed a fine-scale mapping study of a 700 kb region including 441 genotyped and more than 1300 imputed genetic variants in 48,155 cases and 43,612 controls of European descent, 6269 cases and 6624 controls of East Asian descent and 1116 cases and 932 controls of African descent in the Breast Cancer Association Consortium (BCAC; http://bcac.ccge.medschl.cam.ac.uk/ ), and in 15,252 BRCA1 mutation carriers in the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). Stepwise regression analyses were performed to identify independent association signals. Data from the Encyclopedia of DNA Elements project (ENCODE) and the Cancer Genome Atlas (TCGA) were used for functional annotation.

RESULTS

Analysis of data from European descendants found evidence for four independent association signals at 12p11, represented by rs7297051 (odds ratio (OR) = 1.09, 95 % confidence interval (CI) = 1.06-1.12; P = 3 × 10(-9)), rs805510 (OR = 1.08, 95 % CI = 1.04-1.12, P = 2 × 10(-5)), and rs1871152 (OR = 1.04, 95 % CI = 1.02-1.06; P = 2 × 10(-4)) identified in the general populations, and rs113824616 (P = 7 × 10(-5)) identified in the meta-analysis of BCAC ER-negative cases and BRCA1 mutation carriers. SNPs rs7297051, rs805510 and rs113824616 were also associated with breast cancer risk at P < 0.05 in East Asians, but none of the associations were statistically significant in African descendants. Multiple candidate functional variants are located in putative enhancer sequences. Chromatin interaction data suggested that PTHLH was the likely target gene of these enhancers. Of the six variants with the strongest evidence of potential functionality, rs11049453 was statistically significantly associated with the expression of PTHLH and its nearby gene CCDC91 at P < 0.05.

CONCLUSION

This study identified four independent association signals at 12p11 and revealed potentially functional variants, providing additional insights into the underlying biological mechanism(s) for the association observed between variants at 12p11 and breast cancer risk.