ORIGINAL RESEARCH article

Front. Physiol., 21 July 2021

Sec. Medical Physics and Imaging

Volume 12 - 2021 | https://doi.org/10.3389/fphys.2021.663898

Impaired Cerebrovascular Reactivity in Huntington’s Disease

  • 1. Department of Radiology, Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, United States

  • 2. Department of Neurology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States

Abstract

There is increasing evidence that impairments of cerebrovascular function and/or abnormalities of the cerebral vasculature might contribute to early neuronal cell loss in Huntington’s disease (HD). Studies in both healthy individuals as well as in patients with other neurodegenerative disorders have used an exogenous carbon dioxide (CO2) challenge in conjunction with functional magnetic resonance imaging (fMRI) to assess regional cerebrovascular reactivity (CVR). In this study, we explored potential impairments of CVR in HD. Twelve gene expanded HD individuals, including both pre-symptomatic and early symptomatic HD and eleven healthy controls were administered a gas mixture targeting a 4–8 mmHg increase in CO2 relative to the end-tidal partial pressure of CO2 (PETCO2) at rest. A Hilbert Transform analysis was used to compute the cross-correlation between the time series of regional BOLD signal changes (ΔBOLD) and increased PETCO2, and to estimate the response delay of ΔBOLD relative to PETCO2. After correcting for age, we found that the cross-correlation between the time series for regional ΔBOLD and for PETCO2 was weaker in HD subjects than in controls in several subcortical white matter regions, including the corpus callosum, subcortical white matter adjacent to rostral and caudal anterior cingulate, rostral and caudal middle frontal, insular, middle temporal, and posterior cingulate areas. In addition, greater volume of dilated perivascular space (PVS) was observed to overlap, primarily along the periphery, with the areas that showed greater ΔBOLD response delay. Our preliminary findings support that alterations in cerebrovascular function occur in HD and may be an important, not as yet considered, contributor to early neuropathology in HD.

Introduction

Huntington’s Disease (HD) is a devastating fully penetrant autosomal dominant progressive, rare neurological disorder that is characterized by progressive motor dysfunction, emotional disturbances, dementia, and weight loss (Conneally, 1984; Hersch and Rosas, 2008). It is caused by the abnormal expansion of a CAG repeat length in the gene that codes for an ubiquitously expressed protein, huntingtin (mHtt) (MacDonald et al., 1993), which can be found in neurons as well as all major components of the neurovascular unit including the basal lamina, endothelial cell, pericytes, smooth muscle cells (Lin et al., 2013; Drouin-Ouellet et al., 2015), astrocytes, and oligodendrocytes and microglia (Jansen et al., 2017). Chronic mHtt expression has been shown to alter the neurovasculature in an HD transgenic mouse model by upregulating wingless-related integration site (WNT) signaling activity and consequently inhibiting proper endothelial cell differentiation and maturation in the brain (Lim et al., 2017). The upregulation of WNT signaling activity in HD endothelial cells appeared to induce an increase in pericyte number (Stapor et al., 2014) and occurred early in the disease, leading to blood vessel sprouting, elongation, and early maturation. By late stage of disease, prominent aberrant vasculature was present (Padel et al., 2018). Alterations in the neurovasculature could result in increased cerebral blood volume, reduced cerebral blood flow (CBF), increased small vessel density and increased blood-brain barrier (BBB) permeability, as has reported in rodent models of HD and patient post-mortem tissue (Chen et al., 2012; Franciosi et al., 2012; Lin et al., 2013; Hua et al., 2014; Drouin-Ouellet et al., 2015; Hsiao et al., 2015). These studies provided early support for neurovascular dysfunction in HD. Impairments in neurovascular function have also been associated with enlarged perivascular spaces (PVS), which are considered not only a hallmark of small vessel disease (Brown et al., 2018), but of perivascular inflammation and impaired clearance of proteins (Brown et al., 2018). Using a novel automated segmentation algorithm (Chan et al., 2020b), we found that the load of enlarged PVS was greater in patients with HD in a manner that directly correlated with disease severity and caudate atrophy. Together, these findings converge to provide support for early impairment of cerebrovascular function, including prior to neuronal cell loss or to early clinical symptoms including weight loss, behavioral changes (Padel et al., 2018) or motor dysfunction (Lin et al., 2013; Drouin-Ouellet et al., 2015) in HD.

To date, studies have not as yet evaluated early cerebrovascular dysfunction in vivo in HD. In particular, mild physiological stress can be used to depict subtle vascular abnormalities prior to the failure in autoregulation. It is possible that early cerebrovascular dysfunction could contribute to early neuropathological changes in HD. In this study, we used hypercapnic stress, administering an air mixture with a low concentration of carbon dioxide (CO2) as a mild physiological stress, to probe regional cerebrovascular reactivity (CVR) in HD gene-expanded individuals, both motor pre-symptomatic and early symptomatic individuals. We mapped regional cerebrovascular responses to the CO2 challenge using blood oxygenation level dependent (BOLD) signal changes measured with functional magnetic resonance imaging (fMRI). This approach has a higher temporal resolution and higher contrast to noise ratio than arterial spin labeling (ASL), and might therefore be better able to more precisely demonstrate the vascular responses in white matter (Schmithorst et al., 2014). The extent of regional impairment in CVR was evaluated using the cross-correlation between regional BOLD signal changes and increased CO2 level, indicated by end-tidal partial pressure of CO2 (PETCO2), using a Hilbert Transform analysis (Saad et al., 2003). A correlation was also employed with the total functional capacity (TFC), a measure of clinical severity, to explore the potential clinical impact associated with impaired cerebrovascular function. Our findings support that early changes in CVR occur in HD and support the importance of considering neurovascular alterations as a key pathophysiological mechanism that could both contribute to subclinical progression and could also serve as a novel target for intervention.

Materials and Methods

Participants

Huntington’s disease gene-expanded individuals, both motor pre-symptomatic and early symptomatic, were recruited from the Center of Excellence in the Department of Neurology at Massachusetts General Hospital. Healthy volunteers were recruited through the Partners hospital network and were screened to exclude neurological, mental and medical disorders, drug abuse, and contraindications for exposure to high magnetic field. Clinical assessments were conducted by an HD specialist (HDR). All MRI scanning was performed at the Athinoula A. Martinos Center for Biomedical Imaging at the Massachusetts General Hospital. All the experimental procedures were explained to the subjects. Signed informed consent was obtained before any study procedures were conducted. All components of this study were performed in compliance with the Declaration of Helsinki and all procedures were approved by the MGH Human Research Committee (IRB Protocol Number: 2014P002240).

MRI Acquisition

Magnetic resonance imaging brain scanning was performed on a 3-Tesla scanner (Siemens Medical, Erlangen, Germany). The head was immobilized in a standard head coil with foam pads. The following whole brain MRI datasets were acquired: (1) standard high-resolution 3D sagittal images acquired with volumetric T1-weighted 3D-MEMPRAGE sequence (TR = 2530 ms, TE = 1.74, 3.6, 5.46, and 7.32 ms, flip angle = 7°, FOV = 256 × 256 mm, matrix = 256 × 256, slice thickness = 1 mm); (2) T2-SPACE (TR = 3200 ms, TE = 454, flip angle = 120°, FOV = 256 × 256 mm, matrix = 256 × 256, slice thickness = 1 mm); and (3) BOLD-fMRI images acquired with gradient-echo echo planar imaging (EPI) sequence (TR = 1250 ms, TE = 30 ms, flip angle = 90°, FOV = 220 × 220 mm, matrix = 64 × 64, thickness = 5 mm, slice gap = 1 mm) during the hypercapnic challenge.

Hypercapnic Challenge With Externally Administered CO2

Subjects wore a nose-clip and breathed through a mouth-piece using an MRI-compatible circuit designed to maintain the PETCO2 within ±1–2 mmHg of the target PETCO2 (Banzett et al., 2000; McKay et al., 2003). Given the potential for inter-individual variance in the resting end-tidal partial pressure of CO2 (PETCO2) (West, 1992), the resting PETCO2 was assessed for each individual subject via calibrated capnography before the hypercapnic challenge. The fraction of inspired carbon dioxide was adjusted to produce steady-state conditions of normocapnia and mild hypercapnia (4–8 mmHg above the subject’s resting PETCO2). The CO2 challenge paradigm consisted of 2 consecutive phases (normocapnia and mild hypercapnia) repeating six times with three epochs of 4 mmHg increase and three epochs of 8 mmHg increase of PETCO2 (Figure 1). The normocapnia phase lasted 60–90 s, while the mild hypercapnia phase lasted 30 s. The total duration of the hypercapnic challenge lasted 10 min.

FIGURE 1

Physiological Recording

Physiological changes, including the partial pressure of carbon dioxide (PCO2) and oxygen (PO2), and respiratory flow, were measured simultaneously with the MRI acquisition. PCO2 and PO2 were sampled through an air filter connected to the mouthpiece and the sampled gases were measured using gas analyzers (Capstar-100, Oxystar-100, CWE, Inc., PA, United States), after calibrating to the barometric pressure of the day of MRI scanning and correcting for vapor pressure. Respiratory flow was measured using a respiratory flow head (MTL300L, ADInstruments, Inc., CO, United States) on a breathing circuit via calibrated spirometer (FE141, ADInstruments, Inc., CO, United States). The recording of PCO2, PO2, respiratory flow was integrated into Powerlab and LabChart (ADInstruments, Inc., CO, United States). The physiological measurements were synchronized with the MRI images using trigger signals from the scanner. The BOLD-fMRI images and physiological recordings were stored for off-line data analysis.

Data Analysis

Processing of Physiological Data

Physiological data were analyzed using Matlab R2014a (Mathworks, Inc., Natick, MA, United States). Technical delays of PCO2 and PO2 were corrected by cross-correlating the time series of PCO2 and PO2 with the respiratory flow tracing acquired with the spirometer. End inspiration (I) and end expiration (E) were defined on the time series of PO2 and PCO2 and were verified by the inspiratory and expiratory phases of the respiratory flow time series. The breath-to-breath PETCO2 and end-tidal O2 (PETO2) were extracted at the end expiration of PCO2 and PO2 time series, respectively.

Cross-Correlation Analysis Between BOLD-fMRI Data and PETCO2

The BOLD-fMRI data were imported into the Analysis of Functional NeuroImage (AFNI) (Cox, 1996) (National Institute of Mental Health)1 software for data analysis. Details of preprocessing with AFNI are given as follows. The first 12 volumes of each functional dataset that were collected before equilibrium magnetization was reached were discarded. Each functional dataset was corrected for slice timing, motion-corrected and co-registered to the first image of the first functional dataset using a three-dimensional volume registration. In the co-registered functional dataset, the time series of each voxel was de-trended with the fifth order of polynomials to remove the low drift frequency, and the components of motion were removed using orthogonal projection. The clean functional dataset was then normalized to its mean intensity value across the time series. Voxels located within the ventricles and outside the brain defined in the parcellated brain volume using FreeSurfer version 5.3 (Dale et al., 1999; Fischl et al., 1999) (MGH/MIT/HMS Athinoula A. Martinos Center for Biomedical Imaging, Boston)2 were excluded. Individual subject brain volumes with time series of percent BOLD signal changes (ΔBOLD) were derived. Taking advantage of Fast Fourier transforms, a Hilbert Transform analysis (Saad et al., 2003) was used to measure the cross-correlation between the ΔBOLD and the reference time series of PETCO2 in each voxel and the associated delay of BOLD response relative to the changes in PETCO2.

The statistical parametric maps for individual subjects were cluster-corrected using a threshold estimated with a Monte Carlo simulation algorithm. Individual subject brain volumes with cross-correlation coefficients were registered onto each subject’s anatomical scan and transformed to the standardized space of Talairach and Tournoux (Talairach and Tournoux, 1988). In order to protect against type I error, an individual voxel probability threshold of p < 0.005 was held to correct the overall significance level to α < 0.05. Monte Carlo simulation was used to correct for multiple comparisons (Gold et al., 1998). Based upon a Monte Carlo simulation, using 2000 iteration processed with ClustSim program (Ward, 1997), it was estimated that a 226 mm3 contiguous volume would provide the overall corrected threshold of p < 0.05.

Group Analysis of Demographics, CVR, and Disease Severity Between Patients and Controls

Descriptive summaries were computed by diagnostic group. Categorical variables were summarized as frequencies and percentages and group differences were assessed using either the Chi-square test or Fisher’s exact test. Continuous variables were summarized as median (interquartile range = 25th and 75th percentiles) and group differences were assessed using Kruskal-Wallis test.

Cross-correlation coefficients in individual subject maps were transformed to z-scores for the group analysis (z = ln[(1 + r)/(1-r)]/2) (Fisher, 1921). For each subject, these z-scores were averaged in each of the 160 brain regions parcellated using FreeSurfer version 5.3 (Dale et al., 1999; Fischl et al., 1999). A linear regression model was constructed to quantify the association between CVR, as indicated by the average z-score, and disease designation (HD and control) after adjusting for age. Separate models were estimated for each region. For each model, linear combinations of parameter estimates were computed to estimate average z-values by group (fixing age = 50), and their associated 95% confidence intervals. Differences in average z-scores (HD-control) were also estimated, after adjusting for age, and the 95% confidence and unadjusted p-value.

Partial spearman correlation coefficients (Liu et al., 2018) were estimated to summarize the monotonic relationship between z-scores and TFC for all subjects after adjusting for age. False discovery rate adjusted p-values were computed to account for multiple comparisons (Benjamini and Hochberg, 1995). All group analyses were performed using R 4.0.23.

Results

Twenty-three participants (14 males, 9 females, aged from 22 to 62 years) were included. Eleven were healthy volunteers (median = 32.0 years, interquartile range IQR = 27.5–46.0 years; 9M and 2F) and the remaining 12 were HD gene-expanded individuals (median = 51.5 years, interquartile range IQR = 47.3–57.3 years; 5M and 7F), including both pre-symptomatic and early symptomatic HD subjects. The selected clinical characteristics of the gene-expanded subjects are shown in Table 1. Significant differences were observed between HD and control groups in terms of age (p = 0.026) but marginally for gender (p = 0.123).

TABLE 1

ControlHDp
N1112
Age32.00 [27.50, 46.00]51.50 [47.25, 57.25]0.026*
Sex: Male9 (81.8)5 (41.7)0.123
TFC13 [11, 13]

Subject demographics.

Categorical variable (i.e., sex) is summarized as frequency (percentage) and continuous variable (i.e., age) is summarized as median [interquartile range]. TFC, total functional capacity. *significant difference p < 0.05.

Group Comparison of Cross-Correlation of ΔBOLD With PETCO2 Between HD Gene-Expanded Subjects and Controls

Brain maps with the cross-correlation of ΔBOLD with PETCO2 superimposed on the T1 image and the associated ΔBOLD response delay superimposed on the T2-T1 ratio image for individual gene-expanded subjects are shown in Figure 2 and Supplementary Figure 1, respectively.

FIGURE 2

After correcting for age, we found that the cross-correlation between the time series for regional ΔBOLD and for PETCO2 in several subcortical white matter regions was significantly weaker in gene-expanded subjects than in controls (Figure 2 and Table 2). Regions of significance included corpus callosum (mid-posterior segment: HD-Control estimate = −0.23, 95% CI = −0.37, −0.08, p = 0.002, pfdr = 0.041; central segment: HD-Control estimate = −0.24, 95% CI = −0.38, −0.11, p = 0.001, pfdr < 0.029), subcortical white matter adjacent to rostral and caudal anterior cingulate, rostral and caudal middle frontal, insular, middle temporal, and posterior cingulate areas. Greater dilated PVS load was also observed to overlap with and along the periphery of the areas that showed greater ΔBOLD response delays (Supplementary Figure 1). In the symptomatic HD group, no significant correlation was found between TFC and the z-scores indicating the cross correlation of ΔBOLD with PETCO2 after adjusting for age (pfdr > 0.05).

TABLE 2

Corpus callosumControlHDHD-ControlUnadjusted
Estimate (95% CI)Estimate (95% CI)Estimate (95% CI)p-value
Posterior segment0.44 (0.34,0.55)0.40 (0.33,0.48)−0.04 (−0.17, 0.09)0.560
Midposterior segment0.52 (0.41,0.64)0.30 (0.22,0.38)−0.23 (−0.37, −0.08)0.002*
Central segment0.53 (0.42,0.64)0.28 (0.20,0.36)−0.24 (−0.38, −0.11)0.001*
Midanterior segment0.53 (0.38,0.67)0.32 (0.21,0.42)−0.21 (−0.39, −0.03)0.020
Anterior segment0.56 (0.41,0.71)0.32 (0.22,0.43)−0.24 (−0.42, −0.06)0.010

Left HemisphereRight Hemisphere

Brain regionsControl  Estimate (95% CI)HD Estimate (95% CI)HD-Control Estimate (95% CI)Unadjusted p-valueControl  Estimate (95% CI)HD Estimate (95% CI)HD-Control Estimate (95% CI)Unadjusted p-value

Thalamus0.59 (0.42, 0.75)0.57 (0.45, 0.69)−0.02 (−0.23, 0.19)0.8580.58 (0.41, 0.74)0.56 (0.45, 0.68)−0.01 (−0.22, 0.19)0.904
Caudate0.55 (0.38, 0.71)0.51 (0.39, 0.62)−0.04 (−0.24, 0.16)0.6960.56 (0.39, 0.72)0.52 (0.40, 0.64)−0.04 (−0.24, 0.16)0.706
Putamen0.54 (0.36, 0.71)0.52 (0.39, 0.65)−0.02 (−0.24, 0.20)0.8660.57 (0.40, 0.74)0.52 (0.40, 0.64)−0.05 (−0.26, 0.16)0.631
Pallidum0.55 (0.43, 0.66)0.37 (0.29, 0.46)−0.18 (−0.32, −0.03)0.0170.54 (0.42, 0.66)0.37 (0.28, 0.45)−0.18 (−0.33, −0.03)0.022
Hippocampus0.63 (0.49, 0.76)0.47 (0.37, 0.56)−0.16 (−0.33, 0.01)0.0570.64 (0.50, 0.77)0.45 (0.35, 0.54)−0.19 (−0.36, −0.02)0.027
Amygdala0.59 (0.45, 0.73)0.41 (0.30, 0.51)−0.18 (−0.36, −0.01)0.0420.55 (0.42, 0.69)0.38 (0.29, 0.48)−0.17 (−0.34, 0.00)0.045
Accumbens area0.45 (0.30, 0.59)0.39 (0.28, 0.50)−0.05 (−0.24, 0.13)0.5540.46 (0.33, 0.58)0.35 (0.25, 0.44)−0.11 (−0.27, 0.05)0.164
Banks of superior temporal sulcus0.68 (0.50, 0.86)0.57 (0.44, 0.70)−0.11 (−0.33, 0.12)0.3430.69 (0.51, 0.87)0.57 (0.44, 0.70)−0.12 (−0.34, 0.10)0.285
Caudal anterior cingulate0.63 (0.47, 0.79)0.63 (0.52, 0.75)0.00 (−0.19, 0.20)0.9690.63 (0.47, 0.79)0.61 (0.50, 0.73)−0.02 (−0.22, 0.18)0.856
Caudal middle frontal0.71 (0.54, 0.87)0.58 (0.46, 0.70)−0.13 (−0.33, 0.08)0.2200.67 (0.50, 0.85)0.60 (0.47, 0.73)−0.07 (−0.29, 0.15)0.513
Cuneus0.65 (0.45, 0.86)0.61 (0.46, 0.76)−0.04 (−0.30, 0.21)0.7410.61 (0.41, 0.81)0.60 (0.45, 0.75)−0.01 (−0.26, 0.24)0.918
Entorhinal0.39 (0.30, 0.48)0.31 (0.24, 0.37)−0.09 (−0.20, 0.03)0.1380.35 (0.25, 0.46)0.32 (0.25, 0.40)−0.03 (−0.16, 0.10)0.657
Fusiform0.55 (0.39, 0.71)0.56 (0.44, 0.68)0.01 (−0.19, 0.21)0.9050.54 (0.39, 0.69)0.54 (0.43, 0.65)0.00 (−0.18, 0.19)0.985
Inferior parietal0.71 (0.55, 0.87)0.55 (0.43, 0.67)−0.16 (−0.36, 0.04)0.1270.70 (0.53, 0.86)0.55 (0.43, 0.67)−0.15 (−0.36, 0.06)0.163
Inferior temporal0.45 (0.34, 0.56)0.43 (0.35, 0.51)−0.02 (−0.16, 0.11)0.7480.45 (0.33, 0.56)0.43 (0.35, 0.52)−0.02 (−0.16, 0.13)0.834
Isthmus of cingulate0.55 (0.37, 0.72)0.58 (0.46, 0.71)0.04 (−0.18, 0.25)0.7340.56 (0.39, 0.73)0.59 (0.46, 0.71)0.03 (−0.18, 0.24)0.809
Lateral occipital0.63 (0.47, 0.80)0.55 (0.44, 0.67)−0.08 (−0.28, 0.12)0.4360.62 (0.46, 0.78)0.58 (0.46, 0.69)−0.05 (−0.24, 0.15)0.648
Lateral orbitofrontal0.45 (0.33, 0.57)0.38 (0.30, 0.47)−0.07 (−0.21, 0.08)0.3550.43 (0.32, 0.54)0.41 (0.33, 0.49)−0.02 (−0.16, 0.12)0.783
Lingual0.57 (0.38, 0.76)0.62 (0.48, 0.75)0.04 (−0.19, 0.28)0.7250.57 (0.38, 0.76)0.60 (0.46, 0.74)0.04 (−0.20, 0.27)0.759
Medial orbitofrontal0.46 (0.36, 0.56)0.34 (0.27, 0.42)−0.12 (−0.24, 0.01)0.0600.45 (0.35, 0.54)0.35 (0.29, 0.42)−0.09 (−0.21, 0.02)0.118
Middle temporal0.57 (0.43, 0.70)0.48 (0.38, 0.57)−0.09 (−0.25, 0.08)0.2870.58 (0.45, 0.70)0.48 (0.39, 0.58)−0.09 (−0.25, 0.06)0.237
Parahippocampal0.55 (0.42, 0.68)0.47 (0.37, 0.56)−0.09 (−0.25, 0.08)0.3060.54 (0.40, 0.69)0.44 (0.33, 0.54)−0.11 (−0.29, 0.07)0.248
Paracentral0.64 (0.45, 0.83)0.66 (0.52, 0.80)0.02 (−0.22, 0.26)0.8670.63 (0.45, 0.81)0.67 (0.53, 0.80)0.03 (−0.19, 0.26)0.764
Pars opercularis0.65 (0.47, 0.84)0.61 (0.47, 0.74)−0.05 (−0.27, 0.18)0.6760.63 (0.46, 0.80)0.62 (0.50, 0.75)−0.01 (−0.22, 0.20)0.923
Pars orbitalis0.51 (0.35, 0.67)0.49 (0.37, 0.60)−0.02 (−0.22, 0.18)0.8380.49 (0.34, 0.65)0.54 (0.42, 0.65)0.04 (−0.15, 0.23)0.677
Pars triangularis0.59 (0.42, 0.76)0.54 (0.42, 0.67)−0.04 (−0.25, 0.17)0.6840.58 (0.43, 0.74)0.56 (0.45, 0.68)−0.02 (−0.22, 0.18)0.843
Peri-calcarine0.62 (0.40, 0.84)0.62 (0.46, 0.79)0.00 (−0.28, 0.28)0.9840.60 (0.39, 0.81)0.60 (0.44, 0.75)0.00 (−0.26, 0.26)0.983
Postcentral0.67 (0.50, 0.84)0.62 (0.50, 0.75)−0.05 (−0.26, 0.17)0.6720.64 (0.47, 0.81)0.62 (0.49, 0.74)−0.03 (−0.23, 0.18)0.800
Posterior cingulate0.64 (0.46, 0.82)0.61 (0.48, 0.74)−0.03 (−0.25, 0.20)0.8150.61 (0.44, 0.78)0.60 (0.48, 0.73)−0.01 (−0.22, 0.21)0.960
Precentral0.68 (0.51, 0.84)0.59 (0.47, 0.71)−0.08 (−0.29, 0.12)0.4260.65 (0.48, 0.81)0.60 (0.48, 0.72)−0.04 (−0.25, 0.16)0.661
Precuneus0.63 (0.45, 0.81)0.59 (0.46, 0.71)−0.04 (−0.26, 0.17)0.6950.64 (0.46, 0.82)0.58 (0.45, 0.71)−0.06 (−0.28, 0.16)0.602
Rostral anterior cingulate0.57 (0.41, 0.73)0.52 (0.41, 0.64)−0.04 (−0.24, 0.16)0.6790.58 (0.44, 0.72)0.50 (0.40, 0.60)−0.08 (−0.25, 0.09)0.357
Rostral middle frontal0.65 (0.50, 0.81)0.54 (0.43, 0.66)−0.11 (−0.30, 0.08)0.2660.64 (0.48, 0.79)0.56 (0.45, 0.67)−0.08 (−0.26, 0.11)0.419
Superior frontal0.63 (0.47, 0.78)0.58 (0.46, 0.69)−0.05 (−0.25, 0.14)0.6020.63 (0.47, 0.80)0.58 (0.46, 0.70)−0.05 (−0.25, 0.15)0.602
Superior parietal0.68 (0.52, 0.84)0.60 (0.48, 0.71)−0.08 (−0.28, 0.11)0.4060.68 (0.52, 0.84)0.59 (0.47, 0.71)−0.09 (−0.29, 0.11)0.361
Superior temporal0.63 (0.48, 0.78)0.52 (0.40, 0.63)−0.11 (−0.30, 0.07)0.2380.62 (0.47, 0.78)0.54 (0.42, 0.65)−0.09 (−0.28, 0.11)0.381
Supramarginal0.69 (0.52, 0.86)0.62 (0.50, 0.74)−0.07 (−0.28, 0.14)0.5130.65 (0.48, 0.82)0.62 (0.50, 0.74)−0.03 (−0.24, 0.18)0.765
Frontal pole0.36 (0.25, 0.48)0.28 (0.19, 0.36)−0.09 (−0.23, 0.06)0.2430.41 (0.28, 0.54)0.33 (0.23, 0.42)−0.08 (−0.25, 0.08)0.303
Transverse temporal0.63 (0.42, 0.84)0.61 (0.46, 0.77)−0.02 (−0.27, 0.24)0.9080.62 (0.41, 0.82)0.64 (0.49, 0.78)0.02 (−0.23, 0.27)0.873
Insula0.63 (0.45, 0.82)0.54 (0.41, 0.68)−0.09 (−0.32, 0.14)0.4410.62 (0.45, 0.79)0.56 (0.43, 0.68)−0.06 (−0.27, 0.15)0.570
Banks of superior temporal sulcus WM0.61 (0.46, 0.76)0.43 (0.32, 0.54)−0.18 (−0.36, 0.01)0.0590.57 (0.44, 0.71)0.40 (0.31, 0.50)−0.17 (−0.33, −0.01)0.042
Caudal anterior cingulate WM0.65 (0.52, 0.78)0.40 (0.31, 0.50)−0.25 (−0.40, −0.09)0.002*0.66 (0.55, 0.78)0.40 (0.31, 0.48)−0.27 (−0.41, −0.12)< 0.001*
Caudal middle frontal WM0.63 (0.51, 0.75)0.40 (0.31, 0.49)−0.23 (−0.38, −0.07)0.004*0.64 (0.49, 0.78)0.44 (0.33, 0.54)−0.20 (−0.38, −0.02)0.026
Cuneus WM0.68 (0.50, 0.86)0.59 (0.46, 0.72)−0.10 (−0.32, 0.13)0.3970.67 (0.49, 0.84)0.57 (0.44, 0.70)−0.10 (−0.31, 0.12)0.389
Entorhinal WM0.47 (0.36, 0.57)0.31 (0.23, 0.39)−0.16 (−0.29, −0.03)0.0170.49 (0.38, 0.61)0.31 (0.23, 0.40)−0.18 (−0.33, −0.04)0.014
Fusiform WM0.59 (0.45, 0.73)0.46 (0.35, 0.56)−0.13 (−0.30, 0.04)0.1370.57 (0.44, 0.69)0.45 (0.35, 0.54)−0.12 (−0.28, 0.03)0.125
Inferior parietal WM0.64 (0.49, 0.78)0.44 (0.34, 0.54)−0.20 (−0.37, −0.02)0.0290.61 (0.48, 0.75)0.42 (0.33, 0.52)−0.19 (−0.36, −0.02)0.026
Inferior temporal WM0.52 (0.42, 0.63)0.40 (0.32, 0.48)−0.12 (−0.25, 0.01)0.0650.50 (0.40, 0.60)0.40 (0.32, 0.47)−0.10 (−0.23, 0.03)0.119
Isthmus of cingulate WM0.61 (0.49, 0.73)0.44 (0.35, 0.53)−0.17 (−0.32, −0.02)0.0250.60 (0.49, 0.72)0.43 (0.34, 0.51)−0.17 (−0.32, −0.03)0.017
Lateral occipital WM0.63 (0.48, 0.78)0.49 (0.38, 0.60)−0.14 (−0.32, 0.05)0.1390.61 (0.47, 0.74)0.49 (0.39, 0.58)−0.12 (−0.28, 0.05)0.155
Lateral orbitofrontal WM0.53 (0.42, 0.65)0.35 (0.27, 0.43)−0.18 (−0.32, −0.05)0.0090.50 (0.40, 0.60)0.38 (0.31, 0.46)−0.12 (−0.24, 0.01)0.066
Lingual WM0.61 (0.45, 0.78)0.54 (0.42, 0.66)−0.07 (−0.27, 0.13)0.4990.63 (0.47, 0.79)0.52 (0.41, 0.64)−0.11 (−0.30, 0.09)0.284
Medial orbitofrontal WM0.52 (0.42, 0.62)0.35 (0.28, 0.42)−0.17 (−0.30, −0.05)0.0060.51 (0.41, 0.61)0.35 (0.27, 0.42)−0.16 (−0.29, −0.04)0.012
Middle temporal WM0.63 (0.49, 0.77)0.44 (0.34, 0.54)−0.19 (−0.36, −0.02)0.0330.63 (0.52, 0.74)0.42 (0.34, 0.50)−0.21 (−0.35, −0.07)0.003*
Parahippocampal WM0.61 (0.48, 0.74)0.44 (0.34, 0.53)−0.17 (−0.34, −0.01)0.0330.61 (0.49, 0.73)0.42 (0.33, 0.50)−0.20 (−0.35, −0.05)0.010
Paracentral WM0.59 (0.44, 0.74)0.48 (0.36, 0.59)−0.12 (−0.31, 0.07)0.2310.53 (0.39, 0.67)0.46 (0.36, 0.56)−0.07 (−0.24, 0.10)0.415
Pars opercularis WM0.70 (0.56, 0.84)0.47 (0.37, 0.58)−0.22 (−0.40, −0.05)0.0120.67 (0.53, 0.82)0.46 (0.35, 0.56)−0.22 (−0.39, −0.04)0.015
Pars orbitalis WM0.59 (0.43, 0.75)0.47 (0.35, 0.58)−0.12 (−0.31, 0.08)0.2360.57 (0.42, 0.73)0.51 (0.40, 0.62)−0.06 (−0.25, 0.13)0.533
Pars triangularis WM0.63 (0.49, 0.78)0.43 (0.33, 0.54)−0.20 (−0.38, −0.03)0.0240.63 (0.50, 0.76)0.45 (0.36, 0.55)−0.18 (−0.34, −0.01)0.033
Peri-calcarine WM0.65 (0.48, 0.81)0.54 (0.42, 0.66)−0.10 (−0.31, 0.10)0.3270.63 (0.47, 0.79)0.52 (0.41, 0.64)−0.11 (−0.30, 0.09)0.292
Postcentral WM0.71 (0.55, 0.86)0.55 (0.44, 0.67)−0.16 (−0.35, 0.04)0.1120.67 (0.52, 0.82)0.55 (0.44, 0.66)−0.12 (−0.31, 0.07)0.209
Posterior cingulate WM0.65 (0.52, 0.77)0.39 (0.30, 0.48)−0.26 (−0.41, −0.10)0.001*0.60 (0.48, 0.71)0.39 (0.31, 0.47)−0.20 (−0.34, −0.07)0.004*
Precentral WM0.63 (0.50, 0.76)0.43 (0.34, 0.53)−0.20 (−0.36, −0.03)0.0190.60 (0.47, 0.73)0.43 (0.34, 0.53)−0.17 (−0.33, 0.00)0.048
Precuneus WM0.62 (0.49, 0.75)0.46 (0.36, 0.56)−0.16 (−0.32, 0.00)0.0560.59 (0.46, 0.73)0.44 (0.35, 0.54)−0.15 (−0.32, 0.02)0.088
Rostral anterior cingulate WM0.62 (0.50, 0.74)0.39 (0.30, 0.48)−0.23 (−0.39, −0.08)0.003*0.63 (0.51, 0.74)0.36 (0.28, 0.45)−0.26 (−0.40, −0.12)<0.001*
Rostral middle frontal WM0.66 (0.53, 0.79)0.42 (0.32, 0.51)−0.24 (−0.40, −0.08)0.003*0.66 (0.54, 0.78)0.43 (0.35, 0.52)−0.23 (−0.37, −0.08)0.003*
Superior frontal WM0.65 (0.51, 0.78)0.45 (0.35, 0.55)−0.20 (−0.37, −0.03)0.0240.63 (0.50, 0.76)0.44 (0.35, 0.54)−0.19 (−0.35, −0.02)0.027
Superior parietal WM0.63 (0.50, 0.76)0.45 (0.36, 0.55)−0.18 (−0.34, −0.02)0.0270.58 (0.45, 0.72)0.46 (0.36, 0.56)−0.13 (−0.30, 0.05)0.149
Superior temporal WM0.65 (0.52, 0.79)0.43 (0.33, 0.52)−0.23 (−0.39, −0.07)0.0060.68 (0.54, 0.82)0.48 (0.37, 0.58)−0.21 (−0.38, −0.03)0.021
Supramarginal WM0.68 (0.54, 0.82)0.48 (0.38, 0.58)−0.20 (−0.37, −0.03)0.0210.61 (0.48, 0.74)0.46 (0.37, 0.56)−0.15 (−0.31, 0.02)0.084
Frontal pole WM0.37 (0.25, 0.49)0.31 (0.22, 0.40)−0.06 (−0.21, 0.09)0.4530.45 (0.29, 0.60)0.40 (0.28, 0.51)−0.05 (−0.24, 0.14)0.600
Transverse temporal WM0.69 (0.52, 0.86)0.52 (0.40, 0.65)−0.17 (−0.38, 0.04)0.1190.66 (0.49, 0.82)0.56 (0.44, 0.68)−0.10 (−0.31, 0.11)0.356
Insula WM0.67 (0.55, 0.80)0.44 (0.35, 0.53)−0.23 (−0.39, −0.08)0.004*0.68 (0.55, 0.80)0.45 (0.36, 0.54)−0.23 (−0.38, −0.07)0.004*
Deep white matter0.51 (0.40, 0.61)0.32 (0.24, 0.40)−0.19 (−0.32, −0.06)0.0050.49 (0.38, 0.60)0.33 (0.25, 0.41)−0.16 (−0.30, −0.03)0.016

Association between Fisher-transformed cross-correlation values and disease designation (HD vs. Control) after adjusting for age.

Disease designation-specific estimates, and 95% confidence intervals, of Fisher-transformed cross-correlation among 50-years are provided, along with difference between HD subjects and controls (regardless of age, due to the lack of interaction between age and designation), and the unadjusted p-vales (p). * significant difference pfdr < 0.05. WM, white matter.

Discussion

Evidence from transgenic mouse models has suggested that cerebrovascular dysfunction might contribute to the early neuropathology in HD. To our knowledge, our study is the first to confirm regional impairments of cerebrovascular function in vivo using an externally administered CO2 challenge in human subjects with pre-symptomatic and early HD.

In our study, impaired CVR was predominantly found in subcortical white matter regions in HD gene-expanded individuals (Figure 2), including in the body of the corpus callosum, forceps major, superior longitudinal fasciculus, and the cingulum bundle. These correspond to white matter regions that are both known to demonstrate altered microstructural integrity (Rosas et al., 2018) and adjacent to cortical brain areas shown previously with reduced cerebral blood flow (Chen et al., 2012) and/or cortical atrophy (Rosas et al., 2002, 2008; Hobbs et al., 2010, 2011). Increased PVS burden was also found in several subcortical white matter regions including the insula, caudal middle frontal and middle temporal areas (Chan et al., 2020b) in HD. Impaired cerebrovascular function and increased PVS burden were associated with each other in cortical areas also showing reduced cerebral blood flow. These findings, together, suggest that cerebrovascular dysfunction is both early and significant in HD (Brown et al., 2018). Furthermore, they suggest that alterations in cerebrovascular function could potentially independently contribute to the neuropathology of HD.

This was supported by the presence of mutant huntingtin (mHtt) (MacDonald et al., 1993) not only in neurons but in all major components of the neurovascular unit including the basal lamina, endothelial cell, pericytes, smooth muscle cells (Lin et al., 2013; Drouin-Ouellet et al., 2015), astrocytes, oligodendrocytes, and microglia (Jansen et al., 2017). Chronic mHtt expression appears to not only alter the neurovasculature by the upregulation of WNT signaling (Lim et al., 2017) but also by affecting the expression and activity of endothelial nitric oxide (Deckel and Duffy, 2000), a vasoactive agent which could affect CVR. Disrupted cerebrovascular function may lead to reductions in CBF and prolonged relative hypoxia and/or neurovascular coupling (Kisler et al., 2017), which could then result in an increased PVS burden (Brown et al., 2018). This is consistent with our observation in HD gene-expanded subjects that greater dilated PVS load appeared to overlap with and occur along the periphery of the areas showing greater delays in cerebrovascular response (Supplementary Figure 1). Increases in PVS load have been associated with neuroinflammation, pericyte loss and breakdown of the blood-brain barrier (Sweeney et al., 2018; Nation et al., 2019), further exacerbating pathological cascades. Delayed cerebrovascular responses (Marstrand et al., 2002) or negative BOLD responses (Sam et al., 2016) to hypercapnic challenge have been previously reported in brain areas with white matter hyperintensities, compared to normal appearing white matter by as much as more than a year (Sam et al., 2016). Our finding of these changes including during the motor pre-manifest period supports that these alterations could influence not only the neuropathology of HD but contribute to clinical symptoms, including early cognitive deficits (Brown et al., 2018).

Taken together, our results support an important role for altered CVR in the pathophysiology of HD and support future longitudinal studies using CVR challenge to evaluate disease onset and/or disease progression.

Limitations of This Study

Since this is the first study to evaluate CVR with externally administered CO2 to HD patients, we limited our study of gene-expanded subjects who were either motor pre-manifest or in early stages of HD. Therefore, our study is limited by a relatively small sample size and with the exclusion of more advanced patients who might have also demonstrated more profound alterations. Another limitation is the age of the HD subjects and healthy controls was significantly different; this could influence the interpretation of the differences in CVR between these two subject groups. In order to address this, we included age as a regressor of no interest to correct for the potential contribution of aging in the group analyses. Increased arterial blood pressure resulting from moderate hypercapnia can be a confound in the CVR assessment. Although we limited the mild hypercapnic condition to 4–8 mmHg above the subject’s resting PETCO2, it would be ideal to continuously measure and monitor the blood pressure changes while the subjects are under hypercapnic challenge which we did not incorporate in the current study. Variations stemming from gas delivery apparatus, design of the breathing circuit, or of the protocol, to the analysis of the imaging data in CVR MRI studies using externally administered CO2 challenge have been reported previously (Fierstra et al., 2013; Liu et al., 2019; Chan et al., 2020a). These may be due in part to the level of the complexity of the set-up of the CO2 challenge, varying tolerability in different subject or patient populations, variability in the physiological responses, in addition to differences in analytic approaches used to measure the response. It would be ideal to have a standardized vasoactive stimulus available that could be applied in different settings and in different subject/patient populations; this should be an important consideration in future studies.

Conclusion

Our preliminary findings support that alterations in neurovascular function occur in HD and may be an important, not as yet considered, independent contributor to early neuropathology and clinical symptoms in HD.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Ethics statement

The studies involving human participants were reviewed and approved by MGH Human Research Committee. The patients/participants provided their written informed consent to participate in this study.

Author contributions

SC and HR conceived and designed the study. SC and HR performed the experiments. SC and NM analyzed the data. SC, NM, KK, SH, and HR interpreted results of the experiments. SC, NM, and HR edited and revised the manuscript. All authors reviewed and approved the final version of manuscript.

Funding

This work was supported by Dake Family Foundation, and the National Institutes of Health R01NS042861, R01NS114562, and R21AT010955. It was also supported by the resources provided by the Shared Instrumentation Grant S10RR023043.

Acknowledgments

We are very grateful to the individuals who participated in this study, who so generously contributed their time and energy to this work and without whom it would not have been possible. We thank the individuals who helped with the recruitment and scanning of patients, including Lois Kwon, Natalie Connors, Paul Wilkens, Matt Linnehan, and Eugene Hsu.

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2021.663898/full#supplementary-material

Supplementary Figure 1

Maps of delayed cerebrovascular responses to hypercapnic challenge superimposed on T2-T1 ratio maps for individual subjects. Warm colors represent long delay response while cold colors represent short delay. Multiple comparisons are corrected using clustering at an overall p < 0.05.

References

  • 1

    BanzettR. B.GarciaR. T.MoosaviS. H. (2000). Simple contrivance “clamps” end-tidal PCO(2) and PO(2) despite rapid changes in ventilation.J. Appl. Physiol.8815971600. 10.1152/jappl.2000.88.5.1597

  • 2

    BenjaminiY.HochbergY. (1995). Controlling the false discovery rate: a practival and powerful approach to multiple testing.J. R. Statist. Soc.57289300. 10.1111/j.2517-6161.1995.tb02031.x

  • 3

    BrownR.BenvenisteH.BlackS. E.CharpakS.DichgansM.JoutelA.et al (2018). Understanding the role of the perivascular space in cerebral small vessel disease.Cardiovasc. Res.11414621473. 10.1093/cvr/cvy113

  • 4

    ChanS. T.EvansK. C.SongT. Y.SelbJ.van der KouweA.RosenB. R.et al (2020a). Dynamic brain-body coupling of breath-by-breath O2-CO2 exchange ratio with resting state cerebral hemodynamic fluctuations.PLoS One15:e0238946. 10.1371/journal.pone.0238946

  • 5

    ChanS. T.MercaldoN. D.RavinaB.HerschS. M.RosasH. D. (2020b). Association of dilated perivascular spaces and disease severity in patients with Huntington’s disease.Neurology2020:11121. 10.1212/WNL.0000000000011121

  • 6

    ChenJ. J.SalatD. H.RosasH. D. (2012). Complex relationships between cerebral blood flow and brain atrophy in early Huntington’s disease.Neuroimage5910431051. 10.1016/j.neuroimage.2011.08.112

  • 7

    ConneallyP. M. (1984). Huntington disease: genetics and epidemiology.Am. J. Hum. Genet.36506526.

  • 8

    CoxR. W. A. F. N. I. (1996). software for analysis and visualization of functional magnetic resonance neuroimages.Comput. Biomed. Res.29162173. 10.1006/cbmr.1996.0014

  • 9

    DaleA. M.FischlB.SerenoM. I. (1999). Cortical surface-based analysis. I. Segmentation and surface reconstruction.Neuroimage9179194. 10.1006/nimg.1998.0395

  • 10

    DeckelA. W.DuffyJ. D. (2000). Vasomotor hyporeactivity in the anterior cerebral artery during motor activation in Huntington’s disease patients.Brain Res.872258261. 10.1016/s0006-8993(00)02506-3

  • 11

    Drouin-OuelletJ.SawiakS. J.CisbaniG.LagaceM.KuanW. L.Saint-PierreM.et al (2015). Cerebrovascular and blood-brain barrier impairments in Huntington’s disease: Potential implications for its pathophysiology.Ann. Neurol.78160177. 10.1002/ana.24406

  • 12

    FierstraJ.SobczykO.Battisti-CharbonneyA.MandellD. M.PoublancJ.CrawleyA. P.et al (2013). Measuring cerebrovascular reactivity: what stimulus to use?J. Physiol.59158095821. 10.1113/jphysiol.2013.259150

  • 13

    FischlB.SerenoM. I.DaleA. M. (1999). Cortical surface-based analysis. II: Inflation, flattening, and a surface-based coordinate system.Neuroimage9195207. 10.1006/nimg.1998.0396

  • 14

    FisherR. A. (1921). On the ‘probable error’ of a coefficient of correlation deduced from a small sample.Metron1132.

  • 15

    FranciosiS.RyuJ. K.ShimY.HillA.ConnollyC.HaydenM. R.et al (2012). Age-dependent neurovascular abnormalities and altered microglial morphology in the YAC128 mouse model of Huntington disease.Neurobiol. Dis.45438449. 10.1016/j.nbd.2011.09.003

  • 16

    GoldS.ChristianB.ArndtS.ZeienG.CizadloT.JohnsonD. L.et al (1998). Functional MRI statistical software packages: a comparative analysis.Hum. Brain Mapp.67384. 10.1002/(sici)1097-0193(1998)6:2<73::aid-hbm1>3.0.co;2-h

  • 17

    HerschS. M.RosasH. D. (2008). Neuroprotection for Huntington’s disease: ready, set, slow.Neurotherapeutics5226236. 10.1016/j.nurt.2008.01.003

  • 18

    HobbsN. Z.BarnesJ.FrostC.HenleyS. M.WildE. J.MacdonaldK.et al (2010). Onset and progression of pathologic atrophy in Huntington disease: a longitudinal MR imaging study.AJNR Am. J. Neuroradiol.3110361041. 10.3174/ajnr.A2018

  • 19

    HobbsN. Z.PedrickA. V.SayM. J.FrostC.Dar SantosR.ColemanA.et al (2011). The structural involvement of the cingulate cortex in premanifest and early Huntington’s disease.Mov. Disord.2616841690. 10.1002/mds.23747

  • 20

    HsiaoH. Y.ChenY. C.HuangC. H.ChenC. C.HsuY. H.ChenH. M.et al (2015). Aberrant astrocytes impair vascular reactivity in Huntington disease.Ann. Neurol.78178192. 10.1002/ana.24428

  • 21

    HuaJ.UnschuldP. G.MargolisR. L.van ZijlP. C.RossC. A. (2014). Elevated arteriolar cerebral blood volume in prodromal Huntington’s disease.Mov. Disord.29396401. 10.1002/mds.25591

  • 22

    JansenA. H.van HalM.Op den KelderI. C.MeierR. T.de RuiterA. A.SchutM. H.et al (2017). Frequency of nuclear mutant huntingtin inclusion formation in neurons and glia is cell-type-specific.Glia655061. 10.1002/glia.23050

  • 23

    KislerK.NelsonA. R.MontagneA.ZlokovicB. V. (2017). Cerebral blood flow regulation and neurovascular dysfunction in Alzheimer disease.Nat. Rev. Neurosci.18419434. 10.1038/nrn.2017.48

  • 24

    LimR. G.QuanC.Reyes-OrtizA. M.LutzS. E.KedaigleA. J.GipsonT. A.et al (2017). Huntington’s Disease iPSC-Derived Brain Microvascular Endothelial Cells Reveal WNT-Mediated Angiogenic and Blood-Brain Barrier Deficits.Cell Rep.1913651377. 10.1016/j.celrep.2017.04.021

  • 25

    LinC. Y.HsuY. H.LinM. H.YangT. H.ChenH. M.ChenY. C.et al (2013). Neurovascular abnormalities in humans and mice with Huntington’s disease.Exp. Neurol.2502030. 10.1016/j.expneurol.2013.08.019

  • 26

    LiuP.De VisJ. B.LuH. (2019). Cerebrovascular reactivity (CVR) MRI with CO2 challenge: A technical review.Neuroimage187104115. 10.1016/j.neuroimage.2018.03.047

  • 27

    LiuQ.LiC.WangaV.ShepherdB. E. (2018). Covariate-adjusted Spearman’s rank correlation with probability-scale residuals.Biometrics74595605. 10.1111/biom.12812

  • 28

    MarstrandJ. R.GardeE.RostrupE.RingP.RosenbaumS.MortensenE. L.et al (2002). Cerebral perfusion and cerebrovascular reactivity are reduced in white matter hyperintensities.Stroke33972976. 10.1161/01.str.0000012808.81667.4b

  • 29

    McKayL. C.EvansK. C.FrackowiakR. S.CorfieldD. R. (2003). Neural correlates of voluntary breathing in humans.J. Appl. Physiol.9511701178. 10.1152/japplphysiol.00641.2002

  • 30

    NationD. A.SweeneyM. D.MontagneA.SagareA. P.D’OrazioL. M.PachicanoM.et al (2019). Blood-brain barrier breakdown is an early biomarker of human cognitive dysfunction.Nat. Med.25270276. 10.1038/s41591-018-0297-y

  • 31

    MacDonaldM. E.AmbroseC. M.DuyaoM. P.MyersR. H.LinC.SrinidhiL.et al (1993). A novel gene containing a trinucleotide repeat that is expanded and unstable on Huntington’s disease chromosomes. The Huntington’s Disease Collaborative Research Group.Cell72971983. 10.1016/0092-8674(93)90585-e

  • 32

    PadelT.RothM.GacebA.LiJ. Y.BjorkqvistM.PaulG. (2018). Brain pericyte activation occurs early in Huntington’s disease.Exp. Neurol.305139150. 10.1016/j.expneurol.2018.03.015

  • 33

    RosasH. D.LiuA. K.HerschS.GlessnerM.FerranteR. J.SalatD. H.et al (2002). Regional and progressive thinning of the cortical ribbon in Huntington’s disease.Neurology58695701. 10.1212/wnl.58.5.695

  • 34

    RosasH. D.SalatD. H.LeeS. Y.ZaletaA. K.PappuV.FischlB.et al (2008). Cerebral cortex and the clinical expression of Huntington’s disease: complexity and heterogeneity.Brain131(Pt 4), 10571068. 10.1093/brain/awn025

  • 35

    RosasH. D.WilkensP.SalatD. H.MercaldoN. D.VangelM.YendikiA. Y.et al (2018). Complex spatial and temporally defined myelin and axonal degeneration in Huntington disease.Neuroimage Clin.20236242. 10.1016/j.nicl.2018.01.029

  • 36

    SaadZ. S.DeYoeE. A.RopellaK. M. (2003). Estimation of FMRI response delays.Neuroimage18494504.

  • 37

    SamK.CrawleyA. P.ConklinJ.PoublancJ.SobczykO.MandellD. M.et al (2016). Development of White Matter Hyperintensity Is Preceded by Reduced Cerebrovascular Reactivity.Ann. Neurol.80277285. 10.1002/ana.24712

  • 38

    SchmithorstV. J.Hernandez-GarciaL.VannestJ.RajagopalA.LeeG.HollandS. K. (2014). Optimized simultaneous ASL and BOLD functional imaging of the whole brain.J. Magnetic Resonance Imaging JMRI3911041117. 10.1002/jmri.24273

  • 39

    StaporP. C.SweatR. S.DashtiD. C.BetancourtA. M.MurfeeW. L. (2014). Pericyte dynamics during angiogenesis: new insights from new identities.J. Vasc. Res.51163174. 10.1159/000362276

  • 40

    SweeneyM. D.SagareA. P.ZlokovicB. V. (2018). Blood-brain barrier breakdown in Alzheimer disease and other neurodegenerative disorders.Nat. Rev. Neurol.14133150. 10.1038/nrneurol.2017.188

  • 41

    TalairachJ.TournouxP. (1988). Co-Planar Stereotaxic Atlas of the Human Brain.New York, NY: Thieme Medical.

  • 42

    WardL. L. (1997). Simultaneous inference for fMRI data.Wisconsin: Medical College of Wisconsin.

  • 43

    WestJ. B. (1992). Pulmonary Pathophysiology.Baltimore: Williams & Wilkins.

Summary

Keywords

cerebrovascular reactivity, exogenous CO2 challenge, Huntington’s disease, functional magnetic resonance imaging, presymptomatic, perivascular space

Citation

Chan ST, Mercaldo ND, Kwong KK, Hersch SM and Rosas HD (2021) Impaired Cerebrovascular Reactivity in Huntington’s Disease. Front. Physiol. 12:663898. doi: 10.3389/fphys.2021.663898

Received

15 February 2021

Accepted

09 June 2021

Published

21 July 2021

Volume

12 - 2021

Edited by

James Duffin, University of Toronto, Canada

Reviewed by

Lashmi Venkatraghavan, University of Toronto, Canada; Jorn Fierstra, University Hospital Zurich, Switzerland

Updates

Copyright

*Correspondence: Herminia D. Rosas,

These authors have contributed equally to this work and share first authorship

This article was submitted to Medical Physics and Imaging, a section of the journal Frontiers in Physiology

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics