Skip to main content
eLife logoLink to eLife
. 2022 Nov 17;11:e69916. doi: 10.7554/eLife.69916

Unleashing a novel function of Endonuclease G in mitochondrial genome instability

Sumedha Dahal 1, Humaira Siddiqua 1, Shivangi Sharma 1, Ravi K Babu 1, Diksha Rathore 1, Sheetal Sharma 2, Sathees C Raghavan 1,
Editors: Wolf-Dietrich Heyer3, Jessica K Tyler4
PMCID: PMC9711528  PMID: 36394256

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.

Figure 1. Biochemical studies to investigate formation of G-quadruplex structure in Region I of the mitochondrial genome.

(A) Sequence of the mitochondrial Region I with direct repeats and inverted repeats. Blue-shaded arrows represent the direct repeats, and yellow arrows represent the inverted repeats. The sequence in blue is G rich sequence predicted to form a G-quadruplex structure. The sequence after 9 bp deletion is also shown. (B) The oligomeric sequence of the G-rich strand (indicated by ‘G1’) and the complementary C-rich strand (indicated by ‘C1’). These oligomers are derived from the predicted mitochondrial G-quadruplex forming Region I. The stretches of guanines are marked in blue. (C) The G and C-rich strands were incubated in the presence of 100 mM KCl and resolved on 15% native polyacrylamide gels in the absence (-KCl) or presence (+KCl) of KCl (100 mM) in the gel and running buffer. The substrate, intramolecular (Intra G), and intermolecular (Inter G) quadruplex structures are indicated. (D) CD spectra for G1 in the presence (red) or absence of 100 mM KCl (black). CD spectra for G1 (green) and C1 (black) in the presence of 100 mM KCl. In each case where KCl was added, the respective oligomer resuspended in Tris-EDTA buffer was incubated for 1 h at 37 °C and spectra were recorded using JASCO J-810 spectropolarimeter (scan range of 220–300 nm). (E) DMS protection assay for the Region I of mitochondria. The oligomer was allowed to fold into G4 DNA and then treated with DMS, followed by cleavage with piperidine. The products were resolved on a 15% denaturing PAGE. Substrate indicates the gel-purified DNA from the reaction incubated with or without KCl. All the positions of guanines are indicated. (F) A representative 2D model for intramolecular G-quadruplex structure formed at mitochondrial Region I based on the reactivity of guanine to DMS and CD studies for determining the orientation of DNA fold. Refer also Figure 1—figure supplement 1.

Figure 1—source data 1. EMSA and CD studies to show formation of G-quadruplex structure in Region I.
Figure 1—source data 2. Gel shift assay to show abrogation of G-quadruplex structure in mutants of Region I.

Figure 1.

Figure 1—figure supplement 1. Evaluation of G-quadruplex formation at mitochondrial Region I.

Figure 1—figure supplement 1.

Related to Figure 1. (A) Schematic representation of mitochondrial genome showing positions of five G-quadruplex DNA motifs. (B) List of oligomers used for studying effect of mutation of G stretches on G-quadruplex formation. Guanine stretches are indicated in blue, while sequences in red are mutations introduced to G stretches. (C) Gel shift assay in the absence (left panel) and presence (right panel) of KCl showing the impact of mutation of G stretches on intramolecular G4 DNA formation.

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—source data 1. Primer extension studies to show formation of G-quadruplex structure in Region I containing plasmid.

Figure 2.

Figure 2—figure supplement 1. Schematic showing cloning of the mitochondrial Region I and its mutant to generate plasmids, pDI1 and pDI2.

Figure 2—figure supplement 1.

Related to Figure 2. The duplex region containing the G stretches are depicted in blue, while the mutated nucleotides are marked in red. A 68 nucleotide sequence is also shown with inverted repeats marked in orange.

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—source data 1. Sequence of clones after bisulphite treatment in a plasmid containing Region I of mitochondria.
Figure 3—source data 2. Sequence of clones after bisulphite treatment in a region I of mitochondria.

Figure 3.

Figure 3—figure supplement 1. Schematic showing conversion of cytosine to thymine following bisulfite modification assay.

Figure 3—figure supplement 1.

Related to Figure 3. (A, B) Treatment with sodium bisulfite can result in deamination of cytosine, leading to uracil when present on a single-stranded DNA (A). This C to U conversion can lead to C to T change following PCR and sequencing (B). (C) Table showing the summary of clone sequenced from mitochondrial region containing plasmid and mitochondrial genome for bisulfite sequencing for both top and bottom strand.

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—source data 1. Localization of BG4 to mitochondria.
Figure 4—source data 2. Immunofluorescence showing the localization of BG4 to mitochondria.

Figure 4.

Figure 4—figure supplement 1. Evaluation of existence of G-quadruplex in mitochondria of cells.

Figure 4—figure supplement 1.

Related to Figure 4. (A) Immunofluorescence study showing colocalization of BG4 (antibody that binds to G4 DNA) and mitochondrial genome. Nucleus from HeLa cells is stained with DAPI (blue color). Mitochondria are stained with MitoTracker green FM and G4 DNA with BG4 (Alexa-Fluor 568), a merged image is shown. A merge of red and green is depicted by Coste’s mask (colocalization is represented as a white dot). The experiments were done independently in HeLa and HEK 293T cells and the quantitation of colocalization of BG4 with MitoTracker is shown in Figure 2. (B) Increasing concentrations (0.2, 0.5, 1, and 2 µg) of purified BG4 was incubated with G-rich (G1) and C-rich (C1) oligomers and products were resolved on 5% native polyacrylamide gel. In each case a lane without protein (lane 1 and 6) served as control. The substrate and the bound complex are marked.

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—source data 1. BG4 ChIP to show the binding of BG4 to mitochondrial G-quadruplex forming regions.

Figure 5.

Figure 5—figure supplement 1. Evaluation of existence of G-quadruplex in mitochondrial DNA.

Figure 5—figure supplement 1.

Related to Figure 5. Schematic showing the position of primers used for mito IP studies. GR1-GR5 represents primers that can amplify G-quadruplex forming motifs and CR1-CR10 represents the random control regions.

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—source data 1. Gel profiles showing the mitochondrial induced cleavage at mitochondrial region I.
Figure 6—source data 2. Cleavage assay on plasmid bearing wildtype and mutant G4 sequence.

Figure 6.

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 1.

Related to Figure 6. (A) Table showing the positions of different mutations at G4 DNA motif on plasmid bearing Region I. (B) Cleavage assay using different mutations at G4 DNA motif on plasmid bearing Region I.
Figure 6—figure supplement 2. DNA sequences that support formation of different conformation of G-quadruplex DNA structures in Region I.

Figure 6—figure supplement 2.

Related to Figure 6. The presence of five G-stretches (shown in blue) and three GNG stretches (shown in green) allows this region to fold into four different conformations of G-quadruplex structure.

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—source data 1. Gel profiles showing the Endonuclease-G-induced cleavage at mitochondrial region I.
Figure 7—source data 2. Gel profiles showing the purification of different endonucleases.

Figure 7.

Figure 7—figure supplement 1. Overexpression, purification and activity assay of different endonucleases.

Figure 7—figure supplement 1.

Related to Figure 7. (A, B) Overexpression and purification of Endonuclease G protein (wild type and mutant). The purity and identity of the protein was confirmed using SDS-PAGE for wild type (A) and mutant (B). (C) The identity and purity of the protein was confirmed using the western blotting for both wild type and mutant. Fractions 5 and 6 (wild type) and 4 and 5 (ΔH151A) were used for the western blot in both the cases. (D). SDS gel profile showing the overexpression and purification of CtIP. Western blotting (bottom panel) was performed for confirmation of the protein. (E) Activity assay showing resection of the substrate (ssDNA) upon addition of increasing concentrations (30, 60 and 90 ng) of purified CtIP (lanes 2–4). (F) SDS gel profile showing the overexpression and purification of FEN1. Western blotting confirming the identity of the protein is also shown. (G) Activity assay showing the cleavage of the substrate with flap DNA upon the addition of increasing concentrations (50, 100, and 500 ng) of purified FEN1 (lanes 2–5). (H) Silver-stained gel profile showing the overexpression and purification of RAGs from mammalian cells. Western blotting confirming the presence of RAG1 and RAG2. (I) Activity assay showing the cleavage of the RSS substrate upon addition of different fractions of purified RAGs (lanes 2–8). (J) Cleavage assay for different nucleases was performed on a plasmid containing mitochondrial Region I, pDI1. Following incubation with different purified nucleases CtIP, FEN1, and RAGs, primer extension assay was carried out using γ-32P radiolabeled VKK11 primer and resolved on 8% denaturing polyacrylamide gel along with sequencing ladder. Lane 1 is primer alone, lanes 2, 5, 8 are without protein, lane 11 is with mitochondrial extract and lanes 12–15 (A, C, G, T) represents the sequencing ladder. Lanes 3,4 represents (CtIP); lanes 6,7 (FEN1); lanes 9, 10 (RAGs) incubated samples. In each case, 30 and 60 ng of protein was used. ‘M’ is 50 bp ladder.

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—source data 1. Localization of Endonuclease G to mitochondria.
Figure 8—source data 2. Immunofluorescence showing the Localization of Endonuclease G to mitochondria.

Figure 8.

Figure 8—figure supplement 1. Analysis of Endonuclease G localization in mitochondria.

Figure 8—figure supplement 1.

Related to Figure 8.
Representative image showing colocalization of Endonuclease G using MitoTracker Green FM (Mt Green) in HeLa cells. Alexa-568 conjugated secondary antibody was used for the detection of Endonuclease G. DAPI is used to stain the nucleus.

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—source data 1. ChIP assay showing the binding of Endonuclease G with the mitochondrial G-quadruplex regions within cells.
Figure 9—source data 2. ChIP assay showing the binding of Endonuclease G with the mitochondrial G-quadruplex regions when purified Endonuclease G was used.
Figure 9—source data 3. P1 nuclease assay showing the binding of Endonuclease G to mitochondrial G quadruples regions.

Figure 9.

Figure 9—figure supplement 1. Binding of Endonuclease G to G-quadruplex regions of the mitochondrial genome.

Figure 9—figure supplement 1.

Related to Figure 9. (A) Representative immunofluorescence images showing the colocalization of Endonuclease G and BG4. The ‘Merged’ image shown in left is a colocalization of DAPI, Endonuclease G and BG4 foci, while ‘Merged’ in right shows colocalization of Endonuclease and BG4 foci. (B) Schematic showing the binding of purified Endonuclease G to the mitochondrial genome and the pull-down of a bound region using Endonuclease G antibody and protein A/G beads. These regions were then used for semi-quantitative and real-time PCR using appropriate primers. (C) Schematic showing the position of primers used for mito-IP studies. The upper panel shows the primer positions for G-quadruplex forming regions and lower panels for control regions. (D) 5 G-quadruplex forming regions and 10 random regions were used for amplification using mito-IP DNA as template (lower panel). Upper panel shows the amplification for input DNA. (E) Bar graph was plotted for the threshold cycle against different primers following real-time PCR. Input DNA served as template control. ‘Antibody’ control in which no antibody was added served as negative control. Blue bars represent G-quadruplex regions, while green bars represent the control regions. Error bar represents three independent biological repeats.
Figure 9—figure supplement 2. Binding of Endonuclease G to G-quadruplex regions of the mitochondrial genome.

Figure 9—figure supplement 2.

Related to Figure 9. (A) SDS profile and western blotting of the pulldown sample when mitochondrial DNA was incubated with mitochondrial extracts. As described in the methodology, mitochondrial extract was allowed to bind to mitochondrial DNA and Endonuclease G bound DNA was pulled down using Endonuclease G antibody. After the pull down, the extract was loaded into a SDS PAGE with Endonuclease G (lane 3) or IgG antibody (lane 2). For reference, mitochondrial extract was also loaded (lane 4). M is the protein marker. For other details refer main Figure 9. (B) P1 nuclease foot-printing to investigate the binding of Endonuclease G to G-quadruplex forming Region I. Radiolabelled oligomers were incubated with Endonuclease G protein and subjected to P1 nuclease and electrophoresed on 18% denaturing PAGE. In each case, lanes 1, 5, 9, 13, 17, 21 are substrate alone, lanes 2, 6, 10, 14, 18, 22 are Endonuclease G alone. Lanes 3, 7, 11, 15, 19, 23 are P1 nuclease alone treated samples and lanes 4, 8, 12, 16, 20, 24 are Endonuclease G plus P1 nuclease treated samples. C1 is C-strand, G1 is G strand, M1, M2 and M3 are mutants while random sequence (RN) is the oligomer with equal G-C content as G1 in a random manner. 50 ng of purified Endonuclease G and 0.03 U of P1 nuclease was used for the assay.

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.

Figure 10.

(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.

Figure 10—source data 1. Sub localization of Endonuclease G with or without induction of stress.

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.

Figure 11.

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

Transparent reporting form

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

  1. 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]
  2. 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]
  3. 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]
  4. 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]
  5. 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]
  6. 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]
  7. 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]
  8. 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]
  9. 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]
  10. 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]
  11. 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]
  12. 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]
  13. 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]
  14. 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]
  15. 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]
  16. 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]
  17. 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]
  18. 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]
  19. 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]
  20. 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]
  21. 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]
  22. 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]
  23. 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]
  24. 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]
  25. 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]
  26. 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]
  27. 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]
  28. 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]
  29. 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]
  30. García-Lepe UO, Bermúdez-Cruz RM. Mitochondrial Genome Maintenance: Damage and Repair Pathways. UK: IntechOpen; 2019. [Google Scholar]
  31. 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]
  32. 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]
  33. 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]
  34. 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]
  35. 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]
  36. 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]
  37. 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]
  38. 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]
  39. 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]
  40. 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]
  41. 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]
  42. Krishnan KJ, Turnbull DM. Mitochondrial DNA and genetic disease. Essays in Biochemistry. 2010;47:139–151. doi: 10.1042/bse0470139. [DOI] [PubMed] [Google Scholar]
  43. 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]
  44. 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]
  45. 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]
  46. 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]
  47. 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]
  48. 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]
  49. 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]
  50. 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]
  51. 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]
  52. Marni G. Mitochondrial DNA Deletion Syndromes. Seattle): University of Washington; 2003. [PubMed] [Google Scholar]
  53. 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]
  54. 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]
  55. 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]
  56. 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]
  57. 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]
  58. 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]
  59. 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]
  60. 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]
  61. 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]
  62. 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]
  63. 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]
  64. 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]
  65. 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]
  66. 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]
  67. 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]
  68. 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]
  69. 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]
  70. 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]
  71. 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]
  72. 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]
  73. 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]
  74. 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]
  75. 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]
  76. 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]
  77. 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]
  78. 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]
  79. 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]
  80. 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]
  81. Sambrook J, Fritsch EF, Maniatis T. Molecular Cloning A Laboratory Manual. Second edn. New York: Cold Spring Harbor Laboratory Press; 1989. [Google Scholar]
  82. 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]
  83. 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]
  84. 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]
  85. 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]
  86. 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]
  87. 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]
  88. 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]
  89. 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]
  90. 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]
  91. 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]
  92. 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]
  93. 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]
  94. 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]
  95. 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]
  96. 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]
  97. 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]
  98. 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]
  99. 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]
  100. 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]
  101. 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]
  102. 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]
  103. 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]
  104. 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]
  105. 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]
  106. 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]
  107. 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]
  108. 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]
  109. 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]
  110. 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]
  111. 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]
  112. 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]
  113. 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]
  114. 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]
  115. 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]

Editor's evaluation

Wolf-Dietrich Heyer 1

This manuscript is of interest to researchers in the field of mitochondrial genome stability and mitochondrial genetic diseases and reports valuable findings that convincingly demonstrate that endonuclease G preferentially binds to mitochondrial genome regions which have a potential for forming G4 tetraplexes and induces DNA breaks that may lead to a common 9 bp deletion in the mitochondrial genome.

Decision letter

Editor: Wolf-Dietrich Heyer1
Reviewed by: Wolf-Dietrich Heyer2, Sergei Mirkin

Our editorial process produces two outputs: (i) public reviews designed to be posted alongside the preprint for the benefit of readers; (ii) feedback on the manuscript for the authors, including requests for revisions, shown below. We also include an acceptance summary that explains what the editors found interesting or important about the work.

Decision letter after peer review:

Thank you for submitting your article "Unleashing a Novel Function of Endonuclease G in Mitochondrial Genome Instability" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, including Wolf-Dietrich Heyer as the Reviewing Editor and Reviewer #1, and the evaluation has been overseen by Jessica Tyler as the Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Sergei Mirkin (Reviewer #3).

The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.

Essential revisions:

The manuscript makes interesting observations that should advance the field of mitochondrial DNA metabolism and will be off interest for researchers in the field of mitochondrial genome stability and mitochondrial diseases.

The revision that is needed to make this manuscript, as submitted, publishable is too extensive. If the authors decide to address the comments by the additional experimentation needed, they are invited to resubmit a revised manuscript as a new submission that I would be happy to handle as BRE and seek to have the same reviewers evaluate the new submission.

Alternatively, the authors may want to consider tightening the manuscript and omit the reconstitution studies making sure that their interpretations are supported by the data and consider alternative mechanisms, when necessary. This more focused version would still contain sufficient new insights and allow the authors to further develop the reconstitution experiments with the required exhaustive separate incubation controls for a future publication.

1) Independent of specific criticism towards individual experiments or experimental tools, the manuscript is difficult to read, poorly assembled and would need a major overhaul in the text and figures. The text requires careful editing, there are many disagreements, and many cases where referrals such as 'its' or 'this' are unclear. The figures should be labeled that the experiments can be understood independent of the legend or method description. A few examples are given here.

• Figure 1A: The sequences of C1 and G1 are needed here. Giving more is confusing, as they do not relate to the figure. In Figure S1, all oligo sequences should be listed.

• Figure S1A: Define -/+ KCl under gel in legend.

• Figure S1B: Indicate inverted repeat sequences.

• Figure S3C: Where are the gel data? What is shown here? The description in the legend does not allow to understand the figure, what is the difference between the 3 lanes that are labeled as +.

• Figure S3: Unclear what the Westernblots show in D, F, and H.

• Figure S5F: Define the lanes. What is loaded in lane 4?

• Figure 4D: Label upper input and lower pull-down, also Figure 2G.

2) A key problem in the manuscript is that the authors mix their interpretation into the description of the results. A few examples are given here, but the authors should check the entire manuscript and correct this issue. The authors should carefully separate description of the results from their interpretation and refrain from overinterpreting their data (see also below). There is no need to overstate the result; exceptions or unexpected results must be discussed and acknowledged. For example:

• Figure 2F/G: What happened with CR7? Low threshold value but no PCR product? This is not described in the text.

• Figure 3D: The upper cleavage site stays the same, unlike what is described in the text.

• Figure S3A, B: It is unclear which band is EndoG. Is the Westernblot aligned with material loaded in slots 4-7 of the protein-stained gel? If so, why does the Westernblot intensity not reflect the staining intensity?

• Figure 3F: The Mt extract does only show the major EndoG sites but not the minor sites, which add up to a significant proportion of the products.

• Page 10, line 10: The mutant plasmid still has one cleavage site, cleavage is not abrogated as stated in the text.

• Page 10, line 22: Figure 3D is extract not purified EndoG as stated in line 22 of page 10.

• Figure 4E and F: What is the difference between both experiments? E: extract, F: cells? This could be clarified in the figure and text.

• Figure S5G: The description of the P1 results is not precise. There is clear protection of the C strand by EndoG in the upper third of the ladder.

• Figure 7A requires more description and quantitation. Figure 7 B/E: Why is the DAPI signal so variable? Figure 7G: the difference between lanes 1-3 and 4-6 should become apparent from the figure. For Figure 7G and H, only menadione data are shown, there are no IR data as stated in the text.

3) The BG4 antibody is problematic. Page 8, line 22: Why is there a BG4 signal in the cytoplasm in rho (0) cells? Is this also seen in rho (+) cells?

Also, the authors may want to discuss in more details image modification via Coste's mask. I can clearly see yellow dots in merged images for rho (0) cells, but somehow they are not converting to the white dotes when the Coste's mask is applied.

Figure 2 and 4: There are many BG4 signals? Do the signals all represent G4 tetraplex DNA? How can the authors ascertain the specificity and usefulness of this reagent?

The control in Figure 4A/B needs to be defined. What would be expectation from random overlap in this experiment given that both signals are very abundant?

4) The microscopy approach is less informative and overinterpreted. The level of resolution obtained with the confocal microscope approach used does not allow a conclusion on the localization beyond saying that the protein is present in the mitochondria and is not suitable for distinguishing the different mitochondrial compartments. As for the fractionation procedure used, I am not sure that it allows the separation of the matrix fraction from the intermembrane compartment. Furthermore, the description of the results displayed is not clear and there seem to be some mistakes in the labelling of the figures, making it harder to draw clear conclusions.

5) The mutant substrates are an excellent control for testing the G4 tetraplex hypothesis. What is the experimental evidence that these mutant oligonucleotides do not form tetraplexes?

6) The attempt to reconstruct the deletion formation in vitro is welcomed. However, these experiments lack a crucial control. In Figure 5B, there is product in the negative control? This raises the concern that bands labeled MMEJ products are PCR artifacts. The authors must conduct a separate incubation control by incubating both oligos separately with extracts, extract the DNA, mix, and then perform PCR detection to ascertain that bands labeled as products are formed during the extract incubation and not PCR artifacts. The recombination literature is replete with this artifact, and it must be eliminated here.

The separate incubation controls are also needed for the experiments with mt plasmid DNA and mitochondrial genomics DNA in Figure 5F and Figure S6.

In their exploration of the mechanisms leading to deletions, the data presented supports their model invoking MMEJ, at least for the in vitro system as shown by the sequencing of the deletions obtained either with the plasmid or with mitochondrial DNA. When mitochondrial DNA was used, however, it would seem that mutations in that region accumulate but they are not associated with deletions, which are proportionally less abundant than in the plasmid system. This finding merits more discussion.

The authors hypothesize that MMEJ repair of the DNA break caused by the endonuclease G produces the 9 bp-deletion. If this is the case, why just one out of four repeats is always deleted? Further, the authors may want to tone down their MMEJ model, since no direct evidence for this mechanism is provided. Other mechanisms, such as nick-repair can potentially lead to the same outcome.

7) Page 13, lines 20/21: Give the exact numbers for 'few' and 'many others'.

8) The EndoG nuclease experiments and conclusion that EndoG is responsible for the observed cleavages, would be strengthened by purifying a catalytic mutant of EndoG and showing dependance of cleavage on the EndoG active site. This is a standard control for the nuclease field.

9) The discussion is not well organized and should include an evaluation what the possible physiological function of the observation of EndoG release upon oxidative stress and deletion formation might be. Are these off-pathway accidents?

Is there a role of EndoG on G4 RNA tetraplexes? This refers also to the point about the BG4 antibody; does it recognize G4 RNA?

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

Thank you for resubmitting your work entitled "Unleashing a Novel Function of Endonuclease G in Mitochondrial Genome Instability" for further consideration by eLife. Your revised article has been evaluated by Jessica Tyler (Senior Editor) and a Reviewing Editor.

The manuscript has been improved but there are some remaining issues that need to be addressed, as outlined below:

The manuscript makes interesting observations that should advance the field of mitochondrial DNA metabolism and will be off interest for researchers in the field of mitochondrial genome stability and mitochondrial diseases.

The revised manuscript is much improved, although there are number of points that should be addressed in a final revision.

Essential revisions:

1) The manuscript still requires some language editing.

2) Page 6, line 33: Overstatement, there are some pause site w/o KCl; change to significantly fewer sites and less pausing.

3) Figure 2 D: Primer DI12 is too inefficient to be meaningful. A new primer would be needed to make the point. The question is whether these data are needed in light of the mutant control?

4) Page 10, line 24-27 and Figure 6E: Make an additional mutant to test the prediction of the origin of the second remaining band.

5) Some points concerning the localisation of the protein in the mitochondria sub-fractions and its proposed recruitment to the matrix are still not convincing (please see specific comments below). There are enough novel data to justify its publication even without including the above-mentioned aspects, or if they are only presented as preliminary data to propose a hypothesis. The discussion should be amended to reflect that.

• Line 372: The fact that EndoG is present in mitochondria was already documented by several papers. It is not clear where the novelty is in this Results section. Furthermore, the title mentions TFAM as an inner membrane protein, which is not the case. As mentioned in the text, although TFAM is found associated with the membrane, it is a matrix protein. The low level of co-localisation and the limited resolution of the microscope used make do not allow to conclude much beyond the already known presence of Endo G in mitochondria. More detailed localisation studies would require electron microscopy or super resolution microscopy.

• Line 396: Based on the same arguments used above, the results described in the first paragraph of this section do not add much. The resolution of the images does not allow to conclude the presence of Endo G at G-quadruplexes through the co-localisation with the BG4 antibody. However, the PCR experiments provide evidence for that.

• The sub-mitochondrial fractionation described does not correspond to the references cited. In both papers cited there is a first step is a hypotonic shock that would yield a pellet of mitoplasts that can then be separated into inner membrane and matrix fractions. Here the fractionation is directly done in high salt and, to my understanding there is no an intermediate step allowing to get rid of the intermembrane fraction or the outer membrane.

There is all along an ambiguity regarding the localisation of EndoG. As mentioned by the authors the protein is normally found in the intermembrane space, therefore, the majority of it should be found in the soluble fraction with their fractionation protocol or lost if they use the hypotonic one. However, they find a substantial fraction in the membrane (inner or total?) fraction and this fraction does not decrease upon treatment while the soluble does. The question is then, where does the newly soluble protein come from?

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

Thank you for resubmitting your work entitled "Unleashing a Novel Function of Endonuclease G in Mitochondrial Genome Instability" for further consideration by eLife. Your revised article has been evaluated by Jessica Tyler (Senior Editor) and a Reviewing Editor.

The manuscript has been improved but there are some remaining issues that need to be addressed, as outlined below:

1) There is still a contradiction in the section describing the localisation of EndoG within the mitochondria. The title of that section (page 11) is "Endonuclease G is expressed in mitochondria and colocalises with inner mitochondria membrane protein TFAM within the cells". However, based on the experiments described in the text that follows and shown in Figure 8, the authors conclude that it colocalises with TFAM which they define as "a transcription factor known to be present in the mitochondrial matrix".

2) As described the protocol for mitochondrial fractionation still does not provide a separation of the inner membrane from the outer one, nor does it separate the intermembrane space from the matrix (unless the description omitted a centrifugation step after the hypotonic shock). If this is the case the authors should be more careful in their statements and mention in the text on page 13, as they did in Figure 10, that what they are looking at is the partial release from the membrane to the soluble fractions.

eLife. 2022 Nov 17;11:e69916. doi: 10.7554/eLife.69916.sa2

Author response


Essential revisions:

The manuscript makes interesting observations that should advance the field of mitochondrial DNA metabolism and will be off interest for researchers in the field of mitochondrial genome stability and mitochondrial diseases.

The revision that is needed to make this manuscript, as submitted, publishable is too extensive. If the authors decide to address the comments by the additional experimentation needed, they are invited to resubmit a revised manuscript as a new submission that I would be happy to handle as BRE and seek to have the same reviewers evaluate the new submission.

Alternatively, the authors may want to consider tightening the manuscript and omit the reconstitution studies making sure that their interpretations are supported by the data and consider alternative mechanisms, when necessary. This more focused version would still contain sufficient new insights and allow the authors to further develop the reconstitution experiments with the required exhaustive separate incubation controls for a future publication.

We thank the editors and reviewers for their positive comments and suggestions on our manuscript. We have omitted reconstitution experiments as suggested and tightened the manuscript by reorganizing the figures and text to address all the comments, as detailed below.

1) Independent of specific criticism towards individual experiments or experimental tools, the manuscript is difficult to read, poorly assembled and would need a major overhaul in the text and figures. The text requires careful editing, there are many disagreements, and many cases where referrals such as 'its' or 'this' are unclear. The figures should be labeled that the experiments can be understood independent of the legend or method description. A few examples are given here.

• Figure 1A: The sequences of C1 and G1 are needed here. Giving more is confusing, as they do not relate to the figure. In Figure S1, all oligo sequences should be listed.

As suggested by the reviewer, we have removed the oligomer list of mutants from Figure 1 (now named as Figure 1B) and moved the mutants to Supplement Figure 1 (Figure 1—figure supplement 1B and Supplement Table 1).

• Figure S1A: Define -/+ KCl under gel in legend.

The suggested label is mentioned below the gel and added also to Figure legend (new Figure 1C and Figure 1—figure supplement 1C) and (page 33, lines 5-9 as well as page 40, lines 10-12).

• Figure S1B: Indicate inverted repeat sequences.

As suggested by the reviewer, we have now indicated the inverted repeat sequences (new Figure 1A and Figure 2—figure supplement 1).

• Figure S3C: Where are the gel data? What is shown here? The description in the legend does not allow to understand the figure, what is the difference between the 3 lanes that are labeled as +.

We are extremely sorry for the confusion. The gel data and the quantitation are now presented together in the same figure to avoid confusion (New Figure 6E, F). The 3 lanes labelled as + are the triplicate samples, which is now explained in the legend (page 36, lines 5-10).

• Figure S3: Unclear what the Westernblots show in D, F, and H.

The western blots in D, F, H confirm the identity of purified proteins in the fractions used (New Figure 7—figure supplement 1). Details are added to the legend (Page 41, lines 16-26).

• Figure S5F: Define the lanes. What is loaded in lane 4?

In old figure S5F (New Figure 9—figure supplement 2A), Lane 1 is marker, lanes 2 and 3 are pull down samples loaded after incubation of secondary antibody and anti-endonuclease G, respectively with mitochondrial extract. Lane 4 is mitochondrial extract.

• Figure 4D: Label upper input and lower pull-down, also Figure 2G.

As suggested, we have labelled the figures (new Figures 5C, 9D, Figure 9—figure supplement 1D).

2) A key problem in the manuscript is that the authors mix their interpretation into the description of the results. A few examples are given here, but the authors should check the entire manuscript and correct this issue. The authors should carefully separate description of the results from their interpretation and refrain from overinterpreting their data (see also below). There is no need to overstate the result; exceptions or unexpected results must be discussed and acknowledged. For example:

• Figure 2F/G: What happened with CR7? Low threshold value but no PCR product? This is not described in the text.

When sequence of the region spanning CR7 was analysed for potential G4 motifs, we observed sequence that may support formation of 2 plate G4 structure. This could explain the observed binding of BG4 during semiquantitative PCR. However, the threshold value was based on real time PCR which is highly sensitive when compared to semiquantitative PCR. This could be the reason for absence of a band in the CR7. This is also described in the text now (Page 9, lines 11–18).

• Figure 3D: The upper cleavage site stays the same, unlike what is described in the text.

The mitochondrial G-quadruplex region I has the ability to form G-quadruplexes in multiple forms (Dahal et al., 2022). Mutating two of the sequences will still leave with 3-G stretches and 3 GNG sequences which makes the region, a ligand for Endonuclease G binding and cleavage (Figure 6—figure supplement 1). This is described in the main text (Page 11, lines 14-20).

• Figure S3A, B: It is unclear which band is EndoG. Is the Westernblot aligned with material loaded in slots 4-7 of the protein-stained gel? If so, why does the Westernblot intensity not reflect the staining intensity?

We have repeated the experiment and provided a new western blot for showing the presence of Endo G in the purified protein. We have marked the band representing Endo G in the figure. For both Endo G and its mutant, fractions 4 and 5, respectively were used for blotting (New Figure 7—figure supplement 1A-C).

• Figure 3F: The Mt extract does only show the major EndoG sites but not the minor sites, which add up to a significant proportion of the products.

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, there may be additional proteins, which can bind to G4 DNA, thus preventing cleavage at multiple sites. This has been now added in Page 11, lines 14-20.

• Page 10, line 10: The mutant plasmid still has one cleavage site, cleavage is not abrogated as stated in the text.

Although, 2 of the G-stretches are mutated, Region I still can fold into a G-quadruplex, thus one of the cleavage sites is not abrogated. We have corrected the statement in the mentioned paragraph (Page 10, line 22-26).

• Page 10, line 22: Figure 3D is extract not purified EndoG as stated in line 22 of page 10.

We apologize for the mistake. We have now corrected the statement (Page 10, lines 5-6, New figure 6).

• Figure 4E and F: What is the difference between both experiments? E: extract, F: cells? This could be clarified in the figure and text.

We apologise for the confusion. Figure 4E (New Figure 9E) corresponds to the use of mitochondrial extract and 4F (New Figure 9F) refers use of cells. The same is included in the text as well (page 13, lines 9-19).

• Figure S5G: The description of the P1 results is not precise. There is clear protection of the C strand by EndoG in the upper third of the ladder.

As suggested we have edited the text to match with the results in the Figure (New Figure 9—figure supplement 2B) (page 13, lines 20-25).

• Figure 7A requires more description and quantitation. Figure 7 B/E: Why is the DAPI signal so variable? Figure 7G: the difference between lanes 1-3 and 4-6 should become apparent from the figure. For Figure 7G and H, only menadione data are shown, there are no IR data as stated in the text.

We have now added description and quantitated figure 10A-D. We have also changed the images in figure 10B/E and the DAPI signal is now comparable between the panels. We apologise for the mistake. As mentioned, figure 10G and H is only for Menadione and not for IR.

3) The BG4 antibody is problematic. Page 8, line 22: Why is there a BG4 signal in the cytoplasm in rho (0) cells? Is this also seen in rho (+) cells?

BG4 binds to G-quadruplexes of nuclear DNA and RNA in cytoplasm (Biffi et al., 2014; Biffi et al., 2013). All the BG4 outside nucleus is not only because of G4 DNA in mitochondria, but also due to RNA G4 in the cytoplasm (Biffi et al., 2014). We have used mitotracker to indicate the mitochondria. Recently, we have also extensively investigated the specificity of BG4 towards G4 DNA (Javadekar et al., 2020). The BG4 signal in the cytoplasm is also seen in rho (+) cells.

Also, the authors may want to discuss in more details image modification via Coste's mask. I can clearly see yellow dots in merged images for rho (0) cells, but somehow they are not converting to the white dotes when the Coste's mask is applied.

We have now replaced the images to take care of the reviewer’s comments (Figure 4D).

Figure 2 and 4: There are many BG4 signals? Do the signals all represent G4 tetraplex DNA? How can the authors ascertain the specificity and usefulness of this reagent?

The control in Figure 4A/B needs to be defined. What would be expectation from random overlap in this experiment given that both signals are very abundant?

Approximately 7,00,000 canonical G4 forming motifs of 25-35 nt are reported in the human nuclear genome (Chambers et al., 2015). Besides, recently 2 plate G-quadruplexes and involvement of GNG motifs were reported during G-quadruplex structure formation (Biffi et al., 2013; Das et al., 2016). If one considers all these into account, ~2% of the genome may support the formation of G-quadruplexes in human genome. Thus, the number of BG4 signals can be justified in the nucleus. Within mitochondria 5 canonical G4 motifs (Dahal et al., 2022) and >200 number of 2 plate G-quadruplexes (Brázda et al., 2019) have been reported. Considering that several mtDNA are present in one mitochondrion and several mitochondria are present in one cell, impact can be several fold. In addition to that RNA G4 structures have also been reported in cytoplasm (Biffi et al., 2014; Biffi et al., 2013).

The specificity and usefulness of the reagent has been extensively shown in the previous studies (Biffi et al., 2013; Das et al., 2016; Javadekar et al., 2020)

4) The microscopy approach is less informative and overinterpreted. The level of resolution obtained with the confocal microscope approach used does not allow a conclusion on the localization beyond saying that the protein is present in the mitochondria and is not suitable for distinguishing the different mitochondrial compartments.

We agree to the reviewer’s comments that the conclusion on localization based on microscopy approach is not preferred, and thus apart from microscopy we have used the fractionation of mitochondrial compartments to support the data (new Figure 10G, H).

As for the fractionation procedure used, I am not sure that it allows the separation of the matrix fraction from the intermembrane compartment. Furthermore, the description of the results displayed is not clear and there seem to be some mistakes in the labelling of the figures, making it harder to draw clear conclusions.

The fractionation procedure has been used as described in the previous literature (Sinha et al., 2010; Suzuki et al., 2002). We apologise for the mistakes in labelling and now it has been corrected and presented (new Figure 10).

5) The mutant substrates are an excellent control for testing the G4 tetraplex hypothesis. What is the experimental evidence that these mutant oligonucleotides do not form tetraplexes?

The mutant substrates have been used for the study to understand their mobility along with the wild type oligomers in the gel shift assays (Figure 1—figure supplement 1B, C). Apart from this, circular dichroism was performed and shown previously for this region and its mutants that they do not fold into G4 DNA (Dahal et al., 2022). Besides, we have also investigated the same region (with and without mutation) after cloning in a plasmid (new Figure 2).

6) The attempt to reconstruct the deletion formation in vitro is welcomed. However, these experiments lack a crucial control. In Figure 5B, there is product in the negative control? This raises the concern that bands labeled MMEJ products are PCR artifacts. The authors must conduct a separate incubation control by incubating both oligos separately with extracts, extract the DNA, mix, and then perform PCR detection to ascertain that bands labeled as products are formed during the extract incubation and not PCR artifacts. The recombination literature is replete with this artifact, and it must be eliminated here.

The separate incubation controls are also needed for the experiments with mt plasmid DNA and mitochondrial genomics DNA in Figure 5F and Figure S6.

In their exploration of the mechanisms leading to deletions, the data presented supports their model invoking MMEJ, at least for the in vitro system as shown by the sequencing of the deletions obtained either with the plasmid or with mitochondrial DNA. When mitochondrial DNA was used, however, it would seem that mutations in that region accumulate but they are not associated with deletions, which are proportionally less abundant than in the plasmid system. This finding merits more discussion.

The authors hypothesize that MMEJ repair of the DNA break caused by the endonuclease G produces the 9 bp-deletion. If this is the case, why just one out of four repeats is always deleted? Further, the authors may want to tone down their MMEJ model, since no direct evidence for this mechanism is provided. Other mechanisms, such as nick-repair can potentially lead to the same outcome.

As discussed above, based on the suggestion by editors and reviewers we have now removed reconstitution studies from the manuscript. All suggested controls will be performed and a separate manuscript will be submitted at a later timepoint. We sincerely thank the editors for the constructive criticism and for the suggestion.

7) Page 13, lines 20/21: Give the exact numbers for 'few' and 'many others'.

We have now removed this figure from the manuscript completely and edited the text accordingly.

8) The EndoG nuclease experiments and conclusion that EndoG is responsible for the observed cleavages, would be strengthened by purifying a catalytic mutant of EndoG and showing dependance of cleavage on the EndoG active site. This is a standard control for the nuclease field.

We thank the editor and reviewer for the suggestion. We have now purified the catalytic mutant of Endo G and shown the cleavage pattern when the mutant protein is used (new Figure 7D, E).

9) The discussion is not well organized and should include an evaluation what the possible physiological function of the observation of EndoG release upon oxidative stress and deletion formation might be. Are these off-pathway accidents?

­­­­­As suggested, we have now reorganized the discussion. The observation that EndoG is released into matrix upon stress is new in the field. Although our data suggest that it may play a role in deletion formation, further studies are required to demonstrate this (Page 16, lines 13-17; page 17, lines 15-20).

Is there a role of EndoG on G4 RNA tetraplexes? This refers also to the point about the BG4 antibody; does it recognize G4 RNA?

Previously, it has been reported that BG4 binds to G4 RNA (Biffi et al., 2014; Biffi et al., 2013; Brázda et al., 2019). However, affinity of Endo G on G4 RNA tetraplexes is yet to be investigated. A detailed study is required before we can comment on this aspect.

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

The manuscript has been improved but there are some remaining issues that need to be addressed, as outlined below:

The manuscript makes interesting observations that should advance the field of mitochondrial DNA metabolism and will be off interest for researchers in the field of mitochondrial genome stability and mitochondrial diseases.

The revised manuscript is much improved, although there are number of points that should be addressed in a final revision.

We thank the editors for their overall positive comments on our manuscript. Based on their suggestions following changes have been incorporated into the manuscript as detailed below.

Essential revisions:

1) The manuscript still requires some language editing.

We have gone through the manuscript and made additional language corrections as suggested.

2) Page 6, line 33: Overstatement, there are some pause site w/o KCl; change to significantly fewer sites and less pausing.

As suggested, we have modified the text (Page 6, line 33-35).

3) Figure 2 D: Primer DI12 is too inefficient to be meaningful. A new primer would be needed to make the point. The question is whether these data are needed in light of the mutant control?

As rightly pointed out by the reviewer, the efficacy of primer (DI12) binding is inefficient, and this data does not add up significantly to the overall content of Figure 2. Thus, we have removed the panel from figure 2 and reorganised the rest of the figure (Revised Figure 2).

4) Page 10, line 24-27 and Figure 6E: Make an additional mutant to test the prediction of the origin of the second remaining band.

As suggested by the reviewer, we have now generated a new mutant to test the predicted basis of the second band and presented it as Figure 6-supplement figure 1. Results showed that upon generation of additional mutants, the second band, due to cleavage, disappeared. This further confirms the involvement of G4 DNA during the observed nicking (Page 10, lines 23-26).

5) Some points concerning the localisation of the protein in the mitochondria sub-fractions and its proposed recruitment to the matrix are still not convincing (please see specific comments below). There are enough novel data to justify its publication even without including the above-mentioned aspects, or if they are only presented as preliminary data to propose a hypothesis. The discussion should be amended to reflect that.

We agree with the reviewer, and as suggested, the immunolocalisation data showing the recruitment of Endonuclease G to the matrix is now removed from Figure 10. However, we have maintained sub-fractionation data as a preliminary observation (Figure 10).

• Line 372: The fact that EndoG is present in mitochondria was already documented by several papers. It is not clear where the novelty is in this Results section. Furthermore, the title mentions TFAM as an inner membrane protein, which is not the case. As mentioned in the text, although TFAM is found associated with the membrane, it is a matrix protein. The low level of co-localisation and the limited resolution of the microscope used make do not allow to conclude much beyond the already known presence of Endo G in mitochondria. More detailed localisation studies would require electron microscopy or super resolution microscopy.

We agree with the reviewer’s comments that the conclusion on sub-localisation based on the standard microscopy technique may not be sufficient to make the point. Thus, we have removed the immunofluorescence figures related to the localisation of EndoG to different mitochondrial compartments.

• Line 396: Based on the same arguments used above, the results described in the first paragraph of this section do not add much. The resolution of the images does not allow to conclude the presence of Endo G at G-quadruplexes through the co-localisation with the BG4 antibody. However, the PCR experiments provide evidence for that.

We thank the reviewer for their comment. We agree that the results described do not add much and the presence of BG4 cannot be concluded based on only colocalization with BG4. However, as stated, since we have other experiments (PCR-based ChIP) to evaluate the same question, we have maintained the data in the manuscript (Figures 4 and 9A).

• The sub-mitochondrial fractionation described does not correspond to the references cited. In both papers cited there is a first step is a hypotonic shock that would yield a pellet of mitoplasts that can then be separated into inner membrane and matrix fractions. Here the fractionation is directly done in high salt and, to my understanding there is no an intermediate step allowing to get rid of the intermembrane fraction or the outer membrane.

We apologise for not providing a detailed protocol. As mentioned in Sinha et al., we first performed the hypotonic shock to yield mitoplasts and separated them into inner membrane and matrix fractions (Sinha et al., 2010, Hum Mol Genet, vol 19, page 1248-1262; Suzuki et al., 2002, J Cell Sci, vol 115, pages 1895-1905). Accordingly, we have rewritten the section by providing sufficient detail to avoid ambiguity (Page 20, lines 24-32).

There is all along an ambiguity regarding the localisation of EndoG. As mentioned by the authors the protein is normally found in the intermembrane space, therefore, the majority of it should be found in the soluble fraction with their fractionation protocol or lost if they use the hypotonic one. However, they find a substantial fraction in the membrane (inner or total?) fraction and this fraction does not decrease upon treatment while the soluble does. The question is then, where does the newly soluble protein come from?

We agree with the reviewer that there has been ambiguity in the sublocalization of EndoG. While there are few papers suggesting the localization to inner membrane (David et al., 2006 Cell Death Differ;13(7):1147-55,Ishihara and Shimamoto (2006) JBC; 281, (10): 6726 –6733) others suggest its presence in the inner membrane space (Ohsato et al., 2002 Eur. J. Biochem.269, 5765–5770, Zhou et al., 2016 Science 353 (6297):394-9, Davies et al., 2003 Nucleic acids Res. 31(4): 1364–1373). With the protocol that we have used (the hypotonic method), if it was present only in the inner membrane space, Endo G would have been completely lost. However, further investigation is required to comment on this aspect.

We apologise for the ambiguity of quantification. We have now re-quantified the western blots and presented it as Figure 10B.

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

The manuscript has been improved but there are some remaining issues that need to be addressed, as outlined below:

1) There is still a contradiction in the section describing the localisation of EndoG within the mitochondria. The title of that section (page 11) is "Endonuclease G is expressed in mitochondria and colocalises with inner mitochondria membrane protein TFAM within the cells". However, based on the experiments described in the text that follows and shown in Figure 8, the authors conclude that it colocalises with TFAM which they define as "a transcription factor known to be present in the mitochondrial matrix".

We thank the reviewer for pointing out this and apologize for the error. We have now corrected the title section on page 12, lines 1, 2.

2) As described the protocol for mitochondrial fractionation still does not provide a separation of the inner membrane from the outer one, nor does it separate the intermembrane space from the matrix (unless the description omitted a centrifugation step after the hypotonic shock). If this is the case the authors should be more careful in their statements and mention in the text on page 13, as they did in Figure 10, that what they are looking at is the partial release from the membrane to the soluble fractions.

Thank you for the suggestion. We have now edited the text to avoid the discrepancy (page 13, line 31).

Associated Data

    This section collects any data citations, data availability statements, or supplementary materials included in this article.

    Supplementary Materials

    Figure 1—source data 1. EMSA and CD studies to show formation of G-quadruplex structure in Region I.
    Figure 1—source data 2. Gel shift assay to show abrogation of G-quadruplex structure in mutants of Region I.
    Figure 2—source data 1. Primer extension studies to show formation of G-quadruplex structure in Region I containing plasmid.
    Figure 3—source data 1. Sequence of clones after bisulphite treatment in a plasmid containing Region I of mitochondria.
    Figure 3—source data 2. Sequence of clones after bisulphite treatment in a region I of mitochondria.
    Figure 4—source data 1. Localization of BG4 to mitochondria.
    Figure 4—source data 2. Immunofluorescence showing the localization of BG4 to mitochondria.
    Figure 5—source data 1. BG4 ChIP to show the binding of BG4 to mitochondrial G-quadruplex forming regions.
    Figure 6—source data 1. Gel profiles showing the mitochondrial induced cleavage at mitochondrial region I.
    Figure 6—source data 2. Cleavage assay on plasmid bearing wildtype and mutant G4 sequence.
    Figure 7—source data 1. Gel profiles showing the Endonuclease-G-induced cleavage at mitochondrial region I.
    Figure 7—source data 2. Gel profiles showing the purification of different endonucleases.
    Figure 8—source data 1. Localization of Endonuclease G to mitochondria.
    Figure 8—source data 2. Immunofluorescence showing the Localization of Endonuclease G to mitochondria.
    Figure 9—source data 1. ChIP assay showing the binding of Endonuclease G with the mitochondrial G-quadruplex regions within cells.
    Figure 9—source data 2. ChIP assay showing the binding of Endonuclease G with the mitochondrial G-quadruplex regions when purified Endonuclease G was used.
    Figure 9—source data 3. P1 nuclease assay showing the binding of Endonuclease G to mitochondrial G quadruples regions.
    Figure 10—source data 1. Sub localization of Endonuclease G with or without induction of stress.
    Transparent reporting form

    Data Availability Statement

    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.


    Articles from eLife are provided here courtesy of eLife Sciences Publications, Ltd

    RESOURCES