<title>Introduction</title> <p>Cancer is driven by mutations in genes that control the proliferation of cells, their survival and their integrity. Screens aimed at identifying such cancer genes often use chromosomal location and/or functional properties to select candidates genes for subsequent mutation analysis ##REF##8553070##[1]##–##REF##17344846##[4]##. Although many candidate cancer gene loci have been identified, the labor-intensive mutation analysis severely hampers finding the corresponding cancer gene. Other gene search strategies have focused on aberrant gene expression patterns to identify candidates. For example, gene mutants that result in premature termination codons were identified by screening for genes that were specifically expressed following chemical inhibition of nonsense mediated RNA decay ##REF##15300251##[5]##. Furthermore, fusion genes in prostate cancer were identified by screening for outliers in a large cohort of gene-expression profiles ##REF##16254181##[6]##.</p> <p>Human cancer gene mutations frequently result in the skipping of one or several exons from the encoded transcripts ##REF##8033105##[7]##–##REF##12600935##[9]##. Exon-skipping mutations may be caused by nucleotide substitutions within the consensus splice sites or by deletions that span entire exons. In addition, exon-skipping mutations may be caused by relatively small intragenic insertions, deletions or duplications. Even though exon-skipping mutations represent an estimated 10–20% of all cancer-related gene mutations ##REF##17344846##[4]##, ##REF##12600935##[9]##–##REF##16839875##[12]##, no high throughput method has been available to screen for such mutations. Here, we describe Pattern Based Correlation (PAC) as an approach to identify candidate cancer genes by screening for exon-skipping events in a global fashion. Detailed mutation analysis is then restricted only to the PAC-identified outlier exons. As a proof-of-principle, we demonstrate the efficacy of the PAC strategy on previously identified exon-skipping mutations in breast cancer cell lines and in clinical brain tumor samples. We also demonstrate that PAC can identify novel exon skipping events with underlying genetic changes in known cancer genes and in randomly-selected PAC-identified outlier exons.</p>
<title>Materials and Methods</title> <title>Samples</title> <p>Our collection of 41 publicly-available human breast cancer cell lines had been subjected to mutational screens of seven tumor suppressor genes: <italic>BRCA1</italic> (Breast Cancer Susceptibility Gene 1; OMIM 113705), <italic>CDH1</italic> (E-cadherin; OMIM 192090), <italic>MAP2K4</italic> (MAP Kinase Kinase 4, a.k.a. <italic>MKK4</italic>; OMIM 601335), <italic>PTEN</italic> (Phosphatase and Tensin Homolog; OMIM 601728), <italic>p16</italic> (CDK4-inhibitor, a.k.a. <italic>INK4A</italic>, <italic>CDKN2A</italic>; OMIM 600160), <italic>p53</italic> (Tumor Protein p53; OMIM 191170) and <italic>RB1</italic> (Retinoblastoma Susceptibility Gene 1; OMIM 180200) ##REF##16397213##[14]##–##REF##16541312##[18]## (Nagel <italic>et al.</italic> submitted for publication). Mutational analysis involved sequencing the entire coding region of these genes on genomic DNA as well as analysis of the encoded transcript. The twelve breast cancer cell lines used for this study were: CAMA-1, EVSA-T, HCC1937, MDA-MB-134VI, MDA-MB-157, MDA-MB-435s, MDA-MB-436, MDA-MB-453, MDA-MB-468, MPE600, OCUB-F and SK-BR-5. Clinical glioblastoma specimens were frozen in liquid nitrogen immediately upon surgical resection from patients at Erasmus University Medical Center, as described elsewhere ##REF##17575129##[13]##. Pathological review revealed at least 70% tumor nuclei for each specimen. Mutation analysis of the <italic>EGFR</italic> oncogene (Epidermal Growth Factor Receptor; OMIM 131550) in the glioblastomas was performed by conventional RT-PCR and subsequent sequencing of transcripts from samples with <italic>EGFR</italic> amplifications. <italic>EGFR</italic> transcript expression was quantified by Real-Time RT-PCR, using primers that amplified exons 2–3 or exons 22–23 and thus allowed discrimination of wild-type <italic>EGFR</italic> transcripts and the <italic>EGFRvIII</italic> isoform.</p> <title>Exon-level expression profiling</title> <p>Total RNA was isolated using the Qiagen RNeasy kit for the breast cancer cell lines and using Trizol followed by RNeasy for the glioblastoma specimens ##REF##16357140##[23]##. RNA quality was assessed using the Agilent Bioanalyser, requiring RNA integrity &gt;7.0 ##REF##16448564##[24]##. All further processing of the samples was performed according the Affymetrix GeneChip Whole Transcript (WT) Sense Target Labeling Assay. Affymetrix GeneChip Human Exon 1.0 ST Arrays were used to determine the expression level of virtually all exons present in the human genome (1.4 million probe sets covering &gt;1 million exon clusters). For this study, we used expression data of the 290,000 core probe sets that are supported by putative full-length mRNA from e.g. the RefSeq database (Geo dataset accession number GSE9385) . Signal processing was performed after sketch normalization by using Affymetrix ExACT 1.2.1 software and the PLIER algorithm, described in Affymetrix GeneChip Exon Array Whitepaper “Gene Signal Estimates from Exon Arrays” and Technote “Guide to Probe Logarithmic Intensity Error (PLIER) Estimation” (<ext-link ext-link-type="uri" xlink:href="http://www.affymetrix.com/support/technical">www.affymetrix.com/support/technical</ext-link>).</p> <title>PAttern based Correlation (PAC)</title> <p>Predicted exon expression levels were calculated by using the PAC algorithm, described in Whitepaper “Alternative Transcript Analysis Methods for Exon Arrays”, where the predicted expression level of the exon (Exon-pr) equals the overall expression of its transcript in that sample (Transcript-m: the meta probe set expression level) multiplied by the average expression level of that exon among all samples (Exon-ave) and divided by the average overall expression of the transcript among all samples (Transcript-ave), all 2-logarithm transformed. In formula:
</p> <p>PAC values were calculated by subtracting the predicted expression level of the exon in that sample from its measured expression level (Exon-m), again with 2-logarithm transformation:
</p> <p>Meta probe set expression levels were calculated using all core probe sets of a transcript with PLIER signal estimates &gt;50. To enrich for probe sets with significant expression above background, PAC values were calculated using exons and transcripts that had PLIER signal estimates &gt;50 ##REF##17575129##[13]##. Identification of outlier exons was performed without prior knowledge of the mutation data.</p>
<title>Results</title> <title>Outlier exon identification by Pattern Based Correlation (PAC)</title> <p>In this study we have developed a new approach to screen exon-skipping events in human cancer samples. Because mutations in cancer often are highly heterogeneous with respect to their intragenic location, individual tumors often express unique RNA species. Screening for mutations that result in skipping of one or more exons in the encoded transcript therefore requires screening for unique, exon-skipped, transcripts within a specific sample cohort. Briefly, exon-level expression profiles are generated using Affymetrix Human Exon Arrays, which determine the expression level of virtually all exons present in the human genome. The PAttern based Correlation (PAC) algorithm is then used to calculate the predicted expression level of each exon (or probe set). We then identify outlier exons by subtracting the predicted expression level of exons from their measured expression level, with values equaling zero when the measured expression level of an exon was similar to its predicted expression level (formulated in detail under Methods). The PAC algorithm has been used to identify differential splicing between predefined groups ##REF##17575129##[13]##. In this study, we have tested the PAC algorithm for its ability to identify aberrantly expressed exons in single samples of a well defined cohort of cell lines or tumors. PAC effectively normalizes the variability in gene expression levels between samples and, in a single sample, normalizes the variability in signal intensity between probe sets of the same transcript (##FIG##0##Fig. 1##).</p> <title>PAC detects exon-skipping events in breast cancer cell lines</title> <p>We tested the feasibility of the PAC strategy on a panel of 12 human breast cancer cell lines that had been screened for mutations in seven tumor suppressor genes: <italic>BRCA1</italic>, <italic>CDH1</italic>, <italic>MAP2K4</italic>, <italic>PTEN</italic>, <italic>p16</italic>, <italic>p53</italic> and <italic>RB1</italic>
##REF##16397213##[14]##–##REF##16541312##[18]##, and unpublished results). Mutation analysis was performed by sequencing of the complete coding sequences of the genes and analysis of all mutations on both genomic gene fragments and transcripts. Together, the 12 cell lines contained seven gene mutants that should be detectable by PAC, as they resulted in the skipping of eight exons from among four tumor suppressor genes (mutations are detailed in ##TAB##0##Table 1##). We have explored the PAC strategy at different cut-off levels, identifying outlier exons that were expressed less than 16-fold, 8-fold, 4-fold, 2.8-fold and 2.5-fold than their predicted expression level (i.e. PAC values of -4.0, -3.0, -2.0, -1.5 and -1.3, respectively). Outlier exons were identified without prior knowledge of the mutation data.</p> <p>From the total of 3.4 million core probe sets that we assayed for the 12 cell lines (290,000 core probe sets per sample), PAC identified 21,151 (0.6%) outlier probe sets at PAC value −4.0 and 94,590 (2.8%) outlier probe sets at PAC value −1.3 (##FIG##1##Fig. 2A##). All probe sets at PAC values &lt;−2.0 (34,137 probe sets corresponding to 31,357 exons and 10,247 genes) are listed in supplementary data ##SUPPL##0##Table S1##. When all PAC values are plotted in a frequency histogram, a tail towards the negative end is observed (##FIG##1##Fig. 2A##). This skewed distribution gives a rough estimation of the false positive rate at the various PAC levels. PAC of the seven fully characterized tumor suppressor genes in the 12 cell lines involved analysis of 1,200 exons (1,752 probe sets). PAC correctly detected six of the eight skipped exons when using PAC value −4.0, seven skipped exons were detected at PAC value −2.0 and all eight skipped exons were detected at PAC value −1.3 (##FIG##1##Fig. 2C##). Importantly, the number of false positive outlier exons was substantially reduced at PAC value −4.0 as compared to PAC value −1.3, resulting in an increase of the true positive rate from 9% to 24% of the identified outlier exons (##FIG##1##Fig. 2D##). For comparison, random sampling of 24/1200 exons has a &gt;85% probability of not finding any true positive mutation and only a &lt;10<sup>−8</sup> chance of finding 6 or more. For the known cancer genes used in our study, the true positive rate of our approach thus by far exceeds random exon selection. In this respect, the reduction of the number of false positive candidate genes may initially be far more beneficial for a gene search project than accurate identification of all true positive outlier exons. Together, our results show that the PAC strategy is reliable in detecting exon-skipping mutants in cancer cell lines.</p> <title>PAC performance in samples with heterogeneous transcript expression</title> <p>Similar to other genetic screening methods, PAC is most suited to detect homozygous genetic changes. For example, the lowest PAC value when 50% wild-type transcripts are present (as may be the case for a heterozygous genetic change) is −1.0. The somewhat compromised detection of skipped exons at PAC value −4.0 as compared to PAC value −1.3 (i.e. six vs. all eight skipped exons) in our panel of breast cancer cell lines therefore may have been caused by the expression of a second aberrant transcript that still includes (part of) the exon. Indeed, a second <italic>CDH1</italic> transcript length of minor intensity was detected in CAMA-1 (##FIG##2##Fig. 3A##), the splice site mutant that had been detected only at PAC value −1.3.</p> <p>To further asses the performance of PAC in samples with heterogeneous (wild-type and mutant) transcript expression, we performed PAC on 14 clinical glioblastoma specimens (selected to contain &gt;70% tumor nuclei) that had genomic amplifications of the <italic>EGFR</italic> oncogene. Glioblastomas with <italic>EGFR</italic> amplifications frequently carry an intragenic deletion of exons 2 through 7, resulting in expression of the constitutively active <italic>EGFRvIII</italic> isoform ##REF##11302336##[8]##, ##REF##10728703##[19]##. However, glioblastomas expressing the <italic>EGFRvIII</italic> isoform also frequently express wild-type <italic>EGFR</italic> transcripts. This heterogeneous <italic>EGFR</italic> expression is related to amplification of the <italic>EGFR</italic> locus prior to the deletion of exons ##REF##15700833##[20]##, although non-malignant cells in the glioblastoma specimens may also express <italic>EGFR</italic>. Of the fourteen glioblastoma samples used in this study, six expressed <italic>EGFRvIII</italic> (a total of 36 skipped exons) of which five also expressed significant levels of wild-type <italic>EGFR</italic> transcripts as determined by quantitative Real-Time PCR (qPCR) (##FIG##2##Fig. 3B##) (insufficient RNA remained of the sixth sample with <italic>EGFRvIII</italic> expression to perform qPCR).</p> <p>From the total of 4.1 million core probe sets that we assayed for these 14 samples (290,000 core probe sets per sample), PAC identified 1,646 (0.04%) outlier probe sets at PAC value −4.0 and 39,936 (1.0%) outlier probe sets at PAC value −1.3 (##FIG##1##Fig. 2B##). PAC thus identified three to ten-fold less outlier exons in the glioblastomas as compared to the breast cancer cell lines (##FIG##1##Fig. 2A##). All probe sets at PAC values &lt;−2.0 (11,287 probe sets, corresponding to 10,903 exons and 6,264 genes) are listed in supplementary data ##SUPPL##0##Table S1##. This smaller number of outlier exons in the glioblastomas may be related to their homogeneous histopathology and their highly similar gene expression profiles ##REF##17575129##[13]##, ##REF##15378038##[21]##, to the presence of non-neoplastic cells in the tumor samples, or may reflect sampling biases due to small cohort sizes.</p> <p>PAC of the <italic>EGFR</italic> gene in the 14 glioblastomas involved the analysis of 392 exons (434 probe sets). PAC detected two of six <italic>EGFRvIII</italic> expressing tumors (12 of the 36 skipped exons) at PAC values −2.0 and lower (##TAB##0##Table 1## and ##FIG##1##Fig. 2C##). Of the two glioblastomas with <italic>EGFRvIII</italic> that had been detected by PAC, one had significantly (i.e. &gt;5 fold) more mutant than wild-type <italic>EGFR</italic> transcripts. In this tumor, the Ct value difference was &gt;2 between qPCR fragments inside (measuring only wild-type <italic>EGFR</italic> transcripts) and outside (measuring both wild-type and <italic>EGFRvIII</italic> transcripts) the <italic>EGFR</italic> exon 2–7 deletion region (##FIG##2##Fig. 3B##). The other glioblastoma had a similar expression level difference between wild-type and <italic>EGFRvIII</italic> transcripts (Ct value difference of ∼1.5) as the three glioblastomas that had not been detected by PAC, but had lower overall <italic>EGFR</italic> transcript levels. It appears that PAC detection of the <italic>EGFRvIII</italic> isoform is determined by the overall expression level of <italic>EGFR</italic> transcripts in combination with the ratio of <italic>EGFRvIII</italic> and wild-type <italic>EGFR</italic> transcripts, where samples with too high <italic>EGFR</italic> transcript levels may escape PAC detection due to saturation of the probe sets involved. These results show that the PAC strategy can detect exon-skipping mutants in clinical cancer specimens if the ratio mutant/wild-type transcript level is high and when probe sets are within the linear detection range of the microarray.</p> <title>PAC performance in detecting recurrent outlier exons</title> <p>PAC performance can also be challenged by recurrent outlier exons. Such frequently skipped exons will result in an underestimation of the exon/transcript ratio in the PAC algorithm and so increase PAC values. We therefore evaluated the performance of PAC in detecting recurrent outlier exons by reiterated replacement of <italic>EGFRvIII</italic> expressing samples with samples that expressed only wild-type <italic>EGFR</italic> (##FIG##3##Fig. 4A##). When six of 14 samples express <italic>EGFRvIII</italic>, the deletion of exons 2–7 in GBM67 is not detected by PAC. PAC values indeed decreased with decreasing ratios of wild-type versus mutant samples. However, the decrease was relatively small and resulted in the identification of only one of the six deleted exons once the ratio had dropped to one mutant sample among 14 samples. We also simulated PAC detection of recurrent mutations with two breast cancer cell lines, of which HCC1937 had skipped <italic>RB1</italic> exon 22, and we were already able to identify the mutant from among two samples up to even five mutants from among six samples (##FIG##3##Fig. 4B##). These simulation experiments indicate that PAC performs well in identifying recurrent exon-skipping mutations.</p> <title>Detection of nucleotide substitutions and novel genetic changes by PAC</title> <p>The performance of PAC was further evaluated by analysis of outlier exons selected from all candidates at PAC value ≤−2.0 in breast cancer cell lines and clinical glioblastoma samples. In total, 44 and 68 outlier exons were screened in breast cancer cell lines and glioblastoma samples respectively. Sequence analysis of PCR amplified outlier exons identified two novel exon skipping events and two novel genetic base changes in glioblastoma samples, as well as a number of previously reported base changes (homozygous SNPs) in breast cancer cell lines (n = 5) and glioblastomas (n = 16). RT-PCR experiment results are detailed in supplementary data ##SUPPL##1##Table S2##.</p> <p>The majority of genetic changes identified by PAC were single nucleotide changes, both in breast cancer cell lines (five known SNPs) and in glioblastomas (two novel base changes and 16 known SNPs). Moreover, two out of ten previously identified oncogenic point mutations that did not result in exon skipping events were also PAC detected in our cohort of breast cancer cell lines: <italic>MAP2K4</italic> c.551C&gt;G in MDA-MB-134VI and <italic>PTEN</italic> c.274G&gt;C in CAMA-1; ##UREF##0##[16]##, ##REF##17314276##[17]## (##FIG##4##Fig. 5##). Single nucleotide mismatches have been used to define hybridization specificity on other Affymetrix microarray platforms. By analogy, single nucleotide substitutions in cancer may also cause reduced hybridization to the probes on the microarray and thus be detected as outlier exons by PAC. Indeed, all of the PAC detected base changes and SNPs were centrally localized within the probe set selection region and overlap with several of its individual probes (##FIG##4##Fig. 5##).</p> <p>One of the PAC identified novel exon skipping events was predicted to result in a deletion of the four 3′ end exons of <italic>EGFR</italic> (##FIG##5##Fig. 6A##). This exon-skipping event was due to a genomic deletion as determined using semi quantitative PCR on genomic tumor DNA. Compared to the 5′ end of the <italic>EGFR</italic> locus in GBM157, the 3′ end showed less amplification (ΔCt −2.5) whereas other samples showed equal amplification between the 5′ and 3′ end of the gene (ΔCt 0.3±1.9). Similar 3′ deletions in <italic>EGFR</italic> have been observed previously in gliomas ##REF##10728703##[19]##. The second exon-skipping event predicted by PAC would result in a deletion of exon 30 in the <italic>FCGBP</italic> cDNA (##FIG##5##Fig. 6B##). This deletion will cause a frameshift that is predicted to result in a truncated protein. The absence of exon 30 was confirmed by RT-PCR and sequence analysis (##FIG##5##Fig. 6C##). Novel identified single base changes include a single base change 1934C&gt;G (S645C) in the <italic>EGFR</italic> gene, (##FIG##5##Fig. 6A and D##), and a single base change 946G&gt;A (G316R) in the <italic>TLE2</italic> gene (##FIG##5##Fig. 6E##). The G316R (946G&gt;A) mutation in <italic>TLE2</italic> is rendered “pathological” by PMut (mmb2.pcb.ub.es:8080/PMut/) and “not tolerated” by SIFT BLink (blocks.fhcrc.org/sift/SIFT_BLink_submit.html). In summary, the novel exon skipping events and base changes identified by analysis of a selected set of outlier exons confirms the suitability of PAC to identify candidate cancer genes.</p>
<title>Discussion</title> <p>We have developed an approach that uses PAttern based Correlation (PAC) to screen for cancer gene mutations that cause exon skipping in the encoded transcripts. We demonstrate that PAC correctly detected all of seven previously identified exon-skipping mutants in breast cancer cell lines and two of six mutants in clinical glioblastoma samples. The true and false positive rates were determined at various stringency levels. Importantly, PAC identified a number of novel genetic changes, including those affecting splicing, that previously had gone undetected. These novel genetic changes are either in known cancer genes (<italic>EGFR</italic>), result in a frameshift (<italic>FCGBP</italic>) or are rendered “not tolerated” by gene prediction algorithms PMut and SIFT BLink (<italic>TLE2</italic>). Additional experiments are required to determine whether changes in the novel candidate cancer genes (<italic>FCGBP</italic> and <italic>TLE2</italic>) are indeed oncogenic. A significant number of nucleotide substitutions that are located within the probe set selection region are also PAC detected (##FIG##4##Fig. 5##). Our results thus classify PAC as a reliable approach to screen for candidate cancer genes in a global fashion.</p> <p>Gene expression profiling at the level of individual exons has only recently become feasible through the release of exon arrays. Here, we have explored the efficacy of PAC to identify exon-skipping mutants, but the strategy may also be used to deduce the primary structure of gene transcripts ##REF##17575129##[13]##, ##REF##17456239##[22]##. It is important to note that the PAC algorithm, detailed under Materials and Methods, is in essence a simple formula that predicts outlier exons based on fold change differences between measured and predicted exon expression levels. Other approaches can also be used to identify outliers (e.g. &gt;n standard deviations from the mean expression level) but need to account for the non-linearity in gene expression levels between samples (especially for cancer genes) and the limited sample size. Because of the high true positive rates obtained by PAC, we did not further explore alternative statistical approaches.</p> <p>The PAC algorithm is independent of array platform or organism, allowing application of the PAC strategy in a wide variety of biological systems. Several algorithms for exon-level expression profiling are commercially available, including Stratagene ArrayAssist (<ext-link ext-link-type="uri" xlink:href="http://www.stratagene.com">www.stratagene.com</ext-link>), Partek Genomics Suite (<ext-link ext-link-type="uri" xlink:href="http://www.partek.com">www.partek.com</ext-link>) and Genomatix Suite (<ext-link ext-link-type="uri" xlink:href="http://www.genomatix.de">www.genomatix.de</ext-link>). Although each of these software packages is relatively straight-forward, important advantages of PAC are that it allows detection of unique outlier exons without any prior knowledge of the encoding gene or its transcript structure and that it does not require predefined subgroups of samples with differential expression of the outlier exons.</p> <p>As with any global screening strategy, PAC has its preconditions for detecting outlier exons. First and foremost, identification of outlier exons requires their transcript expression level to be within the linear detection range of the exon array, which is determined by their transcript expression level as well as the hybridization efficiency and specificity of the probe sets involved. The constituency of the test samples is another consideration, particularly when both mutant and wild-type transcripts may be expressed. For example, the breast cancer cell line cohort included two splice site mutants that escaped detection by PAC because each had a second transcript length of major intensity that resulted from cryptic splicing (<italic>BRCA1</italic> c.5396+1G&gt;A in MDA-MB-436 ##REF##16397213##[14]## and <italic>p16</italic> c.150+2T&gt;C in MDA-MB-436 (Nagel <italic>et al.</italic>, submitted for publication). Furthermore, PAC detection of the <italic>EGFRvIII</italic> transcript isoform in clinical glioblastomas was determined by the overall expression level of <italic>EGFR</italic> transcripts, that was near the limits of linear detection in all five <italic>EGFRvIII</italic> glioblastomas, but also by the ratio of the <italic>EGFRvIII</italic> isoform versus wild-type <italic>EGFR</italic> transcripts (##FIG##2##Fig. 3B##). A corollary is that PAC performance may be compromised in detecting an outlier exon when wild-type transcripts represent more than one-fourth of all transcripts of that particular gene, which could be the case in tumor samples with less than 75% neoplastic cells. However, expression levels of mutant and wild-type alleles typically are disproportional to their allele frequency and detection by PAC thus again is determined by the (relative) expression level of the outlier transcript. PAC therefore performs best in the absence of wild-type transcript expression. Homozygous transcripts are predominantly found among tumor suppressor genes, where often one allele is mutated accompanied by loss of the other allele.</p> <p>The influence of allele ratios was further stressed in our simulations of recurrent outlier detection by PAC: The <italic>EGFRvIII</italic> isoform in GBM67 was detected only once it was present as a unique outlier among 14 samples, whereas it had not been detected in our original PAC screen that included five other <italic>EGFRvIII</italic> expressing glioblastomas (##FIG##3##Fig. 4A##). However, this sub optimal PAC performance appeared not related to the recurrence of outliers, as recurrent outliers were easily identified among cell lines − even when present in five out of six cell lines (##FIG##3##Fig. 4B##). The simulation experiments also revealed that two cell lines were sufficient to reliably detect outlier exons and that more than eight cell lines did not further improve PAC performance, whereas for clinical tumor samples ten appeared the minimum but twenty would be preferred (##FIG##3##Fig. 4##).</p> <p>How efficient might PAC be in detecting mutations in cancer genomes? From our selection of outlier exons, we identified ∼20% (21/112) SNPs, ∼2% (2/112) novel base changes and ∼2% (2/112) exon skipping events. When including all nucleotide substitutions, the false positive rate in these experiments is ∼76%. By extension, amplification and sequencing 1,763 reactions on a single sample (all outliers at PAC values &lt;−4.0) can be expected to yield as much as 34 novel base changes and 34 exon skipping events. Therefore, our approach can be classified as a highly efficient screening method for candidate cancer genes, especially when compared to random selection of exons. Additional studies should then be performed to determine whether identified changes are causal for the tumor formation and/or progression, for example by screening for additional mutations (e.g. deletions, missense mutations) in other tumor samples or by functional analysis of the identified mutants.</p>
<p>Conceived and designed the experiments: MS PvdS PSS PF. Performed the experiments: FE LB JN ED AH MV MW PF. Analyzed the data: FE JKP PvdS PF. Wrote the paper: MS PF.</p> <title>Background</title> <p>Identification of genes that are causally implicated in oncogenesis is a major goal in cancer research. An estimated 10–20% of cancer-related gene mutations result in skipping of one or more exons in the encoded transcripts. Here we report on a strategy to screen in a global fashion for such exon-skipping events using PAttern based Correlation (PAC). The PAC algorithm has been used previously to identify differentially expressed splice variants between two predefined subgroups. As genetic changes in cancer are sample specific, we tested the ability of PAC to identify aberrantly expressed exons in single samples.</p> <title>Principal Findings</title> <p>As a proof-of-principle, we tested the PAC strategy on human cancer samples of which the complete coding sequence of eight cancer genes had been screened for mutations. PAC detected all seven exon-skipping mutants among 12 cancer cell lines. PAC also identified exon-skipping mutants in clinical cancer specimens although detection was compromised due to heterogeneous (wild-type) transcript expression. PAC reduced the number of candidate genes/exons for subsequent mutational analysis by two to three orders of magnitude and had a substantial true positive rate. Importantly, of 112 randomly selected outlier exons, sequence analysis identified two novel exon skipping events, two novel base changes and 21 previously reported base changes (SNPs).</p> <title>Conclusions</title> <p>The ability of PAC to enrich for mutated transcripts and to identify known and novel genetic changes confirms its suitability as a strategy to identify candidate cancer genes.</p>
<title>Supporting Information</title>
<fig id="pone-0003007-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003007.g001</object-id><label>Figure 1</label><caption><title>PAC detection of an exon-skipping <italic>PTEN</italic> mutant.</title><p>(A) Normalized expression data of all exons within the <italic>PTEN</italic> gene. Each exon probe set is represented by a dot in the solid line; multiple probe sets may be directed against the same exon. (B) PAC normalizes the variability in gene expression levels between samples and, in a single sample, the variability in signal intensity between probe sets of the same transcript. PAC calculation therefore allows rapid detection of skipping of <italic>PTEN</italic> exon 4 in breast cancer cell line MDA-MB-468 due to a <italic>PTEN</italic> c.253+1G&gt;T splice site mutation that we previously had identified ##REF##17314276##[17]##.</p></caption></fig> <fig id="pone-0003007-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003007.g002</object-id><label>Figure 2</label><caption><title>Performance of PAC to detect exon-skipping mutants.</title><p>(A) and (B) Total number of PAC-detected outlier probe sets from among 290,000 core probe sets in 12 breast cancer cell lines and in 14 glioblastomas, respectively. (C) Number of skipped exons detected by PAC as a percentage of all eight skipped exons present in the breast cancer cell lines, or as a percentage of the 36 skipped <italic>EGFR</italic> exons present in the glioblastomas (see ##TAB##0##Table 1##). (D) Total number of outlier exons (true plus false positives) and number of true positive outlier exons detected by PAC among the seven tumor suppressor genes and the <italic>EGFR</italic> oncogene. True positive outlier exons include all PAC detected skipped exons and two missense mutations (<italic>PTEN</italic> c.274G&gt;C in CAMA1, <italic>MAP2K4</italic> c.551C&gt;G in MDA-MB-134VI).</p></caption></fig> <fig id="pone-0003007-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003007.g003</object-id><label>Figure 3</label><caption><title>Compromised detection due to heterogeneous transcript expression.</title><p>Skipping of <italic>CDH1</italic> exon 11 in breast cancer cell line CAMA-1 was only detected at PAC value −1.3, likely due to expression of a second aberrant transcript variant (*) that was detected by conventional RT-PCR. (B) Expression of <italic>EGFR</italic> transcripts was detected in glioblastoma samples by Real-Time RT-PCR, using primers designed to anneal inside the exon 2–7 deletion region of the <italic>EGFRvIII</italic> isoform (gray bars) or outside the deletion region (black bars). Differences in Ct values between the two transcript fragments are indicative for <italic>EGFRvIII</italic> isoform expression levels. All five samples with the <italic>EGFRvIII</italic> isoform also expressed significant amounts of wild-type <italic>EGFR</italic> transcripts, likely compromising outlier detection by PAC (indicated by “detected” and “not detected”). Wild-type, samples with normal transcripts; Controls, non-malignant brain specimens.</p></caption></fig> <fig id="pone-0003007-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003007.g004</object-id><label>Figure 4</label><caption><title>Performance of PAC to detect recurrent outlier exons.</title><p>(A) Simulation experiment to determine PAC performance in detecting recurrent exon-skipping events among clinical glioblastoma samples, where mutant samples express the <italic>EGFRvIII</italic> isoform with deletion of exons 2 through 7. The cohort of 14 glioblastomas included six mutant samples that were replaced by wild-type samples through reiteration, based on their position from left to right in Fig. 3B. Deletion of <italic>EGFR</italic> exon 6 in sample GBM67 was detected only as unique mutant sample. (B) Simulation experiment to determine PAC performance in detecting recurrent exon-skipping events among breast cancer cell lines, using the wild-type cell line CAMA-1 and the <italic>RB1</italic> exon 22 deletion mutant HCC1937. The two cell lines were analyzed under various cohort sizes, with either the wild-type or the mutant cell line as single sample. The mutant sample was still detected at PAC value −2.0 with five recurrent mutants among six samples. The average expression level of <italic>RB1</italic> exon 22 dropped below PLIER 50 when more than five mutants were simulated, precluding PAC analysis (see <xref ref-type="sec" rid="s4">Materials and Methods</xref>).</p></caption></fig> <fig id="pone-0003007-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003007.g005</object-id><label>Figure 5</label><caption><title>Identification of nucleotide substitutions by PAC.</title><p>(A) PAC predicts skipping of the 5′ end of <italic>PTEN</italic> exon 5 in the CAMA-1 breast cancer cell line. This cell line contains a nucleotide substitution within the identified exon. This base change does not induce exon skipping but is centrally located within all three probes of the probe set (B). The central location suggests this mutation causes a reduced affinity to the probes on the exon-array.</p></caption></fig> <fig id="pone-0003007-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003007.g006</object-id><label>Figure 6</label><caption><title>PAC Identification of novel genetic changes.</title><p>(A) PAC detection of novel genetic changes in <italic>EGFR</italic>. PAC predicted skipping of the last four exons of GBM157 and the 5′ end of exon 17 in GBM172. Semi-quantitative PCR on genomic DNA confirmed the deletion in GBM157 (not shown). (B) PAC predicts skipping of exon 30 in the <italic>FCGBP</italic> gene in GBM60. (C) RT-PCR confirmed the <italic>FCGBP</italic> exon skipping event in GBM60; other tumors did not show this exon skipping. (D) Direct sequencing identified a single base change in <italic>EGFR</italic> in GBM172 (as predicted by PAC, see Fig. 6A). (E) Confirmation of a PAC predicted change in the <italic>TLE2</italic> gene in GBM60. The nucleotide substitution overlaps with individual probes of the probe set.</p></caption></fig>
<table-wrap id="pone-0003007-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003007.t001</object-id><label>Table 1</label><caption><title>Detection of exon-skipping mutants by PAC.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Tumor sample</td><td align="left" rowspan="1" colspan="1">Gene</td><td align="left" rowspan="1" colspan="1">Mutation</td><td align="left" rowspan="1" colspan="1">mRNA</td><td align="left" rowspan="1" colspan="1">PAC detection</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">OCUB-F</td><td align="left" rowspan="1" colspan="1">
<italic>CDH1</italic>
</td><td align="left" rowspan="1" colspan="1">c.49_163del</td><td align="left" rowspan="1" colspan="1">r.49_163del115 (Ex2)</td><td align="left" rowspan="1" colspan="1">detected</td></tr><tr><td align="left" rowspan="1" colspan="1">MDA-MB-134VI</td><td align="left" rowspan="1" colspan="1">
<italic>CDH1</italic>
</td><td align="left" rowspan="1" colspan="1">c.688_832del</td><td align="left" rowspan="1" colspan="1">r.688_832del145 (Ex6)</td><td align="left" rowspan="1" colspan="1">detected</td></tr><tr><td align="left" rowspan="1" colspan="1">MPE600</td><td align="left" rowspan="1" colspan="1">
<italic>CDH1</italic>
</td><td align="left" rowspan="1" colspan="1">c.1138_1320del</td><td align="left" rowspan="1" colspan="1">r.1138_1320del183 (Ex9)</td><td align="left" rowspan="1" colspan="1">detected</td></tr><tr><td align="left" rowspan="1" colspan="1">CAMA-1</td><td align="left" rowspan="1" colspan="1">
<italic>CDH1</italic>
</td><td align="left" rowspan="1" colspan="1">c.1712-1G&gt;A</td><td align="left" rowspan="1" colspan="1">r.1566_1712del147 (Ex11)</td><td align="left" rowspan="1" colspan="1">detected</td></tr><tr><td align="left" rowspan="1" colspan="1">MDA-MB-468</td><td align="left" rowspan="1" colspan="1">
<italic>PTEN</italic>
</td><td align="left" rowspan="1" colspan="1">c.253+1G&gt;T</td><td align="left" rowspan="1" colspan="1">r.210_253del44 (Ex4)</td><td align="left" rowspan="1" colspan="1">detected</td></tr><tr><td align="left" rowspan="1" colspan="1">MDA-MB-453</td><td align="left" rowspan="1" colspan="1">
<italic>p53</italic>
</td><td align="left" rowspan="1" colspan="1">c.994_1182del</td><td align="left" rowspan="1" colspan="1">r.994_1182del189 (Ex10-11)</td><td align="left" rowspan="1" colspan="1">detected</td></tr><tr><td align="left" rowspan="1" colspan="1">HCC1937</td><td align="left" rowspan="1" colspan="1">
<italic>RB1</italic>
</td><td align="left" rowspan="1" colspan="1">c.2212_2325del</td><td align="left" rowspan="1" colspan="1">r.2212_2325del114 (Ex22)</td><td align="left" rowspan="1" colspan="1">detected</td></tr><tr><td align="left" rowspan="1" colspan="1">Glioblastoma 67</td><td align="left" rowspan="1" colspan="1">
<italic>EGFR</italic>
</td><td align="left" rowspan="1" colspan="1">c.89_889del</td><td align="left" rowspan="1" colspan="1">r.89_889del801 (Ex2-7)</td><td align="left" rowspan="1" colspan="1">not detected</td></tr><tr><td align="left" rowspan="1" colspan="1">Glioblastoma 96</td><td align="left" rowspan="1" colspan="1">
<italic>EGFR</italic>
</td><td align="left" rowspan="1" colspan="1">c.89_889del</td><td align="left" rowspan="1" colspan="1">r.89_889del801 (Ex2-7)</td><td align="left" rowspan="1" colspan="1">detected</td></tr><tr><td align="left" rowspan="1" colspan="1">Glioblastoma 142</td><td align="left" rowspan="1" colspan="1">
<italic>EGFR</italic>
</td><td align="left" rowspan="1" colspan="1">c.89_889del</td><td align="left" rowspan="1" colspan="1">r.89_889del801 (Ex2-7)</td><td align="left" rowspan="1" colspan="1">detected</td></tr><tr><td align="left" rowspan="1" colspan="1">Glioblastoma 149</td><td align="left" rowspan="1" colspan="1">
<italic>EGFR</italic>
</td><td align="left" rowspan="1" colspan="1">c.89_889del</td><td align="left" rowspan="1" colspan="1">r.89_889del801 (Ex2-7)</td><td align="left" rowspan="1" colspan="1">not detected</td></tr><tr><td align="left" rowspan="1" colspan="1">Glioblastoma 163</td><td align="left" rowspan="1" colspan="1">
<italic>EGFR</italic>
</td><td align="left" rowspan="1" colspan="1">c.89_889del</td><td align="left" rowspan="1" colspan="1">r.89_889del801 (Ex2-7)</td><td align="left" rowspan="1" colspan="1">not detected</td></tr><tr><td align="left" rowspan="1" colspan="1">Glioblastoma 164</td><td align="left" rowspan="1" colspan="1">
<italic>EGFR</italic>
</td><td align="left" rowspan="1" colspan="1">c.89_889del</td><td align="left" rowspan="1" colspan="1">r.89_889del801 (Ex2-7)</td><td align="left" rowspan="1" colspan="1">not detected</td></tr></tbody></table></alternatives></table-wrap>
<disp-formula></disp-formula> <disp-formula></disp-formula>
<supplementary-material content-type="local-data" id="pone.0003007.s001"><label>Table S1</label><caption><p>Complete list of outlier exons at PAC values &lt;−2.0</p><p>(13.64 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0003007.s002"><label>Table S2</label><caption><p>Summary of PCR confirmation results for a randomly selected cohort of outlier exons</p><p>(0.04 MB XLS)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><label/><p>Stated are mutations that result in the expression of an aberrant transcript variant in breast cancer cell lines (top) and clinical glioblastoma samples (bottom). All mutations were determined by sequencing genomic DNA fragments as well as transcripts, where the cDNA sequence c.1 corresponds to the adenosine residue of the ATG initiation codon in the Genbank reference sequence. Genbank accession numbers: Z13009, NM_000314, NM_000546, NM_000321 and NM_005228.3 for <italic>CDH1</italic>, <italic>PTEN</italic>, <italic>p53</italic>, <italic>RB1</italic> and <italic>EGFR</italic>, respectively. Detection of the mutations by PAC was at PAC value −1.3.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>This work was supported by the Susan G. Komen Breast Cancer Foundation (BCTR0601309), the Dutch Cancer Society Koningin Wilhelmina Fonds, DDHK 2002-2687 and Erasmus MC Mrace 2005.</p></fn></fn-group>
<graphic xlink:href="pone.0003007.g001"/> <graphic id="pone-0003007-t001-1" xlink:href="pone.0003007.t001"/> <graphic xlink:href="pone.0003007.g002"/> <graphic xlink:href="pone.0003007.g003"/> <graphic xlink:href="pone.0003007.g004"/> <graphic xlink:href="pone.0003007.g005"/> <graphic xlink:href="pone.0003007.g006"/> <graphic xlink:href="pone.0003007.e001.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0003007.e002.jpg" mimetype="image" position="float"/>
<media xlink:href="pone.0003007.s001.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0003007.s002.xls"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["16"], "element-citation": ["\n"], "surname": ["Su", "Song", "Repasky", "Schutte", "Kern"], "given-names": ["GH", "JJ", "EA", "M", "SE"], "year": ["2002"], "article-title": ["Mutation rate of MAP2K4/MKK4 in breast carcinoma."], "source": ["Hum Mutat"], "volume": ["19"], "fpage": ["81"]}]
PLoS One. 2008 Aug 20; 3(8):e3007