flexibility that would have been an evolution-
ary advantage for Anthropoidea during the
major environmental perturbations of the later
Paleogene, including the global cooling event at
the Eocene-Oligocene boundary.
REFERENCES AND NOTES
- R. Hoffstetter, R. Lavocat,C. R. Acad. Sci. (Paris) 271 , 172– 175
(1970). - J. G. Fleagle,Primate Adaptation and Evolution(Academic
Press, ed. 3, 2013). - A. Houle,Am. J. Phys. Anthropol. 109 , 541–559 (1999).
- J. G. Fleagle, C. C. Gilbert, inPrimate Biogeography: Progress
and Prospects, S. M. Lehman, J. G. Fleagle, Eds. (Springer,
2006), pp. 375–418. - M. F. Tejedor, thesis, Universidad Nacional de la Plata (2000).
- M. Bondet al.,Nature 520 , 538–541 (2015).
- P.-O. Antoineet al.,Proc. Biol. Sci. 279 , 1319–1326 (2012).
- E. R. Seiffertet al.,Hist. Biol. 30 , 204–226 (2018).
- K. E. Campbell Jr.,Nat. Hist. Mus. Los Angeles County Sci. Ser.
40 ,i–vi, 1–163 (2004). - M. Arnal, A. G. Kramarz, M. G. Vucetich, C. D. Frailey,
K. E. Campbell Jr.,Pap. Palaeontol.10.1002/spp2.1264
(2019). - C. D. Frailey, K. E. Campbell Jr.,Nat. Hist. Mus. Los Angeles
County Sci. Ser. 40 , 71–130 (2004). - F. J. Goin, A. M. Candela,Nat. Hist. Mus. Los Angeles County
Sci. Ser. 40 , 15–60 (2004). - E. L. Simons, R. F. Kay,Nature 304 , 624–626 (1983).
- E. L. Simons, E. R. Seiffert, P. S. Chatrath, Y. Attia,Folia
Primatol. 72 , 316–331 (2001). - K. C. Beard, P. M. C. Coster,Am. J. Phys. Anthropol. 159 ,
714 – 721 (2016). - N. Egi, M. Takai, N. Shigehara, T. Tsubamoto,Int. J. Primatol. 25 ,
211 – 236 (2004). - D. Silvestroet al.,Syst. Biol. 68 , 78–92 (2019).
- G. F. Gunnellet al.,Nat. Commun. 9 , 3193 (2018).
- F. Ronquistet al.,Syst. Biol. 61 , 539–542 (2012).
- E. L. Simons, R. F. Kay,Am. J. Primatol. 15 , 337–347 (1988).
- J. G. Fleagle, R. F. Kay,J. Hum. Evol. 16 , 483–532 (1987).
- R. E. Kass, A. E. Raftery,J. Am. Stat. Assoc. 90 , 773–795 (1995).
- K. G. Milleret al.,Geol. Soc. Am. Bull. 120 , 34–53 (2008).
- J. G. Fleagle, E. L. Simons,Am. J. Phys. Anthropol. 97 , 235– 289
(1995). - M. F. Teaford, M. C. Maas, E. L. Simons,Am. J. Phys. Anthropol.
101 , 527–543 (1996). - E. C. Kirk, E. L. Simons,J. Hum. Evol. 40 , 203–229 (2001).
- E. R. Seiffertet al.,Science 310 , 300–304 (2005).
- N. J. Stevens, E. R. Seiffert, E. M. Roberts, P. M. O’Connor,
Am. J. Phys. Anthropol. 156 (suppl.), 296–297 (2015). - L. Marivauxet al.,J. Hum. Evol. 97 , 159–175 (2016).
- J. I. Blochet al.,Nature 533 , 243–246 (2016).
- E. Seiffertet al., A parapithecid stem anthropoid of African
origin in the Paleogene of South America, Version 1,
Dryad (2020); https://doi.org/10.5061/dryad.0vt4b8gtt.
ACKNOWLEDGMENTS
The Instituto Geológico, Minero, y Metalúrgico (INGEMMET), Lima,
and particularly C. Chacaltana Budiel and L. Tejada Medina have
facilitated and supported paleontological and geological work at
the Santa Rosa locality since 1995. M. Vilca of INGEMMET and
S. Hynek participated in fieldwork in 2016. Members of the Santa
Rosa indigenous community assisted with fieldwork in 1998 and
- I. Mawhinney and M. Shazab helped to process the Santa
Rosa matrix. T. Jashashvili scanned theUcayalipithecusspecimens.
M.F.T. is a research associate at the Gothenburg Global Biodiversity
Centre.Funding:J. Wigmore, W. Rhodes, and R. Seaver helped
to fund the 1998 expedition that led to the recovery of the
Ucayalipithecuspartial upper molars. The Leakey Foundation,
Gordon Getty, and A. Stenger supported fieldwork in 2016. Micro–
computed tomography scanning was supported by the Keck School of
Medicine of USC and the U.S. National Science Foundation (BCS-
1231288).Author contributions:K.E.C. ran the Santa Rosa project.
K.E.C. collected matrix from Santa Rosa in 1998. J.G.F., F.M.C., and
K.E.C. collected matrix from Santa Rosa in 2016. M.F.T., N.M.N.,
and M.B. identified the upper molars ofUcayalipithecusin the collection
from the 1998 field season. E.R.S. identified the lower molars of
Ucayalipithecusin the collection from the 2016 field season. E.R.S.,
M.F.T., J.G.F., N.M.N., M.B., and K.E.C. analyzed the fossils. E.R.S.
scored the taxa for phylogenetic analysis, ran the phylogenetic and
biogeographic analyses, and created digital models. J.G.F., K.E.C.,
D.d.V., and E.R.S. picked the Santa Rosa matrix; D.d.V. created digital
models. E.R.S. wrote the first draft and created the figures. All authors
read and edited the paper.Competing interests:The authors
declare no competing interests.Data and materials availability:Data
used for the phylogenetic and biogeographic analyses are available
for download on the Dryad Digital Repository ( 31 ). Digital surface
models of theUcayalipithecusspecimens are available for download
on MorphoSource (project 872). Fossils are deposited in the
Colección Paleontológica del INGEMMET (CPI), Lima, Perú.
SUPPLEMENTARY MATERIALS
science.sciencemag.org/content/368/6487/194/suppl/DC1
Materials and Methods
Supplementary Text
5 November 2019; accepted 19 February 2020
10.1126/science.aba1135
NEUROSCIENCE
Dopaminylation of histone H3 in ventral tegmental
area regulates cocaine seeking
Ashley E. Lepack^1 , Craig T. Werner^2 , Andrew F. Stewart^1 , Sasha L. Fulton^1 , Ping Zhong^3 ,
Lorna A. Farrelly^1 , Alexander C. W. Smith^1 , Aarthi Ramakrishnan^1 , Yang Lyu^1 , Ryan M. Bastle^1 ,
Jennifer A. Martin^2 , Swarup Mitra^2 , Richard M. O’Connor^1 , Zi-Jun Wang^3 , Henrik Molina^4 ,
Gustavo Turecki^5 , Li Shen^1 , Zhen Yan^3 , Erin S. Calipari^6 , David M. Dietz^2 , Paul J. Kenny^1 , Ian Maze1,7*
Vulnerability to relapse during periods of attempted abstinence from cocaine use is
hypothesized to result from the rewiring of brain reward circuitries, particularly ventral
tegmental area (VTA) dopamine neurons. How cocaine exposures act on midbrain dopamine
neurons to precipitate addiction-relevant changes in gene expression is unclear. We found
that histone H3 glutamine 5 dopaminylation (H3Q5dop) plays a critical role in cocaine-induced
transcriptional plasticity in the midbrain. Rats undergoing withdrawal from cocaine showed
an accumulation of H3Q5dop in the VTA. By reducing H3Q5dop in the VTA during withdrawal,
we reversed cocaine-mediated gene expression changes, attenuated dopamine release in the
nucleus accumbens, and reduced cocaine-seeking behavior. These findings establish a
neurotransmission-independent role for nuclear dopamine in relapse-relatedtranscriptional
plasticity in the VTA.
C
ocaine increases dopamine neurotrans-
mission from the ventral tegmental area
(VTA) to reward-relevant brain regions.
This action is central to its addictive
properties. Non-neurotransmission roles
for dopamine in cocaine dependency have not
been considered. Our laboratory has recent-
ly described a role for serotonin (5-HT) in de-
veloping 5-HTergic neurons, whereby 5-HT
located in the nucleus of these neurons was
shown to covalently attach to histone proteins—
specifically on H3 glutamine 5 (H3Q5)—to regu-
late gene expression through a process called
serotonylation ( 1 , 2 ). We had hypothesized
that this mechanism may generalize to other
monoamines in brain, such as dopamine. If
true, this process could potentially play a role
in the addiction-relevant actions of drugs that
stimulate dopaminergic transmission.
Drug addictions are defined by pathological
drug-seeking behavior that persists despite ad-
verse consequences. Prolonged vulnerability to
relapse is hypothesized to reflect the function-
al rewiring of brain reward circuitries ( 3 , 4 ).
This is precipitated, at least in part, by drug-
induced transcriptional plasticity in midbrain
dopamine neurons ( 5 , 6 ). Histone mechanisms
that control chromatin structures, and conse-
quently gene expression, regulate addiction-
relevant behaviors ( 7 , 8 ). Given that histone
H3canbemodifiedbymonoaminesinre-
sponse to fluctuations in intracellular availa-
bility, we assessed whether dopamine, like
5-HT, can be transferred to the H3 N-terminal
tail. We performed targeted, peptide-based liq-
uid chromatography–tandem mass spectrome-
try (LC-MS/MS) after in vitro transglutaminase 2
(TGM2) ( 1 , 2 , 9 ) enzymatic assays with dopa-
mine. Peptide LC-MS/MS analyses (fig. S1, A to
D) revealed Q5 as a reactive substrate for the
dopaminyl mark [H3 glutamine 5 dopaminyla-
tion (H3Q5dop)]. Given that the serotonyl modi-
fication can exist both in isolation (H3Q5ser)
and in combination with H3 lysine 4 trimethyl-
ation (H3K4me3Q5ser), we examined the effect
of K4me3 on TGM2-mediated dopaminylation
in vitro. Unmodified versus methylated mono-
nucleosomes were subjected to TGM2 dopa-
minylation assays. Using an antibody against
H3Q5dop, vide infra, we found that TGM2
equally dopaminylates unmodified and K4me3
SCIENCEsciencemag.org 10 APRIL 2020•VOL 368 ISSUE 6487 197
(^1) Department of Neuroscience, Friedman Brain Institute,
Icahn School of Medicine at Mount Sinai, New York, NY
10029, USA.^2 Department of Pharmacology and Toxicology,
Program in Neuroscience, State University of New York at
Buffalo, Buffalo, NY 14214, USA.^3 Department of Physiology
and Biophysics, Jacobs School of Medicine and Biomedical
Sciences, State University of New York at Buffalo, Buffalo,
NY 14203, USA.^4 Proteomics Resource Center, The
Rockefeller University, New York, NY 10065, USA.
(^5) Department of Psychiatry, McGill University, Montreal, QC
H3A 1A1, Canada.^6 Department of Pharmacology, Center for
Addiction Research, Vanderbilt University School of
Medicine, Nashville, TN 37232, USA.^7 Department of
Pharmacological Sciences, Icahn School of Medicine at
Mount Sinai, New York, NY 10029, USA.
*Corresponding author. Email: [email protected]
RESEARCH | REPORTS