Mice
Wildtype C57Bl/6J mice were obtained from The Jackson Laboratory and bred at the National Cancer Institute-Frederick Cancer Research and Development Center. Irg1−/− mice were obtained from Dr. Michael Diamond (Washington University School of Medicine, St. Louis, MO) and backcrossed 10 generations onto C57Bl/6J background. LysMCre or AlbCre mice on C57BL/6J background at NCI Frederick breeding core were crossed with Irg1fl/fl mice (generously provided by Dr. Diamond) to generate LysMCre x Irg1fl/fl and AlbCre x Irg1fl/fl mice, respectively. The C57Bl/6J background of all mice was confirmed by DartMouse (DartMouse Lab, Lebanon, NH; under NIH grants 1S10OD021667-01 and 5P30CA023108). Male mice were used at 8–10 weeks of age in accordance with an approved NCI Frederick Animal Care and Use Committee protocol. The mouse model of non-alcoholic fatty liver disease was induced by feeding mice ad libitum with RD Western diet (Research Diets Inc., New Brunswick, NJ) supplemented with 10% sucrose in drinking water for 12 weeks as described 15, 16, 43, 44. Some mice were treated 3x/week i.p. with 50 mg/kg 4-Octyl itaconate (Cayman Chemical, Ann Arbor, MI) in 40% cyclodextrin/PBS as described 18, 39 for 12 weeks.
Human tissue
Frozen liver tissue from deidentified 16 patients with non-alcoholic steatohepatitis (NASH) and 10 non-NASH controls were provided through the Clinical Biospecimen Repository and Processing Core of the Pittsburgh Liver Research Center. The non-NASH controls had diagnoses of hemangioma, focal nodular hyperplasia, hepatic adenoma, hepatocellular carcinoma and benign liver mass. Samples were homogenized in lysis buffer or 80% methanol and processed as described below for qPCR and metabolomics, respectively.
Histologic analysis/quantitation
Liver tissue sections were fixed in 4% buffered formalin and embedded in paraffin or snap frozen in OCT compound. Fixed, paraffin embedded tissue sections were incubated with antibody to F4/80 (clone BM4, eBioscience), followed by horseradish peroxidase or alkaline phosphatase conjugated secondary antibodies. Frozen tissues were stained in Oil red O (Sigma Aldrich). At least 5 non-overlapping fields per sample were imaged using Aperio ImageScope software (Leica Biosystems, Deer Park, IL) and quantitation was performed using ImageJ software using thresholding of object sizes to eliminate debris. Average ORO droplet size was calculated by dividing the sum of stained area by the total number of droplets in the same field.
Isolation of F4/80 + macrophages
Liver tissues were homogenized in 200U/mg type 1 collagenase (Worthington Biochemical, Lakewood, NJ), 1mg/ml dispase 2 (Sigma) and 0.5 mg/ml DNase 1 (Sigma) using a GentleMacs tissue dissociator (Miltenyi Biotec). The single cell suspension was filtered through a 70 uM filter. The washed cells were incubated with biotinylated anti-F4/80 antibody (clone BM8; BioLegend, San Diego, CA) followed by magnetically coupled streptavidin microbeads (Miltenyi Biotec). F4/80+ cells were positively selected by magnetic separation (Miltenyi BIotec) to > 97% purity. The flow through fraction, containing the non-macrophages was saved for analysis.
Quantitative PCR
Total RNA was isolated using High Pure RNA Isolation kits (Roche Diagnostics). RNA was reverse transcribed using the High Capacity cDNA Archive kit (Applied Biosystems). Murine and human Irg1 and IL-6 expression was examined using validated gene expression assays (Applied Biosystems). Briefly, 10 ng cDNA was put in a final volume of 20 µl containing 10 µl Taqman Universal PCR mix and 1 µl primer/probe gene expression assay. All samples were run on an ABI 7300 real-time PCR system and analyzed using the ΔΔCT method 45. Gene expression was normalized to the level of the housekeeping gene HPRT or GAPDH.
Metabolomics
For detection of itaconate and TCA intermediates, cell pellets (~ 1x10^6 cells) were washed and resuspended in 80% methanol. Phase separation was achieved by centrifugation at 4°C. The methanol-water phase containing polar metabolites was dried using a vacuum concentrator. Targeted measurements on the resuspended metabolite samples were performed through electrospray ionization mass spectrometry (ESI-LC–MS/MS) analysis as described 46.
For broad metabolomics analysis, cells were washed with cold 150 mM ammonium acetate and resuspended in cold 80% methanol/0.1M formic acid. Phase separation and lyophilization of samples was performed as described above. Lyophilized samples were reconstituted in 50% acetonitrile and analyzed in an Agilent 1290 Infinity II ultra-high performance liquid chromatography (UHPLC) system coupled with an Agilent 6546 quadrupole time-of-flight (Q-TOF) mass spectrometer. Acquired data was analyzed using Agilent MassHunter Profinder 8.0, which performed data processing of > 500 target metabolites derived from an in-house Accurate Mass Retention Time (AMRT) metabolite standard library 47.
Lipidomics
Cell pellets (~ 1x10^6 cells) were resuspended in methanol, extracted with chloroform and centrifuged for 5 min at 4C. Dried lipid extracts were reconstituted with 50 µL methanol/chloroform (9:1, v/v). Analysis was carried out an Agilent 1290 Infinity II ultra-high performance liquid chromatography (UHPLC) system coupled with an Agilent 6546 quadrupole time-of-flight (Q-TOF) mass spectrometer. Final curation of annotated lipids was performed with Agilent MassHunter Profinder (Version 10.0.2) as described 48. Agilent MassHunter Mass Profiler (Professional Version 15.1) was used for statistical analysis.
Acyl-CoA Analysis
Acyl-CoAs were extracted by resuspension of the cell pellets in 10% trichloroacetic acid, sonication, centrifugation at 17,000 x g 10min 4C, and extraction of the supernatant by 80% MetOH (final concentration) containing 25mM ammonium acetate. Supernatant from further centrifugation at 17,000 x g 10 min 4C was dried using a vacuum concentrator and then resuspended in 50 µL of water. Liquid chromatography-high resolution mass spectrometry analysis was conducted as previously described 49 on a Q Exactive Plus operating in positive ion mode coupled to a Vanquish Duo LC. Analysts were blinded to sample identity.
Biochemical analyses on Serum and Tissue Samples
Blood was collected at the indicated times and placed into serum separator tubes which were spun 12,000 x g for 5 min and stored frozen until use. Serum insulin levels were determined using a commercial ELISA (Cayman Chemical). Free fatty acids were determined using an enzyme-based assay kit (Abcam). A triglyceride colorimetric assay kit (Cayman Chemical) was used to determine triglycerides in serum samples or homogenized tissues.
For glucose tolerance testing, mice were fasted 16 hr. Blood glucose levels were tested using a commercially available blood glucose meter and test strips (CVS Health). Baseline glucose levels were read and then mice were challenged i.p. with 2 g/kG sterile glucose in 100 ul saline. For insulin tolerance testing, 5 hr fasted mice were injected i.p. with a single dose of insulin (0.5U/kG). Blood glucose was read every 15 min for 2 hr.
Preparation of hepatocytes
Hepatocytes were isolated from freshly isolated mouse livers (Wildtype mice) as described 42. Unless indicated otherwise, cultured hepatocytes were treated overnight in complete media (DMEM supplemented with 5% FBS; 1x insulin-transferrin-selenium (Thermo Fisher Scientific) with 10 mM sodium itaconate in the presence or absence of a 1:10 dilution of sterile, chemically defined lipid mixture 1 (2 ug/ml arachidonic acid; 10 ug/ml each linoleic, myristic, oleic, palmitic and stearic acids; 0.22 mg/ml cholesterol; Sigma-Aldrich, St. Louis, MO).
Extracellular flux analysis
Isolated hepatocytes were seeded at 1.7×104 cells per 96 well in complete media and treated as described overnight. Plated cells were washed gently with PBS and incubated with Seahorse assay media supplemented with 2mM glutamine and 25mM Glucose for 1h at 37°C with no CO2. Extracellular flux analysis was performed at 37°C with no CO2 using the XF-96 analyzer (Seahorse Bioscience) per the manufacturer’s instructions. Port additions and times were used as indicated in the figures.
Flow cytometric analysis
Cells (1 x 106) were incubated in cell staining buffer (0.1% BSA, 0.1% sodium azide) containing 250 µg/ml 2.4G2 ascites for 15 min. Cells were stained with either fluorescently conjugated antibodies to F4/80 (clone BM8) or Ly6 (clone 1A8), or the fluorescent probe BODIPY 493/503 (Thermo Fisher Scientific) for 20 min. After washing, the labeled cells were analyzed on an LSR II flow cytometer using FACS-DIVA software (Becton Dickinson).
Bulk RNA Sequencing
Total RNA was isolated from treated hepatocytes using a PureLink RNA Mini kit with on-column DNase digestion (Invitrogen). All samples met quality control of RIN > 8.0 using a Bioanalyzer. 500 ng total RNA was used via the stranded mRNA ligation kit (Illumina) as the input to an mRNA capture, cDNA synthesis and standard Illumina library prep prior to sequencing on the NextSeq 2000 instrument (Illumina, Inc). HiSeq Real Time Analysis software (RTA 3.9) was used for processing raw data files. Illumina bcl2fastq2.20 was used to demultiplex and convert binary base calls and qualities to fastq format. Trimmed sequencing reads were mapped to mouse reference genome mm10 and annotated using GENCODE STAR aligner (version 2.7.0f) with two-pass alignment. RSEM (version 1.3.1) was used for gene and transcript quantification based on GENCODE annotation files.
For RNAseq data processing, raw counts were uploaded to the Partek Flow platform and all analysis was carried out using default parameters. Counts were normalized and filtered to remove low expressed features (< 1.0). Batch removal was carried out using a general linear model. Data was visualized using principal component analysis and differential gene expression analysis was carried out using analysis of variance (ANOVA) testing comparing different treatment conditions. Gene set enrichment analysis was done as described 50. Briefly, raw counts were normalized with the median ratio method from DESeq2. Batch effects were removed and visualized using PCA. Normalized and batch corrected counts were imported to the GSEA software version 4.2.3 (Broad Institute). Gene sets were obtained from MSigDB (Broad Institute) or generated de novo in Gene Matrix Transposed (GMT) file format. Analysis was run using default parameters, with the exception of “Min size,” which was reduced to 2.
Statistical analysis
Statistical differences between groups were analyzed using GraphPad Prism software. Significance is indicated by *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001 values.