Collagen molecular organization preservation in human fascia lata and periosteum after tissue engineering

1 Introduction

Bone regeneration is a well-regulated process requiring the recruitment of skeletal stem/progenitor cells (SSPCs) from adjacent tissues. Among them periosteum is a connective tissue that surrounds the external surface of bones and has been widely described as physiologically essential for bone development and healing (Duchamp de Lageneste and Colnot, 2019; Julien et al., 2022; Perrin and Colnot, 2022). This membrane consists of a fibrous outer layer composed of type I collagen and a highly cellular inner layer composed of fibroblasts, osteoprogenitor cells and mesenchymal stem cells (Squier, Ghoneim, and Kremenak, 1990; Dwek, 2010). However, in case of critical bone defect after a massive trauma, a tumoral resection or an osteomyelitis, the surrounding soft tissues, such as periosteum, are usually damaged and the normal physiological healing capacity is insufficient (Vidal et al., 2020).

Current strategies for the reconstruction of large bone defects exploit surgical approaches such as Ilizarov’s (Aktuglu, Erol, and Vahabi, 2019) or Masquelet’s induced-membrane (Masquelet and Begue, 2010) techniques. However, these methods require lengthy and tedious processing. Bone repair therefore remains a major concern for orthopedic surgeons (Evans, Chang, and Knothe Tate, 2013; Delloye et al., 2014; Dalisson et al., 2021) which has led to the emergence of a series of alternative tissue engineering approaches in recent decades (Chen et al., 2015; Lou et al., 2021; Manon et al., 2023). Among these, Manon et al. recently showed the suitability of using human fascia lata matrices as a periosteal-like scaffold for bone regeneration, mainly through their collagen type I composition and appropriate mechanical properties (Manon et al., 2022; Manon et al., 2023). The human fascia lata (HFL) is a dense connective tissue made up of three layers of collagen fiber bundles with little cellularity (Pancheri et al., 2014; Blottner et al., 2019).

Along with bone marrow and muscle, the periosteum is an important source of skeletal stem/progenitor cells (SSPCs) with a high potential for bone regeneration (Zhang et al., 2005; Duchamp de Lageneste et al., 2018). These progenitors reside in a rich environment formed by the periosteal extracellular matrix, where cellular and molecular crosstalk is crucial to achieve bone regeneration (Julien et al., 2022). This three-dimensional scaffold therefore not only provides a structural and mechanical support (Costa-Almeida et al., 2015; Hussey, Dziki, and Badylak, 2018; Sainio and Järveläinen, 2020) but is also a modulator of tissue homeostasis (Theocharis et al., 2016). Consequently, any alteration of this structure can lead to the disruption of tissue function (Theocharis et al., 2016). Collagen constitutes the major ECM component (Shoulders and Raines, 2009; Yamauchi and Sricholpech, 2012) and plays a key role in maintaining biological and architectural integrity of tissues as well as providing physical support (Dong and Lv, 2016). Lin et al. (2018) further emphasized this point by revealing that cell-free collagen scaffold of decellularized periosteum alone efficiently promotes orthotopic remodeling of skeletal defects in vivo (Lin et al., 2018).

Fibril forming collagen type I, the most prevalent one, provides tensile strength to tissues (Theocharis et al., 2016) through its complex hierarchical structure. Encoded by the COL1A1 and COL1A2 genes, collagen precursors undergo series of post-translational modifications (PTMs) to form mature molecules (Gjaltema and Bank, 2017; Rappu et al., 2019), which are stabilized by two main types of intermolecular crosslinks important for their assembly and for collagen functions (Yamauchi and Sricholpech, 2012). First, enzymatic crosslinks, considered a post-translational modification of lysine residues, are mediated by the lysyl hydroxylases and lysyl oxidases through a highly controlled mechanism. Initially formed as immature divalent crosslinking (e.g., dehydrodihydroxylysinonorleucine), some are converted into mature trivalent crosslinks (e.g., pyridinoline) (Thorpe, 2010; Yamauchi and Sricholpech, 2012; Saito and Marumo, 2015; Kamml et al., 2023). Secondly, non-enzymatic crosslinks result from a glycation process forming advanced glycation end products (AGEs), among which pentosidine is the most frequently characterized (Paul and Bailey, 1996; Pritchard and Willett, 2017). All AGEs are formed by a Lysine-Arginine or Lysine-Lysine interaction and predominantly accumulate in collagen because of its low turnover. Their emergence leads to an increase in connective tissue rigidity and insolubility (Paul and Bailey, 1996; Gautieri et al., 2014). Finally, interfibrillar bridges also result from the interactions of non-collagenous components between glycosaminoglycans (GAG)-decorin (DCN) complexes of adjacent fibrils (Scott, 2003; Thorpe, 2010; Costa-Almeida et al., 2015). DCN is a member of class 1 small leucine-rich proteoglycans (SLRPs) essential for the regulation of various tissular processes (Halper, 2014; Gubbiotti et al., 2016) but most notably studied for its involvement in collagen fibrillogenesis and matrix stability. This protein also acts as a reservoir of growth factors and cytokines (Gubbiotti et al., 2016; Zhang et al., 2018).

In the present study harvested human periosteum (HP) and HFL samples were treated with several protocols of decellularization (Manon et al., 2023) as well as by irradiation. Whereas the general microstructure, composition, and biosafety of the ECM-derived scaffold had received considerable attention, a special analysis of collagen per se is required as the main component contributing to matrix-cell crosstalks and, hence, to the global integrity of the periosteum. Therefore, since decellularization and sterilization processes can have a negative impact on the ECM, the aim of the present study was to characterize the collagen scaffold obtained after these engineering processes. Here, we hypothesized and sought to confirm the global persistence of important ECM components of HP and HFL samples after chemical decellularization and sterilization, and therefore a persistence of their functional and structural properties. The ultimate goal of this study is to confirm the maintenance of HFL matrices properties in order to consider their future use as periosteal scaffold substitutes to promote bone regeneration.

2 Materials and methods

2.1 Sample collection and processing

HFL and HP were harvested from deceased donors (3 donors of average 86.81 ± 3.7 years) of the Human Anatomy Department of UCLouvain (Body donation service, IRB00008535, Brussels, Belgium) following the same procedure as previously described (Manon et al., 2022) and were subjected to five chemical decellularization (D) protocols (Table 1) performed under continuous agitation (N = 3 donors x 2 tissues x 5 protocols = 30 analyzes). The soap-based detergent protocols, D1 to D4, consisted of different combinations of Triton, Sodium dodecyl sulfate (SDS) and Deoxyribonuclease (DNase). Triton is a non-ionic detergent that disturbs the lipid-lipid and protein-DNA interactions. SDS is an ionic detergent that denatures protein structure and solubilizes nuclear and cellular membranes. DNase treatment is employed to cleave internucleotide bonds and eliminate undesirable cellular waste. In protocol D1, tissues were maintained for 12 h in SDS, followed by 4 h in Triton. The sequence was reversed for D2, with 12 h in Triton and by 4 h in SDS. For protocols D3 and D4, the sequence was the same as for D1 and D2, respectively, with an increase in the duration of the first treatment step maintained for 24 h instead of 12 h (SDS and Triton, respectively). Each protocol ended with a 12 h DNase treatment. D5 protocol is the one routinely used in the institutional tissue bank (van Steenberghe et al., 2017; Manon et al., 2023). In brief, it includes sequential immersion in pure acetone to degrease, ether to neutralize the previous compounds and 70° ethanol to further degrease. Following a washing step, sample are immersed in a hypertonic solution of NaCl (7%) and NaOH to finally been placed in H₂O₂ (7.5%). Quality controls were performed to confirm the full decellularization and to measure the pH and traces of H₂O₂ and acetone. To assess the efficacy of decellularization protocols, microscopic analyses were performed on DAPI and Hematoxylin and Eosin (HE) stained slices first. Then, 3 random samples of 25 mg were harvested from native and decellularized tissues of all protocols, and freeze-dried. The DNA content was measured using the Quant-iT PicoGreen DNA assay kit (ThermoFisher, Waltham, MA, United States) just after being extracted using the DNeasy^® Blood & Tissue kit (Qiagen, Hilden, Germany). The very precise steps are previously described (Manon et al., 2022; Manon et al., 2023). Samples were then stored at −80°C.

Table 1. Decellularization protocols (D), from D1 to D5.

Samples of decellularized HFL (D1-D5) were then frozen and sterilized by gamma(γ)-irradiation (HFLi, HFL irradiated) under dry ice using a minimal exposure dose of 25 kGy (IBA Mediris, Fleurus, Belgium), the recommended standard sterilization method for clinical applications (Dziedzic-Goclawska et al., 2005; Seto, Gatt, and Dunn, 2008; Harrell et al., 2018).

To perform the various analyses, samples were randomly taken from the entire decellularized scaffold to consider the whole tissue heterogeneity.

2.2 Polarized light microscopy

After paraffin histology processing, sirius red (SR) staining, an acidic hydrophilic stain, was used to highlight and evaluate qualitative modification in collagen organization (n_(HFL) = 4 samples/protocol, n_(HP) = 3 samples/protocol, with 3 sections per native samples and 1 per decellularization protocol). Deparaffinized and rehydrated sections were immersed in 1% Phosphomolybdic acid and then in 0.1% Sirius Red solution for 2 h. Sections were then incubated in hydrochloric acid before washing and mounting steps, then imaged under polarized light. Parallelly bound to collagen fibrils, this dye intensifies the birefringence of the samples which is manifested by shades of color under polarized light, each associated with properties of fiber organization. Thin fibers (diameter ≤0.8 µm) with a more scattered histoarchitecture are characterized by a weak birefringence and are detected in green/yellow while thick (diameter between 1.6–2.4 µm) organized ones appeared as shades of red and orange (Dayan et al., 1989). These parameters are quantified in order to highlight a possible alteration in collagen organization.

SR-stained sections were first imaged using the Axioplan microscope (ZEISS) with a 10-fold magnification. Linear polarization of light was obtained by using two filters that allow light to pass only in one direction of vibration: the polarizer and the analyzer, placed perpendicularly (Rittié, 2017). All slides were examined at the same time to ensure consistency of the experiment and were digitalized using the Nikon Digital sight DS SMC camera with NIS-Element BR 3.0 software. The average of five different acquisitions was assessed for each HFL section to cover the whole tissue surface, and three for HP, the surface of the sections being smaller. Collagen organization was then estimated with ImageJ software (Schneider, Rasband, and Eliceiri 2012) as a surface percentage of each color with a constructed macro (IREC imaging platform, 2IP). The analysis displayed green, yellow, orange, and red polarization colors.

To further validate our findings, whole slide images were then acquired with the ZEISS Axio Scan.Z1 Slide scanner under polarized light in order to approximate the proportion of each fiber type. These acquisitions were analyzed with QuPath (Bankhead et al., 2017) calibrated to detect the different shades of colors. With this approach, unlike with the previous one, only green, orange, and red colors are detected.

2.3 Raman spectroscopy

2.3.1 Spectra acquisition

Molecular composition of fresh frozen and processed HFL, HFLi and HP samples (n = 3/protocol for each type of samples) was investigated using Raman spectroscopy (Renishaw confocal Raman inVia Quontor). Spectra of each sample were obtained as the average of 3 consecutive spectra randomly distributed over the sample surface. Each spectrum was collected using a 785 nm monochromatic laser beam with a power of 50 mW, a 50-fold magnification and the Renishaw Centrus 21KP95 detector. Two types of spectra ranges were recorded.

Acquisitions were first performed with a static mode in a wavelength range of 700–1800 cm⁻¹ with an integration time of 10 s and 15 accumulations. The presence of non-enzymatic crosslinks (pentosidine, 1345 cm⁻¹) and hydroxyproline (872 cm⁻¹) was normalized to the amount of tissue using the proline peak (920 cm⁻¹). Enzymatic crosslink level is assessed by the ratio of mature to immature crosslinks (pyridinoline/dehydrodihydroxylysinonorleucine, 1670/1690 cm⁻¹). A 1657 cm⁻¹ supplementary band was added for the Amide-I peak-fit model of HP. This lipid vibration band was used to correct its potential interference on the measure of enzymatic crosslinks (Mandair et al., 2020).

To estimate the collagen-bound water content, a second wavelength range of 2700–3800 cm⁻¹ was recorded with extended mode, 10 s integration time, and 6 accumulations. Three parameters were considered: organic matrix-related water (3243 cm⁻¹), amine (NH) groups beside organic matrix-related water (3333 cm⁻¹) and hydroxide (OH) groups of hydroxyproline beside collagen-related water (3457 cm⁻¹). Each of these parameters was normalized to the amount of organic matrix (mainly collagen) assessed at a wavelength of 2940 cm⁻¹ (Unal and Akkus, 2015).

2.3.2 Spectra processing

Raman spectra were processed using two different softwares: the Windows®-based Raman Environment (WiRE 3.5, Renishaw) and MATLAB PLS toolbox (PLS_Toolbox 9.0 (23592) Eigenvector Research, Inc., Manson, WA United States 98831). WiRE software was used for peak deconvolution, applying a baseline correction. The latter was used to deconvolve the Amide I subpeaks (1605, 1635, (1657), 1666, 1686 cm⁻¹) and those present in the 2700–3800 cm⁻¹ spectra (2943, 3220, 3333, 3457 cm⁻¹). Models were built by measuring the average centers (intensity wavelength), width and Gaussian percentage of each subpeak derived from the deconvolution. The PLS toolbox automates the measurement of height and area of each peak of interest with a local baseline correction, after the manual creation of a model, according to the mean values of the peaks of interest from all spectra.

2.4 Protein extraction and enzyme-linked immunosorbent assay

Decorin (DCN) content was evaluated using the enzyme-linked immunosorbent assay (ELISA). For this purpose, the protein extraction was performed as previously described (Manon et al., 2022). Samples of each HFL and HP were weighted, cut and lysed with the buffer solution containing the radioimmunoprecipitation assay (RIPA, ThermoFisher Scientific, 89900), the protease inhibitor cocktail (Sigma-Aldrich, P8340-1 ML), and phosphatase inhibitors (PhosSTOP™, Roche, Basel, Switzerland). The samples were subjected to four cycles in the Precellys homogenizer (Bertin Technologies SAS, France) at 7,200 rpm. The protein concentration in the supernatant was determined with the Pierce™ BCA Protein assay kit (ThermoFisher Scientific, 23227) and normalized to the sample weight.

Subsequently, DCN content was evaluated using the RayBio^® Human Decorin ELISA kit (RayBiotech, ELH-Decorin-1,13485) following an adaptation of the manufacturer’s instructions. Briefly, since the protein concentration was low, larger volumes were used to achieve the required concentration per well. Therefore, 145 µL of tissue homogenates and standards were incubated in the pre-coated 96 well microplate at room temperature for 2.5 h. Following the washing steps, biotinylated antibody was added for 1 h. The subsequent washing step removed unbound antibodies before the addition of streptavidin conjugated horseradish peroxidase (HRP) solution. 5′-Tetramethylbenzidine (TMB) substrate was added after four washes and incubated protected from light. After 30 min, the reaction was stopped (Stop solution of 0.2 M sulfuric acid) and the plate was red at 450 nm. Standard curve was obtained using Sigmaplot 13 software (version 20, SPSS, Inc., Chicago, IL, United States) and results were expressed as pg of DCN per mg of total protein (pg/mg).

2.5 GAG quantification

The content of sulfated glycosaminoglycans (GAGs) was quantified in native and decellularized HP and HFL samples (n = 3) using the Blyscan Sulfated Glycosaminoglycan Assay (Biocolor LTD., Carrickfergus, Northern Ireland) and following the manufacturer’s protocol. Sulfated GAGs include chondroitin sulfate, dermatan sulfate, heparan sulfate and keratan sulfate. Three samples of 25 mg of each tissue were harvested. Samples were then first freeze-dried and dry-weighted. Before the quantitative measurement, each sample was placed in a Papain reagent (50 mL Na₂HPO₄/NaH₂PO₄ 0.2 M pH 6.4, 400 mg sodium acetate, 200 mg EDTA and 40 mg cysteine HCL, 250 µL papain) overnight at 65°C for GAG extraction. Following centrifugation (10 min, 10,000 rpm), 100 µL of each sample were taken and 1 mL of Dye Reagent was added for 30 min. The Eppendorfs were centrifuged again at 12,000 rpm for 10 min and the supernatant discarded. Dissociation reagent (500 µL) was added for 10 min and, after a final centrifugation (5 min, 12,000 rpm), 200 µL of each sample were transferred to a 96-well plate. Standards were also prepared according to the manufacturer’s protocol. The plate was read at 630 nm three times and the average of the three readings was used for final quantification. The concentration of GAGs was expressed in µg/mg dry weight.

2.6 Statistics

Statistical analyses were carried out using GraphPad Prism (Version 8.4.0, Inc., San Diego, CA, United States) and SPSS software (Version 27, IBM SPSS, Inc., Chicago, IL, United States). Data were expressed as the mean ± standard error of mean (SEM) and compared with the non-parametric Mann-Whitney or Kruskal-Wallis test because of the small dataset and the subsequent non-respect of normality (tested by Shapiro-Wilk test). Pairwise comparisons were performed using the Dunn’s test. All tests were two-tailed. The significance level was considered at a p-value < 0.05.

3 Results

All the decellularization protocols successfully achieved the Crapo’s criteria (Crapo et al., 2011). These standards include an absence of nuclei visible by histological analysis (DAPI, HE staining), a DNA quantity of less than 50 ng/mg dry weight and a length of less than 200 bp for persistent DNA fragments after treatment. In all the samples, DNA concentration felt below the critical threshold of 50 ng/mg of dry weight and no cellular traces were visible on microscopic analyses.

3.1 Weak alteration of collagen structure

HFL samples appeared to be more disorganized after protocol D5 treatment than the other subgroups, as highlighted by a decrease, although not significant, of thick (1.6–2.4 µm) red hue fibers (ImageJ: −29%; QuPath: −32.6%) and an increase in thin (≤0.8 µm) green ones (ImageJ: +87.9%; QuPath: +53.7%) (Figure 1). HP samples seemed to be more affected by protocol D1. Although not significant, the D1-treated tissues presented the lowest proportion of red fibers (ImageJ: −33.8%; QuPath: −68.3%) and the highest of green ones (ImageJ: +300%, QuPath: +189%) compared to natives. No significant differences in the amount of intermediate fibers (yellow/orange) were recorded.

Figure 1. Qualitative (A) and quantitative (B) evaluations of the collagen organization of sections through native (n = 12) and decellularized HFL (n = 4) and HP (n = 3). (A) Sirius Red staining under bright field (SR, left) and polarized light (PL, right) obtained with ZEISS Axio Scan.Z1 Slide. Same analyses were carried out with polarized light microscope but are not shown. Variation in the polarization of collagen fibers under polarized light reflects differences in fiber thickness, density and packing: from green (thin, ≤0.8 µm, and loosely packed) to red (thick, 1.6–2.4 µm, and tightly packed) (Dayan et al., 1989). (B) Quantitative analyses were performed using ImageJ for the microscopic images and Qupath for scanner images. The quantitative analysis confirmed the qualitative assumptions: D5 and D1 protocols affected the most the HFL and HP samples respectively, with a higher average of thin fibers and lower of thick ones assessed by both ImageJ (microscope) and QuPath (scanner), despite remaining not significant.

3.2 Weak difference in the matrix molecular composition in native tissues

Raman spectroscopy revealed no statistical differences between the molecular composition of native HFL and HP. Both tissues presented similar amounts of hydroxyproline and collagen-bound water, although HFL samples showed a slightly higher proportion of enzymatic crosslinks (+13.6%, p > 0.05) and a slightly lower amount of non-enzymatic crosslinks (−55.8%, p > 0.05) than native HP.

3.3 Weak alteration of the matrix molecular composition in decellularized tissues

3.3.1 Human fascia lata

Increase in tissue hydration was the most marked modification after HFL decellularization (Table 2), at the limit of significance for all three biomarkers of collagen-related water for all D-protocols: the hydroxyproline OH groups beside water bound to collagen (IR 3457/2943, p = 0.054), the NH groups beside organic matrix related water (IR 3333/2943, p = 0.052) and the organic matrix related water parameter (IR 3220/2943, p = 0.061). These increases were then confirmed, by a pairwise Dunn’s comparison, for the three biomarkers in protocols D3 (IR 3220/2943, p = 0.030; IR 3333/2943 p = 0.047 and IR 3457/2943, p = 0.037) and D5 (p = 0.030, p = 0.023 and p = 0.037 respectively) when compared to native tissues. An average loss of 25.2% (p = 0.195) in the amount of enzymatic crosslinks was also observed after the D5 treatment, which was the largest, though not significant, loss among all protocols. Hydroxyproline and pentosidine content were not significantly affected by the chemical treatments as well.

Table 2. Raman spectroscopy data of Human Fascia Lata matrices after decellularization with 5 different protocols (D1-D5) (n = 3/protocol and native). AGEs: Advanced Glycation End Products, NH: amine group, OH: hydroxide. Results are expressed as Mean ± SEM (Kruskal-Wallis).

3.3.2 Human periosteum

HP samples hydration also increased with all decellularization protocols (Table 3) but in a less important way than HFL. D1-HP presented the lowest average of enzymatic crosslinks (−37.4%, p = 0.137) and pentosidine (−46.7%, p = 0.332) among the different decellularization protocols. Once again, the content of hydroxyproline did not significantly vary according to the treatment.

Table 3. Raman spectroscopy data of Human Periosteum matrices after decellularization with 5 different protocols (D1-D5) (n = 3/protocol and native). AGEs: Advanced Glycation End Products, NH: amine group, OH: hydroxide. Results are expressed as Mean ± SEM (Kruskal-Wallis).

3.3.3 Irradiated human fascia lata

By multiple comparison analyses, we highlighted a further loss in the proportion of enzymatic crosslinks in D5i sterilized samples, leading to a significant difference with natives (−29.4%, p = 0.037) (Figure 2). However, hydroxyproline and pentosidine levels were similar to that of native tissues (Table 4). Moreover, while hydration increased after decellularization (D1 to D5) of both tissues (HFL and HP), this proportion decreased after sterilization processes for HFL (D1i to D5i), without returning to basal levels of natives; the data obtained with protocol D5i remained the most different.

Figure 2. Measurement of enzymatic crosslinks by Raman spectroscopy in samples of native HFL and HFL decellularized according to the fifth protocol, followed (D5i) or not (D5) by gamma irradiation. The decellularization chemical treatment impaired the amount of enzymatic crosslinks but without significance (−25.2%, p = 0.195). However, after sterilization, the amount of enzymatic crosslink was further reduced and this reduction became significant (−29.4%, p = 0.037) compared to native.

Table 4. Raman spectroscopy data of Human Fascia Lata matrix after decellularization with 5 different protocols followed by irradiation (D1i-D5i) (n = 3/protocol and native). AGEs: Advanced Glycation End Products, NH: amine group, OH: hydroxide. Results are expressed as Mean ± SEM (Kruskal-Wallis) and significant values are in bold-italic.

3.4 Interfibrillar bridges of decorin

ELISA analyses of protein extracts from HP and HFL samples provide a deeper insight into the effect of chemical treatments on interfibrillar bridges (Figure 3). The DCN levels in HFL were approximately 6.40 ng/mL (±0.23) for all subgroups except D5-tissues, which averaged 2.3 ng/mL (±0.33) (p = 0.036). HP samples of all D-protocols presented a homogenous alteration of DCN levels compared with the native subgroup, with a mean of 2.2 ng/mL (±0.15) versus 6.2 ng/mL (±0.04) for natives. DCN levels in protocol D2 were the lowest, exhibiting a significant decrease compared to natives (p < 0.001), followed by protocol D4 (p = 0.010).

Figure 3. Decorin (DCN) levels (ng/mL) assessed by ELISA in HFL (A) and HP (B) protein extracts according to sample treatments (D1-D5). DCN in HFLs was significantly disrupted by tissue engineering processes in D5-treated tissue (p = 0.036). Concerning HP, DCN level was reduced with all decellularization protocols, but most importantly in D2 (p < 0.001) and D4 (p = 0.010)-treated tissues. Error bars represent the SEM. D: Decellularization protocol.

Since interfibrillar bridges are formed by an association between DCN and sulfated GAG, we wanted to identify a relationship in their persistence within the tissue-engineered samples. In both tissue types, total sulfated GAG levels were similarly and significantly decreased after all D-protocols (Figure 4A). HFL decellularized samples contained average levels of 0.24 µg of GAG per mg of dry weight (±0.20), while native sample levels averaged 3.32 μg/mg (±0.38). For HP samples, a similar ratio was found, with an average of 0.66 μg/mg (±0.19) for decellularized tissues and 3.06 μg/mg (±0.79) for native tissues. No correlation was observed between GAG and DCN levels for both HFL (R² = 0.022) and HP (R² = 0.784) samples (Figure 4B). Changes in DCN and GAG level, resulting from chemical treatments were not related.

Figure 4. Impact of decellularization approaches on GAGs and DCN. (A) Quantitative analyses of GAG content in HFL and HP samples. The levels significantly decreased after all chemical treatments for both HFL (D1 (p = 0.0004), D2 (p = 0.024), D3 (p < 0.0001), D4 (p = 0.012), D5 (p < 0.0001)) and HP (D1 (p = 0.0002), D2 (p = 0.007), D3 (p < 0.0001), D4 (p < 0.0001), D5 (p = 0.007)) samples. (B) Variation of DCN and GAGs resulting from the chemical treatments are not related for both HFL (R² = 0.022) and HP (R² = 0.784). However, all decellurazed samples present similar profiles clearly distinguishable from native tissues. Error bars represent the standard error of the mean (SEM). GAGs: Glycosaminoglycans, D: Decellularization protocol, DCN: decorin.

4 Discussion

None of the five decellularization protocols (D1 to D5, Table 1) significantly altered collagen arrangement. Detergent-free treatment (D5) disorganized HFL samples to a greater extent, significantly reducing DCN and GAGs, and non-significantly reducing collagen cross-links. HP samples, treated with detergent-based treatment D1, showed collagen disorganization, and a reduction in enzymatic cross-links. Overall, decellularization had minimal impact on the matrix, while the addition of gamma irradiation impact leads to a significant negative effect on the integrity of collagen cross-links.

Structural modifications were therefore most pronounced in HFL when treated with the D5 protocol and in HP with the D1 protocol. This was first observed with the assessment of collagen bundle birefringence using PLM to understand the impact of the engineered treatments on the architectural arrangement (Dayan et al., 1989; Ribeiro et al., 2013; Rittié, 2017). The robustness of the results is confirmed by the convergence of birefringence measurements between various analysis tools. Interestingly, each of these subgroups, HFL-D5 and HP-D1, also showed the lowest levels of enzymatic crosslinks as measured by Raman spectroscopy. In the literature, optical density measures under polarized light were already used as indicator of crosslinked collagen organization (Kirby, Heuerman, and Yellon, 2018; Lien, Chen, and Phan, 2022).

Previous studies confirmed collagen persistence in decellularized scaffolds (Duisit et al., 2017; Manon et al., 2023). Being the major ECM component, collagen biological and architectural functions have already been well characterized (Snedeker and Gautieri, 2014; Dong and Lv, 2016). Ensuring the persistence of enzymatic crosslinks in decellularized collagen scaffolds is therefore crucial for maintaining the functional properties of native matrices, contributing to fibril stability and improving tissue biomechanical competence. Therefore, physical and chemical methods are nowadays used to achieve intermolecular collagen crosslinking in engineered scaffolds (Dong and Lv, 2016).

In addition to collagen enzymatic crosslinks, its overall molecular composition has an impact on its structure and biomechanical properties. Indeed, AGEs are known to alter the mechanical strength of tissues (Saito and Marumo, 2010). In addition, as Unal and Akkus (2015) have shown, the hydration level of bones also plays a role in their mechanical properties. Moreover, it is also known that physical training has an impact on collagen crimps, which can be analyzed by measuring birefringence (Mazon et al., 2018). Therefore, by taking all collagen parameters into account, we could in future investigate a relationship between the structural data of our scaffolds and their mechanical properties.

It is important to note, however, that the samples were taken from old subjects, which may also have an impact on the results. Since the rate of crosslinks is age-related, the slight decrease in enzymatic crosslinks after the engineering protocols could be explained by an age-related decrease, itself associated with a decrease in the mechanical capacity of the tissue. These decreases could therefore possibly be significant in tissues from younger subjects, which should be verified in the future. Moreover, it should be noted that these analyses were carried out on a small number of samples randomly taken from the entire decellularized scaffold, allowing potential tissue heterogeneity to be taken into account.

Raman spectra of the HFL samples differed slightly from those of the periosteal ones, particularly in the amide I region (1670/1690 cm⁻¹). We suspected interference from tissues other than periosteum, adding some background in the spectra. The samples were harvested from old subjects and, despite efforts to remove muscle and fat, some residues remained due to the thinness of the HP. Indeed, anomalies in the spectra were observed for all protocols, including native tissue, with the exception of protocol D5, the only one to include a degreasing step. Therefore, in line with the literature, an additional peak (1657 cm⁻¹) was added for deconvolution of this amide region, known to correct for potential interference from lipid vibrations (Mandair et al., 2020). In addition, we observed increased hydration in samples after decellularization, probably because decellularization disrupts the natural barriers regulating hydration. This leads to increased water uptake during abundant rinsing steps. Irradiation, however, leads to a partial loss of this increased hydration, without returning to baseline levels.

Since these decellularized tissues are intended to regenerate in vivo damaged structures, sterilization processes are required to prevent pathogen transmission (Dziedzic-Goclawska et al., 2005; Seto, Gatt, and Dunn, 2008; Harrell et al., 2018). Althoug gamma irradiation is commonly used in bone banks for sterilization, it causes radiolysis of water molecules, releasing free radicals (Nguyen et al., 2007) that damage collagen (Seto, Gatt, and Dunn, 2008; Harrell et al., 2018). Previous studies on bones showed a negative synergic effect on mechanical resistance when irradiation was added to freeze-dried samples (Cornu et al., 2000; Cornu et al., 2011), attributed to a significant reduction in enzymatic crosslinks density (Cornu, 2012). This is consistent with our results revealing a more pronounced alteration of crosslinks post-irradiation. Future studies should explore decellularization procedures under totally sterile conditions to avoid the detrimental effects of subsequent irradiation on tissue properties.

We did not investigate irradiated periosteum since the priority, for the creation of the new regenerative membrane, was the use of irradiated decellularized HFL (Manon et al., 2023).

DCN-GAG complexes are essential for the ECM organization, especially in collagen fibrillogenesis, as evidenced by the alteration of the latter in DCN-deficient mice (Gubbiotti et al., 2016; Zhang et al., 2018). This PG is also involved in stress transfer between collagen fibrils (Thorpe, 2010) and in the bioavailability of growth factors such as TGF-β (Gubbiotti et al., 2016; Zhang et al., 2018). Alteration in decorin content vary according to the tissue type, with HFL samples globally maintaining their DCN concentration after decellularization, except for D5, while HP samples show a decrease in all D-protocols.

Furthermore, it has already been demonstrated that detergent-based decellularization affects the GAG composition of ECM (Duisit et al., 2017; Liguori et al., 2020; Uhl et al., 2020), as confirmed by our quantitative measurements, with potential dysfunction of persistent GAGs (Uhl et al., 2020). However, the maintenance of functional GAGs is necessary. Abundantly distributed on cell surface and in the ECM, they modulate various physiological processes. They provide biochemical and mechanical support in the ECM and influence cell behavior (adhesion, proliferation, differentiation) (Salbach et al., 2012; Menezes et al., 2022).

Despite the lower proportions of DCN and GAG in decellularized samples compared to native ones, no correlation could be demonstrated in their alteration after chemical treatments. Since GAGs are ubiquitous on the cell surface, their removal by the decellularization protocol is expected. Therefore, proportionally, the loss of GAGs could be greater than that of DCN.

The present study confronts different decellularization approaches. Protocols D1 to D4 are based on the use of Triton, SDS and DNase and ensure cell membrane lysis, protein solubilization and residual DNA elimination. Triton and SDS are known for their efficiency in decellularization, but also for altering the structure of the ECM (Crapo et al., 2011). Protocol D5 consists of a succession of solvent- and detergent-based baths, which has also already demonstrated its decellularization potential (Manon et al., 2023), including the monitoring of Crapo’s criteria (Crapo et al., 2011). In this protocol, solvents such as ethanol and acetone are used for cell lysis and lipid removal but are both also known as altering the ECM ultrastructure. The additional advantage of this protocol is that it also neutralizes prions, which the other protocols do not allow.

As qualitatively shown in the Figure 1 on the SR stained sections, more undulating collagen structure are observed in decellularized tissues than in native tissues. The crosslink alteration could modifiy the fibers alignment. In the future, analysis of the collagen distribution could support the observations made in this study. Approaches such as transmission electron microscopy (Adeoye et al., 2022), histology and contrast-enhanced microCT (Maes et al., 2022) could be used and be gathered in order to precisely cartograph the entire HFL structure.

In conclusion, a relationship was observed between the ECM structure and its molecular composition. The decellularization approaches demonstrated no significant alteration of the matrices, which justifies, along with the analysis of other parameters, the future clinical use of these processed matrices for bone regeneration. However, the combination of D5 treatment with subsequent irradiation appears the most deleterious for tissue properties. In the future, an analysis of the entire sample surface would provide more accurate information about potential sample heterogeneity, and precise quantitative analyses of the crosslink content could further validate our results. It would also be interesting to correlate all the observations of this study with the mechanical properties of the tissues.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by the Human Anatomy Department of UCLouvain (Body donation service, IRB00008535, Brussels, Belgium. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

JV: Writing–original draft, Software, Methodology, Investigation, Data curation. JM: Writing–review and editing, Writing–original draft, Supervision, Resources, Methodology, Investigation, Funding acquisition, Data curation, Conceptualization. AC: Writing–review and editing, Software, Data curation. RE: Writing–review and editing. LF: Writing–review and editing, Methodology, Investigation, Data curation. TS: Writing–review and editing. BL: Conceptualization, Writing–review and editing, Supervision. CB: Conceptualization, Writing–review and editing, Supervision. OC: Supervision, Writing–review and editing, Conceptualization.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the F.S.R. Fund (Fonds Spécial pour la Recherche, Belgium) and F.N.R.S. Aspirant Fund granted to JM (Fonds National de la Recherche Scientifique, Application ID 40004991, Belgium), and FRIA grant awarded to JV (Fonds pour la Formation à la Recherche dans l’Industrie et dans l’Agriculture, Application ID 40021733, Belgium).

Acknowledgments

The authors would like to express their sincere gratitude to all the colleagues for their contribution to this work, including Xavier Merny, Olivier Goffin, Laurence Lekeux, Dimitri Roman, Angélique Arts, Shalaw Fawzi, Andreia Teixeira da Cunha, Dimitri Pinto de Carvalho, Anthony Marais, Paulo Ferreira Gomes and all the members of the institutional tissue bank, Guillaume Falgayrac, Caroline Bouzin and the 2IP platform, Christine de Ville de Goyet, Bernard Caelen, Nicolas Charles-Pirlot, Michelle Cougnon, Walter Hudders, Grégoire André. The authors would like to extend special gratitude to the donors used in this study, along with their families.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References