- Research article
- Open access
- Published:
Loss of heterozygosity: what is it good for?
BMC Medical Genomics volume 8, Article number: 45 (2015)
Abstract
Background
Loss of heterozygosity (LOH) is a common genetic event in cancer development, and is known to be involved in the somatic loss of wild-type alleles in many inherited cancer syndromes. The wider involvement of LOH in cancer is assumed to relate to unmasking a somatically mutated tumour suppressor gene through loss of the wild type allele.
Methods
We analysed 86 ovarian carcinomas for mutations in 980 genes selected on the basis of their location in common regions of LOH.
Results
We identified 36 significantly mutated genes, but these could only partly account for the quanta of LOH in the samples. Using our own and TCGA data we then evaluated five possible models to explain the selection for non-random accumulation of LOH in ovarian cancer genomes: 1. Classic two-hit hypothesis: high frequency biallelic genetic inactivation of tumour suppressor genes. 2. Epigenetic two-hit hypothesis: biallelic inactivation through methylation and LOH. 3. Multiple alternate-gene biallelic inactivation: low frequency gene disruption. 4. Haplo-insufficiency: Single copy gene disruption. 5. Modified two-hit hypothesis: reduction to homozygosity of low penetrance germline predisposition alleles. We determined that while high-frequency biallelic gene inactivation under model 1 is rare, regions of LOH (particularly copy-number neutral LOH) are enriched for deleterious mutations and increased promoter methylation, while copy-number loss LOH regions are likely to contain under-expressed genes suggestive of haploinsufficiency. Reduction to homozygosity of cancer predisposition SNPs may also play a minor role.
Conclusion
It is likely that selection for regions of LOH depends on its effect on multiple genes. Selection for copy number neutral LOH may better fit the classic two-hit model whereas selection for copy number loss may be attributed to its effect on multi-gene haploinsufficiency. LOH mapping alone is unlikely to be successful in identifying novel tumour suppressor genes; a combined approach may be more effective.
Background
Cancer cells undergo multiple genetic and epigenetic hits in the development of tumorigenic phenotypes, including somatic point mutations, increases in copy number, gene deletions, gene rearrangements, translocations and promoter hypermethylation [1]. These random events are selected for due to their effect on oncogenes, where the aberration activates the gene to promote tumorigenesis (e.g. KRAS, MYC), and on tumour suppressor genes (TSG), where the genetic or epigenetic aberrations is inactivating (e.g. TP53, PTEN), since the normal function of these genes is to restrict tumorigenic potential.
Loss of heterozygosity (LOH) is a common genetic event in many cancer types, so-called because of the early observations of a change in polymorphic markers from a heterozygous state in the germline to an apparently homozygous state in the tumour DNA [2]. LOH is a general term that encompasses both LOH with copy number losses (CNL-LOH) and copy number neutral LOH (CNN-LOH). In CNL-LOH all or part of a chromosome is deleted. CNN-LOH originates either through a homologous recombination event (“gene conversion”), or because the retained chromosome was duplicated either before or after the LOH event. LOH is strongly associated with loss of the wild-type allele in individuals with an inherited cancer predisposition syndrome and carry a germline mutation in genes such as RB1 in retinoblastoma or BRCA1 in breast and ovarian cancer [2, 3]. This “second hit” hypothesis was initially proposed by Knudson based on his observations of the incidence of familial retinoblastoma [4] and has been widely accepted as a mechanism for the complete inactivation of tumour suppressor genes, both in a germline context and the sporadic cancer context where the first hit is a somatic event, such as mutation of TP53. As a consequence, mapping of common regions of minimal LOH has historically been a popular strategy to pursue the identification of novel TSGs without the need for segregation data from large cancer families. However, such analyses have been generally been unsuccessful leading to speculation that the approach is technically and conceptually flawed [5], and even to whether there is any selective advantage to LOH events. Nonetheless, we previously used SNP mapping arrays to analyse LOH in ovarian carcinomas of diverse histological subtypes, with the rationale that the newer methodology would at least overcome some of the previous technical issues with LOH analyses [6]. We mapped a number of minimal regions of LOH containing tumour suppressor gene candidates, including regions of homozygous deletion encompassing genes such as MAP2K4 [7]. Advances in massively parallel sequencing has enabled the current study where we report targeted sequencing of 980 candidate tumour suppressor genes in 86 ovarian carcinomas, most of which have matched SNP array data enabling the assessment of the importance of LOH in the selection for somatic mutations in ovarian cancer. We evaluated a number of different histological subtypes, since these have different etiologies and causative genes.
Methods
Ethics statement
Accrual and use of patient material for this study was approved by the following Human Research Ethics Committees: Peter MacCallum Cancer Centre Human Research Ethics Committee, Southampton Hospital Human Research Ethics Committee, University of Melbourne Human Research Ethics Committee, Queensland Institute of Medical Research Human Research Ethics Committee, Westmead Hospital Human Research Ethics Committee. All individuals gave written informed consent for the use of their tissue in research. This project was approved by the Peter MacCallum Cancer Centre Human Research Ethics Committee (Approval # 09/29).
Ovarian tumour cohort
A tumour cohort (n = 86) comprising a variety of histological subtypes including serous (n = 45), endometrioid (n = 28), mucinous (n = 7) and clear cell (n = 6) were obtained through the Australia Ovarian Cancer Study, the Peter MacCallum Cancer Centre Tissue Bank, or from patients presenting to hospitals in the south of England [8]. The majority of tumour DNA samples were needle microdissected to ensure greater than 70 % cancer epithelial cell component; other samples were processed from tissue where the reference haematoxylin and eosin stained section showed >70 % tumour epithelial cells. Matching peripheral blood samples were also collected from patients at time of tumour collection and used as a source of germline DNA for somatic mutation detection. Details of the cohort are listed in Additional file 1: Table S1.
Library preparation, target enrichment and sequencing
Library preparation was performed as previously described [9] following the Illumina genomic DNA library preparation protocol (Illumina, San Diego, CA) using an input of 200 ng of tumour or matched normal lymphocyte DNA. Seven custom multiplexing adapters compatible with Illumina single-end sequencing were used and indexed DNA samples were pooled equally prior to PCR enrichment. A boutique exon capture (SureSelect, Agilent Technologies, Santa Clara, CA) was used to enrich for coding exons of candidate tumour suppressor genes (n = 980, Additional file 1: Tables S2 and S3) and known cancer genes (TP53, BRCA1, BRCA2) according to the recommended protocol. Capture probes were designed using default parameters in eArray (Agilent Technologies).
Sequencing of target-enriched DNA libraries were performed using an Illumina GAIIx, generating 75 bp single-end sequence reads. Image analysis and base calling was performed using the Genome Analyser Pipeline v1.5-1.7. Sequence reads were aligned to the human reference genome (GRCh37/hg19 assembly) using BWA [10] and any remaining unmapped reads aligned with Novoalign [11]. The mean coverage for bases within target regions was 70-fold and 92 % had at least 10-fold coverage. This was followed by local realignment with GATK [12]. Point mutations and insertions/deletions (indels) were identified using GATK and Dindel [13] respectively, and annotated according to Ensembl release 56. Sequence variants were called as somatic alterations only when (i) the variant was not called in the matched normal sample or identified as a germline alteration in another tumour/normal pair (ii) the variant was not seen in > =2 independent reads in the matched normal sample following manual inspection of sequence reads using the Integrated Genomics Viewer [14] (iii) the variant was identified in bi-directional sequence reads.
A selection of variants that met the above criteria for a somatic mutation (n = 202) were subjected to validation by conventional PCR amplification and bidirectional capillary electrophoresis on the ABI3130 Genetic Analyser using BigDye Terminator v3.1 sequencing chemistry (Applied Biosystems, Foster City, CA).
SNP arrays and loss of heterozygosity
Affymetrix SNP Mapping array data was obtained for the 86 sequenced cases, 54 by SNP6 arrays (GSE19539, [15]), 26 by 500 K arrays (previously published in [6]), and six previously unreported low-grade endometrioid cases. Affymetrix SNP6 CEL files, HM27 methylation array data (level 3), Agilent expression array data (level 3) and somatic mutation data from 266 tumors generated by The Cancer Genome Atlas (TCGA) were downloaded from the TCGA Data Portal. LOH was detected as described previously in Partek Genomics Suite (Partek, St Louis, MO), using allele-specific copy number that compared the tumour genotype to the matching normal genotype, and evaluated the copy number at heterozygous alleles [6]. The “min” allele had to have a value of <0.5 copies to be called LOH, thus excluding regions of allelic imbalance where at least one copy of both alleles was present. The results published here are in whole or part based upon data generated by The Cancer Genome Atlas pilot project established by the NCI and NHGRI. Information about TCGA and the investigators and institutions who constitute the TCGA research network can be found at http://cancergenome.nih.gov.
Results and discussion
A candidate TSG screen in ovarian cancer – selection of genes from LOH regions
Candidate ovarian tumour suppressor genes (n = 980) were selected for analysis on the basis of their location in frequent regions of LOH or deletion (Additional file 1: Tables S2 and S3) from our previously published SNP array analysis of 122 primary ovarian carcinomas of various histologies [6]. The regions, located on 20 different chromosome arms, met the following three criteria. Firstly, minimal overlapping regions of LOH were included if they were detected in greater than 35 % of all ovarian carcinomas analysed, or secondly, in >35 % of subtype specific minimal overlapping regions of LOH (4 of 9 clear cell carcinomas, 5 of 12 low-grade endometrioid carcinomas (grades 1 and 2), 6 of 16 mucinous carcinomas and 23 of 64 high-grade serous/endometrioid carcinomas (grades 2 and 3 for the serous subtype, grade 3 for endometrioid carcinomas)). Finally, all homozygous deletions within frequent regions of LOH along with the overlapping portion of all recurrent homozygous deletions were included. This gene list included genes with well established roles in cancer such as CDKN2A and PTEN, but for the purposes of this analysis they were included in the “candidate” LOH genes. In addition, the known ovarian cancer genes TP53, BRCA1 and BRCA2 were included despite lying outside the minimal regions of LOH.
A candidate TSG screen in ovarian cancer – correlation of mutations with LOH
A targeted mutation screen was conducted on the 86 ovarian cancer cases including high-grade serous and endometrioid, low-grade endometrioid, clear cell and mucinous subtypes. Somatic coding mutations were detected in both candidate (561 variants in 366 genes, Additional file 2: Table S4) and known cancer genes (58 TP53 mutations in 56 cases and two mutations in BRCA1). Eighty-nine genes had two or more non-synonymous mutations. The classic two-hit hypothesis would predict that driver genes have homozygous, deleterious mutations in samples with LOH. With respect to deleterious mutation status this was certainly true for TP53 and BRCA1 where a high proportion of somatic mutations were truncating (25/58 and 2/2, respectively) compared to an overall truncating mutation frequency of 13 % (72/561). In addition, among the 53 cases with TP53 or BRCA1 somatic mutations where SNP data was available, 50 (94 %) showed LOH of the wild-type allele. This was is sharp contrast with the other candidate genes, where only 181/520 showed LOH of the wild-type allele (35 %); in particular, there was no significant difference in LOH of the wild-type allele between non-synonymous mutations (134/381 with LOH, 35.2 %) and synonymous mutations (47/139 with LOH, 33.8 %). The overall frequency of non-synonymous compared to synonymous mutations was 73 % (411/561) for the candidate TSGs, but 100 % of mutations in known cancer genes were non-synonymous (60/60). This difference in ratio suggests that the majority of mutations in candidate TSGs from LOH regions are likely to be passenger events, since this rate might be expected without any strong positive selection [16]. The lack of difference in LOH between synonymous and non-synonymous also implies that there is limited selection for homozygosity for the majority of gene mutations.
Significance analysis of recurrently mutated gene candidates
Within the list of mutated genes, we applied a number of filters to assess whether any genes could function as tumor suppressors under either a one-hit or two-hit mechanism. Firstly, significantly mutated genes were identified using the MuSiC algorithm [17], which determines the significance of the observed mutation rate of each gene based on the background mutation rate in the sample cohort. Three known ovarian cancer genes (TP53, PTEN and CDKN2A) were identified by all three tests (convolution, likelihood ratio and Fisher’s combined p-value tests) with a false discovery rate (FDR) of less than 0.10. At this FDR the genes DNAH9, LINGO1, MEF2C, SAMD11, STARD5, ZNRF4 and ZNF287 were also identified, although each was supported only by the likelihood ratio test.
Secondly, the 125 genes with recurrent mutations were assessed for the proportion of cases with biallelic mutation, including by homozygous deletion from SNP array data. For genes with mutations in three or more cases, the proportion of biallelic mutations was greater than 80 % of mutated samples for eight genes (AL355987.1, CASK, CDKN2A, MAP2K4, NF1, PTEN, RB1 and TP53) and between 60-80 % for seven genes (FANCA, GRAMD4, GPR98, IL16, MYOCD, SYNE1 and TEX15). Biallelic mutations were detected in 2/2 mutated samples in an additional nine genes (SKG223, APOOL, BRCA1, CDH8, DACH2, EPHX2, FARP1, PNMA3 and RAI1).
Finally, genes recurrently targeted by inactivating mutations were identified. Mutations with overtly deleterious consequences were considered for this analysis, including nonsense and essential splice site mutations, frameshift indels and gene deletions. Although missense amino acid changes and in-frame indels can also negatively impact gene function, interpreting these mutations in the absence of functional validation is challenging. Sixteen genes were identified where more than half of their mutations would be considered clearly deleterious, including seven known ovarian cancer genes (PTEN, CDKN2A, MAP2K4, PIK3R1, RB1, FANCA and BRCA1). These three analyses identified 36 genes as possible tumour suppressors (Table 1), and it was notable that seven well characterised tumour suppressors were identified by at least two of the three methods (BRCA1, TP53, RB1, PTEN, CDKN2A, FANCA, and MAP2K4), although others were only identified by one method (NF1, PIK3R1). In contrast, 22 of the 27 (81 %) novel/less well characterized genes were identified by only 1 method, indicating that regions of LOH are not strongly enriching for novel genes with classic tumor suppressor gene characteristics.
Loss of heterozygosity – what is it good for?
From the data above it appears that we did not identify dominant, very frequently mutated novel genes where selection for a classic two-hit tumour suppressor gene was apparent. So what, if anything, is the LOH for? We considered five possibilities (Fig. 1) and assessed each in turn.
Models of LOH. Boxes = genes; “X” = inactivating mutation; A, B = alternative alleles of a single nucleotide polymorphism. In the top panels, the black line on the graph represents the overall frequency of LOH observed in tumour samples across the chromosome, while the red bars are the frequency of mutation in a particular gene. Thus, for the classic two-hit model, the frequency of mutation is similar to the frequency of LOH, while in the low frequency model, the frequency of LOH is higher than the mutation rate, because each sample is mutated in a different gene. In the bar graphs below, at left, the red bars represent the frequency of the A allele that is retained in samples with LOH at the locus; thus, the risk locus (*) has a higher proportion of the risk allele (A) retained after LOH compared to a non-risk locus, where the A and B alleles are equally retained. At right, the graphs represents the average reduction in expression of a gene in samples with LOH, compared to samples without LOH; genes in LOH regions show a reduction in expression
-
1.
Classic two-hit hypothesis: high frequency biallelic genetic inactivation of TSG
-
2.
Epigenetic two-hit hypothesis: biallelic inactivation through methylation and LOH
-
3.
Multiple alternate gene biallelic inactivation: low frequency gene disruption
-
4.
Haplo-insufficiency: Single copy gene disruption
-
5.
Modified two-hit hypothesis: reduction to homozygosity of predisposition alleles
-
1.
Classic two-hit hypothesis: high frequency biallelic genetic inactivation of TSG
This mechanism is demonstrably true for many known tumour suppressor genes, with TP53 being a clear example of a gene functioning as a classical TSG in ovarian cancer [18, 19]. However, from our data and large published studies such as TCGA, it is clear that novel genes with a high frequency of biallelic mutations are exceedingly rare and can not explain the bulk of the observed LOH. For example, 8p undergoes LOH in >40 % of ovarian carcinomas, but no gene in this region is mutated at frequency higher than 3 % in our or any other study, although homozygous deletion can target, for example, CSMD1 in 11 % of cases [20]. It remains a possibility, however, that genes not represented on our targeted or exome sequencing platforms could still be the target of such LOH, for example long non-coding RNAs.
-
2.
Epigenetic two-hit hypothesis: high frequency biallelic inactivation through methylation and LOH
Somatic gene mutation is not the only mechanism of biallelic inactivation. Some TSGs can be inactivated through a combination of LOH and promoter hypermethylation, for example MLH1. This methylation can be acquired somatically or may be a consequence of imprinting. We assessed this possibility using TCGA ovarian cancer methylation data. Globally, we observed that there was no enrichment for methylation in regions of LOH – in samples with LOH at a locus, on average 12.7 % of genes were strongly methylated (probe value of >0.75), whereas 13.65 % of genes were strongly methylated when there was no LOH, (Fig. 2a). CNL-LOH was less likely to have strongly methylated genes than CNN-LOH (12.3 % vs 12.9 %, p < 0.0001, Chi-squared test). However, when we analysed the X chromosome separately, we found that samples with any LOH were more likely to have low methylation levels (45.3 % of genes had a probe value of <0.25, compared to 35.5 % in samples without LOH, p < 0.0001 Chi-squared test).
Methylation. a Frequency of probes on the HM27 methylation array that have high (value > 0.75), intermediate (0.25-0.75) and low (<0.25) methylation associated with LOH in a sample, comparing all autosomes (no difference between LOH and no LOH) and chromosome X (more low methylation probes in LOH regions). b Considering only probes that were significantly different between LOH and no LOH, the frequency of significant probes where the mean methylation ratio (LOH/no LOH) was increased (>1.5, higher in LOH) or decreased (<0.75, higher in no LOH). c Percentage of significant methylation probes that are located in regions where the majority (>2/3) of LOH is either copy number loss, neutral or neither. Only regions with at least 20 % LOH were included
Detection of methylation is challenging from both technical and biological perspectives. Tumour and cell type heterogeneity may influence the degree of methylation detected, so we also took an alternative approach where we used the methylation array data to test whether there were genes that were more strongly methylated in samples with LOH compared to samples without LOH. Using a multiple testing correction p-value threshold of 2.2x10−6, there were 1584/22374 (7 %) methylation probes that were significantly differentially methylated. Interestingly, 28 % of these significant probes were located on the X chromosome and indeed 51 % of all probes on the X chromosome were significantly differentially methylated, with lower average levels of methylation in samples with LOH compared to samples without LOH. On the autosomes, the outcome was reversed: 50.3 % of the statistically significant probes had a fold-change difference in mean methylation of >1.5, while only 1.1 % had a fold-change difference of <0.75 (Fig. 2b). Thus, for the X chromosome it appears there is selection for retaining the active copy, perhaps because loss of this copy would be cell lethal as an effective homozygous inactivation of the chromosome. In contrast, for the autosomes there appears to be selection for increased methylation by LOH.
We also evaluated whether there was any difference by the type of LOH. For those genes occurring in a region of LOH with at least 20 % frequency, we determined whether a probe was in a CNL-LOH enriched locus (>66 % of samples with LOH also had CN loss) or a CNN-LOH enriched locus (>66 % of samples with LOH were CNN). Of the CNL-enriched probes, 11.3 % were significantly differentially methylated, compared to 21.7 % of CNN-enriched probes (p < 0.0001, Chi-squared test, Fig. 2c). This data would support a model whereby differential methylation is more commonly selected for in regions of CNN-LOH than CNL-LOH.
-
3.
Multiple alternate gene biallelic inactivation: low frequency gene disruption
Another possibility is that particular loci harbour multiple TSGs but individual tumours only require one to be inactivated and the gene targeted can differ from tumour to tumour. If this is the case then locating the TSGs by mapping overlapping regions of LOH would incorrectly flag the interval between two TSG as the likely location of the TSG – in effect then the peak LOH regions may not be the most likely places to find the targeted gene(s). To evaluate this possibility, we used TCGA data to see whether regions of LOH were enriched for somatic mutations on a sample-by-sample basis. Cases with both somatic exome and SNP array data were used (n = 266). There were 13,148 coding somatic mutations, of which 29.7 % were located within a region of LOH in the sample where it was observed. The average overlap of all the genes assayed with regions of LOH per sample was 35.5 %. Thus, somatic mutations are if anything under-represented in regions of LOH (Binomial test p < 0.0001). Given that most of these mutations are likely to be passengers, we evaluated whether this was true for non-synonymous or overtly deleterious mutations (nonsense, frameshift, essential splice site). For deleterious mutations, 38.2 % were in regions of LOH (p = 0.035, Binomial test), whereas only 25.2 % of the non-synonymous mutations were in regions of LOH, similar to the 22.2 % observed for synonymous mutations. The signal for deleterious mutations was substantially reduced if TP53 was excluded (34.9 % of deleterious and 28.5 % of other non-synonymous mutations had LOH). Thus, regions of LOH are slightly enriched for deleterious mutations, but not for other non-synonymous or silent mutations.
We then evaluated whether there was a difference in mutation frequency in CNN versus CNL regions of LOH on a case by case basis as above (Fig. 3). Excluding TP53, there were fewer mutations in regions of CNL-LOH than would have been expected based on the overall percentage of the exome affected (19.8 % of mutations were in CNL-LOH regions, whereas 26.8 % of the exome was affected by CNL-LOH, p < 0.0001, Binomial test). The difference was less striking when considering overtly deleterious mutations only (24.6 % vs 26.8 %, p = 0.09, Binomial test). For CNN-LOH, the overall difference was small (8.8 % of mutations vs 8.7 % of the exome affected by CNN-LOH, p = 0.7, Binomial test), however there were more deleterious mutations than expected (10.3 % vs 8.7 %, p = 0.05, Binomial test). When TP53 was included, both total mutations (9.2 %) and deleterious mutations (11.4 %) in CNN-LOH regions were increased. Silent mutations were the most likely to be underrepresented in CNL-LOH regions (17.2 % vs 26.8 %). It is possible, therefore, that mutations are seen less often in CNL-LOH regions simply as a consequence of decreased DNA dosage. The enrichment of deleterious mutations in CNN-LOH regions, however, suggests the presence of positive selection for mutations in TSGs.
Mutation load. Frequency of mutations of various types in copy number neutral and copy number loss regions, compared to the overall frequency of LOH across the exome (“LOH overall”). Deleterious mutations are enriched in CNN-LOH; all other mutations types are less frequent in copy number loss regions
-
4.
Haplo-insufficiency: Single copy gene disruption
We and others have shown that loss of a single gene copy can reduce gene expression [6, 21]. A recent study showed that regions of copy number loss are enriched for tumour-suppressor genes [22], but that each gene might have a limited effect on its own. Chromosome complementation studies, where all or part of a chromosome is introduced into cell lines with LOH of that chromosome via microcell-mediated monochromosome transfer, have frequently been able to show reduction in tumorigenicity of the cell line thus complemented [23, 24], but have only rarely been able to implicate a single gene responsible [25, 26]. Thus, haplo-insufficiency of multiple genes, each with a small effect, could contribute to the non-random pattern of LOH observed in ovarian cancer, especially for chromosomal regions that are weighted towards CNL-LOH such as 8p and X, rather than CNN-LOH, such as 17.
We previously observed a correlation between the percentage of genes under-expressed and the percentage of cases with CNL-LOH, as opposed to CNN-LOH, in a region-wise comparison of LOH vs. no LOH [6]. In an analysis of TCGA data, we compared the expression of 10,925 genes between cases with and without LOH. 3,780 genes (34.6 %) were significantly differentially expressed (at a multiple testing p-value threshold of 4.56x10−6), and all significant genes were under-expressed in samples with LOH compared to samples without LOH. When comparing CNN-LOH to CNL-LOH in genes with at least 20 % frequency of LOH, only 1/163 genes (0.6 %) at CNN-LOH enriched loci were significantly differentially expressed, compared to 2701/5740 (47 %) CNL-LOH enriched genes (p < 0.0001, Chi-squared test, Fig. 4). This result supports the idea that chromosomal regions with CNL-LOH may contain genes where loss of a single copy results in reduced gene expression and a selective advantage to the cell. In contrast, chromosomal regions with little copy number loss may contain essential genes for which haplo-insufficiency is cell lethal.
-
5.
Modified two-hit hypothesis: reduction to homozygosity of predisposition alleles
In familial cancer predisposition syndromes, it is common for the remaining wild-type allele to be lost by LOH, for example BRCA1 pathogenic variants are usually reduced to homozygosity in breast and ovarian carcinomas [3, 27]. However, common low-penetrance risk alleles could also be targeted by LOH leading to an enhancement of their cancer-promoting role. We assessed this using nine SNP loci identified in the iCOGs study [28, 29] as predisposing to all ovarian cancer types or high-grade serous ovarian cancer. Two of these SNPs were present on the Affymetrix SNP 6.0 array, the remainder were represented by SNPs in linkage disequilibrium (r2 > 0.7, from HapMap [30]). Where possible, up to four linked SNPs were evaluated.
For each SNP, we assessed whether cases were heterozygous in their normal DNA, and what proportion of these with LOH of the region were homozygous for the risk allele in the tumour DNA using TCGA and our own data (n = 364). Interestingly, two SNPs at 10p12 linked to the risk allele rs1232180 were found to be significantly more likely to have lost the non-risk allele than the risk allele (Table 2). Some other SNPs linked to a risk allele also showed significantly non-random loss of the non-risk allele, but the data were not consistent across all SNPs examined at the locus (e.g. 17q21, 3q25 and 9p22). It is not clear whether these discrepancies could be due to technical variation in the SNP calling; alternative methods may be required to assess this possibility. The remainder of SNPs were not significant, however several are uncommon, limiting the power of the analysis. Thus, it is possible that some LOH may be selected for through the phenotypic effect of reduction to homozygosity of predisposition alleles.
Conclusion
The broader relevance of LOH in cancer has been debated for some time [5, 31] although many of the criticisms stemmed from technical issues that are being overcome by newer methodologies. Our initial assumption for this study was that we would detect high-frequency mutated genes in the minimal peak regions of LOH we had defined by LOH mapping using these newer methodologies; i.e. a classic two-hit model. However, the biology of LOH does not support this assumption and with large-scale tumour studies it is now possible to explore the many possibilities for the functional significance of this genetic event as summarised in Table 3. We suggest that the non-random patterns of LOH detected in cancer are a result of multiple different mechanisms operating to affect multiple genes, which may differ from tumour to tumour yet collectively play a role in the development of the tumorigenic phenotype. It is worth noting the differences in CNL-LOH versus CNN- LOH, with the latter appearing more relevant for selection of deleterious mutations and methylation, in contrast to global changes in gene expression. Identifying the specific driver genes targeted in a particular cancer remains a challenge given the multiple possible reasons for selection of an LOH event.
Abbreviations
- AOCS:
-
Australian Ovarian Cancer Study
- LOH:
-
Loss of heterozygosity
- CNL-LOH:
-
Copy number loss LOH
- CNN-LOH:
-
Copy number neutral LOH
- SNP:
-
Single nucleotide polymorphism
- TCGA:
-
The Cancer Genome Atlas
- TSG:
-
Tumour suppressor gene
References
Vogelstein B, Papadopoulos N, Velculescu VE, Zhou S, Diaz Jr LA, Kinzler KW. Cancer genome landscapes. Science. 2013;339(6127):1546–58. doi:10.1126/science.1235122.
Cavenee WK, Dryja TP, Phillips RA, Benedict WF, Godbout R, Gallie BL, et al. Expression of recessive alleles by chromosomal mechanisms in retinoblastoma. Nature. 1983;305(5937):779–84.
Merajver SD, Frank TS, Xu J, Pham TM, Calzone KA, Bennett-Baker P, et al. Germline BRCA1 mutations and loss of the wild-type allele in tumors from families with early onset breast and ovarian cancer. Clin Cancer Res. 1995;1(5):539–44.
Knudson Jr AG. Mutation and cancer: statistical study of retinoblastoma. Proc Natl Acad Sci U S A. 1971;68(4):820–3.
Tomlinson IP, Lambros MB, Roylance RR. Loss of heterozygosity analysis: practically and conceptually flawed? Genes Chromosomes Cancer. 2002;34(4):349–53. doi:10.1002/gcc.10085.
Gorringe KL, Ramakrishna M, Williams LH, Sridhar A, Boyle SE, Bearfoot JL, et al. Are there any more ovarian tumor suppressor genes? A new perspective using ultra high-resolution copy number and loss of heterozygosity analysis. Genes Chromosomes Cancer. 2009;48(10):931–42.
Davis SJ, Choong DY, Ramakrishna M, Ryland GL, Campbell IG, Gorringe KL. Analysis of the mitogen-activated protein kinase kinase 4 (MAP2K4) tumor suppressor gene in ovarian cancer. BMC Cancer. 2011;11:173. doi:10.1186/1471-2407-11-173.
Bryan EJ, Watson RH, Davis M, Hitchcock A, Foulkes WD, Campbell IG. Localization of an ovarian cancer tumor suppressor gene to a 0.5-cM region between D22S284 and CYP2D, on chromosome 22q. Cancer Res. 1996;56(4):719–21.
Ryland GL, Bearfoot JL, Doyle MA, Boyle SE, Choong DY, Rowley SM, et al. MicroRNA genes and their target 3'-untranslated regions are infrequently somatically mutated in ovarian cancers. PLoS One. 2012;7(4):e35805. doi:10.1371/journal.pone.0035805.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14):1754–60. doi:10.1093/bioinformatics/btp324.
Novocraft Technologies. http://www.novocraft.com.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297-303. doi:10.1101/gr.107524.110
Albers CA, Lunter G, Macarthur DG, McVean G, Ouwehand WH, Durbin R. Dindel: Accurate indel calls from short-read data. Genome Res. 2010. doi:10.1101/gr.112326.110
Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G et al. Integrative genomics viewer. Nat Biotechnol. 2011;29(1):24-6. doi:10.1038/nbt.1754
Ramakrishna M, Williams LH, Boyle SE, Bearfoot JL, Sridhar A, Speed TP, et al. Identification of candidate growth promoting genes in ovarian cancer through integrated copy number and expression analysis. PLoS One. 2010;5(4), e9983. doi:10.1371/journal.pone.0009983.
Pleasance ED, Cheetham RK, Stephens PJ, McBride DJ, Humphray SJ, Greenman CD, et al. A comprehensive catalogue of somatic mutations from a human cancer genome. Nature. 2010;463(7278):191–6. doi:10.1038/nature08658.
Dees ND, Zhang Q, Kandoth C, Wendl MC, Schierding W, Koboldt DC, et al. MuSiC: identifying mutational significance in cancer genomes. Genome Res. 2012;22(8):1589–98. doi:10.1101/gr.134635.111.
Ahmed AA, Etemadmoghadam D, Temple J, Lynch AG, Riad M, Sharma R, et al. Driver mutations in TP53 are ubiquitous in high grade serous carcinoma of the ovary. J Pathol. 2010;221(1):49–56. doi:10.1002/path.2696.
McManus DT, Murphy M, Arthur K, Hamilton PW, Russell SE, Toner PG. p53 mutation, allele loss on chromosome 17p, and DNA content in ovarian carcinoma. J Pathol. 1996;179(2):177–82. doi:10.1002/(SICI)1096-9896(199606)179:2<177::AID-PATH561>3.0.CO;2-2.
Cancer Genome Atlas Research N. Integrated genomic analyses of ovarian carcinoma. Nature. 2011;474(7353):609–15. doi:10.1038/nature10166.
Pollack JR, Sorlie T, Perou CM, Rees CA, Jeffrey SS, Lonning PE, et al. Microarray analysis reveals a major direct role of DNA copy number alteration in the transcriptional program of human breast tumors. Proc Natl Acad Sci U S A. 2002;99(20):12963–8.
Davoli T, Xu AW, Mengwasser KE, Sack LM, Yoon JC, Park PJ, et al. Cumulative haploinsufficiency and triplosensitivity drive aneuploidy patterns and shape the cancer genome. Cell. 2013;155(4):948–62. doi:10.1016/j.cell.2013.10.011.
Cody NA, Ouellet V, Manderson EN, Quinn MC, Filali-Mouhim A, Tellis P, et al. Transfer of chromosome 3 fragments suppresses tumorigenicity of an ovarian cancer cell line monoallelic for chromosome 3p. Oncogene. 2007;26(4):618–32. doi:10.1038/sj.onc.1209821.
Stronach EA, Sellar GC, Blenkiron C, Rabiasz GJ, Taylor KJ, Miller EP, et al. Identification of clinically relevant genes on chromosome 11 in a functional model of ovarian cancer tumor suppression. Cancer Res. 2003;63(24):8648–55.
England NL, Cuthbert AP, Trott DA, Jezzard S, Nobori T, Carson DA, et al. Identification of human tumour suppressor genes by monochromosome transfer: rapid growth-arrest response mapped to 9p21 is mediated solely by the cyclin-D-dependent kinase inhibitor gene, CDKN2A (p16INK4A). Carcinogenesis. 1996;17(8):1567–75.
Murakami Y. Functional cloning of a tumor suppressor gene, TSLC1, in human non-small cell lung cancer. Oncogene. 2002;21(45):6936–48. doi:10.1038/sj.onc.1205825.
Berchuck A, Heron KA, Carney ME, Lancaster JM, Fraser EG, Vinson VL, et al. Frequency of germline and somatic BRCA1 mutations in ovarian cancer. Clin Cancer Res. 1998;4(10):2433–7.
Earp MA, Kelemen LE, Magliocco AM, Swenerton KD, Chenevix-Trench G, Australian Cancer S, et al. Genome-wide association study of subtype-specific epithelial ovarian cancer risk alleles using pooled DNA. Hum Genet. 2014;133(5):481–97. doi:10.1007/s00439-013-1383-3.
Pharoah PD, Tsai YY, Ramus SJ, Phelan CM, Goode EL, Lawrenson K, et al. GWAS meta-analysis and replication identifies three new susceptibility loci for ovarian cancer. Nat Genet. 2013;45(4):362–70. doi:10.1038/ng.2564.70e1-2.
International HapMap C. The International HapMap Project. Nature. 2003;426(6968):789–96. doi:10.1038/nature02168.
Presneau N, Manderson EN, Tonin PN. The quest for a tumor suppressor gene phenotype. Curr Mol Med. 2003;3(7):605–29.
Futreal PA, Coin L, Marshall M, Down T, Hubbard T, Wooster R, et al. A census of human cancer genes. Nat Rev Cancer. 2004;4(3):177–83. doi:10.1038/nrc1299.
The Cancer Genome Atlas Research Network. Integrated genomic analyses of ovarian carcinoma. Nature. 2011;474(7353):609-15. doi:10.1038/nature10166
Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA et al. The cBio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics Data. Cancer Discov. 2012;2(5):401-4. doi:10.1158/2159-8290.CD-12-0095
Acknowledgements
We gratefully acknowledge the cooperation of the institutions in Australia participating in the Australian Ovarian Cancer Study (AOCS). We also acknowledge the contribution of the AOCS study nurses, research assistants and all clinical and scientific collaborators and would like to thank all of the women who participated in AOCS. The AOCS management group comprises David Bowtell (Peter MacCallum Cancer Centre, East Melbourne, Victoria, Australia), Georgia Chenevix-Trench (Queensland Institute of Medical Research, Brisbane, Queensland, Australia), Adele Green (Queensland Institute of Medical Research, Brisbane, Queensland, Australia), Penny Webb (Queensland Institute of Medical Research, Brisbane, Queensland, Australia), Anna deFazio (Westmead lnstitute for Cancer Research, Westmead Millennium Institute, Westmead, New South Wales, Australia ), Dorota Gertig (Victorian Cervical Cytology Registry, Carlton South, Victoria, Australia. The scientific and clinical collaborators are: ACT- R Stuart-Harris; NSW- F Kirsten, J Rutovitz, P Clingan, A Glasgow, A Proietto, S Braye, G Otton, J Shannon, T Bonaventura, J Stewart, S Begbie, M Friedlander, D Bell, S Baron-Hay, A Ferrier (dec.), G Gard, D Nevell, N Pavlakis, S Valmadre, B Young, C Camaris, R Crouch, L Edwards, N Hacker, D Marsden, G Robertson, P Beale, J Beith, J Carter, C Dalrymple, R Houghton, P Russell, L Anderson, M Links, J Grygiel, J Hill, A Brand, K Byth, R Jaworski, P Harnett, R Sharma, G Wain; QLD- D Purdie, D Whiteman, B Ward, D Papadimos, A Crandon, M Cummings, K Horwood. A Obermair, L Perrin, D Wyld, J Nicklin; SA- M Davy, MK Oehler, C Hall, T Dodd, T Healy, K Pittman, D Henderson, J Miller, J Pierdes, A Achan; TAS- P Blomfield, D Challis, R McIntosh, A Parker; VIC- B Brown, R Rome, D Allen, P Grant, S Hyde, R Laurie M Robbie, D Healy, T Jobling, T Manolitsas, J McNealage, P Rogers, B Susil, E Sumithran, I Simpson, I Haviv, K Phillips, D Rischin, S Fox, D Johnson, S Lade, P Waring, M Loughrey, N O’Callaghan, B Murray, L Mileshkin, P Allan; V Billson, J Pyman, D Neesham, M Quinn, A Hamilton, O McNally, C Underhill, R Bell, LF Ng, R Blum, V Ganju; WA- I Hammond, A McCartney (dec.), C Stewart, Y Leung, M Buck, N Zeps (WARTN). Further information can be found at http://www.aocstudy.org. This work was supported by the National Health and Medical Research Council of Australia (NHMRC ID #628773) and the Emer Casey Foundation.
Author information
Authors and Affiliations
Consortia
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Author contributions
GLR conducted SNP array, mutation validation and all massively parallel sequencing experiments and analysed sequencing data; MAD and JL performed sequencing alignments and variant calling; DG performed bioinformatics analyses; SEB, DYHC and SMR extracted DNA and undertook Sanger sequencing validation; AOCS and DDLB contributed samples and data; RWT contributed to assay design and optimisation; IGC conceived of the study and participated in its design and coordination; and KLG conceived of the study, participated in its design and coordination, undertook bioinformatics analyses and drafted the manuscript. All authors read and approval the final manuscript.
Author’s information
Ian G. Campbell and Kylie L. Gorringe, co-senior authors.
Additional files
Additional file 1: Table S1.
Cohort information, Table S2. Regions of LOH analysed by targeted sequencing, Table S3. List of the genes in regions of LOH analysed by targeted sequencing. (DOCX 135 kb)
Additional file 2: Table S4.
Excel file with all synonymous and non-synonymous mutations identified by targeted sequencing. (XLSX 141 kb)
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.
The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Ryland, G.L., Doyle, M.A., Goode, D. et al. Loss of heterozygosity: what is it good for?. BMC Med Genomics 8, 45 (2015). https://doi.org/10.1186/s12920-015-0123-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12920-015-0123-z