Biological Screening of Glycyrrhiza glabra L. from Different Origins for Antidiabetic and Anticancer Activity

Background: Geographical variation may affect the phytochemistry as well as the biological activities of Glycyrrhiza glabra (licorice) root. Herein, a series of biological activities were performed to evaluate the impact of geographical origin on the biological potential of eight different licorice samples. Methodology: Cell culture studies were performed for cytotoxicity (MCF7, HCT116, HepG2, and MRC5), glucose uptake assay (HepG2), and glutathione peroxidase activity (HepG2), whereas α-amylase inhibition activity was tested for antidiabetic potential. Results: The Indian sample was observed to be more cytotoxic against MCF7 (22%) and HCT116 (43%) with an IC50 value of 56.10 (±2.38) μg/mL against the MCF7 cell line. The glucose uptake was seen with a mean value of 96 (±2.82) and a range of 92–101%. For glutathione peroxidase activity (GPx), the Syrian (0.31 ± 0.11) and Pakistani samples (0.21 ± 0.08) revealed a significant activity, whereas the Palestinian (70 ± 0.09) and Indian samples (68±0.06) effectively inhibited the α-amylase activity, with the lowest IC50 value (67.11 ± 0.97) μg/mL for the Palestinian sample. The statistical models of PCA (principal component analysis) and K-mean cluster analysis were performed to correlate the geographical origin, extract yield, and biological activities for the eight licorice samples of different origins. Conclusion: The licorice samples exhibited significant cytotoxic, GPx, and α-amylase inhibitory activity. The samples with higher extract yield showed more potential in these biological activities.


Introduction
Glycyrrhiza glabra is commonly known as licorice; it is considered the most important species of the Glycyrrhiza genus [1,2]. It is widely distributed globally, mostly in Asia, but it is also found in Europe, Africa, and Australia [3]. Glycyrrhiza glabra has considerable medicinal, pharmaceutical, and industrial value [4]. The roots of Glycyrrhiza glabra have been particularly well known for their medicinal effects since ancient times [5]. The populations of Rome, Greece, India, and China, use the roots to treat respiratory ailments such as asthma and bronchitis [5,6]. Consequently, pharmaceutical companies use Glycyrrhiza glabra in many cough syrup preparations [3]. Moreover, the strong sweet taste of the roots is of significant interest to industrial bodies. It acts as a sweetener in many tobacco products, chewing gums, ice cream flavors, and candies [3,7]. Notably, only some species are very sweet [7]. This sweet taste of Glycyrrhiza glabra is linked to the level of a triterpenoid saponin known as glycyrrhizin. It is considered the primary phytochemical and most abundant compound (10-25%) of the Glycyrrhiza glabra root extract [1,7]. Additionally, it has been used as a quality marker in the pharmacopeias of countries such as Japan and China [8,9]. In the plant, it is present as glycyrrhizic acid salt [10]. In several studies, geographical, environmental, and genetic factors were found to influence phytochemical levels, including that of glycyrrhizic acid [1]. Zhang et al. tested the glycyrrhizic acid content of five samples from different locations in China and found significant glycyrrhizic acid variations between the samples. The soil components were proposed as the most crucial factor in the abovementioned variation [11]. Similarly, the climate temperature influenced the glycyrrhizic acid levels of various Iranian Glycyrrhiza glabra roots from many areas in samples investigated by Hosseini et al. [12]. A substantial amount of literature has been published on the Glycyrrhiza glabra extracts and their phytochemical cytotoxic activities. For example, the ethanolic root extract of Glycyrrhiza glabra suppresses the growth of the (MDA-MB-231) breast cancer cell line in a dose-dependent manner [13]. Equally, the proliferation of the (T47D, MCF7, MDA-MB-231, and MDA-MB-361) breast cancer cell lines, the SiHa cervical cancer cell line, and the (A2780) ovarian cancer cell line was notably inhibited by the ethyl acetate root extract. A considerable amount of 18-β glycyrrhetic acid was found when a phytochemical analysis was performed [14]. The compound showed anticancer activity against A549 lung cancer cells when examined by Luo et al. [15]. Additionally, glycyrrhizic acid was also found to be cytotoxic in numerous studies. For instance, it showed significant antileukemic activity in K562 chronic myeloid leukemia cells and EL-4 lymphoma-bearing mice [16]. Furthermore, it induced apoptotic cell death in gastric cancerous cells [17]. Similar results were found in a study by Cai et al., when hepatic cancer cells were used [18]. The exact mechanism of action is unclear, but several mechanisms have been proposed for this anticancer activity [3]. Other biological activities such as anti-inflammatory, antioxidant, antimicrobial, and antiviral have been proven for Glycyrrhiza glabra extracts or its derivatives in numerous studies [1,[19][20][21][22]. However, few studies have discussed the antidiabetic effects. For example, Takii et al. studied the effect of glycyrrhizin administration on non-insulin-dependent diabetic rats. After eight weeks, the blood glucose levels were suppressed remarkably in the treated rats compared to controlled non-insulin-dependent diabetic rats (p < 0.05). Likewise, the postprandial glucose level dropped notably after 60 min of glucose administration [23]. Similar findings were obtained by Sen et al. when the streptozotocin-induced diabetic model was used [24]. Glabridin could be suggested as the chemical constituent responsible for antidiabetic activity [25].
To the best of our knowledge, there are no studies on the correlation between the quality and origin of Glycyrrhiza glabra and its cytotoxic and antidiabetic activities. Thus, we attempt to investigate that correlation by using glycyrrhizic acid level as a quality marker.

Results
For descriptive analysis, principal component analysis, and K-mean cluster analysis, SPSS (Statistical Package for the Social Sciences) Version 22.0 was used for data entry and correlation analysis. A multi-view (three-dimensional) representation of the extract yield and GA amount versus the pharmacological activities (cytotoxicity, GPx, and glucose uptake assay) is presented in Figure 1.

Cytotoxicity Assay
The cytotoxicity for these samples was observed with a mean value of 56.25 (±15.41) for MCF7 and 75.50 (±14.60) for HCT116, whereas the ranges (%) observed were 22-74 and 43-90 for MCF7 and HCT116, respectively. The lowest % viability for both cell lines was observed in the Indian-origin sample, i.e., 22% (MCF7) and 43% (HCT116), as shown in Table 1. This sample, upon further MTT assay investigation (six different concentrations for IC50) in MRC5, MCF7, and HCT116 cell lines revealed the lowest IC50 value of 56.10 (±2.38) μg/mL against the MCF7 cell line, as shown in Table 2.   A mean value of 96 (±2.82) with a range (%) of 92.00-101.00 was observed for the glucose utilization assay (HepG2) in these eight different geographical samples. The glucose uptake assay revealed the lack of any significant (p < 0.05) increase in glucose uptake and its utilization in HepG2 cells (100 µg/mL), as shown in Table 1. The standard drug Metformin showed a value of 118% (±0.07).

The Effect on Glutathione Peroxidase Activity (GPx)
The GPx activity (100 µg/mL) revealed a mean value of 0.64 (±0.29) within the range of 0.21-0.93. The study revealed a significant GPx activity for the Syrian (0.31 ± 0.11; p = 0.01) and Pakistani samples (0.21 ± 0.08; p = 0.01) as compared to the standard drug Quercetin (0.36 ± 0.10; p = 0.01). The GPx activity results for all the samples are shown in Table 1.

The Effect of Extracts on α-Amylase Activity
The α-amylase activity was observed with a mean of 63.12 (±0.29) and a range (%) of 0.21-0.93. The initial screening for these samples (500 µg/mL) exhibited a significant inhibition for α-amylase (>50%); hence, they were further investigated at six different concentrations in order to determine the IC 50 value. The samples observed with the highest inhibition during initial screening were the Palestinian (70 ± 0.09) and Indian samples (68 ± 0.06), as shown in Table 1. For IC 50 value determination, the lowest value was observed for the Palestinian (67.11 ± 0.97), followed by the Indian sample (74.87 ± 1.26) µg/mL. The IC 50 values observed for all the samples are reported in Table 3. Table 3. α-Amylase inhibitory activity of the selected extracts (IC 50 ± SD µg/mL).

Principal Component Analysis
The PCA data for Glycyrrhiza glabra suggested three components, as shown in the scree plot Figure 2 and Table 4. The three components showed an individual and cumulative %variability of PC1 29.774 (29.774), PC2 28.607 (58.381), and PC3 23.467 (81.847). PC1, with more %variability, was loaded with cytotoxicity activities (HCT116 and MCF7) only. This reveals a high %variation of cytotoxic activities in relation to other biological activities, as well as a lack of correlation with the geographical origin, extract yield, and GA amount in these samples. As shown in Table 1, all the geographical origin samples exhibited cytotoxic potential with a similar pattern; however, the Indian and Georgian samples were significantly more cytotoxic. Yet, none of these samples had a high extract yield or GA amount either. PC2, with a variability of 28.607%, was loaded with the geographical origin, extract yield, and GA amount along with the α-amylase activity. This shows α-amylase activity to be more correlated with extract yield and GA amount, as all the Glycyrrhiza glabra samples of different geographical origins exhibited significant α-amylase activity. The last component, i.e., PC3 with a variability of 23.467%, was loaded with the geographical origin, glucose uptake, as well as GPx activity. PC3 suggests a higher correlation between glucose uptake and GPx activity, and geographical origin. For glucose uptake, more activity was shown by the geographical samples with either a high extract yield (Morocco, 99 ± 0.05; Palestine, 96 ± 0.09) or GA amount (Syria, 101 ± 0.07; Egypt, 96 ± 0.07). Likewise, GPx activity was observed to be more significant in the Pakistani and Syrian samples (0.21 ± 0.08 and 0.31 ± 0.11) with a high GA amount (p = 0.01). A graphical representation is also shown in Figure 3.
yield or GA amount either. PC2, with a variability of 28.607%, was loaded with the graphical origin, extract yield, and GA amount along with the α-amylase activity. shows α-amylase activity to be more correlated with extract yield and GA amount, a the Glycyrrhiza glabra samples of different geographical origins exhibited significan amylase activity. The last component, i.e., PC3 with a variability of 23.467%, was loa with the geographical origin, glucose uptake, as well as GPx activity. PC3 sugges higher correlation between glucose uptake and GPx activity, and geographical origin. glucose uptake, more activity was shown by the geographical samples with either a h extract yield (Morocco, 99 ± 0.05; Palestine, 96 ± 0.09) or GA amount (Syria, 101 ± 0 Egypt, 96 ± 0.07). Likewise, GPx activity was observed to be more significant in the P stani and Syrian samples (0.21 ± 0.08 and 0.31 ± 0.11) with a high GA amount (p = 0.01 graphical representation is also shown in Figure 3.

K-Mean Cluster Analysis
The K-mean suggested three clusters for Glycyrrhiza glabra samples, where all the samples were classified on the basis of extract yield and GA amount vs. activities. The three clusters consist of cluster (samples) 1(3), 2(4), and 3(1). Cluster one consisted of three Pharmaceuticals 2023, 16, 7 6 of 13 samples with a higher GA amount reported (Syria, Pakistan, and Egypt), whereas cluster 2 consisted of four samples with a higher amount of extract yield (America, Palestine, Georgia, and Morocco). Figure 4 represents the cluster numbers/geographical origin samples with the corresponding activities exhibited by these samples. As evident, cluster 1, i.e., with a higher amount of GA, showed the highest glucose uptake activity, with a minimal activity potential for cytotoxicity (HCT116, MCF7). Cluster 2, i.e., with higher extract yield samples, exhibited almost all the activities except glucose uptake, whereas cluster 3 (one sample only i.e., Indian origin) exhibited a high potential for α-amylase activity. The K-mean cluster further supports the data from PCA analysis, where the GA amount was observed with a negative correlation for most activities. Hence, it may be concluded that the Glycyrrhiza glabra extracts from various geographical origins contain numerous phytochemicals which play an important role in different activities. The presence of GA in high amounts is not a guarantee that a sample will exhibit all the biological activities. The cluster distributions with the samples, F-values, and p-values are shown in Table 5.

K-Mean Cluster Analysis
The K-mean suggested three clusters for Glycyrrhiza glabra samples, where all the samples were classified on the basis of extract yield and GA amount vs. activities. The three clusters consist of cluster (samples) 1(3), 2(4), and 3(1). Cluster one consisted of three samples with a higher GA amount reported (Syria, Pakistan, and Egypt), whereas cluster 2 consisted of four samples with a higher amount of extract yield (America, Palestine, Georgia, and Morocco). Figure 4 represents the cluster numbers/geographical origin samples with the corresponding activities exhibited by these samples. As evident, cluster 1, i.e., with a higher amount of GA, showed the highest glucose uptake activity, with a minimal activity potential for cytotoxicity (HCT116, MCF7). Cluster 2, i.e., with higher extract yield samples, exhibited almost all the activities except glucose uptake, whereas cluster 3 (one sample only i.e., Indian origin) exhibited a high potential for α-amylase activity. The K-mean cluster further supports the data from PCA analysis, where the GA amount was observed with a negative correlation for most activities. Hence, it may be concluded that the Glycyrrhiza glabra extracts from various geographical origins contain numerous phytochemicals which play an important role in different activities. The presence of GA in high amounts is not a guarantee that a sample will exhibit all the biological activities. The cluster distributions with the samples, F-values, and p-values are shown in Table 5.

Discussion
This is the first study to report the phytochemical profile and evaluate the quality of Glycyrrhiza glabra samples of different geographical origins using various biological activities. The phytochemical quantification for GA amount was conducted in our previously published study [26], and the extracts thus quantified were subjected to multiple biologi-

Discussion
This is the first study to report the phytochemical profile and evaluate the quality of Glycyrrhiza glabra samples of different geographical origins using various biological activities. The phytochemical quantification for GA amount was conducted in our previously published study [26], and the extracts thus quantified were subjected to multiple biological assays in order to correlate the Glycyrrhiza glabra phytochemistry with its biological potential. Eight Glycyrrhiza glabra samples of different origins were found and collected from Al-Khobar, an eastern province of Saudi Arabia. The green extraction and quantification were followed by a set of biological activities, including cytotoxicity, GPx, glucose uptake, and α-amylase activity. Statistical models were used to establish the correlation between the extract yield, GA amount, and biological activities.
The initial screening for cytotoxicity (100 µg/mL) revealed a significant potential against HCT116 and MCF7, especially in the samples with a higher GA amount (Syria, Pakistan, and Egypt). However, the sample of Indian origin demonstrated a significantly higher cytotoxic effect against the HCT-116 and MCF7 cell lines. In order to determine the IC 50 value for cytotoxicity, the Indian sample was studied further against the HCT116, MCF7, and MRC5 cell lines with different dose ranges ( Table 2). In the context of GA amount, the Indian sample showed a lower amount as compared to other samples, which suggests that the cytotoxic effects may not be linked to GA content only. A number of other chemical constituents in the extract may play a vital role in cytotoxic and anticancer effects, including licochalcone-A [27][28][29][30] glycyrrhizinic acid [31], 18b-glycyrrhetinic acid [31][32][33][34][35][36], Isoliquiritigenin [37], formononetin [38], and glabridin [39]. 18b-glycyrrhetinic acid, one of the derivatives of glycyrrhetinic acid, is present in a higher concentration in the Glycyrrhiza glabra root as compared to 18α-glycyrrhetinic acid [40,41], being reported with a dose-dependent cytotoxic effect against the colon cancer cell line HCT-116 [32]. Likewise, glycyrrhetinic acid and glycyrrhizic acid have been reported to induce apoptosis in cancer cells [31,42,43] with an established cytotoxic potential of Glycyrrhiza glabra root, as previously reported [44][45][46].
Glucose uptake activity is essential for studying the metabolic diseases of cancer, diabetes, myocardial infarction, etc. [47]. A number of pathways are involved in the regulation of glucose in cancer cells, posing hurdles in the pathway of PI3K signaling or the selective blockage of glucose transporter type1 (GLUT1), thereby resulting in low glucose uptake [48,49]. Most of the available drugs are known to exert the anticancer effect via targeting GLUT1 to deprive tumor cells of glucose [48,50]. Herein, although the licorice samples of Moroccan and Syrian origin revealed a high glucose uptake, the comparative activity was lower than for the standard drug metformin. A dose-dependent stimulation of glucose uptake has been reported for licorice constituents such as glabridin [25] and glycyrrhizin [51].
Glutathione peroxidase (GPx) plays an important role in protection from oxidative damage, where it converts reduced glutathione to oxidized glutathione. An abnormal level of GPx has been reported in free radical-related disorders [52,53]. In our study, liver cell line lysates revealed a significant GPx activity for the Syrian, Pakistani, and Indian-origin samples as compared to the standard drug used. To the best of our knowledge, this assay is performed for the first time using Glycyrrhiza glabra extracts of different geographical origins. GPx is known for its complex effect on the development and progression of cancer due to its vital role in modulating the intracellular ROS [54]. An imbalance in essential ROS may lead to lessened cell growth with the encouragement of apoptotic pathways, whereas a surplus of cellular ROS plays a vital role in facilitating the multifaceted cascade of events leading to cell death via apoptotic pathways. For GPx, its high level may prevent oxidative damage and inflammation, although it may also block apoptotic cell death, hence leading to a higher survival rate for altered cells. This indicates the complex role of GPx, i.e., excess of GPx may allow for tumor survival and growth that may be reduced by the intake of suitable natural alternatives, such as licorice. Several in vivo toxicity studies already confirmed the role of licorice in normalizing or restoring the GPx level in the toxicity studies models [55][56][57][58]. The inhibition of carbohydrate digesting enzyme, i.e., α-amylase, is one of the investigative targets for treating diabetes [59]. The Glycyrrhiza glabra samples from Palestine and India revealed a significant α-amylase inhibitory activity as compared to the standard drug acarbose. Previous studies have reported a more effective inhibition of α-amylase than acarbose [60], which is in line with our current results.
The statistical models of PCA and K-mean cluster analysis of correlation revealed more activities for the licorice samples with higher extract yields. The samples with higher amounts of GA were found to have significant glucose uptake activity only, whereas the remaining tested activities were found to be significantly better in the samples with higher extract yields. This may be due to the presence of various other chemical classes present in these extracts, which imparts potential to these samples as observed in the cytotoxic, GPx, and α-amylase activities. Hence, it may be concluded that the presence of one chemical class or constituent may not be claimed for all the biological or pharmacological activities. The presence of multiple chemical constituents or phytochemical classes may also result in a synergistic effect, hence showing comparatively more significant results when tested in various throughput screenings.

Extraction and Phytochemical Analysis
The Glycyrrhiza glabra samples of different geographical origins (Syria, Egypt, America, Pakistan, India, Palestine, Georgia, and Morocco) were extracted with the help of the ultrasonication technique using water as a green solvent. The green analytical method of UHPLC-MS/MS was applied to quantify the amount of glycyrrhizic acid (GA) in each sample. The details regarding the optimization and validation for green extraction and analysis of GA in these samples were reported in our previous study [26]. For convenience, the extract yield and GA amount/sample are reported in Table 1 herein.

Determination of Cytotoxicity and Selectivity
The cytotoxicity of the extracts was evaluated by MTT assay, as previously reported [61]. The MCF7, HCT116, and MRC5 cell lines were tested to determine the cytotoxicity and selectivity of the extracts, while HepG2 was used to determine the effective concentrations of the extracts. Concentrations that did not significantly influence cell viability were further investigated in other assays. The MCF7 and HCT116 cells were cultured separately in 96-well plates (3 × 10 3 /well) and incubated at 37 • C overnight. The first set of experiments screened the samples for cytotoxic activity at 100 µg/mL concentration (DMSO 0.4%; n = 3). Plates were incubated for 48 h, after which MTT was added to each well and the plates were incubated further for 3 h. The supernatant was removed, and the MTT product was solubilized by adding DMSO to each well. Absorbance was read using a multi-plate reader (BIORAD, PR 4100, Hercules, CA, USA). The optical density of the purple formazan A550 was proportional to the number of viable cells, which was calculated as an inhibition percentage compared to control cells and listed in Table 1. The extracts with the highest percentage of inhibition were selected for the second set of experiments to determine the half-maximal inhibitory concentration (IC 50 ) and selectivity using the same cancer cell lines and the same fibroblast cell line. The range of extract concentrations tested were 500, 250, 100, 50, 10, and 1 µg/mL, and doxorubicin was used as a positive control to determine the IC 50 values using GraphPad Prism (San Diego, CA, USA) ( Table 2).

Glucose Uptake Assay
Glucose uptake was evaluated for the Glycyrrhiza glabra extracts as previously outlined by Odeyemi et al. [62]. After seeding HepG2 cells (5000 cells/well) in a 96-well plate and growing them overnight in a 37 • C humidified incubator with 5% CO 2 , cells were exposed to extract concentrations that had no influence on cell viability and incubated at 37 • C for 48 h, with 0.2 µg/mL metformin as a positive control. The culture medium was then aspirated, and the cells were incubated with an incubation medium containing DMEM diluted with 8 mM glucose, 0.1% bovine serum albumin (BSA), and PBS for 3 h at 37 • C. The assay was carried out by transferring the incubation medium to a new 96-well plate to measure the glucose concentration in the medium, following the manufacturer's protocol (GAGO20, Sigma Aldrich, St Louis, MO, USA). Absorbance was measured at 550 nm on a multi-plate reader (BIORAD, PR 4100, Hercules, CA, USA). We calculated the amount of glucose used by the cells by subtracting the cell-containing wells from the cell-free wells.

Determination of Glutathione Peroxidase Activity
The activity of glutathione peroxidase (GPx) was evaluated using HepG2 cell lysates in a colorimetric assay, performed according to the manufacturer's protocol (ab102530, Abcam, Cambridge, UK). HepG2 cells (5 × 10 5 cells/ well) were seeded in a 12-well plate and treated with extract concentrations that do not influence cell viability. These were incubated for 48 h, with 25 µg/mL of quercetin used as a positive control, after which the cells were collected and washed with cold PBS. The collected cells were depleted of all GSSG by incubating the samples with glutathione reductase (GR) and reduced glutathione (GSH) for 15 min. GPx activity was determined by adding cumene hydroperoxide and incubating for 5 min. The absorbance was measured before and after adding cumene hydroperoxide at 340 nm, using a multi-plate reader (BIORAD, PR 4100, Hercules, CA, USA). The decrease in NADPH (measured at OD = 340 nm) is proportional to GPx activity. The concentration of GPx in the test samples was calculated as recommended by the manufacturer's protocol.

α-Amylase Inhibition Activity
All extracts were initially tested at 500 µg/mL, and those that demonstrated an inhibitory effect on the α-amylase enzyme were further investigated at a range of concentrations (1000, 500, 100, 50, 25, and 10 µg/mL), with acarbose as a positive control. The inhibitory activity of α-amylase was ascertained as described by Quan et al. [63]. An aliquot of 1 mg α-amylase from Aspergillus oryzae in a phosphate buffer was prepared, and 20 µL was added to each well in a 96-well plate along with 20 µL of the extract samples that were diluted in a phosphate buffer. After mixing, the plate was incubated for 10 min at 37 • C, then 30 µL of starch (0.05% in deionized water) was added to each reaction well and incubated for a further 8 min at 37 • C. The reaction was then halted by the addition of 20 µL of hydrochloric acid (1M) followed by 100 µL of iodine reagent (0.25 mM) to each well. Control wells were set up by replacing the enzyme with buffer and adding acarbose to create a positive control. Absorbance was measured using a multi-plate reader (BIORAD, PR 4100, Hercules, CA, USA) at 550 nm for each well, calculating the percentage of inhibition utilizing the equation: % inhibition = (A−C/B−C) × 100, where A = the absorbance of the reaction mixture in the presence of the extract, B = the absorbance of the mixture without the enzyme, and C = the absorbance of the reaction mixture in the absence of any extract.

Statistical Analysis
The results are expressed as the mean ± standard deviation (SD) from at least three independent experiments. All data were analyzed with a p-value of < 0.05 considered as significant in order to find the statistical significance between treated groups and controls using GraphPad Prism 9.2.0 (GraphPad, San Diego, CA, USA).

Conclusions
The tested Glycyrrhiza glabra samples revealed a significant potential in different biological activities, including cytotoxicity, glucose uptake, GPx, and α-amylase activity. The samples were effectively classified into groups consisting of higher extract yield and higher GA amount samples. The statistical analysis established a higher correlation for samples with higher extract yield vs. biological activities, as compared to the samples with a higher GA amount. Hence, it is of utmost importance to further investigate and explore the complete metabolomics profile along with in vitro and in vivo studies for all these Glycyrrhiza glabra samples in order to reveal the presence of phytochemical classes. This may provide a more effective data source for a clear demarcation of the quality of Glycyrrhiza glabra samples based on geographical origin, which is sensitive to various environmental, agricultural, storage, and transport factors.  Institutional Review Board Statement: This study does not require any ethical approval as no human or animal subjects were used.

Informed Consent Statement:
The authors agree to publish the article in its current format.

Data Availability Statement:
Any data generated during this study is presented in this document.