Serial block face scanning electron microscopy for the study of cardiac muscle ultrastructure at nanoscale resolutions

Serial Block Face-Scanning Electron Microscopy as a Burgeoning Technology

Serial block face scanning electron microscopy (SBF-SEM), also referred to as serial block-face electron microscopy, is an advanced ultrastructural imaging technique that enables three-dimensional visualization that provides largerx- and y-axis ranges than other volumetric EM techniques. While SEM is first introduced in the 1930s, SBF-SEM is developed as a novel method to resolve the 3D architecture of neuronal networks across large volumes with nanometer resolution by Denk and Horstmann in 2004. Here, the authors provide an accessible overview of the advantages and challenges associated with SBF-SEM. Beyond this, the applications of SBF-SEM in biochemical domains as well as potential future clinical applications are briefly reviewed. Finally, the alternative forms of artificial intelligence-based segmentation which may contribute to devising a feasible workflow involving SBF-SEM, are also considered.

Quantitative super-resolution microscopy reveals the differences in the nanoscale distribution of nuclear phosphatidylinositol 4,5-bisphosphate in human healthy skin and skin warts

Introduction: Imaging of human clinical formalin-fixed paraffin-embedded (FFPE) tissue sections provides insights into healthy and diseased states and therefore represents a valuable resource for basic research, as well as for diagnostic and clinical purposes. However, conventional light microscopy does not allow to observe the molecular details of tissue and cell architecture due to the diffraction limit of light. Super-resolution microscopy overcomes this limitation and provides access to the nanoscale details of tissue and cell organization. Methods: Here, we used quantitative multicolor stimulated emission depletion (STED) nanoscopy to study the nanoscale distribution of the nuclear phosphatidylinositol 4,5-bisphosphate (nPI(4,5)P2) with respect to the nuclear speckles (NS) marker SON. Results: Increased nPI(4,5)P2 signals were previously linked to human papillomavirus (HPV)-mediated carcinogenesis, while NS-associated PI(4,5)P2 represents the largest pool of nPI(4,5)P2 visualized by staining and microscopy. The implementation of multicolor STED nanoscopy in human clinical FFPE skin and wart sections allowed us to provide here the quantitative evidence for higher levels of NS-associated PI(4,5)P2 in HPV-induced warts compared to control skin. Discussion: These data expand the previous reports of HPV-induced increase of nPI(4,5)P2 levels and reveal for the first time the functional, tissue-specific localization of nPI(4,5)P2 within NS in clinically relevant samples. Moreover, our approach is widely applicable to other human clinical FFPE tissues as an informative addition to the classical histochemistry.

Tissue Clearing and Its Application in the Musculoskeletal System

The musculoskeletal system is an integral part of the human body. Currently, most skeletal muscle research is conducted through conventional histological sections due to technological limitations and the structure of skeletal muscles. For studying and observing bones and muscles, there is an urgent need for three-dimensional, objective imaging technologies. Optical tissue-clearing technologies seem to offer a novel and accessible approach to research of the musculoskeletal system. Using this approach, the components which cause refraction or prevent light from penetrating into the tissue are physically and chemically eliminated; then the liquid in the tissue is replaced with high-refractive-index chemicals. This innovative method, which allows three-dimensional reconstruction at the cellular and subcellular scale, significantly improves imaging depth and resolution. Nonetheless, this technology was not originally developed to image bones or muscles. When compared with brain and nerve organs which have attracted considerable attention in this field, the musculoskeletal system contains fewer lipids and has high levels of hemoglobin, collagen fibers, and inorganic hydroxyapatite crystals. Currently, three-dimensional imaging methods are widely used in the diagnosis and treatment of skeletal and muscular illnesses. In this regard, it is vitally important to review and evaluate the optical tissue-clearing technologies currently employed in the musculoskeletal system, so that researchers may make an informed decision. In the meantime, this study offers guidelines and recommendations for expanding the use of this technology in the musculoskeletal system.

A cell-based framework for modeling cardiac mechanics

Cardiomyocytes are the functional building blocks of the heart-yet most models developed to simulate cardiac mechanics do not represent the individual cells and their surrounding matrix. Instead, they work on a homogenized tissue level, assuming that cellular and subcellular structures and processes scale uniformly. Here we present a mathematical and numerical framework for exploring tissue-level cardiac mechanics on a microscale given an explicit three-dimensional geometrical representation of cells embedded in a matrix. We defined a mathematical model over such a geometry and parametrized our model using publicly available data from tissue stretching and shearing experiments. We then used the model to explore mechanical differences between the extracellular and the intracellular space. Through sensitivity analysis, we found the stiffness in the extracellular matrix to be most important for the intracellular stress values under contraction. Strain and stress values were observed to follow a normal-tangential pattern concentrated along the membrane, with substantial spatial variations both under contraction and stretching. We also examined how it scales to larger size simulations, considering multicellular domains. Our work extends existing continuum models, providing a new geometrical-based framework for exploring complex cell-cell and cell-matrix interactions.

Three-dimensional visualization of the cardiac ryanodine receptor clusters and the molecular-scale fraying of dyads

Clusters of ryanodine receptor calcium channels (RyRs) form the primary molecular machinery of intracellular calcium signalling in cardiomyocytes. While a range of optical super-resolution microscopy techniques have revealed the nanoscale structure of these clusters, the three-dimensional (3D) nanoscale topologies of the clusters have remained mostly unresolved. In this paper, we demonstrate the exploitation of molecular-scale resolution in enhanced expansion microscopy (EExM) along with various 2D and 3D visualization strategies to observe the topological complexities, geometries and molecular sub-domains within the RyR clusters. Notably, we observed sub-domains containing RyR-binding protein junctophilin-2 (JPH2) occupying the central regions of RyR clusters in the deeper interior of the myocytes (including dyads), while the poles were typically devoid of JPH2, lending to a looser RyR arrangement. By contrast, peripheral RyR clusters exhibited variable co-clustering patterns and ratios between RyR and JPH2. EExM images of dyadic RyR clusters in right ventricular (RV) myocytes isolated from rats with monocrotaline-induced RV failure revealed hallmarks of RyR cluster fragmentation accompanied by breaches in the JPH2 sub-domains. Frayed RyR patterns observed adjacent to these constitute new evidence that the destabilization of the RyR arrays inside the JPH2 sub-domains may seed the primordial foci of dyad remodelling observed in heart failure. This article is part of the theme issue 'The cardiomyocyte: new revelations on the interplay between architecture and function in growth, health, and disease'.

A joint alignment and reconstruction algorithm for electron tomography to visualize in-depth cell-to-cell interactions

Electron tomography allows one to obtain 3D reconstructions visualizing a tissue's ultrastructure from a series of 2D projection images. An inherent problem with this imaging technique is that its projection images contain unwanted shifts, which must be corrected for to achieve reliable reconstructions. Commonly, the projection images are aligned with each other by means of fiducial markers prior to the reconstruction procedure. In this work, we propose a joint alignment and reconstruction algorithm that iteratively solves for both the unknown reconstruction and the unintentional shift and does not require any fiducial markers. We evaluate the approach first on synthetic phantom data where the focus is not only on the reconstruction quality but more importantly on the shift correction. Subsequently, we apply the algorithm to healthy C57BL/6J mice and then compare it with non-obese diabetic (NOD) mice, with the aim of visualizing the attack of immune cells on pancreatic beta cells within type 1 diabetic mice at a more profound level through 3D analysis. We empirically demonstrate that the proposed algorithm is able to compute the shift with a remaining error at only the sub-pixel level and yields high-quality reconstructions for the limited-angle inverse problem. By decreasing labour and material costs, the algorithm facilitates further research directed towards investigating the immune system's attacks in pancreata of NOD mice for numerous samples at different stages of type 1 diabetes.

Electron microscopy of cardiac 3D nanodynamics: form, function, future

The 3D nanostructure of the heart, its dynamic deformation during cycles of contraction and relaxation, and the effects of this deformation on cell function remain largely uncharted territory. Over the past decade, the first inroads have been made towards 3D reconstruction of heart cells, with a native resolution of around 1 nm³, and of individual molecules relevant to heart function at a near-atomic scale. These advances have provided access to a new generation of data and have driven the development of increasingly smart, artificial intelligence-based, deep-learning image-analysis algorithms. By high-pressure freezing of cardiomyocytes with millisecond accuracy after initiation of an action potential, pseudodynamic snapshots of contraction-induced deformation of intracellular organelles can now be captured. In combination with functional studies, such as fluorescence imaging, exciting insights into cardiac autoregulatory processes at nano-to-micro scales are starting to emerge. In this Review, we discuss the progress in this fascinating new field to highlight the fundamental scientific insight that has emerged, based on technological breakthroughs in biological sample preparation, 3D imaging and data analysis; to illustrate the potential clinical relevance of understanding 3D cardiac nanodynamics; and to predict further progress that we can reasonably expect to see over the next 10 years.

Image-Driven Modeling of Nanoscopic Cardiac Function: Where Have We Come From, and Where Are We Going?

Complementary developments in microscopy and mathematical modeling have been critical to our understanding of cardiac excitation–contraction coupling. Historically, limitations imposed by the spatial or temporal resolution of imaging methods have been addressed through careful mathematical interrogation. Similarly, limitations imposed by computational power have been addressed by imaging macroscopic function in large subcellular domains or in whole myocytes. As both imaging resolution and computational tractability have improved, the two approaches have nearly merged in terms of the scales that they can each be used to interrogate. With this review we will provide an overview of these advances and their contribution to understanding ventricular myocyte function, including exciting developments over the last decade. We specifically focus on experimental methods that have pushed back limits of either spatial or temporal resolution of nanoscale imaging (e.g., DNA-PAINT), or have permitted high resolution imaging on large cellular volumes (e.g., serial scanning electron microscopy). We also review the progression of computational approaches used to integrate and interrogate these new experimental data sources, and comment on near-term advances that may unify understanding of the underlying biology. Finally, we comment on several outstanding questions in cardiac physiology that stand to benefit from a concerted and complementary application of these new experimental and computational methods.

Comprehensive assessment of myocardial remodeling in ischemic heart disease by synchrotron propagation based X-ray phase contrast imaging

Cardiovascular research is in an ongoing quest for a superior imaging method to integrate gross-anatomical information with microanatomy, combined with quantifiable parameters of cardiac structure. In recent years, synchrotron radiation-based X-ray Phase Contrast Imaging (X-PCI) has been extensively used to characterize soft tissue in detail. The objective was to use X-PCI to comprehensively quantify ischemic remodeling of different myocardial structures, from cell to organ level, in a rat model of myocardial infarction. Myocardial infarction-induced remodeling was recreated in a well-established rodent model. Ex vivo rodent hearts were imaged by propagation based X-PCI using two configurations resulting in 5.8 µm and 0.65 µm effective pixel size images. The acquired datasets were used for a comprehensive assessment of macrostructural changes including the whole heart and vascular tree morphology, and quantification of left ventricular myocardial thickness, mass, volume, and organization. On the meso-scale, tissue characteristics were explored and compared with histopathological methods, while microstructural changes were quantified by segmentation of cardiomyocytes and calculation of cross-sectional areas. Propagation based X-PCI provides detailed visualization and quantification of morphological changes on whole organ, tissue, vascular as well as individual cellular level of the ex vivo heart, with a single, non-destructive 3D imaging modality.

Nano-scale morphology of cardiomyocyte t-tubule/sarcoplasmic reticulum junctions revealed by ultra-rapid high-pressure freezing and electron tomography

Detailed knowledge of the ultrastructure of intracellular compartments is a prerequisite for our understanding of how cells function. In cardiac muscle cells, close apposition of transverse (t)-tubule (TT) and sarcoplasmic reticulum (SR) membranes supports stable high-gain excitation-contraction coupling. Here, the fine structure of this key intracellular element is examined in rabbit and mouse ventricular cardiomyocytes, using ultra-rapid high-pressure freezing (HPF, omitting aldehyde fixation) and electron microscopy. 3D electron tomograms were used to quantify the dimensions of TT, terminal cisternae of the SR, and the space between SR and TT membranes (dyadic cleft). In comparison to conventional aldehyde-based chemical sample fixation, HPF-preserved samples of both species show considerably more voluminous SR terminal cisternae, both in absolute dimensions and in terms of junctional SR to TT volume ratio. In rabbit cardiomyocytes, the average dyadic cleft surface area of HPF and chemically fixed myocytes did not differ, but cleft volume was significantly smaller in HPF samples than in conventionally fixed tissue; in murine cardiomyocytes, the dyadic cleft surface area was higher in HPF samples with no difference in cleft volume. In both species, the apposition of the TT and SR membranes in the dyad was more likely to be closer than 10 nm in HPF samples compared to CFD, presumably resulting from avoidance of sample shrinkage associated with conventional fixation techniques. Overall, we provide a note of caution regarding quantitative interpretation of chemically-fixed ultrastructures, and offer novel insight into cardiac TT and SR ultrastructure with relevance for our understanding of cardiac physiology.

A multifaceted analysis reveals two distinct phases of chloroplast biogenesis during de-etiolation in Arabidopsis

Light triggers chloroplast differentiation whereby the etioplast transforms into a photosynthesizing chloroplast and the thylakoid rapidly emerges. However, the sequence of events during chloroplast differentiation remains poorly understood. Using Serial Block Face Scanning Electron Microscopy (SBF-SEM), we generated a series of chloroplast 3D reconstructions during differentiation, revealing chloroplast number and volume and the extent of envelope and thylakoid membrane surfaces. Furthermore, we used quantitative lipid and whole proteome data to complement the (ultra)structural data, providing a time-resolved, multi-dimensional description of chloroplast differentiation. This showed two distinct phases of chloroplast biogenesis: an initial photosynthesis-enabling 'Structure Establishment Phase' followed by a 'Chloroplast Proliferation Phase' during cell expansion. Moreover, these data detail thylakoid membrane expansion during de-etiolation at the seedling level and the relative contribution and differential regulation of proteins and lipids at each developmental stage. Altogether, we establish a roadmap for chloroplast differentiation, a critical process for plant photoautotrophic growth and survival.

Multiscale cardiac imaging spanning the whole heart and its internal cellular architecture in a small animal model

Cardiac pumping depends on the morphological structure of the heart, but also on its subcellular (ultrastructural) architecture, which enables cardiac contraction. In cases of congenital heart defects, localized ultrastructural disruptions that increase the risk of heart failure are only starting to be discovered. This is in part due to a lack of technologies that can image the three-dimensional (3D) heart structure, to assess malformations; and its ultrastructure, to assess organelle disruptions. We present here a multiscale, correlative imaging procedure that achieves high-resolution images of the whole heart, using 3D micro-computed tomography (micro-CT); and its ultrastructure, using 3D scanning electron microscopy (SEM). In a small animal model (chicken embryo), we achieved uniform fixation and staining of the whole heart, without losing ultrastructural preservation on the same sample, enabling correlative multiscale imaging. Our approach enables multiscale studies in models of congenital heart disease and beyond.

Postnatal cardiovascular adaptation

The heart undergoes rapid transformations in function during the transition to extrauterine life. Our understanding of the adaptive physiology underlying this process is able to inform the clinical management of infants who are struggling to complete this complex transition. Much of our knowledge of the cardiac transition is derived from the preterm infant in whom the preparative adaptations are incomplete and clinical sequelae all too common. This review will re-examine the cardiac transition highlighting the physiology that drives it and suggest appropriate clinical intervention to support the process.

Diffusional and Electrical Properties of T-Tubules Are Governed by Their Constrictions and Dilations

Cardiac t-tubules (TTs) form a network of complex surface membrane invaginations that is essential for proper excitation-contraction coupling. Although electron and optical microscopy studies provided a wealth of important information about the structure of TTs, assessing their functional properties remains a challenge. In this study, we investigated the diffusional accessibility of TTs in intact isolated adult mouse ventricular myocytes using, to our knowledge, a novel fluorescence-based assay. In this approach, a small part of TTs is first locally filled with fluorescent dextran and then its diffusion out of TTs is monitored after rapid removal of extracellular dextran. In normal cells, diffusion of 3 kDa dextran is characterized by an average time constant of 3.9 ± 1.2 s with the data ranging from 1.8 to 10.5 s. The data are consistent with essentially free diffusion of dextran in TTs although measurable contribution of binding is also evident. TT fluorescence is abolished in cells treated with high concentration of formamide or after hyposmotic stress. Importantly, the assay we use allows for quantitative, repetitive measurements of subtle dynamic changes in TT structure of the same cell that are not possible to observe with other approaches. In particular, dextran diffusion rate decreases two-to-threefold during cell swelling, suggesting significant structural remodeling of TTs. Computer modeling shows that diffusional accessibility and electrical properties of TTs are primarily determined by the constrictions and dilations of individual TTs and that, from a functional perspective, TTs cannot be considered as a network of cylinders of the same average diameter. Constriction/dilation model of cardiac TTs is in a quantitative agreement with previous high-resolution microscopy studies of TT structure and alternative measurements of diffusional and electrical time constants of TTs. The data also show that the apparent electrical length constant of cardiac TTs is likely several-fold smaller than that estimated in earlier studies.

Hypertension-induced changes in the rat myocardium during the development of cardiac hypertrophy - a comparison between the left and the right ventricle

The hypertrophy of the cardiac muscle is one of the most significant maladaptive mechanisms activated in response to increased workload. It is associated with histological and ultrastructural alterations, changes in the quantitative parameters and the expression of different enzymes. While the structural and functional consequences of systemic hypertension on the left ventricle have been well evaluated, the right ventricle has received less attention. The aim of the present study was to analyse and compare the changes in the left and right ventricle during the development of cardiac hypertrophy initiated by systemic hypertension in different age groups of spontaneously hypertensive rats. Therefore, we studied the histology and ultrastructure of the cells of the myocardium, evaluated the immunohistochemical expression of the enzyme neuronal nitric oxide synthase and conducted a quantitative analysis of several morphometric parameters. We used three groups of spontaneously hypertensive rats. For the quantitative analysis we also used three age groups of age- and weight-matched control animals (normotensive Wistar rats). In both ventricles, we described cardiomyocytic hypertrophy, focal myocytolysis and increased collagen deposition in the interstitial space. Our observations on the ultrastructural level were associated with changes in the cardiomyocytic nuclei, the arrangement, maturity and organisation of the myofibrils, the localisation and ultrastructure of the mitochondria, the development and maturity of the intercalated discs, as well as changes in the components of the interstitium. The immunohistochemical expression of neuronal nitric oxide synthase in the left ventricle was stronger than that in the right ventricle across all age groups. The comparative quantitative analysis revealed that changes in the studied morphometric parameters in the two ventricles occurred disproportionately. In conclusion, the present study characterised the development of cardiac hypertrophy in response to systemic hypertension in both ventricles and demonstrated the involvement of the right ventricle.

Subcellular connectomic analyses of energy networks in striated muscle

Mapping biological circuit connectivity has revolutionized our understanding of structure-function relationships. Although connectomic analyses have primarily focused on neural systems, electrical connectivity within muscle mitochondrial networks was recently demonstrated to provide a rapid mechanism for cellular energy distribution. However, tools to evaluate organelle connectivity with high spatial fidelity within single cells are currently lacking. Here, we developed a framework to quantitatively assess mitochondrial network connectivity and interactions with cellular sites of energy storage, utilization, and calcium cycling in cardiac, oxidative, and glycolytic muscle. We demonstrate that mitochondrial network configuration, individual mitochondrial size and shape, and the junctions connecting mitochondria within each network are consistent with the differing contraction demands of each muscle type. Moreover, mitochondria-lipid droplet interaction analyses suggest that individual mitochondria within networks may play specialized roles regarding energy distribution and calcium cycling within the cell and reveal the power of connectomic analyses of organelle interactions within single cells.

Three-dimensional electron microscopy techniques for unravelling mitochondrial dysfunction in heart failure and identification of new pharmacological targets

A hallmark of heart failure is mitochondrial dysfunction leading to a bioenergetics imbalance in the myocardium. Consequently, there is much interest in targeting mitochondrial abnormalities to attenuate the pathogenesis of heart failure. This review discusses (i) how electron microscopy (EM) techniques have been fundamental for the current understanding of mitochondrial structure–function, (ii) the paradigm shift in resolutions now achievable by 3‐D EM techniques due to the introduction of direct detection devices and phase plate technology, and (iii) the application of EM for unravelling mitochondrial pathological remodelling in heart failure. We further consider the tremendous potential of multi‐scale EM techniques for the development of therapeutics, structure‐based ligand design and for delineating how a drug elicits nanostructural effects at the molecular, organelle and cellular levels. In conclusion, 3‐D EM techniques have entered a new era of structural biology and are poised to play a pivotal role in discovering new therapies targeting mitochondria for treating heart failure.

Linked Articles

This article is part of a themed section on Mitochondrial Pharmacology: Featured Mechanisms and Approaches for Therapy Translation. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v176.22/issuetoc

Nano-structural analysis of engrafted human induced pluripotent stem cell-derived cardiomyocytes in mouse hearts using a genetic-probe APEX2

Many studies have shown the feasibility of in vivo cardiac transplantation of human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) in animal experiments. However, nano-structural confirmation of the successful incorporation of the engrafted iPSC-CMs including electron microscopy (EM) has not been accomplished, partly because identification of graft cells in EM has proven to be difficult. Using APEX2, an engineered ascorbate peroxidase imaging tag, we successfully localized and analyzed the fine structure of sarcomeres and the excitation contraction machinery of iPSC-CMs 6 months after their engraftment in infarcted mouse hearts. APEX2 made iPSC-CMs visible in multiple imaging modalities including light microscopy, X-ray microscopic tomography, transmission EM, and scanning EM. EM tomography allowed assessment of the differentiation state of APEX2-positive iPSC-CMs and analysis of the fine structure of the sarcomeres including T-tubules and dyads.

First evidence of octacalcium phosphate@osteocalcin nanocomplex as skeletal bone component directing collagen triple-helix nanofibril mineralization

Tibia trabeculae and vertebrae of rats as well as human femur were investigated by high-resolution TEM at the atomic scale in order to reveal snapshots of the morphogenetic processes of local bone ultrastructure formation. By taking into account reflections of hydroxyapatite for Fourier filtering the appearance of individual alpha-chains within the triple-helix clearly shows that bone bears the feature of an intergrowth composite structure extending from the atomic to the nanoscale, thus representing a molecular composite of collagen and apatite. Careful Fourier analysis reveals that the non-collagenous protein osteocalcin is present directly combined with octacalcium phosphate. Besides single spherical specimen of about 2 nm in diameter, osteocalcin is spread between and over collagen fibrils and is often observed as pearl necklace strings. In high-resolution TEM, the three binding sites of the γ-carboxylated glutamic acid groups of the mineralized osteocalcin were successfully imaged, which provide the chemical binding to octacalcium phosphate. Osteocalcin is attached to the collagen structure and interacts with the Ca-sites on the (100) dominated hydroxyapatite platelets with Ca-Ca distances of about 9.5 Å. Thus, osteocalcin takes on the functions of Ca-ion transport and suppression of hydroxyapatite expansion.

A guide to analysis and reconstruction of serial block face scanning electron microscopy data

Serial block face scanning electron microscopy (SBF-SEM) is a relatively new technique that allows the acquisition of serially sectioned, imaged and digitally aligned ultrastructural data. There is a wealth of information that can be obtained from the resulting image stacks but this presents a new challenge for researchers - how to computationally analyse and make best use of the large datasets produced. One approach is to reconstruct structures and features of interest in 3D. However, the software programmes can appear overwhelming, time-consuming and not intuitive for those new to image analysis. There are a limited number of published articles that provide sufficient detail on how to do this type of reconstruction. Therefore, the aim of this paper is to provide a detailed step-by-step protocol, accompanied by tutorial videos, for several types of analysis programmes that can be used on raw SBF-SEM data, although there are more options available than can be covered here. To showcase the programmes, datasets of skeletal muscle from foetal and adult guinea pigs are initially used with procedures subsequently applied to guinea pig cardiac tissue and locust brain. The tissue is processed using the heavy metal protocol developed specifically for SBF-SEM. Trimmed resin blocks are placed into a Zeiss Sigma SEM incorporating the Gatan 3View and the resulting image stacks are analysed in three different programmes, Fiji, Amira and MIB, using a range of tools available for segmentation. The results from the image analysis comparison show that the analysis tools are often more suited to a particular type of structure. For example, larger structures, such as nuclei and cells, can be segmented using interpolation, which speeds up analysis; single contrast structures, such as the nucleolus, can be segmented using the contrast-based thresholding tools. Knowing the nature of the tissue and its specific structures (complexity, contrast, if there are distinct membranes, size) will help to determine the best method for reconstruction and thus maximize informative output from valuable tissue.

Creating a Structurally Realistic Finite Element Geometric Model of a Cardiomyocyte to Study the Role of Cellular Architecture in Cardiomyocyte Systems Biology

With the advent of three-dimensional (3D) imaging technologies such as electron tomography, serial-block-face scanning electron microscopy and confocal microscopy, the scientific community has unprecedented access to large datasets at sub-micrometer resolution that characterize the architectural remodeling that accompanies changes in cardiomyocyte function in health and disease. However, these datasets have been under-utilized for investigating the role of cellular architecture remodeling in cardiomyocyte function. The purpose of this protocol is to outline how to create an accurate finite element model of a cardiomyocyte using high resolution electron microscopy and confocal microscopy images. A detailed and accurate model of cellular architecture has significant potential to provide new insights into cardiomyocyte biology, more than experiments alone can garner. The power of this method lies in its ability to computationally fuse information from two disparate imaging modalities of cardiomyocyte ultrastructure to develop one unified and detailed model of the cardiomyocyte. This protocol outlines steps to integrate electron tomography and confocal microscopy images of adult male Wistar (name for a specific breed of albino rat) rat cardiomyocytes to develop a half-sarcomere finite element model of the cardiomyocyte. The procedure generates a 3D finite element model that contains an accurate, high-resolution depiction (on the order of ~35 nm) of the distribution of mitochondria, myofibrils and ryanodine receptor clusters that release the necessary calcium for cardiomyocyte contraction from the sarcoplasmic reticular network (SR) into the myofibril and cytosolic compartment. The model generated here as an illustration does not incorporate details of the transverse-tubule architecture or the sarcoplasmic reticular network and is therefore a minimal model of the cardiomyocyte. Nevertheless, the model can already be applied in simulation-based investigations into the role of cell structure in calcium signaling and mitochondrial bioenergetics, which is illustrated and discussed using two case studies that are presented following the detailed protocol.

An automated workflow for segmenting single adult cardiac cells from large-volume serial block-face scanning electron microscopy data

This paper presents a new algorithm to automatically segment the myofibrils, mitochondria and nuclei within single adult cardiac cells that are part of a large serial-block-face scanning electron microscopy (SBF-SEM) dataset. The algorithm only requires a set of manually drawn contours that roughly demarcate the cell boundary at routine slice intervals (every 50th, for example). The algorithm correctly classified pixels within the single cell with 97% accuracy when compared to manual segmentations. One entire cell and the partial volumes of two cells were segmented. Analysis of segmentations within these cells showed that myofibrils and mitochondria occupied 47.5% and 51.6% on average respectively, while the nuclei occupy 0.7% of the cell for which the entire volume was captured in the SBF-SEM dataset. Mitochondria clustering increased at the periphery of the nucleus region and branching points of the cardiac cell. The segmentations also showed high area fraction of mitochondria (up to 70% of the 2D image slice) in the sub-sarcolemmal region, whilst it was closer to 50% in the intermyofibrillar space. We finally demonstrate that our segmentations can be turned into 3D finite element meshes for cardiac cell computational physiology studies. We offer our large dataset and MATLAB implementation of the algorithm for research use at www.github.com/CellSMB/sbfsem-cardiac-cell-segmenter/. We anticipate that this timely tool will be of use to cardiac computational and experimental physiologists alike who study cardiac ultrastructure and its role in heart function.

A new look at the heart-novel imaging techniques

The development and successful implementation of cutting-edge imaging technologies to visualise cardiac anatomy and function is a key component of effective diagnostic efforts in cardiology. Here, we describe a number of recent exciting advances in the field of cardiology spanning from macro- to micro- to nano-scales of observation, including magnetic resonance imaging, computed tomography, optical mapping, photoacoustic imaging, and electron tomography. The methodologies discussed are currently making the transition from scientific research to routine clinical use, albeit at different paces. We discuss the most likely trajectory of this transition into clinical research and standard diagnostics, and highlight the key challenges and opportunities associated with each of the methodologies.

A computational model of spatio-temporal cardiac intracellular calcium handling with realistic structure and spatial flux distribution from sarcoplasmic reticulum and t-tubule reconstructions

Intracellular calcium cycling is a vital component of cardiac excitation-contraction coupling. The key structures responsible for controlling calcium dynamics are the cell membrane (comprising the surface sarcolemma and transverse-tubules), the intracellular calcium store (the sarcoplasmic reticulum), and the co-localisation of these two structures to form dyads within which calcium-induced-calcium-release occurs. The organisation of these structures tightly controls intracellular calcium dynamics. In this study, we present a computational model of intracellular calcium cycling in three-dimensions (3-D), which incorporates high resolution reconstructions of these key regulatory structures, attained through imaging of tissue taken from the sheep left ventricle using serial block face scanning electron microscopy. An approach was developed to model the sarcoplasmic reticulum structure at the whole-cell scale, by reducing its full 3-D structure to a 3-D network of one-dimensional strands. The model reproduces intracellular calcium dynamics during control pacing and reveals the high-resolution 3-D spatial structure of calcium gradients and intracellular fluxes in both the cytoplasm and sarcoplasmic reticulum. We also demonstrated the capability of the model to reproduce potentially pro-arrhythmic dynamics under perturbed conditions, pertaining to calcium-transient alternans and spontaneous release events. Comparison with idealised cell models emphasised the importance of structure in determining calcium gradients and controlling the spatial dynamics associated with calcium-transient alternans, wherein the probabilistic nature of dyad activation and recruitment was constrained. The model was further used to highlight the criticality in calcium spark propagation in relation to inter-dyad distances. The model presented provides a powerful tool for future investigation of structure-function relationships underlying physiological and pathophysiological intracellular calcium handling phenomena at the whole-cell. The approach allows for the first time direct integration of high-resolution images of 3-D intracellular structures with models of calcium cycling, presenting the possibility to directly assess the functional impact of structural remodelling at the cellular scale.

The organisation of the membrane and sub-cellular structures of cells in the heart closely controls the coupling between its electrical and mechanical function. Computational models of the cellular calcium handling system, which is responsible for this electro-mechanical coupling, have been developed in recent years to study underlying structure-function relationships. Previous models have been largely idealised in structure; we present a new model which incorporates experimental data describing the high-resolution organisation of the primary structures involved in calcium dynamics. Significantly, the structure of the intracellular calcium store is modelled for the first time. The model is shown to reproduce calcium dynamics in control cells in both normal and abnormal conditions, demonstrating its suitability for future investigation of structure-function relationships. Thus, the model presented provides a powerful tool for the direct integration of experimentally acquired structural data in healthy and diseased cells and assessment of the role of structure in regulating normal and abnormal calcium dynamics.

A computational study of the role of mitochondrial organization on cardiac bioenergetics

All cells in the body have a specific shape and internal organization which is specific to that cell's function. Heart cells are rod-shaped, and contain arrays of contractile protines (myofibrils) and mitochondria (organelles that produce energy) that are aligned along the length of the rod. This arrangement is presumed to allow the cell to generate maximal contractile force for each heartbeat and for energy metabolites to be readily available to generate this force. Heart disease phenotypes, such as diabetic cardiomyopathy and heart failure, exhibit altered organization of mitochondria. However, physiological and computational studies have predominantly investigated the effect of the biochemical changes that accompany the disease alone, such as reduced rates of ATP production by mitochondria. We present a modeling study that examines the effect of mitochondrial organization on energy metabolite distribution during the heartbeat. A 2D micrograph of the cell cross-section was selected from a 3D image stack of structural data of a cardiac cell. The image was used to generate a 2D finite element model, on which mitochondrial oxidative phosphorylation and energy metabolite diffusion was modelled. Results illustrate that mitochondrial density can induce heterogeneity in the distribution of metabolites across the cell area. We discuss the implications of these findings and avenues for future, more indepth studies.

Nanoscale Imaging of Collagen Gels with Focused Ion Beam Milling and Scanning Electron Microscopy

In vitro polymerized type I collagen hydrogels have been used extensively as a model system for three-dimensional (3D) cell and tissue culture, studies of fibrillogenesis, and investigation of multiscale force transmission within connective tissues. The nanoscale organization of collagen fibrils plays an essential role in the mechanics of these gels and emergent cellular behavior in culture, yet quantifying 3D structure with nanoscale resolution to fully characterize fibril organization remains a significant technical challenge. In this study, we demonstrate that a new imaging modality, focused ion beam scanning electron microscopy (FIB-SEM), can be used to generate 3D image datasets for visualizing and quantifying complex nanoscale organization and morphometry in collagen gels. We polymerized gels at a number of concentrations and conditions commonly used for in vitro models, stained and embedded the samples, and performed FIB-SEM imaging. The resulting image data had a voxel size of 25 nm, which is the highest resolution 3D data of a collagen fibril network ever obtained for collagen gels. This resolution was essential for discerning individual fibrils, fibril paths, and their branching and grouping. The resulting volumetric images revealed that polymerization conditions have a significant impact on the complex fibril morphology of the gels. We segmented the fibril network and demonstrated that individual collagen fibrils can be tracked in 3D space, providing quantitative analysis of network descriptors such as fibril diameter distribution, length, branch points, and fibril aggregations. FIB-SEM 3D reconstructions showed considerably less lateral grouping and overlap of fibrils than standard 2D SEM images, likely due to artifacts in SEM introduced by dehydration. This study demonstrates the utility of FIB-SEM for 3D imaging of collagen gels and quantitative analysis of 3D fibril networks. We anticipate that the method will see application in future studies of structure-function relationships in collagen gels as well as native collagenous tissues.

Post-Myocardial Infarction T-tubules Form Enlarged Branched Structures With Dysregulation of Junctophilin-2 and Bridging Integrator 1 (BIN-1)

BACKGROUND

Heart failure is a common secondary complication following a myocardial infarction (MI), characterized by impaired cardiac contraction and t-tubule (t-t) loss. However, post-MI nano-scale morphological changes to the remaining t-ts are poorly understood.

METHOD AND RESULTS

We utilized a porcine model of MI, using a nonlethal microembolization method to generate controlled microinfarcts. Using serial block face scanning electron microscopy, we report that post-MI, after mild left-ventricular dysfunction has developed, t-ts are not only lost in the peri-infarct region, but also the remnant t-ts form enlarged, highly branched disordered structures, containing a dense intricate inner membrane. Biochemical and proteomics analyses showed that the calcium release channel, ryanodine receptor 2 (RyR2), abundance is unchanged, but junctophilin-2 (JP2), important for maintaining t-t trajectory, is depressed (-0.5×) in keeping with the t-ts being disorganized. However, immunolabeling shows that populations of RyR2 and JP2 remain associated with the remodeled t-ts. The bridging integrator 1 protein (BIN-1), a regulator of tubulogensis, is upregulated (+5.4×), consistent with an overdeveloped internal membrane system, a feature not present in control t-ts. Importantly, we have determined that t-ts, in the remote region, are narrowed and also contain dense membrane folds (BIN-1 is up-regulated +3.4×), whereas the t-ts have a radial organization comparable to control JP2 is upregulated +1.7×.

CONCLUSIONS

This study reveals previously unidentified remodeling of the t-t nano-architecture in the post-MI heart that extends to the remote region. Our findings highlight that targeting JP2 may be beneficial for preserving the orientation of the t-ts, attenuating the development of hypocontractility post-MI.

Studying dyadic structure-function relationships: a review of current modeling approaches and new insights into Ca2+ (mis)handling

Excitation–contraction coupling in cardiac myocytes requires calcium influx through L-type calcium channels in the sarcolemma, which gates calcium release through sarcoplasmic reticulum ryanodine receptors in a process known as calcium-induced calcium release, producing a myoplasmic calcium transient and enabling cardiomyocyte contraction. The spatio-temporal dynamics of calcium release, buffering, and reuptake into the sarcoplasmic reticulum play a central role in excitation–contraction coupling in both normal and diseased cardiac myocytes. However, further quantitative understanding of these cells’ calcium machinery and the study of mechanisms that underlie both normal cardiac function and calcium-dependent etiologies in heart disease requires accurate knowledge of cardiac ultrastructure, protein distribution and subcellular function. As current imaging techniques are limited in spatial resolution, limiting insight into changes in calcium handling, computational models of excitation–contraction coupling have been increasingly employed to probe these structure–function relationships. This review will focus on the development of structural models of cardiac calcium dynamics at the subcellular level, orienting the reader broadly towards the development of models of subcellular calcium handling in cardiomyocytes. Specific focus will be given to progress in recent years in terms of multi-scale modeling employing resolved spatial models of subcellular calcium machinery. A review of the state-of-the-art will be followed by a review of emergent insights into calcium-dependent etiologies in heart disease and, finally, we will offer a perspective on future directions for related computational modeling and simulation efforts.

The cellular basis of platelet secretion: Emerging structure/function relationships

Platelet activation has long been known to be accompanied by secretion from at least three types of compartments. These include dense granules, the major source of small molecules; α-granules, the major protein storage organelle; and lysosomes, the site of acid hydrolase storage. Despite ~60 years of research, there are still many unanswered questions about the cell biology of platelet secretion: for example, how are these secretory organelles organized to support cargo release and what are the key routes of cargo release, granule to plasma membrane or granule to canalicular system. Moreover, in recent years, increasing evidence points to the platelet being organized for secretion of the contents from other organelles, namely the dense tubular system (endoplasmic reticulum) and the Golgi apparatus. Conceivably, protein secretion is a widespread property of the platelet and its organelles. In this review, we concentrate on the cell biology of the α-granule and its structure/function relationships. We both review the literature and discuss the wide array of 3-dimensional, high-resolution structural approaches that have emerged in the last few years. These have begun to reveal new and unanticipated outcomes and some of these are discussed. We are hopeful that the next several years will bring rapid advances to this field that will resolve past controversies and be clinically relevant.

Superresolution Imaging of Clinical Formalin Fixed Paraffin Embedded Breast Cancer with Single Molecule Localization Microscopy

Millions of archived formalin-fixed, paraffin-embedded (FFPE) specimens contain valuable molecular insight into healthy and diseased states persevered in their native ultrastructure. To diagnose and treat diseases in tissue on the nanoscopic scale, pathology traditionally employs electron microscopy (EM), but this platform has significant limitations including cost and painstaking sample preparation. The invention of single molecule localization microscopy (SMLM) optically overcame the diffraction limit of light to resolve fluorescently labeled molecules on the nanoscale, leading to many exciting biological discoveries. However, applications of SMLM in preserved tissues has been limited. Through adaptation of the immunofluorescence workflow on FFPE sections milled at histological thickness, cellular architecture can now be visualized on the nanoscale using SMLM including individual mitochondria, undulations in the nuclear lamina, and the HER2 receptor on membrane protrusions in human breast cancer specimens. Using astigmatism imaging, these structures can also be resolved in three dimensions to a depth of ~800 nm. These results demonstrate the utility of SMLM in efficiently uncovering ultrastructural information of archived clinical samples, which may offer molecular insights into the physiopathology of tissues to assist in disease diagnosis and treatment using conventional sample preparation methods.

Electrical defects of the transverse-axial tubular system in cardiac diseases

Electrical excitability is an essential feature of cardiomyocytes and the homogenous propagation of the action potential is guaranteed by a complex network of membrane invaginations called the transverse-axial tubular system (TATS). TATS structural remodelling is a hallmark of cardiac diseases and we demonstrated that this can be accompanied by electrical defects at single T-tubular level. Using a random-access multi-photon (RAMP) microscope, we found that pathological T-tubules can fail to conduct action potentials, which delays local Ca²⁺ release. Although the underlying causes for T-tubular electrical failure are still unknown, our findings suggest that they are likely to be related to local ultrastructural alterations. Here, we first review the experimental approach that allowed us to observe and dissect the consequences of TATS electrical dysfunction and then propose two different strategies to unveil the reasons for T-tubular electrical failures. The first strategy consists in a correlative approach, in which the failing T-tubule identified with the RAMP microscope is then imaged with electron microscopy. The second approach exploits the diffusion of molecules within TATS to gain insights into the local TATS structure, even without a thorough reconstruction of the tubular network. Although challenging, the local electrical failure occurring at single T-tubules is a fundamental question that needs to be addressed and could provide novel insights in cardiac pathophysiology.

Three-dimensional ultrastructural analysis of cells in the periodontal ligament using focused ion beam/scanning electron microscope tomography

The accurate comprehension of normal tissue provides essential data to analyse abnormalities such as disease and regenerative processes. In addition, understanding the proper structure of the target tissue and its microenvironment may facilitate successful novel treatment strategies. Many studies have examined the nature and structure of periodontal ligaments (PDLs); however, the three-dimensional (3D) structure of cells in normal PDLs remains poorly understood. In this study, we used focused ion beam/scanning electron microscope tomography to investigate the whole 3D ultrastructure of PDL cells along with quantitatively analysing their structural properties and ascertaining their orientation to the direction of the collagen fibre. PDL cells were shown to be in contact with each other, forming a widespread mesh-like network between the cementum and the alveolar bone. The volume of the cells in the horizontal fibre area was significantly larger than in other areas, whereas the anisotropy of these cells was lower than in other areas. Furthermore, the orientation of cells to the PDL fibres was not parallel to the PDL fibres in each area. As similar evaluations are recognized as being challenging using conventional two-dimensional methods, these novel 3D findings may contribute necessary knowledge for the comprehensive understanding and analysis of PDLs.

Analysis of Brain Mitochondria Using Serial Block-Face Scanning Electron Microscopy

Human brain is a high energy consuming organ that mainly relies on glucose as a fuel source. Glucose is catabolized by brain mitochondria via glycolysis, tri-carboxylic acid (TCA) cycle and oxidative phosphorylation (OXPHOS) pathways to produce cellular energy in the form of adenosine triphosphate (ATP). Impairment of mitochondrial ATP production causes mitochondrial disorders, which present clinically with prominent neurological and myopathic symptoms. Mitochondrial defects are also present in neurodevelopmental disorders (e.g. autism spectrum disorder) and neurodegenerative disorders (e.g. amyotrophic lateral sclerosis, Alzheimer's and Parkinson's diseases). Thus, there is an increased interest in the field for performing 3D analysis of mitochondrial morphology, structure and distribution under both healthy and disease states. The brain mitochondrial morphology is extremely diverse, with some mitochondria especially those in the synaptic region being in the range of <200 nm diameter, which is below the resolution limit of traditional light microscopy. Expressing a mitochondrially-targeted green fluorescent protein (GFP) in the brain significantly enhances the organellar detection by confocal microscopy. However, it does not overcome the constraints on the sensitivity of detection of relatively small sized mitochondria without oversaturating the images of large sized mitochondria. While serial transmission electron microscopy has been successfully used to characterize mitochondria at the neuronal synapse, this technique is extremely time-consuming especially when comparing multiple samples. The serial block-face scanning electron microscopy (SBFSEM) technique involves an automated process of sectioning, imaging blocks of tissue and data acquisition. Here, we provide a protocol to perform SBFSEM of a defined region from rodent brain to rapidly reconstruct and visualize mitochondrial morphology. This technique could also be used to provide accurate information on mitochondrial number, volume, size and distribution in a defined brain region. Since the obtained image resolution is high (typically under 10 nm) any gross mitochondrial morphological defects may also be detected.

Electron tomography of rabbit cardiomyocyte three-dimensional ultrastructure

The field of cardiovascular research has benefitted from rapid developments in imaging technology over the last few decades. Accordingly, an ever growing number of large, multidimensional data sets have begun to appear, often challenging existing pre-conceptions about structure and function of biological systems. For tissue and cell structure imaging, the move from 2D section-based microscopy to true 3D data collection has been a major driver of new insight. In the sub-cellular domain, electron tomography is a powerful technique for exploration of cellular structures in 3D with unparalleled fidelity at nanometer resolution.

Electron tomography is particularly advantageous for studying highly compartmentalised cells such as cardiomyocytes, where elaborate sub-cellular structures play crucial roles in electrophysiology and mechanics. Although the anatomy of specific ultra-structures, such as dyadic couplons, has been extensively explored using 2D electron microscopy of thin sections, we still lack accurate, quantitative knowledge of true individual shape, volume and surface area of sub-cellular domains, as well as their 3D spatial interrelations; let alone of how these are reshaped during the cycle of contraction and relaxation. Here we discuss and illustrate the utility of ET for identification, visualisation, and analysis of 3D cardiomyocyte ultrastructures such as the T-tubular system, sarcoplasmic reticulum, mitochondria and microtubules.

Nanoscale bio-platforms for living cell interrogation: current status and future perspectives

The living cell is a complex entity that dynamically responds to both intracellular and extracellular environments. Extensive efforts have been devoted to the understanding intracellular functions orchestrated with mRNAs and proteins in investigation of the fate of a single-cell, including proliferation, apoptosis, motility, differentiation and mutations. The rapid development of modern cellular analysis techniques (e.g. PCR, western blotting, immunochemistry, etc.) offers new opportunities in quantitative analysis of RNA/protein expression up to a single cell level. The recent entries of nanoscale platforms that include kinds of methodologies with high spatial and temporal resolution have been widely employed to probe the living cells. In this tutorial review paper, we give insight into background introduction and technical innovation of currently reported nanoscale platforms for living cell interrogation. These highlighted technologies are documented in details within four categories, including nano-biosensors for label-free detection of living cells, nanodevices for living cell probing by intracellular marker delivery, high-throughput platforms towards clinical current, and the progress of microscopic imaging platforms for cell/tissue tracking in vitro and in vivo. Perspectives for system improvement were also discussed to solve the limitations remains in current techniques, for the purpose of clinical use in future.

Examination of the Effects of Heterogeneous Organization of RyR Clusters, Myofibrils and Mitochondria on Ca2+ Release Patterns in Cardiomyocytes

Spatio-temporal dynamics of intracellular calcium, [Ca²⁺]_i, regulate the contractile function of cardiac muscle cells. Measuring [Ca²⁺]_i flux is central to the study of mechanisms that underlie both normal cardiac function and calcium-dependent etiologies in heart disease. However, current imaging techniques are limited in the spatial resolution to which changes in [Ca²⁺]_i can be detected. Using spatial point process statistics techniques we developed a novel method to simulate the spatial distribution of RyR clusters, which act as the major mediators of contractile Ca²⁺ release, upon a physiologically-realistic cellular landscape composed of tightly-packed mitochondria and myofibrils. We applied this method to computationally combine confocal-scale (~ 200 nm) data of RyR clusters with 3D electron microscopy data (~ 30 nm) of myofibrils and mitochondria, both collected from adult rat left ventricular myocytes. Using this hybrid-scale spatial model, we simulated reaction-diffusion of [Ca²⁺]_i during the rising phase of the transient (first 30 ms after initiation). At 30 ms, the average peak of the simulated [Ca²⁺]_i transient and of the simulated fluorescence intensity signal, F/F₀, reached values similar to that found in the literature ([Ca²⁺]_i ≈1 μM; F/F₀≈5.5). However, our model predicted the variation in [Ca²⁺]_i to be between 0.3 and 12.7 μM (~3 to 100 fold from resting value of 0.1 μM) and the corresponding F/F₀ signal ranging from 3 to 9.5. We demonstrate in this study that: (i) heterogeneities in the [Ca²⁺]_i transient are due not only to heterogeneous distribution and clustering of mitochondria; (ii) but also to heterogeneous local densities of RyR clusters. Further, we show that: (iii) these structure-induced heterogeneities in [Ca²⁺]_i can appear in line scan data. Finally, using our unique method for generating RyR cluster distributions, we demonstrate the robustness in the [Ca²⁺]_i transient to differences in RyR cluster distributions measured between rat and human cardiomyocytes.

Calcium (Ca²⁺) acts as a signal for many functions in the heart cell, from its primary role in triggering contractions during the heartbeat to acting as a signal for cell growth. Cellular function is tightly coupled to its ultra-structural organization. Spatially-realistic and biophysics-based computational models can provide quantitative insights into structure-function relationships in Ca²⁺ signaling. We developed a novel computational model of a rat ventricular myocyte that integrates structural information from confocal and electron microscopy datasets that were independently acquired and includes: myofibrils (protein complexes that contract during the heartbeat), mitochondria (organelles that provide energy for contraction), and ryanodine receptors (RyR, ion channels that release the Ca²⁺ required to trigger myofibril contraction from intracellular stores). Using this model, we examined [Ca²⁺]_i dynamics throughout the cell cross-section at a much higher resolution than previously possible. We estimated the size of structural maladaptation that would cause disease-related alterations in [Ca²⁺]_i dynamics. Using our methods for data integration, we also tested whether reducing the density of RyRs in human cardiomyocytes (~40% relative to rat) would have a significant effect on [Ca²⁺]_i. We found that Ca²⁺ release patterns between the two species are similar, suggesting Ca²⁺ dynamics are robust to variations in cell ultrastructure.

Telocytes and putative stem cells in ageing human heart

Tradition considers that mammalian heart consists of about 70% non-myocytes (interstitial cells) and 30% cardiomyocytes (CMs). Anyway, the presence of telocytes (TCs) has been overlooked, since they were described in 2010 (visit http://www.telocytes.com). Also, the number of cardiac stem cells (CSCs) has not accurately estimated in humans during ageing. We used electron microscopy to identify and estimate the number of cells in human atrial myocardium (appendages). Three age-related groups were studied: newborns (17 days–1 year), children (6–17 years) and adults (34–60 years). Morphometry was performed on low-magnification electron microscope images using computer-assisted technology. We found that interstitial area gradually increases with age from 31.3 ± 4.9% in newborns to 41 ± 5.2% in adults. Also, the number of blood capillaries (per mm²) increased with several hundreds in children and adults versus newborns. CMs are the most numerous cells, representing 76% in newborns, 88% in children and 86% in adults. Images of CMs mitoses were seen in the 17-day newborns. Interestingly, no lipofuscin granules were found in CMs of human newborns and children. The percentage of cells that occupy interstitium were (depending on age): endothelial cells 52–62%; vascular smooth muscle cells and pericytes 22–28%, Schwann cells with nerve endings 6–7%, fibroblasts 3–10%, macrophages 1–8%, TCs about 1% and stem cells less than 1%. We cannot confirm the popular belief that cardiac fibroblasts are the most prevalent cell type in the heart and account for about 20% of myocardial volume. Numerically, TCs represent a small fraction of human cardiac interstitial cells, but because of their extensive telopodes, they achieve a 3D network that, for instance, supports CSCs. The myocardial (very) low capability to regenerate may be explained by the number of CSCs, which decreases fivefold by age (from 0.5% to 0.1% in newborns versus adults).