Gene expression profiles in rat brain disclose CNS signature genes and regional patterns of functional specialisation

Gene expression profiles in Malpighian tubules of the vector leafhopper Psammotettix striatus (L.) revealed regional functional diversity and heterogeneity

Background

Many leafhoppers are known as pests and disease vectors of economically important plants. Previous studies of the physiological functions of vector leafhoppers have mainly focused on the salivary glands and the alimentary tract that are deemed to be associated with digestion, host defense and phytoplasma and/or virus transmission. By contrast, the significance of Malpighian tubules (MTs) is less studied. To clarify the physiological function of MTs of the vector leafhopper Psammotettix striatus that transmits phytoplasma triggering the wheat blue dwarf disease, we performed a transcriptome study on P. striatus MTs and compared gene expression profiles among different anatomical regions in the tubules (i.e., MT1+2, the anterior segment together with the sub-anterior segment; MT3, the inflated segment; and MT4, the distal segment).

Results

Transcriptome of P. striatus MTs generate a total of 42,815 high-quality unigenes, among which highly expressed unigenes are mainly involved in organic solute transport, detoxification and immunity in addition to osmoregulation. Region-specific comparative analyses reveal that all these MT regions have functions in osmoregulation, organic solute transport and detoxification, but each region targets different substrates. Differential expression and regional enrichment of immunity-related effector activities and molecules involved in phagocytosis and the biosynthesis of antimicrobial peptides among different regions indicate that MT1+2 and MT4 have the ability to eliminate the invading pathogens. However, in MT3 which secrets brochosomes to the integument and eggs as physical barriers, disulfide-isomerase, acidic ribosomal protein P and many other unigenes were highly expressed, which can be attractive candidate genes for future studies of the biosynthesis and the origin of brochosomes.

Conclusions

Psammotettix striatus MTs perform multiple physiological functions as versatile organs than just excretory organs with osmoregulatory function. Heterogeneity of physiological functions among different MT regions is related to organic solute transport, detoxification, immunity and brochosome biosynthesis in addition to osmoregulation, and each region targets different substrates. These functions may be helpful for P. striatus to resist pathogens from habitats and to utilize a wider range of host plants, which may assist the transmission and spread of phytoplasmas. The results provide potential molecular targets for the exploit of chemical and/or gene-silencing insecticides.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08300-6.

The Expression of Hypoxia-Induced Gene 1 (Higd1a) in the Central Nervous System of Male and Female Rats Differs According to Age

HIGD1A (hypoxia-induced gene domain protein-1a), a mitochondrial inner membrane protein present in various cell types, has been mainly associated with anti-apoptotic processes in response to stressors. Our previous findings have shown that Higd1a mRNA is widely expressed across the central nervous system (CNS), exhibiting an increasing expression in the spinal cord from postnatal day 1 (P1) to 15 (P15) and changes in the distribution pattern from P1 to P90. During the first weeks of postnatal life, the great plasticity of the CNS is accompanied by cell death/survival decisions. So we first describe HIGD1A expression throughout the brain during early postnatal life in female and male pups. Secondly, based on the fact that in some areas this process is influenced by the sex of individuals, we explore HIGD1A expression in the sexual dimorphic nucleus (SDN) of the medial preoptic area, a region that is several folds larger in male than in female rats, partly due to sex differences in the process of apoptosis during this period. Immunohistochemical analysis revealed that HIGD1A is widely but unevenly expressed throughout the brain. Quantitative Western blot analysis of the parietal cortex, diencephalon, and spinal cord from both sexes at P1, P5, P8, and P15 showed that the expression of this protein is predominantly high and changes with age but not sex. Similarly, in the sexual dimorphic nucleus, the expression of HIGD1A varied according to age, but we were not able to detect significant differences in its expression according to sex. Altogether, these results suggest that HIGD1A protein is expressed in several areas of the central nervous system following a pattern that quantitatively changes with age but does not seem to change according to sex.

Genetic and neuroendocrine regulation of the postpartum brain

Changes in expression of hundreds of genes occur during the production and function of the maternal brain that support a wide range of processes. In this review, we synthesize findings from four microarray studies of different maternal brain regions and identify a core group of 700 maternal genes that show significant expression changes across multiple regions. With those maternal genes, we provide new insights into reward-related pathways (maternal bonding), postpartum depression, social behaviors, mental health disorders, and nervous system plasticity/developmental events. We also integrate the new genes into well-studied maternal signaling pathways, including those for prolactin, oxytocin/vasopressin, endogenous opioids, and steroid receptors (estradiol, progesterone, cortisol). A newer transcriptional regulation model for the maternal brain is provided that incorporates recent work on maternal microRNAs. We also compare the top 700 genes with other maternal gene expression studies. Together, we highlight new genes and new directions for studies on the postpartum brain.

MicroRNA expression is altered in lateral septum across reproductive stages

MicroRNAs (miRNAs) inhibit RNA targets and may contribute to postpartum central nervous system (CNS) gene expression changes, although this has never been tested. In the present study, we directly evaluated miRNA levels using RNA sequencing during reproduction in female mice in the lateral septum (LS). We found the reliable and robust changes of miRNAs away from the virgin stage at the three other stages, namely pregnant, day 1 postpartum, and day 8 postpartum. For a given miRNA that was significantly different from the virgin condition in more than one group, the direction of change was always the same. Overall, we identified 32 upregulated miRNAs and 25 downregulated miRNAs that were consistently different from the virgin state. 'Arm switching' occurs for miR-433-3 and miR-7b. Unexpectedly, a third of upregulated miRNAs (relative to virgin) were highly localized within the 12qF1 region of chromosome 12 that includes the Dlk1-Dio3 gene cluster implicated in stem cell and neuronal differentiation. Over 1500 genes were targeted by multiple upregulated miRNAs with about 100 genes targeted by five or more miRNAs. Over 1000 genes were targeted by multiple downregulated miRNAs with about 50 genes targeted by five or more miRNAs. Half of the target genes were regulated by up and downregulated miRNAs, indicating homeostatic regulation. Transcriptional regulation was the most enriched pathway for genes linked to up or down regulated miRNAs. Other enriched pathways included protein kinase activity (e.g., MAP kinase), CNS development, axon guidance, neurotrophin signaling, neuron development/differentiation, and neurogenesis. Previously published postpartum LS gene expression changes were enrichment for LS miRNA targets, as expected. Surprisingly, postpartum gene expression changes from other regions were also enriched against LS miRNA targets, suggesting a core group of miRNAs may act across the CNS during reproduction. Together, we directly examine miRNAs and find significant alterations in the postpartum brain.

The Effects of Different Repetitive Transcranial Magnetic Stimulation (rTMS) Protocols on Cortical Gene Expression in a Rat Model of Cerebral Ischemic-Reperfusion Injury

Although repetitive Transcranial Magnetic Stimulation (rTMS) in treatment of stroke in humans has been explored over the past decade the data remain controversial in terms of optimal stimulation parameters and the mechanisms of rTMS long-term effects. This study aimed to explore the potential of different rTMS protocols to induce changes in gene expression in rat cortices after acute ischemic-reperfusion brain injury. The stroke was induced by middle cerebral artery occlusion (MCAO) with subsequent reperfusion. Changes in the expression of 96 genes were examined using low-density expression arrays after MCAO alone and after MCAO combined with 1Hz, 5Hz, continuous (cTBS) and intermittent (iTBS) theta-burst rTMS. rTMS over the lesioned hemisphere was given for two weeks (with a 2-day pause) in a single daily session and a total of 2400 pulses. MCAO alone induced significant upregulation in the expression of 44 genes and downregulation in 10. Two weeks of iTBS induced significant increase in the expression of 52 genes. There were no downregulated genes. 1Hz and 5Hz had no significant effects on gene expression, while cTBS effects were negligible. Upregulated genes included those involved in angiogenesis, inflammation, injury response and cellular repair, structural remodeling, neuroprotection, neurotransmission and neuronal plasticity. The results show that long-term rTMS in acute ischemic-reperfusion brain injury induces complex changes in gene expression that span multiple pathways, which generally promote the recovery. They also demonstrate that induced changes primarily depend on the rTMS frequency (1Hz and 5Hz vs. iTBS) and pattern (cTBS vs. iTBS). The results further underlines the premise that one of the benefits of rTMS application in stroke may be to prime the brain, enhancing its potential to cope with the injury and to rewire. This could further augment its potential to favorably respond to rehabilitation, and to restore some of the loss functions.

Postnatal neurodevelopmental expression and glutamate-dependent regulation of the ZNF804A rodent homologue

The zinc finger protein ZNF804A rs1344706 variant is a replicated genome-wide significant risk variant for schizophrenia and bipolar disorder. While its association with altered brain structure and cognition in patients and healthy risk allele carriers is well documented, the characteristics and function of the gene in the brain remains poorly understood. Here, we used in situ hybridization to determine mRNA expression levels of the ZNF804A rodent homologue, Zfp804a, across multiple postnatal neurodevelopmental time points in the rat brain. We found changes in Zfp804a expression in the rat hippocampus, frontal cortex, and thalamus across postnatal neurodevelopment. Zfp804a mRNA peaked at postnatal day (P) 21 in hippocampal CA1 and DG regions and was highest in the lower cortical layers of frontal cortex at P1, possibly highlighting a role in developmental migration. Using immunofluorescence, we found that Zfp804a mRNA and ZFP804A co-localized with neurons and not astrocytes. In primary cultured cortical neurons, we found that Zfp804a expression was significantly increased when neurons were exposed to glutamate [20μM], but this increase was blocked by the N-methyl-d-aspartate receptor (NMDAR) antagonist MK-801. Expression of Comt, Pde4b, and Drd2, genes previously shown to be regulated by ZNF804A overexpression, was also significantly changed in an NMDA-dependent manner. Our results describe, for the first time, the unique postnatal neurodevelopmental expression of Zfp804a in the rodent brain and demonstrate that glutamate potentially plays an important role in the regulation of this psychiatric susceptibility gene. These are critical steps toward understanding the biological function of ZNF804A in the mammalian brain.

Integrated annotation and analysis of in situ hybridization images using the ImAnno system: application to the ear and sensory organs of the fetal mouse

An in situ hybridization (ISH) study was performed on 2000 murine genes representing around 10% of the protein-coding genes present in the mouse genome using data generated by the EURExpress consortium. This study was carried out in 25 tissues of late gestation embryos (E14.5), with a special emphasis on the developing ear and on five distinct developing sensory organs, including the cochlea, the vestibular receptors, the sensory retina, the olfactory organ, and the vibrissae follicles. The results obtained from an analysis of more than 11,000 micrographs have been integrated in a newly developed knowledgebase, called ImAnno. In addition to managing the multilevel micrograph annotations performed by human experts, ImAnno provides public access to various integrated databases and tools. Thus, it facilitates the analysis of complex ISH gene expression patterns, as well as functional annotation and interaction of gene sets. It also provides direct links to human pathways and diseases. Hierarchical clustering of expression patterns in the 25 tissues revealed three main branches corresponding to tissues with common functions and/or embryonic origins. To illustrate the integrative power of ImAnno, we explored the expression, function and disease traits of the sensory epithelia of the five presumptive sensory organs. The study identified 623 genes (out of 2000) concomitantly expressed in the five embryonic epithelia, among which many (∼12%) were involved in human disorders. Finally, various multilevel interaction networks were characterized, highlighting differential functional enrichments of directly or indirectly interacting genes. These analyses exemplify an under-represention of "sensory" functions in the sensory gene set suggests that E14.5 is a pivotal stage between the developmental stage and the functional phase that will be fully reached only after birth.

Quantitative proteomic profiling reveals novel region-specific markers in the adult mouse brain

Despite major advances in neuroscience, a comprehensive understanding of the structural and functional components of the adult brain compartments remains to be fully elucidated at a quantitative molecular level. Indeed, over half of the soluble- and membrane-annotated proteins are currently unmapped within online digital brain atlases. In this study, two complementary approaches were used to assess the unique repertoire of proteins enriched within select regions of the adult mouse CNS, including the brain stem, cerebellum, and remaining brain hemispheres. Of the 1200 proteins visualized by 2D-DIGE, approximately 150 (including cytosolic and membrane proteins) were found to exhibit statistically significant changes in relative abundance thus representing putative region-specific brain markers. In addition to using a high-precision (18) O-labeling strategy for the quantitative LC-MS/MS mapping of membrane proteins isolated from myelin-enriched fractions, we have identified over 1000 proteins that have yet to be described in any other mammalian myelin proteome. A comparison of our myelin proteome was made to an existing transcriptome database containing mRNA abundance profiles during oligodendrocyte differentiation and has confirmed statistically significant abundance changes for ∼500 of these newly mapped proteins, thus revealing new roles in oligodendrocyte and myelin biology. These data offer a resource for the neuroscience community studying the molecular basis for specialized neuronal activities in the CNS and myelin-related disorders. The MS proteomics data associated with this manuscript have been deposited to the ProteomeXchange Consortium with the dataset identifier PXD000327 (http://proteomecentral.proteomexchange.org/dataset/PXD000327).

Model of tryptophan metabolism, readily scalable using tissue-specific gene expression data

Tryptophan is utilized in various metabolic routes including protein synthesis, serotonin, and melatonin synthesis and the kynurenine pathway. Perturbations in these pathways have been associated with neurodegenerative diseases and cancer. Here we present a comprehensive kinetic model of the complex network of human tryptophan metabolism based upon existing kinetic data for all enzymatic conversions and transporters. By integrating tissue-specific expression data, modeling tryptophan metabolism in liver and brain returned intermediate metabolite concentrations in the physiological range. Sensitivity and metabolic control analyses identified expected key enzymes to govern fluxes in the branches of the network. Combining tissue-specific models revealed a considerable impact of the kynurenine pathway in liver on the concentrations of neuroactive derivatives in the brain. Moreover, using expression data from a cancer study predicted metabolite changes that resembled the experimental observations. We conclude that the combination of the kinetic model with expression data represents a powerful diagnostic tool to predict alterations in tryptophan metabolism. The model is readily scalable to include more tissues, thereby enabling assessment of organismal tryptophan metabolism in health and disease.

LOC689986, a unique gene showing specific expression in restricted areas of the rodent neocortex

Background

The neocortex is a highly specialised and complex brain structure, involved in numerous tasks, ranging from processing and interpretation of somatosensory information, to control of motor functions. The normal function linked to distinct neocortical areas might involve control of highly specific gene expression, and in order to identify such regionally enriched genes, we previously analysed the global gene expression in three different cortical regions (frontomedial, temporal and occipital cortex) from the adult rat brain. We identified distinct sets of differentially expressed genes. One of these genes, namely the hypothetical protein LOC689986 (LOC689986), was of particular interest, due to an almost exclusive expression in the temporal cortex.

Results

Detailed analysis of LOC689986 in the adult rat brain confirmed the expression in confined areas of parieto-temporal cortex, and revealed highly specific expression in layer 4 of the somatosensory cortex, with sharp borders towards the neighbouring motor cortex. In addition, LOC689986 was found to be translated in vivo, and was detected in the somatosensory cortex and in the Purkinje cells of the cerebellar cortex. The protein was present in neuronal dendrites and also in astrocyte cells. Finally, this unique gene is apparently specific for, and highly conserved in, the vertebrate lineage.

Conclusions

In this study, we have partially characterised the highly conserved LOC689986 gene, which is specific to the vertebrate linage. The gene displays a distinct pattern of expression in layer 4 of the somatosensory cortex, and areas of the parieto-temporal cortex in rodents.

Distinct functional brain regional integration of Casp3, Ascl1 and S100a6 gene expression in spatial memory

Evaluating the brain structural expression of defined genes involved in basic biological processes of neurogenesis, apoptosis or neural plasticity may facilitate the understanding of genetic mechanisms underlying spatial memory. The aim of the present study was to compare Ascl1, Casp3 and S100a6 gene expression in the hippocampus, prefrontal cortex and cerebellum of adult rats in water maze spatial memory performance. After four days training, the mean platform time (<10s) was evidence of stable long-term spatial memory formation. Real time PCR analysis revealed a positive inter-structural correlation for S100a6/Casp gene expression between the prefrontal cortex and the cerebellum but a negative correlation for S100a6/Ascl1 transcribed genes between the prefrontal cortex and hippocampus during swimming in the active controls. However, during spatial memory performance there was only one inter-structural correlation between the prefrontal cortex and cerebellum with respect to Casp3 expression, though there were intra-structural correlations between Casp3/Ascl1 transcriptions within the prefrontal cortex and hippocampus as well as between Ascl1/S100a6 in the cerebellum. In active learners versus naive controls, the transcrption of all genes was augmented in the prefrontal cortex but Casp3 and Ascl1 were also elevated in hippocampus whilst only S100a6 increased in the cerebellum. The findings endorsed the role of the hippocampus in memory acquisition in addition to an integrative relationship with the prefrontal cortex and cerebellum. This structural and molecular configuration is important for creation of novel neural circuitry for consolidation and reconsolidation of memory trace with an involvement of coupled processes of neurogenesis, apoptosis or neural plasticity.

Regional features of the expression of genes involved in neurogenesis and apoptosis in the brain of adult rats

The expression of mRNA of genes involved in neurogenesis and apoptosis (Apaf1, Ascl1, Bax, Bcl2, Casp3, Casp8, Casp9, Dffb, Myh10, Naip2, Napa, Notch2, Numb, Pura, S100a6, Tnf) in the prefrontal cortex, hippocampus, and cerebellum was studied in adult rats. The content of mRNA of these genes (except Apaf1) was several-fold higher in the cerebellum than in the hippocampus and brain cortex. In the hippocampus, the expression of Apaf1 was significantly lower than in the prefrontal cortex, while the expression of Ascl1, Pura, S100b, and Tnf was higher. Regional differences in the direction, strength, and numbers of significant correlations between the expression of the studied genes were detected. Documented differences in gene expression were regarded as validation of the structural and functional cooperation of neurogenesis and apoptosis at the molecular genetic level.

Developmental traumatic brain injury decreased brain derived neurotrophic factor expression late after injury

Pediatric traumatic brain injury (TBI) is a major cause of acquired cognitive dysfunction in children. Hippocampal Brain Derived Neurotrophic Factor (BDNF) is important for normal cognition. Little is known about the effects of TBI on BDNF levels in the developing hippocampus. We used controlled cortical impact (CCI) in the 17 day old rat pup to test the hypothesis that CCI would first increase rat hippocampal BDNF mRNA/protein levels relative to SHAM and Naïve rats by post injury day (PID) 2 and then decrease BDNF mRNA/protein by PID14. Relative to SHAM, CCI did not change BDNF mRNA/protein levels in the injured hippocampus in the first 2 days after injury but did decrease BDNF protein at PID14. Surprisingly, BDNF mRNA decreased at PID 1, 3, 7 and 14, and BDNF protein decreased at PID 2, in SHAM and CCI hippocampi relative to Naïve. In conclusion, TBI decreased BDNF protein in the injured rat pup hippocampus 14 days after injury. BDNF mRNA levels decreased in both CCI and SHAM hippocampi relative to Naïve, suggesting that certain aspects of the experimental paradigm (such as craniotomy, anesthesia, and/or maternal separation) may decrease the expression of BDNF in the developing hippocampus. While BDNF is important for normal cognition, no inferences can be made regarding the cognitive impact of any of these factors. Such findings, however, suggest that meticulous attention to the experimental paradigm, and possible inclusion of a Naïve group, is warranted in studies of BDNF expression in the developing brain after TBI.

Gene-based analysis of regionally enriched cortical genes in GWAS data sets of cognitive traits and psychiatric disorders

BACKGROUND

Despite its estimated high heritability, the genetic architecture leading to differences in cognitive performance remains poorly understood. Different cortical regions play important roles in normal cognitive functioning and impairment. Recently, we reported on sets of regionally enriched genes in three different cortical areas (frontomedial, temporal and occipital cortices) of the adult rat brain. It has been suggested that genes preferentially, or specifically, expressed in one region or organ reflect functional specialisation. Employing a gene-based approach to the analysis, we used the regionally enriched cortical genes to mine a genome-wide association study (GWAS) of the Norwegian Cognitive NeuroGenetics (NCNG) sample of healthy adults for association to nine psychometric tests measures. In addition, we explored GWAS data sets for the serious psychiatric disorders schizophrenia (SCZ) (n = 3 samples) and bipolar affective disorder (BP) (n = 3 samples), to which cognitive impairment is linked.

PRINCIPAL FINDINGS

At the single gene level, the temporal cortex enriched gene RAR-related orphan receptor B (RORB) showed the strongest overall association, namely to a test of verbal intelligence (Vocabulary, P = 7.7E-04). We also applied gene set enrichment analysis (GSEA) to test the candidate genes, as gene sets, for enrichment of association signal in the NCNG GWAS and in GWASs of BP and of SCZ. We found that genes differentially expressed in the temporal cortex showed a significant enrichment of association signal in a test measure of non-verbal intelligence (Reasoning) in the NCNG sample.

CONCLUSION

Our gene-based approach suggests that RORB could be involved in verbal intelligence differences, while the genes enriched in the temporal cortex might be important to intellectual functions as measured by a test of reasoning in the healthy population. These findings warrant further replication in independent samples on cognitive traits.

Large-scale neuroinformatics for in situ hybridization data in the mouse brain

Large-scale databases of the brain are providing content to the neuroscience community through molecular, cellular, functional, and connectomic data. Organization, presentation, and maintenance requirements are substantial given the complexity, diverse modalities, resolution, and scale. In addition to microarrays, magnetic resonance imaging, and RNA sequencing, several in situ hybridization databases have been constructed due to their value in spatially localizing cellular expression. Scalable techniques for processing and presenting these data for maximum utility in viewing and analysis are key for end user value. We describe methods and use cases for the Allen Brain Atlas resources of the adult and developing mouse.

Adaptation of the CHARM DNA methylation platform for the rat genome reveals novel brain region-specific differences

Comprehensive High-throughput Arrays for Relative Methylation (CHARM) was recently developed as an experimental platform and analytic approach to assess DNA methylation (DNAm) at a genome-wide level. Its initial implementation was for human and mouse. We adapted it for rat and sought to examine DNAm differences across tissues and brain regions in this model organism. We extracted DNA from liver, spleen, and three brain regions: cortex, hippocampus, and hypothalamus from adult Sprague Dawley rats. DNA was digested with McrBC, and the resulting methyl-depleted fraction was hybridized to the rat CHARM array along with a mock-treated fraction. Differentially methylated regions (DMRs) between tissue types were detected using normalized methylation log-ratios. In validating 24 of the most significant DMRs by bisulfite pyrosequencing, we detected large mean differences in DNAm, ranging from 33-59%, among the most significant DMRs in the across-tissue comparisons. The comparable figures for the hippocampus vs. hypothalamus DMRs were 14-40%, for the cortex vs. hippocampus DMRs, 12-29%, and for the cortex vs. hypothalamus DMRs, 5-35%, with a correlation of r(2) = 0.92 between the methylation differences in 24 DMRs predicted by CHARM and those validated by bisulfite pyrosequencing. Our adaptation of the CHARM array for the rat genome yielded highly robust results that demonstrate the value of this method in detecting substantial DNAm differences between tissues and across different brain regions. This platform should prove valuable in future studies aimed at examining DNAm differences in particular brain regions of rats exposed to environmental stimuli with potential epigenetic consequences.

Decreased proliferation of adult hippocampal stem cells during cocaine withdrawal: possible role of the cell fate regulator FADD

The current study uses an extended access rat model of cocaine self-administration (5-h session per day, 14 days), which elicits several features manifested during the transition to human addiction, to study the neural adaptations associated with cocaine withdrawal. Given that the hippocampus is thought to have an important role in maintaining addictive behavior and appears to be especially relevant to mechanisms associated with withdrawal, this study attempted to understand how extended access to cocaine impacts the hippocampus at the cellular and molecular levels, and how these alterations change over the course of withdrawal (1, 14, and 28 days). Therefore, at the cellular level, we examined the effects of cocaine withdrawal on cell proliferation (Ki-67+ and NeuroD+ cells) in the DG. At the molecular level, we employed a 'discovery' approach with gene expression profiling in the DG to uncover novel molecules possibly implicated in the neural adaptations that take place during cocaine withdrawal. Our results suggest that decreased hippocampal cell proliferation might participate in the adaptations associated with drug removal and identifies 14 days as a critical time-point of cocaine withdrawal. At the 14-day time-point, gene expression profiling of the DG revealed the dysregulation of several genes associated with cell fate regulation, highlighting two new neurobiological correlates (Ascl-1 and Dnmt3b) that accompany cessation of drug exposure. Moreover, the results point to Fas-Associated protein with Death Domain (FADD), a molecular marker previously associated with the propensity to substance abuse and cocaine sensitization, as a key cell fate regulator during cocaine withdrawal. Identifying molecules that may have a role in the restructuring of the hippocampus following substance abuse provides a better understanding of the adaptations associated with cocaine withdrawal and identifies novel targets for therapeutic intervention.

Brain expression quantitative trait locus mapping informs genetic studies of psychiatric diseases

Genome-wide association study (GWAS) can be used to identify genes that increase the risk of psychiatric diseases. However, much of the disease heritability is still unexplained, suggesting that there are genes to be discovered. Functional annotation of the genetic variants may increase the power of GWAS to identify disease genes, by providing prior information that can be used in Bayesian analysis or in reducing the number of tests. Expression quantitative trait loci (eQTLs) are genomic loci that regulate gene expression. Genetic mapping of eQTLs can help reveal novel functional effects of thousands of single nucleotide polymorphisms (SNPs). The present review mainly focused on the current knowledge on brain eQTL mapping, and discussed some major methodological issues and their possible solutions. The frequently ignored problems of batch effects, covariates, and multiple testing were emphasized, since they can lead to false positives and false negatives. The future application of eQTL data in GWAS analysis was also discussed.

Early and sustained increase in the expression of hippocampal IGF-1, but not EPO, in a developmental rodent model of traumatic brain injury

Pediatric traumatic brain injury (pTBI) is the leading cause of traumatic death and disability in children in the United States. Impaired learning and memory in these young survivors imposes a heavy toll on society. In adult TBI (aTBI) models, cognitive outcome improved after administration of erythropoietin (EPO) or insulin-like growth factor-1 (IGF-1). Little is known about the production of these agents in the hippocampus, a brain region critical for learning and memory, after pTBI. Our objective was to describe hippocampal expression of EPO and IGF-1, together with their receptors (EPOR and IGF-1R, respectively), over time after pTBI in 17-day-old rats. We used the controlled cortical impact (CCI) model and measured hippocampal mRNA levels of EPO, IGF-1, EPOR, IGF-1R, and markers of caspase-dependent apoptosis (bcl2, bax, and p53) at post-injury days (PID) 1, 2, 3, 7, and 14. CCI rats performed poorly on Morris water maze testing of spatial working memory, a hippocampally-based cognitive function. Apoptotic markers were present early and persisted for the duration of the study. EPO in our pTBI model increased much later (PID7) than in aTBI models (12 h), while EPOR and IGF-1 increased at PID1 and PID2, respectively, similar to data from aTBI models. Our data indicate that EPO expression showed a delayed upregulation post-pTBI, while EPOR increased early. We speculate that administration of EPO in the first 1-2 days after pTBI would increase hippocampal neuronal survival and function.

Gene expression in the rat brain: high similarity but unique differences between frontomedial-, temporal- and occipital cortex

BACKGROUND

The six-layered neocortex of the mammalian brain may appear largely homologous, but is in reality a modular structure of anatomically and functionally distinct areas. However, global gene expression seems to be almost identical across the cerebral cortex and only a few genes have so far been reported to show regional enrichment in specific cortical areas.

RESULTS

In the present study on adult rat brain, we have corroborated the strikingly similar gene expression among cortical areas. However, differential expression analysis has allowed for the identification of 30, 24 and 11 genes enriched in frontomedial -, temporal- or occipital cortex, respectively. A large proportion of these 65 genes appear to be involved in signal transduction, including the ion channel Fxyd6, the neuropeptide Grp and the nuclear receptor Rorb. We also find that the majority of these genes display increased expression levels around birth and show distinct preferences for certain cortical layers and cell types in rodents.

CONCLUSIONS

Since specific patterns of expression often are linked to equally specialised biological functions, we propose that these cortex sub-region enriched genes are important for proper functioning of the cortical regions in question.

Upregulation of immunoglobulin-related genes in cortical sections from multiple sclerosis patients

Multiple sclerosis (MS) is a demyelinating disease of the central nervous system (CNS). Microarray-based global gene expression profiling is a promising method, used to study potential genes involved in the pathogenesis of the disease. In the present study, we have examined global gene expression in normal-appearing gray matter and gray matter lesions from the cortex of MS patients, and compared them with cortical gray matter samples from controls. We observed a massive upregulation of immunoglobulin (Ig)-related genes in cortical sections of MS patients. Using immunohistochemistry, the activation of Ig genes seems to occur within plasma cells in the meninges. As synthesis of oligoclonal IgGs has been hypothesized to be caused by the activation of Epstein-Barr virus (EBV)-infected B-cells, we screened the brain samples for the presence of EBV by real-time quantitative polymerase chain reaction (qPCR) and immunohistochemistry, but no evidence of active or latent EBV infection was detected. This study demonstrates that genes involved in the synthesis of Igs are upregulated in MS patients and that this activation is caused by a small number of meningeal plasma cells that are not infected by EBV. The findings indicate that the Ig-producing B-cells found in the cerebrospinal fluid (CSF) of MS patients could have meningeal origin.

Mammalian epoxide hydrolases in xenobiotic metabolism and signalling

Epoxide hydrolases catalyse the hydrolysis of electrophilic--and therefore potentially genotoxic--epoxides to the corresponding less reactive vicinal diols, which explains the classification of epoxide hydrolases as typical detoxifying enzymes. The best example is mammalian microsomal epoxide hydrolase (mEH)-an enzyme prone to detoxification-due to a high expression level in the liver, a broad substrate selectivity, as well as inducibility by foreign compounds. The mEH is capable of inactivating a large number of structurally different, highly reactive epoxides and hence is an important part of the enzymatic defence of our organism against adverse effects of foreign compounds. Furthermore, evidence is accumulating that mammalian epoxide hydrolases play physiological roles other than detoxification, particularly through involvement in signalling processes. This certainly holds true for soluble epoxide hydrolase (sEH) whose main function seems to be the turnover of lipid derived epoxides, which are signalling lipids with diverse functions in regulatory processes, such as control of blood pressure, inflammatory processes, cell proliferation and nociception. In recent years, the sEH has attracted attention as a promising target for pharmacological inhibition to treat hypertension and possibly other diseases. Recently, new hitherto uncharacterised epoxide hydrolases could be identified in mammals by genome analysis. The expression pattern and substrate selectivity of these new epoxide hydrolases suggests their participation in signalling processes rather than a role in detoxification. Taken together, epoxide hydrolases (1) play a central role in the detoxification of genotoxic epoxides and (2) have an important function in the regulation of physiological processes by the control of signalling molecules with an epoxide structure.

Refined anatomical isolation of functional sleep circuits exhibits distinctive regional and circadian gene transcriptional profiles

Powerful new approaches to study molecular variation in distinct neuronal populations have recently been developed enabling a more precise investigation of the control of neural circuits involved in complex behaviors such as wake and sleep. We applied laser capture microdissection (LCM) to isolate precise brain nuclei from rat CNS at opposing circadian time points associated with wake and sleep. Discrete anatomical and temporal analysis was performed to examine the extent of variation in the transcriptional control associated with both identifiable anatomical nuclei and with light/dark cycle. Precise isolation of specific brain nuclei regulating sleep and arousal, including the LC, SCN, TMN, VTA, and VLPO, demonstrated robust changes in gene expression. Many of these differences were not observed in previous studies where whole brain lysates or gross dissections were used to probe for changes in gene expression. The robust and differential profiles of genomic data obtained from the approaches used herein underscore the requirement for careful anatomical refinement in CNS gene expression studies designed to understand genomic control within behaviorally-linked, but functionally isolated brain nuclei.

An 8.9 Mb 19p13 duplication associated with precocious puberty and a sporadic 3.9 Mb 2q23.3q24.1 deletion containing NR4A2 in mentally retarded members of a family with an intrachromosomal 19p-into-19q between-arm insertion

In a 2 and a half-year-old girl with onset of puberty before the age of 5 months, short stature, hand anomalies and severe mental retardation, an 8.9 Mb interstitial 19p13 duplication containing 215 predicted genes was detected. It was initially assumed that the duplication involved the kisspeptin receptor gene, GPR54, known to stimulate induction of puberty, but more refined duplication mapping excluded this possibility. In an attempt to further understand the genotype-phenotype correlation, global gene expression was measured in skin fibroblasts. The overall expression pattern was quite similar to controls, and only about 25% of the duplicated genes had an expression level that was increased by more than 1.3-fold, with no obvious changes that could explain the precocious puberty. The proband's mother carried a balanced between-arm insertion of the duplicated segment that resembled a pericentric inversion. The same insertion was found in several other family members, including one who had lost a daughter with severe mental retardation and menarche at the age of 10 years. Another close relative was severely mentally retarded, but neither dysmorphic nor microcephalic. His phenotype was initially ascribed to a presumed cryptic chromosome 19 imbalance caused by the 19p-into19q insertion, but subsequent array-CGH detected a 3.9-Mb deletion of 2q23.3q24.1. This novel microdeletion involves seven genes, of which FMNL2, a suggested regulator of Rho-GTPases, and NR4A2, an essential gene for differentiation of dopaminergic neurons, may be critical genes for the proposed 2q23q24 microdeletion syndrome.

A comparison of four clustering methods for brain expression microarray data

Background

DNA microarrays, which determine the expression levels of tens of thousands of genes from a sample, are an important research tool. However, the volume of data they produce can be an obstacle to interpretation of the results. Clustering the genes on the basis of similarity of their expression profiles can simplify the data, and potentially provides an important source of biological inference, but these methods have not been tested systematically on datasets from complex human tissues. In this paper, four clustering methods, CRC, k-means, ISA and memISA, are used upon three brain expression datasets. The results are compared on speed, gene coverage and GO enrichment. The effects of combining the clusters produced by each method are also assessed.

Results

k-means outperforms the other methods, with 100% gene coverage and GO enrichments only slightly exceeded by memISA and ISA. Those two methods produce greater GO enrichments on the datasets used, but at the cost of much lower gene coverage, fewer clusters produced, and speed. The clusters they find are largely different to those produced by k-means. Combining clusters produced by k-means and memISA or ISA leads to increased GO enrichment and number of clusters produced (compared to k-means alone), without negatively impacting gene coverage. memISA can also find potentially disease-related clusters. In two independent dorsolateral prefrontal cortex datasets, it finds three overlapping clusters that are either enriched for genes associated with schizophrenia, genes differentially expressed in schizophrenia, or both. Two of these clusters are enriched for genes of the MAP kinase pathway, suggesting a possible role for this pathway in the aetiology of schizophrenia.

Conclusion

Considered alone, k-means clustering is the most effective of the four methods on typical microarray brain expression datasets. However, memISA and ISA can add extra high-quality clusters to the set produced by k-means, so combining these three methods is the method of choice.

Optimisation of region-specific reference gene selection and relative gene expression analysis methods for pre-clinical trials of Huntington's disease

Background

Transcriptional dysregulation is an early, key pathogenic mechanism in Huntington's disease (HD). Therefore, gene expression analyses have biomarker potential for measuring therapeutic efficacy in pre-clinical trials, particularly those aimed at correcting gene expression abnormalities. Housekeeping genes are commonly used as endogenous references in gene expression studies. However, a systematic study comparing the suitability of candidate reference genes for use in HD mouse models has not been performed. To remedy this situation, 12 housekeeping genes were examined to identify suitable reference genes for use in expression assays.

Results

We found that commonly used reference genes are dysregulated at later time points in the R6/2 mouse model of HD. Therefore, in order to reliably measure gene expression changes for use as pre-clinical trial biomarkers, we set out to identify suitable reference genes for use in R6/2 mice. The expression of potential reference genes was examined in striatum, cortex and cerebellum from 15 week old R6/2 and matched wild-type littermates. Expression levels of candidate reference genes varied according to genotype and brain region. GeNorm software was used to identify the three most stably expressed genes for each brain region. Relative quantification methods using the geometric mean of three reference genes for normalisation enables accurate determination of gene expression levels in wild-type and R6/2 mouse brain regions.

Conclusion

Our study has identified a reproducible, reliable method by which we able to accurately determine the relative expression level of target genes in specific brain regions, thus increasing the potential of gene expression analysis as a biomarker in HD pre-clinical trials.

Transcriptome analysis identifies genes with enriched expression in the mouse central extended amygdala

The central extended amygdala (EAc) is an ensemble of highly interconnected limbic structures of the anterior brain, and forms a cellular continuum including the bed nucleus of the stria terminalis (BNST), the central nucleus of the amygdala (CeA) and the nucleus accumbens shell (AcbSh). This neural network is a key site for interactions between brain reward and stress systems, and has been implicated in several aspects of drug abuse. In order to increase our understanding of EAc function at the molecular level, we undertook a genome-wide screen (Affymetrix) to identify genes whose expression is enriched in the mouse EAc. We focused on the less-well known BNST-CeA areas of the EAc, and identified 121 genes that exhibit more than twofold higher expression level in the EAc compared with whole brain. Among these, 43 genes have never been described to be expressed in the EAc. We mapped these genes throughout the brain, using non-radioactive in situ hybridization, and identified eight genes with a unique and distinct rostro-caudal expression pattern along AcbSh, BNST and CeA. Q-PCR analysis performed in brain and peripheral organ tissues indicated that, with the exception of one (Spata13), all these genes are predominantly expressed in brain. These genes encode signaling proteins (Adora2, GPR88, Arpp21 and Rem2), a transcription factor (Limh6) or proteins of unknown function (Rik130, Spata13 and Wfs1). The identification of genes with enriched expression expands our knowledge of EAc at a molecular level, and provides useful information to toward genetic manipulations within the EAc.

Thyroid hormone action in the adult brain: gene expression profiling of the effects of single and multiple doses of triiodo-L-thyronine in the rat striatum

Thyroid hormones have profound effects on mood and behavior, but the molecular basis of thyroid hormone action in the adult brain is relatively unknown. In particular, few thyroid hormone-dependent genes have been identified in the adult brain despite extensive work carried out on the developing brain. In this work we performed global analysis of gene expression in the adult rat striatum in search for genomic changes taking place after administration of T(3) to hypothyroid rats. The hormone was administered in two different schedules: 1) a single, large dose of 25 microg per 100 g body weight (SD) or 2) 1.5 microg per 100 g body weight once daily for 5 d (RD). Twenty-four hours after the single or last of multiple doses, gene expression in the striatum was analyzed using Codelink microarrays. SD caused up-regulation of 149 genes and down-regulation of 88 genes. RD caused up-regulation of 18 genes and down-regulation of one gene. The results were confirmed by hybridization to Affymetrix microarrays and by TaqMan PCR. Among the genes identified are genes involved in circadian regulation and the regulation of signaling pathways in the striatum. These results suggest that thyroid hormone is involved in regulation of striatal physiology at multiple control points. In addition, they may explain the beneficial effects of large doses of thyroid hormone in bipolar disorders.

Identification of a set of genes showing regionally enriched expression in the mouse brain

Background

The Pleiades Promoter Project aims to improve gene therapy by designing human mini-promoters (< 4 kb) that drive gene expression in specific brain regions or cell-types of therapeutic interest. Our goal was to first identify genes displaying regionally enriched expression in the mouse brain so that promoters designed from orthologous human genes can then be tested to drive reporter expression in a similar pattern in the mouse brain.

Results

We have utilized LongSAGE to identify regionally enriched transcripts in the adult mouse brain. As supplemental strategies, we also performed a meta-analysis of published literature and inspected the Allen Brain Atlas in situ hybridization data. From a set of approximately 30,000 mouse genes, 237 were identified as showing specific or enriched expression in 30 target regions of the mouse brain. GO term over-representation among these genes revealed co-involvement in various aspects of central nervous system development and physiology.

Conclusion

Using a multi-faceted expression validation approach, we have identified mouse genes whose human orthologs are good candidates for design of mini-promoters. These mouse genes represent molecular markers in several discrete brain regions/cell-types, which could potentially provide a mechanistic explanation of unique functions performed by each region. This set of markers may also serve as a resource for further studies of gene regulatory elements influencing brain expression.

Auditory hallucinations in schizophrenia: the role of cognitive, brain structural and genetic disturbances in the left temporal lobe

In this article we review research in our laboratory on auditory hallucinations using behavioral and MRI measure. The review consists of both previously published and new data that for the first time is presented together in a cohesive way. Auditory hallucinations are among the most common symptoms in schizophrenia, affecting more than 70% of the patients. We here advance the hypothesis that auditory hallucinations are internally generated speech perceptions that are lateralized to the left temporal lobe, in the peri-Sylvian region. From this we predict that hallucinating patients should have problems identifying a simultaneously presented external speech sound, as measured through performance on the dichotic listening (DL) paradigm with consonant-vowel syllables, since this technique lateralizes the stimulus input. Across a series of behavioral experiments, we have shown that patients with schizophrenia who experience frequent auditory hallucinations fail to demonstrate an expected right ear advantage on the dichotic listening test. Absence of a right ear advantage is indicative of a functional deficit in the left peri-Sylvian region. The results also revealed that patients with ongoing auditory hallucinations were more impaired than patients with previous hallucinations, and that a higher score on the hallucination item in a standard symptom rating scale (BPRS) correlated negatively with number of correct reports for the right ear stimulus. Moreover, we have found that schizophrenia patients fail to shift attention to the left ear stimulus, when explicitly instructed to focus on the right or left ear stimulus only, thus showing a deficit in inhibition of attention and response-inhibition. The behavioral DL data are substantiated in two MR morphometry studies that revealed significant reductions in grey matter density in the left peri-Sylvian region in hallucinating patients, and patients with reduced left temporal lobe grey matter density. Hallucinating patients also failed to show a right ear advantage in the dichotic listening test. Ongoing fMRI studies are focused on the underlying synaptic and molecular mechanisms by investigating the effects of the glutamate antagonist drug memantine on auditory perception and speech lateralization, and examination of temporal cortex-specific gene expression in the left peri-Sylvian region.

Sharing and reusing gene expression profiling data in neuroscience

As public availability of gene expression profiling data increases, it is natural to ask how these data can be used by neuroscientists. Here we review the public availability of high-throughput expression data in neuroscience and how it has been reused, and tools that have been developed to facilitate reuse. There is increasing interest in making expression data reuse a routine part of the neuroscience tool-kit, but there are a number of challenges. Data must become more readily available in public databases; efforts to encourage investigators to make data available are important, as is education on the benefits of public data release. Once released, data must be better-annotated. Techniques and tools for data reuse are also in need of improvement. Integration of expression profiling data with neuroscience-specific resources such as anatomical atlases will further increase the value of expression data.

Synaptic activity-induced global gene expression patterns in the dentate gyrus of adult behaving rats: induction of immunity-linked genes

Gene expression in adult neuronal circuits is dynamically modulated in response to synaptic activity. Persistent changes in synaptic strength, as seen during high-frequency stimulation (HFS)-induced long-term potentiation (LTP), require new gene expression. While modulation of many individual genes has been shown, an understanding of LTP as a complex dynamical response requires elucidation of the global gene expression signature and its impact on biologically meaningful gene sets. In this study, we demonstrate that LTP induction in the dentate gyrus of awake freely moving rats was associated with changes in the expression of genes linked to signal transduction, protein trafficking, cell structure and motility, and other processes consistent with the induction of mechanisms of synaptic reorganization and growth. Interestingly, the most significantly over-represented gene sets were related to immunity and defense, including T-cell-mediated immunity and major histocompatibility complex (MHC) class I-mediated immunity. Real-time PCR confirmed the upregulation of a panel of immune-linked genes including the rt1-a/ce family, and the MHC class II members cd74, rt1-Ba and rt1-Da. These genes were N-methyl-d-aspartate receptor-independent and not induced following HFS-LTP induction in anesthetized rats, indicating a gene response specific to behaving rats. Our data support recent assumptions that immunity-associated processes are functionally linked to adaptive neuronal responses in the brain, although the differential expression of immunity-linked genes could also be related to the HFS per se.