Multi-polygenic scores in psychiatry: From disorder specific to transdiagnostic perspectives.

The dense co-occurrence of psychiatric disorders questions the categorical classification tradition and motivates efforts to establish dimensional constructs with neurobiological foundations that transcend diagnostic boundaries. In this study, we examined the genetic liability for eight major psychiatric disorder phenotypes under both a disorder-specific and a transdiagnostic framework. The study sample (n = 513) was deeply phenotyped, consisting of 452 patients from tertiary care with mood disorders, anxiety disorders (ANX), attention-deficit/hyperactivity disorder (ADHD), autism spectrum disorders, and/or substance use disorders (SUD) and 61 unaffected comparison individuals. We computed subject-specific polygenic risk score (PRS) profiles and assessed their associations with psychiatric diagnoses, comorbidity status, as well as cross-disorder behavioral dimensions derived from a rich battery of psychopathology assessments. High PRSs for depression were unselectively associated with the diagnosis of SUD, ADHD, ANX, and mood disorders (p < 1e-4). In the dimensional approach, four distinct functional domains were uncovered, namely the negative valence, social, cognitive, and regulatory systems, closely matching the major functional domains proposed by the Research Domain Criteria (RDoC) framework. Critically, the genetic predisposition for depression was selectively reflected in the functional aspect of negative valence systems (R2 = 0.041, p = 5e-4) but not others. This study adds evidence to the ongoing discussion about the misalignment between current psychiatric nosology and the underlying psychiatric genetic etiology and underscores the effectiveness of the dimensional approach in both the functional characterization of psychiatric patients and the delineation of the genetic liability for psychiatric disorders.

Predicting ADHD in alcohol dependence using polygenic risk scores for ADHD.

Attention deficit hyperactivity disorder (ADHD) is a neurodevelopmental disorder with a high degree of comorbidity, including substance misuse. We aimed to assess whether ADHD polygenic risk scores (PRS) could predict ADHD diagnosis in alcohol dependence (AD). ADHD PRS were generated for 1223 AD subjects with ADHD diagnosis information and 1818 healthy controls. ADHD PRS distributions were compared to evaluate the differences between healthy controls and AD cases with and without ADHD. We found increased ADHD PRS means in the AD cohort with ADHD (mean 0.30, standard deviation (SD) 0.92; p = 3.9 × 10-6 ); and without ADHD (mean - 0.00, SD 1.00; p = 5.2 × 10-5 ) compared to the healthy control subjects (mean - 0.17, SD 0.99). The ADHD PRS means differed within the AD group with a higher ADHD PRS mean in those with ADHD, odds ratio (OR) 1.34, confidence interval (CI) 1.10 to 1.65; p = 0.002. This study showed a positive relationship between ADHD PRS and risk of ADHD in individuals with co-occurring AD indicating that ADHD PRS may have utility in identifying individuals that are at a higher or lower risk of ADHD. Further larger studies need to be conducted to confirm the reliability of the results before ADHD PRS can be considered as a robust biomarker for diagnosis.

Identification of genetic overlap and novel risk loci for attention-deficit/hyperactivity disorder and bipolar disorder.

Differential diagnosis between childhood onset attention-deficit/hyperactivity disorder (ADHD) and bipolar disorder (BD) remains a challenge, mainly due to overlapping symptoms and high rates of comorbidity. Despite this, genetic correlation reported for these disorders is low and non-significant. Here we aimed to better characterize the genetic architecture of these disorders utilizing recent large genome-wide association studies (GWAS). We analyzed independent GWAS summary statistics for ADHD (19,099 cases and 34,194 controls) and BD (20,352 cases and 31,358 controls) applying the conditional/conjunctional false discovery rate (condFDR/conjFDR) statistical framework that increases the power to detect novel phenotype-specific and shared loci by leveraging the combined power of two GWAS. We observed cross-trait polygenic enrichment for ADHD conditioned on associations with BD, and vice versa. Leveraging this enrichment, we identified 19 novel ADHD risk loci and 40 novel BD risk loci at condFDR <0.05. Further, we identified five loci jointly associated with ADHD and BD (conjFDR < 0.05). Interestingly, these five loci show concordant directions of effect for ADHD and BD. These results highlight a shared underlying genetic risk for ADHD and BD which may help to explain the high comorbidity rates and difficulties in differentiating between ADHD and BD in the clinic. Improving our understanding of the underlying genetic architecture of these disorders may aid in the development of novel stratification tools to help reduce these diagnostic difficulties.

CNVs conferring risk of autism or schizophrenia affect cognition in controls.

In a small fraction of patients with schizophrenia or autism, alleles of copy-number variants (CNVs) in their genomes are probably the strongest factors contributing to the pathogenesis of the disease. These CNVs may provide an entry point for investigations into the mechanisms of brain function and dysfunction alike. They are not fully penetrant and offer an opportunity to study their effects separate from that of manifest disease. Here we show in an Icelandic sample that a few of the CNVs clearly alter fecundity (measured as the number of children by age 45). Furthermore, we use various tests of cognitive function to demonstrate that control subjects carrying the CNVs perform at a level that is between that of schizophrenia patients and population controls. The CNVs do not all affect the same cognitive domains, hence the cognitive deficits that drive or accompany the pathogenesis vary from one CNV to another. Controls carrying the chromosome 15q11.2 deletion between breakpoints 1 and 2 (15q11.2(BP1-BP2) deletion) have a history of dyslexia and dyscalculia, even after adjusting for IQ in the analysis, and the CNV only confers modest effects on other cognitive traits. The 15q11.2(BP1-BP2) deletion affects brain structure in a pattern consistent with both that observed during first-episode psychosis in schizophrenia and that of structural correlates in dyslexia.

Nonsense mutation in the LGR4 gene is associated with several human diseases and other traits.

Low bone mineral density (BMD) is used as a parameter of osteoporosis. Genome-wide association studies of BMD have hitherto focused on BMD as a quantitative trait, yielding common variants of small effects that contribute to the population diversity in BMD. Here we use BMD as a dichotomous trait, searching for variants that may have a direct effect on the risk of pathologically low BMD rather than on the regulation of BMD in the healthy population. Through whole-genome sequencing of Icelandic individuals, we found a rare nonsense mutation within the leucine-rich-repeat-containing G-protein-coupled receptor 4 (LGR4) gene (c.376C>T) that is strongly associated with low BMD, and with osteoporotic fractures. This mutation leads to termination of LGR4 at position 126 and fully disrupts its function. The c.376C>T mutation is also associated with electrolyte imbalance, late onset of menarche and reduced testosterone levels, as well as an increased risk of squamous cell carcinoma of the skin and biliary tract cancer. Interestingly, the phenotype of carriers of the c.376C>T mutation overlaps that of Lgr4 mutant mice.

Rate of de novo mutations and the importance of father's age to disease risk.

Mutations generate sequence diversity and provide a substrate for selection. The rate of de novo mutations is therefore of major importance to evolution. Here we conduct a study of genome-wide mutation rates by sequencing the entire genomes of 78 Icelandic parent-offspring trios at high coverage. We show that in our samples, with an average father's age of 29.7, the average de novo mutation rate is 1.20 × 10(-8) per nucleotide per generation. Most notably, the diversity in mutation rate of single nucleotide polymorphisms is dominated by the age of the father at conception of the child. The effect is an increase of about two mutations per year. An exponential model estimates paternal mutations doubling every 16.5 years. After accounting for random Poisson variation, father's age is estimated to explain nearly all of the remaining variation in the de novo mutation counts. These observations shed light on the importance of the father's age on the risk of diseases such as schizophrenia and autism.

Fine-scale recombination rate differences between sexes, populations and individuals.

Meiotic recombinations contribute to genetic diversity by yielding new combinations of alleles. Recently, high-resolution recombination maps were inferred from high-density single-nucleotide polymorphism (SNP) data using linkage disequilibrium (LD) patterns that capture historical recombination events. The use of these maps has been demonstrated by the identification of recombination hotspots and associated motifs, and the discovery that the PRDM9 gene affects the proportion of recombinations occurring at hotspots. However, these maps provide no information about individual or sex differences. Moreover, locus-specific demographic factors like natural selection can bias LD-based estimates of recombination rate. Existing genetic maps based on family data avoid these shortcomings, but their resolution is limited by relatively few meioses and a low density of markers. Here we used genome-wide SNP data from 15,257 parent-offspring pairs to construct the first recombination maps based on directly observed recombinations with a resolution that is effective down to 10 kilobases (kb). Comparing male and female maps reveals that about 15% of hotspots in one sex are specific to that sex. Although male recombinations result in more shuffling of exons within genes, female recombinations generate more new combinations of nearby genes. We discover novel associations between recombination characteristics of individuals and variants in the PRDM9 gene and we identify new recombination hotspots. Comparisons of our maps with two LD-based maps inferred from data of HapMap populations of Utah residents with ancestry from northern and western Europe (CEU) and Yoruba in Ibadan, Nigeria (YRI) reveal population differences previously masked by noise and map differences at regions previously described as targets of natural selection.

A variant in CDKAL1 influences insulin response and risk of type 2 diabetes.

We conducted a genome-wide association study for type 2 diabetes (T2D) in Icelandic cases and controls, and we found that a previously described variant in the transcription factor 7-like 2 gene (TCF7L2) gene conferred the most significant risk. In addition to confirming two recently identified risk variants, we identified a variant in the CDKAL1 gene that was associated with T2D in individuals of European ancestry (allele-specific odds ratio (OR) = 1.20 (95% confidence interval, 1.13-1.27), P = 7.7 x 10(-9)) and individuals from Hong Kong of Han Chinese ancestry (OR = 1.25 (1.11-1.40), P = 0.00018). The genotype OR of this variant suggested that the effect was substantially stronger in homozygous carriers than in heterozygous carriers. The ORs for homozygotes were 1.50 (1.31-1.72) and 1.55 (1.23-1.95) in the European and Hong Kong groups, respectively. The insulin response for homozygotes was approximately 20% lower than for heterozygotes or noncarriers, suggesting that this variant confers risk of T2D through reduced insulin secretion.

Variant of transcription factor 7-like 2 (TCF7L2) gene confers risk of type 2 diabetes.

We have previously reported suggestive linkage of type 2 diabetes mellitus to chromosome 10q. We genotyped 228 microsatellite markers in Icelandic individuals with type 2 diabetes and controls throughout a 10.5-Mb interval on 10q. A microsatellite, DG10S478, within intron 3 of the transcription factor 7-like 2 gene (TCF7L2; formerly TCF4) was associated with type 2 diabetes (P = 2.1 x 10(-9)). This was replicated in a Danish cohort (P = 4.8 x 10(-3)) and in a US cohort (P = 3.3 x 10(-9)). Compared with non-carriers, heterozygous and homozygous carriers of the at-risk alleles (38% and 7% of the population, respectively) have relative risks of 1.45 and 2.41. This corresponds to a population attributable risk of 21%. The TCF7L2 gene product is a high mobility group box-containing transcription factor previously implicated in blood glucose homeostasis. It is thought to act through regulation of proglucagon gene expression in enteroendocrine cells via the Wnt signaling pathway.

Genetics of gene expression and its effect on disease.

Common human diseases result from the interplay of many genes and environmental factors. Therefore, a more integrative biology approach is needed to unravel the complexity and causes of such diseases. To elucidate the complexity of common human diseases such as obesity, we have analysed the expression of 23,720 transcripts in large population-based blood and adipose tissue cohorts comprehensively assessed for various phenotypes, including traits related to clinical obesity. In contrast to the blood expression profiles, we observed a marked correlation between gene expression in adipose tissue and obesity-related traits. Genome-wide linkage and association mapping revealed a highly significant genetic component to gene expression traits, including a strong genetic effect of proximal (cis) signals, with 50% of the cis signals overlapping between the two tissues profiled. Here we demonstrate an extensive transcriptional network constructed from the human adipose data that exhibits significant overlap with similar network modules constructed from mouse adipose data. A core network module in humans and mice was identified that is enriched for genes involved in the inflammatory and immune response and has been found to be causally associated to obesity-related traits.

GWAS meta-analysis reveals key risk loci in essential tremor pathogenesis.

Essential tremor (ET) is a prevalent neurological disorder with a largely unknown underlying biology. In this genome-wide association study meta-analysis, comprising 16,480 ET cases and 1,936,173 controls from seven datasets, we identify 12 sequence variants at 11 loci. Evaluating mRNA expression, splicing, plasma protein levels, and coding effects, we highlight seven putative causal genes at these loci, including CA3 and CPLX1. CA3 encodes Carbonic Anhydrase III and carbonic anhydrase inhibitors have been shown to decrease tremors. CPLX1, encoding Complexin-1, regulates neurotransmitter release. Through gene-set enrichment analysis, we identify a significant association with specific cell types, including dopaminergic and GABAergic neurons, as well as biological processes like Rho GTPase signaling. Genetic correlation analyses reveals a positive association between ET and Parkinson's disease, depression, and anxiety-related phenotypes. This research uncovers risk loci, enhancing our knowledge of the complex genetics of this common but poorly understood disorder, and highlights CA3 and CPLX1 as potential therapeutic targets.

Sequence variants at CYP1A1-CYP1A2 and AHR associate with coffee consumption.

Coffee is the most commonly used stimulant and caffeine is its main psychoactive ingredient. The heritability of coffee consumption has been estimated at around 50%. We performed a meta-analysis of four genome-wide association studies of coffee consumption among coffee drinkers from Iceland (n = 2680), The Netherlands (n = 2791), the Sorbs Slavonic population isolate in Germany (n = 771) and the USA (n = 369) using both directly genotyped and imputed single nucleotide polymorphisms (SNPs) (2.5 million SNPs). SNPs at the two most significant loci were also genotyped in a sample set from Iceland (n = 2430) and a Danish sample set consisting of pregnant women (n = 1620). Combining all data, two sequence variants significantly associated with increased coffee consumption: rs2472297-T located between CYP1A1 and CYP1A2 at 15q24 (P = 5.4 · 10(-14)) and rs6968865-T near aryl hydrocarbon receptor (AHR) at 7p21 (P = 2.3 · 10(-11)). An effect of ∼0.2 cups a day per allele was observed for both SNPs. CYP1A2 is the main caffeine metabolizing enzyme and is also involved in drug metabolism. AHR detects xenobiotics, such as polycyclic aryl hydrocarbons found in roasted coffee, and induces transcription of CYP1A1 and CYP1A2. The association of these SNPs with coffee consumption was present in both smokers and non-smokers.

Common genetic variants associated with open-angle glaucoma.

Open-angle glaucoma (glaucoma) is a major eye disorder characterized by optic disc pathology. Recent genome-wide association studies identified new loci associated with clinically relevant optic disc parameters, such as the optic disc area and vertical cup-disc ratio (VCDR). We examined to what extent these loci are involved in glaucoma. The loci studied include ATOH7, CDC7/TGFBR3 and SALL1 for optic disc area, and CDKN2B, SIX1, SCYL1/LTBP3, CHEK2, ATOH7 and DCLK1 for VCDR. We performed a meta-analysis using data from six independent studies including: the Rotterdam Study (n= 5736), Genetic Research in Isolated Populations combined with Erasmus Rucphen Family study (n= 1750), Amsterdam Glaucoma Study (n= 296) and cohorts from Erlangen and Tübingen (n= 1363), Southampton (n= 702) and deCODE (n= 36 151) resulting in a total of 3161 glaucoma cases and 42 837 controls. Of the eight loci, we found significant evidence (P= 1.41 × 10(-8)) for the association of CDKN2B with glaucoma [odds ratio (OR) for those homozygous for the risk allele: 0.76; 95% confidence interval (CI): 0.70-0.84], for the role of ATOH7 (OR: 1.28; 95% CI: 1.12-1.47) and for SIX1 (OR: 1.20; 95% CI: 1.10-1.31) when adjusting for the number of tested loci. Furthermore, there was a borderline significant association of CDC7/TGFBR3 and SALL1 (both P= 0.04) with glaucoma. In conclusion, we found consistent evidence for three common variants (CDKN2B, ATOH7 and SIX1) significantly associated with glaucoma. These findings may shed new light on the pathophysiological protein pathways leading to glaucoma, and point to pathways involved in the growth and development of the optic nerve.

DanMAC5: a browser of aggregated sequence variants from 8,671 whole genome sequenced Danish individuals.

OBJECTIVES: Allele counts of sequence variants obtained by whole genome sequencing (WGS) often play a central role in interpreting the results of genetic and genomic research. However, such variant counts are not readily available for individuals in the Danish population. Here, we present a dataset with allele counts for sequence variants (single nucleotide variants (SNVs) and indels) identified from WGS of 8,671 (5,418 females) individuals from the Danish population. The data resource is based on WGS data from three independent research projects aimed at assessing genetic risk factors for cardiovascular, psychiatric, and headache disorders. To enable the sharing of information on sequence variation in Danish individuals, we created summarized statistics on allele counts from anonymized data and made them available through the European Genome-phenome Archive (EGA, https://identifiers.org/ega. DATASET: EGAD00001009756 ) and in a dedicated browser, DanMAC5 (available at www.danmac5.dk ). The summary level data and the DanMAC5 browser provide insight into the allelic spectrum of sequence variants segregating in the Danish population, which is important in variant interpretation. DATA DESCRIPTION: Three WGS datasets with an average coverage of 30x were processed independently using the same quality control pipeline. Subsequently, we summarized, filtered, and merged allele counts to create a high-quality summary level dataset of sequence variants.

Deficit of homozygosity among 1.52 million individuals and genetic causes of recessive lethality.

Genotypes causing pregnancy loss and perinatal mortality are depleted among living individuals and are therefore difficult to find. To explore genetic causes of recessive lethality, we searched for sequence variants with deficit of homozygosity among 1.52 million individuals from six European populations. In this study, we identified 25 genes harboring protein-altering sequence variants with a strong deficit of homozygosity (10% or less of predicted homozygotes). Sequence variants in 12 of the genes cause Mendelian disease under a recessive mode of inheritance, two under a dominant mode, but variants in the remaining 11 have not been reported to cause disease. Sequence variants with a strong deficit of homozygosity are over-represented among genes essential for growth of human cell lines and genes orthologous to mouse genes known to affect viability. The function of these genes gives insight into the genetics of intrauterine lethality. We also identified 1077 genes with homozygous predicted loss-of-function genotypes not previously described, bringing the total set of genes completely knocked out in humans to 4785.

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.

Rare variants with large effects provide functional insights into the pathology of migraine subtypes, with and without aura.

Migraine is a complex neurovascular disease with a range of severity and symptoms, yet mostly studied as one phenotype in genome-wide association studies (GWAS). Here we combine large GWAS datasets from six European populations to study the main migraine subtypes, migraine with aura (MA) and migraine without aura (MO). We identified four new MA-associated variants (in PRRT2, PALMD, ABO and LRRK2) and classified 13 MO-associated variants. Rare variants with large effects highlight three genes. A rare frameshift variant in brain-expressed PRRT2 confers large risk of MA and epilepsy, but not MO. A burden test of rare loss-of-function variants in SCN11A, encoding a neuron-expressed sodium channel with a key role in pain sensation, shows strong protection against migraine. Finally, a rare variant with cis-regulatory effects on KCNK5 confers large protection against migraine and brain aneurysms. Our findings offer new insights with therapeutic potential into the complex biology of migraine and its subtypes.

A genome-wide meta-analysis identifies 50 genetic loci associated with carpal tunnel syndrome.

Carpal tunnel syndrome (CTS) is the most common entrapment neuropathy and has a largely unknown underlying biology. In a genome-wide association study of CTS (48,843 cases and 1,190,837 controls), we found 53 sequence variants at 50 loci associated with the syndrome. The most significant association is with a missense variant (p.Glu366Lys) in SERPINA1 that protects against CTS (P = 2.9 × 10-24, OR = 0.76). Through various functional analyses, we conclude that at least 22 genes mediate CTS risk and highlight the role of 19 CTS variants in the biology of the extracellular matrix. We show that the genetic component to the risk is higher in bilateral/recurrent/persistent cases than nonrecurrent/nonpersistent cases. Anthropometric traits including height and BMI are genetically correlated with CTS, in addition to early hormonal-replacement therapy, osteoarthritis, and restlessness. Our findings suggest that the components of the extracellular matrix play a key role in the pathogenesis of CTS.

Thirty novel sequence variants impacting human intracranial volume.

Intracranial volume, measured through magnetic resonance imaging and/or estimated from head circumference, is heritable and correlates with cognitive traits and several neurological disorders. We performed a genome-wide association study meta-analysis of intracranial volume (n = 79 174) and found 64 associating sequence variants explaining 5.0% of its variance. We used coding variation, transcript and protein levels, to uncover 12 genes likely mediating the effect of these variants, including GLI3 and CDK6 that affect cranial synostosis and microcephaly, respectively. Intracranial volume correlates genetically with volumes of cortical and sub-cortical regions, cognition, learning, neonatal and neurological traits. Parkinson's disease cases have greater and attention deficit hyperactivity disorder cases smaller intracranial volume than controls. Our Mendelian randomization studies indicate that intracranial volume associated variants either increase the risk of Parkinson's disease and decrease the risk of attention deficit hyperactivity disorder and neuroticism or correlate closely with a confounder.