Dose-dependent genotype effects of BDNF Val66Met polymorphism on default mode network in early stage Alzheimer's disease

In humans, brain-derived neurotrophic factor (BDNF) has been shown to play a pivotal role in neurocognition, and its gene contains a functional polymorphism (Val66Met) that may explain individual differences in brain volume and memory-related activity. In this study, we enrolled 186 Alzheimer's disease (AD) patients who underwent 3D T1 magnetic resonance imaging, and explored the gray matter (GM) structural covariance networks (SCN). The patients were divided into three groups according to their genotype: Met/Met (n = 45), Val/Met (n = 86) and Val/Val (n = 55). Seed-based analysis was performed focusing on four SCN networks. Neurobehavioral scores served as the major outcome factor. Only peak cluster volumes of default mode medial temporal lobe network showed significant genotype interactions, of which the interconnected peak clusters showed dose-dependent genotype effects. There were also significant correlations between the cognitive test scores and interconnected-cluster volumes, especially in the orbitofrontal cortex. These findings support the hypothesis that BDNF rs6265 polymorphisms modulate entorhinal cortex-interconnected clusters and the valine allele was associated with stronger structural covariance patterns that determined the cognitive outcomes.


IntroductIon
Neurotrophins are a group of molecules which play a key role in regulating neural survival, development and maintenance. The most abundant neurotrophic factor is brain-derived neurotrophic factor (BDNF) [1,2] and has been reported to be of clinical significance in mediating the hippocampus network function [3,4]. Activity-dependent BDNF secretion has been reported to be required for long-term potentiation and depression [5,6] while impaired patterns of discrimination and learning deficits have been observed in BDNF knockout mice [7]. The salient pathological feature in Alzheimer's disease (AD) is hippocampal network degeneration, and BDNF has been shown to play a protective role in attenuating amyloid-related toxicity [5,8].
In humans, the BDNF gene has been mapped to chromosome 11p14. 1

. A common single nucleotide
Research Paper: Gerotarget (Focus on Aging) www.impactjournals.com/oncotarget polymorphisms (SNP) consisting of a missense change of the coding exon at position 66 (Val66Met,rs6265) has been shown to produce non-conservative amino acid changes (Val to Met). The Met allele carriers have been shown to have poorer episodic memory, and this may be reflected in lower activities of BDNF with impaired neuronal processing and trafficking [9]. A recent metaanalysis study [10] showed no significant associations between the BDNF Val66Met polymorphism and the risk of developing AD in dominant (Met vs. Val/Val), recessive (Met/Met vs. Val) and additive (Met/Met vs. Val/Val) carriers. Another meta-analysis reported that the Met66 allele conferred susceptibility to AD in women but not in men [11].
The influence of BDNF polymorphisms on human cognition has been linked with its impact on changes in brain structure. However, whether the Met allele leads to worse outcomes is controversial. Smaller brain volumes have been reported in Met carriers in the hippocampal and dorsolateral prefrontal areas [12], hippocampus [13], hippocampal, amygdala, thalamus, fusiform gyrus and frontal gyrus [14]. Other studies have reported no significant differences in hippocampal volume, whole brain volume or memory scores between Met carriers and Val homozygotes [15], and no significant genetic effects of BDNF on cognitive performance [16]. A protective role of the Met allele in regional gray matter (GM) volume and neurocognitive performance have been reported [17][18][19]. Given the role of BDNF as a crucial mediator in the maintenance of neuronal function, it is not known whether the Val66Met genotype is related to large-scale network modulation in patients with AD.
Recent research has suggested that highly related regions show covariance in morphometric characteristics, so called structural covariance. Structural covariance networks (SCNs) can be used to test the influence of the genotype with careful control of other factors. Three SCNs have been reported to be relevant to patients with AD: default mode network (DMN) [20][21][22], salience network [23] and executive control network [24,25]. A recent report suggested that the DMN may be comprised of multiple, spatially dissociated but interactive components [26], of which two subsystems are particularly relevant: the "medial temporal lobe subsystem", and the "dorsal medial prefrontal cortex subsystem" (or the midline core subsystem).
The potential mechanisms of genetic-based neurobiology are still under investigation, however a number of studies have highlighted how genetic variations [27][28][29][30], epigenetic [31,32] or metabolic interactions [33] may affect organization of the brain or therapeutic programs [34]. In this study, we hypothesized that BDNF Val66Met functional polymorphisms may modulate the large-scale structural covariance pattern in patients with AD, and that network alterations may also determine the neurobehavioral characteristics.

results demographic data, cognitive data and nPI
The demographic characteristics and neuropsychiatric test results of the three genotype groups are listed in Table 1. The Val66Met genotype distribution was in Hardy-Weinberg equilibrium. There were no significant differences in apolipoprotein E4 allele distribution, MMSE scores, total or individual CASI scores, NPI total scores or levels of cerebrovascular risk biomarkers among the groups.

Patterns of structural associations in the patients and genetic variants
Among the three BDNF genotypes and four seed regions ( Figure 1A), there were no significant differences in the GM volumes of each seed ( Figure 1B). Networks showing structural associations with the seed regions for each genotype are shown in Figure 1C and Supplementary

seed region volumes and relationships with the cognitive scores
We first explored whether each seed region volume was correlated with the selected cognitive www.impactjournals.com/oncotarget test (Supplementary Figure 1). The DMN entorhinal (Supplementary Figure 1A) and PCC (Supplementary Figure 1B) seed volumes both showed significant correlations with the MMSE scores, CASI total scores, short-term memory, orientation and EFT scores, while the PCC seed also showed significant correlations with CASI language and drawing scores. For the salience network, only the CASI short-term memory score showed a significant correlation with the frontoinsular seed volume (Supplementary Figure 1C). For the executive control network, the seed region volume showed significant correlations with all test scores (Supplementary Figure  1D). The correlations between the seed volume and test scores suggested greater clinical significance of the PCC seed and dorsolateral prefrontal seed on the prediction of overall cognitive test scores.

Peak clusters showing significant interactions between genotypes (Val/Val > Met/Met or Val/Met > Met/Met)
For each seed, we further explored the genotypic interactions with regards to the topography showing differences in structural covariance between seed and peak clusters ( Figure 2 and Supplementary Figure 2).
Only peak clusters connected to the right entorhinal seed ( Table 2) showed topographic similarities in structural covariance in genotype group analysis (  Data are presented as mean (standard deviation) or number (percentage; %) Abbreviations: CASI, Cognitive Ability Screening Instrument; Attention, verbal fluency, abstract thinking, and mental manipulation sub-domain scores of the CASI were used to assess executive function. BDNF, brain APOE4 carriers were defined as the presence of one or two APOE4 alleles.  Figure 2D).   There was no significant difference in seed volume among the three genotypes (p > 0.05). Z-statistic maps (p < 0.05, corrected with a false discovery rate with extended cluster voxels > 100). The images are displayed on a standard brain render. (Supplementary Tables 13-15). Of note, in the covariance linking the PCC, frontoinsular and dorsolateral prefrontal seeds, there were no overlapping regions of Val/Val > Val/ Met > Met/Met, suggesting a lack of a genetic dosage effect from these three seeds.

dosage effects (Val/Val > Val/Met > Met/Met) of the bdnF val allele via entorhinal seed-based connectivity covariance strength alterations
To investigate whether there was a genetic dosage effect, the volumes of the five peak clusters that were connected with the entorhinal seed were extracted, as they showed significant interactions between Val/Val > Met/Met or Val/Met > Met/Met (Figure 3). While the  entorhinal seed volume ( Figure 1B) and seed-connected cluster volume ( Figure 3A-3E) were not significantly different among the three genotypes, the genetic dosage effect of the Val allele was via changes in the strength of entorhinal seed-based covariance connectivity (Val/Val > Val/Met > Met/Met) (Figure 3).

Clinical significance of entorhinal seed-related peak clusters showing genotype differences
The clinical significance of each peak cluster linking to the entorhinal seed showing a genotype effect was evaluated by correlation analysis (Figure 4). Among the five peak cluster volumes and in all three genotypes, only the left orbitofrontal cluster volume ( Figure 4A, bold lines) demonstrated significant correlations in predicting MMSE, CASI total scores, CASI orientation and EFT scores.

dIscussIon
The results of this study provide data on the network-specific genetic influence of BDNF Val66Met in the early stages of AD. There were three main findings. First, the seed-based SCN pattern validated the hypothesis that the BDNF functional polymorphism targets largescale brain networks rather than selectively in a focal brain region. Second, a genetic dosage effect of the BDNF Val allele was found only in the DMN medial temporal subsystem with cortical hubs consisting of entorhinal, PCC and prefrontal regions. The presence of the genetic dosage effect suggests that the mechanism related to this functional polymorphism may be via modulation of the strength of structural covariance between spatially scattered but functionally coherent regions rather than having a direct impact on the seed or peak cluster volume. Third, while the volumes of the seed and peak clusters of the DMN medial temporal subsystem correlated significantly with cognitive test scores, only the orbitofrontal peak cluster volume predicted cognitive test scores in the three BDNF genotype groups.

Network-specific genetic influence on DMN medial temporal lobe subsystem
The salient episodic memory impairment in AD has been reported to be strongly associated with hippocampal/ entorhinal volume [35,36], while decreased expressions of BDNF and its receptor, tyrosine receptor kinase B, were also located in the hippocampal and frontal regions [37]. Pathological molecules of the amyloid or tau protein may spread via synaptic connections, and BDNF has been reported to attenuate the neurotoxicity [5,8]. The correlation of structural covariance between regions could therefore reveal the genetic effect of AD at the network level. Taken together with our results, the clinical features, BDNF activity and hippocampal network appear to be highly correlated in the early stages of AD.
While our results suggest that there was no direct effect of BDNF genotype on hippocampal or entorhinal volume, we found that the BDNF genotype may specifically modulate the DMN medial temporal lobe subsystem structural covariance patterns, of which the increased correlations between entorhinal seed and peak clusters demonstrated the genetic dosage effect of the Val allele. While early amyloid burden or GM atrophy has been reported to involve the DMN medial temporal lobe subsystem in the early stage of AD [38], lower BDNF protein levels in the entorhinal cortex seed region have also been reported [39] implying the vulnerability of this region. Therefore, within the entorhinal seed and DMN, BDNF activity can be regarded as an important factor against pathological toxicity.
These alterations can be caused by several factors, and the disease-and genetic-related changes in structural covariance can be difficult to interpret. As the cortical hubs of the DMN are highly functionally anchored, the spatially scattered DMN-related clusters that show structural covariance may validate the notion that BDNF Val66Met polymorphisms only affect activity-dependent pathways [9]. Attenuated correlations between brain regions may suggest disconnectivity related to lower BDNF activity. The explanation of the dosage effect using correlation strength between seed and seed-based connectivity clusters was based on reports that neurons transfected with Met-BDNF show diminished neuronal integrity [9] and less hippocampal dendritic complexity in haploinsufficient BDNF mice [40]. In transgenic mice (BDNF [Met/Met]) [41], normal levels of BDNF have been reported in the brain despite defective secretion from neurons, suggesting the genetic susceptibility of Met homozygotes. As such, the seed region of Met carriers may show less connectivity in mediating the cortical hubs of the DMN network. Alternatively, localized or limited degeneration of a covariance network is also possible in patients with the Met allele [42].
In contrast, stronger correlations modulated by increased BDNF activity may suggest greater regional connectivity and synchronized GM loss in regions targeted by the pathological process. Although the pairwise correlations did not necessarily influence the global organizational properties of structural covariance, as shown by the lack of differences in regional volumes among genotypes, the significant correlations between cognitive outcomes and seed (or peak clusters such as the anterior cingulate, orbitofrontal and dorsolateral prefrontal lobes) volumes within the DMN hubs still address the pathological-genetic alterations of local clustering. The findings of the DMN still reflect segregated and less integrated components in AD networks that could predict cognitive outcomes.

bdnF genotype polymorphisms target different GM degenerative network patterns
Among the four seed-based SCNs, the Val (Val/ Val and Val/Met) carriers had more extended voxels of structural covariance in the DMN medial temporal lobe subsystem. In contrast, the Met (Met/Met and Val/ Met) carriers showed stronger structural associations in the DMN midline core system, salience and executive control networks. As BDNF activity was different among the genotypes, the SCN patterns may reflect different strengths of structural covariance between cortical GM and seed region. The SCN patterns may therefore provide evidence that BDNF genotypes target specific large-scale brain networks in patients with AD that may modulate distinct degenerative patterns. In normal human cortical morphology, Met carriers have been reported to have smaller hippocampal and prefrontal cortical volumes [12,43]. However, the genetic effects of Met carriers on smaller entorhinal or prefrontal cortical volumes were not established in this study.
Studies on the genetic risk of developing AD [10,11] or the direct influence of BDNF Val66Met functional polymorphisms on regional volume have reported inconsistent results [12][13][14][15]44]. As significant apolipoprotein E and BDNF interactions have been reported in normal elderly and AD patients to predict episodic memory performance [45][46][47][48], our study design may ignore the epistatic effect of the apolipoprotein E gene. However, linkage disequilibrium either from other functional polymorphisms in the BDNF gene or from another nearby gene was still possible. Taken together, the strength of connectivity alterations with regards to the interactions of genotypes supports the hypothesis that BDNF functional activity contributes to cognitive outcomes in the early stages of AD, and that this may be mediated by polymorphisms at the Val66Met locus.
Possible explanations why the BDNF genotypes showed inconsistent results on brain structures may be related to epistasis [44], variability in individual hippocampal activity and its neural substrates [49], differences in racial distribution of the Met allele with differences in group strategies for BDNF genotypes [50], and interactions between genotypes and specific intracerebral pathology [51]. In surgically-resected hippocampi from patients with epilepsy, increased levels of BDNF mRNA and protein have been noted, indicating that epileptic activity may up-regulate protein levels via BDNF gene expression [52,53]. In patients with chronic epilepsy, however, decreased BDNF secretion has been associated with a decompensated hippocampus [54].As such, secretion of BDNF can be dynamic depending on the integrity of the hippocampus, disease typs and thus the disease stage. We enrolled AD patients with an early stage of disease and the network influence could be different at different stages of disease.

GM alterations reflected the endophenotype of the bdnF genotype
Despite our observations that BDNF Val66Met polymorphisms may affect structural covariance patterns, there were no significant differences in cognitive measurements in our genotype groups. On the basis of the none-significance in the PCC and dorsolateral prefrontal seed volumes among the three genotype groups, the volumes of both seeds were still significantly correlated with the cognitive scores. The finding supports the hypothesis that brain structure may be considered to be an endophenotype that is more sensitive to a genetic effect than behavioral level. The PCC is an early target of amyloid deposition in both patients with AD [55] and in cognitively normal subjects with a positive family history of AD [56]. The dorsolateral prefrontal cortex has been reported to frequently confer executive function, and activation of the prefrontal neural resources have been reported to compensate for the posterior degenerative process [57].

study limitations
An important limitation of this study is that we did not include a control group. The effects of BDNF Val66Met genetic variations on regional GM volume and functional state in healthy subjects have been reported [16,47,50]. The enrolment of controls may help to understand how the disease interacts with degenerative processes, and whether different genotypes have similar effects on normative structures. Nonetheless, the aim of this study was to determine whether changes in structural correlations between areas of the brain are AD-specific. Our results may provide further evidence reported in healthy elders, that genetic variations of BDNF in AD may mediate SCN patterns and the strength of structural covariance rather than focally. Another potential limitation is that we used seed-based analysis, with an emphasis on the SCNs and reported the genetic dosage effect focusing on the Val allele. The use of independent component analysis [58] or resting state functional MRI data may elucidate whether other networks are involved and also validate the findings of the DMN medial temporal network observed in this study. Finally, as the clinical significance was established in restricted nodes showing peak correlation, whether these nodes were affected by other genetic or biological effects remains an important issue that needs to be tested in further studies. Nonetheless, our genotype groups were carefully matched in major possible confounders, highlighting the unique role of this BDNF functional polymorphism.

conclusIons
In summary, our structural covariance analysis supports the role of BDNF functional polymorphisms in modulating the GM degenerative scaffold in the early stages of AD. Our results suggest that the Val66Met functional polymorphism carries different weightings on the DMN, salience and executive control networks. Within the DMN medial temporal lobe subsystem, the genetic dosage effect of the Val allele on entorhinal seedbased SCN was validated by correlation strength gradients between seed and peak clusters, thus supporting the importance of BDNF in AD degenerative network.

MAterIAls And Methods
This study was conducted in accordance with the Declaration of Helsinki and was approved by the Institutional Review Board of Chang Gung Memorial Hospital. The study participants were treated at the Cognition and Aging Center, Department of General Neurology, Kaohsiung Chang Gung Memorial Hospital. A total of 186 subjects (82 males, 104 females) were included after the consensus of a panel composed of neurologists, neuropsychologists, neuroradiologists and experts in nuclear medicine [59]. AD was diagnosed according to the International Working Group criteria [60] with a clinical diagnosis of typical AD. All of the patients had a clinical dementia rating scores of 0.5 or 1, and all of the patients were in a stable condition under treatment with acetylcholine esterase inhibitors from the time of diagnosis. The exclusion criteria were a history of clinical stroke, a modified Hachinski ischemic score > 4 [61] and depression.

study working scheme
The patients were classified into three groups based on the genotype: Met/Met homozygotes (n = 45), Met/Val (n = 86) and Val/Val homozygotes (n = 55). To avoid possible confounding, the three genotype groups were matched for gender, age, educational level, clinical dementia rating [62] and Mini-Mental State Examination (MMSE) scores [63] ( Table 1). The working scheme was as follows. First, the SCN was established by seed-based correlation analysis. Seed regions anchoring the DMN (medial temporal subsystem or midline core subsystem), salience and executive control networks were selected to generate SCNs. Differences in each seed regional volume and SCN peak cluster volume were then compared among the three genotype groups. Finally, to evaluate the genetic dosage effects of the Val allele on SCN clusters [9], the SCNs showing significant genotype interactions (i.e. Val/ Val > Met/Met or Met/Val > Met/Met) in seed-peak cluster covariance were modeled. Correlations between peak clusters showing a dosage effect with cognitive test scores were assessed to evaluate the clinical significance.

clinical and neurobehavioral assessments
After enrolment, demographic data and a family history of each patient were recorded, and each patient underwent physical and neurological examinations. A trained neuro-psychologist administered the tests. The 30-item MMSE [63] and cognitive ability-screening instrument (CASI) [64] total scores were used as a global assessment of cognitive function. Attention, verbal fluency, abstract thinking, and mental manipulation sub-domain scores of the CASI were used to assess executive function test (EFT) [65], while the non-executive domains included orientation, short-and long-term memory, language ability, and drawing. We used the neuropsychiatric inventory (NPI) to evaluate changes in behavior.

Genotyping
Genomic DNA was extracted from blood samples using a commercial kit (Qiagen, Gentra Puregene Blood Kit), followed by genotyping for BDNF Val66Met polymorphisms using the polymerase chain reaction (PCR)-restriction fragment length polymorphism method. Genotyping was conducted with the operator blinded to the clinical data. The apolipoprotein E4 genotype was also determined using a PCR-restriction fragment length polymorphism assay and restriction enzyme HhaI [66]. Apolipoprotein E4 carriers were defined as those with one or two E4 alleles [59].

Image acquisition
MR images were acquired using a 3.0T MRI scanner (Excite, GE Medical Systems, Milwaukee, WI, USA). Structural images were acquired for structural covariance analysis using the following protocols: a T1weighted, inversion-recovery-prepared, three-dimensional, gradient-recalled acquisition in a steady-state sequence with a repetition time/echo time/inversion time of 8,600 ms/minimal/450 ms, a 256 × 256 mm field of view, and a 1-mm slice sagittal thickness with a resolution of 0.5 × 0.5 × 1 mm 3 .

data analysis
Image preprocessing and statistical analysis were performed using SPM8 (SPM8, Wellcome Trust Centre of Cognitive Neurology, University College London, UK, http://www.fil.ion.ucl.ac.uk/spm/). The T1 images were reoriented, realigned, and normalized using the standard Montreal Neurological Institute (MNI) space. The images were then segmented into GM and white matter. Related tissue segments were used to create a custom template using the diffeomorphic anatomical registration using exponentiated lie algebra (DARTEL) approach. The DARTEL approach is one of the highest ranking registration methods in patients with AD [70]. The modulated and warped images were then smoothed using a Gaussian kernel of 8 mm full width at half maximum. A direct comparison among the modulated segmented GM volumes of Val/Val homozygotes, Met/Val and Met/ Met homozygotes using voxel-based morphometry [71] showed no significant differences with the threshold set at www.impactjournals.com/oncotarget p < 0.05, corrected for a false discovery rate (FDR) with a cluster size > 100 voxels.

statistic analysis
Clinical and laboratory data were expressed as mean ± standard deviation. Analysis of variance with Bonferroni correction for multiple comparisons was used to compare levels of cerebrovascular risk biomarkers or continuous variables among the Val/Val homozygotes, and Met/ Val and Met/Met homozygotes. All statistical analyses were conducted using SPSS software (SPSS version 22 for Windows ® , SPSS Inc., Chicago, IL). Statistical significance was set at p < 0.05.
From the modified GM images, the GM volumes of a 4-mm radius sphere around the seed ROI coordinates were calculated, followed by four separate correlation analyses using the extracted GM volumes as the covariates of interest. The three BDNF genotype groups were modeled separately. For each BDNF genotype group, specific contrasts were set to identify voxels that showed positive correlations for each seed ROI. The results reflected the SCNs of each ROI and the threshold was set at p < 0.01, corrected for FDR with a cluster size > 100 voxels.
In addition, to investigate how genetic variance may interfere with structural covariance patterns, voxels showing significant differences in the regression slopes in each seed-peak cluster correlations were compared, pointing to possible interactions between Val/Val > Met/ Met or Val/Met > Met/Met. The genetic dosage model was based on an in vivo study in which the expression of BDNF was found to be highest in Val/Val, followed by Val/Met and Met/Met [75]. Specific t contrasts were established to map the voxels that expressed significant between-group associations. The threshold for the resulting statistical parametric maps was set at p < 0.001 (uncorrected) with a cluster size > 100 voxels. In addition, for the peak clusters showing significant between-group differences, a 4-mm radius sphere was placed on the peak voxel, and the GM volumes were then calculated for regression analysis. The seed ROI was considered to be the predictive variable for the extracted SCN peak voxel, and the threshold was set at p < 0.05 with multiple corrections. To evaluate the clinical significance of the seed or identified peak voxel, we used a linear regression model with the cognitive test scores serving as the dependent variable. The threshold was set at p < 0.05 with multiple corrections.