Protocol | Open | Open Peer Review | Published:
The relationship between central and peripheral oxytocin concentrations: a systematic review and meta-analysis protocol
Systematic Reviewsvolume 5, Article number: 49 (2016)
Research examining the effects of oxytocin (OT) interventions on psychiatric, social-behavioral, and social-cognitive outcomes regularly collect peripheral levels of OT as markers of central bioavailability. Such inferences rest on the assumption that central and peripheral levels of OT are directly associated. However, conflicting evidence from coordinated sampling of central and peripheral OT question the validity of this assumption. The purpose of this meta-analysis is to evaluate the correlation between central and peripheral OT, as well as to account for potential heterogeneity in the literature.
Studies that report coordinated sampling of central and peripheral OT in humans and animals will be identified. Research investigating concentrations in both basal states and after exogenous administration will be considered. PubMed and Embase databases will be searched, along with citation lists of retrieved articles. Peer-reviewed studies written in English published from 1971 onwards will be included in the meta-analysis. Data will be extracted from eligible studies for a random-effects meta-analysis. For each study, a summary effect size, heterogeneity, risk of bias, publication bias, and the effect of categorical and continuous moderator variables will be determined.
This systematic review and meta-analysis will identify and synthesize evidence to determine if there is an association between central and peripheral OT concentrations. If significant associations are observed, evidence would provide a rationale for future research to use peripheral measures as a proxy for central OT concentrations.
Systematic review registration
Oxytocin (OT) is a neuropeptide that is primarily synthesized in the paraventricular and supraoptical nuclei of the hypothalamus. OT plays a vital role in facilitating a range of physiological functions, such as labor induction and lactation . However, recent research interest has focused more on its effects on social cognitive functions including emotion recognition, trust, and intersubjective selectivity [2–4]. Alongside these efforts determining the role of OT in typical behavior and cognition, research began to investigate the efficacy of this neuropeptide in the treatment of psychiatric disorders—such as autism spectrum disorder  and schizophrenia —and the overall role of the OT system in psychopathology.
Three broad research approaches have been adopted to better understand the role of OT in cognition, behavior, and psychiatric illness: (i) the comparison of a psychological outcome measure or neurobiological effect after OT and placebo treatment; (ii) the comparison of basal OT concentrations between a pathological (e.g., psychiatric) and non-pathological control group; and (iii) the assessment of basal OT concentration covariance with various psychological or other biological phenomena. While central to the latter two methods in this list, peripheral concentrations of OT (e.g., blood plasma, saliva, urine) have been assessed within each of these research traditions, typically in an effort to approximate central OT concentrations and/or central bioavailability after OT administration. Using peripheral OT concentrations to index central concentrations is obviously appealing given the invasive methods required to collect centrally circulating fluids in humans, such as cerebrospinal fluid (CSF). However, the relationship between central and peripheral OT concentrations remains unclear. Although some animal research indicates that central release from the hypothalamus and peripheral release via the posterior pituitary is coordinated [7–9], other research does not support this [10, 11]. Research is also mixed in humans, with some results consistent with related levels of central and peripheral endogenous OT , while others report no significant associations . In relation to exogenous OT, one human study found that peripheral and central levels were not related after intranasal administration; however, as the authors of the study acknowledge, the small sample size suggest caution in interpreting these results . Together, it appears unclear if peripheral OT measures can be used as a proxy for central nervous system (CNS) concentrations and central OT bioavailability.
Meta-analysis enables statistical pooling of effect sizes and is a valuable tool for synthesizing results across studies and increasing power. In the proposed meta-analysis, we will examine studies in which both central and peripheral measures of OT have been sampled. This will enable greater understanding of the dynamics of neuropeptide release and distribution and help determine whether peripheral concentration is a reliable marker for central concentration.
As a primary outcome, an overall summary effect size for the relationship between central and peripheral concentrations of OT will be calculated from the effect sizes of all the individual studies included in the review. Since central and peripheral samples of OT have been obtained simultaneously in several species, the relevant population will not be limited to humans. This summary effect size will also generalize across experimental paradigms and interventions. Additional primary outcomes include summary effect sizes for the relationship between central and peripheral concentrations of OT assessments within species category and within experimental paradigm category (see below for details). The impact of other moderators will be treated as secondary outcomes.
The present study is a meta-analysis of available studies where coordinated central and peripheral measures of OT have been sampled. The aim is to examine the relationship between central and peripheral OT levels, as well as to provide a rationale for discrepancies between the results of independent studies. If a significant overall correlation is observed, then peripheral oxytocin concentrations may represent an index of central oxytocin levels. This protocol is registered with PROSPERO (CRD42015027864) and has been reported here according to the PRISMA-P  guidelines [see Additional file 1]. The procedures described in this protocol follow recent recommendations for the meta-analysis of correlational data .
We will conduct a systematic literature search in three iterations to collect studies that simultaneously measure central and peripheral concentrations of OT. The specificity and sensitivity of different search criteria were examined in a pilot search. In the first iteration, searches will be performed in PubMed and Embase with the following combination of terms: (oxytocin) AND (concentration* OR level*) AND (plasma OR blood OR saliva* OR urin*) AND (central OR csf OR "cerebrospinal fluid"). Further limits for the search will be set in order to focus the search on full text articles that were published between January 1, 1971, and the date in which the search will be performed. In a second iteration, review articles will be consulted for relevant citations. Since any combination of terms is unlikely to exhibit perfect sensitivity to relevant articles, and review articles are unlikely to review a complete set of relevant articles, a third iteration will be performed in which reference lists within studies and citing articles of selected studies in the first and second iterations will be examined for remaining studies that include the critical measures.
Inclusion and exclusion criteria
In this meta-analysis, we will include studies that meet the following criteria: (a) OT sampled in fluids with access to central (e.g., local extracellular fluid, cerebrospinal fluid) and peripheral (e.g., blood plasma, urine, saliva) regions of the body; (b) correlations between central and peripheral concentrations obtainable either through the article, the authors directly, statistical transformation of available data (e.g., Spearman to Pearson correlation), or “data scraping” from available scatterplots; (c) article published in a peer-reviewed journal between 1971 (discovery of ELISA) and June 2015; and (d) written in English.
Studies that provide simultaneous measures of central and peripheral OT concentrations are not only heterogeneous by the correlations they report, they also differ from each other on several potentially critical methodological aspects. These include species examined, participant characteristics, experimental paradigm, sampling type, biochemical analysis of OT concentration, as well as whether or not extraction of neuropeptide from sampled substance was performed. The following potential moderator variables will be examined a priori in this meta-analysis in order to account for heterogeneity in the literature.
Participants and species. Participants include clinical and non-clinical human populations. Animal studies will also be included to collect a broader dataset of studies assessing peripheral and central OT concentrations, with species encompassing non-human primates and other mammals. It is possible that studies in non-human animals may yield results differing from those in human studies given differences in physiology to humans. For example, the olfactory epithelium covers as much as 50 % of the rodent nasal cavity  whereas the proportion in humans is only 3–15 % [18, 19]. As intranasal administration of OT appears to depend on uptake through the olfactory epithelium for direct entry into the brain , such differences in physiology moderate effect sizes. Furthermore, within human studies, we will examine whether effect sizes are moderated by (a) gender of participants, (b) presence or absence of physical illness, and (c) presence or absence of psychiatric diseases.
Experimental paradigms. These range from the examination of baseline levels of central and peripheral OT [13, 21], and levels following some environmental stimulus such as light , and stressors [7, 23], to examinations of concentrations after manipulation of OT levels by osmotic stimulation , drug administration [24–26], or the administration of intranasal [14, 27] or intraperitoneal  OT.
Since the pharmacokinetics of OT seem to exhibit quite complex, time-dependent patterns , there is a possibility that differences in experimental designs account for differences in correlations between studies. For this moderator, (a) baseline samples, (b) samples after exogenous intranasal OT administration, (c) samples after intravenous or intraperitoneal OT administration, (d) samples after intracerebral OT administration, and (e) samples after other experimental manipulations (e.g., osmotic stimulation) will be treated as levels within the experimental paradigms moderator variable. In studies that report both baseline correlation and correlations after manipulation, the baseline sample will be assigned to the baseline level in the moderator analysis, while experimental samples will be assigned to experimental levels.
Sampling types. These include CSF concentrations and concentrations in central extracellular fluids, and measures of concentrations in peripheral blood, urine, and saliva. A number of discrepancies have been seen in reports of OT concentration measures. For example, urine samples have been reported to give OT concentration estimates at a million times above usual values for plasma samples . Furthermore, in a study on OT levels in parents interacting with their infants, Feldman et al.  found no correlation between urine concentrations of OT and concentrations in blood plasma or saliva. In studies that administer OT intranasally, large amounts of the neuropeptide may be drawn through the floor of the nasal vault, missing nose-to-brain targets in the upper posterior regions, down to the oral mucosa and gastrointestinal tract. In these cases, it is difficult to distinguish between exogenous OT that has dripped down into the oral cavity and salivary OT that represent systematically circulating levels. Indeed, saliva OT has been examined for its validity as a biomarker, with negative results reported .
Method for biochemical analysis of OT concentration. Two main methods for biochemical detection of sampled substance are reported in the literature: enzyme immunoassay (EIA) and radioimmunoassay (RIA). In a validation procedure directed at these methods, Christensen et al.  found that RIA was more sensitive to low concentrations of OT and that concentration estimates provided by RIA corresponded more closely to estimates by highly specific two-dimensional liquid chromatography-tandem mass spectrometry assay [32, 33]. Thus, studies that employ EIA might exhibit increased amounts of noise and accordingly weakened correlations between central and peripheral oxytocin. Another related concern regarding biochemical analysis of OT concentration is whether or not neuropeptides are extracted from the sampled substance prior to EIA or RIA.
Extraction of neuropeptides from the sampled substance is not universal among studies [31, 34]. Measures taken from raw substance are more susceptible to cross-sensitivity, resulting in values that reflect variability of unrelated molecules in addition to the true level of oxytocin [34, 35]. Szeto et al.  found estimates of OT concentrations using EIA to be two orders of magnitude higher in non-extracted plasma as compared to extracted, additionally finding no significant association between OT level estimates from extracted and non-extracted plasma. Thus, there is a possibility that studies where substances are extracted prior to immunoassay provide more accurate correlations between central and peripheral oxytocin if collected simultaneously. In addition, Szeto et al.  found that while EIA estimates increased a 100-fold in extracted plasma as compared to non-extracted, extraction did not increase sensitivity of biochemical analysis using RIA. There seems, in other words, to be an interaction effect of analysis method and extraction on OT concentration estimates. Accordingly, there might be an interaction effect of analysis method and extraction on estimates of correlations between central and peripheral OT concentrations. Therefore, in this meta-analysis, the moderator biochemical analysis method contains two variables with two levels each (RIA vs EIA, and extraction prior to immunoassay or not), and we will examine the main effects of both variables as well as interaction effects.
Age. CSF concentrations of OT seem to be related to age in macaque monkeys, decreasing with age in infants while increasing with age in mothers . Changes in the oxytocin system with age [38, 39] may also involve the extent of coordinated release or even the degree to which OT may penetrate the blood brain barrier.
Risk of bias. A risk of bias tool has been developed to quantify potential biases in effect sizes related to study quality (see below).
Year of publication. Early, more preliminary studies, may potentially exhibit different effect sizes, so year of publication is included as an additional moderator to assess bias.
Level of sample coordination. In some studies, central and peripheral samples are collected simultaneously, whereas in others, there is a brief interval between central and peripheral samples . In this meta-analysis, levels of sample coordination will be coded as simultaneous or non-simultaneous.
Quality assessment (risk of bias)
A custom tool will be developed to assess risk of bias via a systematic procedure for determining the quality of eligible studies [see Additional file 2]. The risk of bias tool is adapted from a previously used tool . Using this tool, studies will be given a score of potential bias by two raters. This tool, in combination with publication bias (described below), will be used to assess the quality of overall evidence. The evidence will be categorized as high (no evidence of publication bias and mean quality scores are ≥80 %), moderate (little-to-no publication bias and quality scores are ≥50 and <80 %), low (evidence of publication bias and quality score ≥50 %), or very low (evidence of publication bias and quality score <50 %).
Data extraction and management
A summary effect size for all studies will be computed with effect sizes (correlations transformed to Fisher’s z) and sample sizes from individual studies as input. Effect sizes will be extracted using one out of four possible procedures, depending on whether and in what form effect sizes are reported in a study. First, if effect sizes are reported, these will be transformed to Fisher’s z. Secondly, if effect sizes are not reported but individual results are reported numerically, r will be calculated and transformed to z. Thirdly, in studies published after 2000 (as 15 years is a common time frame for the retention of clinical data) that do not report effect sizes or individual results in numerical form, authors will be contacted and asked to provide the relevant information. Fourth, in studies either published before 2000 or where authors do not respond to our request, and where all individual results are reported in scatterplot graphs, a plot digitizer (http://arohatgi.info/WebPlotDigitizer) will be used for conversion into numerical values, from which r will be calculated and transformed to z. Studies that do not fall into one of these four categories do not provide the data that is necessary for our analysis and will not be included in the meta-analysis.
Two independent reviewers will extract data from all eligible studies using the data extraction form [see Additional file 3]. This form includes (a) general information on studies including authors, title, number of effect sizes, sample size, and effect sizes; (b) information about the participants including species, gender, age, and physical and mental health status; (c) information about the level the study is on other moderator variables, including study type, biochemical analysis, and neuropeptide extraction; as well as (d) information concerning study quality including publication year and the risk of bias measures as defined by the custom risk of bias tool. Any disagreements will be adjudicated by author DSQ.
The data will be entered into an excel version of the data extraction form, with one excel sheet for each reviewer. For each study, the extracted data will be entered in a separate tab of the data extraction excel sheet. After the data extraction phase of this meta-study, the data will be saved as .csv files and treated as data frames in R for analysis. Data from studies initially selected based on title and abstract and articles included in the review will be documented. Reasons for the exclusion of retrieved articles will also be recorded for eventual documentation in a study search and data extraction flow diagram.
Analysis will be performed with R statistical software version 3.1.1. , using the metafor  and robumeta  packages. All effect sizes will be converted to Fisher’s z prior to analysis, and back to Pearson’s r after analysis. Operations needed for the analysis include the following:
Summary effect size
Since the designs of individual studies vary substantially (e.g., differing species, biochemical methods), real effect sizes of individual studies should vary accordingly. Thus, a random effects model is appropriate in the computation of summary effect size . Summary effect size as well as confidence interval and statistical significance of effect size will be computed in R. In the case of studies reporting more than one correlation, Robust variance estimation will be used with adjusted estimators for small samples (n < 40), if necessary.
The Q-statistic is the sum of the products of weights and squared differences between study effect sizes and summary effect sizes. As such, it is an expression of the heterogeneity of effect sizes in the sample. A significant Q-statistic is indicative of significantly different effect sizes between the studies included in the meta-analysis. In this meta-analysis, we follow the convention of an alpha level of .05 for the Q-statistic. I 2 reflects the proportion of total variance that is attributable to real between studies variance. The I 2 statistic is useful for determining whether the amount of real between studies effect size variance is relatively higher than chance variability. The Q-statistic, the significance of the Q-statistic, and I 2 will be computed and reported.
Categorical moderator analysis
A random effects model with separate estimates of between studies variance will be applied for all categorical moderator variables (species, paradigm, sampling type, biochemical method, and extraction). Summary mean effects are calculated for each subgroup, as described in the previous section for summary effect sizes. T 2 is calculated separately for each subgroup. Means are compared using a two-tailed z-test in order to determine the probability for the difference of observed means given equality of true means, assuming normal distribution. A p value below .05 will be seen as indicative of moderating properties of the variable in question.
Continuous moderator analysis
For the continuous moderator variables (age, risk of bias, and year of publication), meta-regression will be performed to estimate an unstandardized regression coefficient along with the coefficient’s significance level.
Publication bias will be assessed using the regression test for funnel plot asymmetry outlined in Egger et al. . If there is evidence of publication bias, we will use the Trim and Fill method, which provides a reasonable approximation of “missing” effect sizes . With this method, the asymmetric right side of the funnel (which is the side expected to be effected by publication bias in the present meta-analysis) is “trimmed”. The remaining symmetric effect sizes are then used to re-estimate the center of the funnel, after which the trimmed effect sizes and their missing counterpart effect sizes on the left side of the funnel are “filled”. An estimate of the summary effect size and variance will then be calculated on the filled funnel plot as a form of sensitivity analysis to examine the potential impact of missing studies.
There is increasing research on the use of peripheral OT concentrations to determine relationships to a range of psychological and neurobiological phenomena. Moreover, levels of peripheral OT are also used to determine the bioavailability of OT after intranasal administration. These measures are often used as a proxy for central levels of OT but it is unknown if these two measures are in fact related, with mixed results reported in the literature across species and experimental designs. Recent debate has also highlighted the need for a more thorough characterization of OT mechanisms [47, 48] which includes measures of central and peripheral OT bioavailability.
This will be the first systematic review and meta-analysis, to the best of our knowledge, to investigate this association. An up-to-date synthesis of the available data will help determine if future research can use peripheral levels of OT as a measure of central activity or if the collection of CSF concentration is necessary to make these conclusions. Alternatively, the proposed meta-analysis may indicate that central and peripheral OT concentrations are only related under certain conditions, such as after exogenous OT administration. Predetermined measures of publication bias and study quality will be used to describe the confidence in the resulting body of evidence.
central nervous system
Burbach JP, Young LJ, Russell J. Oxytocin: synthesis, secretion, and reproductive functions. Knobil Neill’s Physiolc Reprod. 2006;2:3055–128.
Guastella AJ, Einfeld SL, Gray KM, Rinehart NJ, Tonge BJ, Lambert TJ, et al. Intranasal oxytocin improves emotion recognition for youth with autism spectrum disorders. Biol Psychiatry. 2010;67(7):692–4.
Kosfeld M, Heinrichs M, Zak PJ, Fischbacher U, Fehr E. Oxytocin increases trust in humans. Nature. 2005;435(7042):673.
De Dreu CKW. Oxytocin modulates cooperation within and competition between groups: an integrative review and research agenda. Horm Behav. 2012;61:419–28.
Modi ME, Young LJ. The oxytocin system in drug discovery for autism: animal models and novel therapeutic strategies. Horm Behav. 2012;61:340–50.
MacDonald K, Feifel D. Oxytocin in schizophrenia: a review of evidence for its therapeutic effects. Acta Neuropsychiatrica. 2012;24:130–46.
Wotjak CT, Ganster J, Kohl G, Holsboer F, Landgraf R, Engelmann M. Dissociated central and peripheral release of vasopressin, but not oxytocin, in response to repeated swim stress: new insights into the secretory capacities of peptidergic neurons. Neuroscience. 1998;85(4):1209–22.
Ross HE, Cole CD, Smith Y, Neumann ID, Landgraf R, Murphy AZ, et al. Characterization of the oxytocin system regulating affiliative behavior in female prairie voles. Neuroscience. 2009;162(4):892–903.
Landgraf R, Neumann I, Schwarzberg H. Central and peripheral release of vasopressin and oxytocin in the conscious rat after osmotic stimulation. Brain Res. 1988;457(2):219–25.
Amico JA, Challinor SM, Cameron JL. Pattern of oxytocin concentrations in the plasma and cerebrospinal fluid of lactating rhesus monkeys (Macaca mulatta): evidence for functionally independent oxytocinergic pathways in primates. J Clin Endocrinol Metabol. 1990;71(6):1531–5. http://doi.org/10.1210/jcem-71-6-1531.
Robinson IC, Jones PM. Oxytocin and neurophysin in plasma and CSF during suckling in the guinea-pig. Neuroendocrinology. 1982;34(1):59–63.
Carson DS, Berquist SW, Trujillo TH, Garner JP, Hannah SL, Hyde SA, et al. Cerebrospinal fluid and plasma oxytocin concentrations are positively correlated and negatively predict anxiety in children. Mol Psychiatry. 2014;20:1085–90. http://doi.org/10.1038/mp.2014.132.
Kagerbauer SM, Martin J, Schuster T, Blobner M, Kochs EF, Landgraf R. Plasma oxytocin and vasopressin do not predict neuropeptide concentrations in human cerebrospinal fluid. J Neuroendocrinol. 2013;25(7):668–73. http://doi.org/10.1111/jne.12038.
Striepens N, Kendrick KM, Hanking V, Landgraf R, Wüllner U, Maier W, et al. Elevated cerebrospinal fluid and blood concentrations of oxytocin following its intranasal administration in humans. Sci Rep. 2013;3. http://doi.org/10.1038/srep03440.
Shamseer L, Moher D, Clarke M, Ghersi D, Liberati A, Petticrew M, et al. Preferred reporting items for systematic review and meta-analysis protocols (PRISMA-P) 2015: elaboration and explanation. BMJ (Clin Res Ed). 2015;349:7647–7. http://doi.org/10.1136/bmj.g7647.
Quintana DS. From pre-registration to publication: a nontechnical primer for conducting a meta-analysis to synthesize correlational data. Front Psychol. 2015;6:1549. http://doi.org/10.3389/fpsyg.2015.01549.
Illum L. Nasal delivery. The use of animal models to predict performance in man. J. Drug Target. 1996;3:427–42.
Morrison EE, Costanzo RM. Morphology of olfactory epithelium in humans and other vertebrates. Microsc Res Tech. 1992;23:49–61.
Popp JA, Monteiro-Riviere NA. Macroscopic, microscopic, and ultrastructrural anatomy of the nasal cavity, rat. Respir Syst. 1985:3–10
Thorne RG, Pronk GJ, Padmanabhan V, Frey WH. Delivery of insulin-like growth factor-I to the rat brain and spinal cord along olfactory and trigeminal pathways following intranasal aministration. Neuroscience. 2004;127:481–96.
Jokinen J, Chatzittofis A, Hellström C, Nordström P, Uvnäs-Moberg K, Åsberg M. Low CSF oxytocin reflects high intent in suicide attempters. Psychoneuroendocrinology. 2012;37(4):482–90.
Devarajan K, Rusak B. Oxytocin levels in the plasma and cerebrospinal fluid of male rats: effects of circadian phase, light and stress. Neurosci Lett. 2004;367(2):144–7.
Babygirija R, Bülbül M, Yashimoto S, Ludwig K, Takahashi T. Central and peripheral release of oxytocin following chronic homotypic stress in rats. Autonomic Neurosci: Basic Clin. 2012;167(1–2):56–60.
Brown DC, Perkowski SZ, Shofer F, Amico JA. Effect of centrally administered opioid receptor agonists on CSF and plasma oxytocin concentrations in dogs. Am J Vet Res. 2001;62(4):496–9.
Ludwig M, Callahan MF, Morris M. Effects of tetrodotoxin on osmotically stimulated central and peripheral vasopressin and oxytocin release. Neuroendocrinology. 1995;62:619–27.
van Wimersma Greidanus TB, Janssen S, Frankhujzen-Siervogel JC, Maigret C, van de Heijning BJ. Effect of central administration of the kappa-opiate receptor agonist U 69.593 on neurohypophyseal hormone levels in blood and cerebrospinal fluid. Neuropeptides. 1996;30(5):425–5.
Neumann ID, Maloumby R, Beiderbeck DI, Lukas M, Landgraf R. Increased brain and plasma oxytocin after nasal and peripheral administration in rats and mice. Psychoneuroendocrinology. 2013;38(10):1985–93.
Amico J, Ulbrecht J, Robinson A. Clearance studies of oxytocin in humans using radioimmunoassay measurements of the hormone in plasma and urine. J Clin Endocrinol Metab. 1987;64:340–5.
Feldman R, Gordon I, Zagoory-Sharon O. Maternal and paternal plasma, salivary, and urinary oxytocin and parent-infant synchrony: considering stress and affiliation components of human bonding. Dev Sci. 2011;14(4):752–61.
Horvat-Gordon M, Granger DA, Schwartz EB, Nelson VJ, Kivlighan KT. Oxytocin is not a valid biomarker when measured in saliva by immunoassay. Physiol Behav. 2005;84(3):445–8. http://doi.org/10.1016/j.physbeh.2005.01.007.
Christensten JC, Shiyanov PA, Estepp JR, Schlager JJ. Lack of association between human plasma oxytocin and interpersonal trust in a Prisoner’s Dilemma paradigm. PLoS One. 2014;9(12):e116172.
Zhang G, Zhang Y, Fast DM, Lin Z, Steenwyk R. Ultra sensitive quantitation of endogenous oxytocin in rat and human plasma using a two-dimensional liquid chromatography-tandem mass spectrometry assay. Anal Biochem. 2011;416(1):45–5.
Johnsen E, Leknes S, Wilson SR, Lundanes E. Liquid chromatography-mass spectrometry platform for both small neurotransmitters neurotransmitters in blood, with automatic and robust solid phase extraction. Sci Rep. 2015;5:9308.
McCullough ME, Churchland PS, Mendez AJ. Problems with measuring peripheral oxytocin: can the data on oxytocin and human behavior be trusted? Neurosci Biobehav Rev. 2013;37:1485–92.
Robinson KJ, Hazon N, Lonergan M, Pomeroy PP. Validation of enzyme-linked immunoassay (ELISA) for plasma oxytocin in a novel mammal species reveals potential errors induced by sampling procedure. J Neurosci Methods. 2014;226:73–9.
Szeto A, McCabe PM, Nation DA, Tabak BA, Rossetti MA, McCullough ME, et al. Evaluation of enzyme immunoassay and radioimmunoassay methods for the measurement of plasma oxytocin. Psychosom Med. 2011;73(5):393–400.
Parker KJ, Hoffman CL, Hyde SA, Cummings CS, Maestripieri D. Effects of age on cerebrospinal fluid oxytocin levels in free-ranging adult female and infant rhesus macaques. Behav Neurosci. 2010;124(3):428–33.
Goudsmit E, Fliers E, Swaab DF. Vasopressin and oxytocin excretion in the Brown-Norway rat in relation to aging, water metabolism and testosterone. Mech Ageing Dev. 1988;44(3):241–52.
Zbuzek V, Fuchs AR, Zbuzek VK, Wu W. Neurohypophyseal aging: differential changes in oxytocin and vasopressin release, studied in Fischer 344 and Sprague–Dawley rats. Neuroendocrinology. 1988;48(6):619–26.
Alvares GA, Quintana DS, Hickie IB, Guastella AJ. Autonomic nervous system dysfunction in psychiatric disorders and the impact of psychotropic medications: a systematic review and meta-analysis. J Psychiatry Neurosci. 2016;41(2):89–104.
Core Team R. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2014 http://www.R-project.org/.
Viechtbauer W. Conducting meta-analyses in R with the metafor package. J Stat Softw. 2010;36(3):1–48.
Fischer Z, Tipton E. Robumeta: Robust variance meta-regression. R package version 1.6. http://CRAN.R-project.org/package=robumeta. 2015.
Borenstein M, Hedges LV, Higgins JPT, Rothstein HR. Introduction to meta-analysis. UK: John Wiley & Sons; 2009.
Egger M, Smith GD, Scheider M, Minder C. Bias in meta-analysis detected by a simple, graphical test. BMJ. 1997;315:629–34.
Duval S, Tweedie R. Trim and fill: a simple funnel-plot-based method of testing and adjusting for publication bias in meta-analysis. Biometrics. 2000;56(2):455–63.
Leng G, Ludwig M. Intranasal Oxytocin: Myths and Delusions. Biol Psychiatry. 2016;79(3): 243-250. http://doi.org/10.1016/j.biopsych.2015.05.003.
Quintana DS, Woolley JD. Intranasal oxytocin mechanisms can be better understood, but its effects on social cognition and behavior are not to be sniffed at. Biol Psychiatry. 2016;79(8):e49-e50. http://doi.org/10.1016/j.biopsych.2015.06.021.
We thank Hege Kristin Ringnes (University of Oslo Library) for providing guidance on our systematic search strategy, and Maiken Egknud for assisting with the development of the data extraction form. The Research Council of Norway (RCN) and OptiNose AS contributed to funding this review through a BIA project grant (219483) via salary support to DSQ and project support to OAA, LTW, and DSQ. MV received salary support from the RCN via a grant for students in clinical psychology programmes (Studentstipend for profesjonsstudenter i psykologi som skriver hovedoppgave). GAA did not receive funding.
OAA, LTW, and DSQ are investigators in a project studying oxytocin’s effects after intranasal delivery partnered by OptiNose AS (Oslo, Norway) and funded by a User-driven Research based Innovation (BIA) grant (219483) from the Research Council of Norway (RCN). The RCN and partner contributed to funding this review (through salary support to DSQ and project support to OAA, LTW, and DSQ); however, they had no influence in the ideas contained in the manuscript and no role in the writing of the manuscript.
DSQ conceived the study idea. MV, GAA, OAA, LTW, and DSQ contributed to the design of the systematic review. MV, GAA, and DSQ contributed to the data analysis plan. All authors contributed to the write up and editing of the manuscript and approved the final manuscript. DSQ takes responsibility, and acts as guarantor, for the contents of the protocol and review.
PRISMA-P (Preferred Reporting Items for Systematic Review and Meta-Analysis Protocols) 2015 checklist. Recommended items to address in a systematic review protocol. (PDF 165 kb)
Risk of bias tool. Criteria for objective assessment of risk of bias for individual studies included in the meta-analysis. (PDF 37 kb)
Data extraction form. Record of pertinent study characteristics for each included study. (PDF 37 kb)