Identifying canonical and replicable multi-scale intrinsic connectivity networks in 100k+ resting-state fMRI datasets.

Despite the known benefits of data-driven approaches, the lack of approaches for identifying functional neuroimaging patterns that capture both individual variations and inter-subject correspondence limits the clinical utility of rsfMRI and its application to single-subject analyses. Here, using rsfMRI data from over 100k individuals across private and public datasets, we identify replicable multi-spatial-scale canonical intrinsic connectivity network (ICN) templates via the use of multi-model-order independent component analysis (ICA). We also study the feasibility of estimating subject-specific ICNs via spatially constrained ICA. The results show that the subject-level ICN estimations vary as a function of the ICN itself, the data length, and the spatial resolution. In general, large-scale ICNs require less data to achieve specific levels of (within- and between-subject) spatial similarity with their templates. Importantly, increasing data length can reduce an ICN's subject-level specificity, suggesting longer scans may not always be desirable. We also find a positive linear relationship between data length and spatial smoothness (possibly due to averaging over intrinsic dynamics), suggesting studies examining optimized data length should consider spatial smoothness. Finally, consistency in spatial similarity between ICNs estimated using the full data and subsets across different data lengths suggests lower within-subject spatial similarity in shorter data is not wholly defined by lower reliability in ICN estimates, but may be an indication of meaningful brain dynamics which average out as data length increases.

A method for estimating and characterizing explicitly nonlinear dynamic functional network connectivity in resting-state fMRI data.

The past 10 years have seen an explosion of approaches that focus on the study of time-resolved change in functional connectivity (FC). FC characterization among networks at a whole-brain level is frequently termed functional network connectivity (FNC). Time-resolved or dynamic functional network connectivity (dFNC) focuses on the estimation of transient, recurring, whole-brain patterns of FNC. While most approaches in this area have attempted to capture dynamic linear correlation, we are particularly interested in whether explicitly nonlinear relationships, above and beyond linear, are present and contain unique information. This study thus proposes an approach to assess explicitly nonlinear dynamic functional network connectivity (EN dFNC) derived from the relationship among independent component analysis time courses. Linear relationships were removed at each time point to evaluate, typically ignored, explicitly nonlinear dFNC using normalized mutual information (NMI). Simulations showed the proposed method estimated explicitly nonlinearity over time, even within relatively short windows of data. We then, applied our approach on 151 schizophrenia patients, and 163 healthy controls fMRI data and found three unique, highly structured, mostly long-range, functional states that also showed significant group differences. In particular, explicitly nonlinear relationships tend to be more widespread than linear ones. Results also highlighted a state with long range connections to the visual domain, which were significantly reduced in schizophrenia. Overall, this work suggests that quantifying EN dFNC may provide a complementary and potentially valuable tool for studying brain function by exposing relevant variation that is typically ignored.

Moving beyond the 'CAP' of the Iceberg: Intrinsic connectivity networks in fMRI are continuously engaging and overlapping.

Resting-state functional magnetic resonance imaging is currently the mainstay of functional neuroimaging and has allowed researchers to identify intrinsic connectivity networks (aka functional networks) at different spatial scales. However, little is known about the temporal profiles of these networks and whether it is best to model them as continuous phenomena in both space and time or, rather, as a set of temporally discrete events. Both categories have been supported by series of studies with promising findings. However, a critical question is whether focusing only on time points presumed to contain isolated neural events and disregarding the rest of the data is missing important information, potentially leading to misleading conclusions. In this work, we argue that brain networks identified within the spontaneous blood oxygenation level-dependent (BOLD) signal are not limited to temporally sparse burst moments and that these event present time points (EPTs) contain valuable but incomplete information about the underlying functional patterns. We focus on the default mode and show evidence that is consistent with its continuous presence in the BOLD signal, including during the event absent time points (EATs), i.e., time points that exhibit minimum activity and are the least likely to contain an event. Moreover, our findings suggest that EPTs may not contain all the available information about their corresponding networks. We observe distinct default mode connectivity patterns obtained from all time points (AllTPs), EPTs, and EATs. We show evidence of robust relationships with schizophrenia symptoms that are both common and unique to each of the sets of time points (AllTPs, EPTs, EATs), likely related to transient patterns of connectivity. Together, these findings indicate the importance of leveraging the full temporal data in functional studies, including those using event-detection approaches.

Resting-state thalamic dysconnectivity in schizophrenia and relationships with symptoms.

BACKGROUND: Schizophrenia (SZ) is a severe neuropsychiatric disorder associated with disrupted connectivity within the thalamic-cortico-cerebellar network. Resting-state functional connectivity studies have reported thalamic hypoconnectivity with the cerebellum and prefrontal cortex as well as thalamic hyperconnectivity with sensory cortical regions in SZ patients compared with healthy comparison participants (HCs). However, fundamental questions remain regarding the clinical significance of these connectivity abnormalities. METHOD: Resting state seed-based functional connectivity was used to investigate thalamus to whole brain connectivity using multi-site data including 183 SZ patients and 178 matched HCs. Statistical significance was based on a voxel-level FWE-corrected height threshold of p < 0.001. The relationships between positive and negative symptoms of SZ and regions of the brain demonstrating group differences in thalamic connectivity were examined. RESULTS: HC and SZ participants both demonstrated widespread positive connectivity between the thalamus and cortical regions. Compared with HCs, SZ patients had reduced thalamic connectivity with bilateral cerebellum and anterior cingulate cortex. In contrast, SZ patients had greater thalamic connectivity with multiple sensory-motor regions, including bilateral pre- and post-central gyrus, middle/inferior occipital gyrus, and middle/superior temporal gyrus. Thalamus to middle temporal gyrus connectivity was positively correlated with hallucinations and delusions, while thalamus to cerebellar connectivity was negatively correlated with delusions and bizarre behavior. CONCLUSIONS: Thalamic hyperconnectivity with sensory regions and hypoconnectivity with cerebellar regions in combination with their relationship to clinical features of SZ suggest that thalamic dysconnectivity may be a core neurobiological feature of SZ that underpins positive symptoms.

Decreased hemispheric connectivity and decreased intra- and inter- hemisphere asymmetry of resting state functional network connectivity in schizophrenia.

Many studies have shown that schizophrenia patients have aberrant functional network connectivity (FNC) among brain regions, suggesting schizophrenia manifests with significantly diminished (in majority of the cases) connectivity. Schizophrenia is also associated with a lack of hemispheric lateralization. Hoptman et al. (2012) reported lower inter-hemispheric connectivity in schizophrenia patients compared to controls using voxel-mirrored homotopic connectivity. In this study, we merge these two points of views together using a group independent component analysis (gICA)-based approach to generate hemisphere-specific timecourses and calculate intra-hemisphere and inter-hemisphere FNC on a resting state fMRI dataset consisting of age- and gender-balanced 151 schizophrenia patients and 163 healthy controls. We analyzed the group differences between patients and healthy controls in each type of FNC measures along with age and gender effects. The results reveal that FNC in schizophrenia patients shows less hemispheric asymmetry compared to that of the healthy controls. We also found a decrease in connectivity in all FNC types such as intra-left (L_FNC), intra-right (R_FNC) and inter-hemisphere (Inter_FNC) in the schizophrenia patients relative to healthy controls, but general patterns of connectivity were preserved in patients. Analyses of age and gender effects yielded results similar to those reported in whole brain FNC studies.

Prefrontal cortical thinning links to negative symptoms in schizophrenia via the ENIGMA consortium.

BACKGROUND: Our understanding of the complex relationship between schizophrenia symptomatology and etiological factors can be improved by studying brain-based correlates of schizophrenia. Research showed that impairments in value processing and executive functioning, which have been associated with prefrontal brain areas [particularly the medial orbitofrontal cortex (MOFC)], are linked to negative symptoms. Here we tested the hypothesis that MOFC thickness is associated with negative symptom severity. METHODS: This study included 1985 individuals with schizophrenia from 17 research groups around the world contributing to the ENIGMA Schizophrenia Working Group. Cortical thickness values were obtained from T1-weighted structural brain scans using FreeSurfer. A meta-analysis across sites was conducted over effect sizes from a model predicting cortical thickness by negative symptom score (harmonized Scale for the Assessment of Negative Symptoms or Positive and Negative Syndrome Scale scores). RESULTS: Meta-analytical results showed that left, but not right, MOFC thickness was significantly associated with negative symptom severity (ß std = -0.075; p = 0.019) after accounting for age, gender, and site. This effect remained significant (p = 0.036) in a model including overall illness severity. Covarying for duration of illness, age of onset, antipsychotic medication or handedness weakened the association of negative symptoms with left MOFC thickness. As part of a secondary analysis including 10 other prefrontal regions further associations in the left lateral orbitofrontal gyrus and pars opercularis emerged. CONCLUSIONS: Using an unusually large cohort and a meta-analytical approach, our findings point towards a link between prefrontal thinning and negative symptom severity in schizophrenia. This finding provides further insight into the relationship between structural brain abnormalities and negative symptoms in schizophrenia.

Newborn insula gray matter volume is prospectively associated with early life adiposity gain.

BACKGROUND: The importance of energy homeostasis brain circuitry in the context of obesity is well established, however, the developmental ontogeny of this circuitry in humans is currently unknown. Here, we investigate the prospective association between newborn gray matter (GM) volume in the insula, a key brain region underlying energy homeostasis, and change in percent body fat accrual over the first six months of postnatal life, an outcome that represents among the most reliable infant predictors of childhood obesity risk. METHODS: A total of 52 infants (29 male, 23 female, gestational age at birth=39(1.5) weeks) were assessed using structural MRI shortly after birth (postnatal age at MRI scan=25.9(12.2) days), and serial Dual X-Ray Absorptiometry shortly after birth (postnatal age at DXA scan 1=24.6(11.4) days) and at six months of age (postnatal age at DXA scan 2=26.7(3.3) weeks). RESULTS: Insula GM volume was inversely associated with change in percent body fat from birth to six-months postnatal age and accounted for 19% of its variance (ß=-3.6%/S.D., P=0.001). This association was driven by the central-posterior portion of the insula, a region of particular importance for gustation and interoception. The direction of this effect is in concordance with observations in adults, and the results remained statistically significant after adjusting for relevant covariates and potential confounding variables. CONCLUSIONS: Altogether, these findings suggest an underlying neural basis of childhood obesity that precedes the influence of the postnatal environment. The identification of plausible brain-related biomarkers of childhood obesity risk that predate the influence of the postnatal obesogenic environment may contribute to an improved understanding of propensity for obesity, early identification of at-risk individuals, and intervention targets for primary prevention.

Positive symptoms associate with cortical thinning in the superior temporal gyrus via the ENIGMA Schizophrenia consortium.

OBJECTIVE: Based on the role of the superior temporal gyrus (STG) in auditory processing, language comprehension and self-monitoring, this study aimed to investigate the relationship between STG cortical thickness and positive symptom severity in schizophrenia. METHOD: This prospective meta-analysis includes data from 1987 individuals with schizophrenia collected at seventeen centres around the world that contribute to the ENIGMA Schizophrenia Working Group. STG thickness measures were extracted from T1-weighted brain scans using FreeSurfer. The study performed a meta-analysis of effect sizes across sites generated by a model predicting left or right STG thickness with a positive symptom severity score (harmonized SAPS or PANSS-positive scores), while controlling for age, sex and site. Secondary models investigated relationships between antipsychotic medication, duration of illness, overall illness severity, handedness and STG thickness. RESULTS: Positive symptom severity was negatively related to STG thickness in both hemispheres (left: ßstd = -0.052; P = 0.021; right: ßstd = -0.073; P = 0.001) when statistically controlling for age, sex and site. This effect remained stable in models including duration of illness, antipsychotic medication or handedness. CONCLUSION: Our findings further underline the important role of the STG in hallmark symptoms in schizophrenia. These findings can assist in advancing insight into symptom-relevant pathophysiological mechanisms in schizophrenia.

The role of the ITIH3 rs2535629 variant in antipsychotic response.

INTRODUCTION: There is mounting evidence that schizophrenia risk variants influence response to antipsychotic medication. Common single nucleotide polymorphisms (SNPs) in or near the inter-alpha-trypsin inhibitor heavy chain H3 (ITIH3) gene have been repeatedly associated with schizophrenia and related psychiatric disorders in genome-wide association studies. Here, we provide the first study to assess the relevance of the ITIH3 rs2535629 SNP in response to antipsychotic medication. METHODS: The rs2535629 SNP was genotyped in N=256 patients receiving various antipsychotics for up to 26weeks. Treatment response was assessed using the Brief Psychiatric Rating Scale (BPRS) including its positive and negative subscales. Follow-up analyses were performed after stratifying for ethnicity and medication. RESULTS: We found significant association of rs2535629 with improvement of negative symptoms in patients of European ancestry after six months of clozapine treatment (F1,87=8.8, pcorr=0.032). Patients homozygous for the minor A-allele showed the best improvement of negative BPRS scores. However, we observed no association between rs2535629 and changes in total BPRS score in the entire sample or the clozapine-treated subgroup. DISCUSSION: Although there was no association of genotype with overall changes in BPRS scores, the greater improvement of negative symptoms in minor allele carriers indicates that rs2535629 may help to identify a subset of schizophrenia patients with better treatment response to clozapine. Therefore, our findings provide the first suggestive evidence that rs2535629 is relevant in antipsychotic response.

Subcortical brain volume abnormalities in 2028 individuals with schizophrenia and 2540 healthy controls via the ENIGMA consortium.

The profile of brain structural abnormalities in schizophrenia is still not fully understood, despite decades of research using brain scans. To validate a prospective meta-analysis approach to analyzing multicenter neuroimaging data, we analyzed brain MRI scans from 2028 schizophrenia patients and 2540 healthy controls, assessed with standardized methods at 15 centers worldwide. We identified subcortical brain volumes that differentiated patients from controls, and ranked them according to their effect sizes. Compared with healthy controls, patients with schizophrenia had smaller hippocampus (Cohen's d=-0.46), amygdala (d=-0.31), thalamus (d=-0.31), accumbens (d=-0.25) and intracranial volumes (d=-0.12), as well as larger pallidum (d=0.21) and lateral ventricle volumes (d=0.37). Putamen and pallidum volume augmentations were positively associated with duration of illness and hippocampal deficits scaled with the proportion of unmedicated patients. Worldwide cooperative analyses of brain imaging data support a profile of subcortical abnormalities in schizophrenia, which is consistent with that based on traditional meta-analytic approaches. This first ENIGMA Schizophrenia Working Group study validates that collaborative data analyses can readily be used across brain phenotypes and disorders and encourages analysis and data sharing efforts to further our understanding of severe mental illness.

Variable telomere length across post-mortem human brain regions and specific reduction in the hippocampus of major depressive disorder.

Stress can be a predisposing factor to psychiatric disorders and has been associated with decreased neurogenesis and reduced hippocampal volume especially in depression. Similarly, in white blood cells chronic psychological stress has been associated with telomere shortening and with mood disorders and schizophrenia (SZ). However, in previous post-mortem brain studies from occipital cortex and cerebellum, no difference in telomere length was observed in depression. We hypothesized that in psychiatric disorders, stress-driven accelerated cellular aging can be observed in brain regions particularly sensitive to stress. Telomere length was measured by quantitative-PCR in five brain regions (dorsolateral prefrontal cortex, hippocampus (HIPP), amygdala, nucleus accumbens and substantia nigra (SN)) in major depressive disorder (MDD), bipolar disorder, SZ and normal control subjects (N = 40, 10 subjects per group). We observed significant differences in telomere length across brain regions suggesting variable levels of cell aging, with SN and HIPP having the longest telomeres and the dorsolateral prefrontal cortex the shortest. A significant decrease (P < 0.02) in telomere length was observed specifically in the HIPP of MDD subjects even after controlling for age. In the HIPP of MDD subjects, several genes involved in neuroprotection and in stress response (FKBP5, CRH) showed altered levels of mRNA. Our results suggest the presence of hippocampal stress-mediated accelerated cellular aging in depression. Further studies are needed to investigate the cellular specificity of these findings.

Huntington's disease cerebrospinal fluid seeds aggregation of mutant huntingtin.

Huntington's disease (HD), a progressive neurodegenerative disease, is caused by an expanded CAG triplet repeat producing a mutant huntingtin protein (mHTT) with a polyglutamine-repeat expansion. Onset of symptoms in mutant huntingtin gene-carrying individuals remains unpredictable. We report that synthetic polyglutamine oligomers and cerebrospinal fluid (CSF) from BACHD transgenic rats and from human HD subjects can seed mutant huntingtin aggregation in a cell model and its cell lysate. Our studies demonstrate that seeding requires the mutant huntingtin template and may reflect an underlying prion-like protein propagation mechanism. Light and cryo-electron microscopy show that synthetic seeds nucleate and enhance mutant huntingtin aggregation. This seeding assay distinguishes HD subjects from healthy and non-HD dementia controls without overlap (blinded samples). Ultimately, this seeding property in HD patient CSF may form the basis of a molecular biomarker assay to monitor HD and evaluate therapies that target mHTT.

Dynamic functional connectivity analysis reveals transient states of dysconnectivity in schizophrenia.

Schizophrenia is a psychotic disorder characterized by functional dysconnectivity or abnormal integration between distant brain regions. Recent functional imaging studies have implicated large-scale thalamo-cortical connectivity as being disrupted in patients. However, observed connectivity differences in schizophrenia have been inconsistent between studies, with reports of hyperconnectivity and hypoconnectivity between the same brain regions. Using resting state eyes-closed functional imaging and independent component analysis on a multi-site data that included 151 schizophrenia patients and 163 age- and gender matched healthy controls, we decomposed the functional brain data into 100 components and identified 47 as functionally relevant intrinsic connectivity networks. We subsequently evaluated group differences in functional network connectivity, both in a static sense, computed as the pairwise Pearson correlations between the full network time courses (5.4 minutes in length), and a dynamic sense, computed using sliding windows (44 s in length) and k-means clustering to characterize five discrete functional connectivity states. Static connectivity analysis revealed that compared to healthy controls, patients show significantly stronger connectivity, i.e., hyperconnectivity, between the thalamus and sensory networks (auditory, motor and visual), as well as reduced connectivity (hypoconnectivity) between sensory networks from all modalities. Dynamic analysis suggests that (1), on average, schizophrenia patients spend much less time than healthy controls in states typified by strong, large-scale connectivity, and (2), that abnormal connectivity patterns are more pronounced during these connectivity states. In particular, states exhibiting cortical-subcortical antagonism (anti-correlations) and strong positive connectivity between sensory networks are those that show the group differences of thalamic hyperconnectivity and sensory hypoconnectivity. Group differences are weak or absent during other connectivity states. Dynamic analysis also revealed hypoconnectivity between the putamen and sensory networks during the same states of thalamic hyperconnectivity; notably, this finding cannot be observed in the static connectivity analysis. Finally, in post-hoc analyses we observed that the relationships between sub-cortical low frequency power and connectivity with sensory networks is altered in patients, suggesting different functional interactions between sub-cortical nuclei and sensorimotor cortex during specific connectivity states. While important differences between patients with schizophrenia and healthy controls have been identified, one should interpret the results with caution given the history of medication in patients. Taken together, our results support and expand current knowledge regarding dysconnectivity in schizophrenia, and strongly advocate the use of dynamic analyses to better account for and understand functional connectivity differences.

Absence of evidence for bornavirus infection in schizophrenia, bipolar disorder and major depressive disorder.

In 1983, reports of antibodies in subjects with major depressive disorder (MDD) to an as-yet uncharacterized infectious agent associated with meningoencephalitis in horses and sheep led to molecular cloning of the genome of a novel, negative-stranded neurotropic virus, Borna disease virus (BDV). This advance has enabled the development of new diagnostic assays, including in situ hybridization, PCR and serology based on recombinant proteins. Since these assays were first implemented in 1990, more than 80 studies have reported an association between BDV and a wide range of human illnesses that include MDD, bipolar disorder (BD), schizophrenia (SZ), anxiety disorder, chronic fatigue syndrome, multiple sclerosis, amyotrophic lateral sclerosis, dementia and glioblastoma multiforme. However, to date there has been no blinded case-control study of the epidemiology of BDV infection. Here, in a United States-based, multi-center, yoked case-control study with standardized methods for clinical assessment and blinded serological and molecular analysis, we report the absence of association of psychiatric illness with antibodies to BDV or with BDV nucleic acids in serially collected serum and white blood cell samples from 396 subjects, a study population comprised of 198 matched pairs of patients and healthy controls (52 SZ/control pairs, 66 BD/control pairs and 80 MDD/control pairs). Our results argue strongly against a role for BDV in the pathogenesis of these psychiatric disorders.

Association study of cannabinoid receptor 1 (CNR1) gene in tardive dyskinesia.

Tardive dyskinesia (TD) is a severe, debilitating movement disorder observed in 25-30% of the patients treated with typical antipsychotics. Cannabinoid receptor 1 (CNR1) activators tend to inhibit movement, an effect prevented by rimonabant and other selective CNR1 antagonists. Furthermore, CNR1 receptor is downregulated in Huntington's disease and upregulated in Parkinson's disease. Twenty tagSNPs spanning the CNR1 gene were analyzed in schizophrenia patients of European ancestry (n=191; 74 with TD). Significant genotypic (P=0.012) and allelic (P=0.012) association was observed with rs806374 (T>C). Carriers of the CC genotype were more likely to be TD positive (CC vs TT+TC, odds ratio=3.4 (1.5-7.8), P=0.003) and had more severe TD (CC vs TT+TC; 9.52±9.2 vs 5.62±6.9, P=0.046). These results indicate a possible role of CNR1 in the development of TD in our patient population. However, these observations are marginal after correcting for multiple testing and need to be replicated in a larger patient population.

Genome-wide association study of CSF biomarkers Abeta1-42, t-tau, and p-tau181p in the ADNI cohort.

OBJECTIVES: CSF levels of Aß1-42, t-tau, and p-tau181p are potential early diagnostic markers for probable Alzheimer disease (AD). The influence of genetic variation on these markers has been investigated for candidate genes but not on a genome-wide basis. We report a genome-wide association study (GWAS) of CSF biomarkers (Aß1-42, t-tau, p-tau181p, p-tau181p/Aß1-42, and t-tau/Aß1-42). METHODS: A total of 374 non-Hispanic Caucasian participants in the Alzheimer's Disease Neuroimaging Initiative cohort with quality-controlled CSF and genotype data were included in this analysis. The main effect of single nucleotide polymorphisms (SNPs) under an additive genetic model was assessed on each of 5 CSF biomarkers. The p values of all SNPs for each CSF biomarker were adjusted for multiple comparisons by the Bonferroni method. We focused on SNPs with corrected p<0.01 (uncorrected p<3.10×10(-8)) and secondarily examined SNPs with uncorrected p values less than 10(-5) to identify potential candidates. RESULTS: Four SNPs in the regions of the APOE, LOC100129500, TOMM40, and EPC2 genes reached genome-wide significance for associations with one or more CSF biomarkers. SNPs in CCDC134, ABCG2, SREBF2, and NFATC4, although not reaching genome-wide significance, were identified as potential candidates. CONCLUSIONS: In addition to known candidate genes, APOE, TOMM40, and one hypothetical gene LOC100129500 partially overlapping APOE; one novel gene, EPC2, and several other interesting genes were associated with CSF biomarkers that are related to AD. These findings, especially the new EPC2 results, require replication in independent cohorts.

Effect of dopamine D3 receptor gene polymorphisms and clozapine treatment response: exploratory analysis of nine polymorphisms and meta-analysis of the Ser9Gly variant.

D2 blockade has been implicated in having a central role in antipsychotic response. However, treatment refractoriness, in spite of complete D2 blockade, as well as the efficacy of clozapine (CLZ) in a portion of this patient population, indicates the involvement of other factors as well. Several lines of evidence suggest a role for D3. Furthermore, an earlier meta-analysis by Jönsson et al. (2003) (n=233) suggested a role for genetic variation in the D3 gene. Relevant to this study, Jönsson et al. found the Ser allele of the D3 serine-to-glycine substitution at amino acid position 9 (Ser9Gly) polymorphism to be associated with worse CLZ response compared with the Gly allele. In this study, we attempt to validate these findings by performing a meta-analysis in a much larger sample (n=758). Eight other variants were also tested in our own sample to explore the possible effect of other regions of the gene. We report a negative but consistent trend across individual studies in our meta-analysis for the DRD3 Ser allele and poor CLZ response. A possible minor role for this single-nucleotide polymorphism cannot be disregarded, as our sample size may have been insufficient. Other DRD3 variants and haplotypes of possible interest were also identified for replication in future studies.

Association study of tardive dyskinesia and five DRD4 polymorphisms in schizophrenia patients.

Tardive dyskinesia (TD) is a side effect of chronic antipsychotic medication exposure. Abnormalities in dopaminergic activity in the nigro-striatal system have been most often suggested to be involved because the agents that cause TD share in common potent antagonism of dopamine D(2) receptors (DRD2). Thus, a number of studies have focused on the association of dopamine system gene polymorphisms and TD, with the most consistent findings being an association between TD and the Ser9Gly polymorphism of the DRD3 gene and the TaqIA site 3' of the DRD2 gene. The DRD4 gene codes for the third member of the D(2)-like dopamine receptor family, and the variable number tandem-repeat polymorphism in exon 3 of DRD4 has been associated with TD. However, other polymorphisms have not been thoroughly examined. In this study, we investigated five polymorphisms spanning the DRD4 gene and their association with TD in our European Caucasian sample (N=171). Although the exon 3 variable number tandem repeat was not associated with TD, haplotypes consisting of four tag polymorphisms were associated with TD in males. This study suggests that DRD4 may be involved in TD in the Caucasian population, although further replication studies are needed.