Identification of a weight loss-associated causal eQTL in MTIF3 and the effects of MTIF3 deficiency on human adipocyte function

  1. Mi Huang
  2. Daniel Coral
  3. Hamidreza Ardalani
  4. Peter Spegel
  5. Alham Saadat
  6. Melina Claussnitzer
  7. Hindrik Mulder
  8. Paul W Franks  Is a corresponding author
  9. Sebastian Kalamajski  Is a corresponding author
  1. Genetic and Molecular Epidemiology Unit, Department of Clinical Sciences, Clinical Research Centre, Lund University, Sweden
  2. Department of Chemistry, Centre for Analysis and Synthesis, Lund University, Sweden
  3. Metabolism Program, Broad Institute of MIT and Harvard, United States
  4. Unit of Molecular Metabolism, Department of Clinical Sciences, Clinical Research Centre, Lund University, Sweden
  5. Department of Nutrition, Harvard T.H. Chan School of Public Health, United States

Abstract

Genetic variation at the MTIF3 (Mitochondrial Translational Initiation Factor 3) locus has been robustly associated with obesity in humans, but the functional basis behind this association is not known. Here, we applied luciferase reporter assay to map potential functional variants in the haplotype block tagged by rs1885988 and used CRISPR-Cas9 to edit the potential functional variants to confirm the regulatory effects on MTIF3 expression. We further conducted functional studies on MTIF3-deficient differentiated human white adipocyte cell line (hWAs-iCas9), generated through inducible expression of CRISPR-Cas9 combined with delivery of synthetic MTIF3-targeting guide RNA. We demonstrate that rs67785913-centered DNA fragment (in LD with rs1885988, r2 > 0.8) enhances transcription in a luciferase reporter assay, and CRISPR-Cas9-edited rs67785913 CTCT cells show significantly higher MTIF3 expression than rs67785913 CT cells. Perturbed MTIF3 expression led to reduced mitochondrial respiration and endogenous fatty acid oxidation, as well as altered expression of mitochondrial DNA-encoded genes and proteins, and disturbed mitochondrial OXPHOS complex assembly. Furthermore, after glucose restriction, the MTIF3 knockout cells retained more triglycerides than control cells. This study demonstrates an adipocyte function-specific role of MTIF3, which originates in the maintenance of mitochondrial function, providing potential explanations for why MTIF3 genetic variation at rs67785913 is associated with body corpulence and response to weight loss interventions.

Editor's evaluation

In this study, Huang et al. perform detailed functional genomics assays in cultured adipocytes to provide mechanistic insight underlying an important obesity GWAS locus. These studies not only demonstrate allele-specific effects of MITF3 as a potential causal gene for variations in rs67785913 and rs1885988 alleles, but further provide a foundational framework from bridging GWAS associations to actionable pathways. The study strengths include genetic manipulation followed by detailed biochemical characterization to mimic and test the impacts of association, where future studies potentially involving in vivo characterizations could further inform the metabolic consequences of these observations.

https://doi.org/10.7554/eLife.84168.sa0

Introduction

Over 650 million people are obese and often suffer from metabolic abnormalities, including dyslipidemia, type 2 diabetes, and hypertension (Adams et al., 2006; May et al., 2020). It is widely believed that obesity results from an interplay between genetic and environmental factors (Thomas, 2010), but the biological mechanisms behind these interactions are poorly understood.

Genetic variation (rs12016871) at MTIF3 (encoding the Mitochondrial Translation Initiation Factor 3 protein [Kuzmenko et al., 2014]) has been robustly associated with body mass index (BMI) in humans (Locke et al., 2015). Several subsequent studies have linked MTIF3 genetic variation with the response to weight loss interventions, including diet, exercise, and bariatric surgery (Papandonatos et al., 2015; Rasmussen-Torvik et al., 2015), and with weight-related effects of habitual diet (Nettleton et al., 2015). For example, analyses in two of the world’s largest randomized controlled weight loss trials (Diabetes Prevention Program [DPP] and Look AHEAD) found that homozygous minor allele carriers (rs1885988) were slightly more prone to weight gain in the control arm, yet achieved significantly greater weight loss at 12-month post-randomization and retained lost weight longer (18–36 months) than major allele carriers (Papandonatos et al., 2015). Elsewhere, the same locus has been associated with greater and more sustained weight loss following bariatric surgery (Rasmussen-Torvik et al., 2015).

Mtif3 loss in the mouse results in cardiomyopathy owing to impaired translation initiation from mitochondrial mRNAs and uncoordinated assembly of OXPHOS complexes in heart and skeletal muscle (Rudler et al., 2019). In the human hepatocyte-like HepG2 cell line, MTIF3 loss decreases the translation of the mitochondrial-encoded ATP synthase membrane subunit 6 (ATP6) mRNA without affecting cellular proliferation (Chicherin et al., 2020). No human genomic mutations leading to total MTIF3 deficiency have been reported, but the studies outlined above suggest that MTIF3 may influence obesity predisposition and weight loss potential by modulating mitochondrial function; thus, MTIF3 may play a key role in adipose tissue metabolic homeostasis, as adipocyte mitochondria not only provide ATP, but also impact adipocyte-specific biological processes such as adipogenesis, lipid metabolism, thermogenesis, and regulation of whole-body energy homeostasis (Gregoire et al., 1998; Boudina and Graham, 2014).

In this study, we aimed to experimentally dissect the molecular mechanisms that could underlie the correlation between MTIF3 genetic variation and weight loss intervention outcomes. We hypothesized that among the common genetic variants in MTIF3, one (or more) is causal for altered MTIF3 expression. Secondly, we hypothesized that MTIF3 content in human white adipocytes influences adipocyte-specific, obesity-related traits under basal and perturbed metabolic conditions. For the latter, we used glucose restriction to mimic the effects of in vivo lifestyle interventions focused on energy restriction and expenditure.

Results

rs67785913 is a regulatory variant for MTIF3 expression

The MTIF3 rs1885988 C allele is associated with enhanced weight loss and weight retention following lifestyle intervention trials in DPP and Look AHEAD cohorts (Papandonatos et al., 2015). In the GTEx database, the rs1885988 associates with an eQTL in subcutaneous fat (Figure 1A), with C allele carriers having significantly higher MTIF3 expression (normalized effect size: 0.15, p = 0.0000032).

Figure 1 with 1 supplement see all
Identification of rs67785913 as a causal cis-eQTL for MTIF3.

(A) Violin plot of MTIF3 expression in subcutaneous adipose tissue for rs1885988 from Genotype-Tissue Expression (GTEx) Project eQTL. (B) Same as in (A), but for rs67785913. (C) Representative Sanger sequencing traces of rs67785913 CTCT/CTCT and CT/CT clones obtained after CRISPR/Cas9-mediated allele editing and single-cell cloning. (D) Normalized Z-score plot of luciferase reporter assays using vectors carrying different DNA fragments of the MTIF3 gene cloned into pGL4.23 luciferase reporter vector. Hypothesis testing was performed by comparing the transcriptional enhancer activity of each of the 12 vectors (F1–12) to the empty vector (minP). All data were plotted as mean ± standard deviation (SD), n = 4 independent experiments, p values are presented in each graph; ordinary one-way analysis of variance (ANOVA) was used for statistical analysis. (E) Relative MTIF3 expression (mRNA) in rs67785913 allele-edited cells 2 days before, at, or 2 days post-differentiation induction (day −2, 0, and 2, respectively). n = 3 clonal populations for CTCT/CTCT genotype, n = 5 clonal populations for CT/CT genotype, error bars show SD. (F) as in (E), but for GTF3A (mRNA) expression. Two-tailed Student’s t-test was used; p values are presented in each graph.

To experimentally validate and fine map the potential causal DNA variation in the haplotype block tagged by rs1885988, we looked up all tightly linked (r2 > 0.8) single-nucleotide polymorphisms (SNPs) in HaploReg database v4.1 (Ward and Kellis, 2012). We then PCR-amplified and cloned 12 DNA fragments from that haploblock, altogether comprising the linked SNP loci, into luciferase reporter plasmids. As shown in Figure 1D, by comparing the luciferase signals with minimal promoter (minP) construct, only one DNA fragment (F11), encompassing the rs67785913 locus, could enhance luciferase transcription. Coincidentally, the rs67785913 also shows an eQTL effect on MTIF3 expression in subcutaneous adipose tissue in GTEx database, with the major CT allele associated with lower expression than the minor CTCT allele (normalized effect size: −0.16, p = 3.0 × 10−8) (Figure 1B). To demonstrate an allele-specific regulatory effect on MTIF3 expression, we then used CRISPR-Cas9 to substitute the major CT for the minor CTCT allele at the rs67785913 locus in the pre-adipocyte hWAs cell line. Due to rather low CRISPR editing efficiency of that locus, we needed to genotype over 700 single-cell clones to obtain five CT/CT and three CTCT/CTCT clones without random indels, as confirmed by Sanger sequencing (Figure 1C). We then examined MTIF3 expression in these clones at pre- and post-adipogenic differentiation induction, and found rs67785913 CTCT/CTCT to confer higher MTIF3 expression at all time points (Figure 1E), although without apparent change on adipogenic differentiation markers (Figure 1—figure supplement 1). As rs67785913 also correlates with an altered GTF3A expression in other tissues (e.g., muscle, lung), we also detected, but found no apparent difference in GTF3A expression in rs67785913-edited cells (Figure 1F).

Generating inducible Cas9-expressing pre-adipocyte cell line (hWas-iCas9)

Next, we intended to use the rs67785913-edited cells in functional genomics experiments to examine the phenotypic consequences of the eQTL. To conduct meaningful studies of gene × environment interaction, it is desirable to use similarly differentiated cells with comparable baselines (e.g., similar triglyceride or mitochondrial content). Unfortunately, marginally different passage numbers between control and experimental groups can confound adipogenic differentiation. This problem can originate during single-cell cloning to create genetic knockouts/knockins, and became apparent with our rs67785913 allele-edited cells. While the mean values of adipogenic markers were similar in both rs67785913 genotypes (Figure 1—figure supplement 1), the variation between clones of the same genotype precluded the use of these cells in gene × environment studies. To circumvent this, we instead established an inducible Cas9-expressing pre-adipocyte cell line that allowed us to knockout MTIF3 after completed adipogenic differentiation. As illustrated in Figure 2, we co-transfected hWAs with two plasmids: one encoding piggyBac transposase, and the other carrying piggyBac transposon-flanked doxycycline-inducible Cas9 and constitutively expressed puromycin resistance genes. In this setup, piggyBac transposase drives the integration of the piggyBac-transposon-flanked genes, and transgenic cells are then selected and expanded in puromycin-supplemented culture medium. We have thus obtained an hWAs cell line with doxycycline-inducible Cas9 expression and maintained adipogenic differentiation capacity (henceforth called hWAs-iCas9). The Cas9-expressing differentiated cells could then be transfected with relatively low molecular weight synthetic single guide RNAs (sgRNAs) that complex with intracellularly expressed Cas9 and target the gene exon of interest to generate random indels (in essence, gene knockouts). We used this method here to determine the functional role of MTIF3 in adipocyte biology.

The workflow of establishing hWAs-iCas9 cell line and its application in studying MTIF3 and environment interactions in vitro.

Generation of MTIF3 knockout in hWAs-iCas9 mature adipocytes

To investigate the role of MTIF3 in human adipocyte development and energy metabolism we generated stable MTIF3 knockouts in differentiated hWAs-iCas9 adipocytes. We designed Cas9-specific sgRNA to generate random indels in the exon expressed in all three MTIF3 protein-encoding transcripts (Figure 3A) and obtained a >80% reduction in MTIF3 protein levels in every experiment, as assessed by western blotting (Figure 3B–D, Figure 3—figure supplement 1, and Figure 4A, B). To assess off-target effects of CRISPR-Cas9, we also performed T7EI assays on PCR-amplified top 5 predicted off-target sites and did not observe any detectable off-targeting (data are not shown).

Figure 3 with 1 supplement see all
MTIF3 perturbation in mature adipocytes does not affect adipocyte-specific protein expression or total triglyceride content.

(A) An illustration of Cas9-specific single guide RNA (sgRNA)-binding site in the exon expressed in all three MTIF3 protein-encoding transcripts. (B) Representative Sanger sequencing of control and knockout hWAs mature adipocytes. (C) Immunoblots of adipocyte markers in scrambled control and MTIF3 knockout adipocytes, n = 5 independent experiments. (D) Quantitative analysis of MTIF3 band densities in (C). (E) Quantitative analysis of ACC, FABP4, and FAS band densities in (C). (F) Representative Oil-red O staining images of control and MTIF3 knockout in hWAs mature adipocytes. Scale bar is 200 µm. (G) Total triglyceride content in scrambled control (SC) and MTIF3 knockout (KO) cells. n = 3 independent experiments. Error bars show standard deviation in all plots. Statistical analysis was performed using two-tailed Student’s t-test, p values are presented in each graph. Uncropped blot images for (C) and raw.scn data files can be found in Figure 3—source data 1.

MTIF3 perturbation in mature adipocytes disrupts mitochondrial gene expression and OXPHOS complex assembly.

(A) Immunoblots of mitochondrial genome-encoded proteins in scrambled control and MTIF3 knockout adipocytes. (B) Quantitative analysis of band densities in (A). (C) qPCR for mitochondrial gene expression in scrambled control and MTIF3 knockout adipocytes, n = 5 independent experiments. (D) Relative mitochondrial DNA content in scrambled control and MTIF3 knockout adipocytes, n = 5 independent experiments. (E) Immunoblots of mitochondrial OXPHOS complex assembly after Blue Native-PAGE electrophoresis, n = 4 independent experiments. (F) Quantitative analysis of band densities in (E). Error bars show standard deviation in all plots. Statistical analysis was performed using two-tailed Student’s t-test, p values are presented in each graph. Uncropped blot images for (A) and raw.scn data files can be found in Figure 4—source data 1. Uncropped blot images for (E) and raw.scn data files can be found in Figure 4—source data 2.

MTIF3 knockout in mature adipocytes does not affect adipogenic marker or lipid content

Although the MTIF3 knockout in hWAs-iCas9 cells was generated after the cells were differentiated, we wanted to ensure the genetic perturbation did not affect adipogenic markers or triglyceride content, as that could confound results from downstream functional studies. Incidentally, we observed that the quantities of the adipogenic markers, including ACC, FABP4, and FAS were comparable in control and MTIF3 knockout cells (Figure 3C, E; see also Figure 3—figure supplement 1). Similarly, there were no apparent differences in Oil-red O or total triglyceride content (Figure 3F, G).

MTIF3 knockout disrupts mitochondrial DNA-encoding gene and protein expression, mitochondrial content, as well as mitochondrial OXPHOS assembly in hWAs-iCas9 adipocytes

MTIF3 is a mitochondrial translation initiation factor; thus, we examined the effects of MTIF3 ablation on differentiated hWAs adipocyte mitochondrial respiration chain. Assessed by western blotting, the MTIF3 knockout cells had significantly decreased COX II (subunit of OXPHOS complex IV) and ND2 (subunit of OXPHOS complex I), trending decrease of CYTB (subunit of OXPHOS complex III), and unchanged ATP8 (subunit of OXPHOS complex V) content (Figure 4A , B). Moreover, using qPCR, we observed an altered expression of several mitochondrial DNA-encoding genes. Specifically, MTIF3 deficiency led to higher expression of MT-ND1, MT-ND2, a trending increase of MT-ND4, and lower expression of MT-ND3, and MT-CO3 (Figure 4C). In addition, we also found significantly reduced mitochondrial content in MTIF3 knockout adipocytes (Figure 4D). Taken together, the above data suggest MTIF3 knockout disrupts mitochondrial DNA-encoding gene and protein expression.

Next, we hypothesized the above observations could have originated from the insufficient MTIF3 supply during OXPHOS complex assembly (a role previously ascribed to MTIF3 [Rudler et al., 2019]). To test this, we used Blue Native-PAGE to examine OXPHOS complexes in mitochondria isolated from MTIF3 knockout adipocytes. MTIF3 deficiency led to decreased complex III2/IV2 and IV1, and a trending decreased complex V/III2 + IV1 assembly. In contrast, OXPHOS complex II assembly was significantly increased in MTIF3 knockout cells (Figure 4E, F). Interestingly, we also observed faster-migrating undefined bands in MTIF3 knockout adipocytes (Figure 4E), which could be single chain proteins, or mistranslation or degradation products. Lastly, The OXPHOS complex I + III2 + IVn appeared to be less abundant in MTIF3 knockout mitochondria, although the bands appeared more diffuse and could not be quantified (Figure 4E).

MTIF3 knockout affects mitochondrial respiration in hWAs-iCas9 adipocytes

Having established the role of MTIF3 in adipocyte mitochondria OXPHOS complex assembly, and in mitochondrial gene expression, we then investigated the mitochondrial function in MTIF3-ablated differentiated hWAs adipocytes using Seahorse Mito Stress Test. Additionally, to avoid potential cofounders caused by the high glucose content in the differentiation medium, we adapted the cells to 1 g/l glucose growth medium for 3 days before running the assay. As shown in Figure 3, MTIF3 knockout cells exhibited lower basal oxygen consumption rate (OCR), as well as lower ATP-forming capacity, the latter estimated by calculating OCR decrease after blocking ATP synthase with oligomycin (Figure 5A–C). MTIF3 knockout cells also showed a trending decrease in maximal respiration OCR (Figure 5D). Furthermore, both MTIF3 knockout and control cells, had comparable proton leak OCR, non-mitochondrial respiration OCR and coupling efficiency (Figure 5E–G).

Cellular mitochondrial respiration in hWAs adipocytes.

(A) The average oxygen consumption rate (OCR) traces during basal respiration, and after addition of oligomycin, FCCP, and rotenone/antimycin A. (B) Basal respiration OCR, n = 4 different cell passages. (C) ATP production OCR, n = 4 different cell passages. (D) Maximal respiration OCR, n = 4 different cell passages. (E) Proton leak OCR, n = 4 different cell passages. (F) Non-mitochondrial respiration OCR, n = 4 different cell passages. (G) Coupling efficiency, n = 4 different cell passages. Error bars show standard deviation. Statistical analyses were performed using paired Student’s t-test in each condition, p values are presented in each graph.

MTIF3 knockout affects hWAs-iCas9 adipocyte endogenous fatty acid oxidation

Next, we used Seahorse to assess the endogenous fatty acid oxidation in MTIF3 knockout versus control cells treated with etomoxir (an inhibitor of carnitine palmitoyl transferase). We found that MTIF3 ablation mimics the effect of etomoxir on basal endogenous fatty acid oxidation OCR. Furthermore, while etomoxir decreases basal fatty acid oxidation OCR in control cells, it does not markedly decrease it in MTIF3 knockout cells (Figure 6A, B).

Figure 6 with 1 supplement see all
MTIF3 perturbation affects adipocyte fatty acid oxidation.

(A) A representative Seahorse oxygen consumption rate (OCR) trace for endogenous fatty acid oxidation assay. MTIF3 knockout and scrambled control adipocytes were treated with or without etomoxir for 15 min before the assay. Following the basal OCR measurement, oligomycin, FCCP (carbonyl cyanide-p-trifluoromethoxyphenylhydrazone), and rotenone + antimycin A were added sequentially to measure the detection of ATP production OCR, maximal respiration OCR and non-mitochondrial respiration OCR. (B) Basal endogenous fatty acid oxidation OCR in scrambled control (SC) and MTIF3 knockout (KO) adipocytes, n = 4 independent experiments. (C) Upper panel: workflow of glucose restriction in differentiated adipocytes; Lower left panel: total triglyceride content in scrambled control (SC) and MTIF3 knockout (KO) adipocytes in 25 mM glucose conditions; Lower right panel: triglyceride content in adipocytes cultured in glucose-restricted conditions (5, 3, and 1 mM) relative to adipocytes cultured in 25 mM glucose, n = 4 independent experiments. (D) Z-score-normalized data for glycerol release in scrambled control and MTIF3 knockout adipocytes under basal, insulin-stimulated, and isoproterenol-stimulated conditions, n = 4 independent experiments. (E) qPCR for mitochondrial and adipocyte-related gene expression in scrambled control and MTIF3 knockout adipocytes. Error bars show standard deviation in all plots. Statistical analysis was performed using two-tailed Student’s t-test, p values are presented in each graph.

MTIF3 knockout affects hWAs-iCas9 adipocyte triglyceride content after glucose restriction challenge

To mimic the interactions between MTIF3 content and dietary intervention on weight change, we generated hypertrophic control and MTIF3 knockout hWAs-iCas9 adipocytes and then used glucose-limited medium, not supplemented with free fatty acids (FFAs), to mimic energy restriction in vivo (schematic shown in Figure 6C). Triglyceride content decreased both in control and MTIF3 knockout cells after 3 days of different levels of glucose restriction when compared with 25 mM glucose medium. Interestingly, a more extensive decrease in triglyceride content occurred in control cells cultured in 1 mM glucose medium (p = 0.053), and a similar trending decrease, albeit with higher coefficient of variation, occurred in 3 and 5 mM glucose medium (Figure 6C).

MTIF3 knockout does not affect lipolysis-mediated glycerol release in hWAs adipocytes

Owing to the effects of MTIF3 ablation on triglyceride content and on fatty acid oxidation, described above, we then examined the effects of MTIF3 knockout on lipolysis. We measured basal, insulin-attenuated, and isoproterenol-stimulated glycerol release in differentiated hWAs cells. As shown in Figure 6D, in all three conditions, glycerol release in control and MTIF3 knockout cells was comparable. In addition, basal glycerol release in glucose-restricted conditions was similar (Figure 6—figure supplement 1), and significantly reduced in low glucose versus high glucose assay medium.

MTIF3 knockout affects mitochondrial function- and fatty acid oxidation-related gene expression

Next, we examined how MTIF3 ablation affects the gene expression programmes pertinent to mitochondrial function, fatty acid oxidation, lipolysis and lipid catabolism. As shown in Figure 6E, MTIF3 knockout cells had decreased expression of the mitochondria-related MT-CO1, and the fatty acid oxidation-related ACADM and ACAT1, but unchanged expression of other genes involved in mitochondrial function and lipid metabolism (TFAM, TOMM20, PRDM16, CPT1B, ABHD5, PNPLA2, and ACACB).

MTIF3 knockout results in glucose level-depending alterations in metabolism

Considering the observed effect of MTIF3 ablation on mitochondrial function and fatty acid oxidation, we assessed the metabolite profile in MTIF3 knockout cells. Using combined GC/MS and LC/MS metabolite profiling resulted in relative quantification of 110 metabolites. First, we analyzed metabolite profiles at a global level using PCA. The score plot reveals a clear systematic difference in the metabolite profile between cells in 25 mM glucose versus in glucose restriction (Figure 7A). Interestingly, differences between MTIF3 knockout and control cells at 25 mM glucose are observed along principal component 1 (PC1), whereas differences between genotypes at glucose restriction are observed along PC2, suggesting the effect of MTIF3 ablation to depend on the calorie level. Next, to identify alterations in metabolite levels underlying this differential response, we analyzed data using orthogonal projections to latent structures discriminant analysis (OPLS-DA) separately at 25 mM glucose condition (two components R2 = 0.82, Q2 = 0.66) and at glucose restriction (two components, R2 = 0.95, Q2 = 0.52). These analyses revealed systematic differences between genotypes at both growth conditions (Figure 7 B, C). Next, to examine whether the differences between genotype depended on growth condition, we combined the correlations from the two OPLS-DA models into a shared and unique structures plot (Figure 7D). These analyses revealed levels of intermediates in cytosolic metabolic pathways connected to the glycolysis, such as glycerate 3-phosphate, glycerol 2-phosphate, UDP-N-acetylglucosamine, and ribose 5-phosphate, to be lower in MTIF3 knockout cells at both glucose concentrations. Interestingly, levels of fatty acids, ranging from 9 to 17 carbons and including several odd-chain fatty acids, were lower in the MTIF3 knockout cells only at 25 mM glucose condition. At glucose-restricted conditions, levels of both essential and non-essential amino acids were lower in the knockout. Finally, we analyzed data using two-way analysis of variance (ANOVA), incorporating glucose concentration and genotype, thereby providing information on effects at the individual metabolite level. These analyses revealed 18 and 20 significantly different metabolites between control and MTIF3 knockout cells at 25 mM glucose condition and glucose-restricted conditions, respectively (q <0.05). These included ribose 5-phosphate, glycerate 3-phosphate, glycerol 2-phosphate, and glycerol 3-phosphate (Figure 7E).

Mass spectrometry-based metabolomics data for control (SC) and MTIF3 knockout (KO) cells in 25 mM glucose (NF, normal feeding) and 5 mM glucose (GR, glucose-restricted) conditions.

(A) Principal component analysis (PCA) score plot displaying the discrimination between MTIF3 knockout and control cells in normal and glucose-restricted conditions (PC1: 28%, PC2: 19%). (B) Orthogonal projections to latent structures discriminant analysis (OPLS-DA) score plot showing classification of MTIF3 knockout and control cells in 25 mM glucose condition. (C) OPLS-DA score plots showing classification of MTIF3 knockout and control cells in glucose-restricted condition. (D) Shared and unique structures (SUS) plot, based on OPLS-DA models in (B, C), showing glucose concentration-dependent differences between MTIF3 knockout and control cells. (E) Box plots showing the abundance of some of the significantly altered metabolites in normal and MTIF3 knockout cells in normal and glucose-restricted conditions. Statistical analysis was performed using two-way analysis of variance (ANOVA) test, p values are presented in each graph.

Discussion

Excessive weight gain caused by dietary excess, and its effects on adipocyte lipid metabolism, can cause life-threatening disease (Appleton et al., 2013; Denis and Obin, 2013; Chu et al., 2017; Lotta et al., 2018). Findings from clinical trials (Papandonatos et al., 2015), a bariatric surgery case series (Rasmussen-Torvik et al., 2015) and epidemiological cohorts (Nettleton et al., 2015) showed the MTIF3 variation modulates weight loss-promoting exposures on body weight (see also UK Biobank analysis in Supplementary file 1b). Here, we validated MTIF3 rs1885988 C allele correlates with higher MTIF3 expression in subcutaneous fat tissue, and our in vitro luciferase reporter assay and CRISPR-Cas9 genome editing results revealed that the tightly linked rs67785913 variant is likely to be the actual eQTL for MTIF3 expression (Figure 1).

Pre-adipocyte cell lines have been used extensively for studying adipogenic differentiation (Ruiz-Ojeda et al., 2016), but their application for lipid metabolism studies has been limited, especially in the context of gene–environment interaction. This is largely due to the variation of differentiation capacity across cell passages (Poulos et al., 2010) and differential genetic effects on adipocyte differentiation (Kamble et al., 2020), which can alter baseline phenotypes in differentiated cells. Therefore, we established an inducible Cas9-expressing human pre-adipocyte cell line (hWAs-iCas9), which enabled us to generate gene knockout of interest in differentiated adipocytes, thus circumventing these limitations (Figure 2). Using the inducible knockout cell model, we then tested interactions between MTIF3 and environmental changes (lifestyle mimetics). Our data reveal that MTIF3 deficiency mediated disrupted mitochondrial respiration, probably as a consequence of decreased OXPHOS complex assembly. This led to perturbed cellular functions, including reduced fatty acid oxidation and a trending increase in intracellular triglyceride content. These data indicate that MTIF3 plays an important role in lipid metabolism in human adipocytes.

In adipose tissue, mitochondria play an essential role not only by ensuring ATP supply but also by triggering cellular signaling pathways that require reactive oxygen species generated by OXPHOS complexes I and III (Tormos et al., 2011). These, and complexes IV and V, are partially encoded by mitochondrial DNA (Maechler and Wollheim, 2001). The post-transcriptional rate-limiting translation step can be promoted by MTIF3, as it facilitates initiation complex formation on mitochondrial 55S ribosomes in the presence of MTIF2, fMet-tRNA, and poly(A,U,G) (Bhargava and Spremulli, 2005). Previous studies have shown that loss of MTIF3 results in an imbalanced assembly of OXPHOS complexes in muscle and heart in mouse (Rudler et al., 2019) and decreased translation of ATP6 mRNA in hepatocyte-like HepG2 cells (Chicherin et al., 2020). Our data show that in adipocytes MTIF3 deficiency results in lower content of mitochondrial DNA-encoding proteins and altered expression of several mitochondrial DNA-encoding genes. Furthermore, it leads to disassembly of mitochondrial respiration OXPHOS complex, along with impaired mitochondrial respiration rate (Figures 4 and 5). Altogether, our results suggest that MTIF3 vastly affects the mitochondrial electron transport chain.

Mitochondrial dysfunction in white adipose tissue has been frequently associated with obesity (Kaaman et al., 2007), with the presumed mechanisms being decreased fatty acid oxidation and ATP production (reviewed in Heinonen et al., 2020). We have found that one causal link connecting these factors may be the MTIF3 content (Figures 5 and 6). Furthermore, a previous study described an inverse relationship between mitochondrial capacity and weight change (Jokinen et al., 2018); thus, conceivably, any genetic variation-modulated change in MTIF3 expression could influence a dietary intervention outcome. We attempted to test this hypothesis in vitro, exposing hypertrophic adipocytes to glucose restriction, thereby mimicking weight loss-promoting exposures in vivo. We found that MTIF3-deficient adipocytes exposed to glucose restriction challenge responded less through changed triglyceride content, indicating limited capacity for lipid metabolism under glucose restriction (Figure 6). Intriguingly, we did not observe altered lipolysis (glycerol release) in MTIF3 knockouts, either in the presence or absence of insulin and isoproterenol, or during glucose restriction. Here, we speculate that MTIF3 knockouts, having reduced fatty acid oxidation, re-esterify the released fatty acids to triglycerides, and thus retain more triglycerides during glucose restriction.

Mitochondrial function is robustly associated with insulin sensitivity (Böhm et al., 2020; Pietiläinen et al., 2008; Heinonen et al., 2015), both of which can be improved through lifestyle intervention aiming at weight loss (Phielix et al., 2010; Toledo et al., 2007; Larson-Meyer et al., 2006; Civitarese et al., 2007). Furthermore, impaired adipocyte differentiation derived from mitochondrial dysfunction in adipose tissue often associates with insulin resistance in humans and animal models (Sakers et al., 2022; Zhu et al., 2022). In our study, we did not observe impaired adipocyte differentiation due to MTIF3 deficiency or rs67785913 eQTL. These results, however, may not be directly translatable to in vivo conditions, in part because in vitro differentiation protocols employ relatively highly concentrated adipogenic compounds. Further in vivo studies are therefore needed to establish any link between MTIF3 content and insulin resistance. Similarly, long-term in vivo studies on the effect of altered MTIF3 expression on body weight are warranted to illuminate the translatability of the short-term or low effect size in vitro experiments (e.g., Figures 5 and 6). We envisage that the MTIF3 effect on adipocyte metabolism, while not dramatic in some of the presented data, could translate into larger effect size over time in vivo. Here, MTIF3 effect on metabolism in other tissues (e.g., muscle) would further contribute to body weight regulation. Lastly, one should bear in mind that our in vitro data were generated in cells highly depleted of MTIF3, and the extent to which a more moderate MTIF3 deficiency in vivo (e.g., conferred by common genetic variants) influences adipogenic differentiation or long-term diet-induced weight loss is currently unknown.

Finally, the metabolomic analysis added another dimension to the role of MTIF3 in regulating adipocyte metabolism (Figure 7). The lower content of glycolysis intermediates and odd-chain fatty acids in MTIF3 knockout adipocytes could indicate blunted lipogenesis, while the decreased essential and non-essential amino acid level in MTIF3 knockout cells under glucose restriction could be an adaption to the lower energy supply owing to impaired mitochondrial function (Hansson et al., 2004). The metabolite data, and the earlier described lower fatty acid oxidation capacity of MTIF3-deficient cells, suggest that MTIF3 plays a vital role in triglyceride metabolism in adipocytes, and provides further insight into the previously reported role of this protein on weight loss induced by dietary intervention (Papandonatos et al., 2015).

In summary, we experimentally demonstrated that the common genetic variant rs67785913 is a functional polymorphism that causes modulated MTIF3 expression. Our functional genomics studies demonstrate that MTIF3 is essential for mitochondrial electron transport chain complex assembly, mitochondrial function, and lipid metabolism in human adipocyte cell lines. MTIF3 content may also influence adipocyte triglyceride content in conditions that mimic dietary restriction, reflecting a gene x environment interaction. Since our findings show that higher MTIF3 content in adipocytes increases mitochondrial function, this helps explain the previously observed interaction between MTIF3 locus and weight loss interventions. Furthermore, this suggests, although it remains to be demonstrated, that people who carry genetic variants that increase MTIF3 expression (e.g., the rs67785913 CTCT allele) may benefit more from lifestyle interventions targeting weight loss.

Lastly, we established a novel, efficient, and simple method to generate gene knockouts in differentiated adipocyte cell lines, and the method should be suitable for generating other gene knockouts.

Materials and methods

Allele-specific MTIF3 expression in subcutaneous adipose tissue

Request a detailed protocol

Data and plots presented in this manuscript for allele-specific MTIF3 expression in subcutaneous adipose tissue were obtained from GTEx Portal (https://www.gtexportal.org/home/) on 10/5/2022.

Cell culture

Request a detailed protocol

A human white pre-adipocyte cell line (hWAs) was previously isolated and immortalized from human subcutaneous white adipose tissue of a female subject, aged 56 with a BMI of 30.8. The generation and characterization of hWAs were described previously (Xue et al., 2015), and the cell line was kindly shared by Professor Yu-Hua Tseng (Joslin Diabetes Center, Harvard Medical School, USA). hWAS cells. For expansion, cells were cultured in 25 mM Dulbecco’s Modified Eagle Medium (DMEM) with GlutaMAX (10566016, Thermo Fisher Scientific), 10% fetal bovine serum (FBS; HyClone, GE Healthcare, Uppsala, Sweden), and 1% (100 U/ml) penicillin/streptomycin (15140122, Thermo Fisher Scientific). The cells were passaged at 90% confluence, and tested negative for mycoplasma.

DNA isolation and luciferase reporter assays

Request a detailed protocol

Genomic DNA was isolated from hWAs cells using DNeasy Blood and Tissue kit (69506, QIAGEN) according to the manufacturer’s manual. To fine map the transcriptional regulatory regions in the MTIF3 locus, we first identified the common genetic variants which were in tight linkage disequilibrium (r2 ≥ 0.8) with the lead variant rs1885988 in HaploReg v4.1 (Ward and Kellis, 2012). The thus identified 31 SNPs were tiled down into 12 DNA segments of the MTIF3 gene, as shown in Supplementary file 1a. These segments were then PCR amplified from hWAs DNA, all ranging from 700 to 1600 bp in size (depending on PCR primer design constraints), and with all SNP loci located several hundred bp from the ends of each fragment. The PCR primers were also designed to include flanking KpnI and EcoRV sites to allow cloning into the pGL4.23 minimal promoter luc2 luciferase reporter vector (Promega). For the reporter assays, hWAs were seeded into 96-well plates, and on the following day transfected with 95 ng of the pGL4.23 vectors and 5 ng pGL4.75 CMV-Renilla reporter vectors (for normalization), using Lipofectamine 3000 (Thermo Fisher Scientific), in technical duplicates. Two days after transfection, the luc2 and Renilla signals were detected using Dual-Glo Stop&Glo reagents (E2920, Promega). The averages of technical duplicates were used to calculate luc2:Renilla ratios, which were then Z-score normalized to allow statistical evaluation across four independent experiments.

gRNAs and ssDNA design for CRISPR/Cas9 mediated editing of rs67785913 in hWAs cells

Request a detailed protocol

To edit the rs67785913 CT allele to the minor CTCT allele in hWAs cell genome, CRISPR/Cas9 D10A nickase (Alt-R S.p. Cas9 D10A Nickase V3, IDT) and two sgRNAs, and an ssDNA donor template were used. The sgRNA spacer sequences were: 5′-TTCAATAAGAAATTCCTCAA-3′ and 5′-GAAGAAAAAGGGGGGACACG-3′. The ssDNA sequence was 5′-TGTGGACTCGCAGTCTGCCCTTGAGGAATTTCTTATTGAAGAAGAAAAAGAGGGGGGACACGGGGCCCAGACCCCCAGCACCCGGCTTTCGAGCAGGCTC-3′. All oligonucleotides, sgRNAs, and ssDNA were purchased from Integrated DNA Technologies. The transfection was performed using Nucleofector 2b device (program A-033) (Lonza, Sweden) in nucleofector reagent L (Lonza, Sweden) mixed with 5 × 105 hWAs cells, 120 pmol Cas9 nickase, 104 pmol sgRNA, and 300 pmol ssDNA. To increase the homology directed repair (HDR) editing efficiency, cells were incubated at 32°C for 2 days in growth medium containing 30 μM HDR enhancer (Alt-R HDR Enhancer V2, IDT). Subsequently, cells were transferred to 37°C for 3 days. For single-cell cloning, the hWAs cells were seeded at low density (2 cells/well in a 96-well plate) and allowed to expand for 3 weeks. Then the genomic DNA was extracted using QuickExtract DNA Extraction Solution (Lucigen) from the apparent single-cell clonal populations. To identify the allele-edited homozygous clones, PCR was used to amplify the DNA fragment surrounding rs67785913 using primer pairs as below: Forward 5′–3′ GATTTGCAGGTGAGCAGACA, Reverse 5′–3′ ACTTGGAAATGGCCAAGATG; the amplicon was then subjected to Sanger sequencing to confirm the DNA sequence of each clone.

Generation of inducible CRISPR/Cas9-expressing hWAs cell line (hWAs-iCas9)

Request a detailed protocol

hWAs cells were first seeded at 80,000 cells/well in 6-well plates and transfected with 200 ng Super PiggyBac transposase (PB210PA-1, System Biosciences) and 500 ng pPB-rtTA-hCas9-puro-PB plasmid (kind gift from Dr. William Pu) (Wang et al., 2017) using Lipofectamine 3000 (Thermo Fisher Scientific). The plasmid carries a doxycycline-inducible promoter driving the expression of Cas9, and a puromycin resistance gene, all flanked by piggyBac transposon integration sequences. After 2 days, the transfected cells were selected and expanded for 3 weeks in growth medium with 1 μg/ml puromycin, to obtain cells with genomically integrated inducible Cas9 construct.

Differentiation of hWAs-iCas9 pre-adipocytes into mature adipocytes

Request a detailed protocol

hWAs-iCas9 pre-adipocytes were seeded into 24- or 96-well plates at the density of 40,000 or 8000 cells/well, respectively. After 3 days, the cells reached confluency and were then incubated for 12 days with the differentiation cocktail, with medium changes every 3 days, as described before (Shamsi and Tseng, 2017). To increase the accumulation of lipid droplets, 30 μM FFA (Linoleic Acid-Oleic Acid-Albumin) (L9655, Sigma-Aldrich) was added to the differentiation medium (Aprile et al., 2020).

CRISPR/Cas9 guide RNA design and off-targeting check

Request a detailed protocol

To generate MTIF3 knockout adipocytes, guide RNA spacer sequence targeting MTIF3 exon 5, expressed in all known MTIF3 protein-encoding transcripts (as reported at https://www.ensembl.org), was selected. The spacer sequence was 5′-GCAATAGGGGACAACTGTGC-3′, and full-length sgRNA was purchased from IDT. Furthermore, the hWAs genomic sequence surrounding the gRNA-binding site was amplified by PCR using the primers 5′-CCACTTGTCTTGGGGACAGT-3′ and 5′-CTGGGAATGGTGGTTGAATC-3′, then analyzed by Sanger sequencing to ensure sequence match between gRNA spacer and the intended target locus. The potential off-target sites were predicted using CRISPR-Cas9 guide RNA design checker (https://eu.idtdna.com), and the genomic regions surrounding the top 5 off-target sites were PCR amplified from the genomic DNA extracted from MTIF3-knockout and scramble control cells. The amplicons were then analyzed for any heteroduplexes generated by off-targeting using T7EI assay (IDT, Alt-R Genome Editing Detection Kit).

sgRNA transfection and MTIF3 knockout in mature adipocytes

Request a detailed protocol

After 12 days of differentiation, Cas9 expression was induced in mature adipocytes by adding 2 μg/ml doxycycline to the growth medium. On the following day, 30 nM pre-designed sgRNA was delivered into the cells using Lipofectamine RNAiMAX (13778075, Thermo Fisher Scientific) according to the manufacturer’s protocol; in parallel, 30 nM negative control crRNA (1072544, IDT) was used to transfect the scrambled control cells. One day post-transfection, cells were washed with phosphate-buffered saline (PBS) and incubated in normal growth medium for at least 3 days before functional assays carried out.

Oil-red O staining

Request a detailed protocol

After MTIF3 knockout, the differentiated white adipocytes were washed twice with PBS and fixed for 10–20 min with 4% buffered formalin at room temperature. The cells were then stained with Oil-red O solution for 30 min at room temperature, followed by five washes with distilled water. The stained cells were visualized using light microscopy.

Glucose restriction challenge for hWAs-iCas9 adipocytes

Request a detailed protocol

The hWAs-iCas9 adipocytes were firstly differentiated to mature adipocytes as described above (in FFA-supplemented medium), then the MTIF3 knockouts and scrambled controls were generated, also as described above. After 3 days, the mature adipocytes were incubated in DMEM medium (11966025, Thermo Fisher Scientific) without FFA, and supplemented with different glucose concentrations (5, 3, and 1 mM) for the glucose restriction test. The cells were then incubated for 3 days, and the triglyceride content was determined as described above.

RNA isolation and qPCR gene expression assays

Request a detailed protocol

Total RNA was extracted from cells using RNeasy plus Kit (74034, QIAGEN) together with Qiazol reagent (79306, QIAGEN). RNA purity was assessed using Nanodrop (Nanodrop, Wilmington, USA), and cDNA was synthesized using SuperScript IV VILO Master Mix (11756500, Thermo Fisher Scientific). Then, RT-qPCR was performed on ViiA7 qRT-PCR system (PE Applied Biosystems, Foster City, CA, USA), using predesigned Taqman assays following the manufacturer’s instructions. The Taqman assays (Thermo Fisher Scientific, Uppsala, Sweden) were: MTIF3 (Hs00794538_m1), GTF3A (Hs00157851_m1), ADIPOQ (Hs00977214_m1), PPARG (Hs01115513_m1), CEBPA (Hs00269972_s1), SREBF1 (Hs02561944_s1), FASN (Hs01005622_m1), TFAM (Hs01073348_g1), MT-CO1 (Hs02596864_g1), PRDM16 (Hs00223161_m1), TOMM20 (Hs03276810_g1), CPT1B (Hs00189258_m1), ACADM (Hs00936584_m1), ACAT1 (Hs00608002_m1), ABHD5 (Hs01104373_m1), PNP1A2 (Hs00386101_m1), ACACB (Hs01565914_m1), MT-ND1 (Hs02596873_s1), MT-ND2 (Hs02596874_g1), MT-ND3 (Hs02596875_s1), MT-ND4 (Hs02596876_g1), MT-CO2 (Hs02596865_g1), MT-CO3 (Hs02596866_g1), HPRT-1 (Hs99999909_m1), TBP (Hs00427620_m1), and RPL13A (Hs03043885_g1). The relative gene expression was calculated using the delta Ct method, and the target gene expression was normalized to the mean Ct of three reference genes HPRT-1, TBP, and RPL13A.

Western blotting

Request a detailed protocol

Cells were washed twice with ice-cold PBS and lysed in 1% sodium dodecyl sulfate buffer for 10 min, then passed through a QIAshredder (79654, QIAGEN) and centrifuged for 15 min at 14,000 × g. The supernatant was subsequently collected and protein concentration quantified using the BCA assays (23225, Thermo Fisher Scientific). To assess target protein expression, 10 μg lysates were loaded into 4–20% Mini-PROTEAN TGX Stain-Free Protein Gels (Bio-Rad Laboratories AB, Solna, Sweden) and separated, followed by transfer of polyvinylidene difluoride (PVDF) membranes (1704156, Bio-Rad Laboratories AB). After blocking in 5% bovine serum albumin (BSA) solution for 1 hr, the membranes were incubated with primary antibodies against MTIF3 (14219-1-AP, Proteintech), OXPHOS complex (45-8199, Thermo Fisher Scientific), FABP4, ACC, FAS (12589, Cell Signalling Technology), ATP8 (26723-1-AP, Proteintech), ND2 (19704-1-AP, Proteintech), CYTB (55090-1-AP, Proteintech), and corresponding horseradish peroxidase (HRP)-conjugated secondary antibodies (anti-mouse IgG, Cell Signalling Technology; anti-rabbit IgG, Cell Signaling Technology). TBS with 0.1% (vol/vol) Tween-20 was used for membrane washing, and TBS with 2% BSA was used for antibody incubation. To visualize the blots, Clarity western ECL substrate was added to the membrane and a CCD camera used to acquire images and Image Lab software (Bio-Rad Laboratories AB, Solna, Sweden) were used to develop the images. ImageJ software was used to quantify the protein bands. After detection of the protein targets, the membranes were stripped using Restore Western Blot Stripping Buffer (21059, Thermo Fisher Scientific) and blotted using anti-β-Actin antibody (4967, Cell Signaling Technology) or anti-GAPDH antibody (ab37168, Abcam).

Blue Native polyacrylamide gel electrophoresis and immunoblotting

Request a detailed protocol

Differentiated scrambled control and MTIF3 knockout hWAs-iCas9 cells were adapted to 5.5 mM glucose growth medium for 3 days to mimic the physiological glucose concentration. The Blue Native polyacrylamide gel electrophoresis (BN-PAGE) was performed as described previously (Singh and Duchen, 2022). Briefly, mitochondria were isolated using Mitochondria Isolation Kit (89874, Thermo Fisher Scientific). NativePAGE Sample Prep Kit (BN2008, Invitrogen) was then used for mitochondrial protein extraction and BN-PAGE sample preparation. For native gel electrophoresis, 20 µg mitochondrial protein was loaded to precast 3–12% gradient Blue Native gels (BN1001, Invitrogen) and separated according to the manufacturer’s instructions. The proteins were then electroblotted onto PVDF membrane and probed with anti-OXPHOS antibody cocktail (45-8199, Thermo Fisher Scientific) and a corresponding secondary antibody. The blots were then imaged and analyzed as described above.

Relative mitochondrial content measurement

Request a detailed protocol

To examine the effects of MTIF3 knockout on mitochondrial biogenesis in white adipocytes, relative amount of mtDNA was quantified using a qPCR-based method described previously (Ajaz et al., 2015). Briefly, total DNA was extracted and quantified using QIAamp DNA Mini Kit (catalogue number: 56304, QIAGEN) from scrambled control and MTIF3 knockout cells. For qPCR, equal amounts of total DNA from each sample were mixed with SYBR Green master mix (catalogue number: A25742, Thermo Fisher Scientific) and with primers targeting mitochondrial and nuclear genes, then the samples were run on ViiA7 qRT-PCR system (PE Applied Biosystems, Foster City, CA, USA). The relative mtDNA content was calculated as ΔCt (Ct of nuclear target − Ct of mitochondrial target).

Mitochondrial function in MTIF3 knockout adipocytes

Request a detailed protocol

To directly assess the effects of MTIF3 on mitochondrial respiration in adipocytes we used the Seahorse XF (Seahorse Bioscience, North Billerica, MA) to measure cellular respiration OCR under different conditions. hWAs-iCas9 cells were seeded at 8000 cells/well in a Seahorse 24-well plate, then differentiated and induced for MTIF3 knockout or with scrambled control, as described above. Then, cells were adapted in 1 g/l growth medium (31885049, Thermo Fisher Scientific) for 3 days. Mitochondrial function was then assessed using the Seahorse XF-24 instrument according to a protocol optimized for the adipocyte cell line. Briefly, to measure OCR independent of oxidative phosphorylation, 2 μM oligomycin (O4876, Sigma-Aldrich) was added to the cells. Subsequently, 2 μM FCCP (carbonyl cyanide-p-trifluoromethoxyphenylhydrazone) (C2920, Sigma-Aldrich) and 5 μM respiratory chain inhibitors: rotenone (R8875, Sigma-Aldrich) and antimycin A (A8674, Sigma-Aldrich) were added to measure maximal respiration and basal rates of non-mitochondrial respiration. Cells were then frozen at −80°C for at least 4 hr, then the plate was dried, and DNA was extracted with CyQUANT Cell Lysis Buffer (C7027, Thermo Fisher Scientific). Total DNA was then quantified by Quant-iT PicoGreen dsDNA Assay Kit (P7589, Thermo Fisher Scientific) against a lambda DNA-generated standard curve.

Endogenous long chain fatty acid oxidation in adipocytes

Request a detailed protocol

The Seahorse mitochondrial analyzer was used to test the effects of MTIF3 loss on endogenous long chain fatty acid oxidation in adipocytes. Prior to the assay, adipocytes were incubated overnight with substrate-limited medium: DMEM (A14430, Thermo Fisher Scientific); 0.5 mM glucose (103577-100, Angilent); 1.0 mM glutamine (103579-100, Angilent); 0.5 mM carnitine (C0283, Sigma-Aldrich); 1% FBS (SV30160.03, HyClone). On day of the assay, the substrate-limited medium was replaced with FAO assay medium: 1× Krebs-Henseleit Buffer (KHB) was supplemented with 2.5 mM glucose, 0.5 mM carnitine, and 5 mM N-2-hydroxyethylpiperazine-N-2-ethane sulfonic acid (HEPES), and the pH was adjusted to pH 7.4 with NaOH. The cells were then treated for 15 min with either 40 μM etomoxir (E1905, Sigma-Aldrich) or only with the solvent (dimethyl sulfoxide, DMSO). Etomoxir inhibits carnitine palmitoyltransferase (CPT)-1 and diglyceride acyltransferase (DGAT) activity in mitochondria, and thus inhibits mitochondrial fatty acid oxidation (Xu et al., 2003; Griesel et al., 2010). The OCR was then measured as described above.

Mass spectrometry-based metabolite profiling

Request a detailed protocol

The mature hWAs-iCas9 adipocytes were firstly induced for MTIF3 knockout, followed by glucose restriction challenge as described above. The cells were quenched on dry ice and metabolites were extracted using a previously optimized protocol (Danielsson et al., 2010).

For analysis of low molecular weight metabolites, extracts were reconstituted in 100 µl of MeOH/water (8/2, vol/vol) and 60 µl was transferred to new Eppendorf tubes and evaporated to dryness using a miVac concentrator (SP Scientific, NY) for 3 hr at 30°C. Dried samples were methoximated using 20 µl of methoxyamine hydrochloride in pyridine (Thermo Scientific, MA) by shaking at 3000 rpm for 30 min at room temperature (VWR, PA). Afterward, 20 µl of N-methyl-N-(trimethylsilyl) trifluoroacetamide (MSTFA) + 1% trimethylsilyl chloride (Thermo Scientific, MA) was added to each sample and shaken at 3000 rpm at room temperature for 1 hr. Samples were transferred to glass vials and immediately analyzed using an Agilent 6890 gas chromatograph connected to an Agilent 5975CL VL MSD mass spectrometer controlled by MassHunter Workstation software 10.0 (Agilent, Atlanta, GA). One µl sample was injected at 270°C on an HP-5MS column (30 m length, 250 µm ID, 0.25 µm phase thickness). with a helium gas flow rate of 1 ml/min and a temperature gradient starting at 70°C for 2 min, increasing 15°C/min to 320°C and held for 2 min. Data were acquired using electron ionization at 70 eV in either full scan (50–550 m/z) or single ion monitoring mode. The MS-DIAL version 4.7 was used for raw peak extraction, peak alignment, deconvolution, peak annotation, and integration of peaks.

Amino acids and free fatty acids were chemically derivatized and analyzed using a previously described method (Meng et al., 2021). Briefly, 40 µl of the samples were mixed with 20 µl of 3-nitrophenylhydrazine (3-NPH) (Sigma-Aldrich, MO), followed by addition of 20 µl of 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide hydrochloride (EDC) (Thermo Scientific, MA) and shaking at 3000 rpm at room temperature for 1 hr. Samples were analyzed using an Agilent 1260 ultra-performance liquid chromatograph coupled with an Agilent 6495 tandem mass spectrometer and controlled by MassHunter version 8.0 (Agilent Technologies, CA). Three µl sample was injected on an Agilent Eclipse RRHD C18 column (2.1 × 150 mm, 1.8 µm) (Agilent Technologies, CA) with a flow rate of 0.6 ml/min and a column oven temperature of 50°C. The mobile phases A and B were 0.1% formic acid (Fisher Chemical, Prague, Czech Republic) in Milli-Q water (Merck, Millipore, MO) and acetonitrile (VWR, Paris, France), respectively. Gradient elution was performed as follows: held at 5% B from 0 to 1 min, changed linearly to 90% B in 10 min, changed from 90% B to 100% B in 13 min, held at 100% B for 2 min, returned to 5% B (initial condition) in 0.1 min, and held at 5% B for 2 min. Analyses were conducted in negative electrospray ionization mode (ESI) mode with the nebulizer gas pressure set at 20 psi, ion capillary voltage at 2500 V, gas temperature at 150°C, and sheath gas temperature at 250°C. Data were recorded in multiple reaction monitoring (MRM) mode, with two transitions for each analyte.

Lipolysis quantification in differentiated hWAs cells

Request a detailed protocol

Differentiated scrambled control or MTIF3 knockout cells were washed twice with PBS and then incubated with DMEM containing 2% free fatty acid-free BSA for 2 hr. For the insulin or isoproterenol-stimulated lipolysis, 100 nM insulin (I2643, Sigma-Aldrich) or 10 μM isoproterenol (1351005, Sigma-Aldrich) was added in the medium separately. After the incubation, the medium was collected, and the glycerol content was measured using Glycerol-Glo Assay (J3150, Promega).

Total triglyceride measurement

Request a detailed protocol

Triglyceride-Glo Assay kit (J3161, Promega) was used to quantify total triglyceride content in scrambled control or MTIF3 knockout cells cultured either in 25 mM glucose or glucose restriction medium. Briefly, cells were collected in 50 μl kit lysis buffer at room temperature for 1 hr. Then 2 μl lysate was mixed with 8 μl glycerol lysis solution with lipase, and incubated at 37°C for 30 min. Subsequently, 10 μl glycerol solution was mixed with 10 μl glycerol detection solution supplemented with reductase substrate and kinetic enhancer, and transferred into a 384-well plate. After 1-hr incubation at room temperature, the luminescence was detected using CLARIOstar plate reader (BMG Labtech, Germany), and the triglyceride concentration was calculated using a standard curve generated from glycerol standards and normalized to total protein measured using BCA assays (23227, Thermo Fisher Scientific).

Statistics

For each assay, the number of biological and technical replicates, standard deviation and statistical significance are reported in the figure legends. Hypothesis tests were performed using two-tailed Student’s t-test, one-way ANOVA, or paired t-test. A nominal p value of <0.05 was considered statistically significant. All analyses were undertaken using Prism GraphPad 9.0 software (La Jolla California, USA), SIMCA 17.0 (Sartorius Stedim Data Analytics, Malmö, Sweden), Rstudio 1.4, and Microsoft Excel 365. For the metabolome data, ANOVA (aov) was performed in R with genotype and glucose concentration as independent variables with Tukey’s test post hoc (TukeyHSD). Significance was defined as q < 0.05 using multiple testing adjustment according to the false discovery rate (p.adjust).

Appendix 1

Appendix 1—key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Gene
(Homo sapiens)
MTIF3UCSC Genome BrowserGRCh38/hg38
Cell line
(Homo sapiens)
hWAsTseng laboratory at Joslin Diabetes CenterXue et al., 2015
Cell line
(Homo sapiens)
hWAs-iCas9This paperCell line maintained at Lund University
Diabetes Center
AntibodyAnti-MTIF3 (rabbit polyclonal antibody)ProteintechCat: 14219-1-APWB (1:2000)
AntibodyAnti-OXPHOS antibody cocktail (mouse polyclonal antibody)Thermo Fisher ScientificCat: 45-8199WB (1:1000)
AntibodyAnti-FABP4 (rabbit polyclonal antibody)Cell Signalling TechnologyCat: 12589WB (1:1000)
AntibodyAnti-ACC rabbit polyclonal antibodyCell Signalling TechnologyCat: 12589WB (1:1000)
AntibodyAnti-FAS rabbit polyclonal antibodyCell Signalling TechnologyCat: 12589WB (1:1000)
AntibodyAnti-ATP8
(rabbit polyclonal antibody)
ProteintechCat: 26723-1-APWB (1:2000)
AntibodyAnti-ND2 (rabbit polyclonal antibody)ProteintechCat: 19704-1-APWB (1:2000)
AntibodyAnti-CYTB (rabbit polyclonal antibody)ProteintechCat: 55090-1-APWB (1:2000)
AntibodyAnti-β-Actin (rabbit polyclonal antibody)Cell Signaling TechnologyCat: #4967WB (1:10,000)
AntibodyAnti-GAPDH (rabbit polyclonal antibody)AbcamCat: ab37168WB (1:10,000)
AntibodyAnti-rabbit IgG, HRP-linked Antibody (goat polyclonal antibody)Cell Signaling TechnologyCat: #7074WB (1:10,000)
AntibodyAnti-mouse IgG, HRP-linked Antibody (horse polyclonal antibody)Cell Signaling TechnologyCat: #7076WB (1:10,000)
Recombinant DNA reagentSuper PiggyBac transposase (plasmid)System BiosciencesPB210PA-1
Recombinant DNA reagentpGL4.23 vectorsPromegaE8411
Recombinant DNA reagentpGL4.75 CMV-Renilla reporter vectorsPromegaE6931
Recombinant DNA reagentpPB-rtTA-hCas9-puro-PB plasmiddoi:10.1038/nprot.2016.152
Sequence-based reagentPCR primer (Forward) for rs67785913 genotypingIDT5′–3′: GATTTGCAGGTGAGCAGACA
Sequence-based reagentPCR primer (Reverse) for rs67785913 genotypingIDT5′–3′: ACTTGGAAATGGCCAAGATG
Sequence-based reagentsgRNA for rs67785913 editingIDTSpacer sequence: 5′-TTCAATAAGAAATTCCTCAA-3′
Sequence-based reagentsgRNA for rs67785913 editingIDTSpacer sequence: 5′-GAAGAAAAAGGGGGGACACG-3
Sequence-based reagentDonor template for rs67785913 editingIDTssDNA sequence:
5′’TGTGGACTCGCAGTCTGCCCTTGAGGAATTTCTTATTGAAGAAGAAAAAGAGGGGGGACACGGGGCCCAGACCCCCAGCACCCGGCTTTCGAGCAGGCTC-3′
Sequence-based reagentsgRNA against MTIF3IDTDesign ID: Hs.Cas9.MTIF3.1.ABSpacer sequence:
5′-GCAATAGGGGACAA
CTGTGC-3′
Sequence-based reagentTaqman assay for MTIF3Thermo Fisher ScientificHs00794538_m1
Sequence-based reagentTaqman assay for GTF3AThermo Fisher ScientificHs00157851_m1
Sequence-based reagentTaqman assay for ADIPOQThermo Fisher ScientificHs00977214_m1
Sequence-based reagentTaqman assay for PPARGThermo Fisher ScientificHs01115513_m1
Sequence-based reagentTaqman assay for CEBPAThermo Fisher ScientificHs00269972_s1
Sequence-based reagentTaqman assay for SREBF1Thermo Fisher ScientificHs02561944_s1
Sequence-based reagentTaqman assay for FASNThermo Fisher ScientificHs01005622_m1
Sequence-based reagentTaqman assay for TFAMThermo Fisher ScientificHs01073348_g1
Sequence-based reagentTaqman assay for MT-CO1Thermo Fisher ScientificHs02596864_g1
Sequence-based reagentTaqman assay for PRDM16Thermo Fisher ScientificHs00223161_m1
Sequence-based reagentTaqman assay for TOMM20Thermo Fisher ScientificHs03276810_g1
Sequence-based reagentTaqman assay for CPT1BThermo Fisher ScientificHs00189258_m1
Sequence-based reagentTaqman assay for ACADMThermo Fisher ScientificHs00936584_m1
Sequence-based reagentTaqman assay for ACAT1Thermo Fisher ScientificHs00608002_m1
Sequence-based reagentTaqman assay for ABHD5Thermo Fisher ScientificHs01104373_m1
Sequence-based reagentTaqman assay for PNP1A2Thermo Fisher ScientificHs00386101_m1
Sequence-based reagentTaqman assay for ACACBThermo Fisher ScientificHs01565914_m1
Sequence-based reagentTaqman assay for MT-ND1Thermo Fisher ScientificHs02596873_s1
Sequence-based reagentTaqman assay for MT-ND2Thermo Fisher ScientificHs02596874_g1
Sequence-based reagentTaqman assay for MT-ND3Thermo Fisher ScientificHs02596875_s1
Sequence-based reagentTaqman assay for MT-ND4Thermo Fisher ScientificHs02596876_g1
Sequence-based reagentTaqman assay for MT-CO2Thermo Fisher ScientificHs02596865_g1
Sequence-based reagentTaqman assay for MT-CO3Thermo Fisher ScientificHs02596866_g1
Sequence-based reagentTaqman assay for HPRT-1Thermo Fisher ScientificHs99999909_m1
Sequence-based reagentTaqman assay for TBPThermo Fisher ScientificHs00427620_m1
Sequence-based reagentTaqman assay for RPL13AThermo Fisher ScientificHs03043885_g1
Commercial assay or kitDNeasy Blood and Tissue kitQIAGEN69506
Commercial assay or kitDual-Glo Stop&Glo reagentsPromegaE2920
Commercial assay or kitAlt-R S.p. Cas9 D10A Nickase V3IDT1081058
Commercial assay or kitNucleofector reagent LLonzaVCA-1005
Commercial assay or kitAlt-R HDR Enhancer V2IDT10007921
Commercial assay or kitQuickExtract DNA Extraction SolutionLucigenQE09050
Commercial assay or kitAlt-R Genome Editing Detection KitIDT1075932
Commercial assay or kitMitochondrial isolation kitThermo Fisher Scientific89874
Commercial assay or kitNativePAGE Sample Prep KitInvitrogenBN2008
Commercial assay or kitQuant-iT PicoGreen dsDNA Assay KitThermo Fisher ScientificP7589
Commercial assay or kitTriglyceride-Glo Assay kitPromegaJ3161
Commercial assay or kitGlycerol-Glo AssayPromegaJ3150
Chemical compound, drugInsulinSigma-AldrichI2643
Chemical compound, drugIsoproterenolSigma-Aldrich1351005
Chemical compound, drugGlucoseAngilent103577-100
Chemical compound, drugGlutamineAngilent103579-100
Chemical compound, drugCarnitineSigma-AldrichC0283
Chemical compound, drugPyridineThermo Scientific019378.K2
Chemical compound, drugN-Methyl-N-(trimethylsilyl) trifluoroacetamideThermo ScientificA13141.22
Chemical compound, drugTrimethylsilyl chlorideThermo ScientificA12535.30
Chemical compound, drug3-NitrophenylhydrazineSigma-AldrichN21804
Chemical compound, drug1-Ethyl-3-(3-dimethylaminopropyl)carbodiimide hydrochlorideThermo Scientific22980
Chemical compound, drugFormic acidFisher ChemicalA117-50
Chemical compound, drugEtomoxirSigma-AldrichE1905
Chemical compound, drugOligomycinSigma-AldrichO4876
Chemical compound, drugFCCPSigma-AldrichC2920
Chemical compound, drugRotenoneSigma-AldrichR8875
Chemical compound, drugAntimycin ASigma-AldrichA8674

Data availability

All data generated or analyzed in this study are included in the figures and the source data files. Source data files are provided for Figures 3 and 4, and for Figure 3—figure supplement 1.

References

    1. Locke AE
    2. Kahali B
    3. Berndt SI
    4. Justice AE
    5. Pers TH
    6. Day FR
    7. Powell C
    8. Vedantam S
    9. Buchkovich ML
    10. Yang J
    11. Croteau-Chonka DC
    12. Esko T
    13. Fall T
    14. Ferreira T
    15. Gustafsson S
    16. Kutalik Z
    17. Luan J
    18. Mägi R
    19. Randall JC
    20. Winkler TW
    21. Wood AR
    22. Workalemahu T
    23. Faul JD
    24. Smith JA
    25. Zhao JH
    26. Zhao W
    27. Chen J
    28. Fehrmann R
    29. Hedman ÅK
    30. Karjalainen J
    31. Schmidt EM
    32. Absher D
    33. Amin N
    34. Anderson D
    35. Beekman M
    36. Bolton JL
    37. Bragg-Gresham JL
    38. Buyske S
    39. Demirkan A
    40. Deng G
    41. Ehret GB
    42. Feenstra B
    43. Feitosa MF
    44. Fischer K
    45. Goel A
    46. Gong J
    47. Jackson AU
    48. Kanoni S
    49. Kleber ME
    50. Kristiansson K
    51. Lim U
    52. Lotay V
    53. Mangino M
    54. Leach IM
    55. Medina-Gomez C
    56. Medland SE
    57. Nalls MA
    58. Palmer CD
    59. Pasko D
    60. Pechlivanis S
    61. Peters MJ
    62. Prokopenko I
    63. Shungin D
    64. Stančáková A
    65. Strawbridge RJ
    66. Sung YJ
    67. Tanaka T
    68. Teumer A
    69. Trompet S
    70. van der Laan SW
    71. van Setten J
    72. Van Vliet-Ostaptchouk JV
    73. Wang Z
    74. Yengo L
    75. Zhang W
    76. Isaacs A
    77. Albrecht E
    78. Ärnlöv J
    79. Arscott GM
    80. Attwood AP
    81. Bandinelli S
    82. Barrett A
    83. Bas IN
    84. Bellis C
    85. Bennett AJ
    86. Berne C
    87. Blagieva R
    88. Blüher M
    89. Böhringer S
    90. Bonnycastle LL
    91. Böttcher Y
    92. Boyd HA
    93. Bruinenberg M
    94. Caspersen IH
    95. Chen Y-DI
    96. Clarke R
    97. Daw EW
    98. de Craen AJM
    99. Delgado G
    100. Dimitriou M
    101. Doney ASF
    102. Eklund N
    103. Estrada K
    104. Eury E
    105. Folkersen L
    106. Fraser RM
    107. Garcia ME
    108. Geller F
    109. Giedraitis V
    110. Gigante B
    111. Go AS
    112. Golay A
    113. Goodall AH
    114. Gordon SD
    115. Gorski M
    116. Grabe H-J
    117. Grallert H
    118. Grammer TB
    119. Gräßler J
    120. Grönberg H
    121. Groves CJ
    122. Gusto G
    123. Haessler J
    124. Hall P
    125. Haller T
    126. Hallmans G
    127. Hartman CA
    128. Hassinen M
    129. Hayward C
    130. Heard-Costa NL
    131. Helmer Q
    132. Hengstenberg C
    133. Holmen O
    134. Hottenga J-J
    135. James AL
    136. Jeff JM
    137. Johansson Å
    138. Jolley J
    139. Juliusdottir T
    140. Kinnunen L
    141. Koenig W
    142. Koskenvuo M
    143. Kratzer W
    144. Laitinen J
    145. Lamina C
    146. Leander K
    147. Lee NR
    148. Lichtner P
    149. Lind L
    150. Lindström J
    151. Lo KS
    152. Lobbens S
    153. Lorbeer R
    154. Lu Y
    155. Mach F
    156. Magnusson PKE
    157. Mahajan A
    158. McArdle WL
    159. McLachlan S
    160. Menni C
    161. Merger S
    162. Mihailov E
    163. Milani L
    164. Moayyeri A
    165. Monda KL
    166. Morken MA
    167. Mulas A
    168. Müller G
    169. Müller-Nurasyid M
    170. Musk AW
    171. Nagaraja R
    172. Nöthen MM
    173. Nolte IM
    174. Pilz S
    175. Rayner NW
    176. Renstrom F
    177. Rettig R
    178. Ried JS
    179. Ripke S
    180. Robertson NR
    181. Rose LM
    182. Sanna S
    183. Scharnagl H
    184. Scholtens S
    185. Schumacher FR
    186. Scott WR
    187. Seufferlein T
    188. Shi J
    189. Smith AV
    190. Smolonska J
    191. Stanton AV
    192. Steinthorsdottir V
    193. Stirrups K
    194. Stringham HM
    195. Sundström J
    196. Swertz MA
    197. Swift AJ
    198. Syvänen A-C
    199. Tan S-T
    200. Tayo BO
    201. Thorand B
    202. Thorleifsson G
    203. Tyrer JP
    204. Uh H-W
    205. Vandenput L
    206. Verhulst FC
    207. Vermeulen SH
    208. Verweij N
    209. Vonk JM
    210. Waite LL
    211. Warren HR
    212. Waterworth D
    213. Weedon MN
    214. Wilkens LR
    215. Willenborg C
    216. Wilsgaard T
    217. Wojczynski MK
    218. Wong A
    219. Wright AF
    220. Zhang Q
    221. Brennan EP
    222. Choi M
    223. Dastani Z
    224. Drong AW
    225. Eriksson P
    226. Franco-Cereceda A
    227. Gådin JR
    228. Gharavi AG
    229. Goddard ME
    230. Handsaker RE
    231. Huang J
    232. Karpe F
    233. Kathiresan S
    234. Keildson S
    235. Kiryluk K
    236. Kubo M
    237. Lee J-Y
    238. Liang L
    239. Lifton RP
    240. Ma B
    241. McCarroll SA
    242. McKnight AJ
    243. Min JL
    244. Moffatt MF
    245. Montgomery GW
    246. Murabito JM
    247. Nicholson G
    248. Nyholt DR
    249. Okada Y
    250. Perry JRB
    251. Dorajoo R
    252. Reinmaa E
    253. Salem RM
    254. Sandholm N
    255. Scott RA
    256. Stolk L
    257. Takahashi A
    258. Tanaka T
    259. van ’t Hooft FM
    260. Vinkhuyzen AAE
    261. Westra H-J
    262. Zheng W
    263. Zondervan KT
    264. Heath AC
    265. Arveiler D
    266. Bakker SJL
    267. Beilby J
    268. Bergman RN
    269. Blangero J
    270. Bovet P
    271. Campbell H
    272. Caulfield MJ
    273. Cesana G
    274. Chakravarti A
    275. Chasman DI
    276. Chines PS
    277. Collins FS
    278. Crawford DC
    279. Cupples LA
    280. Cusi D
    281. Danesh J
    282. de Faire U
    283. den Ruijter HM
    284. Dominiczak AF
    285. Erbel R
    286. Erdmann J
    287. Eriksson JG
    288. Farrall M
    289. Felix SB
    290. Ferrannini E
    291. Ferrières J
    292. Ford I
    293. Forouhi NG
    294. Forrester T
    295. Franco OH
    296. Gansevoort RT
    297. Gejman PV
    298. Gieger C
    299. Gottesman O
    300. Gudnason V
    301. Gyllensten U
    302. Hall AS
    303. Harris TB
    304. Hattersley AT
    305. Hicks AA
    306. Hindorff LA
    307. Hingorani AD
    308. Hofman A
    309. Homuth G
    310. Hovingh GK
    311. Humphries SE
    312. Hunt SC
    313. Hyppönen E
    314. Illig T
    315. Jacobs KB
    316. Jarvelin M-R
    317. Jöckel K-H
    318. Johansen B
    319. Jousilahti P
    320. Jukema JW
    321. Jula AM
    322. Kaprio J
    323. Kastelein JJP
    324. Keinanen-Kiukaanniemi SM
    325. Kiemeney LA
    326. Knekt P
    327. Kooner JS
    328. Kooperberg C
    329. Kovacs P
    330. Kraja AT
    331. Kumari M
    332. Kuusisto J
    333. Lakka TA
    334. Langenberg C
    335. Marchand LL
    336. Lehtimäki T
    337. Lyssenko V
    338. Männistö S
    339. Marette A
    340. Matise TC
    341. McKenzie CA
    342. McKnight B
    343. Moll FL
    344. Morris AD
    345. Morris AP
    346. Murray JC
    347. Nelis M
    348. Ohlsson C
    349. Oldehinkel AJ
    350. Ong KK
    351. Madden PAF
    352. Pasterkamp G
    353. Peden JF
    354. Peters A
    355. Postma DS
    356. Pramstaller PP
    357. Price JF
    358. Qi L
    359. Raitakari OT
    360. Rankinen T
    361. Rao DC
    362. Rice TK
    363. Ridker PM
    364. Rioux JD
    365. Ritchie MD
    366. Rudan I
    367. Salomaa V
    368. Samani NJ
    369. Saramies J
    370. Sarzynski MA
    371. Schunkert H
    372. Schwarz PEH
    373. Sever P
    374. Shuldiner AR
    375. Sinisalo J
    376. Stolk RP
    377. Strauch K
    378. Tönjes A
    379. Trégouët D-A
    380. Tremblay A
    381. Tremoli E
    382. Virtamo J
    383. Vohl M-C
    384. Völker U
    385. Waeber G
    386. Willemsen G
    387. Witteman JC
    388. Zillikens MC
    389. Adair LS
    390. Amouyel P
    391. Asselbergs FW
    392. Assimes TL
    393. Bochud M
    394. Boehm BO
    395. Boerwinkle E
    396. Bornstein SR
    397. Bottinger EP
    398. Bouchard C
    399. Cauchi S
    400. Chambers JC
    401. Chanock SJ
    402. Cooper RS
    403. de Bakker PIW
    404. Dedoussis G
    405. Ferrucci L
    406. Franks PW
    407. Froguel P
    408. Groop LC
    409. Haiman CA
    410. Hamsten A
    411. Hui J
    412. Hunter DJ
    413. Hveem K
    414. Kaplan RC
    415. Kivimaki M
    416. Kuh D
    417. Laakso M
    418. Liu Y
    419. Martin NG
    420. März W
    421. Melbye M
    422. Metspalu A
    423. Moebus S
    424. Munroe PB
    425. Njølstad I
    426. Oostra BA
    427. Palmer CNA
    428. Pedersen NL
    429. Perola M
    430. Pérusse L
    431. Peters U
    432. Power C
    433. Quertermous T
    434. Rauramaa R
    435. Rivadeneira F
    436. Saaristo TE
    437. Saleheen D
    438. Sattar N
    439. Schadt EE
    440. Schlessinger D
    441. Slagboom PE
    442. Snieder H
    443. Spector TD
    444. Thorsteinsdottir U
    445. Stumvoll M
    446. Tuomilehto J
    447. Uitterlinden AG
    448. Uusitupa M
    449. van der Harst P
    450. Walker M
    451. Wallaschofski H
    452. Wareham NJ
    453. Watkins H
    454. Weir DR
    455. Wichmann H-E
    456. Wilson JF
    457. Zanen P
    458. Borecki IB
    459. Deloukas P
    460. Fox CS
    461. Heid IM
    462. O’Connell JR
    463. Strachan DP
    464. Stefansson K
    465. van Duijn CM
    466. Abecasis GR
    467. Franke L
    468. Frayling TM
    469. McCarthy MI
    470. Visscher PM
    471. Scherag A
    472. Willer CJ
    473. Boehnke M
    474. Mohlke KL
    475. Lindgren CM
    476. Beckmann JS
    477. Barroso I
    478. North KE
    479. Ingelsson E
    480. Hirschhorn JN
    481. Loos RJF
    482. Speliotes EK
    483. LifeLines Cohort Study
    484. ADIPOGen Consortium
    485. AGEN-BMI Working Group
    486. CARDIOGRAMplusC4D Consortium
    487. CKDGen Consortium
    488. GLGC
    489. ICBP
    490. MAGIC Investigators
    491. MuTHER Consortium
    492. MIGen Consortium
    493. PAGE Consortium
    494. ReproGen Consortium
    495. GENIE Consortium
    496. International Endogene Consortium
    (2015) Genetic studies of body mass index yield new insights for obesity biology
    Nature 518:197–206.
    https://doi.org/10.1038/nature14177

Decision letter

  1. Marcus M Seldin
    Reviewing Editor; University of California, Irvine, United States
  2. David E James
    Senior Editor; University of Sydney, Australia
  3. Sean M Hartig
    Reviewer; Baylor College of Medicine, United States

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

Decision letter after peer review:

Thank you for submitting your article "Identification of a weight loss-associated causal eQTL in MTIF3 and the effects of MTIF3 deficiency on human adipocyte function" for consideration by eLife. Your article has been reviewed by 2 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and David James as the Senior Editor. The following individual involved in the review of your submission has agreed to reveal their identity: Sean M. Hartig (Reviewer #1).

The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission. We all remain enthusiastic about the implications of these functional genomics experiments characterizing the locus and most requested revisions surround more detailed evaluation primarily in cell culture models used.

Essential revisions:

As you will see below, the primary suggestions for improving the study lie in more detailed molecular and functional characterization of MTIF3 and the effects of the locus. A few of these include protein quantification and staining (oil red O) for adipocytes where MTIF3 is perturbed. In addition, the reviewers requested a few text-based clarifications regarding the interpretation of the locus, listed below. Obviously, if you have any questions feel, please free to reach out.

Reviewer #1 (Recommendations for the authors):

It is important to point out that adipocyte differentiation and the expression of the associated metabolic programs critically contribute to the maintenance of insulin sensitivity in obesity. Defects in adipocyte differentiation derived from mitochondrial functions that reduce ATP synthesis and storage of surfeit carbons in fat tissue almost always associate with insulin resistance in animal models and people (see Trends Cell Biol 32:351-364, 2022; Cell. 185:419-446, 2022). The in vitro studies can add to understanding the important roles of MTIF3, but the authors should use more caution when extrapolating to in vivo effects of the eQTL.

Reviewer #2 (Recommendations for the authors):

1. The authors should give their hypothesis and interpretation of their results in the manuscript. There are no clear connections between results, which makes it difficult to understand the rationale of their experiments.

2. The labels of Figure 1 do not match the figure legends.

3. In figure 1F, what is the readout of the relative MTIF3expression? Are they based on western blot? If yes, it would be better to include the blog images.

4. In the method, the authors also added extra lipids in the differentiated medium to increase the lipid droplets. I'm curious about the effect of lipid deprivation in KO cells rather than glucose deprivation. The metabolomics showed the change of fatty acids at full medium. There might be some changes in fatty acid utilization in KO cells.

https://doi.org/10.7554/eLife.84168.sa1

Author response

Essential revisions:

As you will see below, the primary suggestions for improving the study lie in more detailed molecular and functional characterization of MTIF3 and the effects of the locus. A few of these include protein quantification and staining (oil red O) for adipocytes where MTIF3 is perturbed. In addition, the reviewers requested a few text-based clarifications regarding the interpretation of the locus, listed below. Obviously, if you have any questions feel, please free to reach out.

We thank for all the comments and suggestions the editors and reviewers made to help us improve our manuscript. We marked all the changes we made in the revised manuscript in yellow, and we also included an Appendix key resources table to list all the reagents we used in this study. Note that we also removed the p values and n (sample size) numbers in the main text, as they are already listed in the figures and figure legends. We also corrected some typos that were missed in the original manuscript. Furthermore, to meet the revision submission guidelines, we have included Key Resources Table as an appendix at the end of the manuscript, reformatted previous supplemental tables as Supplementary files 1a and 1b, and added legends to source files (uncropped blot images, and raw data files).

Reviewer #1 (Recommendations for the authors):

It is important to point out that adipocyte differentiation and the expression of the associated metabolic programs critically contribute to the maintenance of insulin sensitivity in obesity. Defects in adipocyte differentiation derived from mitochondrial functions that reduce ATP synthesis and storage of surfeit carbons in fat tissue almost always associate with insulin resistance in animal models and people (see Trends Cell Biol 32:351-364, 2022; Cell. 185:419-446, 2022). The in vitro studies can add to understanding the important roles of MTIF3, but the authors should use more caution when extrapolating to in vivo effects of the eQTL.

We have now incorporated your suggestion, and expanded on it in the manuscript (lines 294-311).

Reviewer #2 (Recommendations for the authors):

1. The authors should give their hypothesis and interpretation of their results in the manuscript. There are no clear connections between results, which makes it difficult to understand the rationale of their experiments.

Thank you for your suggestions. In the revised manuscript, we have included the hypotheses in the introduction (Lines: 70-76) and added or clarified the rationale behind the experiments in several paragraphs of the results (lines: 104-116, 138-140, 159-161, 170-171, 182-183, 199-200, 215-216). We also added more interpretation of our data to the discussion, among the highlighted parts, and the summary statement in lines 321-331.

2. The labels of Figure 1 do not match the figure legends.

We have corrected this error in the revised manuscript in the figure legend.

3. In figure 1F, what is the readout of the relative MTIF3expression? Are they based on western blot? If yes, it would be better to include the blog images.

In Figure 1F, the relative MTIF3 expression represents mRNA level based on RT-qPCR assay. We have clarified that in the revised figure legend (lines 803, 806).

4. In the method, the authors also added extra lipids in the differentiated medium to increase the lipid droplets. I'm curious about the effect of lipid deprivation in KO cells rather than glucose deprivation. The metabolomics showed the change of fatty acids at full medium. There might be some changes in fatty acid utilization in KO cells.

We might not have been clear on this in the original manuscript – we only use extra lipids during the differentiation process. Then, after generating the knockout and during the restriction experiments, the adipocytes are maintained in medium without extra lipids, to ensure more efficient lipolysis. In other words, your suggestion for lipid deprivation has already been met in our experiments. We have clarified that in the methods (lines: 432-439) and in the results (line 192). As you suggest, the fatty acid utilization is affected by the MTIF3 knockout, as demonstrated by the changed metabolite profile (Figure 7) and fatty acid oxidation (Figure 6).

References

1. Rudler, D.L., et al., Fidelity of translation initiation is required for coordinated respiratory complex assembly. Sci Adv, 2019. 5(12): p. eaay2118.

2. Renes, J., et al., Calorie restriction-induced changes in the secretome of human adipocytes, comparison with resveratrol-induced secretome effects. Biochim Biophys Acta, 2014. 1844(9): p. 1511-22.movie

https://doi.org/10.7554/eLife.84168.sa2

Article and author information

Author details

  1. Mi Huang

    Genetic and Molecular Epidemiology Unit, Department of Clinical Sciences, Clinical Research Centre, Lund University, Malmö, Sweden
    Contribution
    Formal analysis, Investigation, Methodology, Writing – original draft
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-6938-5322
  2. Daniel Coral

    Genetic and Molecular Epidemiology Unit, Department of Clinical Sciences, Clinical Research Centre, Lund University, Malmö, Sweden
    Contribution
    Formal analysis, Investigation, Writing – original draft
    Competing interests
    No competing interests declared
  3. Hamidreza Ardalani

    Department of Chemistry, Centre for Analysis and Synthesis, Lund University, Lund, Sweden
    Contribution
    Investigation, Methodology, Writing – original draft
    Competing interests
    No competing interests declared
  4. Peter Spegel

    Department of Chemistry, Centre for Analysis and Synthesis, Lund University, Lund, Sweden
    Contribution
    Formal analysis, Investigation
    Competing interests
    No competing interests declared
  5. Alham Saadat

    Metabolism Program, Broad Institute of MIT and Harvard, Cambridge, United States
    Contribution
    Conceptualization, Resources, Supervision
    Competing interests
    No competing interests declared
  6. Melina Claussnitzer

    Metabolism Program, Broad Institute of MIT and Harvard, Cambridge, United States
    Contribution
    Conceptualization, Resources, Funding acquisition
    Competing interests
    No competing interests declared
  7. Hindrik Mulder

    Unit of Molecular Metabolism, Department of Clinical Sciences, Clinical Research Centre, Lund University, Malmö, Sweden
    Contribution
    Resources, Funding acquisition
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6593-8417
  8. Paul W Franks

    1. Genetic and Molecular Epidemiology Unit, Department of Clinical Sciences, Clinical Research Centre, Lund University, Malmö, Sweden
    2. Department of Nutrition, Harvard T.H. Chan School of Public Health, Boston, United States
    Contribution
    Conceptualization, Supervision, Funding acquisition, Writing – original draft, Project administration
    Contributed equally with
    Sebastian Kalamajski
    For correspondence
    paul.franks@med.lu.se
    Competing interests
    No competing interests declared
  9. Sebastian Kalamajski

    Genetic and Molecular Epidemiology Unit, Department of Clinical Sciences, Clinical Research Centre, Lund University, Malmö, Sweden
    Contribution
    Conceptualization, Formal analysis, Supervision, Funding acquisition, Investigation, Methodology, Writing – original draft, Project administration
    Contributed equally with
    Paul W Franks
    For correspondence
    sebastian.kalamajski@med.lu.se
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6600-9302

Funding

European Research Council (ERC-2015-CoG -681742 NASCENT)

  • Paul W Franks

Vetenskapsrådet

  • Hindrik Mulder

LUDC-IRC

  • Hindrik Mulder

China Scholarship Council (201708420158)

  • Mi Huang

Albert Påhlsson Foundation

  • Sebastian Kalamajski

The Crafoord Foundation (20210571)

  • Sebastian Kalamajski

The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Jennifer Doudna (UC Berkeley) for initial support with concepts relating to using CRISPR in in vitro studies of gene–environment interactions. We also thank Yu-Hua Tseng in Joslin Diabetes Centre for providing the hWAs cells.

This work was supported by China Scholarship Council (201708420158 to Mi Huang); the European Commission (ERC-2015-CoG-681742 NASCENT to PWF), and Swedish Research Council (Distinguish Young Research Reward in Medicine) (to PWF), LUDC-IRC and The Swedish Research Council (to HM), and by The Albert Påhlsson Foundation (to SK).

Senior Editor

  1. David E James, University of Sydney, Australia

Reviewing Editor

  1. Marcus M Seldin, University of California, Irvine, United States

Reviewer

  1. Sean M Hartig, Baylor College of Medicine, United States

Version history

  1. Received: October 12, 2022
  2. Preprint posted: October 17, 2022 (view preprint)
  3. Accepted: March 5, 2023
  4. Accepted Manuscript published: March 6, 2023 (version 1)
  5. Version of Record published: March 17, 2023 (version 2)

Copyright

© 2023, Huang et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 840
    Page views
  • 88
    Downloads
  • 0
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Mi Huang
  2. Daniel Coral
  3. Hamidreza Ardalani
  4. Peter Spegel
  5. Alham Saadat
  6. Melina Claussnitzer
  7. Hindrik Mulder
  8. Paul W Franks
  9. Sebastian Kalamajski
(2023)
Identification of a weight loss-associated causal eQTL in MTIF3 and the effects of MTIF3 deficiency on human adipocyte function
eLife 12:e84168.
https://doi.org/10.7554/eLife.84168

Share this article

https://doi.org/10.7554/eLife.84168

Further reading

    1. Biochemistry and Chemical Biology
    2. Cell Biology
    Kristian Davidsen, Jonathan S Marvin ... Lucas B Sullivan
    Research Article

    Intracellular levels of the amino acid aspartate are responsive to changes in metabolism in mammalian cells and can correspondingly alter cell function, highlighting the need for robust tools to measure aspartate abundance. However, comprehensive understanding of aspartate metabolism has been limited by the throughput, cost, and static nature of the mass spectrometry (MS)-based measurements that are typically employed to measure aspartate levels. To address these issues, we have developed a green fluorescent protein (GFP)-based sensor of aspartate (jAspSnFR3), where the fluorescence intensity corresponds to aspartate concentration. As a purified protein, the sensor has a 20-fold increase in fluorescence upon aspartate saturation, with dose-dependent fluorescence changes covering a physiologically relevant aspartate concentration range and no significant off target binding. Expressed in mammalian cell lines, sensor intensity correlated with aspartate levels measured by MS and could resolve temporal changes in intracellular aspartate from genetic, pharmacological, and nutritional manipulations. These data demonstrate the utility of jAspSnFR3 and highlight the opportunities it provides for temporally resolved and high-throughput applications of variables that affect aspartate levels.

    1. Cell Biology
    2. Immunology and Inflammation
    Chinky Shiu Chen Liu, Tithi Mandal ... Dipyaman Ganguly
    Research Article

    T cells are crucial for efficient antigen-specific immune responses and thus their migration within the body, to inflamed tissues from circulating blood or to secondary lymphoid organs, plays a very critical role. T cell extravasation in inflamed tissues depends on chemotactic cues and interaction between endothelial adhesion molecules and cellular integrins. A migrating T cell is expected to sense diverse external and membrane-intrinsic mechano-physical cues, but molecular mechanisms of such mechanosensing in cell migration are not established. We explored if the professional mechanosensor Piezo1 plays any role during integrin-dependent chemotaxis of human T cells. We found that deficiency of Piezo1 in human T cells interfered with integrin-dependent cellular motility on ICAM-1-coated surface. Piezo1 recruitment at the leading edge of moving T cells is dependent on and follows focal adhesion formation at the leading edge and local increase in membrane tension upon chemokine receptor activation. Piezo1 recruitment and activation, followed by calcium influx and calpain activation, in turn, are crucial for the integrin LFA1 (CD11a/CD18) recruitment at the leading edge of the chemotactic human T cells. Thus, we find that Piezo1 activation in response to local mechanical cues constitutes a membrane-intrinsic component of the ‘outside-in’ signaling in human T cells, migrating in response to chemokines, that mediates integrin recruitment to the leading edge.