Risk Variants in Three Alzheimer’s Disease Genes Show Association with EEG Endophenotypes
Abstract
Background:
Dementia due to Alzheimer’s disease (AD) is a complex neurodegenerative disorder, which much of heritability remains unexplained. At the clinical level, one of the most common physiological alterations is the slowing of oscillatory brain activity, measurable by electroencephalography (EEG). Relative power (RP) at the conventional frequency bands (i.e., delta, theta, alpha, beta-1, and beta-2) can be considered as AD endophenotypes.
Objective:
The aim of this work is to analyze the association between sixteen genes previously related with AD: APOE, PICALM, CLU, BCHE, CETP, CR1, SLC6A3, GRIN2 β, SORL1, TOMM40, GSK3 β, UNC5C, OPRD1, NAV2, HOMER2, and IL1RAP, and the slowing of the brain activity, assessed by means of RP at the aforementioned frequency bands.
Methods:
An Iberian cohort of 45 elderly controls, 45 individuals with mild cognitive impairment, and 109 AD patients in the three stages of the disease was considered. Genomic information and brain activity of each subject were analyzed.
Results:
The slowing of brain activity was observed in carriers of risk alleles in IL1RAP (rs10212109, rs9823517, rs4687150), UNC5C (rs17024131), and NAV2 (rs1425227, rs862785) genes, regardless of the disease status and situation towards the strongest risk factors: age, sex, and APOE ɛ4 presence.
Conclusion:
Endophenotypes reduce the complexity of the general phenotype and genetic variants with a major effect on those specific traits may be then identified. The found associations in this work are novel and may contribute to the comprehension of AD pathogenesis, each with a different biological role, and influencing multiple factors involved in brain physiology.
INTRODUCTION
Dementia due to Alzheimer’s disease (AD) is a common, age-related, neurodegenerative disorder leading to progressive memory loss and impairments in speech and behavior. The prevalence of AD is higher among females than males, even when individuals at the same age are compared [1–4]. Indeed, AD diagnosis is not straightforward in early stages, and symptoms are frequently dismissed and confused with normal aging. AD progression is often characterized in three main stages: mild (or early-stage), moderate, and severe (or late-stage) AD, depending on the level of functionality of the patient [5]. Mild cognitive impairment (MCI) is characterized by a slight but perceptible decline in cognitive abilities, including memory and thinking skills, and plays an important role in early diagnosis. Individuals with MCI do not necessarily develop AD, but are in greater risk, especially those with amnestic MCI, which is considered a prodromal stage of the disease [6, 7].
Heritability plays an important role in AD [8], APOE having been reported since the early 1990s as the strongest genetic risk for the disease, through allele ɛ4 (reviewed in [9]). From the genetic point of view, two subtypes of the disease are generally considered: 1) familial early-onset AD, and 2) late-onset AD (LOAD). LOAD appears in sporadic cases after 65 years old, and is approximately 20 times more prevalent than the familial early-onset subtype [10]. Heritability of familial AD is often explained by the presence of rare variants in a few genes, which are expected to have a strong effect [11–13]. Contrarily, LOAD is generally associated with common variants, expected to have a small impact when analyzed individually, but a large effect when aggregated [14, 15]. Thus, to reliably identify the individual contribution of a single polymorphism to LOAD risk, huge genetic datasets are needed, which may carry problems of ultrahigh dimensionality when a correlation between genome-wide information and the LOAD phenotype is attempted. The analysis of endophenotypes has proven to be a useful approach to simplify the general AD phenotype, helping to dissect the role of specific genetic variants correlated with the disease [16, 17].
Electroencephalography (EEG) measures the electrical activity of the brain, acquiring voltage fluctuations (derived from synaptic potentials) by electrodes placed on the scalp. In previous works, physiological alterations caused by neurodegeneration have been studied by means of EEG analyses [18–20]. AD is known to be strongly correlated with a general slowing of the EEG, measured through relative power (RP) calculations. Specifically, a progressive increase of RP in low frequency bands (i.e., delta and theta), along with a progressive decrease at high frequency bands (i.e., alpha and beta), have been consistently associated with AD progression [21–23]. These alterations seem to be related with neurophysiological and anatomical disturbances, such as hippocampal atrophy [24], cortical disinhibition or hyper excitability [25, 26]. For comprehensive reviews on this topic see [27–30].
RP calculations at each frequency band are objective and reliable measurements, and thus effective as AD endophenotypes [31–34]. Genetic variants with a major effect on those specific traits may be then identified, as analyzed in previous studies considering four AD risk variants in APOE (rs7412 and rs429358), PICALM (rs3851179), and CLU (rs11136000) genes, and cohorts composed by AD patients, MCI subjects, relatives of AD patients, or cognitively healthy young or elderly controls [35–38]. Two of these studies revealed conflicting results in what concerns the presence of APOE ɛ4 allele in AD patients [35, 38]. In [38], a slowing of EEG oscillatory activity in AD patients carrying APOE ɛ4 alleles (rs7412 and rs429358) was observed through a decreased power in alpha band, the same occurring for AD patients’ relatives along with an increased power in delta and theta bands. Contrarily, in [35], AD ɛ4 non-carrier patients showed a decrease in alpha frequency band and an increase in delta one. On the other hand, presence of CLU (rs11136000) and PICALM (rs3851179) risk alleles in cognitively healthy elders was associated with increased alpha and beta RP values [36, 37].
The main goal of this work is to analyze the correlation between: (a) the EEG relative power in the conventional frequency bands (i.e., delta, theta, alpha, beta-1, and beta-2) measured in resting state, and (b) genetic variants in sixteen candidate genes previously associated with AD in an Iberian cohort composed by AD patients, individuals with MCI and controls. To achieve this, the correlation between EEG measurements and both (c) the strongest risk factors for AD: age, sex, and APOE ɛ4 presence, and (d) the status of the subjects regarding AD, was analyzed and RP values corrected accordingly for further analyses.
MATERIAL AND METHODS
Subjects
We studied an Iberian cohort of: 1) LOAD patients in different stages of the disease: mild (MIL), moderate (MOD), and severe (SEV) AD; 2) individuals with MCI; and 3) cognitively healthy elderly controls (CON). AD patients and MCI individuals were clinically diagnosed following the criteria of the National Institute on Aging and Alzheimer’s Association (NIA-AA) [39]. AD staging was based mainly on the Mini-Mental State Examination (MMSE) [40]. The control group was composed by individuals over 68 years old, with no signs of dementia or history of neurological disease. Neither patients nor controls were diagnosed with any other neurologic and psychiatric diseases (other than AD and MCI), and were not using any drugs that might affect EEG signal.
At the time of sample collection, subjects were residents of the autonomous community of Castile and Leon, northwestern Spain, or of the northern region of Portugal. Saliva and buccal swabs were selected as sources of biological sampling to maintain the process as noninvasive as possible. Biological samples and EEG data were collected from 253 individuals, mostly equally distributed by population: Portugal (PT) and Spain (SP), and subgroup: controls (25 PT + 26 SP), individuals with MCI (26 PT + 25 SP), and patients with mild (25 PT + 26 SP), moderate (25 PT + 25 SP), and severe LOAD (25 PT + 25 SP). This project was approved by the Ethics Committee of the University of Porto (report # 38/CEUP/2018), Portugal, and written informed consent was obtained from all participants, family and/or legal representatives.
After DNA quality control assessment, prior to genotyping (with a minimum quantity of 10 ng/μl in 45μl minimum volume and integrity of 90%of gDNA greater than 10 Kb in size; see below for more details), 54 subjects were removed before microarray processing since their biological samples did not fulfill the minimum requirements. The distribution of the 199 subjects which biological samples passed the quality control procedures are presented in Table 1, along with information on sex, mean age, and MMSE scores. Situation of the cohort regarding the strongest risk factors for AD: sex, age, and APOE ɛ4 allele, was analyzed (Supplementary Tables 1 and 2). A Dunn’s test of independence supported that individuals with MCI are significantly older than both AD patients (p = 1.0e–03) and controls (p = 1.5e–04), and a Chi-square test of independence showed a significant association between gender and AD status (p = 1.4e–02) with more female patients and male controls than expected (Supplementary Tables 1 and 2). Also, the distribution of APOE ɛ4 allele carriers significantly differed among sampling groups (p = 2.8e–02), with more carrier patients and non-carrier MCI individuals and controls (Supplementary Table 2).
Table 1
Number of subjects that qualified for further analyses after genotyping quality control, and respective mean age and MMSE scores
| Subgroups | Cases | Total | |||||||||
| CON | MCI | MIL | MOD | SEV | |||||||
| Population | PT | ES | PT | ES | PT | ES | PT | ES | PT | ES | |
| Females | 11 | 12 | 16 | 15 | 15 | 12 | 15 | 18 | 5 | 17 | 136 |
| Males | 10 | 12 | 7 | 7 | 6 | 12 | 2 | 4 | 1 | 2 | 63 |
| Total | 45 | 45 | 45 | 39 | 25 | 199 | |||||
| Mean age±SD (y) | 79.7±7.3 | 84.8±7.1 | 80.7±6.5 | 81.3±8.1 | 79.9±6.4 | 81.4±7.3 | |||||
| Mean±SD MMSE | 28.9±1.1 | 23.3±3.0 | 21.7±3.3 | 13.0±3.7 | 3.8±4.2 | - | |||||
CON, controls; MCI, individuals with mild cognitive impairment; MIL, patients with mild Alzheimer’s disease; MOD, patients with moderate Alzheimer’s disease; SEV, patients with severe Alzheimer’s disease; PT, subjects residents of northern Portugal; SP, subjects residents of Castile and Leon, Spain.
Genotyping
A sample of saliva or three buccal swabs were collected from each of the initial 253 participants of the study. Preference was given for collecting a saliva sample of 2 ml with the Oragene DNA 500 collection kit (DNAgenotek). Buccal swabs were only used for patients at a more advanced stage of the disease, unable to voluntarily spit. After DNA extraction and quality control assessment, samples were genotyped using the genome wide Axiom Spain Biobank Array (Thermo Fisher Scientific) at the Spanish National Center for Genotyping (CEGEN, Santiago de Compostela, Spain, CEGEN-PRB3-ISCIII; supported by grant PT17/0019, of the PE I + D+i 2013-2016, funded by ISCIII and ERDF).
Variant calling quality control (QC) was performed in accordance with the Affymetrix’s best practices guide, and a widely used protocol was followed for both individual and marker analysis [41]. All the analyses were computed with Affymetrix Power Tools and PLINK [42]. In the variant calling QC, individuals with dish QC or QC call rates below the defined thresholds were not considered for further analysis, as well as those with heterozygosity rate greater than the defined acceptance threshold. Finally, the probes belonging to the recommended calling categories were selected and the corresponding variants annotated according to the Genome Reference Consortium Human Build 37 (GRCh37) single nucleotide polymorphism (SNP) assembly. In per-individual QC analysis, the sex of the individuals was confirmed through the sex chromosomes’ homozygosity rate, and duplicates or related individuals were identified through the estimation of identity by descent (IBD) rates. For the pairwise cases where IBD estimates revealed to be greater than the established, one of the two individuals was removed (preference for remaining in the study was given first to patients, then to females and finally to individuals with higher call rates). Finally, individuals with divergent ancestry were disregarded. For this, principal component (PC) analyses were computed using the software EIGENSOFT [43, 44] and merging the dataset with one, publicly available, from the 1000 Genomes Project (1KGP) [45], containing 12 different populations from four ancestry groups: East Asian, African, European and Admixed American. In per-marker QC analysis, SNPs with a significant deviation (α= 1E–06) from Hardy-Weinberg equilibrium in control samples were eliminated, as well as those with missingness rate greater than 5%. Significant differences in missing genotype rates between cases and controls (α= 1E–05) were also considered, and SNPs with no variance in the sample were disregarded.
Gene selection
Among the plethora of markers associated with AD, sixteen candidate genes with functional relevance in the brain were selected for further analyses: APOE, PICALM, CLU, BCHE, CETP, CR1, SLC6A3, GRIN2 β, SORL1, TOMM40, GSK3 β, UNC5C, OPRD1, NAV2, HOMER2, and IL1RAP. Some of these markers were selected among those showing strongest association with the disease in traditional case-control genome wide studies [14, 15], while others were selected among those associated with the disease via relevant quantitative endophenotypes, as measures of structural imaging or amyloid-β deposition [49, 55–57, 60–62, 64, 65, 67–69]. The second set of genes were found to be interesting to analyze as traditional case-control studies, considering general (and complex) disease phenotypes, may not uncover genetic risk associations with relevant phenotypes, through which silent signals of the disease may appear before diagnosis being possible. The inclusion of these genes in our analysis could improve the knowledge of their role in brain physiology.
The set of sixteen selected candidate genes and the QC procedures described above, resulted in the analysis of 796 common variants (minimum allele frequency ≥5%) in 199 subjects (Supplementary Table 3). PC analysis showed no population substructure, the gain to the understanding of data with each additional PC revealing to be approximately linear (Supplementary Figure 1). At this point is noteworthy that one of the two variants that determine the APOE genotype, rs7412, failed QC and was further genotyped using a standard Sanger sequencing protocol along with rs429358 SNP (this for validation purposes). The remaining three LOAD risk variants previously associated with EEG measurements in LOAD framework (i.e., APOE rs429358, PICALM rs3851179, and CLU rs11136000) [35–38] passed the QC procedures and their genotypes were considered for further analysis and replication purposes.
Electroencephalographic recordings
Five minutes of resting-state EEG activity was acquired for each subject using a 19-channel Nihon Kohden Neurofax JE-921A System at the electrodes: Fp1, Fp2, Fz, F3, F4, F7, F8, Cz, C3, C4, T3, T4, T5, T6, Pz, P3, P4, O1, and O2, following the International System 10–20. Signals were recorded at a sampling frequency of 500 Hz with common average reference. Subjects were asked to remain awake with closed eyes during acquisition. EEG data were preprocessed according to the following steps [18, 70, 71]: i) mean removal; ii) digital filtering using a Hamming window bandpass finite impulse response (FIR) filter in the band of interest (i.e., 1–30 Hz); iii) independent component analysis (ICA) to remove oculographic and cardiographic artefacts; iv) segmentation into 5 s epochs; and v) visual rejection of epochs with artefacts. In this study, conventional EEG frequency bands were considered: delta (1–4 Hz), theta (4–8 Hz), alpha (8–13 Hz), beta-1 (13–19 Hz), and beta-2 (19–30 Hz). Gamma band was not included in the analyses due to the possible interference of muscle artefacts in its frequency range [72, 73]. In order to quantify the relative contribution of the previous frequency bands to the global power spectrum, RP was computed. Specifically, RP was calculated from the normalized power spectral density by summing the contribution of each spectral component in a specific frequency band [74].
Statistical analysis
Power analysis
Considering the set of 796 variants and the cohort of 199 subjects, the statistical power of the study for identifying variants explaining a range of phenotypic variance was quantified in the non-centrality parameter [75], taking into account the interplay between experimental sample size, allele frequency and effect size (Supplementary Table 4). For detecting a SNP explaining 10%of trait variance, the statistical power reached 75.8%, increasing this figure to 97.3%considering a proportion of variance of 15%(Bonferroni-corrected significance level α= 6.28E – 05).
The variance in phenotype Y: Var(Y), explained by genetic variant X, can be decomposed into two components: Var (Y) = β2Var (X) + σ2, where β is the effect size of Var (X) and σ2 quantifies the remaining variance that can be explained by other factors or genetic variants. The first parcel can be estimated by
Correlation between EEG data and covariates: age, sex, APOE ɛ4 presence, and AD status
The normality of the EEG data was assessed (Shapiro-Wilk test), as well as the correlation between EEG data and the strongest risk factors for AD: sex (Kruskal-Wallis’ test), age (Pearson’s correlation test), and APOE ɛ4 presence (Kruskal-Wallis’ test). Independence analyses were also computed considering the AD status of the subjects (Dunn’s test), either assuming AD patients as a unique group or considering the stage of the disease. RP spatial patterns were also explored considering subjects’ status. All the computations were performed using R and PLINK [42, 77], and the significance threshold considered for rejecting the null hypothesis was α= 0.05, corrected accordingly through Bonferroni’s method. Results reaching nominal statistical significance may be presented and highlighted, but only correlations reaching the corrected statistical significance were considered for further analyses.
Before assessing correlation with the genetic data, RP values of each frequency band were corrected for the covariates significantly correlated with each one of them. This was computed by adding them as covariates in the linear model: E (Y) = β0 + βXX + βGG where β0 is the intercept coefficient, βx is the coefficient of the covariates present in matrix X, and βG is the coefficient of genotype G, i.e., the parameter that quantifies the association between the genotype and the expected value of the outcome Y, which is one of the five brainwaves.
Correlation between genetics and EEG measures
Correlation between (corrected) RP values at each EEG frequency band and the genetic variants exhibited by the subjects was assessed through Kruskal-Wallis testing. Information of linkage disequilibrium (LD) between pairs of genetic variants was obtained from LDlink [78], accessed on July 15, 2020, for “Iberian Populations in Spain”. For each one of the 796 variants where differences reached the significance level α= 0.005, allele analysis were performed to assess the possibility of identifying a ‘risk allele’, consistently associated with the slowing of brain activity. Indeed, the results concerning the relationship between the set of 796 genetic variants and the EEG measurements were presented considering α= 0.005 and α= 0.05/796 = 6.28E–05. Aware of the modest size of the sample compared to the number of variants analyzed, the first significance level was considered to identify candidate variants to be further analyzed in the future.
RESULTS
Correlation between EEG data and covariates: age, sex, APOE ɛ4 presence, and AD status
Evidence to refute the null hypothesis regarding normality of the EEG data was found for all brainwaves (delta p = 4.7E–06, theta p = 5.3E–07, alpha p = 2.0E–07, beta-1 p = 9.5E–09, and beta-2 p = 1.8E–10).
RP values showed no association neither with the sex of the subjects nor with the presence of APOE ɛ4 allele in their genotype (α= 1.0E-02, Supplementary Table 5). On the other hand, age showed to be positively correlated with RP theta values (α= 1.0E-02, Supplementary Table 5).
Statistically significant correlations between EEG-based measurements and sampling groups: AD patients, individuals with MCI and controls, were found for the five frequency bands (α= 3.3E-03, Supplementary Table 6) and were in accordance with previous studies. Particularly, when controls were compared with AD patients, statistically significant differences were found for all frequency bands, with patients showing higher RP values for low frequency bands (i.e., delta and theta), while lower for alpha, beta-1, and beta-2. Noteworthy, theta rhythms differentiated between controls and both MCI and AD groups, at a statistically significant level, which may point to the suitability of RP in theta band as case-control biomarker. We also explored the RP spatial patterns for AD, MCI, and control groups. Supplementary Figure 2 shows the previously mentioned slowing process of brain activity associated to AD and MCI. The statistically significant differences between groups mainly appear at parietal and right frontal areas in delta band, at central and parietal areas in theta band, at temporal areas in alpha band, and across all the scalp in beta-1 band. Beta-2 band only showed statistically significant differences at the P3 channel.
When AD patients were analyzed in subgroups according to the disease stage, delta and alpha RP values showed statistically significant differences within AD subgroups (α=1.0E-03, Supplementary Table 7, Supplementary Figure 3).
For further analyses, theta RP values were then corrected for the age of the subjects, RP values of all frequency bands were corrected for the status of the subject: AD patient, individual with MCI and control, and RP values of delta and alpha were also corrected for the stage of the disease. Correlation analyses considering separately the sampling subgroups were also computed.
Correlation between genetics and EEG measures
Some of the analyzed 796 common variants showed to be correlated with the (corrected) EEG-based measures. Considering a significance level α= 0.005, six variants within NAV2 (rs1425227 and rs862785), IL1RAP (rs10212109, rs9823517 and rs4687150), and UNC5C (rs17024131) genes showed to be correlated with RP values (Table 2). One variant, rs4687150, in IL1RAP, reached a significant correlation after Bonferroni correction, with RP in delta (p = 3.37E–05). The variants within NAV2, rs1425227 and rs862785, showed to be in low LD (R2 = 0.006), as well as the pair of variants rs4687150 and rs10212109 (R2 = 0.002) in IL1RAP. On the other hand, variants rs10212109 and rs9823517 in IL1RAP are in non-negligible LD (R2 = 0.143). The smallest p-values were observed for the variants in IL1RAP with the lowest level of LD: rs4687150 and rs10212109 (R2 = 0.002), regarding delta (p = 3.37E–05) and theta (p = 2.65E–04) frequency bands, respectively. Proportions of variance in phenotype explained by a given variant (PVE) were computed for each of the found significant correlations (Table 2). Despite PVEs not being routinely computed, the range of obtained estimates (from 1.95%to 7.96%) are acceptable in the framework of a complex human trait, although admittedly high. Indeed, polygenic combinations, which may include dozens of variants, can explain a great amount of the variance. This is the case, for example, of the total variance in each lipid trait [79]. The power of our study to detect an association like the one of rs4687150 and delta frequency band (PVE = 6.11%) is ∼34%(see the Statistical Analysis section and Supplementary Table 4). Generically, the larger the effect size of a variant, the larger is expected to be the proportion of variance of the phenotype explained by it, and so the more likely is to detect such association.
Table 2
Genetic variants showing statistically significant differences in RP values of at least one frequency band (Kruskal-Wallis’ test, significance level α= 0.005). RP values corrected for: 1AD, MCI, and CON status, 2MIL, MOD, and SEV stage, or 3 the age of the subjects. The risk allele is the one associated with the EEG slowing, reflected in an increase of RP in delta and theta bands, and a decrease of RP in alpha and beta bands. Proportions of variance in phenotype explained by the genetic variants (PVE) are presented in parentheses
| Gene | Ref. SNP | Risk allele | Significant p-values (PVE) | ||||
| Delta1,2 | Theta1,3 | Alpha1,2 | Beta-11 | Beta-21 | |||
| IL1RAP | rs10212109 | C | 0.000265 (5.41%) | ||||
| IL1RAP | rs9823517 | G | 0.000965 (7.96%) | ||||
| IL1RAP | rs4687150 | T | 3.37E–05* (6.11%) | 0.000647 (3.80%) | |||
| UNC5C | rs17024131 | T | 0.000887 (1.95%) | ||||
| NAV2 | rs1425227 | T | 0.001434 (5.57%) | 0.000292 (6.69%) | |||
| NAV2 | rs862785 | G | 0.000280 (5.00%) | ||||
*Significance level reached after Bonferroni’s correction (α= 6.28e–5).
When the sampling groups: AD patients, individuals with MCI, and controls were analyzed separately, statistically Bonferroni corrected significant differences were found for the three IL1RAP variants in some subgroups and frequency bands (Table 3). Indeed, rs4687150 showed to be correlated with delta RP values in AD patients group (p = 5.62e–05), while rs10212109 and rs9823517, showed to be correlated with theta RP values in the subgroup of controls (p = 4.01e–05 and p = 1.72e–06, respectively).
Table 3
p-values regarding pairwise correlation (Kruskal-Wallis’ test, significance level α= 0.005) between the genetic variants presented in Table 2 and the RP values for each EEG frequency band, considering the sampling groups: AD patients, individuals with MCI, and controls. Theta RP values were corrected for the age of the subjects
| Gene | IL1RAP | UNC5C | NAV2 | ||||
| Ref. SNP | rs10212109 | rs9823517 | rs4687150 | rs17024131 | rs1425227 | rs862785 | |
| Delta | Controls | 4.53e–02 | |||||
| MCI Individuals | |||||||
| AD patients | 5.62e–05* | ||||||
| Theta1 | Controls | 4.01e–05* | 1.72e–06* | 3.98e–04 | |||
| MCI Individuals | 1.38e–02 | ||||||
| AD patients | 5.21e–03 | 6.70e–03 | 2.21e–03 | ||||
| Alpha | Controls | 2.82e–02 | |||||
| MCI Individuals | 9.78e–03 | ||||||
| AD patients | 5.12e–03 | ||||||
| Beta-1 | Controls | 4.21e–02 | 1.97e–02 | ||||
| MCI Individuals | 4.55e–03 | ||||||
| AD patients | 1.26e–03 | ||||||
| Beta-2 | Controls | 1.42e–02 | |||||
| MCI Individuals | 2.19e–02 | ||||||
| AD patients | 8.25e–04 | ||||||
1RP values corrected for the age of the individuals; *Significance level reached after Bonferroni’s correction (α= 6.28e–05).
For all the six variants where the significance level α= 0.005 was reached, allele correlation was analyzed and in all the cases a risk allele was associated with the slowing of brain activity (higher delta and theta, and lower alpha and beta mean RP values). It is noteworthy that no evidence of differentiation were found between the frequency of the risk allele among AD patients and controls (Supplementary Table 8). The distribution of genotypes according with AD status is presented in Supplementary Table 9.
Concerning the four variants previously associated with RP values in the context of AD [35–38], only APOE variants rs7412 and rs429358, showed evidence of correlation with RP values in beta-1 frequency band (Supplementary Table 10).
DISCUSSION
Correlation between EEG data and covariates: age, sex, APOE ɛ4 presence, and AD status
When analyzed independently, EEG data provided results in accordance with the literature [80]. As expected, AD patients showed a significant slowing of brain oscillatory activity compared with controls, and within AD subgroup a progressive slowing of EEG was observed along with an increasing of the severity of the disease. Also, theta RP values showed to be correlated with the age of the subjects. From the obtained results is noteworthy that theta rhythms can differentiate between controls and both MCI and AD groups, at a statistically significant level. The suitability of theta band to reflect dementia conditions was already reported in previous studies [81–83]. This agrees with our results, appointing EEG theta band alterations as a potential biomarker to detect hints of neural deterioration, even in pre-clinical states. Indeed, our results supports a likely association between EEG power in theta band and cognitive impairment. Previous studies reported correlation between healthy cognition and reduced tonic theta power [84–86], indicating higher values during resting state as a potential consequence of cognitive impairment. In addition, negative correlation between theta power and hippocampal volume was obtained by means of magnetic resonance imaging [87, 88], which could be related with loss of CA1 hippocampal pyramidal neurons [89]. Another well-known negative correlation in AD is delta and theta power with MMSE [90], which may provide diagnostic value to EEG.
Correlation between genetics and EEG measures
After correlating the genetic information of 796 variants from 16 genes previously associated with LOAD and the RP at each EEG frequency (corrected accordingly for the different covariates), three genes harbored the highest number of variants correlated with EEG-related measures: IL1RAP, UNC5C, and NAV2. Among the genotyped variants are those from APOE (rs7412 and rs429358), P ICALM (rs3851179), and CLU (rs11136000) genes previously associated with RP values [35–38]. Our analysis only showed evidence of correlation between APOE variants and RP values in beta-1 frequency band. Theories considering that hyperactivity and hyperconnectivity in preclinical stages of AD may drive to profound disconnection and hypoactivity in AD patients have been supported [91], which may be the cause of these discrepancies.
Overall, theta RP values showed the highest correlation with the genetic variants tested. Previous studies considered alterations in theta-related EEG activity to be associated with amyloid plaque deposition, which is strongly correlated with AD [92]. Also, correlations between theta band disturbances and known AD biomarkers, such as total-tau (T-tau) and phosphorylated tau (P-tau) presence in cerebrospinal fluid, have been reported [93]. At the physiological level, cerebellar-evoked prefrontal synchronization in theta frequency band showed to be modulated by GABA, and positively associated with working memory performance [94]. This may indicate that a GABAergic dependent set of interneurons have a key role on the cortex, modulating theta-burst stimulation the cortical excitability of distant interconnected cortical areas [95]. Evidence suggests that the cerebellum likely exerts its control on the cortex by a GABAergic dependent set of interneurons and cerebellar theta-burst stimulation modulates cortical excitability of distant interconnected cortical areas [96]. The hippocampal theta oscillations are also functionally relevant in humans since they appear to explain how the hippocampus organizes information from a broad range of neocortical networks, being correlated with behaviors such as memory and spatial navigation [97]. All in all, these observations suggest that theta rhythms are intimately related with complex brain processes that become hindered with neurodegeneration.
Relationship between IL1RAP and RP in delta and theta frequency
A genome-wide association study of longitudinal amyloid accumulation in AD patients implicated the microglial activation of IL1RAP, presenting the rs12053868-G allele association with raised amyloid accumulation and greater declines in temporal cortex thickness (over 2 years). Moreover, MCI carriers of this allele showed to be more likely to progress into dementia due to AD, having a faster 2-year decline in episodic memory performance. Indeed, this gene presents itself as a molecular link between the immune system and synapse formation. In this work, three variants from IL1RAP gene were found to be correlated with RP in at least one frequency band. Variants rs4687150 and rs10212109 are in low LD (R2 = 0.002) and the smallest p-values were found for delta (p = 3.37E–05) and theta (p = 2.65E–04) frequency bands. IL1RAP (Interleukin-1 receptor accessory protein) is a constituent of the IL-1 (interleukin-1), a potent pro-inflammatory cytokine that promotes microglial activation, receptor complex and its downstream signaling pathway, essential components for mediating the immune responses of the IL-1 family of cytokines [98]. Besides neuroinflammation and cell death in neurodegenerative conditions (such as stroke and head injuries), IL-1 signaling is also involved in sleep regulation, learning, memory, and the lipid metabolism in the brain [99–102]. In the brain, IL1RAP exists in two isoforms, IL1RAP and IL1RAPb, with differences only in the C-terminal region [103]. Although this gene is not highly expressed in the brain (https://www.gtexportal.org/), synaptogenic activities of this receptor were found in cultured cortical neurons. Knockdown of both isoforms suppressed synapse formation and knock-out mice showed a decreased of the spine densities of cortical and hippocampal pyramidal neurons, suggesting that these receptors function as cell adhesion molecules and organize synapse formation in the brain [104]. IL-1 is chronically upregulated in dementia due to AD and it is believed to be part of an inflammatory cycle that drives AD pathology. Sustained IL-1β overexpression has been shown to increase tau phosphorylation, despite a substantial reduction in amyloid load [105]. Previous evidence indicates that higher rates of phosphorylation of tau protein is related with neuronal and axonal loss [106]. Also, greater amounts of phosphorylated tau have been found to be associated with alterations in theta rhythms due to lower neuronal excitability [107], probably due to axonal degradation. In addition, higher power in delta band has been attributed to neuronal injury [108], which could be a side effect of the previous point.
Relationship between UNC5C and RP in alpha frequency band
UNC5C has been identified as related to AD in various association studies, suggesting a notable role in the disease [53, 60, 109, 110]. Evidence from cell models show that UNC5C polymorphisms could lead to AD pathogenesis by activating death-associated protein kinase 1 (DAPK1). Meanwhile, DAPK1 itself was shown to be involved in modulating tau protein accumulation, amyloid-β toxicity and neuronal apoptosis/autophagy [53, 111, 112]. This gene belongs to the UNC5H receptor family, including UNC5A, UNC5B, UNC5C, and UNC5D, which are highly expressed in the nervous system (https://www.gtexportal.org/), and seems to be enriched in the hippocampus of AD brain [113]. It is a transmembrane receptor for netrin-1 and can trigger apoptosis under the absence of the netrin-1 ligand. As the receptor of neurin-1, UNC5C plays a crucial role in mediating axon repulsion of neuronal growth cones and cell migration in the developing nervous system [114, 115]. Being UNC5C involved in mechanism during axon development, alterations in this gene could trigger alterations in connections between neuronal groups equivalent to lesions in the cerebral white matter. In this line, a positive correlation between white matter vascular injury and sources of alpha EEG, were found in MCI individuals [116]. Additional evidence showed correlations between alpha power and white matter volume [117]. These findings suggest that alpha activity may be affected by aberrant neuronal pathways caused by altered development mechanisms of the nervous system.
Relationship between NAV2 and RP in theta and beta frequency bands
NAV2 (neuron navigator-2), expressed in the brain, is essential for nervous system development [118]. NAV2 is highly expressed in the hippocampus, cortex, cerebellum, and thalamus [119] and, more specifically, in postmitotic neurons involved in cell migration and neurite outgrowth. There is some evidence that the variant rs1425227 can alter the binding motive of a transcription factor active in bipolar neurons, and thus may have a regulatory effect in NAV2. Indeed, this variant lies within a promoter region, with enhancer histone marks, and within open chromatin in the brain, which suggests active transcription. NAV2 mutant mice embryos have been shown to display a reduction in nerve fiber density, as well as specific defects in cranial nerves, being required for normal cranial nerve development and blood pressure regulation [120]. In this study, NAV2 SNPs are associated mostly with theta and beta band alterations. This gene was previously mentioned to be involved in episodic memory tasks [121]. It turns out that processes of this type are closely related to theta band activity [122–124], which seem reasonable to occur during an EEG recording. Hence, we suggest a greater impact of NAV2 expression in certain cognitive functions under neurodegeneration conditions.
Limitations and future research lines
Some limitations were faced in this work and should be examined in future research. Since this study aimed to ascertain relations between RP values and a reasonable quantity of genetic features, the reliability of the results may be sensitive to database size. This was particularly noticeable when the analyses were computed considering separately the sampling subgroups: AD patients, MCI individuals, and controls. The genetic associations found in this work are novel and future research is important to ensure the correlation with the disease. In order to improve statistical power, enlarging the sample of study should be taken into account. It is also noteworthy that RP was the only feature considered in this study to characterize neurodynamic alterations. Although RP is reliable describing the deterioration process of dementia, connectivity measures (amplitude envelope correlation, weighted phase lag index, etc.) or network parameters (clustering coefficient, characteristic path length and betweenness centrality, among others), could provide additional insights in this regard. In future research, we aim to apply alternative metrics that allow to obtain further information on physiological deterioration associated to the progression of AD. In addition, whole scalp EEG activity was considered, losing spatial-relative data. Several studies pointed particular brain areas to be associated with more marked EEG alterations in AD [125–128], and hence the average activity from all electrodes may diminish the statistical significance of the associations, which in fact may be stronger.
CONCLUSIONS
RP values at the conventional EEG frequency bands (i.e., delta, theta, alpha, beta-1, and beta-2) showed to be correlated with some genetic variants in genes previously associated with AD. Globally, theta frequency band showed the greatest correlation with the analyzed genetic variants and is noteworthy its potential as case-control biomarker.
Novel associations between variants in IL1RAP, UNC5C, and NAV2 genes, and the RP values of some frequency bands (delta and theta, alpha, and theta and beta, respectively) were identified. In all the cases a risk allele was associated with the slowing of brain activity, specifically with the increase of RP values in low frequency bands (i.e., delta and theta), and with a decrease at high frequency bands (i.e., alpha and beta).
These associations may contribute to the comprehension of AD pathogenesis, each with a different biological role, and influencing multiple factors that contribute to brain physiology.
ACKNOWLEDGMENTS
We deeply thank the three anonymous reviewers, whose comments and suggestions deeply contributed to the improvement of this work.
We deeply thank all participants, families and institutions involved, namely: Asociación de Familiares de Enfermos de Alzheimer de Ávila, Ávila; Associação de Pensionistas e Reformados de Viana do Castelo, Viana do Castelo;Casa do Povo de Alvito S.Pedro, Barcelos; Santa Casa da Misericórdia de Vila Nova de Gaia; Obra Social Nossa Senhora da Boa Viagem, Porto; Gero Vida, Villaralbo (Zamora); Asociación de Familiares de Alzheimer de León; Residencia San Raimundo en Coreses; Centro de Dia S. João de Deus, da Santa Casa da Misericórdia do Porto; Lar Santa Rita, da Santa Casa da Misericórdia de Caminha; Centro Social e Cultural de Vila Praia de Âncora; Lar Casa de Magalhães; Freixo, Ponte de Lima; Centro de Dia Memória de Mim; Lavra; Armonía Centro de Día, Zamora.
This work was supported by ‘European Commission’ and ‘European Regional Development Fund’ under the project “Análisis y correlación entre el genoma completo y la actividad cerebral para la ayuda en el diagnóstico de la enfermedad de Alzheimer” (Project 0378_AD_EEGWA_2_P), (Co-operation Programme INTERREG V-A Spain-Portugal POCTEP 2014-2020) and the COMPETE 2020-Operacional Programme for Competitiveness and Internationalisation (POCI), Portugal 2020. Portuguese funds are supporting this work through FCT-Fundação para a Ciência e a Tecnologia/Ministério da Ciência, Tecnologia e Inovação in the framework of the project “Institute for Research and Innovation in Health Sciences” (POCI-01-0145-FEDER-007274). SM, AML, IG, and NP are funded by FCT: CEECIND/00684/2017, IF/01262/2014, and CEECIND/02609/2017, and through the Decreto-Lei nº 57/2016 de 29 de Agosto, respectively. Spanish funds are supporting this work through ‘Ministerio de Ciencia e Innovación – Agencia Estatal de Investigación’ and ‘European Regional Development Fund’ under project PGC2018-098214-A-I00 and by ‘CIBER en Bioingeniería, Biomateriales y Nanomedicina (CIBER-BBN)’ through ‘Instituto de Salud Carlos III’ co-funded with ‘European Regional Development Fund’ funds. MA is funded by the Grant RYC-2015-18241 from the Spanish Government.
Authors’ disclosures available online (https://www.j-alz.com/manuscript-disclosures/20-0963r2).
SUPPLEMENTARY MATERIAL
[1] The supplementary material is available in the electronic version of this article: https://dx.doi.org/10.3233/JAD-200963.
REFERENCES
[1] | Altmann A , Tian L , Henderson VW , Greicius MD , Alzheimer’s Disease Neuroimaging Initiative ((2014) ) Sex modifies the APOE-related risk of developing Alzheimer disease. Ann Neurol 75: , 563–573. |
[2] | Beam CR , Kaneshiro C , Jang JY , Reynolds CA , Pedersen NL , Gatz M ((2018) ) Differences between women and men in incidence rates of dementia and Alzheimer’s disease. J Alzheimers Dis 64: , 1077–1083. |
[3] | Lin KA , Choudhury KR , Rathakrishnan BG , Marks DM , Petrella JR , Doraiswamy PM , Alzheimer’s Disease Neuroimaging Initiative ((2015) ) Marked gender differences in progression of mild cognitive impairment over 8 years. Alzheimers Dement (N Y) 1: , 103–110. |
[4] | Rocca WA , Grossardt BR , Shuster LT ((2014) ) Oophorectomy, estrogen, and dementia: A 2014 update. Mol Cell Endocrinol 389: , 7–12. |
[5] | Reisberg B , Ferris SH , de Leon MJ , Crook T ((1982) ) The Global Deterioration Scale for assessment of primary degenerative dementia. Am J Psychiatry 139: , 1136–1139. |
[6] | Grundman M , Petersen RC , Ferris SH , Thomas RG , Aisen PS , Bennett DA , Foster NL , Jack CR Jr. , Galasko DR , Doody R , Kaye J , Sano M , Mohs R , Gauthier S , Kim HT , Jin S , Schultz AN , Schafer K , Mulnard R , van Dyck CH , Mintzer J , Zamrini EY , Cahn-Weiner D , Thal LJ ((2004) ) Mild cognitive impairment can be distinguished from Alzheimer disease and normal aging for clinical trials. Arch Neurol 61: , 59–66. |
[7] | Petersen RC ((2004) ) Mild cognitive impairment as a diagnostic entity. J Intern Med 256: , 183–194. |
[8] | Brainstorm Consortium, Anttila V , Bulik-Sullivan B , Finucane HK , Walters RK , Bras J , Duncan L , Escott-Price V , Falcone GJ , Gormley P , et al. ((2018) ) Analysis of shared heritability in common disorders of the brain. Science 360: , eaap8757. |
[9] | Belloy ME , Napolioni V , Greicius MD ((2019) ) A quarter century of APOE and Alzheimer’s disease: Progress to date and the path forward. Neuron 101: , 820–838. |
[10] | Zhu X-C , Tan L , Wang H-F , Jiang T , Cao L , Wang C , Wang J , Tan C-C , Meng X-F , Yu J-T ((2015) ) Rate of early onset Alzheimer’s disease: A systematic review and meta-analysis. Ann Transl Med 3: , 38–38. |
[11] | Goate A , Chartier-Harlin MC , Mullan M , Brown J , Crawford F , Fidani L , Giuffra L , Haynes A , Irving N , James L , et al. ((1991) ) Segregation of a missense mutation in the amyloid precursor protein gene with familial Alzheimer’s disease. Nature 349: , 704–706. |
[12] | Sherrington R , Rogaev EI , Liang Y , Rogaeva EA , Levesque G , Ikeda M , Chi H , Lin C , Li G , Holman K , Tsuda T , Mar L , Foncin JF , Bruni AC , Montesi MP , Sorbi S , Rainero I , Pinessi L , Nee L , Chumakov I , Pollen D , Brookes A , Sanseau P , Polinsky RJ , Wasco W , Da Silva HA , Haines JL , Perkicak-Vance MA , Tanzi RE , Roses AD , Fraser PE , Rommens JM , St George-Hyslop PH ((1995) ) Cloning of a gene bearing missense mutations in early-onset familial Alzheimer’s disease. Nature 375: , 754–760. |
[13] | Cacace R , Sleegers K , Van Broeckhoven C ((2016) ) Molecular genetics of early-onset Alzheimer’s disease revisited. Alzheimers Dement 12: , 733–748. |
[14] | Jansen IE , Savage JE , Watanabe K , Bryois J , Williams DM , Steinberg S , Sealock J , Karlsson IK , Hägg S , Athanasiu L , Voyle N , Proitsi P , Witoelar A , Stringer S , Aarsland D , Almdahl IS , Andersen F , Bergh S , Bettella F , Bjornsson S , Brækhus A , Bråthen G , de Leeuw C , Desikan RS , Djurovic S , Dumitrescu L , Fladby T , Hohman TJ , Jonsson PV , Kiddle SJ , Rongve A , Saltvedt I , Sando SB , Selbæk G , Shoai M , Skene NG , Snaedal J , Stordal E , Ulstein ID , Wang Y , White LR , Hardy J , Hjerling-Leffler J , Sullivan PF , van der Flier WM , Dobson R , Davis LK , Stefansson H , Stefansson K , Pedersen NL , Ripke S , Andreassen OA , Posthuma D ((2019) ) Genome-wide meta-analysis identifies new loci and functional pathways influencing Alzheimer’s disease risk. Nat Genet 51: , 404–413. |
[15] | Lambert JC , Ibrahim-Verbaas CA , Harold D , Naj AC , Sims R , Bellenguez C , DeStafano AL , Bis JC , Beecham GW , Grenier-Boley B , et al. ((2013) ) Meta-analysis of 74,046 individuals identifies 11 new susceptibility loci for Alzheimer’s disease. Nat Genet 45: , 1452–1458. |
[16] | Braskie MN , Ringman JM , Thompson PM ((2011) ) Neuroimaging measures as endophenotypes in Alzheimer’s disease. Int J Alzheimers Dis 2011: , 490140. |
[17] | de Geus EJ ((2010) ) From genotype to EEG endophenotype: A route for post-genomic understanding of complex psychiatric disease? Genome Med 2: , 63. |
[18] | Maturana-Candelas A , Gómez C , Poza J , Pinto N , Hornero R ((2019) ) EEG characterization of the Alzheimer’s disease continuum by means of multiscale entropies. Entropy (Basel) 21: , 544. |
[19] | Núñez P , Poza J , Gómez C , Barroso-García V , Maturana-Candelas A , Tola-Arribas MA , Cano M , Hornero R ((2020) ) Characterization of the dynamic behavior of neural activity in Alzheimer’s disease: Exploring the non-stationarity and recurrence structure of EEG resting-state activity. J Neural Eng 17: , 016071. |
[20] | Ruiz-Gómez SJ , Hornero R , Poza J , Maturana-Candelas A , Pinto N , Gómez C ((2019) ) Computational modeling of the effects of EEG volume conduction on functional connectivity metrics. Application to Alzheimer’s disease continuum. J Neural Eng 16: , 066019. |
[21] | Dauwels J , Vialatte F , Cichocki A ((2010) ) Diagnosis of Alzheimer’s disease from EEG signals: Where are we standing? Curr Alzheimer Res 7: , 487–505. |
[22] | Jeong J ((2004) ) EEG dynamics in patients with Alzheimer’s disease. Clin Neurophysiology 115: , 1490–1505. |
[23] | Poza J , Gómez C , García M , Tola-Arribas MA , Carreres A , Cano M , Hornero R ((2017) ) Spatio-temporal fluctuations of neural dynamics in mild cognitive impairment and Alzheimer’s disease. Curr Alzheimer Res 14: , 924–936. |
[24] | Moretti DV , Miniussi C , Frisoni GB , Geroldi C , Zanetti O , Binetti G , Rossini PM ((2007) ) Hippocampal atrophy and EEG markers in subjects with mild cognitive impairment. Clin Neurophysiol 118: , 2716–2729. |
[25] | Fernandez-Mendoza J , Li Y , Vgontzas AN , Fang J , Gaines J , Calhoun SL , Liao D , Bixler EO ((2016) ) Insomnia is associated with cortical hyperarousal as early as adolescence. Sleep 39: , 1029–1036. |
[26] | Jaworska N , Berrigan L , Ahmed AG , Gray J , Bradford J , Korovessis A , Fedoroff P , Knott V ((2012) ) Resting electrocortical activity in adults with dysfunctional anger: A pilot study. Psychiatry Res 203: , 229–236. |
[27] | Cassani R , Estarellas M , San-Martin R , Fraga FJ , Falk TH ((2018) ) Systematic review on resting-state EEG for Alzheimer’s disease diagnosis and progression assessment. Dis Markers 2018: , 5174815. |
[28] | Horvath A , Szucs A , Csukly G , Sakovics A , Stefanics G , Kamondi A ((2018) ) EEG and ERP biomarkers of Alzheimer’s disease: A critical review. Front Biosci (Landmark Ed) 23: , 183–220. |
[29] | Rossini P , Iorio R , Vecchio F , DuBois B ((2020) ) Early diagnosis of Alzheimer’s disease: The role of biomarkers including advanced EEG signal analysis. Report from the IFCN-sponsored panel of experts. Clin Neurophysiol 131: , 1287–1310. |
[30] | Yang S , Bornot JMS , Wong-Lin K , Prasad G ((2019) ) M/EEG-based bio-markers to predict the MCI and Alzheimer’s disease: A review from the ML perspective. IEEE Trans Biomed Eng 66: , 2924–2935. |
[31] | de Geus EJC ((2002) ) Introducing genetic psychophysiology. Biol Psychol 61: , 1–10. |
[32] | Smit DJ , Boersma M , van Beijsterveldt CE , Posthuma D , Boomsma DI , Stam CJ , de Geus EJ ((2010) ) Endophenotypes in a dynamically connected brain. Behav Genet 40: , 167–177. |
[33] | Smit DJ , Posthuma D , Boomsma DI , Geus EJ ((2005) ) Heritability of background EEG across the power spectrum. Psychophysiology 42: , 691–697. |
[34] | Zietsch BP , Hansen JL , Hansell NK , Geffen GM , Martin NG , Wright MJ ((2007) ) Common and specific genetic influences on EEG power bands delta, theta, alpha, and beta. Biol Psychol 75: , 154–164. |
[35] | de Waal H , Stam CJ , de Haan W , van Straaten EC , Blankenstein MA , Scheltens P , van der Flier WM ((2013) ) Alzheimer’s disease patients not carrying the apolipoprotein E ɛ4 allele show more severe slowing of oscillatory brain activity. Neurobiol Aging 34: , 2158–2163. |
[36] | Ponomareva N , Andreeva T , Protasova M , Shagam L , Goltsov A , Malina D , Fokin V , Mitrofanov A , Rogaev E ((2013) ) Age-dependent effect of Alzheimer’s risk variant of CLU on EEG alpha rhythm in non-demented adults. Front Aging Neurosci 5: , 86. |
[37] | Ponomareva NV , Andreeva TV , Protasova MS , Shagam LI , Malina DD , Goltsov AY , Fokin VF , Illarioshkin SN , Rogaev EI ((2017) ) Quantitative EEG during normal aging: Association with the Alzheimer’s disease genetic risk variant in PICALM gene. Neurobiol Aging 51: , 177.e171–177.e178. |
[38] | Ponomareva NV , Korovaitseva GI , Rogaev EI ((2008) ) EEG alterations in non-demented individuals related to apolipoprotein E genotype and to risk of Alzheimer disease. Neurobiol Aging 29: , 819–827. |
[39] | Jack CR Jr. , Bennett DA , Blennow K , Carrillo MC , Dunn B , Haeberlein SB , Holtzman DM , Jagust W , Jessen F , Karlawish J , Liu E , Molinuevo JL , Montine T , Phelps C , Rankin KP , Rowe CC , Scheltens P , Siemers E , Snyder HM , Sperling R , Contributors ((2018) ) NIA-AA Research Framework: Toward a biological definition of Alzheimer’s disease. Alzheimers Dement 14: , 535–562. |
[40] | Folstein MF , Folstein SE , McHugh PR ((1975) ) “Mini-mental state”. A practical method for grading the cognitive state of patients for the clinician. J Psychiatr Res 12: , 189–198. |
[41] | Anderson CA , Pettersson FH , Clarke GM , Cardon LR , Morris AP , Zondervan KT ((2010) ) Data quality control in genetic case-control association studies. Nat Protoc 5: , 1564–1573. |
[42] | Purcell S , Neale B , Todd-Brown K , Thomas L , Ferreira MA , Bender D , Maller J , Sklar P , de Bakker PI , Daly MJ , Sham PC ((2007) ) PLINK: A tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: , 559–575. |
[43] | Patterson N , Price AL , Reich D ((2006) ) Population structure and eigenanalysis. PLoS Genet 2: , e190. |
[44] | Price AL , Patterson NJ , Plenge RM , Weinblatt ME , Shadick NA , Reich D ((2006) ) Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet 38: , 904–909. |
[45] | Genomes Project C , Auton A , Brooks LD , Durbin RM , Garrison EP , Kang HM , Korbel JO , Marchini JL , McCarthy S , McVean GA , Abecasis GR ((2015) ) A global reference for human genetic variation. Nature 526: , 68–74. |
[46] | Arias-Vásquez A , Isaacs A , Aulchenko YS , Hofman A , Oostra BA , Breteler M , van Duijn CM ((2007) ) The cholesteryl ester transfer protein (CETP) gene and the risk of Alzheimer’s disease. Neurogenetics 8: , 189–193. |
[47] | Hooper C , Killick R , Lovestone S ((2008) ) The GSK3 hypothesis of Alzheimer’s disease. J Neurochem 104: , 1433–1439. |
[48] | Ma X-Y , Yu J-T , Wang W , Wang H-F , Liu Q-Y , Zhang W , Tan L ((2013) ) Association of TOMM40 polymorphisms with late-onset Alzheimer’s disease in a Northern Han Chinese population. Neuromolecular Med 15: , 279–287. |
[49] | Murphy EA , Roddey JC , McEvoy LK , Holland D , Hagler DJ Jr. , Dale AM , Brewer JB , Alzheimer’s Disease Neuroimaging Initiative ((2012) ) CETP polymorphisms associate with brain structure, atrophy rate, and Alzheimer’s disease risk in an APOE-dependent manner. Brain Imaging Behav 6: , 16–26. |
[50] | Fehér Á , Juhász A , Pákáski M , Kálmán J , Janka Z ((2014) ) Association between the 9 repeat allele of the dopamine transporter 40bp variable tandem repeat polymorphism and Alzheimer’s disease. Psychiatry Res 220: , 730–731. |
[51] | Lyall DM , Harris SE , Bastin ME , Muñoz Maniega S , Murray C , Lutz MW , Saunders AM , Roses AD , Valdés Hernández Mdel C , Royle NA , Starr JM , Porteous DJ , Wardlaw JM , Deary IJ ((2014) ) Alzheimer’s disease susceptibility genes APOE and TOMM40, and brain white matter integrity in the Lothian Birth Cohort 1936. Neurobiol Aging 35: , 1513.e25–33. |
[52] | Pérez-Palma E , Bustos BI , Villamán CF , Alarcón MA , Avila ME , Ugarte GD , Reyes AE , Opazo C , De Ferrari GV , Alzheimer’s Disease Neuroimaging Initiative; NIA-LOAD/NCRAD Family Study Group ((2014) ) Overrepresentation of glutamate signaling in Alzheimer’s disease: Network-based pathway enrichment using meta-analysis of genome-wide association studies. PloS One 9: , e95413–e95413. |
[53] | Wetzel-Smith MK , Hunkapiller J , Bhangale TR , Srinivasan K , Maloney JA , Atwal JK , Sa SM , Yaylaoglu MB , Foreman O , Ortmann W , Rathore N , Hansen DV , Tessier-Lavigne M , Mayeux R , Pericak-Vance M , Haines J , Farrer LA , Schellenberg GD , Goate A , Behrens TW , Cruchaga C , Watts RJ , Graham RR ((2014) ) A rare mutation in UNC5C predisposes to late-onset Alzheimer’s disease and increases neuronal cell death. Nat Med 20: , 1452–1457. |
[54] | Khondoker M , Newhouse S , Westman E , Muehlboeck JS , Mecocci P , Vellas B , Tsolaki M , Kłoszewska I , Soininen H , Lovestone S , Dobson R , Simmons A ((2015) ) Linking genetics of brain changes to Alzheimer’s disease: Sparse whole genome association scan of regional MRI volumes in the ADNI and AddNeuroMed Cohorts. J Alzheimers Dis 45: , 851–864. |
[55] | Roussotte FF , Gutman BA , Hibar DP , Madsen SK , Narr KL , Thompson PM ((2015) ) Carriers of a common variant in the dopamine transporter gene have greater dementia risk, cognitive decline, and faster ventricular expansion. Alzheimers Dement 11: , 1153–1162. |
[56] | Yan J , Kim S , Nho K , Chen R , Risacher SL , Moore JH , Saykin AJ , Shen L , Alzheimer’s Disease Neuroimaging I ((2015) ) Hippocampal transcriptome-guided genetic analysis of correlated episodic memory phenotypes in Alzheimer’s disease. Front Genet 6: , 117. |
[57] | Darvesh S ((2016) ) Butyrylcholinesterase as a diagnostic and therapeutic target for Alzheimer’s disease. Curr Alzheimer Res 13: , 1173–1177. |
[58] | Hu C , Chen W , Myers SJ , Yuan H , Traynelis SF ((2016) ) Human GRIN2B variants in neurodevelopmental disorders. J Pharmacol Sci 132: , 115–121. |
[59] | Ji H , Wang Y , Liu G , Chang L , Chen Z , Zhou D , Xu X , Cui W , Hong Q , Jiang L , Li J , Zhou X , Li Y , Guo Z , Zha Q , Niu Y , Weng Q , Duan S , Wang Q ((2017) ) Elevated OPRD1 promoter methylation in Alzheimer’s disease patients. PLoS One 12: , e0172335. |
[60] | Sun JH , Wang HF , Zhu XC , Yu WJ , Tan CC , Jiang T , Tan MS , Tan L , Yu JT ((2016) ) The impact of UNC5C genetic variations on neuroimaging in Alzheimer’s disease. Mol Neurobiol 53: , 6759–6767. |
[61] | Macdonald IR , Maxwell SP , Reid GA , Cash MK , DeBay DR , Darvesh S ((2017) ) Quantification of butyrylcholinesterase activity as a sensitive and specific biomarker of Alzheimer’s disease. J Alzheimers Dis 58: , 491–505. |
[62] | Zettergren A , Höglund K , Kern S , Thorvaldsson V , Johan Skoog M , Hansson O , Andreasen N , Bogdanovic N , Blennow K , Skoog I , Zetterberg H ((2019) ) Association of IL1RAP-related genetic variation with cerebrospinal fluid concentration of Alzheimer-associated tau protein. Sci Rep 9: , 2460–2460. |
[63] | Andreoli V , De Marco EV , Trecroci F , Cittadella R , Di Palma G , Gambardella A ((2014) ) Potential involvement of GRIN2B encoding the NMDA receptor subunit NR2B in the spectrum of Alzheimer’s disease. J Neural Transm (Vienna) 121: , 533–542. |
[64] | Ramanan VK , Risacher SL , Nho K , Kim S , Swaminathan S , Shen L , Foroud TM , Hakonarson H , Huentelman MJ , Aisen PS , Petersen RC , Green RC , Jack CR , Koeppe RA , Jagust WJ , Weiner MW , Saykin AJ ((2014) ) APOE and BCHE as modulators of cerebral amyloid deposition: A florbetapir PET genome-wide association study. Mol Psychiatry 19: , 351–357. |
[65] | Warstadt NM , Dennis EL , Jahanshad N , Kohannim O , Nir TM , McMahon KL , de Zubicaray GI , Montgomery GW , Henders AK , Martin NG , Whitfield JB , Jack CR Jr. , Bernstein MA , Weiner MW , Toga AW , Wright MJ , Thompson PM ((2014) ) Serum cholesterol and variant in cholesterol-related gene CETP predict white matter microstructure. Neurobiol Aging 35: , 2504–2513. |
[66] | Haddley K , Vasiliou AS , Ali FR , Paredes UM , Bubb VJ , Quinn JP ((2008) ) Molecular genetics of monoamine transporters: Relevance to brain disorders. Neurochem Res 33: , 652–667. |
[67] | Hohman TJ , Koran MEI , Thornton-Wells TA ((2014) ) Interactions between GSK3β and amyloid genes explain variance in amyloid burden. Neurobiol Aging 35: , 460–465. |
[68] | Roussotte FF , Jahanshad N , Hibar DP , Sowell ER , Kohannim O , Barysheva M , Hansell NK , McMahon KL , de Zubicaray GI , Montgomery GW , Martin NG , Wright MJ , Toga AW , Jack CR Jr. , Weiner MW , Thompson PM ((2014) ) A commonly carried genetic variant in the delta opioid receptor gene, OPRD1, is associated with smaller regional brain volumes: Replication in elderly and young populations. Hum Brain Mapp 35: , 1226–1236. |
[69] | Ramanan VK , Risacher SL , Nho K , Kim S , Shen L , McDonald BC , Yoder KK , Hutchins GD , West JD , Tallman EF , Gao S , Foroud TM , Farlow MR , De Jager PL , Bennett DA , Aisen PS , Petersen RC , Jack CR Jr. , Toga AW , Green RC , Jagust WJ , Weiner MW , Saykin AJ ((2015) ) GWAS of longitudinal amyloid accumulation on 18F-florbetapir PET in Alzheimer’s disease implicates microglial activation gene IL1RAP. Brain 138: , 3076–3088. |
[70] | Núñez P , Poza J , Bachiller A , Gomez-Pilar J , Lubeiro A , Molina V , Hornero R ((2017) ) Exploring non-stationarity patterns in schizophrenia: Neural reorganization abnormalities in the alpha band. J Neural Eng 14: , 046001. |
[71] | Ruiz-Gómez SJ , Gómez C , Poza J , Martínez-Zarzuela M , Tola-Arribas MA , Cano M , Hornero R ((2018) ) Measuring alterations of spontaneous EEG neural coupling in Alzheimer’s disease and mild cognitive impairment by means of cross-entropy metrics. Front Neuroinform 12: , 76. |
[72] | Whitham EM , Lewis T , Pope KJ , Fitzgibbon SP , Clark CR , Loveless S , DeLosAngeles D , Wallace AK , Broberg M , Willoughby JO ((2008) ) Thinking activates EMG in scalp electrical recordings. Clin Neurophysiol 119: , 1166–1175. |
[73] | Whitham EM , Pope KJ , Fitzgibbon SP , Lewis T , Clark CR , Loveless S , Broberg M , Wallace A , DeLosAngeles D , Lillie P , Hardy A , Fronsko R , Pulbrook A , Willoughby JO ((2007) ) Scalp electrical recording during paralysis: Quantitative evidence that EEG frequencies above 20Hz are contaminated by EMG. Clin Neurophysiol 118: , 1877–1888. |
[74] | Nunez P , Poza J , Gomez C , Rodriguez-Gonzalez V , Hillebrand A , Tola-Arribas MA , Cano M , Hornero R ((2019) ) Characterizing the fluctuations of dynamic resting-state electrophysiological functional connectivity: Reduced neuronal coupling variability in mild cognitive impairment and dementia due to Alzheimer’s disease. J Neural Eng 16: , 056030. |
[75] | Visscher PM , Wray NR , Zhang Q , Sklar P , McCarthy MI , Brown MA , Yang J ((2017) ) 10 years of GWAS discovery: Biology, function, and translation. Am J Hum Genet 101: , 5–22. |
[76] | Shim H , Chasman DI , Smith JD , Mora S , Ridker PM , Nickerson DA , Krauss RM , Stephens M ((2015) ) A multivariate genome-wide association analysis of 10 LDL subfractions, and their response to statin treatment, in 1868 caucasians. PLoS One 10: , e0120758. |
[77] | R Development Core Team (2010) R Foundation for Statistical Computing, Vienna Austria. |
[78] | Machiela MJ , Chanock SJ ((2015) ) LDlink: A web-based application for exploring population-specific haplotype structure and linking correlated alleles of possible functional variants. Bioinformatics 31: , 3555–3557. |
[79] | Teslovich TM , Musunuru K , Smith AV , Edmondson AC , Stylianou IM , Koseki M , Pirruccello JP , Ripatti S , Chasman DI , Willer CJ , et al. ((2010) ) Biological, clinical and population relevance of 95 loci for blood lipids. Nature 466: , 707–713. |
[80] | Gomez C , Perez-Macias JM , Poza J , Fernandez A , Hornero R ((2013) ) Spectral changes in spontaneous MEG activity across the lifespan. J Neural Eng 10: , 066006. |
[81] | Bennys K , Rondouin G , Vergnes C , Touchon J ((2001) ) Diagnostic value of quantitative EEG in Alzheimer’s disease. Neurophysiol Clin 31: , 153–160. |
[82] | Huang C , Wahlund L , Dierks T , Julin P , Winblad B , Jelic V ((2000) ) Discrimination of Alzheimer’s disease and mild cognitive impairment by equivalent EEG sources: A cross-sectional and longitudinal study. Clin Neurophysiol 111: , 1961–1967. |
[83] | Musaeus CS , Engedal K , Høgh P , Jelic V , Mørup M , Naik M , Oeksengaard AR , Snaedal J , Wahlund LO , Waldemar G , Andersen BB ((2018) ) EEG theta power is an early marker of cognitive decline in dementia due to Alzheimer’s disease. J Alzheimers Dis 64: , 1359–1371. |
[84] | Caravaglios G , Castro G , Costanzo E , Di Maria G , Mancuso D , Muscoso EG ((2010) ) Theta power responses in mild Alzheimer’s disease during an auditory oddball paradigm: Lack of theta enhancement during stimulus processing. J Neural Transm 117: , 1195–1208. |
[85] | Klimesch W , Doppelmayr M , Yonelinas A , Kroll NEA , Lazzara M , Rohm D , Gruber W ((2001) ) Theta synchronization during episodic retrieval: Neural correlates of conscious awareness. Cogn Brain Res 12: , 33–38. |
[86] | Klimesch W , Vogt F , Doppelmayr M ((1999) ) Interindividual differences in alpha and theta power reflect memory performance. Intelligence 27: , 347–362. |
[87] | Grunwald M , Busse F , Hensel A , Kruggel T , Riedel-Heller S , Wolf M , Arendt T , Gertz HJ ((2001) ) Correlation between cortical theta activity and hippocampal volumes in health, mild cognitive impairment, and mild dementia. J Clin Neurophysiol 18: , 178–184. |
[88] | Grunwald M , Hensel A , Wolf H , Weiss T , Gertz HJ ((2007) ) Does the hippocampal atrophy correlate with the cortical theta power in elderly subjects with a range of cognitive impairment? J Clin Neurophysiol 24: , 22–26. |
[89] | Diehl B , Najm I , Mohamed A , Babb T , Ying Z , Bingaman W ((2002) ) Interictal EEG, hippocampal atrophy, and cell densities in hippocampal sclerosis and hippocampal sclerosis associated with microscopic cortical dysplasia. J Clin Neurophysiol 19: , 157–162. |
[90] | Fonseca LC , Tedrus GMAS , Prandi LR , Almeida AM , Furlanetto DS ((2011) ) Alzheimer’s disease: Relationship between cognitive aspects and power and coherence EEG measures. Arq Neuropsiquiatr 69: , 875–881. |
[91] | Koelewijn L , Lancaster TM , Linden D , Dima DC , Routley BC , Magazzini L , Barawi K , Brindley L , Adams R , Tansey KE , Bompas A , Tales A , Bayer A , Singh K ((2019) ) Oscillatory hyperactivity and hyperconnectivity in young APOE-ɛ4 carriers and hypoconnectivity in Alzheimer’s disease. Elife 8: , e36011. |
[92] | Scott L , Feng J , Kiss T , Needle E , Atchison K , Kawabe TT , Milici AJ , Hajos-Korcsok E , Riddell D , Hajos M ((2012) ) Age-dependent disruption in hippocampal theta oscillation in amyloid-beta overproducing transgenic mice. Neurobiol Aging 33: , 1481.e13–23. |
[93] | Stomrud E , Hansson O , Minthon L , Blennow K , Rosen I , Londos E ((2010) ) Slowing of EEG correlates with CSF biomarkers and reduced cognitive speed in elderly with normal cognition over 4 years. Neurobiol Aging 31: , 215–223. |
[94] | Du XM , Rowland LM , Summerfelt A , Choa FS , Wittenberg GF , Wisner K , Wijtenburg A , Chiappelli J , Kochunov P , Hong LE ((2018) ) Cerebellar-stimulation evoked prefrontal electrical synchrony is modulated by GABA. Cerebellum 17: , 550–563. |
[95] | Casula EP , Pellicciari MC , Ponzo V , Bassi MS , Veniero D , Caltagirone C , Koch G ((2016) ) Cerebellar theta burst stimulation modulates the neural activity of interconnected parietal and motor areas. Sci Rep 6: , 36191. |
[96] | Casula EP , Pellicciari MC , Picazio S , Caltagirone C , Koch G ((2016) ) Spike-timing-dependent plasticity in the human dorso-lateral prefrontal cortex. Neuroimage 143: , 204–213. |
[97] | Zhang H , Jacobs J ((2015) ) Traveling theta waves in the human hippocampus. J Neurosci 35: , 12477–12487. |
[98] | Gabay C , Lamacchia C , Palmer G ((2010) ) IL-1 pathways in inflammation and human diseases. Nat Rev Rheumatol 6: , 232–241. |
[99] | Goshen I , Yirmiya R ((2009) ) Interleukin-1 (IL-1): A central regulator of stress responses. Front Neuroendocrinol 30: , 30–45. |
[100] | Krueger JM ((2008) ) The role of cytokines in sleep regulation. Curr Pharm Des 14: , 3408–3416. |
[101] | Matsuki T , Horai R , Sudo K , Iwakura Y ((2003) ) IL-1 plays an important role in lipid metabolism by regulating insulin levels under physiological conditions. J Exp Med 198: , 877–888. |
[102] | Rothwell NJ , Luheshi GN ((2000) ) Interleukin 1 in the brain Biology,athology and therapeutic target. Trends Neurosci 23: , 618–625 |
[103] | Smith DE , Lipsky BP , Russell C , Ketchem RR , Kirchner J , Hensley K , Huang Y , Friedman WJ , Boissonneault V , Plante MM , Rivest S , Sims JE ((2009) ) A central nervous system-restricted isoform of the interleukin-1 receptor accessory protein modulates neuronal responses to interleukin-1. Immunity 30: , 817–831. |
[104] | Yoshida T , Shiroshima T , Lee SJ , Yasumura M , Uemura T , Chen X , Iwakura Y , Mishina M ((2012) ) Interleukin-1 receptor accessory protein organizes neuronal synaptogenesis as a cell adhesion molecule. J Neurosci 32: , 2588–2600. |
[105] | Ghosh S , Wu MD , Shaftel SS , Kyrkanides S , LaFerla FM , Olschowka JA , O’Banion MK ((2013) ) Sustained interleukin-1β overexpression exacerbates tau pathology despite reduced amyloid burden in an Alzheimer’s mouse model. J Neurosci 33: , 5053–5064. |
[106] | Anderson JM , Hampton DW , Patani R , Pryce G , Crowther RA , Reynolds R , Franklin RJ , Giovannoni G , Compston DA , Baker D , Spillantini MG , Chandran S ((2008) ) Abnormally phosphorylated tau is associated with neuronal and axonal loss in experimental autoimmune encephalomyelitis and multiple sclerosis. Brain 131: , 1736–1748. |
[107] | Mondragón-Rodríguez S , Salas-Gallardo A , González-Pereyra P , Macías M , Ordaz B , Peña-Ortega F , Aguilar-Vázquez A , Orta-Salazar E , Díaz-Cintra S , Perry G , Williams S ((2018) ) Phosphorylation of Tau protein correlates with changes in hippocampal theta oscillations and reduces hippocampal excitability in Alzheimer’s model. J Biol Chem 293: , 8462–8472. |
[108] | Assenza G , Pellegrino G , Tombini M , Di Pino G , Di Lazzaro V ((2015) ) Wakefulness delta waves increase after cortical plasticity induction. Clin Neurophysiol 126: , 1221–1227. |
[109] | Jiao B , Liu X , Tang B , Hou L , Zhou L , Zhang F , Zhou Y , Guo J , Yan X , Shen L ((2014) ) Investigation of TREM2, PLD3, and UNC5C variants in patients with Alzheimer’s disease from mainland China. Neurobiol Aging 35: , 2422 e2429–2422 e2411. |
[110] | Korvatska O , Leverenz JB , Jayadev S , McMillan P , Kurtz I , Guo X , Rumbaugh M , Matsushita M , Girirajan S , Dorschner MO , Kiianitsa K , Yu C-E , Brkanac Z , Garden GA , Raskind WH , Bird TD ((2015) ) R47H variant of TREM2 associated with Alzheimer disease in a large late-onset family: Clinical, genetic, and neuropathological study. JAMA Neurol 72: , 920–927. |
[111] | Hashimoto Y , Toyama Y , Kusakari S , Nawa M , Matsuoka M ((2016) ) An Alzheimer disease-linked rare mutation potentiates netrin receptor uncoordinated-5C-induced signaling that merges with amyloid β precursor protein signaling. J Biol Chem 291: , 12282–12293. |
[112] | Küry S , Garrec C , Airaud F , Breheret F , Guibert V , Frenard C , Jiao S , Bonneau D , Berthet P , Bossard C , Ingster O , Cauchin E , Bezieau S ((2014) ) Evaluation of the colorectal cancer risk conferred by rare UNC5C alleles. World J Gastroenterol 20: , 204–213. |
[113] | Kim D , Ackerman SL ((2011) ) The UNC5C netrin receptor regulates dorsal guidance of mouse hindbrain axons. J Neurosci 31: , 2167–2179. |
[114] | Mehlen P , Guenebeaud C ((2010) ) Netrin-1 and its dependence receptors as original targets for cancer therapy. Curr Opin Oncol 22: , 46–54. |
[115] | Poliak S , Morales D , Croteau LP , Krawchuk D , Palmesino E , Morton S , Cloutier JF , Charron F , Dalva MB , Ackerman SL , Kao TJ , Kania A ((2015) ) Synergistic integration of Netrin and ephrin axon guidance signals by spinal motor neurons. Elife 4: , e10841. |
[116] | Babiloni C , Frisoni GB , Pievani M , Vecchio F , Infarinato F , Geroldi C , Salinari S , Ferri R , Fracassi C , Eusebi F , Rossini PM ((2008) ) White matter vascular lesions are related to parietal-to-frontal coupling of EEG rhythms in mild cognitive impairment. Hum Brain Mapp 29: , 1355–1367. |
[117] | Smit DJ , Boomsma DI , Schnack HG , Hulshoff Pol HE , de Geus EJ ((2012) ) Individual differences in EEG spectral power reflect genetic variance in gray and white matter volumes. Twin Res Hum Genet 15: , 384–392. |
[118] | Maes T , Barcelo A , Buesa C ((2002) ) Neuron navigator: A human gene family with homology to unc-53, a cell guidance gene from Caenorhabditis elegans. Genomics 80: , 21–30. |
[119] | Pook C , Ahrens JM , Clagett-Dame M ((2020) ) Expression pattern of Nav2 in the murine CNS with development. Gene Expr Patterns 35: , 119099. |
[120] | McNeill EM , Roos KP , Moechars D , Clagett-Dame M ((2010) ) Nav2 is necessary for cranial nerve development and blood pressure regulation. Neural Dev 5: , 6. |
[121] | Wang KS , Liu Y , Xu C , Liu X , Luo X ((2017) ) Family-based association analysis of NAV2 gene with the risk and age at onset of Alzheimer’s disease. J Neuroimmunol 310: , 60–65. |
[122] | Klimesch W , Schimke H , Schwaiger J ((1994) ) Episodic and semantic memory: An analysis in the EEG theta and alpha band. Electroencephalogr Clin Neurophysiol 91: , 428–441. |
[123] | Nyhus E , Curran T ((2010) ) Functional role of gamma and theta oscillations in episodic memory. Neurosci Biobehav Rev 34: , 1023–1035. |
[124] | Staudigl T , Hanslmayr S ((2013) ) Theta oscillations at encoding mediate the context-dependent nature of human episodic memory. Curr Biol 23: , 1101–1106. |
[125] | Escudero J , Abasolo D , Hornero R , Espino P , Lopez M ((2006) ) Analysis of electroencephalograms in Alzheimer’s disease patients with multiscale entropy. Physiol Meas 27: , 1091–1106. |
[126] | Fraga FJ , Falk TH , Kanda PAM , Anghinah R ((2013) ) Characterizing Alzheimer’s disease severity via resting-awake EEG amplitude modulation analysis. PloS One 8: , e72240–e72240. |
[127] | Locatelli T , Cursi M , Liberati D , Franceschi M , Comi G ((1998) ) EEG coherence in Alzheimer’s disease. Electroencephalogr Clin Neurophysiol 106: , 229–237. |
[128] | Wang R , Wang J , Li S , Yu H , Deng B , Wei X ((2015) ) Multiple feature extraction and classification of electroencephalograph signal for Alzheimers’ with spectrum and bispectrum. Chaos 25: , 013110. |
