- Research article
- Open Access
Massive interstitial copy-neutral loss-of-heterozygosity as evidence for cancer being a disease of the DNA-damage response
BMC Medical Genomics volume 8, Article number: 42 (2015)
The presence of loss-of-heterozygosity (LOH) mutations in cancer cell genomes is commonly encountered. Moreover, the occurrences of LOHs in tumor suppressor genes play important roles in oncogenesis. However, because the causative mechanisms underlying LOH mutations in cancer cells yet remain to be elucidated, enquiry into the nature of these mechanisms based on a comprehensive examination of the characteristics of LOHs in multiple types of cancers has become a necessity.
We performed next-generation sequencing on inter-Alu sequences of five different types of solid tumors and acute myeloid leukemias, employing the AluScan platform which entailed amplification of such sequences using multiple PCR primers based on the consensus sequences of Alu elements; as well as the whole genome sequences of a lung-to-liver metastatic cancer and a primary liver cancer. Paired-end sequencing reads were aligned to the reference human genome to identify major and minor alleles so that the partition of LOH products between homozygous-major vs. homozygous-minor alleles could be determined at single-base resolution. Strict filtering conditions were employed to avoid false positives. Measurements of LOH occurrences in copy number variation (CNV)-neutral regions were obtained through removal of CNV-associated LOHs.
We found: (a) average occurrence of copy-neutral LOHs amounting to 6.9 % of heterologous loci in the various cancers; (b) the mainly interstitial nature of the LOHs; and (c) preference for formation of homozygous-major over homozygous-minor, and transitional over transversional, LOHs.
The characteristics of the cancer LOHs, observed in both AluScan and whole genome sequencings, point to the formation of LOHs through repair of double-strand breaks by interhomolog recombination, or gene conversion, as the consequence of a defective DNA-damage response, leading to a unified mechanism for generating the mutations required for oncogenesis as well as the progression of cancer cells.
As a common feature of cancer cells, LOHs have been investigated by cytogenetics, fluorescence in situ hybridization, comparative genomic hybridization (CGH), array-CGH, and single nucleotide polymorphism (SNP)-based microarrays [1–4]. With the application of next-generation sequencing, analysis of LOH in cancer can further be conducted at the level of single base resolution [5, 6]. However, owing to the importance of LOHs giving rise to loss of major alleles and inactivation of tumor suppressor genes, hitherto investigations of LOHs in cancers have been focused mainly on LOHs that yield homozygous-minor genotypes. Yet a comprehensive understanding of the properties and origins of LOHs in cancers requires analysis of all types of LOHs in multiple cancers.
Accordingly, in the present study next-generation sequencing was applied to determine at single-base resolution the LOHs in the genomic sequences of various types of cancers, covering not only sequence regions that have undergone loss of heterozygosity but also single nucleotide changes where a heterozygous position has mutated to a homozygous one. Thirty tumor-control pairs of six different types of cancers including glioma (glioblastoma and astroglioma), acute myeloid leukemia, gastric adenocarcinoma, hepatocellular carcinoma, primary lung cancer (pulmonary squamous-cell carcinoma, adenocarcinoma and neuroendocrinal carcinoma), and lung-to-brain metastatic adenocarcinoma were analyzed with the AluScan platform established by our laboratory, based on the capture of ~8–25 Mb/genome of inter-Alu sequences by inter-Alu PCR amplification using multiple consensual Alu sequence-based primers for next-generation sequencing .
The results obtained on both LOH mutations forming homozygous-major genotypes and those forming homozygous-minor genotypes have yielded a comprehensive LOH landscape across different types of cancers that identifies the outstanding characteristics of cancer LOHs: (a) occurrence of massive percentile LOH mutations of heterozygous residues in the cancer genomes, far exceeding the percentile gain-of-heterozygosity (GOH) mutations of homozygous-major residues; (b) cancer LOHs are mainly interstitial ones indicative of gene conversion rather than segmental deletion as the major underlying mechanism for their production; and (c) cancer LOHs display preferences for the production of homozygous-major genotypes over homozygous-minor genotypes, and for transitional over transversional changes. These characteristics of cancer LOHs, determined using the AluScan platform and also confirmed by the whole-genome sequences reported for a lung-to-liver metastatic cancer  and a primary liver cancer , indicate that cancer LOHs are generated mainly by repair of double-strand breaks (DSB) through interhomolog recombination with the homologous chromosome serving as repair template. The massive scale of the interhomolog recombinations called for by the cancer LOHs suggests that a defective DNA-damage response, by weakening cell cycle checkpoints in the cancer cells, allows the entry of DSB-bearing DNA into the S-phase of the cell cycle, thereby enabling interhomolog recombination and production of the LOH and tag-along GOH mutations needed by the cells during their post-oncogenesis as well as pre-oncogenesis phases.
Participation in this study was voluntary and informed consent was obtained from each of the Han Chinese patients. Institutional Ethics Committees approvals for this study were granted by Hong Kong University of Science and Technology, Second Military Medical University of Shanghai, The First Hospital of Nanjing, Jiangsu Cancer Hospital, Chinese University of Hong Kong, and Capital Medical University of Beijing. Of the thirty patients, five were diagnosed for adenocarcinoma of the stomach; five for glioma including three with glioblastoma and two with astroglioma; five for acute myeloid leukemia; five for primary hepatocellular carcinoma; five for lung cancer including two with pulmonary squamous-cell carcinoma, two with pulmonary adenocarcinoma and one with pulmonary neuroendocrinal carcinoma; and five for lung-to-brain metastatic adenocarcinoma. Detailed medical records of the patients are given in Additional file 1: Table S2 footnotes.
DNA samples from normal white blood cells and leukemia cells were prepared using phenol-chloroform extraction; and DNA samples from normal lung tissue and solid tumors were prepared using DNAzol® Reagent from Life Technologies.
Normal white blood cells were employed as controls for the tumor samples, except for Lung 4 and Lung 5 in Table 1 where normal lung tissue served as controls. In the case of leukemias, normal white blood cells were separated from leukemia cells by Ficoll density gradient centrifugation. Sequence data on the whole genome sequences of a lung-to-liver metastatic lung adenocarcinoma and its blood cell control , and a primary hepatitis B positive hepatocellular carcinoma and its normal liver tissue control (case number DD59) , were obtained online from EBI-SRA at www.ebi.ac.uk/ena/data/view/ERP001071 and NCBI Short Reads Archive (accession number SRA076160) respectively.
Inter-Alu PCR and next-generation sequencing
Next-generation sequencing technologies have transformed genetics through their ability to produce giga-bases of sequence information in a single run. However, next-generation sequencing of a subset of the genome captured by inter-Alu PCR with an amplicon range vastly enhanced by the use of both ‘head type’ and ‘tail type’ Alu consensus sequence-based PCR primers of opposing orientations could substantially reduce the amount of sample DNA required as well as data analysis . Such sets of next-generation sequenced inter-Alu PCR amplicons, or AluScans, were employed in this study to provide an expedited scan of the mutations in exons, introns and non-coding regions. For this purpose, a 25-μl PCR reaction mixture contained 2 μl Bioline 10× NH4 buffer (160 mM ammonium sulfate, 670 mM Tris–HCl, pH 8.8, 0.1 % stabilizer; www.bioline.com), 3 mM MgCl2, 0.15 mM dNTP mix, 1 unit Taq polymerase, 0.1 μg DNA sample, and 0.075 μM each of the four following Alu-based PCR primers: AluY278T18 (5′-GAGCGAGACTCCGTCTCA-3′); AluY66H21 (5′-TGGTCTCGATCTCCTGACCTC-3′); R12A/267 (5′-AGCGAGACTCCG-3′) and L12A/8 (5′-TGAGCCACCGCG-3′). DNA denaturation was carried out at 95 °C for 5 min, followed by 30 cycles each of 30 s at 95 °C, 30 s at 50 °C, and 5 min at 72 °C, and finally another 7 min at 72 °C. Amplicons were purified with ethanol precipitation, sequenced on the Illumina-Solexa platform at Beijing Genomics Institute (Shenzhen, China) and mapped to the reference human genome hg19 .
AluScan sequencing data mapping and variant analysis
Paired-end sequencing reads were aligned to the GRCh37.p2 (Feb 2009) reference human genome using BWA (Burrows-Wheeler Aligner, version 0.6.1) with default settings . After BAM format-transfer and sorting by SAMtools (Sequence Alignment/Map, version 0.1.18) , the reads were further recalibrated and locally realigned using GATK (Genome Analysis Tool-Kit, version Lite-2.1-8-gbb7f038)  according to the standard framework .
The module ‘UnifiedGenotyper’ in GATK was employed to perform genotyping, and LOH and GOH callings for each sample were conducted with default settings. An LOH was defined as the conversion of a locus from heterozygosity in control to homozygosity in tumor, whereas a GOH was defined as the conversion of a locus from homozygosity in control to heterozygosity in tumor. Regions of read depths < 8 in either the tumor sample or its paired control would not be analyzed further. For homozygous reference loci, allele frequency must be 100 %. For homozygous non-reference loci, non-reference allele frequency must be 100 % with QD ≥ 20. For heterozygous loci, the non-reference allele must be ≥ 35 % and ≤ 65 % with Quality by Depth (QD) ≥ 4. Based on these conditions, with a ≥ 8 read depth, recognition of any site on the control sequence as a heterozygous site required a minimum of three reads bearing the non-reference allele. Since LOH would be called at this site on the tumor sequence only when the homozygous genotype was observed to be 100 %, an LOH event would be scored only if all these three non-reference allele reads in the control were no longer observed. Strand bias filter was employed to ensure SB values < −0.01 for both heterozygous loci and homozygous non-reference loci. See Additional file 2: Methods for more details on methods.
A variety of algorithms have been designed for CNV calling from whole genome sequencing (WGS) or exome sequencing data. However, the special features of AluScan data rendered difficult the calling of CNV using the calling algorithms designed for WGS or exome sequencing. Accordingly, the AluScanCNV method developed by us to call CNVs from AluScan and other types of sequence data, based on Geary-Hinkley transformation of read-depth ratios between either paired test-control samples or between test samples and a reference template constructed from reference samples , was employed in the present study for the identification of CN-gains and CN-losses.
Genic locations of mutations
The possible genic locations of the called LOH and GOH loci were identified through comparison with the Ensembl gene list from the UCSC database (http://genome.ucsc.edu/cgi-bin/hgTables), the TSGene database (Tumor Suppressor Gene: http://bioinfo.mc.vanderbilt.edu/TSGene/), and the NCG4.0 database (Network of Cancer Genes: http://ncg.kcl.ac.uk/).
Variant analysis of whole genome sequencing data
Raw whole genome sequencing data for a lung-to-liver metastatic cancer from a lung adenocarcinoma , and a hepatitis B positive hepatocellular carcinoma  were aligned to the GRCh37.p2 (Feb 2009) reference human genome by using BWA and GATK to extract all the overlapping sites between blood cell control and tumor tissue in the case of the lung-to-liver metastatic cancer, and between normal liver tissue control and tumor tissue in the case of the primary liver cancer. Genotyping and variant callings were performed for each of these two DNA samples with the ‘UnifiedGenotyper’ module in GATK with default settings as described for the AluScan samples. Regions of read depths <15 were first filtered out. Further filtration was achieved using the following criteria: for the homozygous reference loci, allele frequency must be 100 % and Quality by Depth (QD) ≥ 1; for heterozygous loci, non-reference allele frequency must be ≥ 25 % and ≤ 75 % with QD ≥ 4; for homozygous non-reference loci, non-reference allele frequency must be 100 % with QD ≥20. Phred-scaled p-value using Fisher’s exact test to detect strand bias (FS) was employed to ensure FS value ≤ 12 for both heterozygous loci and homozygous non-reference loci. The percentile rate of LOH (or GOH) was estimated by dividing the number of LOHs (or GOHs) by the number of total Mm sites (or MM or mm sites) that had passed through the filtration steps described above, and multiplying by 100 %.
High rate of LOH occurrence
The genome-wide single base-resolution LOHs detected in thirty primary and secondary cancers using the AluScan platform are summarized in Table 1 and Additional file 3: Table S1. In these tables, any allele at a residue in the genome that corresponded to the allele represented in the reference human genome hg19 was regarded as M, viz. the reference or major allele; on the other hand, any allele that differed from the allele on hg19 was regarded as m, viz. the minor allele. On this basis, the results in the tables showed that the heterozygous Mm residues in the genomes of gastric cancers, gliomas, leukemias, liver cancers and primary and secondary lung cancers all displayed exceptionally high mutation rates leading to an all-sample RMm equal to 9.15 % of all Mm residues analyzed, with a great majority of the mutations giving rise to LOH to yield an RLOH of 8.59 %. In contrast, RGOH-M for gain-of-heterozygosity was merely 1.54 × 10−3 % of all MM residues analyzed.
Accordingly, based on the all-sample averages, the RLOH/RGOH-M ratio for the thirty cancers was 8.59 % / 1.54 × 10−3 % = 5.58 × 103. Because the number of analyzed residues (viz. base pairs that were mapped in both the paired cancer and control samples) varied with the cancer-control pair, aggregate average RLOH was also estimated directly as the ratio [all LOHs detected] / [all Mm residues analyzed] in the thirty cancer-control pairs, and aggregate average RGOH-M as the ratio [all GOH mutations of MM residues detected] / [all MM residues analyzed]: thereby the aggregate RLOH/RGOH-M ratio for the 30 cancers was 6.91 % / 1.30 × 10−3 % = 5.32 × 103. Thus, by either route, RLOH was more than 5000-fold greater than RGOH-M. Given the far greater number of MM residues than Mm residues in genomic sequences, but far smaller RGOH-M than RLOH values, the total number of GOH-M mutations in the 30 cancers were of the same order of magnitude as the total number of LOH mutations (6360 and 8315 respectively) (Additional file 4: Table S3.31). The vastly unequal rates of LOH and GOH occurrences also rendered unlikely that the massive numbers of LOHs arose from technical errors, which would have produced LOHs and GOHs randomly at comparable rates. In fact, based on the use of high density whole genome SNP arrays, occurrence of copy number-neutral LOHs was also found to be frequent in gastrointestinal stromal tumors where contamination of tumor samples with normal cells was generally low, leading to the suggestion that the frequency of copy number-neutral LOHs might tend to be underestimated in solid tumors on account of the low percentages of tumor cells in the samples . In the present study, as indicated in Methods, all MM or mm genotypes arising from Mm genotypes must be 100 % for them to be called as LOHs in order to minimize or obviate the effects of varied percentages of tumor cells in a tumor sample.
Chromosomal distribution of LOH
The LOHs detected in the thirty cancers analyzed by AluScan sequencing (Table 1) were mostly distributed over interstitial sites along the lengths of different chromosomes without extraordinary clustering at the ends of chromosomes (Figs. 1 and 2). This fundamentally interstitial character of cancer LOHs was readily discerned in the leukemias (e.g. Fig. 1c) where the majority of the chromosomes displayed sparse, isolated LOH occurrences. It was equally evident in chromosomal regions with a high density of LOHs, e.g. chr 19, chr 17, chr 1p and chr 22q in Fig. 2a, where dense LOHs were closely interspersed with dense GOHs, indicating that the dense LOHs could not be derived from long stretches of CN-losses which would be incompatible with the co-occurrence of dense GOHs within the same stretches.
Detection of CNV revealed some CNV-dense segments, e.g. CN-losses in chr 1p, chr 9p and chr 9q of Glioma 1 which coincided with frequently observed CN-losses in gliomas [16–18], and CN-gains in chr 1q of both Glioma 1 and Lung-Brain 4 as well as chr 19p of Glioma 4 (Fig. 1b, Additional file 5: Table S6; Additional file 6: Figure S1), but most of the LOHs observed in the cancers apart from the gliomas were copy-neutral ones unassociated with either CN-losses or CN-gains (Additional file 7: Table S4). Even for the glioma samples, the aggregate average RGOH-M and RLOH were only changed from 2.69 × 10−3 to 2.45 × 10−3 % and from 12.29 to 12.79 % respectively after subtraction of CNV-associated GOHs and LOHs (Table 1). Locations of common and rare fragile sites  overlapped some LOH-dense regions such as those in chr 1p, chr 11q, chr 12q, chr 19p, chr19q, and chr 22q, but not all LOH-dense regions (Fig. 2a).
In the lung-to-liver metastatic cancer and primary liver cancer analyzed by WGS sequencing, the LOHs detected were similarly found to be mainly interstitial along the lengths of chromosomes rather than concentrated near the ends of chromosomes (Fig. 2b, c). The GOHs were also thickly interspersed with the LOHs along various chromosomes, again ruling out extended segments of CN-losses as a major cause of the LOHs (Additional file 8: Table S5). While the AluScan results (Fig. 2a) showed some of the chromosomes such as chr 19, chr 17, and chr 22q to be particularly enriched with LOHs and GOHs, the LOHs and GOHs revealed by WGS were more uniformly distributed among different autosomal chromosomes, which is in agreement with the elevated density of Alu elements in chr 19 and chr 17 and therefore enhanced sequence capture from these chromosomes by AluScan (Additional file 4: Table S3.34). Earlier we also found single nucleotide variations to be elevated in the vicinity of Alu elements .
Notably, the RLOH/RGOH-M ratio was ~640-3200 for the two sets of whole genome sequencing (WGS) data obtained by Ju et al.  and Ouyang et al.  (bottom two lines respectively, Table 1), compared to ~5300 for the all sample AluScan-based aggregate results. Therefore these two sequencing platforms were in agreement regarding the far greater percentile LOH mutations in cancer genomes compared to GOH mutations.
Preferences for reference alleles and transitional changes in cancer LOHs
In an LOH event, a heterozygous Mm residue with two different allelic bases on homologous chromosomes is mutated to a homozygous MM residue or a homozygous mm residue. In this regard, any residue in a human genome can be classified into the A-, G-, T- or C-family, depending on whether the reference allele at the same nucleotide position in the hg19 reference was an A, G, T or C. A preference for the reference (or M) allele represented by dark red or blue columns over the non-reference (or m) allele represented by light red or blue columns was observed for the LOHs arising from all 12 different types of Mm residues (Fig. 3a). Furthermore, transitional LOHs that produced the MM genotype (dark red columns), converting Ag to AA, Ga to GG, Tc to TT and Ct to CC, were far more prominent than transversional LOHs that produced the MM genotype (dark blue columns), converting At to AA, Ac to AA, Gt to GG, etc. Since the M-allele preference applied to LOHs arising from Mm residues belonging to all of the A-, G-, T- and C-families, it was distinct from GC-biased gene conversions .
M-allele preference was also displayed by the LOHs in the secondary lung-to-liver metastatic cancer  and the primary liver cancer  analyzed using WGS. Amongst the LOHs giving rise to MM genotypes in Fig. 3b, there was greater prominence of transitional LOHs than transversional LOHs as well. In Fig. 3c, the Ag-AA, Tc-TT and Ct-CC transitions remained prominent, but two of the transversions, viz. Ac-AA and Tg-TT, were likewise pronounced; whether this was related to the use of normal liver tissue as control in this instance remains to be determined. Thus the AluScan and WGS results were in accord with respect to M-allele preference, and both favored a prominence of transitional LOHs. A greater prominence of transitional relative to transversional LOHs is difficult to explain in terms of technical biases but consistent with the known higher transitional mutation rates relative to transversional mutation rates.
Occurrence of LOH and GOH in cancer related genes
A large number of tumor suppressor genes and other cancer related genes in the AluScan sequence reads of the 30 cancer genomes, identified based on the TSGene and NCG4.0 databases, underwent LOH and GOH mutations (Additional file 9: Table S7). Table 2 shows the genes that displayed such mutations in multiple cancer samples, which would represent genes with a high likelihood to participate in the oncogenic pathways of the various cancers.
To facilitate delineation of the genetic basis of human diseases, the bioinformatics tool CIPHER  has been developed to predict and prioritize disease genes based on the concordance between human protein network and disease phenotype network. Figure 4a shows an interaction network module associated with the solid tumor group, where the high-risk genes identified for this group by CIPHER based on database STRING  (Additional file 10: Table S8) interact with high likelihood genes from Table 2 that displayed multiple LOH and/or GOH mutations in the four solid tumor types (shown as color-coded sectors). Notably the genes CTNNA3, DLGAP1, CDH1, CHEK2, SAE1, SP1, RPS6KB1, AXIN1 and DNAJA3 were included in the high-likelihood genes in Table 2 as well as the CIPHER-identified high-risk genes. The network in Fig. 4a illustrates the utility of the identified multi-LOH/GOH genes in combination with CIPHER for analyzing potential protein-protein interactions in oncogenic networks.
Besides the multiple LOH- or GOH-bearing cancer-related genes detected, a host of other genes from the Ensembl database that were unlisted in either TSGene or NCG4.0 also underwent LOH and GOH mutations (Additional file 9: Tables S7E-F), and those unlisted genes that displayed highest occurrences of LOH or GOH, or in the largest numbers of cancer samples, are shown in Tables 2E-F, e.g. FRMD4A with fourteen LOHs, CAMK1D with 12 LOHs and 14 GOHs, AFAP1 with twelve GOHs, and LARGE with eleven GOHs. Based on the multiplicity of LOH and/or GOH mutations displayed by these genes in the various cancer samples, they would represent potential cancer-related genes that merit further investigation: eight of these genes, viz. DLGAP1, SHANK2, EIF2B5, SAE1, INSR, ATXN1 and RBFOX1 are included in Fig. 4a along with 20 tumor-suppressor and cancer related genes from Table 2A-D. These findings underline the usefulness of comprehensive LOH and GOH tracking in multiple cancers for uncovering potential cancer-related genes.
On account of the complexity of cancer cells, genomic studies provide an excellent approach to find surprises . In the present study, a characterization of the landscape of cancer LOHs revealed the surprisingly massive rates of LOH formation in various cancers, far exceeding the rates of GOH-M formation, and these cancer LOHs displayed a number of special properties.
Features of cancer copy-neutral LOHs
Unequal incidences in solid tumors and leukemias
In Table 1, the average leukemia RLOH of 0.61 % Mm residues was significantly lower than the average solid tumor RLOH of 10.2 % Mm residues with p = 0.0003 (Fig. 4b), in accord with previous reports of lower mutation rates in leukemias compared to solid tumors [25, 26]. Furthermore, insofar that the leukemia and solid tumor samples were analyzed using the same procedures, the leukemia RLOH of 0.61 % Mm residues suggests that the maximum technical error in the estimation of both leukemia and solid tumor RLOH incurred by false-positive calling of LOHs would not exceed 0.61 % Mm residues, amounting to only 0.61/10.2, or 6 % error, for the solid tumors. At the other extreme, the exact causes for the extraordinary >20 % RLOH rates of Gastric 1, Glioma 3, Glioma 5, Liver 1, and Lung 5 were undetermined; these cancers could be entering a terminal state of cellular disarray, and contribution from treatment modality-induced chromosomal instability also could not be ruled out.
Among different chromosomes, dense LOHs were present over large portions of chr 19, chr 17, chr 16p, chr 22q, and parts of chr 1p, chr 6p, chr 9q and chr 11q, but were relatively sparse in chr 4, chr 8, chr 13, chr 18 and chr 21 (Fig. 2a). Some but not all of the dense LOH regions overlapped with the locations of known common or rare fragile sites. Association of LOHs with CNVs was evident in some instances, notably in Glioma 1, but such CNV-associated LOHs represented only a minor fraction of the LOHs observed (Fig. 1, Additional file 7: Table S4).
The mainly interstitial character of the cancer LOHs was evident from the well-spaced LOH occurrences in the sparse LOH regions of Fig. 2a and most regions in Fig. 2b and c, as well as the distribution of LOH fragment lengths showing a substantial fraction of fragments that were ≤ 1 Mb in size (Fig. 5; Additional file 11: Table S9). It was also evident in the LOH-dense regions in Fig. 2a, where the crowded interspersion of LOHs and GOHs was incompatible with the dense LOHs being the result of any extended stretches of CN-losses. The interstitial nature of the major fraction of copy-neutral LOHs observed supports gene conversion being an important mechanism in the production of cancer LOHs.
Reduced LOHs on X-chromosomes
RLOH, the rate of LOH occurrence expressed as the percentile conversion of Mm residues analyzed, varied between different chromosomes (Fig. 6a). Among the autosomal chromosomes, this rate ranged from 5.5 % in chr 20 to 8.4 % in chr 13 (Additional file 4: Table S3.31). However, it was only 3.6 % in the X chromosomes in the female samples, which suggests that the inactive Barr-body configuration of one of the X-chromosomes  might constrain its participation in interhomolog recombination and LOH production relative to other chromosomes.
Interestingly, the four chromosomes exhibiting the highest levels of LOH/Mm, viz. chr 13, chr 22, chr 18 and chr 14, are known for their propensity to chromosomal instability as indicated by the frequency of trisomy [28, 29]. Also, the ~10−5 frequency of LOH occurrence on the various chromosomes expressed on a per bp basis (Fig. 6b) was similar to the 10−6 to 5 × 10−5 per bp frequencies displayed by TCGA ovarian cystadenoma and glioblastoma genomes .
Preferences for reference alleles and transitional changes
The partition of Mm genotypes between the MM and mm outcomes of LOH was characterized by a preference for MM over mm genotypes, with transitional LOHs being more prominent than transversional ones in the process (Additional file 12: Figure S2). Plausibly the M-alleles in the human genome were selected in the course of primate and human evolution to result in the adoption of MM genotypes at >99.9 % of the residues in the human genome. This bias in favor of the MM over Mm and mm genotypes suggests that, although the presence of m-alleles in the genome confers beneficial sequence diversity and gene-dosage modulation, functional advantages generally accrue to the M-alleles over m-alleles at most of the base positions in the genome. Accordingly, although key LOHs introducing mm-genotypes into tumor suppressor genes could facilitate oncogenesis , excessive accumulation of mm-genotypes might be detrimental to the cancer cells themselves, leading to selection against them and M-allele preference.
The taller column heights of transitional over transversional LOHs giving rise to MM-genotypes in Fig. 3a, b, and to a lesser extent in Fig. 3c, was in accord with the greater frequencies of transitional compared to transversional mutations in organisms.
Mode of double strand break repair in cancers
Since random point mutations would yield similar rates of mutation of Mm residues to yield LOHs and MM residues to yield GOHs, the vastly higher RLOH than RGOH-M values in Table 1 rule out point mutations as an important mechanism for the generation of cancer LOHs. This together with the limited role of CNV focuses attention on repair of double strand breaks (DSB) as a major source of cancer LOHs. In eukaryotic cells, DSBs are repaired by a spectrum of mechanisms through non-homologous end-joining, and homologous recombinations (HR) that include crossover pathways, break-induced replication, and synthesis-dependent strand-annealing employing a repair template supplied by sister chromatid in inter-sister chromatid recombination (ISR), a homologous chromosome in interhomolog recombination (IHR), or some ectopic sequence to bring about gene conversion .
Break-induced replication would be inconsistent with the largely interstitial character of the LOHs, and deletions of base residues at Mm sites due to non-homologous end-joining would not be called as SNVs (either LOHs or GOHs) in variant analysis (see Methods). In contrast, the HR process is known to produce relatively short patches of new DNA , in accord with the interstitial character of the cancer LOHs. Because ISR enhances both RMM and RLOH comparably whereas IHR enhances RLOH far in excess of RMM, the finding of RLOH > > RGOH-M with the different types of cancers strongly favored IHR over ISR as the underlying mechanism for the massive cancer LOHs. This conclusion was supported by the linear correlations between the levels of LOH, GOH-M and GOH-m, where GOH-M denotes GOHs arising from homozygous MM residues, and GOH-m denotes GOHs arising from homozygous mm residues (Fig. 7). That the occurrence of GOH mutations in MM and mm tagged along in correlation with the occurrence of LOH mutations was to be expected with IHR, which would cause not only the Mm residues involved to undergo LOH mutation due to the use of the allelic template for DSB repair, but also significant mutations of the Mm, MM and mm residues involved owing to the highly error-prone nature of the DNA polymerase employed for invading strand elongation in the course of IHR, increasing the mutation rate by up to 1400-fold . That RGOH-m was also much greater than RGOH-M would be consistent with the presence of extensive selection against mm-genotypes during cancer cell evolution, as suggested by the preference for M-alleles amongst the cancer LOHs.
The choice between ISR and IHR for DSB repair depends on a wide array of regulatory factors . IHR prevails in meiotic cells. In mitotic cells, ISR prevails over IHR on account of the close positioning of sister chromatids secured by cohesin [34–36], but IHR can be increased upon induction of a chromosomal DSB by two to three orders of magnitude up to 1 % of the target sites to yield short tract nonreciprocal gene conversions . Cellular systems are also known to undergo enhanced LOH production under special circumstances. In a Bloom mouse model, Blm-deficient ES cells exhibited elevated mitotic recombination rates with an 18-fold increase in somatic LOH . In aging diploid mother yeast cells, recombination rates could increase to 200-fold the rate in young cells to result in an age-induced switch to a hyper-recombination state . Thus the degenerative states in ageing yeast cells and the various types of tumors analyzed in Table 1 shared the common attribute of hyper-recombination. Interestingly, copy-neutral LOHs in ovarian cancer were found to be more frequent in older patients, suggesting that the effects of cancer and ageing could be additive in this regard .
The rate of DSB occurrence in cells has been estimated at about 50 DSBs per cell per cell cycle  or ten per day , and it is increased by both exogenous agents such as chemicals, ultraviolet and ionizing radiation, and endogenous events such as arrested replication forks, nucleases and reactive oxygen species from cellular metabolism [30, 43]. In the event that both leading and lagging strands of DNA are synthesized by discontinuous synthesis in human cells as has been suggested for Escherichia coli [44, 45], DNA synthesis itself can be a significant source of DSBs . Evidence for oncogene-induced DSBs has been provided by using the presence of p53 binding protein 1 (53BP1) nuclear foci as indicator of DSBs, whereby 10–20 foci per cell could be detected in cancer cell lines but not in proliferating normal cells .
Defective DNA-damage response
To cope with the continuous threat of DSBs, eukaryotic cells possess the capacity to mount a DNA-damage response (DDR) that arrests cell-cycle progression at the G1-S, intra-S and G2-M checkpoints to increase the time available for DNA repair; if the DNA damage cannot be removed, chronic DDR triggers cell death by apoptosis or cellular senescense [47–49]. Analysis of the relationship between DDR and oncogenesis has brought important insight into how oncogene activation-induced DNA hyper-replication could lead to S-phase DNA damage, onset of DDR and abrogation of cell cycle checkpoints, leading to a circumvention of the apoptosis and senescence pathways normally elicited by DDR, and hence oncogenesis [46, 50–52]. Notably, this chain of events, by diminishing or nullifying the action of the cell cycle checkpoints, not only would increase the influx of DSB-bearing DNA into the S-phase during the pre-oncogenesis phase to induce oncogenesis, but also may be expected to continue in the post-oncogenesis phase.
While usage of IHR for DSB repair is suppressed in favor of ISR in mitotic cells, the situation is radically altered upon oncogenesis and relaxation of checkpoints to enable the entry of DSB-bearing DNAs into S-phase. Under these circumstances, when a DSB-bearing chromatid replicates to yield a sister chromatid, the two sister chromatids will be unable to provide a useful repair template to one another. Instead, either a homologous chromosome or a homologous ectopic sequence will have to supply the requisite repair template . Given the reliable presence of the homologous chromosome as template compared to the haphazard availability of a homologous ectopic sequence, DSB repair by IHR will prevail giving rise to hyper-accumulation of LOH mutations in the cancer cells. Although DSB repair by HR in general can be a source of mutations arising from error-prone polymerases or replication forks [32, 54, 55], in the presence of intact cell cycle checkpoints ISR will be strongly preferred over IHR for DSB repair during S-G2, such that the error-prone polymerases would bring about comparable RGOH-M and RLOH rates instead of the vastly higher RLOH than RGOH-M rates found in cancer samples. Therefore only impaired cell cycle checkpoints arising from a defective DDR can lead to the hugely greater RLOH than RGOH-M displayed by different types of cancer cells.
The question has been posed regarding how might the requisite genetic changes, estimated to be about six mutations, be acquired to initiate oncogenesis, and whether mutator phenotypes participate in the process . This question is important not only to oncogenesis, but also to the post-oncogenesis phase with respect to the mutations needed to implement the manifold hallmarks of the neoplastic state, i.e. sustaining proliferative signaling, evading growth suppressors, resisting cell death, enabling replicative immortality, inducing angiogenesis, and activating invasion and metastasis , and complex metabolic reprogramming to support rapid growth even under conditions of fluctuating oxygen tension through enhanced glucose uptake, aerobic glycolysis, decreased conversion of pyruvate to acetyl-CoA etc. . Thus solid tumors compared to leukemias are confronted with irregular vascularization and more fluctuations in oxygen tension and nutrient supplies, which might contribute to the higher percentile LOH and GOH occurrences displayed by solid tumors compared to leukemias (Fig. 4b). Overall, the intimate relationship between defective DDR and cancer is clearly underlined by the numerous human genetic diseases that are associated with both DDR defects and predisposition to cancers including xeroderma pigmentosum, familial breast cancer, Bloom syndrome, MYH-associated polyposis, hereditary nonpolyposis colorectal cancer, etc. .
Accordingly, in the pre-oncogenesis phase, the weakening of cell cycle checkpoints caused by oncogene-induced DDR alterations can usher in mutations to initiate oncogenesis [46, 50–52]. In the post-oncogenesis phase, a defective DDR allows the continued accumulation of LOHs and GOHs, which can be self-amplifying insofar that some of the accumulated mutations can further alter DDR, to result in the wide landscape of mutations including LOHs, GOHs, CNVs, indels and chromosomal instabilities that characterize cancers, thus meeting fully the mutations needed for post-oncogenesis alterations which may be expected to outnumber those required for oncogenesis. Table 2 illustrates some of the LOH and GOH mutations in tumor suppressor and other cancer related genes potentially arising from interhomolog recombination enabled by a defective DDR that could be utilized by various cancers.
In conclusion, analysis of cancer LOHs and GOHs in the present study has revealed evidence for the occurrence of LOHs and tag-along GOHs in cancers brought about by repair of DSBs through interhomolog recombination under conditions of relaxed cell cycle checkpoints due to a defective DDR. On this basis, cancer may be regarded foremost as a disease of the DNA-damage response, where the mutator phenotype arising from DDR derangement provides a unified mechanism for generating interhomolog recombination-induced mutations (Fig. 8) to drive the initiation, development and aggressiveness of the neoplastic state from its oncogenic beginning to its terminal stages of unconstrained growth and proliferation.
Availability of supporting data
All supporting data are included as additional files.
Comparative genome hybridization
Copy number variation
DNA damage response
Double strand break
Fisher’s exact test to detect strand bias
Genome Analysis Tool-Kit
Gain of heterozygosity
Inter-sister chromatid recombination
Loss of heterozygosity
Quality by Depth
Percentile mutation of MM residues in the form of GOH
Percentile mutation of mm residues in the form of GOH
Percentile mutation of Mm residues in the form of LOH
Percentile mutation of Mm residues
Percentile mutation of mm residues
Single nucleotide polymorphism
Whole genome sequencing
Knudson AG. Antioncogenes and human cancer. Proc Natl Acad Sci USA. 1993;90:10914–21.
Zin R, Pham K, Ashleigh M, Ravine D, Waring P, Charles A. SNP-based arrays complement classic cytogenetics in the detection of chromosomal aberrations in Wilms’ tumor. Cancer Genet. 2012;205:80–93.
Abkevich V, Timms KM, Hennessy BT, Potter J, Carey MS, Meyer LA, et al. Patterns of genomic loss of heterozygosity predict homologous recombination repair defects in epithelial ovarian cancer. Br J Cancer. 2012;107:1776–82.
Lourenço N, Hélias-Rodzewicz Z, Bachet JB, Brahimi-Adouane S, Jardin F, van Nhieu JT, et al. Copy-neutral loss of heterozygosity and chromosome gains and losses are frequent in gastrointestinal stromal tumors. Mol Cancer. 2014;13:246.
Ha G, Roth A, Lai D, Bashashati A, Ding J, Goya R, et al. Integrative analysis of genome-wide loss of heterozygosity and monoallelic expression at nucleotide resolution reveals disrupted pathways in triple-negative breast cancer. Genome Res. 2012;22:1995–2007.
Pedersen BS, De S. Loss of heterozygous preferentially occurs in early replicating regions in cancer genomes. Nucl Acid Res. 2013;41:7615–24.
Mei L, Ding X, Tsang SY, Pun FW, Ng SK, Yang J, et al. AluScan: a method for genome-wide scanning of sequence and structure variations in the human genome. BMC Genomics. 2011;12:564–72.
Ju YS, Lee WC, Shin JY, Lee S, Bleazard T, Won JK, et al. A transforming KIF5B and RET gene fusion in lung adenocarcinoma revealed from whole-genome and transcriptome sequencing. Genome Res. 2012;22:436–45.
Ouyang L, Lee J, Park CK, Mao M, Shi Y, Gong Z, et al. Whole-genome sequencing of matched primary and metastatic hepatocellular carcinomas. BMC Med Genomics. 2014;7:2.
Genome Reference Consortium. UCSC Genome Browser, Index of /goldenPath/hg19/chromosomes.http://hgdownload.soe.ucsc.edu/goldenPath/hg19/chromosomes/. Accessed 10 Nov 2014.
Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25:1754–60.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map (SAM) Format and SAMtools. Bioinformatics. 2009;25:2078–9.
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:1297–303.
DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43:491–8.
Yang JF, Ding XF, Chen L, Mat WK, Xu MZ, Chen JF, et al. Copy number variation analysis based on AluScan sequences. J Clin Bioinformatics. 2014;4:15.
Reyes-Botero G, Dehais C, Idbaih A, Martin-Duverneuil N, Lahutte M, Carpentier C, et al. Contrast enhancement in 1p/19q-codeleted anaplastic oligodendrogliomas is associated with 9p loss, genomic instability, and angiogenic gene expression. Neuro Oncol. 2014;16:662–70.
Boots-Sprenger SH, Sijben A, Rijntjes J, Tops BB, Idema AJ, Rivera AL, et al. Significance of complete 1p/19q co-deletion, IDH1 mutation and MGMT promoter methylation in gliomas: use with caution. Mod Pathol. 2013;26:922–9.
Coco S, Valdora F, Bonassi S, Scaruffi P, Stigliani S, Oberthuer A, et al. Chromosome 9q and 16q loss identified by genome-wide pooled-analysis are associated with tumor aggressiveness in patients with classic medulloblastoma. Omics. 2011;15:273–80.
Schwartz M, Zlotorynski E, Kerem B. The molecular basis of common and rare fragile sites. Cancer Lett. 2006;232:13–26.
Ng SK, Xue H. Alu-associated enhancement of single nucleotide polymorphisms in the human genome. Gene. 2006;368:110–6.
Munch K, Mailund T, Dutheil JY, Schierup MH. A fine-scale recombination map of the human-chimpanzee ancestor reveals faster changes in humans than in chimpanzees and a strong impact of GC-baised gene conversion. Genome Res. 2014;24:467–74.
Wu X, Jiang R, Zhang MQ, Li S. Network-based global inference of human disease genes. Mol Systems Biol. 2008;4:189.
Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, et al. STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acid Res. 2013;41:D808–15.
Garraway LA, Lander ES. Lessons from the cancer genome. Cell. 2013;28:17–27.
Downing JR, Wilson RK, Zhang J, Mardis ER, Pui CH, Ding L, et al. The Pediatric Cancer Genome Project. Nat Genet. 2012;44:619–22.
Kandoth C, McLellan MD, Vandin F, Ye K, Niu B, Lu C, et al. Mutational landscape and significance across 12 major cancer types. Nature. 2013;502:333–41.
Walker CL, Cargile CB, Floy KM, Delannoy M, Migeon BR. The Barr body is a looped X chromosome formed by telomere association. Proc Natl Acad Sci USA. 1991;88:6191–5.
Morris JK, Savva GM. The risk of fetal loss following a prenatal diagnosis of trisomy 13 or trisomy 18. Am J Med Genet Part A. 2008;146:827–32.
Kanamori M, Yasuda T, Hori T, Suzuki K. Giant invasive sacral Schwannoma showing chromosomal numerical aberrations [−14, +18, +22]. Asian Spine J. 2013;7:227–31.
Mehta A, Haber JE. Sources of DNA double-strand breaks and models of recombinational repair. Cold Spring Harb Perspect Biol. 2014;6:a016428.
Jain S, Sugawara N, Lydeard J, Vaze M, Le Gac NT, Haber JE. A recombination execution checkpoint regulates the choice of homologous recombination pathway during DNA double-strand break repair. Genes Dev. 2009;23:291–303.
Hicks WM, Kim M, Haber JE. Increased mutagenesis and unique mutation signature associated with mitotic gene conversion. Science. 2010;329:82–5.
Chapman JR, Taylor MRG, Boulton SJ. Playing the end game: DNA double-strand break repair pathway choice. Mol Cell. 2012;47:497–510.
Kadyk LC, Hartwell LH. Sister chromatids are preferred over homologous substrates for recombination repair in Saccharomyces cerevisiae. Genetics. 1992;132:387–402.
Johnson RD, Jasin M. Double-strand-break-induced homologous recombination in mammalian cells. Biochem Soc Trans. 2001;29:196–201.
Sjogren C, Strom L. S-phase and DNA damage activated establishment of sister chromatid cohesion-importance for DNA repair. Exp Cell Res. 2010;316:1445–53.
Taghian DG, Nickoloff JA. Chromosomal double-strand breaks induce gene conversion at high frequency in mammalian cells. Mol Cell Biol. 1997;17:6386–93.
Luo G, Santoro IM, McDaniel LD, Nishijima I, Mills M, Youssoufian H, et al. Cancer predisposition caused by elevated mitotic recombination in Bloom mice. Nat Genet. 2000;26:424–9.
McMurray MA, Gottschling DE. An age-induced switch to a hyper- recombinational state. Science. 2003;301:1908–11.
Pedersen BS, Konstantinopoulos PA, Spillman MA, De S. Copy neutral loss of heterozygosity is more frequent in older ovarian cancer patients. Genes Chromosomes Cancer. 2013;52:794–801.
Vilenchik MM, Knudson AG. Endogenous DNA double-strand breaks: Production, fidelity of repair and induction of cancer. Proc Natl Acad Sci USA. 2003;100:12871–6.
Gospodinov A, Herceg Z. Chromatin structure in double strand break repair. DNA Repair. 2013;12:800–10.
Shrivastav M, De Haro LP, Nickoloff JA. Regulation of DNA double-strand break repair pathway choice. Cell Res. 2008;18:134–47.
Amado L, Kuzminov A. The replication intermediates in Escherichioa coli are not the product of DNA processing or uracil excision. J Biol Chem. 2006;281:22635–46.
Lehmann AR, Fuchs RP. Gaps and forks in DNA replication: rediscovering old models. DNA Repair. 2006;5:1495–8.
Halazonetis TD, Gorgoulis VG, Bartek J. An oncogene-induced DNA damage model for cancer development. Science. 2008;319:1352–5.
Jackson SP. The DNA-damage response: new molecular insights and new approaches to cancer therapy. Biochem Soc Trans. 2009;37:483–94.
Kastan MB, Bartek J. Cell-cycle checkpoints and cancer. Nature. 2004;432:316–23.
Jackson SP, Bartek J. The DNA-damage response in human biology and disease. Nature. 2009;461:1071–8.
Bartkova J, Hořejší Z, Koed K, Krämer A, Tort F, Zieger K, et al. DNA damage response as a candidate anti-cancer barrier in early human tumorigenesis. Nature. 2005;434:864–70.
Gorgoulis VG, Vassiliou LVF, Karakaidos P, Zacharatos P, Kotsinas A, Liloglou T, et al. Activation of the DNA damage checkpoint and genomic instability in human precancerous lesions. Nature. 2005;434:907–13.
Di Micco R, Fumagalli M, Cicalese A, Piccinin S, Gasparini P, Luise C, et al. Oncogene-induced senescence is a DNA damage response triggered by DNA hyper-replication. Nature. 2006;444:638–42.
Delacote F, Lopez BS. Importance of the cell cycle phase for the choice of the appropriate DSB repair pathway for genome stability maintenance. Cell Cycle. 2008;7:33–8.
Michel B, Flores MJ, Viguera E, Grompone G, Seigneur M, Bidnenko V. Rescue of arrested replication forks by homologous recombination. Proceedings of the National Academy of Sciences. 2001;98:8181–8.
Iraqui I, Chekkal Y, Jmari N, Pietrobon V, Fréon K, Costes A, et al. (2012). Recovery of arrested replication forks by homologous recombination is error-prone. PLoS genetics. 2012;8:e1002976.
Beckman RA, Loeb LA. Efficiency of carcinogenesis with and without a mutator mutation. Proc Natl Acad Sci USA. 2006;103:14140–5.
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144:646–74.
Kroemer G, Pouyssegur J. Tumor cell metabolism: cancer’s Achilles’ heel. Cancer Cell. 2008;13:472–82.
Ciccia A, Elledge SJ. The DNA damage response: making it safe to play with knives. Mol Cell. 2010;40:179–204.
Broad Institute TCGA Genome Data Analysis Center (2014): SNP6 Copy number analysis (GISTIC2). Broad Institute of MIT and Harvard. http://gdac.broadinstitute.org/runs/analyses__2014_10_17/reports/cancer/CESC/CopyNumber_Gistic2/nozzle.html. doi:10.7908/C1D21WFF; doi:10.7908/C12Z14C9; doi:10.7908/C1VH5MP9; doi:10.7908/C1XP73SJ; doi:10.7908/C1MP525H; doi:10.7908/C1T152H9; doi:10.7908/C11G0K5G; doi:10.7908/C1N29VT5; doi:10.7908/C1Z60MXC; doi:10.7908/C1X34WCP; doi:10.7908/C1S46QT9; doi:10.7908/C1K64H04; doi:10.7908/C17943KQ; doi:10.7908/C11G0K6X; doi:10.7908/C1QN65PF; doi:10.7908/C1TB15VM; doi:10.7908/C1CZ3633; doi:10.7908/C15B01DV; doi:10.7908/C17W6B4D.
The study was supported by grants to H. Xue from University Grants Council of Hong Kong SAR (VPRDO09/10.SC08, VPRDO14SC01, DG14SC02, SRFI11SC06 and SRFI11SC06PG), grants to S. Li from National Science Foundation of China (91229201 and 81225025), grant to L. Zhang from 863 Program, Ministry of Science and Technology, China (2012AA02A201), and grants to J. F. Chen (National 973 Basic Research Program of China, No. 2013CB911300; National Natural Science Foundation of China, No. 81272469; and Natural Science Foundation of Jiangsu Province special clinical project No. BL2012016). Y. Kumar was recipient of Hong Kong PhD Fellowship from Government of Hong Kong SAR. F.W. Pun was recipient of Research Fellowship from HKUST Jockey Club Institute of Advanced Study. We thank Prof. J. Tze-Fei Wong for valuable discussion.
The authors declare no competing financial or non-financial interests.
HX conceived and initiated the study; LC, ZX, LX, XXH, GT, JMW, YL, WSP, WQW, LWZ, JFC and HYW organized and collected the clinical samples and data; YK, JY, TH, WKM, FWP, PL, THYC, XD, SKN, SYT, PZ and SL analyzed the data; and HX, YK, JY and TH wrote the paper. All authors read and approved the final manuscript.
Yogesh Kumar, Jianfeng Yang and Taobo Hu contributed equally to this work.
Different types of point mutations in cancer samples grouped by genotypes.
Summary of tumor-control pairs analyzed.
Numbers of genotypes and SNV mutations on different chromosomes in various cancer samples.
Complete list of CNV sites.
Chromosomal distributions of LOHs, GOHs and CNVs in cancer samples. Symbols representing the different types of mutations are given in legend of Fig. 1.
Complete list of LOH sites and their CNV status.
Complete list of GOH sites and their CNV status.
LOH- and GOH-bearing genes present in TSGene, NCG4.0 and Ensembl databases.
Proteins identified by CIPHER as high-risk for the solid tumor group. The gene ranks indicated in % for the high-risk genes in the four different types of cancers were predicted by CIPHER. Covered length refers to the average length of gene sequence captured by the 30 AluScans analyzed in Table 2.
LOH fragment lengths.
About this article
Cite this article
Kumar, Y., Yang, J., Hu, T. et al. Massive interstitial copy-neutral loss-of-heterozygosity as evidence for cancer being a disease of the DNA-damage response. BMC Med Genomics 8, 42 (2015). https://doi.org/10.1186/s12920-015-0104-2
- Copy number variation
- Double strand break repair
- Gene conversion
- Inter-homologous recombination