Aggregation tests identify new gene associations with breast cancer in populations with diverse ancestry.

BACKGROUND: Low-frequency variants play an important role in breast cancer (BC) susceptibility. Gene-based methods can increase power by combining multiple variants in the same gene and help identify target genes. METHODS: We evaluated the potential of gene-based aggregation in the Breast Cancer Association Consortium cohorts including 83,471 cases and 59,199 controls. Low-frequency variants were aggregated for individual genes' coding and regulatory regions. Association results in European ancestry samples were compared to single-marker association results in the same cohort. Gene-based associations were also combined in meta-analysis across individuals with European, Asian, African, and Latin American and Hispanic ancestry. RESULTS: In European ancestry samples, 14 genes were significantly associated (q < 0.05) with BC. Of those, two genes, FMNL3 (P = 6.11 × 10-6) and AC058822.1 (P = 1.47 × 10-4), represent new associations. High FMNL3 expression has previously been linked to poor prognosis in several other cancers. Meta-analysis of samples with diverse ancestry discovered further associations including established candidate genes ESR1 and CBLB. Furthermore, literature review and database query found further support for a biologically plausible link with cancer for genes CBLB, FMNL3, FGFR2, LSP1, MAP3K1, and SRGAP2C. CONCLUSIONS: Using extended gene-based aggregation tests including coding and regulatory variation, we report identification of plausible target genes for previously identified single-marker associations with BC as well as the discovery of novel genes implicated in BC development. Including multi ancestral cohorts in this study enabled the identification of otherwise missed disease associations as ESR1 (P = 1.31 × 10-5), demonstrating the importance of diversifying study cohorts.

Genome- and transcriptome-wide association studies of 386,000 Asian and European-ancestry women provide new insights into breast cancer genetics.

By combining data from 160,500 individuals with breast cancer and 226,196 controls of Asian and European ancestry, we conducted genome- and transcriptome-wide association studies of breast cancer. We identified 222 genetic risk loci and 137 genes that were associated with breast cancer risk at a p < 5.0 × 10-8 and a Bonferroni-corrected p < 4.6 × 10-6, respectively. Of them, 32 loci and 15 genes showed a significantly different association between ER-positive and ER-negative breast cancer after Bonferroni correction. Significant ancestral differences in risk variant allele frequencies and their association strengths with breast cancer risk were identified. Of the significant associations identified in this study, 17 loci and 14 genes are located 1Mb away from any of the previously reported breast cancer risk variants. Pathways analyses including 221 putative risk genes identified multiple signaling pathways that may play a significant role in the development of breast cancer. Our study provides a comprehensive understanding of and new biological insights into the genetics of this common malignancy.

Relevance of the MHC region for breast cancer susceptibility in Asians.

BACKGROUND: Human leukocyte antigen (HLA) genes play critical roles in immune surveillance, an important defence against tumors. Imputing HLA genotypes from existing single-nucleotide polymorphism datasets is low-cost and efficient. We investigate the relevance of the major histocompatibility complex region in breast cancer susceptibility, using imputed class I and II HLA alleles, in 25,484 women of Asian ancestry. METHODS: A total of 12,901 breast cancer cases and 12,583 controls from 12 case-control studies were included in our pooled analysis. HLA imputation was performed using SNP2HLA on 10,886 quality-controlled variants within the 15-55 Mb region on chromosome 6. HLA alleles (n = 175) with info scores greater than 0.8 and frequencies greater than 0.01 were included (resolution at two-digit level: 71; four-digit level: 104). We studied the associations between HLA alleles and breast cancer risk using logistic regression, adjusting for population structure and age. Associations between HLA alleles and the risk of subtypes of breast cancer (ER-positive, ER-negative, HER2-positive, HER2-negative, early-stage, and late-stage) were examined. RESULTS: We did not observe associations between any HLA allele and breast cancer risk at P < 5e-8; the smallest p value was observed for HLA-C*12:03 (OR = 1.29, P = 1.08e-3). Ninety-five percent of the effect sizes (OR) observed were between 0.90 and 1.23. Similar results were observed when different subtypes of breast cancer were studied (95% of ORs were between 0.85 and 1.18). CONCLUSIONS: No imputed HLA allele was associated with breast cancer risk in our large Asian study. Direct measurement of HLA gene expressions may be required to further explore the associations between HLA genes and breast cancer risk.

Polygenic risk scores for prediction of breast cancer risk in Asian populations.

PURPOSE: Non-European populations are under-represented in genetics studies, hindering clinical implementation of breast cancer polygenic risk scores (PRSs). We aimed to develop PRSs using the largest available studies of Asian ancestry and to assess the transferability of PRS across ethnic subgroups. METHODS: The development data set comprised 138,309 women from 17 case-control studies. PRSs were generated using a clumping and thresholding method, lasso penalized regression, an Empirical Bayes approach, a Bayesian polygenic prediction approach, or linear combinations of multiple PRSs. These PRSs were evaluated in 89,898 women from 3 prospective studies (1592 incident cases). RESULTS: The best performing PRS (genome-wide set of single-nucleotide variations [formerly single-nucleotide polymorphism]) had a hazard ratio per unit SD of 1.62 (95% CI = 1.46-1.80) and an area under the receiver operating curve of 0.635 (95% CI = 0.622-0.649). Combined Asian and European PRSs (333 single-nucleotide variations) had a hazard ratio per SD of 1.53 (95% CI = 1.37-1.71) and an area under the receiver operating curve of 0.621 (95% CI = 0.608-0.635). The distribution of the latter PRS was different across ethnic subgroups, confirming the importance of population-specific calibration for valid estimation of breast cancer risk. CONCLUSION: PRSs developed in this study, from association data from multiple ancestries, can enhance risk stratification for women of Asian ancestry.

Body mass index and type 2 diabetes and breast cancer survival: a Mendelian randomization study.

The causal relationship between body mass index (BMI) and type 2 diabetes (T2D) and breast cancer prognosis is still ambiguous. The aim of this study was to investigate the prognostic effect of BMI and T2D on breast cancer disease-free survival (DFS) among Asian individuals. In this two-sample Mendelian randomization (MR) study, the instrumental variables (IVs) were identified using a genome-wide association study (GWAS) among 24,000 participants in the Taiwan Biobank. Importantly, the validity of these IVs was confirmed with a previous large-scale GWAS (Biobank Japan Project, BBJ). In this study, we found that a genetic predisposition toward higher BMI (as indicated by BMI IVs, F = 86.88) was associated with poor breast cancer DFS (hazard ratio [HR] = 6.11; P < 0.001). Furthermore, higher level of genetically predicted T2D (as indicated by T2D IVs) was associated with an increased risk of recurrence of and mortality from breast cancer (HR = 1.43; P < 0.001). Sensitivity analyses, including the weighted-median approach, MR-Egger regression, Radial regression and Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO) supported the consistency of our findings. Finally, the causal relationship between BMI and poor breast cancer prognosis was confirmed in a prospective cohort study. Our MR analyses demonstrated the causal relationship between the genetic prediction of elevated BMI and a greater risk of T2D with poor breast cancer prognosis. BMI and T2D have important clinical implications and may be used as prognostic indicators of breast cancer.

Functional annotation of the 2q35 breast cancer risk locus implicates a structural variant in influencing activity of a long-range enhancer element.

A combination of genetic and functional approaches has identified three independent breast cancer risk loci at 2q35. A recent fine-scale mapping analysis to refine these associations resulted in 1 (signal 1), 5 (signal 2), and 42 (signal 3) credible causal variants at these loci. We used publicly available in silico DNase I and ChIP-seq data with in vitro reporter gene and CRISPR assays to annotate signals 2 and 3. We identified putative regulatory elements that enhanced cell-type-specific transcription from the IGFBP5 promoter at both signals (30- to 40-fold increased expression by the putative regulatory element at signal 2, 2- to 3-fold by the putative regulatory element at signal 3). We further identified one of the five credible causal variants at signal 2, a 1.4 kb deletion (esv3594306), as the likely causal variant; the deletion allele of this variant was associated with an average additional increase in IGFBP5 expression of 1.3-fold (MCF-7) and 2.2-fold (T-47D). We propose a model in which the deletion allele of esv3594306 juxtaposes two transcription factor binding regions (annotated by estrogen receptor alpha ChIP-seq peaks) to generate a single extended regulatory element. This regulatory element increases cell-type-specific expression of the tumor suppressor gene IGFBP5 and, thereby, reduces risk of estrogen receptor-positive breast cancer (odds ratio = 0.77, 95% CI 0.74-0.81, p = 3.1 × 10-31).

Breast Cancer Polygenic Risk Score and Contralateral Breast Cancer Risk.

Previous research has shown that polygenic risk scores (PRSs) can be used to stratify women according to their risk of developing primary invasive breast cancer. This study aimed to evaluate the association between a recently validated PRS of 313 germline variants (PRS313) and contralateral breast cancer (CBC) risk. We included 56,068 women of European ancestry diagnosed with first invasive breast cancer from 1990 onward with follow-up from the Breast Cancer Association Consortium. Metachronous CBC risk (N = 1,027) according to the distribution of PRS313 was quantified using Cox regression analyses. We assessed PRS313 interaction with age at first diagnosis, family history, morphology, ER status, PR status, and HER2 status, and (neo)adjuvant therapy. In studies of Asian women, with limited follow-up, CBC risk associated with PRS313 was assessed using logistic regression for 340 women with CBC compared with 12,133 women with unilateral breast cancer. Higher PRS313 was associated with increased CBC risk: hazard ratio per standard deviation (SD) = 1.25 (95%CI = 1.18-1.33) for Europeans, and an OR per SD = 1.15 (95%CI = 1.02-1.29) for Asians. The absolute lifetime risks of CBC, accounting for death as competing risk, were 12.4% for European women at the 10th percentile and 20.5% at the 90th percentile of PRS313. We found no evidence of confounding by or interaction with individual characteristics, characteristics of the primary tumor, or treatment. The C-index for the PRS313 alone was 0.563 (95%CI = 0.547-0.586). In conclusion, PRS313 is an independent factor associated with CBC risk and can be incorporated into CBC risk prediction models to help improve stratification and optimize surveillance and treatment strategies.

European polygenic risk score for prediction of breast cancer shows similar performance in Asian women.

Polygenic risk scores (PRS) have been shown to predict breast cancer risk in European women, but their utility in Asian women is unclear. Here we evaluate the best performing PRSs for European-ancestry women using data from 17,262 breast cancer cases and 17,695 controls of Asian ancestry from 13 case-control studies, and 10,255 Chinese women from a prospective cohort (413 incident breast cancers). Compared to women in the middle quintile of the risk distribution, women in the highest 1% of PRS distribution have a ~2.7-fold risk and women in the lowest 1% of PRS distribution has ~0.4-fold risk of developing breast cancer. There is no evidence of heterogeneity in PRS performance in Chinese, Malay and Indian women. A PRS developed for European-ancestry women is also predictive of breast cancer risk in Asian women and can help in developing risk-stratified screening programmes in Asia.

Identification of novel breast cancer susceptibility loci in meta-analyses conducted among Asian and European descendants.

Known risk variants explain only a small proportion of breast cancer heritability, particularly in Asian women. To search for additional genetic susceptibility loci for breast cancer, here we perform a meta-analysis of data from genome-wide association studies (GWAS) conducted in Asians (24,206 cases and 24,775 controls) and European descendants (122,977 cases and 105,974 controls). We identified 31 potential novel loci with the lead variant showing an association with breast cancer risk at P < 5 × 10-8. The associations for 10 of these loci were replicated in an independent sample of 16,787 cases and 16,680 controls of Asian women (P < 0.05). In addition, we replicated the associations for 78 of the 166 known risk variants at P < 0.05 in Asians. These findings improve our understanding of breast cancer genetics and etiology and extend previous findings from studies of European descendants to Asian women.

Germline Mutation in 1338 BRCA-Negative Chinese Hereditary Breast and/or Ovarian Cancer Patients: Clinical Testing with a Multigene Test Panel.

Differences in the mutation spectrum across ethnicities suggest the importance of identifying genes in addition to common high penetrant genes to estimate the associated breast cancer risk in China. A total of 1338 high-risk breast cancer patients who tested negative for germline BRCA1, BRCA2, TP53, and PTEN mutations between 2007 and 2017 were selected from the Hong Kong Hereditary Breast Cancer Family Registry. Patient samples were subjected to next-generation DNA sequencing using a multigene panel (Color Genomics). All detected pathogenic variants were validated by bidirectional DNA sequencing. The sequencing data were coanalyzed by a bioinformatics pipeline developed in-house. Sixty-one pathogenic variants (4.6%) were identified in this cohort in 11 cancer predisposition genes. Most carriers (77.1%) had early onset of breast cancer (age <45 years), 32.8% had family members with breast cancer, and 11.5% had triple-negative breast cancer. The most common mutated genes were PALB2 (1.4%), RAD51D (0.8%), and ATM (0.8%). A total of 612 variants of unknown significance were identified in 494 patients, and 87.4% of the variants of unknown significance were missense mutations. Pathogenic variants in cancer predisposition genes beyond BRCA1, BRCA2, TP53, and PTEN were detected in an additional 4.6% of patients using the multigene panel. PALB2 (1.4%) and RAD51D (0.8%) were the most commonly mutated genes in patients who tested mutation negative by a four-gene panel.

Association of Genomic Domains in BRCA1 and BRCA2 with Prostate Cancer Risk and Aggressiveness.

Pathogenic sequence variants (PSV) in BRCA1 or BRCA2 (BRCA1/2) are associated with increased risk and severity of prostate cancer. We evaluated whether PSVs in BRCA1/2 were associated with risk of overall prostate cancer or high grade (Gleason 8+) prostate cancer using an international sample of 65 BRCA1 and 171 BRCA2 male PSV carriers with prostate cancer, and 3,388 BRCA1 and 2,880 BRCA2 male PSV carriers without prostate cancer. PSVs in the 3' region of BRCA2 (c.7914+) were significantly associated with elevated risk of prostate cancer compared with reference bin c.1001-c.7913 [HR = 1.78; 95% confidence interval (CI), 1.25-2.52; P = 0.001], as well as elevated risk of Gleason 8+ prostate cancer (HR = 3.11; 95% CI, 1.63-5.95; P = 0.001). c.756-c.1000 was also associated with elevated prostate cancer risk (HR = 2.83; 95% CI, 1.71-4.68; P = 0.00004) and elevated risk of Gleason 8+ prostate cancer (HR = 4.95; 95% CI, 2.12-11.54; P = 0.0002). No genotype-phenotype associations were detected for PSVs in BRCA1. These results demonstrate that specific BRCA2 PSVs may be associated with elevated risk of developing aggressive prostate cancer. SIGNIFICANCE: Aggressive prostate cancer risk in BRCA2 mutation carriers may vary according to the specific BRCA2 mutation inherited by the at-risk individual.

Two truncating variants in FANCC and breast cancer risk.

Fanconi anemia (FA) is a genetically heterogeneous disorder with 22 disease-causing genes reported to date. In some FA genes, monoallelic mutations have been found to be associated with breast cancer risk, while the risk associations of others remain unknown. The gene for FA type C, FANCC, has been proposed as a breast cancer susceptibility gene based on epidemiological and sequencing studies. We used the Oncoarray project to genotype two truncating FANCC variants (p.R185X and p.R548X) in 64,760 breast cancer cases and 49,793 controls of European descent. FANCC mutations were observed in 25 cases (14 with p.R185X, 11 with p.R548X) and 26 controls (18 with p.R185X, 8 with p.R548X). There was no evidence of an association with the risk of breast cancer, neither overall (odds ratio 0.77, 95%CI 0.44-1.33, p = 0.4) nor by histology, hormone receptor status, age or family history. We conclude that the breast cancer risk association of these two FANCC variants, if any, is much smaller than for BRCA1, BRCA2 or PALB2 mutations. If this applies to all truncating variants in FANCC it would suggest there are differences between FA genes in their roles on breast cancer risk and demonstrates the merit of large consortia for clarifying risk associations of rare variants.

DNA methylation at an enhancer of the three prime repair exonuclease 2 gene (TREX2) is linked to gene expression and survival in laryngeal cancer.

BACKGROUND: Genetic aberrations in DNA repair genes are linked to cancer, but less is reported about epigenetic regulation of DNA repair and functional consequences. We investigated the intragenic methylation loss at the three prime repair exonuclease 2 (TREX2) locus in laryngeal (n = 256) and colorectal cancer cases (n = 95) and in pan-cancer data from The Cancer Genome Atlas (TCGA). RESULTS: Significant methylation loss at an intragenic site of TREX2 was a frequent trait in both patient cohorts (p = 0.016 and < 0.001, respectively) and in 15 out of 22 TCGA studies. Methylation loss correlated with immunohistochemically staining for TREX2 (p < 0.0001) in laryngeal tumors and improved overall survival of laryngeal cancer patients (p = 0.045). Chromatin immunoprecipitation, demethylation experiments, and reporter gene assays revealed that the region of methylation loss can function as a CCAAT/enhancer binding protein alpha (CEBPA)-responsive enhancer element regulating TREX2 expression. CONCLUSIONS: The data highlight a regulatory role of TREX2 DNA methylation for gene expression which might affect incidence and survival of laryngeal cancer. Altered TREX2 protein levels in tumors may affect drug-induced DNA damage repair and provide new tailored therapies.

Distinct mutation spectrum, clinical outcome and therapeutic responses of typical complex/monosomy karyotype acute myeloid leukemia carrying TP53 mutations.

The present study aimed to define a subtype of complex/monosomal karyotype (CK/MK) acute myeloid leukemia (AML) by its distinct clinical features, p53 signaling and responses to p53 targeting agents. Ninety-eight young adults (range: 21-60 years; median: 49 years) with CK/MK AML were studied. They received standard induction, consolidation and allogeneic hematopoietic stem cell transplantation from siblings or matched unrelated donors if available. Chromosomal abnormalities most commonly affected chromosome 5 (30%), 7 (22%) and 17 (21%). Next generation sequencing of a 54-myeloid gene panel were available in 76 patients. Tumor protein 53 (TP53) mutations were most common (49%) and associated with the presence of -5/5q- (P < .001) and -17/17p- (P < .001), but not -7/7q- (P = .370). This "typical" CK/MK AML subtype was associated with significantly lower presenting white cell counts, higher number of karyotypic abnormalities, and inferior leukemia-free and overall survivals, compared with CK/MK AML without the typical features. Blood or bone marrow samples from typical CK/MK AML patients showed defective p53 signaling upon induction by etoposide. In vitro drug sensitivity analysis showed that they were sensitive to APR-246 that targeted mutant p53, but resistant to MDM2 antagonist MI-77301. Novel therapeutic strategies targeting TP53 mutations in CK/MK AML should be developed and tested in clinical trials.

Breast and ovarian cancer penetrance of BRCA1/2 mutations among Hong Kong women.

Germline mutations in BRCA1 and BRCA2 (BRCA1/2) are associated with increased risk of breast and ovarian cancer. The penetrance of breast and ovarian cancer in BRCA1/2 mutation carriers has been well characterized in Caucasian but not in Asian. Two studies have investigated the breast cancer risk in Asian women with BRCA1/2 mutations, and no published estimates are available for ovarian cancer. Therefore, we estimated the age-specific cumulative risk of BRCA1/2-associated breast and ovarian cancer in Chinese women. From Jan 2007 to Nov 2015, the Hong Kong Hereditary Breast Cancer Family Registry identified 1635 families with hereditary breast-ovarian cancer. Among probands in these families, 66 had BRCA1 mutations, 84 had BRCA2 mutations, and 1,485 tested negative for BRCA1/2 mutations. Using the female first-degree relatives of these probands, we estimated the risk of breast and ovarian cancer using a modified marginal likelihood approach. Estimates of breast cancer penetrance by age 70 were 53.7% (95% CI 34.5-71.6%) for BRCA1 mutation carriers and 48.3% (95% CI 31.8-68.5%) for BRCA2. The estimated risk of ovarian cancer by age 70 was 21.5% and 7.3% for Chinese women carrying BRCA1 or BRCA2 mutation respectively. A meta-analysis of available studies in Asian women revealed pooled estimates of breast cancer risk by age 70 of 44.8% (95% CI 33-57.2%) and 40.7% (95% CI 31.3-50.9%) for BRCA1 and BRCA2 mutation carriers respectively. These data suggest that BRCA1/2-associated breast cancer risk for Chinese women is similar to that for Caucasian women, although BRCA1/2-associated ovarian cancer risks are lower for Chinese women.

Rapid detection of BRCA1/2 recurrent mutations in Chinese breast and ovarian cancer patients with multiplex SNaPshot genotyping panels.

BRCA1/2 mutations are significant risk factors for hereditary breast and ovarian cancer (HBOC), its mutation frequency in HBOC of Chinese ethnicity is around 9%, in which nearly half are recurrent mutations. In Hong Kong and China, genetic testing and counseling are not as common as in the West. To reduce the barrier of testing, a multiplex SNaPshot genotyping panel that targeted 25 Chinese BRCA1/2 mutation hotspots was developed, and its feasibility was evaluated in a local cohort of 441 breast and 155 ovarian cancer patients. For those who tested negative, they were then subjected to full-gene testing with next-generation sequencing (NGS). BRCA mutation prevalence in this cohort was 8.05% and the yield of the recurrent panel was 3.52%, identifying over 40% of the mutation carriers. Moreover, from 79 Chinese breast cancer cases recruited overseas, 2 recurrent mutations and one novel BRCA2 mutation were detected by the panel and NGS respectively. The developed genotyping panel showed to be an easy-to-perform and more affordable testing tool that can provide important contributions to improve the healthcare of Chinese women with cancer as well as family members that harbor high risk mutations for HBOC.

Detection of Germline Mutation in Hereditary Breast and/or Ovarian Cancers by Next-Generation Sequencing on a Four-Gene Panel.

Mutation in BRCA1/BRCA2 genes accounts for 20% of familial breast cancers, 5% to 10% of which may be due to other less penetrant genes which are still incompletely studied. Herein, a four-gene panel was used to examine the prevalence of BRCA1, BRCA2, TP53, and PTEN in hereditary breast and ovarian cancers in Southern Chinese population. In this cohort, 948 high-risk breast and/or ovarian patients were recruited for genetic screening by next-generation sequencing (NGS). The performance of our NGS pipeline was evaluated with 80 Sanger-validated known mutations and eight negative cases. With appropriate bioinformatics analysis pipeline, the detection sensitivity of NGS is comparable with Sanger sequencing. The prevalence of BRCA1/BRCA2 germline mutations was 9.4% in our Chinese cohort, of which 48.8% of the mutations arose from hotspot mutations. With the use of a tailor-made algorithm, HomopolymerQZ, more mutations were detected compared with single mutation detection algorithm. The frequencies of PTEN and TP53 were 0.21% and 0.53%, respectively, in the Southern Chinese patients with breast and/or ovarian cancers. High-throughput NGS approach allows the incorporation of control cohort that provides an ethnicity-specific data for polymorphic variants. Our data suggest that hotspot mutations screening such as SNaPshot could be an effective preliminary screening alternative adopted in a standard clinical laboratory without NGS setup.

Germline RECQL mutations in high risk Chinese breast cancer patients.

Recently, RECQL was reported as a new breast cancer susceptibility gene. RECQL belongs to the RECQ DNA helicase family which unwinds double strand DNA and involved in the DNA replication stress response, telomere maintenance and DNA repair. RECQL deficient mice cells are prone to spontaneous chromosomal instability and aneuploidy, suggesting a tumor-suppressive role of RECQL in cancer. In this study, RECQL gene mutation screening was performed on 1110 breast cancer patients who were negative for BRCA1, BRCA2, TP53 and PTEN gene mutations and recruited from March 2007 to June 2015 in the Hong Kong Hereditary and High Risk Breast Cancer Program. Four different RECQL pathogenic mutations were identified in six of the 1110 (0.54 %) tested breast cancer patients. The identified mutations include one frame-shift deletion (c.974_977delAAGA), two splicing site mutations (c.394+1G>A, c.867+1G>T) and one nonsense mutation (c.796C>T, p.Gln266Ter). Two of the mutations (c.867+1G>T and p.Gln266Ter) were seen in more than one patients. This study provides the basis for existing of pathogenic RECQL mutations in Southern Chinese breast cancer patients. The significance of rare variants in RECQL gene in the estimation of breast cancer risk warranted further investigation in larger cohort of patients and in other ethnic groups.

Altered regulation of DNA ligase IV activity by aberrant promoter DNA methylation and gene amplification in colorectal cancer.

Colorectal cancer (CRC) presents as a very heterogeneous disease which cannot sufficiently be characterized with the currently known genetic and epigenetic markers. To identify new markers for CRC we scrutinized the methylation status of 231 DNA repair-related genes by methyl-CpG immunoprecipitation followed by global methylation profiling on a CpG island microarray, as altered expression of these genes could drive genomic and chromosomal instability observed in these tumors. We show for the first time hypermethylation of MMP9, DNMT3A and LIG4 in CRC which was confirmed in two CRC patient groups with different ethnicity. DNA ligase IV (LIG4) showed strong differential promoter methylation (up to 60%) which coincided with downregulation of mRNA in 51% of cases. This functional association of LIG4 methylation and gene expression was supported by LIG4 re-expression in 5-aza-2'-deoxycytidine-treated colon cancer cell lines, and reduced ligase IV amounts and end-joining activity in extracts of tumors with hypermethylation. Methylation of LIG4 was not associated with other genetic and epigenetic markers of CRC in our study. As LIG4 is located on chromosome 13 which is frequently amplified in CRC, two loci were tested for gene amplification in a subset of 47 cases. Comparison of amplification, methylation and expression data revealed that, in 30% of samples, the LIG4 gene was amplified and methylated, but expression was not changed. In conclusion, hypermethylation of the LIG4 promoter is a new mechanism to control ligase IV expression. It may represent a new epigenetic marker for CRC independent of known markers.