Mads Thomassen

Gene expression meta-analysis identifies chromosomal regions and candidate genes involved in breast cancer metastasis

Breast cancer cells exhibit complex karyotypic alterations causing deregulation of numerous genes. Some of these genes are probably causal for cancer formation and local growth whereas others are causal for the various steps of metastasis. In a fraction of tumors deregulation of the same genes might be caused by epigenetic modulations, point mutations or the influence of other genes. We have investigated the relation of gene expression and chromosomal position, using eight datasets including more than 1200 breast tumors, to identify chromosomal regions and candidate genes possibly causal for breast cancer metastasis. By use of “Gene Set Enrichment Analysis” we have ranked chromosomal regions according to their relation to metastasis. Overrepresentation analysis identified regions with increased expression for chromosome 1q41–42, 8q24, 12q14, 16q22, 16q24, 17q12–21.2, 17q21–23, 17q25, 20q11, and 20q13 among metastasizing tumors and reduced gene expression at 1p31–21, 8p22–21, and 14q24. By analysis of genes with extremely imbalanced expression in these regions we identified DIRAS3 at 1p31, PSD3, LPL, EPHX2 at 8p21–22, and FOS at 14q24 as candidate metastasis suppressor genes. Potential metastasis promoting genes includes RECQL4 at 8q24, PRMT7 at 16q22, GINS2 at 16q24, and AURKA at 20q13.

On the origin and diffusion of BRCA1 c.5266dupC (5382insC) in European populations

The BRCA1 mutation c.5266dupC was originally described as a founder mutation in the Ashkenazi Jewish (AJ) population. However, this mutation is also present at appreciable frequency in several European countries, which raises intriguing questions about the origins of the mutation. We genotyped 245 carrier families from 14 different population groups (Russian, Latvian, Ukrainian, Czech, Slovak, Polish, Danish, Dutch, French, German, Italian, Greek, Brazilian and AJ) for seven microsatellite markers and confirmed that all mutation carriers share a common haplotype from a single founder individual. Using a maximum likelihood method that allows for both recombination and mutational events of marker loci, we estimated that the mutation arose some 1800 years ago in either Scandinavia or what is now northern Russia and subsequently spread to the various populations we genotyped during the following centuries, including the AJ population. Age estimates and the molecular evolution profile of the most common linked haplotype in the carrier populations studied further suggest that c.5266dupC likely entered the AJ gene pool in Poland approximately 400–500 years ago. Our results illustrate that (1) BRCA1 c.5266dupC originated from a single common ancestor and was a common European mutation long before becoming an AJ founder mutation and (2) the mutation is likely present in many additional European countries where genetic screening of BRCA1 may not yet be common practice.

The OncoArray Consortium: a Network for Understanding the Genetic Architecture of Common Cancers.

Background: Common cancers develop through a multistep process often including inherited susceptibility. Collaboration among multiple institutions, and funding from multiple sources, has allowed the development of an inexpensive genotyping microarray, the OncoArray. The array includes a genome-wide backbone, comprising 230,000 SNPs tagging most common genetic variants, together with dense mapping of known susceptibility regions, rare variants from sequencing experiments, pharmacogenetic markers, and cancer-related traits. Methods: The OncoArray can be genotyped using a novel technology developed by Illumina to facilitate efficient genotyping. The consortium developed standard approaches for selecting SNPs for study, for quality control of markers, and for ancestry analysis. The array was genotyped at selected sites and with prespecified replicate samples to permit evaluation of genotyping accuracy among centers and by ethnic background. Results: The OncoArray consortium genotyped 447,705 samples. A total of 494,763 SNPs passed quality control steps with a sample success rate of 97% of the samples. Participating sites performed ancestry analysis using a common set of markers and a scoring algorithm based on principal components analysis. Conclusions: Results from these analyses will enable researchers to identify new susceptibility loci, perform fine-mapping of new or known loci associated with either single or multiple cancers, assess the degree of overlap in cancer causation and pleiotropic effects of loci that have been identified for disease-specific risk, and jointly model genetic, environmental, and lifestyle-related exposures. Impact: Ongoing analyses will shed light on etiology and risk assessment for many types of cancer. Cancer Epidemiol Biomarkers Prev; 26(1); 126–35. ©2016 AACR.

BRCA1 R1699Q variant displaying ambiguous functional abrogation confers intermediate breast and ovarian cancer risk

Background Clinical classification of rare sequence changes identified in the breast cancer susceptibility genes BRCA1 and BRCA2 is essential for appropriate genetic counselling of individuals carrying these variants. We previously showed that variant BRCA1 c.5096G>A p.Arg1699Gln in the BRCA1 transcriptional transactivation domain demonstrated equivocal results from a series of functional assays, and proposed that this variant may confer low to moderate risk of cancer. Methods Measures of genetic risk (report of family history, segregation) were assessed for 68 BRCA1 c.5096G>A p.Arg1699Gln (R1699Q) families recruited through family cancer clinics, comparing results with 34 families carrying the previously classified pathogenic BRCA1 c.5095C>T p.Arg1699Trp (R1699W) mutation at the same residue, and to 243 breast cancer families with no BRCA1 pathogenic mutation (BRCA-X). Results Comparison of BRCA1 carrier prediction scores of probands using the BOADICEA risk prediction tool revealed that BRCA1 c.5096G>A p.Arg1699Gln variant carriers had family histories that were less ‘BRCA1-like’ than BRCA1 c.5095C>T p.Arg1699Trp mutation carriers (p<0.00001), but more ‘BRCA1-like’ than BRCA-X families (p=0.0004). Further, modified segregation analysis of the subset of 30 families with additional genotyping showed that BRCA1 c.5096G >A p.Arg1699Gln had reduced penetrance compared with the average truncating BRCA1 mutation penetrance (p=0.0002), with estimated cumulative risks to age 70 of breast or ovarian cancer of 24%. Conclusions Our results provide substantial evidence that the BRCA1 c.5096G>A p.Arg1699Gln (R1699Q) variant, demonstrating ambiguous functional deficiency across multiple assays, is associated with intermediate risk of breast and ovarian cancer, highlighting challenges for risk modelling and clinical management of patients of this and other potential moderate-risk variants.

Myeloma cell expression of 10 candidate genes for osteolytic bone disease. Only overexpression of DKK1 correlates with clinical bone involvement at diagnosis

Osteolytic bone disease (OBD) in multiple myeloma (MM) is caused by interactions between MM cells and the bone marrow microenvironment and is characterized by increased osteoclastic bone resorption and decreased osteoblastic bone formation. Recently, the role of osteoblast inhibition has come into focus, especially the possible role of overexpression of DKK1, an inhibitor of the Wnt signalling pathway. Further, CKS2, PSME2 and DHFR have also been reported as candidate genes for OBD. We studied the gene expression by quantitative reverse transcription polymerase chain reaction of TNFSF11 (RANKL), TNFSF11A (RANK), TNFRSF11B (OPG), CCL3 (MIP1A), CCL4 (MIP1B), PTHR1 (PTHrp), DKK1, CKS2, PSME2 and DHFR in purified, immunophenotypic FACS‐sorted plasma cells from 171 newly diagnosed MM patients, 20 patients with monoclonal gammopathy of undetermined significance and 12 controls. The gene expressions of the analysed genes were correlated with radiographically assessed OBD. Only overexpression of DKK1 was correlated to the degree of OBD. Myeloma cells did not express TNFSF11A, TNFSF11, or TNFRSF11B, and very rarely expressed CCL3 and PTHR11. CCL4, CKS2, PSME2 and DHFR were variably expressed, but the expression of these genes showed no correlation with OBD. In contrast, loss of PSME2 expression in MM plasma cells was significantly correlated with OBD.

Gene expression profiling with principal component analysis depicts the biological continuum from essential thrombocythemia over polycythemia vera to myelofibrosis.

The recent discovery of the Janus activating kinase 2 V617F mutation in most patients with polycythemia vera (PV) and half of those with essential thrombocythemia (ET) and primary myelofibrosis (PMF) has favored the hypothesis of a biological continuum from ET over PV to PMF. We performed gene expression profiling of whole blood from control subjects (n = 21) and patients with ET (n = 19), PV (n = 41), and PMF (n = 9) using DNA microarrays. Applying an unsupervised method, principal component analysis, to search for patterns in the data, we demonstrated a separation of the four groups with biological relevant overlaps between the different entities. Moreover, the analysis separates Janus activating kinase 2-negative ET patients from Janus activating kinase 2-positive ET patients. Functional annotation analysis demonstrates that clusters of gene ontology terms related to inflammation, immune system, apoptosis, RNA metabolism, and secretory system were the most significantly deregulated terms in the three different disease groups. Our results yield further support for the hypothesis of a biological continuum originating from ET over PV to PMF. Functional analysis suggests an important implication of these gene ontology clusters in the pathogenesis of these neoplasms and in disease evolution from ET over PV to PMF.

Evaluation of Nine Somatic Variant Callers for Detection of Somatic Mutations in Exome and Targeted Deep Sequencing Data.

Next generation sequencing is extensively applied to catalogue somatic mutations in cancer, in research settings and increasingly in clinical settings for molecular diagnostics, guiding therapy decisions. Somatic variant callers perform paired comparisons of sequencing data from cancer tissue and matched normal tissue in order to detect somatic mutations. The advent of many new somatic variant callers creates a need for comparison and validation of the tools, as no de facto standard for detection of somatic mutations exists and only limited comparisons have been reported. We have performed a comprehensive evaluation using exome sequencing and targeted deep sequencing data of paired tumor-normal samples from five breast cancer patients to evaluate the performance of nine publicly available somatic variant callers: EBCall, Mutect, Seurat, Shimmer, Indelocator, Somatic Sniper, Strelka, VarScan 2 and Virmid for the detection of single nucleotide mutations and small deletions and insertions. We report a large variation in the number of calls from the nine somatic variant callers on the same sequencing data and highly variable agreement. Sequencing depth had markedly diverse impact on individual callers, as for some callers, increased sequencing depth highly improved sensitivity. For SNV calling, we report EBCall, Mutect, Virmid and Strelka to be the most reliable somatic variant callers for both exome sequencing and targeted deep sequencing. For indel calling, EBCall is superior due to high sensitivity and robustness to changes in sequencing depths.

Molecular profiling of peripheral blood cells from patients with polycythemia vera and related neoplasms: Identification of deregulated genes of significance for inflammation and immune surveillance

Essential thrombocythemia (ET), polycythemia vera (PV) and primary myelofibrosis (PMF) are hematopoietic stem cell neoplasms that may be associated with autoimmune or chronic inflammatory disorders. Earlier gene expression profiling studies have demonstrated aberrant expression of genes involved in inflammatory responses, mainly being performed on granulocytes or CD34+ cells. Using gene expression profiling of whole blood from patients with ET (n=16), PV (n=36), and PMF (n=9), several genes involved in inflammation and immune regulation were found to be significantly deregulated. Our findings may reflect chronic inflammation to be of pathogenetic importance for the progression of these neoplasms toward the myelofibrotic end-stage and may also account for the increased frequency of second cancer in these diseases.

Nordic collaborative study of the BARD1 Cys557Ser allele in 3956 patients with cancer: enrichment in familial BRCA1/BRCA2 mutation-negative breast cancer but not in other malignancies

Background: BARD1 was originally identified as a BRCA1-interacting protein but has also been described in tumour-suppressive functions independent of BRCA1. Several studies have indicated that the BARD1 gene is a potential target for germline changes predisposing to breast and ovarian cancer. The C-terminal Cys557Ser change has previously been uncovered to associate with an increased risk of breast cancer and was recently shown to result in defective apoptotic activities. Aim and methods: Conformation-sensitive gel electrophoresis, minisequencing, TaqMan assays, denaturing high-performance liquid chromatography analysis and DNA sequencing were used to investigate the prevalence of the Cys557Ser allele in a large Nordic case–control study cohort consisting of 2906 patients with breast or ovarian cancer, 734 with prostate cancer, 188 with colorectal cancer, 128 men with breast cancer, and 3591 controls from Finland, Iceland, Denmark, Sweden and Norway. Results: The frequency of the BARD1 Cys557Ser variant seemed to increase among patients from families with breast or ovarian cancer lacking BRCA1 or BRCA2 mutations: a significant difference was obtained compared with controls (6.8% v 2.7%; p<0.001; odds ratio (OR) 2.6; 95% confidence interval (CI) 1.7 to 4.0) and with patients from BRCA1/BRCA2 mutation-positive families (6.8% v 2.2%; p = 0.01; OR 3.2; 95% CI 1.2 to 8.3). In contrast, no major association with male breast, ovarian, colorectal or prostate cancer was observed. Additionally, a novel BARD1 allele resulting in Ser558Pro was identified in familial breast cancer cases. Conclusion: These results provide further evidence that BARD1 Cys557Ser confers a slightly increased risk of breast cancer in women.

Increased gene expression of histone deacetylases in patients with Philadelphia-negative chronic myeloproliferative neoplasms

Abstract Myeloproliferation, myeloaccumulation (decreased apoptosis), inflammation, bone marrow fibrosis and angiogenesis are cardinal features of the Philadelphia-negative chronic myeloproliferative neoplasms: essential thrombocythemia (ET), polycythemia vera (PV) and primary myelofibrosis (PMF). Histone deacetylases (HDACs) have a critical role in modulating gene expression and, accordingly, in the control of cell pathobiology and cancer development. HDAC inhibition has been shown to inhibit tumor growth (impaired myeloproliferation), to modulate the balance between pro- and antiapoptotic proteins in favor of apoptosis (enhanced apoptosis) and also to inhibit angiogenesis. Recently, enhanced HDAC enzyme activity has been found in CD34+cells from patients with PMF, enzyme activity levels highly exceeding those recorded in other chronic myeloproliferative neoplasms (CMPNs). The raised levels correlated to the degree of splenomegaly, suggesting that HDAC might be recruited as ET or PV progresses into myelofibrosis or PMF progresses into a more advanced stage. Accordingly, HDAC inhibition is an obvious novel therapeutic approach in these neoplasms. Using global gene expression profiling of whole blood from patients with CMPNs, we have found a pronounced deregulation of HDAC genes, involving significant up-regulation of the HDAC genes 9 and 11, with the highest expression levels being found in patients with ET (HDAC9 and 11), PMF (HDAC9) and CMPNs (both HDAC9 and HDAC11). Furthermore, we have identified that the HDAC6 gene is progressively expressed in patients with ET, PV and PMF, reflecting a steady accumulation of abnormally expressed HDAC6 during disease evolution. Our results lend further support to HDACs as important epigenetic targets in the future treatment of patients with CMPNs. Since the highest expression levels of HDAC genes were recorded in ET, in PMF and in the entire CMPN group, their down-regulation by HDAC inhibitors might be associated with decreased disease activity, including reduction of splenomegaly.