Genome-Wide Analysis Reveals Coating of the Mitochondrial Genome by TFAM

Abstract

Mitochondria contain a 16.6 kb circular genome encoding 13 proteins as well as mitochondrial tRNAs and rRNAs. Copies of the genome are organized into nucleoids containing both DNA and proteins, including the machinery required for mtDNA replication and transcription. The transcription factor TFAM is critical for initiation of transcription and replication of the genome, and is also thought to perform a packaging function. Although specific binding sites required for initiation of transcription have been identified in the D-loop, little is known about the characteristics of TFAM binding in its nonspecific packaging state. In addition, it is unclear whether TFAM also plays a role in the regulation of nuclear gene expression. Here we investigate these questions by using ChIP-seq to directly localize TFAM binding to DNA in human cells. Our results demonstrate that TFAM uniformly coats the whole mitochondrial genome, with no evidence of robust TFAM binding to the nuclear genome. Our study represents the first high-resolution assessment of TFAM binding on a genome-wide scale in human cells.

About the Authors:

Yun E. Wang

Affiliation: Division of Biology, California Institute of Technology, Pasadena, California, United States of America

Georgi K. Marinov

Affiliation: Division of Biology, California Institute of Technology, Pasadena, California, United States of America

Barbara J. Wold

Affiliation: Division of Biology, California Institute of Technology, Pasadena, California, United States of America

David C. Chan

* E-mail: [email protected]

Affiliations Division of Biology, California Institute of Technology, Pasadena, California, United States of America, Howard Hughes Medical Institute, California Institute of Technology, Pasadena, California, United States of America

Introduction

Mitochondria are essential eukaryotic organelles, serving as the epicenter of ATP production in the cell through oxidative phosphorylation. To perform this bioenergetic function, mitochondria utilize gene products encoded by the mitochondrial genome, a circular DNA that is 16.6 kb long. This genome is organized into DNA/protein structures termed nucleoids [1]. Mitochondrial DNA (mtDNA) encodes thirteen components of the electron transport chain, as well as 22 tRNAs and two ribosomal RNA genes. These gene products are essential for the proper function of the respiratory chain, and therefore maintenance of mtDNA levels and sequence fidelity is essential for cellular bioenergetics. In a human cell, there are hundreds to thousands of copies of the mtDNA genome [2,3]. Damage or depletion of mtDNA causes numerous inherited disorders, including Alpers’ Disease, ataxia neuropathy spectrum, and progressive external ophthalmoplegia [4,5]. Furthermore, loss and damage to mtDNA has been implicated in cardiovascular disease [6–9], diabetes [10–12], neurodegenerative disorders such as Alzheimer’s [13,14], and aging [15,16]. Strikingly, increasing mtDNA copy number promotes cell survival or function in many models of disease associated with decreased mtDNA abundance, such as diabetes [12,17], aging [18], Alzheimer’s [19], and Parkinson’s [20,21]. Thus, it is critical to understand how mtDNA copy number and integrity are maintained.

Mitochondrial transcription factor A (TFAM) is a DNA binding protein that plays multiple roles in regulating mtDNA function. As a sequence-specific transcription factor, it binds upstream of the light strand promoter (LSP) and heavy strand promoter 1 (HSP1) to activate initiation of transcription. At these sites, the footprint of TFAM binding is ~22 bp long [22,23]. As a result, TFAM is essential for production of gene products from the mitochondrial genome. In addition, TFAM is required for normal mtDNA copy number, because RNA primers generated from LSP are used to prime mtDNA replication [24,25]. Mice heterozygous for a knockout of TFAM exhibit not only an expected reduction (22%) in mitochondrial transcript levels in the heart and kidney, but also a universal 34% reduction in mtDNA copy number across all assayed tissues. Furthermore, homozygous knockout mice have no detectable levels of mtDNA and die during embryogenesis [26], highlighting the importance of TFAM in maintenance of mtDNA levels and in cellular and organismal viability.

Apart from its sequence-specific functions, TFAM is thought to organize the mtDNA genome by coating it in a nonspecific manner. Although how TFAM packages mtDNA is not well-understood, it is known to bind nonspecifically to DNA [27] and is estimated to be sufficiently abundant to coat the genome completely [28–30]. One model suggests that nonspecific binding radiates from the TFAM LSP binding site, which acts as a nucleation site for subsequent cooperative binding in a phased pattern to yield an inter-genome homogeneous pattern of binding [31,32]. The packaging function of TFAM appears to have important consequences for maintenance of the mtDNA genome. A TFAM variant that is deficient in transcriptional activation but competent in DNA binding is capable of preventing mtDNA depletion [33]. Therefore, as a prominent component of mtDNA nucleoids, TFAM appears to coat the mitochondrial genome, perhaps protecting it from turnover or deleterious damage.

Despite the importance of the associations of TFAM with mtDNA in the maintenance of mtDNA integrity and in cellular viability, these interactions have only been visualized in vivo at low resolution [34]. Therefore, to capture a high-resolution profile of TFAM-mtDNA interactions across the entire mitochondrial genome, we performed chromatin immunoprecipitation followed by massively parallel sequencing (ChIP-seq) for TFAM in human HeLa cells.

Results

Detection of TFAM-DNA interactions using ChIP-seq

To characterize TFAM binding to both the mitochondrial and nuclear genomes in an unbiased manner, we performed ChIP-seq targeting TFAM in HeLa cells. Because ChIP-seq data is highly dependent on the use of high-quality antibodies, we generated two new TFAM monoclonal antibodies (20G2C12 and 20F8A9) that efficiently immunoprecipitated TFAM (Figure 1A). Both of these antibodies gave clean mitochondrial and nucleoid signals in immunofluorescence experiments with cultured HeLa cells (Figure 1C,D). The 20G2C12 antibody also performed well in Western blots of whole-cell lysates, recognizing a single protein band of ~23 kDa (Figure 1B).

[Figure omitted. See PDF.]

Figure 1. Characterization of TFAM monoclonal antibodies.

(A) Immunoprecipitation of TFAM from cell lysates. HeLa cell lysate was applied to sheep anti-mouse Dynabeads conjugated to anti-Myc, 20G2C12 TFAM antibody, 20F8A9 TFAM antibody, or a 50/50 mixture of 20G2C12 and 20F8A9 TFAM antibodies. The labeled bands are: 1) Antibody heavy chain; 2) antibody light chain; 3) TFAM. (B) Western blot using the 20G2C12 antibody detects a ~23kDa band. (C and D) Immunocytochemistry showing TFAM localization. Mitochondria were identified by PPIF staining; mtDNA was identified by anti-DNA staining. There was no evidence for nuclear localization of TFAM using either antibody.

https://doi.org/10.1371/journal.pone.0074513.g001

Given the high efficiency of 20G2C12 in immunoprecipitating TFAM, as well as its high specificity, we used it to capture TFAM-associated DNA fragments for ChIP-seq analysis. DNA was sonicated prior to immunoenrichment and size-selected prior to library building so that the average fragment length of the final library was centered around 200 bp, a fragment distribution allowing for high-resolution deconvolution of binding events. We generated 3 replicates and matching controls. The sequencing depth of all samples was between 18 million and 48 million mappable reads, which is generally sufficient for comprehensive identification of transcription factor binding sites [35].

A common concern with ChIP-seq datasets is the variability of enrichment for true binding events as compared to background. In a typical ChIP-seq experiment, a minority of sequencing reads originates from binding events, with the majority representing random genomic DNA. Even for the same DNA binding factor, large variations in the strength of enrichment can be observed, and therefore it is critical to assess the degree of enrichment before downstream analysis. A number of ChIP-seq quality control metrics have been developed [35] for nuclear transcription factors. However, TFAM is expected to bind to the mitochondrial genome, which has very different characteristics from the nuclear genome. In addition, it is predicted to bind both in the classical localized manner [36] as well as broadly across the mitochondrial genome. As a result, metrics for evaluating nuclear transcription factors are not well-suited for analysis of TFAM binding data. We therefore examined the fraction of sequencing reads in our libraries mapping to the mitochondria as a proxy for the enrichment of TFAM binding events. Strikingly, between 30% and 75% of TFAM ChIP-seq reads mapped to the mitochondrial genome, while less than 2% of reads mapped to the mitochondrial genome in the input samples, indicating that our TFAM ChIP-seq datasets are indeed highly enriched for TFAM binding events (Figure 1B). We note that 75% ChIP enrichment is extremely high (in fact, practically unprecedented) for any transcription factor dataset [35], thus underscoring the high experimental quality of our datasets.

Because partial copies of the mitochondrial genome are also present in the nuclear genome, not all reads originating from mtDNA can be mapped uniquely. Therefore, we characterized TFAM binding to mtDNA and to the nuclear genome separately. We analyzed mitochondrial binding events by aligning sequencing reads to the mitochondrial genome alone (restricting our analysis to reads mapping perfectly without any mismatches to further increase mapping accuracy), and analyzed binding to the nuclear genome by aligning only the reads which did not map to the mitochondrial genome, as outlined in Figure 2A. For a standard nuclear transcription factor, this approach may cause some reads originating from the nuclear genome to artificially map to the mitochondrial genome. However, given that TFAM is known to bind to the mitochondrial genome and the extremely high enrichment for TFAM binding to mtDNA in our TFAM ChIP-seq libraries, this should not be a significant confounding factor.

[Figure omitted. See PDF.]

Figure 2. ChIP-seq analysis of genome-wide TFAM binding.

(A) Overview of computational processing of data. Reads were trimmed to 36 bp and then either mapped against the mitochondrial genome (ChrM), or the complete hg19 version of the genome. After removing multireads and alignments to the mitochondrial genome, peaks in the nuclear genome were called using MACS2. (B) The proportion of sequencing reads mapping to chrM in ChIP and input datasets. All replicates of the ChIP-seq resulted in at least 30% of reads mapping to the mitochondrial genome, much greater than the 0.4-1.9% of reads mapping to mtDNA in the input datasets. Replicates 1-3 were performed using the 20G2C12 antibody, while Replicate 4 was performed using the 20F8A9 antibody.

https://doi.org/10.1371/journal.pone.0074513.g002

TFAM coats the mitochondrial genome

As discussed above, TFAM has not only been proposed to bind specifically to well-defined binding sites in the D-loop, but has also been suggested to play a nonspecific packaging role in the nucleoid that is essential for mtDNA integrity. However, little is known about the pattern of non-specific binding of TFAM to the mitochondrial genome. Localized binding at the D-loop and diffuse binding across the rest of the genome are expected to result in distinct ChIP-seq signal profiles. Localized, “point-source” binding to DNA results in an asymmetric distribution of reads mapping to the forward and reverse strand around the binding site of the protein [36,37], while diffuse binding does not produce such strand asymmetry.

To characterize TFAM binding to mtDNA, we examined the forward and reverse strand read distribution after mapping TFAM ChIP-seq and input library reads to the mitochondrial genome. Strikingly, we did not observe regions of obvious enrichment and strand asymmetry in the D-loop; in particular, we did not see specific binding at the predicted HSP1 and LSP sites. On the whole, the TFAM ChIP-seq signal was broadly distributed over the whole mitochondrial chromosome, and while coverage was not perfectly uniform, the amplitude of the non-uniformity was not significant, and the signal profile closely tracked that of the input sample (Figure 3). The low level of non-uniformity likely results from sequencing biases, which has been documented to skew coverage [38,39]. Because our libraries were carefully size-selected for fragments in the 200 bp range, discrete TFAM binding sites would be expected to yield discrete signal localizations. Therefore, we interpret these results as evidence for the uniform coating of the whole mitochondrial genome by TFAM. We observed one region of apparent localized enrichment exhibiting strand asymmetry in the ND2 ORF near the origin of light strand replication (OL) (Figure 3F), which we discuss in the Discussion section.

[Figure omitted. See PDF.]

Figure 3. Coating of the mitochondrial genome by TFAM in HeLa cells.

Circos plot of plus strand and minus strand TFAM ChIP-seq and input read density signal over chrM. (A, E) Annotation of protein coding (green on forward/heavy strand, red on reverse/light strand), ribosomal RNA (blue) and tRNA (blue on forward/heavy strand, grey on reverse/light strand) transcripts. (B) D-loop (black), LSP promoter (large red tile), known LSP TFAM binding site (small red tile), HSP promoter (large blue tile), known HSP1 TFAM binding site (small blue tile), and origins of heavy strand replication (Ori-b, orange tile; OH, yellow tile). (C) TFAM ChIP-seq signal on forward (red) and reverse (blue) strands. (D) Input signal on forward (red) and reverse (blue) strands. (F) Origin of light strand replication (yellow tile). Note that the input signal is exaggerated 60-fold relative to the ChIP-seq signal in order to visualize coverage irregularities. The signal from the TFAM ChIP-seq largely follows that of the input, indicating generalized binding across the mitochondrial genome.

https://doi.org/10.1371/journal.pone.0074513.g003

To further verify our results, we carried out ChIP-seq against TFAM with a second TFAM monoclonal antibody, 20F8A9. We obtained similar results (Figure S1) and found significant correlation between the 20F8A9 dataset and the three datasets obtained from the 20G2C12 antibody datasets (p < 0.0001).

No evidence for binding to the nuclear genome

Previous studies have suggested that TFAM can be found in the nucleus and that it modulates the transcription of nuclear genes. In rat neonatal cardiac myocytes, TFAM was found to bind to the promoter of SERCA2, the homolog of human sarco(endo) plasmic reticulum calcium-ATPase 2 (ATP2A2), and was implicated in regulating its transcription [40]. Given the extremely high degree of TFAM binding enrichment in our datasets, any robust nuclear TFAM binding events should be readily detectable. To analyze nuclear binding, we excluded all sequencing reads mapping to the mitochondrial genome and used the resulting set of reads to identify putative TFAM binding sites. We first looked for significant global read clustering using cross-correlation between reads mapping to the forward and the reverse DNA strands [35,36]. Cross-correlation plots for input samples and for TFAM ChIP-seq datasets were indistinguishable from each other (Figure 4A,B). Next, we called putative TFAM binding sites using MACS2 [41]. Using default settings (corresponding to a q-value cut-off of 10-2), we identified 72, 137 and 153 sites respectively for the three replicates generated with antibody 20G2C12, and a single site for the 20F8A9 antibody. However, manual inspection of each of the identified sites revealed that all were likely to represent artifacts, mostly associated with repetitive DNA sequences, as none had the expected strand asymmetry of read distribution around a binding site. Instead, the two strand profiles at each site were identical (summarized in Figure 4D, with the classic nuclear transcription factor NRSF shown for comparison in Figure 4C), and numerous unmappable regions and repetitive elements were present in the immediate vicinity of many of the called sites. Inspection of the ATP2A2 gene revealed no TFAM enrichment neither in the promoter region nor anywhere else in the neighborhood of the gene (Figure 4E). Furthermore, we do not detect nuclear localization of TFAM in our cells (Figure 1C). Therefore, in HeLa cells under normal growth conditions, we find no evidence for specific binding of TFAM to nuclear target genes.

[Figure omitted. See PDF.]

Figure 4. Absence of TFAM binding to the nuclear genome.

(A) Cross-correlation plot of input DNA computed over the nuclear genome. (B) Cross-correlation plot of TFAM ChIP-seq computed over the nuclear genome. (C) Distribution of ChIP-seq reads mapping to the plus and minus strand around called binding sites in a ChIP-seq dataset for the NRSF transcription factor [51] in HeLa cells, generated by the ENCODE consortium [52]. (D) Distribution of TFAM ChIP-seq reads mapping to the plus and minus strand around called binding sites indicates lack of real binding sites. (E) No ChIP-seq enrichment around the promoter of the SERCA2/ATP2A2 gene, previously suggested to be a TFAM target.

https://doi.org/10.1371/journal.pone.0074513.g004

Discussion

Previous in vitro studies have suggested that TFAM binds specifically to LSP and HSP1, and that it may also bind nonspecifically in a phased manner. Furthermore, evidence has been presented for its nuclear localization and action as a canonical nuclear transcription factor in rat neonatal cardiac myocytes. However, no direct genome-wide measurements of TFAM binding have been previously reported. Our TFAM ChIP-seq data reveal very high enrichment for reads mapping to the mitochondrial genome, but a binding pattern that largely mirrors the read distribution observed in the input DNA, suggesting broad, non-specific binding to mitochondrial genome. This pattern is highly reproducible, indicating that the average population-wide state of TFAM-mtDNA interactions is stable. We found no correlation between irregularities in TFAM signal distribution and characteristics of the mitochondrial genome such as GC content (data not shown). Thus, we conclude that TFAM binds to the mitochondrial genome nonspecifically and without bias when cells are grown under typical culture conditions. Although we do not observe the synchronized phased binding seen in in vitro studies, we cannot rule out a model where individual mtDNAs have such a pattern of binding initiating from a non-universal nucleation site.

Strikingly, we did not observe localized enrichment of binding at the known LSP and HSP1 TFAM binding sites. Peak patterns mirrored that of the input in these regions, and no ChIP-seq peaks displaying the canonical strand asymmetry in read distribution were observed. This finding can be explained by a model in which the interaction of TFAM with the LSP and HSP1 binding sites is relatively transient and infrequent compared to a more stable non-specific association with the genome in its packaging state.

We did detect one site in the genome exhibiting the characteristics of a specific, localized ChIP-seq peak, centered at 5175 bp in the ND2 ORF. The localized nature of the ChIP signal at this site suggests higher occupancy of TFAM. This peak localizes to 546 bp upstream of the OL. Strikingly, TFAM has previously been localized 520 bp upstream of the OL of rat mtDNA [42–44]. We found no sequence similarity between the rat and human sites, and in general this region of the mtDNA genome shows low homology between the two species. Further work will be required to understand the significance of this putative TFAM binding site.

Finally, analysis of all datasets for TFAM binding to the nuclear genome yielded no hits distinguishable from common ChIP-seq artifacts. Although Watanabe et al. observed regulation of the SERCA2 gene in rat myocytes, we did not detect TFAM binding at the promoter of its ortholog in humans. Previous studies have shown nuclear localization of TFAM in rat hepatoma cells [45], as well as an alternate isoform of TFAM in mouse testis nuclei [46]. We have thus far been unable to detect nuclear TFAM localization in HeLa cells (Figure 1C), suggesting that nuclear localization and transcriptional regulation may be cell type or perhaps species-dependent. ChIP-seq in different cell lines may be able to detect such nuclear interactions.

We demonstrate here the first high-resolution ChIP-seq analysis of TFAM binding to the mitochondrial genome. Aside from generalized, largely non-specific binding across the mitochondrial genome, we detected a putative specific binding site upstream of the origin of light strand replication. We do not observe the expected binding at the known HSP1 and LSP sites, nor did we identify any nuclear binding sites. An area that remains to be explored is the dynamic nature of TFAM-DNA interactions with respect to both the nuclear and mitochondrial genomes. ChIP-chip on the yeast mitochondrial genome has shown that metabolic changes can lead to differential binding of the yeast TFAM homolog, Abf2p [47]. It is possible that such remodeling also occurs in the mammalian system, and further studies will provide insight into the dynamic nature of the mtDNA-protein interactions within the nucleoid that serve to protect its integrity.

Materials and Methods

Cell growth and treatment

HeLaS3 cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM, Invitrogen #11995) containing 10% bovine serum (Invitrogen #16170), penicillin and streptomycin, and additional L-glutamine (2mM). Cells were fed 24 hours before harvest for ChIP-seq, which was performed at 80-90% confluency.

Antibody Production and characterization

Antibodies were produced by the Caltech Monoclonal Antibody Facility and raised against the full-length TFAM protein in mouse. Immunoprecipitation with 20G2C12 and 20F8A9 TFAM antibodies and Myc antibody (Santa Cruz #sc-40) was performed according to established protocols using M-280 sheep anti-mouse Dynabeads (Invitrogen #11201D). Immunoblotting of IP products was performed using a monoclonal TFAM 18G102B2E11 antibody, also custom generated, at 1:2000, with goat anti-mouse HRP antibody (1:10,000, Jackson ImmunoResearch #115-056-003). Immunoblotting of HeLa whole cell lysate with 20G2C12 was performed at a 1:200 dilution and with goat anti-mouse HRP antibody.

Immunocytochemistry

HeLa cells cultured as described above were plated onto poly-lysine coated glass coverslips 48 hours prior to fixation in formaldehyde and permeabilization with 0.1% Triton X-100. For colocalization of TFAM to mitochondria, 20G2C12 or 20F8A9 antibodies were used at 1:10 in conjunction with PPIF at 1:200 (ProteinTech #18466-1-AP). Secondary antibodies were goat anti-mouse AF488 (1:500, Invitrogen #A11001) and donkey anti-rabbit AF546 (1:500, Invitrogen #A10040). Cells were also stained with DAPI to visualize nuclei. Immunocytochemistry to visualize colocalization of mitochondrial nucleoids and TFAM was performed sequentially due to both antibodies being raised in mouse. Sequential immunostaining yielded no background fluorescence due to cross-antibody reactivity (data not shown). Order was as follows: anti-TFAM antibody (1:10); goat anti-mouse AF488 (1:500, Invitrogen #A11001); anti-DNA antibody (1:25, Millipore #CBL186); goat anti-mouse AF555 (1:500, Invitrogen #A21426), DAPI. Images were acquired with a Zeiss LSM 710 confocal microscope with PlanApochromat 63X/1.4 oil objective. Z-stack acquisitions were converted to maximum z-projections using ImageJ software.

Chromatin immunoprecipitation and sequencing

ChIP experiments and preparation of DNA for sequencing were performed following standard procedures [48] with some modifications. Cells were fixed for 10min at RT in 1% formaldehyde, harvested using a cell scraper, washed once in ice-cold PBS, and resuspended in RIPA buffer with protease inhibitor. The sample was then sonicated using a 3.2mm microtip (QSonica Sonicator 4000) at 30s on/30s off intervals and 40% amplitude for 180min while in a -30°C 3:1 isopropanol and water bath containing dry ice. Subsequent steps were performed as per the standard protocol. DNA was size-selected during library building to an average fragment size of 200bp. Libraries were sequenced using Illumina GAIIx and Illumina HiSeq 2000. Sequencing data is available under GEO accession record GSE48176.

Sequencing data processing and analysis

Sequencing reads were trimmed down to 36 bp and then mapped against either the female set of human chromosomes (excluding the Y chromosome and all random chromosomes and haplotypes) or the mitochondrial genome alone, using the hg19 version of the human genome as a reference. Bow tie 0.12.7 [49] was used for aligning reads, not allowing for any mismatches between the reads and the reference. ChIP-seq peaks were called using MACS2 [41] with default settings except for the mfold parameter, which was lowered to (2,30). Circos plots were generated using Circos version 0.60 [50]. Additional data processing was carried out using custom-written python scripts. ENCODE data was downloaded from the UCSC browser (http://hgdownload-test.cse.ucsc.edu/goldenPath/hg19/encodeDCC/wgEncodeHaibTfbs) and its use here complies with its terms of usage. Pearson correlation coefficient, t-test, and p values were calculated using embedded and custom Microsoft Excel functions.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Comparison of profiles of TFAM binding to mitochondrial genome.

Circos plots of TFAM ChIP-seq experiments: (1) 20F8A9 antibody ChIP-Seq; (2) 20G2C12 replicate 1; (3) 20G2C12 replicate 2; (4) 20G2C12 replicate 3. Read profiles are very similar across replicates and antibodies.

https://doi.org/10.1371/journal.pone.0074513.s001

(TIF)

Acknowledgments

We thank Igor Antoshechkin and the Millard and Muriel Jacobs Genetics and Genomics Laboratory for assistance with library building and high-throughput sequencing. We thank Henry Amrhein and Diane Trout for computational assistance, Susan Ou for monoclonal antibody production, and Elizabeth Nelson for help with antibody characterization.

Author Contributions

Conceived and designed the experiments: YEW DCC. Performed the experiments: YEW. Analyzed the data: YEW GKM. Wrote the manuscript: YEW GKM BJW DCC.

Citation: Wang YE, Marinov GK, Wold BJ, Chan DC (2013) Genome-Wide Analysis Reveals Coating of the Mitochondrial Genome by TFAM. PLoS ONE 8(8): e74513. https://doi.org/10.1371/journal.pone.0074513

References

1. Bogenhagen DF, Rousseau D, Burke S (2008) The layered structure of human mitochondrial DNA nucleoids. J Biol Chem 283(6): 3665-3675. PubMed: 18063578.

2. Bogenhagen D, Clayton DA (1974) The number of mitochondrial deoxyribonucleic acid genomes in mouse L and human HeLa cells. Quantitative isolation of mitochondrial deoxyribonucleic acid. J Biol Chem 249(24): 7991-7995. PubMed: 4473454.

3. Satoh M, Kuroiwa T (1991) Organization of multiple nucleoids and DNA molecules in mitochondria of a human cell. Exp Cell Res 196(1): 137-140. doi:10.1016/0014-4827(91)90467-9. PubMed: 1715276.

4. Suomalainen A, Isohanni P (2010) Mitochondrial DNA depletion syndromes - many genes, common mechanisms. Neuromuscul Disord 20(7): 429-437. doi:10.1016/j.nmd.2010.03.017. PubMed: 20444604.

5. Stumpf JD, Saneto RP, Copeland RC (2013) Clinical and Molecular Features of POLG-Related Mitochondrial Disease. Cold Spring Harb Perspect Biol 4(5): a011395. PubMed: 23545419.

6. Sugiyama S, Hattori K, Hayakawa M, Ozawa T (1991) Quantitative analysis of age-associated accumulation of mitochondrial DNA with deletion in human age-associated accumulation of mitochondrial DNA with deletion in human hearts. Biochem Biophys Res Commun 180(2): 894-899. doi:10.1016/S0006-291X(05)81149-0. PubMed: 1953759.

7. Ide T, Tsutsui H, Hayashidani S, Kang D, Suematsu N et al. (2001) Mitochondrial DNA damage and dysfunction associated with oxidative stress in failing hearts after myocardial infarction. Circ Res 88(5): 529-535. doi:10.1161/01.RES.88.5.529. PubMed: 11249877.

8. Karamanlidis G, Nascimben L, Couper GS, Shekar PS, del Monte F et al. (2010) Defective DNA replication impairs mitochondrial biogenesis in human failing hearts. Circ Res 106(9): 1541-1548. doi:10.1161/CIRCRESAHA.109.212753. PubMed: 20339121.

9. Karamanlidis G, Bautista-Hernandez FV, Fynn-Thompson F, Del Nido P, Tian R (2011) Impaired mitochondrial biogenesis precedes heart failure in right ventricular hypertrophy in congenital heart disease. Circ Heart Fail 4(6): 707-713. doi:10.1161/CIRCHEARTFAILURE.111.961474. PubMed: 21840936.

10. Maassen JA, 'T Hart LM, Van Essen E, Heine RJ, Nijpels G et al. (2004) Mitochondrial diabetes: molecular mechanisms and clinical presentation. Diabetes 53 Suppl 1: S103-S109. doi:10.2337/diabetes.53.2007.S103. PubMed: 14749274.

11. Simmons RA, Suponitsky-Kroyter I, Selak MA (2005) Progressive accumulation of mitochondrial DNA mutations and decline in mitochondrial function lead to beta-cell failure. J Biol Chem 280(31): 28785-28791. doi:10.1074/jbc.M505695200. PubMed: 15946949.

12. Gauthier BR, Wiederkehr A, Baquié M, Dai C, Powers AC et al. (2009) PDX1 deficiency causes mitochondrial dysfunction and defective insulin secretion through TFAM suppression. Cell Metab 10(2): 110-118. doi:10.1016/j.cmet.2009.07.002. PubMed: 19656489.

13. Coskun PE, Beal MF, Wallace DC (2004) Alzheimer’s brains harbor somatic mtDNA control-region mutations that suppress mitochondrial transcription and replication. Proc Natl Acad Sci USA 101(29): 10726-10731. doi:10.1073/pnas.0403649101. PubMed: 15247418.

14. Coskun P, Wyrembak J, Schriner SE, Chen HW, Marciniack C et al. (2012) A mitochondrial etiology of Alzheimer and Parkinson disease. Biochim Biophys Acta 1820(5): 553-564. doi:10.1016/j.bbagen.2011.08.008. PubMed: 21871538.

15. Corral-Debrinski M, Shoffner JM, Lott MT, Wallace DC (1992) Association of mitochondrial DNA damage with aging and coronary atherosclerotic heart disease. Mutat Res 275(3-6): 169-180. doi:10.1016/0921-8734(92)90021-G. PubMed: 1383759.

16. Trifunovic A, Larsson NG (2008) Mitochondrial dysfunction as a cause of ageing. J Intern Med 263(2): 167-178. doi:10.1111/j.1365-2796.2007.01905.x. PubMed: 18226094.

17. Suarez J, Hu Y, Makino A, Fricovsky E, Wang H et al. (2008) Alterations in mitochondrial function and cytosolic calcium induced by hyperglycemia are restored by mitochondrial transcription factor A in cardiomyocytes. Am J Physiol Cell Physiol 295(6): 1561-1568. doi:10.1152/ajpcell.00076.2008. PubMed: 19060297.

18. Hayashi Y, Yoshida M, Yamato M, Ide T, Wu Z et al. (2008) Reverse of age-dependent memory impairment and mitochondrial DNA damage in microglia by an overexpression of human mitochondrial transcription factor a in mice. J Neurosci 28(34): 8624-8634. doi:10.1523/JNEUROSCI.1957-08.2008. PubMed: 18716221.

19. Xu S, Zhong M, Zhang L, Wang Y, Zhou Z et al. (2009) Overexpression of TFAM protects mitochondria against beta-amyloid-induced oxidative damage in SH-SY5Y cells. FEBS J 276(14): 3800-3809. doi:10.1111/j.1742-4658.2009.07094.x. PubMed: 19496804.

20. Keeney PM, Quigley CK, Dunham LD, Papageorge CM, Iyer S et al. (2009) Mitochondrial gene therapy augments mitochondrial physiology in a Parkinson's disease cell model. Hum Gene Ther 20(8): 897-907. doi:10.1089/hum.2009.023. PubMed: 19374590.

21. Piao Y, Kim HG, Oh MS, Pak YK (2012) Overexpression of TFAM, NRF-1 and myr-AKT protects the MPP(+)-induced mitochondrial dysfunctions in neuronal cells. Biochim Biophys Acta:1820(5): 577-585. doi:10.1016/j.bbagen.2011.08.007. PubMed: 21856379.

22. Fisher RP, Clayton DA (1988) Purification and characterization of human mitochondrial transcription factor 1. Mol Cell Biol 8(8): 3496-3509. PubMed: 3211148.

23. Ngo HB, Kaiser JT, Chan DC (2011) The mitochondrial transcription and packaging factor TFAM imposes a U-turn on mitochondrial DNA. Nat Struct Mol Biol 18(11): 1290-1296. doi:10.1038/nsmb.2159. PubMed: 22037171.

24. Chang DD, Clayton DA (1984) Precise identification of individual promoters for transcription of each strand of human mitochondrial DNA. Cell 36(3): 635-643. doi:10.1016/0092-8674(84)90343-X. PubMed: 6697390.

25. Chang DD, Clayton DA (1985) Priming of human mitochondrial DNA replication occurs at the light-strand promoter. Proc Natl Acad Sci USA 82(2): 351-355. doi:10.1073/pnas.82.2.351. PubMed: 2982153.

26. Larsson NG, Wang J, Wilhelmsson H, Oldfors A, Rustin P et al. (1998) Mitochondrial transcription factor A is necessary for mtDNA maintenance and embryogenesis in mice. Nat Genet 18(3): 231-236. doi:10.1038/ng0398-231. PubMed: 9500544.

27. Fisher RP, Parisi MA, Clayton DA (1989) Flexible recognition of rapidly evolving promoter sequences by mitochondrial transcription factor 1. Genes Dev 3(12b):2202-17.

28. Alam TI, Kanki T, Muta T, Ukaji K, Abe Y et al. (2003) Human mitochondrial DNA is packaged with TFAM. Nucleic Acids Res 31(6): 1640-1645. doi:10.1093/nar/gkg251. PubMed: 12626705.

29. Ekstrand MI, Falkenberg M, Rantanen A, Park CB, Gaspari M et al. (2004) Mitochondrial transcription factor A regulates mtDNA copy number in mammals. Hum Mol Genet 13(9): 935-944. doi:10.1093/hmg/ddh109. PubMed: 15016765.

30. Kaufman BA, Durisic N, Mativetsky JM, Costantino S, Hancock MA et al. (2007) The mitochondrial transcription factor TFAM coordinates the assembly of multiple DNA molecules into nucleoid-like structures. Mol Biol Cell 18(9): 3225-3236. doi:10.1091/mbc.E07-05-0404. PubMed: 17581862.

31. Fisher RP, Lisowsky T, Parisi MA, Clayton DA (1992) DNA wrapping and bending by a mitochondrial high mobility group-like transcriptional activator protein. J Biol Che 267(5): 3358-3367. PubMed: 1737790.

32. Ghivizzani SC, Madsen CS, Nelen MR, Ammini CV, Hauswirth WW (1994) In organello footprint analysis of human mitochondrial DNA: human mitochondrial transcription factor A interactions at the origin of replication. Mol Cell Biol 14(12): 7717-7730. PubMed: 7969115.

33. Kanki T, Ohgaki K, Gaspari M, Gustafsson CM, Fukuoh A et al. (2004) Architectural role of mitochondrial transcription factor A in maintenance of human mitochondrial DNA. Mol Cell Biol 24(22): 9823-9834. doi:10.1128/MCB.24.22.9823-9834.2004. PubMed: 15509786.

34. Ohgaki K, Kanki T, Fukuoh A, Kurisaki H, Aoki Y et al. (2007) The C-terminal tail of mitochondrial transcription factor a markedly strengthens its general binding to DNA. J Biochem 141(2): 201-211. PubMed: 17167045.

35. Landt SG, Marinov GK, Kundaje A, Kheradpour P, Pauli F et al. (2011) ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia. Genome Res 22(9): 1813-1831. PubMed: 22955991.

36. Kharchenko PV, Tolstorukov MY, Park PJ (2008) Design and analysis of ChIP-seq experiments for DNA-binding proteins. Nat Biotechnol 26: 1351-1359. doi:10.1038/nbt.1508. PubMed: 19029915.

37. Pepke S, Wold B, Mortazavi A (2009) Computation for ChIP-seq and RNA-seq studies. Nat Methods 6(11 Suppl): S22-S32. doi:10.1038/nmeth.1371. PubMed: 19844228.

38. Dohm JC, Lottaz C, Borodina T, Himmelbauer H (2008) Substantial biases in ultra-short read data sets from high-throughput DNA sequencing. Nucleic Acids Res 36(16): e105. doi:10.1093/nar/gkn425. PubMed: 18660515.

39. Ross MG, Russ C, Costello M, Hollinger A, Lennon NJ et al. (2013) Characterizing and measuring bias in sequence data. Genome Biol 14(5): R51. doi:10.1186/gb-2013-14-5-r51. PubMed: 23718773.

40. Watanabe A, Arai M, Koitabashi N, Niwano K, Ohyama Y et al. (2011) Mitochondrial transcription factors TFAM and TFB2M regulate Serca2 gene transcription. Cardiovasc Res 90(1): 57-67. doi:10.1093/cvr/cvq374. PubMed: 21113058.

41. Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS et al. (2008) Model-based analysis of ChIP-Seq (MACS). Genome Biol 9(9): R137. doi:10.1186/gb-2008-9-9-r137. PubMed: 18798982.

42. Gadaleta G, D’Elia D, Capaccio L, Saccone C, Pepe G (1996) Isolation of a 25-kDa protein binding to a curved DNA upstream the origin of the L strand replication in the rat mitochondrial genome. J Biol Chem 271(23): 13537-13541. doi:10.1074/jbc.271.23.13537. PubMed: 8662779.

43. Cingolani G, Capaccio L, D’Elia D, Gadaleta G (1997) In organelle footprinting analysis of rat mitochondrial DNA: protein interaction upstream of the Ori-L. Biochem Biophys Res Commun 231(3): 856-860. doi:10.1006/bbrc.1997.6203. PubMed: 9070910.

44. Pierro P, Capaccio L, Gadaleta G (1999) The 25 kDa protein recognizing the rat curved region upstream of the origin of the L-strand replication is the rat homologue of the human mitochondrial transcription factor A. FEBS Lett 457(3): 307-310. doi:10.1016/S0014-5793(99)01055-8. PubMed: 10471798.

45. Dong X, Ghoshal K, Majumder S, Yadav SP, Jacob ST (2002) Mitochondrial transcription factor A and its downstream targets are up-regulated in a rat hepatoma. J Biol Chem 277(45): 43309-43318. doi:10.1074/jbc.M206958200. PubMed: 12198131.

46. Larsson NG, Garman JD, Oldfors A, Barsh GS, Clayton DA (1996) A single mouse gene encodes the mitochondrial transcription factor A and a testis-specific nuclear HMG-box protein. Nat Genet 13(3): 296-302. doi:10.1038/ng0796-296. PubMed: 8673128.

47. Kucej M, Kucejova B, Subramanian R, Chen XJ, Butow RA (2008) Mitochondrial nucleoids undergo remodeling in response to metabolic cues. J Cell Sci 121(11): 1861-1868. doi:10.1242/jcs.028605. PubMed: 18477605.

48. Johnson DS, Mortazavi A, Myers RM, Wold B (2007) Genome-wide mapping of in vivo protein-DNA interactions. Science 316(5830): 1497-1502. doi:10.1126/science.1141319. PubMed: 17540862.

49. Langmead B, Trapnell C, Pop M, Salzberg SL (2009) Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol 10(3): R25. doi:10.1186/gb-2009-10-3-r25. PubMed: 19261174.

50. Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R et al. (2009) Circos: an information aesthetic for comparative genomics. Genome Res 19(9): 1639-1645. doi:10.1101/gr.092759.109. PubMed: 19541911.

51. Schoenherr CJ, Anderson DJ (1995) The neuron-restrictive silencer factor (NRSF): a coordinate repressor of multiple neuron-specific gene. Science 267(5202): 1360-1363. doi:10.1126/science.7871435. PubMed: 7871435.

52. ENCODE Project Consortium (2011) A user’s guide to the encyclopedia of DNA elements (ENCODE). PLOS Biol 9(4):e1001046. PubMed: 21526222.