Genome-wide Association Analysis in Humans Links Nucleotide Metabolism to Leukocyte Telomere Length

Leukocyte telomere length (LTL) is a heritable biomarker of genomic aging. In this study, we perform a genome-wide meta-analysis of LTL by pooling densely genotyped and imputed association results across large-scale European-descent studies including up to 78,592 individuals. We identify 49 genomic regions at a false dicovery rate (FDR) < 0.05 threshold and prioritize genes at 31, with five highlighting nucleotide metabolism as an important regulator of LTL. We report six genome-wide significant loci in or near SENP7, MOB1B, CARMIL1, PRRC2A, TERF2, and RFWD3, and our results support recently identified PARP1, POT1, ATM, and MPHOSPH6 loci. Phenome-wide analyses in >350,000 UK Biobank participants suggest that genetically shorter telomere length increases the risk of hypothyroidism and decreases the risk of thyroid cancer, lymphoma, and a range of proliferative conditions. Our results replicate previously reported associations with increased risk of coronary artery disease and lower risk for multiple cancer types. Our findings substantially expand current knowledge on genes that regulate LTL and their impact on human health and disease.

Dose response of the 16p11.2 distal copy number variant on intracranial volume and basal ganglia

Carriers of large recurrent copy number variants (CNVs) have a higher risk of developing neurodevelopmental disorders. The 16p11.2 distal CNV predisposes carriers to e.g., autism spectrum disorder and schizophrenia. We compared subcortical brain volumes of 12 16p11.2 distal deletion and 12 duplication carriers to 6882 non-carriers from the large-scale brain Magnetic Resonance Imaging collaboration, ENIGMA-CNV. After stringent CNV calling procedures, and standardized FreeSurfer image analysis, we found negative dose-response associations with copy number on intracranial volume and on regional caudate, pallidum and putamen volumes (β = -0.71 to -1.37; P < 0.0005). In an independent sample, consistent results were obtained, with significant effects in the pallidum (β = -0.95, P = 0.0042). The two data sets combined showed significant negative dose-response for the accumbens, caudate, pallidum, putamen and ICV (P = 0.0032, 8.9 × 10-6, 1.7 × 10-9, 3.5 × 10-12 and 1.0 × 10-4, respectively). Full scale IQ was lower in both deletion and duplication carriers compared to non-carriers. This is the first brain MRI study of the impact of the 16p11.2 distal CNV, and we demonstrate a specific effect on subcortical brain structures, suggesting a neuropathological pattern underlying the neurodevelopmental syndromes.

Measurement and genetic architecture of lifetime depression in the Netherlands as assessed by LIDAS (Lifetime Depression Assessment Self-report)

BACKGROUND

Major depressive disorder (MDD) is a common mood disorder, with a heritability of around 34%. Molecular genetic studies made significant progress and identified genetic markers associated with the risk of MDD; however, progress is slowed down by substantial heterogeneity as MDD is assessed differently across international cohorts. Here, we used a standardized online approach to measure MDD in multiple cohorts in the Netherlands and evaluated whether this approach can be used in epidemiological and genetic association studies of depression.

METHODS

Within the Biobank Netherlands Internet Collaboration (BIONIC) project, we collected MDD data in eight cohorts involving 31 936 participants, using the online Lifetime Depression Assessment Self-report (LIDAS), and estimated the prevalence of current and lifetime MDD in 22 623 unrelated individuals. In a large Netherlands Twin Register (NTR) twin-family dataset (n ≈ 18 000), we estimated the heritability of MDD, and the prediction of MDD in a subset (n = 4782) through Polygenic Risk Score (PRS).

RESULTS

Estimates of current and lifetime MDD prevalence were 6.7% and 18.1%, respectively, in line with population estimates based on validated psychiatric interviews. In the NTR heritability estimates were 0.34/0.30 (s.e. = 0.02/0.02) for current/lifetime MDD, respectively, showing that the LIDAS gives similar heritability rates for MDD as reported in the literature. The PRS predicted risk of MDD (OR 1.23, 95% CI 1.15-1.32, R2 = 1.47%).

CONCLUSIONS

By assessing MDD status in the Netherlands using the LIDAS instrument, we were able to confirm previously reported MDD prevalence and heritability estimates, which suggests that this instrument can be used in epidemiological and genetic association studies of depression.

Genetic architecture of subcortical brain structures in 38,851 individuals

Subcortical brain structures are integral to motion, consciousness, emotions and learning. We identified common genetic variation related to the volumes of the nucleus accumbens, amygdala, brainstem, caudate nucleus, globus pallidus, putamen and thalamus, using genome-wide association analyses in almost 40,000 individuals from CHARGE, ENIGMA and UK Biobank. We show that variability in subcortical volumes is heritable, and identify 48 significantly associated loci (40 novel at the time of analysis). Annotation of these loci by utilizing gene expression, methylation and neuropathological data identified 199 genes putatively implicated in neurodevelopment, synaptic signaling, axonal transport, apoptosis, inflammation/infection and susceptibility to neurological disorders. This set of genes is significantly enriched for Drosophila orthologs associated with neurodevelopmental phenotypes, suggesting evolutionarily conserved mechanisms. Our findings uncover novel biology and potential drug targets underlying brain development and disease.

A characterization of cis- and trans-heritability of RNA-Seq-based gene expression

Insights into individual differences in gene expression and its heritability (h²) can help in understanding pathways from DNA to phenotype. We estimated the heritability of gene expression of 52,844 genes measured in whole blood in the largest twin RNA-Seq sample to date (1497 individuals including 459 monozygotic twin pairs and 150 dizygotic twin pairs) from classical twin modeling and identity-by-state-based approaches. We estimated for each gene h²_total, composed of cis-heritability (h²_cis, the variance explained by single nucleotide polymorphisms in the cis-window of the gene), and trans-heritability (h²_res, the residual variance explained by all other genome-wide variants). Mean h²_total was 0.26, which was significantly higher than heritability estimates earlier found in a microarray-based study using largely overlapping (>60%) RNA samples (mean h² = 0.14, p = 6.15 × 10^-258). Mean h²_cis was 0.06 and strongly correlated with beta of the top cis expression quantitative loci (eQTL, ρ = 0.76, p < 10^-308) and with estimates from earlier RNA-Seq-based studies. Mean h²_res was 0.20 and correlated with the beta of the corresponding trans-eQTL (ρ = 0.04, p < 1.89 × 10^-3) and was significantly higher for genes involved in cytokine-cytokine interactions (p = 4.22 × 10^-15), many other immune system pathways, and genes identified in genome-wide association studies for various traits including behavioral disorders and cancer. This study provides a thorough characterization of cis- and trans-h² estimates of gene expression, which is of value for interpretation of GWAS and gene expression studies.

Testing Familial Transmission of Smoking With Two Different Research Designs

Introduction

Classical twin studies show that smoking is heritable. To determine if shared family environment plays a role in addition to genetic factors, and if they interact (G×E), we use a children-of-twins design. In a second sample, we measure genetic influence with polygenic risk scores (PRS) and environmental influence with a question on exposure to smoking during childhood.

Methods

Data on smoking initiation were available for 723 children of 712 twins from the Netherlands Twin Register (64.9% female, median birth year 1985). Children were grouped in ascending order of risk, based on smoking status and zygosity of their twin-parent and his/her co-twin: never smoking twin-parent with a never smoking co-twin; never smoking twin-parent with a smoking dizygotic co-twin; never smoking twin-parent with a smoking monozygotic co-twin; and smoking twin-parent with a smoking or never smoking co-twin. For 4072 participants from the Netherlands Twin Register (67.3% female, median birth year 1973), PRS for smoking were computed and smoking initiation, smoking heaviness, and exposure to smoking during childhood were available.

Results

Patterns of smoking initiation in the four group children-of-twins design suggested shared familial influences in addition to genetic factors. PRS for ever smoking were associated with smoking initiation in all individuals. PRS for smoking heaviness were associated with smoking heaviness in individuals exposed to smoking during childhood, but not in non-exposed individuals.

Conclusions

Shared family environment influences smoking, over and above genetic factors. Genetic risk of smoking heaviness was only important for individuals exposed to smoking during childhood, versus those not exposed (G×E).

Implications

This study adds to the very few existing children-of-twins (CoT) studies on smoking and combines a CoT design with a second research design that utilizes polygenic risk scores and data on exposure to smoking during childhood. The results show that shared family environment affects smoking behavior over and above genetic factors. There was also evidence for gene-environment interaction (G×E) such that genetic risk of heavy versus light smoking was only important for individuals who were also exposed to (second-hand) smoking during childhood. Together, these findings give additional incentive to recommending parents not to expose their children to cigarette smoking.

Genome-wide Burden of Rare Short Deletions Is Enriched in Major Depressive Disorder in Four Cohorts

BACKGROUND

Major depressive disorder (MDD) is moderately heritable, with a high prevalence and a presumed high heterogeneity. Copy number variants (CNVs) could contribute to the heritable component of risk, but the two previous genome-wide association studies of rare CNVs did not report significant findings.

METHODS

In this meta-analysis of four cohorts (5780 patients and 6626 control subjects), we analyzed the association of MDD to 1) genome-wide burden of rare deletions and duplications, partitioned by length (<100 kb or >100 kb) and other characteristics, and 2) individual rare exonic CNVs and CNV regions.

RESULTS

Patients with MDD carried significantly more short deletions than control subjects (p = .0059) but not long deletions or short or long duplications. The confidence interval for long deletions overlapped with that for short deletions, but long deletions were 70% less frequent genome-wide, reducing the power to detect increased burden. The increased burden of short deletions was primarily in intergenic regions. Short deletions in cases were also modestly enriched for high-confidence enhancer regions. No individual CNV achieved thresholds for suggestive or significant association after genome-wide correction. p values < .01 were observed for 15q11.2 duplications (TUBGCP5, CYFIP1, NIPA1, and NIPA2), deletions in or near PRKN or MSR1, and exonic duplications of ATG5.

CONCLUSIONS

The increased burden of short deletions in patients with MDD suggests that rare CNVs increase the risk of MDD by disrupting regulatory regions. Results for longer deletions were less clear, but no large effects were observed for long multigenic CNVs (as seen in schizophrenia and autism). Further studies with larger sample sizes are warranted.

Genome-wide association meta-analysis of 30,000 samples identifies seven novel loci for quantitative ECG traits

Genome-wide association studies (GWAS) of quantitative electrocardiographic (ECG) traits in large consortia have identified more than 130 loci associated with QT interval, QRS duration, PR interval, and heart rate (RR interval). In the current study, we meta-analyzed genome-wide association results from 30,000 mostly Dutch samples on four ECG traits: PR interval, QRS duration, QT interval, and RR interval. SNP genotype data was imputed using the Genome of the Netherlands reference panel encompassing 19 million SNPs, including millions of rare SNPs (minor allele frequency < 5%). In addition to many known loci, we identified seven novel locus-trait associations: KCND3, NR3C1, and PLN for PR interval, KCNE1, SGIP1, and NFKB1 for QT interval, and ATP2A2 for QRS duration, of which six were successfully replicated. At these seven loci, we performed conditional analyses and annotated significant SNPs (in exons and regulatory regions), demonstrating involvement of cardiac-related pathways and regulation of nearby genes.

Maternal and fetal genetic effects on birth weight and their relevance to cardio-metabolic risk factors

Birth weight variation is influenced by fetal and maternal genetic and non-genetic factors, and has been reproducibly associated with future cardio-metabolic health outcomes. In expanded genome-wide association analyses of own birth weight (n = 321,223) and offspring birth weight (n = 230,069 mothers), we identified 190 independent association signals (129 of which are novel). We used structural equation modeling to decompose the contributions of direct fetal and indirect maternal genetic effects, then applied Mendelian randomization to illuminate causal pathways. For example, both indirect maternal and direct fetal genetic effects drive the observational relationship between lower birth weight and higher later blood pressure: maternal blood pressure-raising alleles reduce offspring birth weight, but only direct fetal effects of these alleles, once inherited, increase later offspring blood pressure. Using maternal birth weight-lowering genotypes to proxy for an adverse intrauterine environment provided no evidence that it causally raises offspring blood pressure, indicating that the inverse birth weight-blood pressure association is attributable to genetic effects, and not to intrauterine programming.

A Genetic Investigation of the Well-Being Spectrum

The interrelations among well-being, neuroticism, and depression can be captured in a so-called well-being spectrum (3-phenotype well-being spectrum, 3-WBS). Several other human traits are likely linked to the 3-WBS. In the present study, we investigate how the 3-WBS can be expanded. First, we constructed polygenic risk scores for the 3-WBS and used this score to predict a series of traits that have been associated with well-being in the literature. We included information on loneliness, big five personality traits, self-rated health, and flourishing. The 3-WBS polygenic score predicted all the original 3-WBS traits and additionally loneliness, self-rated health, and extraversion (R² between 0.62% and 1.58%). Next, using LD score regression, we calculated genetic correlations between the 3-WBS and the traits of interest. From all candidate traits, loneliness and self-rated health were found to have the strongest genetic correlations (r_g = - 0.79, and r_g= 0.64, respectively) with the 3-WBS. Lastly, we use Genomic SEM to investigate the factor structure of the proposed spectrum. The best model fit was obtained for a two-factor model including the 5-WBS traits, with two highly correlated factors representing the negative- and positive end of the spectrum. Based on these analyses we propose to include loneliness and self-rated health in the WBS and use a 5-phenotype well-being spectrum in future studies to gain more insight into the determinants of human well-being.

Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use

Tobacco and alcohol use are leading causes of mortality that influence risk for many complex diseases and disorders1. They are heritable2,3 and etiologically related4,5 behaviors that have been resistant to gene discovery efforts6-11. In sample sizes up to 1.2 million individuals, we discovered 566 genetic variants in 406 loci associated with multiple stages of tobacco use (initiation, cessation, and heaviness) as well as alcohol use, with 150 loci evidencing pleiotropic association. Smoking phenotypes were positively genetically correlated with many health conditions, whereas alcohol use was negatively correlated with these conditions, such that increased genetic risk for alcohol use is associated with lower disease risk. We report evidence for the involvement of many systems in tobacco and alcohol use, including genes involved in nicotinic, dopaminergic, and glutamatergic neurotransmission. The results provide a solid starting point to evaluate the effects of these loci in model organisms and more precise substance use measures.

A Potential Role for the STXBP5-AS1 Gene in Adult ADHD Symptoms

We aimed to detect Attention-deficit/hyperactivity (ADHD) risk-conferring genes in adults. In children, ADHD is characterized by age-inappropriate levels of inattention and/or hyperactivity-impulsivity and may persists into adulthood. Childhood and adulthood ADHD are heritable, and are thought to represent the clinical extreme of a continuous distribution of ADHD symptoms in the general population. We aimed to leverage the power of studies of quantitative ADHD symptoms in adults who were genotyped. Within the SAGA (Study of ADHD trait genetics in adults) consortium, we estimated the single nucleotide polymorphism (SNP)-based heritability of quantitative self-reported ADHD symptoms and carried out a genome-wide association meta-analysis in nine adult population-based and case-only cohorts of adults. A total of n = 14,689 individuals were included. In two of the SAGA cohorts we found a significant SNP-based heritability for self-rated ADHD symptom scores of respectively 15% (n = 3656) and 30% (n = 1841). The top hit of the genome-wide meta-analysis (SNP rs12661753; p-value = 3.02 × 10^-7) was present in the long non-coding RNA gene STXBP5-AS1. This association was also observed in a meta-analysis of childhood ADHD symptom scores in eight population-based pediatric cohorts from the Early Genetics and Lifecourse Epidemiology (EAGLE) ADHD consortium (n = 14,776). Genome-wide meta-analysis of the SAGA and EAGLE data (n = 29,465) increased the strength of the association with the SNP rs12661753. In human HEK293 cells, expression of STXBP5-AS1 enhanced the expression of a reporter construct of STXBP5, a gene known to be involved in "SNAP" (Soluble NSF attachment protein) Receptor" (SNARE) complex formation. In mouse strains featuring different levels of impulsivity, transcript levels in the prefrontal cortex of the mouse ortholog Gm28905 strongly correlated negatively with motor impulsivity as measured in the five choice serial reaction time task (r² = - 0.61; p = 0.004). Our results are consistent with an effect of the STXBP5-AS1 gene on ADHD symptom scores distribution and point to a possible biological mechanism, other than antisense RNA inhibition, involved in ADHD-related impulsivity levels.

Dopaminergic Genetic Variants and Voluntary Externally Paced Exercise Behavior

PURPOSE

Most candidate gene studies on the neurobiology of voluntary exercise behavior have focused on the dopaminergic signaling pathway and its role in the mesolimbic reward system. We hypothesized that dopaminergic candidate genes may influence exercise behavior through additional effects on executive functioning and that these effects are only detected when the types of exercise activity are taken into account.

METHODS

Data on voluntary exercise behavior and at least one single-nucleotide polymorphism/variable number of tandem repeat (VNTR) were available for 12,929 participants of the Netherlands Twin Registry. Exercise activity was classified as externally paced if a high level of executive function skill was required. The total volume of voluntary exercise (minutes per week) as well as the volume specifically spent on externally paced activities were tested for association with nine functional dopaminergic polymorphisms (DRD1: rs265981, DRD2/ANKK1: rs1800497, DRD3: rs6280, DRD4: VNTR 48 bp, DRD5: VNTR 130-166 bp, DBH: rs2519152, DAT1: VNTR 40 bp, COMT: rs4680, MAOA: VNTR 30 bp), a polygenic score (PGS) based on nine alleles leading to lower dopamine responsiveness, and a PGS based on three alleles associated with both higher reward sensitivity and better executive functioning (DRD2/ANKK1: "G" allele, COMT: Met allele, DAT1: 440-bp allele).

RESULTS

No association with total exercise volume or externally paced exercise volume was found for individual alleles or the nine-allele PGS. The volume of externally paced exercise behavior was significantly associated with the reward and executive function congruent PGS. This association was driven by the DAT1 440-bp and COMT Met allele, which acted as increaser alleles for externally paced exercise behavior.

CONCLUSIONS

Taking into account the types of exercise activity may increase the success of identifying genetic variants and unraveling the neurobiology of voluntary exercise behavior.

Multivariate genome-wide analyses of the well-being spectrum

We introduce two novel methods for multivariate genome-wide-association meta-analysis (GWAMA) of related traits that correct for sample overlap. A broad range of simulation scenarios supports the added value of our multivariate methods relative to univariate GWAMA. We applied the novel methods to life satisfaction, positive affect, neuroticism, and depressive symptoms, collectively referred to as the well-being spectrum (N_obs = 2,370,390), and found 304 significant independent signals. Our multivariate approaches resulted in a 26% increase in the number of independent signals relative to the four univariate GWAMAs and in an ~57% increase in the predictive power of polygenic risk scores. Supporting transcriptome- and methylome-wide analyses (TWAS and MWAS, respectively) uncovered an additional 17 and 75 independent loci, respectively. Bioinformatic analyses, based on gene expression in brain tissues and cells, showed that genes differentially expressed in the subiculum and GABAergic interneurons are enriched in their effect on the well-being spectrum.

Genome-wide association study of offspring birth weight in 86 577 women identifies five novel loci and highlights maternal genetic effects that are independent of fetal genetics

Genome-wide association studies of birth weight have focused on fetal genetics, whereas relatively little is known about the role of maternal genetic variation. We aimed to identify maternal genetic variants associated with birth weight that could highlight potentially relevant maternal determinants of fetal growth. We meta-analysed data on up to 8.7 million SNPs in up to 86 577 women of European descent from the Early Growth Genetics (EGG) Consortium and the UK Biobank. We used structural equation modelling (SEM) and analyses of mother–child pairs to quantify the separate maternal and fetal genetic effects. Maternal SNPs at 10 loci (MTNR1B, HMGA2, SH2B3, KCNAB1, L3MBTL3, GCK, EBF1, TCF7L2, ACTL9, CYP3A7) were associated with offspring birth weight at P < 5 × 10⁻⁸. In SEM analyses, at least 7 of the 10 associations were consistent with effects of the maternal genotype acting via the intrauterine environment, rather than via effects of shared alleles with the fetus. Variants, or correlated proxies, at many of the loci had been previously associated with adult traits, including fasting glucose (MTNR1B, GCK and TCF7L2) and sex hormone levels (CYP3A7), and one (EBF1) with gestational duration. The identified associations indicate that genetic effects on maternal glucose, cytochrome P450 activity and gestational duration, and potentially on maternal blood pressure and immune function, are relevant for fetal growth. Further characterization of these associations in mechanistic and causal analyses will enhance understanding of the potentially modifiable maternal determinants of fetal growth, with the goal of reducing the morbidity and mortality associated with low and high birth weights.

The moderating role of SES on genetic differences in educational achievement in the Netherlands

Parental socioeconomic status (SES) is a strong predictor of children's educational achievement (EA), with an increasing effect throughout development. Inequality in educational outcomes between children from different SES backgrounds exists in all Western countries. It has been proposed that a cause of this inequality lies in the interplay between genetic effects and SES on EA, which might depend on society and the equality of the education system. This study adopted two approaches, a classical twin design and polygenic score (PGS) approach, to address the effect of parental SES on EA in a large sample of 12-year-old Dutch twin pairs (2479 MZ and 4450 DZ twin pairs with PGSs for educational attainment available in 2335 children) from the Netherlands Twin Register (NTR). The findings of this study indicated that average EA increased with increasing parental SES. The difference in EA between boys and girls became smaller in the higher SES groups. The classical twin design analyses based on genetic covariance structure modeling pointed to lower genetic, environmental, and thus phenotypic variation in EA at higher SES. Independent from a child's PGS, parental SES predicted EA. However, the strength of the association between PGS and EA did not depend on parental SES. In a within-family design, the twin with a higher PGS scored higher on EA than the co-twin, demonstrating that the effect of the PGS on EA was at least partly independent from parental SES. To conclude, EA depended on SES both directly and indirectly, and SES moderated the additive genetic and environmental components of EA. Adding information from PGS, in addition to parental SES, improved the prediction of children's EA.

Publisher Correction: Genetic analysis of over 1 million people identifies 535 new loci associated with blood pressure traits

In the version of this article originally published, the name of author Martin H. de Borst was coded incorrectly in the XML. The error has now been corrected in the HTML version of the paper.

Genetic analysis of over 1 million people identifies 535 new loci associated with blood pressure traits

High blood pressure is a highly heritable and modifiable risk factor for cardiovascular disease. We report the largest genetic association study of blood pressure traits (systolic, diastolic and pulse pressure) to date in over 1 million people of European ancestry. We identify 535 novel blood pressure loci that not only offer new biological insights into blood pressure regulation but also highlight shared genetic architecture between blood pressure and lifestyle exposures. Our findings identify new biological pathways for blood pressure regulation with potential for improved cardiovascular disease prevention in the future.

Genetic vulnerability to schizophrenia is associated with cannabis use patterns during adolescence

BACKGROUND

Previously reported comorbidity between schizophrenia and substance use may be explained by shared underlying risk factors, such as genetic background. The aim of the present longitudinal study was to investigate how a genetic predisposition to schizophrenia was associated with patterns of substance use (cannabis use, smoking, alcohol use) during adolescence (comparing ages 13-16 with 16-20 years).

METHOD

Using piecewise latent growth curve modelling in a longitudinal adolescent cohort (RADAR-Y study, N = 372), we analyzed the association of polygenic risk scores for schizophrenia (PRS; p-value thresholds (p_t) < 5e-8 to p_t < 0.5) with increase in substance use over the years, including stratified analyses for gender. Significance thresholds were set to adjust for multiple testing using Bonferroni at p ≤ 0.001.

RESULTS

High schizophrenia vulnerability was associated with a stronger increase in cannabis use at age 16-20 (PRS thresholds p_t < 5e-5 and p_t < 5e-4; p_t < 5e-6 was marginally significant), whereas more lenient PRS thresholds (PRS thresholds p_t < 5e-3 to p_t < 0.5) showed the reverse association. For smoking and alcohol, no clear relations were found.

CONCLUSIONS

In conclusion, our findings support a relation between genetic risk to schizophrenia and prospective cannabis use patterns during adolescence. In contrast, no relation between alcohol and smoking was established.

Genome-wide association analysis links multiple psychiatric liability genes to oscillatory brain activity

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 GABA_A 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.

Analysis of shared heritability in common disorders of the brain

Disorders of the brain can exhibit considerable epidemiological comorbidity and often share symptoms, provoking debate about their etiologic overlap. We quantified the genetic sharing of 25 brain disorders from genome-wide association studies of 265,218 patients and 784,643 control participants and assessed their relationship to 17 phenotypes from 1,191,588 individuals. Psychiatric disorders share common variant risk, whereas neurological disorders appear more distinct from one another and from the psychiatric disorders. We also identified significant sharing between disorders and a number of brain phenotypes, including cognitive measures. Further, we conducted simulations to explore how statistical power, diagnostic misclassification, and phenotypic heterogeneity affect genetic correlations. These results highlight the importance of common genetic variation as a risk factor for brain disorders and the value of heritability-based methods in understanding their etiology.

Genome-wide meta-analysis of cognitive empathy: heritability, and correlates with sex, neuropsychiatric conditions and cognition

We conducted a genome-wide meta-analysis of cognitive empathy using the 'Reading the Mind in the Eyes' Test (Eyes Test) in 88,056 research volunteers of European Ancestry (44,574 females and 43,482 males) from 23andMe Inc., and an additional 1497 research volunteers of European Ancestry (891 females and 606 males) from the Brisbane Longitudinal Twin Study. We confirmed a female advantage on the Eyes Test (Cohen's d=0.21, P<2.2 × 10^-16), and identified a locus in 3p26.1 that is associated with scores on the Eyes Test in females (rs7641347, P_meta=1.58 × 10^-8). Common single nucleotide polymorphisms explained 5.8% (95% CI: 4.5%-7.2%; P=1.00 × 10^-17) of the total trait variance in both sexes, and we identified a twin heritability of 28% (95% CI: 13%-42%). Finally, we identified significant genetic correlation between the Eyes Test and anorexia nervosa, openness (NEO-Five Factor Inventory), and different measures of educational attainment and cognitive aptitude.

Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression

Major depressive disorder (MDD) is a common illness accompanied by considerable morbidity, mortality, costs, and heightened risk of suicide. We conducted a genome-wide association meta-analysis based in 135,458 cases and 344,901 controls and identified 44 independent and significant loci. The genetic findings were associated with clinical features of major depression and implicated brain regions exhibiting anatomical differences in cases. Targets of antidepressant medications and genes involved in gene splicing were enriched for smaller association signal. We found important relationships of genetic risk for major depression with educational attainment, body mass, and schizophrenia: lower educational attainment and higher body mass were putatively causal, whereas major depression and schizophrenia reflected a partly shared biological etiology. All humans carry lesser or greater numbers of genetic risk factors for major depression. These findings help refine the basis of major depression and imply that a continuous measure of risk underlies the clinical phenotype.

2SNP heritability and effects of genetic variants for neutrophil-to-lymphocyte and platelet-to-lymphocyte ratio

Neutrophil-to-lymphocyte ratio (NLR) and platelet-to-lymphocyte ratio (PLR) are important biomarkers for disease development and progression. To gain insight into the genetic causes of variance in NLR and PLR in the general population, we conducted genome-wide association (GWA) analyses and estimated SNP heritability in a sample of 5901 related healthy Dutch individuals. GWA analyses identified a new genome-wide significant locus on the HBS1L-MYB intergenic region for PLR, which replicated in a sample of 2538 British twins. For platelet count, we replicated three known genome-wide significant loci in our cohort (at CCDC71L-PIK3CG, BAK1 and ARHGEF3). For neutrophil count, we replicated the PSMD3 locus. For the identified top SNPs, we found significant cis and trans expression quantitative trait loci effects for several loci involved in hematological and immunological pathways. Linkage Disequilibrium score (LD) regression analyses for PLR and NLR confirmed that both traits are heritable, with a polygenetic SNP heritability for PLR of 14.1%, and for NLR of 2.4%. Genetic correlations were present between ratios and the constituent counts, with the genetic correlation (r=0.45) of PLR with platelet count reaching statistical significance. In conclusion, we established that two important biomarkers have a significant heritable SNP component, and identified the first genome-wide locus for PLR.

Genome-wide meta-analysis associates HLA-DQA1/DRB1 and LPA and lifestyle factors with human longevity

Genomic analysis of longevity offers the potential to illuminate the biology of human aging. Here, using genome-wide association meta-analysis of 606,059 parents’ survival, we discover two regions associated with longevity (HLA-DQA1/DRB1 and LPA). We also validate previous suggestions that APOE, CHRNA3/5, CDKN2A/B, SH2B3 and FOXO3A influence longevity. Next we show that giving up smoking, educational attainment, openness to new experience and high-density lipoprotein (HDL) cholesterol levels are most positively genetically correlated with lifespan while susceptibility to coronary artery disease (CAD), cigarettes smoked per day, lung cancer, insulin resistance and body fat are most negatively correlated. We suggest that the effect of education on lifespan is principally mediated through smoking while the effect of obesity appears to act via CAD. Using instrumental variables, we suggest that an increase of one body mass index unit reduces lifespan by 7 months while 1 year of education adds 11 months to expected lifespan.

Variability in human longevity is genetically influenced. Using genetic data of parental lifespan, the authors identify associations at HLA-DQA/DRB1 and LPA and find that genetic variants that increase educational attainment have a positive effect on lifespan whereas increasing BMI negatively affects lifespan.

Impact of common genetic determinants of Hemoglobin A1c on type 2 diabetes risk and diagnosis in ancestrally diverse populations: A transethnic genome-wide meta-analysis

BACKGROUND

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.

METHODS & FINDINGS

Using genome-wide association meta-analyses in up to 159,940 individuals from 82 cohorts of European, African, East Asian, and South Asian ancestry, we identified 60 common genetic variants associated with HbA1c. We classified variants as implicated in glycemic, erythrocytic, or unclassified biology and tested whether additive genetic scores of erythrocytic variants (GS-E) or glycemic variants (GS-G) were associated with higher T2D incidence in multiethnic longitudinal cohorts (N = 33,241). Nineteen glycemic and 22 erythrocytic variants were associated with HbA1c at genome-wide significance. GS-G was associated with higher T2D risk (incidence OR = 1.05, 95% CI 1.04-1.06, per HbA1c-raising allele, p = 3 × 10-29); whereas GS-E was not (OR = 1.00, 95% CI 0.99-1.01, p = 0.60). In Europeans and Asians, erythrocytic variants in aggregate had only modest effects on the diagnostic accuracy of HbA1c. Yet, in African Americans, the X-linked G6PD G202A variant (T-allele frequency 11%) was associated with an absolute decrease in HbA1c of 0.81%-units (95% CI 0.66-0.96) per allele in hemizygous men, and 0.68%-units (95% CI 0.38-0.97) in homozygous women. The G6PD variant may cause approximately 2% (N = 0.65 million, 95% CI 0.55-0.74) of African American adults with T2D to remain undiagnosed when screened with HbA1c. Limitations include the smaller sample sizes for non-European ancestries and the inability to classify approximately one-third of the variants. Further studies in large multiethnic cohorts with HbA1c, glycemic, and erythrocytic traits are required to better determine the biological action of the unclassified variants.

CONCLUSIONS

As G6PD deficiency can be clinically silent until illness strikes, we recommend investigation of the possible benefits of screening for the G6PD genotype along with using HbA1c to diagnose T2D in populations of African ancestry or groups where G6PD deficiency is common. Screening with direct glucose measurements, or genetically-informed HbA1c diagnostic thresholds in people with G6PD deficiency, may be required to avoid missed or delayed diagnoses.

Genome-wide Regional Heritability Mapping Identifies a Locus Within the TOX2 Gene Associated With Major Depressive Disorder

BACKGROUND

Major depressive disorder (MDD) is the second largest cause of global disease burden. It has an estimated heritability of 37%, but published genome-wide association studies have so far identified few risk loci. Haplotype-block-based regional heritability mapping (HRHM) estimates the localized genetic variance explained by common variants within haplotype blocks, integrating the effects of multiple variants, and may be more powerful for identifying MDD-associated genomic regions.

METHODS

We applied HRHM to Generation Scotland: The Scottish Family Health Study, a large family- and population-based Scottish cohort (N = 19,896). Single-single nucleotide polymorphism (SNP) and haplotype-based association tests were used to localize the association signal within the regions identified by HRHM. Functional prediction was used to investigate the effect of MDD-associated SNPs within the regions.

RESULTS

A haplotype block across a 24-kb region within the TOX2 gene reached genome-wide significance in HRHM. Single-SNP- and haplotype-based association tests demonstrated that five of nine genotyped SNPs and two haplotypes within this block were significantly associated with MDD. The expression of TOX2 and a brain-specific long noncoding RNA RP1-269M15.3 in frontal cortex and nucleus accumbens basal ganglia, respectively, were significantly regulated by MDD-associated SNPs within this region. Both the regional heritability and single-SNP associations within this block were replicated in the UK-Ireland group of the most recent release of the Psychiatric Genomics Consortium (PGC), the PGC2-MDD (Major Depression Dataset). The SNP association was also replicated in a depressive symptom sample that shares some individuals with the PGC2-MDD.

CONCLUSIONS

This study highlights the value of HRHM for MDD and provides an important target within TOX2 for further functional studies.

Novel Blood Pressure Locus and Gene Discovery Using Genome-Wide Association Study and Expression Data Sets From Blood and the Kidney

Elevated blood pressure is a major risk factor for cardiovascular disease and has a substantial genetic contribution. Genetic variation influencing blood pressure has the potential to identify new pharmacological targets for the treatment of hypertension. To discover additional novel blood pressure loci, we used 1000 Genomes Project-based imputation in 150 134 European ancestry individuals and sought significant evidence for independent replication in a further 228 245 individuals. We report 6 new signals of association in or near HSPB7, TNXB, LRP12, LOC283335, SEPT9, and AKT2, and provide new replication evidence for a further 2 signals in EBF2 and NFKBIA. Combining large whole-blood gene expression resources totaling 12 607 individuals, we investigated all novel and previously reported signals and identified 48 genes with evidence for involvement in blood pressure regulation that are significant in multiple resources. Three novel kidney-specific signals were also detected. These robustly implicated genes may provide new leads for therapeutic innovation.

Interaction between the FTO gene, body mass index and depression: meta-analysis of 13701 individuals

BackgroundDepression and obesity are highly prevalent, and major impacts on public health frequently co-occur. Recently, we reported that having depression moderates the effect of the FTO gene, suggesting its implication in the association between depression and obesity.AimsTo confirm these findings by investigating the FTO polymorphism rs9939609 in new cohorts, and subsequently in a meta-analysis.MethodThe sample consists of 6902 individuals with depression and 6799 controls from three replication cohorts and two original discovery cohorts. Linear regression models were performed to test for association between rs9939609 and body mass index (BMI), and for the interaction between rs9939609 and depression status for an effect on BMI. Fixed and random effects meta-analyses were performed using METASOFT.ResultsIn the replication cohorts, we observed a significant interaction between FTO, BMI and depression with fixed effects meta-analysis (β = 0.12, P = 2.7 × 10^-4) and with the Han/Eskin random effects method (P = 1.4 × 10^-7) but not with traditional random effects (β = 0.1, P = 0.35). When combined with the discovery cohorts, random effects meta-analysis also supports the interaction (β = 0.12, P = 0.027) being highly significant based on the Han/Eskin model (P = 6.9 × 10^-8). On average, carriers of the risk allele who have depression have a 2.2% higher BMI for each risk allele, over and above the main effect of FTOConclusionsThis meta-analysis provides additional support for a significant interaction between FTO, depression and BMI, indicating that depression increases the effect of FTO on BMI. The findings provide a useful starting point in understanding the biological mechanism involved in the association between obesity and depression.

Using a multivariate model to assess the interactive effects of demographics and lifestyle on the hematological profile

AIM

To assess the extent to which a multivariate approach to modeling interrelated hematological indices provides more informative results than the traditional approach of modeling each index separately.

MATERIALS & METHODS

The effects of demographics and lifestyle on ten hematological indices collected from a Dutch population-based sample (n = 3278) were studied, jointly using multivariate distance matrix regression and separately using linear regression.

RESULTS

The multivariate approach highlighted the main effects of all predictors and several interactions; the traditional approach highlighted only main effects.

CONCLUSION

The multivariate approach provides more power than traditional methods to detect effects on interrelated biomarkers, suggesting that its use in future research may help identify subgroups that benefit from different treatment or prevention measures.

Genetic loci associated with heart rate variability and their effects on cardiac disease risk

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<r_g<-0.55) and blood pressure (-0.35<r_g<-0.20). These findings provide clinically relevant biological insight into heritable variation in vagal heart rhythm regulation, with a key role for genetic variants (GNG11, RGS6) that influence G-protein heterotrimer action in GIRK-channel induced pacemaker membrane hyperpolarization.

Conditional eQTL analysis reveals allelic heterogeneity of gene expression

In recent years, multiple eQTL (expression quantitative trait loci) catalogs have become available that can help understand the functionality of complex trait-related single nucleotide polymorphisms (SNPs). In eQTL catalogs, gene expression is often strongly associated with multiple SNPs, which may reflect either one or multiple independent associations. Conditional eQTL analysis allows a distinction between dependent and independent eQTLs. We performed conditional eQTL analysis in 4,896 peripheral blood microarray gene expression samples. Our analysis showed that 35% of genes with a cis eQTL have at least two independent cis eQTLs; for several genes up to 13 independent cis eQTLs were identified. Also, 12% (671) of the independent cis eQTLs identified in conditional analyses were not significant in unconditional analyses. The number of GWAS catalog SNPs identified as eQTL in the conditional analyses increases with 24% as compared to unconditional analyses. We provide an online conditional cis eQTL mapping catalog for whole blood (https://eqtl.onderzoek.io/), which can be used to lookup eQTLs more accurately than in standard unconditional whole blood eQTL databases.

Short communication: Genetic association between schizophrenia and cannabis use

BACKGROUND AND AIM

Previous studies have shown a relationship between schizophrenia and cannabis use. As both traits are substantially heritable, a shared genetic liability could explain the association. We use two recently developed genomics methods to investigate the genetic overlap between schizophrenia and cannabis use.

METHODS

Firstly, polygenic risk scores for schizophrenia were created based on summary statistics from the largest schizophrenia genome-wide association (GWA) meta-analysis to date. We analysed the association between these schizophrenia polygenic scores and multiple cannabis use phenotypes (lifetime use, regular use, age at initiation, and quantity and frequency of use) in a sample of 6,931 individuals. Secondly, we applied LD-score regression to the GWA summary statistics of schizophrenia and lifetime cannabis use to calculate the genome-wide genetic correlation.

RESULTS

Polygenic risk scores for schizophrenia were significantly (α<0.05) associated with five of the eight cannabis use phenotypes, including lifetime use, regular use, and quantity of use, with risk scores explaining up to 0.5% of the variance. Associations were not significant for age at initiation of use and two measures of frequency of use analyzed in lifetime users only, potentially because of reduced power due to a smaller sample size. The LD-score regression revealed a significant genetic correlation of r_g=0.22 (SE=0.07, p=0.003) between schizophrenia and lifetime cannabis use.

CONCLUSIONS

Common genetic variants underlying schizophrenia and lifetime cannabis use are partly overlapping. Individuals with a stronger genetic predisposition to schizophrenia are more likely to initiate cannabis use, use cannabis more regularly, and consume more cannabis over their lifetime.

Novel genetic loci associated with hippocampal volume

The hippocampal formation is a brain structure integrally involved in episodic memory, spatial navigation, cognition and stress responsiveness. Structural abnormalities in hippocampal volume and shape are found in several common neuropsychiatric disorders. To identify the genetic underpinnings of hippocampal structure here we perform a genome-wide association study (GWAS) of 33,536 individuals and discover six independent loci significantly associated with hippocampal volume, four of them novel. Of the novel loci, three lie within genes (ASTN2, DPP4 and MAST4) and one is found 200 kb upstream of SHH. A hippocampal subfield analysis shows that a locus within the MSRB3 gene shows evidence of a localized effect along the dentate gyrus, subiculum, CA1 and fissure. Further, we show that genetic variants associated with decreased hippocampal volume are also associated with increased risk for Alzheimer's disease (rg=-0.155). Our findings suggest novel biological pathways through which human genetic variation influences hippocampal volume and risk for neuropsychiatric illness.

Novel genetic loci underlying human intracranial volume identified through genome-wide association

Intracranial volume reflects the maximally attained brain size during development, and remains stable with loss of tissue in late life. It is highly heritable, but the underlying genes remain largely undetermined. In a genome-wide association study of 32,438 adults, we discovered five previously unknown loci for intracranial volume and confirmed two known signals. Four of the loci were also associated with adult human stature, but these remained associated with intracranial volume after adjusting for height. We found a high genetic correlation with child head circumference (ρgenetic = 0.748), which indicates a similar genetic background and allowed us to identify four additional loci through meta-analysis (Ncombined = 37,345). Variants for intracranial volume were also related to childhood and adult cognitive function, and Parkinson's disease, and were enriched near genes involved in growth pathways, including PI3K-AKT signaling. These findings identify the biological underpinnings of intracranial volume and their link to physiological and pathological traits.

A principal component meta-analysis on multiple anthropometric traits identifies novel loci for body shape

Large consortia have revealed hundreds of genetic loci associated with anthropometric traits, one trait at a time. We examined whether genetic variants affect body shape as a composite phenotype that is represented by a combination of anthropometric traits. We developed an approach that calculates averaged PCs (AvPCs) representing body shape derived from six anthropometric traits (body mass index, height, weight, waist and hip circumference, waist-to-hip ratio). The first four AvPCs explain >99% of the variability, are heritable, and associate with cardiometabolic outcomes. We performed genome-wide association analyses for each body shape composite phenotype across 65 studies and meta-analysed summary statistics. We identify six novel loci: LEMD2 and CD47 for AvPC1, RPS6KA5/C14orf159 and GANAB for AvPC3, and ARL15 and ANP32 for AvPC4. Our findings highlight the value of using multiple traits to define complex phenotypes for discovery, which are not captured by single-trait analyses, and may shed light onto new pathways.

Heritability of Behavioral Problems in 7-Year Olds Based on Shared and Unique Aspects of Parental Views

In studies of child psychopathology, phenotypes of interest are often obtained by parental ratings. When behavioral ratings are obtained in the context of a twin study, this allows for the decomposition of the phenotypic variance, into a genetic and a non-genetic part. If a phenotype is assessed by a single rater, heritability is based on the child’s behavior as expressed in the presence of that particular rater, whereas heritability based on assessments by multiple raters allows for the estimation of the heritability of the phenotype based on rater agreement, as well as the heritability of the rater specific view of the behavior. The aim of this twin study was to quantify the rater common and rater specific contributions to the variation in children’s behavioral problems. We estimated the heritability of maternal and paternal ratings of the Child Behavior Checklist (CBCL) 6–18 empirical emotional and behavioral problem scales in a large sample of 12,310 7-year old Dutch twin pairs. Between 30 and 59% of variation in the part of the phenotype parents agree upon was explained by genetic effects. Common environmental effects that make children in the same family similar explained less variance, ranging between 0 and 32%. For unique views of their children’s behavioral problems, heritability ranged between 0 and 20% for maternal and between 0 and 22% for paternal views. Between 7 and 24% of the variance was accounted for by common environmental factors specific to mother and father’s views. The proportion of rater shared and rater specific heritability can be translated into genetic correlations between parental views and inform the design and interpretation of results of molecular genetic studies. Genetic correlations were nearly or above 0.7 for all CBCL based psychopathology scales. Such large genetic correlations suggest two practical guidelines for genome-wide association studies (GWAS): when studies have collected data from either fathers or mothers, the shared genetic aetiology in parental ratings indicates that is possible to analyze paternal and maternal assessments in a single GWAS or meta-analysis. Secondly, if a study has collected information from both parents, a gain in statistical power may be realized in GWAS by the simultaneous analysis of the data.

Genome-wide associations for birth weight and correlations with adult disease

Birth weight (BW) has been shown to be influenced by both fetal and maternal factors and in observational studies is reproducibly associated with future risk of adult metabolic diseases including type 2 diabetes (T2D) and cardiovascular disease. These life-course associations have often been attributed to the impact of an adverse early life environment. Here, we performed a multi-ancestry genome-wide association study (GWAS) meta-analysis of BW in 153,781 individuals, identifying 60 loci where fetal genotype was associated with BW (P < 5 × 10-8). Overall, approximately 15% of variance in BW was captured by assays of fetal genetic variation. Using genetic association alone, we found strong inverse genetic correlations between BW and systolic blood pressure (Rg = -0.22, P = 5.5 × 10-13), T2D (Rg = -0.27, P = 1.1 × 10-6) and coronary artery disease (Rg = -0.30, P = 6.5 × 10-9). In addition, using large -cohort datasets, we demonstrated that genetic factors were the major contributor to the negative covariance between BW and future cardiometabolic risk. Pathway analyses indicated that the protein products of genes within BW-associated regions were enriched for diverse processes including insulin signalling, glucose homeostasis, glycogen biosynthesis and chromatin remodelling. There was also enrichment of associations with BW in known imprinted regions (P = 1.9 × 10-4). We demonstrate that life-course associations between early growth phenotypes and adult cardiometabolic disease are in part the result of shared genetic effects and identify some of the pathways through which these causal genetic effects are mediated.

Causes of variation in the neutrophil-lymphocyte and platelet-lymphocyte ratios: a twin-family study

AIM

Neutrophil-lymphocyte ratio (NLR) and platelet-lymphocyte ratio (PLR) are biomarkers for disease development, for whom little is known about causes of variation in the general population.

MATERIALS & METHODS

We estimated the heritability of PLR and NLR and examined their association with gender, demographic, lifestyle and environmental factors in a Dutch nonpatient twin family population (n = 8108).

RESULTS

Heritability was estimated at 64% for PLR and 36% for NLR. Men had on average higher NLR, but lower PLR levels than women. PLR and NLR increased significantly with age, decreased in colder months and showed small but significant sex- and age-specific associations with body composition and smoking.

CONCLUSION

NLR and PLR levels are heritable and influenced by age, sex and environmental factors, such as seasonal conditions and lifestyle.

A Genome-Wide Association Meta-Analysis of Attention-Deficit/Hyperactivity Disorder Symptoms in Population-Based Pediatric Cohorts

OBJECTIVE

The aims of this study were to elucidate the influence of common genetic variants on childhood attention-deficit/hyperactivity disorder (ADHD) symptoms, to identify genetic variants that explain its high heritability, and to investigate the genetic overlap of ADHD symptom scores with ADHD diagnosis.

METHOD

Within the EArly Genetics and Lifecourse Epidemiology (EAGLE) consortium, genome-wide single nucleotide polymorphisms (SNPs) and ADHD symptom scores were available for 17,666 children (<13 years of age) from nine population-based cohorts. SNP-based heritability was estimated in data from the three largest cohorts. Meta-analysis based on genome-wide association (GWA) analyses with SNPs was followed by gene-based association tests, and the overlap in results with a meta-analysis in the Psychiatric Genomics Consortium (PGC) case-control ADHD study was investigated.

RESULTS

SNP-based heritability ranged from 5% to 34%, indicating that variation in common genetic variants influences ADHD symptom scores. The meta-analysis did not detect genome-wide significant SNPs, but three genes, lying close to each other with SNPs in high linkage disequilibrium (LD), showed a gene-wide significant association (p values between 1.46 × 10(-6) and 2.66 × 10(-6)). One gene, WASL, is involved in neuronal development. Both SNP- and gene-based analyses indicated overlap with the PGC meta-analysis results with the genetic correlation estimated at 0.96.

CONCLUSION

The SNP-based heritability for ADHD symptom scores indicates a polygenic architecture, and genes involved in neurite outgrowth are possibly involved. Continuous and dichotomous measures of ADHD appear to assess a genetically common phenotype. A next step is to combine data from population-based and case-control cohorts in genetic association studies to increase sample size and to improve statistical power for identifying genetic variants.

The genetics of blood pressure regulation and its target organs from association studies in 342,415 individuals

To dissect the genetic architecture of blood pressure and assess effects on target organ damage, we analyzed 128,272 SNPs from targeted and genome-wide arrays in 201,529 individuals of European ancestry, and genotypes from an additional 140,886 individuals were used for validation. We identified 66 blood pressure-associated loci, of which 17 were new; 15 harbored multiple distinct association signals. The 66 index SNPs were enriched for cis-regulatory elements, particularly in vascular endothelial cells, consistent with a primary role in blood pressure control through modulation of vascular tone across multiple tissues. The 66 index SNPs combined in a risk score showed comparable effects in 64,421 individuals of non-European descent. The 66-SNP blood pressure risk score was significantly associated with target organ damage in multiple tissues but with minor effects in the kidney. Our findings expand current knowledge of blood pressure-related pathways and highlight tissues beyond the classical renal system in blood pressure regulation.

Detection of gene-environment interaction in pedigree data using genome-wide genotypes

Heritability may be estimated using phenotypic data collected in relatives or in distantly related individuals using genome-wide single nucleotide polymorphism (SNP) data. We combined these approaches by re-parameterizing the model proposed by Zaitlen et al and extended this model to include moderation of (total and SNP-based) genetic and environmental variance components by a measured moderator. By means of data simulation, we demonstrated that the type 1 error rates of the proposed test are correct and parameter estimates are accurate. As an application, we considered the moderation by age or year of birth of variance components associated with body mass index (BMI), height, attention problems (AP), and symptoms of anxiety and depression. The genetic variance of BMI was found to increase with age, but the environmental variance displayed a greater increase with age, resulting in a proportional decrease of the heritability of BMI. Environmental variance of height increased with year of birth. The environmental variance of AP increased with age. These results illustrate the assessment of moderation of environmental and genetic effects, when estimating heritability from combined SNP and family data. The assessment of moderation of genetic and environmental variance will enhance our understanding of the genetic architecture of complex traits.

Correction: The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study

[This corrects the article DOI: 10.1371/journal.pgen.1005378.].

The 3p21.1-p21.3 Hereditary Vascular Retinopathy Locus Increases the Risk for Raynaud's Phenomenon and Migraine

Previously, we described a large Dutch family with hereditary vascular retinopathy (HVR), Raynaud's phenomenon and migraine. A locus for HVR was mapped on chromosome 3p21.1-p21.3, but the gene has not yet been identified. The fact that all three disorders share a vascular aetiology prompted us to study whether the HVR haplotype also contributed to Raynaud's phenomenon and migraine in this family. Whereas the parent-child transmission disequilibrium test (TDT) did not reach significance, the sibling TDT revealed that the HVR haplotype harbours a susceptibility factor for Raynaud's phenomenon and migraine. Identification of the HVR gene will improve the understanding of the pathophysiology of HVR, Raynaud's phenomenon and migraine.

Meta-analysis of 375,000 individuals identifies 38 susceptibility loci for migraine

Migraine is a debilitating neurological disorder affecting around one in seven people worldwide, but its molecular mechanisms remain poorly understood. There is some debate about whether migraine is a disease of vascular dysfunction or a result of neuronal dysfunction with secondary vascular changes. Genome-wide association (GWA) studies have thus far identified 13 independent loci associated with migraine. To identify new susceptibility loci, we carried out a genetic study of migraine on 59,674 affected subjects and 316,078 controls from 22 GWA studies. We identified 44 independent single-nucleotide polymorphisms (SNPs) significantly associated with migraine risk (P < 5 × 10(-8)) that mapped to 38 distinct genomic loci, including 28 loci not previously reported and a locus that to our knowledge is the first to be identified on chromosome X. In subsequent computational analyses, the identified loci showed enrichment for genes expressed in vascular and smooth muscle tissues, consistent with a predominant theory of migraine that highlights vascular etiologies.

International Genome-Wide Association Study Consortium Identifies Novel Loci Associated With Blood Pressure in Children and Adolescents

BACKGROUND

Our aim was to identify genetic variants associated with blood pressure (BP) in childhood and adolescence.

METHODS AND RESULTS

Genome-wide association study data from participating European ancestry cohorts of the Early Genetics and Lifecourse Epidemiology (EAGLE) Consortium was meta-analyzed across 3 epochs; prepuberty (4-7 years), puberty (8-12 years), and postpuberty (13-20 years). Two novel loci were identified as having genome-wide associations with systolic BP across specific age epochs: rs1563894 (ITGA11, located in active H3K27Ac mark and transcription factor chromatin immunoprecipitation and 5'-C-phosphate-G-3' methylation site) during prepuberty (P=2.86×10(-8)) and rs872256 during puberty (P=8.67×10(-9)). Several single-nucleotide polymorphism clusters were also associated with childhood BP at P<5×10(-3). Using a P value threshold of <5×10(-3), we found some overlap in variants across the different age epochs within our study and between several single-nucleotide polymorphisms in any of the 3 epochs and adult BP-related single-nucleotide polymorphisms.

CONCLUSIONS

Our results suggest that genetic determinants of BP act from childhood, develop over the lifecourse, and show some evidence of age-specific effects.

Genetic variants associated with subjective well-being, depressive symptoms, and neuroticism identified through genome-wide analyses

Very few genetic variants have been associated with depression and neuroticism, likely because of limitations on sample size in previous studies. Subjective well-being, a phenotype that is genetically correlated with both of these traits, has not yet been studied with genome-wide data. We conducted genome-wide association studies of three phenotypes: subjective well-being (n = 298,420), depressive symptoms (n = 161,460), and neuroticism (n = 170,911). We identify 3 variants associated with subjective well-being, 2 variants associated with depressive symptoms, and 11 variants associated with neuroticism, including 2 inversion polymorphisms. The two loci associated with depressive symptoms replicate in an independent depression sample. Joint analyses that exploit the high genetic correlations between the phenotypes (|ρ^| ≈ 0.8) strengthen the overall credibility of the findings and allow us to identify additional variants. Across our phenotypes, loci regulating expression in central nervous system and adrenal or pancreas tissues are strongly enriched for association.

Genome-wide association study identifies 74 loci associated with educational attainment

Educational attainment is strongly influenced by social and other environmental factors, but genetic factors are estimated to account for at least 20% of the variation across individuals. Here we report the results of a genome-wide association study (GWAS) for educational attainment that extends our earlier discovery sample of 101,069 individuals to 293,723 individuals, and a replication study in an independent sample of 111,349 individuals from the UK Biobank. We identify 74 genome-wide significant loci associated with the number of years of schooling completed. Single-nucleotide polymorphisms associated with educational attainment are disproportionately found in genomic regions regulating gene expression in the fetal brain. Candidate genes are preferentially expressed in neural tissue, especially during the prenatal period, and enriched for biological pathways involved in neural development. Our findings demonstrate that, even for a behavioural phenotype that is mostly environmentally determined, a well-powered GWAS identifies replicable associated genetic variants that suggest biologically relevant pathways. Because educational attainment is measured in large numbers of individuals, it will continue to be useful as a proxy phenotype in efforts to characterize the genetic influences of related phenotypes, including cognition and neuropsychiatric diseases.