Abstract
The major regulator of the neuroendocrine stress response in the brain is corticotropin releasing factor (CRF), whose transcription is controlled by CREB and its cofactors CRTC2/3 (TORC2/3). Phosphorylated CRTCs are sequestered in the cytoplasm, but rapidly dephosphorylated and translocated into the nucleus following a stressful stimulus. As the stress response is attenuated by oxytocin (OT), we tested whether OT interferes with CRTC translocation and, thereby, Crf expression. OT (1 nmol, i.c.v.) delayed the stress-induced increase of nuclear CRTC3 and Crf hnRNA levels in the paraventricular nucleus of male rats and mice, but did not affect either parameter in the absence of the stressor. The increase in Crf hnRNA levels at later time points was parallel to elevated nuclear CRTC2/3 levels. A direct effect of Thr4 Gly7-OT (TGOT) on CRTC3 translocation and Crf expression was found in rat primary hypothalamic neurons, amygdaloid (Ar-5), hypothalamic (H32), and human neuroblastoma (Be(2)M17) cell lines. CRTC3, but not CRCT2, knockdown using siRNA in Be(2)M17 cells prevented the effect of TGOT on Crf hnRNA levels. Chromatin-immunoprecipitation demonstrated that TGOT reduced CRTC3, but not CRTC2, binding to the Crf promoter after 10 min of forskolin stimulation. Together, the results indicate that OT modulates CRTC3 translocation, the binding of CRTC3 to the Crf promoter and, ultimately, transcription of the Crf gene.
SIGNIFICANCE STATEMENT The neuropeptide oxytocin has been proposed to reduce hypothalamic-pituitary-adrenal (HPA) axis activation during stress. The underlying mechanisms are, however, elusive. In this study we show that activation of the oxytocin receptor in the paraventricular nucleus delays transcription of the gene encoding corticotropin releasing factor (Crf), the main regulator of the stress response. It does so by sequestering the coactivator of the transcription factor CREB, CRTC3, in the cytosol, resulting in reduced binding of CRTC3 to the Crf gene promoter and subsequent Crf gene expression. This novel oxytocin receptor-mediated intracellular mechanism might provide a basis for the treatment of exaggerated stress responses in the future.
Introduction
The neuropeptide corticotropin releasing factor (CRF), synthesized in the hypothalamic paraventricular nucleus (PVN), is the main activator of the hypothalamo-pituitary-adrenal (HPA) axis during stress. Following its release into portal vessels of the median eminence, CRF stimulates the release of ACTH from the pituitary gland. CRF-producing neurons are under negative feedback by glucocorticoids, thus preventing the detrimental effects of long-term excessive HPA axis activity. CRF cells are further regulated by other factors during stress, including the neuropeptide oxytocin (OT; Neumann et al., 2000b; Windle et al., 2004). Known for its role in reproduction in the periphery, central OT via its OT receptor (OTR) promote prosocial behavior and dampen HPA axis activity (Donaldson and Young, 2008; Neumann and Landgraf, 2012). Elevated levels of endogenous OT, as observed during the peripartum period (Van Tol et al., 1988; Neumann et al., 1993; Hillerer et al., 2011), as well as chronic infusion of synthetic OT (Windle et al., 2004), reduce stress-induced ACTH response in rats. However, earlier studies have reported potentiating effects of peripheral OT on stress-induced ACTH release (Gibbs et al., 1984; Petersson et al., 1999; Ondrejcakova et al., 2010). Although there is evidence that OT is released within the PVN in response to acute stress (Engelmann et al., 1999), it is unknown whether such locally released OT directly alters the activity of CRF cells during stress.
One of the early hallmarks of CRF cell activity during stress is increased Crf gene transcription via a cAMP/PKA-dependent mechanism; presumably to replenish protein stores within the cells (Aguilera and Liu, 2012). The downstream binding of phosphorylated cAMP-responsive element binding protein (pCREB) to a cAMP-responsive element (CRE) in the promoter is accompanied by binding of a CREB coactivator, CREB-regulated transcriptional coactivators (CRTC1-3) to CREB. Under basal conditions, CRTC is phosphorylated and bound to the scaffolding protein 14-3-3 in the cytoplasm. Following its dephosphorylation, CRTC translocates to the nucleus, where it binds CREB via its bZIP domain acting as a coactivator for the recruitment of CBP/p300 to the Crf promoter for gene transcription to commence (Conkright et al., 2003; Liu et al., 2011). As OT was shown to phosphorylate CREB in the hippocampus (Tomizawa et al., 2003), we hypothesized that OT regulates Crf expression by controlling the translocation and, therefore, nuclear availability of CRTC2 and CRTC3, which play an important role in the regulation of Crf gene transcription (Liu et al., 2011), at least in vitro. In the present study we found that OT delays stress-induced Crf gene transcription in the PVN through inhibition of CRTC3 translocation, and subsequent CRTC3, but not CRTC2, binding to the Crf promoter.
Materials and Methods
Animals
Adult male and female Wistar rats, adult male and female C57BL/6J mice (Charles River Laboratories, 250–300 g and 25–30 g, respectively, at the beginning of the experiment), and transgenic male and female C57BL/6J OTR-Venus reporter mice (Tohoku University Japan; Yoshida et al., 2009) were housed in separate rooms under standard laboratory conditions in groups of 3–4 (12 h light/dark cycle, 22–24°C, lights on at 06:00 h, food and water ad libitum). All animal experiments were performed between 08:00 and 11:00 A.M., in accordance with the Guide for the Care and Use of Laboratory Animals by the National Institutes of Health, and were approved by the government of the Oberpfalz, Germany, and the Institutional Animal Care and Use Committee of Tohoku University.
Stereotaxic implantation of an intracerebroventricular guide cannula in rats and mice and effects of intracerebroventricular infusion of OT or TGOT on protein phosphorylation and anxiety
Stereotaxic implantation of a guide cannula above the third ventricle (rats) or the lateral ventricle (mice) for subsequent intracerebroventricular infusion into the vicinity of the PVN was performed under isoflurane anesthesia and semisterile conditions. Following surgery, the animals received a subcutaneous injection of antibiotics (0.03 ml enrofloxacin; 100 mg/1 ml Baytril, Bayer) and allowed to recover for at least 6 d. Animals were single-housed after surgery and rats were handled daily to habituate them to the central infusion procedure to avoid nonspecific stress responses during the experiment. Mice were left undisturbed in their home cage.
To determine the effect of OT on hypothalamic Crf expression and CRTC translocation, vehicle, OT (Sigma-Aldrich, rats 1 nmol/2 μl), or Thr4 Gly7-OT (TGOT; Bachem, rats 1 nmol/2 μl; mice 0.5 nmol/2 μl) were infused through an intracerebroventricular infusion cannula (stainless steel, rats: 28 G, 14.7 mm long; mice 25 G, 10.7 mm long) inserted into an indwelling intracerebroventricular guide cannula (rats 23 G, 12 mm long; mice 21 G, 8 mm long), stereotaxically implanted 2 mm above the third (rats) or lateral (mice) ventricle (rats AP: +0.9 mm bregma, ML: −1.2 mm lateral, DV: +5 mm below the surface of the skull, angle: 10°; mice: AP +0.2 mm, ML+1.0 mm, V+1.4 mm; Paxinos and Watson, 1998; Baird et al., 2008; Roy et al., 2011) as described previously (Toth et al., 2012).
The effectiveness of intracerebroventricular OT or TGOT was assessed by measuring the phosphorylation of CaMKII and MEK1/2 in the PVN, which are two known markers of OTR activation (Jurek et al., 2012; and unpublished observation). Although infusion of OT directly into the PVN, or locally released endogenous OT, is anxiolytic, infusion of OT into the lateral ventricle of rats and mice is not (Waldherr and Neumann, 2007; Zoicas et al., 2014). Whether OT is anxiolytic, when infused into the third ventricle, i.e., in the immediate vicinity of the PVN, has not been addressed before. Therefore, we tested whether our infusion procedure affects anxiety-like behavior, using the light/dark box (LDB) test (see below).
Light/dark box
The LDB was performed 10 min after intracerebroventricular infusion of either OT or TGOT as previously described (Slattery and Neumann, 2010; Jurek et al., 2012). Briefly, LDB setup consisted of two compartments; one lit compartment (40 × 50 cm, 350 lux, light box) and one dark compartment (40 × 30 cm, 70 lux). The floors in each compartment were divided into squares (10 × 10 cm) and the compartments were connected via a small opening (7.5 × 7.5 cm). Rats were placed in the light compartment, and measurements of anxiety and locomotor activity (line crosses, time spent in each compartment, rearings, latency to first dark compartment entry, latency to first re-entry in light compartment, and total entries into light compartment) during the 5 min test were assessed online via a camera located above the box, by an observer blind to treatment.
Effects of intracerebroventricular OT or TGOT on stress-induced Crf gene expression and CRTC2/3 translocation
To analyze the influence of OT or TGOT on basal and stress-induced Crf expression and CRTC translocation in the PVN, four groups of rats were used in separate cohorts: (1) intracerebroventricular vehicle (Ringer)/no stress group (Veh-NS), (2) intracerebroventricular OT (TGOT)/no stress group (OT-NS/TGOT-NS), (3) intracerebroventricular vehicle/restraint stress group (Veh-RS), and (4) intracerebroventricular OT (TGOT)/restraint stress group (OT-RS/TGOT-RS). The time course of the experiment is depicted in Figure 1. OT and TGOT are equally effective at the concentration applied (1 nmol in 2 μl; Blume et al., 2008; Lukas et al., 2011; Viviani et al., 2011; Jurek et al., 2012) and have a similar affinity for the OTR (Manning et al., 2008). TGOT was used in some of the experiments, outlined in detail below, to exclude the contribution of the vasopressin receptors to the observed effects as TGOT has minimal affinity for the V1a receptor (Manning et al., 2008). The drugs, or their vehicle, were infused 10 min before the restraint stress procedure, and the animals were returned to their home cage after infusion. Then, rats were restrained for 10, 15, or 30 min in a Plexiglas cylinder (12 cm diameter) with ventilation holes. The nonstressed controls, which remained in their home cage, were time-matched to accommodate for the time of the stressor. Immediately after the respective restraint period, stressed and control rats were decapitated, trunk blood was collected and brains were rapidly removed, frozen on dry-ice, and stored at −80°C until cryosectioning for Crf expression or CRTC translocation analysis. The brains were cut at a thickness of 250 μm, and the PVN was microdissected with a tissue puncher (Fine Science Tools, 1.8 mm diameter; Jurek et al., 2012), placed in 1.5 ml microcentrifuge tubes on dry ice and stored at −80°C until RNA extraction using TriFast Gold (PeqLab) according to the manufacturer's protocol (see below). For CRTC translocation analysis, microdissected PVNs were directly transferred to lysis buffer (see below). In addition, cortical, hippocampal, and septal tissue served as control for brain region-specific effects of stress and OT. The restraint stress protocol for male C57BL/6J mice was similar to that for rats with minor modifications. Mice were restraint using 50 ml Falcon tubes (BD Falcon) with ventilation holes for 10 min, immediately anesthetized, decapitated, and the brains removed for protein isolation from hypothalamic tissue.
Crf/OTR colocalization
Transgenic male OTR-Venus reporter mice were restrained for 30 min, returned to their home cage for 90 min, and then killed by cervical dislocation. This protocol has been previously shown to induce high levels of Crf mRNA (Aguilera and Liu, 2012).
Unstressed and stressed male, as well as pregnant female mice, were deeply anesthetized by tribromoethanol, and the brains from male mice and embryonic day (E)18 embryos from females were immediately dissected out. The brains were postfixed with 4% paraformaldehyde (PFA)/phosphate buffer overnight, and cryoprotected in 30% sucrose. Tissues were embedded in OCT compound (Tissue-Tek) and cut in 30-μm-thick sagittal sections with a cryostat (Leica, CM1950). Sections were individually mounted on slides and processed for in situ hybridization (ISH) to visualize the expression of Crf. Nonradioactive ISH was performed as previously described (Assimacopoulos et al., 2003) with minor modifications. Glass slides with sections were fixed by 4% PFA. Sections were then permeabilized by detergents (1% Nonidet P-40, 1% SDS, 0.5% deoxycholate), and hybridized with digoxigenin-labeled probes (Roche) at 70°C overnight. The cRNA probe sequence was cloned from embryonic day 18 mouse whole brain cDNA by using specific primers (5′-GAGAGAATTCTAGAGCCTGTCTTGTCTGTG-3′ and 5′-GAGACTCGAGAGCATGGGCAATACAAATAA-3′). Excess probes were washed out and sections were blocked with sheep serum and incubated with solution containing alkaline phosphatase-conjugated with digoxigenin-labeled antibodies (Roche). The color was developed with Fastred (Roche) according to the manufacturer's instructions. For double-immunostaining, sections were then incubated with anti-GFP antibody (Santa Cruz Biotechnology, sc-8334, rabbit polyclonal, 1:200), which recognizes the Venus protein, as it is a modified version of GFP (Yoshida et al., 2009), overnight at 4°C. The sections were rinsed and incubated in species-specific secondary antibodies, which were tagged with AlexaFluor 488 (Molecular Probes). Images were taken with a confocal laser microscope (Zeiss LSM780).
Cell culture and stimulation
Primary hypothalamic neurons.
Primary cultures of hypothalamic neurons were prepared from fetal Wistar rats, E18, by collagenase dispersion and plated in six-well plates at 37°C/5% CO2 as previously described (Liu et al., 2008). After 24 h in the presence of serum, cultures were maintained for 8 additional days in neurobasal medium containing B27 supplement (both Life Technologies, Invitrogen). Cytosine arabinoside (Sigma-Aldrich) was added to a final concentration of 5 μm from day 4 onward to prevent glial cell proliferation. On day 10, growth medium was replaced by supplement-free neurobasal medium containing 0.1% BSA. After 1 h of preincubation in this medium, cells were incubated with (FSK+) or without (FSK−) the adenylate cyclase stimulator forskolin (FSK; 1 μm, Sigma-Aldrich) in presence (TGOT+) or in absence (TGOT−) of the specific OTR agonist TGOT (10 nm, Bachem) for the periods of time indicated in Results. Following incubation, cells were harvested and RNA was isolated as described below.
Hypothalamic H32 and amygdaloid Ar5 cells.
The hypothalamic rat cell line H32 (Mugele et al., 1993) and the amygdaloid rat cell line Ar5 (Dalwadi and Uht, 2013) were cultured in DMEM (Life Technologies) containing 10% fetal bovine serum, 10% horse serum, and 1% penicillin/streptomycin (Life Technologies) at 37°C and 5% CO2. Before the experiments, the cells were transferred to 100 mm plates at a density of 3 × 106 cells per plate. Twenty-four hours later, the medium was changed to serum-free medium containing 0.1% BSA for 2 h. To determine the signaling pathways mediating the effects of OT, cells were incubated as described above with (TGOT+) or without (TGOT−) TGOT (10 nm), in the presence (FSK+) or absence (FSK−) of FSK (1 μm H32, 50 μm Ar-5). After incubation for the time periods indicated in the results section and figure legends, cytoplasmic and nuclear proteins were extracted for Western blot analysis for CRTC2/3 and pCREB levels as described below.
Human Be(2)-M17 cells.
The human neuroblastoma cell line Be(2)-M17 (European Collection of Cell Cultures, no. 95011816) was cultured in DMEM/F12 (1:1; Invitrogen) containing l-glutamine, 2.4 g/L sodium bicarbonate, 15% heat-inactivated fetal bovine serum (Invitrogen), 1% nonessential amino acids (Invitrogen), and 0.1 mg/ml gentamycin (Invitrogen 15750-060) at 37°C and 5% CO2. Three days before the experiment, the cells were differentiated to neurons by adding retinoic acid (Sigma-Aldrich) to a final concentration of 5 μm. On the day of experiment, cells were incubated in serum-free DMEM/F12 (+ 0.1% BSA) for 2 h to reduce basal activation of gene transcription initiated by serum components. The TGOT (10 nm), des-Gly-NH2d(CH2)5[Tyr(Me)2Thr4]-OVT (referred to as OTA, 1 μm; Manning et al., 1989), and FSK (50 μm) stimulation and RNA isolation protocols of the Be(2)-M17 cells were similar to those described for H32 cells.
RNA isolation and qPCR
Punched rat PVN tissue from the stress experiments was treated with 1 ml TriFast Gold (PeqLab), vortexed to homogenize tissue, and RNA isolated according to the protocol provided by the manufacturer with some modifications. Briefly, the aqueous supernatant obtained from chloroform precipitation with TriFast was transferred to an RNeasy Mini Kit (Qiagen) column, washed, treated with RNase-free DNase (according to the manufacturer's protocol, Qiagen), washed, and eluted with nuclease-free water. RNA content was determined with the aid of a NanoDrop photospectrometer (PeqLab; Liu and Aguilera, 2009).
To isolate RNA from the stimulated cells, the medium was aspirated off, and 1 ml of TriFast Gold Reagent (PeqLab) was added to the six-well plates (primary cells) or 10 cm culture dishes (Be(2)-M17 cells). The lysed cells were scraped off using a cell scraper and collected in RNase-free 1.5 ml tubes. RNA was isolated as described for brain punches.
Three hundred nanograms of total RNA per sample were used for reverse transcription into cDNA using Super Script III First strand Synthesis System for RT-PCR (Invitrogen). Relative quantification of Crf (NM_000756.1) hnRNA levels was performed using SYBR Green (QuantiFast Qiagen), ribosomal protein L13A (Rpl13A, NR_073024), and cyclophilin A (CycA, NM_021130) as housekeeping genes (Bonefeld et al., 2008; Table 1). Specificity of the qPCR was assured by omitting reverse transcription and by using ddH2O as template. As the results obtained using Rpl13A and CycA yielded similar results, only those for Rpl13A are shown. The PCR protocol consisted of an initial denaturation step of 5 min at 95°C, followed by 50 cycles of denaturation at 95°C for 10 s, and annealing/extension at 60°C for 45 s. At the end of the protocol, a melting curve was generated and PCR products were analyzed by agarose gel electrophoresis to confirm the specificity of the primers. All samples were run in triplicate.
Western blotting
Single hypothalamic punches containing the PVN from the stress experiments were lysed in 500 μl of lysis buffer included in the nuclear extraction kit (Active Motif) that was used to isolate the proteins from both the cytosolic and nuclear fractions. Be(2)-M17 cells were lysed according to Active Motif's protocol for cell lysis. Nuclear extracts from H32 cells and Ar5 cells were prepared using NE-PER Nuclear and Cytoplasmic Extraction Reagent (Pierce) according to the manufacturer's protocol. Western blot analysis was used to detect CRTC2 and 3 translocation and CREB phosphorylation as previously described (Liu et al., 2008), using Lamin A, β-tubulin, Ras, and histone deacetylase 1 as loading controls. For CREB, phosphorylated CREB (pCREB), and CRTC3, 15 μg of cytoplasmic or nuclear extract were loaded and separated in a 6% Tris-glycine gel (Invitrogen). To detect CRTC2 from Be(2)-M17 lysates, 30 μg of protein were loaded onto the gel. Equal amounts of protein from PVN punches were loaded onto SDS gels for analysis (Table 2).
The primary antibodies were diluted in Tris-buffered saline supplemented with 0.1% Tween 20 (TBS-T, Sigma-Aldrich), or 5% fat-free dry milk powder in TBS. After blocking the membrane for 1 h in 5% milk powder or BSA, the blots were incubated with the primary antibodies at the respective dilution overnight at 4°C with gentle agitation. After washing and incubating with the appropriate horseradish peroxidase-labeled secondary antibody, immunoreactive bands were visualized using enhanced chemiluminescence plus detection system and film exposure (Liu et al., 2008), or with the ChemiDoc XRS+ Imager (Bio-Rad). The intensity of the bands on the film was quantified using ImageJ (NIH). Images obtained with the ChemiDoc XRS+ were analyzed with the accompanying software. Following imaging, blots were stripped with a mild stripping solution (Millipore), or 0.2 N NaOH and assayed for HDAC1, Lamin A, total CREB, Ras, or β-tubulin as loading controls.
To control for cross-contamination of cytoplasmic and nuclear fractions, blots were reprobed with nuclear TATA-box binding protein and cytoplasmic anti-Ras antibody. Negligible cross-contamination between fractions was detected.
siRNA transfection
Be(2)-M17 cells (4 × 106 cells) plated in 100 × 20 mm cell culture dishes (Sarstedt) were transfected with 33 nm small interfering RNA (siRNA) and Lipofectamine RNAiMAX Reagent (Invitrogen), diluted in OptiMEM (Invitrogen) to inhibit CRTC2 and CRTC3 synthesis. siRNA specific for CRTC2, CRTC3, and the negative control Silencer no. 5 were purchased from Ambion/Life Technologies. The human oligonucleotide sequences were 5′-CUAUAGUCCUGCCUACUUAtt-3′ (sense) for CRTC2 (NM_181715), and 5′-GCACAUCAAGGUUUCAGCAtt-3′ (sense) for CRTC3 (NM_022769; Table 1). After 12 h, transfection medium was removed and replaced by DMEM/F12 (supplemented with 0.1% BSA, sterile filtered) for 2 h. The cells were then stimulated as described above with (FSK+) or without (FSK−) FSK (50 μm) in the presence (TGOT+) or absence (TGOT−) of TGOT (10 nm) for 1 h, and gene and protein expression were analyzed as described above. Verification of knockdown was performed by RT-qPCR for CRTC2 and CRTC3 mRNA, as well as by Western blot with the respective CRTC antibodies.
Chromatin immunoprecipitation
To investigate whether CRTC3 binds directly to the Crf promoter, Be(2)-M17 cells (7 × 106 cells) were cultured in 75 cm2 flasks (Sarstedt) in growth medium supplemented with 5 μm retinoic acid for 3 d until 80% confluence (∼20 × 106 cells per flask). Cells were then stimulated with FSK and TGOT for 10 min, as described above. Next, the stimulation medium was removed and replaced by 10 ml of 1% formaldehyde fixation solution for 10 min at room temperature. Fixation was stopped by addition of glycine to a final concentration of 125 mm and incubation for 5 min at room temperature. Fixed cells were washed with ice-cold 1× PBS supplemented with 1 mm of the protease inhibitor phenylmethanesulfonylfluoride (PMSF), and harvested with 5 ml of the same solution. The cell suspension was centrifuged for 8 min at 370 × g and 4°C, and the pellet was frozen at −80°C until lysis in cell lysis buffer (10 mm HEPES, 85 mm KCl, 1 mm EDTA, 1% NP-40). Frozen nuclear pellets were thawed, resuspended in 350 μl of nuclear lysis buffer (50 mm Tris/HCl, 1% SDS, 0.5% Empigen BB, 10 mm EDTA, 1 mm PMSF, 1× Protease Inhibitor (Roche) and sonicated five times for 10 s (30 s sample cooling between each sonification step) at output control 2 (Branson Sonifier 250) to produce 0.2–1 kbp DNA fragments of chromatin. Five percent of the sonicated sample was separated and kept as input control. Immunoprecipitation was performed with either 4 μg of CRTC3 antibody (Abcam, ab91654) or rabbit IgG as negative control at 4°C under rotation overnight. DNA-protein-antibody complexes were collected using protein A sepharose beads, washed, eluted, reverse cross-linked with proteinase K, and purified using Qiaquick PCR purification columns (Qiagen). Immunoprecipitated Crf promoter was quantified using quantitative real-time PCR with primers designed to amplify the human Crf promoter region containing the CRE (forward: 5′-AGTCATAAGAAGCCCTTCCA-3′; reverse: 5′-CAACACTGAATCTCACATCCA-3′; Table 1). Fold-change in promoter occupancy is calculated as (percentage input − negative control signal)/(TGOT−/FSK+) group signal.
Statistical analysis
One-way (factor treatment) or two-way (factors treatment × time) ANOVA, followed by the Student Newmann–Keuls test, were performed for statistical analyses of behavioral and molecular experiments. Statistical significance was accepted at p < 0.05. Data are presented as mean ± or + SEM, as indicated in the figure legends. Statistical analyses were performed using SPSS (v19) for windows.
Results
Intracerebroventricular infusion of TGOT induces MEK1/2 and CaMKII phosphorylation in the PVN
Infusion of TGOT (1 nmol/2 μl) into the third ventricle doubled the phosphorylation of both CaMKII (1.9-fold change; p = 0.047) and MEK1/2 (2.1-fold change; p = 0.001) within 10 min in the PVN (Fig. 2A,B), but not in the cortex, hippocampus, nor septum. In confirmation of previous results after intracerebroventricular OT or its antagonist in unstressed male and female rats under basal conditions (Neumann et al., 2000a; Slattery and Neumann, 2010) there were no effects of intracerebroventricular OT or TGOT on anxiety-like behavior in the LDB (data not shown). Thus, intracerebroventricular infusion of OT or TGOT either into the lateral or third ventricle activates the OTR and downstream intracellular signaling pathways in the PVN, without influencing anxiety-related behavior.
OT delays stress-induced Crf gene transcription and plasma ACTH increase
Vehicle-treated/restraint stressed (VEH-RS) rats displayed a rapid increase in Crf hnRNA in the PVN with maximal levels observed between 10 and 15 min of restraint, and returning to near-basal levels by 30 min (Fig. 3A). OT delayed this stress-induced response (OT-RS group), with peak Crf hnRNA values found at 15 min and a slower decrease to basal levels (two-way ANOVA; factors treatment × time F(3,48) = 195.4; p = 0.001; Fig. 3A). After 10 min of stress, OT-RS rats showed a reduced Crf hnRNA level compared with VEH-RS controls (p = 0.001), whereas at the 30 min time point, Crf hnRNA levels of OT-RS rats were higher than both basal levels and those of VEH-RS animals (p = 0.001; Fig. 3A). There was no difference in Crf hnRNA expression within the PVN between OT-NS and VEH-NS rats, indicating that OT has no effect on Crf expression under basal, nonstressed conditions.
As expected, no significant changes in Crf mRNA levels were measured within the time points examined (up to 30 min) in any of the groups (Fig. 3B), consistent with the observation that changes in Crf mRNA levels require at least 1 h to occur (Ma et al., 1997).
We found minor, yet statistically significant, effects of OT treatment on the plasma ACTH response to restraint stress (treatment × time F(3,36) = 3.4; p = 0.032, Fig. 3C). In the VEH-RS group, plasma ACTH concentrations peaked at 10 min and remained elevated until the end of the stressor. In the OT-RS rats, plasma ACTH tended to be lower at 10 min (p = 0.065 vs VEH-RS), and showed a delayed peak 15 min after RS (p = 0.039 vs OT-RS at 10 min), thus, being significantly higher (p = 0.035) than those in VEH-RS animals at 15 min. At 30 min, ACTH levels of Veh-RS and OT-RS-treated rats did not differ (Fig. 3C; p = 0.109). Thus, OT infused intracerebroventricularly before the onset of restraint stress delays the increase of Crf gene transcription following restraint stress, with minor effects on ACTH plasma levels.
OT reduces CRTC3 activation/translocation in the PVN during short exposure (10 min) to restraint stress
To study the mechanism underlying the effects of OT on Crf gene expression in vivo, we examined the translocation of CRTC following 10 min of restraint stress. In VEH-RS rats, stress exposure decreased cytosolic levels of CRTC2 (not significant) and CRTC3 (p = 0.007) compared with VEH-NS rats (Fig. 4A,C). This reduction was reflected in the nuclear compartment, where CRTC2 and CRTC3 levels were both significantly increased (p = 0.03 for both; Fig. 4B,D; n = 7).
Treatment of rats with TGOT before restraint (TGOT-RS) prevented the nuclear translocation of CRTC2 and CRTC3, with CRTC levels in both compartments not being significantly different from those in the Veh-RS control-group (Figs. 4A–D). TGOT infusion in nonstressed animals had no effect either on CRTC2 and CRTC3 levels in the cytosol, or in the nucleus. These data suggest that TGOT inhibits the stress-induced activation, and subsequent translocation, predominantly of CRTC3 and to a lesser extent of CRTC2. Neither restraint stress nor intracerebroventricular TGOT infusion affected CRTC translocation in the cortex, demonstrating the specificity of the effects of stress and TGOT (data not shown; for region-specific CRTC expression levels, see Watts et al., 2011).
As CRTCs are cofactors of CREB signaling, we assessed whether TGOT inhibits the stress-induced CREB phosphorylation. Intracerebroventricular infusion of TGOT did not influence the increase of pCREB following 10 min of restraint stress, suggesting that OT regulates Crf gene expression through modulation of CRTC translocation, rather than CREB phosphorylation (n = 5; Fig. 4E,F).
Together, our in vivo results suggest that OT delays the peak of Crf transcription activity, and that this effect is mediated by the attenuation of CRTC nuclear translocation.
Colocalization of Crf and OTR in mice
To determine whether the effects of TGOT on Crf expression are mediated by OTR expressed in Crf-positive cells, we used OTR-Venus reporter mice (Yoshida et al., 2009). However, the number of GFP/OTR-positive cells, was extremely low in the PVN (as previously reported by Yoshida et al., 2009), making it impossible to study colocalization of Crf and OTR. Nevertheless, in stressed mice there was a substantial augmentation of Crf-expressing cells (Fig. 5A,B), and CRTC3 translocation to the nucleus (Fig. 5J). In contrast, the number of GFP/OTR-positive cells was not increased by the 30 min stress protocol. The number of cells expressing GFP/OTR in the adult amygdala, on the other hand, was high, but these cells did not express Crf (Fig. 5D–F), as previously reported (Gray and Magnuson, 1992; Gray, 1993; Haubensak et al., 2010). We conducted further studies using hypothalamic primary neuron cultures and cell lines shown to coexpress the OTR and Crf to determine whether the observed OT effect on stress-induced Crf expression in vivo is, at least partly, due to direct actions on Crf-positive cells.
Interestingly, a high degree of colocalized Crf and OTR was observed in the amygdala of embryonic OTR-reporter mice, in contrast to adult mice, suggesting developmental changes in the localization and coexpression of Crf and OTR-expressing cells within the amygdala (Fig. 5G–I). Thus, we also included an embryonic amygdala cell line in our investigations.
OT reduces FSK-induced Crf hnRNA and delays CRTC3 translocation in rat hypothalamic and amygdala cells
Stimulation of rat embryonic (E18) primary hypothalamic neurons with FSK (1 μm) and/or TGOT (10 nm) revealed a striking similarity in the pattern of Crf hnRNA expression compared with our in vivo findings from adult rat hypothalamic PVN tissue. FSK induced a rapid increase of Crf hnRNA levels, which peaked at 20–30 min (p < 0.001 for both time points), and returned to near-basal levels (p = 0.061) at 90 min. The peak levels at 20–30 min in the TGOT−/FSK+ group were lowered by the presence of TGOT in the incubation medium (TGOT+/FSK+ group), being significantly reduced to 74 ± 4.6% at 30 min (p = 0.002). Also, TGOT prolonged the duration of FSK-induced Crf gene expression, as Crf hnRNA levels were still elevated at the 90 min time point (Fig. 6A; treatment × time; F(5,84) = 4.6; p < 0.001).
The protein yield of the primary cultures was too low to separate cytosolic and nuclear fractions to analyze CRTC2 and CRTC3 translocation and CREB phosphorylation by Western blotting. Therefore, we made use of the immortalized rat hypothalamic cell line H32 (Mugele et al., 1993), because of its hypothalamic origin and expression of endogenous OTR (Blume et al., 2008). Although these cells have lost Crf gene expression, the CREB–CRTC signaling pathway that is known to regulate Crf gene expression (Liu et al., 2011), and intracellular pathways that are coupled to the OTR, are still intact (Blume et al., 2008). In H32 cells, we found no effect of TGOT on the FSK-induced, transient CREB phosphorylation, peaking at 10 min and returning to baseline 20 min after the incubation started (Fig. 6B). Likewise, the FSK-induced increase of nuclear CRTC2 levels was not affected by TGOT during the first 60 min after administration (Fig. 6C), but a further increase in nuclear CRTC2 was found 90 min after the application of both compounds (p = 0.047).
In sharp contrast, TGOT had profound effects on the FSK-stimulated increase in nuclear CRTC3 levels of H32 cells (Fig. 6D). FSK alone augmented CRTC3 in the nucleus with a peak at 10 min after the onset of FSK stimulation (p < 0.001), and maintained elevated nuclear CRTC3 levels during the 90 min of incubation (p < 0.001). Treatment of H32 cells with TGOT in FSK+ cells reduced the fast peak that was seen in TGOT−/FSK+ cells at 10 min (p = 0.007). Nuclear CRTC3 levels steadily increased until their maximum was reached at 60 min to a level that was ∼1.4-fold higher than the peak levels induced in TGOT−/FSK+-treated cells (treatment × time, F(6,64) = 5.3; p < 0.001). After 60 min, a rapid drop in nuclear CRTC3 levels in TGOT+/FSK+-treated cells annulled the differences between TGOT+/FSK+ and TGOT−/FSK+ cells. Similar effects of TGOT on CRTC3 translocation were observed in the amygdala cell line Ar-5. As experimental evaluation of the optimal dose of FSK in Be(2)M17 cells revealed maximal CRTC3 translocation and CRF transcription at 50 μm FSK (see below; Mulchahey et al., 1999) we chose this concentration also for Ar-5 cells. FSK induced a peak level of nuclear CRTC3 at 30 min, which declined to almost basal levels at 90 min. OT reduced this FSK-induced peak of nuclear CRTC3 significantly (p = 0.009) at 30 min by ∼30% (Fig. 6E; treatment × time, F(12,99) = 3.893; p < 0.001). In all of the experiments, in primary, H32, and Ar-5 cells cells, TGOT+/FSK− was without effect (Fig. 6).
The data from the restraint stress experiments, together with those obtained in cell culture, suggest that OT controls Crf expression through inhibition of nuclear translocation of CRTC3 during stress at early time points, but not under basal (NS or FSK−) conditions.
OT inhibits Crf gene transcription in the human adult cell line Be(2)-M17 via CRTC3
To test the hypothesis that OT exerts its effects on Crf gene transcription via CRTC3, it is necessary to assess Crf hnRNA and nuclear CRTC levels in the same cell type. Therefore, we used the human adult neuroblastoma cell line Be(2)-M17, which differentiates into neurons following exposure to retinoic acid. These cells express Crf and the OTR, as well as the CREB-CRTC signaling pathway, which allowed us to determine effects of TGOT on CRTC translocation, CREB phosphorylation, and Crf gene expression in one system. In initial experiments, lower concentrations (1–10 μm) of FSK, as used for the H32 cells, failed to induce reproducible increases in Crf hnRNA, whereas a 50 μm concentration resulted in a reproducible and constant increase in Crf transcript (for comparable concentrations, see Sala et al., 2000; Heo et al., 2013). Therefore, this concentration was chosen for the rest of the Be(2)-M17 and Ar-5 experiments. Application of TGOT at a concentration of 10 and 100 nm were equally effective in inhibiting the FSK-induced increase of Crf expression, so we used 10 nm throughout the rest of the experiments.
A temporal analysis showed that stimulation of Be(2)-M17 cells with 50 μm of FSK augmented Crf hnRNA levels, which reached their maximum after 60 min and remained elevated until at least 90 min. TGOT affected FSK-induced Crf hnRNA (treatment × time F(3,41) = 4.4; p = 0.01), with lower levels observed at 60 min (p < 0.001) and 90 min (p = 0.005) after the onset of the incubation compared with vehicle/FSK. Incubation with TGOT alone was without effect on Crf hnRNA levels at any of the time points measured (Fig. 7A). The effects of TGOT on FSK-induced Crf transcription were blocked by a specific OTR antagonist (OTA, 1 μm; p < 0.001 vs OTA−/TGOT+/FSK+; Fig. 7B), demonstrating the involvement of OTRs, rather than vasopressin receptors. OTA without FSK or TGOT had no effect on Crf hnRNA levels. Furthermore, the delayed increase of FSK-induced Crf hnRNA levels in the presence of TGOT translated into reduced mRNA synthesis at 90 min (p = 0.008), but increased mRNA levels at time points from 180 min on (p < 0.001; Fig. 7C). Finally, CRTC translocation correlated with the observed effects on Crf transcription. FSK induced a short-lived increase of nuclear CRTC2 levels (p = 0.005), peaking after 10 min of incubation and returning to basal levels shortly thereafter (Fig. 7D), which is in contrast to H32 cells, where the response lasts for up to at least 90 min. In the presence of TGOT, FSK (TGOT+/FSK+) failed to augment nuclear CRTC2, however, no significant difference between the TGOT−/FSK+ and the TGOT+/FSK+ group was found (treatment × time, F(3,31) = 0.47; p = 0.707).
As observed with CRTC2, FSK induced a peak level of CRTC3 at 10 min of treatment with elevated, submaximal levels at 30 and 60 min of treatment (treatment × time; F(4,52) = 2,56; p = 0.05). In contrast, in the TGOT+/FSK+ group this peak of nuclear CRTC3 at 10 min was absent (p = 0.016), however, levels steadily increased until at least 60 min of treatment, being no longer different from the TGOT−/FSK+ group at the 30 and 60 min time points (Fig. 7E).
The FSK-induced phosphorylation of CREB was not affected by TGOT (Fig. 7G), again suggesting that OT does not control Crf expression by modulation of CREB phosphorylation.
Thus, the effects of TGOT on FSK-induced changes in Crf gene expression seen in the human neuroblastoma cell line Be(2)-M17 as measured by hnRNA levels, mRNA levels, CRTC2 and CRTC3 translocation, and CREB phosphorylation are strikingly similar to those observed in the rat hypothalamic and amygdala cells, both in vivo and in vitro, as well as in the mouse PVN in vivo. Therefore, we concluded that the Be(2)-M17 cells are a suitable model to further study the modulatory role of OT on CRTC translocation and the control of Crf gene expression.
CRTC3 is necessary for the inhibition of Crf gene expression by OT
To determine whether CRTC2 and/or CRTC3 are necessary for the inhibition of Crf gene expression by OT, we downregulated the expression of CRTC2 and CRTC3 by transfecting the Be(2)-M17 cells with siRNA oligonucleotides. The CRTC2 siRNA construct downregulated both CRTC2 and CRTC3 mRNA levels by 60% (p = 0.002) and 20% (p = 0.019), respectively (F(3,15) = 20.45; p < 0.001; Table 3). The CRTC2 mRNA downregulation was mirrored at the protein level (64%), whereas CRTC3 protein levels did not change. The CRTC3 siRNA construct specifically reduced both CRTC3 mRNA (F(3,15) = 3.54; p = 0.048) and protein levels by 60% (p = 0.017). The nonspecific control siRNA oligonucleotide was without effect (Table 3).
To assess whether the inhibitory effect of TGOT on FSK-induced Crf transcription depends on CRTC2 or CRTC3, we stimulated the cells for 60 min with FSK, when Crf hnRNA levels were maximal (Fig. 8A). Statistical analysis revealed a significant effect of TGOT+/FSK+ treatment (F(3,27) = 33.93; p < 0.001; Fig. 8A), and of pretreatment with siRNA constructs (F(3,27) = 17.4; p < 0.001). Student–Newman–Keuls post hoc analyses confirmed the significant reduction in FSK+-stimulated Crf hnRNA levels by TGOT in the three groups pretreated with either vehicle (p < 0.001), control oligo (p = 0.018), or siCRTC2 (p = 0.005). Cells that were pretreated with the CRTC3 siRNA construct revealed a less pronounced FSK-induced rise in Crf hnRNA (50% of vehicle or control oligo pretreated cells, p < 0.001). Importantly, TGOT treatment was without effect on Crf transcript levels in CRTC3 siRNA transfected cells. We exclude a floor effect, as levels of both, TGOT−/FSK+ and TGOT+/FSK+, remained three times higher than TGOT−/FSK− controls (this group was set to 1; Fig. 8A, dashed line). This indicates that the effect of TGOT on FSK-induced Crf gene expression depends exclusively on CRTC3. The lack of effect of pretreatment with the CRTC2 siRNA construct might be explained by the observation that nuclear CRTC2 levels have already returned to basal levels after 60 min of FSK incubation (Fig. 7C).
To demonstrate the unique role of CRTC3 in the reduction of FSK-induced Crf expression by TGOT further, we studied the binding of CRTC2 and CRTC3 to the Crf promoter after 10 min of FSK incubation. At this time point, the FSK-induced increase of both CRTC2 and CRTC3 nuclear levels are maximal (Figs. 7C,D). The chromatin immunoprecipitation analyses showed that binding of CRTC2 and CRTC3 to the Crf promoter was induced by FSK (Fig. 8B,C; one-way ANOVA, F(3,15) = 10.7; p = 0.001 for CRTC2; F(3,19) = 4.4; p = 0.02 for CRTC3), and that TGOT prevented CRTC3 (p = 0.028), but not CRTC2 binding to the Crf promoter. Again, TGOT+/FSK− was without effect. This supports the hypothesis that OT mediates the downregulation of Crf transcription during stress via CRTC3, but not CRTC2.
Discussion
The current study describes a novel molecular mechanism recruited by OT to regulate Crf expression during acute stress exposure. Central to this mechanism is the inhibition, or delay, of CRTC3 translocation from the cytosol to the nucleus, resulting in reduced CRTC3 binding to the Crf promoter and Crf gene expression. The control of Crf expression is a completely novel role of CRTC3 in the brain, which is only active during stress, not under basal (stress-free) conditions. The present data show that it operates within the PVN of both rats and mice in vivo, and in cells of human and rat (Wistar and Sprague-Dawley strains) origin, suggesting that this mechanism is conserved across mammalian species. Furthermore, our in vitro data show that prolonged exposure to FSK and TGOT leads to a potentiation of Crf expression. Thus, the early inhibition and the late potentiation of Crf expression by CRTC3 translocation and promoter binding might explain the contradictory results of previous studies describing either inhibitory (Windle et al., 2004) or potentiating effects (Petersson et al., 1999) of OT on ACTH release.
The present in vitro experiments clearly show that OT can directly regulate Crf transcription in Crf-expressing neurons. Importantly, the kinetics of the cellular responses in vitro parallel those we observed in vivo, which is consistent with a direct effect of OT on Crf-expressing cells. However, the inability to demonstrate colocalization of the OTR and Crf in the PVN raises questions as to whether the in vitro findings do indeed apply to the mechanism by which OTR activation, and subsequent CRTC3 translocation, reduce Crf expression in vivo. Lack of visualization of the OTR within the PVN is likely to be the consequence of low expression levels, consistent with a previous report where OTR was detected in the PVN only following infusion of an OTR antagonist, suggesting that OTR expression is under negative control of its ligand (Freund-Mercier et al., 1994). Even the use of the reporter mice in our study could not reliably reveal the expression of OTRs in the PVN as previously described (Yoshida et al., 2009), although it is known from behavioral (this study; Blume et al., 2008; Jurek et al., 2012), physiological (Moos and Richard, 1989; Leng et al., 2008; Jurek et al., 2012; van den Burg et al., 2015), and molecular studies (Freund-Mercier et al., 1994; Dabrowska et al., 2011) that the OTR is expressed in the PVN. Consistent with this, a number of Crf cells are surrounded by OT-positive boutons, suggesting that these Crf cells express the OTR (Dabrowska et al., 2011). Moreover, Crf expression can be induced during stress in magnocellular neurons that are Crf-negative under control conditions (Kresse et al., 2001). This has been shown for stressors of a systemic nature, especially inflammation (Kresse et al., 2001) and adrenalectomy (Swanson et al., 1983). Thus, it is conceivable that the effects observed in vivo are mediated, at least in part, by a small number of OTRs expressed in Crf neurons, undetectable in the present experimental conditions.
On the other hand, it is also possible that the effects in vivo are indirect, mediated by OTR located in afferent projections to Crf neurons in the PVN. It has been shown that the inhibitory effect on Crf expression following intracerebroventricular administration of OT can be blocked by the GABAA receptor blocker bicuculline (Bülbül et al., 2011), suggesting a role for OT-sensitive inhibitory interneurons. However, this is unlikely because none of the OTR-expressing neurons in the PVN appear to be GABAergic (Dabrowska et al., 2013). Thus, although direct and indirect mechanisms may be involved, our data clearly demonstrate that OT administration alters the temporal dynamics of CRTC3 translocation, and subsequent Crf transcription under stress conditions; both in vivo and in vitro. Moreover, the ability of CRTC3 siRNA to block the effects of OT on Crf hnRNA in Be(2)-M17 cells indicates that the modulatory action of OT on Crf transcription involves the CRTC3 pathway.
A rather unexpected finding was that OT only influences the translocation of CRTC3 and not that of CRTC2. Although we did observe nuclear translocation of CRTC2 following 10 min of restraint, OT did not affect nuclear CRTC2 levels. Earlier studies have shown that, unlike CRTC1, both CRTC2 and CRTC3 control the transcription of the Crf gene (Liu et al., 2011, 2012). However, these studies focused on the recovery period of rats following 1 h of restraint stress, a stress-protocol that has been shown to reduce OTR content in the PVN of voles (Smith and Wang, 2013). This finding led us to study the effects of OT on Crf expression during a short acute stress paradigm, when the OTRs are still fully expressed and functional.
The actions of OT on cAMP/CREB-regulated Crf gene expression were limited to CRTC3 regulation, as shown by the lack of effect of OT on stress-induced CREB Ser133 phosphorylation. Both our in vivo and in vitro experiments showed a similar short-lived, stress- or FSK-induced increase of CREB phosphorylation in the absence, as well as in the presence, of OT or the specific OTR agonist TGOT. This is consistent with the notion that the actions of CRTCs are independent from the CREB Ser133 phosphorylation site (Takemori et al., 2007).
The lack of effect of OT on pCREB also makes the involvement of regulatory mechanisms related to CREB activity unlikely, such as competition of the inhibitory factor ICER with pCREB at the Crf promoter. Furthermore, ICER is Gαs-cAMP activated, whereas the OTR is coupled to Gαi and Gαq, with PKC/PLCβ as downstream effectors, rather than cAMP (Gimpl and Fahrenholz, 2001).
The intracellular signaling pathway that couples the OTR to CRTC3 is currently not known. A major protein kinase known to regulate CRTC phosphorylation and nuclear trafficking is salt-inducible kinase (SIK; Sasaki et al., 2011; Clark et al., 2012). SIK is a member of the mammalian AMP-activated protein kinase (AMPK) family. AMPK is activated by OT in skeletal muscle cells (Lee et al., 2008; Florian et al., 2010), and it is thus possible that the OTR also activates SIK. Its two isoforms, SIK1 and SIK2, are present in hypothalamic neurons of the PVN and have been shown to inhibit Crf transcription by impairing CRTC trafficking to the nucleus (Liu et al., 2011). An additional target of SIK1 may be the myocyte enhancer factor 2 (MEF2), which also has a binding site in the Crf promoter and regulates its gene transcription following activation by CaMKII, MEK1/2, ERK5/BMK, and p38 (Zhao et al., 1999; Lu et al., 2000; Flavell et al., 2006). These kinases are coupled to the OTR (this study; Devost et al., 2008; Kim et al., 2015), whereas other target genes of MEF2, such as RGS2 or PACAP (Flavell et al., 2008) have also been implicated in the anxiolytic effect of OT (Park et al., 2002; Okimoto et al., 2012), and regulation of Crf expression (Stroth and Eiden, 2010). Whether the OTR exclusively regulates the CRTC3 pathway via SIK1, or has also implications on the activity of MEF2, again via SIK1, to alter the expression of Crf and Crf-related proteins, remains to be elucidated.
In our study, acute synthetic OT exerts only a minor influence on acute HPA-axis activity, as it had only minimal effects on the increase of plasma ACTH during stress. It is thus unlikely that OT modulates Crf release, leaving the acute stress response intact. The observed actions of OT on Crf gene transcription during the early phase of the stress response might rather serve to reallocate energy away from costly investments like gene expression. Replenishment of depleted CRF stores (Aguilera and Liu, 2012) could be enhanced by OT later, as suggested by increased Crf transcription, at both the hnRNA and mRNA levels, and CRTC3 translocation at later time points in our in vitro studies. A repression of Crf gene transcription as part of a negative feedback mechanism is unlikely, because, in our protocols, OT is already high before the onset of the stressor (such as following successful mating in male rats; Waldherr and Neumann, 2007) and Crf transcription is only delayed, not repressed, over the whole 30 min period of restraint. Furthermore, there is no corticosterone in the cell culture medium that could be implicated in the effects of OT on the CRTC3–Crf expression pathway. However, the lack of negative feedback in vitro could account for the late (4–6 h) potentiation of FSK-induced Crf gene transcription by TGOT in vitro.
Our study has revealed a novel, and specific, role for the CREB cofactor CRTC3 and is the first study that describes a clearly defined physiological role of CRTC3 in the brain; namely regulation of stress-induced Crf transcription. Importantly, CRTC3 is recruited during stress, but, in sharp contrast to CRTC2, its translocation from the cytosol to the nucleus is regulated by the OTR. This particular stress-CRTC3 pathway, which can be modulated by OT, may not be unique to the Crf gene alone, but to other CREB-regulated genes, such as Crfr1/2 and Crf-bp as well. Thus, a further search for stress-related genes that are under the control of the OT-CRTC3 pathway seems warranted and might lead to a better understanding of the regulation of gene transcription during stress.
Footnotes
↵*I.D.N. and E.H.v.d.B. share last authorship.
This work was supported by the German Research Foundation (I.D.N., E.v.d.B.), the Boehringer Ingelheim Foundation (B.J.), the Intramural Research Program, NICHD (G.A.), and SRPBS, Japan (MEXT; Y.H., K.N.). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the paper.
The authors declare no competing financial interests.
- Correspondence should be addressed to Inga D. Neumann, University of Regensburg, Regensburg 93040, Germany. inga.neumann{at}biologie.uni-regensburg.de