Gpr176 is a Gz-linked orphan G-protein-coupled

Full text

Turn on search term navigation

ARTICLE

Received 15 Sep 2015 | Accepted 30 Dec 2015 | Published 17 Feb 2016

Masao Doi1,2,*, Iori Murai1,*, Sumihiro Kunisue1,*, Genzui Setsu1, Naohiro Uchio1, Rina Tanaka1,Sakurako Kobayashi1, Hiroyuki Shimatani1, Hida Hayashi1, Hsu-Wen Chao1,2, Yuuki Nakagawa1, Yukari Takahashi1, Yunhong Hotta1, Jun-ichirou Yasunaga3, Masao Matsuoka3, Michael H. Hastings4, Hiroshi Kiyonari5& Hitoshi Okamura1,2

G-protein-coupled receptors (GPCRs) participate in a broad range of physiological functions. A priority for fundamental and clinical research, therefore, is to decipher the function of over 140 remaining orphan GPCRs. The suprachiasmatic nucleus (SCN), the brains circadian pacemaker, governs daily rhythms in behaviour and physiology. Here we launch the SCN orphan GPCR project to (i) search for murine orphan GPCRs with enriched expression in the SCN, (ii) generate mutant animals decient in candidate GPCRs, and (iii) analyse the impact on circadian rhythms. We thereby identify Gpr176 as an SCN-enriched orphan GPCR that sets the pace of circadian behaviour. Gpr176 is expressed in a circadian manner by SCN neurons, and molecular characterization reveals that it represses cAMP signalling in an agonist-independent manner. Gpr176 acts independently of, and in parallel to, the Vipr2 GPCR, not through the canonical Gi, but via the unique G-protein subclass Gz.

DOI: 10.1038/ncomms10583 OPEN

Gpr176 is a Gz-linked orphan G-protein-coupled receptor that sets the pace of circadian behaviour

1 Department of Systems Biology, Graduate School of Pharmaceutical Sciences, Kyoto University, Sakyo-ku, Kyoto 606-8501, Japan. 2 Core Research for Evolutional Science and Technology, Japan Science and Technology Agency, Saitama 332-0012, Japan. 3 Laboratory of Virus Control, Institute for Virus Research, Kyoto University, Sakyo-ku, Kyoto 606-8507, Japan. 4 Division of Neurobiology, Medical Research Council Laboratory of Molecular Biology, Cambridge CB2 0QH, UK. 5 Animal Resource Development Unit and Genetic Engineering Team, RIKEN Center for Life Science Technologies, Kobe 650-0047, Japan. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed to M.D. (email: mailto:[email protected]

Web End [email protected] ) or to H.O. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 7:10583 | DOI: 10.1038/ncomms10583 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10583

a positive regulator of cAMP22,23 and demonstrated to be necessary for SCN time-keeping24,25. Yet, much less is known about the molecular identity of GPCR that negatively regulates cAMP production in the SCN.cAMP synthesis is positively or negatively regulated by Gs or Gi family members, respectively. While the Gs family contains two subtypes (Gs1 and Gs2), the inhibitory members include three Gi (Gi1, Gi2, and Gi3) and one Gz. All Gi members, except Gz, are substrates of pertussis toxin (PTX). PTX mediates ADP ribosylation at the carboxyl terminal cysteine residue ( 4 position), inhibiting Gi activity. Because Gz lacks this

cysteine residue, it can inhibit adenylyl cyclase activity in a PTX-insensitive manner26. Differing from Gi, Gz is expressed mainly in the brain27,28, but its roles in the brain are not understood.

In the present study we surveyed all known orphan GPCRs expressed in the SCN, and identied three SCN-enriched genes: Gpr176, Gpr19 and Calcr. We generated knockout mice for each of them and demonstrated that Gpr176 is a unique orphan GPCR that can set the pace of circadian behaviour. Gpr176 is expressed mainly in the brain, with prominent expression in the SCN, and its protein abundance uctuates in a circadian fashion. Molecular characterization further revealed that this orphan receptor has an agonist-independent basal activity to repress cAMP production.Notably, the unique G-protein subclass Gz, but not the canonical Gi, is required for the activity of Gpr176. We show that Gpr176 (negative regulator of cAMP) acts independently of, and additively with, the Vipr2.

ResultsThe SCN-orphan GPCR project identies Gpr176. We rst constructed a list of GPCR genes that are potentially enriched in the SCN (Fig. 1a). Using SCN microarray data29 we surveyed all known mouse non-odorant GPCR genes and rank-ordered them based on their relative levels. The list of the top 100 genes (Fig. 1a) identied 23 orphan GPCRs, whose in vivo functions remain unknown. The receptors with known key functions in the SCN are also included in this list, for example, Adcyap1r1 (refs 30,31), Prokr2 (refs 32,33), Avpr1a (refs 34,35) and Vipr2 (refs 24,25), suggesting the relevance of the screen.

The genes highlighted in blue in Fig. 1a were further analysed by in situ hybridization. These include all listed orphan GPCRs (Gpr37l1, Lphn1, Gpr37, Gprc5b, Gpr85, Darc, Gpr176, Lphn3, Lphn2, Gpr48, Gpr56, Gpr19, Gpr123, Gpr83, Ackr3, Gpr125, Gpr22, Gpr153, Gpr45, Gpr68, Gpr6, Cckar and Gpr75) and three selected GPCRs whose roles in the SCN are not known (Ntsr2, Oprl1, and Calcr). To specify their histological distribution, we

Figure 1 | Gpr176 is an SCN-enriched orphan GPCR that sets the pace of circadian timing. (a) The SCN-GPCR project leading to the identication of Gpr176. The bar graph shows the rank order of expression of the top 100 GPCRs (classes A, B, and C) detected in the SCN microarray analysis (GEO accession number: GSE28574). Red bars are the receptors known to be expressed in the SCN. The horizontal line indicates the threshold of statistical signicance of expression. The genes highlighted in blue, which include all listed orphan GPCRs, were characterized further by in situ hybridization using radiolabeled gene-specic probes. Arrows indicate robust positive SCN signals for Gpr176, Calcr and Gpr19. (b) Snake-plot representation of the mouse Gpr176. The residues highlighted in yellow indicate the DRYxxV motif located at the cytoplasmic end of the transmembrane helix III. (c) Northern blotting for Gpr176 with a mouse multiple-tissue blot (Clontech). (d) Representative mouse coronal brain section immunolabeled for Gpr176. Scale bar, 1 mm. (e) Schematic representations of the mouse Gpr176 gene, targeting construct, and the resulting mutant allele. A genomic region downstream of the start codon (ATG) of Gpr176 was deleted. Grey box: probe used for Southern blot. (f) Southern blot of KpnI-digested DNA from Gpr176 / , Gpr176 /

and Gpr176 / mice. Genomic fragments from the wild-type (7.8 kb) and mutant (12.3 kb) alleles are indicated. (g) Immunohistochemical conrmation of Gpr176 deciency in the SCN of Gpr176 / mice. Scale bar, 100 mm. oc, optic chiasm; v, third ventricle. (h) Representative locomotor activity records of

C57BL/6J-backcrossed Gpr176 / and Gpr176 / mice. Mice were housed in LD and then transferred to DD. Periods of darkness are indicated by grey backgrounds. Data are shown in double-plotted format. Each horizontal line represents 48 h; the second 24-h period is plotted to the right and below the rst. (i) Circadian periods of free-running activities in DD. Periods of individual mice are plotted. Bars indicate means.e.m. (n 12, for each genotype).

***Po0.001, Students t-test. (j) Representative Per1-luc bioluminescence records from organotypic SCN slices of Gpr176 / (light and dark blue traces) and Gpr176 / (red and orange traces) mice. (k) Periods of Per1-luc rhythm in Gpr176/ and Gpr176 / SCN slices (meanss.e.m., n 4 for each).

**Po0.01, Students t-test.

In mammals, the principal circadian pacemaker governing daily rhythms in behaviour and physiology resides in the suprachiasmatic nucleus (SCN) of the hypothalamus1,2.

Most tissues outside the SCN also contain local clocks (the so-called peripheral clocks), and their rhythms are synchronized, harmoniously, by an array of direct or indirect signals from the SCN2,3. Thus, the SCN lies at the top of a hierarchical, multioscillator system distributed across the body4,5.

Malfunction of the circadian clock has been linked to the pathogenesis of a wide variety of diseases6, including sleep-wake disorder, tumorigenesis, obesity, diabetes and hypertension. Drug efcacy and toxicity are also under circadian regulation7. These lines of evidence support the potential value of developing drugs that target the circadian clock, and pioneer studies have already identied synthetic compounds that selectively target the key intracellular clock components, cryptochromes (Cry1 and Cry2)8 and REV-ERBa and b9. Because their targets are distributed across the body, these compounds can modulate both the central and peripheral clocks equally8,9. In contrast, the development of drugs that specically target the SCN remains an unfullled opportunity for circadian pharmacology4,7.

G-protein-coupled receptors (GPCRs) constitute the largest family of cell surface receptors, participating in a broad range of physiological functions. It has been appreciated that GPCRs are the most common target of pharmaceutical drugs: more than 30% of clinically marketed drugs target GPCR function10. Intriguingly, there are still 4140 orphan GPCRs whose cognate ligands are not known, and deciphering their physiological function remains a priority for both clinical and fundamental research1115. We speculated that some orphan GPCRs, whose physiological functions have remained unknown, might exist in the SCN and function as potential modulators of the circadian system.

Structurally, GPCRs possess two different conformations, active and inactive. Agonists lock the receptor structure in its active form, antagonists block agonist action, and inverse agonists stabilize the receptor in its inactive form. In the absence of ligands, GPCRs spontaneously interchange between the two conformations; active and inactive, generating agonist-independent baseline activity1618. Although the magnitude of this spontaneous activity differs strikingly between GPCRs, some of the orphan GPCRs exhibit signicant levels of intrinsic activity19,20.

Depending on the type of G-protein to which the GPCR is coupled, a variety of downstream signalling pathways can be activated. Circadian uctuation of cAMP signal is crucial for the maintenance of circadian clock function in the SCN21. In this context, the Vipr2 GPCR for vasoactive intestinal peptide (Vip) is

2 NATURE COMMUNICATIONS | 7:10583 | DOI: 10.1038/ncomms10583 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10583 ARTICLE

labelled mouse brain sections with gene-specic probes, and relatively strong, SCN-specic signals were observed for Gpr176, Gpr19 and Calcr (Fig. 1a). On the other hand, the remaining genes were either only faintly expressed in the SCN or located more broadly in the regions outside of the SCN. Thus, based on

this intensity and regionality, we selected Gpr176, Gpr19 and Calcr as candidates of interest.

We generated knockout mice for Gpr176, Gpr19 and Calcr (Fig. 1e,f, see also Methods). Homozygous deletion of each of them did not cause any gross abnormalities, lethality or infertility,

Gpr176

PACAP receptor

SCN-GPCR project

DNA microarray SCN-GPCR Top100

Log2 GPCR mRNA levels

Prokineticin receptor

SCN microarray

In situ hybridization

Creation of KO mice

Behavioural test

Identification of a new relevant SCN-GPCR

AVP receptor

VIP/PACAP receptor

Ntsr2 Gpr37l1 Oprl1

Lphn1 Gpr37 Gprc5b Darc

Gpr85 Gpr176 Lphn3 Calcr Gpr48

Lphn2

Gpr56 Gpr19 Gpr83

Gpr123 Ackr3 Gpr125 Gpr153

Gpr22 Gpr45 Gpr68 Cckar Gpr75

Gpr6

b c d

Gpr176

Heart

Brain

Spleen

Lung

Liver

Skeletal muscle

Kidney

Testis

Extracellular

9.5

7.5

4.4

2.4

Intracellular

IHC: anti-Gpr176

e f

Gpr176

7.6 kb

Wild-type (+) allele

Mutant () allele

Probe

KpnI

ATG

12.3 kb

+/+

Targeting construct

Exon1

Tau-LacZ

loxP

Neo

7.6 kb

loxP

12.3 kb

Tau-LacZ

loxP

Neo

loxP

KpnI

g h i

Gpr176

Time (h)

Gpr176

Time (h)

0 0

+/+ 12 0

24.5

1 7

21 28 35

***

Free running period (h)

Gpr176 +/+

Days

1 7

21 28 35

24.0

23.5

Gpr176 /

Days

1 7

21 28 35

Days

1 7

21 28 35

0 0

/ 12 0

23.0

Gpr176 +/+

Gpr176/

j k

Gpr176 +/+ Gpr176 +/+

Gpr176 / Gpr176 /

Per1-luc SCN slice culture

25.5

130

25.0

Bioluminescence

(106 photon counts per 20 min)

125

Period (h)

120

24.5

115

110

24.0

105

23.5

100 0 12 24 36 48 60 72 84 108 120

Time in culture (h)

Gpr176+/+

Gpr176/

NATURE COMMUNICATIONS | 7:10583 | DOI: 10.1038/ncomms10583 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10583

allowing for circadian locomotor activity tests on respective adult mutant mice. Animals were entrained on a 12-h light:12-h dark (LD) cycle for 2 weeks and then transferred to constant darkness (DD) to monitor the endogenous SCN-driven locomotor activity rhythm. We found that whereas Gpr176-decient mice exhibited a signicantly shorter circadian period compared with wild-type littermate mice, mice lacking Gpr19 or Calcr did not show any signicant difference in period by genotype (circadian periods(h), means.e.m. (n 6 for each genotype); for Gpr176 / ,

23.780.04; Gpr176 / , 23.360.04 (Po0.01); for Gpr19 / ,23.790.04; Gpr19 / , 23.920.07 (P40.1); for Calcr / ,23.840.03; Calcr / , 23.860.04 (P40.1); Students t-test). Importantly, Gpr176 mutant mice that had been backcrossed to the C57BL/6J background over 10 generations also had a similarly short-period phenotype (Fig. 1h; period (w2 periodogram)s.e.m:

Gpr176 / , 23.390.03 h; Gpr176 / , 23.770.02 h (n 12,

both genotypes); Fig. 1i, Po0.001, Students t-test). On the other hand, as shown in Supplementary Fig. 1, Gpr176-decient mice were normal for light-induced phase resetting; their resetting response to early- and late-night light pulses were indistinguishable from wild-type mice. Thus, Gpr176 is not likely to be a light signal-related receptor for the SCN. Rather, this orphan receptor appears to be involved in the determination of intrinsic period of the SCN.

To complement the data on circadian behaviour, we performed real-time bioluminescence imaging of cultured SCN slices from wild-type and Gpr176 / mice carrying a Per1-promoter-luciferase (Per1-luc) reporter gene36 (Fig. 1j,k).

Crucially, bioluminescence rhythm of the Gpr176 / SCN slice had a shorter period than that of the wild-type control (Gpr176/ 25.110.14 h; Gpr176 / 24.390.16 h;

means.e.m., n 4 for each, Po0.01, t-test); the short-period

phenotype of behaviour is thus ascribable to the SCN. In contrast, the short period was not observed for their lung slices (Gpr176 / 24.760.09 h; Gpr176 / 24.810.31 h; n 4

for each), suggesting that the effect of Gpr176 deciency may be specic to the central oscillator.

Of note, northern blot analysis with different tissues (Fig. 1c; Supplementary Fig. 2a,b) revealed that Gpr176 messenger RNA (mRNA) is expressed chiey in the brain. In situ hybridization was therefore further conducted to map the neural expression from the forebrain to the medulla oblongata (Supplementary Fig. 2c). Notably, the highest expression of Gpr176 was observed within the SCN. Relatively strong and restricted expression was also seen in the subfornical organ (SFO), organum vasculosum of the lamina terminalis (OVLT) and cerebellar occulus (Fl), among other regions, suggesting that Gpr176 in extra-SCN sites may contribute to different brain functions that we have not yet systematically investigated.

Clock-controlled circadian expression of Gpr176 in the SCN. Gpr176 encodes a class A orphan GPCR37(Fig. 1b), whose in vivo function is unexplored. To further investigate its role in the circadian system, we proled temporal expression of Gpr176 in the SCN (Supplementary Fig. 2d). To do this, we performed quantitative in situ hybridization, using mice housed either in LD or DD, their brains being collected at 4-h intervals across 24-h cycles. Notably, in both LD and DD, the SCN had a time-of-dayspecic expression of Gpr176 (Supplementary Fig. 2d). In LD, Gpr176 was highest at night at ZT16 (ZT represents Zeitgeber time; ZT0 denotes lights-on and ZT12 lights-off) and lowest in the early morning at ZT0 (Po0.01, peak versus trough, one-way analysis of variance (ANOVA) with Bonferroni post hoc test). Similarly, in DD, Gpr176 was highest in the subjective night at CT16 (CT represents circadian time; CT0 denotes the beginning of the subjective day and CT12 the beginning of the subjective

night) and lowest in the subjective morning at CT4 (Po0.01, peak versus trough). Furthermore, we found that the subjective-night peak expression of Gpr176 was severely damped in mice decient in the core clock components Cry1 and Cry2, which completely lack a functional circadian clock38,39 (Supplementary Fig. 2e), indicating that the SCN clockwork controls Gpr176 gene expression. A light pulse given at night did not cause any acute change in Gpr176 expression in the SCN (Supplementary Fig. 2f). Thus the SCN appears to direct Gpr176 expression independently of light.

Characteristic distribution of Gpr176 across the whole SCN. We developed an afnity-puried rabbit polyclonal antibody against Gpr176 and performed immunohistochemistry. Coronal brain sections from wild-type mice (Fig. 1d) revealed marked SCN-specic immunoreactivity to Gpr176 (see also Supplementary Fig. 3a), while little or no immunostaining was observed for Gpr176 / SCN (Fig. 1g; Supplementary Fig. 3a).

The SCN is composed of anatomically heterogeneous sub-regions1,40. We thus studied the topographical distribution of Gpr176 by staining serial coronal brain sections covering the whole rostral-caudal extent of the SCN (Supplementary Fig. 3c). Positive immunostaining of Gpr176 was widely observed from the rostral to caudal extremities of the SCN. Moreover, we noticed that the immunoreactivity to Gpr176 was relatively strong in the dorsomedial area of the SCN, a region also referred to as the SCN shell41. This characteristic distribution prole was reminiscent of the Vipr2 receptor for Vip42. Interestingly, anti-Vipr2 immunoreactivity in the SCN, revealed by rabbit polyclonal antiserum, was widespread, extending from the rostral to caudal margins, and more intense in the dorsomedial than the ventrolateral region of the SCN (Supplementary Fig. 3b,c)42.

Immunolocalization of Gpr176 was also compared with that of neuropeptide markers in the SCN. We performed double-label confocal immunohistochemistry (Supplementary Fig. 4) and found that the two non-overlapping SCN populations, Vip- and vasopressin (AVP)-ergic neurons, both appear to express Gpr176. Interestingly, Vipr2 is also expressed in both populations42.

Histological relationship between Gpr176 and Vipr2. The resemblance of immunohistological distribution of Gpr176 to that of Vipr2 attracted our attention. Vipr2 is a class C GPCR that serves as a Vip receptor in the SCN and has already been shown to play a key role in circadian pacemaking24,25. We thus claried whether Gpr176 colocalizes with this important receptor.

To visualize Vipr2 expression in the SCN, we generated anti-Vipr2 chicken polyclonal antibody, for which we conrmed specic SCN staining and its absence in mice lacking Vipr2 (Supplementary Fig. 5). Dual-label immunouorescence microscopy (Fig. 2a) revealed markedly overlapping immunoreactivities for Gpr176 (red) and Vipr2 (green). Both signals tended to be more intense in the dorsomedial than ventrolateral area of the SCN (Fig. 2a). Moreover, high-magnication images showed that both receptors were expressed in almost all individual cells with comparable subcellular locations (Fig. 2b). Dual-labelled cells were also observed when the SCN neurons were dispersed in culture (Fig. 2c), conrming the colocalization of Gpr176 and Vipr2.

Circadian expression proles of Gpr176 and Vipr2 in the SCN were then compared. Brains were collected from mice at six time points in DD. Because of widespread distribution of Gpr176 and Vipr2, expression levels of each receptor were evaluated as a sum of the whole immunoreactivity from the rostral to caudal extremities of the SCN (10 sections per brain; n 6 mice for each

data point). The results revealed that Gpr176 immunoreactivity

4 NATURE COMMUNICATIONS | 7:10583 | DOI: 10.1038/ncomms10583 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10583 ARTICLE

a e

Vip2r /

Gpr176 Vipr2 Merge

Gpr176 +/+

Gpr176 /

0 0

12 0

0 0

12 0

0 0

12 0

0 0

12 0

14 21 28 35 42 49 56

Days Amplitude Amplitude

14 21 28 35 42 49 56

Gpr176

Vipr2

Gpr176 Vipr2

1,500

1,000

500

800 600 400 200

1,000

800 600 400 200

Gpr176 Vipr2 DAPI

16 18 20 22 24 26 28 30 32

Period (h)

16 18 20 22 24 26 28 30 32

Period (h)

DAPI

Period (h)

16 18 20 22 24 26 28 30 32

Period (h)

f g

Vip2r /

Vip2r +/+

Gpr176 /

Gpr176 +/+ Gpr176 +/+ Gpr176 /

12 0

DAPI Vipr2 Gpr176

0 0

12 0

0 0

12 0

0 0

12 0

P.C. Vipr2

Days

14 21 28 35 42 49 56

Gpr176

0 4 8 12 16 20

Gpr176

Vipr2

1,000

800 600 400 200

600

400

2,000 1,500 1,000

1,500

1,000

500

16 18 20 22 24 26 28 30 32

Period (h)

200 16 18 20 22 24 26 28 30 32

Period (h)

500 16 18 20 22 24 26 28 30 32

Period (h)

16 18 20 22 24 26 28 30 32

Period (h)

100

Relative Gpr176 intensity

Relative Vipr2 intensity

***

100

***

Gpr176 Vipr2

Per cent

Free running period (h)

24.0

23.5

23.0

22.5

22.0

0 4 8

Circadian time

12 16 20

+/+

Gpr176

Vipr2

0 +/+

+/+

Gpr176

Vipr2

+/+

/ /

Figure 2 | Characterization of histological and genetic relationships between Gpr176 and Vipr2. (a) Double-label immunouorescence of Gpr176 and Vipr2 in the mouse SCN. Coronal brain sections were immunolabeled with antibodies against Gpr176 (rabbit polyclonal) and Vipr2 (chicken polyclonal). Representative confocal pictures are shown, with a merged image of Gpr176 (red) and Vipr2 (green). Scale bar, 100 mm. (b) Enlargement of the boxed area in (a). Merge shows combined images for Gpr176 (red), Vipr2 (green), and DAPI-based nuclear staining (blue). Scale bar, 20 mm. (c) Immunouorescence and phase contrast (P.C.) images of dispersed SCN neurons. Cells were immunolabelled for Gpr176 and Vipr2. Merge shows combined images of Gpr176 (red), Vipr2 (green), and DAPI (blue). Scale bar, 20 mm. (d) Antiphasic circadian expression proles of Gpr176 and Vipr2 in the SCN. Values (means.e.m.) indicate relative immunoreactivities of Gpr176 and Vipr2 at 6 time points in DD (n 6 brains for each data point). Representative images of

the immunolabeled SCN sections are shown on the top. Scale bar, 200 mm. (eg) Representative double-plotted actograms and w2 periodograms of locomotor activity rhythms of mice carrying wild-type, Gpr176 / and Vipr2 / alleles. Vipr2 / mice exhibit a single (e) or multiple (f) circadian periods in DD, while all wild-type and Gpr176 / mice had a single, stable circadian period (g). Diagonal line on periodogram shows signicance at

Po0.001. See also Supplementary Fig. 6. (h) Percentage of mice expressing a single dominant circadian period (black) or multiple circadian periods (grey). n 8 12 for each genotype. (i) Group data showing individual and combined effects of Gpr176 / and Vipr2 / on circadian period of locomotor

activity rhythms. Values (means.e.m.) indicate free-running periods of mice with a single dominant period in DD (n 5 8). ***Po0.001, one-way

ANOVA with Bonferroni post hoc test.

was highest in the subjective night at CT16 and lowest in the subjective day at CT4 (Po0.01, peak versus trough, one-way

ANOVA with Bonferroni post hoc test; Fig. 2d). This protein cycle is almost in phase with its mRNA expression that peaks at CT16 (Supplementary Fig. 2d,e). In comparison, Vipr2

immunoreactivity displayed an opposite circadian cycle, characterized by a robust decrease in the subjective night at CT16 (Fig. 2d) as reported previously42. These results illustrate that circadian proles of Gpr176 and Vipr2 in the SCN are antiphasic.

NATURE COMMUNICATIONS | 7:10583 | DOI: 10.1038/ncomms10583 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10583

Genetic relationship between Gpr176 and Vipr2. We next sought to understand the relationship between Gpr176 and Vipr2 (Fig. 2ei). To do this, we generated double decient mice (Gpr176 / ; Vipr2 / ). By crossing double heterozygous mice (Gpr176 / ;Vipr2 / ), we obtained genetic background (C57BL/6J)-matched male wild-type (n 8), respective single

(n 8 for Gpr176 / , n 11 for Vipr2 / ), and double

knockout (n 12) mice, and compared their behavioural rhythms

in DD. As observed above, all single Gpr176 knockout mice (Gpr176 / ; Vipr2 / ) had shorter periods than those of control mice (Gpr176 /;Vipr2 / ) (Fig. 2i; Supplementary

Fig. 6a,b). Vipr2 / mice had altered circadian rhythms as documented previously25,43 (Fig. 2eg), displaying either greatly reduced or multiple circadian periods in DD. About half of Vipr2 / mice (6 of 11) had a single, short-circadian period (Fig. 2e), whereas the other half (5 of 11) simultaneously

expressed two or more statistically signicant circadian periods (Fig. 2f; Actograms of all individual mice are available in Supplementary Fig. 6). Concomitant deletion of Gpr176 and Vipr2 (Gpr176 / ;Vipr2 / ) did not change this mixed phenotype (Fig. 2e,f): about half of doubly decient mice (6 of 12) still had a single short-circadian period (Fig. 2h), indicating that Gpr176 is not involved in the penetrance of this phenotype. Under these conditions, we compared period length between genotypes (Fig. 2i) and found that the deletion of Gpr176 could induce further signicant shortening of the circadian period length in Vipr2 / background: the circadian period of double knockout mice (Gpr176 / ;Vipr2 / ) in DD was22.200.03 h (w2 periodogram, means.e.m.), which was signicantly shorter than that of Gpr176 wild-type, Vipr2 knockout mice at 22.630.05 h (Po0.001, one-way ANOVA with Bonferroni post hoc test). Thus, the shortening effect of

a b

Dox ()

Dox (+)

kDa

d e

Flp-In TREx293 Gpr176 (tet-on)/ Vipr2 (constitutive)

100

Dox ()

Dox (+)

Dox ()

Dox (+)

Dox ()

Dox (+)

Vip

Gpr176

Vipr2

75 75

cAMP (pmol per 105cells)

Specific 125 I-Vip binding (c.p.m. 103 )

Vipr2

[cAMP]i

Gpr176

Tet-On

Gpr176

***

Dox (+)

Forskolin (Fsk)

***

Gpr176

DAPI

Dox ()

***

Vipr2

Vehicle

100 nM 10 nM 10 M

1 nM

Fsk

100 nM

10 nM

Vip

125I-Vip

1 nM

f g

Gpr176

(tet-on)

Gpr176RDY (tet-on)

Gpr176V145R (tet-on)

2.5

2.0

Human ** * *

DRY RDY V R

+ +

Mouse

Chicken

Xenopus

Zebrafish

137-AIALDRYYSVLYPLERKIS-155

136-AIALDRYYSVLYPLERKIS-154

154-AIALDRYYSVLYPLERKIT-172

147-AIALDRYYSVLYPLERKIS-165

156-VIALDRYYSVLYPLERKIS-174

Relative cAMP levels

1.5

***

Gpr176 (tet-on):

DRYxxV motif

RDY V145R

1.0

Dox + +

***

Anti-Gpr176

Anti-Vipr2

Anti-Gi

0.5

Vehicle Vip Fsk

PTX-treated Gpr176

(tet-on)

Vehicle Vip Fsk

Figure 3 | Gpr176 basal activity inhibits VipVipr2cAMP signalling. (a) Experimental schemes based on Flp-In TREx293-Gpr176(tet-on)/Vipr2 cells. (b) Immunoblots of Gpr176 and Vipr2 with cell membrane fractions of Flp-In TREx293-Gpr176(tet-on)/Vipr2 cells. Cells were treated with either 1 mg ml 1

Dox or control vehicle for 24 h. Gi was used as a loading control. (c) Confocal images of representative dual-label immunouorescent staining of Gpr176 and Vipr2 in Dox-treated or non-treated Flp-In TREx293-Gpr176(tet-on)/Vipr2 cells. (d) Antagonistic basal activity of Gpr176 on Vip- and forskolin (Fsk)-stimulated cAMP accumulation in Flp-In TREx293-Gpr176(tet-on)/Vipr2 cells. Cells of the same batch were cultured in parallel with or without Dox for 24 h then resuspended in 0.5 mM IBMX (non-selective phosphodiesterase inhibitor)-containing assay buffer for 1 h and stimulated with Vip or Fsk for 15 min. Individual treatments were run in each experiment in triplicate. cAMP levels were determined by enzyme immunoassays. Values are expressed as means.e.m. per 1 105 cells. ***Po0.001 in comparison with the control stimulated group. (e) Unaltered, dose-dependent [125I]-Vip binding activities

of Dox-treated Flp-In TREx293-GPR176(tet-on)/Vipr2 cells. Values are meanss.e.m. of four independent experiments. (f) Immunoblots of mutant cell lines expressing either Gpr176RDY or Gpr176V145R. Point mutations were introduced into the conserved DRYxxV motif residues on Gpr176 in Flp-In TREx293-Gpr176(tet-on)/Vipr2 cells. Western blotting was done as described in b. (g) Diminished basal activities of Gpr176 by DRYxxV motif mutation. Dox-treated or non-treated cells were stimulated with Vip (100 nM), Fsk (10 mM), or vehicle then subjected to cAMP assays. Where specied, cells were treated with PTX (100 ng ml 1) for 16 h before assay. Each experiment was done in triplicate. cAMP values (meanss.e.m.) are plotted relative to those of

Fsk-stimulated uninduced cells. ***Po0.001, **Po0.01 in comparison with the control stimulated group.

6 NATURE COMMUNICATIONS | 7:10583 | DOI: 10.1038/ncomms10583 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10583 ARTICLE

Gpr176 / was not masked by, and acted additively with, the absence of Vipr2.

Gpr176 basal activity antagonizes VipVipr2cAMP signalling. To explore potential functional interaction between Gpr176 and Vipr2 we developed a heterologous expression system (Fig. 3a). Agonist-independent constitutive function has been assigned to a

number of orphan GPCRs. Thus, in the absence of a known ligand, we looked for a constitutive action or inuence of Gpr176 on VipVipr2-mediated signalling (Fig. 3a). We employed the Flp-In TREx293 cell system and established a clonal cell line that expressed Vipr2 constitutively (Fig. 3b,c, middle) and Gpr176 in a doxycycline-dependent, inducible manner (Fig. 3b,c, upper; see Methods for details). Confocal microscopy

Gnaz probe

Anti-sense Sense #2V

siRNA

DN-Gz

kDa

37 kDa

()

(+)

NT NC #1

Gz (DN-Gz)

RGSZ1

Mock Infection

Dox () Dox (+)

Dox () Dox (+) Dox () Dox (+)

1.5

NC siRNA

Gz siRNA#1

Gz siRNA#2

1.5

Mock

RGSZ1

1.5

DN-Gz () DN-Gz (+)

NS NS

Gnaz

Gnai1

Gnai2

Gnai3

Absolute mRNA levels (mmol mol1)

Relative cAMP levels

0.5

***

Relative cAMP levels

0 1 2 3 4 5 30

25 35 40

1.0

***

Relative cAMP levels

1.0

***

0.5

Flp-In TREx293

NIH3T3

Vehicle Fsk Vehicle Fsk Vehicle Fsk Vehicle Fsk Vehicle Fsk

Vehicle Fsk Vehicle Fsk

g h

3T3 Tet-on 3G Gpr176

SCN punches

TREx293

3T3 Tet-on 3G Gpr176

Gz()

Gz(+)

kDa

Dox () Dox (+) Dox () Dox (+)

Gpr176 +/+

***

1.5

1.0

0.5

1.5

1.0

0.5

Gpr176 /

Gz (+)

Gz ()

PTX (+)

Tublin

Gpr176

Dox

Relative cAMP levels

cAMP concentrations

(pmol mg1 SCN protein)

1,000

800

600

400

200

100

***

+ +

***

*** **

Vehicle Fsk

100 M

Vehicle Fsk

100 M

Fsk 10 M

Vehicle Fsk

100 M

Fsk 10 M

CT4 CT16

Figure 4 | Gpr176 couples to Gz. (a) In situ hybridization of the gene encoding Gz (Gnaz) in the mouse SCN using a digoxigenin-labelled gene-specic probe. Scale bar, 100 mm. (b) Expressions of the genes encoding Gz (Gnaz), Gi1 (Gnai1), Gi2 (Gnai2) and Gi3 (Gnai3) in TREx293 (red) and NIH3T3 (blue)

cells. Values (meanss.e.m.) were determined by qRTPCR and normalized to those of the ribosomal phosphoprotein P0 (Rplp0)-encoding gene (n 3).

(c) Interference of Gpr176-mediated cAMP repression by knocking down Gz protein with two different siRNA mixtures. For simplicity, Flp-In TREx293 cells bearing only Gpr176 were used (no concomitant Vipr2 expression in this assay). Cells were transfected with siRNAs for 96 h, with or without Dox for the nal 24 h, and stimulated by Fsk (100 mM). Transfection was also performed with a negative-control siRNA (NC). Immunoblots show efcient knockdown of Gz. There was no off-target effect on Gi. NT, non-transfected cells. cAMP levels are plotted relative to those of Fsk-stimulated uninduced cells. Data are the means.e.m. of at least two independent experiments with three replicates each. ***Po0.001 in comparison with the control stimulated group. NS, not signicant. (d) Inactivation of Gpr176-mediated cAMP reduction by RGSZ1 coexpression. Cells were infected with lentiviruses encoding GFP (mock) or HA-tagged RGSZ1 and treated with or without Dox for 24 h. Western blot shows delivered proteins (anti-HA). cAMP values in Fsk-stimulated (100 mM) and nonstimulated (vehicle) cells are shown as described in c. (e) Interruption of Gpr176-mediated cAMP reduction by coexpression of DN-Gz. DN-Gz was stably introduced into the cells, and Dox-treated or non-treated cells were stimulated with Fsk (100 mM) or vehicle. Western blot (anti-Gz) shows the predominance of DN-Gz over the endogenous protein. cAMP values are expressed as described in c. (f) Generation and western blot conrmation of

NIH3T3-based Tet-inducible Gpr176 cell lines with or without stable expression of Gz. Cells were treated with or without Dox and immunoblotted for indicated proteins. Note that Gz is scant in parental NIH3T3 cells. (g) cAMP repressive activities expressed by Gpr176 in NIH3T3 Gz( ) cells. Dox-treated

or non-treated cells were stimulated with Fsk or vehicle. Where specied, cells were treated with PTX for 16 h before assay. cAMP values are expressed as described in c. (h) cAMP levels in the SCN of Gpr176 / and Gpr176 / mice at CT4 and CT16 in DD. SCN punches taken from mice were subjected to cAMP quantication by using a cAMP-specic enzyme immunoassay kit. Each data point represents a single animal. Bars indicate meanss.e.m. for

Gpr176/ CT4 (n 11), Gpr176 / CT4 (n 12), Gpr176 / CT16 (n 11), and Gpr176 / CT16 (n 11). ***Po0.001, two-way ANOVA with

Bonferroni post hoc test.

NATURE COMMUNICATIONS | 7:10583 | DOI: 10.1038/ncomms10583 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10583

revealed markedly increased plasma membrane-localized immunouorescence of Gpr176 in induced cells. On the other hand, constant membrane localization was observed for Vipr2.

Induction of Gpr176 led to an attenuation of VipVipr2-mediated cAMP signalling (Fig. 3d): Vip-promoted cAMP accumulation was signicantly reduced by doxycycline treatment (Po0.01, versus nontreatment control) at any of the dosages of

Vip used for stimulation (100, 10, and 1 nM). This was not due to a difference in cell viability or cell number between doxycyclineand control-treated groups, the growth rate of the cells being equivalent in the two groups. Furthermore, we assayed the same number of cells (Fig. 3d) resuspending 3 105 viable cells in

serum-free assay buffer for stimulation with Vip.

Neither plasma membrane localization nor total protein amount of Vipr2 changed appreciably with doxycycline treatment (Fig. 3b,c). Moreover, 125I-labelled Vip binding was almost identical between doxycycline-treated and non-treated cells (Fig. 3e). Thus, attenuation of Vip signalling by Gpr176 induction is likely to arise from changes in the downstream pathway. In accordance with this notion, forskolin-stimulated cAMP accumulation was also blunted by the induction of Gpr176 (Fig. 3d). Forskolin is a cell-permeable drug that directly activates adenylyl cyclases. Thus, the mode of action of Gpr176 does not necessarily depend on the integrity of VipVipr2. Consistent with this, Gpr176 exhibited similar basal activity when expressed alone in Flp-In TREx293 cells without concomitant expression of Vipr2 (Supplementary Fig. 7ac).

Contrasting to the overt effect on cAMP, induction of Gpr176 did not bring about any signicant change in inositol phosphate IP1 formation (Supplementary Fig. 7d). Thus, the downstream action of Gpr176 activity appears to be linked specically to cAMP regulation.

Gpr176 contains a conserved AspArgTyrXXVal (DRYxxV) motif at the cytoplasmic end of transmembrane domain III (Figs 1b and 3f). Generally, this motif is important for coupling of GPCRs to the partner G-proteins. We thus wondered whether this motif is required for generation of the basal activity of Gpr176. In other GPCRs, including the related orphan receptor Gpr161 (ref. 19) as well as the other class A non-orphan GPCRs such as Cxcr1 (ref. 44) and Adra1b (ref. 45) a single point mutation of the valine residue located following the DRY triplet sequence led to a drastic loss of the receptor-mediated signalling. We therefore established a mutant Gpr176 inducible cell line where the homologous valine of this protein was mutated to arginine (V145R; Fig. 3f). Notably, regardless of the type of stimuli (Vip or forskolin), mutant Gpr176 did not yield any noticeable reduction of cAMP (Fig. 3g), although the induction of protein levels was similar to that of the wild-type (Fig. 3f), suggesting that the V145R mutation blocks Gpr176 activity. As an alternative mutation, the DRY sequence was changed to RDY (Fig. 3f)46, but this mutation caused a loss of detectable Gpr176 expression and as a result no associated cAMP reduction was observed (Fig. 3f,g). Thus, the DRYxxV motif seems important for both protein activity and stability of Gpr176. We generated additional mutants on the DRY triplet sequence and tested their activities using a cAMP GloSensor assay, and the results further support this conclusion (see Supplementary Fig. 8).

Pertussis toxin (PTX) has been widely used to test the involvement of Gi signalling. This toxin is able to inactivate all Gi/o family members except Gz. Importantly, PTX displayed a strong inhibitory effect on Gi-coupled sphingosine-1-phosphate (S1P) receptor signalling in TREx293 cells47 (Supplementary Fig. 9). However, we could not detect any noticeable effect of PTX on the Gpr176-mediated signalling (Fig. 3g; Supplementary Fig. 7c): Gpr176 still had an essentially unimpaired capacity to reduce cAMP accumulation even after PTX treatment

(100 ng ml 1, 16 h). Thus, PTX-insensitive G-protein might mediate the downstream action of Gpr176.

Gpr176 couples to Gz. Gz is a unique Gi/o subfamily member that can repress adenylyl cyclases in a PTX-insensitive manner26, and its expression is known to be high particularly in the brain27,28 and in several specic tissues or cells in the periphery48. Microarray data indicate that the gene encoding Gz (Gnaz) is also expressed in the human embryonic kidney (HEK)293 cells49, a parental cell line of Flp-In TREx293 cells. We performed in situ hybridization and conrmed that the gene encoding Gz (Gnaz) is expressed in the mouse SCN (Fig. 4a; Supplementary Fig. 10). Gnaz was also found to be expressed in Flp-In TREx293 cells with abundance comparable to that of the other Gi family members (Fig. 4b). By contrast, the mouse embryonic broblast NIH3T3 cells did not express Gz (Fig. 4b,f) despite displaying high-level expression of various Gi. Akin to Gpr176, Gz is conserved among vertebrates (see http://ensembl.org/Multi/GeneTree/

Web End =http://ensembl.org/Multi/GeneTree/ ).

To test the hypothesis that Gz might be a mediator of Gpr176 signalling, we performed siRNA-mediated knockdown of the endogenous Gz protein in Flp-In TREx293 cells (Fig. 4c), using two different siRNA mixtures, both of which reduced the endogenous Gz protein expression levels to o10% of those observed for control siRNA treatment (NC) or nontransfection (NT) control (Fig. 4c). There was no off-target effect on the levels of Gi (Fig. 4c). Notably, both siRNAs against Gz abrogated Gpr176 activity, while negative control (NC) siRNA did not (Fig. 4c). We also used the regulator of G-protein signalling (RGS) protein family member RGSZ1, which selectively inhibits Gz50,51. We observed that lentiviral expression of this protein prevented the suppression of cAMP levels by Gpr176 (Fig. 4d). Similarly, overexpression of a dominant negative Gz (DN-Gz) protein52 resulted in a loss of the Gpr176-mediated cAMP reduction (Fig. 4e). Furthermore, to test this activity in a different cell system, we developed doxycycline-inducible NIH3T3 cell lines expressing Gpr176 (Fig. 4f). Because of the deciency of Gz in NIH3T3 cells (Fig. 4b), Gz was stably introduced into the cells at levels comparable to those of Flp-In TREx293 cells (Fig. 4f, Gz( )). Under these conditions Gpr176 displayed its effect

(Fig. 4g): a signicant reduction of forskolin-stimulated cAMP accumulation was observed when the cells were treated with doxycycline. In agreement with the unique property of Gz, PTX did not inhibit this reduction (Fig. 4g). Moreover, without exogenously expressed Gz, the cells did not elicit any noticeable activities of Gpr176, conrming that Gpr176 requires Gz for its activity (Fig. 4g, Gz( )).

Gpr176 is a negative modulator of cAMP synthesis in the SCN. The data thus far in cell culture suggest that Gpr176 is a negative regulator of cAMP signalling. Finally, we examined cAMP content in the SCN from wild-type and Gpr176 / mice (Fig. 4h).

In agreement with the above in vitro data, we observed that the deletion of Gpr176 leads to increased cAMP content in the SCN (Fig. 4h): The effect of the deletion of Gpr176 is more evident and signicant at CT16 than CT4 (Po0.001 for CT16, two-way

ANOVA with Bonferroni post hoc test), while even at CT4 there is a similar trend towards increased cAMP levels in Gpr176 /

SCN. Based on these data, we conclude that Gpr176 can act as a cAMP suppressor not only in vitro but also in vivo in the SCN.

DiscussionIn the hope of identifying a new GPCR that tunes the central clock, we searched for orphan GPCRs whose expression is enriched in the SCN. Gene knockout studies of candidate genes of interest revealed that Gpr176 is an SCN-enriched orphan GPCR

8 NATURE COMMUNICATIONS | 7:10583 | DOI: 10.1038/ncomms10583 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10583 ARTICLE

required for normal circadian behaviour. Nearly all SCN neurons express Gpr176, and its abundance uctuates in a circadian manner. Molecular characterization revealed that this orphan receptor has an agonist-independent basal activity to repress cAMP production. Notably, the unique G-protein subclass Gz, but not the canonical Gi, is required for the activity of Gpr176.

In the SCN, Gpr176 colocalizes with Vipr2. Given this overlapping expression, together with the ability of Gpr176 to compete with the VipVipr2cAMP signal, we surmise that the nighttime cAMP repression mediated by Gpr176 may serve as a part of the cAMP repressing mechanism that could counteract the Vip-Vipr2 axis in the SCN. As an additional feature of importance, Gpr176 is expressed mainly in the brain, with prominent expression in the SCN. This region specicity suggests that, as a putative target for circadian therapeutics, Gpr176 might possess the advantage of specicity, compared with broadly expressed neurotransmitter receptors, clock gene products and general modulators of second-messenger signalling such as phosphodiesterases. Thus, the results not only reveal a new signalling module, Gpr176/Gz, in the control of SCN circadian time-keeping, but also provide, thereby, a new class of GPCR signalling as a potential drug target to modulate the central clock in the brain.

In the present study, we focused on genes enriched in the SCN. The rationale for this tissue-specic approach is, rst, the SCN is responsible for circadian behaviour, the ultimate and denitive arbiter of daily life53. Second, we considered that specicity for local mechanisms is required because the SCN differs from the peripheral clocks in phasing. In rodents, for example, the molecular clockwork of the SCN phase-leads that of peripheral clocks by 7 to 11 h (ref. 54). Any drug targeting common circadian mechanisms would be able to affect or reset all the clocks simultaneously, and this might be deleterious for keeping the adaptive phasic order between the tissues (unless such drugs were to be delivered to the target tissue selectivelya very demanding regimen). We thus reasoned that pursuing the molecular mechanisms that underpin the specicity of the SCN would be a valuable alternative way to search for potential drug targets in the central clock. Third, and last, the SCN-GENE selection strategy can be biased for druggable targets. In the current study we found Gpr176 by focusing on the GPCR family, but the project readily lends itself to other target categories, for example, ion channels, suggesting the continued relevance of this screening method, in conjunction with others, for the discovery of clock gene modulators in the SCN.

Agonist-independent basal activity has been observed for a number of physiologically relevant orphan GPCRs: for example, Gpr3, an orphan GPCR endowed with constitutive Gs signalling activity, protects oocytes from ageing20 and modulates amyloid-b production in neurons55. The orphan receptor Gpr161, which participates in Sonic Hedgehog (Shh) signalling during neural tube development, also displays constitutive activity19. Agonist-independent intrinsic activity is also implicated in non-orphan GPCRs physiological function. The odorant receptors (ORs), which are GPCRs, possess an intrinsic activity andin the absence of activating odorantregulate axonal projection of olfactory neurons56. Thus, by analogous inference, intrinsic activity of Gpr176 could underlie its effect in vivo. Indeed, there may be no endogenous ligand, in which case the main regulatory control of activity would be via the level of protein expression, which for Gpr176 was highly circadian in the SCN. Consistent with its negative effect on cAMP, Gpr176 protein increased at night, a phase when cAMP levels in the SCN were decreased to the circadian nadir levels21,29. The absence of Gpr176, in turn, led to reduced suppression of cAMP level in the night, supporting the hypothesis. Nevertheless, the presence of

unidentied endogenous ligands is always difcult to exclude for receptors with constitutive action. In this respect, it is worth noting that even without a known natural ligand, surrogate ligands can be developed for the orphan GPCRs11,17, highlighting the druggable feature of this protein family.

Previous studies demonstrated that pharmacological inhibition of cAMP synthesis in the SCN leads to longer circadian period of behavioural rhythm21,29. Compatible with this, the loss of the cAMP suppressor Gpr176 leads to the opposite phenotype, period shortening. In agreement with the reduced suppression of cAMP signal during the night, we observed that circadian rising phase of Per1-luc activity was signicantly accelerated in the Gpr176 /

SCN slices (see Supplementary Fig. 11; waveform analysis of Fig. 1j), implying that the deletion of the suppressive signal from Gpr176 allows early rising of Per1 expression and thereby shortens circadian period. A simple explanation for the possible underlying mechanism may involve cAMP signal-mediated regulation of Per1 transcription through a cAMP-responsive element on its promoter57. However, the mechanism(s) through which the circadian uctuation of cAMP signal is integrated to the core clock machinery is still unclear in the literature and needs further exploration. In addition, a cohort of double decient mice for Gpr176 and Vipr2 (Gpr176 / ; Vipr2 / )

were still rhythmic in DD, albeit with a severely shortened period of circadian locomotor activity rhythms. These data raise the possibility of potential compensatory mechanisms, perhaps through alternative GPCRs in the SCN, as previously suggested58. A complete understanding of circadian regulation of cAMP signalling in the SCN, and of the concerted roles of Gpr176 (repressor) and Vipr2 (activator) will be a challenge of future study.

Unlike Gi, Gz does not serve as a substrate for PTX because the consensus cysteine residue in the fourth position from the C terminus of Gi proteins is replaced with isoleucine. We showed that Gpr176 links preferentially to Gz. It is generally considered that Gz shares the same receptor coupling prole with the Gi subtypes, but there are some exceptions. For example, Smoothened (smo), which is an orphan G-protein-linked seven-transmembrane protein that mediates Hedgehog signalling, has been shown to activate Gz more drastically than Gi59. The mechanism of specicity remains unclear but may involve the different amino-acid sequence of the C-terminal region between Gi and Gz. Because the C-terminal region of G proteins is important for physical interactions with upstream receptors60, slightly higher hydrophobicity of isoleucine (Gz) over cysteine (Gi) may affect this selectivity. Finally, as a critical difference between Gz and Gi, Gz is expressed predominantly in the brain. The brain distribution of Gz is more widespread than Gpr176, implying that additional Gz-linked orphan GPCRs may remain to be identied in the brain.

Gpr176 is an evolutionally conserved, vertebrate class A orphan GPCR, initially cloned by Hata et al.37 from a human brain cDNA library. Amino-acid sequence analysis reveals that Gpr176 contains four putative glycosylation sites at the N-terminal part37. A relatively large, C-terminal domain of about 200 amino acids (the exact sequence is omitted in the snake plot presentation in Fig. 1b) also characterizes Gpr176. Of interest, this C-terminal cytosolic region is highly conserved among Gpr176 genes in different species, yet does not show homology to any other annotated protein sequences, thus implying a unique feature characterizing this GPCR. Within its seven transmembrane domains, Gpr176 does not have an aspartic acid at position 2.50 (BW numbering), a feature also rarely observed for class A GPCRs. Understanding how these structural features affect molecular functions of Gpr176 will be a topic of future study.

NATURE COMMUNICATIONS | 7:10583 | DOI: 10.1038/ncomms10583 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10583

In summary, we described here that Gpr176 is an SCN-enriched orphan GPCR that can set the pace of circadian behaviour. This is the rst report to assign a function of this orphan receptor in physiology. We revealed that Gpr176 is a previously uncharacterized Gz-linked orphan GPCR that bears intrinsic activity to reduce cAMP production. The discovery of the functional orphan GPCR with a novel mode of action within the SCN would be of help to understand the mechanism that underpins the SCN and thereby facilitate searching for a potential specic drug target to modulate the central clock.

Methods

Mouse strains. Gpr176 / mice (Acc. No. CDB0672K: http://www.cdb.riken.jp/arg/mutant&percnt;20mice&percnt;20list.html

Web End =http://www.cdb.riken.jp/ http://www.cdb.riken.jp/arg/mutant&percnt;20mice&percnt;20list.html

Web End =arg/mutant%20mice%20list.html ) were generated in the RIKEN CDB (Kobe,

Japan) and backcrossed to C57BL/6J for 10 generations. Then, Gpr176 / mice were intercrossed to produce homozygous null and wild-type progenies for behavioural tests. Vipr2 / mice24 were bred on C57BL/6J background58.

A cohort of Gpr176 / ;Vipr2 / mice and control siblings was produced by crossing double heterozygotes (Gpr176 / ;Vipr2 / ) using in vitro fertilization.

Cry-null mice (Cry1 / ;Cry2 / ) were bred as described previously38,39,61. Targeted mutant mice for Gpr19 (Gpr19tm1Dgen) and Calcr (Calcrtm1Dgen) were obtained from the Mutant Mouse Regional Resource Center at the University of North Carolina with a mixed genetic background involving 129P2/OlaHsd

C57BL/6J (https://www.mmrrc.org/

Web End =https://www.mmrrc.org/). All animal experiments were performed under protocols approved by the Animal Care and Experimentation Committee of Kyoto University.

Behavioural activity monitoring. Single-caged adult male littermate mice (8- to 15-week old) were housed individually in lighttight, ventilated closets within a temperature- and humidity-controlled facility. The animals were entrained on a 12-h light (B200 lux uorescent light):12-h dark (LD) cycle for at least 2 weeks and then transferred to DD. Locomotor activity was detected with passive (pyroelectric) infrared sensors (FA-05 F5B; Omron) and the data were analysed with ClockLab software (Actimetrics) developed on MatLab (Mathworks)29. Free-running circadian period was determined with w2 periodogram, based on animal behaviours in a 21-day interval taken 3 days after the start of DD condition. Fast Fourier transform (FFT) spectral analysis of the activity records was conducted with MATLAB Signal Processing Toolbox 6.2 (Mathworks). To extract long-term locomotor activity trends, we applied a moving average with a 3.67-h window size three times to the original locomotor activity data collected every20 min. Then, FFT spectrograms were created through specgram command, with window size 512, overlap set to 506, and sampling set to 72 cycles per day. For light-pulse-induced shift experiments, mice put in DD were exposed to a 15-min light pulse at either CT14 or CT22. Phase shifts (delay at CT14, advance at CT22) were quantied as the time difference between regression lines of activity onset before and after the light application, using ClockLab software.

Per1-luc organotypic tissue slice culture. Per1-luc transgenic mice carry a rey luciferase reporter gene linked to a 7.2-kb genomic DNA fragment covering the 50-upstream region of the mouse Per1 gene62. Per1-luc-Gpr176 / mice were generated by crossing Per1-luc mice with Gpr176 / mice. The SCN slices were prepared according to our standard method36 and kept at 35 C in a sealed 35 mm Petri dish with 1 ml of the culture medium containing 1 mM D-luciferin. Bioluminescence from the cultured SCN was measured with a highly sensitive cryogenic CCD camera (800S: Spectral Instruments) equipped with a microscope (Axiovert 200: Carl Zeiss). Recording was performed every 20 min. Observed data of images were ltered through a median ltre to eliminate cosmic-ray-induced background noise using ImageJ (http://imagej.nih.gov/ij

Web End =http://imagej.nih.gov/ij). For period determination, the bioluminescence values of whole SCN in the image sequence were exported into Excel (Microsoft), where values were detrended by subtraction of baseline bioluminescence based on a running average from 12 h before to 12 h after each time point. Then, the baseline-subtracted data were curve tted to a modied damped sine wave in Prism (Graphpad software) using the following equation: Y Amplitude exp( K X) sin ((2 p X/Period) Phase 2

p/Period), where K is the damping constant, and the period was determined based on the best-t results. As for lung explant culture, culture was performed according to a published method63 with slight modications. In brief, lungs taken from 5-day-old pups were sliced into a small piece (B2 2 0.3 mm3), placed on a

Millicell membrane (PICMORG50, Millipore) with 800 ml of DMEM medium (Sigma), supplemented with 10 mM HEPES (pH 7.2), 2% B27 (Invitrogen),25 units ml 1 penicillin, 25 mg ml 1 streptomycin, and 1 mM luciferin, in 35-mm dish, and air sealed. Bioluminescence was continuously monitored without interruption for 45 d immediately upon placement in culture with a dish-type photon countable luminometer (Kronos Dio, ATTO) at 35 C. Period of circadian luminescence was determined as described for the SCN.

Radioisotopic in situ hybridization. In situ hybridization was performed with free-oating brain sections (30-mm thick), using [33P]-labelled cRNA probes for Gpr176 (nucleotides 33213746, NM_201367), Ntsr2 (10011480, NM_008747), Gpr37l1 (13681849, NM_134438), Lphn1 (76178072, NM_181039), Oprl1 (15132006, NM_011012), Gpr37 (20562528, NM_010338), Gprc5b (22692837, NM_022420), Gpr85 (19422440, NM_145066), Darc (6951103, NM_010045), Lphn3 (54515771, NM_198702), Lphn2 (52135697, NM_001081298), Calcr (6701168, NM_007588), Gpr48 (43204803, NM_172671), Gpr56 (18442334, NM_018882), Gpr19 (250748, NM_008157), Gpr123 (39554486, NM_177469), Gpr83 (15182057, NM_010287), Ackr3 (9311392, NM_007722), Gpr125 (39194377, NM_133911), Gpr22 (25013002, NM_175191), Gpr153 (32053689, NM_178406), Gpr45 (11651722, NM_053107), Gpr68 (26703170, NM_175493), Gpr6 (14061747, AK139367), Cckar (21822593, NM_009827), and Gpr75 (22692739, NM_175490). All fragments were sequenced to verify their identity, and antisense riboprobes were generated.

Digoxigenin in situ hybridization. Digoxigenin in situ hybridization was performed according to our standard method64 with two different probes for Gnaz (NM_010311), a 50UTR probe (296 bp, nucleotides 151446) and a 30UTR probe (232 bp, nucleotides 1,9192,150), the sequences of which are divergent from those of the other Gi/o family members.

Northern blotting. Northern blot analysis was performed with the following probes for Gpr176 (NM_201367): 50 probe (289 bp, nucleotides 132420) and30 probe (426 bp, nucleotides 3,3213,746). Both fragments were labelled with [32P] deoxycytidine triphosphate by random priming and hybridized with the mouse MTN blots (Clontech) to which poly(A) RNA fractions from various tissues

were transferred (2 mg for each tissue).

Microarray analysis. Microarray data have been deposited in the Gene Expression Omnibus under accession code GSE28574 (ref. 29). To identify genes that are enriched in the SCN, SCN punches taken from 10 animals at CT2 and 10 animals at CT14 were pooled together and analysed with a GeneChip Mouse Genome 430 2.0 (Affymetrix)29. The data were normalized with the MAS5 (GCOS 1.4) algorithm, using the default analysis settings and global scaling for normalization. For statistical analysis of the microarray data, we transformed the values into log2 format for ease of comparison and data representation. We then obtained values of 13 probes for the seminal proteins (Svs3, Svs5, Svs6, Sva, Svp2, Sval1 and Sval2) that were unlikely to be expressed in the SCN. We determined the mean and s.d. of these values. We took this mean value to represent zero expression and subtracted it from the value for each receptor for which we obtained a signal. The value of two s.d. was then considered to be the threshold or baseline for our receptors of interest as any values at this level or higher would be considered statistically signicant (95% condence that the true mean of the seminal protein expression would fall within this threshold value)49.

Laser microdissection of the SCN. Coronal brain section (30-mm thick) containing the SCN was prepared using a cryostat microtome (CM3050S, Leica) and mounted on POL-membrane slides (Leica). Sections were xed for 3 min in an ice-cold mixture of ethanol and acetic acid (19:1), rinsed briey in ice-cold water, stained for 30 seconds in ice-cold water containing 0.05% toluidine blue, followed by two brief washes in ice-cold water. After wiping off excess water, slides were quickly air dried at room temperature. As soon as moistures in the sections decreased enough for laser-cutting, cells in the SCN were microdissected using a LMD7000 device (Leica) and lysed in Trizol reagent (Invitrogen), and total RNA was puried using the RNeasy micro kit (Qiagen).

qRTPCR. qRTPCR analysis was performed as described previously with a standard curve method65,66. The data were normalized with Rplp0. The primer sets used were following: for mouse, Gpr176 (NM_201367), Fw: 50-CATCTTCATTGG

CTCGCTAC-30, Rv: 50-CGTATAGATCCACCAGCAAC-30; Gnaz (NM_010311), Fw: 50-CAGCCGTGCTTAGAAACATCG-30, Rv: 50-TCTAGTGACACTCCACC TCC-30; Gnai1 (NM_010305), Fw: 50-AAGCTGACTCGCCTTCCCAG-30, Rv: 50-GTAGTTTACAGTTCTCCACACG-30; Gnai2 (NM_008138), Fw: 50-TGCCTT GAGTGTGTCTGCGTG-30, Rv: 50-CTCAGTGACGTTGGCAGTTG-30; Gnai3 (NM_010306), Fw: 50-GTGCAGTCCGTGTACAAGAG-30, Rv: 50-GATGAATGG ATCCGAGCCAC-30; Per1 (NM_011065), Fw: 50-TGGCTCAAGTGGCAATGA GTC-30, Rv: 50-GGCTCGAGCTGACTGTTCACT-30; and Rplp0 (NM_007475), Fw: 50-CTCACTGAGATTCGGGATATG-30, Rv: 50-CTCCCACCTTGTCTCCA GTC-30; and for human, GNAZ (NM_002073), Fw: 50- CTACGAGGATAACCA GAC-30, Rv: 50-TACGTGTTCTGGCCCTTG-30; GNAI1 (NM_002069), Fw: 50-CATCTCTGACCTTGTTTCAGC-30, Rv: 50-CTTCAACCCAGTGACAAC ACG-30; GNAI2 (NM_002070), Fw: 50-ACTCCGTGCCTTGAGTGTG-30, Rv: 50-TTGTCTGGAACAGCCCTTGG-30; GNAI3 (NM_010306), Fw: 50-GGAAAG TTACGTTCACTTCAACC-30, Rv: 50-TTGGACCCCAAAAGGCACTG-30; and RPLP0 (NM_053275), Fw: 50-ATGCAGCAGATCCGCATGT-30, Rv: 50-TTGCGC ATCATGGTGTTCTT-30.

10 NATURE COMMUNICATIONS | 7:10583 | DOI: 10.1038/ncomms10583 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10583 ARTICLE

Antibodies to Gpr176 and Vipr2. Gpr176 antibody was raised in rabbit using a glutathione-S-transferase-fused Gpr176 mouse protein fragment (amino acids (a.a.) 311515). The raised antibodies were afnity-puried using a maltose-binding protein (MBP)-fused Gpr176 fragment (a.a. 311515). Vipr2 antibody was raised in chicken with a keyhole-limpet hemocyanin (KLH)-conjugated synthetic peptide mapping a C-terminal region of the mouse Vipr2 (a.a. 418 437). The antibodies

were afnity-puried using the antigen peptide. Rabbit antiserum against Vipr2 was purchased from Abcam (ab28624).

Immunohistochemistry. Free-oating immunohistochemistry was performed with 30-mm-thick serial coronal brain sections. To minimize technical variations in immunostaining, sections from different CTs were immunolabelled simultaneously. The primary antibodies used were anti-Gpr176 (nal concentration,0.6 mg ml 1), anti-Vipr2 rabbit antiserum (Abcam, 1:1,000), and anti-Vipr2 puried chicken polyclonal antibody (0.17 mg ml 1). Immunoreactivities were visualized with a peroxidase-based Vectorstain Elite ABC kit (Vector Laboratories)

using diaminobenzidine chromogen. For brightness and contrast, photo-micrographs were processed identically with ImageJ. For quantitative analysis, data were normalized with respect to the difference between signal intensities in equal areas of the SCN and the corpus callosum. Normalized values were summed from the rostral to the caudal margins of the SCN (10 sections per brain), and the sum was considered a measure for the amount of protein in the SCN. Valuesare expressed as meanss.e.m. (n 6, for each time point). For dual-label

immunouorescence, free-oating sections were stained with anti-Gpr176 (rabbit polyclonal, nal concentration, 0.6 mg ml 1) together with either anti-

Vipr2 (chicken polyclonal, 0.17 mg ml 1), anti-Vip (Abnova, guinea pig polyclonal, PAB16648, 1:1,000), or anti-Avp-associated neurophysin II (Santa Cruz, goat polyclonal, sc-27093, 0.2 mg ml 1) antibody. We visualized immunoreactivities using Alexa594-conjugated anti-rabbit IgG (1:1,000; Life Technologies) and

Alexa488-conjugated anti-chicken, guinea pig, or goat IgG (1:1,000; Life Technologies). Nuclei were visualized by staining with 40,60-diamino-2-phenylindole (DAPI).

Immunocytochemistry of dispersed SCN neurons. For colocalization analysis of Gpr176 and Vipr2, dissociated SCN neuronal cultures were used. SCN punches of 15 pups (C57Bl/6) of 45 days of age were pooled and incubated for 40 min at 37 C in Ca2 /Mg2 -free Hanks balanced salt solution (HBSS, Life Technologies)

containing 0.06% papain, 0.02% L-cysteine, and 1 kU ml 1 DNaseI (Sigma). Following the addition of fetal bovine serum to the solution, cells were further dissociated through trituration with a re-polished Pasteur pipette. The dispersals were then ltered through a 100-mm nylon cell strainer (BD Falcon) and resuspended in Neurobasal medium (Life Technologies) containing B27 supplement (Life Technologies) with 2 mM glutamine, 8 mM glutamate, 100 U ml 1 penicillin, and 100 mg ml 1 streptomycin. Cell viability was 8595%. Viable cells were plated on polylysine-coated coverslips in 48-well culture plates at a density of 2,000 cells per mm2. One-half of the culture medium was exchanged with fresh medium every 3 days. After 8 days culture, cells were xed with 4% paraformaldehyde and double-labelled for Gpr176 and Vipr2 using published methods67. We stained cultures with anti-Gpr176 (puried rabbit polyclonal, nal concentration 0.15 mg ml 1) and anti-Vipr2 (puried chicken polyclonal,0.15 mg ml 1) followed by Alexa594-conjugated antirabbit IgG (1:1,000; Life Technologies) and Alexa488-conjugated antichicken IgG (1:1,000; Life

Technologies). Nuclei were visualized using DAPI staining.

Stable cell lines. Flp-In TREx293-Gpr176 cells were generated by stable transfection of Flp-In T-Rex-293 cells (Life Technologies) with a pcDNA5/FRT vector (Life Technologies) containing the untagged full-length coding sequence of the mouse Gpr176 (NM_201367). Similarly, Flp-In TREx293-Vipr2 cells were established with the mouse Vipr2 full-length coding sequence (NM_009511).

To develop Flp-In TREx293-Gpr176(tet-on)/Vipr2 cells, we constructed a modied pcDNA5/FRT vector carrying Vipr2 and Gpr176 under different promoters: while Gpr176 was cloned into a proprietary pcDNA5/FRT cloning site for tet-on induction, Vipr2 was cloned separately into a different position of the vector(at a unique PciI site) in conjunction with a tetracycline-insensitive CMV promoter. Point mutations for Gpr176RDY and Gpr176V145R were introduced into the corresponding constructs with a standard sequential PCR method68. For stable expression of DN-Gz mutant (G204A/E246A/A327S)52, we established Flp-In TREx293-Gpr176(tet-on)/DN-Gz cells by using the modied pcDNA5/FRT vector as mentioned above. For infection of RGSZ1, lentiviruses carrying the hemagglutinin (HA)-tagged full-length coding sequence of RGSZ1 (NM_003702) in tandem with IRES-GFP (HA-RGSZ1-IRES-GFP) were generated with pCSII-EFMCS-IRES-hrGFP vector69, and the cells were infected with an empirical titre of virus that resulted in nearly 100% infection as determined by GFP expression. Stable clonal NIH 3T3 cell lines expressing Gpr176 (NIH3T3 Tet-on 3G-Gpr176) were generated via transfection of NIH3T3 Tet-On 3G cells (Clontech) with a pTRE3G vector (Clontech) containing Gpr176. The established cells were further transfected with a pEF vector (Addgene) containing Gz (NM_010311) along with Linear Puromycin Marker (Clontech) to generate NIH3T3 Tet-on 3G-Gpr176 (tet-on)/Gz double stable cell clones. Cells were cultured in DMEM containing 10%

fetal bovine serum with an appropriate mixture of antibiotics that the manufacturers recommend for the maintenance of the cell clones.

Immunoblot. To avoid high-temperature-induced protein aggregation of GPCR, cell lysates were denatured on ice in Laemmli buffer and subjected to SDSPAGE at 4 C. Immunoblotting was performed using our standard method29 with afnity-puried antibodies against Gpr176 (rabbit polyclonal, nal concentration,0.6 mg ml 1) and Vipr2 (chicken polyclonal, 0.4 mg ml 1). To detect Gz,the plasma membrane fractions were lysed in standard RIPA buffer and immunoprecipitated with anti-Gz antibody (Santa Cruz, sc-388, 1 mg per immunoprecipitaion), and immunoblots were probed with the rabbit antisera against Gz (2919) provided by Dr Manning48,70 (1:100 dilution). Commercially available antibodies against Gi (Abcam, ab3522, 1 mg ml 1) and a-tubulin (Sigma,

T6199, 1 mg ml 1) were used as a control. Blot images have been cropped for presentation. Full size images are presented in Supplementary Fig. 12.

Measurements of cAMP and IP1. After 24 h of treatment with Dox (1 mg ml 1) or vehicle, cells were removed from culture dish with Versene solution (Life

Technologies) and dissociated into single cells through gentle trituration. After ltration with a 100-mm cell strainer (BD Falcon), cells were resuspended in HBSS containing 5 mM HEPES (pH 7.5), 0.1% bovine serum albumin, and 0.5 mM 3-isobutyl-1-methylxanthine (IBMX, a non-selective phosphodiesterase inhibitor). Cells in suspension were incubated at 37 C for 1 h in 24-well plates at a density of 3 105 cells per well, followed by stimulation with forskolin (Nacalai Tesque)

or Vip (Peptide Institute Inc.) at the indicated concentrations for 15 min. The reactions were stopped by adding ice-cold perchloric acid (1 N, nal solution) containing 4 mM theophylline (Sigma). After 1 h at 4 C, the mixtures were centrifuged, and the supernatants were neutralized with ice-cold 0.72 M KOH/0.6 M KHCO3. Following removal of salt precipitants, the extracts were assayed for cAMP concentrations with a cAMP-specic enzyme immunoassay kit (Cayman Chemical)29. Assays with PTX-treated cells were also done with the same protocols, except that the cells were cultured in the presence of PTX (100 ng ml 1, Bio

Academia) for 16 h before assay. When specied, 100 mM sphingosine-1-phosphate (Enzo Life Sciences) was added together with forskolin, to conrm whether Gi-mediated signalling was blocked efciently by the PTX treatment. For IP1 assay, cells in suspension were incubated in IP stimulation buffer (Cisbio; 10 mM HEPES (pH 7.4), 1 mM CaCl2, 0.5 mM MgCl2, 4.2 mM KCl, 146 mM NaCl, 5.5 mM glucose and 50 mM LiCl) for 1 h at 37 C in 24-well plate at a density of 4 105 cells

per well. The reaction was stopped by adding Lysis reagent (Cisbio). The concentrations of IP1 were determined by enzyme immunoassay with IP-One

ELISA kit (Cisbio) according to the manufacturers protocol.

[125I]-Vip binding assay. To prepare membrane fractions, Flp-In TREx293-Gpr176(tet-on)/Vipr2 cells were lysed in hypotonic lysis buffer (25 mM HEPES, pH 7.5, 1 mM EDTA, 1 mM DTT) with protease inhibitors, and passed through a 27-gauge needle 10 times. The lysed cells were centrifuged at 700 g for 5 min to remove nuclei and debris. The supernatant was centrifuged at 20,400g for 30 min. Then, the pellet (membrane fraction) was resuspended in HBSS containing 5 mM HEPES (pH 7.5). For the binding reaction, the membranes (10 mg of protein) were incubated with 1, 10 or 100 nM of [125I]-labelled Vip (PerkinElmer; NEX192) in HBSS binding buffer containing 5 mM HEPES (pH 7.5) and 0.1% bovine serum albumin for 60 min at 4 C. Incubation was performed in a 1.5-ml siliconized tube (Sarstedt). At the termination of incubation, membrane-bound [125I]-Vip was separated from free peptide by centrifugation. The membrane pellets were washed three times with HBSS/5 mM HEPES/0.1% bovine serum albumin, and the washed membranes were assayed for [125I] radioactivity with a 1470 automatic gamma counter (PerkinElmer). Specic binding was calculated by subtracting the radioactivity detected for un-transfected (parental) Flp-In TREx293 cells.

siRNAs. To knockdown Gz, two independent pools of Gz-specic silencer select siRNAs (#1 and #2), each containing three different siRNA duplexes directed against the Gz coding sequence (#1: s5898, s5900 and s499632; #2: s499631, s500931 and s500937; Life Technologies), were introduced into B30% conuent

Flp-In TREx293-Gpr176 cells using Lipofectamine 2000 (Life Technologies) according to the manufacturers instructions. As a control, we also transfected the cells with negative-control siRNA (catalog number 4390846, Life Technologies) at the same concentration as the #1 and #2 mixtures (1.8 nmol siRNA per 10-cm dish). Medium was replaced 6 h after transfection. Four days later, cells were treated with Dox (1 mg ml 1) or vehicle for 24 h and removed from dishes with

Versene solution (Life Technologies) to be subject to cAMP assays. For each experiment, we took a fraction of cells and conrmed by qRTPCR that Gz mRNA accumulation was diminished to o5% by RNA interference.

SCN punch for cAMP measurement. The microdissection of the SCN was performed as described with modications. Animals kept in DD were killed by cervical dislocation, and the eyes were removed under a safety red light. The brain was then isolated from the skull under room light and frozen immediately on dry ice. Coronal brain section (300-mm thick) containing the SCN was prepared using a

NATURE COMMUNICATIONS | 7:10583 | DOI: 10.1038/ncomms10583 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 11

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10583

cryostat microtome (CM3050S, Leica) and mounted on a silicon rubber stage at

17 C. Under a magnifying glass, the bilateral SCN was punched out fromthe frozen section using a blunt 20-gauge syringe needle whose edge had been sharpened by ling. The microdissected SCN (one punch per assay) was then immediately sonicated at 4 C (Bioruptor, COSMO BIO) in 0.1 N HCl solution containing 0.5 mM IBMX. Lysates were claried by centrifugation, and the protein content was determined with a Bradford assay (Nakarai). The amount of extracted cAMP was measured using an enzyme immunoassay kit for cAMP (Cayman Chemical). To increase the sensitivity of the assay, samples were acetylated according to the manufacturers instruction (Cayman Chemical).

GloSensor cAMP assay. We generated additional Gpr176 mutants on the DRY triplet sequence (Supplementary Fig. 8; DAY, AAY, and AAA) and tested their basal activities using the cAMP GloSensor system (Promega), which allows transient transfection-based GPCR assay for cAMP signal. Flp-In TREx293 cells were plated on 35-mm dishes at 1.2 106 cells per dish with a CO2-independent DMEM

medium (Sigma) containing 10 mM HEPES (pH 7.2), 10% FBS, 100 units ml 1 penicillin, 100 mg ml 1 streptomycin at 37 C. In the next day, the cells were transfected with a DNA mixture containing 1 mg of pGloSensor-22F plasmid (Promega) and 1.5 mg of pcDNA3.1 expression plasmid (Life Technologies) encoding either the untagged wild-type Gpr176 or its respective DRY mutants, using Lipofectamine LTX/Plus reagent (Life Technologies). Four hours after the transfection, the medium was refreshed to the medium with 1 mM luciferin, and the cells were further cultured at 37 C for 14 h. Following 2-h incubation at 28 C for equilibrium, GloSensor activities in the cells were measured using a dish-type luminometer (Kronos Dio, ATTO) at 28 C. Recording was performed every 1 min with 2 s of integration. After detection of baseline luminescence activities for10 min, Fsk was added to the culture medium at the nal concentration of 10 mM.

Recording was stopped at 30 min after the Fsk stimulation. Then the cells were immediately lysed into Laemmli buffer for western blot analysis.

References

1. Welsh, D. K., Takahashi, J. S. & Kay, S. A. Suprachiasmatic nucleus: cell autonomy and network properties. Annu. Rev. Physiol. 72, 551577 (2010).

2. Mohawk, J. A., Green, C. B. & Takahashi, J. S. Central and peripheral circadian clocks in mammals. Annu. Rev. Neurosci. 35, 445462 (2012).

3. Schibler, U. & Sassone-Corsi, P. A web of circadian pacemakers. Cell 111, 919922 (2002).

4. Liu, A. C., Lewis, W. G. & Kay, S. A. Mammalian circadian signaling networks and therapeutic targets. Nat. Chem. Biol. 3, 630639 (2007).

5. Herzog, E. D. Neurons and networks in daily rhythms. Nat. Rev. Neurosci. 8, 790802 (2007).

6. Bass, J. Circadian topology of metabolism. Nature 491, 348356 (2012).7. Paschos, G. K., Baggs, J. E., Hogenesch, J. B. & FitzGerald, G. A. The role of clock genes in pharmacology. Annu. Rev. Pharmacol. Toxicol. 50, 187214 (2010).

8. Hirota, T. et al. Identication of small molecule activators of cryptochrome. Science 337, 10941097 (2012).

9. Solt, L. A. et al. Regulation of circadian behaviour and metabolism by synthetic REV-ERB agonists. Nature 485, 6268 (2012).

10. Rask-Andersen, M., Almen, M. S. & Schioth, H. B. Trends in the exploitation of novel drug targets. Nat. Rev. Drug. Discov. 10, 579590 (2011).

11. Civelli, O. et al. G protein-coupled receptor deorphanizations. Annu. Rev. Pharmacol. Toxicol. 53, 127146 (2013).

12. Reinscheid, R. K. et al. Orphanin FQ: a neuropeptide that activates an opioidlike G protein-coupled receptor. Science 270, 792794 (1995).

13. Sakurai, T. et al. Orexins and orexin receptors: a family of hypothalamic neuropeptides and G protein-coupled receptors that regulate feeding behavior. Cell 92, 573585 (1998).

14. Hinuma, S. et al. A prolactin-releasing peptide in the brain. Nature 393,

272276 (1998).

15. Kojima, M. et al. Ghrelin is a growth-hormone-releasing acylated peptide from stomach. Nature 402, 656660 (1999).

16. Seifert, R. & Wenzel-Seifert, K. Constitutive activity of G-protein-coupled receptors: cause of disease and common property of wild-type receptors. Naunyn Schmiedebergs Arch. Pharmacol. 366, 381416 (2002).

17. Bond, R. A. & Ijzerman, A. P. Recent developments in constitutive receptor activity and inverse agonism, and their potential for GPCR drug discovery. Trends Pharmacol. Sci. 27, 9296 (2006).

18. Milligan, G. Constitutive activity and inverse agonists of G protein-coupled receptors: a current perspective. Mol. Pharmacol. 64, 12711276 (2003).

19. Mukhopadhyay, S. et al. The ciliary G-protein-coupled receptor Gpr161 negatively regulates the Sonic hedgehog pathway via cAMP signaling. Cell 152, 210223 (2013).

20. Mehlmann, L. M. et al. The Gs-linked receptor GPR3 maintains meiotic arrest in mammalian oocytes. Science 306, 19471950 (2004).

21. ONeill, J. S., Maywood, E. S., Chesham, J. E., Takahashi, J. S. & Hastings, M. H. cAMP-dependent signaling as a core component of the mammalian circadian pacemaker. Science 320, 949953 (2008).

22. An, S., Irwin, R. P., Allen, C. N., Tsai, C. & Herzog, E. D. Vasoactive intestinal polypeptide requires parallel changes in adenylate cyclase and phospholipase C to entrain circadian rhythms to a predictable phase. J. Neurophysiol. 105, 22892296 (2011).

23. Rea, M. A. VIP-stimulated cyclic AMP accumulation in the suprachiasmatic hypothalamus. Brain Res. Bull. 25, 843847 (1990).

24. Harmar, A. J. et al. The VPAC(2) receptor is essential for circadian function in the mouse suprachiasmatic nuclei. Cell 109, 497508 (2002).

25. Aton, S. J., Colwell, C. S., Harmar, A. J., Waschek, J. & Herzog, E. D. Vasoactive intestinal polypeptide mediates circadian rhythmicity and synchrony in mammalian clock neurons. Nat. Neurosci. 8, 476483 (2005).

26. Wong, Y. H., Conklin, B. R. & Bourne, H. R. Gz-mediated hormonal inhibition of cyclic AMP accumulation. Science 255, 339342 (1992).

27. Fong, H. K., Yoshimoto, K. K., Eversole-Cire, P. & Simon, M. I. Identication of a GTP-binding protein alpha subunit that lacks an apparent ADP-ribosylation site for pertussis toxin. Proc. Natl Acad. Sci. USA 85, 30663070 (1988).

28. Matsuoka, M., Itoh, H., Kozasa, T. & Kaziro, Y. Sequence analysis of cDNA and genomic DNA for a putative pertussis toxin-insensitive guanine nucleotide-binding regulatory protein alpha subunit. Proc. Natl Acad. Sci. USA 85, 53845388 (1988).

29. Doi, M. et al. Circadian regulation of intracellular G-protein signalling mediates intercellular synchrony and rhythmicity in the suprachiasmatic nucleus. Nat. Commun. 2, 327 (2011).

30. Hannibal, J., Brabet, P. & Fahrenkrug, J. Mice lacking the PACAP type I receptor have impaired photic entrainment and negative masking. Am. J. Physiol. Regul. Integr. Comp. Physiol. 295, R2050R2058 (2008).

31. Hannibal, J. et al. Dissociation between light-induced phase shift of the circadian rhythm and clock gene expression in mice lacking the pituitary adenylate cyclase activating polypeptide type 1 receptor. J. Neurosci. 21, 48834890 (2001).

32. Cheng, M. Y. et al. Prokineticin 2 transmits the behavioural circadian rhythm of the suprachiasmatic nucleus. Nature 417, 405410 (2002).

33. Prosser, H. M. et al. Prokineticin receptor 2 (Prokr2) is essential for the regulation of circadian behavior by the suprachiasmatic nuclei. Proc. Natl Acad. Sci. USA 104, 648653 (2007).

34. Yamaguchi, Y. et al. Mice genetically decient in vasopressin V1a and V1b receptors are resistant to jet lag. Science 342, 8590 (2013).

35. Li, J. D., Burton, K. J., Zhang, C., Hu, S. B. & Zhou, Q. Y. Vasopressin receptor V1a regulates circadian rhythms of locomotor activity and expression of clock-controlled genes in the suprachiasmatic nuclei. Am. J. Physiol. Regul. Integr. Comp. Physiol. 296, R824R830 (2009).

36. Yamaguchi, S. et al. Synchronization of cellular clocks in the suprachiasmatic nucleus. Science 302, 14081412 (2003).

37. Hata, S., Emi, Y., Iyanagi, T. & Osumi, T. cDNA cloning of a putative G protein-coupled receptor from brain. Biochim. Biophys. Acta. 1261, 121125 (1995).

38. van der Horst, G. T. et al. Mammalian Cry1 and Cry2 are essential for maintenance of circadian rhythms. Nature 398, 627630 (1999).

39. Vitaterna, M. H. et al. Differential regulation of mammalian period genes and circadian rhythmicity by cryptochromes 1 and 2. Proc. Natl Acad. Sci. USA 96, 1211412119 (1999).

40. Aton, S. J. & Herzog, E. D. Come together, right...now: synchronization of rhythms in a mammalian circadian clock. Neuron 48, 531534 (2005).

41. Moore, R. Y. Circadian rhythms: basic neurobiology and clinical applications. Annu. Rev. Med. 48, 253266 (1997).

42. An, S., Tsai, C., Ronecker, J., Bayly, A. & Herzog, E. D. Spatiotemporal distribution of vasoactive intestinal polypeptide receptor 2 in mouse suprachiasmatic nucleus. J. Comp. Neurol. 520, 27302741 (2012).

43. Hughes, A. T. & Piggins, H. D. Behavioral responses of Vipr2 / mice to

light. J. Biol. Rhythms 23, 211219 (2008).44. Damaj, B. B. et al. Identication of G-protein binding sites of the human interleukin-8 receptors by functional mapping of the intracellular loops. FASEBJ. 10, 14261434 (1996).45. Greasley, P. J. et al. Mutational and computational analysis of the alpha(1b)-adrenergic receptor. Involvement of basic and hydrophobic residues in receptor activation and G protein coupling. J. Biol. Chem. 276, 4648546494 (2001).

46. Imai, T., Suzuki, M. & Sakano, H. Odorant receptor-derived cAMP signals direct axonal targeting. Science 314, 657661 (2006).

47. van Koppen, C. et al. Activation of a high afnity Gi protein-coupled plasma membrane receptor by sphingosine-1-phosphate. J. Biol. Chem. 271, 20822087 (1996).

48. Yang, J. et al. Loss of signaling through the G protein, Gz, results in abnormal platelet activation and altered responses to psychoactive drugs. Proc. Natl Acad. Sci. USA 97, 99849989 (2000).

49. Atwood, B. K., Lopez, J., Wager-Miller, J., Mackie, K. & Straiker, A. Expression of G protein-coupled receptors and related proteins in HEK293, AtT20, BV2, and N18 cell lines as revealed by microarray analysis. BMC Genomics 12, 14 (2011).

12 NATURE COMMUNICATIONS | 7:10583 | DOI: 10.1038/ncomms10583 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10583 ARTICLE

69. Kuwata, H. et al. IL-10-inducible Bcl-3 negatively regulates LPS-induced TNF-alpha production in macrophages. Blood 102, 41234129 (2003).

70. Carlson, K. E., Brass, L. F. & Manning, D. R. Thrombin and phorbol esters cause the selective phosphorylation of a G protein other than Gi in human platelets. J. Biol. Chem. 264, 1329813305 (1989).

Acknowledgements

We thank members of the H.O. laboratory, Drs Akira Hirasawa, and Masahiro Matsuo for technical assistance and discussion. We also thank Dr William J. Schwartz (University of Massachusetts Medical School) for stimulating discussions and kind reading of this manuscript. We are grateful to Dr David R. Manning (University of Pennsylvania) for providing an antiserum against Gz (2919) and to Dr Mitsuo Tagaya (Tokyo University of Pharmacy and Life Science) for providing a DN-Gz expression plasmid. This research was supported by Core Research for Evolutional Science and Technology (CREST), Japan Science and Technology Agency (to H.O.), and the Ministry of Education, Culture, Sports, Science and Technology of Japan through a Funding Program for Next Generation World-Leading Researchers (to M.D.), and grants from Takeda Science Foundation and Kobayashi International Scholarship Foundation (to H.O.), Senri Life Science Foundation, Daiichi Sankyo Foundation of Life Science, Inoue Foundation for Science, Brain Science Foundation, Mochida Memorial Foundation for Medical and Pharmaceutical Research, and Kanae Foundation for the Promotion of Medical Science (to M.D.). M.H.H. is supported by Medical Research Council, UK.

50. Glick, J. L., Meigs, T. E., Miron, A. & Casey, P. J. RGSZ1, a Gz-selective regulator of G protein signaling whose action is sensitive to the phosphorylation state of Gzalpha. J. Biol. Chem. 273, 2600826013 (1998).

51. Wang, J. et al. RGSZ1, a Gz-selective RGS protein in brain. Structure, membrane association, regulation by Galphaz phosphorylation, and relationship to a Gz gtpase-activating protein subfamily. J. Biol. Chem. 273, 2601426025 (1998).

52. Nagahama, M. et al. Inactivation of Galpha(z) causes disassembly of the Golgi apparatus. J. Cell Sci. 115, 44834493 (2002).

53. Okamura, H. Suprachiasmatic nucleus clock time in the mammalian circadian system. Cold Spring Harb. Symp. Quant. Biol. 72, 551556 (2007).

54. Yamazaki, S. et al. Resetting central and peripheral circadian oscillators in transgenic rats. Science 288, 682685 (2000).

55. Thathiah, A. et al. The orphan G protein-coupled receptor 3 modulates amyloid-beta peptide generation in neurons. Science 323, 946951 (2009).

56. Nakashima, A. et al. Agonist-independent GPCR activity regulates anterior-posterior targeting of olfactory sensory neurons. Cell 154, 13141325 (2013).

57. Travnickova-Bendova, Z., Cermakian, N., Reppert, S. M. & Sassone-Corsi, P. Bimodal regulation of mPeriod promoters by CREB-dependent signaling and CLOCK/BMAL1 activity. Proc. Natl Acad. Sci. USA 99, 77287733 (2002).

58. Maywood, E. S., Chesham, J. E., OBrien, J. A. & Hastings, M. H. A diversity of paracrine signals sustains molecular circadian cycling in suprachiasmatic nucleus circuits. Proc. Natl Acad. Sci. USA 108, 1430614311 (2011).

59. Riobo, N. A., Saucy, B., Dilizio, C. & Manning, D. R. Activation of heterotrimeric G proteins by Smoothened. Proc. Natl Acad. Sci. USA 103, 1260712612 (2006).

60. Rasmussen, S. G. et al. Crystal structure of the beta2 adrenergic receptor-Gs protein complex. Nature 477, 549555 (2011).

61. Doi, M. et al. Salt-sensitive hypertension in circadian clock-decient Cry-null mice involves dysregulated adrenal Hsd3b6. Nat. Med. 16, 6774 (2010).

62. Yamaguchi, S. et al. The 5 upstream region of mPer1 gene contains two promoters and is responsible for circadian oscillation. Curr. Biol. 10, 873876 (2000).

63. Yamazaki, S. & Takahashi, J. S. Real-time luminescence reporting of circadian gene expression in mammals. Methods Enzymol. 393, 288301 (2005).

64. Shigeyoshi, Y. et al. Light-induced resetting of a mammalian circadian clock is associated with rapid induction of the mPer1 transcript. Cell 91, 10431053 (1997).

65. Hou, Y., Zhang, H., Miranda, L. & Lin, S. Serious overestimation in quantitative PCR by circular (supercoiled) plasmid standard: microalgal pcna as the model gene. PLoS ONE 5, e9545 (2010).

66. Yamamura, K. et al. Immunolocalization of murine type VI 3beta-hydroxysteroid dehydrogenase in the adrenal gland, testis, skin, and placenta. Mol. Cell Endocrinol. 382, 131138 (2014).

67. Doi, M. et al. Isoform-specic monoclonal antibodies against 3beta-hydroxysteroid dehydrogenase/isomerase family provide markers for subclassication of human primary aldosteronism. J. Clin. Endocrinol. Metab. 99, E257E262 (2014).

68. Cormack, B. Directed mutagenesis using the polymerase chain reaction. Curr. Protoc. Mol. Biol. Chapter 8, Unit 8.5 (2001).

Author contributions

M.D. and H.O. conceived the project and designed the research; M.D., I.M. and S.K. are equally contributing rst authors, who performed experiments in collaboration with G.S., N.U., R.T., S.K., H.S., H.H., H.W.C., Y.N., Y.T., Y.H., J.Y., M.M., M.H.H, H.K. and H.O.; M.D. and H.O. drafted the manuscript, supported by I.M., S.K. and M.H.H.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Doi, M. et al. Gpr176 is a Gz-linked orphan G-protein-coupled receptor that sets the pace of circadian behaviour. Nat. Commun. 7:10583doi: 10.1038/ncomms10583 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

NATURE COMMUNICATIONS | 7:10583 | DOI: 10.1038/ncomms10583 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 13

Word count: 13535

Show less

Abstract

Translate

G-protein-coupled receptors (GPCRs) participate in a broad range of physiological functions. A priority for fundamental and clinical research, therefore, is to decipher the function of over 140 remaining orphan GPCRs. The suprachiasmatic nucleus (SCN), the brain's circadian pacemaker, governs daily rhythms in behaviour and physiology. Here we launch the SCN orphan GPCR project to (i) search for murine orphan GPCRs with enriched expression in the SCN, (ii) generate mutant animals deficient in candidate GPCRs, and (iii) analyse the impact on circadian rhythms. We thereby identify Gpr176 as an SCN-enriched orphan GPCR that sets the pace of circadian behaviour. Gpr176 is expressed in a circadian manner by SCN neurons, and molecular characterization reveals that it represses cAMP signalling in an agonist-independent manner. Gpr176 acts independently of, and in parallel to, the Vipr2 GPCR, not through the canonical Gi, but via the unique G-protein subclass Gz.

Details

Title

Gpr176 is a Gz-linked orphan G-protein-coupled receptor that sets the pace of circadian behaviour

Author

Doi, Masao; Murai, Iori; Kunisue, Sumihiro; Setsu, Genzui; Uchio, Naohiro; Tanaka, Rina; Kobayashi, Sakurako; Shimatani, Hiroyuki; Hayashi, Hida; Chao, Hsu-wen; Nakagawa, Yuuki; Takahashi, Yukari; Hotta, Yunhong; Yasunaga, Jun-ichirou; Matsuoka, Masao; Hastings, Michael H; Kiyonari, Hiroshi; Okamura, Hitoshi

Pages

10583

Publication year

2016

Publication date

Feb 2016

Publisher

Nature Publishing Group

e-ISSN

20411723

Source type

Scholarly Journal

Language of publication

English

DOI

https://doi.org/10.1038/ncomms10583

ProQuest document ID

1765575111

Gpr176 is a Gz-linked orphan G-protein-coupled receptor that sets the pace of circadian behaviour

Jump to:

Full text

Abstract

Details

Suggested sources