Neural Correlates of Vestibular Processing During a Spaceflight Analog With Elevated Carbon Dioxide (CO2): A Pilot Study

Astronauts return to Earth from spaceflight missions with impaired mobility and balance; recovery can last weeks postflight. This is due in large part to the altered vestibular signaling and sensory reweighting that occurs in microgravity. The neural mechanisms of spaceflight-induced vestibular changes are not well understood. Head-down-tilt bed rest (HDBR) is a common spaceflight analog environment that allows for study of body unloading, fluid shifts, and other consequences of spaceflight. Subjects in this context still show vestibular changes despite being in Earth’s gravitational environment, potentially due to sensory reweighting. Previously, we found evidence of sensory reweighting and reduced neural efficiency for vestibular processing in subjects who underwent a 70-day HDBR intervention. Here we extend this work by evaluating the impact of HDBR paired with elevated carbon dioxide (CO₂) to mimic International Space Station conditions on vestibular neural processing. Eleven participants (6 males, 34 ± 8 years) completed 30 days of HDBR combined with 0.5% atmospheric CO₂ (HDBR + CO₂). Participants underwent six functional magnetic resonance imaging (fMRI) sessions pre-, during, and post- HDBR + CO₂ while we measured brain activity in response to pneumatic skull taps (a validated method of vestibular stimulation). We also measured mobility and balance performance several times before and after the intervention. We found support for adaptive neural changes within the vestibular system during bed rest that subsequently recovered in several cortical and cerebellar regions. Further, there were multiple brain regions where greater pre- to post- deactivation was associated with reduced pre- to post- balance declines. That is, increased deactivation of certain brain regions associated with better balance post-HDBR + CO₂. We also found that, compared to HDBR alone (n = 13 males; 29 ± 3 years) HDBR + CO₂ is associated with greater increases in activation of multiple frontal, parietal, and temporal regions during vestibular stimulation. This suggests interactive or additive effects of bed rest and elevated CO₂. Finally, we found stronger correlations between pre- to post- HDBR + CO₂ brain changes and dependence on the visual system during balance for subjects who developed signs of Spaceflight-Associated Neuro-ocular Syndrome (SANS). Together, these findings have clear implications for understanding the neural mechanisms of bed rest and spaceflight-related changes in vestibular processing, as well as adaptation to altered sensory inputs.

Introduction

Microgravity exposure poses unique challenges to human physiology: astronauts encounter body unloading, headward fluid shifts, altered vestibular and proprioceptive inputs, inflight and postflight spatial disorientation (Young et al., 1984), and confined quarters with carbon dioxide (CO₂) levels up to more than ten times higher than those on Earth (Law et al., 2014). Upon return to Earth, astronauts present with multi-systemic consequences, such as declining bone (Sibonga, 2013) and muscle mass (LeBlanc et al., 1995; Stein, 2013), cardiovascular changes (Hargens and Richardson, 2009), and mobility and balance difficulties (Mulavara et al., 2010; Cohen et al., 2012; Wood et al., 2015). Here we focus on the neural vestibular consequences of a spaceflight analog environment, as well as the neural mechanisms underlying declines in vestibularly mediated mobility and balance.

Animal studies have demonstrated peripheral vestibular changes with spaceflight; for instance, utricular afferents become hypersensitive to translational accelerations after return to Earth (Boyle et al., 2001). Although the specific mechanisms for these changes remain unknown, one possibility is that the brain reinterprets afferent sensory input during flight due to the lack of a gravitational reference vector for the otoliths. After return to Earth, this re-interpretation is in conflict with Earth’s gravitational environment and results in postflight vestibular dysfunction (e.g., balance difficulties), followed by slow re-adaptation over the days and weeks following spaceflight (Young et al., 1984; Parker et al., 1985; Mulavara et al., 2010). Astronauts also present with decreased skin sensitivity on the soles of the feet following spaceflight (Lowrey et al., 2014), which has been attributed to in-flight sensory reweighting (i.e., the process of adjusting the magnitude of different sensory contributions to motor control) (Assländer and Peterka, 2014) in compensation for unreliable vestibular inputs in microgravity. A single-subject case study (Demertzi et al., 2016) and recent study of 11 cosmonauts (Pechenkova et al., 2019) examining resting-state and task-based functional magnetic resonance imaging (fMRI) connectivity found evidence for vestibular cortex reorganization and multisensory reweighting following long-duration spaceflight. This work provides preliminary evidence of flight-related central vestibular plasticity. Taken together, it is likely that spaceflight factors influence the neural correlates of vestibular processing; however, the precise mechanisms underlying such changes require further study.

Head-down-tilt bed rest (HDBR) is a common spaceflight analog intervention that permits ground-based study of how axial body unloading alters sensory inputs that subsequently impact neural vestibular processing and vestibular system plasticity. Subjects remain in bed with their head tilted down six degrees to mimic a subset of spaceflight consequences including headward fluid shifts, arterial pressure changes, axial body unloading, and reduced somatosensory input. Although gravitational vector input does not change during HDBR, there is evidence that, similar to spaceflight, axial body unloading contributes to sensory reweighting (Moore et al., 2010; Mulavara et al., 2018; Yuan et al., 2018b). Even though HDBR does not directly affect vestibular inputs, sensory reweighting is thought to affect neural vestibular processing during HDBR; more specifically, the vestibular nuclei receive inputs from the vestibular organs, in addition to proprioceptive signals from the limbs (Fredrickson et al., 1966; Rubin et al., 1979; Yates et al., 2000; Jian et al., 2002). If either vestibular or somatosensory inputs appear to be incorrect or abnormal, the central nervous system may use information from the other system to compensate and maintain performance (Bles et al., 1984; Dieringer, 1995; Horak and Hlavacka, 2001; Carriot et al., 2015). Thus, during HDBR, in the absence of normal somatosensory inputs to the foot, vestibular processing appears to be altered, with vestibular cues weighted more heavily (Mulavara et al., 2018). HDBR also results in reduced functional mobility and decreased postural stability, which are both behaviors that depend upon the vestibular system and multisensory integration (Reschke et al., 2009; Mulder et al., 2014; Koppelmans et al., 2015, 2017; Miller et al., 2018; Mulavara et al., 2018). Thus, taken together, HDBR provides an effective environment for studying neural vestibular adaptation and has applications for both space travel and for better understanding plasticity of the vestibular system and multisensory integration.

In recent years, several fMRI-compatible vestibular stimulation methods, including auditory tone bursts and pneumatic skull taps, have been used to map central vestibular processing networks (Schlindwein et al., 2008; Noohi et al., 2017). Two meta-analyses have revealed a diffuse vestibular processing network, including portions of insular cortex, premotor cortex, inferior parietal cortex, cingulate cortex, and the superior temporal gyri (Lopez et al., 2012a; Zu Eulenburg et al., 2012). However, the most commonly activated regions across several different vestibular stimulation methods were the parietal opercular area (“OP2”) and retroinsular cortex; consequently, these regions are sometimes referred to as “vestibular cortex,” and considered to be the core regions responsible for vestibular processing. In the present work, we stimulated the vestibular system during fMRI using pneumatic skull taps, which we have previously validated in young (Noohi et al., 2017) and older adults (Noohi et al., 2019), and successfully employed in our past HDBR work (Yuan et al., 2018b). Pneumatic skull taps elicit both activation in the vestibular cortex and deactivation in cross-modal sensory regions. Across both young and older adults, we have found associations between greater deactivation of certain subcortical and cerebellar regions in response to vestibular stimulation and better static balance performance (Noohi et al., 2019). This suggests the importance of both brain activation and deactivation, potentially reflecting sensory reweighting, for successful vestibular functioning.

In our past HDBR work, using this skull tap technique, we identified longitudinal brain changes suggestive of upregulation of vestibular processing in response to reduced somatosensory input during 70 days of HDBR (Yuan et al., 2018b). We also found associations between increased frontal, parietal, and occipital brain activity and greater HDBR-related mobility and balance declines, suggestive of reduced neural efficiency post-HDBR (Yuan et al., 2018b). Further, we identified post-HDBR increases in resting state connectivity for a network including the vestibular cortex and the cerebellum (Cassady et al., 2016). These findings indicate that the neural correlates of vestibular processing are altered with HDBR and have functional implications for vestibularly mediated behaviors.

No previous work has investigated the neural correlates of vestibular processing during an intervention combining HDBR with elevated CO₂, which would better mimic the actual conditions on the International Space Station (ISS). Among other effects, exposure to heightened CO₂ increases blood flow to the brain (at least initially) due to cerebral vasodilation (Atkinson et al., 1990; Zhou et al., 2008) and mildly impairs visuomotor function (Manzey and Lorenz, 1998). Although reduced blood levels of CO₂ during voluntary hyperventilation have been associated with increased postural sway (Sakellari et al., 1997), it is unknown how elevated atmospheric CO₂ interacts with central or peripheral vestibular processing.

Here we used fMRI to test changes in the neural response to vestibular stimulation with 30 days of HDBR paired with elevated CO₂ (which we refer to as “HDBR + CO₂”). We hypothesized that similar brain changes would emerge compared to our past work (referred to as “HDBR”), including evidence for HDBR-related upregulation of vestibular networks and reduced neural efficiency. Further, we anticipated that the interaction of HDBR and elevated CO₂ would result in additive neural effects.

We addressed three primary aims in this small pilot sample (n = 11): (1) assess the time course of changes in the neural correlates of vestibular processing and recovery patterns with HDBR + CO₂; (2) examine the functional consequences of HDBR + CO₂ by associating brain changes with mobility and balance declines; and (3) characterize how HDBR + CO₂ differentially affects the neural correlates of vestibular processing compared to HDBR alone. We developed an additional, exploratory aim (4) after about half of the HDBR + CO₂ subjects developed signs of Spaceflight-Associated Neuro-ocular Syndrome (SANS) (Laurie et al., 2019), a condition which manifests with symptoms such as optic disk edema and is estimated to affect between approximately 16 and 50 percent of long-duration astronauts (i.e., those who have completed an ISS mission, which typically last for about 6 months) (Stenger et al., 2017). We characterized subgroup differences between those HDBR + CO₂ subjects who did and did not develop signs of SANS.

Materials and Methods

HDBR + CO₂

Participants

Eleven healthy individuals (six males, five females; mean ± SD age = 34 ± 8 years) provided their written informed consent and participated in 30 days of HDBR + CO₂. This intervention was implemented within the larger study, VaPER (Visual impairment intracranial pressure and Psychological:envihab Research), in which separate investigators evaluated other physiological systems. All study procedures were approved by the local ethical commission of the regional medical association, Ärztekammer Nordrhein, as well as the University of Florida and NASA Institutional Review Boards.

Testing Timeline

Subjects were admitted to:envihab at the German Aerospace Center (Deutsches Zentrum für Luft- und Raumfahrt, DLR) in Cologne, Germany 14 days before the start of HDBR + CO₂. During this time, they completed two baseline data collection (BDC) sessions (Figure 1). Subjects then underwent 30 days of six-degree HDBR with approximately 0.5% (partial pressure = 3.8 mmHg) elevated atmospheric CO₂ (HDT), to match average ISS conditions (Law et al., 2014). Subjects kept a “strict” head-down-tilt position at all times, verified by 24/7 video monitoring. Subjects were instructed to always keep at least one shoulder in contact with the mattress. They were not permitted to use a pillow or to raise or stretch their legs aside from standardized physiotherapy sessions. Subjects remained at the facility for 14 days after HDBR + CO₂ and completed two recovery (R) data collection sessions during this time.

FIGURE 1

Figure 1. Testing timeline. Top: testing timeline for the HDBR + CO₂ group, who completed 30 days of head-down-tilt bed rest (HDBR) with 0.5% elevated atmospheric CO₂. Bottom: testing timeline for the HDBR group, who completed 70 days of HDBR with normal atmospheric CO₂ levels. BDC, baseline data collection; HDT, head-down-tilt bed rest; R, recovery. Circles indicate the day for each MRI scan. Circles with asterisks represent the three time points used to create intercept and slope images for between-group comparisons. All HDBR + CO₂ subjects completed MRI scans on exactly the same days with respect to bed rest. There was some variability in testing days for the HDBR group; average day is plotted, with error bars indicating standard deviation. Mobility and balance data were collected at all time points for the HDBR + CO₂ group, with the exception of HDT 7 and HDT 29 (i.e., participants did not complete standing tasks during bed rest). One additional mobility and balance data collection took place on R0 for the HDBR + CO₂ group within ∼3 h of first standing up. Mobility and balance scores at BDC 7 and R 0 were used for brain-behavior correlations with MRI data from BDC 7 and HDT 29.

fMRI images were collected at six time points: two times pre-, two times during, and two times post-HDBR + CO₂ (Figure 1). Subjects completed mobility and balance testing on the same days as fMRI scans, except for the time points during HDBR + CO₂. Subjects completed one additional mobility and balance testing session on the first recovery day (R0). One individual began testing late and thus did not complete BDC 13; however, this individual did complete the second baseline session (BDC 7).

Head-Down-Tilt Bed Rest

Participants

Thirteen healthy individuals (all males; mean ± SD age = 29 ± 3 years) provided their written informed consent and participated in 70 days of HDBR. All study procedures were approved by the University of Michigan, University of Texas Medical Branch, and NASA Institutional Review Boards. These subjects represent a subset of the 18 total HDBR participants who received the same mode of vestibular stimulation caused by pneumatic skull taps as the HDBR + CO₂ cohort. There were no significant age differences between the HDBR + CO₂ and HDBR participants, and both cohorts passed a minimum physical fitness standard (i.e., an Air Force Class III equivalent physical examination) to participate (Lee et al., 2019a).

Testing Timeline

Participants were admitted to the NASA bed rest facility at the University of Texas Medical Branch, Galveston, TX, United States and completed two BDC sessions in the two weeks prior to starting HDBR (Figure 1). Participants underwent 70 days of HDBR with normal atmospheric CO₂ (HDT). They maintained a six-degree head-down-tilt at all times except for 30 min during each meal, when they were allowed to support their head with their hand. Subjects remained at the facility for 14 days after HDBR and completed two recovery (R) data collection sessions during this time. All participants were a part of larger bed rest studies; thus, the timelines for HDBR + CO₂ and HDBR were restricted by NASA and not identically matched between the two studies.

fMRI images were collected at seven time points: two times pre-, three times during, and two times post-HDBR (Figure 1). As we have previously reported neural vestibular changes with this intervention (Yuan et al., 2018b), here we use these fMRI data only for group comparisons with the HDBR + CO₂ group. We examine only the fMRI scans from BDC 8, HDT 8, and HDT 50, as these points fell closest in time to those collected pre- and during bed rest for the HDBR + CO₂ group and allowed us to compare slopes of change over time in the two groups. See Section “HDBR + CO₂ vs. HDBR Group Comparisons” for details on this analysis approach for making between-group comparisons that account for these differing testing timelines.

fMRI Data Collection

fMRI Acquisition

For the HDBR + CO₂ group, fMRI scans were collected on a 3 Tesla Siemens MRI scanner. A gradient echo T2^∗-weighted echo-planar imaging sequence was used to collect the fMRI scans: TR = 2.5 s, TE = 32 ms, flip angle = 90°, FOV = 192 × 192 mm, matrix = 64 × 64, slice thickness = 3.5 mm, voxel size = 3 × 3 × 3.5 mm³, 37 slices, 96 volumes. A T1-weighted gradient-echo pulse sequence was also collected with parameters: TR = 1.9 s, TE = 2.4 ms, flip angle = 9°, FOV = 250 × 250 mm, matrix = 512 × 512, slice thickness = 1.0 mm, voxel size = 0.49 × 0.49 × 1.0 mm³, 192 slices. Participants maintained the head-down-tilt position at all times using a foam wedge in the scanner. In addition, 0.5% CO₂ was continuously supplied during the HDBR + CO₂ intervention phase (through a mask and tank system when subjects were out of the environmentally controlled wing of the building).

For the HDBR group, fMRI scans were collected on a different 3 Tesla Siemens MRI scanner. A gradient echo T2^∗-weighted echo-planar imaging sequence was used to collect the fMRI scans: TR = 3.66 s, TE = 39 ms, flip angle = 90°, FOV = 240 × 240 mm, matrix = 94 × 94, slice thickness = 4 mm, slice gap = 1 mm, voxel size = 2.55 × 2.55 × 5.0 mm³, 36 slices, 66 volumes. A T1-weighted gradient-echo pulse sequence was also collected with parameters: TR = 1.9 s, TE = 2.49 ms, flip angle = 9°, FOV = 270 × 270 mm, matrix = 288 × 288, slice thickness = 0.90 mm, voxel size = 0.94 × 0.94 × 0.90 mm³, 192 slices. The HDBR participants did not maintain the head-down-tilt position in the scanner; they were supine instead.

Vestibular Stimulation

For both the HDBR + CO₂ and the HDBR cohorts, subjects received vestibular stimulation during fMRI. Subjects received skull taps via a pneumatic tactile pulse system [MR-compatible Pn Tactile Pulse System (PnTPS), Engineering Acoustics Inc.] placed over the lateral cheekbones (Noohi et al., 2017; Yuan et al., 2018b). The skull tapper used compressed air (50–55 psi) to power a small piston that delivered low-force taps (0.6 kg) to the cheekbone. We have recently shown that this approach is well tolerated by subjects, it activates vestibular cortical regions, it results in vestibular-evoked myogenic potentials in eye muscles, and it does not cause excessive head motion (Noohi et al., 2017).

Taps were delivered at 1 Hz, and each tapping block contained 24 taps. Both groups completed one fMRI run with five 24-s blocks of taps on the left cheekbone. Each block of taps was preceded and followed by 20-s rest periods. The HDBR group also completed a second run with taps to the right cheekbone; however, here we examine only the HDBR left tap run to enable direct comparison to the HDBR + CO₂ group. Of note, although the vestibular stimulation parameters and total sequence duration were identical between groups, as the HDBR + CO₂ fMRI sequence included a faster TR and more volumes (TR = 2.5 s; 96 volumes) than the HDBR sequence (TR = 3.66 s; 66 volumes), we acquired more data and thus had more statistical power for the HDBR + CO₂ group. This represents a potential limitation of the present work and is discussed further in Section “Limitations.”

For both groups, the force of the taps was sufficiently low that they did not induce head motion that was greater than for other task runs. No subject moved more than 2.1 mm within any run, which is smaller than the size of one voxel.

fMRI Preprocessing and Subject-Level Analyses

Preprocessing

Image preprocessing was completed using Statistical Parametric Mapping 12 (SPM12, version 7219) (Ashburner et al., 2016) with MatLab R2016a, version 9.0. We used a standard SPM preprocessing pipeline for fMRI (Ashburner et al., 2016). All functional images were corrected for slice timing then realigned and resliced to correct for head motion. As an additional quality check, we used the Artifact Detection Tool (ART)¹ with motion threshold = 2.5 mm and global brain signal Z threshold = 9. There were no within-session movement outliers for either group. Only one individual in the HDBR + CO₂ group had a global intensity outlier present in 4 of 96 volumes for one session; we used the subject-level covariate outputted by ART to minimize effects of these volumes on group-level analyses.

After resetting the origins of each T1 image to the anterior commissure, the T1 images were coregistered to the mean functional image with separation of [2, 1 mm]. The T1 images were segmented using the SPM12 Dartel algorithm with a sampling distance of 1 mm. The forward deformation fields from the T1 segmentation were used to normalize the functional images and the T1 to MNI space. We used 7th degree B-spline normalization for optimal performance (Ashburner et al., 2016). The warped images were spatially smoothed with an 8 mm full−width at half−maximum three−dimensional Gaussian kernel.

Subject-Level Whole Brain Statistical Analyses

At the subject level, we calculated brain activity for each participant on a voxel-by-voxel basis for left cheekbone vestibular stimulation versus rest. We set the first level masking threshold to -infinity and masked out non-brain areas using the “mask_ICV.nii” SPM intracranial volume mask. This allowed for inclusion of all voxels in the first level general linear model (GLM), as opposed to the default SPM masking threshold of 0.80, which includes in the GLM only those voxels with a mean value ≥ 80% of the global signal. We included ART-derived head motion parameters as nuisance variables in the subject-level analyses.

Cerebellar Processing

To improve normalization of the cerebellum and avoid over-stretching (Diedrichsen, 2006; Diedrichsen et al., 2009), we applied specialized processing using portions of both the CEREbellum Segmentation (CERES) (Romero et al., 2017) pipeline and the Spatially Unbiased Infratentorial and cerebellar Template (SUIT) (Diedrichsen, 2006; Diedrichsen et al., 2009) pipeline. We used CERES to segment the cerebellum from each person’s structural T1-weighted image. We then reset the origin of each individual’s cerebellum segmentation in native space to fall within the space of the segment. This allowed us to coregister each subject’s native space segmentation to the SUIT.nii template. We created binary gray matter, white matter, and full cerebellar masks from the CERES native space output and then used the suit_normalize_dartel function to obtain the Affine transformation matrix and flowfield needed to normalize these images into SUIT space.

We coregistered all of the slice timing-corrected, realigned/resliced (but not normalized) whole brain images to the T1-weighted whole brain image that was entered into the CERES pipeline and re-ran the subject-level statistical analyses described above on these non-normalized whole brain images. Then, using the Affine transformation and flowfield from normalizing the structural cerebellar segments to SUIT space, as well as each subject’s native space full cerebellar mask, we applied suit_reslice_dartel to the whole brain functional images to reslice all of the images into SUIT space. Given the small size of cerebellar structures, we applied a 2 mm smoothing kernel to the final functional cerebellar images and masked all second-level statistical results with a binary version of the SUIT.nii template, to avoid any spillover off the cerebellum due to the spatial smoothing. We performed all second-level statistical analyses described below twice: once for the whole brain (excluding the cerebellum) and a second time for only the cerebellum.

fMRI Group-Level Statistics

Neural Response to Vestibular Stimulation

To demonstrate that our pneumatic tapper method was eliciting the expected vestibular system response, we first tested the main effect of vestibular stimulation averaged across all sessions for the HDBR + CO₂ participants at peak-level p < 0.0005 (uncorrected), extent threshold = 10 voxels, controlling for age and sex differences.

Time Course of Neural Vestibular Response to HDBR + CO₂

Similar to our past work (Yuan et al., 2016, 2018a,b), we tested for regions of immediate and cumulative change during bed rest followed by both quick and gradual recovery of brain activation patterns during vestibular stimulation across all six time points. We used flexible factorial analysis (SPM’s mixed model equivalent), controlling for age and sex, assuming independence between but not within subjects, and assuming equal variances between and within subjects (Gläscher and Gitelman, 2008; Kurth et al., 2010). We used several contrast vectors as weights for the statistical analyses to test the hypothesized relative level of activation during each session. Cumulative change (Figure 2A) was modeled as a progressive increase in activity across the course of HDBR + CO₂, with a peak at the end of HDBR + CO₂, and gradual restoration after the conclusion of HDBR + CO₂. Immediate change (Figure 2B) was assumed to onset shortly after the start of HDBR + CO₂, to maintain during HDBR + CO₂, and to end shortly after the conclusion of HDBR + CO₂. We hypothesized that recovery would be either quick (Figure 2C), occurring during bed rest (i.e., between HDT 7 and HDT 29), or that recovery would be more gradual, with altered brain activation patterns still evident at HDT 29 (Figures 2A,B). We tested for both increases and decreases in activation with each of these contrast shapes. To better detect within−subject changes with the complex longitudinal models used in this pilot study, the alpha level was set at p < 0.001 (uncorrected). We report clusters that are at least 10 voxels for the whole brain and k = 5 voxels for the cerebellum.

FIGURE 2

Figure 2. Hypothesized changes in neural vestibular response to HDBR + CO₂. We hypothesized three different longitudinal patterns of brain change: (A) Cumulative change, in which brain changes would slowly increase over the course of HDBR + CO₂, followed by recovery after the conclusion of HDBR + CO₂. (B) Instant change, in which brain changes would immediately onset after the start of HDBR + CO₂, followed by recovery post-HDBR + CO₂. (C) Instant change, in which brain changes would immediately onset the start of HDBR + CO₂, but begin recovery during bed rest. We tested each of these hypotheses by using the contrast vectors shown here as weights in our longitudinal statistical model which assessed brain changes over all six time points for the HDBR + CO₂ group. Solid lines depict the positive version of each contrast; dotted lines depict the negative version of each contrast.

Correlations of Brain and Behavioral Changes With HDBR + CO₂

We computed brain activation differences during vestibular stimulation between the final pre- HDBR + CO₂ time point (BDC 7) and the final time point during HDBR + CO₂ (HDT 29). We also computed the change in mobility and balance scores from BDC 7 to the first post-HDBR + CO₂ time point, R0. To examine regions in which HDBR + CO₂ brain changes were associated with changes in mobility and balance performance, we used a one-sample t-test model controlling for age and sex and included the behavioral change score as a covariate of interest. For each model, we used the Statistical Non-Parametric Mapping (SnPM version 13)² (Nichols and Holmes, 2002) toolbox to run non-parametric permutation tests with 15,000 permutations, variance smoothing = 8 mm kernel for the whole brain analyses and 2 mm kernel for the cerebellar analyses, minimum cluster size = 10 voxels, and threshold = non-parametric p < 0.0005 (uncorrected). The SnPM toolbox is recommended for studies with small sample sizes that may not meet assumptions for parametric testing. The SnPM toolbox calculates pseudo t-statistic images and uses non-parametric permutation testing to assess for significance.

HDBR + CO₂ vs. HDBR Group Comparisons

To examine differences in neural response to vestibular stimulation between bed rest with and without elevated CO₂, we compared both baseline (i.e., intercept) differences between the HDBR + CO₂ and HDBR groups, as well as the slope of change in brain activation across bed rest. As each cohort followed a different testing timeline, we compared the three time points that fell the closest together in time between the groups (indicated by asterisks in Figure 1). As in our previous work (Yuan et al., 2016, 2018b), we calculated a regression intercept and slope for each person using the scans from these three time points. The last image collected before the start of HDBR + CO₂ or HDBR was treated as time = 0 days, assuming that pre-bed rest activation was stable. Calculating the regression intercept allowed us to examine baseline differences between groups, and calculating the regression slope allowed us to compare the rate of change in brain activation during vestibular stimulation between groups.

We used two sample t-tests to examine between-group differences in intercept and slope images. For all group comparisons, we used SnPM non-parametric permutation tests with 15,000 permutations, variance smoothing = 8 mm kernel for whole brain analyses and 2 mm kernel for cerebellar analyses, minimum cluster size = 10 voxels, and threshold = non-parametric p < 0.0005 (uncorrected). In each model, we accounted for age and sex differences. We excluded two individuals from the HDBR cohort from group analyses: one individual had severe artifacts in their HDT 50 scan, and another individual had abnormally high contrast values at the single-subject level, possibly also due to artifacts. Thus there were n = 11 subjects per group for group comparisons.

There were several differences between the HDBR + CO₂ and HDBR images. Images were collected on different scanners, HDBR images showed evidence of greater orbitofrontal dropout compared to HDBR + CO₂ images, and HDBR individuals presented with slightly smaller ventricles. To address this and remain conservative in our analyses, we do not report any between-group orbitofrontal results, and we report with caution one between-group result in close proximity to the ventricles. As these two groups represent highly unique cohorts who have undergone a rare, intensive bed rest intervention with nearly identical behavioral and neuroimaging protocols, we feel that it is still valuable to report on group differences between these two cohorts, although the results of these specific analyses should be interpreted with caution. As we previously reported on longitudinal neural vestibular changes and brain-behavior correlations for the HDBR group (Yuan et al., 2018b), the only HDBR results reported here are the group differences between HDBR + CO₂ and HDBR.

SANS Versus No-SANS Group Comparisons

We performed two exploratory analyses to examine group differences between those HDBR + CO₂ subjects who developed signs of SANS (SANS; n = 5; two males, three females) and those who did not (no-SANS; n = 6; four males, two females). First, we tested for differences between the intercept and slope images for each group. We conducted two-sample parametric t-tests with threshold p < 0.0005, k = 10. Non-parametric testing would not have been possible here, as less than 500 permutations exist for this combination of sample sizes.

Next, we tested for regions where the SANS versus no-SANS groups showed differences in the correlation between pre- to post-HDBR + CO₂ brain change and pre- to post- change in the ratio between the balance—eyes open and balance—eyes closed condition. This ratio score was calculated as: (balance—eyes open score/balance—eyes closed score) ^∗100 and provides a metric of the degree to which an individual relies on vision for maintaining quiet upright stance. Each of these balance tasks is described in Section “Balance Testing.” We selected to compare brain-behavior correlations only for this ratio score here because we previously identified significant differences between the two SANS subgroups on this measure, in which SANS individuals showed greater increases from pre- to post-HDBR + CO₂ in their reliance on vision during balance compared to no-SANS individuals (Lee et al., 2019a). One of the five SANS subjects was excluded from this analysis due to outlier values for the balance—eyes closed condition (described in Section “Balance Testing”). Thus there were n = 4 for the SANS group and n = 6 for the no-SANS group.

Mobility and Balance Testing: HDBR + CO₂ Cohort Only

Although the HDBR + CO₂ participants completed a battery of neurocognitive and sensorimotor assessments at each time point, here we focus on only mobility and balance testing, as these tasks were the most directly related to vestibular processing. We have previously published comprehensive behavioral profiles for both the HDBR + CO₂ (Lee et al., 2019a) and HDBR groups (Koppelmans et al., 2015), as well as vestibular brain-behavior correlations for the HDBR group (Yuan et al., 2018b).

Functional Mobility Test (FMT)

The Functional Mobility Test (FMT) is sensitive to the effects of spaceflight (Mulavara et al., 2010) and to the effects of bed rest (Reschke et al., 2009; Koppelmans et al., 2017). The FMT requires subjects to arise from a seated position and walk through a 6-m × 4-m two-part obstacle course consisting of foam hurdles, pylons, and bars. The first part of the course was completed on a hard floor, and the second part was completed on medium-density foam. Participants were instructed to walk through the course as quickly as possible without touching any of the obstacles. Participants repeated the FMT 10 times per session on five different testing days (Figure 1). Here we analyze only the total time needed to complete the course for the first trial of each session, as we have found this measure to be the most sensitive to intervention-related change. We excluded one subject from FMT analyses, as the subject showed substantial pre- to post- slowing (> ± 2.5 standard deviations from the group average pre- to post- change) and exerted considerable influence on group-level statistics. Thus there were n = 10 subjects for analyses involving FMT.

Balance Testing

Participants completed three balance tasks: (1) balance—eyes open; (2) balance—eyes closed; and (3) balance—eyes closed dynamic head tilt. Details of these tasks have been previously described (Mulder et al., 2014). Participants stood on a foam pad on top of a force platform (Leonardo Mechanograph, Novotec Medical GmbH, Pforzheim, Germany). Participants were instructed to maintain a comfortable stance, keep their arms folded across their chest and remain in a stable, upright posture for 30 s. Foot markers on the foam pad were used to ensure consistent foot placement across trials and between subjects. For the first two conditions, participants kept their head erect and eyes either open or closed. For the eyes closed dynamic head tilt condition, participants kept their eyes shut and made head pitch motions of ±20°, synchronized to a 0.33 Hz metronome tone. Participants repeated all conditions three times during each testing session, and the order of conditions was semi-randomized to ensure that identical conditions did not repeat back-to-back. To minimize the effect of outlier trials, for each condition we examined the median score of the three trials. Scores are reported as equilibrium quotients (EQ), where 100% is a perfect score. EQ scores were calculated using instantaneous anterior-to-posterior peak-to-peak center-of-mass sway angle. For the eyes closed condition only, we excluded one subject, as the subject showed pre- to post-HDBR + CO₂ declines > ± 2 standard deviations from the group average and exerted considerable influence on group-level statistics. Thus there were n = 10 subjects for analyses involving balance-eyes closed scores, but there were n = 11 subjects for all other balance analyses.

Statistical Analyses of Behavioral Data

For completeness, we tested pre- to post- HDBR + CO₂ behavioral change, and we tested recovery for the mobility and balance tasks. In R 3.5.1 (R Core Team, 2013), using the last pre- bed rest time point (BDC 7) and the end of bed rest time point (HDT 29), we calculated a slope of pre- to post- performance change for each subject and conducted a one-sample t-test to determine if the group-average slope was different from 0. We also examined post-HDBR + CO₂ recovery trajectories for the three post-bed rest time points using a linear mixed model with restricted maximum likelihood (REML) estimation via the “lme” function. The model included a random intercept for subject (to allow for different starting points for each person) and the fixed effect of time. In each case, we were interested in whether the fixed effect of time was significant; we tested a quadratic fit for time for each measure as well, but the model including the linear effect of time performed better in all cases.

Results

Neural Response to Vestibular Stimulation

Average BOLD signal during vestibular stimulation versus rest across all subjects and all time points is shown in Figure 3 to illustrate the neural response to the skull tap method. In line with previous work (Lopez et al., 2012b; Zu Eulenburg et al., 2012; Noohi et al., 2017; Yuan et al., 2018b) vestibular stimulation resulted in activation of clusters in the right and left insula (Table 1). Also, as anticipated, we observed widespread deactivation of frontal, temporal, occipital, subcortical, and cerebellar regions. These results demonstrate that our skull tap method was able to engage the vestibular system and produce the expected neural response.

FIGURE 3

Figure 3. Neural response to vestibular stimulation. Vestibular stimulation resulted in activation of insular cortex and widespread deactivation including (A) frontal, temporal, occipital, subcortical, and (B) cerebellar regions. Whole brain and cerebellar results are overlaid onto MNI (A) and SUIT (B) standard templates, respectively; p < 0.0005, k = 10; red = regions of activation; blue = regions of deactivation.

TABLE 1

Table 1. Regions of activation or deactivation in response to vestibular stimulation.

Time Course of Neural Vestibular Response to HDBR + CO₂

We identified multiple longitudinal changes in the neural response to vestibular stimulation across HDBR + CO₂, followed by recovery (Figure 4 and Table 2). Several frontal, parietal, and temporal regions (Figure 4A) showed immediate decreases in activation with HDBR + CO₂, followed by recovery during the intervention and complete recovery by the final bed rest time point, HDT 29. For instance, right inferior temporal gyrus showed a conversion from activation to deactivation with the onset of HDBR + CO₂, followed by recovery of activation of this region by HDT 29.

FIGURE 4

Figure 4. Time course of neural vestibular response to HDBR + CO₂. Left panels: Regions showing whole brain (A,B) and cerebellar (C) longitudinal changes in activation during vestibular stimulation across the six time points. Whole brain and cerebellar results are overlaid onto MNI and SUIT standard templates, respectively; p < 0.001, k = 10 for whole brain analyses; k = 5 for cerebellar analyses. Three clusters survived p < 0.0005, k = 10 correction (see Table 2). Right panels: Example contrast values plotted for peak coordinate within the cluster with the largest T value in each case (cluster indicated by red arrow). Points represent group mean contrast values; error bars represent standard error. Dotted lines depict the hypothesized longitudinal contrasts for “instant decrease, during-HDBR + CO₂ recovery” (A), “instant decrease, post-HDBR + CO₂ recovery” (B), and “instant increase, post-HDBR + CO₂ recovery” (B,C).

TABLE 2

Table 2. Regions showing longitudinal increases and decreases in activation during vestibular stimulation across all six time points.

Several other regions showed patterns of fast change, with changes sustaining throughout HDBR + CO₂ and not restoring until after the conclusion of bed rest. Right superior medial gyrus and right cerebellar lobule VI both showed decreases in deactivation and a conversion to activation with the start of bed rest, followed by recovery by 12 days post- bed rest. While neither of these clusters precisely overlaps with the regions that deactivated on average during vestibular stimulation (Figure 3 and Table 1), other nearby parts of the superior medial gyrus and right cerebellar lobule VI did significantly deactivate in response to vestibular stimulation.

One brainstem cluster showed a fast decrease in activation with bed rest, with a conversion to deactivation of this region during HDBR + CO₂, followed by recovery. Of note, we did not find differences between the brain regions that emerged as significant for the “instant” versus “cumulative change post-HDBR + CO₂ recovery” contrasts, so we have reported only results for the instant change post-HDBR + CO₂ recovery contrasts.

Functional Behavioral Implications

Mobility and Balance Changes With HDBR + CO₂

Subjects showed pre- to post-HDBR + CO₂ declines in mobility, followed by a linear recovery pattern (Figure 5 and Table 3); that is, participants were slower to complete the FMT obstacle course post-HDBR + CO₂, but sped back up by 12 days post-bed rest. The slope of decline in balance scores was only significant for the balance—eyes open condition; however, visually (Figure 5), there was a clear trend that HDBR + CO₂ negatively impacted balance across all three tasks. Similarly, only the balance—eyes closed dynamic head tilt condition showed a significant linear recovery pattern post-HDBR + CO₂, but again, visually, a recovery trend was evident post-bed rest for each of the balance tasks.

FIGURE 5

Figure 5. Mobility and Balance Changes with HDBR + CO₂. Mobility and balance performance pre- and post-HDBR + CO₂. In general, participants showed performance declines followed by recovery post-HDBR + CO₂. The slope of behavioral decline was significant for FMT and a trend for the balance-eyes open condition. Participants showed a significant linear recovery trajectory for FMT and balance-eyes closed dynamic head tilt. Error bars represent standard error. ^∗p < 0.10; ^∗∗p < 0.05; ^∗∗∗p < 0.001.

TABLE 3

Table 3. Mobility and balance change with HDBR + CO₂ and recovery.

Brain—Behavior Correlations

We identified several dozen regions for which pre- to post-HDBR + CO₂ change in neural response to vestibular stimulation correlated with pre- to post- change in mobility and balance performance (Figures 6A–D and Table 4). In general, across all tasks and almost all clusters, we found that greater deactivation of various brain regions was associated with reduced decline or even improvement in behavioral measures. For instance, for the balance—eyes open condition, we found that greater deactivation of the right superior temporal gyrus (Figure 6B) associated with less balance decline and even balance improvement for a few individuals. That is, those with the greatest decreases in activation or increases in deactivation of this region had the best post-HDBR + CO₂ balance performance. Similarly, for the eyes closed and eyes closed dynamic head tilt balance conditions, we found that greater pre- to post- deactivation of right cerebellar lobule I-IV and supplementary motor area, respectively, was associated with less balance decline or even balance improvement (Figures 6C,D).

FIGURE 6

Figure 6. Functional behavioral consequences. Brain—behavior correlations for functional mobility test (FMT; A) and three balance tasks (B–D). Brain change was calculated as the difference in brain activation from pre-HDBR + CO₂ (time BDC 7) to the end of HDBR + CO₂ (time HDT 29). Behavior change was calculated as the difference in mobility and balance performance from pre-HDBR + CO₂ (time BDC 7) to post-HDBR + CO₂ (time R0). Whole brain and cerebellar results are overlaid onto MNI and SUIT standard templates, respectively; non-parametric p < 0.0005, k = 10. Warm colors indicate regions of positive correlation between brain change and behavior change. Cool colors indicate regions of negative correlation between brain change and behavior change. Right side correlation plots include contrast values extracted from peak coordinate inside an example cluster (indicated with red arrows), graphed against behavior change score.

TABLE 4

Table 4. Regions showing associations between pre- to post-HDBR + CO₂ differences in behavioral scores and pre- to post-HDBR + CO₂ change in brain activation during vestibular stimulation.

This relationship between greater deactivation and reduced behavioral decline held across the majority of brain regions that showed significant brain—behavior correlations, including sensorimotor cortex (i.e., supplementary motor area and postcentral gyrus), temporal cortex, occipital cortex, brainstem, and cerebellum (particularly, anterior cerebellum and crus I and II). In several cases, increased deactivation was found in regions that typically deactivate during vestibular stimulation (Figure 3 and Table 1); for instance, this was the case for supplementary motor area, postcentral gyrus, and occipital gyrus (i.e., clusters marked with superscript “d” in Table 4).

There were only a few regions where reduced pre- to post- deactivation was more beneficial for post- HDBR + CO₂ mobility and balance performance. For instance, in the case of right cerebellar lobule VIIb (Figure 6A), decreased deactivation from pre- to post- was associated with less FMT slowing (i.e., less decline). Similarly, for the balance—eyes closed dynamic head tilt condition, less deactivation from pre- to post- in left cerebellar crus II was associated with less balance decline.

HDBR + CO₂ vs. HDBR Group Comparisons

Baseline (Intercept) Differences

Only three regions emerged where HDBR + CO₂ and HDBR had baseline differences in neural response to vestibular stimulation: left inferior temporal gyrus, right superior occipital gyrus, and brainstem (Figure 7A and Table 5). That is, both groups produced similar neural responses to vestibular stimulation pre- bed rest, and thus between-group slope differences can likely be attributed to intervention effects.

FIGURE 7

Figure 7. HDBR + CO₂ vs. HDBR group comparisons. (A) Regions of intercept differences between HDBR + CO₂ and HDBR groups. (B) Regions of slope differences between HDBR + CO₂ and HDBR groups. Results overlaid onto MNI standard template; non-parametric p < 0.0005, k = 10. Cool colors indicate regions where the intercept or slope for the HDBR + CO₂ group was numerically less than the intercept or slope for the HDBR group. Warm colors indicate regions where the slope for the HDBR + CO₂ group was numerically greater than the slope for the HDBR group. (C) Example slope values extracted from L Middle Frontal Gyrus (i.e., the cluster with greatest T value for the HDBR + CO₂ > HDBR contrast). Here, the HDBR + CO₂ group showed increased activation of this region across bed rest, while the HDBR group showed increased deactivation of this region across bed rest.

TABLE 5

Table 5. Regions with intercept differences between HDBR + CO₂ and HDBR subjects.

Slope Differences

There were five clusters across frontal, parietal, and temporal cortex where the HDBR + CO₂ group had a numerically greater slope of change in neural response to vestibular stimulation across the course of bed rest (Figure 7B and Table 6). Although no regions here overlapped with brain areas from the main effect analysis (Figure 3 and Table 1), these clusters were located in close proximity to regions that are expected to show deactivation during vestibular stimulation. In general, the HDBR + CO₂ group showed increases in activation of these regions over the course of bed rest, as well as more within-group variability in neural response, compared to the HDBR group who generally showed increases in deactivation of these regions over the course of bed rest. For instance, in the left middle frontal gyrus (Figure 7C), the HDBR + CO₂ subjects showed a switch from deactivation of this region during vestibular stimulation at BDC 7 to activation of this region at HDT 29. The HDBR group showed the opposite pattern, changing from activation to deactivation of this region.

TABLE 6

Table 6. Regions with differences in slope of activation change with bed rest between HDBR + CO₂ and HDBR subjects.

There was one cluster in the thalamus where the HDBR + CO₂ group showed a reduced slope of change compared to the HDBR group. Here the HDBR + CO₂ group exhibited decreasing activation, and ultimately deactivation of this region during vestibular stimulation at HDT 29, whereas the HDBR group exhibited a transition from deactivation of this region to activation.

SANS vs. No-SANS Group Differences

Five of the 11 HDBR + CO₂ participants developed signs of SANS, including optic disc edema. While this phenomenon is commonly reported following spaceflight (Lee et al., 2016), this is the first bed rest study to induce such effects (Laurie et al., 2019), possibly due to the careful testing for SANS symptoms and the strict head-down-tilt conditions, or the addition of elevated CO₂. As SANS was not anticipated a priori but represents a substantial subgroup of the HDBR + CO₂ cohort, we conducted two exploratory analyses of this unique sample.

At p < 0.0005 and k = 10, there were no regions of intercept or slope difference between the SANS and no-SANS participants. We previously found that SANS individuals increased their reliance on visual information during balance from pre- to post-HDBR + CO₂ (Lee et al., 2019a); that is, SANS subjects had greater increases in their ratio of eyes open balance compared to eyes closed balance. Here we identified multiple frontal, parietal, temporal, and occipital regions where the SANS subjects showed stronger correlations with this balance ratio score compared to the no-SANS subjects (Figure 8 and Table 7). For instance, SANS subjects showed a stronger correlation between greater activation of left middle frontal gyrus (Figure 8) and increased pre- to post- balance ratio score. Several of these clusters included regions typically activated during vestibular stimulation (indicated with a superscript “d” in Table 7). There were no regions of stronger correlation for the no-SANS subjects.

FIGURE 8

Figure 8. SANS vs. No-SANS Group Differences. Left: regions in which HDBR + CO₂ participants who developed signs of Spaceflight-Associated Neuroccular Syndrome (SANS) showed a stronger correlation between pre- to post- change in brain activity during vestibular stimulation and pre- to post- change in balance ratio compared to those who did not develop signs of SANS (no-SANS). Balance ratio was calculated as: (balance—eyes open score/balance—eyes closed score) ^∗ 100. Results overlaid onto MNI standard template; non-parametric p < 0.0005, k = 10. One of the five SANS subjects was excluded from analysis here due to outlier values for their eyes closed balance score; thus n = 4 for the SANS group and n = 6 for the no-SANS group. Right: correlation between brain changes and visual dependence balance ratio change for an example cluster; SANS subjects are shown in pink, no-SANS subjects are shown in purple. Regression lines are shown for each group for reference only and do not represent the results of any statistical tests.

TABLE 7

Table 7. Regions where SANS subjects showed greater correlations compared to no-SANS subjects between pre- to post HDBR + CO₂ change in visual dependence balance ratio and pre- to post-HDBR + CO₂ change in brain activation during vestibular stimulation.

Discussion

Key Findings

Here we identified changes in the neural correlates of vestibular processing with 30 days of HDBR + CO₂. We found multiple regions where brain activation during vestibular stimulation changed quickly after participants started HDBR + CO₂ and recovered either during or post-bed rest, providing support for adaptive plasticity of the vestibular system in response to altered sensory inputs. In multiple cases, increased deactivation of cortical and cerebellar areas was associated with less decline in balance from pre- to post-HDBR + CO₂, suggesting that some of the adaptive neural changes during bed rest may benefit post-bed rest performance of vestibularly mediated behaviors. We found several differences for HDBR + CO₂ compared to HDBR subjects, suggesting interactive or additive effects of bed rest and CO₂. Finally, we noted differences in brain—behavior relationships for SANS versus no-SANS subjects, indicating the need for further study of bed rest-induced ocular symptoms.

Time Course of Neural Vestibular Response to HDBR + CO₂

Similar to our past work (Yuan et al., 2018b), we found multiple longitudinal changes in the neural response to vestibular stimulation, including in several areas in close proximity to regions typically involved in processing vestibular information, as well as in several regions that are not normally activated during vestibular stimulation. These responses could represent adaptive plasticity, in which the enhanced demands of neural processing of altered sensory inputs during HDBR + CO₂ are requiring greater neural resources. For instance, the finding of decreased deactivation in right superior medial gyrus and right cerebellar lobule VI (which are anatomically near to regions that, on average, deactivated in response to vestibular stimulation) suggests a compensatory response. That is, functional brain regions that typically deactivate in response to vestibular input are deactivating less during exposure to an altered sensory environment, potentially to allow for additional brain pathways to aid in processing the novel sensory information. Here, more specifically, it could be that down-weighting of somatosensory input during HDBR, paired with upweighting of vestibular input due to vestibular-somatosensory convergence at the vestibular nuclei (Mulavara et al., 2012), is resulting in a higher neural processing demand. The fast recovery of several of these regions during bed rest suggests an ability of the vestibular system to adjust rapidly to such altered sensory conditions.

In contrast to our past work (Yuan et al., 2018b), we did not identify any unique regions of slow, cumulative brain changes; instead, we found clusters that were statistically significant for both the immediate and cumulative models. This suggests that interactive effects of CO₂ with bed rest might accelerate neural vestibular changes. It could also be that CO₂-related increases in cerebral perfusion enhanced the BOLD signal (Corfield et al., 2001) for HDBR + CO₂ participants, making it easier to detect bed rest-related changes earlier during that intervention. Further, given that the HDBR + CO₂ intervention was about half the duration of the HDBR intervention, it could be that some of the slow, cumulative brain changes that we previously identified require longer than 30 days to develop.

Previously we found that HDBR resulted in upregulation of the vestibular system (Yuan et al., 2018b), which we attributed to either increased sensitivity of the vestibular system during HDBR or to reduced neural efficiency, in which greater activation of vestibular cortical regions would be needed to process vestibular information during HDBR. Here, we did not find clear evidence for reduced neural efficiency (Yuan et al., 2018b). That is, we did not identify any HDBR + CO₂-related increases in activation of vestibular cortical regions. This suggests that elevated CO₂ may augment vestibular processing, as this environment did not produce the same longitudinal reduction in neural efficiency as HDBR alone.

Functional Behavioral Implications

We identified predominantly regions for which increased pre- to post- deactivation during vestibular stimulation associated with reduced balance performance decline, or even performance improvement. This may represent adaptive plasticity during the HDBR + CO₂ intervention, in which some individuals show an enhancement of the expected cortical deactivation response, paired with a dampening of activity in other brain regions that could interfere with processing of vestibular information. This adaptive change could then later manifest as superior post-bed rest balance due to underlying increased specificity of activation of vestibular cortex and deactivation of other sensory regions during the balance tasks. More specifically, it could be that the reduced plantar somatosensory input during HDBR + CO₂ results in down-weighting of somatosensory input, but upweighting of vestibular input. While this reweighting is likely modulated at the level of the vestibular nuclei where the somatosensory and vestibular systems converge (Bles et al., 1984; Dieringer, 1995; Horak and Hlavacka, 2001; Carriot et al., 2015), this reorganization could plausibly manifest as increased deactivation of cortical sensorimotor processing regions, with those individuals who had the most successful reweighting processes presenting with enhanced preservation of balance abilities post- HDBR + CO₂. Given these possible mechanisms, it thus makes sense that we found the most numerous brain-behavior correlations for the eyes closed dynamic head tilt condition, as this condition most specifically tasks the vestibular system.

We similarly found evidence of brain-behavioral relationships in our past HDBR work for vestibular processing (Yuan et al., 2018b) and for neural control of foot movement (Yuan et al., 2018a). Further, we previously identified several regions for which better balance (i.e., reduced postural sway while standing on one leg) correlated with greater deactivation of the brainstem, cerebellar lobule VI, and crus I and II across healthy young and older adults at one time point (Noohi et al., 2019). This fits with the present work, as here we found several brainstem and cerebellar regions for which greater post-bed rest deactivation associated with better balance. Together, these findings provide further support for the notion that those with a more refined neural response to vestibular stimulation in the scanner (e.g., including greater deactivation of brainstem and cerebellum) likely also produce a more refined neural response during balance tasks outside of the scanner and therefore perform better.

Similarly, several mobile neuroimaging studies have identified that older adults exhibit poorer balance paired with increased brain activation during vestibular stimulation (Karim et al., 2013; Lin et al., 2017). For instance, using functional near-infrared spectroscopy (fNIRS), Lin et al. (2017) found greater activation in frontal and occipital regions during vestibular stimulation for older compared to middle-aged adults. This suggests compensatory processes in which older adults require greater neural resources to process the same vestibular information. This fits with the present findings, as those who had the largest post-bed rest balance declines also showed bed rest-related increases in activation or reductions in deactivation across various cortical regions, including frontal and occipital regions. It could be that, similar to older adults, these individuals were recruiting extra brain regions to aid in processing vestibular information in the scanner and then engaging similar compensatory over-recruitment mechanisms outside the scanner during the post-bed rest balance assessments, which ultimately resulted in poorer balance performance.

HDBR + CO₂ vs. HDBR Group Comparisons

CO₂-specific effects or interactive effects of bed rest and CO₂ may be contributing to the identified differences in slope of activation change for the HDBR + CO₂ group versus the HDBR group. CO₂ is a strong vasodilator and, among other effects, results in increased blood flow to the brain (Atkinson et al., 1990; Zhou et al., 2008), as well as increased intensity of the blood oxygen level-dependent (BOLD) signal measured by fMRI (Corfield et al., 2001). This increased cerebral perfusion could be contributing to the identified group differences here. Elevated CO₂ selectively favors frontal lobe perfusion (Bhogal et al., 2015); thus increased slope of change for frontal regions (e.g., middle frontal gyrus) among the HDBR + CO₂ group could be particularly related to perfusion effects on the BOLD signal. As brain deactivation is also an active process, the widespread increases in deactivation that we noted (which were generally associated with better behavioral performance) could also be influenced by CO₂-related perfusion effects.

It has been demonstrated that decreased CO₂ (through voluntary hyperventilation) negatively impacts postural sway, resulting in unsteadiness of balance (Sakellari et al., 1997). Although the mechanisms here are unknown, it has been suggested that hyperventilation disrupts vestibular system compensation, including interfering with central and peripheral somatosensory signals from the lower limbs (Sakellari et al., 1997). Here the differential effects between HDBR + CO₂ and HDBR subjects demonstrate that increased CO₂ may also disrupt normal vestibular processing and vestibular compensatory mechanisms.

SANS vs. No-SANS Group Differences

SANS subjects showed stronger correlations between pre- to post-HDBR + CO₂ change in balance ratio score and brain changes. This suggests a relationship between SANS status and visual contributions to balance. Higher ratio scores are associated with more reliance on external visual cues for balance. As those with the greatest activation of regions such as middle frontal gyrus also showed the greatest pre- to post- increase in visual dependence, this could indicate reduced neural efficiency—in which these individuals are recruiting extra brain regions to aid in processing of the same vestibular information. This indicates reduced efficiency of the vestibular system that manifested behaviorally as greater reliance on the visual system during post-bed rest balance. These findings require validation in future studies to more clearly understand implications for the one third of astronauts who develop SANS (Lee et al., 2016).

Limitations

Limitations of the present work include the small pilot sample, lack of ambulatory control group, and differences in scanning timeline and parameters between the two groups. As the HDBR + CO₂ and HDBR groups were each part of separate, larger bed rest campaigns, the respective testing timelines were restricted by NASA and not matched between groups. Further, only males were included in the HDBR group, as these participants were control subjects for a different investigator’s testosterone supplementation study. Although the vestibular stimulation parameters and sequence duration were identical between the HDBR + CO₂ and HDBR groups, these data were collected on two different Siemens scanners with two slightly different fMRI sequences. The HDBR + CO₂ fMRI sequence included a faster TR and more volumes (TR = 2.5 s; 96 volumes) than the HDBR sequence (TR = 3.66 s; 66 volumes). This amounts to greater statistical power for the HDBR + CO₂ group. While these differences represent a limitation of the between-group comparisons that we report here, multisite neuroimaging studies are increasing in popularity; such studies generally indicate that functional neuroimaging data is robust across sites (Costafreda et al., 2007; Biswal et al., 2010; Gountouna et al., 2010; Noble et al., 2017). Further, the analyses in the present work focus on differences in within-subject brain changes. That is, we have tested for between-group differences in within-person intercepts and slopes of change in brain activity. This makes the present results more robust to any introduced variance due to scanner or sequence differences.

Due to the limited pilot sample size, we use uncorrected p-values for the neuroimaging statistical tests to better detect within- and between-subject differences. Nonetheless, here we demonstrate for the first time the feasibility of characterizing vestibular brain changes during multiple weeks of bed rest combined with elevated CO₂. We present compelling preliminary findings based on these highly unique data, which should be validated in future work.

It should also be noted that it is difficult to fully generalize these bed rest findings to spaceflight. HDBR + CO₂ and HDBR mimic only some of the effects of spaceflight, such as body unloading, altered sensory input, and fluid shifts, but these analogs do not include all features of spaceflight that could impact vestibular processing. Finally, as is typical with the pneumatic skull tap method, the stimulus did not induce any vestibular perception, motion sensation, or movement in head position. Instead of using subjective perception of the vestibular stimulus, assessment of vestibular-evoked myogenic potentials in the eye muscles (oVEMPs) outside of the scanner was used to validate successful stimulation of vestibular organs.

Applications to Spaceflight and Future Directions

The present findings document changes in functional vestibular processing with spaceflight analog environments. These findings support that spaceflight factors do likely influence the neural correlates of vestibular processing; however, there is limited past work investigating vestibular processing in astronauts. Two previous studies suggest changes from pre- to post- spaceflight in resting-state (Demertzi et al., 2016) and task-based connectivity (Pechenkova et al., 2019) in brain networks that support vestibular function. Our past work has identified disrupted white matter structural connectivity in several tracts that underlie sensory integration and vestibular processes and associations of these brain changes with balance declines (Lee et al., 2019b). For instance, we found that astronauts with the largest spaceflight-associated balance disruptions also had the greatest white matter declines in the superior longitudinal fasciculus (Lee et al., 2019b), which connects the temporoparietal and prefrontal cortices and is thought to subserve vestibular functions (Spena et al., 2006). However, no studies to date have tested spaceflight-related changes in functional brain activity during processing of vestibular information. To address this critical literature gap, an ongoing prospective study by our group is measuring brain activity with the pneumatic skull tap paradigm at two time points before and four time points after astronauts complete ISS missions. Comparing these results with the present work will help to elucidate how additional microgravity factors not induced by bed rest (e.g., an altered gravitational vector) might affect the neural correlates of vestibular processing.

Conclusion

Here we demonstrate the feasibility of assessing longitudinal neural vestibular changes following 30 days of HDBR + CO₂. We identify support for specific effects of combined HDBR + CO₂ on vestibular processing, adaptive plasticity of the vestibular system during HDBR + CO₂ followed by fast and slow recovery, and relationships between adaptive plasticity and spared behavioral performance post-HDBR + CO₂. We note some differences between neural processing of vestibular information for HDBR + CO₂ versus HDBR subjects, as well as implications for dependence on visual cues during balance for SANs versus no-SANS subjects. Together, these findings contribute to understanding of how the vestibular system adapts to altered sensory inputs and to understanding of how spaceflight may influence the neural correlates of vestibular processing.

Data Availability Statement

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

Ethics Statement

The studies involving human participants were reviewed and approved by the local ethical commission of the regional medical association, Ärztekammer Nordrhein, as well as the University of Florida and NASA Institutional Review Boards. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

KH analyzed the vestibular fMRI data and mobility/balance behavioral data, created the figures and tables, and wrote the manuscript. JL collected and managed the data and participated in manuscript preparation. NG collected and analyzed the data. IK participated in project design and software development. YD collected and analyzed the data. JB, AM, and RS designed the project, secured funding, and led the interpretation and discussion of the results. All authors participated in revision of the manuscript.

Funding

This work was supported by grants to RS, AM, and JB from the National Aeronautics and Space Administration (NASA; 80NSSC17K0021), the National Space Biomedical Research Institute (NSBRI; NCC 9-58 and SA02802), and the National Institutes of Health and National Center for Advancing Translational Sciences (1UL1RR029876-01). Additionally, during completion of this work KH was supported by a National Science Foundation Graduate Research Fellowship under Grant Nos. DGE-1315138 and DGE-1842473, as well as training grant NIH T32-NS082128.

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.

The reviewer RG declared a past co-authorship with several of the authors AM and JB to the handling Editor.

Acknowledgments

The authors wish to gratefully acknowledge those at the: envihab facility in Cologne, Germany, including Dr. Edwin Mulder, the VaPER study staff, and the bed rest participants who volunteered their time, without whom this project would not have been possible.

Footnotes

1. ^ www.nitrc.org/projects/artifact_detect/
2. ^ http://warwick.ac.uk/snpm

References

Ashburner, J., Barnes, G., Chen, C., Daunizeau, J., Flandin, G., Friston, K., et al. (2016). SPM12 Manual. Available at: http://web.mit.edu/spm_v12/manual.pdf (accessed September 2, 2019).

Google Scholar

Assländer, L., and Peterka, R. J. (2014). Sensory reweighting dynamics in human postural control. J. Neurophysiol. 111, 1852–1864. doi: 10.1152/jn.00669.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Atkinson, J. L., Anderson, R. E., and Sundt, T. M. Jr. (1990). The effect of carbon dioxide on the diameter of brain capillaries. Brain Res. 517, 333–340. doi: 10.1016/0006-8993(90)91046-j

PubMed Abstract | CrossRef Full Text | Google Scholar

Bhogal, A. A., Philippens, M. E., Siero, J. C., Fisher, J. A., Petersen, E. T., Luijten, P. R., et al. (2015). Examining the regional and cerebral depth-dependent BOLD cerebrovascular reactivity response at 7 T. Neuroimage 114, 239–248. doi: 10.1016/j.neuroimage.2015.04.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Biswal, B. B., Mennes, M., Zuo, X.-N., Gohel, S., Kelly, C., Smith, S. M., et al. (2010). Toward discovery science of human brain function. Proc. Natl. Acad. Sci. U.S.A. 107, 4734–4739. doi: 10.1073/pnas.0911855107

PubMed Abstract | CrossRef Full Text | Google Scholar

Bles, W., Jong, J. M. V., and Wit, G. D. (1984). Somatosensory compensation for loss of labyrinthine function. Acta Otolaryngol. 97, 213–221. doi: 10.3109/00016488409130982

PubMed Abstract | CrossRef Full Text | Google Scholar

Boyle, R., Mensinger, A. F., Yoshida, K., Usui, S., Intravaia, A., Tricas, T., et al. (2001). Neural readaptation to Earth’s gravity following return from space. J. Neurophysiol. 86, 2118–2122. doi: 10.1152/jn.2001.86.4.2118

PubMed Abstract | CrossRef Full Text | Google Scholar

Carriot, J., Jamali, M., and Cullen, K. E. (2015). Rapid adaptation of multisensory integration in vestibular pathways. Front. Syst. Neurosci. 9:59. doi: 10.3389/fnsys.2015.00059

PubMed Abstract | CrossRef Full Text | Google Scholar

Cassady, K., Koppelmans, V., Reuter-Lorenz, P., De Dios, Y., Gadd, N., Wood, S., et al. (2016). Effects of a spaceflight analog environment on brain connectivity and behavior. Neuroimage 141, 18–30. doi: 10.1016/j.neuroimage.2016.07.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen, H. S., Kimball, K. T., Mulavara, A. P., Bloomberg, J. J., and Paloski, W. H. (2012). Posturography and locomotor tests of dynamic balance after long-duration spaceflight. J. Vestib. Res. 22, 191–196. doi: 10.3233/VES-2012-0456

PubMed Abstract | CrossRef Full Text | Google Scholar

Corfield, D., Murphy, K., Josephs, O., Adams, L., and Turner, R. (2001). Does hypercapnia-induced cerebral vasodilation modulate the hemodynamic response to neural activation? Neuroimage 13, 1207–1211. doi: 10.1006/nimg.2001.0760

PubMed Abstract | CrossRef Full Text | Google Scholar

Costafreda, S. G., Brammer, M. J., Vêncio, R. Z., Mourao, M. L., Portela, L. A., De Castro, C. C., et al. (2007). Multisite fMRI reproducibility of a motor task using identical MR systems. J. Magn. Reson. Imaging 26, 1122–1126. doi: 10.1002/jmri.21118

PubMed Abstract | CrossRef Full Text | Google Scholar

Demertzi, A., Van Ombergen, A., Tomilovskaya, E., Jeurissen, B., Pechenkova, E., Di Perri, C., et al. (2016). Cortical reorganization in an astronaut’s brain after long-duration spaceflight. Brain Struct. Funct. 221, 2873–2876. doi: 10.1007/s00429-015-1054-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Diedrichsen, J. (2006). A spatially unbiased atlas template of the human cerebellum. Neuroimage 33, 127–138. doi: 10.1016/j.neuroimage.2006.05.056

PubMed Abstract | CrossRef Full Text | Google Scholar

Diedrichsen, J., Balsters, J. H., Flavell, J., Cussans, E., and Ramnani, N. (2009). A probabilistic MR atlas of the human cerebellum. Neuroimage 46, 39–46. doi: 10.1016/j.neuroimage.2009.01.045

PubMed Abstract | CrossRef Full Text | Google Scholar

Dieringer, N. (1995). ‘Vestibular compensation’: neural plasticity and its relations to functional recovery after labyrinthine lesions in frogs and other vertebrates. Prog. Neurobiol. 46, 97–129. doi: 10.1016/0301-0082(94)00063-n

CrossRef Full Text | Google Scholar

Fredrickson, J., Schwarz, D., and Kornhuber, H. (1966). Convergence and interaction of vestibular and deep somatic afferents upon neurons in the vestibular nuclei of the cat. Acta Otolaryngol. 61, 168–188. doi: 10.3109/00016486609127054

PubMed Abstract | CrossRef Full Text | Google Scholar

Gläscher, J., and Gitelman, D. (2008). Contrast Weights in Flexible Factorial Design with Multiple Groups of Subjects (London: SPM), 1–12.

Google Scholar

Gountouna, V.-E., Job, D. E., Mcintosh, A. M., Moorhead, T. W. J., Lymer, G. K. L., Whalley, H. C., et al. (2010). Functional magnetic resonance imaging (fMRI) reproducibility and variance components across visits and scanning sites with a finger tapping task. Neuroimage 49, 552–560. doi: 10.1016/j.neuroimage.2009.07.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Hargens, A. R., and Richardson, S. (2009). Cardiovascular adaptations, fluid shifts, and countermeasures related to space flight. Respir. Physiol. Neurobiol. 169, S30–S33. doi: 10.1016/j.resp.2009.07.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Horak, F., and Hlavacka, F. (2001). Somatosensory loss increases vestibulospinal sensitivity. J. Neurophysiol. 86, 575–585. doi: 10.1152/jn.2001.86.2.575

PubMed Abstract | CrossRef Full Text | Google Scholar

Jian, B., Shintani, T., Emanuel, B., and Yates, B. (2002). Convergence of limb, visceral, and vertical semicircular canal or otolith inputs onto vestibular nucleus neurons. Exp. Brain Res. 144, 247–257. doi: 10.1007/s00221-002-1042-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Karim, H., Fuhrman, S. I., Furman, J. M., and Huppert, T. J. (2013). Neuroimaging to detect cortical projection of vestibular response to caloric stimulation in young and older adults using functional near-infrared spectroscopy (fNIRS). Neuroimage 76, 1–10. doi: 10.1016/j.neuroimage.2013.02.061

PubMed Abstract | CrossRef Full Text | Google Scholar

Koppelmans, V., Bloomberg, J. J., De Dios, Y. E., Wood, S. J., Reuter-Lorenz, P. A., Kofman, I. S., et al. (2017). Brain plasticity and sensorimotor deterioration as a function of 70 days head down tilt bed rest. PLoS One 12:e0182236. doi: 10.1371/journal.pone.0182236

PubMed Abstract | CrossRef Full Text | Google Scholar

Koppelmans, V., Mulavara, A. P., Yuan, P., Cassady, K. E., Cooke, K. A., Wood, S. J., et al. (2015). Exercise as potential countermeasure for the effects of 70 days of bed rest on cognitive and sensorimotor performance. Front. Syst. Neurosci. 9:121. doi: 10.3389/fnsys.2015.00121

PubMed Abstract | CrossRef Full Text | Google Scholar

Kurth, F., Luders, E., and Gaser, C. (2010). VBM8 Toolbox Manual. Jena: University of Jena.

Google Scholar

Laurie, S. S., Macias, B. R., Dunn, J. T., Young, M., Stern, C., Lee, S. M., et al. (2019). Optic disc edema after 30 days of strict head-down tilt bed rest. Ophthalmology 126, 467–468. doi: 10.1016/j.ophtha.2018.09.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Law, J., Van Baalen, M., Foy, M., Mason, S. S., Mendez, C., Wear, M. L., et al. (2014). Relationship between carbon dioxide levels and reported headaches on the international space station. J. Occup. Environ. Med. 56, 477–483. doi: 10.1097/JOM.0000000000000158

PubMed Abstract | CrossRef Full Text | Google Scholar

LeBlanc, A., Rowe, R., Schneider, V., Evans, H., and Hedrick, T. (1995). Regional muscle loss after short duration spaceflight. Aviat. Space Environ. Med. 66, 1151–1154.

PubMed Abstract | Google Scholar

Lee, A. G., Tarver, W. J., Mader, T. H., Gibson, C. R., Hart, S. F., and Otto, C. A. (2016). Neuro-ophthalmology of space flight. J. Neuroophthalmol. 36, 85–91. doi: 10.1097/WNO.0000000000000334

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, J. K., De Dios, Y. E., Kofman, I. S., Mulavara, A. P., Bloomberg, J., and Seidler, R. D. (2019a). Head down tilt bed rest plus elevated CO2 as a spaceflight analog: effects on cognitive and sensorimotor performance. Front. Hum. Neurosci. 13:355. doi: 10.3389/fnhum.2019.00355

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, J. K., Koppelmans, V., Riascos, R. F., Hasan, K. M., Pasternak, O., Mulavara, A. P., et al. (2019b). Spaceflight-associated brain white matter microstructural changes and intracranial fluid redistribution. JAMA Neurol. 76, 412–419. doi: 10.1001/jamaneurol.2018.4882

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, C.-C., Barker, J. W., Sparto, P. J., Furman, J. M., and Huppert, T. J. (2017). Functional near-infrared spectroscopy (fNIRS) brain imaging of multi-sensory integration during computerized dynamic posturography in middle-aged and older adults. Exp. Brain Res. 235, 1247–1256. doi: 10.1007/s00221-017-4893-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopez, C., Blanke, O., and Mast, F. (2012a). The human vestibular cortex revealed by coordinate-based activation likelihood estimation meta-analysis. Neuroscience 212, 159–179. doi: 10.1016/j.neuroscience.2012.03.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopez, C., Schreyer, H.-M., Preuss, N., and Mast, F. W. (2012b). Vestibular stimulation modifies the body schema. Neuropsychologia 50, 1830–1837. doi: 10.1016/j.neuropsychologia.2012.04.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Lowrey, C. R., Perry, S. D., Strzalkowski, N. D., Williams, D. R., Wood, S. J., and Bent, L. R. (2014). Selective skin sensitivity changes and sensory reweighting following short-duration space flight. J. Appl. Physiol. 116, 683–692. doi: 10.1152/japplphysiol.01200.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Manzey, D., and Lorenz, B. (1998). Effects of chronically elevated CO2 on mental performance during 26 days of confinement. Aviat. Space Environ. Med. 69, 506–514.

PubMed Abstract | Google Scholar

Miller, C. A., Kofman, I. S., Brady, R. R., May-Phillips, T. R., Batson, C. D., Lawrence, E. L., et al. (2018). Functional task and balance performance in bed rest subjects and astronauts. Aerosp. Med. Hum. Perform. 89, 805–815. doi: 10.3357/AMHP.5039.2018

PubMed Abstract | CrossRef Full Text | Google Scholar

Moore, S. T., Macdougall, H. G., and Paloski, W. H. (2010). Effects of head-down bed rest and artificial gravity on spatial orientation. Exp. Brain Res. 204, 617–622. doi: 10.1007/s00221-010-2317-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Mulavara, A., Peters, B., Miller, C., Kofman, I., Reschke, M., Taylor, L., et al. (2018). Physiological and functional alterations after spaceflight and bed rest. Med. Sci. Sports Exerc. 50, 1961–1980. doi: 10.1249/MSS.0000000000001615

PubMed Abstract | CrossRef Full Text | Google Scholar

Mulavara, A., Ruttley, T., Cohen, H., Peters, B., Miller, C., Brady, R., et al. (2012). Vestibular-somatosensory convergence in head movement control during locomotion after long-duration space flight. J. Vestib. Res. 22, 153–166. doi: 10.3233/ves-2011-0435

PubMed Abstract | CrossRef Full Text | Google Scholar

Mulavara, A. P., Feiveson, A. H., Fiedler, J., Cohen, H., Peters, B. T., Miller, C., et al. (2010). Locomotor function after long-duration space flight: effects and motor learning during recovery. Exp. Brain Res. 202, 649–659. doi: 10.1007/s00221-010-2171-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Mulder, E., Linnarsson, D., Paloski, W., Rittweger, J., Wuyts, F., Zange, J., et al. (2014). Effects of five days of bed rest with and without exercise countermeasure on postural stability and gait. J. Musculoskelet. Neuronal Interact. 14, 359–366.

PubMed Abstract | Google Scholar

Nichols, T. E., and Holmes, A. P. (2002). Nonparametric permutation tests for functional neuroimaging: a primer with examples. Hum. Brain Mapp. 15, 1–25. doi: 10.1002/hbm.1058

PubMed Abstract | CrossRef Full Text | Google Scholar

Noble, S., Scheinost, D., Finn, E. S., Shen, X., Papademetris, X., Mcewen, S. C., et al. (2017). Multisite reliability of MR-based functional connectivity. Neuroimage 146, 959–970. doi: 10.1016/j.neuroimage.2016.10.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Noohi, F., Kinnaird, C., De Dios, Y., Kofman, I., Wood, S. J., Bloomberg, J. J., et al. (2019). Deactivation of somatosensory and visual cortices during vestibular stimulation is associated with older age and poorer balance. PLoS One 14:e0221954. doi: 10.1371/journal.pone.0221954

PubMed Abstract | CrossRef Full Text | Google Scholar

Noohi, F., Kinnaird, C., Dedios, Y., Kofman, I. S., Wood, S., Bloomberg, J., et al. (2017). Functional brain activation in response to a clinical vestibular test correlates with balance. Front. Syst. Neurosci. 11:11. doi: 10.3389/fnsys.2017.00011

PubMed Abstract | CrossRef Full Text | Google Scholar

Parker, D., Reschke, M., Arrott, A., Homick, J., and Lichtenberg, B. (1985). Otolith tilt-translation reinterpretation following prolonged weightlessness: implications for preflight training. Aviat. Space Environ. Med. 56, 601–606.

PubMed Abstract | Google Scholar

Pechenkova, E., Nosikova, I., Rumshiskaya, A., Litvinova, L., Rukavishnikov, I., Mershina, E., et al. (2019). Alterations of functional brain connectivity after long-duration spaceflight as revealed by fMRI. Front. Physiol. 10:761. doi: 10.3389/fphys.2019.00761

PubMed Abstract | CrossRef Full Text | Google Scholar

R Core Team (2013). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.

Google Scholar

Reschke, M. F., Bloomberg, J. J., Paloski, W. H., Mulavara, A. P., Feiveson, A. H., and Harm, D. L. (2009). Postural reflexes, balance control, and functional mobility with long-duration head-down bed rest. Aviat. Space Environ. Med. 80, A45–A54.

PubMed Abstract | Google Scholar

Romero, J. E., Coupé, P., Giraud, R., Ta, V.-T., Fonov, V., Park, M. T. M., et al. (2017). CERES: a new cerebellum lobule segmentation method. Neuroimage 147, 916–924. doi: 10.1016/j.neuroimage.2016.11.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Rubin, A., Liedgren, S., Ödkvist, L., Larsby, B., and Aschan, G. (1979). Limb input to the cat vestibular nuclei. Acta Otolaryngol. 87, 113–122. doi: 10.3109/00016487909126395

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakellari, V., Bronstein, A., Corna, S., Hammon, C., Jones, S., and Wolsley, C. (1997). The effects of hyperventilation on postural control mechanisms. Brain J. Neurol. 120, 1659–1673. doi: 10.1093/brain/120.9.1659

PubMed Abstract | CrossRef Full Text | Google Scholar

Schlindwein, P., Mueller, M., Bauermann, T., Brandt, T., Stoeter, P., and Dieterich, M. (2008). Cortical representation of saccular vestibular stimulation: VEMPs in fMRI. Neuroimage 39, 19–31. doi: 10.1016/j.neuroimage.2007.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Sibonga, J. D. (2013). Spaceflight-induced bone loss: is there an osteoporosis risk? Curr. Osteoporosis Rep. 11, 92–98. doi: 10.1007/s11914-013-0136-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Spena, G., Gatignol, P., Capelle, L., and Duffau, H. (2006). Superior longitudinal fasciculus subserves vestibular network in humans. Neuroreport 17, 1403–1406. doi: 10.1097/01.wnr.0000223385.49919.61

PubMed Abstract | CrossRef Full Text | Google Scholar

Stein, T. (2013). Weight, muscle and bone loss during space flight: another perspective. Eur. J. Appl. Physiol. 113, 2171–2181. doi: 10.1007/s00421-012-2548-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Stenger, M. B., Tarver, W. J., Brunstetter, T., Gibson, C. R., Laurie, S. S., Lee, S., et al. (2017). Evidence Report: Risk of Spaceflight Associated Neuro-Ocular Syndrome (SANS). Available at: https://humanresearchroadmap.nasa.gov/evidence/reports/SANS.pdf?rnd=0.434276635495143 (accessed February 15, 2018).

Google Scholar

Wood, S. J., Paloski, W. H., and Clark, J. B. (2015). Assessing sensorimotor function following ISS with computerized dynamic posturography. Aerosp. Med. Hum. Perform. 86, A45–A53. doi: 10.3357/AMHP.EC07.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Yates, B., Jian, B., Cotter, L., and Cass, S. (2000). Responses of vestibular nucleus neurons to tilt following chronic bilateral removal of vestibular inputs. Exp. Brain Res. 130, 151–158. doi: 10.1007/s002219900238

PubMed Abstract | CrossRef Full Text | Google Scholar

Young, L. R., Oman, C. M., Watt, D., Money, K. E., and Lichtenberg, B. K. (1984). Spatial orientation in weightlessness and readaptation to earth’s gravity. Science 225, 205–208. doi: 10.1126/science.6610215

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, P., Koppelmans, V., Reuter-Lorenz, P., De Dios, Y., Gadd, N., Riascos, R., et al. (2018a). Change of cortical foot activation following 70 days of head-down bed rest. J. Neurophysiol. 119, 2145–2152. doi: 10.1152/jn.00693.2017

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, P., Koppelmans, V., Reuter-Lorenz, P., De Dios, Y., Gadd, N., Wood, S., et al. (2018b). Vestibular brain changes within 70 days of head down bed rest. Hum. Brain Mapp. 39, 2753–2763. doi: 10.1002/hbm.24037

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, P., Koppelmans, V., Reuter-Lorenz, P. A., De Dios, Y. E., Gadd, N. E., Wood, S. J., et al. (2016). Increased brain activation for dual tasking with 70-days head-down bed rest. Front. Syst. Neurosci. 10:71. doi: 10.3389/fnsys.2016.00071

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, H., Saidel, G. M., and LaManna, J. C. (2008). Cerebral blood flow adaptation to chronic hypoxia. Adv Exp Med Biol. 614, 371–377. doi: 10.1007/978-0-387-74911-2_41

PubMed Abstract | CrossRef Full Text | Google Scholar

Zu Eulenburg, P., Caspers, S., Roski, C., and Eickhoff, S. B. (2012). Meta-analytical definition and functional connectivity of the human vestibular cortex. Neuroimage 60, 162–169. doi: 10.1016/j.neuroimage.2011.12.032

PubMed Abstract | CrossRef Full Text | Google Scholar