{'fraction_non_alphanumeric': 0.03714875377043975, 'fraction_numerical': 0.03992697253532307, 'mean_word_length': 5.286926147704591, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Lack of Association between Oxytocin Receptor (OXTR) Gene Polymorphisms and Alexithymia: Evidence from Patients with Obsessive-Compulsive Disorder

Oxytocin receptor gene single nucleotide polymorphisms have been associated with structural and functional alterations in brain regions, which involve social-emotional processing. Therefore, oxytocin receptor gene polymorphisms may contribute to individual differences in alexithymia, which is considered to be a dysfunction of emotional processing. The aim of this study was to evaluate the association between oxytocin receptor gene single nucleotide polymorphisms or haplotypes and alexithymia in patients with obsessive-compulsive disorder. We recruited 355 patients with obsessive-compulsive disorder (234 men, 121 women). Alexithymia was measured by using the Toronto Alexithymia Scale. We performed single-marker and haplotype association analyses with eight single nucleotide polymorphisms (rs237885, rs237887, rs2268490, rs4686301, rs2254298, rs13316193, rs53576, and rs2268498) in the oxytocin receptor gene. There were no significant associations between any of the eight single nucleotide polymorphism of the oxytocin receptor gene and alexithymia. In addition, a six-locus haplotype block (rs237885-rs237887-rs2268490-rs4686301-rs2254298-rs13316193) was not significantly associated with alexithymia. These findings suggest that genetic variations in the oxytocin receptor gene may not explain a significant part of alexithymia in patients with obsessive-compulsive disorder.


Introduction
Alexithymia is characterized by difficulties in recognizing and effectively expressing one's own feelings [1]. It is difficult for highly alexithymic individuals to distinguish or appreciate the emotions of others [2], and this trait results in an unempathic and ineffective emotional response [3]. These individuals feel less distressed on seeing others in pain, and such blunted emotions can be translated into impaired prosociality [4]. Effective emotional processing, Interview for the DSM-IV (SCID) [27] assessed by a trained psychiatrist (S.J. KIM). This group was composed of patients in different stages of the illness and with different degrees of severity. All patients were taking psychotropic medications (mainly selective serotonin reuptake inhibitors and/or low-dose benzodiazepines). All participants were Korean and gave written informed consent prior to the beginning of this study. The study protocol was approved by the Institutional Review Board of Severance hospital, Yonsei University Health System.

Assessment of alexithymia
The degree of alexithymia was measured with the Toronto Alexithymia Scale , which is a 20-item self-report scale with a five-point Likert-type scale [28]. The TAS-20 comprises three subdimensions: 1) difficulty in identifying feelings (DIF, seven items), 2) difficulty in describing feelings (DDF, five items) and 3) externally oriented thinking (EOT, eight items). The TAS-20 was previously proven to be valid and reliable [28] and was validated for the Korean population [29]. All participants completed the Korean version of the 20-item Toronto Alexithymia Scale [30].

Measures of clinical symptoms
The clinical symptoms and the severity of OCD symptoms were evaluated by using the Yale-Brown obsessive-compulsive scale (Y-BOCS) [31]. The Y-BOCS is a reliable and valid 10-item scale administered by a clinician and used to assess the severity of obsessions and compulsions. Levels of depressive symptoms were assessed by using the Montgomery-Åsperg Depression Rating Scale (MADRS) [32], which is a well-known scale widely used by trained psychiatrists.

SNP selection and Genotyping
First, we selected OXTR SNPs which showed significant association with empathy, prosocial behavior, sociability, emotionality, social cognition, and social deficit by using a PubMed search of previous studies [19,[33][34][35][36][37][38][39][40][41][42]. Nine SNPs were selected considering reported minor allele frequencies (MAF) greater than 0.1 from the 1000 Genomes Project database, JPT sample (June 2010 release). Among them, rs2268491 was excluded because it was in complete linkage disequilibrium (LD; D 0 = 1, r 2 = 0.97) with rs2254298 in our sample. Therefore, eight SNPs with an r 2 threshold < 0.8 in 'pair-wise tagging only' mode using the 'Tagger' program in Haploview [43] were included in the final statistical analysis (Table 1; Fig 1). Peripheral blood samples were obtained from each subject, and genomic DNA was extracted from the leukocytes. The SNaPshot assay was performed according to the manufacturer's instructions (ABI PRISM SNaPShot Multiplex kit, Foster City, CA, USA).

Statistical analyses
Hardy-Weinberg equilibrium values of the SNPs were tested using χ 2 tests. Generalized estimating equations (GEE) were used to analyze the relationship between total or three subdimension scores of TAS-20 and SNPs, adjusting for age, total Y-BOCS, and MADRS scores. An additive model was applied by coding genotypes as 0, 1, or 2, depending on carrier status of the minor allele.

Results
Demographic and clinical characteristics of the subjects are presented in Table 2. The mean age of subjects was 30.48 ± 11.13 years old. The mean age of onset of obsessive-compulsive symptoms was 18.82 ± 9.21. The mean scores of Y-BOCS, and MADRS were 22.16 ± 8.70, and 19.07 ± 8.99, respectively. Because the patients' total scores and some subdimension scores on the TAS-20 were significantly correlated with the total Y-BOCS and MADRS scores, and age, although the size of the correlation between the TAS-20 scores (total and subdimension scores) and the Y-BOCS or age were negligible (absolute r = 0.13-0.19, data not present), we used these factors as covariates in subsequent analyses to control for their potential effects on alexithymia.
We used GEE approaches to examine the effect of the OXTR gene polymorphisms on the total or subdimension scores of the TAS-20 for each of eight SNPs. There were no significant associations found between any of the OXTR SNPs and alexithymia. In addition, a six-locus haplotype block (rs237885-rs237887-rs2268490-rs4686301-rs2254298-rs13316193) also was not significantly associated with the severity of alexithymia. The results are shown in Tables 3,  4 and 5.

Discussion
Initially, we hypothesized that OXTR variants might account for the individual differences of alexithymic traits seen in OCD patients. However, contrary to our initial hypothesis, there were no differences of alexithymic traits according to alleles or genotypes of eight SNPs within OXTR. Haplotype analyses also did not show any association between various haplotypes of OXTR and alexithymia. Alexithymic individuals show impaired emotional experience, deficits of emotional interpretation and emotional face recognition, and lack of empathy [44] and are associated with various psychiatric conditions including autistic disorders and psychopathy [45,46]. Although there have been no reports which directly investigated the influence of OXTR genetic variants on alexithymic traits, there are a number of studies of the association between OXTR polymorphisms and those various alexithymia-related traits (e.g., empathy, emotional facial recognition, and so forth) and psychiatric disorders (e.g., autism spectrum disorders, psychopathy, and so forth).
Recently, Laursen et al. found that subjects with the CC genotype at OXTR rs2268498 and AA genotype at OXTR rs53576 showed higher empathic accuracy [47]. Meanwhile, Uzefovsky et al. reported that the A allele of OXTR rs53576 predicted lower emotional empathy [48]. Melcher et al. found that T allele carriers of OXTR rs2268498 showed more accurate facial emotional recognition skill [49]. In addition, subjects with the A allele of OXTR rs2254298 showed deficient deactivation of the dorsal anterior cingulate gyrus during an emotional face matching task [50]. The SNP rs13316193 C allele of the OXTR gene has been associated with empathy [33], whereas the T allele has been linked to decreased expression of oxytocin receptors in the brain, depressive mood, and greater risk for autism spectrum disorder [19]. One recent metaanalysis study found significant association between autism spectrum disorder and the rs7632287, rs237887, rs2268491, and rs2254298 SNPs of OXTR [51]. Psychopathy has a higher incidence of alexithymic traits, and both psychopathy and alexithymia have common core features of a lack of empathy, insight, and introspection [52]. Several neuroimaging studies suggested that the genetic variants of OXTR modulate the activities of limbic circuits including the amygdala, the hypothalamus, and the cingulate gyrus [53], which are also associated with alexithymia [23].
However, not all of the studies consistently reported a positive association between OXTR variants and alexithymia-related traits. Nyffeler et al. did not find any associations between single SNPs of OXTR (rs2301261, rs53576, rs2144298, or rs2268494) and autism [54]. Tansey  [55]. Also, a meta-analyses study consisting of a large number of participants (N 17000 for rs53576 and N for rs2254298) failed to support the impact of these two OXTR gene variants on five domains of human functioning including biology, personality, social behavior, psychopathology, and autism [56]. However, the resulting lack of association in our study should be interpreted with caution because there are several factors, which are required to be considered. First, the ethnic and cultural backgrounds might influence the results of our study. In fact, several inconclusive results about autism have been reported from different ethnic and cultural backgrounds across studies. In Japanese [57] and Chinese [58] populations, the rs2254398 A allele of OXTR has been reported to be linked to autism. However, in Caucasian autism trios, the rs2254395 G allele of OXTR was overtransmitted to probands with autistic disorders [59]. In addition, as mentioned earlier, one study of a Caucasian sample found no influence of various SNPs of the OXTR gene on autism [55], although they did not include the rs2254398. We did not control the potential effects of undetected population stratification. Notwithstanding, the Korean population is considered to be much more ethnically and culturally homogenous than other populations, due to its distinct language and culture. Although the Koreans are assumed to be relatively free from stratification [60], there is still potential risk of biased results from undetected population stratification. Second, the participants in this study were all OCD patients, and the disease status of the sample might influence the results. For example, in schizophrenic patients, the A allele carriers of the rs2254298 SNP of OXTR had higher empathic concern than non-A allele carriers, whereas in healthy controls, this difference was not found [35]. Therefore, although patients with OCD tend to be highly alexithymic with a relatively large variability, which can increase statistical power, the results cannot be generalized to a non-clinical sample or to other psychiatric conditions. Third, we did not consider gene and environment interaction effects. However, there has been evidence suggesting an OXTR gene by environment interaction on various human traits [56] (e.g., emotional dysregulation and attachment style) or psychiatric diseases [61] (e.g., depression). Lastly, while our sample size was relatively larger than those in previous studies examining associations between OXTR SNPs and alexithymia-related conditions, our power is still limited. When we simulated the power of our study (5,000 trial runs) using the JPT + CHB panel of a 1,000 genome database and phenotypic variance-covariance matrix estimated from this study, the power of the sample size in this study was only 0.342, and was sufficient for detecting only an effect size larger than 0.9, as its regards to an association between rs237885 (the lowest GEE p-value in this study) and three subdimension scores of TAS-20. However, the effect size of rs237885 in this study was only 0.454 (standardized regression coefficient, β), suggesting the possibility of type II errors. Therefore, to confirm our reported findings, a study with a much larger number of samples with stronger statistical power is needed. Therefore, it would be necessary to evaluate the effects of OXTR variants on alexithymia with particular consideration towards those confounding factors and limitations.
Supporting Information S1 Table. Dataset of all subjects used in this study.