Distinct biological signature and modifiable risk factors underlie the comorbidity between major depressive disorder and cardiovascular disease.

Major depressive disorder (MDD) and cardiovascular disease (CVD) are often comorbid, resulting in excess morbidity and mortality. Here we show that CVDs share most of their genetic risk factors with MDD. Multivariate genome-wide association analysis of shared genetic liability between MDD and atherosclerotic CVD revealed seven loci and distinct patterns of tissue and brain cell-type enrichments, suggesting the involvement of the thalamus. Part of the genetic overlap was explained by shared inflammatory, metabolic and psychosocial or lifestyle risk factors. Our data indicated causal effects of genetic liability to MDD on CVD risk, but not from most CVDs to MDD, and showed that the causal effects were partly explained by metabolic and psychosocial or lifestyle factors. The distinct signature of MDD-atherosclerotic CVD comorbidity suggests an immunometabolic subtype of MDD that is more strongly associated with CVD than overall MDD. In summary, we identified biological mechanisms underlying MDD-CVD comorbidity and modifiable risk factors for prevention of CVD in individuals with MDD.

Distinct genomic signatures and modifiable risk factors underly the comorbidity between major depressive disorder and cardiovascular disease.

Major depressive disorder (MDD) and cardiovascular disease (CVD) are often comorbid, resulting in excess morbidity and mortality. Using genomic data, this study elucidates biological mechanisms, key risk factors, and causal pathways underlying their comorbidity. We show that CVDs share a large proportion of their genetic risk factors with MDD. Multivariate genome-wide association analysis of the shared genetic liability between MDD and atherosclerotic CVD (ASCVD) revealed seven novel loci and distinct patterns of tissue and brain cell-type enrichments, suggesting a role for the thalamus. Part of the genetic overlap was explained by shared inflammatory, metabolic, and psychosocial/lifestyle risk factors. Finally, we found support for causal effects of genetic liability to MDD on CVD risk, but not from most CVDs to MDD, and demonstrated that the causal effects were partly explained by metabolic and psychosocial/lifestyle factors. The distinct signature of MDD-ASCVD comorbidity aligns with the idea of an immunometabolic sub-type of MDD more strongly associated with CVD than overall MDD. In summary, we identify plausible biological mechanisms underlying MDD-CVD comorbidity, as well as key modifiable risk factors for prevention of CVD in individuals with MDD.

Depression pathophysiology, risk prediction of recurrence and comorbid psychiatric disorders using genome-wide analyses.

Depression is a common psychiatric disorder and a leading cause of disability worldwide. Here we conducted a genome-wide association study meta-analysis of six datasets, including >1.3 million individuals (371,184 with depression) and identified 243 risk loci. Overall, 64 loci were new, including genes encoding glutamate and GABA receptors, which are targets for antidepressant drugs. Intersection with functional genomics data prioritized likely causal genes and revealed new enrichment of prenatal GABAergic neurons, astrocytes and oligodendrocyte lineages. We found depression to be highly polygenic, with ~11,700 variants explaining 90% of the single-nucleotide polymorphism heritability, estimating that >95% of risk variants for other psychiatric disorders (anxiety, schizophrenia, bipolar disorder and attention deficit hyperactivity disorder) were influencing depression risk when both concordant and discordant variants were considered, and nearly all depression risk variants influenced educational attainment. Additionally, depression genetic risk was associated with impaired complex cognition domains. We dissected the genetic and clinical heterogeneity, revealing distinct polygenic architectures across subgroups of depression and demonstrating significantly increased absolute risks for recurrence and psychiatric comorbidity among cases of depression with the highest polygenic burden, with considerable sex differences. The risks were up to 5- and 32-fold higher than cases with the lowest polygenic burden and the background population, respectively. These results deepen the understanding of the biology underlying depression, its disease progression and inform precision medicine approaches to treatment.

Genome-wide association study of school grades identifies genetic overlap between language ability, psychopathology and creativity.

Cognitive functions of individuals with psychiatric disorders differ from that of the general population. Such cognitive differences often manifest early in life as differential school performance and have a strong genetic basis. Here we measured genetic predictors of school performance in 30,982 individuals in English, Danish and mathematics via a genome-wide association study (GWAS) and studied their relationship with risk for six major psychiatric disorders. When decomposing the school performance into math and language-specific performances, we observed phenotypically and genetically a strong negative correlation between math performance and risk for most psychiatric disorders. But language performance correlated positively with risk for certain disorders, especially schizophrenia, which we replicate in an independent sample (n = 4547). We also found that the genetic variants relating to increased risk for schizophrenia and better language performance are overrepresented in individuals involved in creative professions (n = 2953) compared to the general population (n = 164,622). The findings together suggest that language ability, creativity and psychopathology might stem from overlapping genetic roots.

Genome-wide analyses of ADHD identify 27 risk loci, refine the genetic architecture and implicate several cognitive domains.

Attention-deficit hyperactivity disorder (ADHD) is a prevalent neurodevelopmental disorder with a major genetic component. Here, we present a genome-wide association study meta-analysis of ADHD comprising 38,691 individuals with ADHD and 186,843 controls. We identified 27 genome-wide significant loci, highlighting 76 potential risk genes enriched among genes expressed particularly in early brain development. Overall, ADHD genetic risk was associated with several brain-specific neuronal subtypes and midbrain dopaminergic neurons. In exome-sequencing data from 17,896 individuals, we identified an increased load of rare protein-truncating variants in ADHD for a set of risk genes enriched with probable causal common variants, potentially implicating SORCS3 in ADHD by both common and rare variants. Bivariate Gaussian mixture modeling estimated that 84-98% of ADHD-influencing variants are shared with other psychiatric disorders. In addition, common-variant ADHD risk was associated with impaired complex cognition such as verbal reasoning and a range of executive functions, including attention.

Polygenic liability, stressful life events and risk for secondary-treated depression in early life: a nationwide register-based case-cohort study.

BACKGROUND: In this study, we examined the relationship between polygenic liability for depression and number of stressful life events (SLEs) as risk factors for early-onset depression treated in inpatient, outpatient or emergency room settings at psychiatric hospitals in Denmark. METHODS: Data were drawn from the iPSYCH2012 case-cohort sample, a population-based sample of individuals born in Denmark between 1981 and 2005. The sample included 18 532 individuals who were diagnosed with depression by a psychiatrist by age 31 years, and a comparison group of 20 184 individuals. Information on SLEs was obtained from nationwide registers and operationalized as a time-varying count variable. Hazard ratios and cumulative incidence rates were estimated using Cox regressions. RESULTS: Risk for depression increased by 35% with each standard deviation increase in polygenic liability (p < 0.0001), and 36% (p < 0.0001) with each additional SLE. There was a small interaction between polygenic liability and SLEs (ß = -0.04, p = 0.0009). The probability of being diagnosed with depression in a hospital-based setting between ages 15 and 31 years ranged from 1.5% among males in the lowest quartile of polygenic liability with 0 events by age 15, to 18.8% among females in the highest quartile of polygenic liability with 4+ events by age 15. CONCLUSIONS: These findings suggest that although there is minimal interaction between polygenic liability and SLEs as risk factors for hospital-treated depression, combining information on these two important risk factors could potentially be useful for identifying high-risk individuals.

FarGen: Elucidating the distribution of coding variants in the isolated population of the Faroe Islands.

Here we present results from FarGen Phase I exomes. This dataset is based on the FarGen cohort, which consists of 1,541 individuals from the isolated population of the Faroe Islands. The purpose of this cohort is to serve as a reference catalog of coding variants, and to conduct population genetic studies to better understand the genetic contribution to various diseases in the Faroese population. The first whole-exome data set comprise 465 individuals and a total of 148,267 genetic variants were discovered. Principle Component Analysis indicates that the population is isolated and weakly structured. The distribution of variants in various functional classes was compared with populations in the gnomAD dataset; the results indicated that the proportions were consistent across the cohorts, but probably due to a small sample size, the FarGen dataset contained relatively few rare variants. We identified 19 variants that are classified as pathogenic or likely pathogenic in ClinVar; several of these variants are associated with monogenetic diseases with increased prevalence in the Faroe Islands. The results support previous studies, which indicate that the Faroe Islands is an isolated and weakly structured population. Future studies may elucidate the significance of the 19 pathogenic variants that were identified. The FarGen Phase I dataset is an important step for genetic research in the Faroese population, and the next phase of FarGen will increase the sample size and broaden the scope.

Genome-wide by Environment Interaction Study of Stressful Life Events and Hospital-Treated Depression in the iPSYCH2012 Sample.

Background: Researchers have long investigated a hypothesized interaction between genetic risk and stressful life events in the etiology of depression, but studies on the topic have yielded inconsistent results. Methods: We conducted a genome-wide by environment interaction study (GWEIS) in 18,532 patients with depression from hospital-based settings and 20,184 population controls. All individuals were drawn from the iPSYCH2012 case-cohort study, a nationally representative sample identified from Danish national registers. Information on stressful life events including family disruption, serious medical illness, death of a first-degree relative, parental disability, and child maltreatment was identified from the registers and operationalized as a time-varying count variable. Hazard ratios for main and interaction effects were estimated using Cox regressions weighted to accommodate the case-cohort design. Our replication sample included 22,880 depression cases and 50,378 controls from the UK Biobank. Results: The GWEIS in the iPSYCH2012 sample yielded three novel, genome-wide-significant (p < 5 × 10-8) loci located in the ABCC1 gene (rs56076205, p = 3.7 × 10-10), the AKAP6 gene (rs3784187, p = 1.2 × 10-8), and near the MFSD1 gene (rs340315, p = 4.5 × 10-8). No hits replicated in the UK Biobank (rs56076205: p = .87; rs3784187: p = .93; rs340315: p = .71). Conclusions: In this large, population-based GWEIS, we did not find any replicable hits for interaction. Future gene-by-stress research in depression should focus on establishing even larger collaborative GWEISs to attain sufficient power.

Identification of shared and differentiating genetic architecture for autism spectrum disorder, attention-deficit hyperactivity disorder and case subgroups.

Attention-deficit hyperactivity disorder (ADHD) and autism spectrum disorder (ASD) are highly heritable neurodevelopmental conditions, with considerable overlap in their genetic etiology. We dissected their shared and distinct genetic etiology by cross-disorder analyses of large datasets. We identified seven loci shared by the disorders and five loci differentiating them. All five differentiating loci showed opposite allelic directions in the two disorders and significant associations with other traits, including educational attainment, neuroticism and regional brain volume. Integration with brain transcriptome data enabled us to identify and prioritize several significantly associated genes. The shared genomic fraction contributing to both disorders was strongly correlated with other psychiatric phenotypes, whereas the differentiating portion was correlated most strongly with cognitive traits. Additional analyses revealed that individuals diagnosed with both ASD and ADHD were double-loaded with genetic predispositions for both disorders and showed distinctive patterns of genetic association with other traits compared with the ASD-only and ADHD-only subgroups. These results provide insights into the biological foundation of the development of one or both conditions and of the factors driving psychopathology discriminatively toward either ADHD or ASD.

Differences in the genetic architecture of common and rare variants in childhood, persistent and late-diagnosed attention-deficit hyperactivity disorder.

Attention-deficit hyperactivity disorder (ADHD) is a neurodevelopmental disorder with onset in childhood (childhood ADHD); two-thirds of affected individuals continue to have ADHD in adulthood (persistent ADHD), and sometimes ADHD is diagnosed in adulthood (late-diagnosed ADHD). We evaluated genetic differences among childhood (n = 14,878), persistent (n = 1,473) and late-diagnosed (n = 6,961) ADHD cases alongside 38,303 controls, and rare variant differences in 7,650 ADHD cases and 8,649 controls. We identified four genome-wide significant loci for childhood ADHD and one for late-diagnosed ADHD. We found increased polygenic scores for ADHD in persistent ADHD compared with the other two groups. Childhood ADHD had higher genetic overlap with hyperactivity and autism compared with late-diagnosed ADHD and the highest burden of rare protein-truncating variants in evolutionarily constrained genes. Late-diagnosed ADHD had a larger genetic overlap with depression than childhood ADHD and no increased burden in rare protein-truncating variants. Overall, these results suggest a genetic influence on age at first ADHD diagnosis, persistence of ADHD and the different comorbidity patterns among the groups.

Identifying the Common Genetic Basis of Antidepressant Response.

Background: Antidepressants are a first-line treatment for depression. However, only a third of individuals experience remission after the first treatment. Common genetic variation, in part, likely regulates antidepressant response, yet the success of previous genome-wide association studies has been limited by sample size. This study performs the largest genetic analysis of prospectively assessed antidepressant response in major depressive disorder to gain insight into the underlying biology and enable out-of-sample prediction. Methods: Genome-wide analysis of remission (n remit = 1852, n nonremit = 3299) and percentage improvement (n = 5218) was performed. Single nucleotide polymorphism-based heritability was estimated using genome-wide complex trait analysis. Genetic covariance with eight mental health phenotypes was estimated using polygenic scores/AVENGEME. Out-of-sample prediction of antidepressant response polygenic scores was assessed. Gene-level association analysis was performed using MAGMA and transcriptome-wide association study. Tissue, pathway, and drug binding enrichment were estimated using MAGMA. Results: Neither genome-wide association study identified genome-wide significant associations. Single nucleotide polymorphism-based heritability was significantly different from zero for remission (h 2 = 0.132, SE = 0.056) but not for percentage improvement (h 2 = -0.018, SE = 0.032). Better antidepressant response was negatively associated with genetic risk for schizophrenia and positively associated with genetic propensity for educational attainment. Leave-one-out validation of antidepressant response polygenic scores demonstrated significant evidence of out-of-sample prediction, though results varied in external cohorts. Gene-based analyses identified ETV4 and DHX8 as significantly associated with antidepressant response. Conclusions: This study demonstrates that antidepressant response is influenced by common genetic variation, has a genetic overlap schizophrenia and educational attainment, and provides a useful resource for future research. Larger sample sizes are required to attain the potential of genetics for understanding and predicting antidepressant response.

Mapping genomic loci implicates genes and synaptic biology in schizophrenia.

Schizophrenia has a heritability of 60-80%1, much of which is attributable to common risk alleles. Here, in a two-stage genome-wide association study of up to 76,755 individuals with schizophrenia and 243,649 control individuals, we report common variant associations at 287 distinct genomic loci. Associations were concentrated in genes that are expressed in excitatory and inhibitory neurons of the central nervous system, but not in other tissues or cell types. Using fine-mapping and functional genomic data, we identify 120 genes (106 protein-coding) that are likely to underpin associations at some of these loci, including 16 genes with credible causal non-synonymous or untranslated region variation. We also implicate fundamental processes related to neuronal function, including synaptic organization, differentiation and transmission. Fine-mapped candidates were enriched for genes associated with rare disruptive coding variants in people with schizophrenia, including the glutamate receptor subunit GRIN2A and transcription factor SP4, and were also enriched for genes implicated by such variants in neurodevelopmental disorders. We identify biological processes relevant to schizophrenia pathophysiology; show convergence of common and rare variant associations in schizophrenia and neurodevelopmental disorders; and provide a resource of prioritized genes and variants to advance mechanistic studies.

Atlantic herring (Clupea harengus) population structure in the Northeast Atlantic Ocean.

The Atlantic herring Clupea harengus L has a vast geographical distribution and a complex population structure with a few very large migratory units and many small local populations. Each population has its own spawning ground and/or time, thereby maintaining their genetic integrity. Several herring populations migrate between common feeding grounds and over-wintering areas resulting in frequent mixing of populations. Thus, many herring fisheries are based on mixed populations of different demographic status. In order to avoid over-exploitation of weak populations and to conserve biodiversity, understanding the population structure and population mixing is important for maintaining biologically sustainable herring fisheries. The aim of this study was to investigate the genetic population structure of herring in the Faroese and surrounding waters, and to develop genetic markers for distinguishing between four herring management units (often called stocks), namely the Norwegian spring-spawning herring (NSSH), Icelandic summer-spawning herring (ISSH), North Sea autumn-spawning herring (NSAH), and Faroese autumn-spawning herring (FASH). Herring from the four stocks were sequenced at low coverage, and single nucleotide polymorphisms (SNPs) were called and used for population structure analysis and individual assignment. An ancestry-informative SNP panel with 118 SNPs was developed and tested on 240 individuals. The results showed that all four stocks appeared to be genetically differentiated populations, but at lower levels of differentiation between FASH and ISSH than the other two populations. Overall assignment rate with the SNP panel was 80.7%, and agreement between the genetic and traditional visual assignment was 75.5%. The NSAH and NSSH samples had the highest assignment rate (100% and 98.3%, respectively) and highest agreement between traditional and genetic assignment methods (96.6% and 94.9%, respectively). The FASH and ISSH samples had substantially lower assignment rates (72.9% and 51.7%, respectively) and agreement between traditional and genetic methods (39.5% and 48.4%, respectively).

Genome-wide association study of more than 40,000 bipolar disorder cases provides new insights into the underlying biology.

Bipolar disorder is a heritable mental illness with complex etiology. We performed a genome-wide association study of 41,917 bipolar disorder cases and 371,549 controls of European ancestry, which identified 64 associated genomic loci. Bipolar disorder risk alleles were enriched in genes in synaptic signaling pathways and brain-expressed genes, particularly those with high specificity of expression in neurons of the prefrontal cortex and hippocampus. Significant signal enrichment was found in genes encoding targets of antipsychotics, calcium channel blockers, antiepileptics and anesthetics. Integrating expression quantitative trait locus data implicated 15 genes robustly linked to bipolar disorder via gene expression, encoding druggable targets such as HTR6, MCHR1, DCLK3 and FURIN. Analyses of bipolar disorder subtypes indicated high but imperfect genetic correlation between bipolar disorder type I and II and identified additional associated loci. Together, these results advance our understanding of the biological etiology of bipolar disorder, identify novel therapeutic leads and prioritize genes for functional follow-up studies.

Single nucleotide polymorphisms associated with methotrexate-induced nausea in juvenile idiopathic arthritis.

BACKGROUND: Context: Methotrexate (MTX) is a cornerstone in the treatment of juvenile idiopathic arthritis (JIA). MTX treatment is commonly associated with nausea. Large inter-individual variation exists in the level of MTX-induced nausea, possibly due to genetic factors. PURPOSE: To investigate whether MTX-induced nausea was associated with single nucleotide polymorphisms (SNPs) in genes encoding MTX-transporter proteins, a MTX metabolizing enzyme and a nausea receptor. FINDINGS: Methods: Children aged ≥9 years treated with MTX for JIA were eligible. MTX-induced nausea was registered by the children's completion of a nausea diary (min. 7 days) and the parents' completion of the MTX intolerance severity score (MISS). The selected SNPs were: SLCO1B1 (rs4149056; rs4149081), SLCO1B3 (rs2117032), SLC19A1 (rs1051266), ABCC2 (rs2273697; rs3740066; rs717620), ABCB1 (rs2032582; rs1045642), MTHFR (rs1801131, rs1801133), HTR3A (rs1062613; rs1985242; rs1176713) and HTR3B (rs1176744). RESULT: Enrolled were 121 JIA patients (82 girls: 39 boys) with a median age of 13.3 years (IQR: 11.3-15.1). The median MTX dose was 9.7 mg/m2/week (IQR: 9.0-10.9). The median MTX treatment duration prior to enrolment was 340 days (IQR: 142-766). The SNP analysis was available for 119 patients. MTX intolerance was associated with the genotype distribution of rs1801133 (MTHFR) (p = 0.02). There was no additive effect of the minor alleles for any of the selected SNPs, nor any significant haplotype associations. CONCLUSION: Summary: MTX-induced nausea may be influenced by genetic polymorphisms in a MTX metabolizing enzyme (rs1801133; MTHFR). IMPLICATIONS: Further analyses involving inclusion of larger cohorts are needed to understand the impact of SNPs on MTX-induced nausea in JIA.

Genetic analyses identify widespread sex-differential participation bias.

Genetic association results are often interpreted with the assumption that study participation does not affect downstream analyses. Understanding the genetic basis of participation bias is challenging since it requires the genotypes of unseen individuals. Here we demonstrate that it is possible to estimate comparative biases by performing a genome-wide association study contrasting one subgroup versus another. For example, we showed that sex exhibits artifactual autosomal heritability in the presence of sex-differential participation bias. By performing a genome-wide association study of sex in approximately 3.3 million males and females, we identified over 158 autosomal loci spuriously associated with sex and highlighted complex traits underpinning differences in study participation between the sexes. For example, the body mass index-increasing allele at FTO was observed at higher frequency in males compared to females (odds ratio = 1.02, P = 4.4 × 10-36). Finally, we demonstrated how these biases can potentially lead to incorrect inferences in downstream analyses and propose a conceptual framework for addressing such biases. Our findings highlight a new challenge that genetic studies may face as sample sizes continue to grow.

Risk variants and polygenic architecture of disruptive behavior disorders in the context of attention-deficit/hyperactivity disorder.

Attention-Deficit/Hyperactivity Disorder (ADHD) is a childhood psychiatric disorder often comorbid with disruptive behavior disorders (DBDs). Here, we report a GWAS meta-analysis of ADHD comorbid with DBDs (ADHD + DBDs) including 3802 cases and 31,305 controls. We identify three genome-wide significant loci on chromosomes 1, 7, and 11. A meta-analysis including a Chinese cohort supports that the locus on chromosome 11 is a strong risk locus for ADHD + DBDs across European and Chinese ancestries (rs7118422, P = 3.15×10-10, OR = 1.17). We find a higher SNP heritability for ADHD + DBDs (h2SNP = 0.34) when compared to ADHD without DBDs (h2SNP = 0.20), high genetic correlations between ADHD + DBDs and aggressive (rg = 0.81) and anti-social behaviors (rg = 0.82), and an increased burden (polygenic score) of variants associated with ADHD and aggression in ADHD + DBDs compared to ADHD without DBDs. Our results suggest an increased load of common risk variants in ADHD + DBDs compared to ADHD without DBDs, which in part can be explained by variants associated with aggressive behavior.