Next Article in Journal
Visual-Attentional Load Unveils Slowed Processing Speed in Multiple Sclerosis Patients: A Pilot Study with a Tablet-Based Videogame
Next Article in Special Issue
Speech–Brain Frequency Entrainment of Dyslexia with and without Phonological Deficits
Previous Article in Journal
Leukoencephalopathy with Calcifications and Cysts—The First Polish Patient with Labrune Syndrome
Previous Article in Special Issue
Quantitative Identification of Functional Connectivity Disturbances in Neuropsychiatric Lupus Based on Resting-State fMRI: A Robust Machine Learning Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Scan for Five Brain Oscillatory Phenotypes Identifies a New QTL Associated with Theta EEG Band

1
IPATIMUP—Instituto de Patologia e Imunologia Molecular da Universidade do Porto, 4200-135 Porto, Portugal
2
I3S—Instituto de Investigação e Inovação em Saúde, Universidade do Porto, 4200-135 Porto, Portugal
3
Grupo de Ingeniería Biomédica, Universidad de Valladolid, 47011 Valladolid, Spain
4
Centro de Investigación Biomédica en Red en Bioingeniería, Biomateriales y Nanomedicina, (CIBER-BBN), 47011 Valladolid, Spain
5
Instituto de Investigación en Matemáticas (IMUVA), Universidad de Valladolid, 47011 Valladolid, Spain
6
Associação Portuguesa de Familiares e Amigos de Doentes de Alzheimer, Delegação Norte, 4455-301 Lavra, Portugal
7
Asociación de Familiares y Amigos de Enfermos de Alzheimer y otras demencias de Zamora, 49021 Zamora, Spain
8
Department of Biochemistry, Genetics and Immunology, University of Vigo, 36310 Vigo, Spain
9
TellmeGen, 46010 Valencia, Spain
10
Centro de Matemática da, Universidade do Porto, 4169-007 Porto, Portugal
*
Authors to whom correspondence should be addressed.
Brain Sci. 2020, 10(11), 870; https://doi.org/10.3390/brainsci10110870
Submission received: 15 October 2020 / Revised: 6 November 2020 / Accepted: 10 November 2020 / Published: 18 November 2020
(This article belongs to the Special Issue Human Brain Dynamics: Latest Advances and Prospects)

Abstract

:
Brain waves, measured by electroencephalography (EEG), are a powerful tool in the investigation of neurophysiological traits and a noninvasive and cost-effective alternative in the diagnostic of some neurological diseases. In order to identify novel Quantitative Trait Loci (QTLs) for brain wave relative power (RP), we collected resting state EEG data in five frequency bands (δ, θ, α, β1, and β2) and genome-wide data in a cohort of 105 patients with late onset Alzheimer’s disease (LOAD), 41 individuals with mild cognitive impairment and 45 controls from Iberia, correcting for disease status. One novel association was found with an interesting candidate for a role in brain wave biology, CLEC16A (C-type lectin domain family 16), with a variant at this locus passing the adjusted genome-wide significance threshold after Bonferroni correction. This finding reinforces the importance of immune regulation in brain function. Additionally, at a significance cutoff value of 5 × 10−6, 18 independent association signals were detected. These signals comprise brain expression Quantitative Loci (eQTLs) in caudate basal ganglia, spinal cord, anterior cingulate cortex and hypothalamus, as well as chromatin interactions in adult and fetal cortex, neural progenitor cells and hippocampus. Moreover, in the set of genes showing signals of association with brain wave RP in our dataset, there is an overrepresentation of loci previously associated with neurological traits and pathologies, evidencing the pleiotropy of the genetic variation modulating brain function.

Graphical Abstract

1. Introduction

Neurons, as excitable cells with electrical properties, generate coherent electric and magnetic fields when they are synchronously activated, being akin to current dipoles that may be recorded by electrodes at the scalp [1]. Electroencephalography (EEG) captures the summed electrical activities of neuron populations and EEG signals co-vary strongly with different levels of arousal and consciousness [1]. These signals can be viewed as reflections of outputs from the regulation of synchronization/desynchronization and functional coupling/decoupling in neuron populations with effect on vigilance, motivation and cognitive processes [2,3].
Brain waves recorded as electrophysiological signals harbor an impressive amount of information, thus offering a noninvasive and cost-effective alternative into the investigation of neurophysiological mechanisms. Neural networks change with development, age, disease and personal experience, and abnormalities in the oscillatory activity are associated to neurological disorders such as schizophrenia, autism spectrum disorder (ASD), attention deficit hyperactivity disorder (ADHD) and Alzheimer’s disease (AD) [4,5,6,7,8,9,10]. The mechanisms of synchronization/desynchronization of thalamocortical and ascending activity systems can be analyzed by measuring the amplitude or source activity of resting-state eyes-closed cortical EEG rhythms and can reveal the effects of AD on brain function [4]. In this regard, abnormally reduced spectral coherence in α and β rhythms has been associated with AD, as well as an increase of δ and θ rhythms, contributing with valuable information to build a prediction pipeline for the disease [7,8].
All current evidence suggests that oscillations are essential components in neural computation and not just biomarkers of the process. These manifestations are amongst the most heritable traits in humans [11,12], since individual differences across the EEG frequency spectra in adult population are largely determined by genetic factors. Generally, hereditability is highest around the α frequency band and lower in θ and δ bands [13]. These signals may act as an intermediate phenotype between genetics and behavior, representing endophenotypes [14,15].
The use of brain endophenotypes, such as EEG, has been proposed as a valid approach to uncover the contribution of a given variant to disease, even with modest sample sizes (an important advantage of quantitative phenotypes) [14], by reducing the inherent complexity of behavioral and neurological traits [16,17]. Indeed, although genome-wide association studies (GWAS) have played an important role in uncovering many genetic variants for some conditions, this method has not been so successful for psychiatric and neurological disorders [18]. This is probably due to a more complex genetic architecture of these phenotypes, which reflects the intricacy of the underlying brain processes. Nevertheless, as EEG patterns also likely come from a polygenic model of inheritance, alternative approaches to maximize power to detect relevant biological pathways without large sample sizes are needed. Gene-based and expression-based enrichment analyses increase the power to detect genes affecting the phenotype under study [19,20,21].
Only a few studies to date have searched for genetic variation associated to different brain oscillatory traits, using individuals affected by neurological diseases and controls, and the number of genes identified is still small [22]. We hypothesized that additional quantitative trait loci (QTLs) for brain waves at different frequency bands could be uncovered when using data from cohorts with different EEG phenotypic ranges. For this purpose, resting state EEG data from a novel cohort of 105 patients with late onset Alzheimer’s disease (LOAD), 41 individuals with mild cognitive impairment (MCI) and 45 controls from Iberia were analyzed by means of a genome-wide QTL analysis for brain wave relative power (RP) in five conventional EEG frequency bands (δ, θ, α, β1, and β2), corrected for disease status. With this strategy, we aimed to identify new genes associated with EEG endophenotypes. Our findings contribute a number of interesting candidate genes with a role in brain function.

2. Materials and Methods

2.1. Subjects

The current study included 105 LOAD patients from North Portugal (n = 42) and from the Spanish autonomous community of Castile and León (n = 63) with a clinical diagnosis of dementia due to LOAD. AD patients were diagnosed following the criteria of the National Institute on Aging and Alzheimer’s Association (NIA-AA) [23]. The severity of cognitive impairment in AD was assessed by Mini-Mental State Examination (MMSE) test. Patients were classified in the four main stages of AD development: mild (MIL), moderate (MOD), and severe Alzheimer’s disease (SEV).
In addition, 41 MCI subjects and 45 elderly controls with no history of neurological or major psychiatric disorders from both geographical regions were also included in the study. In total, 191 individuals were analyzed. The age of the subjects ranged between 63 and 97 years old, with a mean age of 79 years old for controls, 84 for individuals with MCI and 81 for LOAD patients.
This project has been approved by the Ethics Committee of the University of Porto (CEUP) (report # 38/CEUP/2018), and written informed consents were obtained from all participants, family or legal representatives.

2.2. EEG Recording and Processing

Resting state EEG was acquired using a 19-channel EEG system (Nihon Kohden Neurofax JE-921A, Tokyo, Japan) at a sampling frequency of 500 Hz. EEG signals were recorded at electrodes F3, F4, F7, F8, Fp1, Fp2, T3, T4, T5, T6, C3, C4, P3, P4, O1, O2, Fz, Cz and Pz of the International 10–20 System and then re-referenced by means of common average referencing. Subjects were asked to remain awake with eyes closed during EEG acquisition. Each five-minute EEG recording was processed by independent component analysis (ICA), digital filtering using a Hamming window bandpass finite impulse response (FIR) filter in the band of interest (1–30 Hz) and by selecting 5 s artifact-free epochs by visual inspection [24]. Frequency bins were defined as the conventional EEG frequency bands: δ (1–4 Hz), θ (4–8 Hz), α (8–13 Hz), β1 (13–19 Hz) and β2 (19–30 Hz). The metric used to describe the distribution of the spectral content of the signals and explore associations was RP. RP quantifies the relative contribution of each frequency band to the global power spectrum. It is calculated from the power spectral density (PSD) function by summing the contribution of each spectral component in a specific band [24].

2.3. Genotyping and Genome-Wide Association Analysis

Saliva was collected from participants either using Oragene DNA (OG-500) self-collection kits (DNA Genotek, Ottawa, Canada) or cotton sterile buccal swabs. DNA was then extracted from the liquid saliva samples using the prepIT DNA extraction kit (DNA Genotek, Ottawa, Canada) and using the Citogene extraction kit (Citomed, Odivelas, Portugal) for the buccal swabs following the manufacturer’s protocols. Samples were genotyped with the Axiom Spain Biobank Array (Thermo Fisher, Waltham, Massachusetts, USA), and the genotyping service was carried out at CEGEN-PRB3-ISCIII, Santigo de Compostela, Spain; supported by grant PT17/0019, of the PE I + D + I. Variant calling and quality control pre-analysis were performed using the Affymetrix Power Tools (APT) and PLINK [25]. Individuals with outlying missing genotype or heterozygosity rates (between μ ± 3σ), as well as individuals estimated as related through identical-by-descent measurements and markers with significant deviation to Hardy–Weinberg equilibrium (HWE, α = 10−7) were removed prior to the analysis. Through Principal Component Analysis (PCA, implemented in SMARTPCA from EIGENSOFT) individuals with divergent ancestry were identified and removed [26]. The final total genotyping rate was 0.996. Differences between populations from PCA were not statistically significant (p = 0.53).
Common variants were selected, considering minimum allele frequencies (MAF) above 5%, which resulted in a set of 388,620 variants. For a nominal significance level α = 0.05, this set of variants leads to a Bonferroni genome-wide significance threshold of 1.29 × 10−7, when considering each of the five frequency bands.
To detect association between variants and brain EEG signals, we implemented a multivariable linear regression model in PLINK. The regression was done with 5 continuous variables, corresponding to each of the 5 different EEG frequency bands selected for analysis (δ, θ, α, β1, and β2). The model included 5 covariates to correct for residual population substructure: the first 2 principal components (PC), age (since RP suffer important changes across the lifespan [27]), sex and disease status (encoded as 0: control; 1: MCI; 2: MIL; 3: MOD; 4; SEV):
Y = a + b1.ADD + b2.PC1 + b3.PC2 + b4.AGE + b5.SEX + b6.dstatus + ϵ
being b1 the coefficient for the applied additive genotypic model of association between SNP and EEG frequency band (ADD). PCs were computed using Eigensoft’s smartpca for the 191 individuals [28]. With this model, all SNP effects found in this study will mostly reflect the genetic effects on brain electrophysiology.
The cutoff value used to evaluate associations was the Bonferroni genome-wide significance threshold of 1.29 × 10−7, but a less stringent cutoff of 5 × 10−6 was also applied (since our sample size is small) for gene-based analysis.

2.4. Common Variant QTL and eQTL Enrichment Analysis

Considering the five linear regressions, top candidate variants were selected considering the significance level α = 0.01 (= 0.05/5). This resulted in the analysis of a set of 16,575 variants.
Functional annotation of GWAS results, comprising gene identification and prioritization by positional, expression-QTL (eQTL) and chromatin interaction (CI) mapping, as well as gene-based pathway and tissue enrichment, was carried out using FUMA (https://fuma.ctglab.nl/) to find likely causal relations from the summary statistics provided (rsID, p-value) [21]. SNPs were mapped to genes up to 10 kb apart and the r2 threshold of the linkage disequilibrium (LD) to define independent significant SNPs was set to 0.8 (≥). SNPs were filtered based on chromatin state in the brain (accessibility of genomic regions every 200 bp based on 15 categorical states) for state ≤ 7 (open in given tissue/cell type). Gene eQTLs were mapped to Genotype-Tissue Expression (GTEx) V8 Brain tissues filtered by raw p-values (α = 5 × 10−4) and chromatin state in the brain ≤ 7. Datasets of brain tissues and cell types were selected for 3D CI mapping—the significance of interaction threshold was set to false discovery rate (FDR) ≤ 10−6, as suggested by Schmitt and colleagues [29]. The promoter region window was set to 250 bp upstream and 500 bp downstream, to overlap transcription start site (TSS) of genes to significantly interacted regions with risk loci. Overlapping genes were used for gene mapping. Only SNPs overlapped with enhancers and promoters from brain epigenomes were mapped. More detailed information on the analyses incorporated in FUMA is available at https://fuma.ctglab.nl/.
The 50 genes with lower p-values were selected and tested for gene-based pathway and brain tissue enrichment, for enrichment in curated gene sets from the HGRI GWAS catalog of reported genes, as well as for evaluating protein interactions using STRING v.11.0 (https://string-db.org/) [30,31].

2.5. Statistical Power Computation

The computation of the statistical power to detect an association is of paramount importance. The power to detect a variant-trait association can be quantified through the non-centrality parameter (NPC), which is the expected value of a test statistic to detect association when the null hypothesis is true [32]:
NPC = n × R2/(1 − R2), with R2 = r2 × q2
where n is the sample size, q2 is the proportion of phenotypic variance explained by a specific causal variant in the population, and r2 is the squared LD correlation between the causal variant and the genotyped one. R 2 is the proportion of variance explained by the genotyped SNP in the population. If the genotypes at the causal locus are in HWE [33,34], then
q2 = 2 × MAF × (1 − MAF) × β2
where β is the effect size of an allele on the phenotype. This assumes that the analysis for detecting an association is computed by regression of the phenotype on the genotype count (zero, one, or two minor alleles). R2 can then be computed as the ratio between q2 and the phenotype’s total variance [35]:
R2 = q2/var(Y).

3. Results

3.1. Linear Regression for 5 Oscillatory Phenotypes Identified One New Association

We tested a total of 388,620 common variants (MAF > 5%). The linear regressions for each of the five frequency bands revealed some deviations from the null distribution (Figure 1). This is more clearly noticed for the θ, β1 and β2 frequency bands. At a significance level α = 5 × 10−6, 19 highly correlated SNPs were identified, corresponding to 18 independent signals (since two of the loci are in LD), namely 1 in δ, 8 in θ, 2 in α, 4 in β1 and, 4 in β2 (Figure S1A–D, Table 1). With p-values ranging from 4.98 × 10−6 to 2.64 × 10−8, these can be interpreted as strong associations, considering our modest sample size. None of the associated SNPs were detected in more than one phenotype (frequency band) and, to the best of our knowledge, none of these associations were previously reported.
An ANOVA test was computed to evaluate the capacity of the EEG signals for distinguishing between LOAD patients in the different disease stages, subjects with MCI, and controls with normal ageing. RP values distinguished between healthy and affected individuals with statistical confidence, in particular with p-values ranging between 1.42 × 10−7 and 1.02 × 10−3 (Figure 2). This pattern has been observed before in other AD cohorts: an increase in the relative power of slow oscillations (δ and θ rhythms) and a decrease in relative power of fast oscillations (α and β rhythms) [36].
One variant within a candidate for a role in brain wave biology, CLEC16A (C-type lectin domain family 16), passed the genome-wide significance threshold after Bonferroni correction. Individuals with at least one copy of the minor allele of this variant showed higher θ RP values, a tendency also observed in the LOAD patients when compared with controls (Figure 2 and Figure 3). CLEC16A is highly expressed in cerebellum (Figure S2) and in cerebellar Purkinje cells in mouse (Allen Mouse Brain Atlas, http://mouse.brain-map.org/gene/show/50215).

3.2. SNP-Based Functional Analysis of Most Significant GWAS Hits

For the SNP-based analysis, 16,575 variants from all the five linear regressions were jointly selected (significance level α = 0.01) for functional annotation, gene identification and prioritization by positional, eQTL and CI mapping using FUMA [21]. From this analysis, 15 genomic associated loci corresponding to the lead SNPs that passed the significance thresholds and mapping conditions settings were evaluated, with a total of 18 individually significant SNPs, 3 of them in LD. A total of 30 genes were mapped to these regions (Table 2).
This approach retrieved biologically meaningful results, since all the genomic loci mapped to the most significant SNPs by FUMA, either as eQTL or CI analysis, are brain-related. Brain eQTLs were identified in caudate basal ganglia, spinal cord, anterior cingulate cortex and hypothalamus, while CIs were detected mainly in adult and fetal cortex, neural progenitor cells and hippocampus. Moreover, the genes mapped to the most significant SNPs are highly expressed in various brain tissues (Figure S3).
A few genes particularly relevant for brain function were identified with this analysis. Those more directly linked to the phenotypes under study are highlighted next. FOXP2 (Forkhead box protein P2) was mapped both by eQTL and CI to rs12705973 (p-value = 2.65 × 10−6, β1 regression coefficient = 0.02) (Figure S4). This gene encodes a transcriptional repressor that plays a role in synapse formation by regulating SRPX2 levels and was shown to be involved in neural mechanisms mediating the speech development [37,38]. It is highly expressed in the head and tail of nucleus caudatus and putamen, which play a role in movement regulation and other nonmotor actions such as procedural learning, associative learning and inhibitory control of action [39].
Even though the variant rs12263011 physically maps to an intron of FRMD4A (p-value = 4.98 × 10−6, θ regression coefficient = 0.07), by CI it was mapped to CDNF (Cerebral dopamine neurotrophic factor) and SUV39H2 (Histone-lysine N-methyltransferase), in the adult and fetal cortex, respectively (Figure S5). The variant rs12263011 is embedded in a region with epigenetic promoter marks and thus it may modulate the expression of these genes approximately 0.5 Mb apart, through CI. FRMD4A (FERM domain-containing protein 4A) is a scaffolding protein that regulates epithelial polarity and has been previously identified as a genetic risk factor for LOAD and cognitive decline [40,41]. Even though this gene is an interesting candidate for an EEG endophenotype, the mechanism through which rs12263011 or any of the linked intronic SNPs (Figure S6) may contribute to the phenotype is not clear. There is some evidence supporting this variant may alter a transcription factor binding motif (https://pubs.broadinstitute.org/mammals/haploreg/) and thus have a regulatory effect if it lies within a promoter region with enhancer-like features, regulating distal target genes through chromatin loops [42]. CDNF is a trophic factor for dopamine neurons, preventing the 6-hydroxydopamine (6-OHDA)-induced degeneration of dopaminergic neurons in substantia nigra [43]. Neurotrophic factors influence the survival, differentiation and maintenance of neurons in the developing and adult nervous system. SUV39H2 specifically trimethylates Lys-9 of histone H3, which is a tag for epigenetic transcriptional repression. It is involved in the circadian rhythm by being recruited to the E-box elements of the circadian target genes such as PER2 or PER1 [44,45].
FUMA functionally mapped rs7149295 (an intronic variant within NIN gene that encodes the centrosomal protein ninein; p-value = 4.88 × 10−6, δ regression coefficient = 0.09), to L2HGDH (L-2-hydroxyglutarate dehydrogenase) by CI in fetal cortex (Figure S7). L2hgdh KO mice exhibit white matter abnormalities, extensive gliosis, microglia-mediated neuroinflammation and an expansion of oligodendrocyte progenitor cells (OPCs).
The intergenic variant rs55908084 (p-value = 2.65 × 10−6, δ regression coefficient = 0.06) was mapped directly to CACNG4 (voltage-dependent calcium channel gamma-4 subunit), the gene upstream to it. CACNG4 regulates the activity of L-type calcium channels and the trafficking and gating properties of AMPA-selective glutamate receptors (AMPARs), promoting their targeting to the cell membrane and synapses and modulating their gating properties by slowing their rates of activation, deactivation and desensitization [46].
Finally, rs1893824 (p-value = 1.66 × 10−6, θ regression coefficient = −0.05 “better”) was identified as an eQTL of the gene in its vicinity, GALR1 (Galanin receptor type 1), in the anterior cingulate cortex (Figure S8). GALR1 is a receptor for the hormone galanin with the highest expression level in the adenohypophysis [47], and there is evidence supporting that this receptor modulates impulse control in prefrontal-hippocampal circuitry [48].

3.3. Gene-Based Expression and GO Analysis

Gene expression analysis for the 15 lead SNPs showed a significant enrichment of genes up-regulated in the hippocampus, anterior cingulate cortex, cortex and caudate basal ganglia (Figure S3). The same analysis was carried out at the gene level, for the 50 most significant genes (Table 3). Each “genescore” computed by FUMA presents a contribution of all variants mapped to it and there is a correction for the gene size. This revealed an enrichment of up-regulated genes in all brain tissues plus coronary tissue (Figure 4).
Gene-ontology (GO) analysis on this set of 50 genes retrieved “Neuron development” and “neuron differentiation” as the topmost enriched categories, and included the following genes: CAMK1D, PRKG1, CDH23, TENM4, NTM, OPCML, NRXN3, RUNX1, DSCAM, CNTN4, UNCSC, TENM3, MAGI2, CNTNAP2, PTPRD and KDM4C (neuron differentiation) (Figure 5). All the other significantly enriched categories are related to neuron or head development. By the same token, GO term enrichment for cellular component analysis revealed significantly enriched categories related to neuron parts (Figure S9). As expected, in this set of genes, there is an overrepresentation of associations with brain phenotypes and neurological pathologies such as: chronotype, schizophrenia, Asperger syndrome, brain connectivity, short-term memory, dimensional psychopathology, bipolar disorder and AD (Table 4).

3.4. Protein–Protein Interactions

Finally, we performed a protein interaction analysis with STRING for the top 50 protein-coding genes that resulted in a network with significantly more interactions than expected and an enrichment p-value = 2.35 × 10−6 (Figure 6). Some of the interactions in the network, deserve a closer look. CNTN4 (Contactin-4), mediates cell surface interactions during nervous system development and in conjunction with one of its binding partners, amyloid precursor protein (APP), has been shown to promote target-specific axon arborization, highlighting its importance for the functional development of a behaviorally-relevant parallel visual pathway [49]. CNTN5 (Contactin-5) is a similar molecule but it is exclusively expressed in the central nervous system (CNS), with strong expression in the cortex and hippocampus, and has been associated with ASD [50]. NTM (neurotrimin) is a neural cell adhesion molecule of the immunoglobulin superfamily that appears to regulate the development of neuronal projections and might have a role in mediating estrogen-induced peripheral sympathetic innervation; Ntm-deficient mice have shown a deficit in emotional learning [51,52]. RBFOX1 (RNA binding protein fox-1) regulates alternative splicing in tissue-specific exons, and its cytoplasmic target mRNAs are enriched in genes involved in cortical development and autism [53].

3.5. Links to Brain Phenotypes and Neuropathologies

Overall the 16,575 input SNPs were mapped by FUMA to a total of 4078 protein coding genes. In order to identify a possible link between these hits and brain pathophysiology, we inspected the scores for genes previously associated with AD by meta-analysis (Table S1) [54,55]. CNTNAP2 (Contactin-associated protein-like 2), with 12 mapped SNPs, reveals a strong association (p-value = 8.92 × 10−10). We also inspected genes previously associated with schizophrenia by genome-wide association study and replication (Table S2) [56]. The strongest association was found for MAGI2 (p-value = 2.08 × 10−8), with 10 mapped SNPs.
The set of the 50 most significant genes in our analysis is also enriched for genes associated with different brain physiological traits, behavioral phenotypes and neurological disorders (Table 3), highlighting the overlap between the biological pathways of brain oscillations and these traits and diseases. The most significant gene in the genescore analysis, CSMD1 (CUB and sushi domain-containing protein 1), is highly expressed in the CNS, particularly in the frontal cortex (Figure S10) and has been associated with schizophrenia. From the 50 most significant genes, two (ZBTB7C, p-value = 2.29 × 10−12; and FRMD4A, p-value = 9.68 × 10−11) were already identified by a single highly significant SNP for each gene. They have also been linked to different brain pathologies. From the 50 most significant genes, ZBTB7C (p-value = 2.29 × 10−12) and FRMD4A (p-value = 9.68 × 10−11) have been previously linked to different brain pathologies [54]. FRMD4A has been identified as a genetic risk factor for LOAD and may modulate the disease progression by altering tau [57].

3.6. Statistical Power and Effect Size

We evaluated the theoretical statistical power to detect a causal SNP, such as the one identified at genome-wide significance (rs71381191), by taking advantage of the relationship between experimental sample size, allele frequency and effect size. The effect sizes for the discovered variants were reported as b1Table 1, referring to the coefficient of the linear regression. Generally, in QTL analyses, where this coefficient is derived from a regression on a continuous variable with widely distributed values, b1 is not directly interpretable as an effect size per se. However, in our study we considered brain wave RPs, which are normalized measurements, with continuous distribution between 0 and 1. It follows that, in this case, the b1 coefficient derived from the linear regression is interpretable as the relative effect size. As such, the value of b1 of 0.11, obtained for the SNP showing statistically significant association, represents an 11% increase in the RP for the θ brain wave, for each additional “effect” allele.
Using the expressions exposed in Section 2.5, with MAF and β values for the genotyped genome-wide significant variant (rs71381191), we calculated the proportion of phenotypic variance explained by a such causal variant in our population (q2 = 0.0015) and divided it by the total θ RP variance to get the proportion of phenotypic variance explained by the genotyped SNP (R2 = 0.14). Finally, the statistical power to detect a SNP explaining 14% of trait variance (such as this one) with a sample size of 191 individuals and assuming our genome-wide significance threshold, was estimated in 62%. We conclude that the relatively large proportion of RP variance explained by rs71381191 (or a linked variant) has likely contributed to the identification of this significant association in our dataset, for which MAF = 7%.

4. Discussion

The easiness and affordability of EEG has made it very appealing as a potential diagnostic tool for several neurological/neurodegenerative diseases. In addition, the characterization of the biological pathways underlying its measurements are thus of high importance. Malone and colleagues, in one of the largest studies involving EEG signals so far, did not find genome-wide significant hits after Bonferroni correction [58]. This is not surprising since the effects of single genetic variants on multifactorial phenotypes such as brain electrophysiological signals are expected to be very small. The largest GWAS of oscillatory power during eyes-closed resting EEG to date found through a gene-based approach that GABRA2, a known genetic marker for alcohol abuse disorder and epilepsy, was significantly related to β wave power [22]. Out of twenty-four other genes, three were significantly associated with α power, showing differential expression in two tissues: GLYCTK in the hippocampus and GNL3 and ITIH4 in the frontal cortex. All of these three genes were previously associated with schizophrenia and bipolar disorder.

4.1. Linear Regression for 5 Oscillatory Phenotypes Identified One New θ QTL

We have performed a GWAS and QTL analysis for one brain endophenotype, EEG RP in five frequency bands (δ, θ, α, β1 and β2) in an Iberian cohort of LOAD patients, individuals with MCI and controls, corrected for the disease status, as well as for residual population stratification, age and sex. In spite of the modest sample size, one variant within CLEC16A, a C-type lectin, passed the genome-wide significance threshold after Bonferroni correction. CLEC16A was shown to participate in the molecular machinery of human leukocyte antigen (HLA) late endosome formation and trafficking, serving as a direct regulator of the HLA-II pathway in antigen-presenting cells [59]. HLA-II is expressed in microglial cells, which are a component of the innate immunity. The relevance of the adaptive immunity in neurodegeneration has been increasingly recognized. An increase in the number of microglial cells, their activation and the disruption of their functions have been demonstrated in neurodegenerative pathologies such as amyotrophic lateral sclerosis (ALS) and frontotemporal dementia [60]. This finding reinforces the importance of immune regulation in brain physiology. Moreover, evidence points to a function of CLEC16A in Purkinje cells. This neuronal cell type releases the γ-aminobutyric acid (GABA) neurotransmitter that regulates synaptic plasticity and network oscillations through synaptic inhibition by interneurons that release GABA [61,62]. At the physiological level, cerebellar-evoked prefrontal synchronization in the θ frequency range has been shown to be modulated by GABA, being positively associated with working memory performance [63]. Evidence suggests that the cerebellum likely exerts its control on the cortex by a GABAergic dependent set of interneurons and cerebellar θ-burst stimulation modulates cortical excitability of distant interconnected cortical areas [64]. Recently, the hippocampal expression of a GABA receptor has been associated with β oscillations, supporting the now reported association as being relevant to the electrophysiology of the brain [22]. Considering our present results, another link between GABAergic system and brain electrophysiology surfaces. In addition, the inhibition of CLEC16A protein function also has been shown to lead to motor impairments, Purkinje cell loss and impaired autophagy, pointing to a role in the function and clearance of autolysosomes that culminates in neurodegeneration [65].

4.2. SNP and Gene-Based Functional Analysis of Most Significant GWAS Hits and Their Link to Neurological Traits and Pathologies

The analysis with the most nominally significant SNPs from all five linear regressions retrieved a total of 30 candidates, all with functional relevance in brain tissues, as suggested by expression and epigenetic data. Most of the loci physically or functionally mapped by FUMA to the lead SNPs (SNPs that passed the threshold criteria) are involved in brain function, namely neuronal development and maintenance (CDNF), cognitive function (FOXP2), circadian rhythm (SUV39H2) and one is likely a key player in brain electrophysiology (GALR1) [39,43,45,66]. Due to their relevance to brain circuitry and physiology these genes have also proven or potential roles in several diseases characterized by impaired cognitive function. Indeed, FRMD4A (FERM domain-containing protein 4A) is a scaffolding protein that regulates epithelial polarity and has been previously identified as a genetic risk factor for LOAD [57,67]. This connection to AD is also reinforced by the interaction with PAR3 which regulates CNTNAP2 spatial localization. CNTNAP2 is required for radial and longitudinal organization of myelinated axons and has been associated with AD and other disorders such as epilepsy, seizures, autism and schizophrenia, and common variants in this gene influence early language acquisition [54,55,68,69,70,71]. Neuropathological analyses in mice with mutations of this gene revealed abnormalities in neuron migration, reduced number of interneurons and abnormal neuronal network activity [72]. Moreover, a homozygous mutation in FRMD4A has been linked to a syndrome of congenital microcephaly and intellectual disability [73]. On the other hand, variants in the ZBTB7C (Zinc finger and BTB domain-containing protein 7C) have been associated with ischemic injury susceptibility, maybe by modulating the ischemic response via neuronal apoptosis [74]. Another gene among the top hits, L2HGDH, is involved in L-2-hydroxyglutaric aciduria (L2HGA), a rare autosomal recessive disorder clinically characterized by a mild psychomotor delay followed by progressive cerebellar ataxia and moderate to severe intellectual disability and a tendency to the development of malignant brain tumors [75]. Moreover, L2hgdh deficiency leads to impaired adult hippocampal neurogenesis and late-onset neurodegeneration in mouse brains [76]. It is thus plausible that variants in L2HGDH may increase the risk for neurodegenerative disorders. Finally, GALR1 has been linked to temporal lobe epilepsy (TLE) and galanin agonists inhibit seizures [47].
More indirect evidence for a contribution to neurological diseases exists for other genes in this set. CDNF, a trophic factor for dopamine neurons that prevents the 6-hydroxydopamine (6-OHDA)-induced degeneration of dopaminergic neurons in substantia nigra may have a role in the evolution of Parkinson’s disease [43]. SUV39H, as an essential part of the circadian system, a timing mechanism responsible for orchestrating many physiological processes including behavior and cognition through epigenetic mechanisms, is a candidate gene for autism susceptibility [77].
As previously mentioned, gene-based approaches increase the power to detect genes affecting the phenotype under study [19,21]. The set of the 50 most significant genes in our analysis is also enriched in genes associated with different brain physiological traits, behavioral phenotypes and neurological disorders, highlighting the overlap between the biological pathways of brain oscillations and these traits and diseases. Indeed, EEG has been a valuable tool in the study of chronotype, brain connectivity, substance abuse, Asperger syndrome, the diagnosis and classification of schizophrenia, and even loneliness [78,79,80,81,82,83].
The most significant gene in the genescore analysis was CSMD1 (CUB and sushi domain-containing protein 1), which is highly expressed in frontal cortex. Neurophysiological deficits have been observed in CSMD1 depleted mice, inducing blunted emotional responses, anxiety and depression [84,85], and it has also been linked to schizophrenia. We also inspected genes previously associated with schizophrenia by genome-wide association study and replication [56]. The strongest association was found for MAGI2 (p-value = 2.08 × 10−8), with 10 mapped SNPs. This gene encodes a molecule that serves as a scaffold for proteins assembling synaptic protein complexes, therefore, with an essential role in synaptic development and maintenance [86]. In addition, common variants of this gene were associated with cognitive impairment in individuals with schizophrenia, and it has been validated as a strong candidate by genome-wide association [56,87].

4.3. Limitations and Future Research Lines

Our study has some limitations that should be addressed. On one hand, the EEG analyses have been performed using the grand-averaged values in order to reduce the dimensionality of the results and simplify statistical analyses. As EEG activity may differ depending on the brain region under study, it could be useful to identify specific affected scalp regions in future studies. For future works, it will be also interesting to acquire and analyze the sleep patterns of the subjects during the night previous to the EEG recording, since some studies reported that disturbances in the sleep-wake cycle and circadian rhythms are common symptoms of AD [88]. On the other hand, it is noteworthy the limited number of subjects analyzed, as already acknowledged, as well as the inclusion of only two closely related populations. The multimodal analysis of larger cohorts, from several diverse populations, will be important to replicate our signals and may uncover other associations between genomic and EEG data.

5. Conclusions

One novel association was found with an interesting candidate for a role in brain wave biology, CLEC16A (C-type lectin domain family 16), with a variant at this locus passing the adjusted genome-wide significance threshold after Bonferroni correction, reinforcing the importance of the immune regulation in brain function. Moreover, at a significance cutoff value of 5 × 10−6, 18 independent association signals were detected. These signals comprise brain expression Quantitative Loci (eQTLs) in caudate basal ganglia, spinal cord, anterior cingulate cortex and hypothalamus, as well as chromatin interactions in adult and fetal cortex, neural progenitor cells and hippocampus. At the same time, in the set of genes showing signals of association with brain wave RP in our dataset, there is an overrepresentation of loci previously associated with neurological traits and pathologies, evidencing the pleiotropy of the genetic variation modulating brain function.
Our results corroborate and strengthen previous findings regarding the biological pathways involved in brain electrophysiology, namely, the role of immunity regulation and GABA neurotransmission, through the identification of a novel candidate gene for brain wave RP modulation. In addition, the complexity inherent to brain phenotypes and the pleiotropy of the variants with potential to modulate brain function is evidenced. In fact, even when dissecting isolated oscillatory endophenotypes, the associated genetic variation has the potential to affect the regulation of other traits and disease risk through diverse mechanisms, such as direct regulation of other genes or by modulating the interaction with other proteins.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-3425/10/11/870/s1, Figure S1: Linear regression p-values for each frequency band; Figure S2: Gene expression for CLEC16A in 54 tissues from GTEx v8; Figure S3: Significantly enriched up-regulated DEGs (differentially expressed genes) for the 15 lead SNPs on GTEx v8 54 tissue types; Figure S4: Circos plot for the rs12705973 genomic associated locus (FOXP2); Figure S5: Circos plot for the rs12705973 genomic associated locus (FRMD4A); Figure S6: Linkage proxies for the variant rs12263011, within FRMD4A; Figure S7: Circos plot for the rs7149295 genomic associated locus (L2HGDH); Figure S8: Circos plot for the rs1893824 genomic associated locus (GALR1); Figure S9: GO (gene-ontology) most significant categories for cellular component; Figure S10: Gene expression for CSMD1 in 54 tissues from GTEx v8; Table S1: Genescores of genes previously associated with AD; Table S2: Genescores of genes previously associated with schizophrenia.

Author Contributions

Conceptualization of this work was developed by A.M.L., N.P., C.G., R.H., M.A. and L.Á., and methodology was developed by A.M.L. and M.Â.R. Project administration was led by N.P., C.G., R.H., L.D. and M.F. Subjects were recruited by L.D., P.S., M.F., M.R., and C.P., which also obtained and handled demographic data, collected saliva samples and supported EEG recordings. C.G., J.P., A.M.-C. and S.J.R.-G. recorded and handled E.E.G. data. I.G. and S.M. handled and pre-processed genetic samples. Statistical analyses were developed by M.Â.R., supervised by A.M.L. and N.P. Original draft preparation was led by M.Â.R. and A.M.L. All the authors have contributed to the writing of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This project is 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 1317_AD-EEGWA), (Cooperation 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, NP and IG are funded by FCT: CEECIND/00684/2017, IF/01262/2014, SFRH/BPD/97414/2013 and CEECIND/02609/2017, respectively. MA is funded by the Grant RYC-2015-18241 from the Spanish Government. 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.

Acknowledgments

To 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; Armonía Centro de Día, Zamora.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Da Silva, F.L. EEG: Origin and Measurement. In EEG-fMRI; Mulert, C., Lemieux, L., Eds.; Springer: Berlin/Heidelberg, Germany, 2009; pp. 19–38. ISBN 978-3-540-87918-3. [Google Scholar]
  2. Nunez, P.L.; Srinivasan, R.; Fields, R.D. EEG functional connectivity, axon delays and white matter disease. Clin. Neurophysiol. 2015, 126, 110–120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Voytek, B.; Knight, R.T. Dynamic Network Communication as a Unifying Neural Basis for Cognition, Development, Aging, and Disease. Biol. Psychiatry 2015, 77, 1089–1097. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Babiloni, C.; Blinowska, K.; Bonanni, L.; Cichocki, A.; De Haan, W.; Del Percio, C.; Dubois, B.; Escudero, J.; Fernández, A.; Frisoni, G.; et al. What electrophysiology tells us about Alzheimer’s disease: A window into the synchronization and connectivity of brain neurons. Neurobiol. Aging 2020, 85, 58–73. [Google Scholar] [CrossRef] [PubMed]
  5. Cremer, R.; Zeef, E.J. What Kind of Noise Increases With Age? J. Gerontol. 1987, 42, 515–518. [Google Scholar] [CrossRef]
  6. Holtzheimer, P.E.; Mayberg, H.S. Stuck in a rut: Rethinking depression and its treatment. Trends Neurosci. 2011, 34, 1–9. [Google Scholar] [CrossRef] [Green Version]
  7. Jelic, V.; Johansson, S.-E.; Almkvist, O.; Shigeta, M.; Julin, P.; Nordberg, A.; Winblad, B.; Wahlund, L.-O. Quantitative electroencephalography in mild cognitive impairment: Longitudinal changes and possible prediction of Alzheimer’s disease. Neurobiol. Aging 2000, 21, 533–540. [Google Scholar] [CrossRef]
  8. Jeong, J. EEG dynamics in patients with Alzheimer’s disease. Clin. Neurophysiol. 2004, 115, 1490–1505. [Google Scholar] [CrossRef]
  9. Weinger, P.M.; Zemon, V.; Soorya, L.; Gordon, J. Low-contrast response deficits and increased neural noise in children with autism spectrum disorder. Neuropsychologia 2014, 63, 10–18. [Google Scholar] [CrossRef]
  10. Winterer, G.; Weinberger, D.R. Genes, dopamine and cortical signal-to-noise ratio in schizophrenia. Trends Neurosci. 2004, 27, 683–690. [Google Scholar] [CrossRef]
  11. Tang, Y.; Chorlian, D.B.; Rangaswamy, M.; Porjesz, B.; Bauer, L.; Kuperman, S.; O’Connor, S.; Rohrbaugh, J.; Schuckit, M.; Stimus, A.; et al. Genetic influences on bipolar EEG power spectra. Int. J. Psychophysiol. 2007, 65, 2–9. [Google Scholar] [CrossRef]
  12. Zietsch, B.P.; Hansen, J.L.; Hansell, N.K.; Geffen, G.M.; Martin, N.G.; Wright, M.J. Common and specific genetic influences on EEG power bands delta, theta, alpha, and beta. Biol. Psychol. 2007, 75, 154–164. [Google Scholar] [CrossRef] [PubMed]
  13. Smit, D.J.A.; Posthuma, D.; Boomsma, D.I.; Geus, E.J.C. Heritability of background EEG across the power spectrum. Psychophysiology 2005, 42, 691–697. [Google Scholar] [CrossRef] [PubMed]
  14. De Geus, E.J. From genotype to EEG endophenotype: A route for post-genomic understanding of complex psychiatric disease? Genome Med. 2010, 2, 63. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Loo, S.K.; Lenartowicz, A.; Makeig, S. Research Review: Use of EEG biomarkers in child psychiatry research-current state and future directions. J. Child Psychol. Psychiatry 2016, 57, 4–17. [Google Scholar] [CrossRef]
  16. De Geus, E.J.C.; Boomsma, D.I. A Genetic Neuroscience Approach to Human Cognition. Eur. Psychol. 2001, 6, 241–253. [Google Scholar] [CrossRef]
  17. Gottesman, I.I.; Gould, T.D. The Endophenotype Concept in Psychiatry: Etymology and Strategic Intentions. Am. J. Psychiatry 2003, 160, 636–645. [Google Scholar] [CrossRef]
  18. Psychiatric GWAS Consortium Coordinating Committee. Genomewide Association Studies: History, Rationale, and Prospects for Psychiatric Disorders. Am. J. Psychiatry 2009, 166, 540–556. [Google Scholar] [CrossRef] [Green Version]
  19. Lamparter, D.; Marbach, D.; Rueedi, R.; Kutalik, Z.; Bergmann, S. Fast and Rigorous Computation of Gene and Pathway Scores from SNP-Based Summary Statistics. PLoS Comput. Biol. 2016, 12, e1004714. [Google Scholar] [CrossRef]
  20. De Leeuw, C.A.; Mooij, J.M.; Heskes, T.; Posthuma, D. MAGMA: Generalized Gene-Set Analysis of GWAS Data. PLoS Comput. Biol. 2015, 11, e1004219. [Google Scholar] [CrossRef]
  21. Watanabe, K.; Taskesen, E.; van Bochoven, A.; Posthuma, D. Functional mapping and annotation of genetic associations with FUMA. Nat. Commun. 2017, 8, 1826. [Google Scholar] [CrossRef] [Green Version]
  22. Smit, D.J.A.; Wright, M.J.; Meyers, J.L.; Martin, N.G.; Ho, Y.Y.W.; Malone, S.M.; Zhang, J.; Burwell, S.J.; Chorlian, D.B.; de Geus, E.J.C.; et al. Genome-wide association analysis links multiple psychiatric liability genes to oscillatory brain activity. Hum. Brain Mapp. 2018, 39, 4183–4195. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Jack, C.R.; Bennett, D.A.; Blennow, K.; Carrillo, M.C.; Dunn, B.; Haeberlein, S.B.; Holtzman, D.M.; Jagust, W.; Jessen, F.; Karlawish, J.; et al. NIA-AA Research Framework: Toward a biological definition of Alzheimer’s disease. Alzheimers Dement. 2018, 14, 535–562. [Google Scholar] [CrossRef] [PubMed]
  24. Núñez, P.; Poza, J.; Gómez, C.; Rodríguez-González, V.; Hillebrand, A.; Tola-Arribas, M.A.; Cano, M.; Hornero, R. 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. 2019, 16, 056030. [Google Scholar] [CrossRef]
  25. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.W.; Daly, M.J.; et al. PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Price, A.L.; Patterson, N.J.; Plenge, R.M.; Weinblatt, M.E.; Shadick, N.A.; Reich, D. Principal components analysis corrects for stratification in genome-wide association studies. Nat. Genet. 2006, 38, 904–909. [Google Scholar] [CrossRef] [PubMed]
  27. Gómez, C.; Pérez-Macías, J.M.; Poza, J.; Fernández, A.; Hornero, R. Spectral changes in spontaneous MEG activity across the lifespan. J. Neural Eng. 2013, 10, 066006. [Google Scholar] [CrossRef]
  28. Patterson, N.; Price, A.L.; Reich, D. Population Structure and Eigenanalysis. PLoS Genet. 2006, 2, e190. [Google Scholar] [CrossRef]
  29. Schmitt, A.D.; Hu, M.; Jung, I.; Xu, Z.; Qiu, Y.; Tan, C.L.; Li, Y.; Lin, S.; Lin, Y.; Barr, C.L.; et al. A Compendium of Chromatin Contact Maps Reveals Spatially Active Regions in the Human Genome. Cell Rep. 2016, 17, 2042–2059. [Google Scholar] [CrossRef] [Green Version]
  30. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING v11: Protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019, 47, D607–D613. [Google Scholar] [CrossRef] [Green Version]
  31. Welter, D.; MacArthur, J.; Morales, J.; Burdett, T.; Hall, P.; Junkins, H.; Klemm, A.; Flicek, P.; Manolio, T.; Hindorff, L.; et al. The NHGRI GWAS Catalog, a curated resource of SNP-trait associations. Nucleic Acids Res. 2014, 42, D1001–D1006. [Google Scholar] [CrossRef]
  32. Visscher, P.M.; Wray, N.R.; Zhang, Q.; Sklar, P.; McCarthy, M.I.; Brown, M.A.; Yang, J. 10 Years of GWAS Discovery: Biology, Function, and Translation. Am. J. Hum. Genet. 2017, 101, 5–22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Weinberg, W. Über Vererbungsgesetze beim Menschen. Z. Indukt. Abstamm.-Vererb. 1908, 1, 440–460. [Google Scholar] [CrossRef] [Green Version]
  34. Hardy, G.H. Mendelian proportions in a mixed population. Science 1908, 28, 49–50. [Google Scholar] [CrossRef]
  35. Shim, H.; Chasman, D.I.; Smith, J.D.; Mora, S.; Ridker, P.M.; Nickerson, D.A.; Krauss, R.M.; Stephens, M. A Multivariate Genome-Wide Association Analysis of 10 LDL Subfractions, and Their Response to Statin Treatment, in 1868 Caucasians. PLoS ONE 2015, 10, e0120758. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Hamm, V.; Héraud, C.; Cassel, J.-C.; Mathis, C.; Goutagny, R. Precocious Alterations of Brain Oscillatory Activity in Alzheimer’s Disease: A Window of Opportunity for Early Diagnosis and Treatment. Front. Cell. Neurosci. 2015, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Fisher, S.E.; Scharff, C. FOXP2 as a molecular window into speech and language. Trends Genet. 2009, 25, 166–177. [Google Scholar] [CrossRef] [PubMed]
  38. Konopka, G.; Bomar, J.M.; Winden, K.; Coppola, G.; Jonsson, Z.O.; Gao, F.; Peng, S.; Preuss, T.M.; Wohlschlegel, J.A.; Geschwind, D.H. Human-specific transcriptional regulation of CNS development genes by FOXP2. Nature 2009, 462, 213–217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Teramitsu, I. Parallel FoxP1 and FoxP2 Expression in Songbird and Human Brain Predicts Functional Interaction. J. Neurosci. 2004, 24, 3152–3163. [Google Scholar] [CrossRef] [PubMed]
  40. Andrews, S.J.; Das, D.; Cherbuin, N.; Anstey, K.J.; Easteal, S. Association of genetic risk factors with cognitive decline: The PATH through life project. Neurobiol. Aging 2016, 41, 150–158. [Google Scholar] [CrossRef] [PubMed]
  41. Huovinen, J.; Helisalmi, S.; Paananen, J.; Laiterä, T.; Kojoukhova, M.; Sutela, A.; Vanninen, R.; Laitinen, M.; Rauramaa, T.; Koivisto, A.M.; et al. Alzheimer’s Disease-Related Polymorphisms in Shunt-Responsive Idiopathic Normal Pressure Hydrocephalus. J. Alzheimers Dis. 2017, 60, 1077–1085. [Google Scholar] [CrossRef] [Green Version]
  42. Jung, I.; Schmitt, A.; Diao, Y.; Lee, A.J.; Liu, T.; Yang, D.; Tan, C.; Eom, J.; Chan, M.; Chee, S.; et al. A compendium of promoter-centered long-range chromatin interactions in the human genome. Nat. Genet. 2019, 51, 1442–1449. [Google Scholar] [CrossRef] [PubMed]
  43. Lindholm, P.; Voutilainen, M.H.; Laurén, J.; Peränen, J.; Leppänen, V.-M.; Andressoo, J.-O.; Lindahl, M.; Janhunen, S.; Kalkkinen, N.; Timmusk, T.; et al. Novel neurotrophic factor CDNF protects and rescues midbrain dopamine neurons in vivo. Nature 2007, 448, 73–77. [Google Scholar] [CrossRef] [PubMed]
  44. Ait-Si-Ali, S.; Guasconi, V.; Fritsch, L.; Yahi, H.; Sekhri, R.; Naguibneva, I.; Robin, P.; Cabon, F.; Polesskaya, A.; Harel-Bellan, A. A Suv39h-dependent mechanism for silencing S-phase genes in differentiating but not in cycling cells. EMBO J. 2004, 23, 605–615. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Duong, H.A.; Weitz, C.J. Temporal orchestration of repressive chromatin modifiers by circadian clock Period complexes. Nat. Struct. Mol. Biol. 2014, 21, 126–132. [Google Scholar] [CrossRef]
  46. Kato, A.S.; Gill, M.B.; Ho, M.T.; Yu, H.; Tu, Y.; Siuda, E.R.; Wang, H.; Qian, Y.-W.; Nisenbaum, E.S.; Tomita, S.; et al. Hippocampal AMPA Receptor Gating Controlled by Both TARP and Cornichon Proteins. Neuron 2010, 68, 1082–1096. [Google Scholar] [CrossRef] [Green Version]
  47. Guipponi, M.; Chentouf, A.; Webling, K.E.B.; Freimann, K.; Crespel, A.; Nobile, C.; Lemke, J.R.; Hansen, J.; Dorn, T.; Lesca, G.; et al. Galanin pathogenic mutations in temporal lobe epilepsy. Hum. Mol. Genet. 2015, 24, 3082–3091. [Google Scholar] [CrossRef] [Green Version]
  48. Messanvi, F.; Perkins, A.; du Hoffmann, J.; Chudasama, Y. Fronto-temporal galanin modulates impulse control. Psychopharmacology 2019. [Google Scholar] [CrossRef]
  49. Osterhout, J.A.; Stafford, B.K.; Nguyen, P.L.; Yoshihara, Y.; Huberman, A.D. Contactin-4 Mediates Axon-Target Specificity and Functional Development of the Accessory Optic System. Neuron 2015, 86, 985–999. [Google Scholar] [CrossRef] [Green Version]
  50. Kleijer, K.T.E.; van Nieuwenhuize, D.; Spierenburg, H.A.; Gregorio-Jordan, S.; Kas, M.J.H.; Burbach, J.P.H. Structural abnormalities in the primary somatosensory cortex and a normal behavioral profile in Contactin-5 deficient mice. Cell Adhes. Migr. 2018, 12, 5–18. [Google Scholar] [CrossRef]
  51. Gil, O.D.; Zanazzi, G.; Struyk, A.F.; Salzer, J.L. Neurotrimin Mediates Bifunctional Effects on Neurite Outgrowth via Homophilic and Heterophilic Interactions. J. Neurosci. 1998, 18, 9312–9325. [Google Scholar] [CrossRef] [Green Version]
  52. Krizsan-Agbas, D.; Pedchenko, T.; Smith, P.G. Neurotrimin is an estrogen-regulated determinant of peripheral sympathetic innervation. J. Neurosci. Res. 2008, 86, 3086–3095. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Lee, J.-A.; Damianov, A.; Lin, C.-H.; Fontes, M.; Parikshak, N.N.; Anderson, E.S.; Geschwind, D.H.; Black, D.L.; Martin, K.C. Cytoplasmic Rbfox1 Regulates the Expression of Synaptic and Autism-Related Genes. Neuron 2016, 89, 113–128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. European Alzheimer’s Disease Initiative (EADI); Genetic and Environmental Risk in Alzheimer’s Disease (GERAD); Alzheimer’s Disease Genetic Consortium (ADGC); Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE); Lambert, J.-C.; Ibrahim-Verbaas, C.A.; Harold, D.; Naj, A.C.; Sims, R.; Bellenguez, C.; et al. Meta-analysis of 74,046 individuals identifies 11 new susceptibility loci for Alzheimer’s disease. Nat. Genet. 2013, 45, 1452–1458. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Jansen, I.E.; Savage, J.E.; Watanabe, K.; Bryois, J.; Williams, D.M.; Steinberg, S.; Sealock, J.; Karlsson, I.K.; Hägg, S.; Athanasiu, L.; et al. Genome-wide meta-analysis identifies new loci and functional pathways influencing Alzheimer’s disease risk. Nat. Genet. 2019, 51, 404–413. [Google Scholar] [CrossRef]
  56. Li, Z.; Chen, J.; Yu, H.; He, L.; Xu, Y.; Zhang, D.; Yi, Q.; Li, C.; Li, X.; Shen, J.; et al. Genome-wide association analysis identifies 30 new susceptibility loci for schizophrenia. Nat. Genet. 2017, 49, 1576–1583. [Google Scholar] [CrossRef]
  57. Yan, X.; Nykänen, N.-P.; Brunello, C.A.; Haapasalo, A.; Hiltunen, M.; Uronen, R.-L.; Huttunen, H.J. FRMD4A–cytohesin signaling modulates the cellular release of tau. J. Cell Sci. 2016, 129, 2003–2015. [Google Scholar] [CrossRef] [Green Version]
  58. Malone, S.M.; Burwell, S.J.; Vaidyanathan, U.; Miller, M.B.; Mcgue, M.; Iacono, W.G. Heritability and molecular-genetic basis of resting EEG activity: A genome-wide association study: Genome-wide association study of resting EEG. Psychophysiology 2014, 51, 1225–1245. [Google Scholar] [CrossRef] [Green Version]
  59. Van Luijn, M.M.; Kreft, K.L.; Jongsma, M.L.; Mes, S.W.; Wierenga-Wolf, A.F.; van Meurs, M.; Melief, M.-J.; van der Kant, R.; Janssen, L.; Janssen, H.; et al. Multiple sclerosis-associated CLEC16A controls HLA class II expression via late endosome biogenesis. Brain 2015, 138, 1531–1547. [Google Scholar] [CrossRef] [Green Version]
  60. Radford, R.A.; Morsch, M.; Rayner, S.L.; Cole, N.J.; Pountney, D.L.; Chung, R.S. The established and emerging roles of astrocytes and microglia in amyotrophic lateral sclerosis and frontotemporal dementia. Front. Cell. Neurosci. 2015, 9. [Google Scholar] [CrossRef] [Green Version]
  61. Gong, N.; Li, Y.; Cai, G.-Q.; Niu, R.-F.; Fang, Q.; Wu, K.; Chen, Z.; Lin, L.-N.; Xu, L.; Fei, J.; et al. GABA Transporter-1 Activity Modulates Hippocampal Theta Oscillation and Theta Burst Stimulation-Induced Long-Term Potentiation. J. Neurosci. 2009, 29, 15836–15845. [Google Scholar] [CrossRef] [Green Version]
  62. Van Welie, I.; Smith, I.T.; Watt, A.J. The metamorphosis of the developing cerebellar microcircuit. Curr. Opin. Neurobiol. 2011, 21, 245–253. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Du, X.; Rowland, L.M.; Summerfelt, A.; Choa, F.-S.; Wittenberg, G.F.; Wisner, K.; Wijtenburg, A.; Chiappelli, J.; Kochunov, P.; Hong, L.E. Cerebellar-Stimulation Evoked Prefrontal Electrical Synchrony Is Modulated by GABA. Cerebellum 2018, 17, 550–563. [Google Scholar] [CrossRef] [PubMed]
  64. Casula, E.P.; Pellicciari, M.C.; Ponzo, V.; Stampanoni Bassi, M.; Veniero, D.; Caltagirone, C.; Koch, G. Cerebellar theta burst stimulation modulates the neural activity of interconnected parietal and motor areas. Sci. Rep. 2016, 6, 36191. [Google Scholar] [CrossRef] [Green Version]
  65. Redmann, V.; Lamb, C.A.; Hwang, S.; Orchard, R.C.; Kim, S.; Razi, M.; Milam, A.; Park, S.; Yokoyama, C.C.; Kambal, A.; et al. Clec16a is Critical for Autolysosome Function and Purkinje Cell Survival. Sci. Rep. 2016, 6, 23326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Bruce, H.A.; Margolis, R.L. FOXP2: Novel exons, splice variants, and CAG repeat length stability. Hum. Genet. 2002, 111, 136–144. [Google Scholar] [CrossRef] [PubMed]
  67. Ikenouchi, J.; Umeda, M. FRMD4A regulates epithelial polarity by connecting Arf6 activation with the PAR complex. Proc. Natl. Acad. Sci. USA 2010, 107, 748–753. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Alarcón, M.; Abrahams, B.S.; Stone, J.L.; Duvall, J.A.; Perederiy, J.V.; Bomar, J.M.; Sebat, J.; Wigler, M.; Martin, C.L.; Ledbetter, D.H.; et al. Linkage, Association, and Gene-Expression Analyses Identify CNTNAP2 as an Autism-Susceptibility Gene. Am. J. Hum. Genet. 2008, 82, 150–159. [Google Scholar] [CrossRef] [Green Version]
  69. Friedman, J.I.; Vrijenhoek, T.; Markx, S.; Janssen, I.M.; van der Vliet, W.A.; Faas, B.H.W.; Knoers, N.V.; Cahn, W.; Kahn, R.S.; Edelmann, L.; et al. CNTNAP2 gene dosage variation is associated with schizophrenia and epilepsy. Mol. Psychiatry 2008, 13, 261–266. [Google Scholar] [CrossRef]
  70. Rodenas-Cuadrado, P.M.; Mengede, J.; Baas, L.; Devanna, P.; Schmid, T.A.; Yartsev, M.; Firzlaff, U.; Vernes, S.C. Mapping the distribution of language related genes FoxP1, FoxP2, and CntnaP2 in the brains of vocal learning bat species. J. Comp. Neurol. 2018, 526, 1235–1266. [Google Scholar] [CrossRef] [Green Version]
  71. Whitehouse, A.J.O.; Bishop, D.V.M.; Ang, Q.W.; Pennell, C.E.; Fisher, S.E. CNTNAP2 variants affect early language development in the general population. Genes Brain Behav. 2011, 10, 451–456. [Google Scholar] [CrossRef]
  72. Peñagarikano, O.; Abrahams, B.S.; Herman, E.I.; Winden, K.D.; Gdalyahu, A.; Dong, H.; Sonnenblick, L.I.; Gruver, R.; Almajano, J.; Bragin, A.; et al. Absence of CNTNAP2 Leads to Epilepsy, Neuronal Migration Abnormalities, and Core Autism-Related Deficits. Cell 2011, 147, 235–246. [Google Scholar] [CrossRef] [Green Version]
  73. Fine, D.; Flusser, H.; Markus, B.; Shorer, Z.; Gradstein, L.; Khateeb, S.; Langer, Y.; Narkis, G.; Birk, R.; Galil, A.; et al. A syndrome of congenital microcephaly, intellectual disability and dysmorphism with a homozygous mutation in FRMD4A. Eur. J. Hum. Genet. 2015, 23, 1729–1734. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Du, R.; Zhou, J.; Lorenzano, S.; Liu, W.; Charoenvimolphan, N.; Qian, B.; Xu, J.; Wang, J.; Zhang, X.; Wang, X.; et al. Integrative Mouse and Human Studies Implicate ANGPT1 and ZBTB7C as Susceptibility Genes to Ischemic Injury. Stroke 2015, 46, 3514–3522. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Haliloglu, G.; Jobard, F.; Oguz, K.; Anlar, B.; Akalan, N.; Coskun, T.; Sass, J.; Fischer, J.; Topcu, M. L-2-Hydroxyglutaric Aciduria and Brain Tumors in Children with Mutations in the L2HGDH Gene: Neuroimaging Findings. Neuropediatrics 2008, 39, 119–122. [Google Scholar] [CrossRef] [PubMed]
  76. Ma, S.; Sun, R.; Jiang, B.; Gao, J.; Deng, W.; Liu, P.; He, R.; Cui, J.; Ji, M.; Yi, W.; et al. L2hgdh Deficiency Accumulates l-2-Hydroxyglutarate with Progressive Leukoencephalopathy and Neurodegeneration. Mol. Cell. Biol. 2017, 37, e00492-16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Feng, D.; Lazar, M.A. Clocks, Metabolism, and the Epigenome. Mol. Cell 2012, 47, 158–167. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Cacioppo, S.; Bangee, M.; Balogh, S.; Cardenas-Iniguez, C.; Qualter, P.; Cacioppo, J.T. Loneliness and implicit attention to social threat: A high-performance electrical neuroimaging study. Cogn. Neurosci. 2016, 7, 138–159. [Google Scholar] [CrossRef]
  79. Ceballos, N.A.; Bauer, L.O.; Houston, R.J. Recent EEG and ERP Findings in Substance Abusers. Clin. EEG Neurosci. 2009, 40, 122–128. [Google Scholar] [CrossRef] [Green Version]
  80. Clarke, A.R.; Barry, R.J.; Indraratna, A.; Dupuy, F.E.; McCarthy, R.; Selikowitz, M. EEG activity in children with Asperger’s Syndrome. Clin. Neurophysiol. 2016, 127, 442–451. [Google Scholar] [CrossRef] [Green Version]
  81. Dvey-Aharon, Z.; Fogelson, N.; Peled, A.; Intrator, N. Schizophrenia Detection and Classification by Advanced Analysis of EEG Recordings Using a Single Electrode Approach. PLoS ONE 2015, 10, e0123033. [Google Scholar] [CrossRef] [Green Version]
  82. Jamal, W.; Das, S.; Maharatna, K.; Pan, I.; Kuyucu, D. Brain connectivity analysis from EEG signals using stable phase-synchronized states during face perception tasks. Phys. A Stat. Mech. Appl. 2015, 434, 273–295. [Google Scholar] [CrossRef] [Green Version]
  83. Putilov, A.A.; Donskaya, O.G.; Verevkin, E.G.; Putilov, D.A. Associations of waking EEG structure with chronotype and trototype of 130 sleep deprived individuals. Biol. Rhythm Res. 2010, 41, 113–136. [Google Scholar] [CrossRef]
  84. Håvik, B.; Le Hellard, S.; Rietschel, M.; Lybæk, H.; Djurovic, S.; Mattheisen, M.; Mühleisen, T.W.; Degenhardt, F.; Priebe, L.; Maier, W.; et al. The Complement Control-Related Genes CSMD1 and CSMD2 Associate to Schizophrenia. Biol. Psychiatry 2011, 70, 35–42. [Google Scholar] [CrossRef] [PubMed]
  85. Steen, V.M.; Nepal, C.; Ersland, K.M.; Holdhus, R.; Nævdal, M.; Ratvik, S.M.; Skrede, S.; Håvik, B. Neuropsychological Deficits in Mice Depleted of the Schizophrenia Susceptibility Gene CSMD1. PLoS ONE 2013, 8, e79501. [Google Scholar] [CrossRef] [PubMed]
  86. Bauß, K.; Knapp, B.; Jores, P.; Roepman, R.; Kremer, H.; v. Wijk, E.; Märker, T.; Wolfrum, U. Phosphorylation of the Usher syndrome 1G protein SANS controls Magi2-mediated endocytosis. Hum. Mol. Genet. 2014, 23, 3923–3942. [Google Scholar] [CrossRef] [Green Version]
  87. Koide, T.; Banno, M.; Aleksic, B.; Yamashita, S.; Kikuchi, T.; Kohmura, K.; Adachi, Y.; Kawano, N.; Kushima, I.; Nakamura, Y.; et al. Common Variants in MAGI2 Gene Are Associated with Increased Risk for Cognitive Impairment in Schizophrenic Patients. PLoS ONE 2012, 7, e36836. [Google Scholar] [CrossRef]
  88. Musiek, E.S. Circadian clock disruption in neurodegenerative diseases: Cause and effect? Front. Pharmacol. 2015, 6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Quantile–quantile plot for the five linear regressions. Ordered quantiles of the log-transformed p-values for each of the five linear regressions (δ: blue; θ: yellow; α: green; β1: red; β2: purple) by the theoretical quantiles. A uniform distribution of p-values should follow the 45° line in orange.
Figure 1. Quantile–quantile plot for the five linear regressions. Ordered quantiles of the log-transformed p-values for each of the five linear regressions (δ: blue; θ: yellow; α: green; β1: red; β2: purple) by the theoretical quantiles. A uniform distribution of p-values should follow the 45° line in orange.
Brainsci 10 00870 g001
Figure 2. Box plots for the relative power distributions in each frequency band for each disease status. Individuals with outlier relative power values are marked as diamonds. Each color represents a disease status (CTRL: blue; MCI: orange; MIL: green; MOD: red; SEV: purple).
Figure 2. Box plots for the relative power distributions in each frequency band for each disease status. Individuals with outlier relative power values are marked as diamonds. Each color represents a disease status (CTRL: blue; MCI: orange; MIL: green; MOD: red; SEV: purple).
Brainsci 10 00870 g002
Figure 3. Relative power for each frequency band and genotype for rs71381191, rs12705973, rs12263011, rs7149295, rs55908084 and rs1893824. The number of altered alleles is color coded: blue for individuals with 0 altered alleles (homozygous for the most frequent allele), orange for heterozygous individuals and green for individuals homozygous for the altered (minor) allele.
Figure 3. Relative power for each frequency band and genotype for rs71381191, rs12705973, rs12263011, rs7149295, rs55908084 and rs1893824. The number of altered alleles is color coded: blue for individuals with 0 altered alleles (homozygous for the most frequent allele), orange for heterozygous individuals and green for individuals homozygous for the altered (minor) allele.
Brainsci 10 00870 g003
Figure 4. Significantly enriched DEG (differentially expressed genes) for the 50 top genes on GTEx v8 54 tissue types. Log-transformed p-values for the enrichment of DEG in each tissue. Significantly enriched tissues are displayed in red (adjusted p-value ≤ 0.05), otherwise in blue. The tissues are ordered by their enrichment’s significance on up-regulated genes, and the graph also displays values for down-regulated genes as well as up/down-regulated.
Figure 4. Significantly enriched DEG (differentially expressed genes) for the 50 top genes on GTEx v8 54 tissue types. Log-transformed p-values for the enrichment of DEG in each tissue. Significantly enriched tissues are displayed in red (adjusted p-value ≤ 0.05), otherwise in blue. The tissues are ordered by their enrichment’s significance on up-regulated genes, and the graph also displays values for down-regulated genes as well as up/down-regulated.
Brainsci 10 00870 g004
Figure 5. Gene-ontology (GO) most significant categories for biological function. GO biological functions (at the left) are ordered by enrichment log-transformed p-value, displayed in blue. Only statistically significant results are shown (padj < 0.05). The proportion of overlapping genes (at the bottom) for each GO term (at the left) is displayed in red and for each term the specific genes are identified in yellow.
Figure 5. Gene-ontology (GO) most significant categories for biological function. GO biological functions (at the left) are ordered by enrichment log-transformed p-value, displayed in blue. Only statistically significant results are shown (padj < 0.05). The proportion of overlapping genes (at the bottom) for each GO term (at the left) is displayed in red and for each term the specific genes are identified in yellow.
Brainsci 10 00870 g005
Figure 6. STRING protein interaction analysis for the 50 top protein-coding genes. Inside each node, the 3D structure of each protein is displayed. The edges represent protein–protein interactions. Blue lines represent known interactions from curated databases and pink ones represent experimentally determined interactions. Black lines represent co-expression and purple ones represent homology.
Figure 6. STRING protein interaction analysis for the 50 top protein-coding genes. Inside each node, the 3D structure of each protein is displayed. The edges represent protein–protein interactions. Blue lines represent known interactions from curated databases and pink ones represent experimentally determined interactions. Black lines represent co-expression and purple ones represent homology.
Brainsci 10 00870 g006
Table 1. List of significant variants (rsID) at a p-value threshold of 5 × 10−6. Significant variants (rsID) and their respective p-value and regression coefficient (b1) for each phenotype (p_δ, p_θ, p_α, p_β1 and p_β2).
Table 1. List of significant variants (rsID) at a p-value threshold of 5 × 10−6. Significant variants (rsID) and their respective p-value and regression coefficient (b1) for each phenotype (p_δ, p_θ, p_α, p_β1 and p_β2).
rsIDpppp_β1p_β2b1
rs7516534 2.76 × 10−6 −0.08
rs12720066 1.39 × 10−6 0.10
rs12705973 2.65 × 10−6 0.02
rs10108126 1.97 × 10−6 0.05
rs10104429 1.55 × 10−7 0.05
rs7125249 3.70 × 10−6 0.03
rs6692346 2.32 × 10−60.04
rs4658030 2.04 × 10−6−0.03
rs77599684 9.17 × 10−70.05
rs6106856 2.21 × 10−60.03
rs71492954.88 × 10−6 0.09
rs12263011 4.98 × 10−6 0.07
rs71381191 2.64 × 10−8 1 0.11
rs12443654 4.19 × 10−6 0.07
rs9930193 4.57 × 10−6 0.07
rs55908084 2.65 × 10−6 0.06
rs9960516 2 3.09 × 10−6 0.05
rs72919581 2 1.82 × 10−6 0.06
rs1893824 1.66 × 10−6 −0.05
1 Values that passed the Bonferroni threshold of significance. 2 SNPs in LD (R2 = 0.73).
Table 2. Gene mapping of all brain wave associated SNPs with FUMA. Only SNPs showing maximum p-value of 0.01 in any of the five linear regressions were selected for this analysis. The table displays the strongly associated genomic loci (Locus), the p-value of the lead SNP for each locus (p), the individually significant SNPs for each locus (Sig_SNPs), the mapped genes (Gene), the gene product (type) and brain tissues from Genotype-Tissue Expression (GTEx) V8 mapped by gene expression-QTL (eqtlMapts) and by chromatin interaction (ciMapts) to the genomic loci.
Table 2. Gene mapping of all brain wave associated SNPs with FUMA. Only SNPs showing maximum p-value of 0.01 in any of the five linear regressions were selected for this analysis. The table displays the strongly associated genomic loci (Locus), the p-value of the lead SNP for each locus (p), the individually significant SNPs for each locus (Sig_SNPs), the mapped genes (Gene), the gene product (type) and brain tissues from Genotype-Tissue Expression (GTEx) V8 mapped by gene expression-QTL (eqtlMapts) and by chromatin interaction (ciMapts) to the genomic loci.
LocuspSig_SNPsGeneTypeeqtlMaptsciMapts
1:31318133:C:T2.76 × 10−6rs7516534LAPTM5p_coding Adult_Cortex
RN7SKP91misc
RNA
RP1-65J11.5antisense
1:192427329:C:G2.04 × 10−6rs4658030
7:87169702:A:C1.39 × 10−6rs12720066
7:114313199:A:G2.65 × 10−6rs12705973FOXP2p_codingBrain_Caudate_basal_gangliaFetal_Cortex:
Neural_Progen_Cell
AC073626.2antisense Fetal_Cortex:
Neural_Progen_Cell
MIR3666miRNA
8:74257947:A:G1.55 × 10−7rs10104429
rs10108126
10:13865505:A:C4.98 × 10−6rs12263011PRPF18p_coding Fetal_Cortex
FRMD4Ap_coding
CDNFp_coding Adult_Cortex:
Fetal_Cortex
HSPA14p_coding Adult_Cortex:
Fetal_Cortex
RP11-398C13.6lincRNA Fetal_Cortex
SUV39H2p_coding Fetal_Cortex
11:45062339:G:T3.70 × 10−6rs7125249
14:51283148:C:T4.88 × 10−6rs7149295L2HGDHp_coding Fetal_Cortex
ATP5Sp_coding Fetal_Cortex
NINp_coding Adult_Cortex:
Hippocampus
RP11-286O18.1antisense
PYGLp_coding Adult_Cortex
16:11156812:A:G2.64 × 10−8rs71381191CLEC16Ap_coding
RPL7P46pseudogene
RP11-66H6.3antisense
16:19375297:C:T4.19 × 10−6rs12443654RPS15Ap_codingBrain_Spinal_cord_cervical_c-1
rs9930193CTA-363E6.2lincRNA
17:65034162:C:G2.65 × 10−6rs55908084CACNG4p_coding
AC005544.1p_coding
RP11-74H8.1antisense
18:21849024:A:G9.17 × 10−7rs77599684OSBPL1Ap_coding
RN7SL247Pmisc_RNA
18:45885064:C:T1.82 × 10−6rs72919581ZBTB7Cp_coding
rs9960516
18:74959125:G:T1.66 × 10−6rs1893824RP11-17M16.2antisenseBrain_Hypothalamus
GALR1p_codingBrain_Anterior_cingulate_cortex_BA24
20:24313473:A:G2.21 × 10−6rs6106856
Lead SNPs are underlined, common findings between SNP-based and Gene-based approaches are bolded.
Table 3. Genescores for the top 50 most significant genes. The table displays the number of SNPs mapped to each gene (nSNPs) as well as the computed p-value for each gene.
Table 3. Genescores for the top 50 most significant genes. The table displays the number of SNPs mapped to each gene (nSNPs) as well as the computed p-value for each gene.
GenenSNPsp-Value
CSMD1353.01 × 10−16
CDH13274.85 × 10−14
PTPRD217.89 × 10−14
SORCS2211.23 × 10−13
RBFOX1202.05 × 10−13
MACROD2201.39 × 10−12
ZBTB7C122.29 × 10−12
RBFOX3186.64 × 10−12
RUNX1141.17 × 10−11
WWOX232.39 × 10−11
FRMD4A149.68 × 10−11
PLCB1131.08 × 10−10
LRP1B121.55 × 10−10
OPCML112.74 × 10−10
CAMK1D133.29 × 10−10
FHIT133.97 × 10−10
PDZD2134.84 × 10−10
CNTN4145.03 × 10−10
CDH2387.60 × 10−10
LPP108.21 × 10−10
MYO16118.79 × 10−10
CNTNAP2128.92 × 10−10
KDM4C121.06 × 10−9
NTM111.09 × 10−9
ASIC291.13 × 10−9
RYR3121.22 × 10−9
RPA3-AS1131.39 × 10−9
OFCC1101.41 × 10−9
DSCAM121.46 × 10−9
TMEM132C121.53 × 10−9
MTUS2111.98 × 10−9
SLC9A9102.12 × 10−9
LDLRAD4113.99 × 10−9
FSTL5114.06 × 10−9
RBMS3104.79 × 10−9
TENM4137.66 × 10−9
NELL1101.06 × 10−8
TENM381.06 × 10−8
PRKG1111.06 × 10−8
GRM781.19 × 10−8
HS3ST481.46 × 10−8
CNTN5101.46 × 10−8
TUSC3141.60 × 10−8
NRXN371.86 × 10−8
DPP6111.92 × 10−8
UNC5C92.06 × 10−8
MAGI2102.08 × 10−8
ADARB282.27 × 10−8
RIMBP272.53 × 10−8
GPC572.59 × 10−8
Genes identified by highly significant single SNPs are bolded.
Table 4. Enrichment on genesets of reported genes for brain conditions/phenotypes from the NHGRI (National Human Genome Research Institute) GWAS (genome-wide association studies) catalog.
Table 4. Enrichment on genesets of reported genes for brain conditions/phenotypes from the NHGRI (National Human Genome Research Institute) GWAS (genome-wide association studies) catalog.
GeneSetNnp-ValueAdjusted p-ValueGenes
Chronotype556137.05 × 10−136.39 × 10−10CNTN5, GPC5, MYO16, NRXN3, RBFOX1, ASIC2, MACROD2, CNTN4, GRM7, FSTL5, MAGI2, CSMD1, PTPRD
Intracranial aneurysm7268.85 × 10−103.97 × 10−7RBFOX1, DSCAM, RBMS3, FHIT, PDZD2, PTPRD
Response to amphetamines3351.09 × 10−93.97 × 10−7NELL1, WWOX, CDH13, FHIT, TENM3
Schizophrenia827121.43 × 10−94.31 × 10−7PRKG1, NELL1, NRXN3, RYR3, CDH13, CNTN4, RBMS3, FHIT, TENM3, MAGI2, CSMD1, KDM4C
Amyotrophic lateral sclerosis (sporadic)16473.58 × 10−98.76 × 10−7CNTN5, OPCML, RYR3, ASIC2, MACROD2, DSCAM, CSMD1
Night sleep phenotypes538103.86 × 10−98.76 × 10−7TENM4, CDH13, LRP1B, CNTN4, GRM7, RBMS3, SLC9A9, OFCC1, MAGI2, CSMD1
Cognitive ability, years of educational attainment or schizophrenia (pleiotropy)19763.72 × 10−74.22 × 10−5CAMK1D, NTM, GPC5, CDH13, LRP1B, CNTN4
Loneliness (multivariate analysis)2939.65 × 10−67.00 × 10−4PRKG1, CNTN5, PTPRD
Asperger disorder622.97 × 10−51.54 × 10−3NTM, FHIT
Middle childhood and early adolescence aggressive behavior622.97 × 10−51.54 × 10−3OPCML, CNTN4
Daytime sleep phenotypes25953.38 × 10−51.66 × 10−3WWOX, CDH13, PLCB1, CNTN4, GRM7
Brain connectivity724.15 × 10−51.88 × 10−3MACROD2, CNTN4
Short-term memory (digit-span task)724.15 × 10−51.88 × 10−3CDH13, PTPRD
Dimensional psychopathology (Negative)927.10 × 10−52.53 × 10−3CDH13, MACROD2
Hippocampal sclerosis927.10 × 10−52.53 × 10−3NELL1, SORCS2
Bipolar disorder (body mass index interaction)1028.86 × 10−52.87 × 10−3CDH23, WWOX
Aggressiveness in attention deficit hyperactivity disorder1121.08 × 10−43.33 × 10−3NTM, CSMD1
Alzheimer’s disease (64°)7031.39 × 10−43.93 × 10−3FRMD4A, MYO16, CNTNAP2
Enrichment on genesets of reported genes for brain conditions/phenotypes (GeneSet) from the NHGRI GWAS catalog [31] for the 50 most significant genes of our study. The table displays the background number of genes in each geneset (N), the sample number of genes present in each geneset (n), the respective p-value and adjusted p-value calculated from the comparison between the background and the sample frequency of genes and the sample genes (Genes).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rebelo, M.Â.; Gómez, C.; Gomes, I.; Poza, J.; Martins, S.; Maturana-Candelas, A.; Ruiz-Gómez, S.J.; Durães, L.; Sousa, P.; Figueruelo, M.; et al. Genome-Wide Scan for Five Brain Oscillatory Phenotypes Identifies a New QTL Associated with Theta EEG Band. Brain Sci. 2020, 10, 870. https://doi.org/10.3390/brainsci10110870

AMA Style

Rebelo MÂ, Gómez C, Gomes I, Poza J, Martins S, Maturana-Candelas A, Ruiz-Gómez SJ, Durães L, Sousa P, Figueruelo M, et al. Genome-Wide Scan for Five Brain Oscillatory Phenotypes Identifies a New QTL Associated with Theta EEG Band. Brain Sciences. 2020; 10(11):870. https://doi.org/10.3390/brainsci10110870

Chicago/Turabian Style

Rebelo, Miguel Ângelo, Carlos Gómez, Iva Gomes, Jesús Poza, Sandra Martins, Aarón Maturana-Candelas, Saúl J. Ruiz-Gómez, Luis Durães, Patrícia Sousa, Manuel Figueruelo, and et al. 2020. "Genome-Wide Scan for Five Brain Oscillatory Phenotypes Identifies a New QTL Associated with Theta EEG Band" Brain Sciences 10, no. 11: 870. https://doi.org/10.3390/brainsci10110870

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop