A pseudolikelihood approach for assessing genetic association in case-control studies with unmeasured population structure

The case-control study design is one of the main tools for detecting associations between genetic markers and diseases. It is well known that population substructure can lead to spurious association between disease status and a genetic marker if the prevalence of disease and the marker allele frequency vary across subpopulations. In this paper, we propose a novel statistical method to estimate the association in case-control studies with unmeasured population substructure. The proposed method takes two steps. First, the information on genomic markers and disease status is used to infer the population substructure; second, the association between the disease and the test marker adjusting for the population substructure is modeled and estimated parametrically through polytomous logistic regression. The performance of the proposed method, relative to the existing methods, on bias, coverage probability and computational time, is assessed through simulations. The method is applied to an end-stage renal disease study in African Americans population.

Trichotillomania comorbidity in a sample enriched for familial obsessive-compulsive disorder

BACKGROUND

This study addresses the strength of associations between trichotillomania (TTM) and other DSM-IV Axis I conditions in a large sample (n = 2606) enriched for familial obsessive-compulsive disorder (OCD), to inform TTM classification.

METHODS

We identified participants with TTM in the Johns Hopkins OCD Family Study (153 families) and the OCD Collaborative Genetics Study, a six-site genetic linkage study of OCD (487 families). We used logistic regression (with generalized estimating equations) to assess the strength of associations between TTM and other DSM-IV disorders.

RESULTS

TTM had excess comorbidity with a number of conditions from different DSM-IV chapters, including tic disorders, alcohol dependence, mood disorders, anxiety disorders, impulse-control disorders, and bulimia nervosa. However, association strengths (odds ratios) were highest for kleptomania (6.6), pyromania (5.8), OCD (5.6), skin picking disorder (4.4), bulimia nervosa (3.5), and pathological nail biting (3.4).

CONCLUSIONS

TTM is comorbid with a number of psychiatric conditions besides OCD, and it is strongly associated with other conditions involving impaired impulse control. Though DSM-5 includes TTM as an OCD-related disorder, its comorbidity pattern also emphasizes the impulsive, appetitive aspects of this condition that may be relevant to classification.

Estimation of Genetic Correlation via Linkage Disequilibrium Score Regression and Genomic Restricted Maximum Likelihood

Genetic correlation is a key population parameter that describes the shared genetic architecture of complex traits and diseases. It can be estimated by current state-of-art methods, i.e., linkage disequilibrium score regression (LDSC) and genomic restricted maximum likelihood (GREML). The massively reduced computing burden of LDSC compared to GREML makes it an attractive tool, although the accuracy (i.e., magnitude of standard errors) of LDSC estimates has not been thoroughly studied. In simulation, we show that the accuracy of GREML is generally higher than that of LDSC. When there is genetic heterogeneity between the actual sample and reference data from which LD scores are estimated, the accuracy of LDSC decreases further. In real data analyses estimating the genetic correlation between schizophrenia (SCZ) and body mass index, we show that GREML estimates based on ∼150,000 individuals give a higher accuracy than LDSC estimates based on ∼400,000 individuals (from combined meta-data). A GREML genomic partitioning analysis reveals that the genetic correlation between SCZ and height is significantly negative for regulatory regions, which whole genome or LDSC approach has less power to detect. We conclude that LDSC estimates should be carefully interpreted as there can be uncertainty about homogeneity among combined meta-datasets. We suggest that any interesting findings from massive LDSC analysis for a large number of complex traits should be followed up, where possible, with more detailed analyses with GREML methods, even if sample sizes are lesser.

Oral Clefts, Maternal Smoking, and TGFA: A Meta-Analysis of Gene-Environment Interaction

Objective

A meta-analysis was performed to examine the association among maternal cigarette smoking, infant genotype at the Taq1 site in the transforming growth factor α (TGFA) locus, and risk of nonsyndromic oral clefts, both cleft palate (CP) and cleft lip with or without cleft palate (CL/P).

Design

Five published case-control studies were included in the meta-analyis. Pooled Mantel-Haenszel odds ratios (OR) and 95% confidence intervals (CIs) were computed. Gene-environment interaction was also assessed by using the pooled data in a case-only analysis and polytomous logistic regression.

Results

Among nonsmoking mothers, there was no evidence of any increased risk for CP if the infant carried the TGFA Taq1 C2 allele. If the mother reported smoking, however, there was an overall increased risk for CP if the infant carried the C2 allele (ORsmokers = 1.95; 95% CI = 1.22 to 3.10). TGFA genotype did not increase risk to CL/P, regardless of maternal smoking status. Polytomous logistic regression revealed a significant overall smoking effect for CL/P (OR = 1.64, 95% CI = 1.33 to 2.02) and CP (OR = 1.42, 95% CI = 1.06 to 1.90).

Conclusions

While maternal smoking was a consistent risk factor for both CL/P and CP across all studies, the suggestive evidence for gene-environment interaction between the infant's genotype at the Taq1 marker in TGFA and maternal smoking was limited to CP. Furthermore, evidence for such gene-environment interaction was strongest in a case-control study drawn from a birth defect registry where infants with non-cleft defects served as controls.

Subtypes of physical frailty: Latent class analysis and associations with clinical characteristics and outcomes

Frailty is a well-recognized geriatric syndrome with various definitions and conceptual frameworks. This study aimed to use latent class analysis to discover potential subtypes of pre-frail and frail older people. Data from the I-Lan Longitudinal Aging Study (ILAS), a community-based cohort study was used for analysis. Latent class analysis was applied to characterize classes or subgroups with different frailty phenotypes among ILAS participants targeting older adults aged 65 and above, capable of completing a 6-meter walk, without severe major or life threatening diseases, and not institutionalized. Latent class analysis identified three distinct subgroups with different frailty phenotypes: non-mobility-type (weight loss and exhaustion), mobility-type frailty (slowness and weakness), and low physical activity. Comparing these groups with the robust group, people with mobility-type frailty had poorer body composition, worse bone health, poorer cognitive function, lower survival (hazard ratio: 6.82, p = 0.019), and poorer overall health outcomes (hazard ratio: 1.67, p = 0.040). People in the non-mobility-type group had poorer bone health and more metabolic serum abnormalities. In conclusion, mobility-type frailty was a better predictor of adverse outcomes. However, further investigation is needed to evaluate how these phenotypic subgroups may help in predicting prognosis or in developing interventions.

Testing homogeneity in semiparametric mixture case-control models

Parametric and semiparametric mixture models have been widely used in applications from many areas, and it is often of interest to test homogeneity in these models. However, hypothesis testing is nonstandard due to the fact that several regularity conditions do not hold under the null hypothesis. We consider a semiparametric mixture case-control model, in the sense that the density ratio of two distributions is assumed to be of an exponential form, while the baseline density is unspecified. This model was first considered by Qin and Liang (2011, biometrics), and they proposed a modified score statistic for testing homogeneity. In this paper, we consider alternative testing procedures based on supremum statistics, which could improve power against certain types of alternatives. We demonstrate the connection and comparison among proposed and existing these approaches. In addition, we provide a unified theoretical justification of the supremum test and other existing test statistics from an empirical likelihood perspective. The finite sample performance of the supremum test statistics were evaluated in simulation studies.

Modeling Linkage Disequilibrium Increases Accuracy of Polygenic Risk Scores

Polygenic risk scores have shown great promise in predicting complex disease risk and will become more accurate as training sample sizes increase. The standard approach for calculating risk scores involves linkage disequilibrium (LD)-based marker pruning and applying a p value threshold to association statistics, but this discards information and can reduce predictive accuracy. We introduce LDpred, a method that infers the posterior mean effect size of each marker by using a prior on effect sizes and LD information from an external reference panel. Theory and simulations show that LDpred outperforms the approach of pruning followed by thresholding, particularly at large sample sizes. Accordingly, predicted R(2) increased from 20.1% to 25.3% in a large schizophrenia dataset and from 9.8% to 12.0% in a large multiple sclerosis dataset. A similar relative improvement in accuracy was observed for three additional large disease datasets and for non-European schizophrenia samples. The advantage of LDpred over existing methods will grow as sample sizes increase.

Genome-wide association study in obsessive-compulsive disorder: results from the OCGAS

Obsessive-compulsive disorder (OCD) is a psychiatric condition characterized by intrusive thoughts and urges and repetitive, intentional behaviors that cause significant distress and impair functioning. The OCD Collaborative Genetics Association Study (OCGAS) is comprised of comprehensively assessed OCD patients with an early age of OCD onset. After application of a stringent quality control protocol, a total of 1065 families (containing 1406 patients with OCD), combined with population-based samples (resulting in a total sample of 5061 individuals), were studied. An integrative analyses pipeline was utilized, involving association testing at single-nucleotide polymorphism (SNP) and gene levels (via a hybrid approach that allowed for combined analyses of the family- and population-based data). The smallest P-value was observed for a marker on chromosome 9 (near PTPRD, P=4.13 × 10(-)(7)). Pre-synaptic PTPRD promotes the differentiation of glutamatergic synapses and interacts with SLITRK3. Together, both proteins selectively regulate the development of inhibitory GABAergic synapses. Although no SNPs were identified as associated with OCD at genome-wide significance level, follow-up analyses of genome-wide association study (GWAS) signals from a previously published OCD study identified significant enrichment (P=0.0176). Secondary analyses of high-confidence interaction partners of DLGAP1 and GRIK2 (both showing evidence for association in our follow-up and the original GWAS study) revealed a trend of association (P=0.075) for a set of genes such as NEUROD6, SV2A, GRIA4, SLC1A2 and PTPRD. Analyses at the gene level revealed association of IQCK and C16orf88 (both P<1 × 10(-)(6), experiment-wide significant), as well as OFCC1 (P=6.29 × 10(-)(5)). The suggestive findings in this study await replication in larger samples.

Simultaneous estimation of the locations and effects of multiple disease loci in case-control studies

The genetic basis of complex diseases often involves multiple causative loci. Under such a disease etiology, assuming one disease locus in linkage disequilibrium mapping is likely to induce bias and lead to efficiency loss in disease locus estimation. An approach is needed for simultaneously localizing multiple functional loci within the same region. However, due to the increasing number of parameters accompanying disease loci, these estimates can be computationally infeasible. To circumvent this problem, we propose to estimate the main and two-adjacent-locus joint effects and a nuisance parameter at the disease loci separately through a linear approximation. Estimates of the genetic effects are entered into a generalized estimating equation to estimate disease loci, and the procedure is conducted iteratively until convergence. The proposed method provides estimates and confidence intervals (CIs) for the disease loci, the genetic main effects, and the joint effects of two adjacent disease loci, with the CIs for the disease loci providing useful regions for further fine-mapping. We apply the proposed approach to a data example of case-control studies. Results of the simulations and data example suggest that the developed method performs well in terms of bias, variance, and coverage probability under scenarios with up to three disease loci.

Partitioning heritability of regulatory and cell-type-specific variants across 11 common diseases

Regulatory and coding variants are known to be enriched with associations identified by genome-wide association studies (GWASs) of complex disease, but their contributions to trait heritability are currently unknown. We applied variance-component methods to imputed genotype data for 11 common diseases to partition the heritability explained by genotyped SNPs (hg(2)) across functional categories (while accounting for shared variance due to linkage disequilibrium). Extensive simulations showed that in contrast to current estimates from GWAS summary statistics, the variance-component approach partitions heritability accurately under a wide range of complex-disease architectures. Across the 11 diseases DNaseI hypersensitivity sites (DHSs) from 217 cell types spanned 16% of imputed SNPs (and 24% of genotyped SNPs) but explained an average of 79% (SE = 8%) of hg(2) from imputed SNPs (5.1× enrichment; p = 3.7 × 10(-17)) and 38% (SE = 4%) of hg(2) from genotyped SNPs (1.6× enrichment, p = 1.0 × 10(-4)). Further enrichment was observed at enhancer DHSs and cell-type-specific DHSs. In contrast, coding variants, which span 1% of the genome, explained <10% of hg(2) despite having the highest enrichment. We replicated these findings but found no significant contribution from rare coding variants in independent schizophrenia cohorts genotyped on GWAS and exome chips. Our results highlight the value of analyzing components of heritability to unravel the functional architecture of common disease.

Genome-wide association analysis identifies 13 new risk loci for schizophrenia

Schizophrenia is a heritable disorder with substantial public health impact. We conducted a multi-stage genome-wide association study (GWAS) for schizophrenia beginning with a Swedish national sample (5,001 cases, 6,243 controls) followed by meta-analysis with prior schizophrenia GWAS (8,832 cases, 12,067 controls) and finally by replication of SNPs in 168 genomic regions in independent samples (7,413 cases, 19,762 controls, and 581 trios). In total, 22 regions met genome-wide significance (14 novel and one previously implicated in bipolar disorder). The results strongly implicate calcium signaling in the etiology of schizophrenia, and include genome-wide significant results for CACNA1C and CACNB2 whose protein products interact. We estimate that ∼8,300 independent and predominantly common SNPs contribute to risk for schizophrenia and that these collectively account for most of its heritability. Common genetic variation plays an important role in the etiology of schizophrenia, and larger studies will allow more detailed understanding of this devastating disorder.

Body mass index and risk of incident hypertension over the life course: the Johns Hopkins Precursors Study

BACKGROUND

The obesity-hypertension link over the life course has not been well characterized, although the prevalence of obesity and hypertension is increasing in the United States.

METHODS AND RESULTS

We studied the association of body mass index (BMI) in young adulthood, into middle age, and through late life with risk of developing hypertension in 1132 white men of The Johns Hopkins Precursors Study, a prospective cohort study. Over a median follow-up period of 46 years, 508 men developed hypertension. Obesity (BMI ≥30 kg/m(2)) in young adulthood was strongly associated with incident hypertension (hazard ratio, 4.17; 95% confidence interval, 2.34-7.42). Overweight (BMI 25 to <30 kg/m(2)) also signaled increased risk (hazard ratio, 1.58; 95% confidence interval, 1.28-1.96). Men of normal weight at age 25 years who became overweight or obese at age 45 years were at increased risk compared with men of normal weight at both times (hazard ratio, 1.57; 95% confidence interval, 1.20-2.07), but not men who were overweight or obese at age 25 years who returned to normal weight at age 45 years (hazard ratio, 0.91; 95% confidence interval, 0.43-1.92). After adjustment for time-dependent number of cigarettes smoked, cups of coffee taken, alcohol intake, physical activity, parental premature hypertension, and baseline BMI, the rate of change in BMI over the life course increased the risk of incident hypertension in a dose-response fashion, with the highest risk among men with the greatest increase in BMI (hazard ratio, 2.52; 95% confidence interval, 1.82-3.49).

CONCLUSIONS

Our findings underscore the importance of higher weight and weight gain in increasing the risk of hypertension from young adulthood through middle age and into late life.

Homeobox genes in obsessive-compulsive disorder

BACKGROUND

Despite evidence that obsessive-compulsive disorder (OCD) is a familial neuropsychiatric condition, progress aimed at identifying genetic determinants of the disorder has been slow. The OCD Collaborative Genetics Study (OCGS) has identified several OCD susceptibility loci through linkage analysis.

METHODS

In this study we investigate two regions on chromosomes 15q and 1q by first refining the linkage region using additional short tandem repeat polymorphic (STRP) markers. We then performed association analysis on single nucleotide polymorphisms (SNP) genotyped (markers placed every 2-4 kb) in the linkage regions in the OCGS sample of 376 rigorously phenotyped affected families.

RESULTS

Three SNPs are most strongly associated with OCD: rs11854486 (P = 0.00005 [0.046 after adjustment for multiple tests]; genetic relative risk (GRR) = 11.1 homozygous and 1.6 heterozygous) and rs4625687 [P = 0.00007 (after adjustment = 0.06); GRR = 2.4] on 15q; and rs4387163 (P = 0.0002 (after adjustment = 0.08); GRR = 1.97) on 1q. The first SNP is adjacent to NANOGP8, the second SNP is in MEIS2, and the third is 150 kb between PBX1 and LMX1A.

CONCLUSIONS

All the genes implicated by association signals are homeobox genes and are intimately involved in neurodevelopment. PBX1 and MEIS2 exert their effects by the formation of a heterodimeric complex, which is involved in development of the striatum, a brain region involved in the pathophysiology of OCD. NANOGP8 is a retrogene of NANOG, a homeobox transcription factor known to be involved in regulation of neuronal development. These findings need replication; but support the hypothesis that genes involved in striatal development are implicated in the pathogenesis of OCD.

Cognitive functioning in compulsive hoarding

The aim of this study is to determine whether neurocognitive performance distinguishes individuals with compulsive hoarding (CH) from those with obsessive-compulsive disorder (OCD). Compared to control subjects, OCD patients and CHs scored significantly worse on the Serial Reaction Time Task suggesting disturbed implicit memory in both patient groups. On the Iowa Gambling Task, an overall learning progression difference over time was found between the CHs, OCD group, and control subjects, suggesting differences in decision-making between the groups. The groups did not differ in performance on the Stop Signal Reaction Time Task (motor inhibition). This study found evidence for impaired implicit memory in CHs, but also in OCD patients, albeit less severe. There was evidence that OCD patients learned more slowly on a decision-making task than CHs and control subjects. This latter finding provides some evidence to suggest that CH and OCD have, at least on this one measure, differing cognitive substrates.

Likelihood ratio testing for admixture models with application to genetic linkage analysis

We consider likelihood ratio tests (LRT) and their modifications for homogeneity in admixture models. The admixture model is a two-component mixture model, where one component is indexed by an unknown parameter while the parameter value for the other component is known. This model is widely used in genetic linkage analysis under heterogeneity in which the kernel distribution is binomial. For such models, it is long recognized that testing for homogeneity is nonstandard, and the LRT statistic does not converge to a conventional χ(2)  distribution. In this article, we investigate the asymptotic behavior of the LRT for general admixture models and show that its limiting distribution is equivalent to the supremum of a squared Gaussian process. We also discuss the connection and comparison between LRT and alternative approaches such as modifications of LRT and score tests, including the modified LRT (Fu, Chen, and Kalbfleisch, 2006, Statistica Sinica 16, 805-823). The LRT is an omnibus test that is powerful to detect general alternative hypotheses. In contrast, alternative approaches may be slightly more powerful to detect certain type of alternatives, but much less powerful for others. Our results are illustrated by simulation studies and an application to a genetic linkage study of schizophrenia.

Association between a serotonin transporter promoter polymorphism (5HTTLPR) and personality disorder traits in a community sample

BACKGROUND

The serotonin transporter (SERT) polymorphism (5HTTLPR) has been reported to be associated with several psychiatric conditions. Specific personality disorders could be intermediate factors in the known relationship between 5HTTLPR and psychiatric disorders. This is the first study to test the association between this polymorphism and dimensions of all DSM-IV personality disorders in a community sample.

METHODS

374 white participants were assessed by clinical psychologists using the International Personality Disorder Examination (IPDE). Associations between dimensions of each DSM-IV personality disorder and the long (l) and short (s) alleles of the 5HTTLPR were evaluated using non-parametric tests and regression models.

RESULTS

The s allele of the 5HTTLPR polymorphism was significantly associated with higher avoidant personality trait scores in the whole sample. Males with the s allele had a significantly lower likelihood of higher obsessive-compulsive personality disorder (OCPD) trait scores, whereas females with the s allele were likely to have higher OCPD personality trait scores.

CONCLUSION

This paper provides preliminary data on the relationship between personality disorders and the 5HTTLPR polymorphism. The relationship of the s allele and avoidant PD is consistent with findings of a nonspecific relationship of this polymorphism to anxiety and depressive disorders. Concerning the unusual sexual dimorphic result with OCPD, several hypotheses are presented. These findings need further replication, including a more detailed study of additional variants in SERT.

Hypothesis testing in a mixture case-control model

We consider a problem of testing mixture proportions using two-sample data, one from group one and the other from a mixture of groups one and two with unknown proportion, λ, for being in group two. Various statistical applications, including microarray study, infectious epidemiological studies, case-control studies with contaminated controls, clinical trials allowing "nonresponders," genetic studies for gene mutation, and fishery applications can be formulated in this setup. Under the assumption that the log ratio of probability (density) functions from the two groups is linear in the observations, we propose a generalized score test statistic to test the mixture proportion. Under some regularity conditions, it is shown that this statistic converges to a weighted chi-squared random variable under the null hypothesis of λ= 0, where the weight depends only on the sampling fraction of both groups. The permutation method is used to provide more reliable finite sample approximation. Simulation results and two real data applications are presented.

Comprehensive family-based association study of the glutamate transporter gene SLC1A1 in obsessive-compulsive disorder

SLC1A1 encodes a neuronal glutamate transporter and is a promising candidate gene for obsessive-compulsive disorder (OCD). Several independent research groups have reported significant associations between OCD and single nucleotide polymorphisms (SNPs) in this gene. Previously, we evaluated 13 SNPs in, or near, SLC1A1 and reported a strong association signal with rs301443, a SNP 7.5 kb downstream of the gene [Shugart et al. (2009); Am J Med Genet Part B 150B:886–892]. The aims of the current study were first, to further investigate this finding by saturating the region around rs301443; and second, to explore the entire gene more thoroughly with a dense panel of SNP markers. We genotyped an additional 111 SNPs in or near SLC1A1, covering from 9 kb upstream to 84 kb downstream of the gene at average spacing of 1.7 kb per SNP, and conducted family-based association analyses in 1,576 participants in 377 families.We found that none of the surrounding markers were in linkage disequilibrium with rs301443, nor were any associated with OCD. We also found that SNP rs4740788, located about 8.8 kb upstream of the gene, was associated with OCD in all families (P = 0.003) and in families with male affecteds (P = 0.002). A three-SNP haplotype (rs4740788–rs10491734–rs10491733) was associated with OCD in the total sample (P = 0.00015) and in families with male affecteds (P = 0.0007). Although of nominal statistical significance considering the number of comparisons, these findings provide further support for the involvement of SLC1A1 in the pathogenesis of OCD.

Linkage and association on 8p21.2-p21.1 in schizophrenia

In the past decade, we and others have consistently reported linkage to a schizophrenia (SZ) susceptibility region on chromosome 8p21. Most recently, in the largest SZ linkage sample to date, a multi-site international collaboration performed a SNP-based linkage scan (~6,000 SNPs; 831 pedigrees; 121 from Johns Hopkins (JHU)), that showed the strongest evidence for linkage in a 1 Mb region of chr 8p21 from rs1561817 to rs9797 (Z(max) = 3.22, P = 0.0004) [Holmans et al. 2009. Mol Psychiatry]. We have investigated this 8p21 peak region further in two ways: first by linkage and family-based association in 106 8p-linked European-Caucasian (EUC) JHU pedigrees using 1,402 SNPs across a 4.4 Mb region surrounding the peak; second, by an independent case-control association study in the genetically more homogeneous Ashkenazim (AJ) (709 cases, 1,547 controls) using 970 SNPs in a further narrowed 2.8 Mb region. Family-based association analyses in EUC pedigrees and case-control analyses in AJ samples reveal significant associations for SNPs in and around DPYSL2 and ADRA1A, candidate genes previously associated with SZ in our work and others. Further, several independent gene expression studies have shown that DPYSL2 is differentially expressed in SZ brains [Beasley et al. 2006. Proteomics 6(11):3414–3425; Edgar et al. 2000. Mol Psychiatry 5(1):85–90; Johnston-Wilson et al. 2000. Mol Psychiatry 5(2):142–149] or in response to psychosis-inducing pharmaceuticals [Iwazaki et al. 2007. Proteomics 7(7):1131–1139; Paulson et al. 2004. Proteomics 4(3):819–825]. Taken together, this work further supports DPYSL2 and the surrounding genomic region as a susceptibility locus for SZ.

Separation anxiety disorder in OCD

BACKGROUND

A history of separation anxiety disorder (SAD) is frequently reported by patients with obsessive-compulsive disorder (OCD). The purpose of this study was to determine if there are clinical differences between OCD-affected individuals with, versus without, a history of SAD.

METHODS

Using data collected during the OCD Collaborative Genetic Study, we studied 470 adult OCD participants; 80 had a history of SAD, whereas 390 did not. These two groups were compared as to onset and severity of OCD, lifetime prevalence of Axis I disorders, and number of personality disorder traits.

RESULTS

OCD participants with a history of SAD were significantly younger than the non-SAD group (mean, 34.2 versus 42.2 years; P<.001). They had an earlier age of onset of OCD symptoms (mean, 8.0 versus 10.5 years; P<.003) and more severe OCD, as measured by the Yale-Brown Obsessive Compulsive Scale (mean, 27.5 versus 25.0; P<.005). In addition, those with a history of SAD had a significantly greater lifetime prevalence of agoraphobia (odds ratio (OR) = 2.52, 95% confidence interval (CI) = 1.4-4.6, P<.003), panic disorder (OR = 1.84, CI = 1.03-3.3 P<.04), social phobia (OR = 1.69, CI 1.01-2.8, P<.048), after adjusting for age at interview, age at onset of OCD, and OCD severity in logistic regression models. There was a strong relationship between the number of dependent personality disorder traits and SAD (adjusted OR = 1.42, CI = 1.2-1.6, P<.001).

CONCLUSIONS

A history of SAD is associated with anxiety disorders and dependent personality disorder traits in individuals with OCD.

Familial aggregation of panic disturbances in Parkinson's disease

Panic disorder has an elevated prevalence in Parkinson's disease (PD). To explore the basis for this co-occurrence, the familial aggregation of panic disorder was examined in patients with PD. Probands and relatives of patients with PD and panic disorder (PD-PANIC; N=20, N=115) and control probands with PD and no active psychiatric illness (PD-NA; N=17, N=108) were interviewed by phone, using a structured interview to determine panic status. Lifetime prevalence of panic and "panic-like" disorders was higher in PD-PANIC than in PD-NA relatives. Panic and "panic-like" disorders are familial disorders in PD.

Detection of SNP-SNP interactions in trios of parents with schizophrenic children

Schizophrenia (SZ) is a heritable and complex psychiatric disorder with an estimated worldwide prevalence of about 1%. Research on the risk factors for SZ has thus far yielded few clues to causes, but has pointed to a heterogeneous etiology that likely involves multiple genes and gene-environment interactions. In this manuscript, we apply a novel method (trio logic regression, Li et al., 2009) to case-parent trio data from a SZ candidate gene study conducted on families of Ashkenazi Jewish descent, and demonstrate the method's ability to detect multi-gene models for SZ risk in the family-based design. In particular, we demonstrate how this method revealed a genotype-phenotype association that includes an allele without marginal effect.

Alcohol consumption and domain-specific cognitive function in older adults: longitudinal data from the Johns Hopkins Precursors Study

OBJECTIVES

The association of alcohol consumption with performance in different cognitive domains has not been well studied.

METHODS

The Johns Hopkins Precursors Study was used to examine associations between prospectively collected information about alcohol consumption ascertained on multiple occasions starting at age 55 years on average with domain-specific cognition at age 72 years. Cognitive variables measured phonemic and semantic fluency, attention, verbal memory, and global cognition.

RESULTS

Controlling for age, hypertension, smoking status, sex, and other cognitive variables, higher average weekly quantity and frequency of alcohol consumed in midlife were associated with lower phonemic fluency. There were no associations with four other measures of cognitive function. With respect to frequency of alcohol intake, phonemic fluency was significantly better among those who drank three to four alcoholic beverages per week as compared with daily or almost daily drinkers. A measure of global cognition was not associated with alcohol intake at any point over the follow-up.

DISCUSSION

Results suggest that higher alcohol consumption in midlife may impair some components of executive function in late life.

Reliability of alcohol recall after 15 years and 23 years of follow-up in the Johns Hopkins Precursors Study

OBJECTIVE

Recall of past alcohol intake is used in many studies of chronic disease, but few studies have been able to examine its long-term reliability.

METHOD

We sought to assess the reliability of recalled alcohol intake assessed at an average age of 70 years in 2001, after 15 and 23 years of follow-up, in a prospective study of medical students in classes 1948 to 1964.

RESULTS

Average reported alcohol intake 15 years and 23 years prior were 6.3 and 7.4 drinks per week, respectively. Recall of alcohol intake overestimated the concurrently reported intake after 15 years by a mean of 0.47 (95% CI [0.10, 0.85]) drinks per week and underestimated intake after 23 years by a mean of 0.79 (95% CI [-1.27, -0.30]) drinks per week, mostly driven by differences between concurrently reported and recalled distilled spirits consumption. Characteristics associated with underestimation of alcohol recall were age of 71 years or older in 2001, self-report of memory difficulties, and self-report of difficulties in physical functioning. In multivariate regression analyses combining 15- and 23-year recall, subjects who reported consumption of more than 14 alcoholic drinks per week in 2001 marginally overestimated recall by slightly more than 1 drink per week (M = 1.18 drinks/week, 95% CI [-0.03, 2.40]).

CONCLUSIONS

Although significant differences were detected, recalled alcohol intake after 15 and 23 years of follow-up is remarkably reliable.

Incorporation of covariates in simultaneous localization of two linked loci using affected relative pairs

Background

Many dichotomous traits for complex diseases are often involved more than one locus and/or associated with quantitative biomarkers or environmental factors. Incorporating these quantitative variables into linkage analysis as well as localizing two linked disease loci simultaneously could therefore improve the efficiency in mapping genes. We extended the robust multipoint Identity-by-Descent (IBD) approach with incorporation of covariates developed previously to simultaneously estimate two linked loci using different types of affected relative pairs (ARPs).

Results

We showed that the efficiency was enhanced by incorporating a quantitative covariate parametrically or non-parametrically while localizing two disease loci using ARPs. In addition to its help in identifying factors associated with the disease and in improving the efficiency in estimating disease loci, this extension also allows investigators to account for heterogeneity in risk-ratios for different ARPs. Data released from the collaborative study on the genetics of alcoholism (COGA) for Genetic Analysis Workshop 14 (GAW 14) were used to illustrate the application of this extended method.

Conclusions

The simulation studies and example illustrated that the efficiency in estimating disease loci was demonstratively enhanced by incorporating a quantitative covariate and by using all relative pairs while mapping two linked loci simultaneously.

On the asymptotic behaviour of the pseudolikelihood ratio test statistic with boundary problems

This paper considers the asymptotic distribution of the likelihood ratio statistic T for testing a subset of parameter of interest θ, θ = (γ, η), H(0) : γ = γ(0), based on the pseudolikelihood L(θ, ϕ̂), where ϕ̂ is a consistent estimator of ϕ, the nuisance parameter. We show that the asymptotic distribution of T under H(0) is a weighted sum of independent chi-squared variables. Some sufficient conditions are provided for the limiting distribution to be a chi-squared variable. When the true value of the parameter of interest, θ(0), or the true value of the nuisance parameter, ϕ(0), lies on the boundary of parameter space, the problem is shown to be asymptotically equivalent to the problem of testing the restricted mean of a multivariate normal distribution based on one observation from a multivariate normal distribution with misspecified covariance matrix, or from a mixture of multivariate normal distributions. A variety of examples are provided for which the limiting distributions of T may be mixtures of chi-squared variables. We conducted simulation studies to examine the performance of the likelihood ratio test statistics in variance component models and teratological experiments.

Incorporating covariates into multipoint association mapping in the case-parent design

BACKGROUND/AIMS

To improve the efficiency of disease locus localization in association mapping using case-parent designs and to assess or account for the main covariate effects and gene-covariate interaction effects, while localizing the disease locus.

METHODS

The present study extends a multipoint fine-mapping approach to incorporate covariates into the association mapping of case-parent designs through parametric and non-parametric modeling. This approach is based on the expected preferential-allele-transmission statistics for transmission from either parent to an affected child.

RESULTS

Simulation studies indicate that the efficiency in estimating the disease locus increases considerably when incorporating a covariate associated with the disease. This is especially true when the genetic effect of the disease locus is small. The proposed approach was applied to a young-onset hypertension data sample. The relative efficiency of estimating the locus of young-onset hypertension increases 110-fold after incorporating triglyceride into the association mapping while localizing the disease variant in the lipoprotein lipase gene in the non-parametric model. By incorporating the information of SNP variants into the fine-mapping, the proposed method further assesses the gene-gene interactions between the SNP and the disease locus.

CONCLUSION

With the incorporation of covariates, the proposed method cannot only improve efficiency in estimating disease loci, but can also elucidate the etiology of a complex disease.

A screen of SLC1A1 for OCD-related alleles

SLC1A1, which encodes the neuronal and epithelial glutamate transporter, is a promising candidate gene for obsessive-compulsive disorder (OCD). In this study, we conducted capillary electrophoresis single-strand conformation polymorphism (CE-SSCP) screen for all 12 identified exons, including all coding regions and approximately 50 bp of flanking introns of the human SLC1A1 in 378 OCD-affected individuals. Full sequencing was completed on samples that showed an aberrant SSCP tracing for identification of the underlying sequence variants. Our aim was to determine if there are differences in the frequencies of relatively common alleles, or rare functional alleles, in 378 OCD cases and 281 ethnically matched controls. We identified one nonsynonymous coding SNP (c.490A > G, T164A) and three synonymous coding SNP (c.81G > C, A27A; c.414A > G, T138T; c.1110T > C, T370T) in case samples. We found no statistical differences in genotype and allele frequencies of common cSNPs in SLC1A1 between the OCD cases and controls. The rare variant T164A was found only in one family. Further investigation of this variant is necessary to determine whether and how it is related to OCD. There was no other evidence of significant accumulation of deleterious coding mutations in SLC1A1 in the OCD cases.

A family-based association study of the glutamate transporter gene SLC1A1 in obsessive-compulsive disorder in 378 families

SLC1A encodes the neuronal and epithelial glutamate transporter and was previously tested as a candidate for obsessive-compulsive disorder (OCD) by several research groups. Recently, three independent research groups reported significant association findings between OCD and several genetic variants in SLC1A1. This study reports the results from a family-based association study, which examined the association between 13 single nucleotide polymorphisms (SNPs) within or in proximity to the SLC1A1 gene. Although we did not replicate association findings for those significant SNPs reported by previous studies, our study indicated a strong association signal with the SNP RS301443 (P-value = 0.000067; Bonferroni corrected P-value = 0.0167) under a dominant model, with an estimated odds ratio of 3.5 (confidence interval: 2.66-4.50). Further, we conducted single SNP analysis after stratifying the full data set by the gender status of affected in each family. The P-value for RS301443 in families with the male affected was 0.00027, and the P-value in families with female affected was 0.076. The fact that we identified a signal which was not previously reported by the other research groups may be due to differences in study designs and sample ascertainment. However, it is also possible that this significant SNP may be part of a regulator for SLC1A1, given that it is roughly 7.5 kb away from the boundary of the SLC1A1 gene. However, this potential finding needs to be validated biologically. Further functional studies in this region are planned by this research group.

Obsessive-compulsive disorder: subclassification based on co-morbidity

BACKGROUND

Obsessive-compulsive disorder (OCD) is probably an etiologically heterogeneous condition. Many patients manifest other psychiatric syndromes. This study investigated the relationship between OCD and co-morbid conditions to identify subtypes.

METHOD

Seven hundred and six individuals with OCD were assessed in the OCD Collaborative Genetics Study (OCGS). Multi-level latent class analysis was conducted based on the presence of eight co-morbid psychiatric conditions [generalized anxiety disorder (GAD), major depression, panic disorder (PD), separation anxiety disorder (SAD), tics, mania, somatization disorders (Som) and grooming disorders (GrD)]. The relationship of the derived classes to specific clinical characteristics was investigated.

RESULTS

Two and three classes of OCD syndromes emerge from the analyses. The two-class solution describes lesser and greater co-morbidity classes and the more descriptive three-class solution is characterized by: (1) an OCD simplex class, in which major depressive disorder (MDD) is the most frequent additional disorder; (2) an OCD co-morbid tic-related class, in which tics are prominent and affective syndromes are considerably rarer; and (3) an OCD co-morbid affective-related class in which PD and affective syndromes are highly represented. The OCD co-morbid tic-related class is predominantly male and characterized by high conscientiousness. The OCD co-morbid affective-related class is predominantly female, has a young age at onset, obsessive-compulsive personality disorder (OCPD) features, high scores on the 'taboo' factor of OCD symptoms, and low conscientiousness.

CONCLUSIONS

OCD can be classified into three classes based on co-morbidity. Membership within a class is differentially associated with other clinical characteristics. These classes, if replicated, should have important implications for research and clinical endeavors.

Genomewide linkage scan of schizophrenia in a large multicenter pedigree sample using single nucleotide polymorphisms

A genomewide linkage scan was carried out in eight clinical samples of informative schizophrenia families. After all quality control checks, the analysis of 707 European-ancestry families included 1615 affected and 1602 unaffected genotyped individuals, and the analysis of all 807 families included 1900 affected and 1839 unaffected individuals. Multipoint linkage analysis with correction for marker-marker linkage disequilibrium was carried out with 5861 single nucleotide polymorphisms (SNPs; Illumina version 4.0 linkage map). Suggestive evidence for linkage (European families) was observed on chromosomes 8p21, 8q24.1, 9q34 and 12q24.1 in nonparametric and/or parametric analyses. In a logistic regression allele-sharing analysis of linkage allowing for intersite heterogeneity, genomewide significant evidence for linkage was observed on chromosome 10p12. Significant heterogeneity was also observed on chromosome 22q11.1. Evidence for linkage across family sets and analyses was most consistent on chromosome 8p21, with a one-LOD support interval that does not include the candidate gene NRG1, suggesting that one or more other susceptibility loci might exist in the region. In this era of genomewide association and deep resequencing studies, consensus linkage regions deserve continued attention, given that linkage signals can be produced by many types of genomic variation, including any combination of multiple common or rare SNPs or copy number variants in a region.

Sapap3 and pathological grooming in humans: Results from the OCD collaborative genetics study

SAP90/PSD95-associated protein (SAPAP) family proteins are post-synaptic density (PSD) components that interact with other proteins to form a key scaffolding complex at excitatory (glutamatergic) synapses. A recent study found that mice with a deletion of the Sapap3 gene groomed themselves excessively, exhibited increased anxiety-like behaviors, and had cortico-striatal synaptic defects, all of which were preventable with lentiviral-mediated expression of Sapap3 in the striatum; the behavioral abnormalities were also reversible with fluoxetine. In the current study, we sought to determine whether variation within the human Sapap3 gene was associated with grooming disorders (GDs: pathologic nail biting, pathologic skin picking, and/or trichotillomania) and/or obsessive-compulsive disorder (OCD) in 383 families thoroughly phenotyped for OCD genetic studies. We conducted family-based association analyses using the FBAT and GenAssoc statistical packages. Thirty-two percent of the 1,618 participants met criteria for a GD, and 65% met criteria for OCD. Four of six SNPs were nominally associated (P < 0.05) with at least one GD (genotypic relative risks: 1.6-3.3), and all three haplotypes were nominally associated with at least one GD (permuted P < 0.05). None of the SNPs or haplotypes were significantly associated with OCD itself. We conclude that Sapap3 is a promising functional candidate gene for human GDs, though further work is necessary to confirm this preliminary evidence of association.

Interaction between the serotonin transporter gene (5-HTTLPR), stressful life events, and risk of depression: a meta-analysis

CONTEXT

Substantial resources are being devoted to identify candidate genes for complex mental and behavioral disorders through inclusion of environmental exposures following the report of an interaction between the serotonin transporter linked polymorphic region (5-HTTLPR) and stressful life events on an increased risk of major depression.

OBJECTIVE

To conduct a meta-analysis of the interaction between the serotonin transporter gene and stressful life events on depression using both published data and individual-level original data.

DATA SOURCES

Search of PubMed, EMBASE, and PsycINFO databases through March 2009 yielded 26 studies of which 14 met criteria for the meta-analysis.

STUDY SELECTION

Criteria for studies for the meta-analyses included published data on the association between 5-HTTLPR genotype (SS, SL, or LL), number of stressful life events (0, 1, 2, > or = 3) or equivalent, and a categorical measure of depression defined by the Diagnostic and Statistical Manual of Mental Disorders (Fourth Edition) or the International Statistical Classification of Diseases, 10th Revision (ICD-10) or use of a cut point to define depression from standardized rating scales. To maximize our ability to use a common framework for variable definition, we also requested original data from all studies published prior to 2008 that met inclusion criteria. Of the 14 studies included in the meta-analysis, 10 were also included in a second sex-specific meta-analysis of original individual-level data.

DATA EXTRACTION

Logistic regression was used to estimate the effects of the number of short alleles at 5-HTTLPR, the number of stressful life events, and their interaction on depression. Odds ratios (ORs) and 95% confidence intervals (CIs) were calculated separately for each study and then weighted averages of the individual estimates were obtained using random-effects meta-analysis. Both sex-combined and sex-specific meta-analyses were conducted. Of a total of 14,250 participants, 1769 were classified as having depression; 12,481 as not having depression.

RESULTS

In the meta-analysis of published data, the number of stressful life events was significantly associated with depression (OR, 1.41; 95% CI,1.25-1.57). No association was found between 5-HTTLPR genotype and depression in any of the individual studies nor in the weighted average (OR, 1.05; 95% CI, 0.98-1.13) and no interaction effect between genotype and stressful life events on depression was observed (OR, 1.01; 95% CI, 0.94-1.10). Comparable results were found in the sex-specific meta-analysis of individual-level data.

CONCLUSION

This meta-analysis yielded no evidence that the serotonin transporter genotype alone or in interaction with stressful life events is associated with an elevated risk of depression in men alone, women alone, or in both sexes combined.

Antidepressant-associated mania: soon after switch from fluoxetine to mirtazapine in an elderly woman with mixed depressive features

Mirtazapine augmentation to a serotonin-reuptake inhibitor has been proposed to boost antidepressant effects and more likely to induce manic switch. Such a combined antidepressant therapy strategy should be used carefully if the patient's refractoriness is attributable to mixed depressive features. Mixed depression is more difficult to be treated by antidepressant monotherapy and related to higher risk of manic switch during treatment. We report a case with no previous history of bipolar disorder, whereas developed full-blown psychotic manic symptoms soon after switch from fluoxetine to mirtazapine. The patient's premorbid characters and clinical presentations suggested an implicit bipolarity that predisposed her to a manic switch. Her manic switch was likely to be triggered by a simulated combined effect because of complex drug interactions during shifting from fluoxetine to mirtazapine. For patients in mixed depressive states, mood stabilizers are preferable to antidepressants.

Maternal transmission effects of the PAX genes among cleft case-parent trios from four populations

Isolated cleft lip with or without cleft palate (CL/P) is among the most common human birth defects, with a prevalence of 1 in 700 live births. The paired box (PAX) genes have been suggested as candidate genes for CL/P based largely on mouse models; however, few human studies have focused on this gene family. This study tests for association between markers in four PAX genes and CL/P using a case-parent trio design considering parent-of-origin effects. Trios from four populations (76 from Maryland, 146 from Taiwan, 35 from Singapore, and 40 from Korea) were genotyped for 34 single nucleotide polymorphisms (SNPs) in the PAX3, PAX6, PAX7, and PAX9 genes. We performed the transmission disequilibrium test (TDT) on individual SNPs. Parent-of-origin effects were assessed using the transmission asymmetry test (TAT) and the parent-of-origin likelihood ratio test (PO-LRT). TDT analysis showed one SNP (rs766325) in PAX7 yielding evidence of linkage and association when parent-of-origin was not considered, with an OR(transmission)=1.62 (P=0.003), and five SNPs in PAX6 (including two pairs in near perfect linkage disequilibrium). TAT analysis of all trios revealed two SNPs in PAX7 and four SNPs in PAX3 showing significant excess maternal transmission. For these six SNPs, the maternal OR(transmission) ranged between 1.74 and 2.40, and PO-LRT was also significant (P-values=0.035-0.012). When this analysis was limited to trios with male cases, SNPs in PAX7 showed higher maternal OR(transmission) and greater significance. PAX genes may influence the risk of CL/P through maternal effects, possibly imprinting, which seems to be stronger among male cases.

Gender differences in genetic linkage and association on 11p15 in obsessive-compulsive disorder families

Several clinical and genetic studies have reported gender differences in obsessive-compulsive disorder (OCD). Previously, we conducted a linkage genome scan using multipoint allele-sharing methods to test for linkage in 219 families participating in the OCD Collaborative Genetics Study. When these families were stratified by proband's gender, suggestive linkage to chromosome 11p15 at marker D11S2362 (KAC(all) = 2.92, P = 0.00012) was detected in families with male probands, but not in the ones with female probands. We have since conducted fine mapping with a denser microsatellite marker panel in the region of 11p15, and detected a significant linkage signal at D11S4146 (KAC(all) = 5.08, P < 0.00001) in the families of male probands. Subsequently, 632 SNPs were genotyped spanning a 4.0 Mb region of the 1 LOD unit interval surrounding the linkage peak in the original families and an additional 165 families. Six SNPs were associated with OCD (P < 0.001): two SNPs were identified when all the families were included, and four SNPs only in male proband families. No SNP showed significant association with the OCD phenotype only in the families with a female proband. The results suggest a possible gender effect in the etiology of OCD.

Fine mapping on chromosome 10q22-q23 implicates Neuregulin 3 in schizophrenia

Linkage studies have implicated 10q22-q23 as a schizophrenia (SZ) susceptibility locus in Ashkenazi Jewish (AJ) and Han Chinese from Taiwan populations. To further explore our previous linkage signal in the AJ population (NPL score: 4.27, empirical p = 2 x 10(-5)), we performed a peakwide association fine mapping study by using 1414 SNPs across approximately 12.5 Mb in 10q22-q23. We genotyped 1515 AJ individuals, including 285 parent-child trios, 173 unrelated cases, and 487 unrelated controls. We analyzed the binary diagnostic phenotype of SZ and 9 heritable quantitative traits derived from a principal components factor analysis of 73 items from our consensus diagnostic ratings and direct assessment interviews. Although no marker withstood multiple test correction for association with the binary SZ phenotype, we found strong evidence of association by using the "delusion" factor as the quantitative trait at three SNPs (rs10883866, rs10748842, and rs6584400) located in a 13 kb interval in intron 1 of Neuregulin 3 (NRG3). Our best p value from family-based association analysis was 7.26 x 10(-7). We replicated this association in the collection of 173 unrelated AJ cases (p = 1.55 x 10(-2)), with a combined p value of 2.30 x 10(-7). After performing 10,000 permutations of each of the phenotypes, we estimated the empirical study-wide significance across all 9 factors (90,000 permutations) to be p = 2.7 x 10(-3). NRG3 is primarily expressed in the central nervous system and is one of three paralogs of NRG1, a gene strongly implicated in SZ. These biological properties together with our linkage and association results strongly support NRG3 as a gene involved in SZ.

The association between parental bonding and obsessive compulsive disorder in offspring at high familial risk

BACKGROUND

The aim of the current study is to estimate the association between parenting factors derived from the Parental Bonding Instrument (PBI) and a lifetime DSM-IV diagnosis of OCD.

METHOD

Data were from approximately 1200 adults from 465 families assessed as part of a large family and genetic study of OCD. The association of three parenting factors, for fathers and mothers, with offspring OCD status were examined; analyses were stratified by parental OCD status and family loading for OCD (multiplex versus sporadic).

RESULTS

Three factors were derived by principal components factor analysis of the PBI (maternal and paternal care, overprotection and control). Maternal overprotection was associated with OCD in offspring with familial OCD (familial cases) but only if neither parent was affected with OCD, which suggests independent but additive environmental and genetic risk (OR = 5.9, 95% CI 1.2, 29.9, p = 0.031). Paternal care was a protective factor in those not at high genetic risk (sporadic cases) (OR = 0.2, 95% CI 0.0, 0.8, p = 0.027). Maternal overprotection was also associated with offspring OCD in sporadic families (OR = 2.9, 95% CI 1.3, 6.6, p = 0.012). The finding that parental overprotection and care were not associated with offspring OCD when at least one parent had OCD addressed directly the hypothesis of maternal or paternal OCD adversely impacting parenting.

CONCLUSIONS

This study provides evidence that aspects of parenting may contribute to the development of OCD among offspring. Prospective studies of children at risk for OCD are needed to explore the direction of causality.

Differential parental transmission of markers in RUNX2 among cleft case-parent trios from four populations

Isolated cleft lip with or without cleft palate (CL/P) is among the most common human birth defects, with a prevalence around 1 in 700 live births. The Runt-related transcription factor 2 (RUNX2) gene has been suggested as a candidate gene for CL/P based largely on mouse models; however, no human studies have focused on RUNX2 as a risk factor for CL/P. This study examines the association between markers in RUNX2 and isolated, nonsyndromic CL/P using a case-parent trio design, while considering parent-of-origin effects. Case-parent trios from four populations (77 from Maryland, 146 from Taiwan, 35 from Singapore, and 40 from Korea) were genotyped for 24 single nucleotide polymorphisms (SNPs) in the RUNX2 gene. We performed the transmission disequilibrium test on individual SNPs. Parent-of-origin effects were assessed using the transmission asymmetry test and the parent-of-origin likelihood ratio test (PO-LRT). When all trios were combined, the transmission asymmetry test revealed a block of 11 SNPs showing excess maternal transmission significant at the P<0.01 level, plus one SNP (rs1934328) showing excess paternal transmission (P=0.002). For the 11 SNPs showing excess maternal transmission, odds ratios of being transmitted to the case from the mother ranged between 3.00 and 4.00. The parent-of-origin likelihood ratio tests for equality of maternal and paternal transmission were significant for three individual SNPs (rs910586, rs2819861, and rs1934328). Thus, RUNX2 appears to influence risk of CL/P through a parent-of-origin effect with excess maternal transmission.

Excess maternal transmission of markers in TCOF1 among cleft palate case-parent trios from three populations

Isolated cleft palate is among the most common human birth defects. The TCOF1 gene has been suggested as a candidate gene for cleft palate based on animal models. This study tests for association between markers in TCOF1 and isolated, nonsyndromic cleft palate using a case-parent trio design considering parent-of-origin effects. Case-parent trios from three populations (comprising a total of 81 case-parent trios) were genotyped for single nucleotide polymorphisms (SNPs) in the TCOF1 gene. We used the transmission disequilibrium test and the transmission asymmetry test on individual SNPs. When all trios were combined, the odds ratio for transmission of the minor allele, OR(transmission), was significant for SNP rs15251 (OR = 2.88, P = 0.007), as well as rs2255796 and rs2569062 (OR = 2.08, P = 0.03; OR = 2.43, P = 0.041; respectively) when parent of origin was not considered. The transmission asymmetry test also revealed one SNP (rs15251) showing excess maternal transmission significant at the P = 0.005 level (OR = 6.50). Parent-of-origin effects were assessed using the parent-of-origin likelihood ratio test on both SNPs and haplotypes. While the parent-of-origin likelihood ratio test was only marginally significant for this SNP (P = 0.136), analysis of haplotypes of rs2255796 and rs15251 suggested excess maternal transmission. Therefore, these data suggest TCOF1 may influence risk of cleft palate through a parent-of-origin effect.

Evidence for potential relationship between SLC1A1 and a putative genetic linkage region on chromosome 14q to obsessive-compulsive disorder with compulsive hoarding

Obsessive-compulsive disorder (OCD) is likely a disorder involving complex genetic transmission. This suggests that multiple genetic and environmental factors are involved in its etiology. This is complicated further by the probability of genetic heterogeneity for this phenotype. In this report, we describe a preliminary approach to deal with both complexities. SLC1A1, a glutamate transporter gene on chromosome 9p, was originally proposed to be related to OCD based on two linkage studies, and subsequently association of OCD to the gene has been replicated. Additionally, genetic linkage to a subtype of OCD, compulsive hoarding, has been reported on chromosome 14q. We hypothesized that both genomic regions contribute to OCD in some instances. Using the analytic program GENEFINDER we found that conditioning linkage on chromosome 14q to a marker adjacent to SLC1A1, reduced the size of the linkage region on chromosome 14q and provided evidence for interaction between the regions on chromosomes 9p and 14q.

Further development of YBOCS dimensions in the OCD Collaborative Genetics study: symptoms vs. categories

Despite progress in identifying homogeneous subphenotypes of obsessive-compulsive disorder (OCD) through factor analysis of the Yale Brown Obsessive-Compulsive Scale Symptom Checklist (YBOCS-SC), prior solutions have been limited by a reliance on presupposed symptom categories rather than discrete symptoms. Furthermore, there have been few attempts to evaluate the familiality of OCD symptom dimensions. The purpose of this study was to extend prior work by this collaborative group in category-based dimensions by conducting the first-ever exploratory dichotomous factor analysis using individual OCD symptoms, comparing these results to a refined category-level solution, and testing the familiality of derived factors. Participants were 485 adults in the six-site OCD Collaborative Genetics Study, diagnosed with lifetime OCD using semi-structured interviews. YBOCS-SC data were factor analyzed at both the individual item and symptom category levels. Factor score intraclass correlations were calculated using a subsample of 145 independent affected sib pairs. The item- and category-level factor analyses yielded nearly identical 5-factor solutions. While significant sib-sib associations were found for four of the five factors, Hoarding and Taboo Thoughts were the most robustly familial (r ICC>or=0.2). This report presents considerable converging evidence for a five-factor structural model of OCD symptoms, including separate factor analyses employing individual symptoms and symptom categories, as well as sibling concordance. The results support investigation of this multidimensional model in OCD genetic linkage studies.

Sex-specific clinical correlates of hoarding in obsessive-compulsive disorder

Little is known about whether the clinical correlates of hoarding behavior are different in men and women with obsessive-compulsive disorder (OCD). In the current study, we evaluated the association of hoarding with categories of obsessions and compulsions, psychiatric disorders, personality dimensions, and other clinical characteristics separately in 151 men and 358 women with OCD who were examined during the OCD Collaborative Genetics Study. We found that, among men but not women, hoarding was associated with aggressive, sexual, and religious obsessions and checking compulsions. In men, hoarding was associated with generalized anxiety disorder and tics whereas, among women, hoarding was associated with social phobia, post-traumatic stress disorder, body dysmorphic disorder, nail biting, and skin picking. In women but not men, hoarding was associated with schizotypal and dependent personality disorder dimensions, and with low conscientiousness. These findings indicate that specific clinical correlates of hoarding in OCD are different in men and women and may reflect sex-specific differences in the course, expression, and/or etiology of hoarding behavior in OCD.

Demographic and clinical characteristics associated with treatment status in family members with obsessive-compulsive disorder

This study investigated the demographic and clinical factors that influence treatment status in family members with obsessive compulsive disorder (OCD). Six hundred and two subjects from the OCD Collaborative Genetics Study were interviewed using the Structured Clinical Interview for DSM-IV (SCID) to diagnose Axis I disorders, and the Yale-Brown Obsessive Compulsive Scale (Y-BOCS) for assessment of OCD symptoms. The demographic and clinical data were compared between subjects who had received treatment and those who had not. A precipitous onset of symptoms, severe illness, multiple obsessions and compulsions, and co-morbid affective disorders were all positively associated with receiving treatment. Older age and the presence of obsessive compulsive personality disorder (OCPD) or OCPD traits were negatively associated with treatment. Gender and age at onset of symptoms did not predict treatment history. The mean duration from onset of symptoms to receiving treatment was 13.8+/-SD 11.9 years, but there was a direct relationship between current age and time to treatment, with younger subjects receiving treatment sooner. Clinical factors are predominant in predicting treatment status in family members with OCD. Although the mean duration from onset of symptoms to treatment was long, younger family members appear to receive treatment sooner.