Juice from Fructus Rosae Roxburghii normalizes blood lipids in mice with diet‐induced hyperlipidemia*†

Abstract Fructus Rosae Roxburghii (FRR) as a dietary supplement is considered to possess anti‐atherosclerosis (AS), and hyperlipidemia (HLP) is material basis for AS formation, so the effects and molecular mechanism of FRR on diet‐induced hyperlipidemic mice were explored. In Diet IV2 group, hepatic steatosis was significantly relieved; meanwhile, TC, TG, LDL‐C, HDL‐C, and ASI in serum were regulated to control level. Thirty‐seven DCEG in Diet I, Diet II, and Diet IV2 groups were obtained by RNA‐seq analysis. Relative mRNA levels were further determined by qRT‐PCR, of which 28 genes were matched with those detected by RNA‐seq. Ten DCEP were verified by targeted quantitative proteomic analysis, but expressive patterns of only six proteins were correlated with qRT‐PCR data. These DCEG and DCEP played important roles in regulating the biosynthesis of BAs and steroids, fatty acid metabolism, and LPO production. They might cooperatively regulate the function of HDL or RCT by PPAR signaling pathway under the FRR action. As we know, it is the first time the potential anti‐atherosclerotic mechanism of FRR regulating the blood lipids was explored.


| INTRODUC TI ON
Hyperlipidemia (HLP) is one of the important risk factors leading to atherosclerosis (AS) (Ross & Harker, 1976), so regulating blood lipids and ameliorating lipid metabolism disorders are considered as important strategies to prevent cardiovascular disease (CVD).
HLP induced by lipid metabolism disorders can lead to increase in lipids and lipoproteins in serum (Brischetto, Connor, Connor, & Matarazzo, 1983;Olsson, Pears, McKellar, Mizan, & Raza, 2001), generally characterized by high total cholesterol (TC) or triglyceride (TG), or low high-density lipoprotein (HDL). These lipoproteins play important role in regulating the blood lipid. Hypertriglyceridemia is a prevalent risk factor for CVD. Elevated TG levels are also important markers for several types of atherogenic lipoproteins (Malloy & Kane, 2001). LDL can transport the cholesterol on the surface of cells by LDL receptor and scavenger receptor, which can prevent excessive lipid accumulation (Goldstein & Brown, 2009;Hansson & Libby, 2006). Damaged LDL receptor would prompt the development of AS by accumulating intracellular lipid. LDL can be transformed into oxidized low-density lipoprotein (OX-LDL) under free radical action, which plays a key role for AS formation (Jialal & Devaraj, 1996). HDL can promote the removal of cholesterol from artery walls by reverse cholesterol transportation (RCT). HDL can protect vascular endothelial cells and decrease the intake of smooth muscle cells for LDL, thereby inhibiting the synthesis of glycosaminoglycan (GAGs) and proliferation of smooth muscle cells (Hessler, Robertson, & Chisolm, 1979).
At present, some medicinal and edible plants were considered to regulate blood lipids, such as radix notoginseng powder and ginkgo biloba leaves (Calixto, 2000;Chang et al., 2007). Fructus Rosae Roxburghii (the fruit of Rosa roxiburghii Tratt.f. normalis Rehd. Et Wils., FRR) belonging to Rosaceae family is a special wild resource in Yunnan-Guizhou plateau, which mainly originates from Guizhou Province of China (Li et al., 2016). It has long been considered as healthy fruit improving the digestion, which was recorded for the first time in the Compendium of Materia Medica in 1765 A.D (Liu et al., 2016). Recently, data have shown that FRR has antioxidant, antimutagenic, antiatherogenic, and antitumor effects due to the presence of various phytochemicals, including vitamin C, flavonoids, polyphenol, organic acids, triterpenes, superoxide orgotein dismutase (SOD), and essential oils that have been found in FRR (Burke, Smidt, & Vuong, 2005;Dai & Yang, 2015;Lu & Bao, 2002;Xu, Vidyarthi, Bai, & Pan, 2019). FRR is rich with nutrients and bioactive properties, and can be used as dietary supplement, which can reduce the risk of CVD, especially atherosclerosis (AS) (Zhang et al., 2001).
Some experiments have found that FRR juice can significantly reduce intracellular cholesterol and improve antioxidant capacity of LDL in plasma, thus reducing the incidence of AS (Sotherden, Uto-Kondo, Ayaori, & Ikewaki, 2018).
The anti-atherosclerosis function of FRR is related to its chemical composition, which may be single function or synergistic effect of bioactive substances, but the mechanisms are still unclear. HLP is the material basis for the formation of AS; at present, we analyzed the effects and molecular mechanism of FRR juice on diet-induced hyperlipidemic mice. In recent years, RNA-seq and targeted quantitative protein analysis have become powerful technologies for identification of mRNA and proteins in tissues or cells. The KM mice model with HLP was firstly reproduced; then, expression of key genes and proteins participated in lipid metabolism and related signaling pathways were analyzed in mice, which could explore potential mechanisms of FRR in anti-atherosclerosis. This not only provides theoretical basis for application of FRR in AS, but also suggests a potential beneficial drug for AS.

| The hyperlipidemic mouse model
The healthy mice were divided into four groups with average weight and maintained in different mice cages, with a 12 hr light/dark cycle and free access to purified water. The control mice were fed normal diet (Diet I) starting from 6 weeks of age. The model group was fed high-fat diet (Diet II). The positive drug group (Diet III) and treatment group (Diet IV) were on high-fat diet and were fed separately with Xuezhikang capsules or FRR juice. In Diet IV group, 25% (IV1), 50% (IV2), and 100% (IV3) FRR juice were prepared by mixing purified water, and then were given to mice by gavage (0.25 ml/10 g). On day 30, blood was collected by retro-orbital bleeding; finally, these mice were euthanized by cervical dislocation.

| The liver pathology and blood lipids in mice
After mice were anatomized, their livers were immediately weighed and washed by physiological saline (0.9%). For pathological analysis, some hepatic tissues in mice were fixed in 10% neutral formalin, then processed by rinsing, dehydration, wax immersion, embeddedness, HE staining, and sealing, and finally observed by optical microscope. The rests of hepatic tissue were preserved in liquid nitrogen. The collected blood from mice was allowed to clotting at room temperature and then centrifuged at 10 min (500 g). The upper serum samples were separated and used for analysis of TC, TG, HDL-C, and LDL-C by the BS-800 automatic biochemical analyzer. The hepatic index was calculated by the ratio of liver weight and body weight. The atherosclerotic index (ASI) is usually used for measurement of AS level, which was calculated by the equation as following: ASI = (TC − HDL)/HDL. If ASI is less than 4, it reflects that AS degree is not serious or in reducing. ASI is smaller, AS degree is much lighter, and the risk of CVD is much lower.

| The RNA-seq analysis for hyperlipidemic mice
RNA extraction. The total RNA of hepatic tissues in mice was extracted by TRIzol reagent in accordance with manufacturer's instructions. The genomic DNA was removed from RNA sample using DNase. The concentration, quality, and integrity of RNA were determined using a NanoDrop 1000 spectrophotometer (Thermo Scientific) and an Agilent 2100 Bioanalyzer (Agilent Technologies). RNA samples with 260/280 nm absorbance ratios of approximately 2.0 and RNA integrity numbers >8.0 were processed for RNA-seq.
RNA-seq library construction. The RNA samples from three replicates of nine hepatic tissues of mice fed three experimental diets were prepared 81 libraries. After total RNA was extracted, eukaryotic mRNA was enriched by Oligo (dT) beads, while prokaryotic mRNA was enriched by removing rRNA by Ribo-ZeroTM Magnetic Kit (Epicentre). Then, enriched mRNA was fragmented into short fragments using fragmentation buffer and reverse transcripted into cDNA with random primers. Second-strand cDNA was synthesized by DNA polymerase I, in the presence of RNase H, dNTP, and buffer.
Then, cDNA fragments were purified with QiaQuick PCR extraction kit, end-repaired, poly(A)-added, and ligated to Illumina sequencing adapters. The ligation products were selected by agarose gel electrophoresis, PCR-amplified, and sequenced using Illumina HiSeqTM 4000 (Gene Denovo). After removing adaptor sequences, ambiguous "N" nucleotides (>5% of "N" ratio), and low-quality sequences (<10 of quality score), remaining clean reads were assembled using Trinity software. The clean reads were submitted to NCBI Short Read Archive (SRA) database. For homology annotation, non-redundant sequences were subjected to public databases using BLASTx or BLASTn algorithm (E-value < 10 −5 ), including NCBI non-redundant protein and nucleotide, Gene Ontology (GO), and KEGG.
Bioinformatic Analysis of RNA-Seq Data. Reads obtained from the sequencing machines were filtered by removing reads containing adapters, reads containing more than 10% of unknown nucleotides (N), and low-quality reads containing more than 50% of low-quality (Q-value ≤ 10) bases. The high-quality clean reads were mapped to ribosome RNA (rRNA) to identify residual rRNA reads. The rRNA removed reads were used for further analysis. The gene abundances were calculated and normalized to RPKM (reads per kb per million reads), which were used for comparing the difference of gene expression. Principal component analysis (PCA) performed with R package models (http://www.r-proje ct.org/) was used to reveal the structure/relationship of the samples. To identify differentially expressed genes (DEG) between different treatment groups, the edgeR package (http://www.r-proje ct.org/) was used. The genes with a fold change ≥2 and a false discovery rate (FDR) < 0.05 were identified as DEG, which were then subjected to enrichment analysis of GO functions and KEGG pathways. GO has three ontologies: molecular function, cellular component, and biological process, which were used to describe properties of genes and their products in any organism. All DEG were mapped to GO terms in the GO database (http://www.geneo ntolo gy.org/) and significantly enriched GO terms were defined by hypergeometric test. KEGG as a public pathway-related database could identify significantly enriched metabolic or signal transduction pathways in DEG. The calculating formula is the same as GO analysis.

| The qRT-PCR analysis
In three treatment groups (Diet I, Diet II, and Diet IV2), 37 differentially coexpressed genes (DCEG) were further analyzed to validate the RNA-seq results by qRT-PCR using CFX96 RT-PCR Detection System (Bio-Rad) with SYBR Green Kit (Promega). The extraction and determination of total RNA from the hepatic tissue were performed as in Section 2.4. cDNA was generated from DNase-treated RNA by the FastQuant RT-PCR Kit (with gDNase). The mixture contained 400 ng RNA, 2 μl Fast RT buffer (10×), 1 μl RT Enzyme Mix, 2 μl FQ-RT Primer Mix, and 5 μl RNase-Free ddH 2 O. The reaction conditions were 42°C for 15 min and 95°C for 3 min. The primers were designed according to Illumina sequencing data with Primer Premier 5. The reaction conditions were 95°C for 90 s followed by 40 cycles consisting of 95°C for 5 s, 60°C for 15 s, and 72°C for 20 s.
The fluorescent flux was recorded, and the reaction continued at 72°C for 3 min. The dissolution rate was measured between 65 and 95°C. Each increase of 0.2°C was maintained for 1 s, and fluorescent flux was recorded. A dissociation curve was determined during PCR program to make sure that specific products were obtained in each run. The gene expression was analyzed by 2 −ΔΔCT method and normalized toward the mean of the reference gene (β-actin). The stability of β-actin was checked, which showed that ΔCt values among all treatments were less than 0.5.

| The targeted quantitative proteomics based on mass spectrometry
Protein extraction and trypsin digestion. The sample was grinded in liquid nitrogen into cell powder and then transferred to a 5 ml centrifuge tube. After that, four volumes of lysis buffer (8 M urea, 1% Triton X-100, 10 mM dithiothreitol, and 1% Protease Inhibitor Cocktail) were added to cell powder, followed by sonication three times on ice using a high-intensity ultrasonic processor (Scientz).
The remaining debris was removed by centrifugation at 20,000 g at 4°C for 10 min. Finally, the protein was precipitated with ice-cold 20% TCA for 2 hr at −20°C. After centrifugation at 12,000 g 4°C for 10 min, the supernatant was discarded. The remaining precipitate was washed with cold acetone for three times. The protein was redissolved in 8 M urea, and the protein concentration was determined with BCA kit according to manufacturer's instructions. For digestion, the protein solution was reduced with 5 mM dithiothreitol for 30 min at 56°C and alkylated with 11 mM iodoacetamide for 15 min at room temperature in darkness. The protein sample was then diluted to urea concentration less than 2 M. Finally, trypsin was added at 1:50 trypsin-to-protein mass ratio for the first overnight digestion and 1:100 trypsin-to-protein mass ratio for the second 4 hr digestion.
LC-MS/MS Analysis. The tryptic peptides were dissolved in 0.1% formic acid in water (solvent A), directly loaded onto a homemade reversed-phase analytical column (15 cm length, 75 μm i.d.).
The gradient was comprised of an increase from 6% to 23% solvent B (0.1% formic acid in 98% acetonitrile) over 38 min, 23% to 35% in 14 min and climbing to 80% in 4 min then holding at 80% for the last 4 min, all at a constant flow rate of 400 nl/min on an EASY-nLC 1000 UPLC system. The peptides were subjected to NSI source followed by tandem mass spectrometry (MS/MS) in Q Exactive™ Plus (Thermo) coupled online to the UPLC. The electrospray voltage applied was 2.0 kV. The m/z scan range was from 350 to 1,000 for the full scan, and intact peptides were detected in the Orbitrap at a resolution of 35,000. Peptides were then selected for MS/MS using NCE setting as 27 and the fragments were detected in the Orbitrap at a resolution of 17,500. Cys and oxidation on Met, and max variable modifications were set as 3. Transition settings: precursor charges were set as 2 and 3, ion charges were set as 1 and 2, and ion types were set as b, y, and p. The product ions were set as from ion 3 to last ion, and the ion match tolerance was set as 0.02 Da.

| Statistical analysis
Based on obtained data, the mean and standard deviation in each group were calculated, and one-way analysis of variance was conducted by using SPSS 19.0 software. p < .05 indicates that the two groups are different, and p < .01 indicates that the two groups have significant difference.

| The ameliorating effect of FRR on hepatic pathology
On day 30, the hepatic pathological changes in mice were analyzed by HE staining in different groups. In Diet I group, hepatic lobules were intact, with radially arranged hepatic cords around central vein and no obvious lipid infiltration. In Diet II group, hepatic cords were loosely arranged, and hepatocytes with medium microvesicular steatosis were swollen. In Diet III and Diet IV2 groups, hepatic steatosis was significantly relieved, but no obvious effect in Diet IV1 and Diet IV3 (Figure 1).

| The modulation of FRR on blood lipids
The blood lipids in different treatments were measured on day 30. In Diet II group, TC (2.15 mmol/L), TG (1.12 mmol/L), LDL-C (0.77 mmol/L), hepatic index (4.95 mmol/L), and ASI (1.55 mmol/L) in serum were significantly increased (p < .05) and HDL-C (0.88 mmol/L) was lower (p < .05) compared with Diet I group. In Diet III group, blood lipid level in serum was close to Diet I group, and hepatic index (4.18 mmol/L) was lower than Diet I (p < .05). In Diet IV2 group, TC, TG, LDL-C, HDL-C, and ASI in serum were regulated to control level, and hepatic index (3.83 mmol/L) was significantly lower than Diet I ( Figure 2, p < .05).
A pairwise comparison was performed by using Diet I as the control, and Diet II and Diet IV as the treatments. A total of 1,270 genes were identified as DEG and used for further analysis. The criteria that were used to determine significantly up-or downregulated genes were a twofold or greater change in expression, p-value < .05 and false discovery rate (FDR) < 0.05. As a result, 804 genes (213 upregulated and 591 downregulated) displayed the changes in expression between Diet I and Diet II, 466 genes (184 upregulated and 282 downregulated) between Diet II and Diet IV2 (Figure 3).

| The identification of DCEG
To reveal the molecular mechanism associated with lipid metabolism, 623 genes in the liver for three dietary treatments were evaluated by RNA-seq. Three hundred eighty-eight genes were found in Diet I and Diet II groups; compared with Diet I, 110 were upregulated and 278 were downregulated in Diet II. Two hundred thirty-five genes were found in Diet II and Diet IV2 groups; compared with Diet II, 90 were upregulated and 145 genes were downregulated in Diet IV2. From these 623 genes, only 37 DCEG were found in all three Diet I, Diet II, and Diet IV2 groups (Table 2), 17 genes were upregulated in Diet II and then downregulated in Diet IV2, 16 genes were downregulated in Diet F I G U R E 2 The modulating effects of FRR on blood lipids of KM mice in different treatment groups on day 30 (TC, TG, LDL-C, HDL-C, hepatic index, and ASI) II and then upregulated in Diet IV2, one gene kept the upregulation and two genes downregulation in Diet II and IV2 groups, and one gene was only detected in Diet II and Diet IV2 groups. These DCEG were participated in lipid metabolism directly or indirectly (Table 2).
To characterize the functional consequences of genes in liver, a Gene Ontology (GO) functional enrichment and a KEGG pathway analysis of DEG were performed for each treatment. These genes were mapped to GO databases and classified into three major functional categories such as biological process, cellular component, and molecular function (Figure 4a,b). The GO enrichment analysis showed that a large number of biological processes were significantly altered (qvalue < 0.05) between transcriptomes of three treatment groups (Diet I vs. II and II vs. IV2), respectively. To identify the active biological pathways of liver in different treatments, the annotated coding sequences of differential genes were mapped to KEGG database. Comparing Diet I and Diet II groups, twenty KEGG terms were significantly enriched (qvalue < 0.05, Figure 5a). The twenty KEGG terms between the groups of Diet II and IV2 were also significantly enriched (Figure 5b). The bigger the RichFactor is, the higher the enrichment is. If Q-value is closer to zero, enrichment is more significant. From Figure 6b, it is seen that the pathway of fat digestion and absorption is significantly enriched.

| The qRT-PCR validation of DCEG
To confirm the reliability of the gene expression profiles generated by RNA-seq analysis, qRT-PCR was used to examine the expression levels of 37 DCEG shown in Table 2 in hepatic mRNA samples.
Nine DCEG that had different expression profiles between RNA-seq analysis and qRT-PCR were discarded from further analysis as ambiguous. For 28 of these genes, the expression patterns detected by qRT-PCR were matched with those genes by RNA-seq analysis (Table 3). They might interact with each other to regulate lipid metabolism by specific signaling pathways, which could commonly maintain metabolic homeostasis (Tables 2 and 3).

| Targeted-protein identification
To verify the protein expression of these DCEG, targeted proteins in liver tissues were analyzed by mass spectrometry. In total, 10 DCEP were detected by targeted quantitative proteomic analysis, which were related to lipid metabolism directly or indirectly (Table 4)

| D ISCUSS I ON
In view of rich nutrients and bioactive properties, FRR has been considered as a dietary supplement that could ameliorate the dyslipidemia and finally alleviate the risk of AS (Xu et al., 2019). However, detailed mechanisms are still unclear. It is well known that HLP is primary factor in the pathophysiological process forming AS (Ross & Harker, 1976). At present, we investigate the protection effects and molecular mechanism of FRR juice on diet-induced hyperlipidemic mice and explore the potential anti-atherosclerosis.  (Figure 1). The TC, TG, LDL-C, HDL-C, and ASI in serum were also alleviated to control level; especially, hepatic index was significantly lower than control group (Figure 2 could play an important role in regulating the body metabolism by multifunction or synergistic effect (Qin & Xu, 1998).

Fatty acids metabolism
Transcriptome is defined as complete set of expressed RNA transcript in a cell or tissue, which can provide understanding for cellular activities including the growth, development, and metabolism in organisms, particularly in exploring potential mechanism of nutrients (Ogłuszka et al., 2017;Pertea et al., 2015). RNA sequencing is an effective method that could serve as a platform for future investigations. At present, totally 1,273 genes in liver were generated, of which 623 genes in three dietary treatments (Diet I, Diet II, and Diet IV2) were annotated against KEGG database. To characterize the functions of genes, these DEG in three treatments were analyzed by GO function (Figure 4) and KEGG pathway ( Figure 5). Pathway enrichment analysis indicated that lipid metabolism was significantly changed in mice fed Diet II and Diet IV2, which highlighted the importance of liver with variable rates of lipid metabolism. Further analysis, 37 DCEG were screened out in three treatment groups (Table 2). These genes were involved in lipid metabolism and signal transduction, which showed that biological processes in HLP mice were significantly altered by FRR. To validate the reproducibility and accuracy of RNA-seq analysis, relative mRNA levels of 37 DCEG were determined by qRT-PCR, of which 28 genes were matched with those detected by RNA-seq analysis (Table 3). Their fold changes in genes expression ratios between the RNA-seq and qRT-PCR analyses were positively correlated, which verified the reliability of RNAseq analysis. These results also showed that qRT-PCR expression patterns of genes were consistent with those of RNA-seq. By data analysis, it was found that these DCEG played important roles in regulating lipid metabolism, mainly including biosynthesis of bile acids (BAs) and steroids, fatty acid metabolism, and lipid peroxides (LPOs) production.
The expression of 10 genes was decreased in Diet II group, but returned to control level or obviously increased in Diet IV2 group, Jefcoate & Larsen, 2019). Ppargc1a as a target could also regulate lipid metabolism. It was found that miR-30b-5p might regulate lipid metabolism by targeting the Ppargc1a in hepatocellular carcinoma Huh-7 cells (Zhang et al., 2020). Trib3 and Lepr were both associated with the control of fatty acid synthesis besides for regulating TG, TC, and HDL-C levels in plasma (Angyal & Kiss-Toth, 2012;Domínguez-Reyes et al., 2015). Mat1a plays an important role in amino acid metabolism and formic acid degradation, especially modulating the effect of fatty acids on plasma homocysteine concentrations (Yang et al., 2019;Poursoleiman et al., 2018). Mat1a is also necessary for normal VLDL assembly and lipid homeostasis in plasma (Cano et al., 2011).
The expression of four genes (Tuba4a,Tubb2a,Tubb4b,and Elovl3) was increased in Diet II group, but decreased in Diet IV2 group. The lipid composition in membrane could regulate tubulin interaction with mitochondrial voltage-dependent anion channel (VDAC). The electron transport chain is major cellular source for mitochondrial reactive oxygen species (ROS). αβ-Tubulin heterodimers could inhibit VDAC conductance in lipid bilayers, thus decreasing mitochondrial metabolism (Fang & Maldonado, 2018;Rostovtseva, Gurnev, Chen, & Bezrukov, 2012). It was also found that the agonist pemafibrate could increase the expression of mitochondrial marker Elovl3 involved in fatty acids oxidation (Araki et al., 2018).  (Reue & Brindley, 2008).
The increase in Acot3 and decrease in Chka both appeared in Diet IV2 group. Acot3 could promote the biosynthesis of polyunsaturated fatty acids (PUFAs); ω-3 and ω-6 PUFAs play important role in regulating the blood lipids and preventing the AS (Figure 6). The finding was consistent with previous data. Accumulating evidences also concluded that FRR could attenuate the risk of AS (Goldstein & Brown, 2009;Olsson et al., 2001), because various kinds of bioactive substances in FRR could regulate lipid metabolism by synergistic effect. Chka is the first enzyme in the CDP-choline pathway, which has been implicated in phospholipid metabolism. Its overexpression in several human tumors has been proposed as a novel cancer drug target (Lacal & Campos, 2015). It was seen that lower expression of Chka was beneficial to maintain lipid metabolism under FRR action.
The protein expression corresponding to DCEG was further analyzed by targeted quantitative proteomic analysis. Only 10 proteins were exactly verified, all of them were downregulated in mice fed Diet II and upregulated in Diet IV2 group (Table 4). Among of them, expressive patterns of six proteins (Cyp2C39, Cyp4a10, Cyp7a1, Cyp2C54, Mpo, and Mat1a) were parallel with qRT-PCR.
These data show that mRNA abundance of specific gene was not necessarily a linear relationship with protein expression, because gene expression was regulated at different levels, especially during transcription and translation, which would affect the protein expression (Gautier et al., 2012;Spriggs, Bushell, & Willis, 2010). The degradation of mRNA and proteins, and protein modification may also lead to this inconformity. Furthermore, the abundance of some proteins might be too low to be detected by targeted quantitative proteomics.
In a word, FRR juice could ameliorate the blood lipids and maintain metabolic balance by regulating the expression of these DCEG and DCEP in mice. Despite different expression patterns at molecular level, they jointly regulated lipid metabolism and normalize the blood lipids in mice. In our experiments, HDL level in mice fed Diet IV2 was increased significantly, and expression of genes and proteins involved in PPAR signaling pathway was also altered. It can be seen that FRR juice could regulate the function of HDL or RCT by PPAR signaling pathway, thus maintaining normal level of blood lipids in mice, which might be important drug targets for AS. This would be very beneficial to explore anti-atherosclerotic mechanism of FRR and provide theoretical basis for application of FRR in AS.

| CON CLUS IONS
Under the FRR action, the TC, TG, LDL, HDL, and ASI in serum were regulated to control level. Thirty-seven DCEG in each group were obtained by RNA-seq analysis. Relative mRNA levels of 28 genes determined by qRT-PCR were matched with those detected by RNA-seq analysis. Ten DCEP were verified by targeted quantitative proteomic analysis, but expressive patterns of only six proteins were correlated with qRT-PCR data. These DCEG and DCEP played important roles in regulating lipid metabolism, mainly including biosynthesis of BAs and steroids, fatty acid metabolism, and LPO production. They might cooperatively regulate the function of HDL or RCT by PPAR signaling pathway under the action of FRR juice.

ACK N OWLED G M ENTS
The authors would like to thank Dr Alexander Kuryatov (University

E TH I C A L A PPROVA L
The study's protocols and procedures were ethically reviewed and approved by the Institutional Animal Care and Use Committee (IACUC) of Guizhou Medical University of China. The authors agree with the paper publication. The data and materials in the study are shared and available. The authors declare that they have no any conflict of interests.