Social cognitive deficits, including deficits in emotion processing, precede onset of illness, are stable impairments over time and important predictors for functional disabilities in patients with schizophrenia.Reference Corcoran, Keilp, Kayser, Klim, Butler and Bruder 1 , Reference Green, Horan and Lee 2 Add-on trials with oxytocin administered intranasally in this patient group suggest a treatment effect on these deficits as well as a modulating effect on amygdala response to facial emotion processing.Reference Shin, Park, Jung, Park, Yun and Jang 3 , Reference Feifel, Shilling and MacDonald 4 Trials from healthy subjects show an effect on amygdala activation in a possible sexual dimorphic pattern, but the diverse direction of effects seen are not yet established.Reference Shamay-Tsoory and Abu-Akel 5 – Reference Domes, Lischke, Berger, Grossmann, Hauenstein and Heinrichs 7 The results from oxytocin trials are inconsistent in general,Reference Shin, Park, Jung, Park, Yun and Jang 3 , Reference Feifel, Shilling and MacDonald 4 , Reference Heringa, Begemann, Goverde and Sommer 8 and this inconsistency has been attributed to individual differences, including genetically based differences in the oxytocin system,Reference Bartz, Zaki, Bolger and Ochsner 9 , Reference Bartholomeusz, Ganella, Labuschagne, Bousman and Pantelis 10 and a poor understanding of drug targets and mechanisms.Reference Leng and Ludwig 11 The field is also troubled by small sample sizes with an average sample size in schizophrenia and oxytocin imaging of 16 and 32, respectively.Reference Li, Chan, McAlonan and Gong 12 – Reference Bethlehem, van Honk, Auyeung and Baron-Cohen 15 The endogenous oxytocin system is less explored, and the contribution of oxytocin pathway genes to amygdala activation has yet to be investigated in this patient group. Rs53576, rs2254298 and rs237902 are single nucleotide polymorphisms (SNPs) in the gene coding for the production of the oxytocin receptor (OXTR) and have previously been associated with disease characteristics in patients with schizophrenia, that is, general psychopathology (rs53576 and rs2254298) and negative symptoms (rs237902 and rs53576) as measured by the Positive and Negative Syndrome Scale (PANSS) and empathic concern (rs2254298).Reference Montag, Brockmann, Lehmann, Muller, Rujescu and Gallinat 16 – Reference Haram, Tesli, Bettella, Djurovic, Andreassen and Melle 18 The results have, however, also here been inconsistent.Reference Bartholomeusz, Ganella, Labuschagne, Bousman and Pantelis 10 , Reference Haram, Tesli, Bettella, Djurovic, Andreassen and Melle 18 , Reference Bakermans-Kranenburg and van Ijzendoorn 19 Thus, to get a better understanding of the contribution of oxytocin pathway gene polymorphisms to the amygdala response to fearful/angry faces, we used a functional magnetic resonance imaging (fMRI) paradigm for matching of emotional faces in the largest oxytocin clinical sample size to date,Reference Shilling and Feifel 20 consisting of patients with schizophrenia spectrum (SCZ) and affective spectrum disorders (AD) as well as a community representative sample of healthy participants (HC), specifically targeting the SNPs: rs53576, rs2254298 and rs237902. As emotion stimuli processing is particularly altered in patients with schizophrenia, we here hypothesised that the contribution of the three SNPs would be more pronounced in the SCZ group. Because of the reported sex differences in the response to exogenous oxytocin, we additionally investigated interaction effects between SNPs and sex in the case of significant associations between SNPs and amygdala response.
Method
Sample characteristics
Participants, all of Caucasian origin, were included in the Norwegian multicentre Thematically Organized Psychosis (TOP) research study recruiting patients from in- and out-patient clinics in the greater Oslo area in south-eastern Norway. None of the participants had a history of head injury, neurological disorder, autoimmune or infectious disorders or malignancies, or an IQ below 70. Patient assessments included clinical and physical examinations by a physician, neuropsychological testing by a psychologist, collection of blood samples for somatic screening and DNA analyses, and MRI. The Structured Clinical Interview for DSM-IV Axis 1 Disorders (SCID-1) was used for diagnostic purposes. The sample of 104 patients with SCZ included the following diagnoses: schizophrenia (n=62), schizophreniform disorder (n=6), schizoaffective disorder (n=12) and psychotic disorder not otherwise specified (NOS) (n=24). The AD group (n=100) consisted of patients with bipolar disorder I (BD I, n=50), bipolar disorder II (BD II, n=38), bipolar disorder NOS (BD NOS, n=6) and depressive psychosis (n=6) (Table 1). The healthy controls (HC) (n=142) were invited via stratified random selection from statistical records of persons from the same catchment area as the patient groups and excluded in case of the presence of severe mental disorders, including schizophrenia and AD in the participant and in first-degree relatives, or current illegal substance use.
SCZ n=104 | AD n=100 | HC n=142 | Total N=346 | Test statistic | ||
---|---|---|---|---|---|---|
F/χ2 | P | |||||
Gender (m/f) | 63/41 | 43/57 | 84/58 | 190/156 | χ2=8.1 | 0.02 |
Age, years: mean (s.d.) | 31.5 (8.9) | 33.2 (11.2) | 33.4 (8.4) | 32.8 (9.5) | F=1.44 | 0.24 a |
Left amygdala activation (COPE): mean (s.d.) | 11.6 (13.7) | 12.3 (13.1) | 11.0 (12.9) | 11.6 (12.9) | χ2=0.73 | 0.69 b |
Right amygdala activation (COPE): mean (s.d.) | 12.1 (15.4) | 12.6 (13.6) | 12.3 (11.7) | 12.3 (13.4) | χ2=0.19 | 0.91 b |
Medicated c N (%) | 91 (87.5) | 80 (80) | – | – | – | – |
AP/LIT/AE/AD/HYP d (N) | 88/1/21/32/12 | 47/10/39/31/13 | – | – | – | – |
Rs2254298 (maf, A e ) | 0.11 | 0.19 | 0.10 | 0.10 | – | – |
Rs53576 (maf, A) | 0.29 | 0.28 | 0.32 | 0.30 | – | – |
Rs237902 (maf, A) | 0.31 | 0.42 | 0.35 | 0.36 | – | – |
Behavioural data | ||||||
Response time f (ms): mean (s.d.) | 1192.8 (314.9) | 1197.6 (339) | 1069.9 (230.9) | 1141 (295.3) | χ2=13.6 | 0.001 |
Accuracy rate f (%): mean (s.d.) | 98.7 (0.06) | 98.5 (0.09) | 99.2 (0.03) | 98.9 (0.06) | χ2=0.29 | 0.87 |
SCZ, schizophrenia spectrum disorders; AD, affective spectrum disorders; HC, healthy controls; COPE, contrast parameter estimates; maf, minor allele frequency.
a Brown–Forsythe test on log-transformed age.
b Non-parametric Kruskal–Wallis test.
c Regular use of psychopharmacological medication.
d Regular use of antipsychotics/lithium/anti epileptics/anti-depressives/hypnotics.
e The minor allele frequency of the A allele.
f Complete behavioural data (response time and accuracy rate) were available for 89 SCZ, 94 AD and 140 HC. For the remaining individuals (SCZ 15, AD 6, HC 2) an accuracy rate (i.e. combined rate for negative and shapes) was available and was used to confirm that the participants paid attention to the task (accuracy rate 96.2%, 96.3% and 93.5%, respectively).
The study received approval from the Norwegian Scientific-Ethical Committees and the Norwegian Data Protection Agency, and all participants provided written informed consent.
Genotyping
DNA was extracted from blood and rs237902 was genotyped using the Affymetrix Human SNP Array 6.0 (Affymetrix Inc., Santa Clara, CA, USA), as previously reported.Reference Athanasiu, Mattingsdal, Kahler, Brown, Gustafsson and Agartz 21 , Reference Djurovic, Gustafsson, Mattingsdal, Athanasiu, Bjella and Tesli 22 All chips were subjected to the Birdseed-v2 algorithm developed by Affymetrix Inc. and Broad InstituteReference Korn, Kuruvilla, McCarroll, Wysoker, Nemesh and Cawley 23 and implemented in the software Affymetrix Power Tools (APT v1.10). Samples were excluded through quality control if they had low-yield (call rate below 95%), if they were duplicates of other samples, if they had a sex determined by X chromosome marker homozygosity different from their reported sex, or if they were calculated to have ancestry other than European.Reference Athanasiu, Mattingsdal, Kahler, Brown, Gustafsson and Agartz 21 The markers did not show deviation from Hardy–Weinberg equilibrium (P<0.001), low-yield (<95%) or minor allele frequencies below 0.01. Following the above-mentioned quality control, rs53576 and rs2254298 were imputed with MaCH (http://www.sph.umich.edu/csg/abecasis/MACH/download/1000G-PhaseI-Interim.html) using the European samples in the Phase I release of the 1000 Genomes Project.
Experimental paradigm
A widely used and validated paradigm was employed to elicit amygdala activation.Reference Carre, Fisher, Manuck and Hariri 24 – Reference Tesli, Kauppi, Bettella, Brandt, Kaufmann and Espeseth 27 In this task, participants selected which of two stimuli (displayed at the bottom of the screen) matched a target stimulus (displayed at the top). In the faces-matching task, the images displayed were human faces expressing either anger or fear (‘negative faces’), or happiness (‘positive faces’). In the sensorimotor control task, geometrical shapes were displayed. The negative and positive faces were presented in two separate runs, each with four blocks of faces interleaved by five blocks of geometrical shapes. The order of the runs was counterbalanced among participants. Within blocks of negative faces, fearful and angry faces were intermixed. Each block of faces consisted of six emotion-specific face trios derived from a standard set of facial-affected pictures.Reference Tottenham, Tanaka, Leon, McCarry, Nurse and Hare 28 Each trial (face/shape stimulus) was presented for 5.4 s with no inter-stimulus interval, for a total block length of 32.6 s. Each run lasted 294 s. The E-prime software (version 1.0 Psychology Software Tools, Inc, Pittsburgh, PA, USA) controlled the presentations of the stimuli using VisualSystem (NordicNeuroLab, Bergen, Norway). Measures of reaction times and task accuracy were recorded through MRI-compatible ResponseGrips (NordicNeuroLab, Bergen, Norway).
BOLD fMRI data acquisition
MRI was performed on a 1.5 T Siemens Magnetom Sonata scanner (Siemens Medical Solutions, Erlangen, Germany) supplied with a standard head coil. We collected 152 functional volumes per run using a T2*-weighted echo-planar imaging (EPI) sequence. Each volume consisted of 24 axial slices with a pixel size of 3 mm in the axial plane and a slice thickness of 4 mm with 1 mm gap between slices (TR=2040 ms, TE=50 ms, FA=90°, matrix 64×64, FOV 224 mm). The first seven volumes and the last one were discarded, leaving 144 volumes for analyses. Prior to fMRI, a sagittal T1-weighted 3D Magnetization Prepared Rapid Gradient Echo (MPRAGE) scan was collected (160 slices, TR=2730 ms, TE=3.93 ms, FA=7°, matrix 192×256, FOV 240 mm, voxel size 1.33×0.94×1 mm3), here used for co-registration purposes.
fMRI data analysis
The fMRI data went through an initial quality check procedure, to detect images with poor quality because of excessive head motion, slice dropout, radiofrequency (RF) artefacts or other noise. It was subsequently processed and analysed using FEAT (FMRI Expert Analysis Tool) version 6.0, part of FSL (fMRIB's Software Library http://www.fmrib.ox.ac.uk/fsl).Reference Smith, Jenkinson, Woolrich, Beckmann, Behrens and Johansen-Berg 29 Conventional preprocessing included motion correction,Reference Jenkinson, Bannister, Brady and Smith 30 non-brain removal,Reference Smith 31 spatial smoothing with a 6-mm full-width half-maximum kernelReference Smith and Brady 32 and a high-pass filter with a 90-s window. A first-level general linear model (GLM) analysis was performed for each run, where the onset and duration of the on-blocks (positive faces and negative faces, respectively) was modelled with that of the off-blocks (geometrical shapes) as implicit baseline. The design matrix was filtered and convolved with a haemodynamic response function before the model fit. This analysis resulted in individual contrast parameter estimates (COPEs) reflecting the faces versus shapes contrast within runs. In the current study, we focused on participants’ amygdala responses to negative stimuli (the negative faces>shapes contrast), because of the implication of altered threat-related facial affect perception in psychotic disordersReference Kohler, Turner, Bilker, Brensinger, Siegel and Kanes 33 , Reference Hall, Whalley, McKirdy, Romaniuk, McGonigle and McIntosh 34 and since this contrast has been reported to reliably engage the amygdala.Reference Tesli, Kauppi, Bettella, Brandt, Kaufmann and Espeseth 27 , Reference Costafreda, Brammer, David and Fu 35 , Reference Zald 36 Amygdala regions of interest (ROIs) were defined in accordance with the probabilistic Harvard–Oxford subcortical atlas provided with FSL, and were thresholded at 25% probability. Mean BOLD signal changes (parameter estimates) across all ROI voxels were extracted from FSL and used in statistical analyses.
Statistics
Group analyses of amygdala activation and behavioural data
Potential differences in amygdala activation, response time and accuracy rate between subgroups (SCZ, AD, and HC) were examined. The distributions were not normal within all subgroups, hence Kruskal–Wallis tests were performed. Sex effects for amygdala activation were tested for in the total group with Mann–Whitney tests because of values without normal distribution across gender. The correlation between right and left amygdala activation was examined using bivariate correlation analyses (Spearman's rho). These analyses were performed in SPSS version 22 (IBM, Armonk, NY).
Genetic associations analyses
The additive effects of allele dosage on right and left amygdala activation in the total group were investigated for rs237902 (genotyped SNP), as well as rs53576 and rs2254298 (imputed SNPs) with age, sex and diagnosis as covariates in linear regression models. In the follow-up analyses, the genetic association tests for all the three SNPs were repeated in the three subgroups separately, with age and sex as covariates (PLINK Version 1.07). In the case of significant genetic association with amygdala activation, interaction analyses between SNP and diagnosis as well as between SNP and sex in linear regression models were performed (R Version 3.2.2). All analyses were performed with and without outliers. The significance threshold was set to P=0.05/3=0.017 as a result of Bonferroni corrections for three SNPs. Linkage disequilibrium (LD) analyses across the three SNPs were performed (R Version 3.2.2).
Results
Group analyses of behavioural data and amygdala activation
We did not detect differences in amygdala activation (left hemisphere: χ2=0.73, P=0.69; right hemisphere: χ2=0.19, P=0.91) or accuracy rate (χ2=0.29, P=0.87) between diagnostic subgroups, but in response time (χ2=13.6, P=0.001) (Table 1). There were no main effects of sex on amygdala activation in the total sample (left hemisphere: z=−1.72, P=0.09), (right hemisphere: z=−1.17, P=0.24). COPE values from the right and left amygdala showed a strong correlation (rho=0.67, P<0.01).
Genetic associations analyses
Among patients with SCZ, rs237902 was significantly associated with left amygdala activation, whereas the rs237902G allele displayed increased risk of low amygdala activation (P=0.014, Bonferroni corrected P=0.04) (Table 2, Fig. 1). An interaction analysis of rs237902 and sex was performed in this group, but no effect was detected. No significant associations between oxytocin polymorphisms and amygdala activation were found within patients with AD, HC or the total sample (Table 2). Interaction analyses showed a significant effect of having a diagnosis of SCZ compared with HC on the association between rs237902 and amygdala activation (left hemisphere: b=−8.15, P=0.001, Bonferroni corrected P=0.003; right hemisphere: b=−6.77, P=0.01, Bonferroni corrected P=0.03), while this was not the case for patients with AD (left hemisphere: b=−3.42, P=0.17; right hemisphere: b=−5.58, P=0.03, Bonferroni corrected P=0.09). Twelve individuals with outliers were identified after inspection of boxplots (Fig. 2), and repeating the analyses without them strengthened the results (Table 3). No meaningful LD was found across the three SNPs (r 2<0.25).
SNP | Amygdala hemisphere | SCZ n=104 | AD n=100 | HC n=142 | Total N=346 | ||||
---|---|---|---|---|---|---|---|---|---|
Effect size | P | Effect size | P | Effect size | P | Effect size | P | ||
Rs2254298 (G) a | Left | b=0.12 | 0.97 | b=−1.35 | 0.67 | b=1.92 | 0.40 | b=0.72 | 0.66 |
Right | b=−1.11 | 0.77 | b=4.24 | 0.20 | b=0.45 | 0.84 | b=1.22 | 0.47 | |
Rs53576 (G) | Left | b=4.93 | 0.026 | b=0.37 | 0.87 | b=−3.36 | 0.06 | b=0.35 | 0.77 |
Right | b=3.28 | 0.19 | b=0.69 | 0.77 | b=−1.14 | 0.50 | b=0.77 | 0.54 | |
Rs237902 (G) | Left | b=−4.99 | 0.014 b | b=−0.25 | 0.90 | b=2.81 | 0.08 | b=−0.56 | 0.59 |
Right | b=−3.84 | 0.097 | b=−2.66 | 0.18 | b=2.65 | 0.08 | b=−1.00 | 0.36 |
SNP, single nucleotide polymorphism; SCZ, schizophrenia spectrum disorders; AD, affective spectrum disorders; HC, healthy controls.
a Results for the G allele dosage are displayed.
b Remains significant after Bonferroni corrections (3 independent tests).
SNP | Amygdala hemisphere | SCZ n=97 | AD n=99 | HC n=138 | Total N=334 | ||||
---|---|---|---|---|---|---|---|---|---|
Effect size | P | Effect size | P | Effect size | P | Effect size | P | ||
Rs2254298 (G) a | Left | b=0.52 | 0.85 | b=−1.94 | 0.51 | b=1.12 | 0.60 | b=0.37 | 0.80 |
Right | b=0.66 | 0.81 | b=3.60 | 0.24 | b=−0.13 | 0.95 | b=1.29 | 0.38 | |
Rs53576 (G) a | Left | b=2.45 | 0.16 | b=−0.46 | 0.83 | b=−2.60 | 0.12 | b=−0.39 | 0.71 |
Right | b=0.36 | 0.84 | b=−0.24 | 0.91 | b=−0.51 | 0.76 | b=−0.21 | 0.84 | |
Rs237902 (G) a c | Left | b=−4.56 | 0.004 b | b=−0.15 | 0.93 | b=1.93 | 0.19 | b=−0.66 | 0.48 |
Right | b=−3.69 | 0.03 | b=−2.54 | 0.16 | b=2.09 | 0.16 | b=−1.07 | 0.25 |
SNP, single nucleotide polymorphism; SCZ, schizophrenia spectrum disorders; AD, affective spectrum disorders; HC, healthy controls.
a Results for the G allele dosage are displayed.
b Remains significant after Bonferroni corrections (3 independent tests).
c Interaction analyses in the total sample showed a significant effect of having a diagnosis of SCZ compared to healthy controls on the association between rs237902 and amygdala activation (left hemisphere; b=−6.95, P=0.002, Bonferroni corrected P=0.006, right hemisphere; b=−6.05, P=0.009, Bonferroni corrected P=0.03), while this was not the case for patients with AD (left hemisphere; b=−2.46, P=0.26, right hemisphere; b=−4.87, P=0.03, Bonferroni corrected P=0.09).
Discussion
In this study, we show attenuated amygdala activation in response to emotionally negative faces in rs237902G allele carriers among patients with schizophrenia spectrum disorders (SCZ), while not in patients with affective disorders or healthy control participants. Rs237902 is an OXTR synonymous variant located at exon 3 (https://genome.ucsc.edu/cgi-bin/hgTracks?db=hg19&lastVirtModeType=default&lastVirtModeExtraState=&virtModeType=default&virtMode=0&nonVirtPosition=&position=chr3%3A8808934%2D8809434&hgsid=480554839_Izi7M39PmzcwlaTbAUfn3La2J6O0), thus with potential influence on protein expression, conformation and function, and relative potential functional and clinical consequences.Reference Sauna and Kimchi-Sarfaty 37 Montag and colleagues found rs237902G to be associated with more negative symptoms in schizophrenia,Reference Montag, Brockmann, Bayerl, Rujescu, Muller and Gallinat 17 a clinical characteristic found to be correlated with emotion processing and face recognition in a meta-analysis.Reference Ventura, Wood, Jimenez and Hellemann 38 Thus, the main result of the current study is in line with this previous finding. Rs237902 could also serve as a proxy for rs237915 (LD: r 2=0.82), the single SNP Loth and colleagues found to be associated with activation in the ventral striatum in response to threat perception in 1445 healthy adolescents, investigating 23 tag SNPs in the OXTR.Reference Loth, Poline, Thyreau, Jia, Tao and Lourdusamy 39 The nucleus accumbens in the ventral striatum plays an important role in the mesocorticolimbic system and has been hypothesised to interplay with the amygdala via dopaminergic–oxytocinergic mechanisms.Reference Shamay-Tsoory and Abu-Akel 5 Loth and colleagues did not find any associations between OXTR SNPs and the amygdala, similar to our findings in healthy subjects in the present study. Combined, these findings imply that oxytocin plays an important role in emotion processing and provide additional evidence for alterations in the endogenous oxytocin system in patients with SCZ.Reference Uhrig, Hirth, Broccoli, von Wilmsdorff, Bauer and Sommer 40
The influence of rs53576 and rs2254298, two intronic polymorphisms, has been extensively studied across social features, psychiatric symptoms and ethnic groups; however, the results are mixed in spite of large meta-analyses.Reference Bakermans-Kranenburg and van Ijzendoorn 19 , Reference LoParo and Waldman 41 , Reference Li, Zhao, Li, Broster, Zhou and Yang 42 Tost and colleagues showed associations between rs53576 and rs2254298 and amygdala activation in separate studies in 228 healthy individuals using the same fMRI paradigm as the present study,Reference Tost, Kolachana, Hakimi, Lemaitre, Verchinski and Mattay 43 , Reference Tost, Kolachana, Verchinski, Bilek, Goldman and Mattay 44 findings neither we nor Loth and colleagues were able to replicate. While Loth and colleagues suggest that this may be because of their use of a different emotion processing paradigm and inclusion of adolescents instead of adults, our lack of replication might be because of a smaller sample of healthy individuals and the use of a clinical sample that may not show associations to the same OXTR polymorphisms as healthy individuals.
Another possible explanation for ambiguous genetic associations regarding oxytocin pathway genes is potential interaction effects with sex.Reference Feldman, Monakhov, Pratt and Ebstein 45 The expression of central oxytocin is steroid-dependent,Reference Insel 46 the central regulation differs between men and women,Reference Dumais and Veenema 47 and men and women respond differently to exogenous oxytocin.Reference Shamay-Tsoory and Abu-Akel 5 Additionally, several brain imaging studies investigating oxytocin pathway genes in healthy individuals have detected sex-specific effects.Reference Tost, Kolachana, Verchinski, Bilek, Goldman and Mattay 44 , Reference Wang, Qin, Liu, Zhou, Wang and Zhang 48 , Reference Wang, Qin, Liu, Wang, Zhang and Jiang 49 However, we were not able to detect any interaction effects between sex and rs237902 in patients with SCZ, and all genetic association analyses included sex as a covariate. Hence, our genetic association results presented are not because of sex differences, and the disorder-specific association with rs237902 did not differ between men and women with SCZ.
As no oxytocin-related polymorphisms have been detected as risk variants for SCZ in large genome-wide association studies (http://www.broadinstitute.org/mpg/ricopili/), we chose a hypothesis-based approach for SNP selection. Rs53576, rs2254298 and rs237902 are the only SNPs previously found associated with disease characteristics in Caucasian schizophrenia samples; Reference Bartholomeusz, Ganella, Labuschagne, Bousman and Pantelis 10 , Reference Montag, Brockmann, Lehmann, Muller, Rujescu and Gallinat 16 – Reference Haram, Tesli, Bettella, Djurovic, Andreassen and Melle 18 hence, we were interested in these polymorphisms in particular. However, only a few studies have focused on the oxytocin pathway genes in patients with schizophrenia, and the possibility of publication bias must be taken into account. The clinical consequences and functions of these variants are also as yet unknown, and their influence on brain function should be interpreted with caution.
To conclude, we found a disorder-specific contribution of rs237902 to amygdala activation during emotional face processing in patients with SCZ. This is, to our knowledge, the first fMRI study investigating oxytocin pathway polymorphisms in this patient group. Our findings suggest that the endogenous oxytocin system contributes to biological underpinnings of emotion processing, that this contribution differs between patients with SCZ and HC and that exonic regions in OXTR are of interest for future research. Future research should investigate additional oxytocin pathway polymorphisms as well as involving other brain circuits, to better understand the biological mechanisms underlying emotion processing. We encourage replication of our findings in further imaging studies and meta-analyses of oxytocin pathway genes in this patient group.
Funding
This work was supported by the Kristian Gerhard Jebsen Foundation, the Research Council of Norway (grant numbers 181831, 147787/320, 167153/V50, 204966/F20) and the Regional Health Authority for South-Eastern Norway Health Authority (grant numbers 2014–102, 2014–097, 2010–074, 2006–258). The funders had no role in study design, data collection, analysis and interpretation, writing of the report or the decision to submit the paper for publication.
Acknowledgements
We thank the patients and healthy controls for participating in the study, and the TOP study group members for contributing with data collection. We also thank Eivind Bakken, Monica Aas and Christian Melbø-Jørgensen for administrative and statistical assistance.
eLetters
No eLetters have been published for this article.