Prodromal Parkinson disease signs are predicted by a whole-blood inflammatory transcriptional signature in young Pink1-/- rats.

BACKGROUND: Parkinson disease (PD) is the fastest growing neurodegenerative disease. The molecular pathology of PD in the prodromal phase is poorly understood; as such, there are no specific prognostic or diagnostic tests. A validated Pink1 genetic knockout rat was used to model early-onset and progressive PD. Male Pink1-/- rats exhibit progressive declines in ultrasonic vocalizations as well as hindlimb and forelimb motor deficits by mid-to-late adulthood. Previous RNA-sequencing work identified upregulation of genes involved in disease pathways and inflammation within the brainstem and vocal fold muscle. The purpose of this study was to identify gene pathways within the whole blood of young Pink1-/- rats (3 months of age) and to link gene expression to early acoustical changes. To accomplish this, limb motor testing (open field and cylinder tests) and ultrasonic vocalization data were collected, immediately followed by the collection of whole blood and RNA extraction. Illumina® Total RNA-Seq TruSeq platform was used to profile differential expression of genes. Statistically significant genes were identified and Weighted Gene Co-expression Network Analysis was used to construct co-expression networks and modules from the whole blood gene expression dataset as well as the open field, cylinder, and USV acoustical dataset. ENRICHR was used to identify the top up-regulated biological pathways. RESULTS: The data suggest that inflammation and interferon signaling upregulation in the whole blood is present during early PD. We also identified genes involved in the dysregulation of ribosomal protein and RNA processing gene expression as well as prion protein gene expression. CONCLUSIONS: These data identified several potential blood biomarkers and pathways that may be linked to anxiety and vocalization acoustic parameters and are key candidates for future drug-repurposing work and comparison to human datasets.

Large-scale gene expression changes in APP/PSEN1 and GFAP mutation models exhibit high congruence with Alzheimer's disease.

Alzheimer's disease (AD) is a complex neurodegenerative disorder with both genetic and non-genetic causes. Animal research models are available for a multitude of diseases and conditions affecting the central nervous system (CNS), and large-scale CNS gene expression data exist for many of these. Although there are several models specifically for AD, each recapitulates different aspects of the human disease. In this study we evaluate over 500 animal models to identify those with CNS gene expression patterns matching human AD datasets. Approaches included a hypergeometric based scoring system that rewards congruent gene expression patterns but penalizes discordant gene expression patterns. The top two models identified were APP/PS1 transgenic mice expressing mutant APP and PSEN1, and mice carrying a GFAP mutation that is causative of Alexander disease, a primary disorder of astrocytes in the CNS. The APP/PS1 and GFAP models both matched over 500 genes moving in the same direction as in human AD, and both had elevated GFAP expression and were highly congruent with one another. Also scoring highly were the 5XFAD model (with five mutations in APP and PSEN1) and mice carrying CK-p25, APP, and MAPT mutations. Animals with the APOE3 and 4 mutations combined with traumatic brain injury ranked highly. Bulbectomized rats scored high, suggesting anosmia could be causative of AD-like gene expression. Other matching models included the SOD1G93A strain and knockouts for SNORD116 (Prader-Willi mutation), GRID2, INSM1, XBP1, and CSTB. Many top models demonstrated increased expression of GFAP, and results were similar across multiple human AD datasets. Heatmap and Uniform Manifold Approximation Plot results were consistent with hypergeometric ranking. Finally, some gene manipulation models, including for TYROBP and ATG7, were identified with reversed AD patterns, suggesting possible neuroprotective effects. This study provides insight for the pathobiology of AD and the potential utility of available animal models.

Prodromal Parkinson disease signs are predicted by a whole-blood inflammatory transcriptional signature in young Pink1-/- rats.

Background: Parkinson disease (PD) is the fastest growing neurodegenerative disease. The molecular pathology of PD in the prodromal phase is poorly understood; as such, there are no specific prognostic or diagnostic tests. A validated Pink1 genetic knockout rat was used to model early-onset and progressive PD. Male Pink1-/- rats exhibit progressive declines in ultrasonic vocalizations as well as hindlimb and forelimb motor deficits by mid-to-late adulthood. Previous RNA-sequencing work identified upregulation of genes involved in disease pathways and inflammation within the brainstem and vocal fold muscle. The purpose of this study was to identify gene pathways within the whole blood of young Pink1-/- rats (3 months of age) and to link gene expression to early acoustical changes. To accomplish this, limb motor testing (open field and cylinder tests) and ultrasonic vocalization data were collected, immediately followed by the collection of whole blood and RNA extraction. Illumina® Total RNA-Seq TruSeq platform was used to profile differential expression of genes. Statistically significant genes were identified and Weighted Gene Co-expression Network Analysis was used to construct co-expression networks and modules from the whole blood gene expression dataset as well as the open field, cylinder, and USV acoustical dataset. ENRICHR was used to identify the top up-regulated biological pathways. Results: The data suggest that inflammation and interferon signaling upregulation in the whole blood is present during early PD. We also identified genes involved in the dysregulation of ribosomal protein and RNA processing gene expression as well as prion protein gene expression. Conclusions: These data identified several potential blood biomarkers and pathways that may be linked to anxiety and vocalization acoustic parameters and are key candidates for future drug-repurposing work and comparison to human datasets.

Distinct patterns of gene expression in the medial preoptic area are related to gregarious singing behavior in European starlings (Sturnus vulgaris).

BACKGROUND: Song performed in flocks by European starlings (Sturnus vulgaris), referred to here as gregarious song, is a non-sexual, social behavior performed by adult birds. Gregarious song is thought to be an intrinsically reinforced behavior facilitated by a low-stress, positive affective state that increases social cohesion within a flock. The medial preoptic area (mPOA) is a region known to have a role in the production of gregarious song. However, the neurochemical systems that potentially act within this region to regulate song remain largely unexplored. In this study, we used RNA sequencing to characterize patterns of gene expression in the mPOA of male and female starlings singing gregarious song to identify possibly novel neurotransmitter, neuromodulator, and hormonal pathways that may be involved in the production of gregarious song. RESULTS: Differential gene expression analysis and rank rank hypergeometric analysis indicated that dopaminergic, cholinergic, and GABAergic systems were associated with the production of gregarious song, with multiple receptor genes (e.g., DRD2, DRD5, CHRM4, GABRD) upregulated in the mPOA of starlings who sang at high rates. Additionally, co-expression network analyses identified co-expressing gene clusters of glutamate signaling-related genes associated with song. One of these clusters contained five glutamate receptor genes and two glutamate scaffolding genes and was significantly enriched for genetic pathways involved in neurodevelopmental disorders associated with social deficits in humans. Two of these genes, GRIN1 and SHANK2, were positively correlated with performance of gregarious song. CONCLUSIONS: This work provides new insights into the role of the mPOA in non-sexual, gregarious song in starlings and highlights candidate genes that may play a role in gregarious social interactions across vertebrates. The provided data will also allow other researchers to compare across species to identify conserved systems that regulate social behavior.

Thyroarytenoid Oxidative Metabolism and Synaptic Signaling Dysregulation in the Female Pink1-/- Rat.

OBJECTIVES AND HYPOTHESIS: Vocal dysfunction, including hypophonia, in Parkinson disease (PD) manifests in the prodromal period and significantly impacts an individual's quality of life. Data from human studies suggest that pathology leading to vocal deficits may be structurally related to the larynx and its function. The Pink1-/- rat is a translational model used to study pathogenesis in the context of early-stage mitochondrial dysfunction. The primary objective of this work was to identify differentially expressed genes in the thyroarytenoid muscle and examine the dysregulated biological pathways in the female rat. METHODS: RNA sequencing was used to determine thyroarytenoid (TA) muscle gene expression in adult female Pink1-/- rats compared with controls. A bioinformatic approach and the ENRICHR gene analysis tool were used to compare the sequencing dataset with biological pathways and processes, disease relationships, and drug-repurposing compounds. Weighted Gene Co-expression Network Analysis was used to construct biological network modules. The data were compared with a previously published dataset in male rats. RESULTS: Significant upregulated pathways in female Pink1-/- rats included fatty acid oxidation and muscle contraction, synaptic transmission, and neuromuscular processes. Downregulated pathways included anterograde transsynaptic signaling, chemical synaptic transmission, and ion release. Several drug treatment options including cetuximab, fluoxetine, and resveratrol are hypothesized to reverse observed genetic dysregulation. CONCLUSIONS: Data presented here are useful for identifying biological pathways that may underlie the mechanisms of peripheral dysfunction including neuromuscular synaptic transmission to the TA muscle. These experimental biomarkers have the potential to be targeted as sites for improving the treatment for hypophonia in early-stage PD. LEVEL OF EVIDENCE: NA Laryngoscope, 133:3412-3421, 2023.

Alzheimer's disease large-scale gene expression portrait identifies exercise as the top theoretical treatment.

Alzheimer's disease (AD) is a complex neurodegenerative disorder that affects multiple brain regions and is difficult to treat. In this study we used 22 AD large-scale gene expression datasets to identify a consistent underlying portrait of AD gene expression across multiple brain regions. Then we used the portrait as a platform for identifying treatments that could reverse AD dysregulated expression patterns. Enrichment of dysregulated AD genes included multiple processes, ranging from cell adhesion to CNS development. The three most dysregulated genes in the AD portrait were the inositol trisphosphate kinase, ITPKB (upregulated), the astrocyte specific intermediate filament protein, GFAP (upregulated), and the rho GTPase, RHOQ (upregulated). 41 of the top AD dysregulated genes were also identified in a recent human AD GWAS study, including PNOC, C4B, and BCL11A. 42 transcription factors were identified that were both dysregulated in AD and that in turn affect expression of other AD dysregulated genes. Male and female AD portraits were highly congruent. Out of over 250 treatments, three datasets for exercise or activity were identified as the top three theoretical treatments for AD via reversal of large-scale gene expression patterns. Exercise reversed expression patterns of hundreds of AD genes across multiple categories, including cytoskeleton, blood vessel development, mitochondrion, and interferon-stimulated related genes. Exercise also ranked as the best treatment across a majority of individual region-specific AD datasets and meta-analysis AD datasets. Fluoxetine also scored well and a theoretical combination of fluoxetine and exercise reversed 549 AD genes. Other positive treatments included curcumin. Comparisons of the AD portrait to a recent depression portrait revealed a high congruence of downregulated genes in both. Together, the AD portrait provides a new platform for understanding AD and identifying potential treatments for AD.

Evaluation of animal model congruence to human depression based on large-scale gene expression patterns of the CNS.

Depression is a complex mental health disorder that is difficult to study. A wide range of animal models exist and for many of these data on large-scale gene expression patterns in the CNS are available. The goal of this study was to evaluate how well animal models match human depression by evaluating congruence and discordance of large-scale gene expression patterns in the CNS between almost 300 animal models and a portrait of human depression created from male and female datasets. Multiple approaches were used, including a hypergeometric based scoring system that rewards common gene expression patterns (e.g., up-up or down-down in both model and human depression), but penalizes opposing gene expression patterns. RRHO heat maps, Uniform Manifold Approximation Plot (UMAP), and machine learning were used to evaluate matching of models to depression. The top ranked model was a histone deacetylase (HDAC2) conditional knockout in forebrain neurons. Also highly ranked were various models for Alzheimer's, including APPsa knock-in (2nd overall), APP knockout, and an APP/PS1 humanized double mutant. Other top models were the mitochondrial gene HTRA2 knockout (that is lethal in adulthood), a modified acetylcholinesterase, a Huntington's disease model, and the CRTC1 knockout. Over 30 stress related models were evaluated and while some matched highly with depression, others did not. In most of the top models, a consistent dysregulation of MAP kinase pathway was identified and the genes NR4A1, BDNF, ARC, EGR2, and PDE7B were consistently downregulated as in humans with depression. Separate male and female portraits of depression were also evaluated to identify potential sex specific depression matches with models. Individual human depression datasets were also evaluated to allow for comparisons across the same brain regions. Heatmap, UMAP, and machine learning results supported the hypergeometric ranking findings. Together, this study provides new insights into how large-scale gene expression patterns may be similarly dysregulated in some animals models and humans with depression that may provide new avenues for understanding and treating depression.

Thyroarytenoid Muscle Gene Expression in a Rat Model of Early-Onset Parkinson's Disease.

OBJECTIVES/HYPOTHESIS: Voice disorders in Parkinson's disease (PD) are early-onset, manifest in the preclinical stages of the disease, and negatively impact quality of life. The complete loss of function in the PTEN-induced kinase 1 gene (Pink1) causes a genetic form of early-onset, autosomal recessive PD. Modeled after the human inherited mutation, the Pink1-/- rat demonstrates significant cranial sensorimotor dysfunction including declines in ultrasonic vocalizations. However, the underlying genetics of the vocal fold thyroarytenoid (TA) muscle that may contribute to vocal deficits has not been studied. The aim of this study was to identify differentially expressed genes in the TA muscle of 8-month-old male Pink1-/- rats compared to wildtype controls. STUDY DESIGN: Animal experiment with control. METHODS: High throughput RNA sequencing was used to examine TA muscle gene expression in adult male Pink1-/- rats and wildtype controls. Weighted Gene Co-expression Network Analysis was used to construct co-expression modules to identify biological networks, including where Pink1 was a central node. The ENRICHR tool was used to compare this gene set to existing human gene databases. RESULTS: We identified 134 annotated differentially expressed genes (P < .05 cutoff) and observed enrichment in the following biological pathways: Parkinson's disease (Casp7, Pink1); Parkin-Ubiquitin proteasome degradation (Psmd12, Psmd7); MAPK signaling (Casp7, Ppm1b, Ppp3r1); and inflammatory TNF-α, Nf-κB Signaling (Casp7, Psmd12, Psmd7, Cdc34, Bcl7a, Peg3). CONCLUSIONS: Genes and pathways identified here may be useful for evaluating the specific mechanisms of peripheral dysfunction including within the laryngeal muscle and have potential to be used as experimental biomarkers for treatment development. LEVEL OF EVIDENCE: NA Laryngoscope, 131:E2874-E2879, 2021.

Creation of a gene expression portrait of depression and its application for identifying potential treatments.

Depression is a complex mental health disorder and the goal here was to identify a consistent underlying portrait of expression that ranks all genes from most to least dysregulated and indicates direction of change relative to controls. Using large-scale neural gene expression depression datasets, a combined portrait (for men and women) was created along with one for men and one for women only. The depressed brain was characterized by a "hypo" state, that included downregulation of activity-related genes, including EGR1, FOS, and ARC, and indications of a lower brain temperature and sleep-like state. MAP kinase and BDNF pathways were enriched with overlapping genes. Expression patterns suggested decreased signaling for GABA and for neuropeptides, CRH, SST, and CCK. GWAS depression genes were among depression portrait genes and common genes of interest included SPRY2 and PSEN2. The portraits were used with the drug repurposing approach of signature matching to identify treatments that could reverse depression gene expression patterns. Exercise was identified as the top treatment for depression for the combined and male portraits. Other non-traditional treatments that scored well were: curcumin, creatine, and albiflorin. Fluoxetine scored best among typical antidepressants. The creation of the portraits of depression provides new insights into the complex landscape of depression and a novel platform for evaluating and identifying potential new treatments.

Mu opioid receptors in the medial preoptic area govern social play behavior in adolescent male rats.

Neural systems underlying important behaviors are usually highly conserved across species. The medial preoptic area (MPOA) has been demonstrated to play a crucial role in reward associated with affiliative, nonsexual, social communication in songbirds. However, the role of MPOA in affiliative, rewarding social behaviors (eg, social play behavior) in mammals remains largely unknown. Here we applied our insights from songbirds to rats to determine whether opioids in the MPOA govern social play behavior in rats. Using an immediate early gene (ie, Egr1, early growth response 1) expression approach, we identified increased numbers of Egr1-labeled cells in the MPOA after social play in adolescent male rats. We also demonstrated that cells expressing mu opioid receptors (MORs, gene name Oprm1) in the MPOA displayed increased Egr1 expression when adolescent rats were engaged in social play using double immunofluorescence labeling of MOR and Egr1. Furthermore, using short hairpin RNA-mediated gene silencing we revealed that knockdown of Oprm1 in the MPOA reduced the number of total play bouts and the frequency of pouncing. Last, RNA sequencing differential gene expression analysis identified genes involved in neuronal signaling with altered expression after Oprm1 knockdown, and identified Egr1 as potentially a key modulator for Oprm1 in the regulation of social play behavior. Altogether, these results show that the MPOA is involved in social play behavior in adolescent male rats and support the hypothesis that the MPOA is part of a conserved neural circuit across vertebrates in which opioids act to govern affiliative, intrinsically rewarded social behaviors.

Complex patterns of dopamine-related gene expression in the ventral tegmental area of male zebra finches relate to dyadic interactions with long-term female partners.

Dopaminergic projections from the ventral tegmental area (VTA) to multiple efferent targets are implicated in pair bonding, yet the role of the VTA in the maintenance of long-term pair bonds is not well characterized. Complex interactions between numerous neuromodulators modify activity in the VTA, suggesting that individual differences in patterns of gene expression in this region may explain individual differences in long-term social interactions in bonded pairs. To test this hypothesis we used RNA-seq to measure expression of over 8000 annotated genes in male zebra finches in established male-female pairs. Weighted gene co-expression network analysis identified a gene module that contained numerous dopamine-related genes with TH found to be the most connected gene of the module. Genes in this module related to male agonistic behaviors as well as bonding-related behaviors produced by female partners. Unsupervised learning approaches identified two groups of males that differed with respect to expression of numerous genes. Enrichment analyses showed that many dopamine-related genes and modulators differed between these groups, including dopamine receptors, synthetic and degradative enzymes, the avian dopamine transporter and several GABA- and glutamate-related genes. Many of the bonding-related behaviors closely associated with VTA gene expression in the two male groups were produced by the male's partner, rather than the male himself. Collectively, results highlight numerous candidate genes in the VTA that can be explored in future studies and raise the possibility that the molecular/genetic organization of the VTA may be strongly shaped by a social partner and/or the strength of the pair bond.

The circadian gene Nr1d1 in the mouse nucleus accumbens modulates sociability and anxiety-related behavior.

Nuclear receptor subfamily 1, group D, member 1 (Nr1d1) (also known as Rev-erb alpha) has been linked to circadian rhythm regulation, mood-related behavior, and disorders associated with social deficits. Recent work from our laboratory found striking decreases in Nr1d1 in the nucleus accumbens (NAc) in the maternal condition and indirect evidence that Nr1d1 was interacting with numerous addiction and reward-related genes to modulate social reward. In this study, we applied our insights from the maternal state to non-parental adult mice to determine whether decreases in Nr1d1 expression in the NAc via adeno-associated viral (AAV) vectors and short hairpin RNA (shRNA)-mediated gene knockdown were sufficient to modulate social behaviors and mood-related behaviors. Knockdown of Nr1d1 in the NAc enhanced sociability, reduced anxiety, but did not affect depressive-like traits in female mice. In male mice, Nr1d1 knockdown had no significant behavioral effects. Microarray analysis of Nr1d1 knockdown in females identified changes in circadian rhythm and histone deacetylase genes and suggested possible drugs, including histone deacetylase inhibitors, that could mimic actions of Nr1d1 knockdown. Quantitative real-time PCR (qPCR) analysis confirmed expression upregulation of genes period circadian clock 1 (Per1) and period circadian clock 2 (Per2) with Nr1d1 knockdown. Evidence for roles for opioid-related genes opioid receptor, delta 1 (Oprd1) and preproenkephalin (Penk) was also found. Together, these results suggest that Nr1d1 in the NAc modulates sociability and anxiety-related behavior in a sex-specific manner and circadian, histone deacetylase, and opioid-related genes may be involved in the expression of these behavioral phenotypes. This article is protected by copyright. All rights reserved.

Down-regulation of fatty acid binding protein 7 (Fabp7) is a hallmark of the postpartum brain.

Fatty acid binding protein 7 (Fabp7) is a versatile protein that is linked to glial differentiation and proliferation, neurogenesis, and multiple mental health disorders. Recent microarray studies identified a robust decrease in Fabp7 expression in key brain regions of the postpartum rodents. Given its diverse functions, Fabp7 could play a critical role in sculpting the maternal brain and promoting the maternal phenotype. The present study aimed at investigating the expression profile of Fabp7 across the postpartum CNS. Quantitative real-time PCR (qPCR) analysis showed that Fabp7 mRNA was consistently down-regulated across the postpartum brain. Of the 9 maternal care-related regions tested, seven exhibited significant decreases in Fabp7 in postpartum (relative to virgin) females, including medial prefrontal cortex (mPFC), nucleus accumbens (NA), lateral septum (LS), bed nucleus of stria terminalis dorsal (BnSTd), paraventricular nucleus (PVN), lateral hypothalamus (LH), and basolateral and central amygdala (BLA/CeA). For both ventral tegmental area (VTA) and medial preoptic area (MPOA) levels of Fabp7 were lower in mothers, but levels of changes did not reach significance. Confocal microscopy revealed that protein expression of Fabp7 in the LS paralleled mRNA findings. Specifically, the caudal LS exhibited a significant reduction in Fabp7 immunoreactivity, while decreases in medial LS were just above significance. Double fluorescent immunolabeling confirmed the astrocytic phenotype of Fabp7-expressing cells. Collectively, this research demonstrates a broad and marked reduction in Fabp7 expression in the postpartum brain, suggesting that down-regulation of Fabp7 may serve as a hallmark of the postpartum brain and contribute to the maternal phenotype.

The circadian gene Nr1d1 in the mouse nucleus accumbens modulates sociability and anxiety-related behaviour.

Nuclear receptor subfamily 1, group D, member 1 (Nr1d1) (also known as Rev-erb alpha) has been linked to circadian rhythm regulation, mood-related behaviour and disorders associated with social deficits. Recent work from our laboratory found striking decreases in Nr1d1 in the nucleus accumbens (NAc) in the maternal condition and indirect evidence that Nr1d1 was interacting with numerous addiction and reward-related genes to modulate social reward. In this study, we applied our insights from the maternal state to nonparental adult mice to determine whether decreases in Nr1d1 expression in the NAc via adeno-associated viral (AAV) vectors and short hairpin RNA (shRNA)-mediated gene knockdown were sufficient to modulate social behaviours and mood-related behaviours. Knockdown of Nr1d1 in the NAc enhanced sociability and reduced anxiety, but did not affect depressive-like traits in female mice. In male mice, Nr1d1 knockdown had no significant behavioural effects. Microarray analysis of Nr1d1 knockdown in females identified changes in circadian rhythm and histone deacetylase genes and suggested possible drugs, including histone deacetylase inhibitors, that could mimic actions of Nr1d1 knockdown. Quantitative real-time PCR (qPCR) analysis confirmed expression upregulation of gene period circadian clock 1 (Per1) and period circadian clock 2 (Per2) with Nr1d1 knockdown. The evidence for roles for opioid-related genes opioid receptor, delta 1 (Oprd1) and preproenkephalin (Penk) was also found. Together, these results suggest that Nr1d1 in the NAc modulates sociability and anxiety-related behaviour in a sex-specific manner, and circadian, histone deacetylase and opioid-related genes may be involved in the expression of these behavioural phenotypes.

Genomic variants in an inbred mouse model predict mania-like behaviors.

Contemporary rodent models for bipolar disorders split the bipolar spectrum into complimentary behavioral endophenotypes representing mania and depression. Widely accepted mania models typically utilize single gene transgenics or pharmacological manipulations, but inbred rodent strains show great potential as mania models. Their acceptance is often limited by the lack of genotypic data needed to establish construct validity. In this study, we used a unique strategy to inexpensively explore and confirm population allele differences in naturally occurring candidate variants in a manic rodent strain, the Madison (MSN) mouse strain. Variants were identified using whole exome resequencing on a small population of animals. Interesting candidate variants were confirmed in a larger population with genotyping. We enriched these results with observations of locomotor behavior from a previous study. Resequencing identified 447 structural variants that are mostly fixed in the MSN strain relative to control strains. After filtering and annotation, we found 11 non-synonymous MSN variants that we believe alter protein function. The allele frequencies for 6 of these variants were consistent with explanatory variants for the Madison strain's phenotype. The variants are in the Npas2, Cp, Polr3c, Smarca4, Trpv1, and Slc5a7 genes, and many of these genes' products are in pathways implicated in human bipolar disorders. Variants in Smarca4 and Polr3c together explained over 40% of the variance in locomotor behavior in the Hsd:ICR founder strain. These results enhance the MSN strain's construct validity and implicate altered nucleosome structure and transcriptional regulation as a chief molecular system underpinning behavior.

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.

Metabotropic glutamate receptor 3 is downregulated and its expression is shifted from neurons to astrocytes in the mouse lateral septum during the postpartum period.

The inhibitory metabotropic glutamate receptor 3 (mGluR3) plays diverse and complex roles in brain function, including synaptic plasticity and neurotransmission. We recently found that mGluR3 is downregulated in the lateral septum (LS) of postpartum females using microarray and qPCR analysis. In this study, we used double fluorescence immunohistochemical approaches to characterize mGluR3 changes in LS of the postpartum brain. The number of mGluR3-immunoractive cells was significantly reduced in the dorsal (LSD) and intermediate (LSI) but not ventral (LSV) parts of the LS in postpartum versus virgin females. mGluR3 immunoreactivity in the LS was found predominantly in neurons (~70%), with a smaller portion (~20%-30%) in astrocytes. Colocalization analysis revealed a reduced mGluR3 expression in neurons but an increased astrocytic localization in postpartum LSI. This change in the pattern of expression suggests that mGluR3 expression is shifted from neurons to astrocytes in postpartum LS, and the decrease in mGluR3 is neuron-specific. Because mGluR3 is inhibitory and negatively regulates glutamate and GABA release, decreases in neuronal expression would increase glutamate and GABA signaling. Given our recent finding that ~90% of LS neurons are GABAergic, the present data suggest that decreases in mGluR3 are a mechanism for elevated GABA in LS in the postpartum state.

Glutamate, GABA, and glutamine are synchronously upregulated in the mouse lateral septum during the postpartum period.

Dramatic structural and functional remodeling occurs in the postpartum brain for the establishment of maternal care, which is essential for the growth and development of young offspring. Glutamate and GABA signaling are critically important in modulating multiple behavioral performances. Large scale signaling changes occur in the postpartum brain, but it is still not clear to what extent the neurotransmitters glutamate and GABA change and whether the ratio of glutamate/GABA remains balanced. In this study, we examined the glutamate/GABA-glutamine cycle in the lateral septum (LS) of postpartum female mice. In postpartum females (relative to virgins), tissue levels of glutamate and GABA were elevated in LS and increased mRNA was found for the respective enzymes producing glutamate and GABA, glutaminase (Gls) and glutamate decarboxylase 1 and 2 (Gad1 and Gad2). The common precursor, glutamine, was elevated as was the enzyme that produces it, glutamate-ammonia ligase (Glul). Additionally, glutamate, GABA, and glutamine were positively correlated and the glutamate/GABA ratio was almost identical in the postpartum and virgin females. Collectively, these findings indicate that glutamate and GABA signaling are increased and that the ratio of glutamate/GABA is well balanced in the maternal LS. The postpartum brain may provide a useful model system for understanding how glutamate and GABA are linked despite large signaling changes. Given that some mental health disorders, including depression and schizophrenia display dysregulated glutamate/GABA ratio, and there is increased vulnerability to mental disorders in mothers, it is possible that these postpartum disorders emerge when glutamate and GABA changes are not properly coordinated.

Addiction and reward-related genes show altered expression in the postpartum nucleus accumbens.

Motherhood involves a switch in natural rewards, whereby offspring become highly rewarding. Nucleus accumbens (NAC) is a key CNS region for natural rewards and addictions, but to date no study has evaluated on a large scale the events in NAC that underlie the maternal change in natural rewards. In this study we utilized microarray and bioinformatics approaches to evaluate postpartum NAC gene expression changes in mice. Modular Single-set Enrichment Test (MSET) indicated that postpartum (relative to virgin) NAC gene expression profile was significantly enriched for genes related to addiction and reward in five of five independently curated databases (e.g., Malacards, Phenopedia). Over 100 addiction/reward related genes were identified and these included: Per1, Per2, Arc, Homer2, Creb1, Grm3, Fosb, Gabrb3, Adra2a, Ntrk2, Cry1, Penk, Cartpt, Adcy1, Npy1r, Htr1a, Drd1a, Gria1, and Pdyn. ToppCluster analysis found maternal NAC expression profile to be significantly enriched for genes related to the drug action of nicotine, ketamine, and dronabinol. Pathway analysis indicated postpartum NAC as enriched for RNA processing, CNS development/differentiation, and transcriptional regulation. Weighted Gene Coexpression Network Analysis (WGCNA) identified possible networks for transcription factors, including Nr1d1, Per2, Fosb, Egr1, and Nr4a1. The postpartum state involves increased risk for mental health disorders and MSET analysis indicated postpartum NAC to be enriched for genes related to depression, bipolar disorder (BPD), and schizophrenia. Mental health related genes included: Fabp7, Grm3, Penk, and Nr1d1. We confirmed via quantitative PCR Nr1d1, Per2, Grm3, Penk, Drd1a, and Pdyn. This study indicates for the first time that postpartum NAC involves large scale gene expression alterations linked to addiction and reward. Because the postpartum state also involves decreased response to drugs, the findings could provide insights into how to mitigate addictions.

Medial prefrontal cortex: genes linked to bipolar disorder and schizophrenia have altered expression in the highly social maternal phenotype.

The transition to motherhood involves CNS changes that modify sociability and affective state. However, these changes also put females at risk for post-partum depression and psychosis, which impairs parenting abilities and adversely affects children. Thus, changes in expression and interactions in a core subset of genes may be critical for emergence of a healthy maternal phenotype, but inappropriate changes of the same genes could put women at risk for post-partum disorders. This study evaluated microarray gene expression changes in medial prefrontal cortex (mPFC), a region implicated in both maternal behavior and psychiatric disorders. Post-partum mice were compared to virgin controls housed with females and isolated for identical durations. Using the Modular Single-set Enrichment Test (MSET), we found that the genetic landscape of maternal mPFC bears statistical similarity to gene databases associated with schizophrenia (5 of 5 sets) and bipolar disorder (BPD, 3 of 3 sets). In contrast to previous studies of maternal lateral septum (LS) and medial preoptic area (MPOA), enrichment of autism and depression-linked genes was not significant (2 of 9 sets, 0 of 4 sets). Among genes linked to multiple disorders were fatty acid binding protein 7 (Fabp7), glutamate metabotropic receptor 3 (Grm3), platelet derived growth factor, beta polypeptide (Pdgfrb), and nuclear receptor subfamily 1, group D, member 1 (Nr1d1). RT-qPCR confirmed these gene changes as well as FMS-like tyrosine kinase 1 (Flt1) and proenkephalin (Penk). Systems-level methods revealed involvement of developmental gene networks in establishing the maternal phenotype and indirectly suggested a role for numerous microRNAs and transcription factors in mediating expression changes. Together, this study suggests that a subset of genes involved in shaping the healthy maternal brain may also be dysregulated in mental health disorders and put females at risk for post-partum psychosis with aspects of schizophrenia and BPD.