- 1Lawrence D. Longo, MD Center for Perinatal Biology, Loma Linda University School of Medicine, Loma Linda, CA, United States
- 2Advanced Imaging and Microscopy Core, Loma Linda University School of Medicine, Loma Linda, CA, United States
- 3Department of Biostatistics, School of Medicine, California Northstate University, Elk Grove, CA, United States
Cerebral arterial vasoreactivity is vital to the regulation of cerebral blood flow. Depolarization of arterial myocytes elicits whole-cell Ca2+ oscillations as well as subcellular Ca2+ sparks due to activation of ryanodine receptors on the sarcoplasmic reticulum. Previous evidence illustrates that contraction of cerebral arteries from sheep and underlying Ca2+ signaling pathways are modified by age and that long-term hypoxia (LTH) causes aberrations in Ca2+ signaling pathways and downstream effectors impacting vasoregulation. We hypothesize that age and LTH affect the influence of membrane depolarization on whole-cell intracellular Ca2+ oscillations and sub-cellular Ca2+ spark activity in cerebral arteries. To test this hypothesis, we examined Ca2+ oscillatory and spark activities using confocal fluorescence imaging techniques of Fluo-4 loaded basilar arterial myocytes of low- and high-altitude term fetal (∼145 days of gestation) and adult sheep, where high-altitude pregnant and non-pregnant sheep were placed at 3,801 m for >100 days. Ca2+ oscillations and sparks were recorded using an in situ preparation evaluated in the absence or presence of 30 mM K+ (30K) to depolarize myocytes. Myocytes from adult animals tended to have a lower basal rate of whole-cell Ca2+ oscillatory activity and 30K increased the activity within cells. LTH decreased the ability of myocytes to respond to depolarization independent of age. These observations illustrate that both altitude and age play a role in affecting whole-cell and localized Ca2+ signaling, which are important to arterial vasoreactivity and cerebral blood flow.
Introduction
In an age where cardiovascular health is so heavily scrutinized, we still struggle to reduce deaths from cardiovascular diseases. With the identification of modifiable risk factors such as smoking, diet, exercise, and obesity, we are now able to help prevent a significant portion of the morbidity and mortality associated with cardiovascular disease; however, environmental exposures that we have little or no control over are important to cardiovascular development and health (Yang et al., 2017). What is more, previous research illustrates that many cardiovascular diseases have origins in fetal development and that gestational stressors have complex effects and impacts on vascular health (Ducsay et al., 2018).
Gestational long-term hypoxia (LTH) from living at high altitude or placental insufficiency of various etiologies has a marked effect on fetal cerebral vascular development that can have repercussions throughout an individual’s lifetime (Ducsay et al., 2018). Decreasing oxygenation during the fetal period plays a significant role in vascular development, affecting vessel structure and function. Low oxygen exposures cause significant changes in the regulation of cellular phenotype, which influences vasoreactivity (Adeoye et al., 2015; Ducsay et al., 2018). Cerebral vascular dysfunction can result in an increased risk of cerebral edema, hemorrhage, and ischemic stroke. Furthermore, there are correlations between brain blood flow, developmental neuropathies, and impaired learning in children as well as psychological disorders such as clinical depression (Koenigs and Grafman, 2009) and dementia (Roher et al., 2012). These interrelationships illustrate the importance of considering the role of dysregulated vessel structure and function in neurological disorders.
Traveling to high altitude has significant effects on cerebrovascular blood flow (CBF). Cerebral arteries dilate in response to acute hypoxia, allowing for distribution of a greater portion of cardiac output to the brain in an attempt to maintain oxygen delivery despite reduced oxygen tensions (Hunter et al., 2003; Lee et al., 2009; Giussani, 2016). However, the reduced cerebral arterial contractility impairs the autoregulation of CBF (Tweed et al., 1983, 1986) and increases the risk of hemorrhage of the germinal matrix leading to varying brain disorders including epilepsy, cerebral palsy, and intellectual disability (Volpe, 1997; Nelson, 2003; Ferriero, 2004; Stoll et al., 2010). Vasodilation of cerebral vessels in response to acute hypoxia is partially regulated by an increase in large conductance potassium (BKCa) channel activation, which is responsible for membrane hyperpolarization and consequent vasorelaxation (Gebremedhin et al., 1994; Brenner et al., 2000; Ledoux et al., 2006; Thorpe et al., 2017). Long-term hypoxia, in comparison to acute hypoxia, can lead to compensatory responses that reduce CBF toward control levels. Studies done on adult animals as well as humans show that a large component of this compensation is due to increased ventilation and hematocrit, both of which serve to maintain cerebral O2 delivery (Brugniaux et al., 2007). Prior to these adaptations and acclimatization, there is a correlation between hypoxia and increased basilar blood flow due to increases in arterial diameter (Jansen et al., 2002; Ainslie and Subudhi, 2014). Studies in fetal animals have shown that gestational LTH leads to a sustained increase in the proportion of cardiac output serving the brain (Ducsay et al., 2018). However, there is an impaired ability to increase cerebral blood flow in response to superimposed acute hypoxia (Pena et al., 2007), suggestive of altered cerebral vasoactivity. This is consistent with evidence that chronic hypoxia alters many of the pathways that regulate calcium signaling, and thus vascular tone, in fetal cerebral arteries (Pearce, 2006).
The current series of studies were designed to interrogate the effects of animal age and long-term hypoxia on calcium signals in basilar arterial myocytes under resting conditions and in response to membrane depolarization. We hypothesized that long-term hypoxia due to high altitude exposure would alter calcium signaling in ways that may restrict BKCa channel activation. We theorized that in arteries from low altitude fetuses, myocyte depolarization would increase whole-cell Ca2+ oscillations and localized rapid calcium signals (Ca2+ sparks) but that LTH would blunt these responses. We further hypothesized that post-natal development would magnify the increases in whole-cell and localized Ca2+ responses to membrane depolarization. This was examined by interrogating Ca2+ signals in basilar cerebral arterial myocytes of fetal and adult sheep housed at low- or high-altitude.
Materials and Methods
Experimental Animals
Surgical and experimental procedures were performed in accordance with the regulations of the Animal Welfare Act, the National Institutes of Health’s Guide for the Care and Use of Laboratory Animals, and “The Guiding Principles in the Care and Use of Animals” approved by the Council of the American Physiological Society and by the Institutional Animal Care and Use Committee of Loma Linda University. Pregnant (n = 7) and non-pregnant ewes (n = 7) of a mixed Western breed were divided into low altitude (normoxic) and high-altitude long-term hypoxic (LTH) groups. All ewes were obtained from Nebeker Ranch in Lancaster, CA, United States at an elevation near sea level, 720 m. Normoxic control ewes were maintained near sea level at 720 m for the duration of their gestation. Animals for the LTH groups were held at Nebeker Ranch under normoxic conditions until 30 days gestation at which time the pregnant and non-pregnant ewes were transported to the Barcroft Laboratory, White Mountain Research Station in Bishop, CA, United States at an elevation of 3,801 m. Previous work shows that residing at the Barcroft laboratory results in a maternal arterial PO2 of 60 ± 3 Torr and a fetal arterial PO2 of 19 ± 2 Torr (Kamitomo et al., 1993). Animals were held at elevation for the remaining ∼110 days of gestation for pregnant ewes and acclimatization for the non-pregnant ewes. Following this acclimatization period, ewes were transported (∼6 h drive) to Loma Linda University (LLU) for study at an elevation of 346 m. Once at LLU, LTH ewes were surgically instrumented with arterial and tracheal catheters. Based on frequent arterial blood gas sampling and adjustment of the rate of N2 flow through the tracheal catheter, arterial PO2 level in the adult sheep was maintained at ∼60 Torr for 2–4 days, mimicking the conditions of the effects of high altitude until the day of study (Kamitomo et al., 1993). After induction with thiopental sodium (10 mg/kg iv), the ewes were intubated and anesthesia was maintained via inhalation of 2–3% isoflurane in O2 for the duration of the surgery. The fetuses were delivered via hysterectomy at a male to female ratio of ∼1:1. Following delivery, fetal sheep were euthanized with an overdose of Euthasol (pentobarbital sodium, 100 mg/kg) and phenytoin sodium (10 mg/kg).
Artery Isolation
Non-pregnant females and mixed sex, near term fetal brains were removed and placed in iced Balanced Salt Solution (BSS) of the following composition (mM): 126 NaCl; 5 KCl; 10 HEPES; 1 MgCl2; 2 CaCl2; 10 glucose; pH 7.4 (adjusted with NaOH). Basilar arteries then were quickly dissected under normoxic conditions in BSS. Basilar arteries were selected from the same anatomical locations in both fetal and adult sheep to maintain segments of similar function and embryological origin. Because of this, there was a significant difference in diameter between fetal and adult arteries (∼200 μm vs. 300 μm, respectively) as described previously (Lin et al., 2003; Tao et al., 2015). All experiments were performed under normoxic conditions at room temperature (∼22–24°C).
Cytosolic Ca2+ Imaging
Intracellular Ca2+ of basilar arterial myocytes was measured in situ with a Ca2+ sensitive fluorescent dye (Fluo-4 AM, Cat No F14201, Invitrogen, Carlsbad, CA, United States) using a Zeiss 710 NLO laser scanning confocal imaging work station (Thornwood, NY, United States) with an inverted microscope (ZEISS Axio Observer), using procedures based on previous studies (Hadley et al., 2012; Harraz et al., 2014; Shen et al., 2018). Fluo-4 AM was dissolved in DMSO creating a 1 mM stock solution. Arteries were placed in BSS and exposed to a Fluo-4 concentration of 10 μM with 0.1% Pluronic F-127 Cat No P6867 (Invitrogen) from a 20% w/v stock solution in DMSO for 1 h in the dark at room temperature. These arterial segments were subsequently washed with BSS for 30 min to facilitate dye esterification and were then cut into linear strips for testing. These strips were pinned to Sylgard en face (Ellsworth Adhesives, Germantown, WI, United States) with fine insect dissecting pins and placed into an open bath imaging chamber (Warner Instruments, Hamden, CT, United States). Myocytes were illuminated at 488 nm via a Krypton-argon laser. Emitted light was captured with a photomultiplier tube with a band limited spectral grating of range 493–622 nm in both full frame and line scan imaging studies.
Whole-Cell Ca2+ Recordings
A time series of 300 full frame Fluo-4 fluorescence images of 512 × 512 pixels were made over 234 s (roughly 780 ms/frame). In order to make sure that the smooth muscle intracellular Ca2+ was recorded, the pinhole was set at an imaging depth of 5.4 μm, which is roughly the depth of an individual smooth muscle cell based on examination of both fixed and live preparations (Hadley et al., 2012; Shen et al., 2018). This thicker optical sectioning of samples was performed to mitigate the effects of sample ruffling, allowing for visualization of more myocytes than would otherwise be possible at other depths. The sample was focused just below the internal elastic lamina to center on the myocytes and avoid the significant autofluorescence of this layer when excited at 488 nm. Images were taken using 12-bit sampling and a water immersion 63X Plan Apochromat, 1.4 NA objective. Under basal conditions a time series was made to assess the Ca2+ oscillatory behavior of individual myocytes. After the video recording 30–50 line scans of 18.9 s were collected as detailed below, which were used to measure Ca2+ spark activity. Careful effort was made to not duplicate imaging regions for the line scans, which ensured proper sampling of individual cells in the arterial wall and reduced the potential for photobleaching and laser induced toxicity. Following basal recordings, BSS with 30 mM K (30K) through equimolar replacement of potassium for sodium was added to the tissues to depolarize the plasma membrane of the myocytes (Hadley et al., 2012; Harraz et al., 2014; Shen et al., 2018). 5 min later we performed another time series recording, followed by series of line scan recordings.
Whole-Cell Ca2+ Analysis
Regions of interest for full frame fluorescent imaging were automatically detected using the LCPro plug in for Fiji (ImageJ) (Francis et al., 2012; Schindelin et al., 2012), following 8 bit image conversion, image registration with StackReg using rigid body settings and cropping the images to 490 × 490 pixels to remove any erroneous pixels due to image registration (Thévenaz et al., 1998; Francis et al., 2012; Shen et al., 2018). LCPro was used to calculate fractional fluorescence intensity for automatically detected regions of interest, based on user defined region sizes, which were circles of 1.56 μm (six pixels). This region size was chosen as it provides analysis within individual cells as opposed to larger regions, which often overlap with adjacent cells.
The program performs analysis of the fluorescent intensity increases for each of the detected regions of interest above statistical noise, which was set to a threshold of P < 0.05. The program then analyzes the signal in the time and space domains for several parameters. In the current study, examinations were made of the following parameters: Amplitude (F/F0), which is the maximum amplitude of the signal transient as defined by the global maximum of the F/F0 curve within the time-period of the signal. Duration (s), which is the time interval defined by the period from one-half of the maximal amplitude values before and after the signal peak. Rise Time (s), which is written as “attack” in the program output and is the time interval defined by the period from one-half of the maximal signal amplitude value preceding the signal peak to the signal peak. Decay (s), which is the time interval defined by the time at the maximum amplitude to the time at one-half of the maximal amplitude following the peak. Area under curve (AUC) (F*s/F0), which is the discrete right-sided Riemann sum of the signal area during the period that is above one-half of the maximum signal amplitude (Francis et al., 2012, 2014, 2016; Shen et al., 2018).
Within the recordings from basilar arterial myocytes several different types of calcium events were observed that could be discriminated based on their temporal signaling characteristics, which were also observed in our previous studies performed in pulmonary arteries (Shen et al., 2018). These include rapid events that were between 0 and 4 s and that were scored as calcium sparks, medium duration oscillations that were defined as being between 4 and 40 s, and finally long duration oscillations that were defined both by being 40 + s and by their unique plateau at their peak fluorescence (Shen et al., 2018). Visual analysis of each data record was used to clarify the divisions between groups and to remove false positives. The number of cells with Ca2+ oscillations were based on visual examination of myocytes in 1,002 μm2 regions of interest with replicates performed in three separate regions per recording for each animal (Hadley et al., 2012; Shen et al., 2018).
Spatiotemporal Ca2+ Signaling Analysis
Spatiotemporal characteristics of the whole-cell Ca2+ oscillatory events were analyzed by performing a cross correlation analysis in time and space for the regions of interest detected by LCPro (Shen et al., 2018). The program is utilized to gain a deeper understanding of local and distant signaling networks within smooth muscle cells as well as to better determine the methods of Ca2+ trafficking between cells. This program combines a custom set of physiology analysis tools written within python that enable us to compare and observe correlative Ca2+ oscillations within a data set. A correlation coefficient of r = 0.8 was used for all data sets based on parameter space searches for correlating oscillations utilizing r values ranging from r = 0.1–1.0 (Shen et al., 2018). “Friends” were determined to be calcium transients that had a correlation coefficient greater than r = 0.8. “Neighbors” were defined as calcium transients occurring in nearby myocytes. A radius of 100 pixels (∼23 μm) was chosen as the critical spatial cutoff because this distance allowed for the inclusion of most all neighboring cells.
Line Scan Ca2+ Recordings and Analysis
Fluo-4 AM loaded basilar arterial myocytes were recorded with a Zeiss LSM 710 NLO laser scanning confocal imaging workstation on an inverted platform (Zeiss Axio Observer Z1). After being loaded with Fluo-4, the arteries were prepared and imaged as detailed above. Line scan images were captured at 529 lines per second and recordings were made for a duration of 18.9 s. The lateral pixel size ranged from 0.0148 to 0.0911 μm per pixel and the pinhole was adjusted so that the cells were imaged to a depth of 2.5 μm, which is about 50% of the width of a myocyte based on prior morphological studies in live cells (Hadley et al., 2012; Shen et al., 2018).
Line scan recordings were analyzed for the percentage of cells with Ca2+ sparks, their frequency, amplitude, and other spatiotemporal characteristics via SparkLab 4.3.1 (Harraz et al., 2015; Shen et al., 2018). Before analysis, background fluorescence was subtracted from each recording assuming a homogenous level of background noise, and the fluorescence recordings normalized (Harraz et al., 2015; Shen et al., 2018). The threshold for spark detection was set to 2 standard deviations above this background level of fluorescence. 3–4 animals were used per group for both subcellular Ca2+ spark and whole-cell Ca2+ oscillation analysis. Refer to the figure legends for specific n values associated with each set of data.
Experimental Protocol
Fluo-4 loaded basilar arteries placed en-face in the imaging chamber were handled akin to our previous imaging studies of cerebral, pulmonary and uterine arteries (Hashad et al., 2017; Shen et al., 2018; Hu et al., 2020). Once a region of interest was identified in the control solution a video recording was made followed by line-scan images of individual myocytes in the arterial wall. These recordings required roughly 1 h. Upon completion of the video and line-scan recordings, the imaging chamber was removed from the microscope workstation and placed on the platform of a stereomicroscope. The control bathing solution was then gently removed with a pipette, and replaced with a 30 mM K bathing solution. To ensure that the bathing solution was fully exchanged the procedure was performed a total of three times. The tissues were then allowed to recover for a minimum of 5 min before resuming imaging to ensure that the cytosolic calcium reached a dynamic equilibrium.
Drugs and Chemical Reagents
Unless otherwise indicated, all reagents were purchased from Sigma-Aldrich (St. Louis, MO, United States).
Statistical Methods
Data analysis and the production of graphs was performed using GraphPad Prism 9.1.2 (La Jolla, CA, United States). The graphs present the data as means ± SD. Continuous variable data was tested for normality prior to analysis. These datasets were not normally distributed and as such they were evaluated using non-parametric statistical tests. A Kruskal–Wallis ANOVA with a Dunn’s multiple comparisons test was performed on non-parametric datasets when making examinations between the experimental groups based on animal age, altitude, and treatment condition. The density of cell firing in video recordings was evaluated by drawing 3 boxes of a known size in each recording and counting the number of cells with events inside the box, with each cell only being counted once. For calcium spark analysis the numbers of cells were evaluated for those cells with or without events. The number of cells in videos or lines-scan recordings that had events was summarized. Statistical comparisons of this discrete data was made by performing contingency analysis between groups using Chi-Square analysis to evaluate potential changes in the incidence of Ca2+ oscillatory and spark events in whole-cell and line-scan recordings. A P-value of P < 0.05 was considered significant and were further broken down to P < 0.01 and P < 0.001 where appropriate.
The specific test performed for each data set are depicted in the figure legends. Sample sizes were determined by several different measurements including the number of animals studied, the number of regions of interest showing Ca2+ oscillations, the number of line scans examined for Ca2+ sparks, and number of Ca2+ events. The percentage of cells firing with Ca2+ sparks was determined by the number of line scans containing Ca2+ events through visual observation relative to the total number of line scans observed.
Results
Whole-Cell Ca2+ Oscillations in Basilar Arterial Myocytes
The first series of studies were designed to test the hypotheses that membrane depolarization would increase Ca2+ oscillatory activity and that long-term hypoxia (LTH) would impair oscillatory events in basilar arterial myocytes. Figure 1A shows a maximum intensity projection of Fluo-4 fluorescence in myocytes of the basilar arterial wall from a time series recording. A video corresponding to this image is provided in Supplementary Material Video 1. The image illustrates that the myocytes are spindle shaped and closely associated with one another and consistent with myocytes from other vascular beds and species that we have examined, including pulmonary and uterine arteries of sheep as well as cerebral arteries of rat and mesenteric arteries of mouse (Hadley et al., 2012; Harraz et al., 2014, 2015; Shen et al., 2018; Hu et al., 2020). Figure 1B shows the Fluo-4 fluorescence over time recorded from two regions of interest (ROI’s) in individual myocytes that were automatically detected with LCPro (Francis et al., 2012; Shen et al., 2018), a custom analysis program and subsequently analyzed for spatial and temporal aspects to the Ca2+ signals. These “medium duration” events were the most common in the whole-cell Ca2+ recordings, with the oscillations spreading relatively uniformly through the myocytes and lasting between 4 and 40 s.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/760176/fphys-12-760176-HTML-r1/image_m/fphys-12-760176-g001.jpg)
Figure 1. Representative Ca2+ responses in basilar arterial myocytes from an adult sheep recorded en face under 30K conditions. (A) maximum intensity projection for Fluo-4 fluorescence of recorded cells using laser scanning confocal microscopy. Arrows point to regions of interest in individual myocytes shown in (B); fluorescence intensity tracing showing spontaneous Ca2+ oscillations in two ROIs. Scale bar (white) = 10 μm.
The recordings of Fluo-4 fluorescence (Figure 2) show that there were a variety of different Ca2+ oscillatory events in addition to the medium duration events described in Figure 1, which were akin to those we identified in pulmonary arterial myocytes (Shen et al., 2018). Results from these imaging studies are presented in Figures 2–7. Figure 2 shows that some Ca2+ events were very rapid, and had features that were reminiscent of Ca2+ sparks in that they were localized in subcellular regions (Figure 2Aa) and were fast, being roughly 1–2 s in duration (Figure 2Ab), being at the edge of detectability with the full-frame recording techniques used for whole-cell Ca2+ recordings (Shen et al., 2018). The corresponding Supplementary Video 2 recording for Figure 2Aa illustrates that there are multiple types of Ca2+ events in the arterial myocytes. The inability to record rapid events with enough fidelity to perform analysis using full-frame recording techniques made at 1.28 Hz is exemplified in the exploded view of a single Ca2+ spark event. These rapid Ca2+ signaling events were scored as Ca2+ sparks, which are classically due to coordinated activation of ryanodine receptors on the sarcoplasmic endoplasmic reticulum (SER) that release Ca2+ into the cytosol (Jaggar et al., 1998a,b; Hadley et al., 2012; Hashad et al., 2017; Shen et al., 2018). These data also illustrate the need to use high-speed line-scan recordings to examine the Ca2+ spark events more thoroughly (Shen et al., 2018), results of which are presented in Figures 8–10.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/760176/fphys-12-760176-HTML-r1/image_m/fphys-12-760176-g002.jpg)
Figure 2. Basilar arterial myocytes show additional distinct forms of Ca2+ signals. The figure shows (a) maximum intensity projections with arrows pointing to regions of interest that correspond to the fluorescence intensity tracings shown in (b). This includes (A) short duration Ca2+ “spark” events (0–4 s) and (B) long duration Ca2+ oscillations (40+s). The labeled time-series tracings of F/F0 for these two types of Ca2+ responses were recorded using Fluo-4 and then detected and evaluated post hoc with LCPro. Scale bars (white) in (Aa and Ba) = 10 μm.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/760176/fphys-12-760176-HTML-r1/image_m/fphys-12-760176-g003.jpg)
Figure 3. Long-term hypoxia reduces the percentage of basilar arterial myocytes with depolarization mediated Ca2+ oscillations in fetal and adult sheep. Each bar represents the percentage of cells with Ca2+ oscillations under control (clear) or with treatment of 30K (black) based on an examination of myocytes in 1,002 μm2 regions of interest. Replicates were performed in 3 separate regions per recording for fetal normoxic (3 animals along with 79 control and 71 30K myocytes), fetal hypoxic (4 animals with 91 control and 74 30K myocytes), adult normoxic (3 animals with 69 control and 68 30K myocytes), and adult hypoxic (4 animals with 61 control and 68 30K myocytes) animals. *P < 0.05, **P < 0.01, ***P < 0.001 indicate significance based on a chi-square test.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/760176/fphys-12-760176-HTML-r1/image_m/fphys-12-760176-g004.jpg)
Figure 4. Long-term hypoxia has minimal effect on the magnitude or temporal aspects of medium duration Ca2+ oscillations in fetal and adult sheep. Effects of membrane depolarization with 30 mM K, long-term hypoxia, and animal age on (A) area under the curve, (B) amplitude of the fractional fluorescence, (C) duration of the event, (D) rise time, and (E) decay time in cytosolic fractional fluorescence. Bars represent mean ± SD for each parameter; closed circles specify control conditions, and open circles treatment with 30K. *P < 0.05, **P < 0.01, ***P < 0.001 indicate significance based on a Kruskal–Wallis ANOVA with a Dunn’s multiple comparisons test based on ranks. Control recordings were made in 207/3 FN, 201/4 FH, 225/3 AN, and 325/4 AH ROIs and animals, respectively. 30K recordings were made in 642/3 FN, 394/4 FH, 1204/3 AN, and 312/4 AH ROIs and animals, respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/760176/fphys-12-760176-HTML-r1/image_m/fphys-12-760176-g005.jpg)
Figure 5. Long-term hypoxia has minimal effect on the magnitude or temporal aspects of long duration Ca2+ oscillations in fetal and adult sheep. Effects of artificial membrane depolarization with 30 mM K, long-term hypoxia, and animal age on (A) area under the curve, (B) amplitude of the fractional fluorescence, (C) duration of the fluorescent event, (D) rise time for fractional fluorescence in the cytosol and (E) decay time for cytosolic fractional fluorescence. Bars represent mean ± SD for each parameter; closed circles specify control conditions, and open circles treatment with 30K. *P < 0.05, **P < 0.01, ***P < 0.001 indicate significance based on a Kruskal–Wallis ANOVA with a Dunn’s multiple comparisons test based on ranks. Control recordings were made in 116/3 FN, 28/4 FH, 61/3 AN, and 48/4 AH ROIs and animals, respectively. 30K recordings were made in 52/3 FN, 5/4 FH, 8/3 AN, and 16/4 AH ROIs and animals, respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/760176/fphys-12-760176-HTML-r1/image_m/fphys-12-760176-g006.jpg)
Figure 6. Spatial and temporal correlations of Ca2+ responses in basilar arterial myocytes from an adult sheep recorded en face under basal conditions. (A) Maximum intensity projection for Fluo-4 fluorescence of recorded cells using laser scanning confocal microscopy. Highly correlated events are colored and plotted around the center (reference) region of interest (ROI, white dot). Gray dots show ROIs of spontaneous Ca2+ oscillations that showed less than an 80% temporal correlation with the white ROI. The large white open circle shows a ∼23 μm boundary that was used for deriving spatial correlation measurements shown for the spatial and temporal analysis in Figure 7. (B) Fluorescence intensity tracing showing spontaneous Ca2+ oscillations, with the black line being the average tracing and the orange lines being the individual ROIs. ROIs were plotted using correlation coefficients with greater than 80% temporal correlation plotted surrounding the center ROI (white dot in A). Key shows the degree of correlation of each ROI relative to the reference ROI.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/760176/fphys-12-760176-HTML-r1/image_m/fphys-12-760176-g007.jpg)
Figure 7. Long-term hypoxia mitigates the influence of 30K mediated membrane depolarization on friends and neighbors. (A–E) Arteries were analyzed for number of correlated ROIs (friends), distance between correlated ROIs, percentage of nearby ROIs (neighbors) that were correlated, and percentage of correlated ROIs that were in nearby cells. Bars represent mean ± SD for each condition; closed circles specify control conditions, and open circles treatment with 30K. *P < 0.05, **P < 0.01, ***P < 0.001 indicate significance based on based on a Kruskal–Wallis ANOVA with a Dunn’s multiple comparisons test based on ranks. Control responses were obtained from 178/3 AN, 235/4 AH, 175/3 FN, and 163/4 FH ROIs and animals, respectively. 30K responses were obtained from 844/3 AN, 244/4 AH, 537/3 FN, and 303/4 FH ROIs and animals, respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/760176/fphys-12-760176-HTML-r1/image_m/fphys-12-760176-g008.jpg)
Figure 8. Representative Ca2+ spark tracings for basilar arterial myocytes from fetal and adult sheep. (A–D) Fluo-4 fluorescence tracings recorded from basilar arterial myocyte line scans of fetal and adult sheep under normoxic and hypoxic conditions analyzed with SparkLab 4.3.1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/760176/fphys-12-760176-HTML-r1/image_m/fphys-12-760176-g009.jpg)
Figure 9. The incidence of Ca2+ sparks activated by membrane depolarization is reduced in fetuses and by hypoxia in adults. (A) percentage of cells with or without Ca2+ sparks and (B) Ca2+ spark firing frequency. Bars represent mean ± SD for each condition. Clear bars and closed circles specify control conditions, while black bars and open circles treatment with 30 K. Data in A were analyzed by a chi-square test while data in (B) were analyzed by Kruskal–Wallis one-way ANOVA with Dunn’s multiple comparison test based on ranks for each group *P < 0.05, **P < 0.01, ***P < 0.001. Recordings were made in FN control (161/3) and 30K (172/3), FH control (171/4) and 30K (185/4), AN control (150/3) and 30K (155/3), and AH control (207/4) and 30K (217/4) lines (cells)/animals, respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/760176/fphys-12-760176-HTML-r1/image_m/fphys-12-760176-g010.jpg)
Figure 10. Magnitude and kinetics of Ca2+ sparks were minimally influenced by animal age, LTH, or membrane depolarization. (A–D) amplitude, full width at half-maximum, full duration at half-maximum and tau exposed to control (closed circles) or 30 K (open circles) for Ca2+ spark events of arterial myocytes from adult and fetal sheep under normoxic and hypoxic conditions. Bars represent means ± SD for each parameter. Data were analyzed by a Kruskal–Wallis one-way ANOVA with Dunn’s multiple comparison test based on ranks for each group **P < 0.01, ***P < 0.001. Recordings were made in FN control (41/161/3) and 30K (199/172/3), FH control (30/171/4) and 30K (136/185/4), AN control (173/150/3) and 30K (365/155/3), and AH control (52/207/4) and 30K (177/217/4) Ca2+ sparks/lines/animals, respectively.
Infrequently, we observed a third category of “long-duration” Ca2+ oscillations with large spatial spread through the myocytes (Figure 2Ba) and that lasted over 40 s. Figure 2Bb shows that these events had a unique quality in that they had rapid increases in cytosolic Ca2+ that plateaued and were sustained for an extended period followed by a prompt relaxation of the Ca2+ back to basal levels. The Supplementary Material Video 3 provides the recording associated with Figure 2B.
For analysis purposes, the oscillations were broken into three discrete categories: short-duration Ca2+ sparks, which lasted less than 4 s, medium-duration Ca2+ oscillations of 4–40 s, and long-duration Ca2+ oscillations that lasted more than 40 s (Shen et al., 2018). Acutely, the subcellular Ca2+ spark events are likely to regulate the activity of voltage and calcium activated potassium channels that give rise to spontaneous transient outward currents and subsequent vasodilatory responses (Jaggar et al., 1998b; Hu et al., 2011, 2012; Hashad et al., 2017). The whole-cell Ca2+ oscillation events in comparison are likely responsible for causing arterial contraction as well as regulating a variety of other cellular processes including metabolism or transcription (Wilson et al., 2005; Goyal et al., 2011; Papamatheakis et al., 2011; Ureña et al., 2013; Tykocki et al., 2017).
The percentage of cells with medium-duration Ca2+ oscillation events were impacted by animal age and LTH, with activity being enhanced by treatment with 30 mM potassium (30K), data that is summarized in Figure 3. The potassium in the extracellular bathing solution was raised to 30 mM in order to depolarize the plasma membrane and increase activation of voltage gated calcium channels. Increased calcium influx across the plasma membrane leads to more calcium in the junctional space between the plasma membrane and SER. The elevated junctional calcium subsequently activates clusters of ryanodine receptors on the SER, which gives rise to subcellular events (Ca2+ sparks, see Figures 8–10) and when combined with other sources of calcium there are pronounced calcium oscillations that spread through the myocyte (Jaggar et al., 1998b; Janiak et al., 2001; Tykocki et al., 2017).
The percentage of cells with calcium events was similar in the fetal and adult groups and unaffected by LTH in the control bathing solution with low (5 mM) extracellular potassium. Membrane depolarization with 30K increased calcium oscillations in normoxic fetal and adult basilar arterial myocytes, an effect that was more pronounced in adult as compared to fetal sheep. Following long term hypoxia, however, membrane depolarization failed to increase cellular calcium oscillations in both fetal and adult basilar arterial preparations.
The kinetics and magnitude of the medium-duration Ca2+ oscillations were then quantified, with summarized data shown in Figure 4. Overall, animal age, LTH, and membrane depolarization had only mild effects on the Ca2+ oscillations. Figure 4A shows that under basal conditions the area under the curve (AUC) of Ca2+ oscillations were unaffected by maturation. However, in adult normoxic myocytes the AUC was slightly increased relative to that in the adult hypoxic period. Membrane depolarization caused a modest increase in the AUC in adult hypoxic myocytes. The AUC was unaffected by membrane depolarization in fetal normoxic and hypoxic groups; however, the AUC was lower in myocytes from fetal hypoxic and adult normoxic relative to fetal normoxic animals following membrane depolarization. The amplitudes of the Ca2+ responses are shown in Figure 4B. Event amplitudes recorded from myocytes of the adult hypoxic and fetal normoxic groups were reduced relative to those from normoxic adults under basal conditions. Membrane depolarization increased Ca2+ event amplitude in the fetal normoxic and adult hypoxic groups. Event amplitude was also elevated in adult hypoxic relative to adult normoxic myocytes following membrane depolarization. The duration of the Ca2+ events is shown in Figure 4C and again there were only modest differences between the groups. There was no impact of animal age on event duration under basal conditions, however, LTH caused a shortening of the oscillatory duration in adult myocytes. Membrane depolarization had no direct influence on the oscillatory duration in any of the four groups, though fetal normoxic myocytes had longer durations as compared to events from fetal hypoxic and adult normoxic myocytes. Figure 4D shows that the Ca2+ rise time during the oscillations was unaffected by LTH or animal age under basal conditions. Membrane depolarization shortened the rise time modestly in only the fetal LTH group. Following membrane depolarization, the rise time was longer in fetal normoxic myocytes as compared to fetal hypoxic and adult normoxic myocytes. Figure 4E shows the decay time for the Ca2+ signal. Under basal conditions, LTH modestly shortened the time for the decay of the oscillatory signal in both fetal and adult myocytes. Membrane depolarization shortened the decay in myocytes from adult normoxic animals, which caused the decay duration to become equivalent with that from LTH adult animals. Following membrane depolarization, the time for Ca2+ decay was also shorter in LTH fetal myocytes relative to the normoxic counterparts.
Long-Duration Ca2+ Oscillations
Long-duration Ca2+ oscillations were significantly less common than medium duration oscillations or Ca2+ sparks as illustrated by the reduced number of events. The effects of animal age, LTH, and membrane depolarization on these events are presented in Figure 5. The most notable influences were observed in the event amplitudes. Adult normoxic animals under basal conditions had increased event amplitudes (Figure 5B), decreased rise times (Figure 5D), and increased decay times relative to normoxic fetuses (Figure 5E). Figure 5B also shows that long-term hypoxia increased event amplitude in fetal animals but decreased the amplitude in adults under control conditions. Membrane depolarization decreased event amplitude in fetal hypoxic and adult normoxic animals but increased the amplitude in adult hypoxic animals.
Spatiotemporal Ca2+ Signaling
The interactions of the Ca2+ oscillations among smooth muscle cells were then analyzed for temporal correlations between the Ca2+ transients, or “friendship,” between oscillations at the various ROIs. The larger the correlation coefficient, the greater the level of event “friendship.” Temporally related correlations between the Ca2+ events were then identified by those events that had a correlation coefficient of r ≥ 0.8 (Shen et al., 2018). Results from one representative sample are provided in Figure 6. Figure 6A shows an image projection with overlaid regions of interest that have varied levels of correlations, while Figure 6B shows the average fluorescence (black line) along with the fluorescence intensity values for individual regions of interest that were correlated. Spatial correlations between events are shown in Figure 7B, a relationship characterized as “neighbors” as delineated by those ROIs that were within a 100 pixels (∼23 μm) radius of any individual ROI, as delineated by the larger white circle on the image in Figure 6B (Shen et al., 2018). Animal age, LTH and membrane depolarization had diverse effects on the numbers of friends, with summary data provided in Figure 7. Figure 7 shows that relative to fetuses there was a decrease in the number of friends in adult arterial myocytes under control conditions. LTH decreased the number of friends in fetal arterial myocytes but increased them in adults under control conditions. Membrane depolarization increased the number of friends in adult normoxic myocytes. As illustrated in Figure 7B, the number of neighbors was also influenced by animal age, LTH, and membrane depolarization. Under basal conditions, there was a decrease in the number of neighbors in adults relative to fetuses. Under basal conditions, LTH increased the number of neighbors in adults but this was unaffected in fetuses. Membrane depolarization increased the number of neighbors in normoxic fetuses and adults, but LTH blunted the responses in myocytes of both age groups. Figure 7C delineates the distance of temporally correlated events (friends). The distance between friends was shorter in normoxic adults relative to fetuses under control and depolarized conditions. Long-term hypoxia resulted in a shortening of the distance only in fetuses under control conditions. Figure 7D provides the percentage of friends who are neighbors, illustrating the spatial and temporal relationships between Ca2+ oscillations. Overall, a majority of temporally correlated oscillatory events (friends) are closely associated (neighbors), however, there were some subtle effects of animal age, LTH, and membrane depolarization. There was a slight increase in the percentage of adult normoxic myocyte ROIs who were neighbors relative to those in fetal normoxic myocytes. Although there were fewer Ca2+ oscillatory events in fetal hypoxic myocytes than their normoxic counterparts, a greater proportion were neighbors under control conditions. Membrane depolarization preferentially increased the percentage of friends who were neighbors in normoxic myocytes independent of animal age. Figure 7E presents the percentage of neighbors who were friends, which provides an index of the temporal relationship between Ca2+ oscillatory events that are spatially related. Unlike the tight spatial relationship between temporally related events, events that occur in nearby regions are not necessarily temporally correlated. Even still, there were influences of LTH and membrane depolarization. Long-term hypoxia preferentially reduced the percentage of neighbors who were friends in fetal myocytes under control conditions. Membrane depolarization preferentially reduced the percentage of neighbors who were friends in myocytes of fetal and adult normoxic animals.
Activation of Ca2+ Sparks
The last series of studies examined the impact of membrane depolarization and hypoxia on the activation of Ca2+ sparks in myocytes from fetuses and adults. Figure 8 shows representative Ca2+ spark tracings from line scan recordings for control groups analyzed with SparkLab (Figures 8A–D; Shen et al., 2018). Our lab and others have shown that membrane depolarization with extracellular K+ can enhance Ca2+ spark activity and vascular contractility, which is why we chose to depolarize the membrane with 30 mM K+ (Jaggar et al., 1998a; Hadley et al., 2012; Papamatheakis et al., 2012; Harraz et al., 2014; Shen et al., 2018). Membrane depolarization with 30K increased the prevalence of myocytes with Ca2+ spark activity in all groups and the frequency of activation as shown in Figures 9A,B. Under control conditions the prevalence and frequency of Ca2+ sparks in myocytes from normoxic adults were greater than in the fetal period. LTH did not impact the prevalence or frequency of Ca2+ sparks in fetuses, but reduced spark activity in adult myocytes to fetal levels. Membrane depolarization increased Ca2+ spark activation in adult animals while LTH preferentially blunted spark activation in adults. Overall, these findings illustrate that LTH impairs spark activity in basilar arteries but support the general influence of membrane depolarization on Ca2+ spark activity that we have shown in our previous studies on pulmonary and uterine arteries from sheep as well as in middle cerebral as well as mesenteric arteries from rats and mice (Hadley et al., 2012; Harraz et al., 2014, 2015; Shen et al., 2018; Hu et al., 2019).
The magnitude and kinetics of Ca2+ sparks were quantified in basilar arterial myocytes with the summary results shown in Figures 10A–D. Overall, there was little impact of animal age, LTH, or membrane depolarization on the quality of the Ca2+ sparks. When comparing events in myocytes from fetuses and adults, the full width of half maximum was narrowed slightly in normoxic control conditions. Long-term hypoxia increased the full width of half maximum of Ca2+ sparks in adult myocytes under control conditions. Membrane depolarization modestly increased the full width of half maximum of Ca2+ sparks in adult normoxic myocytes and there was no impact on the exponential decay time constant (Tau) of the Ca2+ spark events.
Discussion
The current study was designed to examine the impact of animal age and high altitude induced LTH on cytosolic Ca2+ oscillations and sparks in basilar arterial myocytes. These studies were performed because of the integral role basilar arteries have in regulating vascular reactivity and providing blood flow to the brain stem, which is critical to autonomic cardiorespiratory functions that are impacted by LTH. Basilar arterial myocytes had robust Ca2+ oscillations and the number of event sites were increased by membrane depolarization. The data also illustrate that animal age and LTH modify the coupling of membrane depolarization to the activation of whole-cell Ca2+ oscillations as well as the generation of localized Ca2+ sparks.
Whole-Cell Ca2+ Oscillations
Restriction of whole-cell intracellular Ca2+ oscillations by LTH is not unique to the cerebral vasculature. Previously, we reported that Ca2+ oscillatory activity was reduced by LTH in pulmonary arterial myocytes (Hadley et al., 2012; Shen et al., 2018). The parallel changes in Ca2+ signals in basilar and pulmonary arterial myocytes due to LTH suggest there may be commonalities in the mechanisms underlying the changes in oscillatory activity following LTH. Myocyte cytosolic Ca2+ oscillations are dependent on a functional SER (Mufti et al., 2010) and InsP3 receptor activation is important to the generation of oscillatory signals (Adebiyi et al., 2010, 2011; Mufti et al., 2015; Tykocki et al., 2017). These relationships lead to the potential that losses in Ca2+ oscillatory activity and other temporal signaling aspects by LTH may be coupled to impairments in InsP3 signaling. Indeed, previous work from our group shows that LTH reduces InsP3 receptor expression in both fetal and adult middle cerebral arteries (Longo et al., 1996; Ueno et al., 1997). Depression in InsP3 signaling may therefore also underlie the suppression in membrane depolarization induced Ca2+ oscillations we report here. With regards to the current findings that LTH restricts cellular Ca2+ responses, our previous data illustrate that the SER Ca2+ stores of sheep pulmonary arteries are largely intact regardless of age or hypoxic stress (Hadley et al., 2012), although there is significant SER stress as evidenced by structural changes in the SER of pulmonary arterial myocytes as well as increases in markers of SER stress (Leslie et al., 2021). Such evidence suggests that the major aberrations in Ca2+ responses are more likely due to dysregulation of pathways critical to regulating Ca2+ signaling that may involve SER stress as opposed to fundamental disruption of the storage capabilities of the intracellular Ca2+ stores.
Membrane depolarization dependent increases in Ca2+ oscillatory activity was expected as this is common in smooth muscle (Iino, 1990; Hadley et al., 2012; Shen et al., 2018). Increases in oscillatory frequency with membrane depolarization may be the result of changes in the speed of Ca2+ release and uptake at the SER as well as increase Ca2+ flux across the plasma membrane (Keizer et al., 1995; Wilson et al., 2002, 2005; Goyal et al., 2009). The increases in the oscillatory amplitude of cells from fetal animals suggests that membrane depolarization may increase oscillatory signals either through increased activation of InsP3 or ryanodine receptors on the SER or through enhanced recruitment and coupling to L-Type Ca2+ channels (Jaggar et al., 1998b; Blood et al., 2002; del Corsso et al., 2006; Ureña et al., 2013) or other Ca2+ permeable plasma membrane ion channels (Jaggar et al., 1998b; Janiak et al., 2001; Earley et al., 2005; Adebiyi et al., 2010, 2011; Dahan et al., 2012; Ureña et al., 2013; Harraz et al., 2014). In this regard, the picture that emerges is that the impact of LTH on basilar arterial Ca2+ signaling is complex and mediated through modification of multiple pathways important to the regulation of intracellular Ca2+.
The finding that whole-cell Ca2+ oscillations were reduced in basilar arterial myocytes from normoxic fetuses compared to adults builds from our previous evidence of age-related changes in Ca2+ signaling of cerebral arteries (Ducsay et al., 2018). One possibility is that these differences are coupled to developmental related effects on the role of SER and Ca2+ influx pathways to contractility. Indeed, contractions of cerebral arteries from fetal sheep have reduced reliance on release of intracellular Ca2+ stores and increased dependence on Ca2+ influx as compared to adults (Long et al., 2000). These changes in arterial contractions were shown to be coupled to greater L-type Ca2+ channel protein expression in fetus and responsiveness to L-type Ca2+ channel activation (Long et al., 2000; Blood et al., 2002). Even still, Ca2+ oscillations are due to sequential filling and release of internal Ca2+ stores in smooth muscle, with modulation through extracellular Ca2+ influx pathways (Janiak et al., 2001; Wilson et al., 2005; Hume et al., 2009). Given that we did not fully assess the role of extracellular Ca2+ entry or release of intracellular Ca2+ stores to the oscillatory activity in basilar arterial myocytes, we still do not know what role age has to the process of Ca2+ handling.
The influence of membrane depolarization on the temporal and spatial aspects to the Ca2+ oscillations in normoxic animals are compelling. When the tissue was depolarized, there were significantly more events during the recordings in the arterial myocytes as demonstrated by the uptick in spatially related events (neighbors), however, there were fewer temporally related events (friends). On closer examination the neighboring events mainly occurred within individual cells as opposed to adjacent cells. While we have not explored this deeply in the current study, the findings suggest that when myocytes are depolarized the Ca2+ oscillations do not propagate from cell to cell, but rather, the number of events increase through recruitment of unrelated cells across the arterial wall. One explanation for this apparent lack of cell-to-cell communication is that membrane depolarization recruits additional L-Type Ca2+ channels on individual cells in a stochastic manner, which then increases Ca2+ responses in separate cells. Further, the data suggest gap junction connections between myocytes are not being activated, which limits cell-to-cell propagation of Ca2+ signals (Hald et al., 2014; Welsh et al., 2018; Zechariah et al., 2020). Such findings are not unfounded as physiological studies and computational models suggest there is poor electromechanical coupling between arterial myocytes but propagation is spread more easily through the vascular endothelium (Hald et al., 2014; Tykocki et al., 2017; Welsh et al., 2018).
Ca2+ Spark Activity
Previous evidence illustrates that the localized rapid release of calcium, coined Ca2+ sparks, in cerebral vascular myocytes is due to coordinated activation of ryanodine receptor clusters (Jaggar et al., 1998a,b). These Ca2+ spark events are important as they are coupled to vasodilation through activation of large conductance potassium channels (BKCa). Depolarization of the plasma membrane is well regarded to enhance Ca2+ spark activity and the ensuing repolarization of the membrane and vasodilation acts as a negative feedback regulator of vasoreactivity. With regards to the effects of age and LTH, basilar arterial myocytes from adults had far greater Ca2+ spark activity relative to those from fetuses. Furthermore, LTH reduced spark activity in adults to fetal levels and suppressed the depolarization mediated increase in Ca2+ spark activity in adult myocytes. The inability of depolarization to increase Ca2+ spark activity in adults following LTH is interesting and reminiscent of the impact of LTH on fetal pulmonary arterial myocytes (Hadley et al., 2012; Shen et al., 2018). Much like the impact on whole-cell Ca2+ oscillations, this parallel effect of LTH on basilar and pulmonary arterial myocytes suggests there may be a common underlying mechanism whereby LTH suppresses depolarization induced activation of Ca2+ sparks. Presumably, this would be due to a loss in communication between plasma membrane channels and ryanodine receptors on the SER, possibly involving SER stress such as we have shown in fetal pulmonary arterial myocytes (Leslie et al., 2021). These losses in communication potentially include L-type Ca2+ channels but may also involve disruption in T-type Ca2+ channels or TRPV4, which are also important to Ca2+ spark activation in vascular myocytes (Earley et al., 2005; Harraz et al., 2014, 2015).
The impact of age and long-term hypoxia on the role of the Ca2+ spark events to feedback activation of BKCa is not fully resolved. Previous evidence from our group, however, illustrates that BKCa from fetal basilar arteries have increased activity in response to changes in the cytosolic Ca2+ (Lin et al., 2003) mediated through enhanced protein kinase G (Lin et al., 2005, 2006). Further, LTH increases membrane expression of BKCa and Ca2+ dependent channel activity in both fetal and adult basilar arterial myocytes (Tao et al., 2015). Together these data suggest that even though fetal basilar myocytes have reduced Ca2+ spark activity relative to adults there may be preservation of BKCa channel activity and vasodilatory capacity through altered channel regulation. Similarly, vasodilation may be maintained following LTH in adults even though there is suppression in Ca2+ spark activity through enhanced Ca2+ dependent activation of BKCa. Unraveling the full impacts of age and LTH on the regulation of basilar arterial reactivity and the influence this has on brain blood flow will require further interrogation.
Perspectives
The impact of long-term hypoxia on cerebral vascular function remains poorly understood, especially with regards to the differing influences on fetuses and adults. The current data begin to shed light on the complexity of the effects of LTH on cellular Ca2+ signals depending on animal age that are important to vascular reactivity. The adaptations in Ca2+ signaling to LTH reported here are potentially juxtaposed by modifications in the regulation of BKCa channels, which are critical components in the feedback regulation of vascular reactivity and brain blood flow. These effects lead us to speculate that LTH-induced dysregulation in vascular function prompts compensatory responses that act to preserve brainstem nutrient and oxygen delivery.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The animal study was reviewed and approved by the Loma Linda University Health Institutional Animal Care and Use Committee. All study procedures adhered to the Animal Welfare Act, the National Institutes of Health Guide for the Care and Use of Laboratory Animals (https://grants.nih.gov/grants/olaw/Guide-for-the-Care-and-use-of-laboratory-animals.pdf), “The Guiding Principles in the Care and Use of Animals” approved by the Council of the American Physiological Society.
Author Contributions
SW, CW, AB, and LZ contributed to conception and design of the study. MR and NO performed laboratory experimentation. CW and JLP designed customized software for data analysis. SBC, MR, NO, and CR performed data and statistical analysis. CR wrote the first draft of the manuscript. SC wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
This work was supported by the National Science Foundation under Major Research Instrumentation, Division of Biological Infrastructure Grant No. 0923559 (SW) and National Institutes of Health Grants HD083132 and HL149608 (LZ).
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.
Acknowledgments
We would like to thank Shane Glasgow for his role in the work and Dr. Michael Francis for clarifying the LCPro signal analysis algorithms. A portion of this material was performed in the Loma Linda University School of Medicine Advanced Imaging and Microscopy Core that is supported by the Loma Linda University School of Medicine.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2021.760176/full#supplementary-material
References
Adebiyi, A., Narayanan, D., and Jaggar, J. H. (2011). Caveolin-1 assembles type 1 inositol 1,4,5-trisphosphate receptors and canonical transient receptor potential 3 channels into a functional signaling complex in arterial smooth muscle cells. J. Biol. Chem. 286, 4341–4348. doi: 10.1074/jbc.M110.179747
Adebiyi, A., Zhao, G., Narayanan, D., Thomas-Gatewood, C. M., Bannister, J. P., and Jaggar, J. H. (2010). Isoform-selective physical coupling of TRPC3 channels to IP3 receptors in smooth muscle cells regulates arterial contractility. Circ. Res. 106, 1603–1612. doi: 10.1161/CIRCRESAHA.110.216804
Adeoye, O. O., Silpanisong, J., Williams, J. M., and Pearce, W. J. (2015). Role of the sympathetic autonomic nervous system in hypoxic remodeling of the fetal cerebral vasculature. J. Cardiovasc. Pharmacol. 65, 308–316. doi: 10.1097/FJC.0000000000000192
Ainslie, P. N., and Subudhi, A. W. (2014). Cerebral blood flow at high altitude. High Alt. Med. Biol. 15, 133–140. doi: 10.1089/ham.2013.1138
Blood, A. B., Zhao, Y., Long, W., Zhang, L., and Longo, L. D. (2002). L-type Ca2+ channels in fetal and adult ovine cerebral arteries. Am. J. Physiol. Regul. Integr. Comp. Physiol. 282, R131–R138. doi: 10.1152/ajpregu.00318.2001
Brenner, R., Peréz, G. J., Bonev, A. D., Eckman, D. M., Kosek, J. C., Wiler, S. W., et al. (2000). Vasoregulation by the beta1 subunit of the calcium-activated potassium channel. Nature 407, 870–876. doi: 10.1038/35038011
Brugniaux, J. V., Hodges, A. N., Hanly, P. J., and Poulin, M. J. (2007). Cerebrovascular responses to altitude. Respir. Physiol. Neurobiol. 158, 212–223. doi: 10.1016/j.resp.2007.04.008
Dahan, D., Ducret, T., Quignard, J. F., Marthan, R., Savineau, J. P., and Estève, E. (2012). Implication of the ryanodine receptor in TRPV4-induced calcium response in pulmonary arterial smooth muscle cells from normoxic and chronically hypoxic rats. Am. J. Physiol. Lung Cell. Mol. Physiol. 303, L824–L833. doi: 10.1152/ajplung.00244.2011
del Corsso, C., Ostrovskaya, O., Mcallister, C. E., Murray, K., Hatton, W. J., Gurney, A. M., et al. (2006). Effects of aging on Ca2+ signaling in murine mesenteric arterial smooth muscle cells. Mech. Ageing Dev. 127, 315–323. doi: 10.1016/j.mad.2005.12.001
Ducsay, C. A., Goyal, R., Pearce, W. J., Wilson, S., Hu, X.-Q., and Zhang, L. (2018). Gestational hypoxia and developmental plasticity. Physiol. Rev. 98, 1241–1334. doi: 10.1152/physrev.00043.2017
Earley, S., Heppner, T. J., Nelson, M. T., and Brayden, J. E. (2005). TRPV4 forms a novel Ca2+ signaling complex with ryanodine receptors and BKCa channels. Circ. Res. 97, 1270–1279. doi: 10.1161/01.RES.0000194321.60300.d6
Ferriero, D. M. (2004). Neonatal brain injury. N. Engl. J. Med. 351, 1985–1995. doi: 10.1056/NEJMra041996
Francis, M., Qian, X., Charbel, C., Ledoux, J., Parker, J. C., and Taylor, M. S. (2012). Automated region of interest analysis of dynamic Ca2 + signals in image sequences. Am. J. Physiol. Cell Physiol. 303, C236–C243. doi: 10.1152/ajpcell.00016.2012
Francis, M., Waldrup, J. R., Qian, X., Solodushko, V., Meriwether, J., and Taylor, M. S. (2016). Functional tuning of intrinsic endothelial Ca2+ dynamics in swine coronary arteries. Circ. Res. 118, 1078–1090. doi: 10.1161/CIRCRESAHA.115.308141
Francis, M., Waldrup, J., Qian, X., and Taylor, M. S. (2014). Automated analysis of dynamic Ca2+ signals in image sequences. J. Vis. Exp. 88:e51560. doi: 10.3791/51560
Gebremedhin, D., Bonnet, P., Greene, A. S., England, S. K., Rusch, N. J., Lombard, J. H., et al. (1994). Hypoxia increases the activity of Ca(2+)-sensitive K+ channels in cat cerebral arterial muscle cell membranes. Pflugers Arch. 428, 621–630. doi: 10.1007/BF00374586
Giussani, D. A. (2016). The fetal brain sparing response to hypoxia: physiological mechanisms. J. Physiol. 594, 1215–1230. doi: 10.1113/JP271099
Goyal, R., Angermann, J. E., Ostrovskaya, O., Buchholz, J. N., Smith, G. D., and Wilson, S. M. (2009). Enhanced capacitative calcium entry and sarcoplasmic-reticulum calcium storage capacity with advanced age in murine mesenteric arterial smooth muscle cells. Exp. Gerontol. 44, 201–207. doi: 10.1016/j.exger.2008.10.007
Goyal, R., Papamatheakis, D. G., Loftin, M., Vrancken, K., Dawson, A. S., Osman, N. J., et al. (2011). Long-term maternal hypoxia: the role of extracellular Ca2+ entry during serotonin-mediated contractility in fetal ovine pulmonary arteries. Reprod. Sci. 18, 948–962. doi: 10.1177/1933719111401660
Hadley, S. R., Blood, Q., Rubalcava, M., Waskel, E., Lumbard, B., Le, P., et al. (2012). Maternal high-altitude hypoxia and suppression of ryanodine receptor-mediated Ca2+ sparks in fetal sheep pulmonary arterial myocytes. Am. J. Physiol. Lung Cell. Mol. Physiol. 303, L799–L813. doi: 10.1152/ajplung.00009.2012
Hald, B. O., Jacobsen, J. C., Sandow, S. L., Holstein-Rathlou, N. H., and Welsh, D. G. (2014). Less is more: minimal expression of myoendothelial gap junctions optimizes cell-cell communication in virtual arterioles. J. Physiol. 592, 3243–3255. doi: 10.1113/jphysiol.2014.272815
Harraz, O. F., Abd El-Rahman, R. R., Bigdely-Shamloo, K., Wilson, S. M., Brett, S. E., Romero, M., et al. (2014). Ca(V)3.2 channels and the induction of negative feedback in cerebral arteries. Circ. Res. 115, 650–661. doi: 10.1161/CIRCRESAHA.114.304056
Harraz, O. F., Brett, S. E., Zechariah, A., Romero, M., Puglisi, J. L., Wilson, S. M., et al. (2015). Genetic ablation of CaV3.2 channels enhances the arterial myogenic response by modulating the RyR-BKCa axis. Arterioscler. Thromb. Vasc. Biol. 35, 1843–1851. doi: 10.1161/ATVBAHA.115.305736
Hashad, A. M., Mazumdar, N., Romero, M., Nygren, A., Bigdely-Shamloo, K., Harraz, O. F., et al. (2017). Interplay among distinct Ca2+ conductances drives Ca2+ sparks/spontaneous transient outward currents in rat cerebral arteries. J. Physiol. 595, 1111–1126. doi: 10.1113/JP273329
Hu, X. Q., Song, R., Romero, M., Dasgupta, C., Huang, X., Holguin, M. A., et al. (2019). Pregnancy increases Ca(2+) sparks/spontaneous transient outward currents and reduces uterine arterial myogenic tone. Hypertension 73, 691–702. doi: 10.1161/HYPERTENSIONAHA.118.12484
Hu, X. Q., Song, R., Romero, M., Dasgupta, C., Min, J., Hatcher, D., et al. (2020). Gestational hypoxia inhibits pregnancy-induced upregulation of Ca(2+) sparks and spontaneous transient outward currents in uterine arteries via heightened endoplasmic reticulum/oxidative stress. Hypertension 76, 930–942. doi: 10.1161/HYPERTENSIONAHA.120.15235
Hu, X. Q., Xiao, D., Zhu, R., Huang, X., Yang, S., Wilson, S. M., et al. (2012). Chronic hypoxia suppresses pregnancy-induced upregulation of large-conductance Ca2+-activated K+ channel activity in uterine arteries. Hypertension 60, 214–222. doi: 10.1161/HYPERTENSIONAHA.112.196097
Hu, X. Q., Xiao, D., Zhu, R., Huang, X., Yang, S., Wilson, S., et al. (2011). Pregnancy upregulates large-conductance Ca(2+)-activated K(+) channel activity and attenuates myogenic tone in uterine arteries. Hypertension 58, 1132–1139. doi: 10.1161/HYPERTENSIONAHA.111.179952
Hume, J. R., Mcallister, C. E., and Wilson, S. M. (2009). Caffeine inhibits InsP3 responses and capacitative calcium entry in canine pulmonary arterial smooth muscle cells. Vascul. Pharmacol. 50, 89–97. doi: 10.1016/j.vph.2008.11.001
Hunter, C. J., Blood, A. B., White, C. R., Pearce, W. J., and Power, G. G. (2003). Role of nitric oxide in hypoxic cerebral vasodilatation in the ovine fetus. J. Physiol. 549, 625–633. doi: 10.1113/jphysiol.2002.038034
Iino, M. (1990). Biphasic Ca2+ dependence of inositol 1,4,5-trisphosphate-induced Ca release in smooth muscle cells of the guinea pig taenia caeci. J. Gen. Physiol. 95, 1103–1122. doi: 10.1085/jgp.95.6.1103
Jaggar, J. H., Stevenson, A. S., and Nelson, M. T. (1998a). Voltage dependence of Ca2+ sparks in intact cerebral arteries. Am. J. Physiol. 274, C1755–C1761. doi: 10.1152/ajpcell.1998.274.6.C1755
Jaggar, J. H., Wellman, G. C., Heppner, T. J., Porter, V. A., Perez, G. J., Gollasch, M., et al. (1998b). Ca2+ channels, ryanodine receptors and Ca(2+)-activated K+ channels: a functional unit for regulating arterial tone. Acta Physiol. Scand. 164, 577–587. doi: 10.1046/j.1365-201X.1998.00462.x
Janiak, R., Wilson, S. M., Montague, S., and Hume, J. R. (2001). Heterogeneity of calcium stores and elementary release events in canine pulmonary arterial smooth muscle cells. Am. J. Physiol. Cell Physiol. 280, C22–C33. doi: 10.1152/ajpcell.2001.280.1.C22
Jansen, G. F., Kagenaar, D. A., Basnyat, B., and Odoom, J. A. (2002). Basilar artery blood flow velocity and the ventilatory response to acute hypoxia in mountaineers. Respir. Physiol. Neurobiol. 133, 65–74. doi: 10.1016/S1569-9048(02)00152-0
Kamitomo, M., Alonso, J. G., Okai, T., Longo, L. D., and Gilbert, R. D. (1993). Effects of long-term, high-altitude hypoxemia on ovine fetal cardiac output and blood flow distribution. Am. J. Obstetr. Gynecol. 169, 701–707. doi: 10.1016/0002-9378(93)90646-Z
Keizer, J., Li, Y. X., Stojilkoviæ, S., and Rinzel, J. (1995). InsP3-induced Ca2+ excitability of the endoplasmic reticulum. Mol. Biol. Cell 6, 945–951. doi: 10.1091/mbc.6.8.945
Koenigs, M., and Grafman, J. (2009). The functional neuroanatomy of depression: distinct roles for ventromedial and dorsolateral prefrontal cortex. Behav. Brain Res. 201, 239–243. doi: 10.1016/j.bbr.2009.03.004
Ledoux, J., Werner, M. E., Brayden, J. E., and Nelson, M. T. (2006). Calcium-activated potassium channels and the regulation of vascular tone. Physiology (Bethesda) 21, 69–78. doi: 10.1152/physiol.00040.2005
Lee, S. J., Hatran, D. P., Tomimatsu, T., Pena, J. P., Mcauley, G., and Longo, L. D. (2009). Fetal cerebral blood flow, electrocorticographic activity, and oxygenation: responses to acute hypoxia. J. Physiol. 587, 2033–2047. doi: 10.1113/jphysiol.2008.166983
Leslie, E., Lopez, V., Anti, N. A. O., Alvarez, R., Kafeero, I., Welsh, D. G., et al. (2021). Gestational long-term hypoxia induces metabolomic reprogramming and phenotypic transformations in fetal sheep pulmonary arteries. Am. J. Physiol. Lung Cell. Mol. Physiol. 320, L770–L784. doi: 10.1152/ajplung.00469.2020
Lin, M. T., Hessinger, D. A., Pearce, W. J., and Longo, L. D. (2003). Developmental differences in Ca2+-activated K+ channel activity in ovine basilar artery. Am. J. Physiol. Heart Circ. Physiol. 285, H701–H709. doi: 10.1152/ajpheart.00138.2003
Lin, M. T., Hessinger, D. A., Pearce, W. J., and Longo, L. D. (2006). Modulation of BK channel calcium affinity by differential phosphorylation in developing ovine basilar artery myocytes. Am. J. Physiol. Heart Circ. Physiol. 291, H732–H740. doi: 10.1152/ajpheart.01357.2005
Lin, M. T., Longo, L. D., Pearce, W. J., and Hessinger, D. A. (2005). Ca2+-activated K+ channel-associated phosphatase and kinase activities during development. Am. J. Physiol. Heart Circ. Physiol. 289, H414–H425. doi: 10.1152/ajpheart.01079.2004
Long, W., Zhang, L., and Longo, L. D. (2000). Cerebral artery sarcoplasmic reticulum Ca(2+) stores and contractility: changes with development. Am. J. Physiol. Regul. Integr. Comp. Physiol. 279, R860–R873. doi: 10.1152/ajpregu.2000.279.3.R860
Longo, L. D., Ueno, N., Zhao, Y., Pearce, W. J., and Zhang, L. (1996). Developmental changes in alpha 1-adrenergic receptors, IP3 responses, and NE-induced contraction in cerebral arteries. Am. J. Physiol. 271, H2313–H2319. doi: 10.1152/ajpheart.1996.271.6.H2313
Mufti, R. E., Brett, S. E., Tran, C. H., Abd El-Rahman, R., Anfinogenova, Y., El-Yazbi, A., et al. (2010). Intravascular pressure augments cerebral arterial constriction by inducing voltage-insensitive Ca2+ waves. J. Physiol. 588, 3983–4005. doi: 10.1113/jphysiol.2010.193300
Mufti, R. E., Zechariah, A., Sancho, M., Mazumdar, N., Brett, S. E., and Welsh, D. G. (2015). Implications of αvβ3 integrin signaling in the regulation of Ca2+ waves and myogenic tone in cerebral arteries. Arterioscler. Thromb. Vasc. Biol. 35, 2571–2578. doi: 10.1161/ATVBAHA.115.305619
Nelson, K. B. (2003). Can we prevent cerebral palsy? N. Engl. J. Med. 349, 1765–1769. doi: 10.1056/NEJMsb035364
Papamatheakis, D. G., Patel, J. J., Blood, Q., Merritt, T. T., Longo, L. D., and Wilson, S. M. (2012). Depolarization-dependent contraction increase after birth and preservation following long-term hypoxia in sheep pulmonary arteries. Pulm Circ. 2, 41–53. doi: 10.4103/2045-8932.94832
Papamatheakis, D. G., Vemulakonda, S., Blood, Q., Goyal, R., Rubalcava, M., Vrancken, K., et al. (2011). Preservation of serotonin-mediated contractility in adult sheep pulmonary arteries following long-term high-altitude hypoxia. High Alt. Med. Biol. 12, 253–264. doi: 10.1089/ham.2010.1076
Pearce, W. (2006). Hypoxic regulation of the fetal cerebral circulation. J. Appl. Physiol. 100, 731–738. doi: 10.1152/japplphysiol.00990.2005
Pena, J. P., Tomimatsu, T., Hatran, D. P., Mcgill, L. L., and Longo, L. D. (2007). Cerebral blood flow and oxygenation in ovine fetus: responses to superimposed hypoxia at both low and high altitude. J. Physiol. 578, 359–370. doi: 10.1113/jphysiol.2006.119925
Roher, A. E., Debbins, J. P., Malek-Ahmadi, M., Chen, K., Pipe, J. G., Maze, S., et al. (2012). Cerebral blood flow in Alzheimer’s disease. Vasc. Health Risk Manag. 8, 599–611. doi: 10.2147/VHRM.S34874
Schindelin, J., Arganda-Carreras, I., Frise, E., Kaynig, V., Longair, M., Pietzsch, T., et al. (2012). Fiji: an open-source platform for biological-image analysis. Nat. Methods 9, 676–682. doi: 10.1038/nmeth.2019
Shen, C. P., Romero, M., Brunelle, A., Wolfe, C., Dobyns, A., Francis, M., et al. (2018). Long-term high-altitude hypoxia influences pulmonary arterial L-type calcium channel-mediated Ca(2+) signals and contraction in fetal and adult sheep. Am. J. Physiol. Regul. Integr. Comp. Physiol. 314, R433–R446. doi: 10.1152/ajpregu.00154.2017
Stoll, B. J., Hansen, N. I., Bell, E. F., Shankaran, S., Laptook, A. R., Walsh, M. C., et al. (2010). Neonatal outcomes of extremely preterm infants from the NICHD Neonatal Research Network. Pediatrics 126, 443–456. doi: 10.1542/peds.2009-2959
Tao, X., Lin, M. T., Thorington, G. U., Wilson, S. M., Longo, L. D., and Hessinger, D. A. (2015). Long-term hypoxia increases calcium affinity of BK channels in ovine fetal and adult cerebral artery smooth muscle. Am. J. Physiol. Heart Circ. Physiol. 308, H707–H722. doi: 10.1152/ajpheart.00564.2014
Thévenaz, P., Ruttimann, U. E., and Unser, M. (1998). A pyramid approach to subpixel registration based on intensity. IEEE Trans. Image Process. 7, 27–41. doi: 10.1109/83.650848
Thorpe, R. B., Hubbell, M. C., Silpanisong, J., Williams, J. M., and Pearce, W. J. (2017). Chronic hypoxia attenuates the vasodilator efficacy of protein kinase G in fetal and adult ovine cerebral arteries. Am. J. Physiol. Heart Circ. Physiol. 313, H207–H219. doi: 10.1152/ajpheart.00480.2016
Tweed, A., Cote, J., Lou, H., Gregory, G., and Wade, J. (1986). Impairment of cerebral blood flow autoregulation in the newborn lamb by hypoxia. Pediatr. Res. 20, 516–519. doi: 10.1203/00006450-198606000-00007
Tweed, W. A., Cote, J., Pash, M., and Lou, H. (1983). Arterial oxygenation determines autoregulation of cerebral blood flow in the fetal lamb. Pediatr. Res. 17, 246–249. doi: 10.1203/00006450-198304000-00002
Tykocki, N. R., Boerman, E. M., and Jackson, W. F. (2017). Smooth muscle ion channels and regulation of vascular tone in resistance arteries and arterioles. Compr. Physiol. 7, 485–581. doi: 10.1002/cphy.c160011
Ueno, N., Zhao, Y., Zhang, L., and Longo, L. D. (1997). High altitude-induced changes in alpha1-adrenergic receptors and Ins(1,4,5)P3 responses in cerebral arteries. Am. J. Physiol. 272, R669–R674. doi: 10.1152/ajpregu.1997.272.2.R669
Ureña, J., Fernández-Tenorio, M., Porras-González, C., González-Rodríguez, P., Castellano, A., and López-Barneo, J. (2013). A new metabotropic role for L-type Ca(2+) channels in vascular smooth muscle contraction. Curr. Vasc. Pharmacol. 11, 490–496. doi: 10.2174/1570161111311040012
Volpe, J. J. (1997). Brain injury in the premature infant. Neuropathology, clinical aspects, pathogenesis, and prevention. Clin. Perinatol. 24, 567–587. doi: 10.1016/S0095-5108(18)30159-3
Welsh, D. G., Tran, C. H. T., Hald, B. O., and Sancho, M. (2018). The conducted vasomotor response: function, biophysical basis, and pharmacological control. Annu. Rev. Pharmacol. Toxicol. 58, 391–410. doi: 10.1146/annurev-pharmtox-010617-052623
Wilson, S. M., Mason, H. S., Ng, L. C., Montague, S., Johnston, L., Nicholson, N., et al. (2005). Role of basal extracellular Ca2+ entry during 5-HT-induced vasoconstriction of canine pulmonary arteries. Br. J. Pharmacol. 144, 252–264. doi: 10.1038/sj.bjp.0706077
Wilson, S. M., Mason, H. S., Smith, G. D., Nicholson, N., Johnston, L., Janiak, R., et al. (2002). Comparative capacitative calcium entry mechanisms in canine pulmonary and renal arterial smooth muscle cells. J. Physiol. 543, 917–931. doi: 10.1113/jphysiol.2002.021998
Yang, Q., Tong, X., Schieb, L., Vaughan, A., Gillespie, C., Wiltz, J. L., et al. (2017). Vital Signs: Recent Trends in Stroke Death Rates - United States, 2000-2015. MMWR Morb. Mortal Wkly. Rep. 66, 933–939. doi: 10.15585/mmwr.mm6635e1
Zechariah, A., Tran, C. H. T., Hald, B. O., Sandow, S. L., Sancho, M., Kim, M. S. M., et al. (2020). Intercellular conduction optimizes arterial network function and conserves blood flow homeostasis during cerebrovascular challenges. Arterioscler. Thromb. Vasc. Biol. 40, 733–750. doi: 10.1161/ATVBAHA.119.313391
Keywords: arterial myocytes, calcium oscillations in living cells, calcium sparks, ontogeny, high altitude
Citation: Reid C, Romero M, Chang SB, Osman N, Puglisi JL, Wilson CG, Blood AB, Zhang L and Wilson SM (2022) Long-Term Hypoxia Negatively Influences Ca2+ Signaling in Basilar Arterial Myocytes of Fetal and Adult Sheep. Front. Physiol. 12:760176. doi: 10.3389/fphys.2021.760176
Received: 17 August 2021; Accepted: 25 October 2021;
Published: 18 January 2022.
Edited by:
Dieter Blottner, Charité – University Medicine Berlin, GermanyReviewed by:
Josh Maxwell, Wake Forest School of Medicine, United StatesJoseph P. Y. Kao, University of Maryland, Baltimore, United States
Michael D. Stern, National Institute on Aging, National Institutes of Health (NIH), United States
Jianjie Ma, The Ohio State University, United States
Copyright © 2022 Reid, Romero, Chang, Osman, Puglisi, Wilson, Blood, Zhang and Wilson. 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: Sean M. Wilson, seanwilson@llu.edu