Preparation of GSE and inosine
The GSE powder was purchased from Changsha Tianwei Biotechnology Co., Ltd (Changsha, China). It was quantified to contain 95.98% proanthocyanidins with a water content < 0.65% and crude ash < 0.27%. Inosine powder was purchased from Sigma-Aldrich (Louis, USA), purity ≥ 99.00% (HPLC). Inosine was dissolved in DMSO to 20 mM and stored at -80℃ until used in vitro when it was diluted with different concentrations of the corresponding culture medium. Inosine used for intramuscular injection was dissolved in saline (ready-to-use) at a concentration of 100 mg/mL.
Animal experiments and sampling
Chinese Holstein cows used in this study received humane care according to the Guide for the Care and Use of Agricultural Animals in Research and Teaching (FASS, 2010). We conducted the GSE supplementation and inosine injection trial on the same farm, an 8,500-cow commercial intensive dairy farm in Inner Mongolia, China. All cows were selected to meet the basic criteria: clinically healthy, 2 ≤ parity < 6, previous 305-d mature equivalent milk yield > 8,000 kg, and with 3.0 ≤ body condition score (BCS) < 4.0. All cows had free access to water and total mixed ration (TMR) diet (Table S1).
For the GSE supplementation trial, 30 multiparous cows with similar expected calving dates (28 ± 3 d prior to calving) were selected and randomly assigned to two groups consisting of a non-GSE control (CON) and a GSE-fed group (0 or 15 g GSE/d per cow). Following a 7-day adaptation period (all cows were fed the same TMR diet), treatment began 21 days before expected calving and ended 21 days after calving. The GSE was fed individually to each cow during the morning while confined in head stalls for a 30 to 45 min until the supplement was consumed. During the treatment period, 8 cows (5 cows in CON and 3 cows in GSE) with less than 21 days of pre-calving feeding were excluded due to the gap between the expected and actual calving date, 2 cows in CON were excluded due to trauma-related lameness and metritis. This left 10 cows in CON and GSE groups that were included in the final analyses. A schematic diagram of the experiment is depicted in Fig. 1A. During the last 3 consecutive days, dry matter intake (DMI) and milk yield for each cow were recorded and the average of these 3 days used as representative for the experimental period. On the last day of the experiment, milk, serum, liver and subcutaneous adipose biopsies were collected and stored until use. The storage methods were as described in our previous work [32, 33].
Ruminal fluid was collected with a gastric tube rumen fluid sampler (A1141K, Anscitech, Winnipeg, CA) on the last day of the experiment, 1 h before the morning feeding. The first 2 tubes of fluid were discarded to avoid contamination, and 150 mL of fluid was then collected from each cow [34]. Ruminal fluid samples were filtered immediately through 4 layers of gauze and divided into 5 parts for the analysis of pH, VFA, ammonia nitrogen (NH3-N), bacteria, and metabolites. The pH was measured immediately, and the remaining 4 parts of filtered rumen fluid were placed in liquid nitrogen immediately to preserve them for use.
Inosine injection trials were conducted in two stages: experiment Ⅰ (Fig. S5A) and experiment Ⅱ (Fig. 6A). For experiment Ⅰ, 15 postpartum cows without clinical diseases and serum BHBA < 1.2 mM (monitored for 3 consecutive days from day 5 to 7 after calving) were selected and randomly assigned to 3 groups: Control (saline), inosine (4 mg/kg), and inosine (8 mg/kg). From day 7 to 13 after calving, cows received an intramuscular injection of inosine (4, 8 mg/kg, dissolved in saline) before the morning feeding. Control cows were treated with an injection of vehicle (normal saline). Serum samples were collected on the next day before the morning feeding and stored for use. For experiment Ⅱ, 20 cows with serum BHBA < 1.2 mM (Healthy, n = 10) or ≥ 1.2 mM (Ketosis, n = 10) for 3 consecutive days were selected and assigned to 4 groups consisting of Healthy + Saline (n = 5), Healthy + Inosine (n = 5), Ketosis + Saline (n = 5) and Ketosis + Inosine (n = 5). From day 7 to 13 after calving, cows received an intramuscular injection of inosine (8 mg/kg, dissolved in saline) before the morning feeding. Control cows were treated with an injection of vehicle (saline). On the 14th day after calving, milk, serum, liver and subcutaneous adipose biopsies were collected and stored until use.
Milk and serum sample analyses
Milk protein, milk fat and SCC were analyzed by mid-infrared spectroscopy using a MilkoScan Minor FT3 (FOSS, Denmark). Serum concentrations of Glu, NEFA, and BHBA, as well as the activities of alanine aminotransferase (ALT) and aspartate aminotransferase (AST) were determined using a Hitachi 7170 autoanalyzer with commercially-available kits (Glu, GL3815; NFFA, FA115; BHB, RB1008; ALT, AL3801; AST, AS3804; Randox Laboratories, Crumlin, UK). Serum concentrations of haptoglobin (HP) and serum amyloid A (SAA) were determined using commercial ELISA kits (HP: GWB-A43096, GenWay Biotech, San Diego, CA; SAA: LS-F12552, LifeSpan BioSciences Inc., Seattle, WA) according to the manufacturer’s instructions. Serum concentrations of inosine, hypoxanthine, cholesterol, spermidine and 4-hydroxyphen were determined by HPLC [35].
Hematoxylin and Eosin (H&E) and Oil red O staining
Liver and subcutaneous adipose tissues were fixed in 4% formaldehyde neutral buffer solution (P1110, Solarbio, China) and embedded in paraffin, floated on a water bath, picked up onto glass slides, and placed in slide racks. Slides were dewaxed with xylene, rehydrated through descending concentrations of alcohol, and stained with H&E (G1120, Solarbio, China). For Oil red O staining, a portion of liver tissue was frozen in OCT compound (4583, Sakura, Japan), sectioned at -18°C and fixed with 75% alcohol at room temperature for 15 min. Then, the slides were stained with Oil red O (O0625, Sigma-Aldrich, USA) and counterstained with hematoxylin.
Quantitative reverse transcription PCR (qRT-PCR)
Total RNA was prepared from collected tissue using RNAiso Plus (TaKaRa Japan) and reverse-transcribed into cDNA using a reverse transcription kit (TaKaRa Japan) according to the manufacturer’s instructions. Quantitative reverse transcription (qRT)-PCR was performed on a 7500 Real-Time PCR System (Applied Biosystems) using the SYBR Green plus reagent kit (TaKaRa Japan) to analyze relative mRNA abundance of target genes. Relative transcription of each target gene was normalized against the geometric mean of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and β-actin (ACTB) and calculated with the 2–ΔΔCT method. The primer pairs used in this study are listed in Supplemental Table S2.
Western blot analysis
Protein samples were separated on 10 or 15% Tris-glycine gels with a known pre-stained protein ladder (Shanghai Epizyme Biomedical Technology Co. Ltd, Shanghai, China). The membranes were blocked in Tris-buffered saline solution with 0.1% Tween-20 (TBST) containing 3% bovine serum albumin (Sigma-Aldrich Co. St. Louis, MO, USA) for 4 h at room temperature. The blocked membranes were incubated overnight at 4°C with primary antibodies against SREBP-1c (1:500; Cat. NB100-2215; Novus Biologicals, USA), PPARα (1:500; Cat. ab126285; Abcam, UK), p-HSL (1:1,000; Cat. 4139; Cell Signaling Technology; USA), HSL (1:1,000; Cat. 4107; Cell Signaling Technology, USA), ATGL (1:1,000; Cat. ab99532; Abcam, UK); p-AKT (1:1,000; Cat. 4060; Cell Signaling Technology, USA), AKT (1:1,000; Cat. 9272; Cell Signaling Technology, USA), p-NF-κB (1:1,000; Cat. 3033; Cell Signaling Technology, USA), NF-κB (1:1,000; Cat. 4764; Cell Signaling Technology, USA), p-IκBα(1:1,000; Cat. 2859; Cell Signaling Technology, USA), IκBα(1:1,000; Cat. 9242; Cell Signaling Technology, USA) and β-actin (1:1,000; Cat. sc-47778; Santa Cruz Biotechnology, USA). The primary and secondary antibodies were diluted in TBST. Removal of excess primary antibody was carried out by washing the membranes in TBST three times for 5 min each. Membranes were incubated with horseradish peroxidase-conjugated anti-rabbit (Proteintech, Chicago, USA) or anti-mouse (Proteintech, Chicago, USA) immunoglobulin at room temperature for 45 min. Excess secondary antibody was removed by washing the membranes in TBST three times for 5 min each. Lastly, immunoreactive bands were visualized via a Tanon Imaging System (Tanon 4600, Shanghai, China) using an enhanced chemiluminescence solution (Millipore, Bedford, USA). All bands were analyzed using Image-Pro Plus 6.0 (Media Cybernetics, Rockville, MD).
Ruminal fermentation parameters
Ruminal fluid pH was measured immediately after filtration using a portable pH meter (REX-PHBJ-260F, Shanghai, China). Ruminal NH3-N was determined using a colorimetric method [36], and the VFAs were measured by gas chromatography-mass spectrometry (GC-MS) at Majorbio Bio-Pharm technology Co., Ltd (Beijing, China).
Ruminal bacterial DNA extraction, 16S rRNA gene amplification, sequencing and data analysis
The bacterial DNA from 20 ruminal fluid samples was extracted using the E.Z.N.A. ®Stool DNA Kit (D4015, Omega, Inc., USA) according to manufacturer ’s instructions. The extraction quality of DNA was assessed by 1% agarose gel electrophoresis, and the concentration and purity of DNA were determined by NanoDrop2000. The V4 region of 16S rRNA gene was amplified by PCR using 515F (5'-GTGYCAGCMGCCGCGGTAA-3') and 806R (5'-GGACTACNVGGGTWTCTAAT-3'). After mixing the PCR products of the same sample, they were recovered by using 2% agarose gel, purified by Axyprep DNA gel extraction kit (Axygen Biosciences, Union City, California, USA), detected by 2% agarose gel electrophoresis, and detected and quantified by Quantum Fluorometer (Promega, USA). The NEXTFLEX Rapid DNA-Seq Kit was used to build the database. High-throughput sequencing was conducted on the Illumina MiseqPE250 platform. After filtering, clean reads were obtained from the raw reads and merged into tags using FLASH. Using DADA2, bacterial tags were clustered into amplicon sequence variants (ASVs) with 100% sequence similarity. Based on the database of Silva 16S rRNA (v138), the bacterial ASVs were classified by the RDP classifier (threshold was set to 0.7).
Alpha diversity was calculated using Mothur, and rarefaction curves were generated by the R software. The Bray-Curtis distance algorithm was used to quantify the distances between samples, visualization was via Principal Coordinate Analysis (PCoA) and evaluation for statistical significance was done via ANOSIM testing. For differential abundance analysis, the Wilcoxon signed rank test was used to analyze paired samples and adjust for multiple comparisons by FDR.
Ruminal untargeted metabolomics and analysis by liquid chromategraphy-mass spectrometry (LC/MS)
After thawing at room temperature, 20 samples were extracted with acetonitrile: methanol = 1:1(v: v), and an internal standard (0.02 mg/mL, L-2-chlorophenylalanine) added. Samples were centrifuged for 15 min (4°C, 13000 g) and the supernatant removed and blown dry under nitrogen gas. The sample was then re-solubilized with solution (acetonitrile: water = 1:1) and extracted by low-temperature ultrasonication for 5 min (5°C, 40 KHz), followed by centrifugation at 13000 g and 4°C for 10 min. The supernatant was transferred to sample vials for LC-MS analysis. A quality control sample (QC) was prepared by mixing equal volumes of 20 samples. The QC samples were used to monitor the stability of the analysis, and were handled in the same manner as the test samples. Metabolites were identified by searching the HMDB (http://www.hmdb.ca/), Metlin (https://metlin.scripps.edu/) and Majorbio databases. The data matrix obtained by searching databases was uploaded to the Majorbio cloud platform (https://cloud.majorbio.com) for data analysis [37]. The variables from QC samples with relative standard deviation (RSD) > 30% were excluded and log10 logarithmicized to obtain the final data matrix for subsequent analysis.
The R package “ropls” (Version 1.6.2) was used to perform PCA and orthogonal least partial squares discriminant analysis (OPLS-DA), and 7-cycle interactive validation evaluating the stability of the model. The metabolites with Variable importance in projection (VIP) > 1, P < 0.05 were deemed as significantly different based on the VIP obtained by the OPLS-DA model and the p-value generated by student’s t test.
Correlation analysis
The correlation analysis was conducted using the R package "psych" (version 1.8.12). The correlations between differential bacterial and metabolites, four types of data (the basic parameters, hepatic triglyceride, serum parameters, and ruminal fermentation indicators) and differential bacterial, as well as significantly differential metabolites and four types of data, were calculated via spearman's correlation coefficient. The correlation coefficients (r) ranged from − 1 to 1, with positive and negative values indicating positive and negative correlations. The absolute value of r indicated the strength of the correlation between variables. Specifically, r values of -1, 0, and 1 indicated a fully negative correlation, no correlation, and a fully positive correlation, respectively. Correlation significance was determined via P-values below 0.05 and 0.01, indicating significant and extremely significant correlations, respectively.
Cell isolation, culture and processing
Primary hepatocytes and pre-adipocytes were isolated using the collagenase (C0130 and V900893, Sigma-Aldrich, USA) digestion method from 5 Holstein calves (1 d old, female, 30–40 kg, fasting, rectal temperature: 38.7–39.7°C, healthy) purchased from a commercial dairy farm (Changchun, China). The culture procedure for both hepatocytes and pre-adipocytes and the method of pre-adipocyte maturation were the same as described previously [38, 39]. To determine the toxicological effect of inosine, hepatocytes and mature adipocytes were treated with different concentrations (0, 1, 2, 5, 10 µM) of inosine for 12 h and 24 h, respectively. To study the function of inosine in the metabolism of hepatocytes, cells were co-treated with 1.2 mM NEFA and different concentrations of inosine (0, 1, 2, 5, 10 µM) for 12 h. The method of NEFA preparation was described previously [40]. To activate the insulin signaling pathway, hepatocytes were treated with 100 nM bovine insulin (I6634, Sigma-Aldrich, USA) for 30 min before sampling. To investigate the effective receptor of inosine in hepatocyte metabolism, hepatocytes were pretreated separately with 100 nM SLV320 (A1R inhibitor, HY-19533, MCE, USA), 1 µM SCH58261 (A2aR inhibitor, HY-14858, MCE, USA), siRNA against A2bR or A3R, and subsequently incubated with inosine and NEFA.
To study the function of inosine in lipolysis of adipocytes, mature adipocytes were incubated with 10 µM ISO (S2566, Selleck, China) for 3 h and then treated with different concentrations (0, 1, 2, 5, 10 µM) of inosine for 12 h. To activate the insulin signaling pathway, adipocytes were treated with 100 nM insulin for 30 min before sampling. To investigate the function of inosine on inflammation in adipocytes, cells were co-treated with 1 µg/mL LPS (L2630, Sigma-Aldrich, USA) and inosine for 12 h. To investigate the effective receptor of inosine in adipocytes, adipocytes were pretreated separately with SLV320, SCH58261, siRNA against A2bR or A3R, and subsequently incubated with inosine and ISO or LPS.
Cell viability, triglyceride and glycerol content assessment
Cell viability was determined with the Cell Counting Kit (CCK)-8 assay according to the instructions (CK04, Dojindo, Japan). Liver tissue, hepatocytes and adipose were homogenized in lysis solution and then heated at 70℃ for 10 min. After cooling, the sample was vortexed and centrifuged at 2,000 × g for 5 min at room temperature. The supernatant was collected and analyzed using an enzymatic kit (E1013, Applygen, China) following the manufacturer’s instructions. Total protein concentration was measured with the BCA method (P1511, Applygen, China). The glycerol content of the adipocyte culture supernatant was measured using a glycerol measurement kit (E1002, Applygen, China) following the manufacturer’s instructions.
Nile red staining
A fluorescent quantitative method for determining Nile red-stained fat droplets was employed to determine the extent of triglyceride synthesis. Briefly, primary hepatocytes were fixed with 4% paraformaldehyde solution and stained with Nile red (1 µg/mL). After incubation for 30 min at 4°C, cells were washed with PBS and counterstained with 4′,6-diamidino-2-phenylindole (DAPI, Solarbio) for nuclear visualization. Then, cellular Nile red–stained lipid droplets were observed using a confocal scanning laser microscope (Leica).
Statistical analysis
GraphPad Prism 9.0 was used for the statistical analysis. All data were tested for normality of distribution with the Shapiro-Wilk test. Significant difference between two groups was evaluated by Mann-Whitney U test (non-parametric) or two-tailed unpaired Student’s t-test (parametric). For more than two groups of comparison, Kruskal-Walli’s test (non-parametric) or one-way analysis of variance (ANOVA) followed by Dunn’s test for non-parametric samples or Tukey test for parametric samples as a post hoc test. Data were expressed as the mean ± SEM, *P < 0.05, **P < 0.01. Other specific statistical analyses were mentioned in each methods section.