Skip to main content

Distinct high resolution genome profiles of early onset and late onset colorectal cancer integrated with gene expression data identify candidate susceptibility loci

Abstract

Background

Estimates suggest that up to 30% of colorectal cancers (CRC) may develop due to an increased genetic risk. The mean age at diagnosis for CRC is about 70 years. Time of disease onset 20 years younger than the mean age is assumed to be indicative of genetic susceptibility. We have compared high resolution tumor genome copy number variation (CNV) (Roche NimbleGen, 385 000 oligo CGH array) in microsatellite stable (MSS) tumors from two age groups, including 23 young at onset patients without known hereditary syndromes and with a median age of 44 years (range: 28-53) and 17 elderly patients with median age 79 years (range: 69-87). Our aim was to identify differences in the tumor genomes between these groups and pinpoint potential susceptibility loci. Integration analysis of CNV and genome wide mRNA expression data, available for the same tumors, was performed to identify a restricted candidate gene list.

Results

The total fraction of the genome with aberrant copy number, the overall genomic profile and the TP53 mutation spectrum were similar between the two age groups. However, both the number of chromosomal aberrations and the number of breakpoints differed significantly between the groups. Gains of 2q35, 10q21.3-22.1, 10q22.3 and 19q13.2-13.31 and losses from 1p31.3, 1q21.1, 2q21.2, 4p16.1-q28.3, 10p11.1 and 19p12, positions that in total contain more than 500 genes, were found significantly more often in the early onset group as compared to the late onset group. Integration analysis revealed a covariation of DNA copy number at these sites and mRNA expression for 107 of the genes. Seven of these genes, CLC, EIF4E, LTBP4, PLA2G12A, PPAT, RG9MTD2, and ZNF574, had significantly different mRNA expression comparing median expression levels across the transcriptome between the two groups.

Conclusions

Ten genomic loci, containing more than 500 protein coding genes, are identified as more often altered in tumors from early onset versus late onset CRC. Integration of genome and transcriptome data identifies seven novel candidate genes with the potential to identify an increased risk for CRC.

Background

Less than five percent of all patients diagnosed with colorectal cancers (CRC) carry known genetic germline alterations that predispose to the disease [1]. However, it has been estimated that up to 30% of all CRC patients may carry a genetic risk as suggested by young age at onset, multiple tumors in the same patient, and an excess of individuals with CRC within a family [2, 3]. Many studies have tried to identify some of these genetic risk factors, and several recent genome-wide association studies (GWAS) have pinpointed SNP loci on chromosome arms 8q, 10p, 11q, 14q, 15q, 16q, 18q, 19q, and 20p to be associated with CRC [4–10]. Furthermore, a study by Mourra et al. [11] showed that microsatellite loci within chromosome arm 14q, known to be deleted in about 30% of all colorectal cancers, were more frequently lost in tumors from early onset patients.

TP53 mutations and genomic copy number alterations, as well as other somatic genetic and epigenetic alterations, have been shown to accumulate with the adenoma-carcinoma development in CRC [12–17]. Copy number alterations are typically identified using cytogenetic techniques as G-banding, chromosome-based comparative genomic hybridization (cCGH) and array-CGH (aCGH) [18]. For the most common chromosomal aberrations, which include gains at 7q, 7p, 8q, 11q, 13q, and 20q and losses from 1p, 4p, 4q, 8p, 14q, 15q, 17p, and 18, the time of occurrence in the adenoma-carcinoma sequence has been suggested [16, 17, 19]. In addition, cCGH has been used to identify DNA sequences that contain predisposing genes, e.g. changes in chromosome 19 found in tumors from patients with Peutz-Jeghers syndrome, led to the identification of STK11 as the predisposition gene [20]. Other susceptibility loci and genes have been suggested for CRC, typically based on linkage analyses and genome wide SNP analyses [10, 21–28].

Array-CGH allows for increased resolution, improves the chromosome dependent method, and thus facilitates detection of small aberrations and fine-tunes the accuracy of breakpoint determination [29]. In order to identify somatic differences and potential susceptibility loci for CRC, we have compared high resolution (385 000 oligo probe array) DNA copy number profile and TP53 mutation status in carcinomas from late onset and early onset patients without known hereditary CRC syndromes. These data have further been integrated with corresponding gene expression data for each patient.

Methods

Patients and tumor samples

Forty patients diagnosed with CRC, were included in the study. Patient gender and age, and tumor stage and location are shown in Table 1 and Additional file 1. Twenty-three patients with early onset CRC were enrolled from 4 different hospitals in the south-eastern region of Norway. HNPCC, FAP and other known syndromes were excluded after a thorough family and medical history. Seventeen patients with late onset CRC, treated in one of the hospitals, were selected after matching for sex, tumor location and stage.

Table 1 Summary of patient data in the early and late onset CRC series.

Tissue samples were taken and preserved in the operation theatre immediately after resection of the specimen. In the early onset group samples were transferred to tubes with RNA-later RNA Stabilization Reagent (Qiagen, Hilden, Germany), stored at room temperature over night, then transferred to an empty tube and frozen at -80°C for long-term storage. In the late onset group samples were frozen in liquid nitrogen and stored at -80°C until use.

Written informed consent was obtained from all subjects included. The research biobanks are registered according to national legislation and the research project is approved by the Regional Committee for Medical Research Ethics (REK South-East: 1.2005.1629; REK South: 2003, S-02126).

Nucleic acids isolation

The tumor tissue was manually ground in liquid N2 using a mortar and pestle before isolation of DNA. DNA was extracted using a semi-automatic phenol-chloroform extraction method followed by ethanol precipitation in a 340A Nucleic Acid Extractor (Applied Biosystems, Foster City, CA, USA). Ensuring good DNA quality, DNA was measured by using NanoDrop ND-1000 (Thermo Fisher Scientific, Waltham, MA, USA). Samples with OD 260 nm/280 nm >1.7 were included. Thereafter, the DNA was diluted in TE-buffer to a final concentration of 250 ng/μl. The integrity of the DNA was visually inspected on a 1% agarose gel.

RNA was isolated using AllPrep DNA/RNA mini kit (Qiagen). The RNA quality and quantity was measured using the NanoDrop ND-1000 and OD 260 nm/280 nm and OD 260 nm/230 nm was carefully evaluated. Degradation was measured using BioAnalyzer (Agilent Technologies, Santa Clara, CA, USA). All samples included had RNA concentrations above 150 ng/μl, and RNA integrity number (RIN) above 8.0.

Microsatellite instability analysis

MSI status was determined for all tumors and positive samples were excluded from further analyses, in order to avoid any potential undetected HNPCC patient. Determination of MSI status was performed as described by Wu et al.[30], using the Bethesda markers. High degree of microsatellite instability (MSI-H) was defined if two or more markers showed aberrant profile after fragment analysis.

Array Comparative Genomic Hybridization (aCGH)

From each sample, 1 μg DNA was included in the analysis alongside with gender-matched human reference DNA. Each of the two reference samples, one female and one male reference, consisted of a pool of DNA from normal lymphocytes from four healthy persons. The applied microarray platform, Roche NimbleGen, Human Whole-Genome Array CGH Analysis v1 (Roche Diagnostics, Mannheim, Germany), provides measurements from 385 000 unique genomic loci. The aCGH experimental procedure, as well as raw data pre-processing and normalization, was performed by NimbleGen Systems Inc. at their facility in Iceland [31, 32]. For mapping of genomic breakpoints the segMNT v1.1 CGH segmentation analysis algorithm were run in NimbleScan™ software v2.4 (Roche, Basel, Switzerland). A window-size of approximately 60 000 base pairs, calculating the average of 10 probes per window, was further used for analysis.

The female reference pool was hybridized against the male reference pool as a technical control. To determine the threshold for scoring of gain and loss, normalized, log2-transformed ratios were used. Based on the variation in autosomal genomic regions, which should not vary between the two reference samples, thresholds for averaged log2 ratio data were set to 0.1 and -0.1 for gains and losses, respectively. For the control experiment this implied a false positive rate per genomic loci of .005 and .007 for gains and losses, respectively. As expected when hybridizing female DNA vs. male DNA, gain was seen for the whole of chromosome × whereas loss was seen for the whole of chromosome Y. Furthermore, two samples were run twice on two different arrays. For each doublet of samples the squared correlation coefficient (R2) was calculated from all measurements if one or both samples differed from normal copy number. The two samples showed correlation coefficient of 0.84 and 0.75, respectively. For pairs of random samples the correlation coefficient was calculated in the same manner. Median correlation coefficient for ten pairs of random samples was 0.16 (range 0.04 - 0.32).

Statistical differences in frequencies of gains and losses, between the groups of early onset and late onset CRC patients, were calculated using Fisher's exact test. Significance values were for each window/segment calculated based on number of patients with gain or loss compared to no-gain and no-loss, respectively. Affected telomere regions and regions spanning centromeres were not considered as these regions have highly repetitive fragments and may reflect changes that are technical in nature. Also, changes in sex chromosomes were not considered. The number of chromosomal aberrations and the number of breakpoints were calculated as follows: A chromosomal aberration spanning an entire chromosome was recorded as one chromosomal aberration and zero breakpoints; a chromosomal aberration inside a chromosome was recorded as one chromosomal aberration and two breakpoints.

To assign the corresponding chromosome band to specified breakpoints, and identify genes included in given regions we used the UCSC Genome Browser (Mar. 2006 (hg18) assembly, URL: http://www.genome.ucsc.edu) and BioMart genome annotation tool (version 0.7, Ensemble version 53, URL: http://www.biomart.org/), respectively.

Integration of DNA copy numbers to gene expression data

Genome-wide measurements of mRNA levels were obtained by AB1700 gene expression microarrays (Applied Biosystems, Foster City, CA, USA). Information on sample handling and preprocessing of the raw data will be published elsewhere (Ã…gesen et al., manuscript in preparation). For the integration analysis, we used quantile normalized data with gene-wise centering on the median of the dataset.

Each measurement on the mRNA expression array was related to the DNA copy number status of the corresponding genomic loci. A visual basic (VBA) script was written to automate the analyses and to perform integrated statistics genome-wide on a gene-by-gene and locus-by-locus manner. For each gene, and for gain, loss and normal copy number independently, the median of the expression levels for the samples with the DNA copy number changes was divided by the median of the expression levels for the samples with normal copy number. To indicate genes with concomitant gain and over-expression, we set a threshold at 1.5, and genes with concomitant loss and under-expression were indicated when the value was below 0.75. Furthermore, the expression levels for genes located in regions with statistically significant difference in gain or loss in the early onset compared to the late onset CRC samples, were evaluated using Independent samples T-test in SPSS. Annotations of genes are used according to AB1700 annotation file (version; 20060930_ab1700_human, Applied Biosystems). Overview of the study is given in Figure 1.

Figure 1
figure 1

Study design comparing the tumor genomes of early onset and late onset colorectal cancer patients. Tumor samples from early onset and late onset CRC patients were analyzed with aCGH (DNA/genome level) and the data were integrated with the mRNA expression (RNA/transcriptome level) of the same samples. About 10 000 genes show corresponding DNA and RNA level. Exclusion of genes located outside chromosomal regions with statistically significant difference between the two patient groups resulted in 107 genes, summarized in Additional file 6. Further, exclusion of genes whose mRNA expression levels were not statistically significant between patients groups resulted in a short list of seven potential predisposing genes, as summarized in Table 3.

TP53 mutation analysis

The total protein coding region of TP53 was amplified in a multiplex PCR reaction containing five distinct PCR fragments by using flanking intronic primers with M13 tails. Multiplex PCR kit was used as recommended by the vendor (Qiagen, Hilden, Germany). If initial multiplex PCR reaction did not succeed, the same fragments were amplified separately using each of the primer pairs. Primer and fragment details are described in Additional file 2. After visual inspection of quality and quantity on a polyacrylamide gel, the product was purified using Exo-SAP-IT (USB Corporation, Cleveland, OH, USA). Sequencing of the purified products was performed in both 5'- and 3' directions using BigDye Terminator v3.1 kit (Applied BioSystems, Foster City, CA, USA). Primers for the sequencing reaction were identical with those used in the initial PCR for fragments amplified in multiplex reactions. M13-primers were used in sequencing reactions of singleplex PCR products. The resulting sequence product was further purified using multiscreen plates, (MilliPore, Billerica, MA, USA) with Sephadexâ„¢ G-50 Superfine (GE Healthcare, Chalfont St. Giles, UK), and subjected to sequencing at a 3730 DNA Analyzer (Applied Biosystems, Foster City, CA, USA). If a mutation was detected, a new independent PCR product was subjected to sequencing to confirm our finding.

Results

Genomic profile of colorectal carcinomas

The profile of DNA copy number gains and losses across all samples is shown in Figure 2A. Overall, across the 40 MSS colorectal samples the median number of chromosomal aberrations was 31.5 (range 11-93), and the median percentage of probes showing normal copy number was 78% (range 48%-99%).

Figure 2
figure 2

Genomic profile across all chromosomes from 40 tumor samples. A) Copy number profile across all 40 colorectal carcinomas, including early onset (n = 23) and late onset (n = 17) tumors. Percentage of samples with gains is shown in red and losses in green. Gains and losses in the short arm of acrocentric chromosomes are not considered and further discussed. The centromere positions are indicated by dashed lines, and changes recorded at the p-arms and q-arms are shown to the left and right of the dashed lines, respectively. B) Copy number profiles of early onset (solid) as compared to late onset (line) colorectal carcinomas.

Gains at 8q, 13q and 20q and losses from 4p, 4q, 8p, 17p, 18p and 18q were most frequent (range 43%-78%). Aberrations in the p-arm of acrocentric chromosomes were not considered as these most probably reflect technical challenges. The most frequent loss, 18q22.1-22.3 affecting ~10 Mbps, was found in 78% of all samples. This region comprises 34 protein coding genes. The region most often gained, found in 70% of the tumors, were a ~7.5 Mbps region at 20q13.31-13.33. This region contains 10 miRNAs and 100 genes. One of the most commonly affected chromosomes, chromosome 8, showed losses from the p-arm in 65% of all tumors and gains at the q-arm in 43%. For the majority of the samples, the breakpoint between gain and loss at chromosome 8 was located in the p-arm, about 5 Mbps away from the centromere, as shown in Additional file 3.

Genomic profiles of early onset versus late onset colorectal carcinomas

The average fraction of the genome with aberrant copy number was similar in the early onset and late onset CRC series, 26% versus 22%, respectively (P = 0.53). However, both the number of chromosomal aberrations and thus the number of breakpoints differed significantly between the groups (P = 0.02 and P = 0.01, respectively). Chromosome regions with aberrations are listed in Table 2. Overview of aberrations from each patient sample is shown in Additional file 4.

Table 2 Genomic loci with statistically significant copy number variation between early and late onset colorectal cancers.

A plot of gains and losses along a genome axis for each of the two patient groups shows distinct differences (Figure 2B). Statistical calculations suggested that 3.2% of the genome was significantly different between the two groups. This includes 16 different regions at 6 unique chromosomes, spans a total of 58.3 Mbps, and affects 107 genes, (Figure 3). Four regions with statistical significant difference were spanning centromere regions and were excluded from further analysis. Only 0.3 Mbps of these changes (one region) were more often observed in late onset cancers than in the early onset group (Table 2). The parts of the genome that exhibited significant difference between the two tumor groups contains, among others, 574 protein coding genes, 44 miRNA, and 47 pseudogenes, as listed in Additional file 5.

Figure 3
figure 3

Genes mapping to the chromosome bands with copy number aberrations associated with early onset CRC. Chromosomal sites gained or lost with statistical significance in the early compared to the late onset tumor group are shown in red or green boxes, respectively. Chromosomes with no aberrations or aberrations in centromere regions only, are shaded. Genes which have expression levels corresponding to chromosomal aberration, and concurrently located within the statistically significant regions are indicated within the boxes. Genes which also have significantly different expression level between early onset and late onset groups are further marked with green or red background.

Overall, the smallest region of overlap (SRO) of losses with the most statistically significant difference between early onset and late onset patients (P = 0.003) was located 0.24 Mbps adjacent to the centromere in chromosome 4. Only one miRNA maps to this position. Fourteen of the 23 (61%) early onset CRC patients, as compared to two of 17 (12%) late onset patients, displayed loss in this region. The most significant SRO of losses containing protein coding genes was located at 1q21.1 (P = 0.01). This position contains one protein coding gene, namely PRKAB2, encoding a regulatory subunit of the AMP-activated protein kinase (AMPK). Ten of 23 (43%) early onset CRC patients, and only one of 17 (6%) late onset patients, displayed loss of this locus.

Copy number differences between the patient groups were found for large parts of chromosome 4 (4p16.1-q28.3) (Figure 2B). In total, the regions on chromosome 4 with significant differences between the groups contained 258 protein coding genes, 36 miRNAs and 27 pseudogenes.

Two sequence stretches within 10q22.3 and 19q13.31 were gained in 8 of 23 patients in the early onset group versus none among the 17 late onset patients (P = 0.01). The 10q region contains four protein coding genes, SFTPA2, ZCCHC24, PPIF and ZMIZ1, whereas the 19q region includes only PSG3 (pregnancy-specific beta-1-glycoprotein 3). At chromosome 19 a 4 Mbps gained region (19q13.2-q13.31) was significantly more frequent in the early onset group. This comprises 103 protein coding genes and three miRNAs. The protein coding genes include, among others, a group of genes in the carcinoembryonic antigen (CEA) family.

Only one region, spanning 0.3 Mbps at 5q13.2, was gained more often in late onset patients (5 of 17) compared to early onset patients (none of 23). The region contains 4 putative protein coding genes, β-clucoronidase-like protein SMA5, β-clucoronidase-like protein SMA4, ENSG00000205565, and ENSG00000197370.

Integrated genome and transcriptome analysis

An integration analysis of the DNA copy number data and mRNA expression data revealed that 37% (~10 900) of the genes had corresponding DNA copy number and RNA expression levels. Among these, 107 genes were located in genomic regions with significantly different frequency of copy number changes between the two tumor groups (Additional file 6). When sorting these 107 genes based on fold-change of expression levels between the groups, CLC, CEACAM7, FCGBP, hGC1641029.5 (Probe ID:218092) and CD79A were top five for gains, and the top five list for losses included LOC391661(Probe ID:152373), ARSJ, ADH6, SPP1 and BANK1. From the 107 genes, those with significantly different mRNA expression levels between the early onset and late onset groups were identified. The resulting strict gene list included seven genes (Table 3 and Figure 3). Of these, three genes (CLC, LTBP4 and ZNF574), all within chromosome band 19q13.2, were concomitantly gained and up-regulated, and four genes (PPAT, EIF4E, RG9MTD2 and PLA2G12A) were lost and down-regulated. The PPAT gene is located in 4q21, EIF4E and RG9MTD2 in 4q23 and PLA2G12A in 4q25. The expression levels for each gene are shown in Figure 4.

Table 3 Genes/Probe IDs showing statistically significant different DNA copy number and mRNA expression in early versus late onset colorectal cancers.
Figure 4
figure 4

mRNA expression levels of a short list of seven candidate genes for early onset CRC. Log transformed mRNA expression values plotted against age. Red bullets indicate copy number gain, green bullets indicate loss and grey bullets indicate normal copy number. Box plots for each of the patient groups are indicated for each gene.

TP53 mutations and genetic complexity

Mutation status of TP53 was evaluated for all protein coding exons. Sixty-five percent of the samples were mutated, with no significant difference being observed in the early onset group compared to the late onset group (Additional file 1). Five percent of the mutations were located outside the DNA binding domain of TP53. A region frequently lost at 17p comprises the p13.1 band where the TP53 gene is encoded. Loss in this region was observed in 55% of the tumor samples, also with no frequency difference between the early onset and the late onset group. Overall, 80% of the samples displayed mutation in TP53 and/or DNA copy number aberration of 17p13.1. Both loss of 17p13.1 and mutation in TP53 were independently associated with an increased number of chromosomal aberrations (P < 0.01 and P < 0.01, respectively, and P = 0.03 for loss and mutation combined). This difference was still significant when comparing mutated (n = 9) and wild type (n = 8) samples with normal copy number in 17p13.1 only (P = 0.01).

The concurrent loss and reduced expression of PPAT, EIF4E, RG9MTD2, and PLA2G12A, were significantly correlated to TP53 mutations (P = 0.006) in the tumors from the early onset patients. The same trend was seen for tumors with concurrent gain and elevated expression of CLC, LTBP4 and ZNF574, though not statistically significant, (P = 0.17).

Discussion

In Norway, around 5% of all patients with CRC are diagnosed before the age of 55. However, the incidence rate increases with increasing age, and the median age at diagnosis is 70 years. Individuals with early onset of CRC may be carriers of gene variants causing an increased risk for disease, but most of these predisposition genes remain to be identified. The present study demonstrates genetic aberrations that are associated with early onset CRC. Compared to the current literature, the number of patients (n = 40) is fairly high, and the number of measurements provided is very high. To our knowledge no previous study has at this resolution level compared the genomic changes between tumors from patients with ~20 year difference in onset. Methods for genome-wide assessment of DNA copy number alterations in most studies published to date are generally of low to medium resolution. The introduction of aCGH has increased the resolution hundreds of times, and with the microarray platform used in this study the total number of probes is 385 000, giving an average probe interval of 6000 bases. In addition, measurements of mRNA expression levels from the same samples have allowed for a short list of candidate genes.

Genomic alterations in colorectal cancers

The overall DNA copy number profile for CRC found in the present study is in line with previous reports [19], suggesting that the series used are representative for microsatellite stable CRCs. The profile includes frequent gains at 8q, 13q, and 20q and losses from 4p, 4q, 8p, 17p, 18p and 18q.

As expected the most common loss and gain were found at 18q and 20q, respectively. The SROs in these chromosome arms contains a number of protein coding genes, including well known critical cancer genes as BCL2 (18q22), DCC (18q21) and AURKA (20q13). The latter is known to be important for normal chromosome segregation, observed as over-expressed in cancers with chromosomal instability (CIN) [33]. Also present in 20q13 are several members of the SERPIN gene family, including SERPINB5, positively regulated by TP53 [34]. Furthermore, one of the 10 miRNAs in this region is the MIR646 (hsa-miR-646), reported to be related to colorectal cancer (The miR-ontology database, http://ferrolab.dmi.unict.it/miro) [35].

In CRC, aberrations at chromosome 8 are commonly reported as loss in 8p and gain in 8q [36]. Our data shows that the breakpoint is located on the p-arm, approximately five Mbps from the centromere (Additional file 3). This is in accordance with findings from others [37, 38]. This clearly shows that the resolution of aCGH is crucial for a more exact evaluation of breakpoint regions. Loss of gross parts of both arms of chromosome 4 is frequently reported in several cancers [39–44], and the location of tumor suppressor genes important for tumorigenesis in this region is plausible [45]. In this study we observed on average 40% of losses in chromosome 4.

Differences in genome profile between early onset and late onset CRC

In the early onset sample series, the number of chromosomal aberrations and the number of breakpoints were significantly higher than in the late onset sample series. However, the percentage of measurements showing aberrations does not differ significantly. This indicates that about the same amount of DNA is affected by copy number changes in the two series, but the number of events is larger in the early onset group compared to the late onset samples, which again suggests increased genome instability in the early onset tumors.

Fifteen loci at five unique chromosomes were statistically more often altered in the tumors of the early onset group compared to the late onset group (Table 2). The gained region at chromosome 19 with statistically significant difference between the two patient groups spanned 4 Mbps. This region comprises more than one hundred protein coding genes, among them a group of genes in the carcinoembryonic antigen (CEA) family. Genes in this family encode pregnancy-specific beta-I-glycoproteins and carcinoembryonic antigen-related cell adhesion molecules. Also, 19q13.2-q13.31 contains three microRNAs, of which one, the MIR641 (hsa-miR-641), is known to be associated with colorectal cancer [35]. Chromosome 19 is notable for the highest gene density of all human chromosomes, large clustered gene families and high GC content, indicating biological and evolutionary significance [46]. In a recent GWAS study, a region on 19q13.1 was identified as a putative susceptibility locus for CRC [10]. This region is 6.4 Mbps proximal to the region we have identified at 19q13.2-q13.31.

In the present study six smaller regions at chromosome 4, comprising ~70 Mbps in total, were identified with statistical significant losses in the early onset group compared to the late onset patients (Table 2). In total ~250 genes and 30 miRNAs are located in the affected regions. Even though deletions in chromosome 4 is a common feature in CRC, and known cancer genes such as KIT, EGF and FGF2 are situated here, none has yet been verified as predisposing for early onset or hereditary CRC.

Chromosomal loci previously suggested to contain genes that may predispose to CRC [4–10] are not the same as the ones found in the present study. This may best be explained by the difference in inclusion criteria; whereas the GWAS-studies have mainly been performed in samples from individuals from high-risk families, the present study selected patients with young age at diagnosis of disease. Mourra et al. reported a statistically higher frequency of deleted 14q loci in patients <50 years at primary diagnosis versus patients >50 years at diagnosis [11], but this could not be confirmed in the present study.

Integration of DNA copy number and mRNA expression data

The integration analysis of genome copy number data with mRNA expression data identified a short list of 7 target genes differing between the early- and late onset groups. The three gained and up-regulated genes, CLC, LTBP4, and ZNF574, were altered in the same tumors. The CLC gene differentiates the tumor groups the most with regard to mRNA expression levels (P = 0.001). The CLC protein (Charcot-Leyden crystal protein/Galectin-10) is a lysophospholipase in the galectin super family of proteins, normally expressed in eosinophils and basophils, associated with inflammation and some myeloid leukemia [47]. However, data on CLC in cancer is sparse.

The LTBP4 protein (latent TGF-β binding protein 4) has structural roles in the extracellular matrix as well as its participation in the TGFβ-pathway. TGFβ-pathway components are commonly altered in cancer in general, and in colorectal cancer in particular. On of them, SMAD4, cause juvenile polyposis when mutated in germline cells. It is localized at chromosome 18q21 together with SMAD2, in a region frequently deleted in colorectal cancers. TGFβRII mutations are found in ~30% of colorectal cancers and are considered the main mechanism of inhibition of the TGFβ-pathway. Although mutations in TGFβRI are rare, an association between the TGFBR1*6A-polymorphism and cancer has been reported [48]. Together with LAPs (latency-associated proteins), LTBP4 it is involved in the assembly, secretion and targeting of TGFβs [49]. TGF-β has an important, complex and somewhat dual role in the normal cell. Downstream effects are as diverse as promotion of migration, adhesion and differentiation on one hand, to inhibition of growth, cell-cycle control and apoptosis on the other [50]. As opposed to the present results, others have reported a suppression of LTBP4 to cause cancer [51, 52]. An increase in LTBP4 expression might potentially cause colorectal cancer development. Elevated levels might increase deposition of TGFβ in extracellular matrix, thus hindering TGFβ to bind to the receptor, and thereby prevent tumor suppressive downstream reactions such as apoptosis, and growth arrest.

The concomitantly lost and down-regulated genes PPAT, EIF4E, RG9MTD2, and PLA2G12A are located on different bands of chromosome 4. All these four candidate genes and potential tumor suppressor genes had a combined genomic loss and mRNA under-expression in the same patients, presumably associated with monosomy. PLA2G12A is a member of the secreted PLA2s (sPLA2) in the phospholipase A2 family of proteins, found to be involved in tumorigenesis [53–55]. The enzymatic activity of PLA2s is to hydrolyze the fatty acid from membrane phospholipids, which are further metabolized and forms eicosanoids and bioactive lipid mediators [56]. PLA2G12A has been reported to be highly expressed in normal and tumor tissues from the colon [53], which suggests that reduced expression, as reported in this study, might contribute to an oncogenic transformation in a subset of early onset CRCs. The protein EIF4E is a translation initiation factor, and mRNA expression is reported to be elevated in many human cancers [57]. mRNA levels are reported to increase during tumor formation and progression of colorectal cancers [58]. In our study the expression of EIF4E was found to be reduced in the early onset group compared to the late onset group, which may indicate that early onset cancers develop in an EIF4E independent manner.

TP53 mutations and genomic complexity

TP53 mutations are associated with CIN tumors. The TP53 gene was sequenced to confirm an equal representation of mutations in the two age groups, to ensure that the series were unbiased. The seemingly high frequency of 65% mutations reported here is most likely due to the fact that all protein coding exons in TP53 were analyzed, and that only microsatellite stable tumors are included in the study. The samples with loss in 17p13.1 had significantly more aberrations throughout the whole genome, compared to the samples with normal copy number or gain. This is as expected, partly because loss of 17p is associated with CIN phenotype, and partly because one expects any sample with aberration in a given position to have more aberrations globally. Furthermore, a statistically significant difference in percentage of aberrations was seen when comparing mutated (n = 27) to wild type (n = 13) samples, irrespective of DNA copy number in the region where TP53 is located. Samples with no aberrations in TP53 showed a lower degree of genetic complexity measured as the copy number level throughout the genome. This clearly indicates that non-functional TP53, regardless of how it is inactivated, plays an important role in tumorigenesis associated with CIN.

Interestingly, loss of 1q21 containing the AMP kinase subunit β2 (PRKAB2), which was recently shown by the Arnold Levine laboratory to be regulated by TP53 [59], are more often lost in tumors with wild type TP53 than in those with mutated TP53 (P = 0.08).

Cancer susceptibility and genetic pathways

Our finding that carcinomas from early onset patients present a high number of chromosomal aberrations and breakpoints is comparable with the concept that they have inherited a germline mutation in a gene relevant for chromosome stability, with loss of the second allele occurring somatically. Relevant examples for such a relationship are the complex pattern of copy number changes observed in breast carcinomas associated with germline BRCA1 mutations (when compared with sporadic breast carcinomas), and the microsatellite instability observed in colorectal carcinomas from Lynch syndrome patients [60–62]. The data we present here highlight possible candidate susceptibility loci that complement other studies using different strategies [5, 6, 10, 21–25, 63–66]. Alternatively, the changes we have found more frequent in tumors from early onset patients, may reflect changes as a result of somatic mutations related to a yet uncovered germline mutation.

Conclusions

This is the first study that has identified distinct tumor genome profiles in early onset and late onset CRC patients at 6 Kbps resolution level, combined with the corresponding mRNA expression profile. Sixteen genomic loci containing more than 500 coding genes were identified as preferentially altered in cancers from early onset patients when compared to late onset patients. Finally, by integration analysis with gene expression data from the same samples, we identified a short list of seven candidate genes as potentially predisposing to early onset CRC. Further studies are warranted to find out if early onset CRC is caused by inherited low penetrance alleles or germline mutations in these or other genes.

References

  1. de la Chapella A: Genetic predisposition to colorectal cancer. Nat Rev Cancer. 2004, 4: 769-780. 10.1038/nrc1453

    Article  Google Scholar 

  2. Rustgi AK: The genetics of hereditary colon cancer. Genes Dev. 2007, 21: 2525-2538. 10.1101/gad.1593107

    Article  CAS  PubMed  Google Scholar 

  3. Lichtenstein P, Holm NV, Verkasalo PK, Iliadou A, Kaprio J, Koskenvuo M, Pukkala E, Skytthe A, Hemminki K: Environmental and heritable factors in the causation of cancer--analyses of cohorts of twins from Sweden, Denmark, and Finland. N Engl J Med. 2000, 343: 78-85. 10.1056/NEJM200007133430201

    Article  CAS  PubMed  Google Scholar 

  4. Jaeger E, Webb E, Howarth K, Carvajal-Carmona L, Rowan A, Broderick P, Walther A, Spain S, Pittman A, Kemp Z, Sullivan K, Heinimann K, Lubbe S, Domingo E, Barclay E, Martin L, Gorman M, Chandler I, Vijayakrishnan J, Wood W, Papaemmanuil E, Penegar S, Qureshi M, Farrington S, Tenesa A, Cazier JB, Kerr D, Gray R, Peto J, Dunlop M: Common genetic variants at the CRAC1 (HMPS) locus on chromosome 15q13.3 influence colorectal cancer risk. Nat Genet. 2008, 40: 26-28. 10.1038/ng.2007.41

    Article  CAS  PubMed  Google Scholar 

  5. Zanke BW, Greenwood CM, Rangrej J, Kustra R, Tenesa A, Farrington SM, Prendergast J, Olschwang S, Chiang T, Crowdy E, Ferretti V, Laflamme P, Sundararajan S, Roumy S, Olivier JF, Robidoux F, Sladek R, Montpetit A, Campbell P, Bezieau S, O'Shea AM, Zogopoulos G, Cotterchio M, Newcomb P, McLaughlin J, Younghusband B, Green R, Green J, Porteous ME, Campbell H: Genome-wide association scan identifies a colorectal cancer susceptibility locus on chromosome 8q24. Nat Genet. 2007, 39: 989-994. 10.1038/ng2089

    Article  CAS  PubMed  Google Scholar 

  6. Tomlinson I, Webb E, Carvajal-Carmona L, Broderick P, Kemp Z, Spain S, Penegar S, Chandler I, Gorman M, Wood W, Barclay E, Lubbe S, Martin L, Sellick G, Jaeger E, Hubner R, Wild R, Rowan A, Fielding S, Howarth K, Silver A, Atkin W, Muir K, Logan R, Kerr D, Johnstone E, Sieber O, Gray R, Thomas H, Peto J: A genome-wide association scan of tag SNPs identifies a susceptibility variant for colorectal cancer at 8q24.21. Nat Genet. 2007, 39: 984-988. 10.1038/ng2085

    Article  CAS  PubMed  Google Scholar 

  7. Tomlinson IP, Webb E, Carvajal-Carmona L, Broderick P, Howarth K, Pittman AM, Spain S, Lubbe S, Walther A, Sullivan K, Jaeger E, Fielding S, Rowan A, Vijayakrishnan J, Domingo E, Chandler I, Kemp Z, Qureshi M, Farrington SM, Tenesa A, Prendergast JG, Barnetson RA, Penegar S, Barclay E, Wood W, Martin L, Gorman M, Thomas H, Peto J, Bishop DT: A genome-wide association study identifies colorectal cancer susceptibility loci on chromosomes 10p14 and 8q23.3. Nat Genet. 2008, 40: 623-630. 10.1038/ng.111

    Article  CAS  PubMed  Google Scholar 

  8. Broderick P, Carvajal-Carmona L, Pittman AM, Webb E, Howarth K, Rowan A, Lubbe S, Spain S, Sullivan K, Fielding S, Jaeger E, Vijayakrishnan J, Kemp Z, Gorman M, Chandler I, Papaemmanuil E, Penegar S, Wood W, Sellick G, Qureshi M, Teixeira A, Domingo E, Barclay E, Martin L, Sieber O, Kerr D, Gray R, Peto J, Cazier JB, Tomlinson I: A genome-wide association study shows that common alleles of SMAD7 influence colorectal cancer risk. Nat Genet. 2007, 39: 1315-1317. 10.1038/ng.2007.18

    Article  CAS  PubMed  Google Scholar 

  9. Tenesa A, Farrington SM, Prendergast JG, Porteous ME, Walker M, Haq N, Barnetson RA, Theodoratou E, Cetnarskyj R, Cartwright N, Semple C, Clark AJ, Reid FJ, Smith LA, Kavoussanakis K, Koessler T, Pharoah PD, Buch S, Schafmayer C, Tepel J, Schreiber S, Volzke H, Schmidt CO, Hampe J, Chang-Claude J, Hoffmeister M, Brenner H, Wilkening S, Canzian F, Capella G: Genome-wide association scan identifies a colorectal cancer susceptibility locus on 11q23 and replicates risk loci at 8q24 and 18q21. Nat Genet. 2008, 40: 631-637. 10.1038/ng.133

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Houlston RS, Webb E, Broderick P, Pittman AM, Di Bernardo MC, Lubbe S, Chandler I, Vijayakrishnan J, Sullivan K, Penegar S, Carvajal-Carmona L, Howarth K, Jaeger E, Spain SL, Walther A, Barclay E, Martin L, Gorman M, Domingo E, Teixeira AS, Kerr D, Cazier JB, Niittymaki I, Tuupanen S, Karhu A, Aaltonen LA, Tomlinson IP, Farrington SM, Tenesa A, Prendergast JG: Meta-analysis of genome-wide association data identifies four new susceptibility loci for colorectal cancer. Nat Genet. 2008, 40: 1426-1435. 10.1038/ng.262

    Article  CAS  PubMed  Google Scholar 

  11. Mourra N, Zeitoun G, Buecher B, Finetti P, Lagarde A, Adelaide J, Birnbaum D, Thomas G, Olschwang S: High frequency of chromosome 14 deletion in early-onset colon cancer. Dis Colon Rectum. 2007, 50: 1881-1886. 10.1007/s10350-007-9040-3

    Article  PubMed  Google Scholar 

  12. Grady WM, Carethers JM: Genomic and epigenetic instability in colorectal cancer pathogenesis. Gastroenterology. 2008, 135: 1079-1099. 10.1053/j.gastro.2008.07.076

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Lengauer C, Kinzler KW, Vogelstein B: Genetic instability in colorectal cancers. Nature. 1997, 386: 623-627. 10.1038/386623a0

    Article  CAS  PubMed  Google Scholar 

  14. Lothe RA: Microsatellite instability in human solid tumors. Mol Med Today. 1997, 3: 61-68. 10.1016/S1357-4310(96)10055-1

    Article  CAS  PubMed  Google Scholar 

  15. Fearon ER, Vogelstein B: A genetic model for colorectal tumorigenesis. Cell. 1990, 61: 759-767. 10.1016/0092-8674(90)90186-I

    Article  CAS  PubMed  Google Scholar 

  16. Hoglund M, Gisselsson D, Hansen GB, Sall T, Mitelman F, Nilbert M: Dissecting karyotypic patterns in colorectal tumors: two distinct but overlapping pathways in the adenoma-carcinoma transition. Cancer Res. 2002, 62: 5939-5946.

    CAS  PubMed  Google Scholar 

  17. Carvalho B, Postma C, Mongera S, Hopmans E, Diskin S, Wiel van de MA, Van CW, Thas O, Matthai A, Cuesta MA, Terhaar Sive Droste JS, Craanen M, Schrock E, Ylstra B, Meijer GA: Multiple putative oncogenes at the chromosome 20q amplicon contribute to colorectal adenoma to carcinoma progression. Gut. 2009, 58: 79-89. 10.1136/gut.2007.143065

    Article  CAS  PubMed  Google Scholar 

  18. Grade M, Becker H, Liersch T, Ried T, Ghadimi BM: Molecular cytogenetics: genomic imbalances in colorectal cancer and their clinical impact. Cell Oncol. 2006, 28: 71-84.

    CAS  PubMed  Google Scholar 

  19. Diep CB, Kleivi K, Ribeiro FR, Teixeira MR, Lindgjaerde OC, Lothe RA: The order of genetic events associated with colorectal cancer progression inferred from meta-analysis of copy number changes. Genes Chromosomes Cancer. 2006, 45: 31-41. 10.1002/gcc.20261

    Article  CAS  PubMed  Google Scholar 

  20. Hemminki A, Tomlinson I, Markie D, Jarvinen H, Sistonen P, Bjorkqvist AM, Knuutila S, Salovaara R, Bodmer W, Shibata D, de la CA, Aaltonen LA: Localization of a susceptibility locus for Peutz-Jeghers syndrome to 19p using comparative genomic hybridization and targeted linkage analysis. Nat Genet. 1997, 15: 87-90. 10.1038/ng0197-87

    Article  CAS  PubMed  Google Scholar 

  21. Poynter JN, Figueiredo JC, Conti DV, Kennedy K, Gallinger S, Siegmund KD, Casey G, Thibodeau SN, Jenkins MA, Hopper JL, Byrnes GB, Baron JA, Goode EL, Tiirikainen M, Lindor N, Grove J, Newcomb P, Jass J, Young J, Potter JD, Haile RW, Duggan DJ, Le ML: Variants on 9p24 and 8q24 are associated with risk of colorectal cancer: results from the Colon Cancer Family Registry. Cancer Res. 2007, 67: 11128-11132. 10.1158/0008-5472.CAN-07-3239

    Article  CAS  PubMed  Google Scholar 

  22. Valle L, Serena-Acedo T, Liyanarachchi S, Hampel H, Comeras I, Li Z, Zeng Q, Zhang HT, Pennison MJ, Sadim M, Pasche B, Tanner SM, de la Chapelle A: Germline allele-specific expression of TGFBR1 confers an increased risk of colorectal cancer. Science. 2008, 321: 1361-1365. 10.1126/science.1159397

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Wiesner GL, Daley D, Lewis S, Ticknor C, Platzer P, Lutterbaugh J, MacMillen M, Baliner B, Willis J, Elston RC, Markowitz SD: A subset of familial colorectal neoplasia kindreds linked to chromosome 9q22.2-31.2. Proc Natl Acad Sci USA. 2003, 100: 12961-12965. 10.1073/pnas.2132286100

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Neklason DW, Kerber RA, Nilson DB, nton-Culver H, Schwartz AG, Griffin CA, Lowery JT, Schildkraut JM, Evans JP, Tomlinson GE, Strong LC, Miller AR, Stopfer JE, Finkelstein DM, Nadkarni PM, Kasten CH, Mineau GP, Burt RW: Common familial colorectal cancer linked to chromosome 7q31: a genome-wide analysis. Cancer Res. 2008, 68: 8993-8997. 10.1158/0008-5472.CAN-08-1376

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Daley D, Lewis S, Platzer P, MacMillen M, Willis J, Elston RC, Markowitz SD, Wiesner GL: Identification of susceptibility genes for cancer in a genome-wide scan: results from the colon neoplasia sibling study. Am J Hum Genet. 2008, 82: 723-736. 10.1016/j.ajhg.2008.01.007

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Tenesa A, Dunlop MG: New insights into the aetiology of colorectal cancer from genome-wide association studies. Nat Rev Genet. 2009, 10: 353-358. 10.1038/nrg2574

    Article  CAS  PubMed  Google Scholar 

  27. Hong Y, Ho KS, Eu KW, Cheah PY: A susceptibility gene set for early onset colorectal cancer that integrates diverse signaling pathways: implication for tumorigenesis. Clin Cancer Res. 2007, 13: 1107-1114. 10.1158/1078-0432.CCR-06-1633

    Article  CAS  PubMed  Google Scholar 

  28. Tomlinson IP, Dunlop M, Campbell H, Zanke B, Gallinger S, Hudson T, Koessler T, Pharoah PD, Niittymakix I, Tuupanenx S, Aaltonen LA, Hemminki K, Lindblom A, Forsti A, Sieber O, Lipton L, van WT, Morreau H, Wijnen JT, Devilee P, Matsuda K, Nakamura Y, Castellvi-Bel S, Ruiz-Ponte C, Castells A, Carracedo A, Ho JW, Sham P, Hofstra RM, Vodicka P: COGENT (COlorectal cancer GENeTics): an international consortium to study the role of polymorphic variation on the risk of colorectal cancer. Br J Cancer. 2009, 102: 447-454. 10.1038/sj.bjc.6605338

    Article  PubMed Central  PubMed  Google Scholar 

  29. Pinkel D, Segraves R, Sudar D, Clark S, Poole I, Kowbel D, Collins C, Kuo WL, Chen C, Zhai Y, Dairkee SH, Ljung BM, Gray JW, Albertson DG: High resolution analysis of DNA copy number variation using comparative genomic hybridization to microarrays. Nat Genet. 1998, 20: 207-211. 10.1038/2524

    Article  CAS  PubMed  Google Scholar 

  30. Wu Q, Lothe RA, Ahlquist T, Silins I, Trope CG, Micci F, Nesland JM, Suo Z, Lind GE: DNA methylation profiling of ovarian carcinomas and their in vitro models identifies HOXA9, HOXB5, SCGB3A1, and CRABP1 as novel targets. Mol Cancer. 2007, 6: 45-55. 10.1186/1476-4598-6-45

    Article  PubMed Central  PubMed  Google Scholar 

  31. Selzer RR, Richmond TA, Pofahl NJ, Green RD, Eis PS, Nair P, Brothman AR, Stallings RL: Analysis of chromosome breakpoints in neuroblastoma at sub-kilobase resolution using fine-tiling oligonucleotide array CGH. Genes Chromosomes Cancer. 2005, 44: 305-319. 10.1002/gcc.20243

    Article  CAS  PubMed  Google Scholar 

  32. McCormick MR, Selzer RR, Richmond TA: Methods in high-resolution, array-based comparative genomic hybridization. Methods Mol Biol. 2007, 381: 189-211. full_text

    CAS  PubMed  Google Scholar 

  33. Baba Y, Nosho K, Shima K, Irahara N, Kure S, Toyoda S, Kirkner GJ, Goel A, Fuchs CS, Ogino S: Aurora-A expression is independently associated with chromosomal instability in colorectal cancer. Neoplasia. 2009, 11: 418-425.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Yu X, Harris SL, Levine AJ: The regulation of exosome secretion: a novel function of the p53 protein. Cancer Res. 2006, 66: 4795-4801. 10.1158/0008-5472.CAN-05-4579

    Article  CAS  PubMed  Google Scholar 

  35. Cummins JM, He Y, Leary RJ, Pagliarini R, Diaz LA, Sjoblom T, Barad O, Bentwich Z, Szafranska AE, Labourier E, Raymond CK, Roberts BS, Juhl H, Kinzler KW, Vogelstein B, Velculescu VE: The colorectal microRNAome. Proc Natl Acad Sci USA. 2006, 103: 3687-3692. 10.1073/pnas.0511155103

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Hermsen M, Postma C, Baak J, Weiss M, Rapallo A, Sciutto A, Roemen G, Arends JW, Williams R, Giaretti W, De GA, Meijer G: Colorectal adenoma to carcinoma progression follows multiple pathways of chromosomal instability. Gastroenterology. 2002, 123: 1109-1119. 10.1053/gast.2002.36051

    Article  CAS  PubMed  Google Scholar 

  37. Camps J, Grade M, Nguyen QT, Hormann P, Becker S, Hummon AB, Rodriguez V, Chandrasekharappa S, Chen Y, Difilippantonio MJ, Becker H, Ghadimi BM, Ried T: Chromosomal breakpoints in primary colon cancer cluster at sites of structural variants in the genome. Cancer Res. 2008, 68: 1284-1295. 10.1158/0008-5472.CAN-07-2864

    Article  CAS  PubMed  Google Scholar 

  38. Nakao K, Mehta KR, Fridlyand J, Moore DH, Jain AN, Lafuente A, Wiencke JW, Terdiman JP, Waldman FM: High-resolution analysis of DNA copy number alterations in colorectal cancer by array-based comparative genomic hybridization. Carcinogenesis. 2004, 25: 1345-1357. 10.1093/carcin/bgh134

    Article  CAS  PubMed  Google Scholar 

  39. Hertz S, Rothamel T, Skawran B, Giere C, Steinemann D, Flemming P, Becker T, Flik J, Wiese B, Soudah B, Kreipe H, Schlegelberger B, Wilkens L: Losses of chromosome arms 4q, 8p, 13q and gain of 8q are correlated with increasing chromosomal instability in hepatocellular carcinoma. Pathobiology. 2008, 75: 312-322. 10.1159/000151712

    Article  PubMed  Google Scholar 

  40. Cetin E, Cengiz B, Gunduz E, Gunduz M, Nagatsuka H, Bekir-Beder L, Fukushima K, Pehlivan D, MO N, Nishizaki K, Shimizu K, Nagai N: Deletion mapping of chromosome 4q22-35 and identification of four frequently deleted regions in head and neck cancers. Neoplasma. 2008, 55: 299-304.

    CAS  PubMed  Google Scholar 

  41. Krona C, Caren H, Sjoberg RM, Sandstedt B, Laureys G, Kogner P, Martinsson T: Analysis of neuroblastoma tumour progression; loss of PHOX2B on 4p13 and 17q gain are early events in neuroblastoma tumourigenesis. Int J Oncol. 2008, 32: 575-583.

    CAS  PubMed  Google Scholar 

  42. Shivapurkar N, Sood S, Wistuba II, Virmani AK, Maitra A, Milchgrub S, Minna JD, Gazdar AF: Multiple regions of chromosome 4 demonstrating allelic losses in breast carcinomas. Cancer Res. 1999, 59: 3576-3580.

    CAS  PubMed  Google Scholar 

  43. Polascik TJ, Cairns P, Chang WY, Schoenberg MP, Sidransky D: Distinct regions of allelic loss on chromosome 4 in human primary bladder carcinoma. Cancer Res. 1995, 55: 5396-5399.

    CAS  PubMed  Google Scholar 

  44. Ried T, Knutzen R, Steinbeck R, Blegen H, Schrock E, Heselmeyer K, du MS, Auer G: Comparative genomic hybridization reveals a specific pattern of chromosomal gains and losses during the genesis of colorectal tumors. Genes Chromosomes Cancer. 1996, 15: 234-245. 10.1002/(SICI)1098-2264(199604)15:4<234::AID-GCC5>3.0.CO;2-2

    Article  CAS  PubMed  Google Scholar 

  45. Jiang LX, Xu J, Wang ZW, Li DP, Peng ZH, Gao JJ, He L, Zheng HT: Tumor suppress genes screening analysis on 4q in sporadic colorectal carcinoma. World J Gastroenterol. 2008, 14: 5606-5611. 10.3748/wjg.14.5606

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  46. Grimwood J, Gordon LA, Olsen A, Terry A, Schmutz J, Lamerdin J, Hellsten U, Goodstein D, Couronne O, Tran-Gyamfi M, Aerts A, Altherr M, Ashworth L, Bajorek E, Black S, Branscomb E, Caenepeel S, Carrano A, Caoile C, Chan YM, Christensen M, Cleland CA, Copeland A, Dalin E, Dehal P, Denys M, Detter JC, Escobar J, Flowers D, Fotopulos D: The DNA sequence and biology of human chromosome 19. Nature. 2004, 428: 529-535. 10.1038/nature02399

    Article  CAS  PubMed  Google Scholar 

  47. Vermeersch P, Zachee P, Brusselmans C: Acute myeloid leukemia with bone marrow necrosis and Charcot Leyden crystals. Am J Hematol. 2007, 82: 1029- 10.1002/ajh.20907

    Article  PubMed  Google Scholar 

  48. Xu Y, Pasche B: TGF-beta signaling alterations and susceptibility to colorectal cancer. Hum Mol Genet. 2007, 16: R14-R20. 10.1093/hmg/ddl486

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Oklu R, Hesketh R: The latent transforming growth factor beta binding protein (LTBP) family. Biochem J. 2000, 352: 601-610. 10.1042/0264-6021:3520601

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Siegel PM, Massague J: Cytostatic and apoptotic actions of TGF-beta in homeostasis and cancer. Nat Rev Cancer. 2003, 3: 807-821. 10.1038/nrc1208

    Article  CAS  PubMed  Google Scholar 

  51. Mauel S, Kruse B, Etschmann B, Schulenburg von der AG, Schaerig M, Stovesand K, Wilcken B, Sterner-Kock A: Latent transforming growth factor binding protein 4 (LTBP-4) is downregulated in human mammary adenocarcinomas in vitro and in vivo. APMIS. 2007, 115: 687-700. 10.1111/j.1600-0463.2007.apm_453.x

    Article  CAS  PubMed  Google Scholar 

  52. Sterner-Kock A, Thorey IS, Koli K, Wempe F, Otte J, Bangsow T, Kuhlmeier K, Kirchner T, Jin S, Keski-Oja J, von MH: Disruption of the gene encoding the latent transforming growth factor-beta binding protein 4 (LTBP-4) causes abnormal lung development, cardiomyopathy, and colorectal cancer. Genes Dev. 2002, 16: 2264-2273. 10.1101/gad.229102

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  53. Mounier CM, Wendum D, Greenspan E, Flejou JF, Rosenberg DW, Lambeau G: Distinct expression pattern of the full set of secreted phospholipases A2 in human colorectal adenocarcinomas: sPLA2-III as a biomarker candidate. Br J Cancer. 2008, 98: 587-595. 10.1038/sj.bjc.6604184

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  54. Laye JP, Gill JH: Phospholipase A2 expression in tumours: a target for therapeutic intervention?. Drug Discov Today. 2003, 8: 710-716. 10.1016/S1359-6446(03)02754-5

    Article  CAS  PubMed  Google Scholar 

  55. Morioka Y, Ikeda M, Saiga A, Fujii N, Ishimoto Y, Arita H, Hanasaki K: Potential role of group × secretory phospholipase A(2) in cyclooxygenase-2-dependent PGE(2) formation during colon tumorigenesis. FEBS Lett. 2000, 487: 262-266. 10.1016/S0014-5793(00)02350-4

    Article  CAS  PubMed  Google Scholar 

  56. Burke JE, Dennis EA: Phospholipase A2 structure/function, mechanism, and signaling. J Lipid Res. 2009, 50 (Suppl): S237-S242. 10.1194/jlr.R800033-JLR200

    PubMed Central  PubMed  Google Scholar 

  57. Fan S, Ramalingam SS, Kauh J, Xu Z, Khuri FR, Sun SY: Phosphorylated eukaryotic translation initiation factor 4 (eIF4E) is elevated in human cancer tissues. Cancer Biol Ther. 2009, 8: 1463-1469.

    Article  PubMed Central  PubMed  Google Scholar 

  58. Martin ME, Perez MI, Redondo C, Alvarez MI, Salinas M, Fando JL: 4E binding protein 1 expression is inversely correlated to the progression of gastrointestinal cancers. Int J Biochem Cell Biol. 2000, 32: 633-642. 10.1016/S1357-2725(00)00007-8

    Article  CAS  PubMed  Google Scholar 

  59. Feng Z, Hu W, de SE, Teresky AK, Jin S, Lowe S, Levine AJ: The regulation of AMPK beta1, TSC2, and PTEN expression by p53: stress, cell and tissue specificity, and the role of these gene products in modulating the IGF-1-AKT-mTOR pathways. Cancer Res. 2007, 67: 3043-3053. 10.1158/0008-5472.CAN-06-4149

    Article  CAS  PubMed  Google Scholar 

  60. Wessels LF, van WT, Hart AA, van't Veer LJ, Reinders MJ, Nederlof PM: Molecular classification of breast carcinomas by comparative genomic hybridization: a specific somatic genetic profile for BRCA1 tumors. Cancer Res. 2002, 62: 7110-7117.

    CAS  PubMed  Google Scholar 

  61. Jonsson G, Naylor TL, Vallon-Christersson J, Staaf J, Huang J, Ward MR, Greshock JD, Luts L, Olsson H, Rahman N, Stratton M, Ringner M, Borg A, Weber BL: Distinct genomic profiles in hereditary breast tumors identified by array-based comparative genomic hybridization. Cancer Res. 2005, 65: 7612-7621.

    PubMed  Google Scholar 

  62. Peltomaki P, Lothe RA, Aaltonen LA, Pylkkanen L, Nystrom-Lahti M, Seruca R, David L, Holm R, Ryberg D, Haugen A: Microsatellite instability is associated with tumors that characterize the hereditary non-polyposis colorectal carcinoma syndrome. Cancer Res. 1993, 53: 5853-5855.

    CAS  PubMed  Google Scholar 

  63. Blaker H, Mechtersheimer G, Sutter C, Hertkorn C, Kern MA, Rieker RJ, Penzel R, Schirmacher P, Kloor M: Recurrent deletions at 6q in early age of onset non-HNPCC- and non-FAP-associated intestinal carcinomas. Evidence for a novel cancer susceptibility locus at 6q14-q22. Genes Chromosomes Cancer. 2008, 47: 159-164. 10.1002/gcc.20516

    Article  PubMed  Google Scholar 

  64. Skoglund J, Djureinovic T, Zhou XL, Vandrovcova J, Renkonen E, Iselius L, Bisgaard ML, Peltomaki P, Lindblom A: Linkage analysis in a large Swedish family supports the presence of a susceptibility locus for adenoma and colorectal cancer on chromosome 9q22.32-31.1. J Med Genet. 2006, 43: e7- 10.1136/jmg.2005.033928

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  65. Gruber SB, Moreno V, Rozek LS, Rennerts HS, Lejbkowicz F, Bonner JD, Greenson JK, Giordano TJ, Fearson ER, Rennert G: Genetic variation in 8q24 associated with risk of colorectal cancer. Cancer Biol Ther. 2007, 6: 1143-1147. 10.1158/1535-7163.MCT-06-0555

    Article  CAS  PubMed  Google Scholar 

  66. Haiman CA, Le ML, Yamamato J, Stram DO, Sheng X, Kolonel LN, Wu AH, Reich D, Henderson BE: A common genetic risk factor for colorectal and prostate cancer. Nat Genet. 2007, 39: 954-956. 10.1038/ng2098

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This study was funded by grants from the Norwegian Cancer Society (95068/RAL) supporting MB and THÅ as research fellows. Kari Almendingen, Arne Bakka, Gunter Bock, Torunn Fetveit, Hans Joacim Haus, Anders Husby, Ragnhild A. Lothe, Tom Mala, Øystein Mathisen, Ingvild Moberg, Arild Nesbakken, Arve Rennesund, Oddvar Sandvik, Espen Thiis-Evensen, Morten H. Vatn are acknowledged as members of the INFAC (individuals with early onset and/or familial cancer) -study group.

Author information

Authors and Affiliations

Authors

Consortia

Corresponding author

Correspondence to Ragnhild A Lothe.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

MB and THÃ… have performed the experimental analyses. MB has drafted the manuscript, made all figures and performed statistical- and computational analysis. ETE, MAM, MHV, AN and the INFAC-study group have provided the tumor series and the patient information. MRT participated in the evaluation of results and in the manuscript preparation. ETE, AN, and RIS participated in the study design. RIS supervised the bioinformatic analyses, participated in evaluation of results and in the manuscript preparation. RAL conceived the study, was main responsible for the study design and study coordination, participated in evaluation of results and in manuscript preparation. All authors have contributed to the writing and have read and approved the final version of the manuscript.

Electronic supplementary material

Additional file 1: Summary of clinical data for all patients included in the study. (XLS 18 KB)

12943_2009_615_MOESM2_ESM.XLS

Additional file 2:Primer details for TP53 mutation analysis. Details of the fragments amplified for TP53 analysis. Both 5'- and 3' primers have M13 tails, indicated in capital letters. (XLS 26 KB)

12943_2009_615_MOESM3_ESM.PDF

Additional file 3:Detailed profile of gains and losses in chromosome 8. Percentage of gains and losses along chromosome 8 for all colorectal carcinomas in the presented study. Centromere position is indicated by solid line, breakpoint region by dashed line. (PDF 138 KB)

12943_2009_615_MOESM4_ESM.XLS

Additional file 4:Overview of aberrant genomic regions in each carcinoma. Detailed listing of regions displaying gain or loss for each patient sample. (XLS 226 KB)

12943_2009_615_MOESM5_ESM.XLS

Additional file 5:Loci within the chromosome bands that show statistically different copy number changes in colorectal carcinomas between patient groups. Protein coding genes, miRNAs and pseudogenes in all regions with statistically significant difference in copy number between patient groups are extracted by using BioMart with Ensemble version 53. (XLS 34 KB)

12943_2009_615_MOESM6_ESM.XLS

Additional file 6:Genes that show corresponding mRNA expression and DNA copy number changes in early onset CRCs. Genes (with respective probe IDs) with mRNA expression levels corresponding to DNA copy number, also situated in regions with statistically significant difference in early onset compared to late onset patients are listed. The expression ratio is calculated as median expression in samples with gain or loss divided by median expression in samples with normal copy number. The number of tumors with copy number aberration in the early onset and late onset groups are indicated, as well as the P-value. The resulting genes are primarily sorted by increasing P-value, based on aCGH-analysis alone, followed by a sort based on the expression-ratio. (XLS 42 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Berg, M., Ã…gesen, T.H., Thiis-Evensen, E. et al. Distinct high resolution genome profiles of early onset and late onset colorectal cancer integrated with gene expression data identify candidate susceptibility loci. Mol Cancer 9, 100 (2010). https://doi.org/10.1186/1476-4598-9-100

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1476-4598-9-100

Keywords