- Research article
- Open Access
Patient survival and tumor characteristics associated with CHEK2:p.I157T – findings from the Breast Cancer Association Consortium
Breast Cancer Research volume 18, Article number: 98 (2016)
P.I157T is a CHEK2 missense mutation associated with a modest increase in breast cancer risk. Previously, another CHEK2 mutation, the protein truncating c.1100delC has been associated with poor prognosis of breast cancer patients. Here, we have investigated patient survival and characteristics of breast tumors of germ line p.I157T carriers.
We included in the analyses 26,801 European female breast cancer patients from 15 studies participating in the Breast Cancer Association Consortium. We analyzed the association between p.I157T and the clinico-pathological breast cancer characteristics by comparing the p.I157T carrier tumors to non-carrier and c.1100delC carrier tumors. Similarly, we investigated the p.I157T associated risk of early death, breast cancer-associated death, distant metastasis, locoregional relapse and second breast cancer using Cox proportional hazards models.
Additionally, we explored the p.I157T-associated genomic gene expression profile using data from breast tumors of 183 Finnish female breast cancer patients (ten p.I157T carriers) (GEO: GSE24450). Differential gene expression analysis was performed using a moderated t test. Functional enrichment was investigated using the DAVID functional annotation tool and gene set enrichment analysis (GSEA). The tumors were classified into molecular subtypes according to the St Gallen 2013 criteria and the PAM50 gene expression signature.
P.I157T was not associated with increased risk of early death, breast cancer-associated death or distant metastasis relapse, and there was a significant difference in prognosis associated with the two CHEK2 mutations, p.I157T and c.1100delC. Furthermore, p.I157T was associated with lobular histological type and clinico-pathological markers of good prognosis, such as ER and PR expression, low TP53 expression and low grade. Gene expression analysis suggested luminal A to be the most common subtype for p.I157T carriers and CDH1 (cadherin 1) target genes to be significantly enriched among genes, whose expression differed between p.I157T and non-carrier tumors.
Our analyses suggest that there are fundamental differences in breast tumors of CHEK2:p.I157T and c.1100delC carriers. The poor prognosis associated with c.1100delC cannot be generalized to other CHEK2 mutations.
Checkpoint kinase 2 (CHEK2) is a moderate penetrance breast cancer risk gene. The two most frequent CHEK2 mutations in European populations are p.I157T and c.1100delC. Truncating CHEK2 founder mutations (c.1100delC, IVS2 + 1G > A, del5395) confer a higher than twofold increase in the risk of breast cancer [1–3], whereas p.I157T (c.470 T > C, rs17879961), a CHEK2 missense mutation is associated with a milder, 1.4-fold elevation in the risk . The c.1100delC carrier frequency is highest in the Netherlands and in Finland (over 1 %), the other two truncating founder mutations are found mainly in Poland , and p.I157T is most frequent in Finland and in Poland (around 5 %) . Additionally, dozens of rare CHEK2 missense mutations have been found in breast cancer patients, but their contribution to disease risk is minor on a population level and causative role in disease development probably varies greatly [6–8].
The consequences of c.1100delC and p.I157T differ on a molecular level, but both have been shown to severely interfere with the CHEK2 protein activity. C.1100delC is a loss-of-function mutation that induces a premature termination codon in the kinase domain in exon 10 (ter381) leading to a nonsense-mediated mRNA decay, which reduces both mutated and overall CHEK2 mRNA level [9, 10]. C.1100delC truncates CHEK2 protein’s C-terminal kinase domain. The truncated protein is unstable and practically undetectable in mutation carrier cells . Isoleucine 157 (p.I157T) is required for several van der Waals interactions at the interface of forkhead-associated (FHA) and kinase domains of dimerizing CHEK2 peptide chains. Its replacement to threonine (p.I157T) has been shown to interfere with these interactions and to severely impede the CHEK2 homodimerization required for its activation . Furthermore, ectopic expression of human CHEK2:p.I157T failed a rad53/sml complementation assay in yeast suggesting an impaired protein function . Thus, p.I157T possibly disturbs CHEK2 function by competing with the wild-type protein in dimer formation in heterozygous cells in a dominant negative manner .
Since both p.I157T and c.1100delC cause increased risk of breast cancer and compromise the activity of the CHEK2 protein, the question remains whether their effects on patient prognosis would be proportional to their risk effects and how similar the breast cancer phenotypes associated with the mutations would be. C.1100delC is associated with bilateral disease and estrogen receptor (ER)-positive tumors [12–14]. However, although tumors from p.I157T carriers are also predominantly ER-positive , tumors from p.I157T and c.1100delC carriers are associated with phenotypically different types of breast cancer. The lobular histological type is overrepresented among p.I157T mutation carrier tumors , whereas the c.1100delC carrier tumors are typically ductal [13, 14].
We have previously reported CHEK2:c.1100delC heterozygosity to be associated with reduced overall and disease-free survival as well as with increased risk of breast cancer-specific death in a Breast Cancer Association Consortium (BCAC) data set combining mutation carriers from multiple European populations . Here, we report a study investigating thoroughly the prognostic associations of CHEK2:p.I157T as well as pathologic characteristics and genomic gene expression profiles of breast tumors from carriers of germ line p.I157T.
Study subjects for survival and pathology analyses
We included in the analyses female invasive breast cancer patients of European ancestry with a first invasive primary breast cancer enrolled in 15 studies participating in the Breast Cancer Association Consortium (BCAC) (Additional file 1: Table S1). In order to be able to stratify the analyses by study, only BCAC studies providing genotype and survival data of about ten CHEK2:p.I157T carriers were included in the analyses (Additional file 1: Table S2). Altogether, the data set consisted of 26,801 study subjects, of which 590 carried germ line p.I157T and 271 carried c.1100delC mutations (Table 1). Individuals carrying both mutations were excluded from the analyses (n = 4).
CHEK2:p.I157T was first genotyped by independent studies using various methods including MassARRAY iPLEX Gold (Sequenom, San Diego, CA, USA), TaqMan (Applied Biosystems, Life Technologies, Carlsbad, CA, USA) and Fluidigm (Fluidigm, San Francisco, CA, USA) as listed in Additional file 1: Table S1. Quality control was implemented as follows: each study performed duplicate measurements of at least two samples from each sample plate and genotyped 93 CEPH control DNAs (HAPMAPPT01, Coriell Institute for Medical Research, Cambden, NJ, USA). If a study reported more than two discordant genotyping results of the CEPH DNAs, all genotype data from that study was excluded. Later, p.I157T was genotyped centrally using a custom Illumina iSelect genotyping array for the Collaborative Oncological Gene-environment Study (COGS) . Discordant genotyping results were clarified with Sanger sequencing. CHEK2:c.1100delC was genotyped by independent studies using mainly TaqMan (Additional file 1: Table S1), as described earlier .
Pathology data was collected from hospital records or from scientific projects within the individual studies, as described previously . Additionally, the TP53 protein expression was measured by individual studies using immunohistochemical staining as described in Additional file 1: Table S3. The pathology data availability and mutation carrier frequencies varied between independent BCAC studies and therefore all analyses were stratified by study. Pathology analyses were performed using R environment for statistical computing version 3.0.2  including packages vcdExtra  and meta . Comparisons were made between CHEK2 mutations carriers (heterozygous or homozygous) and non-carriers, for both p.I157T and c.1100delC, as well as between carriers of p.I157T and c.1100delC (Table 1). Associations between the mutations and clinico-pathological characteristics were tested with study-stratified Cochran-Mantel-Haenszel test (mantelhaen.test for categorical characteristics and CMHtest for ordinal characteristics). The category of missing data was not included in these comparisons. Differences in age at diagnosis were tested by meta-analysis of age distribution in independent studies using a random effects model (metacont).
Survival analyses were performed using the Cox regression  as implemented in Stata (Stata/SE 10.1 for Windows, StataCorp LP, College Station, TX, USA) comparing CHEK2 mutation carriers and non-carriers, as described above. Study subjects were considered to become at risk at the time of their first invasive breast cancer diagnosis. The data did not consist entirely of incident cases. Therefore, in order to avoid bias caused by late enrollment, we implemented a method called left censoring, which has been proven to provide robust survival estimates for data, which includes also prevalent cases . Survival analysis endpoints included death of any cause, breast cancer-associated death, distant metastasis relapse, locoregional relapse and second breast cancer. Patients were censored at the end of their follow-up period or at the latest 15 years after the initial breast cancer diagnosis in analyses of overall survival and second breast cancer, but at the latest 10 years in analyses of locoregional or distant relapse-free survival as well as in analyses of breast cancer-specific survival. Patients presenting with distant metastases at diagnosis were excluded from the analyses of locoregional relapse-free survival. All analyses were stratified by study.
In addition to univariate analyses, we performed multivariate analyses, which were stratified by study and age category (≤50 years; >50 and ≤70; >70), and adjusted for tumor grade (1, 2 or 3, ordinal), tumor size (1: maximum diameter less than or equal to 20 mm; 2: more than 20 mm and less than or equal to 50 mm; 3 over 50 mm, ordinal), tumor spread in axillary lymph nodes (0 = negative, 1 = positive) and progesterone receptor (PR) status (0 = negative, 1 = positive). ER was not included in the model, because of the non-linear relationship between tumor ER status and patient survival during the 10 years following the diagnosis; patients with ER-negative tumors have a higher risk of dying from breast cancer during the first 5 years after the diagnosis, but the difference in risk between ER-positive and ER-negative tumors levels out after that period [17, 25]. However, since several studies have reported an association between the two CHEK2 mutations and ER-positive disease [12–14] (Table 1), we performed the survival analyses in a subgroup of patients with ER-positive tumors. Only cases with complete data on the pathological markers were included in the multivariate analyses. Univariate survival analyses were performed also in a subgroup of breast cancer patients with lobular tumors, because of the association between p.I157T and lobular breast cancer  (Table 1).
Study subjects for gene expression analysis
Gene expression analyses were performed using a data set of 183 breast tumors from the Helsinki University Hospital (GEO: GSE24450). As described previously, the data set consisted of total RNA samples from 151 tumors from unselected cohorts of breast cancer patients and 32 tumors from additional familial cases hybridized on Illumina HumanHT-12 v3 Expression BeadChips (Illumina Inc., San Diego, CA, USA) [10, 26]. The p.I157T carrier status was defined from peripheral blood samples as described earlier for the BCAC study ‘HEBCS’ (Additional file 1: Table S1). Ten patients were germ line p.I157T carriers and 162 were non-carriers, of which six carried germ line c.1100delC. The c.1100delC carrier tumors were included in the analyses as non-I157T carriers. The p.I157T genotype information was not available for 11 study subjects. These were included in the molecular subtype analysis, but not in differential gene expression or gene set enrichment analysis. The clinico-pathologic characteristics of the 183 tumors are provided in Additional file 1: Table S4.
Gene expression analysis
Gene expression data quality control and quantile normalization was performed in the Bioconductor  as described earlier . Data analyses were performed in R version 3.0.2 and Bioconductor packages genefu , limma [29, 30] and geneplotter . Probes not mapping to any current Entrez Gene entities (GRCh38.p2) were excluded, resulting in a filtered data set of 20,145 genes.
For determining the intrinsic molecular subtypes, expression data of the fifty PAM50 signature genes was extracted from the filtered data set, median centered and standardized per gene by dividing with the standard deviation of the gene’s expression values. Intrinsic subtypes were defined by Pearson correlation between tumors and the luminal A, luminal B, human epidermal growth factor receptor 2 (Her2)-enriched, basal-like and normal-like centroids as implemented in the genefu package [28, 32]. Hierarchical clustering was performed using the Ward’s method . As a comparison to the subtype classification by gene expression, we used the surrogate clinico-pathologic markers to define the subtypes following the St Gallen 2013 criteria (luminal A: ER+, PR+, Her2-, Ki-67-; luminal B (three marker combinations): ER+, PR-, Her2- or ER+, Her2-, Ki-67+ or ER+, Her2+; basal: ER-, PR-, Her2-; Her2 overexpressing: ER-, PR-, Her2+) .
For analysis of differential gene expression the data was filtered by including only genes with highest variation in expression levels over the entire data set (st. dev. ≥ 0.75, 1852 genes). The samples from p.I157T carriers were compared to samples from non-carriers with a moderated t test adjusting for ER, tumor protein 53 (TP53) and Ki-67 protein expression (positive/negative), tumor grade (1, 2, 3, ordinal) as well as histological type (lobular/other). The adjusting covariates were selected from features tabulated in Additional file 1: Table S4 as the most significant factors (p < 0.001) explaining variation in the expression of the 1852 genes as summarized by the first four principal components. Additionally, lobular histologic type was included to avoid bias caused by the association between the p.I157T and lobular type. Data on at least one of the adjusting variables was missing for 12 tumor samples and thus the differential gene expression and gene set enrichment analyses were performed with a set of 160 (ten p.I157T and 150 non-carrier) tumor samples and 1852 genes. Genes with p values below 0.01 were considered to be associated with p.I157T. These were included in a functional enrichment analysis performed using the DAVID functional annotation tool . Functional annotations with Benjamini-Hochberg  corrected p values below 0.01 were considered to be significantly enriched.
Gene set enrichment analysis (GSEA) was performed using a java application available at http://software.broadinstitute.org/gsea following the instructions of the user guide . For the GSEA analysis, the 1852 genes were ranked according to a score calculated as the product of log2(fold change) and log10(p value) from comparisons of p.I157T carrier and non-carrier tumors as described above. All gene sets available at the Molecular Signatures Database (MSigDB) v5.0  were included in the analyses. The p values were corrected for false discovery rate for all other gene sets but the gene sets originating from single publications (‘CGP: chemical and genetic perturbations’ database), which were corrected for the family-wise error rate. Gene sets with corrected p value below 0.05 were considered to be significantly enriched in the p.I157T carrier tumors.
Results and discussion
Our findings from extensive analyses of breast tumor phenotypes and patient survival underline a fundamental difference in breast cancers of the carriers of two CHEK2 mutations, p.I157T and c.1100delC. Significant differences were found in tumor grade and histopathological type as well as in patient survival of p.I157T and c.1100delC carriers, whereas no difference was seen in tumor subtypes: ER+, PR+, Her2- disease was the most common type for carriers of both mutations.
Association of p.I157T with clinico-pathological markers
In our analyses p.I157T was associated with low grade as well as several other markers, which have previously been associated with good prognosis (Table 1). Our analyses confirmed the previously reported associations between p.I157T and ER-positive or lobular breast cancer . Also mixed (ductal and lobular) and tubular histological types were more frequent in p.I157T carrier tumors. Both of ER-positive and lobular tumor types are associated with good short-term prognosis, but increased risk of disease progression in the long run [25, 39]. Furthermore, p.I157T was associated with PR-positive and TP53-negative breast cancer. PR expression is a marker for good prognosis for ER-positive breast cancer and it has been suggested as a surrogate marker separating luminal A and luminal B subtypes according to immunohistochemical data [34, 40, 41]. TP53 immunohistochemical staining is considered to be indicative of somatic TP53 mutations. Strong TP53 staining suggests the presence of stabilizing mutations (primarily missense), whereas absence of staining indicates typically a protein-truncating mutation (nonsense or frameshift), and weak staining a wild-type functional TP53. Both strong and completely negative TP53 staining have been associated with poor prognosis in comparison to weak staining [42–44]. The sensitivity of the assays used in this study did not enable differentiation between normal, low or absent TP53 expression. Therefore, we used binary classification of TP53 immunohistochemical data, the positive category corresponding to high expression (strong staining) and negative category to low expression (Additional file 1: Table S3). Noteworthy, the loss-of-function mutations associated with absent TP53 staining are relatively rare in breast cancer: these are seen in less than 5 % of all tumors [45, 46]. Therefore, it is likely that most of the tumors in the category ‘negative’ (Table 1) represented tumors with wild-type TP53. However, compromised CHEK2 function as a result of the p.I157T mutation could be another way for TP53 silencing as CHEK2 is among the key upstream activators of TP53 .
Like p.I157T, also c.1100delC was associated with ER-positive and PR-positive disease in our data set (Table 1). Furthermore, TP53-positive tumors were slightly less often observed in c.1100delC carriers than in non-carriers, even though the difference was not statistically significant. Significant differences in clinico-pathological features associated with the two CHEK2 mutations were seen in grade and histological type, as the c.1100delC carrier tumors resembled more non-carrier tumors (Table 1).
Breast cancer subtypes
We investigated the I1157T-associated molecular breast cancer subtypes by applying the St Gallen 2013 criteria for immunohistochemical markers  on the BCAC data set as well as St Gallen 2013 and the PAM50 classifier  on the gene expression data set of 183 breast tumors. The subtype classification of the BCAC study subjects relied on the available immunohistochemical markers, ER, PR and Her2. We found both p.I157T and c.1100delC carrier tumors to be predominantly ER+, ER+, Her2-, suggestive of good prognosis ER+ tumors or the luminal A subtype (Table 1) . Also the frequency of ER+ subtypes linked to poor prognosis (ER+, PR-, Her2-; ER+, Her2+), referred to as luminal B , were more common for CHEK2 mutation carriers than for non-carriers. This confirmed previous reports with regard to p.I157T , but was not consistent with previous reports on c.1100delC-associated tumor subtypes [47, 48]. However, the difference between our findings and these reports may have arisen from different overall cohort compositions or from differing classification methods, as the guidelines for subtype classification have changed over the years.
Subtype classifications of the 183 tumors according to gene expression data and immunohistochemical markers were partly contradictory (Table 2). Similar inconsistencies between gene expression-based classification and the surrogate immunohistochemical markers have been reported previously for other data sets [40, 49]. Overall, the division between basal and luminal appeared rather consistent: only 17 (9 %) of the 183 tumors were classified differentially across the luminal-basal axis. PAM50  classified three of the p.I157T carrier tumors as luminal A, two as luminal B and two as basal. Three lobular tumors were classified as normal-like (Table 2). This kind of a misclassification has been reported to be typical for lobular tumors due to their infiltrating growth pattern, which causes the tumor sample to consist of unusually high proportion of non-cancerous stromal cells . St Gallen 2013 criteria classified these normal-like tumors as luminal (ER+). Furthermore, in unsupervised hierarchical clustering of the 183 tumor samples based on expression of the PAM50 signature genes (Fig. 1), two of the normal-like p.I157T tumors (HEL_045 and HEL_174) clustered within the luminal A branch suggesting that luminal A could be their true molecular subtype. In summary, luminal A appeared to be the most common subtype for p.I157T carrier tumors in the gene expression data concordantly with the findings in the BCAC data.
P.I157T carriers had better prognosis than the c.1100delC carriers with regard to overall or breast cancer-specific survival (Table 3a and b). This difference was possibly due to the poor survival associated with c.1100delC as reported previously by several studies [13, 14, 17, 50]. No statistically significant difference in overall or breast cancer-specific survival was seen between p.I157T carriers and non-carriers. Hazard ratios in the analyses of subgroups of patients with ER-positive or lobular tumors were comparable to those of the main analyses (Table 3b and c).
Noteworthy, the c.1100delC carriers included here were only a subset of the study subjects included in the previous report by Weischer and colleagues on survival of c.1100delC carriers of the BCAC studies . This was because the individual BCAC studies, which did not provide sufficient number of p.I157T carriers, were excluded from these analyses. Thus, the lack of statistical significance in some comparisons of survival difference between c.1100delC carriers and non-carriers (Table 3a and b) probably only reflected limited power due to low number of c.1100delC carriers, since the hazard ratios were always consistent with the previous report.
The different prognoses associated with p.I157T and c.1100delC possibly reflect their difference in molecular level severity of functional consequences. Therefore, it would be tempting to assume that the prognosis of all carriers of the truncating mutations would be similar to the prognosis of c.1100delC carriers. However, a recent Polish study combining three truncating CHEK2 founder mutations found no difference between mutation carrier and non-carrier survival . Some part of the conflicting findings could be explained by different patient selection: in the Polish study all patients had been diagnosed before 50 years of age, whereas here and in Weischer et al.  also postmenopausal patients were included in the analyses. Another potential explanation could be mutation-specific survival effects. As the Polish study combined in the analyses three different truncating mutations, the c.1100delC specific effects could have been masked, since it is the least common of the three truncating CHEK2 mutations in Polish population . Similarly as here, the Polish study reported no significant difference in survival of the p.I157T carriers and non-carriers .
The hazard ratios for locoregional relapse and second breast cancer (91 % contralateral, 9 % ipsilateral) associated with p.I157T and c.1100delC were close to the mutations’ relative risk estimates of primary breast cancer (Table 3) [5, 17, 53]. The marginally significant increased risk of locoregional relapse associated with p.I157T in the adjusted analyses (hazard ratio 1.62 [0.99 - 2.66], p value 0.056) warrants further studies, but could merely reflect the baseline risk associated with p.I157T: some of the local recurrences could represent new cancers arising during the 10-year follow-up. The risk of locoregional relapse for c.1100delC carriers was elevated in the univariate analysis but leveled out in the adjusted analysis.
P.I157T associated differentially expressed genes
In order to investigate the molecular biology of p.I157T carrier tumors and to identify potential tumor-driving events and pathways, we performed an analysis of differential gene expression and subsequent functional enrichment analysis comparing ten p.I157T to 150 non-carrier tumors. We found 21 genes to be differentially expressed between p.I157T and non-carrier tumors. All of these had higher expression in the p.I157T carrier tumors (Table 4). When the 160 tumor samples were clustered according to expression of these 21 genes, the p.I157T tumors did not form a distinct cluster (Fig. 2), suggesting that high expression of these genes is not exclusive of the p.I157T mutation carrier tumors, but typical for a subgroup of breast tumors including the mutation carrier tumors. Tumors with different intrinsic subtypes appeared to be dispersed across all branches, similarly as the c.1100delC carriers suggesting that in this data set the c.1100delC carrier tumors would not be similar to the p.I157T tumors.
Enrichment of features associated with lobular breast cancer
The list of 21 differentially expressed genes contained seven collagen genes (Table 4), which were a major driver in the functional enrichment analysis. The enriched annotations from DAVID  analysis included characteristics of the collagen family and their related functions such as ‘focal adhesion’, ‘extracellular matrix (ECM) organization’ and ‘ECM-receptor interaction’ (Additional file 1: Table S5). Similar results were obtained from the GSEA analysis (Additional file 1: Table S6, Additional file 2: Figure S1). Since collagens are usually expressed by stromal fibroblasts, the findings may suggest that infiltrating growth pattern, typical for lobular tumors , could be more common also for non-lobular p.I157T carrier tumors than for the non-carrier tumors. Further support for this hypothesis came from the GSEA, which showed cadherin 1 (CDH1) target genes to be significantly enriched among genes, whose expression was lower in p.I157T than in non-carrier tumors (Additional file 1: Table S6). CDH1 silencing is generally considered as a defining characteristic of lobular tumors and it is often caused by somatic mutations targeting the CDH1 gene itself . However, since the differential gene expression analysis, which was also the basis for the ranked gene list used as an input to GSEA, was adjusted for the lobular tumor type, the impact of the diagnosed lobular cancers on these findings should have been minimal. CDH1 gene expression was lower in p.I157T carriers tumors in the adjusted analysis (log2 fold change -1.12, p value 0.03, Fig. 3), but it did not exceed the preset threshold for significance. Previously, we have reported CDH1 mRNA expression to be higher in c.1100delC carrier than in non-carriers tumors . Therefore, CDH1 expression appears to be yet another factor, which is not shared by breast tumors from carriers of the two CHEK2 mutations, p.157 T and c.1100delC, and possibly reflects somatic changes, which have taken place during the clonal evolution of the p.I157T carrier tumors . Taken together, these results suggest that besides the fact that the lobular tumors are more common among p.I157T carriers and non-carriers, the association between p.I157T and lobular features could be even stronger than what is suggested by the diagnosed histological types.
Enrichment of cancer associated gene signatures
In the GSEA analysis, several independent MSigDB  gene signatures related to epithelial-to-mesenchymal transition (EMT) [55–57], stromal stem cells  or invasive behavior [59, 60] were enriched at the top of the gene list with higher expression in p.I157T carrier tumors than in non-carrier tumors (Additional file 1: Table S6). These observations may reflect higher stromal content of the p.I157T carrier tumor samples, as the samples were not prepared at a single cell level. However, to prevent such confounding effects the tumor sample sections were selected by an experienced breast cancer pathologist. Furthermore, the above mentioned MSigDB signatures originated from carefully designed experiments tailored to detect the true signal from cancerous epithelial cells and to escape the effects of non-cancerous stromal cells. The enrichment of these signatures may suggest that the p.I157T carrier tumors have an intrinsically invasive nature. However, this should have been reflected into poor prognosis for the p.I157T carriers, which we did not see in the survival analyses. On the other hand, it is possible that higher state of differentiation of the tumor cells suggested by low grade accompanied with the invasive nature can be seen in the prognosis only in the long run, and within the 10-year follow-up period is only reflected in the slightly elevated risk of local recurrence. All in all, these observations deserve further studies before any definitive conclusions can be made.
In addition to CDH1, tumor suppressor retinoblastoma 1 (RB1) appeared as a potential gene expression regulator, whose activity was reduced in p.I157T carrier tumors in comparison to non-carrier tumors (Additional file 1: Table S6, Additional file 1: Figure S1). RB1 and its direct downstream target E2F-1 are both targets of the CHEK2 protein [61, 62]. Thus, the differential expression of the RB1 target genes possibly reflects compromised CHEK2 function in the p.I157T carrier tumors.
Noteworthy, the two differential gene expression studies on c.1100delC carrier tumors have reported enrichment of genes of WNT and FGF pathways [10, 47], which regulate the growth and differentiation of normal breast epithelium [63–66]. Among the p.I157T-associated differentially expressed genes we did not see enrichment of any growth factor pathway. These notions on differences in gene expression signatures are more descriptive than definitive by nature, but they further emphasize intrinsic biological differences between p.I157T and 1100deC carrier tumors.
Based on our analyses, breast cancers of p.I157T and c.1100delC CHEK2 mutation carriers differ in disease severity as seen especially in differences in tumor grade and patient survival, as well as in intrinsic biological features as seen in differences in histological type and gene expression profiles. Thus, it appears that even though both mutations have been proven to compromise the protein function [6, 9, 11], they have different consequences on the disease phenotype, and prognostic findings based on one mutation cannot be generalized to the other. Furthermore, our results raise a hypothesis that the increased risk of locoregional relapse for p.I157T carriers could be caused by intrinsically invasive nature of the tumor cells. Future studies with longer follow-up are needed to test this hypothesis.
Breast Cancer Association Consortium
Checkpoint kinase 2
Gene set enrichment analysis
Human epidermal growth factor receptor 2
Tumor protein 53
CHEK2 Breast Cancer Case-Control Consortium. CHEK2*1100delC and susceptibility to breast cancer: a collaborative analysis involving 10,860 breast cancer cases and 9,065 controls from 10 studies. Am J Hum Genet. 2004;74(6):1175–82.
Johnson N, Fletcher O, Naceur-Lombardelli C, dos Santos SI, Ashworth A, Peto J. Interaction between CHEK2*1100delC and other low-penetrance breast-cancer susceptibility genes: a familial study. Lancet. 2005;366(9496):1554–7.
Cybulski C, Wokolorczyk D, Jakubowska A, Huzarski T, Byrski T, Gronwald J, et al. Risk of breast cancer in women with a CHEK2 mutation with and without a family history of breast cancer. J Clin Oncol. 2011;29(28):3747–52.
Kilpivaara O, Vahteristo P, Falck J, Syrjakoski K, Eerola H, Easton D, et al. CHEK2 variant I157T may be associated with increased breast cancer risk. Int J Cancer. 2004;111(4):543–7.
Nevanlinna H, Bartek J. The CHEK2 gene and inherited breast cancer susceptibility. Oncogene. 2006;25(43):5912–9.
Roeb W, Higgins J, King MC. Response to DNA damage of CHEK2 missense mutations in familial breast cancer. Hum Mol Genet. 2012;21(12):2738–44.
Le Calvez-Kelm F, Lesueur F, Damiola F, Vallee M, Voegele C, Babikyan D, et al. Rare, evolutionarily unlikely missense substitutions in CHEK2 contribute to breast cancer susceptibility: results from a breast cancer family registry case-control mutation-screening study. Breast Cancer Res. 2011;13(1):R6.
Bell DW, Kim SH, Godwin AK, Schiripo TA, Harris PL, Haserlat SM, et al. Genetic and functional analysis of CHEK2 (CHK2) variants in multiethnic cohorts. Int J Cancer. 2007;121(12):2661–7.
Anczukow O, Ware MD, Buisson M, Zetoune AB, Stoppa-Lyonnet D, Sinilnikova OM, et al. Does the nonsense-mediated mRNA decay mechanism prevent the synthesis of truncated BRCA1, CHK2, and p53 proteins? Hum Mutat. 2008;29(1):65–73.
Muranen TA, Greco D, Fagerholm R, Kilpivaara O, Kampjarvi K, Aittomaki K, et al. Breast tumors from CHEK2 1100delC-mutation carriers: genomic landscape and clinical implications. Breast Cancer Res. 2011;13(5):R90.
Cai Z, Chehab NH, Pavletich NP. Structure and activation mechanism of the CHK2 DNA damage checkpoint kinase. Mol Cell. 2009;35(6):818–29.
Kilpivaara O, Bartkova J, Eerola H, Syrjakoski K, Vahteristo P, Lukas J, et al. Correlation of CHEK2 protein expression and c.1100delC mutation status with tumor characteristics among unselected breast cancer patients. Int J Cancer. 2005;113(4):575–80.
de Bock GH, Schutte M, Krol-Warmerdam EM, Seynaeve C, Blom J, Brekelmans CT, et al. Tumour characteristics and prognosis of breast cancer patients carrying the germline CHEK2*1100delC variant. J Med Genet. 2004;41(10):731–5.
Schmidt MK, Tollenaar RA, de Kemp SR, Broeks A, Cornelisse CJ, Smit VT, et al. Breast cancer survival and tumor characteristics in premenopausal women carrying the CHEK2*1100delC germline mutation. J Clin Oncol. 2007;25(1):64–9.
Domagala P, Wokolorczyk D, Cybulski C, Huzarski T, Lubinski J, Domagala W. Different CHEK2 germline mutations are associated with distinct immunophenotypic molecular subtypes of breast cancer. Breast Cancer Res Treat. 2012;132(3):937–45.
Huzarski T, Cybulski C, Domagala W, Gronwald J, Byrski T, Szwiec M, et al. Pathology of breast cancer in women with constitutional CHEK2 mutations. Breast Cancer Res Treat. 2005;90(2):187–9.
Weischer M, Nordestgaard BG, Pharoah P, Bolla MK, Nevanlinna H, Van’t Veer LJ, et al. CHEK2*1100delC heterozygosity in women with breast cancer associated with early death, breast cancer-specific death, and increased risk of a second breast cancer. J Clin Oncol. 2012;30(35):4308–16.
Michailidou K, Hall P, Gonzalez-Neira A, Ghoussaini M, Dennis J, Milne RL, et al. Large-scale genotyping identifies 41 new loci associated with breast cancer risk. Nat Genet. 2013;45(4):353–61.
Broeks A, Schmidt MK, Sherman ME, Couch FJ, Hopper JL, Dite GS, et al. Low penetrance breast cancer susceptibility loci are associated with specific breast tumor subtypes: findings from the Breast Cancer Association Consortium. Hum Mol Genet. 2011;20(16):3289–303.
R Development Core Team. R: A language and environment for statistical computing. 2012.
Friendly M. vcdExtra: ‘vcd’ extensions and additions. 2015. R package version 0.6-8.
Schwarzer G. meta: General package for meta-analysis. R package version 4.1-0. 2015.
Cox DR. Regression models and life-tables. J Roy Statist Soc Ser B. 1972;34:187–220.
Azzato EM, Greenberg D, Shah M, Blows F, Driver KE, Caporaso NE, et al. Prevalent cases in observational studies of cancer survival: do they bias hazard ratio estimates? Br J Cancer. 2009;100(11):1806–11.
Colzani E, Liljegren A, Johansson AL, Adolfsson J, Hellborg H, Hall PF, et al. Prognosis of patients with breast cancer: causes of death and effects of time since diagnosis, age, and tumor characteristics. J Clin Oncol. 2011;29(30):4014–21.
Heikkinen T, Greco D, Pelttari LM, Tommiska J, Vahteristo P, Heikkila P, et al. Variants on the promoter region of PTEN affect breast cancer progression and patient survival. Breast Cancer Res. 2011;13(6):R130.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):R80.
Haibe-Kains B, Schroeder M, Bontempi G, Sotirou C, Quackenbush J. genefu: Relevant functions for gene expression analysis, especially in breast cancer. 2013. R package version 1.12.0.
Smyth GK. Limma: linear models for microarray data. In: Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W, editors. Bioinformatics and computational biology solutions using R and Bioconductor. New York: Springer; 2005.
Ritchie ME, Silver J, Oshlack A, Holmes M, Diyagama D, Holloway A, et al. A comparison of background correction methods for two-colour microarrays. Bioinformatics. 2007;23(20):2700–7.
Gentleman R, Biocore. geneplotter: Graphics related functions for Bioconductor; R package version 1.40.0.
Parker JS, Mullins M, Cheang MC, Leung S, Voduc D, Vickery T, et al. Supervised risk predictor of breast cancer based on intrinsic subtypes. J Clin Oncol. 2009;27(8):1160–7.
Ward Jr JH. Hierarchical grouping to optimize an objective function. J Am Stat Assoc. 1963;58:236–44.
Curigliano G, Criscitiello C, Andre F, Colleoni M, Di Leo A. Highlights from the 13th St Gallen International Breast Cancer Conference 2013. Access to innovation for patients with breast cancer: how to speed it up? Ecancermedicalscience. 2013;7:299.
Anonymous. DAVID Bioinformatics Database. 2010. http://david.abcc.ncifcrf.gov/home.jsp.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc. 1995;57(1):289–300.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.
Liberzon A, Subramanian A, Pinchback R, Thorvaldsdottir H, Tamayo P, Mesirov JP. Molecular signatures database (MSigDB) 3.0. Bioinformatics. 2011;27(12):1739–40.
McCart Reed AE, Kutasovic JR, Lakhani SR, Simpson PT. Invasive lobular carcinoma of the breast: morphology, biomarkers and ‘omics’. Breast Cancer Res. 2015;17:12.
Prat A, Cheang MC, Martin M, Parker JS, Carrasco E, Caballero R, et al. Prognostic significance of progesterone receptor-positive tumor cells within immunohistochemically defined luminal A breast cancer. J Clin Oncol. 2013;31(2):203–9.
Dunnwald LK, Rossing MA, Li CI. Hormone receptor status, tumor characteristics, and prognosis: a prospective cohort of breast cancer patients. Breast Cancer Res. 2007;9(1):R6.
Yemelyanova A, Vang R, Kshirsagar M, Lu D, Marks MA, Shih I, et al. Immunohistochemical staining patterns of p53 can serve as a surrogate marker for TP53 mutations in ovarian carcinoma: an immunohistochemical and nucleotide sequencing analysis. Mod Pathol. 2011;24(9):1248–53.
Watanabe G, Ishida T, Furuta A, Takahashi S, Watanabe M, Nakata H, et al. Combined immunohistochemistry of PLK1, p21, and p53 for predicting TP53 status: an independent prognostic factor of breast cancer. Am J Surg Pathol. 2015;39(8):1026–34.
Valgardsdottir R, Tryggvadottir L, Steinarsdottir M, Olafsdottir K, Jonasdottir S, Jonasson JG, et al. Genomic instability and poor prognosis associated with abnormal TP53 in breast carcinomas. Molecular and immunohistochemical analysis. APMIS. 1997;105(2):121–30.
Olivier M, Langerod A, Carrieri P, Bergh J, Klaar S, Eyfjord J, et al. The clinical value of somatic TP53 gene mutations in 1,794 patients with breast cancer. Clin Cancer Res. 2006;12(4):1157–67.
Boyle DP, McArt DG, Irwin G, Wilhelm-Benartzi CS, Lioe TF, Sebastian E, et al. The prognostic significance of the aberrant extremes of p53 immunophenotypes in breast cancer. Histopathology. 2014;65(3):340–52.
Nagel JH, Peeters JK, Smid M, Sieuwerts AM, Wasielewski M, de Weerd V, et al. Gene expression profiling assigns CHEK2 1100delC breast cancers to the luminal intrinsic subtypes. Breast Cancer Res Treat. 2012;132(2):439–48.
Sorlie T, Tibshirani R, Parker J, Hastie T, Marron JS, Nobel A, et al. Repeated observation of breast tumor subtypes in independent gene expression data sets. Proc Natl Acad Sci U S A. 2003;100(14):8418–23.
Cheang MC, Martin M, Nielsen TO, Prat A, Voduc D, Rodriguez-Lescure A, et al. Defining breast cancer intrinsic subtypes by quantitative receptor expression. Oncologist. 2015;20(5):474–82.
Kriege M, Hollestelle A, Jager A, Huijts PE, Berns EM, Sieuwerts AM, et al. Survival and contralateral breast cancer in CHEK2 1100delC breast cancer patients: impact of adjuvant chemotherapy. Br J Cancer. 2014;111(5):1004–13.
Huzarski T, Cybulski C, Wokolorczyk D, Jakubowska A, Byrski T, Gronwald J, et al. Survival from breast cancer in patients with CHEK2 mutations. Breast Cancer Res Treat. 2014;144(2):397–403.
Cybulski C, Wokolorczyk D, Huzarski T, Byrski T, Gronwald J, Gorski B, et al. A deletion in CHEK2 of 5,395 bp predisposes to breast cancer in Poland. Breast Cancer Res Treat. 2007;102(1):119–22.
Weischer M, Bojesen SE, Ellervik C, Tybjaerg-Hansen A, Nordestgaard BG. CHEK2*1100delC genotyping for clinical assessment of breast cancer risk: meta-analyses of 26,000 patient cases and 27,000 controls. J Clin Oncol. 2008;26(4):542–8.
Dennis Jr G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, et al. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003;4(5):3.
Anastassiou D, Rumjantseva V, Cheng W, Huang J, Canoll PD, Yamashiro DJ, et al. Human cancer cells express Slug-based epithelial-mesenchymal transition gene expression signature obtained in vivo. BMC Cancer. 2011;11:529.
Charafe-Jauffret E, Ginestier C, Monville F, Finetti P, Adelaide J, Cervera N, et al. Gene expression profiling of breast cell lines identifies potential new basal markers. Oncogene. 2006;25(15):2273–84.
Jechlinger M, Grunert S, Tamir IH, Janda E, Ludemann S, Waerner T, et al. Expression profiling of epithelial plasticity in tumor progression. Oncogene. 2003;22(46):7155–69.
Boquest AC, Shahdadfar A, Fronsdal K, Sigurjonsson O, Tunheim SH, Collas P, et al. Isolation and transcription profiling of purified uncultured human stromal stem cells: alteration of gene expression after in vitro cell culture. Mol Biol Cell. 2005;16(3):1131–41.
Schuetz CS, Bonin M, Clare SE, Nieselt K, Sotlar K, Walter M, et al. Progression-specific genes identified by expression profiling of matched ductal carcinomas in situ and invasive breast tumors, combining laser capture microdissection and oligonucleotide microarray analysis. Cancer Res. 2006;66(10):5278–86.
Sung SY, Hsieh CL, Law A, Zhau HE, Pathak S, Multani AS, et al. Coevolution of prostate cancer and bone stroma in three-dimensional coculture: implications for cancer growth and metastasis. Cancer Res. 2008;68(23):9996–10003.
Inoue Y, Kitagawa M, Taya Y. Phosphorylation of pRB at Ser612 by Chk1/2 leads to a complex between pRB and E2F-1 after DNA damage. EMBO J. 2007;26(8):2083–93.
Stevens C, Smith L, La Thangue NB. Chk2 activates E2F-1 in response to DNA damage. Nat Cell Biol. 2003;5(5):401–9.
Krejci P, Aklian A, Kaucka M, Sevcikova E, Prochazkova J, Masek JK, et al. Receptor tyrosine kinases activate canonical WNT/beta-catenin signaling via MAP kinase/LRP6 pathway and direct beta-catenin phosphorylation. PLoS One. 2012;7(4):e35826.
Incassati A, Chandramouli A, Eelkema R, Cowin P. Key signaling nodes in mammary gland development and cancer: β-catenin. Breast Cancer Res. 2010;12(6):213.
Pond AC, Herschkowitz JI, Schwertfeger KL, Welm B, Zhang Y, York B, et al. Fibroblast growth factor receptor signaling dramatically accelerates tumorigenesis and enhances oncoprotein translation in the mouse mammary tumor virus-Wnt-1 mouse model of breast cancer. Cancer Res. 2010;70(12):4868–79.
Pond AC, Bin X, Batts T, Roarty K, Hilsenbeck S, Rosen JM. Fibroblast growth factor receptor signaling is essential for normal mammary gland development and stem cell function. Stem Cells. 2013;31(1):178–89.
We thank all the individuals who took part in these studies and all the researchers, clinicians, technicians and administrative staff who have enabled this work to be carried out. Especially we would like to acknowledge Kirsimari Aaltonen, Karl von Smitten, Sofia Khan, Tuomas Heikkinen, Irja Erkkilä; The SEARCH and EPIC teams; The Swedish Medical Research Counsel; Petra Seibold, Dieter Flesch-Janys, Judith Heinz, Nadia Obi, Alina Vrieling, Sabine Behrens, Ursula Eilber, Muhabbet Celik, Til Olchers and Stefan Nickels; staff and participants of the Copenhagen General Population Study. For the excellent technical assistance: Dorthe Uldall Andersen, Maria Birna Arnadottir, Anne Bank, Dorthe Kjeldgård Hansen; Louise Brinton, Mark Sherman, Neonila Szeszenia-Dabrowska, Beata Peplonska, Witold Zatonski, Pei Chao, Michael Stagner; Teresa Selander, Nayana Weerasooriya; Meeri Otsukka, Kari Mononen; Eija Myöhänen, Helena Kemiläinen; and Michael Bremer. The iCOGS study would not have been possible without the contributions of the following: Andrew Berchuck (OCAC), Rosalind A. Eeles, Ali Amin Al Olama, Zsofia Kote-Jarai and Sara Benlloch (PRACTICAL), Georgia Chenevix-Trench, Antonis Antoniou, Lesley McGuffog and Ken Offit (CIMBA), Andrew Lee, and Ed Dicks, Craig Luccarini and the staff of the Centre for Genetic Epidemiology Laboratory, Javier Benitez, Anna Gonzalez-Neira and the staff of the CNIO genotyping unit, Jacques Simard and Daniel C. Tessier, Francois Bacot, Daniel Vincent, Sylvie LaBoissière and Frederic Robidoux and the staff of the McGill University and Génome Québec Innovation Centre, Sune F. Nielsen and the staff of the Copenhagen DNA laboratory, and Julie M. Cunningham, Sharon A. Windebank, Christopher A. Hilker, Jeffrey Meyer and the staff of Mayo Clinic Genotyping Core Facility.
BCAC is funded by Cancer Research UK [C1287/A10118, C1287/A12014] and by the European Community’s Seventh Framework Programme under grant agreement number 223175 (grant number HEALTH-F2-2009-223175) (COGS); The work of the BBCC was partly funded by ELAN-Fond of the University Hospital of Erlangen; The HEBCS was financially supported by the Helsinki University Central Hospital Research Fund, Academy of Finland (266528), the Finnish Cancer Society, The Nordic Cancer Union and the Sigrid Juselius Foundation, The work of TAM has been supported by Finnish Cultural Foundation, Orion Research Foundation, Cancer Society of Finland and Ida Montin Foundation; SEARCH is funded by a programme grant from Cancer Research UK [C490/A10124] and supported by the UK National Institute for Health Research Biomedical Research Centre at the University of Cambridge; The pKARMA study was supported by Märit and Hans Rausings Initiative Against Breast Cancer; The MARIE study was supported by the Deutsche Krebshilfe e.V. [70-2892-BR I, 106332, 108253, 108419], the Hamburg Cancer Society, the German Cancer Research Center (DKFZ) and the Federal Ministry of Education and Research (BMBF) Germany [01KH0402]; The CGPS was supported by the Chief Physician Johan Boserup and Lise Boserup Fund, the Danish Medical Research Council and Herlev Hospital; The MCBCS was supported by the NIH grants CA192393, CA116167, CA176785 an NIH Specialized Program of Research Excellence (SPORE) in Breast Cancer [CA116201], and the Breast Cancer Research Foundation; The PBCS was funded by Intramural Research Funds of the National Cancer Institute, Department of Health and Human Services, USA; The SASBAC study was supported by funding from the Agency for Science, Technology and Research of Singapore (A*STAR), the US National Institute of Health (NIH) and the Susan G. Komen Breast Cancer Foundation; The Ontario Familial Breast Cancer Registry (OFBCR) was supported by grant UM1 CA164920 from the National Cancer Institute (USA). The content of this manuscript does not necessarily reflect the views or policies of the National Cancer Institute or any of the collaborating centers in the Breast Cancer Family Registry (BCFR), nor does mention of trade names, commercial products, or organizations imply endorsement by the USA Government or the BCFR; The OBCS was supported by research grants from the Finnish Cancer Foundation, the Academy of Finland (grant number 250083, 122715 and Center of Excellence grant number 251314), the Sigrid Juselius Foundation, the University of Oulu, the University of Oulu Support Foundation and the special Governmental EVO funds for Oulu University Hospital-based research activities; The KBCP was financially supported by the special Government Funding (EVO) of Kuopio University Hospital grants, Cancer Fund of North Savo, the Finnish Cancer Organizations, and by the strategic funding of the University of Eastern Finland; Financial support for KARBAC was provided through the regional agreement on medical training and clinical research (ALF) between Stockholm County Council and Karolinska Institutet, the Swedish Cancer Society, The Gustav V Jubilee foundation and and Bert von Kantzows foundation; The SZBCS was supported by Grant PBZ_KBN_122/P05/2004; The HABCS study was supported by an intramural grant from Hannover Medical School; Funding for the iCOGS infrastructure came from: the European Community’s Seventh Framework Programme under grant agreement n° 223175 (HEALTH-F2-2009-223175) (COGS), Cancer Research UK (C1287/A10118, C1287/A 10710, C12292/A11174, C1281/A12014, C5047/A8384, C5047/A15007, C5047/A10692, C8197/A16565), the National Institutes of Health (CA128978) and Post-Cancer GWAS initiative (1U19 CA148537, 1U19 CA148065 and 1U19 CA148112 - the GAME-ON initiative), the Department of Defence (W81XWH-10-1-0341), the Canadian Institutes of Health Research (CIHR) for the CIHR Team in Familial Risks of Breast Cancer, Komen Foundation for the Cure, the Breast Cancer Research Foundation, and the Ovarian Cancer Research Fund.
Availability of data and materials
The dataset supporting the conclusions of this article is available in the Gene Expression Omnibus (GEO) repository, GEO: GSE24450, http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/geo/query/acc.cgi?acc=GSE24450.
TAM, HN, MKS, CB, TD, AJ and PDP conceived and designed the study, ran the data analyses and interpreted the results. HN, PDP, PH and DFE provided financial support. MKB, QW and DFE provided administrative support. All authors provided study material and patients, participated in collecting and assembling the data as well as in manuscript writing. All authors approved the final version of the manuscript.
Dr. Peter Fasching has had a consultant or advisory relationship with Roche and Novartis, has received honoraria from Novartis and Genomic Health, and has received research funding from Novartis and Amgen.
Other authors declare no conflicts of interest.
Ethics approval and consent to participate
All participating studies were approved by their institutional review committees. All participants provided written informed consent or did not object to the secondary use of their tissue and data following country-specific regulations.
Description, genotyping methods and references for individual studies as reported by the studies. Table S2. Genotype and follow-up data availability in individual studies. Table S3. Sources and scoring of TP53 immunohistochemistry data used in this study. Table S4. Pathological characteristics of 183 breast tumors used in the gene expression analysis. Table S5. Functional annotations enriched in the 21 differentially expressed genes. Table S6. Enriched gene sets at the high and low edges of a list of 1852 genes ranked according to differences between p.I157T and non-carrier breast tumors. (DOCX 72 kb)
Enrichment plots of the high ranking gene sets from the gene set enrichment analyses (GSEA). (PDF 364 kb)
About this article
Cite this article
Muranen, T.A., Blomqvist, C., Dörk, T. et al. Patient survival and tumor characteristics associated with CHEK2:p.I157T – findings from the Breast Cancer Association Consortium. Breast Cancer Res 18, 98 (2016) doi:10.1186/s13058-016-0758-5
- Breast cancer
- Gene expression