Rare tandem repeat expansions associate with genes involved in synaptic and neuronal signaling functions in schizophrenia.

Tandem repeat expansions (TREs) are associated with over 60 monogenic disorders and have recently been implicated in complex disorders such as cancer and autism spectrum disorder. The role of TREs in schizophrenia is now emerging. In this study, we have performed a genome-wide investigation of TREs in schizophrenia. Using genome sequence data from 1154 Swedish schizophrenia cases and 934 ancestry-matched population controls, we have detected genome-wide rare (<0.1% population frequency) TREs that have motifs with a length of 2-20 base pairs. We find that the proportion of individuals carrying rare TREs is significantly higher in the schizophrenia group. There is a significantly higher burden of rare TREs in schizophrenia cases than in controls in genic regions, particularly in postsynaptic genes, in genes overlapping brain expression quantitative trait loci, and in brain-expressed genes that are differentially expressed between schizophrenia cases and controls. We demonstrate that TRE-associated genes are more constrained and primarily impact synaptic and neuronal signaling functions. These results have been replicated in an independent Canadian sample that consisted of 252 schizophrenia cases of European ancestry and 222 ancestry-matched controls. Our results support the involvement of rare TREs in schizophrenia etiology.

The impact of educational attainment, intelligence and intellectual disability on schizophrenia: a Swedish population-based register and genetic study.

Schizophrenia (SCZ) is highly heterogenous and no subtypes characterizing treatment response or longitudinal course well. Cognitive impairment is a core clinical feature of SCZ and a determinant of poorer outcome. Genetic overlap between SCZ and cognitive traits is complex, with limited studies of comprehensive epidemiological and genomic evidence. To examine the relation between SCZ and three cognitive traits, educational attainment (EDU), premorbid cognitive ability, and intellectual disability (ID), we used two Swedish samples: a national cohort (14,230 SCZ cases and 3,816,264 controls) and a subsample with comprehensive genetic data (4992 cases and 6009 controls). Population-based analyses confirmed worse cognition as a risk factor for SCZ, and the pedigree and SNP-based genetic correlations were comparable. In the genotyped cases, those with high EDU and premorbid cognitive ability tended to have higher polygenetic risk scores (PRS) of EDU and intelligence and fewer rare exonic variants. Finally, by applying an empirical clustering method, we dissected SCZ cases into four replicable subgroups characterized by EDU and ID. In particular, the subgroup with higher EDU in the national cohort had fewer adverse outcomes including long hospitalization and death. In the genotyped subsample, this subgroup had higher PRS of EDU and no excess of rare genetic burdens than controls. In conclusion, we found extensive evidence of a robust relation between cognitive traits and SCZ, underscoring the importance of cognition in dissecting the heterogeneity of SCZ.

Elevated common variant genetic risk for tourette syndrome in a densely-affected pedigree.

Tourette syndrome (TS) is a highly heritable neuropsychiatric disorder with complex patterns of genetic inheritance. Recent genetic findings in TS have highlighted both numerous common variants with small effects and a few rare variants with moderate or large effects. Here we searched for genetic causes of TS in a large, densely-affected British pedigree using a systematic genomic approach. This pedigree spans six generations and includes 122 members, 85 of whom were individually interviewed, and 53 of whom were diagnosed as "cases" (consisting of 28 with definite or probable TS, 20 with chronic multiple tics [CMT], and five with obsessive-compulsive behaviors [OCB]). A total of 66 DNA samples were available (25 TS, 15 CMT, 4 OCB cases, and 22 unaffecteds) and all were genotyped using a dense single nucleotide polymorphism (SNP) array to identify shared segments, copy number variants (CNVs), and to calculate genetic risk scores. Eight cases were also whole genome sequenced to test whether any rare variants were shared identical by descent. While we did not identify any notable CNVs, single nucleotide variants, indels or repeat expansions of near-Mendelian effect, the most distinctive feature of this family proved to be an unusually high load of common risk alleles for TS. We found that cases within this family carried a higher load of TS common variant risk similar to that previously found in unrelated TS cases. Thus far, the strongest evidence from genetic data for contribution to TS risk in this family comes from multiple common risk variants rather than one or a few variants of strong effect.

Correction: Common-variant associations with fragile X syndrome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Association test using Copy Number Profile Curves (CONCUR) enhances power in rare copy number variant analysis.

Copy number variants (CNVs) are the gain or loss of DNA segments in the genome that can vary in dosage and length. CNVs comprise a large proportion of variation in human genomes and impact health conditions. To detect rare CNV associations, kernel-based methods have been shown to be a powerful tool due to their flexibility in modeling the aggregate CNV effects, their ability to capture effects from different CNV features, and their accommodation of effect heterogeneity. To perform a kernel association test, a CNV locus needs to be defined so that locus-specific effects can be retained during aggregation. However, CNV loci are arbitrarily defined and different locus definitions can lead to different performance depending on the underlying effect patterns. In this work, we develop a new kernel-based test called CONCUR (i.e., copy number profile curve-based association test) that is free from a definition of locus and evaluates CNV-phenotype associations by comparing individuals' copy number profiles across the genomic regions. CONCUR is built on the proposed concepts of "copy number profile curves" to describe the CNV profile of an individual, and the "common area under the curve (cAUC) kernel" to model the multi-feature CNV effects. The proposed method captures the effects of CNV dosage and length, accounts for the numerical nature of copy numbers, and accommodates between- and within-locus etiological heterogeneity without the need to define artificial CNV loci as required in current kernel methods. In a variety of simulation settings, CONCUR shows comparable or improved power over existing approaches. Real data analyses suggest that CONCUR is well powered to detect CNV effects in the Swedish Schizophrenia Study and the Taiwan Biobank.

Common-variant associations with fragile X syndrome.

Fragile X syndrome is rare but a prominent cause of intellectual disability. It is usually caused by a de novo mutation that occurs on multiple haplotypes and thus would not be expected to be detectible using genome-wide association (GWA). We conducted GWA in 89 male FXS cases and 266 male controls, and detected multiple genome-wide significant signals near FMR1 (odds ratio = 8.10, P = 2.5 × 10-10). These findings withstood robust attempts at falsification. Fine-mapping yielded a minimum P = 1.13 × 10-14, but did not narrow the interval. Comprehensive functional genomic integration did not provide a mechanistic hypothesis. Controls carrying a risk haplotype had significantly longer FMR1 CGG repeats than controls with the protective haplotype (P = 4.75 × 10-5), which may predispose toward increases in CGG number to the premutation range over many generations. This is a salutary reminder of the complexity of even "simple" monogenetic disorders.

A randomized approach to speed up the analysis of large-scale read-count data in the application of CNV detection.

BACKGROUND: The application of high-throughput sequencing in a broad range of quantitative genomic assays (e.g., DNA-seq, ChIP-seq) has created a high demand for the analysis of large-scale read-count data. Typically, the genome is divided into tiling windows and windowed read-count data is generated for the entire genome from which genomic signals are detected (e.g. copy number changes in DNA-seq, enrichment peaks in ChIP-seq). For accurate analysis of read-count data, many state-of-the-art statistical methods use generalized linear models (GLM) coupled with the negative-binomial (NB) distribution by leveraging its ability for simultaneous bias correction and signal detection. However, although statistically powerful, the GLM+NB method has a quadratic computational complexity and therefore suffers from slow running time when applied to large-scale windowed read-count data. In this study, we aimed to speed up substantially the GLM+NB method by using a randomized algorithm and we demonstrate here the utility of our approach in the application of detecting copy number variants (CNVs) using a real example. RESULTS: We propose an efficient estimator, the randomized GLM+NB coefficients estimator (RGE), for speeding up the GLM+NB method. RGE samples the read-count data and solves the estimation problem on a smaller scale. We first theoretically validated the consistency and the variance properties of RGE. We then applied RGE to GENSENG, a GLM+NB based method for detecting CNVs. We named the resulting method as "R-GENSENG". Based on extensive evaluation using both simulated and empirical data, we concluded that R-GENSENG is ten times faster than the original GENSENG while maintaining GENSENG's accuracy in CNV detection. CONCLUSIONS: Our results suggest that RGE strategy developed here could be applied to other GLM+NB based read-count analyses, i.e. ChIP-seq data analysis, to substantially improve their computational efficiency while preserving the analytic power.

A New Method for Detecting Associations with Rare Copy-Number Variants.

Copy number variants (CNVs) play an important role in the etiology of many diseases such as cancers and psychiatric disorders. Due to a modest marginal effect size or the rarity of the CNVs, collapsing rare CNVs together and collectively evaluating their effect serves as a key approach to evaluating the collective effect of rare CNVs on disease risk. While a plethora of powerful collapsing methods are available for sequence variants (e.g., SNPs) in association analysis, these methods cannot be directly applied to rare CNVs due to the CNV-specific challenges, i.e., the multi-faceted nature of CNV polymorphisms (e.g., CNVs vary in size, type, dosage, and details of gene disruption), and etiological heterogeneity (e.g., heterogeneous effects of duplications and deletions that occur within a locus or in different loci). Existing CNV collapsing analysis methods (a.k.a. the burden test) tend to have suboptimal performance due to the fact that these methods often ignore heterogeneity and evaluate only the marginal effects of a CNV feature. We introduce CCRET, a random effects test for collapsing rare CNVs when searching for disease associations. CCRET is applicable to variants measured on a multi-categorical scale, collectively modeling the effects of multiple CNV features, and is robust to etiological heterogeneity. Multiple confounders can be simultaneously corrected. To evaluate the performance of CCRET, we conducted extensive simulations and analyzed large-scale schizophrenia datasets. We show that CCRET has powerful and robust performance under multiple types of etiological heterogeneity, and has performance comparable to or better than existing methods when there is no heterogeneity.

Implication of LRRC4C and DPP6 in neurodevelopmental disorders.

We performed whole-genome sequencing on an individual from a family with variable psychiatric phenotypes that had a sensory processing disorder, apraxia, and autism. The proband harbored a maternally inherited balanced translocation (46,XY,t(11;14)(p12;p12)mat) that disrupted LRRC4C, a member of the highly specialized netrin G family of axon guidance molecules. The proband also inherited a paternally derived chromosomal inversion that disrupted DPP6, a potassium channel interacting protein. Copy Number (CN) analysis in 14,077 cases with neurodevelopmental disorders and 8,960 control subjects revealed that 60% of cases with exonic deletions in LRRC4C had a second clinically recognizable syndrome associated with variable clinical phenotypes, including 16p11.2, 1q44, and 2q33.1 CN syndromes, suggesting LRRC4C deletion variants may be modifiers of neurodevelopmental disorders. In vitro, functional assessments modeling patient deletions in LRRC4C suggest a negative regulatory role of these exons found in the untranslated region of LRRC4C, which has a single, terminal coding exon. These data suggest that the proband's autism may be due to the inheritance of disruptions in both DPP6 and LRRC4C, and may highlight the importance of the netrin G family and potassium channel interacting molecules in neurodevelopmental disorders. © 2016 Wiley Periodicals, Inc.

Allele-specific copy-number discovery from whole-genome and whole-exome sequencing.

Copy-number variants (CNVs) are a major form of genetic variation and a risk factor for various human diseases, so it is crucial to accurately detect and characterize them. It is conceivable that allele-specific reads from high-throughput sequencing data could be leveraged to both enhance CNV detection and produce allele-specific copy number (ASCN) calls. Although statistical methods have been developed to detect CNVs using whole-genome sequence (WGS) and/or whole-exome sequence (WES) data, information from allele-specific read counts has not yet been adequately exploited. In this paper, we develop an integrated method, called AS-GENSENG, which incorporates allele-specific read counts in CNV detection and estimates ASCN using either WGS or WES data. To evaluate the performance of AS-GENSENG, we conducted extensive simulations, generated empirical data using existing WGS and WES data sets and validated predicted CNVs using an independent methodology. We conclude that AS-GENSENG not only predicts accurate ASCN calls but also improves the accuracy of total copy number calls, owing to its unique ability to exploit information from both total and allele-specific read counts while accounting for various experimental biases in sequence data. Our novel, user-friendly and computationally efficient method and a complete analytic protocol is freely available at https://sourceforge.net/projects/asgenseng/.

CNV analysis in a large schizophrenia sample implicates deletions at 16p12.1 and SLC1A1 and duplications at 1p36.33 and CGNL1.

Large and rare copy number variants (CNVs) at several loci have been shown to increase risk for schizophrenia. Aiming to discover novel susceptibility CNV loci, we analyzed 6882 cases and 11 255 controls genotyped on Illumina arrays, most of which have not been used for this purpose before. We identified genes enriched for rare exonic CNVs among cases, and then attempted to replicate the findings in additional 14 568 cases and 15 274 controls. In a combined analysis of all samples, 12 distinct loci were enriched among cases with nominal levels of significance (P < 0.05); however, none would survive correction for multiple testing. These loci include recurrent deletions at 16p12.1, a locus previously associated with neurodevelopmental disorders (P = 0.0084 in the discovery sample and P = 0.023 in the replication sample). Other plausible candidates include non-recurrent deletions at the glutamate transporter gene SLC1A1, a CNV locus recently suggested to be involved in schizophrenia through linkage analysis, and duplications at 1p36.33 and CGNL1. A burden analysis of large (>500 kb), rare CNVs showed a 1.2% excess in cases after excluding known schizophrenia-associated loci, suggesting that additional susceptibility loci exist. However, even larger samples are required for their discovery.

An inherited duplication at the gene p21 Protein-Activated Kinase 7 (PAK7) is a risk factor for psychosis.

Identifying rare, highly penetrant risk mutations may be an important step in dissecting the molecular etiology of schizophrenia. We conducted a gene-based analysis of large (>100 kb), rare copy-number variants (CNVs) in the Wellcome Trust Case Control Consortium 2 (WTCCC2) schizophrenia sample of 1564 cases and 1748 controls all from Ireland, and further extended the analysis to include an additional 5196 UK controls. We found association with duplications at chr20p12.2 (P = 0.007) and evidence of replication in large independent European schizophrenia (P = 0.052) and UK bipolar disorder case-control cohorts (P = 0.047). A combined analysis of Irish/UK subjects including additional psychosis cases (schizophrenia and bipolar disorder) identified 22 carriers in 11 707 cases and 10 carriers in 21 204 controls [meta-analysis Cochran-Mantel-Haenszel P-value = 2 × 10(-4); odds ratio (OR) = 11.3, 95% CI = 3.7, ∞]. Nineteen of the 22 cases and 8 of the 10 controls carried duplications starting at 9.68 Mb with similar breakpoints across samples. By haplotype analysis and sequencing, we identified a tandem ~149 kb duplication overlapping the gene p21 Protein-Activated Kinase 7 (PAK7, also called PAK5) which was in linkage disequilibrium with local haplotypes (P = 2.5 × 10(-21)), indicative of a single ancestral duplication event. We confirmed the breakpoints in 8/8 carriers tested and found co-segregation of the duplication with illness in two additional family members of one of the affected probands. We demonstrate that PAK7 is developmentally co-expressed with another known psychosis risk gene (DISC1) suggesting a potential molecular mechanism involving aberrant synapse development and plasticity.

One CNV Discordance in NRXN1 Observed Upon Genome-wide Screening in 38 Pairs of Adult Healthy Monozygotic Twins.

Monozygotic (MZ) twins stem from the same single fertilized egg and therefore share all their inherited genetic variation. This is one of the unequivocal facts on which genetic epidemiology and twin studies are based. To what extent this also implies that MZ twins share genotypes in adult tissues is not precisely established, but a common pragmatic assumption is that MZ twins are 100% genetically identical also in adult tissues. During the past decade, this view has been challenged by several reports, with observations of differences in post-zygotic copy number variations (CNVs) between members of the same MZ pair. In this study, we performed a systematic search for differences of CNVs within 38 adult MZ pairs who had been misclassified as dizygotic (DZ) twins by questionnaire-based assessment. Initial scoring by PennCNV suggested a total of 967 CNV discordances. The within-pair correlation in number of CNVs detected was strongly dependent on confidence score filtering and reached a plateau of r = 0.8 when restricting to CNVs detected with confidence score larger than 50. The top-ranked discordances were subsequently selected for validation by quantitative polymerase chain reaction (qPCR), from which one single ~120kb deletion in NRXN1 on chromosome 2 (bp 51017111-51136802) was validated. Despite involving an exon, no sign of cognitive/mental consequences was apparent in the affected twin pair, potentially reflecting limited or lack of expression of the transcripts containing this exon in nerve/brain.

Improving detection of copy-number variation by simultaneous bias correction and read-depth segmentation.

Structural variation is an important class of genetic variation in mammals. High-throughput sequencing (HTS) technologies promise to revolutionize copy-number variation (CNV) detection but present substantial analytic challenges. Converging evidence suggests that multiple types of CNV-informative data (e.g. read-depth, read-pair, split-read) need be considered, and that sophisticated methods are needed for more accurate CNV detection. We observed that various sources of experimental biases in HTS confound read-depth estimation, and note that bias correction has not been adequately addressed by existing methods. We present a novel read-depth-based method, GENSENG, which uses a hidden Markov model and negative binomial regression framework to identify regions of discrete copy-number changes while simultaneously accounting for the effects of multiple confounders. Based on extensive calibration using multiple HTS data sets, we conclude that our method outperforms existing read-depth-based CNV detection algorithms. The concept of simultaneous bias correction and CNV detection can serve as a basis for combining read-depth with other types of information such as read-pair or split-read in a single analysis. A user-friendly and computationally efficient implementation of our method is freely available.

DeepGWAS: Enhance GWAS Signals for Neuropsychiatric Disorders via Deep Neural Network.

Genetic dissection of neuropsychiatric disorders can potentially reveal novel therapeutic targets. While genome-wide association studies (GWAS) have tremendously advanced our understanding, we approach a sample size bottleneck (i.e., the number of cases needed to identify >90% of all loci is impractical). Therefore, computationally enhancing GWAS on existing samples may be particularly valuable. Here, we describe DeepGWAS, a deep neural network-based method to enhance GWAS by integrating GWAS results with linkage disequilibrium and brain-related functional annotations. DeepGWAS enhanced schizophrenia (SCZ) loci by ~3X when applied to the largest European GWAS, and 21.3% enhanced loci were validated by the latest multi-ancestry GWAS. Importantly, DeepGWAS models can be transferred to other neuropsychiatric disorders. Transferring SCZ-trained models to Alzheimer's disease and major depressive disorder, we observed 1.3-17.6X detected loci compared to standard GWAS, among which 27-40% were validated by other GWAS studies. We anticipate DeepGWAS to be a powerful tool in GWAS studies.

Increased Prevalence of Rare Copy Number Variants in Treatment-Resistant Psychosis.

BACKGROUND: It remains unknown why ~30% of patients with psychotic disorders fail to respond to treatment. Previous genomic investigations of treatment-resistant psychosis have been inconclusive, but some evidence suggests a possible link between rare disease-associated copy number variants (CNVs) and worse clinical outcomes in schizophrenia. Here, we identified schizophrenia-associated CNVs in patients with treatment-resistant psychotic symptoms and then compared the prevalence of these CNVs to previously published schizophrenia cases not selected for treatment resistance. METHODS: CNVs were identified using chromosomal microarray (CMA) and whole exome sequencing (WES) in 509 patients with treatment-resistant psychosis (a lack of clinical response to ≥3 adequate antipsychotic medication trials over at least 5 years of psychiatric hospitalization). Prevalence of schizophrenia-associated CNVs in this sample was compared to that in a previously published large schizophrenia cohort study. RESULTS: Integrating CMA and WES data, we identified 47 cases (9.2%) with at least one CNV of known or possible neuropsychiatric risk. 4.7% (n = 24) carried a known neurodevelopmental risk CNV. The prevalence of well-replicated schizophrenia-associated CNVs was 4.1%, with duplications of the 16p11.2 and 15q11.2-q13.1 regions, and deletions of the 22q11.2 chromosomal region as the most frequent CNVs. Pairwise loci-based analysis identified duplications of 15q11.2-q13.1 to be independently associated with treatment resistance. CONCLUSIONS: These findings suggest that CNVs may uniquely impact clinical phenotypes beyond increasing risk for schizophrenia and may potentially serve as biological entry points for studying treatment resistance. Further investigation will be necessary to elucidate the spectrum of phenotypic characteristics observed in adult psychiatric patients with disease-associated CNVs.

Leveraging Base Pair Mammalian Constraint to Understand Genetic Variation and Human Disease.

Although thousands of genomic regions have been associated with heritable human diseases, attempts to elucidate biological mechanisms are impeded by a general inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function that is agnostic to cell type or disease mechanism. Here, single base phyloP scores from the whole genome alignment of 240 placental mammals identified 3.5% of the human genome as significantly constrained, and likely functional. We compared these scores to large-scale genome annotation, genome-wide association studies (GWAS), copy number variation, clinical genetics findings, and cancer data sets. Evolutionarily constrained positions are enriched for variants explaining common disease heritability (more than any other functional annotation). Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Leveraging base-pair mammalian constraint to understand genetic variation and human disease.

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Genome-wide association mapping of loci for antipsychotic-induced extrapyramidal symptoms in mice.

Tardive dyskinesia (TD) is a debilitating, unpredictable, and often irreversible side effect resulting from chronic treatment with typical antipsychotic agents such as haloperidol. TD is characterized by repetitive, involuntary, purposeless movements primarily of the orofacial region. In order to investigate genetic susceptibility to TD, we used a validated mouse model for a systems genetics analysis geared toward detecting genetic predictors of TD in human patients. Phenotypic data from 27 inbred strains chronically treated with haloperidol and phenotyped for vacuous chewing movements were subject to a comprehensive genomic analysis involving 426,493 SNPs, 4,047 CNVs, brain gene expression, along with gene network and bioinformatic analysis. Our results identified ~50 genes that we expect to have high prior probabilities for association with haloperidol-induced TD, most of which have never been tested for association with human TD. Among our top candidates were genes regulating the development of brain motor control regions (Zic4 and Nkx6-1), glutamate receptors (Grin1 and Grin2a), and an indirect target of haloperidol (Drd1a) that has not been studied as well as the direct target, Drd2.

A customized and versatile high-density genotyping array for the mouse.

We designed a high-density mouse genotyping array containing 623,124 single-nucleotide polymorphisms that captures the known genetic variation present in the laboratory mouse. The array also contains 916,269 invariant genomic probes targeted to functional elements and regions known to harbor segmental duplications. The array opens the door to the characterization of genetic diversity, copy-number variation, allele-specific gene expression and DNA methylation, and will extend the successes of human genome-wide association studies to the mouse.