A burden of rare copy number variants in obsessive-compulsive disorder.

Current genetic research on obsessive-compulsive disorder (OCD) supports contributions to risk specifically from common single nucleotide variants (SNVs), along with rare coding SNVs and small insertion-deletions (indels). The contribution to OCD risk from rare copy number variants (CNVs), however, has not been formally assessed at a similar scale. Here we describe an analysis of rare CNVs called from genotype array data in 2248 deeply phenotyped OCD cases and 3608 unaffected controls from Sweden and Norway. Cases carry an elevated burden of CNVs ≥30 kb in size (OR = 1.12, P = 1.77 × 10-3). The excess rate of these CNVs in cases versus controls was around 0.07 (95% CI 0.02-0.11, P = 2.58 × 10-3). This signal was largely driven by CNVs overlapping protein-coding regions (OR = 1.19, P = 3.08 × 10-4), particularly deletions impacting loss-of-function intolerant genes (pLI >0.995, OR = 4.12, P = 2.54 × 10-5). We did not identify any specific locus where CNV burden was associated with OCD case status at genome-wide significance, but we noted non-random recurrence of CNV deletions in cases (permutation P = 2.60 × 10-3). In cases where sufficient clinical data were available (n = 1612) we found that carriers of neurodevelopmental duplications were more likely to have comorbid autism (P < 0.001), and that carriers of deletions overlapping neurodevelopmental genes had lower treatment response (P = 0.02). The results demonstrate a contribution of rare CNVs to OCD risk, and suggest that studies of rare coding variation in OCD would have increased power to identify risk genes if this class of variation were incorporated into formal tests.

Genome-wide association study identifies new loci associated with OCD.

To date, four genome-wide association studies (GWAS) of obsessive-compulsive disorder (OCD) have been published, reporting a high single-nucleotide polymorphism (SNP)-heritability of 28% but finding only one significant SNP. A substantial increase in sample size will likely lead to further identification of SNPs, genes, and biological pathways mediating the susceptibility to OCD. We conducted a GWAS meta-analysis with a 2-3-fold increase in case sample size (OCD cases: N = 37,015, controls: N = 948,616) compared to the last OCD GWAS, including six previously published cohorts (OCGAS, IOCDF-GC, IOCDF-GC-trio, NORDiC-nor, NORDiC-swe, and iPSYCH) and unpublished self-report data from 23andMe Inc. We explored the genetic architecture of OCD by conducting gene-based tests, tissue and celltype enrichment analyses, and estimating heritability and genetic correlations with 74 phenotypes. To examine a potential heterogeneity in our data, we conducted multivariable GWASs with MTAG. We found support for 15 independent genome-wide significant loci (14 new) and 79 protein-coding genes. Tissue enrichment analyses implicate multiple cortical regions, the amygdala, and hypothalamus, while cell type analyses yielded 12 cell types linked to OCD (all neurons). The SNP-based heritability of OCD was estimated to be 0.08. Using MTAG we found evidence for specific genetic underpinnings characteristic of different cohort-ascertainment and identified additional significant SNPs. OCD was genetically correlated with 40 disorders or traits-positively with all psychiatric disorders and negatively with BMI, age at first birth and multiple autoimmune diseases. The GWAS meta-analysis identified several biologically informative genes as important contributors to the aetiology of OCD. Overall, we have begun laying the groundwork through which the biology of OCD will be understood and described.

A Burden of Rare Copy Number Variants in Obsessive-Compulsive Disorder.

Current genetic research on obsessive-compulsive disorder (OCD) supports contributions to risk specifically from common single nucleotide variants (SNVs), along with rare coding SNVs and small insertion-deletions (indels). The contribution to OCD risk from large, rare copy number variants (CNVs), however, has not been formally assessed at a similar scale. Here we describe an analysis of rare CNVs called from genotype array data in 2,248 deeply phenotyped OCD cases and 3,608 unaffected controls from Sweden and Norway. We found that in general cases carry an elevated burden of large (>30kb, at least 15 probes) CNVs (OR=1.12, P=1.77×10-3). The excess rate of these CNVs in cases versus controls was around 0.07 (95% CI 0.02-0.11, P=2.58×10-3). This signal was largely driven by CNVs overlapping protein-coding regions (OR=1.19, P=3.08×10-4), particularly deletions impacting loss-of-function intolerant genes (pLI>0.995, OR=4.12, P=2.54×10-5). We did not identify any specific locus where CNV burden was associated with OCD case status at genome-wide significance, but we noted non-random recurrence of CNV deletions in cases (permutation P = 2.60×10-3). In cases where sufficient clinical data were available (n=1612) we found that carriers of neurodevelopmental duplications were more likely to have comorbid autism (P<0.001), and that carriers of deletions overlapping neurodevelopmental genes had lower treatment response (P=0.02). The results demonstrate a contribution of large, rare CNVs to OCD risk, and suggest that studies of rare coding variation in OCD would have increased power to identify risk genes if this class of variation were incorporated into formal tests.

Association of autism diagnosis and polygenic scores with eating disorder severity.

Among individuals with eating disorders (ED), those with co-occurring autism are often considered to have more severe presentations and poorer prognosis. However, previous findings have been contradictory and limited by small sample size and/or cross-sectional assessment of autistic traits. We examine the hypothesis that autism diagnosis and autism polygenic score (PGS) are associated with increased ED severity in a large ED cohort using a broad range of ED severity indicators. Our cohort included 3189 individuals (64 males) born 1977-2000 with current or previous anorexia nervosa who participated in the Anorexia Nervosa Genetics Initiative-Sweden (ANGI-SE) and for whom genotypes and linkage to national registers were available. We identified 134 (4.2%) individuals with registered autism diagnoses. Individuals with confirmed autism diagnosis had significantly more severe ED across three sets of severity indicators. Some of the largest effects were found for the proportion of individuals who attempted suicide and who received tube feeding (higher in autism), and for the time spent in inpatient care (longer in autism). Results for autism PGS were not statistically significant. Adapting ED treatment to the needs of individuals with co-occurring autism is an important research direction to improve treatment outcome in this group.

First genome-wide association study of esophageal atresia identifies three genetic risk loci at CTNNA3, FOXF1/FOXC2/FOXL1, and HNF1B.

Esophageal atresia with or without tracheoesophageal fistula (EA/TEF) is the most common congenital malformation of the upper digestive tract. This study represents the first genome-wide association study (GWAS) to identify risk loci for EA/TEF. We used a European case-control sample comprising 764 EA/TEF patients and 5,778 controls and observed genome-wide significant associations at three loci. On chromosome 10q21 within the gene CTNNA3 (p = 2.11 × 10-8; odds ratio [OR] = 3.94; 95% confidence interval [CI], 3.10-5.00), on chromosome 16q24 next to the FOX gene cluster (p = 2.25 × 10-10; OR = 1.47; 95% CI, 1.38-1.55) and on chromosome 17q12 next to the gene HNF1B (p = 3.35 × 10-16; OR = 1.75; 95% CI, 1.64-1.87). We next carried out an esophageal/tracheal transcriptome profiling in rat embryos at four selected embryonic time points. Based on these data and on already published data, the implicated genes at all three GWAS loci are promising candidates for EA/TEF development. We also analyzed the genetic EA/TEF architecture beyond the single marker level, which revealed an estimated single-nucleotide polymorphism (SNP)-based heritability of around 37% ± 14% standard deviation. In addition, we examined the polygenicity of EA/TEF and found that EA/TEF is less polygenic than other complex genetic diseases. In conclusion, the results of our study contribute to a better understanding on the underlying genetic architecture of ET/TEF with the identification of three risk loci and candidate genes.

The Genetic Architecture of Obsessive-Compulsive Disorder: Contribution of Liability to OCD From Alleles Across the Frequency Spectrum.

OBJECTIVE: Obsessive-compulsive disorder (OCD) is known to be substantially heritable; however, the contribution of genetic variation across the allele frequency spectrum to this heritability remains uncertain. The authors used two new homogeneous cohorts to estimate the heritability of OCD from inherited genetic variation and contrasted the results with those of previous studies. METHODS: The sample consisted of 2,090 Swedish-born individuals diagnosed with OCD and 4,567 control subjects, all genotyped for common genetic variants, specifically >400,000 single-nucleotide polymorphisms (SNPs) with minor allele frequency (MAF) ≥0.01. Using genotypes of these SNPs to estimate distant familial relationships among individuals, the authors estimated the heritability of OCD, both overall and partitioned according to MAF bins. RESULTS: Narrow-sense heritability of OCD was estimated at 29% (SE=4%). The estimate was robust, varying only modestly under different models. Contrary to an earlier study, however, SNPs with MAF between 0.01 and 0.05 accounted for 10% of heritability, and estimated heritability per MAF bin roughly followed expectations based on a simple model for SNP-based heritability. CONCLUSIONS: These results indicate that common inherited risk variation (MAF ≥0.01) accounts for most of the heritable variation in OCD. SNPs with low MAF contribute meaningfully to the heritability of OCD, and the results are consistent with expectation under the "infinitesimal model" (also referred to as the "polygenic model"), where risk is influenced by a large number of loci across the genome and across MAF bins.

Study protocol of comprehensive risk evaluation for anorexia nervosa in twins (CREAT): a study of discordant monozygotic twins with anorexia nervosa.

BACKGROUND: Anorexia nervosa (AN) is a severe disorder, for which genetic evidence suggests psychiatric as well as metabolic origins. AN has high somatic and psychiatric comorbidities, broad impact on quality of life, and elevated mortality. Risk factor studies of AN have focused on differences between acutely ill and recovered individuals. Such comparisons often yield ambiguous conclusions, as alterations could reflect different effects depending on the comparison. Whereas differences found in acutely ill patients could reflect state effects that are due to acute starvation or acute disease-specific factors, they could also reflect underlying traits. Observations in recovered individuals could reflect either an underlying trait or a "scar" due to lasting effects of sustained undernutrition and illness. The co-twin control design (i.e., monozygotic [MZ] twins who are discordant for AN and MZ concordant control twin pairs) affords at least partial disambiguation of these effects. METHODS: Comprehensive Risk Evaluation for Anorexia nervosa in Twins (CREAT) will be the largest and most comprehensive investigation of twins who are discordant for AN to date. CREAT utilizes a co-twin control design that includes endocrinological, neurocognitive, neuroimaging, genomic, and multi-omic approaches coupled with an experimental component that explores the impact of an overnight fast on most measured parameters. DISCUSSION: The multimodal longitudinal twin assessment of the CREAT study will help to disambiguate state, trait, and "scar" effects, and thereby enable a deeper understanding of the contribution of genetics, epigenetics, cognitive functions, brain structure and function, metabolism, endocrinology, microbiology, and immunology to the etiology and maintenance of AN.

Distribution of the voltage-gated sodium channel Na(v)1.7 in the rat: expression in the autonomic and endocrine systems.

It is generally accepted that the voltage-gated, tetrodotoxin-sensitive sodium channel, Na(V)1.7, is selectively expressed in peripheral ganglia. However, global deletion in mice of Na(V)1.7 leads to death shortly after birth (Nassar et al. [2004] Proc. Natl. Acad. Sci. U. S. A. 101:12706-12711), suggesting that this ion channel might be more widely expressed. To understand better the potential physiological function of this ion channel, we examined Na(V)1.7 expression in the rat by in situ hybridization and immunohistochemistry. As expected, highest mRNA expression levels are found in peripheral ganglia, and the protein is expressed within these ganglion cells and on the projections of these neurons in the central nervous system. Importantly, we found that Na(V)1.7 is present in discrete rat brain regions, and the unique distribution pattern implies a central involvement in endocrine and autonomic systems as well as analgesia. In addition, Na(V)1.7 expression was detected in the pituitary and adrenal glands. These results indicate that Na(V)1.7 is not only involved in the processing of sensory information but also participates in the regulation of autonomic and endocrine systems; more specifically, it could be implicated in such vital functions as fluid homeostasis and cardiovascular control.

Patterned gene programs and target remodeling following axotomy at a major site for sensory innervation.

The genetic mechanisms that a target uses to reestablish the connections of regenerating axons were explored using oligonucleotide microarrays and real-time PCR. In normal and denervated mouse vibrissa follicles, patterns of genetic regulation were assessed in adjacent targets that normally receive different types of sensory and autonomic innervation. We show that a target remodeling occurs following axotomy involving reduced hair growth, altered hair follicle integrity and remodeling of the extracellular matrix. Also, we found two orphan receptors putatively involved in hair growth. Our data further demonstrate region-specific regulation of genes putatively involved in target-axon interactions. Thus, this study shows for the first time that major target remodeling occurs following denervation and suggests putative functions for several novel genes.

Diverse dependencies of developing Merkel innervation on the trkA and both full-length and truncated isoforms of trkC.

This study demonstrates that innervation dependent on two different neurotrophin tyrosine kinase (trk) receptors can form the same types of sensory endings (Merkel endings) in the same target (Merkel cells of vibrissa follicles). Some endings transiently express trkA during their initial development, whereas others express trkC throughout their development. Consequently, elimination of kinase domains of either trkA or trkC each result in a partial loss of Merkel endings, whereas absence of kinase domains of both receptors results in a total loss. At the onset of Merkel ending development, at least one kinase-lacking trkC isoform is transiently expressed on all the follicle cells, while neurotrophin 3 is transiently expressed only in the cells at the middle third of the follicle where the Merkel endings and cells develop. This transient non-neuronal expression of truncated trkC is essential for development of any Merkel endings, whereas some Merkel endings and cells still begin to develop in the absence of neurotrophin 3. Therefore, truncated trkC plays a more important role in the development of this innervation than kinase forms of trkA or trkC or of NT3, the only known ligand for trkC receptors.