Skip to main content

ORIGINAL RESEARCH article

Front. Aging Neurosci., 12 November 2021
Sec. Neurocognitive Aging and Behavior
Volume 13 - 2021 | https://doi.org/10.3389/fnagi.2021.666495

Magnetic Resonance Imaging-Visible Perivascular Spaces in the Basal Ganglia Are Associated With the Diabetic Retinopathy Stage and Cognitive Decline in Patients With Type 2 Diabetes

Eun Young Choi1 Yae Won Park2* Minyoung Lee3 Min Kim1 Christopher Seungkyu Lee1 Sung Soo Ahn2 Jinna Kim2 Seung-Koo Lee2
  • 1Department of Ophthalmology, Yonsei University College of Medicine, Seoul, South Korea
  • 2Department of Radiology, Center for Clinical Imaging Data Science, Research Institute of Radiological Science, Yonsei University College of Medicine, Seoul, South Korea
  • 3Department of Endocrinology, Yonsei University College of Medicine, Seoul, South Korea

Purpose: The aim of this study was to evaluate whether perivascular space (PVS) severity and retinal ganglion cell layer (GCL) thickness differed based on the stage of diabetic retinopathy (DR) and the cognitive status in patients with DR.

Methods: A total of 81 patients with DR (51 in the non-proliferative group and 30 in the proliferative group) were included in this retrospective, cross-sectional study. PVS severity was assessed in the basal ganglia (BG) and centrum semiovale using MRI. The total cerebral small vessel disease (SVD) score was determined based on the numbers of lacunes and microbleeds and the severity of white matter hyperintensity. Optical coherence tomography was used to measure foveal and perifoveal GCL thicknesses. Cerebral SVD markers and cognitive function were compared between the groups, and correlations between the BG-PVS severity and the Mini-Mental Status Examination (MMSE) scores and GCL parameters were evaluated.

Results: Patients with proliferative DR had higher BG-PVS severity (P = 0.012), higher total cerebral SVD scores (P = 0.035), reduced GCL thicknesses in the inferior (P = 0.027), superior (P = 0.046), and temporal (P = 0.038) subfields compared to patients with non-proliferative DR. In addition, the BG-PVS severity was negatively correlated with the MMSE score (P = 0.007), and the GCL thickness was negatively correlated with the BG-PVS severity (P-values < 0.05 for inferior, superior, and temporal subfields).

Conclusion: BG-PVS severity and retinal GCL thickness may represent novel imaging biomarkers reflecting the stage of DR and cognitive decline in diabetic patients. Furthermore, these results suggest a possible link between cerebral and retinal neurodegeneration at the clinical level.

Introduction

Diabetic retinopathy (DR) is a major microvascular complication of type 2 diabetes mellitus (T2DM) and a predictor of other end-organ complications (Aminian et al., 2020), with a global prevalence of 27% in patients with T2DM (Thomas et al., 2019). DR can be divided into two major categories, namely, non-proliferative DR in the early stages and proliferative DR in the more advanced stages. While pathological preretinal neovascularization is a hallmark feature of proliferative DR (Stitt et al., 2016), retinal neurodegeneration is reported to occur earlier in DR pathogenesis (Sohn et al., 2016). Retinal ganglion and amacrine cells are believed to be affected first by diabetes-induced apoptosis (Ng et al., 2016), which can be identified clinically by a decrease in the ganglion cell layer (GCL) thickness on optical coherence tomography (OCT) imaging (Ng et al., 2016). Moreover, previous studies have shown that patients with DR have an increased risk of cognitive impairment (Crosby-Nwaobi et al., 2012; Gupta et al., 2019), although the pathophysiological mechanisms underlying this process are not well known.

While the retina is the classic target organ for diabetic microangiopathy, the brain has also recently been identified as another target organ for diabetic microvascular complications (Woerdeman et al., 2014). The retina and the brain share similar anatomical and physiological characteristics, as both are embryologically derived from the same tissue and possess similar structural and functional features, such as the blood-retinal/blood-brain barriers and glial cell connections (London et al., 2013). In addition, similar to the retinal microvasculature, an increasing body of evidence indicates that hyperglycemia can induce dysfunction in the cerebral microvasculature, causing pathological neurovascular remodeling, blood-brain barrier disruption, and cell damage to neurons and glial cells (Ergul et al., 2012). The retinal glymphatic system, a network between the superficial capillary plexus and the inner nuclear layers, has also been reported to reflect the state of neurodegenerative and inflammatory disorders in the brain (Petzold, 2016).

Cerebral small vessel disease (SVD) encompasses a wide spectrum of cerebrovascular diseases with similar clinical manifestations. Using MRI, conventional markers for SVD include lacunes, microbleeds, and white matter hyperintensity (WMH) (Chen et al., 2019). Cerebral perivascular spaces (PVS), also known as Virchow-Robin spaces, are pial-lined, interstitial, fluid-filled spaces surrounding the penetrating vessels. A growing body of evidence suggests that enlarged PVS represent a novel biomarker for cerebral SVD (Duperron et al., 2018). Enlarged PVS have been correlated with cerebrovascular disease (Hurford et al., 2014; Potter et al., 2015) and neurodegenerative diseases such as dementia (Ramirez et al., 2016; Banerjee et al., 2017; Sepehrband et al., 2021). Moreover, PVS have been shown to be a component of the glymphatic system, a recently discovered, macroscopic, waste-clearance system in the brain (Jessen et al., 2015). A recent study has demonstrated the glymphatic system damage in T2DM rat models (Jiang et al., 2017), and the enlarged PVS may reflect this damage (Jessen et al., 2015). However, the relationships between cerebral PVS enlargement, retinal neurodegeneration, DR stage, and cognitive decline have not yet been evaluated. Therefore, the comprehensive assessments of the brain and retinal changes should be performed to elucidate the underlying mechanisms for these processes at the clinical level and to acquire practical insights into the design of appropriate treatments for the prevention of cognitive impairment in patients with DR.

Therefore, in this study, we investigated whether MRI-visible PVS, other cerebral SVD markers, and retinal GCL thickness differed based on the DR stage and the cognitive status of patients with DR.

Materials and Methods

Patient Population

This retrospective study was approved by the Severance Institutional Review Board, and a waiver of informed consent was obtained (IRB approval number: 2020-3812-001). In total, 158 patients with DR who had undergone detailed neurological examinations for the assessment of cognitive impairment from February 2010 to April 2020 were identified. Among them, 130 patients with DR who had undergone both brain MRI and neuropsychological testing within 1 year of ophthalmologic examinations were enrolled. In this prospective study, participants were excluded based on the following criteria: (1) presence of other retinal diseases (e.g., age-related macular degeneration, central serous chorioretinopathy, retinal vascular occlusion, epiretinal membrane, and macular hole) or optic nerve diseases (e.g., optic neuritis, glaucoma, and ischemic optic neuropathy) (n = 19); (2) presence of old infarcts due to large vessel disease or cardioembolic strokes, past intracerebral hemorrhages, or post-traumatic encephalomalacia detected on MRI (n = 17); (3) a previous history of intravitreal or sub-Tenon’s injections, vitrectomy surgery, or laser photocoagulation for the treatment of diabetic macular edema (n = 7); or (4) a current or past history of severe dementia, other neurological or neurodevelopmental disorders, substance-related problems, or head trauma with loss of consciousness (n = 6). After the consideration of these exclusion criteria, a total of 81 consecutive patients with DR were included in the final study population (mean age = 67.3 ± 9.9 years; 49 females and 32 males).

Demographic and clinical characteristics of patients, such as age, sex, vascular risk factors (e.g., hypertension, dyslipidemia, or cardiovascular diseases), and current or remote history of smoking, were collected. The Mini-Mental Status Examination (MMSE) (Han et al., 2008) assessed a range of elements that include time orientation, spatial orientation, memory registration, attention and calculation, memory recall, language, and space-time configuration. The scores of the MMSE test result were divided into normal and abnormal in 15 percentile standards by using the age and education criteria. “Cognitive impairment” was defined as an adjusted MMSE score in the range of 10–23, with less than 10 for “severe dementia” (Creavin et al., 2016). The Seoul Neuropsychological Screening Battery (SNSB) test results were evaluated to determine four cognitive statuses (i.e., memory, language skills, space-time functions, and executive function) and provided percentile score adjusted for age, sex, and education (Ahn et al., 2010). The SNSB test incorporated the MMSE, Clinical Dementia Rating (Juva et al., 1994), and Global Deterioration Scale (Reisberg et al., 1982).

Assessment of Diabetic Retinopathy Stage and Retinal Ganglion Cell Layer Thickness

Images from laser-scanning ophthalmoscopy (Optos Plc., Dunfermline, United Kingdom) and fluorescein angiography (Heidelberg Engineering, Heidelberg, Germany) were reviewed to diagnose DR and to determine the DR stage. Non-proliferative vs. proliferative DR stages were identified based on the presence of preretinal angiogenesis using fundus fluorescein angiography (Group, 1991). The examination of the GCL thickness in the macula was performed using spectral-domain OCT (Spectralis®, Heidelberg Engineering, Heidelberg, Germany), and the retinal-segmentation software accompanying this device was used to identify and quantify the average thickness of each retinal neuronal layer (Reisberg et al., 1982). The GCL failed was defined as the layer between the ganglion cell-inner plexiform layer (Mwanza et al., 2011), and the GCL thickness was measured in the foveal area (a 1-mm zone centered on the fovea) and in the inferior, superior, nasal, and temporal subfields (a 1- to 3-mm perifoveal zone), as defined by the Early Treatment DR Study. The mean GCL thickness of each subfield in both eyes was analyzed. However, if one eye met exclusion criteria, only GCL parameters in the opposite eye were analyzed. Two trained retinal ophthalmologists (with 8 and 18 years of experience, respectively) measured OCT findings, and the averaged values were used for statistical analyses.

MRI Protocol

All scans were acquired with a 3.0-T scanner (Achieva; Philips Healthcare, Best, Netherlands) with a 32-channel head coil. Head motion was minimized with restraining foam pads. The MRI imaging protocol included T2-weighted (repetition time [TR]/echo time [TE]: 2,800–3,000/80–100 ms; field of view (FOV): 230–240 mm; section thickness: 5 mm; matrix: 256 × 256), T2*-weighted gradient-echo images (TR/TE: 500–1,000/15–25 ms; section thickness: 5 mm; matrix: 256 × 256), and fluid-attenuated inversion recovery (FLAIR) images (TR/TE: 9,000–10,000/110–125 ms; FOV: 240 mm; section thickness: 5 mm; matrix: 256 × 256).

Assessment of Cerebral Small Vessel Disease Markers

The PVS was identified and rated on axial, T2-weighted, MR images based on the previously reported criteria (Wardlaw et al., 2013) by two trained neuroradiologists (with 15 and 8 years of experience, respectively) who were blinded to clinical information. A 4-point visual rating scale (0 = absent PVS, 1 = less than 10 PVS, 2 = 11–20 PVS, 3 = 21–40 PVS, and 4 = more than 40 PVS) was applied to the basal ganglia (BG) and the centrum semiovale (CSO) (Wardlaw et al., 2013). The number of MRI-visible PVS was counted in each hemisphere, and the highest score was recorded. The final consensus rating scale results were used for the analysis.

Lacunes were defined as lesions larger than 3 mm and less than 15 mm in the subcortical area with high signal intensity on the T2-weighted images and a perilesional halo on FLAIR imaging (Wardlaw et al., 2013). The number of lacunes was manually counted in the brain stem, BG, thalamus, pons, and cerebral white matter. Microbleeds were defined as being 10 mm or less in diameter on axial sections of the T2-weighted gradient-echo images, using the previously reported criteria (Greenberg et al., 2009). The total number of microbleeds was counted.

The WMH was defined as a hyperintense white matter lesion on FLAIR images based on the STRIVE (STandards for ReportIng Vascular changes on nEuroimaging) criteria (Wardlaw et al., 2013) and was graded according to the Fazekas scale as deep WMH (0 = absent, 1 = punctate, 2 = early confluent, and 3 = confluent) or periventricular WMH (0 = absent, 1 = caps or pencil-thin lining, 2 = smooth halo, and 3 = irregular WMH extending into the deep white matter) (Fazekas et al., 1987; Wardlaw et al., 2013). The total Fazekas score was calculated by adding together the periventricular and deep WMH scores (0–6) (Park et al., 2018).

The total cerebral SVD score was calculated for each patient using an ordinal scale ranging from 0 to 4 by counting the presence of each of the four MRI features of cerebral SVD, with a score of 1 point indicating the presence of high BG-PVS severity (a score ≥ 2), lacunes, microbleeds, and severe WMH (total Fazekas score > 3) (Staals et al., 2014). Lacunes, microbleeds, and WMHs were rated through consensus discussions between the two trained neuroradiologists.

Statistical Analysis

Baseline clinical and imaging characteristics were compared between non-proliferative DR and proliferative DR groups using the chi-squared or Fisher’s exact test for categorical variables or the independent t-test or Mann–Whitney U test for continuous variables, according to normality. Interobserver agreements for imaging evaluations were calculated using the Cohen’s kappa index. In addition, correlations between the BG-PVS severity and the total cerebral SVD score with the MMSE score and GCL parameters were evaluated using the Spearman’s correlation coefficient analysis. The multivariable logistic regression analysis using the enter method was performed to determine the effect of selected cerebro-ophthalmic parameters (with P-values of < 0.10 in previous comparison and correlation tests) on the final outcome, i.e., DR progression to the proliferative stage. An identical process was performed to predict cognitive impairment (adjusted score of < 24 on the MMSE).

All statistical analyses were performed using the Statistical Package for the Social Sciences (SPSS/IBM Corporation; Chicago, IL, United States) version 23.0 software for Windows. A P-value < 0.05 was considered statistically significant.

Results

Among 81 patients, 51 patients were classified as having non-proliferative DR and 30 were classified as having proliferative DR.

Interobserver Agreement of Perivascular Space Evaluations

The inter-rater reliabilities were good for both BG-PVS (κ = 0.76) and CSO-PVS (κ = 0.73).

Clinical, Cerebral Small Vessel Disease, and Retinal Ganglion Cell Layer According to Diabetic Retinopathy Stage

The clinical and neuropsychological characteristics of the non-proliferative and proliferative DR groups are summarized in Table 1. There were no significant differences between the clinical or neuropsychological characteristics of these two groups.

TABLE 1
www.frontiersin.org

Table 1. Clinical characteristics and neuropsychological function according to the diabetic retinopathy (DR) stage.

Brain MRI and retinal OCT imaging findings for patients in the two DR groups are summarized in Table 2. Patients in the proliferative DR group had a significantly higher BG-PVS severity (P = 0.012) and a significantly higher total cerebral SVD score (P = 0.035) than patients in the non-proliferative DR group. Patients with proliferative DR also showed significantly reduced GCL thicknesses in the inferior (P = 0.027), superior (P = 0.046), and temporal (P = 0.038) subfields compared to patients with non-proliferative DR. Decreased GCL thicknesses were also identified in other subfields (i.e., foveal center and nasal) in the proliferative DR group; however, these results did not reach statistical significance. Figures 1, 2 show representative cases at different stages of DR progression.

TABLE 2
www.frontiersin.org

Table 2. Cerebral and ophthalmic imaging findings according to the DR stage.

FIGURE 1
www.frontiersin.org

Figure 1. (A) A 83-year-old female with non-proliferative diabetic retinopathy (DR). Patchy areas of capillary dropout with microaneurysms were observed on a montage fundus fluorescein angiography (FFA) image. The axial T2-weighted MRI images show grade 0 perivascular space in the basal ganglia (BG-PVS). On fluid-attenuated inversion recovery (FLAIR) images, the total Fazekas score was 2. The total cerebral small vessel disease (SVD) score was 0. (B) A 73-year-old female with proliferative DR. Wide-field FFA image shows the presence of neovascularizations elsewhere around the wider areas of non-perfusion. The axial T2-weighted images show grade 4 BG-PVS. On FLAIR images, the total Fazekas score was 3. The total cerebral SVD score was 1.

FIGURE 2
www.frontiersin.org

Figure 2. Optical coherence tomographic (OCT) images of (A) a 62-year-old female with non-proliferative DR and of (B) a 63-year-old female with proliferative DR. A cross-sectional B-scan of foveal center within the 6 × 6-mm macular cubes is shown in the upper row. The segmentation lines on the B-scan represent the inner retinal nerve fiber layer (RNFL) and outer inner plexiform layer (IPL) boundaries of the slab of ganglion cell layer (GCL). The increased irregularity of the GCL boundaries is observed in (B) when compared to (A). En-face OCT image with color topography of each GCL thickness corresponding to the Early Treatment Diabetic Retinopathy Study (ETDRS) circles is presented in the bottom row. A decrease in GCL thickness is presented in panel B compared to panel A, as a significant color change in the 1- to 3-mm perifoveal zone.

Correlation of BG-PVS Severity and Total Cerebral Small Vessel Disease Score With MMSE Score

The results of the correlation analysis for cerebral SVD makers and neuropsychological function are summarized in Supplementary Table 1. The BG-PVS severity was negatively correlated with the MMSE score (P = 0.007), while the total cerebral SVD score did not show a significant correlation with this score (P = 0.23). Figure 3 shows a corresponding scatterplot. In addition, BG-PVS was negatively correlated with total visuospatial score (P = 0.039), while it was positively correlated with CDR (P = 0.043) and GDS (P = 0.004). Total cerebral SVD score presented a negative correlation with frontal executive function (P < 0.001) and a positive correlation with GDS (P = 0.035).

FIGURE 3
www.frontiersin.org

Figure 3. Correlation plot between the BG-PVS severity and the total cerebral SVD score with the Mini-Mental State Examination (MMSE) score.

The CSO-PVS showed a negative correlation with MMSE (P = 0.005). Total WMH was negatively correlated with total attention (P = 0.015), visuospatial (P = 0.023), memory (P = 0.025), and frontal executive (P < 0.001) function scores.

Correlation of Retinal Ganglion Cell Layer Thickness With BG-PVS Severity and Total Cerebral Small Vessel Disease Score

The correlation analysis results between the GCL thickness and the BG-PVS severity and total cerebral SVD score are summarized in Supplementary Table 2. In this analysis, the GCL thickness showed significant negative correlations with the BG-PVS severity (P = 0.041 for inferior; P = 0.006 for superior; and P = 0.038 for temporal GCL, respectively) and the total cerebral SVD score (P < 0.001 for inferior; P = 0.004 for superior; and P = 0.001 for nasal and temporal GCL, respectively). Figure 4 shows corresponding scatterplots.

FIGURE 4
www.frontiersin.org

Figure 4. Correlation plots between (A) BG-PVS severity and GCL thickness and (B) total cerebral SVD score and GCL thickness.

Supplementary Table 2 also summarizes the results of correlation analyses for the thickness of retinal GCL with other cerebral SVD-related markers. GCL thickness showed significant negative correlations with BG-PVS (P = 0.041 for inferior; P = 0.006 for superior; and P = 0.038 for temporal GCL, respectively), total number of lacunes (P < 0.001 for inferior and superior; and P = 0.012 for temporal GCL, respectively), total WMH (P < 0.001 for inferior; P = 0.005 for superior; P = 0.015 for nasal; and P = 0.042 for temporal GCL, respectively), and total CSVD (P < 0.001 for inferior; P = 0.004 for superior; and P = 0.001 for nasal and temporal GCL, respectively). CSO-PVS was not significantly associated with GCL thickness.

Logistic Regression Analysis for Diabetic Retinopathy Stage Progression and Cognitive Impairment

The determinants among cerebro-ophthalmic parameters obtained by the logistic regression analysis are presented in Table 3. Possible predictors (i.e., age, male sex, underlying cardiovascular diseases, BG-PVS, CSO-PVS, total number of lacunes and WMH, total cerebral SVD score, and inferior and superior GCL thickness) that were significant in previous comparison and correlation tests were included in the multivariable logistic regression models with the enter method. According to the developed model, a higher BG-PVS score (OR = 2.306, P = 0.017) was associated with a higher risk of pathological angiogenesis in DR, while old age (OR = 0.138, P = 0.004) was associated with lower risk. Old age (OR = 2.620, P = 0.012) and higher BG-PVS score (OR = 1.716, P = 0.044) were associated with a higher risk of cognitive impairment, while thicker inferior GCL (OR = 0.669, P = 0.011) was associated with a lower risk.

TABLE 3
www.frontiersin.org

Table 3. Multivariable logistic regression models for the determination of significant cerebro-ophthalmic parameters for DR progression (proliferative stage) and for cognitive impairment (adjusted MMSE score < 24).

Discussion

In this study, we assessed the relationships between cerebral SVD severity and retinal GCL parameters based on the DR stage. Higher BG-PVS severity, higher total cerebral SVD score, and reduced retinal GCL thicknesses in the inferior, superior, and temporal subfields were identified in the proliferative DR group compared with the non-proliferative DR group. In addition, the BG-PVS severity was negatively correlated with the MMSE score, while the retinal GCL thickness was negatively correlated with the BG-PVS severity and the total cerebral SVD score. Taken together, these findings suggest that the BG-PVS severity reflects the retinal neurodegeneration that precedes pathological angiogenesis during DR progression. Moreover, the BG-PVS severity may be a potential biomarker of cognitive decline in patients with DR. Finally, our results provide possible clinical evidence of a link between the retina and the brain on the cognitive impairment in patients with DR.

While a relationship between enlarged PVS and cerebral SVD has been noted (Doubal et al., 2010), PVS have also been proposed to be a part of a macroscopic clearance mechanism called the glymphatic system (Jessen et al., 2015), which may be related to the development of neurodegenerative diseases. A previous study has shown that T2DM suppresses the clearance of interstitial fluid, leading to the impairment of glymphatic system (Jiang et al., 2017). Although the mechanism underlying enlarged PVS is not completely understood, previous studies have shown that associations of disease with visible PVS differ based on their location (Charidimou et al., 2017; Chen et al., 2019). The BG-PVS burden has been associated with hypertensive angiopathy, systemic markers of inflammation, lacunar stroke, and vascular cognitive impairment (Hansen et al., 2015; Chen et al., 2019), while the CSO-PVS burden has been associated with cerebral amyloid angiopathy and Alzheimer’s disease (van Veluw et al., 2016; Chen et al., 2019; Sepehrband et al., 2021). In our study, the BG-PVS severity was associated with a higher DR grade, suggesting that pathological angiogenesis in DR is associated with the systemic markers of inflammation and vascular cognitive impairment. Although a previous meta-analysis has shown that enlarged PVS is not associated with T2DM (Francis et al., 2019), our discrepant results may have been related to the fact that we assessed enlarged PVS based on the DR stage, which represents the diabetic neurovascular remodeling process. Proliferative DR is present only in 10–12% of patients with T2DM (Thomas et al., 2015), and previous studies have not evaluated whether the PVS burden is increased with higher DR severity. Our results provide a better understanding of the pathogenesis underlying brain manifestations during the progression of DR.

Previous studies have shown an increased burden of lacunes in patients with T2DM and DR (Inoue et al., 1998; Sanahuja et al., 2016). The relationship between diabetes and microbleeds, however, remains unclear (Umemura et al., 2017), and the observational studies have reported discordant results concerning the relationship between diabetes and WMH (van Elderen et al., 2010; Sanahuja et al., 2016; de Bresser et al., 2018). In our study, these cerebral SVD markers did not show significant differences based on DR grade. The total cerebral SVD score has recently been proposed as a potential scoring system for the DR stage and has been shown to better represent the overall effects of SVD on the brain than one or two features independently (Cuadrado-Godia et al., 2018). Our results also suggest that the total cerebral SVD score may well reflect the DR stage.

Furthermore, our results suggest a potential effect of BG-PVS severity on the early stages of cognitive impairment in patients with DR. Recent studies have identified an association between BG-PVS severity and cognitive impairment in patients with Parkinson’s disease (PD) and stroke (Arba et al., 2018; Park et al., 2019). These findings suggest that the BG-PVS severity may be a novel indicator of the early pathophysiological processes leading to cognitive impairment in the brain. Considering a recent study showing that PVS from other brain regions may also be associated with PD (Donahue et al., 2021), future studies should consider various anatomical regions for PVS evaluation.

Similar to the glymphatic system of the brain, the existence of neurovascular units in the retina has also been proposed (Metea and Newman, 2007; Newman, 2013). Blood-retinal barrier disruption and retinal neurodegeneration (Bek, 2017; Duh et al., 2017) are more significantly activated by low-grade inflammation (Du et al., 2013; Altmann and Schmidt, 2018) and are early events during the pathogenesis of DR. Glial activation and neuronal apoptosis, both being the prominent features of retinal neurodegeneration (Simó and Hernández, 2014), result in decreased retinal GCL thickness. An association between a reduced retinal GCL thickness and the SVD burden in the brain was first identified in this study through a correlation analysis between the BG-PVS severity and the total cerebral SVD score. These findings suggest that DR progression and cognitive impairment in diabetic patients may share specific mechanisms mainly related to systemic microangiopathy. This common mechanism appears to be associated with the changes in GCL thickness and BG-PVS severity. Further longitudinal studies are necessary to validate the sequence of these changes.

It should be noted that there were several limitations to our study. First, this study was based on the retrospective, single-institutional data with a relatively small size without complete control of covariates and external validation potentially limiting the generalizability of our results. However, to minimize the selection bias, we ensured that the basic characteristics such as age and sex and the underlying conditions of the subgroups were comparable. Second, the lack of an age-matched control group without DR prevents the extrapolation of these results to the general diabetic population. Future studies with the inclusion of age-matched control groups are warranted. Third, variables such as brain volume or body mass index, which have shown an association with PVS (Barisano et al., 2021; Liu et al., 2021), were not analyzed in this study. Further studies should also include these variables in a larger dataset.

Conclusion

Our results demonstrate that the BG-PVS severity may be a novel imaging biomarker in patients with DR, reflecting progression in the DR stage and cognitive impairment. Furthermore, our identification of a possible association between a decreased GCL thickness and the BG-PVS severity provides potential clinical evidence of a link between retinal and brain changes during diabetic neurodegeneration.

Data Availability Statement

The corresponding author had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. The raw data supporting the conclusions of this article will be made available by the correspondence, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the Severance Institutional Review Board (IRB approval number: 2020-3812-001). Written informed consent for participation was not required for this study in accordance with the National Legislation and the Institutional Requirements.

Author Contributions

EYC and YWP contributed to the concept, design, and drafting of the manuscript and performed the statistical analysis. EYC, YWP, MK, CSL, SSA, JK, and S-KL contributed to the acquisition, analysis, and interpretation of data. EYC, YWP, and ML contributed to the critical revision of the manuscript for important intellectual content, administrative, technical, material support, and supervision. YWP obtained funding. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2020R1I1A1A0107164811) and by the “Team Science Award” of Yonsei University College of Medicine (6-2021-0009).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

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

Abbreviations

BG, basal ganglia; CDR, clinical dementia rating scale; CSO, centrum semiovale; CSVD, cerebral small vessel disease; DR, diabetic retinopathy; FOV, field of view; FLAIR, fluid-attenuated inversion recovery; GCL, ganglion cell layer; GDS, global deterioration scale; MMSE, Mini-Mental State Examination; OCT, optical coherence tomography; OR, odds ratio; PVS, perivascular space; SNSB, Seoul neuropsychological screening battery; SVD, small vessel disease; TE, echo time; TR, repetition time; T2DM, type 2 diabetes mellitus; WMH, white matter hyperintensity.

References

Ahn, H. J., Chin, J., Park, A., Lee, B. H., Suh, M. K., Seo, S. W., et al. (2010). Seoul Neuropsychological Screening Battery-dementia version (SNSB-D): a useful tool for assessing and monitoring cognitive impairments in dementia patients. J. Korean Med. Sci. 25, 1071–1076. doi: 10.3346/jkms.2010.25.7.1071

PubMed Abstract | CrossRef Full Text | Google Scholar

Altmann, C., and Schmidt, M. H. H. (2018). The Role of Microglia in Diabetic Retinopathy: Inflammation, Microvasculature Defects and Neurodegeneration. Int. J. Mol. Sci. 19:110. doi: 10.3390/ijms19010110

PubMed Abstract | CrossRef Full Text | Google Scholar

Aminian, A., Zajichek, A., Arterburn, D. E., Wolski, K. E., Brethauer, S. A., Schauer, P. R., et al. (2020). Predicting 10-Year Risk of End-Organ Complications of Type 2 Diabetes With and Without Metabolic Surgery: A Machine Learning Approach. Diabetes Care 43, 852–859. doi: 10.2337/dc19-2057

PubMed Abstract | CrossRef Full Text | Google Scholar

Arba, F., Quinn, T. J., Hankey, G. J., Lees, K. R., Wardlaw, J. M., Ali, M., et al. (2018). Enlarged perivascular spaces and cognitive impairment after stroke and transient ischemic attack. Int. J. Stroke 13, 47–56. doi: 10.1177/1747493016666091

PubMed Abstract | CrossRef Full Text | Google Scholar

Banerjee, G., Kim, H. J., Fox, Z., Jäger, H. R., Wilson, D., Charidimou, A., et al. (2017). MRI-visible perivascular space location is associated with Alzheimer’s disease independently of amyloid burden. Brain 140, 1107–1116. doi: 10.1093/brain/awx003

PubMed Abstract | CrossRef Full Text | Google Scholar

Barisano, G., Sheikh-Bahaei, N., Law, M., Toga, A. W., and Sepehrband, F. (2021). Body mass index, time of day and genetics affect perivascular spaces in the white matter. J. Cereb. Blood Flow Metab. 41, 1563–1578. doi: 10.1177/0271678X20972856

PubMed Abstract | CrossRef Full Text | Google Scholar

Bek, T. (2017). Diameter Changes of Retinal Vessels in Diabetic Retinopathy. Curr. Diab. Rep. 17:82. doi: 10.1007/s11892-017-0909-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Charidimou, A., Boulouis, G., Pasi, M., Auriel, E., Van Etten, E. S., Haley, K., et al. (2017). MRI-visible perivascular spaces in cerebral amyloid angiopathy and hypertensive arteriopathy. Neurology 88, 1157–1164. doi: 10.1212/WNL.0000000000003746

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Wang, J., Shan, Y., Cai, W., Liu, S., Hu, M., et al. (2019). Cerebral small vessel disease: neuroimaging markers and clinical implication. J. Neurol. 266, 2347–2362. doi: 10.1007/s00415-018-9077-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Creavin, S. T., Wisniewski, S., Noel-Storr, A. H., Trevelyan, C. M., Hampton, T., Rayment, D., et al. (2016). Mini-Mental State Examination (MMSE) for the detection of dementia in clinically unevaluated people aged 65 and over in community and primary care populations. Cochrane Database Syst. Rev. 1:Cd011145. doi: 10.1002/14651858.CD011145.pub2

PubMed Abstract | CrossRef Full Text | Google Scholar

Crosby-Nwaobi, R., Sivaprasad, S., and Forbes, A. (2012). A systematic review of the association of diabetic retinopathy and cognitive impairment in people with Type 2 diabetes. Diabetes Res. Clin. Pract. 96, 101–110. doi: 10.1016/j.diabres.2011.11.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Cuadrado-Godia, E., Dwivedi, P., Sharma, S., Ois Santiago, A., Roquer Gonzalez, J., Balcells, M., et al. (2018). Cerebral Small Vessel Disease: A Review Focusing on Pathophysiology, Biomarkers, and Machine Learning Strategies. J. Stroke 20, 302–320. doi: 10.5853/jos.2017.02922

PubMed Abstract | CrossRef Full Text | Google Scholar

de Bresser, J., Kuijf, H. J., Zaanen, K., Viergever, M. A., Hendrikse, J., and Biessels, G. J. (2018). White matter hyperintensity shape and location feature analysis on brain MRI; proof of principle study in patients with diabetes. Sci. Rep. 8:1893. doi: 10.1038/s41598-018-20084-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Donahue, E. K., Murdos, A., Jakowec, M. W., Sheikh-Bahaei, N., Toga, A. W., Petzinger, G. M., et al. (2021). Global and Regional Changes in Perivascular Space in Idiopathic and Familial Parkinson’s Disease. Mov. Disord. 36, 1126–1136. doi: 10.1002/mds.28473

PubMed Abstract | CrossRef Full Text | Google Scholar

Doubal, F. N., Maclullich, A. M., Ferguson, K. J., Dennis, M. S., and Wardlaw, J. M. (2010). Enlarged perivascular spaces on MRI are a feature of cerebral small vessel disease. Stroke 41, 450–454. doi: 10.1161/STROKEAHA.109.564914

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, Y., Veenstra, A., Palczewski, K., and Kern, T. S. (2013). Photoreceptor cells are major contributors to diabetes-induced oxidative stress and local inflammation in the retina. Proc. Natl. Acad. Sci. U. S. A. 110, 16586–16591. doi: 10.1073/pnas.1314575110

PubMed Abstract | CrossRef Full Text | Google Scholar

Duh, E. J., Sun, J. K., and Stitt, A. W. (2017). Diabetic retinopathy: current understanding, mechanisms, and treatment strategies. JCI Insight 2:e93751. doi: 10.1172/jci.insight.93751

PubMed Abstract | CrossRef Full Text | Google Scholar

Duperron, M. G., Tzourio, C., Sargurupremraj, M., Mazoyer, B., Soumaré, A., Schilling, S., et al. (2018). Burden of Dilated Perivascular Spaces, an Emerging Marker of Cerebral Small Vessel Disease, Is Highly Heritable. Stroke 49, 282–287. doi: 10.1161/STROKEAHA.117.019309

PubMed Abstract | CrossRef Full Text | Google Scholar

Ergul, A., Kelly-Cobbs, A., Abdalla, M., and Fagan, S. C. (2012). Cerebrovascular complications of diabetes: focus on stroke. Endocr. Metab. Immune Disord. Drug Targets 12, 148–158. doi: 10.2174/187153012800493477

PubMed Abstract | CrossRef Full Text | Google Scholar

Fazekas, F., Chawluk, J. B., Alavi, A., Hurtig, H. I., and Zimmerman, R. A. (1987). MR signal abnormalities at 1.5 T in Alzheimer’s dementia and normal aging. Am. J. Roentgenol. 149, 351–356. doi: 10.2214/ajr.149.2.351

PubMed Abstract | CrossRef Full Text | Google Scholar

Francis, F., Ballerini, L., and Wardlaw, J. M. (2019). Perivascular spaces and their associations with risk factors, clinical disorders and neuroimaging features: A systematic review and meta-analysis. Int. J. Stroke 14, 359–371. doi: 10.1177/1747493019830321

PubMed Abstract | CrossRef Full Text | Google Scholar

Greenberg, S. M., Vernooij, M. W., Cordonnier, C., Viswanathan, A., Al-Shahi Salman, R., Warach, S., et al. (2009). Cerebral microbleeds: a guide to detection and interpretation. Lancet Neurol. 8, 165–174. doi: 10.1016/S1474-4422(09)70013-4

CrossRef Full Text | Google Scholar

Group, E. T. D. R. S. R. (1991). Grading diabetic retinopathy from stereoscopic color fundus photographs–an extension of the modified Airlie House classification. ETDRS report number 10. Early Treatment Diabetic Retinopathy Study Research Group. Ophthalmology 98, 786–806. doi: 10.1016/S0161-6420(13)38012-9

CrossRef Full Text | Google Scholar

Gupta, P., Gan, A. T. L., Man, R. E. K., Fenwick, E. K., Sabanayagam, C., Mitchell, P., et al. (2019). Association between diabetic retinopathy and incident cognitive impairment. Br. J. Ophthalmol. 103, 1605–1609. doi: 10.1136/bjophthalmol-2018-312807

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, C., Jo, S. A., Jo, I., Kim, E., Park, M. H., and Kang, Y. (2008). An adaptation of the Korean mini-mental state examination (K-MMSE) in elderly Koreans: demographic influence and population-based norms (the AGE study). Arch. Gerontol. Geriatr. 47, 302–310. doi: 10.1016/j.archger.2007.08.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Hansen, T., Cain, J., Thomas, O., and Jackson, A. (2015). Dilated perivascular spaces in the Basal Ganglia are a biomarker of small-vessel disease in a very elderly population with dementia. Am. J. Neuroradiol. 36, 893–898. doi: 10.3174/ajnr.A4237

PubMed Abstract | CrossRef Full Text | Google Scholar

Hurford, R., Charidimou, A., Fox, Z., Cipolotti, L., Jager, R., and Werring, D. J. (2014). MRI-visible perivascular spaces: relationship to cognition and small vessel disease MRI markers in ischaemic stroke and TIA. J. Neurol. Neurosurg. Psychiatry 85, 522–525. doi: 10.1136/jnnp-2013-305815

PubMed Abstract | CrossRef Full Text | Google Scholar

Inoue, T., Fushimi, H., Yamada, Y., Kubo, M., Udaka, F., and Kameyama, M. (1998). The changes of lacunar state during a 5-year period in NIDDM. Diabetes Res. Clin. Pract. 42, 155–160. doi: 10.1016/S0168-8227(98)00106-5

CrossRef Full Text | Google Scholar

Jessen, N. A., Munk, A. S., Lundgaard, I., and Nedergaard, M. (2015). The Glymphatic System: A Beginner’s Guide. Neurochem. Res. 40, 2583–2599.

Google Scholar

Jiang, Q., Zhang, L., Ding, G., Davoodi-Bojd, E., Li, Q., Li, L., et al. (2017). Impairment of the glymphatic system after diabetes. J. Cereb. Blood Flow Metab. 37, 1326–1337. doi: 10.1007/s11064-015-1581-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Juva, K., Sulkava, R., Erkinjuntti, T., Ylikoski, R., Valvanne, J., and Tilvis, R. (1994). Staging the severity of dementia: comparison of clinical (CDR, DSM-III-R), functional (ADL, IADL) and cognitive (MMSE) scales. Acta Neurol. Scand. 90, 293–298. doi: 10.1177/0271678X16654702

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, S., Hou, B., You, H., Zhang, Y., Zhu, Y., Ma, C., et al. (2021). The Association Between Perivascular Spaces and Cerebral Blood Flow, Brain Volume, and Cardiovascular Risk. Front. Aging Neurosci. 13:599724. doi: 10.1111/j.1600-0404.1994.tb02724.x

PubMed Abstract | CrossRef Full Text | Google Scholar

London, A., Benhar, I., and Schwartz, M. (2013). The retina as a window to the brain-from eye research to CNS disorders. Nat. Rev. Neurol. 9, 44–53.

Google Scholar

Metea, M. R., and Newman, E. A. (2007). Signalling within the neurovascular unit in the mammalian retina. Exp. Physiol. 92, 635–640. doi: 10.1038/nrneurol.2012.227

PubMed Abstract | CrossRef Full Text | Google Scholar

Mwanza, J. C., Oakley, J. D., Budenz, D. L., Chang, R. T., Knight, O. J., and Feuer, W. J. (2011). Macular ganglion cell-inner plexiform layer: automated detection and thickness reproducibility with spectral domain-optical coherence tomography in glaucoma. Invest. Ophthalmol. Vis. Sci. 52, 8323–8329. doi: 10.1113/expphysiol.2006.036376

PubMed Abstract | CrossRef Full Text | Google Scholar

Newman, E. A. (2013). Functional hyperemia and mechanisms of neurovascular coupling in the retinal vasculature. J. Cereb. Blood Flow Metab. 33, 1685–1695. doi: 10.1167/iovs.11-7962

PubMed Abstract | CrossRef Full Text | Google Scholar

Ng, D. S., Chiang, P. P., Tan, G., Cheung, C. G., Cheng, C. Y., Cheung, C. Y., et al. (2016). Retinal ganglion cell neuronal damage in diabetes and diabetic retinopathy. Clin. Exp. Ophthalmol. 44, 243–250. doi: 10.1038/jcbfm.2013.145

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, M., Moon, Y., Han, S.-H., Kim, H. K., and Moon, W.-J. (2018). Myelin loss in white matter hyperintensities and normal-appearing white matter of cognitively impaired patients: a quantitative synthetic magnetic resonance imaging study. Eur. Radiol. 29, 4914–4921. doi: 10.1111/ceo.12724

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, Y. W., Shin, N. Y., Chung, S. J., Kim, J., Lim, S. M., Lee, P. H., et al. (2019). Magnetic Resonance Imaging-Visible Perivascular Spaces in Basal Ganglia Predict Cognitive Decline in Parkinson’s Disease. Mov. Disord. 34, 1672–1679. doi: 10.1007/s00330-018-5836-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Petzold, A. (2016). Retinal glymphatic system: an explanation for transient retinal layer volume changes? Brain 139, 2816–2819. doi: 10.1002/mds.27798

PubMed Abstract | CrossRef Full Text | Google Scholar

Potter, G. M., Doubal, F. N., Jackson, C. A., Chappell, F. M., Sudlow, C. L., Dennis, M. S., et al. (2015). Enlarged perivascular spaces and cerebral small vessel disease. Int. J. Stroke 10, 376–381. doi: 10.1093/brain/aww239

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramirez, J., Berezuk, C., Mcneely, A. A., Gao, F., Mclaurin, J., Black, S. E. J. C., et al. (2016). Imaging the perivascular space as a potential biomarker of neurovascular and neurodegenerative diseases. Cell. Mol. Neurobiol. 36, 289–299. doi: 10.1111/ijs.12054

PubMed Abstract | CrossRef Full Text | Google Scholar

Reisberg, B., Ferris, S. H., De Leon, M. J., and Crook, T. (1982). The Global Deterioration Scale for assessment of primary degenerative dementia. Am. J. Psychiatry 139, 1136–1139. doi: 10.1007/s10571-016-0343-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanahuja, J., Alonso, N., Diez, J., Ortega, E., Rubinat, E., Traveset, A., et al. (2016). Increased Burden of Cerebral Small Vessel Disease in Patients With Type 2 Diabetes and Retinopathy. Diabetes Care 39, 1614–1620. doi: 10.1176/ajp.139.9.1136

PubMed Abstract | CrossRef Full Text | Google Scholar

Sepehrband, F., Barisano, G., Sheikh-Bahaei, N., Choupan, J., Cabeen, R. P., Lynch, K. M., et al. (2021). Volumetric distribution of perivascular space in relation to mild cognitive impairment. Neurobiol. Aging 99, 28–43. doi: 10.2337/dc15-2671

PubMed Abstract | CrossRef Full Text | Google Scholar

Simó, R., and Hernández, C. (2014). Neurodegeneration in the diabetic eye: new insights and therapeutic perspectives. Trends Endocrinol. Metab. 25, 23–33. doi: 10.1016/j.neurobiolaging.2020.12.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Sohn, E. H., Van Dijk, H. W., Jiao, C., Kok, P. H., Jeong, W., Demirkaya, N., et al. (2016). Retinal neurodegeneration may precede microvascular changes characteristic of diabetic retinopathy in diabetes mellitus. Proc. Natl. Acad. Sci. U. S. A. 113, E2655–E2664. doi: 10.1016/j.tem.2013.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Staals, J., Makin, S. D., Doubal, F. N., Dennis, M. S., and Wardlaw, J. M. (2014). Stroke subtype, vascular risk factors, and total MRI brain small-vessel disease burden. Neurology 83, 1228–1234. doi: 10.1073/pnas.1522014113

PubMed Abstract | CrossRef Full Text | Google Scholar

Stitt, A. W., Curtis, T. M., Chen, M., Medina, R. J., Mckay, G. J., Jenkins, A., et al. (2016). The progress in understanding and treatment of diabetic retinopathy. Prog. Retin. Eye Res. 51, 156–186. doi: 10.1212/WNL.0000000000000837

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomas, R. L., Dunstan, F. D., Luzio, S. D., Chowdhury, S. R., North, R. V., Hale, S. L., et al. (2015). Prevalence of diabetic retinopathy within a national diabetic retinopathy screening service. Br. J. Ophthalmol. 99, 64–68. doi: 10.1016/j.preteyeres.2015.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomas, R. L., Halim, S., Gurudas, S., Sivaprasad, S., and Owens, D. R. (2019). IDF Diabetes Atlas: A review of studies utilising retinal photography on the global prevalence of diabetes related retinopathy between 2015 and 2018. Diabetes Res. Clin. Pract. 157:107840. doi: 10.1136/bjophthalmol-2013-304017

PubMed Abstract | CrossRef Full Text | Google Scholar

Umemura, T., Kawamura, T., and Hotta, N. (2017). Pathogenesis and neuroimaging of cerebral large and small vessel disease in type 2 diabetes: A possible link between cerebral and retinal microvascular abnormalities. J. Diabetes Investig. 8, 134–148. doi: 10.1016/j.diabres.2019.107840

PubMed Abstract | CrossRef Full Text | Google Scholar

van Elderen, S. G., De Roos, A., De Craen, A. J., Westendorp, R. G., Blauw, G. J., Jukema, J. W., et al. (2010). Progression of brain atrophy and cognitive decline in diabetes mellitus: a 3-year follow-up. Neurology 75, 997–1002. doi: 10.1111/jdi.12545

PubMed Abstract | CrossRef Full Text | Google Scholar

van Veluw, S. J., Biessels, G. J., Bouvy, W. H., Spliet, W. G., Zwanenburg, J. J., Luijten, P. R., et al. (2016). Cerebral amyloid angiopathy severity is linked to dilation of juxtacortical perivascular spaces. J. Cereb. Blood Flow Metab. 36, 576–580. doi: 10.1212/WNL.0b013e3181f25f06

PubMed Abstract | CrossRef Full Text | Google Scholar

Wardlaw, J., Smith, E., Biessels, G., Cordonnier, C., Fazekas, F., Frayne, R., et al. (2013). Neuroimaging standards for research into small vessel disease and its contribution to ageing and neurodegeneration. Lancet Neurol. 12, 822–838. doi: 10.1177/0271678X15620434

PubMed Abstract | CrossRef Full Text | Google Scholar

Woerdeman, J., Van Duinkerken, E., Wattjes, M. P., Barkhof, F., Snoek, F. J., Moll, A. C., et al. (2014). Proliferative retinopathy in type 1 diabetes is associated with cerebral microbleeds, which is part of generalized microangiopathy. Diabetes Care 37, 1165–1168. doi: 10.1016/S1474-4422(13)70124-8

CrossRef Full Text | Google Scholar

Keywords: basal ganglia, cerebral small vessel disease, cognitive decline, diabetes mellitus, diabetic retinopathy, ganglion cell layer, perivascular space

Citation: Choi EY, Park YW, Lee M, Kim M, Lee CS, Ahn SS, Kim J and Lee S-K (2021) Magnetic Resonance Imaging-Visible Perivascular Spaces in the Basal Ganglia Are Associated With the Diabetic Retinopathy Stage and Cognitive Decline in Patients With Type 2 Diabetes. Front. Aging Neurosci. 13:666495. doi: 10.3389/fnagi.2021.666495

Received: 10 February 2021; Accepted: 13 October 2021;
Published: 12 November 2021.

Edited by:

Rita Machado De Olivera, New University of Lisbon, Portugal

Reviewed by:

Henri Leinonen, University of Eastern Finland, Finland
Giuseppe Barisano, University of Southern California, United States
Jeiran Choupan, University of California, Los Angeles, United States

Copyright © 2021 Choi, Park, Lee, Kim, Lee, Ahn, Kim and Lee. 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: Yae Won Park, yaewonpark@yuhs.ac

Download