key: cord-326956-oz047qmf
authors: Lu, Yiping; Li, Xuanxuan; Geng, Daoying; Mei, Nan; Wu, Pu-Yeh; Huang, Chu-Chung; Jia, Tianye; Zhao, Yajing; Wang, Dongdong; Xiao, Anling; Yin, Bo
title: Cerebral Micro-Structural Changes in COVID-19 Patients – An MRI-based 3-month Follow-up Study
date: 2020-08-03
journal: EClinicalMedicine
DOI: 10.1016/j.eclinm.2020.100484
sha: 
doc_id: 326956
cord_uid: oz047qmf

BACKGROUND: Increasing evidence supported the possible neuro-invasion potential of SARS-CoV-2. However, no studies were conducted to explore the existence of the micro-structural changes in the central nervous system after infection. We aimed to identify the existence of potential brain micro-structural changes related to SARS-CoV-2. METHODS: In this prospective study, diffusion tensor imaging (DTI) and 3D high-resolution T1WI sequences were acquired in 60 recovered COVID-19 patients (56.67% male; age: 44.10 ± 16.00) and 39 age- and sex-matched non-COVID-19 controls (56.41% male; age: 45.88 ± 13.90). Registered fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD) were quantified for DTI, and an index score system was introduced. Regional volumes derived from Voxel-based Morphometry (VBM) and DTI metrics were compared using analysis of covariance (ANCOVA). Two sample t-test and Spearman correlation were conducted to assess the relationships among imaging indices, index scores and clinical information. FINDINGS: In this follow-up stage, neurological symptoms were presented in 55% COVID-19 patients. COVID-19 patients had statistically significantly higher bilateral gray matter volumes (GMV) in olfactory cortices, hippocampi, insulas, left Rolandic operculum, left Heschl's gyrus and right cingulate gyrus and a general decline of MD, AD, RD accompanied with an increase of FA in white matter, especially AD in the right CR, EC and SFF, and MD in SFF compared with non-COVID-19 volunteers (corrected p value <0.05). Global GMV, GMVs in left Rolandic operculum, right cingulate, bilateral hippocampi, left Heschl's gyrus, and Global MD of WM were found to correlate with memory loss (p value <0.05). GMVs in the right cingulate gyrus and left hippocampus were related to smell loss (p value <0.05). MD-GM score, global GMV, and GMV in right cingulate gyrus were correlated with LDH level (p value <0.05). INTERPRETATION: Study findings revealed possible disruption to micro-structural and functional brain integrity in the recovery stages of COVID-19, suggesting the long-term consequences of SARS-CoV-2. FUNDING: Shanghai Natural Science Foundation, Youth Program of National Natural Science Foundation of China, Shanghai Sailing Program, Shanghai Science and Technology Development, Shanghai Municipal Science and Technology Major Project and ZJ Lab.

We searched PubMed and the China National Knowledge Infrastructure database for articles published up to June 10, 2020, using the keywords "COVID-19", "SARS-CoV-2", and "neurological", "DTI", and "brain". A total of 29 case reports/series reports and 106 reviews, among which 8 were systematic reviews, about the neurological findings in the central nervous system. No data has been reported about the structural changes in the brain from COVID-19 patients by functional MRI in a prospective way.

In this prospective study, diffusion tensor imaging (DTI) and 3D high-resolution T1WI sequences were acquired in 60 recovered COVID-19 patients and 39 age-and sex-matched non-COVID-19 controls. We found that these recovered COVID-19 patients were more likely to have enlarged olfactory cortices, hippocampi, insulas, Heschl's gyrus, Rolandic operculum and cingulate gyrus, and a general decline of Mean Diffusivity (MD), Axial Diffusivity (AD), Radial Diffusivity (RD) accompanied with an increase of Fractional Anisotropy (FA) in white matter, especially AD in the right Coronal Radiata (CR), External Capsule (EC) and Superior Frontal-occipital Fasciculus (SFF), and MD in SFF compared with non-COVID-19 volunteers. Global Gray Matter Volume (GMV), GMVs in left Rolandic operculum, right cingulate, bilateral hippocampi, left Heschl's gyrus, and Global MD of WM were found to correlate with memory loss. GMVs in right cingulate gyrus and left hippocampus were related to smell loss. MD-GM score, global GMV, and GMV in right cingulate gyrus were correlated with Lactate Dehydrogenase (LDH) level.

Our findings revealed possible disruption to micro-structural and functional brain integrity in the recovery stages of COVID-19, suggesting neuro-invasion potential of SARS-CoV-2. This requires attention since even if the patients recover well from the pneumonia condition, the neurological changes may cause a great burden. More research in the mechanism and route of neuro-invasion of SARS-CoV-2 is expected.

Coronavirus Disease 2019 (COVID- 19) , an illness caused by the novel Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2), is an on-going viral pandemic and has spread to the whole world. To date, it has been spreading globally with nearly 4,700,000 active infections and the total death toll was over 560,000 [1] . As a member of the coronavirus family, SARS-CoV-2 shares a 77.2% amino acid identity, 72.8% sequence identity and structural similarity with Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV) [2, 3] .

With high affinity of the receptor-binding domain of Angiotensin Converting Enzyme-2 (ACE-2), SARS-CoV-2 invades human cells in the same way as SARS-CoV [4, 5] . In light of the neurological invasion of SARS-CoV proved by abundant studies, it is plausible to hypothesize that SARS-CoV-2 has the potential to attack the central nervous system (CNS) as well [6] .

Increasing evidence supported the neuro-invading potential of SARS-CoV-2. According to the first-hand evidence from Wuhan, 36 .4% of COVID-19 patients presented neurological symptoms such as dizziness, headache and impaired consciousness during hospitalization. Furthermore, the percentage and extensiveness were higher in severe patients [7] . Except the frequent olfactory and gustatory dysfunctions in mild-to-moderate COVID-19 patients recorded by 12 European hospitals, scattered cases of various neurological diseases including encephalitis, stroke, micro-hemorrhage, hemorrhage posterior reversible encephalopathy and cerebral venous embolism were also reported in hospitalized patients [8À10] . Additionally, it was documented that the specific SARS-CoV-2 RNA was detected in the cerebral-spinal fluid (CSF) of a COVID-19 patient [11] . Although no definite description of the pathological findings has been recorded, all the aforementioned evidence indicated SARS-CoV-2 was neuro-invasive just like SARS-CoV [12] .

Based on previous researches, coronaviruses can cause demyelination, neurodegeneration, and cellular senescence which accelerate brain aging and exacerbate neurodegenerative pathology [4,13À16] . However, only scattered neurological cases in COVID-19 patients were collected to document the neurological changes during the acute infection period and so far, no long-term observation has been conducted to explore possible structural changes in the CNS. Despite satisfactory recovery in the majority of COVID-19 patients, they will have a great burden if there are neurological consequences. Therefore, it was necessary to investigate the long-term impact of SARS-CoV-2 infection on the CNS, especially on the structures easily attacked by virus and the structures highly-expressing ACE-2 [3, 17] .

In contrast to the pathological methods, in-vivo Magnetic Resonance Imaging (MRI) could reflect the cerebral structures non-invasively. The possible micro-structural damage in CNS could be detected by structural MRI and diffusion tensor imaging (DTI). Axial diffusivity (AD), radial diffusivity (RD), mean diffusivity (MD) and fractional anisotropy (FA) can be calculated using track-based spatial statistics (TBSS) [18] . Together with volumetric analysis, DTI is widely used in a large scale of neuro-radiological studies to detect microstructural changes in patients with cerebral viral infections, Human Immunodeficiency Virus (HIV) and Herpes Simplex Virus (HSV), etc [19, 20] . So far, no researches have been found to describe the cerebral changes after SARS-CoV-2 infection.

Therefore, in the current study, we aimed to apply volumetric and diffusion measurements in recovered COVID-19 patients to identify the existence of potential long-term brain structural changes related to SARS-CoV-2, which could provide better insights to understand the impact of SARS-CoV-2 on the CNS.

This was a prospective study, which was approved by the local ethics committee (No.20200616017) and written informed consent was obtained from each participant.

We planned to enroll 155 consecutive recovered COVID-19 patients who were discharged from Fuyang No.2 People`s Hospital (the only designated hospital for infectious diseases in Fuyang, Anhui Province) to take MRI and clinical follow-up. The criteria of diagnosis and discharge of COVID-19 patients were based on the Polymerase Chain Reaction (PCR) result according to World Health Organization (WHO) guidelines [21] . The non-COVID-19 volunteers were recruited through social media. Each volunteer agreed to take the MRI scan and finished a questionnaire. Finally, 60 COVID-19 patients and 39 age-and sex-matched non-COVID-19 volunteers were successfully recruited to undertake three-dimensional T1-weighted imaging (3D-T1WI) and DTI scans. No artefacts were observed in their MR images.

The clinical characteristics of COVID-19 patients were obtained from the hospital records including age, gender, clinical type (mild/ severe/critical type, according to WHO guidelines [21] ), hospitalization days, existence of fever, cough, gastrointestinal symptoms, known contact history, alcohol intake history, smoking history and laboratory tests including white blood cell (WBC) count, lymphocyte count and lactate dehydrogenase (LDH). The laboratory tests were done at admission.

A semi-quantitative scoring system (CT extent score) was used to estimate the severity of pneumonia by assessing the CT scan at admission [22] . Each lung was divided into upper (above the tracheal carina), lower (below the inferior pulmonary vein) and middle (in between) zones, and each zone was scored based on the following criteria: 0, 0%; 1, < 25%;2, 25% -49%; 3, 50% -74%; 4, > 75%. The abnormal extent was determined by the summation of scores (possible range 0-24).

The detailed neurological symptoms both during infection and at this follow-up visit were collected via self-report, including the existence of headache, vision changes, hearing loss, loss of smell, loss of taste, impaired mobility, limb numbness, tremor, fatigue, myalgia, memory loss, mood changes. Handedness was also included in the questionnaire.

For non-COVID-19 volunteers, the basic information about age, gender, handedness, alcohol intake history, smoking history and underlying diseases were recorded.

Neuroimaging was performed using a 3.0 T Magnetom Skyra (Siemens, Germany). The structural 3 Dimensional T1-weighted Images (3D-T1WI) were acquired in sagittal plane using T1-weighted 3D Magnetization Prepared Rapid Gradient Echo (MPRAGE) sequence, with parameter setting: TR/TE=1900/2.84 ms, flip angle=5°, acquisition matrix=256 £ 256, FOV=280 £ 280 mm 2 , number of slices=160 and slice thickness=0.9mm.

DTI sequence was twice-refocused spin-echo sequence based on single-shot echo-planar acquisition. Diffusion sensitizing gradients were applied along 90 orthogonal directions using two b values (0 and 1000 s/mm 2 ) and other DTI parameters were: TR=3700 ms, TE=92 ms, FOV=220 £ 220 mm 2 , matrix=128 £ 128 and slice thickness=4.0 mm (voxel size 1.7 £ 1.7 £ 4.0 mm). The acquisition time per dataset was approximately 10 minutes. The original raw data were anonymized and transferred from the scanner in the DICOM format.

To calculate the quantitative parameters including gray matter volume (GMV), white matter volume (WMV), FA, MD, RD and AD values of different brain regions, we followed an atlas-based image processing approach (Figure 1 ). DTI datasets were first processed using DSI Studio software (http://dsi-studio.labsolver.org/) to generate FA, MD, RD and AD maps for each subject. Second, we performed a rigid [24] for quantitative parameters extraction in GM and WM, respectively. Finally, the brain was parcellated into 65 (AAL-3) plus 36 (JHU DTI-based WM atlas) anatomical regions. Brain regional volumes, and diffusion indices including FA, MD, RD and AD of WM and GM respectively can be extracted by averaging the values from voxels with partial volume exceeding 90% of the corresponding tissue. The estimated global GMV and WMV were further normalized by correction of intracranial volume.

We stratified our findings based on the following rules in order to further investigate the impact of the factors of interest on the global or regional volumetric or diffusive indices that were significantly different between groups. 1) Handedness: The data was compared between righthanded and lefthanded patients. 2) Neurological symptoms during acute stage: The data was compared between patients with and without neurological symptoms during acute stage. 3) Neurological symptoms at follow-up point: The data was compared between patients with and without neurological symptoms at the follow-up point. 4) Clinical severity. We allocated mild patients into non-severe subgroup, and severe subgroup included severe and critical ones and explored the difference.

We used a semi-quantitative scoring system by respectively comparing regional MD values of GM and WM, and FA values of WM between patients and controls: 0 not involved, MD GM or MD WM value in the patient is within 1 standard deviations (SD) of that in controls; 1 positive mild alteration, the value in the patient is 1 SD above but within 2 SD of that in controls; -1 negative mild alteration, the value in the patient is -1 SD below but within -2 SD of that in controls; 2 positive clear alteration, the value in the patient is 2 SD above that in controls; -2 negative clear alteration, the value in the patient is -2 SD below that in controls. The total index scores were calculated as a sum of the score (named MD-GM score) in each of the 65 brain regions according to AAL atlas, and in each of the 36 regions according to JHU atlas (named MD-WM and FA-WM scores). The MD-GM score ranged from -130 to 130 and the MD-WM and FA-WM scores ranged from -72 to 72 [25] .

The assessment of the correlation was performed among different index scores, regional volumes and diffusion indices with significant inter-group differences, and clinical characteristics.

All statistical analyses were performed with R software (version 3.5.3; R Foundation for Statistical Computing, Vienna, Austria). The Shapiro-Wilk test and Levene's test was used to evaluate the distribution type and homogeneity of variance. Volumetric and diffusion indices were compared between COVID-19 and control groups using analysis of covariance (ANCOVA) to adjust the effects of age, gender and whole brain volume. To control for multiple comparisons in brain regional analysis, false discovery rate (FDR) method was applied to correct the p value. Two sample t-test and Chi-squared test/ Fisher exact test were performed to compare the continuous and categorical measurements between two groups. Spearman Correlation coefficients were calculated to explore relationship between volumetric/ diffusion indices and clinical data. A p value of <0.05 was defined as statistical significance.

No study sponsor had any role in study design; collection, analysis, or interpretation of data; or in writing this paper or the decision to submit for publication. The corresponding author had full access to all study data and had final responsibility for the decision to submit for publication.

We recruited 60 COVID-19 patients (mean age § SD, 44. 10 (15, 25 .00%), vision change (13, 21.67%), myalgia (9, 15.00%), impaired mobility (7, 11.67%), memory loss (8, 13.33%), taste loss (4, 6.67%), limb numbness (4, 6.67%), tremor (4, 6.67%), smell loss (2, 3.33%) and hearing loss(1, 1.67%). Other common symptoms included fever (53/60, 88.33%), cough (34/60, 56.67%), and gastrointestinal discomfort (8/60, 13.33%). At this follow-up visit, 33/60 (55%) patients still had neurological symptoms. The numbers of patients feeling fatigue and mood change declined significantly compared with before (p value < 0.05). 47 (78.33%) of the patients were classified as the mild type of COVID-19, 12 (20.00%) patients were severe type, and 1 (1.67%) was grouped as critical type. (Table 1) 

The comparisons of the mean regional GMVs between the COVID-19 group and the control group were reported in Table 2 and those of WMVs were shown in Table 3 . VBM analysis showed the significantly higher GMV in the left Rolandic operculum area (p value = 0.019), bilateral olfactory cortices (p value = 0.024 and 0.043 for left and right), bilateral insulas (p value = 0.013 and 0.013), bilateral hippocampi (p value < 0.001, and = 0.013), right cingulate gyrus (p value = 0.118) and left Heschl's gyrus (p value < 0.001) in COVID-19 patients. No significant differences were seen in other regional GMVs or any regional WMVs (all p values > 0.05).

The comparisons of diffusion indices were performed. Differences in regional MD values of gray matter between two groups were collected in Table 4 , and regional FA, MD, AD and RD values of white matter were illustrated in Table 5 . In the gray matter of COVID-19 patients, the left insula (p value <0.001), bilateral cingulate gyri (p value = 0.039 and 0.305), right precuneus (p value = 0.279), and right thalamus (p value = 0.033) were found to have significantly lower MD values when compared with the control group. In white matter, the mean regional MD, AD and RD values were generally lower in the COVID-19 group, and the mean regional FA values were generally higher compared with the controls. The mean MD values of right 

The differences in global GMV along with diffusion indices of white matter, the AD, RD, MD, and the FA were summarized in Table 6 . Significantly higher mean global GMV (F = 143, p value= 0.02), higher global FA value (F = 100, p value = 0.021), lower global MD (F = 80, p value= 0.001), AD (F = 87, p value= 0.002) and RD (F = 72, p value= 0.002) were found in the COVID-19 group than the control group.

The regions with significant differences in the volumes and diffusion indices between these two groups were shown in Figure 2 .

The results were collected in Supplementary Table 1 . No statistically significant differences were found in the indices of interest between patients that were righthanded (n = 54) and lefthanded (n=6), patients with (n = 41) and without (n = 19) neurological symptoms during acute stage, patients with (n = 33) and without (n = 27) neurological symptoms at the follow-up point, or patients classified as non-severe (n = 47) and severe (n = 13) (all p values> 0.05). Additionally, we were interested in the changes of regions in the olfactory cortex-related regions in patients with smell loss, thus related regional GMVs and MD of GM were compared (Supplementary Table  2 ). No significant differences were found in either comparison (all p values > 0.05), yet the average regional GMVs and MD of GM of patients with smell loss (n = 2) were generally lower than patients without the symptom (n = 58). The assessment of the correlation was performed among different index scores, regional volumes (including global GMV and 9 regional GMVs) and diffusion indices (including 4 global indices of WM, 1 regional MD of WM, 5 regional MD of GM, and 3 regional AD of WM) with significant inter-group differences, and clinical characteristics.

The significant correlations were shown in Figure 3 and Supplementary Table 3 between neurological symptoms during the acute stage and regional WMVs, diffusion indices, and index scores. MD values of GM in right thalamus and AD values of WM in right EC were positively correlated with vision changes (r = 0.336 and 0.285 respectively, p value = 0.027 and 0.046 respectively). GMVs in right cingulate and hippocampus were negatively correlated with smell loss (r = -0.284 and -0.279 respectively, p value = 0.028 and 0.031 respectively).

The correlations were also shown in Figure 3 and Supplementary  Table 3 . FA-WM score was negatively related to tremor (r = -0.262, p value = 0.043), and MD-WM score was positively correlated (r = 0.279, p value =0.031). Global GMV and regional GMV in left Rolandic operculum, right cingulate, bilateral hippocampus, and the left Heschl's gyrus, as well as global MD value of WM were all negatively correlated with memory loss (r = -0.341, -0.369, -0.394, -0.321, -0.277, -0.279 and -0.285 respectively, p value = 0.008, 0.002, 0.009, 0.012, 0.032, 0.031, 0.027). MD values of GM in cingulate gyri were correlated with fatigue (r = 0.305 and 0.320, p value =0.018 and 0.013 respectively for the left and the right) and numbness in the extremities (r = 0.286, p value = 0.027), and MD of GM in right precuneus correlated to numbness (r = 0.309, p value = 0.016). Figure 3 displayed all the correlation (p value < 0.01).

Other important clinical data included handedness, typing of COVID-19, CT extent score, etc. MDs of GM in bilateral cingulate gyri were found to correlate with clinical type (p value = 0.035 and 0.036 for left and right side). Other significant correlations can be found in Figure 4 and Supplementary Table 4 .

MD-GM score (r = -0.317, p value = 0.014), global GMV (r = -0.282, p value = 0.029), and GMV in right cingulate gyrus (r =0.160, p value = 0.041) correlated with LDH. A variety of indices were found to correlate with WBC count and lymphocyte percentage. (See Figure 4 and Supplementary Table 4 ). A plot with regression was produced for the correlation between LDH and global GMV ( Figure 5 ).

In this prospective work, we recorded the detailed cerebral volumetric, DTI metrics 3 months after SARS-CoV-2 infection by applying DTI and 3D T1WI in 60 recovered COVID-19 patients and 39 age-and sex-matched normal volunteers. Overall, these recovered COVID-19 patients were more likely to have enlarged olfactory cortices, hippocampi, insulas, Heschl's gyrus, Rolandic operculum and cingulate gyrus, and a general decline of MD, AD, RD accompanied with an increase of FA in white matter, especially AD in the right CR, EC and SFF, and MD in SFF compared with non-COVID-19 volunteers. Global GMV, GMVs in left Rolandic operculum, right cingulate, bilateral hippocampus, left Heschl's gyrus, and Global MD of WM were found to correlate with memory loss. GMVs in right cingulate gyrus and left hippocampus were related to smell loss. MD-GM score, global GMV, and GMV in right cingulate gyrus were correlated with LDH level.

Given the olfactory and gustatory dysfunction in COVID-19 patients and evidence of olfactory epithelium invasion by SARS-CoV, the olfactory gyrus was thought to be the first functional area in CNS to be infected by SARS-CoV-2 [26] . To avoid cross-infection, patients were unable to undertake MRI scans during the acute phase, but we are still curious about whether any micro-structural changes exist in the recovery stage of COVID-19 and whether any clues were left to suggest the probable intracranial infection route.

According to our results, significant enlarged volumes were observed in the bilateral olfactory cortices, hippocampi, insulas, left Heschl's gyrus, left Rolandic operculum and right cingulate gyrus. All these structures mentioned above belonged to the central olfactory system. Among them, olfactory cortex, also named as piriform cortex, directly receives axonal projections from olfactory bulb (OB), is referred to as a part of 'primary olfactory cortex'. Other structures are the cortical targets of primary olfactory cortex in the bilateral limbic lobe, temporal cortices, which were referred as 'secondary olfactory cortex' [27] . It was reported that the frequent olfactory loss during the course of upper respiratory tract infections (URTI) resulted in loss of stimulation and subsequent volume loss in the acute phase, while after olfactory recovery, the volumes of GM in the central olfactory system got enlarged subsequently [28, 29] . We analyzed the GMV difference between patients with or without olfactory loss and found the GMVs of the central olfactory system were generally smaller in patients with persistent olfactory loss compared with those without olfactory problems which was aligned with previous similar studies (Supplementary Table 2 ) [28] .

Several possible invasion routes of SARS-CoV-2 were raised including hematogenous, lymphatic and neuro retrograde routes, etc., yet the exact route was unknown. Based on our findings, the gray matter volumetric changes in the central olfactory system provided a speculation that SARS-CoV-2 might enter the CNS via an OBmediated neuronal retrograde route. Two reasons were suspected to play a role in these GMV enlargement: the neurogenesis and functional compensation. Firstly, patients were suspected to experience neurogenesis. It is well accepted that neurogenesis in adults are restricted to two regions, the subventricular zone (SVZ) and the sub-granular layer of the dentate gyrus of the hippocampus [30] . The neuroblasts from SVZ migrate along the rostral migratory stream, enter olfactory cortex first and finally replace interneurons (e.g., periglomerular cells, granular cells) in the OB [30] . Therefore, the increasing neurons possibly resulted in the enlargement of the GMV in the olfactory system. Secondly, in order to compensate for impaired olfaction, the increased functional engagement of brain regions would be hypertrophy, which was proved to have enlarged neurons and increasing number of dendritic spines by experimental studies in sensory deprivation models [28] . These two reasons might be an explanation of the enlarged GMV, while further pathological exploration was needed.

In general, lower diffusivity parameters (MD, AD, RD) and higher FA values were recognized in the white matter from COVID-19 cohort. In addition, the diffusivity (MD and AD) of right CR, EC and SFF decreased significantly. The CR, consisting massive bundle of projection fibers, connects the cortex to the brainstem and the thalamus in afferent and efferent manners [31] . The EC and SFF are series of association fibers, connecting frontal, parietal, and temporal cortices. White matter was not the main target for neurotropic virus; however, the connecting fibers could act as the channel for intracranial viral transmission. Different from the elevated MD with compressed volume in fibers in hydrocephalus, the enlargement of white matter fibers, decreased MD values and elevated FA suggested a greater alignment of fibers and limited diffusion freedom, indicating a possible intrinsic re-construction (e.g. remyelination) process that occurred after infection [32] .

It was interesting to notice that all the diffusivity abnormalities in the white matter restricted in the right hemisphere, without asymmetrical symptoms reported by COVID-19 patients. Since the blood flow rate and volume was higher in the hemisphere dominant for handedness, we hypothesized that the predominance of abnormal white matter diffusivities might be related to the difference of blood volume in bilateral hemispheres [33] . However, after comparing these abnormal diffusivity indices between lefthanded and righthanded patients, no difference was found (Supplementary Table 1 ). The asymmetrical phenomenon was detected in other studies as well [28, 34] . The right-side predominance in odorant perception has been indicated by multiple olfactory functional studies which was not fully understood. The diffusivity changes on the right side might be related to the right-side odorant perception in which further exploration was needed [28] .

During SARS-CoV-2 infection, 41/60 (68.33%) patients had neurological symptoms and over 50% recovered patients still had symptoms 3 months later. It was interesting to find the GMV in hippocampi (a key part in the organization of memory) and cingulate gyri (an important part of limbic system) were negatively related to loss of smell during infection and loss of memory 3 month later, which could support our hypothesis of neurogenesis in these regions mentioned above. Tremor was found to be negatively related to FA_WM score both in the acute stage and at the 3-month follow-up point which indicated a destruction of WM fibers in both hemispheres, possibly resulted from SARS-CoV-2 induced cytokine storm [35] . According to WHO guideline, COVID-19 patients were divided into mild, severe and critical types based on their clinical information and laboratory results [21] . No significant difference was observed between severe and non-severe groups. But the clinical types were positively related to MD values in bilateral cingulate gyri (r=0.271 and 0.272, p=0.035 and 0.036 respectively), implying the more severe the case was, the higher the MD value of bilateral cingulate gyri was presented. Cingulate gyrus usually plays an important role in attention, motivation, decision making, learning, and cost-benefit calculation, as well as conflict and error monitoring, which is frequently affected in limbic encephalitis [36] . Therefore, the dysregulated cytokine response was speculated in severe cases [35] .

After exploring the relationship between laboratory data and DTI metrics, the global GMV was significantly but slightly correlated with the LDH concentration in COVID-19 patients. LDH is one of the key enzymes in the glycolytic pathway, highly expressed in cells from kidney, heart, liver and brain [37] . Elevated concentrations of LDH are observed in patients with encephalitis, ischemic stroke and head injuries [37] . Higher concentration of serum LDH always follows tissue breakdown and is closely linked to the deterioration and poor outcome [38] . The decreased global GMV in LDH-elevated patients might indicate an atrophy due to a severe inflammatory response.

Since the hemostatic abnormalities, including disseminated intravascular coagulation (DIC) and severe inflammatory response, were frequently observed in COVID-19 patients, individuals may predispose to cerebral-vascular events caused by infection and treatment [39] . An index score system was introduced into our study to investigate the existence of micro-structural abnormalities due to microvessel diseases. The ischemic changes are known to be accompanied with lower FA value and higher MD value in the ischemic lesions. However, our findings showed that FA-WM score was higher, MD-GM and MD-WM scores were lower in the COVID-19 group compared with the control group, and the differences were insignificant. Thus, no obvious evidence of micro-vessel diseases was found.

It is important to investigate the relationship between abnormal anatomical brain areas and ACE-2 distribution. It is clarified that SARS-CoV-2 enters the host cell by attaching with ACE-2 via Spike (S) glycoprotein. Therefore, the more expression of ACE-2 might bring more severe abnormalities. The distribution of ACE2 was non-equivalent over the brain and was most frequently expressed in substantia nigra, followed by spinal cord, hippocampus, basal ganglia, limbic system and frontal cortex [17] . Our results suggested that various components in the limbic system were affected structures sharing possible high ACE-2 expression, which were partly aligned with the proposed ACE-2-riched regions. Although we were not able to observe the DTI metrics in substantia nigra since it was not included in the brain atlas we used, it was still very hard to support any relationship between ACE2 expression and affected brain areas.

At the time of writing, there was still no research to detect the cerebral micro-structural changes after SARS-CoV-2 infection from imaging or pathological aspects. Our study gave a hint to possible neurological changes after SARS-CoV-2 infection. The limitations of our study were listed as follows: 1) we did not enroll enough patients with neurological dysfunction or olfactory loss, therefore the relationship between GMV/diffusivity changes and olfactory symptoms would be missed; 2) as a single-centered study, a selection bias might result from limited ethnical and regional characteristics of the participants, and possible mutants of SARS-CoV-2 in other countries, and limit the generalization of the study; 3) the atlas we applied did not contain the structure in brainstem, therefore, we failed to obtain the volumetric and DTI information about the nucleus in brainstem, some of which were quite important, especially the solitary nuclei.

In this prospective study, volumetric and micro-structural abnormalities were detected mainly in the central olfactory cortices, partial white matter in the right hemisphere from recovered COVID-19 patients, providing new evidence to the neurological damage of SARS-CoV-2. The abnormalities in these brain areas might cause long-term burden to COVID-19 patients after recovery, which was thus worth public attention.

Yiping Lu, Bo Yin and Daoying Geng conceived and designed the research. Yiping Lu, Anling Xiao, Xuanxuan Li, Nan Mei, Yajing Zhao and Dongdong Wang contributed to data collection and interpretation. Yiping Lu, Xuanxuan Li, Pu-Yeh Wu, Chu-Chung Huang and Tianye Jia contributed to statistical analysis. Yiping Lu, Xuanxuan Li and Nan Mei contributed to draft and revise the manuscript. Anling Xiao and Bo Yin supervised the whole research. Yiping Lu, Xuanxuan Li, Daoying Geng and Nan Mei contributed equally to this study. All authors gave final approval for the version to be published.

According to the informed consent, all the data involved is only used in this research and our institution. Thus, supporting data is not available for the public due to the ethical restriction.

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. This paper has never been published elsewhere.

A new coronavirus associated with human respiratory disease in China

The SARS-CoV-2 exerts a distinctive strategy for interacting with the ACE2 human receptor

Human coronaviruses: Viral and cellular factors involved in neuroinvasiveness and neuropathogenesis

Tissue distribution of ACE2 protein, the functional receptor for SARS coronavirus. A first step in understanding SARS pathogenesis

Detection of Severe Acute Respiratory Syndrome Coronavirus in the Brain: Potential Role of the Chemokine Mig in Pathogenesis

Neurologic Manifestations of Hospitalized Patients With Coronavirus Disease

Updates on What ACS Reported: Emerging Evidences of COVID-19 with Nervous System Involvement

Cerebral Venous Thrombosis associated with COVID-19 infection: causality or coincidence?

A first case of meningitis/encephalitis associated with SARS-Coronavirus-2

Pathological findings of COVID-19 associated with acute respiratory distress syndrome

The pathogenesis of murine coronavirus infection of the central nervous system

Cell senescence is an antiviral defense mechanism

Does SARS-CoV-2 infection cause chronic neurological complications? GeroScience 2020

A review of pathophysiology and neuropsychiatric manifestations of COVID-19

The spatial and cell-type distribution of SARS-CoV-2 receptor ACE2 in human and mouse brain

Statistical Instability of TBSS Analysis Based on DTI Fitting Algorithm

Brain alterations within the first 100 days of HIV infection

Herpes Simplex Virus-1 Encephalitis in Adults: Pathophysiology, Diagnosis, and Management

Clinical management of severe acute respiratory infection when COVID-19 is suspected (v1.2). Who; 2020

Development and Validation of a Prognostic Nomogram Based on Clinical and CT Features for Adverse Outcome Prediction in Patients with COVID-19

Automated anatomical labelling atlas 3

Tract probability maps in stereotaxic spaces: Analyses of white matter anatomy and tract-specific quantification

Diffusionweighted MRI signal abnormality in sporadic CJD increases in extent and intensity with disease duration

Severe Acute Respiratory Syndrome Coronavirus Infection Causes Neuronal Death in the Absence of Encephalitis in Mice Transgenic for Human ACE2

Magnetic Resonance Imaging of Human Olfactory Dysfunction

Congenital olfactory impairment is linked to cortical changes in prefrontal and limbic brain regions

Brain volume changes in hyposmic patients before and after olfactory training

Human neuroblasts migrate to the olfactory bulb via a lateral ventricular extension

Structure of corona radiata and tapetum fibers in ventricular surgery

Axial diffusivity of the corona radiata correlated with ventricular size in adult hydrocephalus

Hand preference is selectively related to common and internal carotid arterial asymmetry

Anosmia leads to a loss of gray matter in cortical brain areas

COVID-19: consider cytokine storm syndromes and immunosuppression

Anterior cingulate cortex involvement in non-paraneoplastic limbic encephalitis

The Regulation and Function of Lactate Dehydrogenase A: Therapeutic Potential in Brain Tumor

Lactate dehydrogenase 5: An old friend and a new hope in the war on cancer

COVID-19 and Thrombotic or Thromboembolic Disease: Implications for Prevention, Antithrombotic Therapy, and Follow-up

We thank Wei Cheng (Ph.D., Institute of Science and Technology for Brain-Inspired Intelligence, Fudan University) for providing statistical consultation and technological support. 

Supplementary material associated with this article can be found in the online version at doi:10.1016/j.eclinm.2020.100484.