Integrating Constituents Absorbed into Blood, Network Pharmacology, and Quantitative Analysis to Reveal the Active Components in Rubus chingii var. suavissimus that Regulate Lipid Metabolism Disorder

Rubus chingii var. suavissimus (S. K. Lee) L. T. Lu (RS)—a sweet plant also known as Tiancha distributed in the south of China where it is used as a beverage—recently gained extensive attention as adjuvant therapy of diabetes and hypertension. Although pharmacological studies indicate that RS has beneficial effects in regulating lipid metabolism disorder characteristics, the active chemicals responsible for this effect remains unclear. The present study aims to predict the effective substances of RS on regulating lipid metabolism disorder through the analysis of the chemical profile of RS, the absorbed prototype components in rat plasma, and network pharmacology. Also, a UPLC method able to quantify the screened potential effective chemicals of RS products was established. First, a total of 69 components—including diterpene, triterpenoids, flavonoids, polyphenols, and lignans—were systematically characterized in RS. Of those, 50 compounds were detected in the plasma of rats administered with RS extract. Through network pharmacology, 9 potential effective components, 71 target genes, and 20 pathways were predicted to be involved in RS-mediated regulation of lipid metabolism disorder. The quantitative analysis suggested that the contents of potential effective components varied among samples from different marketplaces. In conclusion, the presented results provide a chemical basis for further research of Rubus chingii var. suavissimus.


INTRODUCTION
Lipid metabolism disorder is an important pathogenic factor of diseases characterized by an abnormal lipid metabolism, such as atherosclerotic vascular disease, type 2 diabetes, and nonalcoholic fatty liver disease (NAFLD) (Guo et al., 2017;Zhang et al., 2018). While statins are widely used as classical chemical medicines in the treatment of lipid metabolism disorder, they have side effects such as rhabdomyolysis and digestive problems (Jones and Davidson, 2005;Thompson et al., 2016). Recent observations indicating that some ethnomedicines and natural products have lipid-lowering effects with multitarget, multilink, and slight side effect features resulted in the development of combination therapies with western medicines for the treatment of lipid metabolism disorder (Zhang et al., 2018;Li et al., 2020).
Rubus chingii var. suavissimus (S. K. Lee) L. T. Lu (Flora of China, 2003) or Rubus suavissimus S. Lee (RS) (Li, 1981) also known as Tiancha or sweet tea is a sweet plant distributed in the south of China where it is used as a beverage and the adjuvant therapy of diabetes, the treatment of hypertension, and urinary tract infections (Guangxi Food And Drug Administration, 2011). As an ideal sugar substitute and low-calorie, low-toxic healthcare product, it is comparable with Siraitia grosvenorii (Swingle) C. Jeffrey ex Lu et Z. Y. Zhang and Stevia rebaudiana Bertoni (Liang et al., 2003) and is mainly used as a tea drink or sweetener with medicinal characteristics. For example, RS was shown to improve blood lipid regulation due to its richness of a variety of active ingredients (Ezure and Amano, 2011;Koh et al., 2011;Wang et al., 2015). Rubusoside-the sweetness compound of RS accounting for up to 5%-was shown to be able to lower serum total cholesterol (TC), triglycerides (TG) levels, and lipid peroxides (Sun et al., 2001;Tian et al., 2001). Our previous study showed that both RS (Jiang et al., 2021) and rubusoside (Li et al., 2020) could alleviate high-fat-diet-induced lipid metabolism disorder and liver injury in a golden hamster model. Besides, RS-derived polyphenols could reduce abdominal fat, triglycerides (TG), and total cholesterol (TC) in animal models of high-fat diet (Liu et al., 2014;Wang et al., 2015). Moreover, RS-derived flavonoids could regulate lipid metabolism disorder by scavenging free radicals and preventing lipid peroxidation (Liu et al., 2017). However, the systematic chemical basis of RS and the identification of the active compounds active in mediating the lipid metabolism disorder should be further evaluated.
Analysis of compounds migrating into the blood proved to be a very useful approach for the identification of bioavailable components of traditional Chinese medicine (TCM) . Moreover, the combined analysis of compound characterization and blood absorbent compounds is generally accepted for mining the potential biological active constituents of TCM (Bi et al., 2018). With the development of bioinformatics, network pharmacology has become a novel approach to systematically predict and reveal the active components, molecular targets, and action mechanisms of TCM from the molecular level to the pathway level by establishing a multilevel network model of "components−targets−diseases" .
Therefore, this study aims to predict the active components and molecular mechanisms underlying the protective action of RS on regulating lipid metabolism disorder. First, an ultraperfomance liquid chromatography with quadrupole time-offlight tandem mass spectrometry (UPLC-Q/TOF-MS) method was established to characterize the chemical profile of RS as well as the in vivo absorbed components from RS extractadministrated rats. Based on this, a network pharmacology approach was developed to systematically predict the main active components and possible molecular action mechanisms. In addition, the main chemicals predicted to have the potential activity to regulate lipid metabolism disorder were quantified by a validated UPLC method. These results provide a reference for further research and exploration of pharmacodynamics material basis and mechanisms of RS.

Materials and Reagents
Thirteen batches of plant materials from markets were collected (Table 3) and authenticated as the leaves of Rubus chingii var. suavissimus (S. K. Lee) L. T. Lu by Professor Yi Cai, Guangxi University of Chinese Medicine. The voucher specimen is deposited in the herbarium of School of Pharmacy, Guangxi University of Chinese Medicine.

Animals
Male Sprague Dawley (SD) rats (200 ± 20 g) were obtained from Hunan SJA Laboratory Animal Co., Ltd. (Hunan, China). Animal care and procedures were approved by and conducted according to the standards of the Guangxi University of Chinese Medicine (Nanning, China) (Certificate DW20200713-003). The rats were housed in an animal room (22 ± 2°C, 60 ± 10% relative humidity) on a 12 h dark/12 h light cycle with food and water ad libitum. Before experiments, they were fasted for 24 h but with free access to water.

Reference Stock Solutions
Accurately weighted references of caffeic acid, rutin, hyperoside, isoquercitin, quercitrin, quercetin, kaempferol, rubusoside, steviol, isosteviol, ursolic acid, and oleanolic acid were dissolved in methanol, except for gallic acid (in water) and ellagic acid (in DMSO) to get the stock solutions of each reference compound. All reference stock solutions were stored at 4°C.

Chemical Profile of RS by UPLC-Q/TOF-MS
Preparation of Sample Solutions 0.5 g weighed sample powders (RS08) were extracted with 10 mL 60% ethanol at room temperature for 60 min by ultrasonication, then complemented for the weight loss, and centrifuged at 12,000 rpm for 10 min, after which the supernatant fluid was diluted 10 times with 60% ethanol and filtered through a 0.22 μm filter before analysis.

Preparation of Reference Working Solutions
Each reference stock solution was mixed and diluted with methanol to a suitable concentration and then filtered through a 0.22 μm filter before analysis.

Animal Experiment and Plasma Sample Preparation
Preparation of RS Extract for Intragastric Administration to Rats 200.0 g RS powder (RS08) was extracted three times with 1,600 mL 95% ethanol (60 min per time) by ultrasonication, filtered, and recycled the solvent under reduced pressure to obtain the dried extract. It was diluted by 0.5% sodium carboxymethyl cellulose (CMC-Na) to prepare the suspension for animal experiments.

Preparation of Plasma Sample
Ten male SD rats were randomly divided into control group (two) and RS group (eight). Rats in the RS group were gave at a fourfold human equivalent dosage (4 × 2.06 g/kg), and the blood was collected from the orbital veins in heparinized centrifuge tube at 0.17, 0.33, 0.5, 1, 1.5, 2, 4, 6, 8, and 12 h. Rats in the control group received an equal volume of 0.5% CMC-Na orally after which the blood was collected by the abdominal aortic method. The blood samples were centrifuged at 4,000 r/min for 10 min at 4°C to obtain the plasma which was stored at −80°C.
50 μL plasma samples of the above time points were combined (or 500 μL control plasma samples) and loaded onto a preequilibrated solid-phase extraction (SPE) column and successively eluted by 3 ml of water and 3 ml of methanol. The methanol eluent was collected and evaporated to dryness under nitrogen at 40°C. The residue was reconstituted in 100 μL methanol and centrifuged at 12,000 rpm for 10 min to obtain the supernatant fluid for analysis.
The ESI source was operated at an optimized capillary voltage of 3.5 V (ESI -), with a drying gas flow and temperature of 8 L/min at 200°C, a nebulizer pressure of 2.0 bar ( 29 psi), a collision RF of 500 Vpp, a transfer time of 60.0 µs, and a prepulse storage of 5 µs. Automatic MS/MS experiments were performed adjusting the collision energy values as follows: m/z 50, 10 eV; m/z 500, 25 eV; and m/z 1,000, 40 eV. Argon was used as collision gas for CID in the MS E mode. The ion scan ranged from m/z 50 to 1,000. Sodium formate solution was used as a calibration solution. The MS data were processed through DataAnalysis 4.3 software (Bruker Daltonics, Bremen, Germany).

Bioactive Components of RS and Potential Targets
According to the results of the chemical profiling and components absorbed in plasma of RS, the absorbed prototypes in rat plasma were selected as candidate compounds. The molecular targets of RS were identified by PharmMapper (http://www.lilab-ecust.cn/pharmmapper/), STITCH (http://stitch.embl.de/) and SwissTargetPrediction (http:// www.swisstargetprediction.ch/). Beyond this, some known molecular targets of identified compounds were collected for comprehensive analysis through literature mining. All targets were combined and the duplicates were removed.

Screening of Putative Therapeutic Targets
The related targets of lipid metabolism disorder were collected from GeneCards (https://www.genecards.org/). Next, the overlapping genes between RS and lipid metabolism disorder were matched by the Venn diagram.

Network Construction of Interactions Between RS and Lipid Metabolism Disorder
To gain insight into the molecular mechanisms underlying the action of RS on lipid metabolism disorder, a bioactive compound-therapeutic target network was constructed with Cytoscape 3.7.1 to depict the interactions between target genes and bioactive compounds. Next, the degree value was calculated by using the NetworkAnalyzer plug-in. Compounds with a high degree value were deliberated as main bioactive RS-derived compounds potentially involved in the RS-mediated alleviation of lipid metabolism disorder. The protein-protein interaction (PPI) network was supported by STRING (https://stringdb.org/). For the analysis, the study organism was limited to Homo sapiens, while all other settings were kept as default values.

Enrichment of KEGG Pathways
To decipher the possible signaling pathway involved in the intervention of RS on lipid metabolism disorder, the KEGG pathway enrichment analysis was performed with DAVID 6.8 (https://david.ncifcrf.gov/); terms with p ＜ 0.05 were considered significant and captured.

Quantitative Analysis of RS by UPLC
Preparation of Sample Solutions 0.5 g weighed sample powders were extracted with 10 ml 60% ethanol at room temperature for 60 min by ultrasonication, then complemented for the weight loss, and filtered through a 0.22 μm filter before analysis.

Preparation of Reference Working Solution
Gallic acid, caffeic acid, rutin, ellagic acid, rubusoside, and kaempferol stock solutions were diluted into suitable working solutions with corresponding solvents before analysis.

Method Validation
Stock reference solutions of gallic acid, caffeic acid, rutin, ellagic acid, rubusoside, and kaempferol were prepared and diluted to six appropriate concentrations using the corresponding solvent to get the reference working solutions. Six calibration curves were constructed by plotting the content (Y) vs. the peak areas (X) of each analyte. The limit of quantitation (LOD) and limit of detection (LOQ) were determined at a signal-to-noise ratio (S/N) of about 3 and 10, respectively. The intraday precision was examined by six repetitive injections within one day, and the interday precision was carried out for three consecutive days. Repeatability was acquired by analyzing six replicates of the same sample solution (RS08). To determine the stability, the same sample solution (RS08) was injected at 0, 1, 2, 4, 6, 8, 12, and 24 h,  respectively. The accuracy of the method was expressed by the recovery. Approximately 100% amount of each reference standard equivalent to those in 0.25 g extracts were added to the same RS sample (RS08), after which the sample solutions were prepared in six replicates. Recoveries were calculated by the following equation: recovery (%) 100 × (detected amount−original amount)/spiked amount.

Statistical Analysis
All data were expressed as the mean ± standard deviation (SD). T-test was used to compare the two groups. The multigroup comparisons were analyzed by one-way analysis of variance (ANOVA) with SPSS 24.0 (SPSS Inc., Chicago, IL, United States).

Characterization of Chemical Profile of RS by UPLC-Q/TOF-MS
A UPLC-QTOF/MS method was established to characterize the chemical profile. As MS and fragmentations were most informative in the negative ion mode, all analyses were conducted in the negative ion mode. The base peak intensity (BPI) chromatogram of RS in the negative ion mode is shown in Figure 1. The identification data are summarized in Table 1, and the temporarily identified structures are shown in Supplementary Figure S1. The obtained retention time, quasimolecular ions, and the MS n fragments data were compared with reference compounds, a selfbuilt chemical database, and reported literature. In addition, a ClogP value was calculated by ChemDraw 11.0 (CambridgeSoft, United States) to estimate the retention time of the isomers. In total, 69 compounds were identified or tentatively characterized, including 20 diterpenes, 13 triterpenes, 19 flavonoids, 15 polyphenols, and 2 lignans. Compounds 1,9,18,19,20,21,26,32,44,50,64,66,68, and 69 were identified as gallic acid, caffeic acid, ellagic acid, rutin, hyperoside, isoquercitrin, quercitrin, quercetin, kaempferol, rubusoside, steviol, isosteviol, ursolic acid, and oleanolic acid through the comparison of the retention time and quasimolecular ions with reference compounds, respectively. Compounds 4 (corilagin), 19 (rutin), and 50 (rubusoside) were taken as examples for explaining the fragmentation process (Supplementary Figure S2).
Kaurane-type diterpenoids, oleanane-type triterpenoids, ursanetype triterpenoids, flavone, dihydroflavone, anthocyanidins, and flavonol-all widely described in RS-are sometimes present in the form of glycoside with one or several sugar moieties (Hirono et al., 1990;Lv and Wang, 2007;Ohtani et al., 1992;Peng et al., 2019;Wang and Lv, 2008). Upon analysis of the MS/MS spectra of kaurane-type diterpenoids and oleanane, or ursane-type triterpenoids, it was observed that a series of diagnostic ions was caused by the loss of -H 2 O, -OCH 2 , -CO 2, and -C 6 H 9 O 5 (glucose unit). The main characteristic fragmentation behaviors of RS-flavonoids were the loss of glycosyl, elimination of small neutral fragments such as H 2 O and CO, and retro Diels-Alder (RDA) cleavage. RS-derived polyphenols are mainly galloyl-oxygen-diphenyl-type ellagitannins and gallotannins of which the fragmentation patterns could mainly be summarized as successive elimination of the galloyl group and lactonization (Nakano et al., 2017;Schulza and Chim, 2019).

Optimization of the Methods to Collect and Pretreat Plasma
Considering that the time to being absorbed into the blood is different for each RS extract-derived component, plasma samples were first collected at 10 different time points (0.17 ∼ 12 h) after the administration of RS extract and then finally combined. An SPE column was used for protein precipitation and separation.

Identification of Absorbed Components in Rat Plasma Samples after Oral Administration of RS Extract
Based on reference compounds and the chemical profile of RS and by comparing pre-and postdose rat plasma, an in vivo absorbed prototype component profile of the plasma of RS extract-administrated rats was obtained and analyzed.
As shown in Figure 1 and Table 1, a total of 50 compounds-including 19 diterpenoids, 8 triterpenoids, 13 flavonoids, and 10 polyphenols-were identified or tentatively identified in rat plasma samples after the oral administration of RS extract. The detailed identified information, including observed mass values of quasimolecular ions, mass error, and MS n fragments data, is listed in Supplementary Table S1. Most diterpenoids and triterpenoids existed as unchanged prototypes in rat plasma and displayed high peak intensity in the BPI chromatogram. Compared to flavonoid aglycones, more flavone glycosides were detected. The fact that RS-polyphenols displayed low peak intensity or were even hard to detect might be due to their bioavailability and the liver first-pass effect (D'Archivio et al., 2010;Zhang et al., 2016). First, 50 candidate compounds were selected based on the in vivo absorbed components profile of RS. Then, a total of 1,882 molecular targets were obtained from PharmMapper, STITCH, and SwissTargetPrediction. Next, the name of the identified targets was converted into an official gene symbol by UniProt (https://www.uniprot.org/). Finally, a target library with 521 target genes was established after integrating the results and removing the duplicates.

Acquisition of Therapeutic Targets of Lipid Metabolism Disorder
1,112 targets (relevance score>22.01) related to lipid metabolism disorder were gathered from GeneCards. Through Venn diagram matching, 213 overlapping target genes between RS and lipid metabolism disorder-related targets were identified. Of note, a target with a higher relevance score as calculated by GeneCards indicated a higher pertinence to the disease. As such, the top 71 targets ( Table 2) selected by GeneCards were defined as hub genes and retained for subsequent research. Finally, 36 compounds interacting with those previously defined hub genes were identified.

Interaction Network Construction
As shown in Figure 2A, a network linking 36 bioactive compounds and 71 target genes was visualized by Cytoscape 3.7.1. The green nodes represent bioactive compounds and the orange nodes represent putative therapeutic targets. Besides, a protein-protein interaction network ( Figure 2B) illustrates the interconnection of the 71 hub targets. This network suggested the involvement of the 36 compounds and 71 genes in the beneficial effect of RS on lipid metabolism disorder. In network analysis, a high degree value generally indicates the node is closely connected with other target genes. Based on the calculated results from NetworkAnalyzer, compounds with a degree value ≥ 11 were considered as potential bioactive ingredients. As such, nine compounds including ursolic acid (21.0), oleanolic acid (21.0), quercetin (18.0), kaempferol (16.0), gallic acid (14.0), ellagic acid (14.0), caffeic acid (11.0), rubusoside (14.0), and rutin (16.0) might be the main effective RS-derived compounds involved in the beneficial effect of RS on lipid metabolism disorder.

Comprehensive Analysis of the Compounds, Target Genes, and Signaling Pathways Involved in the Beneficial Effect of RS on Lipid Metabolism Disorder
The network pharmacology analysis suggested that nine potential effective components (ursolic acid, oleanolic acid, quercetin, kaempferol, gallic acid, ellagic acid, caffeic acid, rutin, and rubusoside), 71 target genes, and 20 pathways are directly related to lipid metabolism disorder.
Previously, the potential of the above compounds on the regulation of lipid metabolism disorder and antiobesity have Frontiers in Pharmacology | www.frontiersin.org June 2021 | Volume 12 | Article 630198 8 been shown. For example, ursolic acid and oleanolic acid have regulatory effects on the absorption, synthesis, and metabolism of triglycerides and cholesterol (Zhang and Shen, 2015). Also, intake of gallic acid can be beneficial for the suppression of HFDinduced dyslipidaemia, hepatosteatosis, and oxidative stress in rats (Hsu and Yen, 2007). Of note, gallic acid was previously shown to be partially responsible for the antiangiogenic (Liu et al., 2006) and antiobesity (Koh et al., 2011) activities of Rubus leaf. Also, caffeic acid ameliorates hepatic steatosis and reduces endoplasmic reticulum stress in obese mice by regulating autophagy (Kim et al., 2018). Rutin, quercetin, and kaempferol can alleviate high-fat diet-induced obesity and fatty liver insulin resistance (Gao et al., 2013;Bai et al., 2014;Wang et al., 2020). Our previous study indicated that rubusoside effectively alleviated liver injury of golden hamsters on a high-fat diet by reversing the metabolic pathway disorder such as amino acid metabolism and synthesis of ketone bodies (Li et al., 2020).
The interaction network indicated that RS might mediate lipid metabolism disorder through a multicomponents and multitargets approach, suggesting a synergistic effect between these compounds; however, this should be further analyzed in future studies.
The pathway enrichment analysis noted that RS might regulate multiple pathways involved in lipid metabolism disorder, such as insulin resistance and PPAR and AMPK signaling pathways. When lipid metabolism disorder occurs, excess free fatty acids (FFA) are deposited primarily in the liver and β-cells are injured potentially resulting in insulin resistance (IR). Meanwhile, IR could promote the transcription of fatty acid synthetase (FAS) (Meex and Watt, 2017). On the other hand, the PPAR signaling pathway is involved in adipocyte differentiation, fatty acid conversion, and lipid synthesis. PPARc is one of the key transcriptional activators of adipocyte genes such as glucose transporter 4 (Glut4), lipoprotein lipase (LPL), and fatty acid binding protein 2 (aP2). PPARα and AMPK play important roles in regulating the processes of FFA transport and lipid oxidative metabolism (Madrazo and Kelly, 2008;Lopaschuk et al., 2010). Hypoxia-inducible factor-1α (HIF-1α) mainly improves abnormal lipid metabolism by regulating oxidative stress (Ma et al., 2017). Interestingly, LPL, PPARc, and PPARα-genes which are directly related to the PPAR signaling pathway-are among the top 10 hub targets identified in the present study.
Our previous study showed that RS could alleviate lipid metabolism disorder in a high-fat diet golden hamster model and upregulate the expression of the PPAR pathway mediators PPARα and PPARc, resulting in the activation of downstream adipocyte genes aP2, Glut4, and LPL (Jiang et al., 2021). These findings verified that RS might regulate lipid metabolism disorder through mediating the PPAR signaling pathway. Determination of the Content of Compounds with a Potential Effect on Lipid Metabolism Disorder in Different RS Samples by UPLC In the current study, network pharmacology identified nine compounds-including ursolic acid, oleanolic acid, quercetin, kaempferol, gallic acid, ellagic acid, caffeic acid, rubusoside, and rutin-as potential bioactive ingredients of RS involved in mediating lipid metabolism disorder. However, ursolic acid and oleanolic acid could not be detected in the RS samples due to their weak UV absorptions whereas the content of quercetin was below the detection limit. Therefore, six compounds-gallic acid, caffeic acid, rutin, ellagic acid, rubusoside, and kaempferol-were chosen to determine the content in RS samples obtained from multiple marketplaces (RS01-RS13).

Optimization of Extraction and UPLC Chromatographic Conditions
First, several variables including solvent, extraction methods, extraction time, and number of cycles were optimized to obtain the maximal extraction efficiency. Maximal extraction was achieved through ultrasonication with 20 volumes of 60% ethanol for 60 min once or eight volumes of 95% ethanol for three times in case of intragastric administration to rats. Next, the conditions for chromatographic analysis including the type of column, mobile phase gradients, and detection wavelength were optimized. The results showed that optimal separation was obtained by eluting RS samples on an InertSustain C 18 LC column at 40°C using a linear gradient of 0.1% phosphoric acid and acetonitrile within 40 min. Considering that the maximum absorption wavelength of rubusoside situates at 205 nm while that of the other references situates at 254 nm, a wavelength conversion program was used to obtain a better resolution.

Method Validation of UPLC
The regression equation, linearity range, correlation coefficients, and LODs and LOQs of gallic acid, caffeic acid, rutin, ellagic acid, rubusoside, and kaempferol are summarized in Supplementary  Table S3. All calibration curves indicated good linearity (r 2 ≥ 0.9993) within the test ranges while the overall LODs and LOQs were ranged from 0.09 to 1.00 ng and 0.30 to 3.00 ng, respectively. The RSDs of precision, repeatability, stability, and recovery are shown in Supplementary Table S4. For six analytes, the RSDs of intra-and interday precisions of these analytes were less than 1.50 and 1.93%, respectively. The reproducibility RSDs were less than 2.94%. The RS sample solution was stable within 24 h with an RSD of less than 1.89%. The recoveries ranged from 97.58 to 100.29%, and the recovery RSDs ranged from 1.21 to 2.24%. Based on these data, the developed UPLC method was considered reliable, accurate, and suitable for the quantification of those components.

Sample Analysis
The obtained peaks in the chromatograms were identified through the comparison of the retention times and on-line UV spectra with those of the standard constituents ( Figure 4). Retention times of peaks 1-6 were 5. 07, 11.03, 16.93, 17.29, 28.15, and 29.16 min, respectively. The results are summarized in Table 3 and Figure 5. The contents of gallic acid, caffeic acid, rutin, ellagic acid, rubusoside, and kaempferol in the samples ranged from 42.78 to 254.38 μg/g, 31.52 to 195.86 μg/g, 143.60 to 687.29 μg/g, 0.38 to 2.67 mg/g, 36.74 to 68.38 mg/g, and 93.28 to 140.45 μg/g, respectively.
The results showed that the contents of the determined compounds varied in the RS samples obtained from different markets. For example, high levels of ellagic acid and rubusoside were present in all batches. However, even though the samples were purchased from the same origin, the contents could show significant differences. For example, while samples RS06∼ RS08 were all purchased from Laibin of Guangxi, the content of rubusoside and ellagic acid in these samples ranged from 40.15 to 67.92 mg/ g and 0.85 to 2.67 mg/g, respectively. Similarly, the content of rubusoside and ellagic acid in samples RS01∼ RS05-obtained from Guilin-ranged from 36.85 to 68.38 mg/g and 0.38 to 1.15 mg/g, respectively. Importantly, the effectiveness of herbal medicine is highly related to the quality of the raw materials which on itself is affected by a complex set of factors such as the place of origin, growing conditions, medicinal parts, and processing methods. Previously, a higher content of rubusoside in tender tissues of RS (Yin et al., 2008) which was influenced by origin and dosage form (Fan et al., 2012) was shown.
In order to strengthen the quality control of RS, a strict control of the quality of raw materials from the origin and processing procedure is required. Thus, simultaneous determination methods to identify effective bioactive components could be of importance not only to analyze underlying mechanisms but also to establish future quality control of RS.

CONCLUSION
In this study, a total of 69 compounds-including diterpenes, triterpenes, flavonoids, polyphenols, and lignans-were identified or tentatively characterized in Rubus chingii var. suavissimus by UPLC-Q/TOF-MS. Among them, 50 absorbed prototype components were detected in rat plasma. Next, bioactive compounds-therapeutic target network analysis predicted nine active components (ursolic acid, oleanolic acid, quercetin, kaempferol, gallic acid, ellagic acid, caffeic acid, rutin, and rubusoside), 71 target genes, and 20 pathways to be involved in RS-mediated alleviation of lipid metabolism disorder. Finally, UPLC analysis could show significant variability in the contents of six selected potential active compounds in RS samples obtained from different marketplaces. Taken together, the study provided a basis for further exploration of the mechanism on the regulation of lipid metabolism disorder and the development of Rubus chingii var. suavissimus-derived therapeutic drugs and healthcare products.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding authors.