- 1Department of Radiology, Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, United States
- 2Department of Neurology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States
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.
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).
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. Cross-correlation of ΔBOLD with PETCO2. (A) Brain maps of cross-correlation of ΔBOLD with PETCO2 for individual HD subjects. The strength of cross-correlation increases from light color to dark color. Brain regions without color indicate that the ΔBOLD in these areas are not correlated with the changes in PETCO2. Multiple comparisons are corrected using clustering at an overall p < 0.05. (B) Group comparison of Fisher-transformed cross-correlation values between HD subjects and controls after adjusting for age. Cold colors represent the weaker cross-correlation found in HD subjects than in controls. Multiple comparisons were corrected at false discovery rate pfdr < 0.05.
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. Association between Fisher-transformed cross-correlation values and disease designation (HD vs. Control) after adjusting for age.
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.
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.
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.
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.
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.
Footnotes
References
Banzett, R. B., Garcia, R. T., and Moosavi, S. H. (2000). Simple contrivance “clamps” end-tidal PCO(2) and PO(2) despite rapid changes in ventilation. J. Appl. Physiol. 88, 1597–1600. doi: 10.1152/jappl.2000.88.5.1597
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practival and powerful approach to multiple testing. J. R. Statist. Soc. 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
Brown, R., Benveniste, H., Black, S. E., Charpak, S., Dichgans, M., Joutel, A., et al. (2018). Understanding the role of the perivascular space in cerebral small vessel disease. Cardiovasc. Res. 114, 1462–1473. doi: 10.1093/cvr/cvy113
Chan, S. T., Evans, K. C., Song, T. Y., Selb, J., van der Kouwe, A., Rosen, B. R., et al. (2020a). Dynamic brain-body coupling of breath-by-breath O2-CO2 exchange ratio with resting state cerebral hemodynamic fluctuations. PLoS One 15:e0238946. doi: 10.1371/journal.pone.0238946
Chan, S. T., Mercaldo, N. D., Ravina, B., Hersch, S. M., and Rosas, H. D. (2020b). Association of dilated perivascular spaces and disease severity in patients with Huntington’s disease. Neurology 2020:11121. doi: 10.1212/WNL.0000000000011121
Chen, J. J., Salat, D. H., and Rosas, H. D. (2012). Complex relationships between cerebral blood flow and brain atrophy in early Huntington’s disease. Neuroimage 59, 1043–1051. doi: 10.1016/j.neuroimage.2011.08.112
Conneally, P. M. (1984). Huntington disease: genetics and epidemiology. Am. J. Hum. Genet. 36, 506–526.
Cox, R. W. A. F. N. I. (1996). software for analysis and visualization of functional magnetic resonance neuroimages. Comput. Biomed. Res. 29, 162–173. doi: 10.1006/cbmr.1996.0014
Dale, A. M., Fischl, B., and Sereno, M. I. (1999). Cortical surface-based analysis. I. Segmentation and surface reconstruction. Neuroimage 9, 179–194. doi: 10.1006/nimg.1998.0395
Deckel, A. W., and Duffy, J. D. (2000). Vasomotor hyporeactivity in the anterior cerebral artery during motor activation in Huntington’s disease patients. Brain Res. 872, 258–261. doi: 10.1016/s0006-8993(00)02506-3
Drouin-Ouellet, J., Sawiak, S. J., Cisbani, G., Lagace, M., Kuan, W. L., Saint-Pierre, M., et al. (2015). Cerebrovascular and blood-brain barrier impairments in Huntington’s disease: Potential implications for its pathophysiology. Ann. Neurol. 78, 160–177. doi: 10.1002/ana.24406
Fierstra, J., Sobczyk, O., Battisti-Charbonney, A., Mandell, D. M., Poublanc, J., Crawley, A. P., et al. (2013). Measuring cerebrovascular reactivity: what stimulus to use? J. Physiol. 591, 5809–5821. doi: 10.1113/jphysiol.2013.259150
Fischl, B., Sereno, M. I., and Dale, A. M. (1999). Cortical surface-based analysis. II: Inflation, flattening, and a surface-based coordinate system. Neuroimage 9, 195–207. doi: 10.1006/nimg.1998.0396
Fisher, R. A. (1921). On the ‘probable error’ of a coefficient of correlation deduced from a small sample. Metron 1, 1–32.
Franciosi, S., Ryu, J. K., Shim, Y., Hill, A., Connolly, C., Hayden, M. R., et al. (2012). Age-dependent neurovascular abnormalities and altered microglial morphology in the YAC128 mouse model of Huntington disease. Neurobiol. Dis. 45, 438–449. doi: 10.1016/j.nbd.2011.09.003
Gold, S., Christian, B., Arndt, S., Zeien, G., Cizadlo, T., Johnson, D. L., et al. (1998). Functional MRI statistical software packages: a comparative analysis. Hum. Brain Mapp. 6, 73–84. doi: 10.1002/(sici)1097-0193(1998)6:2<73::aid-hbm1>3.0.co;2-h
Hersch, S. M., and Rosas, H. D. (2008). Neuroprotection for Huntington’s disease: ready, set, slow. Neurotherapeutics 5, 226–236. doi: 10.1016/j.nurt.2008.01.003
Hobbs, N. Z., Barnes, J., Frost, C., Henley, S. M., Wild, E. J., Macdonald, K., et al. (2010). Onset and progression of pathologic atrophy in Huntington disease: a longitudinal MR imaging study. AJNR Am. J. Neuroradiol. 31, 1036–1041. doi: 10.3174/ajnr.A2018
Hobbs, N. Z., Pedrick, A. V., Say, M. J., Frost, C., Dar Santos, R., Coleman, A., et al. (2011). The structural involvement of the cingulate cortex in premanifest and early Huntington’s disease. Mov. Disord. 26, 1684–1690. doi: 10.1002/mds.23747
Hsiao, H. Y., Chen, Y. C., Huang, C. H., Chen, C. C., Hsu, Y. H., Chen, H. M., et al. (2015). Aberrant astrocytes impair vascular reactivity in Huntington disease. Ann. Neurol. 78, 178–192. doi: 10.1002/ana.24428
Hua, J., Unschuld, P. G., Margolis, R. L., van Zijl, P. C., and Ross, C. A. (2014). Elevated arteriolar cerebral blood volume in prodromal Huntington’s disease. Mov. Disord. 29, 396–401. doi: 10.1002/mds.25591
Jansen, A. H., van Hal, M., Op den Kelder, I. C., Meier, R. T., de Ruiter, A. A., Schut, M. H., et al. (2017). Frequency of nuclear mutant huntingtin inclusion formation in neurons and glia is cell-type-specific. Glia 65, 50–61. doi: 10.1002/glia.23050
Kisler, K., Nelson, A. R., Montagne, A., and Zlokovic, B. V. (2017). Cerebral blood flow regulation and neurovascular dysfunction in Alzheimer disease. Nat. Rev. Neurosci. 18, 419–434. doi: 10.1038/nrn.2017.48
Lim, R. G., Quan, C., Reyes-Ortiz, A. M., Lutz, S. E., Kedaigle, A. J., Gipson, T. A., et al. (2017). Huntington’s Disease iPSC-Derived Brain Microvascular Endothelial Cells Reveal WNT-Mediated Angiogenic and Blood-Brain Barrier Deficits. Cell Rep. 19, 1365–1377. doi: 10.1016/j.celrep.2017.04.021
Lin, C. Y., Hsu, Y. H., Lin, M. H., Yang, T. H., Chen, H. M., Chen, Y. C., et al. (2013). Neurovascular abnormalities in humans and mice with Huntington’s disease. Exp. Neurol. 250, 20–30. doi: 10.1016/j.expneurol.2013.08.019
Liu, P., De Vis, J. B., and Lu, H. (2019). Cerebrovascular reactivity (CVR) MRI with CO2 challenge: A technical review. Neuroimage 187, 104–115. doi: 10.1016/j.neuroimage.2018.03.047
Liu, Q., Li, C., Wanga, V., and Shepherd, B. E. (2018). Covariate-adjusted Spearman’s rank correlation with probability-scale residuals. Biometrics 74, 595–605. doi: 10.1111/biom.12812
Marstrand, J. R., Garde, E., Rostrup, E., Ring, P., Rosenbaum, S., Mortensen, E. L., et al. (2002). Cerebral perfusion and cerebrovascular reactivity are reduced in white matter hyperintensities. Stroke 33, 972–976. doi: 10.1161/01.str.0000012808.81667.4b
McKay, L. C., Evans, K. C., Frackowiak, R. S., and Corfield, D. R. (2003). Neural correlates of voluntary breathing in humans. J. Appl. Physiol. 95, 1170–1178. doi: 10.1152/japplphysiol.00641.2002
Nation, D. A., Sweeney, M. D., Montagne, A., Sagare, A. P., D’Orazio, L. M., Pachicano, M., et al. (2019). Blood-brain barrier breakdown is an early biomarker of human cognitive dysfunction. Nat. Med. 25, 270–276. doi: 10.1038/s41591-018-0297-y
MacDonald, M. E., Ambrose, C. M., Duyao, M. P., Myers, R. H., Lin, C., Srinidhi, L., 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. Cell 72, 971–983. doi: 10.1016/0092-8674(93)90585-e
Padel, T., Roth, M., Gaceb, A., Li, J. Y., Bjorkqvist, M., and Paul, G. (2018). Brain pericyte activation occurs early in Huntington’s disease. Exp. Neurol. 305, 139–150. doi: 10.1016/j.expneurol.2018.03.015
Rosas, H. D., Liu, A. K., Hersch, S., Glessner, M., Ferrante, R. J., Salat, D. H., et al. (2002). Regional and progressive thinning of the cortical ribbon in Huntington’s disease. Neurology 58, 695–701. doi: 10.1212/wnl.58.5.695
Rosas, H. D., Salat, D. H., Lee, S. Y., Zaleta, A. K., Pappu, V., Fischl, B., et al. (2008). Cerebral cortex and the clinical expression of Huntington’s disease: complexity and heterogeneity. Brain 131(Pt 4), 1057–1068. doi: 10.1093/brain/awn025
Rosas, H. D., Wilkens, P., Salat, D. H., Mercaldo, N. D., Vangel, M., Yendiki, A. Y., et al. (2018). Complex spatial and temporally defined myelin and axonal degeneration in Huntington disease. Neuroimage Clin. 20, 236–242. doi: 10.1016/j.nicl.2018.01.029
Saad, Z. S., DeYoe, E. A., and Ropella, K. M. (2003). Estimation of FMRI response delays. Neuroimage 18, 494–504.
Sam, K., Crawley, A. P., Conklin, J., Poublanc, J., Sobczyk, O., Mandell, D. M., et al. (2016). Development of White Matter Hyperintensity Is Preceded by Reduced Cerebrovascular Reactivity. Ann. Neurol. 80, 277–285. doi: 10.1002/ana.24712
Schmithorst, V. J., Hernandez-Garcia, L., Vannest, J., Rajagopal, A., Lee, G., and Holland, S. K. (2014). Optimized simultaneous ASL and BOLD functional imaging of the whole brain. J. Magnetic Resonance Imaging JMRI 39, 1104–1117. doi: 10.1002/jmri.24273
Stapor, P. C., Sweat, R. S., Dashti, D. C., Betancourt, A. M., and Murfee, W. L. (2014). Pericyte dynamics during angiogenesis: new insights from new identities. J. Vasc. Res. 51, 163–174. doi: 10.1159/000362276
Sweeney, M. D., Sagare, A. P., and Zlokovic, B. V. (2018). Blood-brain barrier breakdown in Alzheimer disease and other neurodegenerative disorders. Nat. Rev. Neurol. 14, 133–150. doi: 10.1038/nrneurol.2017.188
Talairach, J., and Tournoux, P. (1988). Co-Planar Stereotaxic Atlas of the Human Brain. New York, NY: Thieme Medical.
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.
Edited by:
James Duffin, University of Toronto, CanadaReviewed by:
Lashmi Venkatraghavan, University of Toronto, CanadaJorn Fierstra, University Hospital Zurich, Switzerland
Copyright © 2021 Chan, Mercaldo, Kwong, Hersch and Rosas. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Herminia D. Rosas, rosas@helix.mgh.harvard.edu
†These authors have contributed equally to this work and share first authorship