3D Porous Architecture of Stacks of β-TCP Granules Compared with That of Trabecular Bone: A microCT, Vector Analysis, and Compression Study

The 3D arrangement of porous granular biomaterials usable to fill bone defects has received little study. Granular biomaterials occupy 3D space when packed together in a manner that creates a porosity suitable for the invasion of vascular and bone cells. Granules of beta-tricalcium phosphate (β-TCP) were prepared with either 12.5 or 25 g of β-TCP powder in the same volume of slurry. When the granules were placed in a test tube, this produced 3D stacks with a high (HP) or low porosity (LP), respectively. Stacks of granules mimic the filling of a bone defect by a surgeon. The aim of this study was to compare the porosity of stacks of β-TCP granules with that of cores of trabecular bone. Biomechanical compression tests were done on the granules stacks. Bone cylinders were prepared from calf tibia plateau, constituted high-density (HD) blocks. Low-density (LD) blocks were harvested from aged cadaver tibias. Microcomputed tomography was used on the β-TCP granule stacks and the trabecular bone cores to determine porosity and specific surface. A vector-projection algorithm was used to image porosity employing a frontal plane image, which was constructed line by line from all images of a microCT stack. Stacks of HP granules had porosity (75.3 ± 0.4%) and fractal lacunarity (0.043 ± 0.007) intermediate between that of HD (respectively 69.1 ± 6.4%, p < 0.05 and 0.087 ± 0.045, p < 0.05) and LD bones (respectively 88.8 ± 1.57% and 0.037 ± 0.014), but exhibited a higher surface density (5.56 ± 0.11 mm²/mm³ vs. 2.06 ± 0.26 for LD, p < 0.05). LP granular arrangements created large pores coexisting with dense areas of material. Frontal plane analysis evidenced a more regular arrangement of β-TCP granules than bone trabecule. Stacks of HP granules represent a scaffold that resembles trabecular bone in its porous microarchitecture.

Introduction

Bone is a calcified connective tissue, which fulfils a large number of functions in the body (1). Being both rigid and elastic due to the unique composition of its matrix made of collagen type I and hydroxyapatite crystals, it is adapted to gravity and muscle contractions. Bone cells are able to remodel and adapt bone mass and architecture throughout life allowing repair of microdamage (microcracks, trabecular microfractures) or fracture. To reduce the amount of this material within the body, Nature has developed a number of strategies to provide the maximum biomechanical competence with a genetically controlled bone mass. During aging, epigenetic factors, such as gravity and muscle strains, also modify bone mass and microarchitecture (Wolff’s Law). Trabecular microarchitecture is now a well-recognized type of self-organization that fulfils this condition (2, 3).

Bone loss can occur systemically, for example, during aging. However, localized bone loss can also occur, as in some intra-bony defects, such as solitary or aneurysmal cysts, or acquired bone defects as observed after tooth extraction or during hip prosthesis revision. Granular biomaterials are often used to fill these defects, especially in non-bearing bones. Among the different orthophosphates that can be produced in large quantity by industry, beta-tricalcium phosphate (β-TCP) appears the most suitable, particularly in maxillo-facial and dental surgery (4). Granules of 1000–2000 μm diameter represent the most commonly used size for filling alveolar sockets or increasing bone volume by sinus lift (5). However, when the surgeons place granules within a bone defect, it is the voids between the granules that represent the interconnected space available for vascular sprouts and bone cells to invade the grafted area. Calcium phosphate ceramics are known to be too brittle to be placed in load-bearing areas. Ceramic granules allow a perfect filling of the bone defects, but they are limited to grafted areas where no significant strength is needed (alveolar ridge augmentation, sinus lift …) (5). Little is known about the 3D arrangement of a grafted stack of granules, and if the porosity of the stack is similar to that of trabecular bone to allow the accessibility of biological fluids, body fluids, cells, and vascularization within the grafted site (6, 7). In addition, when filling a bone defect, the surgeon has to avoid crushing of the granules to avoid altering the granules’ microarchitecture (8). Hence, the biomechanical resistance of the granules in compression is also an important factor to be taken into account when choosing the grafting material (9).

The aims of the present study were (1) to compare the pore microarchitecture of the two types of β-TCP granules providing stacks of low porosity (LP) and high porosity (HP), respectively, with that of the two types of trabecular bone [high density (HD) and low density (LD)], in order to assess the granular stack that had the pore microarchitecture more similar to trabecular bone; (2) to compare the biomechanical properties in compression of these two types of granules.”

Materials and Methods

Bone Samples

For the bone with HD, five bone cylinders were harvested from the upper tibia epiphysis of young butchered animals (calves, aged around 15 months, 500–600 kg). The calibrated cylinders (10 mm in diameter) were obtained with a diamond-tipped coring tool (Starlite Industries, Rosemont, PA, USA) perpendicular to the tibial plateau. Similarly, five cylinders (10 mm diameter) with low bone density (LD) were prepared from the upper tibial epiphysis from five aged human cadavers from the Anatomical Laboratory (mean age 78.3 ± 5.2 years). For both the HD and LD cylinders, the articular and subchondral bone were removed to provide cylinders 23–25 mm in length containing only the secondary spongiosa. The subjects had given their body to science before death and had completed a form indicating that their corpses were to be used for medical education and research. Subjects were examined upon arrival at the laboratory of anatomy of the Faculty of Medicine. All bone samples were harvested following a strict protocol of the laboratory of anatomy; e.g., it was impossible to know any detail which could lead to the subject’s identification (causes of death, previous medical records).

β-TCP Granules

Beta-tricalcium phosphate granules were obtained from Kasios (L’Union, France). Granules were prepared by polyurethane foam technology and used as received. This technique has been extensively described elsewhere (10, 11). Two types of granules were used and placed in polyethylene test tubes (10 mm in diameter): compact granules were prepared with 25 g of β-TCP grain in a fixed volume of slurry and provided a LP stacks; lighter granules with 12.5 g of β-TCP grain 1 in the same volume of slurry produced the HP stacks (12).

The morphology of the two types of granules was qualitatively evaluated by scanning electron microscopy (SEM, Jeol JSM 6301F) before analysis. The β-TCP granules were glued on an aluminum stub with a double-sided carbon sellotape. They were gold-coated by sputtering with an MED 020 (Bal-Tec, Chatillon sur Cher, France). Images were captured at a 3-kV acceleration voltage in the secondary electron mode.

Biomechanical Analysis

The compression characteristics of the two types of β-TCP granules were determined on an Instron 5942 machine with Bluehill software (Instron, Elancourt, France) (Figure 1A). As the test is destructive, new intact granules were placed into the cylindrical test cell of the machine and compacted with an adapted punch 10 mm in diameter. The maximum applied load was fixed at 400 N with a 2-mm min⁻¹ displacement. For each type of stacks of granule (HP or LP), the load–displacement curve was recorded (N/mm) automatically by the machine. Stiffness was automatically provided by measuring the slope of the tangent to the curve at the beginning of the deformation (N/mm). The maximum displacement of the punch (in mm) was obtained when the load reached 400 N. The work to failure (N⋅mm) was determined as the area under the load–displacement curve from yield until failure and represents the energy absorbed by bone during plastic deformation (Figure 1B) (13). These parameters were calculated according to previously published equations (14, 15). Analysis was repeated five times for each type of granules and the averaged values were considered.

FIGURE 1

Figure 1. Biomechanical analysis of the two types of stacks of granule (HP and LP). (A) Image of the compression machine; granules are placed in the cylindrical test cell (1) and compressed with the punch (2). (B) Theoretical load–displacement curve between displacement and force obtained by compression; the slope of the tangent to the curve is automatically provided (blue line), thus providing stiffness; the maximum displacement of the punch at 400 N is obtained by the red line; the area under the curve (in green) represents the work to failure. (C) The curves obtained for the HP and LP stacks of granules.

MicroCT Analysis

MicroCT analysis was performed on a Skyscan 1172 microcomputed tomograph (Bruker microCT, Kontich, Belgium). Bone samples were placed in polystyrene test tubes (10 mm in diameter) and scanned while in fixative (10% formalin). Granules of β-TCP were placed in similar tubes and gently agitated to allow the granules to settle and optimize their 3D spatial arrangement in the dry state. Tubes were filled with β-TCP granules up to 30 mm in height. For each series of granules, analysis was done in triplicate. The microCT was operated at 80 kV, 100 μA with a 0.5-mm aluminum filter. The pixel size was fixed at 4.95 μm and the angular rotation step at 0.25°. 3D models of bones and granules were obtained with VG Studiomax (Volume Graphics GmbH, Heidelberg, Germany) operating in the volume rendering mode. Microarchitecture parameters were obtained with the CTAn software (Bruker) after global thresholding of the pores (16). They included measures of porosity (Po, in %), mean pore diameter (Po⋅Diam, in μm), and specific surface [representing the interface between bone or biomaterial and the porosity, BS/TV, in mm²/mm³, where BS is the whole pore surface and TV is the volume of interest according to the ASBMR nomenclature (17)]. For the bone cylinders and the granule stacks, the region of interest (ROI) over-imposed on the 2D binarized sections was a 5-mm wide square; the height of analysis of the samples has concerned 15 mm. Pore diameter was determined by the sphere-fitting method in the CTAn software (18). In order to further investigate the porosity, the binarized stacks of 2D microCT sections of each sample were analyzed by software written in Matlab (MathWorks, Natick, MA, USA). The algorithm for evaluation of porosity used a vector projection on a frontal plane (Vectopor) and has been extensively described elsewhere (19). In the present study, the pixels of the ROI which belonged to the same column received the same pseudo-color according to the number of pixels superimposed on pores. A look-up table (LUT) was designed ranging from deep blue (no porosity) to red (high porosity). Intermediate values depend on the value of porosity along the vector according to the LUT. The frontal plane vector projection image was assigned intensity values for each pixel based on the number of pixels, which contained pores in the related column of the microCT binarized image. The frontal plane image was saved with the colorized LUT and analyzed by the FracLac plugin developed for ImageJ (20, 21). The box-plot fractal dimension (D_f) and the lacunarity (λ) of the frontal plane images were determined as previously described (19). The box-plot fractal dimension evaluates the complexity of the material to fill the reference space and is helpful in the study of porous objects, such as bone (22, 23) or materials (24, 25). Lacunarity is another recently described fractal parameter reported to improve the description of a fractal porous object (26, 27). Lacunarity is influenced by the variation of the pores in a given structure; a low lacunarity reflects homogeneity, while a high lacunarity is an indicator of heterogeneity (28, 29).

Statistical Analysis

Statistical analysis was done with Systat 13 (Systat Inc., San José, CA, USA). For each parameter, results are expressed as mean ± SD. Differences for porosity between the four groups (HD, LD, HP, and LP) were sought with the Kruskal–Wallis non-parametric analysis of variance followed by the Conover–Inman test for all pairwise comparisons. For biomechanical parameters, differences between the HP vs. LP groups were analyzed by the non-parametric Mann–Whitney U test. A difference was considered as significant when p < 0.05.

Results

Qualitative Examination of the β-TCP Granules

The difference in morphology of the granules was evident to the naked eye and confirmed by SEM. The 12.5-g preparation of granules was the most porous with thin walls of sintered β-TCP, while the 25-g granule preparation was denser, with numerous concave surfaces and occasional macropores (Figure 2). In all cases, the inner porosity due to sublimation of the polyurethane foam was clearly identified as triangular voids connected by thin channels.

FIGURE 2

Figure 2. SEM analysis of the two types of β-TCP granules. (A) LP granules prepared with 25 g in the slurry, (B) HP granules prepared with 12.5 g. Note the inner porosity due to the sublimation of the polyurethane foam as triangular voids connected by thin channels (arrows).

Biomechanical Analysis of the β-TCP Granules

The load vs. displacement curves of both HP and LP types of granules are shown in Figure 1C. The first domain of the curve corresponds to the collapse of the granules, after which the deformation is due to their irreversible compression. The HP stack had a statistically significant higher maximum displacement compared to the LP stack, as expected (p < 0.05, Table 1). On the contrary, the stiffness and the work to failure were significantly higher in the LP stacks compared to the HP stack (p < 0.05 in both cases).

TABLE 1

Table 1. Morphometric parameters of bone cylinders with a high-density (HD), low-density (LD), and β-TCP granules (compact low porosity granules: LP; high porosity granules: HP).

Morphometric Analysis of Bone and Granules

A 3D-rendering of the two types of bone cores and granule stacks is shown in Figure 3. Morphometric parameters are listed in Table 1 and illustrated in Figure 4. The HD bone cores had a significantly smaller porosity than LD bones. The LP stacks had the lowest porosity among all groups. Conversely, HP stacks presented a porosity that is intermediate between that of a HD bone and that of a low-density bone. The HP formulation produced a ~28% increase of porosity associated with a reduction of the mass of β-TCP. This was associated with a significant increase in BS/TV in the HP group compared to the LP group (+11.6%). The vector analysis of porosity clearly evidenced differences between the two pairs of materials: HD bones appeared inhomogeneous with areas containing large amounts of bony material (blue areas according to the LUT), while LD had a large number of red spots indicating HP (Figures 5A,B). Very similar findings were observed for the LP and HP stacks of granules: LP stacks had large areas of material but large holes were also present. On the contrary, the morphology of the stacks of HP granules had a more regular distribution of the pores and material (Figures 5C,D). Interestingly, the fractal dimension D_f did not differ among the four groups, but lacunarity λ was high in HD bones and even higher in the LP granules (Table 1).

FIGURE 3

Figure 3. Microtomographic analysis of the different types of materials (appearing in white). (A) High-density bone (HD) with numerous trabeculae and low porosity; (B) low-density bone (LD) from aged subjects with reduced trabeculae and increased porosity; (C) of compact granules (LP) of β-TCP; (D) scaffold of HP granules. All images have been prepared with a square ROI of 5 mm side; so the volume of interest is a parallelepiped (height: 15 mm; square section: 5 mm).

FIGURE 4

Figure 4. Porosity measured by microCT in the bone cylinders and stacks of granules. Statistical analysis was done with Kruskal–Wallis analysis of variance followed by the Conover–Inman post hoc test. Error bars indicate SD. *Significantly different with p < 0.02; **with p < 0.003; ***with p < 0.0001.

FIGURE 5

Figure 5. Images of the frontal planes projected from the different types of porous materials: (A) HD bone, (B) LD bone, (C) LP, and (D) HP granules of β-TCP. The regions with maximal porosity are red, minimal porosity in blue according to the LUT (at the bottom of the image). All images have been prepared with a square ROI of 5 mm side length and the scale bar represents 1 mm.

Discussion

In this series of materials comprising bone cores and granular stacks made of β-TCP granules, a quantitative and morphological analysis of porosity was done. Porosity of the LP stacks was noticeably lower than a HD trabecular bone. However, the porosity of the HP stacks (and conversely the amount of material) was in the range of normal bone volume as determined by classical histomorphometry methods (30). This reflects the microarchitecture of the stacks of HP granules which contain more macropores and thin branches than the LP granules. However, as expected, this effect is detrimental to the biomechanical properties of the HP granules whose stacks have lower compressive mechanical properties, as evidenced by lower stiffness values, lower work to failure, and maximum displacement. The surgeon should be aware not to crush the granules in the grafted areas and the use of a bone syringe with a large diameter may help.

It should be noted that the BS/TV of HP granules is significantly higher than that of LP granules or that of HD and LD bones. This is especially important since BS/BV represents the bioactive surfaces of β-TCP that will be in direct contact with the biological fluids and the cells (31). The BS/TV values obtained in this study with β-TCP granules are similar to those reported by other groups with BCP (a mixture of β-TCP and hydroxyapatite) (32). Recent observations have stressed the interest of increasing the specific surface (BS/TV) of orthophosphate materials to promote bone formation (33). So, the HP formulation represents the highest surface available for osteoblast attachment and proliferation (osteoblasts are bone-forming cells). Because β-TCP is a biomaterial colonized by osteoconduction from nearby trabeculae, this parameter is of great importance (34).

MicroCT analysis is a relatively recently introduced tool in the analysis of bone and materials. Values for parameters, such as BV/TV and Po, are highly correlated with histomorphometry on 2D sections (35). In the present study, the difference in porosity between the stacks of granules can be clearly identified as previously reported (12). Different types of parameters have been proposed to evaluate the 3D arrangement of porous materials (e.g., the box-plot fractal dimension) (22, 23) and the interconnectivity of the pores created between the granules (29). The method was successfully applied to evaluate the complexity of different types of granule aggregates: sandstone (36), coal (37), metal grains (38), and soils (39). The complexity of the granule repartition, evaluated by D_f, did not differ between the four groups. However, fractal lacunarity appears as an interesting new tool capable to differentiate the groups but it does not bring morphological information (12, 28, 29, 40). We have recently developed vector analysis for visualization and quantification of the porosity in objects presenting a complex form (19). The method was found useful to analyze pores in materials prepared with various types of porogens and in bones having a particularly complex shape. This is the case with alveolar bone at the mandible, which is very difficult to analyze because of the presence of the tooth roots. Recently, the bone loss induced by paralysis of Mus masseter and Mus temporalis was evidenced by this technique (41). The frontal plane is produced line by line during the vector analysis of each microCT image and represents the projection of porosity. The complexity of the image can be measured by fractal geometry on this projected 2D image. In the present series, there was no difference for D_f among the four groups under study. However, lacunarity is known to be influenced by the uniformity of the spatial distribution; it is a useful parameter when objects do not differ by their fractal dimension (42, 43). The maximal λ value was reached in the LP (25 g) group meaning that very large pores are created in the 3D arrangement of the granules. On the other hand, the lowest λ values were obtained for the HP (12.5 g group) granules and the LD bone, which presented a regular distribution of pore dimensions. Similar relationships have been reported in studies concerned with the porosity of soils (44, 45).

The presence of pores larger than 100 μm has been recognized as a key factor in the vascular invasion of grafted materials into bone (46, 47). In this study, Po⋅Diam values for bone were within the values typically reported for bone (48); Po⋅Diam of both HP and LP are larger than 300 μm and well suited for body fluid and cell invasion (49).

There are some limitations to this study: (1) no compressive analysis was done on the bone cylinders, but synthetic ceramics are well known to be more brittle than bone; (2) bovine bone was used to obtain the HD cylinders, but the biomechanical parameters for trabecular bone are very close (stiffness in compression for bovine bone: 173 MPa vs. human: 139–472 MPa) (50); and (3) the other microarchitectural characteristics of the β-TCP materials having been reported elsewhere are not duplicated in the present study (12).

Conclusion

The present study shows that the two types of β-TCP granules prepared differently did not produce the same interconnected porosity: LP granules provided very large but less numerous pores and HP granules provided a HP with a large interface. The HP stacks of granules presented a porosity similar to trabecular bone, although the granules were physically independent. This study confirms that vector analysis is a suitable method for the evaluation of porosity of complex objects.

Conflict of Interest Statement

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

Acknowledgments

The authors thank Mrs. Laurence Lechat for secretarial assistance. This work was made possible by grants from ANR, program LabCom “NextBone.” The authors thank Kasios, 18 chemin de la Violette, 31240 L’Union – France for providing the different formulations of granules. The analysis software (VECTOPOR) is now licensed by the Agence de Protection des Programmes (APP). Authors also thank Mrs. Nadine Gaborit for her skillful technical assistance with microCT and Phil Salmon for having reviewed the manuscript.

References

1. Baslé MF. Le tissu osseux: ses niveaux d’organisation (the organisation levels of bone). Tissu osseux et biomatériaux en chirurgie dentaire. Paris: Quintessence International (2014). p. 67–80.

Google Scholar

2. Chappard D, Baslé MF, Legrand E, Audran M. New laboratory tools in the assessment of bone quality. Osteoporos Int (2011) 22:2225–40. doi: 10.1007/s00198-011-1573-6

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Parfitt AM, Matthews CHE, Villanueva AR, Kleerekoper M, Frame B, Rao DS. Relationships between surface, volume and thickness of iliac trabecular bone in aging and in osteoporosis. Implications for the microanatomic and cellular mechanisms of bone loss. J Clin Invest (1983) 72:1396–409. doi:10.1172/JCI111096

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Dorozhkin SV. Biphasic, triphasic and multiphasic calcium orthophosphates. Acta Biomater (2012) 8:963–77. doi:10.1016/j.actbio.2011.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Chappard D, Guillaume B, Mallet R, Pascaretti-Grizon F, Baslé MF, Libouban H. Sinus lift augmentation and beta-TCP: a microCT and histologic analysis on human bone biopsies. Micron (2010) 41:321–6. doi:10.1016/j.micron.2009.12.005

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Gauthier O, Bouler J-M, Aguado E, Pilet P, Daculsi G. Macroporous biphasic calcium phosphate ceramics: influence of macropore diameter and, macroporosity, percentage, on, bone, ingrowth. Biomaterials (1998) 19:133–9. doi:10.1016/S0142-9612(97)00180-4

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Holmes RE, Bucholz RW, Mooney V. Porous hydroxyapatite as a bone graft substitute in diaphyseal defects: a histometric study. J Orthop Res (1987) 5:114–21. doi:10.1002/jor.1100050114

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Long M, Martin C, Cooper M. Confined compression of bone substitute granules. Key Eng Mater (2005) 284:361–4. doi:10.4028/www.scientific.net/KEM.284-286.361

CrossRef Full Text | Google Scholar

9. Rezwan K, Chen Q, Blaker J, Boccaccini AR. Biodegradable and bioactive porous polymer/inorganic composite scaffolds for bone tissue engineering. Biomaterials (2006) 27:3413–31. doi:10.1016/j.biomaterials.2006.01.039

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Schwartzwalder K, Somers H, Somers AV, inventors. Method of Making Porous Ceramics. US patent 3 090 094 (1963).

Google Scholar

11. Filmon R, Retailleau-Gaborit N, Brossard G, Grizon-Pascaretti F, Baslé MF, Chappard D. MicroCT and preparation of ß-TCP granular material by the polyurethane foam method. Image Anal Stereol (2009) 28:103–12. doi:10.5566/ias.v28.p103-112

CrossRef Full Text | Google Scholar

12. Ndiaye M, Terranova L, Mallet R, Mabilleau G, Chappard D. Three-dimensional arrangement of beta-tricalcium phosphate granules evaluated by microcomputed tomography and fractal analysis. Acta Biomater (2015) 11:404–11. doi:10.1016/j.actbio.2014.09.015

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Kadiri MS. Compression de poudres pharmaceutiques et interaction avec l’outillage [PhD thesis]. Toulouse: Analyse expérimentale et modélisation numérique Univ. (2004). Available from: http://ethesis.inp-toulouse.fr/archive/00000096/01/kadiri.pdf

Google Scholar

14. Ritchie RO, Koester KJ, Ionova S, Yao W, Lane NE, Ager JW III. Measurement of the toughness of bone: a tutorial with special reference to small animal studies. Bone (2008) 43:798–812. doi:10.1016/j.bone.2008.04.027

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Turner CH, Burr DB. Basic biomechanical measurements of bone: a tutorial. Bone (1993) 14:595–608. doi:10.1016/8756-3282(93)90081-K

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Tassani S, Korfiatis V, Matsopoulos G. Influence of segmentation on micro-CT images of trabecular bone. J Microsc (2014) 256:75–81. doi:10.1111/jmi.12159

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Dempster DW, Compston JE, Drezner MK, Glorieux FH, Kanis JA, Malluche H, et al. Standardized nomenclature, symbols, and units for bone histomorphometry: a 2012 update of the report of the ASBMR Histomorphometry Nomenclature Committee. J Bone Miner Res (2013) 28:2–17. doi:10.1002/jbmr.1805

CrossRef Full Text | Google Scholar

18. Hildebrand T, Ruegsegger P. A new method for the model-independent assessment of thickness in three-dimensional images. J Microsc (1997) 185:67–75. doi:10.1046/j.1365-2818.1997.1340694.x

CrossRef Full Text | Google Scholar

19. Chappard D, Stancu IC. Porosity imaged by a vector projection algorithm correlates with fractal dimension measured on 3D models obtained by microCT. J Microsc (2015) 258:24–30. doi:10.1111/jmi.12212

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Karperien A. FracLac for ImageJ. Charles Sturt University (1999–2013). Available from: http://rsb.info.nih.gov/ij/plugins/fraclac/FLHelp/Introduction.htm

Google Scholar

21. Rasband WS. ImageJ. Bethesda, MD (1997–2014). Available from: http://imagej.nih.gov/ij/

Google Scholar

22. Chappard D, Legrand E, Haettich B, Chales G, Auvinet B, Eschard JP, et al. Fractal dimension of trabecular bone: comparison of three histomorphometric computed techniques for measuring the architectural two-dimensional complexity. J Pathol (2001) 195:515–21. doi:10.1002/path.970

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Fazzalari NL, Parkinson IH. Fractal properties of subchondral cancellous bone in severe osteoarthritis of the hip. J Bone Miner Res (1997) 12:632–40. doi:10.1359/jbmr.1997.12.4.632

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Haire TJ, Hodgskinson R, Ganney PS, Langton CM. A comparison of porosity, fabric and fractal dimension as predictors of the Young’s modulus of equine cancellous bone. Med Eng Phys (1998) 20:588–93. doi:10.1016/S1350-4533(98)00063-0

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Hu BX, Meerschaert MM, Barrash W, Hyndman DW, He C, Li X, et al. Examining the influence of heterogeneous porosity fields on conservative solute transport. J Contam Hydrol (2009) 108:77–88. doi:10.1016/j.jconhyd.2009.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Smith TG Jr, Lange GD, Marks WB. Fractal methods and results in cellular morphology – dimensions, lacunarity and multifractals. J Neurosci Methods (1996) 69:123–36. doi:10.1016/S0165-0270(96)00080-5

PubMed Abstract | CrossRef Full Text | Google Scholar

27. de Melo RHC. Using Fractal Characteristics Such as Fractal Dimension, Lacunarity and Succolarity to Characterize Texture Patterns on Images. Fluminense: Federal University of Fluminense (2007).

Google Scholar

28. de Melo RHC, Vieira EA, Conci A. Comparing two approaches to compute Lacunarity of mammograms. In: 13th Int Conf on Systems, Signal and Image Processing and Semantic Multimodal Analysis of Digital Media (IWSSIP 2006). Budapest: IEEE Signal Processing Society (2006). p. 299–302.

Google Scholar

29. N’Diaye M, Degeratu C, Bouler JM, Chappard D. Biomaterial porosity determined by fractal dimensions, succolarity and lacunarity on microcomputed tomographic images. Mater Sci Eng C Mater Biol Appl (2013) 33:2025–30. doi:10.1016/j.msec.2013.01.020

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Chappard D, Josselin N, Rouge-Maillart C, Legrand E, Baslé MF, Audran M. Bone microarchitecture in males with corticosteroid-induced osteoporosis. Osteoporos Int (2007) 18:487–94. doi:10.1007/s00198-006-0278-8

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Webb K, Hlady V, Tresco PA. Relationships among cell attachment, spreading, cytoskeletal organization, and migration rate for anchorage-dependent cells on model surfaces. J Biomed Mater Res (2000) 49:362–8. doi:10.1002/(SICI)1097-4636(20000305)49:3<362::AID-JBM9>3.0.CO;2-S

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Frenken J, Bouwman W, Bravenboer N, Zijderveld S, Schulten E, Ten Bruggenkate C. The use of Straumann^® Bone Ceramic in a maxillary sinus floor elevation procedure: a clinical, radiological, histological and histomorphometric evaluation with a 6-month healing period. Clin Oral Implants Res (2010) 21:201–8. doi:10.1111/j.1600-0501.2009.01821.x

CrossRef Full Text | Google Scholar

33. Habibovic P, Yuan H, van den Doel M, Sees TM, van Blitterswijk CA, de Groot K. Relevance of osteoinductive biomaterials in critical-sized orthotopic defect. J Orthop Res (2006) 24:867–76. doi:10.1002/jor.20115

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Araujo MG, Liljenberg B, Lindhe J. Beta-tricalcium phosphate in the early phase of socket healing: an experimental study in the dog. Clin Oral Implants Res (2010) 21:445–54. doi:10.1111/j.1600-0501.2009.01876.x

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Chappard D, Retailleau-Gaborit N, Legrand E, Baslé MF, Audran M. Comparison insight bone measurements by histomorphometry and microCT. J Bone Miner Res (2005) 20:1177–84. doi:10.1359/JBMR.050205

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Katz AJ, Thompson A. Fractal sandstone pores: implications for conductivity and pore formation. Phys Rev Lett (1985) 54:1325. doi:10.1103/PhysRevLett.54.1325

CrossRef Full Text | Google Scholar

37. Friesen W, Mikula R. Fractal dimensions of coal particles. J Colloid Interface Sci (1987) 120:263–71. doi:10.1016/0021-9797(87)90348-1

CrossRef Full Text | Google Scholar

38. Khanbareh H, Wu X, van der Zwaag S. Analysis of the fractal dimension of grain boundaries of AA7050 aluminum alloys and its relationship to fracture toughness. J Mater Sci (2012) 47:6246–53. doi:10.1007/s10853-012-6544-6

CrossRef Full Text | Google Scholar

39. Ding W. Fractal dimension features of soil aggregates distribution with the different reclamation years on loess plateau. Sens Lett (2012) 10:555–61. doi:10.1166/sl.2012.1897

CrossRef Full Text | Google Scholar

40. de Melo RHC, Conci A. How succolarity could be used as another fractal measure in image analysis. Telecommun Syst (2013) 52:1643–55. doi:10.1007/s11235-011-9657-3

CrossRef Full Text | Google Scholar

41. Kun-Darbois JD, Libouban H, Chappard D. Botulinum toxin in masticatory muscles of the adult rat induces bone loss at the condyle and alveolar regions of the mandible associated with a bone proliferation at a muscle enthesis. Bone (2015) 77:75–82. doi:10.1016/j.bone.2015.03.023

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Zeng Y, Payton R, Gantzer C, Anderson S. Fractal dimension and lacunarity of bulk density determined with X-ray computed tomography. Soil Sci Soc Am J (1996) 60:1718–24. doi:10.2136/sssaj1996.03615995006000060016x

CrossRef Full Text | Google Scholar

43. Mandelbrot B. The Fractal Geometry of Nature. New York, NY: Freeman, WH (1977).

Google Scholar

44. Chun HC, Giménez D, Yoon SW. Morphology, lacunarity and entropy of intra-aggregate pores: aggregate size and soil management effects. Geoderma (2008) 146:83–93. doi:10.1016/j.geoderma.2008.05.018

CrossRef Full Text | Google Scholar

45. Latief FDE, Fauzi U. Kozeny-Carman and empirical formula for the permeability of computer rock models. Int J Rock Mech Min Sci (2012) 50:117–23. doi:10.1016/j.ijrmms.2011.12.005

CrossRef Full Text | Google Scholar

46. Karageorgiou V, Kaplan D. Porosity of 3D biomaterial scaffolds and osteogenesis. Biomaterials (2005) 26:5474–91. doi:10.1016/j.biomaterials.2005.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Mastrogiacomo M, Scaglione S, Martinetti R, Dolcini L, Beltrame F, Cancedda R, et al. Role of scaffold internal structure on in vivo bone formation in macroporous calcium phosphate bioceramics. Biomaterials (2006) 27:3230–7. doi:10.1016/j.biomaterials.2006.01.031

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Chappard D, Legrand E, Pascaretti C, Baslé MF, Audran M. Comparison of eight histomorphometric methods for measuring trabecular bone architecture by image analysis on histological sections. Microsc Res Tech (1999) 45:303–12. doi:10.1002/(SICI)1097-0029(19990515/01)45:4/5<303::AID-JEMT14>3.3.CO;2-#

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Hutmacher DW. Scaffolds in tissue engineering bone and cartilage. Biomaterials (2000) 21:2529–43. doi:10.1016/S0142-9612(00)00121-6

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Wang X, Nyman JS, Dong X, Leng H, Reyes M. Fundamental biomechanics in bone tissue engineering. Synth Lect Tissue Eng (2010) 2:1–225. doi:10.2200/S00246ED1V01Y200912TIS004

CrossRef Full Text | Google Scholar