Quantification of collagen fiber organization using three-dimensional Fourier transform-second-harmonic generation imaging

Spectroscopic analysis of the sum-frequency response of the carbon-hydrogen stretching modes in collagen type I

We studied the origin of the vibrational signatures in the sum-frequency generation (SFG) spectrum of fibrillar collagen type I in the carbon-hydrogen stretching regime. For this purpose, we developed an all-reflective, laser-scanning SFG microscope with minimum chromatic aberrations and excellent retention of the polarization state of the incident beams. We performed detailed SFG measurements of aligned collagen fibers obtained from rat tail tendon, enabling the characterization of the magnitude and polarization-orientation dependence of individual tensor elements Xijk2 of collagen's nonlinear susceptibility. Using the three-dimensional atomic positions derived from published crystallographic data of collagen type I, we simulated its Xijk2 elements for the methylene stretching vibration and compared the predicted response with the experimental results. Our analysis revealed that the carbon-hydrogen stretching range of the SFG spectrum is dominated by symmetric stretching modes of methylene bridge groups on the pyrrolidine rings of the proline and hydroxyproline residues, giving rise to a dominant peak near 2942 cm-1 and a shoulder at 2917 cm-1. Weak asymmetric stretches of the methylene bridge group of glycine are observed in the region near 2870 cm-1, whereas asymmetric CH2-stretching modes on the pyrrolidine rings are found in the 2980 to 3030 cm-1 range. These findings help predict the protein's nonlinear optical properties from its crystal structure, thus establishing a connection between the protein structure and SFG spectroscopic measurements.

Fourier transform-based method for quantifying the three-dimensional orientation distribution of fibrous units

Several materials and tissues are characterized by a microstructure composed of fibrous units embedded in a ground matrix. In this paper, a novel three-dimensional (3D) Fourier transform-based method for quantifying the distribution of fiber orientations is presented. The method allows for an accurate identification of individual fiber families, their in-plane and out-of-plane dispersion, and showed fast computation times. We validated the method using artificially generated 3D images, in terms of fiber dispersion by considering the error between the standard deviation of the reconstructed and the prescribed distributions of the artificial fibers. In addition, we considered the measured mean orientation angles of the fibers and validated the robustness using a measure of fiber density. Finally, the method is employed to reconstruct a full 3D view of the distribution of collagen fiber orientations based on in vitro second harmonic generation microscopy of collagen fibers in human and mouse skin. The dispersion parameters of the reconstructed fiber network can be used to inform mechanical models of soft fiber-reinforced materials and biological tissues that account for non-symmetrical fiber dispersion.

Evaluation of growth-induced, mechanical stress in solid tumors and spatial association with extracellular matrix content

Mechanical stresses in solid tumors play an important role in tumor progression and treatment efficacy but their quantification is under-investigated. Here, we developed an experimental and computational approach to calculate growth-induced, residual stresses and applied it to the breast (4T1), pancreatic (PAN02), and fibrosarcoma (MCA205) tumor models. Following resection, tumors are embedded in agarose gels and cuts are made in two perpendicular directions to release residual stress. With the use of image processing, the detailed bulging displacement profile is measured and finite elements models of the bulging geometry are developed for the quantification of the stress levels. The mechanical properties of the tumors are measured in vivo prior to resection with shear wave elastography. We find that the average magnitude of residual stresses ranges from 3.31 to 10.88 kPa, and they are non-uniformly distributed within the tissue due to the heterogeneity of the tumor microenvironment. Interestingly, we demonstrate that a second cut can still release a significant amount of stresses. We further find a strong association of spatial hyaluronan and collagen content with the spatial profile of stress for the MCA205 and PAN02 tumors and a partial association for the 4T1. Interestingly the colocalization of hyaluronan and collagen content had a stronger association with the spatial profile of stress for MCA205, PAN02, and 4T1. Finally, measurements of the elastic modulus with shear wave elastography show a nonlinear correlation with tumor volume for the more fibrotic MCA205 and 4T1 tumors. Overall, our results provide insights for a better understanding of the mechanical behavior of tumors.

Shear wave elastography can stratify rectal cancer response to short-course radiation therapy

Rectal cancer is a deadly disease typically treated using neoadjuvant chemoradiotherapy followed by total mesorectal excision surgery. To reduce the occurrence of mesorectal excision surgery for patients whose tumors regress from the neoadjuvant therapy alone, conventional imaging, such as computed tomography (CT) or magnetic resonance imaging (MRI), is used to assess tumor response to neoadjuvant therapy before surgery. In this work, we hypothesize that shear wave elastography offers valuable insights into tumor response to short-course radiation therapy (SCRT)-information that could help distinguish radiation-responsive from radiation-non-responsive tumors and shed light on changes in the tumor microenvironment that may affect radiation response. To test this hypothesis, we performed elastographic imaging on murine rectal tumors (n = 32) on days 6, 10, 12, 16, 18, 20, 23, and 25 post-tumor cell injection. The study revealed that radiation-responsive and non-radiation-responsive tumors had different mechanical properties. Specifically, radiation-non-responsive tumors showed significantly higher shear wave speed SWS (p < 0.01) than radiation-responsive tumors 11 days after SCRT. Furthermore, there was a significant difference in shear wave attenuation (SWA) (p < 0.01) in radiation-non-responsive tumors 16 days after SCRT compared to SWA measured just one day after SCRT. These results demonstrate the potential of shear wave elastography to provide valuable insights into tumor response to SCRT and aid in exploring the underlying biology that drives tumors' responses to radiation.

Mapping the 3D remodeling of the extracellular matrix in human hypertrophic scar by multi-parametric multiphoton imaging using endogenous contrast

The hypertrophic scar is an aberrant form of wound healing process, whose clinical efficacy is limited by a lack of understanding of its pathophysiology. Remodeling of collagen and elastin fibers in the extracellular matrix (ECM) is closely associated with scar progression. Herein, we perform label-free multiphoton microscopy (MPM) of both fiber components from human skin specimens and propose a multi-fiber metrics (MFM) analysis model for mapping the structural remodeling of the ECM in hypertrophic scars in a highly-sensitive, three-dimensional (3D) manner. We find that both fiber components become wavier and more disorganized in scar tissues, while content accumulation is observed from elastin fibers only. The 3D MFM analysis can effectively distinguish normal and scar tissues with better than 95% in accuracy and 0.999 in the area under the curve value of the receiver operating characteristic curve. Further, unique organizational features with orderly alignment of both fibers are observed in scar-normal adjacent regions, and an optimized combination of features from 3D MFM analysis enables successful identification of all the boundaries. This imaging and analysis system uncovers the 3D architecture of the ECM in hypertrophic scars and exhibits great translational potential for evaluating scars in vivo and identifying individualized treatment targets.

Label-free analysis of inflammatory tissue remodeling in murine lung tissue based on multiphoton microscopy, Raman spectroscopy and machine learning

Inflammatory fibrotic tissue remodeling can lead to severe morbidity. Histopathology grading requires extraction of biopsies and elaborate tissue processing. Label-free optical technologies can provide diagnostic readout without preparation and artificial stainings and show potential for in vivo applications. Here, we present an integration of Raman spectroscopy (RS) and multiphoton microscopy for joint investigation of the bio-chemical composition and morphological features related to cellular components and connective tissue. Both modalities show that collagen signatures were significantly increased in a murine fibrosis model. Furthermore, autofluorescence signatures assigned to immune cells show high correlation with disease severity. RS indicates increased levels of elastin and lipids. Further, we investigated the effect of joint data sets on prediction performance in machine learning models. Although binary classification did not benefit from adding more features, multi-class classification was improved by integrated data sets.

Identification of human ovarian cancer relying on collagen fiber coverage features by quantitative second harmonic generation imaging

Ovarian cancer has the highest mortality rate among all gynecological cancers, containing complicated heterogeneous histotypes, each with different treatment plans and prognoses. The lack of screening test makes new perspectives for the biomarker of ovarian cancer of great significance. As the main component of extracellular matrix, collagen fibers undergo dynamic remodeling caused by neoplastic activity. Second harmonic generation (SHG) enables label-free, non-destructive imaging of collagen fibers with submicron resolution and deep sectioning. In this study, we developed a new metric named local coverage to quantify morphologically localized distribution of collagen fibers and combined it with overall density to characterize 3D SHG images of collagen fibers from normal, benign and malignant human ovarian biopsies. An overall diagnosis accuracy of 96.3% in distinguishing these tissue types made local and overall density signatures a sensitive biomarker of tumor progression. Quantitative, multi-parametric SHG imaging might serve as a potential screening test tool for ovarian cancer.

A Nonparametric Approach for Estimating Three-Dimensional Fiber Orientation Distribution Functions (ODFs) in Fibrous Materials

Many biological tissues contain an underlying fibrous microstructure that is optimized to suit a physiological function. The fiber architecture dictates physical characteristics such as stiffness, diffusivity, and electrical conduction. Abnormal deviations of fiber architecture are often associated with disease. Thus, it is useful to characterize fiber network organization from image data in order to better understand pathological mechanisms. We devised a method to quantify distributions of fiber orientations based on the Fourier transform and the Qball algorithm from diffusion MRI. The Fourier transform was used to decompose images into directional components, while the Qball algorithm efficiently converted the directional data from the frequency domain to the orientation domain. The representation in the orientation domain does not require any particular functional representation, and thus the method is nonparametric. The algorithm was verified to demonstrate its reliability and used on datasets from microscopy to show its applicability. This method increases the ability to extract information of microstructural fiber organization from experimental data that will enhance our understanding of structure-function relationships and enable accurate representation of material anisotropy in biological tissues.

Tissue Imaging and Quantification Relying on Endogenous Contrast

Cell-matrix interactions play an important role in regulating a variety of essential processes in multicellular organisms, and are closely associated with numerous diseases. Modified interactions have major effects upon key features of both cells and extracellular matrix (ECM), and a thorough understanding of changes in these features can lead to critically important insights of diseases as well as the identification of effective therapeutic targets. Here, we summarize recent advances in quantitative, optical imaging of cellular metabolism and ECM spatial organization using endogenous sources of contrast. Specifically, we focus on the two-photon excited fluorescence (TPEF) imaging of autofluorescent cellular coenzymes, NAD(P)H and FAD, for the extraction of metabolic information described by optical biomarkers including cellular redox state, NAD(P)H fluorescence lifetime, and mitochondrial clustering. We show representative applications in assessing adipose tissue function and detecting malignant lesions in human skin, and further demonstrate that a combination of these optical metrics can provide complementary insights into the underlying biological mechanisms. In addition, we review the development of quantitative analysis methods to extract spatial orientation and organization metrics of collagen fibers, a major ECM component, and demonstrate applications of these approaches in two and three dimensions in several diseases, including would healing, osteoarthritis and cancer, as well as assessments of matrix remodeling in hormone-regulated engineered breast tissues. Finally, we summarize this chapter and discuss important research directions that we expect will evolve in the near future.

Spatial frequency metrics for analysis of microscopic images of musculoskeletal tissues

Purpose: Imaging-based metrics for analysis of biological tissues are powerful tools that can extract information such as shape, size, periodicity, and many other features to assess the requested qualities of a tissue. Muscular and osseous tissues consist of periodic structures that are directly related to their function, and so analysis of these patterns likely reflects tissue health and regeneration.Methods: A method for assessment of periodic structures is by analyzing them in the spatial frequency domain using the Fourier transform. In this paper, we present two filters which we developed in the spatial frequency domain for the purpose of analyzing musculoskeletal structures. These filters provide information about 1) the angular orientation of the tissues and 2) their periodicity. We explore periodic structural patterns in the mitochondrial network of skeletal muscles that are reflective of muscle metabolism and myogenesis; and patterns of collagen fibers in the bone that are reflective of the organization and health of bone extracellular matrix.Results: We present an analysis of mouse skeletal muscle in healthy and injured muscles. We used a transgenic mouse that ubiquitously expresses fluorescent protein in their mitochondria and performed 2-photon microscopy to image the structures. To acquire the collagen structure of the bone we used non-linear SHG microscopy of mouse flat bone. We analyze and compare juvenile versus adult mice, which have different structural patterns.Conclusions: Our results indicate that these metrics can quantify musculoskeletal tissues during development and regeneration.

Quantitative Classification of 3D Collagen Fiber Organization From Volumetric Images

Collagen fibers in biological tissues have a complex 3D organization containing rich information linked to tissue mechanical properties and are affected by mutations that lead to diseases. Quantitative assessment of this 3D collagen fiber organization could help to develop reliable biomechanical models and understand tissue structure-function relationships, which impact diagnosis and treatment of diseases or injuries. While there are advanced techniques for imaging collagen fibers, published methods for quantifying 3D collagen fiber organization have been sparse and give limited structural information which cannot distinguish a wide range of 3D organizations. In this article, we demonstrate an algorithm for quantitative classification of 3D collagen fiber organization. The algorithm first simulates five groups, or classifications, of fiber organization: unidirectional, crimped, disordered, two-fiber family, and helical. These five groups are widespread in natural tissues and are known to affect the tissue's mechanical properties. We use quantitative metrics based on features such as preferred 3D fiber orientation and spherical variance to differentiate each classification in a repeatable manner. We validate our algorithm by applying it to second-harmonic generation images of collagen fibers in tendon and cervix tissue that has been sectioned in specified orientations, and we find strong agreement between classification from simulated data and the physical fiber organization. Our approach provides insight for interpreting 3D fiber organization directly from volumetric images. This algorithm could be applied to other fiber-like structures that are not necessarily made of collagen.

Three-dimensional computation of fibre orientation, diameter and branching in segmented image stacks of fibrous networks

Fibre topography of the extracellular matrix governs local mechanical properties and cellular behaviour including migration and gene expression. While quantifying properties of the fibrous network provides valuable data that could be used across a breadth of biomedical disciplines, most available techniques are limited to two dimensions and, therefore, do not fully capture the architecture of three-dimensional (3D) tissue. The currently available 3D techniques have limited accuracy and applicability and many are restricted to a specific imaging modality. To address this need, we developed a novel fibre analysis algorithm capable of determining fibre orientation, fibre diameter and fibre branching on a voxel-wise basis in image stacks with distinct fibre populations. The accuracy of the technique is demonstrated on computer-generated phantom image stacks spanning a range of features and complexities, as well as on two-photon microscopy image stacks of elastic fibres in bovine tendon and dermis. Additionally, we propose a measure of axial spherical variance which can be used to define the degree of fibre alignment in a distribution of 3D orientations. This method provides a useful tool to quantify orientation distributions and variance on image stacks with distinguishable fibres or fibre-like structures.

GPU-accelerated ray-casting for 3D fiber orientation analysis

Orientation analysis of fibers is widely applied in the fields of medical, material and life sciences. The orientation information allows predicting properties and behavior of materials to validate and guide a fabrication process of materials with controlled fiber orientation. Meanwhile, development of detector systems for high-resolution non-invasive 3D imaging techniques led to a significant increase in the amount of generated data per a sample up to dozens of gigabytes. Though plenty of 3D orientation estimation algorithms were developed in recent years, neither of them can process large datasets in a reasonable amount of time. This fact complicates the further analysis and makes impossible fast feedback to adjust fabrication parameters. In this work, we present a new method for quantifying the 3D orientation of fibers. The GPU implementation of the proposed method surpasses another popular method for 3D orientation analysis regarding accuracy and speed. The validation of both methods was performed on a synthetic dataset with varying parameters of fibers. Moreover, the proposed method was applied to perform orientation analysis of scaffolds with different fibrous micro-architecture studied with the synchrotron μCT imaging setup. Each acquired dataset of size 600x600x450 voxels was analyzed in less 2 minutes using standard PC equipped with a single GPU.

Noise reduction and quantification of fiber orientations in greyscale images

Quantification of the angular orientation distribution of fibrous tissue structures in scientific images benefits from the Fourier image analysis to obtain quantitative information. Measurement uncertainties represent a major challenge and need to be considered by propagating them in order to determine an adaptive anisotropic Fourier filter. Our adaptive filter method (AF) is based on the maximum relative uncertainty δcut of the power spectrum as well as a weighted radial sum with weighting factor α. We use a Monte-Carlo simulation to obtain realistic greyscale images that include defined variations in fiber thickness, length, and angular dispersion as well as variations in noise. From this simulation the best agreement between predefined and derived angular orientation distribution is found for evaluation parameters δcut = 2.1% and α = 1.5. The resulting cumulative orientation distribution was modeled by a sigmoid function to obtain the mean angle and the fiber dispersion. A comparison to a state-of-the-art band-pass method revealed that the AF method is more suitable for the application on greyscale fiber images, since the error of the fiber dispersion significantly decreased from (33.9 ± 26.5)% to (13.2 ± 12.7)%. Both methods were found to accurately quantify the mean fiber orientation with an error of (1.9 ± 1.5)° and (2.3 ± 2.1)° in case of the AF and the band-pass method, respectively. We demonstrate that the AF method is able to accurately quantify the fiber orientation distribution in in vivo second-harmonic generation images of dermal collagen with a mean fiber orientation error of (6.0 ± 4.0)° and a dispersion error of (9.3 ± 12.1)%.

Demonstration of flat-top beam illumination in widefield multiphoton microscopy

Multiphoton microscopy provides a suitable technique for imaging biological tissues with submicrometer resolution. Usually a Gaussian beam (GB) is used for illumination, leading to a reduced power efficiency in the multiphoton response and vignetting for a square-shaped imaging area. A flat-top beam (FTB) provides a uniform spatial intensity distribution that equalizes the probability of a multiphoton effect across the imaging area. We employ a customized widefield multiphoton microscope to compare the performance of a square-shaped FTB illumination with that based on using a GB, for both two-photon fluorescence (TPF) and second-harmonic generation (SHG) imaging. The variation in signal-to-noise ratio across TPF images of fluorescent dyes spans ∼5.6  dB for the GB and ∼1.2  dB for the FTB illumination, respectively. For the GB modality, TPF images of mouse colon and Convallaria root, and SHG images of chicken tendon and human breast biopsy tissue showcase ∼20  %   area that are not imaged due to either insufficient or lack of illumination. For quantitative analysis that depends on the illuminated area, this effect can potentially lead to inaccuracies. This work emphasizes the applicability of FTB illumination to multiphoton applications.

Second harmonic generation imaging reveals asymmetry in the rotational helicity of collagen lamellae in chicken corneas

High tensile strength and optical clarity are unique properties of the cornea. These features are dictated by the three-dimensional architecture of corneal lamellae. Therefore, understanding the microscopic details of the cornea's structural organization may contribute to the development of artificial cornea for the treatment of corneal diseases. In this study, the combination of forward second harmonic generation (SHG) microcopy and fast Fourier-transform based image analysis was used to characterize the depth-dependent superstructure of chicken corneal stroma. Our results show that from the surface, adjacent lamellae of anterior chicken cornea lamella rotate in a counterclockwise direction, and the same rotational helicity is observed in left and right corneas. Furthermore, the overall average rotational pitch of lamellae is 0.92 ± 0.11 degree/µm which persists for 176 ± 14 µm in the anterior stroma. As depth further increased, the rate of lamellar rotation decreases. Upon reaching posterior stroma, lamellar orientation remains constant. Throughout the stroma, collagen lamellae in chicken rotate a total of 169 ± 21 degrees. The lack of lamellar rotation in posterior stroma suggests that packing efficiency cannot be used to explain the helicity of depth-dependent rotation of anterior stroma. In addition, although the right cornea has a higher rotational pitch (0.95 ± 11 vs 0.90 ± 10 degrees/µm) and thinner anterior stroma (173 ± 13 vs 179 ± 14 µm) than the left cornea, the two effects cancel each other out and result in similar total angular rotation of anterior stroma (161 ± 23 and 165 degrees ± 21). Finally, our observation of a total angular rotation of 169 ± 21 degrees shows that within experimental error, chicken cornea lamellae rotate around 180 degrees or half of a complete turn. Additional studies are needed to arrive at an explanation of chicken superstructure in three dimensions.

High-speed collagen fiber modeling and orientation quantification for optical coherence tomography imaging

Quantifying collagen fiber architecture has clinical and scientific relevance across a variety of tissue types and adds functionality to otherwise largely qualitative imaging modalities. Optical coherence tomography (OCT) is uniquely suited for this task due to its ability to capture the collagen microstructure over larger fields of view than traditional microscopy. Existing image processing techniques for quantifying fiber architecture, while accurate and effective, are very slow for processing large datasets and tend to lack structural specificity. We describe here a computationally efficient method for quantifying and visualizing collagen fiber organization. The algorithm is demonstrated on swine atria, bovine anterior cruciate ligament, and human cervical tissue samples. Additionally, we show an improved performance for images with crimped fiber textures and low signal to noise when compared to similar methods.

3D organizational mapping of collagen fibers elucidates matrix remodeling in a hormone-sensitive 3D breast tissue model

Hormones play an important role in normal and diseased breast tissue development. However, they can also disrupt cell-matrix interactions and their role in extracellular matrix reorganization during epithelial morphogenesis remains poorly understood, partly due to a lack of sensitive approaches for matrix characterization. Here, we assess the hormonal regulation of matrix reorganization in a three-dimensional (3D) breast tissue culture model using a novel metric, i.e., 3D directional variance, to characterize the 3D organization of collagen fibers visualized via high-resolution, second harmonic generation imaging. This metric enables resolving and quantifying patterns of spatial organization throughout the matrix surrounding epithelial structures treated with 17β-estradiol (E2) alone, and E2 in combination with either promegestone, a progestogen, or prolactin. Addition of promegestone results in the most disorganized fibers, while the E2 alone treatment leads to the most organized ones. Location-dependent organization mapping indicates that only the prolactin treatment leads to significant heterogeneities in the regional organization of collagen fibers, with higher levels of alignment observed at the end of the elongated epithelial structures. The observed collagen organization patterns for all groups persist for tens of micrometers. In addition, a comparison between 3D directional variance and typical 2D analysis approaches reveals an improved sensitivity of the 3D metric to identify organizational heterogeneities and differences among treatment groups. These results demonstrate that 3D directional variance is sensitive to subtle changes in the extracellular micro-environment and has the potential to elucidate reciprocal cell-matrix interactions in the context of numerous applications involving the study of normal and diseased tissue morphogenesis.

Quantization of collagen organization in the stroma with a new order coefficient

Many optical and biomechanical properties of the cornea, specifically the transparency of the stroma and its stiffness, can be traced to the degree of order and direction of the constituent collagen fibers. To measure the degree of order inside the cornea, a new metric, the order coefficient, was introduced to quantify the organization of the collagen fibers from images of the stroma produced with a custom-developed second harmonic generation microscope. The order coefficient method gave a quantitative assessment of the differences in stromal collagen arrangement across the cornea depths and between untreated stroma and cross-linked stroma.

Techniques to assess bone ultrastructure organization: orientation and arrangement of mineralized collagen fibrils

Bone's remarkable mechanical properties are a result of its hierarchical structure. The mineralized collagen fibrils, made up of collagen fibrils and crystal platelets, are bone's building blocks at an ultrastructural level. The organization of bone's ultrastructure with respect to the orientation and arrangement of mineralized collagen fibrils has been the matter of numerous studies based on a variety of imaging techniques in the past decades. These techniques either exploit physical principles, such as polarization, diffraction or scattering to examine bone ultrastructure orientation and arrangement, or directly image the fibrils at the sub-micrometre scale. They make use of diverse probes such as visible light, X-rays and electrons at different scales, from centimetres down to nanometres. They allow imaging of bone sections or surfaces in two dimensions or investigating bone tissue truly in three dimensions, in vivo or ex vivo, and sometimes in combination with in situ mechanical experiments. The purpose of this review is to summarize and discuss this broad range of imaging techniques and the different modalities of their use, in order to discuss their advantages and limitations for the assessment of bone ultrastructure organization with respect to the orientation and arrangement of mineralized collagen fibrils.

Fitting-free algorithm for efficient quantification of collagen fiber alignment in SHG imaging applications

Collagen fiber alignment derived from second harmonic generation (SHG) microscopy images can be important for disease diagnostics. Image processing algorithms are needed to robustly quantify the alignment in images with high sensitivity and reliability. Fourier transform (FT) magnitude, 2D power spectrum, and image autocorrelation have previously been used to extract fiber information from images by assuming a certain mathematical model (e.g. Gaussian distribution of the fiber-related parameters) and fitting. The fitting process is slow and fails to converge when the data is not Gaussian. Herein we present an efficient constant-time deterministic algorithm which characterizes the symmetricity of the FT magnitude image in terms of a single parameter, named the fiber alignment anisotropy R ranging from 0 (randomized fibers) to 1 (perfect alignment). This represents an important improvement of the technology and may bring us one step closer to utilizing the technology for various applications in real time. In addition, we present a digital image phantom-based framework for characterizing and validating the algorithm, as well as assessing the robustness of the algorithm against different perturbations.

Second-harmonic patterned polarization-analyzed reflection confocal microscope

We introduce the second-harmonic patterned polarization-analyzed reflection confocal (SPPARC) microscope-a multimodal imaging platform that integrates Mueller matrix polarimetry with reflection confocal and second-harmonic generation (SHG) microscopy. SPPARC microscopy provides label-free three-dimensional (3-D), SHG-patterned confocal images that lend themselves to spatially dependent, linear polarimetric analysis for extraction of rich polarization information based on the Mueller calculus. To demonstrate its capabilities, we use SPPARC microscopy to analyze both porcine tendon and ligament samples and find differences in both circular degree-of-polarization and depolarization parameters. Moreover, using the collagen-generated SHG signal as an endogenous counterstain, we show that the technique can be used to provide 3-D polarimetric information of the surrounding extrafibrillar matrix plus cells or EFMC region. The unique characteristics of SPPARC microscopy holds strong potential for it to more accurately and quantitatively describe microstructural changes in collagen-rich samples in three spatial dimensions.

Application of quantitative second-harmonic generation microscopy to posterior cruciate ligament for crimp analysis studies

We use second-harmonic generation (SHG) microscopy to quantitatively characterize collagen fiber crimping in the posterior cruciate ligament (PCL). The obtained SHG images are utilized to define three distinct categories of crimp organization in the PCL. Using our previously published spatial-frequency analysis, we develop a simple algorithm to quantitatively distinguish the various crimp patterns. In addition, SHG microscopy reveals both the three-dimensional structural variation in some PCL crimp patterns as well as an underlying helicity in these patterns that have mainly been observed using electron microscopy. Our work highlights how SHG microscopy could potentially be used to link the fibrous structural information in the PCL to its mechanical properties.

Automated quantification of three-dimensional organization of fiber-like structures in biological tissues

Fiber-like structures are prevalent in biological tissues, yet quantitative approaches to assess their three-dimensional (3D) organization are lacking. We develop 3D directional variance, as a quantitative biomarker of truly 3D fibrillar organization by extending the directional statistics formalism developed for describing circular data distributions (i.e. when 0° and 360° are equivalent) to axial ones (i.e. when 0° and 180° are equivalent). Significant advantages of this analysis include its time efficiency, sensitivity and ability to provide quantitative readouts of organization over different size scales of a given data set. We establish a broad range of applications for this method by characterizing collagen fibers, neuronal axons and fibroblasts in the context of cancer diagnostics, traumatic brain injury and cell-matrix interactions in developing engineered tissues. This method opens possibilities for unraveling in a sensitive, and quantitative manner the organization of essential fiber-like structures in tissues and ultimately its impact on tissue function.

From nano to macro: studying the hierarchical structure of the corneal extracellular matrix

In this review, we discuss current methods for studying ocular extracellular matrix (ECM) assembly from the 'nano' to the 'macro' levels of hierarchical organization. Since collagen is the major structural protein in the eye, providing mechanical strength and controlling ocular shape, the methods presented focus on understanding the molecular assembly of collagen at the nanometre level using X-ray scattering through to the millimetre to centimetre level using non-linear optical (NLO) imaging of second harmonic generated (SHG) signals. Three-dimensional analysis of ECM structure is also discussed, including electron tomography, serial block face scanning electron microscopy (SBF-SEM) and digital image reconstruction. Techniques to detect non-collagenous structural components of the ECM are also presented, and these include immunoelectron microscopy and staining with cationic dyes. Together, these various approaches are providing new insights into the structural blueprint of the ocular ECM, and in particular that of the cornea, which impacts upon our current understanding of the control of corneal shape, pathogenic mechanisms underlying ectatic disorders of the cornea and the potential for corneal tissue engineering.

Analysis of collagen organization in mouse achilles tendon using high-frequency ultrasound imaging

Achilles tendon ruptures are traumatic injuries, and techniques for assessing repair outcomes rely on patient-based measures of pain and function, which do not directly assess tendon healing. Consequently, there is a need for a quantitative, in vivo measure of tendon properties. Therefore, the purpose of this study was to validate ultrasound imaging for evaluating collagen organization in tendons. In this study, we compared our novel, high-frequency ultrasound (HFUS) imaging and analysis method to a standard measure of collagen organization, crossed polarizer (CP) imaging. Eighteen mouse Achilles tendons were harvested and placed into a testing fixture where HFUS and CP imaging could be performed simultaneously in a controlled loading environment. Two experiments were conducted: (1) effect of loading on collagen alignment and (2) effect of an excisional injury on collagen alignment. As expected, it was found that both the HFUS and CP methods could reliably detect an increase in alignment with increasing load, as well as a decrease in alignment with injury. This HFUS method demonstrates that structural measures of collagen organization in tendon can be determined through ultrasound imaging. This experiment also provides a mechanistic evaluation of tissue structure that could potentially be used to develop a targeted approach to aid in rehabilitation or monitor return to activity after tendon injury.

Determination of collagen fiber orientation in histological slides using Mueller microscopy and validation by second harmonic generation imaging

We studied the azimuthal orientations of collagen fibers in histological slides of uterine cervical tissue by two different microscopy techniques, namely Mueller polarimetry (MP) and Second Harmonic Generation (SHG). SHG provides direct visualization of the fibers with high specificity, which orientations is then obtained by suitable image processing. MP provides images of retardation (among other polarimetric parameters) due to the optical anisotropy of the fibers, which is enhanced by Picrosirius Red staining. The fiber orientations are then assumed to be those of the retardation slow axes. The two methods, though fully different from each other, provide quite similar maps of average fiber orientations. Overall, our results confirm that MP microscopy provides reliable images of dominant fiber orientations at a much lower cost that SHG, which remains the "gold standard" for specific imaging of collagen fibers using optical microscopy.

Graphics processing unit-based quantitative second-harmonic generation imaging

We adapt a graphics processing unit (GPU) to dynamic quantitative second-harmonic generation imaging. We demonstrate the temporal advantage of the GPU-based approach by computing the number of frames analyzed per second from SHG image videos showing varying fiber orientations. In comparison to our previously reported CPU-based approach, our GPU-based image analysis results in ∼10× improvement in computational time. This work can be adapted to other quantitative, nonlinear imaging techniques and provides a significant step toward obtaining quantitative information from fast in vivo biological processes.

Contribution of inertial cavitation in the enhancement of in vitro transscleral drug delivery

In ocular drug delivery, the sclera is a promising pathway for administering drugs to both the anterior and posterior segments of the eye. Due to the low permeability of the sclera, however, efficient drug delivery is challenging. In this study, pulsed ultrasound (US) was investigated as a potential method for enhancing drug delivery to the eye through the sclera. The permeability of rabbit scleral tissue to a model drug compound, sodium fluorescein, was measured after US-irradiation at 1.1 MHz using time-averaged acoustic powers of 0.5-5.4 W (6.8-12.8 MPa peak negative pressure), with a fixed duty cycle of 2.5% for two different pulse repetition frequencies of 100 and 1000 Hz. Acoustic cavitation activity was measured during exposures using a passive cavitation detector and was used to quantify the level of bubble activity. A correlation between the amount of cavitation activity and the enhancement of scleral permeability was demonstrated with a significant enhancement in permeability of US exposed samples compared to controls. Transmission electron microscopy showed no evidence of significant alteration in viability of tissue exposed to US exposures. A pulsed US protocol designed to maximum cavitation activity may therefore be a viable method for enhancing drug delivery to the eye.

Application of quantitative second-harmonic generation microscopy to dynamic conditions

We present a quantitative second-harmonic generation (SHG) imaging technique that quantifies the 2D spatial organization of collagen fiber samples under dynamic conditions, as an image is acquired. The technique is demonstrated for both a well-aligned tendon sample and a randomly aligned, sparsely distributed collagen scaffold sample. For a fixed signal-to-noise ratio, we confirm the applicability of this method for various window sizes (pixel areas) as well as with using a gridded overlay map that allows for correlations of fiber orientations within a given image. This work has direct impact to in vivo biological studies by incorporating simultaneous SHG image acquisition and analysis.

Application of Fourier transform-second-harmonic generation imaging to the rat cervix

We present the application of Fourier transform-second-harmonic generation (FT-SHG) imaging to evaluate the arrangement of collagen fibers in five nonpregnant rat cervices. Tissue slices from the mid-cervix and near the external orifice of the cervix were analyzed in both two-dimensions (2D) and three-dimensions (3D). We validate that the cervical microstructure can be quantitatively assessed in three dimensions using FT-SHG imaging and observe collagen fibers oriented both in and out-of-plane in the outermost and the innermost layers, which cannot be observed using 2D FT-SHG analysis alone. This approach has the potential to be a clinically applicable method for measuring progressive changes in collagen organization during cervical remodeling in humans.

Chiral imaging of collagen by second-harmonic generation circular dichroism

We provide evidence that the chirality of collagen can give rise to strong second-harmonic generation circular dichroism (SHG-CD) responses in nonlinear microscopy. Although chirality is an intrinsic structural property of collagen, most of the previous studies ignore that property. We demonstrate chiral imaging of individual collagen fibers by using a laser scanning microscope and type-I collagen from pig ligaments. 100% contrast level of SHG-CD is achieved with sub-micrometer spatial resolution. As a new contrast mechanism for imaging chiral structures in bio-tissues, this technique provides information about collagen morphology and three-dimensional orientation of collagen molecules.