ABSTRACT
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.
ABSTRACT
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.
ABSTRACT
Antidepressants exhibit a considerable variation in efficacy, and increasing evidence suggests that individual genetics contribute to antidepressant treatment response. Here, we combined data on antidepressant non-response measured using rating scales for depressive symptoms, questionnaires of treatment effect, and data from electronic health records, to increase statistical power to detect genomic loci associated with non-response to antidepressants in a total sample of 135,471 individuals prescribed antidepressants. We performed genome-wide association meta-analyses, leave-one-out polygenic prediction, and bioinformatics analyses for genetically informed drug prioritization. We identified two novel loci associated with non-response to antidepressants and showed significant polygenic prediction in independent samples. In addition, we investigated drugs that target proteins likely involved in mechanisms underlying antidepressant non-response, and shortlisted drugs that warrant further replication and validation of their potential to reduce depressive symptoms in individuals who do not respond to first-line antidepressant medications. These results suggest that meta-analyses of GWAS utilizing real-world measures of treatment outcomes can increase sample sizes to improve the discovery of variants associated with non-response to antidepressants.
ABSTRACT
AIMS/HYPOTHESIS: Metabolic risk factors and plasma biomarkers for diabetes have previously been shown to change prior to a clinical diabetes diagnosis. However, these markers only cover a small subset of molecular biomarkers linked to the disease. In this study, we aimed to profile a more comprehensive set of molecular biomarkers and explore their temporal association with incident diabetes. METHODS: We performed a targeted analysis of 54 proteins and 171 metabolites and lipoprotein particles measured in three sequential samples spanning up to 11 years of follow-up in 324 individuals with incident diabetes and 359 individuals without diabetes in the Danish Blood Donor Study (DBDS) matched for sex and birth year distribution. We used linear mixed-effects models to identify temporal changes before a diabetes diagnosis, either for any incident diabetes diagnosis or for type 1 and type 2 diabetes mellitus diagnoses specifically. We further performed linear and non-linear feature selection, adding 28 polygenic risk scores to the biomarker pool. We tested the time-to-event prediction gain of the biomarkers with the highest variable importance, compared with selected clinical covariates and plasma glucose. RESULTS: We identified two proteins and 16 metabolites and lipoprotein particles whose levels changed temporally before diabetes diagnosis and for which the estimated marginal means were significant after FDR adjustment. Sixteen of these have not previously been described. Additionally, 75 biomarkers were consistently higher or lower in the years before a diabetes diagnosis. We identified a single temporal biomarker for type 1 diabetes, IL-17A/F, a cytokine that is associated with multiple other autoimmune diseases. Inclusion of 12 biomarkers improved the 10-year prediction of a diabetes diagnosis (i.e. the area under the receiver operating curve increased from 0.79 to 0.84), compared with clinical information and plasma glucose alone. CONCLUSIONS/INTERPRETATION: Systemic molecular changes manifest in plasma several years before a diabetes diagnosis. A particular subset of biomarkers shows distinct, time-dependent patterns, offering potential as predictive markers for diabetes onset. Notably, these biomarkers show shared and distinct patterns between type 1 diabetes and type 2 diabetes. After independent replication, our findings may be used to develop new clinical prediction models.
ABSTRACT
Gene promoter and enhancer sequences are bound by transcription factors and are depleted of methylated CpG sites (cytosines preceding guanines in DNA). The absence of methylated CpGs in these sequences typically correlates with increased gene expression, indicating a regulatory role for methylation. We used nanopore sequencing to determine haplotype-specific methylation rates of 15.3 million CpG units in 7,179 whole-blood genomes. We identified 189,178 methylation depleted sequences where three or more proximal CpGs were unmethylated on at least one haplotype. A total of 77,789 methylation depleted sequences (~41%) associated with 80,503 cis-acting sequence variants, which we termed allele-specific methylation quantitative trait loci (ASM-QTLs). RNA sequencing of 896 samples from the same blood draws used to perform nanopore sequencing showed that the ASM-QTL, that is, DNA sequence variability, drives most of the correlation found between gene expression and CpG methylation. ASM-QTLs were enriched 40.2-fold (95% confidence interval 32.2, 49.9) among sequence variants associating with hematological traits, demonstrating that ASM-QTLs are important functional units in the noncoding genome.
Subject(s)
CpG Islands , DNA Methylation , Quantitative Trait Loci , Humans , Promoter Regions, Genetic , Haplotypes , Alleles , Gene Expression Regulation , Genetic Variation , Nanopore Sequencing/methods , Genome, HumanSubject(s)
Alzheimer Disease , Membrane Glycoproteins , Receptors, Immunologic , Aged , Female , Humans , Male , Alzheimer Disease/ethnology , Alzheimer Disease/genetics , Alzheimer Disease/metabolism , Apolipoproteins E/metabolism , Cholesterol, LDL/metabolism , Clusterin/metabolism , Genetic Predisposition to Disease/genetics , Homozygote , Membrane Glycoproteins/genetics , Membrane Glycoproteins/metabolism , Mutation, Missense , Receptors, Immunologic/genetics , Receptors, Immunologic/metabolismABSTRACT
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.
Subject(s)
Essential Tremor , Genetic Predisposition to Disease , Genome-Wide Association Study , Essential Tremor/genetics , Humans , Polymorphism, Single Nucleotide , Genetic LociABSTRACT
Genomic prediction of antipsychotic dose and polypharmacy has been difficult, mainly due to limited access to large cohorts with genetic and drug prescription data. In this proof of principle study, we investigated if genetic liability for schizophrenia is associated with high dose requirements of antipsychotics and antipsychotic polypharmacy, using real-world registry and biobank data from five independent Nordic cohorts of a total of N = 21,572 individuals with psychotic disorders (schizophrenia, bipolar disorder, and other psychosis). Within regression models, a polygenic risk score (PRS) for schizophrenia was studied in relation to standardized antipsychotic dose as well as antipsychotic polypharmacy, defined based on longitudinal prescription registry data as well as health records and self-reported data. Meta-analyses across the five cohorts showed that PRS for schizophrenia was significantly positively associated with prescribed (standardized) antipsychotic dose (beta(SE) = 0.0435(0.009), p = 0.0006) and antipsychotic polypharmacy defined as taking ≥2 antipsychotics (OR = 1.10, CI = 1.05-1.21, p = 0.0073). The direction of effect was similar in all five independent cohorts. These findings indicate that genotypes may aid clinically relevant decisions on individual patients´ antipsychotic treatment. Further, the findings illustrate how real-world data have the potential to generate results needed for future precision medicine approaches in psychiatry.
Subject(s)
Antipsychotic Agents , Biological Specimen Banks , Multifactorial Inheritance , Polypharmacy , Registries , Schizophrenia , Humans , Antipsychotic Agents/administration & dosage , Antipsychotic Agents/therapeutic use , Male , Female , Schizophrenia/drug therapy , Schizophrenia/genetics , Middle Aged , Adult , Psychotic Disorders/drug therapy , Psychotic Disorders/genetics , Cohort Studies , AgedABSTRACT
Precision medicine has the ambition to improve treatment response and clinical outcomes through patient stratification and holds great potential for the treatment of mental disorders. However, several important factors are needed to transform current practice into a precision psychiatry framework. Most important are 1) the generation of accessible large real-world training and test data including genomic data integrated from multiple sources, 2) the development and validation of advanced analytical tools for stratification and prediction, and 3) the development of clinically useful management platforms for patient monitoring that can be integrated into health care systems in real-life settings. This narrative review summarizes strategies for obtaining the key elements-well-powered samples from large biobanks integrated with electronic health records and health registry data using novel artificial intelligence algorithms-to predict outcomes in severe mental disorders and translate these models into clinical management and treatment approaches. Key elements are massive mental health data and novel artificial intelligence algorithms. For the clinical translation of these strategies, we discuss a precision medicine platform for improved management of mental disorders. We use cases to illustrate how precision medicine interventions could be brought into psychiatry to improve the clinical outcomes of mental disorders.
ABSTRACT
BACKGROUND: In 2021, the American College of Medical Genetics and Genomics (ACMG) recommended reporting actionable genotypes in 73 genes associated with diseases for which preventive or therapeutic measures are available. Evaluations of the association of actionable genotypes in these genes with life span are currently lacking. METHODS: We assessed the prevalence of coding and splice variants in genes on the ACMG Secondary Findings, version 3.0 (ACMG SF v3.0), list in the genomes of 57,933 Icelanders. We assigned pathogenicity to all reviewed variants using reported evidence in the ClinVar database, the frequency of variants, and their associations with disease to create a manually curated set of actionable genotypes (variants). We assessed the relationship between these genotypes and life span and further examined the specific causes of death among carriers. RESULTS: Through manual curation of 4405 sequence variants in the ACMG SF v3.0 genes, we identified 235 actionable genotypes in 53 genes. Of the 57,933 participants, 2306 (4.0%) carried at least one actionable genotype. We found shorter median survival among persons carrying actionable genotypes than among noncarriers. Specifically, we found that carrying an actionable genotype in a cancer gene was associated with survival that was 3 years shorter than that among noncarriers, with causes of death among carriers attributed primarily to cancer-related conditions. Furthermore, we found evidence of association between carrying an actionable genotype in certain genes in the cardiovascular disease group and a reduced life span. CONCLUSIONS: On the basis of the ACMG SF v3.0 guidelines, we found that approximately 1 in 25 Icelanders carried an actionable genotype and that carrying such a genotype was associated with a reduced life span. (Funded by deCODE Genetics-Amgen.).
Subject(s)
Disease , Genomics , Longevity , Humans , Alleles , Genetic Testing , Genetic Variation , Genotype , Iceland/epidemiology , Longevity/genetics , Disease/genetics , Cardiovascular Diseases/genetics , Neoplasms/geneticsABSTRACT
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.
ABSTRACT
Background: The anti-seizure medication vigabatrin (VGB) is effective for controlling seizures, especially infantile spasms. However, use is limited by VGB-associated visual field loss (VAVFL). The mechanisms by which VGB causes VAVFL remains unknown. Average peripapillary retinal nerve fibre layer (ppRNFL) thickness correlates with the degree of visual field loss (measured by mean radial degrees). Duration of VGB exposure, maximum daily VGB dose, and male sex are associated with ppRNFL thinning. Here we test the hypothesis that common genetic variation is a predictor of ppRNFL thinning in VGB exposed individuals. Identifying pharmacogenomic predictors of ppRNFL thinning in VGB exposed individuals could potentially enable safe prescribing of VGB and broader use of a highly effective drug. Methods: Optical coherence topography (OCT) and GWAS data were processed from VGB-exposed individuals (n = 71) recruited through the EpiPGX Consortium. We conducted quantitative GWAS analyses for the following OCT measurements: (1) average ppRNFL, (2) inferior quadrant, (3) nasal quadrant, (4) superior quadrant, (5) temporal quadrant, (6) inferior nasal sector, (7) nasal inferior sector, (8) superior nasal sector, and (9) nasal superior sector. Using the summary statistics from the GWAS analyses we conducted gene-based testing using VEGAS2. We conducted nine different PRS analyses using the OCT measurements. To determine if VGB-exposed individuals were predisposed to having a thinner RNFL, we calculated their polygenic burden for retinal thickness. PRS alleles for retinal thickness were calculated using published summary statistics from a large-scale GWAS of inner retinal morphology using the OCT images of UK Biobank participants. Results: The GWAS analyses did not identify a significant association after correction for multiple testing. Similarly, the gene-based and PRS analyses did not reveal a significant association that survived multiple testing. Conclusion: We set out to identify common genetic predictors for VGB induced ppRNFL thinning. Results suggest that large-effect common genetic predictors are unlikely to exist for ppRNFL thinning (as a marker of VAVFL). Sample size was a limitation of this study. However, further recruitment is a challenge as VGB is rarely used today because of this adverse reaction. Rare variants may be predictors of this adverse drug reaction and were not studied here.
ABSTRACT
High-throughput proteomics platforms measuring thousands of proteins in plasma combined with genomic and phenotypic information have the power to bridge the gap between the genome and diseases. Here we performed association studies of Olink Explore 3072 data generated by the UK Biobank Pharma Proteomics Project1 on plasma samples from more than 50,000 UK Biobank participants with phenotypic and genotypic data, stratifying on British or Irish, African and South Asian ancestries. We compared the results with those of a SomaScan v4 study on plasma from 36,000 Icelandic people2, for 1,514 of whom Olink data were also available. We found modest correlation between the two platforms. Although cis protein quantitative trait loci were detected for a similar absolute number of assays on the two platforms (2,101 on Olink versus 2,120 on SomaScan), the proportion of assays with such supporting evidence for assay performance was higher on the Olink platform (72% versus 43%). A considerable number of proteins had genomic associations that differed between the platforms. We provide examples where differences between platforms may influence conclusions drawn from the integration of protein levels with the study of diseases. We demonstrate how leveraging the diverse ancestries of participants in the UK Biobank helps to detect novel associations and refine genomic location. Our results show the value of the information provided by the two most commonly used high-throughput proteomics platforms and demonstrate the differences between them that at times provides useful complementarity.
Subject(s)
Blood Proteins , Disease Susceptibility , Genomics , Genotype , Phenotype , Proteomics , Humans , Africa/ethnology , Asia, Southern/ethnology , Biological Specimen Banks , Blood Proteins/analysis , Blood Proteins/genetics , Datasets as Topic , Genome, Human/genetics , Iceland/ethnology , Ireland/ethnology , Plasma/chemistry , Proteome/analysis , Proteome/genetics , Proteomics/methods , Quantitative Trait Loci , United KingdomABSTRACT
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.
Subject(s)
Attention Deficit Disorder with Hyperactivity , Bipolar Disorder , Schizophrenia , Male , Female , Humans , Genome-Wide Association Study , Depression , Bipolar Disorder/epidemiology , Bipolar Disorder/genetics , Schizophrenia/epidemiology , Schizophrenia/genetics , Attention Deficit Disorder with Hyperactivity/epidemiology , Attention Deficit Disorder with Hyperactivity/genetics , Polymorphism, Single Nucleotide/genetics , Genetic Predisposition to DiseaseABSTRACT
The genetic basis of the human vocal system is largely unknown, as are the sequence variants that give rise to individual differences in voice and speech. Here, we couple data on diversity in the sequence of the genome with voice and vowel acoustics in speech recordings from 12,901 Icelanders. We show how voice pitch and vowel acoustics vary across the life span and correlate with anthropometric, physiological, and cognitive traits. We found that voice pitch and vowel acoustics have a heritable component and discovered correlated common variants in ABCC9 that associate with voice pitch. The ABCC9 variants also associate with adrenal gene expression and cardiovascular traits. By showing that voice and vowel acoustics are influenced by genetics, we have taken important steps toward understanding the genetics and evolution of the human vocal system.
Subject(s)
Speech Acoustics , Voice , Humans , Speech/physiology , AcousticsABSTRACT
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.
Subject(s)
Genome , Polymorphism, Single Nucleotide , Female , Humans , Polymorphism, Single Nucleotide/genetics , Whole Genome Sequencing/methods , Genomics , DenmarkABSTRACT
Tourette Syndrome (TS) is a complex neurodevelopmental disorder characterized by vocal and motor tics lasting more than a year. It is highly polygenic in nature with both rare and common previously associated variants. Epidemiological studies have shown TS to be correlated with other phenotypes, but large-scale phenome wide analyses in biobank level data have not been performed to date. In this study, we used the summary statistics from the latest meta-analysis of TS to calculate the polygenic risk score (PRS) of individuals in the UK Biobank data and applied a Phenome Wide Association Study (PheWAS) approach to determine the association of disease risk with a wide range of phenotypes. A total of 57 traits were found to be significantly associated with TS polygenic risk, including multiple psychosocial factors and mental health conditions such as anxiety disorder and depression. Additional associations were observed with complex non-psychiatric disorders such as Type 2 diabetes, heart palpitations, and respiratory conditions. Cross-disorder comparisons of phenotypic associations with genetic risk for other childhood-onset disorders (e.g.: attention deficit hyperactivity disorder [ADHD], autism spectrum disorder [ASD], and obsessive-compulsive disorder [OCD]) indicated an overlap in associations between TS and these disorders. ADHD and ASD had a similar direction of effect with TS while OCD had an opposite direction of effect for all traits except mental health factors. Sex-specific PheWAS analysis identified differences in the associations with TS genetic risk between males and females. Type 2 diabetes and heart palpitations were significantly associated with TS risk in males but not in females, whereas diseases of the respiratory system were associated with TS risk in females but not in males. This analysis provides further evidence of shared genetic and phenotypic architecture of different complex disorders.
Subject(s)
Attention Deficit Disorder with Hyperactivity , Autism Spectrum Disorder , Diabetes Mellitus, Type 2 , Tourette Syndrome , Male , Female , Humans , Tourette Syndrome/genetics , Autism Spectrum Disorder/genetics , Attention Deficit Disorder with Hyperactivity/genetics , Risk FactorsABSTRACT
BACKGROUND: Tourette syndrome (TS) is a childhood-onset neurodevelopmental disorder of complex genetic architecture and is characterized by multiple motor tics and at least one vocal tic persisting for more than 1 year. METHODS: We performed a genome-wide meta-analysis integrating a novel TS cohort with previously published data, resulting in a sample size of 6133 individuals with TS and 13,565 ancestry-matched control participants. RESULTS: We identified a genome-wide significant locus on chromosome 5q15. Integration of expression quantitative trait locus, Hi-C (high-throughput chromosome conformation capture), and genome-wide association study data implicated the NR2F1 gene and associated long noncoding RNAs within the 5q15 locus. Heritability partitioning identified statistically significant enrichment in brain tissue histone marks, while polygenic risk scoring of brain volume data identified statistically significant associations with right and left thalamus volumes and right putamen volume. CONCLUSIONS: Our work presents novel insights into the neurobiology of TS, thereby opening up new directions for future studies.