13 C 6 Lysine isotope labeling study.
Male Bl6/N mice (12 weeks old) were fed a custom diet (Silantes) containing more than 99% 13C6 lysine. The protocol for labeling was described previously30. After 1, 2, and 8 weeks (protocol 1) or 1, 2, and 3 weeks (protocol 2), mice were sacrificed and perfused with ice-cold PBS. All denominated organs were snap-frozen from the same mouse and stored at -80 °C degrees. Blood was separated in erythrocyte and plasma by centrifugation.
Metabolomics sample preparation.
Frozen organs were kept on dry ice and 10 mg tissue was weighed. Then, 800 ul ice-cold extraction solution containing Acetonitril : Methanol : Water 2:2:1 was added. Samples were subjected to tissue homogenization using a multiplex-bead-beater (Storm) for 30 sec (liver), 1 min (lung), or 45 sec (all other tissues). The supernatant was transferred to a prechilled Eppendorf tube. Then, beads were washed with 200 ul of the extraction solution, and the wash solution was added to the remaining homogenate. The homogenate was incubated at -20 °C degrees for 2 hrs. Then, tissue was spun down at 4 °C degrees for 20 min, and the supernatant was transferred to another vial. The supernatant containing the extract was transferred into a speed vac and dried down. The next morning, the dried-down extract was once resuspended in 100 ul (per 10 mg tissue) of acetonitrile-water 1:1, and the solution was centrifuged at 4 °C degrees. Then, the solution was transferred to autosampler vials and stored at -80 °C until further use.
Untargeted metabolomics analysis and mass spectrometry.
LC-MS/MS analysis for metabolomics was performed as previously described22. Data was annotated with in-source fragments and adducts. Quality controls were run every 5 samples. The mass spectrometer was initially calibrated using NaFormate peaks and in addition post-run. For untargeted metabolomic analysis, we used a UHPLC-MS approach. For fractionation, we used hydrophilic interaction liquid chromatography (HILIC) fractionation and reversed-phase (RP) chromatography as previously described22. We used a quadrupole time-of-flight instrument (Impact II, Bruker, Bremen, Germany) coupled to a ultrahigh-performance liquid chromatography (UHPLC) device (Bruker Elute, Bruker, Billerica, MA), or to an Agilent Infinity 1290 UHPLC device (Agilent, USA). The MS was calibrated using sodium formate (post-run mass calibration). Data were acquired over an m/z range of 50 to 1000 Da in positive ion mode and negative ion mode (HILIC only). Electrospray source conditions were set as follows: end plate offset, 500 V; dry gas temperature, 200 °C; drying gas, 6 liters/min; nebulizer, 1.6 bar; and capillary voltage, 3500 V.
To increase metabolome coverage and minimize ion suppression, we used a dual fractionation strategy. For RP separation, an ACQUITY BEH C18 column (1.0 × 100 mm, 1.7-µm particle size; Waters Corporation, Milford, MA) was used, and for HILIC fractionation, a ACQUITY BEH amide (1.0 × 100 mm, 1.7-µm particle size; Waters Corporation, Milford, MA) column was used. Flow was 150 µl/min, and a binary buffer system consisting of buffer A (0.1% FA) and buffer B (0.1% FA in acetonitrile) was used. The gradient for RP was: 99% A for 1 min, 1% A over 9 min, 35% A over 13 min, 60% A over 3 min, and held at 60% A for an additional 1 min. The gradient for HILIC consisted of 1% A for 1 min, 35% A over 13 min, 60% A over 3 min, and held at 60% A for an additional 1 min. The injection volume was always 2 µl. For molecule identification purposes, putative molecules of interest were fragmented using three different collision energies (10, 20 eV) or ramp collision energies (20 to 50 eV).
Untargeted metabolomics data analysis.
Bruker Raw files (*.d) were transformed into mzml files using the compassxport_converter.py script (Bruker). Then, data files were uploaded to XCMS online, and differential peaks were extracted in positive and negative ion mode42. Feature Detection was performed with the centwave method with the following options: ppm = 10, minimum peak width = 5, maximum peak width = 20, mzdiff = 0.01, Signal/Noise = 6, Integration method = 1, prefilter peaks = 3, prefilter intensity = 100, noise filter = 100. For Retention time correction, we used obiwarp method with profStep = 1. Alignment was perfomed with bw 05, minfrac = 0.5, mzwid = 0.015. mminsamp = 1, and max = 100. Statistical test was performed using an unpaired parametric t-test (Welsh), with post-hoc analysis = TRUE,. Statistical filtering was performed for priorization of features of interest, including a fold change of at least 1.5, and a p-value lower than 0.05, and a corrected p value (q-value) lower than 0.05. The analysis included PCA and multivariate analyses. Metabolites were identified based on 1) unique mass, 2) MS2 spectra comparison to authentic standard 3) coelution with authentic standard, and 4) isotopic pattern. The following adducts were routinely considered : Na, H, for positive mode, and Cl, formate for negative mode. Intensity data were plotted using instantclue32, ggplot2, or circosplot.R package.
Kidney dissection.
Kidney cortices were isolated under manual control (protocol 1, Fig. 1). Kidneys were manually dissected using a stereomicroscope. Based on anatomical criterions, the kidneys were dissected into cortex, inner stripe outer medulla, outer stripe outer medulla and inner medulla (protocol 2, Suppl. Figure 3). The samples were snap-frozen and stored at -80 °C or dry ice.
Proteomics sample preparation.
Kidney samples of 13C6 labeled kidneys were subjected to proteomics analysis using a tryptic in solution digestion protocol followed by nLC-MS/MS analysis. In brief, kidney pieces were minced and homogenized using a glass homogenizer in 8M urea containing 10 mM Ammoniumbicarbonate as well as protease and phosphatase inhibitor cocktail (Thermo). The homogenate was spun down at 6 °C for 20 min, and the supernatant was kept for further analysis. A small aliquot was subjected to protein measurement using BCA assay (Thermo). The proteins were reduced and alkylated using 5 mM DTT (30 min) and 10 mM IAA (1 hr) in the dark. Then, urea concentration was diluted to < 2M and LysC (1:100 w/w) ratio was added, and the mixture was digested for 16 hrs at 37 °C. The next day, the reaction was terminated by addition of 2% formic acid. Peptides were desalted using in-house made stage-tips and analyzed by nLC-MS/MS.
Proteomics analysis.
The peptides were separated by reverse-phase nanoflow-LC-MS/MS analysis and sprayed into a quadrupole-orbitrap tandem mass spectrometer (qExactive plus, thermo scientific) as previously described43. Raw Proteomics data were parsed with MaxQuant v 1.5.3.3.44, using a Uniprot RefSeq reference proteome database from Jan 2017, with using LysC (cuts after each Lysine) as a protease. Multiplexicity of the analysis was 2, with 13C6 labels in proteins as a modification in the second channel. The analysis has been previously described45. The non-normalized ratios (Heavy/13C6 over Light12C6) were used for further analysis using Perseus46 software suite and filtering for ratios in all experiments, as well as for annotation with GO terms.
Dahl salt sensitive (D/SS) rats.
The strain of re-derived Rapp Dahl SS rats used in studies (SS/JrHsdMcwi, RRID:RGD_1579902) has been inbred for more than 50 generations at Medical College of Wisconsin. Male and female animals at the age of 8 weeks were used for experiments. Rats were maintained on AIN-76A custom diet, either low salt (LS; 0.4% NaCl, # 113755, Dyets Inc.) or high salt (HS; 4% or 8% NaCl, # 113756 or # 100078, respectively, Dyets Inc.). Water and food were provided ad libitum.
Dahl SS (D/SS) rat is a widely used model of salt-induced hypertension and CKD. Since the derivation of the D/SS rat in 1962, there have been numerous phenotyping studies demonstrating the importance of the kidney in the regulation of blood pressure. Cowley et al. showed that upon consuming a high salt diet, the D/SS rat rapidly becomes hypertensive and exhibits severe renal damage, yet GFR is not changed until 10 days after consuming a high salt diet47.
Lysine treatment study.
Experimental animals received either vehicle (water) or L-Lysine (17 mg/ml) via drinking water (n = 6 per group). Blood pressure was measured by telemetry. D/SS rats were anesthetized with 2–3% (vol/vol) isoflurane and a blood pressure transmitter (PA-C40; DSI) was surgically implanted subcutaneously, with the catheter tip secured in the abdominal aorta via the femoral artery. After a 3-day recovery period, blood pressure was measured with a DSI system (“telemetry”) in conscious, freely moving SS rats under HS diet protocol, similar to those described previously48–50.
For urine collection, rats were placed in metabolic cages (no. 40615, Laboratory Products) for a 24 hrs urine collection. These urine samples were used to determine electrolytes, microalbumin, and creatinine. Whole blood and urine electrolytes and creatinine were measured with a blood gas and electrolyte analyzer (ABL system 800 Flex, Radiometer, Copenhagen, Denmark)48. Kidney function was determined by measuring albuminuria using a fluorescent assay (Albumin Blue 580 dye, Molecular Probes, Eugene, OR) read by a fluorescent plate reader (FL600, Bio-Tek, Winooski, VT).
Ketogenic diet study.
Male D/SS rats (n = 7 per group) were fed either control high salt diet (4% NaCl, # 113756, Dyets Inc.) or high salt ketogenic diet (Keto Diet 4% NaCl, TD.190564, Teklad Custom Diet) for 28 days. Food composition is depicted in supplementary Fig. 8A. Water and food were provided ad libitum. Blood pressure acquisition, urine and whole blood analyzes were performed as descried above in lysine treatment study.
Pulsed 13C6 Lysine labeling in the SS rat.
Male D/SS rats on HS protocol were administrated with 13C6 L-Lysine-2HCl (#1860969, Thermo Scientific) at day 13 HS diet (8% NaCl; 24 hrs before sacrifice). Intraperitoneal injection of 13C6 L-Lysine-2HCl (340 mM; 200 µl of PBS solution) was performed 24 hrs before sacrifice.
Albumin uptake cell culture studies.
Confocal microscopy was used to detect uptake of fluorescent albumin (AlexaFluor-647 albumin, 40 µg/ml) in confluent monolayers of OK proximal tubule epithelial cells51 (visualized by fluorescent F-actin, 488 nm), grown on transwell filters under orbital shear stress as previously described and coincubated with lysine or glycine52. Cells were incubated for 1 hr in serum-free culture media containing fluorescent albumin ), after pretreatment of 5–50 mM L-Lysine (or other amino acids) for 2 h, or after overnight treatment with 1–10 mM L-Lysine. Treatment was performed from the apical side unless otherwise indicated. Cell-associated albumin was quantified by spectrofluorimetry.
Intravital dual photon microscopy.
All surgical, imaging, rat fluorescent albumin, and image analysis procedures were performed as described previously53,54. Imaging was conducted using an Olympus FV1000 microscope adapted for an intravital two-photon microscopy with high-sensitivity gallium arsenide nondescanned 12-bit detectors. Animals were anesthetized with pentobarbital sodium (50 mg/ml). A jugular venous line was used to introduce fluorescent rat albumin (Texas Red labeled) and high-molecular weight dextran (150 kDa FITC-labeled, TdB Consultancy, Uppsala, Sweden).
Immunohistochemistry.
Rat kidneys were fixed in 10% formalin and processed for paraffin embedding as previously described55. Kidney sections were cut at 4 µm, dried, and deparaffinized for subsequent labeling by streptavidin-biotin immunohistochemistry. After deparaffinization, slides were treated with a citrate buffer (pH 6) for total of 35 min. Slides were blocked with a perioxidase block (Dako, Coppenhagen, Denmark), avidin block (Vector Laboratories, Burlingame, CA), biotin block (Vector Laboratories), and serum-free protein block (Dako). Tissue sections were incubated for 90 min in antibody to Kidney Injury Molecule-1 (Rat KIM-1Ab, 1:300, #AF3689, R&D Systems, Inc) or megalin (lipoprotein-related protein 2 (LRP2) Ab, 1:2500, from Dr. Franziska Theilig University of Kiel, Kiel, Germany). Secondary detection was performed with goat anti-goat or anti-rabbit biotinylated IgG (Biocare, Tempe, AZ) followed by streptavidin-horseradish peroxidase (Biocare) and visualized with diaminobenzidine (Dako). All slides were counterstained with Mayer's hematoxylin (Dako), dehydrated, and mounted with permanent mounting medium (Sakura, Torrance, CA).
Bioinformatic image analysis via convolutional neural net analysis.
Rat kidneys were cleared of blood, formalin fixed, paraffin embedded, sectioned, and mounted on slides as previously described56. Slides were stained with Masson’s trichrome stain. The localization and scoring of glomeruli was performed by a novel and robust application of convolutional neural nets23,57. A cumulative distribution plot was generated (OriginPro 9.0) using glomerular injury scores based on a scale of 0–4 as previously described58, and the probability for a corresponding score interval was calculated (more than 3,500 glomeruli per group). Cortex protein cast analysis was performed using a color deconvolution filter and Analyze Particles mode in the Fiji image application (ImageJ 1.51u, NIH).
Bioinformatic algorithms for mass-difference based isotope selection.
Labeled isotopologues were detected using an in-house isotracker script (see Data Availability and github submission) that operated as follows: first, we search for groups of features that, according to their m/z difference (< 10 ppm), could correspond to an isotopic envelope composed of a light isotope and heavy isotopes. Only features within 2 seconds of retention time difference were allowed to be grouped into a single isotopic envelope. Next, we compared the isotopic envelope between labeled and unlabeled samples. The envelopes presenting at least one heavy isotope with statistically significant higher abundance in labeled samples compared to unlabeled samples were retained for further analysis.
Correlation-based isotope selection approach.
Labeled isotopologues were detected using an in-house script (see Data Availability and github submission) that operated as follows: The custom R script is designed to identify isotopes from stable labeled isotope LC-MS data based on an input list of compounds and sum formulas and is available online including documentation via https://github.com/hpbenton/targeted_isotopes. It uses mzR and MSnBase from the bioconductor repository to open and manipulate the data. Once the raw data is opened each file is independently searched for a list of possible compound hits. These compounds are searched by creating a small bin of a user chosen ppm range around the mass. The vector of data undergoes a smoothing using a Savitzky Golay filter. Any compound that is above a given threshold (default 1000 counts) and is also above the chosen signal to noise is selected. Then, since the formula is known, any and all isotopes are searched within the same ppm range and at that retention time range. If the isotope peak also satisfies the above criteria the two vectors are correlated to help confirm a true positive isotope. Most isotopes are correlated above 0.959, Supplemental Fig. 1. The script used a list of 400 lysine metabolites derived from KEGG and METLIN as input.
Synthesis of Nε-malonyl-Lysine.
For a detailed description of Nε-malonyl-Lysine and Nα-malonyl-lysine isomers, please see the supplementary material and methods. Molecules were synthesized as described in the supplementary materal and methods and characterized by NMR and mass spectrometry. Volatiles were removed under reduced pressure and the crude product was purified by mass directed preparative reversed phase HPLC to give the formic acid salt which was directly used as an analytical standard and fragmented at 10, 20 and 20–50 eV in ESI in positive ion mode.
Malonyl-CoA and lysine in vitro reaction.
1 uM lysine and a 10 uM of malonyl-CoA was incubated together in 10 ul of PBS, pH = 8 at 37 degrees for 1 h. Both the isotope-labeled (13C6) and non-isotope-labeled (12C6) form of lysine were used, in order to exclude unspecific molecule products. Mixtures were analyzed on a QQQ as well as a QTOF machine. Specific transitions were used in order to detect both heavy and light forms of N-e-malonyl-lysine as well as lysine, and malonyl-CoA.
Malonyl-CoA assay.
Malonyl-CoA levels from tissue lysates were determined using a commercial rat ELISA assay (MyBiosource.com) according to the manufacturer’s manual.
Targeted metabolomics.
Targeted metabolomic analysis was performed on a triple-quadrupole (QQQ) mass spectrometer (Agilent Triple Quadrupole 6490, San Diego, CA), and the LC part was coupled to a high-performance liquid chromatography (HPLC) system (1290 Infinity, Agilent Technologies) coupled to ion funnel. For glycolysis and TCA product metabolite, a ZIC-pHILIC (Sequant column; 2.1 × 150 mm) was used for separation. Cycle time was 100 ms. Collision energies and product ions (MS2 or quantifier and qualifier ion transitions) were optimized. Electrospray ionization source conditions were set as follows: gas temperature, 250 °C; gas flow, 12 liters/min; Nebulizer, 20 psi; sheath gas temperature, 350 °C; cap voltage, 2000 V; and nozzle voltage, 1000 V. The gradient consisted of buffer A and buffer B. Buffer A was 95:5 H2O:acetonitrile, 20 mM NH4OAc, 20 mM NH4OH (pH 9.4). Buffer B was acetonitrile. The gradient with A/B ratios were as follows: T0, 10:90; T1.5, 10:90; T20, 60:40; T25, off. Five microliters was injected. For analysis of lysine metabolites, identical column and chromatography conditions as in the “untargeted metabolomics” section were used. In all cases, a standard curve was recorded and integrated using the mass hunter platform (Agilent). The method for the TCA cycle including transitions was previously published 22. The transitions used for malonyl-Lysine were as follows: 233 -> 84; 233 -> 129.09, 233 -> 147.10 (for non-labeled malonyl-Lysine) and 239 -> 89, 239 -> 134, 239 -> 153 for 13C6 labeled malonyl-Lysine).
Immunoblot.
Protein samples in RIPA buffer were loaded onto Novex 4–12% Bis-Tris gels (Life Tech) and were transferred onto nitrocellulose membranes with the Novex semi-dry transfer apparatus (Life Tech). After blocking in 5% milk-TBST for 1 h at room temperature, blots were incubated overnight in 5% BSA-TBST (1:1000 acetylated-lysine CST, 1:1000 malonyl-lysine CST, 1:5000 beta-actin Genscript A00702) at 4 deg Celsius. After washing in TBST, blots were incubated in 1:5000 HRP-conjugated secondary antibodies (mouse anti-rabbit Jackson Immunoresearch 211-032-171, rabbit anti-mouse Jackson Immunoresearch 211-035-109) in 1% milk-TBST. Blots were incubated with ECL Western blotting substrate (Pierce Scientific 32106) and were processed by autoradiography. The used ladder was Biorad Precision Plus All-blue ladder.
Study approval.
All studies using D/SS rats were conducted at Medical College of Wisconsin and protocols were approved by the MCW Animal Care and Use Committees and were performed in accordance with the standards set forth by the NIH Guide for the Care and Use of Laboratory Animals (National Academies Press, 2011). All mice studies were conducted at the MPI for Heart Lung and Blood Research Bad Nauheim as previously described in accordance to the local authorities31,60,61.
Statistical analysis.
Two-tailed t-tests were used for comparison unless otherwise indicated. Analysis of metabolomics and proteomics data is described in the respective sections.