Preparation of a skin equivalent phantom with interior micron-scale vessel structures for optical imaging experiments

A popular alternative of preparing multilayer or microfluidic chip based phantoms could have helped to simulate the subsurface vascular network, but brought inevitable problems. In this work, we describe the preparation method of a single layer skin equivalent tissue phantom containing interior vessel channels, which mimick the superficial microvascular structure. The fabrication method does not disturb the optical properties of the turbiding matrix material. The diameter of the channels reaches a value of 50 μm. The size, as well as the geometry of the generated vessel structures are investigated by using the SD-OCT system. Our preliminary results confirm that fabrication of such a phantom is achievable and reproducible. Prospectively, this phantom is used to calibrate the optical angiographic imaging approaches. © 2014 Optical Society of America OCIS codes: (170.0170) Medical optics and biotechnology (110.7050) Turbid media; (170.3880) Medical and biological imaging; (160.4760) Optical properties; (350.0350) Other areas of optics; References and links 1. M. Firbank and D. T. Delpy, “A design for a stable and reproducible phantom for use in near infra-red imaging and spectroscopy,” Phys. Med. Biol. 38(6), 847 (1993) 2. B. W. Pogue and M. S. Patterson, “Review of tissue simulating phantoms for optical spectroscopy, imaging and dosimetry,” J. Biomed. Opt. 11(4), 041102 (2006) 3. M. Lualdi, A. Colombo, B. Farina, S. Tomatis, and R. Marchesini, “A phantom with tissue-like optical properties in the visible and near infrared for use in photomedicine,” Laser. Surg. Med. 39(3), 237–243 (2001) 4. D. D. Roystonand, R. S. Poston, and S. A. Prahl, “Optical properties of scattering and absorbing materials used in the development of optical phantoms at 1064 nm,” J. Biomed. Opt. 1(1), 110–116 (1996) 5. G. Lamouche, B. F. Kennedy, K. M. Kennedy, C. E. Bisaillon, A. Curatolo, G. Campbell, V. Pazos, and D. D. Sampson, “Review of tissue simulating phantoms with controllable optical, mechanical and structural properties for use in optical coherence tomography,” Biomed. Opt. Express 3(6), 1381–1398 (2012) 6. R. J. Cooper, R. Eames, J. Brunker, L. C. Enfield, A. P. Gibson, and J. C. Hebden, “A tissue equivalent phantom for simultaneous near-infrared optical tomography and EEG,” Biomed. Opt. Express 1(2), 425–430 (2010) 7. A. V. Mudaliar, Development of a Phantom Tissue for Blood Perfusion Measurement and Noninvasive Blood Perfusion Estimation in Living Tissue (Virginia Polytechnic Institute and State University, 2007) #214728 $15.00 USD Received 25 Jun 2014; revised 10 Aug 2014; accepted 15 Aug 2014; published 22 Aug 2014 (C) 2014 OSA 1 September 2014 | Vol. 5, No. 9 | DOI:10.1364/BOE.5.003140 | BIOMEDICAL OPTICS EXPRESS 3140 8. C. L. De Korte, E. I. Cespedes, A. F. W. Van der Steen, B. Norder, and K. te Nijenhuis, “Elastic and acoustic properties of vessel mimicking material for elasticity imaging,” Ultrasonic imaging 19(2), 112–126 (1997) 9. W. Dabrowski, J. Dunmore-Buyze, R. N. Rankin, D. W. Holdsworth, and A. Fenster, “A real vessel phantom for imaging experimentation,” Med. Phys. 24, 687 (1997) 10. A. Surowiec, P. N. Shrivastava, M. Astrahan, and Z. Petrovich, “Utilization of a multilayer polyacrylamide phantom for evaluation of hyperthermia applicators,” Int. J. Hyperther. 8(6), 795–807 (1992) 11. R. B. Saager, C. Kondru, A. U. Kendrew, K. Sry, F. Ayers, and A. J. Durkin, “Multilayer silicone phantoms for the evaluation of quantitative optical techniques in skin imaging,” Proc. SPIE 7567, 756706 (2010) 12. A. V. Bykovand, A. P. Popov, A. V. Priezzhev, and R. Myllylä, “Multilayer tissue phantoms with embedded capillary system for OCT and DOCT imaging,” Proc. SPIE 8091, 80911R (2012) 13. A. V. Bykov, A. P. Popov, M. Kinnunen, T. Prykäri, A. V. Priezzhev, and R. Myllylä, “Skin phantoms with realistic vessel structure for OCT measurements,” Proc. SPIE 7376, 73760F (2010) 14. V. V. Tuchin, A. N. Bashkatov, E. A. Genina, V. I. Kochubey , V. V. Lychagov, S. A. Portnov, N. A. Trunina, D. R. Miller, S. Cho, B. Shim, M. Kim, J. Oh, H. Eum, Y. Ku, D. Kim, and Y. Yang, “Finger tissue model and blood perfused skin tissue phantom,” Proc. SPIE 7898, 78980Z (2011) 15. J. T. Wang, J. C. James, C. P. Liang, N. Woolsey, J. C. Ramella-Roman, Y. Chen, and T. J. Pfefer, “Threedimensional printing of tissue phantoms for biophotonic imaging,” Opt. Lett. 39(10), 3010–3013 (2014) 16. L. Luu, P. A. Roman, S. A. Mathews, and J. C. Ramella-Roman, “Microfluidics based phantoms of superficial vascular network,” Biomed. Opt. Express 3(6), 1350–1364 (2012) 17. N. T. Nguyen and S. T. Wereley, Fundamentals and Applications of Microfluidics (Artech House, 2002) 18. A. N. Bashkatov, E. A. Genina, and V. V. Tuchin, “Optical properties of skin, subcutaneous, and muscle tissues: a review,” Journal of Innovative Optical Health Sciences, 4(1), 9–38 (2011) 19. T. J. Farrell, M. S. Patterson, and B. Wilson, “A diffusion theory model of spatially resolved, steady-state diffuse reflectance for the noninvasive determination of tissue optical properties in vivo,” Med. Phys. 19, 879 (1992) 20. S. J. Madsen, M. S. Patterson, and B. C. Wilson, “The use of India ink as an optical absorber in tissue-simulating phantoms,” Phys. Med. Biol. 37(4), 985 (1992) 21. S. A. Prahl, M. J. C. van Gemert, and A. J. Welch, “Determining the optical properties of turbid mediaby using the adding–doubling method,” Appl. Opt. 32(4), 559–568 (1993) 22. V. V. Tuchin, Tissue Optics: Light Scattering Methods and Instruments for Medical Diagnosis (SPIE press Bellingham, 2012) 23. A. E. Cerussi, R. Warren, B. Hill, D. Roblyer, A. Leproux, A. F. Durkin, T. D. O ’ Sullivan, S. Keene, H. Haghany, T. Quang, W. M. Mantulin, and B. J. Tromberg, “Tissue phantoms in multicenter clinical trials for diffuse optical technologies,” Biomed. Opt. Express 3(5), 966–971 (2012)


Introduction
Optical imaging studies are widely performed in clinical or experimental trials to reconstruct human organs, such as brain, liver, breast and skin.Inversely, these organs are most commonly mimicked to calibrate the optical imaging approaches and to determine the adaptivity of such devices on biological objects [1,2].For this purpose, these mimicked body parts should be modeled with equivalent tissue phantoms.Phantoms for optical studies have already fulfilled the requirements of stability, biological compatibility, and replicated the representative values of specific optical properties of the real tissue [1][2][3][4][5][6].These requirements have now swayed into generating realistic physiological structures, that permits an increasing complexity of optical imaging studies on tissue phantoms [7][8][9].
Mapping the microvascular network under skin promotes the understanding of some specific types of vascular dermatosis or the local angiogenic activities during the tumor progression in skin.To verify the feasibility of applying optical imaging approaches on recovering the superficial microvasculature, skin phantoms with subsurface structures should be particularly prepared for this scenario.In view of skin equivalent phantoms, many researches have been conducted onto the concept of multilayers that imitate the layers of the cutis: epidermis, dermis.Multilayer structure favors to assess some imaging modalities and their associated algorithms [9][10][11][12].It is even practical to emboss vessel structures in multilayer phantoms to include biological / artificial solution [8,13].As a typical example, a blood perfused finger phantom was prepared by folding glass tubes between layers to create an array of vessel structures [14].J. T. Wang et al. presented an additive production process to print a vessel pattern inside a phantom in their pioneer work [15].These embossed structures lack of an exact dimensional analogy to the real microvascular according to the morphology of skin [13][14][15].Fast prototyping of pattern on microfluidic chip is becoming a popular option to simulate the nature of superficial microvasculature and achieve a dimensional similarity [16].
However, experiments on multi-layer phantom [14], 3-D printed phantom [15] or microfluidic chip [16], are not easily duplicated in a few steps or without particular equipments.It means, they have an innate shortcoming of cumbersome to mass manufacture.Traditional fabrication of microfluidic chips using soft lithography, hot embossing could be even more complicated and costly [17].
It is crucial to decrease the diameter of the embedded vessel structures to a value resembling a real microvasculature.Another challenge is, that the optical properties of this phantom replicate those of skin tissue and should not be significantly changed by the fabrication process.This work pertains to prepare a skin equivalent phantom with micron-scale interior hollow vessel channels to utilize biological / artificial inclusion.Ultimately, this phantom will be used to calibrate the optical angiographic approaches, like SD-OCT, diffuse optical imaging or hyperspectral imaging studied in our lab.

Designing and preparing the casting mold
The diameter of a real microvascular structure under skin ranges from 10 to 100 µm, the capillary network is mostly bedded in the first 200 µm of the skin subsurface.To simulate this, we intend to derive thin hollow vessel channels with a diameter of 50 µm at the depth of 150 to 200 µm from the surface.The entire progress is summarized as embedding-and-etching copper wires in the matrix material of the phantom.
An iron cylinder is firstly prepared as the casting mold to regulate the positioning of the vessel channels in phantom.We use a ps-laser (TBWP Fuego ps-laser, Time-Bandwidth Products AG, Switzerland) to drill holes through the wall of this mold.A laser beam spot is located on the defined positions by rotating the mold with a chuck.The focused laser beam scans over the wall in a circular hatched pattern to ablate the metal.By this method, a minimal drilling threshold of 50 µm in diameter as well as a sharp edge around the holes is produced.Conventional drilling machine, like an electrical discharging drilling machine, could also be used in this case.No specific means of machining are necessary.Holes with a diameter of 100 µm are stamped through the wall to fix the copper wires (φ = 50 µm, Conrad Electronic GmbH, Germany).Likewise, holes with a diameter of 400 µm are drilled to fix the 21 gauge cannulas as the mounting brackets of the copper wires.Altogether 8 holes (each 4 for every 10 degrees symmetrically on the opposite of the cylinder mold) for inserting copper wires and 4 holes (respectively at 60°, 140°, 240°, 320°) for inserting canullas are proposed.In this design as shown in Fig. 1(a), 2 mounting brackets perpendicularly cross over 4 copper wires.Copper wires and cannulas are plugged into the corresponding drilling holes afterwards (see Fig. 1(b)).Hollow vessel channels should be left in the phantom after the copper wires and canullas are eliminated from the matrix material.

Preparing the matrix material
For simulated studies on skin via optical imaging approaches, especially with diffuse optical imaging (DOI), the phantom should allow for similar absorption coefficient μ a and reduced scattering coefficient μ s to a real skin tissue in the wavelength dependence.Therefore, we choose the data from ex-vivo caucasian skin tissues [18,19] as the reference.The tissue phantom is constructed with castable aliphatic polyurethanes (WC-781, BJB Entprise Co., US).To replicate the optical properties of the caucasian skin, 0.7 g of titanium oxide powder (Sigma-Aldrich GmbH, Germany) and 3.3 µL of india ink (Pelikan GmbH, Germany) [4,18,20,21] are added into 100 mL of polyurethane component.The additives are dispersed by 20 min stirring, followed with an ultrasound bath (Elmasonic P, Elma Hans Schmidbauer GmbH, Germany) under 60 • C. We also use a vacuum autoclave to extract the air bubbles from the viscous mixture before casting.Casting is excuted at 30 • C in a convection chamber (BINDER GmbH, Germany), where a dehydrated ambience helps to reduce the reactive air bubbles in polyurethane.Solidified phantom is de-molded after 48 h, when most mechanical properties ( e.g.shore hardness = 82 D) are reached.

Etching the retained copper wires
Ferric chloride is used to eliminate the folded copper wires from the fabricated phantom.To avoid the swelling of polyurethane, the concentration of FeCl 3 ( Sigma-Aldrich Co., US) is restricted at 50 g/L, without adding acid adjuvant.We precipitate the de-molded phantom into the ferric chloride solution, keep the ambient temperature at 45 • C and stir the solution every 1 hour to avoid the material stagnation on the local corrosion knot.Ferric chloride solution perfuses into the channels and dissolves the wires without disrupting the wall of the channels, thus produces a finished phantom.
However, the etching progress might alterate the optical properties of the phantom.Possible attenuation mainly stems from the dying of titanium oxide particles on the surface during their contact with the ferric chloride solution.Due to this, we use UV-IR spectrophotometry (UV-3600 seires, SHIMADZU Co., Japan) to characterize the optical properties through the matrix material of the phantom before and after etching, as well as to verify the chemical resistance of polyurethane to ferric chloride.To do this, spectra of absorbance A b , diffuse reflectance R d , diffuse transmittance T d and total transmittance T t are measured.We also investigate the geometry of the vessel structures and evaluate the uniformity of the diameter along the vessel channels by using a SD-OCT (Telesto-II-SP1, Thorlabs GmbH, Germany).

Results and discussion
We prepare a phantom as the object sample containing 4 separate rows of vessel channels with replicated optical properties (see Fig. 2(a)).Besides, a simplified reference phantom with 2 rows of vessel channels is constructed without inserting scatters and absorbers.This transparent phantom is individually prepared to monitor the dissolving of the copper wires in the matrix before and after etching show the same trace with an overlaying error bar in the full wavelength range.The standard deviations of these spectra do not exceed ± 1 %.Values of R d alterate a bit after etching, with a drop of around 4 %.This drop is mostly caused by the induced change of the surface status, as the surface layer of the phantom slab is dyed during etching.
To verify the reproduciblity of these results, same experiments are respectively repeated on 4 phantom slabs with different content of titanium oxide.Table 1 summerizes the attenuation rate of the T t , T d , A b and R d before and after etching, with regarding to the weight of titanium oxide in 100 mL of polyurethane.According to Tbl. 1, a higher content of TiO 2 could have increased the change of optical properties, espcially diffuse reflectance.Nevertheless, This change does not exceed 5 %, which is rather slight and therefore acceptable.By comparing these spectra, it can be verified that the ferric chloride solution does not reveal any significant alteration in the optical properties of the turbid phantom slab.

Optical properties compared with ex-vivo skin tissue
The reference optical properties are reported elsewhere in [18,19], the paradigm of absorption coefficient μ a and reduced scattering coefficient μ s are attained from the whole skin slabs out of 21 caucasian human objects at wavelengths of 500, 600 and 700 nm and converted using a single intergration sphere and inverse adding doubling [21].Likewise, we measure the same coefficients from the tissue phantom at the same wavelengths (see Tbl. 2).The absorption coefficient μ a is close to that of the real ex-vivo skin tissue.μ s is slightly overestimated.This could be further modified by reducing the amount of the titanium oxide particles.The refractive index of phantom is lower than the real value at 500 nm [22].Nevertheless, this is not disturbing, when μ a and μ s approach the values of the ex-vivo tissue [23].Summarizing from the previous results, we have successfully fabricated the object phantom with 4 seperate rows of vessels and a transparent phantom to ensure the feasibility of including liquid into the hollow vessels.The object sample replicates the optical properties of a caucasian skin.Comparing the registered spectra from before and after etching, it is confirmed that the chemical resistance of polyurethane to ferric chloride at 45 • C is a strong hold, the etching progress does not lead to any significant attenuation in the optical properties of the phantom.Such a slight attenuation on the surface of the phantom slab can even be neglected, as if the dyed section is ground off beforehand.Limit stays, that the reducing the attenation on the wall of the generated vessel channels is beyond our ability, although such a slight drop of the diffuse reflectance at 4 -5 % is acceptable.
In contrast, direct engraving on the matrix material via thermal, mechanical or photonic fabrication brings inevitable problems.For instance in the fast prototyping method [16], discussion about the changes of optical properties in the heat affected zone (HAZ), induced by laser micromachining, is missing.Moreover, it is critical to package the microfluidic chip with an adhesive material, whose optical properties attenuate over time.Other engraving methods could mostly disturb the optical properties of the turbiding matrix material [17].Although the nature of the mutilayers of skin cutis can not be simulated, like described in [2,10,14], the diameter of the generated vessels, that we pursued, is more close to that of a real vascular compared to the results from Refs.[12][13][14][15].Further discussion goes onto investigating the geometry of the generated vessel structures.

Investigating the geometry of vessel channels
The diameter and depth of the vessel channels are separately measured on the finished tissue phantom over around the mounting bracket using the optical microscopy (OM) and optical coherence tomography (OCT).The OM graphic (see Fig. 4(a)) demonstrates a single vessel channel in the transparent tissue phantom.Regarding a volume shrinkage ratio of the aliphatic polyurethane, the approximately achievable diameter of the vessel channels should be around 55 µm.The measured value is 59 µm on average, which shows an excellent agreement to the estimated value.
In order to validate the vessel geometry more directly, we refer to OCT for a reconstruction of the vessel structures entirely, or partially around a certain area of interest.Using OCT, we can recover the cross section passing along a vessel channel (see Fig. 4(b)), along the mounting bracket across the vessel channels (see Fig. 4(c)) through a B-scan, and a segmented vessel structures in a 3-D volume view (see Fig. 4(d)).Similarly, the measured value of vessel diameter is about 52 µm on average.The vessels are uniform in diameter along the entire structure, with a tolerance of ± 3.3 µm.Conclusion could be drawn, that the geometric parameters of the vessel channels conform the original intension of design.
The diameter of the generated vessels is decided by the master mold of the inserted copper wires.However, the core diameter of the wire with a value smaller than 50 µm is not off-theshelf.Further decrease of vessel dimension could be implemented by using an ultra-thin-wire.In this case, we attempt to produce a thinner wire by eroding the φ 50 µm copper wires to φ 20 µm with electrolysis.Experiments are repeated with a low erosion rate to get a core diameter close to the aimed value.Geometry of the generated vessel are observed after the eroded wires are dissolved from a transparent phantom.The results are summarized in Tbl. 3. We demonstrate that the embedding-and-etching method can create a hollow vessel with a diameter down to φ = 20 µm.Although we doubt the dimensional uniformity of such a vessel channel, the value of the diameter still reaches the unique feature value, recorded elsewhere in [16].Uniformity could be simply improved by using ultra-thin-wires.

Deformation of the vessel channels
It is apparently indicated in the OCT B-scan image (see Fig. 5(a)), that the vessel channels may deviate from the supposed normal latitude.In this subfigure, the distance from the upper outline of the generated vessel channels to the surface ranges from φ 120 µm to φ 210 µm.This occurs because of the thermometrical force inside the polyurethane in tangential direction, induced by rapid cooling during solidification.It is even more visibly obvious from the OCT top view (see Fig. 5(b)), that the channels bend along the structure.Such a deviation could even reach dozens of micrometers.This phenomena, although random and even inevitable, could be reduced by cooling the polyurethane in a cooling rate at aroud 1 • C / min until 30 • C.

Improving homogeneity of scatter
The homogeneity of scatters influences the light propagation in turbid media.The agglomeration of the titanium oxide particles may exist until the matrix material is solidified, thus forms a scattering cluster.In our case, the diameter of the generated vessel channels encloses to the feret diameter of the scattering clusters.Scattering clusters might be misidentified, or fracture the real target of interest, namely the vessels.Such scattering clusters must be, and could be reduced by a longer ultra-sound bath, which smashs and disperses the powder agglomeration.To do this, a 2 h ultra-sound-bath is executed to substitute the 20 min stirring.To check the consequence of this homogenization, scatter homogeneity is evaluated from the OCT B-scan image by using an image post-processing software (ImageJ, National Institutes of Health, US) respectively before and after homogenization.The scattering clusters are thresholded, and morphologically opened during the post processing of the OCT images (see Fig. 6 from (a) to (c), or from (b) to (d)).Those scattering clusters, whose size exceed a certain value, are contrasted against the background.By comparing Fig. 6(d) to (c), we find that the speckeling spots are mostly vanishing from Fig. 6(d).The claim is, that the amount and the size of scattering clusters decrease after homogenization.Further problems like the residual air bubbles could be solved by inserting anti-foam components.In this way, quality of the tissue phantom could be ensured.

Conclusion
In the framework of preparing a skin equivalent tissue phantom containing thin interior vessel channels, prototype with 4 separate rows of channels are constructed.Optical properties similar to real ex-vivo caucasian skin tissues are replicated and not attenuated by the embedding-andetching progress.With the OCT studies, we assure the feasibility of creating a vessel structure with a diameter of 52 µm at a depth of 150 to 200 µm from the phantom surface.Our fabrication method is easily reproduced, and does not require any specific equipments.Future work will be conducted onto upgrading the geometric complexity of the vascular pattern and varying the optical properties for a phantom databank.The fabrication of this phantom also provides a fresh frame of mind to prepare a multi-layer tissue phantom, with an overlaying microvascular network, and meanwhile, without using any adhesive sealing material.Our upcoming work will be to recover this superficial microvascular pattern with a customized diffuse optical imaging apporach.

Fig. 1 .
Fig. 1.(a) Schematic of the cylinder casting mold design for regulating 4 rows of vessel channels in the tissue phantom and (b) image of a finished casting mold with 4 rows of copper wires and 21 gauge cannulas inserted through the drilling holes

Fig. 2 .Figure 3 Fig. 3 .
Fig. 2. Picture of (a) a finished tissue phantom with 4 rows of vessel channels perpendicularly crossing over 2 mounting brackets and (b) a finished transparent tissue phantom with 2 rows of channels, they are included with india ink to highlight the hollow vessel channels against the matrix material

Fig. 4 .
Fig. 4. (a) OM graphic of a generated vessel channel dyed with india ink in a transparent tissue phantom (b) OCT B-scan image of the cross section passing along a vessel channel (c) OCT B-scan image of the cross section passing along a mounting bracket across the the vessel channels (d) OCT volume view of the segmentated vessel structures, including vessel channels and mounting brackets

Fig. 5 .
Fig. 5. (a) OCT image of the cross section across over 4 vessel channel, showing that they are not located on the same latitude of around 160 µm to the surface (b) OCT top view of vessel channels in tissue phantom : vessel channels bend in tangential direction along the structure due to the inner stress induced during solidification

Fig. 6 .
Fig. 6.Origin OCT B-Scan image of the cross section of (a) a typically inhomogeneous tissue phantom (b) homogenous tissue phantom, subfigures (c) and (d) demonstrate the segmentated scattering clusters corresponding to (a) and (b)

Table 1 .
Attenuation in the values of T t , T d , A b and R d with regarding to different content of titanium oxide in 100 mL of polyurethane

Table 2 .
Optical properties of the ex-vivo caucasian skin tissues and tissue phantom calcu-

Table 3 .
Estimated and real values of diameter of the eroded copper wires and generated vessel channels in a transparent phantom slab