Skip to main content

ORIGINAL RESEARCH article

Front. Integr. Neurosci., 16 December 2021
This article is part of the Research Topic Bridging the Gap: Implementing Building Blocks of Nervous System to Simulate Neuronal Functions at Different Levels View all 5 articles

Silencing of Activity During Hypoxia Improves Functional Outcomes in Motor Neuron Networks in vitro

  • 1Department of Neuromedicine and Movement Science, Faculty of Medicine and Health Sciences, Norwegian University of Science and Technology, Trondheim, Norway
  • 2Department of Neurology, St. Olav’s Hospital, Trondheim University Hospital, Trondheim, Norway
  • 3Department of Pharmacology and Clinical Neurosciences, Division of Neuro, Head, and Neck, Umeå University Hospital, Umeå, Sweden
  • 4Department of Community Medicine and Rehabilitation, Umeå University, Umeå, Sweden

The effects of hypoxia, or reduced oxygen supply, to brain tissue can be disastrous, leading to extensive loss of function. Deoxygenated tissue becomes unable to maintain healthy metabolism, which leads to increased production of reactive oxygen species (ROS) and loss of calcium homoeostasis, with damaging downstream effects. Neurons are a highly energy demanding cell type, and as such they are highly sensitive to reductions in oxygenation and some types of neurons such as motor neurons are even more susceptible to hypoxic damage. In addition to the immediate deleterious effects hypoxia can have on neurons, there can be delayed effects which lead to increased risk of developing neurodegenerative diseases such as amyotrophic lateral sclerosis (ALS), even if no immediate consequences are apparent. Furthermore, impairment of the function of various hypoxia-responsive factors has been shown to increase the risk of developing several neurodegenerative disorders. Longitudinal assessment of electrophysiological network activity is underutilised in assessing the effects of hypoxia on neurons and how their activity and communication change over time following a hypoxic challenge. This study utilised multielectrode arrays and motor neuron networks to study the response to hypoxia and the subsequent development of the neuronal activity over time, as well as the effect of silencing network activity during the hypoxic challenge. We found that motor neuron networks exposed to hypoxic challenge exhibited a delayed fluctuation in multiple network activity parameters compared to normoxic networks. Silencing of activity during the hypoxic challenge leads to maintained bursting activity, suggesting that functional outcomes are better maintained in these networks and that there are activity-dependent mechanisms involved in the network damage following hypoxia.

Introduction

Hypoxia is a condition of reduced oxygen supply to brain tissue that can occur in a variety of circumstances, from low atmospheric oxygen to reduced blood flow. It is distinguished from anoxia, a more adverse condition of no oxygen supply, which is not considered in this study. Low oxygen levels can have severe adverse effects on affected tissue, including enhanced production of reactive oxygen species (ROS), impaired ATP production and loss of calcium homoeostasis (Kristian and Siesjo, 1998), all of which can lead to tissue damage and cell death. Neuronal tissue has a very high energy consumption, which makes it particularly vulnerable to hypoxic insult (Hossmann, 1999), for example following a stroke or respiratory arrest. This vulnerability is not limited to the immediate and drastic brain damage often associated with sustained reduction of oxygen supply to the brain. It can manifest in more subtle ways, notably increased risk of developing neurodegenerative disorders (Sun et al., 2006; Guglielmotto et al., 2009; Bhatia et al., 2017) as well as aggravate pathological processes in patients with existing neurodegenerative disease and animal models of neurodegenerative disease (Zhang et al., 2019). Compounding the link between hypoxia and neurodegenerative disorders is also a growing body of evidence that genetic traits which impair proper function of various hypoxia-responsive elements can increase the risk of developing neurodegenerative disorders, including Parkinson’s disease (PD) and Amyotrophic lateral sclerosis (ALS) (Greenway et al., 2006; van Es et al., 2011; Zou et al., 2012).

Hypoxia causes neurotoxic effects which may be the consequence of different mechanisms, both activity-dependent and activity-independent. Changes in neuronal excitability are widely reported in both dissociated neuronal cell cultures and brain slice preparations, albeit with contradictory findings. For example, some studies have found that hypoxia leads to neuronal hyperexcitability in Schaffer collaterals of rat hippocampal slices (Schiff and Somjen, 1985) and CA1 neurons in rat hippocampal slices (Godukhin et al., 2002), while other studies have reported decreased excitability in CA1 neurons in mouse hippocampal slices (Gu and Haddad, 2001) and dissociated mouse hippocampal neurons (Gavello et al., 2012). There may also be cell type-specific vulnerabilities which contribute to this deviance in findings, such as the reported vulnerability of inhibitory neurons (Turovsky et al., 2013). The influx of calcium into hypoxic neurons has been widely implicated in the neurotoxic process (Kristian and Siesjo, 1998), and NMDA-receptors appear to be an important avenue of Ca2+ influx due to their proximity to mitochondria, which act as intracellular calcium stores (Starkov et al., 2004). Blocking NMDA-receptors has been shown to reduce the formation of ROS (Lafon-Cazal et al., 1993). Activation of GABA-A receptors has also been shown to ameliorate oxidative damage during hypoxia due to Cl influx in rodents (Zhao et al., 2003; Kumar and Goyal, 2008), and inhibition of activity of hippocampal neurons with leptin also leads to improved outcomes following hypoxia (Gavello et al., 2012). Together these results suggest an activity-dependent mechanism for neuronal damage following hypoxia.

Although neuronal tissue is susceptible to hypoxia in general, there are notable differences in how different neuronal populations respond to hypoxia, and some populations are clearly more vulnerable than others. However, most of the research on neuronal susceptibility examines the vulnerable hippocampal neurons, while other vulnerable groups such as motor neurons are less frequently studied. However, not only are motor neurons particularly susceptible to hypoxia (Watzlawik et al., 2015), but repeated exposure to hypoxic conditions also appears to greatly increase the risk of developing motor neuron disorders like ALS (Vanacore et al., 2010). Motor neurons are highly energy demanding, and findings indicate that motor neurons with ALS pathology exhibit reduced mitochondrial respiration and oxidative ATP production, with a shift toward increased glycolysis (Hor et al., 2021) which can contribute to overall energy deficits which contribute to disease progression (Dupuis et al., 2011). This mirrors the effects of hypoxia, where mitochondrial activation of hypoxia inducible factor 1α (HIF-1α) in response to reduced concentration of oxygen leads to inhibition of oxidative phosphorylation and increased production of ATP through glycolysis, resulting in an overall energy deficit compared to normoxic conditions (Solaini et al., 2010). This makes the study of the responses of motor neurons to hypoxic challenge also relevant to studies of ALS, both through shared pathways and as a potential environmental contributor to the disease.

Studies of neuronal activity are underutilised in the study of hypoxia exposure, especially with regard to the effects on broader network dynamics. This study uses in vitro motor neuron networks on multielectrode arrays (MEAs) to examine responses to hypoxic insult in a longitudinal perspective. It considers the effects of network activity during the episode in maintaining healthy network function following hypoxia.

Materials and Methods

In vitro Neural Networks

Human motor neuron progenitors derived from a healthy donor (ax0078) were purchased from Axol Bioscience. The cells were seeded on two Axion Biosystems 48-well Cytoview MEA plates (M768-tMEA-48B) and two 8-well chambered slides (Ibidi 80841). Prior to seeding, all culture wells were coated with 0.05% polyethylenimine (PEI) diluted in HEPES (both Sigma-Aldrich). MEA wells were each coated with 70 μL while slide wells were coated with 250 μL, before being left overnight in an incubator at 37°C and 5% CO2. Water reservoirs in the MEA plates were filled with distilled water to reduce evaporation. The next day, the PEI was removed, and the wells were rinsed four times with distilled water, before they were left to air dry at room temperature overnight. The wells were then coated with natural mouse laminin (Thermo Fisher Scientific), at 20 μg/mL, diluted in PBS. MEA wells were each coated with 70 μL while slide wells were coated with 250 μL. Cells were plated immediately after removing the laminin, without washing. Each MEA well was seeded with 25,000 motor neuron progenitors and 2,500 human astrocytes (Thermo Fisher Scientific Gibco K1884, lot 1948466) in mixed suspension, using Axol Motor Neuron Recovery media (ax0071) supplemented with retinoic acid at 0.1 μM and Y-27623 2HCL at 10 μM according to the supplier’s instructions. Each slide well was seeded with 50,000 motor neuron progenitors and 5,000 human astrocytes in the same media. The ratio of neurons to astrocytes was based on recommendations from Axion Biosystems to ensure trophic support for the neuron progenitors without interfering with electrophysiological recordings. The MEAs and slides were left in the incubator for 1 h before adding more media to reach a final volume of 300 μL in each MEA well and 400 μL in each slide well. The next day, half of the cell media volume was replaced with Axol Motor Neuron Recovery media supplemented with retinoic acid, Y-27623 2HCL and 1% Pen-Strep. The day after, and every other day for the duration of the experiment, half of the media volume was replaced with Axol Motor Neuron Maintenance media (ax0072) supplemented with 10 ng/mL CNTF, 5 ng/mL BDNF, 0.5 μM retinoic acid and 1% Pen Strep. A timeline for coating, seeding and media composition is shown in Figure 1A.

FIGURE 1
www.frontiersin.org

Figure 1. Experiment timeline. (A) The experimental timeline. Cells were seeded on day 0 after coating with polyethyleminine (PEI) and laminin. RI: Rock inhibitor Y-27623 2HCL, PS, Penicillin-Streptomycin antibiotics; ICC, Immunocytochemistry. (B) Detailed timeline of hypoxic challenge. Following day 56 activity was recorded every day, and media was half replaced every other day. (C) Network bursts were identified by bins of high firing rate, with the threshold indicated by the red horizontal line.

Electrophysiological Recordings

Recordings of electrical activity were made using an Axion Maestro acquisition tool, with temperature of 37°C and 5% CO2, recording 1 h of baseline activity at a sampling rate of 12.5 kHz. Each MEA well recorded from 16 electrodes. After being placed in the Axion Maestro, neural networks were left for 1 h to allow activity to stabilise before recordings were started. The first recording was made at 49 DIV and subsequent recordings were made every other day until the neural networks reached 53 DIV, and daily after this. Data analysis was based on Axion spike detection using a dynamic threshold detection of 7 standard deviations. Spike times were exported via Neuroexplorer 5 to Matlab R2019b and analysed using custom scripts, see section “Data Analysis”. Due to a technical fault in the acquisition process, the recording of activity at 58 DIV was not stored for the hypoxic plate (data point omitted), and the recording of the non-hypoxic plate at 58 DIV only stored 23 min of the 1 h planned recording.

Immunocytochemistry and Imaging

To confirm motor neuron identity as well as expression of relevant neurotransmitters, parallel neural networks derived from the same progenitors were tested by immunocytochemistry (ICC) at 49 DIV for the markers in Table 1, as well as nuclear staining with Hoechst (bisBenzimide H 33342 trihydrochloride, 14533, Sigma-Aldrich, 1:10,000). Prior to staining, cells were fixed in 4% Paraformaldehyde in PBS, before blocking was carried out with 5% Goat serum and 0.6% Triton-X in PBS. Primary and secondary antibody staining was done using 2.5% Goat serum and 0.3% Triton-X in PBS at the above indicated antibody dilutions.

TABLE 1
www.frontiersin.org

Table 1. Immunocytochemistry (ICC): Markers, antibodies and concentrations.

Images were acquired with an EVOS M5000 microscope (Invitrogen Thermo Fisher Scientific) with DAPI (AMEP4650), CY5 (AMEP4656), GFP (AMEP4651), and TxRed (AMEP4655) LED light cubes and an Olympus UPLSAP020x lens, 20x/0.75 NA (N1480500). All image processing was done in Fiji/ImageJ.

Hypoxic Challenge

Two multiwell MEA plates; hypoxic condition and normoxic control, plated with 24 and 21 motor neuron networks of healthy motor neurons, respectively, aged 54 DIV were used. Within the hypoxic group and the normoxic control, half the motor neuron networks were inhibited with GABA (12 networks in the hypoxic group and 11 in the normoxic control group) while the remaining networks were uninhibited (12 networks in the hypoxic group and 10 in the normoxic control group), prior to exposure to hypoxic challenge at 1% O2 and 5% CO2 for 1 h, during which activity was recorded. 30 μL of 550 μM GABA was added to the 300 μL of cell media in each culture for a final GABA concentration of 50 μM. Uninhibited motor neuron networks were only supplemented with 30 μL of cell media with no GABA. The experiment therefore had four groups in total: inhibited hypoxic (n = 12), uninhibited hypoxic (n = 12), inhibited normoxic (n = 11), and uninhibited normoxic (n = 10). Prior to the experiment, baseline activity was recorded for all neural networks for 1 h. After the hypoxia/control exposure, half media changes were made to all neural networks. Starting the day after the hypoxic challenge, activity was recorded for 1 h daily for 10 days, and half of the media volume was replenished every other day. The last recordings were made at 64 DIV. A more detailed timeline of the hypoxic challenge is shown in Figure 1B.

Data Analysis

Network bursts were identified by examining a binned spike distribution using 50 ms bins. Every bin that exceeded a threshold of firing rate was identified as potentially part of a burst. The threshold used was mean plus 5 standard deviations, indicated in red in Figure 1C. Then, the number of firing electrodes in each of these bins was checked, and bins with high firing rate but fewer than 20% of active electrodes in the recording fired were considered non-bursting activity. Finally, subsequent time bins with bursting were combined into longer bursts. From this burst train, identifying burst beginnings and ends as well as size in terms of participating electrodes, further measures were derived, such as the proportion of spikes which occurred in network bursts, hereafter referred to as burst propensity. Network synchrony was measured by the coherence index.

Because the data was not normally distributed, all summary data is shown as median ± median absolute deviation and comparisons between groups were assessed by Wilcoxon rank sum test, with Bonferroni correction for multiple comparisons. Relative dispersion of the data was assessed by an alternative to the coefficient of variation based on interquartile range, RCVQ (Arachchige et al., 2020).

Results

Immunocytochemistry

The neural networks were imaged live throughout development until 49 DIV, and matured from a dissociated, even layer of cells into progressively more aggregated clusters connected by fasciculated neurites. After maturation at 49 DIV, immunostaining confirmed expression of GABA and glutamate receptors as well as heavy neurofilament, a mature cytoskeletal marker (Figure 2C) and motor neuron markers Islet 1 (Figure 2A), HB9 and ChAT (Figure 2B).

FIGURE 2
www.frontiersin.org

Figure 2. Immunocytochemistry. After maturation, neurons expressed motor neuron specific markers as well as relevant neurotransmitter receptors and markers of mature cytoskeleton. (A) Overlap of expression of motor neuron specific marker Islet1 and neuronal marker NeuN, alongside heavy neurofilament, indicates the presence of mature motor neurons. Scale bar 50 μm. (B) The co-expression of motor neuron markers HB9 and ChAT further confirms motor neuron identity. Scale bar 50 μm. (C) Expression of receptors for GABA and glutamate confirms the capacity for excitatory signalling within the motor networks, as well as the susceptibility of the networks to the GABA-inhibition. Scale bar 100 μm.

Neural Network Activity and Response to Hypoxic Challenge

By 49 DIV, the motor neuron networks had largely reached a stable state in terms of activity, burst propensity and synchrony. Following the hypoxic challenge at 54 DIV, activity remained fairly stable for all groups until 59–60 DIV. The normoxic networks remained stable after this point and throughout the recording period. The hypoxic networks, however, exhibited fluctuations in activity after this point. Both inhibited and uninhibited hypoxic networks showed consistent changes in firing rate from day to day (Figure 3A), alternating between high (inhibited 47.27 ± 0.90 Hz, uninhibited 21.42 ± 2.80 Hz) and low (inhibited 11.61 ± 0.76 Hz, uninhibited 2.00 ± 0.54 Hz). Burst propensity and synchrony were different between the inhibited and the uninhibited hypoxic networks. The inhibited networks maintained entirely stable burst propensity comparable to control networks, while the uninhibited networks fluctuated in this regard in a similar way as seen in firing rate (Figure 3B). The uninhibited and inhibited hypoxic networks both showed similar fluctuations from low to high synchrony in co-occurrence with the fluctuations in firing rate and burst propensity, although in an opposite manner (Figure 3C).

FIGURE 3
www.frontiersin.org

Figure 3. Activity of neural networks before and after the hypoxic challenge. The activity of the motor neuron networks is described in terms of firing rate (A), fraction of spikes in network bursts, or bursting propensity (B), and coherence index, a measure of network synchrony (C). The stippled line at 54 DIV indicates when the hypoxic challenge was carried out, and the data points presented for this day are baseline recordings pre-hypoxia. Red labels along the Culture Age axis indicates that cell media was replaced these days. Lines and shaded regions indicate median ± median absolute deviation.

To assess the fluctuations of the network activity, the relative dispersion in the longitudinal data for each network was evaluated by RCVQ prior to and after the hypoxic challenge. The results of the statistical tests are shown in Table 2. Pre-hypoxia, there was no significant difference between the longitudinal dispersion of any electrophysiology parameter across the experimental groups, as seen in Figures 4A–C. As shown in Figure 4D, the hypoxic networks subsequently exhibited significant increases in dispersion of firing rate compared to normoxic networks, regardless of inhibition status (p < 0.0083). A similar relationship was observed in dispersion of synchrony following hypoxia, as seen in Figure 4F, which was also seen regardless of inhibition status (p < 0.0083). The dispersion in burst propensity following hypoxia was significantly increased in uninhibited hypoxic networks compared to normoxic control networks (p < 0.0083), seen in Figure 4E. However, there was no significant difference between the dispersion in burst propensity of normoxic control networks and inhibited hypoxic networks. Furthermore, the latter group showed significantly less dispersion in burst propensity compared to uninhibited hypoxic networks, suggesting that GABA inhibition contributed to maintained bursting activity in motor neuron networks.

TABLE 2
www.frontiersin.org

Table 2. Statistical comparison of network activity.

FIGURE 4
www.frontiersin.org

Figure 4. Dispersion in network activity pre and post the hypoxic challenge. The dispersion of the electrophysiology parameters was assessed by the RCVQ of the activity of each network and compared by Wilcoxon rank sum test, prior to and following the hypoxic challenge. The longitudinal dispersion pre-hypoxia showed no significant differences in terms of firing rate (A), bursting propensity (B), or synchrony (C). Post-hypoxia the hypoxic networks showed similar significant increases in longitudinal dispersion in terms of firing rate (D) and synchrony (F). However, only the uninhibited hypoxic networks showed a significant increase in the dispersion of bursting propensity, while the inhibited hypoxic networks showed significantly less longitudinal dispersion, in line with control networks (E). *p < 0.0083, **p < 0.0017, ***p < 0.00017. Bars and error bars indicate median ± median absolute deviation.

In summary, the uninhibited hypoxic networks demonstrated changes in activity which only became evident 5–6 days after the hypoxic challenge, consisting of alternating decreases in firing rate and burst propensity and increases in network synchrony. The inhibited hypoxic networks on the other hand, showed the same fluctuations in firing rate and synchrony, but not in burst propensity. Normoxic control networks, did not exhibit any fluctuations regardless of inhibition status. The fluctuations led to significant increases in longitudinal dispersion of the network activity parameters for the hypoxic networks, with the exception of the dispersion in burst propensity of the inhibited hypoxic networks, which remained similar to the levels seen in normoxic control networks and significantly lower than uninhibited hypoxic networks.

Of note is the timing of the fluctuations in terms of media changes in handling the networks. Points of high firing rate occurred closer to the previous media change (about 24 h), while points of low firing rate coincided with longer time elapsed since the previous media change, indicated by red x-axis labels in Figure 3 (about 48 h).

Hypoxic Challenge

Neuronal activity was recorded during hypoxia, and there were no clear differences observed between the hypoxic and the normoxic networks in terms of firing rate. All the inhibited networks remained quiet, with firing rate near or at zero throughout the recording period (Figure 5A). Furthermore, there were no clear differences between the activity displayed during the baseline observations and the hypoxic challenge, for either normoxic controls or hypoxic networks, in any assessed network activity parameter (Figures 5B–G).

FIGURE 5
www.frontiersin.org

Figure 5. Hypoxic challenge and GABA inhibition. The hypoxia exposed, uninhibited networks showed similar activity to uninhibited networks normoxic networks, while inhibited networks were completely silenced. (A) Data during the hypoxia session was divided into 1-min bins and is shown here throughout the 1-h recording. Uninhibited networks showed stable activity throughout the episode, while inhibited networks were completely silenced. Lines and shaded regions indicate median ± median absolute deviation. (B–G) There were no significant differences between the baseline electrophysiological activity and the activity during the hypoxic exposure in terms of firing rate, burst propensity or synchrony, for normoxic controls or hypoxic networks, assessed by Wilcoxon rank sum test, all p > 0.05. Bars and error bars indicate median ± median absolute deviation.

Discussion

By the time the motor neuron networks in this study were mature at 49 DIV they exhibited stable activity in terms of firing rate, a high degree of burst propensity and stable levels of synchrony as seen up to 54 DIV in Figure 3, which are in line with observations for functionally mature networks (Pasquale et al., 2008). These networks also expressed markers of mature neuron cytoskeletons and relevant receptors including GABAb-1 receptors, shown in Figure 2C, as well as specific markers for motor neurons, shown in Figures 2A,B.

Interestingly, the hypoxic networks showed no changes in activity during the episode itself compared to the normoxic controls. There were no discernible differences in firing rate, burst firing or synchrony in the uninhibited hypoxic networks compared to the pre-hypoxia baseline (Figures 5C,E,G). The same was predictably true for uninhibited normoxic networks (Figures 5B,D,F). Previous reports have indicated that exposure to hypoxia can lead to acute hyperexcitability leading to elevated firing rate in hypothalamic neurons (Shonis and Waldrop, 1995) while others have found that chemical hypoxia leads to reduced excitability (Englund et al., 2004), although investigations of ongoing activity during hypoxia rather than post-hypoxic recovery are rare, especially in motor neurons. Such changes in activity were not observable in the current motor neuron networks, and changes in activity following exposure to hypoxic challenge only became evident later on.

It should be noted that the astrocytes co-cultured in the motor neuron networks may provide neuroprotective effects to the motor neurons in this study, a widely reported phenomenon (Marina et al., 2016). However, since all experimental groups contained the same ratio of neurons to astrocytes, and evidence suggests that astrocytes are more resilient than neurons to hypoxic conditions (Alves et al., 2000), it seems unlikely that potential effects of astrocytes is a major contributing factor to the observed group differences.

This study demonstrated that mature networks of motor neurons exposed to hypoxic challenge exhibit aberrant fluctuations in network activity which were not evident in networks maintained in normoxic conditions, seen clearly in Figure 3. The fluctuations occurred over the course of 59–64 DIV, and involved changes from high to low firing rate, from normoxic levels of burst propensity to low levels, and from normoxic levels of network synchrony to elevated levels. This loss of homoeostasis and stable network function differs significantly from the normoxic networks, as seen in Figure 4, showing the detrimental effects of the hypoxic exposure on network function. Interestingly, this effect did not become apparent until 5–6 days following the hypoxic challenge, as seen in Figure 3. This may suggest that there are compensatory mechanisms or adaptive responses which can maintain function close to that of normoxic counterparts, but that these measures are only effective for a limited time.

Furthermore, the timing of the observed changes in activity follow a very characteristic pattern, with deviations from normoxic-like levels observed on alternate days which correspond to observations immediately prior to cell media replacement, meaning that more of the nutrients in the cell media were depleted compared to the observations where network activity was more similar to that of normoxic networks. It is possible that this is a result of an increased metabolic burden due to effects of the hypoxia and any homeostatic mechanisms which fail to maintain as nutrients are depleted. High metabolic requirements are seen across multiple populations of neurons vulnerable to neurodegenerative disease (Muddapu et al., 2020) and may contribute to their selective vulnerability. Similarly, during hypoxia, due to the uncoupling of oxidative phosphorylation, the production of ATP becomes less efficient, eventually leading to elevated oxidative stress (Jain et al., 2015; Snyder et al., 2017; Chen et al., 2018). The effects of hypoxia can thus be more adverse for neuronal populations with high energy requirements. The results from this study may therefore suggest that motor neuron networks exposed to hypoxia exhibit increased metabolic stress, which leads to loss of functional homoeostasis and could be a predisposing factor for further network degeneration.

Some of the changes in network activity appear to be related to activity during the hypoxic challenge, since inhibition of the network with GABA improves functional outcomes in terms of maintaining network burst propensity. While inhibition of the normoxic networks appeared to have no lasting effect, the hypoxic networks inhibited during the hypoxic challenge exhibited normoxic levels of burst propensity throughout the experiment (Figures 3B, 4E), despite exhibiting similar fluctuations as the uninhibited hypoxic networks in terms of the other activity parameters (Figures 3A,C, 4D,F). Network bursts are an emergent property of activity of networks of neurons and is considered to be more reliable in terms of signalling and can carry more information (Zeldenrust et al., 2018). The maintenance of this parameter in the silenced hypoxic networks indicates improved functional outcomes as a result of the silencing, even though there are still clear consequences of the hypoxic challenge. It is possible that there are activity-dependent mechanisms which occur during hypoxia which lead to a delayed loss of function, either due to plasticity-induced changes during the episode or due to metabolic stress incurred by maintaining activity during hypoxic challenge. Acute hypoxia like in this study has been reported to activate HIF-1 through mitochondria-dependent mechanisms, a function of which is to upregulate anaerobic respiration and inhibiting the oxygen dependent TCA cycle, which places the affected cells in an energy-deficient state (Solaini et al., 2010; Hernandez-Gerez et al., 2019). From the results of our study, it appears that the activity-dependent damage was due to maintenance of normal activity, indicated by the results in Figure 4, rather than due to changes in activity during the hypoxic challenge. Maintaining activity during such circumstances of less energy production can lead to increased production of ROS which have the potential to cause damage to the network. Elevated levels of oxidative stress markers have previously been reported in mouse models of ALS exposed to hypoxic conditions (Kim et al., 2013). It seems likely that inhibition the networks during the hypoxic challenge would reduce the energy requirements of the motor neurons, thus reducing the metabolic stress. Inhibiting the networks contributed to maintaining pre-hypoxic functionality in terms of network burst propensity, although the persistent fluctuations in firing rate and synchrony indicate that the hypoxic challenge had effects on the networks which were not dependent on activity during the hypoxic challenge itself.

Although there are conflicting reports on the immediate effects of hypoxia exposure on network activity, it is evident that the mechanisms involved, at least in part, increased oxidative damage through calcium dysregulation (Kristian and Siesjo, 1998; Starkov et al., 2004). In this context, there is a particular association with the activity of NMDA-receptors and findings suggest that blocking them can reduce the production of ROS during hypoxia (Lafon-Cazal et al., 1993) as well as ameliorate post-hypoxic hyperexcitability (Godukhin et al., 2002). When addressing hypoxia and considering interventions it is vital to know what pathways can be targetted in order to improve outcomes for patients, for example those who require assisted ventilation, which can be associated with acute hypoxic events (Keith and Pierson, 1996; Mahmood et al., 2013). This study suggests that functional outcomes can be preserved by targetting activity-dependent pathways.

The study of hypoxia and the mechanisms by which it damages neural networks bears a direct relevance to neurodegenerative diseases like ALS. The body of evidence, which shows that mutations in hypoxia-responsive elements like angiogenin (ANG) can increase the likelihood of developing ALS (Greenway et al., 2006; van Es et al., 2011; Zou et al., 2012), indicates that some mechanisms of healthy responses to hypoxia are also involved in the ongoing health of motor neurons which fail in ALS. Additionally, SOD1 genetic models of ALS exhibit altered responses to hypoxia, which may exacerbate their susceptibility to hypoxic stress (Cimini et al., 2014), and elevated levels of hypoxia responsive elements including ANG, apoptosis inducing factor and HIF-1α (Zhang et al., 2011) are observed in both ALS patients (Cronin et al., 2006; van Es et al., 2014; Morelli et al., 2020) and mouse models (Xu et al., 2011). Furthermore, exposure to repeated hypoxic episodes appears to increase the risk of developing ALS (Vanacore et al., 2010), and studies in mouse models of ALS have found that stabilising HIF-1α, in effect enhancing the hypoxic response mechanism, reduced degeneration of spinal motor neurons and myofibrils of the lower limbs, as well as increasing the lifespan of the mice (Nomura et al., 2019). Mouse models of ALS exposed to hypoxic conditions also exhibit accelerated disease progression compared to normoxic controls (Kim et al., 2013). These findings suggest a relationship between the mechanisms of resilience to hypoxia and the risk of developing ALS as well as the severity of the disease progression. The results of this study indicate that there are activity-dependent mechanisms which lead to network dysfunction following hypoxia. Not unexpectedly, studies of ALS show that activity-dependent processes, hyperexcitability and formation of ROS is directly linked with motor neuron damage (Kaus and Sareen, 2015; Verma et al., 2018; Marcuzzo et al., 2019). ALS has also been observed to alter cortical hub regions, i.e., highly interconnected structures which are vital for proper network function, early on in the disease process (Verstraete et al., 2011; Ma et al., 2015). Given that cortical hubs are known to be highly metabolically demanding (Bullmore and Sporns, 2012; Collin et al., 2014), it is natural to consider that they might be some of the first parts of the neural network to show deterioration in circumstances where there is reduced capacity to respond healthily to activity-dependent adversity, such as a reduction in functionality of the hypoxia-responsive machinery.

In conclusion, this study examined the responses of motor neuron networks to a hypoxic challenge and how the functionality of the networks develops over time, with and without inhibition of activity during the challenge. We found that although activity during the hypoxic challenge remained similar to both baseline measurements and to normoxic controls, hypoxic networks began to show clear deviations from the behaviour of normoxic networks, involving fluctuations in activity. Inhibition of network activity during the hypoxic challenge contributed to preserved bursting behaviour compared to uninhibited networks, an improved functional outcome. These responses could not have been captured without a longitudinal assessment of the network activity, which is rarely seen in the existing literature. Elucidating patterns of dysfunction following hypoxic challenge and the mechanisms involved has implications beyond the hypoxia itself, such as in diseases like ALS, where the mechanisms of degeneration appear to be at least in part shared with those observed in hypoxia. Assessing these mechanisms can thus contribute to a more precise investigations of these properties in a highly multifactorial neurodegenerative disorder.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

VF conceived of and designed the study, acquired the data, performed the statistical analysis, and drafted the manuscript. AS and IS participated in the design of the study and helped to revise the manuscript. All authors read and approved the final version of the manuscript.

Funding

The authors would like to acknowledge financial support from ALS Norge, the Liaison Committee for Education, Research and Innovation in Central Norway (Samarbeidsorganet HMN- NTNU) Grant Number 2018/42521, and the Joint Research Committee (Felles Forskningsutvalg–St. Olav’s Hospital HF and the Faculty of Medicine, NTNU).

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 Morten Høydal, Department of Circulation and Medical Imaging, NTNU, for facilitating access to the Axion electrophysiology system.

References

Alves, P. M., Fonseca, L. L., Peixoto, C. C., Almeida, A. C., Carrondo, M. J., and Santos, H. (2000). NMR studies on energy metabolism of immobilized primary neurons and astrocytes during hypoxia, ischemia and hypoglycemia. NMR Biomed. 13, 438–448. doi: 10.1002/nbm.665

PubMed Abstract | CrossRef Full Text | Google Scholar

Arachchige, C. N. P. G., Prendergast, L. A., and Staudte, R. G. (2020). Robust analogs to the coefficient of variation. J. Appl. Stat. 1–23. doi: 10.1080/02664763.2020.1808599

CrossRef Full Text | Google Scholar

Bhatia, D., Ardekani, M. S., Shi, Q., and Movafagh, S. (2017). “Hypoxia and its emerging therapeutics in neurodegenerative, inflammatory and renal diseases,” in Hypoxia and Human Diseases, eds J. Zheng and C. Zhou (London: IntechOpen). doi: 10.1016/j.pharmthera.2016.03.018

CrossRef Full Text | Google Scholar

Bullmore, E., and Sporns, O. (2012). The economy of brain network organization. Nat. Rev. Neurosci. 13, 336–349. doi: 10.1038/nrn3214

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, R., Lai, U. H., Zhu, L., Singh, A., Ahmed, M., and Forsyth, N. R. (2018). Reactive oxygen species formation in the brain at different oxygen levels: the role of hypoxia inducible factors. Front. Cell Dev. Biol. 6:132. doi: 10.3389/fcell.2018.00132

CrossRef Full Text | Google Scholar

Cimini, S., Rizzardini, M., Biella, G., and Cantoni, L. (2014). Hypoxia causes autophagic stress and derangement of metabolic adaptation in a cell model of amyotrophic lateral sclerosis. J. Neurochem. 129, 413–425. doi: 10.1111/jnc.12642

PubMed Abstract | CrossRef Full Text | Google Scholar

Collin, G., Sporns, O., Mandl, R. C., and van den Heuvel, M. P. (2014). Structural and functional aspects relating to cost and benefit of rich club organization in the human cerebral cortex. Cereb. Cortex 24, 2258–2267. doi: 10.1093/cercor/bht064

PubMed Abstract | CrossRef Full Text | Google Scholar

Cronin, S., Greenway, M. J., Ennis, S., Kieran, D., Green, A., Prehn, J. H., et al. (2006). Elevated serum angiogenin levels in ALS. Neurology 67, 1833–1836. doi: 10.1212/01.wnl.0000244466.46020.47

CrossRef Full Text | Google Scholar

Dupuis, L., Pradat, P. F., Ludolph, A. C., and Loeffler, J. P. (2011). Energy metabolism in amyotrophic lateral sclerosis. Lancet Neurol. 10, 75–82.

Google Scholar

Englund, M., Bjurling, M., Edin, F., Hyllienmark, L., and Brismar, T. (2004). Hypoxic excitability changes and sodium currents in hippocampus CA1 neurons. Cell Mol. Neurobiol. 24, 685–694. doi: 10.1023/b:cemn.0000036405.53992.78

PubMed Abstract | CrossRef Full Text | Google Scholar

Gavello, D., Rojo-Ruiz, J., Marcantoni, A., Franchino, C., Carbone, E., and Carabelli, V. (2012). Leptin counteracts the hypoxia-induced inhibition of spontaneously firing hippocampal neurons: a microelectrode array study. PLoS One 7:e41530. doi: 10.1371/journal.pone.0041530

CrossRef Full Text | Google Scholar

Godukhin, O., Savin, A., Kalemenev, S., and Levin, S. (2002). Neuronal hyperexcitability induced by repeated brief episodes of hypoxia in rat hippocampal slices: involvement of ionotropic glutamate receptors and L-type Ca(2+) channels. Neuropharmacology 42, 459–466. doi: 10.1016/s0028-3908(02)00005-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Greenway, M. J., Andersen, P. M., Russ, C., Ennis, S., Cashman, S., Donaghy, C., et al. (2006). ANG mutations segregate with familial and ‘sporadic’ amyotrophic lateral sclerosis. Nat. Genet. 38, 411–413. doi: 10.1038/ng1742

PubMed Abstract | CrossRef Full Text | Google Scholar

Gu, X. Q., and Haddad, G. G. (2001). Decreased neuronal excitability in hippocampal neurons of mice exposed to cyclic hypoxia. J. Appl. Physiol. 91, 1245–1250. doi: 10.1152/jappl.2001.91.3.1245

PubMed Abstract | CrossRef Full Text | Google Scholar

Guglielmotto, M., Aragno, M., Autelli, R., Giliberto, L., Novo, E., Colombatto, S., et al. (2009). The up-regulation of BACE1 mediated by hypoxia and ischemic injury: role of oxidative stress and HIF1alpha. J. Neurochem. 108, 1045–1056. doi: 10.1111/j.1471-4159.2008.05858.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hernandez-Gerez, E. I, Fleming, N., and Parson, S. H. (2019). A role for spinal cord hypoxia in neurodegeneration. Cell Death Dis. 10:861. doi: 10.1038/s41419-019-2104-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Hor, J. H., Santosa, M. M., Lim, V. J. W., Ho, B. X., Taylor, A., Khong, Z. J., et al. (2021). ALS motor neurons exhibit hallmark metabolic defects that are rescued by SIRT3 activation. Cell Death Differ. 28, 1379–1397. doi: 10.1038/s41418-020-00664-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Hossmann, K.-A. (1999). “The hypoxic brain,” in Hypoxia: Into the Next Millennium, eds R. C. Roach, P. D. Wagner, and P. H. Hackett (Boston, MA: Springer US), 155–169.

Google Scholar

Jain, K., Prasad, D., Singh, S. B., and Kohli, E. (2015). Hypobaric hypoxia imbalances mitochondrial dynamics in rat brain hippocampus. Neurol. Res. Int. 2015:742059. doi: 10.1155/2015/742059

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaus, A., and Sareen, D. (2015). ALS patient stem cells for unveiling disease signatures of motoneuron susceptibility: perspectives on the deadly mitochondria, er stress and calcium triad. Front. Cell Neurosci. 9:448. doi: 10.3389/fncel.2015.00448

CrossRef Full Text | Google Scholar

Keith, R. L., and Pierson, D. J. (1996). Complications of mechanical ventilation. A bedside approach. Clin. Chest. Med. 17, 439–451. doi: 10.1016/s0272-5231(05)70326-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, S. M., Kim, H., Lee, J. S., Park, K. S., Jeon, G. S., Shon, J., et al. (2013). Intermittent hypoxia can aggravate motor neuronal loss and cognitive dysfunction in ALS mice. PLoS One 8:e81808. doi: 10.1371/journal.pone.0081808

CrossRef Full Text | Google Scholar

Kristian, T., and Siesjo, B. K. (1998). Calcium in ischemic cell death. Stroke 29, 705–718. doi: 10.1161/01.str.29.3.705

CrossRef Full Text | Google Scholar

Kumar, A., and Goyal, R. (2008). Gabapentin attenuates acute hypoxic stress-induced behavioral alterations and oxidative damage in mice: possible involvement of GABAergic mechanism. Indian J. Exp. Biol. 46, 159–163.

PubMed Abstract | Google Scholar

Lafon-Cazal, M., Pietri, S., Culcasi, M., and Bockaert, J. (1993). NMDA-dependent superoxide production and neurotoxicity. Nature 364, 535–537. doi: 10.1038/364535a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, X., Zhang, J., Zhang, Y., Chen, H., Li, R., Wang, J., et al. (2015). Altered cortical hubs in functional brain networks in amyotrophic lateral sclerosis. Neurol Sci. 36, 2097–2104. doi: 10.1007/s10072-015-2319-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Mahmood, N. A., Chaudry, F. A., Azam, H., Ali, M. I., and Khan, M. A. (2013). Frequency of hypoxic events in patients on a mechanical ventilator. Int. J. Crit. Illn. Inj. Sci. 3, 124–129. doi: 10.4103/2229-5151.114272

PubMed Abstract | CrossRef Full Text | Google Scholar

Marcuzzo, S., Terragni, B., Bonanno, S., Isaia, D., Cavalcante, P., Cappelletti, C., et al. (2019). Hyperexcitability in cultured cortical neuron networks from the G93A-SOD1 amyotrophic lateral sclerosis model mouse and its molecular correlates. Neuroscience 416, 88–99. doi: 10.1016/j.neuroscience.2019.07.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Marina, N., Kasymov, V., Ackland, G. L., Kasparov, S., and Gourine, A. V. (2016). Astrocytes and brain hypoxia. Adv. Exp. Med. Biol. 903, 201–207. doi: 10.1007/978-1-4899-7678-9_14

CrossRef Full Text | Google Scholar

Morelli, C., Tiloca, C., Colombrita, C., Zambon, A., Soranna, D., Lafronza, A., et al. (2020). CSF angiogenin levels in amyotrophic lateral sclerosis-frontotemporal dementia spectrum. Amyotroph. Lateral Scler. Frontotemporal. Degener. 21, 63–69. doi: 10.1080/21678421.2019.1704016

PubMed Abstract | CrossRef Full Text | Google Scholar

Muddapu, V. R., Dharshini, S. A. P., Chakravarthy, V. S., and Gromiha, M. M. (2020). Neurodegenerative diseases - is metabolic deficiency the root cause? Front. Neurosci. 14:213. doi: 10.3389/fnins.2020.00213

CrossRef Full Text | Google Scholar

Nomura, E., Ohta, Y., Tadokoro, K., Shang, J., Feng, T., Liu, X., et al. (2019). Imaging hypoxic stress and the treatment of amyotrophic lateral sclerosis with dimethyloxalylglycine in a mice model. Neuroscience 415, 31–43. doi: 10.1016/j.neuroscience.2019.06.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Pasquale, V., Massobrio, P., Bologna, L. L., Chiappalone, M., and Martinoia, S. (2008). Self-organization and neuronal avalanches in networks of dissociated cortical neurons. Neuroscience 153, 1354–1369. doi: 10.1016/j.neuroscience.2008.03.050

PubMed Abstract | CrossRef Full Text | Google Scholar

Schiff, S. J., and Somjen, G. G. (1985). Hyperexcitability following moderate hypoxia in hippocampal tissue slices. Brain Res. 337, 337–340. doi: 10.1016/0006-8993(85)90071-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Shonis, C. A., and Waldrop, T. G. (1995). In vitro effects of GABA and hypoxia on posterior hypothalamic neurons from spontaneously hypertensive and Wistar-Kyoto rats. Brain Res. Bull. 36, 461–466. doi: 10.1016/0361-9230(94)00213-k

PubMed Abstract | CrossRef Full Text | Google Scholar

Snyder, B., Shell, B., Cunningham, J. T., and Cunningham, R. L. (2017). Chronic intermittent hypoxia induces oxidative stress and inflammation in brain regions associated with early-stage neurodegeneration. Physiol. Rep. 5:e13258. doi: 10.14814/phy2.13258

PubMed Abstract | CrossRef Full Text | Google Scholar

Solaini, G., Baracca, A., Lenaz, G., and Sgarbi, G. (2010). Hypoxia and mitochondrial oxidative metabolism. Biochim. Biophys. Acta. 1797, 1171–1177. doi: 10.1016/j.bbabio.2010.02.011

CrossRef Full Text | Google Scholar

Starkov, A. A., Chinopoulos, C., and Fiskum, G. (2004). Mitochondrial calcium and oxidative stress as mediators of ischemic brain injury. Cell Calcium 36, 257–264. doi: 10.1016/j.ceca.2004.02.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, X., He, G., Qing, H., Zhou, W., Dobie, F., Cai, F., et al. (2006). Hypoxia facilitates Alzheimer’s disease pathogenesis by up-regulating BACE1 gene expression. Proc. Natl. Acad. Sci. U.S.A. 103, 18727–18732. doi: 10.1073/pnas.0606298103

PubMed Abstract | CrossRef Full Text | Google Scholar

Turovsky, E. A., Turovskaya, M. V., Kononov, A. V., and Zinchenko, V. P. (2013). Short-term episodes of hypoxia induce posthypoxic hyperexcitability and selective death of GABAergic hippocampal neurons. Exp. Neurol. 250, 1–7. doi: 10.1016/j.expneurol.2013.09.006

PubMed Abstract | CrossRef Full Text | Google Scholar

van Es, M. A., Schelhaas, H. J., van Vught, P. W., Ticozzi, N., Andersen, P. M., Groen, E. J., et al. (2011). Angiogenin variants in Parkinson disease and amyotrophic lateral sclerosis. Ann. Neurol. 70, 964–973.

Google Scholar

van Es, M. A., Veldink, J. H., Schelhaas, H. J., Bloem, B. R., Sodaar, P., van Nuenen, B. F., et al. (2014). Serum angiogenin levels are elevated in ALS, but not Parkinson’s disease. J. Neurol. Neurosurg. Psychiatry 85, 1439–1440. doi: 10.1136/jnnp-2013-307168

PubMed Abstract | CrossRef Full Text | Google Scholar

Vanacore, N., Cocco, P., Fadda, D., and Dosemeci, M. (2010). Job strain, hypoxia and risk of amyotrophic lateral sclerosis: results from a death certificate study. Amyotroph. Lateral Scler. 11, 430–434. doi: 10.3109/17482961003605796

PubMed Abstract | CrossRef Full Text | Google Scholar

Verma, M., Wills, Z., and Chu, C. T. (2018). Excitatory dendritic mitochondrial calcium toxicity: implications for parkinson’s and other neurodegenerative diseases. Front. Neurosci. 12:523. doi: 10.3389/fnins.2018.00523

CrossRef Full Text | Google Scholar

Verstraete, E., Veldink, J. H., Mandl, R. C., van den Berg, L. H., and van den Heuvel, M. P. (2011). Impaired structural motor connectome in amyotrophic lateral sclerosis. PLoS One 6:e24239. doi: 10.1371/journal.pone.0024239

CrossRef Full Text | Google Scholar

Watzlawik, J. O., Kahoud, R. J., O’Toole, R. J., White, K. A., Ogden, A. R., Painter, M. M., et al. (2015). Abbreviated exposure to hypoxia is sufficient to induce CNS dysmyelination, modulate spinal motor neuron composition, and impair motor development in neonatal mice. PLoS One 10:e0128007. doi: 10.1371/journal.pone.0128007

CrossRef Full Text | Google Scholar

Xu, R., Wu, C., Zhang, X., Zhang, Q., Yang, Y., Yi, J., et al. (2011). Linking hypoxic and oxidative insults to cell death mechanisms in models of ALS. Brain Res. 1372, 133–144. doi: 10.1016/j.brainres.2010.11.056

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeldenrust, F., Wadman, W. J., and Englitz, B. (2018). Neural coding with bursts-current state and future perspectives. Front. Comput. Neurosci. 12:48. doi: 10.3389/fncom.2018.00048

CrossRef Full Text | Google Scholar

Zhang, F., Niu, L., Li, S., and Le, W. (2019). Pathological impacts of chronic hypoxia on Alzheimer’s Disease. ACS Chem. Neurosci. 10, 902–909. doi: 10.1021/acschemneuro.8b00442

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Yan, J., Chang, Y., ShiDu Yan, S., and Shi, H. (2011). Hypoxia inducible factor-1 as a target for neurodegenerative diseases. Curr. Med. Chem. 18, 4335–4343. doi: 10.2174/092986711797200426

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, T., Wang, F. Z., Huang, Y. H., Yao, H., Liu, Z. W., and Fan, M. (2003). [Protective effect of GABA against hypoxic injury in rat hippocampal slices]. Zhongguo Ying Yong Sheng Li Xue Za Zhi 19, 16–19.

PubMed Abstract | Google Scholar

Zou, Z. Y., Wang, X. N., Liu, M. S., Sun, Q., Li, X. G., Cui, L. Y., et al. (2012). Identification of a novel missense mutation in angiogenin in a Chinese amyotrophic lateral sclerosis cohort. Amyotroph. Lateral Scler. 13, 270–275. doi: 10.3109/17482968.2011.643900

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: hypoxia, multielectrode array (MEA) recording, motor neuron disease, activity-dependent mechanisms, longitudinal, network activity

Citation: Fiskum V, Sandvig A and Sandvig I (2021) Silencing of Activity During Hypoxia Improves Functional Outcomes in Motor Neuron Networks in vitro. Front. Integr. Neurosci. 15:792863. doi: 10.3389/fnint.2021.792863

Received: 11 October 2021; Accepted: 01 December 2021;
Published: 16 December 2021.

Edited by:

Pierre Pouget, Université Pierre et Marie Curie, France

Reviewed by:

Chiara F. Valori, German Center for Neurodegeneratives, Helmholtz Association of German Research Centers (HZ), Germany
Johannes Boltze, University of Warwick, United Kingdom

Copyright © 2021 Fiskum, Sandvig and Sandvig. 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: Vegard Fiskum, dmVnYXJkLmZpc2t1bUBudG51Lm5v

Disclaimer: 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.