research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775

The use of ethanol as contrast enhancer in synchrotron X-ray phase-contrast imaging leads to heterogeneous myocardial tissue shrinkage: a case report

crossmark logo

aBCN MedTech, Universitat Pompeu Fabra, Barcelona, Spain, bInstitute of neurosciences, Department of Medicine, School of Medicine and Health Sciences, University of Barcelona, Barcelona, Spain, cEuropean Synchrotron Radiation Facility, Grenoble, France, dInstitute of Cardiovascular Science, University College London, London, United Kingdom, eInstitut d'Investigacions Biomèdiques August Pi i Sunyer (IDIBAPS), Barcelona, Spain, fArrhythmia Unit, Department of Cardiology, Hospital Clínic de Barcelona, Barcelona, Spain, gCentre for Biomedical Research on Rare Diseases (CIBER-ER), Hospital Clínic de Barcelona, Barcelona, Spain, hPaul Scherrer Institut, Villigen PSI, Switzerland, and iCardiovascular Resarch Group iCare4Kids, Institut de Recerca Sant Joan de Déu, Esplugues de Llobregat, Spain
*Correspondence e-mail: gabriel.bernardino@upf.edu, p.garciacanadilla@gmail.com

Edited by M. Guizar-Sicairos, Paul Scherrer Institute, Switzerland (Received 27 May 2024; accepted 18 October 2024)

One of the main limitations of conventional absorption-based X-ray micro-computed tomography imaging of biological samples is the low inherent X-ray contrast of soft tissue. To overcome this limitation, the use of ethanol as contrast agent has been proposed to enhance image contrast of soft tissues through dehydration. Some authors have shown that ethanol shrinks and hardens the tissue too much, also causing small tissue ruptures due to fast dehydration. However, the local tissue deformation occurring as a consequence of tissue dehydration and whether tissue shrinkage can modify myocardial architecture has not been quantified yet. The aim of this paper is to quantify the local myocardial tissue deformation due to ethanol dehydration based on 3D non-rigid registration and perform a detailed characterization of its myocardial tissue organization, before and after ethanol dehydration. A rat adult heart was imaged with synchrotron-radiation-based X-ray phase contrast imaging (X-PCI) three times: before, 9 h after and 342 h after ethanol immersion. The total volume shrinkage as well as changes in the left ventricular myocardial thickness were computed. Then, to determine local deformation of the heart caused by ethanol dehydration, the related 3D tomographic datasets were registered by means of a non-rigid registration algorithm. Finally, changes on the orientation and organization of myocytes were assessed. Our results show that the use of ethanol in synchrotron X-PCI can improve image contrast, but the tissue shrinkage is not homogeneous thus changing the local myocardial organization.

1. Introduction

The study of the macro- and micro-structure of cardiac tissue is essential for a better understanding of the cardiac structure and function relationship, both in healthy and under pathological conditions. To date, conventional histology remains the gold standard for the analysis of cardiac tissue microstructure. Although being very powerful in terms of contrast and resolution, it involves several processes such as tissue dehydration, slicing and staining thus altering tissue native structure. In conventional histology, tissue sections are individually visualized under a microscope resulting in two-dimensional (2D) microscopic images. The possibility of reconstructing an entire volume from serial 2D slices has been extensively investigated, but still remains very challenging (Pichat et al., 2018[Pichat, J., Iglesias, J. E., Yousry, T., Ourselin, S. & Modat, M. (2018). Med. Image Anal. 46, 73-105.]).

Three-dimensional (3D) high-resolution techniques such as X-ray microcomputed tomography (micro-CT) have emerged as a powerful tool for 3D imaging of biological samples with high-resolution (<1 µm voxel size) and allowing the production of virtual tissue slices without destroying the sample, and is therefore sometimes referred to as virtual histology (Bournonville et al., 2019[Bournonville, S. de, Vangrunderbeeck, S. & Kerckhofs, G. (2019). Contrast Media Mol. Imaging, 2019, 1-9.]; Busse et al., 2018[Busse, M., Müller, M., Kimm, M. A., Ferstl, S., Allner, S., Achterhold, K., Herzen, J. & Pfeiffer, F. (2018). Proc. Natl Acad. Sci. USA, 115, 2293-2298.]; Albers et al., 2018[Albers, J., Markus, M. A., Alves, F. & Dullin, C. (2018). Sci. Rep. 8, 7712.]). Moreover, with the development of synchrotron facilities, more brilliant X-ray sources are available, thus improving image quality and resolution and significantly reducing the scanning time. However, one of the main limitations of absorption-based X-ray micro-CT imaging of biological samples is the low absorption of soft tissue. To overcome this limitation several different approaches have been proposed, from the use of different contrast agents such as iodine or phospho­tungstic acid (PTA) that bind to the tissue of interest, hence increasing its X-ray attenuation coefficient (Bournonville et al., 2019[Bournonville, S. de, Vangrunderbeeck, S. & Kerckhofs, G. (2019). Contrast Media Mol. Imaging, 2019, 1-9.]; Silva et al., 2015[Silva, J. M. S., Zanette, I., Noël, P. B., Cardoso, M. B., Kimm, M. A. & Pfeiffer, F. (2015). Sci. Rep. 5, 14088.]; Koç et al., 2019[Koç, M. M., Aslan, N., Kao, A. P. & Barber, A. H. (2019). Microsc. Res. Tech. 82, 812-848.]; Metscher, 2009[Metscher, B. D. (2009). BMC Physiol. 9, 11.]; O et al., 2021[O, J., Kwon, H. J., Cho, T. H., Woo, S. H., Rhee, Y. H. & Yang, H. M. (2021). PLoS One, 16, e0254264.]), to the use of other imaging modalities such as X-ray phase-contrast imaging (X-PCI) that uses the refractive properties of brilliant X-rays to enhance soft tissue image contrast (Bravin et al., 2013[Bravin, A., Coan, P. & Suortti, P. (2013). Phys. Med. Biol. 58, R1-R35.]). A variety of X-PCI methods including grating or crystal interferometry, analyzers or propagation-based (PB) are currently available (Bravin et al., 2013[Bravin, A., Coan, P. & Suortti, P. (2013). Phys. Med. Biol. 58, R1-R35.]). Among these methods, PB X-PCI is the simplest solution in terms of experimental setup and has proven to be very powerful for the structural examination of cardiac tissue (Reichardt et al., 2020[Reichardt, M., Töpperwien, M., Khan, A., Alves, F. & Salditt, T. (2020a). J. Med. Imag. 7, 1.]; Reichardt et al., 2020[Reichardt, M., Frohn, J., Khan, A., Alves, F. & Salditt, T. (2020b). Biomed. Opt. Expr. 11, 2633.]; Varray et al., 2017[Varray, F., Mirea, I., Langer, M., Peyrin, F., Fanton, L. & Magnin, I. E. (2017). Med. Image Anal. 38, 117-132.]; Mirea et al., 2015[Mirea, I., Varray, F., Zhu, Y. M., Fanton, L., Langer, M., Jouk, P. S., Michalowicz, G., Usson, Y. & Magnin, I. E. (2015). Functional Imaging and Modeling of the Heart, Vol. 9126 of Lecture Notes in Computer Science, pp. 172-179. Springer.]; Wang et al., 2019[Wang, S., Mirea, I., Varray, F., Liu, W. Y. & Magnin, I. E. (2019). Lecture Notes in Computer Science (Including Subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics). Springer.]; Gonzalez-Tendero et al., 2017[Gonzalez-Tendero, A., Zhang, C., Balicevic, V., Cardenes, R., Loncaric, S., Butakoff, C., Paun, B., Bonnin, A., Garcia-Cañadilla, P., Muñoz-Moreno, E., Gratacos, E., Crispi, F. & Bijnens, B. (2017). Eur Heart J. Cardiovasc. Imaging, 18, 732-741.]; Dejea et al., 2019[Dejea, H., Garcia-Canadilla, P., Cook, A. C., Guasch, E., Zamora, M., Crispi, F., Stampanoni, M., Bijnens, B. & Bonnin, A. (2019). Sci. Rep. 9, 6996.]; Soveral et al., 2020[Soveral, I., Crispi, F., Walter, C., Guirado, L., García-Cañadilla, P., Cook, A., Bonnin, A., Dejea, H., Rovira-Zurriaga, C., Sánchez de Toledo, J., Gratacós, E., Martínez, J. M., Bijnens, B. & Gómez, O. (2020). Ultrasound Obstet. Gyne, 56, 837-849.]; Loncaric et al., 2021[Loncaric, F., Garcia-Canadilla, P., Garcia-Alvarez, A., Sanchis, L., Prat, S., Doltra, A., Quintana, E., Pereda, D., Dejea, H., Bonnin, A., Sitges, M. & Bijnens, B. (2021). Front. Cardiovasc. Med. 8, 670734.]; Garcia-Canadilla et al., 2018[Garcia-Canadilla, P., Dejea, H., Bonnin, A., Balicevic, V., Loncaric, S., Zhang, C., Butakoff, C., Aguado-Sierra, J., Vázquez, M., Jackson, L. H., Stuckey, D. J., Rau, C., Stampanoni, M., Bijnens, B. & Cook, A. C. (2018). Circ. Cardiovasc. Imaging, 11, e007753.]; Planinc et al., 2021[Planinc, I., Garcia-Canadilla, P., Dejea, H., Ilic, I., Guasch, E., Zamora, M., Crispi, F., Stampanoni, M., Milicic, D., Bijnens, B., Bonnin, A. & Cikes, M. (2021). Sci. Rep. 11, 14020.]).

However, although soft tissue contrast is highly increased in X-PCI, several authors have proposed the use of ethanol to increase image contrast even further (Reichardt et al., 2020a[Reichardt, M., Töpperwien, M., Khan, A., Alves, F. & Salditt, T. (2020a). J. Med. Imag. 7, 1.]; Reichardt et al., 2020b[Reichardt, M., Frohn, J., Khan, A., Alves, F. & Salditt, T. (2020b). Biomed. Opt. Expr. 11, 2633.]; Varray et al., 2017[Varray, F., Mirea, I., Langer, M., Peyrin, F., Fanton, L. & Magnin, I. E. (2017). Med. Image Anal. 38, 117-132.]; Mirea et al., 2015[Mirea, I., Varray, F., Zhu, Y. M., Fanton, L., Langer, M., Jouk, P. S., Michalowicz, G., Usson, Y. & Magnin, I. E. (2015). Functional Imaging and Modeling of the Heart, Vol. 9126 of Lecture Notes in Computer Science, pp. 172-179. Springer.]; Patzelt et al., 2019[Patzelt, M., Mrzilkova, J., Dudak, J., Krejci, F., Zemlicka, J., Karch, J., Musil, V., Rosina, J., Sykora, V., Horehledova, B. & Zach, P. (2019). Jpn J. Radiol. 37, 500-510.]; Mrzílková et al., 2019[Mrzílková, J., Patzelt, M., Gallina, P., Wurst, Z., Šeremeta, M., Dudák, J., Krejčí, F., Žemlička, J., Musil, V., Karch, J., Rosina, J. & Zach, P. (2019). Biomed Res. Intl, 2019, 1-7.]; Dudak et al., 2016[Dudak, J., Zemlicka, J., Karch, J., Patzelt, M., Mrzilkova, J., Zach, P., Hermanova, Z., Kvacek, J. & Krejci, F. (2016). Sci. Rep. 6, 30385.]; Shirai et al., 2014[Shirai, R., Kunii, T., Yoneyama, A., Ooizumi, T., Maruyama, H., Lwin, T.-T., Hyodo, K. & Takeda, T. (2014). J. Synchrotron Rad. 21, 795-800.]; Takeda et al., 2014[Takeda, T., Yoneyama, A. & Hyodo, K. (2014). Photon Factory Activity Report 2014. Photon Factory, Tsukuba, Japan.]; Takeda et al., 2013a[Takeda, T., Lwin, T. T., Yoneyama, A. & Hyodo, K. (2013). Photon Factory Activity Report 2013. Photon Factory, Tsukuba, Japan.]; Lwin et al., 2022[Lwin, T.-T., Yoneyama, A., Kokubo, S., Maruyama, H., Hyodo, K. & Takeda, T. (2022). Acta Radiol. 63, 1102-1109.]; Kunii et al., 2013[Kunii, T., Shirai, R., Yoneyama, A., Maruyama, H., Lwin, T. T. & Takeda, T. (2013). Med. Imaging Technol. 31, 132-135.]; Takeda et al., 2013b[Takeda, T., Kunii, T., Sirai, R., Ohizumi, T., Maruyama, H., Hyodo, K., Yoneyama, A. & Ueda, K. (2013). J. Phys. Conf. Ser. 425, 022004.]). Ethanol is commonly used for tissue fixation (Patzelt et al., 2019[Patzelt, M., Mrzilkova, J., Dudak, J., Krejci, F., Zemlicka, J., Karch, J., Musil, V., Rosina, J., Sykora, V., Horehledova, B. & Zach, P. (2019). Jpn J. Radiol. 37, 500-510.]; Haque et al., 2020[Haque, Z., Rahman, A., Khan, Z. I., Hussan, M. T. & Alam, M. (2020). Int. J. Morphol. 38, 1371-1375.]; Essen et al., 2010[Essen, H. F. van, Verdaasdonk, M. A. M., Elshof, S. M., de Weger, R. A. & van Diest, P. J. (2010). J. Clin. Pathol. 63, 1090-1094.]) as well as a contrast agent both in conventional (tube-based) and synchrotron-based micro-CT imaging (Reichardt et al., 2020a[Reichardt, M., Töpperwien, M., Khan, A., Alves, F. & Salditt, T. (2020a). J. Med. Imag. 7, 1.]; Reichardt et al., 2020b[Reichardt, M., Frohn, J., Khan, A., Alves, F. & Salditt, T. (2020b). Biomed. Opt. Expr. 11, 2633.]; Varray et al., 2017[Varray, F., Mirea, I., Langer, M., Peyrin, F., Fanton, L. & Magnin, I. E. (2017). Med. Image Anal. 38, 117-132.]; Mirea et al., 2015[Mirea, I., Varray, F., Zhu, Y. M., Fanton, L., Langer, M., Jouk, P. S., Michalowicz, G., Usson, Y. & Magnin, I. E. (2015). Functional Imaging and Modeling of the Heart, Vol. 9126 of Lecture Notes in Computer Science, pp. 172-179. Springer.]; Shirai et al., 2014[Shirai, R., Kunii, T., Yoneyama, A., Ooizumi, T., Maruyama, H., Lwin, T.-T., Hyodo, K. & Takeda, T. (2014). J. Synchrotron Rad. 21, 795-800.]; Kunii et al., 2013[Kunii, T., Shirai, R., Yoneyama, A., Maruyama, H., Lwin, T. T. & Takeda, T. (2013). Med. Imaging Technol. 31, 132-135.]; Takeda et al., 2013b[Takeda, T., Kunii, T., Sirai, R., Ohizumi, T., Maruyama, H., Hyodo, K., Yoneyama, A. & Ueda, K. (2013). J. Phys. Conf. Ser. 425, 022004.]). However, despite the increase in image contrast, several studies have reported significant tissue shrinkage due to ethanol dehydration as well as small tissue ruptures due to fast dehydration (Patzelt et al., 2019[Patzelt, M., Mrzilkova, J., Dudak, J., Krejci, F., Zemlicka, J., Karch, J., Musil, V., Rosina, J., Sykora, V., Horehledova, B. & Zach, P. (2019). Jpn J. Radiol. 37, 500-510.]; Dudak et al., 2016[Dudak, J., Zemlicka, J., Karch, J., Patzelt, M., Mrzilkova, J., Zach, P., Hermanova, Z., Kvacek, J. & Krejci, F. (2016). Sci. Rep. 6, 30385.]; Hołda et al., 2016[Hołda, M. K., Klimek-Piotrowska, W., Koziej, M., Piątek, K. & Hołda, J. (2016). J. Anat. 229, 334-340.]; Hoshino et al., 2014[Hoshino, M., Uesugi, K., Tsukube, T. & Yagi, N. (2014). J. Synchrotron Rad. 21, 1347-1357.]; Buytaert et al., 2014[Buytaert, J., Goyens, J., De Greef, D., Aerts, P. & Dirckx, J. (2014). Microsc. Microanal. 20, 1208-1217.]; Vickerton et al., 2013[Vickerton, P., Jarvis, J. & Jeffery, N. (2013). J. Anat. 223, 185-193.]). Most of these studies have only quantified tissue shrinkage globally, without evaluating whether the tissue shrinks homogeneously or not, and how this affects the myocardial organization.

Cardiovascular disease is the leading cause of death in Europe, with over 3.9 million deaths (ESC, 2019[ESC (2019). About Cardiovascular Disease in ESC Member Countries, https://www.escardio.org/The-ESC/Press-Office/Fact-sheets.]). It involves cardiac remodelling at the cellular, tissue and organ level. Cardiomyocytes, which are the primary functional muscle cells in the myocardium, are aggregated and aligned in a pre-dominant direction within the myocardium, forming a complex 3D mesh. This complex organization of cardiomyocytes within the myocardium determines the propagation of electrical waves and force development within cardiac tissue. Therefore, a comprehensive structural assessment of cardiac remodelling at different length scales is needed in order to better understand how changes in cardiac structure affect cardiac function in cardiovascular disease. Synchrotron-based X-PCI has proved to be a powerful tool for the study of structural changes of hearts from rabbit fetuses (Gonzalez-Tendero et al., 2017[Gonzalez-Tendero, A., Zhang, C., Balicevic, V., Cardenes, R., Loncaric, S., Butakoff, C., Paun, B., Bonnin, A., Garcia-Cañadilla, P., Muñoz-Moreno, E., Gratacos, E., Crispi, F. & Bijnens, B. (2017). Eur Heart J. Cardiovasc. Imaging, 18, 732-741.]; Garcia-Canadilla et al., 2019[Garcia-Canadilla, P., de Vries, T., Gonzalez-Tendero, A., Bonnin, A., Gratacos, E., Crispi, F., Bijnens, B. & Zhang, C. (2019). PLoS One, 14, e0218192.]), young rats (Planinc et al., 2021[Planinc, I., Garcia-Canadilla, P., Dejea, H., Ilic, I., Guasch, E., Zamora, M., Crispi, F., Stampanoni, M., Milicic, D., Bijnens, B., Bonnin, A. & Cikes, M. (2021). Sci. Rep. 11, 14020.]; Dejea et al., 2019[Dejea, H., Garcia-Canadilla, P., Cook, A. C., Guasch, E., Zamora, M., Crispi, F., Stampanoni, M., Bijnens, B. & Bonnin, A. (2019). Sci. Rep. 9, 6996.]), human fetuses (Garcia-Canadilla et al., 2018[Garcia-Canadilla, P., Dejea, H., Bonnin, A., Balicevic, V., Loncaric, S., Zhang, C., Butakoff, C., Aguado-Sierra, J., Vázquez, M., Jackson, L. H., Stuckey, D. J., Rau, C., Stampanoni, M., Bijnens, B. & Cook, A. C. (2018). Circ. Cardiovasc. Imaging, 11, e007753.]) as well as human endomyocardial biopsies (Planinc et al., 2023[Planinc, I., Ilic, I., Dejea, H., Garcia-Canadilla, P., Gasparovic, H., Jurin, H., Milicic, D., Skoric, B., Stampanoni, M., Bijnens, B., Bonnin, A. & Cikes, M. (2023). Transpl. Int. 36, 11046.]; Loncaric et al., 2021[Loncaric, F., Garcia-Canadilla, P., Garcia-Alvarez, A., Sanchis, L., Prat, S., Doltra, A., Quintana, E., Pereda, D., Dejea, H., Bonnin, A., Sitges, M. & Bijnens, B. (2021). Front. Cardiovasc. Med. 8, 670734.]) at micrometre resolution with sufficient contrast to distinguish myocytes at the cellular level and without the need for destructive sectioning or tissue processing. However, in order to properly quantify myocardial microstructure, especially cardiomyocytes aggregates (often referred to as myofibers and sheets/sheetlets), it is crucial to minimize structural changes in cardiac tissue as a consequence of the imaging process.

The aim of this study is to evaluate the local myocardial tissue deformation due to ethanol dehydration based on 3D non-rigid registration of synchrotron-based X-PCI images of a rat heart and detailed characterization of its myocardial tissue organization, before and after ethanol dehydration.

2. Materials and methods

2.1. Sample preparation

One 1.5-month-old male Wistar rat was obtained from Charles River Laboratories. The rat was housed individually and maintained at 21°C with a 12 h day/night cycle. Food and water were administrated ad libitum. The rat was sacrificed after two weeks, the heart was rapidly excised, rinsed in PBS + 2% heparin solution and immersed in 4% paraformaldehyde. Animal care and experimentation conformed to the European Union (Directive 2010/63/UE) and Spanish guidelines (RD 53/2013) for the use of experimental animals. Approval was obtained from the local animal research ethics committee `Comité de Ética de Experimentación Animal (CEEA)' (CEEA 68/5435).

For the image acquisition, the heart was placed in a cylindrical plastic holder filled with degassed deionized water to hold the heart in place without compression and avoid motion artefacts during tomography acquisition. The heart was first scanned with synchrotron PB X-PCI without further sample preparation. Then, the heart was immersed in successive baths (lasting 30 min) increasing the ethanol concentration from 10% to 70% in 10% steps as described by Varray et al. (2017[Varray, F., Mirea, I., Langer, M., Peyrin, F., Fanton, L. & Magnin, I. E. (2017). Med. Image Anal. 38, 117-132.], 2013[Cibc (2013). Seg3D: Volumetric Image Segmentation and Visualization. Scientific Computing and Imaging Institute (SCI), https://www.mendeley.com/catalogue/a702d9d1-b36d-38f0-9c4b-44f7c552a1b2/?utm_source=desktop&utm_medium=1.17.9&utm_campaign=open_catalog&userDocumentId=%7Bc869e468-ef14-4cd4-988b-46570e2.]). After that, the heart was scanned again. Then the heart was kept in 70% ethanol for 12 days and subsequently the third and last scan was performed. Therefore, the sample was scanned three times: before, 9 h after and 342 h after ethanol immersion. All three scans were performed with the same PB X-PCI setup and acquisition parameters.

2.2. Image acquisition

The heart was scanned with synchrotron PB X-PCI at the TOMCAT beamline (X02DA) of the Swiss Light Source (Paul Scherrer Institute, Switzerland). The plastic tube containing the heart was placed on the rotation stage located 333 cm from the detector. A 20 keV parallel synchrotron X-ray beam (monochromaticity bandwidth of 2%) was used to image the sample with 5.8 µm effective pixel size and field of view (FoV) of 11.83 mm × 3.29 mm. Acquisition was performed by rotating the sample through 360° acquiring 2501 projections (exposure time: 20 ms). X-rays were converted to visible light through a LuAG:Ce 300 µm scintillator and detected by an sCMOS camera (PCO.Edge 4.2). Since the heart was larger than the FoV, several overlapping scans from base to apex were performed to fully cover the full height. Additionally, two series of 50 flat and 20 dark images were acquired for flat-field and dark-field corrections of each acquisition. The acquired projections were reconstructed using the Gridrec algorithm (Marone & Stampanoni, 2012[Marone, F. & Stampanoni, M. (2012). J. Synchrotron Rad. 19, 1029-1037.]). All datasets were reconstructed with the same contrast limits. Finally, the reconstructed datasets were stitched in order to obtain a single full heart dataset.

2.3. Image analysis

2.3.1. Image alignment

A multi-resolution rigid registration between all image datasets was performed in order to correct for differences in sample position during scanning. Briefly, the rigid transformation that minimized the Matte's mutual information (MMI) (Mattes et al., 2001[Mattes, D., Haynor, D. R., Vesselle, H., Lewellyn, T. K. & Eubank, W. (2001). Proc. SPIE, 4322, 1609-1620.]) between moving (9 h and 342 h after ethanol immersion) and reference (before ethanol) image datasets was found. To do that, a multi-resolution framework consisting of three stages was used, in which the resolution of the original image was successively reduced with respective subsampling factors of 4, 8 and 16. The minimization problem was solved via the limited Broyden–Fletcher–Goldfarb–Shanno algorithm (L-BFGS) optimization method (Mattes et al., 2001[Mattes, D., Haynor, D. R., Vesselle, H., Lewellyn, T. K. & Eubank, W. (2001). Proc. SPIE, 4322, 1609-1620.]). Therefore, both after ethanol immersion datasets were transformed to be aligned to the reference dataset.

2.3.2. Cardiac tissue segmentation

Cardiac tissue was semi-automatically segmented using the two-stage pixel classification module from the open-source software Ilastik (Berg et al., 2019[Berg, S., Kutra, D., Kroeger, T., Straehle, C. N., Kausler, B. X., Haubold, C., Schiegg, M., Ales, J., Beier, T., Rudy, M., Eren, K., Cervantes, J. I., Xu, B., Beuttenmueller, F., Wolny, A., Zhang, C., Koethe, U., Hamprecht, F. A. & Kreshuk, A. (2019). Nat. Methods, 16, 1226-1232.]; Sommer et al., 2011[Sommer, C., Sträehle, C., Köthe, U. & Hamprecht, F. A. (2011). 2011 IEEE International  Symposium on Biomedical Imaging, pp. 230-233. IEEE.]). Then, different morphological operations were applied to smooth the segmentation using an in-house algorithm implemented in MATLAB. Finally, blood clots inside the ventricular cavities were manually removed using Seg3D software (Cibc, 2013[Cibc (2013). Seg3D: Volumetric Image Segmentation and Visualization. Scientific Computing and Imaging Institute (SCI), https://www.mendeley.com/catalogue/a702d9d1-b36d-38f0-9c4b-44f7c552a1b2/?utm_source=desktop&utm_medium=1.17.9&utm_campaign=open_catalog&userDocumentId=%7Bc869e468-ef14-4cd4-988b-46570e2.]).

2.3.3. Quantification of image contrast and quality

The parameters used to evaluate the differences in image contrast and quality due to ethanol dehydration were the tissue contrast (C), the contrast-to-noise ratio (CNR) (Brombal et al., 2019[Brombal, L., Arfelli, F., Donato, S., Rigon, L., Golosio, B., Oliva, P., Longo, R., Bonazza, D., Zanconati, F., Contillo, A., Taibi, A., Delogu, P., Mettivier, G. & Tromba, G. (2019). J. Med. Imaging, 6, 031402.], 2018[Brombal, L., Donato, S., Dreossi, D., Arfelli, F., Bonazza, D., Contillo, A., Delogu, P., Di Trapani, V., Golosio, B., Mettivier, G., Oliva, P., Rigon, L., Taibi, A. & Longo, R. (2018). Phys. Med. Biol. 63, aaf2e1.]) and the signal-to-noise ratio (SNR) as follows,

[C\ \left(\%\right) = {{{I_{\rm{T}}}-{I_{\rm{B}}}}\over{{I_{\rm{B}}}}}\times100, \eqno(1)]

[{\rm{CNR}} = {{{I_{\rm{T}}}-{I_{\rm{B}}}}\over{\left[\left(\sigma_{\rm{T}}^2+\sigma_{\rm{B}}^2\right)/2\right]^{1/2}}}, \eqno(2)]

[{\rm{SNR}}\ \left({\rm{dB}}\right) = 20\log\left(I_{\rm{T}}/\sigma_{\rm{B}}\right), \eqno(3)]

where subscripts T and B denote the cardiac tissue and background, respectively, and Ii and σi are the mean and standard deviation voxel intensity values, respectively. These parameters were computed in three different image slices from apex to base (see Fig. 1[link]) for each dataset. Finally, distributions of the voxel intensity values within the cardiac tissue were also computed.

[Figure 1]
Figure 1
Demonstration of the contrast improvement in synchrotron X-PCI and sample shrinkage due to ethanol dehydration between (a) native heart kept in formalin, (b) 9 h after and (b) 342 h after immersion in 70% ethanol solution. (d) Distribution of pixel intensities in the three datasets. All the datasets were reconstructed using the same contrast limits. Dashed lines in panel (a) indicate the three apico-basal slices selected for quantification.
2.3.4. Image non-rigid registration

To evaluate and quantify the tissue shrinkage produced by ethanol dehydration and study whether the deformation was uniform, a 3D non-rigid registration between the before (fixed) and 342 h after ethanol (moving) datasets was performed to obtain the local deformation field that warped the moving towards the fixed dataset. The moving dataset was previously aligned to the fixed images using the rigid registration described previously.

The symmetric force additive demons algorithm, which is a variation of the classical demons algorithm in which the forces are symmetric (Wang et al., 2005[Wang, H., Dong, L., O'Daniel, J., Mohan, R., Garden, A. S., Ang, K. K., Kuban, D. A., Bonnen, M., Chang, J. Y. & Cheung, R. (2005). Phys. Med. Biol. 50, 2887-2905.]; Rogelj & Kovačič, 2006[Rogelj, P. & Kovačič, S. (2006). Med. Image Anal. 10, 484-493.]), was used for computing the transformation that wraps the moving to fixed image. The demons algorithm assumes that the intensity of homologous points of moving and fixed images are equal. However, this assumption is not valid since the image intensity changes among the different datasets due to the use of ethanol. Therefore, we used histogram equalization prior to the non-rigid registration in order to create a more exact match between the intensities of the two datasets. Due to the large size of the original datasets (4080 × 4080 × 3747) we used a multi-resolution pyramids approach, in which the resolution of the original image was successively reduced, consisting of three stages with respective down-sampling factors of 4, 8 and 16. The implementation of the symmetric demons algorithm available in Simple ITK was used (Lowekamp et al., 2013[Lowekamp, B. C., Chen, D. T., Ibáñez, L. & Blezek, D. (2013). Front. Neuroinform. 7, 73735.]; Yaniv et al., 2018[Yaniv, Z., Lowekamp, B. C., Johnson, H. J. & Beare, R. (2018). J. Digit. Imaging, 31, 290-303.]; Beare et al., 2018[Beare, R., Lowekamp, B. & Yaniv, Z. (2018). J. Stat. Soft. 86, 1-35.]). A more detailed explanation of the non-rigid registration method can be found in Section S1 of the supporting information.

Full resolution non-rigid registration in a region of interest. In order to study the local tissue shrinkage produced by ethanol dehydration in more detail, a region of interest (ROI) of 150 × 150 × 150 pixels within the left ventricle (LV) was manually selected, in which the non-rigid registration was refined. To do so, a second non-rigid registration between the fixed (before ethanol) and moving (342 h after ethanol non-rigid registered) ROIs at full resolution was carried out using the same symmetric demons algorithm. The size of the ROI was chosen in order to be able to quantify changes in cardiac tissue at full resolution without prohibitive computational cost.

Evaluation of registration performance. The Dice similarity coefficient between the reference and registered datasets was computed to quantitatively evaluate the quality of the registrations. The Dice coefficient was calculated as

[{\rm{Dice}} = {{ 2\,\left|{X \cap Y}\right| }\over{ \left|X\right|+\left|Y\right| }}, \eqno(4)]

where |XY| represents the common voxels between both datasets, and |X| and |Y| represent the number of voxels in X and Y, respectively.

2.3.5. Quantification of global and local tissue deformation

To quantify global tissue shrinkage due to ethanol dehydration, the myocardial volume of all three aligned datasets was computed. Then, the LV wall thickness was automatically measured in the anterior, posterior septal and lateral walls in three apico-basal image slices (see Fig. 1[link]).

Finally, the non-rigid registration was then used to recover the deformation fields and quantify the local tissue shrinkage produced by ethanol dehydration. To do that, the Jacobian determinant, JD, of the deformation field was calculated. In particular, JD measures how the voxel volume changes after non-rigid registration, where JD > 1.0 denotes volume dilation while JD < 1.0 means volume shrinkage.

2.3.6. Analysis of myocardial organization

Myocardial organization was assessed by means of a structure tensor (ST) based method implemented in MATLAB as described in the literature (Gonzalez-Tendero et al., 2017[Gonzalez-Tendero, A., Zhang, C., Balicevic, V., Cardenes, R., Loncaric, S., Butakoff, C., Paun, B., Bonnin, A., Garcia-Cañadilla, P., Muñoz-Moreno, E., Gratacos, E., Crispi, F. & Bijnens, B. (2017). Eur Heart J. Cardiovasc. Imaging, 18, 732-741.]; Garcia-Canadilla et al., 2018[Garcia-Canadilla, P., Dejea, H., Bonnin, A., Balicevic, V., Loncaric, S., Zhang, C., Butakoff, C., Aguado-Sierra, J., Vázquez, M., Jackson, L. H., Stuckey, D. J., Rau, C., Stampanoni, M., Bijnens, B. & Cook, A. C. (2018). Circ. Cardiovasc. Imaging, 11, e007753.]; Baličević et al., 2015[Baličević, V., Lončarić, S., Cárdenes, R., Gonzalez-Tendero, A., Paun, B., Crispi, F., Butakoff, C. & Bijnens, B. (2015). Functional Imaging and Modeling of the Heart, Vol. 9126 of Lecture Notes in Computer Science, pp. 111-119. Springer.]), in all three datasets: before, 9 h after and 342 h after ethanol immersion. Briefly, the image intensity gradient in the x, y and z directions was obtained for each voxel using a central difference algorithm within the voxel's cubical neighbourhood. The ST was then calculated as the cross-product of gradient vectors. Then, eigen-decomposition of the ST was performed thus obtaining the three eigenvalues (λi) and eigenvectors (vi).

Different parameters based on the relative magnitude of the eigenvalues were calculated to quantify the anisotropy of myocardial tissue. In a manner analogous to diffusion tensor imaging (DTI), fractional anisotropy of the structure tensor, FAST, is defined as

[{\rm{FA}}_{\rm{ST}} = \left({{3\over2}}\right)^{\!\!1/2} \, {{ \left[ \big(\lambda_1-{\hat{\,\,\lambda}}\big)^2 + \big(\lambda_2-{\hat{\,\,\lambda}}\big)^2 + \big(\lambda_3-{\hat{\,\,\lambda}}\big)^2\right]^{1/2} }\over{ \left(\lambda_1^2+\lambda_2^2+\lambda_3^2\right)^{1/2} }}, \eqno(5)]

where [\hat{\,\,\lambda}] = [(\lambda_1+\lambda_2+\lambda_3)/3] is the average of the tensor's eigenvalues. In the same manner, the spherical (Cs), linear (Cl) and planar anisotropy (Cp) are defined (Nelson et al., 2018[Nelson, L. J., Smith, R. A. & Mienczakowski, M. (2018). Composites Part A, 104, 108-119.]),

[C_{\rm{s}} = {{\lambda_3}\over{\lambda_1}}, \qquad C_{\rm{l}} = {{\lambda_2-\lambda_3}\over{\lambda_1}}, \qquad C_{\rm{p}} = {{\lambda_1-\lambda_2}\over{\lambda_1}}, \eqno(6)]

with Cs + Cl + Cp = 1, and 0 ≤ Cs, Cl, Cp ≤ 1.

Then, the eigenvector with the smallest eigenvalue (λ3), known as the tertiary eigenvector (v3), was selected as the vector pointing in the main direction of the myocyte aggregates. In order to assess the orientation of the myocyte aggregates, its helical angle (HA), defined as the angle between the tertiary eigenvector v3 and the transverse plane defined by the local circumferential and radial directions of the cylindrical coordinate system of the heart, was calculated. All the quantitative parameters (FAST, Cs, Cl, Cp and HA) were calculated in the rigidly aligned datasets.

Quantification of myocyte aggregates orientation. The myocyte aggregates orientation was quantified in the same three apico-basal image slices selected to quantify image contrast and quality (see Fig. 1[link]). In order to quantify the changes in myocardial organization induced by ethanol dehydration, the calculated deformation fields of the non-rigid registrations between the after and before ethanol immersion datasets obtained in Section 2.3.4[link] were applied to all parameters calculated (FAST, Cs, Cl, Cp and HA) of the moving datasets. Distributions of FAST, Cs, Cl, Cp and HA within the cardiac tissue at each apico-basal image slice were computed.

Finally, transmural profiles of HA from endo- to epicardium were obtained for four different ROIs within the LV: anterior, lateral, posterior and septal walls, and for each of the three apico-basal image slices. Next, a linear regression fitting, y = β1x + β0, was applied to all transmural HA profiles. The R2 coefficient, of the linear fitting, which characterizes the linearity of the profiles, was also computed.

2.4. Statistical analysis

Normality of data was tested with a Shapiro–Wilks test, adequate for a small number of samples. Normally distributed variables were expressed as mean ± standard deviation while non-normally distributed continuous variables were expressed as medians with interquartile ranges.

3. Results

3.1. Quantification of image contrast and quality

After 9 h of ethanol treatment, the tissue contrast (C) was slightly improved, from 1.08% to 2.80%, also demonstrated by the increase in CNR (see Table 1[link]) without significant changes in image quality. After 342 h of ethanol immersion, the tissue contrast was further improved up to 3.88%, as clearly noticed also in the raw X-PCI images [Figs. 1[link](a)–1(c)], which is also demonstrated by the wider distribution of the voxel intensities [see Fig. 1[link](d)]. However, the image quality was slightly reduced by 3 dB, explained by the higher amount of image artefacts resulting in an increase of the background's standard deviation (σB).

Table 1
Quantification of changes in cardiac tissue contrast and image quality between all three datasets: before, after 9 h and after 342 h of ethanol treatment

CNR: contrast to noise ratio; C: tissue contrast; SNR: signal-to-noise ratio

  Before After 9 h After 342 h
CNR 0.097 ± 0.028 0.195 ± 0.023 0.127 ± 0.015
C (%) 1.084 ± 0.178 2.804 ± 0.094 3.883 ± 0.251
SNR (dB) 27.17 ± 2.85 28.91 ± 1.17 24.15 ± 0.63

3.2. Quantification of global and local tissue shrinkage

We found that long-term immersion (342 h) of the heart in ethanol led to a significant reduction of myocardial volume of −27.7% (from 0.75 cm3 to 0.54 cm3) due to dehydration (Table 2[link], Fig. 2[link]). Moreover, long-term immersion in ethanol also resulted in a heterogeneous reduction of LV myocardial thickness of about −18.33% on average, being more pronounced in the basal and mid-anterior and lateral walls (Table 2[link]).

Table 2
Myocardial volume and left ventricular (LV) wall thickness measurements

  Before After 9 h After 342 h
Value Change (%) Value Change (%)
Myocardial volume (cm3) 0.75 0.75 0 0.54 −27.69
LV wall thickness (mm)
 Basal anterior 3.19 3.15 1.30 2.34 −26.88
 Basal septal 3.16 3.16 0.05 2.53 −19.98
 Basal posterior 3.84 3.84 −0.03 3.13 −18.40
 Basal lateral 3.10 3.05 1.69 2.27 −26.79
 Mid-anterior 4.51 4.49 0.37 3.27 −27.48
 Mid-septal 2.91 2.94 -0.78 2.56 −12.13
 Mid-posterior 2.53 2.52 0.53 2.31 −8.90
 Mid-lateral 3.82 3.78 0.89 2.99 −21.50
 Apical anterior 2.77 2.77 −0.05 2.30 −16.95
 Apical septal 3.07 3.07 0.03 2.80 −8.79
 Apical posterior 2.44 2.43 0.41 1.98 −18.85
 Apical lateral 2.80 2.79 0.28 2.43 −13.34
Mean ± STD 3.18 ± 0.60 3.17 ± 0.60 −0.39 ± 0.66 2.58 ± 0.39 −18.33 ± 6.65
[Figure 2]
Figure 2
Top: three-dimensional reconstruction of the rat heart before and 342 h after 70% ethanol immersion rigidly aligned. Bottom: 2D resliced images along the three orthogonal planes x (blue), y (green) and z (red), clearly showing the good alignment between both datasets and the volume shrinkage of the heart after long-term immersion in ethanol.

Regarding the local deformation, results of the non-rigid registration between the before and 342 h after ethanol immersion datasets can be found in Figs. S1 and S2. A Dice similarity coefficient of 0.97 demonstrates the good matching between reference and registered datasets, which is also illustrated in Fig. S1. When performing a finer registration in a ROI (Fig. S2), the Dice coefficient was 0.99, thus demonstrating an improvement of the previous coarser non-rigid registration.

Looking at Fig. 3[link] one can realize that local deformation induced by ethanol dehydration is not uniform, since there are areas with a high degree of shrinkage (dark blue) and areas where tissue shrinkage is very low (turquoise/green). Values higher than 100% (expansion, since 100% corresponds to no volume changes) are observed along the spaces between myocardium and cavities, indicating that cavity space is increased due to tissue contraction.

[Figure 3]
Figure 3
Quantification of the local myocardial volume changes due to ethanol dehydration. (a) Colour maps of the local Jacobian determinants of the deformation field, JD, at three levels from base to apex, obtained after the non-rigid registration between the heart before and 342 h after ethanol immersion. (b) Local JD at three different image slices obtained after the non-rigid registration at full resolution within the ROI indicated with a dashed red rectangle. Reddish values indicated tissue expansion while blueish values indicated tissue shrinkage.

3.3. Quantification of myocardial organization

Figs. 4[link] and 5[link] show colour maps of the local values of FAST and HA, respectively, for all three datasets: before, 9 h after and 342 h after ethanol immersion. Regarding FAST, there is a significant increase of about 0.16 due to tissue dehydration 342 h after ethanol immersion, while differences are minimal between the heart before and 9 h after ethanol immersion (Fig. 4[link]).

[Figure 4]
Figure 4
Quantification of myocardial organization via fractional anisotropy (FA). Local FA at three levels from base to apex of the rat heart (a) before, (b) 9 h after and (c) 342 h after ethanol immersion. Difference in local FA between the heart before and (d) 9 h after and (e) 342 h after ethanol immersion. The average FA difference ([\bar\varepsilon]) was calculated for each image slice.
[Figure 5]
Figure 5
Quantification of local helical angle (HA) at three levels from base to apex of the rat heart (a) before, (b) 9 h after and (c) 342 h after ethanol immersion. Difference in local HA between the heart before and 9 h after (d) and 342 h after (e) ethanol immersion. The average HA difference ([\bar\varepsilon]) was calculated for each image slice.

For HA, it can be observed that, while there are only small differences between the heart before and 9 h after ethanol immersion, the differences are much bigger (the average HA difference is about 18°; Fig. 5[link]) between the heart before and 342 h after ethanol immersion. There are no differences neither in the HA transmural profiles nor HA distributions between the heart before and 9 h after ethanol immersion (Fig. 6[link]; Table S1). However, while the slope and linearity coefficient of the transmural profile of HA remain unchanged, a significant decrease in the intercept of the transmural profile of HA is noticed in the heart 342 h after ethanol immersion (41.18 ± 9.59 versus 36.93 ± 9.17; p = 0.0015), being more pronounced in the basal septal and apical anterior walls (see Table 3[link] and Fig. 5[link]).

Table 3
Parameters of the linear fitting (y = β1x + β0) of the HA transmural profile within the LV rat heart before and after 342 h of ethanol dehydration

  Slope: β1 (°) Intercept: β0 (°) Linearity (R2)
  Before After 342 h Before After 342 h Before After 342 h
Basal anterior −75.45 −69.27 33.85 31.20 0.97 0.96
Basal septal −94.87 −92.09 46.77 30.21 0.99 0.95
Basal posterior −40.22 −33.67 33.92 31.91 0.83 0.85
Basal lateral −69.48 −59.44 28.34 21.58 0.99 0.98
Mid-anterior −54.87 −57.04 31.64 28.36 0.94 0.91
Mid-septal −97.32 −96.79 43.62 39.51 0.95 0.95
Mid-posterior −40.91 −36.87 36.83 36.88 0.96 0.94
Mid-lateral −69.75 −70.97 36.49 34.63 0.93 0.98
Apical anterior −66.13 −61.68 62.81 48.34 0.92 0.89
Apical septal −82.47 −93.52 47.85 53.45 0.97 0.97
Apical posterior −48.53 −55.82 43.50 45.78 0.98 0.97
Apical lateral −86.29 −84.41 48.59 41.34 0.93 0.96
Mean ± STD −68.8 ± 19.6 −67.63 ± 21.03 41.18 ± 9.59 36.93 ± 9.17 0.94 ± 0.04 0.94 ± 0.04
[Figure 6]
Figure 6
Quantification of myocytes aggregates orientation via helical angle (HA). (a) Histograms of HA in the left ventricle (LV) and in the right ventricle (RV) across the basal, mid-ventricular and apical slices of the rat heart before (blue), 9 h after (red) and 342 h after (green) ethanol immersion. (b) LV transmural profiles of HA across four LV segments (anterior septal, posterior and lateral) in the same three apico-basal slices of the rat heart before (blue), 9 h after (red) and 342 h after (green) ethanol immersion. Solid lines: sample mean. Shadows: ± standard deviation.

Finally, while distributions of Cp, Cs and Cl of the heart before and 9 h after ethanol immersion are very similar, an increase in planar anisotropy Cp together with the subsequent decrease in Cs and Cl is observed in the heart 342 h after ethanol immersion (see Fig. S3). This can be explained by the fact that, since the ventricular intramyocardial space between layers of myocytes, also known as cleavage planes, is increased due to tissue shrinkage, the laminar arrangement of myocytes becomes more evident.

4. Discussion

In the present study we have quantitatively shown that the use of ethanol dehydration in synchrotron X-PCI can improve image contrast at the expense of inhomogeneous tissue shrinkage, which has been shown to disrupt local myocardial organization. This has been demonstrated by the changes in the quantitative parameters commonly used to quantify myocardial organization such as HA and FA.

To our knowledge this is the first study that performs a detailed quantitative analysis of the local myocardial changes produced by ethanol dehydration. First, a rat heart was imaged with synchrotron-based X-PCI, before, 9 h after and 342 h after its immersion in 70% ethanol solution. Then, we quantified image contrast, global volume and myocardial thickness changes due to ethanol dehydration. In order to determine local deformation of the heart caused by ethanol dehydration, the related 3D tomographic datasets were registered by means of a non-rigid registration algorithm. The non-rigid registration was then used to recover the deformation fields, and the local relative volume change was computed. Finally, myocardial organization parameters such as FA and the orientation of myocytes aggregates were also computed and compared in the three tomographic datasets.

The global volume shrinkage was found to be 27.7% after long-term immersion of the heart in 70% ethanol. Moreover, we found a heterogeneous reduction of LV myocardial thickness of about 18.3 ± 6.7% on average, with maximum values of about 27% in the basal anterior and lateral and mid-anterior walls. Similar results were found by Patzelt et al. (2019[Patzelt, M., Mrzilkova, J., Dudak, J., Krejci, F., Zemlicka, J., Karch, J., Musil, V., Rosina, J., Sykora, V., Horehledova, B. & Zach, P. (2019). Jpn J. Radiol. 37, 500-510.]) who evaluated three ethanol fixation protocols for heart and lung micro-CT imaging including fixed ethanol concentrations of 50% and 97% and four consequent ethanol baths of increasing concentrations. They reported a 21% volume shrinkage, high sample stiffness and small tissue ruptures due to fast dehydration in the heart after 168 h immersion in 97% ethanol. The authors concluded that fixation in a series of ascending ethanol concentrations provided the best results in terms of contrast enhancement and tissue preservation avoiding tissue ruptures. Similarly, Dudak et al. (2016[Dudak, J., Zemlicka, J., Karch, J., Patzelt, M., Mrzilkova, J., Zach, P., Hermanova, Z., Kvacek, J. & Krejci, F. (2016). Sci. Rep. 6, 30385.]) also tested three different ethanol fixation techniques (50%, 97% and series of increasing concentration 50–97%) in soft tissue murine organs and found very similar results in terms of contrast enhancement. However, they did not quantify the volume changes due to tissue dehydration. Holda et al. (2016[Hołda, M. K., Klimek-Piotrowska, W., Koziej, M., Piątek, K. & Hołda, J. (2016). J. Anat. 229, 334-340.]) investigated the effect of different types of fixatives and concentrations on the dimensions of cardiac tissue and found that, after 168 h immersion in 85.5% ethanol, heart weight was significantly reduced by 16%. They also found a statistically significant but heterogeneous myocardial wall thickness reduction of up to 16% in the RV. On the other hand, Buytaer et al. (2014[Buytaert, J., Goyens, J., De Greef, D., Aerts, P. & Dirckx, J. (2014). Microsc. Microanal. 20, 1208-1217.]) quantified the total volume shrinkage of three tissue types: bone, muscle and nerve tissue, caused by four common specimen preparations: two water-based and two ethanol-based methods using the same dehydrating alcohol series (30, 60, 90, 100 and 100%, each for at least a day). To do that, they scanned the samples before (fresh) and after specimen preparation for light sheet fluorescence microscopy (LSFM) or micro-CT. They found that LSFM preparation of the samples, regardless of the tissue type, caused the most volume shrinkage compared with the other methods, with an average volume shrinkage for muscle of 56 ± 2%.

Most of the previous studies evaluated the global volume shrinkage produced by ethanol dehydration. The present study, however, determined the local volume changes induced by ethanol within the whole heart by means of a non-rigid registration method. The most important advantage of this technique is the possibility to analyse local shrinkage and deformations of different parts of the heart. Our results clearly show that shrinkage is anisotropic with areas with a high degree of shrinkage and areas where shrinkage is very low, or even where the tissue expands as illustrated in Fig. 3[link]. The maximum expansion values were observed along the spaces between myocardium and cavities, explained by the increase of the cavity space due to tissue contraction. Schulz et al. (2011[Schulz, G., Crooijmans, H. J. A., Germann, M., Scheffler, K., Müller-Gerbl, M. & Müller, B. (2011). J. Neurosci. Methods, 202, 17-27.]) used a similar approach to quantify the local deformation produced by formalin fixation in the brain. They found maximal local volume strain values of 22% and −24% due to formalin fixation. The maximal expansion values were mainly located at the ventricles which can be explained by the enlargement of the ventricles due to the lower pressure of the brain outside the cranium.

Myocardial walls are formed by myocytes that are aggregated together forming a complex 3D network surrounded by a matrix of connective tissue, showing predominant directions within the ventricular myocardium. The complex organization of these aggregates of myocytes within the myocardium determines the propagation of electrical waves as well as the force development within cardiac tissue. However, there is still some controversy regarding the fashion in which the individual cardiomyocytes are aggregated together (Anderson et al., 2019[Anderson, R. H., Niederer, P. F., Sanchez-Quintana, D., Stephenson, R. S. & Agger, P. (2019). J. Anat. 235, 697-705.]). Several studies have suggested that myocytes are also grouped forming flattened structures of approximately four myocytes thick named: `sheets', `sheetlets', `laminar' or `lamellar' (Varray et al., 2017[Varray, F., Mirea, I., Langer, M., Peyrin, F., Fanton, L. & Magnin, I. E. (2017). Med. Image Anal. 38, 117-132.]; Mirea et al., 2015[Mirea, I., Varray, F., Zhu, Y. M., Fanton, L., Langer, M., Jouk, P. S., Michalowicz, G., Usson, Y. & Magnin, I. E. (2015). Functional Imaging and Modeling of the Heart, Vol. 9126 of Lecture Notes in Computer Science, pp. 172-179. Springer.]; Wang et al., 2019[Wang, S., Mirea, I., Varray, F., Liu, W. Y. & Magnin, I. E. (2019). Lecture Notes in Computer Science (Including Subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics). Springer.]; Varray et al., 2013[Varray, F., Wang, L., Fanton, L., Zhu, Y.-M. & Magnin, I. E. (2013). Functional Imaging and Modeling of the Heart, Vol. 7945 of Lecture Notes in Computer Science, pp. 150-157. Berlin, Heidelberg: Springer.]; LeGrice et al., 1995[LeGrice, I. J., Smaill, B. H., Chai, L. Z., Edgar, S. G., Gavin, J. B. & Hunter, P. J. (1995). Am. J. Physiol. Heart Circ. Physiol. 269, H571-H582.], 2005[LeGrice, I., Pope, A. & Smaill, B. (2005). Interstitial Fibrosis in Heart Failure, pp. 3-21. Springer.]; Gilbert et al., 2012[Gilbert, S. H., Benoist, D., Benson, A. P., White, E., Tanner, S. F., Holden, A. V., Dobrzynski, H., Bernus, O. & Radjenovic, A. (2012). Am. J. Physiol. Heart Circ. Physiol. 302, H287-H298.]). In this work we have shown that cleavage planes are more evident after long-term immersion of the heart in 70% ethanol, especially in the papillary muscles and endocardial surface as illustrated in Fig. 1[link] and Figs. S1 and S2. This is also illustrated by the increase in the planar anisotropy Cp, and decrease in the spherical anisotropy Cs of the heart after long-term immersion in ethanol. Interestingly, most of the works that have described the existence of laminar structures within the myocardium have used ethanol as a fixative solution, as part of the tissue processing and/or to improve image contrast, or other contrast agents that also cause tissue shrinkage (Mirea et al., 2015[Mirea, I., Varray, F., Zhu, Y. M., Fanton, L., Langer, M., Jouk, P. S., Michalowicz, G., Usson, Y. & Magnin, I. E. (2015). Functional Imaging and Modeling of the Heart, Vol. 9126 of Lecture Notes in Computer Science, pp. 172-179. Springer.]; Wang et al., 2019[Wang, S., Mirea, I., Varray, F., Liu, W. Y. & Magnin, I. E. (2019). Lecture Notes in Computer Science (Including Subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics). Springer.]; Varray et al., 2017[Varray, F., Mirea, I., Langer, M., Peyrin, F., Fanton, L. & Magnin, I. E. (2017). Med. Image Anal. 38, 117-132.], 2013[Cibc (2013). Seg3D: Volumetric Image Segmentation and Visualization. Scientific Computing and Imaging Institute (SCI), https://www.mendeley.com/catalogue/a702d9d1-b36d-38f0-9c4b-44f7c552a1b2/?utm_source=desktop&utm_medium=1.17.9&utm_campaign=open_catalog&userDocumentId=%7Bc869e468-ef14-4cd4-988b-46570e2.]; LeGrice et al., 1995[LeGrice, I. J., Smaill, B. H., Chai, L. Z., Edgar, S. G., Gavin, J. B. & Hunter, P. J. (1995). Am. J. Physiol. Heart Circ. Physiol. 269, H571-H582.], 2005[LeGrice, I., Pope, A. & Smaill, B. (2005). Interstitial Fibrosis in Heart Failure, pp. 3-21. Springer.]; Gilbert et al., 2012[Gilbert, S. H., Benoist, D., Benson, A. P., White, E., Tanner, S. F., Holden, A. V., Dobrzynski, H., Bernus, O. & Radjenovic, A. (2012). Am. J. Physiol. Heart Circ. Physiol. 302, H287-H298.]).

We have also shown that the heterogeneous deformation occurring as a consequence of tissue dehydration also had an impact on the 3D organization of the myocardium. After long-term immersion of the heart in ethanol we observed a significant increase of the FAST as well as significant changes in the HA. Our results suggest that when performing quantitative studies of myocardial organization it is very important that all the samples are processed in the same way, especially when comparing different groups, and that the possible effects of the use of chemical agents such as the non-uniform deformation of the tissue have to be taken into the account.

4.1. Limitations

The main limitation of the present study is the number of samples, since we have only performed the experiment in a single rat heart. However, we think that one sample was enough as a proof of principle to illustrate the effects of ethanol on myocardial tissue. Moreover, in this study we have also evaluated one protocol for tissue dehydration. Therefore, in order to optimize sample preparation, several ethanol protocols should be evaluated and compared, as well as different tissue types. In this way, study protocols can be optimized to tissue type and size, by maximizing contrast while minimizing tissue alterations.

5. Conclusions

In conclusion, we have quantitatively demonstrated that the use of ethanol as contrast agent in synchrotron X-PCI enhances image contract but the tissue deformation occurring as consequence of dehydration is not uniform, thus altering the 3D myocardial organization.

6. Related literature

The following references, not cited in the main body of the paper, have been cited in the supporting information: Pennec et al (1999[Pennec, X., Cachier, P. & Ayache, N. (1999). Medical Image Computing and Computer-Assisted Intervention (MICCAI'99), Vol. 1679 of Lecture Notes in Computer Science, pp. 597-606. Berlin, Heidelberg: Springer.]); Thirion (1998[Thirion, J. (1998). Med. Image Anal. 2, 243-260.]).

Supporting information


Acknowledgements

We acknowledge the Paul Scherrer Institut, Villigen, Switzerland, for provision of synchrotron radiation beam time at the TOMCAT beamline X02DA of the SLS and would like to thank Anne Bonnin for assistance.

Funding information

The following funding is acknowledged: MICIU/AEI /10.13039/501100011033 and FSE+ (grant No. #RYC2022-035960-I to Gabriel Bernardino); Strategic Focal Area `Personalized Health and Related Technologies (PHRT)' (grant No. #2017-303 to Hector Dejea); Chan Zuckerberg Initiative DAF (grant No. 2022-316777 to Hector Dejea); Fundació Daniel Bravo i Andreu (grant to Patricia Garcia-Canadilla).

References

First citationAlbers, J., Markus, M. A., Alves, F. & Dullin, C. (2018). Sci. Rep. 8, 7712.  Web of Science CrossRef PubMed Google Scholar
First citationAnderson, R. H., Niederer, P. F., Sanchez–Quintana, D., Stephenson, R. S. & Agger, P. (2019). J. Anat. 235, 697–705.  CrossRef PubMed Google Scholar
First citationBaličević, V., Lončarić, S., Cárdenes, R., Gonzalez-Tendero, A., Paun, B., Crispi, F., Butakoff, C. & Bijnens, B. (2015). Functional Imaging and Modeling of the Heart, Vol. 9126 of Lecture Notes in Computer Science, pp. 111–119. Springer.  Google Scholar
First citationBeare, R., Lowekamp, B. & Yaniv, Z. (2018). J. Stat. Soft. 86, 1–35.  CrossRef Google Scholar
First citationBerg, S., Kutra, D., Kroeger, T., Straehle, C. N., Kausler, B. X., Haubold, C., Schiegg, M., Ales, J., Beier, T., Rudy, M., Eren, K., Cervantes, J. I., Xu, B., Beuttenmueller, F., Wolny, A., Zhang, C., Koethe, U., Hamprecht, F. A. & Kreshuk, A. (2019). Nat. Methods, 16, 1226–1232.  Web of Science CrossRef CAS PubMed Google Scholar
First citationBournonville, S. de, Vangrunderbeeck, S. & Kerckhofs, G. (2019). Contrast Media Mol. Imaging, 2019, 1–9.  Google Scholar
First citationBravin, A., Coan, P. & Suortti, P. (2013). Phys. Med. Biol. 58, R1–R35.  Web of Science CrossRef PubMed Google Scholar
First citationBrombal, L., Arfelli, F., Donato, S., Rigon, L., Golosio, B., Oliva, P., Longo, R., Bonazza, D., Zanconati, F., Contillo, A., Taibi, A., Delogu, P., Mettivier, G. & Tromba, G. (2019). J. Med. Imaging, 6, 031402.  Google Scholar
First citationBrombal, L., Donato, S., Dreossi, D., Arfelli, F., Bonazza, D., Contillo, A., Delogu, P., Di Trapani, V., Golosio, B., Mettivier, G., Oliva, P., Rigon, L., Taibi, A. & Longo, R. (2018). Phys. Med. Biol. 63, aaf2e1.  CrossRef Google Scholar
First citationBusse, M., Müller, M., Kimm, M. A., Ferstl, S., Allner, S., Achterhold, K., Herzen, J. & Pfeiffer, F. (2018). Proc. Natl Acad. Sci. USA, 115, 2293–2298.  Web of Science CrossRef CAS PubMed Google Scholar
First citationBuytaert, J., Goyens, J., De Greef, D., Aerts, P. & Dirckx, J. (2014). Microsc. Microanal. 20, 1208–1217.  CrossRef PubMed Google Scholar
First citationCibc (2013). Seg3D: Volumetric Image Segmentation and Visualization. Scientific Computing and Imaging Institute (SCI), https://www.mendeley.com/catalogue/a702d9d1-b36d-38f0-9c4b-44f7c552a1b2/?utm_source=desktop&utm_medium=1.17.9&utm_campaign=open_catalog&userDocumentId=%7Bc869e468-ef14-4cd4-988b-46570e2Google Scholar
First citationDejea, H., Garcia-Canadilla, P., Cook, A. C., Guasch, E., Zamora, M., Crispi, F., Stampanoni, M., Bijnens, B. & Bonnin, A. (2019). Sci. Rep. 9, 6996.  Web of Science CrossRef PubMed Google Scholar
First citationDudak, J., Zemlicka, J., Karch, J., Patzelt, M., Mrzilkova, J., Zach, P., Hermanova, Z., Kvacek, J. & Krejci, F. (2016). Sci. Rep. 6, 30385.  Web of Science CrossRef PubMed Google Scholar
First citationESC (2019). About Cardiovascular Disease in ESC Member Countries, https://www.escardio.org/The-ESC/Press-Office/Fact-sheetsGoogle Scholar
First citationEssen, H. F. van, Verdaasdonk, M. A. M., Elshof, S. M., de Weger, R. A. & van Diest, P. J. (2010). J. Clin. Pathol. 63, 1090–1094.  PubMed Google Scholar
First citationGarcia-Canadilla, P., Dejea, H., Bonnin, A., Balicevic, V., Loncaric, S., Zhang, C., Butakoff, C., Aguado-Sierra, J., Vázquez, M., Jackson, L. H., Stuckey, D. J., Rau, C., Stampanoni, M., Bijnens, B. & Cook, A. C. (2018). Circ. Cardiovasc. Imaging, 11, e007753.  PubMed Google Scholar
First citationGarcia-Canadilla, P., de Vries, T., Gonzalez-Tendero, A., Bonnin, A., Gratacos, E., Crispi, F., Bijnens, B. & Zhang, C. (2019). PLoS One, 14, e0218192.  PubMed Google Scholar
First citationGilbert, S. H., Benoist, D., Benson, A. P., White, E., Tanner, S. F., Holden, A. V., Dobrzynski, H., Bernus, O. & Radjenovic, A. (2012). Am. J. Physiol. Heart Circ. Physiol. 302, H287–H298.  CrossRef PubMed Google Scholar
First citationGonzalez-Tendero, A., Zhang, C., Balicevic, V., Cardenes, R., Loncaric, S., Butakoff, C., Paun, B., Bonnin, A., Garcia-Cañadilla, P., Muñoz-Moreno, E., Gratacos, E., Crispi, F. & Bijnens, B. (2017). Eur Heart J. Cardiovasc. Imaging, 18, 732–741.  PubMed Google Scholar
First citationHaque, Z., Rahman, A., Khan, Z. I., Hussan, M. T. & Alam, M. (2020). Int. J. Morphol. 38, 1371–1375.  CrossRef Google Scholar
First citationHołda, M. K., Klimek–Piotrowska, W., Koziej, M., Piątek, K. & Hołda, J. (2016). J. Anat. 229, 334–340.  PubMed Google Scholar
First citationHoshino, M., Uesugi, K., Tsukube, T. & Yagi, N. (2014). J. Synchrotron Rad. 21, 1347–1357.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationKoç, M. M., Aslan, N., Kao, A. P. & Barber, A. H. (2019). Microsc. Res. Tech. 82, 812–848.  PubMed Google Scholar
First citationKunii, T., Shirai, R., Yoneyama, A., Maruyama, H., Lwin, T. T. & Takeda, T. (2013). Med. Imaging Technol. 31, 132–135.  Google Scholar
First citationO, J., Kwon, H. J., Cho, T. H., Woo, S. H., Rhee, Y. H. & Yang, H. M. (2021). PLoS One, 16, e0254264.  PubMed Google Scholar
First citationLeGrice, I., Pope, A. & Smaill, B. (2005). Interstitial Fibrosis in Heart Failure, pp. 3–21. Springer.  Google Scholar
First citationLeGrice, I. J., Smaill, B. H., Chai, L. Z., Edgar, S. G., Gavin, J. B. & Hunter, P. J. (1995). Am. J. Physiol. Heart Circ. Physiol. 269, H571–H582.  CrossRef Google Scholar
First citationLoncaric, F., Garcia-Canadilla, P., Garcia-Alvarez, A., Sanchis, L., Prat, S., Doltra, A., Quintana, E., Pereda, D., Dejea, H., Bonnin, A., Sitges, M. & Bijnens, B. (2021). Front. Cardiovasc. Med. 8, 670734.  CrossRef PubMed Google Scholar
First citationLowekamp, B. C., Chen, D. T., Ibáñez, L. & Blezek, D. (2013). Front. Neuroinform. 7, 73735.  CrossRef Google Scholar
First citationLwin, T.-T., Yoneyama, A., Kokubo, S., Maruyama, H., Hyodo, K. & Takeda, T. (2022). Acta Radiol. 63, 1102–1109.  CrossRef PubMed Google Scholar
First citationMarone, F. & Stampanoni, M. (2012). J. Synchrotron Rad. 19, 1029–1037.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationMattes, D., Haynor, D. R., Vesselle, H., Lewellyn, T. K. & Eubank, W. (2001). Proc. SPIE, 4322, 1609–1620.  CrossRef Google Scholar
First citationMetscher, B. D. (2009). BMC Physiol. 9, 11.  Google Scholar
First citationMirea, I., Varray, F., Zhu, Y. M., Fanton, L., Langer, M., Jouk, P. S., Michalowicz, G., Usson, Y. & Magnin, I. E. (2015). Functional Imaging and Modeling of the Heart, Vol. 9126 of Lecture Notes in Computer Science, pp. 172–179. Springer.  Google Scholar
First citationMrzílková, J., Patzelt, M., Gallina, P., Wurst, Z., Šeremeta, M., Dudák, J., Krejčí, F., Žemlička, J., Musil, V., Karch, J., Rosina, J. & Zach, P. (2019). Biomed Res. Intl, 2019, 1–7.  Google Scholar
First citationNelson, L. J., Smith, R. A. & Mienczakowski, M. (2018). Composites Part A, 104, 108–119.  CrossRef Google Scholar
First citationPatzelt, M., Mrzilkova, J., Dudak, J., Krejci, F., Zemlicka, J., Karch, J., Musil, V., Rosina, J., Sykora, V., Horehledova, B. & Zach, P. (2019). Jpn J. Radiol. 37, 500–510.  CrossRef PubMed Google Scholar
First citationPennec, X., Cachier, P. & Ayache, N. (1999). Medical Image Computing and Computer-Assisted Intervention (MICCAI'99), Vol. 1679 of Lecture Notes in Computer Science, pp. 597–606. Berlin, Heidelberg: Springer.  Google Scholar
First citationPichat, J., Iglesias, J. E., Yousry, T., Ourselin, S. & Modat, M. (2018). Med. Image Anal. 46, 73–105.  CrossRef PubMed Google Scholar
First citationPlaninc, I., Garcia-Canadilla, P., Dejea, H., Ilic, I., Guasch, E., Zamora, M., Crispi, F., Stampanoni, M., Milicic, D., Bijnens, B., Bonnin, A. & Cikes, M. (2021). Sci. Rep. 11, 14020.  CrossRef PubMed Google Scholar
First citationPlaninc, I., Ilic, I., Dejea, H., Garcia-Canadilla, P., Gasparovic, H., Jurin, H., Milicic, D., Skoric, B., Stampanoni, M., Bijnens, B., Bonnin, A. & Cikes, M. (2023). Transpl. Int. 36, 11046.  CrossRef PubMed Google Scholar
First citationReichardt, M., Frohn, J., Khan, A., Alves, F. & Salditt, T. (2020b). Biomed. Opt. Expr. 11, 2633.  CrossRef Google Scholar
First citationReichardt, M., Töpperwien, M., Khan, A., Alves, F. & Salditt, T. (2020a). J. Med. Imag. 7, 1.  CrossRef Google Scholar
First citationRogelj, P. & Kovačič, S. (2006). Med. Image Anal. 10, 484–493.  CrossRef PubMed Google Scholar
First citationSchulz, G., Crooijmans, H. J. A., Germann, M., Scheffler, K., Müller-Gerbl, M. & Müller, B. (2011). J. Neurosci. Methods, 202, 17–27.  Web of Science CrossRef PubMed Google Scholar
First citationShirai, R., Kunii, T., Yoneyama, A., Ooizumi, T., Maruyama, H., Lwin, T.-T., Hyodo, K. & Takeda, T. (2014). J. Synchrotron Rad. 21, 795–800.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationSilva, J. M. S., Zanette, I., Noël, P. B., Cardoso, M. B., Kimm, M. A. & Pfeiffer, F. (2015). Sci. Rep. 5, 14088.  CrossRef PubMed Google Scholar
First citationSommer, C., Sträehle, C., Köthe, U. & Hamprecht, F. A. (2011). 2011 IEEE International Symposium on Biomedical Imaging, pp. 230–233. IEEE.  Google Scholar
First citationSoveral, I., Crispi, F., Walter, C., Guirado, L., García–Cañadilla, P., Cook, A., Bonnin, A., Dejea, H., Rovira–Zurriaga, C., Sánchez de Toledo, J., Gratacós, E., Martínez, J. M., Bijnens, B. & Gómez, O. (2020). Ultrasound Obstet. Gyne, 56, 837–849.  CrossRef Google Scholar
First citationTakeda, T., Kunii, T., Sirai, R., Ohizumi, T., Maruyama, H., Hyodo, K., Yoneyama, A. & Ueda, K. (2013). J. Phys. Conf. Ser. 425, 022004.  CrossRef Google Scholar
First citationTakeda, T., Lwin, T. T., Yoneyama, A. & Hyodo, K. (2013). Photon Factory Activity Report 2013. Photon Factory, Tsukuba, Japan.  Google Scholar
First citationTakeda, T., Yoneyama, A. & Hyodo, K. (2014). Photon Factory Activity Report 2014. Photon Factory, Tsukuba, Japan.  Google Scholar
First citationThirion, J. (1998). Med. Image Anal. 2, 243–260.  CrossRef PubMed Google Scholar
First citationVarray, F., Mirea, I., Langer, M., Peyrin, F., Fanton, L. & Magnin, I. E. (2017). Med. Image Anal. 38, 117–132.  CrossRef PubMed Google Scholar
First citationVarray, F., Wang, L., Fanton, L., Zhu, Y.-M. & Magnin, I. E. (2013). Functional Imaging and Modeling of the Heart, Vol. 7945 of Lecture Notes in Computer Science, pp. 150–157. Berlin, Heidelberg: Springer.  Google Scholar
First citationVickerton, P., Jarvis, J. & Jeffery, N. (2013). J. Anat. 223, 185–193.  Web of Science CrossRef PubMed Google Scholar
First citationWang, H., Dong, L., O'Daniel, J., Mohan, R., Garden, A. S., Ang, K. K., Kuban, D. A., Bonnen, M., Chang, J. Y. & Cheung, R. (2005). Phys. Med. Biol. 50, 2887–2905.  CrossRef PubMed Google Scholar
First citationWang, S., Mirea, I., Varray, F., Liu, W. Y. & Magnin, I. E. (2019). Lecture Notes in Computer Science (Including Subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics). Springer.  Google Scholar
First citationYaniv, Z., Lowekamp, B. C., Johnson, H. J. & Beare, R. (2018). J. Digit. Imaging, 31, 290–303.  CrossRef PubMed Google Scholar

This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775
Follow J. Synchrotron Rad.
Sign up for e-alerts
Follow J. Synchrotron Rad. on Twitter
Follow us on facebook
Sign up for RSS feeds