Truncating NFKB1 variants cause combined NLRP3 inflammasome activation and type I interferon signaling and predispose to necrotizing fasciitis

In monogenic autoinflammatory diseases, mutations in genes regulating innate immune responses often lead to uncontrolled activation of inflammasome pathways or the type I interferon (IFN-I) response. We describe a mechanism of autoinflammation potentially predisposing patients to life-threatening necrotizing soft tissue inflammation. Six unrelated families are identified in which affected members present with necrotizing fasciitis or severe soft tissue inflammations. Exome sequencing reveals truncating monoallelic loss-of-function variants of nuclear factor κ light-chain enhancer of activated B cells (NFKB1) in affected patients. In patients' macrophages and in NFKB1-variant-bearing THP-1 cells, activation increases both interleukin (IL)-1β secretion and IFN-I signaling. Truncation of NF-κB1 impairs autophagy, accompanied by the accumulation of reactive oxygen species and reduced degradation of inflammasome receptor nucleotide-binding oligomerization domain, leucine-rich repeat-containing protein 3 (NLRP3), and Toll/IL-1 receptor domain-containing adaptor protein inducing IFN-β (TRIF), thus leading to combined excessive inflammasome and IFN-I activity. Many of the patients respond to anti-inflammatory treatment, and targeting IL-1β and/or IFN-I signaling could represent a therapeutic approach for these patients.

Mutational signatures and their association with survival and gene expression in urological carcinomas

Different sources of mutagenesis cause consistently identifiable patterns of mutations and mutational signatures that mirror the various carcinogenetic processes. We used publicly available data from the Cancer Genome Atlas to evaluate the associations between the activity of the mutational signatures and various survival endpoints in six types of urological cancers after adjusting for established prognostic factors. The predictive power of the signatures was evaluated with dynamic area under curve models. In addition, links between mutational signature activities and differences in gene expression patterns were analysed. APOBEC-related signature SBS2 was associated with improved overall survival (OS) and disease-specific survival (DSS) in bladder carcinomas in the multivariate analysis, while clock-like signature SBS1 predicted shortened DSS and progression-free interval (PFI) in clear cell renal cell carcinomas (ccRCC). In papillary renal cell carcinomas (pRCC), SBS45 was a predictor of improved outcomes, and APOBEC-related SBS13 was a predictor of worse outcomes. Gene expression analyses revealed various enriched pathways between the low- and high-signature groups. Interestingly, in both the ccRCC and pRCC cohorts, the genes of several members of the melanoma antigen (MAGE) family were highly upregulated in the signatures, which predicted poor outcomes, and downregulated in signatures, which were associated with improved survival. To summarize, SBS signatures provide substantial prognostic value compared with just the traditional prognostic factors in certain cancer types. APOBEC-related SBS2 and SBS13 seem to provide robust prognostic information for particular urological cancers, maybe driven by the expression of specific groups of genes, including the MAGE gene family.

DeepIFC: Virtual fluorescent labeling of blood cells in imaging flow cytometry data with deep learning

Imaging flow cytometry (IFC) combines flow cytometry with microscopy, allowing rapid characterization of cellular and molecular properties via high-throughput single-cell fluorescent imaging. However, fluorescent labeling is costly and time-consuming. We present a computational method called DeepIFC based on the Inception U-Net neural network architecture, able to generate fluorescent marker images and learn morphological features from IFC brightfield and darkfield images. Furthermore, the DeepIFC workflow identifies cell types from the generated fluorescent images and visualizes the single-cell features generated in a 2D space. We demonstrate that rarer cell types are predicted well when a balanced data set is used to train the model, and the model is able to recognize red blood cells not seen during model training as a distinct entity. In summary, DeepIFC allows accurate cell reconstruction, typing and recognition of unseen cell types from brightfield and darkfield images via virtual fluorescent labeling.

Rare variants in complement system genes associate with endothelial damage after pediatric allogeneic hematopoietic stem cell transplantation

Introduction

Complement system has a postulated role in endothelial problems after hematopoietic stem cell transplantation (HSCT). In this retrospective, singlecenter study we studied genetic complement system variants in patients with documented endotheliopathy. In our previous study among pediatric patients with an allogeneic HSCT (2001-2013) at the Helsinki University Children´s Hospital, Finland, we identified a total of 19/122 (15.6%) patients with vascular complications, fulfilling the criteria of capillary leak syndrome (CLS), venoocclusive disease/sinusoidal obstruction syndrome (VOD/SOS) or thrombotic microangiopathy (TMA).

Methods

We performed whole exome sequencing (WES) on 109 patients having an adequate pre-transplantation DNA for the analysis to define possible variations and mutations potentially predisposing to functional abnormalities of the complement system. In our data analysis, we focused on 41 genes coding for complement components.

Results

50 patients (45.9%) had one or several, nonsynonymous, rare germline variants in complement genes. 21/66 (31.8%) of the variants were in the terminal pathway. Patients with endotheliopathy had variants in different complement genes: in the terminal pathway (C6 and C9), lectin pathway (MASP1) and receptor ITGAM (CD11b, part of CR3). Four had the same rare missense variant (rs183125896; Thr279Ala) in the C9 gene. Two of these patients were diagnosed with endotheliopathy and one with capillary leak syndrome-like problems. The C9 variant Thr279Ala has no previously known disease associations and is classified by the ACMG guidelines as a variant of uncertain significance (VUS). We conducted a gene burden test with gnomAD Finnish (fin) as the reference population. Complement gene variants seen in our patient population were investigated and Total Frequency Testing (TFT) was used for execution of burden tests. The gene variants seen in our patients with endotheliopathy were all significantly (FDR < 0.05) enriched compared to gnomAD. Overall, 14/25 genes coding for components of the complement system had an increased burden of missense variants among the patients when compared to the gnomAD Finnish population (N=10 816).

Discussion

Injury to the vascular endothelium is relatively common after HSCT with different phenotypic appearances suggesting yet unidentified underlying mechanisms. Variants in complement components may be related to endotheliopathy and poor prognosis in these patients.

The clinical picture of ERCC6L2 disease: from bone marrow failure to acute leukemia

Biallelic germ line excision repair cross-complementing 6 like 2 (ERCC6L2) variants strongly predispose to bone marrow failure (BMF) and myeloid malignancies, characterized by somatic TP53-mutated clones and erythroid predominance. We present a series of 52 subjects (35 families) with ERCC6L2 biallelic germ line variants collected retrospectively from 11 centers globally, with a follow-up of 1165 person-years. At initial investigations, 32 individuals were diagnosed with BMF and 15 with a hematological malignancy (HM). The subjects presented with 19 different variants of ERCC6L2, and we identified a founder mutation, c.1424delT, in Finnish patients. The median age of the subjects at baseline was 18 years (range, 2-65 years). Changes in the complete blood count were mild despite severe bone marrow (BM) hypoplasia and somatic TP53 mutations, with no significant difference between subjects with or without HMs. Signs of progressive disease included increasing TP53 variant allele frequency, dysplasia in megakaryocytes and/or erythroid lineage, and erythroid predominance in the BM morphology. The median age at the onset of HM was 37.0 years (95% CI, 31.5-42.5; range, 12-65 years). The overall survival (OS) at 3 years was 95% (95% CI, 85-100) and 19% (95% CI, 0-39) for patients with BMF and HM, respectively. Patients with myelodysplastic syndrome or acute myeloid leukemia with mutated TP53 undergoing hematopoietic stem cell transplantation had a poor outcome with a 3-year OS of 28% (95% CI, 0-61). Our results demonstrated the importance of early recognition and active surveillance in patients with biallelic germ line ERCC6L2 variants.

Germline ERCC excision repair 6 like 2 (ERCC6L2) mutations lead to impaired erythropoiesis and reshaping of the bone marrow microenvironment

Despite the inclusion of inherited myeloid malignancies as a separate entity in the World Health Organization Classification, many established predisposing loci continue to lack functional characterization. While germline mutations in the DNA repair factor ERCC excision repair 6 like 2 (ERCC6L2) give rise to bone marrow failure and acute myeloid leukaemia, their consequences on normal haematopoiesis remain unclear. To functionally characterise the dual impact of germline ERCC6L2 loss on human primary haematopoietic stem/progenitor cells (HSPCs) and mesenchymal stromal cells (MSCs), we challenged ERCC6L2-silenced and patient-derived cells ex vivo. Here, we show for the first time that ERCC6L2-deficiency in HSPCs significantly impedes their clonogenic potential and leads to delayed erythroid differentiation. This observation was confirmed by CIBERSORTx RNA-sequencing deconvolution performed on ERCC6L2-silenced erythroid-committed cells, which demonstrated higher proportions of polychromatic erythroblasts and reduced orthochromatic erythroblasts versus controls. In parallel, we demonstrate that the consequences of ERCC6L2-deficiency are not limited to HSPCs, as we observe a striking phenotype in patient-derived and ERCC6L2-silenced MSCs, which exhibit enhanced osteogenesis and suppressed adipogenesis. Altogether, our study introduces a valuable surrogate model to study the impact of inherited myeloid mutations and highlights the importance of accounting for the influence of germline mutations in HSPCs and their microenvironment.

A germline exome analysis reveals harmful POT1 variants in multiple myeloma patients and families

Observations of inherited susceptibility to multiple myeloma have led to active research in defining predisposing genes to the disease. Here, we analysed 128 plasma cell dyscrasia patients' germline whole-exome sequencing data. Rare dominantly inherited pathogenic or likely pathogenic (P/LP) variant was found in 9.4% of the patients. Among the P/LP variants, CHEK2 (p. Thr410MetfsTer15) was the most prevalent (n = 5, 3.9%). Interestingly, P/LP variants in POT1 were identified in three patients (2.3%). Our findings broaden the spectrum of POT1-related cancers and demonstrate the importance of the germline genetic analysis in hematological malignancies.

Mutational Signatures Associate With Survival in Gastrointestinal Carcinomas

BACKGROUND/AIM

Mutational signatures reflect common patterns based on the counts of mutations and their sequence context. The prognostic value of these signatures, mirroring various carcinogenetic processes of cancers, are unexplored in gastrointestinal cancers. Our aim was to evaluate possible prognostic relevance of mutational signatures in gastrointestinal carcinomas after adjusting with the traditional prognostic factors.

MATERIALS AND METHODS

We used publicly available data from The Cancer Genome Atlas and Pan-Cancer Analysis of Whole Genomes to evaluate the associations between survival endpoints and activity of mutational signatures in seven types of gastrointestinal cancers.

RESULTS

Most strikingly, the high activity of age-related single-base substitution 5 (SBS5) and SBS40 signatures were in rectal adenocarcinomas associated with both improved overall survival (OS) [for SBS5 hazard ratio (HR) 0.130; 95% CI=0.03-0.56, for SBS40 HR=0.072; 95% CI=0.012-0.44, respectively] and similarly also to rectal cancer-specific survival. In patients with left-sided (but not right-sided) colon adenocarcinoma, the high activity of SBS2 signatures, formed due to APOBEC activity, predicted shortened OS. In pancreatic cancer, the high activity of SBS10b, caused by polymerase epsilon exonuclease proofreading defects, was associated both with longer OS (HR=0.44; 95% CI=0.205-0.96) and pancreatic cancer-specific survival (HR=0.32; 95% CI=0.112-0.91).

CONCLUSION

Several mutational signatures seem to have clinically meaningful, cancer-specific associations with prognosis among gastrointestinal cancers.

Enrichment of cancer-predisposing germline variants in adult and pediatric patients with acute lymphoblastic leukemia

Despite recent progress in acute lymphoblastic leukemia (ALL) therapies, a significant subset of adult and pediatric ALL patients has a dismal prognosis. Better understanding of leukemogenesis and recognition of germline genetic changes may provide new tools for treating patients. Given that hematopoietic stem cell transplantation, often from a family member, is a major form of treatment in ALL, acknowledging the possibility of hereditary predisposition is of special importance. Reports of comprehensive germline analyses performed in adult ALL patients are scarce. Aiming at fulfilling this gap of knowledge, we investigated variants in 93 genes predisposing to hematologic malignancies and 70 other cancer-predisposing genes from exome data obtained from 61 adult and 87 pediatric ALL patients. Our results show that pathogenic (P) or likely pathogenic (LP) germline variants in genes associated with predisposition to ALL or other cancers are prevalent in ALL patients: 8% of adults and 11% of children. Comparison of P/LP germline variants in patients to population-matched controls (gnomAD Finns) revealed a 2.6-fold enrichment in ALL cases (CI 95% 1.5–4.2, p = 0.00071). Acknowledging inherited factors is crucial, especially when considering hematopoietic stem cell transplantation and planning post-therapy follow-up. Harmful germline variants may also predispose patients to excessive toxicity potentially compromising the outcome. We propose integrating germline genetics into precise ALL patient care and providing families genetic counseling.

Familial aggregation of early-onset haematological malignancies

Population-based studies on familial aggregation of haematological malignancies (HM) have rarely focused specifically on early-onset HMs. We estimated standardized incidence ratios (SIR) and cumulative risks of relatives with Hodgkin lymphoma (HL), non-Hodgkin lymphomas (NHL), acute lymphoblastic leukaemia/lymphoma (ALL/LBL) and acute myeloid leukaemia (AML) when index persons and relatives were diagnosed with early-onset HM. A total of 8791 patients aged ≤40 years and diagnosed with primary HM in Finland from 1970 to 2012 were identified from the Finnish Cancer Registry and their 75 774 family members were retrieved from the population registry. SIRs for concordant HMs were elevated among first-degree relatives in all of the most common HMs of children and adolescents and young adults (AYA). The risk was highest among siblings with HL (SIR 9·09, 95% confidence interval 5·55-14·04) and AML (8·29, 1·00-29·96). HL also had the highest cumulative risk for siblings at ≤40 years of age (0·92% vs. 0·11% in the population). In conclusion, significantly elevated SIRs indicate a role of shared aetiological factors in some families, which should be noted in the clinical setting when caring for patients with early-onset HMs.

Retrotransposon insertions can initiate colorectal cancer and are associated with poor survival

Genomic instability pathways in colorectal cancer (CRC) have been extensively studied, but the role of retrotransposition in colorectal carcinogenesis remains poorly understood. Although retrotransposons are usually repressed, they become active in several human cancers, in particular those of the gastrointestinal tract. Here we characterize retrotransposon insertions in 202 colorectal tumor whole genomes and investigate their associations with molecular and clinical characteristics. We find highly variable retrotransposon activity among tumors and identify recurrent insertions in 15 known cancer genes. In approximately 1% of the cases we identify insertions in APC, likely to be tumor-initiating events. Insertions are positively associated with the CpG island methylator phenotype and the genomic fraction of allelic imbalance. Clinically, high number of insertions is independently associated with poor disease-specific survival.

Retrotransposons are usually dormant in healthy tissue, but become activated during malignancy. Here, in colorectal cancer, Cajuso et al. show that retrotransposon activity associates with clinical features of the disease.

Comprehensive evaluation of coding region point mutations in microsatellite-unstable colorectal cancer

Microsatellite instability (MSI) leads to accumulation of an excessive number of mutations in the genome, mostly small insertions and deletions. MSI colorectal cancers (CRCs), however, also contain more point mutations than microsatellite‐stable (MSS) tumors, yet they have not been as comprehensively studied. To identify candidate driver genes affected by point mutations in MSI CRC, we ranked genes based on mutation significance while correcting for replication timing and gene expression utilizing an algorithm, MutSigCV. Somatic point mutation data from the exome kit‐targeted area from 24 exome‐sequenced sporadic MSI CRCs and respective normals, and 12 whole‐genome‐sequenced sporadic MSI CRCs and respective normals were utilized. The top 73 genes were validated in 93 additional MSI CRCs. The MutSigCV ranking identified several well‐established MSI CRC driver genes and provided additional evidence for previously proposed CRC candidate genes as well as shortlisted genes that have to our knowledge not been linked to CRC before. Two genes, SMARCB1 and STK38L, were also functionally scrutinized, providing evidence of a tumorigenic role, for SMARCB1 mutations in particular.

Impact of constitutional TET2 haploinsufficiency on molecular and clinical phenotype in humans

Clonal hematopoiesis driven by somatic heterozygous TET2 loss is linked to malignant degeneration via consequent aberrant DNA methylation, and possibly to cardiovascular disease via increased cytokine and chemokine expression as reported in mice. Here, we discover a germline TET2 mutation in a lymphoma family. We observe neither unusual predisposition to atherosclerosis nor abnormal pro-inflammatory cytokine or chemokine expression. The latter finding is confirmed in cells from three additional unrelated TET2 germline mutation carriers. The TET2 defect elevates blood DNA methylation levels, especially at active enhancers and cell-type specific regulatory regions with binding sequences of master transcription factors involved in hematopoiesis. The regions display reduced methylation relative to all open chromatin regions in four DNMT3A germline mutation carriers, potentially due to TET2-mediated oxidation. Our findings provide insight into the interplay between epigenetic modulators and transcription factor activity in hematological neoplasia, but do not confirm the putative role of TET2 in atherosclerosis.

Abstract 2535: Single-cell transcriptome analysis on lymphocytes of GATA2 deficiency patients

GATA2 is an essential transcription factor in the proliferation and maintenance of hematopoietic stem and progenitor cells. Germline GATA2 mutations cause hematological malignancies, especially myelodysplastic syndrome and acute myeloid leukemia. Pre-leukemic phenotype of GATA2 deficient patients varies markedly. However, peripheral blood (PB) lymphopenia is characteristic for almost all the patients and encounters most prominently B and natural killer (NK) cells. We aimed to gain understanding on the effects of GATA2 deficiency in hematopoiesis by comparing expression patterns in lymphocytes in GATA2 wild type and GATA2 deficient individuals.

We conducted single-cell transcriptomic analysis using 10x Genomics Chromium Single Cell 3′ Solution on PB samples from three healthy controls and three GATA2 deficient patients (p.Q328X, p.T354M, and one haploinsufficiency due to unidentified reason). Erythrocytes were lysed prior to washing, filtering and resuspending white blood cells in 0.04% BSA/PBS (1x106 cells/mL). The gel beads in emulsion were generated using Chromium Single Cell 3' v2 Reagent Kit. The sample libraries, prepared according to 10x Genomics instructions, were sequenced in Illumina HiSeq 2500. Cell Ranger v1.3 analysis pipelines (10x Genomics) were used to demultiplex and convert Chromium single cell 3' RNA-sequencing barcode; read data to FASTQ files; and generate align reads and gene-cell matrices. The data were analyzed using Seurat, an R package designed for single cell analysis. Cell type clusters were identified through known markers from Seurat's graph-based clustering.

Significant gene expression differences in lymphocytes (GATA2 deficient patients vs. controls) were observed (Wilcoxon rank sum test, Bonferroni-corrected p &lt;0.05). In NK cells the distinction was found e.g. in KLRB1, CD74, CCL5, STAT1, and LAG3, all linked to immune responses, and BTG1 and DDIT4, regulating cell proliferation and differentiation. Expression changes in immune system-associated genes were also noted in T cells e.g. in STAT1, IL7R, CD48, and CXCR4, as well as in NPM1, TPT1, and TXNIP, genes identified to participate in diverse cellular processes. Also, several genes encoding ribosomal proteins were differentially expressed between patients' and controls' T cells. Comparison of B cells revealed changes in e.g. GZMH, PPBP, and LCN2, which, too, play roles in the immune system.

We observed significant dissimilarities in gene expression pattern in GATA2 deficient patients' vs. controls' lymphocytes. These were reflected in several genes involved in immune system and cellular processes. Connection between our findings and patients' clinical manifestations is intriguing and may have an impact in the development of future patient care. We will extend our material and include also bone marrow sample single-cell transcriptomes and healthy mutation carriers to the analysis.

Citation Format: Elina A. Hirvonen, Inari Natri, Pirkko Mattila, Jenni Mattila, Esa Pitkänen, Kimmo Porkka, Outi Kilpivaara, Ulla Wartiovaara-Kautto. Single-cell transcriptome analysis on lymphocytes of GATA2 deficiency patients [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 2535.

Abstract 4381: The mobile genome of colorectal cancer: Characterization of retrotransposon insertions in 202 colorectal cancer whole genomes

The aim of this study is to characterize the landscape of retrotransposon insertions in 202 colorectal cancer (CRC) patients to elucidate their role in colorectal tumorigenesis. Retrotransposons are genetic sequences that can copy themselves into an mRNA intermediate and insert elsewhere in the genome. Germline retrotranspositions can contribute to genetic variation however, high retrotransposon activity can also lead to genetic instability. Albeit retrotransposons are usually repressed in normal adult tissues, they become highly active in several cancer types, being most active in epithelial cancers such as CRC. Moreover, retrotransposon insertions have been difficult to detect with previous methodological approaches, which has resulted in very few genome-wide studies.

Whole genome sequencing (WGS) was performed in 202 tumors and their corresponding normal tissue on the Illumina HiSeq 2000 platform; with 100 base-pairs paired-end reads. Each sample was sequenced to a minimum of 40x median coverage. To detect retrotransposon insertions we utilized TraFiC (Transposome Finder in Cancer). To identify retrotransposon-mediated transductions that were not identifiable using TraFiC, we utilized DELLY and the European database of L1-HS retrotransposon insertions in humans (euL1db). To characterize somatic insertions and decrease the rate of false positives, all insertions present in any of the normal genomes were filtered. After careful visual inspection of the paired-end read data, we estimated 80% of true somatic calls.

We detected high retrotransposon activity in CRC with remarkable variation among patients. The mean number of insertions per tumor was 25, ranging from 0 to 197 insertions. In addition, we were able to identify few highly active retrotransposons that accounted for 70% of the retrotransposon-mediated transductions. We identified a total of 5065 somatic insertions, 2% of these insertions were located in exons, whereas 45% were in introns. We identified several genes with recurrent insertions, some of these loci being known fragile sites. However, we identified 13 known cancer genes with two or more insertions. Furthermore, we identified two patients with insertions in exon 16 of APC, suggesting that these insertions could be initiating tumorigenic events. To conclude, we found retrotransposon insertions to be a common phenomenon in CRC, contributing to genome instability and colorectal tumorigenesis.

Citation Format: Tatiana Cajuso, Päivi Sulo, Riku Katainen, Ulrika Hänninen, Tomas Tanskanen, Johanna Kondelin, Jiri Hamberg, Niko Välimäki, Kimmo Palin, Outi Kilpivaara, Esa Pitkänen, Lauri A. Aaltonen. The mobile genome of colorectal cancer: Characterization of retrotransposon insertions in 202 colorectal cancer whole genomes [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 4381. doi:10.1158/1538-7445.AM2017-4381

Whole-exome sequencing identifies novel candidate predisposition genes for familial polycythemia vera

Background

Polycythemia vera (PV), characterized by massive production of erythrocytes, is one of the myeloproliferative neoplasms. Most patients carry a somatic gain-of-function mutation in JAK2, c.1849G > T (p.Val617Phe), leading to constitutive activation of JAK-STAT signaling pathway. Familial clustering is also observed occasionally, but high-penetrance predisposition genes to PV have remained unidentified.

Results

We studied the predisposition to PV by exome sequencing (three cases) in a Finnish PV family with four patients. The 12 shared variants (maximum allowed minor allele frequency <0.001 in Finnish population in ExAC database) predicted damaging in silico and absent in an additional control set of over 500 Finns were further validated by Sanger sequencing in a fourth affected family member. Three novel predisposition candidate variants were identified: c.1254C > G (p.Phe418Leu) in ZXDC, c.1931C > G (p.Pro644Arg) in ATN1, and c.701G > A (p.Arg234Gln) in LRRC3. We also observed a rare, predicted benign germline variant c.2912C > G (p.Ala971Gly) in BCORL1 in all four patients. Somatic mutations in BCORL1 have been reported in myeloid malignancies. We further screened the variants in eight PV patients in six other Finnish families, but no other carriers were found.

Conclusions

Exome sequencing provides a powerful tool for the identification of novel variants, and understanding the familial predisposition of diseases. This is the first report on Finnish familial PV cases, and we identified three novel candidate variants that may predispose to the disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s40246-017-0102-x) contains supplementary material, which is available to authorized users.

CTCF/cohesin-binding sites are frequently mutated in cancer

Cohesin is present in almost all active enhancer regions, where it is associated with transcription factors. Cohesin frequently colocalizes with CTCF (CCCTC-binding factor), affecting genomic stability, expression and epigenetic homeostasis. Cohesin subunits are mutated in cancer, but CTCF/cohesin-binding sites (CBSs) in DNA have not been examined for mutations. Here we report frequent mutations at CBSs in cancers displaying a mutational signature where mutations in A•T base pairs predominate. Integration of whole-genome sequencing data from 213 colorectal cancer (CRC) samples and chromatin immunoprecipitation sequencing (ChIP-exo) data identified frequent point mutations at CBSs. In contrast, CRCs showing an ultramutator phenotype caused by defects in the exonuclease domain of DNA polymerase ɛ (POLE) displayed significantly fewer mutations at and adjacent to CBSs. Analysis of public data showed that multiple cancer types accumulate CBS mutations. CBSs are a major mutational hotspot in the noncoding cancer genome.

Frequent L1 retrotranspositions originating from TTC28 in colorectal cancer

L1 element retrotranspositions have been found to alter expression of genes neighboring the insertion sites, potentially involving them in tumorigenesis and tumor progression. In colorectal cancer (CRC), L1 insertions have been found to target genes with a role in tumorigenesis. Structural changes such as L1 insertions are identifiable by whole genome sequencing (WGS). In this study, we observed frequent somatic L1 retrotranspositions originating from TTC28 using deep coverage WGS data from 92 CRC tumor-normal sample pairs. In two cases the event had targeted NOVA1 gene (p=0.025). In addition, a germline retrotransposition event from TTC28 to GABRA4 was found to be a common polymorphism in the Finnish population. Thus while some events may be tumorigenic, others are likely to be neutral. Our data contradict a recent study where a similar signal in TTC28 was interpreted as a common inactivating translocation. While much work remains to be performed to understand the biological significance of retrotranspositions in cancer, accurate identification of these events is a prerequisite for success.

Exome sequencing reveals frequent inactivating mutations in ARID1A, ARID1B, ARID2 and ARID4A in microsatellite unstable colorectal cancer

ARID1A has been identified as a novel tumor suppressor gene in ovarian cancer and subsequently in various other tumor types. ARID1A belongs to the ARID domain containing gene family, which comprises of 15 genes involved, for example, in transcriptional regulation, proliferation and chromatin remodeling. In this study, we used exome sequencing data to analyze the mutation frequency of all the ARID domain containing genes in 25 microsatellite unstable (MSI) colorectal cancers (CRCs) as a first systematic effort to characterize the mutation pattern of the whole ARID gene family. Genes which fulfilled the selection criteria in this discovery set (mutations in at least 4/25 [16%] samples, including at least one nonsense or splice site mutation) were chosen for further analysis in an independent validation set of 21 MSI CRCs. We found that in addition to ARID1A, which was mutated in 39% of the tumors (18/46), also ARID1B (13%, 6/46), ARID2 (13%, 6/46) and ARID4A (20%, 9/46) were frequently mutated. In all these genes, the mutations were distributed along the entire length of the gene, thus distinguishing them from typical MSI target genes previously described. Our results indicate that in addition to ARID1A, other members of the ARID gene family may play a role in MSI CRC.

Identification of candidate oncogenes in human colorectal cancers with microsatellite instability

Microsatellite instability can be found in approximately 15% of all colorectal cancers. To detect new oncogenes we sequenced the exomes of 25 colorectal tumors and respective healthy colon tissue. Potential mutation hot spots were confirmed in 15 genes; ADAR, DCAF12L2, GLT1D1, ITGA7, MAP1B, MRGPRX4, PSRC1, RANBP2, RPS6KL1, SNCAIP, TCEAL6, TUBB6, WBP5, VEGFB, and ZBTB2; these were validated in 86 tumors with microsatellite instability. ZBTB2, RANBP2, and PSRC1 also were found to contain hot spot mutations in the validation set. The form of ZBTB2 associated with colorectal cancer increased cell proliferation. The mutation hot spots might be used to develop personalized tumor profiling and therapy.

MYBL2 is a sub-haploinsufficient tumor suppressor gene in myeloid malignancy

A common deleted region (CDR) in both myelodysplastic syndromes (MDS) and myeloproliferative neoplasms (MPN) affects the long arm of chromosome 20 and has been predicted to harbor a tumor suppressor gene. Here we show that MYBL2, a gene within the 20q CDR, is expressed at sharply reduced levels in CD34+ cells from most MDS cases (65%; n = 26), whether or not they harbor 20q abnormalities. In a murine competitive reconstitution model, Mybl2 knockdown by RNAi to 20-30% of normal levels in multipotent hematopoietic progenitors resulted in clonal dominance of these 'sub-haploinsufficient' cells, which was reflected in all blood cell lineages. By 6 months post-transplantation, the reconstituted mice had developed a clonal myeloproliferative/myelodysplastic disorder originating from the cells with aberrantly reduced Mybl2 expression. We conclude that downregulation of MYBL2 activity below levels predicted by classical haploinsufficiency underlies the clonal expansion of hematopoietic progenitors in a large fraction of human myeloid malignancies. DOI:http://dx.doi.org/10.7554/eLife.00825.001.

Characterization of uterine leiomyomas by whole-genome sequencing

BACKGROUND

Uterine leiomyomas are benign but affect the health of millions of women. A better understanding of the molecular mechanisms involved may provide clues to the prevention and treatment of these lesions.

METHODS

We performed whole-genome sequencing and gene-expression profiling of 38 uterine leiomyomas and the corresponding myometrium from 30 women.

RESULTS

Identical variants observed in some separate tumor nodules suggested that these nodules have a common origin. Complex chromosomal rearrangements resembling chromothripsis were a common feature of leiomyomas. These rearrangements are best explained by a single event of multiple chromosomal breaks and random reassembly. The rearrangements created tissue-specific changes consistent with a role in the initiation of leiomyoma, such as translocations of the HMGA2 and RAD51B loci and aberrations at the COL4A5-COL4A6 locus, and occurred in the presence of normal TP53 alleles. In some cases, separate events had occurred more than once in single tumor-cell lineages.

CONCLUSIONS

Chromosome shattering and reassembly resembling chromothripsis (a single genomic event that results in focal losses and rearrangements in multiple genomic regions) is a major cause of chromosomal abnormalities in uterine leiomyomas; we propose that tumorigenesis occurs when tissue-specific tumor-promoting changes are formed through these events. Chromothripsis has previously been associated with aggressive cancer; its common occurrence in leiomyomas suggests that it also has a role in the genesis and progression of benign tumors. We observed that multiple separate tumors could be seeded from a single lineage of uterine leiomyoma cells. (Funded by the Academy of Finland Center of Excellence program and others.).

Exome sequencing in diagnostic evaluation of colorectal cancer predisposition in young patients

OBJECTIVE

Early-onset colorectal cancer (CRC), defined here as age of onset less than 40 years, develops frequently in genetically predisposed individuals. Next-generation sequencing is an increasingly available option in the diagnostic workup of suspected hereditary susceptibility, but little is known about the practical feasibility and additional diagnostic yield of the technology in this patient group.

MATERIALS AND METHODS

We analyzed 38 young CRC patients derived from a set of 1514 CRC cases. All 38 tumors had been tested in our laboratory for microsatellite instability (MSI), and Sanger sequencing had been used to screen for MLH1 and MSH2 mutations in MSI cases. Also, gastrointestinal polyposis had been diagnosed clinically and molecularly. Family histories were acquired from national registries. If inherited syndromes had not been diagnosed in routine diagnostic efforts (n = 23), normal tissue DNA was analyzed for mutations in a comprehensive set of high-penetrance genes (MLH1, MSH2, MSH6, PMS2, APC, MUTYH, SMAD4, BMPR1A, LKB1/STK11, and PTEN) by exome sequencing.

RESULTS

CRC predisposition syndromes were confirmed in 42% (16/38) of early-onset CRC patients. Hereditary nonpolyposis colorectal cancer was diagnosed in 12 (32%) patients, familial adenomatous polyposis in three (7.9%), and juvenile polyposis in one (2.6%) patient. Exome sequencing revealed one additional MLH1 mutation. Over half of the patients had advanced cancers (Dukes C or D, 61%, 23/38). The majority of nonsyndromic patients had unaffected first-degree relatives and microsatellite-stable tumors.

CONCLUSIONS

Microsatellite instability positivity or gastrointestinal polyposis characterized all patients with unambiguous highly penetrant germline mutations. In our series, exome sequencing produced little added value in diagnosing the underlying predisposition conditions.

Abstract 3156: New candidate oncogenes discovered in microsatellite unstable colorectal cancer

Around 15% of colorectal cancers (CRCs) show microsatellite instability (MSI). A characteristic high passenger mutation load has discouraged systematic mutation screens in MSI CRCs. To systematically search for novel MSI CRC oncogenes, the exomes of 25 MSI CRC-normal pairs were sequenced as the discovery set. Observed hot spots were confirmed by Sanger sequencing and further validated in a set of 86 MSI CRCs. Mutational hot spots were confirmed in 15 genes and three genes showed additional hot spot mutations in the validation set. These three sites were highly conserved across species. The hot spot sites of these three genes were screened in 75 microsatellite stable CRCs and 12 MSI CRC cell lines with negative results. Next, we analyzed the subcellular localization of wild-type and mutant proteins. The findings of this study support the idea that cancer genomes are heterogeneous and characterized by few, frequently mutated “mountains” (e.g. BRAF and KRAS) and numerous, less frequently mutated “hills”. The identified mutational hot spots may prove important in developing personalized tumor profiling and therapy.

Citation Format: Alexandra E. Gylfe, Johanna Kondelin, Mikko Turunen, Heikki Ristolainen, Riku Katainen, Esa Pitkänen, Eevi Kaasinen, Ville Rantanen, Tomas Tanskanen, Heli J. Lehtonen, Kimmo Palin, Minna Taipale, Jussi Taipale, Laura Renkonen-Sinisalo, Heikki Järvinen, Jan Böhm, Jukka-Pekka Mecklin, Ari Ristimäki, Sari Tuupanen, Auli Karhu, Outi Kilpivaara, Pia Vahteristo, Lauri A. Aaltonen. New candidate oncogenes discovered in microsatellite unstable colorectal cancer. [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 3156. doi:10.1158/1538-7445.AM2013-3156

Diagnostic cancer genome sequencing and the contribution of germline variants

Whole-genome sequencing (WGS) is revolutionizing medical research and has the potential to serve as a powerful and cost-effective diagnostic tool in the management of cancer. We review the progress to date in the use of WGS to reveal how germline variants and mutations may be associated with cancer. We use colorectal cancer as an example of how the current level of knowledge can be translated into predictions of predisposition. We also address challenges in the clinical implementation of the variants in germline DNA identified through cancer genome sequencing. We call for the international development of standards to facilitate the clinical use of germline information arising from diagnostic cancer genome sequencing.

Nationwide registry-based analysis of cancer clustering detects strong familial occurrence of Kaposi sarcoma

Many cancer predisposition syndromes are rare or have incomplete penetrance, and traditional epidemiological tools are not well suited for their detection. Here we have used an approach that employs the entire population based data in the Finnish Cancer Registry (FCR) for analyzing familial aggregation of all types of cancer, in order to find evidence for previously unrecognized cancer susceptibility conditions. We performed a systematic clustering of 878,593 patients in FCR based on family name at birth, municipality of birth, and tumor type, diagnosed between years 1952 and 2011. We also estimated the familial occurrence of the tumor types using cluster score that reflects the proportion of patients belonging to the most significant clusters compared to all patients in Finland. The clustering effort identified 25,910 birth name-municipality based clusters representing 183 different tumor types characterized by topography and morphology. We produced information about familial occurrence of hundreds of tumor types, and many of the tumor types with high cluster score represented known cancer syndromes. Unexpectedly, Kaposi sarcoma (KS) also produced a very high score (cluster score 1.91, p-value <0.0001). We verified from population records that many of the KS patients forming the clusters were indeed close relatives, and identified one family with five affected individuals in two generations and several families with two first degree relatives. Our approach is unique in enabling systematic examination of a national epidemiological database to derive evidence of aberrant familial aggregation of all tumor types, both common and rare. It allowed effortless identification of families displaying features of both known as well as potentially novel cancer predisposition conditions, including striking familial aggregation of KS. Further work with high-throughput methods should elucidate the molecular basis of the potentially novel predisposition conditions found in this study.

HOXB13 G84E mutation in Finland: population-based analysis of prostate, breast, and colorectal cancer risk

BACKGROUND

A recently identified germline mutation G84E in HOXB13 was shown to increase the risk of prostate cancer. In a family-based analysis by The International Consortium for Prostate Cancer Genetics (ICPCG), the G84E mutation was most prevalent in families from the Nordic countries of Finland (22.4%) and Sweden (8.2%).

METHODS

To further investigate the importance of G84E in the Finns, we determined its frequency in more than 4,000 prostate cancer cases and 5,000 controls. In addition, 986 breast cancer and 442 colorectal cancer (CRC) cases were studied. Genotyping was conducted using TaqMan, MassARRAY iPLEX, and sequencing. Statistical analyses were conducted using Fisher exact test, and overall survival was analyzed using Cox modeling.

RESULTS

The frequency of the G84E mutation was significantly higher among patients with prostate cancer and highest among patients with a family history of the disease, hereditary prostate cancer [8.4% vs. 1.0% in controls; OR 8.8; 95% confidence interval (CI), 4.9-15.7]. The mutation contributed significantly to younger age (≤55 years) at onset and high prostate-specific antigen (PSA; ≥20 ng/mL) at diagnosis. An association with increased prostate cancer risk in patients with prior benign prostate hyperplasia (BPH) diagnosis was also revealed. No statistically significant evidence for a contribution in CRC risk was detected, but a suggestive role for the mutation was observed in familial BRCA1/2-negative breast cancer.

CONCLUSIONS

These findings confirm an increased cancer risk associated with the G84E mutation in the Finnish population, particularly for early-onset prostate cancer and cases with substantially elevated PSA.

IMPACT

This study confirms the overall importance of the HOXB13 G84E mutation in prostate cancer susceptibility.

Somatic MED12 mutations in uterine leiomyosarcoma and colorectal cancer

Background:

Mediator complex participates in transcriptional regulation by connecting regulatory DNA sequences to the RNA polymerase II initiation complex. Recently, we discovered through exome sequencing that as many as 70% of uterine leiomyomas harbour specific mutations in exon 2 of mediator complex subunit 12 (MED12). In this work, we examined the role of MED12 exon 2 mutations in other tumour types.

Methods:

The frequency of MED12 exon 2 mutations was analysed in altogether 1158 tumours by direct sequencing. The tumour spectrum included mesenchymal tumours (extrauterine leiomyomas, endometrial polyps, lipomas, uterine leiomyosarcomas, other sarcomas, gastro-intestinal stromal tumours), hormone-dependent tumours (breast and ovarian cancers), haematological malignancies (acute myeloid leukaemias, acute lymphoid leukaemias, myeloproliferative neoplasms), and tumours associated with abnormal Wnt-signalling (colorectal cancers (CRC)).

Results:

Five somatic alterations were observed: three in uterine leiomyosarcomas (3/41, 7% Gly44Ser, Ala38_Leu39ins7, Glu35_Leu36delinsVal), and two in CRC (2/392, 0.5% Gly44Cys, Ala67Val).

Conclusion:

Somatic MED12 exon 2 mutations were observed in uterine leiomyosarcomas, suggesting that a subgroup of these malignant tumours may develop from a leiomyoma precursor. Mutations in CRC samples indicate that MED12 may, albeit rarely, contribute to CRC tumorigenesis.

Breast tumors from CHEK2 1100delC-mutation carriers: genomic landscape and clinical implications

INTRODUCTION

Checkpoint kinase 2 (CHEK2) is a moderate penetrance breast cancer risk gene, whose truncating mutation 1100delC increases the risk about twofold. We investigated gene copy-number aberrations and gene-expression profiles that are typical for breast tumors of CHEK2 1100delC-mutation carriers.

METHODS

In total, 126 breast tumor tissue specimens including 32 samples from patients carrying CHEK2 1100delC were studied in array-comparative genomic hybridization (aCGH) and gene-expression (GEX) experiments. After dimensionality reduction with CGHregions R package, CHEK2 1100delC-associated regions in the aCGH data were detected by the Wilcoxon rank-sum test. The linear model was fitted to GEX data with R package limma. Genes whose expression levels were associated with CHEK2 1100delC mutation were detected by the bayesian method.

RESULTS

We discovered four lost and three gained CHEK2 1100delC-related loci. These include losses of 1p13.3-31.3, 8p21.1-2, 8p23.1-2, and 17p12-13.1 as well as gains of 12q13.11-3, 16p13.3, and 19p13.3. Twenty-eight genes located on these regions showed differential expression between CHEK2 1100delC and other tumors, nominating them as candidates for CHEK2 1100delC-associated tumor-progression drivers. These included CLCA1 on 1p22 as well as CALCOCO1, SBEM, and LRP1 on 12q13. Altogether, 188 genes were differentially expressed between CHEK2 1100delC and other tumors. Of these, 144 had elevated and 44, reduced expression levels.Our results suggest the WNT pathway as a driver of tumorigenesis in breast tumors of CHEK2 1100delC-mutation carriers and a role for the olfactory receptor protein family in cancer progression. Differences in the expression of the 188 CHEK2 1100delC-associated genes divided breast tumor samples from three independent datasets into two groups that differed in their relapse-free survival time.

CONCLUSIONS

We have shown that copy-number aberrations of certain genomic regions are associated with CHEK2 mutation 1100delC. On these regions, we identified potential drivers of CHEK2 1100delC-associated tumorigenesis, whose role in cancer progression is worth investigating. Furthermore, poorer survival related to the CHEK2 1100delC gene-expression signature highlights pathways that are likely to have a role in the development of metastatic disease in carriers of the CHEK2 1100delC mutation.

IDH1 and IDH2 mutation studies in 1473 patients with chronic-, fibrotic- or blast-phase essential thrombocythemia, polycythemia vera or myelofibrosis

In a multi-institutional collaborative project, 1473 patients with myeloproliferative neoplasms (MPN) were screened for isocitrate dehydrogenase 1 (IDH1)/IDH2 mutations: 594 essential thrombocythemia (ET), 421 polycythemia vera (PV), 312 primary myelofibrosis (PMF), 95 post-PV/ET MF and 51 blast-phase MPN. A total of 38 IDH mutations (18 IDH1-R132, 19 IDH2-R140 and 1 IDH2-R172) were detected: 5 (0.8%) ET, 8 (1.9%) PV, 13 (4.2%) PMF, 1 (1%) post-PV/ET MF and 11 (21.6%) blast-phase MPN (P<0.01). Mutant IDH was documented in the presence or absence of JAK2, MPL and TET2 mutations, with similar mutational frequencies. However, IDH-mutated patients were more likely to be nullizygous for JAK2 46/1 haplotype, especially in PMF (P=0.04), and less likely to display complex karyotype, in blast-phase disease (P<0.01). In chronic-phase PMF, JAK2 46/1 haplotype nullizygosity (P<0.01; hazard ratio (HR) 2.9, 95% confidence interval (CI) 1.7-5.2), but not IDH mutational status (P=0.55; HR 1.3, 95% CI 0.5-3.4), had an adverse effect on survival. This was confirmed by multivariable analysis. In contrast, in both blast-phase PMF (P=0.04) and blast-phase MPN (P=0.01), the presence of an IDH mutation predicted worse survival. The current study clarifies disease- and stage-specific IDH mutation incidence and prognostic relevance in MPN and provides additional evidence for the biological effect of distinct JAK2 haplotypes.

Genomic subtypes of breast cancer identified by array-comparative genomic hybridization display distinct molecular and clinical characteristics

INTRODUCTION

Breast cancer is a profoundly heterogeneous disease with respect to biologic and clinical behavior. Gene-expression profiling has been used to dissect this complexity and to stratify tumors into intrinsic gene-expression subtypes, associated with distinct biology, patient outcome, and genomic alterations. Additionally, breast tumors occurring in individuals with germline BRCA1 or BRCA2 mutations typically fall into distinct subtypes.

METHODS

We applied global DNA copy number and gene-expression profiling in 359 breast tumors. All tumors were classified according to intrinsic gene-expression subtypes and included cases from genetically predisposed women. The Genomic Identification of Significant Targets in Cancer (GISTIC) algorithm was used to identify significant DNA copy-number aberrations and genomic subgroups of breast cancer.

RESULTS

We identified 31 genomic regions that were highly amplified in > 1% of the 359 breast tumors. Several amplicons were found to co-occur, the 8p12 and 11q13.3 regions being the most frequent combination besides amplicons on the same chromosomal arm. Unsupervised hierarchical clustering with 133 significant GISTIC regions revealed six genomic subtypes, termed 17q12, basal-complex, luminal-simple, luminal-complex, amplifier, and mixed subtypes. Four of them had striking similarity to intrinsic gene-expression subtypes and showed associations to conventional tumor biomarkers and clinical outcome. However, luminal A-classified tumors were distributed in two main genomic subtypes, luminal-simple and luminal-complex, the former group having a better prognosis, whereas the latter group included also luminal B and the majority of BRCA2-mutated tumors. The basal-complex subtype displayed extensive genomic homogeneity and harbored the majority of BRCA1-mutated tumors. The 17q12 subtype comprised mostly HER2-amplified and HER2-enriched subtype tumors and had the worst prognosis. The amplifier and mixed subtypes contained tumors from all gene-expression subtypes, the former being enriched for 8p12-amplified cases, whereas the mixed subtype included many tumors with predominantly DNA copy-number losses and poor prognosis.

CONCLUSIONS

Global DNA copy-number analysis integrated with gene-expression data can be used to dissect the complexity of breast cancer. This revealed six genomic subtypes with different clinical behavior and a striking concordance to the intrinsic subtypes. These genomic subtypes may prove useful for understanding the mechanisms of tumor development and for prognostic and treatment prediction purposes.

A germline JAK2 SNP is associated with predisposition to the development of JAK2(V617F)-positive myeloproliferative neoplasms

Polycythemia vera, essential thrombocythemia and primary myelofibrosis are myeloproliferative neoplasms (MPN) characterized by multilineage clonal hematopoiesis. Given that the identical somatic activating mutation in the JAK2 tyrosine kinase gene (JAK2(V617F)) is observed in most individuals with polycythemia vera, essential thrombocythemia and primary myelofibrosis, there likely are additional genetic events that contribute to the pathogenesis of these phenotypically distinct disorders. Moreover, family members of individuals with MPN are at higher risk for the development of MPN, consistent with the existence of MPN predisposition loci. We hypothesized that germline variation contributes to MPN predisposition and phenotypic pleiotropy. Genome-wide analysis identified an allele in the JAK2 locus (rs10974944) that predisposes to the development of JAK2(V617F)-positive MPN, as well as three previously unknown MPN modifier loci. We found that JAK2(V617F) is preferentially acquired in cis with the predisposition allele. These data suggest that germline variation is an important contributor to MPN phenotype and predisposition.

Family history, genetic testing, and clinical risk prediction: pooled analysis of CHEK2 1100delC in 1,828 bilateral breast cancers and 7,030 controls

If breast cancers arise independently in each breast the odds ratio (OR) for bilateral breast cancer for carriers of CHEK2 1100delC should be approximately 5.5, the square of the reported OR for a first primary (OR, 2.34). In the subset of bilateral cases with one or more affected relatives, the predicted carrier OR should be approximately 9. We have tested these predictions in a pooled set of 1,828 cases with 2 primaries and 7,030 controls from 8 studies. The second primary OR for CHEK2 1100delC carriers was 6.43 (95% confidence interval, 4.33-9.56; P < 0.0001), significantly greater than the published estimate for a first primary (P < 0.001) but consistent with its square. The predicted increase in carrier OR with increasing numbers of affected relatives was seen using bilateral cases from the UK (P(trend) = 0.0003) and Finland (P(trend) = 0.37), although not using those from the Netherlands and Russia (P = 0.001 for heterogeneity between countries). Based on a standard genetic model, we predict lifetime risks for CHEK2 1100delC carrier and noncarrier daughters of bilateral breast cancer cases of 37% and 18%, respectively. Our results imply that clinical management of the daughter of a woman with bilateral breast cancer should depend on her CHEK2 1100delC carrier status. This and other moderate penetrance breast cancer susceptibility alleles, together with family history data, will thus identify increasing numbers of women at potentially very high risk. Before such predictions are accepted by clinical geneticists, however, further population-based evidence is needed on the effect of CHEK2 1100delC and other moderate penetrance alleles in women with a family history of breast cancer.

NAD(P)H:quinone oxidoreductase 1 NQO1*2 genotype (P187S) is a strong prognostic and predictive factor in breast cancer

NQO1 guards against oxidative stress and carcinogenesis and stabilizes p53. We find that a homozygous common missense variant (NQO1(*)2, rs1800566(T), NM_000903.2:c.558C>T) that disables NQO1 strongly predicts poor survival among two independent series of women with breast cancer (P = 0.002, N = 1,005; P = 0.005, N = 1,162), an effect particularly evident after anthracycline-based adjuvant chemotherapy with epirubicin (P = 7.52 x 10(-6)) and in p53-aberrant tumors (P = 6.15 x 10(-5)). Survival after metastasis was reduced among NQO1(*)2 homozygotes, further implicating NQO1 deficiency in cancer progression and treatment resistance. Consistently, response to epirubicin was impaired in NQO1(*)2-homozygous breast carcinoma cells in vitro, reflecting both p53-linked and p53-independent roles of NQO1. We propose a model of defective anthracycline response in NQO1-deficient breast tumors, along with increased genomic instability promoted by elevated reactive oxygen species (ROS), and suggest that the NQO1 genotype is a prognostic and predictive marker for breast cancer.

Prognostic role of HuR in hereditary breast cancer

PURPOSE

HuR is an mRNA-binding protein that enhances the stability of certain transcripts and can regulate their translation. Elevated cytoplasmic expression of HuR protein has been linked to carcinogenesis and is associated with reduced survival in breast, ovarian, and gastric adenocarcinomas.

EXPERIMENTAL DESIGN

Here, we have explored the relevance of HuR in familial breast cancer. Tumor samples were collected from patients with identified BRCA1 (n = 51) or BRCA2 (n = 47) mutations or familial non-BRCA1/2 cases (n = 525), and analyzed by immunohistochemistry.

RESULTS

Among familial non-BRCA1/2 breast cancer patients, cytoplasmic HuR protein expression was present in 39.4% of the cases and was associated with estrogen receptor negativity, progesterone receptor negativity, p53 positivity, high tumor grade, and ductal type of the tumor. In multivariate analysis, cytoplasmic HuR expression was an independent marker of reduced survival in the non-BRCA1/2 group along with tumor size >2 cm, lymph node metastasis, and high histologic grade. In patients with BRCA1 or BRCA2 mutations, cytoplasmic HuR expression was more frequent (62.7% for BRCA1 and 61.7% for BRCA2) than in the non-BRCA1/2 group, but in BRCA-mutated subgroups cytoplasmic HuR expression did not associate with survival.

CONCLUSIONS

Our results show that HuR is an important prognostic factor in familial breast cancer patients and may contribute to carcinogenesis in this disease.

Comprehensive analysis of NuMA variation in breast cancer

Background

A recent genome wide case-control association study identified NuMA region on 11q13 as a candidate locus for breast cancer susceptibility. Specifically, the variant Ala794Gly was suggested to be associated with increased risk of breast cancer.

Methods

In order to evaluate the NuMa gene for breast cancer susceptibility, we have here screened the entire coding region and exon-intron boundaries of NuMa in 92 familial breast cancer patients and constructed haplotypes of the identified variants. Five missense variants were further screened in 341 breast cancer cases with a positive family history and 368 controls. We examined the frequency of Ala794Gly in an extensive series of familial (n = 910) and unselected (n = 884) breast cancer cases and controls (n = 906), with a high power to detect the suggested breast cancer risk. We also tested if the variant is associated with histopathologic features of breast tumors.

Results

Screening of NuMA resulted in identification of 11 exonic variants and 12 variants in introns or untranslated regions. Five missense variants that were further screened in breast cancer cases with a positive family history and controls, were each carried on a unique haplotype. None of the variants, or the haplotypes represented by them, was associated with breast cancer risk although due to low power in this analysis, very low risk alleles may go unrecognized. The NuMA Ala794Gly showed no difference in frequency in the unselected breast cancer case series or familial case series compared to control cases. Furthermore, Ala794Gly did not show any significant association with histopathologic characteristics of the tumors, though Ala794Gly was slightly more frequent among unselected cases with lymph node involvement.

Conclusion

Our results do not support the role of NuMA variants as breast cancer susceptibility alleles.

Basal cytokeratins in breast tumours among BRCA1, BRCA2 and mutation-negative breast cancer families

INTRODUCTION

Finding new immunohistochemical markers that are specific to hereditary breast cancer could help us to select candidates for BRCA1/BRCA2 mutation testing and to understand the biological pathways of tumour development.

METHODS

Using breast cancer tumour microarrays, immunohistochemical expression of cytokeratin (CK)-5/6, CK-14 and CK-17 was evaluated in breast tumours from BRCA1 families (n = 46), BRCA2 families (n = 40), non-BRCA1/BRCA2 families (n = 358) and familial breast cancer patients with one first-degree relative affected by breast or ovarian cancer (n = 270), as well as from patients with sporadic breast cancer (n = 364). Staining for CK-5/6, CK-14 and CK-17 was compared between these groups and correlated with other clinical and histological factors.

RESULTS

CK-5/6, CK-14 and CK-17 were detected mostly among oestrogen receptor (ER)-negative, progesterone receptor (PR)-negative and high-grade tumours. We found the highest percentages of samples positive for these CKs among ER-negative/HER2-negative tumours. In univariate analysis, CK-14 was significantly associated with tumours from BRCA1 (39%; P < 0.0005), BRCA2 (27%; P = 0.011), and non-BRCA1/BRCA2 (21%; P < 0.005) families, as compared with sporadic tumours (10%). However, in multivariate analysis, CKs were not found to be independently associated with BRCA1 or BRCA2 mutation status, and the most effective predictors of BRCA1 mutations were age at onset, HER2 status, and either ER or PR status.

CONCLUSION

Although our study confirms that basal CKs can help to identify BRCA1 mutation carriers, this effect was weaker than previously suggested and CKs did not independently predict BRCA1 mutation either from sporadic or familial breast cancer cases. The most effective, independent predictors of BRCA1 mutations were age at onset, HER2 status, and either ER or PR status, as compared with sporadic or non-BRCA1/BRCA2 cancers.

Combined effect of CCND1 and COMT polymorphisms and increased breast cancer risk

BACKGROUND

Estrogens are crucial tumorigenic hormones, which impact the cell growth and proliferation during breast cancer development. Estrogens are metabolized by a series of enzymes including COMT, which converts catechol estrogens into biologically non-hazardous methoxyestrogens. Several studies have also shown the relationship between estrogen and cell cycle progression through activation of CCND1 transcription.

METHODS

In this study, we have investigated the independent and the combined effects of commonly occurring CCND1 (Pro241Pro, A870G) and COMT (Met108/158Val) polymorphisms to breast cancer risk in two independent Caucasian populations from Ontario (1228 breast cancer cases and 719 population controls) and Finland (728 breast cancer cases and 687 population controls). Both COMT and CCND1 polymorphisms have been previously shown to impact on the enzymatic activity of the coded proteins.

RESULTS

Here, we have shown that the high enzymatic activity genotype of CCND1High (AA) was associated with increased breast cancer risk in both the Ontario [OR: 1.3, 95%CI (1.0-1.69)] and the Finland sample [OR: 1.4, 95%CI (1.01-1.84)]. The heterozygous COMTMedium (MetVal) and the high enzymatic activity of COMTHigh (ValVal) genotype was also associated with breast cancer risk in Ontario cases, [OR: 1.3, 95%CI (1.07-1.68)] and [OR: 1.4, 95%CI (1.07-1.81)], respectively. However, there was neither a statistically significant association nor increased trend of breast cancer risk with COMTHigh (ValVal) genotypes in the Finland cases [OR: 1.0, 95%CI (0.73-1.39)]. In the combined analysis, the higher activity alleles of the COMT and CCND1 is associated with increased breast cancer risk in both Ontario [OR: 2.22, 95%CI (1.49-3.28)] and Finland [OR: 1.73, 95%CI (1.08-2.78)] populations studied. The trend test was statistically significant in both the Ontario and Finland populations across the genotypes associated with increasing enzymatic activity.

CONCLUSION

Using two independent Caucasian populations, we have shown a stronger combined effect of the two commonly occurring CCND1 and COMT genotypes in the context of breast cancer predisposition.

RAD51 135G-->C modifies breast cancer risk among BRCA2 mutation carriers: results from a combined analysis of 19 studies

RAD51 is an important component of double-stranded DNA-repair mechanisms that interacts with both BRCA1 and BRCA2. A single-nucleotide polymorphism (SNP) in the 5' untranslated region (UTR) of RAD51, 135G-->C, has been suggested as a possible modifier of breast cancer risk in BRCA1 and BRCA2 mutation carriers. We pooled genotype data for 8,512 female mutation carriers from 19 studies for the RAD51 135G-->C SNP. We found evidence of an increased breast cancer risk in CC homozygotes (hazard ratio [HR] 1.92 [95% confidence interval {CI} 1.25-2.94) but not in heterozygotes (HR 0.95 [95% CI 0.83-1.07]; P=.002, by heterogeneity test with 2 degrees of freedom [df]). When BRCA1 and BRCA2 mutation carriers were analyzed separately, the increased risk was statistically significant only among BRCA2 mutation carriers, in whom we observed HRs of 1.17 (95% CI 0.91-1.51) among heterozygotes and 3.18 (95% CI 1.39-7.27) among rare homozygotes (P=.0007, by heterogeneity test with 2 df). In addition, we determined that the 135G-->C variant affects RAD51 splicing within the 5' UTR. Thus, 135G-->C may modify the risk of breast cancer in BRCA2 mutation carriers by altering the expression of RAD51. RAD51 is the first gene to be reliably identified as a modifier of risk among BRCA1/2 mutation carriers.

The DNA damage signalling kinase ATM is aberrantly reduced or lost in BRCA1/BRCA2-deficient and ER/PR/ERBB2-triple-negative breast cancer

The ataxia-telangiectasia-mutated (ATM) kinase is a key transducer of DNA damage signals within the genome maintenance machinery and a tumour suppressor whose germline mutations predispose to familial breast cancer. ATM signalling is constitutively activated in early stages of diverse types of human malignancies and cell culture models in response to oncogene-induced DNA damage providing a barrier against tumour progression. As BRCA1 and BRCA2 are also components of the genome maintenance network and their mutations predispose to breast cancer, we have examined the ATM expression in human breast carcinomas of BRCA1/2 mutation carriers, sporadic cases and familial non-BRCA1/2 patients. Our results show that ATM protein expression is aberrantly reduced more frequently among BRCA1 (33%; P=0.0003) and BRCA2 (30%; P=0.0009) tumours than in non-BRCA1/2 tumours (10.7%). Furthermore, the non-BRCA1/2 tumours with reduced ATM expression were more often estrogen receptor (ER) negative (P=0.0002), progesterone receptor (PR) negative (P=0.004) and were of higher grade (P=0.0004). In our series of 1013 non-BRCA1/2 cases, ATM was more commonly deficient (20%; P=0.0006) and p53 was overabundant (47%; P<0.0000000001) among the difficult-to-treat ER/PR/ERBB2-triple-negative subset of tumours compared with cases that expressed at least one of these receptors (10 and 16% of aberrant ATM and p53, respectively). We propose a model of 'conditional haploinsufficiency' for BRCA1/2 under conditions of enhanced DNA damage in precancerous lesions resulting in more robust activation and hence increased selection for inactivation or loss of ATM in tumours of BRCA1/2 mutation carriers, with implications for genomic instability and curability of diverse subsets of human breast cancer.

ATM variants and cancer risk in breast cancer patients from Southern Finland

Background

Individuals heterozygous for germline ATM mutations have been reported to have an increased risk for breast cancer but the role for ATM genetic variants for breast cancer risk has remained unclear. Recently, a common ATM variant, ATMivs38 -8T>C in cis with the ATMex39 5557G>A (D1853N) variant, was suggested to associate with bilateral breast cancer among familial breast cancer patients from Northern Finland. We have here evaluated the 5557G>A and ivs38-8T>C variants in an extensive case-control association analysis. We also aimed to investigate whether there are other ATM mutations or variants contributing to breast cancer risk in our population.

Methods

Two common ATM variants, 5557G>A and ivs38-8T>C, previously suggested to associate with bilateral breast cancer, were genotyped in an extensive set of 786 familial and 884 unselected breast cancer cases as well as 708 healthy controls. We also screened the entire coding region and exon-intron boundaries of the ATM gene in 47 familial breast cancer patients and constructed haplotypes of the patients. The identified variants were also evaluated for increased breast cancer risk among additional breast cancer cases and controls.

Results

Neither of the two common variants, 5557G>A and ivs38-8T>C, nor any haplotype containing them, was significantly associated with breast cancer risk, bilateral breast cancer or multiple primary cancers in any of the patient groups or subgoups. Three rare missense alterations and one intronic change were each found in only one patient of over 250 familial patients studied and not among controls. The fourth missense alteration studied further was found with closely similar frequencies in over 600 familial cases and controls.

Conclusion

Altogether, our results suggest very minor effect, if any, of ATM genetic variants on familial breast cancer in Southern Finland. Our results do not support association of the 5557G>A or ivs38-8T>C variant with increased breast cancer risk or with bilateral breast cancer.

CHEK2 I157T associates with familial and sporadic colorectal cancer

BACKGROUND

Recently, a functionally defective CHEK2 variant I157T has been proposed to associate with an increased risk of several types of cancer. We investigated the CHEK2 I157T variant for colorectal cancer (CRC) predisposition in a large population based study including a significant number of familial CRC cases.

METHODS

We screened the CHEK2 I157T variant in a population based series of 1042 Finnish CRC patients using restriction fragment length polymorphism. Mutation status was studied for correlation with clinical characteristics and family history of CRC and other cancers.

RESULTS

The frequency of CHEK2 I157T was significantly higher in CRC patients (7.8%, 76/972) than in healthy population controls (5.3%, 100/1885) (OR = 1.5, 95% CI 1.1 to 2.1, p = 0.008). The significant association of CHEK2 I157T with CRC was observed among patients with (10.4%, 14/135) and without (7.4%, 62/837) a family history of CRC (OR = 2.1, 95% CI 1.1 to 3.7, p = 0.01; OR = 1.4, 95% CI 1.0 to 2.0, p = 0.03; respectively). A trend towards higher variant frequency was also noted among patients with multiple primary tumours and a family history of any cancer.

CONCLUSIONS

CHEK2 I157T associates with an increased risk of CRC: the association was observed both among familial and sporadic CRC patients. Furthermore, the higher frequency of I157T among patients with multiple primary tumours as well as those with a family history of any cancer supports a role for CHEK2 I157T as a susceptibility allele for multiple cancer types.