TMEM132D, a new candidate for anxiety phenotypes: evidence from human and mouse studies

Genome-wide association study of major anxiety disorders in 122,341 European-ancestry cases identifies 58 loci and highlights GABAergic signaling

The major anxiety disorders (ANX; including generalized anxiety disorder, panic disorder, and phobias) are highly prevalent, often onset early, persist throughout life, and cause substantial global disability. Although distinct in their clinical presentations, they likely represent differential expressions of a dysregulated threat-response system. Here we present a genome-wide association meta-analysis comprising 122,341 European ancestry ANX cases and 729,881 controls. We identified 58 independent genome-wide significant ANX risk variants and 66 genes with robust biological support. In an independent sample of 1,175,012 self-report ANX cases and 1,956,379 controls, 51 of the 58 associated variants were replicated. As predicted by twin studies, we found substantial genetic correlation between ANX and depression, neuroticism, and other internalizing phenotypes. Follow-up analyses demonstrated enrichment in all major brain regions and highlighted GABAergic signaling as one potential mechanism underlying ANX genetic risk. These results advance our understanding of the genetic architecture of ANX and prioritize genes for functional follow-up studies.

Immunoglobulin genes expressed in lymphoblastoid cell lines discern and predict lithium response in bipolar disorder patients

Bipolar disorder (BD) is a neuropsychiatric mood disorder manifested by recurrent episodes of mania and depression. More than half of BD patients are non-responsive to lithium, the first-line treatment drug, complicating BD clinical management. Given its unknown etiology, it is pertinent to understand the genetic signatures that lead to variability in lithium response. We discovered a set of differentially expressed genes (DEGs) from the lymphoblastoid cell lines (LCLs) of 10 controls and 19 BD patients belonging mainly to the immunoglobulin gene family that can be used as potential biomarkers to diagnose and treat BD. Importantly, we trained machine learning algorithms on our datasets that predicted the lithium response of BD subtypes with minimal errors, even when used on a different cohort of 24 BD patients acquired by a different laboratory. This proves the scalability of our methodology for predicting lithium response in BD and for a prompt and suitable decision on therapeutic interventions.

Genetic insights into the neurobiology of anxiety

Anxiety and fear are evolutionarily conserved emotions that increase the likelihood of an organism surviving threatening situations. Anxiety and vigilance states are regulated by neural networks involving multiple brain regions. In anxiety disorders, this intricate regulatory system is disturbed, leading to excessive or prolonged anxiety or fear. Anxiety disorders have both genetic and environmental risk factors. Genetic research has the potential to identify specific genetic variants causally associated with specific phenotypes. In recent decades, genome-wide association studies (GWASs) have revealed variants predisposing to neuropsychiatric disorders, suggesting novel neurobiological pathways in the etiology of these disorders. Here, we review recent human GWASs of anxiety disorders, and genetic studies of anxiety-like behavior in rodent models. These studies are paving the way for a better understanding of the neurobiological mechanisms underlying anxiety disorders.

Using coding and non-coding rare variants to target candidate genes in patients with severe tinnitus

Tinnitus is the phantom percept of an internal non-verbal set of noises and tones. It is reported by 15% of the population and it is usually associated with hearing and/or brain disorders. The role of structural variants (SVs) in coding and non-coding regions has not been investigated in patients with severe tinnitus. In this study, we performed whole-genome sequencing in 97 unrelated Swedish individuals with chronic tinnitus (TIGER cohort). Rare single nucleotide variants (SNV), large structural variants (LSV), and copy number variations (CNV) were retrieved to perform a gene enrichment analysis in TIGER and in a subgroup of patients with severe tinnitus (SEVTIN, n = 34), according to the tinnitus handicap inventory (THI) scores. An independent exome sequencing dataset of 147 Swedish tinnitus patients was used as a replication cohort (JAGUAR cohort) and population-specific datasets from Sweden (SweGen) and Non-Finish Europeans (NFE) from gnomAD were used as control groups. SEVTIN patients showed a higher prevalence of hyperacusis, hearing loss, and anxiety when they were compared to individuals in the TIGER cohort. We found an enrichment of rare missense variants in 6 and 8 high-constraint genes in SEVTIN and TIGER cohorts, respectively. Of note, an enrichment of missense variants was found in the CACNA1E gene in both SEVTIN and TIGER. We replicated the burden of missense variants in 9 high-constrained genes in the JAGUAR cohort, including the gene NAV2, when data were compared with NFE. Moreover, LSVs in constrained regions overlapping CACNA1E, NAV2, and TMEM132D genes were observed in TIGER and SEVTIN.

Systematic review of genome-wide association studies of anxiety disorders and neuroticism

OBJECTIVES

To summarise SNP associations identified by genome-wide association studies (GWASs) of anxiety disorders and neuroticism; to appraise the quality of individual studies, and to assess the ancestral diversity of study participants.

METHODS

We searched PubMed, Scopus, PsychInfo and PubPsych for GWASs of anxiety disorders, non-diagnostic traits (such as anxiety sensitivity), and neuroticism, and extracted all SNPs that surpassed genome-wide significance. We graded study quality using Q-genie scores and reviewed the ancestral diversity of included participants.

RESULTS

32 studies met our inclusion criteria. A total of 563 independent significant variants were identified, of which 29 were replicated nominally in independent samples, and 3 were replicated significantly. The studies had good global quality, but many smaller studies were underpowered. Phenotypic heterogeneity for anxiety (and less so for neuroticism) seemed to reflect the complexity of capturing this trait. Ancestral diversity was poor, with 70% of studies including only populations of European ancestry.

CONCLUSION

The functionality of genes identified by GWASs of anxiety and neuroticism deserves further investigation. Future GWASs should have larger sample sizes, more rigorous phenotyping and include more ancestrally diverse population groups.

Differential expression of the Tmem132 family genes in the developing mouse nervous system

The family of novel transmembrane proteins (TMEM) 132 have been associated with multiple neurological disorders and cancers in humans, but have hardly been studied in vivo. Here we report the expression patterns of the five Tmem132 genes (a, b, c, d and e) in developing mouse nervous system with RNA in situ hybridization in wholemount embryos and tissue sections. Our results reveal differential and partially overlapping expression of multiple Tmem132 family members in both the central and peripheral nervous system, suggesting potential partial redundancy among them.

Effects of stressful life-events on DNA methylation in panic disorder and major depressive disorder

BACKGROUND

Panic disorder (PD) is characterized by recurrent panic attacks and higher affection of women as compared to men. The lifetime prevalence of PD is about 2-3% in the general population leading to tremendous distress and disability. Etiologically, genetic and environmental factors, such as stress, contribute to the onset and relapse of PD. In the present study, we investigated epigenome-wide DNA methylation (DNAm) in respond to a cumulative, stress-weighted life events score (wLE) in patients with PD and its boundary to major depressive disorder (MDD), frequently co-occurring with symptoms of PD.

METHODS

DNAm was assessed by the Illumina HumanMethylation450 BeadChip. In a meta-analytic approach, epigenome-wide DNAm changes in association with wLE were first analyzed in two PD cohorts (with a total sample size of 183 PD patients and 85 healthy controls) and lastly in 102 patients with MDD to identify possible overlapping and opposing effects of wLE on DNAm. Additionally, analysis of differentially methylated regions (DMRs) was conducted to identify regional clusters of association.

RESULTS

Two CpG-sites presented with p-values below 1 × 10-05 in PD: cg09738429 (p = 6.40 × 10-06, located in an intergenic shore region in next proximity of PYROXD1) and cg03341655 (p = 8.14 × 10-06, located in the exonic region of GFOD2). The association of DNAm at cg03341655 and wLE could be replicated in the independent MDD case sample indicating a diagnosis independent effect. Genes mapping to the top hits were significantly upregulated in brain and top hits have been implicated in the metabolic system. Additionally, two significant DMRs were identified for PD only on chromosome 10 and 18, including CpG-sites which have been reported to be associated with anxiety and other psychiatric phenotypes.

CONCLUSION

This first DNAm analysis in PD reveals first evidence of small but significant DNAm changes in PD in association with cumulative stress-weighted life events. Most of the top associated CpG-sites are located in genes implicated in metabolic processes supporting the hypothesis that environmental stress contributes to health damaging changes by affecting a broad spectrum of systems in the body.

Transcriptome-Wide Association Study Provides Insights Into the Genetic Component of Gene Expression in Anxiety

Anxiety disorders are common mental disorders that often result in disability. Recently, large-scale genome-wide association studies (GWASs) have identified several novel risk variants and loci for anxiety disorders (or anxiety traits). Nevertheless, how the reported risk variants confer risk of anxiety remains unknown. To identify genes whose cis-regulated expression levels are associated with risk of anxiety traits, we conducted a transcriptome-wide association study (TWAS) by integrating genome-wide associations from a large-scale GWAS (N = 175,163) (which evaluated anxiety traits based on Generalized Anxiety Disorder 2-item scale (GAD-2) score) and brain expression quantitative trait loci (eQTL) data (from the PsychENCODE and GTEx). We identified 19 and 17 transcriptome-wide significant (TWS) genes in the PsychENCODE and GTEx, respectively. Intriguingly, 10 genes showed significant associations with anxiety in both datasets, strongly suggesting that genetic risk variants may confer risk of anxiety traits by regulating the expression of these genes. Top TWS genes included RNF123, KANSL1-AS1, GLYCTK, CRHR1, DND1P1, MAPT and ARHGAP27. Of note, 25 TWS genes were not implicated in the original GWAS. Our TWAS identified 26 risk genes whose cis-regulated expression were significantly associated with anxiety, providing important insights into the genetic component of gene expression in anxiety disorders/traits and new clues for future drug development.

Genetic contributions to anxiety disorders: where we are and where we are heading

Anxiety disorders are among the most common psychiatric disorders worldwide. They often onset early in life, with symptoms and consequences that can persist for decades. This makes anxiety disorders some of the most debilitating and costly disorders of our time. Although much is known about the synaptic and circuit mechanisms of fear and anxiety, research on the underlying genetics has lagged behind that of other psychiatric disorders. However, alongside the formation of the Psychiatric Genomic Consortium Anxiety workgroup, progress is rapidly advancing, offering opportunities for future research.Here we review current knowledge about the genetics of anxiety across the lifespan from genetically informative designs (i.e. twin studies and molecular genetics). We include studies of specific anxiety disorders (e.g. panic disorder, generalised anxiety disorder) as well as those using dimensional measures of trait anxiety. We particularly address findings from large-scale genome-wide association studies and show how such discoveries may provide opportunities for translation into improved or new therapeutics for affected individuals. Finally, we describe how discoveries in anxiety genetics open the door to numerous new research possibilities, such as the investigation of specific gene-environment interactions and the disentangling of causal associations with related traits and disorders.We discuss how the field of anxiety genetics is expected to move forward. In addition to the obvious need for larger sample sizes in genome-wide studies, we highlight the need for studies among young people, focusing on specific underlying dimensional traits or components of anxiety.

A review of systems biology research of anxiety disorders

The development of "omic" technologies and deep phenotyping may facilitate a systems biology approach to understanding anxiety disorders. Systems biology approaches incorporate data from multiple modalities (e.g., genomic, neuroimaging) with functional analyses (e.g., animal and tissue culture models) and mathematical modeling (e.g., machine learning) to investigate pathological biophysical networks at various scales. Here we review: i) the neurobiology of anxiety disorders; ii) how systems biology approaches have advanced this work; and iii) the clinical implications and future directions of this research. Systems biology approaches have provided an improved functional understanding of candidate biomarkers and have suggested future potential for refining the diagnosis, prognosis, and treatment of anxiety disorders. The systems biology approach for anxiety disorders is, however, in its infancy and in some instances is characterized by insufficient power and replication. The studies reviewed here represent important steps to further untangling the pathophysiology of anxiety disorders.

Genome-wide association study of panic disorder reveals genetic overlap with neuroticism and depression

Panic disorder (PD) has a lifetime prevalence of 2-4% and heritability estimates of 40%. The contributory genetic variants remain largely unknown, with few and inconsistent loci having been reported. The present report describes the largest genome-wide association study (GWAS) of PD to date comprising genome-wide genotype data of 2248 clinically well-characterized PD patients and 7992 ethnically matched controls. The samples originated from four European countries (Denmark, Estonia, Germany, and Sweden). Standard GWAS quality control procedures were conducted on each individual dataset, and imputation was performed using the 1000 Genomes Project reference panel. A meta-analysis was then performed using the Ricopili pipeline. No genome-wide significant locus was identified. Leave-one-out analyses generated highly significant polygenic risk scores (PRS) (explained variance of up to 2.6%). Linkage disequilibrium (LD) score regression analysis of the GWAS data showed that the estimated heritability for PD was 28.0-34.2%. After correction for multiple testing, a significant genetic correlation was found between PD and major depressive disorder, depressive symptoms, and neuroticism. A total of 255 single-nucleotide polymorphisms (SNPs) with p < 1 × 10^-4 were followed up in an independent sample of 2408 PD patients and 228,470 controls from Denmark, Iceland and the Netherlands. In the combined analysis, SNP rs144783209 showed the strongest association with PD (pcomb = 3.10  × 10^-7). Sign tests revealed a significant enrichment of SNPs with a discovery p-value of <0.0001 in the combined follow up cohort (p = 0.048). The present integrative analysis represents a major step towards the elucidation of the genetic susceptibility to PD.

The C. elegans homolog of human panic-disorder risk gene TMEM132D orchestrates neuronal morphogenesis through the WAVE-regulatory complex

TMEM132D is a human gene identified with multiple risk alleles for panic disorders, anxiety and major depressive disorders. Defining a conserved family of transmembrane proteins, TMEM132D and its homologs are still of unknown molecular functions. By generating loss-of-function mutants of the sole TMEM132 ortholog in C. elegans, we identify abnormal morphologic phenotypes in the dopaminergic PDE neurons. Using a yeast two-hybrid screen, we find that NAP1 directly interacts with the cytoplasmic domain of human TMEM132D, and mutations in C. elegans tmem-132 that disrupt interaction with NAP1 cause similar morphologic defects in the PDE neurons. NAP1 is a component of the WAVE regulatory complex (WRC) that controls F-actin cytoskeletal dynamics. Decreasing activity of WRC rescues the PDE defects in tmem-132 mutants, whereas gain-of-function of TMEM132D in mammalian cells inhibits WRC, leading to decreased abundance of select WRC components, impaired actin nucleation and cell motility. We propose that metazoan TMEM132 family proteins play evolutionarily conserved roles in regulating NAP1 protein homologs to restrict inappropriate WRC activity, cytoskeletal and morphologic changes in the cell.

Efficient multiplexed genome engineering with a polycistronic tRNA and CRISPR guide-RNA reveals an important role of detonator in reproduction of Drosophila melanogaster

Genome association studies in human and genetic studies in mouse implicated members of the transmembrane protein 132 (TMEM132) family in multiple conditions including panic disorder, hearing loss, limb and kidney malformation. However, the presence of five TMEM132 paralogs in mammalian genomes makes it extremely challenging to reveal the full requirement for these proteins in vivo. In contrast, there is only one TMEM132 homolog, detonator (dtn), in the genome of fruit fly Drosophila melanogaster, enabling straightforward research into its in vivo function. In the current study, we generate multiple loss-of-function dtn mutant fly strains through a polycistronic tRNA-gRNA approach, and show that most embryos lacking both maternal and paternal dtn fail to hatch into larvae, indicating an essential role of dtn in Drosophila reproduction.

Shared transethnic genetic basis of panic disorder and psychiatric and related intermediate phenotypes

Panic disorder (PD), a common anxiety disorder, is modestly heritable. The genetic basis of anxiety disorders overlaps with that of other psychiatric disorders and their intermediate phenotypes in individuals of European ancestry. Here, we investigated the transethnic polygenetic features shared between Japanese PD patients and European patients with psychiatric disorders and their intermediate phenotypes by conducting polygenic risk score (PRS) analyses. Large-scale European genome-wide association study (GWAS) datasets (n = 7,556-1,131,881) for ten psychiatric disorders and seven intermediate phenotypes were utilized as discovery samples. PRSs derived from these GWASs were calculated for Japanese target subjects [718 PD patients and 1,717 healthy controls (HCs)]. The effects of these PRSs from European GWASs on the risk of PD in Japanese patients were investigated. The PRSs from European studies of anxiety disorders were marginally higher in Japanese PD patients than in HCs (p = 0.013). Regarding other psychiatric disorders, the PRSs for depression in European patients were significantly higher in Japanese PD patients than in HCs (p = 2.31×10^-4), while the PRSs for attention-deficit/hyperactivity disorder in European patients were nominally lower in Japanese PD patients than in HCs (p = 0.024). Regarding health-related, personality-based and cognitive intermediate phenotypes, the PRSs for loneliness (especially isolation) in European individuals were significantly higher in Japanese PD patients than in HCs (p = 9.02×10^-4). Furthermore, Japanese PD patients scored nominally higher than HCs in PRSs for neuroticism in European people (p = 3.37×10^-3), while Japanese PD patients scored nominally lower than HCs in PRSs for tiredness (p = 0.025), educational attainment (p = 0.035) and cognitive function (p = 9.63×10^-3). Our findings suggest that PD shares transethnic genetic etiologies with other psychiatric disorders and related intermediate phenotypes.

Genetic Biomarkers of Panic Disorder: A Systematic Review

(1) Background: Although panic disorder (PD) is one of the most common anxiety disorders severely impacting quality of life, no effective genetic testing exists; known data on possible genetic biomarkers is often scattered and unsystematic which complicates further studies. (2) Methods: We used PathwayStudio 12.3 (Elsevier, The Netherlands) to acquire literature data for further manual review and analysis. 229 articles were extracted, 55 articles reporting associations, and 32 articles reporting no associations were finally selected. (3) Results: We provide exhaustive information on genetic biomarkers associated with PD known in the scientific literature. Data is presented in two tables. Genes COMT and SLC6A4 may be considered the most promising for PD diagnostic to date. (4) Conclusions: This review illustrates current progress in association studies of PD and may indicate possible molecular mechanisms of its pathogenesis. This is a possible basis for data analysis, novel experimental studies, or developing test systems and personalized treatment approaches.

Compositional Features and Codon Usage Pattern of Genes Associated with Anxiety in Human

Codon usage bias (CUB) is the unequal usage of synonymous codon; some codons are more preferred than others. CUB analysis has applications in understanding the molecular organization of genome, genetics, gene expression, and molecular evolution. Bioinformatic approach was used to analyze the protein-coding sequences of genes involved in the anxiety to understand the patterns of codon usage as no work was reported yet. The improved effective number of codons (Nc) values ranged from 43.55 to 55.06, with a mean of 44.57, suggested that the overall CUB was low for genes associated with anxiety. The overall GC and AT content was 54.76 and 45.24, respectively. Relative synonymous codon usage (RSCU) analysis revealed that most frequently used codon ended mostly with C or G. The over-represented codons in genes associated with anxiety were CTG, ATC, GTG, AGC, ACC, and GCC, while under-represented codons were TTA, CTT, CTA, ATA, GTT, GTA, TCG, CCG, GCG, CAA, and CGT. Correlation analysis was performed between overall nucleotide composition and its 3rd codon positions, and observed highly significant (p < 0.01) correlation between them suggested that both mutation pressure and natural selection might affect the pattern of CUB. The highly significant correlation (0.598**, p < 0.01) was also observed between GC12 with GC3 suggested that directional mutation pressure might acted on all codon positions for genes associated with anxiety.

Genome-wide association study of shared liability to anxiety disorders in Army STARRS

Anxiety disorders (ANX), namely generalized anxiety, panic disorder, and phobias, are common, etiologically complex syndromes that show increasing prevalence and comorbidity throughout adolescence and beyond. Few genome-wide association studies (GWAS) examining ANX risk have been published and almost exclusively in individuals of European ancestry. In this study, we phenotyped participants from the Army Study To Assess Risk and Resilience in Servicemembers (STARRS) to approximate DSM-based ANX diagnoses. We factor-analyzed those to create a single dimensional anxiety score for each subject. GWAS were conducted using that score within each of three ancestral groups (EUR, AFR, LAT) and then meta-analyzed across ancestries (NTotal = 16,510). We sought to (a) replicate prior ANX GWAS findings in ANGST; (b) determine whether results extended to other ancestry groups; and (c) meta-analyze with ANGST for increased power to identify novel susceptibility loci. No reliable genome-wide significant SNP associations were detected in STARRS. However, SNPs within the CAMKMT gene located in region 2p21 associated with shared ANX risk in ANGST were replicated in EUR soldiers but not other ancestry groups. Combining EUR STARRS and ANGST (N = 28,950) yielded a more robust 2p21 association signal (p = 9.08x10-11 ). Gene-based analyses supported three genes within 2p21 and LBX1 on chromosome 10. More powerful ANX genetic studies will be required to identify further loci.

Circular RNA expression and association with the clinicopathological characteristics in papillary thyroid carcinoma

Papillary thyroid carcinoma (PTC) is the most common type of thyroid cancer. Circular RNAs (circRNAs) are a novel class of RNAs, with higher stability and tissue specificity, which may be of value as novel clinical markers. High-throughput RNA sequencing was used to profile the expression of circRNAs in 5 pairs of cancer and normal tissues, and reverse transcription-quantitative PCR (RT-qPCR) analysis was employed to verify the results of the RNA sequencing in 45 cases of PTC. The dysregulated circRNA expression and clinicopathological characteristics were assessed and the potential roles of circRNAs in the cellular miRNA and mRNA network were predicted using bioinformatics analysis. The results demonstrated that, compared with normal tissues, a total of 53 circRNAs were dysregulated in tumour tissues, and 8 circRNAs were validated at the mRNA level (P<0.001 and P<0.01). Among those, the expression of chr5:161330882-161336769- (P=0.015), chr9:22046750-22097364+ (P=0.041) and chr8:18765448-18804898- (P=0.036) were obviously associated with the BRAF^V600E mutation, chr12:129699809-129700698- was associated with capsular invasion (P=0.025) and chr5:38523418-38530666- was associated with pT stage (P=0.037) and lymph node metastasis (P=0.002). Therefore, some dysregulated circRNAs were found to be associated with BRAF^V600E mutation, capsular invasion, advanced pT stage and lymph node metastasis of PTC, indicating that circRNAs may be involved in tumourigenesis and cancer progression, and they may be putative biomarkers for the diagnosis and evaluation of progression of PTC.

Genetic and epigenetic analyses of panic disorder in the post-GWAS era

Panic disorder (PD) is a common and debilitating neuropsychiatric disorder characterized by panic attacks coupled with excessive anxiety. Both genetic factors and environmental factors play an important role in PD pathogenesis and response to treatment. However, PD is clinically heterogeneous and genetically complex, and the exact genetic or environmental causes of this disorder remain unclear. Various approaches for detecting disease-causing genes have recently been made available. In particular, genome-wide association studies (GWAS) have attracted attention for the identification of disease-associated loci of multifactorial disorders. This review introduces GWAS of PD, followed by a discussion about the limitations of GWAS and the major challenges facing geneticists in the post-GWAS era. Alternative strategies to address these challenges are then proposed, such as epigenome-wide association studies (EWAS) and rare variant association studies (RVAS) using next-generation sequencing. To date, however, few reports have described these analyses, and the evidence remains insufficient to confidently identify or exclude rare variants or epigenetic changes in PD. Further analyses are therefore required, using sample sizes in the tens of thousands, extensive functional annotations, and highly targeted hypothesis testing.

Shared genetic etiology between anxiety disorders and psychiatric and related intermediate phenotypes

BACKGROUND

Psychiatric disorders and related intermediate phenotypes are highly heritable and have a complex, overlapping polygenic architecture. A large-scale genome-wide association study (GWAS) of anxiety disorders identified genetic variants that are significant on a genome-wide. The current study investigated the genetic etiological overlaps between anxiety disorders and frequently cooccurring psychiatric disorders and intermediate phenotypes.

METHODS

Using case-control and factor score models, we investigated the genetic correlations of anxiety disorders with eight psychiatric disorders and intermediate phenotypes [the volumes of seven subcortical brain regions, childhood cognition, general cognitive ability and personality traits (subjective well-being, loneliness, neuroticism and extraversion)] from large-scale GWASs (n = 7556-298 420) by linkage disequilibrium score regression.

RESULTS

Among psychiatric disorders, the risk of anxiety disorders was positively genetically correlated with the risks of major depressive disorder (MDD) (rg ± standard error = 0.83 ± 0.16, p = 1.97 × 10-7), schizophrenia (SCZ) (0.28 ± 0.09, p = 1.10 × 10-3) and attention-deficit/hyperactivity disorder (ADHD) (0.34 ± 0.13, p = 8.40 × 10-3). Among intermediate phenotypes, significant genetic correlations existed between the risk of anxiety disorders and neuroticism (0.81 ± 0.17, p = 1.30 × 10-6), subjective well-being (-0.73 ± 0.18, p = 4.89 × 10-5), general cognitive ability (-0.23 ± 0.08, p = 4.70 × 10-3) and putamen volume (-0.50 ± 0.18, p = 5.00 × 10-3). No other significant genetic correlations between anxiety disorders and psychiatric or intermediate phenotypes were observed (p > 0.05). The case-control model yielded stronger genetic effect sizes than the factor score model.

CONCLUSIONS

Our findings suggest that common genetic variants underlying the risk of anxiety disorders contribute to elevated risks of MDD, SCZ, ADHD and neuroticism and reduced quality of life, putamen volume and cognitive performance. We suggest that the comorbidity of anxiety disorders is partly explained by common genetic variants.

Oscillotherapeutics - Time-targeted interventions in epilepsy and beyond

Oscillatory brain activities support many physiological functions from motor control to cognition. Disruptions of the normal oscillatory brain activities are commonly observed in neurological and psychiatric disorders including epilepsy, Parkinson's disease, Alzheimer's disease, schizophrenia, anxiety/trauma-related disorders, major depressive disorders, and drug addiction. Therefore, these disorders can be considered as common oscillation defects despite having distinct behavioral manifestations and genetic causes. Recent technical advances of neuronal activity recording and analysis have allowed us to study the pathological oscillations of each disorder as a possible biomarker of symptoms. Furthermore, recent advances in brain stimulation technologies enable time- and space-targeted interventions of the pathological oscillations of both neurological disorders and psychiatric disorders as possible targets for regulating their symptoms.

Genetic Variants Associated With Anxiety and Stress-Related Disorders: A Genome-Wide Association Study and Mouse-Model Study

IMPORTANCE

Anxiety and stress-related disorders are among the most common mental disorders. Although family and twin studies indicate that both genetic and environmental factors play an important role underlying their etiology, the genetic underpinnings of anxiety and stress-related disorders are poorly understood.

OBJECTIVES

To estimate the single-nucleotide polymorphism-based heritability of anxiety and stress-related disorders; to identify novel genetic risk variants, genes, or biological pathways; to test for pleiotropic associations with other psychiatric traits; and to evaluate the association of psychiatric comorbidities with genetic findings.

DESIGN, SETTING, PARTICIPANTS

This genome-wide association study included individuals with various anxiety and stress-related diagnoses and controls derived from the population-based Lundbeck Foundation Initiative for Integrative Psychiatric Research (iPSYCH) study. Lifetime diagnoses of anxiety and stress-related disorders were obtained through the national Danish registers. Genes of interest were further evaluated in mice exposed to chronic social defeat. The study was conducted between June 2016 and November 2018.

MAIN OUTCOMES AND MEASURES

Diagnoses of a relatively broad diagnostic spectrum of anxiety and stress-related disorders.

RESULTS

The study sample included 12 655 individuals with various anxiety and stress-related diagnoses and 19 225 controls. Overall, 17 740 study participants (55.6%) were women. A total of 7308 participants (22.9%) were born between 1981-1985, 8840 (27.7%) between 1986-1990, 8157 (25.6%) between 1991-1995, 5918 (18.6%) between 1996-2000, and 1657 (5.2%) between 2001-2005. Standard association analysis revealed variants in PDE4B to be associated with anxiety and stress-related disorder (rs7528604; P = 5.39 × 10-11; odds ratio = 0.89; 95% CI, 0.86-0.92). A framework of sensitivity analyses adjusting for mental comorbidity supported this result showing consistent association of PDE4B variants with anxiety and stress-related disorder across analytical scenarios. In mouse models, alterations in Pde4b expression were observed in those mice displaying anxiety-like behavior after exposure to chronic stress in the prefrontal cortex (P = .002; t = -3.33) and the hippocampus (P = .001; t = -3.72). We also found a single-nucleotide polymorphism heritability of 28% (standard error = 0.027) and that the genetic signature of anxiety and stress-related overlapped with psychiatric traits, educational outcomes, obesity-related phenotypes, smoking, and reproductive success.

CONCLUSIONS AND RELEVANCE

This study highlights anxiety and stress-related disorders as complex heritable phenotypes with intriguing genetic correlations not only with psychiatric traits, but also with educational outcomes and multiple obesity-related phenotypes. Furthermore, we highlight the candidate gene PDE4B as a robust risk locus pointing to the potential of PDE4B inhibitors in treatment of these disorders.

A genome-wide association meta-analysis of prognostic outcomes following cognitive behavioural therapy in individuals with anxiety and depressive disorders

Major depressive disorder and the anxiety disorders are highly prevalent, disabling and moderately heritable. Depression and anxiety are also highly comorbid and have a strong genetic correlation (rg ≈ 1). Cognitive behavioural therapy is a leading evidence-based treatment but has variable outcomes. Currently, there are no strong predictors of outcome. Therapygenetics research aims to identify genetic predictors of prognosis following therapy. We performed genome-wide association meta-analyses of symptoms following cognitive behavioural therapy in adults with anxiety disorders (n = 972), adults with major depressive disorder (n = 832) and children with anxiety disorders (n = 920; meta-analysis n = 2724). We estimated the variance in therapy outcomes that could be explained by common genetic variants (h2SNP) and polygenic scoring was used to examine genetic associations between therapy outcomes and psychopathology, personality and learning. No single nucleotide polymorphisms were strongly associated with treatment outcomes. No significant estimate of h2SNP could be obtained, suggesting the heritability of therapy outcome is smaller than our analysis was powered to detect. Polygenic scoring failed to detect genetic overlap between therapy outcome and psychopathology, personality or learning. This study is the largest therapygenetics study to date. Results are consistent with previous, similarly powered genome-wide association studies of complex traits.

Genetics of Anxiety Disorders

PURPOSE OF REVIEW

Anxiety disorders are among the most common mental disorders with a lifetime prevalence of over 20%. Clinically, anxiety is not thought of as a homogenous disorder, but is subclassified in generalized, panic, and phobic anxiety disorder. Anxiety disorders are moderately heritable. This review will explore recent genetic and epigenetic approaches to anxiety disorders explaining differential susceptibility risk.

RECENT FINDINGS

A substantial portion of the variance in susceptibility risk can be explained by differential inherited and acquired genetic and epigenetic risk. Available data suggest that anxiety disorders are highly complex and polygenic. Despite the substantial progress in genetic research over the last decade, only few risk loci for anxiety disorders have been identified so far. This review will cover recent findings from large-scale genome-wide association studies as well as newer epigenome-wide studies. Progress in this area will likely require analysis of much larger sample sizes than have been reported to date. We discuss prospects for clinical translation of genetic findings and future directions for research.

Concordance of genetic variation that increases risk for anxiety disorders and posttraumatic stress disorders and that influences their underlying neurocircuitry

BACKGROUND

There have been considerable recent advances in understanding the genetic architecture of anxiety disorders and posttraumatic stress disorder (PTSD), as well as the underlying neurocircuitry of these disorders. However, there is little work on the concordance of genetic variations that increase risk for these conditions, and that influence subcortical brain structures. We undertook a genome-wide investigation of the overlap between the genetic influences from single nucleotide polymorphisms (SNPs) on volumes of subcortical brain structures and genetic risk for anxiety disorders and PTSD.

METHOD

We obtained summary statistics of genome-wide association studies (GWAS) of anxiety disorders (N_cases= 7016, N_controls= 14,745), PTSD (European sample; N_cases= 2424, N_controls= 7113) and of subcortical brain structures (N = 13,171). SNP Effect Concordance Analysis (SECA) and Linkage Disequilibrium (LD) Score Regression were used to examine genetic pleiotropy, concordance, and genome-wide correlations respectively. SECAs conditional false discovery was used to identify specific risk variants associated with anxiety disorders or PTSD when conditioning on brain related traits.

RESULTS

For anxiety disorders, we found evidence of significant concordance between increased anxiety risk variants and variants associated with smaller amygdala volume. Further, by conditioning on brain volume GWAS, we identified novel variants that associate with smaller brain volumes and increase risk for disorders: rs56242606 was found to increase risk for anxiety disorders, while two variants (rs6470292 and rs683250) increase risk for PTSD, when conditioning on the GWAS of putamen volume.

LIMITATIONS

Despite using the largest available GWAS summary statistics, the analyses were limited by sample size.

CONCLUSIONS

These preliminary data indicate that there is genome wide concordance between genetic risk factors for anxiety disorders and those for smaller amygdala volume, which is consistent with research that supports the involvement of the amygdala in anxiety disorders. It is notable that a genetic variant that contributes to both reduced putamen volume and PTSD plays a key role in the glutamatergic system. Further work with GWAS summary statistics from larger samples, and a more extensive look at the genetics underlying brain circuits, is needed to fully delineate the genetic architecture of these disorders and their underlying neurocircuitry.

Epigenetic Mechanisms Within the Cingulate Cortex Regulate Innate Anxiety-Like Behavior

BACKGROUND

Pathological anxiety originates from a complex interplay of genetic predisposition and environmental factors, acting via epigenetic mechanisms. Epigenetic processes that can counteract detrimental genetic risk towards innate high anxiety are not well characterized.

METHODS

We used female mouse lines of selectively bred high (HAB)- vs low (LAB)-innate anxiety-related behavior and performed select environmental and pharmacological manipulations to alter anxiety levels as well as brain-specific manipulations and immunohistochemistry to investigate neuronal mechanisms associated with alterations in anxiety-related behavior.

RESULTS

Inborn hyperanxiety of high anxiety-like phenotypes was effectively reduced by environmental enrichment exposure. c-Fos mapping revealed that hyperanxiety in high anxiety-like phenotypes was associated with blunted challenge-induced neuronal activation in the cingulate-cortex, which was normalized by environmental enrichment. Relating this finding with epigenetic modifications, we found that high anxiety-like phenotypes (compared with low-innate anxiety phenotypes) showed reduced acetylation in the hypoactivated cingulate-cortex neurons following a mild emotional challenge, which again was normalized by environmental enrichment. Paralleling the findings using environmental enrichment, systemic administration of histone-deacetylase-inhibitor MS-275 elicited an anxiolytic-like effect, which was correlated with increased acetylated-histone-3 levels within cingulate-cortex. Finally, as a proof-of-principle, local MS-275 injection into cingulate-cortex rescued enhanced innate anxiety and increased acetylated-histone-3 within the cingulate-cortex, suggesting this epigenetic mark as a biomarker for treatment success.

CONCLUSIONS

Taken together, the present findings provide the first causal evidence that the attenuation of high innate anxiety-like behavior via environmental/pharmacological manipulations is epigenetically mediated via acetylation changes within the cingulate-cortex. Finally, histone-3 specific histone-deacetylase-inhibitor could be of therapeutic importance in anxiety disorders.

Panic disorders: The role of genetics and epigenetics

Panic disorder is characterized by symptoms with abrupt surges of fear with palpitations, sweating, trembling, heat sensations. Considering its disease burden on each individual and on society, understanding its etiology is important. Though no one specific etiology has been known, like other psychiatric disorders, multiple factors such as genetic, environmental, neurobiological, psychopathological factors have been suggested. In this article, we reviewed currently known etiologies and related study results, regarding especially genetic and epigenetic aspects of the panic disorder. Early studies, including twin studies, family studies, adoption studies suggested highly familial trait of panic disorder. Linkage studies, either, found panic disorder is not a single gene disorder but confirmed existence of multiple related genes. Chromosome and candidate gene studies found few related genes, NPY, ADORA2A, COMT, IKBKE. Newer method, genome-wide association studies (GWAS) have been searching for newer genes. No genome-wide significant genes, however, were detected, confirming previously known candidate genes, NPY5R on 4q31.3-32, BDKRB2 on 14q32, instead. Epigenetic modification has also been studied on many different psychiatric disorders. Monoamine oxidase A (MAOA) hypomethylation, taken together with negative life events, showed relation with panic disorder. Glutamate decarbodylases 1 (GAD1) hypomethylation was also specific on panic disorder patients. Relation with noradrenaline transporter (NET) gene SLC6a2 promoter methylation has also been studied. In conclusion, no specific gene or epigenetic pattern can fully explain etiology of panic disorder. Few genes and epigenetic patterns, however, showed strong association with panic disorder compared to healthy controls. Considering its multivariable background, further studies with larger populations can confirm current results and clarify etiologies of panic disorder.

Whole genome association study of brain-wide imaging phenotypes: A study of the ping cohort

Neuropsychological disorders have a biological basis rooted in brain function, and neuroimaging data are expected to better illuminate the complex genetic basis of neuropsychological disorders. Because they are biological measures, neuroimaging data avoid biases arising from clinical diagnostic criteria that are subject to human understanding and interpretation. A challenge with analyzing neuroimaging data is their high dimensionality and complex spatial relationships. To tackle this challenge, we introduced a novel distance covariance tests that can assess the association between genetic markers and multivariate diffusion tensor imaging measurements, and analyzed a genome-wide association study (GWAS) dataset collected by the Pediatric Imaging, Neurocognition, and Genetics (PING) study. We also considered existing approaches as comparisons. Our results showed that, after correcting for multiplicity, distance covariance tests of the multivariate phenotype yield significantly greater power at detecting genetic markers affecting brain structure than standard mass univariate GWAS of individual neuroimaging biomarkers. Our results underscore the usefulness of utilizing the distance covariance to incorporate neuroimaging data in GWAS.

TMEM132: an ancient architecture of cohesin and immunoglobulin domains define a new family of neural adhesion molecules

Summary

The molecular functions of TMEM132 genes remain poorly understood and under-investigated despite their mutations associated with non-syndromic hearing loss, panic disorder and cancer. Here we show the full domain architecture of human TMEM132 family proteins solved using in-depth sequence and structural analysis. We reveal them to be five previously unappreciated cell adhesion molecules whose domain architecture has an early holozoan origin prior to the emergence of choanoflagellates and metazoa. The extra-cellular portions of TMEM132 proteins contain five conserved domains including three tandem immunoglobulin domains, and a cohesin domain homologue, the first such domain found in animals. These findings strongly predict a cellular adhesion function for TMEM132 family, connecting the extracellular medium with the intracellular actin cytoskeleton.

Contact

luis.sanchez-pulido@igmm.ed.ac.uk.

Supplementary information

Supplementary data are available at Bioinformatics online.

Whole-exome sequencing and gene-based rare variant association tests suggest that PLA2G4E might be a risk gene for panic disorder

Panic disorder (PD) is characterized by recurrent and unexpected panic attacks, subsequent anticipatory anxiety, and phobic avoidance. Recent epidemiological and genetic studies have revealed that genetic factors contribute to the pathogenesis of PD. We performed whole-exome sequencing on one Japanese family, including multiple patients with panic disorder, which identified seven rare protein-altering variants. We then screened these genes in a Japanese PD case-control group (384 sporadic PD patients and 571 controls), resulting in the detection of three novel single nucleotide variants as potential candidates for PD (chr15: 42631993, T>C in GANC; chr15: 42342861, G>T in PLA2G4E; chr20: 3641457, G>C in GFRA4). Statistical analyses of these three genes showed that PLA2G4E yielded the lowest p value in gene-based rare variant association tests by Efficient and Parallelizable Association Container Toolbox algorithms; however, the p value did not reach the significance threshold in the Japanese. Likewise, in a German case-control study (96 sporadic PD patients and 96 controls), PLA2G4E showed the lowest p value but again did not reach the significance threshold. In conclusion, we failed to find any significant variants or genes responsible for the development of PD. Nonetheless, our results still leave open the possibility that rare protein-altering variants in PLA2G4E contribute to the risk of PD, considering the function of this gene.

Polymorphism in Tmem132d regulates expression and anxiety-related behavior through binding of RNA polymerase II complex

TMEM132D is a candidate gene, where risk genotypes have been associated with anxiety severity along with higher mRNA expression in the frontal cortex of panic disorder patients. Concurrently, in a high (HAB) and low (LAB) trait anxiety mouse model, Tmem132d was found to show increased expression in the anterior cingulate cortex (aCC) of HAB as compared to LAB mice. To understand the molecular underpinnings underlying the differential expression, we sequenced the gene and found two single-nucleotide polymorphisms (SNPs) in the promoter differing between both lines which could explain the observed mRNA expression profiles using gene reporter assays. In addition, there was no difference in basal DNA methylation in the CpG Island that encompasses the HAB vs. LAB Tmem132d promoter region. Furthermore, we found significantly higher binding of RNA polymerase II (POLR2A) to the proximal HAB-specific SNP (rs233264624) than the corresponding LAB locus in an oligonucleotide pull-down assay, suggesting increased transcription. Virus mediated overexpression of Tmem132d in the aCC of C57BL/6 J mice could confirm its role in mediating an anxiogenic phenotype. To model gene-environmental interactions, HAB mice exposed to enriched environment (HAB-EE) responded with decreased anxiety levels but, had enhanced Tmem132d mRNA expression as compared to standard-housed HAB (HAB-SH) mice. While LAB mice subjected to unpredictable chronic mild stress (LAB-UCMS) exhibited higher anxiety levels and had lower mRNA expression compared to standard-housed LAB (LAB-SH) mice. Chromatin immunoprecipitation revealed significantly higher binding of POLR2A to rs233264624 in HAB-EE, while LAB-UCMS had lower POLR2A binding at this locus, thus explaining the enhanced or attenuated expression of Tmem132d compared to their respective SH controls. To further investigate gene-environment interactions, DNA methylation was assessed using Illumina 450 K BeadChip in 74 panic disorder patients. Significant methylation differences were observed in two CpGs (cg26322591 and cg03283235) located in TMEM132D depending on the number of positive life events supporting the results of an influence of positive environmental cues on regulation of Tmem132d expression in mice.

Anxiety Associated Increased CpG Methylation in the Promoter of Asb1: A Translational Approach Evidenced by Epidemiological and Clinical Studies and a Murine Model

Epigenetic regulation in anxiety is suggested, but evidence from large studies is needed. We conducted an epigenome-wide association study (EWAS) on anxiety in a population-based cohort and validated our finding in a clinical cohort as well as a murine model. In the KORA cohort, participants (n=1522, age 32-72 years) were administered the Generalized Anxiety Disorder (GAD-7) instrument, whole blood DNA methylation was measured (Illumina 450K BeadChip), and circulating levels of hs-CRP and IL-18 were assessed in the association between anxiety and methylation. DNA methylation was measured using the same instrument in a study of patients with anxiety disorders recruited at the Max Planck Institute of Psychiatry (MPIP, 131 non-medicated cases and 169 controls). To expand our mechanistic understanding, these findings were reverse translated in a mouse model of acute social defeat stress. In the KORA study, participants were classified according to mild, moderate, or severe levels of anxiety (29.4%/6.0%/1.5%, respectively). Severe anxiety was associated with 48.5% increased methylation at a single CpG site (cg12701571) located in the promoter of the gene encoding Asb1 (β-coefficient=0.56 standard error (SE)=0.10, p (Bonferroni)=0.005), a protein hypothetically involved in regulation of cytokine signaling. An interaction between IL-18 and severe anxiety with methylation of this CpG cite showed a tendency towards significance in the total population (p=0.083) and a significant interaction among women (p=0.014). Methylation of the same CpG was positively associated with Panic and Agoraphobia scale (PAS) scores (β=0.005, SE=0.002, p=0.021, n=131) among cases in the MPIP study. In a murine model of acute social defeat stress, Asb1 gene expression was significantly upregulated in a tissue-specific manner (p=0.006), which correlated with upregulation of the neuroimmunomodulating cytokine interleukin 1 beta. Our findings suggest epigenetic regulation of the stress-responsive Asb1 gene in anxiety-related phenotypes. Further studies are necessary to elucidate the causal direction of this association and the potential role of Asb1-mediated immune dysregulation in anxiety disorders.

"DNA Methylation signatures in panic disorder"

Panic disorder (PD) affects about four million Europeans, with women affected twice as likely as men, causing substantial suffering and high economic costs. The etiopathogenesis of PD remains largely unknown, but both genetic and environmental factors contribute to risk. An epigenome-wide association study (EWAS) was conducted to compare medication-free PD patients (n = 89) with healthy controls (n = 76) stratified by gender. Replication was sought in an independent sample (131 cases, 169 controls) and functional analyses were conducted in a third sample (N = 71). DNA methylation was assessed in whole blood using the Infinium HumanMethylation450 BeadChip. One genome-wide association surviving FDR of 5% (cg07308824, P = 1.094 × 10-7, P-adj = 0.046) was identified in female PD patients (N = 49) compared to controls (N = 48). The same locus, located in an enhancer region of the HECA gene, was also hypermethylated in female PD patients in the replication sample (P = 0.035) and the significance of the association improved in the meta-analysis (P-adj = 0.004). Methylation at this CpG site was associated with HECA mRNA expression in another independent female sample (N = 71) both at baseline (P = 0.046) and after induction by dexamethasone (P = 0.029). Of 15 candidates, 5 previously reported as associated with PD or anxiety traits also showed differences in DNA methylation after gene-wise correction and included SGK1, FHIT, ADCYAP1, HTR1A, HTR2A. Our study examines epigenome-wide differences in peripheral blood for PD patients. Our results point to possible sex-specific methylation changes in the HECA gene for PD but overall highlight that this disorder is not associated with extensive changes in DNA methylation in peripheral blood.

Association between genes on chromosome 19p13.2 and panic disorder

Panic disorder (PD) is a severe and disabling mental disorder, which is moderately heritable. In a previous study, we carried out a genome-wide association study using patients with PD and control individuals from the isolated population of the Faroe Islands and identified chromosome 19p13.2 as a candidate region. To further investigate this chromosomal region for association with PD, we analysed eight single nucleotide polymorphisms (SNPs) in three candidate genes - small-nuclear RNA activating complex, polypeptide 2 (SNAPC2), mitogen-activated protein kinase kinase 7 (MAP2K7) and leucine-rich repeat containing 8 family, member E (LRRC8E) - these genes have previously been directly or indirectly implicated in other mental disorders. A total of 511 patients with PD and 1029 healthy control individuals from the Faroe Islands, Denmark and Germany were included in the current study. SNPs covering the gene region of SNAPC2, MAP2K7 and LRRC8E were genotyped and tested for association with PD. In the Faroese cohort, rs7788 within SNAPC2 was significantly associated with PD, whereas rs3745383 within LRRC8E was nominally associated. No association was observed between the analysed SNPs and PD in the Danish cohorts. In the German women, we observed a nominal association between rs4804833 within MAP2K7 and PD. We present further evidence that chromosome 19p13.2 may harbour candidate genes that contribute towards the risk of developing PD. Moreover, the implication of the associated genes in other mental disorders may indicate shared genetic susceptibility between mental disorders. We show that associated variants may be sex specific, indicating the importance of carrying out a sex-specific association analysis of PD.

GLRB allelic variation associated with agoraphobic cognitions, increased startle response and fear network activation: a potential neurogenetic pathway to panic disorder

The molecular genetics of panic disorder (PD) with and without agoraphobia (AG) are still largely unknown and progress is hampered by small sample sizes. We therefore performed a genome-wide association study with a dimensional, PD/AG-related anxiety phenotype based on the Agoraphobia Cognition Questionnaire (ACQ) in a sample of 1370 healthy German volunteers of the CRC TRR58 MEGA study wave 1. A genome-wide significant association was found between ACQ and single non-coding nucleotide variants of the GLRB gene (rs78726293, P=3.3 × 10^-8; rs191260602, P=3.9 × 10^-8). We followed up on this finding in a larger dimensional ACQ sample (N=2547) and in independent samples with a dichotomous AG phenotype based on the Symptoms Checklist (SCL-90; N=3845) and a case-control sample with the categorical phenotype PD/AG (N_combined =1012) obtaining highly significant P-values also for GLRB single-nucleotide variants rs17035816 (P=3.8 × 10^-4) and rs7688285 (P=7.6 × 10^-5). GLRB gene expression was found to be modulated by rs7688285 in brain tissue, as well as cell culture. Analyses of intermediate PD/AG phenotypes demonstrated increased startle reflex and increased fear network, as well as general sensory activation by GLRB risk gene variants rs78726293, rs191260602, rs17035816 and rs7688285. Partial Glrb knockout mice demonstrated an agoraphobic phenotype. In conjunction with the clinical observation that rare coding GLRB gene mutations are associated with the neurological disorder hyperekplexia characterized by a generalized startle reaction and agoraphobic behavior, our data provide evidence that non-coding, although functional GLRB gene polymorphisms may predispose to PD by increasing startle response and agoraphobic cognitions.

Modulation of defensive reactivity by GLRB allelic variation: converging evidence from an intermediate phenotype approach

Representing a phylogenetically old and very basic mechanism of inhibitory neurotransmission, glycine receptors have been implicated in the modulation of behavioral components underlying defensive responding toward threat. As one of the first findings being confirmed by genome-wide association studies for the phenotype of panic disorder and agoraphobia, allelic variation in a gene coding for the glycine receptor beta subunit (GLRB) has recently been associated with increased neural fear network activation and enhanced acoustic startle reflexes. On the basis of two independent healthy control samples, we here aimed to further explore the functional significance of the GLRB genotype (rs7688285) by employing an intermediate phenotype approach. We focused on the phenotype of defensive system reactivity across the levels of brain function, structure, and physiology. Converging evidence across both samples was found for increased neurofunctional activation in the (anterior) insular cortex in GLRB risk allele carriers and altered fear conditioning as a function of genotype. The robustness of GLRB effects is demonstrated by consistent findings across different experimental fear conditioning paradigms and recording sites. Altogether, findings provide translational evidence for glycine neurotransmission as a modulator of the brain's evolutionary old dynamic defensive system and provide further support for a strong, biologically plausible candidate intermediate phenotype of defensive reactivity. As such, glycine-dependent neurotransmission may open up new avenues for mechanistic research on the etiopathogenesis of fear and anxiety disorders.

Genome-wide association study for behavior, type traits, and muscular development in Charolais beef cattle

Behavior, type traits, and muscular development are of interest for beef cattle breeding. Genome-wide association studies (GWAS) enable the identification of candidate genes, which enables gene-based selection and provides insight in the genetic architecture of these traits. The objective of the current study was to perform a GWAS for 3 behavior traits, 12 type traits, and muscular development in Charolais cattle. Behavior traits, including aggressiveness at parturition, aggressiveness during gestation period, and maternal care, were scored by farmers. Type traits, including udder conformation, teat, feet and legs, and locomotion, were scored by trained classifiers. Data used in the GWAS consisted of 3,274 cows with phenotypic records and genotyping information for 44,930 SNP. When SNP had a false discovery rate (FDR) smaller than 0.05, they were referred to as significant. When SNP had a FDR between 0.05 and 0.20, they were referred to as suggestive. Four significant and 12 suggestive regions were detected for aggressiveness during gestation, maternal care, udder balance, teat thinness, teat length, foot angle, foot depth, and locomotion. These 4 significant and 12 suggestive regions were not supported by other significant SNP in close proximity. No SNP with major effects were detected for behavior and type traits, and SNP associations for these traits were spread across the genome, suggesting that behavior and type traits were influenced by many genes, each explaining a small part of genetic variance. The GWAS identified 1 region on chromosome 2 significantly associated with muscular development, which included the myostatin gene (), which is known to affect muscularity. No other regions associated with muscular development were found. Results showed that the myostatin region associated with muscular development had pleiotropic effects on udder volume, teat thinness, rear leg, and leg angle.

BDNF Val66Met Genotype Interacts With a History of Simulated Stress Exposure to Regulate Sensorimotor Gating and Startle Reactivity

Reduced expression of Brain-Derived Neurotrophic Factor (BDNF) has been implicated in the pathophysiology of schizophrenia. The BDNF Val66Met polymorphism, which results in deficient activity-dependent secretion of BDNF, is associated with clinical features of schizophrenia. We investigated the effect of this polymorphism on Prepulse Inhibition (PPI), a translational model of sensorimotor gating which is disrupted in schizophrenia. We utilized humanized BDNFVal66Met (hBDNFVal66Met) mice which have been modified to carry the Val66Met polymorphism, as well as express humanized BDNF in vivo. We also studied the long-term effect of chronic corticosterone (CORT) exposure in these animals as a model of history of stress. PPI was assessed at 30ms and 100ms interstimulus intervals (ISI). Analysis of PPI at the commonly used 100ms ISI identified that, irrespective of CORT treatment, the hBDNFVal/Met genotype was associated with significantly reduced PPI. In contrast, PPI was not different between hBDNFMet/Met and hBDNFVal/Val genotype mice. At the 30ms ISI, CORT treatment selectively disrupted sensorimotor gating of hBDNFVal/Met heterozygote mice but not hBDNFVal/Val or hBDNFMet/Met mice. Analysis of startle reactivity revealed that chronic CORT reduced startle reactivity of hBDNFVal/Val male mice by 51%. However, this was independent of the effect of CORT on PPI. In summary, we provide evidence of a distinct BDNFVal66Met heterozygote-specific phenotype using the sensorimotor gating endophenotype of schizophrenia. These data have important implications for clinical studies where, if possible, the BDNFVal/Met heterozygote genotype should be distinguished from the BDNFMet/Met genotype.

Functional neuroanatomy in panic disorder: Status quo of the research

AIM

To provide an overview of the current research in the functional neuroanatomy of panic disorder.

METHODS

Panic disorder (PD) is a frequent psychiatric disease. Gorman et al (1989; 2000) proposed a comprehensive neuroanatomical model of PD, which suggested that fear- and anxiety-related responses are mediated by a so-called “fear network” which is centered in the amygdala and includes the hippocampus, thalamus, hypothalamus, periaqueductal gray region, locus coeruleus and other brainstem sites. We performed a systematic search by the electronic database PubMed. Thereby, the main focus was laid on recent neurofunctional, neurostructural, and neurochemical studies (from the period between January 2012 and April 2016). Within this frame, special attention was given to the emerging field of imaging genetics.

RESULTS

We noted that many neuroimaging studies have reinforced the role of the “fear network” regions in the pathophysiology of panic disorder. However, recent functional studies suggest abnormal activation mainly in an extended fear network comprising brainstem, anterior and midcingulate cortex (ACC and MCC), insula, and lateral as well as medial parts of the prefrontal cortex. Interestingly, differences in the amygdala activation were not as consistently reported as one would predict from the hypothesis of Gorman et al (2000). Indeed, amygdala hyperactivation seems to strongly depend on stimuli and experimental paradigms, sample heterogeneity and size, as well as on limitations of neuroimaging techniques. Advanced neurochemical studies have substantiated the major role of serotonergic, noradrenergic and glutamatergic neurotransmission in the pathophysiology of PD. However, alterations of GABAergic function in PD are still a matter of debate and also their specificity remains questionable. A promising new research approach is “imaging genetics”. Imaging genetic studies are designed to evaluate the impact of genetic variations (polymorphisms) on cerebral function in regions critical for PD. Most recently, imaging genetic studies have not only confirmed the importance of serotonergic and noradrenergic transmission in the etiology of PD but also indicated the significance of neuropeptide S receptor, CRH receptor, human TransMEMbrane protein (TMEM123D), and amiloride-sensitive cation channel 2 (ACCN2) genes.

CONCLUSION

In light of these findings it is conceivable that in the near future this research will lead to the development of clinically useful tools like predictive biomarkers or novel treatment options.

Genome-wide association study of generalized anxiety symptoms in the Hispanic Community Health Study/Study of Latinos

Although generalized anxiety disorder (GAD) is heritable and aggregates in families, no genomic loci associated with GAD have been reported. We aimed to discover potential loci by conducting a genome-wide analysis of GAD symptoms in a large, population-based sample of Hispanic/Latino adults. Data came from 12,282 participants (aged 18-74) in the Hispanic Community Health Study/Study of Latinos. Using a shortened Spielberger Trait Anxiety measure, we analyzed the following: (i) a GAD symptoms score restricted to the three items tapping diagnostic features of GAD as defined by DSM-V; and (ii) a total trait anxiety score based on summing responses to all ten items. We first calculated the heritability due to common variants (h2SNP ) and then conducted a genome-wide association study (GWAS) of GAD symptoms. Replication was attempted in three independent Hispanic cohorts (Multi-Ethnic Study of Atherosclerosis, Women's Health Initiative, Army STARRS). The GAD symptoms score showed evidence of modest heritability (7.2%; P = 0.03), while the total trait anxiety score did not (4.97%; P = 0.20). One genotyped SNP (rs78602344) intronic to thrombospondin 2 (THBS2) was nominally associated (P = 5.28 × 10-8 ) in the primary analysis adjusting for psychiatric medication use and significantly associated with the GAD symptoms score in the analysis excluding medication users (P = 4.18 × 10-8 ). However, meta-analysis of the replication samples did not support this association. Although we identified a genome-wide significant locus in this sample, we were unable to replicate this finding. Evidence for heritability was also only detected for GAD symptoms, and not the trait anxiety measure, suggesting differential genetic influences within the domain of trait anxiety. © 2016 Wiley Periodicals, Inc.

Epigenome-wide association study of DNA methylation in panic disorder

Background

Panic disorder (PD) is considered to be a multifactorial disorder emerging from interactions among multiple genetic and environmental factors. To date, although genetic studies reported several susceptibility genes with PD, few of them were replicated and the pathogenesis of PD remains to be clarified. Epigenetics is considered to play an important role in etiology of complex traits and diseases, and DNA methylation is one of the major forms of epigenetic modifications. In this study, we performed an epigenome-wide association study of PD using DNA methylation arrays so as to investigate the possibility that different levels of DNA methylation might be associated with PD.

Methods

The DNA methylation levels of CpG sites across the genome were examined with genomic DNA samples (PD, N = 48, control, N = 48) extracted from peripheral blood. Methylation arrays were used for the analysis. β values, which represent the levels of DNA methylation, were normalized via an appropriate pipeline. Then, β values were converted to M values via the logit transformation for epigenome-wide association study. The relationship between each DNA methylation site and PD was assessed by linear regression analysis with adjustments for the effects of leukocyte subsets.

Results

Forty CpG sites showed significant association with PD at 5% FDR correction, though the differences of the DNA methylation levels were relatively small. Most of the significant CpG sites (37/40 CpG sites) were located in or around CpG islands. Many of the significant CpG sites (27/40 CpG sites) were located upstream of genes, and all such CpG sites with the exception of two were hypomethylated in PD subjects. A pathway analysis on the genes annotated to the significant CpG sites identified several pathways, including “positive regulation of lymphocyte activation.”

Conclusions

Although future studies with larger number of samples are necessary to confirm the small DNA methylation abnormalities associated with PD, there is a possibility that several CpG sites might be associated, together as a group, with PD.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-016-0307-1) contains supplementary material, which is available to authorized users.

Whole-exome sequencing implicates DGKH as a risk gene for panic disorder in the Faroese population

The demographic history of the isolated population of the Faroe Islands may have induced enrichment of variants rarely seen in outbred European populations, including enrichment of risk variants for panic disorder (PD). PD is a common mental disorder, characterized by recurring and unprovoked panic attacks, and genetic factors have been estimated to explain around 40% of the risk. In this study the potential enrichment of PD risk variants was explored based on whole-exome sequencing of 54 patients with PD and 211 control individuals from the Faroese population. No genome-wide significant associations were found, however several single variants and genes showed strong association with PD, where DGKH was found to be the strongest PD associated gene. Interestingly DGKH has previously demonstrated genome-wide significant association with bipolar disorder as well as evidence of association to other mental disorders. Additionally, we found an enrichment of PD risk variants in the Faroese population; variants with otherwise low frequency in more outbreed European populations. © 2016 Wiley Periodicals, Inc.

Meta-analysis of genome-wide association studies of anxiety disorders

Anxiety disorders (ADs), namely generalized AD, panic disorder and phobias, are common, etiologically complex conditions with a partially genetic basis. Despite differing on diagnostic definitions based on clinical presentation, ADs likely represent various expressions of an underlying common diathesis of abnormal regulation of basic threat-response systems. We conducted genome-wide association analyses in nine samples of European ancestry from seven large, independent studies. To identify genetic variants contributing to genetic susceptibility shared across interview-generated DSM-based ADs, we applied two phenotypic approaches: (1) comparisons between categorical AD cases and supernormal controls, and (2) quantitative phenotypic factor scores (FS) derived from a multivariate analysis combining information across the clinical phenotypes. We used logistic and linear regression, respectively, to analyze the association between these phenotypes and genome-wide single nucleotide polymorphisms. Meta-analysis for each phenotype combined results across the nine samples for over 18 000 unrelated individuals. Each meta-analysis identified a different genome-wide significant region, with the following markers showing the strongest association: for case-control contrasts, rs1709393 located in an uncharacterized non-coding RNA locus on chromosomal band 3q12.3 (P=1.65 × 10(-8)); for FS, rs1067327 within CAMKMT encoding the calmodulin-lysine N-methyltransferase on chromosomal band 2p21 (P=2.86 × 10(-9)). Independent replication and further exploration of these findings are needed to more fully understand the role of these variants in risk and expression of ADs.

A powerful phenotype for gene-finding studies derived from trajectory analyses of symptoms of anxiety and depression between age seven and 18

To investigate the utility of longitudinal data in genetic analyses of symptoms of anxiety and depression, we assessed individual differences between age 7 and 18 using growth mixture models, and investigated the genetic and non-genetic factors contributing to the trajectories. Mothers of 7,706 girl and 7,418 boy twins from the Netherlands Twin Register rated the anxious depression scale (SxAnxDep) of the Child Behavior Check List (CBCL) at age 7, 10, and 12 years. Two thousand seven hundred and six girl and 1,856 boy twins completed the Youth Self Report (YSR) at age 14, 16, and 18. While individual trajectories varied considerably, these differences were largely idiosyncratic and could not be grouped into separate latent classes with class-specific average growth curves. The intercept, which reflects the individuals' baseline level of SxAnxDep across time, explained 55-58% of the total phenotypic variance. The slope factor, which captures a common average trend over time, did not explain variance in the phenotype. This finding also underlines the high level of idiosyncrasy of trajectories that lack a common longitudinal structure. The analyses of twin data showed that the random intercept factor of SxAnxDep during childhood and during adolescence is considerably more heritable than the observations at any single age, namely between 60% and 84%. One explanation is that different factors contribute to the level of symptoms of anxiety and depression at any given time point, including temporary events and emotions. When considering baseline stability, these temporary influences average out, with the result of a more reliable and more heritable phenotype. © 2015 Wiley Periodicals, Inc.

Biological markers for anxiety disorders, OCD and PTSD - a consensus statement. Part I: Neuroimaging and genetics

OBJECTIVES

Biomarkers are defined as anatomical, biochemical or physiological traits that are specific to certain disorders or syndromes. The objective of this paper is to summarise the current knowledge of biomarkers for anxiety disorders, obsessive-compulsive disorder (OCD) and post-traumatic stress disorder (PTSD).

METHODS

Findings in biomarker research were reviewed by a task force of international experts in the field, consisting of members of the World Federation of Societies for Biological Psychiatry Task Force on Biological Markers and of the European College of Neuropsychopharmacology Anxiety Disorders Research Network.

RESULTS

The present article (Part I) summarises findings on potential biomarkers in neuroimaging studies, including structural brain morphology, functional magnetic resonance imaging and techniques for measuring metabolic changes, including positron emission tomography and others. Furthermore, this review reports on the clinical and molecular genetic findings of family, twin, linkage, association and genome-wide association studies. Part II of the review focuses on neurochemistry, neurophysiology and neurocognition.

CONCLUSIONS

Although at present, none of the putative biomarkers is sufficient and specific as a diagnostic tool, an abundance of high-quality research has accumulated that will improve our understanding of the neurobiological causes of anxiety disorders, OCD and PTSD.

Genome-wide significant loci for addiction and anxiety

BACKGROUND

Psychiatric comorbidity is common among individuals with addictive disorders, with patients frequently suffering from anxiety disorders. While the genetic architecture of comorbid addictive and anxiety disorders remains unclear, elucidating the genes involved could provide important insights into the underlying etiology.

METHODS

Here we examine a sample of 1284 Mexican-Americans from randomly selected extended pedigrees. Variance decomposition methods were used to examine the role of genetics in addiction phenotypes (lifetime history of alcohol dependence, drug dependence or chronic smoking) and various forms of clinically relevant anxiety. Genome-wide univariate and bivariate linkage scans were conducted to localize the chromosomal regions influencing these traits.

RESULTS

Addiction phenotypes and anxiety were shown to be heritable and univariate genome-wide linkage scans revealed significant quantitative trait loci for drug dependence (14q13.2-q21.2, LOD=3.322) and a broad anxiety phenotype (12q24.32-q24.33, LOD=2.918). Significant positive genetic correlations were observed between anxiety and each of the addiction subtypes (ρg=0.550-0.655) and further investigation with bivariate linkage analyses identified significant pleiotropic signals for alcohol dependence-anxiety (9q33.1-q33.2, LOD=3.054) and drug dependence-anxiety (18p11.23-p11.22, LOD=3.425).

CONCLUSIONS

This study confirms the shared genetic underpinnings of addiction and anxiety and identifies genomic loci involved in the etiology of these comorbid disorders. The linkage signal for anxiety on 12q24 spans the location of TMEM132D, an emerging gene of interest from previous GWAS of anxiety traits, whilst the bivariate linkage signal identified for anxiety-alcohol on 9q33 peak coincides with a region where rare CNVs have been associated with psychiatric disorders. Other signals identified implicate novel regions of the genome in addiction genetics.

Candidate genes in panic disorder: meta-analyses of 23 common variants in major anxiogenic pathways

The utilization of molecular genetics approaches in examination of panic disorder (PD) has implicated several variants as potential susceptibility factors for panicogenesis. However, the identification of robust PD susceptibility genes has been complicated by phenotypic diversity, underpowered association studies and ancestry-specific effects. In the present study, we performed a succinct review of case-control association studies published prior to April 2015. Meta-analyses were performed for candidate gene variants examined in at least three studies using the Cochrane Mantel-Haenszel fixed-effect model. Secondary analyses were also performed to assess the influences of sex, agoraphobia co-morbidity and ancestry-specific effects on panicogenesis. Meta-analyses were performed on 23 variants in 20 PD candidate genes. Significant associations after correction for multiple testing were observed for three variants, TMEM132D rs7370927 (T allele: odds ratio (OR)=1.27, 95% confidence interval (CI): 1.15-1.40, P=2.49 × 10(-6)), rs11060369 (CC genotype: OR=0.65, 95% CI: 0.53-0.79, P=1.81 × 10(-5)) and COMT rs4680 (Val (G) allele: OR=1.27, 95% CI: 1.14-1.42, P=2.49 × 10(-5)) in studies with samples of European ancestry. Nominal associations that did not survive correction for multiple testing were observed for NPSR1 rs324891 (T allele: OR=1.22, 95% CI: 1.07-1.38, P=0.002), TPH1 rs1800532 (AA genotype: OR=1.46, 95% CI: 1.14-1.89, P=0.003) and HTR2A rs6313 (T allele: OR=1.19, 95% CI: 1.07-1.33, P=0.002) in studies with samples of European ancestry and for MAOA-uVNTR in female PD (low-active alleles: OR=1.21, 95% CI: 1.07-1.38, P=0.004). No significant associations were observed in the secondary analyses considering sex, agoraphobia co-morbidity and studies with samples of Asian ancestry. Although these findings highlight a few associations, PD likely involves genetic variation in a multitude of biological pathways that is diverse among populations. Future studies must incorporate larger sample sizes and genome-wide approaches to further quantify the observed genetic variation among populations and subphenotypes of PD.

Polymorphisms in the TMEM132D region are associated with panic disorder in HLA-DRB1*13:02-negative individuals of a Japanese population

We herein report an association between TMEM132D and panic disorder (PD) in a Japanese population, evaluating the effects of HLA-DRB1*13:02, which we previously reported as a susceptibility genetic factor for PD. SNPs in TMEM132D showed significant associations with PD in subjects without HLA-DRB1*13:02 (rs4759997; P=5.02×10⁻⁶, odds ratio=1.50) but not in those with the HLA allele. TMEM132D might have a role in the development of PD in subjects without HLA-DRB1*13:02.