- 1Department of Clinical and Experimental Cardiology, Heart Center, Amsterdam Cardiovascular Sciences, Amsterdam UMC, University of Amsterdam, Amsterdam, Netherlands
- 2Department of Medical Biology, Amsterdam Cardiovascular Sciences, Amsterdam UMC, University of Amsterdam, Amsterdam, Netherlands
- 3Laboratory Genetic Metabolic Diseases, Amsterdam UMC, University of Amsterdam, Amsterdam Gastroenterology and Metabolism, Amsterdam Cardiovascular Sciences, Amsterdam, Netherlands
- 4Department of Pediatric Metabolic Diseases, Wilhelmina Children’s Hospital, University Medical Center Utrecht, Utrecht, Netherlands
- 5Institute of Pharmacology and Toxicology, Technische Universität Dresden, Dresden, Germany
- 6Department of Pediatric Metabolic Diseases, Emma Children's Hospital, Amsterdam UMC, University of Amsterdam, Amsterdam, Netherlands
Patients with a deficiency in very long-chain acyl-CoA dehydrogenase (VLCAD), an enzyme that is involved in the mitochondrial beta-oxidation of long-chain fatty acids, are at risk for developing cardiac arrhythmias. In human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs), VLCAD deficiency (VLCADD) results in a series of abnormalities, including: 1) accumulation of long-chain acylcarnitines, 2) action potential shortening, 3) higher systolic and diastolic intracellular Ca2+ concentrations, and 4) development of delayed afterdepolarizations. In the fatty acid oxidation process, carnitine is required for bidirectional transport of acyl groups across the mitochondrial membrane. Supplementation has been suggested as potential therapeutic approach in VLCADD, but its benefits are debated. Here, we studied the effects of carnitine supplementation on the long-chain acylcarnitine levels and performed electrophysiological analyses in VLCADD patient-derived hiPSC-CMs with a ACADVL gene mutation (p.Val283Ala/p.Glu381del). Under standard culture conditions, VLCADD hiPSC-CMs showed high concentrations of long-chain acylcarnitines, short action potentials, and high delayed afterdepolarizations occurrence. Incubation of the hiPSC-CMs with 400 µM L-carnitine for 48 h led to increased long-chain acylcarnitine levels both in medium and cells. In addition, carnitine supplementation neither restored abnormal action potential parameters nor the increased occurrence of delayed afterdepolarizations in VLCADD hiPSC-CMs. We conclude that long-chain acylcarnitine accumulation and electrophysiological abnormalities in VLCADD hiPSC-CMs are not normalized by carnitine supplementation, indicating that this treatment is unlikely to be beneficial against cardiac arrhythmias in VLCADD patients.
Introduction
Patients with a deficiency in very long-chain acyl-CoA dehydrogenase (VLCAD; EC 1.3.99.3), the enzyme catalyzing the first step of the mitochondrial beta-oxidation of long-chain fatty acids (Houten and Wanders, 2010; Knottnerus et al., 2018), are at risk for developing liver, skeletal, and heart muscle dysfunction (Ribas and Vargas, 2020; Wanders et al., 2020), including cardiac arrhythmias (Bonnet et al., 1999). Traditionally patients with VLCAD deficiency (VLCADD; OMIM 609575) are treated with dietary restriction of long-chain triglycerides, supplementation of medium-chain triglycerides, and prevention of catabolic state (Yamada and Taketani, 2019). Carnitine supplementation had also been proposed for patients with VLCADD in order to treat secondary carnitine deficiency and to increase the transport of acyl compounds out of the mitochondria (Treem et al., 1991; Winter, 2003).
Carnitine is required for the transport of activated long-chain fatty acids into the mitochondrial matrix, by formation of long-chain acylcarnitines (LCAC) from acyl-CoA esters, and is therefore essential for fatty acid oxidation (Houten and Wanders, 2010). In humans, carnitine is mainly derived from the diet, but can also be synthesized by the liver, kidney and brain (Vaz and Wanders, 2002; Reuter and Evans, 2012; Wanders et al., 2020). Tissues like skeletal muscle and the heart acquire carnitine from the circulation. Its indispensable role in metabolism is illustrated by patients who suffer from primary carnitine deficiency (OMIM 212140) due to excessive urinary carnitine wasting (Longo et al., 2006). So far, the evidence for a beneficial effect of carnitine supplementation in VLCADD patients is not conclusive (Nasser et al., 2012; (Spiekerkoetter et al., 2010; Ribas and Vargas, 2020). A major concern may be the potential increased formation in LCACs (Spiekerkoetter et al., 2009; Ribas and Vargas, 2020), especially since LCACs are associated with the development of cardiac arrhythmia (Corr et al., 1989). On the other hand, extremely low free carnitine levels, as seen in patients with primary carnitine deficiency, are also linked to abnormalities in electrical morphology of the heart (Roussel et al., 2016). Therefore, more insight on effects and safety of carnitine supplementation on electrophysiological abnormalities for VLCADD treatment is needed.
So far, the cardiac effects of carnitine supplementation have been studied in mice. In VLCADD murine models, it was found that supplementation of carnitine can indeed lead to a more pronounced accumulation of LCACs (Liebig et al., 2006; Primassin et al., 2008) without replenishment of free carnitine in cardiac tissue (Liebig et al., 2006). This is in contrast with the findings in a long-chain acyl-CoA dehydrogenase (LCAD) knock out (KO) murine model, in which carnitine supplementation did not lead to increased long-chain acylcarnitine levels and even normalized myocardial triglycerides (Bakermans et al., 2013). However, mice are considerably different from humans with respect to fatty acid oxidation as well as carnitine metabolism (Chegary et al., 2009), and therefore it is uncertain whether findings in these murine models can be directly extrapolated to the human situation and be applied to formulation of a human therapy. Very recently, we generated a human-based cardiac VLCADD model by generation of human induced pluripotent stem cells (hiPSCs) from patients carrying mutations in the VLCAD-encoding ACADVL gene and differentiated them into cardiomyocytes (hiPSC-CMs) (Knottnerus et al., 2020). We found that VLCAD deficiency in this hiPSC-CMs model results in accumulation of LCACs, higher systolic and diastolic intracellular Ca2+ (Ca2+i) concentrations, action potential (AP) shortening, and development of delayed afterdepolarizations (DADs). Furthermore, we found that pre-incubation of VLCADD hiPSC-CMs with either resveratrol or etomoxir–compounds that give rise to an enhanced mitochondrial biogenesis or inhibit fatty acid transport into the mitochondria, respectively (Lopaschuk et al., 2010; Houtkooper et al., 2012)– led to the normalization of the acylcarnitine levels and restored electrophysiological and Ca2+i abnormalities (Knottnerus et al., 2020). Importantly, these findings suggest that treatment with such compounds may be beneficial for VLCADD patients. In addition, these findings also indicate that our hiPSC-CM model is suitable for drug discovery studies of this metabolic disorder. In the present study, we used the VLCADD hiPSC-CM model to study the effects of carnitine supplementation on the LCAC profile and cellular electrophysiology.
Materials and Methods
Human Induced Pluripotent Stem Cell Derived Cardiomyocytes Generation
hiPSC-CMs were generated from the hiPSC line, iVLCADD1, which was derived from skin fibroblasts of a VLCADD-affected woman with mutations (p.Val283Ala/p.Glu381del) in the ACADVL gene (Knottnerus et al., 2020). Differentiation of the hiPSCs into hiPSC-CMs was performed in RPMI 1,640 medium supplemented with B27 (Gibco), initiated with CHIR99021, Activin A and BMP4, followed by Wnt-pathway inhibition by IWP4, differentiation proceeded in presence of RPMI medium supplemented with B27 and insulin until day 30 as we have previously described in detail (Knottnerus et al., 2020). Next, we performed a metabolic selection-based enrichment for hiPSC-CMs by applying glucose-depleted RPMI medium containing 4 mM lactate for six days, thereby removing a large proportion of non-cardiomyocytes (Tohyama et al., 2013). During the whole process, no serum was supplied to the culture medium. Finally, 36 days after the start of the differentiation process, hiPSC-CMs were dissociated to single cells using a two-step procedure starting with 5× TrypLE Select Enzyme (Gibco) to detach the cells from culture flasks, followed by a combination of liberase (Roche Chemicals) and Elastase (Serva) to obtain single cells (Knottnerus et al., 2020). The cells were seeded on matrigel-coated glass coverslips and cultured for 10 days in RPMI medium supplemented with 2% B27, 50 U/ml penicillin, and 50 µg/ml streptomycin and 48 h prior to measurement 400 µM L-carnitine (C0158, Sigma-Aldrich) was added to the culture medium in paired experiments. hiPSC-CMs generated from a control hiPSC line, iCTRL (Dudek et al., 2013), were used as control and were cultured in the absence of carnitine.
Acylcarnitine Profiling
The intracellular acylcarnitine profile of the hiPSC-CMs was assessed after 30 days of differentiation as described previously (Knottnerus et al., 2020). In short, hiPSC-CMs from the iCTRL and VLCADD1 line, here named iCTRL-CMs and iVLACDD1-CMs respectively, were harvested with trypsin and cell pellets or 50 µL of culture medium were extracted in 500 µL acetonitrile with 10 µL internal standards (2H3-C3-carnitine, 2H3-C8 carnitine, 2H3-C16 carnitine) and measured with high performance liquid chromatography mass spectrometry on a Q-ExactiveTM mass spectrometer (Thermo Scientific). Data was normalized to level of protein per well, determined using bicinchonic acid (BCA) assay with human serum albumin (HSA) in 0.4 mmol/L NaOH used as a standard (Knottnerus, et al., 2020).
Patch-Clamp Measurements
Data Acquisition
APs and afterdepolarizations were recorded at 36 ± 0.2°C using the perforated patch-clamp technique and an Axopatch 200B amplifier (Molecular Devices, Sunnyvale, CA, United States). Data acquisition, voltage control, and analysis were accomplished using custom software. Signals were low-pass filtered with a cut-off frequency of 5 kHz and digitized at 40 and 3 kHz for APs and DADs, respectively. Cell membrane capacitance (Cm, in pF) was estimated by dividing the time constant of the decay of the capacitive transient in response to 5 mV hyperpolarizing voltage clamp steps from −40 mV by the series resistance. Cm of hiPSC-CMs was 23.9 ± 1.3 pF (mean ± SEM, n = 41), and did not differed significantly between iCTRL-CMs, untreated iVLCADD1-CMs, and carnitine-treated iVLCADD1-CMs. Patch pipettes with a resistance of 2–3 MΩ were pulled from borosilicate glass (Harvard Apparatus) and filled with solution containing (in mM): 125 K-gluconate, 20 KCl, 5 NaCl, 0.44 Amphotericin-B, 10 HEPES; pH set to 7.2 (KOH). Cells were superfused with modified Tyrode’s solution containing (in mM): 140 NaCl, 5.4 KCl, 1.8 CaCl2, 1.0 MgCl2, 5.5 D-glucose, 5 HEPES; pH set to 7.4 (NaOH). All potentials were corrected for the estimated liquid junction potential of −15 mV (Barry and Lynch, 1991).
Action Potential Recordings
hiPSC-CMs offer a great human-based model for cardiac disease modeling (Hoekstra et al., 2012), and drug discovery and cardiotoxicity screenings (Mordwinkin et al., 2013), but a major limitation for electrophysiological studies is that they have a small or even complete lack of the inward rectifying K+ current (IK1) (Veerman et al., 2015). IK1 plays a major role in maintaining a stable resting membrane potential in cardiomyocytes (Nerbonne and Kass, 2005), and consequent to its absence, hiPSC-CMs have a depolarized maximal diastolic potential and are frequently spontaneously active. To overcome this limitation, we injected an in silico IK1 with kinetics of Kir2.1 channels and a 2 pA/pF current density through dynamic clamp, as previously described in detail (Meijer van Putten et al., 2015). The dynamic clamp injected IK1 resulted in quiescent hiPSC-CMs and we elicited APs at 0.2–4 Hz by 3 ms, ≈1.2× threshold current pulses through the patch pipette. Susceptibility to DADs was tested by applying a 3 Hz pacing episode (10 s) followed by an 8 s pause. After the pause, a single AP was evoked to test the inducibility of early afterdepolarizations (EADs). DADs were defined as depolarizations that occurred after the fast pacing and which were larger than 1 mV. EADs were defined as oscillations in membrane potential which interrupt or retard repolarization of the AP were defined as EADs accordingly to the definition of Cranefield (1977). The AP parameters analyzed were maximum diastolic potential (MDP, in mV), maximal AP amplitude (APA, in mV), maximum upstroke velocity (Vmax, in V/s) AP duration at 20, 50, and 90% repolarization (APD20, APD50, and APD90, respectively, in ms). Parameters from 10 consecutive APs were averaged and number of DADs and EADs was counted and averaged over five fast pacing recording traces.
Statistical Analysis
Statistical analysis was carried out with SigmaStat 3.5 software (Systat Software, Inc., San Jose, CA, United States). Normality and equal variance assumptions were tested with the Kolmogorov-Smirnov and the Levene median test, respectively. Group comparisons were performed with unpaired t-test, One-Way ANOVA followed by Holm-Sidak post hoc test, Two-Way Repeated Measures ANOVA followed by pairwise comparison using the Student–Newman–Keuls test or Fisher Exact Test with Freeman-Halton extension.
In case of non-normally distributed parameters, Kruskal-Wallis tests followed by pairwise comparisons with Dunn’s methods were applied. Data are presented as mean ± SEM or as boxplots, in which boxes represent lower quartile, median and upper quartile, and whiskers 1.5× interquartile range. p < 0.05 was considered statistically significant.
Results
Long-Chain Acylcarnitine Concentrations
First, we measured the acylcarnitine profiles in iCTRL-CMs and iVLCADD1-CMs cultured under standard culture conditions and in the presence of 400 µM carnitine for 48 h. Compared to iCTRL-CMs, LCAC concentrations were higher in iVLCADD1-CMs cultured in the absence of carnitine (Figure 1A). Culturing the iVLCADD1-CMs in the presence of carnitine resulted even in a further increase (Figure 1A). For example, C18:1-acylcarnitine increased from 38.8 ± 4.0 pmol/mg (three replicates) in iVLCADD1-CMs cultured in standard medium without additional carnitine to 138.5 ± 38.9 pmol/mg (three replicates) in iVLCADD1-CMs cultured with carnitine. In iVLCADD1-CMs, carnitine also resulted in a significant increase in LCACs levels in the culture medium with most prominent effects on C14:1-acylcarnitine (Figure 1B). This indicates that the hiPSC-CMs are able to excrete the accumulating LCACs.
FIGURE 1. Box plots of acylcarnitine profiles in cell pellets and culture medium. (A) Box plots of acylcarnitine profiles in cell pellets of iCTRL-CMs and iVLCADD1-CMs cultured in absence or in presence of 400 µM carnitine for 48 h. Asterisks indicate significant differences between groups (One-Way ANOVA followed by Holm-Sidak post hoc test or Kruskal-Wallis tests followed by pairwise comparisons with Dunn’s methods). Data from three differentiations obtained from cell pellets 30 days after differentiation. (B) Box plots of acylcarnitine profiles measured in culture medium of VLCADD1-CMs cultured under standard culture conditions or in presence of 400 µM carnitine for 48 h. Asterisks indicate significant differences between groups (Unpaired t-test; data from three differentiations measured 30 days after differentiation.
Carnitines Incubation Does Not Affect Action Potentials
Next, we assessed the AP parameters of the iCTRL-CMs and compared them to those of the iVLCADD1-CMs cultured in the absence or presence of carnitine. The AP parameters were analyzed as depicted in Figure 2A. Figure 2B shows a typical AP at 1 Hz stimulation of an iCTRL-CMs as well as typical 1 Hz stimulated APs recorded from an iVLCADD1-CMs cultured in the absence or presence of carnitine. Average AP parameters obtained from a total of 16 iCTRL-CMs, 11 iVLCADD1-CMs cultured in standard conditions, and 14 iVLCADD1-CMs cultured in the presence of carnitine are summarized in Figure 2C. The average MDP was around −80 mV and the Vmax was between 50 and 150 V/s, but both MDP and Vmax did not differ significantly between the three hiPSC-CMs groups. In all cell lines, APs overshot the 0 mV level, but the APA was significantly lower in the iVLCADD1-CMs in the absence as well as in the presence of carnitine compared to the iCTRL-CMs. iCTRL-CMs had an APA of 120 ± 2.7 mV (n = 16), while the iVLCADD1-CMs lines showed an APA of 99 ± 3.9 mV (n = 11) and 102 ± 3.4 mV (n = 14) in the absence and presence of carnitine, respectively. APs of both iVLCADD1-CMs groups repolarized earlier and faster than iCTRL-CMs (Figure 2A), resulting in a significantly shorter APD20, APD50, and APD90. For example, APD90 was 190 ± 9.7 mV (n = 16) in iCTRL-CMs and 127 ± 12.2 mV (n = 11) and 126 ± 21.0 mV (n = 14) in the iVLCADD1-CMs in absence and presence of carnitine. Neither APA nor the APDs differ significantly between iVLCADD1-CMs cultured in the presence of 400 µM carnitine and iVLCADD1-CMs cultured under standard conditions.
FIGURE 2. Action potential (AP) characteristics at 1 Hz of single iCTRL-CMs and iVLCADD1-CMs cultured in absence or presence of carnitine. (A) Illustration of the analyzed AP parameters, i.e., maximum diastolic potential (MDP), AP amplitude (APA), maximum upstroke velocity (Vmax), and AP duration at 20, 50, and 90% repolarization (APD20, APD50, and APD90, respectively). (B) Typical APs of iCTRL-CMs and iVLCADD1-CMs cultured in absence or presence of carnitine. (C) Box plots of the APs characteristics. Asterisks indicate significant differences between groups (One-Way ANOVA followed by Holm-Sidak post hoc test or Kruskal-Wallis tests followed by pairwise comparisons with Dunn’s methods). iCTRL-CM and iVLCADD1-CM data are from five to four differentiations, respectively, and were measured 46 and 47 days after differentiation.
Cardiac APs show a frequency dependence in morphology (Carmeliet, 2004) and the iVLCADD1-CM electrical phenotype as well as potential effects of carnitine may be frequency-related. For this reason, we studied the AP characteristics at 0.2–4 Hz stimulation. Figure 3A shows typical APs and Figure 3B summarizes the average AP parameters that were significantly different between iCTRL-CMs and iVLCADD1-CMs at 1 Hz. APA, APD20, APD50, and APD90 of iCTRL-CMs were significantly different from iVLCADD1-CMs cultured without and in the presence of carnitine at all frequencies, with exception of the APDs at 4 Hz in absence of carnitine (Figure 3B). As illustrated in Figure 3, the APs of iCTRL-CMs shows a clear frequency dependency with a decrease in APA and APDs upon both slower and faster pacing frequencies compared to the 1 Hz pacing frequency. This is depicted in more detail in Figure 4A, in which the iCTRL-CMs data presented in Figure 3B are replotted, but now with asterisks which indicate the significant differences between stimulus frequencies. In the VLCADD1 lines on the other hand, frequency dependence, especially those of the APDs, was less prominent as can be seen from the limited number of significant differences (Figures 4B,C). Thus, in iVLCADD1-CMs, APA and APDs are decreased at all frequencies compared to iCTRL-CMs, whereas the frequency dependence of AP parameters is hardly present in iVLCADD1-CMs. In addition, these data demonstrate that carnitine does not prevent these effects.
FIGURE 3. AP parameters at 0.2–4 Hz pacing of iCTRL-CMs and iVLCADD1-CMs cultured in absence or presence of carnitine. (A) Typical APs of iCTRL-CMs and iVLCADD1-CMs cultured in absence or presence of carnitine stimulated at 0.2, 0.5, 1, 2, 3, and 4 Hz. (B) APA, APD20, APD50, and APD90 at 0.2–4 Hz. Asterisks indicate significant differences between groups (Two-way Repeated Measures ANOVA). iCTRL-CM and iVLCADD1-CM data are from five to four differentiations, respectively, and were measured 46 and 47 days after differentiation.
FIGURE 4. (A–C) Frequency dependence of APA, APD20, APD50, and APD90 in iCTRL-CMs (A) and iVLCADD1-CMs cultured in absence (B) or presence of carnitine (C). Asterisks indicate significant differences between stimulus frequencies (Two-way Repeated Measures ANOVA). iCTRL-CM and iVLCADD1-CM data are from five to four differentiations, respectively, and were measured 46 and 47 days after differentiation.
Carnitines Incubation Does Not Affect the Occurrence of Afterdepolarizations
We also assessed the effects of carnitine on the occurrence of DADs and EADs, both important cellular electrophysiological triggers for cardiac arrhythmias (Hoffman and Rosen, 1981; Wit, 2018). Susceptibility to DADs was tested by applying a fast, 3-Hz pacing episode (10-s) followed by an 8-s pause (Figure 5B). After the pause, we evoked a single AP (Figure 5B) to test the inducibility of EADs. In addition, we counted the occurrence of EADs during continuous 0.2 Hz stimulation. The results of DADs and EADS are shown in Figures 5A,B, respectively.
FIGURE 5. Delayed and early afterdepolarization (DAD and EAD, respectively) inducibility of iCTRL-CMs and/or iVLCADD1-CMs cultured in absence or presence of carnitine. (A) Typical DAD examples of an iCTRL-CM and iVLCADD1-CMs cultured in absence or presence of carnitine. The arrows indicate the DADs. (B) Illustration of the protocol to induce DADs. (C) Average occurrence of DADs. (D) APs at 0.2 Hz pacing of an iVLCADD1-CM cultured in absence and an iVLCADD1-CM cultured in presence of carnitine. EADs were only present once every three to five APs and the APs exhibiting an EAD are indicated by dashed lines. (E) EADs occurrence at 0.2 Hz in iCTRL-CMs and iVLCADD1-CMs cultured in absence or presence of carnitine. Asterisks indicate significant differences between groups (Kruskal-Wallis test followed by pairwise comparisons with Dunn’s methods). iCTRL-CM and iVLCADD1-CM data are from five to four differentiations, respectively, and were measured 46 and 47 days after differentiation.
As illustrated in the typical examples of Figures 5A, fast pacing resulted in DADs (indicated with arrows) in the post-pacing pause in both iVLCADD1-CMs groups and not in the iCTRL-CMs. Figure 5C summarizes the occurrence of DADs for a total of 10 iCTRL-CMs, 11 iVLCADD1-CMs cultured in the absence, and 13 iVLCADD1-CMs cultured in the presence of carnitine. As already observed in the typical examples of Figure 5A, and confirming our previous findings (Knottnerus et al., 2020), there is a significant increase of DADs in the iVLCADD1-CMs lines (Figure 5C). However, the number of DADs was not significantly different between the iVLCADD1-CMs cultured in absence and presence of carnitine. EADs were never observed during the single AP evoked after the 8-s pause. In addition to the experiments presented in Figure 5A, we re-analyzed the experiments in which we paced the hiPSC-CMs continuously at a slow 0.2 Hz stimulus frequency for the occurrence of EADs. The data, obtained from a total of 12 iCTRL-CMs, 11 iVLCADD1-CMs cultured in standard conditions, and 13 iVLCADD1-CMs cultured in the presence of carnitine, further demonstrates that the presence of EADs was rare (Figures 5D,E). In iVLCADD1-CMs cultured without as well as in iVLCADD1-CMs cultured in the presence of carnitine, we only observed EADs in one cell (Figures 5D,E) and the EADs occurred once every three to five APs. No EADs were ever observed in the iCTRL-CMs. Thus, these experiments demonstrate a high incidence of DADs in iVLCADD1-CMs, but highlight the lack of protective effect of carnitine supplementation in afterdepolarization occurrence.
Discussion
In recent years, supplementation of L-carnitine to patients with a lcFAO disorder has remained a topic of debate (Spiekerkoetter et al., 2010; Pena et al., 2016; Yamada and Taketani, 2019). The main argument against L-carnitine supplementation is the suspicion that high levels of intrinsic and/or circulating LCACs may cause cardiac arrhythmias (Corr et al., 1989). Here, we studied the effects of L-carnitine supplementation in hiPSC-CMs generated from a VLCADD patient. We observed an increase in intracellular levels of LCACs, as well as an increase in excreted LCACs. There was neither an improvement in AP parameters nor in the occurrence of afterdepolarizations after L-carnitine supplementation. We conclude that LCACs accumulation and electrophysiological abnormalities in VLCADD hiPSC-CMs are not restored by carnitine supplementation, which suggests that carnitine treatment is not beneficial against cardiac arrhythmias in VLCADD patients.
iVLCADD1-CMs cultured under standard conditions exhibit a lower APA and shorter AP at a wide range of stimulus frequencies compared to iCTRL-CMs (Figures 2, 3). The exact mechanism is unknown, but is likely related to the LCAC accumulation (Figure 1A) and the consequent L-type Ca2+ current (ICa,L) decrease (Wu and Corr, 1992) and increased systolic Ca2+i concentrations resulting in enhanced Ca2+-induced ICa,L inactivation and increased K+ currents (Knottnerus et al., 2020). Indeed, reducing LCAC concentrations with either resveratrol or etomoxir restored APA and APDs in VLCADD hiPSC-CMs (Knottnerus et al., 2020). This pharmacological rescue of the metabolic and electrical phenotypes further indicates that these electrical abnormalities were not due to different genetic backgrounds in our used cell lines. The latter is also supported by experiments with a second control hiPSC line (for details, see Veerman et al., 2016), which show almost identical AP parameters (and the absence of DADs and EADs) as the presently used iCTRL-CMs (data not shown). While the iCTRL-CMs showed a clear frequency dependency of APA and APDs, this was less pronounced in the standard cultured iVLCADD1-CMs (Figure 4). The APA and APD decrease in iCTRL-CMs at faster stimulus frequencies is likely due to a combination of incomplete recovery of ICa,L and accumulation of repolarizing delayed rectifier K+ currents (Carmeliet, 2004). The AP shortening at slow stimulus frequencies in iCTRL-CMs, on the other hand, is likely due the repolarizing transient outward K+ current, Ito1. In hiPSC-CMs, Ito1 has an extremely slow recovery from inactivation (Cordeiro et al., 2013) and consequently plays only a role in AP repolarizing at low stimulus frequencies (Ma et al., 2018). The loss of frequency dependency of APA and APDs observed in iVLCADD1-CMs may be related to the already above mentioned decrease in ICa,L function and the obvious fast AP repolarization in iVLCADD1-CMs, which can mask subtle changes in net currents (Gaur et al., 2020). We found a high incidence of DADs in iVLCADD1-CMs cultured under standard conditions (Figures 5A,C), consistent with our previous study (Knottnerus et al., 2020). In human CMs, DADs are due to Ca2+i overload-induced spontaneous Ca2+ release from the sarcoplasmic reticulum (SR), which activates an inward current carried by the Na+-Ca2+ exchanger (Verkerk et al., 2001). Ca2+i overload was previously observed in hiPSC-CMs from VLCADD patients (Knottnerus et al., 2020). EADs in iCTRL-CMs were absent and rare in iVLCADD-CMs (Figures 5D,E). EADs can be divided in the so-called phase 2 and late-phase 3 EADs which are supposed to be due to different mechanisms (Qu et al., 2013). Phase 2 EADs are caused by re-activation of ICa,L which may occur during prolonged APs, while late-phase 3 EADs are likely due to the Na+-Ca2+ exchange current when the AP is substantially shortened such that the Ca2+ transient outlast the AP repolarization (Qu et al., 2013). It is tempting to define the observed EADs as late-phase 3 EADs based on the observed AP shortening (Figure 2) and previously observed increased systolic Ca2+i (Knottnerus et al., 2020) which may result in a larger Na+-Ca2+ exchange current, but further studies are required to address this issue in detail.
The electrophysiological properties of the iVLCADD1-CMs cultured for 48 h in the presence of carnitine supplementation were not different from iVLCADD1-CMs cultured under standard conditions (Figures 2–5), demonstrating that the lack of effect of carnitine in restoring electrophysiological abnormalities. Supplementation of carnitine for 48 h was sufficient to increase the LCAC concentrations (Figure 1A), indicating that the absence of electrophysiological effects is not related to the length of the exposure time to the drugs in this particular hiPSC-CM model. In addition, we have previously observed that 48 h of drug treatment is sufficient to affect mitochondrial function in our iVLCADD1-CMs model, therefore indicating that this hiPSC-CMs model is useful for metabolic disorder drug studies (Knottnerus et al., 2020), even though hiPSC-CMs may have a less developed mitochondrial system for FAO flux compared to adult cardiomyocytes due to their immaturity (Chen et al., 2016). While the carnitine treatment increased the LCAC concentrations, electrophysiological abnormalities were not aggravated. This suggests that there is not a simple linear relationship between AP shortening and incidence of DADs and the LCACs. In addition, it indicates that carnitine treatment can be used for other, non-cardiac VLCADD-induced dysfunctions without increasing the risk for cardiac arrhythmias.
Secondary carnitine deficiency may occur in patients with VLCADD as a result of acylcarnitine formation from high levels of non-degraded acyl-CoAs. The pathophysiological relevance of carnitine depletion for the signs and symptoms observed in this disorder is not clear. Still, carnitine supplementation is often part of the treatment of VLCADD, especially when low (plasma) carnitine concentrations are found in patients (Pena et al., 2016). Evaluation of the efficacy of this intervention is very complex, due to the low numbers of patients and heterogeneous phenotypic presentations. In addition, acylcarnitine profiles measured in plasma do not correlate well with profiles measured in tissue (muscle tissue or white adipose tissue) (Schooneman et al., 2014). We found in iVLCADD1-CMs an increased LCAC accumulation as well as an increased LCAC secretion after carnitine supplementation. So far, there is limited knowledge on the cellular export mechanisms of acylcarnitines (Kim et al., 2017), but the increased LCAC accumulation induced by carnitine supplementation is in agreement with findings in mouse models for VLCADD (Liebig et al., 2006; Primassin et al., 2008; Bakermans et al., 2013). Bakermans and colleagues (Bakermans et al., 2013), however, only observed an increase of C16 and C18:1-carnitine whereas C14:1 levels were not affected in bloodspots after carnitine supplementation.
The heart lacks enzymes for carnitine synthesis and is therefore dependent on carnitine import from the circulation via the sodium dependent organic cation transporter (OCTN2). In a recent study, it was shown that haploinsufficiency of OCTN2 in lcFAO disorder mice results in decreased free carnitine levels and subsequently decrease in tissue and plasma LCAC accumulation but does not significantly affect clinically relevant outcome parameters (hypoglycemia, heart weight or liver weight) (Ranea-Robles et al., 2020). However, it should be noted that the C14:1-acylcarnitine levels were still higher in lcFAO disorder with OCTN2 haploinsufficiency when compared to wild type in heart, liver and plasma. Therefore, toxicity of acylcarnitines cannot be excluded.
In the present study, we have tested the effects of carnitine treatment on single VLCADD hiPSC-CMs and not on iCTRL-CMs. Whereas this approach allows a detailed study of the effects of a well-controlled dose of carnitine on the intrinsic LCAC concentrations, AP configuration and afterdepolarizations of the VLCADD hiPSC-CMs, it does not take into account the metabolic and structural complexity and heterogeneity of the intact myocardium or body. Therefore, caution is warranted when translating our in vitro results to the in vivo situation. Furthermore, although our VLCADD model has a human background thereby avoiding potential species differences, our study was limited to one patient-specific hiPSC line with ACADVL gene mutation. This provided us the opportunity to investigate the effects of carnitine on cardiomyocytes derived from this patient-specific hiPSC line in detail, but further studies are needed to test the carnitine effects on other VLCADD models. In the present study, we focused on LCAC concentrations and electrophysiology in hiPSC-CMs but VLCADD has a broad clinical spectrum (Ribas and Vargas, 2020; Wanders et al., 2020). In future studies, our VLCADD hiPSC line may also be differentiated to a liver and skeletal muscle cell fate (van der Wall et al., 2018; Corbett and Duncan, 2019) to assesses the effects of VLCADD and potential treatment in these cell types.
In conclusion, carnitine enhanced LCAC concentrations in hiPSC-CMs with ACADVL gene mutations and did not improve the electrophysiological abnormalities. This indicates that carnitine treatment is unlikely to be beneficial against cardiac arrhythmias in patients with ACADVL gene mutations.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author Contributions
SK, CB, and RH conceived and designed the study. KG generated the hiPSC lines. IM. maintained the hiPSC lines and differentiated them to cardiomyocytes. IM, SK, SF and JB dissociated cells and performed biochemical experiments. AV and VP performed electrophysiology measurements. SK and AV analyzed and interpreted the data and drafted the manuscript. LIJ, GV, RW, and FW contributed to scientific discussions. All authors critically revised the manuscript and gave final approval of the version to be published.
Funding
KG was supported by the Free State of Saxony and the European Union EFRE (SAB project “PhänoKard”) and by the DFG (GU595/3-1, IRTG2251). CB was supported by the Dutch Heart Foundation (CVON PREDICT2 project), Netherlands Organization for Scientific Research (VICI fellowship, 016.150.610) and Fondation Leducq. RH was supported by a VIDI grant from ZonMw (no. 91715305) and a grant from the Velux Stiftung (no. 1063).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Abbreviations
AP, action potential; APA, maximal action potential amplitude; APD20, action potential duration at 20% repolarization; APD50, action potential duration at 50% repolarization; APD90, action potential duration at 90% repolarization; Ca2+i, intracellular Ca2+; Cm, cell membrane capacitance; DAD, delayed afterdepolarization; EAD, early afterdepolarization; hiPSC-CMs, human induced pluripotent stem cell derived cardiomyocytes; hiPSCs, human induced pluripotent stem cells; ICa,L, L-type Ca2+ current; IK1, inward rectifying K+ current; Ito1, transient outward K+ current; KO, knock out; LCAC, long-chain acylcarnitines; LCAD, long-chain acyl-CoA dehydrogenase; MDP, maximum diastolic potential; VLCAD, very long-chain acyl-CoA dehydrogenase; VLCADD, VLCAD deficiency; Vmax, maximum upstroke velocity.
References
Bakermans, A. J., van Weeghel, M., Denis, S., Nicolay, K., Prompers, J. J., and Houten, S. M. (2013). Carnitine supplementation attenuates myocardial lipid accumulation in long-chain acyl-CoA dehydrogenase knockout mice. J. Inherit. Metab. Dis. 36 (6), 973–981. doi:10.1007/s10545-013-9604-4
Barry, P. H., and Lynch, J. W. (1991). Liquid junction potentials and small cell effects in patch-clamp analysis. J. Membr. Biol. 121, 101–117. doi:10.1007/BF01870526
Bonnet, D., Martin, D., de Lonlay, P., Villain, E., Jouvet, P., Rabier, D., et al. (1999). Arrhythmias and conduction defects as presenting symptoms of fatty acid oxidation disorders in children. Circulation. 100 (22), 2248–2253. doi:10.1161/01.cir.100.22.2248
Carmeliet, E. (2004). Intracellular Ca2+ concentration and rate adaptation of the cardiac action potential. Cell. Calcium. 35 (6), 557–573. doi:10.1016/j.ceca.2004.01.010
Chegary, M., te Brinke, H., Ruiter, J. P., Wijburg, F. A., Stoll, M. S., Minkler, P. E., et al. (2009). Mitochondrial long chain fatty acid beta-oxidation in man and mouse. Biochim. Biophys. Acta. 1791 (8), 806–815. doi:10.1016/j.bbalip.2009.05.006
Corbett, J. L., and Duncan, S. A. (2019). IPSC-derived hepatocytes as a platform for disease modeling and drug discovery. Front. Med. (Lausanne). 6, 265. doi:10.3389/fmed.2019.00265
Cordeiro, J. M., Nesterenko, V. V., Sicouri, S., Goodrow, Jr, R. J., Treat, J. A., Desai, J. A., et al. (2013). Identification and characterization of a transient outward K+ current in human induced pluripotent stem cell-derived cardiomyocytes. J. Mol. Cell Cardiol. 60: 36–46. doi:10.1016/j.yjmcc.2013.03.014
Corr, P. B., Creer, M. H., Yamada, K. A., Saffitz, J. E., and Sobel, B. E. (1989). Prophylaxis of early ventricular fibrillation by inhibition of acylcarnitine accumulation. J. Clin. Invest. 83 (3), 927–936. doi:10.1172/JCI113978
Cranefield, P. F. (1977). Action potentials, afterpotentials, and arrhythmias. Circ. Res. 41 (4), 415–423. doi:10.1161/01.res.41.4.415
Dudek, J., Cheng, I.-F., Balleininger, M., Vaz, F. M., Streckfuss-Bömeke, K., Hübscher, D., et al. (2013). Cardiolipin deficiency affects respiratory chain function and organization in an induced pluripotent stem cell model of Barth syndrome. Stem Cell Res. 11 (2), 806–819. doi:10.1016/j.scr.2013.05.005
Gaur, N., Ortega, F., Verkerk, A. O., Mengarelli, I., Krogh-Madsen, T., Christini, D. J., et al. (2020). Validation of quantitative measure of repolarization reserve as a novel marker of drug induced proarrhythmia. J. Mol. Cell. Cardiol. 145, 122–132. doi:10.1016/j.yjmcc.2020.04.019
Hoekstra, M., Mummery, C. L., Wilde, A. A. M., Bezzina, C. R., and Verkerk, A. O. (2012). Induced pluripotent stem cell derived cardiomyocytes as models for cardiac arrhythmias. Front. Physiol. 3, 346. doi:10.3389/fphys.2012.00346
Hoffman, B. F., and Rosen, M. R. (1981). Cellular mechanisms for cardiac arrhythmias. Circ. Res. 49 (1), 1–15. doi:10.1161/01.res.49.1.1
Houten, I. Y., Matsa, E., and Wu, J. C. (2016). Induced pluripotent stem cells: at the heart of cardiovascular precision medicine. Nat. Rev. Cardiol. 13 (6), 333–349. doi:10.1038/nrcardio.2016.36
Houten, S. M., and Wanders, R. J. A. (2010). A general introduction to the biochemistry of mitochondrial fatty acid β-oxidation. J. Inherit. Metab. Dis. 33 (5), 469–477. doi:10.1007/s10545-010-9061-2
Houtkooper, R. H., Pirinen, E., and Auwerx, J. (2012). Sirtuins as regulators of metabolism and healthspan. Nat. Rev. Mol. Cell Biol. 13 (4), 225–238. doi:10.1038/nrm3293
Kim, H. I., Raffler, J., Lu, W., Lee, J. J., Abbey, D., Saleheen, D., et al. (2017). Fine mapping and functional analysis reveal a role of SLC22A1 in acylcarnitine transport. Am. J. Hum. Genet. 101 (4), 489–502. doi:10.1016/j.ajhg.2017.08.008
Knottnerus, S. J. G., Bleeker, J. C., Wüst, R. C. I., Ferdinandusse, S., IJlst, L., Wijburg, F. A., et al. (2018). Disorders of mitochondrial long-chain fatty acid oxidation and the carnitine shuttle. Rev. Endocr. Metab. Disord. 19 (1), 93–106. doi:10.1007/s11154-018-9448-1
Knottnerus, S. J. G., Mengarelli, I., Wüst, R. C. I., Baartscheer, A., Bleeker, J. C., Coronel, R., et al. (2020). Electrophysiological abnormalities in VLCAD deficient hiPSC-cardiomyocytes can be improved by lowering accumulation of fatty acid oxidation intermediates. Int. J. Mol. Sci. 21 (7), 2589. doi:10.3390/ijms21072589
Liebig, M., Gyenes, M., Brauers, G., Ruiter, J. P., Wendel, U., Mayatepek, E., et al. (2006). Carnitine supplementation induces long-chain acylcarnitine production: studies in the VLCAD-deficient mouse. J. Inherit. Metab. Dis. 29 (2-3), 343–344. doi:10.1007/s10545-006-0249-4
Longo, N., di San Filippo, C. A., and Pasquali, M. (2006). Disorders of carnitine transport and the carnitine cycle. Am. J. Med. Genet. C. 142 (2), 77–85. doi:10.1002/ajmg.c.30087
Lopaschuk, G. D., Ussher, J. R., Folmes, C. D., Jaswal, J. S., and Stanley, W. C. (2010). Myocardial fatty acid metabolism in health and disease. Physiol. Rev. 90 (1), 207–258. doi:10.1152/physrev.00015.2009
Ma, D., Liu, Z., Loh, L. J., Zhao, Y., Li, G., Liew, R., et al. (2018). Identification of an INa-dependent and Ito-mediated proarrhythmic mechanism in cardiomyocytes derived from pluripotent stem cells of a Brugada syndrome patient. Sci. Rep. 8 (1), 11246. doi:10.1038/s41598-018-29574-5
Meijer van Putten, R. M. E., Mengarelli, I., Guan, K., Zegers, J. G., van Ginneken, A. C. G., Verkerk, A. O., et al. (2015). Ion channelopathies in human induced pluripotent stem cell derived cardiomyocytes: a dynamic clamp study with virtual Ik1. Front. Physiol. 6, 7. doi:10.3389/fphys.2015.00007
Mordwinkin, N. M., Burridge, P. W., and Wu, J. C. (2013). A review of human pluripotent stem cell-derived cardiomyocytes for high-throughput drug discovery, cardiotoxicity screening, and publication standards. J. Cardiovasc. Transl. Res. 6 (1), 22–30. doi:10.1007/s12265-012-9423-2
Nasser, M., Javaheri, H., Fedorowicz, Z., and Noorani, Z. (2012). Carnitine supplementation for inborn errors of metabolism. Cochrane Database Syst. Rev. 15 (2), Cd006659. doi:10.1002/14651858.CD006659.pub3
Nerbonne, J. M., and Kass, R. S. (2005). Molecular physiology of cardiac repolarization. Physiol. Rev. 85, 1205–1253. doi:10.1152/physrev.00002.2005
Pena, L. D., van Calcar, S. C., Hansen, J., Edick, M. J., Walsh-Vockley, C., Leslie, N., et al. (2016). Outcomes and genotype-phenotype correlations in 52 individuals with VLCAD deficiency diagnosed by NBS and enrolled in the IBEM-IS database. Mol. Genet. Metab. 118 (4), 272–281. doi:10.1016/j.ymgme.2016.05.007
Primassin, S., Ter Veld, F., Mayatepek, E., and Spiekerkoetter, U. (2008). Carnitine supplementation induces acylcarnitine production in tissues of very long-chain acyl-CoA dehydrogenase-deficient mice, without replenishing low free carnitine. Pediatr. Res. 63 (6), 632–637. doi:10.1203/PDR.0b013e31816ff6f0
Qu, Z., Xie, L. H., Olcese, R., Karagueuzian, H. S., Chen, P. S., Garfinkel, A., et al. (2013). Early afterdepolarizations in cardiac myocytes: beyond reduced repolarization reserve. Cardiovasc. Res. 99 (1), 6–15. doi:10.1093/cvr/cvt104
Ranea-Robles, P., Yu, C., van Vlies, N., Vaz, F. M., and Houten, S. M. (2020). Slc22a5 haploinsufficiency does not aggravate the phenotype of the long-chain acyl-CoA dehydrogenase KO mouse. J. Inherit. Metab. Dis. 43 (3), 486–495. doi:10.1002/jimd.12204
Reuter, S. E., and Evans, A. M. (2012). Carnitine and acylcarnitines: pharmacokinetic, pharmacological and clinical aspects. Clin. Pharmacokinet. 51 (9), 553–572. doi:10.1007/BF03261931
Ribas, G. S., and Vargas, C. R. (2020). Evidence that oxidative disbalance and mitochondrial dysfunction are involved in the pathophysiology of fatty acid oxidation disorders. Cell. Mol. Neurobiol. [in Press]. doi:10.1007/s10571-020-00955-7
Roussel, J., Labarthe, F., Thireau, J., Ferro, F., Farah, C., Roy, J., et al. (2016). Carnitine deficiency induces a short QT syndrome. Heart Rhythm. 13 (1), 165–174. doi:10.1016/j.hrthm.2015.07.027
Schooneman, M. G., Achterkamp, N., Argmann, C. A., Soeters, M. R., and Houten, S. M. (2014). Plasma acylcarnitines inadequately reflect tissue acylcarnitine metabolism. Biochim. Biophys. Acta. 1841 (7), 987–994. doi:10.1016/j.bbalip.2014.04.001
Spiekerkoetter, U., Bastin, J., Gillingham, M., Morris, A., Wijburg, F., and Wilcken, B. (2010). Current issues regarding treatment of mitochondrial fatty acid oxidation disorders. J Inherit. Metab. Dis. 33 (5), 555–561. doi:10.1007/s10545-010-9188-1
Spiekerkoetter, U., Lindner, M., Santer, R., Grotzke, M., Baumgartner, M. R., Boehles, H., et al. (2009). Treatment recommendations in long-chain fatty acid oxidation defects: consensus from a workshop. J. Inherit. Metab. Dis. 32 (4), 498–505. doi:10.1007/s10545-009-1126-8
Tohyama, S., Hattori, F., Sano, M., Hishiki, T., Nagahata, Y., Matsuura, T., et al. (2013). Distinct metabolic flow enables large-scale purification of mouse and human pluripotent stem cell-derived cardiomyocytes. Cell. Stem. Cell. 12, 127–137. doi:10.1016/j.stem.2012.09.013
Treem, W. R., Stanley, C. A., Hale, D. E., Leopold, H. B., and Hyams, J. S. (1991). Hypoglycemia, hypotonia, and cardiomyopathy: the evolving clinical picture of long-chain acyl-CoA dehydrogenase deficiency. Pediatrics. 87 (3), 328–333.
van der Wall, E., Herrero-Hernandez, P., Wan, R., Broeders, M., in’t Groen, S. L. M., van Gestel, T. J. M., et al. (2018). Large-scale expansion of human iPSC-derived skeletal muscle cells for disease modeling and cell-based therapeutic strategies. Stem Cell Rep. 10, 1975–1990. doi:10.1016/j.stemcr.2018.04.002
Vaz, F. M., and Wanders, R. J. A. (2002). Carnitine biosynthesis in mammals. Biochem J. 361 (Pt 3), 417–429. doi:10.1042/0264-6021:3610417
Veerman, C. C., Kosmidis, G., Mummery, C. L., Casini, S., Verkerk, A. O., and Bellin, M. (2015). Immaturity of human stem-cell-derived cardiomyocytes in culture: fatal flaw or soluble problem? Stem Cells Dev. 24 (9), 1035–1052. doi:10.1089/scd.2014.0533
Veerman, C. C., Mengarelli, I., Guan, K., Stauske, M., Barc, J., Tan, H. L., et al. (2016). hiPSC-derived cardiomyocytes from Brugada Syndrome patients without identified mutations do not exhibit clear cellular electrophysiological abnormalities. Sci. Rep. 6, 30967. doi:10.1038/srep30967
Verkerk, A. O., Veldkamp, M. W., Baartscheer, A., Schumacher, C. A., Klöpping, C., van Ginneken, A. C. G., et al. (2001). Ionic mechanism of delayed afterdepolarizations in ventricular cells isolated from human end-stage failing hearts. Circulation. 104 (22), 2728–2733. doi:10.1161/hc4701.099577
Wanders, R. J. A., Visser, G., Ferdinandusse, S., Vaz, F. M., and Houtkooper, R. H. (2020). Mitochondrial fatty acid oxidation disorders: laboratory diagnosis, pathogenesis, and the complicated route to treatment. J. Lipid Atheroscler. 9 (3), 313–333. doi:10.12997/jla.2020.9.3.313
Winter, S. C. (2003). Treatment of carnitine deficiency. J. Inherit. Metab. Dis. 26 (2), 171–180. doi:10.1023/a:1024433100257
Wit, A. L. (2018). Afterdepolarizations and triggered activity as a mechanism for clinical arrhythmias. Pacing Clin. Electrophysiol. 41 (8), 883–896. doi:10.1111/pace.13419
Wu, J., and Corr, P. B. (1992). Influence of long-chain acylcarnitines on voltage-dependent calcium current in adult ventricular myocytes. Am. J. Physiol. 263 (2), H410–H417. doi:10.1152/ajpheart.1992.263.2.H410
Keywords: very long-chain acyl-CoA dehydrogenase, arrhythmia < cardiovascular, acylcarnitines, action potential, human induced pluripotent stem cell derived cardiomyocytes, patients, treatment, carnitine
Citation: Verkerk AO, Knottnerus SJG, Portero V, Bleeker JC, Ferdinandusse S, Guan K, IJlst L, Visser G, Wanders RJA, Wijburg FA, Bezzina CR, Mengarelli I and Houtkooper RH (2021) Electrophysiological Abnormalities in VLCAD Deficient hiPSC-Cardiomyocytes Do not Improve with Carnitine Supplementation. Front. Pharmacol. 11:616834. doi: 10.3389/fphar.2020.616834
Received: 13 October 2020; Accepted: 24 November 2020;
Published: 12 January 2021.
Edited by:
Tamer M. A. Mohamed, University of Louisville, United StatesReviewed by:
Jerry Vockley, University of Pittsburgh, United StatesLukas Cyganek, University Medical Center Göttingen, Germany
Eric Stephen Goetzman, University of Pittsburgh, United States
Copyright © 2021 Verkerk, Knottnerus, Portero, Bleeker, Ferdinandusse, Guan, IJlst, Visser, Wanders, Wijburg, Bezzina, Mengarelli and Houtkooper. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Connie R. Bezzina, c.r.bezzina@amsterdamumc.nl; Riekelt H. Houtkooper, r.h.houtkooper@amsterdamumc.nl
†These authors have contributed equally to this work.