Endogenous oxytocin levels in children with autism: Associations with cortisol levels and oxytocin receptor gene methylation
============================================================================================================================

* Evenepoel Margaux
* Moerkerke Matthijs
* Daniels Nicky
* Chubar Viktoria
* Claes Stephan
* Turner Jonathan
* Vanaudenaerde Bart
* Willems Lynn
* Verhaeghe Johan
* Prinsen Jellina
* Steyaert Jean
* Boets Bart
* Alaerts Kaat

## Abstract

**Background** Alterations in the brain’s oxytocinergic system have been suggested to play an important role in the pathophysiology of autism spectrum disorder (ASD), but insights from pediatric populations are sparse.

**Methods** We examined salivary oxytocin in school-aged children with (n=80) and without (n=40) ASD (boys/girls 4/1), as well as characterizations of DNA methylation (DNAm) of the oxytocin receptor gene (*OXTR*). Cortisol levels were also assessed to examine links between the oxytocinergic system and hypothalamic-pituitary-adrenal (HPA) axis signaling.

**Results** Children with ASD displayed altered (diminished) oxytocin levels in the morning, but not in the afternoon, after a mildly stress-inducing social interaction session. Notably, in the control group, higher oxytocin levels were predictive of lower stress-induced cortisol, likely reflective of a *protective* stress-regulatory mechanism for buffering HPA stress activity. In children with ASD, on the other hand, a more *reactive* stress regulatory mechanism was evident, involving a significant rise in oxytocin levels from the morning to the afternoon upon stress-induced cortisol release, i.e., to reactively cope with heightened HPA activity.

Regarding epigenetic modifications, no overall pattern of *OXTR* hypo- or hypermethylation was evident in ASD. In control children, a notable association between *OXTR* methylation and levels of cortisol was evident, likely indicative of a compensatory downregulation of *OXTR* methylation (higher oxytocin receptor expression) in children with heightened HPA axis activity.

**Conclusion** Together, these observations bear important insights into altered oxytocinergic signaling in ASD, which may aid in establishing relevant biomarkers for diagnostic and/or treatment evaluation purposes targeting the oxytocinergic system in ASD.

## Introduction

Autism spectrum disorder (ASD) is a neurodevelopmental condition characterized by difficulties with social communication and interaction, combined with expressions of restricted and repetitive behaviors and interests. Aside these core ASD symptoms, individuals with ASD often also display a broad range of comorbid symptoms including attachment difficulties 1 and (social) stress and anxiety 2.

While the etiology and underlying neurobiological mechanisms responsible for the heterogeneous clinical presentation of ASD are not yet fully understood, recent lines of evidence suggest that alterations in the brain’s oxytocinergic system play an important role in the pathophysiology of ASD 3. Oxytocin is a nonapeptide produced in the paraventricular and supraoptic nuclei of the hypothalamus 4, from which oxytocinergic neurons project to various regions of the brain’s reward system and social brain circuits (e.g. amygdala) 5. Through its widespread neuromodulatory effects, oxytocin is identified as a key mediator of affiliative and prosocial behaviors, including interpersonal bonding, attachment and social attunement 6–8.

Initial preclinical studies have demonstrated links between aberrant oxytocinergic signaling and ASD symptom expression in rodent models of autism 9–11. However, investigating oxytocin concentrations in biological samples (mostly plasma and saliva) collected from humans with ASD, yielded a more mixed pattern of results (see meta-analysis of Moerkerke et al. (2021) 12). In particular, only young children with ASD (6-to-9 year old) showed reduced oxytocin levels compared to controls in four out of five studies 13–17, whereas none of the studies including 9-to-12 year old children with ASD reported a significant group difference in endogenous oxytocin levels 18–20. Moreover, studies including adolescents (11-to-16 year old) 18–20 and adults 21,22 with ASD revealed a mixed pattern of results, demonstrating no clear evidence for alterations in oxytocin levels in ASD.

Aside assessments of circulating oxytocin, also variations in epigenetic modifications of the oxytocin receptor gene (*OXTR*) have been examined in relation to ASD. The functioning of the circulating oxytocin is dependent on the functioning of the oxytocin receptor 23, which itself is encoded by the *OXTR* 24. DNA methylation (DNAm) is one of the most extensively studied epigenetic mechanisms involved in the regulation of gene transcription, such that increased DNAm frequency is generally associated with decreased gene transcription, and therefore less expression of the respective protein 25. While evidence is sparse, a recent review on *OXTR* DNAm by Moerkerke et al. (2021) 26 pointed towards a pattern of *hyper*methylation in adults with ASD 27,28, indicative of reduced *OXTR* expression. Contradictory, also patterns of *hypo*methylation were identified in children with ASD 29,30, suggesting differential ASD-related alterations in *OXTR* DNAm dependent on developmental stage. However, prior studies not only varied largely in terms of population-dependent characteristics (such as age), but also with regard to the investigated DNAm *OXTR* sites (e.g. intron versus exon CpG sites) as well as other design-related factors (e.g. adopted questionnaires/scales for assessing social/ASD-related traits) 26.

Despite these initial insights into the role of oxytocin hormonal variations and *OXTR* epigenetic modifications in ASD, research directly investigating the link between *OXTR* DNAm and hormonal oxytocin levels remains sparse. One study by Dadds et al. (2014) demonstrated that increased *OXTR* DNAm in peripheral blood was associated with lower plasma oxytocin levels in children/adolescents with oppositional-defiant/conduct disorder (4-to-16 year old) 31. However, another study in adults with psychotic disorders found no significant association between methylation of site −934 and plasma oxytocin 32. Likewise, a later study with female macaques also failed to identify any significant relationship between *OXTR* DNAm of distinct CpG sites and oxytocin levels in cerebrospinal fluid 33.

Finally, several lines of evidence also suggest a significant interplay between the oxytocinergic system and signaling of the hypothalamus-pituitary-adrenal (HPA) axis in stress regulation 21. For example, prior research in animals and humans provided indications of coupled releases of oxytocin and the HPA ‘stress-hormone’ cortisol, particularly in response to stressors, allowing oxytocin to inhibit stress-induced rises in cortisol 34. More research is needed however, to investigate these associations further, particularly in neuropsychiatric conditions such as ASD, in which elevated cortisol levels (see 35 for recent meta-analysis) and altered oxytocinergic signaling are apparent.

To fill these gaps, the present study aimed to examine possible alterations in peripheral hormonal levels of oxytocin and cortisol (assessed from saliva samples) as well as possible alterations in *OXTR* DNAm in a representative sample of school-aged children with and without ASD. Importantly, while endocrinological levels of e.g. cortisol are known to evolve during the course of the day 36 and in response to social encounters or stressful situations 32,37, this possibility of diurnal variations has been largely overlooked in prior research examining hormonal variations in oxytocin in ASD populations. It is however well-established that also oxytocin levels can covary depending on the impact of experimental (social) stimulations 38. The present study will therefore assess, for the first time, the interplay between oxytocin and cortisol hormonal levels in children with ASD, at two time points, either at ‘baseline’ in the morning or in the afternoon, after a (mildly stressful) social interaction session.

## Methods

This study was approved by the Ethical Committee for Biomedical Research at the University of Leuven, KU Leuven (S61358) in accordance with The Code of Ethics of the World Medical Association (Declaration of Helsinki). Written informed consent from the parents and assent from the child were obtained prior to the study.

### 1. Participants

Eighty children with a formal diagnosis of ASD, aged between 8-12 years, were recruited through the Leuven Autism Expertise Centre at the Leuven University Hospital between July 2019 and January 2021. Also, 40 age- and gender-matched typically developing (TD) peers were recruited.

Main inclusion criteria comprised a clinical diagnosis of ASD (only for children with ASD), premenstrual girls, intelligence quotient above 70 and native Dutch speaker. Main exclusion criteria comprised a history of any neurological disorder (stroke, concussion, epilepsy etc.), any significant physical disorder (liver, renal, cardiac pathology) or any neuropsychiatric diagnosis (only for TD children).

In addition, the Autism Diagnostic Observation Schedule (ADOS-2) 39 was acquired (**Table 1**). For all children, estimates of intelligence were acquired using four subtests of the Wechsler Intelligence Scale for Children, Fifth Edition, Dutch version 40 (**Table 1**).

View this table:
[Table 1.](http://medrxiv.org/content/early/2022/12/16/2022.12.15.22283521/T1)

Table 1. Participants’ characteristics

Children were also thoroughly characterized on autism symptom domains, using the parent-reported versions of the *Social Responsiveness Scale*, second edition (SRS-2) 41,42 and the *Repetitive Behavior Scale-Revised* (RBS-R) 43. Also, the parent-rated *Child Behavior Checklist* (CBCL) was obtained to assess behavioral problems and competences 44 (see **Supplementary Methods** for detailed descriptions of the adopted questionnaires).

As outlined in **Table 1**, groups were matched on age and sex, although verbal and performance intelligence quotient were overall higher in the TD group. As expected, children of the ASD group demonstrated significantly higher scores on the parent-rated SRS-2 and RBS-R, indicating more social impairments and more frequent expressions of restricted and repetitive behavior (*p* < .001, **Table 1**). The parent-rated CBCL also showed diagnosis-related effects, indicating more severe deficits in the ASD group, compared to the TD group (*p* < .05, **Table 1**). Note that the (biological) data collected for the current report were part of a larger clinical study including additional neurophysiological assessments (see also next section).

### 2. Oxytocin and cortisol salivary concentrations

For each child, salivary samples were acquired at two time points: (i) an AM sample, acquired at home, in the morning, within 30 minutes after awakening and before breakfast; and (ii) a PM sample, acquired in the afternoon, after completion of an experimental session performed at the Leuven University hospital. Importantly, the PM sample was collected within 30 minutes after finalizing an experimental test session, consisting of a semi-structured social interaction with an unknown experimenter, which could be experienced as moderately stressful, especially to the pediatric participants (see **Supplementary methods)**.

Salivary samples were collected using Salivette cotton swabs (Sarstedt AG & Co.). For the analysis of salivary **oxytocin levels**, the commercial enzyme immunoassay oxytocin ELISA kit of Enzo Life Sciences, Inc. was used, similar to prior studies 45–47. All sample extraction and concentration procedures were conducted in accordance with manufacturer’s instructions. Measurements were performed on undiluted samples (100 μl), and sample concentrations were calculated according to plate-specific standard curves.

Analysis of the **cortisol levels** was performed using the commercial enzyme immunoassay Cortisol ELISA kit of Salimetrics, Europe. Measurements were performed on undiluted samples (25 μl), and sample concentrations were calculated according to plate-specific standard curves. More detailed information regarding the salivary collection procedures and analyses is provided in **Supplementary methods**.

### 3. DNA methylation of the oxytocin receptor gene

Additional salivary samples were obtained via the Oragene DNA sample collection kit (DNA Genotek Inc., Canada) to address epigenetic variations in the level of methylation (DNAm) at two CpG sites (−934 and -924) of *OXTR* (hg19, chr3:8,810,729-8,810,845).

After data collection, 200 ng DNA was extracted from the samples and bisulfite converted following the manufacturer’s protocol (EZ-96 DNAm Kit, Zymo Research, Irvine, CA, USA). Bisulfite converted DNA was stored at −80 °C until further analysis. Next, the levels of methylation at two CpG sites (i.e., -934 and -924) of *OXTR* (hg19, chr3:8,810,729-8,810,845) were determined using Pyrosequencer (Qiagen, Hilden, Germany) and analyzed using Pyromark Q96 software. Laboratory procedures and analyses were conducted in accordance with manufacturer’s protocols and software settings (e.g. for determining unreliable samples) 44. Protocols for the PCR amplification and Pyrosequencing analysis were adapted from Krol et al. (2019) 48. More information regarding the adopted PCR primers can be found in **Supplementary methods**.

### 4. Data handling and statistical procedures

Prior to analysis, hormonal and DNAm data were log-transformed (log10) to deal with skewed data. Next, diagnosis-related differences were assessed by subjecting oxytocin and cortisol hormonal data to repeated-measure analyses-of-variance (ANOVA) with the between-subject factor group (ASD, TD) and the within-subject factor time point (AM, PM) (see **Supplementary table 1-2** for full ANOVA models).

*OXTR* DNAm data were subjected to independent sample t-tests to examine diagnosis-related differences in DNAm at CpG site -924 and -934.

Pearson correlation analyses were performed to assess associations between the two hormonal systems in children with or without ASD and between hormonal levels and *OXTR* DNAm. All statistical analyses were executed with SPSS (version 28.0, IBM).

## Results

### Diagnosis-related differences

#### Oxytocin levels

A repeated-measure ANOVA revealed no main effects of group (*F*(1,113) = 2.45; *p* = .120) or time-point (*F*(1,113) = 0.87; *p* = .354). A significant group x time point interaction (*F*(1,113) = 4.23; *p* = .042) was noted, indicating significantly lower morning oxytocin levels in ASD, compared to TD children (AM sample, independent-samples *t*(115) = - 2.68; *p* = .009, **Figure 1A and Table 1**) but no differential afternoon oxytocin levels in ASD, compared to TD children (PM sample *t*(115) = -0.49; *p* = .629, **Figure 1A and Table 1**). Notably, while oxytocin levels remained overall constant from AM to PM in the TD group (dependent-samples *t*(38) = -0.67; *p* = .507), the ASD group displayed a steep and significant increase in oxytocin levels from AM to PM (*t*(75) = 2.61; *p* =.011) (**Figure 1A**).

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/12/16/2022.12.15.22283521/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2022/12/16/2022.12.15.22283521/F1)

Figure 1. Diagnosis-related differences
Salivary oxytocin (**panel A**) and cortisol hormonal levels (**panel B**), as well as OXTR gene methylation data (**panel C**) are visualized separately for children with ASD and typically developing children (TD).

Hormonal concentrations are visualized separately for the morning sample (AM), and the afternoon sample (PM), acquired after a mildly stress-inducing social interaction session. Methylation frequency of the OXTR gene is visualized separately for CpG site -934 and CpG site -924.

**Panel A**. Oxytocin levels were significantly higher in the TD group, compared to the ASD group, for the morning AM sample, not for the PM sample collected after the social interaction session. While oxytocin levels remained overall constant across AM and PM measurements in the TD group, the ASD group displayed a steep increase in oxytocin levels from the AM measurement to the PM measurement.

**Panel B**. Cortisol levels were not significantly different between the ASD and the TD group. Both groups displayed a similar cortisol peak in the morning and a similar level of cortisol at the PM measurement, after finalizing the social interaction session.

**Panel C**. *OXTR* DNAm at CpG sites -934 and -924 were not significantly different between the ASD group and the TD group.

#### Cortisol levels

A repeated-measure ANOVA revealed no significant main effect of group (*F*(1,116) = .002; *p* = .965) nor group x time point interaction (*F*(1,116) = 1.86; *p* = .175) (**Figure 1B and Table 1**). A significant main effect of time point was evident however (*F*(1,116) = 286.83; *p* < .001), indicating that in both groups, a steep decline in cortisol levels was evident from AM to PM. Accordingly, it appears that - as a group - children with or without ASD display a similar cortisol peak in the morning and a similar pattern of cortisol diminishment during the course of the day.

### *OXTR* DNAm

No significant diagnosis-related differences in *OXTR* DNAm were evident between the ASD and TD group, either at CpG site -934 (*t*(113) = -0.33; *p* = .745) or at CpG site -924 (*t*(113) = 0.44; *p* = .659) (**Figure 1C**).

### Association between oxytocin and cortisol hormonal levels

In the TD group, a tight interaction between oxytocin and cortisol hormonal systems was evident, indicating that high morning oxytocin levels (AM sample) were predictive of lower ‘stress-induced’ cortisol levels (PM sample) (Pearson *r* = -.35; *p* = .027) (**Figure 2A**). Higher oxytocin levels in the afternoon, i.e. after the (mildly stressful) social interaction session, were also significantly associated with reduced cortisol levels in the afternoon (PM sample) (*r* = - .34; *p* = .032) (**Figure 2B**).

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/12/16/2022.12.15.22283521/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2022/12/16/2022.12.15.22283521/F2)

Figure 2. Association between oxytocin and cortisol hormonal levels
Relationships between oxytocin and cortisol hormonal levels are visualized separately for children with ASD and typically developing children (TD).

**Panel A**. In the TD group, oxytocin levels at AM were significantly associated with cortisol levels at PM, indicating that high morning oxytocin was predictive of low ‘stress-induced’ PM cortisol levels. Relationships were only significantly evident in the TD group, not in the ASD group.

**Panel B**. Oxytocin levels at PM were also significantly associated with cortisol levels at PM, indicating that higher afternoon oxytocin levels were associated with lower ‘stress-induced’ cortisol levels. Also here, the relationship was only evident in the TD group, not in the ASD group.

**Panel C**. Finally, in the ASD group, the observed marked increases in oxytocin levels during the day (from the AM to the PM assessment) were significantly associated with higher afternoon ‘stress-induced’ cortisol levels.

Strikingly, these associations were only evident in the TD group, not in the ASD group, suggesting a lack of coupling of these hormonal systems in children with ASD (all, *r* < .13; *p* > .05). **Supplementary table 3** reports a full overview of the assessed associations between oxytocin and cortisol levels, separately for each group, including a direct comparison of correlation values between ASD and control groups.

Notably, considering the marked AM to PM increase in oxytocin levels seen in the ASD group, we additionally explored whether these AM to PM changes in oxytocin levels were potentially related to cortisol levels. It was revealed that indeed, only in the ASD group (*r* = .23; *p* = .041), not in the TD group (*r* = .01; *p* = .964), higher cortisol levels at PM were associated with a significant rise in oxytocin levels during the day (from AM to PM) (**Figure 2C**).

### Associations between hormonal levels and OXTR DNAm

#### CpG site -924

Only in the TD group (*r* = -.41; *p* = .011), not in the ASD group (*r* = .05; *p* = .682), a significant association between *OXTR* DNAm of site -924 and PM cortisol levels was evident, indicating that children with lower methylation frequency of CpG site -924 (higher *OXTR* expression) displayed higher afternoon cortisol levels (**Figure 3**). *OXTR* DNAm of site -924 was not significantly associated with AM cortisol levels or AM/PM oxytocin levels either in children with or without ASD (all, *r* < .07; *p* > .05; see **Supplementary table 3** for a full overview).

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/12/16/2022.12.15.22283521/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2022/12/16/2022.12.15.22283521/F3)

Figure 3. Associations between hormonal levels and OXTR DNAm
Associations between hormonal levels and methylation of the *OXTR*, visualized separately for children with and without ASD.

In the TD group, not in the ASD group, a significant negative association was evident between *OXTR* methylation of CpG site -924 and cortisol levels, indicating that children with lower methylation frequency (higher *OXTR* expression) displayed higher afternoon cortisol levels.

#### CpG site -934

Methylation of CpG site -934 was not significantly associated with AM or PM oxytocin or cortisol levels, either in children with or without ASD (all, *r* < .08; *p* > .05, see **Supplementary table 3**).

## Discussion

In the current study, we investigated alterations in oxytocin and cortisol hormonal systems as well as *OXTR* DNAm in school-aged children with ASD, compared to a group of age and gender matched TD children.

Only in terms of ‘baseline’ oxytocin, children with ASD displayed reduced levels of peripheral oxytocin assessed in the morning, not in terms of afternoon oxytocin levels, assessed after induction of a mildly stressful social interaction session. These findings therefore provide indications that only *‘trait-dependent’* oxytocin levels, but not stress-reactive *‘state-dependent’* levels of oxytocin are altered in ASD. Further, while overall cortisol levels were similar in children with and without ASD, differential couplings between the oxytocin and cortisol hormonal system were evident, indicating that only in the TD group, not in the ASD group, higher oxytocin levels were predictive of dampened afternoon ‘stress-induced’ cortisol levels. Finally, DNAm at two distinct CpG sites of *OXTR* were not altered in children with ASD, although note that only in the TD group, methylation frequency was significantly associated with hormonal reactivity of the cortisol system.

Diagnosis-related differences in oxytocin hormonal levels were evident in children with ASD, albeit exclusively in terms of ‘trait-dependent’ morning oxytocin, not in terms of ‘state-dependent’ afternoon oxytocin levels, measured after an experimental social interaction session. Prior meta-analytic analyses examining altered oxytocin levels in ASD yielded a mixed pattern of results, predominantly depending on the age of the included participants 12. While it is difficult to directly compare studies, the current data provide important indications that also the time of the day as well as the particular experimental and social context (with or without experimental manipulation) are important factors to consider for explaining possible variations in oxytocin levels.

Closer inspection of the recent meta-analysis of Moerkerke et al. (2021) 12 showed that only in 11 out of 18 studies the timepoint or time range of sampling was reported. Of these studies, six reported to have acquired ‘morning’ oxytocin samples (AM) 13,15,17,20,21,49, and four of these six reported similar group-related differences as in our study, indicating lower morning oxytocin levels in ASD compared to the TD group 13,15,21,49. In accordance with our study, the only two studies who indicated to have sampled oxytocin in the afternoon (PM), reported no significant group differences 18,19. Further, in the remaining studies who reported sampling timings, it appeared that acquisition timing was not standardized across individuals (i.e., ranging from morning to afternoon), therefore rendering direct comparisons with our current results was difficult 22,50,51. As the time point and context of sampling may have an important impact on observed levels of oxytocin and correlations with other outcome measures, future studies are urged to report this important design-specific information in more detail.

With regard to cortisol levels, in the current study, no differences in morning or afternoon cortisol levels were identified in children with or without ASD, indicating no overall alterations in HPA axis reactivity, as demonstrated before (21, although see 52). Notably, depending on diagnosis, differential couplings between the oxytocinergic system and the HPA axis (as indexed with cortisol) were evident, indicating that only in the TD children, not in children with ASD, heightened (morning) oxytocin levels were predictive of dampened ‘state-dependent’ release of cortisol. Particularly, this negative coupling was strongest between morning oxytocin and afternoon cortisol levels measured after the mildly stress-inducing social interaction session. These results therefore suggest that higher ‘trait-dependent’ levels of morning oxytocin, seen in TD children, may constitute an important *protective* mechanism to buffer HPA axis activity upon (social) stressors encountered during the day. In children with ASD, on the other hand, this protective oxytocinergic HPA-stress attenuation mechanism may be dysfunctional or absent, considering the overall diminished levels of morning oxytocin in the children with ASD, and the apparent lack of coupling between oxytocin and cortisol systems in this group.

Overall, our observations of oxytocinergic dampening of HPA - cortisol responses are in line with prior studies showing acute, dampening effects of intranasal oxytocin administration on cortisol reactivity to laboratory tasks 53,54. Interestingly, meta-analytic evidence showed that the attenuating effects were most pronounced in response to challenging laboratory tasks that produced a robust stimulation of the HPA axis, as well as in clinical populations with ASD or schizophrenia, compared to neurotypicals 53. This latter observation is important, as the authors interpreted this finding to reflect an increased sensitivity to the exogenously administered oxytocin among those with a clinical diagnosis, presumably because they may present lower baseline oxytocin levels. Albeit indirectly, our current observations corroborate this interpretation by showing that morning levels of salivary oxytocin are indeed significantly diminished in children with ASD. Exogenously administered oxytocin may therefore be particularly effective in ASD, i.e., for restoring the otherwise dysfunctional oxytocinergic HPA stress attenuation.

Notably, while diminished levels of oxytocin were initially noted in children with ASD, a significant increase was observed from AM to PM, reaching a similar level of concentration as noted in TD children. In TD children, levels remained overall stable throughout the day. While speculative, the steep increase in oxytocin release in ‘state-dependent’ afternoon oxytocin, compared to ‘trait-dependent’ morning oxytocin, can be hypothesized to reflect a compensatory mechanism of *reactive* oxytocin release upon the mildly stressful social interaction session. Indeed, while absolute oxytocin and cortisol levels were not significantly associated in children with ASD, exploratory analyses confirmed that, in the ASD group, higher state-dependent, afternoon cortisol levels were paralleled by a steeper morning to afternoon rise in oxytocin levels, likely reflective of a coupled co-release of cortisol and oxytocin upon stress-inducing tasks, as demonstrated before in both animals 55 and humans 34. In this view, the parallel release of oxytocin upon HPA activation may be considered a *reactive* stress-regulatory mechanism, which - over time - may facilitate a dampening and further reduction of the cortisol-induced stress response. As highlighted by Brown et al. (2016), positive associations between oxytocin and cortisol levels upon acute stress induction are thought to be reflective of an initial co-release, for subsequently dampening stress responses and facilitating coping behaviors 34. These and our findings therefore highlight the importance of studying the interaction among oxytocinergic and HPA stress systems from a dynamic, time-varying perspective, rather than studying static levels of a single marker.

To sum up, in children with ASD, the *positive* association between the steep rise in oxytocin levels and stress-induced HPA-cortisol levels may be indicative of a *reactive* stress-regulatory mechanism. In TD children, on the other hand, the identified *negative* association between trait (morning) levels of oxytocin and stress-induced HPA-cortisol levels are likely reflective of a *protective* oxytocinergic stress-regulatory mechanism.

Aside hormonal characterizations, also variations in DNAm frequencies of distinct CpG sites of *OXTR* were assessed, but no significant group differences were identified. These results accord to observations of a recent study by Siu et al. (2021), neither showing significant group differences in the level of methylation of CpG sites -934 and -924 in children and adolescents with and without ASD 30. Despite some variation, a recent review pointed towards differential patterns of hypo-versus hyper-methylation, respectively in children and adults with ASD, although the number of included studies was limited, and overall sample sizes were small to modest 26. Also, considerable design-related variations were noted by the authors, *e*.*g*. in terms of included CpG sites, rendering it difficult to draw conclusive interpretations regarding ASD diagnosis-related alterations in *OXTR* DNAm based on the existing evidence.

In terms of relationships between hormonal levels and *OXTR* DNAm, only the TD group displayed a significant negative association between *OXTR* DNAm (at CpG site -924) and state-dependent (afternoon) cortisol levels, indicating that lower methylation frequencies were associated with higher cortisol levels. This finding accords to a prior study identifying a similar negative relationship between stress-induced cortisol and *OXTR* DNAm frequency in individuals with social anxiety disorder (SAD) 56. In the latter study, the pattern of *OXTR* hypomethylation associated with SAD was interpreted to be reflective of a compensatory mechanism, i.e., facilitating an upregulation of oxytocin receptors to cope with heightened HPA axis activity in SAD. More research is needed however, to corroborate this interpretation and pattern of results in future studies.

While the current study provides important new insights into the oxytocinergic system of children with ASD and its interactions with HPA stress reactivity, the following limitations and recommendations are noted. First, as indicated, the included children reflected a rather homogenous group of high-functioning children within a tight pre-pubertal age range, rendering generalizability of the identified effects to more heterogeneous and/or younger/older children or adolescents uncertain. Also, while the included number of boys and girls in our sample reflected the well-documented four-to-one male bias in ASD prevalence 57, future research is warranted to examine the observed effects also in larger samples of girls with ASD, especially given prior reports of gender-related differences in oxytocinergic function 58. Further, while the inclusion of multiple salivary sampling time points (both in the morning and afternoon) was a strong asset for examining dynamic trait (morning) versus state-dependent (after experimental manipulation) hormonal levels, it is noted that any future design would benefit from amplifying the sampling frequency even further, allowing an even more fine-grained sampling of dynamic changes in hormonal systems interactions 34.

To conclude, this study provides new insights into altered oxytocinergic signaling in children with ASD and its interactions with the cortisol HPA stress system. Results provide important indications that only ‘trait-dependent’ oxytocin levels, but not (social) stress-induced ‘state-dependent’ levels of oxytocin are altered in ASD. These observations suggest that TD children without ASD may rely more on a *protective* stress-regulatory mechanism, involving high circulating trait levels of oxytocin for buffering stress reactivity, whereas children with ASD may rely more heavily on a *reactive* stress regulatory mechanism, involving more pronounced stress-induced, state-dependent releases of oxytocin upon HPA - cortisol activation, i.e., for reactively dampening stress responses and facilitating coping behaviors. While more work is needed, the identified results may aid in establishing relevant biomarkers for diagnostic and/or treatment evaluation purposes, as well as for the development of new (personalized) treatment approaches for targeting the oxytocinergic system in ASD.

## Supporting information

Supplementary material [[supplements/283521_file03.pdf]](pending:yes)

## Data Availability

All data produced in the present study are available upon reasonable request to the authors

## Disclosure

The authors declare no conflict of interest.

## Acknowledgments

This research was supported by internal funding of the KU Leuven (ELG-D2857-C14/17/102), a Doctor Gustave Delport fund of the King Baudouin Foundation, a Branco Weiss fellowship of the Society in Science - ETH Zurich, and an EOS Excellence of Science grant (G0E8718N, HUMVISCAT) granted to KA and BB. ME is supported by an FWO aspirant fundamental fellowship [11N1222N]. JP is supported by an FWO junior postdoctoral fellowship [1257621N].

The funding sources had no further role in study design; in the collection, analysis and interpretation of data; in the writing of the report; and in the decision to submit the paper for publication.

## Footnotes

*   * Joint first authors

*   † Joint senior authors

*   Received December 15, 2022.
*   Revision received December 15, 2022.
*   Accepted December 16, 2022.


*   © 2022, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.Heather, M. Clinical observations of the differences between children on the autism spectrum and those with attachment problems: The Coventry Grid. Good Autusm Practice 11, 46–59 (2010).
    
    

2.  2.McVey, A. J. The neurobiological presentation of anxiety in autism spectrum disorder: A systematic review. Autism Res 12, 346–369 (2019).
    
    

3.  3.Marotta, R. et al. The Neurochemistry of Autism. Brain Sci 10, (2020).
    
    

4.  4.Jones, C., Barrera, I., Brothers, S., Ring, R. & Wahlestedt, C. Oxytocin and social functioning. Dialogues Clin Neurosci 19, 193–201 (2017).
    
    

5.  5.Fineberg, S. K. & Ross, D. A. Oxytocin and the Social Brain. Biol Psychiatry 81, e19 (2017).
    
    

6.  6.Jurek, B. & Neumann, I. D. The Oxytocin Receptor: From Intracellular Signaling to Behavior. Physiol Rev 98, 1805–1908 (2018).
    
    [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%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 

7.  7.Bartz, J. A., Zaki, J., Bolger, N. & Ochsner, K. N. Social effects of oxytocin in humans: context and person matter. Trends Cogn Sci 15, 301–309 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tics.2011.05.002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21696997&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000293436600004&link_type=ISI) 

8.  8.Bakermans-Kranenburg, M. J. & van IJzendoorn, M. H. Sniffing around oxytocin: review and meta-analyses of trials in healthy and clinical groups with implications for pharmacotherapy. Transl Psychiatry 3, 1–14 (2013).
    
    

9.  9.Peñagarikano, O. et al. Exogenous and evoked oxytocin restores social behavior in the Cntnap2 mouse model of autism. Sci Transl Med 7, (2015).
    
    

10. 10.Silverman, J. L. et al. Low stress reactivity and neuroendocrine factors in the BTBR T+tf/J mouse model of autism. Neuroscience 171, 1197–1208 (2010).
    
    

11. 11.Schaller, F. et al. A single postnatal injection of oxytocin rescues the lethal feeding behaviour in mouse newborns deficient for the imprinted Magel2 gene. Hum Mol Genet 19, 4895–4905 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddq424&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20876615&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000284642200012&link_type=ISI) 

12. 12.Moerkerke, M. et al. Endogenous Oxytocin Levels in Autism-A Meta-Analysis. Brain Sci 11, (2021).
    
    

13. 13.Zhang, H. F. et al. Plasma Oxytocin and Arginine-Vasopressin Levels in Children with Autism Spectrum Disorder in China: Associations with Symptoms. Neurosci Bull 32, 423–432 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12264-016-0046-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27342432&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 

14. 14.Feldman, R., Golan, O., Hirschler-Guttenberg, Y., Ostfeld-Etzion, S. & Zagoory-Sharon, O. Parent-child interaction and oxytocin production in pre-schoolers with autism spectrum disorder. Br J Psychiatry 205, 107–112 (2014).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImJqcHJjcHN5Y2giO3M6NToicmVzaWQiO3M6OToiMjA1LzIvMTA3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMTIvMTYvMjAyMi4xMi4xNS4yMjI4MzUyMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

15. 15.Husarova, V. M. et al. Plasma Oxytocin in Children with Autism and Its Correlations with Behavioral Parameters in Children and Parents. Psychiatry Investig 13, 174–183 (2016).
    
    

16. 16.Fujisawa, T. X., Tanaka, S., Saito, D. N., Kosaka, H. & Tomoda, A. Visual attention for social information and salivary oxytocin levels in preschool children with autism spectrum disorders: an eye-tracking study. Front Neurosci 8, (2014).
    
    

17. 17.Jacobson, J. D. et al. Evidence for alterations in stimulatory G proteins and oxytocin levels in children with autism. Psychoneuroendocrinology 40, 159–169 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.psyneuen.2013.11.014&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24485488&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000331921300019&link_type=ISI) 

18. 18.Mariscal, M. G. et al. Blood oxytocin concentration positively predicts contagious yawning behavior in children with autism spectrum disorder. Autism Res 12, 1156–1161 (2019).
    
    

19. 19.Tanaka, S. et al. Participatory Art Activities Increase Salivary Oxytocin Secretion of ASD Children. Brain Sci 10, 1–10 (2020).
    
    

20. 20.Taurines, R. et al. Oxytocin plasma concentrations in children and adolescents with autism spectrum disorder: correlation with autistic symptomatology. Atten Defic Hyperact Disord 6, 231–239 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12402-014-0145-y&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24989441&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 

21. 21.Bakker-Huvenaars, M. J. et al. Saliva oxytocin, cortisol, and testosterone levels in adolescent boys with autism spectrum disorder, oppositional defiant disorder/conduct disorder and typically developing individuals. Eur Neuropsychopharmacol 30, 87–101 (2020).
    
    

22. 22.Miller, M. et al. Oxytocin and vasopressin in children and adolescents with autism spectrum disorders: sex differences and associations with symptoms. Autism Res 6, 91–102 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/aur.1270&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23413037&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000318117500003&link_type=ISI) 

23. 23.Breton, C. & Zingg, H. H. Expression and region-specific regulation of the oxytocin receptor gene in rat brain. Endocrinology 138, 1857–1862 (1997).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/en.138.5.1857&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9112379&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997WU33000013&link_type=ISI) 

24. 24.Kimura, T., Tanizawa, O., Mori, K., Brownstein, M. J. & Okayama, H. Structure and expression of a human oxytocin receptor. Nature 356, 526–529 (1992).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/356526a0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1313946&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1992HN23100060&link_type=ISI) 

25. 25.Kusui, C. et al. DNA methylation of the human oxytocin receptor gene promoter regulates tissue-specific gene suppression. Biochem Biophys Res Commun 289, 681–686 (2001).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/bbrc.2001.6024&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11726201&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000172751600007&link_type=ISI) 

26. 26.Moerkerke, M. et al. Oxytocin receptor gene (OXTR) DNA methylation is associated with autism and related social traits – A systematic review. Res Autism Spectr Disord 85, 101785 (2021).
    
    

27. 27.Gregory, S. G. et al. Genomic and epigenetic evidence for oxytocin receptor deficiency in autism. BMC Med 7, 1–13 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1741-7015-7-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19134183&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 

28. 28.Andari, E. et al. Epigenetic modification of the oxytocin receptor gene: implications for autism symptom severity and brain functional connectivity. Neuropsychopharmacology 2020 45:7 45, 1150–1158 (2020).
    
    

29. 29.Elagoz Yuksel, M., Yuceturk, B., Faruk Karatas, O., Ozen, M. & Dogangun, B. The altered promoter methylation of oxytocin receptor gene in autism. J Neurogenet 30, 280–284 (2016).
    
    

30. 30.Siu, M. T. et al. DNA Methylation of the Oxytocin Receptor Across Neurodevelopmental Disorders. J Autism Dev Disord 51, 3610–3623 (2021).
    
    

31. 31.Dadds, M. R. et al. Methylation of the oxytocin receptor gene and oxytocin blood levels in the development of psychopathy. Dev Psychopathol 26, 33–40 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0954579413000497&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24059811&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 

32. 32.Rubin, L. H. et al. Sex and diagnosis specific associations between DNA methylation of the oxytocin receptor gene with emotion processing and temporal-limbic and prefrontal brain volumes in psychotic disorders. Biol Psychiatry Cogn Neurosci Neuroimaging 1, 141–151 (2016).
    
    

33. 33.de Leon, D. et al. Methylation of OXT and OXTR genes, central oxytocin, and social behavior in female macaques. Horm Behav 126, (2020).
    
    

34. 34.Brown, C. A., Cardoso, C. & Ellenbogen, M. A. A meta-analytic review of the correlation between peripheral oxytocin and cortisol concentrations. Front Neuroendocrinol 43, 19–27 (2016).
    
    

35. 35.Gao, J. et al. Alteration of peripheral cortisol and autism spectrum disorder: A meta-analysis. Front Psychiatry 13, (2022).
    
    

36. 36.Smyth, J. M. et al. Individual differences in the diurnal cycle of cortisol. Psychoneuroendocrinology 22, 89–105 (1997).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0306-4530(96)00039-X&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9149331&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997WW79900003&link_type=ISI) 

37. 37.Miller, R. & Kirschbaum, C. Cultures under stress: A cross-national meta-analysis of cortisol responses to the Trier Social Stress Test and their association with anxiety-related value orientations and internalizing mental disorders. Psychoneuroendocrinology 105, 147–154 (2019).
    
    

38. 38.Badr, H. A. & Zauszniewski, J. A. Kangaroo care and postpartum depression: The role of oxytocin. Int J Nurs Sci 4, 179 (2017).
    
    

39. 39.Lord, C. et al. Autism Diagnostic Observation Schedule: ADOS-2. Western Psychological Services (2012).
    
    

40. 40.Wechsler, D. WISC-V-NL. Wechsler Intelligence Scale for Children, Fifth Edition, Dutch version.. Amsterdam: Pearson Benelux B.V. (2018).
    
    

41. 41.Constantino, J. N. & Gruber, C. P. Social responsiveness scale (2nd. ed.). Manual. Western Psychological Services (2012).
    
    

42. 42.Roeyers, H., Thys, M., Druart, C., de Schryver, M. & Schittekatte, M. SRS-2: Screeningslijst voor autismespectrumstoornissen. Hogrefe Uitgevers (2015).
    
    

43. 43.Bodfish, J. W., Symons, F. J., Parker, D. E. & Lewis, M. H. Varieties of repetitive behavior in autism: comparisons to mental retardation. J Autism Dev Disord 30, 237–243 (2000).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1023/A:1005596502855&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11055459&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000089674100006&link_type=ISI) 

44. 44.Luis Royo, J., Hidalgo, M. & Ruiz, A. Pyrosequencing protocol using a universal biotinylated primer for mutation detection and SNP genotyping. Nat Protoc 2, 1734–1739 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nprot.2007.244&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17641638&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000253139200016&link_type=ISI) 

45. 45.Lebowitz, E. R. et al. Salivary oxytocin in clinically anxious youth: Associations with separation anxiety and family accommodation. Psychoneuroendocrinology 65, 35–43 (2016).
    
    

46. 46.Feldman, R., Gordon, I. & Zagoory-Sharon, O. The cross-generation transmission of oxytocin in humans. Horm Behav 58, 669–676 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.yhbeh.2010.06.005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20558167&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000282469700015&link_type=ISI) 

47. 47.Alaerts, K., Steyaert, J., Vanaudenaerde, B., Wenderoth, N. & Bernaerts, S. Changes in endogenous oxytocin levels after intranasal oxytocin treatment in adult men with autism: An exploratory study with long-term follow-up. Eur Neuropsychopharmacol 43, 147–152 (2021).
    
    

48. 48.Krol, K. M., Puglia, M. H., Morris, J. P., Connelly, J. J. & Grossmann, T. Epigenetic modification of the oxytocin receptor gene is associated with emotion processing in the infant brain. Dev Cogn Neurosci 37, 100648 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.dcn.2019.100648&link_type=DOI) 

49. 49.Abdulamir, H. A., Abdul-Rasheed, O. F. & Abdulghani, E. A. Low oxytocin and melatonin levels and their possible role in the diagnosis and prognosis in Iraqi autistic children. Saudi Med J 37, 29–36 (2016).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic21qIjtzOjU6InJlc2lkIjtzOjc6IjM3LzEvMjkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8xMi8xNi8yMDIyLjEyLjE1LjIyMjgzNTIxLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

50. 50.Yang, C. J. et al. The cortisol, serotonin and oxytocin are associated with repetitive behavior in autism spectrum disorder. Res Autism Spectr Disord 18, 12–20 (2015).
    
    

51. 51.Althaus, M. et al. Do blood plasma levels of oxytocin moderate the effect of nasally administered oxytocin on social orienting in high-functioning male adults with autism spectrum disorder? Psychopharmacology (Berl) 233, 2737–2751 (2016).
    
    

52. 52.Tordjman, S. et al. Altered circadian patterns of salivary cortisol in low-functioning children and adolescents with autism. Psychoneuroendocrinology 50, 227–245 (2014).
    
    

53. 53.Cardoso, C., Kingdon, D. & Ellenbogen, M. A. A meta-analytic review of the impact of intranasal oxytocin administration on cortisol concentrations during laboratory tasks: moderation by method and mental health. Psychoneuroendocrinology 49, 161–170 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.psyneuen.2014.07.014&link_type=DOI) 

54. 54.Cardoso, C., Ellenbogen, M. A., Orlando, M. A., Bacon, S. L. & Joober, R. Intranasal oxytocin attenuates the cortisol response to physical stress: a dose-response study. Psychoneuroendocrinology 38, 399–407 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.psyneuen.2012.07.013&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22889586&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000316775400010&link_type=ISI) 

55. 55.Neumann, I. D. & Landgraf, R. Balance of brain oxytocin and vasopressin: implications for anxiety, depression, and social behaviors. Trends Neurosci 35, 649–659 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tins.2012.08.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22974560&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000311177000001&link_type=ISI) 

56. 56.Ziegler, C. et al. Oxytocin receptor gene methylation: converging multilevel evidence for a role in social anxiety. Neuropsychopharmacology 40, 1528–1538 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/npp.2015.2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25563749&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 

57. 57.Fombonne, E. Epidemiology of pervasive developmental disorders. Pediatr Res 65, 591–598 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1203/PDR.0b013e31819e7203&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19218885&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000266271800001&link_type=ISI) 

58. 58.Carter, C. S. Sex differences in oxytocin and vasopressin: implications for autism spectrum disorders? Behavioural brain research 176, 170–186 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbr.2006.08.025&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17000015&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F16%2F2022.12.15.22283521.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000243456800015&link_type=ISI)