BXD Recombinant Inbred Mice Research Articles

Systems genetics identifies methionine as a high risk factor for Alzheimer's disease.

As a dietary strategy, methionine restriction has been reported to promote longevity and regulate metabolic disorders. However, the role and possible regulatory mechanisms underlying methionine in neurodegenerative diseases such as Alzheimer's disease (AD), remain unexplored. This study utilized the data from BXD recombinant inbred (RI) mice to establish a correlation between the AD phenotype in mice and methionine level. Gene enrichment analysis indicated that the genes associated with the concentration of methionine in the midbrain are involved in the dopaminergic synaptic signaling pathway. Protein interaction network analysis revealed that glycogen synthase kinase 3 beta (GSK-3β) was a key regulator of the dopaminergic synaptic pathway and its expression level was significantly correlated with the AD phenotype. Finally, in vitro experiments demonstrated that methionine deprivation could reduce the expression of Aβ and phosphorylated Tau, suggesting that lowering methionine levels in humans may be a preventive or therapeutic strategy for AD. In conclusion, our findings support that methionine is a high risk factor for AD. These findings predict potential regulatory network, theoretically supporting methionine restriction to prevent AD.

Ninein, a candidate gene for ethanol anxiolysis, shows complex exon-specific expression and alternative splicing differences between C57BL/6J and DBA/2J mice.

Ethanol's anxiolytic actions contribute to increased consumption and the development of Alcohol Use Disorder (AUD). Our laboratory previously identified genetic loci contributing to the anxiolytic-like properties of ethanol in BXD recombinant inbred mice, derived from C57BL/6J (B6) and DBA/2J (D2) progenitor strains. That work identified Ninein (Nin) as a candidate gene underlying ethanol's acute anxiolytic-like properties in BXD mice. Nin has a complex exonic content with known alternative splicing events that alter cellular distribution of the NIN protein. We hypothesize that strain-specific differences in Nin alternative splicing contribute to changes in Nin gene expression and B6/D2 strain differences in ethanol anxiolysis. Using quantitative reverse-transcriptase PCR to target specific Nin splice variants, we identified isoform-specific exon expression differences between B6 and D2 mice in prefrontal cortex, nucleus accumbens and amygdala. We extended this analysis using deep RNA sequencing in B6 and D2 nucleus accumbens samples and found that total Nin expression was significantly higher in D2 mice. Furthermore, exon utilization and alternative splicing analyses identified eight differentially utilized exons and significant exon-skipping events between the strains, including three novel splicing events in the 3' end of the Nin gene that were specific to the D2 strain. Additionally, we document multiple single nucleotide polymorphisms in D2 Nin exons that are predicted to have deleterious effects on protein function. Our studies provide the first in-depth analysis of Nin alternative splicing in brain and identify a potential genetic mechanism altering Nin expression and function between B6 and D2 mice, thus possibly contributing to differences in the anxiolytic-like properties of ethanol between these strains. This work adds novel information to our understanding of genetic differences modulating ethanol actions on anxiety that may contribute to the risk for alcohol use disorder.

The ameliorative effects of choline on ethanol-induced cell death in the neural tube of susceptible BXD strains of mice.

Fetal alcohol spectrum disorders (FASD) are the leading preventable cause of intellectual disability, providing the impetus for evaluating various potential treatments to ameliorate ethanol's teratogenic effects, particularly in the nervous system. One treatment is the dietary supplement choline which has been shown to mitigate at least some of ethanol's teratogenic effects. The present study was designed to investigate the effects of genetics on choline's efficacy in ameliorating cell death in the developing neural tube. Previously, we examined BXD recombinant inbred mice, and their parental C57BL/6 J (B6) and DBA/2 J strains, and identified strains that were sensitive to ethanol's teratogenic actions. Thus, we used these strains to identify response to choline treatment. Timed pregnant mice from 4 strains (B6, BXD51, BXD73, BXD2) were given either ethanol or isocaloric maltose-dextrin (5.8 g/kg in two administrations separated by 2 h) with choline at one of 3 doses: 0, 100 or 250 mg/kg. Subjects were exposed via intragastric gavage on embryonic day 9 and embryos were collected 7 h after the initial ethanol administrations. Cell death was analyzed using TUNEL staining in the developing forebrain and brainstem. Choline ameliorated the ethanol-induced cell death across all 4 strains without causing enhanced cell death in control mice. Choline was effective in both the developing telencephalon and in the brainstem. Both doses diminished cell death, with some differences across strains and brain regions, although the 100 mg/kg dose was most consistent in mitigating ethanol-related cell death. Comparisons across strains showed that there was an effect of strain, particularly in the forebrain at the higher dose. These results show that choline is effective in ameliorating ethanol-induced cell death at this early stage of nervous system development. However, there were some strain differences in its efficacy, especially at the high dose, providing further evidence of the importance of genetics in influencing the ability of choline to protect against prenatal alcohol exposure.

Expression Levels of the Tnni3k Gene in the Heart Are Highly Associated with Cardiac and Glucose Metabolism-Related Phenotypes and Functional Pathways.

Troponin-I interacting kinase encoded by the TNNI3K gene is expressed in nuclei and Z-discs of cardiomyocytes. Mutations in TNNI3K were identified in patients with cardiac conduction diseases, arrhythmias, and cardiomyopathy. We performed cardiac gene expression, whole genome sequencing (WGS), and cardiac function analysis in 40 strains of BXD recombinant inbred mice derived from C57BL/6J (B6) and DBA/2J (D2) strains. Expression quantitative trait loci (eQTLs) mapping and gene enrichment analysis was performed, followed by validation of candidate Tnni3k-regulatory genes. WGS identified compound splicing and missense T659I Tnni3k variants in the D2 parent and some BXD strains (D allele) and these strains had significantly lower Tnni3k expression than those carrying wild-type Tnni3k (B allele). Expression levels of Tnni3k significantly correlated with multiple cardiac (heart rate, wall thickness, PR duration, and T amplitude) and metabolic (glucose levels and insulin resistance) phenotypes in BXDs. A significant cis-eQTL on chromosome 3 was identified for the regulation of Tnni3k expression. Furthermore, Tnni3k-correlated genes were primarily involved in cardiac and glucose metabolism-related functions and pathways. Genes Nodal, Gnas, Nfkb1, Bmpr2, Bmp7, Smad7, Acvr1b, Acvr2b, Chrd, Tgfb3, Irs1, and Ppp1cb were differentially expressed between the B and D alleles. Compound splicing and T659I Tnni3k variants reduce cardiac Tnni3k expression and Tnni3k levels are associated with cardiac and glucose metabolism-related phenotypes.

Abstract 11608: Differences in Genetic Variants and Common Genetic Regulators between Arrhythmogenic Cardiomyopathy and Arrhythmogenic Dilated Cardiomyopathy Identified Using Systems Genetics Approaches

Introduction: Prominent arrhythmias are features of arrhythmogenic cardiomyopathy (ACM) and arrhythmogenic dilated cardiomyopathy (aDCM). Despite clinicopathological differences, these diseases share common genetic causes. Hypothesis: This study aimed to identify variants in ACM and aDCM and map the candidate genetic regulators using a mouse genetic reference population (GRP). Methods: Whole exome sequencing was performed on 226 individuals with ACM and 102 patients with aDCM. Bioinformatics analysis was applied to identify the pathogenic or likely pathogenic variants followed by direct sequencing confirmation. Expression quantitative trait loci (eQTLs) mapping to heart transcriptomes of GRP of BXD recombinant inbred (RI) mice was conducted in GeneNetwork (http://www.genenetwork.org). Results: In the ACM cohort, a total of 82 variants (MAF<0.01) in 24 genes have been confirmed being positive in 42 probands. Depending on the frequency of the identified variants, the top 10 genes (excluding TTN) with the most pathogenic variants were PKP2, MYH6, FLNC, RYR2, MYH7B, PLEC, JPH2, SYNE1, NRAP and POSTN. In the aDCM cohort, 141 pathogenic variants have been identified and 91 patients harbored at least 1 deleterious variant. The top 10 genes were RYR2, FLNC, MYH7, DSP, SCN5A,MYPN, NEBL, CTNNA3, ACAP6 and SLC22A5. We analyzed the distribution of the first principal component (PC1) of mRNA expression levels of those 10 genes (eigentraits) from each cohort in the heart across BXD strains. Interval mapping for PC1 eigentraits for ACM identified a significant eQTL (P<0.05) in the mouse genome at 24-36 Mb of chromosome (Chr) 6. In aDCM, interval mapping for PC1 eigentraits identified 3 suggestive eQTLs at Chr 5 (44-53 Mb), Chr 6 (32-34 Mb) and Chr 15 (82-84 Mb). Further gene exploration of these eQTLs suggested that Mtpn, Lep, Cyb5r3, Tspo, Slit2, and Ppargc1a could be common regulators that interact with those top ACM and aDCM genes. Conclusions: This study identified diverse genes in ACM and aDCM cases. Our joint human and mouse analysis using a systems genetics approach identified genetic regulators that could affected expression of genes with pathogenic variants in ACM and aDCM patients.

Multi-omics analysis reveals expression complexity and functional diversity of mouse kinome.

Protein kinases are a crucial component of signaling pathways involved in a wide range of cellular responses, including growth, proliferation, differentiation, and migration. Systematic investigation of protein kinases is critical to better understand phosphorylation-mediated signaling pathways and may provide insights into the development of potential therapeutic drug targets. Here we perform a systems-level analysis of the mouse kinome by analyzing multi-omics data. We used bulk and single-cell transcriptomic data from the C57BL/6J mouse strain to define tissue- and cell-type-specific expression of protein kinases, followed by investigating variations in sequence and expression between C57BL/6J and DBA/2J strains. We then profiled a deep brain phosphoproteome from C57BL/6J and DBA/2J strains as well as their reciprocal hybrids to infer the activity of the mouse kinome. Finally, we performed phenome-wide association analysis using the BXD recombinant inbred (RI) mice (a cross between C57BL/6J and DBA/2J strains) to identify any associations between variants in protein kinases and phenotypes. Collectively, our study provides a comprehensive analysis of the mouse kinome by investigating genetic sequence variation, tissue-specific expression patterns, and associations with downstream phenotypes.

Hypothalamic Irak4 is a genetically controlled regulator of hypoglycemia-induced glucagon secretion

ObjectivesGlucagon secretion to stimulate hepatic glucose production is the first line of defense against hypoglycemia. This response is triggered by so far incompletely characterized central hypoglycemia-sensing mechanisms, which control autonomous nervous activity and hormone secretion. The objective of this study was to identify novel hypothalamic genes controlling insulin-induced glucagon secretion. MethodsTo obtain new information on the mechanisms of hypothalamic hypoglycemia sensing, we combined genetic and transcriptomic analysis of glucagon response to insulin-induced hypoglycemia in a panel of BXD recombinant inbred mice. ResultsWe identified two QTLs on chromosome 8 and chromosome 15. We further investigated the role of Irak4 and Cpne8, both located in the QTL on chromosome 15, in C57BL/6J and DBA/2J mice, the BXD mouse parental strains. We found that the poor glucagon response of DBA/2J mice was associated with higher hypothalamic expression of Irak4, which encodes a kinase acting downstream of the interleukin-1 receptor (Il-1R), and of Il-ß when compared with C57BL/6J mice. We showed that intracerebroventricular administration of an Il-1R antagonist in DBA/2J mice restored insulin-induced glucagon secretion; this was associated with increased c-fos expression in the arcuate and paraventricular nuclei of the hypothalamus and with higher activation of both branches of the autonomous nervous system. Whole body inactivation of Cpne8, which encodes a Ca++-dependent regulator of membrane trafficking and exocytosis, however, had no impact on insulin-induced glucagon secretion. ConclusionsCollectively, our data identify Irak4 as a genetically controlled regulator of hypoglycemia-activated hypothalamic neurons and glucagon secretion.

A systems genetics approach delineates the role of Bcl2 in leukemia pathogenesis

Leukemia is a group of malignancies of the blood forming tissues, and is characterized by the uncontrolled proliferation of blood cells. In the United States, it accounts for approximately 3.5% and 4% of all cancer-related incidences and mortalities, respectively. The current study aimed to explore the role of Bcl2 and associated genes in leukemia pathogenesis using a systems genetics approach. The transcriptome data from BXD Recombinant Inbred (RI) mice was analyzed to identify the expression of Bcl2 in myeloid cells. eQTL mapping was performed to select the potential chromosomal region and subsequently identify the candidate gene modulating the expression of Bcl2. Furthermore, gene enrichment and protein-protein interaction (PPI) analyses of the Bcl2-coexpressed genes were performed to demonstrate the role of Bcl2 in leukemia pathogenesis. The Bcl2-coexpressed genes were found to be enriched in various hematopoietic system related functions, and multiple pathways related to signaling, immune response, and cancer. The PPI network analysis demonstrated direct interaction of hematopoietic function related genes, such as Bag3, Bak1, Bcl2l11, Bmf, Mapk9, Myc, Ppp2r5c, and Ppp3ca with Bcl2. The eQTL mapping identified a 4.5 Mb genomic region on chromosome 11, potentially regulating the expression of Bcl2. A multi-criteria filtering process identified Top2a, among the genes located in the mapped locus, as the best candidate upstream regulator for Bcl2 expression variation. Hence, the current study provides better insights into the role of Bcl2 in leukemia pathogenesis and demonstrates the significance of our approach in gaining new knowledge on leukemia. Furthermore, our findings from the PPI network analysis and eQTL mapping provide supporting evidence of leukemia-associated genes, which can be further explored for their functional importance in leukemia. Data availabilityThe myeloid cell transcriptomic data of the BXD mice used in this study can be accessed through our GeneNetwork (http://www.genenetwork.org) with the accession number of GN144.

Genetic dissection of glutathione S-transferase omega-1: identification of novel downstream targets and Alzheimer's disease pathways.

Alzheimer's disease (AD) is affected by genetic factors. Polymorphisms in the glutathione S-transferase omega-1 (Gsto1) gene have been shown by genetic correlation analyses performed in different ethnic populations to be genetic risk factors for AD. Gene expression profile data from BXD recombinant inbred mice were used in combination with genetic and bioinformatic analyses to characterize the mechanisms underlying regulation of Gsto1 variation regulation and to identify network members that may contribute to AD risk or progression. Allele-specific assays confirmed that variation in Gsto1 expression is controlled by cis-expression quantitative trait loci. We found that Gsto1 mRNA levels were related to several central nervous system traits, such as glial acidic fibrillary protein levels in the caudate putamen, cortical gray matter volume, and hippocampus mossy fiber pathway volume. We identified 2168 genes whose expression was highly correlated with that of Gsto1. Some genes were enriched for the most common neurodegenerative diseases. Some Gsto1-related genes identified in this study had previously been identified as susceptibility genes for AD, such as APP, Grin2b, Ide, and Psenen. To evaluate the relationships between Gsto1 and candidate network members, we transfected astrocytes with Gsto1 siRNA and assessed the effect on putative downstream effectors. We confirmed that knockdown of Gsto1 had a significant influence on Pa2g4 expression, suggesting that Pa2g4 may be a downstream effector of Gsto1, and that both genes interact with other genes in a network during AD pathogenesis.

BXD Recombinant Inbred Mice as a Model to Study Neurotoxicity.

BXD recombinant inbred (RI) lines represent a genetic reference population derived from a cross between C57BL/6J mice (B6) and DBA/2J mice (D2), which through meiotic recombination events possesses recombinant chromosomes containing B6 or D2 haplotype segments. The quantitative trait loci (QTLs) are the locations of segregating genetic polymorphisms and are fundamental to understanding genetic diversity in human disease susceptibility and severity. QTL mapping represents the typical approach for identifying naturally occurring polymorphisms that influence complex phenotypes. In this process, genotypic values at markers of known genomic locations are associated with phenotypic values measured in a segregating population. Indeed, BXD RI strains provide a powerful tool to study neurotoxicity induced by different substances. In this review, we describe the use of BXD RI lines to understand the underlying mechanisms of neurotoxicity in response to ethanol and cocaine, as well as metals and pesticide exposures.

The BXD21/TyJ recombinant inbred strain as a model for innate inflammatory response in distinct brain regions

Oxidative stress and inflammatory cytokines affect the human brain, increasing the risk for mood and cognitive disorders. Such risk might be selective to brain-specific regions. Here, we determined whether BXD recombinant inbred (RI) mice strains are more suitable than C57BL/6J mice for the understanding of the relationship between antioxidant response and inflammatory responses. We hypothesized that inflammatory responses could be independent of antioxidant response and be inherent to brain-specific regions. This hypothesis will be addressed by the analyses of mRNA expression. We explored, at 7-months-of-age, the innate activation of proinflammatory cytokines (tumor necrosis factor alpha (TNFα) and interleukin 6 (IL-6), as well as Kelch-like ECH-associating protein 1 (Keap1), nuclear factor erythroid 2 related factor 2 (Nrf2) and glutathione peroxidase 1 (Gpx1) mRNA in both male and female BXD84/RwwJ RI, BXD21/TyJ RI and control strain (C57BL/6J mice). We report that: (1) The cerebellum is more sensitive to antioxidant response in the BXD21/TyJ RI strain; (2) The cerebellum, hippocampus and striatum show increased levels of cytokines in the BXD21/TyJ RI strain; (3) The BXD RI strain has lower brain weight relative to control strain (C57BL/6 mice). In conclusion, our novel data show the utility of the BXD21/TyJ RI strain mice in offering mechanistic insight into Nrf2’s role in the inflammatory system.

Cardiac phenotyping in the BXD strains as a genetic reference population for cardiovascular diseases

IntroductionDiscovery of genetic modifiers in cardiomyopathy (CM) has proved difficult in human cohorts because the wide range of clinical manifestations due to joint effects of causal and multiple modifier genes and environment. Therefore, combining human GWAS findings with additional animal model and molecular data using systems genetics is powerful approach to understand complex diseases by focusing on how genes work together in groups rather than singly. The main challenge for a systems genetics analysis is the need for comparatively large sample sizes drawn from a genetic reference population (GRP) and the use of advanced statistical and computational methods and models. The study was aimed to establish a murine GRP for CM systems genetics analysis. The BXD recombinant inbred (RI) strains derived from C57BL/6J and DBA/2J (D2) mice is the largest murine GRP. The parental D2 strain expresses CM key features, which make BXD strains particularly suited for systems genetics analysis for CMs.HypothesisBXD strains can be used to study genetic influences and modifiers in CMs.MethodsCM‐associated traits (heart function, rhythm and mass, cardiomyocyte size, fibrosis) and heart transcriptome were assessed in male and female BXD RI mice (40 lines, N>5) at 5–6 months of age. Echocardiography, ECG and blood pressure was performed followed by cardiac histology and microarray analysis (Affymetrix Mouse Gene 2.0 ST).ResultsEchocardiography revealed more than two‐fold variance in LV end‐systolic and end‐diastolic volumes among tested BXD strains. BXD155 has the largest end‐systolic LV volume (42.8±9.0 μl) and BXD74 has the smallest volume (12.1±4.8 μl). Accordingly, %EF also varied greatly among BXDs. Severe hypertension was found in many BXD strains with highest systolic and diastolic blood pressure (200/159 mmHg) in D2. Atrioventricular (AV) block, bradycardia, ST elevation, and inverted T wave was revealed in several BXD lines. The lowest heart rate was found in BXD171 mice (369±20 beat/min). Heart weight, cardiomyocyte size, and cardiac collagen volume fraction significantly varied among BXD strains as well. D2 strain has the largest hearts (185.1 mg) and B6 has the smallest hearts (129.8 mg), while the average heart weight among BXD strains is 144.7mg. Several CM causal and candidate genes such as Mybpc3, Prkag2, Ankrd1, Myl3, and Mypn have mutation in D2 mouse or they are trans‐regulated by the gene with mutation in D2.ConclusionsCM phenotypes from D2 mouse are segregated among the BXD family of strains and CM‐associated traits vary among BXD strains, demonstrating that the BXD family as an excellent platform to identify genetic factors and modifiers that influence CM phenotypes. Identification of these factors will facilitate improved understanding of gene modifications and phenotypes in human CMs.This abstract is from the Experimental Biology 2019 Meeting. There is no full text article associated with this abstract published in The FASEB Journal.

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

BackgroundGene co-expression studies can provide important insights into molecular and cellular signaling pathways. The GeneNetwork database is a unique resource for co-expression analysis using data from a variety of tissues across genetically distinct inbred mice. However, extraction of biologically meaningful co-expressed gene sets is challenging due to variability in microarray platforms, probe quality, normalization methods, and confounding biological factors. In this study, we tested whether literature derived functional cohesion could be used as an objective metric in lieu of ‘ground truth’ to evaluate the quality of probes and microarray datasets.ResultsWe examined Sirtuin-3 (Sirt3) co-expressed gene sets extracted from either liver or brain tissues of BXD recombinant inbred mice in the GeneNetwork database. Depending on the microarray platform, there were as many as 26 probes that targeted different regions of Sirt3 primary transcript. Co-expressed gene sets (ranging from 100–1000 genes) associated with each Sirt3 probe were evaluated using the previously developed literature-derived cohesion p-value (LPv) and benchmarked against ‘gold standards’ derived from proteomic studies or Gene Ontology classifications. We found that the maximal F-measure was obtained at an average window size of 535 genes. Using set size of 500 genes, the Pearson correlations between LPv and F-measure as well as between LPv and mitochondrial gene enrichment p-values were 0.90 and 0.93, respectively. Importantly, we found that the LPv approach can distinguish high quality Sirt3 probes. Analysis of the most functionally cohesive Sirt3 co-expressed gene set revealed core metabolic pathways that were shared between hippocampus and liver as well as distinct pathways which were unique to each tissue. These results are consistent with other studies that suggest Sirt3 is a key metabolic regulator and has distinct functions in energy-producing vs. energy-demanding tissues.ConclusionsOur results provide proof-of-concept that literature cohesion analysis is useful for evaluating the quality of probes and microarray datasets, particularly when experimentally derived gold standards are unavailable. Our approach would enable researchers to rapidly identify biologically meaningful co-expressed gene sets and facilitate discovery from high throughput genomic data.

Genetic Influences on the Amount of Cell Death in the Neural Tube of BXD Mice Exposed to Acute Ethanol at Midgestation.

Fetal alcohol spectrum disorders (FASD) have a strong genetic component although the genes that underlie this are only beginning to be elucidated. In the present study, one of the most common phenotypes of FASD, cell death within the early developing neural tube, was examined across a genetic reference population in a reverse genetics paradigm with the goal of identifying genetic loci that could influence ethanol (EtOH)-induced apoptosis in the early developing neural tube. BXD recombinant inbred mice as well as the parental strains were used to evaluate genetic differences in EtOH-induced cell death after exposure on embryonic day 9.5. Dams were given either 5.8g/kg EtOH or isocaloric maltose-dextrin in 2 doses via intragastric gavage. Embryos were collected 7hours after the initial exposure and cell death evaluated via TUNEL staining in the brainstem and forebrain. Genetic loci were evaluated using quantitative trait locus (QTL) analysis at GeneNetwork.org. Significant strain differences were observed in the levels of EtOH-induced cell death that were due to genetic effects and not confounding variables such as differences in developmental maturity or cell death kinetics. Comparisons between the 2 regions of the developing neural tube showed little genetic correlation with the QTL maps exhibiting no overlap. Significant QTLs were found on murine mid-chromosome 4 and mid-chromosome 14 only in the brainstem. Within these chromosomal loci, a number of interesting candidate genes were identified that could mediate this differential sensitivity including Nfia (nuclear factor I/A) and Otx2 (orthodenticle homeobox 2). These studies demonstrate that the levels of EtOH-induced cell death occur in strain- and region-dependent manners. Novel QTLs on mouse Chr4 and Chr14 were identified that modulate the differential sensitivity to EtOH-induced apoptosis in the embryonic brainstem. The genes underlying these QTLs could identify novel molecular pathways that are critical in this phenotype.

An Application of the Bayesian Periodicity Test to Identify Diurnal Rhythm Genes in the Brain.

Biological systems are extremely dynamic and many aspects of cellular processes show rhythmic circadian patterns. Extracting such information from large expression data is challenging. In this work, we present a modified application of the Empirical Bayes periodicity test to identify genes with diurnal rhythmic behavior in two brain regions. The hypothalamus and amygdala gene expression data were generated from 100 BXD recombinant inbred mice during the day hours. Brain samples were collected over the course of two days. We first filtered the transcripts based on rank correlation at matched time points between day-1 and day-2. We then applied the proposed test of periodicity to identify diurnal rhythm genes in the full cohort and gender-specific sub-cohorts. In hypothalamus, at a Benjamini-Hochberg false discovery rate (BH-FDR) of 0.01, we identified 15 transcripts with cyclic behavior in the full cohort, none, and 53 transcripts in the female and male cohort, respectively. Similarly, in amygdala, we identified 58 diurnal rhythm genes in the full cohort, and 1 and 28 in the female and male cohorts, respectively. In conclusion, we present a modified version of the empirical Bayes periodicity test to detect periodic expression patterns. Our results demonstrate that this approach can capture cyclic patterns from relatively noisy expression data sets.

Divergence in Morris Water Maze-Based Cognitive Performance under Chronic Stress Is Associated with the Hippocampal Whole Transcriptomic Modification in Mice.

Individual susceptibility determines the magnitude of stress effects on cognitive function. The hippocampus, a brain region of memory consolidation, is vulnerable to stressful environments, and the impact of stress on hippocampus may determine individual variability in cognitive performance. Therefore, the purpose of this study was to define the relationship between the divergence in spatial memory performance under chronically unpredictable stress and an associated transcriptomic alternation in hippocampus, the brain region of spatial memory consolidation. Multiple strains of BXD (B6 × D2) recombinant inbred mice went through a 4-week chronic variable stress (CVS) paradigm, and the Morris water maze (MWM) test was conducted during the last week of CVS to assess hippocampal-dependent spatial memory performance and grouped animals into low and high performing groups based on the cognitive performance. Using hippocampal whole transcriptome RNA-sequencing data, differential expression, PANTHER analysis, WGCNA, Ingenuity's upstream regulator analysis in the Ingenuity Pathway Analysis® and phenotype association analysis were conducted. Our data identified multiple genes and pathways that were significantly associated with chronic stress-associated cognitive modification and the divergence in hippocampal dependent memory performance under chronic stress. Biological pathways associated with memory performance following chronic stress included metabolism, neurotransmitter and receptor regulation, immune response and cellular process. The Ingenuity's upstream regulator analysis identified 247 upstream transcriptional regulators from 16 different molecule types. Transcripts predictive of cognitive performance under high stress included genes that are associated with a high occurrence of Alzheimer's and cognitive impairments (e.g., Ncl, Eno1, Scn9a, Slc19a3, Ncstn, Fos, Eif4h, Copa, etc.). Our results show that the variable effects of chronic stress on the hippocampal transcriptome are related to the ability to complete the MWM task and that the modulations of specific pathways are indicative of hippocampal dependent memory performance. Thus, the divergence in spatial memory performance following chronic stress is related to the unique pattern of gene expression within the hippocampus.

Dissection of Z‐disc myopalladin gene network involved in the development of restrictive cardiomyopathy using system genetics approach

BackgroundGenetic differences mediate individual differences in susceptibility to cardiomyopathy and severity of disease symptoms. Myopalladin (Mypn) is one of the genes associated with familial cardiomyopathies including dilated, hypertrophic and restrictive cardiomyopathy (RCM). Mutant knock‐in mice carrying the human Q529X‐MYPN mutation recapitulate the RCM phenotype the same as seen in the patients. The aim of this study was to investigate the regulation of Mypn and identify gene networks involved in RCM development.MethodsGene expression values were measured in the heart of a large family of BXD recombinant inbred (RI) mice derived from C57BL/6J and DBA/2J. The proteomic data was collected from Mypn knock‐in and knock‐out mice. Expression quantitative trait locus (eQTL) mapping methods and gene enrichment analysis were used to identify Mypn gene pathway and evaluate co‐expression networks.ResultsWe found a wide range of expression variation in Mypn among the large cohort of BXD RI strains. Two upstream genetic loci, modulators of the expression of Mypn, have been identified on chromosomes 1 and 5. Candidate genes within those loci include Ncoa2, Vcpip1, Sgk3, and Lgi2. We also identified 15 sarcomeric genes that interact with Mypn, and constructed the Mypn gene network. Among them, well‐known cardiomyopathy‐associated genes such as Ldb3, Des, Actn2, Fhod3, Tpm2, and Syne1 were found to be members of this network. Several novel genes (Myo18b, Fhod3, Myom1, Bmp10, Myl4, Obscn, and Pdlim5) were also related to the Mypn network. In particular, two members (Syne1 and Myom1) have been confirmed at the protein level as high chance genes working as modifiers to induce cardiomyopathy.ConclusionsThrough a systematic genetic strategy we found several novel genes interacting with Mypn that may play an important role in the development of RCM. This work demonstrates that a systematic genetics approach is a powerful tool for constructing co‐expression networks that define the biological process categories within which similarly regulated genes function.Support or Funding InformationR01HL128350 (Lu L), R01 HL53392 and R01 HL087000 (Towbin JA)

BXD recombinant inbred strains participate in social preference, anxiety and depression behaviors along sex-differences in cytokines and tactile allodynia

Depression and anxiety are the most common psychiatric disorders, representing a major public health concern. Dysregulation of oxidative and inflammatory systems may be associated with psychiatric disorders, such as depression and anxiety. Due to the need to find appropriate animal models to the understanding of such disorders, we queried whether 2 BXD recombinant inbred (RI) mice strains (BXD21/TyJ RI and BXD84/RwwJ RI mice) and C57BL/6 wild-type mice show differential performance in depression and anxiety related behaviors and biomarkers. Specifically, we assessed social preference, elevated plus maze, forced swim, and Von Frey tests at 3–4 months-of-age, as well as activation of cytokines and antioxidant mRNA levels in the cortex at 7 months-of-age. We report that (1) the BXD84/RwwJ RI strain exhibits anxiety disorder and social avoidance-like behavior (2) BXD21/TyJ RI strain shows a resistance to depression illness, and (3) sex-dependent cytokine profiles and allodynia with elevated inflammatory activity were inherent to male BXD21/TyJ RI mice. In conclusion, we provide novel data in favor of the use of BXD recombinant inbred mice to further understand anxiety and depression disorders.

Dissection of Z-disc myopalladin gene network involved in the development of restrictive cardiomyopathy using system genetics approach.

AIMTo investigate the regulation of Myopalladin (Mypn) and identify its gene network involved in restrictive cardiomyopathy (RCM).METHODSGene expression values were measured in the heart of a large family of BXD recombinant inbred (RI) mice derived from C57BL/6J and DBA/2J. The proteomics data were collected from Mypn knock-in and knock-out mice. Expression quantitative trait locus (eQTL) mapping methods and gene enrichment analysis were used to identify Mypn regulation, gene pathway and co-expression networks.RESULTSA wide range of variation was found in expression of Mypn among BXD strains. We identified upstream genetic loci at chromosome 1 and 5 that modulate the expression of Mypn. Candidate genes within these loci include Ncoa2, Vcpip1, Sgk3, and Lgi2. We also identified 15 sarcomeric genes interacting with Mypn and constructed the gene network. Two novel members of this network (Syne1 and Myom1) have been confirmed at the protein level. Several members in this network are already known to relate to cardiomyopathy with some novel genes candidates that could be involved in RCM.CONCLUSIONUsing systematic genetics approach, we constructed Mypn co-expression networks that define the biological process categories within which similarly regulated genes function. Through this strategy we have found several novel genes that interact with Mypn that may play an important role in the development of RCM.

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

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