Twin pair analysis uncovers novel links between DNA methylation, mitochondrial DNA quantity and obesity
=======================================================================================================

* Aino Heikkinen
* Vivienne F C Esser
* Sara Lundgren
* Seung Hyuk T Lee
* Antti Hakkarainen
* Jesper Lundbom
* Juho Kuula
* Per-Henrik Groop
* Sini Heinonen
* Päivi Pajukanta
* Jaakko Kaprio
* Kirsi H Pietiläinen
* Shuai Li
* Miina Ollikainen

## Abstract

Alterations in mitochondrial metabolism in obesity may indicate disrupted communication between mitochondria and nucleus, crucial for adapting to changing metabolic demands. Epigenetic modifications, particularly DNA methylation, may influence this intricate interplay, though the specifics remain poorly understood. Leveraging data from the subcohort of the Finnish Twin Cohort (n=173; 86 full twin pairs) that includes comprehensive measurements of obesity-related outcomes, mitochondrial DNA quantity (mtDNAq) and nuclear DNA methylation levels in adipose and muscle tissue, we identified one locus at *SH3BP4* (cg19998400) significantly associated with mtDNAq in adipose tissue (FDR<0.05). *SH3BP4* methylation correlated with its gene expression. Additionally, 14 out of the 35 obesity-related traits displayed significant associations with both *SH3BP4* methylation and mtDNAq in adipose tissue. Using the method that infers causality from examination of familial confounding (ICE FALCON) our data suggests that mtDNAq, insulin sensitivity and certain body fat measures are causal to *SH3BP4* methylation. The examination of mtDNAq and obesity-related traits suggested causation from mtDNAq to obesity which could not, however, be distinguished from potential unmeasured within-individual confounding. In conclusion, our findings underscore the impact of mtDNAq on DNA methylation and expression of the *SH3BP4* gene within adipose tissue, with potential implications for obesity.

Keywords
*   DNA methylation
*   Mitochondrial DNA
*   Obesity
*   Monozygotic twins
*   ICE FALCON

## Introduction

The significance of mitochondria as a primary energy source for cell growth and survival is indisputable. While mitochondria contain their own circular 16.6 kb genome, most of the genes required for mitochondrial function are encoded in the nuclear DNA. Therefore, cells require synchronized communication between the nucleus and mitochondria, known as mitonuclear communication, to adapt to changing metabolic demands. This intricate interplay can be subject to regulation by epigenetic mechanisms1, including DNA methylation2, which may influence the activity of gene expression. Consequently, the importance of DNA methylation for mitochondrial function has been highlighted in the literature3.

The existing body of literature supports a bidirectional relationship between mitochondrial function and nuclear DNA methylation. On one hand, different mitochondrial haplotypes and variants can lead to differences in DNA methylation4–7 via numerous metabolic pathways, such as the methionine cycle and the production of methyl groups3. On the other hand, nuclear DNA methylation may impact mitochondrial metabolism by regulating mitochondrial-associated gene expression8,9. The term ‘mitochondrial metabolism’ is used here as an umbrella term for any measurable trait associated with mitochondrial metabolism such as mtDNA quantity (mtDNAq), copy number, biogenesis, dynamics, and OXPHOS (oxidative phosphorylation) activity.

The vital role of mitochondria in cellular energy metabolism has placed them in the center of interest in human traits and diseases with metabolic symptoms, including obesity. Obesity is characterized by excess body weight which can lead to a range of systemic metabolic disturbances, and it has been previously associated with compromised mitochondrial biogenesis and OXPHOS capacity10,11. Furthermore, obesity is a highly heterogeneous trait, and the precise molecular phenotypes contributing to the altered mitochondrial metabolism remain unclear. Additionally, the exact tissue-specific roles of mitochondria in obesity are not well documented, although it has been shown that adipose tissue metabolism may be more affected by the acquired weight than muscle tissue11.

The research on the role of DNA methylation underlying obesity-associated decline in mitochondrial metabolism remains scarce. A recent study showed that mtDNAq influences cardiovascular disease and mortality through changes in DNA methylation12. DNA methylation profiles of adipocyte progenitor cells of individuals with obesity have also been linked to mitochondrial metabolism13. Despite these insights, the causal pathways and comprehensive characterization of the orchestrated effects of DNA methylation and mitochondrial decline in obesity and obesity-related phenotypes remain elusive.

Here, we aimed to identify DNA methylation sites associated with differences in mtDNAq in two primary tissues affected by obesity, namely subcutaneous adipose tissue and skeletal muscle. Additionally, we explored whether these methylation sites are linked to a comprehensive range of obesity-related outcomes and other related phenotypes, including several anthropometric and body composition measures, clinical biomarkers, physical activity and biological aging14. Finally, we explored the potential causal relationships between mtDNAq, DNA methylation and obesity-related measures with a method called Inference about Causation from Examination of FAmilial CONfounding (ICE FALCON)15, utilizing the monozygotic (MZ) twin pairs of the study cohort.

## Results

A total of 173 individuals participated in the study, comprising complete 81 MZ twin pairs and 5 DZ twin pairs (Table 1). The age range of the study cohort spanned from 23 to 70 years old, with females accounting for 59% of the cohort (Table 1). The mean BMI and fat percentage were 29.2kg/m2 and 39.6%, respectively, indicating that the sample is predominantly with overweight. In addition, average fasting glucose levels (5.6 mmol/l) and

View this table:
[Table 1.](http://medrxiv.org/content/early/2024/04/03/2024.04.02.24304959/T1)

Table 1. Participant characteristics.

HOMA-IR measures (2.1 units) were slightly elevated which suggests a reduced insulin sensitivity. Twenty-two participants were diagnosed with type 2 diabetes. The main analysis strategy is described in Figure 1, while Supplementary Table 1 summarizes the epigenetic age acceleration estimates of the study participants. Furthermore, the BMI-discordant MZ twin pairs with adipose tissue data (n=71 pairs) used in the causal inference analysis are presented in Supplementary Table 2.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/04/03/2024.04.02.24304959/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2024/04/03/2024.04.02.24304959/F1)

Figure 1. Study flowchart.

### Differential methylation analysis

To investigate whether nuclear DNA methylation was associated with mtDNAq, we performed an epigenome-wide association analysis (EWAS). In adipose tissue (n=153 individuals), we identified one CpG site (*cg19998400*) that was inversely associated with mtDNAq (FDR=0.002) (Fig. 2a). The CpG is in the enhancer region of *SH3BP4* (SH3 domain binding protein 4) gene that codes for a protein involved in intracellular signaling pathways. Another CpG (*cg17468563*) located in the *DHRS3* (Dehydrogenase/Reductase 3) enhancer region showed marginal association with mtDNAq (FDR=0.078) (Fig. 2a). Muscle tissue CpG methylation (n=155 individuals) was not associated with mtDNAq (Fig. 2b). The association profiles of the two tissues did not show similar patterns, as none of the top 100 CpG were found to be in common between adipose and muscle samples. Additionally, there was no strong correlation between the regression effect sizes (Supplementary Fig. 1). The rest of the manuscript focuses explicitly on the findings in adipose tissue, given that we did not detect any statistically significant mtDNAq-linked DNA methylation sites in muscle.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/04/03/2024.04.02.24304959/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2024/04/03/2024.04.02.24304959/F2)

Figure 2. Volcano plot of the epigenome-wide association study on mtDNA quantity
in (a) adipose tissue (n=153 individuals) and (b) skeletal muscle (n=155 individuals). Red dots indicate CpGs with FDR<0.10.

### Causal inference between adipose tissue mitochondrial DNA quantity and identified CpG sites

To investigate the potential evidence for causation underlying the associations between mtDNAq and the two identified CpG sites in adipose tissue, we performed ICE FALCON analysis for the complete MZ twin pairs with data on adipose tissue DNA methylation (n=68 MZ pairs). We analyzed each CpG separately given that we cannot assume them to have identical causal pathways. For instance, the *cg17468563* has been reported to be under genetic control of multiple loci16 whereas no meQTLs have been identified for *cg19998400*. Our data was consistent with the hypothesis that mtDNAq is causally linked to cg*19998400* (Table 2). This was suggested when we set mtDNAq as the independent variable and observed a marginal cross-twin cross-trait association (Model 2: Bcotwin=-0.160, p =0.091) which was attenuated towards null when adjusting for the twin’s own mtDNAq. There was no evidence of change in the within-individual association. Conversely, when using mtDNAq as the dependent variable, we did not observe any cross-twin cross-trait associations in Models 2 or 3. However, the increased absolute value of the regression coefficient in Model 3 is in line with mtDNAq being causal to methylation at *cg19998400* (Table 2). The change was not statistically significant, which may reflect a low sample size and consequently reduced statistical power. The results from ICE FALCON on *cg17468563* were more ambiguous, with the regression coefficients pointing to either causation from mtDNAq to methylation at *cg17468563* or presence of within-individual confounding (Table 2).

View this table:
[Table 2.](http://medrxiv.org/content/early/2024/04/03/2024.04.02.24304959/T2)

Table 2. 
Results from ICE FALCON analysis between mtDNA quantity and the identified CpG sites (n = 68 monozygotic twin pairs) in adipose tissue. Regression models were adjusted for age, sex, smoking, BMI and methylation beadchip and row. P-values < 0.05 are bolded.

### Causal inference between adipose tissue DNA methylation and gene expression

As DNA methylation is a potential mechanism to regulate gene expression, we investigated whether the DNA methylation at *cg19998400* and *cg17468563* were associated with the expression of their respective genes, *SH3BP4* (n=80 individuals) and *DHRS3* (n=80 individuals) in adipose tissue. There was a positive correlation between the expression and methylation of *SH3BP* (rPearson=0.46, p<0.001) (Fig. 3a), where the methylation explained solely 18% of the variation (marginal R-squared) and with familial confounding 58% (conditional R-squared) of the variation in gene expression (Fig. 3c). For the *DHRS3* expression and methylation, we observed a negative correlation (rPearson=-0.66, p<0.001) (Fig. 3b), with a marginal R-squared of 42% and conditional R-squared 56% (Fig 3c).

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/04/03/2024.04.02.24304959/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2024/04/03/2024.04.02.24304959/F3)

Figure 3. Relationship between DNA methylation and gene expression at two genomic loci in adipose tissue.
(a) Correlation between methylation at *cg19998400* and *SH3BP4* expression (n=80 individuals).

(b) Correlation between methylation at *cg17468563* and *DHRS3* expression (n=80 individuals).

(c) Variation in *SH3BP4* and *DHRS3* expression explained by methylation at *cg19998400* and *cg17468563* (yellow bars), respectively, and combined with common familial factors (genetic or environmental) (blue bars) (n=80 individuals).

DNAme = DNA methylation; L2PM = log2-counts per million; r = Pearson correlation coefficient

It is important to acknowledge that changes in gene expression patterns can also directly influence DNA methylation levels. Considering this, we conducted an ICE FALCON analysis to explore the direction of causality between gene expression and DNA methylation. We observed a marginal cross-twin cross-trait association between *SH3BP4* expression and *cg19998400* methylation (Model 2: Bcotwin=-0.112, p =0.078) (Supplementary Table 3) that attenuated towards null when adjusting for twin’s own *cg19998400* levels, whereas the twin’s own regression coefficient remained approximately the same. This observation is consistent with *cg19998400* methylation influencing the *SH3BP4* expression. On the other hand, the possibility for the presence of within-individual confounding cannot be ruled out. In contrast, with *DHRS3* and *cg17468563*, we did not find any compelling evidence for either causation or familial confounding (Supplementary Table 3).

### Obesity-related variables associated with adipose tissue *SH3BP4* methylation and mitochondrial DNA quantity

Considering the well-established connection between mitochondrial dysfunction and obesity, we were interested to study whether the mtDNAq-associated CpG sites would be also associated with different obesity-related variables. We specifically focused on the DNA methylation of *cg19998400* at the *SH3BP4* locus (herein referred to as *SH3BP4* methylation), as our data suggested a potential causal relationship with mtDNAq, in contrast to *cg17468563* methylation at *DHRS3.* Among a set of 35 obesity-related traits, six showed significant association only with mtDNAq, four with *SH3BP4* methylation, and 14 with both mtDNAq and *SH3BP4* methylation (FDR<0.05) (Fig 4., Supplementary Tables 4-5).

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/04/03/2024.04.02.24304959/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2024/04/03/2024.04.02.24304959/F4)

Figure 4. Standardized beta coefficients and standard errors for the associations between obesity-related outcomes and methylation at *cg19998400* (X-axis) or mtDNA quantity (Y-axis) (n=42-142 individuals) in adipose tissue.
Yellow background indicates variables with FDR<0.05 associated with *cg19998400* methylation only, blue indicates variables associated with mtDNA quantity only, and green indicates variables associated with both.

Three distinct measures of epigenetic age acceleration (EAA), namely Horvath, Hannum and PhenoAge, displayed a negative association with mtDNAq but not with *SH3BP4* methylation. In addition to EAA, mtDNAq showed negative associations with fat-free mass and adipsin levels, and positive with sports activity. Obesity-related variables that associated exclusively with *SH3BP4* methylation included elevated blood triglyceride levels, hsCRP, fasting glucose levels and systolic blood pressure. The 14 shared associations between mtDNAq and *SH3BP4* methylation included parameters mainly related to body fat composition, insulin sensitivity and HDL cholesterol levels. Consistent with existing literature, higher mtDNAq were correlated with lower body fat, and higher insulin sensitivity, HDL cholesterol levels and adiponectin levels.

### Causal inference between obesity-related traits, and adipose tissue *SH3BP4* methylation and mitochondrial DNA quantity

To discern the potential causal relationship between the 14 obesity-related variables, and both *SH3BP4* methylation and mtDNAq, we employed an ICE FALCON analysis for the complete MZ twin pairs in the cohort (Supplementary Table 2). Our findings indicate that certain variables related to insulin resistance and ectopic fat may exert causal influence on *SH3BP4* methylation, as suggested by significant cross-twin cross-trait association in Model 2 (body fat percentage Bcotwin=0.170, p=0.074; intra-abdominal fat Bcotwin=0.235, p =0.031; HOMA-IR Bcotwin=0.278, p =0.002; Matsuda Bcotwin=-0.253, p=0.006; fasting insulin Bcotwin=0.235, p =0.004) that attenuated towards null after conditioning on twin’s own corresponding measures (Fig. 5b). The changes in regression coefficients were significant in fat percentage (p=0.047) and marginally significant in Matsuda (p=0.062). Reversing the regression and using *SH3BP4* methylation as predictor variable also suggested causality from these obesity-related traits to methylation, based on the behavior of the regression coefficients (i.e., the cross-twin cross-trait association increased substantially while the Bself remained relatively stable) (Fig. 5c). Results of other obesity-related traits, specifically those measuring body size, liver fat percentage, subcutaneous fat and HDL cholesterol, were less clear. While the behavior of co-twin’s regression coefficients in ICE FALCON may indicate that these variables are a consequence of *SH3BP4* methylation (Fig 5b-c), the opposing signs between Model 1 and Model 2 coefficients, when using *SH3BP4* methylation as predictor variable X, can also suggest the presence of within-individual confounding (Fig 5c). Similarly, the ICE FALCON analysis on the association between mtDNAq and obesity-related variables suggested either mtDNAq being causal to most of the obesity-related variables or being subject to unmeasured within-individual confounding (Fig 5d-e).

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/04/03/2024.04.02.24304959/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2024/04/03/2024.04.02.24304959/F5)

Figure 5. Results from the ICE FALCON analysis between adipose tissue mtDNA quantity, *SH3BP4* methylation and obesity-related outcomes.
(a) Behavior of the ICE FALCON regression coefficients of Models 1-3 in each causal scenario.

(b-e) Standardized regression coefficients and 95% confidence intervals analysis for the ICE FALCON regression coefficients (n=21-71 pairs) of (b) obesity-related outcomes regressed against *SH3BP4* methylation, (c) *SH3BP4* methylation regressed against obesity-related outcomes, (d) obesity-related outcomes regressed against mtDNA quantity and (e) mtDNA quantity regressed against obesity-related outcomes.

‘Self’ represents the association between twin’s own X and Y variables whereas ‘Cotwin’ is the cross-twin cross-trait association i.e. the association between twin’s own X variable with their co-twin’s Y variable. ‘Adjusted’ refers to the regression coefficients derived from the ICE FALCON Model 3 that includes both twin’s own and their cotwin X variables.

HDL = high-density lipoprotein; Subcut. = subcutaneous; Ia. = intra-abdominal

## Discussion

In this study, we identified a significant association between mtDNAq and DNA methylation at *SH3BP4*, which correlated with its gene expression levels. This methylation site was also linked to numerous obesity-related traits, specifically those that measure body fat composition and insulin sensitivity. Leveraging our data on monozygotic twin pairs, we identified potentially causal associations from mtDNAq and obesity-related outcomes to *SH3BP4* methylation in adipose tissue.

The dynamic interplay between mitochondria and the nucleus plays a pivotal role in responding to diverse extracellular signals and metabolic conditions, as suggested by our findings: modifications in mtDNAq may precede changes in *SH3BP4* methylation, which would indicate a retrograde signaling from mitochondria to nuclear DNA methylation. This observation aligns with prior research, although limited, showing the impact of reduced mtDNAq on nuclear DNA methylation in human embryonic kidney cell lines12. One of the main hypotheses for retrograde signaling includes the importance of mitochondria in the methionine cycle, and in the production of S-adenosylmethionine. S-adenosylmethionine is a primary methyl donor in cells, which interacts with DNA methyltransferases (DNMT) and therefore can affect DNA methylation17. However, the targeted nature of the mitochondria-dependent methylation at specific genomic loci, such as *SH3BP4*, requires further investigation, along with its potential functional implications.

Our study indicated that the DNA methylation of *SH3BP4* and *DHRS3* associated with mtDNAq may be reflected at the gene expression level in adipose tissue. Specifically, methylation of *cg19998400* located at the 5’UTR region of *SH3BP4* correlated positively with *SH3BP4* expression whereas methylation at the gene body (*cg17468563*) of *DHRS3* displayed a negative correlation, which is contrary to the previously observed general pattern of positive association between gene body methylation and expression18,19. However, recent research indicates a far more complex relationship between these two factors, heavily influenced by the underlying genomic context20–22. Although our study does not pinpoint the molecular mechanisms underlying the association, we underscore the potential implication of these genes and their methylation status in relation to varying mtDNAq in obesity. Moreover, we demonstrated that DNA methylation of both *SH3BP4* and *DHRS3* accounted for a substantial proportion of the variation in the expression of these genes, alongside shared familial factors within the twin pairs. These factors encompass both genetic elements, such as eQTLs, and environmental factors, like age and lifestyle, which cannot be distinguished apart using MZ twin pairs only. Consequently, future studies with the inclusion of a cohort of dizygotic twins or other family members could provide insights into the relative significance of genes and environment.

The link we identified between the mtDNAq-associated methylation site in *SH3BP4* and various obesity-related outcomes, particularly those assessing insulin sensitivity and body fat composition, implies a potential role of mtDNAq-induced methylation in the etiology of obesity. *SH3BP4* acts as a negative regulator in many signaling pathways such as mTORC1 (mammalian target of rapamycin complex 1), a key promoter of cell growth and proliferation23, and has been observed to promote adipogenesis, possibly by regulating mitochondrial functions24–26. Previous research has demonstrated that mTOR signaling is compromised in obesity, potentially influenced by factors such as diet quality27–29 and oxidative stress30. Still, the precise contribution of *SH3BP4* methylation and gene expression to cellular function and disease development remains unclear.

Our analysis indicates that changes in DNA methylation at *SH3BP4* may result from alterations in body insulin sensitivity and intra-abdominal fat accumulation. Previous studies conducted on blood have also suggested that methylation may be a consequence rather than a cause of obesity31,32. Intriguingly, our findings suggest that only certain obesity-related outcomes, perhaps those more closely tied to metabolic disruptions such as insulin resistance and ectopic fat, are causally associated to *SH3BP4* methylation, while others closely related, such as BMI, are not. This may reflect the high metabolic heterogeneity often observed among people with similar BMI33. Nevertheless, the discrepancy warrants further investigation into the pathways that result in changes in *SH3BP4* methylation in the context of obesity.

The connection between mtDNAq and obesity has been established10,34,35, yet the exploration of causal pathways has been limited. We identified a dualistic relationship between mtDNAq and obesity-related outcomes which pointed to either causation from mtDNAq to several obesity-related outcomes or the presence of unmeasured within-individual confounding. Excessive nutrient intake is a plausible factor influencing this association, as it is known to impair mitochondrial function36–38 and contribute to the excess body weight. Whether nutrient intake serves as a confounder in the association or initiates a pathway mediated by mitochondria leading to obesity remains uncertain. It is also possible that there is a circular relationship between mtDNAq and obesity. Additionally, changes in mtDNAq may manifest only after changes in other mitochondrial parameters, which were not covered in this study, thereby limiting the identification of causal associations between mitochondrial function and obesity, and DNA methylation. Despite this, we demonstrate that both mtDNAq as well as specific obesity-related outcomes may be causal to DNA methylation at *SH3BP4*, via shared or independent molecular pathways.

Aging is widely linked with a decline in mitochondrial metabolism, including reduced mtDNAq39,40. We revealed an association between mtDNAq and EAA, measured with Hannum, Horvath and PhenoAge, in adipose tissue, which can indicate that mitochondrial metabolism is one of the key components in driving biological aging. These clocks have been reported to exhibit similar transcriptional signals with one another41. While many of the epigenetic clocks were originally developed for whole blood (except for Horvath that is a multi-tissue clock), it has been shown that some of the clocks, including PhenoAge and Hannum, work fairly robustly in other tissues too41. Moreover, PhenoAge shows increased age acceleration in cells with depleted mitochondria41.

Our study underscores the importance of investigating diseases-affected tissues beyond readily available blood samples. Exploring two primary tissues affected by excess weight, we discovered that the associations between DNA methylation and mtDNAq in obesity are not uniform. Only adipose tissue DNA methylation, but not muscle, was found to be associated with mtDNAq, which may reflect the different roles of these two tissues in obesity. For instance, alterations in adipose tissue seem to be more profoundly associated with metabolic health than those in muscle in obesity11. It may be that changes in muscle emerge only after systemic changes, such as insulin resistance, that are followed by the lipid accumulation. It is important to note that our cohort comprised mostly healthy participants, and therefore our findings cannot be necessarily generalized to more severe health complications such as metabolic syndrome or type 2 diabetes. Furthermore, mtDNAq as a proxy for mitochondrial biogenesis may not be directly comparable between the two tissues.

This study encompasses several strengths. First, the carefully phenotyped twin cohort for a comprehensive range of obesity-related outcomes enables disentangling the most significant molecular phenotypes of obesity to mtDNAq. Second, the inclusion of adipose and skeletal muscle tissues broadens the examination of the impact of excess weight and other obesity-related outcomes across tissues. In addition, using monozygotic twins, we can use statistical methods such as ICE FALCON to investigate the causality of the observed associations. The ICE FALCON offers a robust approach to explore the causal relationship in the observational data using related individuals, including twins15,42–44 without the need for genetic instrumental variables. To our knowledge, there is a lack of established genetic variants or polygenic scores to estimate mtDNAq specifically in adipose and muscle tissues. This absence of genetic information on mtDNAq prevented us from using statistical methods such as MR-CoD45, which depend on genetic data.

However, there are limitations to consider that include the cross-sectional nature of the study, as well as the modest sample size, which substantially limits statistical power. Nevertheless, the uniqueness of our dataset, to our knowledge, sets it apart as the first to integrate DNA methylation and mitochondrial quantity in individuals with obesity using adipose and muscle tissues.

Overall, we demonstrate a potential causal link from adipose tissue mitochondrial metabolism to DNA methylation and expression of *SH3BP4*. Additionally, this connection holds significance in obesity, where certain outcomes related to insulin sensitivity and intra-abdominal fat were seen to contribute to *SH3BP4* methylation, influenced either by mtDNAq or through alternative pathways. We propose the existence of a complex network interconnecting DNA methylation and mitochondrial metabolism in obesity, contributing to the multifaceted nature of obesity as a phenotype. Comprehensive examination of their interplay with various metabolic parameters holds promise for advancing our understanding of the intricate metabolic landscape in obesity.

## Methods

### Study cohort

The study participants originate from the metabolic substudy from the larger Finnish Twin Cohort (FinnTwin1246, FinnTwin1647 and Older Finnish Twin Cohort48). The present substudy was designed to study obesity-related metabolism, and the participants were initially invited to participate based on their self-reported weight and height. Our data consisted of 173 individuals, comprising those with both subcutaneous adipose and muscle tissue data available (n=141), individuals with only subcutaneous adipose tissue (n=16), and those with only skeletal muscle data (n=16) (Table 1).

### Clinical data

The selected study participants were deeply phenotyped for obesity-related clinical measures as described more in detail in van der Kolk et al. 202111. Briefly, anthropometric and body composition were measured after overnight fasting. *Body mass index* (BMI) was calculated from *weight* and *height* (kg/m2), measured during on-site visits. In addition, *waist circumference* was measured and *waist-to-hip ratio* (WHR) calculated. *Fat mass*, *fat percentage* and *lean mass* were quantified using dual-energy X-ray absorptiometry. *Intra-abdominal* and *subcutaneous fat volumes* were measured using magnetic resonance imaging (MRI), and *liver fat content* using magnetic resonance spectroscopy (MRS). Supine *blood pressure* measurements were also taken (mean of three measurements).

Blood samples were obtained after overnight fasting, and concentrations of *plasma glucose*, *serum insulin*, *plasma total cholesterol*, *low-density lipoprotein* (LDL), *high-density lipoprotein* (HDL), *triglycerides*, *high sensitivity C-reactive protein* (hsCRP), *alanine aminotransferase* (ALAT) and *aspartate aminotransferase* (ASAT) were measured using standard HUSLAB clinical laboratory procedures. *Homeostatic model assessment-insulin resistance index* (HOMA-IR) and the *Matsuda index* for insulin sensitivity were calculated from the standard 4-point oral glucose tolerance test (OGTT).

Levels of different forms of physical activity (*sport, work, leisure and total*) were assessed using the Baecke questionnaire49,50.

### Sample Collection

The samples of this study have been used in previous research11,51,52. Briefly, the adipose and muscle tissue samples were collected from subcutaneous abdominal adipose tissue and vastus lateralis muscle, respectively, under local anesthesia (lidocaine). Adipose tissue samples were taken using a surgical technique or needle biopsy, and muscle samples through Bergström needle biopsy. Determination of *adipocyte volume (dm3)* has been previously described in Lapatto et al. 202352. The tissue samples for DNA/RNA extraction were snap-frozen in liquid nitrogen.

### Mitochondrial DNA quantity

DNA was extracted from adipose and muscle tissue biopsies using AllPrep DNA/RNA/miRNA Universal Kit (Qiagen). The amount of mtDNA was quantified using quantitative PCR (qPCR), targeting for two mitochondrial encoded genes *ND5* (NADH dehydrogenase 5) and *CYTB* (cytochrome b), and normalized to genomic DNA as measured from *APP* (amyloid-beta precursor protein) and *B2M* (beta-2-microglobulin). Data was processed using the 2-ΔΔCt method with qbase+ software (Biogazelle) to obtain calibrated normalized relative quantities (CNRQ) for the mtDNA quantity. In cases of missing data (due to poor sample quality) for either *ND5* or *CYTB*, we applied a stochastic regression method to input missing values, using the other gene as a reference. Subsequently, we calculated the mean of *ND5* and *CYTB*, which served as the metric for mtDNAq.

### DNA methylation data

High molecular weight DNA was extracted from adipose and muscle biopsies with AllPrep DNA/RNA/miRNA Universal Kit (Qiagen) or QIAmp DNA Mini Kit (Qiagen) and bisulfite converted with an EZ DNA Methylation Kit (ZYMO Research) according to the manufacturers’ protocol. DNA methylation was quantified using Infinium HumanMethylation450K (adipose tissue) or HumanMethylationEPIC BeadChip arrays (adipose and muscle tissues).

DNA methylation data was preprocessed and normalized with R package *meffil*53. Due to the modest sample size, adipose tissue 450K and EPIC data were preprocessed together, omitting the platform-specific probes from the analysis. After background and bias correction, we excluded bad quality samples with following criteria: i) Median difference in X and Y chromosome intensities > 3 standard deviations (SDs), ii) Median methylated vs. unmethylated intensity > 3 SDs, iii) unreliable control probes, iv) detection p-value >0.01 in more than 20% of probes and v) bead number < 3 in more than 20% of the probes.

We then applied quantile normalization to reduce technical variation by adjusting for methylation sample slide and control probe PCs. Number of PCs included was estimated from a scree plot separately for adipose and muscle data (Supplementary Fig. 2). Bad quality probes were removed with following criteria: i) Detection p-value > 0.01 in more than 20% samples, ii) bead number < 3 in more than 20% samples, ii) SNP probes and iv) ambiguous mapping probes54,55. After QC, the number of samples and probes were 153/411,585 and 155/765,201 CpG sites for adipose and muscle data, respectively. The data was then beta mixture quantile (BMIQ) normalized to adjust for type2 probe bias.

Because a considerable amount of batch effect remained from 450K and EPIC platforms in adipose tissue, we applied ComBat to minimize the effect of a platform (Supplementary Fig. 3). Methylation M-values, calculated as the log2 ratio between the methylated versus unmethylated probe intensities, were used in the statistical analysis56.

### Epigenetic age acceleration estimates

Epigenetic age for each individual was calculated from the preprocessed DNA methylation data. We used the principal component (PC) versions of the original *Hannum*57*, Horvath*58*, GrimAge*59 and *PhenoAge*60 clocks as that has been shown to remove bias caused by technical variation in certain CpGs61. Other epigenetic clocks applied were DunedinPACE62 which measures a pace of aging, as well as muscle specific epigenetic clock MEAT63 that was calculated for skeletal muscle tissue only. Epigenetic age acceleration measures, defined as the residuals from regressing an epigenetic age estimate on chronological age, were used in the statistical analyses.

### RNA sequencing data

Adipose tissue RNA sequencing data was available for a subset of the twins (n=80 individuals). The generation and preprocessing of the RNA-seq data has been described in detail elsewhere11. Briefly, RNA was extracted using AllPrep DNA/RNA/miRNA Universal Kit (Qiagen) with DNase I (Qiagen) digestion according to manufacturers’ protocol. After the calculation of RNA integrity numbers, the libraries were prepared with Illumina Stranded mRNA preparation and sequenced with Illumina HiSeq2000 platform. The data was mapped against human reference genome 38, the quality was calculated and read counts generated.

### Statistical analysis

#### Differential methylation analysis

To identify individual CpG sites associated with mtDNAq in adipose and muscle tissue, we performed an EWAS using R package *limma*64 that fits a linear model for each probe and computes moderated Bayes t-statistics. The models were adjusted for known biological and behavioral (age, sex, smoking status), and technical (beadchip date and row) covariates, as well as cell type proportions. In the absence of reference-based cell type deconvolution methods developed specifically for either adipose or muscle tissues, we applied EpiDISH65 to adjust for key known cell types: fibroblasts and epithelial cells. In addition, the adipose tissue model was adjusted for the proportion of fat cells derived from EpiDISH. The fraction of immune cells was omitted from the model due to high correlation with other cell types. The relatedness of twins in a pair was accounted for as a blocking factor in the model.

#### Associations between DNA methylation and gene expression in adipose tissue

DNA methylation is known to influence gene expression, due to which we investigated whether the identified DNA methylation sites correlated with the expression of their closest genes. We applied Pearson correlation to see to what extent DNA methylation and the expression of their respective genes were related. In addition, using linear mixed effects modeling between gene expression (outcome) and DNA methylation (predictor), and adjusting for the relatedness of the twins, we investigated how much of the variation in gene expression is explained i) by the changes in DNA methylation solely (marginal R squared) and ii) together with common familial factors shared within the twin pairs (conditional R squared).

#### Associations between obesity-related outcomes, mitochondrial DNA quantity and DNA methylation

Given that there is a well-established link between mitochondria and obesity, we were interested to explore whether the mtDNAq-associated CpG methylation was further linked to various obesity-related measures. We performed a generalized equation estimation (gee) regression using the R package *geepack*66 with ‘exchangeable’ correlation structure to account for within twin pair similarities. The CpG methylation was used as an outcome variable and each obesity-related outcome separately as a predictor variable. Those obesity-related variables that did not follow normal distribution were log10 transformed before the analysis (Supplementary Tables 3-4). The models were adjusted for age, sex and smoking, beadchip date, and row. We verified the associations between the obesity-related outcomes and mtDNAq and chose variables that exhibited statistically significant (FDR<0.05) links with both mtDNAq and CpG methylation.

#### Causal inference

To assess the potential evidence for causal relationship or common familial confounding underlying the identified associations between DNA methylation, mtDNAq and obesity-related outcomes, we applied a statistical method called ICE FALCON (Inference about Causation from Examination of FAmilial CONfounding)15. The method is based on regression models of observational data of related individuals, specifically twins, and assesses the changes in twin’s own and co-twin’s regression coefficient from with and without adjusting the counterparts’ predictor variables (Models 1-3 below). We restricted the ICE FALCON analysis for the complete MZ pairs in our cohort. ![Formula][1]</img>  Briefly, if the observed cross-twin cross-trait association (βcotwin) diminishes after adjusting for within-individual association it suggests a causal effect from the predictor variable to the outcome. Conversely, if the βcotwin appears only when conditioning on the within-individual association, it indicates a causal relationship from the outcome to the predictor. If both the βcotwin and the βself decrease after conditioning on each other, it indicates familial confounding (genetic or environmental) between the predictor and outcome. We examined βcotwin with and without conditioning on within-individual association using generalized equation estimation in R package *geepack*62 using ‘exchangeable’ correlation structure for the twin pairs. All analyses were adjusted for age, sex and smoking. Technical covariates of beadchip date and row were added as covariates to adjust analyses including DNA methylation data.

Additionally, the ICE FALCON between mtDNAq and DNA methylation sites were further adjusted for BMI. We then calculated the changes in regression coefficients, for which we estimated standard errors of changes with non-parametric bootstrapping generating 1,000 datasets having the original sample size. The models were then reversed, e.g. using the previous X variable as Y variable to gain additional evidence on the causal pathway.

## Supporting information

Supplementary Tables and Figures [[supplements/304959_file02.pdf]](pending:yes)

## Author contributions

A.He. and M.O. designed the study. A.He. performed the data analyses, visualized the data and wrote the first draft of the manuscript. K.H.P., S.H., M.O. and J.Ka. collected and generated the data. A.He. and S.H.T.L. preprocessed the omics data. J.Ku., A.Ha. and J.L. performed the MRI and P-H.G. the DEXA imaging. P.P. generated the RNA sequencing data. V.F.C.E. and S.Li. developed and assisted on the ICE FALCON analysis. M.O. and S.Lu. supervised the work. All authors critically commented and edited the final version of the manuscript.

## Competing interests

The authors declare that they have no competing interests.

## Ethics

The Ethics Committee of the Hospital District of Helsinki and Uusimaa approved the protocols of this FTC substudy data collections, and all participants provided their written informed consent. The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Declaration of Helsinki.

## Data availability

RNA sequencing data is part of the ‘Twin Study’ and deposited with the Biobank of the Finnish Institute for Health and Welfare ([https://thl.fi/en/research-and-development/thl-biobank/for-researchers/sample-collections/twin-study](https://thl.fi/en/research-and-development/thl-biobank/for-researchers/sample-collections/twin-study)). For details on accessing the data, see: [https://thl.fi/en/research-and-development/thl-biobank/for-researchers/application-process](https://thl.fi/en/research-and-development/thl-biobank/for-researchers/application-process). The rest of the FTC data is not publicly available due to the restrictions of informed consent. However, the FTC data is available through the Institute for Molecular Medicine Finland (FIMM) Data Access Committee (DAC) (fimmdac{at}helsinki.fi) for authorized researchers who have IRB/ethics approval and an institutionally approved study plan.

## Funding

This study is supported by the following funds: University of Helsinki, Faculty of Medicine, Doctoral School of Population Health (DOCPOP) (A.H.), an Australian Government Research Training Program (RTP) Scholarship (V.F.C.E.), Academy of Finland (#328685, #307339, #297908 and #251316, M.O; # 338417, S.H.; #335443, #314383, #266286, K.H.P.), Academy of Finland Center of Excellence in Complex Disease Genetics (#352792) (J.K.) and Centre of Excellence in Research on Mitochondria, Metabolism and Disease (FinMIT) (#272376) (K.H.P.), Sigrid Juselius Foundation (M.O.), Liv o Hälsa society (M.O.), Minerva Foundation (M.O.), Novo Nordisk Foundation (grants #NNF20OC0060547, #NNF17OC0027232, #NNF10OC1013354, K.H.P.; #NNF23SA0083953, S.H.), National Health and Medical Research Council Investigator Grant (GNT2017373) (S.Li.), Diabetes Research Foundation (S.H., K.H.P.), Paulo Foundation (S.H., K.H.P.), Gyllenberg Foundation (K.H.P.) Finnish Medical Foundation (K.H.P.), University of Helsinki and Helsinki University Hospital (K.H.P., S.H.) and Government Research Funds (K.H.P.).

## Acknowledgements

The authors thank all the study participants. We also acknowledge the computational resources of the Institute for Molecular Medicine Finland (FIMM) Technology Center. S.Lu. was involved in this research while affiliated with the Institute for Molecular Medicine Finland (FIMM) and is currently affiliated with Nightingale Health Plc.

*   Received April 2, 2024.
*   Revision received April 2, 2024.
*   Accepted April 3, 2024.


*   © 2024, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.Matilainen, O., Quirós, P. M. & Auwerx, J. Mitochondria and Epigenetics – Crosstalk in Homeostasis and Stress. Trends Cell Biol. 27, 453–463 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tcb.2017.02.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28274652&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 

2.  2.Oh, S.-Y. et al. Alternative Usages of Multiple Promoters of the Acetyl-CoA Carboxylase β Gene Are Related to Differential Transcriptional Regulation in Human and Rodent Tissues*. J. Biol. Chem. 280, 5909–5916 (2005).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEwOiIyODAvNy81OTA5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDQvMDMvMjAyNC4wNC4wMi4yNDMwNDk1OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

3.  3.F. C. Lopes, A. Mitochondrial metabolism and DNA methylation: a review of the interaction between two genomes. Clin. Epigenetics 12, 182 (2020).
    
    

4.  4.Atilano, S. R. et al. Mitochondrial DNA variants can mediate methylation status of inflammation, angiogenesis and signaling genes. Hum. Mol. Genet. 24, 4491–4503 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddv173&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25964427&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 

5.  5.Lee, W. T. et al. Mitochondrial DNA haplotypes induce differential patterns of DNA methylation that result in differential chromosomal gene expression patterns. Cell Death Discov. 3, 17062 (2017).
    
    

6.  6.Cortés-Pereira, E. et al. Differential Association of Mitochondrial DNA Haplogroups J and H With the Methylation Status of Articular Cartilage: Potential Role in Apoptosis and Metabolic and Developmental Processes. Arthritis Rheumatol. 71, 1191–1200 (2019).
    
    

7.  7.Laaksonen, J. et al. Mitochondrial genome-wide analysis of nuclear DNA methylation quantitative trait loci. Hum. Mol. Genet. 31, 1720–1732 (2022).
    
    

8.  8.Gong, Y. et al. Hypermethylation of Cox5a Promoter Is Associated with Mitochondrial Dysfunction in Skeletal Muscle of High Fat Diet-Induced Insulin Resistant Rats. PLoS One 9, e113784 (2014).
    
    

9.  9.Lee, W. et al. Mitochondrial DNA copy number is regulated by DNA methylation and demethylation of POLGA in stem and cancer cells and their differentiated progeny. Cell Death Dis. 6, e1664–e1664 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/cddis.2015.34&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25719248&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 

10. 10.Heinonen, S. et al. Impaired Mitochondrial Biogenesis in Adipose Tissue in Acquired Obesity. Diabetes 64, 3135–3145 (2015).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6OToiNjQvOS8zMTM1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDQvMDMvMjAyNC4wNC4wMi4yNDMwNDk1OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

11. 11.van der Kolk, B. W. et al. Molecular pathways behind acquired obesity: Adipose tissue and skeletal muscle multiomics in monozygotic twin pairs discordant for BMI. Cell Reports Med. 2, 100–226 (2021).
    
    

12. 12.Castellani, C. A. et al. Mitochondrial DNA copy number can influence mortality and cardiovascular disease via methylation of nuclear DNA CpGs. Genome Med. 12, 84 (2020).
    
    

13. 13.Ejarque, M. et al. Adipose tissue mitochondrial dysfunction in human obesity is linked to a specific DNA methylation signature in adipose-derived stem cells. Int. J. Obes. 43, 1256–1268 (2019).
    
    

14. 14.Lundgren, S. et al. BMI is positively associated with accelerated epigenetic aging in twin pairs discordant for body mass index. J. Intern. Med. 292, 627–640 (2022).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/joim.13528&link_type=DOI) 

15. 15.Li, S., Bui, M. & Hopper, J. L. Inference about causation from examination of familial confounding (ICE FALCON): a model for assessing causation analogous to Mendelian randomization. Int. J. Epidemiol. 49, 1259–1269 (2020).
    
    

16. 16.Min, J. L. et al. Genomic and phenotypic insights from an atlas of genetic effects on DNA methylation. Nat. Genet. 53, 1311–1321 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-021-00923-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 

17. 17.Maddocks, O. D. K., Labuschagne, C. F., Adams, P. D. & Vousden, K. H. Serine Metabolism Supports the Methionine Cycle and DNA/RNA Methylation through De Novo ATP Synthesis in Cancer Cells. Mol. Cell 61, 210–221 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.molcel.2015.12.014&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26774282&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 

18. 18.Lister, R. et al. Human DNA methylomes at base resolution show widespread epigenomic differences . Nature 462, 315–322 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature08514&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19829295&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000271899300037&link_type=ISI) 

19. 19.Ball, M. P. et al. Targeted and genome-scale strategies reveal gene-body methylation signatures in human cells. Nat. Biotechnol. 27, 361–368 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nbt.1533&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19329998&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000264971800022&link_type=ISI) 

20. 20.Shakhbazov, K. et al. Shared genetic control of expression and methylation in peripheral blood. BMC Genomics 17, 278 (2016).
    
    

21. 21.Hemani, G., Tilling, K. & Davey Smith, G. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLOS Genet. 13, e1007081 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1007081&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29149188&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 

22. 22.Kreibich, E., Kleinendorst, R., Barzaghi, G., Kaspar, S. & Krebs, A. R. Single-molecule footprinting identifies context-dependent regulation of enhancers by DNA methylation. Mol. Cell 83, 787–802.e9 (2023).
    
    

23. 23.Kim, Y.-M. et al. SH3BP4 Is a Negative Regulator of Amino Acid-Rag GTPase-mTORC1 Signaling. Mol. Cell 46, 833–846 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.molcel.2012.04.007&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22575674&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000306097000013&link_type=ISI) 

24. 24.Kim, J. E. & Chen, J. Regulation of Peroxisome Proliferator–Activated Receptor-γ Activity by Mammalian Target of Rapamycin and Amino Acids in Adipogenesis. Diabetes 53, 2748–2756 (2004).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6MTA6IjUzLzExLzI3NDgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wNC8wMy8yMDI0LjA0LjAyLjI0MzA0OTU5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

25. 25.Polak, P. et al. Adipose-Specific Knockout of *raptor* Results in Lean Mice with Enhanced Mitochondrial Respiration. Cell Metab. 8, 399–410 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cmet.2008.09.003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19046571&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000260675300007&link_type=ISI) 

26. 26.Huwatibieke, B. et al. Mammalian Target of Rapamycin Signaling Pathway Regulates Mitochondrial Quality Control of Brown Adipocytes in Mice. Frontiers in Physiology vol. 12 (2021).
    
    

27. 27.Hara, K. et al. Amino Acid Sufficiency and mTOR Regulate p70 S6 Kinase and eIF-4E BP1 through a Common Effector Mechanism *. J. Biol. Chem. 273, 14484–14494 (1998).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyNzMvMjMvMTQ0ODQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wNC8wMy8yMDI0LjA0LjAyLjI0MzA0OTU5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

28. 28.Sancak, Y. et al. The Rag GTPases Bind Raptor and Mediate Amino Acid Signaling to mTORC1. Science (80-.). 320, 1496–1501 (2008).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzMjAvNTg4Mi8xNDk2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDQvMDMvMjAyNC4wNC4wMi4yNDMwNDk1OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

29. 29.Kim, E., Goraksha-Hicks, P., Li, L., Neufeld, T. P. & Guan, K.-L. Regulation of TORC1 by Rag GTPases in nutrient response. Nat. Cell Biol. 10, 935–945 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncb1753&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18604198&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000258147100011&link_type=ISI) 

30. 30.DeYoung, M. P., Horak, P., Sofer, A., Sgroi, D. & Ellisen, L. W. Hypoxia regulates TSC1/2–mTOR signaling and tumor suppression through REDD1-mediated 14–3–3 shuttling. Genes Dev. 22, 239–251 (2008).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXNkZXYiO3M6NToicmVzaWQiO3M6ODoiMjIvMi8yMzkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wNC8wMy8yMDI0LjA0LjAyLjI0MzA0OTU5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

31. 31.Mendelson, M. M. et al. Association of Body Mass Index with DNA Methylation and Gene Expression in Blood Cells and Relations to Cardiometabolic Disease: A Mendelian Randomization Approach. PLoS Med. 14, e1002215–e1002215 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1002215&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28095459&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 

32. 32.Wahl, S. et al. Epigenome-wide association study of body mass index, and the adverse outcomes of adiposity. Nature 541, 81–86 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature20784&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28002404&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 

33. 33.Neeland, I. J., Poirier, P. & Després, J.-P. Cardiovascular and Metabolic Heterogeneity of Obesity. Circulation 137, 1391–1406 (2018).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjExOiIxMzcvMTMvMTM5MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA0LzAzLzIwMjQuMDQuMDIuMjQzMDQ5NTkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

34. 34.Kaaman, M. et al. Strong association between mitochondrial DNA copy number and lipogenesis in human white adipose tissue. Diabetologia 50, 2526–2533 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00125-007-0818-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17879081&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000250721400015&link_type=ISI) 

35. 35.Agius, R., Pace, N. P. & Fava, S. Reduced leukocyte mitochondrial copy number in metabolic syndrome and metabolically healthy obesity. Frontiers in Endocrinology vol. 13 (2022).
    
    

36. 36.Liesa, M. & Shirihai, O. S. Mitochondrial Dynamics in the Regulation of Nutrient Utilization and Energy Expenditure. Cell Metab. 17, 491–506 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cmet.2013.03.002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23562075&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000326265600006&link_type=ISI) 

37. 37.Zhang, F. F. et al. Dietary Patterns Are Associated with Levels of Global Genomic DNA Methylation in a Cancer-Free Population12. J. Nutr. 141, 1165–1171 (2011).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToibnV0cml0aW9uIjtzOjU6InJlc2lkIjtzOjEwOiIxNDEvNi8xMTY1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDQvMDMvMjAyNC4wNC4wMi4yNDMwNDk1OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

38. 38.Perfilyev, A. et al. Impact of polyunsaturated and saturated fat overfeeding on the DNA-methylation pattern in human adipose tissue: a randomized controlled trial1, 2, 3. Am. J. Clin. Nutr. 105, 991–1000 (2017).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpjbiI7czo1OiJyZXNpZCI7czo5OiIxMDUvNC85OTEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wNC8wMy8yMDI0LjA0LjAyLjI0MzA0OTU5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

39. 39.Mengel-From, J. et al. Mitochondrial DNA copy number in peripheral blood cells declines with age and is associated with general health among elderly. Hum. Genet. 133, 1149–1159 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00439-014-1458-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24902542&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 

40. 40.Zhang, R., Wang, Y., Ye, K., Picard, M. & Gu, Z. Independent impacts of aging on mitochondrial DNA quantity and quality in humans. BMC Genomics 18, 890 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12864-017-4287-0&link_type=DOI) 

41. 41.Liu, Z. et al. Underlying features of epigenetic aging clocks in vivo and in vitro. Aging Cell 19, e13229 (2020).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32930491&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 

42. 42.Li, S. et al. Causal effect of smoking on DNA methylation in peripheral blood: a twin and family study. Clin. Epigenetics 10, 18 (2018).
    
    

43. 43.Nissen, F. I. et al. Is There a Causal Relationship between Physical Activity and Bone Microarchitecture? A Study of Adult Female Twin Pairs. J. Bone Miner. Res. 38, 951–957 (2023).
    
    

44. 44.Bui, M. et al. Population-based cohort study of oral contraceptive use and risk of depression. Epidemiol. Psychiatr. Sci. 32, e39 (2023).
    
    

45. 45.Minică, C. C., Dolan, C. V, Boomsma, D. I., de Geus, E. & Neale, M. C. Extending Causality Tests with Genetic Instruments: An Integration of Mendelian Randomization with the Classical Twin Design. Behav. Genet. 48, 337–349 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10519-018-9904-4&link_type=DOI) 

46. 46.Rose, R. J. et al. FinnTwin12 Cohort: An Updated Review. Twin Res. Hum. Genet. 22, 302–311 (2019).
    
    

47. 47.Kaidesoja, M. et al. FinnTwin16: A Longitudinal Study from Age 16 of a Population-Based Finnish Twin Cohort. Twin Res. Hum. Genet. 22, 530–539 (2019).
    
    

48. 48.Kaprio, J. et al. The Older Finnish Twin Cohort — 45 Years of Follow-up. Twin Res. Hum. Genet. 22, 240–254 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/thg.2019.54&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31462340&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 

49. 49.Baecke, J. A., Burema, J. & Frijters, J. E. R. A short questionnaire for the measurement of habitual physical activity in epidemiological studies. Am. J. Clin. Nutr. 36, 936–942 (1982).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpjbiI7czo1OiJyZXNpZCI7czo4OiIzNi81LzkzNiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA0LzAzLzIwMjQuMDQuMDIuMjQzMDQ5NTkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

50. 50.Mustelin, L. et al. Genetic Influences on Physical Activity in Young Adults: A Twin Study. Med. Sci. Sport. Exerc. 44, (2012).
    
    

51. 51.Sillanpää, E. et al. Blood and skeletal muscle ageing determined by epigenetic clocks and their associations with physical activity and functioning. Clin. Epigenetics 13, 110 (2021).
    
    

52. 52.Lapatto, H. A. K. et al. Nicotinamide riboside improves muscle mitochondrial biogenesis, satellite cell differentiation, and gut microbiota in a twin study. Sci. Adv. 9, eadd5163 (2023).
    
    

53. 53.Min, J. L., Hemani, G., Davey Smith, G., Relton, C. & Suderman, M. Meffil: efficient normalization and analysis of very large DNA methylation datasets. Bioinformatics 34, 3983–3989 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bty476&link_type=DOI) 

54. 54.Chen, Y. et al. Discovery of cross-reactive probes and polymorphic CpGs in the Illumina Infinium HumanMethylation450 microarray. Epigenetics 8, 203–209 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4161/epi.23470&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23314698&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000314680300009&link_type=ISI) 

55. 55.Zhou, W., Laird, P. W. & Shen, H. Comprehensive characterization, annotation and innovative use of Infinium DNA methylation BeadChip probes. Nucleic Acids Res. 45, e22–e22 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkw967&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 

56. 56.Du, P. et al. Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinformatics 11, 587 (2010).
    
    

57. 57.Hannum, G. et al. Genome-wide Methylation Profiles Reveal Quantitative Views of Human Aging Rates. Mol. Cell 49, 359–367 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.molcel.2012.10.016&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23177740&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000314379400018&link_type=ISI) 

58. 58.Horvath, S. DNA methylation age of human tissues and cell types. Genome Biology vol. 14 [http://genomebiology.com/2013/14/10/R115](http://genomebiology.com/2013/14/10/R115) (2013).
    
    

59. 59.Lu, A. T. et al. DNA methylation GrimAge strongly predicts lifespan and healthspan. Aging (Albany. NY). 11, 303–327 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18632/aging.101684&link_type=DOI) 

60. 60.Levine, M. E. et al. An epigenetic biomarker of aging for lifespan and healthspan. Aging (Albany. NY). 10, 573–591 (2018).
    
    

61. 61.Higgins-Chen, A. T. et al. A computational solution for bolstering reliability of epigenetic clocks: implications for clinical trials and longitudinal tracking. Nat. Aging 2, 644–661 (2022).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s43587-022-00248-2&link_type=DOI) 

62. 62.Belsky, D. W. et al. DunedinPACE, a DNA methylation biomarker of the pace of aging. Elife 11, e73420 (2022).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.73420&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35029144&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 

63. 63.Voisin, S. et al. An epigenetic clock for human skeletal muscle. J. Cachexia. Sarcopenia Muscle 11, 887– 898 (2020).
    
    

64. 64.Ritchie, M. E. et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43, e47–e47 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkv007&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25605792&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F04%2F03%2F2024.04.02.24304959.atom) 

65. 65.Teschendorff, A. E., Breeze, C. E., Zheng, S. C. & Beck, S. A comparison of reference-based algorithms for correcting cell-type heterogeneity in Epigenome-Wide Association Studies. BMC Bioinformatics 18, 105 (2017).
    
    

66. 66.Højsgaard, S., Halekoh, U. & Yan, J. The R Package geepack for Generalized Estimating Equations. J. Stat. Softw. 15, 1–11 (2005).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/jss.v015.i02&link_type=DOI)

 [1]: /embed/graphic-8.gif