Indoor incense burning impacts cognitive functions and brain functional connectivity in community older adults

To investigate (1) the effects of indoor incense burning upon cognition over 3 years; (2) the associations between indoor incense burning with the brain’s structure and functional connectivity of the default mode network (DMN); and (3) the interactions between indoor incense burning and vascular disease markers upon cognitive functions. Community older adults without stroke or dementia were recruited (n = 515). Indoor incense use was self-reported as having burnt incense at home ≥ weekly basis over the past 5 years. Detailed neuropsychological battery was administered at baseline (n = 227) and the Montreal Cognitive Assessment at baseline and year 3 (n = 515). MRI structural measures and functional connectivity of the DMN were recorded at baseline. Demographic and vascular risk factors and levels of outdoor pollutants were treated as covariates. Indoor incense burning was associated with reduced performance across multiple cognitive domains at baseline and year 3 as well as decreased connectivity in the DMN. It interacted with diabetes mellitus, hyperlipidemia and white matter hyperintensities to predict poorer cognitive performance. Indoor incense burning is (1) associated with poorer cognitive performance over 3 years; (2) related to decreased brain connectivity; and (3) it interacts with vascular disease to predispose poor cognitive performance.

The objectives of this study are to investigate, in older adults free of stroke and dementia, 1) the effects of household indoor incense burning upon cognition over 3 years; 2) the effects of indoor incense burning upon structural MRI markers of SVD (WMH and lacune) and medial temporal lobe atrophy (MTLA), a proxy of AD, as well as functional connectivity of the default mode network (DMN); and 3) whether indoor incense burning moderates the relationships between vascular risk factors and structural imaging markers with cognitive functions.

Methods
Study design. Prospective, 3-year longitudinal case-controlled study.

University of Hong Kong RISK Index for Screening Subclinical Brain Lesions in Community-dwelling in Hong
Kong (CU-RISK Study). Inclusion criteria were (1) age ≥ 65 years; (2) community-dwelling and (3) written informed consent given. Exclusion criteria were (1) dementia at baseline, as defined by the locally validated education-adjusted cut-off scores on the Mini-mental State Examination 16 or a clinical diagnosis of dementia; (2) history of stroke; (3) inadequately controlled psychiatric disorders; and (4) physical or sensory impediments hindering participation in cognitive testing. Baseline data collection took place between November 2011 and March 2016. A subset of 515 randomly selected participants received follow-up between October 2014 and December 2017. Flow diagram describing the sample recruitment and cognitive data collection is presented as Fig. 2. The CU-RISK study complies with the relevant guidelines in the Declaration of Helsinki and approval was granted from the Chinese University of Hong Kong -New Territories East Cluster Clinical Research Ethics Committee. This study is part of the Chinese University of Hong Kong "Brain Health Brings Health" programme.
Recording of incense use habits and potential confounders. Incense use habit was recorded based on self-report as current regular practise of incense burning at home on a weekly or more frequent basis over the past 5 years. An apartment of total area between 200 to 800 square feet consisting of a living room and 1-3 bedrooms represents a typical indoor household environment in Hong Kong. www.nature.com/scientificreports www.nature.com/scientificreports/ Outdoor air pollutant levels. The general regional air quality data (AQD) published by Department of Environment Protection of the Government of Hong Kong Special Administrative Region was collected from a publicly accessible website (http://epic.epd.gov.hk/EPICDI/air/station/?lang=en) to measure the levels of 6 outdoor air pollutants: fine suspended particulates, NO 2 , O 3 , SO 2 , nitrogen oxides (NO) and respirable suspended particulates measured hourly over 13 districts of Hong Kong. The yearly averaged readings of the different pollutants measured in the station matched to the residential address of the individual participants were recorded.

Cognitive assessment. Cognitive assessment was conducted by trained research assistants. The Montreal
Cognitive Assessment (Hong Kong version; MoCA) 17 covers cognitive domains including learning and memory, executive and visuospatial functions, language, attention, working memory, abstraction and orientation and is shown to be sensitive to mild cognitive dysfunction and cognitive change over time 18,19 . The total score is used as an index of global cognitive functions. In addition to the total score, domain scores for memory, executive functions, visuospatial functions, language and attention were calculated from a previously published method 20 . A subset of 274 participants also underwent detailed neuropsychological testing covering memory, executive functions/activation, language and visuospatial functions using the National Institute of Neurological Disorders and Stroke and Canadian Stroke Network VCI Harmonization 60-minute neuropsychology protocol at baseline 21,22 . Only the MoCA was administered at year 3 ( Fig. 2).
Data Pre-processing. MRI markers including WMH volume, lacune and MTLA were measured. The WMH was defined as hyperintensity as reflected on FLAIR image, according to the STRIVE standard 23 and automatically segmented with manual correction where appropriate. WMH segmentation was performed based on a validated pipeline -coarse-to-fine detection of WMH using co-registered T1W and FLAIR 24 . Lacunae were defined according to the STRIVE standard and segmented in T1W images 23 . The number was counted as isolated regions automatically. MTLA was visually rated on T1W images using the Schelten's scale 25 . fMRI imaging data were pre-processed using Statistical Parametric Mapping software version 12 (http://www.fil.ion.ucl.ac.uk/spm/). The first 10 volumes of each participant were discarded to allow for T1W equilibration effects. Then the functional images were slice-time corrected for timing offsets between different slices and realigned to the first image to correct for head motion between scans. The high-resolution T1W image was then co-registered to the mean of the corrected functional images. In order to improve the brain tissues segmentation in aging, a multi-channel segmentation approach in SPM 12 was used. The FLAIR image was coregistered to the T1W image, and the coregistered T1W and FLAIR images were combined underwent multi-channel segmentation to extract different tissues including gray matter, white matter and cerebrospinal fluid. A study-specific template was created using the DARTEL toolbox 26 . The functional images were then spatially normalized to the standard MNI space by using the nonlinear normalization parameters estimated by the DARTEL toolbox, resampled to 3 × 3 × 3 mm 3 and spatially smoothed with a 6 mm full-width half-maximum Gaussian kernel. Finally, additional preprocessing steps were implemented to the normalized function images to eliminate the effect of low-frequency drifts and physiological noise, which included removing linear trends, temporally band-pass filtering (0.01-0.1 Hz) and regressing out several nuisance signals (six head-motion profiles, the averaged signals from white matter, cerebrospinal fluid, and the whole brain and the first derivatives thereof) 27 .
Functional connectivity analysis in the DMN. The DMN was identified using a seed-based correlation approach. First, a 5-mm-radius sphere seed centred on the posterior cingulate cortex (MNI coordinates: -7, -43, 33) was selected 28 . The mean BOLD time series of the seed was then extracted and correlated with every voxel of the whole brain using Pearson correlation. The correlation coefficients were then converted to z-scores using Fisher's r-to-z transformation. A one-sample t-test was used to identify the DMN for each group, the p-value threshold was set to p < 0.05 (false discovery rate correction (FDR) for multiple comparisons, voxel size > 30). Then, analysis of covariance (ANCOVA) was used to evaluate the group differences between the incense users and non-users in the DMN with age, sex and years of education as covariates. The significance threshold was set to p < 0.05 (FDR corrected, voxel size > 30).
Statistical analysis. Statistical analysis was performed using IBM SPSS version 21.0. Statistical significance was determined at p < 0.05 for all analyses. Participants were categorized into those who reported having regularly burned incense (users) and who did not (non-users). Group comparisons on demographic and clinical data were conducted using independent sample t tests or Chi-squared tests as appropriate. ANCOVA was used to compare the cognitive performance between the groups at baseline and year 3. The group effect size on cognitive outcomes was measured by the Cohen's d statistic, with 0.2, 0.5, and ≥0.8 indicating small, medium and large effects, respectively 29 . Furthermore, a repeated measure ANCOVA was used and a group x time interaction term was calculated to test whether incense users and non-users experienced a different cognitive trajectory over 3 years as measured by the MoCA Total score. Linear regression models were constructed to examine the associations between incense burning and cognitive functions at baseline and year 3 and the change in MoCA Total score over 3 years. Variables selected as covariates were demographic factors (age and years of education), vascular risk factors (hypertension, diabetes mellitus [DM] and hyperlipidemia, current smoking, current alcohol drinking) and 6 types of outdoor air pollutants levels. Table 1 enlists the cognitive and imaging markers used as outcomes in this study.
Continuous variables were examined first for collinearity. As significant collinearity was found between fine suspended particulates and respirable suspended particulates and between NO 2 and nitrogen oxides, only the results based on fine suspended particulates and NO 2 are presented. Two-way ANCOVA was conducted to examine the interactions between incense burning with vascular risk factors or structural imaging markers upon MoCA Total score while controlling for demographic and other vascular risk factors. For imaging markers, participants were classified into low and high WMH levels by median split at 3.74mm 3 based on data of the whole sample. MTLA was determined as ratings ≥2 on the Schelten's scale 25 .
In view of the group difference in age and year of education in the whole sample, propensity score was used to match incense users and non-users with age and years of education as matching covariates for functional imaging data analysis and as sensitivity analysis for structural imaging data. Each incense user was matched to one non-user without replacement. The distributions of matching covariates between groups were compared with the matched sample to ensure a balanced matching. Inverse probability of treatment weighting (IPTW) using the propensity score was used to create a synthetic sample such that the baseline covariates were independent of grouping.

Ethics approval and consent to participate. The Chinese University of Hong Kong -New Territories
East Cluster Clinical Research. Ethics Committee approved the CU-RISK study. Informed consent was obtained from the study participants.

Results
Five hundred and fifteen participants were recruited. Compared to non-users, users were 1.1 years older, had 2.7 years less education and were less likely to have hypertension. (Table 2). Exposure to outdoor pollutants was similar between groups. The 227 participants with detailed neuropsychological testing were similar in age and sex compared to those without. Participants with detailed testing had slightly lower MoCA Total score (mean difference −1.42 [3.81]) and less education (mean difference in years −1.84 [4.81]) compared to those without detailed testing. The interval between baseline and follow-up was 38.3 (4.1) months. Descriptive statistics and effect size estimates (where applicable) are presented in Table 2.  Table 2).
Year 3. The MoCA was administered at year 3, significant group difference persisted for MoCA Total score, MoCA Visuospatial Domain score and MoCA Memory Domain score ( Table 2). The group x time interaction (2020) 10:7090 | https://doi.org/10.1038/s41598-020-63568-6 www.nature.com/scientificreports www.nature.com/scientificreports/ was not statistically significant (p = 0.732), indicating that the trajectory of change in MoCA Total score between baseline and year 3 did not differ between incense users and non-users.  Table 3 shows the results of the linear regression models. Incense burning was not associated with change in MoCA Total score over 3 years (β = 0.021, p = 0.651).  www.nature.com/scientificreports www.nature.com/scientificreports/ Neuroimaging outcomes. Structural MRI. No group difference was found in WMH volume, lacune count and MTLA. Sensitivity analysis conducted in the propensity matched age-and education-matched sample (n = 200, incense users n = 104, non-users n = 96) yielded essentially the same results (Table 2).

Associations between incense burning and cognitive outcomes at baseline and year 3.
Resting-state fMRI. Functional connectivity of the DMN of the non-users and users are shown in Fig. 3a.b, respectively. Significant DMN clusters of both groups were mainly identified in the superior/medial frontal gyrus, precuneus, and middle temporal gyrus. Incense users had significantly decreased connectivity in the precuneus, medial frontal gyrus, left angular and right middle temporal gyrus (Fig. 3c). Table 2 shows the numerical group comparison of the averaged functional connectivity values in the brain regions with significant group differences.
Interactions between incense burning with vascular risk factors and structural imaging markers upon cognitive performance. Significant interactions were observed between incense burning with DM (p = 0.031), hyperlipidemia (p = 0.036) and WMH volume (p = 0.008) on baseline MoCA Total score, indicating that incense burning had significantly more negative impact upon global cognition in the presence of DM, hyperlipidemia and increasing WMH volume. No interaction was found between incense burning with hypertension, the presence of lacune or significant MTLA on the MoCA Total score (Fig. 4).

Discussion
This study shows that regular indoor incense burning is associated with poorer performance across multiple cognitive domains over 3 years. However, indoor incense burning was not associated with a more rapid decline in cognitive functions as measured by the MoCA Total score. Furthermore, although participants who practised indoor incense burning did not have more structural brain changes in terms of SVD lesions measured as WMH or lacune or MTLA, which is an imaging biomarker for AD, a subtle impact of incense burning upon the brain was evident on the level of functional connectivity of the DMN. The DMN is active at rest and during introspective, self-referential processing. The proper switching between the DMN and active brain networks is responsible for a variety of cognitive processes and its dysregulation is associated with a number of neurological and psychiatric conditions including attention deficit hyperactivity disorder 30 , autism spectrum disorder 31 , schizophrenia 32 , mild cognitive impairment 33 , AD as well as the progression of mild cognitive impairment into dementia 34 . This study suggests that indoor incense burning induces functional changes in the brain that may reduce cognitive resilience through functional connectivity alternation and thereby increases vulnerability for future cognitive decline, a hypothesis to be tested in future longitudinal studies. Moreover, although incense users did not have a higher frequency of DM, hyperlipidemia or WMH, incense burning appeared to interact with these vascular diseases to predispose poor cognitive functioning, highlighting the potential role of incense burning as a risk Figure 3. Results of fMRI seed-based functional connectivity in the DMN. One-sample t-test was used to identify the DMN with p-value threshold set to p < 0.05 with FDR correction (voxel size >30) in non-incense users (a) and incense users (b). Two-sample t-test was used to evaluate the group differences. Compared to nonusers, self-report incense users showed significantly decreased connectivity in the Precuneus, Medial Frontal Gyrus, left Angular and right Middle Temporal Gyrus. The blue indicates decreased connectivity (p < 0.05) in the incense group (c). www.nature.com/scientificreports www.nature.com/scientificreports/ factor for VCI. Given the high prevalence of vascular burden among older adults, indoor incense burning should be practised with caution in this vulnerable population.
The precise mechanisms underlying the relationships between indoor incense burning with cognitive and brain changes need further investigations. Previous studies have reported that a wide range of air pollutants including ultrafine particles, PM, and volatile organic compounds produced from incense burning are detrimental to the vascular health 35 . Studies have demonstrated the capacity of incense smoke to induce oxidative stress and inflammation [36][37][38][39] , which may lead to significantly reduced vascular nitric oxide levels and increased levels of endothelin-1 and inflammatory mediators, such as granulocyte-macrophage-colony stimulating factor. The unregulated levels of these mediators would destruct the vascular vasodilation and constriction 40,41 , thereby promoting vascular dysfunction. It was suggested the water-soluble factions in particulate matter such as the transition metals with redox potential can also accelerate the process in membrane lipid peroxidation, which trigger endothelial cell mutations 42 . Moreover, the greater surface-to-mass ratio of ultrafine particles and PM can further increase the contacts with cells and enhance the damages, which may eventually lead to hemorrhagic events in blood vessels 43 .
To our knowledge, this is the first study to demonstrate the adverse effects of indoor household incense burning upon cognitive functions and brain health. The strengths of this study include a relatively large and well-defined community sample with a combination of cognitive, structural and functional imaging data and longitudinal follow-up. Also, the influence of regionally measured levels of a variety of outdoor air pollutants was taken into account as an important confounder. This study has a number of limitations. First, as incense use habit was evaluated from self-report, its accuracy might have been subjected to reporting bias. Second, although the common practise is to burn incense one to three times a day, each time using one to three incense sticks, details of incense use habits and home environment pertinent to the amount of exposure, such as the frequency of incense burning, the type and number of incense sticks used, indoor airflow variables, the specific location at home where the incense was burned and concomitant sources of indoor air pollutions, were not recorded. Third, past exposure to incense smoke was not assessed. However, it is safe to assume that these older participants have practised incense burning as religious rituals for many decades. Fourth, without real-time indoor air quality monitoring, the contributory roles of the various pollutants contained in incense smoke upon cognitive and neuroimaging outcomes could not be measured precisely. Last but not least, detailed neuropsychology and neuroimaging were not repeated at longitudinal follow up. Despite these limitations, this study has identified indoor incense burning as a novel and easily modifiable risk factor for adverse cognitive and brain health.