LRG1 is an adipokine that promotes insulin sensitivity and suppresses inflammation
Abstract
While dysregulation of adipocyte endocrine function plays a central role in obesity and its complications, the vast majority of adipokines remain uncharacterized. We employed bio-orthogonal non-canonical amino acid tagging (BONCAT) and mass spectrometry to comprehensively characterize the secretome of murine visceral and subcutaneous white and interscapular brown adip ocytes. Over 600 proteins were identified, the majority of which showed cell type-specific enrichment. We here describe a metabolic role for leucine-rich α–2 glycoprotein 1 (LRG1) as an obesity-regulated adipokine secreted by mature adipocytes. LRG1 overexpression significantly improved glucose homeostasis in diet-induced and genetically obese mice. This was associated with markedly reduced white adipose tissue macrophage accumulation and systemic inflammation. Mechanistically, we found LRG1 binds cytochrome c in circulation to dampen its pro-inflammatory effect. These data support a new role for LRG1 as an insulin sensitizer with therapeutic potential given its immunomodulatory function at the nexus of obesity, inflammation, and associated pathology.
Editor's evaluation
This paper presents a fundamental advance by elucidating the function of LRG1 as an adipokine. The authors provide compelling evidence for the biological effects of LRG1 on metabolism and its potential connections to metabolic diseases.
https://doi.org/10.7554/eLife.81559.sa0Introduction
Obesity is a major threat to human health due to its association with serious comorbidities (Angelantonio et al., 2016; Poirier et al., 2006). Now considered an important endocrine organ, adipose tissue secretes a constellation of bioactive peptides, or adipokines, many of which regulate whole-body energy homeostasis and inflammation (Funcke and Scherer, 2019). Dysregulation of adipose tissue endocrine function is a key feature of obesity and a major contributor to its sequelae. However, the vast majority of adipokines remain unstudied. Addressing this knowledge gap calls for a detailed characterization of the adipose secretome to better understand how adipocytes communicate with other cells to coordinate systemic metabolism.
Mammals possess white (WAT) and brown adipose tissues (BAT), with divergent effects on whole-body metabolism. Visceral (Visc) WAT is particularly associated with obesity-related diseases, whereas subcutaneous (SubQ) WAT is comparatively benign (Fox et al., 2007). On the other hand, active BAT is associated with improved cardiometabolic health (Becher et al., 2021; Lee et al., 2014). While Visc fat contains predominantly white adipocytes that efficiently store energy, SubQ depots contain a mixture of white and thermogenic beige adipocytes which, like brown adipocytes, dissipate energy as heat via uncoupled respiration (Chouchani et al., 2019). In addition to these bioenergetic properties, transplantation studies in mice have demonstrated that secretory mediators likely convey the metabolic effects of different adipose tissues (Tran et al., 2008). Moreover, BAT-specific secreted factors, called batokines, are increasingly appreciated as important contributors to the benefits associated with BAT activation (Villarroya et al., 2017).
Transcriptomic profiling coupled with secretion prediction algorithms is often utilized to identify putative secreted proteins, but this approach could omit non-classically secreted factors. Alternatively, direct detection of proteins in conditioned medium (CM) can be performed using mass spectrometry (MS). Such proteomic analyses are commonly performed in fetal bovine serum (FBS)-free conditions because FBS proteins interfere with detection of relatively low abundance secreted proteins in CM. Serum starvation, however, induces a complex and unpredictable response from cultured cells, making physiological interpretation difficult (Pirkmajer and Chibalin, 2011).
To overcome these limitations, a chemoproteomic technology called BONCAT has been applied for secretome profiling (Dieterich et al., 2006; Eichelbaum et al., 2012). L-Azidohomoalanine (AHA) is a non-toxic, non-canonical amino acid with structural similarity to L-methionine (Met) (Figure 1A). AHA can be recognized by cells’ native translational machinery and incorporated into the nascent proteome in place of Met. The azide functional group in AHA does not exist in nature, so peptides containing AHA can be targeted for bio-orthogonal chemical conjugation strategies called click chemistry. Applying this technique to cultured adipocytes allows for selective enrichment of low abundance nascent secreted proteins for proteomic analysis of serum-containing CM.
Here, we employed BONCAT to comprehensively profile the adipose secretome and nascent serum proteome of mice, from which we identified a novel adipokine, leucine-rich α–2 glycoprotein 1 (LRG1). LRG1 is secreted by mature adipocytes and increased in obesity. We demonstrate that LRG1 improves fasting glucose and insulin tolerance and reduces adipose tissue macrophage accumulation. At the molecular level, LRG1 binds extracellular cytochrome c (Cyt c) released from dead/dying cells and dampens Cyt c’s pro-inflammatory effect on macrophages. These data reveal LRG1 as a crucial regulator of metabolic health as an insulin sensitizer and suppressor of systemic inflammation in obesity.
Results
Characterization of primary adipocyte secretome
To profile the secretome of different types of primary adipocytes, stromal vascular fraction (SVF) from epididymal white adipose tissue (eWAT), inguinal white adipose tissue (iWAT), and interscapular brown adipose tissue (BAT) from C57BL/6 J (B6) mice was differentiated in vitro into primary Visc, SubQ, and Brown adipocytes, respectively (Figure 1B). By day 6 of differentiation, all three cell types showed comparable lipid droplet accumulation and expression of mature adipocyte markers such as Fabp4, Pparg2, and Adipoq (Figure 1—figure supplement 1A, B). In vitro primary adipocytes recapitulated characteristic gene expression patterns in vivo, with expression of Prdm16 and Ucp1 highest in Brown, intermediate in SubQ, and lowest in Visc (Figure 1—figure supplement 1C). Adipocytes were pulsed with Met-deficient media containing 0.1 mM AHA and 10% FBS for 24 hours, after which CM was collected. Except for Pparg2, we did not observe a significant effect of AHA pulse on expression of mature adipocyte or thermogenic genes (Figure 1—figure supplement 1D). To visualize azide-labeled secreted proteins, CM was subjected to a copper-catalyzed azide-alkyne cycloaddition (CuAAC) reaction with tetramethylrhodamine (TAMRA)-alkyne. In-gel fluorescence analysis highlighted adipocyte-derived proteins in CM as bands positive for TAMRA fluorescence (Figure 1C).
For identification and quantitative analysis of the secretome across the three adipose cell types, azide-labeled proteins were enriched from CM by conjugating with alkyne-agarose beads. On-bead digested peptides were subjected to label-free MS. We detected a total of 742 proteins, of which 138 were excluded as reverse hits or contaminants (Figure 1—figure supplement 1E). The intensity-based absolute quantification (iBAQ) value estimates a protein’s molar abundance, and summation of iBAQ intensities (Σ iBAQ) of a sample can estimate total moles of proteins detected (Shin et al., 2013). Σ iBAQ values closely correlated with in-gel fluorescence (Figure 1C), both of which showed highest abundance of secreted proteins in SubQ CM (Figure 1D).
We next examined how many of these AHA-labeled proteins are reported or predicted secreted factors. Referencing gene ontology (GO) cellular component terms for 594 genes from 604 detected proteins, a large majority of proteins (413/594, 70%) were annotated to be extracellular (Figure 1E). We performed retrospective secretion prediction analysis on 604 proteins, evaluating the proportion of proteins predicted to be secreted by the classical ER/Golgi pathway (SignalP5.0, TMHMM2.0), sequence-based deep learning method for subcellular localization prediction (DeepLoc1.0), glycosylphosphatidylinositol (GPI)-anchors (PredGPI), and non-classical routes (SecretomeP2.0) (Almagro Armenteros et al., 2017; Almagro Armenteros et al., 2019; Bendtsen et al., 2004; Krogh et al., 2001; Pierleoni et al., 2008). Over two-thirds of the proteins (413/604, 68%) were predicted to be secreted via the classical pathway, while non-classical secretion constituted about 10% (62/604) (Figure 1F). Overall, 79% of detected proteins (479/604) met at least one criterion of secretion prediction. Notably, more than 20% of identified proteins would not have been predicted to be secreted by these algorithms, highlighting the imperfect nature of in silico predictions and the value of directly measuring secreted protein levels.
In-gel fluorescence analysis showed differential band patterns across the three types of adipocytes, suggesting cell type-specific secretory profiles (Figure 1C). Consistent with this observation, principal component analysis (PCA) showed that each cell type formed a distinct cluster (Figure 1—figure supplement 1F). Comparing all nine AHA-pulsed samples against each other showed higher correlations within biological replicates (R>0.97) than those across different cell types (0.74<R < 0.91) (Figure 1—figure supplement 1G). For quantitative analysis, we focused on 424 proteins detected in at least 2/3 replicates (Figure 1—figure supplement 1E) and compared their Log2-transformed label-free quantification (LFQ) intensities, imputing any missing values with a left-shifted Gaussian distribution. One-way ANOVA yielded 348 proteins with significantly different (q<0.01) secretory profiles across cell types, and unbiased clustering was performed based on protein levels in CM (Figure 1G). Clusters 1 and 4 contained proteins enriched in CM of Visc and Brown, respectively. Cluster 2 proteins were most abundant in CM of SubQ followed by Visc, both white adipocytes. Cluster 3 was highly secreted by SubQ and Brown adipocytes with proteins such as SLIT2 (Svensson et al., 2016), suggesting enrichment of the beige/brown secretome, as well as regulators of immune cell migration and angiogenesis. Cluster 4 comprised the largest number of proteins (125/348), consistent with Brown CM containing the greatest number of unique proteins (Figure 1H). Ranking proteins by decreasing order of abundance, we found high levels of many well-described adipocyte-derived factors such as CFD, RETN, ADIPOQ, and RBP4 (Figure 1I). Interestingly, many cluster 2 proteins ranked highly in abundance, while cluster 4 proteins skewed towards lower abundance (Figure 1I). Hence, while SubQ adipocytes demonstrate high secretory capacity, the brown adipocyte secretome is characterized by a diverse array of proteins, many of which are secreted at lower levels.
We next performed pathway analysis on each of the clusters (Figure 1J). Visc CM-enriched cluster 1 was overrepresented with proteins involved in immune response and complement activation, such as C1QA and CFH (Figure 1K). Proteins with a role in collagen fibril organization and extracellular matrix formation, such as COL1A1 and SFRP, were enriched in cluster 2 (Figure 1L). Cluster 3 (beige/brown enriched) showed overrepresentation of angiogenesis regulators (Figure 1M). Finally, cluster 4 Brown CM-specific proteins were enriched for axon guidance factors such as NRP1 and NGF, consistent with the importance of innervation in thermogenesis (Figure 1K; Chi et al., 2018; Chi et al., 2021; Wang et al., 2020; Zeng et al., 2019).
Profiling of the serum proteome in vivo
Previous studies have shown that AHA can be administered in vivo to label tissue proteins (Calve et al., 2016; McClatchy et al., 2015). However, it has not been tested whether this method can label the nascent serum proteome. We administered 0.1 g/kg/day of AHA IP to chow-fed B6 mice for 2 days (Figure 2A) and did not notice any major adverse effects based on body weight (Figure 2B). CuAAC conjugation of serum proteins with TAMRA-alkyne and in-gel fluorescence analysis showed increased signal across most bands in the AHA group along with some AHA group-specific fluorescent bands (Figure 2C). We used alkyne-agarose beads to enrich the azide-labeled nascent proteome and performed MS analysis. Σ iBAQ showed successful enrichment in animals injected with AHA (Figure 2D). Even without dietary Met restriction or depletion of abundant serum proteins, we were able to identify and quantitate 180 proteins, including classical adipokines such as ADIPOQ, adipsin (CFD), and RBP4 (Figure 2E).
Because AHA can be incorporated into the proteome of any tissue, we employed bioinformatic analyses to compare tissue expression levels of genes that encode serum proteins. We cross-referenced 177 genes from 180 detected proteins with publicly available transcriptomic datasets, such as ENCODE (RNA-Seq based) and BioGPS (microarray based) (Davis et al., 2018; Wu et al., 2016). A t-distributed stochastic neighbor embedding (t-SNE) plot was generated with each protein and its tissue mRNA levels from the ENCODE/LICR dataset to visualize each gene’s basal mRNA expression pattern across multiple tissues from lean adult mice (Figure 2F). The majority of proteins were most highly expressed by the liver, which comprised the largest cluster (98/177, 55%). BAT, the only adipose tissue profiled in the dataset, formed a cluster of 14 genes (7.9%). As expected, classical adipokines known for highly adipose-specific expression such as adiponectin (Adipoq) and adipsin (Cfd) belonged to this group, along with a recently described batokine, Kng2 (Peyrou et al., 2020).
BioGPS offers microarray-based transcriptomic data across a much wider variety of mouse tissue and cell types, including eWAT and BAT. We numerically scaled the degree of adipose tissue enrichment for each detected serum protein. Adipose enrichment of a gene was defined as the number of tissues with expression significantly lower than that of eWAT, BAT, or combined. We divided this value by the total number of pairwise comparisons (i.e., total number of tissues - 2 adipose tissues) to obtain the percentage adipose tissue enrichment score. Among the top-enriched genes were well-described adipokines such as Adipoq, where 76 of 86 tissues (88.4%) expressed Adipoq at significantly lower levels than eWAT and BAT combined (Figure 2G). The remaining 10 tissues that did not meet the statistical significance cutoff still expressed Adipoq at lower levels compared to eWAT or BAT and included adipose-embedded tissues such as mammary glands. We also identified genes yet to be described as adipokines, such as H2-Q10 and Cpn2 (Figure 2G). Many genes, such as Rbp4, Agt, and Lrg1, showed high expression in adipose tissues as well as liver. On the t-SNE plot, these genes were grouped with other liver-specific genes, but located closer to the adipose tissue cluster (Figure 2F). Still, their percentage enrichment scores were >80%, as few other tissues express those genes.
LRG1 is secreted by mature adipocytes and increased in obesity
To identify uncharacterized adipokines with a putative role in whole-body metabolism, we prioritized factors (1) detected in adipocyte CM, (2) enriched in adipose tissues with scores >80%, and (3) present in the nascent serum proteome (Figure 3A). Intersection analysis between (1) and (2) identified 23 proteins (Supplementary file 1), of which 9 proteins were identified in the serum proteome (Figure 3A). ADIPOQ, CFD, and RBP4 have already been identified as adipokines, validating our search strategy. Haptoglobin (HP), ceruloplasmin (CP), and angiotensinogen (AGT) have well-characterized biological functions. We therefore focused on LRG1, a protein with relatively unknown metabolic function.
LRG1, or leucine-rich α–2 glycoprotein 1, has been shown to promote angiogenesis by modulating TGF-β signaling (Wang et al., 2013), and this process has been implicated in various complications of diabetes (Hong et al., 2019; Liu et al., 2020; Zhang et al., 2019). In 3T3-L1 adipocytes, Lrg1 expression correlates with adipogenesis (Birsoy et al., 2011). In BAT, Lrg1 expression is dependent on an EHMT1-driven differentiation program (Ohno et al., 2013). Recently, whole-body LRG1 loss of function has been reported to reduce obesity and improve metabolic health by reducing hepatosteatosis (He et al., 2021). However, the role of adipose-derived LRG1 in whole-body metabolism has not been studied. Consistent with our tissue enrichment analysis, qPCR of mouse tissues showed that Lrg1 mRNA was mostly expressed in adipose tissues and liver (Figure 3B). To determine which cell type(s) within adipose tissue express Lrg1, we fractionated eWAT, iWAT, and BAT to separate floating mature adipocytes from the SVF. We confirmed successful enrichment of Fabp4, Pparg2, and Adipoq in the mature adipocyte fraction and Adgre1 and Pecam1 in the SVF fraction (Figure 3—figure supplement 1A). Lrg1 was successfully co-enriched in the adipocyte fraction in all three depots (Figure 3C). Finally, we cultured primary Visc, SubQ, and Brown adipocytes and confirmed that Lrg1 mRNA is induced >70-fold in all three cell types during in vitro adipogenesis (Figure 3D). Consistent with mRNA data, we found robust levels of LRG1 protein in CM of mature adipocytes, whereas preadipocytes did not secrete detectable LRG1 (Figure 3E, quantified in Figure 3—figure supplement 1B). Of note, LRG1 protein levels detected by western blot were consistent with MS results (Figure 3—figure supplement 1C).
Because Lrg1 expression correlated with adipogenesis, we next sought to determine if Lrg1 directly contributes to adipocyte differentiation. LRG1 has been reported to bind cytochrome c (Cyt c), and this property has been used for depletion of LRG1 by adsorption (Codina et al., 2010; Cummings et al., 2006). We depleted bovine LRG1 in FBS by incubating with Cyt c-conjugated agarose beads and collecting the supernatant. We validated that after 2 rounds of depletion, LRG1 in FBS was no longer detectable, and elution of bound proteins from the Cyt c-agarose beads confirmed specific depletion of LRG1 (Figure 3—figure supplement 2A). SVF from adipose tissues of LRG1-KO mice (generation of which is described further below) was able to differentiate into mature adipocytes in LRG1-depleted media, and the degree of adipogenesis was largely unaffected by adding back recombinant human LRG1 (rhLRG1), as determined by lipid accumulation and mature adipocyte gene expression (Figure 3—figure supplement 2B, C). These results demonstrate that LRG1 does not affect adipogenesis.
To assess whether adipose tissues are a significant contributor to circulating LRG1 levels, we collected serum from chow-fed lean mice or diet-induced obese (DIO) mice on high fat diet (HFD) for 4 or 9 weeks. Serum LRG1 protein levels increased with age and obesity (Figure 3F, quantified in Figure 3—figure supplement 2D). qPCR of major Lrg1-expressing tissues showed significant induction of Lrg1 mRNA in iWAT of DIO mice, but not in liver (Figure 3G). With iWAT expansion in obesity, this induction likely contributes to elevated serum LRG1 levels in this state. Obesity is characterized by chronic low-grade inflammation, with elevated circulating inflammatory cytokines (Lackey and Olefsky, 2016). We observed that treatment of primary SubQ adipocytes with recombinant TNFα induced expression of Lrg1 mRNA (Figure 3H) and protein in (Figure 3I, quantified in Figure 3—figure supplement 2E). Taken together, these results demonstrate that LRG1 is an obesity-induced adipokine.
LRG1 overexpression improves glucose homeostasis in diet-induced obesity
To explore whether LRG1 as an adipokine affects whole-body energy homeostasis, we used viral vectors to overexpress LRG1 in vivo. Adenovirus encoding eGFP (Ad-eGFP) or C-terminally FLAG-tagged LRG1 (Ad-LRG1-FL) was administered to obese B6 mice on HFD for 10 weeks (Figure 4A). Plasma western blot confirmed 44-fold overexpression of LRG1 in the Ad-LRG1-FL group 5 days after infection (Figure 4B, C). We observed no difference in body weights between groups (Figure 4D), but an insulin tolerance test (ITT) showed that the Ad-LRG1-FL group had a significantly enhanced insulin response (p=0.045) (Figure 4E).
To study the longitudinal effects of chronic LRG1 overexpression, we used AAV8, which has tropism for liver and adipose tissue. We administered AAV-eGFP or AAV-LRG1-FL to 6-week-old male B6 mice (Figure 4F) and confirmed 12-fold increased plasma LRG1 in the latter (Figure 4G, H). Tissue qPCR showed significant Lrg1 overexpression in the liver, iWAT, and BAT, with the greatest increase in the liver (Figure 4—figure supplement 1A). At the protein level, however, eWAT and iWAT demonstrated the highest degree of overexpression (Figure 4I, quantified in Figure 4—figure supplement 1B). During 3 months of HFD, both groups gained an equivalent amount of weight to around 50 g (Figure 4J), with no difference in tissue weights at sacrifice (Figure 4—figure supplement 1C). Fasting glucose measurements from 4 to 14 weeks on HFD showed that the AAV-LRG1-FL group had significantly lower fasting glucose levels (p=0.0097) with a dampened peak (LRG1: 225.9±6.5 mg/dL vs. eGFP: 263.7±10.8 mg/dL, mean ± SEM) at week 18 (Figure 4K). Glucose and insulin tolerance tests demonstrated that AAV-LRG1-FL mice had markedly improved glucose (p=0.0029) (Figure 4L, AUC shown in Figure 4—figure supplement 1D) and insulin tolerance (p=0.023) (Figure 4M, AOC shown in Figure 4—figure supplement 1E). To determine if insulin secretion is affected by LRG1, we performed a glucose-stimulated insulin secretion (GSIS) assay. Glucose tolerance was significantly improved in the AAV-LRG1-FL group (p=0.020) (Figure 4—figure supplement 1F) without a significant difference in plasma insulin secretion (Figure 4—figure supplement 1G), and plasma non-esterified fatty acid (NEFA) levels were completely suppressed in both conditions (Figure 4—figure supplement 1G, H). The homeostatic model assessment of insulin resistance (HOMA-IR) is a widely used model to estimate insulin resistance (Matthews et al., 1985). HOMA-IR values trended lower in the LRG1 group, consistent with insulin sensitization (Figure 4—figure supplement 1I). To determine whether hepatic gluconeogenesis contributes to the difference in fasting glucose levels, we performed an oral pyruvate tolerance test (PTT) (Figure 4—figure supplement 1J) and qPCR of hepatic gluconeogenic gene expression (Figure 4—figure supplement 1K) and observed no significant difference between the groups. Taken together, these observations suggest that LRG1 overexpression prevents obesity-related dysregulation of glucose homeostasis by insulin sensitization.
LRG1 loss of function elevates fasting blood glucose in diet-induced obesity
To test whether LRG1 loss of function affects glucose homeostasis, we generated whole-body LRG1-KO animals in B6 background using CRISPR-Cas9 targeting exon 2 of Lrg1 (Figure 4—figure supplement 2A). This led to a frameshift mutation in Lrg1 (Figure 4—figure supplement 2B) and absence of LRG1 protein in plasma (Figure 4—figure supplement 2C). The mouse line was subsequently backcrossed to the parental B6 strain for at least 5 generations to reduce any potential off-target effects. LRG1-KO (Lrg1-/-) and WT (Lrg1+/+) littermate controls on HFD showed no difference in body weight (Figure 4—figure supplement 2D, E) or tissue weights at sacrifice (Figure 4—figure supplement 2F). Interestingly, male LRG1-KO animals demonstrated significantly higher fasting glucose levels compared to WT littermates throughout the HFD challenge (p=0.010) with a higher peak (KO: 261.6±11.3 mg/dL vs. eGFP: 225.0±10.9 mg/dL, mean ± SEM) at week 12 on HFD (Figure 4—figure supplement 2G). Therefore, fasting blood glucose levels are reciprocally regulated by LRG1 gain and loss of function. At the time points tested, however, GTT and ITT did not show significant differences between the genotypes (Figure 4—figure supplement 2H–K).
LRG1 overexpression delays diabetic phenotype and promotes WAT expansion in db/db mice
While B6 mice develop severe obesity upon HFD feeding, they demonstrate only transient and mild hyperglycemia with moderate insulin resistance (Winzell and Ahrén, 2004; Kleinert et al., 2018). Because LRG1 overexpression in B6 animals mitigated hyperglycemia, we explored whether LRG1 can improve glucose homeostasis in C57BLKS/J-Leprdb/db (db/db) mice, a more extreme model of obesity-related type 2 diabetes. Due to a leptin receptor (Lepr) mutation and genetic background, db/db animals demonstrate hyperphagia and early-onset obesity, along with profound hyperglycemia and hyperinsulinemia (Kleinert et al., 2018). We confirmed that obesity and hyperglycemia over 400 mg/dL develop in db/db animals as early as 7 weeks of age, while littermate misty mice with a WT Lepr gene (m/m) maintain fasting glucose levels below 200 mg/dL (Figure 5—figure supplement 1A, B). Similar to B6 DIO mice, db/db mice showed higher circulating LRG1 levels compared to lean littermates (Figure 5—figure supplement 1C, D).
We administered AAV-eGFP or AAV-LRG1-FL to db/db mice at 4 weeks of age, before development of severe hyperglycemia (Figure 5A). In LRG1-overexpressing db/db mice, we confirmed 4.3-fold increased circulating LRG1 (Figure 5B, C) and that the adipose tissues were the major sites of Lrg1 overexpression (Figure 5—figure supplement 1E, F). Starting 2 weeks post-injection, LRG1-overexpressing animals demonstrated accelerated weight gain, such that by 10 weeks of age the LRG1 group weighed 19.6% more (p=0.004) than eGFP controls (Figure 5D). Concomitantly, the AAV-LRG1-FL group showed delayed onset of hyperglycemia. At week 6, we observed frank hyperglycemia in eGFP animals, whereas glucose levels in the LRG1 group were 33.6% lower (eGFP: 399.5±31.7 mg/dL vs. LRG1: 265.4±19.6 mg/dL, mean ± SEM; p=0.0020) (Figure 5E). Fasting plasma insulin levels were almost halved in the LRG1 group (eGFP: 12.5±1.3 ng/mL vs. LRG1: 6.3±0.6 ng/mL, mean ± SEM; p=0.0008) (Figure 5F). HOMA-IR was significantly lower in LRG1-overexpressing mice (Figure 5—figure supplement 1G). Consistent with these findings, ITT profiles significantly differed between the two groups (p=0.019) (Figure 5G), largely due to differences in basal glucose levels (Figure 5—figure supplement 1H). By 8 weeks of age, body weights continued to diverge between the groups, but we no longer observed significant differences in fasting blood glucose or plasma insulin concentrations (Figure 5D–F).
To further characterize the metabolic effects of LRG1 overexpression in db/db animals, we performed indirect calorimetry between 6 and 7 weeks of age (Figure 5A). Measurement of total energy expenditure (EE) across 5 days following acclimation (Figure 5H) and regression analysis against body weight using a generalized linear model using CalR (Mina et al., 2018) showed that there is a significant difference in EE between the two groups (p=0.0010) as well as a significant interaction effect between EE and body weight (p=0.0014) (Figure 5I). This suggests that LRG1 overexpression leads to a greater increase in EE per unit increase in body weight. Respiratory exchange ratio (RER) was higher in the LRG1 group compared to eGFP controls (eGFP: 0.89±0.013 vs. LRG1: 0.93±0.014, mean ± SEM; p=0.046), suggesting a fuel preference for carbohydrates (Figure 5J, K). The eGFP group showed greater locomotor activity than the LRG1 group (eGFP: 234±17.5 vs. LRG1: 185±9.27, mean ± SEM; p=0.027), especially during the dark cycle (Figure 5L, M). Finally, LRG1-overexpressing animals had significantly greater food intake, consuming approximately 14% more food each day (p=0.039) (Figure 5N).
Tissue weight measurements revealed that LRG1-overexpressing animals had accelerated gain of eWAT (p<0.0001) and iWAT (p=0.0001) mass, such that by week 10, their eWAT and iWAT were 53% and 40% heavier, respectively, than those of eGFP controls (Figure 5O, P). We observed no difference in weights of BAT, liver, or gastrocnemius muscle (Figure 5O). To further determine if LRG1 affects the insulin signaling pathway, we measured IRS-1 protein levels in eWAT and liver. The LRG1 group showed increased IRS-1 protein, especially in eWAT (Figure 5—figure supplement 1I and J). qPCR analysis of hepatic gluconeogenic gene expression showed no significant difference except for Pck1, which was higher in the LRG1 group (Figure 5—figure supplement 1K). These results reveal that LRG1 overexpression elicits a constellation of metabolic effects in db/db animals that result in improved insulin sensitivity and WAT expansion-driven weight gain.
LRG1 suppresses obesity-associated systemic inflammation
Based on the accelerated eWAT and iWAT expansion in LRG1-overexpressing db/db animals, we hypothesized that WAT could be a target of LRG1 action. We analyzed paraffin-embedded, hematoxylin and eosin (H&E) stained tissue sections and found that eWAT from B6 animals sacrificed at 21 weeks of age (13 weeks on HFD) was characterized by accumulation of macrophages forming crown-like structures (CLS) (Figure 6A). The distal portion of eWAT was especially susceptible to CLS formation in eGFP controls, while the same region in LRG1-overexpressing animals displayed an 82% reduction in CLS number (p<0.0001) (Figure 6A, B). B6 iWAT contained fewer CLS compared to eWAT and did not show major morphological differences between groups (Figure 6—figure supplement 1A, B). In db/db eWAT, we observed significantly reduced CLS in the LRG1 group by 81% at week 7 (p=0.0001) and 85% at week 10 (p<0.0001) (Figure 6C, D). CLS in db/db iWAT from LRG1-overexpressing animals were reduced by 83% at week 7 (p=0.0043) and 88% at week 10 (p<0.0001) (Figure 6C, E). Quantification of adipocyte cross-sectional areas in db/db eWAT slides showed significantly larger adipocyte size in the LRG1 group, especially at week 7 (p=0.0052) (Figure 6—figure supplement 1C).
Obesity is associated with non-alcoholic fatty liver disease (NAFLD), characterized by hepatic steatosis with or without inflammation (Farrell et al., 2019). While HFD-fed B6 mice rarely demonstrate liver injury or inflammation, mild necroinflammation can be observed in db/db liver as early as 1 month of age (Trak-Smayra et al., 2011). Liver sections from the B6 cohort showed a similar degree of steatosis in both groups, without any inflammatory lesions (Figure 6—figure supplement 1D). In the db/db cohort, both eGFP and LRG1 showed no statistically significant differences in the degree of hepatosteatosis (Figure 6—figure supplement 1E), but inflammatory foci were found only in eGFP-expressing animals (Figure 6C).
Inflammation is a key link between obesity and insulin resistance (Saltiel and Olefsky, 2017). We hypothesized that LRG1 may mediate its insulin sensitizing effect via attenuation of inflammation in susceptible organs. We performed RNA-Seq analysis on eWAT from AAV-treated db/db animals harvested at week 7 (midpoint) and 10 (endpoint). We performed differential gene expression analysis between the eGFP and LRG1 groups at midpoint and subjected the list of significant genes (p<0.01) to GO gene-set enrichment analysis (GSEA). All of the top 20 differentially regulated pathways showed highly significant enrichment (p=0.0067) and were immune-related, including leukocyte activation, innate immune response, and cytokine production (Figure 6F). The enrichment score for each of these pathways was negative in the LRG1 group, indicating down-regulation of inflammatory processes in these animals. GSEA of significantly regulated genes at the 10-week endpoint yielded similar down-regulation of immune-related pathways in the LRG1 group (Figure 6—figure supplement 1F). To visualize which genes are most significantly differentially regulated, we plotted a heatmap of 68 genes that showed significant difference (adjusted P<0.01) with Log2 fold-change of >2 or <-2 between the two groups at both time points. Consistent with the pathway analysis, the LRG1 group showed significant down-regulation of cytokines and chemokines such as Ccl3, Ccl4, and Cxcl9; metalloproteinases such as Mmp12 and Mmp13 known to be highly expressed by macrophages; and various immunoglobulin subunit genes (Figure 6G).
Our histological analysis demonstrated significantly reduced CLS in eWAT of LRG1-overexpressing animals. Deconvolution algorithms such as CIBERSORTx allow estimation of cell populations from bulk RNA-Seq datasets (Newman et al., 2019). We performed CIBERSORTx analysis on the expression dataset to gain further insight into differences in immune cell populations. CIBERSORTx estimated that the LRG1 group contains fewer immune cells (Figure 6H). Macrophages were predicted to constitute the majority of immune cells in every eWAT sample analyzed, and the LRG1 group demonstrated a lower absolute quantity of macrophages (Figure 6H), without affecting their relative proportions (Figure 6—figure supplement 1G). Many of the differentially regulated genes between eGFP- and LRG1-overexpressing animals encode chemokines and cytokines. We performed a multiplex cytokine assay to assess whether these differences are reflected in serum levels, and at midpoint (week 7), the LRG1 group showed 28–86% reduction of circulating chemokines such as MCP-1, MIP-1α, MIG, and IP-10 and cytokines such as TNFα and IFNγ (Figure 6I). Many of these differences subsided by week 10, mostly due to a reduction of cytokine levels in the eGFP group (Figure 6—figure supplement 1H). Taken together, these results suggest that LRG1 overexpression attenuated pro-inflammatory processes associated with obesity.
LRG1 binds extracellular cytochrome c and blocks its pro-inflammatory effect on macrophages
LRG1 contains leucine-rich repeat (LRR) domains, which form a structural framework for protein-protein interactions. We hypothesized that LRG1’s immunomodulatory function is mediated by protein-protein interactions. LRG1 has been reported to bind cytochrome c (Cyt c), a mitochondrial protein (Cummings et al., 2006). In addition to its role in the respiratory chain and intrinsic apoptosis pathway, Cyt c is released into the extracellular space following cell death (Jemmerson et al., 2002; Renz et al., 2001) and mediates pro-inflammatory signals as a damage-associated molecular pattern (DAMP) (Grazioli and Pugin, 2018). Binding between LRG1 and extracellular Cyt c has been shown to reduce apoptosis of lymphocytes and cancer cells in vitro (Codina et al., 2010; Jemmerson et al., 2021). However, the in vivo biological context and effect of this interaction have remained unknown. Adipocyte death in obesity is a key event promoting macrophage infiltration and WAT inflammation (Cinti et al., 2005), but the exact triggers of metabolic inflammation remain unidentified. We explored whether extracellular Cyt c released by dead/dying adipocytes could be a key mediator of macrophage recruitment/activation, and if LRG1 in turn modulates the pro-inflammatory action of Cyt c.
We first examined whether LRG1 binds extracellular Cyt c in the circulation. We utilized the C-terminal FLAG-tag of overexpressed LRG1 and α-FLAG antibody-conjugated beads to co-immunoprecipitate (IP) LRG1-FL and interacting partners. In the serum of db/db animals transduced with AAV-LRG1-FL, we were able to successfully co-IP both LRG1-FL and Cyt c from the serum, indicating the two proteins indeed circulate as a complex (Figure 7A). While elevated Cyt c levels in circulation have been observed in conditions involving cell death or inflammation (Adachi et al., 2004; Alleyne et al., 2001; Barczyk et al., 2005; Ben-Ari et al., 2003), it is not known whether obesity is associated with increased serum Cyt c levels. Western blot analysis of serum revealed that Cyt c is increased in db/db animals compared to lean littermates (Figure 7B, quantified in Figure 7—figure supplement 1A). Interestingly, extracellular Cyt c levels trended higher at 7 than at 10 weeks of age (Figure 7B, quantified in Figure 7—figure supplement 1A), correlating with the time point when LRG1 overexpression has a potent physiological effect. In B6 mice, HFD led to increased serum Cyt c levels (Figure 7C, quantified in Figure 7—figure supplement 1B), and the trend continued at 15 weeks on HFD (Figure 7—figure supplement 1C, D). To test if dying adipocytes contribute to circulating Cyt c, we subjected primary SubQ adipocytes to conditions that induce apoptosis or necrosis. In CM of SubQ adipocytes treated with staurosporine (STS) for 24 hr, we detected a robust increase in Cyt c compared to vehicle controls (Figure 7D, quantified in Figure 7—figure supplement 1E). CM collected from adipocytes incubated in hypoxic conditions or treated with TNFα also showed increased Cyt c compared to normoxic or vehicle controls, respectively (Figure 7D, quantified in Figure 7—figure supplement 1E). We next collected interstitial fluid (ISF) from eWAT and iWAT of B6 mice and found higher Cyt c levels in HFD-fed mice (Figure 7E, quantified in Figure 7—figure supplement 1F), suggesting adipose tissues are a source of elevated circulating Cyt c in DIO mice.
Extracellular Cyt c has been shown to exert a pro-inflammatory effect by acting on the toll-like receptor 4 (TLR4)-mediated innate immune signaling pathway in astrocytes (Wenzel et al., 2019). Because circulating monocytes contribute to adipose tissue macrophages in obesity (Oh et al., 2012), we tested whether bone-marrow derived macrophages (BMDMs) respond to extracellular Cyt c by upregulating pro-inflammatory genes. Treating BMDMs with horse Cyt c led to >30-fold induction of lipopolysaccharide (LPS)-responsive pro-inflammatory genes such as Il1b, Cxcl10, Il6, and Nos2 (p<0.0001), while pre-treatment of macrophages with a small molecule TLR4 inhibitor, TAK-242, prevented this induction (Figure 7F). These data confirm that Cyt c in the extracellular space acts as a DAMP that activates innate immune signaling and polarizes macrophages into a more pro-inflammatory state. We then explored how co-treatment of Cyt c and recombinant human LRG1 (rhLRG1) affects Cyt c’s pro-inflammatory effect on macrophages. Prior to treatment, media containing Cyt c, rhLRG1, or both were incubated for one hour to allow protein-protein interactions to occur. Compared to single treatment, co-treatment with rhLRG1 significantly dampened Cyt c-mediated induction of Il1b (P<0.0001), Cxcl10 (P=0.024), and Il6 (P=0.0008), while Nos2 induction was not affected by co-treatment with rhLRG1 (Figure 7G). Interestingly, pre-treatment of BMDMs with rhLRG1 followed by addition of Cyt c did not attenuate the latter’s pro-inflammatory effect, suggesting that Cyt c-LRG1 complex formation is necessary for LRG1’s modulatory effect (Figure 7—figure supplement 1G). Expression of alternatively activated (M2) macrophage markers was not affected by co-treatment with Cyt c and rhLRG1 (Figure 7—figure supplement 1H). We next explored if this immunomodulatory role of LRG1 can in turn reduce adipose tissue damage. Western blot of db/db eWAT showed reduction of cleaved Caspase-3 in the LRG1 group (Figure 7H), suggesting reduction of apoptotic cells with LRG1 overexpression. We also observed decreased release of Cyt c into the ISF of eWAT and liver (Figure 7I), consistent with lower circulating Cyt c detected in the serum (Figure 7A). Taken together, these results suggest LRG1 can modulate the pro-inflammatory gene expression responses induced by Cyt c and subsequently reduce tissue damage and cell death.
Cyt c is a small (12.5 kDa), globular protein with a short half-life in circulation (Radhakrishnan et al., 2007). Small proteins such as Cyt c are expected to be rapidly cleared via glomerular filtration in the kidney. With the tight binding interaction between LRG1 and Cyt c, we hypothesized that much of extracellular Cyt c circulates as a complex with LRG1. To assess whether this interaction affects Cyt c clearance, we performed retroorbital injection of horse Cyt c into WT and LRG1-KO mice and collected blood at time points from 5 min to 6 hr post-injection. Plasma western blot showed that LRG1-KO animals almost completely cleared Cyt c by 2 hr, whereas in WT animals, an initial decrease in Cyt c was followed by nearly constant levels between 1 and 6 hr post-injection (Figure 7—figure supplement 1I). Fitting one-phase decay functions to relative Cyt c intensities showed that while both WT and LRG1-KO curves have similar half-lives (WT: 16.10 min vs. KO: 17.32 min), the WT trace plateaus at 0.2791, whereas the LRG1-KO curve approaches zero (Figure 7—figure supplement 1J, K). Concomitantly, appearance of Cyt c in the urine was observed as soon as 30 min after the injection, demonstrating that clearance of excess Cyt c is indeed through the kidneys (Figure 7—figure supplement 1L, M). Urine Ponceau S staining showed that detected proteins were all below 25 kDa, the apparent size limit for urinary excretion (Figure 7—figure supplement 1L). We reasoned that in LRG1-KO animals, Cyt c circulates as an unbound, free form that can be rapidly excreted by glomerular filtration, whereas the presence of LRG1 in WT animals prevents a portion of Cyt c from excretion due to formation of a complex larger than the size cutoff for glomerular filtration. Indeed, endogenous Cyt c was virtually undetectable in the serum of HFD-fed LRG1-KO mice (Figure 7—figure supplement 1N). Taken together, these data directly demonstrate that Cyt c is a native ligand of LRG1, and that this interaction suppresses the pro-inflammatory effect of extracellular Cyt c and cell death (Figure 7J).
Discussion
We applied a chemoproteomic technology to comprehensively profile the secretome of three major types of murine primary adipocytes, revealing unique secretory profiles of each cell type with differential functional enrichment. This method was also applied in vivo to profile the nascent serum proteome, and bioinformatic analysis suggested that adipose tissue is an important contributor to the serum proteome. The intersection of the adipocyte CM and serum proteome included classical adipokines such as adiponectin (Wang and Scherer, 2016), adipsin (Lo et al., 2014; Gómez-Banoy et al., 2019), and RBP4 (Yang et al., 2005), all of which have important roles in whole-body energy homeostasis. Leptin was not detected in our dataset, as has been the case in most other MS-based approaches (Ali Khan et al., 2018). This is most likely due to its very low expression in in vitro-differentiated adipocytes (MacDougald et al., 1995; O’Connor et al., 2021). LRG1 shares a similar expression signature with these adipokines, but its metabolic function has not been characterized. We demonstrate, using two types of viral vectors and two different mouse models of obesity, a novel metabolic role for adipose-derived LRG1 as a regulator of glucose homeostasis by promoting insulin sensitization. This insulin-sensitizing effect in LRG1 gain of function is associated with a dramatic reduction in systemic inflammation and with LRG1’s ability to bind Cyt c and modulate its pro-inflammatory effect as a DAMP.
There have been various other studies that used MS to profile the secretome of adipocytes (Ali Khan et al., 2018; Deshmukh et al., 2019; Lehr et al., 2012; Roca-Rivada et al., 2015; Rosenow et al., 2010; Zhou et al., 2009). Comparison of our dataset with those from two most recent publications (Ali Khan et al., 2018; Deshmukh et al., 2019) highlights key advances with our approach as well as potential caveats (Supplementary file 2). Our dataset stands out for its use of all three major types of adipocytes including Visc cells, which have important implications for cardiometabolic complications of obesity (Fox et al., 2007). While the total number of secreted proteins in our dataset may be lower than the other two studies (604 vs 1337 and 1866), the number of known secreted proteins is higher (472 vs 453 and 471). The low proportion of intracellular proteins detected in our dataset demonstrates the advantage BONCAT provides by enabling FBS supplementation, which is essential for cell viability. It must be noted, however, the undefined nature and batch variability of FBS can make the study less generalizable, and that FBS contamination is still present, in up to 138/742 of proteins (18.6%). In addition, there could be a potential bias for Met-rich peptide sequences. Finally, although AHA pulse did not affect expression of adipocyte and thermogenic genes tested, its full biological effect is unknown and may include protein misfolding and deficiency in single-carbon metabolism. The remaining 20% of proteins not predicted or annotated to be secreted include intracellular proteins such as tubulin and actin. However, many are non-classically secreted proteins such as FABP4/5 (Hotamisligil and Bernlohr, 2015). Additionally, many proteins (such as GAPDH) that were initially thought to be leakage products turn out to be secreted in a regulated fashion (Takenouchi et al., 2015), and our dataset includes other novel factors we have independently confirmed as bona fide secreted factors. Our findings can therefore greatly expand the scope of the adipose secretome that includes non-classical modes of secretion, but each factor will need to be carefully studied to validate whether its presence in the media is from a regulated secretion process or from non-specific leakage.
We have also demonstrated that BONCAT can label nascent serum proteins in vivo. This technique can be further utilized to study the nascent serum proteome in response to stimuli difficult to model in vitro, such as obesity and changes in ambient temperature. Novel technologies for in vivo profiling of cell type-specific proteomes using non-canonical amino acids (Alvarez-Castelao et al., 2017) or proximity labeling (Branon et al., 2018; Wei et al., 2021) are on the horizon. This study provides a baseline profile of the nascent serum proteome in lean adult mice from which cell type enrichment analyses can be performed, and bioinformatic pipelines implemented here to compare tissue expression would be useful in cross-referencing with transcriptomic datasets.
AAV-mediated LRG1 overexpression in the adipose tissues led to a constellation of metabolic effects in both B6 and db/db mice, resulting in improvement of glucose/insulin tolerance as well as a reduction of obesity-associated inflammation. While endogenous LRG1 levels increase with obesity, further induction of LRG1 expression via AAV transduction protected mice from obesity-related complications. These findings suggest that while LRG1 plays a protective role in obesity, the degree of LRG1 induction is not sufficient, and that AAV-mediated overexpression could be augmenting this physiologic response. LRG1 overexpression had a more dramatic effect in the db/db model, likely due to the severity of inflammation and metabolic derangements. In both mouse models, LRG1 overexpression led to an improvement in fasting glucose levels, which also contributed to differences in GTT and ITT profiles. We believe it is the improvement in insulin sensitivity that drives both findings, given increased IRS-1 protein levels in eWAT and lack of difference in hepatic gluconeogenic gene expression or pyruvate tolerance.
Adipose tissue expansion in mice is characterized by periods of adipocyte turnover, during which there is a sharp increase in adipocyte death and monocyte/macrophage accumulation (Rosen and Spiegelman, 2014; Weisberg et al., 2003). In B6 eWAT, adipocyte death and CLS formation peak around 12–16 weeks on HFD (Strissel et al., 2007), which corresponds to when a parallel surge in fasting blood glucose levels is observed. In db/db mice, the remodeling seems to be initiated earlier, as CLS in WAT emerge at 7 weeks of age, and chemokine/cytokine expression and circulating levels are more highly induced at week 7 than week 10. Periods when LRG1’s insulin-sensitizing and anti-inflammatory effects are particularly effective coincide with these times of adipocyte death and turnover. In db/db mice, we also observed accelerated WAT expansion with LRG1 overexpression. Adipose tissue hypertrophy likely drives this process, given larger adipocyte cross-sectional areas and LRG1’s lack of adipogenic effect. As inflammation is a key mechanistic link between obesity and insulin resistance (Saltiel and Olefsky, 2017), these findings strongly suggest LRG1 promotes insulin sensitization via modulation of inflammation and cell death during adipose tissue remodeling.
Recently, whole-body LRG1 loss of function has been reported to reduce obesity and improve insulin sensitivity by reduction of hepatosteatosis (He et al., 2021). In the LRG1-KO mice we generated, we have not observed any effect on body weight or glucose/insulin tolerance profiles. This discrepancy is potentially due to differences in how the knockout lines were generated (gene cassette vs. CRISPR-Cas9). The work presented here sheds light on the biology of adipose-derived LRG1 with AAV overexpression, which has not been studied before. Our work also expands our understanding of LRG1 function in its native tissue environment, the adipose tissue, where it reduces inflammation and cell death.
Our studies build on our understanding of the LRG1-Cyt c interaction that has thus far been studied in cell culture models (Codina et al., 2010; Jemmerson et al., 2021). We propose LRG1 functions as a buffer against deleterious effects of extracellular Cyt c released from dying/dead cells. During obesity-driven adipose tissue remodeling, increase in adipocyte death contributes to increased circulating Cyt c, which could be an important pro-inflammatory signal triggering monocyte/macrophage recruitment. The pro-inflammatory setting of obesity also induces LRG1 in other adipocytes, and this increase in circulating LRG1 modulates inflammation by directly binding Cyt c. Of note, both TLR4 and LRG1 contain LRR domains, so the Cyt c-neutralizing effect by LRG1 could involve steric interference as the two proteins compete for binding Cyt c. It is interesting to note that LRG1 overexpression led to reduction in circulating Cyt c levels. As our data suggest LRG1-Cyt c complex formation inhibits renal excretion of Cyt c, the reduction of serum Cyt c in the context of LRG1 overexpression is most consistent with decreased Cyt c release from tissues. Together with our data showing reduction of Cyt c in the ISF of eWAT and liver from LRG1-overexpressing mice, LRG1’s immunomodulatory role in these metabolic tissues may contribute to decreased Cyt c release by reducing tissue damage. Many questions remain on how LRG1 and Cyt c regulate inflammation and cell death. Whether LRG1-Cyt c interaction affects tissue-resident macrophages or also impacts immune cell recruitment remain to be studied. Mechanistically, as extracellular Cyt c is pro-apoptotic, it is possible that LRG1-Cyt c binding directly prevents apoptosis either by steric hindrance or by simultaneous binding of LRG1 to Cyt c and TGF-β1 (Jemmerson, 2021). The scope of LRG1’s metabolic function beyond its interaction with Cyt c will also continue to be explored in our future work.
While the current study focused on LRG1’s effect on adipose tissue, it is important to emphasize that LRG1-mediated anti-inflammatory effects were systemic, as observed by the decrease in inflammatory foci in LRG1-overexpressing db/db liver. As hepatocyte apoptosis contributes to progression of NASH (non-alcoholic steatohepatitis) in NAFLD (Feldstein et al., 2003), further investigation may uncover a hepato-protective role for LRG1 in NAFLD. Multiple biomarker studies in humans have shown that increased LRG1 is associated with not only obesity (Pek et al., 2018) but also conditions involving cell death and inflammation, such as appendicitis (Kentsis et al., 2010), solid tumors (Belczacka et al., 2019) and autoimmune diseases (Naka and Fujimoto, 2018). Our study suggests that induction of LRG1 in these settings could be a compensatory mechanism to modulate inflammation.
Enhanced renal clearance of Cyt c in LRG1-KO mice provides further evidence that LRG1 complexes with Cyt c in circulation and highlights the need to measure LRG1 levels in future studies utilizing extracellular Cyt c as a biomarker. Comparable serum Cyt c half-lives between WT and LRG1-KO suggest LRG1 loss of function does not affect kidney function. Divergence in Cyt c clearance profiles occurs towards the plateau, and we believe this difference represents the portion of Cyt c that stays bound to LRG1 in WT animals. Inflammation in obesity, especially in B6 DIO mice, is chronic and relatively low grade. Given our observation that HFD-fed LRG1-KO mice have reduced Cyt c levels in circulation, enhanced clearance of Cyt c in LRG1-KO mice could serve to counteract Cyt c release and protect from its deleterious effects. It is possible therefore that the effect of LRG1 loss of function may be more pronounced in conditions involving acute or high-grade inflammation. It is interesting to note that in streptozotocin (STZ)-induced diabetic mouse models, LRG1 loss of function had a protective role in development of diabetic nephropathy (Hong et al., 2019). In addition to the reported LRG1-mediated aberrant glomerular angiogenesis as a mechanism for this finding, glomerular deposition of LRG1-Cyt c and subsequent immune response could be another potential pathogenic mechanism.
In summary, we identified LRG1 as a novel adipokine and discovered its metabolic role as an insulin sensitizer and suppressor of inflammation. LRG1’s in vivo function is associated with Cyt c binding, and this interaction sheds light on previously unappreciated molecular players of adipocyte-macrophage interaction during adipose tissue expansion and remodeling. LRG1 or targeting extracellular Cyt c could be an attractive therapeutic approach for treatment of not only obesity but a variety of inflammatory conditions.
Methods
Animals
C57BL/6 J mice were purchased from Jackson Laboratories. C57BLKS-Leprdb homozygote males and C57BLKS-Dock7m homozygote males were purchased from Jackson Laboratories. LRG1-KO (Lrg1-/-) line was generated in C57BL/6 J background by the CRISPR and Gene Editing Center at the Rockefeller University and backcrossed with C57BL/6 J mice for at least five generations to minimize off-target effects. Cohorts of LRG1-KO and WT littermates for characterization studies were obtained by intercrossing male and female heterozygotes. Animals were maintained at the Rockefeller University Comparative Biosciences Center, housed at 23 °C and maintained at 12 hr light:dark cycles. The animals were group-housed with ad libitum access to food and water except during metabolic characterization studies. C57BL/6 J mice were fed standard chow diet (LabDiet 5053), and where specified, started on 60% high fat diet (HFD, Research Diets D12492) feeding at 6 or 8 weeks of age. C57BLKS-Leprdb and C57BLKS-Dock7m were fed standard chow diet. For virus injections, each cage containing 4 mice was randomly assigned to eGFP or LRG1 group. Injections were performed in designated ABSL-2 housing rooms and the transduced mice were quarantined for 72 hr post-injection before transferred to regular housing rooms. Experiments involving adenoviral and AAV8 vectors were performed in accordance with the institutional ABSL-2 guidelines. Except in the case of mortality, no mice were excluded from study analysis. AAV experiments consisted of 2 B6 cohorts and 3 db/db cohorts. All animal studies were performed in accordance with the institutional guidelines of the Rockefeller University Institutional Animal Care and Use Committee (IACUC).
Metabolic characterization of mice
For studies involving diet-induced obesity, mouse body weights were monitored once a week, providing fresh 60% high fat diet (Research Diets) at least once a week. For fasting blood glucose measurements, mice were single housed in the morning in cages with fresh bedding and access to water but without food. Mice were kept in a procedure room free of noise or vibration throughout the experiment. After 6 hr, blood was collected from the tail vein, and glucose levels were measured using a glucose meter. For plasma insulin ELISA, blood was also collected in EDTA-coated capillary tubes, which were centrifuged at 2000×g for 15 min at 4 °C to collect plasma. Insulin ELISA was performed using Ultra Sensitive Mouse Insulin ELISA Kit (Crystal Chem). Glucose (GTT), pyruvate (PTT) and insulin tolerance tests (ITT) were performed following a similar 6 hr fasting procedure and started by administering at time 0 indicated doses of glucose or Novolin R human insulin (Novo Nordisk) by intraperitoneal injection or sodium pyruvate by oral gavage. Following injection/gavage, blood glucose measurements were taken from the tail vein at indicated timepoints. During an ITT procedure, mice with glucose measurements below 20 mg/dL or showing signs of hypoglycemia were rescued by 1 g/kg glucose IP injection and excluded from the study. For glucose-stimulated insulin secretion (GSIS) assay, mice were fasted for 6 hr, followed by intraperitoneal injection of glucose and subsequent blood glucose measurement and blood collection at indicated intervals. Plasma insulin levels were measured as described above. Plasma non-esterified fatty acid (NEFA) concentrations were measured using the HR Series NEFA-HR(2) assay (Wako/Fujifilm).
Indirect calorimetry
Phenomaster automated home cage phenotyping system (TSE Systems, Bad Homburg, Germany) was used for the study. Mice were housed individually in controlled environmental chambers maintained at 22 °C, 12 hr light:dark cycles, and 40% humidity. Mice were given ad libitum access to food and water, and bedding was changed every 5 days. For O2 and CO2 measurements, gas was sampled at 15-min intervals with a settle time of 3 min, at 0.25 L/min sample flow rate. Raw values were analyzed using CalR (Mina et al., 2018). Daily food intake was measured directly by weighing remaining food pellets in each cage.
Generation of LRG1-KO mice
CRISPR guide RNAs were designed using CRISPOR.org (Concordet and Haeussler, 2018) and were used as two-part synthetic crRNA and tracrRNA (Alt-RTM CRISPR guide RNA, Integrated DNA Technologies, Inc). Cas9 protein, crRNA, and tracrRNA were assembled to ctRNP using protocols described previously (Shola et al., 2021). Two crRNAs were assembled to ctRNPs and electroporated to one-cell-stage mouse embryos to assess their efficiency in generating indels on the Exon 2 of Lrg1 gene. To prepare for the microinjection mix, crRNA-B which binds to genomic target sequence “AATCTCGGTGGGACCATGGCAGG” was selected for its high on-target efficiency and low off-target potential. The final injection mix was made of 0.6 µM of guide RNA (crRNA +tracrRNA) and 0.3 µM of Cas9 protein according to protocols described previously (Shola et al., 2021). The injection mix was then delivered to 0.5 days of fertilized C57BL/6 J mouse embryos using well-established pronuclear injection and surgical protocols (Shola et al., 2021).
Cells
Mouse primary stromal vascular fraction (SVF) cells were obtained from adipose tissues of 6- to 8-week-old male mice by collagenase digestion and plated on collagen I-coated dishes. SVF cells from epididymal white adipose tissue (eWAT) were grown in ITS media containing 1.5:1 mixture of low-glucose DMEM:MCDB201 supplemented with 2% FBS (Gemini), 1% ITS premix (Corning), 0.1 mM L-ascorbic acid 2-phosphate (Sigma), 10 ng/mL bFGF (Thermo), 0.5% penicillin/streptomycin (P/S, Gibco), and 0.2% primocin (InvivoGen). SVF cells from inguinal white adipose tissue (iWAT) and interscapular brown adipose tissue (BAT) were grown in DMEM/F-12 GlutaMAX medium (Gibco) containing 10% FBS and 1% P/S. Once grown to confluence, differentiation and maintenance of primary adipocytes across cell types were done using DMEM/F-12 GlutaMAX medium containing 10% FBS and 1% P/S. eWAT and iWAT SVF cells were induced to differentiate with an adipogenic cocktail (0.5 mM IBMX, 1 µM dexamethasone, 1 µM rosiglitazone, and 850 nM insulin) for the first 2 days, followed by 2 days of 1 µM rosiglitazone and 850 nM insulin, after which the cells were maintained in 850 nM insulin for additional 2–4 days. SVF cells from BAT were differentiated as above but with 17 nM insulin. Experiments with primary adipocytes were performed between days 6 and 8 of differentiation. All cultured primary adipocytes were checked for lipid accumulation under a phase-contrast microscope before studies. Primary SVF and adipocytes were maintained at 37 °C with 10% CO2.
Bone marrow-derived macrophages (BMDMs) were obtained from 8- to 10-week-old males. Femurs and tibias were dissected, cleaned, and sterilized with ethanol before flushed of bone marrow cells, which were plated onto petri dishes. Bone marrow cells were differentiated in RPMI-1640 medium (Gibco) supplemented with 20% heat-inactivated FBS (Sigma), 1% P/S, and 100 ng/mL M-CSF (Biolegend) for 6–7 days, changing media every 2–3 days. Differentiated BMDMs were washed, trypsinized, and plated onto TC-treated culture plates for overnight before studies were performed. Experiments with BMDMs were performed between days 6 and 7 of differentiation, in X-VIVO 15 medium (Lonza). BMDMs were differentiated and maintained at 37 °C with 5% CO2.
HEK293A cells were purchased from Invitrogen and grown using 4.5 g/L glucose DMEM (Gibco) supplemented with 10% FBS and 1% P/S at 37 °C with 5% CO2. Cells with passage number under 20 were used for adenovirus production. Cells were validated to be mycoplasma free.
Oil Red O
Cells were washed with PBS and fixed with 10% formalin in PBS. Cells were washed with distilled water and 100% propylene glycol. 0.5% Oil Red O staining solution in propylene glycol was added and cells and incubated overnight at room temperature. Following staining, cells were washed with 85% propylene glycol solution followed by washes with distilled water. Stained cells were stored in 50% of glycerol.
Proteomic analysis using BONCAT
Conditioned medium generation
On day 6 of differentiation, primary adipocytes on collagen-coated 6-well plates were washed twice with warm PBS and pulsed with 1 mL/well of Met-free DMEM containing 10% dialyzed FBS, 1% P/S, 17 nM or 850 nM insulin, and either 0.1 mM AHA or 0.1 mM Met. Following 24 hr incubation at 37 °C with 10% CO2, conditioned media (CM) from six wells (1 plate) were collected and pooled, filtered through a 0.22 µm PES membrane syringe filter unit, and supplemented with ½ tabs of EDTA-free cOmplete mini protease inhibitor cocktail (Roche) and PhosSTOP (Roche). CM was concentrated using a 3 kDa centrifugal filter unit (Millipore).
AHA administration and serum collection
Mice were injected with 0.1 g/kg/day AHA or PBS IP for 2 consecutive days and sacrificed 24 hr following the second injection. Following decapitation, truncal blood was collected, allowed to clot for 15 min at room temperature, and centrifuged at 2000×g for 15 min at 4 °C to collect serum.
In-gel fluorescence analysis
Concentrated CM or serum was dialyzed with phosphate-buffered RIPA (10 mM phosphate buffer pH 7.2, 1% Triton X-100, 0.1% Na deoxycholate, 0.1% SDS, 140 mM NaCl) supplemented with EDTA-free cOmplete mini protease inhibitor cocktail (Roche) and PhosSTOP (Roche) using a 3 kDa centrifugal filter unit (Millipore) and protein concentration was determined using Pierce BCA Protein Assay Kit (Thermo Scientific) using a dilution series of bovine serum albumin as protein standards. Copper(I)-catalyzed azide-alkyne cycloaddition reaction with TAMRA-alkyne (Invitrogen) was performed by mixing 200 µg of CM proteins with 0.1 mM TAMRA-alkyne, 1 mM TCEP, 0.1 mM TBTA, and 1 mM of CuSO4 in phosphate-buffered RIPA and rotated end-over-end for 1 hr at room temperature under protection from light. Following methanol/chloroform precipitation, the dried protein pellet was dissolved in Laemmli loading buffer. Following polyacrylamide gel electrophoresis, the gel was briefly washed with distilled H2O and imaged with a Typhoon 5400 imager (GE Healthcare) using 532 nm excitation and a 580 nm detection filter.
Enrichment of labeled proteins
Azide-labeled nascent protein in the concentrated CM was enriched using Click-iT Protein Enrichment Kit (Invitrogen). Serum was diluted 1:1 with the lysis buffer provided with the kit and subjected to enrichment. Enrichment and resin wash was performed following the protocol from Eichelbaum and Krijgsveld, 2014.
On-bead digestion
Extensively washed beads were incubated with Lys-C endopeptidase (Wako) in 4 M urea and 0.14 M NH4HCO3 by shaking at 1400 rpm for 6 hr at room temperature. The resin mixture was further digested by adding trypsin (Promega) in 2 M urea and 0.14 M NH4HCO3 and incubated by shaking at 1400 rpm overnight at room temperature. The following day, the digestion reaction was quenched by adding trifluoroacetic acid.
LC-MS/MS
Tryptic peptides were desalted (Rappsilber et al., 2007) and separated by reverse phase nano-LC-MS/MS (column: 12 cm/75um C18 built-in-emitter column, Nikkyo Technos Co., Ltd. Japan, EasyLC 1200, Thermo Scientific) using a 70-min analytical gradient, increasing from 2% B/98% A to 38%B/62% A (A: 0.1% formic acid, B: 80% Acetonitrile/0.1% formic acid) at 300 nL/min. The mass spectrometer (Fusion Lumos, Thermo Scientific) was operated in high/high mode (120,000 and 30,000 for MS1 and MS2, respectively). Auto Gain Control was set at 50,000 for MS2. MS1 scan range was set to m/z 375–1500 and m/z 110 was set as lowest recorded mass in MS2. One-point lock mass calibration was used. All data were quantified and searched against a Uniprot mouse database using MaxQuant (v. v. 1.6.0.13) (Cox et al., 2014). Oxidation of methionine and protein N-terminal acetylation were allowed as variable modifications, cysteine carbamidomethyl was set as a fixed modification, and two missed cleavages were allowed. The ‘match between runs’ option was enabled, and false discovery rates for proteins and peptides were set to 1%. Protein abundances measured using label free quantitation (Tyanova et al., 2016).
Proteomic data analysis
Proteomic datasets were analyzed using Perseus v1.6.14.0 (Tyanova et al., 2016). Of the detected proteins, those flagged as reverse, only identified by site, and potential contaminants were excluded from the analysis. For quantitative analysis, LFQ or iBAQ intensities were employed as indicated; LFQ intensities were used for comparisons across samples, while iBAQ intensities were used to compare abundances across different proteins. Imputation of undetected data points for Log2(LFQ) intensities was performed by assigning values from a normal distribution of 0.3 width and 1.8 down shift. Principal component analysis (PCA) was performed with imputed Log2(LFQ) intensities. Scatterplot representation of Log2(LFQ) intensities was generated without imputation. Differentially secreted proteins were identified with ANOVA using permutation-based FDR, with FDR set at 0.01 and number of randomizations at 250.
GO cell component analysis
Gene symbols from detected proteins were submitted to Retrieve/ID mapping tool on the UniProt website (https://uniprot.org). List of genes that are annotated with the following gene ontology cell component terms were obtained: extracellular region (5576), extracellular space (5615), extracellular matrix (31012), plasma membrane (5886), cytosol (5829), nucleus (5634), mitochondrion (5739), endoplasmic reticulum (5783), and Golgi apparatus (5794).
Secretion prediction analysis
UniProt accession IDs of the detected proteins were submitted to Retrieve/ID mapping tool on the UniProt website (https://uniprot.org) to obtain the FASTA sequences, which were used as inputs for various secretion prediction algorithms using the web-based query system. We defined classically secreted proteins as having SignalP5.0 score >0.5 and 0 or 1 predicted transmembrane domains by TMHMM2.0. Subcellular localization prediction analysis was performed using DeepLoc1.0 and searched for proteins whose predicted location is extracellular. PredGPI specificity score >99% was used as the threshold to determine if a protein is expected to be GPI-anchored. P roteins with SecretomeP2.0 score >0.6 and SignalP5.0 score ≤0.5 were considered non-classically secreted.
Cluster analysis and functional annotation
Hierarchical clustering was performed on z-score-transformed Log2(LFQ) values using the complete-linkage method and split into 4 clusters by dendrogram. DAVID v6.8 was used to generate functional annotation of clusters (Huang et al., 2009). List of genes encoding the proteins of each of the 4 clusters were compared to a background gene list of total detected proteins in the proteomic dataset. Per cluster, top 4 overrepresented pathways in the gene ontology biological process terms were reported.
Adipose tissue enrichment analysis
To identify candidate genes enriched within brown and white adipose relative other tissue types as well to each other we used the BioGPS datasets (Su et al., 2004). The Mouse GNF1M Gene Atlas datasets (GSE1133) were downloaded from BioGPS portal (Su et al., 2004) and imported into Limma Bioconductor package (Ritchie et al., 2015) for Log2 transformation and differential expression analysis. All pair-wise comparisons for both brown adipose and white adipose tissues against all other tissue types were performed using limma as well as the direct comparison between brown and white adipose tissues. Genes with a Log2 fold change greater than 4 and a Benjamini-Hochberg-corrected FDR of 0.05 within pair-wise comparisons were considered significantly enriched. Genes were further scored by the total number of pair-wise comparisons where genes were found to be enriched in both adipose tissues, brown adipose tissue or white adipose tissue compared to other tissues in the tissue atlas.
RNA isolation, cDNA synthesis, and RT-qPCR
Total RNA was extracted from cultured cells using RLT buffer (Qiagen) and from tissues using TRIzol (Invitrogen) and purified using RNeasy Mini Kit (Qiagen). cDNA was synthesized from 1 µg of RNA using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosciences). Power SYBR Green (Life Technologies) was used for RT-qPCR reactions performed with QuantStudio 6 Flex Real-Time PCR System (Thermo Scientific) in a 384-well format. Relative fold changes of mRNA levels were calculated using the ΔΔCT method with 18 S rRNA as loading control. qPCR primers are provided in Supplementary file 3.
Adipose tissue fractionation
Adipose tissues from 8-week-old C57BL/6 J WT male mice were dissected and minced. eWAT and iWAT were digested in a buffer containing 10 mg/mL collagenase D (Roche), 2.4 mg/mL Dispase II (Roche), and 10 mM CaCl2 in PBS. For BAT, 2 x BAT digestion buffer containing 125 mM NaCl, 5 mM KCl, 1.3 mM CaCl2, 5 mM glucose, 1% P/S, and 4% BSA was prepared, which was diluted 1:1 with PBS and used to dissolve collagenase B (Roche) at a final concentration of 1.5 mg/mL. Following collagenase digestion of the tissues in a 37 °C water bath, the mature adipocyte fraction was separated from the SVF pellet by centrifugation at 500×g for 10 min at 4 °C. The two fractions were transferred to two separate tubes, washed with DMEM/F-12 GlutaMAX containing 10% FBS and 1% P/S, and vortexed in TRIzol for RNA extraction.
RNA-Sequencing and immune cell deconvolution
Extracted RNA samples were analyzed for RNA integrity number (RIN) using the Bioanalyzer (Agilent) and sequenced using Illumina NovaSeq at the Rockefeller University Genomics Resource Center. Reads were trimmed with Cutadapt, aligned to mm10 reference genome using STAR, and quantified using featureCounts. Differential gene expression analysis was performed using DESeq2 (Love et al., 2014). Pathway analysis was performed using clusterProfiler (Yu et al., 2012). Deconvolution analysis was performed with CIBERSORTx (Newman et al., 2019) using ImmuCC signature matrix (Chen et al., 2017).
Tissue expression pattern analysis
Mouse tissue mRNA sequencing data from ENCODE was downloaded from GSE36026. Reads were mapped and quantified as above and gene expression was normalized using DESeq2. 177 genes encoding 180 proteins detected in mouse nascent serum were selected for t-SNE analysis. Briefly, average normalized expression in a tissue was divided by summed expression across tissues. Tissue with the highest relative expression was designated as the highest expressing tissue for a gene. t-SNE analysis was performed on relative expression values with R package Rtsne version 0.16 (Krijthe, 2015) using a perplexity of 30 and maximum iteration of 1000.
Interstitial fluid collection
Mice were sacrificed, then carefully perfused with approximately 20 mL dPBS. Whole tissues were carefully excised with minimal tissue disturbance and laid flat onto a 70 µm nylon mesh over a tube containing 2 µL of 10 x protease inhibitor solution (Roche). Tubes were sealed and centrifugated at 300 g at 4 °C. Following collection, aqueous infranatant was collected for downstream processing.
Immunoblot
Upon collection, conditioned medium (CM) was filtered using a 0.22 µm PES membrane syringe filter unit to remove cell debris and concentrated using a 3 kDa centrifugal filter units (Millipore). Tissue protein lysates were collected by homogenizing flash-frozen tissue in RIPA buffer (10 mM Tris pH 8.0, 1 mM EDTA, 1% Triton X-100, 0.1% Na deoxycholate, 0.1% SDS, and 140 mM NaCl) containing cOmplete mini protease inhibitor cocktail (Roche) and PhosSTOP (Roche). Protein concentration was determined using Pierce BCA Protein Assay Kit (Thermo Scientific) using a dilution series of bovine serum albumin as protein standards. Mouse serum samples were loaded at equal volume. Pre-cast polyacrylamide gels were used for electrophoresis, after which protein was transferred to PVDF membrane using standard techniques. Immunoblots were incubated with indicated primary antibodies and developed using Western Lightning Plus-ECL (PerkinElmer) and imaged on an autoradiographic film or using a Bio-Rad Gel Doc system.
Generation of LRG1-depleted FBS
Affi-Gel10 (Bio-Rad) was washed with cold distilled water and resuspended in binding buffer (0.2 M NaHCO3 and 0.5 M NaCl). Equine Cyt c was dissolved in binding buffer at 40 mg/mL concentration and conjugated with Affi-Gel10 for 1 hr at room temperature. Reaction was quenched with buffer containing 0.5 M ethanolamine and 0.5 M NaCl pH 8.3 for 1 hr. Conjugated beads were washed with copious amounts of PBS and elution buffer (0.1 M acetate and 0.5 M NaCl pH 4.0). For depletion of LRG1 from FBS, 2 mL of 50% Cyt c-agarose bead suspension was added to 50 mL of FBS and incubated overnight at 4 °C and repeated once more with fresh conjugated beads. FBS was sterile-filtered before use. Specific binding of LRG1 to Cyt c-agarose beads was confirmed by incubating the beads used for depletion with elution buffer for 10 min at room temperature and subjecting the supernatant to SDS-PAGE electrophoresis followed by Coomassie stain or western blot against LRG1.
Co-immunoprecipitation
Serum was diluted 1:1 with PBS containing 0.02% Tween-20 (PBS-T) and incubated with α-FLAG M2 beads (Sigma) overnight at 4 °C. Following wash with PBS-T, bound proteins were eluted by heating the beads at 70 °C in Laemmli buffer containing 50 mM glycine buffer pH 2.8 and 9% (v/v) β-mercaptoethanol.
Adenovirus production and purification
Adenoviral vectors were created using the AdEasy system (Luo et al., 2007). C-terminally 3xFLAG-tagged murine LRG1 was cloned into pAdTrack-CMV (AddGene) linearized with XhoI (NEB) and HindIII (NEB) using In-Fusion HD Cloning Kit (Takara). pAdTrack-CMV and pAdTrack-CMV-LRG1-FL plasmids were linearized with PmeI (NEB) and recombined into pAdEasy-1 vector via electrophoretic transformation of recombination-competent BJ5183-AD-1 cells (Agilent) with the linearized product and selection for kanamycin-resistant clones. Plasmids from validated clones were transformed into recombination-deficient XL-10 Gold ultracompetent cells (Agilent), which were used to generate pAd-eGFP and pAd-LRG1-FL plasmids and purified using Plasmid Maxi Kit (Qiagen).
Crude adenovirus was produced by transfecting PacI (NEB)-linearized pAd vectors into HEK293A cells (Invitrogen), which were incubated at 37 °C with 5% CO2 for 10–14 days with media supplementation every 3–5 days until most cells showed cytopathic effect/detachment. Both cells and the culture medium were collected, lysed by 3 cycles of freeze-thaw between dry ice-ethanol and room-temperature water baths, and centrifuged at 3500×g for 15 min at 4 °C to obtain the supernatant crude virus. Round 1 amplification product was obtained by transducing HEK293A cells with the crude virus and repeating the above collection, lysis, and centrifugation steps.
To obtain round 2 amplification product, twelve 15 cm plates of HEK293A cells were transduced with round 1 adenovirus and incubated at 37 °C with 5% CO2 until most cells demonstrated cytopathic effect. As with previous rounds, cells and media were collected, lysed by freeze-thaw cycles, and centrifuged to obtain the supernatant. The supernatant was treated with benzonase, and adenoviral particles were purified from the crude mixture using the Vivapure AdenoPACK 100 kit (Sartorius). Purified virus was dialyzed with buffer containing 20 mM Tris pH 8, 25 mM NaCl, and 2.5% (w/v) glycerol and concentrated using a 100 kDa centrifugal filter unit provided with the kit. Titer of the adenovirus was determined using Adeno-X Rapid Titer Kit (Takara).
AAV8 vector preparation
C-terminally 3xFLAG-tagged LRG1 was cloned into pENN.AAV.CB7.CI.eGFP.WPRE.rBG (Addgene) linearized with EcoRI (NEB) and BglII (NEB) using In-Fusion HD Cloning Kit (Takara). The original eGFP-expressing and cloned LRG1-FL plasmids were transformed into Stable Competent E. coli (NEB), purified using Plasmid Maxi Kit (Qiagen), and shipped to Penn Vector Core (PA, USA) for AAV8 production.
In vivo adenovirus/AAV8 transduction
In vivo adenoviral transduction studies were performed using purified adenovirus from second round of amplification. Adenovirus was injected at a dose of 1010 pfu/mouse. AAV8 was injected at 1011 GC/mouse. The mice were briefly anesthetized with isoflurane for virus injection via the retroorbital route. Following injection, the mice were quarantined in an ABSL-2 housing room for 72 hr before transferred back to regular housing conditions.
Histology
Dissected tissues were fixed in 10% neutral buffered formalin for 3 days at room temperature and transferred to 70% ethanol. Paraffin embedding, sectioning, and H&E staining was performed by the Memorial Sloan Kettering Cancer Center Laboratory of Comparative Pathology. H&E-stained sections were imaged using a wide-field fluorescence/brightfield/DIC microscope (Zeiss) at the Rockefeller University Bioimaging Resource Center. Crown-like structures (CLS) were identified as any adipocyte in a field of view with cellular infiltrates indicated by nuclear staining surrounding a majority of adipocyte perimeter. Objectives were used as indicated and chosen based on CLS enumerability (<30 CLS per field). 5 fields from 3 animals were quantified per group. Adipocyte area was measured from a minimum of 289 adipocytes from 2 to 4 animals per group using Cellpose 2.0 for image segmentation (Stringer et al., 2021).
Liver Triglyceride Quantification
Per 100 mg of flash-frozen liver, 1 mL of 5% NP-40 lysis buffer was added to generate a homogenate. The homogenized tissue samples were then subject to two cycles of heating at 80–100°C for 5 minutes followed by cooling at room temperature. The resulting mixture was centrifugated, and supernatant was collected. Triglyceride concentration of the supernatant was measured using Triglyceride Quantification Colorimetric/Fluorometric Kit (Sigma, MAK266) following manufacturer instructions.
Multiplex cytokine panel
Mouse serum samples were diluted 1:1 with PBS, snap frozen using liquid N2 and shipped to Eve Technologies (Alberta, Canada) on dry ice. The Mouse Cytokine Array / Chemokine Array 31-Plex (MD31) panel was used to quantify the levels of cytokines and chemokines.
Cytochrome c clearance assay
12-week-old chow-fed LRG1-KO and WT littermate males were injected retro-orbitally with 40 mg/kg equine cytochrome c (Cyt c) in PBS. Blood was collected from the tail vein immediately prior to injection (0 min) and 5 min, 15 min, 30 min, 1 h, 2 h, 4 h, and 6 h post-injection into EDTA-coated capillary tubes and kept in ice until further processing. During blood collection, each mouse was placed on a metal grating above a clean plastic wrap to allow collection of excreted urine, if any. Plasma was isolated via centrifugation at 2000×g for 15 min at 4 °C. Immunoblot against Cyt c was performed with WT and KO plasma samples run pairwise to enable relative quantification of Cyt c signal. Quantification was performed using ImageJ (Schneider et al., 2012).
Statistical analysis
Unless otherwise noted, data are presented as mean ± SEM, with n number specified in the figure legends. Statistical analyses were performed with GraphPad Prism 9. Binary comparisons were performed with Welch’s t-test to account for possible difference in variance. Statistical analysis of data involving 3 or more conditions (levels) of a single variable was performed using one-way ANOVA followed by post hoc tests using the Dunnett method to compare every mean with a control mean or the Tukey method to compare every mean with every other mean. Data measured across multiple time points as in GTT, ITT, and PTT were analyzed with repeated measures two-way ANOVA, reporting group factor p-values. Analysis of data from a two-factor experimental setup was performed with two-way ANOVA or two-way mixed effects ANOVA in the case of an uneven n number, reporting group factor p-values. For post hoc tests, the Tukey method was used when comparing every mean with every other mean and the Šídák method was employed when a selected set of means were compared.
Materials availability
The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD035318. RNA-Seq data have been deposited to GEO with the identifier GSE208219. Reagents including unique biological materials are available from the corresponding author upon request.
Code availability
Publicly available codes are available in the relevant references.
Appendix 1
Data availability
The proteomic dataset has been deposited to Proteomexchange PRIDE under accession PXD035318. RNA-Seq data have been deposited to GEO under accession GSE208219. All original gels and blots are available as Source Data Files.
-
PRIDEID PXD035318. LRG1 is an adipokine that promotes insulin sensitivity and suppresses inflammation.
-
NCBI Gene Expression OmnibusID GSE208219. AAV-mediated LRG1 overexpression in diabetic mice.
-
NCBI Gene Expression OmnibusID GSE36026. RNA-seq from ENCODE/LICR.
-
NCBI Gene Expression OmnibusID GSE1133. tissue-specific pattern of mRNA expression.
References
-
Serum cytochrome c level as a prognostic indicator in patients with systemic inflammatory response syndromeClinica Chimica Acta; International Journal of Clinical Chemistry 342:127–136.https://doi.org/10.1016/j.cccn.2003.12.011
-
Comparative secretome analyses of primary murine white and brown adipocytes reveal novel adipokinesMolecular & Cellular Proteomics 17:2358–2370.https://doi.org/10.1074/mcp.RA118.000704
-
Cytochrome-C detection: a diagnostic marker for myocardial infarctionApplied Biochemistry and Biotechnology 90:97–105.https://doi.org/10.1385/abab:90:2:97
-
SignalP 5.0 improves signal peptide predictions using deep neural networksNature Biotechnology 37:420–423.https://doi.org/10.1038/s41587-019-0036-z
-
Cell-type-specific metabolic labeling of nascent proteomes in vivoNature Biotechnology 35:1196–1201.https://doi.org/10.1038/nbt.4016
-
Serum cytochrome c indicates in vivo apoptosis and can serve as a prognostic marker during cancer therapyInternational Journal of Cancer 116:167–173.https://doi.org/10.1002/ijc.21037
-
Proteomics biomarkers for solid tumors: current status and future prospectsMass Spectrometry Reviews 38:49–78.https://doi.org/10.1002/mas.21572
-
Circulating soluble cytochrome c in liver disease as a marker of apoptosisJournal of Internal Medicine 254:168–175.https://doi.org/10.1046/j.1365-2796.2003.01171.x
-
Feature-based prediction of non-classical and leaderless protein secretionProtein Engineering, Design & Selection 17:349–356.https://doi.org/10.1093/protein/gzh037
-
Efficient proximity labeling in living cells and organisms with turboidNature Biotechnology 36:880–887.https://doi.org/10.1038/nbt.4201
-
Adipocyte death defines macrophage localization and function in adipose tissue of obese mice and humansJournal of Lipid Research 46:2347–2355.https://doi.org/10.1194/jlr.M500294-JLR200
-
CRISPOR: intuitive guide selection for CRISPR/cas9 genome editing experiments and screensNucleic Acids Research 46:W242–W245.https://doi.org/10.1093/nar/gky354
-
Accurate proteome-wide label-free quantification by delayed normalization and maximal peptide ratio extraction, termed maxlfqMolecular & Cellular Proteomics 13:2513–2526.https://doi.org/10.1074/mcp.M113.031591
-
The encyclopedia of DNA elements (encode): data portal updateNucleic Acids Research 46:D794–D801.https://doi.org/10.1093/nar/gkx1081
-
Selective enrichment of newly synthesized proteins for quantitative secretome analysisNature Biotechnology 30:984–990.https://doi.org/10.1038/nbt.2356
-
Combining pulsed SILAC labeling and click-chemistry for quantitative secretome analysisMethods in Molecular Biology 1174:101–114.https://doi.org/10.1007/978-1-4939-0944-5_7
-
Beyond adiponectin and leptin: adipose tissue-derived mediators of inter-organ communicationJournal of Lipid Research 60:1648–1697.https://doi.org/10.1194/jlr.R094060
-
LRG1 is an adipokine that mediates obesity-induced hepatosteatosis and insulin resistanceThe Journal of Clinical Investigation 131:e148545.https://doi.org/10.1172/JCI148545
-
LRG1 promotes diabetic kidney disease progression by enhancing TGF-β-induced angiogenesisJournal of the American Society of Nephrology 30:546–562.https://doi.org/10.1681/ASN.2018060599
-
Metabolic functions of fabps—mechanisms and therapeutic implicationsNature Reviews Endocrinology 11:592–605.https://doi.org/10.1038/nrendo.2015.122
-
Release of intact, monomeric cytochrome c from apoptotic and necrotic cellsCell Death and Differentiation 9:538–548.https://doi.org/10.1038/sj.cdd.4400981
-
Paradoxical roles of leucine-rich α2-glycoprotein-1 in cell death and survival modulated by transforming growth factor-beta 1 and cytochrome cFrontiers in Cell and Developmental Biology 9:744908.https://doi.org/10.3389/fcell.2021.744908
-
Animal models of obesity and diabetes mellitusNature Reviews. Endocrinology 14:140–162.https://doi.org/10.1038/nrendo.2017.161
-
Predicting transmembrane protein topology with a hidden markov model: application to complete genomesJournal of Molecular Biology 305:567–580.https://doi.org/10.1006/jmbi.2000.4315
-
Regulation of metabolism by the innate immune systemNature Reviews. Endocrinology 12:15–28.https://doi.org/10.1038/nrendo.2015.189
-
Identification and validation of novel adipokines released from primary human adipocytesMolecular & Cellular Proteomics 11:M111.https://doi.org/10.1074/mcp.M111.010504
-
Pulsed azidohomoalanine labeling in mammals (PALM) detects changes in liver-specific LKB1 knockout miceJournal of Proteome Research 14:4815–4822.https://doi.org/10.1021/acs.jproteome.5b00653
-
Association of circulating proinflammatory marker, leucine-rich-α2-glycoprotein (LRG1), following metabolic/bariatric surgeryDiabetes/Metabolism Research and Reviews 34:e3029.https://doi.org/10.1002/dmrr.3029
-
Serum starvation: caveat emptorAmerican Journal of Physiology. Cell Physiology 301:C272–C279.https://doi.org/10.1152/ajpcell.00091.2011
-
Obesity and cardiovascular disease: pathophysiology, evaluation, and effect of weight loss: an update of the 1997 american heart association scientific statement on obesity and heart disease from the obesity committee of the council on nutrition, physical activity, and metabolismCirculation 113:898–918.https://doi.org/10.1161/CIRCULATIONAHA.106.171016
-
Circulating levels of cytochrome c after resuscitation from cardiac arrest: a marker of mitochondrial injury and predictor of survivalAmerican Journal of Physiology-Heart and Circulatory Physiology 292:H767–H775.https://doi.org/10.1152/ajpheart.00468.2006
-
Identification of novel human adipocyte secreted proteins by using SGBS cellsJournal of Proteome Research 9:5389–5401.https://doi.org/10.1021/pr100621g
-
Inflammatory mechanisms linking obesity and metabolic diseaseThe Journal of Clinical Investigation 127:1–4.https://doi.org/10.1172/JCI92035
-
Nih image to imagej: 25 years of image analysisNature Methods 9:671–675.https://doi.org/10.1038/nmeth.2089
-
Molecular architecture of the chick vestibular hair bundleNature Neuroscience 16:365–374.https://doi.org/10.1038/nn.3312
-
Generation of mouse model (KI and CKO) via easi-CRISPRMethods in Molecular Biology 2224:1–27.https://doi.org/10.1007/978-1-0716-1008-4_1
-
Cellpose: a generalist algorithm for cellular segmentationNature Methods 18:100–106.https://doi.org/10.1038/s41592-020-01018-x
-
Pathology of the liver in obese and diabetic ob/ob and db/db mice fed a standard or high-calorie dietInternational Journal of Experimental Pathology 92:413–421.https://doi.org/10.1111/j.1365-2613.2011.00793.x
-
Brown adipose tissue as a secretory organNature Reviews. Endocrinology 13:26–35.https://doi.org/10.1038/nrendo.2016.136
-
Adiponectin, the past two decadesJournal of Molecular Cell Biology 8:93–100.https://doi.org/10.1093/jmcb/mjw011
-
Cell type-selective secretome profiling in vivoNature Chemical Biology 17:326–334.https://doi.org/10.1038/s41589-020-00698-y
-
Obesity is associated with macrophage accumulation in adipose tissueThe Journal of Clinical Investigation 112:1796–1808.https://doi.org/10.1172/JCI19246
-
Cytochrome c can be released into extracellular space and modulate functions of human astrocytes in a toll-like receptor 4-dependent mannerBiochimica et Biophysica Acta. General Subjects 1863:129400.https://doi.org/10.1016/j.bbagen.2019.07.009
-
BioGPS: building your own mash-up of gene annotations and expression profilesNucleic Acids Research 44:D313–D316.https://doi.org/10.1093/nar/gkv1104
-
Leucine-rich α-2-glycoprotein predicts proliferative diabetic retinopathy in type 2 diabetesJournal of Diabetes and Its Complications 33:651–656.https://doi.org/10.1016/j.jdiacomp.2019.05.021
-
Quantitative analysis of secretome from adipocytes regulated by insulinActa Biochimica et Biophysica Sinica 41:910–921.https://doi.org/10.1093/abbs/gmp085
Article and author information
Author details
Funding
American Diabetes Association (1-17-ACE-17)
- Paul Cohen
National Institutes of Health (RC2 DK129961)
- Paul Cohen
National Institute of General Medical Sciences (T32GM007739)
- Chan Hee J Choi
- Sarah K Szwed
Sarnoff Cardiovascular Research Foundation
- Samir Zaman
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Jeffrey M Friedman, James C Lo, Howard C Hang, Jingyi Chi, and Kaja Plucińska for feedback and discussions. We also thank Chingwen Yang of the CRISPR and Genome Editing Center at the Rockefeller University for generating the LRG1-KO line. We are grateful to Tao Tong and Alison North of the Bio-Imaging Resource Center at the Rockefeller University for assistance with microscopy. CHJC and SKS were supported by a Medical Scientist Training Program grant from the National Institute of General Medical Sciences of the National Institutes of Health under award number T32GM007739 to the Weill Cornell/Rockefeller/Sloan Kettering Tri-Institutional MD-PhD Program. SZ was supported by the Sarnoff Cardiovascular Research Fellowship Program. PC was supported by the American Diabetes Association Pathway to Stop Diabetes Grant 1–17-ACE-17 and by NIH grant RC2 DK129961. Data was generated by the Proteomics Resource Center at The Rockefeller University (RRID:SCR_017797) using instrumentation funded by the Sohn Conferences Foundation and the Leona M and Harry B Helmsley Charitable Trust. Bio-Render was used to make the schematic diagrams. The content of this study is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Ethics
All animal studies were performed in accordance with the institutional guidelines of the Rockefeller University Institutional Animal Care and Use Committee (IACUC) protocol (21010-H). Experiments involving adenoviral and AAV8 vectors were performed under general anesthesia using isoflurane, in accordance with the institutional ABSL-2 guidelines.
Copyright
© 2022, Choi 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
-
- 2,406
- views
-
- 496
- downloads
-
- 12
- 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
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.
-
- 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.