Histone Modifications in a Mouse Model of Early Adversities and Panic Disorder: Role for Asic1 and Neurodevelopmental Genes

Davide Cittaro, Valentina Lampis, Alessandra Luchetti, Roberto Coccurello, Alessandro Guffanti, Armando Felsani, Anna Moles, Elia Stupka, Francesca R D' Amato, Marco Battaglia, Davide Cittaro, Valentina Lampis, Alessandra Luchetti, Roberto Coccurello, Alessandro Guffanti, Armando Felsani, Anna Moles, Elia Stupka, Francesca R D' Amato, Marco Battaglia

Abstract

Hyperventilation following transient, CO2-induced acidosis is ubiquitous in mammals and heritable. In humans, respiratory and emotional hypersensitivity to CO2 marks separation anxiety and panic disorders, and is enhanced by early-life adversities. Mice exposed to the repeated cross-fostering paradigm (RCF) of interference with maternal environment show heightened separation anxiety and hyperventilation to 6% CO2-enriched air. Gene-environment interactions affect CO2 hypersensitivity in both humans and mice. We therefore hypothesised that epigenetic modifications and increased expression of genes involved in pH-detection could explain these relationships. Medullae oblongata of RCF- and normally-reared female outbred mice were assessed by ChIP-seq for H3Ac, H3K4me3, H3K27me3 histone modifications, and by SAGE for differential gene expression. Integration of multiple experiments by network analysis revealed an active component of 148 genes pointing to the mTOR signalling pathway and nociception. Among these genes, Asic1 showed heightened mRNA expression, coherent with RCF-mice's respiratory hypersensitivity to CO2 and altered nociception. Functional enrichment and mRNA transcript analyses yielded a consistent picture of enhancement for several genes affecting chemoception, neurodevelopment, and emotionality. Particularly, results with Asic1 support recent human findings with panic and CO2 responses, and provide new perspectives on how early adversities and genes interplay to affect key components of panic and related disorders.

Figures

Figure 1
Figure 1
(A) Profile heatmap around TSS of RefSeq genes. Read counts were extracted for all ChIP-seq experiments within a region spanning ±5 kb around TSS. The gradient blue-to-red color indicates high-to-low counts in the corresponding region. (B) Clustered heatmap of ChIP-seq profiles. Correlations were evaluated over the genome-wide signal of all experiments, including genomic DNA, and shown as a diagonal matrix. Rows and columns are clustered according to inter-sample correlation. Color intensity is proportional to the correlation value reported in each cell. Markers for actively transcribed regions (H3K4me3 and H3Ac) cluster together in the upper left region. H3K27me3 clusters with input DNA, probably due to its broad enrichment profile.
Figure 2
Figure 2
(A) Active circuit of 145 genes identified by network analysis (see methods). Each node is colored by the fold change found by SAGE when comparing RCF vs. CT mice. (B) Functional analysis of the active circuit. Scores are combined for EnrichR analysis on Reactome Pathways and MGI Mouse Phenotype gene sets.
Figure 3. Asic1 gene Expression.
Figure 3. Asic1 gene Expression.
Histograms on the left represent mean ± SD levels of expression of Asic1 gene found by SAGE analysis in Set 2 MOs of RCF and CT mice, as described in the main text: mean log 2(E) ± SD Fold Change respectively: 6.88 ± 0.37 and 5.89 ± 0.17, p = 0.047. Histograms on the right represent mean ± SD levels of expression of Asic1 gene found by real time PCR in the MO of 5 RCF and 4 CT mice: mean ± SD of Delta Cycle Threshold (dCT) as referred to Atcb housekeeping gene, respectively: 9.52 ± 1.16 and 12.10 ± 1.38, Student t = 2.64, p = 0.033. Asterisks over the RCF bars denote significant difference in comparison to the CT group.
Figure 4. Respiratory responses.
Figure 4. Respiratory responses.
Mean ± SD values of the percentage of increment (delta) of tidal volume from baseline (air) to 6% CO2-enriched air, among 5 RCF female (26.88 ± 8.10) vs. 5 CT female (16.63 ± 4.12) adult (PND = 80–90) mice (Student t = 2.54, DF = 8, p = 0.03), and 5 RCF male (30.46 + 6.34) vs. 3 CT male (20.68 + 2.02) adult (PND = 80–90) mice (Student t = 2.53, DF = 6, p = 0.05). Asterisk indicates a significance of p ≤ 0.05 in the comparisons of RCF to CT animals divided by sex. The results of an ANOVA with females and males pooled together are reported in the main text.
Figure 5. Nociception.
Figure 5. Nociception.
Results of the formalin test performed in 6 RCF vs. 7 CT adult (PND = 80–90) male mice. Responses were recorded continuously for 40 minutes and calculated in blocks of consecutive 5 minute periods. The figure shows the time spent licking the injected paw (Mean ± SE in seconds) for each time interval; ‘*’ and ‘**’ indicate significance of respectively p 

References

    1. Guyenet P. G. et al.. Central CO2 chemoreception and integrated neural mechanisms of cardiovascular and respiratory control. J Appl Physiol 108, 995–1002 (2010).
    1. Brannan S. et al.. Neuroimaging of cerebral activations and deactivations associated with hypercapnia and hunger for air. Proc Natl Acad Sci USA 98, 2029–2034 (2001).
    1. Nattie E. CO2 Brainstem chemoreceptors and breathing. Prog Neurobiol 59, 299–331 (1999).
    1. Ziemann A. E. et al.. The amygdala is a chemosensor that detects carbon dioxide and acidosis to elicit fear behavior. Cell 139, 1012–1 (2009).
    1. Han F. & Strohl K. P. Inheritance of ventilatory behavior in rodent models. Respir Physiol 121, 247–256 (2000).
    1. Battaglia M. et al.. A genetic study of the acute anxious response to carbon dioxide stimulation in man. J Psychiatr Res 41, 906–917 (2007).
    1. Battaglia M. et al.. A genetically informed study of the association between childhood separation anxiety, sensitivity to CO2, panic disorder, and the effect of childhood parental loss. Arch Gen Psychiatry 66, 64–71 (2009).
    1. Klein D. F. False suffocation alarms, spontaneous panics, and related conditions. An integrative hypothesis. Arch Gen Psychiatry 50, 306–317 (1993).
    1. Battaglia M., Bertella S., Ogliari A., Bellodi L. & Smeraldi E. Modulation by muscarinic antagonists of the response to carbon dioxide challenge in panic disorder. Arch Gen Psychiatry 58, 114–119 (2001).
    1. Esquivel G., Schruers K. R., Maddock R. J., Colasanti A. & Griez E. J. Acids in the brain, a factor in panic? J Psychopharmacol 24, 639–647 (2010).
    1. Magnotta V. A. et al.. Detecting activity-evoked pH changes in human brain. Proc Natl Acad Sci USA 109, 8270–8273 (2012).
    1. Vollmer L. L., Strawn J. R. & Sah R. Acid–base dysregulation and chemosensory mechanisms in panic disorder, a translational update. Transl Psychiatry 5, e572 (2015).
    1. Spatola C. A. et al.. Gene- environment interactions in panic disorder and CO sensitivity, Effects of events occurring early in life. Am J Med Genet B Neuropsychiatr Genet 156, 79–88 (2011).
    1. Battaglia M., Ogliari A., D’Amato F. & Kinkead R. Early-life risk factors for panic and separation anxiety disorder, Insights and outstanding questions arising from human and animal studies of CO2 sensitivity. Neurosci Biobehav Rev 46, 455–464 (2014).
    1. D’Amato F. R. et al.. Unstable maternal environment, separation anxiety, and heightened CO2 sensitivity induced by gene-by-environment interplay. PLoS One 6, 18637 (2011).
    1. Luchetti A. et al.. Early manipulations effects in mice, differences and similarities between postnatal handling and repeated cross-fostering. Front Behav Neurosci 9, 93 (2015).
    1. Ventura R. et al.. Postnatal aversive experience impairs sensitivity to natural rewards and increases susceptibility to negative events in adult life. Cereb Cortex 23, 1606–1617 (2013).
    1. Feil R. & Fraga M. F. Epigenetics and the environment, emerging patterns and implications. Nat Rev Genet 13, 97–109 (2012).
    1. Xie L., Korkmaz K. S., Braun K. & Bock J. Early life stress-induced histone acetylations correlate with activation of the synaptic plasticity genes Arc and Egr1 in the mouse hippocampus. J Neurochem 125, 457–464 (2013).
    1. Dai H. Histone Modification Patterns and Their Responses to Environment. Curr Environ Health Rep 1, 11–21 (2014).
    1. Zhou V. W., Goren A. & Bernstein B. E. Charting histone modifications and the functional organization of mammalian genomes. Nat Rev Genet 12, 7–18 (2011).
    1. Goossens L. et al.. Brainstem response to hypercapnia, A symptom provocation study into the pathophysiology of panic disorder. J Psychopharmacol, 28, 449–456 (2014).
    1. Barski A. et al.. High-resolution profiling of histone methylations in the human genome. Cell 129, 823–837 (2007).
    1. Lyu D. et al.. The mTOR signaling pathway regulates pain-related synaptic plasticity in rat entorhinal-hippocampal pathways. Mol pain 9, 64 (2013).
    1. Wemmie J. A., Taugher R. J. & Kreple C. J. Acid-sensing ion channels in pain and disease. Nat Rev Neurosci 14, 461–471 (2013).
    1. Smoller J. W. et al.. The Human Ortholog of Acid-Sensing Ion Channel Gene ASIC1a Is Associated With Panic Disorder and Amygdala Structure and Function. Biol Psychiatry 76, 902–910 (2014).
    1. Savage J. E. et al.. Validation of Candidate Anxiety Disorder Genes Using a Carbon Dioxide Challenge Task. Biol Psychol 109, 61–66 (2015).
    1. Dickinson A. L., Leach M. C. & Flecknell P. A. Influence of early neonatal experience on nociceptive responses and analgesic effects in rats. Lab anim 43, 11–16 (2009).
    1. McWilliams L. A., Goodwin R. D. & Cox B. J. Depression and anxiety associated with three pain conditions, results from a nationally representative sample. Pain 111, 77–83 (2004).
    1. Gadermann A. M., Alonso J., Vilagut G., Zaslavsky A. M. & Kessler R. C. Comorbidity and disease burden in the national comorbidity survey replication (NCS-R). Depress anxiety 29, 797–806 (2012).
    1. Yuzaki M. Cbln1 and its family proteins in synapse formation and maintenance. Curr Opin Neurobiol 21, 215–220 (2011).
    1. Caronia G., Wilcoxon J., Feldman P. & Grove E. A. Bone morphogenetic protein signaling in the developing telencephalon controls formation of the hippocampal dentate gyrus and modifies fear-related behavior. J Neurosci 30, 6291–6301 (2010).
    1. Samanta J. et al.. BMPR1a signaling determines numbers of oligodendrocytes and calbindin-expressing interneurons in the cortex. J Neurosci 27, 7397–7407 (2007).
    1. Canterini S., Mangia F. & Fiorenza M. T. Thg‐1pit gene expression in granule cells of the developing mouse brain and in their synaptic targets, mature Purkinje, and mitral cells. Dev dyn 234, 689–697 (2005).
    1. Walton N. M. et al.. Gastrin‐Releasing Peptide Contributes to the Regulation of Adult Hippocampal Neurogenesis and Neuronal Development. Stem Cells 32, 2454–2466 (2014).
    1. Roesler R., Kent P., Luft T., Schwartsmann G. & Merali Z. Gastrin-releasing peptide receptor signaling in the integration of stress and memory. Neurobiol Learn Memory 112, 44–52 (2014).
    1. Garcia V. A. et al.. Neonatal gastrin-releasing peptide receptor blockade reduces maternal odor preference in rats. Behav brain res 214, 456–459 (2010).
    1. Traka M. et al.. WDR81 Is Necessary for Purkinje and Photoreceptor Cell Survival. J Neurosci 33, 6834–6844 (2013).
    1. Liang H. et al.. Neural development is dependent on the function of specificity protein 2 in cell cycle progression. Development 140, 552–561 (2013).
    1. Taugher R. J. et al.. The bed nucleus of the stria terminalis is critical for anxiety-related behavior evoked by CO2 and acidosis. J Neurosci 34, 10247–10255 (2014).
    1. Richey J. A., Schmidt N. B., Hofmann S. G. & Timpano K. R. Temporal and structural dynamics of anxiety sensitivity in predicting fearful responding to a 35% CO2 challenge. J Anxiety Disord 24, 423–432 (2010).
    1. Roberson-Nay R. et al.. Evidence for Distinct Genetic Effects Associated with Response to 35% CO2. Depress Anxiety 30, 259–266 (2013).
    1. Davies M. N. et al.. Functional annotation of the human brain methylome identifies tissue-specific epigenetic variation across brain and blood. Genome Biol 13, R43 (2012).
    1. Huang H. S., Matevossian A., Jiang Y. & Akbarian S. Chromatin immunoprecipitation in postmortem brain. J Neurosci Methods 156, 284–92 (2006).
    1. Zhang Y. et al.. Model- based analysis of ChIP-Seq (MACS). Genome Biol 9, R137 (2008).
    1. Li Q., Brown J. B., Huang H. & Bickel P. J. Measuring reproducibility of high-throughput experiments. Ann Appl Stat 5, 1752–1779 (2011).
    1. Robinson M. D., McCarthy D. J. & Smyth G. K. EdgeR, a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140 (2009).
    1. Pruitt K. D. et al.. RefSeq, an update on mammalian reference sequences. Nucleic Acids Res 42, 756–763 (2014).
    1. Ritchie M. E. et al.. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acid Res 43, e47 (2015).
    1. Zuberi K. et al.. GeneMANIA prediction server 2013 update. Nucleic Acids Res 41, W115–W122 (2013).
    1. Mitra K., Carvunis A. R., Ramesh S. K. & Ideker T. Integrative approaches forfinding modular structure inbiological networks. Nat Rev Genet 14, 719–732 (2013)
    1. Chen E. Y. et al.. Enrichr, interactive and collaborative HTML5 gene list enrichment analysis tool. BMC bioinformatics 14, 128 (2013).

Source: PubMed

3
Se inscrever