A randomized trial shows dose-frequency and genotype may determine the therapeutic efficacy of intranasal oxytocin
==================================================================================================================

* Juan Kou
* Yingying Zhang
* Feng Zhou
* Cornelia Sindermann
* Christian Montag
* Benjamin Becker
* Keith M Kendrick

## Abstract

**Background** The neuropeptide oxytocin is proposed as a promising therapy for social dysfunction by modulating amygdala-mediated social-emotional behavior. Although clinical trials report some benefits of chronic treatment it is unclear whether efficacy may be influenced by dose frequency or genotype.

**Methods** In a randomized, double blind, placebo-controlled pharmaco-fMRI trial (150 male subjects) we investigated acute and different chronic (every day or on alternate days for 5 days) intranasal oxytocin (24IU) effects and oxytocin receptor genotype-mediated treatment sensitivity on amygdala responses to face emotions. We also investigated similar effects on resting state functional connectivity between the amygdala and prefrontal cortex.

**Results** A single dose of oxytocin reduced amygdala responses to all face emotions but for threatening (fear and anger) and happy faces this effect was abolished after daily doses for 5 days but maintained by doses given every other day. The latter dose regime also enhanced associated anxious-arousal attenuation for fear faces. Oxytocin effects on reducing amygdala responses to face emotions only occurred in AA homozygotes of rs53576 and A carriers of rs2254298. The effects of oxytocin on resting state functional connectivity were not influenced by either dose-frequency or receptor genotype.

**Conclusions** Infrequent chronic oxytocin administration may be therapeutically most efficient and its anxiolytic neural and behavioral actions are highly genotype-dependent in males.

## Introduction

Intranasal oxytocin (OT) has been proposed as novel treatment to attenuate social dysfunctions and anxiety in autism spectrum disorder, social anxiety and schizophrenia (Kendrick et al., 2017; Meyer-Lindenberg et al., 2011; Young and Barrett, 2015) as well as an augmentative strategy to facilitate fear extinction (Eckstein et al., 2015) and working memory performance (Zhao et al., 2018). Following an initial landmark study in humans demonstrating that OT reduces amygdala responses to threatening (fear and angry) face stimuli (Kirsch et al., 2005), subsequent preclinical research has consistently reported that its effects on amygdala function (attenuated fear reactivity and increased amygdala-prefrontal intrinsic connectivity - Domes et al., 2007; Eckstein et al., 2017; Kendrick et al., 2017; Meyer-Lindenberg et al., 2011; Quintana et al., 2016; Spengler et al., 2017; Sripada et al., 2013; Wang et al., 2017) may represent its primary therapeutic-relevant neural mechanism of action.

The proposed neural mechanisms of OT have been established using acute single dose administration protocols, with a dose-response study demonstrating 24 International Units (IU) as optimal for reducing amygdala responses to fear stimuli using conventional intranasal spray devices (Spengler et al., 2017). To date however, no human studies have demonstrated whether chronic treatment protocols can produce enhanced functional effects and initial clinical trials using chronic (twice daily) administration over periods of a week or longer have reported inconsistent, or at best modest therapeutic efficacy on social-emotional dysfunctions (Guastella et al., 2015; Halverson et al., 2019; Jarskog et al., 2017; Kendrick et al., 2017; Watanabe et al., 2015; Yamasue et al., 2018).

G-protein-coupled neuropeptide receptors typically exhibit internalization following even a relatively short period of constant exposure to their target peptides and cannot subsequently respond (desensitized) until they are recycled back the surface of the cell membrane (Pierce et al., 2002; Lohse and Hofmann, 2015). Growing evidence from both in vitro and animal studies suggest there may be extensive internalization and subsequent desensitization of OT receptors (OXTR) following repeated OT administration (Smith et al., 2006; Stoop, 2012). In terms of neural OXTR receptors those in the amygdala may be particularly susceptible to desensitization (Terenzi and Ingram, 2005) and this region critically mediates many of the effects of intranasal OT on social cognition and anxiety (Kendrick et al., 2017). Moreover, rodent models have reported that in contrast to single doses, chronic administration of OT can actually produce social impairment (Bales et al., 2013; Du et al., 2017; Huang et al., 2014) in the context of reduced receptor expression in the amygdala and nucleus accumbens (Du et al., 2017). Furthermore, an initial study reported that chronic versus single administration of OT results in divergent neurochemical changes in both the rodent and human frontal cortex (Benner et al., 2018), emphasizing high translational relevance of the preclinical animal models for clinical trials employing chronic treatment protocols. Thus, empirically evaluated optimal treatment protocols for chronic intranasal OT administration in humans are urgently needed to determine its therapeutic potential.

A further unresolved issue is that different OT receptor (OXTR) genotypes have been associated with social behaviors and may influence sensitivity of behavioral and neural responses to intranasal OT. In particular, OXTR polymorphisms rs53576 and rs2254298 have been associated with autism (Cataldo et al., 2018), deficits in social and emotional processing and anxiety (Jurek and Neumann, 2018; Parker et al., 2014; Yang et al., 2017) and individual variations in behavioral and neural responses to intranasal OT (Chen et al., 2015; Feng et al., 2015). Indeed, variants of rs2254298 may represent a trans-diagnostic biomarker for social dysfunctions (Brüne, 2012). It is therefore imperative to establish whether OXTR polymorphisms could influence sensitivity to effects of intranasal OT.

While single doses of OT can produce numerous functional effects, two of key importance are attenuated amygdala responses to face emotions, particularly fear, and increased resting state functional connectivity between the amygdala and prefrontal cortex (Domes et al., 2007; Eckstein et al., 2017; Kirsch et al., 2005; Koch et al., 2016; Quintana et al., 2016; Spengler et al., 2017; Sripada et al., 2013). While most task studies have focused on the effects of OT on amygdala responses to fear faces using a variety of different task paradigms some studies have reported valence independent effects of OT (Domes et al., 2007; Koch et al., 2016; Quintana et al., 2016). Both these task and resting state neural markers have been primarily associated with attenuated anxiety in terms of reduced responses to social threat and enhanced top-down control of emotion (Zhao et al., 2019), although effects of OT on the amygdala may influence a range of social cognition domains (Kendrick et al., 2017; Meyer-Lindenberg et al., 2011; Young and Barrett, 2015).

To help determine optimal treatment protocols for chronic intranasal OT administration the current pre-registered double-blind, randomized between subject placebo (PLC)-controlled pharmacological neuroimaging trial therefore aimed to replicate previous findings for effects of a single dose of OT (24IU) on reduced amygdala responses to fear or other face emotions (Domes et al., 2007; Kirsch et al., 2005; Koch et al., 2016; Quintana et al., 2016; Spengler et al., 2017; Wang et al., 2017) in a large sample at the whole brain level and on increasing resting-state functional connectivity between the amygdala and prefrontal cortex (Eckstein et al., 2017; Sripada et al., 2013). Next, the primary outcome of the trial was to establish whether different repeat dosing protocols (daily, OT5 group, or every other day, OT3 group, for 5 days) modulate amygdala-centered neural and behavioral responses and resting state functional connectivity. As a secondary outcome, we further aimed to establish if effects of OT on primary outcome measures were influenced by OXTR genotype, focusing on the rs53576 and rs2254298 polymorphisms most frequently associated with social function and individual variations in intranasal OT sensitivity.

## Material and Methods

The main study objectives were to investigate both the effects of acute (single dose) and repeated doses (daily or every other day for 5 days) of 24IU OT versus PLC on two primary outcomes: (1) amygdala and behavioral responses to face emotions (particularly fear faces) and (2) resting state functional connectivity between the amygdala and frontal cortex. We investigated modulatory influences of OXTR rs53576 and rs2254298 polymorphisms on the sensitivity to intranasal OT as a secondary outcome.

### Participants

A total of 150, right-handed healthy adult male subjects were enrolled according to common inclusion and exclusion criteria for human OT-administration studies (see **SI**). Given that methods for estimating power for whole brain task-based fMRI analyses are limited (Poldrack et al., 2017), subject number was initially chosen on the basis of recent large-scale studies investigating dose effects of OT (Spengler et al., 2017 n = 111 males; Lieberz et al., 2019, 90 females). A priori power estimates for behavioral analyses were calculated using Webpower (Zhang and Yuan, 2018) giving a range of 60-80% power for ANOVA main and interaction effects and post-hoc tests. For full details see SI. Calculating power for fMRI experiments is more problematic and likely to be much lower, although since the amygdala was our only pre-registered ROI of interest, and the majority of fMRI post-hoc analyses used a single amygdala ROI, this avoids the problem of significantly reduced power when undirected whole brain analyses are performed (see Cremers et al., 2017). The current chosen sample size is also comparable to the previous largest fMRI studies investigating task-related effects of OT (Spengler et al., 2017; Lieberz et al., 2019). Three subjects failed to meet inclusion criteria (were currently taking medications) and 9 subjects were excluded due to failure to complete the study or excessive head movement (see **Fig. S1**). Subjects were randomly (computer-generated randomization) assigned to repeated intranasal treatment (single daily dose on five consecutive days) of (1) placebo (PLC; n = 46, M ± SD, 22.46± 2.3 years), (2) oxytocin (OT5; n = 49, M ± SD, 21.78 ± 2.3 years), or interleaved OT and PLC (OT on days 1, 3, 5, PLC on days 2, 4; OT3; n = 43, M ± SD, 21.02 ± 2.0 years) (see **SI**). To ensure compliance all subjects were supervised during self-administration of all nasal sprays. Groups did not differ on a number of potential confounders (see **SI** and **Table S4**). The study was approved by the local ethics committee (Institutional Review Board, University of Electronic Science and Technology of China) and subjects provided written informed consent. The study was in accordance with the latest revision of the Declaration of Helsinki, pre-registered at Clinical Trials.gov ([NCT03610919](http://medrxiv.org/lookup/external-ref?link_type=CLINTRIALGOV&access_num=NCT03610919&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) - [https://clinicaltrials.gov/ct2/show/NCT03610919](https://clinicaltrials.gov/ct2/show/NCT03610919)) and in line with recommendations for trials in psychological experiments (Guidi et al., 2018) (see **Fig. S1** for Consort flow diagram). Subjects received 500 RMB monetary compensation for participation.

### Experimental procedures

The study employed a double-blind, randomized, placebo-controlled, between-subject design. The OT and PLC sprays used in the 3 groups were supplied by Sichuan Meike Pharmaceutical Co. Ltd, Sichuan, China in identical dispenser bottles containing identical ingredients (glycerine and sodium chloride) other than OT. In line with recommended guidelines experiments started 45 minutes after intranasal administration (Guastella et al., 2013). Recent studies have confirmed that intranasal OT enters the brain (Lee et al 2018), increases cerebrospinal fluid and blood concentrations (Striepens et al., 2013) and produces activity changes in brain regions expressing OXTRs in humans at around 45 minutes after administration (Paloyelis et al., 2016; Quintana et al 2019). In post-treatment interviews on days 1 and 5 subjects could not identify better than chance whether they had received OT or PLC (χ2 < 0.1, ps > 0.2, **Table S4**). For OXTR genotyping subjects provided buccal swaps on the 1st day for analysis of OXTR rs2254298, rs53576 SNPs (see **SI** and (Montag et al., 2017)).

For the implicit face-emotion processing task 208 grayscale facial stimuli displaying happy, neutral, angry or fearful facial expressions (n = 26 per category, 50% female) were used (two different stimulus sets - see **SI** for details).

### Primary outcomes and analysis plan

Face emotion-related amygdala responses were assessed using the event-related implicit face processing fMRI paradigm on treatment days 1 and 5. Intrinsic amygdala connectivity was assessed by means of resting state fMRI assessments before the task-paradigm. Valence, arousal and intensity ratings (scale: 1-9) for the facial stimuli were collected immediately after MRI acquisition as additional behavioral outcomes (further details see **SI**). In line with the main aim of the study changes in amygdala face emotion reactivity and amygdala intrinsic connectivity between the 1st and 5th day served as primary outcome measures. Changes in valence, arousal and intensity ratings for the different face emotions represented an associated behavioral outcome measure.

### Imaging acquisition and analysis

Blood oxygenation level-dependent contrast functional images were acquired using standard sequences on a 3T GE MR750 system. High resolution T1-weighted structural MRI data were acquired on both 1st and 5th days to improve normalization. MRI data was preprocessed using validated procedures in SPM12 (Friston et al., 1994) (Statistical Parametric Mapping; [http://www.fil.ion.ucl.ac.uk/spm](http://www.fil.ion.ucl.ac.uk/spm)) and Data Processing Assistant or Resting-State fMRI (Yan, 2010) (DPARSFA; [http://rfmri.org/DPARSF](http://rfmri.org/DPARSF)) (for details see **SI**). First level General Linear Models (GLM) for the task-related fMRI data included separate regressors for the four emotional conditions, gender identity rating period and 6 movement parameters and appropriate contrasts were subjected to a second level random effects analysis.

### Primary outcomes: The effects of OT on amygdala and behavioral face emotion reactivity and intrinsic connectivity

Differences in amygdala and behavioral responses to face emotions were examined using mixed ANOVAs with treatment group (PLC, OT3, OT5) as between-subject and time point (1st and 5th days) and face emotion (angry, fear, happy and neutral) as within-subject factors. Dependent variables were extracted target region amygdala responses using the template of left and right amygdala masks from the Brainnetome atlas (Fan et al., 2016) to faces (parameter estimates extracted using Marsbar, [http://marsbar.sourceforge.net](http://marsbar.sourceforge.net)) and valence, intensity and arousal ratings. Post-hoc comparisons explored significant interactions for amygdala and behavioral ratings. Associations between amygdala responses and behavioral ratings were assessed using Pearson correlation.

Differences in resting state amygdala functional connectivity were determined using a voxel-wise seed-to-whole-brain ANOVA as implemented in the SPM flexible factorial design with treatment (PLC, OT3, OT5) as between-subject and time point (1st and 5th days) as within-subject factors and amygdala connectivity maps as dependent variable. (See data analysis flow chart in **SI**)

### Secondary outcome measures: influence of OXTR genotype

For rs2254298 subjects were divided into A carriers (AA and AG) and A non-carriers (GG) and rs53576 into G carriers (GG and GA) and G non-carriers (AA). To explore effects of genotype as a secondary outcome measure, OXTR group was included as an additional between-subject factor in the corresponding ANOVAs and amygdala reactivity to face emotions, amygdala-prefrontal cortex intrinsic functional connectivity and behavior rating response served as dependent variables.

### Control for treatment effects on brain structure

To control for potential confounding effects of single- and repeated OT-administration on brain structure, a voxel-based morphometry (VBM) analysis was conducted on the T1-weighted images acquired on both testing days using SPM12 standard procedures (Ashburner, 2007; Ashburner and Friston, 2005)(see **SI**).

### Thresholding and statistical implementation

In line with our regional hypothesis the analysis of the implicit emotional face task focused on the amygdala. To account for the complex design of the trial we employed two analytic approaches: (1) A region of interest analysis (ROI) on the extracted beta estimates from atlas-based masks for the left and right amygdala (Brainnetome Atlas, Fan et al., 2016) respectively. To this end the extracted estimates were subjected to three-way ANOVAs (treatment ×time point ×face emotion) examining effects on the left and right amygdala separately. The p-value for the analysis was adopted to p<0.025 (Bonferroni corrected multiple comparison 0.05/2). (2) To facilitate a more specific localization of the effects and to determine a more specific seed mask for the seed-to-whole brain resting state functional connectivity analysis a voxel-wise small volume corrected (SVC) analysis for the entire bilateral amygdala was additionally employed. To this end a peak-level Family-wise error (FWE) small volume correction was applied to a single mask encompassing the bilateral amygdala as provided in the Brainnetome Atlas (Fan et al., 2016) (total volume of the mask, 4572.27mm3) (FWE p < .05). For the seed-to-whole brain resting state data analysis a whole-brain cluster-based FWE-correction with an initial cluster forming threshold of p < .001 and FWE p < .05 was employed. Behavioral and neural (extracted parameter estimates) indices were further analyzed using SPSS 22.0 with appropriate ANOVA models and Bonferroni corrected post-hoc tests. P < 0.05 two-tailed was considered significant. Partial eta squared (*η*2p) and Cohen’s *d* were computed for behavioral results as measures of effect size and SPSS outlier analysis was used to confirm the absence of statistical outliers. Full details of the post-hoc and correlation analyses are provided in the SI (Tables S7 - S11).

## Results

### Primary outcome measure: ROI analysis of acute and repeated doses on OT-evoked changes in amygdala responses to face emotions

Examination of the extracted parameter estimates from the amygdala (a priori ROI approach in accordance with trial pre-registration) revealed a significant three-way ANOVA for the right amygdala (F(6,405)=2.563; p<0.019), but not the left amygdala (F(6,405)=1.860; p=0.087). Post-hoc analyses on the extracted beta values from the right amygdala mask on the 1st day revealed that both OT groups compared to PLC exhibited a suppression of right amygdala reactivity (fear: ps < 0.015, angry: ps < 0.016; happy: ps < 0.041: neutral: ps < 0.017). Analysis on the 5th day revealed that the OT3 group exhibited a suppression of right amygdala reactivity to fear and angry faces relative to both the PLC (fear: p = 0.003, angry: p = 0.008) and OT5 groups (fear: p = 0.021; angry: p = 0.082) but the OT5 group did not (ps > 0.315, relative to PLC group). The OT3 group also exhibited reduced right amygdala reactivity to happy faces on the 5th day compared with the PLC group (p = 0.028) but there was no significant difference between OT3 and OT5 groups (p = 0.272). There were no significant effects of OT on the 5th day for neutral faces (ps>0.203, although OT3 vs PLC was marginal, p = 0.052). While amygdala reactivity to fear faces did not differ in the PLC group between days 1 and 5 (p = 0.356), its reduction in the OT3 group on day 5 was equivalent to day 1 (p = 0.538) whereas in the OT5 group the reduction on day 5 was significantly less than on day 1 (p = 0.035). Amygdala reactivity to angry faces decreased significantly more in the OT3 group on day 5 relative to day 1 (p = 0.006), although it also did so in the PLC group for both angry (p = 0.007) and neutral (p = 0.028) (**Fig. 1b-e**).

![Fig. 1](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/08/25/2020.08.23.20180489/F1.medium.gif)

[Fig. 1](http://medrxiv.org/content/early/2020/08/25/2020.08.23.20180489/F1)

Fig. 1 
The effect of intranasal oxytocin (OT) treatment on neural responses to faces. (a) The mask of right amygdala from the Brainnetome atlas (Fan et al., 2016)(b-e) Post hoc analysis on extracted parameter estimates from the right amygdala mask for all face emotion categories on the 1st and 5th days.# p < 0.1, * p < 0.05, ** p < 0.01, two-tailed. Bars depict standard errors.

### Primary outcome measure: behavioral rating analysis of acute and repeated doses to face emotions

Mixed ANOVAs on the behavioral data revealed a significant main effect of treatment group on arousal ratings (arousal: F2, 135 = 5.06, p = 0.008, *η*2p = 0.07) due to ratings generally being lower in the OT3 group. There was also a marginal significant treatment group x time point x face emotion interaction effect for arousal (F6,405 = 1.83, p = 0.093, *η*2p = 0.03) and a significant one for intensity (F6,405 = 3.327, p = 0.004, *η*2p = 0.05). Post-hoc analysis showed that intensity ratings were significantly reduced for fear faces in the OT3 group on the 5th day (fear: p = 0.004 versus PLC, *d* = 0.62, 95% CI, −1.299 to −0.247, p = 0.088 versus OT5, *d* = 0.33, 95% CI, −0.968 to 0.068). Exploratory post-hoc analysis of arousal ratings revealed a similar reduction to fear faces on day 5 (p = 0.001 versus PLC, *d* = 0.77, 95% CI, −1.544 to −0.420, p = 0.033 versus OT5, *d* = 0.42, 95% CI, −1.158 to −0.051). Arousal ratings to fear faces were also decreased in the OT3 group on the 1st day compared to the PLC group but not the OT5 group (arousal: p = 0.015 versus PLC, *d* = 0.60, 95% CI, 0.129 to 1.175, p = 0.355 versus OT5, 95% CI, −0.268 to 0.744). Arousal and intensity ratings to angry faces were reduced on day 5 compared with day 1 in the OT3 (arousal: p = 0.001, *d* = 0.30, 95% CI, 0.164 to 0.619, intensity: p = 0.001, *d* = 0.32, 95% CI, 0.173 to 0.641) and OT5 groups (arousal: p = 0.023, *d* = 0.19, 95% CI, 0.032 to 0.461, intensity: p = 0.006, *d* = 0.26, 95% CI, 0.088 to 0.526). There were no significant differences for angry, neutral and happy faces between OT and PLC groups on the 1st day (ps >0.194, although OT5 vs PLC was marginal, p=0.051 for neutral faces) and 5th day (ps >0.671, although OT3 vs PLC p=0.053 was marginal for angry faces). There was however a significant group difference between OT3 and OT5 on the 1st and 5th day for happy and neutral faces (ps<0.007) (**Fig.2**). While there was a significant 3-way interaction for valence ratings (F6, 405 = 2.63, p = 0.016, *η*2p = 0.04), post-hoc analysis revealed no significant differences between OT and PLC groups on the 1st (ps > 0.200) or 5th days (ps> 0.173).

![Fig. 2](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/08/25/2020.08.23.20180489/F2.medium.gif)

[Fig. 2](http://medrxiv.org/content/early/2020/08/25/2020.08.23.20180489/F2)

Fig. 2 
The effects of oxytocin (OT) on arousal and intensity ratings for face emotions.

(a) OT decreased arousal and intensity ratings of fear faces on both the 1st and the 5 days of treatment in the group receiving OT on alternate days (OT3) compared with placebo (PLC) and ratings were significantly reduced on day 5 relative to day 1. The group receiving OT every day (OT5) did not show significant effects on either day 1 or 5. No effects of OT on ratings for (b) anger, (c) happy and (d) happy face emotions. # p < 0.1, * p < 0.05, ** p < 0.01 two-tailed, Bars depict standard errors.

The whole brain voxel based morphometry analysis revealed no significant effects of acute or repeated doses of OT on gray matter volumes.

### Associations between amygdala responses and behavioral ratings to emotional faces

Associations between right amygdala activation (extracted parameter estimates from the right amygdala mask as provided by the Brainnetome atlas) and emotional arousal and intensity scores to fear faces were observed in the PLC group on the 1st day (arousal: r = 0.295, p = 0.023; intensity: r = 0.455, p = 0.001) demonstrating that greater amygdala activation in response to fear faces was associated with increased anxiety. This association was absent in both OT groups on the 1st day (arousal: OT3, r = 0.041, p = 0.398; OT5, r = 0.050, p = 0.366; intensity: OT3, r = 0.098, p = 0.265; OT5, r = 0.164, p = 0.130). A similar trend occurred on the 5th day (arousal: PLC, r = 0.306, p = 0.019; OT3, r = −0.011, p = 0.471; OT5 group, r = −0.120, p = 0.206; intensity: PLC, r = 0.212, p = 0.079; OT3, r = −0.108, p = 0.246; OT5 group, r =-0.043, p = 0.385) (**Fig.3**). There was a similar relationship with responses to angry faces for intensity but not arousal ratings in the PLC group on the 1st day. (arousal: r = 0.107, p = 0.240; intensity: r = 0.297, p = 0.023). No associations were found for responses to neutral or happy faces.

![Fig. 3](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/08/25/2020.08.23.20180489/F3.medium.gif)

[Fig. 3](http://medrxiv.org/content/early/2020/08/25/2020.08.23.20180489/F3)

Fig. 3 
Associations between intensity and arousal ratings of fear faces and amygdala activation in the placebo (PLC) and two OT treatment groups (OT3 and OT) on days 1 and 5. The magnitude of right amygdala (from the Brainnetome atlas) responses was positively correlated with arousal and intensity ratings in the PLC group on day 1 but not in the OT groups. There was a similar trend on day 5. *p< 0.05. **p<0.01

### Primary outcome measure: effects of acute and repeated oxytocin doses on resting-state functional connectivity

To facilitate a more precise localization of the coordinates of amygdala for the resting state analysis a voxel-wise mixed ANOVA with treatment group (PLC, OT3, OT5) as between-subject, time point (1st and 5th days) and face emotion (happy, neutral, angry and fear) as within-subject factor was computed on the task-based data. This analysis revealed a significant three-way interaction effect in the right amygdala after small volume correction for the bilateral mask encompassing the entire amygdala (k = 7, pFWE = 0.021, F6,405=4.47, peak MNI co-ordinate: x=18, y=-7, z=-16) (**Fig.4a**) confirming the findings from the ROI analysis which revealed a significant interaction effect for the right amygdala only. An amygdala seed-to-whole brain fMRI resting-state analysis employing this region as seed (6 radius sphere centered at x=18, y=-7, z=-16) by mixed-effect ANOVA revealed a main effect of treatment (OT3, OT5 and PLC) on amygdala functional coupling with the anterior prefrontal cortex (peak MNI coordinate x = −3, y =53, z = −19, F2, 135 = 16.55, pFWE = 0.001, k = 88) (see Fig. 4a) but no treatment group x time point interaction. Right amygdala intrinsic connectivity with the prefrontal cortex was stronger in both the OT3 (k = 74, pFWE = 0.014, x = −3, y = 53, z = −19) and OT5 groups (k = 127, pFWE = 0.001, x = −3, y = 53, z = −19) relative to PLC. The OT3 and OT5 groups did not differ significantly on the 1st (p = 0.147) or 5th (p = 0.107) days (see Fig. 4b). Correlation analysis revealed no significant association between resting state and amygdala task-related changes (r=-0.08, p=0.59 in PLC group) or behavioral ratings of emotional faces (−0.085<rs<0.055, ps > 0.57).

![Fig. 4](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/08/25/2020.08.23.20180489/F4.medium.gif)

[Fig. 4](http://medrxiv.org/content/early/2020/08/25/2020.08.23.20180489/F4)

Fig. 4 
The effects of oxytocin (OT) on resting state functional connectivity (a) schematic visualization of the identified functional pathway between the right amygdala and anterior prefrontal cortex (replicating Eckstein et al., 2017). (b) Seed to whole-brain analysis of functional connectivity (FC, seed region right amygdala, 6mm sphere, x=18, y=-7, z=-16) revealed significantly increased functional connectivity between right amygdala and anterior prefrontal cortex (pFWE =0.039, x = −36, y =59, z = −7) in both OT treatment groups (OT3 and OT5) on both, days 1 and 5 relative to placebo (PLC). Line graphs depict the extracted parameter estimates from right amygdala connectivity with anterior prefrontal cortex (M ± SE). * p < 0.05, ** p < 0.01 two-tailed. Bars indicate standard errors.

### Secondary outcome measure: Associations between OT-effects and OXTR genotype

The number of G-carriers and G-non-carriers of rs53576 and A-carriers and A-non-carriers of rs2254298 did not differ between the three groups and both SNPs satisfied Hardy Weinberg Equilibrium (see **Tables S2 and S3**). A four-way mixed ANOVA with treatment group and genotype as between- and time point and face emotion as within-subject factors revealed significant treatment x genotype interactions for both rs2254298 and rs53576 (rs 2254298: F2, 105 = 3.17, p = 0.046, rs53576: F2, 105 = 3.23, p = 0.044) for amygdala responses to emotional faces, although they would not survive a Bonferroni correction for 4 genotypes (i.e. p<0.0125). Post-hoc analysis showed that reduced amygdala reactivity to emotional faces in both OT groups occurred in A-carriers of rs2254298 (ps < 0.009) but not A-non-carriers (ps>0.587) and G-non-carriers of rs 53576 (ps < 0.002) but not G-carriers (ps>0.196) (see Fig. 5a-b). For behavioral ratings there were significant treatment group x genotype interactions for both arousal and intensity ratings for rs2254298 (arousal: F2, 105= 6.60, p = 0.002, *η*2p = 0.11; intensity: F2, 105= 3.85, p = 0.024, *η*2p = 0.07). Post-hoc analysis showed that decreased ratings for faces in the OT3 group only occurred in A-carriers (arousal: versus PLC: p = 0.002, *d* = 4.68, 95% CI, 0.363 to 1.524 and OT5 p < 0.001, *d* = 6.56, 95% CI, 0.743 to 1.904; intensity: PLC: p = 0.021, *d* = 3.40, 95% CI, 0.109 to 1.301; OT5: p = 0.003, *d* = 4.40, 95% CI, 0.314 to 1.507). No significant main or interaction effects involving treatment and genotype were found for resting state functional connectivity between amygdala and prefrontal cortex (rs 2254298 all ps > 0.123; rs53576 all ps > 0.566).

![Fig. 5](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/08/25/2020.08.23.20180489/F5.medium.gif)

[Fig. 5](http://medrxiv.org/content/early/2020/08/25/2020.08.23.20180489/F5)

Fig. 5 
Influence of oxytocin (OT) receptor genotype on the effects of intranasal OT on right amygdala and behavioral responses to face emotions (n = 111 subjects). (a) for rs 53576 only G-non carriers (i.e. AA) showed reduced amygdala responses to faces (PLC, nG+=15, nG-=21; OT3, nG+=17, nG-=24; OT5, nG+=20, nG-=14); (b) for rs2254298 only A-carriers (i.e. AA and AG) in the groups with intranasal oxytocin (OT) treatment showed reduced amygdala responses to faces (PLC, nA+=20, nA-=16; OT3, nA+=23, nA-=18; OT5, nA+=20, nA-=14); (c) for rs2254298 only A-carriers (i.e. AA and AG) in the groups with intranasal oxytocin (OT) treatment showed reduced emotional arousal ratings ** p < 0.01, two-tailed t-test (Bonferroni corrected significance threshold of p = 0.0125). Bars indicate standard errors.

## Discussion

Overall, our findings validate OT-reductions in amygdala responses to face emotions and associated emotional intensity and arousal ratings for fear faces, together with increased resting state functional connectivity between the amygdala and prefrontal cortex as robust markers for its putative therapeutic mechanisms of action. While, these task-dependent and resting-state effects of OT show a different sensitivity to repeated intranasal doses and OXTR genotype our findings indicate that to achieve maximal effects of chronic treatment on neural and behavioral responses to face emotion stimuli and resting-state changes the optimal protocol may be to administer 24IU doses every other day rather than daily.

Our findings are highly consistent with preclinical animal models demonstrating OXTR desensitization following repeated doses of OT in some brain regions (Terenzi and Ingram, 2005; Smith et al., 2005; Stoop, 2012) and that chronic administration can reduce brain OXTR expression (Du et al., 2017) and alter patterns of neurochemical release (Benner et al., 2018). Importantly, chronic doses of OT in rodents fail to produce anxiolytic effects normally seen with single doses (Du et al., 2017) which mirrors our present observations. The apparent long-lasting desensitization effects of daily OT administration on amygdala responses to threatening and happy faces may be contributed to by dose magnitude and so possibly lower daily doses might be less problematic. We also observed some evidence for reduced amygdala and behavioral responses to repeated emotional face presentations in the PLC group in line with previous observations (Fischer et al., 2003; Wright et al., 2001), even though different sets of face stimuli were used at the two time-points. This may have served to reduce time-dependent effects of OT to some extent.

The amygdala is one of the main neural substrates mediating OT’s functional effects and its attenuation of reactivity to threatening stimuli is considered as a putative primary therapeutic mechanism of action (see Jurek and Neumann, 2018; Kendrick et al., 2017; Neumann and Slattery, 2016). However, in the current study we demonstrated similar effects of OT in reducing right amygdala responses to all emotional faces, although desensitization effects were only significant for threatening (fear and anger) and happy faces. While some previous studies have found effects of OT only for fear faces (Spengler et al., 2017) or for fear and angry ones (Kirsch et al., 2005), our current findings are in agreement with a number of others reporting that it can reduce amygdala responses to both positive and negative (Domes et al., 2007) or all face emotions (Koch et al., 2016; Quintana et al., 2016). A possible explanation for these differences could be that studies, including our own, using presentation strategies where subjects are not directly required to identify or discriminate between face emotions (Domes et al., 2007; Koch et al., 2016; Quintana et al., 2016) have tended to report valence-independent results. Importantly however, only the magnitude of amygdala responses to fear faces was positively associated with arousal and intensity ratings in the PLC group and OT exerted an anxiolytic action by both reducing these ratings and abolishing their correlation with amygdala activation. Dysregulations in amygdala and behavioral responses to threatening stimuli have been observed across major psychiatric disorders, particularly those characterized by marked social impairments and anxiety (Hennessey et al., 2018; Neumann and Slattery, 2016). There was some evidence that behavioral anxiolytic effects of OT were greater after repeated compared to single doses in the group receiving OT every other day for angry faces but not in terms of amygdala responses. Possibly the short 5 day period for OT-treatment was not sufficient to demonstrate robust enhancement effects of repeated doses and time-dependent enhancement effects of longer duration treatment need to be investigated in studies on clinical populations.

A major finding from our secondary outcome measure is that OT’s attenuations of neural and behavioral responses to face emotions are highly dependent on OXTR genotype. Only G non-carriers of rs53576 and A-carriers of rs2254298 showed significant reductions in amygdala responses to face emotions under OT although only A-carriers of rs2254298 showed evidence for significant reductions in arousal and intensity ratings. A-carriers of both SNPs have frequently been associated with social dysfunction in autism and social anxiety disorders (Cataldo et al., 2018; Jurek and Neumann, 2018). A recent haplotype-based analysis of OXTR SNPs including rs53576 and rs2254298 also indicated an association with sensitivity to OT-effects on face recognition (Chen et al., 2015).

Interestingly, task-related and intrinsic network changes produced by OT showed a strikingly different sensitivity to repeated doses and OXTR genotype and there were no associations between resting state connectivity and either amygdala or behavioral responses to face emotions. However, for the resting state effects of OT there was no advantage of giving OT daily as opposed to every other day. Intrinsic networks may therefore be less sensitive to the effects of repeated OT administration or OXTR genotype than those engaged by external social emotional stimuli, although arguably the effects of OT on neural processing during actual processing of social cues should be of greatest therapeutic relevance.

There are several limitations in the current study. Firstly, we only included male subjects to increase relevance to potential therapeutic use in autism and to avoid menstrual cycle effects. There is however increasing evidence for sex-dependent effects of OT, notably for amygdala responses (Gao et al., 2016; Jurek and Neumann, 2018; Lieberz et al., 2019; Luo et al., 2017) and also for rs53576 (Feng et al., 2015). Thus, females may show different responses to repeated doses of OT and associations with receptor genotype. Secondly, for the genotype analysis, subject numbers are relatively low for establishing robust associations. While treatment x genotype interactions only passed Bonferroni correction for arousal ratings, post-hoc tests for both behavioral and amygdala effects were highly significant.. The AA allele of rs53576 also occurs more frequently in Asian and Caucasian populations (Butovskaya et al., 2016). Finally, we only examined the effects of one dose and time course for the effects of OT and there may be dose-dependent and temporal variations. The 24IU dose and 45 min post-treatment time course used are the same as in many previous studies demonstrating both task and resting-state effects of OT (Kendrick et al., 2017) and maximal effects on reduced amygdala responses to fear faces occur with a 24IU dose from 45-75 min (Spengler et al., 2017). However, higher doses may produce different effects (Spengler et al., 2017) and a recent regional cerebral blood flow study has reported reduced resting amygdala activity at shorter time points following 40IU intranasal OT (Martins et al., 2020).

In conclusion, the current study provides the first evidence for an important influence of OT dose frequency and receptor genotype on its modulation of neural and behavioral responses to emotional faces in healthy human subjects. Dose frequency therefore requires further empirical evaluation in patient populations given that it can critically determine treatment efficacy in clinical trials employing chronic administration.

## Data Availability

All data referred to in the manuscript will be available

## Funding

This project was supported by National Natural Science Foundation of China (NSFC) grant numbers 31530032 (KMK) and 91632117 (BB).

## Author contributions

JK and KMK designed the experiment. JK, YZ and CS carried out the experiment. JK, KMK, CM, FZ and BB analyzed the experiment and JK, KMK, CM and BB wrote the paper. All authors contributed to the conception of the study and approved the paper.

## Competing interests

The authors declare that they have no competing interests.

## Acknowledgement

We thank the reviewers and editor for their comments and suggestions for improving the manuscript.

*   Received August 23, 2020.
*   Revision received August 23, 2020.
*   Accepted August 25, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  Bales, K. L., Perkeybile, A. M., Conley, O. G., Lee, M. H., Guoynes, C. D., Downing, G. M., … & Mendoza, S. P. (2013). Chronic intranasal oxytocin causes long-term impairments in partner preference formation in male prairie voles. Biological Psychiatry, 74(3), 180–188.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2012.08.025&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23079235&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000321443100008&link_type=ISI) 

2.  Benner, S., Aoki, Y., Watanabe, T., Endo, N., Abe, O., Kuroda, M., … & Kasai, K. (2018). Neurochemical evidence for differential effects of acute and repeated oxytocin administration. Molecular Psychiatry, 1–11.
    
    

3.  Brüne, M. (2012). Does the oxytocin receptor polymorphism (rs2254298) confer ‘vulnerability’ for psychopathology or ‘differential susceptibility’? Insights from evolution. BMC Medicine, 10(1), 38.
    
    

4.  Butovskaya, P. R., Lazebny, O. E., Sukhodolskaya, E. M., Vasiliev, V. A., Dronova, D. A., Fedenok, J. N., … & Butovskaya, M. L. (2016). Polymorphisms of two loci at the oxytocin receptor gene in populations of Africa, Asia and South Europe. BMC Genetics, 17(1), 17.
    
    

5.  Cataldo, I., Azhari, A., & Esposito, G. (2018). A review of oxytocin and arginine-vasopressin receptors and their modulation of autism spectrum disorder. Frontiers in Molecular Neuroscience, 11, 27.
    
    

6.  Chen, F. S., Kumsta, R., Dvorak, F., Domes, G., Yim, O. S., Ebstein, R. P., & Heinrichs, M. (2015). Genetic modulation of oxytocin sensitivity: a pharmacogenetic approach. Translational Psychiatry, 5(10), e664–e664.
    
    

7.  Cremers, H. R., Wager, T. D., & Yarkoni, T. (2017). The relation between statistical power and inference in fMRI. PloS one, 12(11), e0184923.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0184923)&link_type=DOI) 

8.  Domes, G., Heinrichs, M., Gläscher, J., Büchel, C., Braus, D. F., & Herpertz, S. C. (2007). Oxytocin attenuates amygdala responses to emotional faces regardless of valence. Biological Psychiatry, 62(10), 1187–1190.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2007.03.025&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17617382&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000250905800019&link_type=ISI) 

9.  Du, P., He, Z., Cai, Z., Hao, X., Dong, N., Yuan, W., … & Tai, F. (2017). Chronic central oxytocin infusion impairs sociability in mandarin voles. Pharmacology Biochemistry and Behavior, 161, 38–46.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pbb.2017.09.006&link_type=DOI) 

10. Eckstein, M., Becker, B., Scheele, D., Scholz, C., Preckel, K., Schlaepfer, T. E., … & Hurlemann, R. (2015). Oxytocin facilitates the extinction of conditioned fear in humans. Biological Psychiatry, 78(3), 194–202.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2014.10.015&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25542304&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 

11. Eckstein, M., Markett, S., Kendrick, K. M., Ditzen, B., Liu, F., Hurlemann, R., & Becker, B. (2017). Oxytocin differentially alters resting state functional connectivity between amygdala subregions and emotional control networks: Inverse correlation with depressive traits. Neuroimage, 149, 458–467.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2017.01.078&link_type=DOI) 

12. Eklund, A., Nichols, T. E., & Knutsson, H. (2016). Cluster failure: Why fMRI inferences for spatial extent have inflated false-positive rates. Proceedings of the National Academy of Sciences, 113(28), 7900–7905.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTEzLzI4Lzc5MDAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOC8yNS8yMDIwLjA4LjIzLjIwMTgwNDg5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

13. Fan, L., Li, H., Zhuo, J., Zhang, Y., Wang, J., Chen, L., … & Fox, P. T. (2016). The human brainnetome atlas: a new brain atlas based on connectional architecture. Cerebral Cortex, 26(8), 3508–3526.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cercor/bhw157&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27230218&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 

14. Feng, C., Lori, A., Waldman, I. D., Binder, E. B., Haroon, E., & Rilling, J. K. (2015). A common oxytocin receptor gene (OXTR) polymorphism modulates intranasal oxytocin effects on the neural response to social cooperation in humans. Genes, Brain and Behavior, 14(7), 516–525.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/gbb.12234&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26178189&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 

15. Fischer, H., Wright, C.I., Whalen, P.J., McInerney, S.C., Shine, L.M., & Rauch, S.L. (2003). Brain habituation during repeated exposure to fearful and neutral faces. Brain Research Bulletin 59(5), 387–392.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0361-9230(02)00940-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12507690&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000180428200008&link_type=ISI) 

16. Friston, K. J., Holmes, A. P., Worsley, K. J., Poline, J. P., Frith, C. D., & Frackowiak, R. S. (1994). Statistical parametric maps in functional imaging: a general linear approach. Human Brain Mapping, 2(4), 189–210.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hbm.460020402&link_type=DOI) 

17. Gao, S., Becker, B., Luo, L., Geng, Y., Zhao, W., Yin, Y., … & Kendrick, K.M. (2016). Oxytocin, the peptide that bonds the sexes also divides them. Proceedings of the National Academy of Sciences, 113(27), 7650–7654.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTEzLzI3Lzc2NTAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOC8yNS8yMDIwLjA4LjIzLjIwMTgwNDg5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

18. Guastella, A. J., Gray, K. M., Rinehart, N. J., Alvares, G. A., Tonge, B. J., Hickie, I. B., … & Einfeld, S. L. (2015). The effects of a course of intranasal oxytocin on social behaviors in youth diagnosed with autism spectrum disorders: a randomized controlled trial. Journal of Child Psychology and Psychiatry, 56(4), 444–452.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/jcpp.12305&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25087908&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 

19. Guastella, A. J., Hickie, I. B., McGuinness, M. M., Otis, M., Woods, E. A., Disinger, H. M., … & Banati, R. B. (2013). Recommendations for the standardisation of oxytocin nasal administration and guidelines for its reporting in human research. Psychoneuroendocrinology, 38(5), 612–625.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.psyneuen.2012.11.019&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23265311&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000318465500002&link_type=ISI) 

20. Guidi, J., Brakemeier, E. L., Bockting, C. L., Cosci, F., Cuijpers, P., Jarrett, R. B., … & Rief, W. (2018). Methodological recommendations for trials of psychological interventions. Psychotherapy and Psychosomatics, 87(5), 276–284.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000490574&link_type=DOI) 

21. Halverson, T., Jarskog, L. F., Pedersen, C., & Penn, D. (2019). Effects of oxytocin on empathy, introspective accuracy, and social symptoms in schizophrenia: A 12-week twice-daily randomized controlled trial. Schizophrenia Research, 204, 178–182.
    
    

22. Hennessey, T., Andari, E., & Rainnie, D. G. (2018). RDoC-based categorization of amygdala functions and its implications in autism. Neuroscience & Biobehavioral Reviews, 90, 115–129.
    
    

23. Huang, H., Michetti, C., Busnelli, M., Manago, F., Sannino, S., Scheggia, D., … & Scattoni, M. L. (2014). Chronic and acute intranasal oxytocin produce divergent social effects in mice. Neuropsychopharmacology, 39(5), 1102–1114.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/npp.2013.310&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24190025&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000332918100007&link_type=ISI) 

24. Jarskog, L. F., Pedersen, C. A., Johnson, J. L., Hamer, R. M., Rau, S. W., Elliott, T., & Penn, D. L. (2017). A 12-week randomized controlled trial of twice-daily intranasal oxytocin for social cognitive deficits in people with schizophrenia. Schizophrenia Research, 185, 88–95.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.schres.2017.01.008&link_type=DOI) 

25. Jurek, B., & Neumann, I. D. (2018). The oxytocin receptor: from intracellular signaling to behavior. Physiological Reviews, 98(3), 1805–1908.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/physrev.00031.2017&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29897293&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 

26. Kendrick, K. M., Guastella, A. J., & Becker, B. (2017). Overview of human oxytocin research. In Behavioral Pharmacology of Neuropeptides: Oxytocin (pp. 321-348). Springer, Cham. pp. 321-348.
    
    

27. Kirsch, P., Esslinger, C., Chen, Q., Mier, D., Lis, S., Siddhanti, S., … & Meyer-Lindenberg, A. (2005). Oxytocin modulates neural circuitry for social cognition and fear in humans. Journal of Neuroscience, 25(49), 11489–11493.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjExOiIyNS80OS8xMTQ4OSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA4LzI1LzIwMjAuMDguMjMuMjAxODA0ODkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

28. Koch, S.B.J., van Zuiden, M., Nawijn, L., Frijling, J.L., Veltman, D.J. (2016). Intranasal oxytocin administration dampens amygdala reactivity towards emotional faces in male and female PTSD patients. Neuropsychopharmacology, 41, 1495–1504.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/npp.2015.299&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26404844&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 

29. Lieberz, J., Scheele, D., Spengler, F. B., Matheisen, T., Schneider, L., Stoffel-Wagner, B., … & Hurlemann, R. (2019). Kinetics of oxytocin effects on amygdala and striatal reactivity vary between women and men. Neuropsychopharmacology, 1-7.
    
    

30. Lee, M. R., Scheidweiler, K. B., Diao, X. X., Akhlaghi, F., Cummins, A., Huestis, M. A., … & Averbeck, B. B. (2018). Oxytocin by intranasal and intravenous routes reaches the cerebrospinal fluid in rhesus macaques: determination using a novel oxytocin assay. Molecular Psychiatry, 23(1), 115–122.
    
    

31. Lohse, M. J., & Hofmann, K. P. (2015). Spatial and Temporal Aspects of Signaling by G-Protein-Coupled Receptors. Molecular Pharmacology, 88(3), 572–578.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoibW9scGhhcm0iO3M6NToicmVzaWQiO3M6ODoiODgvMy81NzIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOC8yNS8yMDIwLjA4LjIzLjIwMTgwNDg5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

32. Luo, L., Becker, B., Geng, Y., Zhao, Z., Gao, S., Zhao, W., … & Hu, J. (2017). Sex-dependent neural effect of oxytocin during subliminal processing of negative emotion faces. Neuroimage, 162, 127–137.
    
    

33. Martins, D.A., Mazibuko, N, Zelaya, F., Vasilakopoulou, S., Loveridge, J., Oates, U., … & Paloyelis, Y. (2020). Effects of route of administration on oxytocin-induced changes in regional cerebral blood flow in humans. Nature Communications doi: 10.1038/s41467-020-14845-5
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-020-14845-5&link_type=DOI) 

34. Meyer-Lindenberg, A., Domes, G., Kirsch, P., & Heinrichs, M. (2011). Oxytocin and vasopressin in the human brain: social neuropeptides for translational medicine. Nature Reviews Neuroscience, 12(9), 524–538.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrn3044&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21852800&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 

35. Montag, C., Sindermann, C., Melchers, M., Jung, S., Luo, R., Becker, B., … & Kendrick, K. M. (2017). A functional polymorphism of the OXTR gene is associated with autistic traits in Caucasian and Asian populations. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics, 174(8), 808–816.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ajmg.b.32596&link_type=DOI) 

36. Mueller, K., Lepsien, J., Möller, H. E., & Lohmann, G. (2017). Commentary: Cluster failure: Why fMRI inferences for spatial extent have inflated false-positive rates. Frontiers in Human Neuroscience, 11, 345.
    
    

37. Neumann, I. D., & Slattery, D. A. (2016). Oxytocin in general anxiety and social fear: a translational approach. Biological Psychiatry, 79(3), 213–221.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2015.06.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26208744&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 

38. Paloyelis, Y., Doyle, O. M., Zelaya, F. O., Maltezos, S., Williams, S. C., Fotopoulou, A., & Howard, M. A. (2016). A spatiotemporal profile of in vivo cerebral blood flow changes following intranasal oxytocin in humans. Biological Psychiatry, 79(8), 693–705.
    
    

39. Parker, K. J., Garner, J. P., Libove, R. A., Hyde, S. A., Hornbeak, K. B., Carson, D. S., … & Hardan, A. Y. (2014). Plasma oxytocin concentrations and OXTR polymorphisms predict social impairments in children with and without autism spectrum disorder. Proceedings of the National Academy of Sciences, 111(33), 12258–12263.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTExLzMzLzEyMjU4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMjUvMjAyMC4wOC4yMy4yMDE4MDQ4OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

40. Pierce, K. L., Premont, R. T., & Lefkowitz, R. J. (2002). Seven-transmembrane receptors. Nature Reviews Molecular Cell Biology, 3(9), 639–650.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrm908&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12209124&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000177787300013&link_type=ISI) 

41. Poldrack, R. A., Baker, C. I., Durnez, J., Gorgolewski, K. J., Matthews, P. M., Munafò, M. R., … & Yarkoni, T. (2017). Scanning the horizon: towards transparent and reproducible neuroimaging research. Nature Reviews Neuroscience, 18(2), 115.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrn.2016.167&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28053326&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 

42. Quintana, D. S., Rokicki, J., van der Meer, D., Alnæs, D., Kaufmann, T., Córdova-Palomera, A., … & Westlye, L. T. (2019). Oxytocin pathway gene networks in the human brain. Nature Communications, 10(1), 1–12.
    
    

43. Quintana, D.S., Westlye, L.T., Alnaes, D., Rustan, Ø.G., Kaufmann, T.,…& Andreassen, O.A. (2016). Low dose intranasal oxytocin delivered with Breath Powered device dampens anygdala response to emotional stimuli: A peripheral effect-controlled within-subjects randomized dose-response fMRI trial. Psychoneuroendocrinology, 69, 180–188.
    
    

44. Russell, J. A. (2018). 50 years of advances in neuroendocrinology. Brain and Neuroscience Advances, 2, 1–20.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/2398212818812012&link_type=DOI) 

45. Smith, M. P., Ayad, V. J., Mundell, S. J., McArdle, C. A., Kelly, E., & López Bernal, A. (2006). Internalization and desensitization of the oxytocin receptor is inhibited by Dynamin and clathrin mutants in human embryonic kidney 293 cells. Molecular Endocrinology, 20(2), 379–388.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/me.2005-0031&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16179383&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000235129600012&link_type=ISI) 

46. Spengler, F. B., Schultz, J., Scheele, D., Essel, M., Maier, W., Heinrichs, M., & Hurlemann, R. (2017). Kinetics and dose dependency of intranasal oxytocin effects on amygdala reactivity. Biological Psychiatry, 82(12), 885–894.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2017.04.015&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28629540&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 

47. Sripada, C. S., Phan, K. L., Labuschagne, I., Welsh, R., Nathan, P. J., & Wood, A. G. (2013). Oxytocin enhances resting-state connectivity between amygdala and medial frontal cortex. International Journal of Neuropsychopharmacology, 16(2), 255–260.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S1461145712000533&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22647521&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 

48. Stoop, R. (2012). Neuromodulation by oxytocin and vasopressin. Neuron, 76(1), 142–159.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuron.2012.09.025&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23040812&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000309799000011&link_type=ISI) 

49. Striepens, N., Kendrick, K. M., Hanking, V., Landgraf, R., Wüllner, U., Maier, W., & Hurlemann, R. (2013). Elevated cerebrospinal fluid and blood concentrations of oxytocin following its intranasal administration in humans. Scientific Reports, 3, 3440.
    
    

50. Terenzi, M. G., & Ingram, C. D. (2005). Oxytocin-induced excitation of neurones in the rat central and medial amygdaloid nuclei. Neuroscience, 134(1), 345–354.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroscience.2005.04.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15961240&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 

51. Wang, D., Yan, X., Li, M., & Ma, Y. (2017). Neural substrates underlying the effects of oxytocin: a quantitative meta-analysis of pharmaco-imaging studies. Social Cognitive and Affective Neuroscience, 12(10), 1565–1573.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/scan/nsx085&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29048602&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 

52. Watanabe, T., Kuroda, M., Kuwabara, H., Aoki, Y., Iwashiro, N., Tatsunobu, N., … & Kasai, K. (2015). Clinical and neural effects of six-week administration of oxytocin on core symptoms of autism. Brain, 138(11), 3400–3412.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awv249&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26336909&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 

53. Wright, C. I., Fischer, H., Whalen, P. J., McInerney, S. C., Shin, L. M., Rauch, S. L. (2001). Differential prefrontal cortex and amygdala habituation to repeatedly presented emotional stimuli. Neuroreport, 12,379–383.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00001756-200102120-00039&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11209954&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F25%2F2020.08.23.20180489.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000166702100039&link_type=ISI) 

54. Yamasue, H., Okada, T., Munesue, T., Kuroda, M., Fujioka, T., Uno, Y., … & Yoshimura, Y. (2018). Effect of intranasal oxytocin on the core social symptoms of autism spectrum disorder: a randomized clinical trial. Molecular Psychiatry, 1-10.
    
    

55. Chao-Gan, Y., & Yu-Feng, Z. (2010). DPARSF: a MATLAB toolbox for “pipeline” data analysis of resting-state fMRI. Frontiers in Systems Neuroscience, 4.
    
    

56. Yang, S., Dong, X., Guo, X., Han, Y., Song, H., Gao, L., … & Zhang, X. (2017). Serum oxytocin levels and an oxytocin receptor gene polymorphism (rs2254298) indicate social deficits in children and adolescents with autism spectrum disorders. Frontiers in Neuroscience, 11, 221.
    
    

57. Young, L. J., & Barrett, C. E. (2015). Can oxytocin treat autism? Science, 347(6224), 825–826.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNDcvNjIyNC84MjUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOC8yNS8yMDIwLjA4LjIzLjIwMTgwNDg5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

58. Zhang, Z., & Yuan, K.H. (2018). Practical Statistical Power Analysis Using Webpower and R (Eds). Granger, IN: ISDSA Press.
    
    

59. Zhao, W., Becker, B., Yao, S., Ma, X., Kou, J., & Kendrick, K. M. (2019). Oxytocin enhancement of the placebo effect may be a novel therapy for working memory impairments. Psychotherapy and Psychosomatics, 88(2), 125–126.
    
    

60. Zhao, Z., Yao, S., Li, K., Sindermann, C., Zhou, F., Zhao, W., … & Becker, B. (2019). Real-time functional connectivity-informed neurofeedback of amygdala-frontal pathways reduces anxiety. Psychotherapy and Psychosomatics, 88(1), 5–15.