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.

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.

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.

Echocardiography phenotyping in murine genetic reference population of BXD strains reveals significant QTLs associated with cardiac function and morphology.

The genetic reference population of recombinant inbred BXD mice has been derived from crosses between C57BL/6J and DBA/2J strains. The DBA/2J parent exhibits cardiomyopathy phenotypes, whereas C57BL/6J has normal heart. BXD mice are sequenced for studying genetic interactions in cardiomyopathies. The study aimed to assess cardiomyopathy traits in BXDs and investigate the quantitative genetic architecture of those traits. Echocardiography, blood pressure, and cardiomyocyte size parameters obtained from 44 strains of BXD family (n > 5/sex) at 4-5 mo of age were associated with heart transcriptomes and expression quantitative trait loci (eQTL) mapping was performed. More than twofold variance in ejection fraction (EF%), fractional shortening (FS%), left ventricular volumes (LVVols), internal dimensions (LVIDs), mass (LVM), and posterior wall (LVPW) thickness was found among BXDs. In male BXDs, eQTL mapping identified Ndrg4 on chromosome 8 QTL to be positively correlated with LVVol and LVID and negatively associated with cardiomyocyte diameter. In female BXDs, significant QTLs were found on chromosomes 7 and 3 to be associated with LVPW and EF% and FS%, respectively, and Josd2, Dap3, and Tpm3 were predicted as strong candidate genes. Our study found variable cardiovascular traits among BXD strains and identified multiple associated QTLs, suggesting an influence of genetic background on expression of echocardiographic and cardiomyocyte diameter traits. Increased LVVol and reduced EF% and FS% represented dilated cardiomyopathy, whereas increased LV mass and wall thickness indicated hypertrophic cardiomyopathy traits. The BXD family is ideal for identifying candidate genes, causal and modifier, that influence cardiovascular phenotypes.NEW & NOTEWORTHY This study aimed to establish a cardiac phenotype-genotype correlation in murine genetic reference population of BXD RI strains by phenotyping the echocardiography, blood pressure, and cardiomyocyte diameter traits and associating each collected phenotype with genetic background. Our study identified several QTLs and candidate genes that have significant association with cardiac hypertrophy, ventricular dilation, and function including systolic hyperfunction and dysfunction.

HOXD8 suppresses renal cell carcinoma growth by upregulating SHMT1 expression.

Amplification of amino acids synthesis is reported to promote tumorigenesis. The serine/glycine biosynthesis pathway is a reversible conversion of serine and glycine catalyzed by cytoplasmic serine hydroxymethyltransferase (SHMT)1 and mitochondrial SHMT2; however, the role of SHTM1 in renal cell carcinoma (RCC) is still unclear. We found that low SHMT1 expression is correlated with poor survival of RCC patients. The in vitro study showed that overexpression of SHMT1 suppressed RCC proliferation and migration. In the mouse tumor model, SHMT1 significantly retarded RCC tumor growth. Furthermore, by gene network analysis, we found several SHMT1-related genes, among which homeobox D8 (HOXD8) was identified as the SHMT1 regulator. Knockdown of HOXD8 decreased SHMT1 expression, resulting in faster RCC growth, and rescued the SHMT1 overexpression-induced cell migration defects. Additionally, ChIP assay found the binding site of HOXD8 to SHMT1 promoter was at the -456~-254 bp region. Taken together, SHMT1 functions as a tumor suppressor in RCC. The transcription factor HOXD8 can promote SHMT1 expression and suppress RCC cell proliferation and migration, which provides new mechanisms of SHMT1 in RCC tumor growth and might be used as a potential therapeutic target candidate for clinical treatment.

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

BACKGROUND: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: Compound splicing and T659I Tnni3k variants reduce cardiac Tnni3k expression and Tnni3k levels are associated with cardiac and glucose metabolism-related phenotypes.

Autoreactive B cells in SLE, villains or innocent bystanders?

The current concepts for development of autoreactive B cells in SLE (systemic lupus erythematosus) focus on extrinsic stimuli and factors that provoke B cells into tolerance loss. Traditionally, major tolerance loss pathways are thought to be regulated by factors outside the B cell including autoantigen engagement of the B-cell receptor (BCR) with simultaneous type I interferon (IFN) produced by dendritic cells, especially plasmacytoid dendritic cells (pDCs). Later, in autoreactive follicles, B-cells encounter T-follicular helper cells (Tfh) that produce interleukin (IL)-21, IL-4 and pathogenic cytokines, IL-17 and IFN gamma (IFNÉ£). This review discusses these mechanisms and also highlights recent advances pointing to the peripheral transitional B-cell stage as a major juncture where transient autocrine IFNß expression by developing B-cells imprints a heightened susceptibility to external factors favoring differentiation into autoantibody-producing plasmablasts. Recent studies highlight transitional B-cell heterogeneity as a determinant of intrinsic resistance or susceptibility to tolerance loss through the shaping of B-cell responsiveness to cytokines and other environment factors.

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.

Systems genetics analysis defines importance of TMEM43/LUMA for cardiac- and metabolic-related pathways.

Broad cellular functions and diseases including muscular dystrophy, arrhythmogenic right ventricular cardiomyopathy (ARVC5) and cancer are associated with transmembrane protein43 (TMEM43/LUMA). The study aimed to investigate biological roles of TMEM43 through genetic regulation, gene pathways and gene networks, candidate interacting genes, and up- or downstream regulators. Cardiac transcriptomes from 40 strains of recombinant inbred BXD mice and two parental strains representing murine genetic reference population (GRP) were applied for genetic correlation, functional enrichment, and coexpression network analysis using systems genetics approach. The results were validated in a newly created knock-in Tmem43-S358L mutation mouse model (Tmem43S358L) that displayed signs of cardiac dysfunction, resembling ARVC5 phenotype seen in humans. We found high Tmem43 levels among BXDs with broad variability in expression. Expression of Tmem43 highly negatively correlated with heart mass and heart rate among BXDs, whereas levels of Tmem43 highly positively correlated with plasma high-density lipoproteins (HDL). Through finding differentially expressed genes (DEGs) between Tmem43S358L mutant and wild-type (Tmem43WT) lines, 18 pathways (out of 42 found in BXDs GRP) that are involved in ARVC, hypertrophic cardiomyopathy, dilated cardiomyopathy, nonalcoholic fatty liver disease, Alzheimer's disease, Parkinson's disease, and Huntington's disease were verified. We further constructed Tmem43-mediated gene network, in which Ctnna1, Adcy6, Gnas, Ndufs6, and Uqcrc2 were significantly altered in Tmem43S358L mice versus Tmem43WT controls. Our study defined the importance of Tmem43 for cardiac- and metabolism-related pathways, suggesting that cardiovascular disease-relevant risk factors may also increase risk of metabolic and neurodegenerative diseases via TMEM43-mediated pathways.

Identification of cyclin D1 as a major modulator of 3-nitropropionic acid-induced striatal neurodegeneration.

Mitochondria dysfunction occurs in the aging brain as well as in several neurodegenerative disorders and predisposes neuronal cells to enhanced sensitivity to neurotoxins. 3-nitropropionic acid (3-NP) is a naturally occurring plant and fungal neurotoxin that causes neurodegeneration predominantly in the striatum by irreversibly inhibiting the tricarboxylic acid respiratory chain enzyme, succinate dehydrogenase (SDH), the main constituent of the mitochondria respiratory chain complex II. Significantly, although 3-NP-induced inhibition of SDH occurs in all brain regions, neurodegeneration occurs primarily and almost exclusively in the striatum for reasons still not understood. In rodents, 3-NP-induced striatal neurodegeneration depends on the strain background suggesting that genetic differences among genotypes modulate toxicant variability and mechanisms that underlie 3-NP-induced neuronal cell death. Using the large BXD family of recombinant inbred (RI) strains we demonstrate that variants in Ccnd1 - the gene encoding cyclin D1 - of the DBA/2 J parent underlie the resistance to 3-NP-induced striatal neurodegeneration. In contrast, the Ccnd1 variant inherited from the widely used C57BL/6 J parental strain confers sensitivity. Given that cellular stress triggers induction of cyclin D1 expression followed by cell-cycle re-entry and consequent neuronal cell death, we sought to determine if the C57BL/6 J and DBA/2 J Ccnd1 variants are differentially modulated in response to 3-NP. We confirm that 3-NP induces cyclin D1 expression in striatal neuronal cells of C57BL/6 J, but this response is blunted in the DBA/2 J. We further show that striatal-specific alternative processing of a highly conserved 3'UTR negative regulatory region of Ccnd1 co-segregates with the C57BL/6 J parental Ccnd1 allele in BXD strains and that its differential processing accounts for sensitivity or resistance to 3-NP. Our results indicate that naturally occurring Ccnd1 variants may play a role in the variability observed in neurodegenerative disorders involving mitochondria complex II dysfunction and point to cyclin D1 as a possible therapeutic target.

New insights into the transcriptional regulation of benzoxazinoid biosynthesis in wheat.

This article comments on: Batyrshina ZS, Shavit R, Yaakov B, Bocobza S, Tzin V. 2022. The transcription factor gene TaMYB31 regulates the benzoxazinoid biosynthetic pathway in wheat. Journal of Experimental Botany73, 5634-5649. Benzoxazinoids (BXDs) are abundant indole-derived specialized metabolites in several monocot crop species, such as wheat, maize, and rye. They function in plant immunity against herbivorous arthropods and fungal pathogens, but also as iron chelators, in metal tolerance, and as allelochemicals. Although BXD biosynthetic pathways have been studied extensively and are well described, information about the transcriptional regulation of BXD biosynthesis is scarce. In the current issue of JXB, Batyrshina et al. (2022) identified the transcription factor gene TaMYB31 in the tetraploid wheat Triticum turgidum and verified its function as a component of BXD metabolism in the hexaploid wheat Triticum aestivum, where it regulates constitutive and stress-inducible BXD biosynthesis.

Boundaries potentiate polycomb response element-mediated silencing.

BACKGROUND: Epigenetic memory plays a critical role in the establishment and maintenance of cell identities in multicellular organisms. Polycomb and trithorax group (PcG and TrxG) proteins are responsible for epigenetic memory, and in flies, they are recruited to specialized DNA regulatory elements termed polycomb response elements (PREs). Previous transgene studies have shown that PREs can silence reporter genes outside of their normal context, often by pairing sensitive (PSS) mechanism; however, their silencing activity is non-autonomous and depends upon the surrounding chromatin context. It is not known why PRE activity depends on the local environment or what outside factors can induce silencing. RESULTS: Using an attP system in Drosophila, we find that the so-called neutral chromatin environments vary substantially in their ability to support the silencing activity of the well-characterized bxdPRE. In refractory chromosomal contexts, factors required for PcG-silencing are unable to gain access to the PRE. Silencing activity can be rescued by linking the bxdPRE to a boundary element (insulator). When placed next to the PRE, the boundaries induce an alteration in chromatin structure enabling factors critical for PcG silencing to gain access to the bxdPRE. When placed at a distance from the bxdPRE, boundaries induce PSS by bringing the bxdPREs on each homolog in close proximity. CONCLUSION: This proof-of-concept study demonstrates that the repressing activity of PREs can be induced or enhanced by nearby boundary elements.

Genetic differences in ethanol consumption: effects on iron, copper, and zinc regulation in mouse hippocampus.

One common characteristic of neurodegenerative diseases is dysregulation of iron, usually with observed increases in its concentration in various regions. Heavy alcohol consumption is believed to contribute to such iron dysregulation in the brain with accompanying dementia. To examine this effect and related genetic-based individual differences in an animal model, we subjected female mice from 12 BXD recombinant inbred strains to 16 weeks of alcohol consumption using the drinking in the dark (DID) method. Daily consumption was recorded and at the end of 16 weeks hippocampus tissues harvested. Concentrations of iron, copper and zinc were measured using X-ray fluorescence technology. The results showed that, DID increased iron overall across all strains, ranging from 3 to 68%. Copper and Zinc both decreased, ranging from 0.4-42 and 5-35% respectively. Analysis of variance revealed significant strain by treatment interactions for all three metals. Additionally, in the DID group, we observed strain differences in reduction of hippocampus mass. These findings are particularly interesting to us because high alcohol consumption in humans has been associated with neurodegeneration and dementia related to disruption of iron regulation. The findings of alcohol consumption associated decreases in copper and zinc are novel. The role of copper regulation and neurological function related to alcohol consumption is as yet largely unexplored. The role of zinc is better known as a neuromodulator in the hippocampus and appears to be protective against neurological damage. It would seem then, that the alcohol-related decrease in zinc in the hippocampus would be of concern and warrants further study.

Genome-wide transcriptome architecture in a mouse model of Gulf War Illness.

Gulf War Illness (GWI) is thought to be a chronic neuroimmune disorder caused by in-theater exposure during the 1990-1991 Gulf War. There is a consensus that the illness is caused by exposure to insecticides and nerve agent toxicants. However, the heterogeneity in both development of disease and clinical outcomes strongly suggests a genetic contribution. Here, we modeled GWI in 30 BXD recombinant inbred mouse strains with a combined treatment of corticosterone (CORT) and diisopropyl fluorophosphate (DFP). We quantified transcriptomes from 409 prefrontal cortex samples. Compared to the untreated and DFP treated controls, the combined treatment significantly activated pathways such as cytokine-cytokine receptor interaction and TNF signaling pathway. Protein-protein interaction analysis defined 6 subnetworks for CORT + DFP, with the key regulators being Cxcl1, Il6, Ccnb1, Tnf, Agt, and Itgam. We also identified 21 differentially expressed genes having significant QTLs related to CORT + DFP, but without evidence for untreated and DFP treated controls, suggesting regions of the genome specifically involved in the response to CORT + DFP. We identified Adamts9 as a potential contributor to response to CORT + DFP and found links to symptoms of GWI. Furthermore, we observed a significant effect of CORT + DFP treatment on the relative proportion of myelinating oligodendrocytes, with a QTL on Chromosome 5. We highlight three candidates, Magi2, Sema3c, and Gnai1, based on their high expression in the brain and oligodendrocyte. In summary, our results show significant genetic effects of the CORT + DFP treatment, which mirrors gene and protein expression changes seen in GWI sufferers, providing insight into the disease and a testbed for future interventions.

Systems Genetics Approaches in Mouse Models of Group A Streptococcal Necrotizing Soft-Tissue Infections.

Mouse models are invaluable resources for studying the pathogenesis and preclinical evaluation of therapeutics and vaccines against many human pathogens. Infections caused by group A streptococcus (GAS, Streptococcus pyogenes) are heterogeneous ranging from mild pharyngitis to severe invasive necrotizing fasciitis, a subgroup of necrotizing soft-tissue infections (NSTIs). While several strains of mice including BALB/c, C3H/HeN, CBA/J, and C57BL/10 offered significant insights, the human specificity and the interindividual variations on susceptibility or resistance to GAS infections limit their ability to mirror responses as seen in humans. In this chapter, we discuss the advanced recombinant inbred (ARI) BXD mouse model that mimics the genetic diversity as seen in humans and underpins the feasibility to map multiple genes (genetic loci) modulating GAS NSTI. GAS produces a myriad of virulence factors, including superantigens (SAg). Superantigens are potent immune toxins that activate T cells by cross-linking T cell receptors with human leukocyte antigen class-II (HLA-II) molecules expressed on antigen-presenting cells. This leads to a pro-inflammatory cytokine storm and the subsequent multiple organ damage and shock. Inbred mice are innately refractive to SAg-mediated responses. In this chapter, we discuss the versatility of the HLA-II transgenic mouse model that allowed the biological validation of known genetic associations to GAS NSTI. The combined utility of ARI-BXD and HLA-II mice as complementary approaches that offer clinically translatable insights into pathomechanisms driven by complex traits and host genetic context and novel means to evaluate the in vivo efficiency of therapies to improve outcomes of GAS NSTI are also discussed.

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.

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

BACKGROUND: 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. METHODS: 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.8 g/kg EtOH or isocaloric maltose-dextrin in 2 doses via intragastric gavage. Embryos were collected 7 hours 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. RESULTS: 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). CONCLUSIONS: 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.

A novel heat shock protein alpha 8 (Hspa8) molecular network mediating responses to stress- and ethanol-related behaviors.

Genetic differences mediate individual differences in susceptibility and responses to stress and ethanol, although, the specific molecular pathways that control these responses are not fully understood. Heat shock protein alpha 8 (Hspa8) is a molecular chaperone and member of the heat shock protein family that plays an integral role in the stress response and that has been implicated as an ethanol-responsive gene. Therefore, we assessed its role in mediating responses to stress and ethanol across varying genetic backgrounds. The hippocampus is an important mediator of these responses, and thus, was examined in the BXD family of mice in this study. We conducted bioinformatic analyses to dissect genetic factors modulating Hspa8 expression, identify downstream targets of Hspa8, and examined its role. Hspa8 is trans-regulated by a gene or genes on chromosome 14 and is part of a molecular network that regulates stress- and ethanol-related behaviors. To determine additional components of this network, we identified direct or indirect targets of Hspa8 and show that these genes, as predicted, participate in processes such as protein folding and organic substance metabolic processes. Two phenotypes that map to the Hspa8 locus are anxiety-related and numerous other anxiety- and/or ethanol-related behaviors significantly correlate with Hspa8 expression. To more directly assay this relationship, we examined differences in gene expression following exposure to stress or alcohol and showed treatment-related differential expression of Hspa8 and a subset of the members of its network. Our findings suggest that Hspa8 plays a vital role in genetic differences in responses to stress and ethanol and their interactions.

Identification of candidate genes that underlie the QTL on chromosome 1 that mediates genetic differences in stress-ethanol interactions.

Alcoholism, stress, and anxiety are strongly interacting heritable, polygenetic traits. In a previous study, we identified a quantitative trait locus (QTL) on murine chromosome (Chr) 1 between 23.0 and 31.5 Mb that modulates genetic differences in the effects of ethanol on anxiety-related phenotypes. The goal of the present study was to extend the analysis of this locus with a focus on identifying candidate genes using newly available data and tools. Anxiety-like behavior was evaluated with an elevated zero maze following saline or ethanol injections (1.8 g/kg) in C57BL/6J, DBA2J, and 72 BXD strains. We detected significant effects of strain and treatment and their interaction on anxiety-related behaviors, although surprisingly, sex was not a significant factor. The Chr1 QTL is specific to the ethanol-treated cohort. Candidate genes in this locus were evaluated using now standard bioinformatic criteria. Collagen 19a1 (Col19a1) and family sequence 135a (Fam135a) met most criteria but have lower expression levels and lacked biological verification and, therefore, were considered less likely candidates. In contrast, two other genes, the prenylated protein tyrosine phosphate family member Ptp4a1 (protein tyrosine phosphate 4a1) and the zinc finger protein Phf3 (plant homeoDomain finger protein 3) met each of our bioinformatic criteria and are thus strong candidates. These findings are also of translational relevance because both Ptp4a1 and Phf3 have been nominated as candidates genes for alcohol dependence in a human genome-wide association study. Our findings support the hypothesis that variants in one or both of these genes modulate heritable differences in the effects of ethanol on anxiety-related behaviors.

Identification of ethanol analgesia quantitative trait loci and candidate genes in BXD recombinant inbred mouse lines.

Alcohol consumption produces acute analgesic effects, and people experiencing pain conditions may drink alcohol to alleviate discomfort. However, tolerance to the analgesic properties of alcohol could prompt escalating consumption and dependence. Both nociception and alcohol-induced analgesia are under significant genetic control. Understanding the genetic architecture of these processes could inform better treatment options for people with pain conditions. This study aims to identify quantitative trait loci (QTL) driving variation in ethanol-induced analgesia across BXD recombinant inbred mouse lines. Male and female mice from 62 BXD strains received ethanol or saline oral gavage for five days and were tested for hot plate (HP) latency at baseline, Day 1, and Day 5. QTL mapping of HP phenotypes identified a significant provisional QTL on chromosome 17 for Day 1 HP latency in mice receiving ethanol. An additional highly suggestive QTL was present on chromosome 9 for the difference in pre- and post-ethanol thermal nociception. Candidate genes within QTL support intervals were provisionally identified using HP phenotypic correlations to transcriptomic database, expression QTL analysis, and other bioinformatics inquiries. The combined behavioral and bioinformatic analyses yielded strong ethanol analgesia candidate genes, specifically Myo6. Thus, the results of this genetic study of ethanol-induced analgesia in BXD mouse strains may contribute significantly to our understanding of the molecular basis for individual variation in the analgesic response to acute ethanol.