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.

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.

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.

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/.

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.

The recombinational anatomy of a mouse chromosome.

Among mammals, genetic recombination occurs at highly delimited sites known as recombination hotspots. They are typically 1-2 kb long and vary as much as a 1,000-fold or more in recombination activity. Although much is known about the molecular details of the recombination process itself, the factors determining the location and relative activity of hotspots are poorly understood. To further our understanding, we have collected and mapped the locations of 5,472 crossover events along mouse Chromosome 1 arising in 6,028 meioses of male and female reciprocal F1 hybrids of C57BL/6J and CAST/EiJ mice. Crossovers were mapped to a minimum resolution of 225 kb, and those in the telomere-proximal 24.7 Mb were further mapped to resolve individual hotspots. Recombination rates were evolutionarily conserved on a regional scale, but not at the local level. There was a clear negative-exponential relationship between the relative activity and abundance of hotspot activity classes, such that a small number of the most active hotspots account for the majority of recombination. Females had 1.2x higher overall recombination than males did, although the sex ratio showed considerable regional variation. Locally, entirely sex-specific hotspots were rare. The initiation of recombination at the most active hotspot was regulated independently on the two parental chromatids, and analysis of reciprocal crosses indicated that parental imprinting has subtle effects on recombination rates. It appears that the regulation of mammalian recombination is a complex, dynamic process involving multiple factors reflecting species, sex, individual variation within species, and the properties of individual hotspots.

Association of candidate genes with phenotypic traits relevant to anorexia nervosa.

This analysis is a follow-up to an earlier investigation of 182 genes selected as likely candidate genetic variations conferring susceptibility to anorexia nervosa (AN). As those initial case-control results revealed no statistically significant differences in single nucleotide polymorphisms, herein, we investigate alternative phenotypes associated with AN. In 1762 females, using regression analyses, we examined the following: (i) lowest illness-related attained body mass index; (ii) age at menarche; (iii) drive for thinness; (iv) body dissatisfaction; (v) trait anxiety; (vi) concern over mistakes; and (vii) the anticipatory worry and pessimism versus uninhibited optimism subscale of the harm avoidance scale. After controlling for multiple comparisons, no statistically significant results emerged. Although results must be viewed in the context of limitations of statistical power, the approach illustrates a means of potentially identifying genetic variants conferring susceptibility to AN because less complex phenotypes associated with AN are more proximal to the genotype and may be influenced by fewer genes.

Crossover interference underlies sex differences in recombination rates.

In many organisms, recombination rates differ between the two sexes. Here we show that in mice, this is because of a shorter genomic interference distance in females than in males, measured in Mb. However, the interference distance is the same in terms of bivalent length. We propose a model in which the interference distance in the two sexes reflects the compaction of chromosomes at the pachytene stage of meiosis.

A survey of airway responsiveness in 36 inbred mouse strains facilitates gene mapping studies and identification of quantitative trait loci.

Airway hyper-responsiveness (AHR) is a critical phenotype of human asthma and animal models of asthma. Other studies have measured AHR in nine mouse strains, but only six strains have been used to identify genetic loci underlying AHR. Our goals were to increase the genetic diversity of available strains by surveying 27 additional strains, to apply haplotype association mapping to the 36-strain survey, and to identify new genetic determinants for AHR. We derived AHR from the increase in airway resistance in females subjected to increasing levels of methacholine concentrations. We used haplotype association mapping to identify associations between AHR and haplotypes on chromosomes 3, 5, 8, 12, 13, and 14. And we used bioinformatics techniques to narrow the identified region on chromosome 13, reducing the region to 29 candidate genes, with 11 of considerable interest. Our combined use of haplotype association mapping with bioinformatics tools is the first study of its kind for AHR on these 36 strains of mice. Our analyses have narrowed the possible QTL genes and will facilitate the discovery of novel genes that regulate AHR in mice.

Increased burden of ultra-rare structural variants localizing to boundaries of topologically associated domains in schizophrenia.

Despite considerable progress in schizophrenia genetics, most findings have been for large rare structural variants and common variants in well-imputed regions with few genes implicated from exome sequencing. Whole genome sequencing (WGS) can potentially provide a more complete enumeration of etiological genetic variation apart from the exome and regions of high linkage disequilibrium. We analyze high-coverage WGS data from 1162 Swedish schizophrenia cases and 936 ancestry-matched population controls. Our main objective is to evaluate the contribution to schizophrenia etiology from a variety of genetic variants accessible to WGS but not by previous technologies. Our results suggest that ultra-rare structural variants that affect the boundaries of topologically associated domains (TADs) increase risk for schizophrenia. Alterations in TAD boundaries may lead to dysregulation of gene expression. Future mechanistic studies will be needed to determine the precise functional effects of these variants on biology.

Characterization of Single Gene Copy Number Variants in Schizophrenia.

BACKGROUND: Genetic studies of schizophrenia have implicated numerous risk loci including several copy number variants (CNVs) of large effect and hundreds of loci of small effect. In only a few cases has a specific gene been clearly identified. Rare CNVs affecting a single gene offer a potential avenue to discovering schizophrenia risk genes. METHODS: CNVs were generated from exome sequencing of 4913 schizophrenia cases and 6188 control subjects from Sweden. We integrated two CNV calling methods (XHMM and ExomeDepth) to expand our set of single-gene CNVs and leveraged two different approaches for validating these variants (quantitative polymerase chain reaction and NanoString). RESULTS: We found a significant excess of all rare CNVs (deletions: p = .0004, duplications: p = .0006) and single-gene CNVs (deletions: p = .04, duplications: p = .03) in schizophrenia cases compared with control subjects. An expanded set of CNVs generated from integrating multiple approaches showed a significant burden of deletions in 11 of 21 gene sets previously implicated in schizophrenia and across all genes in those sets (p = .008), although no tests survived correction. We performed an extensive validation of all deletions in the significant set of voltage-gated calcium channels among CNVs called from both exome sequencing and genotyping arrays. In total, 4 exonic, single-gene deletions were validated in schizophrenia cases and none in control subjects (p = .039), of which all were identified by exome sequencing. CONCLUSIONS: These results point to the potential contribution of single-gene CNVs to schizophrenia, indicate that the utility of exome sequencing for CNV calling has yet to be maximized, and note that single-gene CNVs should be included in gene-focused studies using other classes of variation.

Otitis media: a genome-wide linkage scan with evidence of susceptibility loci within the 17q12 and 10q22.3 regions.

BACKGROUND: Otitis media (OM) is a common worldwide pediatric health care problem that is known to be influenced by genetics. The objective of our study was to use linkage analysis to map possible OM susceptibility genes. METHODS: Using a stringent diagnostic model in which only those who underwent tympanostomy tube insertion at least once for recurrent/persistent OM are considered affected, we have carried out a genome-wide linkage scan using the 10K Affymetrix SNP panel. We genotyped 403 Caucasian families containing 1,431 genotyped individuals and 377 genotyped affected sib pairs, and 26 African American families containing 75 genotyped individuals and 27 genotyped affected sib pairs. After careful quality control, non-parametric linkage analysis was carried out using 8,802 SNPs. RESULTS: In the Caucasian-only data set, our most significant linkage peak is on chromosome 17q12 at rs226088 with a p-value of 0.00007. Other peaks of potential interest are on 10q22.3 (0.00181 at rs1878001), 7q33 (0.00105 at rs958408), 6p25.1 (0.00261 at rs554653), and 4p15.2 (0.00301 at rs2133507). In the combined Caucasian and African American dataset, the 10q22.3 peak becomes more significant, with a minimal p-value of 0.00026 at rs719871. Family-based association testing reveals signals near previously implicated genes: 513 kb from SFTPA2 (10q22.3), 48 kb from IFNG (12q14), and 870 kb from TNF (6p21.3). CONCLUSION: Our scan does not provide evidence for linkage in the previously reported regions of 10q26.3 and 19q13.43. Our best-supported linkage regions may contain susceptibility genes that influence the risk for recurrent/persistent OM. Plausible candidates in 17q12 include AP2B1, CCL5, and a cluster of other CCL genes, and in 10q22.3, SFTPA2.

Exploration of large, rare copy number variants associated with psychiatric and neurodevelopmental disorders in individuals with anorexia nervosa.

Anorexia nervosa (AN) is a serious and heritable psychiatric disorder. To date, studies of copy number variants (CNVs) have been limited and inconclusive because of small sample sizes. We conducted a case-only genome-wide CNV survey in 1983 female AN cases included in the Genetic Consortium for Anorexia Nervosa. Following stringent quality control procedures, we investigated whether pathogenic CNVs in regions previously implicated in psychiatric and neurodevelopmental disorders were present in AN cases. We observed two instances of the well-established pathogenic CNVs in AN cases. In addition, one case had a deletion in the 13q12 region, overlapping with a deletion reported previously in two AN cases. As a secondary aim, we also examined our sample for CNVs over 1 Mbp in size. Out of the 40 instances of such large CNVs that were not implicated previously for AN or neuropsychiatric phenotypes, two of them contained genes with previous neuropsychiatric associations, and only five of them had no associated reports in public CNV databases. Although ours is the largest study of its kind in AN, larger datasets are needed to comprehensively assess the role of CNVs in the etiology of AN.

Defective carbohydrate metabolism in mice homozygous for the tubby mutation.

Tub is a member of a small gene family, the tubby-like proteins (TULPs), with predominant expression in neurons. Mice carrying a mutation in Tub develop retinal and cochlear degeneration as well as late-onset obesity with insulin resistance. During behavioral and metabolic testing, we found that homozygous C57BL/6J-Tub(tub) mice have a lower respiratory quotient than C57BL/6J controls before the onset of obesity, indicating that tubby homozygotes fail to activate carbohydrate metabolism and instead rely on fat metabolism for energy needs. In concordance with this, tubby mice show higher excretion of ketone bodies and accumulation of glycogen in the liver. Quantitation of liver mRNA levels shows that, during the transition from light to dark period, tubby mice fail to induce glucose-6-phosphate dehydrogenase (G6pdh), the rate-limiting enzyme in the pentose phosphate pathway that normally supplies NADPH for de novo fatty acid synthesis and glutathione reduction. Reduced G6PDH protein levels and enzymatic activity in tubby mice lead accordingly to lower levels of NADPH and reduced glutathione (GSH), respectively. mRNA levels for the lipolytic enzymes acetyl-CoA synthetase and carnitine palmitoyltransferase are increased during the dark cycle and decreased during the light period, and several citric acid cycle genes are dysregulated in tubby mice. Examination of hypothalamic gene expression showed high levels of preproorexin mRNA leading to accumulation of orexin peptide in the lateral hypothalamus. We hypothesize that abnormal hypothalamic orexin expression leads to changes in liver carbohydrate metabolism and may contribute to the moderate obesity observed in tubby mice.

Clozapine-induced agranulocytosis is associated with rare HLA-DQB1 and HLA-B alleles.

Clozapine is a particularly effective antipsychotic medication but its use is curtailed by the risk of clozapine-induced agranulocytosis/granulocytopenia (CIAG), a severe adverse drug reaction occurring in up to 1% of treated individuals. Identifying genetic risk factors for CIAG could enable safer and more widespread use of clozapine. Here we perform the largest and most comprehensive genetic study of CIAG to date by interrogating 163 cases using genome-wide genotyping and whole-exome sequencing. We find that two loci in the major histocompatibility complex are independently associated with CIAG: a single amino acid in HLA-DQB1 (126Q) (P=4.7 × 10(-14), odds ratio (OR)=0.19, 95% confidence interval (CI)=0.12-0.29) and an amino acid change in the extracellular binding pocket of HLA-B (158T) (P=6.4 × 10(-10), OR=3.3, 95% CI=2.3-4.9). These associations dovetail with the roles of these genes in immunogenetic phenotypes and adverse drug responses for other medications, and provide insight into the pathophysiology of CIAG.

Endogenous retinoids in the pathogenesis of alopecia areata.

Alopecia areata (AA) is an autoimmune disease that attacks anagen hair follicles. Gene array in graft-induced C3H/HeJ mice revealed that genes involved in retinoic acid (RA) synthesis were increased, whereas RA degradation genes were decreased in AA compared with sham controls. This was confirmed by immunohistochemistry in biopsies from patients with AA and both mouse and rat AA models. RA levels were also increased in C3H/HeJ mice with AA. C3H/HeJ mice were fed a purified diet containing one of the four levels of dietary vitamin A or an unpurified diet 2 weeks before grafting and disease progression followed. High vitamin A accelerated AA, whereas mice that were not fed vitamin A had more severe disease by the end of the study. More hair follicles were in anagen in mice fed high vitamin A. Both the number and localization of granzyme B-positive cells were altered by vitamin A. IFNγ was also the lowest and IL13 highest in mice fed high vitamin A. Other cytokines were reduced and chemokines increased as the disease progressed, but no additional effects of vitamin A were seen. Combined, these results suggest that vitamin A regulates both the hair cycle and immune response to alter the progression of AA.

CGDSNPdb: a database resource for error-checked and imputed mouse SNPs.

The Center for Genome Dynamics Single Nucleotide Polymorphism Database (CGDSNPdb) is an open-source value-added database with more than nine million mouse single nucleotide polymorphisms (SNPs), drawn from multiple sources, with genotypes assigned to multiple inbred strains of laboratory mice. All SNPs are checked for accuracy and annotated for properties specific to the SNP as well as those implied by changes to overlapping protein-coding genes. CGDSNPdb serves as the primary interface to two unique data sets, the 'imputed genotype resource' in which a Hidden Markov Model was used to assess local haplotypes and the most probable base assignment at several million genomic loci in tens of strains of mice, and the Affymetrix Mouse Diversity Genotyping Array, a high density microarray with over 600,000 SNPs and over 900,000 invariant genomic probes. CGDSNPdb is accessible online through either a web-based query tool or a MySQL public login. Database URL: http://cgd.jax.org/cgdsnpdb/

Patterns of recombination activity on mouse chromosome 11 revealed by high resolution mapping.

The success of high resolution genetic mapping of disease predisposition and quantitative trait loci in humans and experimental animals depends on the positions of key crossover events around the gene of interest. In mammals, the majority of recombination occurs at highly delimited 1-2 kb long sites known as recombination hotspots, whose locations and activities are distributed unevenly along the chromosomes and are tightly regulated in a sex specific manner. The factors determining the location of hotspots started to emerge with the finding of PRDM9 as a major hotspot regulator in mammals, however, additional factors modulating hotspot activity and sex specificity are yet to be defined. To address this limitation, we have collected and mapped the locations of 4829 crossover events occurring on mouse chromosome 11 in 5858 meioses of male and female reciprocal F1 hybrids of C57BL/6J and CAST/EiJ mice. This chromosome was chosen for its medium size and high gene density and provided a comparison with our previous analysis of recombination on the longest mouse chromosome 1. Crossovers were mapped to an average resolution of 127 kb, and thirteen hotspots were mapped to <8 kb. Most crossovers occurred in a small number of the most active hotspots. Females had higher recombination rate than males as a consequence of differences in crossover interference and regional variation of sex specific rates along the chromosome. Comparison with chromosome 1 showed that recombination events tend to be positioned in similar fashion along the centromere-telomere axis but independently of the local gene density. It appears that mammalian recombination is regulated on at least three levels, chromosome-wide, regional, and at individual hotspots, and these regulation levels are influenced by sex and genetic background but not by gene content.