1. Introduction
The developed world has seen a sharp increase in the rate and prevalence of diseases associated with the localized and general loss of bone mass. Numerous studies have reported that many entities and medications are associated with this loss [
1,
2,
3,
4,
5,
6], which together with an increasingly sedentary lifestyle have a profound impact on bone metabolism [
7]. Serious bone lesions are difficult to redress through autologous processes, which means that therapeutic procedures need to be used instead, but these are still limited [
8].
Tissue engineering and cell therapies provide alternatives for inducing and promoting bone regeneration, focusing on personalized and precision regenerative medicine. Among other aspects, this reduces the morbidity of other treatments, such as surgery. The availability of biomaterials with an appropriate physical structure, adhesion potential, and cell proliferation, plays an influential role in the regeneration and repair of bone tissue [
9,
10,
11], with a crucial factor being the advanced state of knowledge on stem cells [
12,
13,
14] and on the growth factors that favor the differentiation and acceleration of these processes [
15,
16].
Given the scarce bioavailability of autologous grafts, the low osteoinductive capacity of xenografts appears to be the ideal solution. Cells, scaffolds (three-dimensional matrices or supports), and growth factors are the key factors [
17], besides an appropriate physical-chemical environment. This requires the introduction of standardized protocols, from the gathering of cell material through to the establishment of clinical criteria, and this is the path taken by numerous studies.
This study assesses an experimental model of bone regeneration involving a mandibular defect in rats, using mesenchymal stem cells (MSCs) obtained from pulp material for finding a suitable biological substitute for regenerating critical bone defects and increasing the potential of MSCs through biological supports with a high concentration of osteogenic factors.
2. Results
2.1. Obtaining Dental Pulp Mesenchymal Stem Cells (DP-MSC)
The DP-MSC were successfully obtained from the molar pulp chamber (Image 1).
Image 1.
Extraction of pulp.
Image 1.
Extraction of pulp.
2.2. Cell Culture and Expansion
The results show that MSCs can be isolated by applying enzymatic digestion methods to dental pulp tissue. The cells obtained from the culture had a long and flat fibroblast-type morphology and were attached to the plastic surface on the culture, as essential characteristics of MSCs (Image 2).
Image 2.
Confluence of MSCs under optical microscopy.
Image 2.
Confluence of MSCs under optical microscopy.
2.3. Characterisation. Flow Cytometry
The cultures from each specimen were investigated for the percentage of MSCs by flow cytometry, using the following panel: CD49e PE (HMa5-1, BioLegend® cat. 103905), 7-AAD (BioLegend® cat. 420403), CD11b/c PE-Cy7 (OX-42, BioLegend® cat. 201817), CD29 APC (HMb1-1, Miltenyi Biotec cat. 130-123-829), CD90-1 APC-/Fire™ 750 (Thy1.1, BioLegend® cat. 202543), and CD45 Pacific Blue™ (OX-1, BioLegend® cat. 202225). The samples were acquired in the BD FACSCanto™ II System and the data were analysed using FlowJo™ Software (v. 10, ©Becton Dickinson & Company).
The gating strategy used for analyzing the MSC population is illustrated in
Figure 1. The percentage of MSCs was calculated from live cells that did not express lineage markers (CD45 and CD11b/c) and which expressed CD90 and CD29 simultaneously.
In total, from six rats, this population of MSCs was identified in 94.58% ± 3.62 of the isolates. From step two all the samples already presented 90% of the population with the characteristics of the MSCs (Lineage-CD90-1+CD29+) and were maintained until step 10.
2.4. Verifying the Cell Viability of the Implants
Following the surgery, the remains of the vials were seeded on plates. Within a few days’ cells appeared attached to the plate, and were therefore viable and proliferating.
2.5. In Vivo Results in the Experimental Model
An essential step in our assay involved creating a critical bone defect on the mandible that was large enough to avoid spontaneous regeneration. A critical defect is defined as one that regenerates less than 10% of the bone during the animal’s lifetime [
18], and in a rat mandible this will mean an area of more than 4 mm in diameter [
19,
20,
21].
None of the specimens died in the post-operatory period, nor recorded fractures, with a mandibular defect of 5 mm in diameter in all cases.
The overall condition of the specimens was good in all the groups, with some swelling of soft tissue due to the osteotomy, which disappeared spontaneously.
The post-surgical clinical monitoring of the specimens treated did not reveal any rejections or degenerative changes due to the dental lesion. The membrane recorded bioresorption, low immunogenicity, stability, and prolonged protection of the covered graft.
2.5.1. CONTROL Group
Study with microcomputed tomography (microCT)
Image 3 shows the defect on the mandible in the Control group, prior to its resection at three (A) and six (B) months of the study.
Image 3.
MicroCT image of a defect in Group C at 3 (A) and 6 months (B).
Image 3.
MicroCT image of a defect in Group C at 3 (A) and 6 months (B).
2.5.2. G+M Group
Study with microcomputed tomography (microCT)
At three months, we observed a significant collapse of the soft tissues and a minor restitution of the defect. The biomaterial has not diminished in volume. The radiopacity of the lesion is highly heterogenous in the center of the defect (Image 4 C).
At six months (Image 4 D), a highlight is the geometric morphology of the defect, which has become smaller, acquiring a denser appearance and with significant irregularity in its perimeter and a significantly greater radiological density than in the preceding state and greater uniformity, although the granulated aspect remained. There is an attempted formation of trabeculae between the biomaterial and the cortical, with the hazy appearance of dense bone pathways in the defect treated with the membrane + scaffold, compared to the situation at three months (Image 4 C).
Image 4.
MicroCT image of a defect in the G+M Group at 3 months (C). CT image of a defect in the G+M Group at 6 months (D).
Image 4.
MicroCT image of a defect in the G+M Group at 3 months (C). CT image of a defect in the G+M Group at 6 months (D).
2.5.3. G+M+S Group
Study with microcomputed tomography (microCT)
At three months, uniform levels of radiopacity were recorded, with significant consistency in the perimeter and partial structural restoration with new bone (Image 5 E). At six months, there was almost complete radiological repair of the defects in all cases (80-100%), revealing a high consistency with the bone perimeters and a very uniform distribution of the radiological density (Image 5 F).
Image 5.
MicroCT image of a defect in the G+M+S group at 3 months (E). Radiological image of a defect in the G+M+S group at 6 months (F).
Image 5.
MicroCT image of a defect in the G+M+S group at 3 months (E). Radiological image of a defect in the G+M+S group at 6 months (F).
2.6. Descriptive Histological Analysis
The untreated Control group recorded a fibrous cicatricial reaction, sometimes overrunning the edges of the wound, with no signs of bone regeneration in either of the two timeframes studied (Image 6).
Image 6.
Fibrous reaction of an untreated lesion. The fibrous cicatricial reaction overran the edges of the lesion. There is no sign of any reactive bone regeneration.
Image 6.
Fibrous reaction of an untreated lesion. The fibrous cicatricial reaction overran the edges of the lesion. There is no sign of any reactive bone regeneration.
At three months, one of the treated groups, G+M (A), showed few signs of bone regeneration, and the edges of the lesion were readily observable. Particles of biomaterial could be seen together with chronic inflammatory infiltration with the presence of multinucleated giant cells (Image 7).
For the same timeframe, the other treated group G+M+S (B) revealed major cellularity, with a considerable decrease in inflammatory reaction bone-forming activity. It recorded osteoclastogenesis and neoformed bone trabeculae surrounding the biomaterial. There are active osteoblasts and, in some cases, neoformed bone trabeculae around the biomaterial. There are areas with osteoblast precursor cells and osteoblasts partially included in the osteoid, which they are calcifying. There is also abundant immature reticular bone together with osteocytes in the lines of bone development (Image 7).
Image 7.
Haematoxylin/eosin of the bone defect at three months. G+M group. Various particles of biomaterial and chronic inflammatory component (A). G+M+S group. Osteogenesis. Black arrows: osteoblasts included in the osteoid (B).
Image 7.
Haematoxylin/eosin of the bone defect at three months. G+M group. Various particles of biomaterial and chronic inflammatory component (A). G+M+S group. Osteogenesis. Black arrows: osteoblasts included in the osteoid (B).
At six months, the treated G+M group continued to have particles of biomaterial encapsulated in a fibrotic matrix. The only signs of new bone were on the edges of the lesion, with a chronic inflammatory component. By contrast, the group treated with DP-MSC records full bone regeneration with well-organized, vascularized bone with a lamellar structure surrounded by haversian canals. The edges of the lesion are not distinguishable (Image 8).
Image 8.
Haematoxylin/eosin of the bone defect at six months. G+M group. Chronic inflammation (A). G+M+S group. Osteogenesis (B).
Image 8.
Haematoxylin/eosin of the bone defect at six months. G+M group. Chronic inflammation (A). G+M+S group. Osteogenesis (B).
2.7. Molecular Studies
At three and six months, an assessment was made of the TGF-β growth factor and the levels of type 1 procollagen, the serum concentrations of the calciotropic hormone calcitonin and the parathyroid hormone (PTH), and the endoglin expression.
As the neovascularization marker, we studied the endoglin crucial for TGF-β signaling in the endothelial cells. The results for the G+M+S group at three months reveal a significant increase in both endoglin expression and the production of growth factors such as TGFβ1. No changes were found in either the serum levels of endoglin or TGFβ1 in any of the groups at six months (
Figure 2 and
Figure 3).
The serum levels of procollagen l were higher in the treated groups at three months compared to the Sham and Control groups, with no differences between them. The G+M group at six months continued to record higher serum levels than the Sham and Control groups, while in the G+M+S groups these levels were similar to those in the Sham group (
Figure 4).
Following the lesion, we observed higher levels of PTH in all the groups compared to the Sham group, except for the Control group at three months and the G+M+S group at 6 months. At the same time, the treated groups maintained high serum levels compared to the Control group, with the exception of the G+M+S group at six months (
Figure 5).
No changes were found in the serum levels of calcitonin in the Sham, Control, and G+M groups. The G+M+S group, however, recorded a significant decrease at three months, attaining the Sham group’s levels after six months (
Figure 6).
3. Discussion
Our study confirms once more the accessibility of dental pulp to efficiently extract DP-MSC. The experimental model used is ideal for studying the efficacy of DP-MSC in osteogenesis, adjudging it to be financially viable with a high performance.
Although the characterization for human derived-MSCs has a well-established criterion, rat MSCs still lack the same consistency, mostly regarding surface markers. Some markers, such as CD29 and CD90 were more frequent, while others (e.g., CD73, CD105, CD44, and CD49e) were less so [
22,
23,
24,
25,
26,
27]. Isolates were evaluated here for the presence of the CD29, CD90 and CD49e expression and the absence of the expression of lineage markers, CD11b/c and CD45. The presence of the MSC population was confirmed in all the samples tested by the Lineage-CD29+CD90+ profile. Interestingly, the samples did not express the CD49e marker. Known as α5 integrins, this molecule associates non-covalently with CD29, a β1 integrin, to form the fibronectin receptor [
28]. The relatively low percentage of the CD49e+ population, ranging from 1-30%, might be explained by the sensitivity of the integrins α5β1 to trypsin [
29,
30]. The use of this proteolytic enzyme could therefore alter the antibody binding site, diminishing its recognition by the clone of the antibody used. Nevertheless, at the same time, the expression of CD90 and CD29 and the absence of lineage markers identify the MSC population.
The results from the radiological and histological studies confirm that there is no evidence of osteogenic activity in either Group C or in the G+M group, confirming that these conditions require enough cell material and a high osteogenic capability. A major characteristic of DP-MSC that makes them a promising osteogenic candidate is their co-differentiation in osteoblasts and endotheliocytes, fully integrating the blood vessels within the ossification lines, prompting the formation of vascularized bone tissue [
31,
32]. A clinical trial conducted in 2019 [
33] has found that these cells produce extracellular vesicles (EVs) that contain pro-angiogenic factors (e.g., VEGF, angiopoietin, or angiogenin). The Gen-Os® biomaterial’s angiogenic potential is also greater than other synthetic materials, as it favors the secretion of VEGF [
34].
DP-MSC are therefore the mainstay of the osteogenic process ending at six months, with full bone formation. DP-MSC contain immunomodulating and anti-inflammatory factors such as IL-6, IL-8, transforming growth factor-beta (TGF-β), hepatocyte growth factor (HGF), and indoleamine 2,3-dioxygenase (IDO) [
35,
36], all of which play an important part in regulating the equilibrium between the anti- and pro-inflammatory properties of cytokines [
36,
37,
38,
39], suppressing the activation of T cells, the proliferation of peripheral blood mononuclear cells, and even allogeneic immune responses [
35,
38,
39,
40].
In this study, the endoglin in the group treated with DP-MSC at three months, and like TGF-β, recorded high levels of expression, with this protein being part of the receptor complex of TGF-β that is highly expressed in several kinds of proliferating cells and stem cells in differentiation during the initial stages of wound healing. Angiogenesis is a major step in bone healing, and endoglin is highly expressed in proliferating endothelial cells and pericytes during angiogenesis. This means it is not only a proliferation marker, as it also regulates cell functions such as adhesion, migration, and cell permeability. Endoglin has a significant role to play in the regulation of the movement of stem cells from a specific niche to the damaged region [
41]. Besides high levels of TGF-β1, our results have revealed an increase in collagen I content in the G+M+S group at three months, confirming that the TGF-β1 factor is involved in cell differentiation and collagen synthesis [
42]. As type I collagen is the main component of the organic matrix, this increase in our results will render the bone tissue flexible and resistant to impact loads [
43].
TGF-β regulates the proliferation and differentiation of chondrocytes and osteoblasts [
44,
45], although mainly, its osteoinductive properties render it a major growth factor, performing a key role in the different stages of the bone healing process [
46] by regulating the production of the extracellular matrix (ECM), proteases, protease inhibitors, migration, chemotaxis, and the proliferation of different kinds of cells, including stem cells, which regulate angiogenesis and the formation of granulation tissue, remodeling of the matrix, and scar maturation.
The implementation of DP-MSC on the biomaterials also increases the serum levels of PTH at three months, favoring bone regeneration, as it has osteoclastic and osteogenic activity, whereby the bone formation process promoted by the PTH is very similar to the same process in physiological conditions. Bone repair is not a unique osteogenic process and requires angiogenesis and the participation of osteoclasts [
47]. Bone repair without osteoclastic activity may lead to excessive bone development [
48]. PTH has been shown to regulate the differentiation pathway of bone marrow mesenchymal stem/stromal cells (BMSC), increase their osteogenic differentiation, and inhibit their adipogenic differentiation [
49,
50,
51,
52,
53].
The study has also found a reduction in the serum concentration of CT in the group treated with DP-MSC. The changes observed in the levels of calciotropic hormones in rats could increase bone resorption and promote mineral metabolism in the DPSC group [
54,
55].
The alloplastic matrices take a backseat as the mere facilitator of cell functions, although it is essential to use the most appropriate ones to provide cell support and facilitate regeneration with osteoinductive factors. Each bone substitute has different properties, and Gen-Os® (OsteoBiol®, Italy) appears to provide the best guarantees for osteogenesis. Indeed, an assay [
56] shows that the proliferation and recruitment of MSC is higher with Gen-Os® (OsteoBiol®, Italy) than with other alloplast grafts such as Bio-Oss®. It also produces a smaller inflammatory reaction than other synthetic biomaterials [
57].
It seems fair to say that the properties of the different elements used together favour osteogenesis, acting as an ideal bone substitute and enabling the bone to regenerate while upholding all its characteristics.
4. Materials and Methods
The sample consisted of 46 WISTAR Han™ (RccHan™:WIST) adult male rats (8 months - 350 g) supplied by Harlan Laboratories Models, S.L.
All the experiments were conducted following the guidelines established by Spain’s Royal Decree 53/013 (Approval by the Ethics Committee at Salamanca University, filed under number 361).
4.1. Experimental Design, Specimens, and Groups (Figure 7)
The molars of rats were used for the extraction of DP-MSC (two pulps per extraction), as described in due course (N = 6).
The 40 remaining rats were distributed into the following groups:
Sham (Sham): Subject to the same anesthetic/surgical technique as the other groups except for the bone lesion and its treatment, as indicated below. The specimens remained caged for three months, whereupon the samples were gathered and they were euthanized (N = 4).
Control (C): The bone lesion was effected but not treated, as indicated below. The specimens remained caged for three and six months (N = 6 in each case), whereupon the samples were gathered and they were euthanised (N = 12).
Gen-Os® + Evolution® (G+M): The bone lesion was effected and then treated with the substitute bone biomaterial and covered with resorbable membrane, as indicated below The specimens remained caged for three and six months (N = 6 in each case), whereupon the samples were gathered and they were euthanized (N = 12).
Gen-Os® + Evolution® + DPSC (G+M+S): The bone lesion was effected and then treated with the substitute bone biomaterial plus DP-MSC and covered with resorbable membrane, as indicated below. The specimens remained caged for three and six months (N = 6 in each case), whereupon the samples were gathered and they were euthanized (N = 12).
Figure 7.
Flowchart scheme of the study’s experimental design.
Figure 7.
Flowchart scheme of the study’s experimental design.
4.2. Biomaterials
The Evolution® (OsteoBiol®, Italy) resorbable membrane used involved dense collagen fibres that covered the induced bone defect.
The Gen-Os® (OsteoBiol®, Italy) 250-1000 µm granulated bone substitute used was of biological origin (pig) based on hydroxyapatite (HA) and type I collagen. Its composition is 80% spongy bone and 20% cortical with high osteoconductivity. This three-dimensional matrix provided the support for the MSCs from the dental pulp.
4.3. Obtaining the DP-MSC from Pulp Tissue
Rat molars were extracted under general intraperitoneal anesthetic, and the specimens were then euthanized with an anesthetic overdose.
The molars were then crushed with a gouge forceps, obtaining pulp tissue with a Hemingway® bone curette (Model 1145/0).
The samples were placed in a falcon tube (50 ml) with 30 ml of DMEM (Dulbecco´s Modified Eagle Medium®, Invitrogen), supplemented with penicillin (100 U/ml), streptomycin (100 μg/ml), and fetal bovine serum (FBS) at 10% (Invitrogen), and cold stored for a maximum of six hours.
For enzymatic disaggregation, the tissue was incubated for 70 minutes at 37ºC and shaken with collagenase I (at 0.2%) (Worthington Biochemicals Corporation®) and dispase II (at 0.4%) (Sanko Junyaku Co. Ltd, Tokyo, Japan). It was centrifuged for 10 minutes at 1200 rpm, and the sediment was spread on 9 cm2 plates and incubated at 37°C.
Cell culture and growth was performed by incubating the cells at 37ºC in a humid environment with 5% CO2 in a DMEM (Dulbecco's Modified Eagle Medium®), supplemented with FBS (10%) (Equitech-Bio Inc., Kerrville, TX) and with penicillin (100 U/ml) and streptomycin (100 μg/ml). The culture medium was replaced twice per week, and any unattached cells were discarded.
The culture was maintained until a cell confluence of at least 80-90% was reached, whereupon the culture medium was eliminated, washed with sterile phosphate buffered saline - PBS (Invitrogen Corporation, Paisely, UK®) and incubated for five minutes at 37ºC with trypsin 1 (0.05%) (Invitrogen Corporation, Paisely, UK®).
Trypsin 1 was neutralized before the addition of a culture medium for MSCs. Cellularity was increased by again seeding a concentration of 5,000 cells/cm2 in larger beakers, up to step 10 (S10).
Cells were isolated for their immunophenotypic characterisation and induction to differentiation in several steps (1, 3, 6, and 10) to verify their osteoblastic capacity. Two parallel samples were cultivated under the same conditions (at 37ºC, 5% CO2, and > 90% humidity): one for osteoblast differentiation and the other as control.
The medium for osteogenic differentiation (NH OsteoDiff Medium, Miltenyi Biotec, Germany) contained β-glycerol-phosphate, ascorbic acid-2-phosphate, dexamethasone, and FBS. The cultures were maintained for 10 days and the medium was changed every 3-4 days. They were then washed with PBS cold 70% ethanol fraction (Merck KGaA, Darmstadt, Germany) and fixed for 10 minutes.
Osteoblast differentiation was confirmed by observing alkaline phosphatase activity by staining for 20-30 minutes with NBT/BCIP (Nitroblue tetrazolium chloride 5-bromo-4-chloro-3-indolyl-phosphate, Roche, Basel, Switzerland), and contrasting for two minutes with hematoxylin (1 mL) (Merck KGaA, Darmstadt, Germany).
4.4. Surgical Procedure
All the surgical procedures were performed under strictly sterile conditions in a laminar flow cabinet.
Following intraperitoneal anesthetization with a mixture of xylazine, ketamine, and saline solution (ratio 2:3:3, at a dose of 0.2 mL/100 g body weight), the surgical area was shaved and an antiseptic solution of povidone-iodine (Betadine, Asta Medica®) was applied.
A 15 mm longitudinal incision was made 2 mm above the lower edge of the mandibular body, providing submandibular entry, accessing the ascendant branch and angle of the right jaw, where a 5 mm diameter circular bone lesion was made with an electric micromotor and 5 mm trephine drill, and considered to be a critical bone defect (Image 9). The bone defect was subsequently covered as described for each group and the incisions were closed in layers. The specimens were caged in a state of good health until three or six months had elapsed after surgery, with direct ad libitum access to food and water.
Image 9.
Access to the insertion crest. Circular osteotomy with trephine drill.
Image 9.
Access to the insertion crest. Circular osteotomy with trephine drill.
The post-surgical clinical monitoring involved assessing the specimens’ overall condition, the appearance of the wound and surgical area, bleeding or discharge from the wound, rejection of the biomaterials, and any degenerative changes due to the dental lesion.
Once the monitoring period had ended, the next step involved gathering the samples and euthanising the specimens:
Total blood by aortic puncture, centrifuged (20 minutes at 4500 rpm and 4ºC), extraction of plasma, aliquoted, and frozen at -80ºC.
Perilesional tissue (bone and muscle) placed in liquid nitrogen and stored at -80ºC.
Affected mandibulae removed in bloc, removal of the tissue on the bone surface, with some being immersed in formaldehyde 3.7-4.0% w/v buffered to pH = 7, and others stored at -80°C.
4.5. Verification of Cell Viability
The DP-MSC inserted in the biomaterial were viable, with a small volume of the implanted cells being retained and returned to the laboratory for re-seeding.
4.6. Macroscopic Study
Prior to euthanizing, a descriptive assessment was made of the following: anatomic and tissue formation; signs of infection; existence of fractures; presence or disappearance of biomaterials; degenerative dental transformations; bone sequestrations; surface morphology, and solidity of the defect.
4.7. Scanning with Micro-Computed Tomography (micro-CT)
Once the period of monitoring had elapsed, and prior to euthanizing, two specimens for each period of evolution and group were subject to an x-ray evaluation through micro-computed tomography (mCT) with SuperArgus equipment from SEDECAL Medical Systems.
The evaluation addressed the regeneration of the defect in the four groups and the different timeframes, together with the changes in the morphology of the biomaterial following its implanting.
It also considered the biomaterial’s loss of homogeneity over the study period as a process of integration in the host bone and remodeling.
Bone trabeculae were also detected within the biomaterial.
The biomaterial’s dispersion within the medullary cavity was studied, together with the appearance on the periphery of the implanted area of bone trabeculae forming “bridges” between the biomaterial and the adjacent bone cortical.
4.8. Histomorphometric Evaluation
The hemimandibles were obtained and fixed in 4% formalin for 48 hours before being immersed in 10% EDTA for 30 days. The bone tissues were then embedded and the paraffin blocks were cut into 4 μm sections. The sections from each sample were stained with hematoxylin-eosin (H&E) for their comparative histological study with optical microscopy (Olympus, Tokyo, Japan).
4.9. ELISA
This technique was used to determine the concentrations of calciotropic hormones such as calcitonin (CT) and parathormone (PTH), the growth factor TGF-β, the levels of type I procollagen, and the endoglin expression, in the serum of the rats in each group. This involved commercial kits, and the procedures followed the manufacturer’s instructions. Parathyroid Hormone (PTH): Enzyme-linked Immunosorbent Assay kit ABK1-E1599. Calcitonin: Enzyme-linked Immunosorbent Assay kit ABK1-E1354. Collagen Type 1 Alpha 2(COL1α2): Enzyme-linked Immunosorbent Assay kit ABK1-744. Transforming Growth Factor Beta 1 (TGFβ1): Enzyme-linked Immunosorbent Assay kit ABK1-E1239. Endoglin (ENG): Enzyme-linked Immunosorbent Assay kit ABK1-E2467.
4.10. Statistical Study
The quantitative data have been represented as X±SD (mean ± standard deviation). A value of p < 0.05 was accepted as a significant result. The statistical software used was NCSS 2007 and Gess 2006 - Version: 07.1.21 - Released 1 June 2011 (Dr Jerry L. Hintze, Kaysville, Utah, USA).
Author Contributions
Conceptualization, F.J.G.C., H.H.A., and M.B.G.C.; methodology, F.J.G.C., H.H.A., J.B.S., and M.B.G.C.; validation, F.J.G.C., H.H.A., and M.B.G.C.; formal analysis, F.J.G.C.; investigation, H.H.A., M.B.G.C., F.P.M. and S.M.O.; resources, J.B.S., M.B.G.C. and S.M.O.; writing—original draft preparation, F.J.G.C., F.S.G., S.M.O., and M.B.G.C.; writing—review and editing, F.J.G.C., F.S.G., S.M.O., and M.B.G.C.; visualization, H.H.A. S.M.O., and M.B.G.C.; supervision, F.J.G.C., and F.S.G.; project administration, F.J.G.C., and MBGC; funding acquisition, F.J.G.C., and H.H.A. All authors have read and agreed to the published version of the manuscript.