AbstractMitochondrial DNA copy number is strictly regulated during development as naive cells differentiate into mature cells to ensure that specific cell types have sufficient copies of mitochondrial DNA to perform their specialised functions. Mitochondrial DNA haplotypes are defined as specific regions of mitochondrial DNA that cluster with other mitochondrial sequences to show the phylogenetic origins of maternal lineages. Mitochondrial DNA haplotypes are associated with a range of phenotypes and disease. To understand how mitochondrial DNA haplotypes induce these characteristics, we used four embryonic stem cell lines that have the same set of chromosomes but possess different mitochondrial DNA haplotypes. We show that mitochondrial DNA haplotypes influence changes in chromosomal gene expression and affinity for nuclear-encoded mitochondrial DNA replication factors to modulate mitochondrial DNA copy number, two events that act synchronously during differentiation. Global DNA methylation analysis showed that each haplotype induces distinct DNA methylation patterns, which, when modulated by DNA demethylation agents, resulted in skewed gene expression patterns that highlight the effectiveness of the new DNA methylation patterns established by each haplotype. The haplotypes differentially regulate α-ketoglutarate, a metabolite from the TCA cycle that modulates the TET family of proteins, which catalyse the transition from 5-methylcytosine, indicative of DNA methylation, to 5-hydroxymethylcytosine, indicative of DNA demethylation. Our outcomes show that mitochondrial DNA haplotypes differentially modulate chromosomal gene expression patterns of naive and differentiating cells by establishing mitochondrial DNA haplotype-specific DNA methylation patterns.