The nondestructive characterization of engineered tissues has generated enormous interest recently. Frequently, tissue engineers are interested in monitoring the local microenvironment in an engineered tissue and assessing parameters, such as cell number, cell function, protein expression, and mechanical properties. Most current techniques that characterize these parameters are either destructive, require adding (often expensive) labels to the sample, or making macroscopic measurements with no possibility of assessing local variations.1, 2, 3, 4 On the microscopic level, engineered tissue constructs often have local heterogeneities, and the patency of implanted constructs can depend on assessing these local variations.
In general, optical methods have shown great promise for nondestructively characterizing engineered tissues. One particularly promising method is optical coherence tomography (OCT), which offers several unique technological advantages for characterizing engineered tissues5. Commercial OCT systems can achieve an axial spatial resolution of and penetrate up to in highly scattering media, imaging at video-rate speeds. As a fiber-based technology, OCT can be coupled into flexible handheld probes, allowing them to be easily incorporated into incubators.
Several recent studies have used OCT to characterize engineered tissues. Initially, Tan used OCT to visualize how fibroblasts change the architecture of a 3D chitosan scaffold,6 and then to monitor cell migration, proliferation, and detachment off a calcium phosphate scaffold by trypsinization.7 Rey visualized chemotaxis-induced cell migration in 3-D agarose substrates.8 OCT has also been used to visualize how cells change the porosity in poly(l-lactic acid) scaffolds9 and chitosan scaffolds.10
These papers clearly demonstrate the ability of OCT to characterize features such as cells (or a cluster of cells) on a weakly scattering background of an imaged tissue construct. However, as the tissue construct develops, gel compaction, cell proliferation, and matrix remodeling change the size and density of scatterers in the sample. At this stage, the construct becomes a highly scattering medium, in which the cells can no longer be resolved. Properly characterizing mature tissue constructs just prior to implantation will be important for ensuring successful outcomes, but given the highly scattering nature of mature constructs, making such assessments is quite challenging.
Moreover, the bulk of the work in characterizing engineered tissue constructs using OCT has been qualitative image analysis, based primarily on identifying features of interest (i.e., a cell) by a subjective operator. Given the qualitative nature of these methods, systematically surveying large volumes of data ( ) imaged at micron resolution is both cumbersome and impractical. However, OCT data contain quantitative information on the optical scattering properties of the imaged sample. The optical properties of interest here are the scattering coefficient [ ] and anisotropy factor [dimensionless],11 which provide a measure of the density and size of scattering particles in the medium, respectively. With proper theoretical modeling, it is possible to describe the depth-dependent OCT signal as a function of and , as well as the system parameters.12 By fitting a depth-dependent OCT signal that is representative of local nanoliter volumes in the sample to a theoretical model, it is possible to determine the local scattering properties of the sample. This approach was used to evaluate the optical properties of one of the most common types of engineered tissues, collagen gels, and characterize their development when seeded with cells.
Collagen gels are one of the simplest types of engineered tissues, in which cells are suspended in a 3-D collagen I matrix. Over time, the cells pull on the collagen fibrils causing gel compaction, wherein the gel becomes smaller, with decreased water content, and visibly more opaque. Some mesenchymal cells, such as fibroblasts and smooth muscle cells (SMCs), also contract and remodel the collagen matrix.13 An example of SMC-seeded collagen gels at and five days is shown in Fig. 1 . Note that the gels are large and nearly transparent, while the five-day gels have contracted and have a cloudy appearance.
In this paper, collagen gels are imaged with different cell-seeding densities over a five-day period and the data are used to measure their optical properties. Over time, the reflectivity of the gels increased while the attenuation did not change, which degraded the ability to visualize structures such as cells. This reflectivity increase was due to a decrease in anisotropy , which caused increased backscatter of light. No such changes were observed in acellular gels. In each set of gels imaged at five days, reflectivity values followed a bimodal distribution, with the major population of gel sites showing a doubling of reflectivity and a sizeable subpopulation showing a roughly tenfold increase in reflectivity relative to the first . Time dynamics of the gels showed that these changes in optical properties were gradual. A careful analysis of the optical properties measured during the earliest time point showed that a correlation existed between cell density and , but that cell density and were uncorrelated. Finally, mechanisms that could be responsible for the observed changes in optical properties are discussed.
The theory of light propagation in low-coherence imaging systems used to process the data in this paper was derived from the inverse Monte Carlo method.12 For a homogeneous turbid medium characterized by a scattering coefficient , an anisotropy factor , an absorption coefficient , and a refractive index , the depth-dependent OCT reflectance signal can be described as an exponential decayis the fraction of light delivered to the focus that is backscattered by the focal volume of tissue within the solid angle of collection of the microscope’s objective lens, and the attenuation describes the signal attenuation as light propagates to and from the focal volume. Specifically, is the axial resolution of the imaging system and is a geometry factor accounting for the extra pathlength caused by off-axis propagation during delivery and collection by the objective lens. Analysis has shown that for an objective lens with numerical aperture (NA), the parameter can be approximated by . For low-NA imaging systems frequently used in OCT, is close to 1.0 (the exact value of for high NA systems has been discussed12). The product is the number of scattering events that occur within the focal volume, and the parameter represents the fraction of light scattered within the focal volume that is backscattered into solid angle of collection of the objective lens. For a medium with scattering phase function , is given by
In our analysis, a Henyey–Greenstein phase function14 is assumed. One of the terms in Eq. 3 is , , that decreases the effectiveness of scattering when the scattering is forward directed such that photons, despite multiple scattering, still manage to reach the focus and cannot be distinguished from the unscattered light. For highly anisotropic media, has the effect of reducing the observed attenuation in the OCT signal. The function can be approximated by a numerical expression determined from Monte Carlo simulations12 The coherence gating in OCT improves on RCM by rejecting photons whose path length in the tissue does not match the round-trip path length to/from the focus, which yields about a 1000-fold noise reduction that enables deeper imaging.12 Preliminary data show that this model is valid approximately up to a (one-way) optical depth of , before multiple scattering begins to considerably interfere with the signal.
For this model to be valid, it is necessary to convert the interferometric OCT envelope (in arbitrary units) into dimensionless units of reflectance , i.e., the fraction of delivered light that is collected and detected by the microscope. To do this, the signal is calibrated against the peak signal from a known interface, such as an oil-glass interface, or against the reflectivity of an optical scattering phantom with known optical properties. This calibration is crucial in order to map the observations and into the optical properties and .
Materials and Methods
All cell culture materials were purchased from GIBCO (Invitrogen, Carlsbad, California) unless otherwise noted. The experiments used primary SMCs isolated from a single baboon carotid artery that was chemically digested. Passage 9 SMCs were used in preparing collagen gels. The cells were fed with SMC growth medium (SGM), consisting of minimum essential medium (MEM) supplemented with 10% fetal bovine serum, 1% L-glutamine, and 1% penicillin-streptomycin.
Collagen Gel Preparation
All collagen gels were prepared using the same basic protocol. Soluble calf skin collagen (part no. 150026, MP Biomedicals, Solon, Ohio) was dissolved in acetic acid at a concentration of . In preparing the collagen gels, eight parts collagen were mixed with two parts 5× MEM and one part NaOH (Sigma, Sigma-Aldrich, St. Louis, Missouri), with NaOH titrated until the pH was neutral as indicated by a color change in phenol red in the solution. Immediately thereafter, five parts SGM containing the predetermined number of cells were added to the mixture and aliquoted as samples into a 12-well plate (Corning, Lowell, Massachusetts). The weakly scattering suspension thermally gelled for at . Thereafter, SGM was added to feed each gel. Overall, 24 gels were prepared at each cell concentration ( , , and ), which is denoted below as a “set” of collagen gels. Additionally, a set of acellular gels was also prepared and imaged. Acellular collagen gels were prepared in an identical manner to the SMC collagen gels, except that cells were absent from the mixed solutions.
Gels that were cultured for longer than were detached from the well at with a spatula and allowed to further compact. These gels were fed SGM every two days. At each imaging time point, three gels were removed from their wells and imaged by OCT. After OCT imaging, each sample was fixed for subsequent histological examination.
In addition to the four sets of collagen gels described above, optical property data acquired with a different OCT system are also presented. The two other sets of gels (one seeded with SMCs and one acellular) were imaged with a different OCT system at the and five-day time points. There were minor differences between these two sets and the four above, namely, the NA was 0.08 and the initial collagen density was . The OCT images, gel compaction data, and histology images for these two sets were nearly identical to the data presented here and are thus not shown. Altogether, imaging of the six sets of collagen gels was done over a span of .
A time-domain OCT imaging system was setup as previously described.15 Briefly, a broadband light source operating at (B&W Tek, Newark, Delaware) with a full width at half maximum (FWHM) bandwidth was coupled into a fiber-based Michelson interferometer. A rapid scanning optical delay line16 varied the delay in the reference arm. A 10× microscope objective (Newport, Irvine, California) focused the light on the sample. Both scanning mirrors in the reference and sample arms ( scanning) were from GSI Lumonics (Bedford, Massachusetts). The interference signal between the backreflected light from the sample and reference arms was measured using an InGaAs photodiode (New Focus, Santa Clara, California). The OCT system was controlled through a data acquisition card (National Instruments, Austin, Texas) using software written in Matlab (MathWorks, Natick, Massachusetts). As the system scanned axially, the interference of the sample and reference arms yielded an oscillating signal. The Hilbert transform was performed in software postacquisition using Matlab to specify the envelope of this oscillating signal, which was proportional to the reflected electric field. Squaring this envelope yielded the signal proportional to reflected intensity.
During OCT imaging, the sample was mounted on a glass slide and immersed in MEM to maintain hydration throughout imaging. Each gel was imaged at five sites, with the sample repeatedly imaged 10–12 times at each imaging site. The sample orientation was randomized between imaging sites. In the 10–12 images acquired at each site, the translation stage on which the sample was mounted was manually moved toward the objective lens by to a position . An axial OCT scan versus depth was obtained at each lens position, yielding a signal . In separate experiments using an oil-glass interface that could be positioned at different depths, a focus function was specified, , which accounted for the variation in response as the position of a structure of interest deviated from the central focus position . showed that the portions of each image in the range experienced only minor defocusing to roughly 81% of the maximum value. Moreover, was chosen as the focus shift distance since is of the depth of focus in our OCT system . The portions of each image in the range provided the strongest signals and were stitched together to yield a final image. Defocusing effects were negligible in comparison to the speckle fluctuations in the signal, which spanned an order of magnitude. In this manner, a quasi-focus tracking image was prepared, appropriate for the theoretical OCT model.
OCT Image Processing Algorithm
The OCT data from all samples were processed in a similar manner. A single, Matlab-based algorithm for evaluation of optical properties from the data was organized in the following separate modular steps:
1. The OCT data at each imaging site were read into Matlab and stitched together to form a quasi-focus tracked image. Using user-defined landmarks as points of reference (typically the air-sample interface or the corners of the sample), images taken at a single site were stitched together.
2. The data from a stitched image were then calibrated against the signal from a known reference. Initially, the data were calibrated against a mineral oil-glass interface . However, polystyrene sphere-agarose gel phantom data acquired during the same period of collagen gel imaging revealed that sphere phantoms are a more stable calibration method, and thus all the relevant collagen gel images were cross-calibrated against sphere data from 210 regions of interest (ROIs) from 30 imaging sites in six phantoms.
3. To align OCT images, the most superficial pixel of the gel was found at each lateral position, and the image was aligned by shifting each A-scan axially. For samples with the top surface located at the air-MEM interface, this edge was found by fitting a parabola to the interface pixels. For samples submerged in MEM with an irregular surface, the user manually drew a line along the gel surface using the computer mouse. Once the sample surface was identified, the A-scans were shifted accordingly.
4. The user subdivided each calibrated, stitched, aligned image into five ROIs along the direction for fitting. The dimensions of each ROI were approximately . The average signal over at each position yielded .
5. For each ROI, versus was fit to Eq. 1 using a linear fit. A weighted-least-squares fitting routine was used, with the weight at each pixel being inversely proportional to the variance of the signal through all the pixels at that depth. In addition to the fitted parameters and , and their respective confidence intervals, goodness-of-fit statistics—the mean pixel error and the adjusted regression coefficient —were documented for each ROI.
Collagen gel compaction was monitored in the gels at eight time points and expressed as a percent of the initial gel volume . The volume of the gels was indirectly measured via the volume of the media in the well that could be aspirated after compaction, i.e., gel minus the aspirated volume, with the coming from a gel fed with SGM. The gel compaction was measured on SMC gels with three cell seeding densities and an initial collagen density of . Calibration experiments showed that this method overestimated the collagen gel volume by .
After imaging, samples were fixed overnight in 10% neutral buffered formalin and stored in 70% ethanol at . The samples were embedded in paraffin, and sections were cut and stained with hematoxylin and eosin (H&E) on glass slides. Sections of the selected samples were cut from positions starting at the radius of the sample toward the center, in increments. The histological sections were cut so they would expose a cross section of the gels, but they did not exactly match the imaging planes in the OCT data since the gel orientation was randomized between the imaging sites. The histological sections were imaged by phase-contrast microscopy (Optiphot, Nikon, Tokyo, Japan) at 20×.
Data Analysis and Statistics
Each linear fit of data at each ROI resulted in a pair. The data were plotted as a scatter plot against a grid of iso- and iso- contours generated by the model of Eqs. 1, 2, 3. The contours map experiment to theory. Alternatively, the data were assessed by looking at the distributions of fitted values. For both and , the probability density functions (PDFs), and , were computed by binning and values from an individual group into histograms made of Gaussian shaped kernels. Treating and as statistically independent random variables, a joint PDF was computed using . Finally, the peaks of the joint PDFs were mapped back to with the a priori knowledge of the basic system parameters and NA using Eq. 2, 3.
To determine the accuracy and precision of the algorithm used in this paper, preliminary experiments were conducted to statistically characterize the fitting results. Specifically, intra-ROI, inter-ROI, inter-site, and inter-sample coefficients of variation were calculated for sphere phantoms, collagen gels, and five-day collagen gels (Table 1 ). Data in Table 1 are expressed as deviation over samples, and are in units of percent. The intra-ROI variability describes the reproducibility of 10 measurements of and on the same ROI. The inter-ROI variability describes the variation of the mean and of each ROI between five ROIs on the same site on a gel. The inter-site variability describes the variation of the mean data of five ROIs on a site between five sites on the same phantom/collagen gel. The inter-sample variability describes the variation of the mean data of five sites on one phantom between three different phantoms (or collagen gels). Because of this averaging, the inter-site and inter-sample variability in Table 1 are lower than the inter-ROI measurements. The sphere data in Table 1 illustrate the accuracy and reproducibility of the measurements, the collagen gel data illustrate the inherent variations in the gel preparation, and the five-day collagen gel data illustrate the inherent variations in gel development. The resolution of the algorithm for determining optical properties was specified as the FWHM of the PDF for either or values. The FWHM resolution in was finite at . For , the resolution scaled with the fitted results, with the FWHM being of the mean. In terms of accuracy, the algorithm output was shown to overestimate the theoretically predicted values for both and by .
Statistics describing the accuracy and reproducibility of the algorithm. Listed are the coefficients of variation for intra-ROI, inter-ROI (intra-site), inter-site (intra-sample), and inter-sample measurement in both μ and ρ , which are expressed in percent as means±standard deviations over n samples. Data are given for sphere phantoms, as well as the key time points of 1 and 5d during collagen gel development.
|Sampletype||Sphere phantoms||Collagen gels, 1d||Collagen gels, 5d|
Outliers that are > 3 standard deviations from the mean (corresponded to <4% of data) were excluded.
Ability to Visualize Cells in OCT Images of Collagen Gels Decreases with Time
Representative images of collagen gels with different cell seeding densities (acellular, , , and ) at , , and five days are shown in Figs. 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 , with corresponding histology for the SMC gels in Figs. 2, 2, 2, 2, 2, 2, 2, 2, 2. The bright line on the top of some of the OCT images represents the air-MEM interface. Overall, the OCT images appeared rather spatially uniform at any given time point, with the only visualized structures being bright spots that presumably represent individual cells or clusters of cells. Such bright spots on a weakly scattering background can be seen in the images of SMC gels [Figs. 2, 2, 2]. The density of these small “bright” spots increases with increasing cell seeding density. The bright spots were in size, comparable to SMCs. At [Figs. 2, 2, 2, 2] a few more of these bright spots can be seen in the gels [Fig. 2], and to a lesser extent in the and gels [Figs. 2 and 2–2(G)]. At [Figs. 2, 2, 2, 2], the images look different. In the gels [Fig. 2], there are a few bright spots that are difficult to resolve. The and gels [Figs. 2 and 2] show a superficial high signal region that decays. For both gel conditions, it is not possible to resolve any cells in the foreground, as the bright spots have blended into the signal rich background. Note that because of the high superficial signal and the attenuation, the images appeared featureless and uniform. These gels have entered a scattering regime in which qualitative image visualization does not yield much useful information and in which quantitative image processing techniques are needed. In all the acellular samples [Figs. 2, 2, 2], there were larger, modestly brighter, and darker regions in the background of all the images caused by heterogeneity in scattering from the collagen, which did not appear to change much over time. Histologically, the H&E sections in Figs. 2, 2, 2, 2, 2, 2, 2, 2, 2 show an increase in cell density over time for each of the three different seeding densities.
Collagen Gels Exhibit an Increase in Reflectivity/Decrease in Anisotropy over Time
The optical properties of the developing collagen gels from a single representative set are shown on a scatter plot in Fig. 3 . Also plotted on the scatter plot is a grid that maps theoretical values to the observed pairs. Here, the initial collagen density was and the initial cell density was . Figure 3 compares SMC and acellular gels at two key time points of one and five days, with each datum point representing a pair from a single ROI. In assessing the SMC and acellular gels at , the relative proximity of these two groups suggests that the scattering is dominated by the collagen fibrils, not the cells, which make sense given the small volume fraction the cells occupy. Also, note that the reflectivity in the SMC gels more than doubled from one to five days, whereas the corresponding acellular gels had only a modest increase in reflectivity. There was a large variance in attenuation values in all the samples. This feature has been observed consistently over multiple sets of collagen gels and results in poorer ability to resolve in collagen gels, in contrast to sphere phantoms (Table 1). Moreover, neither the SMC nor acellular gels experienced a major shift in attenuation over five days.
A summary of the optical property results from multiple sets of collagen gels is detailed in Table 2 , which lists the peak values and range ( or of peak value) for each and of the joint PDFs. These ranges provide an estimate of how and values were distributed around the joint PDF peaks. In comparing multiple sets of gels, the optical properties follow a general trend of starting out at low reflectivity and low attenuation at times , and gradually increasing reflectivity over five days, approximately twofold, while attenuation remained relatively constant. These observed ( , ) pairs corresponded to values that varied greatly, from , while values were consistently closely centered around 0.90 in the first . At five days, values of SMC gels consistently decreased 30–50%, while decreased to . Moreover, at five days, there was also a consistent subpopulation of data with higher reflectivity values, which were approximately tenfold the reflectivity observed in the first . Thus, at five days the reflectivity of collagen gels followed a bimodal distribution as can be seen in the reported peak ranges in Table 2. Note that the peak value is consistently closer to the lower boundary than the upper boundary, which suggests that although most of the data were centered around the peak, there was always a sizeable subpopulation with even higher values. In contrast to the SMC gels, acellular gels from set 4 and set 6 had very minor shifts toward either higher or lower reflectivity, respectively. Nonetheless, despite a large variance in , the increase in of collagen gels over time was both consistent and reproducible. Thus, although the spatial heterogeneity of collagen structures in the acellular gels remained constant over time, the subpopulation of higher values across the various sets of five-day SMC gels suggests that some spatial heterogeneity developed in the SMC gels over time.
Summary of optical properties from multiple sets of collagen gels, from the earliest time point ( 6h in sets 1–3, 24h in sets 4 and 5) to the final time point ( 5d for all sets). For both μ and ρ , we report both the peak of the joint PDF and the range over which the joint PDF > 1∕e* maximum. Square brackets denote concentration. Collagen density units are mg/ml; cell density units are cells/ml; μ and μs units are cm−1 ; ρ and g units are dimensionless.
|Set||No. gels||[collagen]t=0||[cell]t=0||Cell passage no.||μearly||μfinal||μearly pk. range||μfinal pk. range|
|61||6||2 or 53||0||NA||38.46||60.26||9.12–92.26||26.00–89.13|
|Set||pk. range||pk. range|
Note that sets 5 and 6 were imaged with a different OCT system.
Note 3 samples were imaged at eight different time points.
The initial collagen concentration in the gels imaged at 1d and 5d were 2 and 5mg∕mL , respectively.
Time Dynamics of Collagen Gels
Although Fig. 3 shows the overall shift in the local optical properties of the collagen gels, it is also worthwhile to visualize the time dynamics of the optical properties. However, because the shift in the data over time is slow and many data points lie on top of one another, the actual 2-D distribution of and values is presented as the joint PDF. A representative set of the calculated joint PDFs at , , three days, and five days are plotted superimposed on a grid in Figs. 4, 4, 4, 4 , showing the gradual increase in reflectivity over five days. Figures 4 and 4 shows that there are not many changes in the optical properties over the first . However, from one to five days [Figs. 4, 4, 4] the center of the PDF shifted to the right toward higher reflectivity, without changing much in attenuation. There was also a substantial subpopulation in the five-day time point that had an even higher reflectivity (lower anisotropy), similar to the other data sets (Table 2).
Additionally, the peaks of the joint PDFs shown in Figs. 4, 4, 4, 4 were mapped from space back to space. Plots of and versus time for eight time points spanning five days are displayed in Figs. 4 and 4. Error bars represent the width of the joint PDF at the iso- and iso- lines of the PDF peak. The versus time curve [Fig. 4] shows that over one to five days, there was a moderate but steady decrease in anisotropy. A similar moderate but steady decrease is also seen in the versus time plot [Fig. 4]. Note the nonlinear mapping of to that occurs at values , such that reflectivity can increase and attenuation can show no change while and both decrease.
Macroscopically, the SMC collagen gels showed a moderate but consistent rate of gel compaction, i.e., the rate of decrease in collagen gel volume as the cells pull on the fibrils was steady [Fig. 4]. Gel compaction was assessed at the eight imaging time points by quantifying the volume of media in the gel well just before imaging. The volume of gels from each of the three SMC groups shrunk to about 20–30% of their original volume. These data agreed with end point measurements by Cummings 17
Cell Density Correlates with Scattering Coefficient of Collagen Gels
When analyzing the optical property results, the data were closest to the zero time point and thus had minimal change in both the cells and the collagen. In plotting the peaks of the joint PDFs on a grid [Fig. 5 ], these pairs were colocated on an iso- line, and that their value increased with cell concentration. The corresponding and values for the four points were backtracked, as in Figs. 4 and 4, and plotted against the cell seeding density [Figs. 5 and 5]. Included in Figs. 5 and 5 are linear fits of both and against cell seeding density. The fit results showed that was highly correlated with cell concentration and increased linearly with cell density at a rate of plus a offset that represents scattering from collagen alone [Fig. 5]. Moreover, the measured from acellular collagen corresponds to a increase of of collagen density. A similarly high correlation between cell density and was obtained by correlating cell-like structures (i.e., “bright spots”) counted manually in the gel images (data not shown). In contrast to , remained flat relative to the concentration [Fig. 5] and was uncorrelated . Such an increase in without a change in is expected for a static system with an increasing scattering particle (i.e., cell) density. However, the imaged collagen gels are dynamic over five days, and no such clear cut correlations existed between either or and cell concentration at any of the later time points. This suggests that, over five days, the contraction process involved a fundamental change in the overall collagen gel architecture.
In this paper, OCT was used to image developing collagen gels at different time points and to measure the optical properties from the images by fitting the OCT signal to a theoretical model. The results showed that over time, there was a consistent twofold increase in reflectivity with little change in the attenuation in four separate sets of experiments with various cell seeding densities (Table 2, Figs. 3 and 4). Also, consistently present in each set of gels imaged at five days was a sizeable subpopulation of sites on the gel with a higher reflectivity, roughly tenfold greater than the first . These trends in attenuation and reflectivity correspond to a decrease in anisotropy and a moderate decrease in scattering coefficient , suggesting an overall decrease in the size of the scattering particles.
Several models have been proposed for measuring optical properties from OCT signals. The first model is the single-scattering (SS) model, in which the signal is assumed to follow an exponential decay , since for most tissues at near-infrared wavelengths. Despite similarities between the SS and models used here, our results illustrate several flaws in using the SS model for measuring optical properties. First, the model fits two parameters that map back to two unknowns, and , whereas the SS model only fits one parameter . However, our results repeatedly showed that the shift with time in optical properties of the collagen gels was more clearly observed in , not . Second, the main difference in the attenuation term of the two models is the factor , which describes how scattered light can still reach the focus despite multiple scattering. It is that makes the mapping from to nonlinear (and hence nontrivial) at high values, to which the collagen gel data mapped.
Most of the work on fitting optical properties from OCT imaging has focused on differentiating different types of atherosclerotic plaques.18, 19, 20 Each of these studies used a different theoretical model. Van der Meer 19 used the SS model and showed differences in between normal arteries and atherosclerotic plaques. Levitz 18 used the extended Huygens–Fresnel (EHF) model21 and showed differences in anisotropy between normal arteries and atherosclerotic plaques. The EHF model has the added advantage of accounting for multiple scattering and defocusing; that is, it is capable of evaluating both and a measure of anisotropy similar to without the need for the (time consuming) quasi-focus tracking scheme implemented in our study. However, the EHF model is highly nonlinear and thus requires a computationally burdensome algorithm that does not always converge to a solution. Alternatively, there are other models that have been used to measure optical properties from highly scattering samples. The model used by Xu 20 is very similar to the model used here, but with key differences. First, the is missing from this model, and second, is represented as the backscattering coefficient , with no nonlinear mapping back to .
Similar to our work, Xu 22 characterized cytochrome oxidase activity in astrocytes suspended in agarose gels by measuring changes in via the fitted attenuation, but using the SS model with its inherent limitations. Nonetheless, it is also worth noting that in their cell-agarose samples, a negligible amount of scattering and absorption from the agarose contributed to the signal; thus, all the scattering and absorption measured could be attributed to the cells. Figure 5 shows that cells are mild scatterers and that a very high density is needed to obtain a measureable change in attenuation. The attenuation values measured by Xu 22 required very high cell seeding densities , which is in general agreement with our results.
The experiments presented here highlight several advantages to measuring local optical properties of samples using OCT. The method is nondestructive, which makes continuous and time-lapse monitoring of the tissue gels feasible. Because this OCT method depends on endogenous contrast, no labels are needed, which has the added benefits of both preserving the sample unperturbed in its native environment and reducing the costs associated with fluorescent labels. Moreover, the visual information contained in the OCT images complements the quantitative optical property data. And although these experiments were destructive in nature—samples were prepared, removed from the sterile incubator and imaged, and then fixed for histological analysis—one of OCT’s primary advantages is that it is a fiber-based technology that can be sterilized and integrated into the incubator for continuous, time-lapse monitoring.6 The adaptation of confocal mosaicing technology23 to OCT will enable rapid surveying of large sample volumes. These features, if successfully implemented, could enable automated nondestructive evaluation of large volumes of engineered tissues.
In this study, samples were imaged and images were further divided into volume regions of about , or . This resulted in a rather large number of ROIs with volumes on the nanoliter scale, which partly affected the analyses. Implicit in the definition of and is the assumption of averaging over some volume. The nanoliter volumes assessed here approach the size limit over which these properties are defined, but nonetheless, it is still possible to define optical properties within these volumes. These small volumes, from which optical properties are measured, are a direct consequence of the spatial resolution and field of view of OCT images. Indeed, ROI volumes on the nanoliter scale have been used by others when measuring optical properties from OCT data.18, 19, 20, 22
One consequence of these small sampling volumes is an increase in the variance of the optical properties, which may simply reflect the intrinsic heterogeneity of the samples on this size scale. The variance in optical properties of acellular gels was much larger than those in SMC gels (Figs. 3 and 4, Table 2), paradoxically suggesting that acellular gels were less homogeneous than SMC gels. However, this is not the complete picture. A careful examination of the acellular gels with the naked eye showed that some regions are cloudier than others, confirming millimeter-scale heterogeneity, seen in the OCT images and observed in the optical property data. Moreover, the histological sections [Figs. 2, 2, 2, 2, 2, 2, 2, 2, 2] showed that within the field of view of an OCT or histology image, the cells were generally uniformly distributed. Noting that cells are mild scatterers but with a higher reflectivity than the background collagen, the scattering in the SMC gels from the “uniformly” distributed cells contrasted with the collagen background and shifted the optical properties accordingly. However, in the acellular gels, all the scattering was due to the slowly varying background, which explains the increased variance observed in the acellular gels. Indeed, perhaps the large spatial variability in collagen distribution may offer clues as to why some regions in tissue-engineered implants are at increased risk of mechanical failure.24
The optical properties at five days showed a consistent increase in with little change in , which correspond to a decrease in anisotropy with a modest decrease in . Such changes in optical properties suggest a decrease in the overall size of the scattering particles, not an increase in particle density, which would increase , and hence increase both and . In attempting to identify which microstructural changes caused this shift in optical properties, we (initially) proposed four potential reasons: increased collagen density, increased cell density, SMC contraction (i.e., a change in cell morphology), and collagen remodeling. Further analyses of some of the data presented here rules out increases in collagen and cell density, and preliminary results (not shown) suggest that collagen remodeling is responsible for the changes in optical properties.
To test if increased collagen density (i.e., gel compaction) caused the twofold increase in reflectivity, the collagen density was plotted against both and for the different collagen gel samples and the correlation coefficient was calculated. There was little to no correlation between collagen density and optical properties: for versus collagen density, ; for versus collagen density, (data not shown). These poor correlations suggest that collagen density was not the primary cause of the measured shifts in optical properties. Furthermore, the changes in optical properties caused by changes in collagen density were very subtle , and are well below the algorithm resolution for at . Such an increase in with little change in does not explain the shift in optical properties over time observed in our data, namely, a decrease in a modest decrease in . Also, data from set 6 (Table 2), in which higher density acellular gels were imaged at five days, suggests that gel compaction did not cause the change in optical properties, because the characteristic increase in reflectivity (with no change in attenuation) that was consistently observed in the SMC gels was absent. Of note, we do not expect that extracellular matrix laid down by SMCs over the five-day incubation period contributed significantly to the change in optical properties, because previous work has shown the matrix synthesized by SMCs embedded in both collagen and fibrin gels to make up of the total collagen content.25, 26
Additionally, the correlation between the optical properties and cell density was checked. An increase in cell density may be caused either by cell proliferation and/or compounded by gel compaction. However, an alamar blue assay on three extra sets of gels did not show an increase in cell metabolism over five days, suggesting no proliferation took place (data not shown). The result is in agreement with previously published results.27 Nonetheless, to confirm whether cell density correlates with optical properties, cells were manually counted in 15 histological H&E sections for each gel (three analogous sites in five sections separated by ), which were normalized by the area of the collagen in that section. These cell density counts were compared to both the and values for that gel, and the correlation coefficient was calculated. The correlations between the cell density and both and were rather weak. For versus cell density, , whereas for versus cell density (data not shown). These weak correlations, together with the data showing how cell density correlates with but not (Fig. 5) suggest that although an increase in cell density contributed partly to an increase in , it was not the primary reason for the twofold increase in shown in the data.
The data showed that over time, there was an overall decrease in anisotropy factor of the collagen gels, as well as a moderate decrease in scattering coefficient . The decrease in seems a bit counterintuitive because it is expected for there to be an overall increase in the density of scattering particles. However, treating the scattering particles in the gels as discretely sized microspheres, Mie theory provides some additional insight.11 Mie theory states that the relationship between the particle size and scattering cross section is highly nonlinear. Thus, consider two scattering solutions, one made with small microspheres and the other with large microspheres, but where in both cases the spheres occupy the same volume fraction. It is obvious that the small sphere solution will have a smaller value than the large sphere solution. However, these two scattering media will also have different values, with the small sphere solution having a smaller than the large sphere solution. Thus, if the increase in the number of small scattering particles (as indicated by the higher reflectivity) was caused by a breakdown of larger particles, one would expect a moderate decrease in the along with a decrease in . Such a scenario supports the hypothesis that remodeling of the collagen matrix into small fibril fragments caused the observed changes in optical properties.
In this study, OCT was used to image development in several sets of collagen gels with different cell seeding densities, as well as acellular gels, over a five-day period. After the first , changes in the collagen gels made structures (i.e., cells) impossible to resolve. To characterize the samples quantitatively, the optical properties of the gels were determined using an image-processing algorithm that fit the OCT signal from a region of interest to a theoretical model. In doing so, the observed fit results (attenuation and reflectivity ) were mapped back to the optical properties, the scattering coefficient and the scattering anisotropy factor . Within the first , the optical properties of the SMC and acellular gels were similar (Figs. 3 and 4), suggesting that the optical properties were dominated by scattering from collagen fibrils. Over time, collagen gels with SMCs consistently exhibited a twofold increase in reflectivity, while acellular gels did not (Figs. 3 and 4, Table 2). The shift in optical properties toward higher reflectivity was gradual. A sizeable subpopulation of the SMC gels at five days had reflectivity values roughly tenfold higher than those in the first were consistently observed in five sets of gels. Such a shift corresponded to a decrease in and an overall decrease in the size of the scattering particles. Additionally, at the earliest imaging time point ( , Fig. 5), cell seeding density was strongly correlated to , but not to . All in all, the data suggest that SMCs change the architecture of the collagen gel and that OCT can visualize and quantify these changes.