Unraveling the molecular mechanisms of Ace2-mediated post-COVID-19 cognitive dysfunction through systems genetics approach.

The dysregulation of Angiotensin-converting enzyme 2 (ACE2) in central nervous system is believed associates with COVID-19 induced cognitive dysfunction. However, the detailed mechanism remains largely unknown. In this study, we performed a comprehensive system genetics analysis on hippocampal ACE2 based on BXD mice panel. Expression quantitative trait loci (eQTLs) mapping showed that Ace2 was strongly trans-regulated, and the elevation of Ace2 expression level was significantly correlated with impaired cognitive functions. Further Gene co-expression analysis showed that Ace2 may be correlated with the membrane proteins in Calcium signaling pathway. Further, qRT-PCR confirmed that SARS-CoV-2 spike S1 protein upregulated ACE2 expression together with eight membrane proteins in Calcium Signaling pathway. Moreover, such elevation can be attenuated by recombinant ACE2. Collectively, our findings revealed a potential mechanism of Ace2 in cognitive dysfunction, which could be beneficial for COVID-19-induced cognitive dysfunction prevention and potential treatment.

Systems genetics of influenza A virus-infected mice identifies TRIM21 as a critical regulator of pulmonary innate immune response.

Tripartite motif 21 (TRIM21) is a cytosolic Fc receptor that targets antibody-bound, internalized pathogens for destruction. Apart from this intrinsic defense role, TRIM21 is implicated in autoimmune diseases, inflammation, and autophagy. Whether TRIM21 participates in host interactions with influenza A virus (IAV), however, is unknown. By computational modeling of body weight and lung transcriptome data from the BXD parents (C57BL/6 J (B6) and DBA/2 J (D2)) and 41 BXD mouse strains challenged by IAV, we reveal that a Trim21-associated gene network modulates the early host responses to IAV infection. Trim21 transcripts were significantly upregulated in infected mice of both B6 and D2 backgrounds. Its expression was significantly higher in infected D2 than in infected B6 early after infection and significantly correlated with body weight loss. We identified significant trans-eQTL on chromosome 14 that regulates Trim21 expression. Nr1d2 and Il3ra were among the strongest candidate genes. Pathway analysis found Trim21 to be involved in inflammation and immunity related pathways, such as inflammation signaling pathways (TNF, IL-17, and NF-κB), viral detection signaling pathways (NOD-like and RIG-I-like), influenza, and other respiratory viral infections. Knockdown of TRIM21 in human lung epithelial A549 cells significantly augmented IAV-induced expression of IFNB1, IFNL1, CCL5, CXCL10, and IFN-stimulated genes including DDX58 and IFIH1, among others. Our data suggest that a TRIM21-associated gene network is involved in several aspects of inflammation and viral detection mechanisms during IAV infection. We identify and validate TRIM21 as a critical regulator of innate immune responses to IAV in human lung epithelial cells.

TEAD1, MYO7A and NDUFC2 are novel functional genes associated with glucose metabolism in BXD recombinant inbred population.

AIM: The liver is an important metabolic organ that governs glucolipid metabolism, and its dysfunction may cause non-alcoholic fatty liver disease, type 2 diabetes mellitus, dyslipidaemia, etc. We aimed to systematic investigate the key factors related to hepatic glucose metabolism, which may be beneficial for understanding the underlying pathogenic mechanisms for obesity and diabetes mellitus. MATERIALS AND METHODS: Oral glucose tolerance test (OGTT) phenotypes and liver transcriptomes of BXD mice under chow and high-fat diet conditions were collected from GeneNetwork. QTL mapping was conducted to pinpoint genomic regions associated with glucose homeostasis. Candidate genes were further nominated using a multi-criteria approach and validated to confirm their functional relevance in vitro. RESULTS: Our results demonstrated that plasma glucose levels in OGTT were significantly affected by both diet and genetic background, with six genetic regulating loci were mapped on chromosomes 1, 4, and 7. Moreover, TEAD1, MYO7A and NDUFC2 were identified as the candidate genes. Functionally, siRNA-mediated TEAD1, MYO7A and NDUFC2 knockdown significantly decreased the glucose uptake and inhibited the transcription of genes related to insulin and glucose metabolism pathways. CONCLUSIONS: Our study contributes novel insights to the understanding of hepatic glucose metabolism, demonstrating the impact of TEAD1, MYO7A and NDUFC2 on mitochondrial function in the liver and their regulatory role in maintaining in glucose homeostasis.

Systems genetics analysis reveals the common genetic basis for pain sensitivity and cognitive function.

BACKGROUND: There is growing evidence of a strong correlation between pain sensitivity and cognitive function under both physiological and pathological conditions. However, the detailed mechanisms remain largely unknown. In the current study, we sought to explore candidate genes and common molecular mechanisms underlying pain sensitivity and cognitive function with a transcriptome-wide association study using recombinant inbred mice from the BXD family. METHODS: The pain sensitivity determined by Hargreaves' paw withdrawal test and cognition-related phenotypes were systematically analyzed in 60 strains of BXD mice and correlated with hippocampus transcriptomes, followed by quantitative trait locus (QTL) mapping and systems genetics analysis. RESULTS: The pain sensitivity showed significant variability across the BXD strains and co-varies with cognitive traits. Pain sensitivity correlated hippocampual genes showed a significant involvement in cognition-related pathways, including glutamatergic synapse, and PI3K-Akt signaling pathway. Moreover, QTL mapping identified a genomic region on chromosome 4, potentially regulating the variation of pain sensitivity. Integrative analysis of expression QTL mapping, correlation analysis, and Bayesian network modeling identified Ring finger protein 20 (Rnf20) as the best candidate. Further pathway analysis indicated that Rnf20 may regulate the expression of pain sensitivity and cognitive function through the PI3K-Akt signaling pathway, particularly through interactions with genes Ppp2r2b, Ppp2r5c, Col9a3, Met, Rps6, Tnc, and Kras. CONCLUSIONS: Our study demonstrated that pain sensitivity is associated with genetic background and Rnf20-mediated PI3K-Akt signaling may involve in the regulation of pain sensitivity and cognitive functions.

Novel Pathogenic Mutation of P209L in TRPC6 Gene Causes Adult Focal Segmental Glomerulosclerosis.

Focal segmental glomerulosclerosis (FSGS) is a leading kidney disease, clinically associated with proteinuria and progressive renal failure. The occurrence of this disease is partly related to gene mutations. We describe a single affected family member who presented with FSGS. We used high-throughput sequencing, sanger sequencing to identify the pathogenic mutations, and a systems genetics analysis in the BXD mice was conducted to explore the genetic regulatory mechanisms of pathogenic genes in the development of FSGS. We identified high urinary protein (++++) and creatinine levels (149 µmol/L) in a 29-year-old male diagnosed with a 5-year history of grade 2 hypertension. Histopathology of the kidney biopsy showed stromal hyperplasia at the glomerular segmental sclerosis and endothelial cell vacuolation degeneration. Whole-exome sequencing followed by Sanger sequencing revealed a heterozygous missense mutation (c.643C > T) in exon 2 of TRPC6, leading to the substitution of arginine with tryptophan at position 215 (p.Arg215Trp). Systems genetics analysis of the 53 BXD mice kidney transcriptomes identified Pygm as the upstream regulator of Trpc6. Those two genes are jointly involved in the regulation of FSGS mainly via Wnt and Hippo signaling pathways. We present a novel variant in the TRPC6 gene that causes FSGS. Moreover, our data suggested TRPC6 works with PYGM, as well as Wnt and Hippo signaling pathways to regulate renal function, which could guide future clinical prevention and targeted treatment for FSGS outcomes.

When diet meets genetics.

Gene expression profiling of a diverse mouse population helps to decipher how a fat-rich diet contributes to inflammatory bowel disease.

Genetic and dietary modulators of the inflammatory response in the gastrointestinal tract of the BXD mouse genetic reference population.

Inflammatory gut disorders, including inflammatory bowel disease (IBD), can be impacted by dietary, environmental, and genetic factors. While the incidence of IBD is increasing worldwide, we still lack a complete understanding of the gene-by-environment interactions underlying inflammation and IBD. Here, we profiled the colon transcriptome of 52 BXD mouse strains fed with a chow or high-fat diet (HFD) and identified a subset of BXD strains that exhibit an IBD-like transcriptome signature on HFD, indicating that an interplay of genetics and diet can significantly affect intestinal inflammation. Using gene co-expression analyses, we identified modules that are enriched for IBD-dysregulated genes and found that these IBD-related modules share cis-regulatory elements that are responsive to the STAT2, SMAD3, and REL transcription factors. We used module quantitative trait locus analyses to identify genetic loci associated with the expression of these modules. Through a prioritization scheme involving systems genetics in the mouse and integration with external human datasets, we identified Muc4 and Epha6 as the top candidates mediating differences in HFD-driven intestinal inflammation. This work provides insights into the contribution of genetics and diet to IBD risk and identifies two candidate genes, MUC4 and EPHA6, that may mediate IBD susceptibility in humans.

Identification of the regulatory mechanism of ACE2 in COVID-19-induced kidney damage with systems genetics approach.

Studies showed that SARS-CoV-2 can directly target the kidney and induce renal damage. As the cell surface receptor for SARS-CoV-2 infection, the angiotensin-converting enzyme 2 (ACE2) plays a pivotal role for renal physiology and function. Thus, it is important to understand ACE2 through which pathway influences the pathogenesis of renal damage induced by COVID-19. In this study, we first performed an eQTL mapping for Ace2 in kidney tissues in 53 BXD mice strains. Results demonstrated that Ace2 is highly expressed and strongly controlled by a genetic locus on chromosome 16 in the kidney, with six genes (Dnase1, Vasn, Usp7, Abat, Mgrn1, and Rbfox1) dominated as the upstream modulator, as they are highly correlated with Ace2 expression. Gene co-expression analysis showed that Ace2 co-variates are significantly involved in the renin-angiotensin system (RAS) pathway which acts as a reno-protector. Importantly, we also found that Ace2 is positively correlated with Pdgf family members, particularly Pdgfc, which showed the most association among the 76 investigated growth factors. Mammalian Phenotype Ontology enrichment indicated that the cognate transcripts for both Ace2 and Pdgfc were mainly involved in regulating renal physiology and morphology. Among which, Cd44, Egfr, Met, Smad3, and Stat3 were identified as hub genes through protein-protein interaction analysis. Finally, in aligning with our systems genetics findings, we found ACE2, pdgf family members, and RAS genes decreased significantly in the CAKI-1 kidney cancer cells treated with S protein and receptor binding domain structural protein. Collectively, our data suggested that ACE2 work with RAS, PDGFC, as well as their cognate hub genes to regulate renal function, which could guide for future clinical prevention and targeted treatment for COVID-19-induced renal damage outcomes. KEY MESSAGES: • Ace2 is highly expressed and strongly controlled by a genetic locus on chromosome 16 in the kidney. • Ace2 co-variates are enriched in the RAS pathway. • Ace2 is strongly correlated with the growth factor Pdgfc. • Ace2 and Pdgfc co-expressed genes involved in the regulation of renal physiology and morphology. • SARS-CoV-2 spike glycoprotein induces down-regulation of Ace2, RAS, and Pdgfc.

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.

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 AVAILABILITY: The 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.

Transcriptome-wide association study reveals cholesterol metabolism gene Lpl is a key regulator of cognitive dysfunction.

Cholesterol metabolism in the brain plays a crucial role in normal physiological function, and its aberrations are associated with cognitive dysfunction. The present study aimed to determine which cholesterol-related genes play a vital role in cognitive dysfunction and to dissect its underlying molecular mechanisms using a systems genetics approach in the BXD mice family. We first systematically analyzed the association of expression of 280 hippocampal genes related to cholesterol metabolism with cognition-related traits and identified lipoprotein lipase (Lpl) as a critical regulator. This was further confirmed by phenome-wide association studies that indicate Lpl associated with hippocampus volume residuals and anxiety-related traits. By performing expression quantitative trait locus mapping, we demonstrate that Lpl is strongly cis-regulated in the BXD hippocampus. We also identified ∼3,300 genes significantly (p < 0.05) correlated with the Lpl expression. Those genes are mainly involved in the regulation of neuron-related traits through the MAPK signaling pathway, axon guidance, synaptic vesicle cycle, and NF-kappa B signaling pathway. Furthermore, a protein-protein interaction network analysis identified several direct interactors of Lpl, including Rab3a, Akt1, Igf1, Crp, and Lrp1, which indicates that Lpl involves in the regulation of cognitive dysfunction through Rab3a-mediated synaptic vesicle cycle and Akt1/Igf1/Crp/Lrp1-mediated MAPK signaling pathway. Our findings demonstrate the importance of the Lpl, among the cholesterol-related genes, in regulating cognitive dysfunction and highlighting the potential signaling pathways, which may serve as novel therapeutic targets for the treatment of cognitive dysfunction.

Genetic control of the pluripotency epigenome determines differentiation bias in mouse embryonic stem cells.

Genetically diverse pluripotent stem cells display varied, heritable responses to differentiation cues. Here, we harnessed these disparities through derivation of mouse embryonic stem cells from the BXD genetic reference panel, along with C57BL/6J (B6) and DBA/2J (D2) parental strains, to identify loci regulating cell state transitions. Upon transition to formative pluripotency, B6 stem cells quickly dissolved naïve networks adopting gene expression modules indicative of neuroectoderm lineages, whereas D2 retained aspects of naïve pluripotency. Spontaneous formation of embryoid bodies identified divergent differentiation where B6 showed a propensity toward neuroectoderm and D2 toward definitive endoderm. Genetic mapping identified major trans-acting loci co-regulating chromatin accessibility and gene expression in both naïve and formative pluripotency. These loci distally modulated occupancy of pluripotency factors at hundreds of regulatory elements. One trans-acting locus on Chr 12 primarily impacted chromatin accessibility in embryonic stem cells, while in epiblast-like cells, the same locus subsequently influenced expression of genes enriched for neurogenesis, suggesting early chromatin priming. These results demonstrate genetically determined biases in lineage commitment and identify major regulators of the pluripotency epigenome.

Hippocampal Transcriptome-Wide Association Study Reveals Correlations Between Impaired Glutamatergic Synapse Pathway and Age-Related Hearing Loss in BXD-Recombinant Inbred Mice.

Age-related hearing loss (ARHL) is associated with cognitive dysfunction; however, the detailed underlying mechanisms remain unclear. The aim of this study is to investigate the potential underlying mechanism with a system genetics approach. A transcriptome-wide association study was performed on aged (12-32 months old) BXD mice strains. The hippocampus gene expression was obtained from 56 BXD strains, and the hearing acuity was assessed from 54 BXD strains. Further correlation analysis identified a total of 1,435 hearing-related genes in the hippocampus (p < 0.05). Pathway analysis of these genes indicated that the impaired glutamatergic synapse pathway is involved in ARHL (p = 0.0038). Further gene co-expression analysis showed that the expression level of glutamine synthetase (Gls), which is significantly correlated with ARHL (n = 26, r = -0.46, p = 0.0193), is a crucial regulator in glutamatergic synapse pathway and associated with learning and memory behavior. In this study, we present the first systematic evaluation of hippocampus gene expression pattern associated with ARHL, learning, and memory behavior. Our results provide novel potential molecular mechanisms involved in ARHL and cognitive dysfunction association.

Ace2 and Tmprss2 Expressions Are Regulated by Dhx32 and Influence the Gastrointestinal Symptoms Caused by SARS-CoV-2.

Studies showed that the gastrointestinal (GI) tract is one of the most important pathways for SARS-CoV-2 infection and coronavirus disease 2019 (COVID-19). As SARS-CoV-2 cellular entry depends on the ACE2 receptor and TMPRSS2 priming of the spike protein, it is important to understand the molecular mechanisms through which these two proteins and their cognate transcripts interact and influence the pathogenesis of COVID-19. In this study, we quantified the expression, associations, genetic modulators, and molecular pathways for Tmprss2 and Ace2 mRNA expressions in GI tissues using a systems genetics approach and the expanded family of highly diverse BXD mouse strains. The results showed that both Tmprss2 and Ace2 are highly expressed in GI tissues with significant covariation. We identified a significant expression quantitative trait locus on chromosome 7 that controls the expression of both Tmprss2 and Ace2. Dhx32 was found to be the strongest candidate in this interval. Co-expression network analysis demonstrated that both Tmprss2 and Ace2 were located at the same module that is significantly associated with other GI-related traits. Protein-protein interaction analysis indicated that hub genes in this module are linked to circadian rhythms. Collectively, our data suggested that genes with circadian rhythms of expression may have an impact on COVID-19 disease, with implications related to the timing and treatment of COVID-19.

Lifespan-Associated Gene Expression Signatures of Recombinant BXD Mice Implicates Coro7 and Set in Longevity.

Although genetic approaches have identified key genes and pathways that promote longevity, systems-level approaches are less utilized. Here, we took advantage of the wealth of omics data characterizing the BXD family of mice. We associated transcript and peptide levels across five tissues from both female and male BXD isogenic lines with their median lifespan. We identified over 5000 genes that showed a longevity correlation in a given tissue. Surprisingly, we found less than 1% overlap among longevity-correlating genes across tissues and sex. These 1% shared genes consist of 51 genes, of which 13 have been shown to alter lifespan. Only two genes -Coro7 and Set- showed a longevity correlation in all tissues and in both sexes. While differential regulation of aging across tissues and sex has been reported, our systems-level analysis reveals two unique genes that may promote healthy aging in unique sex- and tissue-agnostic manner.

Paraquat Toxicogenetics: Strain-Related Reduction of Tyrosine Hydroxylase Staining in Substantia Nigra in Mice.

Paraquat (PQ) is a putative risk factor for the development of sporadic Parkinson's disease. To model a possible genetic basis for individual differences in susceptibility to exposure to PQ, we recently examined the effects of paraquat on tyrosine hydroxylase (TH)-containing neurons in the substantia nigra pars compacta (SNc) of six members of the BXD family of mice (n = 2-6 per strain). We injected males with 5 mg/kg paraquat weekly three times. The density of TH+ neurons counted by immunocytochemistry at 200x in eight or more sections through the SNc is reduced in five of the six strains relative to control (N = 4 ± 2 mice per strain). TH+ loss ranged from 0 to 20% with an SEM of 1%. The heritability was estimated using standard ANOVA and jackknife resampling and is 0.37 ± 0.05 in untreated animals and 0.47 ± 0.04 in treated animals. These results demonstrate genetic modulation and GxE variation in susceptibility to PQ exposure and the loss of TH staining in the substantia nigra.

The Genetic Dissection of Ace2 Expression Variation in the Heart of Murine Genetic Reference Population.

Background: A high inflammatory and cytokine burden that induces vascular inflammation, myocarditis, cardiac arrhythmias, and myocardial injury is associated with a lethal outcome in COVID-19. The SARS-CoV-2 virus utilizes the ACE2 receptor for cell entry in a similar way to SARS-CoV. This study investigates the regulation, gene network, and associated pathways of ACE2 that may be involved in inflammatory and cardiovascular complications of COVID-19. Methods: Cardiovascular traits were determined in the one of the largest mouse genetic reference populations: BXD recombinant inbred strains using blood pressure, electrocardiography, and echocardiography measurements. Expression quantitative trait locus (eQTL) mapping, genetic correlation, and functional enrichment analysis were used to identify Ace2 regulation, gene pathway, and co-expression networks. Results: A wide range of variation was found in expression of Ace2 among the BXD strains. Levels of Ace2 expression are negatively correlated with cardiovascular traits, including systolic and diastolic blood pressure and P wave duration and amplitude. Ace2 co-expressed genes are significantly involved in cardiac- and inflammatory-related pathways. The eQTL mapping revealed that Cyld is a candidate upstream regulator for Ace2. Moreover, the protein-protein interaction (PPI) network analysis inferred several potential key regulators (Cul3, Atf2, Vcp, Jun, Ppp1cc, Npm1, Mapk8, Set, Dlg1, Mapk14, and Hspa1b) for Ace2 co-expressed genes in the heart. Conclusions: Ace2 is associated with blood pressure, atrial morphology, and sinoatrial conduction in BXD mice. Ace2 co-varies with Atf2, Cyld, Jun, Mapk8, and Mapk14 and is enriched in the RAS, TGFß, TNFα, and p38α signaling pathways, involved in inflammation and cardiac damage. We suggest that all these novel Ace2-associated genes and pathways may be targeted for preventive, diagnostic, and therapeutic purposes in cardiovascular damage in patients with systemic inflammation, including COVID-19 patients.

Modeling the Genetic Basis of Individual Differences in Susceptibility to Gulf War Illness.

Between 25% and 30% of the nearly one million military personnel who participated in the 1991 Persian Gulf War became ill with chronic symptoms ranging from gastrointestinal to nervous system dysfunction. This disorder is now referred to as Gulf War Illness (GWI) and the underlying pathophysiology has been linked to exposure-based neuroinflammation caused by organophosphorous (OP) compounds coupled with high circulating glucocorticoids. In a mouse model of GWI we developed, corticosterone was shown to act synergistically with an OP (diisopropylflurophosphate) to dramatically increase proinflammatory cytokine gene expression in the brain. Because not all Gulf War participants became sick, the question arises as to whether differential genetic constitution might underlie individual differences in susceptibility. To address this question of genetic liability, we tested the impact of OP and glucocorticoid exposure in a genetic reference population of 30 inbred mouse strains. We also studied both sexes. The results showed wide differences among strains and overall that females were less sensitive to the combined treatment than males. Furthermore, we identified one OP-glucocorticoid locus and nominated a candidate gene-Spon1-that may underlie the marked differences in response.

Identification of a Functional Non-coding Variant in the GABA A Receptor α2 Subunit of the C57BL/6J Mouse Reference Genome: Major Implications for Neuroscience Research.

GABA type-A (GABA-A) receptors containing the α2 subunit (GABRA2) are expressed in most brain regions and are critical in modulating inhibitory synaptic function. Genetic variation at the GABRA2 locus has been implicated in epilepsy, affective and psychiatric disorders, alcoholism and drug abuse. Gabra2 expression varies as a function of genotype and is modulated by sequence variants in several brain structures and populations, including F2 crosses originating from C57BL/6J (B6J) and the BXD recombinant inbred family derived from B6J and DBA/2J. Here we demonstrate a global reduction of GABRA2 brain protein and mRNA in the B6J strain relative to other inbred strains, and identify and validate the causal mutation in B6J. The mutation is a single base pair deletion located in an intron adjacent to a splice acceptor site that only occurs in the B6J reference genome. The deletion became fixed in B6J between 1976 and 1991 and is now pervasive in many engineered lines, BXD strains generated after 1991, the Collaborative Cross, and the majority of consomic lines. Repair of the deletion using CRISPR-Cas9-mediated gene editing on a B6J genetic background completely restored brain levels of GABRA2 protein and mRNA. Comparison of transcript expression in hippocampus, cortex, and striatum between B6J and repaired genotypes revealed alterations in GABA-A receptor subunit expression, especially in striatum. These results suggest that naturally occurring variation in GABRA2 levels between B6J and other substrains or inbred strains may also explain strain differences in anxiety-like or alcohol and drug response traits related to striatal function. Characterization of the B6J private mutation in the Gabra2 gene is of critical importance to molecular genetic studies in neurobiological research because this strain is widely used to generate genetically engineered mice and murine genetic populations, and is the most widely utilized strain for evaluation of anxiety-like, depression-like, pain, epilepsy, and drug response traits that may be partly modulated by GABRA2 function.

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

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