- Research article
- Open Access
Genetical genomic determinants of alcohol consumption in rats and humans
BMC Biologyvolume 7, Article number: 70 (2009)
We have used a genetical genomic approach, in conjunction with phenotypic analysis of alcohol consumption, to identify candidate genes that predispose to varying levels of alcohol intake by HXB/BXH recombinant inbred rat strains. In addition, in two populations of humans, we assessed genetic polymorphisms associated with alcohol consumption using a custom genotyping array for 1,350 single nucleotide polymorphisms (SNPs). Our goal was to ascertain whether our approach, which relies on statistical and informatics techniques, and non-human animal models of alcohol drinking behavior, could inform interpretation of genetic association studies with human populations.
In the HXB/BXH recombinant inbred (RI) rats, correlation analysis of brain gene expression levels with alcohol consumption in a two-bottle choice paradigm, and filtering based on behavioral and gene expression quantitative trait locus (QTL) analyses, generated a list of candidate genes. A literature-based, functional analysis of the interactions of the products of these candidate genes defined pathways linked to presynaptic GABA release, activation of dopamine neurons, and postsynaptic GABA receptor trafficking, in brain regions including the hypothalamus, ventral tegmentum and amygdala. The analysis also implicated energy metabolism and caloric intake control as potential influences on alcohol consumption by the recombinant inbred rats. In the human populations, polymorphisms in genes associated with GABA synthesis and GABA receptors, as well as genes related to dopaminergic transmission, were associated with alcohol consumption.
Our results emphasize the importance of the signaling pathways identified using the non-human animal models, rather than single gene products, in identifying factors responsible for complex traits such as alcohol consumption. The results suggest cross-species similarities in pathways that influence predisposition to consume alcohol by rats and humans. The importance of a well-defined phenotype is also illustrated. Our results also suggest that different genetic factors predispose alcohol dependence versus the phenotype of alcohol consumption.
The term "genetical genomics" is entering the common parlance of researchers, denoting the combined use of genetic marker information and transcriptome analysis. Complex trait phenotyping can be fruitfully combined with genetical genomic analysis to ascertain the candidate genes and gene product interaction pathways which significantly influence the variation in expression of a phenotype of interest. We and others have utilized the genetical genomics and phenomic approaches to identify genes and pathways important in genetically influenced complex traits such as obesity, respiratory function and "addictive" behavior [1–3]. In the area of addictive behavior, our laboratory has focused on certain endophenotypes, including the genetic contributors to acute functional tolerance to ethanol [4, 5] and alcohol preference in mice . These studies produced results that implicated the protein products of genes important in learning and memory as contributors to acute functional tolerance to ethanol, and proteins important in orosensory systems and information processing as being important in alcohol preference in mice.
The genus/species Rattus norvegicus has been used extensively in studies of addictive behavior, not only with regard to ethanol, but also to cocaine, opiates, cannabinoids, nicotine, etc. [6, 7]. Although several research groups have suggested that similar or identical biochemical systems in particular brain areas mediate self-administration of addictive drugs [6, 8–10], such proposals have not been tested by applying the unbiased approach of genetical genomics, combined with phenotyping, to studies with rats. In the present study, we have used the HXB/BXH recombinant inbred rat strains , which represent a unique resource for alcohol research. These rats have previously been used for quantitative trait locus (QTL) analysis of various cardiovascular phenotypes and metabolic and behavioral traits . We have generated brain transcriptome data for 28 recombinant inbred strains of HXB/BXH rats , and have combined these data with genetic marker data for these animals . We have also tested the rats in a standard procedure for measuring alcohol preference  and these data have allowed us to generate insight into a complex of genes and their protein products which have significant association with the trait of alcohol preference/consumption in rats.
We have also taken the opportunity to directly compare the candidate genes and pathways for alcohol consumption that we identified using the rodent model, to candidate genes in human populations. We used the "Addictions Array"  to assess the relationships among a panel of genetic markers (SNPs) for "alcoholism and addiction" candidate genes, and the phenotype of alcohol consumption, in two populations of humans that had been characterized in the WHO/ISBRA Study on State and Trait Markers of Alcohol Dependence . The results of this analysis demonstrate a convergence of the human and animal results and provide a differentiation between the genetic polymorphisms associated with predisposition to alcohol drinking and the susceptibility to alcohol dependence.
Male rats from well-characterized HXB/BXH recombinant inbred (RI) rat strains  were used for these studies. The rats were rederived into and maintained in a colony at the University of California, San Diego. These rats were developed from an intercross between two inbred strains, the Wistar origin spontaneously hypertensive rat (SHR/Ola) and Brown Norway congenic (BN-Lx/Cub) . The RI strains replicate the F2 generation of the intercross and have been brother-sister mated for more than 90 generations. Development of the RI strain set utilized gender reciprocal crossing, providing two strain sets that differ in the source of mitochondrial DNA and the Y chromosome. The HXB rats carry mitochondrial DNA from the SHR rats and the Y chromosome from the BN-Lx rats, while the BXH strains are the reverse.
Data on alcohol consumption were gathered on 23 HXB/BXH strains and the two progenitor strains. The number of rats per strain ranged from 9 to 12, with 242 total rats being utilized to measure alcohol consumption. In the first week (week 0) of treatment, rats were given 10% ethanol as their only choice of fluid. For the next seven weeks (week 1 - week 7), the rats were given a choice of two bottles, one with water and one with a 10% (v/v) ethanol solution . Volumes of water and alcohol solution consumed were measured by weight on Monday, Wednesday and Friday of each week. The placement of the two bottles was varied after each measurement to avoid a place preference effect. For the current study we used alcohol consumption data from the second week of the two-bottle choice paradigm. These data were chosen to reflect a stable level of alcohol consumption/"preference" and to facilitate comparisons with alcohol consumption/preference in rats selectively bred for this trait , which use a similar paradigm for selection. The stability of this measure was confirmed by comparison of the mean levels of alcohol consumption for each strain during week 2 and week 7 of the two-bottle choice procedure. A correlation analysis resulted in an r2 = 0.62 (P < 0.001).
Behavioral QTL Analysis
Behavioral QTLs (bQTLs) were calculated for the alcohol consumption phenotype using strain means for average daily alcohol consumption in grams per kilogram during week 2 of the two-bottle choice period. Individual values were not included if they were more than two standard deviations from the strain mean. The recently developed STAR Consortium SNP set of genetic markers was used in this QTL analysis . There were 21 recombinant inbred strains plus the two progenitor strains that had both alcohol consumption data and SNP data. Using these 23 strains, 962 unique strain distribution patterns were identified in the STAR Consortium SNP set. Strain means for alcohol consumption were used in a marker regression QTL analysis. QTLs with empirical P-values less than 0.05 [based on permutation, ] were considered significant, and markers with logarithm of odds (LOD) scores above 2.0 were considered "suggestive" and were included as potential bQTL for alcohol consumption. The 20 Mb region around the significant and suggestive markers was used as the bQTL interval. bQTL analyses were conducted using the R/qtl package in R .
Naive (non alcohol-exposed) male HXB/BXH RI rats were used for microarray analysis. Rats were group-housed and individual rats were quickly anesthetized with isofluorane/air and decapitated according to a protocol approved by the UCSD IACUC. Brains were rapidly removed, sectioned sagittally into two hemispheres, and frozen on dry ice or in liquid nitrogen. The right half-hemispheres were promptly shipped on dry ice to the University of Colorado, where they were kept at -80°C until used for RNA extraction. CodeLink Rat Whole Genome Bioarrays were obtained from G.E. Healthcare/Amersham Biosciences (Piscataway, NJ).
Total RNA Extraction
RNA from the right halves of brains of five to seven rats per strain (26 RI strains and the progenitor strains) was used for these experiments (our prior studies and those of others  showed no statistically significant differences in gene expression between left and right halves of rodent brain). Brains were homogenized in lysis buffer using a Polytron (Brinkmann Instruments, Westbury, NY, USA), and total RNA was isolated using the RNeasy Midi Kit (Qiagen) (Valencia, CA, USA) following the protocol supplied by the manufacturer. An additional clean-up of total RNA was carried out using the RNeasy Mini Kit (Qiagen). Quality of the samples was assessed on an Agilent Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA).
Expression Analysis Using CodeLink Rat Bioarrays
Each CodeLink Rat Whole Genome Bioarray used in these experiments has ~34,000 probes for rat transcripts and expressed sequence tags (ESTs), plus a number of positive and negative control probes. Using the protocol supplied by the manufacturer, double-stranded cDNA was synthesized from the total RNA (4 μg of total RNA from an individual rat brain was used for each array) and was used to obtain biotin-labeled cRNA by an in vitro transcription reaction. Biotin-labeled cRNA was recovered using the RNeasy kit (Qiagen) and concentration and purity were assessed spectrophotometrically. Biotin-labeled cRNA from each rat was then fragmented and hybridized with an individual CodeLink Bioarray. The Bioarrays were subsequently stained for 30 min with Cy5-streptavidin and washed before scanning. Arrays were scanned using a GenePix 4000B Scanner [Axon Instruments (Molecular Devices Corp), Union City, CA, USA] and the images were quantitated with CodeLink Expression Analysis (v 4.1) software (Amersham Biosciences).
Raw intensity values were obtained from the CodeLink processing software. Entire arrays were examined for quality control purposes both before and after normalization. Before normalization, arrays were examined by boxplots and coefficient of variation (CV) plots for consistency, along with background levels and proportion of probes within each category for the CodeLink quality flags. After normalization, samples were examined using hierarchical clustering and pair-wise scatter plots to identify samples that severely deviated from other samples within the same strain (such samples were eliminated from further analysis).
In preparation for normalization, probes were removed from the datasets if they were one of the negative or positive controls placed on the array by the manufacturer. Next, individual values were eliminated based on the quality flags assigned by the CodeLink Expression Analysis Software. Values were eliminated if they were flagged as M (spot was identified to be defective through image inspection at manufacturing), C (spot has a high level of background contamination), I (spot has an irregular shape), or S (spot has a high number of saturated pixels). Values were retained if they were flagged G (spot is good) or L (spot is below local background noise). In addition, to be able to take the log base 2 transformation of the background adjusted intensity values, all background adjusted intensity values below zero were replaced with the value 0.00001. The data were then normalized using a cyclic LOESS (locally weighted scatterplot smoothing) procedure executed in R which accounted for the missing intensity values. Normalization and quality control were executed in R using the affy, codelink, and limma packages .
Identification of Candidate Genes for Alcohol Consumption
After normalization of the expression data, several filters were applied in order to identify candidate genes for alcohol preference in the HXB/BXH RI strains.
Probes were eliminated from further analysis if they did not represent a valid gene. Each probe is represented by a 30-mer base sequence. These 30-mers were compared to the most recent (November 2004) assembly of the rat genome using the BLAST-like alignment tool (BLAT) available from the UCSC website . A probe was retained if at least 28 of the bases matched a region in the genome with less than a two base gap within the sequence. This region also had to fall within an established exon for a RefSeq rat gene, or within an exon sequence for a RefSeq mouse or human gene whose entire sequence matched a region in the rat genome.
Expression QTL analysis
Expression QTL (eQTL) were established using the HXB/BXH RI brain gene expression data. For this analysis, 26 strains (2 progenitor strains and 24 RI strains) were represented among 124 samples (data from two strains for which genotyping was not available were not included). The quality control measures and normalization procedures described above were utilized. The STAR marker set described in the behavioral QTL methods section was also utilized in this analysis. Because more strains were included in this analysis than in the bQTL analysis, 1,184 unique strain distribution patterns were represented in this marker set.
To calculate eQTLs, mean expression levels within strains were used as phenotypic values in a QTL analysis implemented in QTLReaper, which is written in C and compiled as a Python module. A weighted marker regression analysis was used within this program to calculate likelihood ratio statistic (LRS) scores for each marker. LRS scores were transformed to LOD scores for convenience by dividing by 4.61. The regression is weighted to account for the different number of arrays within strains used to calculate strain means. The weight is based on the repeatability of the transcript intensity and number of arrays used to calculate the strain mean . The empirical P-value with respect to the maximum LOD score was calculated for each transcript by permutation . The number of permutations per transcript was increased until the maximum LOD score from the true data was no longer in the top ten of LOD scores from the permutation, or until 1,000,000 permutations were calculated. For eQTLs with empirical P-values less than 0.10, 95% confidence limits for location were calculated using 1,000 bootstrap samples . Probes were retained if they had a significant (P < 0.05) or suggestive (P < 0.10) eQTL.
eQTL/bQTL overlap filter
Probes whose eQTL 95% confidence interval for genomic location overlapped the chromosomal location of any of the four behavioral QTLs were retained for further consideration.
A broad-sense heritability was calculated for each probe in the expression data used in the correlation analysis. Probes were retained if they had a broad sense heritability ≥ 50%.
Detection Limit Filter
Detection limits are calculated by the CodeLink software as the median intensity of the local background plus 1.5 local background standard deviations. Probes were eliminated if values were below detection limits for ≥ 50% of the samples.
The correlation of alcohol consumption and gene expression was modeled using a joint mixed model for both the expression and alcohol consumption data, Y ijk = μ k + s ik + ε ijk for i = 1,...,25 strains, j = 1,...nik rats, and k = 1,2 where k = 1 indicates a gene expression measure and k = 2 indicates an alcohol consumption measure. The εijk are independent and identically distributed from a normal distribution with a 0 mean and a measure-specific variance (σk 2). The sik are random strain-measurement effects. The pair, [si1, si2], are independent and identically distributed as a bivariate normal with an unstructured covariance where a separate between-strain variance is calculated for expression and alcohol consumption plus a covariance between expression and alcohol consumption. This model was executed in SAS version 9.1 using a maximum likelihood model in the linear mixed model function called PROC MIXED. A likelihood ratio statistic was used to determine if the covariance between gene expression and alcohol consumption is significantly different from zero. Probes which demonstrated an overlap of their eQTL location and the genomic location of one of the four bQTLs and which passed the other filters listed above, were used in our correlation analysis. Probes were retained if they had a significant covariance with alcohol consumption (P < 0.05).
Proportion of Variance Explained and Model Building
Once a list of candidate genes was established, the proportion of genetic variance in alcohol consumption among strains that was explained by gene expression was examined. A multivariate genetic model for alcohol consumption was constructed using a linear regression with strain means for both gene expression and alcohol consumption. Genes were selected for inclusion using a forward stepwise selection process with a significance criterion of 0.01 for entrance into the model. A forward selection method was chosen over a backward selection method in this instance because of the large number of potential covariates (candidate genes). The proportion of genetic variance explained was determined for each candidate gene individually and for the combination of transcripts in the final multivariate genetic model which described the strain variation in alcohol consumption.
Functional Analysis of Candidate Gene Products
In order to determine the functions of the candidate gene products, and assemble them into pathways, a literature search was performed using PubMed. The gene symbols (NCBI) and gene product names were used as key words.
In order to replicate differences in gene expression associated with alcohol consumption, gene expression profiles were analyzed using Affymetrix GeneChip Rat Gene 1.0 ST arrays (Affymetrix, Santa Clara, CA, USA), which utilize probesets that span the coding region of each gene. RNA from brains of five male rats from HXB RI strain 23 (low alcohol consumption) and five male rats from HXB RI strain 26 (high alcohol consumption) was used for these studies. These particular strains were chosen because the expression levels of all of the candidate genes (assayed on the CodeLink arrays (Amersham Biosciences)) showed an appropriate direction of differential expression, compared to the correlation analysis, although not all of the differences in the candidate gene expression were statistically significant between these particular strains. RNA (0.3 μg per rat) was processed and hybridized to the Affymetrix arrays (cDNA from one rat to each individual array) according to the manufacturer's protocol. Arrays were labeled and scanned as previously described . For data analysis, raw data and the normalized probeset data were assessed for quality using the 'aroma.affymetrix' package in R  and Expression Console from Affymetrix . All ten arrays passed quality control. Raw intensity values for individual probes were corrected for background using the RMA correction, probe values were quantile normalized, and probes were summarized into probesets  in aroma.affymetrix. The probe intensities were compared between the two strains using a two-sample t-test assuming unequal variance on both the probeset and the probe level. Probe level data were only examined for those probes that fell within the same exon, upstream, or downstream region as the original CodeLink probe.
Human Association Study
Subjects for the human genetics association study were recruited as part of the WHO/ISBRA Study on State and Trait Markers of Alcohol Use and Dependence . More specifically, phenotypic data from Caucasian subjects from Montreal, Canada and Sydney, Australia were analyzed since these subjects were also genotyped using the Addiction Array . The phenotypic data were collected using a structured interview instrument that was developed for the WHO/ISBRA Study. This instrument was developed in concert with the development of the Alcohol Use Disorders and Associated Disabilities Interview Schedule (AUDADIS), and included eight major sections [see ]. Of particular relevance for the current study, the interview included a section that gathered information on beverage-specific frequency and quantity of alcohol drinking during the past 30 days, and the symptoms experienced with various levels of drinking. Another section gathered information that allowed for International Classification of Diseases -10 (ICD-10) and Diagnostic and Statistical Manual -IV (DSM-IV) diagnosis of alcohol abuse, dependence and withdrawal syndromes. In addition, family history of alcohol problems or alcohol dependence, major depression, drug use and abuse and other psychiatric problems were assessed (see Additional file 1, Table S2). Additional information on the demographic characteristics of the human subjects is contained in Additional file 1, Table S3. For more details on the original study and data collection, see Glanz et al. .
Genotype data were obtained from a custom array using the Illumina Goldengate SNP technology assay platform. This array includes a panel of markers designed to extract full haplotype information for 130 candidate genes associated with alcoholism, other drug addictions, and mood and anxiety disorders. The design and performance of this array are described in detail by Hodgkinson et al. . Briefly, the array contains 1,350 tagSNPs representing the 130 addiction-related genes, and an additional 186 markers identified as highly informative for ancestry (AIMs).
The main phenotypic outcome measure of this analysis was alcohol consumption reported in grams per kilogram of body weight per day. The total amount of alcohol consumed in the last 30 days was calculated based on information that was collected in the WHO/ISBRA questionnaire on beverage-specific frequency and quantity of drinking and this total was divided by 30 to obtain the daily average.
Prior to the univariate analyses, the SNP dataset was subjected to strict quality control standards. Individual SNPs were eliminated if they were not in Hardy-Weinberg equilibrium (FDR<0.01) or if they had minor allele frequency less than 5%. Subjects were eliminated if their genotype was not determined for at least 80% of the SNPs according to standards set forth in Hodgkinson et al. . The ancestry informative markers were not considered in any of the association analyses with the exception of the tests for stratification.
Since the subjects came from two distinct geographical locations, it is reasonable to consider that there may be some genetic differences between these populations. Population stratification with respect to the outcome, alcohol consumption, was tested against the AIMs using the concept outlined in Pritchard and Rosenberg . A Fisher's combined probability test was used to calculate significance, due to the continuous nature of the outcome. Differences in allele frequencies between the two populations for the AIMs were also tested using Pritchard and Rosenberg's  method of combining chi-square statistics.
The association between alcohol consumption and each individual SNP was tested in an ANOVA using a genotype model. P-values were adjusted using a false discovery rate (FDR). SNPs with an FDR < 0.05 were considered to show a significant association.
Haplotype blocks were determined in Haploview  for genes with significant SNPs, using the Gabriel et al.  method for identifying block boundaries. Haplotype block structure was examined in the WHO/ISBRA population that had the significant SNP association, and in the CEPH (Utah Residents with Northern and Western European Ancestry) population in the HapMap project  for comparison. The haplotype block structure in the HapMap population was used to construct haplotypes for individual subjects. Subjects were assigned haplotype pairs and the posterior probabilities of those haplotype pairs were determined using haplo.em in the haplo.stats package of R. If haplotypes had less than five counts, they were recoded to the closest haplotype. If the closest haplotype differed by more than one SNP, data related to the minor haplotype were deleted. All further analyses were weighted by the posterior probability of that particular haplotype pair. Additive, recessive, and dominant representations for all haplotypes were examined in a forward selection model building procedure where the criterion for entrance into the model was P < 0.01.
Multivariate Genetic Model
A backwards selection linear regression was used to create a multivariate model that potentially could include effects for all genes that were significant in the univariate analysis. To start the model selection procedure, the best fit univariate model for each gene, either haplotype-based or SNP-based, according to the Bayesian Information Criterion [BIC, 30], was entered into the model. An exit criterion of P > 0.01 was used in the model building process.
Multivariate Model With Covariates
A multivariate model was also created using the candidate genes and the covariates shown in Additional file 1, Table S2. A backwards selection model building process was used, with a more conservative exit criterion of P > 0.001.
Association Between Covariates in Final Multivariate Model and Candidate Genes
Each candidate gene was tested for association with each of the covariates that remained in the final multivariate model. Each gene by covariate pair was tested independently. Genes were represented using their best fit univariate model as described above. For the binary covariates, a Fisher's exact test was used for association. For continuous covariates, an ANOVA model was used for association. All P-values were adjusted (within population) for multiple comparisons using an FDR .
Alcohol Consumption by HXB/BXH RI Rat Strains
Figure 1 shows the average daily alcohol consumption for each of the HXB/BXH RI strains and the two progenitor strains. Alcohol consumption varied among strains, and it was calculated that 49% of the variance could be attributed to strain (genetic variance).
Four genomic regions were identified as having significant or suggestive associations with alcohol consumption in the HXB/BXH rats (Figure 2). The most significant marker was on chromosome 1 at 224.7 Mb, and it had a LOD score of 4.02 (empirical P-value = 0.05). This particular QTL and two other regions had multiple adjacent markers with a LOD score above 2. The bQTL region in each case was extended 10 Mb on either side of the adjacent markers which had LOD scores above 2.0 (Figure 2).
Gene Expression Data and Filtering of Gene Expression Data from HXB/BXH RI Strains
Quality Control and Normalization
After quality control http://Phenogen.UCHSC.edu, 121 arrays remained from the 25 strains used for the alcohol consumption and gene expression correlation analysis. We eliminated 497 control probes from further analysis along with a total of 31,379 (<1%) individual expression values due to quality control flags. In addition 45,722 expression values (1% of remaining values) were changed to 0.00001 because original (local background-corrected) spectroscopically derived values were less than 0. A total of 33,849 probes on each of the 121 arrays were included in the further analysis (filtering and correlation analysis).
As noted in methods, a number of filters were used to identify a list of potential candidate genes for alcohol preference in the HXB/BXH RI rats. Starting with the quality controlled, normalized dataset, we initially determined probes that could be matched with exonic regions of the rat genome (or with a homologous human or mouse sequence). There are 14,622 unique rat RefSeq DNA IDs specified by the UCSC Genome Browser. Of these, 8,433 genes (10,162 probes) are included on the CodeLink Whole Genome Rat Array. We identified 4,966 additional probes on the CodeLink array that were associated with either a human or mouse homolog gene (RefSeq DNA ID) whose sequence matched an area of the rat genome. Of these 15,128 probes, 12,466 (82%) were present (by our criteria) in the brain samples that we analyzed. Significant or suggestive eQTLs could be identified for 2,164 probes (17%), and, of these, 378 probes had an eQTL that overlapped one of the bQTLs for alcohol consumption in these rats. Transcripts with eQTLs that overlap bQTLs were considered to represent likely candidate genes for the alcohol consumption behavior . The heritability filter was used to insure that the expression levels of genes in the candidate gene list had a high heritability, and 308 probes met the criterion of a broad sense heritability ≥ 0.5. We also applied a filter based on detectable levels of expression in the rat brains. Only transcripts that were expressed in brains of ≥ 50% of the individual samples were included. Although this choice is arbitrary, we did not wish to include genes with low or undetectable levels of brain expression in a large number of strains, as such transcripts could produce unreliable results. This filter was passed by 304 probes. The final filter was correlation of probe expression with alcohol consumption. By initially applying all of the previous filters, the risk of false positives in this correlation step was reduced. Of the 304 starting probes, 20 had a significant (P < 0.05) correlation with alcohol consumption in the HXB/BXH RI rat strains. These probes are listed in Table 1.
The chromosomal locations, as well as the location of the bQTLs and eQTLs for the 20 probes, are shown in Table 1. Nine of the candidate genes were posited to be cis-regulated, i.e., the chromosomal location of the gene is within 5 Mb of the eQTL location. The expression of the rest of the candidate genes was posited to be regulated from an area in the genome distant from the physical location of the gene (trans-regulation). The expression of seven of the trans-regulated genes was controlled from an eQTL on chromosome 6. This finding suggests the possibility that there is, for example, a polymorphic gene coding for a transcription factor within this eQTL region, and that this transcription factor contributes to the regulation of the expression of all seven of the candidate genes. Although 18 putative transcription factor genes could be identified in the eQTL region on chromosome 6 (see Additional file 1, Table S1), there is not sufficient information available on the exact nature of the binding sites for these transcription factors to determine whether such binding sites are overrepresented in the promoter regions of the trans-regulated genes. Table 1 also shows the correlation coefficients for the correlation with alcohol consumption, and the results of the heritability and expression level filters. The expression levels of 13 of these probes were positively correlated with alcohol consumption, and expression levels of seven were negatively correlated.
The transcript that explains the greatest proportion of genetic variance in alcohol consumption in the HXB/BXH RI strains is P2rx4 (r2 = 0.39, model P = .0009). When we performed a forward stepwise regression analysis to determine a multivariate model for alcohol consumption, the final model contained P2rx4 and Tmem2. The model explained 57% of the genetic variance among strain means for alcohol consumption.
Replication Study of Candidate Gene Expression in HXB/BXH Rat Strains
The results obtained using the CodeLink arrays were assessed using a different gene expression measurement platform. The analysis of the Affymetrix Array data for the expression levels of the candidate genes in brains of rats from two HXB/BXH strains are shown in Table 2. These two strains were chosen based on their differences in voluntary alcohol consumption (see Figure 1) and their differences in brain candidate gene expression levels determined using the CodeLink arrays. Table 2 includes data from the CodeLink and Affymetrix arrays for these two strains to allow comparisons. In evaluating the replicability of results, it is necessary to keep in mind that even though the candidate genes displayed a significant correlation with alcohol consumption when all RI strains were considered, not all of the candidate genes (based on CodeLink results) displayed significantly different expression levels in these two particular strains. However, the expression differences in these strains were all in the same direction as the overall correlation. Therefore, in comparing to the results for the Affymetrix arrays, we took into account both the significance and direction of expression differences. Three candidate genes displayed significant differences in expression, in the same direction, using both arrays: Afap1l1, Cckbr and P2rx4. Another six candidate genes showed significant differences using CodeLink arrays, and the same direction of differences on the Affymetrix arrays (Tnk2, Tubb6, Ptpre, Mogat2, Nek3 and Tmem2). Three of the candidate genes did not show significant differences in expression using the CodeLink arrays, but either did show a significant difference in expression levels in the same direction (Mc4r), or showed non-significant differences in expression in the same direction using the Affymetrix arrays (F2rl1, Abat). Based on this analysis, we placed more confidence in these 12 of the 20 original candidate genes. The rest of the candidate genes did not show significantly different levels of expression using the Affymetrix arrays, and showed differences in expression in opposite directions on the two platforms. Since the expression values obtained with CodeLink arrays are based on one probe, while the expression values from the Affymetrix arrays represent a summary from multiple probes, spread throughout each gene, we also compared results for probes on the Affymetrix arrays that correspond in location to the probes on the CodeLink arrays. In most cases, the results for Affymetrix individual probes and probesets were the same. For Tnk2, the result from the individual probe on the Affymetrix array indicated a significantly different expression level between strains in the same direction as had been determined using the CodeLink array between the tested strains. We have limited our discussion of candidate genes to those that were found (same direction of differential expression) using both CodeLink and Affymetrix arrays (Table 2).
Human Association Study
Of the original 606 subjects in the WHO/ISBRA data from Montreal, 545 were self-reported Caucasian. Alcohol consumption in the Montreal population ranged from abstainers (0 g/kg/day) to heavy drinkers (maximum, 8.4 g/kg/day; see Additional file 1, Figure S1). Likewise, in the Sydney population, 242 of the original 285 subjects were Caucasian and alcohol consumption ranged from abstainers (0 g/kg/day) to heavy drinkers (maximum, 7.3 g/kg/day; see Additional file 1, Figure S1). Other covariates, and demographic characteristics, were comparable between the two populations with the exception of gender because only male subjects were recruited in Sydney (Additional file 1, Tables S2 and S3). For the purpose of comparability, we report results from only the male Caucasian subjects (n = 280) from the Montreal population.
Out of the 1,350 gene-related SNPs, 98 SNPs were not genotyped in any of the Montreal subjects (98 in Sydney subjects) and 88 were not informative across the Montreal sample (49 in the Sydney sample). Further evaluation of the remaining SNPs resulted in 121 SNPs being eliminated because they were not in Hardy-Weinberg equilibrium (HWE) in the Montreal sample (156 SNPs in the Sydney sample) (FDR<0.01). An additional 203 SNPs were eliminated in the Montreal sample (203 in Sydney) that had a minor allele frequency (MAF) that was less than 5% (46 SNPs in the Montreal population (79 in Sydney) did not meet HWE and had a MAF<5%). In addition, 110 subjects from Montreal and 61 subjects from Sydney were eliminated for having a call rate that was less than 80%, based on the criteria described in Hodgkinson et al. . These exclusion criteria left 840 SNPs on 435 subjects for analysis in the Montreal population (220 males) and 845 SNPs on 181 subjects in the Sydney population.
There was no evidence for population stratification at each site with respect to the outcome, i.e., alcohol consumption (chi-square = 344, df = 352, P-value = 0.62). However, there was substantial evidence for genetic difference between the population gathered in Sydney vs the population from Montreal (chi-square = 318, df = 176, P-value = 3.1 × 10-10). Because there was such strong evidence for genetic differences between the populations, the two populations were analyzed separately.
Genetic and Phenotypic Associations with Alcohol Consumption
Table 3 shows the SNPs that were significantly associated with alcohol consumption in the univariate analysis. In the Montreal population, SNPs representing two unique genes (GAD1, MPDZ) were significantly associated with alcohol consumption. Although the significant SNP in MPDZ only reached statistical significance in the dataset that combined both males and females, its FDR was 0.09 in Montreal males. For consistency, all further analyses in the Montreal population were only conducted on male subjects. In the Sydney population, SNPs representing four unique genes (CHRM5, GABRB2, MAPK1, PPP1R1B) were associated with alcohol consumption.
We determined the haplotype block that contained each significant SNP. For GAD1 in the Montreal sample, three SNPs significantly associated with alcohol consumption were located within the same known haplotype block (Additional file 1, Figure S2). Haplotypes and individual SNPs were evaluated for their association with alcohol consumption to determine the best fit model for each gene.
The final genetic multivariate model for Montreal male subjects included the recessive effect of haplotype 2 (see Additional file 1, Table S4) in MPDZ and the SNP genotype effect for rs701492 in GAD1 (Table 4). In combination, these two genetic factors explained 16% of the total variance in alcohol consumption. For the subjects from Sydney, three of the candidate genes were represented by SNPs in the final genetic multivariate model (Table 4). For the fourth candidate gene, MAPK1, the best fit to the data was a model with both an additive and dominant component for haplotype 3 (see Additional file 1, Table S4); however, this gene was not included in the final genetic multivariate model.
We also analyzed the relationship of phenotypic characteristics with alcohol consumption levels by males in Sydney and Montreal. As expected, our data showed that alcohol consumption was positively correlated with alcohol dependence in both populations (Additional file 1, Figure S3). The multivariate model with covariates for the Sydney population included alcohol abuse and the diagnosis of alcohol dependence (DSM-IV) (Table 4). In the population from Montreal, the multivariate model with covariates included a diagnosis of major depression and a family history of depression, as well as alcohol abuse and dependence (Table 4). In both populations the models including genetic and phenotypic characteristics accounted for 62% of the variance in alcohol consumption. Our further analysis indicated that none of the genotypes or haplotypes that were found to be significantly associated with alcohol consumption were significantly associated with alcohol abuse or alcohol dependence in either population (FDR >0.05). In fact, we found none of the SNPs represented on the "addiction array" to be significantly associated with DSM-IV alcohol dependence in either the Sydney or Montreal populations (FDR >0.05). A caveat to keep in mind is that the sample sizes may not have provided enough power to allow us to detect associations using the binary outcomes of alcohol dependence or abuse.
The reasons for differential intake of ethanol by individuals in a population, be it humans or other animals, have been the subject of an immense amount of research. The current dogma posits that both environmental and genetic factors contribute to individual differences in alcohol consumption [32, 33]. What is sometimes confusing about the literature on alcohol consumption is the lack of discrimination between alcohol consumption in non-dependent vs alcohol-dependent individuals . The DSM IV and ICD-10 criteria for alcohol dependence focus intensively on dichotomizing alcohol drinking behavior in dependent and nondependent individuals [35, 36] and several recent reviews by authors of this current work , and others , have detailed evidence and hypotheses regarding the progression of events which transform the non-dependent alcohol intake phenotype to the alcohol-dependent alcohol intake phenotype.
Epidemiologic studies do indicate that the higher the levels of alcohol intake by an individual, the higher is the propensity for an individual to become alcohol dependent . In this regard, high levels of alcohol intake become a risk factor for transition to alcohol dependence, but there is no a priori reason to assume that the genetic factors that may be responsible for modulating an individual's non-dependent levels of alcohol consumption are the same factors that predispose or protect an individual from becoming dependent. In other words, a genetic relationship between the propensity for high alcohol intake and the propensity for becoming alcohol dependent has not been demonstrated. In fact, one can interpret some of the data collected with mice to show a dissociation between a propensity for high alcohol preference and propensity for physical dependence. C57BL/6 mice are a genetically inbred strain that has been demonstrated to have one of the highest levels of alcohol preference, but these mice have a low propensity to develop alcohol dependence (exemplified by signs of CNS withdrawal hyperexcitability) [2, 39, 40]. It should thus be emphasized that the phenotype we are investigating in this paper is non-dependent alcohol intake.
What emanates as a conclusion from our genetical genomic/phenomic approach for the search for "candidate genes" that influence non-dependent alcohol intake is that a group of genes/gene products that can be linked to systems which, in the rat (and human), control appetite and satiety, play an important role in variations in the non-dependent alcohol intake phenotype (see 5). It is of interest that our earlier work, which involved a meta-analysis of data from recombinant inbred mice, selectively bred mice and a large panel of inbred mice, indicated that orosensory systems in the mice are critically involved in alcohol selection in a two-bottle choice experimental design . Our two studies, taken together with a number of prior studies , indicate that alcohol selection by rodents involves the sensory and caloric information transduction systems dealing with the recognition, and the satiating and rewarding properties of foodstuffs. Ethanol, unlike other psychoactive drugs, has significant caloric value (seven cal/gm) and these calories have been clearly shown to have relevance to an organism's energy status .
Figure 3 illustrates that most of the candidate genes that we identified can be related within a neurobiological pathway that was constructed based on analysis of the literature. Many of the candidate genes identified in this study code for products (protein) that have been implicated in regulation of feeding and energy metabolism. For example, the cholecystokinin 2 (CCK2) receptor is expressed in hypothalamic regions such as the paraventricular nucleus (PVN), which contains the circuitry necessary for maintaining energy homeostasis (food intake vs energy expenditure). CCK released from the gut acts as a "satiety hormone" , and can activate neurons in various hypothalamic regions . The melanocortin 4 receptor (Mc4r) is also expressed in the hypothalamic regions that regulate energy metabolism, and activation of these receptors by the endogenous agonist, MSH, also inhibits food intake. It has been suggested that these actions involve changes in GABA release from hypothalamic interneurons [45, 46]. Both Mc4r and CCK2R are expressed at higher levels in animals that consume higher amounts of alcohol, suggesting that these animals may have an inherent propensity for reduced food intake and lower energy expenditure.
However, animals and humans ingest food not only for nourishment, but also for the rewarding properties of food, and motivational mechanisms are important for generating responses needed for food-seeking and consumption behavior. The neuronal pathways and neurotransmitter substances that regulate food intake for energy homeostasis, and for the rewarding properties of food, are intimately connected but not identical. Particularly relevant to our current studies are the hypothalamic nuclei with connections to the ventral tegmental area (VTA) that generate and carry information on hunger and satiety, and initiate the cascade of events resulting in "wanting/liking" food or other "rewards" [47, 48]. Much of the hypothalamic input to the VTA which provides information on the energy status of the mammalian organism is initiated in the arcuate nucleus of the hypothalamus and transmitted via the lateral hypothalamus to the VTA dopaminergic neurons. The transmitter substance that has been associated with stimulatory input from the lateral hypothalamus is orexin, which activates the dopaminergic neurons of the VTA . GABA neurons provide the major direct inhibitory input to the dopaminergic neurons of the VTA, but GABA neurons arising in the arcuate nucleus also provide inhibitory input to the orexin neurons residing in the lateral hypothalamus . Thus, GABAergic activity related to food/energy requirements can generate inhibitory influences on the VTA dopaminergic neuron firing either directly or through inhibition of orexin signaling to the VTA.
The products of many of the candidate genes identified in this study can affect GABAergic neuronal activity. One of the most studied systems in this regard, and one suggested to link the homeostatic and reward pathways associated with food intake , is the endocannabinoid system. The endogenous agonists anandamide and 2-AG (2-arachidonoyl glycerol) act primarily through presynaptic cannabinoid 1 (CB1) receptors to inhibit GABA release . The endocannabinoid system in the hypothalamus has been implicated in regulation of both food and ethanol intake [52, 53]. 2-AG can act as an anterograde or retrograde messenger to inhibit GABA release , and can then be inactivated by hydrolysis by monoacylglycerol lipase or fatty acid amide hydrolase . More recently, however, another group of enzymes has been proposed to inactivate 2-AG, not by hydrolysis, but by conversion of 2-AG to diacylglycerol. The enzyme carrying out this reaction is the product of Mogat2 (also referred to as Mgat2)). Mogat2 mRNA was found to be expressed at lower levels in the higher alcohol-drinking rats, possibly resulting in a reduced rate of metabolism of 2-AG, leading to greater inhibition of GABA release.
The cholecystokinin 2 (CCK2) receptor is expressed not only in the hypothalamus, but also in other brain regions, such as the VTA or amygdala, where it is localized presynaptically on GABA neurons [57, 58]. The activation of the CCK2 receptor has been shown to promote action potential-induced GABA release in certain brain areas (e.g., hippocampus) , but more recent data suggest that the initial increase in GABA release produced by CCK via the CCK2 receptors is followed by a reduction in GABA release . Our findings of higher levels of the mRNA for the CCK2 receptor in high alcohol consuming rats could be construed as another differential means of modulating GABA release in low and high-drinking rats. Another candidate gene product, the P2X4 purinergic receptor (the product of P2rx4), is also localized presynaptically on GABA neuron terminals . Activation of P2X4 receptors on GABA terminals may increase GABA release . Since P2X4 receptor mRNA is expressed at lower levels in the high alcohol-consuming animals, this again would point to diminished levels of GABA release in such rats.
Another identified candidate gene may further regulate the excitability of GABA neurons and GABA release. Ptpre encodes the receptor type protein tyrosine phosphatase epsilon, which has been reported to dephosphorylate and decrease the activity of the voltage-gated potassium channel, Kv2.1 . The Kv2.1 channel is found on somata and proximal dendrites of inhibitory (GABAergic) interneurons [see ]. The expression of Ptpre is lower in animals that consume higher amounts of alcohol, and if this results in a lesser level of phosphatase activity and higher levels of the channel activity, one would expect lower excitability of GABAergic interneurons in these animals.
On the other hand, another of the candidate genes, Abat, codes for the enzyme 4-aminobutyrate (GABA) aminotransferase, which degrades GABA. This gene is expressed at lower levels in animals with higher ethanol consumption, and this situation could contribute to enhanced GABA levels in GABAergic synapses in the high alcohol-consuming rats.
The effects on GABAergic activity described above could contribute to a decreased hypothalamic GABAergic inhibitory effect on VTA dopaminergic activity in the high alcohol-consuming rats. Furthermore, the products of a number of the candidate genes can more directly affect the activity of midbrain dopamine neurons. MSH, administered into the VTA, increases dopamine release in the nucleus accumbens via actions on VTA melanocortin 4 receptors . Interestingly, administration of a melanocortin 4 receptor agonist to mice (icv) reduced alcohol drinking , and administration of a melanocortin 3/4 receptor agonist to alcohol-preferring rats reduced alcohol intake in association with reduced food intake . The CCK2 receptor is also found in the VTA  and likely responds to CCK release from GABAergic neurons synapsing directly with dopaminergic cell bodies in this region. Administration of CCK into the VTA results in increased dopamine release in the nucleus accumbens, mediated by activation of the CCK2 receptors in the VTA [57, 66]. Taken together with the inherent differences in expression of genes whose products influence GABA release, the greater expression of the CCK2 and MC4 receptors, if it occurs in the VTA of the high alcohol consuming animals, generates a picture of lower GABAergic inhibitory tone and an increased propensity for direct dopaminergic neuron activation in the VTA by MSH and CCK.
Recent work has shown that the quantities of ethanol consumed by rats in a two-bottle choice paradigm, can, in fact, enhance the electrophysiological activity of VTA dopaminergic neurons . What is unclear about such studies is whether these effects of ethanol are direct or indirect, through modulation of transmitters influencing dopamine neuron firing. For example, acutely, ethanol inhibits the function of the P2X4 receptor , and this action may reduce the presynaptic function of the P2X4 receptor and reduce GABA release, thus disinhibiting dopaminergic activity. Xiao and Ye  have also reported that ethanol can acutely inhibit GABAergic neuron activity in the VTA via other indirect mechanisms, while enhancing such activity in other brain areas. These actions of ethanol in the VTA would tend to disinhibit dopaminergic neuron activity and synergize with the inherently lower GABAergic tone suggested by our studies.
Another intriguing aspect of our candidate gene list is the presence of a number of genes coding for products that can affect postsynaptic aspects of GABA receptor trafficking. These candidate genes and gene products are related to rho GTPases, which affect the cytoskeleton, membrane trafficking and cell adhesion through their regulation of actin dynamics. For example, Cdc42 is a member of the rho GTPase family that can stimulate actin polymerization, affecting specific steps of vesicle trafficking such as those involved in endocytosis . Cdc42 can be activated by a protein called collybistin, which is a gephyrin-binding protein . Gephyrin is a scaffolding protein for GABAA and glycine receptors , and these interactions provide support for a specific role for Cdc42 in clustering and trafficking of GABAA receptors.
Two of the candidate genes that we identified are Cdc42-associated proteins. Tnk2 codes for a cytosolic non-receptor tyrosine kinase, also known as Ack (activated Cdc42-associated tyrosine kinase). This kinase has been identified as a downstream effector of Cdc42 which is important for regulation of receptor degradation via endocytic mechanisms . It is widely expressed in brain, including various nuclei of the mesencephalon . F2rl1 codes for a protease-activated receptor known as Par2. This G protein-coupled receptor is also widely expressed in brain, including relatively high expression in amygdala, and is found in the soma and dendrites of GABAergic and glutamatergic neurons [74, 75]. Par2 agonists activate several cell signaling molecules, including phosphoinositide 3-kinase. This latter activity requires Cdc42 activation by Par2 . The negative association of the expression of Ack with alcohol preference, and the positive association of Par2 expression, could enhance the endocytic recycling process associated with Cdc42-mediated actin polymerization in the animals that consume higher amounts of alcohol, leading to more rapid GABA receptor desensitization.
The products of some other candidate genes are associated with actin organization and microtubule activity, which is also important for the trafficking of receptors and other proteins. Afap1l1 codes for an actin filament associated protein. The expression of this gene is positively correlated with alcohol consumption. Microtubules are involved not only in receptor transport, but also in the movement of gephyrin to and from the synaptic plasma membrane . Dynein is a microtubule motor that participates in axonal transport of neurofilaments , and that interacts with beta tubulin , the protein product of another of the candidate genes, Tubb6. The expression of all of the above-mentioned candidate genes is associated positively with alcohol consumption in the RI rat strains, supporting the postulate that there may be enhanced GABAA receptor trafficking (endocytosis), in the brains of animals that have a propensity to consume more alcohol.
Overall, our analysis of the functional pathways defined by the identification of the candidate genes in the HXB/BXH rat strains focuses attention on GABAergic and dopaminergic activity that may set a tone that predisposes to (or against) voluntary alcohol consumption. We are proposing that, in rats, a lower inherent GABAergic tone generated by reduction in presynaptic release and a more responsive GABA receptor desensitization system predisposes to higher alcohol consumption in a free choice experimental paradigm. There are numerous studies showing that administration of agonists or antagonists of GABA or dopamine receptors can alter alcohol consumption or self-administration by rats [e.g., [80–90]].
Although we have focused attention on neuronal systems which mediate the animals' recognition of the energy status of the body and the rewarding properties of caloric substances, it should be clearly stated that the candidate genes we identified, and their products, do play important roles in other anatomically-defined neural systems. A case in point may be that dopaminergic neurons which innervate the dorsal striatum (cell bodies in the substantia nigra) may have as much to do with appetitive behavior as the systems we describe innervating the nucleus accumbens from the A10 nucleus (VTA) . Another area of brain where a number of the products of the candidate genes that we identified interact is the amygdala. For instance, Mc4r is expressed in several nuclei of the amygdala and in the bed nucleus of the stria terminalis [92, 93]. Studies with melanocortin 4 receptor null mutant mice showed that restoration of the receptor in a population of amygdala neurons reduced food intake , similar to the effect of activation of this receptor in hypothalamus. Cckbr is expressed in the amygdala as well, and CCK increases GABA release in amygdala via the CCK2 receptor . Higher levels of GABA release in the amygdala may set a tone that predisposes to higher alcohol consumption, given the finding that delivery of GABAA receptor antagonists to the amygdala reduced ethanol self-administration . The endocannabinoid system can also interact with the CCK system in the amygdala, as well as modulating GABA release in this area of brain. CB1 receptor agonists inhibit the release of CCK from GABAergic neurons and Chhatwal et al.  have shown that the effects of cannabinoids on the extinction of conditioned fear responses are mediated through CCK2 receptors. Clearly the neuronal systems that affect conditioned fear (anxiety) responses are also involved with appetitive and reward pathways [97, 98].
The ultimate goal of genetic studies using animal models of alcohol consumption is to identify candidate genes that influence a human's level of alcohol consumption. It should be stressed at this point, that the phenotype we utilized for our genetic association studies with humans was also the quantitative measure of alcohol intake. This phenotype was chosen to allow proper comparison with our studies with the rats. We also performed a separate analysis of the genetic association of the phenotype of alcohol dependence defined by both DSM IV and ICD-10 criteria with the panel of 1,350 SNP marker included on the "Addiction Array" . The polymorphisms that we found to be associated with the quantitative alcohol consumption phenotype in human populations, identified genes whose products are involved in the pathways determined from our studies of differential levels of alcohol intake by the panel of recombinant inbred rats. Particularly evident were components of the GABAergic signaling pathway, including the β2 subunit of the GABAA receptor in the Sydney population, and the GABA synthetic enzyme, GAD1, as well as the MPDZ protein, which can act as a scaffolding protein for the GABAB receptor , in the Montreal population. MPDZ has previously been linked to alcohol withdrawal seizure susceptibility in mice , but not with levels of alcohol intake. The other genes identified in the Sydney population have also been linked to mesolimbic dopaminergic activity. CHRM5, the M5 muscarinic cholinergic receptor subtype, localized on VTA dopamine neurons, is thought to contribute to tonic excitatory regulation of dopamine transmission . As noted earlier, our studies do indicate that the products of the candidate genes predisposing high alcohol drinking by rats can be linked to the modulation of VTA dopamine neuron function. Additionally, muscarinic cholinergic receptors in the VTA have been shown to play an important role in alcohol selection by rats selectively bred (P rats) for high alcohol intake . PPP1R1B, also known as DARPP-32, as well as MAP kinase (the product of MAPK1 gene), are downstream targets of dopamine D1 receptors in the medium spiny neurons of the shell of the nucleus accumbens and other brain areas receiving input from the VTA. Although the "Addiction Array" SNPs did not allow us to assess polymorphisms in a number of the candidate genes identified in our studies with rats (e.g., Mc4r), it has been demonstrated that activation of melanocortin 4 receptors in the nucleus accumbens could enhance dopamine D1 receptor-mediated cyclic AMP production , and modulate the enzymatic activity of DARPP-32. The Mc4r mRNA, which we found to be elevated in brains of the alcohol-preferring rats, is reported to be found in the nucleus accumbens, apparently in medium spiny neurons, as well as in the VTA and hypothalamus of humans [92, 103].
The overall impression that is generated by cursorily examining the results of our studies with rats and humans, and our prior studies with mice, is that little evidence may have been produced to indicate that identical genes or gene products predispose free choice alcohol intake in rodents or humans. What may be missed is the fact that certain identical neurobiologic pathways have been identified in all of these investigations. In rodents, one can posit that neurobiologic systems that participate in sensing and transducing information about the rewarding or aversive properties of foodstuffs play an important role in oral consumption of ethanol in a free-choice paradigm. Polymorphisms in the loci of genes involved in such an appetitive pathway are also associated with quantitative measures of alcohol intake in humans [104, 105]. The proposition that appetitive sensory systems and/or caloric qualities of ethanol contribute to ethanol drinking by rodents is not novel [41, 106–108], but the identification of such a relationship by an unbiased genetical genomic/phenomic technique allows one to progress to examination of specific pathways and genetic mechanisms in future studies.
Our results also cast some light on the influence of genetic polymorphisms on levels of alcohol drinking vs alcohol dependence in humans. Epidemiologic studies suggest a strong correlation between levels of ethanol consumed and the diagnosis of alcohol dependence , and alcohol intake has been taken by some as a surrogate (endophenotype) for alcohol dependence . In our study, the polymorphisms that we found to be associated with the level of alcohol consumption by humans continued to have an influence in the multivariate model, even when current alcohol dependence and alcohol abuse were accounted for. This result, and the lack of association of the candidate genes for levels of alcohol consumption, with alcohol dependence/abuse in the same individuals, suggest that the genetic factors that we identified as predisposing factors for alcohol consumption may have little direct influence on alcohol dependence or abuse, but would be important for generating the major risk factor for dependence (i.e., high levels of alcohol consumption).
A recent genome-wide association study of alcohol dependence in humans used genes identified by differential mRNA expression in alcohol-consuming rats as a means of attempting to generate more credence in the candidate genes with modest statistical support for association with alcohol dependence in the analysis of the human data . Our studies would caution against such an approach because genetic determinants for alcohol consumption by animals or humans may not be identical with genetic determinants of ethanol dependence. This contention is also supported by QTL analysis of free-choice alcohol consumption  and alcohol dependence/withdrawal  in animals, and alcohol consumption/alcohol dependence in humans [113–115]. The most recent and instructive study in this regard is a study with humans performed by Hansell et al. . This study utilized quantity/frequency measures of alcohol consumption as a phenotype for a linkage study and found two types of QTLs. The LOD scores for certain QTLs were diminished if the primarily high level alcohol consumers were utilized for analysis, and other QTLs were enhanced if primarily high level consumers were used. What was interesting is that the QTLs which were enhanced when individuals with lower levels of alcohol consumption were dropped from the analysis were those which, in other studies, were related to "alcoholism" or alcohol dependence. QTLs identified as significant and having maximal LOD scores using the full range of alcohol consumption in this population, were not in areas of the human genome previously associated with alcohol dependence. However, at least one such QTL was previously found by de Andrade et al. to be associated with the number of drinks (on average) consumed by the subjects collected in the COGA study .
The genetical genomics approach, in combination with phenomics, is a powerful method for determination of candidate genes that contribute to the predisposition to alcohol consumption in rats. Informatics-based analysis of the function of candidate gene products led to the consideration that GABAergic function, and particularly GABA release modulated by peptidergic, purinergic and endogenous cannabinoid systems, as well as GABA receptor trafficking, are important components of the genetic/biochemical pathways that contribute to alcohol drinking. Our results highlight the need to identify neuronal signaling networks based on candidate genes, rather than focusing on individual genes and gene products, when attempting to understand the genetic basis of complex behavioral traits. The comparison of rat and human genetic contributors to the trait of alcohol consumption suggested that one can extrapolate from pathways - not necessarily specific genes - found in animals to begin to elucidate cross-species similarities in the genetic basis of behavior. Our results also emphasize the importance of carefully defining phenotypes for genetical genomic approaches; in this case, although high levels of alcohol drinking are phenotypically correlated with alcohol dependence, the genetic factors that contribute to the full range of alcohol consumption versus alcohol dependence in humans are distinct.
Drake TA, Schadt EE, Lusis AJ: Integrating genetic and gene expression data: application to cardiovascular and metabolic traits in mice. Mamm Genome. 2006, 17 (6): 466-479.
Tabakoff B, Saba L, Kechris K, Hu W, Bhave SV, Finn DA, Grahame NJ, Hoffman PL: The genomic determinants of alcohol preference in mice. Mamm Genome. 2008, 19 (5): 352-365.
Ganguly K, Stoeger T, Wesselkamper SC, Reinhard C, Sartor MA, Medvedovic M, Tomlinson CR, Bolle I, Mason JM, Leikauf GD, et al: Candidate genes controlling pulmonary function in mice: transcript profiling and predicted protein structure. Physiol Genomics. 2007, 31 (3): 410-421.
Tabakoff B, Bhave SV, Hoffman PL: Selective breeding, quantitative trait locus analysis, and gene arrays identify candidate genes for complex drug-related behaviors. J Neurosci. 2003, 23 (11): 4491-4498.
Hu W, Saba L, Kechris K, Bhave SV, Hoffman PL, Tabakoff B: Genomic insights into acute alcohol tolerance. J Pharmacol Exp Ther. 2008, 326 (3): 792-800.
Feltenstein MW, See RE: The neurocircuitry of addiction: an overview. Br J Pharmacol. 2008, 154 (2): 261-274.
Bell RL, Rodd ZA, Lumeng L, Murphy JM, McBride WJ: The alcohol-preferring P rat and animal models of excessive alcohol drinking. Addict Biol. 2006, 11 (3-4): 270-288.
Koob GF: Drugs of abuse: anatomy, pharmacology and function of reward pathways. Trends Pharmacol Sci. 1992, 13 (5): 177-184.
Kalivas PW, O'Brien C: Drug addiction as a pathology of staged neuroplasticity. Neuropsychopharmacology. 2008, 33 (1): 166-180.
Kelley AE: Memory and addiction: shared neural circuitry and molecular mechanisms. Neuron. 2004, 44 (1): 161-179.
Printz MP, Jirout M, Jaworski R, Alemayehu A, Kren V: Genetic Models in Applied Physiology. HXB/BXH rat recombinant inbred strain platform: a newly enhanced tool for cardiovascular, behavioral, and developmental genetics and genomics. J Appl Physiol. 2003, 94 (6): 2510-2522.
The STAR Consortium: SNP and haplotype mapping for genetic analysis in the rat. Nature Genetics. 2008, 40: 560-566.
Murphy JM, Stewart RB, Bell RL, Badia-Elder NE, Carr LG, McBride WJ, Lumeng L, Li TK: Phenotypic and genotypic characterization of the Indiana University rat lines selectively bred for high and low alcohol preference. Behav Genet. 2002, 32 (5): 363-388.
Hodgkinson CA, Yuan Q, Xu K, Shen PH, Heinz E, Lobos EA, Binder EB, Cubells J, Ehlers CL, Gelernter J, et al: Addictions biology: haplotype-based analysis for 130 candidate genes on a single array. Alcohol Alcohol. 2008, 43 (5): 505-515.
Glanz J, Grant B, Monteiro M, Tabakoff B: WHO/ISBRA Study on State and Trait Markers of Alcohol Use and Dependence: analysis of demographic, behavioral, physiologic, and drinking variables that contribute to dependence and seeking treatment. International Society on Biomedical Research on Alcoholism. Alcohol Clin Exp Res. 2002, 26 (7): 1047-1061.
Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping. Genetics. 1994, 138 (3): 963-971.
Broman KW, Wu H, Sen Ś, Churchill GA: R/qtl: QTL mapping in experimental crosses. Bioinformatics. 2003, 19 (7): 889-890.
Brown VM, Ossadtchi A, Khan AH, Yee S, Lacan G, Melega WP, Cherry SR, Leahy RM, Smith DJ: Multiplex three-dimensional brain gene expression mapping in a mouse model of Parkinson's disease. Genome Res. 2002, 12 (6): 868-884.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, et al: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-
Kent WJ: BLAT--the BLAST-like alignment tool. Genome Res. 2002, 12 (4): 656-664.
Carlborg O, De Koning DJ, Manly KF, Chesler E, Williams RW, Haley CS: Methodological aspects of the genetic dissection of gene expression. Bioinformatics. 2005, 21 (10): 2383-2393.
Visscher PM, Thompson R, Haley CS: Confidence intervals in QTL mapping by bootstrapping. Genetics. 1996, 143 (2): 1013-1020.
Bengtsson H, Simpson K, Bullard J, Hansen K: aroma.affymetrix: A generic framework in R for analyzing small to very large Affymetrix data sets in bounded memory. Tech Report #745. 2008, Berkeley: Department of Statistics, University of California
Affymetrix: Quality assessment of exon and gene arrays. Affymetrix GeneChip® Gene and Exon Array Whitepaper Collection. 2003
Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31 (4): e15-
Pritchard JK, Rosenberg NA: Use of unlinked genetic markers to detect population stratification in association studies. Am J Hum Genet. 1999, 65 (1): 220-228.
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265.
Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, et al: The structure of haplotype blocks in the human genome. Science. 2002, 296 (5576): 2225-2229.
The International HapMap Consortium: The International HapMap Project. Nature. 2003, 426: 789-796.
Schwarz G: Estimating the dimension of a model. Annals of Statistics. 1978, 6: 461-464.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. JR Statist Soc B. 1995, 57: 289-300.
Rhee SH, Hewitt JK, Young SE, Corley RP, Crowley TJ, Stallings MC: Genetic and environmental influences on substance initiation, use, and problem use in adolescents. Arch Gen Psychiatry. 2003, 60 (12): 1256-1264.
Kendler KS, Schmitt E, Aggen SH, Prescott CA: Genetic and environmental influences on alcohol, caffeine, cannabis, and nicotine use from early adolescence to middle adulthood. Arch Gen Psychiatry. 2008, 65 (6): 674-682.
Kuo PH, Aggen SH, Prescott CA, Kendler KS, Neale MC: Using a factor mixture modeling approach in alcohol dependence in a general population sample. Drug Alcohol Depend. 2008, 98 (1-2): 105-114.
Hasin D, Hatzenbuehler ML, Keyes K, Ogburn E: Substance use disorders: diagnostic and Statistical Manual of Mental Disorders, fourth edition (DSM-IV) and International Classification of Diseases, tenth edition (ICD-10). Addiction. 2006, 101 (Suppl 1): 59-75.
Moss HB, Chen CM, Yi HY: DSM-IV criteria endorsement patterns in alcohol dependence: relationship to severity. Alcohol Clin Exp Res. 2008, 32 (2): 306-313.
Koob GF: Neurobiological substrates for the dark side of compulsivity in addiction. Neuropharmacology. 2009, 56 (Suppl 1): 18-31.
Saha TD, Chou SP, Grant BF: Toward an alcohol use disorder continuum using item response theory: results from the National Epidemiologic Survey on Alcohol and Related Conditions. Psychol Med. 2006, 36 (7): 931-941.
Kozell LB, Hitzemann R, Buck KJ: Acute alcohol withdrawal is associated with c-Fos expression in the basal ganglia and associated circuitry: C57BL/6J and DBA/2J inbred mouse strain analyses. Alcohol Clin Exp Res. 2005, 29 (11): 1939-1948.
Yoneyama N, Crabbe JC, Ford MM, Murillo A, Finn DA: Voluntary ethanol consumption in 22 inbred mouse strains. Alcohol. 2008, 42 (3): 149-160.
Bachmanov AA, Kiefer SW, Molina JC, Tordoff MG, Duffy VB, Bartoshuk LM, Mennella JA: Chemosensory factors influencing alcohol perception, preferences, and consumption. Alcohol Clin Exp Res. 2003, 27 (2): 220-231.
Buemann B, Astrup A: How does the body deal with energy from alcohol?. Nutrition. 2001, 17 (7-8): 638-641.
Gao Q, Horvath TL: Neuronal control of energy homeostasis. FEBS Lett. 2008, 582 (1): 132-141.
Chen ZW, Chang CS, Leil TA, Olcese R, Olsen RW: GABAA receptor-associated protein regulates GABAA receptor cell-surface number in Xenopus laevis oocytes. Mol Pharmacol. 2005, 68 (1): 152-159.
Cowley MA, Pronchuk N, Fan W, Dinulescu DM, Colmers WF, Cone RD: Integration of NPY, AGRP, and melanocortin signals in the hypothalamic paraventricular nucleus: evidence of a cellular basis for the adipostat. Neuron. 1999, 24 (1): 155-163.
Barnes S, Whistler HL, Hughes J, Woodruff GN, Hunter JC: Effect of cholecystokinin octapeptide on endogenous amino acid release from the rat ventromedial nucleus of the hypothalamus and striatum. J Neurochem. 1991, 56 (4): 1409-1416.
Berridge KC: Food reward: brain substrates of wanting and liking. Neurosci Biobehav Rev. 1996, 20 (1): 1-25.
Kelley AE: Ventral striatal control of appetitive motivation: role in ingestive behavior and reward-related learning. Neurosci Biobehav Rev. 2004, 27 (8): 765-776.
Palmiter RD: Is dopamine a physiologically relevant mediator of feeding behavior?. Trends Neurosci. 2007, 30 (8): 375-381.
Di Marzo V, Ligresti A, Cristino L: The endocannabinoid system as a link between homoeostatic and hedonic pathways involved in energy balance regulation. Int J Obes (Lond). 2009, 33 (Suppl 2): S18-24.
Lovinger DM: Presynaptic modulation by endocannabinoids. Handb Exp Pharmacol. 2008, 435-477. 184
Viveros MP, de Fonseca FR, Bermudez-Silva FJ, McPartland JM: Critical role of the endocannabinoid system in the regulation of food intake and energy metabolism, with phylogenetic, developmental, and pathophysiological implications. Endocr Metab Immune Disord Drug Targets. 2008, 8 (3): 220-230.
Kola B, Farkas I, Christ-Crain M, Wittmann G, Lolli F, Amin F, Harvey-White J, Liposits Z, Kunos G, Grossman AB, et al: The orexigenic effect of ghrelin is mediated through central activation of the endogenous cannabinoid system. PLoS One. 2008, 3 (3): e1797-
Hentges ST, Low MJ, Williams JT: Differential regulation of synaptic inputs by constitutively released endocannabinoids and exogenous cannabinoids. J Neurosci. 2005, 25 (42): 9746-9751.
Matias I, Bisogno T, Di Marzo V: Endogenous cannabinoids in the brain and peripheral tissues: regulation of their levels and control of food intake. Int J Obes (Lond). 2006, 30 (Suppl 1): S7-S12.
Shi Y, Cheng D: Beyond triglyceride synthesis: the dynamic functional roles of MGAT and DGAT enzymes in energy metabolism. Am J Physiol Endocrinol Metab. 2009, 297 (1): E10-18.
Mercer LD, Le VQ, Nunan J, Jones NM, Beart PM: Direct visualization of cholecystokinin subtype2 receptors in rat central nervous system using anti-peptide antibodies. Neurosci Lett. 2000, 293 (3): 167-170.
Ghijsen WE, Leenders AG, Wiegant VM: Regulation of cholecystokinin release from central nerve terminals. Peptides. 2001, 22 (8): 1213-1221.
Miller KK, Hoffer A, Svoboda KR, Lupica CR: Cholecystokinin increases GABA release by inhibiting a resting K+ conductance in hippocampal interneurons. J Neurosci. 1997, 17 (13): 4994-5003.
Deng PY, Lei S: Bidirectional modulation of GABAergic transmission by cholecystokinin in hippocampal dentate gyrus granule cells of juvenile rats. J Physiol. 2006, 572 (Pt 2): 425-442.
Hugel S, Schlichter R: Presynaptic P2X receptors facilitate inhibitory GABAergic transmission between cultured rat spinal cord dorsal horn neurons. J Neurosci. 2000, 20 (6): 2121-2130.
Tiran Z, Peretz A, Attali B, Elson A: Phosphorylation-dependent regulation of Kv2.1 Channel activity at tyrosine 124 by Src and by protein-tyrosine phosphatase epsilon. J Biol Chem. 2003, 278 (19): 17509-17514.
Lindblom J, Opmane B, Mutulis F, Mutule I, Petrovska R, Klusa V, Bergstrom L, Wikberg JE: The MC4 receptor mediates alpha-MSH induced release of nucleus accumbens dopamine. Neuroreport. 2001, 12 (10): 2155-2158.
Navarro M, Cubero I, Chen AS, Chen HY, Knapp DJ, Breese GR, Marsh DJ, Thiele TE: Effects of melanocortin receptor activation and blockade on ethanol intake: a possible role for the melanocortin-4 receptor. Alcohol Clin Exp Res. 2005, 29 (6): 949-957.
Polidori C, Geary N, Massi M: Effect of the melanocortin receptor stimulation or inhibition on ethanol intake in alcohol-preferring rats. Peptides. 2006, 27 (1): 144-149.
Reum T, Schafer U, Marsden CA, Fink H, Morgenstern R: Cholecystokinin increases extracellular dopamine overflow in the anterior nucleus accumbens via CCK(B) receptors in the VTA assessed by in vivo voltammetry. Neuropeptides. 1997, 31 (1): 82-88.
Stuber GD, Hopf FW, Hahn J, Cho SL, Guillory A, Bonci A: Voluntary ethanol intake enhances excitatory synaptic strength in the ventral tegmental area. Alcohol Clin Exp Res. 2008, 32 (10): 1714-1720.
Davies DL, Machu TK, Guo Y, Alkana RL: Ethanol sensitivity in ATP-gated P2X receptors is subunit dependent. Alcohol Clin Exp Res. 2002, 26 (6): 773-778.
Xiao C, Ye JH: Ethanol dually modulates GABAergic synaptic transmission onto dopaminergic neurons in ventral tegmental area: role of mu-opioid receptors. Neuroscience. 2008, 153 (1): 240-248.
Ridley AJ: Rho GTPases and actin dynamics in membrane protrusions and vesicle trafficking. Trends Cell Biol. 2006, 16 (10): 522-529.
Kneussel M, Betz H: Clustering of inhibitory neurotransmitter receptors at developing postsynaptic sites: the membrane activation model. Trends Neurosci. 2000, 23 (9): 429-435.
Shen F, Lin Q, Gu Y, Childress C, Yang W: Activated Cdc42-associated kinase 1 is a component of EGF receptor signaling complex and regulates EGF receptor degradation. Mol Biol Cell. 2007, 18 (3): 732-742.
Urena JM, La Torre A, Martinez A, Lowenstein E, Franco N, Winsky-Sommerer R, Fontana X, Casaroli-Marano R, Ibanez-Sabio MA, Pascual M, et al: Expression, synaptic localization, and developmental regulation of Ack1/Pyk1, a cytoplasmic tyrosine kinase highly expressed in the developing and adult brain. J Comp Neurol. 2005, 490 (2): 119-132.
Striggow F, Riek-Burchardt M, Kiesel A, Schmidt W, Henrich-Noack P, Breder J, Krug M, Reymann KG, Reiser G: Four different types of protease-activated receptors are widely expressed in the brain and are up-regulated in hippocampus by severe ischemia. Eur J Neurosci. 2001, 14 (4): 595-608.
Luo W, Wang Y, Reiser G: Protease-activated receptors in the brain: receptor expression, activation, and functions in neurodegeneration and neuroprotection. Brain Res Rev. 2007, 56 (2): 331-345.
Wang P, Kumar P, Wang C, Defea KA: Differential regulation of class IA phosphoinositide 3-kinase catalytic subunits p110 alpha and beta by protease-activated receptor 2 and beta-arrestins. Biochem J. 2007, 408 (2): 221-230.
Jacob TC, Moss SJ, Jurd R: GABA(A) receptor trafficking and its role in the dynamic modulation of neuronal inhibition. Nat Rev Neurosci. 2008, 9 (5): 331-343.
Motil J, Dubey M, Chan WK, Shea TB: Inhibition of dynein but not kinesin induces aberrant focal accumulation of neurofilaments within axonal neurites. Brain Res. 2007, 1164: 125-131.
Raff EC, Hoyle HD, Popodi EM, Turner FR: Axoneme beta-tubulin sequence determines attachment of outer dynein arms. Curr Biol. 2008, 18 (12): 911-914.
Samson HH, Tolliver GA, Schwarz-Stevens K: Oral ethanol self-administration: a behavioral pharmacological approach to CNS control mechanisms. Alcohol. 1990, 7 (3): 187-191.
Hyytiä P, Koob GF: GABAA receptor antagonism in the extended amygdala decreases ethanol self-administration in rats. Eur J Pharmacol. 1995, 283 (1-3): 151-159.
Di Chiara G, Acquas E, Tanda G: Ethanol as a neurochemical surrogate of conventional reinforcers: the dopamine-opioid link. Alcohol. 1996, 13 (1): 13-17.
Hodge CW, Haraguchi M, Chappelle AM, Samson HH: Effects of ventral tegmental microinjections of the GABAA agonist muscimol on self-administration of ethanol and sucrose. Pharmacol Biochem Behav. 1996, 53 (4): 971-977.
Samson HH, Chappell A: Muscimol injected into the medial prefrontal cortex of the rat alters ethanol self-administration. Physiol Behav. 2001, 74 (4-5): 581-587.
Samson HH, Chappell A: Injected muscimol in pedunculopontine tegmental nucleus alters ethanol self-administration. Alcohol. 2001, 23 (1): 41-48.
Tabakoff B, Hoffman PL: Alcohol addiction: an enigma among us. Neuron. 1996, 16 (5): 909-912.
Roberto M, Madamba SG, Moore SD, Tallent MK, Siggins GR: Ethanol increases GABAergic transmission at both pre- and postsynaptic sites in rat central amygdala neurons. Proc Natl Acad Sci USA. 2003, 100 (4): 2053-2058.
Nie Z, Schweitzer P, Roberts AJ, Madamba SG, Moore SD, Siggins GR: Ethanol augments GABAergic transmission in the central amygdala via CRF1 receptors. Science. 2004, 303 (5663): 1512-1514.
Lobo IA, Harris RA: GABA(A) receptors and alcohol. Pharmacol Biochem Behav. 2008, 90 (1): 90-94.
Hodge CW, Chappelle AM, Samson HH: GABAergic transmission in the nucleus accumbens is involved in the termination of ethanol self-administration in rats. Alcohol Clin Exp Res. 1995, 19 (6): 1486-1493.
Hnasko TS, Perez FA, Scouras AD, Stoll EA, Gale SD, Luquet S, Phillips PE, Kremer EJ, Palmiter RD: Cre recombinase-mediated restoration of nigrostriatal dopamine in dopamine-deficient mice reverses hypophagia and bradykinesia. Proc Natl Acad Sci USA. 2006, 103 (23): 8858-8863.
Mountjoy KG, Mortrud MT, Low MJ, Simerly RB, Cone RD: Localization of the melanocortin-4 receptor (MC4-R) in neuroendocrine and autonomic control circuits in the brain. Mol Endocrinol. 1994, 8 (10): 1298-1308.
Kishi T, Aschkenasi CJ, Lee CE, Mountjoy KG, Saper CB, Elmquist JK: Expression of melanocortin 4 receptor mRNA in the central nervous system of the rat. J Comp Neurol. 2003, 457 (3): 213-235.
Balthasar N, Dalgaard LT, Lee CE, Yu J, Funahashi H, Williams T, Ferreira M, Tang V, McGovern RA, Kenny CD, et al: Divergence of melanocortin pathways in the control of food intake and energy expenditure. Cell. 2005, 123 (3): 493-505.
Chung L, Moore SD: Cholecystokinin enhances GABAergic inhibitory transmission in basolateral amygdala. Neuropeptides. 2007, 41 (6): 453-463.
Chhatwal JP, Gutman AR, Maguschak KA, Bowser ME, Yang Y, Davis M, Ressler KJ: Functional interactions between endocannabinoid and CCK neurotransmitter systems may be critical for extinction learning. Neuropsychopharmacology. 2009, 34 (2): 509-521.
Jhou TC, Fields HL, Baxter MG, Saper CB, Holland PC: The rostromedial tegmental nucleus (RMTg), a GABAergic afferent to midbrain dopamine neurons, encodes aversive stimuli and inhibits motor responses. Neuron. 2009, 61 (5): 786-800.
Hernández-Rabaza V, Hontecillas-Prieto L, Velazquez-Sanchez C, Ferragud A, Perez-Villaba A, Arcusa A, Barcia JA, Trejo JL, Canales JJ: The hippocampal dentate gyrus is essential for generating contextual memories of fear and drug-induced reward. Neurobiol Learn Mem. 2008, 90 (3): 553-559.
Balasubramanian S, Fam SR, Hall RA: GABAB receptor association with the PDZ scaffold Mupp1 alters receptor stability and function. J Biol Chem. 2007, 282 (6): 4162-4171.
Fehr C, Shirley RL, Belknap JK, Crabbe JC, Buck KJ: Congenic mapping of alcohol and pentobarbital withdrawal liability loci to a <1 centimorgan interval of murine chromosome 4: identification of Mpdz as a candidate gene. J Neurosci. 2002, 22 (9): 3730-3738.
Miller AD, Blaha CD: Midbrain muscarinic receptor mechanisms underlying regulation of mesoaccumbens and nigrostriatal dopaminergic transmission in the rat. Eur J Neurosci. 2005, 21 (7): 1837-1846.
Katner SN, McBride WJ, Lumeng L, Li TK, Murphy JM: Alcohol intake of P rats is regulated by muscarinic receptors in the pedunculopontine nucleus and VTA. Pharmacol Biochem Behav. 1997, 58 (2): 497-504.
Hsu R, Taylor JR, Newton SS, Alvaro JD, Haile C, Han G, Hruby VJ, Nestler EJ, Duman RS: Blockade of melanocortin transmission inhibits cocaine reward. Eur J Neurosci. 2005, 21 (8): 2233-2242.
Landgren S, Jerlhag E, Zetterberg H, Gonzalez-Quintela A, Campos J, Olofsson U, Nilsson S, Blennow K, Engel JA: Association of pro-ghrelin and GHS-R1A gene polymorphisms and haplotypes with heavy alcohol use and body mass. Alcohol Clin Exp Res. 2008, 32 (12): 2054-2061.
Hansell NK, Agrawal A, Whitfield JB, Morley KI, Gordon SD, Lind PA, Pergadia ML, Montgomery GW, Madden PA, Todd RD, et al: Can we identify genes for alcohol consumption in samples ascertained for heterogeneous purposes?. Alcohol Clin Exp Res. 2009, 33 (4): 729-739.
Belknap JK, Belknap ND, Berg JH, Coleman R: Preabsorptive vs. postabsorptive control of ethanol intake in C57BL/6J and DBA/2J mice. Behav Genet. 1977, 7 (6): 413-425.
Gentry RT, Dole VP: Why does a sucrose choice reduce the consumption of alcohol in C57BL/6J mice?. Life Sci. 1987, 40 (22): 2191-2194.
Jerlhag E, Egecioglu E, Landgren S, Salome N, Heilig M, Moechars D, Datta R, Perrissoud D, Dickson SL, Engel JA: Requirement of central ghrelin signaling for alcohol reward. Proc Natl Acad Sci USA. 2009, 106 (27): 11318-11323.
de Andrade M, Olswold CL, Slusser JP, Tordsen LA, Atkinson EJ, Rabe KG, Slager SL: Identification of genes involved in alcohol consumption and cigarettes smoking. BMC Genet. 2005, 6 (Suppl 1): S112-
Treutlein J, Cichon S, Ridinger M, Wodarz N, Soyka M, Zill P, Maier W, Moessner R, Gaebel W, Dahmen N, et al: Genome-wide association study of alcohol dependence. Arch Gen Psychiatry. 2009, 66 (7): 773-784.
Bice P, Valdar W, Zhang L, Liu L, Lai D, Grahame N, Flint J, Li TK, Lumeng L, Foroud T: Genomewide SNP screen to detect quantitative trait loci for alcohol preference in the high alcohol preferring and low alcohol preferring mice. Alcohol Clin Exp Res. 2009, 33 (3): 531-537.
Bergeson SE, Kyle Warren R, Crabbe JC, Metten P, Gene Erwin V, Belknap JK: Chromosomal loci influencing chronic alcohol withdrawal severity. Mamm Genome. 2003, 14 (7): 454-463.
Agrawal A, Hinrichs AL, Dunn G, Bertelsen S, Dick DM, Saccone SF, Saccone NL, Grucza RA, Wang JC, Cloninger CR, et al: Linkage scan for quantitative traits identifies new regions of interest for substance dependence in the Collaborative Study on the Genetics of Alcoholism (COGA) sample. Drug Alcohol Depend. 2008, 93 (1-2): 12-20.
Guerrini I, Cook CC, Kest W, Devitgh A, McQuillin A, Curtis D, Gurling HM: Genetic linkage analysis supports the presence of two susceptibility loci for alcoholism and heavy drinking on chromosome 1p22.1-11.2 and 1q21.3-24.2. BMC Genet. 2005, 6 (1): 11-
Bergen AW, Yang XR, Bai Y, Beerman MB, Goldstein AM, Goldin LR: Genomic regions linked to alcohol consumption in the Framingham Heart Study. BMC Genet. 2003, 4 (Suppl 1): S101-
Wan S, Browning KN, Coleman FH, Sutton G, Zheng H, Butler A, Berthoud HR, Travagli RA: Presynaptic melanocortin-4 receptors on vagal afferent fibers modulate the excitability of rat nucleus tractus solitarius neurons. J Neurosci. 2008, 28 (19): 4957-4966.
Kawaguchi Y, Kondo S: Parvalbumin, somatostatin and cholecystokinin as chemical markers for specific GABAergic interneuron types in the rat frontal cortex. J Neurocytol. 2002, 31 (3-5): 277-287.
Xiao C, Zhou C, Li K, Davies DL, Ye JH: Purinergic type 2 receptors at GABAergic synapses on ventral tegmental area dopamine neurons are targets for ethanol action. J Pharmacol Exp Ther. 2008, 327 (1): 196-205.
The authors thank Janet Hopkins for performance of assays for measuring mRNA levels; Ron Smith, Laura Breen and Joseph Gatewood for technical assistance with behavioral data collection and Dr. Samuel P. Hazen for expert advice on database management. WHO/ISBRA Study on State and Trait Markers of Alcoholism is an abbreviation for World Health Organization/International Society for Biomedical Research on Alcoholism Study on State and Trait Markers of Alcohol Use and Dependence Investigators. This work was supported in part by NIAAA, NIH (AA013162 (BT); AA006420 (GK); AA016649-INIA Project (PH); AA016663-INIA Project (BT); AA013522-INIA Project (RLB); AA013517-INIA Project (GK); AA-U01-Developmental Grant - INIA Project (M Printz); AA16922 (KK)); NHLBI, NIH (HL35018 (M Printz)); Banbury Fund (BT); Pearson Center for Alcoholism and Addiction Research (GK); National Health and Medical Research Council of Australia (WHO/ISBRA Sydney); La Fondation Jean Lapointe (WHO/ISBRA Montreal); National Health Research Development Program, Canada (WHO/ISBRA Montreal); Wellcome Trust Program for Cardiovascular Functional Genomics (JM); Ministry of Education of the Czech Republic (1M6837805002(M Pravenec)) and Howard Hughes Medical Institute (M Pravenec).
BT planned and supervised all experiments, analyzed data and drafted the manuscript. LS carried out eQTL and bQTL analyses, correlation and all statistical analyses. MPrintz provided data and tissue from the HXB/BXH rats, participated in alcohol consumption studies, and analyzed alcohol consumption data. PF participated in HXB/BXH rat alcohol consumption experiments and analyzed data. CH performed studies using the addiction array. DG supervised addiction array studies and analyzed data. GK conceived, organized and supervised the HXB/BXH alcohol consumption studies. HNR organized and supervised the HXB/BXH rat alcohol consumption studies and analyzed data. KK performed transcription factor and transcription factor binding site analysis. RLB provided data and tissue from P and NP rats. NH performed genotyping of HXB/BXH RI strains. MH performed genotyping of HXB/BXH RI strains. MPravenec provided tissue from HXB/BXH rats. JM provided genotype data from HXB/BXH RI rats. LL recruited subjects from Montreal and provided data and DNA. MD supervised studies and recruited subjects in Montreal. KMC recruited subjects in Australia and provided data and DNA. JW supervised human studies and recruited subjects in Australia. JS recruited subjects in Australia and analyzed data. BG created questionnaire for WHO/ISBRA study and generated algorithms for data analysis. PLH analyzed microarray data and drafted manuscript. WHO/ISBRA Investigators are included as authors on all studies using data from the WHO/ISBRA study, by formal agreement.