Excitation–contraction coupling in vertebrate hearts is underpinned by calcium (Ca2+) release from Ca2+ release units (CRUs). CRUs are formed by clusters of channels called ryanodine receptors on the sarcoplasmic reticulum (SR) within the cardiomyocyte. Distances between CRUs influence the diffusion of Ca2+, thus influencing the rate and strength of excitation–contraction coupling. Avian myocytes lack T-tubules, so Ca2+ from surface CRUs (peripheral couplings, PCs) must diffuse to internal CRU sites of the corbular SR (cSR) during centripetal propagation. Despite this, avian hearts achieve higher contractile rates and develop greater contractile strength than many mammalian hearts, which have T-tubules to provide simultaneous activation of the Ca2+ signal through the myocyte. We used 3D electron tomography to test the hypothesis that the intracellular distribution of CRUs in the avian heart permits faster and stronger contractions despite the absence of T-tubules. Nearest edge–edge distances between PCs and cSR, and geometric information including surface area and volume of individual cSR, were obtained for each cardiac chamber of the white leghorn chicken. Computational modelling was then used to establish a relationship between CRU distance and cell activation time in the avian heart. Our data suggest that cSR clustered close together along the Z-line is vital for rapid propagation of the Ca2+ signal from the cell periphery to the cell centre, which would aid in the strong and fast contractions of the avian heart.
The four-chambered avian heart possesses many similarities to its mammalian counterpart. Both groups independently evolved a fully divided ventricle, capable of producing fast contractile rates with robust pressure development. However, despite similarities in cardiac performance and gross organ morphology, there are substantial differences on a cellular and subcellular level within the cardiomyocytes. Avian atrial and ventricular myocytes are long (>100 µm) and thin (3–9 µm), with a small cross-sectional area (∼56 µm2) and cell volume (∼10 pl); these features provide a large surface area to volume ratio (Akester, 1981; Bogdanov et al., 1995; Dzialowski and Crossley, 2015). It is likely that these morphological features allow excitation–contraction coupling to work within the necessary timescale without sarcolemmal invaginations, known as T-tubules (Sommer and Johnson, 1969). T-tubules are found in all mammalian ventricular myocytes and in atrial myocytes of large mammals, and are important for synchronised depolarisation in these larger cells (Dibb et al., 2009). The lack of T-tubules, characteristic of ectothermic vertebrates (fish, amphibians and non-avian reptiles), has been associated with slower heart rates and less robust contractile function (Shiels and Galli, 2014). We hypothesise that the subcellular organisation of calcium (Ca2+) release units (CRUs) within the avian myocyte may reconcile this apparent enigma.
Contraction and relaxation of the heart is underpinned by Ca2+ cycling in the cardiomyocytes. Following membrane depolarisation, Ca2+ enters the cell via voltage-gated L-type Ca2+ channels (LTCC). The initial Ca2+ that enters via LTCCs induces further Ca2+ release from the intracellular stores of the sarcoplasmic reticulum (SR), in a process called Ca2+-induced Ca2+ release (CICR). SR Ca2+ is released by the ryanodine receptor (RyR) channels, which cluster in structures defined as CRUs. In birds, there are two types of CRU. The first are peripheral couplings (PCs), which are clusters of RyRs on the surface of junctional SR (jSR), directly adjacent to sarcolemmal LTCCs. The second type of CRU is formed by RyR clusters which are not associated with the surface sarcolemma, known as corbular SR (cSR). cSR are shorter, rounder segments of SR than those found at the PCs. Despite being positioned micrometres from the sarcolemma (Sommer, 1995), these internal CRUs contribute to global Ca2+ release and are necessary for robust excitation–contraction coupling (Franzini-Armstrong et al., 2005). Ca2+ released at PCs must diffuse centripetally to the cSR/EjSR. Ca2+ is taken up by contracting myofilaments and adjacent mitochondria, or can be buffered in the cytosol, and thus [Ca2+] falls rapidly the further it travels from a CRU (Sobie et al., 2006). The distribution of CRUs, as well as the depletion of SR Ca2+ stores, ensures that excitation–contraction coupling is not endlessly regenerative and unstable, a situation that would be pro-arrhythmic and deleterious to the myocyte (Cannell et al., 1995).
There is a limited amount of functional data on Ca2+ flux pathways in avian cardiomyocytes. [3H]ryanodine binding studies in pigeon and ﬁnch left ventricle (LV) indicate that the density and Ca2+ sensitivity of RyRs are similar to those of mammalian RyRs (Junker et al., 1994). An early study showed a high density of LTCC current and a prominent T-type Ca2+ channel current in whole-cell voltage-clamped finch ventricular myocytes (Bogdanov et al., 1995), which the authors suggested could aid in CICR in the absence of T-tubules. This same study reported faster inactivation kinetics of ICa in finch compared with rat myocytes and suggested that this may facilitate fast heart rates.
We are unaware of studies that quantify the dynamics of avian intracellular Ca2+ movement directly (Kim et al., 2000), but it is possible to speculate on this by examining Ca2+ cycling in other cells that are similar in structure. Ectotherm cardiomyocytes, mammalian neonatal cardiomyocytes and mammalian nodal and Purkinje cells also lack T-tubules and tend to be thinner in diameter than mammalian ventricular cardiomyocytes. Studies from these cell types show that Ca2+ levels rise rapidly at the periphery, followed by a time-dependent rise in the cell centre (Boyden et al., 2000; Woo et al., 2003; Shiels and White, 2005; Stuyvers et al., 2005; Louch et al., 2015).
calcium-induced calcium release
calcium release unit
distance between PCs
distance between Z-lines
distance between cSR along a Z-line
extended junctional SR
L-type calcium channel
transmission electron microscopy
Our current understanding of the structural organisation of the Ca2+ release system in birds discussed above has been characterised using 2D transmission electron microscopy (TEM) (Franzini-Armstrong et al., 1999; Perni et al., 2012). However, a more realistic representation of the CRU distribution in space can be achieved with 3D electron microscopy. In this study, we characterised the distribution of CRUs in myocytes from each of the four chambers of the heart of the white leghorn chicken (Gallus gallus domesticus) using electron tomography (ET). This technique enables the reconstruction and visualisation of subcellular organisation of the SR network in a detailed 3D structure. Tomograms from both atria and both ventricles were reconstructed, and the structures of interest were segmented to (1) study the nearest edge–edge distances between PCs and between the cSR; and (2) obtain cSR volume and surface area to determine Ca2+ capacity and RyR cluster size. These inter-CRU distances were used to inform a computer model of Ca2+ wave dynamics which tested the effect of inter-CRU distance on whole-cell Ca2+ activation.
MATERIALS AND METHODS
Tissue samples and specimen preparation
Three adult chickens [white leghorn variety, Gallus gallus domesticus (Linnaeus 1758); 1.5–2 kg in body mass] were acquired from Hinchliffe's Farm, Huddersfield, UK. White leghorns are a slow-growing bird, bred for egg production with no pre-disposition to cardiac dysfunction (Mirsalimi et al., 1993). Chickens were transported to the University of Manchester Biological Services Facility, where they were humanely killed with pentobarbital, followed by dislocation of the neck in accordance with the Animals (Scientific Procedures) Act 1986. Hearts were excised, and tissue from each of the four chambers (left atrium, LA; right atrium, RA; left ventricle, LV; right ventricle, RV) was cut into 1 mm3 samples and placed into Karnovsky fixative with CaCl2 (2% formaldehyde, 2.5% glutaraldehyde, 50 mmol l−1 CaCl2, 0.1 mol l−1 Hepes). Tissue preparation was according to the Ellisman protocol (Deerinck et al., 2010). All reagents were obtained from Sigma.
Sectioning was carried out on a Reichert-Jung UltraCut E ultramicrotome using a glass knife. Sections of ∼100 nm thickness (gold coloured ribbons) were produced for observation using a FEI Tecnai 12 Bio Twin transmission electron microscope operated at 100 kV to assess tissue preparation prior to ET.
ET was performed using an FEI Tecnai G2 Polara transmission electron microscope at the University of Manchester Electron Microscopy Facility, operated at 300 kV using magnifications of approximately ×10,000, ×12,000 and ×15,500. For ET, sections of ∼400 nm thickness (green/purple coloured ribbons) were collected on 200 mesh grids; 10 nm gold fiducial markers were added to the sections to aid tilt series collection, alignment and dataset reconstruction. Tilt series were taken at areas where the myofibrils were longitudinal, as this aids the identification of the cSR. Single-axis tilt series were acquired at 1 deg intervals from −60 to +60 deg, or as close to this maximal angle range as possible. Tilt series were obtained from multiple cells within each of the atria and ventricles for all individuals. Alignment, reconstruction and segmentation were performed in the open-source software package IMOD (Kremer et al., 1996). eTomo was used for alignment and reconstruction, and 3dmod was used to view the tomograms, segment the structures of interest and obtain the distance, volume and surface area measurements from the 3D reconstructions.
Following manual segmentation, the edge–edge distances between the CRUs were measured. PCs were identified as specialised regions of SR in close apposition (10–15 nm) to the sarcolemma (Junker et al., 1994), while cSR were identified as spherical structures (∼100 nm in diameter) attached to the SR network, predominantly found along Z-lines (Asghari et al., 2009). Nearest edge–edge distances between PCs and between cSR were obtained from IMOD by drawing a line in the tomogram volume between the segmented CRUs. Nearest edge–edge distances are important for understanding cellular Ca2+ dynamics as they indicate the minimum distance that must be crossed by Ca2+ released at one site in order to act on an adjacent site (Perni et al., 2012). As cSR are distributed around the Z-lines, we measured distances between cSR distributed around the same Z-line. Geometric data regarding individual cSR were obtained by fully segmenting each cSR in every slice of each tomogram.
A 2D Ca2+ wave model representing an avian cardiomyocyte was constructed to assess the relationship between inter-CRU distance and whole-cell Ca2+ activation time. The myocyte was constructed as 8 μm wide (y-axis) and 136 μm long (x-axis) (Kim et al., 2000). PCs were placed at regular intervals along the cell membrane in the x-direction; Z-lines were placed at regular intervals along the x-axis and extended across the y-axis; and cSR were placed at regular intervals on the Z-lines (see Fig. 5A). The model geometry was discretised at a space step of 0.05 µm and thus PCs were placed at a depth from the cell membrane of 50 nm. This is greater than the biological distance (∼10–15 nm) but we were required to accept this limitation because of the computational power required to provide finer discretisation. CRUs were separated by distances d1 (distance between PCs), d2 (distance between Z-lines) and d3 (distance between cSR along a Z-line). These distances were varied in a systematic way to encompass the range of values measured in the tomogram datasets. Distance d1 was tested at 0.25, 0.5 and 0.75 µm, d3 at 0.1, 0.2, 0.4 and 0.6 µm, and d2 was varied to represent Z-line spacing indicative of a myocyte at rest (1.7–1.9 µm), as well as encompassing values that could be achieved during myocardial stretch (2.1 µm) and myocardial contraction (1.5 µm). The distance ranges are tabulated in Table S1.
Ca2+ dynamics model
The Ca2+ dynamics model was adapted from previous models (Cheng et al., 1993; Smith et al., 1998; Izu et al., 2001). CRU activation began at the surface sarcolemma and moved towards the centre of the cell. Cell membrane activation initiated propagation of the Ca2+ wave in the y-direction. At t=0, [Ca2+] was raised to 100 μmol l−1 at y=0 μm and y=8 μm for 1 ms. Ca2+ dynamics were then permitted to evolve for the duration of one heartbeat (225 ms), according to the equations given below (see ‘Model equations for Ca2+ wave propagation simulations’). Whole-cell activation was achieved when all grid locations reached a [Ca2+] value of 100 µmol l−1 or more at least once (as per Smith et al., 1998; Izu et al., 2001). Ca2+ release was simulated as a stochastic process; release at a CRU occurred when local [Ca2+] exceeded 0.1 µmol l−1 as described previously (Smith et al., 1998). The probability of a given CRU being open was assumed to be proportional to Jpump, the Ca2+ resequestration flux by Ca2+ ATPases located on the CRUs (Smith et al., 1998) (see below). Once a CRU was open, it was permitted to release Ca2+ for 10 ms (Smith et al., 1998). The spatial network of CRUs were diffusively coupled to permit simulation of Ca2+ waves (Izu et al., 2001).
The effect of mobile buffers was omitted. This may affect the absolute slope of the relationship between Ca2+ activation time and the distance between CRUs, but is not expected to impact the relative effect of changing distance (d1, d2 and d3) on activation time.
Robust implicit finite difference solvers developed previously (Kharche et al., 2017) were used to solve the reaction–diffusion equations using a maximum temporal step of 0.02 ms (Table S1). Ten simulations were performed for each combination of distances (d1, d2, d3), and an average activation time from the 10 simulations was taken to form each data point. Parallelisation was implemented using message passing interface (MPI). Each simulation required 24 processors 4 h to complete.
Model equations for Ca2+ wave propagation simulations
Statistics were performed using GraphPad Prism with unpaired t-tests and one-way ANOVA, with Tukey post hoc analysis as specified in the figure legends. The threshold for statistical significance was P<0.05. The results are expressed as mean±s.e.m., with the number of measurements, from the number of tissue sections and the number of animals, provided in each legend. The spread of data is shown using scatter plots of individual data points with mean values overlaid.
Segmented model for calculating distance
Each tomogram is a reconstruction of an approximately 400 nm thick section of avian cardiac tissue (Fig. 1; Movie 1). When the tissue is sectioned longitudinally with respect to the axis of the cells, it is apparent that cSR (yellow) are localised at the Z-lines and thus are separated at roughly the length of the sarcomere (Fig. 1A,B). The 3D structural model (Fig. 1C; Movie 2) gives a representative display of cSR spread along Z-lines, as well as the entire SR network adjacent to the sarcolemma.
PCs are defined as clusters of RyRs on the surface of the jSR which are directly adjacent (10–15 nm) to LTCCs on the sarcolemma (Fig. 2A,B). The mean (±s.e.m.) nearest edge–edge distances between PCs in each of the four chambers of the heart were as follows: LA 377±19 nm (n=380), RA 347±22 nm (n=228), LV 334±26 nm (n=193), RV 462±30 nm (n=195) (n=number of individual distances measured; Fig. 2C). Distances between PCs in the RV were significantly greater (one-way ANOVA, P<0.05) than those measured in the other cardiac chambers.
The internal CRUs in avian hearts are the cSR, which are spherical structures of approximately 100 nm diameter, typically found staggered along the Z-line (Fig. 3A,B; Movie 3). cSR are notably larger in size than the network tubules of free SR (fSR); however, they are highly polymorphous in their geometry, varying in shape and size (described below). The mean (±s.e.m.) nearest edge–edge distances between cSR at the same Z-line (Fig. 3C) in each of the four chambers of the heart were as follows: LA 423±16 nm (n=286), RA 501±20 nm (n=189), LV 485±24 nm (n=133), RV 465±23 nm (n=194). Distances measured in the LA were significantly shorter than those in the RA but not different from those of the two ventricles (one-way ANOVA, P<0.05).
Geometric models for individual cSR
ET enables the 3D rendering of structures of interest (Fig. 4A,B). Thus, diameter (Fig. 4C), surface area (Fig. 4D) and volume (Fig. 4E) of cSR in each chamber of an individual bird heart were measured to provide insight into their role in Ca2+ release (Table S2). The scatterplots in Fig. 4C–E highlight the polymorphic nature of the geometry of bird cSR. The LA cSR were smaller (one-way ANOVA, P<0.05) in volume and surface area than those of the other chambers.
Computational modelling Ca2+ dynamics
To understand how inter-CRU distance may impact on Ca2+ activation time in an avian myocyte devoid of T-tubules, we constructed a spatially extended 2D Ca2+ wave model. In the model, CRUs were placed at locations separated by distances d1 (distance between PCs), d2 (distance between Z-lines) and d3 (distance between cSR along a Z-line) (Fig. 5A). Frames from a simulation of Ca2+ waves show the mode of diffusion across the cell, with [Ca2+] initially increasing at the periphery before diffusion into the interior along the Z-lines (Fig. 5B).
CRU distance in the model was varied in a systematic way to encompass the range of values measured in the tomograms, in order to deduce the effect on whole-cell activation time (Fig. 6). This analysis suggests that distance between PCs (d1) has little effect on whole-cell Ca2+ activation time (Fig. 6A). In contrast, when the distance between cSR (d3) is increased from 0.1 µm to 0.2 µm, activation time is delayed by approximately 3.5 ms; as d3 is increased to 0.4 µm, activation time is more than 10 ms longer than that at 0.1 µm (Fig. 6B). At larger d3 distances (beyond those measured between cSR within a Z-line in the current study), activation time is slowed even further. Under these conditions, widely separated PCs (d1) compound the effect and activation time requires approximately 20 ms. In general, activation time is unaffected as the distance between Z-lines is varied (d2) in a manner anticipated during sarcomeric stretch and contraction (Fig. 6). Indeed, because distance d3 was considerably smaller than distance d2, the CRUs along the Z-lines activated prior to Ca2+ passively diffusing to neighbouring Z-lines that were not proximal to the PCs.
Aves occupy a unique position in vertebrate evolution. They possess a four-chambered heart with a fully divided ventricle, which arose independently from that in mammals. However, despite gross structural similarities to mammalian hearts, at a cellular level, the avian cardiomyocyte more closely resembles that of (non-avian) reptiles. This study is the first to use ET to investigate the subcellular distribution of the intracellular Ca2+ release machinery in the avian heart. We measured distances between PCs and between cSR along a Z-line, and found the PCs to be closer together and the cSR further apart than those reported previously using 2D TEM (Franzini-Armstrong et al., 1999; Perni et al., 2012). Our study has also revealed the diversity of cSR volumes and the varied distribution of these structures in all four cardiac chambers. We used a computational approach to test how varying the distance between these subcellular CRUs affects Ca2+ activation time. Our model highlights the importance of the distance between cSR along a Z-line (d3). Together, our findings suggest that cSR clustered close together along the Z-line are vital for rapid propagation of the Ca2+ signal from the cell periphery to the cell centre, and facilitate strong and fast contractions.
CRUs in the form of PCs have been identified in the hearts of all vertebrates studied, except in the frog ventricle, and are most often associated with the Z-lines (Shiels and Galli, 2014). The mean distance between nearest neighbour PCs in the white leghorn chicken ranged from 334 nm in the LV to 462 nm in the RV, which is comparable to previous values measured using 2D TEM in the LV [472 nm (Franzini-Armstrong et al., 1999) and 567 nm (Perni et al., 2012)]. We found similar distances between PCs in the LV, RA and LA but slightly longer distances in the RV. The functional significance of this pattern is unclear. Size and frequency of peripheral CRUs vary between vertebrate species and between cardiac chambers within a species, and in some studies this variability has been related to the efficacy of excitation–contraction coupling (Perni et al., 2012; Shiels and Galli, 2014). Indeed, animals with high heart rates (i.e. finch and rat, resting heart rate ∼300–350 beats min−1) have closer PCs than animals with slower heart rates (Perni et al., 2012) [i.e. chicken, resting heart rate ∼200 beats min−1; lizard and fish, heart rate dependent on temperature but generally below 120 beats min−1 (Farrell, 1991)]. The distances measured in our study support previous work in chicken (Perni et al., 2012), which suggests PCs are activated by adjacent LTCCs in the sarcolemmal membrane directly, with little spread of activating Ca2+ longitudinally between PCs. Indeed, our simulations suggest that changing the distance between PCs between 250 and 750 nm has little impact on whole-cell activation, particularly when the distances between cSR are narrow (Fig. 6). This differs from mammalian atrial myocytes and from birds with faster heart rates and/or contractile force (i.e. finch; Perni et al., 2012), where extensive and propagative CICR between neighbouring CRUs at the cell periphery has been observed (Chen-Izu et al., 2006). In the finch heart, closely packed PCs are further aided by a higher density of LTCCs (Bogdanov et al., 1995) to achieve fast heart rates.
The CRUs of the avian cSR can be thought of as Ca2+ release relay stations, carrying the wave of Ca2+ from the myocyte periphery to the cell centre in the absence of T-tubules. We found the average distance between cSR at the same Z-line to range from 422 nm in the LA, to 500 nm in the RA, which are significantly larger than those reported by Franzini-Armstrong's group for the LV [148 nm (Franzini-Armstrong et al., 1999) and 235 nm (Perni et al., 2012)] using 2D TEM images taken in a transverse section. Franzini-Armstrong et al. (1999) state that ‘larger distances were ignored’ in their measurements, so although the exact threshold used for excluding large distances is not specified, it may explain the discrepancy between studies. The heterogeneous distribution of organelles (e.g. mitochondria), combined with changes in the distance between organelles during the contraction–relaxation cycle are both features that help prevent circular, endless propagation of a Ca2+ wave. Indeed, our computational model indicates that the distance between cSR has a dramatic effect on the activation time. As cSR distance is varied from those measured previously (∼150–235 nm) to those measured in our study (∼400–500 nm), activation time more than doubles (Fig. 6). Thus, chains of closely distributed cSR at a Z-line increase the probability that Ca2+ will be able to diffuse and activate neighbouring units.
We provide the first 3D geometrical description of cSR in chicken hearts. A similar diverse and polymorphic description has been detailed for the jSR and the T-tubule lattice in mouse cardiomyocytes (Hayashi et al., 2009). Average diameters for cSR from each of the four chambers ranged from 120 to 130 nm, a parameter that aided cSR identification. The average volume of a cSR in each of the four chambers ranged from 446,000 to 600,000 nm3, but the largest and smallest ranged from 120,000 to 1,410,000 nm3. The average surface area from each of the four chambers ranged from 22,000 to 28,000 nm2 but the extremes were 9000 nm2 and 60,000 nm2. Thus, there is extreme heterogeneity in the avian cardiomyocyte's Ca2+ release structures.
We were unable to resolve individual RyR positions on the cSR in this study. However, given the average surface area of cSR was between 22,000 and 28,000 nm2, and knowing that a single RyR is 29×29 nm2, or 841 nm2 (Chen-Izu et al., 2006; Baddeley et al., 2009), we calculate that the average chicken cSR could support between 26 and 33 RyR molecules. RyRs probably do not occupy the entire cSR surface; for instance, the dyad is only partially filled, and RyR clusters often present a checkerboard appearance (Baddeley et al., 2009; Asghari et al., 2014). We therefore estimate the average cSR possesses approximately 13–16 RyRs.
Calculations concerning jSR Ca2+ capacity (whereby jSR of volume 7×10−4 µm3 contains 21,000 Ca2+ ions; Sobie et al., 2002), lead us to suggest the average avian cSR would contain between 13,380 and 18,000 Ca2+ ions. The minimum and maximum volumes of cSR found here would thus contain 3600 and 42,300 Ca2+ ions, respectively. A 1 pA current releases around 3000 Ca2+ ions per millisecond (Sobie et al., 2002), and thus the total capacity of the cSR analysed in this study permits Ca2+ release for ∼5 ms in total before the SR is depleted of Ca2+. SR Ca2+ available for release at these CRUs is replenished by SERCA Ca2+ uptake along the fSR network, which diffuses back towards the junctional area containing the RyRs. Thus, the total capacity of the cSR does not finitely dictate the duration of Ca2+ release. It is likely that the total capacity can be released in the order of tens of milliseconds, correlating with previous measurements for known spark durations (Sobie et al., 2002; Stern et al., 2013) although the actual spark duration in avian cardiomyocytes is not known. These estimates correlate with the Ca2+ activation time of our model, which ranges between 5 and 25 ms, with a mean of ∼15 ms, for the whole cell to be activated, depending on CRU spacing.
Ca2+ dynamics model
We used a spatially extended 2D cell model to improve the interpretation of our experimental data. The model showed that the variation of distances between cSR within a Z-line supersedes that of distances between PCs at the cell membrane for setting Ca2+ activation time of the cell. Our model has limitations that need to be considered when interpreting our results. First, stimulation was assumed to be generated by a rapidly propagating electrical pulse along the length of the cardiomyocyte membrane. Therefore, the initiated Ca2+ wave had an inward propagation. Thus, barring any numerical and boundary condition effects, it may be expected that PCs played a limited role in the output of whole-cell Ca2+ activation. Our model also incorporated local stochastic Ca2+ release, or sparks. The Ca2+ released by the sparks diffused in the cytoplasm uniformly in all directions; thus, the Ca2+ diffusing from a cSR would stimulate release from the directly adjacent CRU along the same Z-line prior to stimulating other CRUs. In previous studies, inter-CRU distance was assumed to be uniform in both the x- and y-direction, thus the work here extends these earlier models (Izu et al., 2001, 2013). However, in our model the locations of the CRUs were assumed to be symmetric, which aligns with the positioning along the Z-lines that we observed experimentally in the bird heart. However, stochasticity in the spatial distribution of CRUs has been observed (Qu et al., 2014), and will be included in future studies. Lastly, because of the limited functional studies on Ca2+ dynamics in avian myocytes (Kim et al., 2000), the model was parameterised based on mammalian data. Further functional studies are need to accurately measure Ca2+ flux in avian cardiomyocytes. Improved numerical schemes and solvers are also required to permit exploration of the experimental data in more depth. However, peak systolic Ca2+ was reached in a similar time frame and at a similar level in our model and that in mammals (Bers, 2002), despite the difference in cell structure between species.
Perspectives for avian cardiac performance
Domestic chickens have an average heartbeat of 218 beats min−1 (Prosheva et al., 2015). A single cardiac cycle is ∼0.275 s and composed of the following electrocardiogram components: P wave (atrial depolarisation) ∼30 ms, PQ interval ∼60 ms (AV conduction), QRS (ventricular depolarisation) ∼30 ms, and QT interval (duration of the ventricular action potential) ∼140 ms (Prosheva et al., 2015). Ca2+ release follows the depolarising wave and thus our mean Ca2+ activation time of 12–20 ms is within a reasonable range. In the scenario where cSR distance is at or beyond the level we observed in our tomographic data, activation time would be closer to 25–30 ms. However, in vivo, the rate of Ca2+ removal and SR replenishment will also influence the speed of the calcium cycle, as an increase/decrease in heart rate is predominantly dictated by the duration of the interval between T and P waves (Dzialowski and Crossley, 2015). Functional studies are required in avian myocytes to reveal the complexities of Ca2+ signalling in this cell type; for instance, the speed of Ca2+ extrusion and specific channel currents.
In birds with faster heart rates, such as the finch or hummingbird, rapid cardiac contraction relies on even faster relaying of the Ca2+ signal. This is brought about by closer spacing of CRUs, and in the case of the finch heart, the existence of EjSR in place of cSR (Perni et al., 2012). EjSR are more elongated CRUs, which might offer a greater capacity for Ca2+ storage, increased numbers of RyRs per cluster for Ca2+ release, or a better geometry for the spread of signal to neighbouring sites.
Our measurements of distance between PCs and between cSR, coupled with computational simulations, suggest that in chicken cardiomyocytes the Ca2+ transient would be initiated at the periphery, be large in concentration around this sub-sarcolemmal space, and diffuse toward the interior along chains of cSR. The resulting effect would be spatial inhomogeneity and non-synchronous spread of Ca2+ across the whole myocyte. However, upon activation of many CRUs, given that the Ca2+ flux is of large enough magnitude, diffusion between cSR positioned at adjacent Z-lines is possible, becoming more probable the closer they are to one another. This pattern is similar to that found in neonatal mammalian myocytes (Louch et al., 2015). As noted by others previously (Perni et al., 2012), the extensive and almost exclusive location of cSR along the Z-lines in bird myocytes induces a high degree of refractoriness to longitudinal Ca2+ wave propagation, with activation initiated independently at PCs cascading through the cSR. Thus, the organisation of CRUs and the short diffusional distance for Ca2+ transport in narrow cells allows for strong and fast contractions in avian myocytes and reinforces the connection between structural organisation of the myocyte, the CRUs, and the strength and rate of cardiac contraction across vertebrate classes.
The authors wish to thank the staff in the Faculty EM Facility for their assistance, in particular Tobias Starborg and Samantha Forbes, and the Wellcome Trust for equipment grant support to the EM Facility. We also thank Ben Newman for tissue samples, Dr Jonathan Codd for assistance with the birds, and Drs Katharine Dibb and Ashraf Kitmitto for useful discussions.
Conceptualization: C.P., H.A.S.; Methodology: T.M.D.S., S.K., C.P., H.A.S.; Formal analysis: T.M.D.S., S.K., C.P.; Resources: S.K., H.A.S.; Data curation: S.K.; Writing - original draft: T.M.D.S., H.A.S.; Writing - review & editing: T.M.D.S., S.K., C.P., H.A.S.; Visualization: T.M.D.S., C.P., H.A.S.; Supervision: C.P., H.A.S.; Project administration: H.A.S.; Funding acquisition: H.A.S.
H.A.S., C.P. and T.M.D.S. were supported by The University of Manchester. We thank Compute Canada and SharcNET Canada for providing computing resources to S.R.K.
The authors declare no competing or financial interests.