Genetic influences on eight psychiatric disorders based on family data of 4 408 646 full and half-siblings, and genetic data of 333 748 cases and controls.

BACKGROUND: Most studies underline the contribution of heritable factors for psychiatric disorders. However, heritability estimates depend on the population under study, diagnostic instruments, and study designs that each has its inherent assumptions, strengths, and biases. We aim to test the homogeneity in heritability estimates between two powerful, and state of the art study designs for eight psychiatric disorders. METHODS: We assessed heritability based on data of Swedish siblings (N = 4 408 646 full and maternal half-siblings), and based on summary data of eight samples with measured genotypes (N = 125 533 cases and 208 215 controls). All data were based on standard diagnostic criteria. Eight psychiatric disorders were studied: (1) alcohol dependence (AD), (2) anorexia nervosa, (3) attention deficit/hyperactivity disorder (ADHD), (4) autism spectrum disorder, (5) bipolar disorder, (6) major depressive disorder, (7) obsessive-compulsive disorder (OCD), and (8) schizophrenia. RESULTS: Heritability estimates from sibling data varied from 0.30 for Major Depression to 0.80 for ADHD. The estimates based on the measured genotypes were lower, ranging from 0.10 for AD to 0.28 for OCD, but were significant, and correlated positively (0.19) with national sibling-based estimates. When removing OCD from the data the correlation increased to 0.50. CONCLUSIONS: Given the unique character of each study design, the convergent findings for these eight psychiatric conditions suggest that heritability estimates are robust across different methods. The findings also highlight large differences in genetic and environmental influences between psychiatric disorders, providing future directions for etiological psychiatric research.

Genome-wide common and rare variant analysis provides novel insights into clozapine-associated neutropenia.

This corrects the article DOI: 10.1038/mp.2016.97.

ASD and schizophrenia show distinct developmental profiles in common genetic overlap with population-based social communication difficulties.

Difficulties in social communication are part of the phenotypic overlap between autism spectrum disorders (ASD) and schizophrenia. Both conditions follow, however, distinct developmental patterns. Symptoms of ASD typically occur during early childhood, whereas most symptoms characteristic of schizophrenia do not appear before early adulthood. We investigated whether overlap in common genetic influences between these clinical conditions and impairments in social communication depends on the developmental stage of the assessed trait. Social communication difficulties were measured in typically-developing youth (Avon Longitudinal Study of Parents and Children, N⩽5553, longitudinal assessments at 8, 11, 14 and 17 years) using the Social Communication Disorder Checklist. Data on clinical ASD (PGC-ASD: 5305 cases, 5305 pseudo-controls; iPSYCH-ASD: 7783 cases, 11 359 controls) and schizophrenia (PGC-SCZ2: 34 241 cases, 45 604 controls, 1235 trios) were either obtained through the Psychiatric Genomics Consortium (PGC) or the Danish iPSYCH project. Overlap in genetic influences between ASD and social communication difficulties during development decreased with age, both in the PGC-ASD and the iPSYCH-ASD sample. Genetic overlap between schizophrenia and social communication difficulties, by contrast, persisted across age, as observed within two independent PGC-SCZ2 subsamples, and showed an increase in magnitude for traits assessed during later adolescence. ASD- and schizophrenia-related polygenic effects were unrelated to each other and changes in trait-disorder links reflect the heterogeneity of genetic factors influencing social communication difficulties during childhood versus later adolescence. Thus, both clinical ASD and schizophrenia share some genetic influences with impairments in social communication, but reveal distinct developmental profiles in their genetic links, consistent with the onset of clinical symptoms.

Polygenic risk for schizophrenia and neurocognitive performance in patients with schizophrenia.

Both neurocognitive deficits and schizophrenia are highly heritable. Genetic overlap between neurocognitive deficits and schizophrenia has been observed in both the general population and in the clinical samples. This study aimed to examine if the polygenic architecture of susceptibility to schizophrenia modified neurocognitive performance in schizophrenia patients. Schizophrenia polygenic risk scores (PRSs) were first derived from the Psychiatric Genomics Consortium (PGC) on schizophrenia, and then the scores were calculated in our independent sample of 1130 schizophrenia trios, who had PsychChip data and were part of the Schizophrenia Families from Taiwan project. Pseudocontrols generated from the nontransmitted parental alleles of the parents in these trios were compared with alleles in schizophrenia patients in assessing the replicability of PGC-derived susceptibility variants. Schizophrenia PRS at the P-value threshold (PT) of 0.1 explained 0.2% in the variance of disease status in this Han-Taiwanese samples, and the score itself had a P-value 0.05 for the association test with the disorder. Each patient underwent neurocognitive evaluation on sustained attention using the continuous performance test and executive function using the Wisconsin Card Sorting Test. We applied a structural equation model to construct the neurocognitive latent variable estimated from multiple measured indices in these 2 tests, and then tested the association between the PRS and the neurocognitive latent variable. Higher schizophrenia PRS generated at the PT of 0.1 was significantly associated with poorer neurocognitive performance with explained variance 0.5%. Our findings indicated that schizophrenia susceptibility variants modify the neurocognitive performance in schizophrenia patients.

The iPSYCH2012 case-cohort sample: new directions for unravelling genetic and environmental architectures of severe mental disorders.

The Integrative Psychiatric Research (iPSYCH) consortium has established a large Danish population-based Case-Cohort sample (iPSYCH2012) aimed at unravelling the genetic and environmental architecture of severe mental disorders. The iPSYCH2012 sample is nested within the entire Danish population born between 1981 and 2005, including 1 472 762 persons. This paper introduces the iPSYCH2012 sample and outlines key future research directions. Cases were identified as persons with schizophrenia (N=3540), autism (N=16 146), attention-deficit/hyperactivity disorder (N=18 726) and affective disorder (N=26 380), of which 1928 had bipolar affective disorder. Controls were randomly sampled individuals (N=30 000). Within the sample of 86 189 individuals, a total of 57 377 individuals had at least one major mental disorder. DNA was extracted from the neonatal dried blood spot samples obtained from the Danish Neonatal Screening Biobank and genotyped using the Illumina PsychChip. Genotyping was successful for 90% of the sample. The assessments of exome sequencing, methylation profiling, metabolome profiling, vitamin-D, inflammatory and neurotrophic factors are in progress. For each individual, the iPSYCH2012 sample also includes longitudinal information on health, prescribed medicine, social and socioeconomic information, and analogous information among relatives. To the best of our knowledge, the iPSYCH2012 sample is the largest and most comprehensive data source for the combined study of genetic and environmental aetiologies of severe mental disorders.

Genome-wide common and rare variant analysis provides novel insights into clozapine-associated neutropenia.

The antipsychotic clozapine is uniquely effective in the management of schizophrenia; however, its use is limited by its potential to induce agranulocytosis. The causes of this, and of its precursor neutropenia, are largely unknown, although genetic factors have an important role. We sought risk alleles for clozapine-associated neutropenia in a sample of 66 cases and 5583 clozapine-treated controls, through a genome-wide association study (GWAS), imputed human leukocyte antigen (HLA) alleles, exome array and copy-number variation (CNV) analyses. We then combined associated variants in a meta-analysis with data from the Clozapine-Induced Agranulocytosis Consortium (up to 163 cases and 7970 controls). In the largest combined sample to date, we identified a novel association with rs149104283 (odds ratio (OR)=4.32, P=1.79 × 10-8), intronic to transcripts of SLCO1B3 and SLCO1B7, members of a family of hepatic transporter genes previously implicated in adverse drug reactions including simvastatin-induced myopathy and docetaxel-induced neutropenia. Exome array analysis identified gene-wide associations of uncommon non-synonymous variants within UBAP2 and STARD9. We additionally provide independent replication of a previously identified variant in HLA-DQB1 (OR=15.6, P=0.015, positive predictive value=35.1%). These results implicate biological pathways through which clozapine may act to cause this serious adverse effect.

Genome-wide common and rare variant analysis provides novel insights into clozapine-associated neutropenia.

This corrects the article DOI: 10.1038/mp.2016.97.

Exome arrays capture polygenic rare variant contributions to schizophrenia.

Schizophrenia is a highly heritable disorder. Genome-wide association studies based largely on common alleles have identified over 100 schizophrenia risk loci, but it is also evident from studies of copy number variants (CNVs) and from exome-sequencing studies that rare alleles are also involved. Full characterization of the contribution of rare alleles to the disorder awaits the deployment of sequencing technology in very large sample sizes, meanwhile, as an interim measure, exome arrays allow rare non-synonymous variants to be sampled at a fraction of the cost. In an analysis of exome array data from 13 688 individuals (5585 cases and 8103 controls) from the UK, we found that rare (minor allele frequency < 0.1%) variant association signal was enriched among genes that map to autosomal loci that are genome-wide significant (GWS) in common variant studies of schizophrenia genome-wide association study (PGWAS = 0.01) as well as gene sets known to be enriched for rare variants in sequencing studies (PRARE = 0.026). We also identified the gene-wise equivalent of GWS support for WDR88 (WD repeat-containing protein 88), a gene of unknown function (P = 6.5 × 10(-7)). Rare alleles represented on exome chip arrays contribute to the genetic architecture of schizophrenia, but as is the case for GWAS, very large studies are required to reveal additional susceptibility alleles for the disorder.

Genome-wide autozygosity is associated with lower general cognitive ability.

Inbreeding depression refers to lower fitness among offspring of genetic relatives. This reduced fitness is caused by the inheritance of two identical chromosomal segments (autozygosity) across the genome, which may expose the effects of (partially) recessive deleterious mutations. Even among outbred populations, autozygosity can occur to varying degrees due to cryptic relatedness between parents. Using dense genome-wide single-nucleotide polymorphism (SNP) data, we examined the degree to which autozygosity associated with measured cognitive ability in an unselected sample of 4854 participants of European ancestry. We used runs of homozygosity-multiple homozygous SNPs in a row-to estimate autozygous tracts across the genome. We found that increased levels of autozygosity predicted lower general cognitive ability, and estimate a drop of 0.6 s.d. among the offspring of first cousins (P=0.003-0.02 depending on the model). This effect came predominantly from long and rare autozygous tracts, which theory predicts as more likely to be deleterious than short and common tracts. Association mapping of autozygous tracts did not reveal any specific regions that were predictive beyond chance after correcting for multiple testing genome wide. The observed effect size is consistent with studies of cognitive decline among offspring of known consanguineous relationships. These findings suggest a role for multiple recessive or partially recessive alleles in general cognitive ability, and that alleles decreasing general cognitive ability have been selected against over evolutionary time.

Analysis of exome sequence in 604 trios for recessive genotypes in schizophrenia.

Genetic associations involving both rare and common alleles have been reported for schizophrenia but there have been no systematic scans for rare recessive genotypes using fully phased trio data. Here, we use exome sequencing in 604 schizophrenia proband-parent trios to investigate the role of recessive (homozygous or compound heterozygous) nonsynonymous genotypes in the disorder. The burden of recessive genotypes was not significantly increased in probands at either a genome-wide level or in any individual gene after adjustment for multiple testing. At a system level, probands had an excess of nonsynonymous compound heterozygous genotypes (minor allele frequency, MAF ⩽ 1%) in voltage-gated sodium channels (VGSCs; eight in probands and none in parents, P = 1.5 × 10(-)(4)). Previous findings of multiple de novo loss-of-function mutations in this gene family, particularly SCN2A, in autism and intellectual disability provide biological and genetic plausibility for this finding. Pointing further to the involvement of VGSCs in schizophrenia, we found that these genes were enriched for nonsynonymous mutations (MAF ⩽ 0.1%) in cases genotyped using an exome array, (5585 schizophrenia cases and 8103 controls), and that in the trios data, synaptic proteins interacting with VGSCs were also enriched for both compound heterozygosity (P = 0.018) and de novo mutations (P = 0.04). However, we were unable to replicate the specific association with compound heterozygosity at VGSCs in an independent sample of Taiwanese schizophrenia trios (N = 614). We conclude that recessive genotypes do not appear to make a substantial contribution to schizophrenia at a genome-wide level. Although multiple lines of evidence, including several from this study, suggest that rare mutations in VGSCs contribute to the disorder, in the absence of replication of the original findings regarding compound heterozygosity, this conclusion requires evaluation in a larger sample of trios.

The genetic architecture of pediatric cognitive abilities in the Philadelphia Neurodevelopmental Cohort.

The objective of this analysis was to examine the genetic architecture of diverse cognitive abilities in children and adolescents, including the magnitude of common genetic effects and patterns of shared and unique genetic influences. Subjects included 3689 members of the Philadelphia Neurodevelopmental Cohort, a general population sample comprising those aged 8-21 years who completed an extensive battery of cognitive tests. We used genome-wide complex trait analysis to estimate the SNP-based heritability of each domain, as well as the genetic correlation between all domains that showed significant genetic influence. Several of the individual domains suggested strong influence of common genetic variants (for example, reading ability, h(2)g=0.43, P=4e-06; emotion identification, h(2)g=0.36, P=1e-05; verbal memory, h(2)g=0.24, P=0.005). The genetic correlations highlighted trait domains that are candidates for joint interrogation in future genetic studies (for example, language reasoning and spatial reasoning, r(g)=0.72, P=0.007). These results can be used to structure future genetic and neuropsychiatric investigations of diverse cognitive abilities.

Copy number variation in schizophrenia in Sweden.

Schizophrenia (SCZ) is a highly heritable neuropsychiatric disorder of complex genetic etiology. Previous genome-wide surveys have revealed a greater burden of large, rare copy number variations (CNVs) in SCZ cases and identified multiple rare recurrent CNVs that increase risk of SCZ although with incomplete penetrance and pleiotropic effects. Identification of additional recurrent CNVs and biological pathways enriched for SCZ CNVs requires greater sample sizes. We conducted a genome-wide survey for CNVs associated with SCZ using a Swedish national sample (4719 cases and 5917 controls). High-confidence CNV calls were generated using genotyping array intensity data, and their effect on risk of SCZ was measured. Our data confirm increased burden of large, rare CNVs in SCZ cases as well as significant associations for recurrent 16p11.2 duplications, 22q11.2 deletions and 3q29 deletions. We report a novel association for 17q12 duplications (odds ratio=4.16, P=0.018), previously associated with autism and mental retardation but not SCZ. Intriguingly, gene set association analyses implicate biological pathways previously associated with SCZ through common variation and exome sequencing (calcium channel signaling and binding partners of the fragile X mental retardation protein). We found significantly increased burden of the largest CNVs (>500 kb) in genes present in the postsynaptic density, in genomic regions implicated via SCZ genome-wide association studies and in gene products localized to mitochondria and cytoplasm. Our findings suggest that multiple lines of genomic inquiry--genome-wide screens for CNVs, common variation and exonic variation--are converging on similar sets of pathways and/or genes.

Detecting large copy number variants using exome genotyping arrays in a large Swedish schizophrenia sample.

Although copy number variants (CNVs) are important in genomic medicine, CNVs have not been systematically assessed for many complex traits. Several large rare CNVs increase risk for schizophrenia (SCZ) and autism and often demonstrate pleiotropic effects; however, their frequencies in the general population and other complex traits are unknown. Genotyping large numbers of samples is essential for progress. Large cohorts from many different diseases are being genotyped using exome-focused arrays designed to detect uncommon or rare protein-altering sequence variation. Although these arrays were not designed for CNV detection, the hybridization intensity data generated in each experiment could, in principle, be used for gene-focused CNV analysis. Our goal was to evaluate the extent to which CNVs can be detected using data from one particular exome array (the Illumina Human Exome Bead Chip). We genotyped 9100 Swedish subjects (3962 cases with SCZ and 5138 controls) using both standard genome-wide association study (GWAS) and exome arrays. In comparison with CNVs detected using GWAS arrays, we observed high sensitivity and specificity for detecting genic CNVs 400 kb including known pathogenic CNVs along with replicating the literature finding that cases with SCZ had greater enrichment for genic CNVs. Our data confirm the association of SCZ with 16p11.2 duplications and 22q11.2 deletions, and suggest a novel association with deletions at 11q12.2. Our results suggest the utility of exome-focused arrays in surveying large genic CNVs in very large samples; and thereby open the door for new opportunities such as conducting well-powered CNV assessment and comparisons between different diseases. The use of a single platform also minimizes potential confounding factors that could impact accurate detection.

Genome-wide association study of obsessive-compulsive disorder.

Obsessive-compulsive disorder (OCD) is a common, debilitating neuropsychiatric illness with complex genetic etiology. The International OCD Foundation Genetics Collaborative (IOCDF-GC) is a multi-national collaboration established to discover the genetic variation predisposing to OCD. A set of individuals affected with DSM-IV OCD, a subset of their parents, and unselected controls, were genotyped with several different Illumina SNP microarrays. After extensive data cleaning, 1465 cases, 5557 ancestry-matched controls and 400 complete trios remained, with a common set of 469,410 autosomal and 9657 X-chromosome single nucleotide polymorphisms (SNPs). Ancestry-stratified case-control association analyses were conducted for three genetically-defined subpopulations and combined in two meta-analyses, with and without the trio-based analysis. In the case-control analysis, the lowest two P-values were located within DLGAP1 (P=2.49 × 10(-6) and P=3.44 × 10(-6)), a member of the neuronal postsynaptic density complex. In the trio analysis, rs6131295, near BTBD3, exceeded the genome-wide significance threshold with a P-value=3.84 × 10(-8). However, when trios were meta-analyzed with the case-control samples, the P-value for this variant was 3.62 × 10(-5), losing genome-wide significance. Although no SNPs were identified to be associated with OCD at a genome-wide significant level in the combined trio-case-control sample, a significant enrichment of methylation QTLs (P<0.001) and frontal lobe expression quantitative trait loci (eQTLs) (P=0.001) was observed within the top-ranked SNPs (P<0.01) from the trio-case-control analysis, suggesting these top signals may have a broad role in gene expression in the brain, and possibly in the etiology of OCD.

Genome-wide association study of Tourette's syndrome.

Tourette's syndrome (TS) is a developmental disorder that has one of the highest familial recurrence rates among neuropsychiatric diseases with complex inheritance. However, the identification of definitive TS susceptibility genes remains elusive. Here, we report the first genome-wide association study (GWAS) of TS in 1285 cases and 4964 ancestry-matched controls of European ancestry, including two European-derived population isolates, Ashkenazi Jews from North America and Israel and French Canadians from Quebec, Canada. In a primary meta-analysis of GWAS data from these European ancestry samples, no markers achieved a genome-wide threshold of significance (P<5 × 10(-8)); the top signal was found in rs7868992 on chromosome 9q32 within COL27A1 (P=1.85 × 10(-6)). A secondary analysis including an additional 211 cases and 285 controls from two closely related Latin American population isolates from the Central Valley of Costa Rica and Antioquia, Colombia also identified rs7868992 as the top signal (P=3.6 × 10(-7) for the combined sample of 1496 cases and 5249 controls following imputation with 1000 Genomes data). This study lays the groundwork for the eventual identification of common TS susceptibility variants in larger cohorts and helps to provide a more complete understanding of the full genetic architecture of this disorder.

Common genetic variants differentially influence the transition from clinically defined states of fasting glucose metabolism.

AIMS/HYPOTHESIS: Common genetic variants have been associated with type 2 diabetes. We hypothesised that a subset of these variants may have different effects on the transition from normal fasting glucose (NFG) to impaired fasting glucose (IFG) than on that from IFG to diabetes. METHODS: We identified 16 type 2 diabetes risk variants from the Illumina Broad Candidate-gene Association Resource (CARe) array genotyped in 26,576 CARe participants. Participants were categorised at baseline as NFG, IFG or type 2 diabetic (n = 16,465, 8,017 or 2,291, respectively). Using Cox proportional hazards and likelihood ratio tests (LRTs), we compared rates of progression by genotype for 4,909 (NFG to IFG) and 1,518 (IFG to type 2 diabetes) individuals, respectively. We then performed multinomial regression analyses at baseline, comparing the risk of assignment to the NFG, IFG or diabetes groups by genotype. RESULTS: The rate of progression from NFG to IFG was significantly greater in participants carrying the risk allele at MTNR1B (p = 1 × 10(-4)), nominally greater at GCK and SLC30A8 (p < 0.05) and nominally smaller at IGF2BP2 (p = 0.01) than the rate of progression from IFG to diabetes by the LRT. Results of the baseline, multinomial regression model were consistent with these findings. CONCLUSIONS/INTERPRETATION: Common genetic risk variants at GCK, SLC30A8, IGF2BP2 and MTNR1B influence to different extents the development of IFG and the transition from IFG to type 2 diabetes. Our findings may have implications for understanding the genetic contribution of these variants to the development of IFG and type 2 diabetes.

Risk-based fluconazole prophylaxis of Candida bloodstream infection in a medical intensive care unit.

Candida bloodstream infection (CBSI) accounted for 50% of bloodstream infections in our medical intensive care unit (MICU) in 2004. Our objective was to evaluate a risk-based fluconazole prophylaxis program. CBSI incidence, patient demographics, and unit metrics were retrospectively reviewed for 2004. Starting on January 2005, patients meeting pre-specified criteria were placed on risk-based fluconazole prophylaxis and their outcomes, adverse events, and unit metrics were prospectively collected. The inclusion criteria were based on a clinical prediction rule and included an MICU stay greater than 72 h, broad-spectrum antibiotics, and central venous catheter, along with at least two of the following: mechanical ventilation for at least 48 h, any type of dialysis, parenteral nutrition, pancreatitis, systemic steroids, or other systemic immunosuppressive agents. For 2004, the unit had nine CBSI, corresponding to a rate of 3.4 CBSI/1,000 line-days. Four cases were caused by C. albicans, four by C. glabrata, and one by C. tropicalis. The mean +/- standard deviation (SD) APACHE II score for these patients was 25 +/- 9. In 2005, a total of 36 patients (2.6% of all unit admissions) received prophylaxis and the unit had two CBSI, corresponding to a rate of 0.79 CBSI/1,000 line-days. One patient had C. albicans and the other had C. tropicalis. The mean +/- SD APACHE II score for these patients was 21 +/- 8. The mean +/- SD duration of fluconazole prophylaxis was 8 +/- 6 days. Fluconazole was discontinued in two patients due to non-severe adverse events (acute eosinophilia, elevated transaminases). The attributable cost of CBSI in the unit in 2004 was $63,000 per episode. The total cost for the 36 courses of fluconazole was $6,000. When comparing the 2004 CBSI patients and the 2005 prophylaxis patients, we found similar acuity, demographics, and risk factors, with no differences in MICU or hospital mortality or length of stay. Risk-based fluconazole prophylaxis in an MICU with a high incidence of CBSI was safe and cost-effective when applied to a limited number of patients and produced a significant decrease in the incidence of this disease.

A high-density SNP linkage scan with 142 combined subtype ADHD sib pairs identifies linkage regions on chromosomes 9 and 16.

As part of the International Multi-centre ADHD Genetics project we completed an affected sibling pair study of 142 narrowly defined Diagnostic and Statistical Manual of Mental Disorders, fourth edition combined type attention deficit hyperactivity disorder (ADHD) proband-sibling pairs. No linkage was observed on the most established ADHD-linked genomic regions of 5p and 17p. We found suggestive linkage signals on chromosomes 9 and 16, respectively, with the highest multipoint nonparametric linkage signal on chromosome 16q23 at 99 cM (log of the odds, LOD=3.1) overlapping data published from the previous UCLA (University of California, Los Angeles) (LOD>1, approximately 95 cM) and Dutch (LOD>1, approximately 100 cM) studies. The second highest peak in this study was on chromosome 9q22 at 90 cM (LOD=2.13); both the previous UCLA and German studies also found some evidence of linkage at almost the same location (UCLA LOD=1.45 at 93 cM; German LOD=0.68 at 100 cM). The overlap of these two main peaks with previous findings suggests that loci linked to ADHD may lie within these regions. Meta-analysis or reanalysis of the raw data of all the available ADHD linkage scan data may help to clarify whether these represent true linked loci.