Skip to main content

Bilirubin ameliorates murine atherosclerosis through inhibiting cholesterol synthesis and reshaping the immune system

Abstract

Atherosclerosis is a chronic inflammatory disease caused mainly by lipid accumulation and excessive inflammatory immune response. Although the lipid-lowering and cardioprotective properties of bilirubin, as well as the negative relationship between bilirubin and atherosclerosis, were well documented, it is not yet clear whether bilirubin can attenuate atherosclerosis in vivo. In this study, we investigated the role of bilirubin in improving atherosclerosis. We found that mildly elevated bilirubin significantly reduced the risk factors of atherosclerosis, such as plasma glucose, total cholesterol, and low-density lipoprotein cholesterol, and the formation of atherosclerotic plaques, liver total cholesterol, and cholesterol ester concentration in apolipoprotein E-deficient (ApoE−/−) mice fed a western-type (high fat) diet. It was further found that bilirubin could promote the degradation of 3-Hydroxy-3-Methylglutaryl-CoA Reductase (HMGCR), a rate-limiting enzyme for endogenous cholesterol synthesis. Using mass cytometry-based high dimensional single cell analysis, we observed a decrease of natural killer cells and an increase of dendritic cells and myeloid-derived suppressor cells, which all are closely associated with atherosclerosis risk factors and contribute to the improvement of atherosclerosis, in ApoE−/− mice treated with bilirubin. By in-depth analysis, modulation of multiple spleen or peripheral blood T cell clusters exhibiting either positive or negative correlations with total cholesterol or low-density lipoprotein cholesterol was detected after bilirubin treatment. In this study, we demonstrate that bilirubin serves as a negative regulator of atherosclerosis and reduces atherosclerosis by inhibiting cholesterol synthesis and modulating the immune system.

Introduction

Cardiovascular disease (CVD) has become the leading cause of death and morbidity worldwide [1]. Atherosclerosis is the principal cause of CVDs, cerebral infarctions, and peripheral vascular diseases. Atherosclerosis is defined as a chronic inflammatory disease of the arteries. Plaque decomposition or plaque rupture can block arteries, causing severe acute cardiovascular events such as myocardial infarction, stroke, and sudden death [2,3,4,5].

Bilirubin is a final product of heme catabolism. In microsomes, heme is converted into carbon dioxide, iron, and biliverdin IX (biliverdin) and thereafter biliverdin reductase transforms biliverdin into bilirubin [6]. The physiological concentration of bilirubin is beneficial for oxidative stress-related pathology, a key event in the development of atherosclerosis [7,8,9]. Bilirubin levels have been shown to be inversely correlated with CVD mortality and severity [10,11,12,13,14,15,16]. Numerous studies have shown that bilirubin offers cardioprotective properties in Gilbert syndrome (GS) patients, who suffer sustained unconjugated hyperbilirubinemia. A low prevalence (2%) of ischemic heart disease (IHD) was found in GS patients, as compared to a higher prevalence (12.1%) among the general population [17]. Hyperbilirubinemia in GS subjects or Gunn rats causes a decrease in circulating cholesterol, triacylglycerol, and oxidized low-density lipoprotein (LDL) by increasing the serum anti-oxidative functions [18, 19]. Bilirubin has the potential to lower lipid, which is associated with cardiovascular protection in GS, as well as dyslipidemia, one of the leading risk factors for atherosclerosis [20]. Oxidative stress plays a critical role in the pathogenesis of atherosclerosis. Bilirubin, a natural endogenous antioxidant, protects vascular cells from oxygen radical damage, inhibits oxidative modification of LDL and improves lipoprotein composition [21, 22]. Bilirubin also inhibits inflammation by downregulating the expression of endothelial cell adhesion molecules, preventing leukocyte adhesion, rolling, and infiltration into vessels [23,24,25]. Additionally, bilirubin reduces neoplastic intima formation by inhibiting smooth muscle cell proliferation and migration. Because of these properties, bilirubin is involved in every aspects of atherosclerosis, inhibiting plaque formation and promoting plaque stabilization [26, 27]. These findings greatly demonstrate the importance of bilirubin in the regulation of lipids and improvement of atherosclerosis. To date, the underlying mechanisms of hypocholesterolemia in hyperbilirubinemia have not been fully elucidated.

Atherosclerosis is a disease associated with inflammation that relies heavily on the immune system for its development and modulation [28]. New high-throughput single cell technologies, such as mass cytometry and single-cell RNA sequencing, enable a comprehensive analysis of the immune system during atherosclerosis progression through uncovering the diversity of > 15 immune cell populations [29,30,31]. The dissection of immune modulation has greatly contributed to the understanding of atherosclerosis mechanisms and has provided novel insights into its treatment via immune therapy. The hypolipidemic effect of bilirubin led us to examine its effect on the development of atherosclerosis in an apolipoprotein E-deficient (ApoE−/−) mouse model, as well as the mechanisms involved from the perspective of lipid regulation and immunomodulation using mass cytometry-based single cell analysis. In the present study, we found that bilirubin negatively regulates atherosclerosis through inhibiting cholesterol synthesis and modulating the immune system.

Methods

Experimental animals

Animal experiments in our study were approved by the Institutional Committee for the Use and Care of Laboratory Animals of Guangzhou Medical University (2017-014) and performed in accordance with the guidelines from Directive 2010/63/EU of the European Parliament on the protection of animals used for scientific purposes or the NIH guidelines. 8-week-old male ApoE−/− mice (C57BL/6.129P2-APOE/J, Beijing Vital River Laboratory Animal Technology Co., Ltd) were divided into two subgroups. One group (Bilirubin, n = 14) received an intraperitoneal injection of bilirubin (20 mg/kg/d, B4126, Sigma) every other day for 10 weeks. Mice injected with the vehicle served as a control (Control, n = 8). Both groups started the western-type diet supplemented with 21% fat (wt/wt) and 0.15% cholesterol (wt/wt), simultaneously receiving the injection of bilirubin or vehicle. All mice were maintained on these diets for 10 weeks. The body weight, diet, and drinking of mice were monitored for a continuous week. The mice were euthanized with an overdose of pentobarbital sodium (200 mg/Kg) for tissue collection.

Biochemical assays

Peripheral blood was collected into a heparin-coated tube by cardiac puncture and centrifuged to isolate plasma. Total bilirubin (T-BIL) in the plasma was determined by the vanadic acid oxidation method. Liver enzymes, including alanine aminotransferase (ALT) and aspartate aminotransferase (AST), were measured by the ultra-violet lactate and malate dehydrogenase methods. Total cholesterol and triacylglycerol levels were measured by an oxidase method. High-density cholesterol (HDL-C) and low-density cholesterol were measured by a direct method. The total bilirubin assay kit was from FUJIFILM Wako pure Chemical Corporation, Japan. The kits for glucose, liver enzymes, lipids, and lipoproteins were from Kehua Bio-Engineering, China. All the spectrophotometric methods were performed according to the manufacture’s instructions using an autoanalyzer.

Histological and morphometric analysis

The heart was perfused with phosphate-buffered saline. Whole aorta from the heart, extending 5–10 mm after bifurcation of the iliac, including the brachiocephalic trunk, left common carotid artery, and subclavian arteries, were dissected free of fat tissues and opened longitudinally before being fixed in formalin solution at room temperature (RT) for 24 h. After staining with Oil Red O (ORO), aortae were pinned on a black wax plate and photographed under stereoscopy at standardized magnification and illumination. The ORO positive area and entire aorta tree area were quantified using Image J software. Aorta lesion areas were expressed as a percentage of the whole aorta tree area. The aortic root and brachiocephalic trunk were serially sectioned into 5 μm cryosections. The first section of the aortic sinus was harvested and collected in pairs on glass slides when all three aortic valves became visible in the lumen of the aorta. Neutral lipids deposition in the lesion was visualized by hematoxylin and eosin (H&E) or ORO staining, respectively. The collagen of the lesion was detected by Masson's Trichrome staining. The liver was dissected and weighed, expressed as a ratio of liver weight to body weight. Liver sections were collected and performed H&E or ORO staining for histological and lipid deposition analysis.

Immunofluorescence staining

The contents of macrophages and smooth muscle cells in aortic sinus lesions were determined by staining with anti–Galectin 3(ab53082) and anti-α-SMA (ab5694) and detecting with Alexa-conjugated antibodies (ab150073 and ab150080). Nuclei were counterstained with DAPI. All antibodies were purchased from Abcam.

Total cholesterol, cholesterol ester, and free cholesterol assay

Hepatic total cholesterol and free cholesterol concentration were measured using a tissue total cholesterol assay kit (E1015; Applygen Technologies, Inc.) and free cholesterol assay kit (E1016; Applygen Technologies, Inc.) following the manufacturer's instructions. Cholesteryl ester was quantified by subtracting the free cholesterol values from the total cholesterol value.

Cell culture and treatment

Human normal hepatocyte LO2 (ATCC) was cultured in Dulbecco's Modified Eagle's Medium (DMEM) containing 10% fetal bovine serum (FBS; Gibco), 1 mM L-glutamine, 100 U/ml penicillin, and 100 μg/ml streptomycin. Hepatocellular carcinoma cell line HepG2 (ATCC) was cultured in RPMI-1640 medium (Gibco, Carlsbad, CA, USA) supplemented with 10% FBS, 1 mM L-glutamine, and 100 U/ml penicillin, and 100 μg/ml streptomycin. Both LO2 and HepG2 were maintained at 37 °C in a humidified atmosphere of 5% CO2. For the dose-dependent assay, cells subcultured in the 6-well plate were treated with various concentrations of bilirubin (0, 3, 6 and 12 μM) for 24 h. For the time-dependent studies, cells were treated with 12 μM bilirubin for 0, 6, 12 and 24 h.

Real-time PCR

Total RNA from cells was extracted with TRIzol reagent (Invitrogen) and reverse-transcribed using PrimeScript™ RT reagent Kit with gDNA Eraser (Takara) according to the manufacturer’s instructions. Real-time PCR was performed on a StemOne™ Real-Tim PCR System (Applied Biosystem, USA) following TB Green™ Premix Ex Taq ™ (Tli RNaseH Plus, Takara) protocol. The primer pairs were used: 18S, forward 5′-CCCAGTAAGTGCGGGTCATAA-3′, reverse 5′- CCGAGGGCCTCACTAAACC-3′; mouse HMGCR forward 5′-TGATTGACCTTTCCAGAGCAAG-3′, reverse 5′- CTAAAATTGCCATTCCACGAGC-3′. Gene expression was calculated based on the −ΔΔCt method. The relative level of mRNA expression was determined using as reference the geometric mean of human/mice 18S (ThermoFisher Scientific).

Western blot

Cells or tissues were collected and homogenized in ice-cold RIPA buffer. Lysates were maintained constant agitation for 2 h at 4 °C. Gently aspirate the supernatant to a prechilled tube after the centrifugation at 13,200 rpm for 10 min at 4 °C. The protein level of lysates was determined by the BCA protein assay. An aliquot of proteins was separated by SDS-PAGE and then transferred to the immunoblot PVDF membrane. The immunoblots were incubated with primary antibodies at 4 °C overnight, followed by HRP-conjugated secondary antibodies (1:5000) at RT for 1 h. The signals of targeted proteins were visualized by employing the ECL detection system (Amersham). Primary antibodies used for immunoblots were rabbit polyclonal antibodies against HMGCR (ab174830, diluted 1:1000) and mouse monoclonal antibody against GAPDH (ab8245, diluted 1:1000). All antibodies were purchased from Abcam.

Co-immunoprecipitation (Co-IP)

Cells were collected and treated with IP lysis buffer containing protease inhibitor cocktails at 4 °C for 30 min. Cell lysates were centrifugated at 13,200 rpm for 10 min at 4 °C to remove cell aggregates. 50 μl of protein extract were dissolved in a DTT-containing loading buffer as an input sample. The remaining protein extract was adjusted to a concentration of 2 μg/μl of Co-IP dilution buffer. 500 μl of protein extract were immunoprecipitated with magnetic beads conjugated with indicated antibodies at RT for 1–2 h in a roller. The beads were washed with prechilled Co-IP washing buffer three times and then heated in 2 × loading buffer for 10 min at 70 °C. After centrifugation at 12,000 rpm for 10 min at 4 °C, the beads were collected by inserting the protein samples in a magnetic stack for 10 min at RT. The well-prepared protein samples were analyzed by western blot. The immunoblots were probed with rabbit polyclonal antibody against HMGCR (ab174830, diluted 1:1000), rabbit polyclonal antibody against ubiquitin (ab179434, diluted 1:5000), and mouse monoclonal antibody against GAPDH (ab8245, diluted 1:1000). All antibodies were purchased from Abcam.

Protein turn-over assay

To evaluate protein half-life, we treated LO2 and HepG2 cells with or without bilirubin (12 μM) in the presence of cycloheximide at a final concentration of 50 µg/ml to block protein synthesis for the indicated time (0, 6, 12, and 24 h). Protein samples were subjected to western blot to analyze the expression of HMGCR.

Peripheral blood and spleen cell preparation

The peripheral blood was collected from the heart after the injection of pentobarbital sodium. After mice were sacrificed, the spleen was isolated and spleen cells suspension was separated by gently crushing the spleen. Cells from peripheral blood or spleen were fixed with Fix I buffer (Fluidigm) for 10 min at RT. After washing with cold PBS, red blood cell lysis buffer was used to remove erythrocytes from peripheral blood or spleen cells.

Cell immunostaining

Peripheral blood or spleen cells were barcoded separately with premade combinations of six different palladium isotopes using the Cell-ID 20-Plex Pd Barcoding Kit (Fluidigm) for 30 min at RT. The samples were washed twice with cell staining buffer and then pooled together. The pooled samples were incubated with anti-CD16/32 (FcR III/II, Biolegend, CA, USA) for 10 min at RT to lower nonspecific antibody binding. These cells were washed twice with cell staining buffer and stained with a cocktail of 26 metal isotope-conjugated antibodies (Additional file 1: Table S1) for 30 min at RT. Stained cells were washed and stained with Intercalator-Ir (Fluidigm) at 4 °C overnight. After washing with cell staining buffer and ddH2O, cells were resuspended in ultrapure water supplemented with EQ Four Element Calibration Beads (Fluidigm) and analyzed using Helios mass cytometer (Fluidigm).

Mass cytometry data processing and analysis

The resulting flow cytometry standard (FCS) files generated by mass cytometry were normalized, randomized and debarcoded using the CyTOF software (Fluidigm). The debarcoded files were uploaded to cytobank.cn, an online platform for cytometry-based single cell analysis. The single cell data without debris and doublets were used for high-dimensional analyses, such as viSNE and FlowSOM.

Statistical analysis

The Shapiro–Wilk test was used to determine whether the data was normally distributed. A parametric test (Student’s t test) was used to determine the statistical significance between 2 groups if the data is normally distributed, otherwise, a nonparametric test (Mann–Whitney test) was used. Correlation analyses were performed using Person’s correlation coefficient. Error bars represent the mean ± standard deviation (SD). p < 0.05 was regarded as statistically significant.

Results

Bilirubin lowers the risk of atherogenesis in ApoE−/− mice

To investigate whether bilirubin can protect atherogenesis in vivo, we used a murine model in which ApoE−/− mice were fed on a western-type (high fat) diet. Mice received an intraperitoneal injection of bilirubin or vehicle every other day for 10 weeks. There was no significant difference between the average body weight of the two groups of mice throughout the 10 weeks of intervention. (Fig. 1A). In comparison with a vehicle control group, exogenous bilirubin administered resulted in a fourfold increase in total plasma bilirubin (Fig. 1B). Treatment with bilirubin reduced plasma glucose, total cholesterol, and LDL cholesterol levels in mice. (Fig. 1C and D), whereas both groups had similar levels of triglycerides and high-density lipoprotein (HDL) cholesterol. Furthermore, multiple correlations were found between glucose, total cholesterol, LDL, HDL, triglyceride, and total bilirubin. (Fig. 1E). For example, total bilirubin levels were negatively and significantly correlated with the levels of glucose and LDL levels in the peripheral blood. The LDL level was positively and significantly correlated with the total cholesterol level. Mice with higher levels of glucose also had higher levels of total cholesterol and LDL. There were significant and positive correlations between the concentration of triglyceride and total cholesterol or HDL (Fig. 1E). Considering the long-term exposure, we also investigated the bilirubin toxicity to the liver through histological and biochemical analysis. Bilirubin-treated mice displayed a downward trend in a lower hepatic enzyme activity (AST and ALT) than control mice (Additional file 1: Fig. S1B). But both groups showed no change in liver weight/body weight. (Additional file 1: Fig. S1A). These findings illustrate that a mildly elevated bilirubin reduced the risks of developing atherosclerosis without causing liver toxicity, suggesting a protective role of bilirubin in atherogenesis.

Fig. 1
figure 1

Bilirubin lowers the risk of atherogenesis in ApoE−/− mice. 22 male ApoE−/− mice (8 weeks) fed with a western-type diet were received the intraperitoneal injection of bilirubin (bilirubin group, 20 mg/Kg, n = 14) and vehicle (control group, n = 8) once every other day for ten weeks. A The body weight of mice was monitored weekly. B–D After 10 weeks of treatment, peripheral blood plasma was collected and the levels of total bilirubin, glucose, total cholesterol (TCHO), low-density lipoprotein cholesterol (LDL), high-density lipoprotein cholesterol (HDL), and total triglyceride (TG) were measured. E A heatmap showing the Pearson correlation coefficients for relationships between the concentrations of total bilirubin (TBIL), glucose (GLU), LDL, HDL, TG, and TCHO in peripheral blood. Dot plots (n = 22) showing the Pearson correlation coefficients for relationships between the concentrations of total bilirubin (TBIL), glucose (GLU), LDL, HDL, TG, and TCHO. Error bars represent mean ± standard deviation (SD). Dots represent individual samples, n (control) = 8, n (bilirubin) = 14. *P < 0.05, **P < 0.01 and ***P < 0.001. Differences were determined by an unpaired t-test. Correlations were determined by a Pearson test

Bilirubin reduces atherosclerosis in ApoE−/− mice

In vivo, bilirubin demonstrated a risk-lowering effect, thus we investigated whether it affects atherosclerotic lesion formation. Morphometric assays were used to measure the aortae and the aortic sinus lesion formation after 10 weeks of intervention. In both groups of mice, lesions were present throughout the aorta, but lesions were reduced in mice treated with bilirubin, as determined by en face ORO staining. (Fig. 2A). The bilirubin treatment significantly decreased the size and lipid load of the aortic sinus plaque and the aortic sinus lesion area. (Fig. 2B). Mice treated with bilirubin also showed reduced brachiocephalic arteries (Fig. 2C). Additionally, bilirubin treatment promoted collagen deposition (Fig. 2B) and α-SMA ( +) smooth muscle cell accumulation without affecting the distribution of Galectin-3( +) macrophages (Fig. 2D). The liver plays a crucial role in lipid metabolism. In the liver, cholesterol is derived from dietary, de novo synthesis, and transport from extrahepatic tissues. In mice fed a western diet for a long time, excessive lipid accumulates in the liver, increasing the risk of atherosclerosis and liver steatosis. H&E and ORO staining showed that bilirubin-treated mice had less hepatic steatosis than vehicle-treated mice (Fig. 2E). Next, The liver was then analyzed using an enzymatic method to determine the content of total cholesterol, cholesterol ester, and free cholesterol. The bilirubin treatment decreased liver total cholesterol and cholesterol ester concentration, but did not affect free cholesterol (Fig. 2F). As a result, bilirubin significantly inhibits the formation of atherosclerotic plaques and favors the phenotype of stable plaque.

Fig. 2
figure 2

Bilirubin inhibits atherosclerotic lesions formation and hepatic lipid accumulation in ApoE−/− mice. ApoE−/− mice fed with a western-type diet were intraperitoneally injected with (bilirubin, n = 14) or without (control, n = 8) bilirubin. After ten weeks, the animals were sacrificed and the aorta, heart, brachiocephalic artery, and liver tissue were collected for morphometry. A Aortas fixed with formalin was performed en face Oil Red O (ORO) staining and the lesion area was quantified. B Series of frozen sections of aortic sinus were subjected to hematoxylin–eosin (H&E), ORO, and Masson staining. Aortic sinus lesion area was quantified in the context of H&E staining. C Histological analysis of brachiocephalic arteries after H&E and ORO staining. The neointima/media ratio was calculated. D Immunofluorescence staining was applied to analyze the distribution of Galectin-3 ( +) (macrophages, red) and α-SMA ( +) (smooth muscle cells, green) in aortic valve lesions. The area of Galectin-3 ( +) and α-SMA ( +) was quantified. E Histological analysis of liver after H&E and ORO staining. Area with oil red color are steatosis area. The black arrows indicate the representative steatosis area. The red arrow indicates the normal area. F The amount of total cholesterol and free cholesterol in liver lysates was measured. Cholesteryl esters were calculated by subtracting the free cholesterol values from the total cholesterol value. Error bars represent mean ± SD. n (control) = 8, n (bilirubin) = 14. *P < 0.05, **P < 0.01 and ***P < 0.001. Differences were determined by an unpaired t-test

Bilirubin promotes HMGCR protein degradation

The liver is the primary organ that synthesizes cholesterol and around 20–25% of total daily cholesterol production occurs here. Since HMGCR is the key rate-limiting enzyme for endogenous cholesterol synthesis, we examined HMGCR expression in vivo and in vitro. HMGCR expression in the liver was significantly reduced in ApoE−/− mice treated with bilirubin compared to vehicle-treated mice (Fig. 3A and Additional file 1: S2A). In some cases, two bands of HMGCR have been detected in our and other studies [32, 33]. The upper one is the HMGCR band and the lower one is the non-specific band. The expression of HMGCR in the liver was also decreased in BALB/c mice treated with bilirubin for three consecutive days (Additional file 1: Fig. S2B). In two hepatic cell lines, bilirubin did not decrease the mRNA expression of HMGCR in vitro (Fig. 3B), but decreased the protein level of HMGCR in a time- and dose-dependent manner (Fig. 3C and D), suggesting that bilirubin may affect HMGCR expression at the translational or post-translational level. After treating two hepatic cell lines with cycloheximide (CHX), we found that bilirubin treatment decreased protein synthesis (Fig. 3E). Considering that HMGCR is degraded through ubiquitination at lysine 89 and 248 [34], we further evaluated the effect of bilirubin on the ubiquitination of HMGCR. Bilirubin increased HMGCR ubiquitination in HepG2 hepatic cells treated with MG132, a proteasome inhibitor (Fig. 3F). Furthermore, bortezomib, a potent protease inhibitor, reversed the decrease in HMGCR levels induced by bilirubin in two hepatic cell lines (Fig. 3G). These results indicate that bilirubin reduces HMGCR level mainly by promoting its ubiquitination and degradation, consequently reducing endogenous cholesterol synthesis.

Fig. 3
figure 3

Bilirubin downregulates HMGCR expression in liver and human hepatocytes by promoting the ubiquitination and degradation of HMGCR. A Western blot analysis of hepatic HMGCR expression in ApoE−/− mice treated with bilirubin (n = 3) or vehicle (n = 3). B Human normal hepatocytes (LO2 or hepatoma cell line (HepG2) were treated with different concentrations of bilirubin (BIL, 0, 3, 6, and 12 μM) for 24 h and the mRNA expression of HMGCR was determined by RT-qPCR. C LO2 or HepG2 cells were treated with different concentrations of bilirubin (0, 3, 6, and 12 μM) for 24 h and the protein level of HMGCR was measured by Western blot. D LO2 or HepG2 cells were treated with 12 μM of bilirubin for 0, 6, 12, or 24 h and the protein level of HMGCR was measured by Western blot. E HepG2 or LO2 cells were incubated with BIL (12 μM) or vehicle in the presence of 50 μg/ml of cycloheximide for the indicated time and the protein level of HMGCR was measured by Western blot. F HepG2 cells were treated with or without BIL (12 μM) in the presence of protease inhibitor MG132 (10 μM) for 24 h. Cell lysates were immunoprecipitated with anti-HMGCR antibody and then probed with an antibody against ubiquitin or HMGCR, respectively. G LO2 or HepG2 cells were incubated with or without BIL (12 μM) in the presence or absence of protease inhibitor bortezomib (Btz, 50 nM) for 24 h. The expression of HMGCR was determined by Western blot. Representative western blot images of three independent experiments (n = 3) with similar results are shown. Asterisks indicate the correct HMGCR bands. All the bands were quantified and the expression of HMGCR was normalized to GAPDH. Normalized values are below the bands

Bilirubin affects immune system in the spleen

Since bilirubin significantly improves atherosclerosis which is closely associated with the immune response, it is likely that bilirubin also has an effect on the immune system. Thus, we used a high-throughput single cell technology, mass cytometry [31, 35], to simultaneously measure the expression of 26 markers on immune cells from the spleen or peripheral blood at the single cell level. ViSNE analysis of CD45+ cells from the spleen was performed to visualize these high-dimensional mass cytometry data in two dimensions. According to the expression of 13 surface markers on the cells displayed on the viSNE map (Fig. 4A), 18 major immune cell populations were identified and gated (Fig. 4B). In each gated population, the expression of 13 markers was the same as the phenotypes of the indicated types of immune cells (Fig. 4C). Next, we compared the proportions of all these cell populations and found that after bilirubin treatment, the percentages of CD11blow dendritic cells (CD11blowDCs) significantly increased, while the proportion CD11bhighDCs markedly decreased (Fig. 4D). In addition, the percentage of CD4 T cells in the spleen T cells was significantly reduced, while the proportion of CD8 T cells increased (Fig. 4D), suggesting that bilirubin induced a change in spleen T cell composition. Interestingly, the proportion of CD4 T cells among spleen T cells was strongly, significantly, and positively related to the LDL and (total cholesterol) TCHO concentration in the blood, while CD8T was strongly, significantly, and negatively related to the LDL and TCHO concentration in blood bilirubin (Fig. 4E), indicating that spleen CD4 T cells are positively associated with atherosclerosis. The proportion of CD11bhighDC that decreased by bilirubin showed a strong and negatively relationship with total bilirubin (Fig. 4E). Moreover, the proportion of CD8 T cells in spleen T cells, as well as the percentages of CD11blowDC, were significantly and negatively associated with the LDL and TCHO concentrations in the blood (Fig. 4E), implicating that these immune cells may favor the improvement of atherosclerosis.

Fig. 4
figure 4

Bilirubin affects immune cells in the spleen. ApoE−/− mice fed with a western-type diet were intraperitoneally injected with (bilirubin, n = 14) or without (control, n = 8) bilirubin. The animals were sacrificed after ten weeks. The spleen cells were collected and simultaneously stained with 26 metal isotope-labeled antibodies. The expression of 26 surface markers was then determined by mass cytometry. A viSNE was used to visualize the distribution of the spleen CD45+ immune cells from both groups. Cells on the viSNE map were colored by normalized expression of indicated surface markers. B 18 cell populations were identified and colored on the viSNE map. C A heatmap showing the normalized expression of 13 indicated markers in 18 cell populations. D Bar plots showing the frequencies of indicated populations in spleen CD45+ or T cells obtained from ApoE−/− mice treated with or without bilirubin. E Dot plots (n = 22) showing the Pearson correlation coefficients for relationships between the concentrations of total bilirubin LDL, or TCHO in peripheral blood and the frequencies of indicated cell populations in the spleen. Error bars represent mean ± SD. Dots represent individual samples, n (control) = 8, n (bilirubin) = 14. *P < 0.05, **P < 0.01. Differences were determined by an unpaired t-test. Correlations were determined by a Pearson test

Bilirubin alters the composition of T cells in the spleen

T cells can influence the process of atherosclerosis by regulating their anti-inflammatory and pro-inflammatory functions and their functions are mainly determined by the direction of their differentiation. T cells are highly heterogeneous and contain various subtypes of differentiated T cells. Thus, we used 12 T cell-related markers to systematically analyze the spleen T cell. After viSNE analysis, heterogeneity in markers was assessed at the single-cell level. From the viSNE map of T cells, we also observed a huge heterogeneity of the T-cell compartments (Additional file 1: Fig. S3A). FlowSOM, an unsupervised cluster analysis for visualizing and interpreting cytometry data, was used on these T cells to map T cell phenotypes exhaustively and to further define any exclusive or significantly altered T cell clusters. 30 T cell clusters containing cells with similar phenotypes were identified and visualized on the viSNE map (Fig. 5A). The expression profiles of all the T cell clusters were visualized in a heatmap (Fig. 5B) and this approach led to the identification of 12 CD4+ clusters, one double-positive cluster, eight CD8+ clusters, five double-negative phenotypes, one natural killer (NK) T cluster and three TCRβ− clusters (Fig. 5B). Among these 30 T cell clusters, proportions of T5 (CD44+CD62L−CXCR3+CD4T), which mainly contains effector memory (em) Th1 cells, and T16 (Ly6C+CD44+CD62L+CD8T), which mainly contains central memory Ly6C+CD8 T (Tcm) cells, were significantly increased after bilirubin treatment (Fig. 5C). In addition, these two clusters, together with T23 which was increased with a trend close to significance, had a mutual, significant, and negative relationship with the concentrations of LDL and TCHO in the blood (Fig. 5D and E). Furthermore, T6 and T13 proportions decreased with a trend close to significance after bilirubin treatment (Fig. 5C). Bilirubin did not affect the proportion of T1 (CD3+CD4+CD25+CD127−Ly6C+) and T2 (CD3+CD4+CD25+CD127−Ly6C−) clusters which mainly contain regulatory T (Treg) cells (Additional file 1: Fig. S3B). These data indicate that bilirubin also can regulate the composition of spleen T cells which are related to the improvement of atherosclerosis.

Fig. 5
figure 5

Bilirubin alters the composition of T cells in the spleen. After determining the expression of 26 surface markers on the spleen cells using mass cytometry, the T cells expressing CD3 were gated and analyzed using FlowSOM. 30 T cell clusters were identified using FlowSOM. A viSNE was used to visualize the distribution of the spleen T cells. 30 T cell clusters were identified and colored on the viSNE map. B A heatmap showing the normalized expression of 12 indicated markers in 30 T cell clusters. C Bar plots showing the frequencies of T5, T16, T6, T13, or T23 in spleen T cells obtained from ApoE−/− mice treated with or without bilirubin. D A heatmap showing the Pearson correlation coefficients for relationships between the concentrations of total bilirubin (TBIL), LDL, HDL, TG, or TCHO in peripheral blood and the frequencies of indicated spleen T cell clusters. E Dot plots (n = 22) showing the Pearson correlation coefficients for relationships between the concentrations of LDL or TCHO in peripheral blood and the frequencies of T5, T16, or T23 clusters. Error bars represent mean ± SD. Dots represent individual samples, n (control) = 8, n (bilirubin) = 14. *P < 0.05. Differences were determined by an unpaired t-test. Correlations were determined by a Pearson test

Bilirubin regulates myeloid-derived suppressor cells (MDSCs), DCs and NK cells in the peripheral blood

After bilirubin administration, the immune cells in the peripheral blood were also comprehensively analyzed using mass cytometry. viSNE was also used to identify the cell populations (Fig. 6A). 19 immune cell populations were obtained and the expression of their phenotypic markers was confirmed (Fig. 6B and C). The percentages of these 19 cell populations were then compared. Proportions of DC, Ly6C+NK, and all NK cells were significantly reduced in the peripheral blood of mice treated with bilirubin (Fig. 6D). The percentages of polymorphonuclear (PMN)-MDSC, neutrophils, and all MDSC significantly increased (Fig. 6D) and were strongly, significantly, and positively correlated with bilirubin concentration (Fig. 6E and F), implicating that bilirubin treatment may directly increase these cells. In addition, the proportion of all MDSC was significantly and negatively correlated with the concentrations of LDL and TCHO in the blood (Additional file 1: Fig. S4A). Considering that higher levels of LDL and TCHO are closely associated with atherogenesis, MDSCs increase by bilirubin may contribute to the improvement of atherogenesis though regulating the levels of LDL and TCHO. Moreover, proportions of Ly6C+NK and all NK cells showed a significant negative correlation with blood bilirubin level (Fig. 6G), suggesting that the reduction of these NK cells may be directly caused by bilirubin administration. The proportion of Ly6C+NK cells correlated positively and significantly correlated with the level of LDL in the blood (Additional file 1: Fig. S4B). Furthermore, B cells were negatively correlated with the concentration of bilirubin and LDL, and positively associated with TCHO concentration (Fig. 6F and Additional file 1: Fig. S4C). These data propose that bilirubin affects the MDSC and NK cells which are closely associated with the improved atherosclerotic blood lipid profiles.

Fig. 6
figure 6

Bilirubin reshapes the immune system in the peripheral blood. ApoE−/− mice fed with a western-type diet were intraperitoneally injected with (bilirubin, n = 14) or without (control, n = 8) bilirubin. The animals were sacrificed after ten weeks. The peripheral white blood (PWB) cells were collected and simultaneously stained with 26 metal isotope-labeled antibodies. The expression of 26 surface markers was then determined by mass cytometry. A viSNE was used to visualize the distribution of the PWB cells from both groups. Cells on the viSNE map were colored by normalized expression of indicated surface markers. B 19 cell populations were identified and colored on the viSNE map. C A heatmap showing the normalized expression of 13 indicated markers in 19 cell populations. D Bar plots showing the frequencies of indicated populations in PWB cells obtained from ApoE−/− mice treated with or without bilirubin. E A heatmap showing the Pearson correlation coefficients for relationships between the concentrations of TBIL, LDL, HDL, TG, or TCHO in peripheral blood and the frequencies of indicated cell populations. F Dot plots (n = 22) showing the Pearson correlation coefficients for relationships between the concentrations of total bilirubin in peripheral blood and the frequencies of indicated cell populations in the PWB cells. G Dot plots showing the Pearson correlation coefficients for relationships between the concentrations of total bilirubin in peripheral blood and the frequencies of NK cells. Error bars represent mean ± SD. Dots represent individual samples, n (control) = 8, n (bilirubin) = 14. *P < 0.05, **P < 0.01. Differences were determined by an unpaired t-test. Correlations were determined by a Pearson test

Bilirubin mildly affects T cell composition in the peripheral blood

Peripheral blood T cell composition after bilirubin treatment was also analyzed using viSNE (Additional file 1: Fig. S5A). 29 T cell clusters, including 13 CD4+ clusters, eight CD8+ clusters, four NKT clusters, two double-negative phenotypes, and two TCRβ− clusters were identified (Fig. 7A and B). Treg cell clusters, T1 and T2, were not changed by bilirubin treatment (Additional file 1: Fig. S5B). T7 (CD44+CCR4+CD4 T), T9 (CD44+CCR4+CCR6+CD4 T), T13 (CD44−CD62L+CD4 T), and T17 (Ly6C+CD44−CD62L+CD8 T) were significantly decreased, whereas T20 (Ly6C−CD44−CD62L−CD8 T) was increased after bilirubin treatment (Fig. 7C). Additionally, proportions of two TCRβ− clusters, T28 and T29, were also changed by bilirubin with a trend close to significance (Additional file 1: Fig. S5C). The percentages of T9 and T28 showed a negative and significant correlation with the total bilirubin level, whereas T20 and T29 showed a positive correlation with total bilirubin (Fig. 7D and E and Additional file 1: S5C). There was a positive and significant correlation between the proportions of T7, T9, T13, and T17 and the amount of LDL in the blood. (Fig. 7F). Except for T17, the other clusters were also positively correlated with TCHO concentration (Additional file 1: Fig. S5D). These results suggest that bilirubin can induce changes of T cell subclusters and that these changed T cells also have a strong correlation with the improvement of atherosclerosis.

Fig. 7
figure 7

Bilirubin modifies the composition of T cells in the peripheral blood. After determining the expression of 26 surface markers on the PWB cells using mass cytometry, the T cells expressing CD3 were gated and analyzed using FlowSOM. 29 T cell clusters were identified. A viSNE was used to visualize the distribution of the peripheral blood T cells. 29 T cell clusters were identified and colored on the viSNE map. B A heatmap showing the normalized expression of 12 indicated markers in 29 T cell clusters. C Bar plots showing the frequencies of T7, T9, T13, T17, or T20 in peripheral blood T cells obtained from ApoE−/− mice treated with or without bilirubin (left panel). D A heatmap showing the Pearson correlation coefficients for relationships between the concentrations of TBIL, LDL, HDL, TG, or TCHO in peripheral blood and the frequencies of peripheral blood T cell clusters. E Dot plots (n = 22) showing the Pearson correlation coefficients for relationships between the concentrations of total bilirubin in peripheral blood and the frequencies of T9 or T20 in peripheral blood T cells. F Dot plots showing the Pearson correlation coefficients for relationships between the concentrations of LDL in peripheral blood and the frequencies of T7, T9, T13, or T17 in peripheral blood T cells. Error bars represent mean ± SD. Dots represent individual samples, n (control) = 8, n (bilirubin) = 14. *P < 0.05, **P < 0.01. Differences were determined by an unpaired t-test. Correlations were determined by a Pearson test

Discussion

Accumulating epidemiological evidence has demonstrated the negative relationship between bilirubin and CVD. Bilirubin, a potential predictor and protective marker of CVD, might become a new therapeutic target in the clinical setting. However, little is known about the direct cardioprotective effects of bilirubin. In the present studies, exogenous bilirubin inhibited lesion formation and lowered circulating cholesterol, LDL cholesterol, glucose concentrations, and hepatic lipid accumulation in ApoE−/− mice receiving western-type diet. We confirmed that bilirubin prompted HMGCR degradation via the proteasome pathway, suggesting that the cardioprotective effects of bilirubin may also be governed by its lipid-lowering properties. Furthermore, bilirubin significantly modulated the immune system, which resulted in a cardioprotective immune profile. These findings highlight the negative regulatory role of bilirubin in atherosclerosis while suggesting bilirubin as a potential treatment option for atherosclerosis.

Dyslipidemia and diabetes are common risk factors for atherosclerosis. A higher level of circulating LDL cholesterol or total cholesterol is independently and positively correlated with the prevalence of atherosclerotic peripheral vessel diseases [36]. In our study, there was a significant reduction of total cholesterol (− 17%), LDL-C (− 28%) and blood glucose (− 26%) in mice treated with bilirubin, confirming the lipid-lowering properties of bilirubin. Studies in animals (Gunn rats) or humans (GS subjects) of hyperbilirubinemia have consistently shown that bilirubin protected against IHD and atherosclerosis in dyslipidemia or diabetes high-risk population [17, 19, 37, 38]. In light of this circumstantial evidence, along with our finding showing that bilirubin attenuates atherosclerotic lesion formation and promotes a stable atherosclerotic plaque, further confirm that bilirubin is a negative regulator in atherosclerosis.

Bilirubin protects against atherosclerosis due to its anti-oxidative and anti-inflammatory properties, as well as the regulation of leukocyte migration mediated by endothelial adhesive molecules [7]. In atherosclerosis, foam cells are formed by the uptake of oxidatively modified lipids as oxLDL by macrophages [39]. Bilirubin has been shown to inhibit lipoprotein oxidation and therefore protects endothelial cells from being damaged by oxidative stress [22]. McNamara et al. showed that the serum remnant-like particle cholesterol (RLP-c), an independent risk factor for CVD, was higher in the male patients in the high bilirubin group than those in the low bilirubin group [40]. It is showed that bilirubin regulates cholesterol and lipid metabolism via activating the aryl hydrocarbon receptor (AhR) signaling pathway, which has been shown to protect against cardiogenesis and oxidative stress [41, 42]. Notably, bilirubin, as a novel PPARα agonist, suppressed lipid accumulation in 3T3-L1 adipocytes and PPARα KO mice [43]. According to a recent study, increased unconjugated bilirubin promotes the clearance of excess cholesterol out of the body by promoting transintestinal cholesterol secretion [18]. These previous studies have shown that bilirubin can regulate LDL via multiple pathways, which is also in accordance with our findings.

The liver plays a critical role in maintaining cholesterol homeostasis by regulating multiple mechanisms, including exogenous cholesterol uptake via LDLR and cholesterol biosynthesis via HMGCR activity, esterification for storage, and reverse cholesterol transport (RCT) [44]. A disruption in any one of these pathways can lead to dyslipidemia. Newly synthesized cholesterol assembles into very low-density lipoprotein and triglycerides, which is then excreted into the bloodstream or stored in lipid droplets. Ectopic lipid accumulation and the deposition of lipid within non-adipose tissues (e.g., muscle and liver) contribute to the pathogenesis of metabolic diseases, such as diabetes and atherosclerosis [45, 46]. Hepatic steatosis, an indicator of atherosclerosis, can be caused by excessive lipid deposits in the liver. We found that ApoE−/− mice treated with bilirubin displayed a decrease in lipid droplets accumulation in the liver. In addition, bilirubin also inhibited cholesterol and cholesterol esters production, suggesting that bilirubin attenuates hepatic steatosis and atherosclerosis by interrupting cholesterol biosynthesis. Our findings are consistent with a previous study, which demonstrated mild hyperbilirubinemia conditions significantly reduced the intracellular lipid accumulation in C2C12 skeletal mouse muscle and HepG2 human liver cells [47].

HMGCR, a key rate-limiting enzyme, mediates liver de novo biosynthesis of cholesterol. Statins can significantly suppress cholesterol biosynthesis through inhibiting HMGCR activity and thus are widely used in treating hypercholesterolemia. Two negative feedback pathways are available to regulate the expression of target genes involved in cholesterol metabolism: the processing of sterol regulatory element-binding proteins 2 (SREBP2), a crucial transcriptional regulator of cholesterol biosynthesis, and sterol-induced HMGCR degradation [48]. Here, bilirubin did not reduce the mRNA expression of HMGCR, but significantly decreased its protein level, suggesting the involvement of a post-translational mechanism. Song et al. reported that the mTORC1-USP20-HMGCR axis, a post-transcriptional mechanism, participated in the turnover of HMGCR and that the inhibition of USP20 might be beneficial for a variety of metabolic diseases, such as diabetes, hyperlipidemia, and CVD [48]. Our findings demonstrated that bilirubin promotes the ubiquitination and degradation of HMGCR, thus contributing to a reduction of cholesterol biosynthesis and the improvement of atherosclerosis.

Atherosclerosis is considered a chronic inflammatory disease with autoimmune responses. Both the innate and adaptive immune responses have been implicated in atherosclerosis progression. Here, we found that bilirubin displays atheroprotective properties in ApoE−/− mice and acts as a regulatory role on both innate and adaptive immune cells. DCs present antigens to T cells and thus provide an important bridge between innate and adaptive immunity. DCs involves atherosclerosis through multiple mechanisms, inducing lipid uptake, antigen presentation, and cytokine production [49, 50]. There are three DC subsets, including conventional DC (cDC) 1, cDC2, and plasmacytoid DC [51]. cDC1, which express low levels of CD11b, play atheroprotective function through inducing Treg cells [52]. We showed that bilirubin increases spleen CD11blowDCs, which correlates negatively with LDL and TCHO levels. Thus, bilirubin likely exerts atheroprotective effects by regulating the functions of these DCs.

CD4+ T cells are important mediators of the adaptive immune response and their subtypes, Th1 cells, have pro-atherogenic functions mainly through secreting pro-inflammatory cytokines, such as IFN-γ, IL-2, IL-3, and tumor necrosis factor [53, 54]. In the present study, CD4 T cells were significantly reduced by bilirubin treatment and this decrease is positively associated with LDL and TCHO, suggesting the involvement of these CD4 T cells in atheroprotective role of bilirubin. Moreover, after in-depth single cell analysis, a decrease of T6, which mainly are Th1 Tem and accounts for above 15% of spleen T cells, contributes to the overall reduction of CD4 T cells. Considering the pro-atherogenic functions of Th1 cells, decrease of effector memory Th1 cell may be one of the causes of bilirubin-induced atherosclerosis improvement.

NK cells, a critical component of the innate immune system, have been demonstrated as an atherogenic factor. NK cells accumulate in the atherosclerotic plaques, and contribute to the expansion of necrotic cores and atherosclerotic lesion development through producing perforin and granzyme B [55]. Antibodies-mediated depletion of NK cells in ApoE−/− mice greatly attenuated atherosclerosis [55]. Another study has also proved the pro-atherogenic role of NK cells by showing that depletion of NK function decreases atherosclerosis in Ldlr−/− mice [56]. Here, bilirubin-induced decrease of NK cells, especially Ly6C+NK cells, was observed. Due to the pro-atherogenic functions of NK cells, the reduction of NK cells after bilirubin treatment may benefit the remission of atherosclerosis.

MDSCs comprise immature myeloid cells that consist mainly of progenitors and precursors of granulocytes, macrophages, and dendritic cells [57]. MDSCs primarily mediate immunosuppressive function through inhibiting T cell activation and proliferation, promoting the development of regulatory T cells, blocking antigen recognition, and suppressing NK cell cytotoxicity [57, 58]. Based on MDSC phenotypic and morphological features, two major subsets, named polymorphonuclear (PMN) and monocytic (M)-MDSC, were distinguished. In mice, the expression of CD11b, Ly6G, and Ly6C were introduced to identify these two MDSC subpopulations: PMN-MDSC (CD11b+Ly6G+Ly6Clow) and monocytic (M)-MDSC (CD11b+Ly6G−Ly6Chigh) [59]. MDSCs have been extensively investigated in tumor immunology, and accumulating evidence has demonstrated their involvement in obesity and atherosclerosis [60, 61]. By suppressing T cells, MDSCs are able to reduce atherosclerotic lesion development in Ldlr−/− mice [61]. In our study, bilirubin significantly attenuated the atherosclerosis and caused a significant increase in all MDSC, especially PMN-MDSCs. Moreover, MDSCs are also significantly and negatively correlated with LDL and TCHO, two important markers of atherosclerosis. Thus, bilirubin-induced MDSC increase may also improve atherosclerosis.

Overall, our data clearly show the lipid-lowering and cardioprotective properties of bilirubin in atherosclerosis and that the detailed mechanisms involve inhibition of cholesterol synthesis through promoting HMGCR degradation and regulation of the immune systems. Notably, mass cytometry-based single cell analysis was introduced to comprehensively dissect the regulatory effects of bilirubin on different immune cell lineages and T cell subsets. Further, these comprehensive data suggest that bilirubin negatively regulates atherosclerosis not only by affecting cholesterol synthesis but also by reshaping the immune system. In addition to the methodological novelty of our work, our results directly propose bilirubin as a novel therapeutic strategy for improving atherosclerosis treatment.

Availability of data and materials

All datasets generated for this study are included in the article. The original data that support the findings of this study are available from the corresponding author upon reasonable request.

References

  1. Barquera S, Pedroza-Tobias A, Medina C, Hernandez-Barrera L, Bibbins-Domingo K, Lozano R, Moran AE. Global overview of the epidemiology of atherosclerotic cardiovascular disease. Arch Med Res. 2015;46:328–38.

    PubMed  Google Scholar 

  2. Thomas H, Diamond J, Vieco A, Chaudhuri S, Shinnar E, Cromer S, Perel P, Mensah GA, Narula J, Johnson CO, et al. Global atlas of cardiovascular disease 2000–2016: the path to prevention and control. Glob Heart. 2018;13:143–63.

    PubMed  Google Scholar 

  3. Dalager-Pedersen S, Ravn HB, Falk E. Atherosclerosis and acute coronary events. Am J Cardiol. 1998;82:37T-40T.

    CAS  PubMed  Google Scholar 

  4. Agnello L, Bellia C, Scazzone C, Bivona G, Iacolino G, Gambino CM, Muratore M, Lo Sasso B, Ciaccio M. Establishing the 99(th) percentile for high sensitivity cardiac troponin I in healthy blood donors from Southern Italy. Biochem Med (Zagreb). 2019;29: 020901.

    Google Scholar 

  5. Agnello L, Bellia C, Lo Sasso B, Pivetti A, Muratore M, Scazzone C, Bivona G, Lippi G, Ciaccio M. Establishing the upper reference limit of Galectin-3 in healthy blood donors. Biochem Med (Zagreb). 2017;27: 030709.

    Google Scholar 

  6. Vitek L, Ostrow JD. Bilirubin chemistry and metabolism; harmful and protective aspects. Curr Pharm Des. 2009;15:2869–83.

    CAS  PubMed  Google Scholar 

  7. Kang SJ, Lee C, Kruzliak P. Effects of serum bilirubin on atherosclerotic processes. Ann Med. 2014;46:138–47.

    CAS  PubMed  Google Scholar 

  8. Ozturk C, Ozturk A. The relationship between bilirubin levels and atherosclerosis. Angiology. 2015;66:96.

    PubMed  Google Scholar 

  9. Vitek L. Bilirubin and atherosclerotic diseases. Physiol Res. 2017;66:S11–20.

    CAS  PubMed  Google Scholar 

  10. Schwertner HA, Jackson WG, Tolan G. Association of low serum concentration of bilirubin with increased risk of coronary artery disease. Clin Chem. 1994;40:18–23.

    CAS  PubMed  Google Scholar 

  11. Djousse L, Levy D, Cupples LA, Evans JC, D’Agostino RB, Ellison RC. Total serum bilirubin and risk of cardiovascular disease in the Framingham offspring study. Am J Cardiol. 2001;87:1196–200.

    CAS  PubMed  Google Scholar 

  12. Novotny L, Vitek L. Inverse relationship between serum bilirubin and atherosclerosis in men: a meta-analysis of published studies. Exp Biol Med (Maywood). 2003;228:568–71.

    CAS  Google Scholar 

  13. Yang XF, Chen YZ, Su JL, Wang FY, Wang LX. Relationship between serum bilirubin and carotid atherosclerosis in hypertensive patients. Intern Med. 2009;48:1595–9.

    PubMed  Google Scholar 

  14. Kang SJ, Kim D, Park HE, Chung GE, Choi SH, Choi SY, Lee W, Kim JS, Cho SH. Elevated serum bilirubin levels are inversely associated with coronary artery atherosclerosis. Atherosclerosis. 2013;230:242–8.

    CAS  PubMed  Google Scholar 

  15. Kawamoto R, Ninomiya D, Hasegawa Y, Kasai Y, Kusunoki T, Ohtsuka N, Kumagi T, Abe M. Mildly elevated serum total bilirubin levels are negatively associated with carotid atherosclerosis among elderly persons with type 2 diabetes. Clin Exp Hypertens. 2016;38:107–12.

    CAS  PubMed  Google Scholar 

  16. Suh S, Cho YR, Park MK, Kim DK, Cho NH, Lee MK. Relationship between serum bilirubin levels and cardiovascular disease. PLoS ONE. 2018;13: e0193041.

    PubMed  PubMed Central  Google Scholar 

  17. Vitek L, Jirsa M, Brodanova M, Kalab M, Marecek Z, Danzig V, Novotny L, Kotal P. Gilbert syndrome and ischemic heart disease: a protective effect of elevated bilirubin levels. Atherosclerosis. 2002;160:449–56.

    CAS  PubMed  Google Scholar 

  18. Bulmer AC, Verkade HJ, Wagner KH. Bilirubin and beyond: a review of lipid status in Gilbert’s syndrome and its relevance to cardiovascular disease protection. Prog Lipid Res. 2013;52:193–205.

    CAS  PubMed  Google Scholar 

  19. Boon AC, Hawkins CL, Bisht K, Coombes JS, Bakrania B, Wagner KH, Bulmer AC. Reduced circulating oxidized LDL is associated with hypocholesterolemia and enhanced thiol status in Gilbert syndrome. Free Radic Biol Med. 2012;52:2120–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Goff DC Jr, Bertoni AG, Kramer H, Bonds D, Blumenthal RS, Tsai MY, Psaty BM. Dyslipidemia prevalence, treatment, and control in the multi-ethnic study of atherosclerosis (MESA): gender, ethnicity, and coronary artery calcium. Circulation. 2006;113:647–56.

    CAS  PubMed  Google Scholar 

  21. Stocker R, Yamamoto Y, McDonagh AF, Glazer AN, Ames BN. Bilirubin is an antioxidant of possible physiological importance. Science. 1987;235:1043–6.

    CAS  PubMed  Google Scholar 

  22. Ziberna L, Martelanc M, Franko M, Passamonti S. Bilirubin is an endogenous antioxidant in human vascular endothelial cells. Sci Rep. 2016;6:29240.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Vogel ME, Idelman G, Konaniah ES, Zucker SD. Bilirubin prevents atherosclerotic lesion formation in low-density lipoprotein receptor-deficient mice by inhibiting endothelial VCAM-1 and ICAM-1 signaling. J Am Heart Assoc. 2017. https://0-doi-org.brum.beds.ac.uk/10.1161/JAHA.116.004820.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Mazzone GL, Rigato I, Ostrow JD, Bossi F, Bortoluzzi A, Sukowati CH, Tedesco F, Tiribelli C. Bilirubin inhibits the TNFalpha-related induction of three endothelial adhesion molecules. Biochem Biophys Res Commun. 2009;386:338–44.

    CAS  PubMed  Google Scholar 

  25. Keshavan P, Deem TL, Schwemberger SJ, Babcock GF, Cook-Mills JM, Zucker SD. Unconjugated bilirubin inhibits VCAM-1-mediated transendothelial leukocyte migration. J Immunol. 2005;174:3709–18.

    CAS  PubMed  Google Scholar 

  26. Peyton KJ, Shebib AR, Azam MA, Liu XM, Tulis DA, Durante W. Bilirubin inhibits neointima formation and vascular smooth muscle cell proliferation and migration. Front Pharmacol. 2012;3:48.

    PubMed  PubMed Central  Google Scholar 

  27. Ollinger R, Bilban M, Erat A, Froio A, McDaid J, Tyagi S, Csizmadia E, Graca-Souza AV, Liloia A, Soares MP, et al. Bilirubin: a natural inhibitor of vascular smooth muscle cell proliferation. Circulation. 2005;112:1030–9.

    PubMed  Google Scholar 

  28. Hill CA, Fernandez DM, Giannarelli C. Single cell analyses to understand the immune continuum in atherosclerosis. Atherosclerosis. 2021. https://0-doi-org.brum.beds.ac.uk/10.1016/j.atherosclerosis.2021.04.003.

    Article  PubMed  Google Scholar 

  29. Winkels H, Ehinger E, Ghosheh Y, Wolf D, Ley K. Atherosclerosis in the single-cell era. Curr Opin Lipidol. 2018;29:389–96.

    CAS  PubMed  PubMed Central  Google Scholar 

  30. Wang J, Tu C, Zhang H, Huo Y, Menu E, Liu J. Single-cell analysis at the protein level delineates intracellular signaling dynamic during hematopoiesis. BMC Biol. 2021;19:201.

    CAS  PubMed  PubMed Central  Google Scholar 

  31. Wang J, Zheng Y, Tu C, Zhang H, Vanderkerken K, Menu E, Liu J. Identification of the immune checkpoint signature of multiple myeloma using mass cytometry-based single-cell analysis. Clin Transl Immunol. 2020;9: e01132.

    Google Scholar 

  32. Leng E, Xiao Y, Mo Z, Li Y, Zhang Y, Deng X, Zhou M, Zhou C, He Z, He J, et al. Synergistic effect of phytochemicals on cholesterol metabolism and lipid accumulation in HepG2 cells. BMC Complement Altern Med. 2018;18:122.

    PubMed  PubMed Central  Google Scholar 

  33. Khan AA, Agarwal H, Reddy SS, Arige V, Natarajan B, Gupta V, Kalyani A, Barthwal MK, Mahapatra NR. MicroRNA 27a is a key modulator of cholesterol biosynthesis. Mol Cell Biol. 2020. https://0-doi-org.brum.beds.ac.uk/10.1128/MCB.00470-19.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Sever N, Song BL, Yabe D, Goldstein JL, Brown MS, DeBose-Boyd RA. Insig-dependent ubiquitination and degradation of mammalian 3-hydroxy-3-methylglutaryl-CoA reductase stimulated by sterols and geranylgeraniol. J Biol Chem. 2003;278:52479–90.

    CAS  PubMed  Google Scholar 

  35. Wang J, Tu C, Zhang H, Zhang J, Feng Y, Deng Y, Huo Y, Xie M, Yang B, Zhou M, Liu J. Loading of metal isotope-containing intercalators for mass cytometry-based high-throughput quantitation of exosome uptake at the single-cell level. Biomaterials. 2020;255: 120152.

    CAS  PubMed  Google Scholar 

  36. Frohlich J, Dobiasova M, Lear S, Lee KW. The role of risk factors in the development of atherosclerosis. Crit Rev Clin Lab Sci. 2001;38:401–40.

    CAS  PubMed  Google Scholar 

  37. Gurel E, Tigen MK. Protective effect of elevated bilirubin levels on cardiovascular disease in patients with Gilbert syndrome. Turk Kardiyol Dern Ars. 2015;43:591–3.

    PubMed  Google Scholar 

  38. Wang D, Tosevska A, Heiss EH, Ladurner A, Molzer C, Wallner M, Bulmer A, Wagner KH, Dirsch VM, Atanasov AG. Bilirubin decreases macrophage cholesterol efflux and ATP-binding cassette transporter A1 protein expression. J Am Heart Assoc. 2017. https://0-doi-org.brum.beds.ac.uk/10.1161/JAHA.117.005520.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Yan J, Horng T. Lipid metabolism in regulation of macrophage functions. Trends Cell Biol. 2020;30:979–89.

    CAS  PubMed  Google Scholar 

  40. McNamara JR, Shah PK, Nakajima K, Cupples LA, Wilson PW, Ordovas JM, Schaefer EJ. Remnant-like particle (RLP) cholesterol is an independent cardiovascular disease risk factor in women: results from the Framingham Heart Study. Atherosclerosis. 2001;154:229–36.

    CAS  PubMed  Google Scholar 

  41. Dietrich C. Antioxidant functions of the aryl hydrocarbon receptor. Stem Cells Int. 2016;2016:7943495.

    PubMed  PubMed Central  Google Scholar 

  42. Phelan D, Winter GM, Rogers WJ, Lam JC, Denison MS. Activation of the Ah receptor signal transduction pathway by bilirubin and biliverdin. Arch Biochem Biophys. 1998;357:155–63.

    CAS  PubMed  Google Scholar 

  43. Stec DE, John K, Trabbic CJ, Luniwal A, Hankins MW, Baum J, Hinds TD Jr. Bilirubin binding to PPARalpha inhibits lipid accumulation. PLoS ONE. 2016;11: e0153427.

    PubMed  PubMed Central  Google Scholar 

  44. Rogers AJ, Leigh J, Berry G, Ferguson DA, Mulder HB, Ackad M. Relationship between lung asbestos fiber type and concentration and relative risk of mesothelioma. A case-control study. Cancer. 1991;67:1912–20.

    CAS  PubMed  Google Scholar 

  45. Shulman GI. Cellular mechanisms of insulin resistance. J Clin Invest. 2000;106:171–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  46. Goodpaster BH, Wolf D. Skeletal muscle lipid accumulation in obesity, insulin resistance, and type 2 diabetes. Pediatr Diabetes. 2004;5:219–26.

    PubMed  Google Scholar 

  47. Hana CA, Klebermass EM, Balber T, Mitterhauser M, Quint R, Hirtl Y, Klimpke A, Somloi S, Hutz J, Sperr E, et al. Inhibition of lipid accumulation in skeletal muscle and liver cells: a protective mechanism of bilirubin against diabetes mellitus type 2. Front Pharmacol. 2020;11: 636533.

    CAS  PubMed  Google Scholar 

  48. DeBose-Boyd RA. Feedback regulation of cholesterol synthesis: sterol-accelerated ubiquitination and degradation of HMG CoA reductase. Cell Res. 2008;18:609–21.

    CAS  PubMed  Google Scholar 

  49. Jongstra-Bilen J, Haidari M, Zhu SN, Chen M, Guha D, Cybulsky MI. Low-grade chronic inflammation in regions of the normal mouse arterial intima predisposed to atherosclerosis. J Exp Med. 2006;203:2073–83.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. Liu P, Yu YR, Spencer JA, Johnson AE, Vallanat CT, Fong AM, Patterson C, Patel DD. CX3CR1 deficiency impairs dendritic cell accumulation in arterial intima and reduces atherosclerotic burden. Arterioscler Thromb Vasc Biol. 2008;28:243–50.

    CAS  PubMed  Google Scholar 

  51. Guilliams M, Ginhoux F, Jakubzick C, Naik SH, Onai N, Schraml BU, Segura E, Tussiwand R, Yona S. Dendritic cells, monocytes and macrophages: a unified nomenclature based on ontogeny. Nat Rev Immunol. 2014;14:571–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  52. Choi JH, Cheong C, Dandamudi DB, Park CG, Rodriguez A, Mehandru S, Velinzon K, Jung IH, Yoo JY, Oh GT, Steinman RM. Flt3 signaling-dependent dendritic cells protect against atherosclerosis. Immunity. 2011;35:819–31.

    CAS  PubMed  Google Scholar 

  53. Buono C, Come CE, Stavrakis G, Maguire GF, Connelly PW, Lichtman AH. Influence of interferon-gamma on the extent and phenotype of diet-induced atherosclerosis in the LDLR-deficient mouse. Arterioscler Thromb Vasc Biol. 2003;23:454–60.

    CAS  PubMed  Google Scholar 

  54. Gupta S, Pablo AM, Jiang X, Wang N, Tall AR, Schindler C. IFN-gamma potentiates atherosclerosis in ApoE knock-out mice. J Clin Invest. 1997;99:2752–61.

    CAS  PubMed  PubMed Central  Google Scholar 

  55. Selathurai A, Deswaerte V, Kanellakis P, Tipping P, Toh BH, Bobik A, Kyaw T. Natural killer (NK) cells augment atherosclerosis by cytotoxic-dependent mechanisms. Cardiovasc Res. 2014;102:128–37.

    CAS  PubMed  Google Scholar 

  56. Whitman SC, Rateri DL, Szilvassy SJ, Yokoyama W, Daugherty A. Depletion of natural killer cell function decreases atherosclerosis in low-density lipoprotein receptor null mice. Arterioscler Thromb Vasc Biol. 2004;24:1049–54.

    CAS  PubMed  Google Scholar 

  57. Gabrilovich DI, Nagaraj S. Myeloid-derived suppressor cells as regulators of the immune system. Nat Rev Immunol. 2009;9:162–74.

    CAS  PubMed  PubMed Central  Google Scholar 

  58. De Veirman K, Van Valckenborgh E, Lahmar Q, Geeraerts X, De Bruyne E, Menu E, Van Riet I, Vanderkerken K, Van Ginderachter JA. Myeloid-derived suppressor cells as therapeutic target in hematological malignancies. Front Oncol. 2014;4:349.

    PubMed  PubMed Central  Google Scholar 

  59. Bronte V, Brandau S, Chen SH, Colombo MP, Frey AB, Greten TF, Mandruzzato S, Murray PJ, Ochoa A, Ostrand-Rosenberg S, et al. Recommendations for myeloid-derived suppressor cell nomenclature and characterization standards. Nat Commun. 2016;7:12150.

    CAS  PubMed  PubMed Central  Google Scholar 

  60. Xia S, Sha H, Yang L, Ji Y, Ostrand-Rosenberg S, Qi L. Gr-1 + CD11b+ myeloid-derived suppressor cells suppress inflammation and promote insulin sensitivity in obesity. J Biol Chem. 2011;286:23591–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  61. Foks AC, Van Puijvelde GH, Wolbert J, Kroner MJ, Frodermann V, Van Der Heijden T, Van Santbrink PJ, Boon L, Bot I, Kuiper J. CD11b + Gr-1+ myeloid-derived suppressor cells reduce atherosclerotic lesion development in LDLr deficient mice. Cardiovasc Res. 2016;111:252–61.

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

None.

Funding

This work was supported by the National Natural Science Foundation of China (81900365); The National Funds for Developing Local Colleges and Universities (B16056001 to J.L.); Natural Science Foundation research team of Guangdong Province (2018B030312001 to J.L.); and Natural Science Foundation of Guangdong Province (2018A0303130024 to G.W., 2017A030313662 to J.L.).

Author information

Authors and Affiliations

Authors

Contributions

JL designed the experiments. GW, LY, and YH performed the experiments and analyzed data. JW performed mass cytometry and analyzed the data. JL and JW wrote the manuscript. All authors have read and approved the final manuscript.

Corresponding authors

Correspondence to Jinheng Wang or Jinbao Liu.

Ethics declarations

Ethics approval and consent to participate

Animal experiments in our study were approved by the Institutional Committee for the Use and Care of Laboratory Animals of Guangzhou Medical University (2017-014).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Figure S1. (A) ApoE−/− mice fed with a western-type diet were intraperitoneally injected with or without bilirubin. After 10 weeks, mice were scarified, and the body and liver weight were measured. The relative liver weight was expressed as liver weight to body weight ratio. (B) Concentrations of liver enzymes were determined by biochemical analysis. Error bars represent mean ± SD. n (control) =8, n (bilirubin) =14. *P < 0.05. Differences were determined by an unpaired t-test. Figure S2. (A) Quantification of hepatic HMGCR protein expression in ApoE−/− mice treated with bilirubin (n=3) or vehicle (control, n=3). (B) BALB/c mice were received intraperitoneal injection of bilirubin (20 mg/Kg/day, n=4) or vehicle (n=4) once a day for 3 days and sacrificed at day 4. The expression levels of HMGCR in the livers were determined by western blot. Bar plots showing the relative HMGCR protein expression in the livers (right panel). *P < 0.05. Figure S3. After determining the expression of 27 surface markers on the spleen cells using mass cytometry, the T cells expressing CD3 were gated and analyzed using viSNE. viSNE was used to visualize the distribution of the spleen T cells from both groups. Cells on the viSNE map were colored by normalized expression of indicated surface markers. (B) Bar plots showing the frequencies of T1 or T2 in spleen T cells obtained from ApoE−/− mice treated with or without bilirubin. Error bars represent mean ± SD. n (control) =8, n (bilirubin) =14. Figure S4. (A) Dot plots showing the Pearson correlation coefficients for relationships between the concentrations of LDL or TCHO in peripheral blood and frequencies of PMN-MDSCs or All MDSCs. (B) Dot plots showing the Pearson correlation coefficients for relationships between the concentrations of LDL in peripheral blood and the frequencies of Ly6C+NK cells. (C) Dot plots (n=22) showing the Pearson correlation coefficients for relationships between the concentrations of LDL or TCHO in peripheral blood and the frequencies of peripheral blood B cells. Correlations were determined by a Pearson test. Figure S5. After determining the expression of 27 surface markers on the PWB cells using mass cytometry, the T cells expressing CD3 were gated and analyzed using viSNE. (A) viSNE was used to visualize the distribution of the peripheral blood T cells from both groups. Cells on the viSNE map were colored by normalized expression of indicated surface markers. (B) Bar plots showing the frequencies of T1 or T2 in peripheral blood T cells obtained from ApoE−/− mice treated with or without bilirubin. (C) Bar plots showing the frequencies of T28 or T29 in peripheral blood T cells obtained from ApoE−/− mice treated with or without bilirubin (left panel). Dot plots showing the Pearson correlation coefficients for relationships between the concentrations of total bilirubin in peripheral blood and frequencies of T28 or T29 in peripheral blood T cells. (D) Dot plots (n=22) showing the Pearson correlation coefficients for relationships between the concentrations of TCHO in peripheral blood and the frequencies of T7, T9, T13 in peripheral blood T cells. Error bars represent mean ± SD. n (control) =8, n (bilirubin) =14. Differences were determined by an unpaired t-test. Correlations were determined by a Pearson test. Table S1. Mass cytometry antibody reagents.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wen, G., Yao, L., Hao, Y. et al. Bilirubin ameliorates murine atherosclerosis through inhibiting cholesterol synthesis and reshaping the immune system. J Transl Med 20, 1 (2022). https://0-doi-org.brum.beds.ac.uk/10.1186/s12967-021-03207-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12967-021-03207-4

Keywords