Identification of a weight loss-associated causal eQTL in MTIF3 and the effects of MTIF3 deficiency on human adipocyte function
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.sa0Introduction
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).
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.
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).
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).
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).
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).
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 protocolData 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 protocolA 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 protocolGenomic 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 protocolTo 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 protocolhWAs 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 protocolhWAs-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 protocolTo 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 protocolAfter 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 protocolAfter 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 protocolThe 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 protocolTotal 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 protocolCells 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 protocolDifferentiated 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 protocolTo 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 protocolTo 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 protocolThe 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 protocolThe 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 protocolDifferentiated 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 protocolTriglyceride-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
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
-
Overweight, obesity, and mortality in a large prospective cohort of persons 50 to 71 years oldThe New England Journal of Medicine 355:763–778.https://doi.org/10.1056/NEJMoa055643
-
Accurate measurement of circulating mitochondrial DNA content from human blood samples using real-time quantitative PCRMethods in Molecular Biology 1264:117–131.https://doi.org/10.1007/978-1-4939-2257-4_12
-
Mitochondrial function/dysfunction in white adipose tissueExperimental Physiology 99:1168–1178.https://doi.org/10.1113/expphysiol.2014.081414
-
Development and optimization of a metabolomic method for analysis of adherent cell culturesAnalytical Biochemistry 404:30–39.https://doi.org/10.1016/j.ab.2010.04.013
-
“ metabolically healthy obesity ”: origins and implicationsMolecular Aspects of Medicine 34:59–70.https://doi.org/10.1016/j.mam.2012.10.004
-
Understanding adipocyte differentiationPhysiological Reviews 78:783–809.https://doi.org/10.1152/physrev.1998.78.3.783
-
Adipose tissue mitochondrial capacity associates with long-term weight loss successInternational Journal of Obesity 42:817–825.https://doi.org/10.1038/ijo.2017.299
-
Gene × dietary pattern interactions in obesity: analysis of up to 68 317 adults of European ancestryHuman Molecular Genetics 24:4728–4738.https://doi.org/10.1093/hmg/ddv186
-
Cell line models for differentiation: preadipocytes and adipocytesExperimental Biology and Medicine 235:1185–1193.https://doi.org/10.1258/ebm.2010.010063
-
Cell models and their application for studying adipogenic differentiation in relation to obesity: a reviewInternational Journal of Molecular Sciences 17:1040.https://doi.org/10.3390/ijms17071040
-
Protocols for generation of immortalized human brown and white preadipocyte cell linesMethods in Molecular Biology 1566:77–85.https://doi.org/10.1007/978-1-4939-6820-6_8
-
Gene -- environment-wide association studies: emerging approachesNature Reviews. Genetics 11:259–272.https://doi.org/10.1038/nrg2764
-
Mitochondrial regulation and white adipose tissue homeostasisTrends in Cell Biology 32:351–364.https://doi.org/10.1016/j.tcb.2021.10.008
Article and author information
Author details
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).
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
-
- 1,004
- views
-
- 112
- downloads
-
- 0
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cell Biology
Anionic lipid molecules, including phosphatidylinositol-4,5-bisphosphate (PI(4,5)P2), are implicated in the regulation of epidermal growth factor receptor (EGFR). However, the role of the spatiotemporal dynamics of PI(4,5)P2 in the regulation of EGFR activity in living cells is not fully understood, as it is difficult to visualize the local lipid domains around EGFR. Here, we visualized both EGFR and PI(4,5)P2 nanodomains in the plasma membrane of HeLa cells using super-resolution single-molecule microscopy. The EGFR and PI(4,5)P2 nanodomains aggregated before stimulation with epidermal growth factor (EGF) through transient visits of EGFR to the PI(4,5)P2 nanodomains. The degree of coaggregation decreased after EGF stimulation and depended on phospholipase Cγ, the EGFR effector hydrolyzing PI(4,5)P2. Artificial reduction in the PI(4,5)P2 content of the plasma membrane reduced both the dimerization and autophosphorylation of EGFR after stimulation with EGF. Inhibition of PI(4,5)P2 hydrolysis after EGF stimulation decreased phosphorylation of EGFR-Thr654. Thus, EGFR kinase activity and the density of PI(4,5)P2 around EGFR molecules were found to be mutually regulated.
-
- Cell Biology
Classical G-protein-coupled receptor (GPCR) signaling takes place in response to extracellular stimuli and involves receptors and heterotrimeric G proteins located at the plasma membrane. It has recently been established that GPCR signaling can also take place from intracellular membrane compartments, including endosomes that contain internalized receptors and ligands. While the mechanisms of GPCR endocytosis are well understood, it is not clear how well internalized receptors are supplied with G proteins. To address this gap, we use gene editing, confocal microscopy, and bioluminescence resonance energy transfer to study the distribution and trafficking of endogenous G proteins. We show here that constitutive endocytosis is sufficient to supply newly internalized endocytic vesicles with 20–30% of the G protein density found at the plasma membrane. We find that G proteins are present on early, late, and recycling endosomes, are abundant on lysosomes, but are virtually undetectable on the endoplasmic reticulum, mitochondria, and the medial-trans Golgi apparatus. Receptor activation does not change heterotrimer abundance on endosomes. Our findings provide a subcellular map of endogenous G protein distribution, suggest that G proteins may be partially excluded from nascent endocytic vesicles, and are likely to have implications for GPCR signaling from endosomes and other intracellular compartments.