Complex control of GABA(A) receptor subunit mRNA expression: variation, covariation, and genetic regulation

The de novo missense mutation F224S in GABRB2, identified in epileptic encephalopathy and developmental delay, impairs GABAAR function

Genetic variants in genes encoding subunits of the γ-aminobutyric acid-A receptor (GABAAR) have been found to cause neurodevelopmental disorders and epileptic encephalopathy. In a patient with epilepsy and developmental delay, a de novo heterozygous missense mutation c.671 T > C (p.F224S) was discovered in the GABRB2 gene, which encodes the β2 subunit of GABAAR. Based on previous studies on GABRB2 variants, this new GABRB2 variant (F224S) would be pathogenic. To confirm and investigate the effects of this GABRB2 mutation on GABAAR channel function, we conducted transient expression experiments using GABAAR subunits in HEK293T cells. The GABAARs containing mutant β2 (F224S) subunit showed poor trafficking to the cell membrane, while the expression and distribution of the normal α1 and γ2 subunits were unaffected. Furthermore, the peak current amplitude of the GABAAR containing the β2 (F224S) subunit was significantly smaller compared to the wild type GABAAR. We propose that GABRB2 variant F224S is pathogenic and GABAARs containing this β2 mutant reduce response to GABA under physiological conditions, which could potentially disrupt the excitation/inhibition balance in the brain, leading to epilepsy.

Alzheimer's disease rewires gene coexpression networks coupling different brain regions

Connectome studies have shown how Alzheimer's disease (AD) disrupts functional and structural connectivity among brain regions. But the molecular basis of such disruptions is less studied, with most genomic/transcriptomic studies performing within-brain-region analyses. To inspect how AD rewires the correlation structure among genes in different brain regions, we performed an Inter-brain-region Differential Correlation (Inter-DC) analysis of RNA-seq data from Mount Sinai Brain Bank on four brain regions (frontal pole, superior temporal gyrus, parahippocampal gyrus and inferior frontal gyrus, comprising 264 AD and 372 control human post-mortem samples). An Inter-DC network was assembled from all pairs of genes across two brain regions that gained (or lost) correlation strength in the AD group relative to controls at FDR 1%. The differentially correlated (DC) genes in this network complemented known differentially expressed genes in AD, and likely reflects cell-intrinsic changes since we adjusted for cell compositional effects. Each brain region used a distinctive set of DC genes when coupling with other regions, with parahippocampal gyrus showing the most rewiring, consistent with its known vulnerability to AD. The Inter-DC network revealed master dysregulation hubs in AD (at genes ZKSCAN1, SLC5A3, RCC1, IL17RB, PLK4, etc.), inter-region gene modules enriched for known AD pathways (synaptic signaling, endocytosis, etc.), and candidate signaling molecules that could mediate region-region communication. The Inter-DC network generated in this study is a valuable resource of gene pairs, pathways and signaling molecules whose inter-brain-region functional coupling is disrupted in AD, thereby offering a new perspective of AD etiology.

Spatial and single-nucleus transcriptomics decoding the molecular landscape and cellular organization of avian optic tectum

The avian optic tectum (OT) has been studied for its diverse functions, yet a comprehensive molecular landscape at the cellular level has been lacking. In this study, we applied spatial transcriptome sequencing and single-nucleus RNA sequencing (snRNA-seq) to explore the cellular organization and molecular characteristics of the avian OT from two species: Columba livia and Taeniopygia guttata. We identified precise layer structures and provided comprehensive layer-specific signatures of avian OT. Furthermore, we elucidated diverse functions in different layers, with the stratum griseum periventriculare (SGP) potentially playing a key role in advanced functions of OT, like fear response and associative learning. We characterized detailed neuronal subtypes and identified a population of FOXG1+ excitatory neurons, resembling those found in the mouse neocortex, potentially involved in neocortex-related functions and expansion of avian OT. These findings could contribute to our understanding of the architecture of OT, shedding light on visual perception and multifunctional association.

Methamphetamine-induced region-specific transcriptomic and epigenetic changes in the brain of male rats

Psychostimulant methamphetamine (METH) is neurotoxic to the brain and, therefore, its misuse leads to neurological and psychiatric disorders. The gene regulatory network (GRN) response to neurotoxic METH binge remains unclear in most brain regions. Here we examined the effects of binge METH on the GRN in the nucleus accumbens, dentate gyrus, Ammon's horn, and subventricular zone in male rats. At 24 h after METH, ~16% of genes displayed altered expression and over a quarter of previously open chromatin regions - parts of the genome where genes are typically active - showed shifts in their accessibility. Intriguingly, most changes were unique to each area studied, and independent regulation between transcriptome and chromatin accessibility was observed. Unexpectedly, METH differentially impacted gene activity and chromatin accessibility within the dentate gyrus and Ammon's horn. Around 70% of the affected chromatin-accessible regions in the rat brain have conserved DNA sequences in the human genome. These regions frequently act as enhancers, ramping up the activity of nearby genes, and contain mutations linked to various neurological conditions. By sketching out the gene regulatory networks associated with binge METH in specific brain regions, our study offers fresh insights into how METH can trigger profound, region-specific molecular shifts.

Association between the SLC6A11 rs2304725 and GABRG2 rs211037 polymorphisms and drug-resistant epilepsy: a meta-analysis

Background: Previous studies have shown that SLC6A11 and GABRG2 are linked to drug-resistant epilepsy (DRE), although there have been conflicting results in the literature. In this study, we systematically assessed the relationship between DRE and these two genes. Methods: We systematically searched the PubMed, Embase, Cochrane Library, Web of Science, Google Scholar, Wanfang Data, CNKI, and VIP databases. To clarify whether heterogeneity existed between studies, tools such as the Q-test and I ² statistic were selected. According to study heterogeneity, we chose fixed- or random-effects models for analysis. We then used the chi-squared ratio to evaluate any bias of the experimental data. Results: In total, 11 trials and 3,813 patients were selected. To investigate the relationship with DRE, we performed model tests on the two genes separately. The results showed that SLC6A11 rs2304725 had no significant correlation with DRE risk in the allele, dominant, recessive, and additive models in a pooled population. However, for the over-dominant model, DRE was correlated with rs2304725 (OR = 1.08, 95% CI: 0.92-1.27, p = 0.33) in a pooled population. Similarly, rs211037 was weakly significantly correlated with DRE for the dominant, recessive, over-dominant, and additive models in a pooled population. The subgroup analysis results showed that rs211037 expressed a genetic risk of DRE in allele (OR = 1.01, 95% CI: 0.76-1.35, p = 0.94), dominant (OR = 1.08, 95% CI: 0.77-1.50, p = 0.65), and additive models (OR = 1.14, 95% CI: 0.62-2.09, p = 0.67) in an Asian population. Conclusion: In this meta-analysis, our results showed that SLC6A11 rs2304725 and GABRG2 rs211037 are not significantly correlated with DRE. However, in the over-dominant model, rs2304725 was significantly correlated with DRE. Likewise, rs211037 conveyed a genetic risk for DRE in an Asian population in the allele, dominant, and additive models.

The likelihood approach for potential role of "GABRG2 (C588T, C315T) gene polymorphisms" on the poor response to carbamazepine therapy in Pakhtun population of Pakistan

BACKGROUND

Gamma-aminobutyric acid A receptor, gamma 2 gene (GABRG2) encode the GABAA receptor which is responsible for fast neuronal inhibition. Polymorphisms in GABGR2 gene affect the clinical response of anti-epileptic drugs (AEDs). Therefore, we carried out an updated study to find the association GABRG2 gene polymorphisms with carbamazepine (CBZ) non-responsive therapy in the Pakhtun population.

METHODS

A clinical prospective cohort study was conducted in 79 CBZ treated patients upon consent after the approval of Khyber Medical University Advanced Study and Research Board. Blood sample were taken at optimal dose of CBZ at base line, third and sixth months of the treatment. Blood level of CBZ was measure through reverse phase high performance liquid chromatography (HPLC). Restriction fragment length polymorphisms techniques were used to genotype GABRG2 gene in these patients. CBZ responses were evaluated on three and six months of study by measuring the decrease in frequency of seizure per week.

RESULTS

The average maximum dose of CBZ was 455 ± 133 mg/day at baseline, 479 ± 142 mg/day at third month and 495 ± 133 mg/day at sixth month of the treatment. CBZ level was found within therapeutic range (4-12 mg/L) without any significant (P > .5) variations among the CC, CT and TT genotypes of GABRG2 (C588T and C315T) gene. But the poor clinical response during CBZ treatment was linked (P < .05) with CT and TT genotypes of GABRG2 (C588T and C315T) gene in Pakhtun Population.

CONCLUSION

A poor response to CBZ was found in variant genotypes (CT and TT) of GABRG2 (C588T and C315T) gene in Pakhtun Population.

Regulation of protein abundance in genetically diverse mouse populations

Genetically diverse mouse populations are powerful tools for characterizing the regulation of the proteome and its relationship to whole-organism phenotypes. We used mass spectrometry to profile and quantify the abundance of 6,798 proteins in liver tissue from mice of both sexes across 58 Collaborative Cross (CC) inbred strains. We previously collected liver proteomics data from the related Diversity Outbred (DO) mice and their founder strains. We show concordance across the proteomics datasets despite being generated from separate experiments, allowing comparative analysis. We map protein abundance quantitative trait loci (pQTLs), identifying 1,087 local and 285 distal in the CC mice and 1,706 local and 414 distal in the DO mice. We find that regulatory effects on individual proteins are conserved across the mouse populations, in particular for local genetic variation and sex differences. In comparison, proteins that form complexes are often co-regulated, displaying varying genetic architectures, and overall show lower heritability and map fewer pQTLs. We have made this resource publicly available to enable quantitative analyses of the regulation of the proteome.

Gabra2 is a genetic modifier of Dravet syndrome in mice

Pathogenic variants in epilepsy genes result in a spectrum of clinical severity. One source of phenotypic heterogeneity is modifier genes that affect expressivity of a primary pathogenic variant. Mouse epilepsy models also display varying degrees of clinical severity on different genetic backgrounds. Mice with heterozygous deletion of Scn1a (Scn1a^+/−) model Dravet syndrome, a severe epilepsy most often caused by SCN1A haploinsufficiency. Scn1a^+/− mice recapitulate features of Dravet syndrome, including spontaneous seizures, sudden death, and cognitive/behavioral deficits. Scn1a^+/− mice maintained on the 129S6/SvEvTac (129) strain have normal lifespan and no spontaneous seizures. In contrast, admixture with C57BL/6J (B6) results in epilepsy and premature lethality. We previously mapped Dravet Survival Modifier loci (Dsm1-Dsm5) responsible for strain-dependent differences in survival. Gabra2, encoding the GABA_A α2 subunit, was nominated as a candidate modifier at Dsm1. Direct measurement of GABA_A receptors found lower abundance of α2-containing receptors in hippocampal synapses of B6 mice relative to 129. We also identified a B6-specific single nucleotide deletion within Gabra2 that lowers mRNA and protein by nearly 50%. Repair of this deletion reestablished normal levels of Gabra2 expression. In this study, we used B6 mice with a repaired Gabra2 allele to evaluate Gabra2 as a genetic modifier of severity in Scn1a^+/− mice. Gabra2 repair restored transcript and protein expression, increased abundance of α2-containing GABA_A receptors in hippocampal synapses, and rescued epilepsy phenotypes of Scn1a^+/− mice. These findings validate Gabra2 as a genetic modifier of Dravet syndrome, and support enhancing function of α₂-containing GABA_A receptors as treatment strategy for Dravet syndrome.

Exome sequencing reveals predominantly de novo variants in disorders with intellectual disability (ID) in the founder population of Finland

The genetics of autosomal recessive intellectual disability (ARID) has mainly been studied in consanguineous families, however, founder populations may also be of interest to study intellectual disability (ID) and the contribution of ARID. Here, we used a genotype-driven approach to study the genetic landscape of ID in the founder population of Finland. A total of 39 families with syndromic and non-syndromic ID were analyzed using exome sequencing, which revealed a variant in a known ID gene in 27 families. Notably, 75% of these variants in known ID genes were de novo or suspected de novo (64% autosomal dominant; 11% X-linked) and 25% were inherited (14% autosomal recessive; 7% X-linked; and 4% autosomal dominant). A dual molecular diagnosis was suggested in two families (5%). Via additional analysis and molecular testing, we identified three cases with an abnormal molecular karyotype, including chr21q22.12q22.2 uniparental disomy with a mosaic interstitial 2.7 Mb deletion covering DYRK1A and KCNJ6. Overall, a pathogenic or likely pathogenic variant was identified in 64% (25/39) of the families. Last, we report an alternate inheritance model for 3 known ID genes (UBA7, DDX47, DHX58) and discuss potential candidate genes for ID, including SYPL1 and ERGIC3 with homozygous founder variants and de novo variants in POLR2F and DNAH3. In summary, similar to other European populations, de novo variants were the most common variants underlying ID in the studied Finnish population, with limited contribution of ARID to ID etiology, though mainly driven by founder and potential founder variation in the latter case.

Association between GABRG2 rs211037 polymorphism and febrile seizures: a meta-analysis

Background

Emerging evidence has implied that the GABRG2 gene play a role in the mechanism of febrile seizure (FS), however, the relationship between GABRG2 rs211037 polymorphism and the risk of FS remains controversial. This meta-analysis was conducted to investigate the relationship of GABRG2 rs211037 polymorphism with the susceptibility to FS.

Methods

MEDLINE, Embase, Cochrane Library and CNKI databases were searched (until April 6, 2019) for eligible studies on the relationship between GABRG2 rs211037 polymorphism and FS. We calculated the odds ratios (ORs) by a fixed or random model with the STATA 15.0 software. Subgroup analyses for the ethnicity, the source of the control, and age and sex matching of controls were conducted.

Results

A total of 8 studies consisting of 775 FS patients and 5162 controls were included in this study. Based on the overall data, he GABRG2 rs211037 polymorphism was not significantly associated with the risk of FS (TT + CT vs CC: OR = 0.95, 95%CI 0.64–1.41, P = 0.80). Notably, the GABRG2 rs211037 variant was significantly associated with decreased risk of FS in Asian populations (TT vs CT + CC: OR = 0.63, 95%CI 0.45–0.88, P = 0.006), but increased risk in Caucasian populations (CT vs CC: OR = 1.56, 95%CI 1.14–2.15, P = 0.006). Significant associations were also detected when healthy controls out of the whole controls were employed for comparison (TT vs CT + CC: OR = 0.59, 95% CI 0.45–0.77, P < 0.001) and when data from studies with age- and sex-matched controls were used (TT + CT vs CC: OR = 0.60, 95% CI 0.43–0.86, P = 0.001).

Conclusion

The GABRG2 rs211037 polymorphism may decrease the risk of FS in Asian populations, while increasing the risk in Caucasian populations. Further well-designed studies with large sample sizes are essential to verify the conclusions in other ethnicities.

A high-resolution in vivo atlas of the human brain's benzodiazepine binding site of GABAA receptors

Gamma-aminobutyric acid (GABA) is the main inhibitory neurotransmitter in the human brain and plays a key role in several brain functions and neuropsychiatric disorders such as anxiety, epilepsy, and depression. For decades, several in vivo and ex vivo techniques have been used to highlight the mechanisms of the GABA system, however, no studies have currently combined the techniques to create a high-resolution multimodal view of the GABA system.

Here, we present a quantitative high-resolution in vivo atlas of the human brain benzodiazepine receptor sites (BZR) located on postsynaptic ionotropic GABA_A receptors (GABA_A Rs), generated on the basis of in vivo [¹¹C]flumazenil Positron Emission Tomography (PET) data. Next, based on ex vivo autoradiography data, we transform the PET-generated atlas from binding values into BZR protein density. Finally, we examine the brain regional association between BZR protein density and ex vivo mRNA expression for the 19 subunits in the GABA_A R, including an estimation of the minimally required expression of mRNA levels for each subunit to translate into BZR protein.

This represents the first publicly available quantitative high-resolution in vivo atlas of the spatial distribution of BZR densities in the healthy human brain. The atlas provides a unique neuroscientific tool as well as novel insights into the association between mRNA expression for individual subunits in the GABA_A R and the BZR density at each location in the brain.

Highlights from the Era of Open Source Web-Based Tools

High digital connectivity and a focus on reproducibility are contributing to an open science revolution in neuroscience. Repositories and platforms have emerged across the whole spectrum of subdisciplines, paving the way for a paradigm shift in the way we share, analyze, and reuse vast amounts of data collected across many laboratories. Here, we describe how open access web-based tools are changing the landscape and culture of neuroscience, highlighting six free resources that span subdisciplines from behavior to whole-brain mapping, circuits, neurons, and gene variants.

Content and Performance of the MiniMUGA Genotyping Array: A New Tool To Improve Rigor and Reproducibility in Mouse Research

The laboratory mouse is the most widely used animal model for biomedical research, due in part to its well-annotated genome, wealth of genetic resources, and the ability to precisely manipulate its genome. Despite the importance of genetics for mouse research, genetic quality control (QC) is not standardized, in part due to the lack of cost-effective, informative, and robust platforms. Genotyping arrays are standard tools for mouse research and remain an attractive alternative even in the era of high-throughput whole-genome sequencing. Here, we describe the content and performance of a new iteration of the Mouse Universal Genotyping Array (MUGA), MiniMUGA, an array-based genetic QC platform with over 11,000 probes. In addition to robust discrimination between most classical and wild-derived laboratory strains, MiniMUGA was designed to contain features not available in other platforms: (1) chromosomal sex determination, (2) discrimination between substrains from multiple commercial vendors, (3) diagnostic SNPs for popular laboratory strains, (4) detection of constructs used in genetically engineered mice, and (5) an easy-to-interpret report summarizing these results. In-depth annotation of all probes should facilitate custom analyses by individual researchers. To determine the performance of MiniMUGA, we genotyped 6899 samples from a wide variety of genetic backgrounds. The performance of MiniMUGA compares favorably with three previous iterations of the MUGA family of arrays, both in discrimination capabilities and robustness. We have generated publicly available consensus genotypes for 241 inbred strains including classical, wild-derived, and recombinant inbred lines. Here, we also report the detection of a substantial number of XO and XXY individuals across a variety of sample types, new markers that expand the utility of reduced complexity crosses to genetic backgrounds other than C57BL/6, and the robust detection of 17 genetic constructs. We provide preliminary evidence that the array can be used to identify both partial sex chromosome duplication and mosaicism, and that diagnostic SNPs can be used to determine how long inbred mice have been bred independently from the relevant main stock. We conclude that MiniMUGA is a valuable platform for genetic QC, and an important new tool to increase the rigor and reproducibility of mouse research.

Minor C allele of the SNP rs7873784 associated with rheumatoid arthritis and type-2 diabetes mellitus binds PU.1 and enhances TLR4 expression

Toll-like receptor 4 (TLR4) is an innate immunity receptor predominantly expressed on myeloid cells and involved in the development of various diseases, many of them with complex genetics. Here we present data on functionality of single nucleotide polymorphism rs7873784 located in the 3'-untranslated region (3'-UTR) of TLR4 gene and associated with various pathologies involving chronic inflammation. We demonstrate that TLR4 3'-UTR strongly enhanced the activity of TLR4 promoter in U937 human monocytic cell line while minor rs7873784(C) allele created a binding site for transcription factor PU.1 (encoded by SPI1 gene), a known regulator of TLR4 expression. Increased binding of PU.1 further augmented the TLR4 transcription while PU.1 knockdown or complete disruption of the PU.1 binding site abrogated the effect. We hypothesize that additional functional PU.1 site may increase TLR4 expression in individuals carrying minor C variant of rs7873784 and modulate the development of certain pathologies, such as rheumatoid arthritis and type-2 diabetes mellitus.

Potentiating α2 subunit containing perisomatic GABAA receptors protects against seizures in a mouse model of Dravet syndrome

KEY POINTS

Dravet syndrome mice (Scn1a^+/- ) demonstrate a marked strain dependence for the severity of seizures which is correlated with GABA_A receptor α₂ subunit expression. The α₂ /α₃ subunit selective positive allosteric modulator (PAM) AZD7325 potentiates inhibitory postsynaptic currents (IPSCs) specifically in perisomatic synapses. AZD7325 demonstrates stronger effects on IPSCs in the seizure resistant mouse strain, consistent with higher α₂ subunit expression. AZD7325 demonstrates seizure protective effects in Scn1a^+/- mice without apparent sedative effects in vivo.

ABSTRACT

GABA_A receptor potentiators are commonly used for the treatment of epilepsy, but it is not clear whether targeting distinct GABA_A receptor subtypes will have disproportionate benefits over adverse effects. Here we demonstrate that the α₂ /α₃ selective positive allosteric modulator (PAM) AZD7325 preferentially potentiates hippocampal inhibitory responses at synapses proximal to the soma of CA1 neurons. The effect of AZD7325 on synaptic responses was more prominent in mice on the 129S6/SvEvTac background strain, which have been demonstrated to be seizure resistant in the model of Dravet syndrome (Scn1a^+/- ), and in which the α₂ GABA_A receptor subunits are expressed at higher levels relative to in the seizure prone C57BL/6J background strain. Consistent with this, treatment of Scn1a^+/- mice with AZD7325 elevated the temperature threshold for hyperthermia-induced seizures without apparent sedative effects. Our results in a model system indicate that selectively targeting α₂ is a potential therapeutic option for Dravet syndrome.

Human brain transcriptome analysis finds region- and subject-specific expression signatures of GABAAR subunits

Altered expression of GABA receptors (GABAARs) has been implicated in neurological and psychiatric disorders, but limited information about region-specific GABAAR subunit expression in healthy human brains, heteromeric assembly of major isoforms, and their collective organization across healthy individuals, are major roadblocks to understanding their role in non-physiological states. Here, by using microarray and RNA-Seq datasets-from single cell nuclei to global brain expression-from the Allen Institute, we find that transcriptional expression of GABAAR subunits is anatomically organized according to their neurodevelopmental origin. The data show a combination of complementary and mutually-exclusive expression patterns that delineate major isoforms, and which is highly stereotypical across brains from control donors. We summarize the region-specific signature of GABAR subunits per subject and its variability in a control population sample that can be used as a reference for remodeling changes during homeostatic rearrangements of GABAAR subunits after physiological, pharmacological or pathological challenges.

Identification of a Functional Non-coding Variant in the GABA A Receptor α2 Subunit of the C57BL/6J Mouse Reference Genome: Major Implications for Neuroscience Research

GABA type-A (GABA-A) receptors containing the α2 subunit (GABRA2) are expressed in most brain regions and are critical in modulating inhibitory synaptic function. Genetic variation at the GABRA2 locus has been implicated in epilepsy, affective and psychiatric disorders, alcoholism and drug abuse. Gabra2 expression varies as a function of genotype and is modulated by sequence variants in several brain structures and populations, including F2 crosses originating from C57BL/6J (B6J) and the BXD recombinant inbred family derived from B6J and DBA/2J. Here we demonstrate a global reduction of GABRA2 brain protein and mRNA in the B6J strain relative to other inbred strains, and identify and validate the causal mutation in B6J. The mutation is a single base pair deletion located in an intron adjacent to a splice acceptor site that only occurs in the B6J reference genome. The deletion became fixed in B6J between 1976 and 1991 and is now pervasive in many engineered lines, BXD strains generated after 1991, the Collaborative Cross, and the majority of consomic lines. Repair of the deletion using CRISPR-Cas9-mediated gene editing on a B6J genetic background completely restored brain levels of GABRA2 protein and mRNA. Comparison of transcript expression in hippocampus, cortex, and striatum between B6J and repaired genotypes revealed alterations in GABA-A receptor subunit expression, especially in striatum. These results suggest that naturally occurring variation in GABRA2 levels between B6J and other substrains or inbred strains may also explain strain differences in anxiety-like or alcohol and drug response traits related to striatal function. Characterization of the B6J private mutation in the Gabra2 gene is of critical importance to molecular genetic studies in neurobiological research because this strain is widely used to generate genetically engineered mice and murine genetic populations, and is the most widely utilized strain for evaluation of anxiety-like, depression-like, pain, epilepsy, and drug response traits that may be partly modulated by GABRA2 function.

A computational strategy for finding novel targets and therapeutic compounds for opioid dependence

Opioids are widely used for treating different types of pains, but overuse and abuse of prescription opioids have led to opioid epidemic in the United States. Besides analgesic effects, chronic use of opioid can also cause tolerance, dependence, and even addiction. Effective treatment of opioid addiction remains a big challenge today. Studies on addictive effects of opioids focus on striatum, a main component in the brain responsible for drug dependence and addiction. Some transcription regulators have been associated with opioid addiction, but relationship between analgesic effects of opioids and dependence behaviors mediated by them at the molecular level has not been thoroughly investigated. In this paper, we developed a new computational strategy that identifies novel targets and potential therapeutic molecular compounds for opioid dependence and addiction. We employed several statistical and machine learning techniques and identified differentially expressed genes over time which were associated with dependence-related behaviors after exposure to either morphine or heroin, as well as potential transcription regulators that regulate these genes, using time course gene expression data from mouse striatum. Moreover, our findings revealed that some of these dependence-associated genes and transcription regulators are known to play key roles in opioid-mediated analgesia and tolerance, suggesting that an intricate relationship between opioid-induce pain-related pathways and dependence may develop at an early stage during opioid exposure. Finally, we determined small compounds that can potentially target the dependence-associated genes and transcription regulators. These compounds may facilitate development of effective therapy for opioid dependence and addiction. We also built a database (http://daportals.org) for all opioid-induced dependence-associated genes and transcription regulators that we discovered, as well as the small compounds that target those genes and transcription regulators.

Cross-species molecular dissection across alcohol behavioral domains

This review summarizes the proceedings of a symposium presented at the "Alcoholism and Stress: A Framework for Future Treatment Strategies" conference held in Volterra, Italy on May 9-12, 2017. Psychiatric diseases, including alcohol-use disorders (AUDs), are influenced through complex interactions of genes, neurobiological pathways, and environmental influences. A better understanding of the common neurobiological mechanisms underlying an AUD necessitates an integrative approach, involving a systematic assessment of diverse species and phenotype measures. As part of the World Congress on Stress and Alcoholism, this symposium provided a detailed account of current strategies to identify mechanisms underlying the development and progression of AUDs. Dr. Sean Farris discussed the integration and organization of transcriptome and postmortem human brain data to identify brain regional- and cell type-specific differences related to excessive alcohol consumption that are conserved across species. Dr. Brien Riley presented the results of a genome-wide association study of DSM-IV alcohol dependence; although replication of genetic associations with alcohol phenotypes in humans remains challenging, model organism studies show that COL6A3, KLF12, and RYR3 affect behavioral responses to ethanol, and provide substantial evidence for their role in human alcohol-related traits. Dr. Rob Williams expanded upon the systematic characterization of extensive genetic-genomic resources for quantifying and clarifying phenotypes across species that are relevant to precision medicine in human disease. The symposium concluded with Dr. Robert Hitzemann's description of transcriptome studies in a mouse model selectively bred for high alcohol ("binge-like") consumption and a non-human primate model of long-term alcohol consumption. Together, the different components of this session provided an overview of systems-based approaches that are pioneering the experimental prioritization and validation of novel genes and gene networks linked with a range of behavioral phenotypes associated with stress and AUDs.

Genome-wide association for testis weight in the diversity outbred mouse population

Testis weight is a genetically mediated trait associated with reproductive efficiency across numerous species. We sought to evaluate the genetically diverse, highly recombinant Diversity Outbred (DO) mouse population as a tool to identify and map quantitative trait loci (QTLs) associated with testis weight. Testis weights were recorded for 502 male DO mice and the mice were genotyped on the GIGAMuga array at ~ 143,000 SNPs. We performed a genome-wide association analysis and identified one significant and two suggestive QTLs associated with testis weight. Using bioinformatic approaches, we developed a list of candidate genes and identified those with known roles in testicular size and development. Candidates of particular interest include the RNA demethylase gene Alkbh5, the cyclin-dependent kinase inhibitor gene Cdkn2c, the dynein axonemal heavy chain gene Dnah11, the phospholipase D gene Pld6, the trans-acting transcription factor gene Sp4, and the spermatogenesis-associated gene Spata6, each of which has a human ortholog. Our results demonstrate the utility of DO mice in high-resolution genetic mapping of complex traits, enabling us to identify developmentally important genes in adult mice. Understanding how genetic variation in these genes influence testis weight could aid in the understanding of mechanisms of mammalian reproductive function.

GABRG2 C588T gene polymorphisms might be a predictive genetic marker of febrile seizures and generalized recurrent seizures: a case-control study in a Romanian pediatric population

INTRODUCTION

This case-control study aimed to assess two single nucleotide polymorphisms of the gene encoding the GABRG2 protein - GABRG2 (3145 G>A) and GABRG2 rs 211037 Asn196Asn (C588T) - in a cohort of pediatric patients from Romania, and evaluate their possible impact on drug-resistant forms of generalized epilepsy and recurrent febrile seizures.

MATERIAL AND METHODS

One hundred and fourteen children with idiopathic generalized epilepsy (group 1) or febrile seizures (group 2) were compared to 153 controls. Peripheral blood samples were assessed using polymerase chain reaction-restriction fragment length polymorphism analysis, with results interpreted based on the disappearance of a restriction site in the C allele (122 bp) compared to the T allele (100 bp + 22 bp).

RESULTS

A significant association was found with the TT homozygous genotype and T allele for both febrile seizures and epilepsy for the C588T locus, while GABRG2 G>A 3145 showed no significant association with any type of seizure. The TT homozygous genotype of GABRG2 Asn196Asn polymorphism was more frequent in patients with a history of febrile seizures (p = 0.0001), without a significant association identified for GABRG2-G>A 3145. Composite analysis showed associations with epilepsy for CC-AG (p = 0.02) and CT-AG (p = 0.007) with the CC-AA combination as reference.

CONCLUSIONS

C588T polymorphism of the GABRG2 gene might be a predictive genetic marker in triggering febrile convulsions. GABRG2 rs211037 TT homozygotes and T allele variants have an increased risk for developing febrile seizures. Recurrent crises and repeated episodes of seizures are more frequent in the GABRG2 Asn196Asn TT genotype polymorphism, with a 45 and 8 times higher risk of developing idiopathic generalized epilepsy and recurrent febrile seizures, respectively.

ZKSCAN1 gene and its related circular RNA (circZKSCAN1) both inhibit hepatocellular carcinoma cell growth, migration, and invasion but through different signaling pathways

There is increasing evidence that circular RNA (circRNA) are involved in cancer development, but the regulation and function of human circRNA remain largely unknown. In this study, we demonstrated that ZKSCAN1, a zinc finger family gene, is expressed in both linear and circular (circZKSCAN1) forms of RNA in human hepatocellular carcinoma (HCC) tissues and cell lines. Here, we analyzed a cohort of 102 patients and found that expression of both ZKSCAN1mRNA and circZKSCAN1 was significantly lower (P < 0.05) in the HCC samples compared with that in matched adjacent nontumorous tissues by reverse transcription PCR (RT-PCR). The low expression level of ZKSCAN1 was only associated with tumor size (P = 0.032), while the cirZKSCAN1 levels varied in patients with different tumor numbers (P < 0.01), cirrhosis (P = 0.031), vascular invasion (P = 0.002), or microscopic vascular invasion (P = 0.002), as well as with the tumor grade (P < 0.001). Silencing both ZKSCAN1mRNA and circZKSCAN1 promoted cell proliferation, migration, and invasion. In contrast, overexpression of both forms of RNA repressed HCC progression in vivo and in vitro. Silencing or overexpression of both forms of RNA did not interfere with each other. RNA-seq revealed a very different molecular basis for the observed effects; ZKSCAN1mRNA mainly regulated cellular metabolism, while circZKSCAN1 mediated several cancer-related signaling pathways, suggesting a nonredundant role for ZKSCAN1mRNA and circRNA. In conclusion, our results revealed two post-translational products (ZKSCAN1mRNA and circZKSCAN1) that cooperated closely with one another to inhibit growth, migration, and invasion of HCC. cirZKSCAN1 might be a useful marker for the diagnosis of HCC.

Analyses of differentially expressed genes after exposure to acute stress, acute ethanol, or a combination of both in mice

Alcohol abuse is a complex disorder, which is confounded by other factors, including stress. In the present study, we examined gene expression in the hippocampus of BXD recombinant inbred mice after exposure to ethanol (NOE), stress (RSS), and the combination of both (RSE). Mice were given an intraperitoneal (i.p.) injection of 1.8 g/kg ethanol or saline, and subsets of both groups were exposed to acute restraint stress for 15 min or controls. Gene expression in the hippocampus was examined using microarray analysis. Genes that were significantly (p < 0.05, q < 0.1) differentially expressed were further evaluated. Bioinformatic analyses were predominantly performed using tools available at GeneNetwork.org, and included gene ontology, presence of cis-regulation or polymorphisms, phenotype correlations, and principal component analyses. Comparisons of differential gene expression between groups showed little overlap. Gene Ontology demonstrated distinct biological processes in each group with the combined exposure (RSE) being unique from either the ethanol (NOE) or stress (RSS) group, suggesting that the interaction between these variables is mediated through diverse molecular pathways. This supports the hypothesis that exposure to stress alters ethanol-induced gene expression changes and that exposure to alcohol alters stress-induced gene expression changes. Behavior was profiled in all groups following treatment, and many of the differentially expressed genes are correlated with behavioral variation within experimental groups. Interestingly, in each group several genes were correlated with the same phenotype, suggesting that these genes are the potential origins of significant genetic networks. The distinct sets of differentially expressed genes within each group provide the basis for identifying molecular networks that may aid in understanding the complex interactions between stress and ethanol, and potentially provide relevant therapeutic targets. Using Ptp4a1, a candidate gene underlying the quantitative trait locus for several of these phenotypes, and network analyses, we show that a large group of differentially expressed genes in the NOE group are highly interrelated, some of which have previously been linked to alcohol addiction or alcohol-related phenotypes.

Systems Genetics Analysis to Identify the Genetic Modulation of a Glaucoma-Associated Gene

Loss of retinal ganglion cells (RGCs) is one of the hallmarks of retinal neurodegenerative diseases, glaucoma being one of the most common. Recently, γ-synuclein (SNCG) was shown to be highly expressed in the somas and axons of RGCs. In various mouse models of glaucoma, downregulation of Sncg gene expression correlates with RGC loss. To investigate the regulation of Sncg in RGCs, we used a systems genetics approach to identify a gene that modulates the expression of Sncg, followed by confirmatory studies in both healthy and diseased retinas. We found that chromosome 1 harbors an eQTL that modulates the expression of Sncg in the mouse retina and identified Pfdn2 as the candidate upstream modulator of Sncg expression. Downregulation of Pfdn2 in enriched RGCs causes a concomitant reduction in Sncg. In this chapter, we describe our strategy and methods for identifying and confirming a genetic modulation of a glaucoma-associated gene. A similar method can be applied to other genes expressed in other tissues.

Genetic and immunohistochemical analysis of HSPA5 in mouse and human retinas

PURPOSE

Photoreceptor degenerative diseases are among the leading causes of vision loss. Although the causative genetic mutations are often known, mechanisms leading to photoreceptor degeneration remain poorly defined. We have previously demonstrated that the photoreceptor membrane-associated protein XAP-1 antigen is a product of the HSPA5 gene. In this study, we used systems genetic methods, statistical modeling, and immunostaining to identify and analyze candidate genes that modulate Hspa5 expression in the retina.

METHODS

Quantitative trait locus (QTL) mapping was used to map the genomic region that regulates Hspa5 in the cross between C57BL/6J X DBA/2J mice (BXD) genetic reference panel. The stepwise refinement of candidate genes was based on expression QTL mapping, gene expression correlation analyses (direct and partial), and analysis of regional sequence variants. The subcellular localization of candidate proteins and HSPA5 in mouse and human retinas was evaluated by immunohistochemistry. Differences in the localization of extracellular HSPA5 were assessed between healthy human donor and atrophic age-related macular degeneration (AMD) donor eyes.

RESULTS

In the eyes of healthy mice, extracellular HSPA5 was confined to the area around the cone photoreceptor outer segments. Mapping variation in Hspa5 mRNA expression levels in the retina revealed a statistically significant trans-acting expression QTL (eQTL) on Chromosome 2 (Chr 2) and a suggestive locus on Chr 15. Sulf2 on Chr 2 was the strongest candidate gene based on partial correlation analysis, Pearson correlation with Hspa5, expression levels in the retina, a missense variant in exon 14, and its reported function in the extracellular matrix and interphotoreceptor matrix. SULF2 is localized to the rod and cone photoreceptors in both human and mouse retinas. In human retinas with no pathology, extracellular HSPA5 was localized around many cones within the macular area. In contrast, fewer HSPA5-immunopositive cones were observed in the retinas from AMD donors.

CONCLUSIONS

We identified Sulf2 as a candidate gene modulating the Hspa5 expression in the retina. The preferential loss of HSPA5 in the interphotoreceptor matrix around cone photoreceptors in atrophic AMD retinas opens up new avenues for exploring the changes in interphotoreceptor matrix (IPM) that are associated with macular disease.

Fine Mapping of a Dravet Syndrome Modifier Locus on Mouse Chromosome 5 and Candidate Gene Analysis by RNA-Seq

A substantial number of mutations have been identified in voltage-gated sodium channel genes that result in various forms of human epilepsy. SCN1A mutations result in a spectrum of severity ranging from mild febrile seizures to Dravet syndrome, an infant-onset epileptic encephalopathy. Dravet syndrome patients experience multiple seizures types that are often refractory to treatment, developmental delays, and elevated risk for SUDEP. The same sodium channel mutation can produce epilepsy phenotypes of varying clinical severity. This suggests that other factors, including genetic, modify the primary mutation and change disease severity. Mouse models provide a useful tool in studying the genetic basis of epilepsy. The mouse strain background can alter phenotype severity, supporting a contribution of genetic modifiers in epilepsy. The Scn1a^+/- mouse model has a strain-dependent epilepsy phenotype. Scn1a^+/- mice on the 129S6/SvEvTac (129) strain have a normal phenotype and lifespan, while [129xC57BL/6J]F1-Scn1a^+/- mice experience spontaneous seizures, hyperthermia-induced seizures and high rates of premature death. We hypothesize the phenotypic differences are due to strain-specific genetic modifiers that influence expressivity of the Scn1a^+/- phenotype. Low resolution mapping of Scn1a^+/- identified several Dravet syndrome modifier (Dsm) loci responsible for the strain-dependent difference in survival. One locus of interest, Dsm1 located on chromosome 5, was fine mapped to a 9 Mb region using interval specific congenics. RNA-Seq was then utilized to identify candidate modifier genes within this narrowed region. Three genes with significant total gene expression differences between 129S6/SvEvTac and [129xC57BL/6J]F1 were identified, including the GABA_A receptor subunit, Gabra2. Further analysis of Gabra2 demonstrated allele-specific expression. Pharmological manipulation by clobazam, a common anticonvulsant with preferential affinity for the GABRA2 receptor, revealed dose-dependent protection against hyperthermia-induced seizures in Scn1a^+/- mice. These findings support Gabra2 as a genetic modifier of the Scn1a^+/- mouse model of Dravet syndrome.

Epilepsy is a neurological disorder affecting approximately 3 million Americans and 1% of the worldwide population. Approximately 70% of patients diagnosed with epilepsy have a genetic basis for their disease. The same genetic mutation can result in epilepsy with varying clinical severity in some individuals. This suggests that additional factors modify the effect of the primary mutation, resulting in the differences observed. Mouse models of epilepsy also exhibit variable severity depending on the inbred strain background. This is a useful model system that enables us to determine other genetic factors that influence phenotype presentation. This study investigates a mouse model of Dravet syndrome. Similar to Dravet syndrome patients, the mice experience a variety of seizure types and have increased risk for premature death. Interestingly, the Dravet syndrome phenotype is completely masked on one mouse strain background. Through a series of genetic and pharmacological techniques we were able to identify a gene that is likely responsible for the strain-dependent phenotype difference of Dravet mice. Identifying genes that modify diseases can help predict patient outcomes and suggest new therapies for the treatment of epilepsy.

The clobazam metabolite N-desmethyl clobazam is an α2 preferring benzodiazepine with an improved therapeutic window for antihyperalgesia

Data from genetically modified mice suggest that benzodiazepine (BDZ)-site agonists with improved selectivity for α2-subtype GABAA receptors (α2GABAAR) are potentially useful for the treatment of neuropathic pain. Subtype-selective compounds available for preclinical tests in rodents support this concept but have not been approved for human use, hindering proof-of-concept studies in patients. We recently proposed that N-desmethyl clobazam (NDMC), the main metabolite of the licensed BDZ clobazam (CBZ), is responsible for most of the antihyperalgesia observed in mice after CBZ administration. In order to assess a potentially favorable pharmacological profile of NDMC, we analyzed differences in the GABAAR subtype specificity of CBZ, NDMC and diazepam (DZP) in recombinant receptors. DZP and CBZ potentiated sedating α1GABAARs and antihyperalgesic α2GABAARs with similar efficacies, whereas NDMC preferred α2GABAARs over α1GABAARs across a wide concentration range. In vivo, DZP and NDMC reduced neuropathic pain at doses between 3 and 30 mg/kg. At these doses, DZP had strong locomotor sedating effects while NDMC caused no or only weak sedation. Sedative effects of NDMC became apparent when the action of NDMC was restricted to α1GABAARs. However, when GABAAR point-mutated mice were studied that allow the analysis of antihyperalgesia and sedation in isolation, we found that, compared to DZP, NDMC had a significantly improved therapeutic window, consistent with its more favorable α2/α1 in vitro activity ratio. Given that NDMC should share the safety profile of its parent compound CBZ, it should be well-suited for proof-of-concept studies in human volunteers or patients.

The abundance of cis-acting loci leading to differential allele expression in F1 mice and their relationship to loci harboring genes affecting complex traits

Background

Genome-wide surveys have detected cis-acting quantitative trait loci altering levels of RNA transcripts (RNA-eQTLs) by associating SNV alleles to transcript levels. However, the sensitivity and specificity of detection of cis- expression quantitative trait loci (eQTLs) by genetic approaches, reliant as it is on measurements of transcript levels in recombinant inbred strains or offspring from arranged crosses, is unknown, as is their relationship to QTL’s for complex phenotypes.

Results

We used transcriptome-wide differential allele expression (DAE) to detect cis-eQTLs in forebrain and kidney from reciprocal crosses between three mouse inbred strains, 129S1/SvlmJ, DBA/2J, and CAST/EiJ and C57BL/6 J. Two of these crosses were previously characterized for cis-eQTLs and QTLs for various complex phenotypes by genetic analysis of recombinant inbred (RI) strains. 5.4 %, 1.9 % and 1.5 % of genes assayed in forebrain of B6/129SF1, B6/DBAF1, and B6/CASTF1 mice, respectively, showed differential allelic expression, indicative of cis-acting alleles at these genes. Moreover, the majority of DAE QTLs were observed to be tissue-specific with only a small fraction showing cis-effects in both tissues. Comparing DAE QTLs in F1 mice to cis-eQTLs previously mapped in RI strains we observed that many of the cis-eQTLs were not confirmed by DAE. Additionally several novel DAE-QTLs not identified as cis-eQTLs were identified suggesting that there are differences in sensitivity and specificity for QTL detection between the two methodologies. Strain specific DAE QTLs in B6/DBAF1 mice were located in excess at candidate genes for alcohol use disorders, seizures, and angiogenesis previously implicated by genetic linkage in C57BL/6J × DBA/2JF2 mice or BXD RI strains.

Conclusions

Via a survey for differential allele expression in F1 mice, a substantial proportion of genes were found to have alleles altering expression in cis-acting fashion. Comparing forebrain and kidney, many or most of these alleles were tissue-specific in action. The identification of strain specific DAE QTLs, can assist in assessment of candidate genes located within the large intervals associated with trait QTLs.

Electronic supplementary material

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

Mutations in GABRG2 receptor gene are not a major factor in the pathogenesis of mesial temporal lobe epilepsy in Indian population

AIM

This study is focused on GABRG2 gene sequence variations in patients with mesial temporal lobe epilepsy (mTLE). The GABAAreceptor is a heteropentameric receptor and alpha-1 beta-2 gamma-2 subunits combination is most abundant and present in almost all regions of the brain. The gamma-2 subunit (GABRG2) gene mutations have been reported in different epilepsy pathologies. In the present study we have looked for GABRG2 gene sequence variations in patients with mTLE.

MATERIALS AND METHODS

Twenty patients (12 females and eight males, age 4.6-38 years) with MTLE were recruited for this investigation. Patients were recommended for epilepsy surgery after all clinical investigations as per the epilepsy protocol. Ethnically matched glioma or meningioma patients were considered as nonepileptic controls. During temporal lobectomy of amygdalohippocampectomy, hippocampal brain tissue samples were resected guided by intraoperative electrocorticography (ECoG) activity. All 11 exons of GABRG2 gene with their flanking intronic regions were amplified by polymerase chain reaction (PCR) and screened by DNA sequencing analysis for sequence variations.

STATISTICAL ANALYSIS USED

Comparison of allele frequencies between patient and control groups was determined using a c(2) test.

RESULTS AND CONCLUSIONS

Total five DNA sequence variations were identified, three in exonic regions (c.643A > G, rs211035), (c.T > A, rs424740), and (c.C > T, rs418210) and two in intronic regions (c.751 + 41A > G, rs211034) and (c.751 + 52G > A, rs 34281163). Allele frequencies of variants identified in this study did not differ between patients and normal controls. Thus, we conclude that GABRG2 gene may not be playing significant role in the development of epilepsy or as a susceptibility gene in patients with MTLE in Indian population.

Cross-Species Integrative Functional Genomics in GeneWeaver Reveals a Role for Pafah1b1 in Altered Response to Alcohol

Identifying the biological substrates of complex neurobehavioral traits such as alcohol dependency pose a tremendous challenge given the diverse model systems and phenotypic assessments used. To address this problem we have developed a platform for integrated analysis of high-throughput or genome-wide functional genomics studies. A wealth of such data exists, but it is often found in disparate, non-computable forms. Our interactive web-based software system, Gene Weaver (http://www.geneweaver.org), couples curated results from genomic studies to graph-theoretical tools for combinatorial analysis. Using this system we identified a gene underlying multiple alcohol-related phenotypes in four species. A search of over 60,000 gene sets in GeneWeaver's database revealed alcohol-related experimental results including genes identified in mouse genetic mapping studies, alcohol selected Drosophila lines, Rattus differential expression, and human alcoholic brains. We identified highly connected genes and compared these to genes currently annotated to alcohol-related behaviors and processes. The most highly connected gene not annotated to alcohol was Pafah1b1. Experimental validation using a Pafah1b1 conditional knock-out mouse confirmed that this gene is associated with an increased preference for alcohol and an altered thermoregulatory response to alcohol. Although this gene has not been previously implicated in alcohol-related behaviors, its function in various neural mechanisms makes a role in alcohol-related phenomena plausible. By making diverse cross-species functional genomics data readily computable, we were able to identify and confirm a novel alcohol-related gene that may have implications for alcohol use disorders and other effects of alcohol.

Multipronged approach to identify and validate a novel upstream regulator of Sncg in mouse retinal ganglion cells

Loss of retinal ganglion cells (RGCs) is one of the hallmarks of retinal neurodegenerative diseases, glaucoma being one of the most common. Mechanistic studies on RGCs are hindered by the lack of sufficient primary cells and consensus regarding their signature markers. Recently, γ-synuclein (SNCG) has been shown to be highly expressed in the somas and axons of RGCs. In various mouse models of glaucoma, downregulation of Sncg gene expression correlates with RGC loss. To investigate the role of Sncg in RGCs, we used a novel systems genetics approach to identify a gene that modulates Sncg expression, followed by confirmatory studies in both healthy and diseased retinae. We found that chromosome 1 harbors an expression quantitative trait locus that modulates Sncg expression in the mouse retina, and identified the prefoldin-2 (PFDN2) gene as the candidate upstream modulator of Sncg expression. Our immunohistochemical analyses revealed similar expression patterns in both mouse and human healthy retinae, with PFDN2 colocalizing with SNCG in RGCs and their axons. In contrast, in retinae from glaucoma subjects, SNCG levels were significantly reduced, although PFDN2 levels were maintained. Using a novel flow cytometry-based RGC isolation method, we obtained viable populations of murine RGCs. Knocking down Pfdn2 expression in primary murine RGCs significantly reduced Sncg expression, confirming that Pfdn2 regulates Sncg expression in murine RGCs. Gene Ontology analysis indicated shared mitochondrial function associated with Sncg and Pfdn2. These data solidify the relationship between Sncg and Pfdn2 in RGCs, and provide a novel mechanism for maintaining RGC health.

A Murine Hypertrophic Cardiomyopathy Model: The DBA/2J Strain

Familial hypertrophic cardiomyopathy (HCM) is attributed to mutations in genes that encode for the sarcomere proteins, especially Mybpc3 and Myh7. Genotype-phenotype correlation studies show significant variability in HCM phenotypes among affected individuals with identical causal mutations. Morphological changes and clinical expression of HCM are the result of interactions with modifier genes. With the exceptions of angiotensin converting enzyme, these modifiers have not been identified. Although mouse models have been used to investigate the genetics of many complex diseases, natural murine models for HCM are still lacking. In this study we show that the DBA/2J (D2) strain of mouse has sequence variants in Mybpc3 and Myh7, relative to widely used C57BL/6J (B6) reference strain and the key features of human HCM. Four-month-old of male D2 mice exhibit hallmarks of HCM including increased heart weight and cardiomyocyte size relative to B6 mice, as well as elevated markers for cardiac hypertrophy including β-myosin heavy chain (MHC), atrial natriuretic peptide (ANP), brain natriuretic peptide (BNP), and skeletal muscle alpha actin (α1-actin). Furthermore, cardiac interstitial fibrosis, another feature of HCM, is also evident in the D2 strain, and is accompanied by up-regulation of type I collagen and α-smooth muscle actin (SMA)-markers of fibrosis. Of great interest, blood pressure and cardiac function are within the normal range in the D2 strain, demonstrating that cardiac hypertrophy and fibrosis are not secondary to hypertension, myocardial infarction, or heart failure. Because D2 and B6 strains have been used to generate a large family of recombinant inbred strains, the BXD cohort, the D2 model can be effectively exploited for in-depth genetic analysis of HCM susceptibility and modifier screens.

Abnormal subcellular localization of GABAA receptor subunits in schizophrenia brain

Inhibitory neurotransmission is primarily mediated by γ-aminobutyric acid (GABA) activating synaptic GABA type A receptors (GABA(A)R). In schizophrenia, presynaptic GABAergic signaling deficits are among the most replicated findings; however, postsynaptic GABAergic deficits are less well characterized. Our lab has previously demonstrated that although there is no difference in total protein expression of the α1-6, β1-3 or γ2 GABA(A)R subunits in the superior temporal gyrus (STG) in schizophrenia, the α1, β1 and β2 GABA(A)R subunits are abnormally N-glycosylated. N-glycosylation is a posttranslational modification that has important functional roles in protein folding, multimer assembly and forward trafficking. To investigate the impact that altered N-glycosylation has on the assembly and trafficking of GABA(A)Rs in schizophrenia, this study used western blot analysis to measure the expression of α1, α2, β1, β2 and γ2 GABA(A)R subunits in subcellular fractions enriched for endoplasmic reticulum (ER) and synapses (SYN) from STG of schizophrenia (N = 16) and comparison (N = 14) subjects and found evidence of abnormal localization of the β1 and β2 GABA(A)R subunits and subunit isoforms in schizophrenia. The β2 subunit is expressed as three isoforms at 52 kDa (β2(52 kDa)), 50 kDa (β2(50 kDa)) and 48 kDa (β2(48 kDa)). In the ER, we found increased total β2 GABA(A)R subunit (β2(ALL)) expression driven by increased β2(50 kDa), a decreased ratio of β(248 kDa):β2(ALL) and an increased ratio of β2(50 kDa):β2(48 kDa). Decreased ratios of β1:β2(ALL) and β1:β2(50 kDa) in both the ER and SYN fractions and an increased ratio of β2(52 kDa):β(248 kDa) at the synapse were also identified in schizophrenia. Taken together, these findings provide evidence that alterations of N-glycosylation may contribute to GABAergic signaling deficits in schizophrenia by disrupting the assembly and trafficking of GABA(A)Rs.

RETRACTED: Transcutaneous vagus nerve stimulation (tVNS) enhances response selection during action cascading processes

This article has been retracted: please see Elsevier Policy on Article Withdrawal (https://www.elsevier.com/about/our-business/policies/article-withdrawal). This retraction comes after a thorough and time-consuming analysis of the published article and the data it reports, following an investigation from the University of Leiden. On date 15th November 2021 European Neuropsychopharmacology received advice from University of Leiden, containing allegations of misconduct concerning the data reported in this article. Due to the nature of the concerns raised about this paper, the Editor-in-Chief examined all documents provided, all aspects of the peer review and editorial process and requested the view of each of the authors. The Editor-in-Chief wishes to acknowledge the co-operation of all parties. The Editor-in-Chief considers that there is a legitimate cause for concern regarding deletion of participants from datasets. The Editor-in-Chief commends the authors for their willingness in participating in this dialogue. However, given the concerns raised, the suggestion of the University of Leiden, and the view expressed by some of the authors, Retraction of this paper is inevitable. The Editor-in-Chief and the entire editorial board will use this case as a reminder to be as diligent as possible in the peer review process and to support transparency and openness in data sharing.

Aberrant DNA methylation of imprinted loci in hepatocellular carcinoma and after in vitro exposure to common risk factors

BACKGROUND

Hepatocellular carcinoma (HCC) is among the most frequent human malignancies and a major cause of cancer-related death worldwide. It is characterized by late detection and fast progression, and it is believed that epigenetic disruption may be one of the molecular mechanisms leading to hepatocarcinogenesis. Previous studies from our group revealed that HCC tumors exhibit specific DNA methylation signatures associated with major risk factors and tumor progression. Imprinted genes are mono-allelically expressed in a parent-of-origin-dependent manner and have been suggested to be more susceptible to deregulation in cancer. To test this notion, we performed a targeted analysis of DNA methylation in known imprinted genes, using HCC samples and in vitro models of carcinogenic exposure.

RESULTS

Analysis of HCC DNA methylation in two independent datasets showed that differentially methylated loci are significantly enriched in imprinted genes. Most of the promoters of imprinted genes were found hypomethylated in HCC tumors compared to surrounding tissues, contrasting with the frequent promoter hypermethylation observed in tumors. We next investigated the status of methylation of the imprinting control region (ICR) of different imprinted clusters and found that the 15q11-13 ICR was significantly hypomethylated in tumors relative to their surrounding tissues. In addition, expression of imprinted genes within this cluster was frequently deregulated in a gene-specific manner, suggesting distinct mechanisms of regulation in this region. Finally, primary human hepatocytes and hepatocyte-like HepaRG cells displayed higher methylation variability in certain imprinted loci after natural hepatitis B virus (HBV) infection and after lipid accumulation, respectively.

CONCLUSION

The methylation status of a large panel of imprinted genes was found deregulated in HCC, suggesting a major role of this mechanism during hepatocarcinogenesis. In vitro models support the hypothesis of imprinted gene methylation as a potential marker of environmental exposures.

Joint genetic analysis of hippocampal size in mouse and human identifies a novel gene linked to neurodegenerative disease

BACKGROUND

Variation in hippocampal volume has been linked to significant differences in memory, behavior, and cognition among individuals. To identify genetic variants underlying such differences and associated disease phenotypes, multinational consortia such as ENIGMA have used large magnetic resonance imaging (MRI) data sets in human GWAS studies. In addition, mapping studies in mouse model systems have identified genetic variants for brain structure variation with great power. A key challenge is to understand how genetically based differences in brain structure lead to the propensity to develop specific neurological disorders.

RESULTS

We combine the largest human GWAS of brain structure with the largest mammalian model system, the BXD recombinant inbred mouse population, to identify novel genetic targets influencing brain structure variation that are linked to increased risk for neurological disorders. We first use a novel cross-species, comparative analysis using mouse and human genetic data to identify a candidate gene, MGST3, associated with adult hippocampus size in both systems. We then establish the coregulation and function of this gene in a comprehensive systems-analysis.

CONCLUSIONS

We find that MGST3 is associated with hippocampus size and is linked to a group of neurodegenerative disorders, such as Alzheimer's.

Short intronic repeat sequences facilitate circular RNA production

Recent deep sequencing studies have revealed thousands of circular noncoding RNAs generated from protein-coding genes. The mechanism by which the spliceosome selects only certain exons to circularize is largely unknown. Liang and Wilusz used extensive mutagenesis of expression plasmids to show that miniature introns containing the splice sites along with short (∼30- to 40-nt) inverted repeats are sufficient to allow the intervening exons to circularize in cells. The intronic repeats and exonic sequences must collaborate with one another, and a functional 3′ end processing signal is required, suggesting that circularization may occur post-transcriptionally.

Recent deep sequencing studies have revealed thousands of circular noncoding RNAs generated from protein-coding genes. These RNAs are produced when the precursor messenger RNA (pre-mRNA) splicing machinery “backsplices” and covalently joins, for example, the two ends of a single exon. However, the mechanism by which the spliceosome selects only certain exons to circularize is largely unknown. Using extensive mutagenesis of expression plasmids, we show that miniature introns containing the splice sites along with short (∼30- to 40-nucleotide) inverted repeats, such as Alu elements, are sufficient to allow the intervening exons to circularize in cells. The intronic repeats must base-pair to one another, thereby bringing the splice sites into close proximity to each other. More than simple thermodynamics is clearly at play, however, as not all repeats support circularization, and increasing the stability of the hairpin between the repeats can sometimes inhibit circular RNA biogenesis. The intronic repeats and exonic sequences must collaborate with one another, and a functional 3′ end processing signal is required, suggesting that circularization may occur post-transcriptionally. These results suggest detailed and generalizable models that explain how the splicing machinery determines whether to produce a circular noncoding RNA or a linear mRNA.

Effects of L-theanine on posttraumatic stress disorder induced changes in rat brain gene expression

Posttraumatic stress disorder (PTSD) is characterized by the occurrence of a traumatic event that is beyond the normal range of human experience. The future of PTSD treatment may specifically target the molecular mechanisms of PTSD. In the US, approximately 20% of adults report taking herbal products to treat medical illnesses. L-theanine is the amino acid in green tea primarily responsible for relaxation effects. No studies have evaluated the potential therapeutic properties of herbal medications on gene expression in PTSD. We evaluated gene expression in PTSD-induced changes in the amygdala and hippocampus of Sprague-Dawley rats. The rats were assigned to PTSD-stressed and nonstressed groups that received either saline, midazolam, L-theanine, or L-theanine + midazolam. Amygdala and hippocampus tissue samples were analyzed for changes in gene expression. One-way ANOVA was used to detect significant difference between groups in the amygdala and hippocampus. Of 88 genes examined, 17 had a large effect size greater than 0.138. Of these, 3 genes in the hippocampus and 5 genes in the amygdala were considered significant (P < 0.05) between the groups. RT-PCR analysis revealed significant changes between groups in several genes implicated in a variety of disorders ranging from PTSD, anxiety, mood disorders, and substance dependence.

Identification of a QTL in Mus musculus for alcohol preference, withdrawal, and Ap3m2 expression using integrative functional genomics and precision genetics

Extensive genetic and genomic studies of the relationship between alcohol drinking preference and withdrawal severity have been performed using animal models. Data from multiple such publications and public data resources have been incorporated in the GeneWeaver database with >60,000 gene sets including 285 alcohol withdrawal and preference-related gene sets. Among these are evidence for positional candidates regulating these behaviors in overlapping quantitative trait loci (QTL) mapped in distinct mouse populations. Combinatorial integration of functional genomics experimental results revealed a single QTL positional candidate gene in one of the loci common to both preference and withdrawal. Functional validation studies in Ap3m2 knockout mice confirmed these relationships. Genetic validation involves confirming the existence of segregating polymorphisms that could account for the phenotypic effect. By exploiting recent advances in mouse genotyping, sequence, epigenetics, and phylogeny resources, we confirmed that Ap3m2 resides in an appropriately segregating genomic region. We have demonstrated genetic and alcohol-induced regulation of Ap3m2 expression. Although sequence analysis revealed no polymorphisms in the Ap3m2-coding region that could account for all phenotypic differences, there are several upstream SNPs that could. We have identified one of these to be an H3K4me3 site that exhibits strain differences in methylation. Thus, by making cross-species functional genomics readily computable we identified a common QTL candidate for two related bio-behavioral processes via functional evidence and demonstrate sufficiency of the genetic locus as a source of variation underlying two traits.

Genetics of gene expression in CNS

Transcriptome studies have revealed a surprisingly high level of variation among individuals in expression of key genes in the CNS under both normal and experimental conditions. Ten-fold variation is common, yet the specific causes and consequences of this variation are largely unknown. By combining classic gene mapping methods-family linkage studies and genomewide association-with high-throughput genomics, it is now possible to define quantitative trait loci (QTLs), single-gene variants, and even single SNPs and indels that control gene expression in different brain regions and cells. This review considers some of the major technical and conceptual challenges in analyzing variation in expression in the CNS with a focus on mRNAs, rather than noncoding RNAs or proteins. At one level of analysis, this work has been highly successful, and we finally have techniques that can be used to track down small numbers of loci that control expression in the CNS. But at a higher level of analysis, we still do not understand the genetic architecture of gene expression in brain, the consequences of expression QTLs on protein levels or on cell function, or the combined impact of expression differences on behavior and disease risk. These important gaps are likely to be bridged over the next several decades using (1) much larger sample sizes, (2) more powerful RNA sequencing and proteomic methods, and (3) novel statistical and computational models to predict genome-to-phenome relations.

Out of the bottleneck: the Diversity Outcross and Collaborative Cross mouse populations in behavioral genetics research

The historical origins of classical laboratory mouse strains have led to a relatively limited range of genetic and phenotypic variation, particularly for the study of behavior. Many recent efforts have resulted in improved diversity and precision of mouse genetic resources for behavioral research, including the Collaborative Cross and Diversity Outcross population. These two populations, derived from an eight way cross of common and wild-derived strains, have high precision and allelic diversity. Behavioral variation in the population is expanded, both qualitatively and quantitatively. Variation that had once been canalized among the various inbred lines has been made amenable to genetic dissection. The genetic attributes of these complementary populations, along with advances in genetic and genomic technologies, makes a systems genetic analyses of behavior more readily tractable, enabling discovery of a greater range of neurobiological phenomena underlying behavioral variation.

Expression, covariation, and genetic regulation of miRNA Biogenesis genes in brain supports their role in addiction, psychiatric disorders, and disease

The role of miRNA and miRNA biogenesis genes in the adult brain is just beginning to be explored. In this study we have performed a comprehensive analysis of the expression, genetic regulation, and co-expression of major components of the miRNA biogenesis pathway using human and mouse data sets and resources available on the GeneNetwork web site (genenetwork.org). We found a wide range of variation in expression in both species for key components of the pathway—Drosha, Pasha, and Dicer. Across species, tissues, and expression platforms all three genes are generally well-correlated. No single genetic locus exerts a strong and consistent influence on the expression of these key genes across murine brain regions. However, in mouse striatum, many members of the miRNA pathway are correlated—including Dicer, Drosha, Pasha, Ars2 (Srrt), Eif2c1 (Ago1), Eif2c2 (Ago2), Zcchc11, and Snip1. The expression of these genes may be partly influenced by a locus on Chromosome 9 (105.67–106.32 Mb). We explored ~1500 brain phenotypes available for the C57BL/6J × DBA/2J (BXD) genetic mouse population in order to identify miRNA biogenesis genes correlated with traits related to addiction and psychiatric disorders. We found a significant association between expression of Dicer and Drosha in several brain regions and the response to many drugs of abuse, including ethanol, cocaine, and methamphetamine. Expression of Dicer, Drosha, and Pasha in most of the brain regions explored is strongly correlated with the expression of key members of the dopamine system. Drosha, Pasha, and Dicer expression is also correlated with the expression of behavioral traits measuring depression and sensorimotor gating, impulsivity, and anxiety, respectively. Our study provides a global survey of the expression and regulation of key miRNA biogenesis genes in brain and provides preliminary support for the involvement of these genes and their product miRNAs in addiction and psychiatric disease processes.

Genetic and molecular network analysis of behavior

This chapter provides an introduction into the genetic control and analysis of behavioral variation using powerful online resources. We introduce you to the new field of systems genetics using "case studies" drawn from the world of behavioral genetics that exploit populations of genetically diverse lines of mice. These lines differ very widely in patterns of gene and protein expression in the brain and in patterns of behavior. In this chapter, we address the following set of related questions: (1) Can we combine massive genomic data sets with large aggregates of precise quantitative data on behavior? (2) Can we map causal relations between gene variants and behavioral differences? (3) Can we simultaneously use these highly coherent data sets to understand more about the underlying molecular and cellular basis of behavior?