Pleiotropic contribution of rbfox1 to psychiatric and neurodevelopmental phenotypes in two zebrafish models

RBFOX1 is a highly pleiotropic gene that contributes to several psychiatric and neurodevelopmental disorders. Both rare and common variants in RBFOX1 have been associated with several psychiatric conditions, but the mechanisms underlying the pleiotropic effects of RBFOX1 are not yet understood. Here we found that, in zebrafish, rbfox1 is expressed in spinal cord, mid- and hindbrain during developmental stages. In adults, expression is restricted to specific areas of the brain, including telencephalic and diencephalic regions with an important role in receiving and processing sensory information and in directing behaviour. To investigate the contribution of rbfox1 to behaviour, we used rbfox1sa15940, a zebrafish mutant line with TL background. We found that rbfox1sa15940 mutants present hyperactivity, thigmotaxis, decreased freezing behaviour and altered social behaviour. We repeated these behavioural tests in a second rbfox1 mutant line with a different genetic background (TU), rbfox1del19, and found that rbfox1 deficiency affects behaviour similarly in this line, although there were some differences. rbfox1del19 mutants present similar thigmotaxis, but stronger alterations in social behaviour and lower levels of hyperactivity than rbfox1sa15940 fish. Taken together, these results suggest that mutations in rbfox1 lead to multiple behavioural changes in zebrafish that might be modulated by environmental, epigenetic and genetic background effects, and that resemble phenotypic alterations present in Rbfox1-deficient mice and in patients with different psychiatric conditions. Our study, thus, highlights the evolutionary conservation of rbfox1 function in behaviour and paves the way to further investigate the mechanisms underlying rbfox1 pleiotropy on the onset of neurodevelopmental and psychiatric disorders.

Decreased Brain Serotonin in rbfox1 Mutant Zebrafish and Partial Reversion of Behavioural Alterations by the SSRI Fluoxetine

RBFOX1 functions as a master regulator of thousands of genes, exerting a pleiotropic effect on numerous neurodevelopmental and psychiatric disorders. A potential mechanism by which RBFOX1 may impact these disorders is through its modulation of serotonergic neurotransmission, a common target for pharmacological intervention in psychiatric conditions linked to RBFOX1. However, the precise effects of RBFOX1 on the serotonergic system remain largely unexplored. Here we show that homozygous rbfox1sa15940 zebrafish, which express a shorter, aberrant rbfox1 mRNA, have significantly reduced serotonin levels in telencephalon and diencephalon. We observed that the acute administration of fluoxetine partially reverses the associated behavioural alterations. The hyperactive phenotype and altered shoaling behaviour of the rbfox1sa15940/sa15940 zebrafish could be reversed with acute fluoxetine exposure in the Open Field and the Shoaling test, respectively. However, in the other paradigms, hyperactivity was not diminished, suggesting a distinct intrinsic motivation for locomotion in the different paradigms. Acute fluoxetine exposure did not reverse the alterations observed in the aggression and social novelty tests, suggesting the involvement of other neurological mechanisms in these behaviours. These findings underscore the importance of investigating the intricate working mechanisms of RBFOX1 in neurodevelopmental and psychiatric disorders to gain a better understanding of the associated disorders along with their pharmacological treatment.

Exploring small non-coding RNAs as blood-based biomarkers to predict Alzheimer's disease

BACKGROUND

Alzheimer's disease (AD) diagnosis relies on clinical symptoms complemented with biological biomarkers, the Amyloid Tau Neurodegeneration (ATN) framework. Small non-coding RNA (sncRNA) in the blood have emerged as potential predictors of AD. We identified sncRNA signatures specific to ATN and AD, and evaluated both their contribution to improving AD conversion prediction beyond ATN alone.

METHODS

This nested case-control study was conducted within the ACE cohort and included MCI patients matched by sex. Patients free of type 2 diabetes underwent cerebrospinal fluid (CSF) and plasma collection and were followed-up for a median of 2.45-years. Plasma sncRNAs were profiled using small RNA-sequencing. Conditional logistic and Cox regression analyses with elastic net penalties were performed to identify sncRNA signatures for A+(T|N)+ and AD. Weighted scores were computed using cross-validation, and the association of these scores with AD risk was assessed using multivariable Cox regression models. Gene ontology (GO) and Kyoto encyclopaedia of genes and genomes (KEGG) enrichment analysis of the identified signatures were performed.

RESULTS

The study sample consisted of 192 patients, including 96 A+(T|N)+ and 96 A-T-N- patients. We constructed a classification model based on a 6-miRNAs signature for ATN. The model could classify MCI patients into A-T-N- and A+(T|N)+ groups with an area under the curve of 0.7335 (95% CI, 0.7327 to 0.7342). However, the addition of the model to conventional risk factors did not improve the prediction of AD beyond the conventional model plus ATN status (C-statistic: 0.805 [95% CI, 0.758 to 0.852] compared to 0.829 [95% CI, 0.786, 0.872]). The AD-related 15-sncRNAs signature exhibited better predictive performance than the conventional model plus ATN status (C-statistic: 0.849 [95% CI, 0.808 to 0.890]). When ATN was included in this model, the prediction further improved to 0.875 (95% CI, 0.840 to 0.910). The miRNA-target interaction network and functional analysis, including GO and KEGG pathway enrichment analysis, suggested that the miRNAs in both signatures are involved in neuronal pathways associated with AD.

CONCLUSIONS

The AD-related sncRNA signature holds promise in predicting AD conversion, providing insights into early AD development and potential targets for prevention.

Three-Photon Infrared Stimulation of Endogenous Neuroreceptors in Vivo

To interrogate neural circuits and crack their codes, in vivo brain activity imaging must be combined with spatiotemporally precise stimulation in three dimensions using genetic or pharmacological specificity. This challenge requires deep penetration and focusing as provided by infrared light and multiphoton excitation, and has promoted two-photon photopharmacology and optogenetics. However, three-photon brain stimulation in vivo remains to be demonstrated. We report the regulation of neuronal activity in zebrafish larvae by three-photon excitation of a photoswitchable muscarinic agonist at 50 pM, a billion-fold lower concentration than used for uncaging, and with mid-infrared light of 1560 nm, the longest reported photoswitch wavelength. Robust, physiologically relevant photoresponses allow modulating brain activity in wild-type animals with spatiotemporal and pharmacological precision. Computational calculations predict that azobenzene-based ligands have high three-photon absorption cross-section and can be used directly with pulsed infrared light. The expansion of three-photon pharmacology will deeply impact basic neurobiology and neuromodulation phototherapies.

Pleiotropic contribution of rbfox1 to psychiatric and neurodevelopmental phenotypes in a zebrafish model

RBFOX1 is a highly pleiotropic gene that contributes to several psychiatric and neurodevelopmental disorders. Both rare and common variants in RBFOX1 have been associated with several psychiatric conditions, but the mechanisms underlying the pleiotropic effects of RBFOX1 are not yet understood. Here we found that, in zebrafish, rbfox1 is expressed in spinal cord, mid- and hindbrain during developmental stages. In adults, expression is restricted to specific areas of the brain, including telencephalic and diencephalic regions with an important role in receiving and processing sensory information and in directing behaviour. To investigate the effect of rbfox1 deficiency on behaviour, we used rbfox1sa15940, a rbfox1 loss-of-function line. We found that rbfox1sa15940 mutants present hyperactivity, thigmotaxis, decreased freezing behaviour and altered social behaviour. We repeated these behavioural tests in a second rbfox1 loss-of-function line with a different genetic background, rbfox1del19, and found that rbfox1 deficiency affects behaviour similarly in this line, although there were some differences. rbfox1del19 mutants present similar thigmotaxis, but stronger alterations in social behaviour and lower levels of hyperactivity than rbfox1sa15940 fish. Taken together, these results suggest that rbfox1 deficiency leads to multiple behavioural changes in zebrafish that might be modulated by environmental, epigenetic and genetic background effects, and that resemble phenotypic alterations present in Rbfox1-deficient mice and in patients with different psychiatric conditions. Our study thus highlights the evolutionary conservation of rbfox1 function in behaviour and paves the way to further investigate the mechanisms underlying rbfox1 pleiotropy on the onset of neurodevelopmental and psychiatric disorders.

The pleiotropic contribution of genes in dopaminergic and serotonergic pathways to addiction and related behavioral traits

Introduction

Co-occurrence of substance use disorders (SUD) and other behavioral conditions, such as stress-related, aggressive or risk-taking behaviors, in the same individual has been frequently described. As dopamine (DA) and serotonin (5-HT) have been previously identified as key neurotransmitters for some of these phenotypes, we explored the genetic contribution of these pathways to SUD and these comorbid phenotypes in order to better understand the genetic relationship between them.

Methods

We tested the association of 275 dopaminergic genes and 176 serotonergic genes with these phenotypes by performing gene-based, gene-set and transcriptome-wide association studies in 11 genome-wide association studies (GWAS) datasets on SUD and related behaviors.

Results

At the gene-wide level, 68 DA and 27 5-HT genes were found to be associated with at least one GWAS on SUD or related behavior. Among them, six genes had a pleiotropic effect, being associated with at least three phenotypes: ADH1C, ARNTL, CHRNA3, HPRT1, HTR1B and DRD2. Additionally, we found nominal associations between the DA gene sets and SUD, opioid use disorder, antisocial behavior, irritability and neuroticism, and between the 5-HT-core gene set and neuroticism. Predicted gene expression correlates in brain were also found for 19 DA or 5-HT genes.

Discussion

Our study shows a pleiotropic contribution of dopaminergic and serotonergic genes to addiction and related behaviors such as anxiety, irritability, neuroticism and risk-taking behavior, highlighting a role for DA genes, which could explain, in part, the co-occurrence of these phenotypes.

Genome-wide DNA methylation analysis in an antimigraine-treated preclinical model of cortical spreading depolarization

BACKGROUND

Cortical spreading depolarization, the cause of migraine aura, is a short-lasting depolarization wave that moves across the brain cortex, transiently suppressing neuronal activity. Prophylactic treatments for migraine, such as topiramate or valproate, reduce the number of cortical spreading depression events in rodents.

OBJECTIVE

To investigate whether cortical spreading depolarization with and without chronic treatment with topiramate or valproate affect the DNA methylation of the cortex.

METHODS

Sprague-Dawley rats were intraperitoneally injected with saline, topiramate or valproate for four weeks when cortical spreading depolarization were induced and genome-wide DNA methylation was performed in the cortex of six rats per group.

RESULTS

The DNA methylation profile of the cortex was significantly modified after cortical spreading depolarization, with and without topiramate or valproate. Interestingly, topiramate reduced by almost 50% the number of differentially methylated regions, whereas valproate increased them by 17%, when comparing to the non-treated group after cortical spreading depolarization induction. The majority of the differentially methylated regions lay within intragenic regions, and the analyses of functional group over-representation retrieved several enriched functions, including functions related to protein processing in the cortical spreading depolarization without treatment group; functions related to metabolic processes in the cortical spreading depolarization with topiramate group; and functions related to synapse and ErbB, MAPK or retrograde endocannabinoid signaling in the cortical spreading depolarization with valproate group.

CONCLUSIONS

Our results may provide insights into the underlying physiological mechanisms of migraine with aura and emphasize the role of epigenetics in migraine susceptibility.

The translational genetics of ADHD and related phenotypes in model organisms

Attention-deficit/hyperactivity disorder (ADHD) is a highly prevalent neurodevelopmental disorder resulting from the interaction between genetic and environmental risk factors. It is well known that ADHD co-occurs frequently with other psychiatric disorders due, in part, to shared genetics factors. Although many studies have contributed to delineate the genetic landscape of psychiatric disorders, their specific molecular underpinnings are still not fully understood. The use of animal models can help us to understand the role of specific genes and environmental stimuli-induced epigenetic modifications in the pathogenesis of ADHD and its comorbidities. The aim of this review is to provide an overview on the functional work performed in rodents, zebrafish and fruit fly and highlight the generated insights into the biology of ADHD, with a special focus on genetics and epigenetics. We also describe the behavioral tests that are available to study ADHD-relevant phenotypes and comorbid traits in these models. Furthermore, we have searched for new models to study ADHD and its comorbidities, which can be useful to test potential pharmacological treatments.

Identification of genetic variants influencing methylation in brain with pleiotropic effects on psychiatric disorders

Psychiatric disorders affect 29% of the global population at least once in the lifespan, and genetic studies have proved a shared genetic basis among them, although the underlying molecular mechanisms remain largely unknown. DNA methylation plays an important role in complex disorders and, remarkably, enrichment of common genetic variants influencing allele-specific methylation (ASM) has been reported among variants associated with specific psychiatric disorders. In the present study we assessed the contribution of ASM to a set of eight psychiatric disorders by combining genetic, epigenetic and expression data. We interrogated a list of 3896 ASM tagSNPs in the brain in the summary statistics of a cross-disorder GWAS meta-analysis of eight psychiatric disorders from the Psychiatric Genomics Consortium, including more than 162,000 cases and 276,000 controls. We identified 80 SNPs with pleiotropic effects on psychiatric disorders that show an opposite directional effect on methylation and gene expression. These SNPs converge on eight candidate genes: ZSCAN29, ZSCAN31, BTN3A2, DDAH2, HAPLN4, ARTN, FAM109B and NAGA. ZSCAN29 shows the broadest pleiotropic effects, showing associations with five out of eight psychiatric disorders considered, followed by ZSCAN31 and BTN3A2, associated with three disorders. All these genes overlap with CNVs related to cognitive phenotypes and psychiatric traits, they are expressed in the brain, and seven of them have previously been associated with specific psychiatric disorders, supporting our results. To sum up, our integrative functional genomics analysis identified eight psychiatric disease risk genes that impact a broad list of disorders and highlight an etiologic role of SNPs that influence DNA methylation and gene expression in the brain.

Comprehensive exploration of the genetic contribution of the dopaminergic and serotonergic pathways to psychiatric disorders

Psychiatric disorders are highly prevalent and display considerable clinical and genetic overlap. Dopaminergic and serotonergic neurotransmission have been shown to play an important role in many psychiatric disorders. Here we aim to assess the genetic contribution of these systems to eight psychiatric disorders (attention-deficit hyperactivity disorder (ADHD), anorexia nervosa (ANO), autism spectrum disorder (ASD), bipolar disorder (BIP), major depression (MD), obsessive-compulsive disorder (OCD), schizophrenia (SCZ) and Tourette's syndrome (TS)) using publicly available GWAS analyses performed by the Psychiatric Genomics Consortium that include more than 160,000 cases and 275,000 controls. To do so, we elaborated four different gene sets: two 'wide' selections for dopamine (DA) and for serotonin (SERT) using the Gene Ontology and KEGG pathways tools, and two'core' selections for the same systems, manually curated. At the gene level, we found 67 genes from the DA and/or SERT gene sets significantly associated with one of the studied disorders, and 12 of them were associated with two different disorders. Gene-set analysis revealed significant associations for ADHD and ASD with the wide DA gene set, for BIP with the wide SERT gene set, and for MD with the core SERT set. Interestingly, interrogation of a cross-disorder GWAS meta-analysis of the eight psychiatric conditions displayed association with the wide DA gene set. To our knowledge, this is the first systematic examination of genes encoding proteins essential to the function of these two neurotransmitter systems in these disorders. Our results support a pleiotropic contribution of the dopaminergic and serotonergic systems in several psychiatric conditions.

Molecular genetics of cocaine use disorders in humans

Drug addiction, one of the major health problems worldwide, is characterized by the loss of control in drug intake, craving, and withdrawal. At the individual level, drugs of abuse produce serious consequences on health and have a negative impact on the family environment and on interpersonal and work relationships. At a wider scale, they have significant socio-economic and public health consequences and they cause delinquency and citizen insecurity. Cocaine, a psychostimulant substance, is one of the most used illicit drugs, especially in America, Western Europe, and Australia. Cocaine use disorders (CUD) are complex multifactorial conditions driven by both genetic and environmental influences. Importantly, not all people who use cocaine develop CUD, and this is due, at least in part, to biological factors that are encoded in the genome of individuals. Acute and repeated use of cocaine induces epigenetic and gene expression changes responsible for the neuronal adaptations and the remodeling of brain circuits that lead to the transition from use to abuse or dependence. The purpose of this review is to delineate such factors, which should eventually help to understand the inter-individual variability in the susceptibility to cocaine addiction. Heritability estimates for CUD are high and genetic risk factors for cocaine addiction have been investigated by candidate gene association studies (CGAS) and genome-wide association studies (GWAS), reviewed here. Also, the high comorbidity that exists between CUD and several other psychiatric disorders is well known and includes phenotypes like schizophrenia, aggression, antisocial or risk-taking behaviors. Such comorbidities are associated with a worse lifetime trajectory, and here we report shared genetic factors that may contribute to them. Gene expression changes and epigenetic modifications induced by cocaine use and chronic abuse in humans are addressed by reviewing transcriptomic studies performed on neuronal cells and on postmortem brains. We report some genes which expression is altered by cocaine that also bear genetic risk variants for the disorder. Finally, we have a glance to the pharmacogenetics of CUD treatments, still in early stages. A better understanding of the genetic underpinnings of CUD will foster the search of effective treatments and help to move forward to personalized medicine.

Behavioural and functional evidence revealing the role of RBFOX1 variation in multiple psychiatric disorders and traits

Common variation in the gene encoding the neuron-specific RNA splicing factor RNA Binding Fox-1 Homolog 1 (RBFOX1) has been identified as a risk factor for several psychiatric conditions, and rare genetic variants have been found causal for autism spectrum disorder (ASD). Here, we explored the genetic landscape of RBFOX1 more deeply, integrating evidence from existing and new human studies as well as studies in Rbfox1 knockout mice. Mining existing data from large-scale studies of human common genetic variants, we confirmed gene-based and genome-wide association of RBFOX1 with risk tolerance, major depressive disorder and schizophrenia. Data on six mental disorders revealed copy number losses and gains to be more frequent in ASD cases than in controls. Consistently, RBFOX1 expression appeared decreased in post-mortem frontal and temporal cortices of individuals with ASD and prefrontal cortex of individuals with schizophrenia. Brain-functional MRI studies demonstrated that carriers of a common RBFOX1 variant, rs6500744, displayed increased neural reactivity to emotional stimuli, reduced prefrontal processing during cognitive control, and enhanced fear expression after fear conditioning, going along with increased avoidance behaviour. Investigating Rbfox1 neuron-specific knockout mice allowed us to further specify the role of this gene in behaviour. The model was characterised by pronounced hyperactivity, stereotyped behaviour, impairments in fear acquisition and extinction, reduced social interest, and lack of aggression; it provides excellent construct and face validity as an animal model of ASD. In conclusion, convergent translational evidence shows that common variants in RBFOX1 are associated with a broad spectrum of psychiatric traits and disorders, while rare genetic variation seems to expose to early-onset neurodevelopmental psychiatric disorders with and without developmental delay like ASD, in particular. Studying the pleiotropic nature of RBFOX1 can profoundly enhance our understanding of mental disorder vulnerability.

Exploring the Contribution to ADHD of Genes Involved in Mendelian Disorders Presenting with Hyperactivity and/or Inattention

Attention-deficit hyperactivity disorder (ADHD) is a complex neurodevelopmental disorder characterized by hyperactivity, impulsivity, and/or inattention, which are symptoms also observed in many rare genetic disorders. We searched for genes involved in Mendelian disorders presenting with ADHD symptoms in the Online Mendelian Inheritance in Man (OMIM) database, to curate a list of new candidate risk genes for ADHD. We explored the enrichment of functions and pathways in this gene list, and tested whether rare or common variants in these genes are associated with ADHD or with its comorbidities. We identified 139 genes, causal for 137 rare disorders, mainly related to neurodevelopmental and brain function. Most of these Mendelian disorders also present with other psychiatric traits that are often comorbid with ADHD. Using whole exome sequencing (WES) data from 668 ADHD cases, we found rare variants associated with the dimension of the severity of inattention symptoms in three genes: KIF11, WAC, and CRBN. Then, we focused on common variants and identified six genes associated with ADHD (in 19,099 cases and 34,194 controls): MANBA, UQCC2, HIVEP2, FOPX1, KANSL1, and AUH. Furthermore, HIVEP2, FOXP1, and KANSL1 were nominally associated with autism spectrum disorder (ASD) (18,382 cases and 27,969 controls), as well as HIVEP2 with anxiety (7016 cases and 14,475 controls), and FOXP1 with aggression (18,988 individuals), which is in line with the symptomatology of the rare disorders they are responsible for. In conclusion, inspecting Mendelian disorders and the genes responsible for them constitutes a valuable approach for identifying new risk genes and the mechanisms of complex disorders.

Reduced cue-induced reinstatement of cocaine-seeking behavior in Plcb1 +/- mice

Cocaine addiction causes serious health problems, and no effective treatment is available yet. We previously identified a genetic risk variant for cocaine addiction in the PLCB1 gene and found this gene upregulated in postmortem brains of cocaine abusers and in human dopaminergic neuron-like cells after an acute cocaine exposure. Here, we functionally tested the contribution of the PLCB1 gene to cocaine addictive properties using Plcb1+/- mice. First, we performed a general phenotypic characterization and found that Plcb1+/- mice showed normal behavior, although they had increased anxiety and impaired short-term memory. Subsequently, mice were trained for operant conditioning, self-administered cocaine for 10 days, and were tested for cocaine motivation. After extinction, we found a reduction in the cue-induced reinstatement of cocaine-seeking behavior in Plcb1+/- mice. After reinstatement, we identified transcriptomic alterations in the medial prefrontal cortex of Plcb1+/- mice, mostly related to pathways relevant to addiction like the dopaminergic synapse and long-term potentiation. To conclude, we found that heterozygous deletion of the Plcb1 gene decreases cue-induced reinstatement of cocaine-seeking, pointing at PLCB1 as a possible therapeutic target for preventing relapse and treating cocaine addiction.

Exploring allele specific methylation in drug dependence susceptibility

Drug dependence is a neuropsychiatric condition that involves genetic, epigenetic and environmental factors. Allele-specific methylation (ASM) is a common and stable epigenetic mechanism that involves genetic variants correlating with differential levels of methylation at CpG sites. We selected 182 single-nucleotide polymorphisms (SNPs) described to influence cis ASM in human brain regions to evaluate their possible contribution to drug dependence susceptibility. We performed a case-control association study in a discovery sample of 578 drug-dependent patients (including 428 cocaine-dependent subjects) and 656 controls from Spain, and then, we followed-up the significant associations in an independent sample of 1119 cases (including 589 cocaine-dependent subjects) and 1092 controls. In the discovery sample, we identified five nominal associations, one of them replicated in the follow-up sample (rs6020251). The pooled analysis revealed an association between drug dependence and rs6020251 but also rs11585570, both overcoming the Bonferroni correction for multiple testing. We performed the same analysis considering only cocaine-dependent patients and obtained similar results. The rs6020251 variant correlates with differential methylation levels of cg17974185 and lies in the first intron of the CTNNBL1 gene, in a genomic region with multiple histone marks related to enhancer and promoter regions in brain. Rs11585570 is an eQTL in brain and blood for the SCP2 and ECHDC2 genes and correlates with differential methylation of cg27535305 and cg13461509, located in the promoter regions of both genes. To conclude, using an approach that combines genetic and epigenetic data, we highlighted the CTNNBL1, SCP2 and ECHDC2 genes as potential contributors to drug dependence susceptibility.

Characterization of an eutherian gene cluster generated after transposon domestication identifies Bex3 as relevant for advanced neurological functions

BACKGROUND

One of the most unusual sources of phylogenetically restricted genes is the molecular domestication of transposable elements into a host genome as functional genes. Although these kinds of events are sometimes at the core of key macroevolutionary changes, their origin and organismal function are generally poorly understood.

RESULTS

Here, we identify several previously unreported transposable element domestication events in the human and mouse genomes. Among them, we find a remarkable molecular domestication that gave rise to a multigenic family in placental mammals, the Bex/Tceal gene cluster. These genes, which act as hub proteins within diverse signaling pathways, have been associated with neurological features of human patients carrying genomic microdeletions in chromosome X. The Bex/Tceal genes display neural-enriched patterns and are differentially expressed in human neurological disorders, such as autism and schizophrenia. Two different murine alleles of the cluster member Bex3 display morphological and physiopathological brain modifications, such as reduced interneuron number and hippocampal electrophysiological imbalance, alterations that translate into distinct behavioral phenotypes.

CONCLUSIONS

We provide an in-depth understanding of the emergence of a gene cluster that originated by transposon domestication and gene duplication at the origin of placental mammals, an evolutionary process that transformed a non-functional transposon sequence into novel components of the eutherian genome. These genes were integrated into existing signaling pathways involved in the development, maintenance, and function of the CNS in eutherians. At least one of its members, Bex3, is relevant for higher brain functions in placental mammals and may be involved in human neurological disorders.

Involvement of the 14-3-3 Gene Family in Autism Spectrum Disorder and Schizophrenia: Genetics, Transcriptomics and Functional Analyses

The 14-3-3 protein family are molecular chaperones involved in several biological functions and neurological diseases. We previously pinpointed YWHAZ (encoding 14-3-3ζ) as a candidate gene for autism spectrum disorder (ASD) through a whole-exome sequencing study, which identified a frameshift variant within the gene (c.659-660insT, p.L220Ffs*18). Here, we explored the contribution of the seven human 14-3-3 family members in ASD and other psychiatric disorders by investigating the: (i) functional impact of the 14-3-3ζ mutation p.L220Ffs*18 by assessing solubility, target binding and dimerization; (ii) contribution of common risk variants in 14-3-3 genes to ASD and additional psychiatric disorders; (iii) burden of rare variants in ASD and schizophrenia; and iv) 14-3-3 gene expression using ASD and schizophrenia transcriptomic data. We found that the mutant 14-3-3ζ protein had decreased solubility and lost its ability to form heterodimers and bind to its target tyrosine hydroxylase. Gene-based analyses using publicly available datasets revealed that common variants in YWHAE contribute to schizophrenia (p = 6.6 × 10-7), whereas ultra-rare variants were found enriched in ASD across the 14-3-3 genes (p = 0.017) and in schizophrenia for YWHAZ (meta-p = 0.017). Furthermore, expression of 14-3-3 genes was altered in post-mortem brains of ASD and schizophrenia patients. Our study supports a role for the 14-3-3 family in ASD and schizophrenia.

Variants of the Aggression-Related RBFOX1 Gene in a Population Representative Birth Cohort Study: Aggressiveness, Personality, and Alcohol Use Disorder

Background: Recently, RBFOX1, a gene encoding an RNA binding protein, has consistently been associated with aggressive and antisocial behavior. Several loci in the gene have been nominally associated with aggression in genome-wide association studies, the risk alleles being more frequent in the general population. We have hence examined the association of four RBFOX1 single nucleotide polymorphisms, previously found related to aggressive traits, with aggressiveness, personality, and alcohol use disorder in birth cohort representative samples. Methods: We used both birth cohorts of the Estonian Children Personality Behavior and Health Study (ECPBHS; original n = 1,238). Aggressiveness was assessed using the Buss-Perry Aggression Questionnaire and the Lifetime History of Aggressiveness structured interview at age 25 (younger cohort) or 33 (older cohort). Big Five personality at age 25 was measured with self-reports and the lifetime occurrence of alcohol use disorder assessed with the MINI interview. RBFOX1 polymorphisms rs809682, rs8062784, rs12921846, and rs6500744 were genotyped in all participants. Given the restricted size of the sample, correction for multiple comparisons was not applied. Results: Aggressiveness was not significantly associated with the RBFOX1 genotype. RBFOX1 rs8062784 was associated with neuroticism and rs809682 with extraversion. Two out of four analyzed RBFOX1 variants, rs8062784 and rs12921846, were associated with the occurrence of alcohol use disorder. Conclusions: In the birth cohort representative sample of the ECPBHS, no association of RBFOX1 with aggressiveness was found, but RBFOX1 variants affected basic personality traits and the prevalence of alcohol use disorder. Future studies on RBFOX1 should consider the moderating role of personality and alcohol use patterns in aggressiveness.

RBFOX1, encoding a splicing regulator, is a candidate gene for aggressive behavior

The RBFOX1 gene (or A2BP1) encodes a splicing factor important for neuronal development that has been related to autism spectrum disorder and other neurodevelopmental phenotypes. Evidence from complementary sources suggests that this gene contributes to aggressive behavior. Suggestive associations with RBFOX1 have been identified in genome-wide association studies (GWAS) of anger, conduct disorder, and aggressive behavior. Nominal association signals in RBFOX1 were also found in an epigenome-wide association study (EWAS) of aggressive behavior. Also, variants in this gene affect temporal lobe volume, a brain area that is altered in several aggression-related phenotypes. In animals, this gene has been shown to modulate aggressive behavior in Drosophila. RBFOX1 has also been associated with canine aggression and is upregulated in mice that show increased aggression after frustration of an expected reward. Associated common genetic variants as well as rare duplications and deletions affecting RBFOX1 have been identified in several psychiatric and neurodevelopmental disorders that are often comorbid with aggressive behaviors. In this paper, we comprehensively review the cumulative evidence linking RBFOX1 to aggression behavior and provide new results implicating RBFOX1 in this phenotype. Most of these studies (genetic and epigenetic analyses in humans, neuroimaging genetics, gene expression and animal models) are hypothesis-free, which strengthens the validity of the findings, although all the evidence is nominal and should therefore be taken with caution. Further studies are required to clarify in detail the role of this gene in this complex phenotype.

Exploring genetic variation that influences brain methylation in attention-deficit/hyperactivity disorder

Attention-deficit/hyperactivity disorder (ADHD) is a neurodevelopmental disorder caused by an interplay of genetic and environmental factors. Epigenetics is crucial to lasting changes in gene expression in the brain. Recent studies suggest a role for DNA methylation in ADHD. We explored the contribution to ADHD of allele-specific methylation (ASM), an epigenetic mechanism that involves SNPs correlating with differential levels of DNA methylation at CpG sites. We selected 3896 tagSNPs reported to influence methylation in human brain regions and performed a case-control association study using the summary statistics from the largest GWAS meta-analysis of ADHD, comprising 20,183 cases and 35,191 controls. We observed that genetic risk variants for ADHD are enriched in ASM SNPs and identified associations with eight tagSNPs that were significant at a 5% false discovery rate (FDR). These SNPs correlated with methylation of CpG sites lying in the promoter regions of six genes. Since methylation may affect gene expression, we inspected these ASM SNPs together with 52 ASM SNPs in high LD with them for eQTLs in brain tissues and observed that the expression of three of those genes was affected by them. ADHD risk alleles correlated with increased expression (and decreased methylation) of ARTN and PIDD1 and with a decreased expression (and increased methylation) of C2orf82. Furthermore, these three genes were predicted to have altered expression in ADHD, and genetic variants in C2orf82 correlated with brain volumes. In summary, we followed a systematic approach to identify risk variants for ADHD that correlated with differential cis-methylation, identifying three novel genes contributing to the disorder.

Truncating variant burden in high-functioning autism and pleiotropic effects of LRP1 across psychiatric phenotypes

BACKGROUND

Previous research has implicated de novo and inherited truncating mutations in autism-spectrum disorder. We aim to investigate whether the load of inherited truncating mutations contributes similarly to high-functioning autism, and to characterize genes that harbour de novo variants in high-functioning autism.

METHODS

We performed whole-exome sequencing in 20 high-functioning autism families (average IQ = 100).

RESULTS

We observed no difference in the number of transmitted versus nontransmitted truncating alleles for high-functioning autism (117 v. 130, p = 0.78). Transmitted truncating and de novo variants in high-functioning autism were not enriched in gene ontology (GO) or Kyoto Encyclopedia of Genes and Genomes (KEGG) categories, or in autism-related gene sets. However, in a patient with high-functioning autism we identified a de novo variant in a canonical splice site of LRP1, a postsynaptic density gene that is a target for fragile X mental retardation protein (FRMP). This de novo variant leads to in-frame skipping of exon 29, removing 2 of 6 blades of the β-propeller domain 4 of LRP1, with putative functional consequences. Large data sets implicate LRP1 across a number of psychiatric disorders: de novo variants are associated with autism-spectrum disorder (p = 0.039) and schizophrenia (p = 0.008) from combined sequencing projects; common variants using genome-wide association study data sets from the Psychiatric Genomics Consortium show gene-based association in schizophrenia (p = 6.6 × E−07) and in a meta-analysis across 7 psychiatric disorders (p = 2.3 × E−03); and the burden of ultra-rare pathogenic variants has been shown to be higher in autism-spectrum disorder (p = 1.2 × E−05), using whole-exome sequencing from 6135 patients with schizophrenia, 1778 patients with autism-spectrum disorder and 7875 controls.

LIMITATIONS

We had a limited sample of patients with high-functioning autism, related to difficulty in recruiting probands with high cognitive performance and no family history of psychiatric disorders.

CONCLUSION

Previous studies and ours suggest an effect of truncating mutations restricted to severe autism-spectrum disorder phenotypes that are associated with intellectual disability. We provide evidence for pleiotropic effects of common and rare variants in the LRP1 gene across psychiatric phenotypes.

Genome-wide association meta-analysis of cocaine dependence: Shared genetics with comorbid conditions

Cocaine dependence is a complex psychiatric disorder that is highly comorbid with other psychiatric traits. Twin and adoption studies suggest that genetic variants contribute substantially to cocaine dependence susceptibility, which has an estimated heritability of 65-79%. Here we performed a meta-analysis of genome-wide association studies of cocaine dependence using four datasets from the dbGaP repository (2085 cases and 4293 controls, all of them selected by their European ancestry). Although no genome-wide significant hits were found in the SNP-based analysis, the gene-based analysis identified HIST1H2BD as associated with cocaine-dependence (10% FDR). This gene is located in a region on chromosome 6 enriched in histone-related genes, previously associated with schizophrenia (SCZ). Furthermore, we performed LD Score regression analysis with comorbid conditions and found significant genetic correlations between cocaine dependence and SCZ, ADHD, major depressive disorder (MDD) and risk taking. We also found, through polygenic risk score analysis, that all tested phenotypes are significantly associated with cocaine dependence status: SCZ (R² = 2.28%; P = 1.21e-26), ADHD (R² = 1.39%; P = 4.5e-17), risk taking (R² = 0.60%; P = 2.7e-08), MDD (R² = 1.21%; P = 4.35e-15), children's aggressive behavior (R² = 0.3%; P = 8.8e-05) and antisocial behavior (R² = 1.33%; P = 2.2e-16). To our knowledge, this is the largest reported cocaine dependence GWAS meta-analysis in European-ancestry individuals. We identified suggestive associations in regions that may be related to cocaine dependence and found evidence for shared genetic risk factors between cocaine dependence and several comorbid psychiatric traits. However, the sample size is limited and further studies are needed to confirm these results.

An integrated analysis of genes and functional pathways for aggression in human and rodent models

Human genome-wide association studies (GWAS), transcriptome analyses of animal models, and candidate gene studies have advanced our understanding of the genetic architecture of aggressive behaviors. However, each of these methods presents unique limitations. To generate a more confident and comprehensive view of the complex genetics underlying aggression, we undertook an integrated, cross-species approach. We focused on human and rodent models to derive eight gene lists from three main categories of genetic evidence: two sets of genes identified in GWAS studies, four sets implicated by transcriptome-wide studies of rodent models, and two sets of genes with causal evidence from online Mendelian inheritance in man (OMIM) and knockout (KO) mice reports. These gene sets were evaluated for overlap and pathway enrichment to extract their similarities and differences. We identified enriched common pathways such as the G-protein coupled receptor (GPCR) signaling pathway, axon guidance, reelin signaling in neurons, and ERK/MAPK signaling. Also, individual genes were ranked based on their cumulative weights to quantify their importance as risk factors for aggressive behavior, which resulted in 40 top-ranked and highly interconnected genes. The results of our cross-species and integrated approach provide insights into the genetic etiology of aggression.

MiR-9, miR-153 and miR-124 are down-regulated by acute exposure to cocaine in a dopaminergic cell model and may contribute to cocaine dependence

Cocaine is one of the most used psychostimulant drugs worldwide. MicroRNAs are post-transcriptional regulators of gene expression that are highly expressed in brain, and several studies have shown that cocaine can alter their expression. In a previous study, we identified several protein-coding genes that are differentially expressed in a dopaminergic neuron-like model after an acute exposure to cocaine. Now, we used the prediction tool WebGestalt to identify miRNA molecules potentially involved in the regulation of these genes. Using the same cellular model, we found that seven of these miRNAs are down-regulated by cocaine: miR-124-3p, miR-124-5p, miR-137, miR-101-3p, miR-9-5p, miR-369-3p and miR-153-3p, the last three not previously related to cocaine. Furthermore, we found that three of the miRNA genes that are differentially expressed in our model (hsa-miR-9-1, hsa-miR-153-1 and hsa-miR-124-3) are nominally associated with cocaine dependence in a case-control study (2,085 cases and 4,293 controls). In summary, we highlighted novel miRNAs that may be involved in those cocaine-induced changes of gene expression that underlie addiction. Moreover, we identified genetic variants that contribute to cocaine dependence in three of these miRNA genes, supporting the idea that genes differentially expressed under cocaine may play an important role in the susceptibility to cocaine dependence.

Evaluation of previous substance dependence genome-wide significant findings in a Spanish sample

BACKGROUND

Substance dependence is a chronic and relapsing disorder explained by genetic and environmental risk factors. The aim of our study is to replicate previous genome-wide significant (GWS) hits identified in substance dependence in general or in cocaine dependence in particular using an independent sample from Spain.

METHODS

We evaluated, in a Spanish sample of 1711 subjects with substance dependence (1011 of them cocaine dependent) and 1719 control individuals, three SNPs identified as GWS in previous studies: rs1868152 and rs2952621 (located near LINC02052 and LINC01854, respectively), associated with substance dependence, and rs2629540 (in the first intron of FAM53B), associated with cocaine dependence.

RESULTS

We replicated the association between rs2952621 and substance dependence under the dominant model (P = 0.020), with the risk allele (T) being the same in our sample and in those two reported previously. We then performed a meta-analysis of the two samples used in the original study that reported the association of rs2952621 with substance dependence (Collaborative Studies on Genetics of Alcoholism (COGA) and Study of Addiction: Genetics and Environment (SAGE)) together with our Spanish sample. The meta-analysis of 3747 cases and 4043 controls confirmed the association (OR = 1.26, 95% CI = 1.15-1.39).

CONCLUSIONS

The rs2952621 variant, located downstream from the yet uncharacterized gene LINC01854, is associated with substance dependence in our Spanish sample. Further research is needed to understand its contribution to the susceptibility to substance dependence.

Association of the PLCB1 gene with drug dependence

Genetic factors involved in the susceptibility to drug addiction still remain largely unknown. MiRNAs seem to play key roles in the drug-induced plasticity of the brain that likely drives the emergence of addiction. In this work we explored the role of miRNAs in drug addiction. With this aim, we selected 62 SNPs located in the 3’UTR of target genes that are predicted to alter the binding of miRNA molecules and performed a case-control association study in a Spanish sample of 735 cases (mainly cocaine-dependent subjects with multiple drug dependencies) and 739 controls. We found an association between rs1047383 in the PLCB1 gene and drug dependence that was replicated in an independent sample (663 cases and 667 controls). Then we selected 9 miRNAs predicted to bind the rs1047383 region, but none of them showed any effect on PLCB1 expression. We also assessed two miRNAs binding a region that contains a SNP in linkage disequilibrium with rs1047383, but although one of them, hsa-miR-582, was found to downregulate PLCB1, no differences were observed between alleles. Finally, we explored the possibility that PLCB1 expression is altered by cocaine and we observed a significant upregulation of the gene in the nucleus accumbens of cocaine abusers and in human dopaminergic-like neurons after cocaine treatment. Our results, together with previous studies, suggest that PLCB1 participates in the susceptibility to drug dependence.

Aggressive behavior in humans: Genes and pathways identified through association studies

Aggressive behavior has both genetic and environmental components. Many association studies have been performed to identify genetic factors underlying aggressive behaviors in humans. In this review we summarize the previous work performed in this field, considering both candidate gene (CGAS) and genome-wide association studies (GWAS), excluding those performed in samples where the primary diagnosis is a psychiatric or neurological disorder other than an aggression-related phenotype. Subsequently, we have studied the enrichment of pathways and functions in GWAS data. The results of our searches show that most CGAS have identified associations with genes involved in dopaminergic and serotonergic neurotransmission and in hormone regulation. On the other hand, GWAS have not yet identified genome-wide significant associations, but top nominal findings are related to several signaling pathways, such as axon guidance or estrogen receptor signaling, and also to neurodevelopmental processes and synaptic plasticity. Future studies should use larger samples, homogeneous phenotypes and standardized measurements to identify genes that underlie aggressive behaviors in humans. © 2016 Wiley Periodicals, Inc.

Genetics of aggressive behavior: An overview

The Research Domain Criteria (RDoC) address three types of aggression: frustrative non-reward, defensive aggression and offensive/proactive aggression. This review sought to present the evidence for genetic underpinnings of aggression and to determine to what degree prior studies have examined phenotypes that fit into the RDoC framework. Although the constructs of defensive and offensive aggression have been widely used in the animal genetics literature, the human literature is mostly agnostic with regard to all the RDoC constructs. We know from twin studies that about half the variance in behavior may be explained by genetic risk factors. This is true for both dimensional, trait-like, measures of aggression and categorical definitions of psychopathology. The non-shared environment seems to have a moderate influence with the effects of shared environment being unclear. Human molecular genetic studies of aggression are in an early stage. The most promising candidates are in the dopaminergic and serotonergic systems along with hormonal regulators. Genome-wide association studies have not yet achieved genome-wide significance, but current samples are too small to detect variants having the small effects one would expect for a complex disorder. The strongest molecular evidence for a genetic basis for aggression comes from animal models comparing aggressive and non-aggressive strains or documenting the effects of gene knockouts. Although we have learned much from these prior studies, future studies should improve the measurement of aggression by using a systematic method of measurement such as that proposed by the RDoC initiative.

Lack of association between the LPR and VNTR polymorphisms of the serotonin transporter gene and cocaine dependence in a Spanish sample

We genotyped the LPR and VNTR polymorphisms of the serotonin transporter gene in 504 cocaine-dependent patients and 508 controls. No association was detected with either polymorphism or with any haplotype combination. This study provided no evidence that these polymorphisms confer susceptibility to cocaine dependence in our sample.

Candidate pathway association study in cocaine dependence: the control of neurotransmitter release

OBJECTIVES

Cocaine is the second most used illegal drug in Europe. The transition from use to dependence involves both genetic and environmental factors. Genetic variation in neurotransmitter systems is involved in the susceptibility to cocaine dependence. We examined the possible contribution to cocaine dependence of 16 genes involved in the cellular machinery that controls neurotransmitter release: genes encoding proteins of the SNARE complex (STX1A, SNAP25, VAMP1 and VAMP2), fusion control elements (SYT1, SYT2, CPLX1, CPLX2, CPLX3 and CPLX4) and regulatory elements (STXBP1, SYP, SNPH, NSF, NAPA and RAB3A).

METHODS

We genotyped 121 SNPs, selected according to genetic coverage criteria, in 360 cocaine-dependent patients and 360 controls from Spain.

RESULTS

Single and multiple-marker analyses revealed a strong association between cocaine dependence and the NSF gene, encoding the N-ethylmaleimide-sensitive factor (P = 5.1e-04, OR = 2.44 (1.45-4.00) and P = 0.001, OR = 1.82 (1.28-2.59), respectively). The presence and absence of psychotic symptoms were also studied. Interestingly, when we considered the time between initial consumption and the onset of cocaine dependence, we observed that the association was mainly restricted to the group of patients that rapidly developed drug dependence (≤ 2 years; P = 2.98e-06, OR = 1.33 (1.20-1.47)).

CONCLUSIONS

Our data show preliminary evidence that NSF may predispose not only to cocaine dependence, but also to an early onset of the dependence.

Active and passive MDMA ('ecstasy') intake induces differential transcriptional changes in the mouse brain

3,4-Methylenedioxymethamphetamine (MDMA, 'ecstasy') is a recreational drug widely used by adolescents and young adults. Although its rewarding effects are well established, there is controversy on its addictive potential. We aimed to compare the consequences of active and passive MDMA administration on gene expression in the mouse brain since all previous studies were based on passive MDMA administration. We used a yoked-control operant intravenous self-administration paradigm combined with microarray technology. Transcriptomic profiles of ventral striatum, frontal cortex, dorsal raphe nucleus and hippocampus were analysed in mice divided in contingent MDMA, yoked MDMA and yoked saline groups, and several changes were validated by quantitative reverse transcription polymerase chain reaction (qRT-PCR). The comparison of contingent MDMA and yoked MDMA vs. yoked saline mice allowed the identification of differential expression in several genes, most of them with immunological and inflammatory functions, but others being involved in neuroadaptation. In the comparison of contingent MDMA vs. yoked MDMA administration, hippocampus and the dorsal raphe nucleus showed statistically significant changes. The altered expression of several genes involved in neuroadaptative changes and synapse function, which may be related to learning self-administration behaviour, could be validated in these two brain structures. In conclusion, our study shows a strong effect of MDMA administration on the expression of immunological and inflammatory genes in all the four brain regions studied. In addition, experiments on MDMA self-administration suggest that the dorsal raphe nucleus and hippocampus may be involved in active MDMA-seeking behaviour, and show specific alterations on gene expression that support the addictive potential of this drug.

Late-onset episodic ataxia type 2 associated with a novel loss-of-function mutation in the CACNA1A gene

We report a patient with typical features of episodic ataxia type 2 (EA2) but with onset in the sixth decade and associated interictal hand dystonia. He was found to bear the novel heterozygous missense mutation p.Gly638Asp (c.1913G>A) in the CACNA1A gene. Functional analysis of the mutation on P/Q channels expressed in HEK 293 cells revealed a reduction of Ca(2+) current densities, a left-shift in the apparent reversal potential, the slowing of inactivation kinetics and the increase in the rate of current recovery from inactivation. These results are consistent with a decrease in Ca(2+) permeability through mutant P/Q channels. To our knowledge, this is just the second patient with late onset EA2 linked to a CACNA1A mutation and the first to carry a loss-of-function missense mutation.

Genetic Analysis of 27 Spanish Patients with Hemiplegic Migraine, Basilar-Type Migraine and Childhood Periodic Syndromes

Familial hemiplegic migraine (FHM) is a rare type of migraine with aura. Mutations in three genes have been described in FHM patients: CACNA1A (FHM1), ATP1A2 (FHM2) and SCN1A (FHM3). We screened 27 Spanish patients with hemiplegic migraine (HM), basilar-type migraine or childhood periodic syndromes (CPS) for mutations in these three genes. Two novel CACNA1A variants, p.Val581Met and p.Tyr1245Cys, and a previously annotated change, p.Cys1534Ser, were identified in individuals with HM, although they have not yet been proven to be pathogenic. Interestingly, p.Tyr1245Cys was detected in a patient displaying a changing, age-specific phenotype that began as benign paroxysmal torticollis of infancy, evolving into benign paroxysmal vertigo of childhood and later becoming HM. This is the first instance of a specific non-synonymous base change being described in a subject affected with CPS. The fact that the molecular screen identified non-synonymous changes in< 15± of our HM patients further stresses the genetic heterogeneity underlying the presumably monogenic forms of migraine.