MAOA EXPRESSION PREDICTS VULNERABILITY FOR ALCOHOL USE
Abstract
The role of the monoamines dopamine (DA) and serotonin (5HT) and the monoamine-metabolizing enzyme monoamine oxidase A (MAOA) have been repeatedly implicated in studies of alcohol use and dependence. Genetic investigations of MAOA have yielded conflicting associations between a common polymorphism (MAOA-LPR) and risk for alcohol abuse. The present study provides direct comparison of tissue-specific MAOA expression and the level of alcohol consumption. We analyzed rhesus macaque MAOA (rhMAOA) expression in blood from males before and after 12-months of alcohol self-administration. In addition, nucleus accumbens core (NAc core) and cerebrospinal fluid (CSF) were collected from alcohol-access and control (no alcohol access) subjects at the 12-month time point for comparison. The rhMAOA expression level in the blood of alcohol-naïve subjects was negatively correlated with subsequent alcohol consumption level. The mRNA expression was independent of rhMAOA-LPR genotype and global promoter methylation. After 12 months of alcohol use, blood rhMAOA expression had decreased in an alcohol dose-dependent manner. Also after 12 months, rhMAOA expression in the NAc core was significantly lower in the heavy drinkers, as compared to control subjects. The CSF measured higher levels of DA and lower DOPAC/DA ratios amongst the heavy drinkers at the same time point. These results provide novel evidence that blood MAOA expression predicts alcohol consumption and that heavy alcohol use is linked to low MAOA expression in both the blood and NAc core. Together, the findings suggest a mechanistic link between dampened MAOA expression, elevated DA and alcohol abuse.
INTRODUCTION
Sustained alcohol use can lead to abuse and dependence, collectively known as alcohol use disorders (AUDs), and ultimately to life-threatening medical, psychological and social consequences1. Risk for hazardous drinking has been ascribed to an array of genetic2, 3 and environmental factors4, 5, while the transition to dependence involves neuronal adaptations to alcohol exposure that evolve over time6, 7. These adaptations translate to behavioral changes that persist long after alcohol abstinence, contributing to the high rates of relapse. Ethanol-induced release of serotonin (5HT) and dopamine (DA) in the nucleus accumbens (NAc)8, 9 and ventral tegmental area of the mesoaccumbens reward pathway have been implicated in the mechanisms underlying development of alcohol dependence10–12 and relapse13.
The monoamine oxidase A (MAOA) enzyme metabolizes the monoamines 5HT, DA and norepinephrine, and therefore its activity has direct implications on the function of the reward system. A promoter-linked length polymorphism upstream of the MAOA gene, MAOA-LPR14, has been widely associated with numerous neuropsychiatric disorders, including alcoholism15–18. It has been suggested that the MAOA-LPR promoter alleles may modulate gene expression14, and correspondingly, the downstream levels of these key neurotransmitters. In addition, Caspi and colleagues19 highlighted the critical role of the environment when they reported that an association between the low activity MAOA-LPR allele with impulsive behaviors was contingent upon exposure to developmental stressors. However, years of research have yielded conflicting results as to the specific MAOA-LPR alleles that contribute to risk for psychiatric disorders20–22 and have also reinforced the possibility of an interaction between genotype and adverse life history20–23. Furthermore, recent studies directly exploring MAOA abundance have detected a lack of association between the MAOA-LPR genotype and postmortem24 or in vivo brain MAOA protein levels25, 26. These findings suggest that MAOA-LPR genotype has little effect in modulating brain MAOA levels or that its effect is dependent on other factors.
Epigenetic regulatory mechanisms enable the adaptive and precise modulation of gene expression. Among them, DNA methylation involves the addition of a methyl group to the cytosine residue within CpG dinucleotides. Hypermethylation in regulatory regions can attenuate gene expression, while DNA methylation within the gene-body has been linked to tissue- and cell-specific gene regulation27. Several studies have identified DNA methylation as a potential molecular link between adverse life history events and risk for neuropsychiatric disorders, including addiction28–33. A study of MAOA promoter methylation in women identified an association between methylation level and alcohol dependence16. In addition, recent studies in healthy males demonstrated that variability in MAOA methylation detected in the blood was correlated with MAOA protein activity in the brain indicating that blood-based methylation levels may serve as a useful biomarker for neural MAOA measures25, 26.
The discovery of human genetic and epigenetic liabilities influencing risk for alcohol dependence relies on the use of retrospective and self-reported data, which may be biased or incorrect. Furthermore, significant findings obtained from retrospective studies of alcoholic and healthy populations are limited in their ability to distinguish preexisting conditions from alcohol-induced effects. These limitations can be addressed by using animal models for the study of long-term alcohol self-administration, providing highly accurate alcohol consumption data and the ability to collect blood before and after alcohol use. A well-established model of oral alcohol self-administration in macaques demonstrates a wide range of voluntary alcohol consumption34–36 and recapitulates a variety of alcohol seeking patterns observed among humans. In addition, macaques have the advantage of sharing similar genetic and epigenetic composition, anatomy, physiology, and endocrine characteristics with humans37, 38. Some functionally similar gene polymorphisms, including MAOA-LPR, are also shared in both species. As with humans, specific rhesus MAOA-LPR alleles have been associated with high levels of alcohol consumption37. Thus the NHP alcohol self-administration model provides a relevant and unique opportunity to investigate the genetic, epigenetic and environmental factors contributing to the vulnerability, progression and pathogenesis of AUDs39.
In this study we test the hypothesis that rhesus MAOA (rhMAOA) expression level in the blood of alcohol-naïve subjects is a predictive biomarker of risk for alcohol abuse. We investigated potential genetic and epigenetic mechanisms that might contribute to the initial population variance in gene expression. We then determined the within-individual changes in blood rhMAOA expression and DNA methylation following chronic alcohol use. To explore the association between neural rhMAOA expression levels and alcohol consumption we compared the mRNA levels in the NAc core of control subjects (no-alcohol access) and animals that consumed alcohol for 12 months. Finally, we considered how the NAc core rhMAOA expression relates to monoamine levels by measuring DA, 5HT and their metabolites in the cerebrospinal fluid (CSF) after 12 months of alcohol self-administration, as the CSF temporarily stores these monoamines and their metabolites40.
MATERIALS AND METHODS
Subjects
Male rhesus macaques (n=38, Macaca mulatta) differing slightly in age, late adolescents and young adults (3.0–5.0 and 5.1–7.0 years at the start of the 12 month open access period, respectively; n=19 per age group), were included in this study. All of the monkeys were born and reared at the Oregon National Primate Research Center (ONPRC) with their mothers until 2–3 years of age. All subjects were initially selected to minimize relatedness. The average kinship coefficient of all subjects in the study was 0.003. Ten of the animals were assigned as controls. The control subjects were selected by matching on age and weight, and then kinship. The mean kinship coefficient did not exceed 0.003 in either the control or ethanol access groups. Controls were housed in the same housing rooms as the experimental subjects and underwent the same training for awake blood draws, medical check-ups and MRI imaging. Controls were also placed on the same diet, timing and order of experimental phases and had equal experience with the research technicians. However, control monkeys did not receive access to alcohol. Instead, at baseline, controls were yoked to a future ethanol monkey based on weight, and received a quantity of 10% maltose-dextrin solution matched in calories to the previous day’s intake of their yoked ethanol monkey. Monkeys were individually housed and ethanol self-administration was induced using schedule-induced polydipsia as previously described34, 35. During the following 12 months, subjects had open access to 4% alcohol and water (ethanol subjects) or water only (control subjects). The average g/kg body weight per day (g/kg/day) was based upon all ethanol or water consumed during the 12 months of open access. For the ethanol self-administration data, technicians were not blind to the condition (ethanol or control) of the monkeys, however the intake data were collected and recorded in an automated fashion by computer and analyzed by individuals who did not interact with or know the drinking status of the monkeys. All of the animal procedures used in this study were approved by the ONPRC IACUC and were performed in accordance with the NIH and the National Resource Council’s Guide for the Care and Use of Laboratory Animals.
Genomic DNA isolation
Blood was collected from awake monkeys34 prior to (0 months) and after 12 months of open alcohol access. Genomic DNA (gDNA) was extracted using the 5 PRIME ArchivePure DNA purification kit (Thermo Scientific, Pittsburgh, PA) following the manufacturer’s protocol.
After 12 months of open access to alcohol or water, a detailed necropsy protocol was used to systematically collect tissues from all subjects41. Animals were sedated (ketamine 15mg/kg), CSF was collected, and then the animals were brought into a surgical plane of anesthesia with intravenous administration of sodium pentobarbital (30–50 mg/kg, i.v.). Following extraction, the entire brain was blocked or sectioned into slides that were fresh frozen and stored at −80° C. The NAc core was later dissected from frozen brain slices. The typical postmortem interval for brain extractions was < 5 minutes. Genomic DNA and RNA were extracted using the All Prep DNA/RNA/miRNA Universal kit (Qiagen Sciences Inc, Germantown, MD) following the manufacturer’s instructions.
Genotyping
The rhMAOA-LPR fragment size analysis was performed as previously described42. The rhMAOA-LPR primers (Supplementary Table 1) generated 206 bp, 224 bp and 242 bp products for the 5, 6 and 7 copy number alleles, respectively (Life Technologies, Carlsbad, CA). Since the rhMAOA gene is located on the X chromosome, there was only a single copy of rhMAOA detected in each of the male subjects.
Bisulfite sequencing
Correction for bisulfite-converted PCR bias was carried out as described by Moskalev et al.43. Methylated and unmethylated human gDNAs (Zymo Research, Irvine, CA) were bisulfite-converted using the EZ DNA Methylation-Gold kit (Zymo Research, Irvine, CA), according to the manufacturer’s instructions. The bisulfite converted DNAs were combined to create a series of methylation rate standards ranging from 100% to 0% methylated. gDNA (500ng) from NHP blood and NAc core was bisulfite converted following the same protocol. Primers were designed to amplify a 223 bp region of the MAOA promoter in both humans and rhesus macaques (Supplementary Table 1). Amplification was carried out using a C1000 Thermal Cycler (Bio-Rad, Hercules, CA) with 20 ng of bisulfite-treated DNA per PCR reaction. The amplification conditions were as follows: Phase 1: 10 cycles of 94°C for 30sec., 68°C for 1min., with a touchdown decrease of 1°C per cycle. Phase 2: 28 cycles of 94°C for 30sec. and 58°C for 45sec. Sequencing libraries were prepared using the NETflex DNA Sequencing Kit (BIOO Scientific, Austin, TX) according to manufacturer’s instructions. The libraries were evaluated using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA) and were normalized to 2nM with 10 mM Tris-HCl. The libraries were then pooled and sequenced on a MiSeq (Illumina, Inc. San Diego, CA) by the Molecular & Cellular Biology Core, (ONPRC, Beaverton, OR) to generate 250base, paired-end reads. The reads were trimmed using Trim Galore44 and aligned to the reference genome using Bismark45. M-bias plots were then generated46, and reads were trimmed further as needed. Bismark alignment data was converted to CpG methylation rate using the Bismark methylation extractor45 and custom scripts.
To assess the linearity of the PCR amplification we plotted the observed versus expected methylation rates obtained from the methylated DNA standard dilution series43, producing a good linear fit (R2=0.98).
RNA isolation and reverse transcription
RNA was extracted from blood collected in PAXgene Blood RNA tubes (PreAnalytix GmbH, Qiagen, Germantown, MD) and stored at −80°C. Total RNA was extracted using the PAXgene Blood mRNA Kit (PreAnalytix GmbH, Qiagen Sciences Inc, Germantown, MD) following the manufacturer’s protocol.
The RNA from the NAc core was extracted using the All Prep DNA/RNA/miRNA Universal kit (Qiagen Sciences, Inc., Germantown, MD) following the manufacturer’s instructions.
RNA quantity and quality was evaluated on a 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA). SuperScript III First-strand Synthesis System (Life Technologies, Carlsbad, CA) was used to reverse-transcribe 500 ng of each RNA sample following the manufacturer’s instructions.
Relative quantitative real-time PCR (qRT-PCR)
The β-Actin gene was selected as the endogenous control since it was equally expressed in all sample types. Standard curves for each gene (rhβ-Actin and rhMAOA), were prepared using serial dilutions of cDNA (from 25 ng to 0.125 ng). Real-time PCR was performed using a QuantStudio 12K Flex (Life Technologies, Carlsbad, CA) with the TaqMan Gene Expression Kit. Each reaction (10μl) contained 5μl of the 2x TaqMan Gene Expression Master Mix (Life Technologies, Carlsbad, CA) and 1μl (12.5ng) of cDNA. For rhβ-Actin, we used 300nM amplification primers with 250nM of rhβ-Actin probe. For rhMAOA, 0.5μl of the custom 20x TaqMan Gene Expression reagents (Supplementary Table 2) was added. The reactions were incubated at 50°C for 2 min., and 95°C for 10 min., followed by 40 cycles of 15 sec at 95°C and 1 min. at 60°C. The QuantStudio 12K Flex system software was used to calculate Ct values. Data from the dilution series were used to generate standard curves for both the rhβ-Actin control and target. Target values were normalized using the endogenous control. Sample 10086 was used as the calibrator to calculate the fold change in expression of other samples. The normalized quotient at 12 months was divided by the quotient at 0 months to determine the fold change. For each experimental sample, the amplification curve for rhβ-Actin and rhMAOA fell within the detectable range of the control titration curves. Two sets of triplicate assays were performed for each sample, and two investigators reviewed the results.
CSF content of monoamines and metabolites
The NAc core was unavailable to analyze for tissue content of monoamines and their metabolites. As the NAc core is one of the main targets of DA and 5HT within the brain, and the CSF acts as a temporary reservoir of monoamines and their metabolites40, we used the CSF to investigate the effects of long-term (12 months) alcohol consumption on monoamine levels in the brain. CSF samples were collected at the same time point as the NAc core tissue, after 12 months of open access to alcohol (ethanol subjects) or water only (control subjects). No CSF sample was available from early time points for comparison.
The CSF content of monoamines and their metabolites was measured as previously described47 with minor modifications. Briefly, the CSF was centrifuged and the supernatants were analyzed for DA, 5HT and their metabolites 3,4-dihydroxyphenylacetic acid (DOPAC), homovanillic acid (HVA) and 5-hydroxyindoleacetic acid (5HIAA) using high-performance liquid chromatography (HPLC) coupled to electrochemical detection at +220 mV (ESA Inc., Chelmsford, MA) and separated on a Luna 100 × 3.0 mm C18 3 μm HPLC column (Phenomenex, Torrance, CA). The mobile phase consisted of 50 mM citric acid, 90 mM sodium dihydrogen phosphate, 1.7–2.0 mM 1-octanesulfonic acid, 50 μM EDTA, 10–12% acetonitrile and 0.3% triethylamine in a volume of 1 L (pH 3.0). Analytes were quantified using Chromeleon v 6.8 software (eDAQ Inc, Colorado Spring, CO) and a calibration curve. Samples were assayed in triplicate and were reviewed by two investigators.
Statistical Analysis
The exploratory nature of the present study limited our ability to estimate a priori the sample size needed to detect the effects investigated in this study. Therefore, we used both a longitudinal study design (blood) and a group comparison based on alcohol consumption level (NAc core, CSF) to maximize the opportunity to detect significant effects. Not all tissues were available from all study subjects at all-time points, and thus the number of subjects per tissue comparison varies slightly. Only animals that had at least two tissues available (blood, NAc or CSF) were included in this study. The maximum number of samples available was used for each analysis.
The Shapiro-Wilk test (appropriate for small sample sizes) was used to assess the normality of the continuous variables (Ave. Ethanol intake, global CpG methylation (Blood 0m, 12m, NAc core 12m), mRNA expression levels (Blood 0m, 12m, NAc core 12m), DA, DOPAC, HVA, 5HT, 5HIAA, DOPAC/DA, HVA/DA, 5HIAA/5HT). All of them followed a normal distribution.
Repeated measures of ANOVA was used to compare longitudinal changes in the relative rhMAOA expression in blood. The ethanol group (LMD, HD) was used as between group factor, and the time point (0m, 12m) as within group factor. Unstructured covariance was assumed to be a within subject covariance/correlation structure. Tukey-Kramer correction for the multiple comparisons was used to correct the overall type I error rate. In addition, the mRNA expression change between 0 and 12 months was compared using a contrast test.
One-way ANOVA was used compare the difference in the NAc core mRNA relative expression and the 5HIAA/5HT between controls and ethanol drinkers (LMD and HD). Prior to applying one-way ANOVA, Levene’s test was used to test homogeneous variance assumption for parametric methods (NAc core mRNA: p=0.1492; 5HIAA/5HT: p=0.3911). Because of the heterogeneous variance difference between controls and ethanol drinkers (LMD and HD; p=0.0014, Levene’s test for equal variance) we used the nonparametric Kruskal-Wallis test to assess the DOPAC/DA levels between controls and low-moderate alcohol drinkers (LMD) and heavy alcohol drinkers (HD). Bonferroni correction for the multiple comparisons was used to correct the overall type I error rate.
The unpaired Student’s t-test for independent samples was used to compare the NAc core rhMAOA expression levels between controls and ethanol drinkers.
Pearson correlation analysis was used to explore associations between rhMAOA mRNA (blood and NAc core), global average methylation (blood and NAc core), average ethanol or water (g/kg/day) consumption and CSF neurotransmitter levels.
A multiple regression model with stepwise selection was used to test whether the inclusion of other independent variables (age (continuous), MAOA-LPR genotype (discrete)) along with regional methylation values improved the accuracy in predicting values for the dependent variables (rhMAOA mRNA levels, average daily ethanol intake).
All analyses were performed using IBM SPSS Statistics 22 (New York, NY) and SAS 9.4 (Cary, NC) with two sided p-values <0.05.
RESULTS
Ethanol drinking patterns
As previously described36, subjects are classified as HD if they consume at least 3 g/kg/day for 20% or more of the open-access days; they are defined as LMD if they do not meet the 20% day threshold (Supplementary Figure 1A). The overall average g/kg/day for heavy drinking subjects ranges from 2.25 g/kg/day to 4.4g/kg/day (Supplementary Figure 1B). Both groups showed distinct ethanol consumption patterns over the 12 months of open access. The LMD group consumed an average of 1.5 g/kg/day during the first month of open access, and increased intake modestly over the 12-month period to 2.25 g/kg/day. In contrast the HD males began at a higher consumption rate at an average of 2.04 g/kg/day, and increased to 3.42 g/kg/day by the end of the study (Supplementary Figure 2). To evaluate whether the range of drinking behaviors reflect general differences in thirst, we compared water and ethanol intake levels. The results demonstrate that ethanol consumption was independent of water consumption level (r=0.105, p=0.633, n=23; Supplementary Figure 3A).
Blood rhMAOA expression levels before and after alcohol consumption
To determine the relationship between rhMAOA expression and risk for hazardous alcohol use, we first evaluated rhMAOA expression in the blood of alcohol naïve subjects, and then recorded the amount of alcohol they subsequently consumed over 12 months of open access to alcohol. There was a significant negative association between pre-alcohol blood rhMAOA mRNA level and the average amount of ethanol they later consumed (r= −0.589, p=0.016, n=14; Figure 1A). Interestingly, the predictive potential of the rhMAOA mRNA levels was observed as early as the 1st month of ethanol access (r= −0.586, p=0.017; Supplementary Table 3) and was stable after 7 months.
The rhMAOA mRNA levels were not associated with the amount of water consumed (r= −0.076, p=0.778, n=16; Supplementary Figure 3B). In addition, the rhMAOA expression and the average alcohol consumption levels were not associated with the age of the subjects (mRNA0m: p(Age)=0.5042, Alcohol12m: p(Age)=0.2416; multiple regression with stepwise selection).
To determine if rhMAOA expression changed following chronic alcohol use, blood from the same subjects was collected after the 12-month open access period. A significant reduction in the rhMAOA mRNA was detected after alcohol use in both groups of ethanol drinkers (LMD: Adjusted p<0.0001, HD: Adjusted p=0.0200; repeated measures ANOVA with Tukey-Kramer correction for the multiple comparisons; Figure 1C). In addition, the decrease in rhMAOA expression was negatively associated with the amount of alcohol consumed (r= −0.679, p=0.008, n=14; Figure 1B). Lighter drinkers that initially had higher levels of mRNA experienced a greater decrease in gene expression (LMD: −0.049 vs HD: −0.019, p=0.0015, repeated measures ANOVA with contrast test), and consequently, all subjects measured similar expression levels after 12 months of alcohol use (LMD: 0.973 vs HD 0.984: Adjusted p=0.6095, repeated measures ANOVA with Tukey-Kramer correction for the multiple comparisons).
rhMAOA expression in the nucleus accumbens core after chronic alcohol use
Based on the significant relationship between rhMAOA expression in the blood and alcohol intake, we next investigated rhMAOA expression in the NAc core, a central component in the reward system. This neural tissue is only accessible at the terminal time point, and thus comparisons were made between rhMAOA gene expression in control, LMD and HD subjects following 12 months of open access to alcohol or water. There was a significant negative association between the level of rhMAOA expression and the amount of ethanol consumed (r= −0.637, p=0.014, n=14; Figure 2A), even after accounting for age (p(Age)=0.3122; multiple regression with stepwise selection). The rhMAOA expression levels were significantly lower in alcohol drinkers as compared to the alcohol naïve control subjects (p=0.023, ncontrol= 10, nethanol= 14; independent t-test). Interestingly, our results indicate that these low levels of expression were driven by the HD subjects (p=0.008, nLMD=7, nHD=7; one-way ANOVA with Bonferroni correction for multiple comparisons, Figure 2B), while there was no significant difference in expression between HD and LMD groups (p=0.122; one-way ANOVA with Bonferroni correction for multiple comparisons, Figure 2B).
Effect of the LPR genotype and promoter methylation on rhMAOA expression and ethanol use
To investigate potential mechanisms regulating rhMAOA expression we considered the contribution of both rhMAOA-LPR genotype and the methylation status within the promoter. Both the 7 copy (“low” activity) and 5–6 copy (“high” activity) alleles were represented with similar frequency (nLPR_7= 12, nLPR_5/6= 16) in the study. The rhMAOA-LPR genotype was not associated with blood and NAc core rhMAOA expression level (Blood0m: p=0.3559; Blood12m: p=0.3198; NAc core: p=0.0505; multiple regression with stepwise selection).
Recent evidence suggested that the CpG methylation levels of a human MAOA promoter region may contribute to the gene regulation26, therefore we focused on the orthologous region of the CGI_1060 rhMAOA promoter in rhesus macaques. Eleven of the 12 CpGs within the region are conserved in humans and rhesus macaques (Supplementary Figure 4). The relative NHP blood methylation levels were very similar to that reported for humans26, while they were lower in the NAc core (Supplementary Figure 4). Neither the global methylation levels of the rhMAOA promoter in the blood or NAc core were significantly associated with ethanol consumption level (Supplementary Table 4). The methylation rate in blood and NAc core was also not associated with mRNA expression in the corresponding tissues (Supplementary Table 4).
Monoamine levels in the cerebrospinal fluid after long-term ethanol consumption
Given the central role of MAOA in catalyzing monoamine metabolism, we considered whether low rhMAOA expression was also correlated with elevated monoamines. Thus monoamines and their metabolites were measured in CSF collected in the rhesus macaques after 12 months of open alcohol or water access. CSF was not collected at earlier time points. A significant positive correlation was detected between the average ethanol consumption and the levels of DA (r=0.557, p=0.038, n=14; Figure 3A). The metabolic rate of DA, measured by the DOPAC/DA ratio, was correlated with the amount of alcohol consumed (r= −0.670, p=0.012, n= 13; Supplementary Table 5) and the NAc core rhMAOA expression level (r=0.722, p=0.043, n=8, Supplementary Table 5). Thus, the more alcohol consumed, the less rhMAOA mRNA expression, and the less DA was metabolized into DOPAC. In addition, the amount of ethanol consumed established a clear separation between the DOPAC/DA levels (Adjusted p(LMD vs HD)= 0.004, Kruskal-Wallis with Bonferroni correction for multiple comparisons; Figure 3B). Interestingly, the DOPAC/DA levels and the DA levels in the control subjects spanned the full range of values observed in the alcohol drinkers (Figure 3B and Supplementary Figure 5, respectively). These results suggest the pre-alcohol existence of variability in DA levels within the population.
Contrary to the DA behavior, there was a negative correlation between the average ethanol consumption and the levels of 5HT (r= −0.708, p=0.033, n=9; Figure 3C) and a positive correlation between the ethanol and 5HIAA levels (r=0.514, p=0.042, n=16; Supplementary Table 5). There were no significant differences in the 5HT metabolic rate (5HIAA/5HT) among ethanol groups (p=0.05, one-way ANOVA with Bonferroni correction for multiple comparisons; Figure 3D).
DISCUSSION
The historic link between MAOA and risk for neuropsychiatric disorders, and the variability in MAOA activity recently described in humans25, 26 led us to investigate its potential contribution to hazardous alcohol use. Our results identified, for the first time, a predictive relationship between blood rhMAOA mRNA level in alcohol-naïve individuals and the subsequent amount of alcohol consumed over the course of a 12-month self-administration study. Individuals that had lower levels of rhMAOA expression in their blood prior to alcohol access went on to drink larger amounts over the next 12 months. Moreover, the chronic alcohol consumption caused a significant reduction in rhMAOA expression in blood that was alcohol-dose dependent, with a smaller change in expression measured in the heavier drinkers. The apparent rhMAOA expression plateau reached by both low and heavy drinking subjects after long-term alcohol might reflect a minimum, or a “floor effect”, since lower expression levels could have been detected. Taken together, these findings suggest that rhMAOA levels in alcohol-naïve individuals may be an informative biomarker to identify risk for alcohol abuse48, 49. However, long-term alcohol consumption effectively establishes similar expression levels among all subjects. Additional studies are needed to understand how rapidly rhMAOA expression changes with alcohol use, and whether mRNA levels rebound during periods of abstinence.
It is notable that the rhMAOA expression levels detected in the blood prior to alcohol exposure are similar to the expression patterns measured in the NAc core after 12 months of alcohol use; rhMAOA levels in the NAc core are lowest among the heaviest alcohol drinkers. Although we do not know the pre-alcohol levels of rhMAOA mRNA in the NAc core we speculate that, as in blood, the pre-alcohol rhMAOA expression levels were lower amongst the individuals that subsequently consumed the most alcohol. Supporting this hypothesis, the rhMAOA expression levels measured in the control subjects spanned the full range of expression values detected in the LMD and some of the HD subjects (Figure 2B). As with blood, the NAc rhMAOA expression may also decrease in response to chronic alcohol use. Consistent with that notion, some of the HD subjects have lower mRNA levels than those detected in the controls.
It is known that alcohol produces its rewarding/reinforcing effects by increasing extracellular DA in the NAc50. The low levels of rhMAOA mRNA present in the NAc of heavy drinking subjects are in agreement with reduced DA metabolism contributing to the maintenance of high DA levels. Accordingly, the low levels of MAOA may enhance the sensitivity of the NAc to the reinforcing effects of alcohol51. Thus, in alcohol vulnerable subjects alcohol encounters the perfect molecular scenario to mount a potent reward response, promoting alcohol dependence. In contrast, subjects with higher amounts of MAOA (LMD) are able to regulate the surge of alcohol-induced DA and trigger a moderate reinforcement response (Figure 4). A recent study showed increases in DA release in response to alcohol occurring preferentially in human social-drinkers at risk for alcohol dependence (based on family-history of alcohol use disorders)52. The authors suggest that this response reflects a neurobiological vulnerability that alters behavioral responses toward alcohol and other rewards in high-risk individuals52. The range of DA detected in control subjects provides evidence of preexisting variability in the levels of DA in the population. In addition, our finding of variability in rhMAOA expression levels in alcohol-naïve subjects suggests a molecular mechanism contributing to the initial variability in DA levels and alcohol-vulnerability (Figure 4).
In contrast to DA, the CSF level of 5HT was negatively correlated with the amount of alcohol consumed. Our results are in agreement with previous reports of decreased extracellular 5HT in alcoholic individuals53, 54. In addition, our finding suggests that the net balance of 5HT amongst heavy drinkers may be more prominently influenced by factors other than MAOA, such as tissue-specific decreases in 5HT biosynthesis55, decreased serotonin transporter density56, 57 or degeneration of serotonergic fibers58. We did not detect a significant difference in 5HT metabolic rate between low and heavy alcoholic subjects, perhaps due to the limited number of LMD subjects for which both 5HT and 5HIAA levels were measurable. In that light, it would be interesting to evaluate additional subjects in future studies.
Our study also provides the first direct analysis of MAOA-LPR genotype and in vivo gene expression level. Although gene-reporter studies of the MAOA-LPR alleles in cell culture detected a relationship between MAOA-LPR genotype and gene expression14, 59–61, our in vivo studies do not support those conclusions. While it is possible that MAOA-LPR genotype influences MAOA expression in tissues other than blood and brain, our studies are consistent with recent reports finding a lack of genotype association25, 26.
Although we did not detect a significant association between rhMAOA promoter methylation and gene expression level in the blood or NAc core, we cannot exclude the possibility that DNA methylation in other regulatory regions contribute to MAOA transcriptional regulation. It will also be of interest to determine how CpG methylation varies among female subjects. A previous report that promoter methylation was associated with alcohol dependence in females16 and that MAOA-LPR genotype was only associated with increased alcoholism among women that also had a history of abuse62 suggests that epigenetic regulation could be a contributing factor in females.
In summary, the results obtained in this study provide novel evidence for the role of rhMAOA in alcohol vulnerability and dependence. First, we show that blood rhMAOA expression level in alcohol-naïve subjects predicts risk for alcohol consumption. In addition, our results support the proposition that rhMAOA is responsive to environmental stressors, such as chronic alcohol exposure. Finally, our findings suggest that a decrease in rhMAOA expression, and an ensuing increase in DA metabolism, may provide a feed-forward mechanism linking this key enzyme with the alcohol addiction process. Future studies will be needed to clarify whether the alcohol-associated changes in rhMAOA expression are reversed during periods of abstinence.
Supplementary Material
rhMAOA (Monoamine oxidase A) mRNA expression rates in blood before and after 12 months of alcohol consumption. A) Correlation analysis of the average daily amount of ethanol (g/kg/day) consumed over 12 months of self-administration and the blood mRNA levels prior to ethanol exposure. B) Correlation analysis of the average daily amount of ethanol (g/kg/day) consumed over 12 months of self-administration and the difference in mRNA levels before and after ethanol exposure. C) Comparison of the total average rhMAOA expression at 0 and 12 months of alcohol use between low-moderate drinkers (LMD) and heavy drinkers (HD). Whiskers represent standard error of means. Repeated measures of ANOVA with Tukey-Kramer correction for multiple comparisons, **** Adjusted p<0.0001; * Adjusted p<0.05.

rhMAOA (Monoamine oxidase A) mRNA expression rates in the nucleus accumbens (NAc) core after 12 months of alcohol consumption. A) Correlation analysis of the average daily amount of ethanol (g/kg/day) consumed over 12 months of self-administration and NAc core mRNA level after ethanol exposure. B) Comparison of the average mRNA level between control subjects, low-moderate drinkers (LMD) and heavy drinkers (HD). Middle line represents the mean. One-way ANOVA with Bonferroni correction for multiple comparisons; ** p<0.01.

The CSF levels of monoamines after 12 months of alcohol consumption. A) Correlation analysis of the average dopamine (DA) measured in the CSF versus average g/kg/day ethanol consumed during 12 months of self-administration. B) Distribution of the DA metabolic activity measured in CSF in control subjects, low-moderate drinkers (LMD) and heavy drinkers (HD). Middle line represents the mean. Kruskal-Wallis test, with Bonferroni correction for multiple comparisons; ** p<0.01. C) Correlation analysis of the average serotonin (5HT) measured in the CSF versus average g/kg/day ethanol consumed during 12 months of self-administration. D) Distribution of 5HT metabolic activity measured in CSF in control subjects, LMD and HD. Middle line represents the mean. One-way ANOVA with Bonferroni correction for multiple comparisons. ** p<0.01.

Model for the contributing role of MAOA to alcohol vulnerability (alcohol naive) and dependence (chronic alcohol). Subjects with increased vulnerability for heavy alcohol consumption initially have lower levels of MAOA and higher levels of dopamine (DA). Following chronic alcohol consumption, MAOA levels decrease further, resulting in further elevation of DA and a more potent reward response. In contrast, individuals with lower risk for alcohol consumption, initially have higher levels of MAOA, and following chronic alcohol consumption, have a more modest increase in DA.

Footnotes
AUTHORS CONTRIBUTION
BF, SRJ, KAG and RCJ were responsible for the study concept and design. KAG provided the animal samples. CH performed the tissue necropsies. SG provided the alcohol and water intake data. RCJ performed DNA, RNA isolation, real-time PCR analysis and bisulfite sequencing. RL performed the initial bisulfite pyrosequencing analysis. LJW performed the informatics analysis of the bisulfite data. BP contributed statistical analysis. RCJ, JL, LJW, RL and BF assisted with data analysis and interpretation of findings. RCJ and BF drafted the manuscript. SRJ, RL, GW, KAG and BP provided critical revision of the manuscript for important intellectual content. All authors critically reviewed content and approved the final version for publication.
FINANCIAL DISCLOSURE
All of the authors declare they do not have conflicts of interest.
Supplementary information is available at Molecular Psychiatry’s Website.
Acknowledgments
This study was supported by NIH grants U01AA020928 (BF), U01AA013510 (KAG), U01AA020890 (GW), U01AA014091 (SRJ), R24AA019431 (KAG), OD011092 (ONPRC) and the Margaret Ann Price Investigator Fund (RL), James Wah Mood Disorders Scholar Fund via Charles T. Bauer Foundation (RL) and The Kenneth Lattman Foundation (GW). We are grateful to the Oregon National Primate Research Center’s Department of Comparative Medicine and Primate Genetics Core Services for their assistance in this project. The authors thank Dr. Julie Hollister-Smith for her technical support and helpful comments. We also appreciate the very helpful critiques provided by the anonymous reviewers.
References
- 1. National Institute of Alcohol Abuse and Addiction2014Alcohol facts and statisticsAvailable at: www.niaaa.nih.gov/alcohol-health/overview-alcohol-consumption/alcohol-facts-and-statistics
- 2. Genetics and alcoholismNature reviews Gastroenterology & hepatology2013108487494[PubMed][Google Scholar]
- 3. Genetic influences on the development of alcoholismCurrent psychiatry reports20131511412[PubMed][Google Scholar]
- 4. Interaction between early life stress and alcohol dependence on neural stress reactivityAddiction biology2014[Google Scholar]
- 5. Interactive effects of childhood maltreatment and recent stressful life events on alcohol consumption in adulthoodJournal of studies on alcohol and drugs2012734559569[PubMed][Google Scholar]
- 6. New insights on neurobiological mechanisms underlying alcohol addictionNeuropharmacology201367223232[PubMed][Google Scholar]
- 7. The neurobiology of alcohol consumption and alcoholism: an integrative historyPharmacology, biochemistry, and behavior20131132037[Google Scholar]
- 8. How adaptation of the brain to alcohol leads to dependence: a pharmacological perspectiveAlcohol research & health: the journal of the National Institute on Alcohol Abuse and Alcoholism2008314310339[PubMed][Google Scholar]
- 9. Regulation of nucleus accumbens dopamine release by the dorsal raphe nucleus in the ratNeurochemical research1992175401407[PubMed][Google Scholar]
- 10. Neurocircuitry of addictionNeuropsychopharmacology: official publication of the American College of Neuropsychopharmacology2010351217238[PubMed][Google Scholar]
- 11. Drugs abused by humans preferentially increase synaptic dopamine concentrations in the mesolimbic system of freely moving ratsProceedings of the National Academy of Sciences of the United States of America1988851452745278[PubMed][Google Scholar]
- 12. Serotonin2C receptors in the nucleus accumbens are involved in enhanced alcohol-drinking behaviorThe European journal of neuroscience201235813681380[PubMed][Google Scholar]
- 13. Microdialysis of dopamine in the nucleus accumbens of alcohol-preferring (P) rats during anticipation and operant self-administration of ethanolAlcoholism, clinical and experimental research2002263318325[Google Scholar]
- 14. A functional polymorphism in the monoamine oxidase A gene promoterHuman genetics19981033273279[PubMed][Google Scholar]
- 15. Neural mechanisms of genetic risk for impulsivity and violence in humansProceedings of the National Academy of Sciences of the United States of America20061031662696274[PubMed][Google Scholar]
- 16. MAOA methylation is associated with nicotine and alcohol dependence in womenAmerican journal of medical genetics Part B, Neuropsychiatric genetics: the official publication of the International Society of Psychiatric Genetics2008147B5565570[Google Scholar]
- 17. Monoamine oxidase in neuropsychiatry and behaviorAmerican journal of human genetics1999653593598[PubMed][Google Scholar]
- 18. MAOA alters the effects of heavy drinking and childhood physical abuse on risk for severe impulsive acts of violence among alcoholic violent offendersAlcoholism, clinical and experimental research2010345853860[Google Scholar]
- 19. Role of genotype in the cycle of violence in maltreated childrenScience20022975582851854[PubMed][Google Scholar]
- 20. MAOA, childhood maltreatment, and antisocial behavior: meta-analysis of a gene-environment interactionBiological psychiatry2014751917[PubMed][Google Scholar]
- 21. MAOA genotype, childhood maltreatment, and their interaction in the etiology of adult antisocial behaviorsBiological psychiatry20147512530[PubMed][Google Scholar]
- 22. Meta-analysis of gene-environment interactions in developmental psychopathologyDevelopment and psychopathology200719410291037[PubMed][Google Scholar]
- 23. Meta-analysis of six genes (BDNF, DRD1, DRD3, DRD4, GRIN2B and MAOA) involved in neuroplasticity and the risk for alcohol dependenceDrug and alcohol dependence2015149259263[PubMed][Google Scholar]
- 24. Investigation of the functional effect of monoamine oxidase polymorphisms in human brainHuman genetics2002110117[PubMed][Google Scholar]
- 25. Evidence that brain MAO A activity does not correspond to MAO A genotype in healthy male subjectsBiological psychiatry2007624355358[PubMed][Google Scholar]
- 26. Evidence that the methylation state of the monoamine oxidase A (MAOA) gene predicts brain activity of MAO A enzyme in healthy menEpigenetics: official journal of the DNA Methylation Society201271011511160[Google Scholar]
- 27. Conserved role of intragenic DNA methylation in regulating alternative promotersNature20104667303253257[PubMed][Google Scholar]
- 28. Epigenetic regulation of the glucocorticoid receptor in human brain associates with childhood abuseNature neuroscience2009123342348[PubMed][Google Scholar]
- 29. Annual Research Review: Epigenetic mechanisms and environmental shaping of the brain during sensitive periods of developmentJournal of child psychology and psychiatry, and allied disciplines2011524398408[Google Scholar]
- 30. Role of CpG islands in the up-regulation of NMDA receptor NR2B gene expression following chronic ethanol treatment of cultured cortical neurons of miceNeurochemistry international2005464313327[PubMed][Google Scholar]
- 31. The site specific demethylation in the 5′-regulatory area of NMDA receptor 2B subunit gene associated with CIE-induced up-regulation of transcriptionPloS one201051e8798[PubMed][Google Scholar]
- 32. Epigenetics-beyond the genome in alcoholismAlcohol research: current reviews2012343293305[PubMed][Google Scholar]
- 33. DNA methylation of the Monoamine Oxidases A and B genes in postmortem brains of subjects with schizophreniaOpen Journal of Psychiatry201224A374383[Google Scholar]
- 34. Drinking typography established by scheduled induction predicts chronic heavy drinking in a monkey model of ethanol self-administrationAlcoholism, clinical and experimental research2008321018241838[Google Scholar]
- 35. The effects of age at the onset of drinking to intoxication and chronic ethanol self-administration in male rhesus macaquesPsychopharmacology2014[Google Scholar]
- 36. Chronic alcohol self-administration in monkeys shows long-term quantity/frequency categorical stabilityAlcoholism, clinical and experimental research2014381128352843[Google Scholar]
- 37. Non-human primate models of alcohol-related phenotypes: the influence of genetic and environmental factorsCurrent topics in behavioral neurosciences201313223249[PubMed][Google Scholar]
- 38. Advances in nonhuman primate alcohol abuse and alcoholism researchPharmacology & therapeutics20031003235255[PubMed][Google Scholar]
- 39. Non-human primate models of inheritance vulnerability to alcohol use disordersAddiction biology2006113–4374385[PubMed][Google Scholar]
- 40. Ventricular cerebrospinal fluid monoamine transmitter and metabolite concentrations reflect human brain neurochemistry in autopsy casesJournal of neurochemistry199054411481156[PubMed][Google Scholar]
- 41. Standardized method for the harvest of nonhuman primate tissue optimized for multiple modes of analysesCell and tissue banking201415199110[PubMed][Google Scholar]
- 42. Genetic load is associated with hypothalamic-pituitary-adrenal axis dysregulation in macaquesGenes, brain, and behavior2012[Google Scholar]
- 43. Correction of PCR-bias in quantitative DNA methylation studies by means of cubic polynomial regressionNucleic acids research20113911e77[PubMed][Google Scholar]
- 45. Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applicationsBioinformatics2011271115711572[PubMed][Google Scholar]
- 46. BSmooth: from whole genome bisulfite sequencing reads to differentially methylated regionsGenome biology20121310R83[PubMed][Google Scholar]
- 47. Methylphenidate and cocaine self-administration produce distinct dopamine terminal alterationsAddiction biology2014192145155[PubMed][Google Scholar]
- 48. Identification of blood biomarkers for psychosis using convergent functional genomicsMolecular psychiatry20111613758[PubMed][Google Scholar]
- 49. Identifying blood biomarkers for mood disorders using convergent functional genomicsMolecular psychiatry2009142156174[PubMed][Google Scholar]
- 50. The mesolimbic dopamine system: the final common pathway for the reinforcing effect of drugs of abuse?Neuroscience and biobehavioral reviews2006302215238[PubMed][Google Scholar]
- 51. Ethanol is self-administered into the nucleus accumbens shell, but not the core: evidence of genetic sensitivityAlcoholism, clinical and experimental research2009331221622171[Google Scholar]
- 52. Differential striatal dopamine responses following oral alcohol in individuals at varying risk for dependenceAlcoholism, clinical and experimental research2014381126134[Google Scholar]
- 53. 5-Hydroxyindoleacetic acid in cerebrospinal fluid in alcoholic patients under different clinical conditionsAlcohol198523415418[PubMed][Google Scholar]
- 54. Ethanol impairs tryptophan transport into the brain and depresses serotoninLife sciences1981292627512755[PubMed][Google Scholar]
- 55. Alterations in brain serotonin synthesis in male alcoholics measured using positron emission tomographyAlcoholism, clinical and experimental research2009332233239[Google Scholar]
- 56. The effects of chronic ethanol self-administration on hippocampal 5-HT1A receptors in monkeysDrug and alcohol dependence2014136135142[PubMed][Google Scholar]
- 57. Serotonin transporter distribution and density in the cerebral cortex of alcoholic and nonalcoholic comparison subjects: a whole-hemisphere autoradiography studyThe American journal of psychiatry20021594599606[PubMed][Google Scholar]
- 58. Brainstem serotonergic neurons in chronic alcoholics with and without the memory impairment of Korsakoff’s psychosisJournal of neuropathology and experimental neurology1993526567579[PubMed][Google Scholar]
- 59. Excess of high activity monoamine oxidase A gene promoter alleles in female patients with panic disorderHuman molecular genetics199984621624[PubMed][Google Scholar]
- 60. Association between monoamine oxidase A activity in human male skin fibroblasts and genotype of the MAOA promoter-associated variable number tandem repeatHuman genetics19991056542551[PubMed][Google Scholar]
- 61. Monoamine oxidase A gene promoter variation and rearing experience influences aggressive behavior in rhesus monkeysBiological psychiatry2005572167172[PubMed][Google Scholar]
- 62. Interaction between a functional MAOA locus and childhood sexual abuse predicts alcoholism and antisocial personality disorder in adult womenMolecular psychiatry2008133334347[PubMed][Google Scholar]