Search

Bleeding in early pregnancy and postpartum hemorrhage (PPH) bear substantial risks, with the former closely associated with pregnancy loss and the latter being the foremost cause of maternal death, underscoring the severe impact on maternal-fetal health. We identified five genetic loci linked to PPH in a meta-analysis. Functional annotation analysis indicated candidate genes HAND2, TBX3 and RAP2C/FRMD7 at three loci and showed that at each locus, associated variants were located within binding sites for progesterone receptors. There were strong genetic correlations with birth weight, gestational duration and uterine fibroids. Bleeding in early pregnancy yielded no genome-wide association signals but showed strong genetic correlation with various human traits, suggesting a potentially complex, polygenic etiology. Our results suggest that PPH is related to progesterone signaling dysregulation, whereas early bleeding is a complex trait associated with underlying health and possibly socioeconomic status and may include genetic factors that have not yet been identified., (© 2024. The Author(s).)

Gene promoter and enhancer sequences are bound by transcription factors and are depleted of methylated CpG sites (cytosines preceding guanines in DNA). The absence of methylated CpGs in these sequences typically correlates with increased gene expression, indicating a regulatory role for methylation. We used nanopore sequencing to determine haplotype-specific methylation rates of 15.3 million CpG units in 7,179 whole-blood genomes. We identified 189,178 methylation depleted sequences where three or more proximal CpGs were unmethylated on at least one haplotype. A total of 77,789 methylation depleted sequences (~41%) associated with 80,503 cis-acting sequence variants, which we termed allele-specific methylation quantitative trait loci (ASM-QTLs). RNA sequencing of 896 samples from the same blood draws used to perform nanopore sequencing showed that the ASM-QTL, that is, DNA sequence variability, drives most of the correlation found between gene expression and CpG methylation. ASM-QTLs were enriched 40.2-fold (95% confidence interval 32.2, 49.9) among sequence variants associating with hematological traits, demonstrating that ASM-QTLs are important functional units in the noncoding genome., (© 2024. The Author(s).)

Essential tremor (ET) is a prevalent neurological disorder with a largely unknown underlying biology. In this genome-wide association study meta-analysis, comprising 16,480 ET cases and 1,936,173 controls from seven datasets, we identify 12 sequence variants at 11 loci. Evaluating mRNA expression, splicing, plasma protein levels, and coding effects, we highlight seven putative causal genes at these loci, including CA3 and CPLX1. CA3 encodes Carbonic Anhydrase III and carbonic anhydrase inhibitors have been shown to decrease tremors. CPLX1, encoding Complexin-1, regulates neurotransmitter release. Through gene-set enrichment analysis, we identify a significant association with specific cell types, including dopaminergic and GABAergic neurons, as well as biological processes like Rho GTPase signaling. Genetic correlation analyses reveals a positive association between ET and Parkinson's disease, depression, and anxiety-related phenotypes. This research uncovers risk loci, enhancing our knowledge of the complex genetics of this common but poorly understood disorder, and highlights CA3 and CPLX1 as potential therapeutic targets., (© 2024. The Author(s).)

Background: Long-read sequencing can enable the detection of base modifications, such as CpG methylation, in single molecules of DNA. The most commonly used methods for long-read sequencing are nanopore developed by Oxford Nanopore Technologies (ONT) and single molecule real-time (SMRT) sequencing developed by Pacific Bioscience (PacBio). In this study, we systematically compare the performance of CpG methylation detection from long-read sequencing., Results: We demonstrate that CpG methylation detection from 7179 nanopore-sequenced DNA samples is highly accurate and consistent with 132 oxidative bisulfite-sequenced (oxBS) samples, isolated from the same blood draws. We introduce quality filters for CpGs that further enhance the accuracy of CpG methylation detection from nanopore-sequenced DNA, while removing at most 30% of CpGs. We evaluate the per-site performance of CpG methylation detection across different genomic features and CpG methylation rates and demonstrate how the latest R10.4 flowcell chemistry and base-calling algorithms improve methylation detection from nanopore sequencing. Additionally, we show how the methylation detection of 50 SMRT-sequenced genomes compares to nanopore sequencing and oxBS., Conclusions: This study provides the first systematic comparison of CpG methylation detection tools for long-read sequencing methods. We compare two commonly used computational methods for the detection of CpG methylation in a large number of nanopore genomes, including samples sequenced using the latest R10.4 nanopore flowcell chemistry and 50 SMRT sequenced samples. We provide insights into the strengths and limitations of each sequencing method as well as recommendations for standardization and evaluation of tools designed for genome-scale modified base detection using long-read sequencing., (© 2024. The Author(s).)

Importance: Recurrent pericarditis is a treatment challenge and often a debilitating condition. Drugs inhibiting interleukin 1 cytokines are a promising new treatment option, but their use is based on scarce biological evidence and clinical trials of modest sizes, and the contributions of innate and adaptive immune processes to the pathophysiology are incompletely understood., Objective: To use human genomics, transcriptomics, and proteomics to shed light on the pathogenesis of pericarditis., Design, Setting, and Participants: This was a meta-analysis of genome-wide association studies of pericarditis from 5 countries. Associations were examined between the pericarditis-associated variants and pericarditis subtypes (including recurrent pericarditis) and secondary phenotypes. To explore mechanisms, associations with messenger RNA expression (cis-eQTL), plasma protein levels (pQTL), and CpG methylation of DNA (ASM-QTL) were assessed. Data from Iceland (deCODE genetics, 1983-2020), Denmark (Copenhagen Hospital Biobank/Danish Blood Donor Study, 1977-2022), the UK (UK Biobank, 1953-2021), the US (Intermountain, 1996-2022), and Finland (FinnGen, 1970-2022) were included. Data were analyzed from September 2022 to August 2023., Exposure: Genotype., Main Outcomes and Measures: Pericarditis., Results: In this genome-wide association study of 4894 individuals with pericarditis (mean [SD] age at diagnosis, 51.4 [17.9] years, 2734 [67.6%] male, excluding the FinnGen cohort), associations were identified with 2 independent common intergenic variants at the interleukin 1 locus on chromosome 2q14. The lead variant was rs12992780 (T) (effect allele frequency [EAF], 31%-40%; odds ratio [OR], 0.83; 95% CI, 0.79-0.87; P = 6.67 × 10-16), downstream of IL1B and the secondary variant rs7575402 (A or T) (EAF, 45%-55%; adjusted OR, 0.89; 95% CI, 0.85-0.93; adjusted P = 9.6 × 10-8). The lead variant rs12992780 had a smaller odds ratio for recurrent pericarditis (0.76) than the acute form (0.86) (P for heterogeneity = .03) and rs7575402 was associated with CpG methylation overlapping binding sites of 4 transcription factors known to regulate interleukin 1 production: PU.1 (encoded by SPI1), STAT1, STAT3, and CCAAT/enhancer-binding protein β (encoded by CEBPB)., Conclusions and Relevance: This study found an association between pericarditis and 2 independent sequence variants at the interleukin 1 gene locus. This finding has the potential to contribute to development of more targeted and personalized therapy of pericarditis with interleukin 1-blocking drugs.

Migraine is a complex neurovascular disease with a range of severity and symptoms, yet mostly studied as one phenotype in genome-wide association studies (GWAS). Here we combine large GWAS datasets from six European populations to study the main migraine subtypes, migraine with aura (MA) and migraine without aura (MO). We identified four new MA-associated variants (in PRRT2, PALMD, ABO and LRRK2) and classified 13 MO-associated variants. Rare variants with large effects highlight three genes. A rare frameshift variant in brain-expressed PRRT2 confers large risk of MA and epilepsy, but not MO. A burden test of rare loss-of-function variants in SCN11A, encoding a neuron-expressed sodium channel with a key role in pain sensation, shows strong protection against migraine. Finally, a rare variant with cis-regulatory effects on KCNK5 confers large protection against migraine and brain aneurysms. Our findings offer new insights with therapeutic potential into the complex biology of migraine and its subtypes., (© 2023. The Author(s).)

Microsatellites are polymorphic tracts of short tandem repeats with one to six base-pair (bp) motifs and are some of the most polymorphic variants in the genome. Using 6084 Icelandic parent-offspring trios we estimate 63.7 (95% CI: 61.9-65.4) microsatellite de novo mutations (mDNMs) per offspring per generation, excluding one bp repeats motifs (homopolymers) the estimate is 48.2 mDNMs (95% CI: 46.7-49.6). Paternal mDNMs occur at longer repeats than maternal ones, which are in turn larger with a mean size of 3.4 bp vs 3.1 bp for paternal ones. mDNMs increase by 0.97 (95% CI: 0.90-1.04) and 0.31 (95% CI: 0.25-0.37) per year of father's and mother's age at conception, respectively. Here, we find two independent coding variants that associate with the number of mDNMs transmitted to offspring; The minor allele of a missense variant (allele frequency (AF) = 1.9%) in MSH2, a mismatch repair gene, increases transmitted mDNMs from both parents (effect: 13.1 paternal and 7.8 maternal mDNMs). A synonymous variant (AF = 20.3%) in NEIL2, a DNA damage repair gene, increases paternally transmitted mDNMs (effect: 4.4 mDNMs). Thus, the microsatellite mutation rate in humans is in part under genetic control., (© 2023. The Author(s).)

Genotypes causing pregnancy loss and perinatal mortality are depleted among living individuals and are therefore difficult to find. To explore genetic causes of recessive lethality, we searched for sequence variants with deficit of homozygosity among 1.52 million individuals from six European populations. In this study, we identified 25 genes harboring protein-altering sequence variants with a strong deficit of homozygosity (10% or less of predicted homozygotes). Sequence variants in 12 of the genes cause Mendelian disease under a recessive mode of inheritance, two under a dominant mode, but variants in the remaining 11 have not been reported to cause disease. Sequence variants with a strong deficit of homozygosity are over-represented among genes essential for growth of human cell lines and genes orthologous to mouse genes known to affect viability. The function of these genes gives insight into the genetics of intrauterine lethality. We also identified 1077 genes with homozygous predicted loss-of-function genotypes not previously described, bringing the total set of genes completely knocked out in humans to 4785., (© 2023. The Author(s).)

Objectives: Erosive hand osteoarthritis (EHOA) is a severe subset of hand osteoarthritis (OA). It is unclear if EHOA is genetically different from other forms of OA. Sequence variants at ten loci have been associated with hand OA but none with EHOA., Methods: We performed meta-analysis of EHOA in 1484 cases and 550 680 controls, from 5 populations. To identify causal genes, we performed eQTL and plasma pQTL analyses, and developed one zebrafish mutant. We analysed associations of variants with other traits and estimated shared genetics between EHOA and other traits., Results: Four common sequence variants associated with EHOA, all with relatively high effect. Rs17013495 ( SPP1 / MEPE , OR=1.40, p=8.4×10 -14 ) and rs11243284 (6p24.3, OR=1.35, p=4.2×10 -11 ) have not been associated with OA, whereas rs11631127 ( ALDH1A2 , OR=1.46, p=7.1×10 -18 ), and rs1800801 ( MGP , OR=1.37, p=3.6×10 -13 ) have previously been associated with hand OA. The association of rs1800801 ( MGP ) was consistent with a recessive mode of inheritance in contrast to its additive association with hand OA (OR homozygotes vs non-carriers=2.01, 95% CI 1.71 to 2.37). All four variants associated nominally with finger OA, although with substantially lower effect. We found shared genetic components between EHOA and other OA measures, grip strength, urate levels and gout, but not rheumatoid arthritis. We identified ALDH1A2 , MGP and BMP6 as causal genes for EHOA, with loss-of-function Bmp6 zebrafish mutants displaying EHOA-like phenotypes., Conclusions: We report on significant genetic associations with EHOA. The results support the view of EHOA as a form of severe hand OA and partly separate it from OA in larger joints., Competing Interests: Competing interests: The authors US, LS, GT, OAS, SS, SHL, DFG, UT and KS are affiliated with deCODE genetics/Amgen and declare competing interests as employees. The author MK reports consultancy/lecture fees outside the submitted work from Pfizer, Novartis, UCB, Galapagos, Flexion, Kiniksa, Jansen, Abbvie, CHDR, all paid to institution. Royalties from Wolters Kluwer and Springer Verlag, paid to institution., (© Author(s) (or their employer(s)) 2023. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.)

Aims: Syncope is a common and clinically challenging condition. In this study, the genetics of syncope were investigated to seek knowledge about its pathophysiology and prognostic implications., Methods and Results: This genome-wide association meta-analysis included 56 071 syncope cases and 890 790 controls from deCODE genetics (Iceland), UK Biobank (United Kingdom), and Copenhagen Hospital Biobank Cardiovascular Study/Danish Blood Donor Study (Denmark), with a follow-up assessment of variants in 22 412 cases and 286 003 controls from Intermountain (Utah, USA) and FinnGen (Finland). The study yielded 18 independent syncope variants, 17 of which were novel. One of the variants, p.Ser140Thr in PTPRN2, affected syncope only when maternally inherited. Another variant associated with a vasovagal reaction during blood donation and five others with heart rate and/or blood pressure regulation, with variable directions of effects. None of the 18 associations could be attributed to cardiovascular or other disorders. Annotation with regard to regulatory elements indicated that the syncope variants were preferentially located in neural-specific regulatory regions. Mendelian randomization analysis supported a causal effect of coronary artery disease on syncope. A polygenic score (PGS) for syncope captured genetic correlation with cardiovascular disorders, diabetes, depression, and shortened lifespan. However, a score based solely on the 18 syncope variants performed similarly to the PGS in detecting syncope risk but did not associate with other disorders., Conclusion: The results demonstrate that syncope has a distinct genetic architecture that implicates neural regulatory processes and a complex relationship with heart rate and blood pressure regulation. A shared genetic background with poor cardiovascular health was observed, supporting the importance of a thorough assessment of individuals presenting with syncope., Competing Interests: Conflicts of interest: The following authors are employees of deCODE genetics/Amgen, Inc.: Hildur M Aegisdottir, Rosa B Thorolfsdottir, Gardar Sveinbjornsson, Olafur A Stefansson, Bjarni Gunnarsson, Vinicius Tragante, Lilja Stefansdottir, Thorgeir E Thorgeirsson, Egil Ferkingstad, Gudmar Thorleifsson, Michael L Frigge, Kristjan E Hjorleifsson, Erna V Ivarsdottir, Anna Helgadottir, Solveig Gretarsdottir, Valgerdur Steinthorsdottir, Asmundur Oddsson, Hannes P Eggertsson, Gisli H Halldorsson, Patrick Sulem, Gudmundur Norddahl, Gudrun Rutsdottir, Gudmundur Thorgeirsson, David O Arnar, Unnur Thorsteinsdottir, Daniel F Gudbjartsson, Hilma Holm, and Kari Stefansson. Karina Banasik and Søren Brunak acknowledge the Novo Nordisk Foundation (grants NNF17OC0027594 and NNF14CC0001)., (© The Author(s) 2023. Published by Oxford University Press on behalf of the European Society of Cardiology. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Detailed knowledge of how diversity in the sequence of the human genome affects phenotypic diversity depends on a comprehensive and reliable characterization of both sequences and phenotypic variation. Over the past decade, insights into this relationship have been obtained from whole-exome sequencing or whole-genome sequencing of large cohorts with rich phenotypic data 1,2 . Here we describe the analysis of whole-genome sequencing of 150,119 individuals from the UK Biobank 3 . This constitutes a set of high-quality variants, including 585,040,410 single-nucleotide polymorphisms, representing 7.0% of all possible human single-nucleotide polymorphisms, and 58,707,036 indels. This large set of variants allows us to characterize selection based on sequence variation within a population through a depletion rank score of windows along the genome. Depletion rank analysis shows that coding exons represent a small fraction of regions in the genome subject to strong sequence conservation. We define three cohorts within the UK Biobank: a large British Irish cohort, a smaller African cohort and a South Asian cohort. A haplotype reference panel is provided that allows reliable imputation of most variants carried by three or more sequenced individuals. We identified 895,055 structural variants and 2,536,688 microsatellites, groups of variants typically excluded from large-scale whole-genome sequencing studies. Using this formidable new resource, we provide several examples of trait associations for rare variants with large effects not found previously through studies based on whole-exome sequencing and/or imputation., (© 2022. The Author(s).)

Long-read sequencing (LRS) promises to improve the characterization of structural variants (SVs). We generated LRS data from 3,622 Icelanders and identified a median of 22,636 SVs per individual (a median of 13,353 insertions and 9,474 deletions). We discovered a set of 133,886 reliably genotyped SV alleles and imputed them into 166,281 individuals to explore their effects on diseases and other traits. We discovered an association of a rare deletion in PCSK9 with lower low-density lipoprotein (LDL) cholesterol levels, compared to the population average. We also discovered an association of a multiallelic SV in ACAN with height; we found 11 alleles that differed in the number of a 57-bp-motif repeat and observed a linear relationship between the number of repeats carried and height. These results show that SVs can be accurately characterized at the population scale using LRS data in a genome-wide non-targeted approach and demonstrate how SVs impact phenotypes.

Despite the important role that monozygotic twins have played in genetics research, little is known about their genomic differences. Here we show that monozygotic twins differ on average by 5.2 early developmental mutations and that approximately 15% of monozygotic twins have a substantial number of these early developmental mutations specific to one of them. Using the parents and offspring of twins, we identified pre-twinning mutations. We observed instances where a twin was formed from a single cell lineage in the pre-twinning cell mass and instances where a twin was formed from several cell lineages. CpG>TpG mutations increased in frequency with embryonic development, coinciding with an increase in DNA methylation. Our results indicate that allocations of cells during development shapes genomic differences between monozygotic twins.

Preeclampsia is a serious complication of pregnancy, affecting both maternal and fetal health. In genome-wide association meta-analysis of European and Central Asian mothers, we identify sequence variants that associate with preeclampsia in the maternal genome at ZNF831/20q13 and FTO/16q12. These are previously established variants for blood pressure (BP) and the FTO variant has also been associated with body mass index (BMI). Further analysis of BP variants establishes that variants at MECOM/3q26, FGF5/4q21 and SH2B3/12q24 also associate with preeclampsia through the maternal genome. We further show that a polygenic risk score for hypertension associates with preeclampsia. However, comparison with gestational hypertension indicates that additional factors modify the risk of preeclampsia.

Pelvic organ prolapse (POP) is a downward descent of one or more of the pelvic organs, resulting in a protrusion of the vaginal wall and/or uterus. We performed a genome-wide association study of POP using data from Iceland and the UK Biobank, a total of 15,010 cases with hospital-based diagnosis code and 340,734 female controls, and found eight sequence variants at seven loci associating with POP (P < 5 × 10 -8 ); seven common (minor allele frequency >5%) and one with minor allele frequency of 4.87%. Some of the variants associating with POP also associated with traits of similar pathophysiology. Of these, rs3820282, which may alter the estrogen-based regulation of WNT4, also associates with leiomyoma of uterus, gestational duration and endometriosis. Rs3791675 at EFEMP1, a gene involved in connective tissue homeostasis, also associates with hernias and carpal tunnel syndrome. Our results highlight the role of connective tissue metabolism and estrogen exposure in the etiology of POP.

Asthma is one of the most common chronic diseases affecting both children and adults. We report a genome-wide association meta-analysis of 69,189 cases and 702,199 controls from Iceland and UK biobank. We find 88 asthma risk variants at 56 loci, 19 previously unreported, and evaluate their effect on other asthma and allergic phenotypes. Of special interest are two low frequency variants associated with protection against asthma; a missense variant in TNFRSF8 and 3' UTR variant in TGFBR1. Functional studies show that the TNFRSF8 variant reduces TNFRSF8 expression both on cell surface and in soluble form, acting as loss of function. eQTL analysis suggests that the TGFBR1 variant acts through gain of function and together with an intronic variant in a downstream gene, SMAD3, points to defective TGFβR1 signaling as one of the biological perturbations increasing asthma risk. Our results increase the number of asthma variants and implicate genes with known role in T cell regulation, inflammation and airway remodeling in asthma pathogenesis.

Background: Breast Cancer 1 gene ( BRCA1 ) is known to be inactivated in breast tumors by promoter methylation. Tumor cells in patients carrying a germline mutation in BRCA1 are sensitive to cytotoxic drugs that cause DNA double strand breaks. However, very little is known on whether patients with BRCA1 promoter methylated tumors are similarly sensitive to cytotoxic drugs. In this study, we address this by making use of extensive follow-up data on patients treated with cyclophosphamide, methotrexate, and fluorouracil in Iceland between 1976 and 2007., Methods: We analyzed BRCA1 promoter methylation by pyrosequencing DNA from tumor samples from 1031 patients with primary breast cancer. Of those, 965 were sporadic cases, 61 were BRCA2 , and five were BRCA1 germline mutation carriers. All cases were examined with respect to clinicopathological parameters and breast cancer-specific survival in patients treated with cytotoxic drugs. Information on chemotherapy treatment in noncarriers was available for 26 BRCA1 methylated tumors and 857 unmethylated tumors., Results: BRCA1 was promoter methylated in 29 sporadic tumors or in 3.0% of cases (29 of 965), whereas none of the tumors derived from BRCA germline mutation carriers were promoter methylated. Important to note, patients with BRCA1 promoter methylation receiving chemotherapeutic drug treatment show highly improved breast cancer-specific survival compared with unmethylated controls (hazard ratio = 0.10, 95% confidence interval = 0.01 to 0.75, two-sided P  = .02)., Conclusions: BRCA1 promoter methylation is predictive of improved disease outcome in patients receiving cyclophosphamide, methotrexate, and fluorouracil drug treatment. Our results support the use of markers indicative of "BRCAness" in sporadic breast cancers to identify patients that are likely to benefit from the use of DNA-damaging agents., (© The Author(s) 2020. Published by Oxford University Press.)

Bone area is one measure of bone size that is easily derived from dual-energy X-ray absorptiometry (DXA) scans. In a GWA study of DXA bone area of the hip and lumbar spine (N ≥ 28,954), we find thirteen independent association signals at twelve loci that replicate in samples of European and East Asian descent (N = 13,608 - 21,277). Eight DXA area loci associate with osteoarthritis, including rs143384 in GDF5 and a missense variant in COL11A1 (rs3753841). The strongest DXA area association is with rs11614913[T] in the microRNA MIR196A2 gene that associates with lumbar spine area (P = 2.3 × 10 -42 , β = -0.090) and confers risk of hip fracture (P = 1.0 × 10 -8 , OR = 1.11). We demonstrate that the risk allele is less efficient in repressing miR-196a-5p target genes. We also show that the DXA area measure contributes to the risk of hip fracture independent of bone density.

In the Methods section of this Article, 'greater than' should have been 'less than' in the sentence 'Putative regions of clustered rearrangements were identified as having an average inter-rearrangement distance that was at least 10 times greater than the whole-genome average for the individual sample. '. The Article has not been corrected.

Nasal polyps (NP) are lesions on the nasal and paranasal sinus mucosa and are a risk factor for chronic rhinosinusitis (CRS). We performed genome-wide association studies on NP and CRS in Iceland and the UK (using UK Biobank data) with 4,366 NP cases, 5,608 CRS cases, and >700,000 controls. We found 10 markers associated with NP and 2 with CRS. We also tested 210 markers reported to associate with eosinophil count, yielding 17 additional NP associations. Of the 27 NP signals, 7 associate with CRS and 13 with asthma. Most notably, a missense variant in ALOX15 that causes a p.Thr560Met alteration in arachidonate 15-lipoxygenase (15-LO) confers large genome-wide significant protection against NP (P = 8.0 × 10 -27 , odds ratio = 0.32; 95% confidence interval = 0.26, 0.39) and CRS (P = 1.1 × 10 -8 , odds ratio = 0.64; 95% confidence interval = 0.55, 0.75). p.Thr560Met, carried by around 1 in 20 Europeans, was previously shown to cause near total loss of 15-LO enzymatic activity. Our findings identify 15-LO as a potential target for therapeutic intervention in NP and CRS.

Genetic diversity arises from recombination and de novo mutation (DNM). Using a combination of microarray genotype and whole-genome sequence data on parent-child pairs, we identified 4,531,535 crossover recombinations and 200,435 DNMs. The resulting genetic map has a resolution of 682 base pairs. Crossovers exhibit a mutagenic effect, with overrepresentation of DNMs within 1 kilobase of crossovers in males and females. In females, a higher mutation rate is observed up to 40 kilobases from crossovers, particularly for complex crossovers, which increase with maternal age. We identified 35 loci associated with the recombination rate or the location of crossovers, demonstrating extensive genetic control of meiotic recombination, and our results highlight genes linked to the formation of the synaptonemal complex as determinants of crossovers., (Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Gallstones are responsible for one of the most common diseases in the Western world and are commonly treated with cholecystectomy. We perform a meta-analysis of two genome-wide association studies of gallstone disease in Iceland and the UK, totaling 27,174 cases and 736,838 controls, uncovering 21 novel gallstone-associated variants at 20 loci. Two distinct low frequency missense variants in SLC10A2, encoding the apical sodium-dependent bile acid transporter (ASBT), associate with an increased risk of gallstone disease (Pro290Ser: OR = 1.36 [1.25-1.49], P = 2.1 × 10 -12 , MAF = 1%; Val98Ile: OR = 1.15 [1.10-1.20], P = 1.8 × 10 -10 , MAF = 4%). We demonstrate that lower bile acid transport by ASBT is accompanied by greater risk of gallstone disease and highlight the role of the intestinal compartment of the enterohepatic circulation of bile acids in gallstone disease susceptibility. Additionally, two low frequency missense variants in SERPINA1 and HNF4A and 17 common variants represent novel associations with gallstone disease.

Benign prostatic hyperplasia and associated lower urinary tract symptoms (BPH/LUTS) are common conditions affecting the majority of elderly males. Here we report the results of a genome-wide association study of symptomatic BPH/LUTS in 20,621 patients and 280,541 controls of European ancestry, from Iceland and the UK. We discovered 23 genome-wide significant variants, located at 14 loci. There is little or no overlap between the BPH/LUTS variants and published prostate cancer risk variants. However, 15 of the variants reported here also associate with serum levels of prostate specific antigen (PSA) (at a Bonferroni corrected P < 0.0022). Furthermore, there is a strong genetic correlation, r g  = 0.77 (P = 2.6 × 10 -11 ), between PSA and BPH/LUTS, and one standard deviation increase in a polygenic risk score (PRS) for BPH/LUTS increases PSA levels by 12.9% (P = 1.6×10 -55 ). These results shed a light on the genetic background of BPH/LUTS and its substantial influence on PSA levels.

Uterine leiomyomas are common benign tumors of the myometrium. We performed a meta-analysis of two genome-wide association studies of leiomyoma in European women (16,595 cases and 523,330 controls), uncovering 21 variants at 16 loci that associate with the disease. Five variants were previously reported to confer risk of various malignant or benign tumors (rs78378222 in TP53, rs10069690 in TERT, rs1800057 and rs1801516 in ATM, and rs7907606 at OBFC1) and four signals are located at established risk loci for hormone-related traits (endometriosis and breast cancer) at 1q36.12 (CDC42/WNT4), 2p25.1 (GREB1), 20p12.3 (MCM8), and 6q26.2 (SYNE1/ESR1). Polygenic score for leiomyoma, computed using UKB data, is significantly correlated with risk of cancer in the Icelandic population. Functional annotation suggests that the non-coding risk variants affect multiple genes, including ESR1. Our results provide insights into the genetic background of leiomyoma that are shared by other benign and malignant tumors and highlight the role of hormones in leiomyoma growth.

Aortic valve stenosis (AS) is the most common valvular heart disease, and valve replacement is the only definitive treatment. Here we report a large genome-wide association (GWA) study of 2,457 Icelandic AS cases and 349,342 controls with a follow-up in up to 4,850 cases and 451,731 controls of European ancestry. We identify two new AS loci, on chromosome 1p21 near PALMD (rs7543130; odds ratio (OR) = 1.20, P = 1.2 × 10 -22 ) and on chromosome 2q22 in TEX41 (rs1830321; OR = 1.15, P = 1.8 × 10 -13 ). Rs7543130 also associates with bicuspid aortic valve (BAV) (OR = 1.28, P = 6.6 × 10 -10 ) and aortic root diameter (P = 1.30 × 10 -8 ), and rs1830321 associates with BAV (OR = 1.12, P = 5.3 × 10 -3 ) and coronary artery disease (OR = 1.05, P = 9.3 × 10 -5 ). The results implicate both cardiac developmental abnormalities and atherosclerosis-like processes in the pathogenesis of AS. We show that several pathways are shared by CAD and AS. Causal analysis suggests that the shared risk factors of Lp(a) and non-high-density lipoprotein cholesterol contribute substantially to the frequent co-occurence of these diseases.

Background: Germline BRCA2 mutations increase risk of breast cancer and other malignancies. BRCA2 has been shown to play a role in telomere protection and maintenance. Telomere length (TL) has been studied as a modifying factor for various diseases, including breast cancer. Previous research on TL in BRCA mutation carriers has produced contradicting results. Methods: We measured blood TL, using a high-throughput monochrome multiplex qPCR method, in a well-defined Icelandic cohort of female BRCA2 mutation carriers ( n = 169), sporadic breast cancer patients ( n = 561), and healthy controls ( n = 537). Results: Breast cancer cases had significantly shorter TL than unaffected women ( P < 0.0001), both BRCA2 mutation carriers ( P = 0.0097) and noncarriers ( P = 0.00006). Using exclusively samples acquired before breast cancer diagnosis, we found that shorter telomeres were significantly associated with increased breast cancer risk in BRCA2 mutation carriers [HR, 3.60; 95% confidence interval (CI), 1.17-11.28; P , 0.025] but not in non-carriers (HR,1.40; 95% CI, 0.89-2.22; P , 0.15). We found no association between TL and breast cancer-specific survival. Conclusions: Blood TL is predictive of breast cancer risk in BRCA2 mutation carriers. Breast cancer cases have significantly shorter TL than unaffected women, regardless of BRCA2 status, indicating that samples taken after breast cancer diagnosis should not be included in evaluations of TL and breast cancer risk. Impact: Our study is built on a well-defined cohort, highly accurate methods, and long follow-up and can therefore help to clarify some previously published, contradictory results. Our findings also suggest that BRCA2 has an important role in telomere maintenance, even in normal blood cells. Cancer Epidemiol Biomarkers Prev; 26(8); 1248-54. ©2017 AACR ., (©2017 American Association for Cancer Research.)

Background: DNA repair of alkylation damage is defective in various cancers. This occurs through somatically acquired inactivation of the MGMT gene in various cancer types, including breast cancers. In addition to MGMT, the two E. coli AlkB homologs ALKBH2 and ALKBH3 have also been linked to direct reversal of alkylation damage. However, it is currently unknown whether ALKBH2 or ALKBH3 are found inactivated in cancer., Methods: Methylome datasets (GSE52865, GSE20713, GSE69914), available through Omnibus, were used to determine whether ALKBH2 or ALKBH3 are found inactivated by CpG promoter methylation. TCGA dataset enabled us to then assess the impact of CpG promoter methylation on mRNA expression for both ALKBH2 and ALKBH3. DNA methylation analysis for the ALKBH3 promoter region was carried out by pyrosequencing (PyroMark Q24) in 265 primary breast tumours and 30 proximal normal breast tissue samples along with 8 breast-derived cell lines. ALKBH3 mRNA and protein expression were analysed in cell lines using RT-PCR and Western blotting, respectively. DNA alkylation damage assay was carried out in cell lines based on immunofluorescence and confocal imaging. Data on clinical parameters and survival outcomes in patients were obtained and assessed in relation to ALKBH3 promoter methylation., Results: The ALKBH3 gene, but not ALKBH2, undergoes CpG promoter methylation and transcriptional silencing in breast cancer. We developed a quantitative alkylation DNA damage assay based on immunofluorescence and confocal imaging revealing higher levels of alkylation damage in association with epigenetic inactivation of the ALKBH3 gene (P = 0.029). In our cohort of 265 primary breast cancer, we found 72 cases showing aberrantly high CpG promoter methylation over the ALKBH3 promoter (27%; 72 out of 265). We further show that increasingly higher degree of ALKBH3 promoter methylation is associated with reduced breast-cancer specific survival times in patients. In this analysis, ALKBH3 promoter methylation at >20% CpG methylation was found to be statistically significantly associated with reduced survival (HR = 2.3; P = 0.012). By thresholding at the clinically relevant CpG methylation level (>20%), we find the incidence of ALKBH3 promoter methylation to be 5% (13 out of 265)., Conclusions: ALKBH3 is a novel addition to the catalogue of DNA repair genes found inactivated in breast cancer. Our results underscore a link between defective alkylation repair and breast cancer which, additionally, is found in association with poor disease outcome.

Obesity is a high risk factor for breast cancer. This relationship could be marked by a specific methylome. The current work was aimed to explore the impact of obesity and menopausal status on variation in breast cancer methylomes. Data from Infinium 450K array-based methylomes of 64 breast tumors were coupled with information on BMI and menopausal status. Additionally, DNA methylation results were validated in 18 non-tumor and 81 tumor breast samples. Breast tumors arising in either pre- or postmenopausal women stratified by BMI or menopausal status alone were not associated with a specific DNA methylation pattern. Intriguingly, the DNA methylation pattern identified in association with the high-risk group (postmenopausal women with high BMI (>25) and premenopausal women with normal or low BMI < 25) exclusively characterized by hypermethylation of 1287 CpG sites as compared with the low-risk group. These CpG sites included the promoter region of fourteen protein-coding genes of which CpG methylation over the ZNF577 promoter region represents the top scoring associated event. In an independent cohort, the ZNF577 promoter methylation remained statistically significant in association with the high-risk group. Additionally, the impact of ZNF577 promoter methylation on mRNA expression levels was demonstrated in breast cancer cell lines after treatment with a demethylating agent (5-azacytidine). In conclusion, the epigenome of breast tumors is affected by a complex interaction between BMI and menopausal status. The ZNF577 methylation quantification is clearly relevant for the development of novel biomarkers of precision therapy in breast cancer., (© 2017 Society for Endocrinology.)

Background: The impact of an inherited BRCA2 mutation on the prognosis of women with breast cancer has not been well documented. We studied the effects of oestrogen receptor (ER) status, other prognostic factors and treatments on survival in a large cohort of BRCA2 mutation carriers., Methods: We identified 285 breast cancer patients with a 999del5 BRCA2 mutation and matched them with 570 non-carrier patients. Clinical information was abstracted from patient charts and pathology records and supplemented by evaluation of tumour grade and ER status using archived tissue specimens. Univariate and multivariate hazard ratios (HR) were estimated for breast cancer-specific survival using Cox regression. The effects of various therapies were studied in patients treated from 1980 to 2012., Results: Among mutation carriers, positive ER status was associated with higher risk of death than negative ER status (HR=1.94; 95% CI=1.22-3.07, P=0.005). The reverse association was seen for non-carriers (HR=0.71; 95% CI: 0.51-0.97; P=0.03)., Conclusions: Among BRCA2 carriers, ER-positive status is an adverse prognostic factor. BRCA2 carrier status should be known at the time when treatment decisions are made.

We analysed whole-genome sequences of 560 breast cancers to advance understanding of the driver mutations conferring clonal advantage and the mutational processes generating somatic mutations. We found that 93 protein-coding cancer genes carried probable driver mutations. Some non-coding regions exhibited high mutation frequencies, but most have distinctive structural features probably causing elevated mutation rates and do not contain driver mutations. Mutational signature analysis was extended to genome rearrangements and revealed twelve base substitution and six rearrangement signatures. Three rearrangement signatures, characterized by tandem duplications or deletions, appear associated with defective homologous-recombination-based DNA repair: one with deficient BRCA1 function, another with deficient BRCA1 or BRCA2 function, the cause of the third is unknown. This analysis of all classes of somatic mutation across exons, introns and intergenic regions highlights the repertoire of cancer genes and mutational processes operating, and progresses towards a comprehensive account of the somatic genetic basis of breast cancer.

In cancer, epigenetic states are deregulated and thought to be of significance in cancer development and progression. We explored DNA methylation-based signatures in association with breast cancer subtypes to assess their impact on clinical presentation and patient prognosis. DNA methylation was analyzed using Infinium 450K arrays in 40 tumors and 17 normal breast samples, together with DNA copy number changes and subtype-specific markers by tissue microarrays. The identified methylation signatures were validated against a cohort of 212 tumors annotated for breast cancer subtypes by the PAM50 method (The Cancer Genome Atlas). Selected markers were pyrosequenced in an independent validation cohort of 310 tumors and analyzed with respect to survival, clinical stage and grade. The results demonstrate that DNA methylation patterns linked to the luminal-B subtype are characterized by CpG island promoter methylation events. In contrast, a large fraction of basal-like tumors are characterized by hypomethylation events occurring within the gene body. Based on these hallmark signatures, we defined two DNA methylation-based subtypes, Epi-LumB and Epi-Basal, and show that they are associated with unfavorable clinical parameters and reduced survival. Our data show that distinct mechanisms leading to changes in CpG methylation states are operative in different breast cancer subtypes. Importantly, we show that a few selected proxy markers can be used to detect the distinct DNA methylation-based subtypes thereby providing valuable information on disease prognosis., (Copyright © 2014 The Authors. Published by Elsevier B.V. All rights reserved.)

Overexpression of the Aurora A kinase has been shown to have prognostic value in breast cancer. Previously, we showed a significant association between AURKA gene amplification and BRCA2 mutation in breast cancer. The aim of this study was to assess the prognostic impact of Aurora A overexpression on breast cancer arising in BRCA2 mutation carriers. Aurora A expression was evaluated by immunohistochemistry on breast tumour tissue microarrays from 107 BRCA2 999del5 mutation carriers and 284 of sporadic origin. Prognostic value of Aurora A nuclear staining was estimated in relation to clinical markers and adjuvant treatment, using multivariate Cox's proportional hazards ratio regression model. BRCA2 wild-type allele loss was measured by TaqMan in BRCA2 mutated tumour samples. All statistical tests were two sided. Multivariate analysis of breast cancer-specific survival, including proliferative markers and treatment, indicated independent prognostic value of Aurora A nuclear staining for BRCA2 mutation carriers (hazards ratio = 7.06; 95% confidence interval = 1.23-40.6; p = 0.028). Poor breast cancer-specific survival of BRCA2 mutation carriers was found to be significantly associated with combined Aurora A nuclear expression and BRCA2 wild type allele loss in tumours (p < 0.001). Multivariate analysis indicated independent prognostic value of both positive Aurora A nuclear staining (hazards ratio = 10.09; 95% confidence interval = 1.19-85.4, p = 0.034) and BRCA2 wild type allele loss (hazards ratio = 9.63; 95% confidence interval = 1.81-51.0, p = 0.008) for BRCA2 mutation carriers. Aurora A nuclear expression was found to be a significant prognostic marker for BRCA2 mutation carriers, independent of clinical parameters and adjuvant treatment. Our conclusion is that treatment benefits for BRCA2 mutation carriers and sporadic breast cancer patients with Aurora A positive tumours may be enhanced by giving attention to Aurora A targeted treatment.

Epithelial-to-mesenchymal transition (EMT) is a plastic process in which fully differentiated epithelial cells are converted into poorly differentiated, migratory and invasive mesenchymal cells, and it has been related to the metastasis potential of tumors. This is a reversible process and cells can also eventually undergo mesenchymal-to-epithelial transition. The existence of a dynamic EMT process suggests the involvement of epigenetic shifts in the phenotype. Herein, we obtained the DNA methylomes at single-base resolution of Madin-Darby canine kidney cells undergoing EMT and translated the identified differentially methylated regions to human breast cancer cells undergoing a gain of migratory and invasive capabilities associated with the EMT phenotype. We noticed dynamic and reversible changes of DNA methylation, both on promoter sequences and gene-bodies in association with transcription regulation of EMT-related genes. Most importantly, the identified DNA methylation markers of EMT were present in primary mammary tumors in association with the epithelial or the mesenchymal phenotype of the studied breast cancer samples., (©2014 American Association for Cancer Research.)

In a recent article, Wang and colleagues reported the discovery of a mechanism by which CARM1 regulates the genomic localization of BAF155 (a SWI/SNF subunit involved in chromatin remodeling) through post-translational methylation at R1064 arginine residues. This modification leads to the relocalization of BAF155-containing SWI/SNF complexes to regions containing genes involved in the Myc oncogenic pathway. The results presented are evidence that these interactions constitute a mechanism by which the BAF155 chromatin remodeling factor contributes to cancer.

Epigenetic regulation and, in particular, DNA methylation have been linked to the underlying genetic sequence. DNA methylation quantitative trait loci (meQTL) have been identified through significant associations between the genetic and epigenetic codes in physiological and pathological contexts. We propose that interrogating the interplay between polymorphic alleles and DNA methylation is a powerful method for improving our interpretation of risk alleles identified in genome-wide association studies that otherwise lack mechanistic explanation. We integrated patient cancer risk genotype data and genome-scale DNA methylation profiles of 3,649 primary human tumors, representing 13 solid cancer types. We provide a comprehensive meQTL catalog containing DNA methylation associations for 21% of interrogated cancer risk polymorphisms. Differentially methylated loci harbor previously reported and as-yet-unidentified cancer genes. We suggest that such regulation at the DNA level can provide a considerable amount of new information about the biology of cancer-risk alleles., (Copyright © 2014 The Authors. Published by Elsevier Inc. All rights reserved.)

In cancer, the overall patterns of epigenetic marks are severely distorted from the corresponding normal cell type. It is now well established that these changes can contribute to cancer development through inactivation of tumor suppressor genes and, conversely, through activation of oncogenes. Recent technological advances have enabled epigenome-wide analyses of cancers that are yielding unexpected findings. The study of cancer epigenetics holds great promise for expanding the range of therapeutic opportunities for personalized medicine. Here, we focus on DNA methylation in breast cancer and the potential implications for clinical management of patients., (Copyright © 2013 American Society for Investigative Pathology. Published by Elsevier Inc. All rights reserved.)

It is not well known to what extent carrying a BRCA2 mutation affects the survival of women with breast cancer and prognostic factors among BRCA2-positive women warrant investigation. Using a record linkage approach we compared the long-term survival in carriers and noncarriers of an inherited BRCA2 founder mutation (999del5), and sought to identify prognostic factors among the BRCA2 mutation-positive subset, including markers of genetic instability (aneuploidy) and mitotic activity (S-phase fraction). We established the genetic status of 2,967 Icelandic breast cancer patients (215 mutation carriers and 2,752 noncarriers) diagnosed from 1955 to 2004, representing 72 % of all cases diagnosed in the country during this period. Tumour ploidy and S-phase fraction were assessed on tumour cells by DNA flow cytometry. Prognostic factors were assessed blindly with respect to mutation status. Univariate and multivariate hazard ratios (HR) were estimated for breast cancer-specific survival by BRCA2 status, using Cox regression. After a median follow-up of 9.5 years, BRCA2 mutation carriers had a higher risk of death from breast cancer than noncarriers (HR 1.64, 95 % CI 1.24-2.16, p < 0.001). The risk increase was restricted to women with diploid tumours (HR 3.03, 95 % CI 1.91-4.79, p < 0.001). Among breast cancer patients with aneuploid tumours, survival of carriers was similar to that of noncarriers (HR 0.76, 95 % CI 0.41-1.41, p = 0.38). Increased tumour size and a positive nodal status predicted worse prognosis in all patients, whereas the highly correlated prognostic factors diploidy, low proliferative activity and a positive estrogen receptor status had reverse effects in mutation carriers and noncarriers. Breast cancer patients who carry the Icelandic founder BRCA2 mutation have inferior long-term survival than noncarriers, but the adverse prognosis is restricted to mutation carriers with diploid, slowly proliferating tumours.

Germline mutations in the BRCA1 or BRCA2 genes are associated with an increased risk of breast and ovarian cancer development. Both genes are involved in DNA repair, and tumors harboring genetic defects in them are thought to be more sensitive to DNA-damaging agents used in chemotherapy. However, as only a minority of breast and ovarian cancer patients carry BRCA1 or BRCA2 mutations, few patients are likely to benefit from these pharmacogenetic biomarkers. Herein, we show that, in cancer cell lines and xenografted tumors, BRCA1 CpG island promoter hypermethylation-associated silencing also predicts enhanced sensitivity to platinum-derived drugs to the same extent as BRCA1 mutations. Most importantly, BRCA1 hypermethylation proves to be a predictor of longer time to relapse and improved overall survival in ovarian cancer patients undergoing chemotherapy with cisplatin.

BRCA1 is a tumor suppressor gene known to be implicated in the development of a subset of breast and ovarian cancers. The tumor suppressor properties of BRCA1 are generally thought to be linked to the gene's critical roles in the network of DNA damage response. In a recent report, BRCA1-mediated epigenetic repression at the promoter region of miRNA-155 was identified as a novel mechanism by which BRCA1 carries out its tumor suppressor functions.

Tetraploidy and aneuploidy can be caused by cell division errors and are frequently observed in many human carcinomas. We have recently reported delayed cytokinesis in primary human fibroblasts from BRCA2 mutation carriers, implying a function for the BRCA2 tumour suppressor in completion of cell division. Here, we address ploidy aberrations in breast tumours derived from BRCA2 germline mutation carriers. Ploidy aberrations were evaluated from flow cytometry histograms on selected breast tumour samples (n=236), previously screened for local BRCA mutations. The ploidy between BRCA2-mutated (n=71) and matched sporadic (n=165) cancers was compared. Differences in ploidy distribution were examined with respect to molecular tumour subtypes, previously defined by immunohistochemistry on tissue microarray sections. Tetraploidy was significantly 3 times more common in BRCA2 breast cancers than sporadic. However, no differences were found in the overall ploidy distribution between BRCA2-mutation carriers and non-carriers. In BRCA2 cancers, tetraploidy was associated with luminal characteristics. The increased frequency of tetraploidy in BRCA2 associated cancers may be linked to cell division errors, particularly cytokinesis. Additionally, tetraploidy emerges predominantly in BRCA2 breast cancers displaying luminal rather than triple-negative phenotypes., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Background: Inherited mutations in the BRCA2 gene greatly increase the risk of developing breast cancer. Consistent with an important role for BRCA2 in error-free DNA repair, complex genomic changes are frequently observed in tumors derived from BRCA2 mutation carriers. Here, we explore the impact of DNA copy-number changes in BRCA2 tumors with respect to phenotype and clinical staging of the disease., Methods: Breast tumors (n = 33) derived from BRCA2 999del5 mutation carriers were examined in terms of copy-number changes with high-resolution aCGH (array comparative genomic hybridization) containing 385 thousand probes (about one for each 7 kbp) and expression of phenotypic markers on TMAs (tissue microarrays). The data were examined with respect to clinical parameters including TNM staging, histologic grade, S phase, and ploidy., Results: Tumors from BRCA2 carriers of luminal and basal/triple-negative phenotypes (TNPs) differ with respect to patterns of DNA copy-number changes. The basal/TNP subtype was characterized by lack of pRb (RB1) coupled with high/intense expression of p16 (CDKN2A) gene products. We found increased proportions of Ki-67-positive cells to be significantly associated with loss of the wild-type (wt) BRCA2 allele in luminal types, whereas BRCA2wt loss was less frequent in BRCA2 tumors displaying basal/TNP phenotypes. Furthermore, we show that deletions at 13q13.1, involving the BRCA2wt allele, represents a part of a larger network of co-occurring genetic changes, including deletions at 6q22.32-q22.33, 11q14.2-q24.1, and gains at 17q24.1. Importantly, copy-number changes at these BRCA2-linked networking regions coincide with those associated with advanced progression, involving the capacity to metastasize to the nodes or more-distant sites at diagnosis., Conclusions: The results presented here demonstrate divergent paths of tumor evolution in BRCA2 carriers and that deletion of the wild-type BRCA2 allele, together with co-occurring changes at 6 q, 11 q, and 17 q, are important events in progression toward advanced disease.

Members of the Polycomb-group (PcG) family of proteins, including EZH2 (enhancer of zeste homolog 2), are involved in establishing epigenetic silencing of developmental genes in adult and embryonic stem cells, and their deregulation has been implicated in cancer. In a recent report, EZH2-mediated epigenetic repression of DNA damage repair in breast tumor initiating cells (BTICs) was identified as a mechanism that could promote expansion of BTICs, and may contribute to cancer progression.

Triple-negative breast cancer (TNBC) occurs in approximately 15% of all breast cancer patients, and the incidence of TNBC is greatly increased in BRCA1 mutation carriers. This study aimed to assess the impact of BRCA1 promoter methylation with respect to breast cancer subtypes in sporadic disease. Tissue microarrays (TMAs) were constructed representing tumors from 303 patients previously screened for BRCA1 germline mutations, of which a subset of 111 sporadic tumors had previously been analyzed with respect to BRCA1 methylation. Additionally, a set of eight tumors from BRCA1 mutation carriers were included on the TMAs. Expression analysis was performed on TMAs by immunohistochemistry (IHC) for BRCA1, pRb, p16, p53, PTEN, ER, PR, HER2, CK5/6, EGFR, MUC1 and Ki-67. Data on BRCA1 aberrations and IHC expression was examined with respect to breast cancer-specific survival. The results demonstrate that CpG island hypermethylation of BRCA1 significantly associates with the basal/triple-negative subtype. Low expression of pRb, and high/intense p16, were associated with BRCA1 promoter hypermethylation, and the same effects were seen in BRCA1 mutated tumors. The expression patterns of BRCA1, pRb, p16 and PTEN were highly correlated, and define a subgroup of TNBCs characterized by BRCA1 aberrations, high Ki-67 (≥ 40%) and favorable disease outcome. In conclusion, our findings demonstrate that epigenetic inactivation of the BRCA1 gene associates with RB/p16 dysfunction in promoting TNBCs. The clinical implications relate to the potential use of targeted treatment based on PARP inhibitors in sporadic TNBCs, wherein CpG island hypermethylation of BRCA1 represents a potential marker of therapeutic response.