Cell Reports
Cell Reports
Cell Reports
IN PHENOTYPING
'HWHUPLQHWKHIXQFWLRQDOLPSDFWRI JHQHVRQPHWDEROLVPLQDQKRXU
Mitochondrial Respiration
Energetic
Aerobic
Metab
olic P
otenti
al
Stressed
Phenotype
Baseline
Phenotype
Quiescent
Glycolytic
Glycolysis
Bethyl Laboratories, Inc. has been dedicated to supporting scientific discovery through
its qualified antibody products and custom antibody services since its founding in 1972.
Every antibody that Bethyl sells has been manufactured to exacting standards at its sole
location in Montgomery, Texas, and has been validated in-house by Bethyls team of
scientists. Antibodies are tested across a range of applications including western blot,
immunoprecipitation, immunohistochemistry, immunocytochemistry, ChIP, proximity
ligation assay and ELISA.
Currently, Bethyls portfolio consists of over 7,150 catalog products; offering close to
5,700 primary antibodies targeting over 2,700 proteins and 1,450 secondary antibodies
raised againstimmunoglobulins from over 25 species. Trial sizes are available for over
4,000 antibodies targeting more than 2,350 protein targets. They are conveniently
priced at $50 and serve as an opportunity to discover for yourself why our antibodies
arereally good.
For really good antibodies, visit bethyl.com/trialsize
Terms & Conditions: $50 pricing for US customers only; international customers please
contact your distributor for details. Trial sizes (catalog # ending in -T) are not available
for all antibodies and existing promotions or discounts do not apply.
2015 Bethyl Laboratories, Inc. All rights reserved.
SAY GOOD-BYE
TO THE DARKROOM
Foreword
We are so pleased to present the Best of Cell Reports 2015. In year 4 we saw some major
changes, including moving to weekly publication. From tales of very long-lived whales
to an analysis of what might have improved last years influenza vaccine, we published
top-notch content from across the life sciences. That made it all the more difficult to
put together this issue, and there were some truly tough choices to make. With the final
content in this supplement, we try to represent the breadth of biology covered by Cell
Reports and have based our choices on reader downloads, altmetric scores, and citation
information. We would urge you (once you have perused the Best of of course!) to browse
our table of contents to explore the full range of Open Access papers we publish in your
field and beyond.
We are really looking forward to 2016. In the last year, Cell Reports has been globetrotting.
We organized the Ubiquitous Ubiquitin Signaling Lablinks in Cambridge, UK and the RNA
in the Nervous System meeting in New York, USA. In November 2015, we cohosted a
major symposium in Singapore on human genomics and organized a day-long cancer and
metabolism meeting in Philadelphia. We plan to be out and about in 2016 and will hopefully
visit your neighborhood, too!
Last, but certainly not least, we want to thank the researchers who are the authors,
reviewers, readers, and editorial and advisory board members that make Cell Reports
possible. With an expanding editorial and production team internally and editorial board
advisors externally, we are so excited about 2016. But for now, kick back, grab a coffee,
tea, or beverage, and enjoy the following selection of papers.
Finally, we are grateful for the generosity of our sponsors, who helped to make this reprint
collection possible.
Jonathan Christison
Program Director, Best of Cell Press
e: [email protected]
p: 617-397-2893
t: @cellPressBiz
Robust performance
Excellent price
A suspension of single cells from pooled C57BL/6 mouse spleen and lymph nodes was prepared to isolate CD4+ T cells using
the MojoSort Mouse CD4 T Cell Isolation Kit. Cells were stained with PE anti-mouse CD4 (clone RM4-4), APC anti-mouse CD3
(145-2C11), and 7-AAD. Grateful Dead cells were excluded from the analysis.
biolegend.com
08-0052-23
Best of 2015
Reports
Leukemia-Associated Somatic Mutations Drive
Distinct Patterns of Age-Related Clonal Hemopoiesis
(continued)
Articles
Comparative Hi-C Reveals that CTCF Underlies
Evolution of Chromosomal Domain Architecture
Resources
Mapping Social Behavior-Induced Brain Activation
at Cellular Resolution in the Mouse
IN VIVO
IMAGING
PROBES
Figure 1
Figure 2
Figure 2
Figure 1: Apoptotic cells are labeled green
using CAS-MAP probes in an in vivo study
of postnatal development. Figure 2: Apoptotic COLO205 tumor cells are labeled in
response to treatment.
rndsystems.com/
Immunoassays
International Corporation
Selective detection of
CPVKIGPURGEKE6EGNNU
%/8RQUKVKXGUCORNG
5RGEKE6%4
6GVTCOGT2'
DOWVCVKQP
6GVTCOGT2'
%&
*.#RGRVKFG%QORNGZ
%&(+6%
%QPVCEV/$.+PVGTPCVKQPCNVQFC[VQNGCTPOQTGCDQWV
QWTRTQFWEVUCPFJQYYGECPJGNR[QWFKUEQXGTYC[U
VQOCMG[QWTUEKGPEGOQTGGHEKGPVCPFGHHGEVKXG
%&(+6%
SIMPLE. EFFECTIVE.
AFFORDABLE.
CELL-BASED
VIABILITY
REPORTER
TOXICITY
METABOLISM
PROLIFERATION
OXYGEN CONSUMPTION
PROTEIN SYNTHESIS
PHAGOCYTOSIS & AUTOPHAGY
CLARIOstar
PHERAstar FSX
Omega series
SPECTROstar Nano
Upgradeable single to
multi-mode lter-based
microplate readers.
The Omega series offers
a combination of exibility
and performance for any
life science application.
EXQUISITE
SENSITIVITY.
2015 Beckman Coulter Life Sciences. All rights reserved. CytoFLEX is for research use only. Class 1 Laser Product. Not for use in diagnostic
purposes. CytoFLEX and CytExpert are trademarks of Xitogen Technologies (Suzhou), Inc., a Beckman Coulter company. Beckman Coulter,
the stylized logo and fast track to success are trademarks of Beckman Coulter, Inc. Beckman Coulter and the stylized logo are registered in
the USPTO. All other trademarks are the property of their respective owners.
Are you
still injecting?
Focus on your research instead, and let
ALZET Osmotic Pumps do the dosing for you.
ALZET pumps are a superior alternative to repetitive
injections and other dosing methods that require
frequent animal handling. These fully implantable
pumps provide continuous and precise administration,
for up to 6 weeks with a single pump, to unrestrained
lab animals as small as mice. ALZET pumps are
economical and easy to use by research personnel.
Connection to a catheter enables direct delivery to
vessels, cerebral ventricles, and other target sites.
Learn more at alzet.com.
Now available: iPRECIO Pumps
Programmable Rellable
Implantable
Small size for mice and rats
Learn more at www.alzet.com/iprecio
HUMANIZED
TM
NSG MICE
START YOUR
STUDY TODAY!
Expect delivery of your order
within 7-10 business days.
VISIT SLAS2016.ORG
FOR FULL PROGRAM
INFORMATION.
SAN DIEGO
CONVENTION CENTER
JANUARY
23-27, 2016
SLAS2016.ORG
Recycling plastic
is fantastic!
Rainin TerraRack
reduces tip rack waste.
Request your free
sample today!
mt.com/rainin-tr
TALK TO HENRIK
MADE IN SWEDEN
PrecisA Monoclonals are precise, accurate and targeted. The stringent production process and
characterization procedure ensures premium performance in approved applications, with dened
specicity, multiplexing opportunities, secured continuity and stable supply.
If youre in need of an antibody to target a protein we dont list in our online catalog you might be
interested in our antibody co-development program. Learn more today at atlasantibodies.com/talktohenrikj
Cell Reports
Report
Leukemia-Associated Somatic Mutations Drive
Distinct Patterns of Age-Related
Clonal Hemopoiesis
Thomas McKerrell,1,13 Naomi Park,2,13 Thaidy Moreno,3 Carolyn S. Grove,1 Hannes Ponstingl,1 Jonathan Stephens,4,5
Understanding Society Scientic Group,6 Charles Crawley,7 Jenny Craig,7 Mike A. Scott,7 Clare Hodkinson,4,8
Joanna Baxter,4,8 Roland Rad,9,10 Duncan R. Forsyth,11 Michael A. Quail,2 Eleftheria Zeggini,12 Willem Ouwehand,4,5,12
Ignacio Varela,3 and George S. Vassiliou1,4,7,*
1Haematological
SUMMARY
Target Codon
DNMT3A
R882
JAK2
V617
NPM1
L287
SRSF2
P95
SF3B1
K666
SF3B1
K700
IDH1
R132
IDH2
R140
IDH2
R172
KRAS
G12
NRAS
G12
NRAS
Q61
KIT
D816
FLT3
D835
FLT3
N676
The important ndings of these studies were based on analysis of exome-sequencing data sets that were generated for
the study of constitutional genomes, thus trading genomewide coverage for reduced sensitivity for detecting small subclonal events. We used the different approach of targeted resequencing of selected leukemia-associated mutation hot spots
in blood DNA from more than 4,000 individuals unselected for
blood disorders. In addition to increasing the sensitivity for detecting subclonal mutations, this approach enabled us to prospectively select and study a large number of elderly individuals.
Our results show that clonal hemopoiesis is signicantly more
common than anticipated, give new insights into the distinct
age-distribution and biological behavior of clonal hemopoiesis
driven by different mutations, and help explain the increased
incidence of myelodysplastic syndromes (MDSs) with advancing
age.
RESULTS
DISCUSSION
To investigate the incidence, target genes, and age distribution
of age-related clonal hemopoiesis (ARCH), we performed targeted re-sequencing for hot spot mutations at 15 gene loci recurrently mutated in myeloid malignancies (Table 1) using blood
DNA from 3,067 blood donors aged 1770 (Wellcome Trust
Case Control Consortium [WTCCC]) and 1,152 unselected
individuals aged 6098 years (United Kingdom Household
Longitudinal Study [UKHLS]; see Figure S1 for detailed age distributions). To do this, we developed and validated a robust
methodology, employing barcoded multiplex PCR of mutational
hot spots followed by next-generation sequencing (MiSeq) and
bioinformatic analysis, to extract read counts and allelic fractions
for reference and non-reference nucleotides. This reliably detected mutation-associated circulating blood cell clones with a
variant allele fraction (VAF) R 0.008 (0.8%; see Supplemental
Experimental Procedures and Figure S2).
1240 Cell Reports 10, 12391245, March 3, 2015 2015 The Authors
Hematological malignancies develop through the serial acquisition of somatic mutations in a process that can take many years
or even decades (Ford et al., 1998; Kyle et al., 2002). Also, it is
clear that the presence of hemopoietic cells carrying leukemiaassociated mutations is only followed by the onset of hematological malignancies in a minority of cases (Busque et al., 2012;
Genovese et al., 2014; Jacobs et al., 2012; Jaiswal et al., 2014;
Laurie et al., 2012; Xie et al., 2014). In order to understand the
incidence and clonal dynamics of pre-leukemic clonal hemopoiesis, we interrogated 15 leukemia-associated mutation hot spots
using a highly sensitive methodology able to detect small clones
with mutations.
We show that clonal hemopoiesis is rare in the young but becomes common with advancing age. In particular, we observed
that ARCH driven by the mutations studied here doubled in
Table 2. Amino Acid Consequences and VAFs of the 112 Clonal Mutations Identied in This Study
Mutation
Hot Spot
DNMT3A R882
Codon
p.R882H
VAF (%)
4.14
Age
Mutation
Hot Spot
25
Codon
VAF (%)
Age
p.R882H
32.02
81
Mutation
Hot Spot
Codon
VAF (%)
Age
IDH1 R132
p.R132H
42.13
84
p.R132C
0.92
92
p.R140Q
6.67
76
p.R882C
2.33
35
p.R882H
1.14
81
p.R882H
3.80
42
p.R882H
3.06
81
IDH2 R140
SRSF2 P95
p.R882H
4.00
42
p.R882H
2.17
81
p.R882H
1.25
43
p.R882H
1.13
82
p.P95R
4.46
70
p.P95L
3.35
72
p.R882H
19.00
48
p.R882H
1.46
82
p.P95H
0.86
73
p.R882H
1.18
49
p.R882C
2.62
82
p.P95H
0.84
77
p.R882S
1.74
49
p.R882C
6.15
89
p.P95L
0.97
79y
p.R882H
9.87
50
p.R882C
2.00
94
p.P95L
0.85
80yy
p.R882H
0.83
51
p.V617F
1.56
34
p.P95H
6.67
80yy
p.R882C
1.10
51
p.V617F
4.91
42
p.P95L
0.96
81
p.R882C
12.50
52
p.V617F
7.72
45
p.P95H
6.40
82
JAK2V617F
p.R882C
1.28
53
p.V617F
0.85
62
p.P95L
2.74
85
p.R882C
2.47
54
p.V617F
25.44
64
p.P95R
7.52
87
p.R882H
1.95
55
p.V617F
7.41
65
p.P95L
5.84
88**
p.R882C
30.22
55
p.V617F
1.03
67
p.P95H
10.48
88**
p.R882C
1.22
56
p.V617F
0.88
71
p.P95R
2.71
88
p.R882H
0.91
58
p.V617F
3.75
71
p.P95R
17.05
90z
p.R882H
4.17
60
p.V617F
1.16
75
p.K700E
1.04
76
SF3B1 K700
p.R882H
5.90
60
p.V617F
2.30
77
p.K700E
6.63
81
p.R882H
9.60
60
p.V617F
1.92
78
p.K700E
0.79
82
p.R882H
2.73
60
p.V617F
2.26
80*
p.K700E
12.59
p.R882C
9.33
60
p.V617F
4.25
80
p.K700E
8.77
83
83zz
p.R882H
7.03
61
p.V617F
1.92
80
p.K700E
1.02
84
p.R882C
1.21
61
p.V617F
3.71
80
p.K700E
0.85
90z
p.R882H
0.86
63
p.V617F
15.48
81
p.R882H
2.54
64
p.V617F
1.21
82
SF3B1 K666
p.K700E
1.37
90
p.K666N
1.33
70
p.R882H
3.19
67
p.V617F
1.62
85
p.K666N
5.01
79
p.R882H
2.74
70
p.V617F
0.83
85
p.K666N
13.36
79y
p.R882H
4.27
74
p.V617F
1.98
86
p.K666N
15.43
80*
p.R882H
0.85
74
p.V617F
25.94
88
p.K666N
4.60
81
p.R882H
0.85
75
p.V617F
10.88
88**
p.K666E
1.09
83zz
p.R882C
1.12
77
p.V617F
2.94
90
p.K666N
35.11
86
p.R882C
1.15
78
p.V617F
1.23
90
p.K666N
19.70
86
p.R882H
1.26
79
p.R882H
16.66
80
p.R882C
4.28
80
p.R882C
3.66
80
KRAS G12
NRAS G12
p.G12 R
0.94
55
p.K666N
16.55
86
p.G12S
2.78
78
p.K666E
3.34
95
p.G12S
1.50
61
p.G12D
0.96
62
Mutations identied in the same sample are highlighted with the same symbol (*, **, y, yy, z, and zz).
bers et al., 2007; Rossi et al., 2005) and humans (Rube et al.,
2011; Taraldsrud et al., 2009), and it is also clear that aging
has a profound effect on the hemopoietic niche, reducing its
ability to sustain polyclonal hemopoiesis, favoring oligo- or
monoclonality instead (Vas et al., 2012). These and many other
observations provide strong evidence that changes in the hemopoietic system subject HSCs to changing pressures during
normal aging, driving clonal selection (Rossi et al., 2008).
A striking example of such selection was described in a 115year-old woman whose peripheral white blood cells were shown
to be primarily the offspring of only two related HSC clones,
whose cargo of approximately 450 somatic mutations did not
include known leukemogenic mutations (Holstege et al., 2014).
In the absence of somatic driver mutations, it is probable that
such selection is driven by well-demonstrated epigenetic differences between individual HSCs (Fraga et al., 2005) or by stochastic events. Furthermore, clonal hemopoiesis in the absence
of a known leukemia-driver mutation was also well documented
recently (Genovese et al., 2014), and whereas unknown or undetected drivers may be responsible for many cases of this phenomenon, it is also highly plausible that a stochastic process
of clonal selection or loss may operate in others. Our study provides evidence that spliceosome gene mutations offer a means
to exploit age-related changes in hemopoiesis to drive clonal hemopoiesis in advanced old age, an observation that blurs the
boundary between driver and passenger mutations. Such
a context dependency is not a surprising attribute for the effects
of spliceosome mutations, which have not, so far, been shown to
impart a primary proliferative advantage to normal hemopoietic
stem and progenitor cells (Matsunawa et al., 2014; Visconte
et al., 2012).
A nal important nding of our study was the almost complete
absence of canonical NPM1 mutations in our collection of more
than 4,000 people, despite the use of a highly sensitive assay for
their detection, designed specically for this study. Among more
than 10 million mapped reads covering this mutation hot spot,
we identied only four reads in a single sample reporting a canonical mutation (mutation A; TCTG duplication). Given their frequency in myeloid leukemia (Cancer Genome Atlas Research
Network, 2013) and the fact that they are not late mutations
(Kronke et al., 2013; Shlush et al., 2014), this observation frames
NPM1 mutations as gatekeepers of leukemogenesis, i.e., their
acquisition appears to be closely associated with the development of frank leukemia. In this light, the frequent co-occurrence
of DNMT3A and NPM1 mutations suggests that the former
behave as rafts that enable NPM1 mutant clones to be
founded and expanded, thus facilitating onward evolution toward acute myeloid leukemia.
We used a highly sensitive method to search for evidence of
clonal hemopoiesis driven by 15 recurrent leukemogenic mutations in more than 4,000 individuals. Our results demonstrate
that the incidence of clonal hemopoiesis is much higher than
suggested by exome-sequencing studies, that spliceosome
gene mutations drive clonal outgrowth primarily in the context
of an aging hemopoietic compartment, and that NPM1 mutations do not drive ARCH, indicating that their acquisition is
closely associated with frank leukemia.
EXPERIMENTAL PROCEDURES
Patient Samples
Samples were obtained with written informed consent and in accordance with
the Declaration of Helsinki and appropriate ethics committee approvals from
all participants (approval reference numbers 10/H0604/02, 07/MRE05/44,
and 05/Q0106/74). Maternal consent was obtained for the use of cord blood
samples. Samples were obtained from 3,067 blood donors aged 1770 years
(WTCCC; UK Blood Services 1 [UKBS1] and UKBS2 common controls),
1,152 unselected individuals aged 6098 years (UKHLS; https://www.
understandingsociety.ac.uk/), 32 patients that had undergone a hemopoietic
stem cell transplant (12 autologous and 20 allogeneic; Tables S3 and S4)
1 month to 14 years previously, and 18 cord blood samples. Age distribution
of the WTCCC and UKHLS cohorts/samples is shown in Figure S1. Hemoglobin concentrations were available for a total of 3,587 of the 4,067 samples from
which adequate sequencing data were obtained for analysis, including 102 of
105 samples with mutations. Full blood count results were available for 2,952
WTCCC samples. The average blood donation frequency for WTCCC donors
was 1.6 donations of one unit per year. Details of donations by individual participants were not available.
Targeted Sequencing
Genomic DNA was used to simultaneously amplify several gene loci using
multiplex PCR, in order to capture and analyze 15 mutational hot spots enriched for, but not exclusive to, targets of mutations thought to arise early in
leukemogenesis (Table 1). We used three multiplex primer combinations
(Plex1-3), guided by our ndings, to capture the targeted mutational hot spots
(Table S1). Primers were designed using the Hi-Plex PCR-MPS (massively
parallel sequencing) strategy (Nguyen-Dumont et al., 2013), except for JAK2
V617 and Plex2 primers, which were designed using MPRIMER (Shen
et al., 2010). These and additional primer sequences used in each Plex and details of PCR- and DNA-sequencing protocols are detailed in Supplemental
Experimental Procedures. Methodological validation experiments are shown
in Figure S2.
Bioinformatic Analysis
Sequencing data were aligned to the human reference genome (hg19) using
BWA. Subsequently, the SAMTOOLS pileup command was used to generate
pileup les from the generated bam les (version 0.1.8; http://samtools.
sourceforge.net; Li et al., 2009). A exible in-house Perl script generated by
our group, MIDAS (Conte et al., 2013), was modied in order to interrogate
only the hot spot nucleotide positions of interest (those with reported mutations in the COSMIC database; Forbes et al., 2015) on the pileup le, considering only those reads with a sequence quality higher than 25 and a mapping
quality higher than 15. For each sample, the numbers of reads reporting the
reference and variant alleles at each position were extracted. VAFs were
derived by dividing the number of reads reporting the most-frequent variant
nucleotide to the total. In order to detect NPM1 mutations with high sensitivity,
Cell Reports 10, 12391245, March 3, 2015 2015 The Authors 1243
SUPPLEMENTAL INFORMATION
Chambers, S.M., Shaw, C.A., Gatza, C., Fisk, C.J., Donehower, L.A., and
Goodell, M.A. (2007). Aging hematopoietic stem cells decline in function and
exhibit epigenetic dysregulation. PLoS Biol. 5, e201.
Conte, N., Varela, I., Grove, C., Manes, N., Yusa, K., Moreno, T., Segonds-Pichon, A., Bench, A., Gudgin, E., Herman, B., et al. (2013). Detailed molecular
characterisation of acute myeloid leukaemia with a normal karyotype using targeted DNA capture. Leukemia 27, 18201825.
AUTHOR CONTRIBUTIONS
G.S.V. conceived and designed the study. G.S.V. and T. McKerrell supervised
the study, analyzed data, and wrote the manuscript. N.P. and T. McKerrell
performed experimental procedures. I.V. and T. Moreno wrote scripts and performed bioinformatics analysis. H.P., T. McKerrell, and G.S.V. performed statistical analyses. E.Z., C.S.G., M.A.Q., and R.R. contributed to study strategy
and to technical and analytical aspects. U.S.S.G., E.Z., W.O, J.C., C.C., J.B.,
J.S., C.H., M.A.S., and D.R.F. contributed to sample acquisition and subject
recruitment.
ACKNOWLEDGMENTS
This project was funded by a Wellcome Trust Clinician Scientist Fellowship
(100678/Z/12/Z; to T. McKerrell) and by the Wellcome Trust Sanger Institute
(grant number WT098051). G.S.V. is funded by a Wellcome Trust Senior
Fellowship in Clinical Science (WT095663MA), and work in his laboratory is
also funded by Leukaemia Lymphoma Research and the Kay Kendal
Leukaemia Fund. I.V. is funded by Spanish Ministerio de Economa y Competitividad subprograma Ramon y Cajal. C.S.G. is funded by a Leukaemia Lymphoma Research Clinical Research Training Fellowship. We thank Servicio
Santander Supercomputacion for their support. We acknowledge use of
DNA from The UK Blood Services Collection of Common Controls (UKBS
collection), funded by the Wellcome Trust grant 076113/C/04/Z, by the Juvenile Diabetes Research Foundation grant WT061858, and by the National Institute of Health Research of England. The collection was established as part of
the Wellcome Trust Case-Control Consortium. We also gratefully acknowledge use of blood DNA samples and data from participants of the UK Household Longitudinal Study (https://www.understandingsociety.ac.uk/), collected
by NatCen and the Institute for Social and Economic Research, University of
Essex, and funded by the Economic and Social Research Council, UK. We
thank the Cambridge Blood and Stem Cell Biobank and the Cancer Molecular
Diagnosis Laboratory, Cambridge Biomedical Research Centre (National Institute for Health Research, UK) for help with sample collection and processing.
Finally, we thank Nathalie Smerdon, Richard Rance, Lucy Hildyard, Ben Softly,
and Britt Killian for help with sample management, DNA sequencing, and data
processing. G.S.V. is a consultant for KYMAB and receives an educational
grant from Celgene.
1244 Cell Reports 10, 12391245, March 3, 2015 2015 The Authors
Ding, L., Ley, T.J., Larson, D.E., Miller, C.A., Koboldt, D.C., Welch, J.S.,
Ritchey, J.K., Young, M.A., Lamprecht, T., McLellan, M.D., et al. (2012). Clonal
evolution in relapsed acute myeloid leukaemia revealed by whole-genome
sequencing. Nature 481, 506510.
Forbes, S.A., Beare, D., Gunasekaran, P., Leung, K., Bindal, N., Boutselakis,
H., Ding, M., Bamford, S., Cole, C., Ward, S., et al. (2015). COSMIC: exploring
the worlds knowledge of somatic mutations in human cancer. Nucleic Acids
Res. 43, D805D811.
Ford, A.M., Bennett, C.A., Price, C.M., Bruin, M.C., Van Wering, E.R., and
Greaves, M. (1998). Fetal origins of the TEL-AML1 fusion gene in identical
twins with leukemia. Proc. Natl. Acad. Sci. USA 95, 45844588.
Fraga, M.F., Ballestar, E., Paz, M.F., Ropero, S., Setien, F., Ballestar, M.L.,
Heine-Suner, D., Cigudosa, J.C., Urioste, M., Benitez, J., et al. (2005). Epigenetic differences arise during the lifetime of monozygotic twins. Proc. Natl.
Acad. Sci. USA 102, 1060410609.
Genovese, G., Kahler, A.K., Handsaker, R.E., Lindberg, J., Rose, S.A., Bakhoum, S.F., Chambert, K., Mick, E., Neale, B.M., Fromer, M., et al. (2014).
Clonal hematopoiesis and blood-cancer risk inferred from blood DNA
sequence. N. Engl. J. Med. 371, 24772487.
Greaves, M., and Maley, C.C. (2012). Clonal evolution in cancer. Nature 481,
306313.
Haferlach, T., Nagata, Y., Grossmann, V., Okuno, Y., Bacher, U., Nagae, G.,
Schnittger, S., Sanada, M., Kon, A., Alpermann, T., et al. (2014). Landscape
of genetic lesions in 944 patients with myelodysplastic syndromes. Leukemia
28, 241247.
Holstege, H., Pfeiffer, W., Sie, D., Hulsman, M., Nicholas, T.J., Lee, C.C., Ross,
T., Lin, J., Miller, M.A., Ylstra, B., et al. (2014). Somatic mutations found in the
healthy blood compartment of a 115-yr-old woman demonstrate oligoclonal
hematopoiesis. Genome Res. 24, 733742.
Jacobs, K.B., Yeager, M., Zhou, W., Wacholder, S., Wang, Z., Rodriguez-Santiago, B., Hutchinson, A., Deng, X., Liu, C., Horner, M.J., et al. (2012). Detectable clonal mosaicism and its relationship to aging and cancer. Nat. Genet. 44,
651658.
Jaiswal, S., Fontanillas, P., Flannick, J., Manning, A., Grauman, P.V., Mar,
B.G., Lindsley, R.C., Mermel, C.H., Burtt, N., Chavez, A., et al. (2014). Agerelated clonal hematopoiesis associated with adverse outcomes. N. Engl. J.
Med. 371, 24882498.
Kronke, J., Bullinger, L., Teleanu, V., Tschurtz, F., Gaidzik, V.I., Kuhn, M.W.,
Rucker, F.G., Holzmann, K., Paschka, P., Kapp-Schworer, S., et al. (2013).
Clonal evolution in relapsed NPM1-mutated acute myeloid leukemia. Blood
122, 100108.
Kyle, R.A., Therneau, T.M., Rajkumar, S.V., Offord, J.R., Larson, D.R., Plevak,
M.F., and Melton, L.J., 3rd. (2002). A long-term study of prognosis in monoclonal gammopathy of undetermined signicance. N. Engl. J. Med. 346,
564569.
Laurie, C.C., Laurie, C.A., Rice, K., Doheny, K.F., Zelnick, L.R., McHugh, C.P.,
Ling, H., Hetrick, K.N., Pugh, E.W., Amos, C., et al. (2012). Detectable clonal
mosaicism from birth to old age and its relationship to cancer. Nat. Genet.
44, 642650.
Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., Marth, G.,
Abecasis, G., and Durbin, R.; 1000 Genome Project Data Processing Subgroup (2009). The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 20782079.
Lin, C.C., Hou, H.A., Chou, W.C., Kuo, Y.Y., Wu, S.J., Liu, C.Y., Chen, C.Y.,
Tseng, M.H., Huang, C.F., Lee, F.Y., et al. (2014). SF3B1 mutations in patients
with myelodysplastic syndromes: the mutation is stable during disease evolution. Am. J. Hematol. 89, E109E115.
Matsunawa, M., Yamamoto, R., Sanada, M., Sato-Otsubo, A., Shiozawa, Y.,
Yoshida, K., Otsu, M., Shiraishi, Y., Miyano, S., Isono, K., et al. (2014). Haploinsufciency of Sf3b1 leads to compromised stem cell function but not to myelodysplasia. Leukemia 28, 18441850.
Nguyen-Dumont, T., Pope, B.J., Hammet, F., Southey, M.C., and Park, D.J.
(2013). A high-plex PCR approach for massively parallel sequencing. Biotechniques 55, 6974.
Nowell, P.C. (1976). The clonal evolution of tumor cell populations. Science
194, 2328.
Rossi, D.J., Jamieson, C.H., and Weissman, I.L. (2008). Stems cells and the
pathways to aging and cancer. Cell 132, 681696.
Rube, C.E., Fricke, A., Widmann, T.A., Furst, T., Madry, H., Pfreundschuh, M.,
and Rube, C. (2011). Accumulation of DNA damage in hematopoietic stem and
progenitor cells during human aging. PLoS ONE 6, e17487.
Shen, Z., Qu, W., Wang, W., Lu, Y., Wu, Y., Li, Z., Hang, X., Wang, X., Zhao, D.,
and Zhang, C. (2010). MPprimer: a program for reliable multiplex PCR primer
design. BMC Bioinformatics 11, 143.
Shlush, L.I., Zandi, S., Mitchell, A., Chen, W.C., Brandwein, J.M., Gupta, V.,
Kennedy, J.A., Schimmer, A.D., Schuh, A.C., Yee, K.W., et al.; HALT Pan-Leukemia Gene Panel Consortium (2014). Identication of pre-leukaemic haematopoietic stem cells in acute leukaemia. Nature 506, 328333.
Taraldsrud, E., Grgaard, H.K., Solheim, S., Lunde, K., Flisand, Y., Arnesen,
H., Seljeot, I., and Egeland, T. (2009). Age and stress related phenotypical
changes in bone marrow CD34+ cells. Scand. J. Clin. Lab. Invest. 69, 7984.
Vas, V., Senger, K., Dorr, K., Niebel, A., and Geiger, H. (2012). Aging of the
microenvironment inuences clonality in hematopoiesis. PLoS ONE 7, e42080.
Visconte, V., Rogers, H.J., Singh, J., Barnard, J., Bupathi, M., Traina, F.,
McMahon, J., Makishima, H., Szpurka, H., Jankowska, A., et al. (2012).
SF3B1 haploinsufciency leads to formation of ring sideroblasts in myelodysplastic syndromes. Blood 120, 31733186.
Woolthuis, C.M., de Haan, G., and Huls, G. (2011). Aging of hematopoietic
stem cells: Intrinsic changes or micro-environmental effects? Curr. Opin. Immunol. 23, 512517.
Wu, S.J., Kuo, Y.Y., Hou, H.A., Li, L.Y., Tseng, M.H., Huang, C.F., Lee, F.Y.,
Liu, M.C., Liu, C.W., Lin, C.T., et al. (2012). The clinical implication of SRSF2
mutation in patients with myelodysplastic syndrome and its stability during disease evolution. Blood 120, 31063111.
Papaemmanuil, E., Gerstung, M., Malcovati, L., Tauro, S., Gundem, G., Van
Loo, P., Yoon, C.J., Ellis, P., Wedge, D.C., Pellagatti, A., et al.; Chronic Myeloid
Disorders Working Group of the International Cancer Genome Consortium
(2013). Clinical and biological implications of driver mutations in myelodysplastic syndromes. Blood 122, 36163627, quiz 3699.
Xie, M., Lu, C., Wang, J., McLellan, M.D., Johnson, K.J., Wendl, M.C.,
McMichael, J.F., Schmidt, H.K., Yellapantula, V., Miller, C.A., et al. (2014).
Age-related mutations associated with clonal hematopoietic expansion and
malignancies. Nat. Med. 20, 14721478.
Rossi, D.J., Bryder, D., Zahn, J.M., Ahlenius, H., Sonu, R., Wagers, A.J., and
Weissman, I.L. (2005). Cell intrinsic alterations underlie hematopoietic stem
cell aging. Proc. Natl. Acad. Sci. USA 102, 91949199.
Zhang, J., Niu, C., Ye, L., Huang, H., He, X., Tong, W.G., Ross, J., Haug, J.,
Johnson, T., Feng, J.Q., et al. (2003). Identication of the haematopoietic
stem cell niche and control of the niche size. Nature 425, 836841.
Cell Reports 10, 12391245, March 3, 2015 2015 The Authors 1245
Cell Reports
Report
Accelerating Novel Candidate Gene Discovery in
Neurogenetic Disorders via Whole-Exome Sequencing
of Prescreened Multiplex Consanguineous Families
Anas M. Alazami,1,23 Nisha Patel,1,23 Hanan E. Shamseldin,1,23 Shamsa Anazi,1 Mohammed S. Al-Dosari,2
Fatema Alzahrani,1 Hadia Hijazi,1 Muneera Alshammari,3 Mohammed A. Aldahmesh,1 Mustafa A. Salih,3 Eissa Faqeih,4
Amal Alhashem,5,6 Fahad A. Bashiri,3 Mohammed Al-Owain,5,7 Amal Y. Kentab,3 Sameera Sogaty,8 Saeed Al Tala,9
Mohamad-Hani Temsah,3 Maha Tulbah,10 Rasha F. Aljelaify,11 Saad A. Alshahwan,6 Mohammed Zain Seidahmed,12
Adnan A. Alhadid,3 Hesham Aldhalaan,13 Fatema AlQallaf,13 Wesam Kurdi,10 Majid Alfadhel,14 Zainab Babay,15
Mohammad Alsogheer,16 Namik Kaya,1 Zuhair N. Al-Hassnan,5,7 Ghada M.H. Abdel-Salam,17 Nouriya Al-Sannaa,18
Fuad Al Mutairi,14 Heba Y. El Khashab,3,19 Saeed Bohlega,13 Xiaofei Jia,20 Henry C. Nguyen,20 Rakad Hammami,1
Nouran Adly,1 Jawahir Y. Mohamed,1 Firdous Abdulwahab,1 Niema Ibrahim,1 Ewa A. Naim,1,21 Banan Al-Younes,1,21
Brian F. Meyer,1,21 Mais Hashem,1 Ranad Shaheen,1 Yong Xiong,20 Mohamed Abouelhoda,1,21
Abdulrahman A. Aldeeri,1,22 Dorota M. Monies,1,21 and Fowzan S. Alkuraya1,5,21,*
1Department
of Genetics, King Faisal Specialist Hospital and Research Center, Riyadh 11211, Saudi Arabia
of Pharmacognosy, College of Pharmacy, King Saud University, Riyadh 11451, Saudi Arabia
3Department of Pediatrics, King Khalid University Hospital and College of Medicine, King Saud University, Riyadh 11451, Saudi Arabia
4Department of Pediatrics, King Fahad Medical City, Riyadh 11525, Saudi Arabia
5Department of Anatomy and Cell Biology, College of Medicine, Alfaisal University, Riyadh 11533, Saudi Arabia
6Department of Pediatrics, Prince Sultan Military Medical City, Riyadh 11159, Saudi Arabia
7Department of Medical Genetics, King Faisal Specialist Hospital and Research Center, Riyadh 11211, Saudi Arabia
8Department of Pediatrics, King Fahad General Hospital, Jeddah 23325, Saudi Arabia
9Department of Pediatrics, Armed Forces Hospital, Khamis Mushayt 62413, Saudi Arabia
10Department of Obstetrics & Gynecology, King Faisal Specialist Hospital, Riyadh 11211, Saudi Arabia
11Center of Excellence for Genomics, King Abdulaziz City for Science and Technology, Riyadh 11442, Saudi Arabia
12Department of Pediatrics, Security Forces Hospital, Riyadh 12625, Saudi Arabia
13Department of Neurosciences, King Faisal Specialist Hospital and Research Center, Riyadh 11211, Saudi Arabia
14Division of Genetics, Department of Pediatrics, King Saud bin Abdulaziz University for Health Sciences, King Abdulaziz Medical City, Riyadh
14611, Saudi Arabia
15Department of Obstetrics and Gynecology, College of Medicine, King Saud University, Riyadh 11451, Saudi Arabia
16Department of Psychiatry, College of Medicine, King Saud University, Riyadh 11451, Saudi Arabia
17Department of Clinical Genetics, Human Genetics and Genome Research Division, National Research Centre, Cairo 12345, Egypt
18Department of Pediatrics, Johns Hopkins Aramco Healthcare, Dhahran 34465, Saudi Arabia
19Department of Pediatrics, Childrens Hospital, Ain Shams University, Cairo 01234, Egypt
20Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT 06520, USA
21Saudi Human Genome Program, King Abdulaziz City for Science and Technology, Riyadh 11442, Saudi Arabia
22Department of Internal Medicine, College of Medicine, King Saud University, Riyadh 11451, Saudi Arabia
23Co-rst author
*Correspondence: [email protected]
http://dx.doi.org/10.1016/j.celrep.2014.12.015
This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/3.0/).
2Department
SUMMARY
Our knowledge of disease genes in neurological disorders is incomplete. With the aim of closing this gap,
we performed whole-exome sequencing on 143
multiplex consanguineous families in whom known
disease genes had been excluded by autozygosity
mapping and candidate gene analysis. This prescreening step led to the identication of 69 recessive genes not previously associated with disease,
of which 33 are here described (SPDL1, TUBA3E,
INO80, NID1, TSEN15, DMBX1, CLHC1, C12orf4,
WDR93, ST7, MATN4, SEC24D, PCDHB4, PTPN23,
TAF6, TBCK, FAM177A1, KIAA1109, MTSS1L,
148 Cell Reports 10, 148161, January 13, 2015 2015 The Authors
clinical presentations that range from the common e.g., intellectual disability (>1%) to the very rare, e.g., neurodegeneration
with brain iron accumulation (one to three per 106) (Kalman
et al., 2012; Maulik et al., 2011). The highly prevalent involvement
of the nervous system in many Mendelian disorders coincides
with the observation that >80% of all human genes are expressed at some stage of brain development (Hawrylycz et al.,
2012) and suggests that the brain is one of the most vulnerable
organs to genetic perturbation. In fact high-resolution microarray
analysis of the human genome reveals that intellectual disability
is the common phenotypic denominator of genomic disorders
that involve losses or gains of genes (Coe et al., 2012).
Variances in clinical presentation are a major obstacle in
establishing a working molecular classication of neurological
disease, because even where the clinical presentation is highly
specic, genetic heterogeneity is the rule. In the setting of
autosomal recessive neurogenetic disorders where parents are
related, a homozygosity scan can serve as a guide to the underlying genetic cause even when the phenotype is atypical (Alkuraya, 2010). Another major challenge in assigning a molecular
classication is that many neurological disease genes have not
been identied yet.
Novel disease gene discovery in this eld has been tremendously abetted by next-generation sequencing, a tool with the
capacity to, theoretically, unravel the genetic cause of all neurological diseases. This full theoretical potential has not yet been
reached unfortunately, although the technology continues to
evolve. For example, two large studies on the genetics of intellectual disability using whole-exome sequencing (WES) provided
a yield of 16%55%, and even though the collective sample size
was >150, only seven novel genes were identied (de Ligt et al.,
2012; Rauch et al., 2012). In these studies, samples could not be
enriched for novel gene discovery, and simultaneously the anticipated mutations were heterozygous, detection of which poses a
challenge for the currently available sequencing technology
(especially regarding insertions/deletions) (Harismendy et al.,
2009). The presence of these two obstacles likely hindered the
authors ability to obtain a higher yield. Thus, alternative/complementary approaches are required to facilitate the discovery of
novel neurogenetic disease genes. In this study, we show that
the analysis of the entire set of autozygous intervals per individual (the autozygome) in multiplex consanguineous patients, as a
prescreen, can markedly increase the yield of WES to identify
candidate genes not previously associated with disease. Even
when known genes were identied using this approach, the
phenotype was often sufciently different to explain why the
gene had been missed by the autozygosity lter. The 33 candidate disease genes we showcase in this study will augment
the global hunt for the genetics of brain development and
function.
RESULTS
Clinical Report
In total, 143 multiplex families met our inclusion criteria (a neurogenetic diagnosis, positive family history, consanguineous
parents, and no candidates identied by autozygosity mapping).
Intellectual disability was the most common clinical feature.
150 Cell Reports 10, 148161, January 13, 2015 2015 The Authors
SPG20
CYP27A1
NPC2
ARFGEF2 NM_006420:c.656_657insC:p.P219fs
PNKP
RYR1
FBN2
12DG0096
12DG1571
10DG0672
11DG1951
11DG1510
12DG0975
12DG0104
EVC2
DDHD2
WDR81
ZNF526
10DG1721
13DG0583
13DG0010
12DG0685
SBF1
AP4M1
09DG0301
11DG1767
NM_004722:c.C952T:p.R318*
NM_020247:c.1744 dup:p.S582Kfs*148
09DG00930 ADCK3
NM_000059.3:c.9152 delC:p.P3051Hfs*11
08DG00385 BRCA2
NM_007254:c.1250_1251insAACGGGTCGCCATCGAC:p.R418Tfs*55
NM_001142294:c.1450_1451insA:p.T484fs
GM2A
09DG0057
Mutation
Gene
ID
Charcot-Marie-Tooth disease
type 4B3
nonsyndromic intellectual
disability
spastic paraplegia
Ellis-van-Creveld syndrome
Fanconi anemia
congenital contractural
arachnodactyly
minicore myopathy
progressive microcephaly,
infantile-onset seizures,
and developmental delay
periventricular heterotopia
with microcephaly
Nieman-Pick disease
cerebrotendinous
xanthomatosis
spastic paraplegia
GM2-gangliosidosis,
AB variant
Published Phenotype
this study
this study
this study
this study
Shaheen
et al. (2013)
Shaheen
et al. (2014)
this study
this study
this study
this study
this study
this study
this study
this study
Reference
this study
Meckel-Gruber syndrome
primordial dwarsm
Observed Phenotype
Table 1. Cases with Mutations in Known Disease Genes following WES, Where the Patient Phenotype Diverged from the Established Literature
152 Cell Reports 10, 148161, January 13, 2015 2015 The Authors
Phenotype
primary microcephaly
and neonatal death
microlissencephaly
and global
developmental delay
primary microcephaly
and global
developmental delay
hydrocephalus,
muscle weakness
and global
developmental delay
primary microcephaly
and global
developmental delay
ID
12DG1528
11DG0443
10DG1705
08DG00041
13DG0167
TSEN15
NID1
INO80
TUBA3E
CCDC99
(SPDL1)
Gene
NM_001127394:
c.226T > G:p.W76G
NM_002508.2:
c.3385+1G > A
INO80:NM_017553:
c.1501T > C:p.S501P,
INO80:NM_017553:
c.3737G > A:p.R1246Q
NM_207312:c.643C >
T:p.R215C
NM_017785:c.1724_
1747 del:p.S575_T582
del
Mutation
missense
splicing
(in-frame
insertion
conrmed
by RTPCR)
missense
missense
in-frame
deletion
Mutation
Type
Reference
PMID:
18711368
PMID:
12480912
and 23674478
PMID:
21947284,
19829069,
24029917,
and 20118933
PMID:
24860126
PMID:
20427577
and 24875059
tRNA splicing
Mutations in other family members,
endonuclease 15 e.g., TSEN54, have been linked to
pontocerebellar hypoplasia. Only
surviving variant and segregates
in family.
Nidogen 1
INO80 complex
subunit E
tubulin, alpha 3e
spindle
apparatus
coiled-coil
protein 1
Table 2. Cases with Mutations in Candidate Genes, along with Available Evidence from the Literature to Support Candidacy
Cell Reports 10, 148161, January 13, 2015 2015 The Authors 153
Phenotype
global developmental
delay, epilepsy and
poor weight gain
myopathy
global developmental
delay
autism spectrum
disorder
global developmental
delay and brain
atrophy
holoprosencephaly
intellectual disability
and epilepsy
primary microcephaly
and global
developmental delay
12DG0929
09DG0405
09DG00102
11DG1513
11DG2479
12DG1901
12DG2051
10DG1069
08DG-00322
Continued
ID
Table 2.
PTPN23
PCDHB4
SEC24D
MATN4
ST7
WDR93
C12orf4
C2orf63
(CLHC1)
DMBX1
Gene
NM_015466:c.3995G >
T:p.R1332L
NM_018938.2:c.915
del:p.K305Nfs*12
NM_014822:c.697G >
C:p.G233R
NM_030590:c.515G >
C:p.G172A
NM_021908:c.489T >
G:p.Y163X
NM_020212:c.280T >
C:p.Y94H
NM_020374:exon6:
c.637_638insAAAC:
p.K213fs
NM_001135598:
c.779G > A:p.R260Q
NM_147192:c.367C >
T:p.R123W
Mutation
missense
frameshift
missense
missense
stopgain
missense
frameshift
missense
missense
Mutation
Type
protein tyrosine
phosphatase,
nonreceptor
type 23
protocadherin
beta 4
SEC24 family,
member D
(S. cerevisiae)
Matrilin 4
suppression of
tumorigenicity 7
WD repeat
domain 93
chromosome 12
open reading
frame 4
Reference
PMID:
19378249
PMID:
22495309
and 22765916
PMID:
23596517
PMID:
11549321
PMID:
10889047
PMID:
23664116
PMID:
15314164
and 17873059
154 Cell Reports 10, 148161, January 13, 2015 2015 The Authors
Phenotype
global developmental
delay and
dysmorphism
global developmental
delay, epilepsy,
dysmorphism,
hypotonia, and VSD
macrocephaly, ID,
dolichocephaly,
and mild obesity
Dandy-Walker
malformation,
hydrocephalus,
exed deformity,
club feet,
micrognathia,
and pleural effusion
neurodegeneration
and brain iron
accumulation
11DG0932
10DG1670
13DG1472
13DG1900
10DG0264
Continued
ID
Table 2.
TAF6 RNA
polymerase II,
TATA box binding
protein (TBP)associated factor
NM_138383:c.1790C >
T:p.T597M
MTSS1L
missense
KIAA1109
frameshift
metastasis
suppressor
1-like
KIAA1109
family with
sequence
similarity 177,
member A1
PMID:
17956977
PMID:
19640479
PMID:
23977024,
19963289
Reference
splicing
TBC1 domain
TBCK knockdown signicantly
(frameshift) containing kinase suppresses mTOR signaling,
which plays a critical role in
multiple neurological disorders.
Only surviving variant and
segregates in family.
missense
Mutation
Type
NM_001079519:c.297_
298insA:p.L99fs
NM_033115:
c.1708+1G > A
NM_005641.3:c.212T >
C:p.I71T
Mutation
FAM177A1
TBCK
TAF6
Gene
Cell Reports 10, 148161, January 13, 2015 2015 The Authors 155
severe psychomotor
retardation, seizure,
and cerebellar
hypoplasia
global developmental
delay and ADHD
neurodegenerative
disease
13DG2274
13DG0832
08DGRC00077
neurodegeneration
with marked white
matter changes with
high lactate peak in
the brain, consistent
with mitochondrial
encephalopathy
primary microcephaly
13DG1935
14DG0152
Phenotype
Continued
ID
Table 2.
ISCA2
ARV1
NM_194279:c.229G >
A:p.G77S
NM_022786:c.565G >
A:p.G189R
NM_005441:c.496A >
G:p.I166V
NM_016121:c.1036_
1073 del:p.P346Tfs*4
KCTD3
CHAF1B
NM_194293:c.4495G >
A:p.E1499K
Mutation
XIRP1
Gene
missense
missense
missense
frameshift
missense
Mutation
Type
iron-sulfur
cluster
assembly 2
homolog
ARV1 homolog
(S. cerevisiae)
chromatin
assembly
factor 1,
subunit B
potassium
channel
tetramerization
domain
containing 3
Reference
PMID:
12145310
PMID:
22177093
PMID:
19623214
PMID:
22366181
PMID:
Maps to the only shared haplotype
22323289
in the extended family. ISCA2 is a
member of mitochondrial iron-sulfur
cluster (ISC) assembly machinery.
Depletion of ISCA2 results in
massively swollen mitochondria that
are devoid of cristae membranes
indicating that it is required for normal
mitochondrial biogenesis. Same
mutation was identied in other
families with identical presentation
(Z.N.A.-H., M. Al-Dosary, M. Alfadhel,
E.F., M. Alsagob, R. Kenana,
R. Almas, O.S. Al-Harazi, H. Al-Hindi,
O.I. Malibari, F.B. Almutari,
T. Al-Sheddi, S. Tulbah, F. Alhadeq,
R. Alamro, A, AlAsmari, M. Almuntashri,
H. Alshaalan, F.A. Al-Mohanna, D. Colak,
N.K., unpublished data). Only surviving
variant and segregates in family.
156 Cell Reports 10, 148161, January 13, 2015 2015 The Authors
cerebellar vermis
hypoplasia, DandyWalker malformation,
hydrocephalus,
developmental delay
10DG0934
11DG0417
14DG0221
cerebellar atrophy,
hydrocephalus, and
global developmental
(cognitive, speech,
and motor) delay
global developmental
delay, typical
Joubert syndrome,
MRI ndings
13DG0274
DPH1
PDPR
NM_001168215:
c.95+3A > G
NM_020401:c.303G >
A:p.M101I
NM_001383.3:c.701T >
C:p.L234P
NM_017990:c.1360G >
T:p.G454C
NM_001146312:
c.1252A > G:p.I418V
MYOCD
intellectual disability
and epilepsy
NM_016077.3:c.254A >
C:p.Q85P
13DG1549
PTRH2
Mutation
NM_015721:c.2452T >
C:p.W818R
global developmental
delay, hearing loss,
and ataxia
13DG0215
Gene
GEMIN4
08DG0048510DG070313DG1542 (08DG00485) global
developmental delay,
severe dystonia, and
congenital cataract
(10DG0703) global
developmental
delay and congenital
cataract(13DG1542)
global developmental
delay, congenital
cataract, tubulopathy,
and severe osteopenia
Phenotype
Continued
ID
Table 2.
splice site
splice site
missense
missense
missense
missense
missense
Mutation
Type
transmembrane
protein 92
nucleoporin
107 kDa
diphthamide
biosynthesis 1
pyruvate
dehydrogenase
phosphatase
regulatory
subunit
myocardin
gem (nuclear
organelle)
associated
protein 4
peptidyl-tRNA
hydrolase 2
PMID:
18765564
PMID:
12867591
PMID:
11914277
PMID:
18218778
Reference
Diphthamide is a unique
posttranslationally modied
histidine found only in translation
elongation factor-2 (eEF2), which
is linked to spinocerebellar ataxia.
Diphthamide modication of eEF2
is essential for normal mouse
development. Only surviving variant
and segregates in family.
missense
family with
sequence
similarity 120A
opposite strand
Reference
Only surviving variant and segregates
in family.
erythrocyte
membrane
protein band
4.1 like 4A
missense
NM_022140:c.1298C >
T:p.S433L
EPB41L4A
spastic paraplegia,
failure to thrive.
10DG0840
Mutation
Type
Gene
Mutation
Phenotype
ID
Continued
Table 2.
Figure 1. Pedigrees of Families with Novel Candidate Genes following WES Analysis
The family ID is presented above each pedigree. A red box indicates the affected family member who was submitted for WES, and asterisks denote all the family
members we had access to for conrming segregation.
Recruitment
Epilepsy
Ataxia
Ausm
Intellectual
Disability
Primary
Microcephaly/
Brain
malformaons
Enrichment step
Global
developmental
delay
WES
Autozygome-guided
sequencing of known
disease genes
is negative
n=143 families
Solved cases
(105)
Unsolved: no
remaining
variants (27)
Exome Sequencing
Exome capture was performed using the TruSeq Exome Enrichment kit
(Illumina). See Supplemental Information for more details. A summary of the
quality control data for exome sequencing is provided in Table S2.
Supplemental Information includes two gures and two tables and can be
found with this article online at http://dx.doi.org/10.1016/j.celrep.2014.12.015.
AUTHOR CONTRIBUTIONS
A.M.A. collected and analyzed data and wrote the manuscript, N.P. collected
and analyzed data and wrote the manuscript, H.E.S. collected and analyzed
data and wrote the manuscript.
ACCESSION NUMBERS
ACKNOWLEDGMENTS
All variants within the 143 exomes in this study can be accessed through the
following link (part of the Saudi Variome Database): http://shgp.kfshrc.edu.
We thank all the families for their enthusiastic participation. This work was
funded in part by KACST 13-BIO1113-20 (F.S.A.). We acknowledge the Saudi
Cell Reports 10, 148161, January 13, 2015 2015 The Authors 159
tion and the rst prenatal screening of cathepsin D deciency (CLN10). Acta
Neuropathol. 117, 201208.
Frosk, P., Weiler, T., Nylen, E., Sudha, T., Greenberg, C.R., Morgan, K.,
Fujiwara, T.M., and Wrogemann, K. (2002). Limb-girdle muscular dystrophy
type 2H associated with mutation in TRIM32, a putative E3-ubiquitin-ligase
gene. Am. J. Hum. Genet. 70, 663672.
Greten-Harrison, B., Polydoro, M., Morimoto-Tomita, M., Diao, L., Williams,
A.M., Nie, E.H., Makani, S., Tian, N., Castillo, P.E., Buchman, V.L., and Chandra, S.S. (2010). abg-Synuclein triple knockout mice reveal age-dependent
neuronal dysfunction. Proc. Natl. Acad. Sci. USA 107, 1957319578.
Harismendy, O., Ng, P.C., Strausberg, R.L., Wang, X., Stockwell, T.B., Beeson, K.Y., Schork, N.J., Murray, S.S., Topol, E.J., Levy, S., and Frazer, K.A.
(2009). Evaluation of next generation sequencing platforms for population
targeted sequencing studies. Genome Biol. 10, R32.
Hawrylycz, M.J., Lein, E.S., Guillozet-Bongaarts, A.L., Shen, E.H., Ng, L.,
Miller, J.A., van de Lagemaat, L.N., Smith, K.A., Ebbert, A., Riley, Z.L., et al.
(2012). An anatomically comprehensive atlas of the adult human brain transcriptome. Nature 489, 391399.
Jaglin, X.H., Poirier, K., Saillour, Y., Buhler, E., Tian, G., Bahi-Buisson, N.,
Fallet-Bianco, C., Phan-Dinh-Tuy, F., Kong, X.P., Bomont, P., et al. (2009).
Mutations in the b-tubulin gene TUBB2B result in asymmetrical polymicrogyria. Nat. Genet. 41, 746752.
Kalman, B., Lautenschlaeger, R., Kohlmayer, F., Buchner, B., Kmiec, T., Klopstock, T., and Kuhn, K.A. (2012). An international registry for neurodegeneration with brain iron accumulation. Orphanet J. Rare Dis. 7, 66.
Keays, D.A., Tian, G., Poirier, K., Huang, G.-J., Siebold, C., Cleak, J., Oliver,
P.L., Fray, M., Harvey, R.J., Molnar, Z., et al. (2007). Mutations in a-tubulin
cause abnormal neuronal migration in mice and lissencephaly in humans.
Cell 128, 4557.
Kelley, L.A., and Sternberg, M.J. (2009). Protein structure prediction on the
Web: a case study using the Phyre server. Nat. Protoc. 4, 363371.
Carr, I.M., Bhaskar, S., OSullivan, J., Aldahmesh, M.A., Shamseldin, H.E.,
Markham, A.F., Bonthron, D.T., Black, G., and Alkuraya, F.S. (2013).
Autozygosity mapping with exome sequence data. Hum. Mutat. 34, 5056.
Mallaret, M., Synofzik, M., Lee, J., Sagum, C.A., Mahajnah, M., Sharkia, R.,
Drouot, N., Renaud, M., Klein, F.A., and Anheim, M. (2014). The tumour suppressor gene WWOX is mutated in autosomal recessive cerebellar ataxia
with epilepsy and mental retardation. Brain 137, 411419.
Cassandrini, D., Biancheri, R., Tessa, A., Di Rocco, M., Di Capua, M., Bruno,
C., Denora, P.S., Sartori, S., Rossi, A., Nozza, P., et al. (2010). Pontocerebellar
hypoplasia: clinical, pathologic, and genetic studies. Neurology 75, 1459
1464.
Martin, H.C., Kim, G.E., Pagnamenta, A.T., Murakami, Y., Carvill, G.L., Meyer, E.,
Copley, R.R., Rimmer, A., Barcia, G., Fleming, M.R., et al.; WGS500 Consortium
(2014). Clinical whole-genome sequencing in severe early-onset epilepsy reveals
new genes and improves molecular diagnosis. Hum. Mol. Genet. 23, 32003211.
Chen, J.-F., Zhang, Y., Wilde, J., Hansen, K.C., Lai, F., and Niswander, L.
(2014). Microcephaly disease gene Wdr62 regulates mitotic progression of
embryonic neural stem cells and brain size. Nat. Commun. 5, 3885.
Maulik, P.K., Mascarenhas, M.N., Mathers, C.D., Dua, T., and Saxena, S.
(2011). Prevalence of intellectual disability: a meta-analysis of populationbased studies. Res. Dev. Disabil. 32, 419436.
Coe, B.P., Girirajan, S., and Eichler, E.E. (2012). A genetic model for neurodevelopmental disease. Curr. Opin. Neurobiol. 22, 829836.
Miller, D.T., Adam, M.P., Aradhya, S., Biesecker, L.G., Brothman, A.R., Carter,
N.P., Church, D.M., Crolla, J.A., Eichler, E.E., Epstein, C.J., et al. (2010).
Consensus statement: chromosomal microarray is a rst-tier clinical diagnostic test for individuals with developmental disabilities or congenital anomalies. Am. J. Hum. Genet. 86, 749764.
Cormier-Daire, V., Amiel, J., Vuillaumier-Barrot, S., Tan, J., Durand, G.,
Munnich, A., Le Merrer, M., and Seta, N. (2000). Congenital disorders of glycosylation IIa cause growth retardation, mental retardation, and facial dysmorphism. J. Med. Genet. 37, 875877.
de Brouwer, A.P., Yntema, H.G., Kleefstra, T., Lugtenberg, D., Oudakker, A.R.,
de Vries, B.B., van Bokhoven, H., Van Esch, H., Frints, S.G., Froyen, G., et al.
(2007). Mutation frequencies of X-linked mental retardation genes in families
from the EuroMRX consortium. Hum. Mutat. 28, 207208.
de Ligt, J., Willemsen, M.H., van Bon, B.W., Kleefstra, T., Yntema, H.G., Kroes,
T., Vulto-van Silfhout, A.T., Koolen, D.A., de Vries, P., Gilissen, C., et al. (2012).
Diagnostic exome sequencing in persons with severe intellectual disability.
N. Engl. J. Med. 367, 19211929.
Ercan-Sencicek, A.G., Jambi, S., Franjic, D., Nishimura, S., Li, M., El-Fishawy,
P., Morgan, T.M., Sanders, S.J., Bilguvar, K., Suri, M., et al. (2014). Homozygous loss of DIAPH1 is a novel cause of microcephaly in humans. Eur. J.
Hum. Genet.
Fritchie, K., Siintola, E., Armao, D., Lehesjoki, A.-E., Marino, T., Powell, C.,
Tennison, M., Booker, J.M., Koch, S., Partanen, S., et al. (2009). Novel muta-
160 Cell Reports 10, 148161, January 13, 2015 2015 The Authors
Morrow, E.M., Yoo, S.-Y., Flavell, S.W., Kim, T.-K., Lin, Y., Hill, R.S., Mukaddes,
N.M., Balkhy, S., Gascon, G., Hashmi, A., et al. (2008). Identifying autism loci
and genes by tracing recent shared ancestry. Science 321, 218223.
Najmabadi, H., Hu, H., Garshasbi, M., Zemojtel, T., Abedini, S.S., Chen, W.,
Hosseini, M., Behjati, F., Haas, S., Jamali, P., et al. (2011). Deep sequencing
reveals 50 novel genes for recessive cognitive disorders. Nature 478, 5763.
Namavar, Y., Barth, P.G., Kasher, P.R., van Ruissen, F., Brockmann, K.,
Bernert, G., Writzl, K., Ventura, K., Cheng, E.Y., Ferriero, D.M., et al.; PCH
Consortium (2011). Clinical, neuroradiological and genetic ndings in pontocerebellar hypoplasia. Brain 134, 143156.
Poirier, K., Saillour, Y., Bahi-Buisson, N., Jaglin, X.H., Fallet-Bianco, C., Nabbout, R., Castelnau-Ptakhine, L., Roubertie, A., Attie-Bitach, T., Desguerre, I.,
et al. (2010). Mutations in the neuronal -tubulin subunit TUBB3 result in
malformation of cortical development and neuronal migration defects. Hum.
Mol. Genet. 19, 44624473.
Quaderi, N.A., Schweiger, S., Gaudenz, K., Franco, B., Rugarli, E.I., Berger,
W., Feldman, G.J., Volta, M., Andol, G., Gilgenkrantz, S., et al. (1997). Opitz
G/BBB syndrome, a defect of midline development, is due to mutations in a
new RING nger gene on Xp22. Nat. Genet. 17, 285291.
Rauch, A., Wieczorek, D., Graf, E., Wieland, T., Endele, S., Schwarzmayr, T.,
Albrecht, B., Bartholdi, D., Beygo, J., Di Donato, N., et al. (2012). Range of
genetic mutations associated with severe non-syndromic sporadic intellectual
disability: an exome sequencing study. Lancet 380, 16741682.
Shaheen, R., Faqeih, E., Alshammari, M.J., Swaid, A., Al-Gazali, L., Mardawi,
E., Ansari, S., Sogaty, S., Seidahmed, M.Z., AlMotairi, M.I., et al. (2013).
Genomic analysis of Meckel-Gruber syndrome in Arabs reveals marked
genetic heterogeneity and novel candidate genes. Eur. J. Hum. Genet. 21,
762768.
Shaheen, R., Faqeih, E., Ansari, S., Abdel-Salam, G., Al-Hassnan, Z.N., AlShidi, T., Alomar, R., Sogaty, S., and Alkuraya, F.S. (2014). Genomic analysis
of primordial dwarsm reveals novel disease genes. Genome Res. 24, 291299.
Sheftel, A.D., Wilbrecht, C., Stehling, O., Niggemeyer, B., Elsasser, H.-P.,
Muhlenhoff, U., and Lill, R. (2012). The human mitochondrial ISCA1, ISCA2,
and IBA57 proteins are required for [4Fe-4S] protein maturation. Mol. Biol.
Cell 23, 11571166.
Song, J., and Markley, J.L. (2007). Three-dimensional structure determined for
a subunit of human tRNA splicing endonuclease (Sen15) reveals a novel
dimeric fold. J. Mol. Biol. 366, 155164.
Tan, D., Rice, W.J., and Sosa, H. (2008). Structure of the kinesin13-microtubule ring complex. Structure 16, 17321739.
Tarpey, P.S., Smith, R., Pleasance, E., Whibley, A., Edkins, S., Hardy, C.,
OMeara, S., Latimer, C., Dicks, E., Menzies, A., et al. (2009). A systematic,
large-scale resequencing screen of X-chromosome coding exons in mental
retardation. Nat. Genet. 41, 535543.
Thevenon, J., Milh, M., Feillet, F., St-Onge, J., Duffourd, Y., Juge, C., Roubertie, A., Heron, D., Mignot, C., Raffo, E., et al. (2014). Mutations in SLC13A5
cause autosomal-recessive epileptic encephalopathy with seizure onset in
the rst days of life. Am. J. Hum. Genet. 95, 113120.
Cell Reports 10, 148161, January 13, 2015 2015 The Authors 161
Cell Reports
Report
Analysis of Intron Sequences Reveals Hallmarks
of Circular RNA Biogenesis in Animals
Andranik Ivanov,1 Sebastian Memczak,1,5 Emanuel Wyler,2,5 Francesca Torti,1,5 Hagit T. Porath,3 Marta R. Orejuela,1
Michael Piechotta,4 Erez Y. Levanon,3 Markus Landthaler,2 Christoph Dieterich,4 and Nikolaus Rajewsky1,*
1Laboratory
for Systems Biology of Gene Regulatory Elements, Berlin Institute for Medical Systems Biology, Max Delbruck Center for
Molecular Medicine, Robert Roessle Strae 10, 13125 Berlin-Buch, Germany
2Laboratory for RNA Biology and Posttranscriptional Regulation, Berlin Institute for Medical Systems Biology, Max Delbruck Center for
Molecular Medicine, Robert Roessle Strae 10, 13125 Berlin-Buch, Germany
3The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat-Gan 5290002, Israel
4Max Planck Institute for Biology of Ageing, Cologne, Joseph Stelzmann Strae 9B, 50931 Ko
ln, Germany
5These authors contributed equally to this work
*Correspondence: [email protected]
http://dx.doi.org/10.1016/j.celrep.2014.12.019
This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/3.0/).
SUMMARY
The resulting symmetric matrix (Figure 2A; Supplemental Experimental Procedures) shows that the intron pairs anking circularized exons stood out from all intron pairs. This was not the case
for the intron length-matched control transcripts (Figure 2A).
We then ranked, genome wide, all intron pairs of all transcripts
annotated in C. elegans by H and calculated the cumulative
fraction of known circRNAs as a function of intron-pair rank (Figure 2B). The top 4,200 (top 1%) of C. elegans intron pairs precisely bracketed 430 (38%) of the already annotated circRNAs,
a highly statistically signicant enrichment (hypergeometric test,
p value < 2.2 3 1016). Together, these results suggest that
RCMs provide highly signicantly improved accuracy in predicting circRNAs compared with intron length. However, we note
that (for C. elegans) the top BLAST score of RCMs in a pair of
introns yielded similar results.
To determine whether the presence of RCMs is a conserved
feature of circRNA formation, we next asked whether we could
predict circRNAs in human by intron sequence analysis.
The Presence of RCMs in Introns Flanking circRNAs Is a
Conserved Feature of circRNA Biogenesis
For human circRNAs, we used 1,067 exons overlapping circRNAs (Memczak et al., 2013). As in C. elegans, RCMs were highly
signicantly enriched (Fishers exact test, p value < 4 3 106;
172 Cell Reports 10, 170177, January 13, 2015 2015 The Authors
double-stranded ALU repeats (Athanasiadis et al., 2004; Levanon et al., 2004; Nishikura, 2010; Ramaswami et al., 2012). We
compared A-to-I conversions (Ramaswami and Li, 2014) in
1,500 bp regions anking circRNA splice sites with (1) other
splice sites in transcripts that produce circRNAs, and (2)
length-matched introns (Figure 4A). A-to-I conversions nearby
circRNA splice sites were enriched compared with controls. In
general, ALU repeats in circRNAs anking introns were edited
signicantly higher compared with expression- and lengthmatched controls (Figures S2A and S2B). Since A-to-I editing
is a hallmark of basepaired RNA, we asked whether RCMs are
preferentially located at sites of editing. We compared the position of the RCMs (dened as the nearest RCMs that match between the pair of introns bracketing circRNA) with the same controls as before (Figure 4B). These results suggest that indeed Ato-I editing preferentially occurs at regions that are basepaired
and proximal (upstream and downstream 200600 nt) to the
splice sites of circularized exons.
To test whether ADAR proteins are involved in circRNA
biogenesis, we codepleted ADAR1 and ADAR2 in HEK293 cells
using RNAi (Supplemental Experimental Procedures). We used
two controls: untreated total RNA and codepletion of three
proteins of the APOBEC family (APOBEC3B, APOBEC3C, and
APOBEC3F), which are known to bind mRNAs (Baltz et al.,
2012). APOBEC enzymes are known to edit single-stranded
DNA or RNA (Vasudevan et al., 2013), but not double-stranded
RNA. The efcacy of the different knockdowns (KDs) was validated by western blotting and qRT-PCR (Figure S2C). Total
RNA extracted from the different experiments was depleted
from rRNAs and sequenced (Supplemental Experimental Procedures). We reproducibly observed that ADAR depletion resulted
in signicantly higher (p value < 2.2 3 1016) circRNA expression
compared with controls (Figure S2D), whereas the linear host
transcripts were less strongly affected (Figure 4C). For example,
84 circRNAs and 11 linear hosts were upregulated more than 2fold (Fishers exact test, p value 5.7 3 1013). This effect was
seen in independent biological replicates, as well as in a comparison of ADAR1/2 KD with APOBEC3 KD (Figures S2E and S2F).
We set out to validate these observations in independently
carried out ADAR1 and ADAR2 KD experiments (using the previous and an independent small interfering RNAs [siRNAs]
against ADAR1) followed by qRT-PCR assays (Experimental
Procedures; Figures 4D, 4E, and S2G). CircRNA candidates
were selected based on the increased fold changes observed
in sequencing data sets. circRNA expression was compared
with the expression of the respective linear host gene. Four
out of eight circRNAs (UBAC2, SMARCA, SPECC1, and
HIPK2) were upregulated upon ADAR1 depletion, whereas
the linear host RNAs did not show consistent expression
changes. PUM1 and CREBBP circRNA were upregulated to
the same level as their linear host transcripts. Consistent
with sequencing data, CDR1as was expressed 4-fold higher
in ADAR1 KD samples compared with control. Two circRNAs
(GAPVD1 and PDS5B) did not show the expression changes
observed in the sequencing data (Figures 4E and S2G). Based
on these ndings, we conclude that ADAR1 depletion can
induce upregulation of circRNAs independently of the expression level of the linear host circRNA.
Cell Reports 10, 170177, January 13, 2015 2015 The Authors 173
To explore whether circRNA anking introns undergo extensive hyperediting (Carmi et al., 2011), we applied a computational pipeline that detects hyperediting events in the RNA
sequencing data sets used by Memczak et al. (2013) to predict
circRNAs, as well as in ADAR KD and control data sets. This
analysis (for details, see Porath et al. 2014) identied 165,000
unique hyperediting sites in the human genome (72.000,
49.000, 45.000, and 11.000 in Memczak et al. [2013] and
two control and ADAR KD samples, respectively). Notably, the
number of hyperedited sites identied in the ADAR KD sample
was 4- to 5-fold smaller compared with controls. We found
that 25% (19%) of circRNA upstream (downstream) introns
had at least one conversion, whereas only 6% of other introns
from the same genes had conversions. The number of conversions per base was also found to be elevated for the upstream
174 Cell Reports 10, 170177, January 13, 2015 2015 The Authors
(downstream) introns (Figures S3AS3C). A similar highly signicant enrichment of hyperediting events was detected for
C. elegans circRNA upstream introns (Supplemental Experimental Procedures; Figure S3D). Together, our data suggest
that A-to-I editing by ADAR is a universal hallmark of circRNA
biogenesis in animals.
DISCUSSION
Our analyses support a model in which RCMs between introns
that bracket an exon promote the circularization of that exon.
We note that this model is powerful enough to enable us to successfully predict and experimentally validate circRNAs. Our data
suggest that this model of circRNA biogenesis is conserved
across animals.
Figure 4. ADAR
Expression
Antagonizes
circRNA
(2013) with the additional ltering step of requiring unique alignments for
both of the read anchors prior to extension. Thereafter, we retained circRNAs
that overlap annotated internal splice sites.
Intron Alignments
We carried out intron alignments using BLAST (Altschul et al., 1990) with
the parameters -task blastn -word_size 6 for C. elegans and -task
blastn -word_size 11 for human. For further analysis, we considered
only alignments that exceeded BLAST score cutoffs of 20 and 100 for
C. elegans and human, respectively. For each intron pair, we calculated
circularization H as a top BLAST score multiplied by the probability of
forming at least one stem around the exon (Supplemental Experimental
Procedures).
Nonrepetitive RCMs were dened as BLAST matches that did not overlap
with any sequence present in UCSC RepeatMasker. Putative circRNAs with
nonrepetitive RCMs are shown in Table S3.
Conservation Analysis
To determine homologous exon groups, we used UCSC liftOver with minMatch = 0.1 option. We found that 71 out of 1,067 human circRNAs (Memczak
et al., 2013) overlapped annotated splice sites and were present in circular
form in mouse. We compared H score distributions using the Mann-Whitney
two-sided test.
qRT-PCR
We performed qPCR using the Maxima SYBR-Green/ROX qPCR master
mix (Thermo Scientic) and a StepOnePlus PCR system (Applied Biosystems). To detect putative head-to-tail junctions, we designed divergent
primers for each circRNA candidate. Ct values for mock/RNase R-treated
circRNAs were normalized to C. elegans spike-in RNA (for standard curves,
see Figures S5AS5C and Table S2). Amplicons were gel or bead puried
(Zymoclean gel DNA recovery kit [Zymo Research]; Agencourt AMPure
XP [Beckman Coulter]) and subjected to Sanger sequencing by LGC Genomics. Conrmed head-to-tail junctions are available in Figure S5D. A list of
the oligos is given in Table S1.
For more details, see the Supplemental Experimental Procedures.
ACCESSION NUMBERS
The data have been deposited in the NCBI Gene Expression Omnibus and are
available under accession numbers GSE63823 and SRP050149. circRNAs
identied in this study, as well as H scores, are available at http://www.
circbase.org.
SUPPLEMENTAL INFORMATION
Supplemental Information includes Supplemental Experimental Procedures,
ve gures, and three tables and can be found with this article online at
http://dx.doi.org/10.1016/j.celrep.2014.12.019.
CONCLUSIONS
AUTHOR CONTRIBUTIONS
A.I. performed all computational analyses except hyperediting. S.M. contributed experimental validation assays. E.W. designed and carried out ADAR
RNAi, RNA-seq, and qPCR experiments and initial analysis of the RNA-seq
data. H.P. and E.L. performed hyperediting detection. F.T. collected and prepared worm samples, carried out initial ADAR KDs, and helped in the design of
the validation experiments. M.O. annotated worm circRNAs. M.P. performed
an initial analysis of the RNA-seq data and was supervised by C.D. A.I. and
N.R. analyzed the data and wrote the paper with input from E.L.
EXPERIMENTAL PROCEDURES
ACKNOWLEDGMENTS
A.I. and N.R. thank Sebastian Kadener (Hebrew University), Albrecht Bindereif
(University of Giessen), and Marvin Jens (N.R. lab) for helpful discussions. We
thank all members of the N.R. lab for discussions and support. This work was
176 Cell Reports 10, 170177, January 13, 2015 2015 The Authors
supported by German-Israeli-Foundation for Scientic Research and Development (G.I.F) and German Ministry for Education and Research (SatNet
program).
Levanon, E.Y., Eisenberg, E., Yelin, R., Nemzer, S., Hallegger, M., Shemesh,
R., Fligelman, Z.Y., Shoshan, A., Pollock, S.R., Sztybel, D., et al. (2004). Systematic identication of abundant A-to-I editing sites in the human transcriptome. Nat. Biotechnol. 22, 10011005.
Liang, D., and Wilusz, J.E. (2014). Short intronic repeat sequences facilitate
circular RNA production. Genes Dev. 28, 22332247.
REFERENCES
Altschul, S.F., Gish, W., Miller, W., Myers, E.W., and Lipman, D.J. (1990). Basic
local alignment search tool. J. Mol. Biol. 215, 403410.
Ashwal-Fluss, R., Meyer, M., Pamudurti, N.R., Ivanov, A., Bartok, O., Hanan,
M., Evantal, N., Memczak, S., Rajewsky, N., and Kadener, S. (2014). circRNA
biogenesis competes with pre-mRNA splicing. Mol. Cell 56, 5566.
Athanasiadis, A., Rich, A., and Maas, S. (2004). Widespread A-to-I RNA editing
of Alu-containing mRNAs in the human transcriptome. PLoS Biol. 2, e391.
Baltz, A.G., Munschauer, M., Schwanhausser, B., Vasile, A., Murakawa, Y.,
Schueler, M., Youngs, N., Penfold-Brown, D., Drew, K., Milek, M., et al.
(2012). The mRNA-bound proteome and its global occupancy prole on protein-coding transcripts. Mol. Cell 46, 674690.
Braun, S., Domdey, H., and Wiebauer, K. (1996). Inverse splicing of a discontinuous pre-mRNA intron generates a circular exon in a HeLa cell nuclear
extract. Nucleic Acids Res. 24, 41524157.
Memczak, S., Jens, M., Elefsinioti, A., Torti, F., Krueger, J., Rybak, A., Maier,
L., Mackowiak, S.D., Gregersen, L.H., Munschauer, M., et al. (2013). Circular
RNAs are a large class of animal RNAs with regulatory potency. Nature 495,
333338.
Morse, D.P., Aruscavage, P.J., and Bass, B.L. (2002). RNA hairpins in noncoding regions of human brain and Caenorhabditis elegans mRNA are edited by
adenosine deaminases that act on RNA, Volume 99 (USA: Proc. Natl. Acad.
Sci), pp. 79067911.
Nishikura, K. (2010). Functions and regulation of RNA editing by ADAR deaminases. Annu. Rev. Biochem. 79, 321349.
Osenberg, S., Paz Yaacov, N., Safran, M., Moshkovitz, S., Shtrichman, R.,
Sherf, O., Jacob-Hirsch, J., Keshet, G., Amariglio, N., Itskovitz-Eldor, J., and
Rechavi, G. (2010). Alu sequences in undifferentiated human embryonic
stem cells display high levels of A-to-I RNA editing. PLoS ONE 5, e11173.
Pasman, Z., Been, M.D., and Garcia-Blanco, M.A. (1996). Exon circularization
in mammalian nuclear extracts. RNA 2, 603610.
Porath, H.T., Carmi, S., and Levanon, E.Y. (2014). A genome-wide map of hyper-edited RNA reveals numerous new sites. Nat. Commun. 5, 4726.
Capel, B., Swain, A., Nicolis, S., Hacker, A., Walter, M., Koopman, P., Goodfellow, P., and Lovell-Badge, R. (1993). Circular transcripts of the testis-determining gene Sry in adult mouse testis. Cell 73, 10191030.
Ramaswami, G., and Li, J.B. (2014). RADAR: a rigorously annotated database
of A-to-I RNA editing. Nucleic Acids Res. 42, D109D113.
Carmi, S., Borukhov, I., and Levanon, E.Y. (2011). Identication of widespread
ultra-edited human RNAs. PLoS Genet. 7, e1002317.
Ramaswami, G., Lin, W., Piskol, R., Tan, M.H., Davis, C., and Li, J.B. (2012).
Accurate identication of human Alu and non-Alu RNA editing sites. Nat.
Methods 9, 579581.
Dubin, R.A., Kazmi, M.A., and Ostrer, H. (1995). Inverted repeats are necessary
for circularization of the mouse testis Sry transcript. Gene 167, 245248.
Glazar, P., Papavasileiou, P., and Rajewsky, N. (2014). circBase: a database
for circular RNAs. RNA 20, 16661670.
Salzman, J., Gawad, C., Wang, P.L., Lacayo, N., and Brown, P.O. (2012). Circular RNAs are the predominant transcript isoform from hundreds of human
genes in diverse cell types. PLoS ONE 7, e30733.
Hansen, T.B., Wiklund, E.D., Bramsen, J.B., Villadsen, S.B., Statham, A.L.,
Clark, S.J., and Kjems, J. (2011). miRNA-dependent gene silencing involving
Ago2-mediated cleavage of a circular antisense RNA. EMBO J 30, 44144422.
Starke, S., Jost, I., Rossbach, O., Schneider, T., Schreiner, S., Hung, L.-H., and
Bindereif, A. (2014). Exon circularization requires canonical splice signals. Cel
Rep 10, this issue.
Hansen, T.B., Jensen, T.I., Clausen, B.H., Bramsen, J.B., Finsen, B., Damgaard, C.K., and Kjems, J. (2013). Natural RNA circles function as efcient
microRNA sponges. Nature 495, 384388.
Vasudevan, A.A.J., Smits, S.H.J., Hoppner, A., Haussinger, D., Koenig, B.W.,
and Munk, C. (2013). Structural features of antiviral DNA cytidine deaminases.
Biol. Chem. 394, 13571370.
Jeck, W.R., and Sharpless, N.E. (2014). Detecting and characterizing circular
RNAs. Nat. Biotechnol. 32, 453461.
Wang, P.L., Bao, Y., Yee, M.-C., Barrett, S.P., Hogan, G.J., Olsen, M.N., Dinneny, J.R., Brown, P.O., and Salzman, J. (2014). Circular RNA is expressed
across the eukaryotic tree of life. PLoS ONE 9, e90859.
Jeck, W.R., Sorrentino, J.A., Wang, K., Slevin, M.K., Burd, C.E., Liu, J., Marzluff, W.F., and Sharpless, N.E. (2013). Circular RNAs are abundant,
conserved, and associated with ALU repeats. RNA 19, 141157.
Zhang, X.-O., Wang, H.-B., Zhang, Y., Lu, X., Chen, L.-L., and Yang, L. (2014).
Complementary sequence-mediated exon circularization. Cell 159, 134147.
Cell Reports 10, 170177, January 13, 2015 2015 The Authors 177
Cell Reports
Report
ThermoMouse: An In Vivo Model
to Identify Modulators of UCP1 Expression
in Brown Adipose Tissue
Andrea Galmozzi,2,3 Si B. Sonne,1,3,4 Svetlana Altshuler-Keylin,1 Yutaka Hasegawa,1 Kosaku Shinoda,1
Ineke H.N. Luijten,1,5 Jae Won Chang,2 Louis Z. Sharp,1 Benjamin F. Cravatt,2 Enrique Saez,2,* and Shingo Kajimura1,*
1UCSF
Diabetes Center, Department of Cell and Tissue Biology, University of California, San Francisco, 35 Medical Center Way,
San Francisco, CA 94143, USA
2Department of Chemical Physiology and The Skaggs Institute for Chemical Biology, The Scripps Research Institute, La Jolla, CA 92037, USA
3Co-rst author
4Present address: Department of Biology, University of Copenhagen, Copenhagen 2200, Denmark
5Present address: The Wenner-Gren Institute, The Arrhenius Laboratories, Stockholm University, Stockholm 106 91, Sweden
*Correspondence: [email protected] (E.S.), [email protected] (S.K.)
http://dx.doi.org/10.1016/j.celrep.2014.10.066
This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/3.0/).
SUMMARY
1584 Cell Reports 9, 15841593, December 11, 2014 2014 The Authors
Cell Reports 9, 15841593, December 11, 2014 2014 The Authors 1585
1586 Cell Reports 9, 15841593, December 11, 2014 2014 The Authors
two serine hydrolases expressed in adipocytes, carboxylesterase 3 (Ces3 or Ces1d) and Ces1f (CesML1). However, the ability
of WWL113 to induce UCP1 expression in brown adipocytes
does not appear to be mediated by Ces3/1f inhibition, for a
structurally distinct Ces3 inhibitor (WWL229) failed to have the
same effect on UCP1 expression, whereas the urea version of
WWL113 (WWL113U), which does not inhibit Ces3, retained
the ability to induce UCP1 (Figure S3A). Nonetheless, we
concluded that WWL113 could be a valuable chemical probe
to further validate our reporter systems.
Induction of UCP1 Expression by WWL113 Relies on
PPARa Signaling
To characterize the effect of WWL113 on endogenous UCP1
expression, we treated cultured brown adipocytes with several
doses of WWL113. At a dose of 1 mM and higher, WWL113 signicantly increased UCP1 protein expression (Figure 4A). The
effect of WWL113 on UCP1 protein levels was largely due to
activation of Ucp1 transcription, because WWL113 powerfully
increased expression of Ucp1 mRNA (Figure 4B). WWL113 treatment also stimulated expression of mRNAs for other thermogenic genes, such as Cidea, Pgc1a, and Cox7a1 (Figure 4B).
WWL113 treatment had no effect on expression of the adipogenic marker Adiponectin. WWL113 treatment for 24 hr was sufcient to activate endogenous Ucp1 mRNA expression without
affecting expression of multiple adipogenic markers (Figures
4C and S3B), indicating that WWL113 enhanced the BAT-selective thermogenic gene program in a cell-autonomous manner
without affecting adipogenesis per se. To test the functional consequences of WWL113 treatment, we examined the extent to
which WWL113 could sensitize brown adipocytes to physiologic
activators of BAT such as norepinephrine. WWL113 pretreatment enhanced the increase in Ucp1 mRNA expression normally
induced by norepinephrine (Figure 4D). More importantly,
WWL113-treated cells showed greater total and uncoupled
respiration in response to norepinephrine and greater uncoupled
basal respiration (Figures 4E and S3C).
Next, we explored the mechanism by which WWL113
increased Ucp1 transcription. Because PPARa is a critical regulator of thermogenic gene expression in BAT (Barbera et al.,
2001), we hypothesized that the action of WWL113 in BAT
could be mediated via PPARa. To test this notion, primary brown
adipocytes were treated with a selective PPARa antagonist
(GW6471) and/or a PPARg antagonist (GW9662) in the presence
or absence of WWL113. The capacity of WWL113 to increase
Ucp1 mRNA expression was largely blunted by treatment with
GW6471, but not with GW9662 (Figure 4F). The inhibitory effect
of GW6471 on WWL113-induced Ucp1 mRNA expression was
not altered when cells were cotreated with GW9662. These
results indicate that the ability of WWL113 to enhance Ucp1
expression is principally dependent on PPARa, but not PPARg,
activity. Because WWL113 is not a direct PPARa activator
(Dominguez et al., 2014), we tested the extent to which cotreatment of differentiated brown adipocytes with WWL113 and a
PPARa agonist (GW9578) would further boost Ucp1 mRNA
levels. Cotreatment with WWL113 and GW9578 modestly but
signicantly increased expression of Ucp1 relative to treatment
with either compound alone (Figure 4G), suggesting that
WWL113 cooperates with the PPARa pathway to enhance
UCP1 expression.
WWL113 Increases UCP1 Expression and the
Thermogenic Response In Vivo
We next tested if WWL113 could enhance UCP1 expression
in vivo. Ucp1 luciferase mice were treated with vehicle or
WWL113 (50 mg/kg once daily) for 5 days. This dose has been
shown to be effective in mice (Dominguez et al., 2014). A robust
and signicant increase (5-fold) in Ucp1-driven luciferase
expression was detected in the interscapular BAT depots of
transgenic mice treated with WWL113 (Figures 5A and 5B).
To conrm these ndings, we examined the ability of WWL113
to enhance thermogenic gene expression in C57BL/6J mice
treated with the compound for 5 days. WWL113 treatment
induced signicant increases in mRNA expression of Ucp1 and
thermogenic genes, such as Pgc1a and Dio2, in the BAT of
wild-type mice (Figure 5C). No change in the general marker
Pparg was seen. Importantly, UCP1 protein expression was
highly induced in vivo by WWL113 (Figure 5D). In contrast, no difference in mRNA expression of Ucp1, Pgc1a, Dio2, and Pparg
was observed in inguinal WAT of WWL113-treated mice (Figure S5A), indicating that the effects of WWL113 on the thermogenic gene program may be specic to brown fat.
Finally, we examined the effects of WWL113 on whole-body
energy expenditure. Mice treated with WWL113 for 7 days
showed no differences in basal energy expenditure, but upon
adrenergic stimulation (CL-316,243 injection), they responded
with a more robust increase in energy expenditure than controls
(Figure 5E). WWL113 did not affect locomotor activity (Figure 5F),
food intake (Figure 5G), or heart rate (Figure 5H). These data indicate that WWL113 increased UCP1 levels in vivo to a functionally
meaningful degree, increasing the adaptive thermogenic capacity of treated mice.
Cell Reports 9, 15841593, December 11, 2014 2014 The Authors 1587
1588 Cell Reports 9, 15841593, December 11, 2014 2014 The Authors
DISCUSSION
We have developed an in vivo monitoring system that allows to
quantitatively track sequential changes in UCP1 expression
within the same individual. Because UCP1 expression shows a
high level of interindividual variation (Boeuf et al., 2002) and
changes dynamically during circadian oscillation (Gerhart-Hines
et al., 2013), seasonal changes (Au-Yong et al., 2009), and
aging (Rogers et al., 2012), this model may prove of wide utility.
18
uoro-labeled 2-deoxy-glucose positron emission tomography (18FDG-PET) scanning has been applied to assess BAT activity in rodents and humans (Cypess et al., 2009; Saito et al.,
2009; van Marken Lichtenbelt et al., 2009; Virtanen et al.,
2009), but detection of 18FDG-PET signals in BAT depends
entirely on glucose uptake. In contrast, the level of UCP1 expression in BAT is a more-direct measure of the thermogenic capacity of this tissue (Nedergaard and Cannon, 2010). Hence, the
transgenic UCP1 reporter we describe provides an opportunity
to identify signaling pathways and transcriptional events that
control thermogenic capacity in brown adipocytes in vivo. Unlike
prior UCP1 models (Cassard-Doulcier et al., 1993), it also en-
Cell Reports 9, 15841593, December 11, 2014 2014 The Authors 1589
differences in locomotor activity, food intake, or heartbeat, indicating that the compound does not enhance basal sympathetic
drive. Together with the nding that WWL113 increases UCP1
expression in brown fat cells in a cell-autonomous manner, these
results suggest that WWL113 boosts energy expenditure primarily
by increasing the content of UCP1 in BAT that can be activated by
physiologic stimuli such as cold.
We previously reported that chronic administration (2 months)
of WWL113 to obese-diabetic mice reduced body weight gain,
improved systemic glucose and lipid homeostasis, and cleared
hepatic steatosis (Dominguez et al., 2014). We ascribed the
effects of WWL113 to inhibition of Ces3 in WAT and liver. Using
a distinct phenotypic screen as a starting point, in this study, we
have found that WWL113 can have Ces3-independent effects in
1590 Cell Reports 9, 15841593, December 11, 2014 2014 The Authors
serve as a tool to monitor the effects of compounds and biological factors on beige cells.
Although they share the thermogenic ability of brown adipocytes, beige adipocytes have a distinct, heterogeneous developmental origin that is not fully understood. Beige adipocytes in
multiple WAT depots can arise from Myf5-positive and negative
cells (Sanchez-Gurmaches and Guertin, 2014), and a subset of
inguinal beige adipocytes originates from a smooth-muscle
lineage (Long et al., 2014). Prior work has shown that regulation
of Ucp1 is distinct in the two types of thermogenic adipocytes
(Guerra et al., 1998; Koza et al., 2000; Xue et al., 2007). We found
that WWL113 could activate UCP1 expression in BAT, but
not in inguinal WAT, implying that WWL113-initiated signaling
events that regulate UCP1 expression may be unique to brown
adipocytes.
BAT and liver are the major organs responsible for triglyceride
uptake (Bartelt et al., 2011), and emerging evidence points to a
close connection between BAT activity, liver function, and systemic lipid homeostasis. For example, defects in thermogenesis
caused by BAT-specic knockout of the enzyme euchromatic
histone-lysine N-methyltransferase 1 resulted in hepatic steatosis and insulin resistance, even when weight-matched mice
Cell Reports 9, 15841593, December 11, 2014 2014 The Authors 1591
2005). Ucp1 luciferase BAC DNA was microinjected into single-cell FVB
embryos and transgenic founders and their offspring identied by PCR
(primers provided in Table S1). Segregation patterns indicate that the transgene inserted into the Y chromosome. Transgenics display no decrease in
fertility or any other abnormalities.
In Vivo Luciferase Imaging
Luciferase activity was monitored using an IVIS Spectrum Instrument (Caliper
Life Sciences). For 3D reconstruction, six images (exposure time: 180 s;
binning: L; F/Stop:1; emission lters: 560660 nm; eld of view: C) were
collected starting 8 min after injection of 150 mg/kg luciferin (Goldbio). In other
experiments, one image per scan (exposure time: 300 s; binning: M; F/Stop:1;
emission lters: open; eld of view: D) was acquired 15 min after luciferin injection. 3D reconstructions and luciferase activity were calculated using Living
Image Software (Caliper Life Sciences).
Immortalized Ucp1 Luciferase Adipocyte Lines
The stromal vascular fraction from interscapular BAT of 3-week-old male
UCP1 luciferase transgenic mice was isolated (Ohno et al., 2012) and cells
infected with a retrovirus expressing large T antigen (pBabe SV40 Large T antigen; Addgene) and selected in puromycin (2 mg/ml). Immortalized preadipocytes were cultured in Dulbeccos modied Eagles medium (DMEM), 10%
fetal bovine serum (FBS), penicillin, and streptomycin. Upon 100% conuence
(day 0), differentiation was induced with medium containing 10% FBS, 5 mg/ml
insulin, 1 nM T3, 0.125 mM indomethacin, 2 mg/ml dexamethasone, and
0.125 mM 3-isobutyl-1-methylxanthine for 2 days. From day 2 on, cells were
cultured only in the presence of insulin and T3.
Phenotypic Screen
A library of approximately 300 compounds (Adibekian et al., 2011; Bachovchin
et al., 2010) was screened. Immortalized Ucp1 luciferase brown preadipocytes
were seeded in 96-well plates and differentiated as described above. At day 8,
cells were exposed to compounds (5 mM for carbamates, 2.5 mM for triazole
ureas) in serum-free DMEM. Rosiglitazone (0.5 mM) served as positive control.
After 16 hr, media was replaced with Glo Lysis Buffer (Promega) and luciferase
activity quantied in a PHERAstar reader (BMG Labtech). Activity was normalized relative to the signal in DMSO-treated cells in each plate. Assays were
performed in triplicate. Compounds inducing >50% increase or decrease in
luciferase activity were selected for follow-up.
Luciferase Assays
Cells or tissues were lysed in Cell Culture Lysis Reagent (Promega) and luciferase activity quantied using the Luciferase Assay System (Promega). Activity
was measured in an Optocomp I reader (MGM Instruments) and normalized to
total protein content.
Gene Expression
Total RNA was isolated using RiboZol reagents (AMRESCO) and reversed
transcribed using an iScript cDNA synthesis kit (Bio-Rad) and quantitative
real-time PCR performed using SYBR green and an ABI ViiA7 machine. Relative mRNA expression was determined by the DD-Ct method with TATAbinding protein as a normalization control. Primer sequences are provided in
Table S1.
Metabolic Studies
Whole-body energy expenditure was measured at ambient temperature using
a Comprehensive Lab Animal Monitoring System (Columbus Instruments) at
the University of California, San Francisco (UCSF) animal metabolic core facility. Wild-type mice were treated daily with WWL113 (50 mg/kg) or vehicle for
7 days (n = 6). The mice were injected intraperitoneally with a b3-adrenergic
receptor-specic agonist CL-316,243 at a dose of 1 mg/kg to examine the
response to adrenergic stimulation. VO2 was normalized by body weight.
Western Analysis
Proteomes were separated by SDS-PAGE. Antibodies used: UCP1 (Abcam;
ab10983); Akt (Cell Signaling Technology; 9272); a-tubulin (Sigma; T8203);
and b-actin (Sigma-Aldrich; AC-15).
Cellular Respiration
Oxygen consumption rate was measured in a MT200A Cell Respirometer
(Strathkelvin), as previously described (Kajimura et al., 2009). Briey, differentiated brown adipocytes treated with DMSO or WWL113 (10 mM) for 48 hr were
trypsinized and incubated in serum-free medium in the presence or absence of
norepinephrine. Uncoupled and nonmitochondrial cellular respiration were
measured using oligomycin (1 mM) and antimycin A (1 mM).
Fat Transplantation
Preadipocytes were implanted as described (Kajimura et al., 2009). Briey,
cultured immortalized Ucp1 luciferase preadipocytes were trypsinized,
washed, and resuspended in PBS. Preadipocytes in a volume of 300 ml (4
3 107 cells) were injected subcutaneously into NCr nude mice (Taconic).
Six days after transplantation, mice were injected with either saline (n = 4) or
rosiglitazone (n = 6; 10 mg/kg) twice daily for 7 days. Luciferase activity was
monitored on days 0, 1, 4, and 7.
Histology
Tissues were xed in 4% paraformaldehyde and embedded in parafn. Sections (7 mm) were analyzed with hematoxylin and eosin staining and immunouorescence to detect UCP1 (ab10983 1:1,000 in 5% goat serum) or GFP
(GFP-1020 1:500 in 5% goat serum) as described (Sharp et al., 2012).
Statistics
Statistical signicance was dened as p < 0.05 and determined by two-tailed
Student t tests, Wilcoxon, or ANOVA, with Dunnes multiple comparison post
hoc analysis.
SUPPLEMENTAL INFORMATION
Supplemental Information includes four gures and one table and can be
found with this article online at http://dx.doi.org/10.1016/j.celrep.2014.10.066.
AUTHOR CONTRIBUTIONS
A.G., S.B.S., E.S., and S.K. designed experiments. J.W.C. and B.F.C. provided
libraries and synthesized compounds. S.K. and S.B.S. generated transgenic
model and cell lines. A.G. performed screen and follow-up. A.G., S.B.S.,
S.A.-K., Y.H., K.S., and S.K. performed in vivo and vitro experiments. L.Z.S.
and I.H.N.L. provided technical help. A.G., S.B.S., E.S., and S.K. wrote the
manuscript.
ACKNOWLEDGMENTS
We thank Haemin Hong, Kathleen Jay, Christophe Paillart, and Denis Glenn for
assistance. Work was supported by NIH grants DK087853 and DK97441 to
S.K. and DK099810 and CA179489 to E.S. S.K. acknowledges support from
the DERC center grant (DK63720), UCSF PBBR program, the Pew Charitable
Trust, and PRESTO from Japan Science and Technology Agency. S.B.S. was
supported by a fellowship from the Alfred Benzon Foundation and A.G. by
fellowship 14POST18200019 from the American Heart Association.
Received: April 2, 2014
Revised: October 13, 2014
Accepted: October 24, 2014
Published: November 26, 2014
REFERENCES
Adibekian, A., Martin, B.R., Wang, C., Hsu, K.L., Bachovchin, D.A., Niessen,
S., Hoover, H., and Cravatt, B.F. (2011). Click-generated triazole ureas as
ultrapotent in vivo-active serine hydrolase inhibitors. Nat. Chem. Biol. 7,
469478.
Au-Yong, I.T., Thorn, N., Ganatra, R., Perkins, A.C., and Symonds, M.E. (2009).
Brown adipose tissue and seasonal variation in humans. Diabetes 58, 2583
2587.
1592 Cell Reports 9, 15841593, December 11, 2014 2014 The Authors
Bachovchin, D.A., Ji, T., Li, W., Simon, G.M., Blankman, J.L., Adibekian, A.,
Hoover, H., Niessen, S., and Cravatt, B.F. (2010). Superfamily-wide portrait
of serine hydrolase inhibition achieved by library-versus-library screening.
Proc. Natl. Acad. Sci. USA 107, 2094120946.
Barbera, M.J., Schluter, A., Pedraza, N., Iglesias, R., Villarroya, F., and Giralt,
M. (2001). Peroxisome proliferator-activated receptor alpha activates transcription of the brown fat uncoupling protein-1 gene. A link between regulation
of the thermogenic and lipid oxidation pathways in the brown fat cell. J. Biol.
Chem. 276, 14861493.
Kopecky, J., Clarke, G., Enerback, S., Spiegelman, B., and Kozak, L.P. (1995).
Expression of the mitochondrial uncoupling protein gene from the aP2 gene
promoter prevents genetic obesity. J. Clin. Invest. 96, 29142923.
Koza, R.A., Hohmann, S.M., Guerra, C., Rossmeisl, M., and Kozak, L.P. (2000).
Synergistic gene interactions control the induction of the mitochondrial uncoupling protein (Ucp1) gene in white fat tissue. J. Biol. Chem. 275, 3448634492.
Lidell, M.E., Betz, M.J., Dahlqvist Leinhard, O., Heglind, M., Elander, L., Slawik,
M., Mussack, T., Nilsson, D., Romu, T., Nuutila, P., et al. (2013). Evidence for
two types of brown adipose tissue in humans. Nat. Med. 19, 631634.
Bartelt, A., Bruns, O.T., Reimer, R., Hohenberg, H., Ittrich, H., Peldschus, K.,
Kaul, M.G., Tromsdorf, U.I., Weller, H., Waurisch, C., et al. (2011). Brown adipose tissue activity controls triglyceride clearance. Nat. Med. 17, 200205.
Long, J.Z., Svensson, K.J., Tsai, L., Zeng, X., Roh, H.C., Kong, X., Rao, R.R.,
Lou, J., Lokurkar, I., Baur, W., et al. (2014). A smooth muscle-like origin for
beige adipocytes. Cell Metab. 19, 810820.
Boeuf, S., Keijer, J., Franssen-Van Hal, N.L., and Klaus, S. (2002). Individual
variation of adipose gene expression and identication of covariated genes
by cDNA microarrays. Physiol. Genomics 11, 3136.
Nedergaard, J., and Cannon, B. (2010). The changed metabolic world with human brown adipose tissue: therapeutic visions. Cell Metab. 11, 268272.
Cassard-Doulcier, A.M., Gelly, C., Fox, N., Schrementi, J., Raimbault, S.,
Klaus, S., Forest, C., Bouillaud, F., and Ricquier, D. (1993). Tissue-specic
and b-adrenergic regulation of the mitochondrial uncoupling protein gene:
control by cis-acting elements in the 50 -anking region. Mol. Endocrinol. 7,
497506.
Ohno, H., Shinoda, K., Spiegelman, B.M., and Kajimura, S. (2012). PPARg agonists induce a white-to-brown fat conversion through stabilization of
PRDM16 protein. Cell Metab. 15, 395404.
Ohno, H., Shinoda, K., Ohyama, K., Sharp, L.Z., and Kajimura, S. (2013).
EHMT1 controls brown adipose cell fate and thermogenesis through the
PRDM16 complex. Nature 504, 163167.
Cypess, A.M., Lehman, S., Williams, G., Tal, I., Rodman, D., Goldne, A.B.,
Kuo, F.C., Palmer, E.L., Tseng, Y.H., Doria, A., et al. (2009). Identication
and importance of brown adipose tissue in adult humans. N. Engl. J. Med.
360, 15091517.
Ortega-Molina, A., Efeyan, A., Lopez-Guadamillas, E., Munoz-Martin, M., Gomez-Lopez, G., Canamero, M., Mulero, F., Pastor, J., Martinez, S., Romanos,
E., et al. (2012). Pten positively regulates brown adipose function, energy
expenditure, and longevity. Cell Metab. 15, 382394.
Cypess, A.M., White, A.P., Vernochet, C., Schulz, T.J., Xue, R., Sass, C.A.,
Huang, T.L., Roberts-Toler, C., Weiner, L.S., Sze, C., et al. (2013). Anatomical
localization, gene expression proling and functional characterization of adult
human neck brown fat. Nat. Med. 19, 635639.
Rogers, N.H., Landa, A., Park, S., and Smith, R.G. (2012). Aging leads to a programmed loss of brown adipocytes in murine subcutaneous white adipose tissue. Aging Cell 11, 10741083.
Dominguez, E., Galmozzi, A., Chang, J.W., Hsu, K.L., Pawlak, J., Li, W., Godio,
C., Thomas, J., Partida, D., Niessen, S., et al. (2014). Integrated phenotypic
and activity-based proling links Ces3 to obesity and diabetes. Nat. Chem.
Biol. 10, 113121.
Enerback, S., Jacobsson, A., Simpson, E.M., Guerra, C., Yamashita, H.,
Harper, M.E., and Kozak, L.P. (1997). Mice lacking mitochondrial uncoupling
protein are cold-sensitive but not obese. Nature 387, 9094.
Fedorenko, A., Lishko, P.V., and Kirichok, Y. (2012). Mechanism of fatty-aciddependent UCP1 uncoupling in brown fat mitochondria. Cell 151, 400413.
Feldmann, H.M., Golozoubova, V., Cannon, B., and Nedergaard, J. (2009).
UCP1 ablation induces obesity and abolishes diet-induced thermogenesis in
mice exempt from thermal stress by living at thermoneutrality. Cell Metab. 9,
203209.
Gerhart-Hines, Z., Feng, D., Emmett, M.J., Everett, L.J., Loro, E., Briggs, E.R.,
Bugge, A., Hou, C., Ferrara, C., Seale, P., et al. (2013). The nuclear receptor
Rev-erba controls circadian thermogenic plasticity. Nature 503, 410413.
Golozoubova, V., Hohtola, E., Matthias, A., Jacobsson, A., Cannon, B., and
Nedergaard, J. (2001). Only UCP1 can mediate adaptive nonshivering thermogenesis in the cold. FASEB J. 15, 20482050.
Grundlingh, J., Dargan, P.I., El-Zanfaly, M., and Wood, D.M. (2011). 2,4-dinitrophenol (DNP): a weight loss agent with signicant acute toxicity and risk
of death. J. Med. Toxicol. 7, 205212.
Saito, M., Okamatsu-Ogura, Y., Matsushita, M., Watanabe, K., Yoneshiro, T.,
Nio-Kobayashi, J., Iwanaga, T., Miyagawa, M., Kameya, T., Nakada, K., et al.
(2009). High incidence of metabolically active brown adipose tissue in healthy
adult humans: effects of cold exposure and adiposity. Diabetes 58, 1526
1531.
Sanchez-Gurmaches, J., and Guertin, D.A. (2014). Adipocyte lineages: tracing
back the origins of fat. Biochim. Biophys. Acta 1842, 340351.
Sears, I.B., MacGinnitie, M.A., Kovacs, L.G., and Graves, R.A. (1996). Differentiation-dependent expression of the brown adipocyte uncoupling protein
gene: regulation by peroxisome proliferator-activated receptor gamma. Mol.
Cell. Biol. 16, 34103419.
Sharp, L.Z., Shinoda, K., Ohno, H., Scheel, D.W., Tomoda, E., Ruiz, L., Hu, H.,
Wang, L., Pavlova, Z., Gilsanz, V., and Kajimura, S. (2012). Human BAT possesses molecular signatures that resemble beige/brite cells. PLoS ONE 7,
e49452.
van Marken Lichtenbelt, W.D., Vanhommerig, J.W., Smulders, N.M., Drossaerts, J.M., Kemerink, G.J., Bouvy, N.D., Schrauwen, P., and Teule, G.J.
(2009). Cold-activated brown adipose tissue in healthy men. N. Engl. J. Med.
360, 15001508.
Virtanen, K.A., Lidell, M.E., Orava, J., Heglind, M., Westergren, R., Niemi, T.,
Taittonen, M., Laine, J., Savisto, N.J., Enerback, S., and Nuutila, P. (2009).
Functional brown adipose tissue in healthy adults. N. Engl. J. Med. 360,
15181525.
Guerra, C., Koza, R.A., Yamashita, H., Walsh, K., and Kozak, L.P. (1998).
Emergence of brown adipocytes in white fat in mice is under genetic control.
Effects on body weight and adiposity. J. Clin. Invest. 102, 412420.
Waki, H., Park, K.W., Mitro, N., Pei, L., Damoiseaux, R., Wilpitz, D.C., Reue, K.,
Saez, E., and Tontonoz, P. (2007). The small molecule harmine is an antidiabetic cell-type-specic regulator of PPARgamma expression. Cell Metab. 5,
357370.
Kajimura, S., and Saito, M. (2014). A new era in brown adipose tissue biology:
molecular control of brown fat development and energy homeostasis. Annu.
Rev. Physiol. 76, 225249.
Warming, S., Costantino, N., Court, D.L., Jenkins, N.A., and Copeland, N.G.
(2005). Simple and highly efcient BAC recombineering using galK selection.
Nucleic Acids Res. 33, e36.
Kajimura, S., Seale, P., Kubota, K., Lunsford, E., Frangioni, J.V., Gygi, S.P.,
and Spiegelman, B.M. (2009). Initiation of myoblast to brown fat switch by a
PRDM16-C/EBP-beta transcriptional complex. Nature 460, 11541158.
Wu, J., Bostrom, P., Sparks, L.M., Ye, L., Choi, J.H., Giang, A.H., Khandekar,
M., Virtanen, K.A., Nuutila, P., Schaart, G., et al. (2012). Beige adipocytes are a
distinct type of thermogenic fat cell in mouse and human. Cell 150, 366376.
Kim, G.W., Lin, J.E., Blomain, E.S., and Waldman, S.A. (2014). Antiobesity
pharmacotherapy: new drugs and emerging targets. Clin. Pharmacol. Ther.
95, 5366.
Xue, B., Rim, J.S., Hogan, J.C., Coulter, A.A., Koza, R.A., and Kozak, L.P.
(2007). Genetic variability affects the development of brown adipocytes in
white fat but not in interscapular brown fat. J. Lipid Res. 48, 4151.
Cell Reports 9, 15841593, December 11, 2014 2014 The Authors 1593
Cell Reports
Report
TRIM28 Represses Transcription of Endogenous
Retroviruses in Neural Progenitor Cells
Liana Fasching,1 Adamandia Kapopoulou,2 Rohit Sachdeva,1 Rebecca Petri,1 Marie E. Jonsson,1 Christian Manne,1
Priscilla Turelli,2 Patric Jern,3 Florence Cammas,4 Didier Trono,2 and Johan Jakobsson1,*
1Laboratory
of Molecular Neurogenetics, Department of Experimental Medical Science, Wallenberg Neuroscience Center and Lund
Stem Cell Center, Lund University, 221 84 Lund, Sweden
2School of Life Sciences, Ecole Polytechnique Fe
derale de Lausanne (EPFL), 1015 Lausanne, Switzerland
3Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, 751 23 Uppsala, Sweden
4IRCM, Institut de Recherche en Cance
rologie de Montpellier, INSERM, U896, Universite Montpellier; Institut Regional du Cancer Montpellier,
Montpellier 34298, France
*Correspondence: [email protected]
http://dx.doi.org/10.1016/j.celrep.2014.12.004
This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/3.0/).
SUMMARY
TRIM28 is a corepressor that mediates transcriptional silencing by establishing local heterochromatin. Here, we show that deletion of TRIM28 in neural
progenitor cells (NPCs) results in high-level expression of two groups of endogenous retroviruses
(ERVs): IAP1 and MMERVK10C. We nd that NPCs
use TRIM28-mediated histone modications to
dynamically regulate transcription and silencing of
ERVs, which is in contrast to other somatic cell types
using DNA methylation. We also show that derepression of ERVs inuences transcriptional dynamics in
NPCs through the activation of nearby genes and
the expression of long noncoding RNAs. These ndings demonstrate a unique dynamic transcriptional
regulation of ERVs in NPCs. Our results warrant
future studies on the role of ERVs in the healthy and
diseased brain.
INTRODUCTION
The mammalian brain is an extremely complex organ harboring
more than a thousand different types of neurons that serve a
wide variety of functions. How this complexity is achieved
remains largely unknown. However, epigenetic mechanisms
such as DNA methylation, histone modication, and noncoding
RNAs are thought to be important in establishing a high diversity
of gene expression from the same template, leading to a spatial
pattern of transcription. How distinct transcriptional programs
are established in different neuronal populations remains poorly
understood, but one interesting recently proposed hypothesis
suggests transposable elements (TEs) to be involved in this process (Muotri et al., 2007; Reilly et al., 2013). TEs are repetitive
mobile genetic elements that were originally considered to be
parasitic DNA without any function, popularly termed junk
DNA. Today, it is becoming increasingly clear that TEs can
act as gene regulatory elements by serving as hubs for chromatin
20 Cell Reports 10, 2028, January 6, 2015 2015 The Authors
fully methylated in DNA extracted from mouse tail (7% unmethylated CpGs, Figure 2F, Fishers exact test one-sided p < 0.05).
Moreover, we found no difference in the level of CpG methylation
between wild-type and Trim28/ NPCs. In summary, these data
suggest that a proportion of the MMERVK10C elements are
spared from undergoing DNA methylation specically in NPCs
during early development.
Increased Expression of IAP1 Results in ERV-Derived
Protein Expression
IAP1 elements, which lose H3K9me3 marks and were also
highly upregulated in Trim28/ NPCs (Figures 3A and 3B), are
Figure
2. Analysis
of
MMERVK10C Provirus
the
Putative
SUPPLEMENTAL INFORMATION
Supplemental Information includes Supplemental Experimental Procedures,
three gures, and three tables and can be found with this article online at
http://dx.doi.org/10.1016/j.celrep.2014.12.004.
AUTHOR CONTRIBUTIONS
L.F., A.K., R.S., R.P., M.E.J., and C.M. designed and performed research and
analyzed data. P.J., P.T., and D.T. designed research and analyzed data.
F.C. contributed reagents. J.J. designed and coordinated the project and
analyzed data. L.F. and J.J. wrote the paper, and all authors reviewed the
manuscript.
ACKNOWLEDGMENTS
We are grateful to A. Bjorklund, S. Quenneville, and all members of the J.J. and
Parmar laboratories for stimulating discussions. We thank U. Jarl, A. Josefsson, C. Isaksson, I. Nilsson, A.-K. Olden, E. Ling, S. Smiljanic, M. Sparrenius,
and E. Tjon for technical assistance. This study was supported by grants from
Swedish Research Council (J.J.), Formas (P.J.), the Swedish Cancer Foundation (J.J.), the Lundqvist, Jeansson, and Crafoord foundations (J.J.), the
Swedish Government Initiative for Strategic Research Areas MultiPark (J.J.),
the French government, CNRS and INSERM (F.C.), and the French Agence
Nationale pour la Recherche (F.C.).
Received: June 10, 2014
Revised: October 28, 2014
Accepted: December 1, 2014
Published: December 24, 2014
REFERENCES
Baillie, J.K., Barnett, M.W., Upton, K.R., Gerhardt, D.J., Richmond, T.A., De
Sapio, F., Brennan, P.M., Rizzu, P., Smith, S., Fell, M., et al. (2011). Somatic
retrotransposition alters the genetic landscape of the human brain. Nature
479, 534537.
Belshaw, R., Dawson, A.L., Woolven-Allen, J., Redding, J., Burt, A., and
Tristem, M. (2005). Genomewide screening reveals high levels of insertional
polymorphism in the human endogenous retrovirus family HERV-K(HML2):
implications for present-day activity. J. Virol. 79, 1250712514.
Bojkowska, K., Aloisio, F., Cassano, M., Kapopoulou, A., Santoni de Sio, F.,
Zangger, N., Offner, S., Cartoni, C., Thomas, C., Quenneville, S., et al.
(2012). Liver-specic ablation of Kruppel-associated box-associated protein
1 in mice leads to male-predominant hepatosteatosis and development of liver
adenoma. Hepatology 56, 12791290.
Cammas, F., Mark, M., Dolle, P., Dierich, A., Chambon, P., and Losson, R.
(2000). Mice lacking the transcriptional corepressor TIF1beta are defective in
early postimplantation development. Development 127, 29552963.
Chuong, E.B., Rumi, M.A., Soares, M.J., and Baker, J.C. (2013).
Endogenous retroviruses function as species-specic enhancer elements in
the placenta. Nat. Genet. 45, 325329.
Conti, L., Pollard, S.M., Gorba, T., Reitano, E., Toselli, M., Biella, G., Sun, Y.,
Sanzone, S., Ying, Q.L., Cattaneo, E., and Smith, A. (2005). Nicheindependent symmetrical self-renewal of a mammalian tissue stem cell.
PLoS Biol. 3, e283.
Coufal, N.G., Garcia-Perez, J.L., Peng, G.E., Yeo, G.W., Mu, Y., Lovci, M.T.,
Morell, M., OShea, K.S., Moran, J.V., and Gage, F.H. (2009). L1 retrotransposition in human neural progenitor cells. Nature 460, 11271131.
ACCESSION NUMBERS
Cowley, M., and Oakey, R.J. (2013). Transposable elements re-wire and
ne-tune the transcriptome. PLoS Genet. 9, e1003234.
The RNA-seq data were deposited in the NCBI Gene Expression Omnibus and
are available under accession number GSE45930.
Douville, R., Liu, J., Rothstein, J., and Nath, A. (2011). Identication of active
loci of a human endogenous retrovirus in neurons of patients with amyotrophic
lateral sclerosis. Ann. Neurol. 69, 141151.
Muotri, A.R., Marchetto, M.C., Coufal, N.G., Oefner, R., Yeo, G., Nakashima,
K., and Gage, F.H. (2010). L1 retrotransposition in neurons is modulated by
MeCP2. Nature 468, 443446.
Perrat, P.N., DasGupta, S., Wang, J., Theurkauf, W., Weng, Z., Rosbash, M.,
and Waddell, S. (2013). Transposition-driven genomic heterogeneity in the
Drosophila brain. Science 340, 9195.
Qin, C., Wang, Z., Shang, J., Bekkari, K., Liu, R., Pacchione, S., McNulty, K.A.,
Ng, A., Barnum, J.E., and Storer, R.D. (2010). Intracisternal A particle genes:
Distribution in the mouse genome, active subtypes, and potential roles as species-specic mediators of susceptibility to cancer. Mol. Carcinog. 49, 5467.
Reichmann, J., Crichton, J.H., Madej, M.J., Taggart, M., Gautier, P., GarciaPerez, J.L., Meehan, R.R., and Adams, I.R. (2012). Microarray analysis of
LTR retrotransposon silencing identies Hdac1 as a regulator of retrotransposon expression in mouse embryonic stem cells. PLoS Comput. Biol. 8,
e1002486.
Reilly, M.T., Faulkner, G.J., Dubnau, J., Ponomarev, I., and Gage, F.H. (2013).
The role of transposable elements in health and diseases of the central nervous
system. J. Neurosci. 33, 1757717586.
Ribet, D., Harper, F., Dupressoir, A., Dewannieux, M., Pierron, G., and
Heidmann, T. (2008). An infectious progenitor for the murine IAP retrotransposon: emergence of an intracellular genetic parasite from an ancient retrovirus.
Genome Res. 18, 597609.
Rowe, H.M., and Trono, D. (2011). Dynamic control of endogenous retroviruses during development. Virology 411, 273287.
Rowe, H.M., Jakobsson, J., Mesnard, D., Rougemont, J., Reynard, S., Aktas,
T., Maillard, P.V., Layard-Liesching, H., Verp, S., Marquis, J., et al. (2010).
KAP1 controls endogenous retroviruses in embryonic stem cells. Nature
463, 237240.
Jurka, J., Kapitonov, V.V., Pavlicek, A., Klonowski, P., Kohany, O., and
Walichiewicz, J. (2005). Repbase Update, a database of eukaryotic repetitive
elements. Cytogenet. Genome Res. 110, 462467.
Rowe, H.M., Kapopoulou, A., Corsinotti, A., Fasching, L., Macfarlan, T.S.,
Tarabay, Y., Viville, S., Jakobsson, J., Pfaff, S.L., and Trono, D. (2013).
TRIM28 repression of retrotransposon-based enhancers is necessary to preserve transcriptional dynamics in embryonic stem cells. Genome Res. 23,
452461.
Karlsson, H., Bachmann, S., Schroder, J., McArthur, J., Torrey, E.F., and
Yolken, R.H. (2001). Retroviral RNA identied in the cerebrospinal uids and
brains of individuals with schizophrenia. Proc. Natl. Acad. Sci. USA 98,
46344639.
Sachdeva, R., Jonsson, M.E., Nelander, J., Kirkeby, A., Guibentif, C., Gentner,
B., Naldini, L., Bjorklund, A., Parmar, M., and Jakobsson, J. (2010). Tracking
differentiating neural progenitors in pluripotent cultures using microRNAregulated lentiviral vectors. Proc. Natl. Acad. Sci. USA 107, 1160211607.
Kumaki, Y., Oda, M., and Okano, M. (2008). QUMA: quantication tool for
methylation analysis. Nucleic Acids Res. 36 (Web Server issue), W170-5.
Santoni de Sio, F.R., Barde, I., Offner, S., Kapopoulou, A., Corsinotti, A.,
Bojkowska, K., Genolet, R., Thomas, J.H., Luescher, I.F., Pinschewer, D.,
et al. (2012a). KAP1 regulates gene networks controlling T-cell development
and responsiveness. FASEB J. 26, 45614575.
Jern, P., and Cofn, J.M. (2008). Effects of retroviruses on host genome function. Annu. Rev. Genet. 42, 709732.
Kunarso, G., Chia, N.Y., Jeyakani, J., Hwang, C., Lu, X., Chan, Y.S., Ng, H.H.,
and Bourque, G. (2010). Transposable elements have rewired the core regulatory network of human embryonic stem cells. Nat. Genet. 42, 631634.
Langmead, B., Trapnell, C., Pop, M., and Salzberg, S.L. (2009). Ultrafast and
memory-efcient alignment of short DNA sequences to the human genome.
Genome Biol. 10, R25.
Li, J., Akagi, K., Hu, Y., Trivett, A.L., Hlynialuk, C.J., Swing, D.A., Volfovsky, N.,
Morgan, T.C., Golubeva, Y., Stephens, R.M., et al. (2012). Mouse endogenous
retroviruses can trigger premature transcriptional termination at a distance.
Genome Res. 22, 870884.
Li, W., Prazak, L., Chatterjee, N., Gruninger, S., Krug, L., Theodorou, D., and
Dubnau, J. (2013). Activation of transposable elements during aging and
neuronal decline in Drosophila. Nat. Neurosci. 16, 529531.
Santoni de Sio, F.R., Massacand, J., Barde, I., Offner, S., Corsinotti, A., Kapopoulou, A., Bojkowska, K., Dagklis, A., Fernandez, M., Ghia, P., et al. (2012b).
KAP1 regulates gene networks controlling mouse B-lymphoid cell differentiation and function. Blood 119, 46754685.
Schultz, D.C., Ayyanathan, K., Negorev, D., Maul, G.G., and Rauscher, F.J.,
3rd. (2002). SETDB1: a novel KAP-1-associated histone H3, lysine 9-specic
methyltransferase that contributes to HP1-mediated silencing of euchromatic
genes by KRAB zinc-nger proteins. Genes Dev. 16, 919932.
Sperber, G.O., Airola, T., Jern, P., and Blomberg, J. (2007). Automated recognition of retroviral sequences in genomic dataRetroTector. Nucleic Acids
Res. 35, 49644976.
Matsui, T., Leung, D., Miyashita, H., Maksakova, I.A., Miyachi, H., Kimura, H.,
Tachibana, M., Lorincz, M.C., and Shinkai, Y. (2010). Proviral silencing in
embryonic stem cells requires the histone methyltransferase ESET. Nature
464, 927931.
Sripathy, S.P., Stevens, J., and Schultz, D.C. (2006). The KAP1 corepressor
functions to coordinate the assembly of de novo HP1-demarcated microenvironments of heterochromatin required for KRAB zinc nger protein-mediated
transcriptional repression. Mol. Cell. Biol. 26, 86238638.
Muotri, A.R., Chu, V.T., Marchetto, M.C., Deng, W., Moran, J.V., and Gage,
F.H. (2005). Somatic mosaicism in neuronal precursor cells mediated by L1
retrotransposition. Nature 435, 903910.
Muotri, A.R., Marchetto, M.C., Coufal, N.G., and Gage, F.H. (2007). The necessary junk: new functions for transposable elements. Hum. Mol. Genet. 16 Spec
No. 2, R159R167.
Whitelaw, N.C., Chong, S., Morgan, D.K., Nestor, C., Bruxner, T.J., Ashe, A.,
Lambley, E., Meehan, R., and Whitelaw, E. (2010). Reduced levels of two
modiers of epigenetic gene silencing, Dnmt3a and Trim28, cause increased
phenotypic noise. Genome Biol. 11, R111.
Wiznerowicz, M., Jakobsson, J., Szulc, J., Liao, S., Quazzola, A., Beermann,
F., Aebischer, P., and Trono, D. (2007). The Kruppel-associated box repressor
domain can trigger de novo promoter methylation during mouse early embryogenesis. J. Biol. Chem. 282, 3453534541.
Wolf, D., and Goff, S.P. (2009). Embryonic stem cells use ZFP809 to silence
retroviral DNAs. Nature 458, 12011204.
Cell Reports
Article
Comparative Hi-C Reveals that CTCF Underlies
Evolution of Chromosomal Domain Architecture
Matteo Vietri Rudan,1 Christopher Barrington,1 Stephen Henderson,1 Christina Ernst,2 Duncan T. Odom,2 Amos Tanay,3
and Suzana Hadjur1,*
1Research
Department of Cancer Biology, Cancer Institute, University College London, 72 Huntley Street, London WC1E 6BT, UK
Research UK, Cambridge Research Institute, Li Ka Shing Centre, Robinson Way, Cambridge CB2 0RE, UK
3Department of Computer Science and Applied Mathematics, Department of Biological Regulation, Weizmann Institute, Rehovot 76100,
Israel
*Correspondence: [email protected]
http://dx.doi.org/10.1016/j.celrep.2015.02.004
This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/3.0/).
2Cancer
SUMMARY
genomes suggests its evolutionary dynamics are context dependent, and conservation can be interrupted by mobile element
activity (Schmidt et al., 2012). Despite these observations, a
link between the evolutionary dynamics of CTCF binding and
the evolution of chromosomal domain organization is yet to be
explored.
Studies that have tracked the evolution of different transcription factor (TF) binding patterns have shown that sequence
evolution alone is incapable of fully explaining the evolutionary
dynamics of TF binding landscapes (Dermitzakis and Clark,
2001; Birney et al., 2007; Borneman et al., 2007; Schmidt
et al., 2010). TF binding landscapes and large-scale chromosomal organization may function cooperatively to drive the evolution of genome regulation. These observations highlight the
importance of multi-species comparative chromosomal structure analysis and its integration with insulator binding proles
across evolution. If the binding patterns of trans-factors such
as CTCF are indeed strong drivers of domain organization,
then their evolutionary dynamics should drive evolutionary conservation and divergence of chromosome domains.
With this in mind, we performed comparative Hi-C in noncycling primary liver cells and analyzed the data together with
CTCF binding proles from the same species and tissue. Analysis of four mammalian Hi-C maps allowed us to explore how
the evolution of CTCF binding proles correlates, and in some
cases likely drives, the evolution of chromosomal topologies.
We nd that the large-scale chromosomal domain structure is
highly conserved between species, in a way that is correlated
with the conservation of both the CTCF binding site and the
orientation of its motif, resulting in directional long-range interactions that demarcate conserved domains. On the other hand, internal domain structure is observed to be more dynamic, and we
discover remarkable correlation between evolutionary dynamics
of CTCF sites and divergence of local insulation structure. Since
the evolution of CTCF binding proles is strongly driven at the
nucleotide level within cis elements, our data suggest that internal domain structure can be modulated exibly through local
sequence evolution. Conversely, we show that interruption of
large-scale domain structure is rare, and we suggest that instead
of local sequence divergence, evolutionary manipulation of
global chromosomal topologies is driven by processes involving
duplications or rearrangements such as inversions, insertions/
deletions, and translocations. We demonstrate this by charting
cases of evolutionary domain shufing in mouse and dog.
RESULTS
Sequence-Driven Evolution of CTCF Binding Proles
CTCF binding is strongly correlated with the topological architecture of mammalian chromosomes and participates in longrange chromatin loops, thereby underlying global contact insulation. We analyzed mouse (Mus musculus [Mmus]), dog (Canis
familiaris [Cfam]), and macaque (Macaca mulatta [Mmul]) CTCF
chromatin immunoprecipitation sequencing (ChIP-seq) proles
from primary liver cells (Schmidt et al., 2012), aiming to dene
how conservation and divergence of the insulator binding landscape co-evolve with chromosomal topology. Pairwise CTCF
ChIP-seq analysis identied conserved or divergent CTCF bind1298 Cell Reports 10, 12971309, March 3, 2015 2015 The Authors
250 Kb
Mmus
CTCF
Cfam
CTCF
Mmul
CTCF
4
6
2
5
4
0
3
-2
Mmus
1
0
0
-4
-4
-6
-2
Mmus
Mmus CTCF
Mmus+/Cfam+
Mmus+/Cfam-
Mmus+/Cfam+
Mmus+/Cfam-
Normalized count
2.5
Contact enrichment
(log2 scale)
52.5 16
14
11
52 -
7
2
1.5
51.5 0
20
40
60
80
100
Band
51 -
Cell Reports 10, 12971309, March 3, 2015 2015 The Authors 1299
Ocun
Mmus
Crossover 160kb
Crossover 80kb
Orth. genes
Clusters
96 42 98 chr 12
chr 10
40 -
38 -
100 -
36 104 -
Cfam
Mmul
Crossover 160kb
Crossover 80kb
Orth. genes
16
14
11
7
2
Clusters
156 -
64 -
66 chr 12
chr 4
154 -
152 -
68 -
150 -
70 -
148 -
Mmus
Prk
Ocun
Hs3st5
Hdac2
Mmul
Hs3st5
Hdac2
Prk
102
150
36.5
102.5
149.5
36
103
35.5
Hs3st5
Hdac2
Prk
37
Cfam
Hs3st5
Hdac2
Prk
70
70.5
149
103.5
148.5
104
148
71
35
34.5
71.5
10
10
12
Mmus interactions
14
12
Cfam interactions
12
Mmul interactions
Ocun interactions
10
10
12
Mmus interactions
14
12
10
10
12
14
Mmus interactions
(legend on next page)
1300 Cell Reports 10, 12971309, March 3, 2015 2015 The Authors
across all species (Figure 2A). The maps also revealed evidence
of intra-domain differences between species (Figure 2B). We
quantied the extent of structural conservation genome-wide
using a computational approach that allowed us to comprehensively describe domain structure at multiple scales. This pairwise
approach revealed extensive genome-wide interspecies conservation of chromosome structure (Figures 2C and S3). A systematic analysis of paired domains in mouse and dog revealed that
conserved domains are smaller in size compared to other domains and are classied as both active and passive clusters (Figure S4). Together, these data facilitated extensive analysis of the
evolution of chromosomal topologies within regions that did not
go through substantial genome rearrangement, allowing examination of the evolution of both large-scale domain borders and
the insulation structure within domains.
Divergent CTCF Binding Drives Local Structural Change
within Domains
Hi-C maps from liver cells of different species allowed us to ask
how the evolutionary dynamics of CTCF correlate with conservation or divergence of domain structure. Analysis of specic loci
showed that conserved CTCF sites were typically located at
the borders of large-scale chromosomal domains that were
themselves conserved between mouse and dog (Figure 3A). To
test these observations globally, we computed the contact insulation proles from either the mouse or dog Hi-C maps around
conserved CTCF sites, showing that these sites indeed globally
served as conserved insulation points (Figure 3B). Similar results
were derived using a comparison of mouse and macaque (Figure S5). We also observed that conserved CTCF sites were
strongly enriched for Rad21 in mouse (79% of conserved sites
compared to 51% mouse-divergent sites co-localize with
Rad21) and that CTCF/cohesin co-occupied sites exhibited
strong contact insulation in all three species (data not shown).
In contrast to these highly stable sites, our data showed that
divergent CTCF sites were located primarily within domains
and exhibited local contact insulation. Comparative analysis of
contact insulation at divergent CTCF sites revealed that indeed
these sites correlated with divergent contact insulation proles.
For example, dog-divergent CTCF sites (Mmus/Cfam+) exhibited local contact insulation specically in the dog genome,
whereas these same sites exhibited background levels of contact insulation when examined in the mouse Hi-C data (Figures
3C and S5). Importantly, the change in insulation following
CTCF binding site evolution was stronger at the local (20-kb)
scale, but was not signicant at the higher (80-kb) scale (Fig-
Cell Reports 10, 12971309, March 3, 2015 2015 The Authors 1301
Cfam Hi-C
Mmus Hi-C
Mmus CTCF
Cfam CTCF
Mmus+/Cfam+
Mmus+/Cfam+
Mmus-/Cfam+
Mmus+/Cfam-
63 90 63.5 89.5 64 -
16
14
11
7
2
89 64.5 -
chr 8
chr 2
Hi-C
Mmus
Mmus+
Cfam+
Cfam
80K
40K
20K
10K
5K
-200
-100
100
200
-200
-100
100
200
CTCF ChIP
C
Mmus+
Cfam-
80K
40K
20K
10K
5K
-200
MmusCfam+
-100
100
200
-200
-100
100
200
80K
40K
20K
10K
5K
-200
-100
100
200
-200
-100
100
Contact enrichment
(log2 scale)
CTCF ChIP
0.3
0.1
0
-0.1
-0.3
-0.5
200
Mmus - Cfam
contact enrichment (80kb)
3
2
1
0
-1
-2
-3
2
1
0
-1
-2
-3
Mmus+
Cfam-
Mmus+
Cfam+
MmusCfam+
Mmus+
Cfam-
1302 Cell Reports 10, 12971309, March 3, 2015 2015 The Authors
of CTCFs long-range interactions. Consistent with this, we proled the genome-wide relative position within chromosomal
domains of Mmus+/Cfam+ conserved CTCF sites grouped according to the orientation of their binding motif as above. We
observed that the conserved CTCF binding sites that are enriched at the edges of conserved domains (Figure 1E) have a
specic orientation of their motif relative to chromosomal domains (Figure 4B). These observations were replicated when
we compared mouse and macaque.
To characterize the contact relationship between evolutionarily stable or exible CTCF sites and to further understand
how they contribute to the evolution of chromosome domain
structure, we performed a high-resolution high-throughput
circular chromosome conformation capture (4C-seq) study. We
designed four 4C-seq viewpoints to a series of neighboring
conserved CTCF binding sites bordering conserved domains in
the mouse and dog as well as to a mouse-specic site. The results showed that each conserved CTCF site engages in very
strong and directional interactions with neighboring conserved
CTCF sites (Figure 4C). Remarkably, the specic interactions
mediated by conserved sites in the mouse genome were themselves precisely conserved in the dog genome (Figure 4D) and
dene the underlying domain structure. In each case, the longrange interaction was anchored by a pair of conserved CTCF
sites whereby one CTCF site had an orientation on the + strand
and the other on the strand and could provide the basis for
the observed directionality of CTCF-mediated interactions.
Moreover, a viewpoint designed to a mouse-divergent site exhibited weak interactions within the mouse domain, analogous
to the local insulation behavior observed in Figure 3B (Figure 4C).
Importantly, the mouse-divergent viewpoint had no prominent
interactions in the dog genome, conrming the specicity of its
interaction network.
Global analysis of Hi-C contacts between pairs of CTCF binding sites stratied according to genomic distances in cis (Sofueva et al., 2013) conrmed the 4C-seq observation systematically (Figure 4E). Consistent with the high-resolution 4C-seq
proles, Hi-C trends showed that conserved CTCF sites strongly
contacted one another within the same domain. Divergent CTCF
sites engaged in signicantly weaker contacts with other divergent sites, even when stratifying thoroughly for genomic distances. Importantly, little to no contact was observed in the
mouse genome between dog-divergent sites. These results
show that evolutionarily stable CTCF sites are engaged in strong
contacts with one another and suggest that in so doing, they
create an interaction network that may support the conservation
Figure 3. Conserved and Divergent CTCF Sites Show Differential Contact Insulation Behavior
(A) Representative Hi-C contact maps for a 2-Mb syntenic region in mouse (left panel) and dog (right panel). Also shown are the CTCF ChIP-seq tracks in each
species as well as the conserved (Mmus+/Cfam+, blue) and divergent (Mmus+/Cfam, gray or Mmus/Cfam+, green) CTCF sites.
(B) Average contact insulation analysis at conserved (Mmus+/Cfam+) CTCF sites in mouse (leftmost panel) and dog (rightmost panel) Hi-C datasets. Arrows
indicate liftOver of sites.
(C) Same as for (B), but for mouse-divergent (Mmus+/Cfam) and dog-divergent (Mmus/Cfam+) sites in both genomes. Divergent sites appear to mediate
weaker, shorter-range insulation that disappears in the species where the site is not bound.
(D) Distribution of the difference in contact insulation between mouse and dog at conserved or divergent CTCF sites at 2060 kb (left panel, 20-kb band) or 80240
kb (right panel, 80-kb band) scales. Divergent sites exhibit a signicant (Kolmogorov-Smirnov test, p < 1 3 1012, marked by an asterisk) shift in their distributions
in the 20-kb band, but not in the 80-kb band, compatible with them mediating insulation at a local scale, but not at higher ranges.
See also Figure S5.
Cell Reports 10, 12971309, March 3, 2015 2015 The Authors 1303
1304 Cell Reports 10, 12971309, March 3, 2015 2015 The Authors
somal domains to the exibility of continuous genomic adaptation. CTCF and cohesin complexes are deeply evolutionarily
conserved, and the data here show that their role in mediating
chromosome topologies and, even more remarkably, the
large-scale building blocks of such topologies are also highly
conserved. Our data suggest that the orientation of the CTCF
motif may underlie the observed directionality of CTCF/cohesin-mediated long-range contacts and provide a rationale by
which specic sites are assembled to dene topological domains. Given that CTCF binding is strongly inuenced by its
consensus sequence, our data suggest that the assembly of
domain structure is hardwired in the genome. This also
has implications for further understanding the nature of the
relationship between CTCF and cohesin, since biochemical
studies have revealed that cohesin subunits interact with
CTCF primarily through its C-terminal tail (Xiao et al., 2011),
placing cohesin on a particular side of the chromosomal
domain.
Interestingly, while we were able to observe cases whereby
local sequence evolution perturbed CTCF binding and disrupted
chromosomal looping, the structures that were affected due to
this insulator divergence were primarily local loops. Cases of
large-scale topological domains that were split or fused due to
insulator divergence were not observed. We hypothesize that
this stability is achieved by a combination of both local purifying
selection on key CTCF binding sites and by buffering of major topological loops by additional factors. Strikingly, the cases of
large-scale domain divergence that we were able to characterize
were all linked with evolutionary genome rearrangements and revealed a mechanism that can reshufe whole domains such that
the rearranged chromosomal modules are aligned with existing
domain borders. It is still, however, formally possible that rearrangements take place between CTCF sites that are mediating
strong interactions.
In addition to the importance of topological domain and insulator conservation described here, the evolutionary dynamics
that couple intra-domain CTCF divergence with changes in the
local domain structure emerge as potentially fundamental for
genome regulation. Loops contained within domains link enhancers (and their bound trans-factors) to target gene promoters. While it is still unclear how such targeting is regulated
and how evolution can manipulate it, based on our data, we hypothesize that exible CTCF binding sites within domains can
inuence looping from the promoter or enhancer as well by
Figure 4. Conserved CTCF Sites Engage in Strong, Directional Interactions with Other Conserved Sites
(A) Average contact insulation analysis in both Mmus and Cfam genomes for conserved (Mmus+/Cfam+) CTCF binding sites grouped according to the orientation
of their binding motif. The consensus motif shown was generated from mouse binding sites.
(B) Genome-wide relative position within chromosomal domains of Mmus+/Cfam+ conserved CTCF sites grouped according to the orientation of their binding
motif as above.
(C and D) (C) 4C-seq analysis and Hi-C maps of a 1.2-Mb syntenic region in the mouse and (D) dog genomes. Shown are 4C-seq viewpoints designed to four
conserved (Mmus+/Cfam+) CTCF binding sites proximal to Hi-C domain borders and one 4C-seq viewpoint located at a mouse-divergent (Mmus+/Cfam) CTCF
site. The symbol above each 4C-seq bait indicates the strand (and orientation) of each viewpoint. Each 4C-seq experiment is represented by the median
normalized 4C-seq coverage in a sliding window of 5 kb (top) and a multi-scale domainogram indicating normalized mean coverage in windows ranging between
2 kb and 50 kb.
(E) Relative intra-domain contact enrichment between CTCF sites (solid lines) as a function of distance when the two sites are <5 kb away from a CTCF binding
site (solid lines) or where only one site is <5 kb away from a CTCF site (dotted lines). Shown are the relative contact enrichments within repressive domains for
conserved (Mmus+/Cfam+, blue), mouse-divergent (Mmus+/Cfam, gray) and dog-divergent (Mmus/Cfam+, green) CTCF sites in mouse Hi-C (left) and dog
Hi-C (right).
Cell Reports 10, 12971309, March 3, 2015 2015 The Authors 1305
Cfam
chr15
Mmus
chr4
108
112
114.2 116Mb
12.6
14.1Mb
Skint gene
family
insertion
B
Cfam
Mmus
Mmus
CTCF
Mmus+
Cfam+
Mmus+
Cfam-
Cfam
CTCF
Mmus+
Cfam+
MmusCfam+
Genes
Genes
chr 15
Contact enrichment
(log2 scale)
2
11 14 16
EXPERIMENTAL PROCEDURES
Liver Homogenization and Fixation
Fresh or frozen liver from mouse, rabbit, macaque, and dog were processed
for Hi-C or 4C-seq libraries. With the exception of mouse, the samples used
for the Hi-C libraries were the same as the material used for CTCF ChIP-seq
(Schmidt et al., 2012). Livers were xed in 10% formalin for 20 min, and
1 g was cut and processed with a Dounce homogenizer (ten strokes with a
CTCF Peak
Reading Primer
Non-reading Primer
chr10:94609250
50 -CCATCTGTTTGAACAAGATC-30
chr10:94623583
chr10:94958324
chr10:94991353
5 -AGTCAGATGGAATGCAGATC-3
0
5 -ATTGCTTTCTCTGGTTGATC-3
50 -CAAGAGAGAGTGGAAACAGG-30
0
50 -CTAGATACAGCAATCAGCCC-30
50 -AGTCACTCCTGCTCCTGTAA-30
50 -AAGCATTGTCCTACGTGATT-30
50 -CCCTTCCCTTCTATGTTTCT-30
5 -GTTTCTGTTGGTTCACGATC-3
chr10:95218005
5 -CTACTCTGGCTTCTATGATC-3
chr15:34606369
50 -GCTCTTGCTCTAAACTGATC-30
50 -TGGACCTCACCTCTCCTA-30
50 -GTCGCATCACTTACTGGG-30
chr15:34596229
5 -TGAGGTCCAGCAGAGATC-3
chr15:34269944
50 -CTCCACTGAGCATTAAGATC-30
50 -GCGGGATAGTTCTTTTCTCT-30
chr15:34244336
50 -CTTATGTGCTCCTCCAGATC-30
50 -AATCATATGCCTCCTCCTCT-30
chr15:33989305
50 -AAAGTAATCCCACCCAGATC-30
50 -CTGAAGGAAACAACAATGTCA-30
loose pestle followed by ten strokes with the tight pestle). After ltration
through a 70-mm nylon cell strainer, the sample was washed twice with PBS,
spinning down at 852 rcf for 5 min at 4 C to collect the cells between washes.
15 3 107 liver cells were then xed for a second time in xation buffer (1%
formaldehyde, 750 mg/ml BSA in DMEM/Hams F12 [Invitrogen]) for 10
30 min at room temperature. The xation reaction was quenched using
0.125 M glycine for 5 min at room temperature. Samples were washed twice
with 10 ml PBS, pelleted into 1 3 107 cells aliquots, and stored at 80 C.
Mouse Hi-C libraries were prepared from fresh liver samples of biological replicates (9-week-old C57/BL6 mouse and the pooled livers from 2- to 4-weekold outbred mice. The libraries for the other three organisms were technical
replicates.
Propidium Iodide Staining of Hepatocytes
Formaldehyde-xed liver cells were lysed on ice in a hypotonic buffer (10 mM
Tris-HCl [pH 8], 10 mM NaCl, 0.2% Igepal CA-640, EDTA-free protease inhibitors) for 30 min. Nuclei were stained with a propidium iodide (PI) staining buffer
(100 mg/ml PI, 50 mg/ml RNase A, 0.05% Triton X-100) for 60 min on ice. Samples were analyzed on a MoFlo cell sorter (Beckman Coulter).
High-Throughput Mapping of Chromatin Interactions via Hi-C
The Hi-C method previously used (Sofueva et al., 2013) was modied to
accommodate primary liver samples. Hepatocytes were lysed in Hi-C lysis
buffer (10 mM Tris-HCl [pH 8], 10 mM NaCl, 0.2% Igepal CA-640, EDTA-free
protease inhibitors) for 30 min on ice. The sample was transferred to Protein
LoBind tubes (Eppendorf) and the nuclei were permeabilized by incubation
with 0.1%0.6% SDS for 1 hr at 37 C with 800 rpm shaking. The reaction
was quenched with 0.67%4% Triton X-100, 1 hr at 37 C, 800 rpm shaking.
Nuclei were digested in 500 ml 1X NEBuffer 2 with 1500 U HindIII (New England
Biolabs) and monitored for maximal digestion of the chromatin template, thus
digestion times ranged from 24-72 hr. All other parts of the Hi-C protocol,
including library preparation were performed as previously described. 75 bp
paired-end sequencing was performed for each library according to manufacturers conditions using the Illumina Hi-seq platform.
Hi-C Interaction Matrix Generation and Domain Calling
Sequencing reads were aligned to the mouse (mm10), rabbit (oryCun2), macaque (rheMac2), and dog (canFam3) genome assemblies using Bowtie
0.12.8 (Langmead et al., 2009). The parameters used for the alignment allowed
a maximum of three mismatches and strictly one alignment per read. Processing of the aligned reads and normalization of the interaction matrices were
performed as previously described (Yaffe and Tanay, 2011; Sofueva et al.,
2013). The pipeline produced normalized matrices of interactions binning the
genome at different resolutions. Interaction matrices for each library were
generated displaying seven different resolutions simultaneously (12,500,
25,000, 50,000, 100,000, 250,000, 500,000, and 1,000,000 bp). Domains
were identied and clustered as described (Sexton et al., 2012) with the modication that scaling factors were inferred using fends 100400 kb apart, to
account for the lower resolution of the mouse map compared to the Drosophila
map. Domain borders were called using the 95% percentile of the scaling
track. A domain-level map was partitioned into two clusters, and clusters
were assigned as passive/active according to Lamin B mouse embryonic
broblast (MEF) data, as before. For the rabbit, macaque, and dog genome,
the Lamin B MEF track for mouse was lifted over to the corresponding genome
to label domain clusters. Domain calls in mouse and dog are available in
Table S1.
ChIP-Seq Analysis
We used previously published ChIP-seq data for CTCF from mouse, macaque,
and dog livers (Schmidt et al., 2012) and for Rad21 for mouse liver (Faure et al.,
2012). Rad21 ChIP-seq data for macaque and dog was prepared as for CTCF.
Mouse, macaque, and dog ChIP-seq reads were mapped using bowtie. Alignment was followed by extension of sequenced tags to 300-bp fragments and
pileup into 50-bp bins. We normalized ChIP-seq coverage by computing the
distribution of pile-up coverage on 50-bp bins and transforming each coverage
value v into log10 (1-quantile(v)). To dene binding sites, we used a simple
threshold on the sum of values from two biological replicates for each CTCF
dataset and for the macaque Rad21 data. Rad21 ChIP data from mouse and
dog were done in single, and the data were thresholded. Thresholds used
were as follows: mouse CTCF = 2.2, macaque CTCF = 2.4, dog CTCF = 2.2,
mouse Rad21 = 2.3, macaque Rad21 = 2.5, dog Rad21 = 3. Different thresholds did not change the results. Binding site width was standardized at 200 bp,
and the ChIP-seq intensity for each site was calculated as the maximum value
across the 200 bp. The relative distribution of CTCF within topological domains
(Figures 1E and S5) was calculated as the distance of each CTCF site from the
center of its domain. Half the size of the domain was added to convert it to a
measure of distance from the edge of the domain, and this number was then
divided by the size of the domain.
Interspecies Comparison of CTCF Sites
Macaque and dog CTCF ChIP-seq libraries were converted to mouse genome
coordinates using the liftOver tool from UCSC. To reduce the chance of inaccurate liftOver, a number of lters were implemented: sites within low-mappability regions, repeats, or windows of 100 kb with insufcient synteny were
excluded. To estimate mappability, each genome was broken into 50-bp reads
and the whole-genome sequence was split into articial reads and then mapped back to the genome. For each 50-bp bin, the mappability score was then
dened to be the portion of articial reads mapped uniquely to that bin. To estimate the level of synteny in the 100 kb around a CTCF site, the mappability
tracks for macaque and dog were converted to the mouse genome using liftOver and all bins for which liftOver was not possible were converted to zeroes.
The converted tracks were subsequently smoothed over 100 kb, and CTCF
Cell Reports 10, 12971309, March 3, 2015 2015 The Authors 1307
sites falling in regions below the top quartile of such smoothed tracks were
excluded from all subsequent analysis. Divergent CTCF sites in mouse and
dog are available in Table S1.
CTCF Binding Energy Function
A CTCF DNA-binding energy function from the Cortex CTCF binding sites
(ENCODE Cortex CTCF mouse, GSM769019; Shen et al., 2012) was used to
prole all genomes for their similarity to the CTCF consensus motif. The
consensus motif is very highly conserved across all species (Schmidt et al.,
2012). Given a set of genomic sites, we compute for each site the maximal energy value within a 200-bp window centered on the point.
Motif Orientation Analysis
Orientation of the motifs underneath conserved CTCF peaks was obtained using MEME (http://meme.nbcr.net/meme/), (Bailey and Elkan, 1994) with the
parameters -revcomp -dna -nmotifs 1 -w 20 -mod zoops -maxsize 100,000.
Crossover Analysis
Crossover analysis was performed as described previously (Sofueva et al.,
2013). The bands used were 57.5, 7.511.25, 1015, 1522.5, 2030, 30
45, 4060, 6090, and 80120 kb.
Distal Contact Analysis
To calculate the average interaction proles for a group of genomic landmarks,
HindIII fragment ends were grouped into classes by associating each end with
a genomic element located within 5 kb and then grouping all fragment ends
associated with an element of the same class. For the mouse, macaque,
and dog genomes, three classes of CTCF sites (conserved, divergent present,
divergent absent) and TSS sites were dened. These classes were further
divided to sites within active or passive Hi-C domains. The remaining fragment
end (not classied given other landmarks) was dened as the background.
4C-Seq
Preparation of 4C-seq samples, libraries, sequencing analysis, and normalization were all performed as previously described (Sofueva et al., 2013). Primer
sequences were chosen to viewpoint sites that were as close as possible to
CTCF ChIP-seq peaks (Table 1). Mouse primers were designed according to
the genome-wide 4C-seq primer database from (van de Werken et al.,
2012). For dog primers, a similar database was generated for the regions of
interest.
ACCESSION NUMBERS
The data analyzed in this study have been deposited in the GEO database with
the accession number GSE65126.
SUPPLEMENTAL INFORMATION
Supplemental Information includes eight gures and one table and can be
found with this article online at http://dx.doi.org/10.1016/j.celrep.2015.02.004.
AUTHOR CONTRIBUTIONS
M.V.R., D.T.O., and S.H. initiated the project. M.V.R. performed the Hi-C and
4C-Seq experiments, including library preparations. D.T.O. and C.E. provided
the liver samples for all species (except mouse) and sequenced the Hi-C libraries. M.V.R., C.B., and A.T. processed and statistically analyzed the data.
M.V.R., A.T., and S.H. wrote the manuscript, with contributions from all
authors.
ACKNOWLEDGMENTS
The authors wish to thank Sevil Sofueva for help with Hi-C and 4C-seq library
preparations; Wen-Ching Chan for analysis; Christine Feig for Rad21 ChIP-seq
data in macaque and dog; Bianca Schmidt and the Cancer Research UK CI
Genomic Core facility for technical assistance with Hi-C sequencing and Pe-
1308 Cell Reports 10, 12971309, March 3, 2015 2015 The Authors
dro Olivares for advice on analysis and data manipulation. We would also
like to acknowledge all members of the Hadjur group for discussions. This
work was supported by the Medical Research Council UK (G0900491/1 and
G1001649) (S.H.), the EPIGENESYS EU NoE (S.H. and A.T.), and Cancer
Research UK (studentship to M.V.R.). D.T.O. is supported by Cancer
Research UK.
Received: September 27, 2014
Revised: December 9, 2014
Accepted: January 29, 2015
Published: February 26, 2015
REFERENCES
Bailey, T.L., and Elkan, C. (1994). Fitting a mixture model by expectation maximization to discover motifs in bipolymers. Proc. Int. Conf. Intell. Syst. Mol. Biol.
2, 2836.
Bell, A.C., West, A.G., and Felsenfeld, G. (1999). The protein CTCF is required
for the enhancer blocking activity of vertebrate insulators. Cell 98, 387396.
Birney, E., Stamatoyannopoulos, J.A., Dutta, A., Guigo, R., Gingeras, T.R.,
Margulies, E.H., Weng, Z., Snyder, M., Dermitzakis, E.T., Stamatoyannopoulos, J.A., et al. (2007). Identication and analysis of functional elements in
1% of the human genome by the ENCODE pilot project. Nature 447, 799
816. http://dx.doi.org/10.1038/nature05874.
Borneman, A.R., Gianoulis, T.A., Zhang, Z.D., Yu, H., Rozowsky, J., Seringhaus, M.R., Wang, L.Y., Gerstein, M., and Snyder, M. (2007). Divergence of
transcription factor binding sites across related yeast species. Science 317,
815819. http://dx.doi.org/10.1126/science.1140748.
Boyden, L.M., Lewis, J.M., Barbee, S.D., Bas, A., Girardi, M., Hayday, A.C., Tigelaar, R.E., and Lifton, R.P. (2008). Skint1, the prototype of a newly identied
immunoglobulin superfamily gene cluster, positively selects epidermal gammadelta T cells. Nat. Genet. 40, 656662. http://dx.doi.org/10.1038/ng.108.
Dermitzakis, E.T., and Clark, A.G. (2001). Differential selection after duplication
in mammalian developmental genes. Mol. Biol. Evol. 18, 557562.
Dixon, J.R., Selvaraj, S., Yue, F., Kim, A., Li, Y., Shen, Y., Hu, M., Liu, J.S., and
Ren, B. (2012). Topological domains in Mamm. Genomes identied by analysis
of chromatin interactions. Nature 485, 376380. http://dx.doi.org/10.1038/nature11082.
Dong, X., Han, S., Zylka, M.J., Simon, M.I., and Anderson, D.J. (2001). A
diverse family of GPCRs expressed in specic subsets of nociceptive sensory
neurons. Cell 106, 619632.
Dowen, J.M., Fan, Z.P., Hnisz, D., Ren, G., Abraham, B.J., Zhang, L.N., Weintraub, A.S., Schuijers, J., Lee, T.I., Zhao, K., and Young, R.A. (2014). Control of
cell identity genes occurs in insulated neighborhoods in mammalian chromosomes. Cell 159, 374387. http://dx.doi.org/10.1016/j.cell.2014.09.030.
Faure, A.J., Schmidt, D., Watt, S., Schwalie, P.C., Wilson, M.D., Xu, H., Ramsay, R.G., Odom, D.T., and Flicek, P. (2012). Cohesin regulates tissue-specic
expression by stabilizing highly occupied cis-regulatory modules. Genome
Res. 22, 21632175. http://dx.doi.org/10.1101/gr.136507.111.
Filippova, G.N., Fagerlie, S., Klenova, E.M., Myers, C., Dehner, Y., Goodwin,
G., Neiman, P.E., Collins, S.J., and Lobanenkov, V.V. (1996). An exceptionally
conserved transcriptional repressor, CTCF, employs different combinations of
zinc ngers to bind diverged promoter sequences of avian and mammalian cmyc oncogenes. Mol. Cell. Biol. 16, 28022813.
Guacci, V., Koshland, D., and Strunnikov, A. (1997). A direct link between sister
chromatid cohesion and chromosome condensation revealed through the
analysis of MCD1 in S. cerevisiae. Cell 91, 4757.
Hadjur, S., Williams, L.M., Ryan, N.K., Cobb, B.S., Sexton, T., Fraser, P.,
Fisher, A.G., and Merkenschlager, M. (2009). Cohesins form chromosomal
cis-interactions at the developmentally regulated IFNG locus. Nature 460,
410413. http://dx.doi.org/10.1038/nature08079.
Hou, C., Li, L., Qin, Z.S., and Corces, V.G. (2012). Gene density, transcription,
and insulators contribute to the partition of the Drosophila genome into
Kim, T.H., Abdullaev, Z.K., Smith, A.D., Ching, K.A., Loukinov, D.I., Green,
R.D., Zhang, M.Q., Lobanenkov, V.V., and Ren, B. (2007). Analysis of the vertebrate insulator protein CTCF-binding sites in the human genome. Cell 128,
12311245. http://dx.doi.org/10.1016/j.cell.2006.12.048.
Seitan, V.C., Hao, B., Tachibana-Konwalski, K., Lavagnolli, T., Mira-Bontenbal, H., Brown, K.E., Teng, G., Carroll, T., Terry, A., Horan, K., et al. (2011). A
role for cohesin in T-cell-receptor rearrangement and thymocyte differentiation. Nature 476, 467471. http://dx.doi.org/10.1038/nature10312.
Langmead, B., Trapnell, C., Pop, M., and Salzberg, S.L. (2009). Ultrafast and
memory-efcient alignment of short DNA sequences to the human genome.
Genome Biol. 10, R25. http://dx.doi.org/10.1186/gb-2009-10-3-r25.
Michaelis, C., Ciosk, R., and Nasmyth, K. (1997). Cohesins: chromosomal proteins that prevent premature separation of sister chromatids. Cell 91, 3545.
Seitan, V.C., Faure, A.J., Zhan, Y., McCord, R.P., Lajoie, B.R., Ing-Simmons,
E., Lenhard, B., Giorgetti, L., Heard, E., Fisher, A.G., et al. (2013). Cohesinbased chromatin interactions enable regulated gene expression within preexisting architectural compartments. Genome Res. 23, 20662077. http://dx.doi.
org/10.1101/gr.161620.113.
Mishiro, T., Ishihara, K., Hino, S., Tsutsumi, S., Aburatani, H., Shirahige, K., Kinoshita, Y., and Nakao, M. (2009). Architectural roles of multiple chromatin insulators at the human apolipoprotein gene cluster. EMBO J. 28, 12341245.
http://dx.doi.org/10.1038/emboj.2009.81.
Sexton, T., Yaffe, E., Kenigsberg, E., Bantignies, F., Leblanc, B., Hoichman,
M., Parrinello, H., Tanay, A., and Cavalli, G. (2012). Three-dimensional folding
and functional organization principles of the Drosophila genome. Cell 148,
458472. http://dx.doi.org/10.1016/j.cell.2012.01.010.
Nativio, R., Wendt, K.S., Ito, Y., Huddleston, J.E., Uribe-Lewis, S., Woodne,
K., Krueger, C., Reik, W., Peters, J.-M., and Murrell, A. (2009). Cohesin is
required for higher-order chromatin conformation at the imprinted IGF2-H19
locus. PLoS Genet. 5, e1000739. http://dx.doi.org/10.1371/journal.pgen.
1000739.
Shen, Y., Yue, F., McCleary, D.F., Ye, Z., Edsall, L., Kuan, S., Wagner, U.,
Dixon, J., Lee, L., Lobanenkov, V.V., and Ren, B. (2012). A map of the cis-regulatory sequences in the mouse genome. Nature 488, 116120. http://dx.doi.
org/10.1038/nature11243.
Nora, E.P., Lajoie, B.R., Schulz, E.G., Giorgetti, L., Okamoto, I., Servant, N.,
Piolot, T., van Berkum, N.L., Meisig, J., Sedat, J., et al. (2012). Spatial partitioning of the regulatory landscape of the X-inactivation centre. Nature 485, 381
385. http://dx.doi.org/10.1038/nature11049.
Parelho, V., Hadjur, S., Spivakov, M., Leleu, M., Sauer, S., Gregson, H.C., Jarmuz, A., Canzonetta, C., Webster, Z., Nesterova, T., et al. (2008). Cohesins
functionally associate with CTCF on mammalian chromosome arms. Cell
132, 422433. http://dx.doi.org/10.1016/j.cell.2008.01.011.
Pauli, A., Althoff, F., Oliveira, R.A., Heidmann, S., Schuldiner, O., Lehner, C.F.,
Dickson, B.J., and Nasmyth, K. (2008). Cell-type-specic TEV protease cleavage reveals cohesin functions in Drosophila neurons. Dev. Cell 14, 239251.
http://dx.doi.org/10.1016/j.devcel.2007.12.009.
Phillips-Cremins, J.E., Sauria, M.E.G., Sanyal, A., Gerasimova, T.I., Lajoie,
B.R., Bell, J.S.K., Ong, C.-T., Hookway, T.A., Guo, C., Sun, Y., et al. (2013).
Architectural protein subclasses shape 3D organization of genomes during
lineage commitment. Cell 153, 12811295. http://dx.doi.org/10.1016/j.cell.
2013.04.053.
Rollins, R.A., Morcillo, P., and Dorsett, D. (1999). Nipped-B, a Drosophila homologue of chromosomal adherins, participates in activation by remote enhancers in the cut and Ultrabithorax genes. Genetics 152, 577593.
Schmidt, D., Wilson, M.D., Ballester, B., Schwalie, P.C., Brown, G.D.,
Marshall, A., Kutter, C., Watt, S., Martinez-Jimenez, C.P., Mackay, S., et al.
(2010). Five-vertebrate ChIP-seq reveals the evolutionary dynamics of transcription factor binding. Science 328, 10361040. http://dx.doi.org/10.1126/
science.1186176.
Schmidt, D., Schwalie, P.C., Wilson, M.D., Ballester, B., Goncalves, A., Kutter,
C., Brown, G.D., Marshall, A., Flicek, P., and Odom, D.T. (2012). Waves of retrotransposon expansion remodel genome organization and CTCF binding in
Sofueva, S., Yaffe, E., Chan, W.-C., Georgopoulou, D., Vietri Rudan, M., MiraBontenbal, H., Pollard, S.M., Schroth, G.P., Tanay, A., and Hadjur, S. (2013).
Cohesin-mediated interactions organize chromosomal domain architecture.
EMBO J. 32, 31193129. http://dx.doi.org/10.1038/emboj.2013.237.
van de Werken, H.J.G., Landan, G., Holwerda, S.J.B., Hoichman, M., Klous,
P., Chachik, R., Splinter, E., Valdes-Quezada, C., Oz, Y., Bouwman, B.A.M.,
et al. (2012). Robust 4C-seq data analysis to screen for regulatory DNA interactions. Nat. Methods 9, 969972. http://dx.doi.org/10.1038/nmeth.2173.
Wendt, K.S., Yoshida, K., Itoh, T., Bando, M., Koch, B., Schirghuber, E., Tsutsumi, S., Nagae, G., Ishihara, K., Mishiro, T., et al. (2008). Cohesin mediates
transcriptional insulation by CCCTC-binding factor. Nature 451, 796801.
http://dx.doi.org/10.1038/nature06634.
Xiao, T., Wallace, J., and Felsenfeld, G. (2011). Specic sites in the C terminus
of CTCF interact with the SA2 subunit of the cohesin complex and are required
for cohesin-dependent insulation activity. Mol. Cell. Biol. 31, 21742183.
http://dx.doi.org/10.1128/MCB.05093-11.
Yaffe, E., and Tanay, A. (2011). Probabilistic modeling of Hi-C contact maps
eliminates systematic biases to characterize global chromosomal architecture. Nat. Genet. 43, 10591065. http://dx.doi.org/10.1038/ng.947.
Yaffe, E., Farkash-Amar, S., Polten, A., Yakhini, Z., Tanay, A., and Simon, I.
(2010). Comparative analysis of DNA replication timing reveals conserved
large-scale chromosomal architecture. PLoS Genet. 6, e1001011. http://dx.
doi.org/10.1371/journal.pgen.1001011.s016.
Zuin, J., Dixon, J.R., van der Reijden, M.I.J.A., Ye, Z., Kolovos, P., Brouwer,
R.W.W., van de Corput, M.P.C., van de Werken, H.J.G., Knoch, T.A., van
IJcken, W.F.J., et al. (2014). Cohesin and CTCF differentially affect chromatin
architecture and gene expression in human cells. Proc. Natl. Acad. Sci. USA
111, 9961001. http://dx.doi.org/10.1073/pnas.1317788111.
Cell Reports 10, 12971309, March 3, 2015 2015 The Authors 1309
Cell Reports
Article
Heterogeneities in Nanog Expression
Drive Stable Commitment to Pluripotency
in the Mouse Blastocyst
Panagiotis Xenopoulos,1,5 Minjung Kang,1,2,5 Alberto Puliato,3 Stefano Di Talia,4 and Anna-Katerina Hadjantonakis1,*
1Developmental
SUMMARY
1508 Cell Reports 10, 15081520, March 10, 2015 2015 The Authors
1B, S1E, and S1F). Moreover, we observed an increased correlation of reporter activity with NANOG protein for the Nanog:H2BGFP transgene, compared to the targeted Nanog transcriptional
reporter in the heterozygous TNGA ESCs (Chambers et al., 2007)
(Figures S1E and S1F). These data suggest that the BAC transgenic reporters we constructed faithfully marked the pluripotent
state in ESC cultures and could be used to probe Nanog expression dynamics at single-cell resolution.
Next, we generated Nanog:H2B-GFP transgenic mice for
single-cell resolution quantitative visualization of the EPI lineage
in vivo. We used live imaging to analyze the distribution and
validate the reporter in embryos (Figure 1C). We rst observed
reporter activity in embryos at the 816 cell stage (Plusa et al.,
2008). Cells displaying high and low GFP levels were rst
observed at the mid blastocyst stage (70100 cells), as a saltand-pepper distribution of EPI and PrE precursors was established within the ICM (Chazaud et al., 2006; Plusa et al., 2008).
In late blastocysts (>100 cells), where the EPI and PrE lineages
have sorted into distinct layers, cells within the ICM forming
the EPI exhibited signicantly elevated levels of GFP compared
to PrE cells located on the surface of the ICM. Thereafter, GFP
was detected within the EPI, albeit at reduced levels in implanting E4.5 and post-implantation E5.5 embryos consistent with the
downregulation of NANOG observed at periimplantation (Chambers et al., 2003). Notably, the reporter was strongly expressed in
TE cells of early and mid blastocysts (corresponding to 3290
cell stage), possibly resulting from robust NANOG localization in
ICM and TE cells at these stages (Figure S1I) (Dietrich and Hiiragi, 2007; Messerschmidt and Kemler, 2010; Morgani et al.,
2013). However, GFP expression in the TE was signicantly
reduced in late (>100 cells) and in implanting blastocysts
(Figures 1C and S1G). TE localization was also noted for other
Nanog-based reporters analyzed at comparable embryonic
stages (Figure S1G). Collectively, these data lead us to conclude
that the Nanog:H2B-GFP reporter faithfully marks the pluripotent
state both in vitro in ESC cultures and in vivo in the emerging EPI
lineage of the mouse blastocyst.
Derivation of ESCs from Mouse Blastocysts Results in
Highly Variable Pluripotency-Associated Gene
Expression
The EPI lineage of the blastocyst represents the in vivo counterpart to ESC cultures propagated in vitro (Boroviak et al., 2014).
We therefore sought to investigate the prole of Nanog:H2BGFP reporter activity in transgenic embryo-derived ESCs. We
used an ESC derivation protocol in which an ICM outgrowth
emerges in conditions promoting ground state pluripotency in
serum-free medium containing 2i (Czechanski et al., 2014).
Indeed, under these conditions GFP was strongly expressed
by all cells of outgrowths from Nanog:H2B-GFPTg/+ blastocysts
(Figure 1D). However, both reporter and NANOG expression
became heterogeneous when the derived ESCs were propagated under standard serum+LIF conditions, either in the presence or absence of mouse embryonic feeders (MEFs). Even
in the presence of MEFs, GFP-low/NANOG-low cells were
observed within ESC colonies (Figure 1E). Fluorescence-activated cell sorting (FACS) and quantitative immunouorescence
analyses conrmed the presence of a highly heterogeneous
Cell Reports 10, 15081520, March 10, 2015 2015 The Authors 1509
Figure 1. BAC-Based Nanog Transcriptional Reporters Faithfully Mark the Pluripotent State in ESCs and Embryos
(A) Immunouorescence images of Nanog:GFPTg/+ and Nanog:H2B-GFPTg/+ ESCs grown in MEF-free conditions including serum+LIF, 2i+LIF, and serum-LIF for
three passages. Orange arrowheads identify GFP-low and NANOG-low cells within ESC colonies.
(B) Quantitative immunouorescence analysis after nuclear segmentation of Nanog:H2B-GFPTg/+ ESCs. GFP (y axis) and NANOG (x axis) uorescence values
plotted for individual ESCs propagated for ve passages in MEF-free serum+LIF conditions then grown for 4 days in various culture conditions.
(C) Reporter expression in live embryos stained with membrane marker FM4-64. GFP-hi cells, white arrowheads; GFP-low cells, orange arrowheads; TE cells
expressing GFP, yellow arrowheads. Dashed line depicts boundary between epiblast (EPI) and visceral endoderm (VE) layers of an E5.5 embryo. Cell number was
determined by staining with Hoechst.
(D) Schematic of ESC derivation. After 20 days, Nanog:H2B-GFPTg/+ ESCs were established in the presence of MEF feeders in serum+LIF conditions and then
propagated in the presence or absence of MEFs.
(E) Immunostaining and analysis of derived Nanog:H2B-GFPTg/+ ESCs. Orange arrowheads mark GFP-low/NANOG-low/OCT4+ cells.
(F) FACS analysis of derived Nanog:H2B-GFPTg/+ ESCs grown in serum+LIF conditions in the presence or absence of MEFs. Numbers in FACS histograms
indicate percentage of GFP (left) and GFP+ (right) populations. r = Pearson correlation coefcient. Scale bar represents 20 mm.
1510 Cell Reports 10, 15081520, March 10, 2015 2015 The Authors
Figure
2. Differential
Expression
of
Nanog:H2B-GFP Reporter in Segregating
EPI and PrE Cells
(A) MINS image analysis pipeline (see Experimental Procedures). GFP and NANOG expression
plotted for individual cells after automated nuclear
segmentation and uorescence intensity normalization (best-t regression curve).
(B) Immunouorescence images of xed Nanog:
H2B-GFPTg/+ embryos. EPI cells, white arrowheads; PrE cells, orange arrowheads; TE cells expressing GFP, yellow arrowheads.
(C) Quantitative immunouorescence analyses.
(D) Reporter, GATA6, and NANOG expression
in Nanogb-geo/b-geo blastocysts carrying the
Nanog:H2B-GFP reporter. NANOG staining was
absent in Nanog mutant transgenic embryos.
r = Pearson correlation coefcient. Scale bar represents 20 mm.
Cell Reports 10, 15081520, March 10, 2015 2015 The Authors 1511
cells emerging within the ICM and provides an accurate quantitative readout of
Nanog expression.
1512 Cell Reports 10, 15081520, March 10, 2015 2015 The Authors
Figure 4. EPI Cells Emerge In Vivo with an Accompanying Increase in Nanog Expression
(A) Schematic of mouse blastocyst stage embryo development.
(B) Quantication of GFP intensity in individual nuclei of a living Nanog:H2B-GFPTg/+ blastocyst recovered at E3.5 (65 cells) and corresponding images of single
time points. Single-cell intensity traces were obtained by tracking single cells in time-lapse experiments. Mitotic and apoptotic cells were also tracked, resulting in
abrupt peaks in uorescence intensity. Cells belonging to the ICM are depicted in purple, while TE cells are depicted in green; cells having undergone apoptosis
are depicted as black lines. The developmental timing and the time spanned by tracks analyzed in this panel are illustrated by an arrow in (A).
(C) Top panel, quantication of relative GFP intensity in daughters/sisters versus mother cells upon division in cells of TE (green bars) or ICM (purple bars) origin.
Daughter cells retained expression levels after mitosis, with most changes within 20%.
Data were obtained from the analysis of time-lapse movies of ve or more embryos. Bottom panel, cell divisions of EPI progenitors occurring at late blastocyst
stages shown in (B). Red arrowheads identify apoptotic events, and yellow arrowheads mark cell divisions. Scale bar represents 20 mm.
Using this methodology, we rst observed that in early blastocysts (until 60 cell stage) all cells (ICM and TE) displayed highly
heterogeneous reporter activity (0400 min in Figure 4B), consistent with previous observations (Dietrich and Hiiragi, 2007;
Ohnishi et al., 2014; Plusa et al., 2008). Following this phase,
TE cells downregulated reporter expression, whereas ICM cells
retained or increased their expression, presumably coinciding
with the establishment of a lineage bias (400 min to end of movie;
Figure 4B). We also noted apoptotic events, occurring during the
process of lineage specication in randomly positioned GFP-hi
or GFP-low ICM cells (red arrowheads; Figure 4B). As described
previously, apoptotic events occur during blastocyst development and are not a consequence of in vitro culture or phototoxicity (Artus et al., 2013; Plusa et al., 2008). Furthermore, we noted
several cell divisions occurring in GFP-hi ICM cells (yellow arrowheads; Figure 4B). Heritability of Nanog:H2B-GFP levels, and
likely cell fate, was observed in all ICM cells after division (Figure 4C). Of note, an increase in GFP uorescence was routinely
observed during mitosis, due to chromosome condensation;
thus any resulting rapid oscillations (of the order of 150 min) in reporter activity were excluded from the analysis, as they were not
Cell Reports 10, 15081520, March 10, 2015 2015 The Authors 1513
1514 Cell Reports 10, 15081520, March 10, 2015 2015 The Authors
Figure 5. State Reversals Can Occur toward a Pluripotent Identity during Lineage Specication in Mid Blastocysts, but Cells Do Not Change
Fate in Late Embryos
(A) Quantication of GFP intensity in single nuclei of a living Nanog:H2B-GFPTg/+ at E3.5 (54 cells) and corresponding snapshots from a time-lapse movie. A
subset of tracks are detailed in the lower plot, corresponding to cells highlighted in the images in the panel on the right. Cell highlighted with an orange dot and red
outline represents a GFP-low cell that upregulated reporter expression and contributed to EPI.
(B) Quantication of GFP intensity in single nuclei of a living PdgfraH2B-GFP/+ at E3.5 (60 cells) and corresponding snapshot of uorescence channel overlapped
with bright-eld images from a time-lapse movie. Only PrE-biased cells express H2B-GFP.
(C) The same analysis applied in (A) was repeated in a later Nanog:H2B-GFPTg/+ embryo at E3.75 (85 cells). ICM cells segregated toward EPI (red) and PrE
(blue) lineages. Cells that underwent apoptosis depicted with black lines. Red arrowheads mark apoptotic events and yellow arrowheads mark cell divisions.
Tracked nuclei are highlighted by dots; the outline of each dot depicts the lineage choice; red outline depicts EPI progenitor cells; blue outline depicts PrE
progenitor cells. Scale bar represents 20 mm.
Cell Reports 10, 15081520, March 10, 2015 2015 The Authors 1515
Figure 6. Cell-State Reversals toward a Pluripotent Identity Are Associated with Changes in Position within ICM
Fluorescence levels of single ICM cells (A, left, and B, lower-left panels) and two-dimensional projections of representative space trajectories (A, central, and B,
lower-central panels) are plotted for the duration of time-lapse movies of Nanog:H2B-GFPTg/+ embryo (A) corresponding to embryo in Figure 5A and (B) in the
presence of ERKi. The distance of each cell from the barycenter of the embryo is plotted versus GFP intensity (A, right, and B, lower right).
(A) EPI cells, green and blue; PrE cell, red trajectories. Blue trajectory represents a cell that switched state acquiring an EPI identity, moving from an initial position
on the surface of ICM inward while increasing Nanog expression.
(B) Top, quantication of GFP uorescence intensities and corresponding snapshots. Bottom, all six trajectories depict the behavior of ICM cells that acquire
an EPI fate increasing reporter activity. All cells change their absolute position due to embryo growth, but there are no events of abrupt spatial change in
position of individual cells within the cohort. Lines represent cell trajectories. Displayed trajectories for mitotic cells were excluded. Unique colors identify
individual cells. The outline of the embryo is drawn for plotting relative position of nuclei. Shades of gray in the trajectory depict time: earlier (light) to later
(darker).
1516 Cell Reports 10, 15081520, March 10, 2015 2015 The Authors
Cell Reports 10, 15081520, March 10, 2015 2015 The Authors 1517
Figure 7. Cell Behaviors during the Emergence of the Pluripotent EPI Lineage
Schematic representation of embryo development from mid-to-late blastocyst. Prior to lineage specication (around 6480 cells): (1) apoptosis occurs randomly
in EPI or PrE progenitors as they segregate to their appropriate layers, (2) segregation is linked to spatial movements within the ICM, and (3) few GFP-low cells
might acquire a pluripotent identity and migrate inward the ICM. After lineage specication (>8090 cells): (1) cells do not uctuate between EPI and PrE states,
and (2) cell divisions occur in EPI-specied population while relocating to the interior of the ICM.
1518 Cell Reports 10, 15081520, March 10, 2015 2015 The Authors
from a niche, namely, neighboring cells and extra-cellular components. Another could be due to the differential timescales.
ICM cells commit to PrE and EPI fates in less than 24 hr, and,
since development proceeds unidirectionally, even though
there may be non-differentiating (or symmetric) cell divisions,
neither cell population arising within the ICM self-renews. By
contrast, ESCs can be maintained indenitely in vitro under
conditions of self-renewal, a timescale that presents an unrestricted period for conversion between alternative states.
minimum in the distance between the cell centroid and the centroids of all cells
in the subsequent frame. If such minimum does not exist, the algorithm uses
cell segmentation to look for maximum overlap between cells at the subsequent time point and the cell being tracked. As additional criteria, changes
in nuclear shape and overall distance are required to lie within control ranges.
Such criteria generate unambiguous cell assignments and provide reliable
tracking of a large fraction of cells (60%70% in this study). The algorithm
was validated by visual inspection, and each track was obtained by matching
forward (prospective) and reverse (retrospective) tracking data (details provided in the Supplemental Experimental Procedures).
EXPERIMENTAL PROCEDURES
SUPPLEMENTAL INFORMATION
Mouse Husbandry
All animal experiments were approved by the Institutional Animal Care and
Use Committee at the Memorial Sloan Kettering Cancer Center. Mice were
maintained under a 12-hr light-dark cycle. Mouse lines used in this study
were Nanog:H2B-GFPTg/+, Nanogb-geo/+ (Mitsui et al., 2003), and NanogGFP/+
(Hatano et al., 2005). Alleles are schematized in Figure S1H.
AUTHOR CONTRIBUTIONS
P.X. and A.-K.H. conceived the project and designed the experiments. P.X.
generated the Nanog:H2B-GFP mouse line and performed ESC and xed
embryo imaging experiments. M.K. performed time-lapse embryo imaging experiments. A.P. and S.D.T. carried out quantitative analyses of experimental
data sets. P.X., M.K., and A.-K.H. wrote the manuscript with input from A.P.
and S.D.T.
ACKNOWLEDGMENTS
We thank X. Lou for developing the algorithm used to calculate regression curves for uorescence intensity normalizations; J. Nichols for the
Nanogb-geo/+ mouse strain; A. Martinez-Arias for TNGA ESCs; S. Nowotschin,
N. Saiz, and N. Schrode for discussions and comments on the manuscript;
the Memorial Sloan Kettering Molecular Cytology and Rockefeller University
Bio-Imaging Core Facilities for use of their instruments for live embryo imaging. Work in A.-K.Hs. laboratory is supported by the NIH (R01-HD052115
and R01-DK084391) and NYSTEM (N13G-236). S.D.T. is supported by NIH
(R00-HD074670). A.P. is supported by Finalized Research and Founding for
Investments in Basic Research (RBAP11BYNP-Newton).
Received: July 14, 2014
Revised: January 4, 2015
Accepted: January 31, 2015
Published: March 5, 2015
REFERENCES
Abranches, E., Bekman, E., and Henrique, D. (2013). Generation and characterization of a novel mouse embryonic stem cell line with a dynamic reporter of
Nanog expression. PLoS ONE 8, e59928.
Artus, J., Kang, M., Cohen-Tannoudji, M., and Hadjantonakis, A.K. (2013).
PDGF signaling is required for primitive endoderm cell survival in the inner
cell mass of the mouse blastocyst. Stem Cells 31, 19321941.
Bessonnard, S., De Mot, L., Gonze, D., Barriol, M., Dennis, C., Goldbeter, A.,
Dupont, G., and Chazaud, C. (2014). Gata6, Nanog and Erk signaling control
cell fate in the inner cell mass through a tristable regulatory network. Development 141, 36373648.
Boroviak, T., Loos, R., Bertone, P., Smith, A., and Nichols, J. (2014). The ability
of inner-cell-mass cells to self-renew as embryonic stem cells is acquired
following epiblast specication. Nat. Cell Biol. 16, 516528.
Chambers, I., Colby, D., Robertson, M., Nichols, J., Lee, S., Tweedie, S., and
Smith, A. (2003). Functional expression cloning of Nanog, a pluripotency
sustaining factor in embryonic stem cells. Cell 113, 643655.
Chambers, I., Silva, J., Colby, D., Nichols, J., Nijmeijer, B., Robertson, M.,
Vrana, J., Jones, K., Grotewold, L., and Smith, A. (2007). Nanog safeguards
pluripotency and mediates germline development. Nature 450, 12301234.
Cell Reports 10, 15081520, March 10, 2015 2015 The Authors 1519
Chazaud, C., Yamanaka, Y., Pawson, T., and Rossant, J. (2006). Early lineage
segregation between epiblast and primitive endoderm in mouse blastocysts
through the Grb2-MAPK pathway. Dev. Cell 10, 615624.
Czechanski, A., Byers, C., Greenstein, I., Schrode, N., Donahue, L.R., Hadjantonakis, A.K., and Reinholdt, L.G. (2014). Derivation and characterization of
mouse embryonic stem cells from permissive and nonpermissive strains.
Nat. Protoc. 9, 559574.
Dietrich, J.E., and Hiiragi, T. (2007). Stochastic patterning in the mouse preimplantation embryo. Development 134, 42194231.
Eckersley-Maslin, M.A., Thybert, D., Bergmann, J.H., Marioni, J.C., Flicek, P.,
and Spector, D.L. (2014). Random monoallelic gene expression increases
upon embryonic stem cell differentiation. Dev. Cell 28, 351365.
Faddah, D.A., Wang, H., Cheng, A.W., Katz, Y., Buganim, Y., and Jaenisch, R.
(2013). Single-cell analysis reveals that expression of nanog is biallelic and
equally variable as that of other pluripotency factors in mouse ESCs. Cell
Stem Cell 13, 2329.
Filipczyk, A., Gkatzis, K., Fu, J., Hoppe, P.S., Lickert, H., Anastassiadis, K., and
Schroeder, T. (2013). Biallelic expression of nanog protein in mouse embryonic
stem cells. Cell Stem Cell 13, 1213.
Frankenberg, S., Gerbe, F., Bessonnard, S., Belville, C., Pouchin, P., Bardot,
O., and Chazaud, C. (2011). Primitive endoderm differentiates via a threestep mechanism involving Nanog and RTK signaling. Dev. Cell 21, 10051013.
Gendrel, A.V., Attia, M., Chen, C.J., Diabangouaya, P., Servant, N., Barillot, E.,
and Heard, E. (2014). Developmental dynamics and disease potential of
random monoallelic gene expression. Dev. Cell 28, 366380.
ska, K., Saiz, N., Piliszek, A., Frankenberg, S., Nichols,
Grabarek, J.B., Zyzyn
J., Hadjantonakis, A.K., and Plusa, B. (2012). Differential plasticity of epiblast
and primitive endoderm precursors within the ICM of the early mouse embryo.
Development 139, 129139.
Guo, G., Huss, M., Tong, G.Q., Wang, C., Li Sun, L., Clarke, N.D., and Robson,
P. (2010). Resolution of cell fate decisions revealed by single-cell gene expression analysis from zygote to blastocyst. Dev. Cell 18, 675685.
Hatano, S.Y., Tada, M., Kimura, H., Yamaguchi, S., Kono, T., Nakano, T.,
Suemori, H., Nakatsuji, N., and Tada, T. (2005). Pluripotential competence of
cells associated with Nanog activity. Mech. Dev. 122, 6779.
Kalmar, T., Lim, C., Hayward, P., Munoz-Descalzo, S., Nichols, J., GarciaOjalvo, J., and Martinez Arias, A. (2009). Regulated uctuations in nanog
expression mediate cell fate decisions in embryonic stem cells. PLoS Biol. 7,
e1000149.
Kang, M., Piliszek, A., Artus, J., and Hadjantonakis, A.K. (2013). FGF4 is
required for lineage restriction and salt-and-pepper distribution of primitive
endoderm factors but not their initial expression in the mouse. Development
140, 267279.
Kurimoto, K., Yabuta, Y., Ohinata, Y., Ono, Y., Uno, K.D., Yamada, R.G., Ueda,
H.R., and Saitou, M. (2006). An improved single-cell cDNA amplication
method for efcient high-density oligonucleotide microarray analysis. Nucleic
Acids Res. 34, e42.
Lou, X., Kang, M., Xenopoulos, P., Munoz-Descalzo, S., and Hadjantonakis,
A.K. (2014). A rapid and efcient 2D/3D nuclear segmentation method for analysis of early mouse embryo and stem cell image data. Stem Cell Reports 2,
382397.
Mitsui, K., Tokuzawa, Y., Itoh, H., Segawa, K., Murakami, M., Takahashi, K.,
Maruyama, M., Maeda, M., and Yamanaka, S. (2003). The homeoprotein
Nanog is required for maintenance of pluripotency in mouse epiblast and ES
cells. Cell 113, 631642.
Morgani, S.M., Canham, M.A., Nichols, J., Sharov, A.A., Migueles, R.P., Ko,
M.S., and Brickman, J.M. (2013). Totipotent embryonic stem cells arise in
ground-state culture conditions. Cell Rep. 3, 19451957.
Munoz Descalzo, S., Rue, P., Garcia-Ojalvo, J., and Martinez Arias, A. (2012).
Correlations between the levels of Oct4 and Nanog as a signature for nave pluripotency in mouse embryonic stem cells. Stem Cells 30, 26832691.
Navarro, P., Festuccia, N., Colby, D., Gagliardi, A., Mullin, N.P., Zhang, W.,
Karwacki-Neisius, V., Osorno, R., Kelly, D., Robertson, M., and Chambers, I.
(2012). OCT4/SOX2-independent Nanog autorepression modulates heterogeneous Nanog gene expression in mouse ES cells. EMBO J. 31, 45474562.
Nichols, J., and Smith, A. (2012). Pluripotency in the embryo and in culture.
Cold Spring Harb. Perspect. Biol. 4, a008128.
Nichols, J., Silva, J., Roode, M., and Smith, A. (2009). Suppression of Erk signalling promotes ground state pluripotency in the mouse embryo. Development 136, 32153222.
Ohnishi, Y., Huber, W., Tsumura, A., Kang, M., Xenopoulos, P., Kurimoto, K.,
Oles, A.K., Arauzo-Bravo, M.J., Saitou, M., Hadjantonakis, A.K., and Hiiragi, T.
(2014). Cell-to-cell expression variability followed by signal reinforcement progressively segregates early mouse lineages. Nat. Cell Biol. 16, 2737.
Okita, K., Ichisaka, T., and Yamanaka, S. (2007). Generation of germlinecompetent induced pluripotent stem cells. Nature 448, 313317.
Plusa, B., Piliszek, A., Frankenberg, S., Artus, J., and Hadjantonakis, A.K.
(2008). Distinct sequential cell behaviours direct primitive endoderm formation
in the mouse blastocyst. Development 135, 30813091.
Schrode, N., Xenopoulos, P., Piliszek, A., Frankenberg, S., Plusa, B., and Hadjantonakis, A.K. (2013). Anatomy of a blastocyst: cell behaviors driving cell fate
choice and morphogenesis in the early mouse embryo. Genesis 51, 219233.
Schrode, N., Saiz, N., Di Talia, S., and Hadjantonakis, A.K. (2014). GATA6
levels modulate primitive endoderm cell fate choice and timing in the mouse
blastocyst. Dev. Cell 29, 454467.
Smith, A. (2013). Nanog heterogeneity: tilting at windmills? Cell Stem Cell
13, 67.
Takahashi, K., and Yamanaka, S. (2006). Induction of pluripotent stem cells
from mouse embryonic and adult broblast cultures by dened factors. Cell
126, 663676.
Tan, F.E., and Elowitz, M.B. (2014). Brf1 posttranscriptionally regulates pluripotency and differentiation responses downstream of Erk MAP kinase. Proc.
Natl. Acad. Sci. USA 111, E1740E1748.
Torres-Padilla, M.E., and Chambers, I. (2014). Transcription factor heterogeneity in pluripotent stem cells: a stochastic advantage. Development 141,
21732181.
Xenopoulos, P., Kang, M., and Hadjantonakis, A.K. (2012). Cell lineage allocation within the inner cell mass of the mouse blastocyst. Results Probl. Cell
Differ. 55, 185202.
MacArthur, B.D., Sevilla, A., Lenz, M., Muller, F.J., Schuldt, B.M., Schuppert,
A.A., Ridden, S.J., Stumpf, P.S., Fidalgo, M., Maayan, A., et al. (2012). Nanogdependent feedback loops regulate murine embryonic stem cell heterogeneity. Nat. Cell Biol. 14, 11391147.
Ying, Q.L., Wray, J., Nichols, J., Batlle-Morera, L., Doble, B., Woodgett, J.,
Cohen, P., and Smith, A. (2008). The ground state of embryonic stem cell
self-renewal. Nature 453, 519523.
1520 Cell Reports 10, 15081520, March 10, 2015 2015 The Authors
7\a^W`Sg]c`^SS`a
eWbV1SZZ@S^]`ba
0`SOYW\U\SeU`]c\R-4`][[]ZSQcZO`US\SbWQab]RSdSZ]^[S\bOZ
\Sc`]PW]Z]Ug1SZZ@S^]`ba^cPZWaVSabV]cUVb^`]d]YW\UQcbbW\USRUS
`SaSO`QVa^O\\W\UbVSa^SQb`c[]TZWTSaQWS\QSa
1SZZ@S^]`baWaO\]^S\OQQSaaX]c`\OZbVOb]TTS`abVS_cOZWbg`WU]`O\R
dWaWPWZWbgg]cdSQ][Sb]Sf^SQbT`][1SZZ>`Saa
2]g]cVOdSO\SePW]Z]UWQOZW\aWUVb-AcP[Wbg]c`^O^S`b]1SZZ@S^]`ba
4]`[]`SW\T]`[ObW]\dWaWbeeeQSZZQ][QSZZ`S^]`ba
/\]^S\OQQSaaX]c`\OZeWbVW[^OQb
4`][1SZZ>`Saa
Cell Reports
Article
Direct Activation of STING in the Tumor
Microenvironment Leads to Potent
and Systemic Tumor Regression and Immunity
Leticia Corrales,1,3 Laura Hix Glickman,2,3 Sarah M. McWhirter,2 David B. Kanne,2 Kelsey E. Sivick,2 George E. Katibah,2
Seng-Ryong Woo,1 Edward Lemmens,2 Tamara Banda,2 Justin J. Leong,2 Ken Metchette,2 Thomas W. Dubensky, Jr.,2,4,*
and Thomas F. Gajewski1,4,*
1Department
of Pathology, The University of Chicago, 929 E57th Street GCIS 3H, Chicago, IL 60637, USA
BioTech, Inc., 626 Bancroft Way, 3C, Berkeley, CA 94710, USA
3Co-rst author
4Co-senior author
*Correspondence: [email protected] (T.W.D.), [email protected] (T.F.G.)
http://dx.doi.org/10.1016/j.celrep.2015.04.031
This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).
2Aduro
SUMMARY
Spontaneous tumor-initiated T cell priming is dependent on IFN-b production by tumor-resident dendritic cells. On the basis of recent observations indicating that IFN-b expression was dependent upon
activation of the host STING pathway, we hypothesized that direct engagement of STING through intratumoral (IT) administration of specic agonists would
result in effective anti-tumor therapy. After proof-ofprinciple studies using the mouse STING agonist
DMXAA showed a potent therapeutic effect, we
generated synthetic cyclic dinucleotide (CDN) derivatives that activated all human STING alleles as
well as murine STING. IT injection of STING agonists
induced profound regression of established tumors
in mice and generated substantial systemic immune
responses capable of rejecting distant metastases
and providing long-lived immunologic memory. Synthetic CDNs have high translational potential as a
cancer therapeutic.
INTRODUCTION
The responsiveness of tumors to immunotherapy depends, at
least in part, on the immunophenotype of the tumor microenvironment (TME) (Gajewski et al., 2013). Substantial evidence indicates that tumor-inltrating lymphocytes (TILs) are correlated
with favorable prognosis in diverse malignancies (Galon et al.,
2012) and predicts a positive clinical outcome in response to
several immunotherapy strategies (Postow et al., 2012; Wolchok
et al., 2013). Understanding the underlying mechanisms that
promote spontaneous T cell inltration is critical toward developing new therapeutic strategies that can be used to effectively
promote an immune-responsive TME.
Innate immune sensing in the TME is a critical step in promoting spontaneous tumor-initiated T cell priming and subsequent
1018 Cell Reports 11, 10181030, May 19, 2015 2015 The Authors
TIL inltration (Fuertes et al., 2011). Transcriptional proling analyses of melanoma patients have revealed that tumors containing inltrating activated T cells are characterized by a type I IFN
transcriptional signature (Harlin et al., 2009). Studies in mice
have demonstrated that type I IFN signaling plays a critical role
in tumor-initiated T cell priming (Diamond et al., 2011; Fuertes
et al., 2011). Mice lacking the IFN-a/b receptor in DCs cannot
reject immunogenic tumors, and CD8a+ DCs from these mice
are defective in antigen cross-presentation to CD8+ T cells.
Furthermore, Baft3/ mice that lack the CD8a+ DC lineage
lose the capacity to spontaneously prime tumor-specic CD8+
T cells (Fuertes et al., 2011; Hildner et al., 2008). These ndings
in humans and in mice indicate that the tumor-resident antigenpresenting cell (APC) compartment is defective in non-T-cell-inamed tumors. Thus, strategies to induce type I IFN signaling
and APC activation in the TME to bridge the innate and adaptive
immune responses may have therapeutic utility.
Recent work has demonstrated that activation of the STING
pathway in tumor-resident host APCs is required for induction
of a spontaneous CD8+ T cell response against tumor-derived
antigens in vivo (Woo et al., 2014). In addition, activation of this
pathway and the subsequent production of IFN-b contributes
to the anti-tumor effect of radiation (Deng et al., 2014), which
can be potentiated with co-administration of a natural STING
agonist. STING (stimulator of interferon genes, also known as
TMEM173, MITA, ERIS, and MPYS) is a transmembrane protein
localized to the ER that undergoes a conformational change in
response to direct binding of cyclic dinucleotides (CDNs), resulting in a downstream signaling cascade involving TBK1 activation, IRF-3 phosphorylation, and production of IFN-b and other
cytokines (Burdette et al., 2011; Burdette and Vance, 2013; Ishikawa and Barber, 2008). IFN-b is the signature cytokine induced
in response to activating STING, by either exogenous CDNs produced by bacterial infection or through binding of a structurally
distinct endogenous CDN produced by a host cyclic GMPAMP synthetase (cGAS) in response to sensing cytosolic double-stranded DNA (dsDNA) (Ablasser et al., 2013; Diner et al.,
2013; McWhirter et al., 2009; Sun et al., 2013; Woodward
et al., 2010; Zhang et al., 2013). These observations suggested
that direct activation of the STING pathway in the TME by intratumoral (IT) injection of specic agonists might be an effective
therapeutic strategy to promote broad tumor-initiated T cell
priming against an individuals tumor antigen repertoire.
To test this therapeutic approach, we began with 5,6-dimethylxanthenone-4-acetic acid (DMXAA), a dened avonoid
compound known as a vascular disrupting agent that was shown
to have anti-tumor activity in mouse models (Baguley and Ching,
1997). This drug ultimately failed in humans when combined with
standard-of-care chemotherapy in a phase 3 efcacy trial in nonsmall-cell lung cancer (Lara et al., 2011). Interestingly, recent
structure-function studies of mouse STING (mSTING) and human STING (hSTING) demonstrated that DMXAA is a direct
ligand for mSTING (Conlon et al., 2013; Gao et al., 2013a; Kim
et al., 2013; Prantner et al., 2012). However, detailed analysis revealed that polymorphisms in hSTING rendered it unable to bind
DMXAA, therefore abrogating its activity in human cells. These
ndings provide a mechanistic insight for the lack of DMXAA efcacy in humans as well as the rationale for the development of
new pharmacologic compounds that potently activate hSTING.
Whereas STING agonists are being developed as vaccine adjuvants (Dubensky et al., 2013; Ebensen et al., 2011; Gray et al.,
2012), whether STING agonists could have direct anti-tumor
therapeutic effects has been under-explored and the lack of
dened agonists that could activate all known hSTING alleles
has been lacking.
In the current report, we conrm that DMXAA is a strong
agonist of the mSTING pathway in vitro and in vivo. We
show that IT injection of DMXAA effectively primes CD8+ T cell
responses to promote rejection of established tumors in a
STING-dependent fashion. Based on these proof-of-concept results, we synthesized a large panel of CDNs and selected compounds capable of activating all known hSTING alleles. Unlike
DMXAA, selected compounds indeed stimulated human PBMCs
to produce IFN-b. Like DMXAA, these STING agonists exhibit
signicant anti-tumor efcacy in several mouse tumor models,
without signicant local or systemic toxicity. Strikingly, direct
IT injection of selected CDNs into established B16 melanoma,
CT26 colon, and 4T1 breast carcinomas resulted in rapid and
profound tumor regression and promoted lasting systemic antigen-specic T cell immunity. These effects were entirely STINGdependent and resulted in regression of non-injected tumors in
the same hosts. We selected dithio-(RP, RP)-[cyclic[A(20 ,50 )
pA(30 ,50 )p]], (ML RR-S2 CDA) as the lead molecule for continued
development. This agent has high translational potential as a
therapeutic intervention strategy to induce activation of the
TME in multiple tumor types, with the mechanistic goal of generating effective tumor-initiated CD8+ T cell priming and lasting
anti-tumor efcacy.
RESULTS
DMXAA Stimulates the STING Pathway In Vitro
We rst conrmed that DMXAA was a functional agonist of the
STING pathway using mouse macrophages in vitro. STING aggregation was assessed using STING/ macrophages expressing mSTING-HA. Control macrophages presented a diffuse
pattern of STING in the cytoplasm, but after 1 hr of incubation
Figure 1. DMXAA Activates the STING Pathway and Promotes the Activation of APCs
(A) STING/ mouse bone-marrow-derived macrophages (BMM) transduced to express STING-HA tag were stimulated for 1 hr with 50 mg/ml DMXAA and stained
with specic antibodies against HA-tag, CD11b, and DAPI. Single cell images were acquired in the ImageStream, and data were analyzed with the IDEAS
software (Amnis; Millipore). The data in the graph represent average of percentage of cells with aggregates from three independent experiments.
(B) WT or STING/ BMM were stimulated with 50 mg/ml of DMXAA for the indicated time points. The amount of pTBK1, total TBK1, pIRF3, total IRF3, STING, and
GAPDH was measured by western blot.
(C) WT or STING/ BMM were stimulated with 50 mg/ml of DMXAA for 12 hr. The amount of secreted IFN-b was measured by ELISA.
1022 Cell Reports 11, 10181030, May 19, 2015 2015 The Authors
Figure 3. Modied CDNs Potently Activate STING and Signal through All Human STING Allelic Variants
(A) Puried human STING binding to CDNs were analyzed by thermal shift assay. Temperature curves are the average from a representative experiment of three
independent experiments performed in duplicate. Tm shift values are mean values SEM.
(B) Domain structure of hSTING is shown with the positions of the amino acid variations (bottom). The allelic frequencies of the hSTING isoforms shown on the left
hand column were obtained from the 1000 Genome Project database as previously described (Yi et al., 2013).
(C) HEK293T cells were stably transfected with the indicated STING alleles. Whole-cell lysates from HEK293T cells stably expressing the indicated full-length
STING-HA proteins were analyzed by western blot with anti-HA antibodies.
(D) HEK293T cells expressing the indicated STING alleles were transfected with an IFN-b-luciferase reporter construct. After 24 hr, cells were stimulated for 6 hr
with the indicated CDN compound (10 mM) and assessed for IFN-b-reporter activity.
(E and F) Human PBMCs from donors with the indicated STING alleles were stimulated with 10 mM of the indicated CDN or 100 mg/ml DMXAA (E), or human
PBMCs from a donor homozygous for the reference variant (STINGREF/REF) were stimulated with 10 mM and 50 mM of the indicated CDN or 100 mg/ml DMXAA (F).
After a 6-hr stimulation, fold induction of IFN-b was measured by qRT-PCR and relative normalized expression was determined by comparison with untreated
controls. Results are representative of at least two independent experiments.
Cell Reports 11, 10181030, May 19, 2015 2015 The Authors 1023
EXPERIMENTAL PROCEDURES
Cells and Cell Isolations
The cells used for the in vivo experiments were: the C57BL/6-derived
melanoma cell lines B16.F10 and B16.F10.SIY (henceforth referred to as
B16.SIY), the breast cancer 4T1 cell line, and the colon cancer CT26 cell lines,
all originally purchased from ATCC. All cells were maintained at 37 C with 5%
CO2 in DMEM supplemented with 10% heat-inactivated FCS, penicillin,
streptomycin, L-arginine, L-glutamine, folic acid, and L-asparagine.
Immortalized WT and STING/ macrophages were obtained as described
in Roberson and Walker (1988). The WT macrophages were obtained from Dr.
K. Fitzgerald (University of Massachusetts). Non-immortalized macrophages
were derived from the bone marrow of WT (C57BL/6) or STING/ mice and
cultured in BMM media (RPMI media with 5% CSF, 5% FBS, 13 L-glutamine,
and 13 pen/strep) for 7 days prior to use. Human PBMCs were isolated by
density-gradient centrifugation using Ficoll-Paque Plus (GE Healthcare).
For stable overexpression of HA-STING in STING/ macrophages, the fulllength mSting-HA DNA sequence was generated by PCR. Sequence encoding
full-length mSTING was amplied from pUNOI-mSTING plasmid (Invivogen)
using a 50 primer containing an EcoRI site, gcagacGAATTCATGCCATACTC
CAACCTGCATCCAGCCATCCCACGGCCCAGAGGTCACCGCTCCAAATAT
GTAGCCCTCATCTTTCTGGTGGCCAG, and a 30 primer containing the HAtag nucleotide sequence, followed by a TGA stop codon and a NotI site, tca
catGCGGCCGCTCAGGCGTAGTCAGGCACGTCGTAAGGATAGATGAGGTC
AGTGCGGAGTGGGAGAGGCTGATCC. The mSTING-HA PCR product was
gel puried and double digested with EcoRI and NotI and then cloned into
the multiple cloning site of pMXS-IRES-GFP with Quick ligation kit (New
England Biolabs).
Cell Reports 11, 10181030, May 19, 2015 2015 The Authors 1025
1026 Cell Reports 11, 10181030, May 19, 2015 2015 The Authors
(D) WT BALB/c mice were implanted with 105 of CT26 tumor cells on both anks. On the days indicated, mice were treated in one ank only with ML RR-S2 CDA
(50 mg) or HBSS vehicle control (n = 8).
(E) WT C57BL/6 were inoculated with 5 3 104 B16.F10 melanoma cells on the right ank at day 0 and implanted i.v. with 105 cells on day 7. Naive mice were
implanted with cells i.v. only as a control. Flank tumors were treated on the days indicated with ML RR-S2 CDA (50 mg), DMXAA (150 mg), or HBSS control (n = 8).
On day 28, lungs were harvested and lung tumor nodules counted. The histogram depicts total numbers of nodules in the ML RR-S2 CDA, DMXAA, or HBSScontrol-treated mice, compared to the untreated i.v.-only tumor implanted mice. The images depict the ML RR-S2 CDA and HBSS-control-treated mice. Data are
representative of at least two independent experiments. Results are shown as mean SEM. **p < 0.01; ***p < 0.001.
Cell Reports 11, 10181030, May 19, 2015 2015 The Authors 1027
were used as controls. For the contralateral experiments, mice were implanted
on both anks and only one tumor was treated. For the B16 melanoma lung
metastasis experiments, mice were implanted on the ank with 5 3 104 cells
B16.F10 on day 0 and then injected intravenously with 1 3 105 cells on
day 7. Lungs were harvested on day 28. Administration of compounds, measurements of tumors, and counting of lung tumors were performed in a blinded
fashion.
IFN-g ELISPOT and SIY-Pentamer Staining
Splenocytes were analyzed 5 days after the rst IT injection of DMXAA or ML
RR-S2 CDA. For the ELISPOTs, 106 splenocytes were plated per well and stimulated overnight with SIY peptide (160 nM) or AH1 (1 mM) peptide, with PMA
(50 ng/ml) plus ionomycin (0.5 mM) as a positive control or medium as negative
control. Spots were developed using the BD mouse IFN-g kit according to the
manufacturers instructions, and the number of spots was measured using an
Immunospot Series 3 Analyzer and analyzed using ImmunoSpot software
(Cellular Technology). For SIY-pentamer staining, splenocytes were preincubated for 15 min with anti-CD16/32 monoclonal antibody (93) to block potential
nonspecic binding and labeled with PE-MHC class I pentamer (Proimmune)
consisting of murine H-2Kb complexed to SIYRYYGL (SIY) peptide, antiTCRb-AF700 (H57-597), anti-CD8-Pacic Blue (53-6.7), anti-CD4-Pacic
Orange (RM4-5) (all antibodies from BioLegend), and the Fixable Viability
Dye eFluor 450 (eBioscience). Stained cells were analyzed using LSR II cytometer with FACSDiva software (BD Biosciences). Data analysis was conducted
with FlowJo software (Tree Star).
Luciferase Assay
104 HEK293T cells were seeded in 96-well plates and transiently transfected
(Lipofectamine 2000) with human IFN-b rey reporter plasmid (Fitzgerald
et al., 2003) and TK-Renilla luciferase reporter for normalization. The following
day, cells were stimulated with 10 mM of each CDN or 100 mg/ml DMXAA using
digitonin permeabilization (50 mM HEPES, 100 mM KCL, 3 mM MgCl2, 0.1 mM
DTT, 85 mM sucrose, 0.2% BSA, 1 mM ATP, 0.1 mM GTP, and 10 mg/ml digitonin) to ensure uniform uptake. After 20 min, stimulation mixtures were
removed and normal media was added. After a total of 6 hr, cell lysates
were prepared and reporter gene activity measured using the Dual Luciferase
Assay System (Promega) on a Spectramax M3 luminometer.
Statistical Analysis
Students paired t test was used to calculate two-tailed p values to estimate
statistical signicance of differences between two treatment groups using
Prism 6 software. The number of mice per group and the statistically signicant
p values are labeled in the gures and/or legends with asterisks.
E.L, T.B., J.J.L., K.M., and T.W.D. are all paid employees of Aduro BioTech,
hold stock in the company, and may be inventors on patent applications
that apply to the CDN molecules described in the manuscript.
Received: November 26, 2014
Revised: February 24, 2015
Accepted: April 14, 2015
Published: May 7, 2015
REFERENCES
Ablasser, A., Goldeck, M., Cavlar, T., Deimling, T., Witte, G., Rohl, I., Hopfner,
K.P., Ludwig, J., and Hornung, V. (2013). cGAS produces a 20 -50 -linked cyclic
dinucleotide second messenger that activates STING. Nature 498, 380384.
Baguley, B.C., and Ching, L.-M. (1997). Immunomodulatory actions of xanthenone anticancer agents. BioDrugs 8, 119127.
Blank, C., Brown, I., Peterson, A.C., Spiotto, M., Iwai, Y., Honjo, T., and Gajewski, T.F. (2004). PD-L1/B7H-1 inhibits the effector phase of tumor rejection by
T cell receptor (TCR) transgenic CD8+ T cells. Cancer Res. 64, 11401145.
Brody, J.D., Ai, W.Z., Czerwinski, D.K., Torchia, J.A., Levy, M., Advani, R.H.,
Kim, Y.H., Hoppe, R.T., Knox, S.J., Shin, L.K., et al. (2010). In situ vaccination
with a TLR9 agonist induces systemic lymphoma regression: a phase I/II study.
J. Clin. Oncol. 28, 43244332.
Burdette, D.L., and Vance, R.E. (2013). STING and the innate immune
response to nucleic acids in the cytosol. Nat. Immunol. 14, 1926.
Burdette, D.L., Monroe, K.M., Sotelo-Troha, K., Iwig, J.S., Eckert, B., Hyodo,
M., Hayakawa, Y., and Vance, R.E. (2011). STING is a direct innate immune
sensor of cyclic di-GMP. Nature 478, 515518.
Cavlar, T., Deimling, T., Ablasser, A., Hopfner, K.P., and Hornung, V. (2013).
Species-specic detection of the antiviral small-molecule compound CMA
by STING. EMBO J. 32, 14401450.
Chandra, D., Quispe-Tintaya, W., Jahangir, A., Asafu-Adjei, D., Ramos, I., Sintim, H.O., Zhou, J., Hayakawa, Y., Karaolis, D.K., and Gravekamp, C. (2014).
STING ligand c-di-GMP improves cancer vaccination against metastatic
breast cancer. Cancer Immunol Res 2, 901910.
Conlon, J., Burdette, D.L., Sharma, S., Bhat, N., Thompson, M., Jiang, Z., Rathinam, V.A.K., Monks, B., Jin, T., Xiao, T.S., et al. (2013). Mouse, but not human STING, binds and signals in response to the vascular disrupting agent
5,6-dimethylxanthenone-4-acetic acid. J. Immunol. 190, 52165225.
SUPPLEMENTAL INFORMATION
Deng, L., Liang, H., Xu, M., Yang, X., Burnette, B., Arina, A., Li, X.-D., Mauceri,
H., Beckett, M., Darga, T., et al. (2014). STING-dependent cytosolic DNA
sensing promotes radiation-induced type I interferon-dependent antitumor immunity in immunogenic tumors. Immunity 41, 843852.
Diamond, M.S., Kinder, M., Matsushita, H., Mashayekhi, M., Dunn, G.P., Archambault, J.M., Lee, H., Arthur, C.D., White, J.M., Kalinke, U., et al. (2011). Type
I interferon is selectively required by dendritic cells for immune rejection of tumors. J. Exp. Med. 208, 19892003.
AUTHOR CONTRIBUTIONS
L.C., L.H.G., T.W.D., and T.F.G. conceived the research and conducted the
experiments with S.M.M., D.B.K., K.E.S., E.L., J.J.L., G.E.K., T.B., S.-R.W.,
and K.M. T.W.D., and T.F.G. conceived and supervised the entire project
and wrote and revised the manuscript with L.C., L.H.G., K.E.S., and S.M.M.
ACKNOWLEDGMENTS
We thank Meredith Leong, Pete Lauer, and Russell Vance for discussions;
Jake Bruml, Ryan Duggan, Michael Y.K. Leung, and Diana Ranoa for technical
assistance; and Elie Diner for retroviral constructs. We thank Hector Nolla at
the UC Berkeley Cancer Research Laboratory (CRL) for help with cell sorting.
L.C. was supported by a post-doctoral fellowship from the Ramon Areces
Foundation. This work was supported by P01 CA97296 and R01CA181160
from the National Cancer Institute. L.H.G., S.M.M., D.B.K., K.E.S., G.E.K.,
1028 Cell Reports 11, 10181030, May 19, 2015 2015 The Authors
Diner, E.J., Burdette, D.L., Wilson, S.C., Monroe, K.M., Kellenberger, C.A., Hyodo, M., Hayakawa, Y., Hammond, M.C., and Vance, R.E. (2013). The innate
immune DNA sensor cGAS produces a noncanonical cyclic dinucleotide that
activates human STING. Cell Rep. 3, 13551361.
Dubensky, T.W., Jr., Kanne, D.B., and Leong, M.L. (2013). Rationale, progress
and development of vaccines utilizing STING-activating cyclic dinucleotide
adjuvants. Ther Adv Vaccines 1, 131143.
Ebensen, T., Schulze, K., Riese, P., Link, C., Morr, M., and Guzman, C.A.
(2007a). The bacterial second messenger cyclic diGMP exhibits potent adjuvant properties. Vaccine 25, 14641469.
Ebensen, T., Schulze, K., Riese, P., Morr, M., and Guzman, C.A. (2007b). The
bacterial second messenger cdiGMP exhibits promising activity as a mucosal
adjuvant. Clin. Vaccine Immunol. 14, 952958.
Ebensen, T., Libanova, R., Schulze, K., Yevsa, T., Morr, M., and Guzman, C.A.
(2011). Bis-(30 ,50 )-cyclic dimeric adenosine monophosphate: strong Th1/Th2/
Th17 promoting mucosal adjuvant. Vaccine 29, 52105220.
Fitzgerald, K.A., McWhirter, S.M., Faia, K.L., Rowe, D.C., Latz, E., Golenbock,
D.T., Coyle, A.J., Liao, S.M., and Maniatis, T. (2003). IKKepsilon and TBK1 are
essential components of the IRF3 signaling pathway. Nat. Immunol. 4,
491496.
Fuertes, M.B., Kacha, A.K., Kline, J., Woo, S.R., Kranz, D.M., Murphy, K.M.,
and Gajewski, T.F. (2011). Host type I IFN signals are required for antitumor
CD8+ T cell responses through CD8alpha+ dendritic cells. J. Exp. Med. 208,
20052016.
Gaffney, B.L., Veliath, E., Zhao, J., and Jones, R.A. (2010). One-ask syntheses of c-di-GMP and the [Rp,Rp] and [Rp,Sp] thiophosphate analogues. Org.
Lett. 12, 32693271.
Gajewski, T.F., Woo, S.R., Zha, Y., Spaapen, R., Zheng, Y., Corrales, L., and
Spranger, S. (2013). Cancer immunotherapy strategies based on overcoming
barriers within the tumor microenvironment. Curr. Opin. Immunol. 25, 268
276.
Galon, J., Page`s, F., Marincola, F.M., Angell, H.K., Thurin, M., Lugli, A., Zlobec,
I., Berger, A., Bifulco, C., Botti, G., et al. (2012). Cancer classication using the
Immunoscore: a worldwide task force. J. Transl. Med. 10, 205.
Gao, P., Ascano, M., Zillinger, T., Wang, W., Dai, P., Serganov, A.A., Gaffney,
B.L., Shuman, S., Jones, R.A., Deng, L., et al. (2013a). Structure-function analysis of STING activation by c[G(20 ,50 )pA(30 ,50 )p] and targeting by antiviral
DMXAA. Cell 154, 748762.
Gao, P., Ascano, M., Wu, Y., Barchet, W., Gaffney, B.L., Zillinger, T., Serganov,
A.A., Liu, Y., Jones, R.A., Hartmann, G., et al. (2013b). Cyclic [G(20 ,50 )pA(30 ,50 )
p] is the metazoan second messenger produced by DNA-activated cyclic
GMP-AMP synthase. Cell 153, 10941107.
Goins, W.F., Huang, S., Cohen, J.B., and Glorioso, J.C. (2014). Engineering
HSV-1 vectors for gene therapy. Methods Mol. Biol. 1144, 6379.
Gray, P.M., Forrest, G., Wisniewski, T., Porter, G., Freed, D.C., DeMartino,
J.A., Zaller, D.M., Guo, Z., Leone, J., Fu, T.M., and Vora, K.A. (2012). Evidence
for cyclic diguanylate as a vaccine adjuvant with novel immunostimulatory activities. Cell. Immunol. 278, 113119.
Harlin, H., Meng, Y., Peterson, A.C., Zha, Y., Tretiakova, M., Slingluff, C.,
McKee, M., and Gajewski, T.F. (2009). Chemokine expression in melanoma
metastases associated with CD8+ T-cell recruitment. Cancer Res. 69, 3077
3085.
Hildner, K., Edelson, B.T., Purtha, W.E., Diamond, M., Matsushita, H., Kohyama, M., Calderon, B., Schraml, B.U., Unanue, E.R., Diamond, M.S., et al.
(2008). Batf3 deciency reveals a critical role for CD8alpha+ dendritic cells
in cytotoxic T cell immunity. Science 322, 10971100.
Ishikawa, H., and Barber, G.N. (2008). STING is an endoplasmic reticulum
adaptor that facilitates innate immune signalling. Nature 455, 674678.
Ishikawa, H., Ma, Z., and Barber, G.N. (2009). STING regulates intracellular
DNA-mediated, type I interferon-dependent innate immunity. Nature 461,
788792.
Jassar, A.S., Suzuki, E., Kapoor, V., Sun, J., Silverberg, M.B., Cheung, L., Burdick, M.D., Strieter, R.M., Ching, L.M., Kaiser, L.R., and Albelda, S.M. (2005).
Activation of tumor-associated macrophages by the vascular disrupting agent
5,6-dimethylxanthenone-4-acetic acid induces an effective CD8+ T-cell-mediated antitumor immune response in murine models of lung cancer and mesothelioma. Cancer Res. 65, 1175211761.
Jin, L., Xu, L.G., Yang, I.V., Davidson, E.J., Schwartz, D.A., Wurfel, M.M., and
Cambier, J.C. (2011). Identication and characterization of a loss-of-function
human MPYS variant. Genes Immun. 12, 263269.
Joseph, W.R., Cao, Z., Mountjoy, K.G., Marshall, E.S., Baguley, B.C., and
Ching, L.M. (1999). Stimulation of tumors to synthesize tumor necrosis factor-alpha in situ using 5,6-dimethylxanthenone-4-acetic acid: a novel
approach to cancer therapy. Cancer Res. 59, 633638.
Kawarada, Y., Ganss, R., Garbi, N., Sacher, T., Arnold, B., and Hammerling, G.J.
(2001). NK- and CD8(+) T cell-mediated eradication of established tumors by
peritumoral injection of CpG-containing oligodeoxynucleotides. J. Immunol.
167, 52475253.
Kim, S., Li, L., Maliga, Z., Yin, Q., Wu, H., and Mitchison, T.J. (2013). Anticancer
avonoids are mouse-selective STING agonists. ACS Chem. Biol. 8, 1396
1401.
Konno, H., Konno, K., and Barber, G.N. (2013). Cyclic dinucleotides trigger
ULK1 (ATG1) phosphorylation of STING to prevent sustained innate immune
signaling. Cell 155, 688698.
Lara, P.N., Jr., Douillard, J.Y., Nakagawa, K., von Pawel, J., McKeage, M.J.,
Albert, I., Losonczy, G., Reck, M., Heo, D.S., Fan, X., et al. (2011). Randomized
phase III placebo-controlled trial of carboplatin and paclitaxel with or without
the vascular disrupting agent vadimezan (ASA404) in advanced non-small-cell
lung cancer. J. Clin. Oncol. 29, 29652971.
Li, L., Yin, Q., Kuss, P., Maliga, Z., Millan, J.L., Wu, H., and Mitchison, T.J.
(2014). Hydrolysis of 230 -cGAMP by ENPP1 and design of nonhydrolyzable
analogs. Nat. Chem. Biol. 10, 10431048.
McWhirter, S.M., Barbalat, R., Monroe, K.M., Fontana, M.F., Hyodo, M.,
Joncker, N.T., Ishii, K.J., Akira, S., Colonna, M., Chen, Z.J., et al. (2009). A
host type I interferon response is induced by cytosolic sensing of the bacterial
second messenger cyclic-di-GMP. J. Exp. Med. 206, 18991911.
Niesen, F.H., Berglund, H., and Vedadi, M. (2007). The use of differential scanning uorimetry to detect ligand interactions that promote protein stability.
Nat. Protoc. 2, 22122221.
Postow, M.A., Harding, J., and Wolchok, J.D. (2012). Targeting immune
checkpoints: releasing the restraints on anti-tumor immunity for patients
with melanoma. Cancer J. 18, 153159.
Prantner, D., Perkins, D.J., Lai, W., Williams, M.S., Sharma, S., Fitzgerald, K.A.,
and Vogel, S.N. (2012). 5,6-Dimethylxanthenone-4-acetic acid (DMXAA) activates stimulator of interferon gene (STING)-dependent innate immune pathways and is regulated by mitochondrial membrane potential. J. Biol. Chem.
287, 3977639788.
Roberson, S.M., and Walker, W.S. (1988). Immortalization of cloned mouse
splenic macrophages with a retrovirus containing the v-raf/mil and v-myc oncogenes. Cell. Immunol. 116, 341351.
Romling, U., Galperin, M.Y., and Gomelsky, M. (2013). Cyclic di-GMP: the rst
25 years of a universal bacterial second messenger. Microbiol. Mol. Biol. Rev.
77, 152.
Sauer, J.D., Sotelo-Troha, K., von Moltke, J., Monroe, K.M., Rae, C.S., Brubaker, S.W., Hyodo, M., Hayakawa, Y., Woodward, J.J., Portnoy, D.A., and
Vance, R.E. (2011). The N-ethyl-N-nitrosourea-induced Goldenticket mouse
mutant reveals an essential function of Sting in the in vivo interferon response
to Listeria monocytogenes and cyclic dinucleotides. Infect. Immun. 79, 688
694.
Slansky, J.E., Rattis, F.M., Boyd, L.F., Fahmy, T., Jaffee, E.M., Schneck, J.P.,
Margulies, D.H., and Pardoll, D.M. (2000). Enhanced antigen-specic antitumor immunity with altered peptide ligands that stabilize the MHC-peptideTCR complex. Immunity 13, 529538.
Sun, L., Wu, J., Du, F., Chen, X., and Chen, Z.J. (2013). Cyclic GMP-AMP synthase is a cytosolic DNA sensor that activates the type I interferon pathway.
Science 339, 786791.
Wallace, A., LaRosa, D.F., Kapoor, V., Sun, J., Cheng, G., Jassar, A., Blouin,
A., Ching, L.M., and Albelda, S.M. (2007). The vascular disrupting agent,
DMXAA, directly activates dendritic cells through a MyD88-independent
mechanism and generates antitumor cytotoxic T lymphocytes. Cancer Res.
67, 70117019.
Wolchok, J.D., Kluger, H., Callahan, M.K., Postow, M.A., Rizvi, N.A., Lesokhin,
A.M., Segal, N.H., Ariyan, C.E., Gordon, R.A., Reed, K., et al. (2013). Nivolumab plus ipilimumab in advanced melanoma. N. Engl. J. Med. 369, 122133.
Woo, S.R., Fuertes, M.B., Corrales, L., Spranger, S., Furdyna, M.J., Leung,
M.Y., Duggan, R., Wang, Y., Barber, G.N., Fitzgerald, K.A., et al. (2014).
STING-dependent cytosolic DNA sensing mediates innate immune recognition of immunogenic tumors. Immunity 41, 830842.
Woodward, J.J., Iavarone, A.T., and Portnoy, D.A. (2010). c-di-AMP secreted
by intracellular Listeria monocytogenes activates a host type I interferon
response. Science 328, 17031705.
Cell Reports 11, 10181030, May 19, 2015 2015 The Authors 1029
Yan, H., Wang, X., KuoLee, R., and Chen, W. (2008). Synthesis and immunostimulatory properties of the phosphorothioate analogues of cdiGMP. Bioorg.
Med. Chem. Lett. 18, 56315634.
Zhang, X., Shi, H., Wu, J., Zhang, X., Sun, L., Chen, C., and Chen, Z.J. (2013).
Cyclic GMP-AMP containing mixed phosphodiester linkages is an endogenous high-afnity ligand for STING. Mol. Cell 51, 226235.
Yi, G., Brendel, V.P., Shu, C., Li, P., Palanathan, S., and Cheng Kao, C. (2013).
Single nucleotide polymorphisms of human STING can affect innate immune
response to cyclic dinucleotides. PLoS ONE 8, e77846.
Zhao, L., Ching, L.M., Kestell, P., and Baguley, B.C. (2002). The antitumour activity of 5,6-dimethylxanthenone-4-acetic acid (DMXAA) in TNF receptor-1
knockout mice. Br. J. Cancer 87, 465470.
1030 Cell Reports 11, 10181030, May 19, 2015 2015 The Authors
Cell Reports
Article
Manipulation of the Quorum Sensing Signal AI-2
Affects the Antibiotic-Treated Gut Microbiota
Jessica Ann Thompson,1,5 Rita Almeida Oliveira,1,5 Ana Djukovic,2 Carles Ubeda,2,3 and Karina Bivar Xavier1,4,*
2Departamento
SUMMARY
Cell Reports 10, 18611871, March 24, 2015 2015 The Authors 1861
phenomenon is likely to also contribute to the interactions between bacteria inhabiting the mammalian gut.
Many quorum sensing signals are species specic; however,
production of and responses to one molecule, autoinducer-2
(AI-2), are observed throughout the bacterial kingdom (Chen
et al., 2002; Miller et al., 2004; Pereira et al., 2013). As AI-2 produced by one species can inuence gene expression in another,
this signal can foster interspecies communication and enable
bacteria to modify behaviors such as virulence, luminescence,
and biolm formation across different species (Armbruster
et al., 2010; Cuadra-Saenz et al., 2012; Pereira et al., 2008; Xavier and Bassler, 2005a). This feature makes AI-2 an excellent
candidate for mediating cell-cell interactions in the mammalian
gut, where hundreds of bacterial species co-exist and interact.
Multiple gut-associated bacteria that encode the AI-2 synthase,
LuxS, or produce AI-2 have already been identied (Antunes
et al., 2005; Hsiao et al., 2014; Lukas et al., 2008; Schauder
et al., 2001). Furthermore, the human commensal bacterium
Ruminococcus obeum was recently reported to inhibit colonization of the mouse gut by V. cholerae, partially through AI-2
signaling, highlighting how AI-2 produced by commensal bacteria can affect invading pathogens (Hsiao et al., 2014). We hypothesized that signaling through AI-2 might also occur between
commensal members of the microbiota and asked whether
AI-2 can shape the species composition of this community under
conditions of dysbiosis.
To investigate this hypothesis, we made use of the natural
signal production and depletion capabilities of the enteric bacterium, Escherichia coli. This organism secretes large amounts of
AI-2 into the environment; it also harbors a highly efcient mechanism for signal uptake and degradation, known as the Lsr transport system (Pereira et al., 2012; Xavier and Bassler, 2005b). By
internalizing and processing AI-2 produced by itself as well as
from other species, Lsr system-expressing bacteria can disrupt
the ability of neighboring species to correctly determine population density and regulate AI-2-dependent behavior appropriately, as shown in vitro in mixed cultures of E. coli and Vibrio
spp. (Xavier and Bassler, 2005a). As a result, this system has
been explored as a potential means for AI-2 interspecies quorum
quenching (Roy et al., 2010). Given that E. coli can also stably
colonize the mouse gut following streptomycin treatment (Conway et al., 2004), we used this bacterium as a tool to manipulate
AI-2 signaling in vivo and demonstrated the accumulation and
depletion of AI-2 in the intestinal tract of gnotobiotic mice. Streptomycin induces gut dysbiosis and has been used extensively to
study Salmonella pathogenesis and E. coli physiology following
the disruption of colonization resistance (Barthel et al., 2003;
Spees et al., 2013). We characterized the changes induced by
streptomycin upon the composition of the microbiota and then
determined the effect of E. coli-mediated AI-2 manipulation
upon this antibiotic-treated community.
RESULTS
Engineered E. coli Mutants Manipulate AI-2 Availability
in the Mouse Gut
To manipulate AI-2 levels in the mouse gut, we constructed a
combination of E. coli mutants affected in the Lsr system that
1862 Cell Reports 10, 18611871, March 24, 2015 2015 The Authors
Figure 1. E. coli Accumulate and Deplete AI-2 in the Gut of Monocolonized Mice
AI-2 activity in cecal extracts harvested from germ-free mice 5 days after
gavage with PBS, WT, DlsrK or DlsrRDluxS; or a 1:1 mix of DlsrK and
DlsrRDluxS; or DlsrK and DlsrKDluxS E. coli, as measured by Vibrio harveyi
bioluminescence. Data shown are the mean, and the error bars correspond to
SD; n = 3. See also Figure S1.
Cell Reports 10, 18611871, March 24, 2015 2015 The Authors 1863
antibiotic treatment, followed by a decrease on day 28 (Figure S2D). This variability also affected the phylogenetic distance
of the communities, as the unweighted UniFrac distance similarly increased after exposure to streptomycin and then
lowered in the bacterial communities present on day 28 (Fig-
1864 Cell Reports 10, 18611871, March 24, 2015 2015 The Authors
Cell Reports 10, 18611871, March 24, 2015 2015 The Authors 1865
Figure 4. Microbiota Composition of Streptomycin-Treated Mice Differs in the Presence of DlsrK Mutant E. coli
Intestinal microbiota composition was analyzed in samples collected from the same mice as in Figure 3 after 28 days of antibiotic exposure (corresponding to
26 days of E. coli colonization; n = 9 per group).
(A and B) Analysis of the overall microbiota of the different E. coli groups by PCoA plots of Jaccard and unweighted Unifrac distances, respectively. The rst two
coordinates are shown. Each group is labeled with a different color, as indicated. Ellipses centered on the categorical averages of the metric distances with a 95%
condence interval for the rst two coordinates of each group were drawn on the associated PCoA.
(CH) Relative abundance of individual OTUs that exhibited a signicant difference between the group colonized with DlsrK and both the other two groups are
shown. Data shown are the median, and error bars show the interquartile range. Data were analyzed with the Wilcoxon test using the Benjamini-Hochberg
correction (*q < 0.1; **q < 0.05; ns, not signicant).
(I) Heatmap showing the fold change in relative abundance of OTUs in mice colonized with DlsrK mutant bacteria divided by the mean of the same OTU in mice
colonized with either the DlsrRDluxS or the control group, DlsrKDluxS, E. coli. All the OTUs that exhibited a signicant difference between the different groups are
shown.
See also Figure S4.
colonized with either DlsrRDluxS or DlsrKDluxS bacteria (Figure 5A). In contrast, the Firmicutes were positively affected
by the presence of AI-2-producing DlsrK mutant bacteria, un-
dergoing 3- to 6-fold increase to 0.8% from a relative abundance of 0.24% and 0.13% in mice colonized with DlsrKDluxS
or DlsrRDluxS mutant bacteria, respectively (Figure 5B). This
1866 Cell Reports 10, 18611871, March 24, 2015 2015 The Authors
corresponding to these two phyla, respectively, encode homologs to the AI-2 synthase (Figure 6).
In conclusion, the presence of DlsrK mutant bacteria, which
accumulate high levels of AI-2, changed the abundance of the
Bacteroidetes and Firmicutes. The latter group were favored
despite the continued presence of streptomycin, providing
evidence that the bacterial communication molecule, AI-2, can
modulate the composition of gut microbiota.
DISCUSSION
resulted in an increase in the ratio of Firmicutes to Bacteroidetes. Though this increase was signicant (Figure 5C), the ratio in these DlsrK mutant-colonized mice remained much lower
than that observed in the microbiota of untreated animals
(0.897). Interestingly, the Bacteroidetes and Firmicutes are predicted to have very different AI-2 production capabilities: 17%
and 83% of currently sequenced genomes (KEGG database)
Cell Reports 10, 18611871, March 24, 2015 2015 The Authors 1867
Figure 6. Higher Prevalence of LuxS Orthologs in the Complete Genomes of Bacteria Belonging to the Firmicutes
Percentage of full genome sequences corresponding to members of the Bacteroidetes and Firmicutes containing protein orthologs of LuxS.
1868 Cell Reports 10, 18611871, March 24, 2015 2015 The Authors
EXPERIMENTAL PROCEDURES
Bacterial Strains and Culture Conditions
All E. coli strains and primers used in this study are listed in supplemental
tables. These are all E. coli K-12 MG1655 derivatives with a streptomycinresistant mutation in rpsL-1(K43N) that express CFP or YFP constitutively.
For details of culture conditions, genetic manipulation, and strain construction,
see Supplemental Experimental Procedures.
Animal Studies
All experiments involving mice were approved by the Institutional Ethics
^ncia and the Portuguese National
Committee at the Instituto Gulbenkian de Cie
Entity (Direcao Geral de Alimentacao e Veterinaria; ref. no. 008957, approval
date 19/03/2013) following the Portuguese legislation (PORT 1005/92), which
complies with the European Directive 86/609/EEC of the European Council.
To demonstrate AI-2 production in the gut, germ-free mice were gavaged
with sterile PBS or 105 CFUs of wild-type, DlsrK, or DlsrRDluxS E. coli strains.
To demonstrate removal of AI-2 from the gut by the E. coli DlsrRDluxS mutant
strain, germ-free mice were gavaged with a 1:1 mix of 105 CFUs of DlsrK and
DlsrKDluxS or DlsrK and DlsrRDluxS mutant bacteria labeled with either CFP
or YFP. Fecal samples were collected 5 days after colonization and plated
to determine bacterial load; cecal contents were harvested for analysis of
AI-2 levels.
6- to 8-week-old male C57BL/6J mice conventionally raised under specic
pathogen-free (SPF) conditions were used to analyze the effects of streptomycin and colonization by E. coli mutant strains upon the intestinal microbiota.
To determine the effects of streptomycin treatment on the gut microbiota
composition, four non-littermate mice were housed individually and maintained
under 5 g/l streptomycin ad libitum in the drinking water (Conway et al. 2004).
Fresh fecal samples were collected prior to antibiotic administration (day 0)
and 2, 4, 7, 12, and 28 days during treatment for subsequent DNA extraction.
To assess the effects of different E. coli mutants upon the gastrointestinal ora,
ve groups of mice (n = 6) were treated and maintained under streptomycin, as
described above. On day 2 of treatment, two mice per group were gavaged with
100 ml PBS containing 108 colony-forming units (CFU) of either ARO071 (DlsrK),
ARO093 (DlsrKDluxS), or ARO081 (DlsrRDluxS) and individually caged (so that
n = 10 per treatment). Fecal samples were collected, part was homogenized in
1 ml sterile PBS and plated to determine colonization levels (CFU/g feces), and
the remainder of each sample was used for subsequent DNA extraction.
Detection of AI-2 Activity
AI-2 activity was measured as previously described (Taga and Xavier, 2011)
using the V. harveyi AI-2 reporter strain TL26 (DluxN DluxS DcqsS; Long
et al., 2009). To determine AI-2 activity in mouse cecal extracts, the cecal contents were homogenized at a 10% weight/volume concentration in 0.1 M
MOPS (pH 7). Samples were centrifuged and ltered, and then an equal volume of methanol was added to precipitate further debris. Supernatants were
vacuum-dried, resuspended at 50% weight/volume in sterile water, and
analyzed by bioassay as above. Enumeration of V. harveyi CFUs demonstrated
no signicant differences in growth of the reporter strain across the samples.
Sample Collection and DNA Extraction
DNA was extracted from fecal samples using the QIAamp DNA Stool Mini Kit
(QIAGEN) following the manufacturers instructions plus an additional membrane disruption step using 0.1-mm glass beads and high-speed shaking.
Samples were stored at 20 C. Total DNA obtained was quantied with Qubit
dsDNA BR Assay Kit (Invitrogen); samples with more than 10 ng/ml were
sequenced and further analyzed (n = 9 per treatment).
Quantication of Bacterial Load by qPCR
Quantitative PCR (qPCR) was performed on DNA extracted from fecal samples
using 16S rRNA universal primers and YFP (CFP)-specic primers to determine total bacterial and YFP-expressing E. coli loads/g feces, respectively.
16S rRNA Gene Amplication, Pyrosequencing, and Analysis
For each sample, the V1V3 region of the 16S rRNA gene was amplied by
PCR and sequenced using a 454 GS FLX Titanium platform following Roche
Cell Reports 10, 18611871, March 24, 2015 2015 The Authors 1869
Bufe, C.G., Jarchum, I., Equinda, M., Lipuma, L., Gobourne, A., Viale, A.,
Ubeda, C., Xavier, J., and Pamer, E.G. (2012). Profound alterations of intestinal
microbiota following a single dose of clindamycin results in sustained susceptibility to Clostridium difcile-induced colitis. Infect. Immun. 80, 6273.
Leatham, M.P., Banerjee, S., Autieri, S.M., Mercado-Lubo, R., Conway, T., and
Cohen, P.S. (2009). Precolonized human commensal Escherichia coli strains
serve as a barrier to E. coli O157:H7 growth in the streptomycin-treated mouse
intestine. Infect. Immun. 77, 28762886.
Bufe, C.G., Bucci, V., Stein, R.R., McKenney, P.T., Ling, L., Gobourne, A., No,
D., Liu, H., Kinnebrew, M., Viale, A., et al. (2015). Precision microbiome reconstitution restores bile acid mediated resistance to Clostridium difcile. Nature
517, 205208.
Ley, R.E., Turnbaugh, P.J., Klein, S., and Gordon, J.I. (2006). Microbial ecology: human gut microbes associated with obesity. Nature 444, 10221023.
Chang, P.V., Hao, L., Offermanns, S., and Medzhitov, R. (2014). The microbial
metabolite butyrate regulates intestinal macrophage function via histone
deacetylase inhibition. Proc. Natl. Acad. Sci. USA 111, 22472252.
Chen, X., Schauder, S., Potier, N., Van Dorsselaer, A., Pelczer, I., Bassler, B.L.,
and Hughson, F.M. (2002). Structural identication of a bacterial quorumsensing signal containing boron. Nature 415, 545549.
Conway, T., Krogfelt, K., and Cohen, P. (2004). The life of commensal Escherichia coli in the mammalian intestine. EcoSal Plus 2004 http://dx.doi.org/10.
1128/ecosalplus.8.3.1.2.
Cuadra-Saenz, G., Rao, D.L., Underwood, A.J., Belapure, S.A., Campagna,
S.R., Sun, Z., Tammariello, S., and Rickard, A.H. (2012). Autoinducer-2 inuences interactions amongst pioneer colonizing streptococci in oral biolms.
Microbiology 158, 17831795.
Long, T., Tu, K.C., Wang, Y., Mehta, P., Ong, N.P., Bassler, B.L., and Wingreen, N.S. (2009). Quantifying the integration of quorum-sensing signals with
single-cell resolution. PLoS Biol. 7, e68.
Louis, P., Hold, G.L., and Flint, H.J. (2014). The gut microbiota, bacterial
metabolites and colorectal cancer. Nat. Rev. Microbiol. 12, 661672.
Lukas, F., Gorenc, G., and Kopecny, J. (2008). Detection of possible AI-2mediated quorum sensing system in commensal intestinal bacteria. Folia
Microbiol. (Praha) 53, 221224.
Maslowski, K.M., Vieira, A.T., Ng, A., Kranich, J., Sierro, F., Yu, D., Schilter,
H.C., Rolph, M.S., Mackay, F., Artis, D., et al. (2009). Regulation of inammatory responses by gut microbiota and chemoattractant receptor GPR43.
Nature 461, 12821286.
McNeil, N.I. (1984). The contribution of the large intestine to energy supplies
in man. Am. J. Clin. Nutr. 39, 338342.
Fabich, A.J., Jones, S.A., Chowdhury, F.Z., Cernosek, A., Anderson, A.,
Smalley, D., McHargue, J.W., Hightower, G.A., Smith, J.T., Autieri, S.M.,
et al. (2008). Comparison of carbon nutrition for pathogenic and commensal
Escherichia coli strains in the mouse intestine. Infect. Immun. 76, 11431152.
Miller, S.T., Xavier, K.B., Campagna, S.R., Taga, M.E., Semmelhack, M.F.,
Bassler, B.L., and Hughson, F.M. (2004). Salmonella typhimurium recognizes
a chemically distinct form of the bacterial quorum-sensing signal AI-2. Mol.
Cell 15, 677687.
Finegold, S.M., Dowd, S.E., Gontcharova, V., Liu, C., Henley, K.E., Wolcott,
R.D., Youn, E., Summanen, P.H., Granpeesheh, D., Dixon, D., et al. (2010).
Pyrosequencing study of fecal microora of autistic and control children.
Anaerobe 16, 444453.
Ng, K.M., Ferreyra, J.A., Higginbottom, S.K., Lynch, J.B., Kashyap, P.C.,
Gopinath, S., Naidu, N., Choudhury, B., Weimer, B.C., Monack, D.M., and
Sonnenburg, J.L. (2013). Microbiota-liberated host sugars facilitate post-antibiotic expansion of enteric pathogens. Nature 502, 9699.
Flint, H.J., Duncan, S.H., Scott, K.P., and Louis, P. (2007). Interactions and
competition within the microbial community of the human colon: links between
diet and health. Environ. Microbiol. 9, 11011111.
Pereira, C.S., McAuley, J.R., Taga, M.E., Xavier, K.B., and Miller, S.T. (2008).
Sinorhizobium meliloti, a bacterium lacking the autoinducer-2 (AI-2) synthase,
responds to AI-2 supplied by other bacteria. Mol. Microbiol. 70, 12231235.
Frank, D.N., St Amand, A.L., Feldman, R.A., Boedeker, E.C., Harpaz, N., and
Pace, N.R. (2007). Molecular-phylogenetic characterization of microbial
community imbalances in human inammatory bowel diseases. Proc. Natl.
Acad. Sci. USA 104, 1378013785.
Pereira, C.S., Santos, A.J., Bejerano-Sagie, M., Correia, P.B., Marques, J.C.,
and Xavier, K.B. (2012). Phosphoenolpyruvate phosphotransferase system
regulates detection and processing of the quorum sensing signal autoinducer-2. Mol. Microbiol. 84, 93104.
Furusawa, Y., Obata, Y., Fukuda, S., Endo, T.A., Nakato, G., Takahashi, D.,
Nakanishi, Y., Uetake, C., Kato, K., Kato, T., et al. (2013). Commensal
microbe-derived butyrate induces the differentiation of colonic regulatory
T cells. Nature 504, 446450.
Pereira, C.S., Thompson, J.A., and Xavier, K.B. (2013). AI-2-mediated signalling in bacteria. FEMS Microbiol. Rev. 37, 156181.
Galley, J.D., Nelson, M.C., Yu, Z., Dowd, S.E., Walter, J., Kumar, P.S., Lyte, M.,
and Bailey, M.T. (2014). Exposure to a social stressor disrupts the community
structure of the colonic mucosa-associated microbiota. BMC Microbiol. 14,
189.
Qin, J., Li, Y., Cai, Z., Li, S., Zhu, J., Zhang, F., Liang, S., Zhang, W., Guan, Y.,
Shen, D., et al. (2012). A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature 490, 5560.
Rakoff-Nahoum, S., Paglino, J., Eslami-Varzaneh, F., Edberg, S., and Medzhitov, R. (2004). Recognition of commensal microora by toll-like receptors is
required for intestinal homeostasis. Cell 118, 229241.
Garner, C.D., Antonopoulos, D.A., Wagner, B., Duhamel, G.E., Keresztes, I.,
Ross, D.A., Young, V.B., and Altier, C. (2009). Perturbation of the small
intestine microbial ecology by streptomycin alters pathology in a Salmonella
enterica serovar typhimurium murine model of infection. Infect. Immun. 77,
26912702.
Reeves, A.E., Koenigsknecht, M.J., Bergin, I.L., and Young, V.B. (2012).
Suppression of Clostridium difcile in the gastrointestinal tracts of germfree
mice inoculated with a murine isolate from the family Lachnospiraceae. Infect.
Immun. 80, 37863794.
Roy, V., Fernandes, R., Tsao, C.-Y., and Bentley, W.E. (2010). Cross species
quorum quenching using a native AI-2 processing enzyme. ACS Chem. Biol.
5, 223232.
Hsiao, A., Ahmed, A.M.S., Subramanian, S., Grifn, N.W., Drewry, L.L., Petri,
W.A., Jr., Haque, R., Ahmed, T., and Gordon, J.I. (2014). Members of the human gut microbiota involved in recovery from Vibrio cholerae infection. Nature
515, 423426.
Kamada, N., Kim, Y.-G., Sham, H.P., Vallance, B.A., Puente, J.L., Martens,
E.C., and Nunez, G. (2012). Regulated virulence controls the ability of a pathogen to compete with the gut microbiota. Science 336, 13251329.
Ruby, E.G. (2008). Symbiotic conversations are revealed under genetic interrogation. Nat. Rev. Microbiol. 6, 752762.
Rutherford, S.T., and Bassler, B.L. (2012). Bacterial quorum sensing: its role in
virulence and possibilities for its control. Cold Spring Harb. Perspect. Med. 2,
pii: a012427.
Kamada, N., Seo, S.-U., Chen, G.Y., and Nunez, G. (2013). Role of the
gut microbiota in immunity and inammatory disease. Nat. Rev. Immunol.
13, 321335.
Salonen, A., Lahti, L., Salojarvi, J., Holtrop, G., Korpela, K., Duncan, S.H., Date,
P., Farquharson, F., Johnstone, A.M., Lobley, G.E., et al. (2014). Impact of diet
and individual variation on intestinal microbiota composition and fermentation
products in obese men. ISME J. 8, 22182230.
Savage, D.C. (1977). Microbial ecology of the gastrointestinal tract. Annu. Rev.
Microbiol. 31, 107133.
1870 Cell Reports 10, 18611871, March 24, 2015 2015 The Authors
Schauder, S., Shokat, K., Surette, M.G., and Bassler, B.L. (2001). The LuxS
family of bacterial autoinducers: biosynthesis of a novel quorum-sensing
signal molecule. Mol. Microbiol. 41, 463476.
Schloss, P.D., Westcott, S.L., Ryabin, T., Hall, J.R., Hartmann, M., Hollister,
E.B., Lesniewski, R.A., Oakley, B.B., Parks, D.H., Robinson, C.J., et al.
(2009). Introducing mothur: open-source, platform-independent, communitysupported software for describing and comparing microbial communities.
Appl. Environ. Microbiol. 75, 75377541.
Sekirov, I., Tam, N.M., Jogova, M., Robertson, M.L., Li, Y., Lupp, C., and Finlay, B.B. (2008). Antibiotic-induced perturbations of the intestinal microbiota
alter host susceptibility to enteric infection. Infect. Immun. 76, 47264736.
Taga, M.E., and Xavier, K.B. (2011). Methods for analysis of bacterial autoinducer-2 production. Curr. Protoc. Microbiol. Chapter 1, Unit 1C.1.
Turnbaugh, P.J., Ley, R.E., Mahowald, M.A., Magrini, V., Mardis, E.R., and
Gordon, J.I. (2006). An obesity-associated gut microbiome with increased
capacity for energy harvest. Nature 444, 10271031.
Turnbaugh, P.J., Quince, C., Faith, J.J., McHardy, A.C., Yatsunenko, T., Niazi,
F., Affourtit, J., Egholm, M., Henrissat, B., Knight, R., and Gordon, J.I. (2010).
Organismal, genetic, and transcriptional variation in the deeply sequenced gut
microbiomes of identical twins. Proc. Natl. Acad. Sci. USA 107, 75037508.
Smith, P.M., Howitt, M.R., Panikov, N., Michaud, M., Gallini, C.A., Bohlooly-Y,
M., Glickman, J.N., and Garrett, W.S. (2013). The microbial metabolites,
short-chain fatty acids, regulate colonic Treg cell homeostasis. Science 341,
569573.
Ubeda, C., Taur, Y., Jenq, R.R., Equinda, M.J., Son, T., Samstein, M., Viale, A.,
Socci, N.D., van den Brink, M.R.M., Kamboj, M., and Pamer, E.G. (2010).
Vancomycin-resistant Enterococcus domination of intestinal microbiota is
enabled by antibiotic treatment in mice and precedes bloodstream invasion
in humans. J. Clin. Invest. 120, 43324341.
Spees, A.M., Wangdi, T., Lopez, C.A., Kingsbury, D.D., Xavier, M.N., Winter,
S.E., Tsolis, R.M., and Baumler, A.J. (2013). Streptomycin-induced inammation enhances Escherichia coli gut colonization through nitrate respiration.
mBio. 4, e00430-13.
Ubeda, C., Bucci, V., Caballero, S., Djukovic, A., Toussaint, N.C., Equinda, M.,
Lipuma, L., Ling, L., Gobourne, A., No, D., et al. (2013). Intestinal microbiota
containing Barnesiella species cures vancomycin-resistant Enterococcus
faecium colonization. Infect. Immun. 81, 965973.
Stecher, B., Robbiani, R., Walker, A.W., Westendorf, A.M., Barthel, M.,
Kremer, M., Chaffron, S., Macpherson, A.J., Buer, J., Parkhill, J., et al.
(2007). Salmonella enterica serovar typhimurium exploits inammation to
compete with the intestinal microbiota. PLoS Biol. 5, 21772189.
Wang, T., Cai, G., Qiu, Y., Fei, N., Zhang, M., Pang, X., Jia, W., Cai, S., and
Zhao, L. (2012). Structural segregation of gut microbiota between colorectal
cancer patients and healthy volunteers. ISME J. 6, 320329.
Stecher, B., Chaffron, S., Kappeli, R., Hapfelmeier, S., Freedrich, S., Weber,
T.C., Kirundi, J., Suar, M., McCoy, K.D., von Mering, C., et al. (2010). Like
will to like: abundances of closely related species can predict susceptibility
to intestinal colonization by pathogenic and commensal bacteria. PLoS
Pathog. 6, e1000711.
Xavier, K.B., and Bassler, B.L. (2005a). Interference with AI-2-mediated bacterial cell-cell communication. Nature 437, 750753.
Xavier, K.B., and Bassler, B.L. (2005b). Regulation of uptake and processing
of the quorum-sensing autoinducer AI-2 in Escherichia coli. J. Bacteriol.
187, 238248.
Cell Reports 10, 18611871, March 24, 2015 2015 The Authors 1871
Cell Reports
Resource
Mapping Social Behavior-Induced Brain
Activation at Cellular Resolution in the Mouse
Yongsoo Kim,1 Kannan Umadevi Venkataraju,1 Kith Pradhan,1 Carolin Mende,1 Julian Taranda,1 Srinivas C. Turaga,2
Ignacio Arganda-Carreras,2 Lydia Ng,3 Michael J. Hawrylycz,3 Kathleen S. Rockland,1,4 H. Sebastian Seung,2
and Pavel Osten1,*
1Cold
SUMMARY
Understanding how brain activation mediates behaviors is a central goal of systems neuroscience. Here,
we apply an automated method for mapping brain
activation in the mouse in order to probe how sexspecic social behaviors are represented in the male
brain. Our method uses the immediate-early-gene
c-fos, a marker of neuronal activation, visualized by
serial two-photon tomography: the c-fos-GFP+ neurons are computationally detected, their distribution
is registered to a reference brain and a brain atlas,
and their numbers are analyzed by statistical tests.
Our results reveal distinct and shared female and
male interaction-evoked patterns of male brain activation representing sex discrimination and social
recognition. We also identify brain regions whose degree of activity correlates to specic features of social
behaviors and estimate the total numbers and the
densities of activated neurons per brain areas. Our
study opens the door to automated screening of
behavior-evoked brain activation in the mouse.
INTRODUCTION
Central to the understanding of brain functions is insight into the
distribution of neuronal activity that drives behavior. Local measurements of brain activity in behaving mice can be made with
electrodes and uorescent calcium indicators (Buzsaki, 2004;
Grewe and Helmchen, 2009), but such approaches provide information regarding only a very small fraction of the 70 million
neurons that comprise the mouse brain. The detection of
elevated levels of the immediate-early genes (IEGs) linked to
recent neuronal activity (Clayton, 2000; Guzowski et al., 2005)
is a more spatially comprehensive technique. While it lacks the
time resolution of electrophysiological recordings or calcium imaging, it does have the potential of providing a complete view of
recent whole-brain activity. Once determined, the whole-brain
IEG-based map can be used to generate structure-function hy292 Cell Reports 10, 292305, January 13, 2015 2015 The Authors
Figure 3. Social Behavior-Activated Areas: The MOB and Its Direct Downstream Circuitry
(A and B) ROI analysis:.The male-female (A) and male-male groups (B) are compared to the handling group and signicantly activated ROIs downstream of the
MOB are displayed. Most of the regions were activated by both stimuli. Heatmap in (A) represents FDR corrected statistical signicance. Numbers in (B) represent
bregma A/P coordinates. See Table S2 for ROI full names.
(CF) Voxel-based analysis revealed activation pattern selective for the female stimulus (red), the male stimulus (green), and shared by both stimuli (yellow). (C)
Both male and female stimuli induced dorsal activation in the MOB. (DF) Dorsoventral separation was detected between the male- and female-evoked activation
in the PIR (D and F) and ENT (E and F).
See also Movie S3 for the full data set.
prefrontal cortices in the male-male interaction was limited to supercial cortical layers (Figure 5C; Movie S3). We also observed
a focal activation in the DR at a specic A/P bregma location in
the male-female data set (Figure 5E; Movie S3).
296 Cell Reports 10, 292305, January 13, 2015 2015 The Authors
Figure 4. Social Behavior-Activated Areas: the AOB and Its Direct Downstream Circuitry
(A and B) ROI analysis. The male-female (A) and male-male groups (B) are compared to the handling group, and signicantly activated ROIs downstream of the
AOB are displayed. The female stimulus activated all AOB downstream regions, while the male stimulus induced only a partial activation of these areas. Heatmap
in (A) represents FDR-corrected statistical signicance. Numbers in (A) represent bregma A/P coordinates. See Table S2 for ROI full names.
(CE) Voxel-based analysis revealed a largely overlapping activation pattern (yellow) in the coactivated AOB (C), MEAad and MEAav (D), and MEApd (E), and
selective female-evoked activation in the MEApv and COApm (E).
See also Movie S3 for the full data set.
distinct and focal LSr activation at A/P coordinates between +0.345 and 0.145 (Figure 6C; Movie S3). The activation in the VMHvl, which was previously shown to play a role
in both sexual and aggressive behaviors (Lin et al., 2011),
was highly overlapping between the male-female and malemale data sets (Figure 6D; Movie S3). In addition, only a medial
part of the PMv was activated in the male-male data set, suggesting a functional subdivision within this structure (Figure 6E;
Movie S3).
Finally, among additional brain areas, the claustrum (CLA), basomedial amygdala (BMA), and intercalated amygdala (IA) were
activated by both the female and male interactions; the capsular
central amygdala (CEAc), basolateral amygdala (BLA), and
thalamic parataenial nucleus (PT) were activated only in
response to the female stimulus; and the temporal associational,
perirhinal, and ectorhinal (TEa, PERI, and ECT) cortical areas
were activated only in response to the male stimulus (Table
S2). The activation of the hippocampal CA2 region linked to
Cell Reports 10, 292305, January 13, 2015 2015 The Authors 297
Table 1.
ROIs
ROIs
Anogentical Snifng
Close Following
Isocortex
Anogentical Snifng
Close Following
Midbrain
ILA
DR
PL
++
ORBm
++
TT
++
Olfactory area
DP
++
AOBgr
++
COAa
COApl
+++
COApm
+++
PAA
+++
TR
Hippocampal formation
ENTmv
Continued
++
CA3
Cortical subplate
EP
BLAa
++
BLAp
BLAv
++
BMAa
++
BMAp
+++
PA
+++
Cerebral nuclei
ACBsh
OT
AAA
LSr
++
+
CEAc
IA
++
MEAad
+++
MEAav
+++
MEApd
+++
MEApv
+++
SI
BSTmg
++
BSTv
++
BSTp
++
Thalamus
MDm
PT
Hypothalamus
AVPV
++
MPN
++
PMv
+++
VMHvl
TU
++
The correlation analysis can be used to add functional signicance to activated regions that were not previously known to be
involved in social behaviors. For example, the activation of the
amygdalar IA and CEAc nuclei was correlated to the anogenital
snifng time, while the PT thalamus activation was correlated
to following. Since both IA and CEAc can inhibit the medial central amygdala (CEAm), which is the output fear pathway (Pitkanen et al., 1997), these data suggest that the IA and CEAc are
activated by chemosensory cues and may act to modulate fear
behaviors during social exploration. The PT, a part of the dorsal
group of thalamic nuclei, projects to the ACB (Kelley and Stinus,
1984) and may play a role in motivational modulation of the malefemale social behavior.
Finally, the quantication of the numbers of c-fos-GFP+ cells
per brain area can provide information about the approximate
percentage of neurons behaviorally recruited in the identied
brain areas. For example, we observed on average 2-fold increase (3,600 c-fos-GFP+ cells per cubic mm) in the prefrontal
cortical areas in the male-female data sets, compared to the
handling control (Figure S7). Since neuronal density in the mouse
cortex is estimated at 80,000100,000 per cubic mm (Herculano-Houzel et al., 2006; Keller and Carlson, 1999; Meyer et al.,
2010), these data suggest less than 5% of neurons is recruited
in response to the female social stimulus. Further, as most brain
areas showed similar c-fos-GFP+ densities, the behavioral
recruitment of a few percent of neurons is likely a general feature
of c-fos activation. This may represent c-fos induction occurring
only in the most strongly activated cells, and such sparse c-fos
induction may be relevant for the proposed sparse coding of
sensory inputs (Olshausen and Field, 2004).
Caveats of the Current Study
There are several caveats associated with our study. First, while
the behavior is limited to 90 s, our assay cannot determine
whether the observed c-fos-GFP induction occurred entirely
during this brief time period or whether some downstream activation occurred during a longer time interval. Second, the behavioral paradigm includes both the introduction and removal of the
stimulus animal from the home cage of the c-fos-GFP male, and
some of the observed activation pattern thus may reect stress
induced by these manipulations. Third, the use of the OVX females in our study restricts the interpretation of the male-female
activation data to social exploration that lacks the effects of
estrous hormones. Thus, male-female interaction with, for
example, estradiol-induced OVX mice may be expected to
induce brain activation partially distinct from the one described
in the current study. Fourth, since the c-fos-GFP reporter labels
60% of all c-fos+ cells detected by immunostaining, it is
possible that some areas with native c-fos activation were
missed in our assay. Finally, fth, c-fos is a member of a family
of IEGs regulated by neuronal activity, and the detection of other
IEGs (such as Arc, homer-1A, or zif-268) can be expected to
reveal partially overlapping activation maps compared to the
c-fos-GFP map identied in our paper. Because neuronal activation in some brain areas may induce other IEGs but fail to induce
c-fos, the c-fos-GFP-based network of brain areas described
here should not be interpreted as a complete brain activation
map evoked by social behavior.
CONCLUSIONS
Our method of c-fos-GFP-based screening generates cellularresolution maps of behaviorally evoked whole-brain activation
in the mouse. The patterns of female and male interactionevoked brain activation revealed clear separation between the
two stimuli, including at the level of brain structures downstream
of both volatile and nonvolatile chemosensory signaling. These
activation patterns were also markedly different from the activation pattern evoked during nonsocial olfactory-enhanced exploratory behavior. These ndings demonstrate that our method
can be used for screening behavior-evoked whole-brain activation, and we envision that future experiments will yield brainmap-like descriptions for other innate behaviors, such as
aggression and defensive behaviors, or cognitive behaviors,
such as attention and decision making. Further, the same
method can be applied to genetic mouse models of neurodevelopmental disorders with the aim of identifying circuit decits
underlying changes in social, cognitive, and other higher-order
brain functions.
were killed and perfused 1 hr later and the brains were xed overnight in 4%
paraformaldehyde, then cut as 50 mm coronal sections. For immunohistochemistry, sections were exposed to rabbit anti-c-fos antibody (1:10,000,
Santa Cruz SC052) and labeled by DAB solution. FIJI (ImageJ) and Volocity
(Perkin-Elmer) were used for cell counting.
Statistics
We ran statistical comparisons between different behavioral groups based on
either ROIs or evenly spaced voxels. Voxels were overlapping 3D spheres with
100 mm diameter each and spaced 20 mm apart from each other. The cell count
of each voxel was calculated as the number of nuclei found within 100 mm from
the center of the voxel in all 3D. To account for multiple comparisons across all
voxel/ROI locations, we thresholded the p values and reported FDRs. For correlation between c-fos-GFP cell counts and social behavior, Pearson correlation R values were calculated between c-fos-GFP cell counts and time spent in
social behaviors. See the Supplemental Experimental Procedures for more
details.
SUPPLEMENTAL INFORMATION
Supplemental Information includes Supplemental Experimental Procedures,
seven gures, three tables, and three movies and can be found with this article
online at http://dx.doi.org/10.1016/j.celrep.2014.12.014.
EXPERIMENTAL PROCEDURES
ACKNOWLEDGMENTS
Animals
Animal procedures were approved by the Cold Spring Harbor Laboratory Animal Care and Use Committee. The c-fos-GFP mice, Tg(Fos-tTA,Fos-EGFP)
line, were obtained from The Jackson Laboratory. In our study, we used the
direct c-fos-GFP signal, whereas several other studies used the tTA protein
to drive other reporter molecules (Garner et al., 2012; Liu et al., 2012; Matsuo
et al., 2008; Reijmers et al., 2007).
Cell Reports 10, 292305, January 13, 2015 2015 The Authors 303
Dong, H.W., and Swanson, L.W. (2004). Projections from bed nuclei of
the stria terminalis, posterior division: implications for cerebral hemisphere
regulation of defensive and reproductive behaviors. J. Comp. Neurol. 471,
396433.
Dong, H.W., and Swanson, L.W. (2006). Projections from bed nuclei of the stria
terminalis, magnocellular nucleus: implications for cerebral hemisphere regulation of micturition, defecation, and penile erection. J. Comp. Neurol. 494,
108141.
Dragunow, M., and Robertson, H.A. (1987). Kindling stimulation induces c-fos
protein(s) in granule cells of the rat dentate gyrus. Nature 329, 441442.
Fenno, L., Yizhar, O., and Deisseroth, K. (2011). The development and application of optogenetics. Annu. Rev. Neurosci. 34, 389412.
Ferguson, J.N., Young, L.J., Hearn, E.F., Matzuk, M.M., Insel, T.R., and Winslow, J.T. (2000). Social amnesia in mice lacking the oxytocin gene. Nat.
Genet. 25, 284288.
Ferguson, J.N., Aldag, J.M., Insel, T.R., and Young, L.J. (2001). Oxytocin in the
medial amygdala is essential for social recognition in the mouse. J. Neurosci.
21, 82788285.
Garner, A.R., Rowland, D.C., Hwang, S.Y., Baumgaertel, K., Roth, B.L.,
Kentros, C., and Mayford, M. (2012). Generation of a synthetic memory trace.
Science 335, 15131516.
Ghosh, S., Larson, S.D., Hefzi, H., Marnoy, Z., Cutforth, T., Dokka, K., and
Baldwin, K.K. (2011). Sensory maps in the olfactory cortex dened by longrange viral tracing of single neurons. Nature 472, 217220.
Grewe, B.F., and Helmchen, F. (2009). Optical probing of neuronal ensemble
activity. Curr. Opin. Neurobiol. 19, 520529.
Matsuo, N., Reijmers, L., and Mayford, M. (2008). Spine-type-specic recruitment of newly synthesized AMPA receptors with learning. Science 319, 1104
1107.
Mattes, D., Haynor, D.R., Vesselle, H., Lewellen, T.K., and Eubank, W. (2003).
PET-CT image registration in the chest using free-form deformations. IEEE
Trans. Med. Imaging 22, 120128.
Meyer, H.S., Wimmer, V.C., Oberlaender, M., de Kock, C.P., Sakmann, B., and
Helmstaedter, M. (2010). Number and laminar distribution of neurons in a thalamocortical projection column of rat vibrissal cortex. Cereb. Cortex 20, 2277
2286.
Morgan, J.I., Cohen, D.R., Hempstead, J.L., and Curran, T. (1987). Mapping
patterns of c-fos expression in the central nervous system after seizure. Science 237, 192197.
Olshausen, B.A., and Field, D.J. (2004). Sparse coding of sensory inputs. Curr.
Opin. Neurobiol. 14, 481487.
Pankevich, D.E., Baum, M.J., and Cherry, J.A. (2004). Olfactory sex discrimination persists, whereas the preference for urinary odorants from estrous females disappears in male mice after vomeronasal organ removal. J. Neurosci.
24, 94519457.
Pfaus, J.G., and Heeb, M.M. (1997). Implications of immediate-early gene induction in the brain following sexual stimulation of female and male rodents.
Brain Res. Bull. 44, 397407.
Pitkanen, A., Savander, V., and LeDoux, J.E. (1997). Organization of intraamygdaloid circuitries in the rat: an emerging framework for understanding
functions of the amygdala. Trends Neurosci. 20, 517523.
Guzowski, J.F., Timlin, J.A., Roysam, B., McNaughton, B.L., Worley, P.F., and
Barnes, C.A. (2005). Mapping behaviorally relevant neural circuits with immediate-early gene expression. Curr. Opin. Neurobiol. 15, 599606.
Ragan, T., Kadiri, L.R., Venkataraju, K.U., Bahlmann, K., Sutin, J., Taranda, J.,
Arganda-Carreras, I., Kim, Y., Seung, H.S., and Osten, P. (2012). Serial twophoton tomography for automated ex vivo mouse brain imaging. Nat. Methods
9, 255258.
Herculano-Houzel, S., Mota, B., and Lent, R. (2006). Cellular scaling rules for
rodent brains. Proc. Natl. Acad. Sci. USA 103, 1213812143.
Reijmers, L.G., Perkins, B.L., Matsuo, N., and Mayford, M. (2007). Localization
of a stable neural correlate of associative memory. Science 317, 12301233.
Hitti, F.L., and Siegelbaum, S.A. (2014). The hippocampal CA2 region is essential for social memory. Nature 508, 8892.
Saper, C.B., Lu, J., Chou, T.C., and Gooley, J. (2005). The hypothalamic integrator for circadian rhythms. Trends Neurosci. 28, 152157.
Sesack, S.R., and Grace, A.A. (2010). Cortico-Basal Ganglia reward network:
microcircuitry. Neuropsychopharmacology 35, 2747.
Keller, A., and Carlson, G.C. (1999). Neonatal whisker clipping alters intracortical, but not thalamocortical projections, in rat barrel cortex. J. Comp. Neurol.
412, 8394.
Keller, M., Douhard, Q., Baum, M.J., and Bakker, J. (2006). Sexual experience
does not compensate for the disruptive effects of zinc sulfatelesioning of the
main olfactory epithelium on sexual behavior in male mice. Chem. Senses 31,
753762.
Kelley, A.E., and Stinus, L. (1984). The distribution of the projection from the
parataenial nucleus of the thalamus to the nucleus accumbens in the rat: an
autoradiographic study. Exp. Brain Res. 54, 499512.
Sosulski, D.L., Bloom, M.L., Cutforth, T., Axel, R., and Datta, S.R. (2011).
Distinct representations of olfactory information in different cortical centres.
Nature 472, 213216.
Kopec, C.D., Bowers, A.C., Pai, S., and Brody, C.D. (2011). Semi-automated
atlas-based analysis of brain histological sections. J. Neurosci. Methods
196, 1219.
Lee, H.M., Giguere, P.M., and Roth, B.L. (2014). DREADDs: novel tools for
drug discovery and development. Drug Discov. Today 19, 469473.
Lin, D., Boyle, M.P., Dollar, P., Lee, H., Lein, E.S., Perona, P., and Anderson,
D.J. (2011). Functional identication of an aggression locus in the mouse hypothalamus. Nature 470, 221226.
Liu, X., Ramirez, S., Pang, P.T., Puryear, C.B., Govindarajan, A., Deisseroth,
K., and Tonegawa, S. (2012). Optogenetic stimulation of a hippocampal
engram activates fear memory recall. Nature 484, 381385.
Mandiyan, V.S., Coats, J.K., and Shah, N.M. (2005). Decits in sexual
and aggressive behaviors in Cnga2 mutant mice. Nat. Neurosci. 8,
16601662.
304 Cell Reports 10, 292305, January 13, 2015 2015 The Authors
Stowers, L., Holy, T.E., Meister, M., Dulac, C., and Koentges, G. (2002). Loss of
sex discrimination and male-male aggression in mice decient for TRP2. Science 295, 14931500.
Sunkin, S.M., Ng, L., Lau, C., Dolbeare, T., Gilbert, T.L., Thompson, C.L., Hawrylycz, M., and Dang, C. (2013). Allen Brain Atlas: an integrated spatio-temporal portal for exploring the central nervous system. Nucleic Acids Res. 41,
D996D1008.
Swanson, L.W. (2000). Cerebral hemisphere regulation of motivated behavior.
Brain Res. 886, 113164.
Taniguchi, H., He, M., Wu, P., Kim, S., Paik, R., Sugino, K., Kvitsiani, D., Fu, Y.,
Lu, J., Lin, Y., et al. (2011). A resource of Cre driver lines for genetic targeting of
GABAergic neurons in cerebral cortex. Neuron 71, 9951013.
Turaga, S.C., Murray, J.F., Jain, V., Roth, F., Helmstaedter, M., Briggman, K.,
Denk, W., and Seung, H.S. (2010). Convolutional networks can learn to
generate afnity graphs for image segmentation. Neural Comput. 22, 511538.
Veening, J.G., Coolen, L.M., de Jong, T.R., Joosten, H.W., de Boer, S.F.,
Koolhaas, J.M., and Olivier, B. (2005). Do similar neural systems subserve
aggressive and sexual behaviour in male rats? Insights from c-Fos and pharmacological studies. Eur. J. Pharmacol. 526, 226239.
Winslow, J.T. (2003). Mouse social recognition and preference. Curr. Protoc.
Neurosci. Chapter 8, 16.
Vousden, D.A., Epp, J., Okuno, H., Nieman, B.J., van Eede, M., Dazai, J., Ragan, T., Bito, H., Frankland, P.W., Lerch, J.P., and Henkelman, R.M. (2014).
Whole-brain mapping of behaviourally induced neural activation in mice. Brain
Struct. Funct.
Yang, C.F., Chiang, M.C., Gray, D.C., Prabhakaran, M., Alvarado, M., Juntti,
S.A., Unger, E.K., Wells, J.A., and Shah, N.M. (2013). Sexually dimorphic neurons in the ventromedial hypothalamus govern mating in both sexes and
aggression in males. Cell 153, 896909.
Yang, C.F., and Shah, N.M. (2014). Representing sex in the brain, one module
at a time. Neuron 82, 261278.
Cell Reports 10, 292305, January 13, 2015 2015 The Authors 305
Cell Reports
Resource
Insights into the Evolution of Longevity
from the Bowhead Whale Genome
Michael Keane,1,18 Jeremy Semeiks,2,18 Andrew E. Webb,3,18 Yang I. Li,4,18,19 Vctor Quesada,5,18 Thomas Craig,1
Lone Bruhn Madsen,6 Sipko van Dam,1 David Brawand,4 Patrcia I. Marques,5 Pawel Michalak,7 Lin Kang,7 Jong Bhak,8
Hyung-Soon Yim,9 Nick V. Grishin,2 Nynne Hjort Nielsen,10 Mads Peter Heide-Jrgensen,10 Elias M. Oziolor,11
Cole W. Matson,11 George M. Church,12 Gary W. Stuart,13 John C. Patton,14 J. Craig George,15 Robert Suydam,15
Knud Larsen,6 Carlos Lopez-Otn,5 Mary J. OConnell,3 John W. Bickham,16,17 Bo Thomsen,6
and Joao Pedro de Magalhaes1,*
1Integrative
Genomics of Ageing Group, Institute of Integrative Biology, University of Liverpool, Liverpool L69 7ZB, UK
Hughes Medical Institute and Departments of Biophysics and Biochemistry, University of Texas Southwestern Medical Center,
Dallas, TX 75390-9050, USA
3Bioinformatics and Molecular Evolution Group, School of Biotechnology, Dublin City University, Glasnevin, Dublin 9, Ireland
4MRC Functional Genomics Unit, University of Oxford, Oxford OX1 3QX, UK
5Departamento de Bioqumica y Biologa Molecular, Instituto Universitario de Oncologa (IUOPA), Universidad de Oviedo, 33006 Oviedo,
Spain
6Department of Molecular Biology and Genetics, Aarhus University, 8830 Tjele, Denmark
7Virginia Bioinformatics Institute, Virginia Tech, Blacksburg, VA 24061, USA
8Personal Genomics Institute, Genome Research Foundation, Suwon 443-270, Republic of Korea
9KIOST, Korea Institute of Ocean Science and Technology, Ansan 426744, Republic of Korea
10Greenland Institute of Natural Resources, 3900 Nuuk, Greenland
11Department of Environmental Science, Center for Reservoir and Aquatic Systems Research (CRASR) and Institute for Biomedical Studies,
Baylor University, Waco, TX 76798, USA
12Department of Genetics, Harvard Medical School, Boston, MA 02115, USA
13The Center for Genomic Advocacy (TCGA) and Department of Biology, Indiana State University, Terre Haute, IN 47809, USA
14Department of Forestry and Natural Resources, Purdue University, West Lafayette, IN 47907, USA
15North Slope Borough, Department of Wildlife Management, Barrow, AK 99723, USA
16Battelle Memorial Institute, Houston, TX 77079, USA
17Department of Wildlife and Fisheries Sciences, Texas A&M University, College Station, TX 77843, USA
18Co-rst author
19Present address: Department of Genetics, Stanford University, Stanford, CA 94305, USA
*Correspondence: [email protected]
http://dx.doi.org/10.1016/j.celrep.2014.12.008
This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/3.0/).
2Howard
SUMMARY
The bowhead whale (Balaena mysticetus) is estimated to live over 200 years and is possibly
the longest-living mammal. These animals should
possess protective molecular adaptations relevant
to age-related diseases, particularly cancer. Here,
we report the sequencing and comparative analysis
of the bowhead whale genome and two transcriptomes from different populations. Our analysis
identies genes under positive selection and bowhead-specic mutations in genes linked to cancer
and aging. In addition, we identify gene gain and
loss involving genes associated with DNA repair,
cell-cycle regulation, cancer, and aging. Our results
expand our understanding of the evolution of
mammalian longevity and suggest possible players
involved in adaptive genetic changes conferring cancer resistance. We also found potentially relevant
changes in genes related to additional processes,
112 Cell Reports 10, 112122, January 6, 2015 2015 The Authors
Libraries
200 bp paired-end
149.1
500 bp paired-end
141.7
48.73
3 kb mate-paired
57.3
19.73
5 kb mate-paired
72.5
24.93
10 kb mate-paired
28.5
9.83
449.1
154.33
Total
51.23
N50 (kb)
34.8
877
Number
113,673
2.1
7,227
2.3
this key resource for studying the bowhead whale and its various
traits, including its exceptional longevity and resistance to
diseases.
EXPERIMENTAL PROCEDURES
DNA and RNA Sampling in Greenland
Bowhead (Balaena mysticetus) DNA used for genome sequencing was isolated from muscle tissue sampled from a 51-year-old female (ID no. 325)
caught in the Disko Bay, West Greenland in 2009 (Heide-Jrgensen et al.,
2012). Tissue samples were stored at 20 C immediately after collection.
Age estimation was performed using the aspartic acid racemization technique
(Garde et al., 2007). CITES no. 12GL1003387 was used for transfer of biological material. Bowhead RNA used for RNA-seq and small RNA analysis
was isolated from two different individuals: kidney samples were from a
44-year-old female (ID no. 500) and muscle samples were isolated from a
44-year-old male (ID no. 322). For more details of the individual whales, see
Heide-Jrgensen et al. (2012).
Genome Sequencing
DNA was extracted following standard protocols, quantied using Qubit and
run on an agarose gel to ensure no degradation had occurred. We then generated 1503 coverage of the genome using the Illumina HiSeq 2000 platform
with 100 bp reads, sequencing paired-end libraries, and mate-paired libraries
with insert sizes of 3, 5, and 10 kb (Table 1). Sequencing was performed at
the Liverpool Centre for Genomic Research (CGR; http://www.liv.ac.uk/
genomic-research/).
Genome Assembly
Libraries were preprocessed in-house by the CGR to remove adaptor sequences. The raw fastq les were trimmed for the presence of the Illumina
adaptor sequence using Cutadapt and then subjected to window-based quality trimming using Sickle with a minimum window quality score of 20. A minimum read-length lter of 10 bp was also applied. Libraries were then assembled with ALLPATHS-LG (Gnerre et al., 2011), which performed all assembly
steps including read error correction, initial read alignment, and scaffolding.
ALLPATHS-LG build 43762 was used with the default input parameters,
including K = 96. Several build parameters were automatically determined by
the software at run time per its standard algorithm. Of 2.88 3 109 paired fragment reads and 1.87 3 109 paired jumping reads, 0.015% were removed as
poly(A) and 1.5% were removed due to low-frequency kmers; 54% of jumping
read pairs were error-corrected, and overall 33% of jumping pairs were redundant. In total, we used 216 Gbp for the 2.3 Gb assembly, meaning that coverage
retained for the assembly was 953. Full assembly and read usage data are
shown in Supplemental Folder 2. Assembly completeness was assayed with
CEGMA by searching for 248 core eukaryotic genes (Parra et al., 2007).
Genome Size Determination
To determine the genome size for bowhead whale, spleen tissues were
acquired from one male (10B17) and one female (10B18). Both whales were
harvested in 2010 as part of the native subsistence hunt in Barrow, Alaska.
Sample processing and staining followed the methods of Vindelv and Christensen (1994). Instrument description and additional methodological details
are provided in Oziolor et al. (2014). Briey, ow cytometric genome size determination is based on propidium iodide uorescent staining of nuclear DNA.
Mean uorescence is calculated for cells in the G0 and G1 phases of the cell
cycle. This method requires direct comparison to known standards to convert
measured uorescence to pg of DNA. The primary standard used in this study
was the domestic chicken (Gallus gallus domesticus). Chicken red blood cells
are widely used as a genome size standard, with an accepted genome size of
C = 1.25 pg. Chicken whole blood was purchased from Innovative Research.
Mouse (Mus musculus) and rat (Rattus norvegicus) were included as internal
checks, with estimates for both falling within 3% of previously published
genome size estimates (Vinogradov, 1998). Spleen tissues from three male
129/SvEvTac laboratory mice and a single male Harlan SD Sprague-Dawley
laboratory rat were used.
118 Cell Reports 10, 112122, January 6, 2015 2015 The Authors
Cell Reports 10, 112122, January 6, 2015 2015 The Authors 119
AUTHOR CONTRIBUTIONS
G.M.C., J.C.G., R.S., J.W.B., B.T., and J.P.M. conceived and designed the
study; L.B.M., E.M.O., and C.W.M. performed the experiments; M.K., J.S.,
A.E.W., Y.I.L., V.Q., L.B.M., S.v.D., D.B., P.I.M., P.M., L.K., J.B., H.-S.Y.,
G.W.S., J.C.P., C.L.-O., M.J.O., J.W.C., and J.P.M. analyzed the data; T.C.,
N.V.G., N.H.N., M.P.H.-J., R.S., and K.L. contributed reagents/materials/analysis tools; and M.K. and J.P.M. wrote the paper.
ACKNOWLEDGMENTS
This project was supported by grants from the Life Extension Foundation and
the Methuselah Foundation to J.P.M. We thank the Inuit whaling captains of
the Alaska Eskimo Whaling Commission and the Barrow Whaling Captains
Association for allowing us to sample their whales and their willingness to support the transcriptome study. This study was also partly funded by the Augustinus Foundation to K.L. T.C. was supported by a Wellcome Trust grant
(WT094386MA) to J.P.M. and M.K. was supported by a studentship from the
University of Liverpools Faculty of Health and Life Sciences. J.S. was supported by grants from the NIH (GM094575 to N.V.G.) and the Welch
Foundation (I-1505 to N.V.G.). Y.I.L. was supported by a Nufeld Department
of Medicine Prize studentship from the University of Oxford. C.L.-O. is an
Investigator of the Botin Foundation also supported by grants from Ministerio
de Economa y Competitividad-Spain and Instituto de Salud Carlos III (RTICC)Spain. M.J.O. and A.E.W. are funded by Science Foundation Ireland Research
Frontiers Programme Grant (EOB2763) to M.J.O. M.J.O. would also like to
acknowledge the Fulbright Commission for the Fulbright Scholar Award
2012-2013. M.J.O. and A.E.W. thank the SFI/HEA Irish Centre for High-End
Computing (ICHEC) for processor time and technical support. This work was
also supported by the Korea Institute of Ocean Science and Technology
(KIOST) in-house program (PE99212). Further thanks to Prof. Chris Ponting
and Dr. Andreas Heger for hosting gene orthology predictions from OPTIC,
the University of Liverpool High Performance Computing facilities for processor time, Eric de Sousa for help with RNA-seq data analysis, and to Louise
Crompton for assistance in compiling and formatting the bibliography. Last,
we are grateful to the staff at the Liverpool Centre for Genomic Research for
advice during this project.
Received: September 7, 2014
Revised: November 21, 2014
Accepted: December 3, 2014
Published: December 24, 2014
REFERENCES
Altschul, S.F., Gish, W., Miller, W., Myers, E.W., and Lipman, D.J. (1990). Basic
local alignment search tool. J. Mol. Biol. 215, 403410.
Austad, S.N. (2010). Methusalehs Zoo: how nature provides us with clues for
extending human health span. J. Comp. Pathol. 142 (Suppl 1), S10S21.
Burki, F., and Kaessmann, H. (2004). Birth and adaptive evolution of a
hominoid gene that supports high neurotransmitter ux. Nat. Genet. 36,
10611063.
Caulin, A.F., and Maley, C.C. (2011). Petos Paradox: evolutions prescription
for cancer prevention. Trends Ecol. Evol. 26, 175182.
Bovine Genome Sequencing and Analysis Consortium, Elsik, C.G., Tellam,
R.L., Worley, K.C., Gibbs, R.A., Muzny, D.M., Weinstock, G.M., Adelson,
D.L., Eichler, E.E., Elnitski, L., et al. (2009). The genome sequence of taurine
cattle: a window to ruminant biology and evolution. Science 324, 522528.
Cornu, M., Albert, V., and Hall, M.N. (2013). mTOR in aging, metabolism, and
cancer. Curr. Opin. Genet. Dev. 23, 5362.
SUPPLEMENTAL INFORMATION
de Magalhaes, J.P. (2013). How ageing processes inuence cancer. Nat. Rev.
Cancer 13, 357365.
120 Cell Reports 10, 112122, January 6, 2015 2015 The Authors
Dobin, A., Davis, C.A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., Batut,
P., Chaisson, M., and Gingeras, T.R. (2013). STAR: ultrafast universal RNA-seq
aligner. Bioinformatics 29, 1521.
Holt, C., and Yandell, M. (2011). MAKER2: an annotation pipeline and genomedatabase management tool for second-generation genome projects. BMC
Bioinformatics 12, 491.
Dupressoir, A., Vernochet, C., Bawa, O., Harper, F., Pierron, G., Opolon, P.,
and Heidmann, T. (2009). Syncytin-A knockout mice demonstrate the critical
role in placentation of a fusogenic, endogenous retrovirus-derived, envelope
gene. Proc. Natl. Acad. Sci. USA 106, 1212712132.
Huang, W., Sherman, B.T., and Lempicki, R.A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat.
Protoc. 4, 4457.
Finch, C. (1990). Longevity, Senescence, and the Genome (Chicago: University of Chicago Press).
Fischer, H., Koenig, U., Eckhart, L., and Tschachler, E. (2002). Human caspase
12 has acquired deleterious mutations. Biochem. Biophys. Res. Commun.
293, 722726.
Flicek, P., Ahmed, I., Amode, M.R., Barrell, D., Beal, K., Brent, S., CarvalhoSilva, D., Clapham, P., Coates, G., Fairley, S., et al. (2013). Ensembl 2013.
Nucleic Acids Res. 41, D48D55.
Friedlander, M.R., Chen, W., Adamidi, C., Maaskola, J., Einspanier, R., Knespel, S., and Rajewsky, N. (2008). Discovering microRNAs from deep
sequencing data using miRDeep. Nat. Biotechnol. 26, 407415.
Kim, E.B., Fang, X., Fushan, A.A., Huang, Z., Lobanov, A.V., Han, L., Marino,
S.M., Sun, X., Turanov, A.A., Yang, P., et al. (2011). Genome sequencing reveals insights into physiology and longevity of the naked mole rat. Nature
479, 223227.
Friedlander, M.R., Mackowiak, S.D., Li, N., Chen, W., and Rajewsky, N. (2012).
miRDeep2 accurately identies known and hundreds of novel microRNA
genes in seven animal clades. Nucleic Acids Res. 40, 3752.
Kurihara, L.J., Kikuchi, T., Wada, K., and Tilghman, S.M. (2001). Loss of UchL1 and Uch-L3 leads to neurodegeneration, posterior paralysis and dysphagia.
Hum. Mol. Genet. 10, 19631970.
Garde, E., Heide-Jorgensen, M.P., Hansen, S.H., Nachman, G., and Forchhammer, M.C. (2007). Age-specic growth and remarkable longevity in narwhals (Monodon monoceros) from West Greenland as estimated by aspartic
acid racemization. J. Mammal. 88, 4958.
Lee, D.Y., Hayes, J.J., Pruss, D., and Wolffe, A.P. (1993). A positive role for histone acetylation in transcription factor access to nucleosomal DNA. Cell 72,
7384.
Gatesy, J., Geisler, J.H., Chang, J., Buell, C., Berta, A., Meredith, R.W.,
Springer, M.S., and McGowen, M.R. (2013). A phylogenetic blueprint for a
modern whale. Mol. Phylogenet. Evol. 66, 479506.
George, J.C., Bada, J., Zeh, J., Scott, L., Brown, S.E., OHara, T., and Suydam,
R. (1999). Age and growth estimates of bowhead whales (Balaena mysticetus)
via aspartic acid racemization. Can. J. Zool. 77, 571580.
Gillet, L.C.J., and Scharer, O.D. (2006). Molecular mechanisms of mammalian
global genome nucleotide excision repair. Chem. Rev. 106, 253276.
Gnerre, S., Maccallum, I., Przybylski, D., Ribeiro, F.J., Burton, J.N., Walker,
B.J., Sharpe, T., Hall, G., Shea, T.P., Sykes, S., et al. (2011). High-quality draft
assemblies of mammalian genomes from massively parallel sequence data.
Proc. Natl. Acad. Sci. USA 108, 15131518.
Goodstadt, L., and Ponting, C.P. (2006). Phylogenetic reconstruction of orthology, paralogy, and conserved synteny for dog and human. PLoS Comput. Biol.
2, e133.
Gori, F., Friedman, L.G., and Demay, M.B. (2006). Wdr5, a WD-40 protein, regulates osteoblast differentiation during embryonic bone development. Dev.
Biol. 295, 498506.
Grabherr, M.G., Haas, B.J., Yassour, M., Levin, J.Z., Thompson, D.A., Amit, I.,
Adiconis, X., Fan, L., Raychowdhury, R., Zeng, Q., et al. (2011). Full-length
transcriptome assembly from RNA-seq data without a reference genome.
Nat. Biotechnol. 29, 644652.
Li, Y., and de Magalhaes, J.P. (2013). Accelerated protein evolution analysis
reveals genes and pathways associated with the evolution of mammalian
longevity. Age (Dordr.) 35, 301314.
Lohse, M., Bolger, A.M., Nagel, A., Fernie, A.R., Lunn, J.E., Stitt, M., and Usadel, B. (2012). RobiNA: a user-friendly, integrated software solution for RNAseq-based transcriptomics. Nucleic Acids Res. 40, W622W627.
Lopez-Otn, C., Blasco, M.A., Partridge, L., Serrano, M., and Kroemer, G.
(2013). The hallmarks of aging. Cell 153, 11941217.
Loytynoja, A., and Goldman, N. (2008). A model of evolution and structure for
multiple sequence alignment. Philos. Trans. R. Soc. Lond. B Biol. Sci. 363,
39133919.
Morgan, C.C., Foster, P.G., Webb, A.E., Pisani, D., McInerney, J.O., and
OConnell, M.J. (2013). Heterogeneous models place the root of the placental
mammal phylogeny. Mol. Biol. Evol. 30, 21452156.
Oziolor, E.M., Bigorgne, E., Aguilar, L., Usenko, S., and Matson, C.W. (2014).
Evolved resistance to PCB- and PAH-induced cardiac teratogenesis, and
reduced CYP1A activity in Gulf killish (Fundulus grandis) populations from
the Houston Ship Channel, Texas. Aquat. Toxicol. 150, 210219.
Parra, G., Bradnam, K., and Korf, I. (2007). CEGMA: a pipeline to accurately
annotate core genes in eukaryotic genomes. Bioinformatics 23, 10611067.
Philo, L.M., Shotts, E.B., and George, J.C. (1993). Morbidity and mortality. In
The Bowhead Whale, J.J. Burns, J.J. Montague, and C.J. Cowles, eds. (Lawrence, Kansas: Allen Press), pp. 275312.
Hedges, S.B., Dudley, J., and Kumar, S. (2006). TimeTree: a public knowledgebase of divergence times among organisms. Bioinformatics 22, 29712972.
Quesada, V., Ordonez, G.R., Sanchez, L.M., Puente, X.S., and Lopez-Otn, C.
(2009). The Degradome database: mammalian proteases and diseases of proteolysis. Nucleic Acids Res. 37, D239D243.
Heger, A., and Ponting, C.P. (2007). Evolutionary rate analyses of orthologs
and paralogs from 12 Drosophila genomes. Genome Res. 17, 18371849.
Heide-Jrgensen, M.P., Garde, E., Nielsen, N.H., Andersen, O.N., and Hansen, S.H. (2012). A note on biological data from the hunt of bowhead whales
in West Greenland 2009-2011. J. Cetacean Res. Manag. 12, 329333.
Ross, P.L., Cheng, I., Liu, X., Cicek, M.S., Carroll, P.R., Casey, G., and Witte,
J.S. (2009). Carboxypeptidase 4 gene variants and early-onset intermediateto-high risk prostate cancer. BMC Cancer 9, 69.
Hoege, C., Pfander, B., Moldovan, G.L., Pyrowolakis, G., and Jentsch, S.
(2002). RAD6-dependent DNA repair is linked to modication of PCNA by ubiquitin and SUMO. Nature 419, 135141.
Seim, I., Ma, S., Zhou, X., Gerashchenko, M.V., Lee, S.G., Suydam, R.,
George, J.C., Bickham, J.W., and Gladyshev, V.N. (2014). The transcriptome
of the bowhead whale Balaena mysticetus reveals adaptations of the
longest-lived mammal. Aging (Albany, N.Y. Online) 6, 879899.
Tacutu, R., Craig, T., Budovsky, A., Wuttke, D., Lehmann, G., Taranukha, D.,
Costa, J., Fraifeld, V.E., and de Magalhaes, J.P. (2013). Human Ageing
Cell Reports 10, 112122, January 6, 2015 2015 The Authors 121
Genomic Resources: integrated databases and tools for the biology and genetics of ageing. Nucleic Acids Res. 41, D1027D1033.
Vinogradov, A.E. (1998). Genome size and GC-percent in vertebrates as determined by ow cytometry: the triangular relationship. Cytometry 31, 100109.
Tanno, M., Ogihara, M., and Taguchi, T. (1996). Age-related changes in proliferating cell nuclear antigen levels. Mech. Ageing Dev. 92, 5366.
Weeda, G., Donker, I., de Wit, J., Morreau, H., Janssens, R., Vissers, C.J.,
Nigg, A., van Steeg, H., Bootsma, D., and Hoeijmakers, J.H.J. (1997). Disruption of mouse ERCC1 results in a novel repair syndrome with growth failure,
nuclear abnormalities and senescence. Curr. Biol. 7, 427439.
Tervo, O.M., Christoffersen, M.F., Parks, S.E., Kristensen, R.M., and Madsen,
P.T. (2011). Evidence for simultaneous sound production in the bowhead
whale (Balaena mysticetus). J. Acoust. Soc. Am. 130, 22572262.
Tian, X., Azpurua, J., Hine, C., Vaidya, A., Myakishev-Rempel, M., Ablaeva, J.,
Mao, Z., Nevo, E., Gorbunova, V., and Seluanov, A. (2013). High-molecularmass hyaluronan mediates the cancer resistance of the naked mole rat. Nature
499, 346349.
West, G.B., Woodruff, W.H., and Brown, J.H. (2002). Allometric scaling of
metabolic rate from molecules and mitochondria to cells and mammals.
Proc. Natl. Acad. Sci. USA 99 (Suppl 1), 24732478.
Yang, Z. (2007). PAML 4: phylogenetic analysis by maximum likelihood. Mol.
Biol. Evol. 24, 15861591.
van Beekum, O., Gao, Y., Berger, R., Koppen, A., and Kalkhoven, E. (2012). A
novel RNAi lethality rescue screen to identify regulators of adipogenesis. PLoS
ONE 7, e37680.
Yang, Z., and dos Reis, M. (2011). Statistical properties of the branch-site test
of positive selection. Mol. Biol. Evol. 28, 12171228.
van Dam, S., Craig, T., and de Magalhaes, J.P. (2015). GeneFriends: a human
RNA-seq-based gene and transcript co-expression database. Nucleic Acids
Res. http://dx.doi.org/10.1093/nar/gku1042
Yang, Z., Wong, W.S.W., and Nielsen, R. (2005). Bayes empirical bayes inference of amino acid sites under positive selection. Mol. Biol. Evol. 22, 1107
1118.
Vervoort, V.S., Viljoen, D., Smart, R., Suthers, G., DuPont, B.R., Abbott, A., and
Schwartz, C.E. (2002). Sorting nexin 3 (SNX3) is disrupted in a patient with a
translocation t(6;13)(q21;q12) and microcephaly, microphthalmia, ectrodactyly, prognathism (MMEP) phenotype. J. Med. Genet. 39, 893899.
Yim, H.S., Cho, Y.S., Guang, X., Kang, S.G., Jeong, J.Y., Cha, S.S., Oh, H.M.,
Lee, J.H., Yang, E.C., Kwon, K.K., et al. (2014). Minke whale genome and
aquatic adaptation in cetaceans. Nat. Genet. 46, 8892.
Vilella, A.J., Severin, J., Ureta-Vidal, A., Heng, L., Durbin, R., and Birney, E.
(2009). EnsemblCompara GeneTrees: Complete, duplication-aware phylogenetic trees in vertebrates. Genome Res. 19, 327335.
Zhang, J., Zhang, Y.P., and Rosenberg, H.F. (2002). Adaptive evolution of a
duplicated pancreatic ribonuclease gene in a leaf-eating monkey. Nat. Genet.
30, 411415.
Vindelv, L.L., and Christensen, I.J. (1994). Detergent and proteolytic enzymebased techniques for nuclear isolation and DNA content analysis. Methods
Cell Biol. 41, 219229.
Zhang, J., Nielsen, R., and Yang, Z. (2005). Evaluation of an improved branchsite likelihood method for detecting positive selection at the molecular level.
Mol. Biol. Evol. 22, 24722479.
122 Cell Reports 10, 112122, January 6, 2015 2015 The Authors
933>
G=C@
47<53@
=<
B63
>C:A3
EWbV1SZZ>`Saa@SdWSea
:]]YW\UT]`OcbV]`WbObWdS`SdWSea]\bVS
T]`ST`]\b]TaQWS\QS-Bc`\b]1SZZ>`Saa
=c`SRWb]`aabOgOP`SOab]TbVSZObSab
RSdSZ]^[S\baa]bVSgQO\Q][[WaaW]\
Sf^S`bQcbbW\USRUS`SdWSea
4W\Rg]c`eOg
eWbV1SZZ>`Saa@SdWSea
eeeQSZZQ][`SdWSea
B]^`]^SZg]c``SaSO`QVT]`eO`RTOabS`bc`\b]
1SZZ>`Saa@SdWSea=c`W\aWUVbTcZOcbV]`WbObWdS
`SdWSea^cPZWaVSROQ`]aabVSZWTSaQWS\QSaW\
]c`^`W[O`g`SaSO`QVO\RB`S\RaX]c`\OZaU]
PSg]\Rag\bVSaWaO\R]TTS`O^]W\b]TdWSe
Traumatic Brain Injury: Clinical, Pathological & Translational Mechanisms (J3) 0RGHUQ3KHQRW\SLF'UXJ'LVFRYHU\'HQLQJWKH3DWK)RUZDUG'
joint with Axons: From Cell Biology to Pathology (J4)
Apr 26, 2016 | Big Sky, Montana | USA
Jan 2427, 2016 | Santa Fe, New Mexico | USA
1HXURORJLFDO'LVRUGHUVRI,QWUDFHOOXODU7UDIFNLQJ$
The Cancer Genome (Q1) joint with Genomics & Personalized Medicine (Q2)
New Approaches to Vaccines for Human & Veterinary Tropical Diseases (M1)
1HZ)URQWLHUVLQ8QGHUVWDQGLQJ7XPRU0HWDEROLVP4
joint with ,PPXQRPHWDEROLVPLQ,PPXQH)XQFWLRQ ,QDPPDWRU\'LVHDVH4
2016
SCIENTIFIC
CONFERENCES
Presenting the most signicant research on cancer etiology, prevention, diagnosis, and treatment
*1.5*%!#+
The Function of Tumor Microenvironment in
Cancer Progression
Conference Co-Chairpersons: Raghu Kalluri, Robert A.
Weinberg, Douglas Hanahan, and Morag Park
*1.5*%!#+
Patient-Derived Cancer Models: Present and Future
Applications from Basic Science to the Clinic
Conference Co-Chairpersons: Manuel Hidalgo,
Hans Clevers, S. Gail Eckhardt, and Joan Seoane
!.1.5!3.(!*/
10th AACR-JCA Joint Conference on
Breakthroughs in Cancer Research:
From Biology to Therapeutics
Co-Chairpersons: Frank McCormick and Tetsuo Noda
!.1.5
1%
AACR Precision Medicine Series:
Cancer Cell Cycle-Tumor Progression
and Therapeutic Response
Conference Chairperson: Julien Sage
Co-Chairpersons: Karen E. Knudsen and
J. Alan Diehl
!.1.5
.$.(* +
AACR Annual Meeting 2016
Program Committee Chairperson: Scott A. Armstrong
,.%(!3.(!*/
Pancreatic Cancer
Conference Co-Chairpersons: Manuel Hidalgo,
Christine Iacobuzio-Donahue, and
Robert H. Vonderheide
5.(* +
AACR Precision Medicine Series:
Targeting the Vulnerabilities of Cancer
Conference Co-Chairpersons: Stephen W. Fesik,
Jeffrey E. Settleman, and Paul Workman
5
%)%
Engineering and Physical Sciences in Oncology
Conference Co-Chairpersons: Rakesh Jain,
Robert Langer, and Joan Brugge
1*!+/0+*
Translational Control of Cancer:
A New Frontier in Cancer Biology and Therapy
Conference Co-Chairpersons: Jennifer A. Doudna,
Frank McCormick, Davide Ruggero, and
Nahum Sonenberg
0+!.*.*%/+
DNA Repair: Tumor Development and
Therapeutic Response
Conference Co-Chairpersons: Robert G. Bristow,
Theodore S. Lawrence, and Maria Jasin
+2!)!.
+*0.!(1!!*
EORTC-NCI-AACR Molecular Targets and
Cancer Therapeutics Symposium
+2!)!.!!)!.
1*%$!.)*5