Demographic characteristics and genotypic data
The demographic characteristics of the subjects are shown in Table 1. One-way ANOVA indicated that there was a significant age difference (F (1,781) =32.83, p<0.01), but not in education between patients and healthy controls. In addition, there was a significant sex difference in age in the healthy control group (F (1,417) =6.59, p<0.05), and significant sex difference in educations in both healthy control and patients (F (1,413) =4.17, p<0.05; F (1,365) =5.35, p<0.05). There was no significant difference in the age of onset between male patients and female patients.
The COMT Val158Met genotype distribution was consistent with Hardy–Weinberg equilibrium in healthy controls (c2=2.83, df=1, p=0.24) and patients (c2=0.03, df=1, p=0.98). The distribution of COMT Val158Met genotype and allele were summarized in Table 1.There was no significant difference in COMT Val158Met genotype (c2=2.27, df=2, p>0.05) and allele distribution (c2=0.52, df=1, p>0.05) between the healthy controls and patients. However, there were significant sex differences in the COMT Val158Met genotype (c2=10.3, df=2, p<0.01) and allele frequency in patients (c2=4.16, df=1, p<0.05), showing that male patients had a higher proportion of Met alleles. However, there was no such a sex difference in healthy controls (c2=0.82, df=2, p>0.05; c2=0.03, df=1, p>0.05, respectively).
Interaction of diagnosis, sex and COMT Val158Met genotype on the cognitive function
As shown in Table 2, three-way MANCOVA indicated the interaction analysis results of the diagnosis, sex, and genotype. The main effects of diagnosis were significant in all cognitive indexes (all p or FDR corrected p<0.05; Cohen’s f=0.54, 0.23, 0.49, 0.30, 0.70,0.47 respectively in immediate memory, visuospatial/constructional, language, attention, delayed memory, and total score). There was significant sex effect in delayed memory (F(1,631)=5.18, p<0.05, FDR corrected p=0.02, Cohen’s f=0.18). There was significant COMT genotype effect on language (F(1,631)=4.92, p<0.05, FDR corrected p=0.03, Cohen’s f=0.18). There was significant diagnosis×sex interaction in delayed memory (F(1,631)=5.58, p<0.05, FDR corrected p=0.09, Cohen’s f=0.19), and diagnosis×COMT interaction in language (F(1,631)=4.77, p<0.05, FDR corrected p=0.10, Cohen’s f=0.17), and sex×COMT interaction (F(1,631)=4.54, p<0.05, FDR corrected p=0.09, Cohen’s f=0.19) in immediate memory (F(1,631)=4.54, p<0.05, FDR corrected p=0.11, Cohen’s f=0.17), but those differences did not pass the correction.
Interaction of sex and COMT Val158Met genotype on the clinical characteristics of schizophrenia patients
The interaction of sex and COMT Val158Met polymorphism on the clinical characteristics of patients with schizophrenia were summarized in Table 3. After adjusting for age and education, two-way MANCOVA showed that there was a significant sex × genotype interaction effect on the positive subscale score (F(1,238)=7.29, p<0.05, FDR corrected p=0.04, Cohen’s f=0.35), but not on the negative subscale, general psychopathology subscale or PANSS total scores (all p>0.05). In addition, there were no main effects of sex or genotype on the PANSS subscales and total scores (all p>0.05). One-way MANCOVA showed that in male patients, there was a significant genotype effect on the positive subscale (F(1,173)=5.42, p<0.05, FDR corrected p=0.03, Cohen’s f=0.50). However, there was no genotype effect in female patients (p>0.05). Further post hoc analysis (Fig 1A) indicated that male patients with Met homozygote and heterozygote had a significantly lower score than female patients with Met homozygote and heterozygote genotypes (p<0.01 ). Male patients with Val/Val genotype had more positive symptoms than Met carriers (Val/Met vs. Met carrier: 13.3±4.8 vs. 11.6±4.6, p<0.05), while in female patients, Met carriers exhibited more positive symptoms slightly (12.3±5.4 vs. 14.7±7.3, p=0.07).
Interaction of sex and COMT Val158Met genotype on cognitive performance in healthy control and schizophrenia patients
In patients, two-way MANCOVA with age, age of onset, and education as covariates indicated that there were significant sex effect on immediate memory (F(1,235)=3.98, p<0.05, FDR corrected p=0.06, Cohen’s f=0.26), delayed memory (F(1,235)=6.91, p<0.01, FDR corrected p=0.02, Cohen’s f=0.34), and total score (F(1,235)=4.33, p<0.05, FDR corrected p=0.06, Cohen’s f=0.27). Moreover, there was a significant genotype effect on language (F(1,235)=6.57, p<0.05, FDR corrected p=0.01, Cohen’s f=0.33). There was a significant sex × genotype interaction on immediate memory (F(1,235)=4.44, p<0.05, FDR corrected p=0.07, Cohen’s f=0.33), attention (F(1,235)=4.92, p<0.05, FDR corrected p=0.05, Cohen’s f=0.29), and total score (F(1,235)=4.21, p<0.05, FDR corrected p=0.04, Cohen’s f=0.27).
One-way MANCOVA showed that in male patients, there was significant genotype effect on all the cognitive indexes: immediate memory (F(1,175)= 12.9, p<0.01, FDR corrected p=0.002, Cohen’s f=0.54), visuospatial/constructional (F(1,175)= 5.9, p<0.05, FDR corrected p=0.03, Cohen’s f=0.37), language (F(1,175)= 4.4, p<0.05, FDR corrected p=0.05, Cohen’s f=0.32), attention (F(1,175)=4.6, p<0.01, FDR corrected p= 0.05, Cohen’s f=0.32), delayed memory (F(1,175)=4.3, p<0.01, FDR corrected p= 0.05, Cohen’s f=0.31), and total score (F(1,175)=9.0, p<0.01, FDR corrected p=0.01, Cohen’s f=0.45). However, in female patients, there were no significant differences between genotypes in any of the cognitive indexes. Further post hoc analysis showed that male patients with Val homozygotes had the lowest score in nearly all the cognitive indexes (Fig1 B,C,D,E,G,F).
In the healthy controls, we did not find any significant sex effects, genotype effects, or sex × genotype interactions on cognitive performance (all p>0.05). Further post hoc analysis showed that in either male or female healthy controls, there was no significant genotype effect on any cognitive indexes (all p>0.05).