Excessive alcohol consumption has a vast societal impact with far-reaching consequences for drinkers and the people around them. In the United States, the cost of alcohol misuse from lost workplace productivity, health-care expenses, criminal-justice expense, and motor-vehicle crashes was estimated to be $249.0 billion (Sacks, Gonzales, Bouchery, Tomedi, & Brewer, 2015). Due to the vast cost associated with alcohol abuse, a variety of approaches have been developed to identify the biological underpinnings that drive alcohol consumption.

One such approach in preclinical research is the development of selectively bred rat lines that voluntarily demonstrate excessive alcohol consumption (Murphy et al., 2002). Examples of such lines include alcohol preferring (P) rats and high-alcohol-drinking (HAD-1) rats, as well as their counterparts the alcohol nonpreferring (NP) rats and low-alcohol-drinking (LAD-1) rats. These lines are widely used by laboratories around the world and fit many of the criteria for an animal model of alcohol use disorders (AUDs; Bell et al., 2016; Rodd, Bell, Sable, Murphy, & McBride, 2004). Research on these lines has led to the identification of therapeutic targets that are also effective in nonselective outbred strains, such as the Long-Evans and Wistar rats (Bell et al., 2016; Bell et al., 2012; McBride, Rodd, Bell, Lumeng, & Li, 2014).

Converging lines of evidence from the selectively bred rat lines and human literature suggest that genetic variation in dopaminergic and cholinergic neurotransmitter systems is associated with alcohol use and abuse (Enoch & Goldman, 2001; McBride, Kimpel, McClintick, Ding, Hyytia, et al., 2013b; Morzorati & Marunde, 2006; Stewart & Li, 1997; Tawa, Hall, & Lohoff, 2016). One approach to studying the role of cholinergic and dopaminergic neurotransmission in AUDs is to use biomarkers that are associated with activity in these neurotransmitter systems.

Ultrasonic vocalizations (USVs) are an established, noninvasive, real-time, functional biomarker of dopaminergic and cholinergic transmission in the rodent limbic system (Brudzynski, 2015). USVs are separated into two common subtypes: 50–55-kHz frequency modulated (FM) USVs, which are emitted in response to rewarding stimuli and reflect positive affective states (Ahrens, Ma, Maier, Duvauchelle, & Schallert, 2009; Avvisati et al., 2016; Brenes et al., 2016; Buck, Vendruscolo, Koob, & George, 2014; Burgdorf & Panksepp, 2006; Ciucci et al., 2009; Maier, Abdalla, Ahrens, Schallert, & Duvauchelle, 2012; Popik et al., 2014; Roccaro-Waldmeyer, Babalian, Müller, & Celio, 2016; Williams & Undieh, 2016) and 22–28-kHz USVs, which are emitted in response to aversive stimuli and reflect negative affect (Bardin et al., 2010; Berger, Williams, McGinnis, & Walker, 2013; Brudzynski, Bihari, Ociepa, & Fu, 1993; Chen, Li, Li, & Kirouac, 2012; Furlanetti, Coenen, Aranda, & Döbrössy, 2015; Lindquist, Jarrard, & Brown, 2004; Parsana, Moran, & Brown, 2012; Wang, Liang, Burgdorf, Wess, & Yeomans, 2008).

As real-time, functional biomarkers of activity in the mesolimbic dopaminergic and cholinergic neurotransmitter systems (Brudzynski, 2015; Wöhr & Schwarting, 2013), USVs can serve as important tools for understanding the role of these systems in AUDs as well as a number of other neurological disorders (Basken, Connor, & Ciucci, 2012; Brenes et al., 2016; Ciucci et al., 2009; Machold, 2013; Pan et al., 2014; Scattoni, Crawley, & Ricceri, 2009; Steele, Whitehouse, Aday, & Prus, 2017). However, due to the time-intensive and labor-intensive nature of USV analyses, this methodology is vastly underutilized. In fact, nearly all published USV studies report only short durations of USV activity (e.g., less than 30 minutes) in counts alone (Armstrong et al., 2001; Brudzynski, 2015; Burgdorf, Wood, Kroes, Moskal, & Panksepp, 2007; Johnson, Grant, Schallert, & Ciucci, 2015; Kaniuga et al., 2016; Portfors, 2007; Scattoni et al., 2009) but do not provide a full USV profile that could be revealed by repeated, long-duration assessment periods. Our development of WAAVES, a highly accurate algorithm that automates the analysis of USVs (Reno & Duvauchelle, 2014; Reno, Marker, Cormack, Schallert, & Duvauchelle, 2013) has allowed us to analyze data from experiments with relatively long time lines, something that was previously not feasible using manual techniques. In addition, WAAVES also allows an examination of each USV in much greater detail. Specifically, in addition to counting the number of USVs, the WAAVES analyses quantify several acoustic characteristics associated with each USV. These characteristics include the USV mean frequency in kilohertz, duration in milliseconds, bandwidth in kilohertz, and relative power in decibels. The result is a rich, multidimensional representation of each unique USV. These characteristics can be combined into large data sets for analyses using advanced statistical models and predictive analytics to examine the intricate relationships between rat USV characteristics and behavioral phenotypes, such as alcohol seeking and consumption. This is the approach taken in some of our recent work (Mittal et al., 2017; Reno et al., 2017) and the current report.

Overview of the present study

The overriding aim of the present work is to provide the first empirical test of the hypothesis that rodent ultrasonic vocalizations emitted by alcohol-naïve rats can be used to predict their future alcohol consumption. To test this hypothesis, we conducted a new study using Long-Evans rats and merged these results with our published data obtained from selectively bred high-alcohol (P and HAD-1) and low-alcohol (NP and LAD-1) drinking lines. The result is an integrative analysis of USV data collected from five different rat lines that highlight a novel and powerful role of USVs as direct predictors of excessive alcohol consumption and biomarkers of neurotransmission in these rat lines.

Instead of focusing on a single strain of rat, in the present study, we take the more challenging, but more broadly informative approach of testing our hypothesis across a number of rat types including selectively bred lines and a classically bred strain. With this goal in mind, we conducted a new experiment in Long Evans (LE) rats using the same procedures used in two previous studies in our lab, one conducted using P versus NP rats (Reno et al., 2017) and the other using HAD-1 versus LAD-1 rats (Mittal et al., 2017).

The rest of this manuscript is organized as follows. In the next section, we introduce the new study using Long-Evans rats. This is followed by a description of the predictive analytic approach and results from the combined sample that includes the new Long-Evans data and the P/NP (Reno et al., 2017) and HAD-1/LAD-1 (Mittal et al., 2017) studies. We conclude with a summary and integration of these results with the existing literature.

Method

Animals

Sixteen male Long-Evans rats were ordered from Harlan (Harlan Laboratories, Indianapolis, Indiana) at approximately 4 weeks of age. The average weight of the rats upon arrival was ~180 grams. Rats were handled 20 minutes per day, 5 days per week for 3 weeks to habituate them to the experimenters. The animals were housed in a 12-hour reverse-light cycle with lights out from 10:00 a.m. to 10:00 p.m. All USV recording sessions and the subsequent EtOH availability sessions began at the beginning of the dark cycle at 10:00 a.m. Rats had ad libitum access to food and water throughout all experiments. Data from the previously published studies consists of six HAD-1, six LAD-1 (Mittal et al., 2017), 10 P, and four NP rats (Reno et al., 2017). All experiments were conducted in accordance with Institutional Animal Care and Use Committee guidelines.

Ultrasonic vocalization recordings

USVs were recorded from alcohol-naïve rats in 4-hour sessions for 3 days/week (e.g., Monday, Wednesday, Friday) for 2 weeks. CM16 microphones were used with an UltraSound Gate interface (Avisoft Bioacoustics) to record USVs at a 250-kHz sampling rate and a 16-bit resolution. The recording sessions started at the beginning of the dark cycle. The animals were weighed and transported to the recording room where they were transferred into recording cages, which were identical to the home cages but used only for the recording sessions. Since novelty or conspecific scents could affect rodent behaviors (Wöhr, Houx, Schwarting, & Spruijt, 2008), each animal was assigned its own recording cage. The animals were given access to food and water in the cage. Based on rat and chamber size, we approximate the distance between the animal’s head and the centered microphone to range from 5 cm to 28.4 cm. Following the recording sessions, the animals were transported back to the vivarium.

Analysis of USVs

USV recordings were analyzed using the WAAVES program developed in our laboratory (Reno & Duvauchelle, 2014; Reno et al., 2013). This program reads USV containing audio files and produces a frequency spectrogram, which is subsequently analyzed using MATLAB’s Image Processing Toolbox (MathWorks, Inc., Natick, MA). The program utilizes a variety of filters to separate USV calls from background noise, and several measurements of interest are quantified and extracted from each call. For 50–55-kHz FM USVs, WAAVES identifies sound objects with a minimum duration of 5 ms occurring in a range of 30 to 120 kHz. An intercall interval of 10 ms was used to discriminate between individual calls and avoid counting call fragments as separate calls. FM USVs were defined as calls that varied more than 5 kHz over the entire duration of the call, and 22–28-kHz calls were identified as sound objects occurring in a frequency range of 20 to 30 kHz with a minimum duration of 200 ms. An intercall interval of 100 ms was used to separate individual calls. Once the calls are identified, several measurements of interest are extracted from each USV call and stored for subsequent analysis. In the present study, the mean frequency, duration, bandwidth, and power (e.g., in decibel measurements relative to sound pressure at human detection levels) for 50–55-kHz FM and 22–28-kHz USV calls were extracted for statistical analysis.

Validation process for WAAVES automation

Validation of WAAVES-generated USV data requires correspondence with human-derived analyses. Experimenters manually analyzed subsets of USV data recorded during the experiment to compare human assessment with WAAVES output. A subset of 10-min USV files recorded from LE rats (22–28-kHz calls: 40 files; 50–55-kHz calls: 40 files) were used for manual validation. The total number of calls identified via manual analysis was correlated with the total number of calls identified by the automated WAAVES program. The correlation coefficients are reported in the results.

Ethanol availability sessions

Following the USV recording rats were given access to 24-hour chronic intermittent ethanol sessions 3 days/week for 4 weeks as previously described (Mittal et al., 2017; Reno et al., 2017). The sessions began at the beginning of the dark cycle. The rats were weighed at the beginning of each session and received three bottles containing either water, 15% EtOH, or 30% EtOH. The bottles were removed 24 hours later and weighed to measure ethanol or water consumption. Any potential loss of fluids due to spillage were insignificant, and therefore not accounted for in the present study. The dose of EtOH in g/kg was measured using the following formula:

$$ \mathsf{EtOH}\kern0.5em \mathsf{Dose}\kern0.5em \left(\mathsf{g}/\mathsf{kg}\right)={\left(\left({\mathsf{0.3}}^{\ast}\kern0.5em \mathsf{Change}\kern0.5em \mathsf{in}\kern0.5em \mathsf{weight}\kern0.5em \mathsf{of}\kern0.5em \mathsf{30}\%\kern0.5em \mathsf{bottle}\right)+\left({\mathsf{0.15}}^{\ast}\kern0.5em \mathsf{Change}\kern0.5em \mathsf{in}\kern0.5em \mathsf{weight}\kern0.5em \mathsf{of}\kern0.5em \mathsf{15}\%\kern0.5em \mathsf{bottle}\right)\right)}^{\ast}\kern0.5em {\mathsf{0.8}}^{\ast}\kern0.5em \left(\mathsf{1}/\left({\mathsf{0.001}}^{\ast}\kern0.5em \mathsf{Weight}\kern0.5em \mathsf{of}\kern0.5em \mathsf{rat}\kern0.5em \mathsf{in}\kern0.5em \mathsf{grams}\right)\right). $$

Advanced statistical methods for testing the USV–alcohol consumption relationship

Traditional USV studies have focused primarily on differences in total call counts in the 22–28-kHz and 50–55-kHz FM categories, across experimental groups. This approach has proven to be a reliable index of underlying neurotransmission and is sensitive to experimental and pharmacological manipulations (Brudzynski, 2013; Ciucci et al., 2009; Knutson, Burgdorf, & Panksepp, 2002; Pan et al., 2014). Using this approach our lab has shown that USVs are responsive to a wide range of drugs of abuse and can serve as a useful model for investigating the neural basis of drug reward and the development of drug craving, behavioral sensitization and drug reward tolerance (Ahrens et al., 2009; Ahrens et al., 2013; Ma, Maier, Ahrens, Schallert, & Duvauchelle, 2010; Maier et al., 2012; Reno et al., 2015; Thakore et al., 2016). Many of these studies, including those from our lab, have used standard parametric statistical approaches (e.g., repeated-measures or two-way ANOVAs) to analyze the USV count data. However, these methods are not sufficiently powerful to analyze the USV acoustic characteristic data captured with WAAVES.

Predictive analytic approaches, such as linear discriminant analysis (LDA), are required when the goal is to determine whether a linear combination of the USV acoustic characteristics can successfully discriminate between rat lines (e.g., P vs. NP or HAD-1 vs. LAD-1 rats). Using LDA on the USV data produces a linear equation that takes the following form:

$$ \mathsf{LDA}\kern0.5em \mathsf{score}={\beta_{\mathsf{Mean}\kern0.3em \mathsf{Frequency}}}^{\ast}\kern0.5em \mathsf{Mean}\kern0.5em \mathsf{Frequency}+{\beta_{\mathsf{Duration}}}^{\ast}\kern0.5em \mathsf{Duration}+{\beta_{\mathsf{Bandwidth}}}^{\ast}\kern0.5em \mathsf{Bandwidth}+{\beta_{\mathsf{Power}}}^{\ast}\kern0.5em \mathsf{Power}. $$

This method computes the β coefficients that can be used to combine the data from all four acoustic characteristics into a singular LDA score that maximizes the percent of USV calls classified into the appropriate experimental classes (e.g., P vs. NP or HAD-1 vs. LAD-1 rats). For example, if we have a dataset consisting of 22–28-kHz calls made by P and NP rats, then the LDA method will produce β coefficients that maximize the difference between the LDA scores of the P and the NP rats respectively, such that we can take the average LDA scores of each group and construct a decision threshold halfway between the two averages to make predictions about unclassified data. That is, if the LDA score (calculated using the equation above) of an unknown animal falls below the decision threshold, then we predict that the rat is a P rat and if it is above the threshold then we predict that the rat is an NP rat (or vice-versa). An additional advantage of this method is that we are able to reduce the four USV characteristics into a one-dimensional LDA score, which can be easily plotted. The LDA method can be implemented using the freely available R package “MASS” (Venables & Ripley, 2002).

Results

Using LDA to accurately separate selectively-bred high-drinking and low-drinking rat lines

Our lab has successfully used linear discriminant analysis to classify HAD-1 and LAD-1 rats with 100% accuracy solely on the basis of the 50–55-kHz FM USV data (Mittal et al., 2017). Similarly, we have also been able to use 22–28-kHz USV data to discriminate alcohol-naïve HAD-1 rats from LAD-1 rats with a classification accuracy of 81.96% (Mittal et al., 2017) and alcohol-naïve P rats from NP rats with 100% accuracy (Reno et al., 2017). Our work on the differences in USV characteristics between the HAD-1/LAD-1 and the P/NP lines suggests that traits associated with alcohol consumption in these lines may also affect USV production. Consequently, USVs may serve as a potential biomarker for excessive alcohol consumption. However, it is also possible that the differences in USV characteristics between these lines may simply have been a byproduct of genetic drift stemming from potential inbreeding within each selectively-bred line (Crabbe, Phillips, Kosobud, & Belknap, 1990). The P and NP rats were derived from a closed colony of Wistar rats, while the HAD-1 and LAD-1 rats were developed using the N/Nih heterogenous stock (Bell et al., 2012; Gongwer, Murphy, McBride, Lumeng, & Li, 1989). Therefore, it is unlikely that the similar differences in USV characteristics observed across these different high and low alcohol consuming lines, derived from distinct breeder stocks, are simply a result of random genetic drift during the selective-breeding process.

Using LDA to predict alcohol consumption in a combined HAD-1, LE, LAD-1 sample

Long-Evans rats, derived from mating a Wistar female with a wild grey rat (Weisbroth, 1969), are commonly used by many labs to study alcohol use and abuse. Although derived from a different breeder stock, HAD-1 and LAD-1 rats show stratification in alcohol consumption levels. Therefore, we first examined alcohol consumption in the combined sample as shown in Fig. 1a. We saw that while the HAD-1 rats consume high levels of alcohol and the LAD-1 rats consumed low levels of alcohol, the LE rats show moderate levels of consumption ranging from HAD-1 levels to LAD-1 levels (see Fig. 1a). Second, we assessed whether similar stratification could be observed in the USV data amongst the HAD-1, LAD-1, and LE rat lines. All of the USVs emitted during the 2-week period were used for this analysis without any references to recording days. Using LDA, we estimated the linear discriminant scores for the three rat lines. When the LDA method is used to discriminate between two classes of data a single linear equation is produced to generate a decision threshold that separates the two classes. However, when dealing with three classes two decision thresholds are needed. Therefore, the LDA method produced two linear discriminant scores for each of the 50–55-kHz FM (see Fig. 2a) and 22–28-kHz USV (see Fig. 2b) categories.

Fig. 1
figure 1

Daily alcohol consumption of HAD-1–LE–LAD-1 and P–LE–NP rats across 12 drinking sessions. Individual symbols represent individual rats on each of the 12 EtOH session, with a slight offset on the x-axis used to minimize overlap between symbols. a Daily alcohol consumption ranged from 2.45 to 17.53 g/kg for HAD-1 rats, from 0.55 to 12.00 g/kg for LE rats, and from 0.28 to 5.34 g/kg for LAD-1 rats. b Daily alcohol consumption ranged from 1.06 to 11.86 g/kg for P rats, from 0.55 to 12.00 g/kg for LE rats, and from 0.36 to 6.39 g/kg for NP rats

Fig. 2
figure 2

Maximal separation between HAD-1–LE–LAD-1 and P–LE–NP rats achieved via linear discriminant analysis using 22–28-kHz and 50–55-kHz FM USV data. a Linear discriminant analysis provided 90.91% separation accuracy between HAD-1 and LE rats, and 100% separation accuracy between HAD-1 and LAD-1 rats based on 50–55-kHz FM USV data. No clear separation was observed between LE and LAD-1 rats. b Linear discriminant analysis based on 22–28-kHz USV data showed 100% separation between HAD-1 and LE rats, 95.45% separation between LAD-1 and LE rats, and 100% separation accuracy between HAD-1 and LAD-1 rats. c Linear discriminant analysis based on 50–55-kHz FM USV data did not show any clear separation between the P, Len, and NP rat lines. d Linear discriminant analysis based on 22–28-kHz USV data showed 92.31% separation between P and LE rats, 95.00% separation between NP and LE rats, and 100% separation accuracy between P and NP rats

We found that the linear discriminant scores from 50–55-kHz FM USVs showed some degree of separation between the HAD-1 and LE rats (90.91% separation), and the HAD-1 and LAD-1 rats (100% separation; see Fig. 2a), though no clear distinction existed between the LDA scores from the LE and LAD-1 rats (54.54% separation). On the other hand, the LDA scores from the 22–28-kHz USV data showed clear clustering across the HAD-1, LAD-1, and LE rat lines (HAD-1 vs. LE separation: 100%; LAD-1 vs. LE separation: 95.45%; HAD-1 vs. LAD-1 separation: 100%; see Fig. 2b). Because 50–55-kHz FM and 22–28-kHz USV emissions are mediated through the mesolimbic dopaminergic and cholinergic neural pathways (Brudzynski, 2015; Wöhr & Schwarting, 2013), these results suggest that basal differences in these transmitter systems may exist between the selectively bred HAD/LAD and the outbred LE rat lines. However, as no study to date has directly explored these differences across these three rat lines, this hypothesis requires further work for serious consideration.

Linear discriminant scores predict future alcohol consumption

Equipped with the linear discriminant scores based on the USV data from alcohol-naïve rats and the future alcohol consumption data, we tested the hypothesis that alcohol consumption levels in each rat can be predicted from their USV acoustic characteristics. To test this hypothesis, we generated a multivariate regression model that used the four linear discriminant values described above (i.e., from Fig. 2a–b) to retrospectively predict total alcohol consumption for each rat during the subsequent alcohol availability sessions. Multivariate linear regression modeling strongly supported our hypothesis and revealed that the linear discriminant scores could be used to predict future alcohol consumption in these rat lines (multiple r 2 = .7532), F(4, 23) = 17.55, p < .0001. Critically, to ensure that this linear model did not overfit the data, we used the leave-one-out cross-validation method. This is a commonly used validation technique that yields a conservative estimate of the accuracy of the regression model (Kearns & Ron, 1999). Using this method, for every rat, a model is trained excluding that rat and the model is then used to make predictions on the data from the excluded rat. Thus, the model predicted values are generated on unseen data that are not used in building the model. As shown in Fig. 3a, the model predicted consumption values showed a high degree of correlation with the total alcohol consumption of the HAD-1, LAD-1, and LE rats during the 4-week alcohol access period (r 2 = .5999), F(1, 26) = 38.99, p < .0001 (se Fig. 3a). This finding provides the first direct evidence that USV acoustic characteristics can predict future alcohol consumption in alcohol-naïve rats. Importantly, this finding is not focused on a single strain of rats, but rather holds across three strains that are commonly used in alcohol research. The breadth of this finding across strains attests to its generality.

Fig. 3
figure 3

Comparison of total alcohol consumption and model predicted alcohol consumption of HAD-1–LE–LAD-1 and P–LE–NP rats during 4 weeks of 24-hour alcohol availability sessions. a Multivariate linear regression models developed using linear discriminant scores from HAD-1, LE, and LAD-1 rats showed a high degree of correlation between the model predicted alcohol consumption and the actual consumption of these rats (r 2 = .5999), F(1, 26) = 38.99, p < .0001. b Multivariate linear regression models developed using linear discriminant scores from P, LE, and NP rats also showed a high degree of correlation between the model predicted alcohol consumption and the actual consumption of these rats (r 2 = .5112), F(1, 28) = 29.29, p < .0001

One might argue that by including selectively bred high-alcohol and low-alcohol consuming lines in the multivariate regression we increased the likelihood of finding that USVs predict future alcohol consumption. This reasoning is problematic because it assumes that USV profiles differ across these lines in a way that is systematically related to alcohol consumption. This is an empirical question that is tested in the current report. Even so, we conducted a secondary analysis focused just on the LE rats. In this LE subsample, we continue to find that alcohol consumption is strongly predicted from the USV profiles, F(1, 14) = 12.90, p < .01.

Using LDA to predict alcohol consumption in a combined P, LE, NP sample

Given the shared ancestry between the LE and the P/NP rat lines, we also sought to investigate potential similarities or differences in USV acoustic characteristics among these three lines. As shown in Fig. 1b, the alcohol consumption of LE rats also ranges from NP rat consumption levels to P rat consumption levels. Therefore, we also tested whether linear discriminant scores could be used to accurately separate the LE rats from the P and the NP rats. As seen in Fig. 2a, the linear discriminant scores from the 50–55-kHz FM USV data did not show any clear distinction between the three rat lines (P vs. LE separation: 80.77%; NP vs. LE separation: 65.00%; P vs. NP separation 71.43%). However, the 22–28-kHz USV associated LDA values showed clear clustering across the P, NP, and LE rats (P vs. LE separation: 92.31%; NP vs. LE separation: 95.00%; P vs. NP separation 100%; see Fig. 2b). Since 22–28-kHz USVs are associated with tegmental cholinergic transmission, these results suggest that there may be differences in basal cholinergic activity in the limbic regions across these three lines. No study has directly compared these neurotransmitter systems between P and LE rats. Together with the HAD/LAD results, these findings highlight a need for future studies to directly explore these systems and compare neurotransmitter activity across these rat lines.

Linear discriminant scores predict future alcohol consumption

Next, we tested whether the linear discriminant scores based on the USV data from alcohol-naïve rats could also predict alcohol consumption in these rats.Footnote 1 Multivariate linear regression models showed that these LDA score could also be used to predict alcohol consumption in these rats (multiple r 2 = .6995), F(4, 25) = 14.55, p < .0001. The model predicted consumption values, generated using the leave-one-out cross-validation method described above, showed a significant correlation with the alcohol consumption of the P, NP, and LE rats (r 2 = .5112), F(1, 28) = 29.29, p < .0001. Figure 3b provides a representation of this relationship using the assigned and predicted alcohol consumption values from one of the 10,000 simulations. These findings complement the HAD-LAD-LE data described above and provide strong evidence that USV data from selectively bred lines can be generalized to Long-Evans rats. Together, these results provide clear and direct evidence that USV data can be used to predict alcohol consumption in rats.

As with the HAD/LAD analysis, one might argue that by including selectively bred high-alcohol and low-alcohol consuming lines in the multivariate regression, we increased the likelihood of finding that USVs predict future alcohol consumption. To address this potential concern, we conducted a secondary analysis focused just on the LE rats. In this LE subsample, we continue to find that alcohol consumption is strongly predicted from the USV profiles, F(1, 14) = 30.77, p < .0001.

Validation of USV analysis

A subset of USV files were manually analyzed using Avisoft SASLab Lite (Avisoft Bioacoustics) to validate WAAVES output. WAAVES-automated analysis and manual analysis results were highly correlated for both the 22–28-kHz (r = .996, p < .0001) and 50–55-kHz calls (r = .970, p < .0001).

Discussion

The present study provided the first known test of the hypothesis that rodent ultrasonic vocalizations obtained while animals are alcohol naïve serve as predictors for future alcohol consumption. Instead of focusing on a single rat strain, we took the more challenging, but more general, approach of testing this hypothesis across five rat lines: P, NP, HAD-1, LAD-1, and Long-Evans. We report the results from a new study with Long-Evans rats that uses the same procedure used in two other studies conducted in our lab: one with P and NP rats (Reno et al., 2017), and the other with HAD-1 and LAD-1 rats (Mittal et al., 2017). In all three studies, we recorded spontaneously emitted USVs while the rats were alcohol naïve and followed this with free alcohol access sessions in which alcohol consumption was measured. We combined these three data sets and utilized predictive analytic techniques to test this hypothesis.

Here, we show that differences in USV expression, which are thought to be mediated through mesolimbic dopaminergic and cholinergic neurotransmitter systems, are linked to the propensity for alcohol consumption in selectively bred and nonselectively bred rat lines. We showed that USVs could serve as real-time functional markers of activity in these systems. Two major findings emerged. First, we showed that linear discriminant scores calculated using the USV data could be used to discriminate Long-Evans rats from rat lines selectively bred for high or low levels of alcohol consumption. These results point to potential differences in neurotransmission across these five rat lines and highlight the need for direct investigation of these systems. Second, using multivariate linear regression, we found that future alcohol consumption could be predicted from USV profiles obtained while animals were alcohol naïve. Thus, these scores provide information about basal neurotransmission in these rat lines, unencumbered by any potential disruptions from alcohol exposure and highlight a novel and potentially invaluable use of USV data as a predictor of alcohol consumption in these five rat lines.

Both pairs of rat lines tested in this study are known to have basal differences in dopaminergic and cholinergic transmission. For instance, P and HAD-1 rats have lower tissue levels of DA and its metabolites (DOPAC and HVA) in the nucleus accumbens and the anterior striatum as compared to their NP and LAD-1 counterparts (Gongwer et al., 1989; Murphy, McBride, Lumeng, & Li, 1982, 1987). This may be due to a potential deficit in the VTA DA projections in these rat lines. Indeed, P rats have reduced DA projections from the VTA to the nucleus accumbens (Zhou, Zhang, Lumeng, & Li, 1995) as well as lower levels of DA D2 receptor binding sites (McBride, Chernet, Dyr, Lumeng, & Li, 1993), as compared to the NP rats. In order to compensate for the lower levels of DA, P rats show increased burst firing of VTA DA neurons, though a similar increase was not seen in the HAD rats (Morzorati, 1998; Morzorati & Marunde, 2006). Alcohol consumption results in reduced D2 auto-receptor function (Engleman, McBride, Li, Lumeng, & Murphy, 2003; Engleman et al., 2000), and enhanced DA efflux (Thielen et al., 2004; Weiss, Lorang, Bloom, & Koob, 1993) and reuptake (Sahr, Thielen, Lumeng, Li, & McBride, 2004).

The potential differences in cholinergic transmission between the selectively bred rat lines remain largely unexplored. However, the few studies that have been conducted suggest that the accumbal cholinergic system may be more active in P rats than NP rats (Bell et al., 2016). For instance, P rats have lower expression levels of striatal nicotinic acetylcholine receptor than NP rats (Tizabi, Getachew, Davila-Garcia, & Taylor, 2001). Studies have also identified higher expression levels of genes encoding choline acetyltransferase, muscarininc acetylcholine receptor 3, and channels responsible for the uptake and vesicular transport of acetylcholine during synthesis in the nucleus accumbens shell of P rats as compared to the NP rats (McBride, Kimpel, McClintick, Ding, Hauser, et al., 2013a; McBride, Kimpel, McClintick, Ding, Hyytia, et al., 2013b). Furthermore, administration of nicotinic antagonists or desensitizing agents can reduce alcohol consumption in P rats (Rezvani et al., 2010; Srisontiyakul, Kastman, Krstew, Govitrapong, & Lawrence, 2016). Although we did not directly explore these systems in our study, these studies provide support to the notion that the USV differences across these lines are a reflection of underlying differences in these neurotransmitter systems.

Although USV recordings are noninvasive and USV recording experiments are relatively easy to set up, USVs have remained relatively underutilized as important preclinical tools for the investigation of neural substrates underlying alcohol and drug abuse, primarily due to the manually intensive nature of the USV analysis process. The development of WAAVES (Reno & Duvauchelle, 2014; Reno et al., 2013), and other automated USV detection methods (Barker, Herrera, & West, 2014; Barker & Johnson, 2017) can help these problems and should encourage more labs to integrate USV recordings to their respective batteries of behavioral assays. Widespread use of USVs can unlock the true potential of this behavioral metric as a noninvasive correlate of underlying neural activity.

Conclusion

This present study provides the first known test of the hypothesis that rodent ultrasonic vocalizations obtained while animals are alcohol naïve serve as a biomarker for and thus predict future alcohol consumption. We find support for this hypothesis in a combined sample of five rat types commonly used in alcohol research. These results highlight a novel role of USVs as biomarkers of neurotransmission and the propensity for high levels of alcohol consumption in rodents.