Introduction

For several decades, cell therapies and tissue engineering have been increasingly considered as new therapeutic strategies to treat large bone defects or non-union fractions.[13] Recreating the complexity of native tissue through three-dimensional (3D) constructs, cell-material biografts, or biomaterials in general remains incredibly challenging.[4] Importantly, cell adhesion, proliferation, differentiation, and subsequenttissue formation, as well as clinical outcome, are strongly dependent on the choice of scaffold design, the ability of cells to interact with the material, and the successful creation of an artificial microenvironment. Ideally scaffolds should recreate the architectural, chemical, mechanical, and electrical properties of the host tissue.[5] A variety of tissues, including bone, are susceptible to electrical stimulation, and the inherent piezoelectric properties of bone tissue has led to the hypothesis that electroactive materials could be suitable for promoting in vivo bone repair and osteogenesis.[6,7] To this end, conjugated organic polymers such as polypyrrole (PPy), polyaniline (PANi), or poly(3,4-ethylenedioxythiophene) (PEDOT) have been increasingly considered for the development of novel electroactive scaffolds due to their reported cytocompatibility.[810] After implantation of such materials in vivo, cell adhesion to a conjugated polymer membrane will be in part supported and maintained by the extracellular matrix (ECM) deposited by cells in the vicinity, and proteins present in the cell media serum. The initial adsorbance of proteins to a material is therefore critical in determining and controlling cell adhesion, immune response, and overall performance of an implant, and is therefore a key consideration during biomaterial design.[11] This is particularly the case when using conjugated polymer-based materials, where the oxidation state of the polymer strongly determines material properties and is responsive to the local environment. Upon doping of a conjugated polymer, a change in electronic configuration or redox state not only leads to the conduction of electrical current, but furthermore alters intermolecular interactions with proteins and significantly influences protein adsorption.[12,13] For example, fibronectin (Fn), a common protein of the ECM which promotes cellular adhesion, has been demonstrated to interact with cells in a protein conformation dependent manner.[14] As such, the redox state-dependent conformation and adhesion of Fn to a conjugated substrate is key in determining interactions at the biological interface.[15] Hence, by measuring the protein adsorption profile we can begin to elucidate the effect of conjugated polymers on the promotion of cellular adhesion. The ability to tune redox states by choice of dopant or synthesis strategy thereby renders conjugated polymers promising candidates for studying controlled protein adsorption and the resultant cell response.

Despite these promising properties of conjugated polymers, their application in tissue engineering remains hindered by many severe drawbacks associated with their use. For example, their low solubility in both aqueous and organic solvents, reduced processability, high heterogeneity when synthesized via traditional routes, lack of cell adhesion sites, and nonbiodegradability are severely limiting. As a result, recent research has focused on the design of biodegradable conductive scaffolds by restricting conjugation to oligomers of 4 to 8 repeating units.[16,17] Oligomers provide a versatile route for customized, end-functionalized units with tuneable solubility and processability[18] that allow for significantly enhanced freedom and flexibility with respect to scaffold design and integration into higher order architectures and biocomposites, while maintaining a significant degree of electroactivity. For example, a tetrameric oligomer of aniline exhibits redox and doping chemistries similar to poly(aniline) and has been shown to present even higher stability during charging-discharging cycling processes due to its shorter chain length and lack of long-chain swelling.[19] Further, oligomers can be synthesized at defined molecular weights, leading to monodisperse products that will result in more concise and reproducible material design.[18] Substantial efforts have gone into the development of aniline tetramer[2022] or aniline pentamer[7] based scaffolds, with conjugated units being grafted to gelatine, poly(lactic acid) (PLA), poly(ethylene glycol), or poly(caprolactone) (PCL).

Preliminary investigations into the cytocompatibility, biodegradation or in vivo response of such materials have been reported, but open questions with respect to the application of such scaffolds in bone tissue engineering still remain.[2325] The aim of the current study therefore was to produce electrospun membranes based on aniline tetramer and evaluate these membranes in vitro. Our study in particular addressed the question of long-term studies and whether these biodegradable, electrically active membranes can successfully be used as substrates for bone tissue engineering. Our study further underpins the challenges encountered in the field and emphasizes on the need of in depth characterization and defined conditions to meet the demand for highly complex in vitro model systems.

Material and methods

Unless stated otherwise, all chemicals were purchased from Sigma-Aldrich, UK and used without further purification. Cell culture consumables and supplements were from Life Technologies, Thermo Fisher, UK, sterile, and used as provided.

Aniline tetramer-co-poly(caprolactone) synthesis and characterization

Aniline tetramer-co-poly(caprolactone) (AT–PCL) in the emeraldine base (EB) state was prepared in a two-step synthesis, as previously reported[22] and detailed in supporting information (SI). AT-PCL was obtained at a molecular weight of ~1300 g/mol [estimated based on nuclear magnetic resonance (NMR) spectra as described in the SI]. The reaction scheme and NMR spectra are displayed in Supplementary Figs. S1 and S2. Optoelectronic activity was assessed by absorption spectroscopy. AT–PCL was dissolved in dimethyl sulfoxide (DMSO) and doped with phytic acid at increasing vol% (0.005–0.1 M). Spectra were recorded in the wavelength range of 200–1000 nm at a step size of 10 nm (Mithras2 Plate reader, Berthold Technologies, Germany).

Electrospinning

AT-PCL was blended with high molecular weight PCL (Mw ~ 70,000–90,000 Da) at w/w ratios of 25%, 50%, or 75% AT–PCL to PCL, and dissolved in a mixture of CHCl3 :CH3OH (9 : 1 vol%). Solutions were doped with phytic acid to a final concentration of 0.02 M. Polymer solutions were processed into fibrous non-wovens on a custom built electrospinning device, as detailed in SI.[26] All scaffolds were spun at 10/ −5 kV, 20 μL/min flow rate, and a needle to collector distance of 15 cm. Fiber diameters were measured with ImageJ software, based on scanning electron micrographs recorded at three random spots on the membrane (nine images, 100 fibers per condition). Electrical characterization of electrospun membranes was carried out using a Keithley 2400 source meter, controlled by customized LabView software. The sheet resistance of dry, fibrous mats, doped in 0.1 M hydrochloric acid (HCl) was measured by standard four-point probe techniques using a custom made probe setup.

Protein–membrane interaction

Atomic force microscopy

Protein-binding force was quantified by atomic force microscopy (AFM). Fn conjugated AFM probes were prepared and characterized as previously described.[15] Briefly, the probes, (PNP-DB, Nanoworld, Switzerland) with a nominal spring constant of 0.06 N/m, were O2 plasma treated for 3 min, and then placed into 1 vol% (ethoxydimethylsilyl)propylamine in toluene for 1 h at room temperature. The probes were rinsed with toluene, followed by PBS, and immersed in 2.5 vol% glutaraldehyde in PBS for 1 h at room temperature. The probes were then rinsed in PBS and immersed in a 10 mg/mL Fn solution from human plasma for 1 h. The probes were finally rinsed and refrigerated in PBS solution until use. AFM force spectroscopy of Fn adhesion was performed using a Keysight Technologies 5500 AFM (USA). The force spectroscopy measurements were performed across a 40 × 40 μm2 area with 64 individual force curves, on three separate locations across the sample. The measurements were repeated with three Fn functionalized probes. Force spectroscopy results were analyzed using PUNIAS software (punias.free.fr). Specifically, the protein unfolding length is measured as the length between two consecutive peaks within the AFM force curve, where the peaks represent a detachment event between the protein and the surface; the total protein unfolding length is calculated as the sum of all the unfolding lengths within each force curve minus the unfolding length between the first two peaks.

Protein adsorption on the membranes

Protein adsorption was further quantified through the use of fluorescein conjugated fibronectin (FITC-Fn). Scaffolds of 6 mm in diameter were placed in black flat bottom 96-well plates, and incubated with 100 μL FITC-Fn solution (10 μg/mL in PBS) for 20 h. Samples were washed three times in PBS and fluorescence intensity measured at λex = 475 nm and λem = 512 nm (Perkin Elmer EnSpire Plate Reader, USA).

In vitro evaluation

Electrospun membranes were prepared and fixed in cell culture plates as previously reported.[27] In brief, scaffolds of 6 mm diameter were cut out with biopsy punches, and fixed with insect pins on cured poly(dimethyl siloxane) (PDMS). Membranes were sterilized under UV overnight. Mouse pre-osteogenic precursor cells (MC3T3-E1, ATCC, UK) were expanded in minimum essential medium Eagle - alpha modification (α-MEM) supplemented with 10 vol% FBS and 1 vol% penicillin/streptomycin. MC3T3-E1 were seeded at a density of 110,000 cells per scaffold and cultured for 3 days in basal media. Osteogenic differentiation was initiated by addition of α-MEM supplemented with 10vol% HyClone FBS (Thermo Scientific, UK), 1 vol% penicillin/streptomycin, 5 μg/mL L-ascorbic acid and 2 mM β-glycerophosphate disodium salt. Media was changed every 2 days over a differentiation period of 21 days. In vitro evaluation was accomplished with n = 3 individual experiments and three individual scaffolds per experiment per condition.

Cell proliferation was assessed based on DNA quantification. Briefly, constructs were washed in 1x PBS, cells lysed in distilled water, and three repeated freeze thawing cycles. DNA was quantified with Hoechst DNA-binding agent. Fluorescence was measured at λex = 360 nm and λem = 460 nm (Perkin Elmer EnSpire Plate Reader, USA). Osteogenic differentiation was assessed after 21 days in culture. For qPCR, scaffolds were lysed in RTL buffer, supplemented with 1 vol% β-mercaptoethanol. RNA was isolated with an RNA isolation kit (RNeasy Mini Kit, Qiagen, Germany) and cDNA synthesis was accomplished with a reverse transcriptase kit (Reverse Transcriptase Kit, Qiagen, Hilden, Germany). Samples were analyzed for the expression of ALPL, COL1A1, and RUNX2, normalized to the housekeeping gene 18S (all primers from Life Technologies, Thermo Fisher, UK). Sequences of forward and reverse primer are found in Supplementary Table S1. qPCR was run on a QuantStudio6 real-time PCR System (Applied Biosystems, USA). Samples were run with Express SYBR green super mix. Fold expression was calculated as 2−ΔΔCt normalized to 18S and gene expression on day 0. Immunohistochemistry was accomplished on formaldehyde-fixed samples. Samples were incubated with anti-osteocalcin antibody (rabbit anti mouse, aa 59 74, MERCK, Millipore, Germany) followed by a fluorescently labeled secondary antibody (anti-rabbit IgG, Alexa 555, Life Technologies, Thermo Fisher, UK). Samples were mounted with glycergel mounting media (Dako, Agilent technologies, Ely, UK). Images were acquired on a confocal microscope (LSM780, Carl Zeiss AG, Switzerland).

Statistical analysis

Significant differences were assessed by one way analysis of variance, followed by pairwise comparison with Bonferroni correction. Results were accepted as significantly different for P <0.05.

Results

Synthesis and electrical characterization of AT-PCL and AT-PCL/PCL membranes

Followingapreviouslypublished synthetic route(Supplementary Fig. S1A), AD-PCL with a molecular weight of -1000 g/mol (7 caprolactone units) was synthesized.[22] Subsequently, a second AD unit was oxidatively coupled to AD-PCL, yielding AT-PCL with a molecular weight of ~1300 g/mol. Due to the acidic reaction conditions, AT PCL was obtained in its doped, emeraldine salt (ES) state and was de-doped in 0.1 M NH4OH, yielding a dark blue/purple product in its emeraldine base (EB) state (Supplementary Fig. S1).

To investigate electrical conductivity, four-point probe measurements were accomplished on HCl-doped electrospun membranes (detailed below). At increasing ratios ofAT PCL to PCL from 25% to 75%, a reduced sheet resistance was observed (24.7 ± 5.6 × 106, 7.1±1.7 × 106, and 3.0 ± 0.6 × 106 Ω/sq, respectively), indicating increased electrical conductivity [Fig. 1(a)]. However, the instability ofHCl doping led to immediate de-doping when immersed in physiological buffer (PBS). This prevents its useas dopantforcellculturemembranes, where stable doping is required. Phytic acid has been reported as a stable dopant for poly(aniline),[28] and we therefore hypothesized that this non-toxic acid would work comparably in the herewith designed system. The optoelectronic activityofAT PCL doped with phytic acid was confirmed by UV–VIS absorbance spectrometry (Fig. 1). The insulating EB state has a characteristic absorption maximum at ~600 nm, attributed to the excitonic transition πBπQ of benzenoid (B) to quinoid (Q). Upon the addition of phytic acid AT PCL was doped to the ES state, resulting in the emergence of an additional adsorption maximum at ~400 nm, characteristic of bipolarons. Interestingly, however, these changes in oxidation and doping were not reflected in four-point probe electrical measurements. Electrospun membranes, as described below, of AT PCL blended with PCL and doped with 0.02 M phytic acid did not conduct any measurable current, and sheet resistance could not be determined at the applied voltage, attributed to the weak acidic properties of phytic acid and its minor doping capabilities compared with stronger acids such as HCl.[28]

FIG. 1.
figure 1

Electrical and optoelectronic characterization of AT–PCL. (a) Sheet resistance of electrospun membranes, produced at different AT–PCL to PCL ratios and doped with 0.1 M HCl. (b) Spectra of AT–PCL in DMSO doped with phytic acid at different molarities.

Aniline tetramer has previously been shown to exhibit similar redox properties to PANi.[19] As a co-polymer with PCL these properties were maintained and AT PCL was shown to retain the characteristic properties of the aniline oligomer. At phytic acid concentrations as low as 0.005 M, we could demonstrate the optoelectronic activity of AT PCL. This warranted further in vitro evaluation and the study of protein adsorption on membranes as a function of AT–PCL/PCL ratios and doping state.

Electrospinning

Fibrous scaffolds prepared by electrospinning have proven particularly versatile for the design of biografts due to their close architectural resemblance to the fibrous ECM, and they have therefore been widely applied in tissue engineering. However, this method cannot typically be applied to fully conjugated systems, due to their poor solubility in organic solvents. Conjugated oligomers, and in particular end functionalized oligomers or block-co-polymers offer opportunities to address this problem, whereby the solubility of an insulating partner can be utilized to modify the system properties. However, electrospinning of aniline tetramer or pentamer based materials for bone tissue engineering has been restricted to very few publications, e.g. Liu et al.[7] in which the conductive unit was grafted to a gelatin-PLLA polymer. As an alternative to chemical grafting, blending polymer systems have emerged as a versatile method to combine distinct material properties during the process of electrospinning. To this end, AT–PCL was blended with high molecular weight PCL (70,000–90,000 Da) and doped by adding phytic acid to the polymer solution prior to electrospinning (final concentration of 0.02 M in CHCl3 :CH3OH). Doped solutions turned green immediately, indicative of the ES state of AT–PCL, whereas non-doped solutions retained the dark purple color characteristic of the EB state. Electrospinning of AT–PCL/PCL blends in a 9 : 1 mixture of CHCl3 :CH3OH resulted in the deposition of fibrous, non-woven, free standing mats, of blue or green color depending on the doping state. Fiber diameters were measured on scanning electron microscopy (SEM) images and the mean fiber diameter was calculated based on 100 individual measurements. For all blending ratios, these were in the sub-micrometer range. For pristine membranes, mean fiber diameters were 0.4 ±0.1 μm (25% AT–PCL); 0.3 ±0.1 μm (50% AT–PCL); 0.3 ±0.1 μm (75% AT–PCL), while doped membranes displayed values of 0.3 ±0.1 μm (25% AT–PCL); 0.4 ±0.1 μm (50% AT–PCL); 0.3 ± 0.1 μm (75% AT–PCL) (Fig. 2). Electrospun membranes retained their color when immersed in PBS or cell culture media for several hours, indicating stable doping, which significantly surpassed that of HCl-doped membranes. AT–PCL/PCL ratios of 50% or 75% resulted in the occasional formation of beaded membranes, due to lower solution viscosity. In the doped state, however, increased electrical conductivity resulted in the deposition of homogenous fibers.

FIG. 2
figure 2

SEM images of fibrous membranes. AT–PCL was blended at different w/w ratios (25%, 50%, and 75%) with high molecular weight PCLand processed by electrospinning. For all AT–PCL to PCL ratios, free standing membranes were produced. (a) Pristine AT–PCL membranes, (b) phytic acid-doped AT–PCL membranes.

After in vivo implantation or in vitro cell culture, a protein layer is adsorbed on the material within the first minutes or hours, then regulates further downstream signals for cell adhesion, cell proliferation, and/or differentiation.[11] We propose that this stable doping for several hours is sufficient to significantly alter protein adsorption and concomitant cell response, and therefore represents an important strategy for the design of electrically active 3D scaffolds for tissue engineering.

Protein-membrane interaction

Controlling protein adsorption via material properties is a powerful means to tune bio-interfaces and steer material-cell interactions. We were therefore interested in protein adsorption on doped versus pristine substrates, following previously published reports that the oxidative state of a conjugated polymer may induce significant changes in both protein adsorption levels and conformation.[12,13,15] As a major constituent of native ECM, we investigated Fn-membrane interactions with atomic force spectroscopy. The required force and work to pull off the protein from the membrane, and the protein’s length at release are displayed in Fig. 3. On doped samples of 50% or 75% AT–PCL, significantly higher forces were required to pull Fn when compared with pristine membranes at the same concentration. Similarly, significantly longer Fn length was observed on pristine 50% AT–PCL compared with doped 50% AT–PCL. Protein length was also shown to be significantly longer on membranes containing 25% AT–PCL compared with 50% or 75%. These results indicate reduced protein-membrane surface interactions on doped substrates, implying a change in protein conformation on the surfaces of doped substrates compared with pristine ones. The observed protein lengths are longer than a single Fn protein (160 nm), which can be attributed to the 3D nature of the electrospun polymer fibers. As the entire AFM probe has been functionalised with Fn, multiple proteins on the apex and side of the probe may bind to the fibers, contributing to longer unbinding lengths.

FIG. 3
figure 3

(a) Fibronectin-membrane interactions recorded with atomic force microscopy (AFM). Force recorded to remove Fn from the substrate, calculated length of the folded/unfolded protein, and work required to remove Fn from the substrate. Results are presented as mean ± SE with * P < 0.05; #P < 0.05 compared with 50% AT–PCL; @P< 0.05 compared with 75% AT–PCL. (b) Schematic representation of the results obtained in (a). (i) Pristine samples: an increased length is indicative for an extended Fn conformation, with larger contact area with the substrate, whereas in (ii) (doped samples), a shorter length indicates a more folded, compact Fn conformation, with reduced contact area. (c) Quantification of FITC-Fn adsorption on the respective substrates. Results are presented as mean ± SD.

We therefore hypothesize that the electrochemical properties of the substrate have a significant effect on protein conformation, with Fn being in a more relaxed, extended conformation on pristine samples, exposing a larger portion of protein binding sites to the membrane, resulting in longer binding length.[12,13] Reduction of the conjugated system to the ES state results in different protein adhesion, as a result of reduced protein-membrane contact (length) induced by a more coil-like folding of Fn. The amount of adsorbed Fn was, however, similar on all substrates, as indicated by the fluorescence intensity of bound FITC-Fn (Fig. 3). We therefore conclude that while the oxidative state of conjugated membranes does not control absolute protein adsorption concentration, it does have a significant effect on conformation. This change in protein conformation could have a significant effect on cell adhesion and further downstream events.[14]

In vitro evaluation

As a model system for osteogenesis in vitro, we assessed osteogenic differentiation of MC3T3-E1 osteogenic precursor cells on the different membranes. Despite phytic acid doping not resulting in membranes of measurable conductivity, we hypothesized the influence of doping on protein conformation would induce a distinct response on MC3T3-E1 differentiation. Intercellular signaling or charge transport might furthermore be influenced by the material properties, warranting the evaluation of our membranes in vitro. DNA content, based on fluorescent spectroscopy after Hoechst staining, remained constant over a period of 21 days and no significant differences in DNA content were observed between different substrates (Supplementary Fig. S3). Our results therefore demonstrate that AT–PCL membranes are suitable cell culture substrates for long-term studies of up to 21 days under differentiation conditions, without any major decrease in cell number. This suggests that even at higher relative concentrations of 75% AT–PCL, the conjugated oligomer does not have a cytotoxic effect. Furthermore, these long-term studies show that phytic acid can be used as a non-toxic dopant.

Gene expression analysis of key osteogenic markers was undertaken by qPCR (Fig. 4). ALPL, COL1A1, and RUNX2 expression was increased on day 21 compared with day 3 on all substrates. Of note, on doped 50% AT–PCL and on pristine 25% AT–PCL, significantly higher COL1A1 expression was found on day 21 compared with day 3. Additionally, RUNX2 expression on pristine 25% AT–PCL was significantly higher than any other condition, in particular when compared with doped 25% AT–PCL or pristine 75% AT–PCL. On all other samples, differences between days 21 and 3 were not statistically significant. Osteocalcin deposition, as assessed by immunohistochemistry and confocal microscopy (Fig. 5), supported the trend of increased osteogenic differentiation on 25% AT–PCL. Increased positive staining for osteocalcin was observed on pristine membranes, with enhanced positive staining on 25% AT–PCL compared with 50% or 75% AT–PCL.

FIG. 4
figure 4

Gene expression of osteogenic genes ALPL, COL1A1, and RUNX2. Results are presented as mean ± SE with *P< 0.05.

FIG. 5
figure 5

Confocal microscopy images of MC3T3-E1 cultured on electrospun membranes of AT–PCL for 21 days under osteogenic differentiation conditions. Cells were stained for osteocalcin (red), actin (green), and nuclei (blue). (a) Phytic acid-doped membranes, (b) pristine membranes.

The in vitro evaluation suggests that 25% AT–PCL membranes are most promising candidates for differentiation of MC3T3-E1 osteogenic precursor cells. The underlying mechanisms are still elusive, yet an altered protein conformation, as indicated by an increased Fn length upon removal (which is indicative for a more elongated, harmonica-like protein conformation), could have affected initial cell attachment, further ECM secretion of MC3T3-E1 and concomitant differentiation. Overall, these results indicate that osteogenic differentiation on membranes of higher AT–PCL content is not fully supported, and no osteo-inductive properties can be attributed to AT–PCL membranes. Limited reports on in vitro evaluation of scaffolds based on conjugated aniline oligomers lie in strong contrast with the vast amount of publications on synthesis routes, and our results can only marginally be related to previous work. While Liu et al. reported on the design of electrospun scaffolds containing aniline pentamer, with enhanced MC3T3-E1 proliferation compared with PLLA scaffolds alone, no qPCR was performed and it remains open to which extent conductive scaffolds induced enhanced gene expression of osteogenic markers.[7]

In summary, our AT PCL scaffolds did allow for cell attachment and supported culture of MC3T3-E1 over 21 days, allowing us to demonstrate cytocompatibility and long-term performance of AT PCL scaffolds in vitro. The cytotoxicity and biodegradation of aniline oligomers is still the subject of controversial discussion and the importance of functional endgroups has been emphasized to be influential in this regard.[17,29,30] Our study further underpins the necessity for sophisticated and highly reliable in vitro models for assessing the cytocompatibility and biodegradability of aniline oligomer-based scaffolds. The high variability observed among experimental repeats during this and prior research highlights the need for accurately controlled synthetic routes to conjugated materials, as a vital step towards the design of scaffolds that allow experimental reproducibility and consistency. Furthermore, analytical methods to assess optoelectronic properties and electrical conductivities need to be optimized to be applied to 3D porous culture substrates. In this respect, oligomeric materials, in contrast to polymers, are particularly attractive due to their potential to be synthesized under highly controlled conditions, resulting in standardized and reproducible chain length and chemical properties.[18] As such, we envisage that such areas of research will find increasing interest in the near future.

Conclusions

In summary, we created electropsun scaffolds of distinct AT PCL content and showed their cytocompatibility in MC3T3-E1 pre-osteogenic cell culture. In a blend system, we demonstrated that AT PCL is spinnable, allowing for the production of electrically active scaffolds of otherwise highly difficult to process conjugated polymers. In this respect, the use of oligomers, and blending with high molecular weight PCL was required to enable electrospinning. Depending on the content and redox state of the conjugated block, a distinct protein adsorption profile was observed, indicating conformational changes of Fn upon adhesion to doped or pristine (non-doped) samples. Our results further highlight the importance of highly controlled chemical synthesis, scaffold design and reliable in vitro models for studying cell material interactions and design 3D cell culture substrates.

Conflict of interest

The authors declare no competing financial interest.

Author contributions

A. G. G. designed the study, conducted experiments, performed data analysis and wrote the manuscript. C. D. S. supported study design, performed data analysis and revised the manuscript. A. G., A. A., and E. J. H. helped with experiments and data analysis. C. M. T. and S. E. H. provided valuable inputs to the study design and revised the manuscript. M. M. S. supervised the study, participated in scientific discussions, and revised the manuscript.