Summary

Quantifying Fibrillar Collagen Organization with Curvelet Transform-Based Tools

Published: November 11, 2020
doi:

Summary

Here, we present a protocol to use a curvelet transform-based, open-source MATLAB software tool for quantifying fibrillar collagen organization in the extracellular matrix of both normal and diseased tissues. This tool can be applied to images with collagen fibers or other types of line-like structures.

Abstract

Fibrillar collagens are prominent extracellular matrix (ECM) components, and their topology changes have been shown to be associated with the progression of a wide range of diseases including breast, ovarian, kidney, and pancreatic cancers. Freely available fiber quantification software tools are mainly focused on the calculation of fiber alignment or orientation, and they are subject to limitations such as the requirement of manual steps, inaccuracy in detection of the fiber edge in noisy background, or lack of localized feature characterization. The collagen fiber quantitation tool described in this protocol is characterized by using an optimal multiscale image representation enabled by curvelet transform (CT). This algorithmic approach allows for the removal of noise from fibrillar collagen images and the enhancement of fiber edges to provide location and orientation information directly from a fiber, rather than using the indirect pixel-wise or window-wise information obtained from other tools. This CT-based framework contains two separate, but linked, packages named “CT-FIRE” and “CurveAlign” that can quantify fiber organization on a global, region of interest (ROI), or individual fiber basis. This quantification framework has been developed for more than ten years and has now evolved into a comprehensive and user-driven collagen quantification platform. Using this platform, one can measure up to about thirty fiber features including individual fiber properties such as length, angle, width, and straightness, as well as bulk measurements such as density and alignment. Additionally, the user can measure fiber angle relative to manually or automatically segmented boundaries. This platform also provides several additional modules including ones for ROI analysis, automatic boundary creation, and post-processing. Using this platform does not require prior experience of programming or image processing, and it can handle large datasets including hundreds or thousands of images, enabling efficient quantification of collagen fiber organization for biological or biomedical applications.

Introduction

Fibrillar collagens are prominent, structural ECM components. Their organization changes impact tissue function and are likely associated with the progression of many diseases ranging from osteogenesis imperfecta1, cardiac dysfunction2, and wound healing3 to different types of cancer including breast4,5,6, ovarian7,8, kidney9, and pancreatic cancers10. Many established imaging modalities can be used to visualize fibrillar collagen such as second harmonic generation microscopy11, stains or dyes in conjunction with bright field or fluorescence microscopy or polarized light microscopy12, liquid crystal-based polarization microscopy (LC-PolScope)13, and electron microscopy14. As the importance of fibrillar collagen organization has become clearer, and the use of these methods has increased, the need for improved collagen fiber analysis approaches has also grown.

There have been many efforts to develop computational methods for automated measurement of fibrillar collagen. Freely available software tools are mainly focused on the calculation of fiber alignment or orientation by adopting either first derivative or structure tensor for pixels15,16, or Fourier transform-based spectrum analysis for image tiles17. All these tools are subject to limitations such as the requirement of manual steps, inaccuracy in detection of the fiber edge in noisy background, or lack of localized feature characterization.

Compared to other freely available open-source free software tools, the methods described in this protocol use CT—an optimal, multiscale, directional image representation method—to remove noise from fibrillar collagen images and enhance or track fiber edges. Information about location and orientation can be provided directly from a fiber rather than by using the indirect pixel-wise or window-wise information to infer the metrics of fiber organization. This CT-based framework18,19,20,21 can quantify fiber organization on a global, ROI, or fiber basis, mainly via two separate, but linked, packages named “CT-FIRE”18,21 and “CurveAlign”19,21. As far as the implementation of the software is concerned, in CT-FIRE, CT coefficients on multiple scales can be used to reconstruct an image that enhances edges and reduces noise. Then, an individual fiber extraction algorithm is applied to the CT-reconstructed image to track fibers for finding their representative center points, extending fiber branches from the center points, and linking fiber branches to form a fiber network. In CurveAlign, CT coefficients on a user-specified scale can be used to track local fiber orientation, where the orientation and locations of curvelets are extracted and grouped to estimate the fiber orientation at the corresponding locations. This resulting quantification framework has been developed for more than ten years and has evolved greatly in many aspects such as functionality, user interface, and modularity. For instance, this tool can visualize local fiber orientation and allows the user to measure up to thirty fiber features including individual fiber properties such as length, angle, width, and straightness, as well as bulk measurements such as density and alignment. Additionally, the user can measure fiber angle relative to manually or automatically segmented boundaries, which, for example, plays an important role in image-based biomarker development in breast cancer22 and pancreatic cancer studies10. This platform provides several feature modules including ones for ROI analysis, automatic boundary creation, and post-processing. The ROI module can be used to annotate different shapes of ROI and conduct corresponding ROI analysis. As an application example, the automatic boundary creation module can be used to register hematoxylin and eosin (H&E) bright field images with second harmonic generation (SHG) images and generate the image mask of tumor boundaries from the registered H&E images. The post-processing module can help facilitate the processing and integration of output data files from individual images for possible statistical analysis.

This quantification platform does not require prior experience of programming or image processing and can handle large datasets including hundreds or thousands of images, enabling efficient quantification of collagen organization for biological or biomedical applications. It has been widely used in different research fields by many researchers all over the world, including ourselves. There are four main publications on CT-FIRE and CurveAlign18,19,20,21, out of which the first three have been cited 272 times (as of 2020-05-04 according to Google Scholar). A review of the publications that cited this platform (CT-FIRE or CurveAlign) indicates that there are approximately 110 journal papers that directly used it for their analysis, in which approximately 35 publications were collaborative with our group, and the others (~ 75) were written by other groups. For instance, this platform was used for the following studies: breast cancer22,23,24, pancreatic cancer10,25, kidney cancer9,26, wound healing3,27,28,29,30, ovarian cancer8,31,7, uterosacral ligament32, hypophosphatemic dentin33, basal cell carcinoma34, hypoxic sarcoma35, cartilage tissue36, cardiac dysfunction37, neurons38, glioblastoma39, lymphatic contractions40, fibrous cacffolds41, gastric cancer42, microtubule43, and bladder fibrosis44. Figure 1 demonstrates the cancer imaging application of CurveAlign to find the tumor-associated collagen signatures of breast cancer19 from the SHG image. Figure 2 describes a typical schematic workflow of this platform. Although these tools have been reviewed technically18,19,21, and a regular protocol20 for alignment analysis with CurveAlign is also available, a visual protocol that demonstrates all the essential features could be useful. A visualized protocol, as presented here, will facilitate the learning process of using this platform as well as more efficiently address concerns and questions that users might have.

Protocol

NOTE: This protocol describes the use of CT-FIRE and CurveAlign for collagen quantification. These two tools have complementary, but different, main goals and are linked together to some extent. CT-FIRE can be launched from the CurveAlign interface to conduct most operations except for advanced post-processing and ROI analysis. For a full operation of CT-FIRE, it should be launched separately. 1. Image collection and image requirement NOTE: The tool can process any im…

Representative Results

These methods have been successfully applied in numerous studies. Some typical applications include: 1) Conklin et al.22 used CurveAlign to calculate tumor-associated collagen signatures, and found that collagen fibers were more frequently aligned perpendicularly to the duct perimeter in ductal carcinoma in situ (DCIS) lesions; 2) Drifka et al.10 used the CT-FIRE mode in CurveAlign to quantify the stromal collagen alignment for pancreatic ductal adenocarcinoma and normal/ch…

Discussion

This protocol describes the use of CT-FIRE and CurveAlign for fibrillar collagen quantification and can be applied to any image with collagen fibers or other line-like or fiber-like elongated structures suitable for analysis by CT-FIRE or CurveAlign. For example, elastin or elastic fibers could be processed in a similar way on this platform. We have tested both tools on computationally generated synthetic fibers21. Depending on the application, users should choose the analysis mode that is most ap…

Divulgaciones

The authors have nothing to disclose.

Acknowledgements

We thank many contributors and users to CT-FIRE and CurveAlign over the years, including Dr. Rob Nowak, Dr. Carolyn Pehlke, Dr. Jeremy Bredfeldt, Guneet Mehta, Andrew Leicht, Dr. Adib Keikhosravi, Dr. Matt Conklin, Dr. Jayne Squirrell, Dr. Paolo Provenzano, Dr. Brenda Ogle, Dr. Patricia Keely, Dr. Joseph Szulczewski, Dr. Suzanne Ponik and additional technical contributions from Swati Anand and Curtis Rueden. This work was supported by funding from Semiconductor Research Corporation, Morgridge Institute for Research, and NIH grants R01CA199996, R01CA181385 and U54CA210190 to K.W.E.

Materials

CT-FIRE Univerity of Wisconsin-Madison N/A open source software available from https://eliceirilab.org/software/ctfire/
CurveAlign University of Wisconsin-Madison N/A open source software available from https://eliceirilab.org/software/curvealign/

Referencias

  1. Nadiarnykh, O., et al. Second harmonic generation imaging microscopy studies of osteogenesis imperfecta. Journal of Biomedical Optics. 12 (5), 051805 (2007).
  2. Kouris, N. A., et al. A nondenatured, noncrosslinked collagen matrix to deliver stem cells to the heart. Regenerative Medicine. 6 (5), 569-582 (2011).
  3. LeBert, D. C., et al. Matrix metalloproteinase 9 modulates collagen matrices and wound repair. Development. 142 (12), 2136-2146 (2015).
  4. Provenzano, P. P., et al. Collagen reorganization at the tumor-stromal interface facilitates local invasion. BMC Medicine. 4 (1), 38 (2006).
  5. Provenzano, P. P., et al. Collagen density promotes mammary tumor initiation and progression. BMC Medicine. 6 (1), 1 (2008).
  6. Conklin, M. W., et al. Aligned collagen is a prognostic signature for survival in human breast carcinoma. The American Journal of Pathology. 178 (3), 1221-1232 (2011).
  7. Alkmin, S., et al. Migration dynamics of ovarian epithelial cells on micro-fabricated image-based models of normal and malignant stroma. Acta Biomaterialia. 100, 92-104 (2019).
  8. Campbell, K. R., Campagnola, P. J. Assessing local stromal alterations in human ovarian cancer subtypes via second harmonic generation microscopy and analysis. Journal of Biomedical Optics. 22 (11), 116008 (2017).
  9. Best, S. L., et al. Collagen organization of renal cell carcinoma differs between low and high grade tumors. BMC Cancer. 19 (1), 490 (2019).
  10. Drifka, C. R., et al. Periductal stromal collagen topology of pancreatic ductal adenocarcinoma differs from that of normal and chronic pancreatitis. Modern Pathology. 28 (11), 1470-1480 (2015).
  11. Campagnola, P. J., et al. Three-dimensional high-resolution second-harmonic generation imaging of endogenous structural proteins in biological tissues. Biophysical Journal. 82 (1), 493-508 (2002).
  12. Arun Gopinathan, P., et al. Study of collagen birefringence in different grades of oral squamous cell carcinoma using picrosirius red and polarized light microscopy. Scientifica. 2015, 1-7 (2015).
  13. Keikhosravi, A., et al. Quantification of collagen organization in histopathology samples using liquid crystal based polarization microscopy. Biomedical Optics Express. 8 (9), 4243-4256 (2017).
  14. Quan, B. D., Sone, E. D. Cryo-TEM analysis of collagen fibrillar structure. Methods in Enzymology. 532, 189-205 (2013).
  15. Boudaoud, A., et al. FibrilTool, an ImageJ plug-in to quantify fibrillar structures in raw microscopy images. Nature Protocols. 9 (2), 457-463 (2014).
  16. Rezakhaniha, R., et al. Experimental investigation of collagen waviness and orientation in the arterial adventitia using confocal laser scanning microscopy. Biomechanics and Modeling in Mechanobiology. 11 (3-4), 461-473 (2012).
  17. Kartasalo, K., et al. CytoSpectre: a tool for spectral analysis of oriented structures on cellular and subcellular levels. BMC Bioinformatics. 16 (1), 1 (2015).
  18. Bredfeldt, J. S., et al. Computational segmentation of collagen fibers from second-harmonic generation images of breast cancer. Journal of Biomedical Optics. 19 (1), 016007 (2014).
  19. Bredfeldt, J. S., et al. Automated quantification of aligned collagen for human breast carcinoma prognosis. Journal of Pathology Informatics. 5 (1), 28 (2014).
  20. Liu, Y., Keikhosravi, A., Mehta, G. S., Drifka, C. R., Eliceiri, K. W. Methods for quantifying fibrillar collagen alignment. Methods in Molecular Biology. 1627, 429-451 (2017).
  21. Liu, Y., et al. Fibrillar collagen quantification with curvelet transform based computational methods. Frontiers in Bioengineering and Biotechnology. 8, 198 (2020).
  22. Conklin, M. W., et al. Collagen alignment as a predictor of recurrence after ductal carcinoma in situ. Cancer Epidemiology and Prevention Biomarkers. 27 (2), 138-145 (2018).
  23. Jallow, F., et al. Dynamic interactions between the extracellular matrix and estrogen activity in progression of ER+ breast cancer. Oncogene. 38 (43), 6913-6925 (2019).
  24. Smirnova, T., et al. Serpin E2 promotes breast cancer metastasis by remodeling the tumor matrix and polarizing tumor associated macrophages. Oncotarget. 7 (50), 82289 (2016).
  25. Fanous, M., Keikhosravi, A., Kajdacsy-Balla, A., Eliceiri, K. W., Popescu, G. Quantitative phase imaging of stromal prognostic markers in pancreatic ductal adenocarcinoma. Biomedical Optics Express. 11 (3), 1354-1364 (2020).
  26. Jiménez-Torres, J. A., Virumbrales-Muñoz, M., Sung, K. E., Lee, M. H., Abel, E. J., Beebe, D. J. Patient-specific organotypic blood vessels as an in vitro model for anti-angiogenic drug response testing in renal cell carcinoma. EBioMedicine. 42, 408-419 (2019).
  27. Govindaraju, P., Todd, L., Shetye, S., Monslow, J., Puré, E. CD44-dependent inflammation, fibrogenesis, and collagenolysis regulates extracellular matrix remodeling and tensile strength during cutaneous wound healing. Matrix Biology. 75, 314-330 (2019).
  28. Henn, D., et al. Cryopreserved human skin allografts promote angiogenesis and dermal regeneration in a murine model. International Wound Journal. 17 (4), 925-936 (2020).
  29. Rico-Jimenez, J., et al. Non-invasive monitoring of pharmacodynamics during the skin wound healing process using multimodal optical microscopy. BMJ Open Diabetes Research and Care. 8 (1), 000974 (2020).
  30. Israel, J. S., et al. Quantification of collagen organization after nerve repair. Plastic and Reconstructive Surgery Global Open. 5 (12), (2017).
  31. Rentchler, E. C., Gant, K. L., Drapkin, R., Patankar, M., Campagnola, P. J. Imaging collagen alterations in STICs and high grade ovarian cancers in the fallopian tubes by second harmonic generation microscopy. Cancers. 11 (11), 1805 (2019).
  32. Hu, C., et al. Imaging collagen properties in the uterosacral ligaments of women with pelvic organ prolapse using spatial light interference microscopy (SLIM). Frontiers in Physics. 7, 72 (2019).
  33. Guirado, E., et al. Disrupted protein expression and altered proteolytic events in hypophosphatemic dentin can be rescued by dentin matrix protein 1. Frontiers in Physiology. 11, 82 (2020).
  34. Kiss, N., et al. Quantitative analysis on ex vivo nonlinear microscopy images of basal cell carcinoma samples in comparison to healthy skin. Pathology & Oncology Research. 25 (3), 1015-1021 (2019).
  35. Lewis, D. M., et al. Collagen fiber architecture regulates hypoxic sarcoma cell migration. ACS Biomaterials Science & Engineering. 4 (2), 400-409 (2018).
  36. Moura, C. C., Bourdakos, K. N., Tare, R. S., Oreffo, R. O., Mahajan, S. Live-imaging of Bioengineered Cartilage Tissue using Multimodal Non-linear Molecular Imaging. Scientific Reports. 9 (1), 1-9 (2019).
  37. Murtada, S. I., et al. Paradoxical aortic stiffening and subsequent cardiac dysfunction in Hutchinson-Gilford progeria syndrome. Journal of The Royal Society Interface. 17 (166), 0066 (2020).
  38. Nichol, R. H., Catlett, T. S., Onesto, M. M., Hollender, D., Gómez, T. M. Environmental elasticity regulates cell-type specific RHOA signaling and neuritogenesis of human neurons. Stem Cell Reports. 13 (6), 1006-1021 (2019).
  39. Pointer, K. B., et al. Association of collagen architecture with glioblastoma patient survival. Journal of Neurosurgery. 126 (6), 1812-1821 (2016).
  40. Razavi, M. S., Leonard-Duke, J., Hardie, B., Dixon, J. B., Gleason, R. L. Axial stretch regulates rat tail collecting lymphatic vessel contractions. Scientific Reports. 10 (1), 1-11 (2020).
  41. Xue, Y., et al. Valve leaflet-inspired elastomeric scaffolds with tunable and anisotropic mechanical properties. Polymers for Advanced Technologies. 31 (1), 94-106 (2020).
  42. Zhou, Z. H., et al. Reorganized collagen in the tumor microenvironment of gastric cancer and its association with prognosis. Journal of Cancer. 8 (8), 1466 (2017).
  43. Zinn, A., et al. The small GTPase RhoG regulates microtubule-mediated focal adhesion disassembly. Scientific Reports. 9 (1), 1-15 (2019).
  44. Zwaans, B. M., et al. Radiation cystitis modeling: A comparative study of bladder fibrosis radio-sensitivity in C57BL/6, C3H, and BALB/c mice. Physiological Reports. 8 (4), (2020).
  45. Devine, E. E., Liu, Y., Keikhosravi, A., Eliceiri, K. W., Jiang, J. J. Quantitative second harmonic generation imaging of leporine, canine, and porcine vocal fold collagen. The Laryngoscope. 129 (11), 2549-2556 (2019).
  46. Schindelin, J., et al. Fiji: an open-source platform for biological-image analysis. Nature Methods. 9 (7), 676-682 (2012).
  47. Zeitoune, A. A., et al. Epithelial ovarian cancer diagnosis of second-harmonic generation images: A semiautomatic collagen fibers quantification protocol. Cancer Informatics. 16, (2017).
  48. Wershof, E., et al. Matrix feedback enables diverse higher-order patterning of the extracellular matrix. PLoS Computational Biology. 15 (10), 1007251 (2019).
  49. Mostaço-Guidolin, L. B., et al. Fractal dimension and directional analysis of elastic and collagen fiber arrangement in unsectioned arterial tissues affected by atherosclerosis and aging. Journal of Applied Physiology. 126 (3), 638-646 (2019).
  50. Thain, D., Tannenbaum, T., Livny, M. Distributed computing in practice: the Condor experience. Concurrency and Computation: Practice and Experience. 17 (2-4), 323-356 (2005).
check_url/es/61931?article_type=t

Play Video

Citar este artículo
Liu, Y., Eliceiri, K. W. Quantifying Fibrillar Collagen Organization with Curvelet Transform-Based Tools. J. Vis. Exp. (165), e61931, doi:10.3791/61931 (2020).

View Video