Biochemical and Molecular Investigation of the Effect of Saponins and Terpenoids Derived from Leaves of Ilex aquifolium on Lipid Metabolism of Obese Zucker Rats

Ilex paraguariensis, the holly tree, is a plant with recognized biological properties, whose aqueous infusions are known as “Yerba mate”, that regulate lipid metabolism, reduce obesity, and improve brain stimulation. In the present study, the effect of standardized saponin and terpenoid fractions of a European taxon, Ilex aquifolium, on blood biochemical parameters in a rat model of metabolic disorder, (fa/fa) Zucker, are presented. The profiles of the volatile fractions of two species and six European varieties of Ilex were investigated. After selecting the best variety, the saponin and terpenoid fractions were isolated and standardized, and animals were fed 10 mg kg−1 b.w. for 8 weeks. A statistically significant decrease in liver adiposity was observed, confirmed by histology and quantitative identification (gas chromatography–mass spectrometry analyses of hepatic lipids. RT-qPCR analysis of gene expression in the aorta revealed that the administration of the terpenoid fraction downregulated LOX-1, suggesting a reduction in atherosclerotic stimuli. In addition, a statistically significant reduction (p < 0.05) in PPARγ for the saponin fraction was observed in the liver. The expression of the ACAT-1 gene in the liver, responsible for the formation of cholesterol esters, increased significantly in the group receiving the terpenoid fraction compared to the control, which was also confirmed by the analysis of individual blood biochemical parameters. The opposite effect was observed for saponins. Taking the above into account, it is shown for the first time that Ilex aquifolium can be a source of compounds that positively influence lipid metabolism.


Introduction
Metabolic disorders such as hepatic steatosis represent a growing problem in modern populations. In addition to classical pharmacotherapy, so-called herbal medicines or dietary supplements that regulate the lipid metabolism of the organism may constitute an alternative. Currently, the most widely used herbal drugs are extracts containing flavonolignan silymarin from Silybum marianum and the derivative hydroxycinnamic acid cynarin of Cynara scolymus, which are used in the treatment of liver damage; they protect the liver from toxic compounds and improve its regeneration. These compounds also exert spasmolytic and cholagogic effects [1]. According to Speroni et al., artichoke leaf extract shows protective effects against oxidative stress in hepatocytes [2]. Lipid-regulating herbs also include those rich in saponin fractions, e.g., Momordica charantia, Gymnema sylvestre, Gypsophila oldhamiana, or Achyranthes aspera. Their mechanisms of action range from appetite regulation, and the inhibition of pancreatic lipase to effects on triglyceride absorption, retrosynthesis, and adipogenesis [3].
In recent years, interest in the traditional South American beverage of "Yerba mate", the water macerate of Ilex paraguariensis, has grown tremendously, which has resulted in a significant increase in research and literature on its biological properties. The active components responsible for its biological activities are: polyphenols, saponins, triterpenoids, and finally caffeine. These have exhibited effects such as vasodilation and lipid lowering, anti-glycation, and anti-obesity effects [4]. In addition, I. paraquariensis has mutagenic or antimutagenic properties, depending on the model. It affects the stimulation of the nervous system and has diuretic properties. "Yerba mate" was also found to protect DNA from oxidation and low-density lipoprotein lipoperoxidation in vitro [5]. The methanolic extracts of the leaves present high antioxidant activity and phenolic compound content such as chlorogenic acid derivatives, flavonoids, and saponins [6]. In 2020, Zapata [7] published results (obtained on a Wistar model with a high-fat and high-sugar diet) indicating that only caffeine is responsible for suppressing lipid accumulation, proven by its modulation of gene expression. Caffeine-free, supercritical CO 2 extracts of I. paraguariensis did not affect fatty acid synthase (Fasn), pyruvate kinase or the microsomal triglyceride transfer protein.
The European varieties of Ilex aquifolium and Ilex meserveae have not yet been as thoroughly characterized as the South American I. paraguariensis for their biochemicals and activities. Our previous experience indicates that, despite geographical separation, the chemical composition of the European and South American Ilex varieties is similar. The I. aquifolium and I. meserveae cultivars are, as with I. paraguariensis, rich in polyphenols, saponins, and triterpenoid fractions (with derivatives of ursolic acid, oleanolic acid, etc.). The main difference in chemotype is the complete absence of the xanthine fraction, i.e., caffeine or theobromine. In the context of the publication by Zapata et al. [7], where it was suggested that for the paraguariensis, the only active fraction (for regulation of lipid metabolism) is xanthines, it seems reasonable to test decaffeinated European varieties. Our preliminary study confirmed the renal protective effect for a high-cholesterol diet of aqueous extracts and saponin and terpenoid fractions of the I. aquifolium variety on a Wistar model [8]. Therefore, it seems reasonable to test whether decaffeinated I. aquifolium extracts, as well as saponin and terpenoid fractions, will have a protective effect on the metabolism of the Zucker knockout rat. The purpose of this study was to determine whether the fractions of I. aquifolium would alter the histopathological and biochemical picture of the Zucker rat, as well as to determine changes in the expression of selected genes responsible for normal body metabolism.

Aroma Profile of Ilex
This is the first paper to present a comparison of the volatile compound profiles in Argentine I. paraguariensis and European I. aquifolium and I. meserveae. The expected results of the study of the volatile fractions of European holly species were the appearance of rich monoterpene and monoterpenoid profiles, which have not yet undergone oxidation reactions during heat treatment. The analysis of volatile compounds in different varieties of Ilex revealed a total of 81 compounds, of which the ten main ones are presented in Table 1 (in addition, Table S1 showing all compounds is available in the Supplementary Data ). The results obtained show a definite difference in the volatile compound profile between Argentine I. paraguariensis and European cultivars of I. aquifolium and I. meserveae. This difference was due to the fact that I. paraguariensis was subjected to a prior drying process that contributed to the formation of more oxidation intermediates, characterized by a lower boiling point. The profiles presented by the different groups can represent the characteristic distribution of the chemotype within a given cultivar and act as specific marker indicators to assist in the identification of individual holly species. A divergence was also observed among the profiles of European cultivars. In I. paraguariensis, the main volatile compounds were found to be 6-methyl-5-heptene-2-one (25.73 µg g −1 ), (E,E)-3,5octadien-2-one (22.82 µg g −1 ), (E,E)-2,4-heptadienal (14.83 µg g −1 ), hexanal (13.80 µg g −1 ) and 2-methyl-2-pentenal (11.77 µg g −1 ). The two main compounds found in I. paraguariensis (6-methyl-5-heptene-2-one and (E,E)-3,5-octadien-2-one) were absent in the European varieties. In I. aquifolium, (E,E)-2,4-heptadienal was also absent. and p-cymene (58.21-82.62 µg g −1 ), α-phellandrene (3.68-5.35 µg g −1 ) and α-pinene (1.98-3.94 µg g −1 ) were the main compounds. In I. meserveae, p-cymene (20.20-32.75 µg g −1 ) was also among the main compounds, but the other compounds were (E)-2-hexenal (12.28-20.97 µg g −1 ) and 3-hexen-1-ol (6.34-11.85 µg g −1 ). Regarding the profile of volatile compounds in I. paraguariensis, this has already been reported in several works. A study by Dallago et al. [9] compared the profile of volatile compounds contained in I. paraguariensis both before the drying process (fresh green leaves) and after the drying stage. The results presented differences in the number of Maillard reaction products formed and their percentage ratios. It was observed that the effect of a higher temperature and oxygen access favored the formation of more secondary products such as ketones, aldehydes, and carboxylic acids. Fresh leaf samples were found to be dominated by compounds from the alcohol group, while dry leaf samples were dominated by aldehydes and ketones. It was also found that long hydrocarbon chains are degraded first, followed by the oxidation of individual functional groups. One paper that shows values similar to those of our team in the profile of volatile compounds of I. paraguariensis is the publication by Araujo et al. [10]. They showed that the main aromatic compounds in mate are (E,E)-2,4-heptadienal (7.8%), (E,Z)-2,4-heptadienal (7.0%) and (E,Z)-3,5-octadien-2-one (5.2%). However, they obtained a much lower value for 6-methyl-5-hepten-2-one (1.4%), which in our study, proved to be the main component of the volatile profile. In a study by Marquez et al. [11], they observed that geranyl acetate (9.38 µg g −1 ), linalool (4.78 µg g −1 ) and isomers (E,E)-and (E,Z)-2,4-heptadienal (3.60 and 1.21 µg g −1 ) were the most abundant in Yerba mate leaves. Conversely, Martins et al. [12] conducted a study on the volatile compounds released in relation to mechanical damage and herbivore exposure. In the control group, which was I. paraguariensis, the main volatile compounds reported were decanal, nonanal, and limonene.

Terpenoid Profile
The analysis of terpenoids from two cultivars of European Ilex, i.e., I. aquifolium and I. meserveae, showed the presence of 12 compounds, of which the seven main ones are presented in Table 2 (in addition, Table S2 showing all compounds is available in the Supplementary Data). The study showed similarities in the terpenoid profiles of the European cultivars. The dominant compounds in both I. aquifolium and I. meserveae were ursolic acid (7.29-16.11 mg g −1 ), oleanolic acid (2.03-6.62 mg g −1 ), α-amyrin (0.36-3.22 mg g −1 ), lupeol (0.42-2.53 mg g −1 ) and uvaol (0.37-2.48 mg g −1 ). For comparison, the reference sample was I. paraguariensis, whose terpenoid profile was similar to those of the European varieties, with the difference being that a significantly lower content of ursolic acid (1.23 mg g −1 ) was observed in Yerba mate leaves. Previous studies have also reported the presence of tritepenes in other Ilex cultivars. Lupeol, betulonic acid, uvaol, ursolic acid, and α-amyrin have been identified in the leaves of I. cornuta and I. latifolia [13]. Compounds from the terpenoid group were also documented in I. centrochinensis and I. macropoda. In the first species, the presence of lupeol and oleanolic acid was revealed, and in the second, betulin was isolated in addition to lupeol [14]. Plants with a high content of pentacyclic tritepenes are often used in phytotherapy due to their valuable medicinal properties. Additionally, they are widely distributed in the plant world and are the subject of phytochemical and pharmacological research. One of them is Ficus carica, a traditional plant in folk medicine used to combat pneumonia, diarrhea, inflammation, and indigestion. Ivanov et al. [15] began to identify the constituents of the nonpolar fraction using GC-MS, and the results of their study show that pentacyclic triterpenoids α-Amyrin, β-amyrin and lupeol are present in ficus leaves. Wolbiś et al. [16], in their study of the quats and leaves of Prunus spinosa, which is known as a tea substitute, proved that more triterpene acids are present, and these are mainly ursolic acid (9.5 mg g −1 ) and oleanolic acid (2.3 mg g −1 ). Conversely, Kowalski [17] obtained similar values for ursolic acid (14.98 mg g −1 ) and oleanolic acid (approximately 5 mg g −1 ) [17] in the leaves of Silphium integrifolium, whose infusion was used to treat urinary tract disorders. Cheun and Zhang [18], in Prunella vulgaris spikes, identified mainly betulinic acid, ursolic acid, and oleanolic acid, whose amounts were 4.7, 4.0 and 0.9 mg g −1 , respectively. A terpenoid fraction was also isolated from the herb Hieracium pilosella [19], used, in the treatment of skin diseases, among others, due to its astringent, antiseptic, and anti-inflammatory activities. This fraction consists of β-amyrin, lupeol, and α-amyrin. In conclusion, terpenoids are a widely distributed group among medicinal plants that show promise for the development of new, multifunctional bioactive agents.

Saponin Profile
Following our previous UHPLC-MS/MS analyses, the saponin profiles of I. aquifolium and I. paraguariensis were compared under the conditions described above [8]. By comparing the saponin profiles of holly and commercial mate, we found that the saponin pattern in I. aquifolium is less complicated than in I. paraguariensis, as shown in Table S3 (available in Supplementary Materials). Therefore, we have tentatively designated 29 and 53 compound respectively, as triterpene glycosides or glycoside esters.
Qualitatively, the profiles were quite different: 16 compounds were found in both species, while 50 were unique to one or the second. In I. paraguariensis, one compound (911, 12.54 min; only in mate), tentatively assigned as matesaponin 1, was the dominant one, while another compound (911, 13.98 min; in both species) appeared to be the main saponin in I. aquifolium. However, numerous already-recognized saponins in the Ilex genus match the MS fragmentation of the last saponin, for example, ilexsaponin B2, latifoloside A, latifoloside B, or latifoloside D; thus, we cannot speculate on dereplication.
Considering the compounds classified as the same in both holly and mate, only six of sixteen had quite small differences below 10 RA% between species (911, 6

Fatty Acid Profile in Ilex Leaves
Analysis of fatty acid methyl esters (FAMEs) from European Ilex cultivars revealed the presence of 17 compounds, and seven of the main FAMEs are summarized in Table 3 (in addition, Table S4 showing all compounds is available in Supplementary Data). It was found that the fatty acid profiles of the Argentine and European varieties were similar to one another. The main compounds identified in the European varieties were palmitic acid (12.05-15.03 mg g −1 ), linoleic acid (LA) (5.52-10.40 mg g −1 ) and α-linolenic acid (ALA) (9.09-22.16 mg g −1 ). Regarding the European cultivars, greater profile similarity was observed between I. aquifolium Ferox Argentea and I. meserveae Blue Angel due to similar amounts of LA (5.84 and 5.63 mg g −1 ) and ALA (9.72 and 9.09 mg g −1 ) acids. However, the other cultivars, i.e., I. aquifolium Alaska, I. aquifolium Rubricaulis Aurea, I. meserveae Golden Girl, and I. meserveae Blue Boy, showed higher amounts of these acids, which were as follows: LA (8.40-10.08 mg g −1 ) and ALA (19.28-22.16 mg g −1 ). I. paraguriensis as a reference sample showed the highest amounts of palmitic (18.08 mg g −1 ) and oleic acid (9.03 mg g −1 ). Stearic and α-linolenic acids were identified at similar levels, namely 6.39 and 5.92 mg g −1 respectively. It was also observed that the level of linoleic acid was lower in the Argentinean cultivar compared to the European cultivars. The fatty acid profile in I. paraguriensis has also been reported by other researchers. In a paper published by Souza et al. [6], the acid that was significantly prominent in terms of the amount in the leaves of I. paraguariensis was α-linolenic acid, and the second acid that was most prominent, but in a much lower amount, was palmitic acid. Compared to the above analyses, the levels of palmitic and LA acids were lower. In addition, in the work of Reis et al. [20] the fatty acid profile of I. paraguariensis was demonstrated, where the main acids were also palmitic acid and α-linolenic acid, but their amounts were similar. As in the previous study, Reis et al. also reported lower contents of stearic acid, oleic acid, and linoleic acid. Chóez-Guaranda et al. [21], in their study on the oxidative activity of fractions of I. guayusa leaves, proved the presence of the following fatty acids in the hexane fraction: palmitic acid, oleic acid, and stearic acid.
Similar fatty acid relationships have been identified in medicinal plants in previous studies. Guil-Guerrero [22], in a study on Plantago major, determined the highest amount of α-linolenic acid (40.04%), followed by palmitic acid (16.59%) and linoleic acid (13.77%). The relationships among fatty acid profiles in medicinal plants were also recorded by Rutto [23], who explored the dietary value of raw and processed Urtica dioica. Stinging nettle showed the highest α-linolenic acid content (49.55%). The opposite situation to the previous case was also observed here, since the second highest acid content was attributed to linoleic acid (23.30%), followed by palmitic acid (17.06%). However, in the fatty acid profile of Melissa officinalis leaves, linoleic acid was found to be the dominant acid (74.08%), while palmitic acid was next (15.77%) [24]. In the leaves of Salvia officinalis, Taarit [25] identified those with 18 carbons in the molecule as the main fatty acids. The acid content decreases with the decrease in the number of double bonds in the acid (or with the decrease in the degree of unsaturation) and is as follows: α-linolenic acid (45.80%), linoleic acid (14.00%) and oleic acid (9.41%).
The polyunsaturated fatty acid/saturated fatty acid (PUFA/SFA) ratio is used when estimating the effect of diet on the cardiovascular system. PUFAs are assumed to lower serum low-density cholesterol (LDL), as opposed to SFAs, which contribute to high serum cholesterol [26]. It follows that the higher the PUFA/SFA ratio, the more positive the effect on the cardiovascular system. The most beneficial plants were I. meserveae Golden Girl (1.90), I. aquifolium Alaska (1.96) and I. aquifolium Rubricaulis Aurea (1.70), showing the highest indices; thus, compared to I. paraguariensis (0.57), they can be considered more beneficial to the circulatory system. However, the PUFA/SFA ratio is too general and does not allow for an assessment of atherogenicity. Therefore, an atherogenicity index (IA) is used, which reveals the ratio of the sum of SFAs and the sum of unsaturated acids (UFAs). SFAs (C12:0, C14:0, C16:0) promote lipid adhesion to the cells of the circulatory and immune systems. However, UFAs inhibit the accumulation of atherosclerotic plaques and influence the reduction of phospholipids, cholesterol, and esterified fatty acids, which is why they are considered anti-atherosclerotic [27]. According to this rule, products with lower IA may  This suggests that European cultivars may have a more beneficial effect than I. paraguariensis in lowering plasma cholesterol levels. An indicator that shows the tendency to form thrombi in blood vessels is the thrombogenicity index (IT). It characterizes the relationship between FAs that are considered prothrombogenic (C12:0, C14:0 and C16:0) and antithrombogenic (n-3, n-6 and MUFA). Therefore, a diet containing products with a lower IT ratio is more beneficial for the cardiovascular system [28]. is similar to that of I. paraguariensis (0.69). In summary, both of these indices (IA and IT) can be used as a means to assess the effect of the composition of Fas on the cardiovascular system. Products with lower IA and IT compositions are associated with a reduced risk of coronary heart disease and better nutritional value.
Another indicator that can better illustrate the effect of the composition of Fas on cardiovascular disease than PUFA/SFA is the hypocholesterolemic/hypercholesterolemic index (HH). It describes the relationship between hypocholesterolemic fatty acids (cis-C18:1 and PUFA) and hypercholesterolemic Fas (C12:0, C14:0 and C16:0). As the value of the h/H index increases, the beneficial effect on the circulatory system also increases [29]. Here, the highest values were also found in I. aquifolium (2.68), I. meserveae Golden Girl (2.63), I. aquifolium Rubricaulis Aurea (2.41) and I. meserveae Blue Boy (2.39). They were twice as high as those of the reference samples, I. paraguariensis (1.24) and I. aquifolium Ferox Argentea (1.17).
Unsaturated fatty acids show different weights in relation to the unsaturation index (UI). This reflects the complex proportions of FA, with different degrees of unsaturation compared to the total FA composition. This index can be used to evaluate the suitability of the products as alternative sources of high-quality PUFAs [30,31]. Regarding the Argentine I. paraguariensis (110.07), all the European varieties tested obtained higher UI indices. The variety most similar to the reference sample was I. meserveae Blue Angel (126.82). The highest UI values were obtained for I. aquifolium (174.73) and I. meserveae Golden Girl (175.87).
The results obtained indicate that European Ilex varieties show potentially positive effects on human health through lowering cholesterol levels, reducing thrombus formation, and reducing the risk of cardiovascular disease.

Feed Consumption, Body Weight and Liver Weight
The initial body weight, the final body weight, and the food intake (Table 4) did not change between the groups. However, a tendency toward decreased body weight after the application of the terpenoid fraction was observed only at the final stage of the experiment. The dynamics of the changes in body weight during the supplementation period are presented in Figure 1. Feed intake in rats with metabolic syndrome affected liver weight, which was highest in the TERP group and lowest in SAP. These differences were not statistically confirmed. In another study, a hypercholesterolemic diet was shown to affect body weight [32]. Conversely, changes in adipose tissue stores may affect the expression of genes responsible for the synthesis of lipolytic hormones and trigger mechanisms related to cholesterol and free fatty acid levels. This has also been confirmed in our own studies. Obese Zucker rats had significantly higher body and organ weights and plasma levels of TNF-α, insulin, and leptin than lean animals [33]. Our own research suggests that these changes may be partially ameliorated by the use of saponins, especially with regard to body weight. The results obtained indicate that European Ilex varieties show potentially positive effects on human health through lowering cholesterol levels, reducing thrombus formation, and reducing the risk of cardiovascular disease.

Feed Consumption, Body Weight and Liver Weight
The initial body weight, the final body weight, and the food intake (Table 4) did not change between the groups. However, a tendency toward decreased body weight after the application of the terpenoid fraction was observed only at the final stage of the experiment. The dynamics of the changes in body weight during the supplementation period are presented in Figure 1. Feed intake in rats with metabolic syndrome affected liver weight, which was highest in the TERP group and lowest in SAP. These differences were not statistically confirmed. In another study, a hypercholesterolemic diet was shown to affect body weight [32]. Conversely, changes in adipose tissue stores may affect the expression of genes responsible for the synthesis of lipolytic hormones and trigger mechanisms related to cholesterol and free fatty acid levels. This has also been confirmed in our own studies. Obese Zucker rats had significantly higher body and organ weights and plasma levels of TNF-α, insulin, and leptin than lean animals [33]. Our own research suggests that these changes may be partially ameliorated by the use of saponins, especially with regard to body weight.

Hematological and Biochemical Parameters
The experimental group supplemented with terpenoids showed a significant increase (p < 0.05) in WBCs (Table 5) in the blood after 8 weeks compared to the control group. Similar changes were also found in the RBC (p < 0.01). In the group that received terpenoids, an increase in WBCs and RBCs was also observed, but no statistical differences  4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Weight

Hematological and Biochemical Parameters
The experimental group supplemented with terpenoids showed a significant increase (p < 0.05) in WBCs (Table 5) in the blood after 8 weeks compared to the control group. Similar changes were also found in the RBC (p < 0.01). In the group that received terpenoids, an increase in WBCs and RBCs was also observed, but no statistical differences were observed. The supplementation used also resulted in HGB and HCT (Table 5). No clear differences were found between the different red cell indices. Conversely, platelets reached the highest concentration in the TERP group, and were slightly lower for SAP, while PLT levels were lowest in the control group. Other studies indicate significantly higher plasma levels of RBC, WBC, monocytes and platelets, and leptin in obese rats [33]. The greatest changes in lipid parameters were found after terpenoid treatment (Table 6). There was a significant increase (p < 0.057) in TG concentration in the TERP group. Increases in the concentrations of total cholesterol, LDL and HDL fractions were also evident in the TERP group compared to the control group. However, these changes were not statistically significant. The use of saponins caused a slight decrease in the concentrations of cholesterol and LDL, and HDL. Supplementation with saponins caused a marked increase in TG concentrations compared to the control group. There was an increase in lipolysis in the TERP group, and the statistical differences compared to the SAP group were clear. The serum NEFA concentration was lowest in the SAP group. The administration of terpenoids or saponins to Zucker rats had no effect on AST activity. Genetically obese Zucker (fa/fa) rats overexpress leptin in the stomach. Gastric leptin expression is regulated by feeding conditions in lean but not obese rats. This would indicate the impaired regulation of leptin expression in the stomachs of obese Zucker rats (as in adipose tissue), possibly related to the lack of the leptin receptor [34]. In a study in mice with an obesity-inducing diet, Yerba extract supplementation was shown to alleviate hyperglycemia and improve insulin sensitivity and plasma lipids [35]. Zucker rats receiving saponins and terpenoids had reduced serum insulin and glucose levels. This indicates their potential to reduce insulin resistance in obese Zucker rats. Our previous study showed that I. aquifolium and its terpenoids had an impact on serum lipid levels in an animal model of hyperlipidemia [32]. I. aquifolium improved insulin sensitivity and decreased blood concentrations of total cholesterol and LDL. Supplementation with terpenoids may be an effective means to reduce oxidative damage to lipids, as it decreases the concentration of MDA and leads to increased blood TAS. However, terpenoids also appeared to mitigate liver damage [32]. In general, a decrease in antioxidant capacity leads to endothelial dysfunction and is characterized by a decrease in NO bioavailability of vasodilator nitric oxide and an increase in endothelium-derived contractility factors, which lead to atherosclerosis [36]. In our study, age-related increases in terpenoids caused an increase in hyperlipidemia and total cholesterol levels, with a trend toward an increase in the LDL and HDL fractions. However, the most intense changes were observed in the TERP group.

Oxidative Status
The mean values of the antioxidant status parameters are summarized in Table 7. Long-term administration of saponins or terpenoids to Zucker rats resulted in a significant increase (p < 0.01) in TAS compared to control rats. At the same time, the GR content decreased in the SAP group (p < 0.01). Terpenoid supplementation also caused a decrease in GR. The applied treatment did not have an effect on IL-6. IL-10 concentrations increased in both experimental groups, with the highest values reached in the TERP group. These changes were not statistically significant. Similarly, in the TERP and SAP groups, the increase in MCP-1 concentration was not statistically confirmed. Terpenoids did not have an effect on ROS values, while saponins caused a significant decrease (p < 0.01).  4 values are mean ± SD; different letters within a row indicate significant differences between groups (p < 0.05; ANOVA followed by the post hoc test); 5 TAS, total antioxidant status; 6 GR, glutathione reductase; 7 IL-6, interleukin-6; 8 IL-10, interleukin-10; 9 MCP-1, blood monocyte chemotactic protein-1; 10 ROS, reactive oxygen species; 11 NO, nitric oxide.
Previous studies have shown that metabolic syndrome characterized by obesity in Zucker rats decreased superoxide dismutase (SOD) activity in the plasma and heart, associated with cardiomyocyte hypertrophy [37]. Supplementation with saponins or terpenoids in our study increased TAS and decreased ROS, with no effect on blood NO concentration. As a consequence of obesity, the metabolic function of GPx was not affected to the same extent as the metabolic activity of SOD [37]. Both compounds tested caused a reduction in GR activity, but the effect of terpenoids was minor, and the effect of saponins was statistically significant (p < 0.01). Conversely, TNF-α and IL-6 are associated with a tendency to develop ischemia or atherosclerotic events [38]. In our study, there was no significant effect of saponins and terpenoids on blood concentrations of TNF-α and IL-6, and terpenoids caused an increase in IL-10. Interleukin 10 is considered an anti-inflammatory cytokine that inhibits the production of pro-inflammatory cytokines such as IL-2, IL-3, and TNF-α, for example. In an earlier study conducted on rats without metabolic syndrome (Wistar), TOS, GPx, and MDA were reduced after supplementation with terpenoids alone. No such effect was seen with water extracts of Ilex spp. [32].

Fatty Acid Profile in Zucker Rat Livers
The chromatogram obtained for the rat liver differed from the chromatogram of fatty acids identified in the leaves of European holly species, primarily in terms of the presence of long-chain acids. Analyses showed that palmitic, elaidic, stearic, arachidonic, and linoleic acids were the most abundant in all test groups (Table 8). Supplementation with the terpenoid fraction and saponins had a significant effect on the change in fatty acid content in the rat liver. A decrease in elaidic acid was observed in the TERP group and an increase in the SAP group; its high content is not desirable due to a decrease in HDL cholesterol. For linoleic acid, no differences were observed between the CON and TERP groups, but a decrease was observed in the SAP group. With respect to arachidonic acid, an increase was observed in the TERP group and a decrease in the SAP group. An analogous trend was also observed for docosahexaenoic acid, which may be due to the conversion of α-linolenic acid, the content of which decreases significantly in both study groups. Comparing the content of saturated and unsaturated acids, SFA acids predominated in all groups. However, in the SAP group, the content of MUFA acids was higher compared to the TERP group. Conversely, the content of PUFA acids in the group fed with the saponin fraction was clearly lower. In total, PUFAω6 and PUFAω3 acids in the SAP group were also observed in a lower proportion than for the CON and TERP groups. The hepatic lipid profile in Zucker rats has also been investigated in other works. González-Torres et al. [39], in their study's control group, showed the proportion of palmitic acid (16.31%), palmitoleic acid (20.40%), oleic acid (31.96%), arachidonic acid (9.78%) and docosahexanoic acid (3.45%) which together accounted for approximately 80% of the total lipid profile. Overall, the proportion of SFA and PUFA acids was similar, but the proportion of MUFA acids was twice as high. Conversely, Fiebig, who studied the effects of training on biological parameters in rats, found that in the group of untrained obese rats, the content of palmitic acid was 35.00%, with palmitoleic acid at 8.32%, oleic acid at 13%, and arachidonic acid at 6.78% [40]. In the study by Aguirre et al., feeding rats with pterostilbene caused an increase in total PUFAs, mainly through an increase in docosohexanoic acid and linoleic and arachidonic acids. The pterostilbene treatment decreased elaidic acid and the sum of MUFA [41].

Transcription of Genes in the Aorta and Liver
Increased TAS and decreased GR in the blood in both the SAP and TERP study groups indicate the plants' antioxidant effects, and increased IL-10 in the blood in both study groups and decreased ROS in the SAP group indicate their anti-inflammatory effects. These findings are correlated with the increased expression of LOX1 mRNA in the aorta (Figure 2). The expression of the receptor for oxidized LDL increases in the case of oxidative stress and inflammatory processes and is one of the most important markers of the increased oxidation of LDL particles with their subsequent incorporation into the arterial walls and formation of atherosclerotic plaques. Atherogenesis stimulating factors, TNF-α, oxLDL, and shear stress increase the expression of LOX1 [42]. In our study, we found an increase in the expression of NOS3 (eNOS) in the aorta in the SAP group (statistically significant) and in the TERP group (not statistically significant), arguing for increased endothelial nitric oxide synthase activity, leading to an increase in vascular endothelial NO, the main vasodilatory factor. Additionally correlated with these results was a decrease, although not statistically significant, in the expression of the metalloproteinase MMP1 mRNA, evident in the TERP group. The reduction, although not statistically significant, in the expression of mRMA NOX1 in the SAP group, indicates a reduction in oxidative stress in the aorta. There was also a reduction, although not statistically significant, in ACAT1 in both groups; its increased expression was described in atherosclerotic plaques. Finally, although not statistically significant, PPAR-α and PPAR-γ mRNA expression can be interpreted as positive. Rats have different metabolic pathways than humans or rabbits. Many authors interpret the reduction in PPAR expression in rats as a positive indicator associated, for example, in the case of PPAR-γ, with reduced lipogenesis. Conversely, PPAR-α plays a pivotal role in controlling fatty acid oxidation. Synthetic PPAR-α agonists are currently used in the clinic as lipid-lowering and anti-atherosclerotic drugs. However, the potencies of these compounds at PPAR-α are lower than those of, e.g., monounsaturated FAE oleoylethanolamide and consequently do not significantly affect food intake [43]. mans or rabbits. Many authors interpret the reduction in PPAR expression in rats as a positive indicator associated, for example, in the case of PPAR-γ, with reduced lipogenesis. Conversely, PPAR-α plays a pivotal role in controlling fatty acid oxidation. Synthetic PPAR-α agonists are currently used in the clinic as lipid-lowering and anti-atherosclerotic drugs. However, the potencies of these compounds at PPAR-α are lower than those of, e.g., monounsaturated FAE oleoylethanolamide and consequently do not significantly affect food intake [43].

Figure 2.
Expression of genes after treatment of saponin and terpenoid fraction in Zucker rat aorta. CON, control group; SAP, rats fed a standard diet with the addition of saponin fraction; TERP, rats fed a standard diet with the addition of terpenoid fraction. Values were normalized against Gapdh expression and presented as a fold-change in relation to gene expression in a control group. Specific comparisons: * p < 0.05 vs. CON; # p < 0.05 vs. SAP (ANOVA followed by the post hoc test). TNFα, tumor necrosis factor alpha; PPARα, peroxisome proliferator-activated receptor alpha; PPARγ, peroxisome proliferator-activated receptor gamma; LOX1, lectin-type oxidized LDL receptor; MMP1, matrix metalloproteinase-1; NOS3, nitric oxide synthase 3; NOX1, NADPH oxidase 1; ACAT1, acetyl-CoA acetyltransferase 1.
In the livers (Figure 3) of the TERP group, expression of the transcription factor LXR1, which regulates cholesterol and fatty acid metabolism and reduces inflammatory processes, was increased [44]. LXR1 regulates fatty acid homeostasis through, among others, SREBP-1c, and this was increased, although not statistically significantly, in the TERP group. Conversely, mRNA expression of HMGR, a gene encoding a key enzyme in the Figure 2. Expression of genes after treatment of saponin and terpenoid fraction in Zucker rat aorta. CON, control group; SAP, rats fed a standard diet with the addition of saponin fraction; TERP, rats fed a standard diet with the addition of terpenoid fraction. Values were normalized against Gapdh expression and presented as a fold-change in relation to gene expression in a control group. Specific comparisons: * p < 0.05 vs. CON; # p < 0.05 vs. SAP (ANOVA followed by the post hoc test). TNFα, tumor necrosis factor alpha; PPARα, peroxisome proliferator-activated receptor alpha; PPARγ, peroxisome proliferator-activated receptor gamma; LOX1, lectin-type oxidized LDL receptor; MMP1, matrix metalloproteinase-1; NOS3, nitric oxide synthase 3; NOX1, NADPH oxidase 1; ACAT1, acetyl-CoA acetyltransferase 1.
In the livers (Figure 3) of the TERP group, expression of the transcription factor LXR1, which regulates cholesterol and fatty acid metabolism and reduces inflammatory processes, was increased [44]. LXR1 regulates fatty acid homeostasis through, among others, SREBP-1c, and this was increased, although not statistically significantly, in the TERP group. Conversely, mRNA expression of HMGR, a gene encoding a key enzyme in the endogenous synthesis of cholesterol in the liver, decreased in the SAP group, although not statistically significantly (statins also act in this mechanism, although they have pleiotropic effects, also positively influencing many other functions). In the SAP group, there was a decrease in PPAR-α and PPAR-γ mRNA expression. The interpretation of these changes is analogous to that for the aorta. The expression of ACAT1 in the liver, responsible for the formation of cholesterol esters, decreased, but not statistically significantly, relative to the control group in the SAP group. endogenous synthesis of cholesterol in the liver, decreased in the SAP group, although not statistically significantly (statins also act in this mechanism, although they have pleiotropic effects, also positively influencing many other functions). In the SAP group, there was a decrease in PPAR-α and PPAR-γ mRNA expression. The interpretation of these changes is analogous to that for the aorta. The expression of ACAT1 in the liver, responsible for the formation of cholesterol esters, decreased, but not statistically significantly, relative to the control group in the SAP group.  Figure 3. Expression of genes after treatment of saponin and terpenoid fraction in Zucker rat liver. CON, control group; SAP, rats fed a standard diet with the addition of saponin fraction; TERP, rats fed a standard diet with the addition of terpenoid fraction. Values were normalized against Gapdh expression and presented as a fold-change in relation to gene expression in a control group. Specific comparisons: * p < 0.05 vs. CON; # p < 0.05 vs. SAP (ANOVA followed by the post hoc test). HMGR, 3-hydroxy-3-methyl-glutaryl-coenzyme A reductase; MMP9, matrix metallopeptidase 9; PPARalpha, peroxisome proliferator-activated receptor alpha; PPARγ, peroxisome proliferator-activated receptor gamma; ACAT1, acetyl-CoA acetyltransferase 1; ACC1, acetyl-CoA carboxylase 1; ACC2, acetyl-CoA carboxylase 2; LXR1, liver X receptor; SREB1c, sterol regulatory element-binding protein 1.

Histopathological Examination
In the CON group ( Figure 4A), histological analysis revealed deep liver steatosis in examined rats. The structure of the organ was blurred, and the portal bile spaces were invisible. The blood vessels were moderately filled with blood. There were cells of the white blood cell system within the larger blood vessels and in the periosteal spaces. Numerous apoptotic corpuscles were found in the periosteal spaces. There were no features Figure 3. Expression of genes after treatment of saponin and terpenoid fraction in Zucker rat liver. CON, control group; SAP, rats fed a standard diet with the addition of saponin fraction; TERP, rats fed a standard diet with the addition of terpenoid fraction. Values were normalized against Gapdh expression and presented as a fold-change in relation to gene expression in a control group. Specific comparisons: * p < 0.05 vs. CON; # p < 0.05 vs. SAP (ANOVA followed by the post hoc test). HMGR, 3-hydroxy-3-methyl-glutaryl-coenzyme A reductase; MMP9, matrix metallopeptidase 9; PPARalpha, peroxisome proliferator-activated receptor alpha; PPARγ, peroxisome proliferator-activated receptor gamma; ACAT1, acetyl-CoA acetyltransferase 1; ACC1, acetyl-CoA carboxylase 1; ACC2, acetyl-CoA carboxylase 2; LXR1, liver X receptor; SREB1c, sterol regulatory element-binding protein 1.

Histopathological Examination
In the CON group ( Figure 4A), histological analysis revealed deep liver steatosis in examined rats. The structure of the organ was blurred, and the portal bile spaces were invisible. The blood vessels were moderately filled with blood. There were cells of the white blood cell system within the larger blood vessels and in the periosteal spaces. Numerous apoptotic corpuscles were found in the periosteal spaces. There were no features associated with the severity of fibrinogenesis processes or fibroblast proliferation. Many hepatocytes showed double cell nuclei, which indicates an attempt to compensate for the processes associated with cell damage. The main change observed in the hepatocytes was the presence of accumulated lipids, which occupied most of the cell. These processes were observed most strongly in the peripheral areas of the lobules, while, in the area of the central vein, this process was much weaker and the hepatocytes contained a greater amount of cytoplasm. vessels were accompanied by leukocytic infiltrates.
In the TERP group ( Figure 4C), the presence of a reduced amount of lipids was observed. No dividing cells or blood steatosis were found. The size of the hepatocytes decreased compared to that of the control group. In the cytoplasm, small lipid droplets were observed.
In both experimental groups ( Figure 4E,F), an increase in the number of peroxisomes was observed compared to the CON group ( Figure 4D). Usually, this increase is not observed in the entire lobuli, but as numerous spots in some neighboring cells ( Figure 4E,D).

Plant Material
Fresh plant material in the form of leaves of I. aquifolium and I. meserveae was obtained from the University of Life Sciences' own cultivation (Vegetable and Ornamental Plant Changes within the SAP group ( Figure 4B) left a picture of fatty liver, but of lower intensity. Less blood vessels were filled with blood, and less swelling of the organ was observed. At the same time, the presence of portal bile spaces was observed between the hepatocytes. In the subcapsular areas, increased fibroblast activity was observed. Blood vessels were accompanied by leukocytic infiltrates.
In the TERP group ( Figure 4C), the presence of a reduced amount of lipids was observed. No dividing cells or blood steatosis were found. The size of the hepatocytes decreased compared to that of the control group. In the cytoplasm, small lipid droplets were observed.
In both experimental groups ( Figure 4E,F), an increase in the number of peroxisomes was observed compared to the CON group ( Figure 4D). Usually, this increase is not observed in the entire lobuli, but as numerous spots in some neighboring cells ( Figure 4E,D).

Plant Material
Fresh plant material in the form of leaves of I. aquifolium and I. meserveae was obtained from the University of Life Sciences' own cultivation (Vegetable and Ornamental Plant Research and Education Station, Psary, Poland). The holly shrubs were grown at the Department of Horticulture in the ornamental plant collection. After harvesting, the cleaned leaves were freeze-dried for 24 h (Lyovac GT 2 apparatus). The lyophilized leaves were crushed to obtain homogeneous material. The reference material was dried leaves of I. paraguariensis, purchased from a local distributor.

Solid-Phase Micro Extraction (SPME) Analysis
Fresh 1 g of material was ground in mortar, transferred to a headspace vial, and stored in a water bath at 40 • C for 15 min. Additionally, 2.5 µg of 2-undecanone was added as an internal standard. The volatile compound fraction was analyzed using a 2 cm fiber (DVB/CAR/PDMS, Supelco, Bellefonte, PA, USA); 15 min exposure and analyte desorption at 220 • C for 3 min. Analysis was performed on a Varian CP-3800/Saturn 2000 (Varian, Walnut Creek, CA, USA) equipped with a Zebron ZB-5 MSI column (30 m, 0.25 mm, 0.25 µm; Phenomenex, Torrance, CA, USA). The instrument was operated under the following conditions: injector 220 • C, gas flow as 1 mL · min −1 with a split ratio of 20. The program was set as follows: 40 • C for 3 min, a rate of 5.0 • C × min −1 from 40 to 110 • C and a rate of 20 • C × min −1 from 110 to 270 • C.

Isolation and Analysis of Terpenoid Fraction
The fractional analysis of natural compounds from Ilex leaves was performed using a gas chromatograph coupled to a mass spectrometer (Shimadzu GC-MS QP 2020 singlequadrupole gas chromatograph mass spectrometer, Shimadzu, Kyoto, Japan). Analyses were performed in triplicate.
To isolate the terpenoid fraction, 100 mg of Ilex leaves was weighed and poured into 5 mL of dichloromethane. An internal standard of 0.5 mg of cholesterol was also added to the sample. The sample was centrifuged on a benchtop centrifuge at 14,000× g rpm for 2 min. The supernatant was collected and washed twice with 4 mL of distilled water, which was then removed. The organic fraction was filtered through a celite filter, transferred to a round bottom flask, and evaporated to dryness on a vacuum evaporator under reduced pressure. Then, 400 µL of pyridine and 50 µL of N,O-Bis (trimethylsilyl)trifluoroacetamide (BSTFA) were added, and the mixture was transferred to a vial. The derivatization process was carried out in GC-MS whisks on a heating panel at 70 • C for 20 min. The prepared samples were stored at −18 • C until GC-MS analyzes were performed. Analyses were performed in triplicate.
Separation of the compounds that occur in the terpenoid fraction was achieved using a Zebron ZB-5 capillary column (30 m, 0.25 mm, 0.25 µm; Phenomenex, Torrance, CA, USA). GC-MS analysis was performed according to the following parameters. Scans were performed from 40 to 600 m/z using electron impact ionization (EI) at 70 eV, in a mode of 10 scans × s −1 . The carrier gas used was helium, with a flow rate of 23 mL × min −1 , with a split ratio of 1:20. The program setup for the analysis was as follows: 180 • C for 1 min, a rate of 5.0 • C × min −1 from 180 to 300 • C for 15 min. The injector was maintained at 280 • C.

Isolation and Analysis of Saponins Fraction
Extraction of the saponin fraction was performed using a modified methodology developed by Włodarczyk et al. [45]. Lyophilized leaves of I. aquifolium and I. paraguariensis (60 g) were macerated for 24 h at room temperature using 600 mL of cold 70% methanol (performed twice). The next step was to combine the extracts, which was filtrated, and a solution of 30 g of lead (II) acetate trihydrate in 70% methanol was added to precipitate the polar components such as phenols, peptides, and polysaccharides. The extract was centrifuged for 15 min at 3000 rpm. The residue containing the precipitated ballast was removed. Residual lead ions in the saponin fraction were removed with the addition of diluted sulfuric acid. The supernatant was diluted with distilled water to a concentration of 40% methanol and then centrifuged. Subsequently, the supernatant was applied to a 20 g acetadecyl SPE column (J.T. Baker), and the saponin-rich fraction was recovered with a minimal volume of pure methanol. The resulting fraction was concentrated on a vacuum evaporator, and the residue was lyophilized.

Isolation and Analysis of Fatty Acids in Leaves and Rat Liver
The Folch method was used to extract lipids from Ilex leaves and rat liver samples [46]. Samples of crushed leaves and freeze-dried liver in the amounts of 0.5 and 0.1 g, respectively, were poured over 10 mL of the solvent mixture chloroform:methanol (in ratio 4:1) and heptadecanoic acid as an internal standard. Samples were allowed to stand for 1 h under a cover at 25 • C (the operation was performed twice), and then the organic phase was filtered and transferred to a round-bottom flask to evaporate the solvents to dryness under reduced pressure on a rotary evaporator (Heidolph-VAP Core). The contents of the flask were hydrolyzed using 2 mL KOH/MeOH, and the entire flask was heated on a heating panel at 70 • C for 5 min. Then, 1.5 mL of 14% (V/V) BF 3 /MeOH was added as a catalyst for the fatty acid esterification reaction. After the extract was cooled, 1 mL NaCl and 4 mL hexane were added to the flask for phase separation. The mixture was shaken for 5 min by hand, and then the hexane fraction was collected. The organic fraction that already contained methyl acid esters was dried on anhydrous magnesium sulphate, filtered on a silica celite filter and transferred to round bottom flasks to concentrate the sample (to 1 mL) on a rotary evaporator. The thus-prepared sample was then transferred to a GC-MS vial and stored at −18 • C until GC-MS analysis.
Separation of fatty acid methyl esters isolated from Ilex leaves and rat liver was achieved on a Zebron ZB-FAME capillary column (60 m, 0.20 mm, 0.20 µm; Phenomenex, Torrance, CA, USA). The GC-MS apparatus was operated under the conditions listed below: scanning was set from 40 to 400 m/z in electron ionization (EI) at 70 • C, the carrier gas type was helium with a flow rate of 0.98 mL min −1 at ratio 1:20, the temperature program was as follows: 80 • C for 2 min; rate of 3.0 • C × min −1 from 80 • C to 180 • C; rate of 8 • C × min −1 from 180 • C to 240 • C for 4 min. The injector was maintained at 280 • C.

Analysis of GC-MS Data
The identification of all obtained components on GC-MS was performed by two independent methods: (1) comparison of the calculated retention indices (RI) by the retention index calculator (against n-alkanes or saturated fatty acids authentic standards) with the RI contained in the NIST 20 database or Lipids Library 1.0; (2) comparison of the obtained spectra with the NIST 20 databases. During data analysis, the following programs were used: AMDIS (v. 2.73) and GCMS Solution (v. 4.20).

Animal and Experimental Design
Six-week-old Zucker (Zucker fatty rats fa/fa) male rats (n = 24) were purchased (Charles River, Madison, WI, USA) and fed a standard food diet for two weeks to acclimate them to their new environment. Animals were housed in pairs under controlled conditions (22 ± 1.6 • C were maintained on a 12 h light:12 h darkness schedule with standard laboratory diet and water available ad libitum). At eight weeks of age, 24 Wistar rats were randomly assigned to 3 groups (n = 8 animals per group): rats fed standard feed (CON group), rats receiving 10 mg/kg body weight of the saponin fraction extracted from I. aquifolium (SAP group), and rats receiving 10 mg/kg body weight of supplemented terpenoid fraction (TERP group). The extracts were administered individually to animals per os. The duration of the experiment was 8 weeks. Rats were fed commercial pellets (containing 16% protein, 2.5% fat, and 7.5% maximum fiber and other nutritional additives, according to the AOAC methods, (2005)). Food intake was monitored, and rats were weighed twice per week. After 8 weeks of feeding, the rats were sacrificed one at a time from each of the four groups and blood and tissue samples collected.
All procedures regarding Zucker rats were in agreement with the Local Ethics Committee (decision no. 09/2019/P1, Wrocław, Poland) and associated guidelines from European Communities Council Directive (no. 86/609/CEE) governing animal protection. The experimental design adopted was due to the limitations imposed by the 3Rs rule.
The following parameters of antioxidant status were also determined. Total antioxidant capacity (TAS) in serum by the colorimetric method based on ABTS (2,2-azine-di-[3ethylbenzothiazoline sulfate]ta) with peroxidase, and glutathione reductase (GR) in whole blood by enzymatic method. Glutathione reductase catalyzes the reduction of glutathione (GSSG) in the presence of NADPH, which is oxidized to NADP+. The decrease in absorbance was measured at 340 nm. Rat Monocyte Chemotactic Protein 1 (MCP-1)-Immunoenzymatic test-sandwich ELISA principle using an antibody specific for rat MCP-1. Kit for the determination of reactive oxygen species (ROS) from MyBioSource Enzyme Immunoassay (San Diego, CA, USA), enzyme immunoassay based on ROS-ROS antigen interactions (immunosorbtion), and the HRP colorimetric detection system for the detection of ROS target antigens in samples. The ELISA kit is designed to detect native and non-recombinant ROS. Nitric oxide (NO) ready kit for the determination of nitric oxide from MyBioSource (San Diego, CA, USA). Colorimetric test based on the quantification of total NO2-/NO3-as a measure of the NO level. The reaction products can be measured with a colorimetric reading at 550 nm. These measurements were conducted using a Synergy fluorescence, luminescence, and absorbance reader from BioTek Instruments (Winooski, VT, USA).

Gene Symbol
Forward Primer Reverse Primer Technical replicates were averaged, and the geometric mean of all Cq values in all samples was calculated and subtracted from the individual sample Cq, yielding ∆Cq. Then, ∆Cq values were linearized by 2 ∆Cq conversion, followed by a normalization to GADPH expression. The obtained normalized relative quantity (NRQ) [50] was subjected to statistical analysis.

Histopathological Examination of Tissues
Samples were fixed in a 4% buffered formalin solution with a pH of 7.2-7.4 for 24 h. Then, the material was rinsed in running water for 24 h, dehydrated in a series of alcohols, and embedded in paraffin. Slides 7 µm thick were routinely stained with hematoxylin and eosin.
Slides after deparaffinization in xylene were rehydrated, permeabilized with 0.1% Triton™ X-100 for 15 min, and blocked with 2% BSA for 1 h at room temperature. The cells were labeled with PMP70 Rabbit Polyclonal Antibody (Invitrogen, Waltham, MA, USA) at 5 µg/mL in 0.1% BSA, incubated at 4 • C overnight and then labeled with Goat anti-Rabbit IgG (H + L), Superclonal™ Recombinant Secondary Antibody, Alexa Fluor 488 (Invitrogen, Paisley, UK). The material was analyzed on a Zeiss Axio Observer D1 inverted phase contrast fluorescence microscope.

Statistical Analysis
Statistical analyses were performed with Statistica version 13.3 (StatSoft, Kraków, Poland). In vivo samples were analyzed in duplicate; the values for each group represent the mean of eight animals (n = 8). All results were expressed as means and standard deviation, and all variables were normalized using the Shapiro-Wilk test. The data obtained were subjected to analysis of variance (ANOVA) using a linear model procedure. The Institutional Review Board Statement: The animal study protocol was approved by the Local Ethics Committee (decision no. 09/2019/P1, Wrocław, Poland) and associated guidelines from European Communities Council Directive (n. 86/609/CEE) governing animal protection.