- 1Department of Anesthesiology and Perioperative Medicine, University of Pittsburgh, Pittsburgh, PA, United States
- 2Division of Cardiology, Department of Medicine, University of Pittsburgh Medical Center, Pittsburgh, PA, United States
- 3Department of Anesthesiology and Perioperative Medicine, University of Pittsburgh Medical Center, Pittsburgh, PA, United States
- 4Department of Physiology, University of Oklahoma Health Sciences Center, Oklahoma City, OK, United States
- 5Cardiac Arrhythmia Center and Neurocardiology Research Program of Excellence, David Geffen School of Medicine at UCLA, Los Angeles, CA, United States
Introduction: Myocardial ischemia disrupts the cardio-spinal neural network that controls the cardiac sympathetic preganglionic neurons, leading to sympathoexcitation and ventricular tachyarrhythmias (VTs). Spinal cord stimulation (SCS) is capable of suppressing the sympathoexcitation caused by myocardial ischemia. However, how SCS modulates the spinal neural network is not fully known.
Methods: In this pre-clinical study, we investigated the impact of SCS on the spinal neural network in mitigating myocardial ischemia-induced sympathoexcitation and arrhythmogenicity. Ten Yorkshire pigs with left circumflex coronary artery (LCX) occlusion-induced chronic myocardial infarction (MI) were anesthetized and underwent laminectomy and a sternotomy at 4−5 weeks post-MI. The activation recovery interval (ARI) and dispersion of repolarization (DOR) were analyzed to evaluate the extent of sympathoexcitation and arrhythmogenicity during the left anterior descending coronary artery (LAD) ischemia. Extracellular in vivo and in situ spinal dorsal horn (DH) and intermediolateral column (IML) neural recordings were performed using a multichannel microelectrode array inserted at the T2-T3 segment of the spinal cord. SCS was performed for 30 min at 1 kHz, 0.03 ms, 90% motor threshold. LAD ischemia was induced pre- and 1 min post-SCS to investigate how SCS modulates spinal neural network processing of myocardial ischemia. DH and IML neural interactions, including neuronal synchrony as well as cardiac sympathoexcitation and arrhythmogenicity markers were evaluated during myocardial ischemia pre- vs. post-SCS.
Results: ARI shortening in the ischemic region and global DOR augmentation due to LAD ischemia was mitigated by SCS. Neural firing response of ischemia-sensitive neurons during LAD ischemia and reperfusion was blunted by SCS. Further, SCS showed a similar effect in suppressing the firing response of IML and DH neurons during LAD ischemia. SCS exhibited a similar suppressive impact on the mechanical, nociceptive and multimodal ischemia sensitive neurons. The LAD ischemia and reperfusion-induced augmentation in neuronal synchrony between DH-DH and DH-IML pairs of neurons were mitigated by the SCS.
Discussion: These results suggest that SCS is decreasing the sympathoexcitation and arrhythmogenicity by suppressing the interactions between the spinal DH and IML neurons and activity of IML preganglionic sympathetic neurons.
1. Introduction
The autonomic nervous system plays a critical role in regulating cardiac function (Armour, 2008). Myocardial ischemia by triggering the cardiac sensory neurons initiates the disruption in the cardiac autonomic nervous system (Zipes, 1990; Salavatian et al., 2019). This persistent afferent signaling by cardiac nociceptive afferents during myocardial ischemia causes an autonomic imbalance, manifested by sympathoexcitation, which leads to cardiac dysfunction and ventricular arrhythmias associated with sudden cardiac death (Zipes, 1990; Rubart and Zipes, 2005; Luqman et al., 2007; Fukuda et al., 2015; Rajendran et al., 2016).
During myocardial ischemic injury, triggered ischemia-sensitive cardiac afferent neurons transmit excitatory information to the dorsal horn of the spinal cord activating the spinal neural network (Dale et al., 2020; Omura et al., 2021). Activation of the ischemia-sensitive dorsal horn neurons leads to increased firing of the connected IML preganglionic sympathetic neuronal population. The excitatory signals from IML neurons to post-ganglionic neurons in the stellate and middle cervical ganglia results in increased sympathetic activity in the heart. This efferent sympathetic output to the heart can be exaggerated with myocardial injury, resulting in malignant ventricular tachyarrhythmias and potential for sudden cardiac death (Bourke et al., 2010; Fukuda et al., 2015). In chronic myocardial infarction (MI), significant neural remodeling at different levels of the cardiac neuraxis alters the global pattern of functional sympathetic innervation and increases cardiac arrhythmogenicity (Armour, 2004; Shen and Zipes, 2014; Vaseghi et al., 2014; Rajendran et al., 2016; Salavatian et al., 2022).
Cardiac neuromodulation therapies engage the underlying afferent and efferent neuronal projections for integrated reflex control of the heart. Neuromodulation via thoracic spinal cord stimulation (SCS) mitigates intrathoracic (intrinsic cardiac and intrathoracic extracardiac) reflexes while helping to preserve cardiac function (Odenstedt et al., 2011; Ardell, 2016; Howard-Quijano et al., 2017). SCS is also known to modify the capacity of cardiac sensory neurons to transduce the ischemic myocardium to the spinal cord dorsal horn neurons, thereby modulating spinal sympathetic efferent neuronal inputs to the heart (Ardell et al., 2019; Salavatian et al., 2019).
We have previously shown that SCS decreases myocardial ischemia-induced cardiac sympathoexcitation and ventricular arrhythmias (Howard-Quijano et al., 2017). However, it is not fully known how the SCS provides a therapeutic effect in the setting of myocardial infarction and myocardial ischemia to reduce cardiac sympathoexcitation and ventricular arrhythmias.
The primary purpose of this study was to determine how SCS modulates processing in the DH and IML spinal cord neuronal network during myocardial ischemia to mitigate cardiac arrhythmogenicity. We used a clinically relevant porcine model to mimic patients with prior MI who experience episodes of ventricular arrhythmias during the new onset of coronary artery ischemia. We hypothesized that SCS decreases myocardial ischemia-induced sympathoexcitation by suppressing intermediolateral nucleus hyperactivity and spinal neural synchrony.
2. Materials and methods
The Institutional Animal Care and Use Committee (IACUC) at the University of Pittsburgh approved the study protocol. All experiments were performed in compliance with the National Institution of Health Guide for the Care and Use of Laboratory Animals.
2.1. Left circumflex coronary artery MI creation
Ten Yorkshire pigs (four males and six females, mean ± SD, 31.2 ± 2.4 kg) were used. All animals were allowed to get acclimated for at least 5 days before survival surgery in the animal facility. Animals were sedated (tiletamine-zolazepam, 4 mg/kg IM + Xylazine, 2 mg/kg IM), intubated, and ventilated under general anesthesia using 1−3% of isoflurane with hemodynamic monitoring. Left circumflex coronary artery (LCX) MI was created as previously described (Salavatian et al., 2022). Under fluoroscopic guidance, an 8F AL2 guide catheter (Boston Scientific, Marlborough, MA, USA) was advanced from the femoral artery to the ostium of the left main coronary artery. A 3 mm percutaneous transluminal angioplasty balloon catheter (Abbott, Chicago, IL, USA) was advanced over a balance middleweight universal guidewire (Abbott, Chicago, IL, USA) into the LCX and positioned at the distal LCX. The balloon was inflated, and 3 mL of polystyrene microspheres (Polybead 90.0 μm, Polysciences, Warrington, PA, USA) followed by 1−2 mL of normal saline were injected through the lumen of the catheter. The balloon was then deflated and removed. MI was confirmed by the presence of ST-elevation or T-wave inversion on electrocardiogram (ECG) and coronary angiography showing a lack of flow in the artery (Figure 1). After the procedure, ECG and blood pressure were monitored until the pig remained hemodynamically stable for at least 20 min before extubation. Immediate external defibrillation was performed if the animal developed sustained ventricular tachycardia or fibrillation. After extubating, all animals were further monitored until they could ambulate.
Figure 1. Experimental protocol and methodology. (A) Chronic myocardial infarction was induced by injection of the microsphere in the left circumflex coronary artery (LCX). The scar zone was measured at the terminal experiments 4–5 weeks later. (B) Following sternotomy, acute myocardial ischemia was performed by ligating the left descending anterior coronary artery (LAD) and was confirmed by the ST-elevation on the epicardial electrograms (EGM). RV/LV: right/left ventricle. (C) 56-electrodes epicardial sock was used to record the electrogram and measure the activation recovery interval and dispersion of the repolarization. (D) A representative image of the location of the neural probe and the spinal cord stimulation (SCS). (E) A Multichannel multi-shank microelectrode array was used to record the activity of spinal cord neurons in the dorsal horn (DH) and intermediolateral column (IML) regions. (F) Spike sorting was performed to identify and classify the spiking activity of each neuron. Activities from 3 neurons (purple, red, pink) and one neuron (green) were identified based on their action potential shapes on channel 1 (CH1) and channel 2 (CH2), respectively. Parts of this figure were created with BioRender.com.
2.2. Terminal experiment – animal preparation
After 4−5 weeks of MI, the terminal experiment was performed and animal preparations were conducted as previously reported (Salavatian et al., 2019; Omura et al., 2021). Animals (n = 10, 50.7 ± 6.5 kg) were sedated with Telazol (4 mg/kg, intramuscular) and Xylazine (2 mg/kg, intramuscular), and mechanically ventilated by tracheal intubation. General anesthesia was induced and maintained with inhaled 2−5% of isoflurane during surgical preparation. Heart rate and ECG were monitored throughout the experiment using a Prucka CardioLab system (GE Healthcare, Fairfield, CT, USA). The carotid artery was catheterized for blood pressure monitoring. In addition, jugular veins were cannulated for intravenous saline infusion (10 ml/kg) and drug administration. Arterial blood gas was checked every hour with adjustment of ventilation if needed to maintain acid-base equilibrium. Body temperature was maintained by external warmers. Dorsal spinal laminectomy was done in a prone position to expose the spinal cord, then median sternotomy was then performed supine to expose the heart. General anesthesia was transitioned to intravenous alpha-chloralose (50 mg/kg initial bolus followed by a 20 mg/kg/h continuous infusion), which limits the impact on the autonomic nervous system and cardiac electrophysiology (Luckl et al., 2008; Omura et al., 2021). The experimental protocol was conducted in the left lateral position to record spinal cord neural activity and cardiac electrophysiology with acute myocardial ischemia. The depth of anesthesia was assessed throughout the experiments by monitoring corneal reflexes, jaw tone, and hemodynamic indices. All animals were euthanized by inducing ventricular fibrillation via potassium chloride administration under deep anesthesia. The heart was excised and chronic infarct size was measured with triphenyl tetrazolium chloride staining as previously reported (Takagawa et al., 2007).
2.3. Recording of spinal cord neuronal activity
Laminectomy was performed at T1-T4 levels of the spinal cord to expose the T2-T3 level and provide enough space for the insertion of the neural probe. The animal was placed in the right lateral position and following a small opening in the dura, a high-density two-dimensional penetrating microarray (64 electrode recording sites; NeuroNexus, Ann Arbor, MI, USA) probe, mounted on a micromanipulator, was inserted in the spinal cord at the T2-T3 level (Figure 1). The neural probe was inserted into the spinal cord to record the neuronal activity of the DH and IML neurons (and not the axons due to the high impedance of the headstage and filtering) in vivo and in situ. Recorded signals were amplified, digitized, and filtered (300 Hz – 3 kHz) using the SmartBox acquisition system (NeuroNexus, Ann Arbor, MI, USA).
2.4. Acute myocardial ischemia
We created acute myocardial ischemia as previously described (Howard-Quijano et al., 2021). Briefly, a Prolene suture was placed around the left anterior descending coronary artery (LAD) below the second diagonal branch of the LAD. The suture was led through a short polyethylene tubing segment, which was then used to ligate the coronary artery to induce cardiac ischemia. LAD ischemia was confirmed with the existence of ST elevations (Figure 1). Short LAD ischemia (3 min) was induced before and 1 min after SCS.
2.5. Cardiac mechanical and nociceptive stimuli
To identify the mechanosensitive and nociceptive neurons, a set of mechanical and nociceptive stimuli were performed. For the mechanical stimuli, the right and left anterior ventricular epicardium were touched for 15 s. To identify nociceptive neurons, ventricular (anterior) epicardial nociceptive stimuli were applied for 1 min using a gauze soaked in bradykinin (10 μM) or capsaicin (1 μM). The neurons that responded to both mechanical and nociceptive stimuli, were identified as multimodal neurons. In this study, we only analyzed the mechanical, nociceptive, and multimodal neurons that also responded to the myocardial ischemia (i.e., ischemia-sensitive mechanical, nociceptive, or multimodal neurons).
2.6. Hemodynamic measurements
Electrocardiogram was continuously recorded on the Prucka CardioLab system (GE Healthcare, Fairfield, CT, USA). A 5 French Millar Mikro-Tip pressure transducer catheter (SPR-350) was inserted into the left ventricle via the left carotid artery and the continuous data was recorded in an MPVS Ultra Pressure-Volume Loop System (Millar Instruments, Houston, TX, USA) throughout the experiment. Left ventricular systolic function was evaluated by the maximum rate of pressure change (dp/dt max).
2.7. Cardiac electrophysiological measurements
Epicardial 56-electrode nylon mesh sock electrodes were placed around the heart and unipolar electrograms (0.05−500 Hz) were measured using a Prucka CardioLab electrophysiology mapping system (GE Healthcare, Fairfield, CT, USA) (Figure 1). We assessed the activation recovery interval (ARI), which has been shown as a surrogate of local action potential duration (Millar et al., 1985). ARIs were calculated with customized software (iScalDyn, University of Utah, Salt Lake City, UT, USA), as previously described (Howard-Quijano et al., 2021). Ischemic regions were defined based on the electrodes with ST elevations and each electrogram with ST-segment changes was analyzed by semiautomated accepted software as well as manually following the guidelines described by Haws and Lux (1990). Global epicardial dispersion of repolarization (DOR) which is associated with a heterogeneity of repolarization time and an increased risk for ventricular arrhythmias, was also analyzed (Vaseghi et al., 2012).
2.8. Spinal cord stimulation
A spinal cord stimulating lead with 8 electrodes was inserted in the epidural space with the leads located at the left T1-T4 spinal cord level and the most cranial pole of the lead at T1 (Figure 1). Current controlled stimulation (model S88 stimulator; Grass Instruments, Quincy, MA, USA) was delivered at 1 kHz, 0.03 ms pulse width for 30 min. Stimulation currents were set at 90% of a motor threshold, which was determined by increasing stimulus intensity at 2 Hz and 0.4 ms pulse width until muscle contractions were observed in the shoulder.
2.9. Neuronal activity analysis
2.9.1. Spike sorting
Extracellular action potentials were recorded using a multichannel microelectrode array. Each electrode is capable of recording the activity of multiple neurons simultaneously. To be able to follow the activity of each neuron and its response to various stimuli and interventions, spike sorting was performed to identify and classify the neurons (Figure 1). Action potentials with a signal-to-noise ratio of greater than 2:1 were detected and were classified using the principal component analysis and cluster on measurements techniques in the Spike2 software program (Cambridge Electronics Design, United Kingdom). Artifacts were identified as simultaneously occurring waveforms on more than 2 adjacent electrodes as based on the specification of the multichannel electrode, it is not possible to record the neural activity from a single neuron on multiple electrodes. The artifacts were removed before the waveform classification.
2.9.2. Identification of intermediolateral neurons
Sympathetic control of the heart is initiated by the preganglionic sympathetic neurons in the IML region of the spinal cord. These preganglionic neurons project their axons into the paravertebral chain to post-ganglionic neurons contained with Stellate and Middle Cervical Ganglia. Stimulation of the paravertebral chain activates the IML neurons antidromically. To identify the IML neurons at the T2 level of the spinal cord, the T2 paravertebral chain was stimulated using a Grass S88 Stimulator (Grass Instruments, Quincy, MA, USA). To set the stimulation current intensity, first, a current threshold to induce a 10% increase in the blood pressure by stimulation at 4 Hz and 4 ms was performed. The actual stimulation to identify the IML neurons was delivered for 1 min at 1 Hz, 1 ms, and at the current threshold defined above. The neurons that had more than 60% one-to-one firing response within the 50 ms of the T2 stimulation pulse and with consistent delay were classified as the IML neurons (Dale et al., 2020; Omura et al., 2021). These neurons were recorded from an electrode zone that was in the IML region. All other recorded neurons were labeled as DH neurons.
2.9.3. Assessment of spinal neural response to myocardial ischemia
To assess the response of spinal neurons to the myocardial ischemia prior to and after SCS, the firing rate of each neuron was compared between the 1 min pre-LAD vs. during the 3 min LAD. Neurons that had a significant increase in the firing rate during the LAD ischemia prior to SCS, were identified as ischemic-sensitive neurons. To evaluate the efficacy of SCS in reducing the firing response of these neurons to the LAD ischemia, the activity of these neurons was followed during the LAD post-SCS to investigate their response change to LAD ischemia pre- vs. post-SCS.
2.9.4. Neural network interactions
In this study, we have used the measure of synchrony method that was previously developed to analyze the neural network interactions in the intrinsic cardiac nervous system (Longpre et al., 2014). In this method, we have calculated the normalized jitter-based synchrony index (SI), (Agmon, 2012) among the pair of neurons to assess the level and statistical significance of synchrony between spinal neurons. Neuron synchronization was analyzed by pair of neurons (one as a reference and one as a target). Coincidence firing between a reference and a target neuron occurred when both neurons fire within a time window of 40 ms. Using the above coincidence firing definition, some neurons with a high firing rate would cause false positive coincidence. These random coincidences are identified and the final SI is normalized based on these random coincidences that might happen by chance.
2.10. Statistical analysis
The sample size of the study was selected based on our previous published work with a similar neural recording protocol. To identify a statistically significant change in the firing rate of each neuron at different time intervals, a statistical test derived from the Skellam distribution was used (Shin et al., 2010) which was previously used for spinal and peripheral neural firing analysis (Beaumont et al., 2013; Salavatian et al., 2019; Omura et al., 2021). A statistically significant change in the activation recovery interval and dispersion of repolarization between baseline and during LAD was assessed by the paired t-test. For the statistical analysis of the neural firing data, since the data were not normally distributed, the Friedman test with Dunn’s multiple comparisons test was performed to investigate the change in firing rate during LAD and the reperfusion. The Kruskal–Wallis test was used to compare the firing rate of mechanosensitive, nociceptive, and multimodal neurons. Wilcoxon test was used to compare the non-paramteric paired synchrony changes. Statistical analyses were performed using Prism (version 8, GraphPad Software Inc., San Diego, CA, USA) and a P-value of 0.05 was considered statistically significant. Data are reported as mean ± SE.
3. Results
We have induced cardiac insult in the different regions of the heart, where chronic MI was percutaneously induced in the left lateral region by injecting microspheres into the LCX coronary artery to create a patchy infarct, and acute myocardial ischemia was induced by ligating the LAD coronary artery. All ten pigs with the MI survived for 4−5 weeks until the terminal experiment and the MI size at the terminal experiment was 10.8 ± 5.5% within the left ventricle. Acute ischemic insult measured by the number of ST elevation leads was comparable between ischemia at pre- and post-SCS. The hemodynamic changes during SCS were not significant (Table 1). No difference was observed between the hemodynamic response of ischemia and reperfusion pre-SCS vs. post-SCS (Table 2).
3.1. Effect of myocardial ischemia on activation recovery interval and dispersion of repolarization
The cardiac electrophysiological data of nine pigs were analyzed and the data from one pig was removed due to the poor quality of the cardiac electrophysiological signal. The ARI is a surrogate for the action potential duration and in this study, we have measured the ARI at baseline and during the LAD to evaluate the magnitude of sympathoexcitation caused by the myocardial ischemia. Sympathetic activation is associated with shortened activation recovery interval duration and increased dispersion of repolarization (Yamakawa et al., 2016). ARI in the ischemic zone was shortened during the myocardial ischemia (baseline:314.3 ± 13.6 ms, LAD: 210.1 ± 7.2 ms, P < 0.0001) indicating the sympathoexcitation (Figure 2). The level of arrhythmogenicity was assessed by measuring the DOR at baseline and during myocardial ischemia. The DOR was increased during myocardial ischemia (baseline: 443.9 ± 93.05 ms2 to LAD: 1941 ± 213.8 ms2, P < 0.0001) showing increased arrhythmogenicity during the LAD ischemia.
Figure 2. Effect of myocardial ischemia on activation recovery interval and dispersion of repolarization. (A) Representative activation recovery interval (ARI) map during the left anterior descending coronary artery (LAD) ischemia pre- and post-spinal cord stimulation (SCS). (B) ARI shortening during LAD ischemia was reduced after SCS, indicating less sympathoexcitation during myocardial ischemia post-SCS (N = 9 pigs). (C) SCS is decreasing the arrhythmogenicity by reducing the augmentation of dispersion of repolarization (DOR) during myocardial ischemia (N = 9 pigs). *Denotes P < 0.05.
3.2. Spinal cord stimulation effect on sympathoexcitation and arrhythmogenicity
To evaluate the effect of SCS on sympathoexcitation and arrhythmogenicity, the extent of ARI shortening and DOR augmentation were compared during the LAD ischemia before and after 30 min of SCS. SCS suppressed the sympathoexcitation during LAD post-SCS (ARI shortening: pre-SCS: −104.10 ± 8.04 ms; post-SCS: −87.80 ± 8.75 ms, P = 0.028) (Figure 2). The arrhythmogenicity was also reduced by SCS, manifested by less augmentation in the DOR during LAD post-SCS (DOR change: pre-SCS: 1497 ± 201.9 ms2; post-SCS: 1260 ± 164.7 ms2, P = 0.024 vs. pre-SCS) (Figure 2).
3.3. Effect of spinal cord stimulation on spinal neural response to myocardial ischemia
Using the multichannel microelectrode array, 2,758 spinal neurons were recorded. Spinal neurons (n = 257) that were activated by myocardial ischemia were labeled as ischemia-sensitive neurons. IML neurons (n = 29) and DH neurons (n = 228) were identified, and their firing activities were followed during the whole protocol. The activity of ischemia-sensitive spinal neurons (n = 257) was evaluated during the myocardial ischemia post-SCS (Figure 3). The ischemia-sensitive neurons’ response to the LAD ischemia and reperfusion was reduced post-SCS (pre-SCS: LAD: 2.45 ± 0.25 Hz, reperfusion: 2.34 ± 0.402 Hz; post-SCS: LAD: −0.07 ± 0.11 Hz, P < 0.0001 vs. LAD pre-SCS, reperfusion: −0.14 ± 0.17 Hz, P < 0.0001 vs. reperfusion pre-SCS). The response of IML neurons (Pre-SCS: LAD: 3.50 ± 0.53 Hz, reperfusion: 4.19 ± 1.24 Hz; Post-SCS: LAD: 0.705 ± 0.67 Hz, P = 0.0027 vs. LAD pre-SCS, reperfusion: 0.29 ± 0.96 Hz, P = 0.0027 vs. reperfusion pre-SCS) and DH neurons (Pre-SCS: LAD: 2.27 ± 0.27 Hz, reperfusion: 2.29 ± 0.39 Hz; Post-SCS: LAD: −0.16 ± 0.104 Hz, P < 0.0001 vs. LAD pre-SCS, reperfusion: −0.20 ± 0.16 Hz, P < 0.0001 vs. reperfusion pre-SCS) to myocardial ischemia and reperfusion was also blunted after the SCS.
Figure 3. Effect of spinal cord stimulation on firing response of spinal neurons to myocardial ischemia. (A) Representative neural activity of five spinal neurons during the left anterior descending coronary artery (LAD) ischemia pre- and post-spinal cord stimulation (SCS). The activity of these five neurons increases during the ischemia, however, they have a different temporal response to the myocardial ischemia and reperfusion. A 30-s averaging window was used for the firing rate data in this panel. (B) The firing rate change (firing rate during LAD or Rep. - firing rate at baseline) of spinal neurons (n = 257) to LAD and reperfusion (Rep.) is shown pre and post-SCS. Negative values, indicate the decreasing firing rate. Data is presented as mean ± SE. ****Denotes adjusted P-value < 0.0001 after Dunn’s multiple comparisons tests.
3.4. Response of effect of spinal cord stimulation on different types of DH neurons
Mechanosensitive (n = 49), nociceptive (n = 37), and multimodal (n = 59) ischemia-sensitive DH neurons were characterized. To investigate if the response of mechanosensitive, nociceptive, and multimodal neurons to LAD ischemia is affected differently by the SCS, we have evaluated the effect of SCS on the firing response of each of these groups to the LAD ischemia and reperfusion (Figure 4). SCS mitigated the response of mechanosensitive (Pre-SCS: LAD: 1.91 ± 0.29 Hz, reperfusion: 2.51 ± 0.60 Hz; Post-SCS: LAD: −0.32 ± 0.26 Hz, P < 0.0001 vs. LAD pre-SCS, reperfusion: −0.25 ± 0.26 Hz, P < 0.0001 vs. reperfusion pre-SCS), nociceptive (Pre-SCS: LAD: 1.29 ± 0.30 Hz, reperfusion: 1.55 ± 0.38 Hz; Post-SCS: LAD: 0.50 ± 0.28 Hz, P = 0.009 vs. LAD pre-SCS, reperfusion: 0.72 ± 0.47 Hz, P = 0.02 vs. reperfusion pre-SCS) and multimodal (Pre-SCS: LAD: 3.62 ± 0.75 Hz, reperfusion: 4.70 ± 1.34 Hz; Post-SCS: LAD: −0.72 ± 0.27 Hz, P < 0.0001 vs. LAD pre-SCS, reperfusion: −1.15 ± 0.44 Hz, P < 0.0001 vs. reperfusion pre-SCS) neurons to the myocardial ischemia.
Figure 4. Effect of spinal cord stimulation on the intermediolateral and dorsal horn neuronal hyperactivity during myocardial ischemia. Spinal cord stimulation (SCS) suppressed the response of panel (A) intermediolateral column (IML) neurons (n = 28, the activity of 1 neuron was removed based on the outlier criteria by ROUT method and Q = 0.1%), (B) dorsal horn (DH) neurons (n = 228), and (C) ischemia-sensitive mechanosensitive (n = 49), nociceptive (n = 37) and multimodal neurons (n = 59) to the left anterior descending coronary artery (LAD) ischemia and reperfusion (Rep.). Data is presented as mean ± SE and *, **, ****denotes P < 0.05, P < 0.01, and p < 0.0001, respectively.
3.5. Neural synchrony response alterations among DH-DH and DH-IML neurons to SCS
The neural synchrony between the DH and IML neurons was evaluated during LAD ischemia, and reperfusion pre- and post-SCS (Figure 5). The LAD ischemia-induced augmentation in neuronal synchrony between DH-DH (pre-SCS:985.30 ± 351.39, vs. post-SCS: 39.60 ± 17.26, p = 0.004) and DH-IML (pre-SCS:165.60 ± 116.76, vs. post-SCS: 11.50 ± 9.45, p = 0.02) pairs of neurons were mitigated by the SCS. SCS also suppressed the reperfusion-induced augmentation in neuronal synchrony between DH-DH (pre-SCS: 1001.60 ± 393.85, vs. post-SCS: 31.30 ± 15.28, p = 0.008) and DH-IML (pre-SCS:175.20 ± 128.41, vs. post-SCS: 13.40 ± 11.49, p = 0.02).
Figure 5. Effect of spinal cord stimulation on spinal neuronal synchrony during myocardial ischemia. (A) Representative activity of two neurons when they have synchronized and unsynchronized activity. (B) The augmentation in the neuronal synchrony among DH-DH neurons (n = 2668, n = 2668) during left anterior descending coronary artery (LAD) ischemia and reperfusion (Rep.) was suppressed by the spinal cord stimulation (SCS). (C) The SCS showed a similar effect on suppressing the synchrony among DH-IML pair of neurons (DH: n = 2668, IML: n = 90) synchrony during LAD, and reperfusion. *Denotes P < 0.05, **denotes P < 0.01.
4. Discussion
We have developed a pre-clinical large animal model that mimics a clinical situation in patients with prior myocardial infarction who experiences episodes of ventricular arrhythmias due to the new onset of acute coronary artery ischemia. This is the first study, to our knowledge, that shows the therapeutic effect of high-frequency SCS on cardiac sympathoexcitation and arrhythmogenicity during acute myocardial ischemia in the presence of chronic MI.
Our major findings are as follows: in a chronic MI heart (1) 3 min of acute myocardial ischemia shortened ARI, which was mitigated by preemptive SCS therapy, (2) acute myocardial ischemia increased DOR, which was mitigated by preemptive SCS therapy. In addition, (3) IML sympathetic preganglionic neural firing rate increased with acute myocardial ischemia and reperfusion injury, which was blunted with preemptive SCS (4) the firing response of the mechanosensitive, nociceptive, and multimodal neurons to myocardial ischemia and reperfusion were mitigated by preemptive SCS therapy, (5) myocardial ischemia and reperfusion increased the synchrony among the DH-DH and DH-IML neurons and this effect was blocked by SCS.
4.1. Effect of SCS on cardiac electrophysiology
We and others have previously shown that SCS can decrease cardiac sympathoexcitation and arrhythmogenicity induced by myocardial ischemia in a healthy heart (Lopshire et al., 2009; Odenstedt et al., 2011; Howard-Quijano et al., 2017, 2021). While the underlying mechanisms are not exactly elucidated yet, we have found that SCS increases an inhibitory neurotransmitter GABA within the thoracic spinal cord, which may decrease an efferent outflow to the heart, thus reducing cardiac sympathoexcitation and ventricular arrhythmias (Howard-Quijano et al., 2023). In this study, we have shown the effectiveness of the high-frequency SCS in reducing acute myocardial ischemia-induced sympathoexcitation and arrhythmogenicity in subjects with prior chronic myocardial infarction. Although we used the low-frequency (50 Hz) SCS in our previous studies and showed the efficacy of SCS to prevent arrhythmogenicity in pre-clinical models, a high-frequency of SCS such as 1 kHz may provide more superior effects on cardiac pathophysiology as a randomized controlled trial comparing the high-frequency and conventional stimulation demonstrated the superiority of the high-frequency stimulation for the treatment in the pain field (Kapural et al., 2015).
4.2. Impact of SCS on spinal neural firing rate during myocardial ischemia
Spinal cord stimulation as a bioelectric neuromodulation engages central neuronal reflexes (Ardell, 2016) to influence ischemia-sensitive spinal cord neurons with regard to both pain perception (Qin et al., 2008) and modulating central-mediated sympathoexcitation (Ding et al., 2008). SCS can mitigate myocardial ischemia-induced reflex activation of cardiac neural networks (Foreman et al., 2000; Ardell et al., 2009). These reflexes are known to exhibit memory, with neural network protective effects extending up to 1 h after SCS (Armour et al., 2002). We have previously shown that acute myocardial ischemia activated cardiac-related DRG neurons and SCS decreased basal activity of such responsiveness to transient myocardial ischemia in the thoracic DRG in a healthy heart (Salavatian et al., 2019).
In this study, we showed that ischemia-activated neurons in the dorsal horn and IML regions of the spinal cord were suppressed by high-frequency SCS, suggesting that SCS mitigated cardiac efferent output from the spinal cord to the heart, which would contribute to lower sympathetic excitation and reduce the risk of ventricular arrhythmias. Interestingly, this mitigation of neuronal firing was also seen during the reperfusion period, which causes tremendous pathophysiological damage to the heart in various ways (Soares et al., 2019), thus indicating that SCS has a memory effect to reduce cardiac sympathoexcitation even after cessation of SCS.
Moreover, we have shown the impact of SCS on blunting the LAD ischemia-induced neuronal hyperexcitation is not significantly different among the ischemia-sensitive mechanical, nociceptive, and multimodal neurons, suggesting that high-frequency SCS is not selective in suppressing just the neural activity of ischemia-sensitive neurons.
4.3. Impact of SCS on spinal neural synchrony during myocardial ischemia
We have previously shown that myocardial ischemia increases the synchrony among the dorsal horn and IML pair of neurons (Dale et al., 2020; Omura et al., 2021). We hypothesize that this increase in synchrony and augmented neuronal communication occurs because of the higher processing needs of the spinal cord during myocardial ischemia. In this study, we showed that this augmented synchrony during both myocardial ischemia and reperfusion is decreased by the SCS. It is important to mention that in this study, we have used a normalized synchrony method. This method normalizes the synchrony index in such a way that the random firing of two high-firing neurons does not necessarily result in a high synchrony index among these neurons.
Myocardial ischemia activates cardiac afferent inputs to the thoracic spinal cord through dorsal root ganglia (Malliani et al., 1981; Fukuda et al., 2015), and it associates with increased coordination within dorsal horn neurons and between dorsal horn to sympathetic preganglionic neurons in the thoracic spinal cord (Dale et al., 2020). In this study, we have also observed the augmentation in the synchronized pairs of neurons in the spinal cord afferent processing sites (between DH and DH neurons) and between the spinal afferent and efferent processing sites (between DH and IML neurons). This suggests that during the myocardial ischemia-reperfusion injury, there is a high level of spinal processing which is manifested by a high firing rate in the DH and high synchrony among the DH-DH neurons and high synchrony among the DH-IML neurons, which leads to the increased firing rate of IML neurons and sympathoexcitation.
4.4. Clinical implications
The clinical trials of the SCS for heart diseases have been equivocal. In the heart failure patients, while the DEFEAT-HF trial showed no improvement in left ventricular end-systolic volume index, peak oxygen uptake, N-terminal pro-B-type natriuretic peptide [NT-proBNP] endpoints (Zipes et al., 2016), the SCS HEART trial demonstrated improvement in New York Heart Association (NYHA) classification, quality of life, left ventricular ejection fraction, left ventricular volume, and peak oxygen consumption (Tse et al., 2015). This equivocal outcome could be a result of not selecting optimized stimulation parameters due to the lack of understanding of the mechanism of the SCS in the modulation of the spinal neural network. In this study, we have evaluated the efficacy of SCS in mitigating the myocardial ischemia-induced sympathoexcitation and arrhythmogenicity in the subjects with myocardial infarction to evaluate if the SCS therapy after a chronic myocardial infarction can provide cardioprotection against the new myocardial ischemic events. The candidates for this therapy would be the patients who have developed a chronic MI and are seeking to find a therapy to decrease the chance of reoccurrence of the fatal ischemic-induced ventricular arrhythmias.
4.5. Limitations
We induced chronic MI into the LCX coronary artery in order to allow for a second acute myocardial ischemia into the LAD coronary artery. It is possible that chronic myocardial infarction into the LAD region creates a different arrhythmogenicity and neuronal remodeling in the cardio-spinal neural pathway. General anesthesia can suppress neuronal activity but is necessary to conduct the animal experiment. Isoflurane was switched to alpha-chloralose to minimize the effects of isoflurane on spinal neural activity. Since this study was performed in an anesthetized animal at rest, in which the sympathetic tone is low, we could not evaluate the effect of SCS on the sympathetic function without the myocardial ischemia. Further studies in normal awake animals are needed to evaluate the impact of the SCS on normal sympathetic function. Another limitation of this study is that since the memory effect of 30 min- SCS was not known, sham ischemia was always performed before the ischemia + SCS.
5. Conclusion
Myocardial ischemia induces sympathoexcitation by increasing the IML activity which is caused by the activation of DH neurons and increased spinal neural synchrony between DH and IML neurons. SCS mitigates the adverse effect of myocardial ischemia by suppressing spinal neural synchrony and IML hyperactivity.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
This animal study was reviewed and approved by the Institutional Animal Care and Use Committee (IACUC).
Author contributions
SS, JAA, JLA, RF, and AM conceived and designed the research. SS, YK, BW, and JF conducted the experiments and analyzed the data. SS, YK, JAA, JLA, RF, KH-Q, and AM interpreted the results of the experiments. SS and YK drafted the manuscript and prepared the figures. All authors reviewed and approved the final version of the manuscript.
Funding
SS was supported by the Competitive Medical Research Fund by the University of Pittsburgh. AM was supported by the NIH R01 HL136836 and NIH R44 DA049630.
Acknowledgments
We thank Eevanna Lundquist, Alyssa Falcione, Stacy Cashman, and Madison Fisher for assistance with the animal experiments. Parts of Figure 1 were created with BioRender.com.
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Agmon, A. (2012). A novel, jitter-based method for detecting and measuring spike synchrony and quantifying temporal firing precision. Neural Syst Circuits 2:5. doi: 10.1186/2042-1001-2-5
Ardell, J. L. (2016). Heart failure: mechanisms of spinal cord neuromodulation for heart disease. Nat. Rev. Cardiol. 13, 127–128.
Ardell, J. L., Cardinal, R., Vermeulen, M., and Armour, J. A. (2009). Dorsal spinal cord stimulation obtunds the capacity of intrathoracic extracardiac neurons to transduce myocardial ischemia. Am. J. Physiol. Regul. Integr. Comp. Physiol. 297, R470–R477. doi: 10.1152/ajpregu.90821.2008
Ardell, J. L., Foreman, R. D., Armour, J. A., and Shivkumar, K. (2019). Cardiac sympathectomy and spinal cord stimulation attenuate reflex-mediated norepinephrine release during ischemia preventing ventricular fibrillation. JCI Insight 4:e131648. doi: 10.1172/jci.insight.131648
Armour, J. (2008). Potential clinical relevance of the ‘little brain’ on the mammalian heart. Exp. Physiol. 93, 165–176. doi: 10.1113/expphysiol.2007.041178
Armour, J. A. (2004). Cardiac neuronal hierarchy in health and disease. Am. J. Physiol. Regul. Integr. Comp. Physiol. 287, R262–R271. doi: 10.1152/ajpregu.00183.2004
Armour, J. A., Linderoth, B., Arora, R. C., Dejongste, M. J. L., Ardell, J. L., Kingma, J. G., et al. (2002). Long-term modulation of the intrinsic cardiac nervous system by spinal cord neurons in normal and ischaemic hearts. Auton. Neurosci. 95, 71–79. doi: 10.1016/s1566-0702(01)00377-0
Beaumont, E., Salavatian, S., Southerland, E. M., Vinet, A., Jacquemet, V., Armour, J. A., et al. (2013). Network interactions within the canine intrinsic cardiac nervous system: implications for reflex control of regional cardiac function. J. Physiol. 591, 4515–4533. doi: 10.1113/jphysiol.2013.259382
Bourke, T., Vaseghi, M., Michowitz, Y., Sankhla, V., Shah, M., Swapna, N., et al. (2010). Neuraxial modulation for refractory ventricular arrhythmias: value of thoracic epidural anesthesia and surgical left cardiac sympathetic denervation. Circulation 121, 2255–2262. doi: 10.1161/CIRCULATIONAHA.109.929703
Dale, E. A., Kipke, J., Kubo, Y., Sunshine, M. D., Castro, P. A., Ardell, J. L., et al. (2020). Spinal cord neural network interactions: implications for sympathetic control of the porcine heart. Am. J. Physiol. Heart Circ. Physiol. 318, H830–H839. doi: 10.1152/ajpheart.00635.2019
Ding, X., Ardell, J. L., Hua, F., Mcauley, R. J., Sutherly, K., Daniel, J. J., et al. (2008). Modulation of cardiac ischemia-sensitive afferent neuron signaling by preemptive C2 spinal cord stimulation: effect on substance P release from rat spinal cord. Am J Physiol Regul Integr Comp Physiol. 294, R93–R101. doi: 10.1152/ajpregu.00544.2007
Foreman, R. D., Linderoth, B., Ardell, J. L., Barron, K. W., Chandler, M. J., Hull, S. S. Jr., et al. (2000). Modulation of intrinsic cardiac neurons by spinal cord stimulation: implications for its therapeutic use in angina pectoris. Cardiovasc Res. 47, 367–375. doi: 10.1016/s0008-6363(00)00095-x
Fukuda, K., Kanazawa, H., Aizawa, Y., Ardell, J. L., and Shivkumar, K. (2015). Cardiac innervation and sudden cardiac death. Circ. Res. 116, 2005–2019. doi: 10.1161/CIRCRESAHA.116.304679
Haws, C. W., and Lux, R. L. (1990). Correlation between in vivo transmembrane action potential durations and activation-recovery intervals from electrograms. Effects of interventions that alter repolarization time. Circulation 81, 281–288. doi: 10.1161/01.cir.81.1.281
Howard-Quijano, K., Kuwabara, Y., Yamaguchi, T., Roman, K., Salavatian, S., Taylor, B., et al. (2023). GABAergic signaling during spinal cord stimulation reduces cardiac arrhythmias in a porcine model. Anesthesiology 138, 372–387. doi: 10.1097/ALN.0000000000004516
Howard-Quijano, K., Takamiya, T., Dale, E. A., Kipke, J., Kubo, Y., Grogan, T., et al. (2017). Spinal cord stimulation reduces ventricular arrhythmias during acute ischemia by attenuation of regional myocardial excitability. Am. J. Physiol. Heart Circ. Physiol. 313, H421–H431. doi: 10.1152/ajpheart.00129.2017
Howard-Quijano, K., Yamaguchi, T., Gao, F., Kuwabara, Y., Puig, S., Lundquist, E., et al. (2021). Spinal cord stimulation reduces ventricular arrhythmias by attenuating reactive gliosis and activation of spinal interneurons. JACC Clin. Electrophysiol. 7, 1211–1225. doi: 10.1016/j.jacep.2021.05.016
Kapural, L., Yu, C., Doust, M. W., Gliner, B. E., Vallejo, R., Sitzman, B. T., et al. (2015). Novel 10-kHz high-frequency therapy (HF10 Therapy) is superior to traditional low-frequency spinal cord stimulation for the treatment of chronic back and leg pain. Anesthesiology 123, 851–860.
Longpre, J. P., Salavatian, S., Beaumont, E., Armour, J. A., Ardell, J. L., and Jacquemet, V. (2014). Measure of synchrony in the activity of intrinsic cardiac neurons. Physiol. Meas. 35, 549–566. doi: 10.1088/0967-3334/35/4/549
Lopshire, J. C., Zhou, X., Dusa, C., Ueyama, T., Rosenberger, J., Courtney, N., et al. (2009). Spinal cord stimulation improves ventricular function and reduces ventricular arrhythmias in a canine postinfarction heart failure model. Circulation 120, 286–294. doi: 10.1161/CIRCULATIONAHA.108.812412
Luckl, J., Keating, J., and Greenberg, J. H. (2008). Alpha-chloralose is a suitable anesthetic for chronic focal cerebral ischemia studies in the rat: a comparative study. Brain Res. 1191, 157–167. doi: 10.1016/j.brainres.2007.11.037
Luqman, N., Sung, R. J., Wang, C. L., and Kuo, C. T. (2007). Myocardial ischemia and ventricular fibrillation: pathophysiology and clinical implications. Int. J. Cardiol. 119, 283–290. doi: 10.1016/j.ijcard.2006.09.016
Malliani, A., Lombardi, F., and Pagani, M. (1981). Functions of afferents in cardiovascular sympathetic nerves. J. Auton. Nerv. Syst. 3, 231–236. doi: 10.1016/0165-1838(81)90065-5
Millar, C. K., Kralios, F. A., and Lux, R. L. (1985). Correlation between refractory periods and activation-recovery intervals from electrograms: effects of rate and adrenergic interventions. Circulation 72, 1372–1379. doi: 10.1161/01.cir.72.6.1372
Odenstedt, J., Linderoth, B., Bergfeldt, L., Ekre, O., Grip, L., Mannheimer, C., et al. (2011). Spinal cord stimulation effects on myocardial ischemia, infarct size, ventricular arrhythmia, and noninvasive electrophysiology in a porcine ischemia-reperfusion model. Heart Rhythm 8, 892–898. doi: 10.1016/j.hrthm.2011.01.029
Omura, Y., Kipke, J. P., Salavatian, S., Afyouni, A. S., Wooten, C., Herkenham, R. F., et al. (2021). Spinal anesthesia reduces myocardial ischemia-triggered ventricular arrhythmias by suppressing spinal cord neuronal network interactions in pigs. Anesthesiology 134, 405–420. doi: 10.1097/ALN.0000000000003662
Qin, C., Farber, J. P., Linderoth, B., Shahid, A., and Foreman, R. D. (2008). Neuromodulation of thoracic intraspinal visceroreceptive transmission by electrical stimulation of spinal dorsal column and somatic afferents in rats. J. Pain 9, 71–78. doi: 10.1016/j.jpain.2007.08.007
Rajendran, P. S., Nakamura, K., Ajijola, O. A., Vaseghi, M., Armour, J. A., Ardell, J. L., et al. (2016). Myocardial infarction induces structural and functional remodelling of the intrinsic cardiac nervous system. J. Physiol. 594, 321–341.
Rubart, M., and Zipes, D. P. (2005). Mechanisms of sudden cardiac death. J. Clin. Invest. 115, 2305–2315.
Salavatian, S., Ardell, S. M., Hammer, M., Gibbons, D., Armour, J. A., and Ardell, J. L. (2019). Thoracic spinal cord neuromodulation obtunds dorsal root ganglion afferent neuronal transduction of the ischemic ventricle. Am. J. Physiol. Heart Circ. Physiol. 317, H1134–H1141. doi: 10.1152/ajpheart.00257.2019
Salavatian, S., Hoang, J. D., Yamaguchi, N., Lokhandwala, Z. A., Swid, M. A., Armour, J. A., et al. (2022). Myocardial infarction reduces cardiac nociceptive neurotransmission through the vagal ganglia. JCI Insight 7:e155747. doi: 10.1172/jci.insight.155747
Shen, M. J., and Zipes, D. P. (2014). Role of the autonomic nervous system in modulating cardiac arrhythmias. Circ. Res. 114, 1004–1021.
Shin, H. C., Aggarwal, V., Acharya, S., Schieber, M. H., and Thakor, N. V. (2010). Neural decoding of finger movements using Skellam-based maximum-likelihood decoding. IEEE Trans. Biomed. Eng. 57, 754–760. doi: 10.1109/TBME.2009.2020791
Soares, R. O. S., Losada, D. M., Jordani, M. C., Évora, P., and Castro, E. S. O. (2019). Ischemia/reperfusion injury revisited: an overview of the latest pharmacological strategies. Int. J. Mol. Sci. 20:5034. doi: 10.3390/ijms20205034
Takagawa, J., Zhang, Y., Wong, M. L., Sievers, R. E., Kapasi, N. K., Wang, Y., et al. (2007). Myocardial infarct size measurement in the mouse chronic infarction model: comparison of area- and length-based approaches. J. Appl. Physiol. 102, 2104–2111. doi: 10.1152/japplphysiol.00033.2007
Tse, H. F., Turner, S., Sanders, P., Okuyama, Y., Fujiu, K., Cheung, C. W., et al. (2015). Thoracic spinal cord stimulation for heart failure as a restorative treatment (SCS HEART study): first-in-man experience. Heart Rhythm 12, 588–595. doi: 10.1016/j.hrthm.2014.12.014
Vaseghi, M., Gima, J., Kanaan, C., Ajijola, O. A., Marmureanu, A., Mahajan, A., et al. (2014). Cardiac sympathetic denervation in patients with refractory ventricular arrhythmias or electrical storm: intermediate and long-term follow-up. Heart Rhythm 11, 360–366. doi: 10.1016/j.hrthm.2013.11.028
Vaseghi, M., Lux, R. L., Mahajan, A., and Shivkumar, K. (2012). Sympathetic stimulation increases dispersion of repolarization in humans with myocardial infarction. Am. J. Physiol. Heart Circ. Physiol. 302, H1838–H1846. doi: 10.1152/ajpheart.01106.2011
Yamakawa, K., Howard-Quijano, K., Zhou, W., Rajendran, P., Yagishita, D., Vaseghi, M., et al. (2016). Central vs. peripheral neuraxial sympathetic control of porcine ventricular electrophysiology. Am. J. Physiol. Regul. Integr. Comp. Physiol. 310, R414–R421. doi: 10.1152/ajpregu.00252.2015
Zipes, D. P. (1990). Influence of myocardial ischemia and infarction on autonomic innervation of heart. Circulation 82, 1095–1105. doi: 10.1161/01.CIR.82.4.1095
Keywords: spinal cord stimulation, neuromodulation, myocardial ischemia, autonomic nervous system, arrhythmias, spinal neural network
Citation: Salavatian S, Kuwabara Y, Wong B, Fritz JR, Howard-Quijano K, Foreman RD, Armour JA, Ardell JL and Mahajan A (2023) Spinal neuromodulation mitigates myocardial ischemia-induced sympathoexcitation by suppressing the intermediolateral nucleus hyperactivity and spinal neural synchrony. Front. Neurosci. 17:1180294. doi: 10.3389/fnins.2023.1180294
Received: 06 March 2023; Accepted: 16 May 2023;
Published: 02 June 2023.
Edited by:
Leah Reznikov, University of Florida, United StatesReviewed by:
Prasanth Puthanveetil, Midwestern University, United StatesAlexander L. Green, University of Oxford, United Kingdom
Copyright © 2023 Salavatian, Kuwabara, Wong, Fritz, Howard-Quijano, Foreman, Armour, Ardell and Mahajan. 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: Aman Mahajan, mahajana@upmc.edu