- Primary research
- Open Access
Oxidative stress-induced chromosome breaks within the ABL gene: a model for chromosome rearrangement in nasopharyngeal carcinoma
Human Genomics volume 12, Article number: 29 (2018)
The mechanism underlying chromosome rearrangement in nasopharyngeal carcinoma (NPC) remains elusive. It is known that most of the aetiological factors of NPC trigger oxidative stress. Oxidative stress is a potent apoptotic inducer. During apoptosis, chromatin cleavage and DNA fragmentation occur. However, cells may undergo DNA repair and survive apoptosis. Non-homologous end joining (NHEJ) pathway has been known as the primary DNA repair system in human cells. The NHEJ process may repair DNA ends without any homology, although region of microhomology (a few nucleotides) is usually utilised by this DNA repair system. Cells that evade apoptosis via erroneous DNA repair may carry chromosomal aberration. Apoptotic nuclease was found to be associated with nuclear matrix during apoptosis. Matrix association region/scaffold attachment region (MAR/SAR) is the binding site of the chromosomal DNA loop structure to the nuclear matrix. When apoptotic nuclease is associated with nuclear matrix during apoptosis, it potentially cleaves at MAR/SAR. Cells that survive apoptosis via compromised DNA repair may carry chromosome rearrangement contributing to NPC tumourigenesis. The Abelson murine leukaemia (ABL) gene at 9q34 was targeted in this study as 9q34 is a common region of loss in NPC. This study aimed to identify the chromosome breakages and/or rearrangements in the ABL gene in cells undergoing oxidative stress-induced apoptosis.
In the present study, in silico prediction of MAR/SAR was performed in the ABL gene. More than 80% of the predicted MAR/SAR sites are closely associated with previously reported patient breakpoint cluster regions (BCR). By using inverse polymerase chain reaction (IPCR), we demonstrated that hydrogen peroxide (H2O2)-induced apoptosis in normal nasopharyngeal epithelial and NPC cells led to chromosomal breakages within the ABL BCR that contains a MAR/SAR. Intriguingly, we detected two translocations in H2O2-treated cells. Region of microhomology was found at the translocation junctions. This observation is consistent with the operation of microhomology-mediated NHEJ.
Our findings suggested that oxidative stress-induced apoptosis may participate in chromosome rearrangements of NPC. A revised model for oxidative stress-induced apoptosis mediating chromosome rearrangement in NPC is proposed.
Nasopharyngeal carcinoma (NPC) is a malignant neoplasm derived from mucosal epithelium of the nasopharynx. According to the World Health Organization (WHO), NPC can be classified into three subtypes according to the degree of epithelial differentiation, namely keratinising squamous cell carcinoma (Type I), non-keratinising squamous cell carcinoma (Type II) and undifferentiated or poorly differentiated carcinoma (Type III) .
NPC is a rare malignancy in most parts of the world; the incidence rates are below one per 100,000 persons per year [2, 3]. However, there are a few well-known notable exceptions . The intermediate rates were reported in South-Eastern Asia, Northern Africa, the Middle East and Arctic Region [3, 4]. The highest incidence rate was observed among Southern Chinese living in central Guangdong province. The annual incidence rates for males and females in central Guangdong province are 23.3 per 100,000 and 8.9 per 100,000, respectively . The NPC incidence rates are generally increasing from Northern China to Southern China [3, 4]. In addition, an exceptionally high incidence rate has been reported among the Bidayuh people, the second biggest ethnic group in Sarawak, Malaysia. The age-adjusted rate of Sarawak residents is 13.5 per 100,000 and 6.2 per 100,000 in males and females, respectively. Although the average rate in Sarawak is intermediate, the incidence rate for Bidayuh people is about 50% higher than that in Hong Kong (a part of the Cantonese region of Guangdong province) .
NPC is strongly associated with Epstein-Barr Virus (EBV) infection [6,7,8] as well as dietary [9,10,11], environmental  and genetic factors [13, 14]. Several genetic aberrations have been reported to be related to the development of NPC, suggesting that NPC tumourigenesis involves multiple genetic changes. These include chromosomal gains or losses [15,16,17,18,19], loss of heterozygosity (LOH) [20,21,22,23], homozygous deletions [24,25,26,27], promoter hypermethylation of tumour suppressor genes [28,29,30,31] and shortening of chromosome telomeres [32, 33].
Although the consistent chromosome rearrangements have long been identified in NPC, the molecular mechanism underlying the chromosome rearrangements of NPC remains poorly understood. In addition to EBV infection, long-term exposures to nitrosamines, formaldehyde, cigarette smoke and wood dust have all been found to be associated with NPC [12, 34,35,36]. More recently, much concern has been raised about the association between chronic inflammation of sinonasal tract and NPC [37, 38]. It is remarkable that all these aetiological factors may trigger oxidative stress [39,40,41,42,43]. Oxidative stress is an imbalance of pro-oxidants and antioxidants resulting in a disruption of redox signalling and control. Pro-oxidants induce oxidative stress either through excessive production of reactive oxygen species (ROS) or inhibition of antioxidant systems . ROS are chemically-reactive molecules containing oxygen which include peroxyl RO·, hydroxyl radical OH·, superoxide O2·- and hydrogen peroxide H2O2 . ROS cause several kinds of DNA damages, including strand cleavage, base modification and DNA-protein cross-linkage [45, 46]. Importantly, formaldehyde and acrolein, a component of cigarette smoke, are reactive aldehydes. In addition, reactive aldehydes may also be produced endogenously during oxidative stress. Aldehydes may cause adduct formation that impairs the function of DNA, RNA and proteins via electrophile-nucleophile interaction. Exposure to environmental aldehydes has been shown to be associated with the onset and development of human diseases that involve oxidative stress. It has been suggested that environmental and endogenous aldehydes may interact additively and exacerbate the cellular oxidative damage .
An evaluation of the levels of 8-hydroxy-2′-deoxyguanosine (8-OHdG), a biomarker of oxidative DNA damage, had been done among NPC patients. The tissue and serum levels of 8-OHdG in NPC patients have been found to be significantly higher than those in control patients . Oxidative stress was suggested to play an important role in carcinogenesis . Since there is a strong link between the aetiological factors of NPC and oxidative stress, it is intriguing to investigate the role of oxidative stress in the molecular mechanisms underlying chromosome rearrangements of NPC.
Oxidative stress may induce apoptosis [50, 51]. H2O2 has been well known as an apoptotic inducer for various human cell types, including osteoblasts [52, 53], sarcoma cells , osteosarcoma cells , hepatoma cells , astrocytoma cells , Jurkat T lymphocytes  and Fanconi’s anaemia cells . Apoptosis or programmed cell death was first described by Kerr et al. (1972). Apoptosis is a type of genetically controlled cell suicide which occurs naturally in multicellular organisms in order to eliminate poisonous cells. Apoptosis is morphologically characterised by condensation of chromatin, fragmentation of nuclei, compaction of cytoplasmic organelles, cell shrinkage and cytoplasmic membrane blebbing [60,61,62]. Apoptosis is related to several biochemical events, including externalisation of phosphatidylserine (PS) on cell membrane, alteration in mitochondrial membrane potential (MMP), release of cytochrome c (cyt c) from mitochondria, caspase activation and internucleosomal cleavage of DNA .
The alteration of nuclear chromatin during apoptosis is often associated with fragmentation of the genomic DNA into high-molecular-weight (HMW) DNA of 30 to 50 and 200 to 300 kbp [64, 65]. These sizes of fragments have been suggested to be derived from the release of loops (50 kbp) or rosettes (300 kbp) of chromatin, probably when they become detached from their binding sites on the nuclear scaffold . Further degradation of the HMW DNA produces the internucleosomal DNA fragments of 180 to 200 bp [67, 68].
Cells undergoing apoptosis may recover from the execution phase of apoptosis upon DNA repair [69, 70]. There are two major double-strand breaks (DSBs) repair pathways, namely homologous recombination (HR) and non-homologous end joining (NHEJ) [71, 72]. Chromosomal DSB repair by HR is predominant during late S/G2 phases of the cell cycle. NHEJ is the more frequently used pathway that can repair a DSB at any time during the cell cycle. These two pathways have different degree of requirement for DNA homology. The HR pathway requires sufficient homology, usually more than 100 bp. Given that the HR DNA repair system ligates two DNA ends with homologous sequences, it gives rise to precise DNA repair. The NHEJ DNA repair system joins two DNA ends without intensive requirement of sequence homology. This pathway joins two DNA terminals with microhomology of a few base pairs [71, 73, 74]. NHEJ pathway has been shown to be prone to cause erroneous repair of DSBs. This may in turn lead to chromosomal aberrations . It has been suggested that interaction of the NHEJ DNA repair system with apoptosis can act as a mechanism leading to translocation in leukaemia .
Chromosomal breakage takes place in the initial stage of chromosome rearrangement and apoptotic DNA fragmentation. It has been observed that chromosome breaks do not randomly occur throughout a gene. Rather, chromosome breaks normally fall within certain regions that contain specific chromatin structures, such as matrix association region/scaffold attachment region (MAR/SAR) [76, 77]. MAR/SAR are DNA sequences where DNA loop structure attaches to nuclear scaffold/matrix proteins . There are two breakpoint cluster regions (BCR) identified in the AF9 gene. These two BCRs are bordered by two experimentally isolated MAR/SARs . The BCR of the mixed lineage leukaemia (MLL) gene has also been found to contain two MAR/SAR sequences . In addition, the most crucial apoptotic nuclease CAD has been reported to associate with the nuclear matrix of apoptotic cells .
We previously demonstrated that in normal nasopharyngeal epithelial and NPC cells, oxidative stress-induced apoptosis resulted in chromosome breaks in the AF9 gene located on chromosome 9p22. We further demonstrated that caspase-activated DNase (CAD) may be a major player in mediating the oxidative stress-induced chromosomal cleavages. A few chromosome breaks were identified within the AF9 region that was previously reported to participate in translocation in an acute lymphoblastic leukaemia (ALL) patient. These findings suggested that oxidative stress-induced apoptosis may play an important role in mediating chromosome rearrangements in NPC . In the present study, we further investigated the potential role of oxidative stress-induced apoptosis by targeting the Abelson murine leukaemia viral oncogene homologue 1 (ABL) gene located on chromosome 9q34. This study targeted the ABL gene because 9q34 is a common region of loss in NPC .
The ABL gene is a proto-oncogene which encodes a 150 kDa nonreceptor protein tyrosine kinase. It was first recognised as the cellular homologue of the v-abl oncogene product of the Abelson murine leukaemia virus [81, 82]. The ABL protein has a complex structure that contains many domains. These domains are found in proteins which are involved in the formation of complexes in signal transduction pathway. It has been demonstrated that overexpression of ABL in fibroblast resulted in growth arrest . The product of ABL-BCR fusion appears to be an abnormal kinase that stimulates the proliferation of myeloid cells leading to chronic myelogenous leukaemia (CML) . The ABL gene is 173,795 bp in length and it consists of 11 exons [Ensembl:ENSG00000097007]. The description of exons and introns in the ABL gene is shown in Additional file 1.
By using MAR/SAR recognition signature (MRS), we predicted 12 possible MAR/SAR sites in the ABL gene. We demonstrated that oxidative stress-induced apoptosis resulted in chromosome breaks in the ABL BCR which contains a MAR/SAR site. We detected shift translocations in H2O2-treated normal nasopharyngeal epithelial cells. Interestingly, we found region of microhomology at the breakpoint junctions. This observation suggests a role for NHEJ DNA repair system in mediating the translocation. At last, we illustrated the possible role of oxidative stress-induced apoptosis in mediating chromosome rearrangements in NPC via NHEJ DNA repair system.
In silico prediction of MAR/SAR by using MAR/SAR recognition signature (MRS)
Potential MAR/SAR sites in the ABL gene were predicted by using MRS. MRS is a bipartite sequence that is strongly associated with MAR/SAR . This bipartite sequence consists of 16 bp nucleotide motif (AWWRTAANNWWGNNNC) within a distance of 200 bp of the 8 bp nucleotide motif (AATAAYAA). However, for our preliminary results in the ABL gene, we only found one MRS (MAR/SAR 9 in Table 1) in the biochemically identified SAR1 . The distance between the 8 bp sequence element and the 16 bp sequence element was 248 bp. Therefore, in this study, we set the maximal distance between the two sequence elements at 250 bp.
By using MRS, we predicted 12 potential MAR/SAR sites in the ABL gene. The nucleotide positions of the MRSs with their sequence composition, relative orientation, distance between the two sequence elements and location of the MRSs in the exon or intron of the ABL gene are shown in Table 1. Out of the 12 predicted MAR/SAR sites, 9 were identified in intron 1 which is the largest intron (approximately 140 kb in length) in the ABL gene (approximately 175 kb in length) (MAR/SAR 1–9 in Table 1). One potential MAR/SAR site was separately found in intron 3 (MAR/SAR 10 in Table 1) and intron 10 (MAR/SAR 11 in Table 1). The distribution of the predicted MAR/SAR sites in the ABL gene is shown in Fig. 1. One MAR/SAR site (MAR/SAR 9 in Table 1) was predicted within the biochemically defined SAR1 which is located in the second intron 1 .
NP69 cells were either left untreated or treated with 100 μM of H2O2 for 16 and 24 h while HK1 cells were either left untreated or treated with 50 μM for 4 and 8 h. Cells treated with CPT was included as a positive control. The cells were then subjected to flow cytometric analyses of PS externalisation and MMP loss.
Phosphatidylserine (PS) externalisation
As shown in Fig. 2a i, the percentages of apoptotic cells detected in NP69 treated with 100 μM of H2O2 for 16 and 24 h were 2.82-fold (p = 0.000170) and 2.87-fold (p = 3.4346E−8) higher than that detected in the untreated control, respectively. The percentages of apoptotic cells detected in HK1 treated with 50 μM of H2O2 for 4 and 8 h were 1.48-fold (p = 0.005735) and 1.92-fold (p = 0.000477) higher than that detected in the untreated control, respectively (Fig. 2b i). Figure 2a i and b ii are the representative dot plot diagrams showing the apoptotic population of H2O2-treated NP69 and HK1 cells, respectively.
Mitochondrial membrane potential (MMP) loss
As shown in Fig. 3a i, the percentages of apoptotic cells detected in NP69 treated with 100 μM of H2O2 for 16 and 24 h were 2.45-fold (p = 0.006) and 2.25-fold (p = 0.002) higher than that detected in the untreated control, respectively. The percentages of apoptotic cells detected in HK1 treated with 50 μM of H2O2 for 4 and 8 h were 1.68-fold (p = 0.009) and 2.18-fold (p = 0.007) higher than that detected in the untreated control, respectively (Fig. 3b i). Figure 3a ii and b ii are the representative contour plot diagrams showing the apoptotic population of H2O2-treated NP69 and HK1 cells, respectively.
IPCR detection of chromosome breaks within the ABL gene mediated by stress-induced apoptosis
NP69 cells at confluency of 30–40% were treated with 10, 50 and 100 μM of H2O2 for 16 and 24 h while HK1 cells at optimal density were treated with 1, 10 or 50 μM of H2O2 for 2, 4, 6 and 8 h. For each cell line, an untreated sample was included to serve as a cell control. Nested IPCR was employed to identify chromosome breaks mediated by stress-induced apoptosis. The IPCR bands representing the ABL cleaved fragments detected were isolated, purified and sequenced.
Figures 4 and 5 show the IPCR results for H2O2-treated NP69 and HK1 cells, respectively. In the manipulation for nested IPCR, Age I (RE2 in Fig. 12) was used to linearise the cyclised DNA. If there is no breakage within the ABL gene, the IPCR product will be approximately 3 kb. On the contrary, if there is any breakage within the ABL gene, it should produce IPCR products which are smaller than 3 kb. As shown in Figures 4a and 5a, the ABL intact fragment of 3 kb is present in all of the samples. This amplification could serve as an internal control by proposing an optimal IPCR condition for the ABL gene. Besides, numerous IPCR bands of less than 3 kb were also obtained. However, these bands are less intense as compared with the intact fragment. This could be due to the competition between the intact fragments and the cleaved fragments for the amplification process. The intact fragments are usually more abundant as compared with the cleaved fragments. Consequently, the amplification of the cleaved fragments would be less efficient in the presence of the intact fragments. Therefore, double digestion with Age I and BsaA I or Age I and EcoR I (RE3 in Fig. 12) was used to eliminate competition from the intact fragments for the nested IPCR reaction. These two different digestions gave rise to the detection of chromosome breaks within different regions. With double digestion of Age I and EcoR I, numerous distinct IPCR bands of less than 3 kb which represent the cleaved ABL fragment were detected in H2O2-treated NP69 (Fig. 4b, lanes 4–9) and H2O2-treated HK1 (Fig. 5b, lanes 4, 7, 9, 10, 11, 12 and 13) cells. A few cleavage bands were detected in the untreated NP69 cells (Fig. 4b, lane 1) which might be due to endogenous DNA breaks in the minority of untreated cells which were unhealthy. However, in general, there were more cleaved fragments detected in the treated samples compared with the untreated sample.
As shown in Fig. 4c, the intact fragment of 3 kb was still detected upon double digestion with Age I and BsaA I. Most probably, this was due to incomplete digestion by Age I and BsaA I. Regardless of the incomplete digestion, several cleavage bands of different sizes were detected in NP69 cells treated with various concentrations of H2O2 for different time points (Fig. 4c, lanes 4–9).
DNA breakages were detected in cell samples treated with various concentrations of H2O2 at various time points. Based on the microscopic analysis and flow cytometric analyses, the optimal concentration and time point were determined. These optimal concentration and time point were used to repeat the experiments in NP69 and HK1 cells. For NP69 cells, we selected concentration of 100 μM with exposure time of 16 and 24 h. The microscopic analysis on H2O2-treated NP69 cells showed that cytoplasmic shrinkage was only observed in cells treated with 100 μM for 16 and 24 h (Additional file 2). In addition, we performed flow cytometric analyses of PS externalisation and MMP loss on NP69 cells treated with 100 μM for 16 and 24 h. In these two flow cytometric assays, both samples showed significantly higher percentage of apoptosis as compared with that of the untreated control (Figs. 2a i and 3a i).
For HK1 cells, we chose concentration of 50 μM with exposure time of 8 h. Cell blebbing, which is one of the major morphological changes, was relatively more prominent in HK1 cells treated with 50 μM of H2O2 for 8 h compared with those treated with lower concentration (1 and 10 μM) and shorter exposure time (2, 4 and 6 h) (Additional file 3). Moreover, the flow cytometric analysis of PS externalisation on sample treated with 50 μM of H2O2 for 8 h showed a more promising result compared with the sample treated with 50 μM of H2O2 for 4 h (Fig. 2b i).
H2O2 treatment was thus repeated in NP69 and HK1 cells with the selected concentration and time point. NP69 cells at confluency of 30–40% were either untreated or treated with 100 μM of H2O2 for 16 h while HK1 cells at confluency of 60–70% were either untreated or treated with 50 μM of H2O2 for 8 h. The cells were then harvested for gDNA extraction and nested IPCR. In the manipulation for nested IPCR, all the samples were subjected to double digestion with Age I and EcoR I (RE3 in Fig. 12).
Fig. 6a shows that numerous IPCR bands of less than 3 kb which represent the cleaved ABL gene detected in NP69 cells treated with H2O2 for 16 h (lanes 9, 10, 12 and 13) and 24 h (lanes 14–19). Three cleavage bands were identified in the untreated control (lanes 2 and 5). This might be due to spontaneous cell death of untreated cells as detected in our flow cytometric analyses. As shown in the bar chart in Fig. 6b, the cleavage frequencies of the ABL gene detected in NP69 cells treated with H2O2 for 16 and 24 h are 1.4-fold (p = 0.004966) and 1.8-fold (p = 0.000009) higher than that of the untreated control, respectively.
Similar findings were obtained from the H2O2 treatment of HK1 cells. Figure 7a is a representative gel picture showing the IPCR result obtained from this experiment. The untreated HK1 cells show a few cleavage bands (lanes 2–7) which most likely due to spontaneous cell death. In contrast, there were numerous cleavage bands identified in the H2O2-treated sample (lanes 8–13). The chart in Fig. 7b shows that the cleavage frequency of the ABL gene detected in H2O2-treated HK1 cells is 1.7-fold higher than that of untreated HK1 cells (p = 0.000197).
In order to confirm that these fragments were derived from the cleaved ABL gene, some of the cleavage IPCR bands were extracted and sequenced. The sequencing results show that they were all derived from the cleaved ABL gene. Table 2 shows the breakpoints identified within the ABL gene in the H2O2-treated cells. A map illustrating the positions of chromosome breaks in HK1 and NP69 cells relative to the MAR/SAR sequences within the ABL gene is shown in Fig. 8.
Intriguingly, we detected two shift translocations in H2O2-treated NP69 cells. Translocation is rarely observed in NPC, compared with deletion and addition. The first shift translocation was identified in NP69 cells exposed to 100 μM of H2O2 for 16 h. As shown in Fig. 9a, the translocated segment was derived from the human lipoma HMGIC fusion partner-like 3 (LHFPL3) gene which is located on chromosome 7. The LHFPL3 gene consisting of three exons is 578,576 bp in length. The description of exons and introns in the LHFPL3 gene is shown in Additional file 4. The translocated segment (228 bp) of the LHFPL3 gene is corresponding to coordinates 108006–108234 [Ensembl:ENSG00000187416]. The breakpoints (108,006 and 108,234) of the LHFPL3 gene were mapped within its first intron. Moreover, region of microhomology (TGCC) was found at the breakpoint junctions. The second shift translocation was identified in NP69 cells exposed to 10 μM of H2O2 for 24 h. The segment translocated to the ABL gene is derived from chromosome 5. The disabled homologue 2 (DAB) gene is 1,263,556 bp at the 5′ end of this segment while a gene encoding for hypothetical protein is 22,122 bp at the 3′ end (Fig. 9b).
Oxidative stress increases genomic instability  which in turn contributes to carcinogenesis [87, 88]. Excessive production of ROS has been associated with mutation and alteration of gene expression . Most of the aetiological factors of NPC were known to generate ROS. These aetiological factors include exposures to nitrosamines, cigarette smoke, formaldehyde and wood dust. EBV infection as well as chronic inflammation of sinonasal tract [39, 41,42,43].
In addition, formaldehyde and acrolein, a component of cigarette smoke, are reactive aldehydes which may impair the function of DNA, RNA and proteins through adduct formation. It has been suggested that the combined interactions of environmental aldehydes and endogenous aldehydes, which are produced during oxidative stress, may exacerbate the cellular oxidative damage .
Although the consistent chromosomal aberrations, such as deletion and addition, have long been identified in NPC, the underlying molecular mechanism requires further investigation. Apoptosis was suggested to participate in the chromosomal translocation process of leukaemia . Given that there is a strong association between the aetiological factors of NPC and oxidative stress, we intended to investigate the role of oxidative stress-induced apoptosis in mediating the chromosome rearrangements of NPC.
We demonstrated that hydrogen peroxide (H2O2), a strong oxidising agent, was able to induce apoptosis in normal nasopharyngeal epithelial cells (NP69) and NPC cells (HK1). Both exposure of PS and disruption of MMP are key events of apoptosis [89, 90]. By using flow cytometric analyses of PS externalisation and MMP loss, we detected significantly higher percentages of apoptosis in H2O2-treated NP69 and HK1 cells as compared with the untreated controls. As compared with NPC cells (HK1), longer exposure time and higher dosage of H2O2 were needed to trigger apoptosis in normal nasopharyngeal epithelial cells (NP69). There may be several possibilities that lead to this variation. The intrachromosomal instability in cancer cells is usually higher than that in normal cells [91, 92]. Furthermore, cancer cells may have a defective DNA repair system which is unable to restore the genomic integrity . These factors imply that cancer cells are more susceptible to DNA damage. When apoptosis is triggered by oxidative stress, DNA fragmentation occurs. Cells try to survive apoptosis through DNA repair. Therefore, as compared with normal cells, cancer cells which have a higher intrachromosomal instability or a defective DNA repair system are usually more vulnerable to apoptosis.
In our previous report, we identified chromosomal breakages within the AF9 gene in H2O2-treated NP69 and HK1 cells. In addition, inhibiting caspase-3 by caspase-3 inhibitor has abolished the AF9 gene cleavages mediated by H2O2-induced apoptosis. Given that caspase-3 is the main activator of CAD-mediated DNA fragmentation in apoptosis, our findings suggested that CAD might be the major player which mediated the chromosomal breakages in H2O2-induced apoptosis . It has been observed that CAD binds to the nuclear matrix during apoptosis . Due to the fact that MAR/SAR sequences are the sites where DNA interacts with the nuclear matrix , it is likely that CAD cleaves the DNA at MAR/SAR sequences when it associates with nuclear matrix. Intriguingly, our previous report demonstrated that oxidative stress-induced apoptosis caused chromosomal breakages within the AF9 BCR which is bordered by two MAR/SARs .
The present study targeted the ABL gene which is located on chromosome 9q34. This gene was targeted because 9q33-34 is one of the common deletion regions in NPC . The ABL gene is the most common fusion partner gene with the breakpoint cluster region (BCR) gene which is located on chromosome 22q11 . The reciprocal translocation t(9;22)(q34;q11) in CML was the first consistent chromosome rearrangement found in malignancy. The ABL-BCR fusion gene was named as the Philadelphia chromosome . This reciprocal translocation was found in approximately 92% of CML patients. Thus, the ABL-BCR fusion gene is recognised as the cytogenetic hallmark of patients suffering from this disease [94, 96]. The presence of Philadelphia chromosome was also reported in 20 to 55% of adults and 2 to 10% of children with acute lymphoblastic leukaemia (ALL)  and rarely (1 to 2%) in acute non-lymphoblastic leukaemia (ANLL) . There are three BCRs found within the ABL gene. The first BCR (BCRA) and the second BCR (BCRB) are located in intron 1b, whereas the third BCR (BCRC) spans through parts of introns 1b to 3. BCRC is the largest BCR of the ABL gene [77, 99].
One biochemically defined MAR/SAR has been previously identified within the BCRC of the ABL gene. This MAR/SAR was designated as SAR1. SAR1 was found within intron 1a . In the present study, we predicted MAR/SAR sites within the ABL gene by using MRS which was proposed to be strongly associated with MAR/SAR . It has been found that the two sequence elements of the MRS exist at a position near the dyad axis of the nucleosome. The wrapping of the DNA around the histone protein complex causes the two sequence elements of the MRS to be physically close together even if they are non-adjacent on the linear DNA. The close proximity between the two sequence elements of the MRS on the positioned nucleosome allows them to generate a protein-binding site in MAR/SAR .
The variation in the distance between the two sequence elements suggested a relation of the MRS to nucleosome organisation. In the Drosophila histone cluster, there was a MAR/SAR identified between the histone H1 and H3 genes. This MAR/SAR was found to contain a few nucleosomes and two MRSs. It was observed that the position of the two MRSs on their respective nucleosomes is similar. The first MRS, where the two sequence elements are overlapping, is found on the dyad axis of a nucleosome. The second MRS, where the two sequence elements are 145 bp apart, is located near the entry and exit sites of a nucleosome. Although the two sequence elements of the MRS are spatially distant, they are brought close together when the DNA is turned around the histone core .
A nucleosome comprises a nucleosome core and a ‘linker’ DNA. The nucleosome core contains 145–147 bp of DNA wrapped around a core histone octamer. The histone octamer consists of two molecules each of the four core histones, namely, H2A, H2B, H3 and H4 . It has been known that the length of ‘linker’ DNA ranges from 15 to 100 bp, depending on the cell types. The ‘linker’ DNA connects one nucleosome to the other (reviewed in ). The nucleosome repeat length (NRL) refers to the length of nucleosomal DNA (145–147 bp) plus the length of linker DNA (15–100 bp) [94, 95]. Using the micrococcal nuclease assay, the NRL has been reported to range from 160 to 240 bp [96, 97].
In the studies by van Drunen et al. (1999), the distance between the two sequence elements of the MRS has been suggested to be within 200 bp . SAR prediction/SAR prediction presently performed in the ABL gene, there was only one MAR/SAR site (MAR/SAR 9 in Table 1) predicted in the experimentally isolated SAR1. The distance between the 8 bp sequence element and the 16 bp sequence element was found to be 248 bp. Given that the NRL may exceed 200 bp, for the mapping of MRS in the present study, the maximal distance between the 8 bp sequence element and the 16 bp sequence element was set at 250 bp. Besides, it is also possible that the two sequence elements, which are 248 bp apart, are located separately on two adjacent nucleosomes. In the positioned nucleosomes, interaction between two adjacent nucleosomes may happen. Thus, it seems possible that even if the two sequence elements are individually located on two adjacent nucleosomes, the wrapping of DNA around the histone protein complex may still cause them to be physically close together and enable them to generate a protein binding site.
We predicted 12 potential MAR/SAR sites within the ABL gene. One MAR/SAR site was predicted within the biochemically defined SAR1. Interestingly, 10 out of these 12 (> 80%) potential MAR/SAR sites are closely associated with the BCRs of the ABL gene (Fig. 1). MAR/SARs 1 and 2 were predicted next to BCRA. MAR/SAR 3 was found within BCRB. MAR/SARs 4 to 10 were identified within BCRC.
By using IPCR, we identified chromosome breaks in H2O2-treated NP69 and HK1 cells. The cleavage frequency of the ABL gene in H2O2-treated cells was significantly higher than that in the untreated control cells. This is true for both NP69 and HK1 cell lines. These results reaffirm our previous findings which demonstrated that oxidative stress-induced apoptosis resulted in chromosomal breakages in normal nasopharyngeal epithelial and NPC cells . Taken together, our findings are consistent with other studies which discovered that H2O2 induced apoptotic DNA fragmentation. It has been demonstrated that H2O2 induced excision of chromosomal DNA loops mediated by topoisomerase II in U937 leukaemic cells . The production of these HMW DNA fragments (50–100 kb loop-sized DNA fragments) is an initial event of apoptosis . It has also shown that in caspase-3-expressing MCF-7 breast carcinoma cells, H2O2 activated DNA fragmentation with nucleosomal intervals . The fragmentation of nuclear DNA into nucleosomal DNA ladders is another hallmark of apoptosis .
Our sequencing results have confirmed that the IPCR bands were derived from the cleaved ABL gene. All of the breakpoints were mapped within the biochemically defined SAR1 of the ABL gene. SAR1 is located in BCRC, the largest BCR of the ABL gene . MAR/SAR is thought to be one of the common chromatin structures within BCRs. The BCRs of AF9, MLL and AF4 genes have all been found to associate with MAR/SAR, suggesting a role for MAR/SAR in non-homologous recombination (NHR) [76,77,78]. MAR/SAR sequences were found to possess DNA unwinding properties [103, 104]. These properties allow them to facilitate the entry of protein factors that take part in chromosome condensation, apoptosis, transcription and replication [104, 105]. However, these unwinding properties also cause MAR/SAR sequences to be more prone to DNA breakage [103, 104].
In addition, two shift translocations were detected in H2O2-treated NP69 cells. One of the translocated segments was derived from the LHFPL3 gene which locates at chromosome 7.
The other translocated segment was derived from chromosome 5. The disabled homologue 2 (DAB) gene is 1,263,556 bp at the 5′ end of this translocated segment while a gene encoding for hypothetical protein is 22,122 bp at the 3′ end.
The LHFPL3 gene is one of the family members of LHFP-like genes. This gene family consists of six family members. All of the family members have been implicated in human diseases. Members of this family are transmembrane proteins which play important roles in extracellular matrix formation, differentiation and proliferation. Most of them have been associated with tumours . The first member, LHFP on chromosome 13q12, was identified, for the first time, as a translocation partner of HMGIC gene on chromosome 12q15 in human lipoma with t(12;13)(q15;q12). Thus, it was annotated as lipoma HMGIC fusion partner (LHFP) gene . The LHFPL1 gene on chromosome Xq23 has been implicated in liver tumour . The LHFPL2 gene on chromosome 5q14.1 was found to be highly expressed in the novel subgroup of ALL [109, 110] and in patients who succumbed fatally to serous epithelial ovarian cancers (SEOC) . The LHFPL4 gene at 3p25.3 was identified as a novel methylation target specific for cervical cancer . Mutation in the LHFPL5 gene (on chromosome 6p21.31) which is also known as tetraspan membrane protein of hair cell stereocilia (TMHS) gene has been found to cause autosomal recessive nonsyndromic deafness .
The LHFPL3 gene is located on chromosome 7q22.1. Deletions involving chromosome 7q22 are commonly observed in uterine leiomyoma (UL). Four distinct deletion intervals have been identified. One of the microdeletions contains the LHFPL3 gene. The single deleted marker in the microdeletion was mapped within the first intron of the LHFPL3 gene. These findings suggested that the LHFPL3 gene is a candidate tumour suppressor gene (TSG) for UL . Deletion of 7q22 has also been associated with leukaemia. A commonly deleted segment of chromosome 7q22 has been identified in patients with a malignant myeloid disease. The LHFPL3 gene is one of the candidate TSGs residing in this deletion interval . More recently, the alteration of LHFPL3 gene has been suggested to be a hallmark of primary glioblastoma .
Intriguingly, region of microhomology (four nucleotides) was found at the breakpoint junctions. This observation suggested that the shift translocation of the LHFPL3 gene might be mediated by NHEJ DNA repair pathway. Based on the analysis of our sequencing data, we illustrated the potential model for the shift translocation of the LHFPL3 gene (Fig. 10). As proposed by Betti and colleagues (2001), the interaction of the NHEJ DNA repair pathway with apoptosis can act as a mechanism leading to translocation in leukaemia. They found that translocation junctions between the MLL gene and the partner DNA contain regions of microhomology consistent with the operation of NHEJ repair process . In addition, it has been found that cells that survive apoptosis may contain rearranged chromosomes that contribute to leukaemogenesis . Taken together, the findings of ours and of others support the notion that the interaction of the NHEJ DNA repair system with oxidative stress-induced apoptosis may be a possible mechanism leading to chromosome rearrangements in NPC.
We previously proposed a potential model for oxidative stress-induced chromosome rearrangements in NPC involving the AF9 gene . Based on the findings of the ABL gene in the present study and additional findings from the literature, we proposed a revised model (Fig. 11). The revision of this model enables us to further elucidate the potential role of oxidative stress-induced apoptosis in mediating chromosome rearrangements in NPC. We propose that oxidative stress plays an essential role in NPC aetiological factors. These include exposure to nitrosamine, wood dust, formaldehyde and cigarette smoke. EBV infection as well as chronic inflammation of sinonasal tract. Oxidative stress-induced apoptosis is initiated by apoptotic signalling. This includes PS externalisation and MMP loss. The apoptotic signalling may in turn result in the activation of the main effector caspase, caspase-3. Caspase-3 cleaves ICAD that contains two caspase-3 cleavage sites. Subsequently, CAD is being released from its chaperone, ICAD. Chromosomal DNA is cleaved by the activated CAD, presumably at MAR/SAR sites. Double strand breaks are primarily repaired through NHEJ pathway which is prone to cause erroneous DNA repair. Cells that evade apoptosis may harbour chromosome rearrangements such as translocation, deletion, addition and inversion. Repeated exposure to these aetiological factors that provoke oxidative stress may therefore contribute to tumourigenesis of NPC.
In the present study, we only focused on the SAR region of the ABL gene. It is difficult to draw a solid conclusion on the role of MAR/SAR in defining the positions of the chromosome breakages. Therefore, for the future work, comparison in the cleavage frequency between the ABL SAR region and non-SAR region may be carried out. This may allow a further elucidation of the potential role of MAR/SAR in mediating the chromosome breakages and rearrangements in oxidative stress-induced apoptosis.
Our findings demonstrated that oxidative stress-induced apoptosis may be a potential mechanism that leads to chromosome rearrangements in NPC. Our results also suggested that NHEJ system is potentially involved in DNA repair in cells undergoing oxidative stress-induced apoptosis. The interaction between NHEJ DNA repair system and oxidative stress-induced apoptosis may lead to chromosome rearrangements in surviving cells. A revised model for oxidative stress-induced apoptosis in mediating chromosome rearrangement in NPC is proposed.
NP69 normal nasopharyngeal epithelial cell line and HK1 NPC cell line were kindly provided by Prof. Tsao Sai Wah (The University of Hong Kong, Hong Kong, China) and Prof. Lo Kwok Wai (The Chinese University of Hong Kong, Hong Kong, China). NP69 is an immortalised nasopharyngeal epithelial cell line which was established by transfection with SV40 large T oncogene. It retains some characteristics of normal nasopharyngeal epithelial cells and is non-tumourigenic. This cell line may provide potential nasopharyngeal epithelial cell model for investigating mechanisms involved in NPC tumourigenesis . HK1 was derived from a Chinese male patient with recurrent squamous NPC 17 ½ years after radiation therapy .
Hydrogen peroxide (H2O2) was bought from MP Biomedicals, USA. Keratinocyte-SFM medium, RPMI 1640 medium, penicillin, streptomycin, fetal bovine serum and l-glutamine were purchased from GIBCO, Invitrogen, USA. Annexin V-Fluorescein isothiocyanate (FITC) Apoptosis Detection Kit I (BD Pharmingen™) and Flow Cytometry Mitochondrial Membrane Potential Detection Kit were bought from BD™ MitoScreen, Becton–Dickinson Biosciences, USA. Camptothecin (CPT) was purchased from Santa Cruz Biotechnology, CA, USA. Ammonium acetate was bought from Merck, Germany. Chloroform was bought from R&M Chemicals, UK. Phenol and Sodium dodecyl sulfate (SDS) were procured from Amresco, USA. Isoamyl alchohol was purchased from Fluka, Switzerland. Phusion High-Fidelity DNA Polymerase was procured from Finnzymes, Finland. PCR primers were from First Base Laboratories. QIAquick Gel Extraction Kit and QIAquick Nucleotide Removal Kit were bought from QIAGEN, Germany. DNA Polymerase I Large (Klenow) Fragment, restriction enzymes and T4 DNA Ligase were obtained from New England Biolabs (NEB), USA. dNTP mix was purchased from Promega, USA.
NP69 cells were grown in Keratinocyte-SFM medium supplemented with 100 μg/ml streptomycin, 100 U/ml penicillin, 40–50 μg/ml Bovine Pituitary Extract (BPE), 4–5 ng/ml recombinant Epidermal Growth Factor (rEGF) and 2% (v/v) heat-inactivated fetal bovine serum. HK1 cells were cultured in RPMI 1640 medium supplemented with 100 μg/ml streptomycin, 100 U/ml penicillin, 2 mM l-glutamine and 10% (v/v) heat-inactivated fetal bovine serum. Cells were cultured at 37 °C with 5% CO2.
In silico prediction of MAR/SAR
The whole sequence of the ABL gene was retrieved from Ensembl (http://www.ensembl.org/index.html) database [Ensembl:ENSG00000097007]. The location of the experimentally defined MAR/SAR was determined from the previous report . By using DNASTAR software (Lasergene, USA), we predicted the possible MAR/SAR sites within the ABL gene. The prediction of MAR/SAR site was performed by searching MRS which comprises two nucleotide motifs. The first nucleotide motif is an 8 bp degenerate sequence, AATAAYAA, where Y = C or T. The second nucleotide motif is a 16 bp degenerate sequence, AWWRTAANNWWGNNNC, where N = A, C, G or T; R = A or G; W = A or T. One mismatch is allowed in the 16 bp degenerate sequence. The 8 bp degenerate sequence has to be exactly matched. The two sequence elements of the MRS should be found within 200 bp apart. The two sequence elements can be present on either Watson or Crick strand and in either order. The two sequence elements may also be overlapping. When there are more than one motif of either 8 or 16 bp found within a distance of 200 bp, they are considered as a single MRS. In addition, when there is more than one MRS identified within close proximity, they are regarded as a single potential MAR/SAR site .
Phosphatidylserine (PS) externalisation
NP69 cells (1.5 × 105) were plated in 150-mm culture dishes containing 15 ml of complete media. When NP69 cells reached confluency of 30–40% on the third day, NP69 cells were either left untreated or treated with 100 μM of H2O2 for 16 and 24 h. HK1 cells (5.5 × 105) were seeded in 150-mm culture dishes containing 15 ml of complete media. When HK1 cells reached confluency of 60–70% on the fourth day, HK1 cells were incubated with 50 μM of H2O2 for 4 and 8 h. NP69 and HK1 cells treated with camptothecin (CPT) were included as positive controls. After incubation, the cells were collected by using StemPro ACCUTASE Cell Dissociation Reagent. The percentage of apoptotic cells was determined by using Annexin V-FITC Apoptosis Detection Kit I as previously described .
Mitochondrial membrane potential (MMP) loss
NP69 and HK1 cells were treated and collected as described above. The percentage of MMP loss in the harvested cells was determined by using Flow Cytometry Mitochondrial Membrane Potential Detection Kit as previously described .
IPCR detection of chromosome breaks within the ABL gene
Induction of apoptosis in normal nasopharyngeal epithelial and NPC cells with H2O2
NP69 (2 × 104) and HK1 (8 × 104) cells were seeded in 60-mm culture plates containing 4 ml of complete media. When NP69 cells reached confluency of 30–40% on the third day, NP69 cells were either left untreated or treated with 10, 50 or 100 μM for 16 and 24 h. When HK1 cells reached confluency of 60–70% on the fourth day, HK1 cells were either left untreated or treated with 1, 10 or 50 μM of H2O2 for 2, 4, 6 and 8 h.
Genomic DNA extraction
At the end of the indicated exposure times, the used medium was discarded. The cells were washed once with cold 1× phosphate-buffered saline (PBS). Genomic DNA extraction was performed as previously described .
Manipulation of the extracted gDNA for nested IPCR
The extracted gDNA was manipulated as described previously  with minor modifications. Figure 12 shows the manipulation steps. Digestion of the gDNA was performed at 37 °C for 16 h with 100 U of Xba I (RE1 in Fig. 12). The staggered four base pairs (CTAG) 5′ overhang was produced by Xba I digestion. The blunt ends were generated by the apoptotic nuclease such as CAD . After Xba I digestion, both ends of the intact targeted DNA fragment were Xba I sites with staggered overhangs. As for the cleaved targeted DNA fragment, one end was the blunt end produced by the apoptotic nuclease, and the other end was the staggered overhang generated by Xba I. To produce blunt-ended fragments, Klenow fill-in was performed with two μg of DNA template, two units of DNA Polymerase I Large (Klenow) Fragment and 33 μM of dNTP mix at 25 °C for 15 min. Cyclisation was then performed with 2000 U of T4 DNA ligase at 16 °C for 16 h. Ethanol precipitation was carried out with 3 M sodium acetate (NaAc) (one volume), pH 5.2 and ice cold absolute ethanol (2.5 volumes). Seventy percent ethanol was used to wash the DNA pellet. The DNA pellet was then air-dried and dissolved in TE, pH 8.0. The DNA sample was divided into three. The DNA samples of tubes 1, 2 and 3 were subjected to digestion with 10 U of Age I (RE2 in Fig. 12), double digestion with 10 U of each Age I and BsaA I (RE3 in Fig. 12), and double digestion with 10 U of each Age I and EcoR I (RE3 in Fig. 12), respectively. These RE digestions were performed at 37 °C for 16 h. Digestion with Age I was used to linearise the cyclised DNA. Double digestion with Age I and BsaA I or Age I and EcoR I was used to eliminate competition from the intact fragments during IPCR. The double digestion with Age I and BsaA I enabled the detection of DNA cleavages occurred within the amplified region towards the 3′ end. The double digestion with Age I and EcoR I enabled the detection of DNA cleavages occurred within the amplified region towards the 5′ end. According to the manufacturer’s protocol, QIAquick Nucleotide Removal Kit (QIAGEN) was used to purify the digested DNA.
The optical density (O.D.) of the purified DNA sample was measured by using an ultraviolet-visible micro-volume spectrophotometer (ND-1000, NanoDrop, USA). Nested IPCR was performed with 1× of HF buffer (containing 1.5 mM of MgCl2), 0.5 μM of each reverse primer and forward primer, 200 μM of dNTP mix, 0.4 U of Phusion High-Fidelity DNA Polymerase and 200 ng of DNA template. To serve as a negative control, sterile ultrapure water was used to replace the DNA template. Cycle condition used in the first round was: 30 s of 98 °C for 1 cycle (initial denaturation), followed by 30 cycles of 98 °C for 10 s (denaturation), 64 °C for 30 s (annealing), 72 °C for 55 s (extension), followed by 1 cycle of 72 °C for 10 min (final extension). Similar cycle condition was used in the second round of IPCR, except that the extension time was 50 s. Two microlitres of 5-fold diluted IPCR product of the first round was used as DNA template. The primers used in the first round of IPCR were 5’-GGTACCTGGTGTCTGTCTCTATC-3′ (reverse) and 5′-AGAAGGTTTATGGGAGATGG-3′ (forward), whereas the primers used in the second round were 5′-TCTCTCATATCTCAGAGCCTTC-3′ (reverse) and 5′-CTTCAGGAGCTCAGACTTTTAC-3′ (forward). The IPCR assays were done by using a Veriti 96 Well Thermal Cycler (Applied Biosystems, USA).
Agarose gel electrophoresis and DNA sequencing
The PCR products were analysed on 1% agarose gel. The agarose gel electrophoresis was performed at 90 V for 1 h and 30 min. The agarose gel was briefly stained with ethidium bromide (0.5 μg/ml) and destained with distilled water. This was followed by visualisation of the gel on an ultraviolet (UV) transilluminator (Vilber Lourmat). The gel image was captured and analysed using a gel documentation (gel doc) and image analysis system (Syngene). The IPCR bands representing cleaved DNA fragments of the ABL gene were purified by using QIAquick Gel Extraction Kit (QIAGEN) according to the manufacturer’s protocol and sequenced. By blasting the human genome database (Genomic BLAST, https://blast.ncbi.nlm.nih.gov/Blast.cgi), the sequencing data obtained was annotated. To identify the breakpoints of the cleaved fragments, the sequencing data was analysed and aligned with the published ABL gene sequence [Ensembl:ENSG00000097007] by using Seqman DNASTAR software (Lasergene, USA). The positions of DNA breaks identified were compared with the location of the MAR/SAR sequence isolated experimentally in the previous study  and the MRS identified in the present study. A genomic map was constructed to depict the positions of the detected DNA breaks relative to the location of the MAR/SAR.
Quantification of gene cleavage frequency
In each set of IPCR, four to seven IPCR replicates were prepared per cell sample. Each experiment consisted of one to three sets of IPCR. The number of IPCR bands representing the ABL cleaved fragments was counted. Gene cleavage frequency expresses the average number of ABL cleaved fragments detected in three independent experiments.
Experiments were repeated three to five times. The significance of differences in the gene cleavage frequency detected by nested IPCR was evaluated by Student’s t test. Data for IPCR are expressed as mean and standard deviation (SD). Differences were considered statistically significant at p value < 0.05. All statistical tests are two sided.
Acute lymphoblastic leukaemia
Breakpoint cluster region
Chronic myelogenous leukaemia
- DAB :
Disabled homologue 2
DNA double-strand breaks
- H2O2 :
Inhibitor of caspase-activated deoxyribonuclease
Inverse polymerase chain reaction
- LHFPL3 :
Human lipoma HMGIC fusion partner-like 3
Matrix association region/scaffold attachment region
- MLL :
Mixed lineage leukaemia
Mitochondrial membrane potential
Non-homologous end joining
Reactive oxygen species
Shanmugaratnam K. Histological typing of nasopharyngeal carcinoma. IARC Sci Publ. 1978;20:3–12.
Hirayama T. Descriptive and analytical epidemiology of nasopharyngeal cancer. IARC Sci Publ. 1978;20:167–89.
Chang ET, Adami HO. The enigmatic epidemiology of nasopharyngeal carcinoma. Cancer Epidemiol Biomark Prev. 2006;15:1765–77.
Huang DP, Lo KW. Aetiological factors and pathogenesis. In: van Hasselt CA, Gibb AG, editors. Nasopharyngeal carcinoma. Hong Kong: The Chinese University Press; 1999. p. 31–60.
Devi BC, Pisani P, Tang TS, Parkin DM. High incidence of nasopharyngeal carcinoma in native people of Sarawak, Borneo Island. Cancer Epidemiol Biomarkers Prev. 2004;13:482–6.
Henle W, Henle G. Evidence for an etiologic relation of the Epstein-Barr virus to human malignancies. Laryngoscope. 1977;87:467–73.
Dickens P, Srivastava G, Loke SL, Chan CW, Liu YT. Epstein-Barr virus DNA in nasopharyngeal carcinomas from Chinese patients in Hong Kong. J Clin Pathol. 1992;45:396–7.
Feng FT, Cui Q, Liu WS, Guo YM, Feng QS, Chen LZ, et al. A single nucleotide polymorphism in the Epstein-Barr virus genome is strongly associated with a high risk of nasopharyngeal carcinoma. Chin J Cancer. 2015;34:563–72.
Fong LY, Ho JH, Huang DP. Preserved foods as possible cancer hazards: WA rats fed salted fish have mutagenic urine. Int J Cancer. 1979;23:542–6.
Armstrong RW, Armstrong MJ, Yu MC, Henderson BE. Salted fish and inhalants as risk factors for nasopharyngeal carcinoma in Malaysian Chinese. Cancer Res. 1983;43:2967–70.
Yong SK, Ha TC, Yeo MCR, Gaborieau V, McKay JD, Wee J. Associations of lifestyle and diet with the risk of nasopharyngeal carcinoma in Singapore: a case–control study. Chin J Cancer. 2017;36:3.
Armstrong RW, Imrey PB, Lye MS, Armstrong MJ, Yu MC, Sani S. Nasopharyngeal carcinoma in Malaysian Chinese: occupational exposures to particles, formaldehyde and heat. Int J Epidemiol. 2000;29:991–8.
Lo KW, Huang DP. Genetic and epigenetic changes in nasopharyngeal carcinoma. Semin Cancer Biol. 2002;12:451–62.
Guo XC, Scott K, Liu Y, Dean M, David V, Nelson GW, et al. Genetic factors leading to chronic Epstein-Barr virus infection and nasopharyngeal carcinoma in South East China: study design, methods and feasibility. Hum Genomics. 2006;2:365–75.
Chen YJ, Ko JY, Chen PJ, Shu CH, Hsu MT, Tsai SF, et al. Chromosomal aberrations in nasopharyngeal carcinoma analyzed by comparative genomic hybridization. Genes Chromosomes Cancer. 1999;25:169–75.
Hui AB, Lo KW, Leung SF, Teo P, Fung MK, To KF, et al. Detection of recurrent chromosomal gains and losses in primary nasopharyngeal carcinoma by comparative genomic hybridisation. Int J Cancer. 1999;82:498–503.
Fang Y, Guan X, Guo Y, Sham J, Deng M, Liang Q, et al. Analysis of genetic alterations in primary nasopharyngeal carcinoma by comparative genomic hybridization. Genes Chromosomes Cancer. 2001;30:254–60.
Yan J, Fang Y, Liang Q, Huang Y, Zeng Y. Novel chromosomal alterations detected in primary nasopharyngeal carcinoma by comparative genomic hybridization. Chin Med J. 2001;114:418–21.
Li X, Wang E, Zhao YD, Ren JQ, Jin P, Yao KT, et al. Chromosomal imbalances in nasopharyngeal carcinoma: a meta-analysis of comparative genomic hybridization results. J Transl Med. 2006;4:4.
Hui AB, Lo KW, Leung SF, Choi PH, Fong Y, Lee JC, et al. Loss of heterozygosity on the long arm of chromosome 11 in nasopharyngeal carcinoma. Cancer Res. 1996;56:3225–9.
Mutirangura A, Pornthanakasem W, Sriuranpong V, Supiyaphun P, Voravud N. Loss of heterozygosity on chromosome 14 in nasopharyngeal carcinoma. Int J Cancer. 1998;78:153–6.
Lo KW, Teo PM, Hui AB, To KF, Tsang YS, Chan SY, et al. High resolution allelotype of microdissected primary nasopharyngeal carcinoma. Cancer Res. 2000;60:3348–53.
Shao JY, Wang HY, Huang XM, Feng QS, Huang P, Feng BJ, et al. Genome-wide allelotype analysis of sporadic primary nasopharyngeal carcinoma from southern China. Int J Oncol. 2000;17:1267–75.
Huang DP, Lo KW, van Hasselt CA, Woo JK, Choi PH, Leung SF, et al. A region of homozygous deletion on chromosome 9p21-22 in primary nasopharyngeal carcinoma. Cancer Res. 1994;54:4003–6.
Lo KW, Huang DP, Lau KM. p16 gene alterations in nasopharyngeal carcinoma. Cancer Res. 1995;55:2039–43.
Hu LF, Eiriksdottir G, Lebedeva T, Kholodniouk I, Alimov A, Chen F, et al. Loss of heterozygosity on chromosome arm 3p in nasopharyngeal carcinoma. Genes Chromosomes Cancer. 1996;17:118–26.
Hui AB, Or YY, Takano H, Tsang RK, To KF, Guan XY, et al. Array-based comparative genomic hybridization analysis identified cyclin D1 as a target oncogene at 11q13.3 in nasopharyngeal carcinoma. Cancer Res. 2005;65:8125–33.
Lo KW, Cheung ST, Leung SF, van Hasselt A, Tsang YS, Mak KF, et al. Hypermethylation of the p16 gene in nasopharyngeal carcinoma. Cancer Res. 1996;56:2721–5.
Kwong J, Lo KW, To KF, Teo PM, Johnson PJ, Huang DP. Promoter hypermethylation of multiple genes in nasopharyngeal carcinoma. Clin Cancer Res. 2002;8:131–7.
Tan SH, Ida H, Goh BC, Hsieh W, Loh M, Ito Y. Analyses of promoter hypermethylation for RUNX3 and other tumor suppressor genes in nasopharyngeal carcinoma. Anticancer Res. 2006;26:4287–92.
Yanatatsaneejit P, Chalermchai T, Kerekhanjanarong V, Shotelersuk K, Supiyaphun P, Mutirangura A, et al. Promoter hypermethylation of CCNA1, RARRES1, and HRASLS3 in nasopharyngeal carcinoma. Oral Oncol. 2008;44:400–6.
Xiao L, Zhou H, Yan J, Liu Y, Song Y. Changes of telomere length in human nasopharyngeal carcinoma. Hua Xi Yi Ke Da Xue Xue Bao. 1999;30:262–4.
Wen Z, Xiao JY, Guo MH. Telomere shortening in the pathogenesis of nasopharyngeal carcinoma. Di Yi Jun Yi Da Xue Xue Bao. 2002;22:329–30.
Poirier S, Ohshima H, de The G, Hubert A, Bourgade MC, Bartsch H. Volatile nitrosamine levels in common foods from Tunisia, South China and Greenland, high-risk areas for nasopharyngeal carcinoma (NPC). Int J Cancer. 1987;39:293–6.
Cheng YJ, Hildesheim A, Hsu MM, Chen IH, Brinton LA, Levine PH, et al. Cigarette smoking, alcohol consumption and risk of nasopharyngeal carcinoma in Taiwan. Cancer Causes Control. 1999;10:201–7.
Vaughan TL, Stewart PA, Teschke K, Lynch CF, Swanson GM, Lyon JL, et al. Occupational exposure to formaldehyde and wood dust and nasopharyngeal carcinoma. Occup Environ Med. 2000;57:376–84.
Hung SH, Chen PY, Lin HC, Ting J, Chung SD. Association of rhinosinusitis with nasopharyngeal carcinoma: a population-based study. Laryngoscope. 2014;124:1515–20.
Tsou YA, Lin CC, Tai CJ, Tsai MH, Tsai TC, Chen CM. Chronic rhinosinusitis and the risk of nasopharyngeal cancer in a Taiwanese health study. Am J Rhinol Allergy. 2014;28:168–72.
Carnevali S, Petruzzelli S, Longoni B, Vanacore R, Barale R, Cipollini M, et al. Cigarette smoke extract induces oxidative stress and apoptosis in human lung fibroblasts. Am J Physiol Lung Cell Mol Physiol. 2003;284:L955–63.
Uslu C, Taysi S, Bakan N. Lipid peroxidation and antioxidant enzyme activities in experimental maxillary sinusitis. Ann Clin Lab Sci. 2003;33:18–22.
van d V, Postma DS, Timens W, ten Hacken NH. Acute effects of cigarette smoke on inflammation and oxidative stress: a review. Thorax. 2004;59:713–21.
Kum C, Kiral F, Sekkin S, Seyrek K, Boyacioglu M. Effects of xylene and formaldehyde inhalations on oxidative stress in adult and developing rats livers. Exp Anim. 2007;56:35–42.
Pylkkanen L, Stockmann-Juvala H, Alenius H, Husgafvel-Pursiainen K, Savolainen K. Wood dusts induce the production of reactive oxygen species and caspase-3 activity in human bronchial epithelial cells. Toxicology. 2009;262:265–70.
Jones DP. Redefining oxidative stress. Antioxid Redox Signal. 2006;8:1865–79.
Halliwell B, Cross CE. Oxygen-derived species: their relation to human disease and environmental stress. Environ Health Perspect. 1994;102(Suppl 10):5–12.
Wiseman H, Halliwell B. Damage to DNA by reactive oxygen and nitrogen species: role in inflammatory disease and progression to cancer. Biochem J. 1996;313(Pt 1):17–29.
LoPachin RM, Gavin T. Molecular mechanisms of aldehyde toxicity: a chemical perspective. Chem Res Toxicol. 2014;27:1081–91.
Huang YJ, Zhang BB, Ma N, Murata M, Tang AZ, Huang GW. Nitrative and oxidative DNA damage as potential survival biomarkers for nasopharyngeal carcinoma. Med Oncol. 2011;28:377–84.
Klaunig JE, Kamendulis LM. The role of oxidative stress in carcinogenesis. Annu Rev Pharmacol Toxicol. 2003;44:239–67.
Koul M, Kumar A, Deshidi R, Sharma V, Singh RD, Singh J, et al. Cladosporol A triggers apoptosis sensitivity by ROS-mediated autophagic flux in human breast cancer cells. BMC Cell Biol. 2017;18:26.
Sun YY, Wang C, Wang LG, Dai ZB, Yang KB. Arsenic trioxide induces apoptosis and the formation of reactive oxygen species in rat glioma cells. Cellular & Molecular Biology Letters. 2018;23
Park BG, Yoo CI, Kim HT, Kwon CH, Kim YK. Role of mitogen-activated protein kinases in hydrogen peroxide-induced cell death in osteoblastic cells. Toxicology. 2005;215:115–25.
Senkoylu A, Yilmaz A, Ergun MA, Ilhan MN, Simsek A, Altun N, et al. Effect of strontium ranelate on hydrogen peroxide-induced apoptosis of CRL-11372 cells. Biochem Genet. 2008;46:197–205.
Yoon SO, Kim MM, Park SJ, Kim D, Chung J, Chung AS. Selenite suppresses hydrogen peroxide-induced cell apoptosis through inhibition of ASK1/JNK and activation of PI3-K/Akt pathways. FASEB J. 2002;16:111–3.
Ogawa Y, Takahashi T, Kobayashi T, Kariya S, Nishioka A, Mizobuchi H, et al. Mechanism of hydrogen peroxide-induced apoptosis of the human osteosarcoma cell line HS-Os-1. Int J Mol Med. 2003;12:459–63.
Li J, Huang CY, Zheng RL, Cui KR, Li JF. Hydrogen peroxide induces apoptosis in human hepatoma cells and alters cell redox status. Cell Biol Int. 2000;24:9–23.
Kwon D, Choi C, Lee J, Kim KO, Kim JD, Kim SJ, et al. Hydrogen peroxide triggers the expression of Fas/FasL in astrocytoma cell lines and augments apoptosis. J Neuroimmunol. 2001;113:1–9.
Stridh H, Kimland M, Jones DP, Orrenius S, Hampton MB. Cytochrome c release and caspase activation in hydrogen peroxide- and tributyltin-induced apoptosis. FEBS Lett. 1998;429:351–5.
Zunino A, Degan P, Vigo T, Abbondandolo A. Hydrogen peroxide: effects on DNA, chromosomes, cell cycle and apoptosis induction in Fanconi's anemia cell lines. Mutagenesis. 2001;16:283–8.
Wyllie AH, Kerr JF, Currie AR. Cell death: the significance of apoptosis. Int Rev Cytol. 1980;68:251–306.
Coleman ML, Sahai EA, Yeo M, Bosch M, Dewar A, Olson MF. Membrane blebbing during apoptosis results from caspase-mediated activation of ROCK I. Nat Cell Biol. 2001;3:339–45.
Shimizu T, Maeno E, Okada Y. Prerequisite role of persistent cell shrinkage in apoptosis of human epithelial cells. Sheng Li Xue Bao. 2007;59:512–6.
Iguchi K, Hirano K, Hamatake M, Ishida R. Phosphatidylserine induces apoptosis in adherent cells. Apoptosis. 2001;6:263–8.
Brown DG, Sun XM, Cohen GM. Dexamethasone-induced apoptosis involves cleavage of DNA to large fragments prior to internucleosomal fragmentation. J Biol Chem. 1993;268:3037–9.
Oberhammer F, Wilson JW, Dive C, Morris ID, Hickman JA, Wakeling AE, et al. Apoptotic death in epithelial cells: cleavage of DNA to 300 and/or 50 kb fragments prior to or in the absence of internucleosomal fragmentation. EMBO J. 1993;12:3679–84.
Filipski J, Leblanc J, Youdale T, Sikorska M, Walker PR. Periodicity of DNA folding in higher order chromatin structures. EMBO J. 1990;9:1319–27.
Compton MM. A biochemical hallmark of apoptosis: internucleosomal degradation of the genome. Cancer Metastasis Rev. 1992;11:105–19.
Bortner CD, Oldenburg NB, Cidlowski JA. The role of DNA fragmentation in apoptosis. Trends Cell Biol. 1995;5:21–6.
Vaughan AT, Betti CJ, Villalobos MJ. Surviving apoptosis. Apoptosis. 2002;7:173–7.
Betti CJ, Villalobos MJ, Diaz MO, Vaughan AT. Apoptotic triggers initiate translocations within the MLL gene involving the nonhomologous end joining repair system. Cancer Res. 2001;61:4550–5.
Paques F, Haber JE. Multiple pathways of recombination induced by double-strand breaks in Saccharomyces cerevisiae. Microbiol Mol Biol Rev. 1999;63:349–404.
Lieber MR, Gu J, Lu H, Shimazaki N, Tsai AG. Nonhomologous DNA end joining (NHEJ) and chromosomal translocations in humans. Subcell Biochem. 2010;50:279–96.
Weaver DT. What to do at an end: DNA double-strand-break repair. Trends Genet. 1995;11:388–92.
Takata M, Sasaki MS, Sonoda E, Morrison C, Hashimoto M, Utsumi H, et al. Homologous recombination and non-homologous end-joining pathways of DNA double-strand break repair have overlapping roles in the maintenance of chromosomal integrity in vertebrate cells. EMBO J. 1998;17:5497–508.
Weinstock DM, Richardson CA, Elliott B, Jasin M. Modeling oncogenic translocations: distinct roles for double-strand break repair pathways in translocation formation in mammalian cells. DNA Repair (Amst). 2006;5:1065–74.
Strissel PL, Strick R, Tomek RJ, Roe BA, Rowley JD, Zeleznik L. DNA structural properties of AF9 are similar to MLL and could act as recombination hot spots resulting in MLL/AF9 translocations and leukemogenesis. Hum Mol Genet. 2000;9:1671–9.
Strick R, Zhang Y, Emmanuel N, Strissel PL. Common chromatin structures at breakpoint cluster regions may lead to chromosomal translocations found in chronic and acute leukemias. Hum Genet. 2006;119:479–95.
Broeker PL, Super HG, Thirman MJ, Pomykala H, Yonebayashi Y, Tanabe S, et al. Distribution of 11q23 breakpoints within the MLL breakpoint cluster region in de novo acute leukemia and in treatment-related acute myeloid leukemia: correlation with scaffold attachment regions and topoisomerase II consensus binding sites. Blood. 1996;87:1912–22.
Lechardeur D, Xu M, Lukacs GL. Contrasting nuclear dynamics of the caspase-activated DNase (CAD) in dividing and apoptotic cells. J Cell Biol. 2004;167:851–62.
Tan SN, Sim SP, Khoo AS. Potential role of oxidative stress-induced apoptosis in mediating chromosomal rearrangements in nasopharyngeal carcinoma. Cell Biosci. 2016;6:35.
Goff SP, Gilboa E, Witte ON, Baltimore D. Structure of the Abelson murine leukemia virus genome and the homologous cellular gene: studies with cloned viral DNA. Cell. 1980;22:777–85.
Wang JY, Ledley F, Goff S, Lee R, Groner Y, Baltimore D. The mouse c-abl locus: molecular cloning and characterization. Cell. 1984;36:349–56.
Sawyers CL, McLaughlin J, Goga A, Havlik M, Witte O. The nuclear tyrosine kinase c-Abl negatively regulates cell growth. Cell. 1994;77:121–31.
Lugo TG, Pendergast AM, Muller AJ, Witte ON. Tyrosine kinase activity and transformation potency of bcr-abl oncogene products. Science. 1990;247:1079–82.
van Drunen CM, Sewalt RG, Oosterling RW, Weisbeek PJ, Smeekens SC, van Driel R. A bipartite sequence element associated with matrix/scaffold attachment regions. Nucleic Acids Res. 1999;27:2924–30.
Sallmyr A, Fan J, Rassool FV. Genomic instability in myeloid malignancies: increased reactive oxygen species (ROS), DNA double strand breaks (DSBs) and error-prone repair. Cancer Lett. 2008;270:1–9.
Strunnikov AV. One-hit wonders of genomic instability. Cell Div. 2010;5
Giam M, Rancati G. Aneuploidy and chromosomal instability in cancer: a jackpot to chaos. Cell Div. 2015;10
Boichuk SV, Minnebaev MM, Mustafin IG. Key role of mitochondria in apoptosis of lymphocytes. Bull Exp Biol Med. 2001;132:1166–8.
Fadok VA, Voelker DR, Campbell PA, Cohen JJ, Bratton DL, Henson PM. Exposure of phosphatidylserine on the surface of apoptotic lymphocytes triggers specific recognition and removal by macrophages. J Immunol. 1992;148:2207–16.
Anderson GR. Genomic instability in cancer. Curr Sci. 2001;81:501–7.
Shen Z. Genomic instability and cancer: an introduction. J Mol Cell Biol. 2011;3:1–3.
Cockerill PN, Garrard WT. Chromosomal loop anchorage of the kappa immunoglobulin gene occurs next to the enhancer in a region containing topoisomerase II sites. Cell. 1986;44:273–82.
Bartram CR, de Klein A, Hagemeijer A, Grosveld G, Heisterkamp N, Groffen J. Localization of the human c-sis oncogene in Ph1-positive and Ph1-negative chronic myelocytic leukemia by in situ hybridization. Blood. 1984;63:223–5.
Nowell PC, Hungerford DA. Chromosome studies on normal and leukemic human leukocytes. J Natl Cancer Inst. 1960;25:85–109.
Shtivelman E, Gale RP, Dreazen O, Berrebi A, Zaizov R, Kubonishi I, et al. bcr-abl RNA in patients with chronic myelogenous leukemia. Blood. 1987;69:971–3.
Melo JV, Gordon DE, Tuszynski A, Dhut S, Young BD, Goldman JM. Expression of the ABL-BCR fusion gene in Philadelphia-positive acute lymphoblastic leukemia. Blood. 1993;81:2488–91.
Chen SJ, Chen Z, Font MP, d'Auriol L, Larsen CJ, Berger R. Structural alterations of the BCR and ABL genes in Ph1 positive acute leukemias with rearrangements in the BCR gene first intron: further evidence implicating Alu sequences in the chromosome translocation. Nucleic Acids Res. 1989;17:7631–42.
Jiang XY, Trujillo JM, Liang JC. Chromosomal breakpoints within the first intron of the ABL gene are nonrandom in patients with chronic myelogenous leukemia. Blood. 1990;76:597–601.
Li TK, Chen AY, Yu C, Mao Y, Wang H, Liu LF. Activation of topoisomerase II-mediated excision of chromosomal DNA loops during oxidative stress. Genes Dev. 1999;13:1553–60.
Kim DK, Cho ES, Um HD. Caspase-dependent and -independent events in apoptosis induced by hydrogen peroxide. Exp Cell Res. 2000;257:82–8.
Sakahira H, Enari M, Nagata S. Cleavage of CAD inhibitor in CAD activation and DNA degradation during apoptosis. Nature. 1998;391:96–9.
Legault J, Tremblay A, Mirault ME. Preferential localization of DNA damage induced by depurination and bleomycin in a plasmid containing a scaffold-associated region. Biochem Cell Biol. 1997;75:369–75.
Bode J, Benham C, Ernst E, Knopp A, Marschalek R, Strick R, et al. Fatal connections: when DNA ends meet on the nuclear matrix. J Cell Biochem Suppl. 2000;Suppl 35:3–22.
Gohring F, Schwab BL, Nicotera P, Leist M, Fackelmayer FO. The novel SAR-binding domain of scaffold attachment factor A (SAF-A) is a target in apoptotic nuclear breakdown. EMBO J. 1997;16:7361–71.
Ptacek T, Song C, Walker CL, Sell SM. Physical mapping of distinct 7q22 deletions in uterine leiomyoma and analysis of a recently annotated 7q22 candidate gene. Cancer Genet Cytogenet. 2007;174:116–20.
Petit MM, Schoenmakers EF, Huysmans C, Geurts JM, Mandahl N, Van de Ven WJ. LHFP, a novel translocation partner gene of HMGIC in a lipoma, is a member of a new family of LHFP-like genes. Genomics. 1999;57:438–41.
Huang C, Guo J, Liu S, Shan Y, Wu S, Cai Y, et al. Isolation, tissue distribution and prokaryotic expression of a novel human X-linked gene LHFPL1. DNA Seq. 2004;15:299–302.
Yeoh EJ, Ross ME, Shurtleff SA, Williams WK, Patel D, Mahfouz R, et al. Classification, subtype discovery, and prediction of outcome in pediatric acute lymphoblastic leukemia by gene expression profiling. Cancer Cell. 2002;1:133–43.
Schwaenen C, Wessendorf S, Kestler HA, Dohner H, Lichter P, Bentz M. DNA microarray analysis in malignant lymphomas. Ann Hematol. 2003;82:323–32.
Collins Y, Tan DF, Pejovic T, Mor G, Qian F, Rutherford T, et al. Identification of differentially expressed genes in clinically distinct groups of serous ovarian carcinomas using cDNA microarray. Int J Mol Med. 2004;14:43–53.
Wang SS, Smiraglia DJ, Wu YZ, Ghosh S, Rader JS, Cho KR, et al. Identification of novel methylation markers in cervical cancer using restriction landmark genomic scanning. Cancer Res. 2008;68:2489–97.
Kalay E, Li Y, Uzumcu A, Uyguner O, Collin RW, Caylan R, et al. Mutations in the lipoma HMGIC fusion partner-like 5 (LHFPL5) gene cause autosomal recessive nonsyndromic hearing loss. Hum Mutat. 2006;27:633–9.
Curtissa NP, Bonifasa JM, Lauchlea JO, Balkmana JD, Kratza CP, Emerlinga BM, et al. Isolation and analysis of candidate myeloid tumor suppressor genes from a commonly deleted segment of 7q22. Genomics. 2005;85:600–7.
Milinkovic V, Bankovic J, Rakic M, Stankovic T, Skender-Gazibara M, Ruzdijic S, et al. Identification of novel genetic alterations in samples of malignant glioma patients. PLoS One. 2013;8:e82108.
Tsao SW, Wang X, Liu Y, Cheung YC, Feng H, Zheng Z, et al. Establishment of two immortalized nasopharyngeal epithelial cell lines using SV40 large T and HPV16E6/E7 viral oncogenes. Biochim Biophys Acta. 2002;1590:150–8.
Huang DP, Ho JH, Poon YF, Chew EC, Saw D, Lui M, et al. Establishment of a cell line (NPC/HK1) from a differentiated squamous carcinoma of the nasopharynx. Int J Cancer. 1980;26:127–32.
Widlak P, Garrard WT. Discovery, regulation, and action of the major apoptotic nucleases DFF40/CAD and endonuclease G. J Cell Biochem. 2005;94:1078–87.
We would like to thank the Director-General of Health Malaysia for his permission to publish this article and the Director of Institute for Medical Research Malaysia for her support. We gratefully acknowledge Prof. Tsao Sai Wah (The University of Hong Kong, Hong Kong, China) and Prof. Lo Kwok Wai (The Chinese University of Hong Kong, Hong Kong, China) for providing the normal nasopharyngeal epithelial cell line (NP69) and NPC cell line (HK1). We are especially grateful to Dr. Cornelis M. van Drunen (Academic Medical Center, Amsterdam, The Netherlands) for his precious opinions. We express our gratitude to Universiti Malaysia Sarawak for all infrastructure support. This is also another project by the members of Malaysian Nasopharyngeal Carcinoma Group.
This project was funded by Ministry of Health Malaysia (grant number: 06-065), and the consent for publication has been granted.
Availability of data and materials
All data generated or analysed during this study are included in this published article.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Description of exons and introns in the ABL gene. (PDF 65 kb)
Microscopic images of NP69 cells after treatment with H2O2. (PDF 107 kb)
Microscopic images of HK1 cells after treatment with H2O2. (PDF 150 kb)
Description of exons and introns in the LHFPL3 gene. (PDF 56 kb)