Unleashing a novel function of Endonuclease G in mitochondrial genome instability
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.
Editor's evaluation
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.
https://doi.org/10.7554/eLife.69916.sa0Introduction
The mitochondria are semiautonomous organelles with their genome. The mammalian mitochondrial genome encodes for 13 respiratory chain proteins, 22 tRNAs and 2 rRNAs (Chen and Butow, 2005; Clayton, 1984; Holthöfer et al., 1999). Due to its vital role in oxidative phosphorylation, mitochondrial DNA is frequently exposed to reactive oxygen species (ROS) and is more prone to DNA damage than nuclear DNA, leading to the accumulation of 50 times more mutations, including deletions (Chen et al., 2011; Hudson et al., 1998; Michikawa et al., 1999; Pakendorf and Stoneking, 2005; Yakes and Van Houten, 1997). Mitochondrial mutations have been associated with ageing and several disease conditions like myopathies, dystonia, cancer etc. (Chen et al., 2011; Taylor and Turnbull, 2005). These mutations include point mutations, mismatches, or deletions, which affect the coding of essential proteins involved in oxidative phosphorylation and the respiratory chain (Penta et al., 2001). Mitochondrial mutations, especially large-scale deletions of 200–5000 bp, were associated with the progression of breast, colorectal, renal, and gastric cancers (Bianchi et al., 2001; Carew and Huang, 2002; Modica-Napolitano et al., 2007; Penta et al., 2001). However, the correlation between mtDNA deletions and cancer remains unsettled. The well-studied phenotypes of mtDNA deletions are Kearns-Sayre Syndrome (KSS), Pearson syndrome, and progressive external ophthalmoplegia (PEO; Marni, 2003). Large stretches of deletions seen in different diseased patients suggest fragile sites in the mitochondrial genome similar to the fragile sites (FRA) sites of the nuclear DNA (Richards and Macaulay, 2001).
The FRA regions and other frequent chromosomal translocations have been observed in late replicating sites in the nuclear DNA, which are prone to replication fork stalling due to possible, stable secondary structure formation or inadvertent action of structure-specific nucleases (Dillon et al., 2010; Lin et al., 2013; Raghavan et al., 2004; Sun et al., 2001). Mitochondrial DNA has its replication machinery that is not governed by cell cycle checkpoints. Therefore, it is likely that the formation of a secondary structure may account for observed mitochondrial fragility.
The presence of alternate DNA structures like G-quadruplex, triplex, and cruciform structures in the mitochondrial genome has been reported along with their association with different mitochondrial diseases (Bharti et al., 2014; Damas et al., 2012; Dong et al., 2014; Oliveira et al., 2013). Breakpoint regions reported in the case of PEO patients were in the proximity of G-quadruplex motifs, further suggesting the role of non-B DNA structures in mitochondrial disorders (Bharti et al., 2014). G4 forming motifs were also located adjacent to mitochondrial deletion breakpoints associated with KSS, a clinical subgroup of mitochondrial encephalomyopathies (Van Goethem et al., 2003; Zeviani et al., 1988). Further, a G4 ligand RHPS4 preferentially localises to mitochondria and affects mitochondrial genome content, transcription elongation, and respiratory function (Falabella et al., 2019a; Oliveira et al., 2013). Moreover, in vitro results suggested that a single-base mutation in mtDNA (mt10251) was sufficient to form G4 DNA (Chu et al., 2019). Importantly, ATP-dependent G4 resolving helicase Pif1 could unfold non-B DNA structures and continue the synthesis by mitochondrial DNA polymerases, Pol γ and PrimPol (Butler et al., 2020). However, further studies are required to confirm the presence of these structures in vivo and establish their contribution towards the mechanism of breakage and deletions.
Mitochondrial DNA (mtDNA) 9 bp deletion caused by loss of one copy of the 9 bp repeat sequence (CCCCCTCTA) in the intergenic region of cytochrome c oxidase II (COII)/tRNALys is one of the most common deletions in mitochondria (Redd et al., 1995; Thomas et al., 1998). It has been associated with various diseases, including hepatocellular carcinoma (Jin et al., 2012; Komandur et al., 2011; Krishnan and Turnbull, 2010; Zhuo et al., 2010). The 9 bp deletion is also widely used as a phylogenetic marker to study evolutionary trends and migration of populations (Soodyall et al., 1996; Yao et al., 2000). The occurrence of deletions in such high frequency indicates the fragility of the region; however, its mechanism is unclear.
Out of the large-scale mtDNA deletions reported, one of the well-studied examples is mt4977, and accumulation of it results in different mitochondrial disorders, including cancer (Yusoff et al., 2019). mt4977 eliminates mitochondrial sequences between 8470 and 13,447 bp, which include the removal of genes required for normal oxidative phosphorylation (OXPHOS) (Lee et al., 1994; Yang et al., 1994; Yen et al., 1991). Recently the role of replication slippage and DSB repair has been shown to play a role in the generation of mtDNA4977 in mammalian cells (Phillips et al., 2017). Previous studies have shown that the presence of a 13 bp direct repeat plays a major role in mt4977, as only one of the repeats was retained following joining, indicating the involvement of microhomology-mediated end joining (Samuels et al., 2004; Srivastava and Moraes, 2005). When restriction endonucleases were used to introduce DSBs, there was the accumulation of deletions, which further implicates the imperfect repair and deletion in the mitochondrial genome (Srivastava and Moraes, 2005). A previous study from our laboratory suggested the existence of microhomology-mediated end joining (MMEJ) in mitochondria, which operate in a microhomology-dependent manner, explaining the frequently seen deletions in mitochondrial DNA (Tadi et al., 2016).
Here, we investigate the mechanism of fragility associated with the most common deletion seen in the mitochondrial genome. The putative role of non-B DNA structures and the mechanism of DNA breakage associated with the formation of mtDNA deletions are extensively studied. We report that the region associated with 9 bp deletion can fold into a G-quadruplex DNA structure using various biochemical and ex vivo methods. For the first time, we show that Endonuclease G can specifically bind to the G-quadruplex structure and cleave the mtDNA. Finally, the joining of the broken DNA could lead to the formation of ‘9 bp deletion’ seen in human mitochondria.
Results
Previously, different non-B DNA motifs, including G4 DNA motifs have been reported in the mitochondrial genome (Dahal et al., 2022; Damas et al., 2012; Dong et al., 2014). It was also suggested that alternate forms of DNA may play a role in the generation of deletions in mitochondria (Dahal et al., 2022; Damas et al., 2012; Dong et al., 2014). In silico studies along with biochemical analyses revealed that 5 classical G4 DNA motifs are present in the mitochondrial genome, which could fold into G4 DNA (Cer et al., 2013; Dahal et al., 2022; Figure 1—figure supplement 1A). Since ‘9 bp deletion’ (8271–8281) is the most frequent mitochondrial deletion and occurs next to one of the 5 G4 motifs with the coordinate position 8252–8295 (designated as mitochondrial Region I) (Figure 1A and B; Dahal et al., 2022), we have investigated the putative role of G4 DNA structure in DNA breakage and formation of the mtDNA deletions. Notably, the chosen G-quadruplex motif is also present in two other large-scale deletions associated with mitochondria.
G-quadruplex motifs from Region I of mitochondrial DNA, when present on shorter DNA, can fold into a G-quadruplex structure
To investigate whether mitochondrial Region I can indeed form G-quadruplex structures, oligomeric DNA containing the G-rich region referred to as ‘G1’ and complementary C strand (C1) were synthesised (Figure 1B). Radiolabelled oligomeric DNA was then subjected to electrophoretic mobility gel shift assay in the absence or presence of KCl. Faster mobility of the G-rich strand compared to its complementary C-rich strand suggested the formation of both intra and intermolecular G4 DNA structure in a KCl-dependent manner (Figure 1C). Importantly, the formation of intramolecular G4 DNA was abrogated when G-stretches were mutated (Figure 1—figure supplement 1B, C). Further, circular dichroism (CD) studies revealed that the spectra exhibited the characteristic of a parallel G-quadruplex DNA with absorption maxima of 260–270 nm and absorption minima of 240 nm as opposed to the C-rich strand, which showed the maxima similar to that of a typical single-stranded DNA (Figure 1D).
DMS protection assay was performed to understand the precise base pairing of guanine nucleotides involved in the G4 plates compared to the ones within connecting loops that facilitate DNA strand folding and, thus, G-quadruplex formation (Nambiar et al., 2011; Nambiar et al., 2013). DMS methylates the N7 position of the guanine base, making it susceptible to subsequent piperidine cleavage. Since the N7 position of G residues forming the G4 plate is engaged upon Hoogsteen base pairing, these residues are protected from DMS treatment and piperidine cleavage. When the G-quadruplex sequences from mitochondrial Region I (G1) were subjected to DMS protection assay, following structure formation, we observed that most G residues in this region showed a significant reduction in cleavage, suggesting protection of these nucleotides from DMS. Based on the observed protection of guanine from DMS, we conclude that these guanines are indeed involved in Hoogsteen bonding (Figure 1E, lane 2). The results suggest that the guanine residues from the G stretches of the mitochondrial Region I motif participate in G-quadruplex formation. Based on CD studies and gel shift assays, a parallel, intramolecular G-quadruplex structure was modelled (Figure 1F). Thus, our study reinforces that the Region I of mitochondrial DNA can fold into an intramolecular G-quadruplex (Figure 1E and F).
The G4 motif present in Region I of mtDNA can fold into a G-quadruplex structure when present on a plasmid DNA
In the context of double-stranded DNA, G-quadruplexes can form only when the Hoogsteen base pairing between the guanine residues is stronger than the Watson and Crick H-bonding required for duplex B-DNA formation. Therefore, the G-quadruplex forming ability of the mitochondrial Region I was tested after cloning into a plasmid, pDI1. Interestingly, an inverted repeat sequence (coordinate positions 8295–8363), which could fold into a cruciform DNA, was also present adjacent to the G4 motif. Thus, that was also part of pDI1 (Figure 2—figure supplement 1). Further, to delineate the role of G4 DNA in the plasmid context, a mutant plasmid with a scrambled sequence corresponding to the G4 motif but with intact inverted repeats (pDI2) was also constructed (Figure 2A, Figure 2—figure supplement 1). Previous studies have shown that the formation of G4 DNA can lead to a pause during replication (Kumari et al., 2015; Nambiar et al., 2013; Voineagu et al., 2008; Wells, 2009). To assess whether the G4 DNA formed at Region I of the mitochondrial genome can act as a replication pause, a primer extension assay was performed using pDI1 (Kumari et al., 2015; Figure 2B and C). A sequencing ladder was also generated for mapping the location of the polymerase pause sites on the plasmid (Figure 2C). Results revealed multiple polymerase pause sites, and the intensity of the pauses increased with increasing concentrations of KCl (25, 50 and 100 mM; Figure 2C). Interestingly, there was significantly fewer sites and less pausing observed in the absence of KCl, indicating the absence of the formation of intramolecular G-quadruplex species in these cases (Figure 2C, lane 2). Considering that the formation of cruciform DNA occurs irrespective of the presence of KCl, these results also suggest that the observed pause was due to the formation of G4 DNA, but not cruciform DNA. Besides, the observed polymerase pauses mapped to the G-quadruplex motif (Figure 2C, lanes 7–9). Abrogation of G4 DNA structure by mutation of G4 DNA motif (pDI2) resulted in the disappearance of the observed pause site, further confirming the above results (Figure 2D, lanes 3, 4 and lane 6, 7). It was also observed that the appearance of strong polymerase pause sites were dependent on K+ (Figure 2E, lanes 4–5; F), but not on other monovalent cations such as Li+ (with no or little G-quadruplex stabilising properties) or Na+ (with limited stabilising properties) (Figure 2E, lanes 6–9; F) suggesting that the non-B DNA responsible for polymerase arrest was G4 DNA.
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.
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).
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).
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.
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).
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).
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—source data 1
- https://cdn.elifesciences.org/articles/69916/elife-69916-fig9-data1-v3.zip
-
Figure 9—source data 2
- https://cdn.elifesciences.org/articles/69916/elife-69916-fig9-data2-v3.zip
-
Figure 9—source data 3
- https://cdn.elifesciences.org/articles/69916/elife-69916-fig9-data3-v3.zip
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.
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).
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
Request a detailed protocolChemicals 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
Request a detailed protocolThe 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.
Cell lines
Request a detailed protocolHeLa (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
Request a detailed protocolshRNA 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
Request a detailed protocolMitochondrial 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
Request a detailed protocolThe 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
Request a detailed protocolFor 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
Request a detailed protocolMitochondria 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
Request a detailed protocolMitochondrial 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
Request a detailed protocolA 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
Request a detailed protocolThe 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
Request a detailed protocolGel 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)
Request a detailed protocolCircular 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
Request a detailed protocolThe 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
Request a detailed protocolPlasmid 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
Request a detailed protocolWild-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
Request a detailed protocolSequencing 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
Request a detailed protocolpDI1, 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
Request a detailed protocolFor 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
Request a detailed protocolMBP 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
Request a detailed protocolFEN1 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
Request a detailed protocolCtIP 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
Request a detailed protocolThe 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
Request a detailed protocolFor 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)
Request a detailed protocolIP 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
Request a detailed protocolHeLa 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
Request a detailed protocol5 µ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)
Request a detailed protocol4X107 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
Request a detailed protocolImmunofluorescence 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
Request a detailed protocolEvery 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.
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
-
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 helicaseThe Journal of Biological Chemistry 289:29975–29993.https://doi.org/10.1074/jbc.M114.567073
-
Mitochondrial genome instability in human cancersMutation Research 488:9–23.https://doi.org/10.1016/s1383-5742(00)00063-6
-
Mitochondrial genetic variation is enriched in G-quadruplex regions that stall DNA synthesis in vitroHuman Molecular Genetics 29:1292–1309.https://doi.org/10.1093/hmg/ddaa043
-
Chromatin immunoprecipitation (CHIP)Cold Spring Harbor Protocols 2009:db.https://doi.org/10.1101/pdb.prot5279
-
Non-B DB v2.0: a database of predicted non-B DNA-forming motifs and its associated toolsNucleic Acids Research 41:D94–D100.https://doi.org/10.1093/nar/gks955
-
High-throughput sequencing of DNA G-quadruplex structures in the human genomeNature Biotechnology 33:877–881.https://doi.org/10.1038/nbt.3295
-
The organization and inheritance of the mitochondrial genomeNature Reviews. Genetics 6:815–825.https://doi.org/10.1038/nrg1708
-
Time-Dependent predominance of nonhomologous DNA end-joining pathways during embryonic development in miceJournal of Molecular Biology 417:197–211.https://doi.org/10.1016/j.jmb.2012.01.029
-
G-Quadruplex formation by single-base mutation or deletion of mitochondrial DNA sequencesBiochimica et Biophysica Acta. General Subjects 1863:418–425.https://doi.org/10.1016/j.bbagen.2018.11.009
-
Transcription of the mammalian mitochondrial genomeAnnual Review of Biochemistry 53:573–594.https://doi.org/10.1146/annurev.bi.53.070184.003041
-
Recognition of (dG) N. (dC) N sequences by endonuclease G. characterization of the calf thymus nucleaseThe Journal of Biological Chemistry 264:3301–3310.
-
Homologous recombination-mediated repair of DNA double-strand breaks operates in mammalian mitochondriaCellular and Molecular Life Sciences 75:1641–1655.https://doi.org/10.1007/s00018-017-2702-y
-
Mitochondrial DNA deletions are associated with non-B DNA conformationsNucleic Acids Research 40:7606–7621.https://doi.org/10.1093/nar/gks500
-
Dna instability at chromosomal fragile sites in cancerCurrent Genomics 11:326–337.https://doi.org/10.2174/138920210791616699
-
High base pair opening rates in tracts of GC base pairsThe Journal of Biological Chemistry 274:6957–6962.https://doi.org/10.1074/jbc.274.11.6957
-
Potential roles for G-quadruplexes in mitochondriaCurrent Medicinal Chemistry 26:2918–2932.https://doi.org/10.2174/0929867325666180228165527
-
G4 resolvase 1 tightly binds and unwinds unimolecular G4-DNANucleic Acids Research 39:7161–7178.https://doi.org/10.1093/nar/gkr234
-
Functional analysis of human FEN1 in Saccharomyces cerevisiae and its role in genome stabilityHuman Molecular Genetics 8:2263–2273.https://doi.org/10.1093/hmg/8.12.2263
-
Age-Associated change in mitochondrial DNA damageFree Radical Research 29:573–579.https://doi.org/10.1080/10715769800300611
-
The mitochondrial DNA 9-bp deletion polymorphism is a risk factor for hepatocellular carcinoma in the Chinese populationGenetic Testing and Molecular Biomarkers 16:330–334.https://doi.org/10.1089/gtmb.2011.0208
-
Mitochondrial insertion-deletion polymorphism: role in disease pathologyGenetic Testing and Molecular Biomarkers 15:361–364.https://doi.org/10.1089/gtmb.2010.0205
-
Mitochondrial DNA and genetic diseaseEssays in Biochemistry 47:139–151.https://doi.org/10.1042/bse0470139
-
Anti-Apoptotic protein BCL2 down-regulates DNA end joining in cancer cellsThe Journal of Biological Chemistry 285:32657–32670.https://doi.org/10.1074/jbc.M110.140350
-
Differential accumulations of 4,977 bp deletion in mitochondrial DNA of various tissues in human ageingBiochimica et Biophysica Acta 1226:37–43.https://doi.org/10.1016/0925-4439(94)90056-6
-
Removal of oxidative DNA damage via FEN1-dependent long-patch base excision repair in human cell mitochondriaMolecular and Cellular Biology 28:4975–4987.https://doi.org/10.1128/MCB.00457-08
-
Transient mammalian cell transfection with polyethylenimine (PEI)Methods in Enzymology 529:227–240.https://doi.org/10.1016/B978-0-12-418687-3.00018-5
-
Processing of G4 DNA by Dna2 helicase/nuclease and replication protein A (RPA) provides insights into the mechanism of Dna2/RPA substrate recognitionThe Journal of Biological Chemistry 283:24359–24373.https://doi.org/10.1074/jbc.M802244200
-
Role of endonuclease G in exogenous DNA stability in HeLa cellsBiochemistry. Biokhimiia 81:163–175.https://doi.org/10.1134/S0006297916020103
-
Mitochondria and human cancerCurrent Molecular Medicine 7:121–131.https://doi.org/10.2174/156652407779940495
-
Mechanism of fragility at BCL2 gene minor breakpoint cluster region during T (14; 18) chromosomal translocationThe Journal of Biological Chemistry 287:8688–8701.https://doi.org/10.1074/jbc.M111.307363
-
Mammalian mitochondrial endonuclease G. digestion of R-loops and localization in intermembrane spaceEuropean Journal of Biochemistry 269:5765–5770.https://doi.org/10.1046/j.1432-1033.2002.03238.x
-
Mitochondrial DNA and human evolutionAnnual Review of Genomics and Human Genetics 6:165–183.https://doi.org/10.1146/annurev.genom.6.080604.162249
-
Mitochondrial DNA in human malignancyMutation Research 488:119–133.https://doi.org/10.1016/s1383-5742(01)00053-9
-
A guide to computational methods for G-quadruplex predictionNucleic Acids Research 48:1–15.https://doi.org/10.1093/nar/gkz1097
-
Evidence for a triplex DNA conformation at the bcl-2 major breakpoint region of the T (14; 18) translocationJournal of Biological Chemistry 280:22749–22760.https://doi.org/10.1074/jbc.M502952200
-
Analysis of non-B DNA structure at chromosomal sites in the mammalian genomeMethods in Enzymology 409:301–316.https://doi.org/10.1016/S0076-6879(05)09017-8
-
The mitochondrial gene tree comes of ageAmerican Journal of Human Genetics 68:1315–1320.https://doi.org/10.1086/320615
-
BookMolecular Cloning A Laboratory Manual (Second edn)New York: Cold Spring Harbor Laboratory Press.
-
Two direct repeats cause most human mtDNA deletionsTrends in Genetics 20:393–398.https://doi.org/10.1016/j.tig.2004.07.003
-
Rna processing in human mitochondriaCell Cycle 10:2904–2916.https://doi.org/10.4161/cc.10.17.17060
-
Efficiency of nonhomologous DNA end joining varies among somatic tissues, despite similarity in mechanismCellular and Molecular Life Sciences 68:661–676.https://doi.org/10.1007/s00018-010-0472-x
-
Role of Magmas in protein transport and human mitochondria biogenesisHuman Molecular Genetics 19:1248–1262.https://doi.org/10.1093/hmg/ddq002
-
Mtdna control-region sequence variation suggests multiple independent origins of an “ asian-specific ” 9-bp deletion in sub-Saharan AfricansAmerican Journal of Human Genetics 58:595–608.
-
Characterization of rat Tom70 as a receptor of the preprotein translocase of the mitochondrial outer membraneJournal of Cell Science 115:1895–1905.https://doi.org/10.1242/jcs.115.9.1895
-
Mitochondrial DNA mutations in human diseaseNature Reviews. Genetics 6:389–402.https://doi.org/10.1038/nrg1606
-
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 repeatsPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 353:955–965.https://doi.org/10.1098/rstb.1998.0260
-
Endonuclease G: a mitochondrial protein released in apoptosis and involved in caspase-independent DNA degradationCell Death and Differentiation 8:1136–1142.https://doi.org/10.1038/sj.cdd.4400944
-
Discovery of the role of non-B DNA structures in mutagenesis and human genomic disordersThe Journal of Biological Chemistry 284:8997–9009.https://doi.org/10.1074/jbc.X800010200
-
DNA polymerase θ increases mutational rates in mitochondrial DNAACS Chemical Biology 13:900–908.https://doi.org/10.1021/acschembio.8b00072
-
A specific 4977-bp deletion of mitochondrial DNA in human ageing skinArchives of Dermatological Research 286:386–390.https://doi.org/10.1007/BF00371798
-
Structural insights into DNA degradation by human mitochondrial nuclease MGME1Nucleic Acids Research 46:11075–11088.https://doi.org/10.1093/nar/gky855
-
Ageing-Associated 5 kb deletion in human liver mitochondrial DNABiochemical and Biophysical Research Communications 178:124–131.https://doi.org/10.1016/0006-291x(91)91788-e
-
Nuclear localization and cell cycle-specific expression of CtIP, a protein that associates with the BRCA1 tumor suppressorThe Journal of Biological Chemistry 275:18541–18549.https://doi.org/10.1074/jbc.M909494199
Article and author information
Author details
Funding
Council of Scientific and Industrial Research, India (37(1692)/17/EMR-11)
- Sathees C Raghavan
Department of Atomic Energy, Government of India (21/01/2016-BRNS/35074)
- Sathees C Raghavan
Department of Biotechnology, Ministry of Science and Technology, India (BT/PR/3458/COE/34/33/2015)
- Sathees C Raghavan
IISc-DBT partnership programme (BT/PR27952-INF/22/212/2018)
- Sathees C Raghavan
Indian Institute of Science
- Sumedha Dahal
- Humaira Siddiqua
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
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.
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).
Copyright
© 2022, Dahal et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,360
- views
-
- 200
- downloads
-
- 13
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cell Biology
- Physics of Living Systems
The endoplasmic reticulum (ER), the largest cellular compartment, harbours the machinery for the biogenesis of secretory proteins and lipids, calcium storage/mobilisation, and detoxification. It is shaped as layered membranous sheets interconnected with a network of tubules extending throughout the cell. Understanding the influence of the ER morphology dynamics on molecular transport may offer clues to rationalising neuro-pathologies caused by ER morphogen mutations. It remains unclear, however, how the ER facilitates its intra-luminal mobility and homogenises its content. It has been recently proposed that intra-luminal transport may be enabled by active contractions of ER tubules. To surmount the barriers to empirical studies of the minuscule spatial and temporal scales relevant to ER nanofluidics, here we exploit the principles of viscous fluid dynamics to generate a theoretical physical model emulating in silico the content motion in actively contracting nanoscopic tubular networks. The computational model reveals the luminal particle speeds, and their impact in facilitating active transport, of the active contractile behaviour of the different ER components along various time–space parameters. The results of the model indicate that reproducing transport with velocities similar to those reported experimentally in single-particle tracking would require unrealistically high values of tubule contraction site length and rate. Considering further nanofluidic scenarios, we show that width contractions of the ER’s flat domains (perinuclear sheets) generate local flows with only a short-range effect on luminal transport. Only contractions of peripheral sheets can reproduce experimental measurements, provided they are able to contract fast enough.
-
- Cell Biology
Cell migration towards stiff substrates has been coined as durotaxis and implicated in development, wound healing, and cancer, where complex interplays between immune and non-immune cells are present. Compared to the emerging mechanisms underlying the strongly adhesive mesenchymal durotaxis, little is known about whether immune cells - migrating in amoeboid mode - could follow mechanical cues. Here, we develop an imaging-based confined migration device with a stiffness gradient. By tracking live cell trajectory and analyzing the directionality of T cells and neutrophils, we observe that amoeboid cells can durotax. We further delineate the underlying mechanism to involve non-muscle myosin IIA (NMIIA) polarization towards the soft-matrix-side but may not require differential actin flow up- or down-stiffness gradient. Using the protista Dictyostelium, we demonstrate the evolutionary conservation of amoeboid durotaxis. Finally, these experimental phenomena are theoretically captured by an active gel model capable of mechanosensing. Collectively, these results may shed new lights on immune surveillance and recently identified confined migration of cancer cells, within the mechanically inhomogeneous tumor microenvironment or the inflamed fibrotic tissues.