# Genetic Contribution to Initial and Progressive Alcohol Intake Among Recombinant Inbred Strains of Mice

Megan K. Mulligan, Wenyuan Zhao, Morgan Dickerson, Danny Arends, Pjotr

Prins, Sonia A. Cavigelli, Elena Terenina, Pierre Mormède, Lu Lu, Byron C.

Jones

► To cite this version:

Megan K. Mulligan, Wenyuan Zhao, Morgan Dickerson, Danny Arends, Pjotr Prins, et al.. Genetic Contribution to Initial and Progressive Alcohol Intake Among Recombinant Inbred Strains of Mice.

Frontiers in Genetics, 2018, 9, Non paginé. 10.3389/fgene.2018.00370. hal-02625991

HAL Id: hal-02625991

https://hal.inrae.fr/hal-02625991

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Genetic Contribution to Initial and Progressive Alcohol Intake Among Recombinant Inbred Strains of Mice

Megan K. Mulligan¹*, Wenyuan Zhao¹, Morgan Dickerson¹, Danny Arends², Pjotr Prins³, Sonia A. Cavigelli⁴, Elena Terenina⁵, Pierre Mormede⁵, Lu Lu¹ and Byron C. Jones¹*

¹Department of Genetics, Genomics, and Informatics, The University of Tennessee Health Science Center, Memphis, TN, United States, ²Albrecht Daniel Thaer-Institut für Agrar- und Gartenbauwissenschaften, Humboldt-Universität zu Berlin, Berlin, Germany, ³Biomedical Genetics, University Medical Center Utrecht, Utrecht, Netherlands, ⁴Department of BioBehavioral Health, The Pennsylvania State University, University Park, PA, United States, ⁵GenPhySE, INRA, ENVT, Université de Toulouse, Castanet-Tolosan, France

We profiled individual differences in alcohol consumption upon initial exposure and during 5 weeks of voluntary alcohol intake in female mice from 39 BXD recombinant inbred strains and parents using the drinking in the dark (DID) method. In this paradigm, a single bottle of 20% (v/v) alcohol was presented as the sole liquid source for 2 or 4 h starting 3 h into the dark cycle. For 3 consecutive days mice had access to alcohol for 2 h followed by a 4th day of 4 h access and 3 intervening days where alcohol was not offered. We followed this regime for 5 weeks. For most strains, 2 or 4 h alcohol intake increased over the 5-week period, with some strains demonstrating greatly increased intake. There was considerable and heritable genetic variation in alcohol consumption upon initial early and sustained weekly exposure. Two different mapping algorithms were used to identify QTLs associated with alcohol intake and only QTLs detected by both methods were considered further. Multiple suggestive QTLs for alcohol intake on chromosomes (Chrs) 2, 6, and 12 were identified for the first 4 h exposure. Suggestive QTLs for sustained intake during later weeks were identified on Chrs 4 and 8. Thirty high priority candidate genes, including *Entpd2*, *Per3*, and *Fto* were nominated for early and sustained alcohol intake QTLs. In addition, a suggestive QTL on Chr 15 was detected for change in 2 h alcohol intake over the duration of the study and *Adcy8* was identified as a strong candidate gene. Bioinformatic analyses revealed that early and sustained alcohol intake is likely driven by genes and pathways involved in signaling, and/or immune and metabolic function, while a combination of epigenetic factors related to alcohol experience and genetic factors likely drives progressive alcohol intake.

**Keywords:** BXD, DID, QTL, alcoholism, alcohol intake, genetic variation, B6, D2

## INTRODUCTION

According to Cloninger (1987) and Babor et al. (1992) there are multiple types of alcohol use disorders and likely different genetic contributions to each type. For each, though there is an initiating event, usually voluntary consumption and subsequent developments from avoidance to steady-state or ever increasing consumption and associated problems.

Genetic populations of animals have been used to model individual differences in the propensity to consume alcohol. One murine population that has contributed greatly to

**OPEN ACCESS**

**Edited by:**

Feng C. Zhou,
Indiana University Bloomington,
United States

**Reviewed by:**

Mark Z. Kos,
University of Texas Rio Grande Valley
Edinburg, United States
Steven Hicks,
Penn State Hershey Children's
Hospital, United States

*Correspondence:

Megan K. Mulligan
mkmulligan@uthsc.edu;
mmulliga@uthsc.edu
Byron C. Jones
bjone129@uthsc.edu

**Specialty section:**

This article was submitted to
*Behavioral and Psychiatric Genetics*,
a section of the journal
*Frontiers in Genetics*

**Received:** 25 May 2018

**Accepted:** 22 August 2018

**Published:** 25 September 2018

**Citation:**

Mulligan MK, Zhao W, Dickerson M, Arends D, Prins P, Cavigelli SA, Terenina E, Mormede P, Lu L and Jones BC (2018) Genetic Contribution to Initial and Progressive Alcohol Intake Among Recombinant Inbred Strains of Mice. *Front. Genet.* 9:370.
doi: 10.3389/fgene.2018.00370
alcohol and addiction research is the BXD recombinant inbred family derived from C57BL/6J (B6) and DBA/2J (D2) inbred strains. Mice from these parental strains were crossbred to produce an F₁ generation and then this generation was interbred for several generations to produce genetically segregating stocks. Next, families were selected for inbreeding to fix the alleles. The approach is described in more detail by Peirce et al. (2004). The result is a large number of inbred strains in which the alleles from B6 and D2 were recombined and redistributed. The BXD strains are segregating over five million variants that distinguish the parental B6 and D2 strains and all of the strains have been densely genotyped. Over 10 brain regions have been subjected to microarray analysis of gene expression and there is a freely available database<fnref n="1" /> consisting of more than 5,000 phenotypes contributed by many laboratories. This includes over 15 alcohol-related data sets in which the BXD family has been used to measure alcohol acceptance, consumption and preference (Phillips et al., 1994; Rodriguez et al., 1994; Gill et al., 1996; Fernandez et al., 1999); *metabolism* (Browman and Crabbe, 2000; Grisel et al., 2002; Philip et al., 2010); *hypothermia*, withdrawal, tolerance, and sensitivity (Belknap et al., 1993; Roberts et al., 1995; Crabbe et al., 1996; Phillips et al., 1996; Buck et al., 1997; Crabbe, 1998; Browman and Crabbe, 2000; Philip et al., 2010); *locomotor response* (Phillips et al., 1995; Browman and Crabbe, 2000; Philip et al., 2010); ethanol induced conditioned taste aversion (Risinger and Cunningham, 1998); and ethanol conditioned place preference (Cunningham, 1995). This family also contributed to the most detailed meta-analysis of genes that contribute to the predisposition for high alcohol consumption (Mulligan et al., 2006).

The research that we present here is a continuation of these studies in which we report on initial early intake and sustained alcohol consumption over a 5-week period in 39 BXD strains using the “drinking in the dark” (DID) protocol described by Rhodes et al. (2005). This protocol is used to elicit high alcohol consumption over a short time span and has been used to model binge-like alcohol consumption. There was considerable genetic variation among inbred strains of mice in alcohol intake using the DID paradigm (Rhodes et al., 2007; Crabbe et al., 2014). However, the DID procedure has typically been used to measure intake during a single week and has never been used to measure alcohol intake over multiple weeks in a recombinant inbred population. Here, we leverage the BXD family and accompanying legacy molecular, alcohol-trait, and other existing phenotypes to provide a systems genetics analysis of the factors driving alcohol consumption in the BXD family.

# MATERIALS AND METHODS

## Animals

The subjects were female mice from 39 BXD recombinant inbred mouse strains and the two parental strains for the BXDs, B6 and D2 (within strain replicates ranged from 1 to 16, **Supplementary Table S1**). The animals were 60–80 days old at

the start of the study. The mice were individually housed and fed a standard laboratory diet (Harlan Teklad 7912) with food and water available *ad libitum* except during exposure to alcohol (see below). The light cycle was 23:00 h lights on and 11:00 h lights off. This light cycle facilitated alcohol administration and measurement of alcohol intake. Mice were weighed weekly, and all procedures included here were approved by the UTHSC Institutional Animal Care and Use Committee.

## Drinking in the Dark

Alcohol consumption was evaluated using the DID method (Rhodes et al., 2005). The protocol calls for 4 consecutive days of testing. Each day, starting on a Tuesday 3 h after lights were turned off, the water bottles were removed from the cages and replaced with 15 ml centrifuge tubes filled with 20% (v/v) ethanol from 95% USP ethanol. On days 1–3 (Tuesday through Thursday) the length of exposure was 2 h and on the 4th day (Friday) the exposure was 4 h. No alcohol was offered in the intervening period (Saturday through Monday). This protocol was repeated weekly for 5 weeks. Tubes were weighed immediately before and after the exposure period. Volume consumed was converted to g/kg body weight of ethanol. Subsequent weekly 2 or 4 h measurements were averaged by strain and used as the dependent variable for data analysis and QTL mapping. The data were deposited in Gene Network (GN)<fnref n="1" /> and are available as traits 20010 through 20014 and 20077 through 20082 in the BXD Published Phenotypes database.

## Data Analysis

Statistical evaluation of daily alcohol intake was performed by ANOVA testing in R using the *lm* function for a 1 between-subjects variable (strain). Heritability at each time point was estimated from the ANOVA results by $ss_{strain}/ss_{total}$ (where ss = sum of squares, Belknap, 1998). Average weekly 2 h or 4 h intake was used to compute slopes and intercepts for each strain over the 5-week observation period.

Genetic correlational analyses among the DID phenotypes, between the DID phenotypes and BXD legacy phenotype data, and traditional quantitative trait loci (QTL) analysis were performed using a combination of R and GN software (Sloan et al., 2016)<fnref n="1" />. Traditional interval mapping in GN was performed using a simple regression method (Haley-Knott or HK) to compute QTL probability given strain genotypes and alcohol intake averaged by strain (Chesler et al., 2005; Mulligan et al., 2017). For traditional interval mapping, genome-wide suggestive (adjusted *p* < 0.63) and significant (adjusted *p* < 0.05) thresholds were determined based on 1,000 permutations of the trait data for each phenotype (GN default). The suggestive threshold is very permissive (see GN glossary of terms and features at www.genenetwork.org/glossary.html for more details) but strikes a balance between detection of false positives and highlighting loci that might be worth further investigation.

QTL analysis using genome-wide efficient mixed model association (GEMMA; Zhou and Stephens, 2012) was performed as a secondary method to traditional HK mapping in GN version 2 to assess the effect of population structure (kinship) between individuals. Full GEMMA Linear Mixed Model (LMM)

<footnote n="1">¹www.GeneNetwork.org</footnote>
support with the optional leave one chromosome out (LOCO) method was recently added to GN. GEMMA software is a computationally efficient LMM method for QTL mapping while explicitly accounting for genetic non-independence within each sample (Zhou and Stephens, 2012). Even for small sample sizes, after running Haley-Knott QTL mapping, GEMMA potentially allows for explorative fine-tuning of results at the SNP level. At each time point, and for each of 7,320 SNPs and phenotypes, we fitted the LMM using a kinship matrix K computed over the SNP genotypes. We also estimated significance thresholds with GEMMA using a permutation approach computing results 1,000 times while shuffling the phenotypes but keeping the genotypes and K the same (Churchill and Doerge, 1994). From every permutation we stored the highest Wald-test *p*-value in an ordered set and set the significance threshold at the 95th percentile and the suggestive threshold at the 67th percentile. For our data set, the average GEMMA significance threshold was LOD 4.1 and the average suggestive threshold was at LOD 3.3.

GEMMA was included as a mapping method for two main reasons. The first reason is that additional power can be gained when using mixed-model association methods. The increase in power results both from accounting for phenotypic covariance due to genetic similarity and by conditioning on associated markers as opposed to a single candidate locus (Yang et al., 2014). The second reason is that the impact of background structure on QTL detection has not been rigorously evaluated in BXD data sets. Thus, GEMMA can be applied in addition to traditional HK mapping to minimize detection of false positive QTLs resulting from kinship. QTLs identified by both approaches are unlikely to result from genetic similarity. Using both approaches has advantages over using a single mapping model because genetic relatedness is not addressed by the HK model and our data set is small (e.g., underpowered) for GEMMA.

For both HK and GEMMA QTL mapping methods, a 1.5 LOD drop from the top marker was used to define QTL confidence interval regions. Loci detected using traditional HK mapping (suggestive threshold or above) that were also detected using the secondary GEMMA method (LOD > 3) were considered for further analysis.

Enrichment analysis was performed using tools available at Enrichr (Kuleshov et al., 2016). Default settings were used for Enrichr.

B6 and D2 polymorphic genes with SNPs and/or small insertions/deletions (InDels) were identified using tools available at the Sanger Institute Mouse Genomes Project<fnref n="2" /> (Keane et al., 2011). The Sanger website provides annotations for variants and, for our analysis, a SNP or InDel was considered to be of potentially high impact if it was annotated as a “coding sequence variant,” “feature elongation,” “feature truncation,” “incomplete terminal codon variant,” “initiator codon variant,” “mature miRNA variant,” “missense variant,” “NMD transcript variant,” “regulatory region ablation,” “regulatory region amplification,” “regulatory region variant,” “splice acceptor variant,” “splice donor variant,” “splice region variant,” “stop gained,” “stop lost,”

“TF binding site variant,” “TFBS ablation,” “TFBS amplification,” “transcript ablation,” or “transcript amplification.”

A large BXD database of hippocampal gene expression profiles generated from 67 naïve BXD strains [Hippocampus Consortium M430v2 (Jun 06); GN110; (Overall et al., 2009)] was used to prioritize candidate genes based on correlation with DID week 1 and 4 h phenotype data and cis expression QTL (eQTL) mapping.

For all candidate genes, literature associations between the term “alcohol” and each candidate gene were mined using the Chilibot<fnref n="3" /> website (Chen and Sharp, 2004).

## RESULTS

### Initial Alcohol Consumption Is Variable and Heritable in the BXD Panel and Increases Over Time

Intake of 20% alcohol for 2 or 4 h during the first and 5th week was variable in female B6, D2, and 39 BXD strains (Figure 1). Alcohol intake during the first 2 h exposure ranged from 1 to 3.55 g/kg compared to the last 2 h (week 5 day 3 or W5D3) exposure, which ranged from 1.67 to 5.01 g/kg. On the first 4 h exposure, alcohol intake ranged from 2.42 to 6.33 g/kg compared to intake on the last exposure, which ranged from 2.77 to 7.41 g/kg. ANOVA revealed a significant strain effect (all *p* < 0.0001) on intake for each time point with associated heritability estimates of 0.3 or greater, except for the first exposure on week 1 day 1 (W1D1) and W4D4 (Table 1). As expected, heritability on the first exposure to alcohol is lower than on subsequent exposure. Also as expected, heritability of 2 h weekly intake averaged for each strain over 3 days is much higher than strain-averaged 4 h intake resulting from a single measurement.

Average 2 or 4 h intake across all strains and the overall distribution of intake changed between week 1 and later weeks reflecting changes in drinking patterns that resulted from increased consumption over time (Figure 1). For 2 h intake, week was significantly associated (*p* < 0.0001; $R^2$ = 0.039) with a 0.14 g/kg increase in alcohol intake per week (~0.6 g/kg average increase in intake between weeks 1 and 5). For 4 h intake, week was significantly associated (*p* < 0.0001; $R^2$ = 0.025) with a 0.2 g/kg increase in alcohol intake per week (~0.8 g/kg average increase in intake between weeks 1 and 5). Pairwise correlations between initial daily, 2 h weekly average, and 4 h weekly alcohol intake tended to be higher between adjacent weeks and degrade as the interval between each week increases (Supplementary Figure S1).

### Genetic Differences in Alcohol Intake Over Time

Most strains demonstrated a gradual increase in alcohol intake over time (Table 2) during 2 or 4 h exposure. Strains with significantly increased 2 h intake over the 5-week period included BXD 40, 24, 49, 34, 50, 77, and 83. Of these strains, BXD 34, 50, 77, and 83 also demonstrated significantly elevated 4 h intake. Of

<footnote n="2">²http://www.sanger.ac.uk/sanger/Mouse_SnpViewer/</footnote>

<footnote n="3">³www.chilibot.net</footnote>
<figure><img src="image_3.png" /></figure>

**FIGURE 1** | Strain distribution of 2 or 4 h alcohol intake over 5 weeks. (A) Alcohol intake of 20% (v/v) over a 2 h access period is shown in g/kg on the Y-axis by week and the X-axis lists the BXD strain numbers with B6 and D2 parental strains shown to the far right (bar graph, each color represents an inbred strain). Population mean and standard deviation are shown for each week by the solid and dashed horizontal line, respectively. Weekly 2 h intake was variable among BXD strains and increased over time. (B) Distribution (histogram) reflects week-averaged intake from individual B6, D2, and 39 BXD strains recorded over three consecutive days of 2 h exposure. The distribution of alcohol intake changed between weeks 1 and 5 indicating increased alcohol intake. (C) Alcohol intake of 20% (v/v) over a 4 h access period is shown in g/kg on the Y-axis by week and the X-axis lists the BXD strain numbers with B6 and D2 parental strains shown to the far right (bar graph). Alcohol intake of 20% alcohol during the 4 h access period (bar graph) was variable among BXD strains and increased over time. (D) Distribution (histogram) reflects intake from B6, D2, and 39 BXD strains recorded over a single 4 h exposure. The distribution of alcohol intake changed between weeks 1 and 5 indicating increased alcohol intake. For the 2 and 4 h exposure, alcohol intake varied by strain both within and across weeks. For both exposures, mean intake increased between weeks 1 and 5.

<table><caption>**TABLE 1** | Anova results table.</caption><thead><tr><th rowspan="2">Strain Effect</th><th colspan="6">2 h access</th><th colspan="5">4 h access</th></tr><tr><th>W1D1</th><th>W1avg</th><th>W2avg</th><th>W3avg</th><th>W4avg</th><th>W5avg</th><th>W1D4</th><th>W2D4</th><th>W3D4</th><th>W4D4</th><th>W5D4</th></tr></thead><tbody><tr><td>R-value</td><td>1.20E-09</td><td>&lt;2.20E-16</td><td>&lt;2.20E-16</td><td>&lt;2.20E-16</td><td>&lt;2.20E-16</td><td>&lt;2.20E-16</td><td>1.22E-11</td><td>1.69E-13</td><td>2.89E-15</td><td>3.80E-08</td><td>2.06E-12</td></tr><tr><td>F-value</td><td>3.36</td><td>5.19</td><td>4.91</td><td>7.89</td><td>5.48</td><td>6.69</td><td>3.79</td><td>4.19</td><td>4.6</td><td>3.0</td><td>3.962</td></tr><tr><td>Df</td><td>39</td><td>40</td><td>40</td><td>40</td><td>40</td><td>40</td><td>40</td><td>40</td><td>40</td><td>40</td><td>40</td></tr><tr><td>h2</td><td>0.29</td><td>0.39</td><td>0.37</td><td>0.49</td><td>0.40</td><td>0.45</td><td>0.32</td><td>0.34</td><td>0.36</td><td>0.27</td><td>0.33</td></tr></tbody></table>

note, some strains with high genetic similarity (~80%, denoted by letters following strain name), such as BXD48 and BXD48a, showed highly divergent intake patterns over time with BXD48

showing significantly increased intake and BXD48a displaying stable or slightly decreased intake over time. The parental B6 strain demonstrated increased 2 h intake over the 5-week period
<table><caption>**TABLE 2** | Change in 2 and 4 h alcohol intake over 5-week.</caption><thead><tr><th rowspan="2">Strain</th><th colspan="3">2 h</th><th colspan="3">4 h</th></tr><tr><th>Slope</th><th>Intercept</th><th>SE</th><th>Slope</th><th>Intercept</th><th>SE</th></tr></thead><tbody><tr><td>BXD40</td><td><strong>0.329</strong></td><td>1.54</td><td>0.07</td><td>0.541</td><td>3.05</td><td>0.54</td></tr><tr><td>BXD24</td><td><strong>0.297</strong></td><td>1.51</td><td>0.07</td><td>0.545</td><td>3.20</td><td>0.55</td></tr><tr><td>C57BL/6J</td><td>0.278</td><td>3.20</td><td>0.11</td><td>-0.077</td><td>6.60</td><td>-0.08</td></tr><tr><td>BXD49</td><td><strong>0.246</strong></td><td>2.48</td><td>0.03</td><td>-0.455</td><td>6.01</td><td>-0.46</td></tr><tr><td>BXD34</td><td><strong>0.238</strong></td><td>1.16</td><td>0.04</td><td><strong>0.408</strong></td><td>2.46</td><td>0.41</td></tr><tr><td>BXD68</td><td>0.218</td><td>2.09</td><td>0.07</td><td>-0.002</td><td>4.96</td><td>0.00</td></tr><tr><td>BXD50</td><td><strong>0.211</strong></td><td>1.18</td><td>0.01</td><td><strong>0.232</strong></td><td>2.83</td><td>0.23</td></tr><tr><td>BXD77</td><td><strong>0.199</strong></td><td>1.22</td><td>0.02</td><td><strong>0.374</strong></td><td>2.47</td><td>0.37</td></tr><tr><td>BXD83</td><td><strong>0.197</strong></td><td>1.50</td><td>0.03</td><td><strong>0.408</strong></td><td>3.02</td><td>0.41</td></tr><tr><td>BXD62</td><td><strong>0.183</strong></td><td>1.57</td><td>0.03</td><td>0.153</td><td>3.59</td><td>0.15</td></tr><tr><td>BXD90</td><td><strong>0.169</strong></td><td>1.93</td><td>0.03</td><td>0.188</td><td>4.05</td><td>0.19</td></tr><tr><td>BXD73a</td><td><strong>0.168</strong></td><td>1.95</td><td>0.01</td><td><strong>0.362</strong></td><td>3.80</td><td>0.36</td></tr><tr><td>BXD75</td><td><strong>0.166</strong></td><td>1.21</td><td>0.04</td><td>0.198</td><td>2.62</td><td>0.20</td></tr><tr><td>BXD100</td><td><strong>0.161</strong></td><td>1.57</td><td>0.04</td><td>0.139</td><td>3.54</td><td>0.14</td></tr><tr><td>BXD70</td><td>0.160</td><td>2.23</td><td>0.06</td><td>0.132</td><td>4.88</td><td>0.13</td></tr><tr><td>BXD98</td><td>0.150</td><td>4.28</td><td>0.30</td><td>-0.318</td><td>6.87</td><td>-0.32</td></tr><tr><td>BXD73b</td><td><strong>0.146</strong></td><td>1.58</td><td>0.01</td><td><strong>0.319</strong></td><td>3.04</td><td>0.32</td></tr><tr><td>BXD78</td><td><strong>0.146</strong></td><td>2.61</td><td>0.03</td><td><strong>0.195</strong></td><td>5.72</td><td>0.19</td></tr><tr><td>BXD71</td><td><strong>0.144</strong></td><td>1.91</td><td>0.03</td><td>0.272</td><td>4.37</td><td>0.27</td></tr><tr><td>BXD60</td><td><strong>0.134</strong></td><td>1.75</td><td>0.02</td><td><strong>0.150</strong></td><td>3.81</td><td>0.15</td></tr><tr><td>BXD65b</td><td><strong>0.134</strong></td><td>1.96</td><td>0.04</td><td>-0.109</td><td>4.42</td><td>-0.11</td></tr><tr><td>BXD48</td><td><strong>0.134</strong></td><td>2.76</td><td>0.02</td><td><strong>0.542</strong></td><td>5.02</td><td>0.54</td></tr><tr><td>BXD44</td><td>0.127</td><td>1.74</td><td>0.06</td><td>0.296</td><td>2.91</td><td>0.30</td></tr><tr><td>BXD101</td><td>0.125</td><td>2.19</td><td>0.06</td><td>0.355</td><td>3.85</td><td>0.35</td></tr><tr><td>BXD87</td><td>0.120</td><td>1.58</td><td>0.04</td><td><strong>0.219</strong></td><td>3.23</td><td>0.22</td></tr><tr><td>BXD73</td><td>0.118</td><td>1.94</td><td>0.04</td><td>0.060</td><td>4.28</td><td>0.06</td></tr><tr><td>BXD65</td><td><strong>0.116</strong></td><td>1.73</td><td>0.03</td><td>0.232</td><td>3.36</td><td>0.23</td></tr><tr><td>BXD61</td><td>0.109</td><td>1.21</td><td>0.08</td><td>-0.372</td><td>3.90</td><td>-0.37</td></tr><tr><td>BXD43</td><td><strong>0.107</strong></td><td>1.83</td><td>0.03</td><td><strong>0.288</strong></td><td>3.39</td><td>0.29</td></tr><tr><td>BXD102</td><td>0.103</td><td>1.14</td><td>0.19</td><td>0.183</td><td>1.84</td><td>0.18</td></tr><tr><td>BXD86</td><td><strong>0.103</strong></td><td>2.09</td><td>0.03</td><td>0.218</td><td>4.13</td><td>0.22</td></tr><tr><td>BXD69</td><td><strong>0.101</strong></td><td>2.03</td><td>0.03</td><td>0.197</td><td>3.58</td><td>0.20</td></tr><tr><td>BXD55</td><td>0.088</td><td>1.79</td><td>0.04</td><td>0.122</td><td>3.17</td><td>0.12</td></tr><tr><td>BXD51</td><td>0.080</td><td>2.36</td><td>0.03</td><td>0.379</td><td>4.08</td><td>0.38</td></tr><tr><td>BXD66</td><td>0.047</td><td>2.47</td><td>0.03</td><td>0.085</td><td>4.33</td><td>0.09</td></tr><tr><td>BXD32</td><td>0.018</td><td>1.61</td><td>0.01</td><td>0.075</td><td>2.73</td><td>0.07</td></tr><tr><td>BXD4S</td><td>0.007</td><td>2.56</td><td>0.21</td><td>-0.356</td><td>5.57</td><td>-0.36</td></tr><tr><td>BXD79</td><td>-0.017</td><td>2.61</td><td>0.03</td><td>-0.080</td><td>5.25</td><td>-0.08</td></tr><tr><td>BXD63</td><td>-0.026</td><td>2.25</td><td>0.05</td><td>-0.125</td><td>4.52</td><td>-0.13</td></tr><tr><td>DBA/2J</td><td>-0.041</td><td>2.17</td><td>0.09</td><td>-0.040</td><td>4.32</td><td>-0.04</td></tr><tr><td>BXD48a</td><td>-0.046</td><td>2.89</td><td>0.05</td><td>-0.264</td><td>5.72</td><td>-0.26</td></tr></tbody></table>

Slope (g/kg per week) and associated y-intercept (g/kg) and standard error (SE) are shown for the 2 or 4 h exposure period by strain. Slopes sorted by 2 h exposure period. All slopes calculated from strain-averaged intake (4 h g/kg) or average weekly intake by strain (2 h g/kg). Bold italic slope values indicate significant (p < 0.05) effect of week on g/kg alcohol intake after linear regression. Intake is variable but generally increases over time for both the 2 and 4 h exposures in the BXD strains. However, slope at 2 h is not always identical to slope at 4 h and both traits are modestly correlated (r = 0.41). Parental strains are highlighted.

and slightly decreased 4 h intake in contrast to the D2 strain, which showed relatively stable or slightly decreased intake for both the 2 and 4 h intake period. However, these changes in

the parental strains were not statistically significant. Although modestly correlated (r = 0.41), change in 2 and 4 h intake over 5 weeks was not identical within strain. We observed a strong, negative correlation (−0.62) between slope (alcohol intake over time) and the y-intercept (baseline intake) for 4 h intake. As demonstrated in **Table 2**, strains with high initial intake tended to evince more stable or slightly decreased drinking over time. This negative relationship was present, albeit much weaker, for the 2 h time point (r = −0.25).

## Identification of QTLs for Alcohol Intake

We first assessed whether there were QTLs associated with the heritable variation in initial (week 1) or sustained weekly alcohol intake using two different mapping methods (traditional HK and GEMMA to correct for kinship). As expected, based on the inclusion of 39 strains, no significant QTLs were detected following multiple test correction (Supplementary Table S2). Using a suggestive threshold, seven QTLs were identified for HK and GEMMA (Table 3). QTLs detected by both methods are expected to be more robust than QTLs detected by HK alone, because these QTLs do not arise as a result of genetic relatedness. We prioritized the overlapping QTLs detected by both methods for further analysis. In all seven cases, inheritance of the B6 parental allele was associated with higher 2 or 4 h alcohol intake. QTLs on Chr 2, 6, and 12 were only identified for 4 h intake on week 1. In contrast, overlapping QTLs on Chr 4 were only identified for 2 h intake on weeks 2 and 3, and the QTL on Chr 8 was identified for 4 h intake on week 3 and 2 h intake on week 5.

Next, we identified QTLs for progressive alcohol intake using slope analysis of 2 or 4 h intake over 5-weeks. No significant or suggestive QTLs were identified for progressive 4 h intake. However, a suggestive QTL on Chr 15 was identified for progressive 2 h alcohol intake using both HK and GEMMA mapping methods (Table 3). In contrast to the other seven QTLs, the 2 h progressive intake QTL was associated with higher intake and inheritance of the D2 parental allele.

## Identification of Candidate Genes Driving Initial Alcohol Intake

The majority of suggestive overlap QTLs detected were associated with variation in initial 4 h alcohol intake on week 1 (W1.4 h, Table 3). Thus we prioritized these QTLs on Chr 2 (9–30.4 Mb), 6 (89–93.6 Mb), and 12 (9.8–15.8 Mb) for identification of candidate genes driving variation in initial early alcohol intake. Hereafter, all QTLs are referred to by their QTL (Q) Chr number. The number of candidate genes located within each confidence interval varies [Q2 = 325, Q6 = 54, Q12 = 25; based on mouse assembly GRCm38/mm10 and UCSC RefSeq (refGene) Table Browser annotations, Supplementary Table S3]. Using genome sequence data generated for the parental strains we identified all genes (excluding gene models, predicted genes, non-coding and pseudo-genes) within each interval that were polymorphic and overlapped by higher impact sequence variants (see Methods). This reduced the number of candidate genes within each interval to Q2 = 60, Q6 = 13, and Q12 = 4 (Supplemental Table S4). These genes were further prioritized using naïve hippocampal
<table><caption>**TABLE 3** | Suggestive QTL detected by multiple mapping methods.</caption><thead><tr><th>Time</th><th>Method</th><th>Chr</th><th>Locus</th><th>Mb</th><th>Additive Effect</th><th>LOD</th><th>CIMb (left)</th><th>CIMb (right)</th></tr></thead><tbody><tr><td>W1.4 h</td><td>HK</td><td>2</td><td>rs13476358</td><td>16.306</td><td>-0.459</td><td>2.43</td><td>9.336</td><td>42.674</td></tr><tr><td>W1.4 h</td><td>GEMMA</td><td>2</td><td>Affy_PC2_15</td><td>15.000</td><td>-</td><td>3.07</td><td>9.000</td><td>30.441</td></tr><tr><td>W2.2 h</td><td>HK</td><td>4</td><td>rs32939068</td><td>153.934</td><td>-0.247</td><td>2.46</td><td>151.094</td><td>156.101</td></tr><tr><td>W2.2 h</td><td>GEMMA</td><td>4</td><td>rs32939068</td><td>153.934</td><td>-</td><td>3.44</td><td>153.372</td><td>155.226</td></tr><tr><td>W3.2 h</td><td>HK</td><td>4</td><td>rs32939068</td><td>153.934</td><td>-0.303</td><td>2.81</td><td>151.094</td><td>156.101</td></tr><tr><td>W3.2 h</td><td>GEMMA</td><td>4</td><td>rs32939068</td><td>153.934</td><td>-</td><td>4.08</td><td>153.372</td><td>155.493</td></tr><tr><td>W1.4 h</td><td>HK</td><td>6</td><td>rs30422489</td><td>91.707</td><td>-0.508</td><td>3.08</td><td>89.011</td><td>93.600</td></tr><tr><td>W1.4 h</td><td>GEMMA</td><td>6</td><td>rs30422489</td><td>91.707</td><td>-</td><td>3.23</td><td>89.011</td><td>93.600</td></tr><tr><td>W3.4 h</td><td>HK</td><td>8</td><td>rs49907965</td><td>89.094</td><td>-0.573</td><td>2.67</td><td>80.460</td><td>95.747</td></tr><tr><td>W3.4 h</td><td>GEMMA</td><td>8</td><td>rs49907965</td><td>89.094</td><td>-</td><td>3.16</td><td>86.380</td><td>95.747</td></tr><tr><td>W5.2 h</td><td>HK</td><td>8</td><td>rs49907965</td><td>89.094</td><td>-0.325</td><td>2.77</td><td>82.870</td><td>95.736</td></tr><tr><td>W5.2 h</td><td>GEMMA</td><td>8</td><td>rs49907965</td><td>89.094</td><td>-</td><td>3.23</td><td>82.870</td><td>95.747</td></tr><tr><td>W1.4 h</td><td>HK</td><td>12</td><td>rs3717933</td><td>11.076</td><td>-0.488</td><td>2.76</td><td>9.864</td><td>15.820</td></tr><tr><td>W1.4 h</td><td>GEMMA</td><td>12</td><td>rs49972008</td><td>12.603</td><td>-</td><td>3.07</td><td>7.991</td><td>26.035</td></tr><tr><td>Slope. 2 h</td><td>HK</td><td>15</td><td>rs31691968</td><td>63.392</td><td>0.045</td><td>3.19</td><td>58.000</td><td>67.990</td></tr><tr><td>Slope. 2 h</td><td>GEMMA</td><td>15</td><td>rs31691968</td><td>63.392</td><td>-</td><td>3.69</td><td>58.000</td><td>67.990</td></tr></tbody></table>

QTLs were mapped for each week (W) of 2 or 4 h alcohol intake or for progressive 2 or 4 h alcohol intake based on slope analysis over five weeks of the DID paradigm. Additive effect, based on HK, mapping, is the contribution of parental alleles to mean trait expression. Negative additive values are associated with higher expression of the B6 allele. CI indicates the 1.5 LOD drop confidence interval. Text color indicates whether higher intake is driven by the inheritance of the B6 (blue) or D2 (red) parental allele at the locus. All QTL were detected at the suggestive threshold.

expression data generated for the BXD family (Table 4 and Supplemental Table S5). For Q2, Arhgap21, Gpr158, Mrpl41, Myo3a, Entpd4, Lhx3, Rapgef1, and Nup188 were all significantly correlated with alcohol intake during the first 4 h exposure in week 1, and the expression of each candidate was also regulated by local sequence variants (cis eQTL). Top candidates for Q6 included Adamts9 and Fgd5. For Q12, Nt5c1b was the only candidate that contained putative high impact variants and whose expression covaried with 4 h alcohol intake on week 1. Products of all candidate genes in Q2, Q6, and Q12 played a role in signal transduction, metabolism, development, or endothelial cell response.

## Identification of Candidate Genes Driving Variation in Sustained Weekly Alcohol Intake

The same strategy used to identify candidates driving variation in initial consumption was used to identify candidates on Q4 (151–155 Mb) and Q8 (86.4–95.7 Mb) associated with variation in 2 h intake on weeks 2 and 3, and variation in 4 h intake on week 3 and 2 h intake on week 5, respectively. The number of candidates within each interval was 51 for Q4 and 93 for Q8 (Supplementary Table S3). The number of candidate genes overlapping higher impact variants within each locus was 36 for Q8 and 16 for Q4 (Supplementary Table S4). For the Q4 interval, the highest priority candidates based on genetic cis-regulation and co-variation analysis using hippocampal expression data

were Per3 and Prdm16 (Table 4 and Supplementary Table S5). Top candidates for the Q8 interval included Cyld, Aktip, Fto, Crnde, Mmp2, and Kifc3. All of these high priority candidate genes were involved in circadian signaling, metabolism, or immune response.

## Identification of Candidate Genes and Pathways Driving Variation in Progressive Alcohol Intake

Variation in progressive 2 h alcohol intake was associated with a suggestive QTL on Q15 (58–68 Mb) that contained 54 candidate genes (Supplementary Table S3). Only four of these genes (Adcy8, Tg, Ndrg1, and Wisp) were overlapped by variants of predicted higher impact (Supplementary Table S4). None of these candidates were modulated by cis eQTLS in naïve BXD hippocampus or correlated with the slope for 2 h alcohol intake over 5 weeks (Supplementary Table S5).

## Trait Covariation With Initial Alcohol Intake

Trait data have been collected for the BXDs since the generation of the first BXD cohort in the late 1970's and much of this data is available in the BXD Published Phenotypes database available on GN. We queried this database to retrieve traits that were significantly correlated with W1.4 h, a heritable and QTL modulated trait associated with alcohol intake at early exposure. Forty-one traits were significantly correlated ($p < 0.005$) with the
<table><caption>TABLE 4 | Prioritized QTL candidate genes.</caption><thead><tr><th>Symbol</th><th>Description</th><th>Gene Wiki</th><th>QTL</th><th>HIP:Mean Expr</th><th>HIP:ciseQTL</th><th>HIP:Trait Correlation</th></tr></thead><tbody><tr><td>Dnajc1</td><td>DnaJ (Hsp40) homolog. subfamily C, member 1</td><td>GRP78 partner in ER; translation, beneficial effects in ER stress induced metabolic dysfunction</td><td>Q2: Initial</td><td>8.74</td><td>Sug</td><td>-</td></tr><tr><td>Armc3</td><td>armadillo repeat containing 3</td><td>?</td><td>Q2: Initial</td><td><strong>5.64</strong></td><td>Sig</td><td>-</td></tr><tr><td><strong>Arhgap21</strong></td><td>Rho GTPase activating protein 21</td><td>Signaling</td><td>Q2: Initial</td><td>11.76</td><td>Sig</td><td>Sig</td></tr><tr><td><strong>Gpr158</strong></td><td>G protein-coupled receptor 158</td><td>Regulator of RGS complexes and Rgs7 signaling in brain; human variants may influence energy expenditure and metabolism</td><td>Q2: Initial</td><td>6.51</td><td>Sig</td><td>Sig</td></tr><tr><td>Myo3a</td><td>Myosin IIIA</td><td>Acin based motor protein with kinase activity</td><td>Q2: Initial</td><td>6.35</td><td>-</td><td>Sig</td></tr><tr><td><strong>Mrpl41</strong></td><td>Mitochondrial ribosomal protein L41</td><td>Regulation of cell death</td><td>Q2: Initial</td><td>10.29</td><td>Sig</td><td>Sig</td></tr><tr><td>Uap1l1</td><td>UDP-N-acetylglucosamine pyrophorylase 1 like 1</td><td>?</td><td>Q2: Initial</td><td>9.01</td><td>Sig</td><td>-</td></tr><tr><td><strong>Entpd2</strong></td><td>Ectonucleoside triphosphate diphosphohydrolase 2</td><td>Metabolism</td><td>Q2: Initial</td><td>9.72</td><td>Sig</td><td>Sig</td></tr><tr><td>Ptgds</td><td>Prostaglandin D2 synthase</td><td>Glucose and insulin metabolism; regulated by estradiol: inflammatory response</td><td>Q2: Initial</td><td>13.90</td><td>Sig</td><td>-</td></tr><tr><td>C8g</td><td>Complement component 8, gamma subunit</td><td>Immune function</td><td>Q2: Initial</td><td>7.13</td><td>Sig</td><td>-</td></tr><tr><td>Kcnt1</td><td>Potassium channel, subfamily T, member 1 (slack, low threshold slowly adapting)</td><td>Localized to the postsynaptic density; involved in learning and memory and initial response to novel situations and environments</td><td>Q2: Initial</td><td>8.98</td><td>Sig</td><td>-</td></tr><tr><td>Lhx3</td><td>LIM homeobox protein 3</td><td>Transcription factor activity, nervous system development</td><td>Q2: Initial</td><td>6.23</td><td>-</td><td>Sig</td></tr><tr><td>Egfl7</td><td>EGF-like domain 7</td><td>Development</td><td>Q2: Initial</td><td>7.81</td><td>Sug</td><td>-</td></tr><tr><td>Rapgef1</td><td>Rap guanine nucleotide exchange factor (GEF) 1</td><td>Signaling; GABAergic neuronal development</td><td>Q2: Initial</td><td>7.41</td><td>Sug</td><td>Sig</td></tr><tr><td>Nup188</td><td>Nucleoporin 188</td><td>Development</td><td>Q2: Initial</td><td>8.50</td><td>-</td><td>Sig</td></tr><tr><td>Plxna1</td><td>plexin A1</td><td>Development; signaling</td><td>Q6: Initial</td><td>9.77</td><td>Sig</td><td>-</td></tr><tr><td>Fgd5</td><td>FYVE, RhoGEF and PH domain containing 5</td><td>Endothelial specific gene</td><td>Q6: Initial</td><td>8.00</td><td>-</td><td>Sig</td></tr><tr><td>Adamts9</td><td>A disintegrin-like and metalloprotease (reprolysin type) with thrombospondin type 1 motif, 9</td><td>Endothelial and vascular response</td><td>Q6: Initial</td><td>9.38</td><td>Sug</td><td>Sig</td></tr><tr><td>Nt5c1b</td><td>5'-nucleotidase, cytosolic 1B</td><td>Metabolism</td><td>Q12: Initial</td><td>7.18</td><td>Sug</td><td>Sig</td></tr><tr><td><strong>Per3</strong></td><td>Period 3</td><td>Circadian signaling</td><td>Q4: Sustained</td><td>7.11</td><td>Sig</td><td>Sig</td></tr><tr><td><strong>Prdm16</strong></td><td>PR domain containing 16</td><td>Metabolism; brown versus white fat differentiation</td><td>Q4: Sustained</td><td>6.52</td><td>Sig</td><td>Sig</td></tr><tr><td>Tnfrsf14</td><td>Tumor necrosis factor receptor superfamily, member 14 (herpesvirus entry mediator)</td><td>Inflammation and immune response</td><td>Q4: Sustained</td><td><strong>5.20</strong></td><td>Sig</td><td>-</td></tr><tr><td>Cyld</td><td>Cylindromatosis (turban tumor syndrome)</td><td>Immune; regulator of NF-kappaB signaling</td><td>Q8: Sustained</td><td>10.41</td><td>Sig</td><td>Sig</td></tr><tr><td>Chd9</td><td>Chromodomain helicase DNA binding protein 9</td><td>Expressed in osteoprogenitors</td><td>Q8: Sustained</td><td>9.55</td><td>Sig</td><td>-</td></tr><tr><td><strong>Aktip</strong></td><td>AKT interacting protein</td><td>Ft1 protein; telomere maintenance; development</td><td>Q8: Sustained</td><td>7.49</td><td>Sig</td><td>Sig</td></tr></tbody></table>

(Continued)
<table><caption>**TABLE 4** | Continued</caption><thead><tr><th>Symbol</th><th>Description</th><th>Gene Wiki</th><th>QTL</th><th>HIP:Mean Expr</th><th>HIP:cis eQTL</th><th>HIP:Trait Correlation</th></tr></thead><tbody><tr><td>Fto</td><td>Fat mass and obesity associated (alpha-ketoglutarate-dependent dioxygenase FTO)</td><td>Metabolism; energy homeostasis; regulator of adipogenesis</td><td>Q8: Sustained</td><td>6.94</td><td>Sig</td><td>Sig</td></tr><tr><td>Crnde</td><td>Colorectal neoplasia differentially expressed (non-protein coding)</td><td>Metabolism; control of glucose and lipid metabolism</td><td>Q8: Sustained</td><td>6.63</td><td>Sig</td><td>Sig</td></tr><tr><td>Mmp2</td><td>Matrix metaloproteinase 2</td><td>Immune system; inflammatory response</td><td>Q8: Sustained</td><td>6.66</td><td>Sig</td><td>Sig</td></tr><tr><td>Coq9</td><td>Coenzyme QB</td><td>Metabolism: deficiency impacts mitochondrial function</td><td>Q8: Sustained</td><td>10.49</td><td>Sig</td><td>Sig</td></tr><tr><td>Kifc3</td><td>Kinesin family member C3</td><td>Peroxisomal transport</td><td>Q8: Sustained</td><td>10.01</td><td>Sig</td><td>Sig</td></tr></tbody></table>

High priority candidate genes selected based on regulation of expression by local sequence variants (cis eQTL) and co-expression with DID intake traits (Slope2 hWeek, 2 h.W3avg, 2 h.W5avg, and W1.4 h). See **Supplementary Table S5** for full results. Gene expression signatures for each gene generated on the Affymetrix M430 platform for BXD hippocampus (HIP, naïve to treatment). Genes in bold have hippocampal gene expression that is significantly regulated by a cis eQTL (Sig) and significantly (p < 0.05) coexpressed (Sig; based on Pearson's r-value) with DID alcohol intake traits. Bold text in the HIP:MeanExpr column used to indicate genes with very low expression in hippocampus. In the hippocampal data set mean expression is 8 ± 2 log₂ units. Function summarized using GN Gene Wiki tool.

first 4 h alcohol intake trait. Representative scatterplots are shown in **Figure 2**. Correlated traits were subdivided into six categories: behavior; blood and brain chemistry, and hematology; drug response (morphine and alcohol); immune function; metabolism, and microbiome (**Supplementary Table S6**).

## DISCUSSION

Here, we report large variability in 2 and 4 h alcohol consumption in B6, D2, and 39 BXD strains using a DID protocol over a 5-week time period. BXD strains demonstrated high, low, or intermediate alcohol intake relative to parental strains (**Figure 1**). Regression analysis revealed that, in general, alcohol intake increased from weeks 1 to 5 regardless of session duration. Although it is important to note that some strains demonstrated greatly increased intake while others remained relatively stable or decreased slightly over the 5-week period. Our study both replicated previous findings on DID alcohol intake and generated novel insight into progressive drinking. In contrast to our study, most of the previous DID studies in male and female inbred strains focused primarily on the first few weeks of DID. These studies reported higher intake in B6 compared to D2, and that alcohol consumption “stabilized” by the second day of access (Rhodes et al., 2005). When these studies were extended for 12 consecutive days in male B6 mice, intake levels were reported to be constant (Rhodes et al., 2005). We also reported stabilization of alcohol intake early in week 1 in that heritability was reduced at first exposure compared to subsequent exposures. This was a result of greater within strain variation upon first experience with alcohol. Accordingly, we also observed much higher intake of 20% alcohol at 2 and 4 h in the B6 parental strain relative to D2. By lengthening the DID protocol to 5 weeks, our work extends the findings of previous studies and we reported increased consumption and, at least for 2 h intake, higher heritability (less within strain variation) during later weeks (3, 4, and 5). Similar increases in intake have been reported for B6 males over a

14-day DID access period (Linsenbardt and Boehm, 2014). Taken together, these results provide evidence that extended access DID protocols of 2 weeks or longer result in increased alcohol intake in B6J. For the first time, we report that this trend of enhanced consumption during extended DID is also apparent in the BXD family of strains. Patterns of initial weekly and progressive intake during 2 or 4 h sessions over 5 weeks varied between B6, D2, and the BXD strains. This variability was heritable and we were able to identify several suggestive loci for variation in initial intake during week 1, and variation in sustained and progressive intake.

The DID paradigm typically employs both a 2 h and a longer 4 h access period. In our study, some strains displayed marked differences in alcohol intake between access periods and identified QTLs were not always consistently detected for both periods. There are several explanations for the lack of complete congruency between the 2 and 4 h measures. First, the 2 h time point is a repeated measure and is thus a much more stable measure of weekly progressive intake. In addition, the extended 4 h access time period probably reflects a larger behavioral repertoire related to differences in learning and consummatory behavior. For example, male B6 mice allowed to drink alcohol in a modified DID paradigm in which they were given a 2 h access period over 14 consecutive days not only increased their intake over sessions, but also learned to “front load” their intake, and consumed the majority of their alcohol on the last session within the first few minutes of the 2 h access period (Linsenbardt and Boehm, 2014). In our study using female mice, we observed increased intake over the 5-week period in B6 females for the 2 h time point (**Table 2**), but the same trend was not observed for the 4 h time point. However, several strains, such as BXD 40, 29, and 34, showed increased intake over 5 weeks for both the 2 and 4 h sessions. It is unknown whether B6 females display front loading, but differences in this behavior, or a ceiling effect on intake in the 4 h access period, may be an additional source of variation that contributes to differences between the 2 and 4 h access period between and within strain in our study. For genetic mapping, the 2 h trait has higher heritability and is more amenable to mapping.
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

**FIGURE 2** | Representative scatterplots representing correlations between BXD legacy trait data and alcohol intake during the first 4 h exposure on week 1. (A) A measure of depression (immobility in the tail suspension test) is negatively correlated with alcohol intake. (B) Fasting sphingomyelin levels are also negatively correlated with intake. (C) A measure of immune function, infection survival time is positively correlated with intake. (D) Measures of one member of the gut microbiome are negatively correlated with intake. (E) A metabolic/immune trait, liver injury following insult, is negatively correlated with intake.

However, it is also possible that different genetic factors control each trait. Temporal dissection of intake over the 2 and 4 h access period for strains with similar and dissimilar 2 and 4 h intake would be needed to reconcile some of these differences.

Our study included 39 BXDs strains with a variable number of replicates within strain (1–16). Traditional QTL mapping incorporates the mean of the dependent measure for each BXD strain. Power to detect QTLs is derived largely from the number of strains (Belknap, 1998; Andreux et al., 2012). However, the traditional mapping method used for BXD data does not account

for population structure and may result in false positives. To account for this, we mapped QTLs using both a simple regression method (HK) in GN that does not account for kinship, and another mapping algorithm (GEMMA) available in GN, version 2, that does. Even though GEMMA is not typically used for smaller data sets, several QTLs were replicated using both methods, albeit at a suggestive level. Both methods used the same marker panel and dependent variables for alcohol intake. QTLs detected by both methods are expected to be more robust against identification of false positives due to genetic relatedness
compared to using either model alone (see methods). With 39
strains and a variable number of replicates within strain, our
study was only powered to detect QTLs of large effect. As
expected, we did not identify any significant loci. However,
using both mapping methods we identified: suggestive QTLs
associated with variation in alcohol intake during early exposure
(W1.4 h) on Chr 2, 6, and 12, and suggestive QTLs associated
with variation in sustained intake during later weeks on Chr 4
and 8. We were also able to identify a QTL on Chr 15 associated
with progressive alcohol intake. To assess whether these QTLs
are likely to reach significance if more strains are added we
winsorized (transformation of the data by limiting extreme
values) the existing 39 strain data to better fit assumptions
of normality assumed when using HK and GEMMA mapping
methods. Winsorizing the data resulted in higher LOD scores
for all QTLs for both mapping methods. This suggests that,
assuming a normal trait distribution in the BXD population,
testing a larger number of BXD strains will result in significant
LOD scores. Thus, the provisional QTLs identified in our study
are worth investigating further and may contain genes that drive
higher or lower initial or sustained alcohol consumption. As such,
we used multiple lines of evidence based on overlapping gene
variants and gene expression to nominate candidate genes for
suggestive QTLs identified in this study. An important caveat
of this approach is that mutations that impact protein function
without modulating gene expression will be missed. However,
this is a powerful and integrative approach that nominates the
best candidates based on all available data. Candidates associated
with genetic predisposition for high or low intake during the 1st
week of exposure to alcohol included *Arhgap21*, *Gpr158*, *Mrpl41*,
*Myo3a*, *Entpd2*, *Lhx3*, *Rapgef1*, and *Nup188* (Q2); *Adamts9* and
*Fgd5* (Q6); and *Nt5c1b* (Q12). All candidates were overlapped by
variants that may impact gene regulation or function, and their
expression in hippocampus was cis-modulated and correlated
with 4 h intake during week 1. Therefore, pre-existing patterns
of expression in these genes due to the presence of functional
gene variants could alter initial alcohol intake behavior. Only
one of these genes, Entpd2, had been directly implicated in
alcohol-related responses in animal model systems (Rico et al.,
2008). Rico and colleagues found that acute alcohol exposure led
to alterations in brain nucleotide (specifically, ATP hydrolysis)
metabolism *in vivo* and a reduction in the mRNA level of
the zebrafish homolog of ENTPD2 (NTPDase2). Although not
directly associated with alcohol related traits in human or
animal model systems, *Nt5c1b* was also involved in nucleotide
metabolism (specifically, adenosine formation following ATP
hydrolysis; Sala-Newby and Newby, 2001). Other candidates were
involved in intracellular signaling pathways (*Arhgap21*, *Gpr158*,
*Rapgef1*, *Adamts9*, *Fgd5*, *Mrpl41*, *Nt5c1b*) that mediated diverse
biological processes related to metabolism, G-protein coupled
receptor signaling, and endothelial response (Table 4). Taken
together, these results suggest that alterations in ATP metabolism
and other signaling pathways could drive variation in initial
alcohol intake in the DID paradigm.

Using a similar approach we also identified *Per3* and *Prdm16*
(Q4) and *Cyld*, *Aktip*, *Fto*, *Crnde*, *Mmp2*, and *Kifc3* (Q8) as
high priority candidate genes for variation (higher or lower) in
sustained intake during later weeks. Variation in 2 h intake on

weeks 2 and 3 was associated with Q4, and variation in 4 h
intake on week 3 and 2 h intake on week 5 was associated
with Q8. Importantly, the Q4 and Q8 intervals were detected
at multiple time points over the course of the 5-week DID
study. These concordant and reproducible QTLs are less likely
to represent false positives. Candidates for the Q4 and Q8
interval play key roles in circadian signaling (*Per3*), metabolic
(*Prdm16*, *Fto*, *Crnde*), and immune response (*Mmp2*, *Cyld*).
Alterations in metabolism (e.g., mitochondrial dysfunction,
insulin/glucose dysregulation), circadian cycle, sleep disruptions,
and inflammation have long been associated with chronic alcohol
use, abuse, and alcoholism. In particular, several of these genes,
*Per3* (Q4) and *Fto*, and *Mmp2* (Q8) have all been previously
associated with ethanol response. In human populations, *PER3*
variants were associated with insomnia severity in alcohol-
dependent patients (Brower et al., 2012). In the BXD panel,
a promoter mutation in *Per3* was associated with expression
variation in naïve animals (cis-modulation) as well as variation
in alcohol-related traits and expression changes following stress
and alcohol exposure (Wang et al., 2012). Altered circadian
expression of *Per3* was also associated with administration of a
liquid alcohol diet in rats (Chen et al., 2004). Another strong
candidate for sustained intake was *Fto*, a gene regulating fat
mass, located on Chr 8 and associated with variation in intake
during week 5. *Fto* variants were associated with frequency
of alcohol consumption (Sobczyk-Kopciol et al., 2011; Young
et al., 2016) and alcohol dependence (Wang et al., 2013) in
human populations. *Mmp2*, also located on Chr 8, is known
to be regulated by alcohol exposure in numerous human and
animal systems and in many different cell and tissue types
(Partridge et al., 1999; Burnham et al., 2007; Fiotti et al., 2008;
Peng et al., 2013), however, a role for this gene in regulating
alcohol intake is unclear. Here we report *Per3* and *Fto* as a
strong candidate genes located on Chr 4 and Chr 8, respectively,
that may underly variation in sustained alcohol intake during
DID in the BXD family. Both were associated with alcohol-
related traits in human and rodent model systems, and their
expression was correlated with 2 h intake during week 3 (negative
correlation, *Per3*) or week 5 (positive correlation, *Fto*) of DID.
In addition, *Per3* contains functional polymorphic regulatory
variants known to impact expression in the BXD panel and *Fto*
is a cis-modulated gene containing multiple missense variants
between B6 and D2.

By extending the DID paradigm to a longer 5-week period
we were also able to identify a suggestive QTL on Chr 15
for variation (higher or lower) in progressive alcohol intake.
None of the genes located within the QTL confidence interval
was cis-modulated in BXD hippocampus. This was an expected
negative result based on the use of an expression data set
generated from naïve animals. However, there were only a
handful of genes (*Adcy8*, *Tg*, *Ndgr1*, and *Wisp*) within this
interval that contained variants likely to impact gene regulation
or function. Of these, only *Adcy8* was previously associated
with alcohol-related traits (*Procopio et al., 2013*). *Procopio* and
colleagues reported that a haplotype within the coding region
of ADCY8 was associated with alcohol dependence co-occurring
with depression in females. Given that progressive alcohol intake
likely reflects changing alcohol consumption over time, it is
reasonable to expect that the underlying gene variants may interact with alcohol exposure to produce variation in intake over time. Of interest, genes located within the Q15 interval were significantly enriched (adjusted *p* < 0.05) for HDAC2 and PPARD binding based on consensus genome-wide ChIP-X data from ENCODE and integrated ChIP Enrichment Analysis (ChEA). None of the other QTL intervals (Q2, Q4, Q6, Q8, or Q12) were enriched for transcription factor binding or histone modifications based on Enrichr analysis, suggesting that gene variants in the Q15 region may exert their effect through alcohol experience dependent epigenetic modulation.

This is the first work of its kind to elucidate genetic influence on ethanol consumption at first opportunity and repeated opportunities over several weeks in a large genetic reference population of mice. As in other studies, the DID paradigm has been applied successfully to reveal genetic influences on alcohol consumption and the impact of alcohol intake on brain gene transcription (Rhodes et al., 2005, 2007; Mulligan et al., 2011; Thiele and Navarro, 2014; Thiele et al., 2014). The DID paradigm allows the animal to achieve significant blood ethanol concentrations, i.e., binge drinking. Indeed, we found that some strains consume large amounts of alcohol and that most of the strains studied increased their consumption by more than 0.5 g/kg. However, a caveat of our study is that we did not measure blood alcohol concentrations (BACs) weekly. We did not attempt to obtain BACs because of the added stressor of blood collection. However, this precludes an analysis of metabolic effects that may have contributed to changes in alcohol intake among strains.

One of the great advantages of a systems approach in the BXD panel is the ability to find associations between a trait of interest and other phenotypes generated in-house, and by others, in the extensive BXD trait database available at GN. Not only were we able to nominate candidate genes for suggestive QTL related to initial, sustained, and progressive alcohol intake, but we were also able to leverage 1000s of legacy trait records available for the BXD strains to associate variation in initial alcohol consumption with variation in metabolic, metagenomic, immune, and behavioral traits in the BXDs (Supplementary Table S4). Some of these associations were expected (e.g., correlations with depression- and anxiety-related traits and alcohol intake; and covariation with morphine and other alcohol response traits), while others were novel and may uncover new areas of interest. For example, and perhaps most surprising, a negative correlation between components of the gut microbiome and initial 4 h alcohol intake (Figure 2D). Of interest, Prevotellaceae was found to be enriched in fecal samples from alcohol cirrhosis patients relative to control samples (Chen et al., 2011). Oral levels of Lactobacillales and Prevotellaceae were also altered following alcohol intake in humans (Fan et al., 2018). Taken together these novel findings could indicate the presence of latent

relationships between alcohol-related diseases and members of the microbial community. This opens the research area to what these associations mean in the general physiology and behavior of mice and other species, by hypothesis generation. In conclusion, we have shown here that there is a large genetic component to both initial, sustained, and progressive consumption of alcohol in the DID paradigm among the BXD family of mice. Moreover, through forward genetic mapping and systems biology, we have identified candidate genes and pathways likely to mediate high alcohol intake.

## AUTHOR CONTRIBUTIONS

MKM, MD, and BCJ wrote the manuscript. MKM, DA, PP, WZ, and BCJ analyzed the data. LL, BCJ, WZ, MKM, ET, SAC, and PM designed and implemented the study.

## ACKNOWLEDGMENTS

We would like to acknowledge Dr. Rob Williams for his intellectual support, contribution of BXD strains, and for GN development. We also acknowledge Arthur Centeno for his assistance with data integration into GN. Partial support for this study was provided by USPHS (NIH) Grants U01AA014425 and R01AA021951.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2018.00370/full#supplementary-material

**FIGURE S1** | Bivariate correlations of alcohol intake at initial exposure and across 2 or 4 h weekly intake. Correlations based on weekly (W) and daily (D) intake (g/kg) for individuals from B6, D2, and 39 BXD strains. Data averaged by strain for the first day of 2 h alcohol exposure (2 h:W1D1) and for each 4 h exposure. Data averaged by strain and week for subsequent 2 h exposures. All correlations were positive. Intake at first exposure for 2 h on W1D1 is highly correlated with subsequent average 2 and 4 h weekly intake. Correlations are higher between adjacent weeks and generally decay slightly over time. Compared to correlations between 4 h intake on the first and last exposure, correlations between intake at first (W1D1 or W1avg) and last exposure are higher for 2 h intake.

**TABLE S1** | Strain and replicate summary.

**TABLE S2** | QTL comparison table.

**TABLE S3** | Candidate genes located within QTL confidence intervals.

**TABLE S4** | Candidate genes within QTL confidence intervals overlapped by higher impact sequence variants.

**TABLE S5** | Full list of high priority candidate genes and evidence for selection.

**TABLE S6** | BXD trait correlations.

## REFERENCES

Andreux, P. A., Williams, E. G., Koutnikova, H., Houtkooper, R. H., Champy, M. F., Henry, H., et al. (2012). Systems genetics of metabolism: the use of the BXD

murine reference panel for multiscalar integration of traits. *Cell* 150, 1287–1299.
doi: 10.1016/j.cell.2012.08.012

Babor, T. F., Dolinsky, Z. S., Meyer, R. E., Hesselbrock, M., Hofmann, M., and Tennen, H. (1992). Types of alcoholics: concurrent and predictive validity of
some common classification schemes. *Br. J. Addict.* 87, 1415-1431. doi: 10.1111/j.1360-0443.1992.tb01921.x

Belknap, J. K. (1998). Effect of within-strain sample size on QTL detection and mapping using recombinant inbred mouse strains. *Behav. Genet.* 28, 29-38. doi: 10.1023/A:1021404714631

Belknap, J. K., Metten, P., Helms, M. L., O'Toole, L. A., Angeli-Gade, S., Crabbe, J. C., et al. (1993). Quantitative trait loci (QTL) applications to substances of abuse: physical dependence studies with nitrous oxide and ethanol in BXD mice. *Behav. Genet.* 23, 213-222. doi: 10.1007/BF01067426

Brower, K. J., Wojnar, M., Sliwerska, E., Armitage, R., and Burmeister, M. (2012). PER3 polymorphism and insomnia severity in alcohol dependence. *Sleep* 35, 571-577. doi: 10.5665/sleep.1748

Browman, K. E., and Crabbe, J. C. (2000). Quantitative trait loci affecting ethanol sensitivity in BXD recombinant inbred mice. *Alcohol. Clin. Exp. Res.* 24, 17-23. doi: 10.1111/j.1530-0277.2000.tb04547.x

Buck, K. J., Metten, P., Belknap, J. K., and Crabbe, J. C. (1997). Quantitative trait loci involved in genetic predisposition to acute alcohol withdrawal in mice. *J. Neurosci.* 17, 3946-3955. doi: 10.1523/JNEUROSCI.17-10-03946.1997

Burnham, E. L., Moss, M., Ritzenthaler, J. D., and Roman, J. (2007). Increased fibronectin expression in lung in the setting of chronic alcohol abuse. *Alcohol. Clin. Exp. Res.* 31, 675-683. doi: 10.1111/j.1530-0277.2007.00352.x

Chen, C. P., Kuhn, P., Advis, J. P., and Sarkar, D. K. (2004). Chronic ethanol consumption impairs the circadian rhythm of pro-opiomelanocortin and period genes mRNA expression in the hypothalamus of the male rat. *J. Neurochem.* 88, 1547-1554. doi: 10.1046/j.1471-4159.2003.02300.x

Chen, H., and Sharp, B. M. (2004). Content-rich biological network constructed by mining PubMed abstracts. *BMC Bioinformatics* 5:147. doi: 10.1186/1471-2105-5-147

Chen, Y., Yang, F., Lu, H., Wang, B., Chen, Y., Lei, D., et al. (2011). Characterization of fecal microbial communities in patients with liver cirrhosis. *Hepatology* 54, 562-572. doi: 10.1002/hep.24423

Chesler, E. J., Lu, L., Shou, S., Qu, Y., Gu, J., Wang, J., et al. (2005). Complex trait analysis of gene expression uncovers polygenic and pleiotropic networks that modulate nervous system function. *Nat. Genet.* 37, 233-242. doi: 10.1038/ng1518

Churchill, G. A., and Doerge, R. W. (1994). Empirical threshold values for quantitative trait mapping. *Genetics* 138, 963-971.

Cloninger, C. R. (1987). Neurogenetic adaptive mechanisms in alcoholism. *Science* 236, 410-416. doi: 10.1126/science.2882604

Crabbe, J. C. (1998). Provisional mapping of quantitative trait loci for chronic ethanol withdrawal severity in BXD recombinant inbred mice. *J. Pharmacol. Exp. Ther.* 286, 263-271.

Crabbe, J. C., Metten, P., Belknap, J. K., Spence, S. E., Cameron, A. J., Schlumbergh, J. P., et al. (2014). Progress in a replicated selection for elevated blood ethanol concentrations in HDID mice. *Genes Brain Behav.* 13, 236-246. doi: 10.1111/gbb.12105

Crabbe, J. C., Phillips, T. J., Gallaher, E. J., Crawshaw, L. I., and Mitchell, S. R. (1996). Common genetic determinants of the ataxic and hypothermic effects of ethanol in BXD/Ty recombinant inbred mice: genetic correlations and quantitative trait loci. *J. Pharmacol. Exp. Ther.* 277, 624-632.

Cunningham, C. L. (1995). Localization of genes influencing ethanol-induced conditioned place preference and locomotor activity in BXD recombinant inbred mice. *Psychopharmacology* 120, 28-41. doi: 10.1007/BF02246142

Fan, X., Peters, B. A., Jacobs, E. J., Gapstur, S. M., Purdue, M. P., Freedman, N. D., et al. (2018). Drinking alcohol is associated with variation in the human oral microbiome in a large study of American adults. *Microbiome* 6:59. doi: 10.1186/s40168-018-0448-x

Fernandez, J. R., Vogler, G. P., Tarantino, L. M., Vignetti, S., Plomin, R., and McClearn, G. E. (1999). Sex-exclusive quantitative trait loci influences in alcohol-related phenotypes. *Am. J. Med. Genet.* 88, 647-652. doi: 10.1002/(SICI) 1096-8628(19991215)88:6<647::AID-AJMG13>3.0.CO;2-6

Fiotti, N., Tubaro, F., Altamura, N., Grassi, G., Moretti, M., Dapas, B., et al. (2008). Alcohol reduces MMP-2 in humans and isolated smooth muscle cells. *Alcohol* 42, 389-395. doi: 10.1016/j.alcohol.2008.02.001

Gill, K., Liu, Y., and Deitrich, R. A. (1996). Voluntary alcohol consumption in BXD recombinant inbred mice: relationship to alcohol metabolism. *Alcohol. Clin. Exp. Res.* 20, 185-190. doi: 10.1111/j.1530-0277.1996.tb01063.x

Grisel, J. E., Metten, P., Wenger, C. D., Merrill, C. M., and Crabbe, J. C. (2002). Mapping of quantitative trait loci underlying ethanol metabolism in BXD recombinant inbred mouse strains. *Alcohol. Clin. Exp. Res.* 26, 610-616. doi: 10.1111/j.1530-0277.2002.tb02582.x

Keane, T. M., Goodstadt, L., Danecek, P., White, M. A., Wong, K., Yalcin, B., et al. (2011). Mouse genomic variation and its effect on phenotypes and gene regulation. *Nature* 477, 289-294. doi: 10.1038/nature10413

Kuleshov, M. V., Jones, M. R., Rouillard, A. D., Fernandez, N. F., Duan, Q., Wang, Z., et al. (2016). Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. *Nucleic Acids Res.* 44, W90-W97. doi: 10.1093/nar/gkw377

Linsenbardt, D. N., and Boehm, S. L. II (2014). Alterations in the rate of binge ethanol consumption: implications for preclinical studies in mice. *Addict. Biol.* 19, 812-825. doi: 10.1111/adb.12052

Mulligan, M. K., Mozhui, K., Prins, P., and Williams, R. W. (2017). GeneNetwork: a toolbox for systems genetics. *Methods Mol. Biol.* 1488, 75-120. doi: 10.1007/978-1-4939-6427-7_4

Mulligan, M. K., Ponomarev, I., Hitzemann, R. J., Belknap, J. K., Tabakoff, B., Harris, R. A., et al. (2006). Toward understanding the genetics of alcohol drinking through transcriptome meta-analysis. *Proc. Natl. Acad. Sci. U.S.A.* 103, 6368-6373. doi: 10.1073/pnas.051088103

Mulligan, M. K., Rhodes, J. S., Crabbe, J. C., Mayfield, R. D., Harris, R. A., and Ponomarev, I. (2011). Molecular profiles of drinking alcohol to intoxication in C57BL/6J mice. *Alcohol. Clin. Exp. Res.* 35, 659-670. doi: 10.1111/j.1530-0277.2010.01384.x

Overall, R. W., Kempermann, G., Peirce, J., Lu, L., Goldowitz, D., Gage, F. H., et al. (2009). Genetics of the hippocampal transcriptome in mouse: a systematic survey and online neurogenomics resource. *Front. Neurosci.* 3:55. doi: 10.3389/ neuro.15.003.2009

Partridge, C. R., Sampson, H. W., and Forough, R. (1999). Long-term alcohol consumption increases matrix metalloproteinase-2 activity in rat aorta. *Life Sci.* 65, 1395-1402. doi: 10.1016/S0024-3205(99)00381-1

Peirce, J. L., Lu, L., Gu, J., Silver, L. M., and Williams, R. W. (2004). A new set of BXD recombinant inbred lines from advanced intercross populations in mice. *BMC Genet.* 5:7. doi: 10.1186/1471-2156-5-7

Peng, C., Li, W. A., Fu, P., Chakraborty, T., Hussain, M., Guthikonda, M., et al. (2013). At low doses ethanol maintains blood-brain barrier (BBB) integrity after hypoxia and reoxygenation: a brain slice study. *Neurol. Res.* 35, 790-797. doi: 10.1179/1743132813Y.0000000198

Philip, V. M., Duvvuru, S., Gomero, B., Ansah, T. A., Blaha, C. D., Cook, M. N., et al. (2010). High-throughput behavioral phenotyping in the expanded panel of BXD recombinant inbred strains. *Genes Brain Behav.* 9, 129-159. doi: 10. 1111/j.1601-183X.2009.00540.x

Phillips, T. J., Crabbe, J. C., Metten, P., and Belknap, J. K. (1994). Localization of genes affecting alcohol drinking in mice. *Alcohol. Clin. Exp. Res.* 18, 931-941. doi: 10.1111/j.1530-0277.1994.tb00062.x

Phillips, T. J., Huson, M., Gwiazdon, C., Burkhart-Kasch, S., and Shen, E. H. (1995). Effects of acute and repeated ethanol exposures on the locomotor activity of BXD recombinant inbred mice. *Alcohol. Clin. Exp. Res.* 19, 269-278. doi: 10.1111/j.1530-0277.1995.tb01502.x

Phillips, T. J., Lessov, C. N., Harland, R. D., and Mitchell, S. R. (1996). Evaluation of potential genetic associations between ethanol tolerance and sensitization in BXD/Ty recombinant inbred mice. *J. Pharmacol. Exp. Ther.* 277, 613-623.

Procopio, D. O., Saba, L. M., Walter, H., Lesch, O., Skala, K., Schlaff, G., et al. (2013). Genetic markers of comorbid depression and alcoholism in women. *Alcohol. Clin. Exp. Res.* 37, 896-904. doi: 10.1111/acer.12060

Rhodes, J. S., Best, K., Belknap, J. K., Finn, D. A., and Crabbe, J. C. (2005). Evaluation of a simple model of ethanol drinking to intoxication in C57BL/6J mice. *Physiol. Behav.* 84, 53-63. doi: 10.1016/j.physbeh.2004. 10.007

Rhodes, J. S., Ford, M. M., Yu, C. H., Brown, L. L., Finn, D. A., Garland, T., et al. (2007). Mouse inbred strain differences in ethanol drinking to intoxication. *Genes Brain Behav.* 6, 1-18. doi: 10.1111/j.1601-183X.2006. 00210.x

Rico, E. P., Roseberg, D. B., Senger, M. R., de Bem Arizi, M., Dias, R. D., Souto, A. A., et al. (2008). Ethanol and acetaldehyde alter NTPDase and 5'-nucleotidase from zebrafish brain membranes. *Neurochem. Int.* 52, 290-296. doi: 10.1016/j.neuint.2007.06.034
Risinger, F. O., and Cunningham, C. L. (1998). Ethanol-induced conditioned taste aversion in BXD recombinant inbred mice. *Alcohol. Clin. Exp. Res.* 22, 1234-1244. doi: 10.1111/j.1530-0277.1998.tb03904.x

Roberts, A. J., Phillips, T. J., Belknap, J. K., Finn, D. A., and Keith, L. D. (1995). Genetic analysis of the corticosterone response to ethanol in BXD recombinant inbred mice. *Behav. Neurosci.* 109, 1199-1208. doi: 10.1037/0735-7044.109.6. 1199

Rodriguez, L. A., Plomin, R., Blizard, D. A., Jones, B. C., and McClearn, G. E. (1994). Alcohol acceptance, preference, and sensitivity in mice. *I. Quantitative genetic analysis using BXD recombinant inbred strains. Alcohol. Clin. Exp. Res.* 18, 1416-1422. doi: 10.1111/j.1530-0277.1994.tb01444.x

Sala-Newby, G. B., and Newby, A. C. (2001). Cloning of a mouse cytosolic 5'-nucleotidase-I identifies a new gene related to human autoimmune infertility-related protein. *Biochim. Biophys. Acta* 1521, 12-18. doi: 10.1016/S0167-4781(01)00278-0

Sloan, Z., Arends, D., Broman, K. W., Centeno, A., Furlotte, N. A., Nijveen, H., et al. (2016). GeneNetwork: framework for web-based genetics. *J. Open Source Softw.* 1:25. doi: 10.21105/joss.00025

Sobczyk-Kopciol, A., Broda, G., Wojnar, M., Kurjata, P., Jakubczyk, A., Klimkiewicz, A., et al. (2011). Inverse association of the obesity predisposing FTO rs9939609 genotype with alcohol consumption and risk for alcohol dependence. *Addiction* 106, 739-748. doi: 10.1111/j.1360-0443.2010.03248.x

Thiele, T. E., Crabbe, J. C., and Boehm, S. L. II (2014). "Drinking in the Dark" (DID): a simple mouse model of binge-like alcohol intake. *Curr. Protoc. Neurosci.* 68, 9.49.1-9.49.12. doi: 10.1002/0471142301.ns0949s68

Thiele, T. E., and Navarro, M. (2014). "Drinking in the dark" (DID) procedures: a model of binge-like ethanol drinking in non-dependent mice. *Alcohol* 48, 235-241. doi: 10.1016/j.alcohol.2013.08.005

Wang, L., Liu, X., Luo, X., Zeng, M., Zuo, L., and Wang, K. S. (2013). Genetic variants in the fat mass- and obesity-associated (FTO) gene are associated with alcohol dependence. *J. Mol. Neurosci.* 51, 416-424. doi: 10.1007/s12031-013-0044-2

Wang, X., Mozhui, K., Li, Z., Mulligan, M. K., Ingels, J. F., Zhou, X., et al. (2012). A promoter polymorphism in the Per3 gene is associated with alcohol and stress response. *Transl. Psychiatry* 2:e73. doi: 10.1038/tp.2011.71

Yang, J., Zaitlen, N. A., Goddard, M. E., Visscher, P. M., and Price, A. L. (2014). Advantages and pitfalls in the application of mixed-model association methods. *Nat. Genet.* 46, 100-106. doi: 10.1038/ng.2876

Young, A. I., Wauthier, F., and Donnelly, P. (2016). Multiple novel gene-by-environment interactions modify the effect of FTO variants on body mass index. *Nat. Commun.* 7:12724. doi: 10.1038/ncomms12724

Zhou, X., and Stephens, M. (2012). Genome-wide efficient mixed-model analysis for association studies. *Nat. Genet.* 44, 821-824. doi: 10.1038/ng.2310

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2018 Mulligan, Zhao, Dickerson, Arends, Prins, Cavigelli, Terenina, Mormede, Lu and Jones. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.