Genome-wide association study identifies 30 obsessive-compulsive disorder associated loci

Obsessive-compulsive disorder (OCD) affects ~1% of the population and exhibits a high SNP-heritability, yet previous genome-wide association studies (GWAS) have provided limited information on the genetic etiology and underlying biological mechanisms of the disorder. We conducted a GWAS meta-analysis combining 53,660 OCD cases and 2,044,417 controls from 28 European-ancestry cohorts revealing 30 independent genome-wide significant SNPs and a SNP-based heritability of 6.7%. Separate GWAS for clinical, biobank, comorbid, and self-report sub-groups found no evidence of sample ascertainment impacting our results. Functional and positional QTL gene-based approaches identified 249 significant candidate risk genes for OCD, of which 25 were identified as putatively causal, highlighting WDR6, DALRD3, CTNND1 and genes in the MHC region. Tissue and single-cell enrichment analyses highlighted hippocampal and cortical excitatory neurons, along with D1- and D2-type dopamine receptor-containing medium spiny neurons, as playing a role in OCD risk. OCD displayed significant genetic correlations with 65 out of 112 examined phenotypes. Notably, it showed positive genetic correlations with all included psychiatric phenotypes, in particular anxiety, depression, anorexia nervosa, and Tourette syndrome, and negative correlations with a subset of the included autoimmune disorders, educational attainment, and body mass index.. This study marks a significant step toward unraveling its genetic landscape and advances understanding of OCD genetics, providing a foundation for future interventions to address this debilitating disorder.

Genome-wide association study identifies new loci associated with OCD

To date, four genome-wide association studies (GWAS) of obsessive-compulsive disorder (OCD) have been published, reporting a high single-nucleotide polymorphism (SNP)-heritability of 28% but finding only one significant SNP. A substantial increase in sample size will likely lead to further identification of SNPs, genes, and biological pathways mediating the susceptibility to OCD. We conducted a GWAS meta-analysis with a 2-3-fold increase in case sample size (OCD cases: N = 37,015, controls: N = 948,616) compared to the last OCD GWAS, including six previously published cohorts (OCGAS, IOCDF-GC, IOCDF-GC-trio, NORDiC-nor, NORDiC-swe, and iPSYCH) and unpublished self-report data from 23andMe Inc. We explored the genetic architecture of OCD by conducting gene-based tests, tissue and celltype enrichment analyses, and estimating heritability and genetic correlations with 74 phenotypes. To examine a potential heterogeneity in our data, we conducted multivariable GWASs with MTAG. We found support for 15 independent genome-wide significant loci (14 new) and 79 protein-coding genes. Tissue enrichment analyses implicate multiple cortical regions, the amygdala, and hypothalamus, while cell type analyses yielded 12 cell types linked to OCD (all neurons). The SNP-based heritability of OCD was estimated to be 0.08. Using MTAG we found evidence for specific genetic underpinnings characteristic of different cohort-ascertainment and identified additional significant SNPs. OCD was genetically correlated with 40 disorders or traits-positively with all psychiatric disorders and negatively with BMI, age at first birth and multiple autoimmune diseases. The GWAS meta-analysis identified several biologically informative genes as important contributors to the aetiology of OCD. Overall, we have begun laying the groundwork through which the biology of OCD will be understood and described.

Polygenic profiles define aspects of clinical heterogeneity in attention deficit hyperactivity disorder

Attention deficit hyperactivity disorder (ADHD) is a complex disorder that manifests variability in long-term outcomes and clinical presentations. The genetic contributions to such heterogeneity are not well understood. Here we show several genetic links to clinical heterogeneity in ADHD in a case-only study of 14,084 diagnosed individuals. First, we identify one genome-wide significant locus by comparing cases with ADHD and autism spectrum disorder (ASD) to cases with ADHD but not ASD. Second, we show that cases with ASD and ADHD, substance use disorder and ADHD, or first diagnosed with ADHD in adulthood have unique polygenic score (PGS) profiles that distinguish them from complementary case subgroups and controls. Finally, a PGS for an ASD diagnosis in ADHD cases predicted cognitive performance in an independent developmental cohort. Our approach uncovered evidence of genetic heterogeneity in ADHD, helping us to understand its etiology and providing a model for studies of other disorders.

Smoking during pregnancy and risk of multiple sclerosis in offspring and mother: A Danish nationwide register-based cohort study

BACKGROUND

The association between intra-uterine exposure to maternal smoking and risk of multiple sclerosis (MS) has been little studied and with conflicting results.

OBJECTIVE

To examine the risk of MS in offspring exposed intra-uterine to maternal smoking. In addition, to re-examine prior observations of an elevated risk of MS among smokers, assuming that self-reported smoking during pregnancy reflects the woman's general smoking habits.

METHODS

The study cohort included all Danish women, pregnant in the period 1991-2018, (n = 789,299) and singletons from these pregnancies (n = 879,135). Nationwide information on maternal smoking during pregnancy and MS cases in the study cohort were obtained from the Medical Birth Register and the National Patient Register. Cox regression analysis was used to estimate hazard ratios (HRs) for the association between smoking and MS risk.

RESULTS

Women who smoked during pregnancy had a 42% increased risk of developing MS compared with non-smoking women (HR = 1.42 (1.32-1.52), n = 1,296). The risk of MS among singletons of women who smoked during pregnancy was 38% higher than that among singletons born to non-smoking women (HR = 1.38 (1.08-1.76), n = 110).

CONCLUSION

Our observations add further to the evidence implicating smoking in the development of MS and suggest that intra-uterine exposure to tobacco smoke may increase MS risk.

Metabolic signature of the pathogenic 22q11.2 deletion identifies carriers and provides insight into systemic dysregulation

Large deletions at chromosome 22q11.2 are known to cause severe clinical conditions collectively known as 22q11.2 deletion syndrome. Notwithstanding the pathogenicity of these deletions, affected individuals are typically diagnosed in late childhood or early adolescence, and little is known of the molecular signaling cascades and biological consequences immediately downstream of the deleted genes. Here, we used targeted metabolomics to compare neonatal dried blood spot samples from 203 individuals clinically identified as carriers of a deletion at chromosome 22q11.2 with 203 unaffected individuals. A total of 173 metabolites were successfully identified and used to inform on systemic dysregulation caused by the genomic lesion and to discriminate carriers from non-carriers. We found 84 metabolites to be differentially abundant between carriers and non-carriers of the 22q11.2 deletion. A predictive model based on all 173 metabolites achieved high Accuracy (89%), Area Under the Curve (93%), F1 (88%), Positive Predictive Value (94%), and Negative Predictive Value (84%) with tyrosine and proline having the highest individual contributions to the model as well as the highest interaction strength. Targeted metabolomics provides insight into the molecular consequences possibly contributing to the pathology underlying the clinical manifestations of the 22q11 deletion and is an easily applicable approach to first-pass screening for carrier status of the 22q11 to prompt subsequent verification of the genomic diagnosis.

The correlates of neonatal complement component 3 and 4 protein concentrations with a focus on psychiatric and autoimmune disorders

Complement components have been linked to schizophrenia and autoimmune disorders. We examined the association between neonatal circulating C3 and C4 protein concentrations in 68,768 neonates and the risk of six mental disorders. We completed genome-wide association studies (GWASs) for C3 and C4 and applied the summary statistics in Mendelian randomization and phenome-wide association studies related to mental and autoimmune disorders. The GWASs for C3 and C4 protein concentrations identified 15 and 36 independent loci, respectively. We found no associations between neonatal C3 and C4 concentrations and mental disorders in the total sample (both sexes combined); however, post-hoc analyses found that a higher C3 concentration was associated with a reduced risk of schizophrenia in females. Mendelian randomization based on C4 summary statistics found an altered risk of five types of autoimmune disorders. Our study adds to our understanding of the associations between C3 and C4 concentrations and subsequent mental and autoimmune disorders.

Multi-ancestry genome-wide association study of cannabis use disorder yields insight into disease biology and public health implications

As recreational use of cannabis is being decriminalized in many places and medical use widely sanctioned, there are growing concerns about increases in cannabis use disorder (CanUD), which is associated with numerous medical comorbidities. Here we performed a genome-wide association study of CanUD in the Million Veteran Program (MVP), followed by meta-analysis in 1,054,365 individuals (ncases = 64,314) from four broad ancestries designated by the reference panel used for assignment (European n = 886,025, African n = 123,208, admixed American n = 38,289 and East Asian n = 6,843). Population-specific methods were applied to calculate single nucleotide polymorphism-based heritability within each ancestry. Statistically significant single nucleotide polymorphism-based heritability for CanUD was observed in all but the smallest population (East Asian). We discovered genome-wide significant loci unique to each ancestry: 22 in European, 2 each in African and East Asian, and 1 in admixed American ancestries. A genetically informed causal relationship analysis indicated a possible effect of genetic liability for CanUD on lung cancer risk, suggesting potential unanticipated future medical and psychiatric public health consequences that require further study to disentangle from other known risk factors such as cigarette smoking.

Multi-ancestry study of the genetics of problematic alcohol use in over 1 million individuals

Problematic alcohol use (PAU), a trait that combines alcohol use disorder and alcohol-related problems assessed with a questionnaire, is a leading cause of death and morbidity worldwide. Here we conducted a large cross-ancestry meta-analysis of PAU in 1,079,947 individuals (European, N = 903,147; African, N = 122,571; Latin American, N = 38,962; East Asian, N = 13,551; and South Asian, N = 1,716 ancestries). We observed a high degree of cross-ancestral similarity in the genetic architecture of PAU and identified 110 independent risk variants in within- and cross-ancestry analyses. Cross-ancestry fine mapping improved the identification of likely causal variants. Prioritizing genes through gene expression and chromatin interaction in brain tissues identified multiple genes associated with PAU. We identified existing medications for potential pharmacological studies by a computational drug repurposing analysis. Cross-ancestry polygenic risk scores showed better performance of association in independent samples than single-ancestry polygenic risk scores. Genetic correlations between PAU and other traits were observed in multiple ancestries, with other substance use traits having the highest correlations. This study advances our knowledge of the genetic etiology of PAU, and these findings may bring possible clinical applicability of genetics insights-together with neuroscience, biology and data science-closer.

Postpartum and non-postpartum depression: a population-based matched case-control study comparing polygenic risk scores for severe mental disorders

It remains inconclusive whether postpartum depression (PPD) and depression with onset outside the postpartum period (MDD) are genetically distinct disorders. We aimed to investigate whether polygenic risk scores (PGSs) for major mental disorders differ between PPD cases and MDD cases in a nested case-control study of 50,057 women born from 1981 to 1997 in the iPSYCH2015 sample in Demark. We identified 333 women with first-onset postpartum depression (PPD group), who were matched with 993 women with first-onset depression diagnosed outside of postpartum (MDD group), and 999 female population controls. Data on genetics and depressive disorders were retrieved from neonatal biobanks and the Psychiatric Central Research Register. PGSs were calculated from both individual-level genetic data and meta-analysis summary statistics from the Psychiatric Genomics Consortium. Conditional logistic regression was used to calculate the odds ratio (OR), accounting for the selection-related reproductive behavior. After adjustment for covariates, higher PGSs for severe mental disorders were associated with increased ORs of both PPD and MDD. Compared with MDD cases, MDD PGS and attention-deficit/hyperactivity disorder PGS were marginally but not statistically higher for PPD cases, with the OR of PPD versus MDD being 1.12 (95% CI: 0 .97-1.29) and 1.11 (0.97-1.27) per-standard deviation increase, respectively. The ORs of PPD versus MDD did not statistically differ by PGSs of bipolar disorder, schizophrenia, or autism spectrum disorder. Our findings suggest that relying on PGS data, there was no clear evidence of distinct genetic make-up of women with depression occurring during or outside postpartum, after taking the selection-related reproductive behavior into account.

Genome-wide association study of placental weight identifies distinct and shared genetic influences between placental and fetal growth

A well-functioning placenta is essential for fetal and maternal health throughout pregnancy. Using placental weight as a proxy for placental growth, we report genome-wide association analyses in the fetal (n = 65,405), maternal (n = 61,228) and paternal (n = 52,392) genomes, yielding 40 independent association signals. Twenty-six signals are classified as fetal, four maternal and three fetal and maternal. A maternal parent-of-origin effect is seen near KCNQ1. Genetic correlation and colocalization analyses reveal overlap with birth weight genetics, but 12 loci are classified as predominantly or only affecting placental weight, with connections to placental development and morphology, and transport of antibodies and amino acids. Mendelian randomization analyses indicate that fetal genetically mediated higher placental weight is causally associated with preeclampsia risk and shorter gestational duration. Moreover, these analyses support the role of fetal insulin in regulating placental weight, providing a key link between fetal and placental growth.

Identification of women at high risk of postpartum psychiatric episodes: A population-based study quantifying relative and absolute risks following exposure to selected risk factors and genetic liability

BACKGROUND

We quantified relative and absolute risks of postpartum psychiatric episodes (PPE) following risk factors: Young age, past personal or family history of psychiatric disorders, and genetic liability.

METHODS

We conducted a register-based study using the iPSYCH2012 case-cohort sample. Exposures were personal history of psychiatric episodes prior to childbirth, being a young mother (giving birth before the age of 21.5 years), having a family history of psychiatric disorders, and a high (highest quartile) polygenic score (PGS) for major depression. PPE was defined within 12 months postpartum by prescription of psychotropic medication or in- and outpatient contact to a psychiatric facility. We included primiparous women born 1981-1999, giving birth before January 1st, 2016. We conducted Cox regression to calculate hazard ratios (HRs) of PPE, absolute risks were calculated using cumulative incidence functions.

RESULTS

We included 8174 primiparous women, and the estimated baseline PPE risk was 6.9% (95% CI 6.0%-7.8%, number of PPE cases: 2169). For young mothers with a personal and family history of psychiatric disorders, the absolute risk of PPE was 21.6% (95% CI 15.9%-27.8%). Adding information on high genetic liability to depression, the risk increased to 29.2% (95% CI 21.3%-38.4%) for PPE.

CONCLUSIONS

Information on prior personal and family psychiatric episodes as well as age may assist in estimating a personalized risk of PPE. Furthermore, additional information on genetic liability could add even further to this risk assessment.

GWAS Meta-Analysis of Suicide Attempt: Identification of 12 Genome-Wide Significant Loci and Implication of Genetic Risks for Specific Health Factors

OBJECTIVE

Suicidal behavior is heritable and is a major cause of death worldwide. Two large-scale genome-wide association studies (GWASs) recently discovered and cross-validated genome-wide significant (GWS) loci for suicide attempt (SA). The present study leveraged the genetic cohorts from both studies to conduct the largest GWAS meta-analysis of SA to date. Multi-ancestry and admixture-specific meta-analyses were conducted within groups of significant African, East Asian, and European ancestry admixtures.

METHODS

This study comprised 22 cohorts, including 43,871 SA cases and 915,025 ancestry-matched controls. Analytical methods across multi-ancestry and individual ancestry admixtures included inverse variance-weighted fixed-effects meta-analyses, followed by gene, gene-set, tissue-set, and drug-target enrichment, as well as summary-data-based Mendelian randomization with brain expression quantitative trait loci data, phenome-wide genetic correlation, and genetic causal proportion analyses.

RESULTS

Multi-ancestry and European ancestry admixture GWAS meta-analyses identified 12 risk loci at p values <5×10-8. These loci were mostly intergenic and implicated DRD2, SLC6A9, FURIN, NLGN1, SOX5, PDE4B, and CACNG2. The multi-ancestry SNP-based heritability estimate of SA was 5.7% on the liability scale (SE=0.003, p=5.7×10-80). Significant brain tissue gene expression and drug set enrichment were observed. There was shared genetic variation of SA with attention deficit hyperactivity disorder, smoking, and risk tolerance after conditioning SA on both major depressive disorder and posttraumatic stress disorder. Genetic causal proportion analyses implicated shared genetic risk for specific health factors.

CONCLUSIONS

This multi-ancestry analysis of suicide attempt identified several loci contributing to risk and establishes significant shared genetic covariation with clinical phenotypes. These findings provide insight into genetic factors associated with suicide attempt across ancestry admixture populations, in veteran and civilian populations, and in attempt versus death.

ADuLT: An efficient and robust time-to-event GWAS

Proportional hazards models have been proposed to analyse time-to-event phenotypes in genome-wide association studies (GWAS). However, little is known about the ability of proportional hazards models to identify genetic associations under different generative models and when ascertainment is present. Here we propose the age-dependent liability threshold (ADuLT) model as an alternative to a Cox regression based GWAS, here represented by SPACox. We compare ADuLT, SPACox, and standard case-control GWAS in simulations under two generative models and with varying degrees of ascertainment as well as in the iPSYCH cohort. We find Cox regression GWAS to be underpowered when cases are strongly ascertained (cases are oversampled by a factor 5), regardless of the generative model used. ADuLT is robust to ascertainment in all simulated scenarios. Then, we analyse four psychiatric disorders in iPSYCH, ADHD, Autism, Depression, and Schizophrenia, with a strong case-ascertainment. Across these psychiatric disorders, ADuLT identifies 20 independent genome-wide significant associations, case-control GWAS finds 17, and SPACox finds 8, which is consistent with simulation results. As more genetic data are being linked to electronic health records, robust GWAS methods that can make use of age-of-onset information will help increase power in analyses for common health outcomes.

Population-based Risk of Psychiatric Disorders Associated with Recurrent CNVs

Recurrent copy number variants (rCNVs) are associated with increased risk of neuropsychiatric disorders but their pathogenic population-level impact is unknown. We provide population-based estimates of rCNV-associated risk of neuropsychiatric disorders for 34 rCNVs in the iPSYCH2015 case-cohort sample (n=120,247). Most observed significant increases in rCNV-associated risk for ADHD, autism or schizophrenia were moderate (HR:1.42-5.00), and risk estimates were highly correlated across these disorders, the most notable exception being high autism-associated risk with Prader-Willi/Angelman Syndrome duplications (HR=20.8). No rCNV was associated with significant increase in depression risk. Also, rCNV-associated risk was positively correlated with locus size and gene constraint. Comparison with published rCNV studies suggests that prevalence of some rCNVs is higher, and risk of psychiatric disorders lower, than previously estimated. In an era where genetics is increasingly being clinically applied, our results highlight the importance of population-based risk estimates for genetics-based predictions.

Multi-PGS enhances polygenic prediction by combining 937 polygenic scores

The predictive performance of polygenic scores (PGS) is largely dependent on the number of samples available to train the PGS. Increasing the sample size for a specific phenotype is expensive and takes time, but this sample size can be effectively increased by using genetically correlated phenotypes. We propose a framework to generate multi-PGS from thousands of publicly available genome-wide association studies (GWAS) with no need to individually select the most relevant ones. In this study, the multi-PGS framework increases prediction accuracy over single PGS for all included psychiatric disorders and other available outcomes, with prediction R2 increases of up to 9-fold for attention-deficit/hyperactivity disorder compared to a single PGS. We also generate multi-PGS for phenotypes without an existing GWAS and for case-case predictions. We benchmark the multi-PGS framework against other methods and highlight its potential application to new emerging biobanks.

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.

Effects of Long-Term Storage on the Biobanked Neonatal Dried Blood Spot Metabolome

Over 2.5 million neonatal dried blood spots (DBS) are stored at the Danish National Biobank. These samples offer extraordinary possibilities for metabolomics research, including prediction of disease and understanding of underlying molecular mechanisms of disease development. Nevertheless, Danish neonatal DBS have been little explored in metabolomics studies. One question that remains underinvestigated is the long-term stability of the large number of metabolites typically assessed in untargeted metabolomics over long time periods of storage. Here, we investigate temporal trends of metabolites measured in 200 neonatal DBS collected over a time course of 10 years, using an untargeted liquid chromatography tandem mass spectrometry (LC-MS/MS) based metabolomics protocol. We found that a majority (71%) of the metabolome was stable during 10 years of storage at -20 °C. However, we found decreasing trends for lipid-related metabolites, such as glycerophosphocholines and acylcarnitines. A few metabolites, including glutathione and methionine, may be strongly influenced by storage, with changes in metabolite levels up to 0.1-0.2 standard deviation units per year. Our findings indicate that untargeted metabolomics of DBS samples, with long-term storage in biobanks, is suitable for retrospective epidemiological studies. We identify metabolites whose stability in DBS should be closely monitored in future studies of DBS samples with long-term storage.

A Review of Major Danish Biobanks: Advantages and Possibilities of Health Research in Denmark

Biobank research may lead to an improved understanding of disease etiology and advance personalized medicine. Denmark (population ~5.9 million) provides a unique setting for population-based health research. The country is a rich source of biobanks and the universal, tax-funded healthcare system delivers routinely collected data to numerous registries and databases. By virtue of the civil registration number (assigned uniquely to all Danish citizens), biological specimens stored in biobanks can be combined with clinical and demographic data from these population-based health registries and databases. In this review, we aim to provide an understanding of advantages and possibilities of biobank research in Denmark. As knowledge about the Danish setting is needed to grasp the full potential, we first introduce the Danish healthcare system, the Civil Registration System, the population-based registries, and the interface with biobanks. We then describe the biobank infrastructures, comprising the Danish National Biobank Initiative, the Bio- and Genome Bank Denmark, and the Danish National Genome Center. Further, we briefly provide an overview of fourteen selected biobanks, including: The Danish Newborn Screening Biobank; The Danish National Birth Cohort; The Danish Twin Registry Biobank; Diet, Cancer and Health; Diet, Cancer and Health - Next generations; Danish Centre for Strategic Research in Type 2 Diabetes; Vejle Diabetes Biobank; The Copenhagen Hospital Biobank; The Copenhagen City Heart Study; The Copenhagen General Population Study; The Danish Cancer Biobank; The Danish Rheumatological Biobank; The Danish Blood Donor Study; and The Danish Pathology Databank. Last, we inform on practical aspects, such as data access, and discuss future implications.

Genetic correlates of vitamin D-binding protein and 25-hydroxyvitamin D in neonatal dried blood spots

The vitamin D binding protein (DBP), encoded by the group-specific component (GC) gene, is a component of the vitamin D system. In a genome-wide association study of DBP concentration in 65,589 neonates we identify 26 independent loci, 17 of which are in or close to the GC gene, with fine-mapping identifying 2 missense variants on chromosomes 12 and 17 (within SH2B3 and GSDMA, respectively). When adjusted for GC haplotypes, we find 15 independent loci distributed over 10 chromosomes. Mendelian randomization analyses identify a unidirectional effect of higher DBP concentration and (a) higher 25-hydroxyvitamin D concentration, and (b) a reduced risk of multiple sclerosis and rheumatoid arthritis. A phenome-wide association study confirms that higher DBP concentration is associated with a reduced risk of vitamin D deficiency. Our findings provide valuable insights into the influence of DBP on vitamin D status and a range of health outcomes.

Associations of psychiatric disorders with sex chromosome aneuploidies in the Danish iPSYCH2015 dataset: a case-cohort study

BACKGROUND

Increased prevalence of mental illness has been reported in clinical studies of sex chromosome aneuploidies, but accurate population-based estimates of the prevalence and clinical detection rate of sex chromosome aneuploidies and the associated risks of psychiatric disorders are needed. In this study, we provide such estimates, valid for children and young adults of the contemporary Danish population.

METHODS

We used the iPSYCH2015 case-cohort dataset, which is based on a source population of single-born individuals born in Denmark between May 1, 1981, and Dec 31, 2008. The case sample comprises all individuals from the source population with a diagnosis of any index psychiatric disorder (schizophrenia spectrum disorder, bipolar disorder, major depressive disorder, autism spectrum disorder, or ADHD) by the end of follow-up (Dec 31, 2015), registered in the hospital-based Danish Psychiatric Central Research Register. The cohort consists of individuals randomly selected from the source population, and overlaps with the case sample. Biobanked blood samples for individuals in the case and cohort samples underwent genotyping and quality-control filtering, after which we analysed microarray data to detect sex chromosome aneuploidy karyotypes (45,X, 47,XXX, 47,XXY, and 47,XYY). We estimated the population-valid prevalence of these karyotypes from the cohort sample. Weighted Cox proportional hazards models were used to estimate the risks of each index psychiatric disorder associated with each sex chromosome aneuploidy karyotype, by use of date of first hospitalisation with the index disorder in the respective case group and the cohort as outcome. The clinical detection rate was determined by comparing records of clinical diagnoses of genetic conditions from the Danish National Patient Register with sex chromosome aneuploidy karyotype determined by our study.

FINDINGS

The assessed sample comprised 119 481 individuals (78 726 in the case sample and 43 326 in the cohort) who had genotyped and quality-control-filtered blood samples, including 64 533 (54%) people of gonadal male sex and 54 948 (46%) of gonadal female sex. Age during follow-up ranged from 0 to 34·7 years (mean 10·9 years [SD 3·5 years]). Information on ethnicity was not available. We identified 387 (0·3%) individuals as carriers of sex chromosome aneuploidies. The overall prevalence of sex chromosome aneuploidies was 1·5 per 1000 individuals. Each sex chromosome aneuploidy karyotype was associated with an increased risk of at least one index psychiatric disorder, with hazard ratios (HRs) of 2·20 (95% CI 1·42-3·39) for 47,XXY; 2·73 (1·25-6·00) for 47,XXX; 3·56 (1·01-12·53) for 45,X; and 4·30 (2·48-7·55) for 47,XYY. All karyotypes were associated with an increased risk of ADHD (HRs ranging from 1·99 [1·24-3·19] to 6·15 [1·63-23·19]), autism spectrum disorder (2·72 [1·72-4·32] to 8·45 [2·49-28·61]), and schizophrenia spectrum disorder (1·80 [1·15-2·80] to 4·60 [1·57-13·51]). Increased risk of major depressive disorder was found for individuals with 47,XXY (1·88 [1·07-3·33]) and 47,XYY (2·65 [1·12-5·90]), and of bipolar disorder for those with 47,XXX (4·32 [1·12-16·62]). The proportion of sex chromosome aneuploidy carriers who had been clinically diagnosed was 93% for 45,X, but lower for 47,XXY (22%), 47,XXX (15%), and 47,XYY (15%). Among carriers, the risk of diagnosis of at least one index psychiatric disorder did not significantly differ between those who had and had not been clinically diagnosed with sex chromosome aneuploidies (p=0·65).

INTERPRETATION

Increased risks of psychiatric disorders associated with sex chromosome aneuploidies, combined with low rates of clinical diagnosis of sex chromosome aneuploidies, compromise the adequate provision of necessary health care and counselling to affected individuals and their families, which might be helped by increased application of genetic testing in clinical settings.

FUNDING

Lundbeck Foundation and National Institutes of Health.

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.

Deep Learning for Cross-Diagnostic Prediction of Mental Disorder Diagnosis and Prognosis Using Danish Nationwide Register and Genetic Data

Importance

Diagnoses and treatment of mental disorders are hampered by the current lack of objective markers needed to provide a more precise diagnosis and treatment strategy.

Objective

To develop deep learning models to predict mental disorder diagnosis and severity spanning multiple diagnoses using nationwide register data, family and patient-specific diagnostic history, birth-related measurement, and genetics.

Design, Setting, and Participants

This study was conducted from May 1, 1981, to December 31, 2016. For the analysis, which used a Danish population-based case-cohort sample of individuals born between 1981 and 2005, genotype data and matched longitudinal health register data were taken from the longitudinal Danish population-based Integrative Psychiatric Research Consortium 2012 case-cohort study. Included were individuals with mental disorders (attention-deficit/hyperactivity disorder [ADHD]), autism spectrum disorder (ASD), major depressive disorder (MDD), bipolar disorder (BD), schizophrenia spectrum disorders (SCZ), and population controls. Data were analyzed from February 1, 2021, to January 24, 2022.

Exposure

At least 1 hospital contact with diagnosis of ADHD, ASD, MDD, BD, or SCZ.

Main Outcomes and Measures

The predictability of (1) mental disorder diagnosis and (2) severity trajectories (measured by future outpatient hospital contacts, admissions, and suicide attempts) were investigated using both a cross-diagnostic and single-disorder setup. Predictive power was measured by AUC, accuracy, and Matthews correlation coefficient (MCC), including an estimate of feature importance.

Results

A total of 63 535 individuals (mean [SD] age, 23 [7] years; 34 944 male [55%]; 28 591 female [45%]) were included in the model. Based on data prior to diagnosis, the specific diagnosis was predicted in a multidiagnostic prediction model including the background population with an overall area under the curve (AUC) of 0.81 and MCC of 0.28, whereas the single-disorder models gave AUCs/MCCs of 0.84/0.54 for SCZ, 0.79/0.41 for BD, 0.77/0.39 for ASD, 0.74/0.38, for ADHD, and 0.74/0.38 for MDD. The most important data sets for multidiagnostic prediction were previous mental disorders and age (11%-23% reduction in prediction accuracy when removed) followed by family diagnoses, birth-related measurements, and genetic data (3%-5% reduction in prediction accuracy when removed). Furthermore, when predicting subsequent disease trajectories of the disorder, the most severe cases were the most easily predictable, with an AUC of 0.72.

Conclusions and Relevance

Results of this diagnostic study suggest the possibility of combining genetics and registry data to predict both mental disorder diagnosis and disorder progression in a clinically relevant, cross-diagnostic setting prior to clinical assessment.

Identification of Novel, Replicable Genetic Risk Loci for Suicidal Thoughts and Behaviors Among US Military Veterans

Importance

Suicide is a leading cause of death; however, the molecular genetic basis of suicidal thoughts and behaviors (SITB) remains unknown.

Objective

To identify novel, replicable genomic risk loci for SITB.

Design, Setting, and Participants

This genome-wide association study included 633 778 US military veterans with and without SITB, as identified through electronic health records. GWAS was performed separately by ancestry, controlling for sex, age, and genetic substructure. Cross-ancestry risk loci were identified through meta-analysis. Study enrollment began in 2011 and is ongoing. Data were analyzed from November 2021 to August 2022.

Main Outcome and Measures

SITB.

Results

A total of 633 778 US military veterans were included in the analysis (57 152 [9%] female; 121 118 [19.1%] African ancestry, 8285 [1.3%] Asian ancestry, 452 767 [71.4%] European ancestry, and 51 608 [8.1%] Hispanic ancestry), including 121 211 individuals with SITB (19.1%). Meta-analysis identified more than 200 GWS (P < 5 × 10-8) cross-ancestry risk single-nucleotide variants for SITB concentrated in 7 regions on chromosomes 2, 6, 9, 11, 14, 16, and 18. Top single-nucleotide variants were largely intronic in nature; 5 were independently replicated in ISGC, including rs6557168 in ESR1, rs12808482 in DRD2, rs77641763 in EXD3, rs10671545 in DCC, and rs36006172 in TRAF3. Associations for FBXL19 and AC018880.2 were not replicated. Gene-based analyses implicated 24 additional GWS cross-ancestry risk genes, including FURIN, TSNARE1, and the NCAM1-TTC12-ANKK1-DRD2 gene cluster. Cross-ancestry enrichment analyses revealed significant enrichment for expression in brain and pituitary tissue, synapse and ubiquitination processes, amphetamine addiction, parathyroid hormone synthesis, axon guidance, and dopaminergic pathways. Seven other unique European ancestry-specific GWS loci were identified, 2 of which (POM121L2 and METTL15/LINC02758) were replicated. Two additional GWS ancestry-specific loci were identified within the African ancestry (PET112/GATB) and Hispanic ancestry (intergenic locus on chromosome 4) subsets, both of which were replicated. No GWS loci were identified within the Asian ancestry subset; however, significant enrichment was observed for axon guidance, cyclic adenosine monophosphate signaling, focal adhesion, glutamatergic synapse, and oxytocin signaling pathways across all ancestries. Within the European ancestry subset, genetic correlations (r > 0.75) were observed between the SITB phenotype and a suicide attempt-only phenotype, depression, and posttraumatic stress disorder. Additionally, polygenic risk score analyses revealed that the Million Veteran Program polygenic risk score had nominally significant main effects in 2 independent samples of veterans of European and African ancestry.

Conclusions and Relevance

The findings of this analysis may advance understanding of the molecular genetic basis of SITB and provide evidence for ESR1, DRD2, TRAF3, and DCC as cross-ancestry candidate risk genes. More work is needed to replicate these findings and to determine if and how these genes might impact clinical care.

Multi-ancestry study of the genetics of problematic alcohol use in >1 million individuals

Problematic alcohol use (PAU) is a leading cause of death and disability worldwide. To improve our understanding of the genetics of PAU, we conducted a large cross-ancestry meta-analysis of PAU in 1,079,947 individuals. We observed a high degree of cross-ancestral similarity in the genetic architecture of PAU and identified 110 independent risk variants in within- and cross-ancestry analyses. Cross-ancestry fine-mapping improved the identification of likely causal variants. Prioritizing genes through gene expression and/or chromatin interaction in brain tissues identified multiple genes associated with PAU. We identified existing medications for potential pharmacological studies by drug repurposing analysis. Cross-ancestry polygenic risk scores (PRS) showed better performance in independent sample than single-ancestry PRS. Genetic correlations between PAU and other traits were observed in multiple ancestries, with other substance use traits having the highest correlations. The analysis of diverse ancestries contributed significantly to the findings, and fills an important gap in the literature.

Accuracy of haplotype estimation and whole genome imputation affects complex trait analyses in complex biobanks

Sample recruitment for research consortia, biobanks, and personal genomics companies span years, necessitating genotyping in batches, using different technologies. As marker content on genotyping arrays varies, integrating such datasets is non-trivial and its impact on haplotype estimation (phasing) and whole genome imputation, necessary steps for complex trait analysis, remains under-evaluated. Using the iPSYCH dataset, comprising 130,438 individuals, genotyped in two stages, on different arrays, we evaluated phasing and imputation performance across multiple phasing methods and data integration protocols. While phasing accuracy varied by choice of method and data integration protocol, imputation accuracy varied mostly between data integration protocols. We demonstrate an attenuation in imputation accuracy within samples of non-European origin, highlighting challenges to studying complex traits in diverse populations. Finally, imputation errors can bias association tests, reduce predictive utility of polygenic scores. Carefully optimized data integration strategies enhance accuracy and replicability of complex trait analyses in complex biobanks.

The genetic background of hydrocephalus in a population-based cohort: implication of ciliary involvement

Hydrocephalus is one of the most common congenital disorders of the central nervous system and often displays psychiatric co-morbidities, in particular autism spectrum disorder. The disease mechanisms behind hydrocephalus are complex and not well understood, but some association with dysfunctional cilia in the brain ventricles and subarachnoid space has been indicated. A better understanding of the genetic aetiology of hydrocephalus, including the role of ciliopathies, may bring insights into a potentially shared genetic aetiology. In this population-based case-cohort study, we, for the first time, investigated variants of postulated hydrocephalus candidate genes. Using these data, we aimed to investigate potential involvement of the ciliome in hydrocephalus and describe genotype-phenotype associations with an autism spectrum disorder. One-hundred and twenty-one hydrocephalus candidate genes were screened in a whole-exome-sequenced sub-cohort of the Lundbeck Foundation Initiative for Integrative Psychiatric Research study, comprising 72 hydrocephalus patients and 4181 background population controls. Candidate genes containing high-impact variants of interest were systematically evaluated for their involvement in ciliary function and an autism spectrum disorder. The median age at diagnosis for the hydrocephalus patients was 0 years (range 0-27 years), the median age at analysis was 22 years (11-35 years), and 70.5% were males. The median age for controls was 18 years (range 11-26 years) and 53.3% were males. Fifty-two putative hydrocephalus-associated variants in 34 genes were identified in 42 patients (58.3%). In hydrocephalus cases, we found increased, but not significant, enrichment of high-impact protein altering variants (odds ratio 1.51, 95% confidence interval 0.92-2.51, P = 0.096), which was driven by a significant enrichment of rare protein truncating variants (odds ratio 2.71, 95% confidence interval 1.17-5.58, P = 0.011). Fourteen of the genes with high-impact variants are part of the ciliome, whereas another six genes affect cilia-dependent processes during neurogenesis. Furthermore, 15 of the 34 genes with high-impact variants and three of eight genes with protein truncating variants were associated with an autism spectrum disorder. Because symptoms of other diseases may be neglected or masked by the hydrocephalus-associated symptoms, we suggest that patients with congenital hydrocephalus undergo clinical genetic assessment with respect to ciliopathies and an autism spectrum disorder. Our results point to the significance of hydrocephalus as a ciliary disease in some cases. Future studies in brain ciliopathies may not only reveal new insights into hydrocephalus but also, brain disease in the broadest sense, given the essential role of cilia in neurodevelopment.

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.

Interplay of ADHD Polygenic Liability With Birth-Related, Somatic, and Psychosocial Factors in ADHD: A Nationwide Study

OBJECTIVE

Attention deficit hyperactivity disorder (ADHD) is a multifactorial neurodevelopmental disorder, yet the interplay between ADHD polygenic risk scores (PRSs) and other risk factors remains relatively unexplored. The authors investigated associations, confounding, and interactions of ADHD PRS with birth-related, somatic, and psychosocial factors previously associated with ADHD.

METHODS

Participants included a random general population sample (N=21,578) and individuals diagnosed with ADHD (N=13,697) from the genotyped Danish iPSYCH2012 case cohort, born between 1981 and 2005. The authors derived ADHD PRSs and identified 24 factors previously associated with ADHD using national registers. Logistic regression was used to estimate associations of ADHD PRS with each risk factor in the general population. Cox models were used to evaluate confounding of risk factor associations with ADHD diagnosis by ADHD PRS and parental psychiatric history, and interactions between ADHD PRS and each risk factor.

RESULTS

ADHD PRS was associated with 12 of 24 risk factors (odds ratio range, 1.03-1.30), namely, small gestational age, infections, traumatic brain injury, and most psychosocial risk factors. Nineteen risk factors were associated with ADHD diagnosis (odds ratio range, 1.20-3.68), and adjusting for ADHD PRS and parental psychiatric history led to only minor attenuations. Only the interaction between ADHD PRS and maternal autoimmune disease survived correction for multiple testing.

CONCLUSIONS

Higher ADHD PRS in the general population is associated with small increases in risk for certain birth-related and somatic ADHD risk factors, and broadly to psychosocial adversity. Evidence of gene-environment interaction was limited, as was confounding by ADHD PRS and family psychiatric history on ADHD risk factor associations. This suggests that the majority of the investigated ADHD risk factors act largely independently of current ADHD PRS to increase risk of ADHD.

The impact of the COVID-19 lockdown on birthweight among singleton term birth in Denmark

In Denmark, a nationwide COVID-19 lockdown was implemented on March 12, 2020 and eased on April 14, 2020. The COVID-19 lockdown featured reduced prevalence of extremely preterm or extremely low birthweight births. This study aims to explore the impact of this COVID-19 lockdown on term birthweights in Denmark. We conducted a nationwide register-based cohort study on 27,870 live singleton infants, born at term (weeks 37-41), between March 12 and April 14, 2015-2020, using data from the Danish Neonatal Screening Biobank. Primary outcomes, corrected for confounders, were birthweight, small-for-gestational-age (SGA), and large-for-gestational-age (LGA), comparing the COVID-19 lockdown to the previous five years. Data were analysed using linear regression to assess associations with birthweight. Multinomial logistic regression was used to assess associations with relative-size-for-gestational-age (xGA) categories. Adjusted mean birthweight was significantly increased by 16.9 g (95% CI = 4.1-31.3) during the lockdown period. A dip in mean birthweight was found in gestational weeks 37 and 38 balanced by an increase in weeks 40 and 41. The 2020 lockdown period was associated with an increased LGA prevalence (aOR 1.13, 95% CI = 1.05-1.21). No significant changes in proportions of xGA groups were found between 2015 and 2019. The nationwide COVID-19 lockdown resulted in a small but significant increase in birthweight and proportion of LGA infants, driven by an increase in birthweight in gestational weeks 40 and 41.

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.

Statistical and functional convergence of common and rare genetic influences on autism at chromosome 16p

The canonical paradigm for converting genetic association to mechanism involves iteratively mapping individual associations to the proximal genes through which they act. In contrast, in the present study we demonstrate the feasibility of extracting biological insights from a very large region of the genome and leverage this strategy to study the genetic influences on autism. Using a new statistical approach, we identified the 33-Mb p-arm of chromosome 16 (16p) as harboring the greatest excess of autism's common polygenic influences. The region also includes the mechanistically cryptic and autism-associated 16p11.2 copy number variant. Analysis of RNA-sequencing data revealed that both the common polygenic influences within 16p and the 16p11.2 deletion were associated with decreased average gene expression across 16p. The transcriptional effects of the rare deletion and diffuse common variation were correlated at the level of individual genes and analysis of Hi-C data revealed patterns of chromatin contact that may explain this transcriptional convergence. These results reflect a new approach for extracting biological insight from genetic association data and suggest convergence of common and rare genetic influences on autism at 16p.

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.

Genome-wide study of early and severe childhood asthma identifies interaction between CDHR3 and GSDMB

BACKGROUND

Asthma with severe exacerbation is one of the most common causes of hospitalization among young children. Exacerbations are typically triggered by respiratory infections, but the host factors causing recurrent infections and exacerbations in some children are poorly understood. As a result, current treatment options and preventive measures are inadequate.

OBJECTIVE

We sought to identify genetic interaction associated with the development of childhood asthma.

METHODS

We performed an exhaustive search for pairwise interaction between genetic single nucleotide polymorphisms using 1204 cases of a specific phenotype of early childhood asthma with severe exacerbations in patients aged 2 to 6 years combined with 5328 nonasthmatic controls. Replication was attempted in 3 independent populations, and potential underlying immune mechanisms were investigated in the COPSAC2010 and COPSAC2000 birth cohorts.

RESULTS

We found evidence of interaction, including replication in independent populations, between the known childhood asthma loci CDHR3 and GSDMB. The effect of CDHR3 was dependent on the GSDMB genotype, and this interaction was more pronounced for severe and early onset of disease. Blood immune analyses suggested a mechanism related to increased IL-17A production after viral stimulation.

CONCLUSIONS

We found evidence of interaction between CDHR3 and GSDMB in development of early childhood asthma, possibly related to increased IL-17A response to viral infections. This study demonstrates the importance of focusing on specific disease subtypes for understanding the genetic mechanisms of asthma.

Maternal pregnancy-related infections and autism spectrum disorder-the genetic perspective

Autism spectrum disorder (ASD) refers to a group of neurodevelopmental disorders which include deficits in behavior, social interaction and communication. ASD has a complex genetic architecture, and it is also influenced by certain environmental exposures. Both types of predisposing factors may be related to immunological mechanisms, involving, for example, immune system genes and infections. Past studies have shown an association between infections occurring during the pregnancy in the mother and increased risk of ASD in the child, an observation which has received recent support from experimental animal studies of ASD-like behavior. The aim of this study was to study the genetic contribution to this effect. We employed genetic correlation analyses across potential ASD subtypes stratified on the basis of maternal pregnancy-related infections within the iPSYCH ASD case-cohort sample, as well as a case-case GWAS. We validated the trends of the genetic correlation analyses observed in our sample using GWAS summary statistics from the PGC ASD study (excluding iPSYCH). The genetic correlation between ASD with a history of maternal pregnancy-related infections and ASD without a history of maternal infections in iPSYCH was r_g = 0.3811. We obtained a similar estimate between the former and the PGC ASD phenotype (r_g = 0.3997). Both estimates are lower compared to the genetic correlation between ASD without a history of maternal infections and the PGC ASD phenotype (r_g = 0.6735), and between ASD with a history of maternal infections occurring only more than 2 months following childbirth and the PGC ASD phenotype (r_g = 0.6293). Additionally, we observed genetic variance between the two main ASD phenotypes using summary statistics from the case-case GWAS in iPSYCH (h²_cc = 0.1059), indicating genome-wide differences between the phenotypes. Our results suggest potentially different etiologies of ASD based on a history of maternal pregnancy-related infections, which may, in part, be genetic. This highlights the relevance of maternal pregnancy-related infections to genetic studies of ASD and provides new insights into the molecular underpinnings of 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.

Deep learning-based integration of genetics with registry data for stratification of schizophrenia and depression

Currently, psychiatric diagnoses are, in contrast to most other medical fields, based on subjective symptoms and observable signs and call for new and improved diagnostics to provide the most optimal care. On the basis of a deep learning approach, we performed unsupervised patient stratification of 19,636 patients with depression [major depressive disorder (MDD)] and/or schizophrenia (SCZ) and 22,467 population controls from the iPSYCH2012 case cohort. We integrated data of disorder severity, history of mental disorders and disease comorbidities, genetics, and medical birth data. From this, we stratified the individuals in six and seven unique clusters for MDD and SCZ, respectively. When censoring data until diagnosis, we could predict MDD clusters with areas under the curve (AUCs) of 0.54 to 0.80 and SCZ clusters with AUCs of 0.71 to 0.86. Overall cases and controls could be predicted with an AUC of 0.81, illustrating the utility of data-driven subgrouping in psychiatry.

Women With Turner Syndrome Are Both Estrogen and Androgen Deficient: The Impact of Hormone Replacement Therapy

CONTEXT

Women with Turner syndrome (TS) suffer from hypergonadotropic hypogonadism, causing a deficit in gonadal hormone secretion. As a consequence, these women are treated with estrogen from the age of 12 years, and later in combination with progesterone. However, androgens have been given less attention.

OBJECTIVE

To assess sex hormone levels in women with TS, both those treated and those nontreated with hormone replacement therapy (HRT), and investigate the impact of HRT on sex hormone levels.

METHODS

At Aarhus University Hospital, 99 women with TS were followed 3 times from August 2003 to February 2010. Seventeen were lost during follow-up. Control group 1 consisted of 68 healthy age-matched control women seen once during this period. Control group 2 consisted of 28 young, eumenorrheic women sampled 9 times throughout the same menstrual cycle. Serum concentrations of follicle-stimulating hormone (FSH), luteinizing hormone (LH), 17β-estradiol, estrone sulfate, DHEAS, testosterone, free androgen index, androstenedione, 17-OH progesterone, and sex hormone-binding globulin (SHBG) were analyzed.

RESULTS

All androgens, 17-OH progesterone, and sex hormone-binding globulin (SHBG) were 30% to 50% lower in TS compared with controls (P < 0.01). FSH, LH, and estrone sulfate were more than doubled in women with TS compared with controls (P < 0.02). Using principal component analysis, we describe a positive correlation between women with TS receiving HRT, elevated levels of SHBG, and decreased levels of androgens.

CONCLUSION

The sex hormone profile in TS reveals a picture of androgen deficiency, aggravated further by HRT. Conventional HRT does not normalize estradiol levels in TS.

The female protective effect against autism spectrum disorder

Autism spectrum disorder (ASD) is diagnosed three to four times more frequently in males than in females. Genetic studies of rare variants support a female protective effect (FPE) against ASD. However, sex differences in common inherited genetic risk for ASD are less studied, particularly within families. Leveraging the Danish iPSYCH resource, we found siblings of female ASD cases (n = 1,707) had higher rates of ASD than siblings of male ASD cases (n = 6,270; p < 1.0 × 10-10). In the Simons Simplex and SPARK collections, mothers of ASD cases (n = 7,436) carried more polygenic risk for ASD than fathers of ASD cases (n = 5,926; 0.08 polygenic risk score [PRS] SD; p = 7.0 × 10-7). Further, male unaffected siblings under-inherited polygenic risk (n = 1,519; p = 0.03). Using both epidemiologic and genetic approaches, our findings strongly support an FPE against ASD's common inherited influences.

Blood Biomarkers and 6- to 7-Year Childhood Outcomes Following Neonatal Encephalopathy

OBJECTIVE

This study aimed to profile the cytokine/chemokine response from day 0 to 7 in infants (≥36 weeks of gestational age) with neonatal encephalopathy (NE) and to explore the association with long-term outcomes.

STUDY DESIGN

This was a secondary study of the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD) Neonatal Research Network randomized controlled trial of whole body hypothermia for NE. Eligible infants with moderate-severe NE were randomized to cooling or normothermia. Blood spots were collected on days 0 to 1, 2 to 4, and 6 to 7. Twenty-four cytokines/chemokines were measured using a multiplex platform. Surviving infants underwent neurodevelopmental assessment at 6 to 7 years. Primary outcome was death or moderate-severe impairment defined by any of the following: intelligence quotient <70, moderate-severe cerebral palsy (CP), blindness, hearing impairment, or epilepsy.

RESULTS

Cytokine blood spots were collected from 109 participants. In total 99 of 109 (91%) were assessed at 6 to 7 years; 54 of 99 (55%) developed death/impairment. Neonates who died or were impaired had lower early regulated upon activation normal T cell expressed and secreted (RANTES) and higher day 7 monocyte chemotactic protein (MCP)-1 levels than neonates who survived without impairment. Though TNF-α levels had no association with death/impairment, higher day 0 to 1 levels were observed among neonates who died/developed CP. On multiple regression analysis adjusted for center, treatment group, sex, race, and level of hypoxic ischemic encephalopathy, higher RANTES was inversely associated with death/impairment (odds ratio (OR): 0.31, 95% confidence interval [CI]: 0.13-0.74), while day seven MCP-1 level was directly associated with death/impairment (OR: 3.70, 95% CI: 1.42-9.61). Targeted cytokine/chemokine levels demonstrated little variation with hypothermia treatment.

CONCLUSION

RANTES and MCP-1 levels in the first week of life may provide potential targets for future therapies among neonates with encephalopathy.

KEY POINTS

· Elevation of specific cytokines and chemokines in neonates with encephalopathy has been noted along with increased risk of neurodevelopmental impairment in infancy.. · Cytokine/chemokines at <7 days were assessed among neonates in a trial of hypothermia for HIE.. · Neonates who died or were impaired at 6 to 7 years following hypoxic-ischemic encephalopathy had lower RANTES and higher MCP-1 levels than those who survived without impairment..

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.

Rare coding variants in ten genes confer substantial risk for schizophrenia

Rare coding variation has historically provided the most direct connections between gene function and disease pathogenesis. By meta-analysing the whole exomes of 24,248 schizophrenia cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in 10 genes as conferring substantial risk for schizophrenia (odds ratios of 3-50, P < 2.14 × 10-6) and 32 genes at a false discovery rate of <5%. These genes have the greatest expression in central nervous system neurons and have diverse molecular functions that include the formation, structure and function of the synapse. The associations of the NMDA (N-methyl-D-aspartate) receptor subunit GRIN2A and AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptor subunit GRIA3 provide support for dysfunction of the glutamatergic system as a mechanistic hypothesis in the pathogenesis of schizophrenia. We observe an overlap of rare variant risk among schizophrenia, autism spectrum disorders1, epilepsy and severe neurodevelopmental disorders2, although different mutation types are implicated in some shared genes. Most genes described here, however, are not implicated in neurodevelopment. We demonstrate that genes prioritized from common variant analyses of schizophrenia are enriched in rare variant risk3, suggesting that common and rare genetic risk factors converge at least partially on the same underlying pathogenic biological processes. Even after excluding significantly associated genes, schizophrenia cases still carry a substantial excess of URVs, which indicates that more risk genes await discovery using this approach.

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.

Accounting for age of onset and family history improves power in genome-wide association studies

Genome-wide association studies (GWASs) have revolutionized human genetics, allowing researchers to identify thousands of disease-related genes and possible drug targets. However, case-control status does not account for the fact that not all controls may have lived through their period of risk for the disorder of interest. This can be quantified by examining the age-of-onset distribution and the age of the controls or the age of onset for cases. The age-of-onset distribution may also depend on information such as sex and birth year. In addition, family history is not routinely included in the assessment of control status. Here, we present LT-FH++, an extension of the liability threshold model conditioned on family history (LT-FH), which jointly accounts for age of onset and sex as well as family history. Using simulations, we show that, when family history and the age-of-onset distribution are available, the proposed approach yields statistically significant power gains over LT-FH and large power gains over genome-wide association study by proxy (GWAX). We applied our method to four psychiatric disorders available in the iPSYCH data and to mortality in the UK Biobank and found 20 genome-wide significant associations with LT-FH++, compared to ten for LT-FH and eight for a standard case-control GWAS. As more genetic data with linked electronic health records become available to researchers, we expect methods that account for additional health information, such as LT-FH++, to become even more beneficial.

Dissecting the Shared Genetic Architecture of Suicide Attempt, Psychiatric Disorders, and Known Risk Factors

BACKGROUND

Suicide is a leading cause of death worldwide, and nonfatal suicide attempts, which occur far more frequently, are a major source of disability and social and economic burden. Both have substantial genetic etiology, which is partially shared and partially distinct from that of related psychiatric disorders.

METHODS

We conducted a genome-wide association study (GWAS) of 29,782 suicide attempt (SA) cases and 519,961 controls in the International Suicide Genetics Consortium (ISGC). The GWAS of SA was conditioned on psychiatric disorders using GWAS summary statistics via multitrait-based conditional and joint analysis, to remove genetic effects on SA mediated by psychiatric disorders. We investigated the shared and divergent genetic architectures of SA, psychiatric disorders, and other known risk factors.

RESULTS

Two loci reached genome-wide significance for SA: the major histocompatibility complex and an intergenic locus on chromosome 7, the latter of which remained associated with SA after conditioning on psychiatric disorders and replicated in an independent cohort from the Million Veteran Program. This locus has been implicated in risk-taking behavior, smoking, and insomnia. SA showed strong genetic correlation with psychiatric disorders, particularly major depression, and also with smoking, pain, risk-taking behavior, sleep disturbances, lower educational attainment, reproductive traits, lower socioeconomic status, and poorer general health. After conditioning on psychiatric disorders, the genetic correlations between SA and psychiatric disorders decreased, whereas those with nonpsychiatric traits remained largely unchanged.

CONCLUSIONS

Our results identify a risk locus that contributes more strongly to SA than other phenotypes and suggest a shared underlying biology between SA and known risk factors that is not mediated by psychiatric disorders.

Genome-wide association study of febrile seizures implicates fever response and neuronal excitability genes

Febrile seizures represent the most common type of pathological brain activity in young children and are influenced by genetic, environmental and developmental factors. In a minority of cases, febrile seizures precede later development of epilepsy.

We conducted a genome-wide association study of febrile seizures in 7635 cases and 83 966 controls identifying and replicating seven new loci, all with P < 5 × 10⁻¹⁰.

Variants at two loci were functionally related to altered expression of the fever response genes PTGER3 and IL10, and four other loci harboured genes (BSN, ERC2, GABRG2, HERC1) influencing neuronal excitability by regulating neurotransmitter release and binding, vesicular transport or membrane trafficking at the synapse. Four previously reported loci (SCN1A, SCN2A, ANO3 and 12q21.33) were all confirmed. Collectively, the seven novel and four previously reported loci explained 2.8% of the variance in liability to febrile seizures, and the single nucleotide polymorphism heritability based on all common autosomal single nucleotide polymorphisms was 10.8%. GABRG2, SCN1A and SCN2A are well-established epilepsy genes and, overall, we found positive genetic correlations with epilepsies (r_g = 0.39, P = 1.68 × 10⁻⁴). Further, we found that higher polygenic risk scores for febrile seizures were associated with epilepsy and with history of hospital admission for febrile seizures. Finally, we found that polygenic risk of febrile seizures was lower in febrile seizure patients with neuropsychiatric disease compared to febrile seizure patients in a general population sample.

In conclusion, this largest genetic investigation of febrile seizures to date implicates central fever response genes as well as genes affecting neuronal excitability, including several known epilepsy genes. Further functional and genetic studies based on these findings will provide important insights into the complex pathophysiological processes of seizures with and without fever.

Comparing Copy Number Variations in a Danish Case Cohort of Individuals With Psychiatric Disorders

Importance

Although the association between several recurrent genomic copy number variants (CNVs) and mental disorders has been studied for more than a decade, unbiased, population-based estimates of the prevalence, disease risks and trajectories, fertility, and mortality to contrast chromosomal abnormalities and advance precision health care are lacking.

Objective

To generate unbiased, population-based estimates of prevalence, disease risks and trajectories, fertility, and mortality of CNVs implicated in neuropsychiatric disorders.

Design, Setting, and Participants

In a population-based case-cohort study, using the Lundbeck Foundation Initiative for Integrative Psychiatric Research (iPSYCH) 2012 database, individuals born between May 1, 1981, and December 31, 2005, and followed up until December 31, 2012, were analyzed. All individuals (n = 57 377) with attention-deficit/hyperactivity disorder (ADHD), major depressive disorder (MDD), schizophrenia (SCZ), autism spectrum disorder (ASD), or bipolar disorder (BPD) were included, as well as 30 000 individuals randomly drawn from the database. Data analysis was conducted from July 1, 2017, to September 7, 2021.

Exposures

Copy number variants at 6 genomic loci (1q21.1, 15q11.2, 15q13.3, 16p11.2, 17p12, and 17q12).

Main Outcomes and Measures

Population-unbiased hazard ratio (HR) and survival estimates of CNV associations with the 5 ascertained psychiatric disorders, epilepsy, intellectual disability, selected somatic disorders, fertility, and mortality.

Results

Participants' age ranged from 1 to 32 years (mean, 12.0 [IQR, 6.9] years) during follow-up, and 38 662 were male (52.3%). Copy number variants broadly associated with an increased risk of autism spectrum disorder and ADHD, whereas risk estimates of SCZ for most CNVs were lower than previously reported. Comparison with previous studies suggests that the lower risk estimates are associated with a higher CNV prevalence in the general population than in control samples of most case-control studies. Significant risk of major depressive disorder (HR, 5.8; 95% CI, 1.5-22.2) and sex-specific risk of bipolar disorder (HR, 17; 95% CI, 1.5-189.3, in men only) were noted for the 1q21.1 deletion. Although CNVs at 1q21.1 and 15q13.3 were associated with increased risk across most diagnoses, the 17p12 deletion consistently conferred less risk of psychiatric disorders (HR 0.4-0.8), although none of the estimates differed significantly from the general population. Trajectory analyses noted that, although diagnostic risk profiles differed across loci, they were similar for deletions and duplications within each locus. Sex-stratified analyses suggest that pathogenicity of many CNVs may be modulated by sex.

Conclusions and Relevance

The findings of this study suggest that the iPSYCH population case cohort reveals broad disease risk for some of the studied CNVs and narrower risk for others, in addition to sex differential liability. This finding on genomic risk variants at the level of a population may be important for health care planning and clinical decision making, and thus the advancement of precision health care.

Preterm birth, stillbirth and early neonatal mortality during the Danish COVID-19 lockdown

Using provisional or opportunistic data, three nationwide studies (The Netherlands, the USA and Denmark) have identified a reduction in preterm or extremely preterm births during periods of COVID-19 restrictions. However, none of the studies accounted for perinatal deaths. To determine whether the reduction in extremely preterm births, observed in Denmark during the COVID-19 lockdown, could be the result of an increase in perinatal deaths and to assess the impact of extended COVID-19 restrictions, we performed a nationwide Danish register-based prevalence proportion study. We examined all singleton pregnancies delivered in Denmark during the COVID-19 strict lockdown calendar periods (March 12-April 14, 2015-2020, N = 31,164 births) and the extended calendar periods of COVID-19 restrictions (February 27-September 30, 2015-2020, N = 214,862 births). The extremely preterm birth rate was reduced (OR 0.27, 95% CI 0.07 to 0.86) during the strict lockdown period in 2020, while perinatal mortality was not significantly different. During the extended period of restrictions in 2020, the extremely preterm birth rate was marginally reduced, and a significant reduction in the stillbirth rate (OR 0.69, 0.50 to 0.95) was observed. No changes in early neonatal mortality rates were found.Conclusion: Stillbirth and extremely preterm birth rates were reduced in Denmark during the period of COVID-19 restrictions and lockdown, respectively, suggesting that aspects of these containment and control measures confer an element of protection. The present observational study does not allow for causal inference; however, the results support the design of studies to ascertain whether behavioural or social changes for pregnant women may improve pregnancy outcomes. What is Known: • The aetiologies of preterm birth and stillbirth are multifaceted and linked to a wide range of socio-demographic, medical, obstetric, foetal, psychosocial and environmental factors. • The COVID-19 lockdown saw a reduction in extremely preterm births in Denmark and other high-income countries. An urgent question is whether this reduction can be explained by increased perinatal mortality. What is New: • The reduction in extremely preterm births during the Danish COVID-19 lockdown was not a consequence of increased perinatal mortality, which remained unchanged during this period. • The stillbirth rate was reduced throughout the extended period of COVID-19 restrictions.