Genetics of gene expression in CNS

Acute intrahippocampal administration of melanin-concentrating hormone impairs memory consolidation and decreases the expression of MCHR-1 and TrkB receptors

Interest in the role of melanin-concentrating hormone (MCH) in memory processes has increased in recent years, with some studies reporting memory-enhancing effects, while others report deleterious effects. Due to these discrepancies, this study seeks to provide new evidence about the role of MCH in memory consolidation and its relation with BDNF/TrkB system. To this end, in the first experiment, increased doses of MCH were acutely administered in both hippocampi to groups of male rats (25, 50, 200, and 500 ng). Microinjections were carried out immediately after finishing the sample trial of two hippocampal-dependent behavioral tasks: the Novel Object Recognition Test (NORT) and the modified Elevated Plus Maze (mEPM) test. Results indicated that a dose of 200 ng of MCH or higher impaired memory consolidation in both tasks. A second experiment was performed in which a dose of 200 ng of MCH was administered alone or co-administered with the MCHR-1 antagonist ATC-0175 at the end of the sample trial in the NORT. Results showed that MCH impaired memory consolidation, while the co-administration with ATC-0175 reverted this detrimental effect. Moreover, MCH induced a significant decrease in hippocampal MCHR-1 and TrkB expression with no modification in the expression of BDNF and NMDA receptor subunits NR1, NR2A, and NR2B. These results suggest that MCH in vivo elicits pro-amnesic effects in the rat hippocampus by decreasing the availability of its receptor and TrkB receptors, thus linking both endogenous systems to memory processes.

The most common European HINT1 neuropathy variant phenotype and its case studies

HINT1 is an ubiquitous homodimeric purine phosphoramidase belonging to the histidine-triad superfamily. In neurons, HINT1 stabilizes the interaction of different receptors and regulates the effects of their signaling disturbances. Changes in HINT1 gene are associated with autosomal recessive axonal neuropathy with neuromyotonia. Aim of the study was detailed description of patients' phenotype with HINT1 homozygous NM_005340.7: c.110G>C (p.Arg37Pro) variant. Seven homozygous and three compound heterozygous patients were recruited and evaluated using standardized tests for CMT patients, in four patients' nerve ultrasonography was performed. The median age of symptom onset was 10 years (range 1-20), with initial complaints being distal lower limb weakness with gait impairment, combined with muscle stiffness, more pronounced in the hands than in the legs and worsened by cold. Arm muscles became involved later, presenting with distal weakness and hypotrophy. Neuromyotonia was present in all reported patients and is thus a diagnostic hallmark. Electrophysiological studies demonstrated axonal polyneuropathy. Impaired mental performance was observed in six out of ten cases. In all patients with HINT1 neuropathy, ultrasound examination showed significantly reduced muscle volume as well as spontaneous fasciculations and fibrillations. The nerve cross-sectional areas of the median and ulnar nerves were closer to the lower limits of the normal values. None of the investigated nerves had structural changes. Our findings broaden the phenotype of HINT1-neuropathy and have implications for diagnostics and ultrasonographic evaluation of HINT1-neuropathy patients.

Natural Mutations Affect Structure and Function of gC1q Domain of Otolin-1

Otolin-1 is a scaffold protein of otoliths and otoconia, calcium carbonate biominerals from the inner ear. It contains a gC1q domain responsible for trimerization and binding of Ca²⁺. Knowledge of a structure-function relationship of gC1q domain of otolin-1 is crucial for understanding the biology of balance sensing. Here, we show how natural variants alter the structure of gC1q otolin-1 and how Ca²⁺ are able to revert some effects of the mutations. We discovered that natural substitutions: R339S, R342W and R402P negatively affect the stability of apo-gC1q otolin-1, and that Q426R has a stabilizing effect. In the presence of Ca²⁺, R342W and Q426R were stabilized at higher Ca²⁺ concentrations than the wild-type form, and R402P was completely insensitive to Ca²⁺. The mutations affected the self-association of gC1q otolin-1 by inducing detrimental aggregation (R342W) or disabling the trimerization (R402P) of the protein. Our results indicate that the natural variants of gC1q otolin-1 may have a potential to cause pathological changes in otoconia and otoconial membrane, which could affect sensing of balance and increase the probability of occurrence of benign paroxysmal positional vertigo (BPPV).

Systems Genetics Approaches in Mouse Models of Group A Streptococcal Necrotizing Soft-Tissue Infections

Mouse models are invaluable resources for studying the pathogenesis and preclinical evaluation of therapeutics and vaccines against many human pathogens. Infections caused by group A streptococcus (GAS, Streptococcus pyogenes) are heterogeneous ranging from mild pharyngitis to severe invasive necrotizing fasciitis, a subgroup of necrotizing soft-tissue infections (NSTIs). While several strains of mice including BALB/c, C3H/HeN, CBA/J, and C57BL/10 offered significant insights, the human specificity and the interindividual variations on susceptibility or resistance to GAS infections limit their ability to mirror responses as seen in humans. In this chapter, we discuss the advanced recombinant inbred (ARI) BXD mouse model that mimics the genetic diversity as seen in humans and underpins the feasibility to map multiple genes (genetic loci) modulating GAS NSTI. GAS produces a myriad of virulence factors, including superantigens (SAg). Superantigens are potent immune toxins that activate T cells by cross-linking T cell receptors with human leukocyte antigen class-II (HLA-II) molecules expressed on antigen-presenting cells. This leads to a pro-inflammatory cytokine storm and the subsequent multiple organ damage and shock. Inbred mice are innately refractive to SAg-mediated responses. In this chapter, we discuss the versatility of the HLA-II transgenic mouse model that allowed the biological validation of known genetic associations to GAS NSTI. The combined utility of ARI-BXD and HLA-II mice as complementary approaches that offer clinically translatable insights into pathomechanisms driven by complex traits and host genetic context and novel means to evaluate the in vivo efficiency of therapies to improve outcomes of GAS NSTI are also discussed.

RNA sequencing profiling of the retina in C57BL/6J and DBA/2J mice: Enhancing the retinal microarray data sets from GeneNetwork

Purpose

The goal of the present study is to provide an independent assessment of the retinal transcriptome signatures of C57BL/6J (B6) and DBA/2J (D2) mice, and to enhance existing microarray data sets for accurately defining the allelic differences in the BXD recombinant inbred strains.

Methods

Retinas from B6 and D2 mice (three of each) were used for the RNA sequencing (RNA-seq) analysis. Transcriptome features were examined for both strains. Differentially expressed genes between the two strains were identified, and bioinformatic analysis was performed to analyze the transcriptome differences between the B6 and D2 strains, including Gene Ontology (GO) analysis, Phenotype and Reactome enrichment, and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. The RNA-seq data were then directly compared with one of the microarray data sets (Department of Defense [DoD] Retina Normal Affy MoGene 2.0 ST RMA Gene Level Microarray Database) hosted on GeneNetwork.

Results

RNA-seq provided an in-depth analysis of the transcriptome of the B6 and D2 retinas with a total of more than 30,000,000 reads per sample. More than 70% of the reads were uniquely mapped, resulting in a total of 18,100 gene counts for all six samples. A total of 1,665 genes were differentially expressed, with 858 of these more highly expressed in the B6 retinas and 807 more highly expressed in the D2 retinas. Several molecular pathways were differentially active between the two strains, including the retinoic acid metabolic process, endoplasmic reticulum lumen, extracellular matrix (ECM) organization, and the PI3K-Akt signaling pathway. The most enriched KEGG pathways were the pentose and glucuronate interconversions pathway, the cytochrome P450 pathway, the protein digestion and absorption pathway, and the ECM-receptor interaction pathway. Each of these pathways had a more than fourfold enrichment. The DoD Normal Retina Microarray Database provided expression profiling for 26,191 annotated transcripts for B6 mouse, D2 mouse, and 53 BXD strains. A total of 13,793 genes in this microarray data set were comparable to the RNA-seq data set. For the B6 and D2 retinas, the RNA-seq data and the microarray data were highly correlated with each other (Pearson's r=0.780 for the B6 mice and 0.784 for D2 mice). These results suggest that the microarray data set can reliably detect differentially expressed genes between the B6 and D2 retinas, with an overall accuracy of 91.1%. Examples of true positive and false positive genes are provided.

Conclusions

Retinal transcriptome features of B6 and D2 mouse strains provide a useful reference for a better understanding of the mouse retina. Generally, the microarray database presented on GeneNetwork shows good agreement with the RNA-seq data, but we note that any allelic difference between B6 and D2 mice should be verified with the latter.

Evaluation of Sirtuin-3 probe quality and co-expressed genes using literature cohesion

BACKGROUND

Gene co-expression studies can provide important insights into molecular and cellular signaling pathways. The GeneNetwork database is a unique resource for co-expression analysis using data from a variety of tissues across genetically distinct inbred mice. However, extraction of biologically meaningful co-expressed gene sets is challenging due to variability in microarray platforms, probe quality, normalization methods, and confounding biological factors. In this study, we tested whether literature derived functional cohesion could be used as an objective metric in lieu of 'ground truth' to evaluate the quality of probes and microarray datasets.

RESULTS

We examined Sirtuin-3 (Sirt3) co-expressed gene sets extracted from either liver or brain tissues of BXD recombinant inbred mice in the GeneNetwork database. Depending on the microarray platform, there were as many as 26 probes that targeted different regions of Sirt3 primary transcript. Co-expressed gene sets (ranging from 100-1000 genes) associated with each Sirt3 probe were evaluated using the previously developed literature-derived cohesion p-value (LPv) and benchmarked against 'gold standards' derived from proteomic studies or Gene Ontology classifications. We found that the maximal F-measure was obtained at an average window size of 535 genes. Using set size of 500 genes, the Pearson correlations between LPv and F-measure as well as between LPv and mitochondrial gene enrichment p-values were 0.90 and 0.93, respectively. Importantly, we found that the LPv approach can distinguish high quality Sirt3 probes. Analysis of the most functionally cohesive Sirt3 co-expressed gene set revealed core metabolic pathways that were shared between hippocampus and liver as well as distinct pathways which were unique to each tissue. These results are consistent with other studies that suggest Sirt3 is a key metabolic regulator and has distinct functions in energy-producing vs. energy-demanding tissues.

CONCLUSIONS

Our results provide proof-of-concept that literature cohesion analysis is useful for evaluating the quality of probes and microarray datasets, particularly when experimentally derived gold standards are unavailable. Our approach would enable researchers to rapidly identify biologically meaningful co-expressed gene sets and facilitate discovery from high throughput genomic data.

C57BL/6 substrain differences in inflammatory and neuropathic nociception and genetic mapping of a major quantitative trait locus underlying acute thermal nociception

Sensitivity to different pain modalities has a genetic basis that remains largely unknown. Employing closely related inbred mouse substrains can facilitate gene mapping of nociceptive behaviors in preclinical pain models. We previously reported enhanced sensitivity to acute thermal nociception in C57BL/6J (B6J) versus C57BL/6N (B6N) substrains. Here, we expanded on nociceptive phenotypes and observed an increase in formalin-induced inflammatory nociceptive behaviors and paw diameter in B6J versus B6N mice (Charles River Laboratories). No strain differences were observed in mechanical or thermal hypersensitivity or in edema following the Complete Freund's Adjuvant model of inflammatory pain, indicating specificity in the inflammatory nociceptive stimulus. In the chronic constrictive nerve injury, a model of neuropathic pain, no strain differences were observed in baseline mechanical threshold or in mechanical hypersensitivity up to one month post-chronic constrictive nerve injury. We replicated the enhanced thermal nociception in the 52.5°C hot plate test in B6J versus B6N mice from The Jackson Laboratory. Using a B6J × B6N-F2 cross (N = 164), we mapped a major quantitative trait locus underlying hot plate sensitivity to chromosome 7 that peaked at 26 Mb (log of the odds [LOD] = 3.81, p < 0.01; 8.74 Mb-36.50 Mb) that was more pronounced in males. Genes containing expression quantitative trait loci associated with the peak nociceptive marker that are implicated in pain and inflammation include Ryr1, Cyp2a5, Pou2f2, Clip3, Sirt2, Actn4, and Ltbp4 (false discovery rate < 0.05). Future studies involving positional cloning and gene editing will determine the quantitative trait gene(s) and potential pleiotropy of this locus across pain modalities.

Mapping Molecular Datasets Back to the Brain Regions They are Extracted from: Remembering the Native Countries of Hypothalamic Expatriates and Refugees

This article focuses on approaches to link transcriptomic, proteomic, and peptidomic datasets mined from brain tissue to the original locations within the brain that they are derived from using digital atlas mapping techniques. We use, as an example, the transcriptomic, proteomic and peptidomic analyses conducted in the mammalian hypothalamus. Following a brief historical overview, we highlight studies that have mined biochemical and molecular information from the hypothalamus and then lay out a strategy for how these data can be linked spatially to the mapped locations in a canonical brain atlas where the data come from, thereby allowing researchers to integrate these data with other datasets across multiple scales. A key methodology that enables atlas-based mapping of extracted datasets-laser-capture microdissection-is discussed in detail, with a view of how this technology is a bridge between systems biology and systems neuroscience.

Unraveling long non-coding RNAs through analysis of high-throughput RNA-sequencing data

Extensive genome-wide transcriptome study mediated by high throughput sequencing technique has revolutionized the study of genetics and epigenetic at unprecedented resolution. The research has revealed that besides protein-coding RNAs, large proportions of mammalian transcriptome includes a heap of regulatory non protein-coding RNAs, the number encoded within human genome is enigmatic. Many taboos developed in the past categorized these non-coding RNAs as ''dark matter" and "junks". Breaking the myth, RNA-seq-- a recently developed experimental technique is widely being used for studying non-coding RNAs which has acquired the limelight due to their physiological and pathological significance. The longest member of the ncRNA family-- long non-coding RNAs, acts as stable and functional part of a genome, guiding towards the important clues about the varied biological events like cellular-, structural- processes governing the complexity of an organism. Here, we review the most recent and influential computational approach developed to identify and quantify the long non-coding RNAs serving as an assistant for the users to choose appropriate tools for their specific research.

Alterations in the expression of a neurodevelopmental gene exert long-lasting effects on cognitive-emotional phenotypes and functional brain networks: translational evidence from the stress-resilient Ahi1 knockout mouse

Many psychiatric disorders are highly heritable and may represent the clinical outcome of early aberrations in the formation of neural networks. The placement of brain connectivity as an 'intermediate phenotype' renders it an attractive target for exploring its interaction with genomics and behavior. Given the complexity of genetic make up and phenotypic heterogeneity in humans, translational studies are indicated. Recently, we demonstrated that a mouse model with heterozygous knockout of the key neurodevelopmental gene Ahi1 displays a consistent stress-resilient phenotype. Extending these data, the current research describes our multi-faceted effort to link early variations in Ahi1 expression with long-term consequences for functional brain networks and cognitive-emotional phenotypes. By combining behavioral paradigms with graph-based analysis of whole-brain functional networks, and then cross-validating the data with robust neuroinformatic data sets, our research suggests that physiological variation in gene expression during neurodevelopment is eventually translated into a continuum of global network metrics that serve as intermediate phenotypes. Within this framework, we suggest that organization of functional brain networks may result, in part, from an adaptive trade-off between efficiency and resilience, ultimately culminating in a phenotypic diversity that encompasses dimensions such as emotional regulation and cognitive function.

Identifying genes for neurobehavioural traits in rodents: progress and pitfalls

Identifying genes and pathways that contribute to differences in neurobehavioural traits is a key goal in psychiatric research. Despite considerable success in identifying quantitative trait loci (QTLs) associated with behaviour in laboratory rodents, pinpointing the causal variants and genes is more challenging. For a long time, the main obstacle was the size of QTLs, which could encompass tens if not hundreds of genes. However, recent studies have exploited mouse and rat resources that allow mapping of phenotypes to narrow intervals, encompassing only a few genes. Here, we review these studies, showcase the rodent resources they have used and highlight the insights into neurobehavioural traits provided to date. We discuss what we see as the biggest challenge in the field - translating QTLs into biological knowledge by experimentally validating and functionally characterizing candidate genes - and propose that the CRISPR/Cas genome-editing system holds the key to overcoming this obstacle. Finally, we challenge traditional views on inbred versus outbred resources in the light of recent resource and technology developments.

Orbitofrontal Neuroadaptations and Cross-Species Synaptic Biomarkers in Heavy-Drinking Macaques

Cognitive impairments, uncontrolled drinking, and neuropathological cortical changes characterize alcohol use disorder. Dysfunction of the orbitofrontal cortex (OFC), a critical cortical subregion that controls learning, decision-making, and prediction of reward outcomes, contributes to executive cognitive function deficits in alcoholic individuals. Electrophysiological and quantitative synaptomics techniques were used to test the hypothesis that heavy drinking produces neuroadaptations in the macaque OFC. Integrative bioinformatics and reverse genetic approaches were used to identify and validate synaptic proteins with novel links to heavy drinking in BXD mice. In drinking monkeys, evoked firing of OFC pyramidal neurons was reduced, whereas the amplitude and frequency of postsynaptic currents were enhanced compared with controls. Bath application of alcohol reduced evoked firing in neurons from control monkeys, but not drinking monkeys. Profiling of the OFC synaptome identified alcohol-sensitive proteins that control glutamate release (e.g., SV2A, synaptogyrin-1) and postsynaptic signaling (e.g., GluA1, PRRT2) with no changes in synaptic GABAergic proteins. Western blot analysis confirmed the increase in GluA1 expression in drinking monkeys. An exploratory analysis of the OFC synaptome found cross-species genetic links to alcohol intake in discrete proteins (e.g., C2CD2L, DIRAS2) that discriminated between low- and heavy-drinking monkeys. Validation studies revealed that BXD mouse strains with the D allele at the C2cd2l interval drank less alcohol than B allele strains. Thus, by profiling of the OFC synaptome, we identified changes in proteins controlling glutamate release and postsynaptic signaling and discovered several proteins related to heavy drinking that have potential as novel targets for treating alcohol use disorder.SIGNIFICANCE STATEMENT Clinical research identified cognitive deficits in alcoholic individuals as a risk factor for relapse, and alcoholic individuals display deficits on cognitive tasks that are dependent upon the orbitofrontal cortex (OFC). To identify neurobiological mechanisms that underpin OFC dysfunction, this study used electrophysiology and integrative synaptomics in a translational nonhuman primate model of heavy alcohol consumption. We found adaptations in synaptic proteins that control glutamatergic signaling in chronically drinking monkeys. Our functional genomic exploratory analyses identified proteins with genetic links to alcohol and cocaine intake across mice, monkeys, and humans. Future work is necessary to determine whether targeting these novel targets reduces excessive and harmful levels of alcohol drinking.

Resources for Systems Genetics

A key characteristic of systems genetics is its reliance on populations that vary to a greater or lesser degree in genetic complexity-from highly admixed populations such as the Collaborative Cross and Diversity Outcross to relatively simple crosses such as sets of consomic strains and reduced complexity crosses. This protocol is intended to help investigators make more informed decisions about choices of resources given different types of questions. We consider factors such as costs, availability, and ease of breeding for common scenarios. In general, we recommend using complementary resources and minimizing depth of resampling of any given genome or strain.

LRRK2 Pathways Leading to Neurodegeneration

Mutations in LRRK2 are associated with inherited Parkinson's disease (PD) in a large number of families, and the genetic locus containing the LRRK2 gene contains a risk factor for sporadic PD. The LRRK2 protein contains several domains that suggest a role in cellular signaling, including a kinase domain. It is also clear that LRRK2 interacts, either physically or genetically, with several other important proteins implicated in PD, suggesting that LRRK2 may be a central player in the pathways that underlie parkinsonism. As such, LRRK2 has been proposed to be a plausible target for therapeutic intervention, with kinase inhibition being pursued most actively. However, there are still several fundamental aspects of LRRK2 biology and function that remain unresolved at this time. This review will focus on the key questions of normal function of LRRK2 and how this might be related to the pathophysiology of PD.

The Genetic Architecture of Murine Glutathione Transferases

Glutathione S-transferase (GST) genes play a protective role against oxidative stress and may influence disease risk and drug pharmacokinetics. In this study, massive multiscalar trait profiling across a large population of mice derived from a cross between C57BL/6J (B6) and DBA2/J (D2)—the BXD family—was combined with linkage and bioinformatic analyses to characterize mechanisms controlling GST expression and to identify downstream consequences of this variation. Similar to humans, mice show a wide range in expression of GST family members. Variation in the expression of Gsta4, Gstt2, Gstz1, Gsto1, and Mgst3 is modulated by local expression QTLs (eQTLs) in several tissues. Higher expression of Gsto1 in brain and liver of BXD strains is strongly associated (P < 0.01) with inheritance of the B6 parental allele whereas higher expression of Gsta4 and Mgst3 in brain and liver, and Gstt2 and Gstz1 in brain is strongly associated with inheritance of the D2 parental allele. Allele-specific assays confirmed that expression of Gsto1, Gsta4, and Mgst3 are modulated by sequence variants within or near each gene locus. We exploited this endogenous variation to identify coexpression networks and downstream targets in mouse and human. Through a combined systems genetics approach, we provide new insight into the biological role of naturally occurring variants in GST genes.

RNA sequencing of transcriptomes in human brain regions: protein-coding and non-coding RNAs, isoforms and alleles

Background

We used RNA sequencing to analyze transcript profiles of ten autopsy brain regions from ten subjects. RNA sequencing techniques were designed to detect both coding and non-coding RNA, splice isoform composition, and allelic expression. Brain regions were selected from five subjects with a documented history of smoking and five non-smokers. Paired-end RNA sequencing was performed on SOLiD instruments to a depth of >40 million reads, using linearly amplified, ribosomally depleted RNA. Sequencing libraries were prepared with both poly-dT and random hexamer primers to detect all RNA classes, including long non-coding (lncRNA), intronic and intergenic transcripts, and transcripts lacking poly-A tails, providing additional data not previously available. The study was designed to generate a database of the complete transcriptomes in brain region for gene network analyses and discovery of regulatory variants.

Results

Of 20,318 protein coding and 18,080 lncRNA genes annotated from GENCODE and lncipedia, 12 thousand protein coding and 2 thousand lncRNA transcripts were detectable at a conservative threshold. Of the aligned reads, 52 % were exonic, 34 % intronic and 14 % intergenic. A majority of protein coding genes (65 %) was expressed in all regions, whereas ncRNAs displayed a more restricted distribution. Profiles of RNA isoforms varied across brain regions and subjects at multiple gene loci, with neurexin 3 (NRXN3) a prominent example. Allelic RNA ratios deviating from unity were identified in > 400 genes, detectable in both protein-coding and non-coding genes, indicating the presence of cis-acting regulatory variants. Mathematical modeling was used to identify RNAs stably expressed in all brain regions (serving as potential markers for normalizing expression levels), linked to basic cellular functions. An initial analysis of differential expression analysis between smokers and nonsmokers implicated a number of genes, several previously associated with nicotine exposure.

Conclusions

RNA sequencing identifies distinct and consistent differences in gene expression between brain regions, with non-coding RNA displaying greater diversity between brain regions than mRNAs. Numerous RNAs exhibit robust allele selective expression, proving a means for discovery of cis-acting regulatory factors with potential clinical relevance.

Electronic supplementary material

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

Segregation of a spontaneous Klrd1 (CD94) mutation in DBA/2 mouse substrains

Current model DBA/2J (D2J) mice lack CD94 expression due to a deletion spanning the last coding exon of the Klrd1 gene that occurred in the mid- to late 1980s. In contrast, DBA/2JRj (D2Rj) mice, crosses derived from DBA/2J before 1984, and C57BL/6J (B6) mice lack the deletion and have normal CD94 expression. For example, BXD lines (BXD1–32) generated in the 1970s by crossing B6 and D2J do not segregate for the exonic deletion and have high expression, whereas BXD lines 33 and greater were generated after 1990 are segregating for the deletion and have highly variable Klrd1 expression. We performed quantitative trait locus analysis of Klrd1 expression by using BXD lines with different generation times and found that the expression difference in Klrd1 in the later BXD set is driven by a strong cis-acting expression quantitative trait locus. Although the Klrd1/CD94 locus is essential for mousepox resistance, the genetic variation among D2 substrains and the later set of BXD strains is not associated with susceptibility to the Influenza A virus PR8 strain. Substrains with nearly identical genetic backgrounds that are segregating functional variants such as the Klrd1 deletion are useful genetic tools to investigate biological function.