S Purcell

Integrated genotype calling and association analysis of SNPs, common copy number polymorphisms and rare CNVs.

Accurate and complete measurement of single nucleotide (SNP) and copy number (CNV) variants, both common and rare, will be required to understand the role of genetic variation in disease. We present Birdsuite, a four-stage analytical framework instantiated in software for deriving integrated and mutually consistent copy number and SNP genotypes. The method sequentially assigns copy number across regions of common copy number polymorphisms (CNPs), calls genotypes of SNPs, identifies rare CNVs via a hidden Markov model (HMM), and generates an integrated sequence and copy number genotype at every locus (for example, including genotypes such as A-null, AAB and BBB in addition to AA, AB and BB calls). Such genotypes more accurately depict the underlying sequence of each individual, reducing the rate of apparent mendelian inconsistencies. The Birdsuite software is applied here to data from the Affymetrix SNP 6.0 array. Additionally, we describe a method, implemented in PLINK, to utilize these combined SNP and CNV genotypes for association testing with a phenotype.

Support for involvement of neuregulin 1 in schizophrenia pathophysiology

Schizophrenia is a common, multigenic psychiatric disorder. Linkage studies, including a recent meta-analysis of genome scans, have repeatedly implicated chromosome 8p12-p23.1 in schizophrenia susceptibility. More recently, significant association with a candidate gene on 8p12, neuregulin 1 (NRG1), has been reported in several European and Chinese samples. We investigated NRG1 for association in schizophrenia patients of Portuguese descent to determine whether this gene is a risk factor in this population. We tested NRG1 markers and haplotypes for association in 111 parent-proband trios, 321 unrelated cases, and 242 control individuals. Associations were found with a haplotype that overlaps the risk haplotype originally reported in the Icelandic population (‘HapICE’), and two haplotypes located in the 3′ end of NRG1 (all P<0.05). However, association was not detected with HapICE itself. Comparison of NRG1 transcript expression in peripheral leukocytes from schizophrenia patients and unaffected siblings identified 3.8-fold higher levels of the SMDF variant in patients (P=0.039). Significant positive correlations (P<0.001) were found between SMDF and HRG-beta 2 expression and between HRG-gamma and ndf43 expression, suggesting common transcriptional regulation of NRG1 variants. In summary, our results suggest that haplotypes across NRG1 and multiple NRG1 variants are involved in schizophrenia.

Parental Monitoring Moderates the Importance of Genetic and Environmental Influences on Adolescent Smoking

Although there is a substantial literature on the role of parenting in adolescent substance use, most parenting effects have been small in magnitude and studied outside the context of genetically informative designs, raising debate and controversy about the influence that parents have on their children (D. C. Rowe, 1994). Using a genetically informative twin-family design, the authors studied the role of parental monitoring on adolescent smoking at age 14. Although monitoring had only small main effects, consistent with the literature, there were dramatic moderation effects associated with parental monitoring: At high levels of parental monitoring, environmental influences were predominant in the etiology of adolescent smoking, but at low levels of parental monitoring, genetic influences assumed far greater importance. These analyses demonstrate that the etiology of adolescent smoking varies dramatically as a function of parenting.

Genome wide significant associations in schizophrenia to ITIH3/4, CACNA1C and SDCCAG8, and extensive replication of associations reported by the Schizophrenia PGC

The Schizophrenia Psychiatric Genome-Wide Association Study Consortium (PGC) highlighted 81 single-nucleotide polymorphisms (SNPs) with moderate evidence for association to schizophrenia. After follow-up in independent samples, seven loci attained genome-wide significance (GWS), but multi-locus tests suggested some SNPs that did not do so represented true associations. We tested 78 of the 81 SNPs in 2640 individuals with a clinical diagnosis of schizophrenia attending a clozapine clinic (CLOZUK), 2504 cases with a research diagnosis of bipolar disorder, and 2878 controls. In CLOZUK, we obtained significant replication to the PGC-associated allele for no fewer than 37 (47%) of the SNPs, including many prior GWS major histocompatibility complex (MHC) SNPs as well as 3/6 non-MHC SNPs for which we had data that were reported as GWS by the PGC. After combining the new schizophrenia data with those of the PGC, variants at three loci (ITIH3/4, CACNA1C and SDCCAG8) that had not previously been GWS in schizophrenia attained that level of support. In bipolar disorder, we also obtained significant evidence for association for 21% of the alleles that had been associated with schizophrenia in the PGC. Our study independently confirms association to three loci previously reported to be GWS in schizophrenia, and identifies the first GWS evidence in schizophrenia for a further three loci. Given the number of independent replications and the power of our sample, we estimate 98% (confidence interval (CI) 78–100%) of the original set of 78 SNPs represent true associations. We also provide strong evidence for overlap in genetic risk between schizophrenia and bipolar disorder.

Influence of RGS2 on Anxiety-Related Temperament, Personality, and Brain Function

CONTEXTnAlthough anxiety disorders are heritable, their genetic and phenotypic complexity has made the identification of susceptibility genes difficult. Well-validated animal models and intermediate phenotypes provide crucial tools for genetic dissection of anxiety. The gene encoding regulator of G protein signaling 2 (Rgs2) is a quantitative trait gene that influences mouse anxiety behavior, making its human ortholog (RGS2) a compelling candidate gene for human anxiety phenotypes.nnnOBJECTIVEnTo examine whether variation in RGS2 is associated with intermediate phenotypes for human anxiety disorders.nnnDESIGNnFamily-based and case-control association analysis of single-nucleotide polymorphisms at the RGS2 locus in 3 independent samples.nnnSETTINGnMassachusetts General Hospital, University of California, San Diego, and San Diego State University.nnnPARTICIPANTSnStudy participants included a family-based sample (n = 119 families) of children who underwent laboratory-based assessments of temperament (behavioral inhibition), a sample of 744 unrelated adults who completed assessments of extraversion and introversion, and 55 unrelated adults who underwent functional magnetic resonance imaging measures of response to emotional faces.nnnMAIN OUTCOME MEASURESnLaboratory-based behavioral measures of childhood temperament, self-report measure of personality, and functional magnetic resonance imaging response to emotion processing.nnnRESULTSnMarkers spanning RGS2 were associated with childhood behavioral inhibition, a temperamental precursor of social anxiety disorder (haplotype P = 3 x 10(-5); odds ratio, 2.99 in complete trios). In independent samples, RGS2 markers, including rs4606, which has previously been associated with RGS2 expression, were also associated with introversion (a core personality trait in social anxiety disorder) and with increased limbic activation (insular cortex and amygdala) during emotion processing (brain phenotypes correlated with social anxiety). The genotype at rs4606 explained 10% to 15% of the variance in amygdala and insular cortex activation to emotional faces.nnnCONCLUSIONSnThese results provide the first evidence that a gene that influences anxiety in mice is associated with intermediate phenotypes for human anxiety disorders across multiple levels of assessment, including childhood temperament, adult personality, and brain function. This translational research suggests that some genetic influences on anxiety are evolutionarily conserved and that pharmacologic modulation of RGS2 function may provide a novel therapeutic approach for anxiety disorders.

Replication of bipolar disorder susceptibility alleles and identification of two novel genome-wide significant associations in a new bipolar disorder case-control sample.

We have conducted a genotyping study using a custom Illumina Infinium HD genotyping array, the ImmunoChip, in a new UK sample of 1218 bipolar disorder (BD) cases and 2913 controls that have not been used in any studies previously reported independently or in meta-analyses. The ImmunoChip was designed before the publication of the Psychiatric Genome-Wide Association Study Consortium Bipolar Disorder Working Group (PGC-BD) meta-analysis data. As such 3106 single-nucleotide polymorphisms (SNPs) with a P-value <1 × 10−3 from the BD meta-analysis by Ferreira et al. were genotyped. We report support for two of the three most strongly associated chromosomal regions in the Ferreira study, CACNA1C (rs1006737, P=4.09 × 10−4) and 15q14 (rs2172835, P=0.043) but not ANK3 (rs10994336, P=0.912). We have combined our ImmunoChip data (569 quasi-independent SNPs from the 3016 SNPs genotyped) with the recently published PGC-BD meta-analysis data, using either the PGC-BD combined discovery and replication data where available or just the discovery data where the SNP was not typed in a replication sample in PGC-BD. Our data provide support for two regions, at ODZ4 and CACNA1C, with prior evidence for genome-wide significant (GWS) association in PGC-BD meta-analysis. In addition, the combined analysis shows two novel GWS associations. First, rs7296288 (P=8.97 × 10−9, odds ratio (OR)=0.9), an intergenic polymorphism on chromosome 12 located between RHEBL1 and DHH. Second, rs3818253 (P=3.88 × 10−8, OR=1.16), an intronic SNP on chromosome 20q11.2 in the gene TRPC4AP, which lies in a high linkage disequilibrium region along with the genes GSS and MYH7B.

GENESiS: creating a composite index of the vulnerability to anxiety and depression in a community-based sample of siblings

There is considerable evidence for a unitary and dimensional view of the genetic vulnerability to symptoms of anxiety and depression. The GENESiS (Genetic Environmental-Nature of Emotional States in Siblings) Study aims to use a multivariate approach to detect genetic loci that contribute to individual differences in this vulnerability dimension. The study used the UK General Practice Research Framework to generate a community-based sample of siblings. Questionnaire measures of anxiety/depression included the short form of the neuroticism scale from the revised Eysenck Personality Questionnaire (EPQ-N), the General Health Questionnaire (GHQ-12), and the anxious arousal and high positive affect subscales from the Mood and Anxiety Symptoms Questionnaire (MASQ-AA and MASQ-HPA). Genetic model-fitting of 2658 unselected sibships provided evidence for a single common genetic (familial) factor that accounted for a substantial proportion of the genetic variances and covariances of these four measures. Using the parameter estimates of this model, we constructed a composite index of this common genetic factor. This index, which has a sib correlation of 0.22, will be used as a quantitative phenotype in the molecular genetic phase of GENESiS.

The serotonin transporter gene as a QTL for ADHD.

Molecular studies of attention deficit hyperactivity disorder (ADHD) have identified susceptibility genes for the categorically diagnosed disorder using operational diagnostic criteria. Here, we take a QTL approach to mapping genes for ADHD using a composite continuous index of ADHD behavior in a large epidemiological sample. Previous studies of clinical ADHD suggest that two functional polymorphisms in the serotonin transporter gene (SLC6A4), one in the 5′‐regulatory region of the gene (5‐HTTLPR) and the other a VNTR (5‐HTTVNTR) in the second intron, as well as a single nucleotide polymorphism in the 3′‐untranslated region (3′‐UTR SNP), may be associated with the disorder. Here, we investigate these polymorphisms as well as an additional ten SNPs spread across the gene. We found significant association with the long (L) allele of the 5‐HTTLPR; Pu2009=u20090.019, but neither the 5‐HTTVNTR nor the 3′‐UTR SNP were significantly associated. Significant associations (Pu2009<u20090.05) were found for a further 5 the 10 other markers tested. We found evidence for two haplotype blocks spanning the region. We found strong evidence for association with the first haplotype block (comprised of four markers), with the significance of a combined primary and secondary test of association reaching an empirical P valueu2009=u20090.0054 for the global test and an empirical P valueu2009=u20090.00081 for the largest local test. Thus, we show here that SLC6A4, which has a major influence on brain serotonin availability, may be a QTL for ADHD.

Parental Phenotypes in Family-Based Association Analysis

Family-based association designs are popular, because they offer inherent control of population stratification based on age, sex, ethnicity, and environmental exposure. However, the efficiency of these designs is hampered by current analytic strategies that consider only offspring phenotypes. Here, we describe the incorporation of parental phenotypes and, specifically, the inclusion of parental genotype-phenotype correlation terms in association tests, providing a series of tests that effectively span an efficiency-robustness spectrum. The model is based on the between-within-sibship association model presented in 1999 by Fulker and colleagues for quantitative traits and extended here to nuclear families. By use of a liability-threshold-model approach, standard dichotomous and/or qualitative disease phenotypes can be analyzed (and can include appropriate corrections for phenotypically ascertained samples), which allows for the application of this model to analysis of the commonly used affected-proband trio design. We show that the incorporation of parental phenotypes can considerably increase power, as compared with the standard transmission/disequilibrium test and equivalent quantitative tests, while providing both significant protection against stratification and a means of evaluating the contribution of stratification to positive results. This methodology enables the extraction of more information from existing family-based collections that are currently being genotyped and analyzed by use of standard approaches.

Association at SYNE1 in both bipolar disorder and recurrent major depression

Genome-wide association studies (GWAS) have identified a number of loci that have strong support for their association with bipolar disorder (BD). The Psychiatric Genome-Wide Association Study (GWAS) Consortium Bipolar Disorder Working Group (PGC-BD) meta-analysis of BD GWAS data sets and replication samples identified evidence (P=6.7 × 10−7, odds ratio (OR)=1.147) of association with the risk of BD at the polymorphism rs9371601 within SYNE1, a gene which encodes nesprin-1. Here we have tested this polymorphism in an independent BD case (n=1527) and control (n=1579) samples, and find evidence for association (P=0.0095) with similar effect sizes to those previously observed in BD (allelic OR=1.148). In a combined (meta) analysis of PGC-BD data (both primary and replication data) and our independent BD samples, we found genome-wide significant evidence for association (P=2.9 × 10−8, OR=1.104). We have also examined the polymorphism in our recurrent unipolar depression cases (n=1159) and control (n=2592) sample, and found that the risk allele was associated with risk for recurrent major depression (P=0.032, OR=1.118). Our findings add to the evidence that association at this locus influences susceptibility to bipolar and unipolar mood disorders.