Polygenic dissection of major depression clinical heterogeneity

The molecular mechanisms underlying major depressive disorder (MDD) are largely unknown. Limited success of previous genetics studies may be attributable to heterogeneity of MDD, aggregating biologically different subtypes. We examined the polygenic features of MDD and two common clinical subtypes (typical and atypical) defined by symptom profiles in a large sample of adults with established diagnoses. Data were from 1530 patients of the Netherlands Study of Depression and Anxiety (NESDA) and 1700 controls mainly from the Netherlands Twin Register (NTR). Diagnoses of MDD and its subtypes were based on DSM-IV symptoms. Genetic overlap of MDD and subtypes with psychiatric (MDD, bipolar disorder, schizophrenia) and metabolic (body mass index (BMI), C-reactive protein, triglycerides) traits was evaluated via genomic profile risk scores (GPRS) generated from meta-analysis results of large international consortia. Single nucleotide polymorphism (SNP)-heritability of MDD and subtypes was also estimated. MDD was associated with psychiatric GPRS, while no association was found for GPRS of metabolic traits. MDD subtypes had differential polygenic signatures: typical was strongly associated with schizophrenia GPRS (odds ratio (OR)=1.54, P=7.8e-9), while atypical was additionally associated with BMI (OR=1.29, P=2.7e-4) and triglycerides (OR=1.21, P=0.006) GPRS. Similar results were found when only the highly discriminatory symptoms of appetite/weight were used to define subtypes. SNP-heritability was 32% for MDD, 38% and 43% for subtypes with, respectively, decreased (typical) and increased (atypical) appetite/weight. In conclusion, MDD subtypes are characterized by partially distinct polygenic liabilities and may represent more homogeneous phenotypes. Disentangling MDD heterogeneity may help the psychiatric field moving forward in the search for molecular roots of depression.

Genome-wide association study of lifetime cannabis use based on a large meta-analytic sample of 32 330 subjects from the International Cannabis Consortium

Cannabis is the most widely produced and consumed illicit psychoactive substance worldwide. Occasional cannabis use can progress to frequent use, abuse and dependence with all known adverse physical, psychological and social consequences. Individual differences in cannabis initiation are heritable (40-48%). The International Cannabis Consortium was established with the aim to identify genetic risk variants of cannabis use. We conducted a meta-analysis of genome-wide association data of 13 cohorts (N=32 330) and four replication samples (N=5627). In addition, we performed a gene-based test of association, estimated single-nucleotide polymorphism (SNP)-based heritability and explored the genetic correlation between lifetime cannabis use and cigarette use using LD score regression. No individual SNPs reached genome-wide significance. Nonetheless, gene-based tests identified four genes significantly associated with lifetime cannabis use: NCAM1, CADM2, SCOC and KCNT2. Previous studies reported associations of NCAM1 with cigarette smoking and other substance use, and those of CADM2 with body mass index, processing speed and autism disorders, which are phenotypes previously reported to be associated with cannabis use. Furthermore, we showed that, combined across the genome, all common SNPs explained 13-20% (P<0.001) of the liability of lifetime cannabis use. Finally, there was a strong genetic correlation (rg=0.83; P=1.85 × 10(-8)) between lifetime cannabis use and lifetime cigarette smoking implying that the SNP effect sizes of the two traits are highly correlated. This is the largest meta-analysis of cannabis GWA studies to date, revealing important new insights into the genetic pathways of lifetime cannabis use. Future functional studies should explore the impact of the identified genes on the biological mechanisms of cannabis use.

Genome-wide study for circulating metabolites identifies 62 loci and reveals novel systemic effects of LPA

Genome-wide association studies have identified numerous loci linked with complex diseases, for which the molecular mechanisms remain largely unclear. Comprehensive molecular profiling of circulating metabolites captures highly heritable traits, which can help to uncover metabolic pathophysiology underlying established disease variants. We conduct an extended genome-wide association study of genetic influences on 123 circulating metabolic traits quantified by nuclear magnetic resonance metabolomics from up to 24,925 individuals and identify eight novel loci for amino acids, pyruvate and fatty acids. The LPA locus link with cardiovascular risk exemplifies how detailed metabolic profiling may inform underlying aetiology via extensive associations with very-low-density lipoprotein and triglyceride metabolism. Genetic fine mapping and Mendelian randomization uncover wide-spread causal effects of lipoprotein(a) on overall lipoprotein metabolism and we assess potential pleiotropic consequences of genetically elevated lipoprotein(a) on diverse morbidities via electronic health-care records. Our findings strengthen the argument for safe LPA-targeted intervention to reduce cardiovascular risk.

Circulating metabolites reflect human health and disease. Here, Kettunen et al. perform a genome-wide association study on 123 circulating metabolic traits and identify novel genetic loci influencing systemic metabolism. They also link new molecular pathways with a known cardiovascular risk factor Lp(a).

Genetic Evidence for Causal Relationships Between Maternal Obesity-Related Traits and Birth Weight

IMPORTANCE

Neonates born to overweight or obese women are larger and at higher risk of birth complications. Many maternal obesity-related traits are observationally associated with birth weight, but the causal nature of these associations is uncertain.

OBJECTIVE

To test for genetic evidence of causal associations of maternal body mass index (BMI) and related traits with birth weight.

DESIGN, SETTING, AND PARTICIPANTS

Mendelian randomization to test whether maternal BMI and obesity-related traits are potentially causally related to offspring birth weight. Data from 30,487 women in 18 studies were analyzed. Participants were of European ancestry from population- or community-based studies in Europe, North America, or Australia and were part of the Early Growth Genetics Consortium. Live, term, singleton offspring born between 1929 and 2013 were included.

EXPOSURES

Genetic scores for BMI, fasting glucose level, type 2 diabetes, systolic blood pressure (SBP), triglyceride level, high-density lipoprotein cholesterol (HDL-C) level, vitamin D status, and adiponectin level.

MAIN OUTCOME AND MEASURE

Offspring birth weight from 18 studies.

RESULTS

Among the 30,487 newborns the mean birth weight in the various cohorts ranged from 3325 g to 3679 g. The maternal genetic score for BMI was associated with a 2-g (95% CI, 0 to 3 g) higher offspring birth weight per maternal BMI-raising allele (P = .008). The maternal genetic scores for fasting glucose and SBP were also associated with birth weight with effect sizes of 8 g (95% CI, 6 to 10 g) per glucose-raising allele (P = 7 × 10(-14)) and -4 g (95% CI, -6 to -2 g) per SBP-raising allele (P = 1×10(-5)), respectively. A 1-SD ( ≈ 4 points) genetically higher maternal BMI was associated with a 55-g higher offspring birth weight (95% CI, 17 to 93 g). A 1-SD ( ≈ 7.2 mg/dL) genetically higher maternal fasting glucose concentration was associated with 114-g higher offspring birth weight (95% CI, 80 to 147 g). However, a 1-SD ( ≈ 10 mm Hg) genetically higher maternal SBP was associated with a 208-g lower offspring birth weight (95% CI, -394 to -21 g). For BMI and fasting glucose, genetic associations were consistent with the observational associations, but for systolic blood pressure, the genetic and observational associations were in opposite directions.

CONCLUSIONS AND RELEVANCE

In this mendelian randomization study, genetically elevated maternal BMI and blood glucose levels were potentially causally associated with higher offspring birth weight, whereas genetically elevated maternal SBP was potentially causally related to lower birth weight. If replicated, these findings may have implications for counseling and managing pregnancies to avoid adverse weight-related birth outcomes.

Meta-analysis of Genome-Wide Association Studies for Extraversion: Findings from the Genetics of Personality Consortium

Extraversion is a relatively stable and heritable personality trait associated with numerous psychosocial, lifestyle and health outcomes. Despite its substantial heritability, no genetic variants have been detected in previous genome-wide association (GWA) studies, which may be due to relatively small sample sizes of those studies. Here, we report on a large meta-analysis of GWA studies for extraversion in 63,030 subjects in 29 cohorts. Extraversion item data from multiple personality inventories were harmonized across inventories and cohorts. No genome-wide significant associations were found at the single nucleotide polymorphism (SNP) level but there was one significant hit at the gene level for a long non-coding RNA site (LOC101928162). Genome-wide complex trait analysis in two large cohorts showed that the additive variance explained by common SNPs was not significantly different from zero, but polygenic risk scores, weighted using linkage information, significantly predicted extraversion scores in an independent cohort. These results show that extraversion is a highly polygenic personality trait, with an architecture possibly different from other complex human traits, including other personality traits. Future studies are required to further determine which genetic variants, by what modes of gene action, constitute the heritable nature of extraversion.

Gene expression in major depressive disorder

The search for genetic variants underlying major depressive disorder (MDD) has not yet provided firm leads to its underlying molecular biology. A complementary approach is to study gene expression in relation to MDD. We measured gene expression in peripheral blood from 1848 subjects from The Netherlands Study of Depression and Anxiety. Subjects were divided into current MDD (N=882), remitted MDD (N=635) and control (N=331) groups. MDD status and gene expression were measured again 2 years later in 414 subjects. The strongest gene expression differences were between the current MDD and control groups (129 genes at false-discovery rate, FDR<0.1). Gene expression differences across MDD status were largely unrelated to antidepressant use, inflammatory status and blood cell counts. Genes associated with MDD were enriched for interleukin-6 (IL-6)-signaling and natural killer (NK) cell pathways. We identified 13 gene expression clusters with specific clusters enriched for genes involved in NK cell activation (downregulated in current MDD, FDR=5.8 × 10(-5)) and IL-6 pathways (upregulated in current MDD, FDR=3.2 × 10(-3)). Longitudinal analyses largely confirmed results observed in the cross-sectional data. Comparisons of gene expression results to the Psychiatric Genomics Consortium (PGC) MDD genome-wide association study results revealed overlap with DVL3. In conclusion, multiple gene expression associations with MDD were identified and suggest a measurable impact of current MDD state on gene expression. Identified genes and gene clusters are enriched with immune pathways previously associated with the etiology of MDD, in line with the immune suppression and immune activation hypothesis of MDD.

The genetics of alcohol dependence: Twin and SNP-based heritability, and genome-wide association study based on AUDIT scores

Alcohol dependence (AD) is among the most common and costly public health problems contributing to morbidity and mortality throughout the world. In this study, we investigate the genetic basis of AD in a Dutch population using data from the Netherlands Twin Register (NTR) and the Netherlands Study of Depression and Anxiety (NESDA). The presence of AD was ascertained via the Alcohol Use Disorders Identification Test (AUDIT) applying cut-offs with good specificity and sensitivity in identifying those at risk for AD. Twin-based heritability of AD-AUDIT was estimated using structural equation modeling of data in 7,694 MZ and DZ twin pairs. Variance in AD-AUDIT explained by all SNPs was estimated with genome-wide complex trait analysis (GCTA). A genome-wide association study (GWAS) was performed in 7,842 subjects. GWAS SNP effect concordance analysis was performed between our GWAS and a recent AD GWAS using DSM-IV diagnosis. The twin-based heritability of AD-AUDIT was estimated at 60% (55-69%). GCTA showed that common SNPs jointly capture 33% (SE = 0.12, P = 0.002) of this heritability. In the GWAS, the top hits were positioned within four regions (4q31.1, 2p16.1, 6q25.1, 7p14.1) with the strongest association detected for rs55768019 (P = 7.58 × 10(-7) ). This first GWAS of AD using the AUDIT measure found results consistent with previous genetic studies using DSM diagnosis: concordance in heritability estimates and direction of SNPs effect and overlap with top hits from previous GWAS. Thus, the use of appropriate questionnaires may represent cost-effective strategies to phenotype samples in large-scale biobanks or other population-based datasets.

Fine mapping the CETP region reveals a common intronic insertion associated to HDL-C

BACKGROUND

Individuals with exceptional longevity and their offspring have significantly larger high-density lipoprotein concentrations (HDL-C) particle sizes due to the increased homozygosity for the I405V variant in the cholesteryl ester transfer protein (CETP) gene. In this study, we investigate the association of CETP and HDL-C further to identify novel, independent CETP variants associated with HDL-C in humans.

METHODS

We performed a meta-analysis of HDL-C within the CETP region using 59,432 individuals imputed with 1000 Genomes data. We performed replication in an independent sample of 47,866 individuals and validation was done by Sanger sequencing.

RESULTS

The meta-analysis of HDL-C within the CETP region identified five independent variants, including an exonic variant and a common intronic insertion. We replicated these 5 variants significantly in an independent sample of 47,866 individuals. Sanger sequencing of the insertion within a single family confirmed segregation of this variant. The strongest reported association between HDL-C and CETP variants, was rs3764261; however, after conditioning on the five novel variants we identified the support for rs3764261 was highly reduced (β_unadjusted=3.179 mg/dl (P value=5.25×10^-509), β_adjusted=0.859 mg/dl (P value=9.51×10^-25)), and this finding suggests that these five novel variants may partly explain the association of CETP with HDL-C. Indeed, three of the five novel variants (rs34065661, rs5817082, rs7499892) are independent of rs3764261.

CONCLUSIONS

The causal variants in CETP that account for the association with HDL-C remain unknown. We used studies imputed to the 1000 Genomes reference panel for fine mapping of the CETP region. We identified and validated five variants within this region that may partly account for the association of the known variant (rs3764261), as well as other sources of genetic contribution to HDL-C.

Large-scale genomic analyses link reproductive aging to hypothalamic signaling, breast cancer susceptibility and BRCA1-mediated DNA repair

Menopause timing has a substantial impact on infertility and risk of disease, including breast cancer, but the underlying mechanisms are poorly understood. We report a dual strategy in ∼70,000 women to identify common and low-frequency protein-coding variation associated with age at natural menopause (ANM). We identified 44 regions with common variants, including two regions harboring additional rare missense alleles of large effect. We found enrichment of signals in or near genes involved in delayed puberty, highlighting the first molecular links between the onset and end of reproductive lifespan. Pathway analyses identified major association with DNA damage response (DDR) genes, including the first common coding variant in BRCA1 associated with any complex trait. Mendelian randomization analyses supported a causal effect of later ANM on breast cancer risk (∼6% increase in risk per year; P = 3 × 10(-14)), likely mediated by prolonged sex hormone exposure rather than DDR mechanisms.

Novel loci associated with usual sleep duration: the CHARGE Consortium Genome-Wide Association Study

Usual sleep duration is a heritable trait correlated with psychiatric morbidity, cardiometabolic disease and mortality, although little is known about the genetic variants influencing this trait. A genome-wide association study (GWAS) of usual sleep duration was conducted using 18 population-based cohorts totaling 47 180 individuals of European ancestry. Genome-wide significant association was identified at two loci. The strongest is located on chromosome 2, in an intergenic region 35- to 80-kb upstream from the thyroid-specific transcription factor PAX8 (lowest P=1.1 × 10(-9)). This finding was replicated in an African-American sample of 4771 individuals (lowest P=9.3 × 10(-4)). The strongest combined association was at rs1823125 (P=1.5 × 10(-10), minor allele frequency 0.26 in the discovery sample, 0.12 in the replication sample), with each copy of the minor allele associated with a sleep duration 3.1 min longer per night. The alleles associated with longer sleep duration were associated in previous GWAS with a more favorable metabolic profile and a lower risk of attention deficit hyperactivity disorder. Understanding the mechanisms underlying these associations may help elucidate biological mechanisms influencing sleep duration and its association with psychiatric, metabolic and cardiovascular disease.

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

Genome-wide association studies (GWAS) have identified more than 100 genetic variants contributing to BMI, a measure of body size, or waist-to-hip ratio (adjusted for BMI, WHR_adjBMI), a measure of body shape. Body size and shape change as people grow older and these changes differ substantially between men and women. To systematically screen for age- and/or sex-specific effects of genetic variants on BMI and WHR_adjBMI, we performed meta-analyses of 114 studies (up to 320,485 individuals of European descent) with genome-wide chip and/or Metabochip data by the Genetic Investigation of Anthropometric Traits (GIANT) Consortium. Each study tested the association of up to ~2.8M SNPs with BMI and WHR_adjBMI in four strata (men ≤50y, men >50y, women ≤50y, women >50y) and summary statistics were combined in stratum-specific meta-analyses. We then screened for variants that showed age-specific effects (G x AGE), sex-specific effects (G x SEX) or age-specific effects that differed between men and women (G x AGE x SEX). For BMI, we identified 15 loci (11 previously established for main effects, four novel) that showed significant (FDR<5%) age-specific effects, of which 11 had larger effects in younger (<50y) than in older adults (≥50y). No sex-dependent effects were identified for BMI. For WHR_adjBMI, we identified 44 loci (27 previously established for main effects, 17 novel) with sex-specific effects, of which 28 showed larger effects in women than in men, five showed larger effects in men than in women, and 11 showed opposite effects between sexes. No age-dependent effects were identified for WHR_adjBMI. This is the first genome-wide interaction meta-analysis to report convincing evidence of age-dependent genetic effects on BMI. In addition, we confirm the sex-specificity of genetic effects on WHR_adjBMI. These results may provide further insights into the biology that underlies weight change with age or the sexually dimorphism of body shape.

Adult body size and body shape differ substantially between men and women and change over time. More than 100 genetic variants that influence body mass index (measure of body size) or waist-to-hip ratio (measure of body shape) have been identified. While there is evidence that some genetic loci affect body shape differently in men than in women, little is known about whether genetic effects differ in older compared to younger adults, and whether such changes differ between men and women. Therefore, we conducted a systematic genome-wide search, including 114 studies (>320,000 individuals), to specifically identify genetic loci with age- and or sex-dependent effects on body size and shape. We identified 15 loci of which the effect on BMI was different in older compared to younger adults, whereas we found no evidence for loci with different effects in men compared to women. The opposite was seen for body shape as we identified 44 loci of which the effect on waist-to-hip ratio differed between men and women, but no difference between younger and older adults were observed. Our observations may provide new insights into the biology that underlies weight change with age or the sexual dimorphism of body shape.

Single Nucleotide Polymorphism Heritability of Behavior Problems in Childhood: Genome-Wide Complex Trait Analysis

OBJECTIVE

Genetic factors contribute to individual differences in behavior problems. In children, genome-wide association studies (GWAS) have yielded the first suggestive results when aiming to identify genetic variants that explain heritability, but the proportion of genetic variance that can be attributed to common single nucleotide polymorphisms (SNPs) remains to be determined, as only a few studies have estimated SNP heritability, with diverging results.

METHOD

Genomic-relationship-matrix restricted maximum likelihood (GREML) as implemented in the software Genome-Wide Complex Trait Analysis (GCTA) was used to estimate SNP heritability (SNP h(2)) for multiple phenotypes within 4 broad domains of children's behavioral problems (attention-deficit/hyperactivity symptoms, internalizing, externalizing, and pervasive developmental problems) and cognitive function. We combined phenotype and genotype data from 2 independent, population-based Dutch cohorts, yielding a total number of 1,495 to 3,175 of 3-, 7-, and 9-year-old children.

RESULTS

Significant SNP heritability estimates were found for attention-deficit/hyperactivity symptoms (SNP h(2) = 0.37-0.71), externalizing problems (SNP h(2) = 0.44), and total problems (SNP h(2) = 0.18), rated by mother or teacher. Sensitivity analyses with exclusion of extreme cases and quantile normalization of the phenotype data decreased SNP h(2) as expected under genetic inheritance, but they remained statistically significant for most phenotypes.

CONCLUSION

We provide evidence of the influence of common SNPs on child behavior problems in an ethnically homogenous sample. These results support the continuation of large GWAS collaborative efforts to unravel the genetic basis of complex child behaviors.

Harmonising and linking biomedical and clinical data across disparate data archives to enable integrative cross-biobank research

A wealth of biospecimen samples are stored in modern globally distributed biobanks. Biomedical researchers worldwide need to be able to combine the available resources to improve the power of large-scale studies. A prerequisite for this effort is to be able to search and access phenotypic, clinical and other information about samples that are currently stored at biobanks in an integrated manner. However, privacy issues together with heterogeneous information systems and the lack of agreed-upon vocabularies have made specimen searching across multiple biobanks extremely challenging. We describe three case studies where we have linked samples and sample descriptions in order to facilitate global searching of available samples for research. The use cases include the ENGAGE (European Network for Genetic and Genomic Epidemiology) consortium comprising at least 39 cohorts, the SUMMIT (surrogate markers for micro- and macro-vascular hard endpoints for innovative diabetes tools) consortium and a pilot for data integration between a Swedish clinical health registry and a biobank. We used the Sample avAILability (SAIL) method for data linking: first, created harmonised variables and then annotated and made searchable information on the number of specimens available in individual biobanks for various phenotypic categories. By operating on this categorised availability data we sidestep many obstacles related to privacy that arise when handling real values and show that harmonised and annotated records about data availability across disparate biomedical archives provide a key methodological advance in pre-analysis exchange of information between biobanks, that is, during the project planning phase.

Directional dominance on stature and cognition in diverse human populations

Homozygosity has long been associated with rare, often devastating, Mendelian disorders, and Darwin was one of the first to recognize that inbreeding reduces evolutionary fitness. However, the effect of the more distant parental relatedness that is common in modern human populations is less well understood. Genomic data now allow us to investigate the effects of homozygosity on traits of public health importance by observing contiguous homozygous segments (runs of homozygosity), which are inferred to be homozygous along their complete length. Given the low levels of genome-wide homozygosity prevalent in most human populations, information is required on very large numbers of people to provide sufficient power. Here we use runs of homozygosity to study 16 health-related quantitative traits in 354,224 individuals from 102 cohorts, and find statistically significant associations between summed runs of homozygosity and four complex traits: height, forced expiratory lung volume in one second, general cognitive ability and educational attainment (P < 1 × 10(-300), 2.1 × 10(-6), 2.5 × 10(-10) and 1.8 × 10(-10), respectively). In each case, increased homozygosity was associated with decreased trait value, equivalent to the offspring of first cousins being 1.2 cm shorter and having 10 months' less education. Similar effect sizes were found across four continental groups and populations with different degrees of genome-wide homozygosity, providing evidence that homozygosity, rather than confounding, directly contributes to phenotypic variance. Contrary to earlier reports in substantially smaller samples, no evidence was seen of an influence of genome-wide homozygosity on blood pressure and low density lipoprotein cholesterol, or ten other cardio-metabolic traits. Since directional dominance is predicted for traits under directional evolutionary selection, this study provides evidence that increased stature and cognitive function have been positively selected in human evolution, whereas many important risk factors for late-onset complex diseases may not have been.

Heritability and Genome-Wide Association Studies for Hair Color in a Dutch Twin Family Based Sample

Hair color is one of the most visible and heritable traits in humans. Here, we estimated heritability by structural equation modeling (N = 20,142), and performed a genome wide association (GWA) analysis (N = 7091) and a GCTA study (N = 3340) on hair color within a large cohort of twins, their parents and siblings from the Netherlands Twin Register (NTR). Self-reported hair color was analyzed as five binary phenotypes, namely “blond versus non-blond”, “red versus non-red”, “brown versus non-brown”, “black versus non-black”, and “light versus dark”. The broad-sense heritability of hair color was estimated between 73% and 99% and the genetic component included non-additive genetic variance. Assortative mating for hair color was significant, except for red and black hair color. From GCTA analyses, at most 24.6% of the additive genetic variance in hair color was explained by 1000G well-imputed SNPs. Genome-wide association analysis for each hair color showed that SNPs in the MC1R region were significantly associated with red, brown and black hair, and also with light versus dark hair color. Five other known genes (HERC2, TPCN2, SLC24A4, IRF4, and KITLG) gave genome-wide significant hits for blond, brown and light versus dark hair color. We did not find and replicate any new loci for hair color.

Meta-analysis of Genome-wide Association Studies for Neuroticism, and the Polygenic Association With Major Depressive Disorder

IMPORTANCE

Neuroticism is a pervasive risk factor for psychiatric conditions. It genetically overlaps with major depressive disorder (MDD) and is therefore an important phenotype for psychiatric genetics. The Genetics of Personality Consortium has created a resource for genome-wide association analyses of personality traits in more than 63,000 participants (including MDD cases).

OBJECTIVES

To identify genetic variants associated with neuroticism by performing a meta-analysis of genome-wide association results based on 1000 Genomes imputation; to evaluate whether common genetic variants as assessed by single-nucleotide polymorphisms (SNPs) explain variation in neuroticism by estimating SNP-based heritability; and to examine whether SNPs that predict neuroticism also predict MDD.

DESIGN, SETTING, AND PARTICIPANTS

Genome-wide association meta-analysis of 30 cohorts with genome-wide genotype, personality, and MDD data from the Genetics of Personality Consortium. The study included 63,661 participants from 29 discovery cohorts and 9786 participants from a replication cohort. Participants came from Europe, the United States, or Australia. Analyses were conducted between 2012 and 2014.

MAIN OUTCOMES AND MEASURES

Neuroticism scores harmonized across all 29 discovery cohorts by item response theory analysis, and clinical MDD case-control status in 2 of the cohorts.

RESULTS

A genome-wide significant SNP was found on 3p14 in MAGI1 (rs35855737; P = 9.26 × 10-9 in the discovery meta-analysis). This association was not replicated (P = .32), but the SNP was still genome-wide significant in the meta-analysis of all 30 cohorts (P = 2.38 × 10-8). Common genetic variants explain 15% of the variance in neuroticism. Polygenic scores based on the meta-analysis of neuroticism in 27 cohorts significantly predicted neuroticism (1.09 × 10-12 < P < .05) and MDD (4.02 × 10-9 < P < .05) in the 2 other cohorts.

CONCLUSIONS AND RELEVANCE

This study identifies a novel locus for neuroticism. The variant is located in a known gene that has been associated with bipolar disorder and schizophrenia in previous studies. In addition, the study shows that neuroticism is influenced by many genetic variants of small effect that are either common or tagged by common variants. These genetic variants also influence MDD. Future studies should confirm the role of the MAGI1 locus for neuroticism and further investigate the association of MAGI1 and the polygenic association to a range of other psychiatric disorders that are phenotypically correlated with neuroticism.

Discovery and Fine-Mapping of Glycaemic and Obesity-Related Trait Loci Using High-Density Imputation

Reference panels from the 1000 Genomes (1000G) Project Consortium provide near complete coverage of common and low-frequency genetic variation with minor allele frequency ≥0.5% across European ancestry populations. Within the European Network for Genetic and Genomic Epidemiology (ENGAGE) Consortium, we have undertaken the first large-scale meta-analysis of genome-wide association studies (GWAS), supplemented by 1000G imputation, for four quantitative glycaemic and obesity-related traits, in up to 87,048 individuals of European ancestry. We identified two loci for body mass index (BMI) at genome-wide significance, and two for fasting glucose (FG), none of which has been previously reported in larger meta-analysis efforts to combine GWAS of European ancestry. Through conditional analysis, we also detected multiple distinct signals of association mapping to established loci for waist-hip ratio adjusted for BMI (RSPO3) and FG (GCK and G6PC2). The index variant for one association signal at the G6PC2 locus is a low-frequency coding allele, H177Y, which has recently been demonstrated to have a functional role in glucose regulation. Fine-mapping analyses revealed that the non-coding variants most likely to drive association signals at established and novel loci were enriched for overlap with enhancer elements, which for FG mapped to promoter and transcription factor binding sites in pancreatic islets, in particular. Our study demonstrates that 1000G imputation and genetic fine-mapping of common and low-frequency variant association signals at GWAS loci, integrated with genomic annotation in relevant tissues, can provide insight into the functional and regulatory mechanisms through which their effects on glycaemic and obesity-related traits are mediated.

Human genetic studies have demonstrated that quantitative human anthropometric and metabolic traits, including body mass index, waist-hip ratio, and plasma concentrations of glucose and insulin, are highly heritable, and are established risk factors for type 2 diabetes and cardiovascular diseases. Although many regions of the genome have been associated with these traits, the specific genes responsible have not yet been identified. By making use of advanced statistical “imputation” techniques applied to more than 87,000 individuals of European ancestry, and publicly available “reference panels” of more than 37 million genetic variants, we have been able to identify novel regions of the genome associated with these glycaemic and obesity-related traits and localise genes within these regions that are most likely to be causal. This improved understanding of the biological mechanisms underlying glycaemic and obesity-related traits is extremely important because it may advance drug development for downstream disease endpoints, ultimately leading to public health benefits.

The association between lower educational attainment and depression owing to shared genetic effects? Results in ~25,000 subjects

An association between lower educational attainment (EA) and an increased risk for depression has been confirmed in various western countries. This study examines whether pleiotropic genetic effects contribute to this association. Therefore, data were analyzed from a total of 9662 major depressive disorder (MDD) cases and 14,949 controls (with no lifetime MDD diagnosis) from the Psychiatric Genomics Consortium with additional Dutch and Estonian data. The association of EA and MDD was assessed with logistic regression in 15,138 individuals indicating a significantly negative association in our sample with an odds ratio for MDD 0.78 (0.75-0.82) per standard deviation increase in EA. With data of 884,105 autosomal common single-nucleotide polymorphisms (SNPs), three methods were applied to test for pleiotropy between MDD and EA: (i) genetic profile risk scores (GPRS) derived from training data for EA (independent meta-analysis on ~120,000 subjects) and MDD (using a 10-fold leave-one-out procedure in the current sample), (ii) bivariate genomic-relationship-matrix restricted maximum likelihood (GREML) and (iii) SNP effect concordance analysis (SECA). With these methods, we found (i) that the EA-GPRS did not predict MDD status, and MDD-GPRS did not predict EA, (ii) a weak negative genetic correlation with bivariate GREML analyses, but this correlation was not consistently significant, (iii) no evidence for concordance of MDD and EA SNP effects with SECA analysis. To conclude, our study confirms an association of lower EA and MDD risk, but this association was not because of measurable pleiotropic genetic effects, which suggests that environmental factors could be involved, for example, socioeconomic status.

Heritability, SNP- and Gene-Based Analyses of Cannabis Use Initiation and Age at Onset

Prior searches for genetic variants (GVs) implicated in initiation of cannabis use have been limited to common single nucleotide polymorphisms (SNPs) typed in HapMap samples. Denser SNPs are now available with the completion of the 1000 Genomes and the Genome of the Netherlands projects. More densely distributed SNPs are expected to track the causal variants better. Therefore we extend the search for variants implicated in early stages of cannabis use to previously untagged common and low-frequency variants. We run heritability, SNP and gene-based analyses of initiation and age at onset. This is the first genome-wide study of age at onset to date. Using GCTA and a sample of distantly related individuals from the Netherlands Twin Register, we estimated that the currently measured (and tagged) SNPs collectively explain 25 % of the variance in initiation (SE = 0.088; P = 0.0016). Chromosomes 4 and 18, previously linked with cannabis use and other addiction phenotypes, account for the largest amount of variance in initiation (6.8 %, SE = 0.025, P = 0.002 and 3.6 %, SE = 0.01, P = 0.012, respectively). No individual SNP- or gene-based test reached genomewide significance in the initiation or age at onset analyses. Our study detected association signal in the currently measured SNPs. A comparison with prior SNP-heritability estimates suggests that at least part of the signal is likely coming from previously untyped common and low frequency variants. Our results do not rule out the contribution of rare variants of larger effect—a plausible source of the difference between the twin-based heritability estimate and that from GCTA. The causal variants are likely of very small effect (i.e., <1 % explained variance) and are uniformly distributed over the genome in proportion to chromosomes’ length. Similar to other complex traits and diseases, detecting such small effects is to be expected in sufficiently large samples.

Genetic analyses benefit from using less heterogeneous phenotypes: an illustration with the hospital anxiety and depression scale (HADS)

Phenotypic heterogeneity of depression has been cited as one of the causes of the limited success to detect genetic variants in genome-wide studies. The 7-item Hospital Anxiety and Depression Scale (HADS-D) was developed to detect depression in individuals with physical health problems. An initial psychometric analysis showed that a short version ("HADS-4") is less heterogeneous and hence more reliable than the full scale, and correlates equally strong with a DSM-oriented depression scale. We compared the HADS-D and the HADS-4 to assess the benefits of using less heterogeneous phenotype measures in genetic analyses. We compared HADS-D and HADS-4 in three separate analyses: (1) twin- and family-based heritability estimation, (2) SNP-based heritability estimation using the software GCTA, and (3) a genome-wide association study (GWAS). The twin study resulted in heritability estimates between 18% and 25%, with additive genetic variance being the largest component. There was also evidence for assortative mating and a dominance component of genetic variance, with HADS-4 having slightly lower estimates of assortment. Importantly, when estimating heritability from SNPs, the HADS-D did not show a significant genetic variance component, while for the HADS-4, a statistically significant amount of heritability was estimated. Moreover, the HADS-4 had substantially more SNPs with small P-values in the GWAS analysis than did the HADS-D. Our results underline the benefits of using more homogeneous phenotypes in psychiatric genetic analyses. Homogeneity can be increased by focusing on core symptoms of disorders, thus reducing the noise in aggregate phenotypes caused by substantially different symptom profiles.

Genome of The Netherlands population-specific imputations identify an ABCA6 variant associated with cholesterol levels

Variants associated with blood lipid levels may be population-specific. To identify low-frequency variants associated with this phenotype, population-specific reference panels may be used. Here we impute nine large Dutch biobanks (~35,000 samples) with the population-specific reference panel created by the Genome of the Netherlands Project and perform association testing with blood lipid levels. We report the discovery of five novel associations at four loci (P value <6.61 × 10⁻⁴), including a rare missense variant in ABCA6 (rs77542162, p.Cys1359Arg, frequency 0.034), which is predicted to be deleterious. The frequency of this ABCA6 variant is 3.65-fold increased in the Dutch and its effect (β_LDL-C=0.135, β_TC=0.140) is estimated to be very similar to those observed for single variants in well-known lipid genes, such as LDLR.

Frequencies of rare variants fluctuate over populations, hampering gene discovery. Here the authors use a population-specific reference panel, the Genome of the Netherlands, to discover four novel loci involved in lipid metabolism, including an exonic variant in ABCA6.

In Silico Post Genome-Wide Association Studies Analysis of C-Reactive Protein Loci Suggests an Important Role for Interferons

BACKGROUND

Genome-wide association studies (GWASs) have successfully identified several single nucleotide polymorphisms (SNPs) associated with serum levels of C-reactive protein (CRP). An important limitation of GWASs is that the identified variants merely flag the nearby genomic region and do not necessarily provide a direct link to the biological mechanisms underlying their corresponding phenotype. Here we apply a bioinformatics-based approach to uncover the functional characteristics of the 18 SNPs that had previously been associated with CRP at a genome-wide significant level.

METHODS AND RESULTS

In the first phase of in silico sequencing, we explore the vicinity of GWAS SNPs to identify all linked variants. In the second phase of expression quantitative trait loci analysis, we attempt to identify all nearby genes whose expression levels are associated with the corresponding GWAS SNPs. These 2 phases generate several relevant genes that serve as input to the next phase of functional network analysis. Our in silico sequencing analysis using 1000 Genomes Project data identified 7 nonsynonymous SNPs, which are in moderate to high linkage disequilibrium (r(2)>0.5) with the GWAS SNPs. Our expression quantitative trait loci analysis, which was based on one of the largest single data sets of genome-wide expression probes (n>5000) identified 23 significantly associated expression probes belonging to 15 genes (false discovery rate <0.01). The final phase of functional network analysis revealed 93 significantly enriched biological processes (false discovery rate <0.01).

CONCLUSIONS

Our post-GWAS analysis of CRP GWAS SNPs confirmed the previously known overlap between CRP and lipids biology. Additionally, it suggested an important role for interferons in the metabolism of CRP.

Educational attainment influences levels of homozygosity through migration and assortative mating

Individuals with a higher education are more likely to migrate, increasing the chance of meeting a spouse with a different ancestral background. In this context, the presence of strong educational assortment can result in greater ancestry differences within more educated spouse pairs, while less educated individuals are more likely to mate with someone with whom they share more ancestry. We examined the association between educational attainment and F_roh (= the proportion of the genome consisting of runs of homozygosity [ROHs]) in ~2,000 subjects of Dutch ancestry. The subjects’ own educational attainment showed a nominally significant negative association with F_roh (p = .045), while the contribution of parental education to offspring F_roh was highly significant (father: p < 10^-5; mother: p = 9×10^-5), with more educated parents having offspring with fewer ROHs. This association was significantly and fully mediated by the physical distance between parental birthplaces (paternal education: p_mediation = 2.4 × 10^-4; maternal education: p_mediation = 2.3 × 10^-4), which itself was also significantly associated with F_roh (p = 9 × 10^-5). Ancestry-informative principal components from the offspring showed a significantly decreasing association with geography as parental education increased, consistent with the significantly higher migration rates among more educated parents. Parental education also showed a high spouse correlation (Spearman’s ρ = .66, p = 3 × 10^-262). We show that less educated parents are less likely to mate with the more mobile parents with a higher education, creating systematic differences in homozygosity due to ancestry differences not directly captured by ancestry-informative principal components (PCs). Understanding how behaviors influence the genomic structure of a population is highly valuable for studies on the genetic etiology of behavioral, cognitive, and social traits.

Intelligence: shared genetic basis between Mendelian disorders and a polygenic trait

Multiple inquiries into the genetic etiology of human traits indicated an overlap between genes underlying monogenic disorders (eg, skeletal growth defects) and those affecting continuous variability of related quantitative traits (eg, height). Extending the idea of a shared genetic basis between a Mendelian disorder and a classic polygenic trait, we performed an association study to examine the effect of 43 genes implicated in autosomal recessive cognitive disorders on intelligence in an unselected Dutch population (N=1316). Using both single-nucleotide polymorphism (SNP)- and gene-based association testing, we detected an association between intelligence and the genes of interest, with genes ELP2, TMEM135, PRMT10, and RGS7 showing the strongest associations. This is a demonstration of the relevance of genes implicated in monogenic disorders of intelligence to normal-range intelligence, and a corroboration of the utility of employing knowledge on monogenic disorders in identifying the genetic variability underlying complex traits.

Joint analysis of psychiatric disorders increases accuracy of risk prediction for schizophrenia, bipolar disorder, and major depressive disorder

Genetic risk prediction has several potential applications in medical research and clinical practice and could be used, for example, to stratify a heterogeneous population of patients by their predicted genetic risk. However, for polygenic traits, such as psychiatric disorders, the accuracy of risk prediction is low. Here we use a multivariate linear mixed model and apply multi-trait genomic best linear unbiased prediction for genetic risk prediction. This method exploits correlations between disorders and simultaneously evaluates individual risk for each disorder. We show that the multivariate approach significantly increases the prediction accuracy for schizophrenia, bipolar disorder, and major depressive disorder in the discovery as well as in independent validation datasets. By grouping SNPs based on genome annotation and fitting multiple random effects, we show that the prediction accuracy could be further improved. The gain in prediction accuracy of the multivariate approach is equivalent to an increase in sample size of 34% for schizophrenia, 68% for bipolar disorder, and 76% for major depressive disorders using single trait models. Because our approach can be readily applied to any number of GWAS datasets of correlated traits, it is a flexible and powerful tool to maximize prediction accuracy. With current sample size, risk predictors are not useful in a clinical setting but already are a valuable research tool, for example in experimental designs comparing cases with high and low polygenic risk.

Common genetic variants influence human subcortical brain structures

The highly complex structure of the human brain is strongly shaped by genetic influences. Subcortical brain regions form circuits with cortical areas to coordinate movement, learning, memory and motivation, and altered circuits can lead to abnormal behaviour and disease. To investigate how common genetic variants affect the structure of these brain regions, here we conduct genome-wide association studies of the volumes of seven subcortical regions and the intracranial volume derived from magnetic resonance images of 30,717 individuals from 50 cohorts. We identify five novel genetic variants influencing the volumes of the putamen and caudate nucleus. We also find stronger evidence for three loci with previously established influences on hippocampal volume and intracranial volume. These variants show specific volumetric effects on brain structures rather than global effects across structures. The strongest effects were found for the putamen, where a novel intergenic locus with replicable influence on volume (rs945270; P = 1.08 × 10(-33); 0.52% variance explained) showed evidence of altering the expression of the KTN1 gene in both brain and blood tissue. Variants influencing putamen volume clustered near developmental genes that regulate apoptosis, axon guidance and vesicle transport. Identification of these genetic variants provides insight into the causes of variability in human brain development, and may help to determine mechanisms of neuropsychiatric dysfunction.

No evidence for genome-wide interactions on plasma fibrinogen by smoking, alcohol consumption and body mass index: results from meta-analyses of 80,607 subjects

Plasma fibrinogen is an acute phase protein playing an important role in the blood coagulation cascade having strong associations with smoking, alcohol consumption and body mass index (BMI). Genome-wide association studies (GWAS) have identified a variety of gene regions associated with elevated plasma fibrinogen concentrations. However, little is yet known about how associations between environmental factors and fibrinogen might be modified by genetic variation. Therefore, we conducted large-scale meta-analyses of genome-wide interaction studies to identify possible interactions of genetic variants and smoking status, alcohol consumption or BMI on fibrinogen concentration. The present study included 80,607 subjects of European ancestry from 22 studies. Genome-wide interaction analyses were performed separately in each study for about 2.6 million single nucleotide polymorphisms (SNPs) across the 22 autosomal chromosomes. For each SNP and risk factor, we performed a linear regression under an additive genetic model including an interaction term between SNP and risk factor. Interaction estimates were meta-analysed using a fixed-effects model. No genome-wide significant interaction with smoking status, alcohol consumption or BMI was observed in the meta-analyses. The most suggestive interaction was found for smoking and rs10519203, located in the LOC123688 region on chromosome 15, with a p value of 6.2×10⁻⁸. This large genome-wide interaction study including 80,607 participants found no strong evidence of interaction between genetic variants and smoking status, alcohol consumption or BMI on fibrinogen concentrations. Further studies are needed to yield deeper insight in the interplay between environmental factors and gene variants on the regulation of fibrinogen concentrations.

Identifying genetic variants for heart rate variability in the acetylcholine pathway

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.

Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche

Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10(-8)) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition.

The ENIGMA Consortium: large-scale collaborative analyses of neuroimaging and genetic data

The Enhancing NeuroImaging Genetics through Meta-Analysis (ENIGMA) Consortium is a collaborative network of researchers working together on a range of large-scale studies that integrate data from 70 institutions worldwide. Organized into Working Groups that tackle questions in neuroscience, genetics, and medicine, ENIGMA studies have analyzed neuroimaging data from over 12,826 subjects. In addition, data from 12,171 individuals were provided by the CHARGE consortium for replication of findings, in a total of 24,997 subjects. By meta-analyzing results from many sites, ENIGMA has detected factors that affect the brain that no individual site could detect on its own, and that require larger numbers of subjects than any individual neuroimaging study has currently collected. ENIGMA's first project was a genome-wide association study identifying common variants in the genome associated with hippocampal volume or intracranial volume. Continuing work is exploring genetic associations with subcortical volumes (ENIGMA2) and white matter microstructure (ENIGMA-DTI). Working groups also focus on understanding how schizophrenia, bipolar illness, major depression and attention deficit/hyperactivity disorder (ADHD) affect the brain. We review the current progress of the ENIGMA Consortium, along with challenges and unexpected discoveries made on the way.

The contribution of the functional IL6R polymorphism rs2228145, eQTLs and other genome-wide SNPs to the heritability of plasma sIL-6R levels

The non-synonymous SNP rs2228145 in the IL6R gene on chromosome 1q21.3 is associated with a wide range of common diseases, including asthma, rheumatoid arthritis, type 1 diabetes and coronary heart disease. We examined the contribution of this functional IL6R gene polymorphism rs2228145 versus other genome-wide SNPs to the variance of sIL-6R levels in blood plasma in a large population-based sample (N ~5,000), and conducted an expression QTL analysis to identify SNPs associated with IL6R gene expression. Based on data from 2,360 twin families, the broad heritability of sIL-6R was estimated at 72 and 51% of the total variance was explained by the functional SNP rs2228145. Converging findings from GWAS, linkage, and GCTA analyses indicate that additional variance of sIL-6R levels can be explained by other variants in the IL6R region, including variants at the 3'-end of IL6R tagged by rs60760897 that are associated with IL6R RNA expression.

Genome-wide association meta-analysis of human longevity identifies a novel locus conferring survival beyond 90 years of age

The genetic contribution to the variation in human lifespan is ∼25%. Despite the large number of identified disease-susceptibility loci, it is not known which loci influence population mortality. We performed a genome-wide association meta-analysis of 7729 long-lived individuals of European descent (≥85 years) and 16 121 younger controls (<65 years) followed by replication in an additional set of 13 060 long-lived individuals and 61 156 controls. In addition, we performed a subset analysis in cases aged ≥90 years. We observed genome-wide significant association with longevity, as reflected by survival to ages beyond 90 years, at a novel locus, rs2149954, on chromosome 5q33.3 (OR = 1.10, P = 1.74 × 10⁻⁸). We also confirmed association of rs4420638 on chromosome 19q13.32 (OR = 0.72, P = 3.40 × 10⁻³⁶), representing the TOMM40/APOE/APOC1 locus. In a prospective meta-analysis (n = 34 103), the minor allele of rs2149954 (T) on chromosome 5q33.3 associates with increased survival (HR = 0.95, P = 0.003). This allele has previously been reported to associate with low blood pressure in middle age. Interestingly, the minor allele (T) associates with decreased cardiovascular mortality risk, independent of blood pressure. We report on the first GWAS-identified longevity locus on chromosome 5q33.3 influencing survival in the general European population. The minor allele of this locus associates with low blood pressure in middle age, although the contribution of this allele to survival may be less dependent on blood pressure. Hence, the pleiotropic mechanisms by which this intragenic variation contributes to lifespan regulation have to be elucidated.

The dopaminergic reward system and leisure time exercise behavior: a candidate allele study

PURPOSE

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.

MATERIALS AND METHODS

Data on exercise activities and at least one SNP/VNTR were available for 8,768 individuals aged 7 to 50 years old that were part of the Netherlands Twin Register. Exercise behavior was quantified as weekly metabolic equivalents of task (MET) spent on exercise activities. Mixed models were fitted in SPSS with genetic relatedness as a random effect.

RESULTS

None of the genetic variants were associated with exercise behavior (P>.02), despite sufficient power to detect small effects.

DISCUSSION AND CONCLUSIONS

We did not confirm that allelic variants involved in dopaminergic function play a role in creating individual differences in exercise behavior. A plea is made for large genome-wide association studies to unravel the genetic pathways that affect this health-enhancing behavior.

Genetic risk score analysis indicates migraine with and without comorbid depression are genetically different disorders

Migraine and major depressive disorder (MDD) are comorbid, moderately heritable and to some extent influenced by the same genes. In a previous paper, we suggested the possibility of causality (one trait causing the other) underlying this comorbidity. We present a new application of polygenic (genetic risk) score analysis to investigate the mechanisms underlying the genetic overlap of migraine and MDD. Genetic risk scores were constructed based on data from two discovery samples in which genome-wide association analyses (GWA) were performed for migraine and MDD, respectively. The Australian Twin Migraine GWA study (N = 6,350) included 2,825 migraine cases and 3,525 controls, 805 of whom met the diagnostic criteria for MDD. The RADIANT GWA study (N = 3,230) included 1,636 MDD cases and 1,594 controls. Genetic risk scores for migraine and for MDD were used to predict pure and comorbid forms of migraine and MDD in an independent Dutch target sample (NTR-NESDA, N = 2,966), which included 1,476 MDD cases and 1,058 migraine cases (723 of these individuals had both disorders concurrently). The observed patterns of prediction suggest that the 'pure' forms of migraine and MDD are genetically distinct disorders. The subgroup of individuals with comorbid MDD and migraine were genetically most similar to MDD patients. These results indicate that in at least a subset of migraine patients with MDD, migraine may be a symptom or consequence of MDD.

Sex differences in the human peripheral blood transcriptome

BACKGROUND

Genomes of men and women differ in only a limited number of genes located on the sex chromosomes, whereas the transcriptome is far more sex-specific. Identification of sex-biased gene expression will contribute to understanding the molecular basis of sex-differences in complex traits and common diseases.

RESULTS

Sex differences in the human peripheral blood transcriptome were characterized using microarrays in 5,241 subjects, accounting for menopause status and hormonal contraceptive use. Sex-specific expression was observed for 582 autosomal genes, of which 57.7% was upregulated in women (female-biased genes). Female-biased genes were enriched for several immune system GO categories, genes linked to rheumatoid arthritis (16%) and genes regulated by estrogen (18%). Male-biased genes were enriched for genes linked to renal cancer (9%). Sex-differences in gene expression were smaller in postmenopausal women, larger in women using hormonal contraceptives and not caused by sex-specific eQTLs, confirming the role of estrogen in regulating sex-biased genes.

CONCLUSIONS

This study indicates that sex-bias in gene expression is extensive and may underlie sex-differences in the prevalence of common diseases.

Further confirmation of the association between anxiety and CTNND2: replication in humans

The rat genome sequencing and mapping consortium found evidence for an association between the catenin-δ2 gene (CTNND2) and anxious behaviour. We replicated these results in humans by carrying out a genetic association test in patients with panic disorder, social phobia, generalized anxiety disorder and/or agoraphobia (N = 1714) and controls (N = 4125). We further explored the association between CTNND2 and other psychiatric disorders based on publicly available genome-wide association results. A gene-based test showed that single nucleotide polymorphisms (SNPs) in CTNND2 have a significantly increased signal (P < 1e(-5) ) and decreased P-values. Single nucleotide polymorphism rs1012176 showed the strongest association with any anxiety disorder (odds ratio: 0.8128, SE = 0.063, P = 0.00099), but this effect was not significant after correction for multiple testing. In available genome-wide association results from the Psychiatric Genomics Consortium we found that SNPs in CTNND2 collectively showed an increased signal for schizophrenia (P < 1e(-5) ) and major depressive disorder (P < 1e(-5) ), but not for bipolar disorder. These signals remained significant after correction for potential confounders. The association between CTNND2 and anxiety was not strong enough to be picked up in the current generation of human genome-wide analyses, indicating the usefulness of and need for animal genetic studies to identify candidate genes for further study in human samples.

Inference of the genetic architecture underlying BMI and height with the use of 20,240 sibling pairs

Evidence that complex traits are highly polygenic has been presented by population-based genome-wide association studies (GWASs) through the identification of many significant variants, as well as by family-based de novo sequencing studies indicating that several traits have a large mutational target size. Here, using a third study design, we show results consistent with extreme polygenicity for body mass index (BMI) and height. On a sample of 20,240 siblings (from 9,570 nuclear families), we used a within-family method to obtain narrow-sense heritability estimates of 0.42 (SE = 0.17, p = 0.01) and 0.69 (SE = 0.14, p = 6 × 10(-)(7)) for BMI and height, respectively, after adjusting for covariates. The genomic inflation factors from locus-specific linkage analysis were 1.69 (SE = 0.21, p = 0.04) for BMI and 2.18 (SE = 0.21, p = 2 × 10(-10)) for height. This inflation is free of confounding and congruent with polygenicity, consistent with observations of ever-increasing genomic-inflation factors from GWASs with large sample sizes, implying that those signals are due to true genetic signals across the genome rather than population stratification. We also demonstrate that the distribution of the observed test statistics is consistent with both rare and common variants underlying a polygenic architecture and that previous reports of linkage signals in complex traits are probably a consequence of polygenic architecture rather than the segregation of variants with large effects. The convergent empirical evidence from GWASs, de novo studies, and within-family segregation implies that family-based sequencing studies for complex traits require very large sample sizes because the effects of causal variants are small on average.

Meta-analysis of gene-level associations for rare variants based on single-variant statistics

Meta-analysis of genome-wide association studies (GWASs) has led to the discoveries of many common variants associated with complex human diseases. There is a growing recognition that identifying "causal" rare variants also requires large-scale meta-analysis. The fact that association tests with rare variants are performed at the gene level rather than at the variant level poses unprecedented challenges in the meta-analysis. First, different studies may adopt different gene-level tests, so the results are not compatible. Second, gene-level tests require multivariate statistics (i.e., components of the test statistic and their covariance matrix), which are difficult to obtain. To overcome these challenges, we propose to perform gene-level tests for rare variants by combining the results of single-variant analysis (i.e., p values of association tests and effect estimates) from participating studies. This simple strategy is possible because of an insight that multivariate statistics can be recovered from single-variant statistics, together with the correlation matrix of the single-variant test statistics, which can be estimated from one of the participating studies or from a publicly available database. We show both theoretically and numerically that the proposed meta-analysis approach provides accurate control of the type I error and is as powerful as joint analysis of individual participant data. This approach accommodates any disease phenotype and any study design and produces all commonly used gene-level tests. An application to the GWAS summary results of the Genetic Investigation of ANthropometric Traits (GIANT) consortium reveals rare and low-frequency variants associated with human height. The relevant software is freely available.

Association of adiposity genetic variants with menarche timing in 92,105 women of European descent

Obesity is of global health concern. There are well-described inverse relationships between female pubertal timing and obesity. Recent genome-wide association studies of age at menarche identified several obesity-related variants. Using data from the ReproGen Consortium, we employed meta-analytical techniques to estimate the associations of 95 a priori and recently identified obesity-related (body mass index (weight (kg)/height (m)(2)), waist circumference, and waist:hip ratio) single-nucleotide polymorphisms (SNPs) with age at menarche in 92,116 women of European descent from 38 studies (1970-2010), in order to estimate associations between genetic variants associated with central or overall adiposity and pubertal timing in girls. Investigators in each study performed a separate analysis of associations between the selected SNPs and age at menarche (ages 9-17 years) using linear regression models and adjusting for birth year, site (as appropriate), and population stratification. Heterogeneity of effect-measure estimates was investigated using meta-regression. Six novel associations of body mass index loci with age at menarche were identified, and 11 adiposity loci previously reported to be associated with age at menarche were confirmed, but none of the central adiposity variants individually showed significant associations. These findings suggest complex genetic relationships between menarche and overall obesity, and to a lesser extent central obesity, in normal processes of growth and development.

Genome-wide meta-analysis identifies new susceptibility loci for migraine

Migraine is the most common brain disorder, affecting approximately 14% of the adult population, but its molecular mechanisms are poorly understood. We report the results of a meta-analysis across 29 genome-wide association studies, including a total of 23,285 individuals with migraine (cases) and 95,425 population-matched controls. We identified 12 loci associated with migraine susceptibility (P<5×10(-8)). Five loci are new: near AJAP1 at 1p36, near TSPAN2 at 1p13, within FHL5 at 6q16, within C7orf10 at 7p14 and near MMP16 at 8q21. Three of these loci were identified in disease subgroup analyses. Brain tissue expression quantitative trait locus analysis suggests potential functional candidate genes at four loci: APOA1BP, TBC1D7, FUT9, STAT6 and ATP5B.

GWAS of 126,559 individuals identifies genetic variants associated with educational attainment

A genome-wide association study (GWAS) of educational attainment was conducted in a discovery sample of 101,069 individuals and a replication sample of 25,490. Three independent single-nucleotide polymorphisms (SNPs) are genome-wide significant (rs9320913, rs11584700, rs4851266), and all three replicate. Estimated effects sizes are small (coefficient of determination R(2) ≈ 0.02%), approximately 1 month of schooling per allele. A linear polygenic score from all measured SNPs accounts for ≈2% of the variance in both educational attainment and cognitive function. Genes in the region of the loci have previously been associated with health, cognitive, and central nervous system phenotypes, and bioinformatics analyses suggest the involvement of the anterior caudate nucleus. These findings provide promising candidate SNPs for follow-up work, and our effect size estimates can anchor power analyses in social-science genetics.

Sex-stratified genome-wide association studies including 270,000 individuals show sexual dimorphism in genetic loci for anthropometric traits

Given the anthropometric differences between men and women and previous evidence of sex-difference in genetic effects, we conducted a genome-wide search for sexually dimorphic associations with height, weight, body mass index, waist circumference, hip circumference, and waist-to-hip-ratio (133,723 individuals) and took forward 348 SNPs into follow-up (additional 137,052 individuals) in a total of 94 studies. Seven loci displayed significant sex-difference (FDR<5%), including four previously established (near GRB14/COBLL1, LYPLAL1/SLC30A10, VEGFA, ADAMTS9) and three novel anthropometric trait loci (near MAP3K1, HSD17B4, PPARG), all of which were genome-wide significant in women (P<5×10(-8)), but not in men. Sex-differences were apparent only for waist phenotypes, not for height, weight, BMI, or hip circumference. Moreover, we found no evidence for genetic effects with opposite directions in men versus women. The PPARG locus is of specific interest due to its role in diabetes genetics and therapy. Our results demonstrate the value of sex-specific GWAS to unravel the sexually dimorphic genetic underpinning of complex traits.

Genome-wide meta-analysis identifies 11 new loci for anthropometric traits and provides insights into genetic architecture

Approaches exploiting trait distribution extremes may be used to identify loci associated with common traits, but it is unknown whether these loci are generalizable to the broader population. In a genome-wide search for loci associated with the upper versus the lower 5th percentiles of body mass index, height and waist-to-hip ratio, as well as clinical classes of obesity, including up to 263,407 individuals of European ancestry, we identified 4 new loci (IGFBP4, H6PD, RSRC1 and PPP2R2A) influencing height detected in the distribution tails and 7 new loci (HNF4G, RPTOR, GNAT2, MRPS33P4, ADCY9, HS6ST3 and ZZZ3) for clinical classes of obesity. Further, we find a large overlap in genetic structure and the distribution of variants between traits based on extremes and the general population and little etiological heterogeneity between obesity subgroups.

Identification of heart rate-associated loci and their effects on cardiac conduction and rhythm disorders

Elevated resting heart rate is associated with greater risk of cardiovascular disease and mortality. In a 2-stage meta-analysis of genome-wide association studies in up to 181,171 individuals, we identified 14 new loci associated with heart rate and confirmed associations with all 7 previously established loci. Experimental downregulation of gene expression in Drosophila melanogaster and Danio rerio identified 20 genes at 11 loci that are relevant for heart rate regulation and highlight a role for genes involved in signal transmission, embryonic cardiac development and the pathophysiology of dilated cardiomyopathy, congenital heart failure and/or sudden cardiac death. In addition, genetic susceptibility to increased heart rate is associated with altered cardiac conduction and reduced risk of sick sinus syndrome, and both heart rate-increasing and heart rate-decreasing variants associate with risk of atrial fibrillation. Our findings provide fresh insights into the mechanisms regulating heart rate and identify new therapeutic targets.

The molecular genetic architecture of self-employment

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.

The use of imputed sibling genotypes in sibship-based association analysis: on modeling alternatives, power and model misspecification

When phenotypic, but no genotypic data are available for relatives of participants in genetic association studies, previous research has shown that family-based imputed genotypes can boost the statistical power when included in such studies. Here, using simulations, we compared the performance of two statistical approaches suitable to model imputed genotype data: the mixture approach, which involves the full distribution of the imputed genotypes and the dosage approach, where the mean of the conditional distribution features as the imputed genotype. Simulations were run by varying sibship size, size of the phenotypic correlations among siblings, imputation accuracy and minor allele frequency of the causal SNP. Furthermore, as imputing sibling data and extending the model to include sibships of size two or greater requires modeling the familial covariance matrix, we inquired whether model misspecification affects power. Finally, the results obtained via simulations were empirically verified in two datasets with continuous phenotype data (height) and with a dichotomous phenotype (smoking initiation). Across the settings considered, the mixture and the dosage approach are equally powerful and both produce unbiased parameter estimates. In addition, the likelihood-ratio test in the linear mixed model appears to be robust to the considered misspecification in the background covariance structure, given low to moderate phenotypic correlations among siblings. Empirical results show that the inclusion in association analysis of imputed sibling genotypes does not always result in larger test statistic. The actual test statistic may drop in value due to small effect sizes. That is, if the power benefit is small, that the change in distribution of the test statistic under the alternative is relatively small, the probability is greater of obtaining a smaller test statistic. As the genetic effects are typically hypothesized to be small, in practice, the decision on whether family-based imputation could be used as a means to increase power should be informed by prior power calculations and by the consideration of the background correlation.

Strong effects of environmental factors on prevalence and course of major depressive disorder are not moderated by 5-HTTLPR polymorphisms in a large Dutch sample

BACKGROUND

There is ongoing interest in the possible interaction of the serotonin-transporter-linked polymorphic region (5-HTTLPR) with environmental factors in determining Major Depressive Disorder (MDD). The current study contributes to this research area by comprehensively examining the interaction-effects and direct-effects of 5-HTTLPR and five environmental factors on MDD prevalence and course in a well-characterized longitudinal sample.

METHODS

The sample consisted of 1625 patients with a CIDI-confirmed diagnosis of MDD and 1698 screened controls from the Netherlands. Four MDD outcomes were studied as dependent variables: one main MDD prevalence-outcome (all MDD), two more severe MDD prevalence-outcomes (suicidal and chronic MDD), and one MDD course outcome (chronic versus non-chronic MDD). Because SNP rs25531 modifies the effect of 5-HTTLPR, haplotypes of 5-HTTLPR and rs25531 were measured. For the four MDD outcome measures, we examined the direct effects of 5-HTTLPR/rs25531-haplotypes, five environmental factors (lifetime and recent stressful life-events, sexual abuse, low educational attainment, and childhood trauma) and their interaction in logistic regression models.

RESULTS

The environmental factors had large and consistent effects on all four MDD outcomes, including course of MDD. The 5-HTTLPR/rs25531-haplotype had a suggestive effect on course of MDD, but not on presence of MDD. Gene-by-environment interaction was significant (<0.05) for one of the 20 tests performed, which is not more than expected by chance.

LIMITATIONS

Environmental factors were not assessed before the onset of MDD.

CONCLUSIONS

Environmental factors had a strong impact on the presence and course of MDD, but no evidence for gene-by-environment interaction was found.

Refining genome-wide linkage intervals using a meta-analysis of genome-wide association studies identifies loci influencing personality dimensions

Personality traits are complex phenotypes related to psychosomatic health. Individually, various gene finding methods have not achieved much success in finding genetic variants associated with personality traits. We performed a meta-analysis of four genome-wide linkage scans (N=6149 subjects) of five basic personality traits assessed with the NEO Five-Factor Inventory. We compared the significant regions from the meta-analysis of linkage scans with the results of a meta-analysis of genome-wide association studies (GWAS) (N∼17 000). We found significant evidence of linkage of neuroticism to chromosome 3p14 (rs1490265, LOD=4.67) and to chromosome 19q13 (rs628604, LOD=3.55); of extraversion to 14q32 (ATGG002, LOD=3.3); and of agreeableness to 3p25 (rs709160, LOD=3.67) and to two adjacent regions on chromosome 15, including 15q13 (rs970408, LOD=4.07) and 15q14 (rs1055356, LOD=3.52) in the individual scans. In the meta-analysis, we found strong evidence of linkage of extraversion to 4q34, 9q34, 10q24 and 11q22, openness to 2p25, 3q26, 9p21, 11q24, 15q26 and 19q13 and agreeableness to 4q34 and 19p13. Significant evidence of association in the GWAS was detected between openness and rs677035 at 11q24 (P-value=2.6 × 10(-06), KCNJ1). The findings of our linkage meta-analysis and those of the GWAS suggest that 11q24 is a susceptible locus for openness, with KCNJ1 as the possible candidate gene.