We have been using the Inbred Long- and Short-Sleep mouse strains (ILS, ISS) and a recombinant inbred panel derived from them, the LXS, to investigate the genetic underpinnings of acute ethanol tolerance which is considered to be a risk factor for alcohol use disorders (AUDs). Here, we have used RNA-seq to examine the transcriptome of whole brain in 40 of the LXS strains 8 hours after a saline or ethanol "pretreatment" as in previous behavioral studies. Approximately 1/3 of the 14,184 expressed genes were significantly heritable and many were unique to the pretreatment. Several thousand cis- and trans-eQTLs were mapped; a portion of these also were unique to pretreatment. Ethanol pretreatment caused differential expression (DE) of 1,230 genes. Gene Ontology (GO) enrichment analysis suggested involvement in numerous biological processes including astrocyte differentiation, histone acetylation, mRNA splicing, and neuron projection development. Genetic correlation analysis identified hundreds of genes that were correlated to the behaviors. GO analysis indicated that these genes are involved in gene expression, chromosome organization, and protein transport, among others. The expression profiles of the DE genes and genes correlated to AFT in the ethanol pretreatment group (AFT-Et) were found to be similar to profiles of HDAC inhibitors. Hdac1, a cis-regulated gene that is located at the peak of a previously mapped QTL for AFT-Et, was correlated to 437 genes, most of which were also correlated to AFT-Et. GO analysis of these genes identified several enriched biological process terms including neuron-neuron synaptic transmission and potassium transport. In summary, the results suggest widespread genetic effects on gene expression, including effects that are pretreatment-specific. A number of candidate genes and biological functions were identified that could be mediating the behavioral responses. The most prominent of these was Hdac1 which may be regulating genes associated with glutamatergic signaling and potassium conductance.
Genetic reference panels are widely used to map complex, quantitative traits in model organisms. We have generated new high-resolution genetic maps of 259 mouse inbred strains from recombinant inbred strain panels (C57BL/6J × DBA/2J, ILS/IbgTejJ × ISS/IbgTejJ, and C57BL/6J × A/J) and chromosome substitution strain panels (C57BL/6J-Chr#, C57BL/6J-Chr#, and C57BL/6J-Chr#). We genotyped all samples using the Affymetrix Mouse Diversity Array with an average intermarker spacing of 4.3 kb. The new genetic maps provide increased precision in the localization of recombination breakpoints compared to the previous maps. Although the strains were presumed to be fully inbred, we found residual heterozygosity in 40% of individual mice from five of the six panels. We also identified de novo deletions and duplications, in homozygous or heterozygous state, ranging in size from 21 kb to 8.4 Mb. Almost two-thirds (46 out of 76) of these deletions overlap exons of protein coding genes and may have phenotypic consequences. Twenty-nine putative gene conversions were identified in the chromosome substitution strains. We find that gene conversions are more likely to occur in regions where the homologous chromosomes are more similar. The raw genotyping data and genetic maps of these strain panels are available at http://churchill-lab.jax.org/website/MDA.
The Inbred Long- and Short-Sleep (ILS, ISS) mouse lines were selected for differences in acute ethanol sensitivity using the loss of righting response (LORR) as the selection trait. The lines show an over tenfold difference in LORR and, along with a recombinant inbred panel derived from them (the LXS), have been widely used to dissect the genetic underpinnings of acute ethanol sensitivity. Here we have sequenced the genomes of the ILS and ISS to investigate the DNA variants that contribute to their sensitivity difference. We identified ~2.7 million high-confidence SNPs and small indels and ~7000 structural variants between the lines; variants were found to occur in 6382 annotated genes. Using a hidden Markov model, we were able to reconstruct the genome-wide ancestry patterns of the eight inbred progenitor strains from which the ILS and ISS were derived, and found that quantitative trait loci that have been mapped for LORR were slightly enriched for DNA variants. Finally, by mapping and quantifying RNA-seq reads from the ILS and ISS to their strain-specific genomes rather than to the reference genome, we found a substantial improvement in a differential expression analysis between the lines. This work will help in identifying and characterizing the DNA sequence variants that contribute to the difference in ethanol sensitivity between the ILS and ISS and will also aid in accurate quantification of RNA-seq data generated from the LXS RIs.
The ILSXISS (LXS) recombinant inbred (RI) panel of mice is a valuable resource for genetic mapping studies of complex traits, due to its genetic diversity and large number of strains. Male and female mice from this panel were used to investigate genetic influences on alcohol consumption in the "drinking in the dark" (DID) model. Male mice (38 strains) and female mice (36 strains) were given access to 20% ethanol during the early phase of their circadian dark cycle for four consecutive days. The first principal component of alcohol consumption measures on days 2, 3, and 4 was used as a phenotype (DID phenotype) to calculate QTLs, using a SNP marker set for the LXS RI panel. Five QTLs were identified, three of which included a significant genotype by sex interaction, i.e., a significant genotype effect in males and not females. To investigate candidate genes associated with the DID phenotype, data from brain microarray analysis (Affymetrix Mouse Exon 1.0 ST Arrays) of male LXS RI strains were combined with RNA-Seq data (mouse brain transcriptome reconstruction) from the parental ILS and ISS strains in order to identify expressed mouse brain transcripts. Candidate genes were determined based on common eQTL and DID phenotype QTL regions and correlation of transcript expression levels with the DID phenotype. The resulting candidate genes (in particular, Arntl/Bmal1) focused attention on the influence of circadian regulation on the variation in the DID phenotype in this population of mice.
BACKGROUND:We previously reported that acute functional tolerance (AFT) to the hypnotic effects of alcohol was significantly correlated with drinking in the dark (DID) in the LXS recombinant inbred panel, but only in mice that had been pretreated with alcohol. Here, we have conducted quantitative trait locus (QTL) mapping for AFT. DNA sequencing of the progenitor ILS and ISS strains and microarray analyses were also conducted to identify candidate genes and functional correlates.METHODS:LXS mice were given either saline or alcohol (5 g/kg) on day 1 and then tested for loss of righting reflex AFT on day 2. QTLs were mapped using standard procedures. Two microarray analyses from brain were conducted: (i) naïve LXS mice and (ii) an alcohol treatment time course in the ILS and ISS. The full genomes of the ILS and ISS were sequenced to a depth of approximately 30×.RESULTS:A significant QTL for AFT in the alcohol pretreatment group was mapped to distal chromosome 4; numerous suggestive QTLs were also mapped. Preference drinking and DID have previously been mapped to the chromosome 4 locus. The credible interval of the significant chromosome 4 QTL spanned 23 Mb and included 716 annotated genes of which 150 had at least 1 nonsynonymous single nucleotide polymorphism or small indel that differed between the ILS and ISS; expression of 48 of the genes was cis-regulated. Enrichment analysis indicated broad functional categories underlying AFT, including proteolysis, transcription regulation, chromatin modification, protein kinase activity, and apoptosis.CONCLUSIONS:The chromosome 4 QTL is a key region containing possibly pleiotropic genes for AFT and drinking behavior. Given that the region contains many viable candidates and a large number of the genes in the interval fall into 1 or more of the enriched functional categories, we postulate that many genes of varying effect size contribute to the observed QTL effect.
Background Many quantitative trait loci (QTL) influencing mouse model phenotypes for alcoholism have been mapped genetically. However, the gene(s) comprising the QTL (QTG) are largely unknown. In previous work, Bennett and colleagues created congenic strains carrying the DBA/2IBG (D2) region for alcohol preference (AP) on chromosome 2, on a C57BL/6IBG (B6) background [1]. Subsequently, interval specific congenic recombinant strains (ISCRS), in which the full D2 QTL region was broken into smaller, partially overlapping regions of introgression, were generated and tested. With information from two ISCRS, the QTL has been mapped onto mouse chromosome 2 (Chr2) in a region of 3.4Mb by using C57BL/6J (B6) x DBA/2J (D2) recombinant inbred (RI) strains as well as by using F2 populations. Several candidate genes, Gad1, Atp5g3, Atf2, Sp3 and Sp9, have been evaluated but none of them has been confirmed for a definitive role in the regulation of the QTL of AP on Chr2 [2,3].
It has been well documented that genetic factors can influence predisposition to develop alcoholism. While the underlying genomic changes may be of several types, two of the most common and disease associated are copy number variations (CNVs) and sequence alterations of protein coding regions. The goal of this study was to identify CNVs and single-nucleotide polymorphisms that occur in gene coding regions that may play a role in influencing the risk of an individual developing alcoholism. Toward this end, two mouse strains were used that have been selectively bred based on their differential sensitivity to alcohol: the Inbred long sleep (ILS) and Inbred short sleep (ISS) mouse strains. Differences in initial response to alcohol have been linked to risk for alcoholism, and the ILS/ISS strains are used to investigate the genetics of initial sensitivity to alcohol. Array comparative genomic hybridization (arrayCGH) and exome sequencing were conducted to identify CNVs and gene coding sequence differences, respectively, between ILS and ISS mice. Mouse arrayCGH was performed using catalog Agilent 1 × 244 k mouse arrays. Subsequently, exome sequencing was carried out using an Illumina HiSeq 2000 instrument. ArrayCGH detected 74 CNVs that were strain-specific (38 ILS/36 ISS), including several ISS-specific deletions that contained genes implicated in brain function and neurotransmitter release. Among several interesting coding variations detected by exome sequencing was the gain of a premature stop codon in the alpha-amylase 2B (AMY2B) gene specifically in the ILS strain. In total, exome sequencing detected 2,597 and 1,768 strain-specific exonic gene variants in the ILS and ISS mice, respectively. This study represents the most comprehensive and detailed genomic comparison of ILS and ISS mouse strains to date. The two complementary genome-wide approaches identified strain-specific CNVs and gene coding sequence variations that should provide strong candidates to contribute to the alcohol-related phenotypic differences associated with these strains.
The purpose of this study is to investigate whether expression profiles of alcoholism-relevant genes in different parts of the brain are correlated differently with those in the liver. Four experiments were conducted. First, we used gene expression profiles from five parts of the brain (striatum, prefrontal cortex, nucleus accumbens, hippocampus, and cerebellum) and from liver in a population of recombinant inbred mouse strains to examine the expression association of 10 alcoholism-relevant genes. Second, we conducted the same association analysis between brain structures and the lung. Third, using five randomly selected, nonalcoholism-relevant genes, we conducted the association analysis between brain and liver. Finally, we compared the expression of 10 alcoholism-relevant genes in hippocampus and cerebellum between an alcohol preference strain and a wild-type control. We observed a difference in correlation patterns in expression levels of 10 alcoholism-relevant genes between different parts of the brain with those of liver. We then examined the association of gene expression between alcohol dehydrogenases ( Adh1 , Adh2 , Adh5 , and Adh7 ) and different parts of the brain. The results were similar to those of the 10 genes. Then, we found that the association of those genes between brain structures and lung was different from that of liver. Next, we found that the association patterns of five alcoholism-nonrelevant genes were different from those of 10 alcoholism-relevant genes. Finally, we found that the expression level of 10 alcohol-relevant genes is influenced more in hippocampus than in cerebellum in the alcohol preference strain. Our results show that the expression of alcoholism-relevant genes in liver is differently associated with the expression of genes in different parts of the brain. Because different structural changes in different parts of the brain in alcoholism have been reported, it is important to investigate whether those structural differences in the brains of those with alcoholism are due to the difference in the associations of gene expression between genes in liver and in different parts of the brain.
BACKGROUND:We hypothesized that rapid tolerance (1-day tolerance) for the duration of the loss of righting reflex ("sleep time" [ST]) was mediated by an increase in acute functional tolerance (AFT). We also hypothesized that increased AFT would correspond to increased drinking. These questions were addressed using the LXS recombinant inbred mouse strain panel.METHODS:Mice were given a pretreatment dose of either saline or 5 g/kg alcohol on day 1. On day 2, mice were tested for ST (4.1 g/kg) using a method with which it is possible to accurately assess AFT. Genetic correlation analysis was conducted among the ST-related variables and also with "drinking in the dark" (DID) which was previously measured by Saba and colleagues (2011).RESULTS:Saline-pretreated mice showed a continuous distribution of ST ranging from ~40 minutes to over 3 hours. Of the 43 strains tested, 9 showed significantly decreased ST after alcohol pretreatment, while in 3 strains, ST was significantly increased. AFT scores ranged from 0 to over 200 mg% in the saline group, and in the alcohol group, 8 strains showed a significant increase in AFT and 2 strains showed significant decrease in AFT. In the saline group, AFT was significantly correlated with ST (r = -0.47), but not in the alcohol group (r = -0.22). DID was significantly correlated with only AFT in the alcohol pretreated group (r = 0.64).CONCLUSIONS:The results suggest that AFT is an important component of the overall ST response, but that the alcohol pretreatment-induced change in AFT does not contribute to rapid ST tolerance. The significant correlation between DID and AFT in the alcohol group suggests that AFT may be a more relevant predictor of drinking behavior than the static measurement of ST. Moreover, preexposure to alcohol seems to change AFT in a way that makes it an even stronger predictor of drinking behavior.
In prior work, congenic strains carrying the DBA/2Igb (D2) region of chromosome 2 (Chr2) for alcohol preference were bred onto a C57BL/6Ibg (B6) background and as predicted were found to reduce voluntary consumption.Subsequently, interval-specific congenic recombinant strains (ISCRS) were generated and also tested.These ISCRS strains reduced the quantitative trait loci (QTL) interval to a comparatively small 3.4 Mb region.Here, we have exploited an integrative approach using both murine and human populations to critically evaluate candidate genes within this region.First, we used bioinformatics tools to search for genes relevant to alcohol preference within the QTL region.Second, we searched for single nucleotide polymorphisms (SNPs) within exons of every gene in this region.Third, we conducted follow-up microarray analyses to identify differentially expressed genes between the B6 and ISCRS strains in mice from each group.Fourth, we analyzed correlations between the expression level of candidate genes and phenotypes of alcohol preference in a large family of BXD recombinant inbred strains derived from B6 and D2.Finally, we evaluated SNP segregation in both BXD mouse strains and in humans who were heavy alcohol drinkers or non-drinkers.Among several potential candidate genes in this region, we identified activating transcription factor 2 (Atf2) as the most plausible gene that would influence alcohol preference.However, the candidacy of Atf2 was only weakly supported when we used a genetic network approach and by focused reanalysis of genome-wide association study data from European-American and African-American populations.Thus, we cannot conclude that Atf2 plays a role in the regulation of the QTL of mouse Chr2.
The mitochondrial ATP synthase, subunit c, isoform 3 gene (Atp5g3) encodes subunit 9, the subunit of the multisubunit enzyme that catalyzes ATP synthesis during oxidative phosphorylation in mitochondria. According to the Ensembl database, Atp5g3 in mice is located on chromosome 2 between 73746504 and 73749383 bp, within the genomic regions of two sets of quantitative trait loci - alcohol preference and body weight. Both of those traits are more influenced by epigenetic factors than many other traits are. Using currently available phenotype and gene expression profiles from the GeneNetwork database, we obtained correlations between Atp5g3 and alcoholism- and obesity-relevant phenotypes. The correlation in expression levels between Atp5g3 and each of its 12 partner genes in the molecular interaction are different in various tissues and genes. Transcriptome mapping indicated that Atp5g3 is differentially regulated in the hippocampus, cerebellum, and liver. Owing to a lack of known polymorphisms of Atp5g3 among three relevant mouse strains, C57BL/6J (B6), DBA/2J (D2), and BALB/cJ, the molecular mechanism for the connection between Atp5g3 and alcoholism and body weight requires further investigation.
Genetic influences on the predisposition to complex behavioral or physiological traits can reflect genetic polymorphisms that lead to altered gene product function, and/or variations in gene expression levels. We have explored quantitative variations in an animal's alcohol consumption, using a genetical genomic/phenomic approach. In our studies, gene expression is correlated with amount of alcohol consumed, and genomic regions that regulate the alcohol consumption behavior and the quantitative levels of gene expression (behavioral and expression quantitative trait loci [QTL]) are determined and used as a filter to identify candidate genes predisposing the behavior. We determined QTLs for alcohol consumption using the LXS panel of recombinant inbred mice. We then identified genes that were: 1) differentially expressed between five high and five low alcohol-consuming lines or strains of mice; and 2) were physically located in, or had an expression QTL (eQTL) within the alcohol consumption QTLs. Comparison of mRNA and protein levels in brains of high and low alcohol consuming mice led us to a bioinformatic examination of potential regulation by microRNAs of an identified candidate transcript, Gnb1 (G protein beta subunit 1). We combined our current analysis with our earlier work identifying candidate genes for the alcohol consumption trait in mice, rats and humans. Our overall analysis leads us to postulate that the activity of the GABAergic system, and in particular GABA release and GABA receptor trafficking and signaling, which involves G protein function, contributes significantly to genetic variation in the predisposition to varying levels of alcohol consumption. This article is part of a Special Issue entitled 'Trends in Neuropharmacology: In Memory of Erminio Costa'.
Results By using the keyword alcohol, we identified a total of 39 genetic elements in the region between markers D2Mit56 and D2Mit10. Among these genetic elements, we found seven with potential function in alcohol preference (Table 1). We then examined the SNPs, insertions and deletions, and gene expression levels of those seven genes. Conclusions Our current data suggest that the Atf2 and Titin genes are potentially the most alcohol relevant genes. However, further experiments and examination are still needed to confirm their candidacy. Several other candidate genes are also in the process of being identified.
Our laboratory has developed an online interactive resource called PhenoGen ( http://phenogen.ucdenver.edu ) which provides an archive of brain and other organ gene expression data from a panel of 20 common inbred mouse strains, and three recombinant inbred (RI) panels (two mouse and one rat). DNA microarray data can also be uploaded to the site where numerous analytical tools can be implemented. An important advantage to the archived data is that each array represents data from a single animal and each strain was sampled 4–7 times, providing an estimate of genetic variance (heritability) of individual transcript levels. These panels also allow genetic mapping of expression QTLs. Overlap of eQTLs with phenotypic QTLs provides a powerful approach to candidate gene identification. These methods are briefly described here and we encourage the use of our site for both scientific discovery and as a teaching tool in quantitative genetics.
Progressive sensorineural hearing loss is the most common form of acquired hearing impairment in the human population. It is also highly prevalent in inbred strains of mice, providing an experimental avenue to systematically map genetic risk factors and to dissect the molecular pathways that orchestrate hearing in peripheral sensory hair cells. Therefore, we ascertained hearing function in the inbred long sleep (ILS) and inbred short sleep (ISS) strains. Using auditory-evoked brain stem response (ABR) and distortion product otoacoustic emission (DPOAE) measurements, we found that ISS mice developed a high-frequency hearing loss at twelve weeks of age that progressed to lower frequencies by 26 weeks of age in the presence of normal endocochlear potentials and unremarkable inner ear histology. ILS mice exhibited milder hearing loss, showing elevated thresholds and reduced DPOAEs at the higher frequencies by 26 weeks of age. To map the genetic variants that underlie this hearing loss we computed ABR thresholds of 63 recombinant inbred stains derived from the ISS and ILS founder strains. A single locus was linked to markers associated with ISS alleles on chromosome 10 with a highly significant logarithm of odds (LOD) score of 15.8. The 2-LOD confidence interval spans approximately 4 Megabases located at position 54-60 Mb. This locus, termed sensorineural hearing loss 1 (Snhl1), accounts for approximately 82% of the phenotypic variation. In summary, this study identifies a novel hearing loss locus on chromosome 10 and attests to the prevalence and genetic heterogeneity of progressive hearing loss in common mouse strains.