Abstract
Cognitive impairments, uncontrolled drinking, and neuropathological cortical changes characterize alcohol use disorder. Dysfunction of the orbitofrontal cortex (OFC), a critical cortical subregion that controls learning, decision-making, and prediction of reward outcomes, contributes to executive cognitive function deficits in alcoholic individuals. Electrophysiological and quantitative synaptomics techniques were used to test the hypothesis that heavy drinking produces neuroadaptations in the macaque OFC. Integrative bioinformatics and reverse genetic approaches were used to identify and validate synaptic proteins with novel links to heavy drinking in BXD mice. In drinking monkeys, evoked firing of OFC pyramidal neurons was reduced, whereas the amplitude and frequency of postsynaptic currents were enhanced compared with controls. Bath application of alcohol reduced evoked firing in neurons from control monkeys, but not drinking monkeys. Profiling of the OFC synaptome identified alcohol-sensitive proteins that control glutamate release (e.g., SV2A, synaptogyrin-1) and postsynaptic signaling (e.g., GluA1, PRRT2) with no changes in synaptic GABAergic proteins. Western blot analysis confirmed the increase in GluA1 expression in drinking monkeys. An exploratory analysis of the OFC synaptome found cross-species genetic links to alcohol intake in discrete proteins (e.g., C2CD2L, DIRAS2) that discriminated between low- and heavy-drinking monkeys. Validation studies revealed that BXD mouse strains with the D allele at the C2cd2l interval drank less alcohol than B allele strains. Thus, by profiling of the OFC synaptome, we identified changes in proteins controlling glutamate release and postsynaptic signaling and discovered several proteins related to heavy drinking that have potential as novel targets for treating alcohol use disorder.
SIGNIFICANCE STATEMENT Clinical research identified cognitive deficits in alcoholic individuals as a risk factor for relapse, and alcoholic individuals display deficits on cognitive tasks that are dependent upon the orbitofrontal cortex (OFC). To identify neurobiological mechanisms that underpin OFC dysfunction, this study used electrophysiology and integrative synaptomics in a translational nonhuman primate model of heavy alcohol consumption. We found adaptations in synaptic proteins that control glutamatergic signaling in chronically drinking monkeys. Our functional genomic exploratory analyses identified proteins with genetic links to alcohol and cocaine intake across mice, monkeys, and humans. Future work is necessary to determine whether targeting these novel targets reduces excessive and harmful levels of alcohol drinking.
Introduction
Alcohol abuse disorder (AUD) is a detrimental worldwide public health problem with increased prevalence over the past 20 years (Whiteford et al., 2013). Accordingly, alcohol ranks higher than 20 abused substances, including heroin, cocaine, and methamphetamine, as the most harmful drug (Nutt et al., 2010). AUD is a chronic relapsing disease characterized by functional and metabolic tolerance, a distinctive withdrawal syndrome, recurrent and heavier alcohol use, and cognitive impairments. Importantly, compelling evidence from studies that have examined the neurobiology of cognitive deficits in alcoholic individuals has demonstrated that poor cognitive function is a risk factor for increased vulnerability to relapse (Charlet et al., 2014), suggesting that improving cognitive performance during abstinence may prevent relapse. However, current pharmacotherapies for treating AUD do not target alcohol-induced cognitive impairments.
Subregions of the prefrontal cortex (PFC) are key regulators of executive functions that are often critically impaired following heavy alcohol use. Among these, the orbitofrontal cortex (OFC) contributes to outcome-based decision-making related to the expected rewards and the pleasurable properties of rewarding stimuli (Stalnaker et al., 2015). Imaging and behavioral studies show that patients with AUD have greater activation of the OFC and reduced performance when completing OFC-dependent tasks compared with healthy control subjects, demonstrating that prolonged drinking can impair OFC function (Fortier et al., 2008, 2009). Consistent with deficits reported in individuals with AUD, rodent models of alcohol dependence induce OFC-dependent behavioral deficits and adaptations in synaptic physiology and intrinsic excitability in OFC neurons (Crews and Boettiger, 2009; Badanich et al., 2011; McGuier et al., 2015; Nimitvilai et al., 2016). Although rodent studies on the effects of alcohol on cortical function are valuable, they are limited by the reduced complexity of the rodent cortex present in primates (Carmichael and Price, 1994; Ongür and Price, 2000). To circumvent these limitations, a nonhuman primate alcohol drinking model was used that closely mimics the range of consumption patterns reported in humans and generates a proportion of monkeys self-administering alcohol in a chronic, heavy manner (Baker et al., 2014). Using this model, we combined comprehensive functional measurements of synaptic physiology with integrated synaptic profiling to identify neuroadaptations and biomarkers in the OFC of heavy-drinking monkeys. Because the nonhuman primate self-administration model is highly translational but is limited by small cohort sizes, additional drinking and integrative genomic experiments were performed using BXD recombinant inbred (RI) strains of mice to validate the novel targets identified from the monkey studies.
Materials and Methods
Alcohol-self administration.
Monkey care and all procedures were approved by the Oregon National Primate Research Center Animal Care and Use Committee at Oregon Health & Science University and performed according to the National Institutes of Health (NIH) Guidelines for the Care and Use of Mammals in Neuroscience and Behavioral Research. Two cohorts of adult male cynomolgus monkeys (n = 21; Macaca fascicularis; age, 67–84 months; weight, 6.02–9.12 kg; World Wide Primates) were used for these studies, and 16 of the monkeys were trained to self-administer alcohol using a schedule-induced polydipsia technique as described previously (Vivian et al., 2001; Grant et al., 2008). Following the 3-month induction phase, monkeys had open access to 4% alcohol (v/v) for 22 h/d for 6 months. Blood samples were obtained by femoral blood sampling at 7 h into the daily 22 h access period every 4–7 d. At the completion of the drinking paradigm, control and drinking monkeys were anesthetized at a time when they would normally begin a drinking session. Brains were rapidly removed following routine procedures (Davenport et al., 2014), and tissue blocks of area 13L were prepared for electrophysiology or protein analysis. Ten monkeys (2 controls and 8 drinkers; MATRR.org INIA Cohort 13) were used for functional assays, and 11 monkeys (3 controls and 8 drinkers; MATRR.org INIA Cohort 9) were used for quantitative proteomics and Western blot analyses. Because of the individual differences in average daily intake in this voluntary alcohol-drinking procedure (Baker et al., 2014, 2017) that may preclude identifying functional or proteomic neuroadaptations, alcohol-drinking monkeys were oversampled in these studies.
Preparation of brain slices.
Brain slices containing area 13L were prepared for whole-cell patch-clamp electrophysiology experiments with the experimenter blind to the treatment conditions. Following the necropsy procedure and rapid removal of the brain, the tissue was blocked coronally for the frontal cortex and mounted in a Vibroslicer (Leica) containing ice-cold oxygenated (95% O2, 5% CO2) sucrose containing buffer, and coronal sections (300 μm) were cut. Slices containing area 13L were immediately placed in a holding chamber containing oxygenated artificial CSF (aCSF) at 34°C for 45 min and kept at room temperature for at least 30 min before transferring to the recording chamber. The glutamate NMDA receptor blocker AP5 (50 μm) was added to the aCSF during incubation to prolong the viability of neurons, and slices were washed with regular oxygenated aCSF for at least 20 min before the whole-cell recordings were performed. Although AP5 is reported to wash out of acute slice preparations (Dozmorov et al., 2004; Herman et al., 2011), application of AP5 for ≥12 h enhances evoked firing (Ishikawa et al., 2009; Lee and Chung, 2014), suggesting that there may be a compensatory increase in intrinsic excitability under these recording conditions. However, all groups were treated identically, so this is unlikely to explain a difference in firing between treatment groups.
The composition of the sucrose-containing cutting solution used was as follows (in mm): 194 sucrose, 30 NaCl, 4.5 KCl, 1.2 NaH2PO4, 1 MgCl2, 10 glucose, and 26 NaHCO3, adjusted to 305–315 mOsm. The composition of the aCSF was as follows (in mm): 125 NaCl, 2.5 KCl, 1.25 NaH2PO4, 1.3 MgCl2, 2.0 CaCl2, 0.4 ascorbate, 10 glucose, and 25 NaHCO3, adjusted to 290–310 mOsm. Both solutions were saturated with 95% O2/5% CO2, pH 7.4.
Whole-cell patch-clamp electrophysiology.
An individual slice was placed in the recording chamber fixed to the stage of an upright microscope (Axio scope examiner D1, Zeiss) and perfused with 29–31°C oxygenated aCSF. Recordings were localized to area 13L of the OFC under infrared optics using a 40× water-immersion objective. Thin-wall borosilicate glass electrodes (outer diameter = 1.5 mm; inner diameter = 1.17 mm) were pulled on a Sutter Instrument P97 Micropipette Puller and had tip resistances ranging from 1.9 to 5.5 MΩ. Patch pipettes filled with an internal solution were slowly lowered onto the layer V/VI pyramidal neurons to obtained a seal (>1 GΩ) followed by breakthrough to gain whole-cell access. All the whole-cell recordings were performed using an Axon MultiClamp 700B amplifier (Molecular Devices) and analyzed with pClamp software (Molecular Devices). Events were filtered at 2 kHz and digitized at a sampling rate of 10 kHz. To measure the intrinsic excitability of layer V/VI area 13L pyramidal neurons, current-clamp recordings were performed. Spike firing was induced by direct current injection (−40 to 240 pA) through patch pipettes filled with a potassium gluconate internal solution [composition as follows (in mm): 120 K-gluconate, 10 KCl, 10 HEPES, 2 MgCl2, 1 EGTA, 2 NaATP, and 0.3 NaGTP]. Recordings were analyzed off-line in AxographX (Axograph) for the number of spikes in response to each current step, resting membrane potential (in millivolts), action potential threshold (in millivolts), height (in millivolts), half-width (in milliseconds), rise time (in milliseconds), interspike interval (ISI; ms), frequency (in hertz), and after-hyperpolarization (AHP; in millivolts). The spike frequency adaptation ratio was calculated as first ISI/last ISI. To monitor spontaneous postsynaptic currents (sPSCs), OFC pyramidal neurons were voltage clamped at a membrane potential of −70 mV with patch pipettes filled with a cesium methanesulfonate internal solution [with composition as follows (in mm): 125 CsMeSO3, 10 CsCl, 5 NaCl, 10 HEPES, 1 EGTA, 2 MgCl2, 5 MgATP, and 0.3 NaGTP]. To maximize the chances of observing changes in synaptic activity in the limited number of OFC slices available, no blockers of synaptic activity were used in these recordings. However, under the recording conditions used, the estimated reversal potentials are biased in favor of detecting AMPA-mediated conductances (reversal potential, ∼0 mV) over those generated by chloride-permeable GABAA receptors (reversal potential, −51.3 mV). Differences in the number of spikes were assessed with two-way ANOVA, and differences in other electrophysiological characteristics of OFC neurons were analyzed with unpaired t test. Spontaneous events were analyzed off-line with AxographX software using a sliding template-matching algorithm. Differences in sPSC parameters were assessed with unpaired t test using GraphPad Prism (GraphPad Software). Comparisons were considered significantly different at p < 0.05, and the experimenter was blind to the treatment conditions for data acquisition of the electrophysiology experiments.
Preparation of iTRAQ-labeled peptides.
For proteomic analysis, Triton X-100 insoluble membrane fractions that are enriched in postsynaptic density (PSD) proteins (Uys et al., 2016) were prepared from 11 monkeys (3 controls and 8 drinkers) following our routine methods (Mulholland et al., 2011; Uys et al., 2016). Due to the limited number of isobaric tags, 8 (3 controls and 5 drinkers) of the samples were used for mass spectrometry (MS), while all 11 samples were used for Western blotting, as described below. Isobaric Tag for Relative and Absolute Quantitation (iTRAQ; AB SCIEX) uses covalent attachment of isobaric tags to individual samples, allowing for qualitative and quantitative analysis of eight protein samples in a single solution by liquid chromatography-tandem MS (LC-MS/MS; Ross et al., 2004). An aliquot of each sample was taken for determination of protein concentration by the bicinchoninic acid assay (Pierce Biotechnology). Each sample (50 μg of protein in equal volumes) was reduced in 5 mm tris-(2-carboxyethyl)phosphine for 1 h at 60°C, alkylated in 9 mm methyl methane-thiosulfonate for 30 min at room temperature, diluted to 0.1% in PPS Silent surfactant (Agilent Technologies) with 50 mm triethylammonium bicarbonate, and digested overnight in trypsin (1:10, enzyme/protein ratio; Promega Corporation) at 37°C. Tryptic peptides were labeled with iTRAQ 8-plex reagents for 2 h at room temperature, and the contents of all digested samples were combined after quenching the labeling reaction with 50 mm ammonium bicarbonate. Control monkeys were labeled with tags 113, 114, and 116, whereas monkeys with a history of alcohol self-administration were labeled with tags 115, 117, 118, 119, and 121.
Strong cationic exchange of labeled peptides and LC-MS/MS.
Lyophilized labeled peptides were solubilized, fractionated by strong cation exchange, and mass analyzed following our previously published methods (Uys et al., 2016). Briefly, eluting peptides were mass analyzed by data-dependent acquisition on the Orbitrap Elite mass spectrometer (Thermo Scientific) with Xcalibur 2.2 software (Thermo Scientific). The top seven most intense ions in the Fourier Transform Mass Spectrometry survey scan were selected for fragmentation by alternating higher-energy collisional dissociation (HCD) at a normalized collision energy of 40% and collision-induced dissociation (CID) at 35%. Full-scan MS spectra were acquired at a target value of 1 × 106 and a resolution of 60,000, and CID MS/MS spectra were acquired in the ion trap at a target value of 1 × 104; the HCD MS/MS spectra were recorded at a target value of 5 × 104 with a resolution of 15,000. Ions with a +1 charge were excluded from selection. Dynamic exclusion was enabled with a repeat count of 2, a duration of 30 s, an exclusion list size of 500, and an exclusion duration of 80 s.
Database searching, peptide identification, quantitation, and functional annotation.
A recent study demonstrated that the human Uniprot database is an excellent alternative strategy to replace deficient nonhuman primate databases for proteomic profiling of monkey tissue samples (Lee et al., 2015). Thus, the raw files were searched using Mascot and SequestHT in Proteome Discoverer 2.1 (Thermo Scientific) against a human database (Swiss-Prot; 20,186 reviewed sequences containing isoforms; downloaded in February 2016). Parameters for peptide identification were as follows: precursor mass tolerance of 10 ppm, fragment mass tolerance was set at 0.8 Da for the CID spectra and 0.1 Da for the HCD spectra, fully tryptic peptides with a maximum of two missed cleavages, static modifications of peptide N termini, and lysines with the iTRAQ reagent and cysteine alkylation with methylthio; methionine oxidation was included as a variable modification. The search results were filtered using Percolator 2.04 (Käll et al., 2007) to yield peptides with a false discovery rate of <1%. The raw data for unique peptides were log transformed and normalized (central mean tendency), and fold changes were calculated using InfernoRDN software (Polpitiya et al., 2008). The log-transformed, normalized data were exported for statistical evaluation using mixed linear modeling (PROC MIXED, SAS software version 9.4). This analytical method has been used previously to analyze MS data (Wildburger et al., 2015) and was selected because of the capacity to handle unbalanced and complex repeated-measures data and the ability to model the variance and correlation structure of repeated-measures experimental designs (Littell et al., 1998). The spectra were nested at the peptide level, and the peptide data were further nested at the protein level. Proteins that were identified by a single peptide were excluded, and a protein was considered differentially expressed if the p value was <0.05. Mascot scores and the percentage of coverage were calculated in Proteome Discover 2.1 software. Although only differentially expressed proteins are reported here, the full raw and normalized datasets are archived in the Monkey Alcohol Tissue Research Resource (www.MATRR.com). Functional annotation of the synaptomics data was performed using Ingenuity Pathway Analysis (QIAGEN) following our previously reported methods (Uys et al., 2016).
Western blotting.
Western blotting procedures for GluA1 followed our routine methods (Padula et al., 2015). Briefly, PSD-enriched fractions from three controls and all eight drinking monkeys were diluted with NuPAGE 4× LDS sample loading buffer (Invitrogen; pH 8.5) containing 50 mm dithiothreitol, and samples were denatured for 10 min at 70°C. Ten micrograms of each sample (controls, n = 3; drinkers, n = 8) was separated using the Bis-Tris (375 mm resolving buffer and 125 mm stacking buffer, pH 6.4; 7.5% acrylamide) discontinuous buffer system with MOPS (3-(N-morpholino)propanesulfonic acid) electrophoresis buffer (50 mm MOPS, 50 mm Tris, 0.1% SDS, 1 mm EDTA, pH 7.7). Protein was then transferred to Immobilon-P PVDF membranes (Millipore) using a semidry transfer apparatus (Bio-Rad). After transfer, blots were washed with PBS containing 0.1% Tween 20 (PBST) and then blocked with PBST containing 5% nonfat dried milk (NFDM) for 1 h at room temperature with agitation. The membranes were then incubated overnight at 4°C with primary antibody against GluA1 (catalog #AB1504, Millipore; RRID: AB_2113602) diluted 1:4000 in PBST containing 0.5% NFDM and washed in PBST before a 1 h incubation at room temperature with horseradish peroxidase-conjugated secondary antibody diluted 1:2000 in PBST. Membranes received a final wash in PBST, and the antigen–antibody complex was detected by enhanced chemiluminescence using a ChemiDoc MP Imaging system (Bio-Rad). Before running the experimental samples, Western blots were performed using different titrations of sample and antibody to establish the linear range for GluA1 in cynomolgus macaque tissue. Bands were quantified by mean optical density using Image Lab software (version 4.0.1; Bio-Rad) and were normalized to a total protein stain (Gürtler et al., 2013). Normalized protein expression data were analyzed with a one-tailed t test, and standard linear regression was used to analyze protein expression levels with optical density values and average daily alcohol intake. The specificity of this antibody for Western blots has been confirmed in GluA1 knock-out mice (Zamanillo et al., 1999) and characterized using human, macaque, and rodent tissue from frontal cortex (Tucholski et al., 2014).
Partial least-squares discriminant analysis regression.
To identify proteins that best discriminate between control and drinking monkeys, we used a partial least-squares discriminant analysis (PLS-DA) regression approach (SOLO, Eigenvector Research). PLS-DA is a multivariate technique that maximizes the separation of the between-groups covariance matrix in large datasets and has been validated as an exploratory approach to identify potential biomarkers in MS data (Jonsson et al., 2005; Rajalahti et al., 2009a). PLS is an extension of multiple regression that is used as a means of exploratory analysis and variable selection when faced with a large number of independent variables and often a relatively small number of dependent variables. For this analysis, monkeys were classified according to their drinking history on the y matrix. The x matrix consisted of log-transformed, normalized (mean centered) spectra data that were rolled up to the protein level using Inferno software. We used the selectivity ratio, calculated as the ratio between explained and residual variance, to determine the proteins with the highest proportion of variance that best discriminate phenotypes, as previously described for identifying biomarkers in MS profiles (Rajalahti et al., 2009a,b). Proteins are reported that have selectivity ratios above the 95% confidence limit that is based on the Hotelling T2 limit, scaled by the eigenvalue of the given principle components.
Drinking-in-the-dark mouse model.
Fifteen female BXD RI strains were obtained from the University of Tennessee Health Science Center and were individually housed in temperature- and humidity-controlled environments. Adult mice (n = 73 total mice; 8–10 weeks old at the start of the experiment) were kept on a 12 h light/dark cycle, and food and water were available ad libitum during the experiment except when alcohol was presented in their home cages. The Pennsylvania State University Institutional Animal Care and Use Committee approved all procedures in accordance with NIH guidelines for the humane care and use of laboratory animals. Binge-like ethanol consumption was induced using a standard 4 d drinking-in-the-dark (DID) protocol (Rhodes et al., 2005; Rinker et al., 2016a). During test week 1, water bottles were replaced with 20% alcohol (v/v) for a 2 h period beginning 3 h into the dark cycle for 3 consecutive days. On the fourth day, alcohol was available for an extended 4 h period. This paradigm was repeated for two additional weeks. Due to an unanticipated increase in the temperature (>30°C) and humidity of the colony room during test week 1 that caused a spike in alcohol drinking, these data were excluded from analysis. Drinking values on day 4 were averaged for each strain and then averaged across test weeks 2 and 3. A two-tailed, unpaired t test was used for analysis in GraphPad software. Strains were genotyped at the C2cd2l interval in GeneNetwork, and all drinking data are available in GeneNetwork.
Bioinformatics analyses.
To provide additional evidence for the alcohol-sensitive proteins and the proteins that best discriminate between low- and heavy-drinking monkeys, we performed two additional open source integrative bioinformatics analyses following our previously reported methods (Padula et al., 2015; McGuier et al., 2016; Rinker et al., 2017). First, we queried the GeneWeaver software system (www.geneweaver.org), a database containing major curated repositories as well as functional genomics results obtained from experiments across nine species (Baker et al., 2012), to retrieve data that implicate the genes that encode the identified proteins in alcohol-related phenomena. Next, we performed targeted analyses using existing genetic and phenotypic data in GeneNetwork (www.genenetwork.org). We correlated PFC robust multi-array average levels in male and female alcohol-naive BXD RI strains of mice [GeneNetwork dataset: VCU BXD PFC Sal M430 2.0 (Dec06) RMA] for C2cd2l (also known as Tmem24; probe set: 1458713_at and 1428095_a_at), Diras2 (1455436_at), and Pycr2 (1448315_a_at) with published phenotypic data, including data for intravenous cocaine self-administration and cocaine-induced conditioned place preference in BXD RI strains (Philip et al., 2010; Dickson et al., 2016). To match the monkey-drinking paradigm, only significant Spearman correlations related to alcohol consumption or preference in two-bottle choice (alcohol vs water) voluntary models were reported. All datasets generated on alcohol and cocaine for this article can be queried on GeneWeaver and GeneNetwork.
Results
Functional studies of OFC pyramidal neurons from heavy drinking macaques
As an initial measure of functional alterations, we recorded passive and active membrane properties of OFC pyramidal neurons in brain slices from control and drinking monkeys. The majority of monkeys used for electrophysiological studies consumed ∼3 g/kg/d (equivalent to ≥12 daily human drinks) and achieved blood alcohol concentrations (BACs) >80 mg/dl that correlated with alcohol intake (Fig. 1A–C). Consistent with known membrane properties of cortical pyramidal neurons in primates (Povysheva et al., 2006), the spike frequency adaptation ratio of deep-layer area 13L pyramidal neurons in controls was 0.45 ± 0.04 (n = 5 cells from two monkeys). All measures of passive membrane (i.e., cell capacitance, resting membrane potential, input resistance) and biophysical properties (i.e., threshold, peak, rise time, half-width, decay, adaptation ratio, AHP amplitude) of the evoked action potentials (APs) were similar between controls and drinkers (Table 1). However, there was a significant reduction in the number of evoked APs (two-way repeated-measures ANOVA, F(14,308) = 2.49, p = 0.002; Sidak's post hoc test, *p < 0.05 vs controls; Fig. 1F,G) and a corresponding increase in rheobase current in OFC neurons from drinking monkeys [controls, n = 5/2 (cells/monkeys); drinkers, n = 19/8; unpaired two-tailed t test, t(22) = 2.42, *p = 0.0239 vs controls; Fig. 1H]. Bath application of 22 (∼100 mg/dl) and 66 mm (∼300 mg/dl) alcohol significantly reduced evoked AP firing in pyramidal neurons in controls (two-way repeated-measures ANOVA, F(28,56) = 1.70, p = 0.046; Sidak's post hoc test, *p < 0.05; n = 3/1; Fig. 1I), but not in drinking monkeys (two-way repeated-measures ANOVA, F(2,16) = 0.131; p = 0.88, n = 9/7; Fig. 1J).
Long-term alcohol drinking reduces evoked spiking in area 13L pyramidal neurons. A, B, Alcohol consumption levels (A; mean ± SD) and BACs (B; mean ± SD) during the 6 months of 22 h access to alcohol self-administration in the INIA Cyno 13 cohort of monkeys used for electrophysiological analysis. C, Representative correlations between alcohol intake and BACs. D, Example image of the ventral aspect of a cynomolgus macaque brain showing prominent anatomical markers and areas 10, 12, 13L, 13M, and 14. E, A representative example of a coronal section of the OFC-containing area 13L that was used for functional and biochemical analyses. F, Representative traces showing reduced evoked spiking in a control and a drinking monkey. G, Number of action potentials at current injections ranging from −40 to −240 pA in control and drinking monkeys. H, Rheobase was increased in the drinking monkeys. I, J, Representative traces and average number of evoked firing in the absence and presence of 22 mm alcohol bath application in control (I) and drinking (J) monkeys. Electrophysiological data are expressed as the mean ± SEM.
Measures of passive membrane and biophysical properties of the evoked APs in deep-layer OFC pyramidal neurons from control and long-term drinking monkeys
Because plasticity of intrinsic excitability coincides with postsynaptic adaptations during learning or alcohol and drug exposure (Sehgal et al., 2013; Kourrich et al., 2015; Nimitvilai et al., 2016), we next determined whether there were functional synaptic adaptations in drinking monkeys. Representative traces of sPSCs obtained from a drinking (gray) and a control (black) monkey are shown in Figure 2, A and B. The average frequency, amplitude, and rise time of sPSCs in OFC pyramidal neurons from control monkeys (Fig. 2C) are consistent with postsynaptic events recorded from pyramidal neurons in monkey prefrontal cortex (Amatrudo et al., 2012; Medalla and Luebke, 2015). Analysis of the sPSCs shows that chronic alcohol self-administration significantly increased amplitude and frequency without altering half-width, rise time, or decay time constant [amplitude: t test, t(20) = 2.057, *p = 0.027; frequency: t test, t(20) = 2.24, *p = 0.037; half-width: t test, t(20) = 1.24, p = 0.11; rise time: t test, t(20) = 1.12, p = 0.45; decay: t test, t(20) = 1.27, p = 0.11; controls, n = 3/2 (cells/monkeys); drinkers, n = 19/5; Fig. 2C]. Frequency histograms for binned (two-way repeated-measures ANOVA, F(13,260) = 5.46, p < 0.0001; Fisher LSD post hoc, *p < 0.007; Fig. 2D) and cumulative (Kolmogorov–Smirnov test, *p < 0.0001; Fig. 2E) distribution revealed rightward shifts of sPSCs amplitudes in pyramidal neurons recorded from drinkers. Together, these functional analyses of OFC pyramidal neurons indicate that alcohol self-administration produces opposing effects on evoked action potential firing and postsynaptic signaling.
Long-term alcohol self-administration enhances the amplitude and frequency of sPSCs in deep-layer area 13L pyramidal neurons. A, B, Representative traces of sPSCs from a control (black trace) and a drinking (gray trace) monkey. C, Average amplitude, frequency, half-width, rise time, and decay in controls and drinkers. Data are expressed as the mean ± SEM. D, E, Frequency histograms for binned and cumulative sPSC amplitude distributions in controls and drinking monkeys.
Synaptomic analysis in drinking monkeys
An unbiased quantitative synaptomics approach was then used to identify synaptic proteins that were altered by drinking. A separate cohort of cynomolgus macaques was used for this study, and their drinking data are shown in Figure 3A–C. After 6 months of free access to alcohol, brains were extracted before the start of a drinking session and PSD-enriched fractions were prepared from area 13L for proteomic analysis. LC/MS-MS analysis identified 976 proteins with isoforms (879 unique proteins and 3809 unique peptides). A normality test of the normalized data fit the expression ratio frequency distribution (R2 = 0.988; Fig. 3D,E). Functional annotation revealed that 33.4% of the unique proteins localize to the plasma membrane or were identified as ion channels or membrane-bound transporters. Analysis showed that long-term drinking significantly altered the expression of a diverse group of 57 proteins (Table 2), of which 29 are transmembrane proteins and 47 were significantly upregulated. We then analyzed all differentially expressed OFC proteins to identify prominent biological networks that are affected by long-term drinking. Proteins linked to cell-to-cell signaling and interaction, cell signaling, and nucleic acid metabolism were the top three biological function annotations of the alcohol-sensitive proteins (Fig. 3F). Data were filtered using a conservative threshold of p < 0.01 (−log(B-H corrected p value) >2). Because Triton X-100 fractions are enriched in excitatory synaptic proteins (Uys et al., 2016), it is possible that these three annotations are also highly ranked when considering the entire dataset. However, analysis showed that these biological annotations were ranked 1st, 6th, and 23rd, respectively, demonstrating that alcohol self-administration affects subsets of OFC proteins. A closer inspection of the significantly dysregulated proteins reveals apparent coordinated adaptations in 23 proteins that control presynaptic glutamate release and postsynaptic glutamatergic signaling (Fig. 3G, Table 2). In addition to excitatory synaptic proteins, GABAA receptors (Centanni et al., 2014) and other key scaffolding and signaling proteins expressed in GABAergic inhibitory synapses (Sassoè-Pognetto et al., 2011) are also present in the Triton X-100 insoluble fraction. However, in contrast with adaptations in excitatory proteins, GABAergic inhibitory synaptic proteins present in the dataset were unaltered by long-term alcohol consumption (Table 3).
Synaptomic profiling of PSD-enriched fractions from control and drinking cynomolgus macaques identifies novel alcohol-sensitive proteins (controls, n = 3; drinkers, n = 5). A, B, Average daily alcohol intake (A) and BACs (B) during the 6 months of self-administration. The monkeys used for MS/MS analysis are shown in the hatched bars. Two of the monkeys (10201 and 10207) reached criteria for heavy drinking (≥3 g/kg on 48% and 38% of the open-access drinking days), and the remaining monkeys fit the criteria for low drinkers achieving ≤3 g/kg intake on <1% of the open access days. Data are the mean ± SD. C, Correlation analysis of average BAC and alcohol intake (g/kg) at 7 h into the daily session for all eight drinking monkeys. D, E, Frequency distribution of the normalized log-transformed median expression ratios of the 976 proteins. The histogram was fit by a normal distribution. Frequency distribution of the ratios is shown in blue bars, and a red line shows the results of the normal fit. F, IPA Core analysis identified highly ranked biological functions of the Triton X-100 insoluble proteins in the area 13L that were differentially expressed by alcohol self-administration. The red line denotes threshold, and the number above each bar shows the number of proteins contained within that annotation. G, Presynaptic and postsynaptic proteins related to glutamatergic signaling that were significantly different between control and drinking monkeys. H, GeneWeaver search results of alcohol-related gene sets (GS) containing GPM6A and PRRT2.
Long-term drinking induces expression changes in proteins in the macaque OFC
Long-term alcohol self-administration does not affect expression levels of key proteins found in GABAergic inhibitory synapses
The increase in GluA1 expression found in the synaptomics screen (Table 2) and demonstration of enhanced postsynaptic currents from the electrophysiology experiments suggests an increase in postsynaptic AMPA receptor function (Fig. 2) in drinking monkeys. Because AMPA receptor expression is increased in the lateral OFC of alcohol-dependent mice (Nimitvilai et al., 2016), GluA1 expression levels in this cohort of monkeys were measured using Western blotting to validate cross-species effect of alcohol on AMPA receptors. Characterization of GluA1 immunoreactivity in monkey OFC samples revealed a linear dynamic range across twofold dilutions between 1.25 and 40 μg of protein (R2 = 0.9956; Fig. 4A,B). We also detected significant differences in total protein intensity and the raw GluA1 signal when there was an intentional 17% difference in control sample load (total protein stain: t test, t(14) = 2.85, *p = 0.013, n = 8/group; anti-GluA1: t test, t(14) = 2.32, *p = 0.036; n = 8/group; Fig. 4C), further demonstrating that this quantitative approach is reliable and sensitive (Gürtler et al., 2013). Consistent with results from the proteomics and functional studies, GluA1 level was significantly elevated in drinking monkeys compared with nondrinking controls (one-tailed t test, t(9) = 2.00, *p = 0.039; controls, n = 3; drinkers, n = 8; Fig. 4D).
Western blot characterization and analysis of synaptic GluA1 expression in the OFC of drinking monkeys. A, Characterization of GluA1 Western blot in macaque cortical tissue (protein loading range, 1.25–40 μg). B, Positive correlation between the amount of protein loaded and GluA1 optical density values. C, Representative images and quantitation of an intentional 17% decrease in total protein load. D, Representative GluA1 blot and quantitation of normalized GluA1 expression in controls and drinkers. Data are the mean ± SEM.
Although these data validate cross-species effects of alcohol on GluA1 expression, a primary goal of synaptic profiling was to identify novel proteins that underpin drinking-induced neuroadaptations. To our knowledge, this is the first study to link alcohol with two synaptic plasma membrane proteins: proline-rich transmembrane protein 2 (PRRT2) and neuronal membrane glycoprotein M6-a (GPM6A). Since these membrane proteins showed the largest fold change, we further explored their relationship with drinking by querying the GeneWeaver.org database (Baker et al., 2012) for functional genomic experiments related to alcohol intake and other alcohol-related behaviors across species. GPM6A is found in 11 published quantitative trait loci (QTLs) for alcohol drinking and other alcohol-related behaviors in humans, rats, and mice (Fig. 3H), while PRRT2 is found in four alcohol-related gene sets (Fig. 3H). PRRT2 is an outer core AMPA receptor complex protein (Schwenk et al., 2012), and GRIA1, the gene encoding the GluA1 subunit of AMPA receptors, is found in 19 alcohol-related gene sets that cross multiple species (GeneWeaver.org analysis not shown). Thus, the synaptomic and functional genomic analyses identified novel alcohol-sensitive synaptic proteins with cross-species genetic links to alcohol consumption phenotypes.
Identification and validation of heavy-drinking synaptic biomarkers
While a recent study has identified a biomarker in blood and a subcortical brain region that predicted heavy alcohol consumption in monkeys (Cervera-Juanes et al., 2016), the identity of synaptic proteins that may drive excessive drinking remains largely unknown. We thus conducted a final exploratory analysis of the synaptomic data using partial least squares discriminant analysis (PLS-DA) to determine which key proteins or groups of proteins best discriminate between control and drinking monkeys. Despite marked changes in functional measures of synaptic excitability in drinking monkeys, PLS-DA only identified two proteins [myosin light chain 6B (MYL6B) and neuronal cell adhesion molecule (NRCAM)] that were above the selectivity ratio criteria that discriminates between controls and drinkers. The low number of discriminatory proteins identified in this analysis suggests unexplained variance in the proteomic dataset, possibly stemming from the allelic diversity reported in alcohol-naive cynomolgus macaques (Kanthaswamy et al., 2013). An equally plausible hypothesis is that individual differences in alcohol intake may contribute to the unexplained variance. To account for this variance, monkeys were classified into low, binge, or heavy drinkers using previously described criteria (Baker et al., 2014). Two of the monkeys (Fig. 3) reached the criteria for heavy drinking (≥3 g/kg on 48% and 38% of the open-access drinking days), and the remaining monkeys fit the criteria for low drinkers achieving ≤3 g/kg intake on <1% of the open access days. PLS-DA of these two groups identified 13 proteins that best discriminated between low and high drinking phenotypes (Fig. 5A). Subsets of the family of 13 proteins fall into functional classes, as 4 proteins have biological processes involving actin cycling (ADD1, ADD3, WASF1, and SPTBN4) and 3 proteins control clathrin-mediated endocytosis (AP2S1, AP2A1, and EPS15L1). Two of the five top proteins are addiction-related proteins in the adducin family (Jung et al., 2013; Han et al., 2015), thus supporting the use of this approach to identify novel brain biomarkers.
Exploratory analysis of the orbitofrontal cortex synaptome identified proteins that best discriminate between low- and high-drinking cynomolgus monkeys and genes with cross-species links to alcohol- and cocaine-related phenotypes. A, Selectivity ratio plot for all proteins identified in the synaptomics screen. Proteins with selectivity ratios above the 95% confidence limit are identified by gene name. B, GeneWeaver search results of alcohol-related gene sets (GS) containing C2CD2L, DIRAS2, and PYCR2. C, D, Correlations between C2cd2l transcript levels in the PFC and cocaine self-administration and cocaine-conditioned place preference in BXD RI strains.
Interestingly, genes encoding the remaining three top proteins (C2CD2L, DIRAS2, PYCR2) are located within published human and rodent QTLs for alcohol drinking and other alcohol-related behaviors (Fig. 5B). To further validate these initial genetic findings, we experimentally tested associations between PFC transcript levels for these genes and voluntary alcohol drinking in BXD RI strains of mice using the GeneNetwork.org software system (Baker et al., 2012). BXD RI strains generated by ethanol-preferring C57BL/6J (B) and ethanol-avoiding DBA/2J (D) mice are valuable for studying genetic drivers of variation in alcohol consumption (Philip et al., 2010) and our previous work used this approach to identify candidate genes as pharmacological targets for reducing heavy drinking (Padula et al., 2015; Rinker et al., 2017). In alcohol-naive BXD strains, PFC transcript levels for C2cd2l, Diras2, and Pycr2 significantly correlate with alcohol consumption in a standard two-bottle choice (alcohol vs water) drinking model (Table 4). C2cd2l was negatively correlated with 12 datasets for alcohol consumption and preference in male and female BXD RI strains across a range of alcohol concentrations (3–15%, v/v). Because of the relationship between C2CD2L and drinking in monkeys and mice, we then exploited the DBA/2J sequence (Keane et al., 2011; Wang et al., 2016) to extract all sequence variants in and around C2cd2l (e.g., rs32609282). While there are no high-impact C2cd2l variants (Wang et al., 2016), C2cd2l is associated with cis-acting expression QTLs (logarithm of odds scores >3; Pandey and Williams, 2014) in multiple addiction-related brain regions. As neurocircuitry and mechanisms of drug-seeking are shared across drug classes (Mulholland et al., 2016), PFC transcript levels for C2cd2l were then correlated with behaviors related to cocaine experience (Philip et al., 2010; Dickson et al., 2016). Comparable to the relationship between C2cd2l and alcohol intake, there were negative correlations between C2cd2l and cocaine self-administration, as well as cocaine-conditioned place preference (Fig. 5C,D). C2cd2l is also found in QTLs for alcohol-, cocaine-, methamphetamine-, morphine-, and nicotine-related phenotypes (GeneWeaver.org analysis not shown), suggesting that C2CD2L may be associated with genetic vulnerability associated with multiple abused substances.
Correlations between PFC transcript levels in male and female alcohol-naive BXD recombinant inbred strains of mice and alcohol consumption and preference in two-bottle choice (alcohol vs water) drinking models
As a final validation for the role of C2cd2l in alcohol drinking, we used a reverse genetic approach (Wang et al., 2016) that can evaluate the spectrum of phenotypes directly linked to the C2cd2l interval and compared alcohol intake in BXD strains that inherited the D allele in comparison with the B allele. Fifteen strains of female BXD mice were allowed to consume alcohol in a binge-like DID model that promotes high levels of alcohol intake. As expected, average alcohol intake in BXDs during the 4 h test drinking sessions varied across strains (range, 2.4–9.5 g/kg/4 h; n = 73 mice; Fig. 6A), and strains that inherited the D allele at C2cd2l on average drank markedly less than strains that inherited the B allele (GeneNetwork record IDs: 13568 and 13571; t test: t(13) = 4.75, ***p < 0.001; Fig. 6B). Consistent with these results, examination of the C2cd2l interval using additional existing datasets in GeneNetwork.org shows that male and female strains that inherited the D allele at C2cd2l consumed significantly lower amounts of alcohol in standard two-bottle choice drinking models [24 h model, 10% alcohol: GeneNetwork record IDs: 10476 (t test, t(19) = 3.58, **p < 0.002), 10477 (t test, t(19) = 2.49, *p < 0.05), and 10582 (t test, t(21) = 2.15, *p = 0.043); 2 h model, 15% alcohol: GeneNetwork record ID 12745 (t test, t(31) = 2.33, *p < 0.027); Fig. 6C,D). Together, these converging analyses corroborate C2CD2L as a high-priority candidate for differential alcohol consumption across species.
Alcohol consumption in BXD RI strains of mice that inherited the B or D allele at the C2cd2l interval. A, Alcohol (20%, v/v) intake across 15 BXD RI strains in a single-bottle, 4 h DID model. The number within each bar for A shows the number of BXD RI mice within each strain. B–D, BXD mice that inherited the D allele at C2cd2l consumed significantly less alcohol in the DID model (B) and two-bottle choice (C) 24 and 2 h access (D) models. The number within each bar shows the number of strains with the B or D alleles for the C2cd2l interval for B–D.
Discussion
A major finding of this study is that long-term drinking produces functional adaptations in OFC pyramidal neurons of macaques. The electrophysiological analysis revealed functional alterations, while the synaptome screen identified novel proteins altered by alcohol self-administration including presynaptic and postsynaptic proteins that regulate glutamatergic, but not GABAergic, signaling. Because tissue samples for all of these studies were acquired before the start of a daily drinking session, these neuroadaptations may reflect an abnormal state of OFC pyramidal neurons that is present in nonabstinent, actively drinking monkeys. Exploratory integrative bioinformatic approaches unveiled a group of proteins that best discriminate between low and heavy alcohol-drinking subjects. The characterization of these proteins and their contribution to the neurophysiological mechanisms of heavy alcohol drinking is an important step in the development of potentially clinically effective therapeutic strategies to target individuals with AUD.
Long-term drinking and functional adaptations
In this study, we show opposing actions of long-term drinking on evoked action potential firing and spontaneous synaptic activity of monkey OFC pyramidal neuron physiology. Drinking reduced evoked firing and markedly enhanced the amplitude and frequency of sPSCs recorded from OFC pyramidal neurons. These findings are in opposition to intrinsic and synaptic plasticity in models of learning that report learning-related enhanced dendritic and somatic excitability as an early mechanism that facilitates synaptic plasticity (Sehgal et al., 2013). However, consistent with our findings, opposing changes in intrinsic and synaptic plasticity have been observed in the nucleus accumbens shell of rats that self-administered cocaine and may reflect a homeostatic response to normalize global neuronal activity or a shift in the signal-to-noise ratio of behaviorally relevant inputs (Kourrich et al., 2015). Regardless of which hypothesis is correct, these findings suggest that abused substances produce aberrant plasticity of somatic and dendritic excitability of neurons in the corticostriatal circuitry. While enhanced postsynaptic excitability and loss of acute alcohol inhibition of firing are consistent with our previously published study in mice (Nimitvilai et al., 2016), OFC neurons in alcohol-dependent mice show enhanced firing following withdrawal. This apparent discrepancy may reflect the use of different withdrawal times between the mouse (3–14 d) and monkey (<6 h) studies. To directly examine this question, we measured evoked firing of OFC neurons prepared from alcohol-dependent mice within 2 h of their last ethanol exposure, similar to the time of the monkey OFC recordings. Compared with air-exposed control mice, there was no significant effect of long-term intermittent alcohol exposure on OFC neuron firing (two-way repeated-measures ANOVA, F(1,29) = 0.624, p = 0.449). These results are consistent with those results from other mouse studies showing withdrawal time-dependent differences in markers of cortical neuron excitability (Holmes et al., 2012; Kroener et al., 2012; McGuier et al., 2015; Pleil et al., 2015; Nimitvilai et al., 2016). With regard to the effects of firing rate changes on behavior, in vivo recordings in mice (Bissonette et al., 2015) and monkeys (Chang et al., 2005) show that both low and high firing rates are related to poor performance on OFC- and PFC-dependent tasks, suggesting that disruption of optimal firing of OFC pyramidal neurons by long-term alcohol consumption may underlie cognitive decline in alcoholic individuals.
Long-term alcohol consumption and AMPA receptors
A notable result from profiling the OFC synaptome was the finding that long-term drinking enhances GluA1 expression, a result confirmed by Western blotting and consistent with the increase in the amplitude of putative AMPA-mediated sPSCs in deep-layer OFC pyramidal neurons observed in drinking monkeys. While we cannot rule out that increases in sPSC amplitude reflect changes in GABA-mediated signaling, GRIA1 levels are also increased in the superior frontal cortex of alcoholic individuals (Lewohl et al., 2000), and increases in AMPA receptor function and expression have been reported in the lateral OFC, but not prelimbic PFC, of alcohol-dependent mice (Kroener et al., 2012; Hu et al., 2015; Nimitvilai et al., 2016). Our recent chemogenetic and lesion study has implicated the OFC in controlling escalated drinking in dependent mice (den Hartog et al., 2016), and there is evidence from studies with rats showing that AMPA receptors affect operant self-administration of alcohol (Wang et al., 2012; Cannady et al., 2016). Interestingly, sPSC frequency was also elevated in the drinking monkeys, and the proteomics screen identified multiple upregulated presynaptic proteins that control glutamate release from synaptic vesicles. A recent study has reported that SV2A, an alcohol-responsive presynaptic protein, regulates dendritic spine density and postsynaptic AMPA receptor expression (Cohen et al., 2011). Although we did not examine dendritic spine changes in drinking monkeys, alcohol dependence enhanced dendritic spine density in the mouse OFC (McGuier et al., 2015). Together, this converging evidence reveals a cross-species effect of long-term alcohol consumption on the GluA1 subunit of AMPA receptors in the OFC and suggests that presynaptic adaptations produced by drinking may be molecular mechanisms that drive functional and morphological postsynaptic adaptations.
Novel alcohol-sensitive synaptic proteins
An advantage of our approach over previous proteomic studies in brain tissue is the ability to solubilize and quantify transmembrane synaptic proteins. Our profiling of the macaque OFC synaptome identified novel alcohol-sensitive transmembrane proteins, and the two synaptic transmembrane proteins that were most altered by long-term drinking (i.e., GPM6A and PRRT2) localize to presynaptic sites where they control glutamate release (Roussel et al., 1998; Valente et al., 2016). In addition, postsynaptic GPM6A has been reported to promote dendritic spine formation (Alfonso et al., 2005), whereas PRRT2 functions in spines to modulate AMPA receptor gating properties and surface expression (Schwenk et al., 2012; Li et al., 2013). Recent studies suggest a central role for GPM6A and PRRT2 in cognitive impairments (Najmabadi et al., 2011; Gregor et al., 2014), and our bioinformatic analysis links the genes encoding these proteins to alcohol drinking and other alcohol-related behaviors. Loss-of-function PRRT2 mutations are a causative factor for familial paroxysmal kinesigenic dyskinesia (PKD; Heron and Dibbens, 2013) and PKD patients with an insertion mutation in PRRT2 show enhanced spontaneous low-frequency fluctuations in corticostriatal circuitry (Luo et al., 2013). Neuroimaging studies in alcoholic individuals reveal aberrant executive control network activity and altered brain connectivity (Müller-Oehring et al., 2015; Zhu et al., 2015) that could reflect disruptions in GPM6A/PRRT2 expression or function. Future studies are needed to validate the importance of these proteins in mediating the altered cognition, glutamatergic signaling, and corticostriatal circuitry observed in individuals with AUD.
Our final query of the OFC synaptome was an exploratory analysis that identified functional groups of proteins that are discriminators of a heavy-drinking phenotype. The subset of identified proteins that best discriminated between high and low drinkers is the known regulators of actin cycling within dendritic spines and has been previously implicated in drug addiction (Jung et al., 2013; Han et al., 2015). Adaptations in dendritic spines on neurons within the corticostriatal circuitry are a shared feature of drug-induced plasticity (Mulholland et al., 2016). Additional bioinformatic analysis revealed that C2CD2L, DIRAS2, and PYCR2 are found within multiple published QTLs for alcohol drinking and dependence in human and rodent studies, suggesting that they may be relevant brain biomarkers for a heavy-drinking phenotype. In support of this, PFC transcript levels for these three genes in BXD strains correlate with voluntary drinking. Like Diras2, several other genes encoding proteins in the RAS family are associated with excessive drinking (Repunte-Canonigo et al., 2010; Ben Hamida et al., 2012; Stacey et al., 2012; Ojelade et al., 2015). Finally, strains with the C2cd2l D allele drink markedly less alcohol across multiple paradigms, including a binge-like drinking model. Thus, these findings suggest that DIRAS2 and C2CD2L are brain biomarkers of heavy drinking and should be explored and validated as potential targets in preclinical models of excessive drinking.
Although the findings reported here uncovered adaptations and potential pharmacotherapeutic targets for the treatment of AUD, there are some limitations to consider. First, sample sizes for the functional measures were unequal due to a limited number of monkeys assigned to the control condition. During these studies, the experimenter was blind to the treatment groups so no adjustments could be made to account for this difference. While the limited number of control monkeys is potentially a concern, the robust effects of drinking on sPSCs and the congruence of the proteomic changes in presynaptic and postsynaptic glutamatergic signaling machinery help to mitigate this limitation. Second, the overrepresentation of high drinkers used for these studies precludes correlations between functional adaptations and individual differences in alcohol intake. Nonetheless, the key proteins that discriminated between low and high male drinkers were validated in multiple alcohol-drinking models using both male and female mice. In addition, as previously discussed, the recording conditions used to monitor sPSCs, although biased in favor of AMPA events, do not rule out the possibility of changes in inhibitory synaptic transmission. Last, because data were collected at a single time point from long-term drinking monkeys, it is unknown at what stage the observed functional neuroadaptations occurred and whether they would persist. These are important considerations as such changes may be long-lasting and could contribute to relapse to alcohol-seeking behaviors even after long periods of abstinence. Future studies are needed to address these questions.
In summary, long-term alcohol self-administration is associated with complex neuroadaptations in deep-layer OFC pyramidal neurons in heavy-drinking macaques. The congruence of alcohol effects across species suggests that aberrant adaptations in glutamatergic signaling are key mechanisms underpinning alcohol-induced OFC dysfunction. Finally, the novel proteins (e.g., C2CD2L) that show links to a heavy-drinking phenotype in rodents, monkeys, and humans have potential as targets for the treatment of AUD. Further clinical and preclinical studies are necessary to determine whether targeting these biomarkers prevents excessive intake and reduces the incidence of relapse.
Footnotes
This work was supported by National Institutes of Health (NIH) Grants AA-020930, AA-023288, RR-024485, AA-024426, AA-019431, AA-013541, AA-016662, AA-013499, AA-010761, AA-009986, AA-021951, and OD-11092. The Medical University of South Carolina Mass Spectrometry Facility receives support from the SC COBRE in Oxidants, Redox Balance, and Stress Signaling (Grant P20 GM103542) and the Office of the Provost. The Orbitrap Elite Mass Spectrometer was acquired through NIH/National Center for Research Resources Grant S10-OD-010731. We thank Susana Comte-Walters. We also thank Dr. James B. Daunais for the macaque frontal cortex images.
The authors declare no competing financial interests.
- Correspondence should be addressed to Dr. Patrick J. Mulholland, Associate Professor, Medical University of South Carolina, Charleston Alcohol Research Center, Department of Neuroscience and Psychiatry and Behavioral Sciences, 67 President Street, MSC 861/IOP 462N, Charleston, SC 29425-8610. mulholl{at}musc.edu.