GWAS meta-analysis reveals key risk loci in essential tremor pathogenesis

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.

A comparison of methods for detecting DNA methylation from long-read sequencing of human genomes

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.

Variants at the Interleukin 1 Gene Locus and Pericarditis

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.

Rare variants with large effects provide functional insights into the pathology of migraine subtypes, with and without aura

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.

Pregnancy-Associated Bleeding and Genetics: Five Sequence Variants in the Myometrium and Progesterone Signaling Pathway are associated with postpartum hemorrhage

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 severity of these complications in maternal-fetal health. Here, we investigated the genetic variation underlying aspects of pregnancy-associated bleeding and identified five loci associated with PPH through a meta-analysis of 21,512 cases and 259,500 controls. 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 (PGR). Furthermore, there were strong genetic correlations with birth weight, gestational duration, and uterine fibroids. Early bleeding during pregnancy (28,898 cases and 302,894 controls) yielded no genome-wide association signals, but showed strong genetic correlation with a variety of human traits, indicative of polygenic and pleiotropic effects. Our results suggest that postpartum bleeding is related to myometrium dysregulation, whereas early bleeding is a complex trait related to underlying health and possibly socioeconomic status.

Sequence variants affecting the genome-wide rate of germline microsatellite mutations

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.

Deficit of homozygosity among 1.52 million individuals and genetic causes of recessive lethality

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.

Genetic variants associated with syncope implicate neural and autonomic processes

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.

Meta-analysis of erosive hand osteoarthritis identifies four common variants that associate with relatively large effect

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.

The sequences of 150,119 genomes in the UK Biobank

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

Genetic predisposition to hypertension is associated with preeclampsia in European and Central Asian women

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.

Genome-wide association identifies seven loci for pelvic organ prolapse in Iceland and the UK Biobank

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.

Eighty-eight variants highlight the role of T cell regulation and airway remodeling in asthma pathogenesis

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.

Asthma is a common allergic airway disease with significant inter-individual heterogeneity. Here, Olafsdottir et al. report a genome-wide meta-analysis of two large population-based cohorts to identify sequence variants that associate with asthma risk and perform follow-up functional analyses on a protective loss-of-function variant in TNFRSF8.

BRCA1 Promoter Methylation Status in 1031 Primary Breast Cancers Predicts Favorable Outcomes Following Chemotherapy

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.

GWAS of bone size yields twelve loci that also affect height, BMD, osteoarthritis or fractures

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.

Characterizing mutagenic effects of recombination through a sequence-level genetic map

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.

A loss-of-function variant in ALOX15 protects against nasal polyps and chronic rhinosinusitis

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.

Genome-wide association meta-analysis yields 20 loci associated with gallstone disease

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.

Variants associating with uterine leiomyoma highlight genetic background shared by various cancers and hormone-related traits

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.

Uterine leiomyomas are common benign tumors. Here, a meta-analysis of two European leiomyoma GWAS uncovers 21 leiomyoma risk variants at 16 loci, providing evidence of genetic overlap between leiomyoma and various benign and malignant tumors and highlighting the role of estrogen in tumor growth.

Abstract 3307: BRCA1 promoter methylation analyzed in 1032 primary breast cancers predicts favorable response to chemotherapy independently of p53 mutations

Individuals with germline mutations in the breast cancer susceptibility gene, BRCA1, are at increased risk of developing breast- and ovarian cancer. BRCA1 is involved in DNA double strand break repair (DSB) by homologous recombination (HR), an error-free DNA repair pathway that uses an intact sister chromatid for the repair. BRCA1 mutated tumors are therefore HR deficient and sensitive to drugs that stall DNA replication forks, resulting in DSBs. Tumor cells lacking the protein are characterized by genomic instability and a pattern of specific genome-wide mutation signature as a consequence of dependence on alternative error prone repair pathways. BRCA1 can also be inactivated epigenetically by CpG promoter hypermethylation. These BRCA1 epigenetically silenced tumors are more likely to be of the basal-like subtype and are strongly associated with this pattern of mutational signature characteristic for individuals with germline mutations in BRCA1. This suggests that the tumorigenesis in the sporadic cancer cases is similar to tumors where germline mutations are present, highlighting the importance of epigenetic silencing of BRCA1. The aim of the study was to examine the clinical significance of epigenetic silencing of BRCA1 in primary breast cancer, focusing on the response to chemotherapy. DNA methylation analysis of the BRCA1 promoter was performed by pyrosequencing tumor samples from 1032 patients with primary breast cancer of which 966 were sporadic, 61 were BRCA2 germline mutation carriers and 5 were BRCA1 mutation carriers. The p53 mutational status was also studied in the same primary tumors by DNA sequencing. We analysed the methylation status of BRCA1 in tumors derived from 1032 patients and found BRCA1 to be promoter hypermethylated in 29 sporadic tumors (29 out of 966, 3.0 %) whereas none of the tumors derived from BRCA mutation carriers were BRCA1 promoter hypermethylated. Patients with BRCA1 promoter hypermethylation receiving chemotherapeutic drug treatment show significantly better breast cancer specific survival (P=0.0192). The trend is opposite in patients not receiving chemotherapeutic treatment, were patients with BRCA1 hypermethylated tumors show worse survival (P=0.074). The impact of BRCA1 methylation with respect to survival in chemotherapy treated patients was found to be independent of p53 mutational status, ER status, tumorsize, lymph node status, year- and age at diagnosis (HR = 0.108; 95% CI = 0.015-0.794; P = 0.0288). BRCA1 promoter hypermethylation predicts improved disease outcome in anti-cancer drug treated patients. Therefore, BRCA1 promoter methylation status in primary breast cancer is an important factor to be considered when deciding treatment options.

Citation Format: Olafur A. Stefansson, Holmfridur Hilmarsdottir, Kristrun Olafsdottir, Laufey Tryggvadottir, Asgerdur Sverrisdottir, Jon G. Jonasson, Jorunn E. Eyfjord, Stefan Sigurdsson. BRCA1 promoter methylation analyzed in 1032 primary breast cancers predicts favorable response to chemotherapy independently of p53 mutations [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3307.

Genome-wide analysis yields new loci associating with aortic valve stenosis

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⁻²²) and on chromosome 2q22 in TEX41 (rs1830321; OR = 1.15, P = 1.8 × 10⁻¹³). Rs7543130 also associates with bicuspid aortic valve (BAV) (OR = 1.28, P = 6.6 × 10⁻¹⁰) and aortic root diameter (P = 1.30 × 10⁻⁸), and rs1830321 associates with BAV (OR = 1.12, P = 5.3 × 10⁻³) and coronary artery disease (OR = 1.05, P = 9.3 × 10⁻⁵). 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.

Aortic valve stenosis (AS) is the most common valvular heart disease. Here the authors identify two new AS loci that also associate with bicuspid aortic valve, aortic root diameter and/or coronary artery disease implicating both developmental abnormalities and atherosclerosis-like processes in AS.

Obesity and menopause modify the epigenomic profile of breast cancer

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.

Telomere Length Is Predictive of Breast Cancer Risk in BRCA2 Mutation Carriers

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.

Landscape of somatic mutations in 560 breast cancer whole-genome sequences

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.

Aurora A is a prognostic marker for breast cancer arising in BRCA2 mutation carriers

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.

A DNA methylation-based definition of biologically distinct breast cancer subtypes

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.

We describe distinct signatures associated with luminal‐B and basal‐like subtypes.

The signatures identified show differences in genes, but also in the CpG context.

A selected set of proxy markers for each signature revealed their clinical relevance.

A comprehensive DNA methylation profile of epithelial-to-mesenchymal transition

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.

CARM1 and BAF155: an example of how chromatin remodeling factors can be relocalized and contribute to cancer

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.

Tumour diploidy and survival in breast cancer patients with BRCA2 mutations

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.

Abstract 2018: Genome-wide analysis of DNA methylation and copy number changes in breast cancers

Epigenetic modifications are known to be critical in regulating normal developmental processes and cellular differentiation. DNA methylation of CpG dinucleotides is a well-characterized epigenetic marker frequently associated with transcriptional silencing. CpG dinucleotides are non-randomly distributed throughout the genome with regions enriched for CpGs being known as CpG islands. It is now well established that epigenetic inactivation of tumor suppressor genes can occur by CpG island promoter methylation. For example, BRCA1 promoter methylation occurs in approximately 10 - 15% of all breast cancers and is associated with gene silencing. Recent technological advances have enabled DNA methylation to be profiled genome-wide. Using the Infinium 450 thousand probe arrays, we have studied DNA methylation and copy number changes to discover new markers of potential clinical relevance in breast cancer. Expression analyses by immunohistochemisty on tissue arrays for subtype-specific markers (ER, PR, HER2, Ki-67, EGFR, CK5/6, CK8, CK18, MUC1 and Claudin-3) carried out previously were used to establish DNA methylation markers of breast cancer subtypes. This integrated approach is presented in clinical context with available data on histological grade, disease stage and breast cancer-specific survival outcomes.

Citation Format: Olafur A. Stefansson, Sebastian Moran, Antonio Gomez, Sergi Sayols Puig, Jorunn Eyfjord, Manel Esteller. Genome-wide analysis of DNA methylation and copy number changes in breast cancers. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 2018. doi:10.1158/1538-7445.AM2013-2018

BRCA1 epigenetic inactivation predicts sensitivity to platinum-based chemotherapy in breast and ovarian cancer

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.

DNA methylation profiling in breast cancer discordant identical twins identifies DOK7 as novel epigenetic biomarker

Using whole blood from 15 twin pairs discordant for breast cancer and high-resolution (450K) DNA methylation analysis, we identified 403 differentially methylated CpG sites including known and novel potential breast cancer genes. Confirming the results in an independent validation cohort of 21 twin pairs determined the docking protein DOK7 as a candidate for blood-based cancer diagnosis. DNA hypermethylation of the promoter region was also seen in primary breast cancer tissues and cancer cell lines. Hypermethylation of DOK7 occurs years before tumor diagnosis, suggesting a role as a powerful epigenetic blood-based biomarker as well as providing insights into breast cancer pathogenesis.

Genomic and phenotypic analysis of BRCA2 mutated breast cancers reveals co-occurring changes linked to progression

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.

Epigenetic silencing and deletion of the BRCA1 gene in sporadic breast cancer

INTRODUCTION

BRCA1 or BRCA2 germline mutations increase the risk of developing breast cancer. Tumour cells from germline mutation carriers have frequently lost the wild-type allele. This is predicted to result in genomic instability where cell survival depends upon dysfunctional checkpoint mechanisms. Tumorigenic potential could then be acquired through further genomic alterations. Surprisingly, somatic BRCA mutations are not found in sporadic breast tumours. BRCA1 methylation has been shown to occur in sporadic breast tumours and to be associated with reduced gene expression. We examined the frequency of BRCA1 methylation in 143 primary sporadic breast tumours along with BRCA1 copy number alterations and tumour phenotype.

METHODS

Primary sporadic breast tumours were analysed for BRCA1alpha promoter methylation by methylation specific PCR and for allelic imbalance (AI) at BRCA1 and BRCA2 loci by microsatellite analysis and TP53 (also known as p53) mutations by constant denaturing gel electrophoresis. The BRCA1 methylated tumours were analysed for BRCA1 copy alterations by fluorescence in situ hybridisation and BRCA1 expression by immunostaining.

RESULTS

BRCA1 methylation was found in 13/143 (9.1%) sporadic breast tumours. The BRCA1 methylated tumours were significantly associated with estrogen receptor (ER) negativity (P = 0.0475) and displayed a trend for BRCA1 AI (P = 0.0731) as well as young-age at diagnosis (< or = 55; P = 0.0898). BRCA1 methylation was not associated with BRCA2 AI (P = 0.5420), although a significant association was found between BRCA1 AI and BRCA2 AI (P < 0.0001).Absent/markedly reduced BRCA1 expression was observed in 9/13 BRCA1 methylated tumours, most of which had BRCA1 deletion. An elevated TP53 mutation frequency was found among BRCA1 methylated tumours (38.5%) compared with non-methylated tumours (17.2%). The BRCA1 methylated tumours were mainly of tumour grade 3 (7/13) and infiltrating ductal type (12/13). Only one methylated tumour was of grade 1.

CONCLUSION

BRCA1 methylation is frequent in primary sporadic breast tumours. We found an indication for BRCA1 methylation to be associated with AI at the BRCA1 locus. Almost all BRCA1 methylated tumours with absent/markedly reduced BRCA1 expression (8/9) displayed BRCA1 deletion. Thus, epigenetic silencing and deletion of the BRCA1 gene might serve as Knudson's two 'hits' in sporadic breast tumorigenesis. We observed phenotypic similarities between BRCA1 methylated and familial BRCA1 tumours, based on BRCA1 deletion, TP53 mutations, ER status, young age at diagnosis and tumour grade.