Abstract
Having its genome makes the mitochondrion a unique and semiautonomous organelle within cells. Mammalian mitochondrial DNA (mtDNA) is a double-stranded closed circular molecule of about 16 kb coding for 37 genes. Mutations, including deletions in the mitochondrial genome, can culminate in different human diseases. Mapping the deletion junctions suggests that the breakpoints are generally seen at hotspots. ‘9 bp deletion’ (8271–8281), seen in the intergenic region of cytochrome c oxidase II/tRNALys, is the most common mitochondrial deletion. While it is associated with several diseases like myopathy, dystonia, and hepatocellular carcinoma, it has also been used as an evolutionary marker. However, the mechanism responsible for its fragility is unclear. In the current study, we show that Endonuclease G, a mitochondrial nuclease responsible for nonspecific cleavage of nuclear DNA during apoptosis, can induce breaks at sequences associated with ‘9 bp deletion’ when it is present on a plasmid or in the mitochondrial genome. Through a series of in vitro and intracellular studies, we show that Endonuclease G binds to G-quadruplex structures formed at the hotspot and induces DNA breaks. Therefore, we uncover a new role for Endonuclease G in generating mtDNA deletions, which depends on the formation of G4 DNA within the mitochondrial genome. In summary, we identify a novel property of Endonuclease G, besides its role in apoptosis and the recently described ‘elimination of paternal mitochondria during fertilisation.
Research organism: E. coli, Human, Rat
Introduction
The mitochondria are semiautonomous organelles with their genome. The mammalian mitochondrial genome encodes for 13 respiratory chain proteins, 22 tRNAs and 2 rRNAs (Chen and Butow, 2005; Clayton, 1984; Holthöfer et al., 1999). Due to its vital role in oxidative phosphorylation, mitochondrial DNA is frequently exposed to reactive oxygen species (ROS) and is more prone to DNA damage than nuclear DNA, leading to the accumulation of 50 times more mutations, including deletions (Chen et al., 2011; Hudson et al., 1998; Michikawa et al., 1999; Pakendorf and Stoneking, 2005; Yakes and Van Houten, 1997). Mitochondrial mutations have been associated with ageing and several disease conditions like myopathies, dystonia, cancer etc. (Chen et al., 2011; Taylor and Turnbull, 2005). These mutations include point mutations, mismatches, or deletions, which affect the coding of essential proteins involved in oxidative phosphorylation and the respiratory chain (Penta et al., 2001). Mitochondrial mutations, especially large-scale deletions of 200–5000 bp, were associated with the progression of breast, colorectal, renal, and gastric cancers (Bianchi et al., 2001; Carew and Huang, 2002; Modica-Napolitano et al., 2007; Penta et al., 2001). However, the correlation between mtDNA deletions and cancer remains unsettled. The well-studied phenotypes of mtDNA deletions are Kearns-Sayre Syndrome (KSS), Pearson syndrome, and progressive external ophthalmoplegia (PEO; Marni, 2003). Large stretches of deletions seen in different diseased patients suggest fragile sites in the mitochondrial genome similar to the fragile sites (FRA) sites of the nuclear DNA (Richards and Macaulay, 2001).
The FRA regions and other frequent chromosomal translocations have been observed in late replicating sites in the nuclear DNA, which are prone to replication fork stalling due to possible, stable secondary structure formation or inadvertent action of structure-specific nucleases (Dillon et al., 2010; Lin et al., 2013; Raghavan et al., 2004; Sun et al., 2001). Mitochondrial DNA has its replication machinery that is not governed by cell cycle checkpoints. Therefore, it is likely that the formation of a secondary structure may account for observed mitochondrial fragility.
The presence of alternate DNA structures like G-quadruplex, triplex, and cruciform structures in the mitochondrial genome has been reported along with their association with different mitochondrial diseases (Bharti et al., 2014; Damas et al., 2012; Dong et al., 2014; Oliveira et al., 2013). Breakpoint regions reported in the case of PEO patients were in the proximity of G-quadruplex motifs, further suggesting the role of non-B DNA structures in mitochondrial disorders (Bharti et al., 2014). G4 forming motifs were also located adjacent to mitochondrial deletion breakpoints associated with KSS, a clinical subgroup of mitochondrial encephalomyopathies (Van Goethem et al., 2003; Zeviani et al., 1988). Further, a G4 ligand RHPS4 preferentially localises to mitochondria and affects mitochondrial genome content, transcription elongation, and respiratory function (Falabella et al., 2019a; Oliveira et al., 2013). Moreover, in vitro results suggested that a single-base mutation in mtDNA (mt10251) was sufficient to form G4 DNA (Chu et al., 2019). Importantly, ATP-dependent G4 resolving helicase Pif1 could unfold non-B DNA structures and continue the synthesis by mitochondrial DNA polymerases, Pol γ and PrimPol (Butler et al., 2020). However, further studies are required to confirm the presence of these structures in vivo and establish their contribution towards the mechanism of breakage and deletions.
Mitochondrial DNA (mtDNA) 9 bp deletion caused by loss of one copy of the 9 bp repeat sequence (CCCCCTCTA) in the intergenic region of cytochrome c oxidase II (COII)/tRNALys is one of the most common deletions in mitochondria (Redd et al., 1995; Thomas et al., 1998). It has been associated with various diseases, including hepatocellular carcinoma (Jin et al., 2012; Komandur et al., 2011; Krishnan and Turnbull, 2010; Zhuo et al., 2010). The 9 bp deletion is also widely used as a phylogenetic marker to study evolutionary trends and migration of populations (Soodyall et al., 1996; Yao et al., 2000). The occurrence of deletions in such high frequency indicates the fragility of the region; however, its mechanism is unclear.
Out of the large-scale mtDNA deletions reported, one of the well-studied examples is mt4977, and accumulation of it results in different mitochondrial disorders, including cancer (Yusoff et al., 2019). mt4977 eliminates mitochondrial sequences between 8470 and 13,447 bp, which include the removal of genes required for normal oxidative phosphorylation (OXPHOS) (Lee et al., 1994; Yang et al., 1994; Yen et al., 1991). Recently the role of replication slippage and DSB repair has been shown to play a role in the generation of mtDNA4977 in mammalian cells (Phillips et al., 2017). Previous studies have shown that the presence of a 13 bp direct repeat plays a major role in mt4977, as only one of the repeats was retained following joining, indicating the involvement of microhomology-mediated end joining (Samuels et al., 2004; Srivastava and Moraes, 2005). When restriction endonucleases were used to introduce DSBs, there was the accumulation of deletions, which further implicates the imperfect repair and deletion in the mitochondrial genome (Srivastava and Moraes, 2005). A previous study from our laboratory suggested the existence of microhomology-mediated end joining (MMEJ) in mitochondria, which operate in a microhomology-dependent manner, explaining the frequently seen deletions in mitochondrial DNA (Tadi et al., 2016).
Here, we investigate the mechanism of fragility associated with the most common deletion seen in the mitochondrial genome. The putative role of non-B DNA structures and the mechanism of DNA breakage associated with the formation of mtDNA deletions are extensively studied. We report that the region associated with 9 bp deletion can fold into a G-quadruplex DNA structure using various biochemical and ex vivo methods. For the first time, we show that Endonuclease G can specifically bind to the G-quadruplex structure and cleave the mtDNA. Finally, the joining of the broken DNA could lead to the formation of ‘9 bp deletion’ seen in human mitochondria.
Results
Previously, different non-B DNA motifs, including G4 DNA motifs have been reported in the mitochondrial genome (Dahal et al., 2022; Damas et al., 2012; Dong et al., 2014). It was also suggested that alternate forms of DNA may play a role in the generation of deletions in mitochondria (Dahal et al., 2022; Damas et al., 2012; Dong et al., 2014). In silico studies along with biochemical analyses revealed that 5 classical G4 DNA motifs are present in the mitochondrial genome, which could fold into G4 DNA (Cer et al., 2013; Dahal et al., 2022; Figure 1—figure supplement 1A). Since ‘9 bp deletion’ (8271–8281) is the most frequent mitochondrial deletion and occurs next to one of the 5 G4 motifs with the coordinate position 8252–8295 (designated as mitochondrial Region I) (Figure 1A and B; Dahal et al., 2022), we have investigated the putative role of G4 DNA structure in DNA breakage and formation of the mtDNA deletions. Notably, the chosen G-quadruplex motif is also present in two other large-scale deletions associated with mitochondria.
G-quadruplex motifs from Region I of mitochondrial DNA, when present on shorter DNA, can fold into a G-quadruplex structure
To investigate whether mitochondrial Region I can indeed form G-quadruplex structures, oligomeric DNA containing the G-rich region referred to as ‘G1’ and complementary C strand (C1) were synthesised (Figure 1B). Radiolabelled oligomeric DNA was then subjected to electrophoretic mobility gel shift assay in the absence or presence of KCl. Faster mobility of the G-rich strand compared to its complementary C-rich strand suggested the formation of both intra and intermolecular G4 DNA structure in a KCl-dependent manner (Figure 1C). Importantly, the formation of intramolecular G4 DNA was abrogated when G-stretches were mutated (Figure 1—figure supplement 1B, C). Further, circular dichroism (CD) studies revealed that the spectra exhibited the characteristic of a parallel G-quadruplex DNA with absorption maxima of 260–270 nm and absorption minima of 240 nm as opposed to the C-rich strand, which showed the maxima similar to that of a typical single-stranded DNA (Figure 1D).
DMS protection assay was performed to understand the precise base pairing of guanine nucleotides involved in the G4 plates compared to the ones within connecting loops that facilitate DNA strand folding and, thus, G-quadruplex formation (Nambiar et al., 2011; Nambiar et al., 2013). DMS methylates the N7 position of the guanine base, making it susceptible to subsequent piperidine cleavage. Since the N7 position of G residues forming the G4 plate is engaged upon Hoogsteen base pairing, these residues are protected from DMS treatment and piperidine cleavage. When the G-quadruplex sequences from mitochondrial Region I (G1) were subjected to DMS protection assay, following structure formation, we observed that most G residues in this region showed a significant reduction in cleavage, suggesting protection of these nucleotides from DMS. Based on the observed protection of guanine from DMS, we conclude that these guanines are indeed involved in Hoogsteen bonding (Figure 1E, lane 2). The results suggest that the guanine residues from the G stretches of the mitochondrial Region I motif participate in G-quadruplex formation. Based on CD studies and gel shift assays, a parallel, intramolecular G-quadruplex structure was modelled (Figure 1F). Thus, our study reinforces that the Region I of mitochondrial DNA can fold into an intramolecular G-quadruplex (Figure 1E and F).
The G4 motif present in Region I of mtDNA can fold into a G-quadruplex structure when present on a plasmid DNA
In the context of double-stranded DNA, G-quadruplexes can form only when the Hoogsteen base pairing between the guanine residues is stronger than the Watson and Crick H-bonding required for duplex B-DNA formation. Therefore, the G-quadruplex forming ability of the mitochondrial Region I was tested after cloning into a plasmid, pDI1. Interestingly, an inverted repeat sequence (coordinate positions 8295–8363), which could fold into a cruciform DNA, was also present adjacent to the G4 motif. Thus, that was also part of pDI1 (Figure 2—figure supplement 1). Further, to delineate the role of G4 DNA in the plasmid context, a mutant plasmid with a scrambled sequence corresponding to the G4 motif but with intact inverted repeats (pDI2) was also constructed (Figure 2A, Figure 2—figure supplement 1). Previous studies have shown that the formation of G4 DNA can lead to a pause during replication (Kumari et al., 2015; Nambiar et al., 2013; Voineagu et al., 2008; Wells, 2009). To assess whether the G4 DNA formed at Region I of the mitochondrial genome can act as a replication pause, a primer extension assay was performed using pDI1 (Kumari et al., 2015; Figure 2B and C). A sequencing ladder was also generated for mapping the location of the polymerase pause sites on the plasmid (Figure 2C). Results revealed multiple polymerase pause sites, and the intensity of the pauses increased with increasing concentrations of KCl (25, 50 and 100 mM; Figure 2C). Interestingly, there was significantly fewer sites and less pausing observed in the absence of KCl, indicating the absence of the formation of intramolecular G-quadruplex species in these cases (Figure 2C, lane 2). Considering that the formation of cruciform DNA occurs irrespective of the presence of KCl, these results also suggest that the observed pause was due to the formation of G4 DNA, but not cruciform DNA. Besides, the observed polymerase pauses mapped to the G-quadruplex motif (Figure 2C, lanes 7–9). Abrogation of G4 DNA structure by mutation of G4 DNA motif (pDI2) resulted in the disappearance of the observed pause site, further confirming the above results (Figure 2D, lanes 3, 4 and lane 6, 7). It was also observed that the appearance of strong polymerase pause sites were dependent on K+ (Figure 2E, lanes 4–5; F), but not on other monovalent cations such as Li+ (with no or little G-quadruplex stabilising properties) or Na+ (with limited stabilising properties) (Figure 2E, lanes 6–9; F) suggesting that the non-B DNA responsible for polymerase arrest was G4 DNA.
Figure 2. Evaluation of formation of G-quadruplex DNA in Region I of the mitochondrial DNA when cloned into a plasmid.
(A) Schematic showing cloning of the mitochondrial Region I and its mutant to generate plasmids pDI1 and pDI2, respectively. The duplex region containing the G stretches is depicted in blue, while the mutated nucleotides are marked in red. (B) Schematic showing the primer extension across plasmid DNA containing mitochondrial Region I. Positions of primer, VKK11 is indicated. The radiolabeled primer (indicated with an asterisk) binds to one of the strands of the plasmid DNA upon heat denaturation and extends till it encounters a non-B DNA, as it blocks the progression of the polymerase. The products (dotted lines) are then resolved on a denaturing PAGE. (C) The plasmid, pDI1, containing the mitochondrial Region I, was used for primer extension studies using radiolabeled primers for the G-rich (VKK11) in a KCl (25, 50, 100 mM) dependent manner. A sequencing ladder was prepared using VKK11 by the chain termination method of sequencing. Pause sites are shown with dotted rectangles in turquoise. The sequence corresponding to the pause site is boxed, and sequence complementary to G-quadruplex forming motifs are indicated in pink. Sequencing ladder was used as marker. (D) Plasmids, pDI1 and its mutant pDI2, containing mutation in G4 motif, were used for primer extension studies using radiolabelled primer (VKK11). Pause sites are shown with dotted rectangles (turquoise). (E) Evaluation of the effect of different cations (Na+, K+, Li+) on non-B DNA formation at mitochondrial Region I when present on a plasmid (pDI1) by using primer VKK11 in a primer extension assay. The pause sites are indicated in dotted rectangle (turquoise). In panels C-F, ‘M’ denotes the 50 nt ladder. (F) Bar diagram showing the effect of ions in primer extension assay. The values are expressed in PSLU (photo-stimulated luminescence units) representing the extent of pause. Refer also Figure 2—figure supplement 1.
Figure 2—figure supplement 1. Schematic showing cloning of the mitochondrial Region I and its mutant to generate plasmids, pDI1 and pDI2.
Mitochondrial Region I can exist as a G-quadruplex structure within the mitochondrial genome and, when present on a plasmid DNA
Formation of G4 DNA on a plasmid or mitochondrial genome could result in single-strandedness at complementary and adjacent regions corresponding to the G4 motif. Sodium bisulphite modification assay was performed on pDI1 to check the single-strandedness in the complementary C strand of Region I. Sodium bisulphite deaminates cytosine resulting in uracil, when present on single-stranded DNA (Figure 3—figure supplement 1A). This change can be read as a C→T conversion after PCR and sequencing (Figure 3—figure supplement 1A, B). Of the 40 clones sequenced, 38 molecules showed C to T conversion to varying extents in either top or bottom strands (Figure 3A; Figure 3—figure supplement 1C). Interestingly, the cytosines complementary to the G-quadruplex motif exhibited single-stranded nature (Figure 3A), suggesting the formation of a G-quadruplex structure in the complementary G-rich strand. Importantly, ~15% of molecules showed continuous conversion at the region corresponding to the G4 motif. These results suggest that Region I of mitochondrial DNA can fold into a non-B DNA even when cloned into a plasmid, consistent with the result obtained following primer extension.
Figure 3. Evaluation of G-quadruplex formation at mitochondrial Region I.
(A) Bisulphite modification assay on plasmid (pDI1) containing mitochondrial Region I. Vertical bar represents the number of times the respective cytosine in the top and bottom strands of mitochondrial Region I is converted to thymine after deamination when treated with sodium bisulphite, followed by PCR. (B, C) Bisulphite sequencing on the mitochondrial genome for determining the formation of G4 DNA at Region I. Each vertical bar represents the conversion of a cytosine to uracil in the top strand or bottom strand. A total of 59 clones were sequenced from both the top and bottom strands (B). Single strandedness was observed in a DNA fragment of 198 nt containing Region I of the mitochondrial genome following bisulphite modification assay (C). Each row represents cytosines present in a DNA molecule. Each dark circle represents the conversion of cytosine to thymine on the top strand after deamination upon treatment with sodium bisulphite, followed by PCR and DNA sequencing. Of the 59 clones sequenced, the most reactive 25 molecules are shown (C). In all the panels, sequences corresponding to the G-quadruplex forming motif are indicated in a dotted rectangular box (mustard yellow). Refer also Figure 3—figure supplement 1.
Figure 3—figure supplement 1. Schematic showing conversion of cytosine to thymine following bisulfite modification assay.
Further, to investigate the existence of G4 DNA structures in the mitochondrial genome, a sodium bisulphite modification assay was performed after extracting the mitochondrial DNA from human cells (Nalm6) under non-denaturing conditions. Results showed that out of 59 DNA molecules analysed following cloning and sequencing, 52 showed C→T conversions to varying extents (Figure 3B and C; Figure 3—figure supplement 1C). Cumulative analysis of converted cytosines on a fragment of 281 bp revealed single-strandedness at Region I corresponding to the complementary region of G4 motifs. When the mitochondrial genome was analysed, regions adjacent to the G4 motif also showed certain levels of single-strandedness (Figure 3B). However, the C→T conversion rate was significantly lower when upstream and downstream sequences were analysed for conversion. Further, among the 50 clones from the top strand, 14 (28%) had continuous conversions of 11 or more cytosines when a stretch of 20 cytosines modifications was considered on a single DNA molecule level (Figure 3C; Figure 3—figure supplement 1C). This suggests that the formation of G4 DNA at Region I may not occur in all DNA molecules of the mitochondrial genome at a given time. However, based on these results in conjunction with biochemical assays, it is evident that the G-quadruplex motif at Region I can indeed fold into a G-quadruplex structure in mitochondria.
BG4, a G-quadruplex binding antibody, binds to G-quadruplex structures in the mitochondrial genome
The previous studies established that BG4 antibody can bind to G-quadruplexes formed in DNA and RNA (Chambers et al., 2015; Das et al., 2016; Javadekar et al., 2020). To study whether BG4 indeed binds to the mitochondrial genome inside cells, immunofluorescence in HeLa and HEK 293T cells was performed after testing the specificity of the antibody (Javadekar et al., 2020). Mito-tracker Deep Red or Mito-tracker Green FM was used for staining mitochondria, while the nucleus was stained with DAPI. BG4 was stained with Alexa fluor-tagged secondary antibodies. Merging of green and red foci resulting in yellow foci was observed in the case of both HeLa and HEK 293T cells, which revealed localisation of BG4 to mitochondria (Figure 4A, Figure 4—figure supplement 1A). The colocalisation was further analysed using the JaCoP plugin in ImageJ software. In both cases, an image with Coste’s mask is presented, showing the extent of localisation in the form of a white dot after merging red and green foci (Figure 4A–C). Colocalisation analyses show Mander’s coefficient value between 0.25–0.4, irrespective of the stain used as MitoTracker, in the case of HeLa and 293T cells indicating that G-quadruplex structures can be detected in the mitochondrial genome in cellular context (Figure 4B and C). To further confirm that the colocalisation was indeed between the MitoTracker and BG4 antibody, immunofluorescence was performed in Rho(0) cells (human osteosarcoma cells with depleted mtDNA) (Dong et al., 2017; Tan et al., 2015). Results suggested that although BG4 foci were observed in the cytoplasm of Rho(0) cells, they do not seem to colocalize with mitochondria (Figure 4D and E).
Figure 4. Evaluation of G-quadruplex structure in mitochondria of cells.
(A) Representative images of HeLa and HEK293T cells showing colocalisation of BG4, the G4 binding antibody to mitochondria following immunofluorescence assay. The nucleus is stained with DAPI (blue colour), mitochondria with MitoTracker DR (red) and BG4 with Alexa-Fluor 488 (green). A merged image is shown with a merge of red and green, as depicted by Coste’s mask (colocalisation is represented as a white dot). (B, C) Quantitation showing colocalisation of BG4 with MitoTracker indicated as dot plots. The colocalisation was quantified using Mander’s colocalisation coefficient (ImageJ software) analysing a minimum of 100 cells as red over green (B) and green over red (C). (D) Representative image of Rho(0) cells showing localization of BG4 to mitochondria following immunofluorescence as investigated in panel (A). (E) Quantitation showing comparison of colocalization of BG4 between HeLa cells and Rho(0) cells shown as dot plots. The colocalization was quantified using Mander’s colocalization coefficient analyzing a minimum of 50 cells as red over green. Refer also Figure 4—figure supplement 1.
Figure 4—figure supplement 1. Evaluation of existence of G-quadruplex in mitochondria of cells.
To examine whether BG4 can indeed bind to mitochondrial G-quadruplexes, in vitro binding of purified BG4 to oligomers harbouring the G-quadruplex region was performed (Figure 4—figure supplement 1B). Results indicated that the bound complex was seen when G-rich oligomer (G1) was incubated with BG4, while such binding was absent when incubated with C-rich oligomer (C1) (Figure 4—figure supplement 1B).
To investigate whether BG4 can bind to G-quadruplex DNA present in the mitochondrial genome, the mito IP-qPCR strategy was employed (Figure 5A, Figure 5—figure supplement 1). Briefly, following crosslinking, mitochondria were isolated, and the region binding to BG4 was immunoprecipitated. Five regions corresponding to G4 motifs and ten random regions without G4 motifs were selected for analysis by using appropriate PCR primers (Figure 5—figure supplement 1). Results showed that all five G-quadruplex forming regions amplified between 10–18 cycles when qPCR was performed (Figure 5B). In contrast, 9 of 10 random regions did not show any amplification following the BG4 pull-down (Figure 5B). Sequence analysis showed that the control region (CR7) amplified after the BG4 pull-down had high GC content. Besides, when CR7 was analysed for potential non-B DNA formation, we noted the presence of two-plate G4 DNA motif. CR7 may likely fold into a 2-plate G4 DNA, which could explain its binding with BG4 (Figure 5B). However, further studies are required to confirm this hypothesis. Samples with input DNA served as a positive control, and as expected, all the regions were amplified. Samples where no antibody was added, did not result in amplification when G4-specific or random region primers were used, revealing the specificity of BG4 towards G4 DNA (Figure 5B). Similar results were also obtained when BG4 pull-down samples were evaluated using semiquantitative PCR (Figure 5C). However, none of the regions amplified when samples from secondary control were used, including CR7 (Figure 5C). Sequence analysis showed that the control region amplified following BG4 pull down had high GC content and predicted for two plates G4 structure, which requires further analysis (Figure 5B and C).
Figure 5. Evaluation of BG4 binding to G-quadruplex structure in mitochondrial genome.
(A) Schematic showing the experimental strategy used for mito-IP using anti-BG4. Briefly, cells were crosslinked and then mitochondria were isolated and sonicated to obtain the small fragments of mitochondrial DNA. Purified BG4 antibody was used along with protein A/G agarose beads to pull down the BG4 bound regions. (B) BG4 bound mtDNA was purified after reverse crosslinking and used for real-time PCR using primers derived from different regions of the mitochondrial genome, which include 5 G-quadruplex forming regions and 10 random regions. Input DNA served as template control. No antibody control was also used. Bars in blue (first 5) are for G-quadruplex forming regions, while green (last 10) are for random regions. Y-axis depicts threshold Ct value obtained following real time PCR for each primer. Error bar represents mean ± SEM. (C) Agarose gel profile showing the amplification of Input DNA (left panel) and BG4 pull down DNA (right panel). ‘M’ denotes 100 bp ladder. Refer also Figure 5—figure supplement 1.
Figure 5—figure supplement 1. Evaluation of existence of G-quadruplex in mitochondrial DNA.
Thus, our results suggest that Region I and other G4 motifs in the human mitochondrial genome could fold into G-quadruplex DNA. However, based on bisulphite DNA sequencing from Region I, the structure may not be formed in all DNA molecules, and the frequency may be <25%.
Mitochondrial extract induces cleavage at or near the G-quadruplex structure at Region I
To elucidate the mechanism responsible for cleavage at mitochondrial Region I, the involvement of mitochondrial nucleases was considered. To test this, the mitochondrial extract was isolated from the rat tissues (testes and spleen) and used for cleavage assay. The extract was confirmed to be mitochondrial in nature and was tested for cross-contamination by western blotting using the anti-cytochrome C (mitochondrial marker), anti-PCNA (nuclear or cytosolic marker), and anti-Actin (loading control) (Figure 6A). Plasmid harbouring mitochondrial Region I (pDI1) was treated with the mitochondrial extract (0.25, 0.5, 1, and 2 μg) and resolved on an agarose gel following deproteinisation. Interestingly, a concentration-dependent decrease in supercoiled form and increase in the nicked form of plasmid DNA was observed when treated with mitochondrial extracts (Figure 6B, lanes 1–5). Further, there was also a concentration-dependent increase in the linear form of the plasmid DNA, indicating presence of two different nicks at proximity leading to a DSB (Figure 6B, lanes 1–5, C). Although there is a possibility that more than one nick could be present in the pDNA, since the 6 bp mutation in the G4 motif (Figure 6—figure supplement 1A) resulted in abrogation of the conversion of supercoiled to nick DNA, we assume that there will be only a single nick present per strand of the plasmid DNA. To test whether the cleavage observed was due to the presence of G4 DNA, the plasmid bearing mutation in the G4 motif, but bearing intact adjacent inverted repeats (pDI2), was treated with mitochondrial extract and analysed. Interestingly, there was a significant reduction in cleavage when two stretches of G’s were mutated (Figure 6B, lanes 6–10, 6 C; Figure 6—figure supplement 1A), indicating the cleavage was indeed due to the presence of G4 DNA.
Figure 6. Evaluation of mitochondrial extract induced cleavage at G4 DNA formed at Region I of the mitochondrial genome.
(A) Purification and characterization of mitochondrial extracts. Mitochondrial extract was prepared from rat testis and spleen, and its purity was evaluated using specific markers by western blotting. Purity of CE (cytosolic extract) and ME (mitochondrial extract) were determined using antibodies against PCNA (nuclear and cytoplasmic marker), cytochrome C (mitochondrial marker) and Actin (loading control). (B) In vitro nicking assay on a plasmid containing wild type (WT pDNA) mitochondrial Region I (pDI1) and mutant (Mt pDNA) of Region I (pDI2). Both wild type and mutant plasmids were incubated with mitochondrial extract from rat testes (ME) and resolved on a 0.8% agarose gel. Lanes 1 and 6 served as the control without any extract in the reaction, whereas lanes 2–5 and lanes 7–10 are with increasing concentration (0.25, 0.5, 1, and 2 µg) of mitochondrial extract for pDI1 and pDI2, respectively. ‘OC’ is open circular, ‘LIN’ is linear, and ‘SC’ is supercoiled. (C) Quantification showing the efficiency of mitochondrial extracts (Mt extracts) mediated cleavage when a plasmid containing mutant and wild type G4 motif derived from Region I was compared. The values in Y-axis are expressed in PSLU (photo-stimulated luminescence units) representing the extent of pause. (D) Primer extension assay on plasmid containing mitochondrial Region I (pDI1) to determine the position and location of cleavage. pDI1 was incubated with either mitochondrial extract prepared from testes or spleen in an increasing concentration (0.25, 0.5, 1, and 2 μg), reaction products were purified and used for primer extension assay. Pause sites are indicated, which correspond to G4 DNA motif. Sequencing ladder was used as marker. Complementary sequence corresponding to G-quadruplex forming motif is indicated in blue. ‘M’ is 50 nt ladder. (E) Primer extension assay on pDI1 and pDI2 (G4 mutant) following incubation with mitochondrial extract prepared from rat testis (1 μg). ‘M’ is 50 nt ladder. ‘+’ indicates addition of 1 μg mitochondrial extract shown in triplicates. (F) Quantification showing the cleavage efficiency of mitochondrial extracts (testes) when a plasmid containing mutant and wild type G4 motif derived from Region I was incubated. Cleavage intensity is shown in Y-axis. Error bar represents ± SEM. Refer also Figure 6—figure supplements 1 and 2.
Figure 6—figure supplement 1. Mitochondrial extract mediated cleavage assay on plasmids bearing different mutations at G4 DNA motif containing plasmid.
Figure 6—figure supplement 2. DNA sequences that support formation of different conformation of G-quadruplex DNA structures in Region I.
For deciphering whether the position of cleavage observed was indeed at Region I, primer extension assay was carried out on pDI1 and its mutant (pDI2 and pDR4; Figure 6—figure supplement 1A) following cleavage by mitochondrial extracts isolated from testes and spleen. Results revealed two specific nicking sites when incubated with mitochondrial extracts from both testes and spleen. These were mapped to the position corresponding to the G-quadruplex forming region when the DNA sequencing ladder was used as a marker (Figure 6D). Interestingly, mutation of the G4 motif resulted in a significant reduction of cleavage activity, as one of the two cleavage sites disappeared when 2 G stretches were mutated (pDI2) from the G4 motif (Figure 6E and F; Figure 6—figure supplement 1). However, the upper cleavage site remains unchanged with the mutation. This may be due to G4 structure formation in alternate conformations, as Region I can fold into multiple forms (Figure 6—figure supplement 2; Dahal et al., 2022). Interestingly, generation of mutation in third stretch of guanines in addition to the above-described mutation (pDR4) resulted in remarkable reduction in the cleavage activity and led to disappearance of both the bands (Figure 6—figure supplement 1). These results confirm the involvement of G4 DNA in the observed cleavage at Region I.
Endonuclease G mediates structure-selective cleavage at mitochondrial Region I
Among different mitochondrial endonucleases, Endonuclease G was shown to bind preferentially to the guanine-rich strands by a previous study (Ruiz-Carrillo and Renaud, 1987). Since we have observed the formation of G-quadruplex DNA structures with multiple conformations (Figure 6—figure supplement 2; Dahal et al., 2022) at Region I of the mitochondrial genome, we considered the possibility of Endonuclease G as the protein responsible for the observed cleavage. In a different study, specific cleavage at kinked DNA by Endonuclease G has also been reported (Ohsato et al., 2002). Based on these observations, we wondered whether Endonuclease G could cleave G4 DNA present at Region I on pDI1. For this, Endonuclease G and its catalytic dead mutant (ΔH151A) were overexpressed in bacteria, purified, and identity was confirmed by western blotting (Figure 7—figure supplement 1A–C). Purified Endonuclease G was subjected to a plasmid nicking assay and analysed on an agarose gel. Results showed that while the wild-type plasmid (pDI1) showed a conformation change from supercoiled to open circular, mutant plasmid (pDI2) retained its supercoiled form (Figure 7A and B). Besides, the linear form of the plasmid DNA was also observed when a higher concentration (60, 90, 120 ng) of Endonuclease G was used for the cleavage in the case of pDI1 (Figure 7A).
Figure 7. Studies to identify mitochondrial nuclease responsible for cleavage at G4 DNA formed at Region I of the mitochondrial genome.
(A) In vitro nicking assay using purified Endonuclease G on wild type and mutant plasmids containing mitochondrial Region I (pDI1 and pDI2). Both wild type and mutant plasmids were treated with increasing concentration of purified Endonuclease G and resolved on a 0.8% agarose gel. Lanes 1 and 6 served as the control without any protein in the reaction whereas lanes 2–5 and lanes 7–10 are with increasing concentration (30, 60, 90, 120 ng) of purified Endonuclease G. (B) Quantification showing the efficiency of Endonuclease-G-mediated cleavage when a plasmid containing mutant and wild type G4 motif derived from Region I was compared. (C) Cleavage assay was performed on a plasmid containing mitochondrial Region I, pDI1 following incubation with purified Endonuclease G. Primer extension assay was carried out using [γ-32P] radiolabeled VKK11 primer and resolved on 8% denaturing polyacrylamide gel. Lanes 2 and 3 represent 30 and 60 ng of Endonuclease G incubated samples. Lane 1 is without protein, lane 4 is with mitochondrial extract (ME). and lanes 5–8 are A, C, G, T represents the sequencing ladder. ‘M’ is 50 bp ladder. Marked regions represent the specific cleavage products. (D) Cleavage assay was performed on a plasmid containing mitochondrial Region I, pDI1 following incubation either with purified Endonuclease G or mutant Endonuclease G (ΔH151A). Primer extension assay was carried out using γ-32P radiolabeled VKK11 primer and resolved on 8% denaturing polyacrylamide gel. Lane 2 is primer alone, Lanes 3 and 6 are without protein, lanes 4 and 5 are with 30 and 60 ng of Endonuclease G, lanes 7 and 8 are with 30 and 60 ng of mutant Endonuclease G incubated samples. ‘M’ is 50 bp ladder. (E) Quantification showing the efficiency of wild type and mutant Endonuclease G (ΔH151A) mediated cleavage when a plasmid containing wild type G4 motif derived from Region I was compared. (F) Western blotting showing immunodepletion of Endonuclease G from rat testicular mitochondrial extracts. Protein A/G beads were incubated with anti-Endonuclease G and then with the extracts. Actin served as a loading control. (G) Immunodepletion of another endonuclease present in mitochondria, MGME1 from rat testicular mitochondrial extracts as described in panel F. (H) Endonuclease G or MGME1-depleted extract was incubated with pDI1 and used for the primer extension using VKK11 primer. Extract without the addition of antibody served as bead control (lane 4). Lane 3 is no protein control, lanes 5 and 6 corresponds to increasing concentrations of MGME1 depleted extracts, while in lanes 7 and 8, increasing concentrations of Endonuclease G depleted extracts were added. ‘M’ is 50 nt ladder. Cleavage positions are indicated by arrow and boxed (red). (I) Bar diagram depicting quantitation showing the impact of immunoprecipitation of Endonuclease G based on multiple experiments. (J) Reconstitution assay was performed by the addition of purified Endonuclease G following its immunodepletion. Lane 2 represents the primer alone; Lane 3 represents the beads control; lane 4 represents the Endonuclease G depleted extract while lanes 5 and 6 represent the addition of purified Endonuclease G to the depleted extracts (performed in duplicate reaction). (K) Bar diagram representing the cleavage intensity after reconstitution assay as shown in panel J. In panels, C, D, H and J, ‘M’ represents 50 nt ladder. In panels E, I and K, quantitation is based on three biological repeats and data is shown with error bar calculated as mean ± SEM (ns not significant, *p<0.05, **p<0.005, ***p<0.0001). Refer also Figure 7—figure supplement 1.
Figure 7—figure supplement 1. Overexpression, purification and activity assay of different endonucleases.
Interestingly, primer extension studies revealed that Endonuclease G cleaved mitochondrial DNA at multiple sites corresponding to a region of G-quadruplex DNA formation at Region I when it was present on pDI1 (Figure 7C, lanes 1–3, 5–8). Importantly, the primary two cleavage positions were comparable to the one observed when the mitochondrial extract was used for the study (Figure 7C, lane 2–4). However, unlike mitochondrial extract, incubation with purified Endonuclease G resulted in cleavage at multiple positions across the G4 motif. It may be that multiple molecules of purified Endonuclease G may bind to different forms of G4 DNA resulting in cleavage at different positions. In the case of mitochondrial extract, additional proteins may bind to G4 DNA, thus preventing cleavage at multiple sites. More importantly, when a protein with a mutation at the nuclease domain of Endonuclease G (ΔH151A) was used, the cleavage activity was abrogated (Figure 7D and E; Figure 7—figure supplement 1B, C). Further, when endonucleases such as CtIP, FEN1 and RAGs were used for cleavage assay, results showed that, unlike Endonuclease G, none of the purified endonucleases exhibited any significant cleavage (Figure 7—figure supplement 1J), although the presence of CtIP (Tadi et al., 2016) and FEN1 (Kalifa et al., 2009; Klungland and Lindahl, 1997; Tadi et al., 2016) has been shown previously in mitochondria. In all cases, the identity and activity of the endonucleases used were confirmed before use (Figure 7—figure supplement 1D–I).
Immunodepletion of Endonuclease G from mitochondrial extracts prepared from rat testes showed a significant reduction in cleavage efficiency at the mitochondrial G-quadruplex Region I, compared to beads alone control (Figure 7F–I). However, immunodepletion of MGME1, another important nuclease known to be present in mitochondria, did not result in any significant difference in the cleavage efficacy (Figure 7G–I; Phillips et al., 2017; Yang et al., 2018). Importantly, the activity was restored when purified Endonuclease G was added back following its immunodepletion from the extract, confirming the role of Endonuclease G in the cleavage (Figure 7J and K).
Endonuclease G is expressed in mitochondria and colocalises with mitochondrial matrix protein TFAM within the cells
Considering that Endonuclease G is one of the essential mitochondrial nucleases (McDermott-Roe et al., 2011; Misic et al., 2016; Wiehe et al., 2018; Zhou et al., 2016), the presence of Endonuclease G was examined in mitochondria, using mitochondrial specific marker, ‘MitoTracker Deep Red or MitoTracker green FM’ and using immunofluorescence with Endonuclease G antibody in three cell lines of different origin (Figure 8A; Figure 8—figure supplement 1). Merge of MitoTracker with tagged Endonuclease G revealed localization of Endonuclease G in mitochondria in HeLa (human cervical cancer), HEK 293T (human embryonic kidney epithelial) and MEF (mouse embryonic fibroblast) cell lines (Figure 8A and B; Figure 8—figure supplement 1). Further, to check if Endonuclease G is present in the mitochondrial matrix, colocalisation analysis was performed along with TFAM, a transcription factor known to be present in the mitochondrial matrix by conducting immunofluorescence studies in HeLa cells (Figure 8C). Results revealed low but distinct levels of colocalisation of Endonuclease G and TFAM, indicating the presence of Endonuclease G in the mitochondrial matrix (Figure 8C and D). Western blotting further confirmed this after fractionation of the mitochondria (see below).
Figure 8. Evaluation of expression of Endonuclease G within different mammalian cells and Endonuclease-G-mediated cleavage at mtDNA following shRNA mediated knockdown within cells.
(A) Localization of Endonuclease G in mitochondria in different cell lines. Representative images of localization of Endonuclease G to mitochondria in HeLa, MEF and HEK293T cells. FITC-conjugated secondary antibodies were used for detecting Endonuclease G proteins. MtDR is Mitotracker Red, a mitochondrial marker. DAPI is used as nuclear stain. (B) Colocalization analyses of Endonuclease G and Mitotracker signals using JaCoP in ImageJ software based on immunofluorescence studies performed in multiple cell lines (see Figure 8—figure supplement 1). Minimum of 50 cells were used for analysis of colocalization of red and green signals and plotted as the colocalization value of green overlapping red. Y-axis depicts the Mander’s colocalization coefficient value calculated for green over red and plotted in the form of dot plot. The significance was calculated using GraphPad Prism 5.0 with respect to secondary control and shown as mean ± SEM (ns not significant, *p<0.05, **p<0.005, ***p<0.0001). (C) Representative images of colocalization of Endonuclease G to mitochondrial matrix (TFAM) in HeLa cells. Conjugated secondary antibodies were used for detecting Endonuclease G (Alexa Fluor 488) and mitochondrial matrix protein, TFAM (Alexa Fluor 568). DAPI is used as nuclear stain. (D) Colocalization analyses of Endonuclease G and TFAM signals using JaCoP in ImageJ software. Minimum of 50 cells were used for analysis of colocalization of red and green signals and plotted. Y-axis depicts the Mander’s colocalization coefficient value calculated for green over red and plotted in the form of dot plot. Control represents the panel where only one of the primary antibodies was used. (E) Knockdown of Endonuclease G from HeLa cell using PEI mediated transfection. shRNA against Endonuclease G cloned plasmid was used for transfection. Cells were harvested after 48 h and mitochondrial extracts were prepared. Western blotting was performed to confirm the knockdown of Endonuclease G from the HeLa cells. Actin served as a loading control. (F) The knockdown extract was incubated with the plasmid and used for the primer extension using VKK11 primer (‘I’ and ‘II’ are two biological repeats). Extract prepared from the sample transfected with scrambled plasmid served as a control (SCR control). Lanes 3, 4, 7, and 8 serve as scrambled controls while lanes 5, 6, 9, and 10 are for knockdown extracts. I and II represent two independent batches of experiments. ‘M’ is a 50 nt ladder. (G) Bar diagram representing the cleavage intensity of the extracts prepared after transfection with scrambled plasmid and shEndo G plasmid. In panels F and G, a minimum of three biological repeats were performed and the data is shown with the error bar calculated as SEM (ns: not significant, *p<0.05, **p<0.005, ***p<0.0001). Refer also Figure 8—figure supplement 1.
Figure 8—figure supplement 1. Analysis of Endonuclease G localization in mitochondria.
Further, shRNA-mediated knockdown of Endonuclease G was performed in HeLa cells. Knockdown efficiency was assessed by western blotting using Actin as the loading control (Figure 8E). HeLa cells transfected with scrambled plasmid served as transfection control. Cleavage assay using mitochondrial extracts following knockdown of Endonuclease G revealed a significant reduction in cleavage efficiency at Region I, as that of scrambled control, suggesting its role in mitochondria in the context of cells (Figure 8F and G).
Endonuclease G binds to the G-quadruplex forming regions of mitochondria inside cells
To investigate if Endonuclease G binds to G-quadruplex forming regions of mitochondria, G-quadruplex specific antibody, BG4 was used for colocalisation studies along with anti-Endonuclease G (Figure 9A and B; Figure 9—figure supplement 1A). Yellow foci in the merged image suggested colocalisation of BG4 to Endonuclease G. Approximately 100 cells were analysed for colocalisation of anti-BG4 and anti-Endonuclease G using the Mander’s colocalisation coefficient, for red over green, was ~0.32, while green over red was ~0.21 indicating the binding of Endonuclease G to G-quadruplex DNA structure (Figure 9B).
Figure 9. Investigation of binding efficacy of Endonuclease G to G4 DNA at Region I of mitochondrial genome.
(A) Representative image showing colocalization of Endonuclease G with BG4 in HeLa cells. Alexa Fluor 568 and Alexa Fluor 488 conjugated secondary antibodies were used for detection of Endonuclease G and BG4 proteins, respectively. DAPI was used as nuclear stain. (B) The quantitation showing colocalization of Endonuclease G and BG4. The colocalization was quantified using Mander’s colocalization coefficient (ImageJ software) by analyzing a minimum of 100 cells and presented as a dot plot. Red plot represents the overlapping of Endonuclease G over BG4 while green plot represents the overlapping of BG4 over Endonuclease G. (C) Schematic showing the pull-down assay used for evaluation of binding of Endonuclease G present in the rat testicular mitochondrial extracts to the mitochondrial genome. Bound regions were pulled out using anti-Endonuclease G and protein A/G beads. Regions of interest were detected by either semi-quantitative PCR or real-time PCR using appropriate primers. (D) Agarose gel profile showing the amplification through semi-quantitative PCR of Input DNA (upper panel) and Endonuclease G pull down DNA (lower panel). Primers specific to 5 G-quadruplex forming regions (GR1-GR5) and 10 random regions (CR1-CR10) were also used for the amplification. (E) Real-time PCR of 5 G-quadruplex forming regions (blue) and 10 random regions (green) following pull-down assay. Input DNA served as template control. Antibody control served as a negative control. Error bar represents three independent biological repeats. (F) Evaluation of binding of Endonuclease G to different regions of the mitochondrial genome within cells by mito IP. Cells were crosslinked and then mitochondria were isolated. Endonuclease G bound DNA was obtained and was amplified for different regions of mitochondria, as explained in panel E. Graph is plotted for the Ct values obtained following real-time PCR as described above. The error bar represents three independent biological repeats. Refer also Figure 9—figure supplement 1, Figure 2.
Figure 9—figure supplement 1. Binding of Endonuclease G to G-quadruplex regions of the mitochondrial genome.
Figure 9—figure supplement 2. Binding of Endonuclease G to G-quadruplex regions of the mitochondrial genome.
An antibody binding assay was performed to examine the binding of Endonuclease G to mitochondrial G quadruplex regions. Purified Endonuclease G was allowed to bind to mitochondrial DNA isolated from Nalm6 cells (Figure 9—figure supplement 1B) and cross-linked. Sonicated mitochondrial DNA was incubated with anti-Endonuclease G and Protein A/G beads. The resulting bound DNA was reverse crosslinked, purified and used for both real-time and semiquantitative PCR (Figure 9—figure supplement 1C, D). The threshold values presented as a bar diagram suggested that Endonuclease G can bind to all G quadruplex regions (GR1- GR5), including the one present at mitochondrial Region I, while the control regions (CR1-CR10) (not known to form canonical G quadruplex) did not bind to Endonuclease G (Figure 9—figure supplement 1C–E). These results reveal the specific binding of Endonuclease G to mitochondrial G quadruplex regions.
Further, the mitochondrial extract was incubated with the mitochondrial DNA, and a pull-down assay was done using Endonuclease G antibody (Figure 9C). Following the pull-down, beads were separated and loaded for protein analysis and western to confirm the pull-down (Figure 9—figure supplement 2A). The remaining was used for DNA isolation, followed by real-time and semiquantitative PCR. Results suggested the binding of Endonuclease G to G-quadruplex regions (Figure 9C–E). Endonuclease G binding assay inside cells was performed to assess the binding of Endonuclease G to mitochondrial G4 DNA (Figure 9F, Figure 9—figure supplement 1C). The antibody-bound mitochondrial DNA was used for real-time PCR to amplify the G quadruplex regions and control regions. Results showed specific amplification of G quadruplex regions and no amplification of control regions, suggesting the specificity of Endonuclease G to G quadruplex regions within the cells (Figure 9F).
P1 nuclease assay was also performed to determine the binding footprint of Endonuclease G to G-quadruplex forming Region I by incubating the wild type and the mutant oligomers of Region I in the presence and absence of Endonuclease G. Interestingly, binding of Endonuclease G to G-rich strand (G1) partially rescued the sensitivity of P1 nuclease (Figure 9—figure supplement 2B, lanes 5–8). In contrast, in the case of C-rich strands, mutants and random oligomers, such protection was absent (Figure 9—figure supplement 2B).
Although Endonuclease G was thought to be localised within the mitochondrial intermembrane space, its proximity to the mitochondrial matrix was supported by multiple experimental evidence suggesting a direct interaction of Endonuclease G with mtDNA (Duguay and Smiley, 2013; McDermott-Roe et al., 2011; Wiehe et al., 2018). Therefore, as a preliminary investigation, we were interested in testing the conditions in which Endonuclease G could be released to the mitochondrial soluble fraction from the inner mitochondrial membrane. To do this, HeLa cells were exposed to Menadione, a mitochondria-specific ROS inducer (25 µM, for 2 hr). Cells were harvested, and fractionation of mitochondria was performed. Fractionated extracts were equalized and used for western blotting studies (Figure 10A). Results showed that exposure to stress conditions resulted in elevated levels of Endonuclease G in the matrix (Figure 10A and B). However, this warrants further detailed investigation. Thus, we hypothesize that different stress conditions may regulate the release of Endonuclease G to mitochondrial matrix.
Figure 10. Investigation of stress conditions that favor the transport of Endonuclease G to mitochondrial matrix.
(A) Western blot showing the presence of Endonuclease G, Cytochrome C and TFAM in either total or supernatant (Sup) and pellet fraction with or without menadione treatment (25 µM) following sub-fractionation of mitochondrial compartments. Lanes 1–3 are for control samples and Lanes 4–6 are menadione treated samples. (B) Bar graph showing the quantitation of presence of Endonuclease G, Cytochrome C and TFAM in either total mitochondria or supernatant (Sup) and pellet fraction with or without menadione treatment (25 µM). Quantitation is based on three biological repeats and the data is shown with the error bar calculated as SEM. The values in Y-axis are expressed in PSLU (photo-stimulated luminescence units) representing the expression.
Based on the present and previous studies (García-Lepe and Bermúdez-Cruz, 2019; Tadi et al., 2016; Wisnovsky et al., 2018), we propose that under stress conditions, Endonuclease G may be released to the matrix from the inner membrane space where it binds to mitochondrial DNA. This could help the induce DNA breaks by Endonuclease G at Region I of the mitochondrial genome. DNA modifying enzymes such as CtIP, MRN, and Exonuclease G may help in generating 3’ overhangs, following which ligation of broken ends by an unknown mechanism leads to the ‘9 bp deletion’ at the Region I of the mitochondrial genome (Figure 11).
Figure 11. Model depicting mechanism of generation of ‘9 bp deletion’ seen in the mitochondrial genome.
When mitochondria are under stress, Endonuclease G releases into matrix from inner membrane space. In the matrix, Endonuclease G binds and induces cleavage at single-double stranded junctions of G4 DNA as indicated by arrows. Exonuclease action (CtIP/MRE11/Exo G) exposes the direct repeats/microhomology region, which is then paired with the help of unknown proteins, most likely using direct repeats. Ligation of broken ends can result in a ‘9 bp deletion’ at the Region I of the mitochondrial genome.
Discussion
The present study in conjunction with previous studies, suggest the existence of five G-quadruplex DNA structures in the mitochondrial genome (Dahal et al., 2022; Damas et al., 2012; Dong et al., 2014), all with three-plate conformation that follows the general empirical formula for G-quadruplexes d(G3+N1-7G3+N1-7G3+N1-7G3+). G-quadruplexes with three-plate are considered energetically stable (Puig Lombardi and Londoño-Vallejo, 2020). A close analysis of patient breakpoint regions revealed the presence of high-frequency mitochondrial deletion junctions (9 bp deletion) adjacent to one of the G-quadruplex structures at Region I.
The presence of G-quadruplex motifs proximal to mitochondrial deletions has been reported (Falabella et al., 2019a; Falabella et al., 2019b; Oliveira et al., 2013). Computational analysis (Bharti et al., 2014) revealed that G4 sequences proximal to mtDNA deletions in several genetic diseases like Kearns-Sayre syndrome, Pearson marrow-pancreas syndrome, Mitochondrial myopathy, Progressive external ophthalmoplegia, etc. G4 motifs were also located adjacent to mitochondrial deletion breakpoints associated with KSS, a clinical subgroup of mitochondrial encephalomyopathies associated with these diseases (Van Goethem et al., 2003; Zeviani et al., 1988). In the present study, we establish that the region containing ‘9 bp deletion’ can fold into G-quadruplex structures, as shown by EMSA and CD studies in a K+dependent manner. Using bisulphite modification assay, we showed that Region I could fold into a G4 DNA in the mitochondrial genome. Immunofluorescence studies and antibody pull-down assays using G4-specific antibody, BG4, provide further support for the occurrence of G4 quadruplex DNA structures within the cells. Thus, with the combinatorial use of in silico studies, and biophysical, biochemical and ex vivo techniques, we demonstrate the formation of G quadruplexes at the region corresponding to a 9 bp deletion in the mitochondrial genome.
The formation of intramolecular parallel G-quadruplex structure in the mitochondrial genome could make the region fragile as there is an increased probability of replication fork slippage or transcription arrest at this region. The primer extension assay revealed multiple pause sites around the G4 motif, suggesting that polymerase arrest could occur at these structures. However, the deletions seen were precise and occurred in a particular sequence. Thus, an alternate possibility for generating a break at these structures could be due to recognition by structure-selective nucleases followed by cleavage at G4 DNA.
Studies have shown that a few enzymes possess the ability to cleave G-quadruplex structures. The action of RAGs at the G-quadruplex site in BCL2 MBR (Raghavan et al., 2004), MRE11 and DNA2 helicase on either side of the G-quadruplex structure in vitro (Masuda-Sasa et al., 2008) being some of the examples. In the current study, we observed that the structure-selective nuclease activity of mitochondrial endonuclease G could induce cleavage at G4 DNA formed at Region I of the mitochondrial genome using both biochemical (primer extension with purified nucleases, mito-IP with purified Endonuclease G and isolated mitochondria, reconstitution assay with plasmid and purified protein) and ex vivo (immunofluorescence, mito-IP, reconstitution assay in mitochondria) assays. Previous studies have reported various other properties of Endonuclease G, which include the generation of primers during mitochondrial DNA replication by DNA polymerase γ, initiation of genomic inversion in herpes simplex type-1 virus (HSV-1), and its role during apoptosis (Côté et al., 1989; Parrish et al., 2001; Ruiz-Carrillo and Renaud, 1987).
Mammalian Endonuclease G is encoded in the nucleus (Gannavaram et al., 2008; Ohsato et al., 2002; Wiehe et al., 2018; Zhou et al., 2016) and upon translocation to mitochondria, the mitochondrial targeting sequence is cleaved off such that a mature nuclease is released (Li et al., 2001; van Loo et al., 2001; Wiehe et al., 2018). Generally, Endonuclease G is known to reside in the intermembrane space of mitochondria, although a low expression has been reported in the matrix (Duguay and Smiley, 2013; McDermott-Roe et al., 2011; Wiehe et al., 2018). Our preliminary results suggest that Endonuclease G may be further released to the matrix when cells are under stress (Figure 10).
In a recent study, it has been shown that Endonuclease G also cleaved mtDNA in the case of unrepaired oxidative lesions (Wiehe et al., 2018). This was helpful in compensatory mtDNA replication through TWINKLE and mtSSB (Wiehe et al., 2018). However, the binding of Endonuclease G to mtDNA and cleavage might be an event dependent upon various factors, including oxidative stress. Therefore, it is possible that in a typical cellular condition, the level of Endonuclease G in the matrix is low; however, upon stress, it is released into the matrix leading to genome fragility. Hence, the stress-induced release may act as a regulatory mechanism concerning the generation of mitochondrial deletions. A recent study showed the complete degradation property of Endonuclease G at pH 6, which could delineate its role during apoptosis, since apoptotic cells show cytosolic acidification (Schäfer et al., 2004). Interestingly, at physiological pH, the enzyme showed the property of nuclease (converting supercoiled DNA to open circular and linear; Schäfer et al., 2004). CPS-6, a mitochondrial endonuclease G in Caenorhabditis elegans, acts with maternal autophagy and proteasome machinery to promote paternal mitochondrial elimination. It relocates from the intermembrane space of paternal mitochondria to the matrix after fertilisation to degrade mitochondrial DNA (Zhou et al., 2016). Endog-/- mice revealed that Endonuclease G could act as a novel determinant of maladaptive cardiac hypertrophy, which is associated with mitochondrial dysfunction and depletion (McDermott-Roe et al., 2011).
Our study suggests that Endonuclease G can bind to G-quadruplexes and nick the DNA, resulting in single-strand breaks, which can be converted as double-strand breaks during replication or other DNA transactions. However, in low frequency, we could also see DSBs directly induced by Endonuclease G. We observed that the binding was specific to G4 DNA, as we did not observe nonspecific binding with regions where GC content was high in the mitochondrial genome. Previous studies have shown that although Endonuclease G showed random nuclease activity at higher concentrations, at optimal concentrations, it was specific to kink DNA (Ohsato et al., 2002). The first study on Endonuclease G showed preference towards the guanine-rich strands, which falls in line with our study as the formation of G quadruplexes requires a G-rich region (Ruiz-Carrillo and Renaud, 1987). Further, Endonuclease G cleaved mtDNA in unrepaired oxidative lesions (Wiehe et al., 2018). It is likely that binding and cleavage by Endonuclease G on mtDNA be an event dependent on multiple factors, including oxidative stress. Based on our data, it is evident that in the presence of stress, levels of Endonuclease G go up in the matrix, which enables the protein to be present proximal to the mitochondrial genome. The findings from mito IP studies suggest that Endonuclease G only binds to the G-quadruplex forming regions and not to the other areas of mtDNA. This provides additional regulation and ensures the mitochondrial genome stability in most instances.
DNA repair in mitochondria is poorly understood. Base excision repair (BER) is the most well-characterized DNA repair pathway in mitochondria, while it is generally believed that nucleotide excision repair is absent (Akbari et al., 2008; de Souza-Pinto et al., 2009; Liu et al., 2008). Homologous recombination repair helps repair DSBs and, thus stability of the genome (Dahal et al., 2018), though classical NHEJ was undetectable (Tadi et al., 2016). In contrast, microhomology-mediated end joining uses flanking direct repeats to fix DSBs and delete the sequences between the repeats (Tadi et al., 2016). Interestingly, the G quadruplex forming Region I is close to the deletion junction 8271:8281 (Redd et al., 1995; Yao et al., 2000) also had a 9 bp direct repeat sequence at the flank site.
In a recent study, a replication slippage-mediated mechanism was attributed to a common 4977 deletion where in frequent replication fork stalling at the junction of the common deletion was observed (Phillips et al., 2017). While this is a replication-dependent process, what we observed may not be immediately connected to replication. However, it is possible that replication can indeed facilitate the formation of G-quadruplexes.
A balance between the formation and resolution of G4 DNA within the mitochondrial genome is essential. The mitochondrial genes are transcribed as a single polycistronic mRNA (Clayton, 1984; Sanchez et al., 2011) and therefore undergo a high transcription rate, which may promote G-quadruplex formation. Besides, random ‘breathing’ of DNA (a process of transient melting of the duplex structure due to thermal fluctuations) may also occur in the GC-rich sequences of the mitochondrial DNA promoting G-quadruplex formation (Dornberger et al., 1999; Jose et al., 2009). Once the structure is formed, G4 resolvases also play a critical role in its resolution. In the nuclear DNA context, resolvases like WRN and BLM helicase are considered as the factors that negatively regulate G4 DNA (Giri et al., 2011; Paeschke et al., 2013). However, such a role for these nucleases in mitochondria is not reported. TWINKLE helicase is known to unwind mitochondrial duplex DNA and regulate replication, although its role in structural resolution is yet to be studied (Phillips et al., 2017).
Hence, our study suggests a novel mechanism that explains fragility in mitochondria, dependent on the structure-selective nuclease activity of Endonuclease G on G-quadruplex DNA. The process could be controlled by cellular stress conditions, which may regulate the release of Endonuclease G to the mitochondrial matrix; however, this requires further extensive investigation. Although we have focused on a single example of ‘9 bp deletion’ often seen in mitochondria, it is very much possible that it may be a general mechanism and may help in explaining several other deletions seen in mitochondria-associated human disorders.
Materials and methods
Enzymes, chemicals, and reagents
Chemicals and reagents used in the present study were purchased from Millipore-Sigma (USA), SRL (India), Himedia (India), and Amresco (USA). Restriction enzymes and other DNA-modifying enzymes were obtained from New England Biolabs (USA). Culture media was from Sera Laboratory International Limited (UK), Lonza (UK). Fetal bovine serum and PenStrep were from Gibco BRL (USA). Radioisotope-labelled nucleotides were from BRIT (India). Antibodies were purchased from Abcam (UK), Santa Cruz Biotechnology (USA), BD (USA), Cell Signaling Technology (USA) and Calbiochem (USA). Oligomeric DNA was purchased from Juniper Life Sciences (India) and Medauxin (India).
Oligomeric DNA
The oligomers used in this study are listed in Table 1. Oligomeric DNA was purified on 12–15% denaturing PAGE when needed (Nambiar and Raghavan, 2012; Raghavan et al., 2005b; Tadi et al., 2016). The details of oligomers used for different experiments are provided in appropriate sections. A more simplified nomenclature is used for certain oligomers while explaining the results of respective experiments.
Table 1. Oligomers used in the study.
Oligomer Name | Sequence | Region |
---|---|---|
VKK11 | 5’- GCTGTGTCGACTACTACGGTCAATGCTCTG –3’ | GR1 |
VKK12 | 5’- CTGAGGTCGACTGGGTGATGAGGAATAGTG - 3’ | |
RBK46 | 5’- TAATCAACACCCTCCTAGCC –3’ | GR2 |
VKK14 | 5’- GATAGTGTCGACGGCTCATGGTAGGGGTAA –3’ | |
SD10 | 5'- TTCGCTGACGCCATAAAACT –3' | GR3 |
SD11 | 5'- ATCAGGGCGTAGTTTGA –3' | |
SD12 | 5'- GCTCACAAGAACTGCTAA –3' | GR4 |
SD13 | 5'- TGGATGCGACAATGGAT –3' | |
SD14 | 5'- TCTTGCACGAAACGGGAT –3' | GR5 |
SD15 | 5'- TAGGATGAGGATGGATAGT –3' | |
RBK15 | 5’- CTACTCCTGCTCGCATCTGC –3’ | CR2 |
RBK16 | 5’- GAAGGTGGTGTTGAGGTTGC –3’ | |
RBK17 | 5’- GCATTGTTCGTTACATGGTCC –3’ | CR3 |
RBK18 | 5’- GTGGAAGCGGATGAGTAAGAAG –3’ | |
RBK19 | 5’- CTCACCACTACAATCTTCCTAG –3’ | CR4 |
RBK20 | 5’- CAAAGATGGTAGAGTAGATGACG –3’ | |
RBK21 | 5’- CTAACCATCTTCTCCTTACACCTAG –3’ | CR5 |
RBK22 | 5’- GTTTGCTAATACAATGCCAGTCAGG –3’ | |
RBK23 | 5’- CGAAGGTGGATTTAGCAGTAAACTG –3’ | CR6 |
RBK24 | 5’- CGGTACTATATCTATTGCGCCAGG –3’ | |
RBK41 | 5’- GTATCATCAACTGATGAGCAAG –3’ | CR1 |
SD2 | 5’- TCAGCAAACCCTGATGAA –3’ | CR7 |
SD3 | 5'- CACTCTACTCTCAGTTTACT –3' | |
SD4 | 5'- ACATCGAATACGCCGCA –3' | CR8 |
SD5 | 5'- AGTTGGTCGTAGCGGAATCG –3' | |
SD6 | 5'- TAGGGTTTATCGTGTGAG –3' | CR9 |
SD7 | 5'- AGTGTGGCGAGTCAGCT –3' | |
SD8 | 5'- TACTCACTCTCACTGCCCAA –3' | CR10 |
SD9 | 5'- TGTTTGTCGTAGGCAGAT-3' | |
VKK21 | 5'- GGATCCATGCGGGCGCTGCGG –3' | |
VKK22 | 5'- GCGGCCGCTCACTTACTGCCCG –3' | |
DI12 | 5'- GCAAACCACAGTTTCATGCCCATC –3' | |
DI13 | 5'- GCCTATAATCACTGCGCCCGCTC –3' | |
VKK1 | 5’- CCCGTATTTACCCTATAGCACCCCCTCTACCCCC –3’ | C1 |
VKK2 | 5’- GGGGGTAGAGGGGGTGCTATAGGGTAAATACGGG –3’ | G1 |
VKK5 | 5’- GTCAGTAGAGGGGGTGCTATAGGGTAAATACGGG –3’ | M1 |
VKK6 | 5’- GTCAGTAGAGAATGTGCTATAGGGTAAATACGGG –3’ | M2 |
VKK7 | 5’- GTCAGTAGAGGGGGTGCTATATCATAAATACGGG –3’ | M3 |
SD 54 | 5'- GGCCAGGGCCCCGCGGTCGAAGCCACTGCC-3' | |
SD 57 | 5'- TCACCTGGCCGCCGCCGCCAACCAC-3' | |
DK27 | 5’-TGGGCTCTAGAGGACATAGAGTAAGTGCT-3’ | |
DK28 | 5’-AGCACTTACTCTATGTCCTCTAGAGCCCA-3’ | |
KD14 | 5’-CAAGCTCGAAATTAACCCTCAC-3’ | |
KD13 | 5’-CCCAGTCACGACGTTGTAAAAC-3’ | |
DI8 | 5’-CTTACAGTGGGCTCTAGAGGGGGTAGATAATAT GCTATAGGGTAAATACTCACTAAAAATCTTTGAA-ATAGGG –3’ |
|
DI9 | 5’-CTAAAAATCTTTGAAATAGGGTGAGTATTTA CCCTATAGCATATTATCTACCCCCTCTAGAGCCCA-CTGTAAG –3’ |
Cell lines
HeLa (human cervical cancer), HCT116 (human colon cancer), MEF (mouse embryonic fibroblast) and HEK 293T (human embryonic kidney epithelial cell line) were purchased from National Centre for Cell Science, Pune, India. Nalm6 and Reh cells were from Dr M. R. Lieber (USA), and Rho(0) cells were from J. Neuzil (Australia). The identity of all these cell lines was confirmed by performing STR profiling (DNA labs India, Hyderabad, India). All the cell lines were found to be free of mycoplasma contamination. The Mycoplasma test was conducted in the laboratory using the MYCOseq kit (Thermo Fisher Scientific, Waltham, Massachusetts) and further validated by DNA labs India, Hyderabad, India. Cells were cultured in RPMI1640 or MEM medium supplemented with 10–15% fetal bovine serum (FBS), 100 μg/ml Penicillin, and 100 μg/ml streptomycin and incubated at 37 °C in a humidified atmosphere containing 5% CO2 as described before (Ghosh et al., 2022; Kumari et al., 2021). Rho(0) cells of B16 cells origin were cultured in DMEM medium supplemented with 15% fetal bovine serum (FBS), 100 μg/ml Penicillin, 100 μg/ml streptomycin, 50 mg/ml Uridine and 1 mM Sodium pyruvate and incubated as described above (Dong et al., 2017; Tan et al., 2015).
Animals
Male Wistar rats (Rattus norvegicus) 4–6 weeks old were purchased from the Central animal facility, Indian Institute of Science (IISc), Bangalore, India, and maintained as per the guidelines of the animal ethical committee in accordance with Indian National Law on animal care and use (CAF/Ethics/526/2016).
shRNA
shRNA used in the study (TRCN0000039643) was purchased from the Resource Center at Division of Biological Science, IISc, Bangalore (funded by DBT: BT/PR4982/AGR/36/718/2012) Indian Institute of Science, Bangalore (India).
Plasmid constructs
Plasmid construction for mitochondrial G-quadruplex assays
Mitochondrial DNA isolated from Nalm6 cells was amplified using primers specific for Region I (VKK11 and VKK12). The PCR products were purified, blunted and ligated into the EcoRV site of pBlueScript SK(+) to obtain the wild-type plasmid for mitochondrial Region I (pDI1). Mutant plasmids for mitochondrial Region I were generated by PCR-mediated site-directed mutagenesis using pDI1 as a template described before (Raghavan et al., 2005c) using the primers VKK11, DI8, VKK12 and DI9 followed by cloning into pBlueScript SK(+). The resulting plasmid, pDI2, had a mutation in two G stretches of Region I. A second mutant plasmid, pDR4 with mutations in three G stretches, was generated by site-directed mutagenesis using pDI2 as a template using the primers DK27, KD13, KD14, and DK28 followed by cloning into pBlueScript SK(+).
Construction of Endonuclease G containing vector for purification of the protein
The endonuclease G coding sequence was amplified using primers VKK21 and VKK22, which were designed with BamHI and NotI restriction sites, respectively, at flanking regions. PCR product was gel purified and cloned into the EcoRV site of pBS SK +cloning vector. For the generation of pRBK2, the coding DNA sequence was digested from the cloning vector using BamHI and NotI enzymes and cloned into the same site (BamHI and NotI) of pET28a+. Mutant plasmid for Endonuclease G was generated by PCR mediated site-directed mutagenesis as described before (Raghavan et al., 2005c) using the primers VKK21, SD54, VKK22 and SD57 followed by cloning into BamHI and NotI site of pRBK2. The resulting plasmid, pDR1, had a mutation in the nuclease domain of Endonuclease G (Histidine to Alanine at position 141).
Plasmid isolation and purification
For each plasmid, after transformation, Escherichia coli was cultured in 500 ml of Luria broth (HiMedia, USA) for 18 h at 37 °C. Isolation of plasmid DNA was performed by the standard alkaline lysis method. It was then purified and precipitated, as described previously (Sambrook et al., 1989). The pellet was dissolved in TE (pH 8.0).
Isolation of mitochondrial DNA
Mitochondria were isolated from Nalm6 cells as described before (Dahal et al., 2018; Tadi et al., 2016). Cells were first homogenized in the buffer [70 mM sucrose, 200 mM mannitol, 1 mM EDTA, and 10 mM HEPES (pH 7.4), 0.5% BSA per g tissue] in a Dounce-type homogenizer (20–30 strokes). Homogenate was centrifuged (3000 rpm, 10 min at 4 °C) to remove nuclei and cellular debris from the supernatant that contains the cytosol and mitochondria. The supernatant fraction was centrifuged at 12,000 rpm (30 min at 4 °C) to pellet the mitochondria. The mitochondrial pellet was washed in suspension buffer [10 mM of Tris–HCl (pH 6.7), 0.15 mM of MgCl2, 0.25 mM of sucrose, 1 mM PMSF and 1 mM DTT] twice, and the mitochondrial pellet was incubated in the buffer containing 75 mM NaCl, 50 mM EDTA (pH 8.0), 1% SDS and 0.5 mg/ml Proteinase K to remove denatured proteins. The resulting lysate was deproteinised by phenol-chloroform extraction. The mitochondrial DNA was precipitated with 3 M Sodium acetate (pH 5.2) and chilled with 100% ethanol. The purity of the mitochondrial DNA preparation was evaluated by PCR amplification with nuclear and mitochondrial DNA-specific primers.
Preparation of mitochondrial extracts
Mitochondrial protein extracts were prepared following the isolation of mitochondria based on differential centrifugation as described (Dahal et al., 2018; Tadi et al., 2016). Testes and spleen from 4- to 6-week-old male Wistar rats were minced on ice and lysed in mitochondrial lysis buffer [50 mM of Tris–HCl (pH 7.5), 100 mM of NaCl, 10 mM of MgCl2, 0.2% Triton X-100, 2 mM of EGTA, 2 mM of EDTA, 1 mM of DTT and 10% glycerol] by mixing (30 min at 4 °C) along with protease inhibitors PMSF (1 mM), aprotinin (1 μg/ml), pepstatin (1 μg/ml) and leupeptin (1 μg/ml). The mitochondrial extract was centrifuged at 12,000 rpm (5 min), and the supernatant (mitochondrial fraction) was aliquoted, snap-frozen and stored at −80 °C till use.
Mitochondrial protein extracts were also prepared from rat tissues as per the manufacturer’s instructions using a mitochondrial extraction kit (Imgenex, USA). The purity of the mitochondrial extracts prepared was checked by immunoblotting using the nuclear and mitochondrial specific markers.
Fractionation of mitochondria
A total of 4×107 HeLa cells were treated with 25 µM Menadione and incubated for 2 hr at 37 °C. The mitochondrial pellet was first sonicated after incubation in hypotonic buffer 10 mM HEPES-KOH pH7.6, 100 mM NaCl, 10 mM Magnesium acetate and 1 mM phenylmethyl sulfonyl fluoride (PMSF) for 15 min. The same buffer containing 500 mM NaCl was added and incubated for another 15 min and then sonicated for 10 min at 37% duty per cycle. The samples were then separated into membrane and soluble fractions by ultracentrifugation at 100,000 g (TLA 110 rotor) for 30 min at 4 °C in OptimaTM TLX table-top ultracentrifuge (Beckman-Coulter) (Sinha et al., 2010; Suzuki et al., 2002). The supernatant fraction was transferred to a new tube. To the pellet, 1 X PBS was added, and both fractions were used for western blotting analysis as described before (Chiruvella et al., 2012).
5’ end labelling of oligomeric substrates
The 5’ end labelling of oligomeric substrates was done using γ[32P]-ATP in the presence of 1 U of T4 Polynucleotide kinase in a buffer containing 20 mM Tris-acetate [pH 7.9], 10 mM magnesium acetate, 50 mM potassium acetate and 1 mM DTT at 37 °C for 1 h. The radiolabelled oligomers were purified on a Sephadex G25 column and stored at –20 °C until further use (Kumar et al., 2010; Sharma et al., 2011; Tadi et al., 2016).
Gel mobility shift assays
Gel mobility shift assays were performed using 5’ radiolabelled oligomeric DNA substrates (4 nM) predicted to form G-quadruplex structures, its complementary oligos and the mutants by incubating in a buffer containing 10 mM Tris-HCl (pH 8.0), 1 mM EDTA, with or without 100 mM KCl at 37 °C and resolved on 15% native polyacrylamide gels (Dahal et al., 2022; Nambiar et al., 2011; Nambiar et al., 2013). In conditions where the effect of KCl was studied, 100 mM KCl was added while preparing gels and the running buffer (1 X TBE). Electrophoresis was carried out at 150 V at RT. The dried gels were exposed, and signals were detected by phosphorImager FLA9000 (Fuji, Japan).
For BG4 binding assay, 5’ radiolabelled oligomeric DNA substrates (4 nM) predicted to form G-quadruplex structures and its complementary oligomer was incubated in a buffer containing 10 mM Tris-HCl (pH 8.0), 1 mM EDTA, with 100 mM KCl at 37 °C for 1 h (Das et al., 2016; Javadekar et al., 2020). It was then incubated along with increasing concentrations of purified BG4 protein (100, 200, 400, and 800 ng) in a buffer containing 250 mM Tris (pH 8.0), 1 mM EDTA, 0.5% Triton X-100, 500 μg/ml BSA and 20 mM DTT for 1 hr at 4 °C. The complex was then resolved on 5% native polyacrylamide gels in the presence of 100 mM KCl in gel and the running buffer (1 X TBE). Electrophoresis was carried out at 100 V at 4 °C. The dried gel was exposed to a PI cassette, and the radioactive signals were detected by phosphorImager FLA9000 (Fuji, Japan).
Circular dichroism (CD)
Circular dichroism studies for analysis of mitochondrial G-quadruplex formation were performed with oligonucleotide substrates (2 µM) corresponding to G-rich sequences and complementary C-rich sequences in a buffer containing 10 mM Tris-HCl (pH 8.0) and 1 mM EDTA. The experiments were set up either in the presence or absence of 100 mM KCl. The spectra were recorded between wavelengths 200 and 300 nm (5 cycles, scan speed of 50 nm/s, RT) using JASCO J-810 spectropolarimeter (Nambiar and Raghavan, 2012; Raghavan et al., 2005a) and analysed using SpectraManager (JASCO J-810 spectropolarimeter).
DMS protection assay
The radiolabelled oligomer VKK19 was incubated with Dimethyl Sulphate (1:250 dilution) in buffer containing 10 mM Tris-HCl (pH 8.0) and 1 mM EDTA either in the presence or absence of 100 mM KCl, at RT for 15 min (Kumari et al., 2019; Nambiar et al., 2013). An equal volume of 10% piperidine was added, and the reactions were incubated at 90 °C for 30 min. The reaction mixtures were diluted 2-fold with double distilled water and vacuum dried in a Speed Vac concentrator. The resulting pellet was again washed with water and vacuum dried. This procedure was repeated thrice. The pellet was finally resuspended in 10 μl TE buffer (10 mM Tris-HCl, (pH 8.0) and 1 mM EDTA), mixed with formamide containing dye, and resolved on a 15% denaturing PAGE.
Polymerase stop assays
Plasmid substrates pDI1, pDI2 (mitochondrial Region I), Vent Exo(-) polymerase was used for primer extension assays (Kumari et al., 2015; Nambiar et al., 2013) in a buffer containing 20 mM Tris-HCl (pH 8.8), 10 mM (NH4)2SO4,10 mM KCl, 2 mM MgSO4, 0.1% Triton X-100 with additional supplementation with 75 mM KCl, LiCl or NaCl where required. After the addition of 200 μM dNTPs, 0.2 U of Vent Exo(-) and 5’ radiolabelled primers (VKK11, VKK12, DI12), the reaction was carried out in a one-step PCR-mediated primer extension assay for mitochondrial Region I (95 °C for 10 min, 55–65°C for 3 min as specified below for each primer and 75 °C for 20 min, as a single cycle) or multi-cycle PCR extension (Denaturation at 95 °C for 5 min, followed by 95 °C for 45 s, 55 °C for 45 s and 72 °C for 45 s for 15 cycles and a final extension at 72 °C for 5 min). The reaction products were loaded on 8% denaturing PAGE, and the products were visualised as described above.
Endonuclease cleavage assay
Wild-type plasmid (pDI1) or mutant plasmid (pDI2) or mitochondrial DNA was incubated with mitochondrial extract or purified Endonucleases (CtIP, FEN1, RAG, Endonuclease G [wildtype or mutant]) in a buffer containing 25 mM MOPS (pH 7.0), 30 mM KCl, 30 mM potassium glutamate and 5 mM MgCl2 at 37 °C for 1 hr. In the control reaction, only DNA and buffer were added. Following termination of the reaction, purified samples were loaded on 0.8% agarose gel or subjected to primer extension assay using a radiolabeled primer (VKK11 or VKK12) and resolved on 8% denaturing PAGE.
Preparation of sequencing ladders
Sequencing ladders for the primers VKK11 and VKK12 were prepared using the cycle sequencing method with a dNTP:ddNTP ratio of 1:20, 1:40, 1:30 and 1:10 for C, T, A and G ladders, respectively. The dNTP:ddNTP mix were separately provided in reaction mixtures containing 10 nM plasmid template, 0.5 μM 5’ radiolabelled primers, 20 mM Tris-HCl (pH 8.8), 10 mM (NH4)2SO4, 10 mM KCl, 2 mM MgSO4, 0.1% Triton X-100 and 1 U of Vent Exo (-) polymerase. The PCR was carried out using the following conditions: Step 1 at 95 °C for 30 s, 60 °C for 30 s and 72 °C for 1 min (25 cycles), followed by Step 2 at 95 °C for 30 s and 72 °C for 2 min (10 cycles).
Sodium bisulphite modification assay
pDI1, plasmid containing mitochondrial Region I or mitochondrial DNA isolated from Nalm6 cells, were treated with sodium bisulfite as described earlier (Raghavan et al., 2004; Raghavan et al., 2006). Briefly, approximately 2 μg of mitochondrial DNA was treated with 12.5 μl of 20 mM hydroquinone and 458 μl of 2.5 M sodium bisulfite (pH 5.2) at 37 °C for 16 hr. The bisulphite-treated DNA was purified using Wizard DNA Clean-Up Kit (Promega, Madison, WI) and desulfonated by treating with 0.3 M NaOH (15 min at 37 °C). The DNA was ethanol-precipitated and resuspended in a 20 μl TE buffer. The mitochondrial Region I was PCR amplified, resolved on 1% agarose gels, purified, and TA cloned. Clones were sequenced to analyse for conversions. The experiment was repeated three independent times, and the cumulative data is presented.
Overexpression and purification of endonucleases
Endonuclease G and its mutant
For protein expression, E. coli Rosetta cells were transformed with pRBK2 or pDR1, grown until OD reaches 1.0 and then induced with 1 mM IPTG at 16 °C for 16 hr. Cells were harvested, and the extract was prepared using extraction buffer (20 mM Tris-HCl (pH 8.0), 0.5 M KCl, 20 mM imidazole (pH 7.0), 20 mM Mercaptoethanol, 10% glycerol, 0.2% Tween 20, 1 mM PMSF) by sonication followed by loading on to Ni-IDA column (Macherey-Nagel, Germany). Endonuclease G or the mutant was eluted in gradient imidazole concentration (100–500 mM) and pure fractions were pooled, and the identity of the protein was confirmed by immunoblotting and used for the assays.
RAGs
MBP cRAGs (RAG1, amino acids 384–1040; RAG2, amino acids 1–383) were purified using a method as described previously (Naik et al., 2010; Raghavan et al., 2005b). Briefly, 293T cells were transfected with 10 μg of plasmid by the calcium phosphate method. After 48 hr of transfection, cells were harvested, and proteins were purified using an amylose resin column (New England Biolabs). Fractions were eluted and checked by silver staining. The activity was and studied by nicking assay on standard recombination signal sequence substrate (AKN1/2; Nishana and Raghavan, 2012).
FEN1
FEN1 was purified using the expression plasmid pET-FEN1 CH described previously (Greene et al., 1999). Briefly, culture was grown until OD reached ~0.5, induced with 0.5 mM IPTG at 37 °C for 3 hr. Cells were harvested, lysed, and the lysate was loaded onto the Ni-NTA column. FEN1 was eluted in increasing imidazole concentrations (100–500 mM), and pure fractions were pooled, and the identity of the protein was confirmed by immunoblotting and used for the assays.
CtIP
CtIP was purified using the expression plasmid pET14b-CtIP as described previously (Yu and Baer, 2000). Briefly, culture was grown till OD reached 0.6, induced with 0.5 mM IPTG at 37 °C for 4 hr. Cells were harvested, lysed, and the lysate was loaded on to Ni-NTA column. FEN1 was eluted in increasing imidazole concentrations (100–800 mM), pure fractions were pooled, and the protein’s identity, was confirmed by immunoblotting and used for the assays.
Overexpression and purification of BG4
The plasmid expressing BG4 protein, pSANG10-3F-BG4, was a gift from Shankar Balasubramanian (Addgene plasmid # 55756). The plasmid was transformed into E. coli, BL21 (DE3), and the culture was expanded by incubating at 30 °C till the O.D. reached up to 0.6 (Biffi et al., 2013; Das et al., 2016; Javadekar et al., 2020). The cells were then induced with 1 mM IPTG for 16 hr at 16 °C, harvested, and resuspended in lysis buffer (20 mM Tris-HCl [pH 8.0], 50 mM NaCl, 5% glycerol, 1% Triton X-100 and 1 mM PMSF). The cells were lysed by sonication, centrifuged, and the supernatant was then loaded onto a Ni-NTA chromatography column (Novagen, USA). BG4 was eluted using increasing concentrations of Imidazole (100–400 mM). BG4 enriched fractions were dialysed against dialysis buffer (PBS containing 0.05% Triton X-100, 1 mM Mercaptoethanol, 5% glycerol and 0.1 mM PMSF) overnight at 4 °C. The identity of the protein was confirmed by immunoblotting using an anti-FLAG antibody (Calbiochem, USA) (Kumari et al., 2019). Activity assay for each batch of BG4 was checked by performing a binding assay as described previously (Javadekar et al., 2020).
Immunoblot analysis
For immunoblotting analysis, approximately 20–30 μg protein was resolved on 8–10% SDS-PAGE (Chiruvella et al., 2008; Kumar et al., 2010; Thomas et al., 2016). Following electrophoresis, proteins were transferred to the PVDF membrane (Millipore, USA), blocked using 5% non-fat milk or BSA in PBS with 0.1% Tween-20. Proteins were detected with appropriate primary antibodies against Cytochrome-C (ab90529), PCNA (SC56), Tubulin (SC5286), CtIP (SC22838), FEN1 (SC13051), RAG1 (SC5599), Endonuclease G (#4969), TFAM (ab176558) and appropriate secondary antibodies as per standard protocol. The blots were developed using a chemiluminescent substrate (Immobilon western, Millipore, USA) and scanned by a gel documentation system (LAS 3000, FUJI, Japan).
Immunoprecipitation (IP)
IP experiments were performed with modifications (Chiruvella et al., 2012; Sharma et al., 2015; Totaro et al., 2011; Yeretssian et al., 2011). Protein A agarose beads (Sigma) were activated by immersing in water and incubated in IP buffer (300 mM NaCl, 20 mM Tris-HCl (pH 8.0), 0.1% NP40, 2 mM EDTA, 2 mM EGTA, and 10% glycerol) for 30 min on ice following which the beads were conjugated with the appropriate antibody at 4 °C for overnight to generate antibody-bead conjugate. The antibody-bead conjugates were then separated by centrifugation and incubated with rat tissue mitochondrial extracts at 4 °C overnight. Then the conjugate bound to the target proteins was separated and washed. Immunoblot analysis confirmed the protein depletion in the resulting supernatant and quantified using Multi Gauge (V3.0). This immunodepleted extract was used for cleavage assay.
Knockdown of Endonuclease G within cells
HeLa cells (10X105) were seeded in a culture petridish. 10 μg of Endonuclease G shRNA plasmid was transfected using linear/branched PEI (Polyethylenimine) polymer (Sigma, 1 mg/ml) (Kumari et al., 2021; Longo et al., 2013; Raymond et al., 2011; Srivastava et al., 2012). Post transfection (48 hr), cells were harvested, and mitochondrial extracts were prepared by differential centrifugation as described above. The knockdown was confirmed using western blotting. Cleavage assay on pDI1 following primer extension was performed using the knockdown extracts. The scrambled plasmid was also used for transfection, which acted as a control for the experiment.
Binding of Endonuclease G to mitochondrial DNA
5 µg of mitochondrial DNA was incubated with either 5 µg purified Endonuclease G or 50 µg of mitochondrial extract in a 100 µl reaction in PBS. DNA-protein interaction was crosslinked by using 0.1% formaldehyde at 37 °C for 10 min, following quenching the reaction with glycine (125 mM) for 5 min. The complex was sonicated (VC 750, Ultrasonic processor) for 10 min with the pulse of 10 sec on/ 10 sec off and incubated at –80 °C for 8 h. 75 ng of Endonuclease G antibody (CST: 4969) was added and incubated at 4 °C (end to end rotation). A total of 15 µl protein A/G beads (Santa cruz, sc-2003) were added to the chromatin-antibody complex and incubated for 2 h. The supernatant was removed entirely, and beads were washed and eluted with high salt immune complex wash buffer (50 mM HEPES (pH 7.9), 500 mM NaCl, 1 mM EDTA, 0.1% SDS, 1% Triton X-100, 0.1% deoxycholate) and incubated at 55 °C for 2 h. DNA protein interactions were reverse crosslinked by incubating the samples overnight at 65 °C. DNA was purified with phenol: chloroform extraction followed by ethanol precipitation (Carey et al., 2009). Endonuclease G pulled-down mitochondrial DNA and was subjected to semiquantitative and real-time PCR of different mitochondrial regions (five G-quadruplex forming regions and ten random control regions). Three independent reactions were performed, and the bar diagram represented threshold values.
Mitochondrial IP (mito IP)
4X107 Nalm6 cells were crosslinked with 1% formaldehyde at 37 oC for 15 min and quenched with 100 µl/ml of 1.375 M glycine. Mitochondria were then isolated from the cells as described above and lysed using buffer (5 mM PIPES, 85 mM KCl, 0.5% NP40). Lysed samples were then sonicated (Diagenode Bioruptor, Belgium) with 30 s on/ 45 s off pulse for 30 cycles. Sonicated samples were stored at −80 °C for at least 8 hr (Carey et al., 2009). For immunoprecipitation, the samples were centrifuged at high speed for 15 min, and the supernatant was divided for input, secondary control and experimental. Relevant antibodies (Endonuclease G or BG4) were added to the sample and allowed to bind for 8–10 hr. Protein A/G-agarose beads (Sigma, USA) were added to the models and incubated for 2 hr. Samples were washed and eluted using high salt wash buffer (50 mM HEPES (pH 7.9), 500 mM NaCl, 1 mM EDTA, 0.1% SDS, 1% Triton X-100, 0.1% deoxycholate) in repeated cycles of centrifugation and then reverse cross-linked by incubating overnight at 65 °C. Finally, DNA was purified following phenol: chloroform extraction and precipitation. The purified DNA was used for real time PCR amplification of different mitochondrial regions that either support formation of G-quadruplex structure or do not support any secondary structure formation.
Immunofluorescence
Immunofluorescence studies were performed as described before (Dahal et al., 2018; Kumari et al., 2019; Ray et al., 2020; Tadi et al., 2016). Approximately 50,000 cells (HeLa, HEK293T, Rho(0), or MEF) were grown in the media for 24 hr, as described before. Cells were grown in chamber slides in MEM medium supplemented with 10% FBS and 1% penicillin-streptomycin (Sigma) for 24 hr. The cells were stained with either 100 ng/ml of MitoTracker Red 580 (Invitrogen) or 100 ng/ml MitoTracker Green FM (Invitrogen) at 37 °C in a CO2 incubator for 30 min (Tadi et al., 2016). Cells were washed twice with 1 X PBS, fixed in 2% paraformaldehyde (20 min) and permeabilised with 0.1% Triton X-100 (5 min) at room temperature. BSA (0.1%) was used for blocking (30 min), and subsequently the cells were incubated with the appropriate primary antibody at room temperature (4 hr). Appropriate FITC or Alexa Fluor conjugated secondary antibodies were used to detect the signal. After washing, the cells were stained with DAPI, mounted with DABCO (Sigma) and imaged under a confocal laser scanning microscope (Zeiss LSM 880 with ×63 magnifications or Olympus FLUOVIEW FV3000 with ×100 magnifications). The images were processed using either Zen Lite or FV31S-SW software.
Statistical analysis
Every experiment was repeated multiple times with at least three biological repeats and presented. Based on the data obtained from repeats, Student’s t-test (two-tailed) was performed, and the error bars were calculated for the bar diagrams representing were the fold change value using GraphPad Prism ver 5.0. Error bars are shown depicting mean ± SEM (ns: not significant, *p<0.05, **p<0.005, ***p<0.0001). Colocalisation analysis in immunofluorescence assay was done using JaCoP in ImageJ software. The values were plotted in GraphPad Prism 5.0, and the significance was calculated using the same. The values were obtained for at least 50–100 cells per sample.
Acknowledgements
We thank Bibha Choudhary, Mridula Nambiar, Urbi Roy, Nitu Kumari, and other members of the SCR laboratory for critical reading and comments on the manuscript. We thank Monica Pandey, Rupa Kumari, Mahesh Hegde, Meghana Manjunath and Supriya Vartak for their technical help. We also thank J Neuzil (Australia), A Agarwal (India) and N Bhatraju for support with Rho(0) cells, S Balasubramanian (UK) for pSANG10-3F-BG4 plasmid, Michael A Resnick (USA) for FEN1 plasmid, Richard Baer (USA) for CtIP-pET14b plasmid, Patrick Swanson (USA) for MBP-RAG constructs, Kumar Somasundaram (India) for pIRES2-EGFP and Arun Kumar (India) for pcDNA3.1RFP. We thank Confocal, FACS and Central Animal Facilities, shRNA core facility of IISc, for their help. This work was supported by grants from CSIR (37 (1692)/17/EMR-11), DAE (21/01/2016-BRNS/35074), DBT-COE (BT/PR/3458/COE/34/33/2015), IISc-DBT partnership program [BT/PR27952-INF/22/212/2018] to SCR. SD is supported by a Senior Research Fellowship (SRF) from IISc, and HS is supported by a Junior Research Fellowship (JRF) from CSIR.
Funding Statement
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Contributor Information
Sathees C Raghavan, Email: sathees@iisc.ac.in.
Wolf-Dietrich Heyer, University of California, Davis, United States.
Jessica K Tyler, Weill Cornell Medicine, United States.
Funding Information
This paper was supported by the following grants:
Council of Scientific and Industrial Research, India 37(1692)/17/EMR-11 to Sathees C Raghavan.
Department of Atomic Energy, Government of India 21/01/2016-BRNS/35074 to Sathees C Raghavan.
Department of Biotechnology, Ministry of Science and Technology, India BT/PR/3458/COE/34/33/2015 to Sathees C Raghavan.
IISc-DBT partnership programme BT/PR27952-INF/22/212/2018 to Sathees C Raghavan.
Indian Institute of Science to Sumedha Dahal, Humaira Siddiqua.
Additional information
Competing interests
No competing interests declared.
No competing interests declared.
Author contributions
Conceptualization, Formal analysis, Validation, Investigation, Methodology, Writing - original draft, Writing – review and editing.
Formal analysis, Validation, Investigation, Methodology, Writing - original draft.
Investigation, Methodology.
Investigation, Methodology.
Conceptualization, Resources, Supervision, Funding acquisition, Validation, Investigation, Methodology, Project administration, Writing – review and editing.
Investigation, Methodology.
Conceptualization, Resources, Supervision, Funding acquisition, Validation, Investigation, Methodology, Project administration, Writing – review and editing.
Ethics
This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the Indian National Law on animal care and use. All of the animals were handled according to approved institutional animal care and use committee protocols (CAF-SOP) of the Indian Institute of Science, Bangalore. The protocol was approved by the Committee on the Ethics of Animal Experiments of the Central Animal Facility (CAF/Ethics/526/2016).
Additional files
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files. Source data for each data is provided along with the figures.
References
- Akbari M, Visnes T, Krokan HE, Otterlei M. Mitochondrial base excision repair of uracil and AP sites takes place by single-nucleotide insertion and long-patch DNA synthesis. DNA Repair. 2008;7:605–616. doi: 10.1016/j.dnarep.2008.01.002. [DOI] [PubMed] [Google Scholar]
- Bharti SK, Sommers JA, Zhou J, Kaplan DL, Spelbrink JN, Mergny JL, Brosh RM. DNA sequences proximal to human mitochondrial DNA deletion breakpoints prevalent in human disease form G-quadruplexes, a class of DNA structures inefficiently unwound by the mitochondrial replicative twinkle helicase. The Journal of Biological Chemistry. 2014;289:29975–29993. doi: 10.1074/jbc.M114.567073. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Bianchi NO, Bianchi MS, Richard SM. Mitochondrial genome instability in human cancers. Mutation Research. 2001;488:9–23. doi: 10.1016/s1383-5742(00)00063-6. [DOI] [PubMed] [Google Scholar]
- Biffi G, Tannahill D, McCafferty J, Balasubramanian S. Quantitative visualization of DNA G-quadruplex structures in human cells. Nature Chemistry. 2013;5:182–186. doi: 10.1038/nchem.1548. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Butler TJ, Estep KN, Sommers JA, Maul RW, Moore AZ, Bandinelli S, Cucca F, Tuke MA, Wood AR, Bharti SK, Bogenhagen DF, Yakubovskaya E, Garcia-Diaz M, Guilliam TA, Byrd AK, Raney KD, Doherty AJ, Ferrucci L, Schlessinger D, Ding J, Brosh RM. Mitochondrial genetic variation is enriched in G-quadruplex regions that stall DNA synthesis in vitro. Human Molecular Genetics. 2020;29:1292–1309. doi: 10.1093/hmg/ddaa043. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Carew JS, Huang P. Mitochondrial defects in cancer. Molecular Cancer. 2002;1:e9. doi: 10.1186/1476-4598-1-9. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Carey MF, Peterson CL, Smale ST. Chromatin immunoprecipitation (CHIP) Cold Spring Harbor Protocols. 2009;2009:db. doi: 10.1101/pdb.prot5279. [DOI] [PubMed] [Google Scholar]
- Cer RZ, Donohue DE, Mudunuri US, Temiz NA, Loss MA, Starner NJ, Halusa GN, Volfovsky N, Yi M, Luke BT, Bacolla A, Collins JR, Stephens RM. Non-B DB v2.0: a database of predicted non-B DNA-forming motifs and its associated tools. Nucleic Acids Research. 2013;41:D94–D100. doi: 10.1093/nar/gks955. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Chambers VS, Marsico G, Boutell JM, Di Antonio M, Smith GP, Balasubramanian S. High-throughput sequencing of DNA G-quadruplex structures in the human genome. Nature Biotechnology. 2015;33:877–881. doi: 10.1038/nbt.3295. [DOI] [PubMed] [Google Scholar]
- Chen XJ, Butow RA. The organization and inheritance of the mitochondrial genome. Nature Reviews. Genetics. 2005;6:815–825. doi: 10.1038/nrg1708. [DOI] [PubMed] [Google Scholar]
- Chen T, He J, Shen L, Fang H, Nie H, Jin T, Wei X, Xin Y, Jiang Y, Li H, Chen G, Lu J, Bai Y. The mitochondrial DNA 4,977-bp deletion and its implication in copy number alteration in colorectal cancer. BMC Medical Genetics. 2011;12:8. doi: 10.1186/1471-2350-12-8. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Chiruvella KK, Kari V, Choudhary B, Nambiar M, Ghanta RG, Raghavan SC. Methyl angolensate, a natural tetranortriterpenoid induces intrinsic apoptotic pathway in leukemic cells. FEBS Letters. 2008;582:4066–4076. doi: 10.1016/j.febslet.2008.11.001. [DOI] [PubMed] [Google Scholar]
- Chiruvella KK, Sebastian R, Sharma S, Karande AA, Choudhary B, Raghavan SC. Time-Dependent predominance of nonhomologous DNA end-joining pathways during embryonic development in mice. Journal of Molecular Biology. 2012;417:197–211. doi: 10.1016/j.jmb.2012.01.029. [DOI] [PubMed] [Google Scholar]
- Chu IT, Wu CC, Chang TC. G-Quadruplex formation by single-base mutation or deletion of mitochondrial DNA sequences. Biochimica et Biophysica Acta. General Subjects. 2019;1863:418–425. doi: 10.1016/j.bbagen.2018.11.009. [DOI] [PubMed] [Google Scholar]
- Clayton DA. Transcription of the mammalian mitochondrial genome. Annual Review of Biochemistry. 1984;53:573–594. doi: 10.1146/annurev.bi.53.070184.003041. [DOI] [PubMed] [Google Scholar]
- Côté J, Renaud J, Ruiz-Carrillo A. Recognition of (dG) N. (dC) N sequences by endonuclease G. characterization of the calf thymus nuclease. The Journal of Biological Chemistry. 1989;264:3301–3310. [PubMed] [Google Scholar]
- Dahal S, Dubey S, Raghavan SC. Homologous recombination-mediated repair of DNA double-strand breaks operates in mammalian mitochondria. Cellular and Molecular Life Sciences. 2018;75:1641–1655. doi: 10.1007/s00018-017-2702-y. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Dahal S, Siddiqua H, Katapadi VK, Iyer D, Raghavan SC. Characterization of G4 DNA formation in mitochondrial DNA and their potential role in mitochondrial genome instability. The FEBS Journal. 2022;289:163–182. doi: 10.1111/febs.16113. [DOI] [PubMed] [Google Scholar]
- Damas J, Carneiro J, Gonçalves J, Stewart JB, Samuels DC, Amorim A, Pereira F. Mitochondrial DNA deletions are associated with non-B DNA conformations. Nucleic Acids Research. 2012;40:7606–7621. doi: 10.1093/nar/gks500. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Das K, Srivastava M, Raghavan SC. GNG motifs can replace a GGG stretch during G-quadruplex formation in a context dependent manner. PLOS ONE. 2016;11:e0158794. doi: 10.1371/journal.pone.0158794. [DOI] [PMC free article] [PubMed] [Google Scholar]
- de Souza-Pinto NC, Mason PA, Hashiguchi K, Weissman L, Tian J, Guay D, Lebel M, Stevnsner TV, Rasmussen LJ, Bohr VA. Novel DNA mismatch-repair activity involving YB-1 in human mitochondria. DNA Repair. 2009;8:704–719. doi: 10.1016/j.dnarep.2009.01.021. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Dillon LW, Burrow AA, Wang YH. Dna instability at chromosomal fragile sites in cancer. Current Genomics. 2010;11:326–337. doi: 10.2174/138920210791616699. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Dong DW, Pereira F, Barrett SP, Kolesar JE, Cao K, Damas J, Yatsunyk LA, Johnson FB, Kaufman BA. Association of G-quadruplex forming sequences with human mtDNA deletion breakpoints. BMC Genomics. 2014;15:677. doi: 10.1186/1471-2164-15-677. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Dong LF, Kovarova J, Bajzikova M, Bezawork-Geleta A, Svec D, Endaya B, Sachaphibulkij K, Coelho AR, Sebkova N, Ruzickova A. Horizontal transfer of whole mitochondria restores tumorigenic potential in mitochondrial DNA-deficient cancer cells. eLife. 2017;6:e22187. doi: 10.7554/eLife.22187. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Dornberger U, Leijon M, Fritzsche H. High base pair opening rates in tracts of GC base pairs. The Journal of Biological Chemistry. 1999;274:6957–6962. doi: 10.1074/jbc.274.11.6957. [DOI] [PubMed] [Google Scholar]
- Duguay BA, Smiley JR. Mitochondrial nucleases ENDOG and EXOG participate in mitochondrial DNA depletion initiated by herpes simplex virus 1 UL12.5. Journal of Virology. 2013;87:11787–11797. doi: 10.1128/JVI.02306-13. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Falabella M, Fernandez RJ, Johnson FB, Kaufman BA. Potential roles for G-quadruplexes in mitochondria. Current Medicinal Chemistry. 2019a;26:2918–2932. doi: 10.2174/0929867325666180228165527. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Falabella M, Kolesar JE, Wallace C, de Jesus D, Sun L, Taguchi YV, Wang C, Wang T, Xiang IM, Alder JK, Maheshan R, Horne W, Turek-Herman J, Pagano PJ, St Croix CM, Sondheimer N, Yatsunyk LA, Johnson FB, Kaufman BA. G-quadruplex dynamics contribute to regulation of mitochondrial gene expression. Scientific Reports. 2019b;9:5605. doi: 10.1038/s41598-019-41464-y. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Gannavaram S, Vedvyas C, Debrabant A. Conservation of the pro-apoptotic nuclease activity of endonuclease G in unicellular trypanosomatid parasites. Journal of Cell Science. 2008;121:99–109. doi: 10.1242/jcs.014050. [DOI] [PubMed] [Google Scholar]
- García-Lepe UO, Bermúdez-Cruz RM. Mitochondrial Genome Maintenance: Damage and Repair Pathways. UK: IntechOpen; 2019. [Google Scholar]
- Ghosh D, Kumari S, Raghavan SC. Depletion of rnaseh2 activity leads to accumulation of DNA double-strand breaks and reduced cellular survivability in T cell leukemia. Journal of Molecular Biology. 2022;434:167617. doi: 10.1016/j.jmb.2022.167617. [DOI] [PubMed] [Google Scholar]
- Giri B, Smaldino PJ, Thys RG, Creacy SD, Routh ED, Hantgan RR, Lattmann S, Nagamine Y, Akman SA, Vaughn JP. G4 resolvase 1 tightly binds and unwinds unimolecular G4-DNA. Nucleic Acids Research. 2011;39:7161–7178. doi: 10.1093/nar/gkr234. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Greene AL, Snipe JR, Gordenin DA, Resnick MA. Functional analysis of human FEN1 in Saccharomyces cerevisiae and its role in genome stability. Human Molecular Genetics. 1999;8:2263–2273. doi: 10.1093/hmg/8.12.2263. [DOI] [PubMed] [Google Scholar]
- Holthöfer H, Kretzler M, Haltia A, Solin ML, Taanman JW, Schägger H, Kriz W, Kerjaschki D, Schlöndorff D. Altered gene expression and functions of mitochondria in human nephrotic syndrome. FASEB Journal. 1999;13:523–532. doi: 10.1096/fasebj.13.3.523. [DOI] [PubMed] [Google Scholar]
- Hudson EK, Hogue BA, Souza-Pinto NC, Croteau DL, Anson RM, Bohr VA, Hansford RG. Age-Associated change in mitochondrial DNA damage. Free Radical Research. 1998;29:573–579. doi: 10.1080/10715769800300611. [DOI] [PubMed] [Google Scholar]
- Javadekar SM, Nilavar NM, Paranjape A, Das K, Raghavan SC. Characterization of G-quadruplex antibody reveals differential specificity for G4 DNA forms. DNA Research. 2020;27:dsaa024. doi: 10.1093/dnares/dsaa024. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Jin Y, Yu Q, Zhou D, Chen L, Huang X, Xu G, Huang J, Gao X, Gao Y, Shen L. The mitochondrial DNA 9-bp deletion polymorphism is a risk factor for hepatocellular carcinoma in the Chinese population. Genetic Testing and Molecular Biomarkers. 2012;16:330–334. doi: 10.1089/gtmb.2011.0208. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Jose D, Datta K, Johnson NP, von Hippel PH. Spectroscopic studies of position-specific DNA “breathing” fluctuations at replication forks and primer-template junctions. PNAS. 2009;106:4231–4236. doi: 10.1073/pnas.0900803106. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Kalifa L, Beutner G, Phadnis N, Sheu SS, Sia EA. Evidence for a role of FEN1 in maintaining mitochondrial DNA integrity. DNA Repair. 2009;8:1242–1249. doi: 10.1016/j.dnarep.2009.07.008. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Klungland A, Lindahl T. Second pathway for completion of human DNA base excision-repair: reconstitution with purified proteins and requirement for DNase IV (FEN1) The EMBO Journal. 1997;16:3341–3348. doi: 10.1093/emboj/16.11.3341. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Komandur S, Venkatasubramanian S, Alluri RV, Rao P, Rao P, Hasan Q. Mitochondrial insertion-deletion polymorphism: role in disease pathology. Genetic Testing and Molecular Biomarkers. 2011;15:361–364. doi: 10.1089/gtmb.2010.0205. [DOI] [PubMed] [Google Scholar]
- Krishnan KJ, Turnbull DM. Mitochondrial DNA and genetic disease. Essays in Biochemistry. 2010;47:139–151. doi: 10.1042/bse0470139. [DOI] [PubMed] [Google Scholar]
- Kumar TS, Kari V, Choudhary B, Nambiar M, Akila TS, Raghavan SC. Anti-Apoptotic protein BCL2 down-regulates DNA end joining in cancer cells. The Journal of Biological Chemistry. 2010;285:32657–32670. doi: 10.1074/jbc.M110.140350. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Kumari R, Nambiar M, Shanbagh S, Raghavan SC. Detection of G-quadruplex DNA using primer extension as a tool. PLOS ONE. 2015;10:e0119722. doi: 10.1371/journal.pone.0119722. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Kumari N, Vartak SV, Dahal S, Kumari S, Desai SS, Gopalakrishnan V, Choudhary B, Raghavan SC. G-Quadruplex structures contribute to differential radiosensitivity of the human genome. IScience. 2019;21:288–307. doi: 10.1016/j.isci.2019.10.033. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Kumari R, Roy U, Desai S, Nilavar NM, Van Nieuwenhuijze A, Paranjape A, Radha G, Bawa P, Srivastava M, Nambiar M, Balaji KN, Liston A, Choudhary B, Raghavan SC. Microrna mir-29c regulates RAG1 expression and modulates V (D) J recombination during B cell development. Cell Reports. 2021;36:109390. doi: 10.1016/j.celrep.2021.109390. [DOI] [PubMed] [Google Scholar]
- Lee HC, Pang CY, Hsu HS, Wei YH. Differential accumulations of 4,977 bp deletion in mitochondrial DNA of various tissues in human ageing. Biochimica et Biophysica Acta. 1994;1226:37–43. doi: 10.1016/0925-4439(94)90056-6. [DOI] [PubMed] [Google Scholar]
- Li LY, Luo X, Wang X. Endonuclease G is an apoptotic DNase when released from mitochondria. Nature. 2001;412:95–99. doi: 10.1038/35083620. [DOI] [PubMed] [Google Scholar]
- Lin W, Sampathi S, Dai H, Liu C, Zhou M, Hu J, Huang Q, Campbell J, Shin-Ya K, Zheng L, Chai W, Shen B. Mammalian Dna2 helicase/nuclease cleaves G-quadruplex DNA and is required for telomere integrity. The EMBO Journal. 2013;32:1425–1439. doi: 10.1038/emboj.2013.88. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Liu P, Qian L, Sung JS, de Souza-Pinto NC, Zheng L, Bogenhagen DF, Bohr VA, Wilson DM, Shen B, Demple B. Removal of oxidative DNA damage via FEN1-dependent long-patch base excision repair in human cell mitochondria. Molecular and Cellular Biology. 2008;28:4975–4987. doi: 10.1128/MCB.00457-08. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Longo PA, Kavran JM, Kim MS, Leahy DJ. Transient mammalian cell transfection with polyethylenimine (PEI) Methods in Enzymology. 2013;529:227–240. doi: 10.1016/B978-0-12-418687-3.00018-5. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Marni G. Mitochondrial DNA Deletion Syndromes. Seattle): University of Washington; 2003. [PubMed] [Google Scholar]
- Masuda-Sasa T, Polaczek P, Peng XP, Chen L, Campbell JL. Processing of G4 DNA by Dna2 helicase/nuclease and replication protein A (RPA) provides insights into the mechanism of Dna2/RPA substrate recognition. The Journal of Biological Chemistry. 2008;283:24359–24373. doi: 10.1074/jbc.M802244200. [DOI] [PMC free article] [PubMed] [Google Scholar]
- McDermott-Roe C, Ye J, Ahmed R, Sun X-M, Serafín A, Ware J, Bottolo L, Muckett P, Cañas X, Zhang J, Rowe GC, Buchan R, Lu H, Braithwaite A, Mancini M, Hauton D, Martí R, García-Arumí E, Hubner N, Jacob H, Serikawa T, Zidek V, Papousek F, Kolar F, Cardona M, Ruiz-Meana M, García-Dorado D, Comella JX, Felkin LE, Barton PJR, Arany Z, Pravenec M, Petretto E, Sanchis D, Cook SA. Endonuclease G is a novel determinant of cardiac hypertrophy and mitochondrial function. Nature. 2011;478:114–118. doi: 10.1038/nature10490. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Michikawa Y, Mazzucchelli F, Bresolin N, Scarlato G, Attardi G. Aging-dependent large accumulation of point mutations in the human mtDNA control region for replication. Science. 1999;286:774–779. doi: 10.1126/science.286.5440.774. [DOI] [PubMed] [Google Scholar]
- Misic V, El-Mogy M, Haj-Ahmad Y. Role of endonuclease G in exogenous DNA stability in HeLa cells. Biochemistry. Biokhimiia. 2016;81:163–175. doi: 10.1134/S0006297916020103. [DOI] [PubMed] [Google Scholar]
- Modica-Napolitano JS, Kulawiec M, Singh KK. Mitochondria and human cancer. Current Molecular Medicine. 2007;7:121–131. doi: 10.2174/156652407779940495. [DOI] [PubMed] [Google Scholar]
- Naik AK, Lieber MR, Raghavan SC. Cytosines, but not purines, determine recombination activating gene (RAG) -induced breaks on heteroduplex DNA structures: implications for genomic instability. The Journal of Biological Chemistry. 2010;285:7587–7597. doi: 10.1074/jbc.M109.089631. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Nambiar M, Goldsmith G, Moorthy BT, Lieber MR, Joshi MV, Choudhary B, Hosur RV, Raghavan SC. Formation of a G-quadruplex at the BCL2 major breakpoint region of the T (14; 18) translocation in follicular lymphoma. Nucleic Acids Research. 2011;39:936–948. doi: 10.1093/nar/gkq824. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Nambiar M, Raghavan SC. Mechanism of fragility at BCL2 gene minor breakpoint cluster region during T (14; 18) chromosomal translocation. The Journal of Biological Chemistry. 2012;287:8688–8701. doi: 10.1074/jbc.M111.307363. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Nambiar M, Srivastava M, Gopalakrishnan V, Sankaran SK, Raghavan SC. G-Quadruplex structures formed at the HOX11 breakpoint region contribute to its fragility during T (10; 14) translocation in T-cell leukemia. Molecular and Cellular Biology. 2013;33:4266–4281. doi: 10.1128/MCB.00540-13. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Nishana M, Raghavan SC. A non-B DNA can replace heptamer of V (D) J recombination when present along with a nonamer: implications in chromosomal translocations and cancer. The Biochemical Journal. 2012;448:115–125. doi: 10.1042/BJ20121031. [DOI] [PubMed] [Google Scholar]
- Ohsato T, Ishihara N, Muta T, Umeda S, Ikeda S, Mihara K, Hamasaki N, Kang D. Mammalian mitochondrial endonuclease G. digestion of R-loops and localization in intermembrane space. European Journal of Biochemistry. 2002;269:5765–5770. doi: 10.1046/j.1432-1033.2002.03238.x. [DOI] [PubMed] [Google Scholar]
- Oliveira PH, da Silva CL, Cabral JMS. An appraisal of human mitochondrial DNA instability: new insights into the role of non-canonical DNA structures and sequence motifs. PLOS ONE. 2013;8:e59907. doi: 10.1371/journal.pone.0059907. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Paeschke K, Bochman ML, Garcia PD, Cejka P, Friedman KL, Kowalczykowski SC, Zakian VA. Pif1 family helicases suppress genome instability at G-quadruplex motifs. Nature. 2013;497:458–462. doi: 10.1038/nature12149. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Pakendorf B, Stoneking M. Mitochondrial DNA and human evolution. Annual Review of Genomics and Human Genetics. 2005;6:165–183. doi: 10.1146/annurev.genom.6.080604.162249. [DOI] [PubMed] [Google Scholar]
- Parrish J, Li L, Klotz K, Ledwich D, Wang X, Xue D. Mitochondrial endonuclease G is important for apoptosis in C. elegans. Nature. 2001;412:90–94. doi: 10.1038/35083608. [DOI] [PubMed] [Google Scholar]
- Penta JS, Johnson FM, Wachsman JT, Copeland WC. Mitochondrial DNA in human malignancy. Mutation Research. 2001;488:119–133. doi: 10.1016/s1383-5742(01)00053-9. [DOI] [PubMed] [Google Scholar]
- Phillips AF, Millet AR, Tigano M, Dubois SM, Crimmins H, Babin L, Charpentier M, Piganeau M, Brunet E, Sfeir A. Single-Molecule analysis of mtDNA replication uncovers the basis of the common deletion. Molecular Cell. 2017;65:527–538. doi: 10.1016/j.molcel.2016.12.014. [DOI] [PubMed] [Google Scholar]
- Puig Lombardi E, Londoño-Vallejo A. A guide to computational methods for G-quadruplex prediction. Nucleic Acids Research. 2020;48:1–15. doi: 10.1093/nar/gkz1097. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Raghavan SC, Swanson PC, Wu X, Hsieh CL, Lieber MR. A non-B-DNA structure at the Bcl-2 major breakpoint region is cleaved by the RAG complex. Nature. 2004;428:88–93. doi: 10.1038/nature02355. [DOI] [PubMed] [Google Scholar]
- Raghavan SC, Chastain P, Lee JS, Hegde BG, Houston S, Langen R, Hsieh CL, Haworth IS, Lieber MR. Evidence for a triplex DNA conformation at the bcl-2 major breakpoint region of the T (14; 18) translocation. Journal of Biological Chemistry. 2005a;280:22749–22760. doi: 10.1074/jbc.M502952200. [DOI] [PubMed] [Google Scholar]
- Raghavan SC, Hsieh CL, Lieber MR. Both V (D) J coding ends but neither signal end can recombine at the bcl-2 major breakpoint region, and the rejoining is ligase IV dependent. Molecular and Cellular Biology. 2005b;25:6475–6484. doi: 10.1128/MCB.25.15.6475-6484.2005. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Raghavan SC, Swanson PC, Ma Y, Lieber MR. Double-strand break formation by the RAG complex at the bcl-2 major breakpoint region and at other non-B DNA structures in vitro. Molecular and Cellular Biology. 2005c;25:5904–5919. doi: 10.1128/MCB.25.14.5904-5919.2005. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Raghavan SC, Tsai A, Hsieh CL, Lieber MR. Analysis of non-B DNA structure at chromosomal sites in the mammalian genome. Methods in Enzymology. 2006;409:301–316. doi: 10.1016/S0076-6879(05)09017-8. [DOI] [PubMed] [Google Scholar]
- Ray U, Raul SK, Gopinatha VK, Ghosh D, Rangappa KS, Mantelingu K, Raghavan SC. Identification and characterization of novel SCR7-based small-molecule inhibitor of DNA end-joining, SCR130 and its relevance in cancer therapeutics. Molecular Carcinogenesis. 2020;59:618–628. doi: 10.1002/mc.23186. [DOI] [PubMed] [Google Scholar]
- Raymond C, Tom R, Perret S, Moussouami P, L’Abbé D, St-Laurent G, Durocher Y. A simplified polyethylenimine-mediated transfection process for large-scale and high-throughput applications. Methods. 2011;55:44–51. doi: 10.1016/j.ymeth.2011.04.002. [DOI] [PubMed] [Google Scholar]
- Redd AJ, Takezaki N, Sherry ST, McGarvey ST, Sofro AS, Stoneking M. Evolutionary history of the COII/trnalys intergenic 9 base pair deletion in human mitochondrial DNAs from the Pacific. Molecular Biology and Evolution. 1995;12:604–615. doi: 10.1093/oxfordjournals.molbev.a040240. [DOI] [PubMed] [Google Scholar]
- Richards M, Macaulay V. The mitochondrial gene tree comes of age. American Journal of Human Genetics. 2001;68:1315–1320. doi: 10.1086/320615. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Ruiz-Carrillo A, Renaud J. Endonuclease G: a (dG) N X (dC) n-specific DNase from higher eukaryotes. The EMBO Journal. 1987;6:401–407. doi: 10.1002/j.1460-2075.1987.tb04769.x. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Sambrook J, Fritsch EF, Maniatis T. Molecular Cloning A Laboratory Manual. Second edn. New York: Cold Spring Harbor Laboratory Press; 1989. [Google Scholar]
- Samuels DC, Schon EA, Chinnery PF. Two direct repeats cause most human mtDNA deletions. Trends in Genetics. 2004;20:393–398. doi: 10.1016/j.tig.2004.07.003. [DOI] [PubMed] [Google Scholar]
- Sanchez M, Mercer TR, Davies SMK, Shearwood AMJ, Nygård KKA, Richman TR, Mattick JS, Rackham O, Filipovska A. Rna processing in human mitochondria. Cell Cycle. 2011;10:2904–2916. doi: 10.4161/cc.10.17.17060. [DOI] [PubMed] [Google Scholar]
- Schäfer P, Scholz SR, Gimadutdinow O, Cymerman IA, Bujnicki JM, Ruiz-Carrillo A, Pingoud A, Meiss G. Structural and functional characterization of mitochondrial endog, a sugar non-specific nuclease which plays an important role during apoptosis. Journal of Molecular Biology. 2004;338:217–228. doi: 10.1016/j.jmb.2004.02.069. [DOI] [PubMed] [Google Scholar]
- Sharma S, Choudhary B, Raghavan SC. Efficiency of nonhomologous DNA end joining varies among somatic tissues, despite similarity in mechanism. Cellular and Molecular Life Sciences. 2011;68:661–676. doi: 10.1007/s00018-010-0472-x. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Sharma S, Javadekar SM, Pandey M, Srivastava M, Kumari R, Raghavan SC. Homology and enzymatic requirements of microhomology-dependent alternative end joining. Cell Death & Disease. 2015;6:e1697. doi: 10.1038/cddis.2015.58. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Sinha D, Joshi N, Chittoor B, Samji P, D’Silva P. Role of Magmas in protein transport and human mitochondria biogenesis. Human Molecular Genetics. 2010;19:1248–1262. doi: 10.1093/hmg/ddq002. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Soodyall H, Vigilant L, Hill AV, Stoneking M, Jenkins T. Mtdna control-region sequence variation suggests multiple independent origins of an “ asian-specific ” 9-bp deletion in sub-Saharan Africans. American Journal of Human Genetics. 1996;58:595–608. [PMC free article] [PubMed] [Google Scholar]
- Srivastava S, Moraes CT. Double-Strand breaks of mouse muscle mtDNA promote large deletions similar to multiple mtDNA deletions in humans. Human Molecular Genetics. 2005;14:893–902. doi: 10.1093/hmg/ddi082. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Srivastava M, Nambiar M, Sharma S, Karki SS, Goldsmith G, Hegde M, Kumar S, Pandey M, Singh RK, Ray P, Natarajan R, Kelkar M, De A, Choudhary B, Raghavan SC. An inhibitor of nonhomologous end-joining abrogates double-strand break repair and impedes cancer progression. Cell. 2012;151:1474–1487. doi: 10.1016/j.cell.2012.11.054. [DOI] [PubMed] [Google Scholar]
- Sun H, Yabuki A, Maizels N. A human nuclease specific for G4 DNA. PNAS. 2001;98:12444–12449. doi: 10.1073/pnas.231479198. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Suzuki H, Maeda M, Mihara K. Characterization of rat Tom70 as a receptor of the preprotein translocase of the mitochondrial outer membrane. Journal of Cell Science. 2002;115:1895–1905. doi: 10.1242/jcs.115.9.1895. [DOI] [PubMed] [Google Scholar]
- Tadi SK, Sebastian R, Dahal S, Babu RK, Choudhary B, Raghavan SC. Microhomology-Mediated end joining is the principal mediator of double-strand break repair during mitochondrial DNA lesions. Molecular Biology of the Cell. 2016;27:223–235. doi: 10.1091/mbc.E15-05-0260. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Tan AS, Baty JW, Dong L-F, Bezawork-Geleta A, Endaya B, Goodwin J, Bajzikova M, Kovarova J, Peterka M, Yan B, Pesdar EA, Sobol M, Filimonenko A, Stuart S, Vondrusova M, Kluckova K, Sachaphibulkij K, Rohlena J, Hozak P, Truksa J, Eccles D, Haupt LM, Griffiths LR, Neuzil J, Berridge MV. Mitochondrial genome acquisition restores respiratory function and tumorigenic potential of cancer cells without mitochondrial DNA. Cell Metabolism. 2015;21:81–94. doi: 10.1016/j.cmet.2014.12.003. [DOI] [PubMed] [Google Scholar]
- Taylor RW, Turnbull DM. Mitochondrial DNA mutations in human disease. Nature Reviews. Genetics. 2005;6:389–402. doi: 10.1038/nrg1606. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Thomas MG, Cook CE, Miller KW, Waring MJ, Hagelberg E. Molecular instability in the COII-trna (Lys) intergenic region of the human mitochondrial genome: multiple origins of the 9-bp deletion and heteroplasmy for expanded repeats. Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences. 1998;353:955–965. doi: 10.1098/rstb.1998.0260. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Thomas E, Gopalakrishnan V, Hegde M, Kumar S, Karki SS, Raghavan SC, Choudhary B. A novel resveratrol based tubulin inhibitor induces mitotic arrest and activates apoptosis in cancer cells. Scientific Reports. 2016;6:34653. doi: 10.1038/srep34653. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Totaro A, Renzi F, La Fata G, Mattioli C, Raabe M, Urlaub H, Achsel T. The human Pat1b protein: a novel mRNA deadenylation factor identified by a new immunoprecipitation technique. Nucleic Acids Research. 2011;39:635–647. doi: 10.1093/nar/gkq797. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Van Goethem G, Martin JJ, Van Broeckhoven C. Progressive external ophthalmoplegia characterized by multiple deletions of mitochondrial DNA: unraveling the pathogenesis of human mitochondrial DNA instability and the initiation of a genetic classification. Neuromolecular Medicine. 2003;3:129–146. doi: 10.1385/NMM:3:3:129. [DOI] [PubMed] [Google Scholar]
- van Loo G, Schotte P, van Gurp M, Demol H, Hoorelbeke B, Gevaert K, Rodriguez I, Ruiz-Carrillo A, Vandekerckhove J, Declercq W, Beyaert R, Vandenabeele P. Endonuclease G: a mitochondrial protein released in apoptosis and involved in caspase-independent DNA degradation. Cell Death and Differentiation. 2001;8:1136–1142. doi: 10.1038/sj.cdd.4400944. [DOI] [PubMed] [Google Scholar]
- Voineagu I, Narayanan V, Lobachev KS, Mirkin SM. Replication stalling at unstable inverted repeats: interplay between DNA hairpins and fork stabilizing proteins. PNAS. 2008;105:9936–9941. doi: 10.1073/pnas.0804510105. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Wells RD. Discovery of the role of non-B DNA structures in mutagenesis and human genomic disorders. The Journal of Biological Chemistry. 2009;284:8997–9009. doi: 10.1074/jbc.X800010200. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Wiehe RS, Gole B, Chatre L, Walther P, Calzia E, Ricchetti M, Wiesmüller L. Endonuclease G promotes mitochondrial genome cleavage and replication. Oncotarget. 2018;9:18309–18326. doi: 10.18632/oncotarget.24822. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Wisnovsky S, Sack T, Pagliarini DJ, Laposa RR, Kelley SO. DNA polymerase θ increases mutational rates in mitochondrial DNA. ACS Chemical Biology. 2018;13:900–908. doi: 10.1021/acschembio.8b00072. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Yakes FM, Van Houten B. Mitochondrial DNA damage is more extensive and persists longer than nuclear DNA damage in human cells following oxidative stress. PNAS. 1997;94:514–519. doi: 10.1073/pnas.94.2.514. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Yang JH, Lee HC, Lin KJ, Wei YH. A specific 4977-bp deletion of mitochondrial DNA in human ageing skin. Archives of Dermatological Research. 1994;286:386–390. doi: 10.1007/BF00371798. [DOI] [PubMed] [Google Scholar]
- Yang C, Wu R, Liu H, Chen Y, Gao Y, Chen X, Li Y, Ma J, Li J, Gan J. Structural insights into DNA degradation by human mitochondrial nuclease MGME1. Nucleic Acids Research. 2018;46:11075–11088. doi: 10.1093/nar/gky855. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Yao YG, Watkins WS, Zhang YP. Evolutionary history of the mtDNA 9-bp deletion in Chinese populations and its relevance to the peopling of East and Southeast Asia. Human Genetics. 2000;107:504–512. doi: 10.1007/s004390000403. [DOI] [PubMed] [Google Scholar]
- Yen TC, Su JH, King KL, Wei YH. Ageing-Associated 5 kb deletion in human liver mitochondrial DNA. Biochemical and Biophysical Research Communications. 1991;178:124–131. doi: 10.1016/0006-291x(91)91788-e. [DOI] [PubMed] [Google Scholar]
- Yeretssian G, Correa RG, Doiron K, Fitzgerald P, Dillon CP, Green DR, Reed JC, Saleh M. Non-Apoptotic role of Bid in inflammation and innate immunity. Nature. 2011;474:96–99. doi: 10.1038/nature09982. [DOI] [PubMed] [Google Scholar]
- Yu X, Baer R. Nuclear localization and cell cycle-specific expression of CtIP, a protein that associates with the BRCA1 tumor suppressor. The Journal of Biological Chemistry. 2000;275:18541–18549. doi: 10.1074/jbc.M909494199. [DOI] [PubMed] [Google Scholar]
- Yusoff AAM, Abdullah WSW, Khair SZNM, Radzak SMA. A comprehensive overview of mitochondrial DNA 4977-bp deletion in cancer studies. Oncology Reviews. 2019;13:409. doi: 10.4081/oncol.2019.409. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Zeviani M, Moraes CT, DiMauro S, Nakase H, Bonilla E, Schon EA, Rowland LP. Deletions of mitochondrial DNA in Kearns-Sayre syndrome. Neurology. 1988;38:1339–1346. doi: 10.1212/wnl.38.9.1339. [DOI] [PubMed] [Google Scholar]
- Zhou Q, Li H, Li H, Nakagawa A, Lin JLJ, Lee ES, Harry BL, Skeen-Gaar RR, Suehiro Y, William D, Mitani S, Yuan HS, Kang BH, Xue D. Mitochondrial endonuclease G mediates breakdown of paternal mitochondria upon fertilization. Science. 2016;353:394–399. doi: 10.1126/science.aaf4777. [DOI] [PMC free article] [PubMed] [Google Scholar]
- Zhuo G, Feng G, Leng J, Yu L, Jiang Y. A 9-bp deletion homoplasmy in women with polycystic ovary syndrome revealed by mitochondrial genome-mutation screen. Biochemical Genetics. 2010;48:157–163. doi: 10.1007/s10528-009-9308-5. [DOI] [PubMed] [Google Scholar]