Exploring the Genetic Risk of Childhood Daytime Urinary Incontinence: A Genome-Wide Association Study.

PURPOSE: Childhood incontinence is stigmatized and underprioritized, and a basic understanding of its pathogenesis is missing. Our goal was to identify risk-conferring genetic variants in daytime urinary incontinence (DUI). MATERIALS AND METHODS: We conducted a genome-wide association study in the Danish iPSYCH2015 cohort. Cases (3024) were identified through DUI diagnosis codes and redeemed prescriptions for DUI medication in patients aged 5 to 20 years. Controls (30,240), selected from the same sample, were matched to cases on sex and psychiatric diagnoses, if any, and down-sampled to a 1:10 case:control ratio. Replication was performed in the Icelandic deCODE cohort (5475 cases/287,773 controls). Single-nucleotide polymorphism heritability was calculated using the genome-based restricted maximum likelihood method. Cross-trait genetic correlation was estimated using linkage disequilibrium score regression. Polygenic risk scores generated with LDpred2-auto and BOLT-LMM were assessed for association. RESULTS: Variants on chromosome 6 (rs12210989, odds ratio [OR] 1.24, 95% CI 1.17-1.32, P = 3.21 × 10-12) and 20 (rs4809801, OR 1.18, 95% CI 1.11-1.25, P = 3.66 × 10-8) reached genome-wide significance and implicated the PRDM13 and RIPOR3 genes. Chromosome 6 findings were replicated (P = .024, OR 1.09, 95% CI 1.01-1.16). Liability scale heritability ranged from 10.20% (95% CI 6.40%-14.00%) to 15.30% (95% CI 9.66%-20.94%). DUI and nocturnal enuresis showed positive genetic correlation (rg = 1.28 ± 0.38, P = .0007). DUI was associated with attention-deficit/hyperactivity disorder (OR 1.098, 95% CI 1.046-1.152, P < .0001) and BMI (OR 1.129, 95% CI 1.081-1.178, P < .0001) polygenic risk. CONCLUSIONS: Common genetic variants contribute to the risk of childhood DUI, and genes important in neuronal development and detrusor smooth muscle activity were implicated. These findings may help guide identification of new treatment targets.

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.

Loss-of-function variants in ITSN1 confer high risk of Parkinson's disease.

Parkinson's disease (PD) is a debilitating neurodegenerative disorder and its rising global incidence highlights the need for the identification of modifiable risk factors. In a gene-based burden test of rare variants (8647 PD cases and 777,693 controls) we discovered a novel association between loss-of-function variants in ITSN1 and PD. This association was further supported with burden data from the Neurodegenerative Disease Knowledge Portal and the Accelerating Medicines Partnership Parkinson's Disease Knowledge Platform. Our findings show that Rho GTPases and disruptions in synaptic vesicle transport may be involved in the pathogenesis of PD, pointing to the possibility of novel therapeutic approaches.

Homozygosity for a stop-gain variant in CCDC201 causes primary ovarian insufficiency.

Age at menopause (AOM) has a substantial impact on fertility and disease risk. While many loci with variants that associate with AOM have been identified through genome-wide association studies (GWAS) under an additive model, other genetic models are rarely considered1. Here through GWAS meta-analysis under the recessive model of 174,329 postmenopausal women from Iceland, Denmark, the United Kingdom (UK; UK Biobank) and Norway, we study low-frequency variants with a large effect on AOM. We discovered that women homozygous for the stop-gain variant rs117316434 (A) in CCDC201 (p.(Arg162Ter), minor allele frequency ~1%) reached menopause 9 years earlier than other women (P = 1.3 × 10-15). The genotype is present in one in 10,000 northern European women and leads to primary ovarian insufficiency in close to half of them. Consequently, homozygotes have fewer children, and the age at last childbirth is 5 years earlier (P = 3.8 × 10-5). The CCDC201 gene was only found in humans in 2022 and is highly expressed in oocytes. Homozygosity for CCDC201 loss-of-function has a substantial impact on female reproductive health, and homozygotes would benefit from reproductive counseling and treatment for symptoms of early menopause.

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.

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.