Migraine is a common neurological disorder with a genetically complex background. This paper describes a meta-analysis of genome-wide association (GWA) studies on migraine, performed by the Dutch-Icelandic migraine genetics (DICE) consortium, which brings together six population-based European migraine cohorts with a total sample size of 10,980 individuals (2446 cases and 8534 controls). A total of 32 SNPs showed marginal evidence for association at a P-value<10(-5). The best result was obtained for SNP rs9908234, which had a P-value of 8.00 × 10(-8). This top SNP is located in the nerve growth factor receptor (NGFR) gene. However, this SNP did not replicate in three cohorts from the Netherlands and Australia. Of the other 31 SNPs, 18 SNPs were tested in two replication cohorts, but none replicated. In addition, we explored previously identified candidate genes in the meta-analysis data set. This revealed a modest gene-based significant association between migraine and the metadherin (MTDH) gene, previously identified in the first clinic-based GWA study (GWAS) for migraine (Bonferroni-corrected gene-based P-value=0.026). This finding is consistent with the involvement of the glutamate pathway in migraine. Additional research is necessary to further confirm the involvement of glutamate.

Twin studies provide evidence that genetic influences contribute strongly to individual differences in exercise behavior. We hypothesize that part of this heritability is explained by genetic variation in the dopaminergic reward system. Eight single nucleotide polymorphisms (SNPs in DRD1: rs265981, DRD2: rs6275, rs1800497, DRD3: rs6280, DRD4: rs1800955, DBH: rs1611115, rs2519152, and in COMT: rs4680) and three variable number of tandem repeats (VNTRs in DRD4, upstream of DRD5, and in DAT1) were investigated for an association with regular leisure time exercise behavior.

Economic variables such as income, education, and occupation are known to affect mortality and morbidity, such as cardiovascular disease, and have also been shown to be partly heritable. However, very little is known about which genes influence economic variables, although these genes may have both a direct and an indirect effect on health. We report results from the first large-scale collaboration that studies the molecular genetic architecture of an economic variable-entrepreneurship-that was operationalized using self-employment, a widely-available proxy. Our results suggest that common SNPs when considered jointly explain about half of the narrow-sense heritability of self-employment estimated in twin data (σ(g)(2)/σ(P)(2) = 25%, h(2) = 55%). However, a meta-analysis of genome-wide association studies across sixteen studies comprising 50,627 participants did not identify genome-wide significant SNPs. 58 SNPs with p<10(-5) were tested in a replication sample (n = 3,271), but none replicated. Furthermore, a gene-based test shows that none of the genes that were previously suggested in the literature to influence entrepreneurship reveal significant associations. Finally, SNP-based genetic scores that use results from the meta-analysis capture less than 0.2% of the variance in self-employment in an independent sample (p≥0.039). Our results are consistent with a highly polygenic molecular genetic architecture of self-employment, with many genetic variants of small effect. Although self-employment is a multi-faceted, heavily environmentally influenced, and biologically distal trait, our results are similar to those for other genetically complex and biologically more proximate outcomes, such as height, intelligence, personality, and several diseases.

Glycated hemoglobin (HbA1c) is used to diagnose type 2 diabetes (T2D) and assess glycemic control in patients with diabetes. Previous genome-wide association studies (GWAS) have identified 18 HbA1c-associated genetic variants. These variants proved to be classifiable by their likely biological action as erythrocytic (also associated with erythrocyte traits) or glycemic (associated with other glucose-related traits). In this study, we tested the hypotheses that, in a very large scale GWAS, we would identify more genetic variants associated with HbA1c and that HbA1c variants implicated in erythrocytic biology would affect the diagnostic accuracy of HbA1c. We therefore expanded the number of HbA1c-associated loci and tested the effect of genetic risk-scores comprised of erythrocytic or glycemic variants on incident diabetes prediction and on prevalent diabetes screening performance. Throughout this multiancestry study, we kept a focus on interancestry differences in HbA1c genetics performance that might influence race-ancestry differences in health outcomes.

Telomere shortening has been associated with multiple age-related diseases such as cardiovascular disease, diabetes, and dementia. However, the biological mechanisms responsible for these associations remain largely unknown. In order to gain insight into the metabolic processes driving the association of leukocyte telomere length (LTL) with age-related diseases, we investigated the association between LTL and serum metabolite levels in 7,853 individuals from seven independent cohorts. LTL was determined by quantitative polymerase chain reaction and the levels of 131 serum metabolites were measured with mass spectrometry in biological samples from the same blood draw. With partial correlation analysis, we identified six metabolites that were significantly associated with LTL after adjustment for multiple testing: lysophosphatidylcholine acyl C17:0 (lysoPC a C17:0, p-value = 7.1 × 10-6), methionine (p-value = 9.2 × 10-5), tyrosine (p-value = 2.1 × 10-4), phosphatidylcholine diacyl C32:1 (PC aa C32:1, p-value = 2.4 × 10-4), hydroxypropionylcarnitine (C3-OH, p-value = 2.6 × 10-4), and phosphatidylcholine acyl-alkyl C38:4 (PC ae C38:4, p-value = 9.0 × 10-4). Pathway analysis showed that the three phosphatidylcholines and methionine are involved in homocysteine metabolism and we found supporting evidence for an association of lipid metabolism with LTL. In conclusion, we found longer LTL associated with higher levels of lysoPC a C17:0 and PC ae C38:4, and with lower levels of methionine, tyrosine, PC aa C32:1, and C3-OH. These metabolites have been implicated in inflammation, oxidative stress, homocysteine metabolism, and in cardiovascular disease and diabetes, two major drivers of morbidity and mortality.

We investigated whether obsessive-compulsive (OC) symptoms from a population-based sample could be analyzed to detect genetic variants influencing obsessive-compulsive disorder (OCD). We performed a genome-wide association studies (GWAS) on the obsession (rumination and impulsions) and compulsion (checking, washing, and ordering/precision) subscales of an abbreviated version of the Padua Inventory (N = 8,267 with genome-wide genotyping and phenotyping). The compulsion subscale showed a substantial and significant positive genetic correlation with an OCD case-control GWAS (r G = 0.61, p = .017) previously published by the Psychiatric Genomics Consortium (PGC-OCD). The obsession subscale and the total Padua score showed no significant genetic correlations (r G = -0.02 and r G = 0.42, respectively). A meta-analysis of the compulsive symptoms GWAS with the PGC-OCD revealed no genome-wide significant Single-Nucleotide Polymorphisms (SNPs combined N = 17,992, indicating that the power is still low for individual SNP effects). A gene-based association analysis, however, yielded two novel genes (WDR7 and ADCK1). The top 250 genes in the gene-based test also showed a significant increase in enrichment for psychiatric and brain-expressed genes. S-Predixcan testing showed that for genes expressed in hippocampus, amygdala, and caudate nucleus significance increased in the meta-analysis with compulsive symptoms compared to the original PGC-OCD GWAS. Thus, the inclusion of dimensional symptom data in genome-wide association on clinical case-control GWAS of OCD may be useful to find genes for OCD if the data are based on quantitative indices of compulsive behavior. SNP-level power increases were limited, but aggregate, gene-level analyses showed increased enrichment for brain-expressed genes related to psychiatric disorders, and increased association with gene expression in brain tissues with known emotional, reward processing, memory, and fear-formation functions.

Oscillatory activity is crucial for information processing in the brain, and has a long history as a biomarker for psychopathology. Variation in oscillatory activity is highly heritable, but current understanding of specific genetic influences remains limited. We performed the largest genome-wide association study to date of oscillatory power during eyes-closed resting electroencephalogram (EEG) across a range of frequencies (delta 1-3.75 Hz, theta 4-7.75 Hz, alpha 8-12.75 Hz, and beta 13-30 Hz) in 8,425 subjects. Additionally, we performed KGG positional gene-based analysis and brain-expression analyses. GABRA2-a known genetic marker for alcohol use disorder and epilepsy-significantly affected beta power, consistent with the known relation between GABAA interneuron activity and beta oscillations. Tissue-specific SNP-based imputation of gene-expression levels based on the GTEx database revealed that hippocampal GABRA2 expression may mediate this effect. Twenty-four genes at 3p21.1 were significant for alpha power (FDR q < .05). SNPs in this region were linked to expression of GLYCTK in hippocampal tissue, and GNL3 and ITIH4 in the frontal cortex-genes that were previously implicated in schizophrenia and bipolar disorder. In sum, we identified several novel genetic variants associated with oscillatory brain activity; furthermore, we replicated and advanced understanding of previously known genes associated with psychopathology (i.e., schizophrenia and alcohol use disorders). Importantly, these psychopathological liability genes affect brain functioning, linking the genes' expression to specific cortical/subcortical brain regions.

Epidemiological and genetic studies on COVID-19 are currently hindered by inconsistent and limited testing policies to confirm SARS-CoV-2 infection. Recently, it was shown that it is possible to predict COVID-19 cases using cross-sectional self-reported disease-related symptoms. Here, we demonstrate that this COVID-19 prediction model has reasonable and consistent performance across multiple independent cohorts and that our attempt to improve upon this model did not result in improved predictions. Using the existing COVID-19 prediction model, we then conducted a GWAS on the predicted phenotype using a total of 1,865 predicted cases and 29,174 controls. While we did not find any common, large-effect variants that reached genome-wide significance, we do observe suggestive genetic associations at two SNPs (rs11844522, p = 1.9x10-7; rs5798227, p = 2.2x10-7). Explorative analyses furthermore suggest that genetic variants associated with other viral infectious diseases do not overlap with COVID-19 susceptibility and that severity of COVID-19 may have a different genetic architecture compared to COVID-19 susceptibility. This study represents a first effort that uses a symptom-based predicted phenotype as a proxy for COVID-19 in our pursuit of understanding the genetic susceptibility of the disease. We conclude that the inclusion of symptom-based predicted cases could be a useful strategy in a scenario of limited testing, either during the current COVID-19 pandemic or any future viral outbreak.

Risky behaviors, such as substance use and unprotected sex, are associated with various physical and mental health problems. Recent genome-wide association studies indicated that variation in the cell adhesion molecule 2 (CADM2) gene plays a role in risky behaviors and self-control. In this phenome-wide scan for risky behavior, it was tested if underlying common vulnerability could be (partly) explained by pleiotropic effects of this gene and how large the effects were. Single nucleotide polymorphism (SNP)-level and gene-level association tests within four samples (25 and Up, Spit for Science, Netherlands Twin Register, and UK Biobank and meta-analyses over all samples (combined sample of 362,018 participants) were conducted to test associations between CADM2, substance- and sex-related risk behaviors, and various measures related to self-control. We found significant associations between the CADM2 gene, various risky behaviors, and different measures of self-control. The largest effect sizes were found for cannabis use, sensation seeking, and disinhibition. Effect sizes ranged from 0.01% to 0.26% for single top SNPs and from 0.07% to 3.02% for independent top SNPs together, with sufficient power observed only in the larger samples and meta-analyses. In the largest cohort, we found indications that risk-taking proneness mediated the association between CADM2 and latent factors for lifetime smoking and regular alcohol use. This study extends earlier findings that CADM2 plays a role in risky behaviors and self-control. It also provides insight into gene-level effect sizes and demonstrates the feasibility of testing mediation. These findings present a good starting point for investigating biological etiological pathways underlying risky behaviors.

Resilience and well-being are strongly related. People with higher levels of well-being are more resilient after stressful life events or trauma and vice versa. Less is known about the underlying sources of overlap and causality between the constructs. In a sample of 11.304 twins and 2.572 siblings from the Netherlands Twin Register, we investigated the overlap and possible direction of causation between resilience (i.e. the absence of psychiatric symptoms despite negative life events) and well-being (i.e. satisfaction with life) using polygenic score (PGS) prediction, twin-sibling modelling, and the Mendelian Randomization Direction of Causality (MR-DoC) model. Longitudinal twin-sibling models showed significant phenotypic correlations between resilience and well-being (.41/.51 at time 1 and 2). Well-being PGS were predictive for both well-being and resilience, indicating that genetic factors influencing well-being also predict resilience. Twin-sibling modeling confirmed this genetic correlation (0.71) and showed a strong environmental correlation (0.93). In line with causality, both genetic (51%) and environmental (49%) factors contributed significantly to the covariance between resilience and well-being. Furthermore, the results of within-subject and MZ twin differences analyses were in line with bidirectional causality. Additionally, we used the MR-DoC model combining both molecular and twin data to test causality, while correcting for pleiotropy. We confirmed the causal effect from well-being to resilience, with the direct effect of well-being explaining 11% (T1) and 20% (T2) of the variance in resilience. Data limitations prevented us to test the directional effect from resilience to well-being with the MR-DoC model. To conclude, we showed a strong relation between well-being and resilience. A first attempt to quantify the direction of this relationship points towards a bidirectional causal effect. If replicated, the potential mutual effects can have implications for interventions to lower psychopathology vulnerability, as resilience and well-being are both negatively related to psychopathology.

Low-frequency 1q21.1 distal deletion and duplication copy number variant (CNV) carriers are predisposed to multiple neurodevelopmental disorders, including schizophrenia, autism and intellectual disability. Human carriers display a high prevalence of micro- and macrocephaly in deletion and duplication carriers, respectively. The underlying brain structural diversity remains largely unknown. We systematically called CNVs in 38 cohorts from the large-scale ENIGMA-CNV collaboration and the UK Biobank and identified 28 1q21.1 distal deletion and 22 duplication carriers and 37,088 non-carriers (48% male) derived from 15 distinct magnetic resonance imaging scanner sites. With standardized methods, we compared subcortical and cortical brain measures (all) and cognitive performance (UK Biobank only) between carrier groups also testing for mediation of brain structure on cognition. We identified positive dosage effects of copy number on intracranial volume (ICV) and total cortical surface area, with the largest effects in frontal and cingulate cortices, and negative dosage effects on caudate and hippocampal volumes. The carriers displayed distinct cognitive deficit profiles in cognitive tasks from the UK Biobank with intermediate decreases in duplication carriers and somewhat larger in deletion carriers-the latter potentially mediated by ICV or cortical surface area. These results shed light on pathobiological mechanisms of neurodevelopmental disorders, by demonstrating gene dose effect on specific brain structures and effect on cognitive function.

We performed an EEG graph analysis on data from 31 typical readers (22.27 ± 2.53 y/o) and 24 dyslexics (22.99 ± 2.29 y/o), recorded while they were engaged in an audiovisual task and during resting-state. The task simulates reading acquisition as participants learned new letter-sound mappings via feedback. EEG data was filtered for the delta (0.5-4 Hz), theta (4-8 Hz), alpha (8-13 Hz), and beta (13-30 Hz) bands. We computed the Phase Lag Index (PLI) to provide an estimate of the functional connectivity between all pairs of electrodes per band. Then, networks were constructed using a Minimum Spanning Tree (MST), a unique sub-graph connecting all nodes (electrodes) without loops, aimed at minimizing bias in between groups and conditions comparisons. Both groups showed a comparable accuracy increase during task blocks, indicating that they correctly learned the new associations. The EEG results revealed lower task-specific theta connectivity, and lower theta degree correlation over both rest and task recordings, indicating less network integration in dyslexics compared to typical readers. This pattern suggests a role of theta oscillations in dyslexia and may reflect differences in task engagement between the groups, although robust correlations between MST metrics and performance indices were lacking.

Resting heart rate is associated with cardiovascular diseases and mortality in observational and Mendelian randomization studies. The aims of this study are to extend the number of resting heart rate associated genetic variants and to obtain further insights in resting heart rate biology and its clinical consequences. A genome-wide meta-analysis of 100 studies in up to 835,465 individuals reveals 493 independent genetic variants in 352 loci, including 68 genetic variants outside previously identified resting heart rate associated loci. We prioritize 670 genes and in silico annotations point to their enrichment in cardiomyocytes and provide insights in their ECG signature. Two-sample Mendelian randomization analyses indicate that higher genetically predicted resting heart rate increases risk of dilated cardiomyopathy, but decreases risk of developing atrial fibrillation, ischemic stroke, and cardio-embolic stroke. We do not find evidence for a linear or non-linear genetic association between resting heart rate and all-cause mortality in contrast to our previous Mendelian randomization study. Systematic alteration of key differences between the current and previous Mendelian randomization study indicates that the most likely cause of the discrepancy between these studies arises from false positive findings in previous one-sample MR analyses caused by weak-instrument bias at lower P-value thresholds. The results extend our understanding of resting heart rate biology and give additional insights in its role in cardiovascular disease development.

Through genome-wide association meta-analyses of up to 133,010 individuals of European ancestry without diabetes, including individuals newly genotyped using the Metabochip, we have increased the number of confirmed loci influencing glycemic traits to 53, of which 33 also increase type 2 diabetes risk (q < 0.05). Loci influencing fasting insulin concentration showed association with lipid levels and fat distribution, suggesting impact on insulin resistance. Gene-based analyses identified further biologically plausible loci, suggesting that additional loci beyond those reaching genome-wide significance are likely to represent real associations. This conclusion is supported by an excess of directionally consistent and nominally significant signals between discovery and follow-up studies. Functional analysis of these newly discovered loci will further improve our understanding of glycemic control.

The incretin hormone glucagon-like peptide 1 (GLP-1) promotes glucose homeostasis and enhances β-cell function. GLP-1 receptor agonists (GLP-1 RAs) and dipeptidyl peptidase-4 (DPP-4) inhibitors, which inhibit the physiological inactivation of endogenous GLP-1, are used for the treatment of type 2 diabetes. Using the Metabochip, we identified three novel genetic loci with large effects (30-40%) on GLP-1-stimulated insulin secretion during hyperglycemic clamps in nondiabetic Caucasian individuals (TMEM114; CHST3 and CTRB1/2; n = 232; all P ≤ 8.8 × 10(-7)). rs7202877 near CTRB1/2, a known diabetes risk locus, also associated with an absolute 0.51 ± 0.16% (5.6 ± 1.7 mmol/mol) lower A1C response to DPP-4 inhibitor treatment in G-allele carriers, but there was no effect on GLP-1 RA treatment in type 2 diabetic patients (n = 527). Furthermore, in pancreatic tissue, we show that rs7202877 acts as expression quantitative trait locus for CTRB1 and CTRB2, encoding chymotrypsinogen, and increases fecal chymotrypsin activity in healthy carriers. Chymotrypsin is one of the most abundant digestive enzymes in the gut where it cleaves food proteins into smaller peptide fragments. Our data identify chymotrypsin in the regulation of the incretin pathway, development of diabetes, and response to DPP-4 inhibitor treatment.

The degree to which genetic factors influence brain connectivity is beginning to be understood. Large-scale efforts are underway to map the profile of genetic effects in various brain regions. The NIH-funded Human Connectome Project (HCP) is providing data valuable for analyzing the degree of genetic influence underlying brain connectivity revealed by state-of-the-art neuroimaging methods. We calculated the heritability of the fractional anisotropy (FA) measure derived from diffusion tensor imaging (DTI) reconstruction in 481 HCP subjects (194/287 M/F) consisting of 57/60 pairs of mono- and dizygotic twins, and 246 siblings. FA measurements were derived using (Enhancing NeuroImaging Genetics through Meta-Analysis) ENIGMA DTI protocols and heritability estimates were calculated using the SOLAR-Eclipse imaging genetic analysis package. We compared heritability estimates derived from HCP data to those publicly available through the ENIGMA-DTI consortium, which were pooled together from five-family based studies across the US, Europe, and Australia. FA measurements from the HCP cohort for eleven major white matter tracts were highly heritable (h(2)=0.53-0.90, p<10(-5)), and were significantly correlated with the joint-analytical estimates from the ENIGMA cohort on the tract and voxel-wise levels. The similarity in regional heritability suggests that the additive genetic contribution to white matter microstructure is consistent across populations and imaging acquisition parameters. It also suggests that the overarching genetic influence provides an opportunity to define a common genetic search space for future gene-discovery studies. Uniquely, the measurements of additive genetic contribution performed in this study can be repeated using online genetic analysis tools provided by the HCP ConnectomeDB web application.

Heart rate variability is an important risk factor for cardiovascular disease and all-cause mortality. The acetylcholine pathway plays a key role in explaining heart rate variability in humans. We assessed whether 443 genotyped and imputed common genetic variants in eight key genes (CHAT, SLC18A3, SLC5A7, CHRNB4, CHRNA3, CHRNA, CHRM2 and ACHE) of the acetylcholine pathway were associated with variation in an established measure of heart rate variability reflecting parasympathetic control of the heart rhythm, the root mean square of successive differences (RMSSD) of normal RR intervals. The association was studied in a two stage design in individuals of European descent. First, analyses were performed in a discovery sample of four cohorts (n = 3429, discovery stage). Second, findings were replicated in three independent cohorts (n = 3311, replication stage), and finally the two stages were combined in a meta-analysis (n = 6740). RMSSD data were obtained under resting conditions. After correction for multiple testing, none of the SNPs showed an association with RMSSD. In conclusion, no common genetic variants for heart rate variability were identified in the largest and most comprehensive candidate gene study on the acetylcholine pathway to date. Future gene finding efforts for RMSSD may want to focus on hypothesis free approaches such as the genome-wide association study.

Reduced cardiac vagal control reflected in low heart rate variability (HRV) is associated with greater risks for cardiac morbidity and mortality. In two-stage meta-analyses of genome-wide association studies for three HRV traits in up to 53,174 individuals of European ancestry, we detect 17 genome-wide significant SNPs in eight loci. HRV SNPs tag non-synonymous SNPs (in NDUFA11 and KIAA1755), expression quantitative trait loci (eQTLs) (influencing GNG11, RGS6 and NEO1), or are located in genes preferentially expressed in the sinoatrial node (GNG11, RGS6 and HCN4). Genetic risk scores account for 0.9 to 2.6% of the HRV variance. Significant genetic correlation is found for HRV with heart rate (-0.74

Circulating metabolites have been shown to reflect metabolic changes during the development of type 2 diabetes. In this study we examined the association of metabolite levels and pairwise metabolite ratios with insulin responses after glucose, glucagon-like peptide-1 (GLP-1) and arginine stimulation. We then investigated if the identified metabolite ratios were associated with measures of OGTT-derived beta cell function and with prevalent and incident type 2 diabetes.

Objectives: Cardiac repolarization may be affected by psychiatric disorders and/or antidepressant use, but evidence for this is inconclusive. This study examined the relationship between depressive and anxiety disorder and use of antidepressants with T-wave amplitude (TWA) and QT-interval. Methods: Data was obtained from the Netherlands Study of Depression and Anxiety (n = 1,383). Depression/anxiety was diagnosed with the DSM-IV based Composite International Diagnostic Interview. The use of tricyclic antidepressants (TCAs), selective serotonin and noradrenalin reuptake inhibitors (SNRIs), and selective serotonin reuptake inhibitors (SSRIs) was established. T-wave amplitude and QT-interval corrected for heart rate (QTc) were obtained from an ECG measured in a type II axis configuration. Results: Compared to controls, persons with depression or anxiety disorders did not show a significantly different TWA (p = 0.58; Cohen's d = 0.046) or QTc (p = 0.48; Cohen's d = -0.057). In spite of known sympathomimetic effects, TCA use (p = 0.26; Cohen's d = -0.162) and SNRI use (p = 0.70; Cohen's d = -0.055) were not significantly associated with a lower TWA. TCA use (p = 0.12; Cohen's d = 0.225) and SNRI use (p = 0.11; Cohen's d = 0.227) were also not significantly associated with a prolonged QTc. Conclusion: We did not find evidence that either depressive/anxiety disorder or antidepressant use is associated with abnormalities in TWA or QTc. Earlier found sympathomimetic effects of TCAs and SNRIs are not evident in these measures of cardiac repolarization.