Next Article in Journal
Consumption of Coffee and Tea Is Associated with Macular Retinal Nerve Fiber Layer Thickness: Results from the UK Biobank
Previous Article in Journal
Pre-Therapeutic Sarcopenia among Cancer Patients: An Up-to-Date Meta-Analysis of Prevalence and Predictive Value during Cancer Treatment
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Nutrients
Volume 15
Issue 5
10.3390/nu15051195
nutrients-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Association between Arachidonic Acid and the Risk of Schizophrenia: A Cross-National Study and Mendelian Randomization Analysis

by
Yan Gao
1,†,
Xiaowen Hu
1,†,
Dandan Wang
1,
Jie Jiang
1,
Minghui Li
1,
Ying Qing
1,
Xuhan Yang
1,
Juan Zhang
1,
Yue Zhang
2,* and
Chunling Wan
1,*
1
Bio-X Institutes, Key Laboratory for the Genetics of Developmental and Neuropsychiatric Disorders (Ministry of Education), Shanghai Jiao Tong University, Shanghai 200030, China
2
Department of Bioinformatics and Biostatistics, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Nutrients 2023, 15(5), 1195; https://doi.org/10.3390/nu15051195
Submission received: 10 January 2023 / Revised: 20 February 2023 / Accepted: 24 February 2023 / Published: 27 February 2023
(This article belongs to the Section Nutrition and Public Health)
Browse Figures
Versions Notes

Abstract

:
Polyunsaturated fatty acids (PUFAs), especially long-chain PUFAs (LCPUFAs), are crucial for both the structural and functional integrity of cells. PUFAs have been reported to be insufficient in schizophrenia, and the resulting cell membrane impairments have been hypothesized as an etiological mechanism. However, the impact of PUFA deficiencies on the onset of schizophrenia remain uncertain. We investigated the associations between PUFAs consumption and schizophrenia incidence rates through correlational analyses and conducted Mendelian randomization analyses to reveal the causal effects. Using dietary PUFA consumption and national schizophrenia incidence rates in 24 countries, we found that incidence rates of schizophrenia were inversely correlated with arachidonic acid (AA) and ω-6 LCPUFA consumption (rAA = −0.577, p < 0.01; rω-6 LCPUFA = −0.626, p < 0.001). Moreover, Mendelian randomization analyses revealed that genetically predicted AA and gamma-linolenic acid (GLA) were protective factors against schizophrenia (ORAA = 0.986, ORGLA = 0.148). In addition, no significant relationships were observed between schizophrenia and docosahexaenoic acid (DHA) or other ω-3 PUFAs. These findings show that the deficiencies of ω-6 LCPUFAs, especially AA, are associated with schizophrenia risk, which sheds novel insight into the etiology of schizophrenia and a promising diet supplementation for the prevention and treatment of schizophrenia.

1. Introduction

Schizophrenia is a complex psychotic disorder characterized by a bundle of neuropsychiatric manifestations, including cognitive impairment, apathy and social withdrawal [1]. The onset of schizophrenia typically occurs in young adulthood and can result in life-long impairment [1]. The median annual incidence rate of schizophrenia is 15.2 per 100,000 with substantial variations across countries and cultural groups [2]. Both genetic and environmental factors that impair neurodevelopment are considered risk factors for schizophrenia, including genetic predisposition provided by mutant neurodevelopmental genes and nutritional factors in prenatal life [1,3]. The indication of maternal malnutrition as a causative factor in schizophrenia was evidenced by higher incidence rates of schizophrenia among offspring during famines that occurred in the Netherlands during 1944–1945 [4] and in China during 1959–1961 [5].
An especially important nutrient group for the course of brain development is polyunsaturated fatty acids (PUFAs). PUFAs are essential fatty acids that are mainly obtained from the diet. The dietary intake of PUFAs varies widely due to geographic and cultural differences across countries [6]. Linoleic acid (LA, C18:2 ω-6) and alpha-linolenic acid (ALA, C18:3 ω-3) are the most abundant dietary PUFAs [6], and can be endogenously converted to long-chain PUFAs (LCPUFAs) to a limited degree. Arachidonic acids (AAs, C20:4 ω-6) and docosahexaenoic acid (DHA, C22:6 ω-3) are the major LCPUFAs in diets in most counties. LCPUFAs are an important nutrient group for brain development. In infancy and early childhood, both ω-6 and ω-3 LCPUFAs, especially AA and DHA, are critical for supporting brain growth and maturation [7]. Several studies on children have revealed improved performance in cognition and motor skills after LCPUFA supplementation [8,9,10], and LCPUFA deficiency can lead to psychotic symptoms in those with neurodevelopmental disorders [11]. Moreover, meta-analyses have shown that PUFA concentrations are insufficient in schizophrenia [12,13], and a similar tendency was found in individuals at an ultrahigh risk of developing psychosis [14].
Our study aims to investigate the associations between PUFAs and schizophrenia risk. Correlations between PUFA dietary consumption and national incidence rates of schizophrenia are evaluated from an epidemiological perspective, and Mendelian randomization (MR) analysis is conducted to address the causative relationships and evade residual confounding. We specifically hypothesize that lower LCPUFAs would be associated with an increased schizophrenia risk.

2. Materials and Methods

Our study consists of cross-national correlation analysis and MR analysis. A schematic diagram of our study design is shown as below (Figure 1).

2.1. Cross-National Correlation between PUFA Consumption and Schizophrenia Incidence Rates

2.1.1. Data Acquisition on the Incidence Rates of Schizophrenia

The incidence rates of schizophrenia were obtained from published epidemiological studies. We included studies published between January 1975 and December 2021, which were identified through the search terms (incidence (Title/Abstract) AND schizophrenia (Title/Abstract)) applied in MEDLINE, PsychINFO and EMBASE, as well as the literature referred to in systematic reviews [2,15]. As the study design had a large impact on estimating the incidence rate of schizophrenia [16], the criteria for including studies in this analysis are listed as follows: Studies were included if they reported primary data on the incidence of schizophrenia according to the diagnostic criteria of the International Classification of Diseases (ICD) for schizophrenia. Studies within special populations only (e.g., twins, young adults) or within birth cohorts were excluded. When multiple studies on schizophrenia incidence in the same country were available, the study that reported the investigated population at risk as the total population or the one that sampled from the most recent period was preferred. Incidence was recorded as a rate per 100,000 population. Ethical approval was obtained in all original studies.

2.1.2. Data Acquisition on the Dietary Consumption of PUFAs

Dietary consumption data of PUFAs were obtained from the Food and Agriculture Organization of the United Nations and the Australian Food Composition Database (previously called NUTTAB). The method of calculations followed the methodological process of establishing a global database of nutrients [17]. Briefly, food consumption data were obtained from the Food and Agriculture Organization, which annually compiles Supply and Utilization Accounts and provides internally consistent information about up to 394 food and agricultural commodities across nations. We extracted the available information about food supply quantity per capita (g/day) for countries with available incidence data (Table S1). To provide a more accurate estimate of the consumption of individual PUFAs, we introduced a wastage index to adjust for retail and consumption loss across countries [18] as well as the refuse factor to adjust for inedible parts [17]. Then, we matched these individual food items with the food items in the Australian Food Composition Database, a database that provides complete information on known dietary sources for PUFA intake. Finally, we added up the contributions of individual food items to each PUFA. Therefore, individual consumption of each PUFA within each country per year was calculated in aggregate. PUFA groups (total ω-3 PUFA, total ω-6 PUFA, total ω-3 LCPUFA, and total ω-6 LCPUFA) consumption was indirectly measured by summing up.

2.1.3. Correlation Adjustment

A twenty-year lag interval was introduced to evaluate the causative relationships between PUFA intake in early life and schizophrenia risks in consideration of an adequate time span for the onset of psychiatric symptoms. Since PUFA consumption is correlated with economic development, we used the GDP per capita to account for the socioeconomic differences, paralleling previous studies [19,20,21,22].

2.1.4. Statistical Analysis

In the cross-national association study, seventeen individual PUFAs and four summary PUFA groups were included for analyses. The analysis involved simple Pearson’s product moment correlation and partial correlation analyses adjusting for economic factors. These statistical procedures were conducted in SPSS 20.

2.2. Mendelian Randomization Analysis for PUFAs and Schizophrenia

To investigate the potential causal associations between PUFAs and schizophrenia, we applied a two-sample MR approach by using summary data from genome-wide association studies (GWASs). Ethical approval was obtained in all original studies.

2.2.1. PUFA Exposure Data Acquisition

We obtained published association results for genetic instruments of ω-6 [23] and ω-3 PUFAs [24] from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, which are the most common datasets of ω-6 and ω-3 PUFA levels used for MR analysis [25,26,27,28]. The summary statistics were based on a meta-analysis of GWAS of individuals of European ancestry, which are 8631 European individuals for ω-6 PUFAs and 8866 European individuals for ω-3 PUFAs (Table 1). Participants for the meta-analysis were drawn from 5 cohort studies (Table S2), including the Atherosclerosis Risk in Communities Study (ARIC, n = 3269), the Cardiovascular Health Study (CHS, n = 2404), the Coronary Artery Risk Development in Young Adults Study (CARDIA, n = 1507), the Invecchiare in Chianti Study (InCHIANTI, n = 1075), and the Multi-Ethnic Study of Atherosclerosis (MESA, n = 707). In these studies, fatty acids were measured using gas chromatography technique. Levels of all individual fatty acids, including LA (average range, 19.96~24.78%), gamma-linolenic acid ((GLA) average range, 0.11~0.12%), dihomo-γ-LA ((DGLA) average range, 3.13~3.33%), AA (average range, 8.00~12.10%), ALA (average range, 0.14~0.44%), DHA (average range, 2.29~3.66%), eicosapentaenoic acid ((EPA) average range, 0.56~0.88%) and docosapentaenoic acid ((DPA) average range, 0.83~0.95%), were expressed as a percentage of total fatty acid (Table 1 and Table S2).

2.2.2. Schizophrenia Outcome Data Acquisition

We obtained the summary statistics data on genetic associations for schizophrenia from Psychiatry Genomics Consortium (PGC) [29]. To mitigate population stratification bias, we extracted genetic data from the schizophrenia GWAS meta-analysis with the largest proportion of European population that included 36,989 cases and 113,075 controls, of which 34,241 sporadic cases and 45,604 ancestry-matched controls were used in the subsequent MR analysis (Table 1). To avoid bias due to overlapping sample sets in MR analyses, we checked the 49 individual participating studies of schizophrenia GWAS and found no overlap with the fatty acid GWASs. All genetic statistical data were downloaded from the PGC database (http://pgc.unc.edu, accessed on 28 January 2021).

2.2.3. Univariable MR Analysis

In this MR analysis, four ω-6 PUFAs (LA, GLA, DGLA, and AA) and four ω-3 PUFAs (ALA, EPA, DHA, and DPA) were included. All included SNPs were significantly associated with PUFA levels at the genome-wide significance level (p < 5 × 10−8). Bias in effect estimates of MR analysis can be induced by correlation between SNPs. In order to minimize this bias, our genetic instruments used in each PUFA were limited to independent SNPs without linkage disequilibrium (R2 < 0.1). We calculated F statistics to assess the strength of each instrument. If F > 10, there is sufficient strength to avoid weak instrument bias in the MR analysis. The F statistics were computed by the formula F = (N-k-1) × R2/(1 − R2), in which R2 = 2 × MAF × (1 − MAF) × β2. N and k refer to the sample size and numbers of the instrument separately. All SNPs used in univariable MR analysis satisfied the MR criteria (F-statistic > 10) [30]. In addition, the instrumental SNPs are biologically relevant to PUFAs. The instruments of single PUFAs with adjustments for their precursors, including GLA, DGLA, and AA, were also obtained in ω-6 PUFAs GWAS from the CHARGE consortium [23].
We harmonized the summary-level data to ensure that the allele for each SNP corresponded between each PUFA and schizophrenia. In the main analysis, we used the fixed-effects inverse-variance-weighted (IVW) method to assess the causal effects between exposure and outcome. The odds ratio (OR) of schizophrenia was calculated per 1 SD increment in genetically predicted plasma fatty acid levels. We applied p < 0.05 as statistical significance criteria to classify SNPs as potentially influential. Two-sample MR analysis was implemented by using the “MendelianRandomization” package in R 4.0.

2.2.4. Multivariable MR Analysis

We also conducted two-sample multivariable MR analysis to simultaneously estimate the direct effect of ω-3 and ω-6 PUFAs on schizophrenia in consideration of their nexus in metabolism. Considering the tight relationships between ω-3 and ω-6 PUFAs of sharing the same metabolic enzymes, multivariable MR analysis focused on estimating LA (18:2 ω-6) and ALA (18:3 ω-3) effects in comparison to each other, as well as AA (C20:4 ω-6) and EPA (C20:5 ω-3) effects in comparison to each other. Links regarding other PUFAs in equivalent positions in the synthesis pathway were not evaluated due to the lack of GWAS summary data from the same GAWS study. A total of eight SNPs related to ω-3 and ω-6 PUFAs were extracted from the CHARGE consortium. The final instrument set contained 5 SNPs due to linkage disequilibrium clumping. The schizophrenia outcome data were previously described. Multivariable MR was conducted using the R package “Mendelian Randomization”. The F-statistic was calculated to assess instrument strength, which was performed using the “MVMR” package.

2.2.5. Pleiotropic Associations

We also considered potential biological and socioeconomic confounders that may affect the associations between PUFAs and schizophrenia risk. We explored the effect of some potential confounders such as C-reactive protein (CRP), educational attainment, and vitamin D levels [31,32,33]. A previous study found that instrumental variables of PUFAs were not causally related to education attainment and vitamin D [34], which illustrate that it was unlikely that PUFA instruments would affect the schizophrenia risk through these traits. Additionally, to address potential confounder CRP, we conducted linkage disequilibrium (LD) analyses [35] to investigate whether the inflammatory cytokine would affect the association between AA and schizophrenia risk. Genetic associations of CRP were obtained from a recent GWAS based on a large sample of >200,000 European individuals, which identified 58 genome-wide significant loci [36]. Four significant SNPs (rs10832027 and rs1582763 in chromosome 11, rs10521222 and rs1558902 in chromosome 16) in CRP GWAS are in the linkage distance to our instruments rs174547 and rs16966952 separately, which were assessed in the subsequent LD analyses. Furthermore, we searched PhenoScanner [37] for other potential pleiotropic associations of the instrumental variables with risk factors for schizophrenia.

3. Results

3.1. Cross-National Correlation between Dietary PUFA Intake and Schizophrenia Incidence Rates

Initially, 41 potentially related papers were collected based on the strict inclusion criteria described above to obtain national incidence rates of schizophrenia [38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78] (Table S3). Eight countries had more than one study reporting the incidence rate of schizophrenia. With further exclusions on studies decades ago or case findings based on the whole population rather than who was at an age known to be at risk, a unique incidence rate of each country was included for subsequent analysis (Table 2). The incidence rates of schizophrenia varied from 4.1 per 100,000 population in the United Kingdom to 58.5 per 100,000 population in India.
Eating habits vary across different countries and cultures, leading to a global difference in the major dietary PUFA categories and a wide range of PUFA intake. The NUTTAB included 17 dietary PUFAs, and the scales of the dietary PUFAs differed by two orders of magnitude; on average, LA accounted for 88.3% of the total PUFA intake, whereas adrenic acid (AdrA, C22:4 ω-6) accounted for only 0.08%. AA and DHA were the major individual LCPUFAs among most countries (Figure 2). PUFA compositions of daily intake across countries and a potential confounder, GDP per capita, are displayed in Table 2.
Correlation analyses showed that a higher consumption of AA was correlated with a lower schizophrenia incidence rate (rAA = −0.577, p < 0.01) (Figure 3A). Moreover, its precursor (eicosadienoic acid, C20:2 ω-6) and product (AdrA) were related to schizophrenia incidence rates in similar relationships. Dietary consumption of other individual PUFAs, including DHA, did not show statistically significant correlations (Figure 3B). Given the potential socioeconomic effects on the consumption of food, we introduced the GDP per capita across different countries to adjust economic levels by partial correlation analysis. As a result, AA consumption remained inversely correlated with the incidence rates of schizophrenia.
We further investigated schizophrenia incidence rates against different PUFA families. Similarly, the results showed that higher ω-6 LCPUFA consumption was correlated with a lower incidence rate of schizophrenia (rω-6 LCPUFA = −0.626, p < 0.001) (Figure 3C), whereas no significant relationship was observed between ω-3 LCPUFA consumption and schizophrenia incidence rates (Figure 3D). In addition to LCPUFAs, an association analysis between total ω-6 PUFA consumption and schizophrenia was conducted and found no significant relationship. Total ω-3 PUFA consumption showed a weak negative association with schizophrenia incidence rates, which was no longer significant after adjustments for economic factors. In addition, there were no correlational relationships between the incidence rates of bipolar disorder or depressive disorder and dietary PUFAs within the limited available dataset [79,80,81,82] (Table S4).
Taken together, these results suggest that the high intake of ω-6 LCPUFAs, especially AA, would decrease the risk of schizophrenia onset.

3.2. Mendelian Randomization Analysis of PUFA and Schizophrenia

To further address the causal relationship between PUFAs and schizophrenia, we performed MR analyses. The variants used as genetic instrumental variables of the eight different PUFAs are listed in Table 3. These SNPs regarding each PUFA were independent without linkage disequilibrium and were not directly associated with schizophrenia at the genome-wide significance threshold. For all instruments related to each PUFA, F-statistics > 10, suggesting that the analyses were unlikely to be affected by weak instrument bias (Table 3). The results of the univariable MR analyses were based on the inverse-variance-weighted method. For ω-6 PUFAs, genetic predisposition to higher levels of AA or GLA was significantly associated with lower schizophrenia risk (ORAA = 0.986, p = 0.030; ORGLA = 0.148, p = 0.004), whereas ω-3 LCPUFAs showed only protective tendencies without statistical significance (Figure 4A). In contrast, the precursors of LCPUFAs, LA and ALA, showed the opposite effects on schizophrenia (ORLA = 1.008, p = 0.264; ORALA = 3.601, p = 0.072) (Figure 4A). Single-SNP analyses between AA and schizophrenia showed that the direction of the causal relationship driven by rs174547 was the same as the direction driven by rs16966952. The results of single-SNP analyses between GLA and schizophrenia were similar (Figure 4B,C).
These associations were driven by the SNP rs174547 in FADS1 and rs16966952 in NTAN1/PDXDC1. These genes are strongly connected with fatty acid metabolism. FADS1 is involved in the desaturation of LA and ALA, contributing to the biosynthesis of LCPUFAs such as AA and DHA. PDXDC1 expression level appeared to be changed in mice fed with a high-fat diet.
The pleiotropy assessment found that rs174547 was most strongly associated with AA and other multiple metabolic traits that are closely related to PUFAs, such as lipid metabolism phenotypes. The SNP rs16966952 was predominantly associated with DGLA and PUFA-related metabolic traits (Table S5). Given the instrumental variables related, most of the traits are downstream of the main exposure of interest. These metabolic relationships are considered as vertical pleiotropy and have no effects on the consequence of MR analysis. Other traits include pulse rate, peripheral blood cell count, height, asthma, hair, or balding pattern, none of which are considered as potential risk factors for schizophrenia. Moreover, neither of these SNPs was directly related to the other risk factors for schizophrenia at the genome-wide significance threshold (Table S5). As for some potential confounders such as CRP and educational attainment and vitamin D levels, there are some evidences to manifest AA-related SNPs would not affect the schizophrenia risk through these traits. Our LD analysis results show that all of the significant SNPs of CRP were independent with our instruments rs174547 and rs16966952 (R2 < 0.1), which implied that CRP is not genetically related to AA (Figure S1). A previous study reported that AA is not causally related to educational attainment and vitamin D [34]. These results suggest it is unlikely that the genetic instruments of AA and GLA affect schizophrenia via these confounders.
Considering the close relationships between PUFAs, we applied further MR analyses. Given the sharing of the elongating and desaturating enzymes by ω-6 and ω-3 PUFAs, multivariable MR analyses were applied to estimate the direct effects of one fatty acid in consideration of another. In the associations between schizophrenia and two PUFAs, AA and EPA, which were assessed together, there is evidence for a causal effect of AA as a protective factor of schizophrenia (ORAA = 0.981, p = 0.009), whereas no statistically significant evidence was found with EPA (Table 4). In the associations between schizophrenia and LA and ALA, which were assessed together, effect estimates for LA and ALA aligned well with the univariable MR estimates (ORLA = 1.005, p = 0.805; ORALA = 3.747, p = 0.512) (Table 4). Given the metabolic pathway in each PUFA group, further univariable MR analyses were applied to estimate the effect of a single PUFA within consideration of its preceding fatty acid. In this MR analysis, controlling for precursor PUFA, low schizophrenia risk was consistently observed with higher levels of AA and GLA (ORAA = 0.986, p = 0.021; ORGLA =0.240, p = 0.002) (Table S6).
Overall, the MR analyses indicate that AA and GLA are protective factors against schizophrenia. Consistent with our cross-national correlation results, insufficient AA may be a potential causal factor leading to the onset of schizophrenia.

4. Discussion

In this study, we aimed to investigate the relationships between PUFAs and schizophrenia risk using epidemiological methodology and Mendelian randomization analysis. Our results showed that ω-6 LCPUFAs, especially AA, are potential protective factors against schizophrenia. A higher intake of AA and ω-6 LCPUFAs, according to measures of food consumption, was significantly correlated with a lower incidence rate of schizophrenia across countries. Consistently, a genetically predicted low level of AA was revealed to increase schizophrenia risk by the MR analysis. In contrast, ω-3 LCPUFAs, including DHA and EPA, showed only a slight tendency to be protective against schizophrenia risk, which was an effect without statistical significance.
Both ω-3 and ω-6 LCPUFAs are rich in seafood and animal food. Previous studies have focused on the correlational relationships between the prevalence of affective disorders and seafood consumption [83,84] and found that lower seafood consumption was robustly associated with higher prevalence rates. Our study focused on the relationships between PUFAs and schizophrenia, and indicates that a high intake of ω-6 LCPUFAs, especially AA, was associated with a low risk of schizophrenia. Careful consideration was given to potential confounding factors. The relationship between each PUFA and schizophrenia was adjusted for based on country socioeconomic levels assessed by GDP per capita. Correspondingly, the results of our two-sample MR study, which is less prone to confounding than cross-national correlation studies, showed that AA was a protective factor against schizophrenia risk. Our MR results were based on large GWASs and the same race, which minimized the possibility of population stratification bias. The genes involved in the MR analysis, the FADS1 gene in LCPUFA synthesis, was reported to influence IQ performance, illustrating the crucial effects of PUFAs on early brain development [85]. NTAN1 is an N-terminal asparagine amidase that is related to social behavior and memory [86]. The results conducted by the MR analyses reveal that these variants may not directly lead to disease but instead affect schizophrenia risk by altering LCPUFAs, especially AA, and their subsequent traits. Therefore, these results both indicate that AA may play an important protective role against schizophrenia onset.
Although both ω-6 and ω-3 LCPUFAs are critical in brain function and development, our results show that ω-6 LCPUFAs, but not ω-3 LCPUFAs, function as protective factors against developing schizophrenia. Our findings provide supportive evidence for the hypothesis that the onset of schizophrenia might be determined by insufficient AA, a major ω-6 LCPUFA in the brain; this is released for the production of its eicosanoid metabolites to support adequate signal transduction [87], and when present at abnormal levels, mental activities due to the structural enrichment of AA and its vital functions in the brain will be disrupted.
AA deficiency impairs normal brain function [88] and contributes to schizophrenia-like phenotypes [89], and supplementation with AA can alleviate psychotic manifestations to some extent [9,10,90,91,92,93]. Lee et al. found that Lpiat1−/− mice, a model of AA deficiency in phosphatidylinositol synthesis, died within a month after birth and showed atrophy of the cerebral cortex and hippocampus [88]. M. Maekawa et al. found that gestational and early postnatal dietary deprivation of AA in mouse offspring elicited schizophrenia-like phenotypes in adulthood [89]. Meanwhile, M. Maekawa et al. [90] found that additional AA in breastmilk increases the total number of neurons in not only normal infant rats but also schizophrenia-like infant rats with information-processing problems, whereas few effects of DHA administration were observed. Other findings showed that high levels of AA helped improve cognitive development in children [9,10] and alleviate cognitive dysfunction in elderly individuals [91,92,93]. A recent study further confirmed this finding by observing that higher ω-6 LCPUFA levels in childhood, but not ω-3 LCPUFA, reduced the risk of psychotic experiences or psychotic disorder in adulthood [11]. These findings suggested that AA plays an essential role in neurodevelopment, which manifests as both a disruption in cognitive performance and psychotic symptoms during the onset of schizophrenia. Our findings provide supportive evidence from an epidemiological perspective and an MR-analysis-based causal inference.
In contrast to AA and ω-6 LCPUFAs, we found that ω-3 LCPUFAs were not significantly related to schizophrenia risk. Several cross-national studies reported strong relationships between fish consumption and mood disorder [83,84,94], which indicates that ω-3 LCPUFAs may have an impact on affective disorders. In contrast, there were no associations between seafood consumption and schizophrenia prevalence [84]. Furthermore, M. Maekawa et al. found that mice supplemented with DHA showed better behaviors in sustaining motivation or the propensity to respond, while mice supplemented with AA showed better performance in a cognitive processing task [89]. Meta-analyses have indicated that ω-3 LCPUFA supplementation may have beneficial effects on affective symptoms in major depressive disorder [95] and bipolar depression [96], whereas it failed to reveal therapeutic benefits on psychotic symptoms in schizophrenia [97]. Taken together, these findings also indicated that ω-3 LCPUFA showed limited influence on the onset of schizophrenia.

5. Conclusions

In conclusion, our findings, combining the results of diet-dependent epidemiological analyses and gene-dependent MR analyses, show that AA and ω-6 LCPUFAs are protective factors against schizophrenia risk. These results imply diet enrichment might help prevent the onset of schizophrenia. Studies on AA supplementation for individuals with ultra-high risk for psychosis would help further characterize the nature of the AA and schizophrenia associations.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/nu15051195/s1. Table S1. Food supply quantity (g per capita per day) in different country; Table S2. Characteristics of participants included in the PUFA meta-analysis GWAS; Table S3. Overall enrolled studies about incidence of schizophrenia in each country; Table S4. Enrolled studies about incidence of affective disorder and dietary PUFAs across countries; Table S5. Pleiotropy assessment of rs174547 in FDAS1 and rs16966952 in NTAN1/PDXDC1; Table S6. The instruments adjusted for the preceding fatty acid levels of three PUFAs and schizophrenia in MR analyses.; Figure S1. Results of linkage disequilibrium analysis between (A) rs174547 and rs10832027, rs1582763 in chromosome 11; (B) rs16966952 and rs10521222, rs1558902 in chromosome 16. Refs. [23,24,35,36,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82] cited in Supplementary Materials.

Author Contributions

Conceptualization, C.W., X.H., Y.G. and Y.Z.; methodology, Y.G., X.H. and Y.Z.; validation, X.H.; formal analysis, Y.G. and Y.Z.; investigation, Y.G.; writing—original draft preparation, Y.G., X.H. and D.W.; writing—review and editing, Y.G., X.H., D.W., J.J., M.L., Y.Q., X.Y., Y.Z. and C.W.; visualization, Y.G. and X.H.; supervision, C.W. and Y.Z.; project administration, J.Z.; funding acquisition, C.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Chinese National Programs for Brain Science and Brain-like Intelligence Technology (grant number: 2021ZD0200800); the National Natural Science Foundation of China (grant number: 81971254, 81771440); Shanghai Municipal Science and Technology Major Project (grant number: 2017SHZDZX01); the National Key R&D Program of China (grant number: 2016YFC1306800, 2016YFC1306900); and the Interdisciplinary Program of Shanghai Jiao Tong University (grant number: ZH2018ZDA40).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. All data are publicly available, no patient contact was made, and no individual identifiers were required. Therefore, ethical approval for the study was not required.

Data Availability Statement

Food consumption data were obtained from the Food and Agriculture Organization of the United Nations, and the nutrition information of food was obtained from the Australian Food Composition Database (previously called NUTTAB). GDP per capita in each country was obtained from the World Bank. The incidence rates of schizophrenia and PUFA related datasets used in Mendelian randomization analyses are from published studies, which are cited in the chapter “References”. The genetic statistical data of schizophrenia can be downloaded from the Psychiatry Genomics Consortium database.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mueser, K.T.; McGurk, S.R. Schizophrenia. Lancet 2004, 363, 2063–2072. [Google Scholar] [CrossRef]
  2. McGrath, J.; Saha, S.; Welham, J.; El Saadi, O.; MacCauley, C.; Chant, D. A systematic review of the incidence of schizophrenia: The distribution of rates and the influence of sex, urbanicity, migrant status and methodology. BMC Med. 2004, 2, 13. [Google Scholar] [CrossRef] [Green Version]
  3. Lewis, L. Mood disorders: Diagnosis, treatment, and support from a patient perspective. Psychopharmacol. Bull. 2001, 35, 186–196. [Google Scholar]
  4. Susser, E.; Neugebauer, R.; Hoek, H.W.; Brown, A.S.; Lin, S.; Labovitz, D.; Gorman, J.M. Schizophrenia after prenatal famine. Further evidence. Arch. Gen. Psychiatry 1996, 53, 25–31. [Google Scholar] [CrossRef] [PubMed]
  5. St Clair, D.; Xu, M.; Wang, P.; Yu, Y.; Fang, Y.; Zhang, F.; Zheng, X.; Gu, N.; Feng, G.; Sham, P.; et al. Rates of adult schizophrenia following prenatal exposure to the Chinese famine of 1959–1961. JAMA 2005, 294, 557–562. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Muskiet, F.A.J. Pathophysiology and Evolutionary Aspects of Dietary Fats and Long-Chain Polyunsaturated Fatty Acids across the Life Cycle. In Fat Detection: Taste, Texture, and Post Ingestive Effects; Frontiers in Neuroscience; Montmayeur, J.P., le Coutre, J., Eds.; CRC Press/Taylor & Francis: Boca Raton, FL, USA, 2010. [Google Scholar]
  7. Janssen, C.I.; Kiliaan, A.J. Long-chain polyunsaturated fatty acids (LCPUFA) from genesis to senescence: The influence of LCPUFA on neural development, aging, and neurodegeneration. Prog. Lipid. Res. 2014, 53, 1–17. [Google Scholar] [CrossRef] [PubMed]
  8. Lepping, R.J.; Honea, R.A.; Martin, L.E.; Liao, K.; Choi, I.Y.; Lee, P.; Papa, V.B.; Brooks, W.M.; Shaddy, D.J.; Carlson, S.E.; et al. Long-chain polyunsaturated fatty acid supplementation in the first year of life affects brain function, structure, and metabolism at age nine years. Dev. Psychobiol. 2019, 61, 5–16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Birch, E.E.; Garfield, S.; Castaneda, Y.; Hughbanks-Wheaton, D.; Uauy, R.; Hoffman, D. Visual acuity and cognitive outcomes at 4 years of age in a double-blind, randomized trial of long-chain polyunsaturated fatty acid-supplemented infant formula. Early Hum. Dev. 2007, 83, 279–284. [Google Scholar] [CrossRef]
  10. Willatts, P.; Forsyth, S.; Agostoni, C.; Casaer, P.; Riva, E.; Boehm, G. Effects of long-chain PUFA supplementation in infant formula on cognitive function in later childhood. Am. J. Clin. Nutr. 2013, 98, 536S–542S. [Google Scholar] [CrossRef] [Green Version]
  11. Thompson, A.D.; Jones, H.J.; Heron, J.; Hibbeln, J.; Sullivan, S.; Zammit, S. Omega-3 and Omega-6 fatty acids and risk of psychotic outcomes in the ALSPAC birth cohort. Schizophr. Res. 2020, 224, 108–115. [Google Scholar] [CrossRef]
  12. Hoen, W.P.; Lijmer, J.G.; Duran, M.; Wanders, R.J.; van Beveren, N.J.; de Haan, L. Red blood cell polyunsaturated fatty acids measured in red blood cells and schizophrenia: A meta-analysis. Psychiatry Res. 2013, 207, 1–12. [Google Scholar] [CrossRef] [PubMed]
  13. Van der Kemp, W.J.; Klomp, D.W.; Kahn, R.S.; Luijten, P.R.; Hulshoff Pol, H.E. A meta-analysis of the polyunsaturated fatty acid composition of erythrocyte membranes in schizophrenia. Schizophr. Res. 2012, 141, 153–161. [Google Scholar] [CrossRef] [PubMed]
  14. Alqarni, A.; Mitchell, T.W.; McGorry, P.D.; Nelson, B.; Markulev, C.; Yuen, H.P.; Schafer, M.R.; Berger, M.; Mossaheb, N.; Schlogelhofer, M.; et al. Comparison of erythrocyte omega-3 index, fatty acids and molecular phospholipid species in people at ultra-high risk of developing psychosis and healthy people. Schizophr. Res. 2020, 226, 44–51. [Google Scholar] [CrossRef] [PubMed]
  15. Charlson, F.J.; Ferrari, A.J.; Santomauro, D.F.; Diminic, S.; Stockings, E.; Scott, J.G.; McGrath, J.J.; Whiteford, H.A. Global Epidemiology and Burden of Schizophrenia: Findings From the Global Burden of Disease Study 2016. Schizophr. Bull. 2018, 44, 1195–1203. [Google Scholar] [CrossRef] [PubMed]
  16. Hogerzeil, S.J.; van Hemert, A.M. Design choices when estimating the treated incidence of schizophrenia. Psychol. Med. 2019, 49, 2635–2636. [Google Scholar] [CrossRef] [Green Version]
  17. Schmidhuber, J.; Sur, P.; Fay, K.; Huntley, B.; Salama, J.; Lee, A.; Cornaby, L.; Horino, M.; Murray, C.; Afshin, A. The Global Nutrient Database: Availability of macronutrients and micronutrients in 195 countries from 1980 to 2013. Lancet Planet Health 2018, 2, e353–e368. [Google Scholar] [CrossRef] [Green Version]
  18. Forsyth, S.; Gautier, S.; Salem, N., Jr. Global Estimates of Dietary Intake of Docosahexaenoic Acid and Arachidonic Acid in Developing and Developed Countries. Ann. Nutr. Metab. 2016, 68, 258–267. [Google Scholar] [CrossRef] [Green Version]
  19. Lassek, W.D.; Gaulin, S.J. Linoleic and docosahexaenoic acids in human milk have opposite relationships with cognitive test performance in a sample of 28 countries. Prostaglandins Leukot. Essent. Fat. Acids 2014, 91, 195–201. [Google Scholar] [CrossRef]
  20. Basu, S.; McKee, M.; Galea, G.; Stuckler, D. Relationship of soft drink consumption to global overweight, obesity, and diabetes: A cross-national analysis of 75 countries. Am. J. Public Health 2013, 103, 2071–2077. [Google Scholar] [CrossRef]
  21. Fox, A.; Feng, W.; Asal, V. What is driving global obesity trends? Globalization or “modernization”? Glob. Health 2019, 15, 32. [Google Scholar] [CrossRef] [Green Version]
  22. Malik, V.S.; Willett, W.C.; Hu, F.B. Global obesity: Trends, risk factors and policy implications. Nat. Rev. Endocrinol. 2013, 9, 13–27. [Google Scholar] [CrossRef]
  23. Guan, W.; Steffen, B.T.; Lemaitre, R.N.; Wu, J.H.Y.; Tanaka, T.; Manichaikul, A.; Foy, M.; Rich, S.S.; Wang, L.; Nettleton, J.A.; et al. Genome-wide association study of plasma N6 polyunsaturated fatty acids within the cohorts for heart and aging research in genomic epidemiology consortium. Circ. Cardiovasc. Genet. 2014, 7, 321–331. [Google Scholar] [CrossRef] [Green Version]
  24. Lemaitre, R.N.; Tanaka, T.; Tang, W.; Manichaikul, A.; Foy, M.; Kabagambe, E.K.; Nettleton, J.A.; King, I.B.; Weng, L.C.; Bhattacharya, S.; et al. Genetic loci associated with plasma phospholipid n-3 fatty acids: A meta-analysis of genome-wide association studies from the CHARGE Consortium. PLoS Genet. 2011, 7, e1002193. [Google Scholar] [CrossRef] [Green Version]
  25. Ghoneim, D.H.; Zhu, J.; Zheng, W.; Long, J.; Murff, H.J.; Ye, F.; Setiawan, V.W.; Wilkens, L.R.; Khankari, N.K.; Haycock, P.; et al. Mendelian Randomization Analysis of n-6 Polyunsaturated Fatty Acid Levels and Pancreatic Cancer Risk. Cancer Epidemiol. Biomark. Prev. 2020, 29, 2735–2739. [Google Scholar] [CrossRef]
  26. Yuan, S.; Larsson, S.C. Association of genetic variants related to plasma fatty acids with type 2 diabetes mellitus and glycaemic traits: A Mendelian randomisation study. Diabetologia 2020, 63, 116–123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Khankari, N.K.; Banbury, B.L.; Borges, M.C.; Haycock, P.; Albanes, D.; Arndt, V.; Berndt, S.I.; Bezieau, S.; Brenner, H.; Campbell, P.T.; et al. Mendelian Randomization of Circulating Polyunsaturated Fatty Acids and Colorectal Cancer Risk. Cancer Epidemiol. Biomark. Prev. 2020, 29, 860–870. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Tomata, Y.; Larsson, S.C.; Hagg, S. Polyunsaturated fatty acids and risk of Alzheimer’s disease: A Mendelian randomization study. Eur. J. Nutr. 2020, 59, 1763–1766. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Schizophrenia Working Group of the Psychiatric Genomics Consortium. Biological insights from 108 schizophrenia-associated genetic loci. Nature 2014, 511, 421–427. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Burgess, S.; Thompson, S.G.; Collaboration, C.C.G. Avoiding bias from weak instruments in Mendelian randomization studies. Int. J. Epidemiol. 2011, 40, 755–764. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Fernandes, B.S.; Steiner, J.; Bernstein, H.G.; Dodd, S.; Pasco, J.A.; Dean, O.M.; Nardin, P.; Goncalves, C.A.; Berk, M. C-reactive protein is increased in schizophrenia but is not altered by antipsychotics: Meta-analysis and implications. Mol. Psychiatry 2016, 21, 554–564. [Google Scholar] [CrossRef] [PubMed]
  32. Dickson, H.; Hedges, E.P.; Ma, S.Y.; Cullen, A.E.; MacCabe, J.H.; Kempton, M.J.; Downs, J.; Laurens, K.R. Academic achievement and schizophrenia: A systematic meta-analysis. Psychol. Med. 2020, 50, 1949–1965. [Google Scholar] [CrossRef] [PubMed]
  33. Radua, J.; Ramella-Cravaro, V.; Ioannidis, J.P.A.; Reichenberg, A.; Phiphopthatsanee, N.; Amir, T.; Yenn Thoo, H.; Oliver, D.; Davies, C.; Morgan, C.; et al. What causes psychosis? An umbrella review of risk and protective factors. World Psychiatry 2018, 17, 49–66. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Seviiri, M.; Law, M.H.; Ong, J.S.; Gharahkhani, P.; Nyholt, D.R.; Olsen, C.M.; Whiteman, D.C.; MacGregor, S. Polyunsaturated Fatty Acid Levels and the Risk of Keratinocyte Cancer: A Mendelian Randomization Analysis. Cancer Epidemiol. Biomark. Prev. 2021, 30, 1591–1598. [Google Scholar] [CrossRef] [PubMed]
  35. Machiela, M.J.; Chanock, S.J. LDlink: A web-based application for exploring population-specific haplotype structure and linking correlated alleles of possible functional variants. Bioinformatics 2015, 31, 3555–3557. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Ligthart, S.; Vaez, A.; Vosa, U.; Stathopoulou, M.G.; de Vries, P.S.; Prins, B.P.; Van der Most, P.J.; Tanaka, T.; Naderi, E.; Rose, L.M.; et al. Genome Analyses of >200,000 Individuals Identify 58 Loci for Chronic Inflammation and Highlight Pathways that Link Inflammation and Complex Disorders. Am. J. Hum. Genet. 2018, 103, 691–706. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Kamat, M.A.; Blackshaw, J.A.; Young, R.; Surendran, P.; Burgess, S.; Danesh, J.; Butterworth, A.S.; Staley, J.R. PhenoScanner V2: An expanded tool for searching human genotype-phenotype associations. Bioinformatics 2019, 35, 4851–4853. [Google Scholar] [CrossRef] [Green Version]
  38. Messias, E.; Sampaio, J.J.; Messias, N.C.; Kirkpatrick, B. Epidemiology of schizophrenia in northeast Brazil. J. Nerv. Ment. Dis. 2000, 188, 118–120. [Google Scholar] [CrossRef]
  39. Whitehorn, D.; Richard, J.C.; Kopala, L.C. Hospitalization in the first year of treatment for schizophrenia. Can. J. Psychiatry 2004, 49, 635–638. [Google Scholar] [CrossRef] [Green Version]
  40. Larach, V.; Varela, L.F.; Vergara, G.; Rojas, I.; Hernández, V.; Casanueva, R.; Irarrázaval, M. Universal healthcare coverage for first episode of schizophrenia-spectrum disorders in Chile: Analysis of the administrative database. Eur. Arch. Psychiatry Clin. Neurosci. 2022, 272, 129–138. [Google Scholar] [CrossRef]
  41. Handal, N.; Dodds, J.H. Statistics for initial admission for schizophrenia in hospitals in Costa Rica. Rev. Panam Salud Publica 1997, 1, 426–434. [Google Scholar]
  42. Folnegović, Z.; Folnegović-Šmalc, V.; Kulčar, Ž. The incidence of schizophrenia in Croatia. Br. J. Psychiatry 1990, 156, 363–365. [Google Scholar] [CrossRef]
  43. Castagnini, A.; Foldager, L. Variations in incidence and age of onset of acute and transient psychotic disorders. Soc. Psychiatry Psychiatr. Epidemiol. 2013, 48, 1917–1922. [Google Scholar] [CrossRef] [PubMed]
  44. Salokangas, R.K.R.; Helminen, M.; Koivisto, A.-M.; Rantanen, H.; Oja, H.; Pirkola, S.; Wahlbeck, K.; Joukamaa, M. Incidence of hospitalised schizophrenia in Finland since 1980: Decreasing and increasing again. Soc. Psychiatry Psychiatr. Epidemiol. 2011, 46, 343–350. [Google Scholar] [CrossRef] [PubMed]
  45. Löffler, W.; Häfner, H. Ecological pattern of first admitted schizophrenics in two German cities over 25 years. Soc. Sci. Med. 1999, 49, 93–108. [Google Scholar] [CrossRef] [PubMed]
  46. Helgason, L. Psychiatric services and mental illness in Iceland. Incidence study (1966–1967) with 6–7 year follow-up. Acta Psychiatr. Scand. Suppl. 1977, 268, 1–140. [Google Scholar]
  47. Rajkumar, S.; Padmavathi, R.; Thara, R.; Menon, M.S. Incidence of schizophrenia in an urban community in madras. Indian J. Psychiatry 1993, 35, 18–21. [Google Scholar] [PubMed]
  48. Blake, B.; Halpenny, J.V.; ni Nuallain, M.; O’Brien, P.F.; O’Hare, A.; Walsh, D. The incidence of mental illness in Ireland—Patients contacting psychiatric services in three Irish Counties. Ir. J. Psychiatry 1984, 5, 23–29. [Google Scholar]
  49. Mulè, A.; Sideli, L.; Capuccio, V.; Fearon, P.; Ferraro, L.; Kirkbride, J.B.; La Cascia, C.; Sartorio, C.; Seminerio, F.; Tripoli, G.; et al. Low incidence of psychosis in Italy: Confirmation from the first epidemiological study in Sicily. Soc. Psychiatry Psychiatr. Epidemiol. 2017, 52, 155–162. [Google Scholar] [CrossRef] [Green Version]
  50. Ohta, Y.; Nakane, Y.; Nishihara, J.; Takemoto, T. Ecological structure and incidence rates of schizophrenia in Nagasaki City. Acta Psychiatr. Scand. 1992, 86, 113–120. [Google Scholar] [CrossRef]
  51. Van Os, J.; Driessen, G.; Gunther, N.; Delespaul, P. Neighbourhood variation in incidence of schizophrenia. Evidence for person-environment interaction. Br. J. Psychiatry 2000, 176, 243–248. [Google Scholar] [CrossRef]
  52. Tapsell, R.; Hallett, C.; Mellsop, G. The rate of mental health service use in New Zealand as analysed by ethnicity. Australas Psychiatry 2018, 26, 290–293. [Google Scholar] [CrossRef] [PubMed]
  53. Johannessen, J.O. Schizofreni-forekomst i et norsk fylke, Rogaland = The occurrence of schizophrenia in the Norwegian county of Rogaland. Nord. Psykiatr. Tidsskr. 1985, 39, 217–223. [Google Scholar] [CrossRef]
  54. Jung, Y.S.; Kim, Y.E.; Go, D.S.; Yoon, S.J. The prevalence, incidence, and admission rate of diagnosed schizophrenia spectrum disorders in Korea, 2008–2017: A nationwide population-based study using claims big data analysis. PLoS ONE 2021, 16, e0256221. [Google Scholar] [CrossRef] [PubMed]
  55. Kazakovtsev, B.A.; Sidoryuk, O.V.; Zrazhevskaya, I.A.; Ovsyannikov, S.A. Regional differences in the incidence of affective disorders. Zh. Nevrol. Psikhiatr. Im S S Korsakova 2019, 119, 10–16. [Google Scholar] [CrossRef]
  56. Vázquez-Barquero, J.L.; Cuesta Nuñez, M.J.; de la Varga, M.; Herrera Castanedo, S.; Gaite, L.; Arenal, A. The Cantabria first episode schizophrenia study: A summary of general findings. Acta Psychiatr. Scand. 1995, 91, 156–162. [Google Scholar] [CrossRef]
  57. Jörgensen, L.; Ahlbom, A.; Allebeck, P.; Dalman, C. The Stockholm non-affective psychoses study (snaps): The importance of including out-patient data in incidence studies. Acta Psychiatr. Scand. 2010, 121, 389–392. [Google Scholar] [CrossRef]
  58. Ajdacic-Gross, V.; Lauber, C.; Warnke, I.; Haker, H.; Murray, R.M.; Rössler, W. Changing incidence of psychotic disorders among the young in Zurich. Schizophr. Res. 2007, 95, 9–18. [Google Scholar] [CrossRef]
  59. Loganovsky, K.N.; Loganovskaja, T.K. At issue: Schizophrenia spectrum disorders in persons exposed to ionizing radiation as a result of the Chernobyl accident. Schizophr. Bull. 2000, 26, 751–773. [Google Scholar] [CrossRef] [Green Version]
  60. Reay, R.; Mitford, E.; McCabe, K.; Paxton, R.; Turkington, D. Incidence and diagnostic diversity in first-episode psychosis. Acta Psychiatr. Scand. 2010, 121, 315–319. [Google Scholar] [CrossRef]
  61. Simon, G.E.; Coleman, K.J.; Yarborough, B.J.H.; Operskalski, B.; Stewart, C.; Hunkeler, E.M.; Lynch, F.; Carrell, D.; Beck, A. First Presentation With Psychotic Symptoms in a Population-Based Sample. Psychiatr. Serv. 2017, 68, 456–461. [Google Scholar] [CrossRef] [Green Version]
  62. Gynther, B.; Charlson, F.; Obrecht, K.; Waller, M.; Santomauro, D.; Whiteford, H.; Hunter, E. The Epidemiology of Psychosis in Indigenous Populations in Cape York and the Torres Strait. EClinicalMedicine 2019, 10, 68–77. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Iacono, W.G.; Beiser, M. Are males more likely than females to develop schizophrenia? Am. J. Psychiatry 1992, 149, 1070–1074. [Google Scholar] [CrossRef] [PubMed]
  64. Nicole, L.; Lesage, A.; Lalonde, P. Lower incidence and increased male:female ratio in schizophrenia. Br. J. Psychiatry 1992, 161, 556–557. [Google Scholar] [CrossRef] [PubMed]
  65. Salokangas, R.K.R. First-contact rate for schizophrenia in community psychiatric care. Consideration of the oestrogen hypothesis. Eur. Arch. Psychiatry Clin. Neurosci. 1993, 242, 337–346. [Google Scholar] [CrossRef]
  66. Suokas, K.; Koivisto, A.M.; Hakulinen, C.; Kaltiala, R.; Sund, R.; Lumme, S.; Kampman, O.; Pirkola, S. Association of Income With the Incidence Rates of First Psychiatric Hospital Admissions in Finland, 1996–2014. JAMA Psychiatry 2020, 77, 274–284. [Google Scholar] [CrossRef]
  67. De Salvia, D.; Barbato, A.; Salvo, P.; Zadro, F. Prevalence and incidence of schizophrenic disorders in Portogruaro. An Italian case register study. J. Nerv. Ment. Dis. 1993, 181, 275–282. [Google Scholar] [CrossRef]
  68. Preti, A.; Miotto, P. Increase in first admissions for schizophrenia and other major psychoses in Italy. Psychiatry Res. 2000, 94, 139–152. [Google Scholar] [CrossRef]
  69. Lasalvia, A.; Bonetto, C.; Tosato, S.; Zanatta, G.; Cristofalo, D.; Salazzari, D.; Lazzarotto, L.; Bertani, M.; Bissoli, S.; De Santi, K.; et al. First-contact incidence of psychosis in north-eastern Italy: Influence of age, gender, immigration and socioeconomic deprivation. Br. J. Psychiatry 2014, 205, 127–134. [Google Scholar] [CrossRef] [Green Version]
  70. Tarricone, I.; Mimmi, S.; Paparelli, A.; Rossi, E.; Mori, E.; Panigada, S.; Carchia, G.; Bandieri, V.; Michetti, R.; Minenna, G.; et al. First-episode psychosis at the West Bologna Community Mental Health Centre: Results of an 8-year prospective study. Psychol. Med. 2012, 42, 2255–2264. [Google Scholar] [CrossRef] [Green Version]
  71. Park Jae-Sup, C.W.-J.; Dong-Wook, K.; Jung-Eun, S.; Seul-Bi, L.; San-Hyeop, L. Status of Treatment for First Episode Psychosis in Korea and Need for Expansion of Coverage; Report from 2017 NHIS Ilsan Hospital; National Health Insurance Service Ilsan Hospital: Goyang, Republic of Korea, 2017; p. 113. [Google Scholar]
  72. Orrico-Sánchez, A.; López-Lacort, M.; Munõz-Quiles, C.; Sanfélix-Gimeno, G.; Diéz-Domingo, J. Epidemiology of schizophrenia and its management over 8-years period using real-world data in Spain. BMC Psychiatry 2020, 20, 149. [Google Scholar] [CrossRef] [Green Version]
  73. Osby, U.; Hammar, N.; Brandt, L.; Wicks, S.; Thinsz, Z.; Ekbom, A.; Sparén, P. Time trends in first admissions for schizophrenia and paranoid psychosis in Stockholm County, Sweden. Schizophr. Res. 2001, 47, 247–254. [Google Scholar] [CrossRef] [PubMed]
  74. Silva, B.; Golay, P.; Boubaker, K.; Bonsack, C.; Morandi, S. Community treatment orders in Western Switzerland: A retrospective epidemiological study. Int. J. Law Psychiatry 2019, 67, 101509. [Google Scholar] [CrossRef] [PubMed]
  75. Brewin, J.; Cantwell, R.; Dalkin, T.; Fox, R.; Medley, I.; Glazebrook, C.; Kwiecinski, R.; Harrison, G. Incidence of schizophrenia in Nottingham: A comparison of two cohorts, 1978–80 and 1992–94. Br. J. Psychiatry 1997, 171, 140–144. [Google Scholar] [CrossRef] [PubMed]
  76. Goldacre, M.; Shiwach, R.; Yeates, D. Estimating incidence and prevalence of treated psychiatric disorders from routine statistics: The example of schizophrenia in Oxfordshire. J. Epidemiol. Community Health 1994, 48, 318–322. [Google Scholar] [CrossRef] [Green Version]
  77. Lee, S.C.; DelPozo-Banos, M.; Lloyd, K.; Jones, I.; Walters, J.T.R.; Owen, M.J.; O’Donovan, M.; John, A. Area deprivation, urbanicity, severe mental illness and social drift—A population-based linkage study using routinely collected primary and secondary care data. Schizophr. Res. 2020, 220, 130–140. [Google Scholar] [CrossRef]
  78. Singh, S.P.; Burns, T.; Amin, S.; Jones, P.B.; Harrison, G. Acute and transient psychotic disorders: Precursors, epidemiology, course and outcome. Br. J. Psychiatry 2004, 185, 452–459. [Google Scholar] [CrossRef]
  79. Jensen, C.M.; Steinhausen, H.C. Time Trends in Lifetime Incidence Rates of First-Time Diagnosed Bipolar and Depressive Disorders Across 16 Years in Danish Psychiatric Hospitals: A Nationwide Study. J. Clin. Psychiatry 2016, 77, e1570–e1575. [Google Scholar] [CrossRef] [Green Version]
  80. Markkula, N.; Lehti, V.; Gissler, M.; Suvisaari, J. Incidence and prevalence of mental disorders among immigrants and native Finns: A register-based study. Soc. Psychiatry Psychiatr. Epidemiol. 2017, 52, 1523–1540. [Google Scholar] [CrossRef] [Green Version]
  81. Selten, J.P.; van Os, J.; Nolen, W.A. First admissions for mood disorders in immigrants to the Netherlands. Soc. Psychiatry Psychiatr. Epidemiol. 2003, 38, 547–550. [Google Scholar] [CrossRef]
  82. Carlborg, A.; Ferntoft, L.; Thuresson, M.; Bodegard, J. Population study of disease burden, management, and treatment of bipolar disorder in Sweden: A retrospective observational registry study. Bipolar Disord. 2015, 17, 76–85. [Google Scholar] [CrossRef]
  83. Hibbeln, J.R. Fish consumption and major depression. Lancet 1998, 351, 1213. [Google Scholar] [CrossRef]
  84. Noaghiul, S.; Hibbeln, J.R. Cross-national comparisons of seafood consumption and rates of bipolar disorders. Am. J. Psychiatry 2003, 160, 2222–2227. [Google Scholar] [CrossRef] [PubMed]
  85. Caspi, A.; Williams, B.; Kim-Cohen, J.; Craig, I.W.; Milne, B.J.; Poulton, R.; Schalkwyk, L.C.; Taylor, A.; Werts, H.; Moffitt, T.E. Moderation of breastfeeding effects on the IQ by genetic variation in fatty acid metabolism. Proc. Natl. Acad. Sci. USA 2007, 104, 18860–18865. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Ingason, A.; Rujescu, D.; Cichon, S.; Sigurdsson, E.; Sigmundsson, T.; Pietilainen, O.P.; Buizer-Voskamp, J.E.; Strengman, E.; Francks, C.; Muglia, P.; et al. Copy number variations of chromosome 16p13.1 region associated with schizophrenia. Mol. Psychiatry 2011, 16, 17–25. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Horrobin, D.F. The membrane phospholipid hypothesis as a biochemical basis for the neurodevelopmental concept of schizophrenia. Schizophr. Res. 1998, 30, 193–208. [Google Scholar] [CrossRef] [PubMed]
  88. Lee, H.C.; Inoue, T.; Sasaki, J.; Kubo, T.; Matsuda, S.; Nakasaki, Y.; Hattori, M.; Tanaka, F.; Udagawa, O.; Kono, N.; et al. LPIAT1 regulates arachidonic acid content in phosphatidylinositol and is required for cortical lamination in mice. Mol. Biol. Cell 2012, 23, 4689–4700. [Google Scholar] [CrossRef] [PubMed]
  89. Maekawa, M.; Watanabe, A.; Iwayama, Y.; Kimura, T.; Hamazaki, K.; Balan, S.; Ohba, H.; Hisano, Y.; Nozaki, Y.; Ohnishi, T.; et al. Polyunsaturated fatty acid deficiency during neurodevelopment in mice models the prodromal state of schizophrenia through epigenetic changes in nuclear receptor genes. Transl. Psychiatry 2017, 7, e1229. [Google Scholar] [CrossRef]
  90. Maekawa, M.; Takashima, N.; Matsumata, M.; Ikegami, S.; Kontani, M.; Hara, Y.; Kawashima, H.; Owada, Y.; Kiso, Y.; Yoshikawa, T.; et al. Arachidonic acid drives postnatal neurogenesis and elicits a beneficial effect on prepulse inhibition, a biological trait of psychiatric illnesses. PLoS ONE 2009, 4, e5085. [Google Scholar] [CrossRef] [Green Version]
  91. Kotani, S.; Sakaguchi, E.; Warashina, S.; Matsukawa, N.; Ishikura, Y.; Kiso, Y.; Sakakibara, M.; Yoshimoto, T.; Guo, J.; Yamashima, T. Dietary supplementation of arachidonic and docosahexaenoic acids improves cognitive dysfunction. Neurosci. Res. 2006, 56, 159–164. [Google Scholar] [CrossRef]
  92. Tokuda, H.; Ito, M.; Sueyasu, T.; Sasaki, H.; Morita, S.; Kaneda, Y.; Rogi, T.; Kondo, S.; Kouzaki, M.; Tsukiura, T.; et al. Effects of combining exercise with long-chain polyunsaturated fatty acid supplementation on cognitive function in the elderly: A randomised controlled trial. Sci. Rep. 2020, 10, 12906. [Google Scholar] [CrossRef]
  93. MacDonald-Wicks, L.; McEvoy, M.; Magennis, E.; Schofield, P.W.; Patterson, A.J.; Zacharia, K. Dietary Long-Chain Fatty Acids and Cognitive Performance in Older Australian Adults. Nutrients 2019, 11, 711. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Hibbeln, J.R. Seafood consumption, the DHA content of mothers’ milk and prevalence rates of postpartum depression: A cross-national, ecological analysis. J. Affect. Disord. 2002, 69, 15–29. [Google Scholar] [CrossRef] [PubMed]
  95. Grosso, G.; Pajak, A.; Marventano, S.; Castellano, S.; Galvano, F.; Bucolo, C.; Drago, F.; Caraci, F. Role of omega-3 fatty acids in the treatment of depressive disorders: A comprehensive meta-analysis of randomized clinical trials. PLoS ONE 2014, 9, e96905. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Sarris, J.; Mischoulon, D.; Schweitzer, I. Omega-3 for bipolar disorder: Meta-analyses of use in mania and bipolar depression. J. Clin. Psychiatry 2012, 73, 81–86. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  97. Fenton, W.S.; Dickerson, F.; Boronow, J.; Hibbeln, J.R.; Knable, M. A placebo-controlled trial of omega-3 fatty acid (ethyl eicosapentaenoic acid) supplementation for residual symptoms and cognitive impairment in schizophrenia. Am. J. Psychiatry 2001, 158, 2071–2074. [Google Scholar] [CrossRef]
Nutrients 15 01195 g001 550
Figure 1. Schematic diagram of our study design. PUFA—polyunsaturated fatty acid; FAO—Food and Agriculture Organization of the United Nations; NUTTAB—Australian Food Composition Database (previously called NUTTAB); GWAS—genome-wide association study; SNP—single-nucleotide polymorphism; CHARGE—Cohorts for Heart and Aging Research in Genomic Epidemiology; LA—linoleic acid; GLA—γ-linolenic acid; DGLA—dihomo-γ-linolenic acid; AA—arachidonic acid; ALA—α-linolenic acid; EPA—eicosapentaenoic acid; DPA—docosapentaenoic acid; DHA—docosahexaenoic acid; and PGC—Psychiatry Genomics Consortium.
Figure 1. Schematic diagram of our study design. PUFA—polyunsaturated fatty acid; FAO—Food and Agriculture Organization of the United Nations; NUTTAB—Australian Food Composition Database (previously called NUTTAB); GWAS—genome-wide association study; SNP—single-nucleotide polymorphism; CHARGE—Cohorts for Heart and Aging Research in Genomic Epidemiology; LA—linoleic acid; GLA—γ-linolenic acid; DGLA—dihomo-γ-linolenic acid; AA—arachidonic acid; ALA—α-linolenic acid; EPA—eicosapentaenoic acid; DPA—docosapentaenoic acid; DHA—docosahexaenoic acid; and PGC—Psychiatry Genomics Consortium.
Nutrients 15 01195 g001
Nutrients 15 01195 g002 550
Figure 2. Incidence rates and estimated daily LCPUFA intakes for countries. Darker blue signifies higher incidence rates of schizophrenia and lighter blue signify lower rates. Grey tones indicate that the incidence rates of schizophrenia in these countries were not available in this study. The size of the pie chart is proportional to the dietary intake of total LCPUFAs. The NUTTAB contained 11 LCPUFAs in total, which green colors represent ω-3 LCPUFAs and orange colors represent ω-6 LCPUFAs. The area outlined by the dark-blue box zoomed into the lower left corner of this figure. BRA—Brazil; CAN—Canada; CRI—Costa Rica; HRV—Croatia; DNK—Denmark; FIN—Finland; DEU—Germany; ISL—Iceland; IND—India; IRL—Ireland; ITA—Italy; JPN—Japan; NLD—Netherlands; NZL—New Zealand; NOR—Norway; RUS—Russian Federation; ESP—Spain; SWE—Sweden; CHE—Switzerland; KOR—the Republic of Korea; UKR—Ukraine; GBR—United Kingdom; and USA—United States.
Figure 2. Incidence rates and estimated daily LCPUFA intakes for countries. Darker blue signifies higher incidence rates of schizophrenia and lighter blue signify lower rates. Grey tones indicate that the incidence rates of schizophrenia in these countries were not available in this study. The size of the pie chart is proportional to the dietary intake of total LCPUFAs. The NUTTAB contained 11 LCPUFAs in total, which green colors represent ω-3 LCPUFAs and orange colors represent ω-6 LCPUFAs. The area outlined by the dark-blue box zoomed into the lower left corner of this figure. BRA—Brazil; CAN—Canada; CRI—Costa Rica; HRV—Croatia; DNK—Denmark; FIN—Finland; DEU—Germany; ISL—Iceland; IND—India; IRL—Ireland; ITA—Italy; JPN—Japan; NLD—Netherlands; NZL—New Zealand; NOR—Norway; RUS—Russian Federation; ESP—Spain; SWE—Sweden; CHE—Switzerland; KOR—the Republic of Korea; UKR—Ukraine; GBR—United Kingdom; and USA—United States.
Nutrients 15 01195 g002
Nutrients 15 01195 g003 550
Figure 3. Correlations between schizophrenia incidence rates and (A) AA, (B) DHA, (C) total ω-6 LCPUFA, and (D) total ω-3 LCPUFA. Correlations described the relationship differences for AA (rAA = −0.577, p < 0.01), DHA (rDHA = −0.169, p > 0.05), total ω-6 LCPUFA (rω-6 LCPUFA = −0.626, p < 0.001), and total ω-3 LCPUFA (rω-3 LCPUFA = −0.104, p > 0.05).
Figure 3. Correlations between schizophrenia incidence rates and (A) AA, (B) DHA, (C) total ω-6 LCPUFA, and (D) total ω-3 LCPUFA. Correlations described the relationship differences for AA (rAA = −0.577, p < 0.01), DHA (rDHA = −0.169, p > 0.05), total ω-6 LCPUFA (rω-6 LCPUFA = −0.626, p < 0.001), and total ω-3 LCPUFA (rω-3 LCPUFA = −0.104, p > 0.05).
Nutrients 15 01195 g003
Nutrients 15 01195 g004 550
Figure 4. Associations of genetically predicted PUFA levels with schizophrenia. (A) This figure shows the Mendelian randomization inverse-variance-weighted (IVW) results for the relationship between the eight PUFAs and the risk for schizophrenia. (B,C) Forest plots of the individual SNP effects for genetically predicted GLA and AA levels on schizophrenia risk. These figures illustrate the contribution of individual instrumental variables used for the association between (B) GLA, (C) AA, and schizophrenia risk. PUFA—polyunsaturated fatty acid; OR—odds ratio per SD; and 95% CI—95% confidence interval. The middle line represents the null (OR = 1.000), and the error bars represent 95% CI.
Figure 4. Associations of genetically predicted PUFA levels with schizophrenia. (A) This figure shows the Mendelian randomization inverse-variance-weighted (IVW) results for the relationship between the eight PUFAs and the risk for schizophrenia. (B,C) Forest plots of the individual SNP effects for genetically predicted GLA and AA levels on schizophrenia risk. These figures illustrate the contribution of individual instrumental variables used for the association between (B) GLA, (C) AA, and schizophrenia risk. PUFA—polyunsaturated fatty acid; OR—odds ratio per SD; and 95% CI—95% confidence interval. The middle line represents the null (OR = 1.000), and the error bars represent 95% CI.
Nutrients 15 01195 g004
Table
Table 1. Information of data on PUFAs and schizophrenia of MR analysis.
Table 1. Information of data on PUFAs and schizophrenia of MR analysis.
PhenotypeSample SizeAncestryWomen%PUFA Concentration, % of Total Fatty Acid
ω-6 PUFAs a8631Europe51~62%LA (19.96~24.78%), GLA (0.11~0.12%), DGLA (3.13~3.33%), AA (8.00~12.10%)
ω-3 PUFAs b8866Europe51~61%ALA (0.14~0.44%), DHA (2.29~3.66%), EPA (0.56~0.88%), DPA (0.83~0.95%)
Schizophrenia c34,241 cases,
45,604 controls		93.5% European
6.5% East Asian		8~73%Not applicable
a Summary statistic for GWAS of ω-6 PUFA from “Genome-wide association study of plasma N6 polyunsaturated fatty acids within the cohorts for heart and aging research in genomic epidemiology consortium” [23]. b Summary statistics for GWAS of ω-3 PUFA from “Genetic loci associated with plasma phospholipid n-3 fatty acids: a meta-analysis of genome-wide association studies from the CHARGE Consortium” [24]. c Summary statistics for GWAS of schizophrenia from “Biological insights from 108 schizophrenia-associated genetic loci” [29]. PUFA—polyunsaturated fatty acid; LA—linoleic acid; GLA—γ-linolenic acid; DGLA—dihomo-γ-linolenic acid; AA—arachidonic acid; ALA—α-linolenic acid; EPA—eicosapentaenoic acid; DPA—docosapentaenoic acid; and DHA—docosahexaenoic acid.
Table
Table 2. The incidence rates of schizophrenia, PUFA consumption, and economic development across different countries.
Table 2. The incidence rates of schizophrenia, PUFA consumption, and economic development across different countries.
Country CodeCountryIncidence of SchizophreniaPUFA Consumption aEconomic Variable
CoveragePeriodIncidence (Per 100,000)AA (mg/Day)DHA (mg/Day)ω-6 LCPUFA (mg/Day)ω-3 LCPUFA (mg/Day)GDP Per Capita b
BRABrazilSubnational1965–199447.2 [38]71.9765.35106.99158.03431.3
CANCanadaSubnational1995–199833.0 [39]231.99180.74356.23443.358590.8
CHLChileNational2004–201713.4 [40]113.12126.64164.82255.503076.6
CRICosta RicaNational1979–198148.2 [41]59.0049.2984.85129.73355.0
HRVCroatiaNational1980–198526.0 [42]100.1454.09165.14126.674852.2
DNKDenmarkNational1995–20089.2 [43]201.39172.04352.87366.5613,333.0
FINFinlandNational200337.8 [44]179.83283.62291.77592.4410,497.5
DEUGermanyNational1987–198919.0 [45]187.32131.59325.37314.242761.2
ISLIcelandNational196727.0 [46]297.87561.37378.271278.491418.1
INDIndiaSubnational1987–198858.5 [47]14.7029.2620.4165.7898.1
IRLIrelandSubnational1974–197727.5 [48]185.71120.95285.16355.25739.3
ITAItalySubnational2008–20119.6 [49]220.95214.39338.11468.6418,728.4
JPNJapanSubnational1979–198021.0 [50]89.70384.11120.07682.04563.6
NLDNetherlandsSubnational1981–199722.3 [51]164.43123.68267.52297.93662.2
NZLNew ZealandNational20146.3 [52]294.28242.54410.68630.0717,400.4
NORNorwaySubnational1982–19835.5 [53]154.97318.68232.38658.061721.4
KORRepublic of KoreaNational2008–201750.0 [54]137.70339.12198.58978.2619057.6
RUSRussian FederationNational20158.0 [55]162.29194.33248.96421.682665.8
ESPSpainSubnational1989–199019.0 [56]157.49263.78234.63525.391145.0
SWESwedenSubnational1997–200613.0 [57]188.92229.14305.55503.1414,272.5
CHESwitzerlandSubnational1977–200511.8 [58]188.47137.66308.20350.498953.3
UKRUkraineSubnational1986–199711.0 [59]149.60114.19248.40261.741417.9
GBRUnited KingdomSubnational1998–20054.1 [60]201.51165.01311.62403.327937.9
USAUnited StatesSubnational2007–20137.5 [61]274.59186.38406.67444.3723,478.6
a Each PUFA consumption is the average value of the consumption for the 20 years before the study period. b GDP per capita is the average value of the GDP per capita for the 20 years before the study period. PUFA—polyunsaturated fatty acid; AA—arachidonic acid; DHA—docosahexaenoic acid; LCPUFA—long-chain PUFA; and GDP—gross domestic product.
Table
Table 3. The instruments of eight PUFAs and schizophrenia in MR analyses.
Table 3. The instruments of eight PUFAs and schizophrenia in MR analyses.
Effect Size Estimates for PUFAs aEffect Size Estimates for Schizophrenia b
PUFASNPChrNearby GeneEffect AlleleF StatisticsβSEpln (OR)SEp
Linoleic acid
(LA, 18:2 ω-6)		rs1074011810JMJD1CG32.910.2480.0438.08 × 10−90.00930.01100.399
rs17454711FADS1C1218.481.4740.0424.98 × 10−2740.02050.01140.072
rs1696695216NTAN1, PDXDC1G62.960.3510.0441.23 × 10−15−0.03940.01170.001
γ-linolenic acid
(GLA, 18:3 ω-6)		rs17454711FADS1T253.260.0160.0012.29 × 10−72−0.02050.01140.072
rs1696695216NTAN1, PDXDC1G36.810.0060.0015.0 × 10−11−0.03940.01170.001
Dihomo-γ-linolenic acid
(DGLA, 20:3 ω-6)		rs17454711FADS1T1211.88−0.3500.0102.63 × 10−151−0.02050.01140.072
rs1696695216NTAN1, PDXDC1G119.700.2200.0207.55 × 10−65−0.03940.01170.001
Arachidonic acid
(AA, 20:4 ω-6)		rs17454711FADS1T4526.151.6910.0253.00 × 10−971−0.02050.01140.072
rs1696695216NTAN1, PDXDC1G40.770.1990.0312.43 × 10−10−0.03940.01170.001
α-linolenic acid
(ALA, 18:3 ω-3)		rs17454711FADS1C256.000.0160.0013.47 × 10−640.02050.01140.072
Eicosapentaenoic acid
(EPA, 20:5ω3)		rs37987136ELOVL2C49.000.0350.0051.93 × 10−120.00810.01100.459
rs17453811FADS1G275.560.0830.0055.37 × 10−58−0.02110.01170.072
Docosapentaenoic acid
(DPA, 22:5 ω-3)		rs7800942GCKRT32.110.0170.0039.04 × 10−90.00070.01110.948
rs37343986ELOVL2C177.780.0400.0039.61 × 10−440.00880.01100.424
rs17454711FADS1T625.000.0750.0033.79 × 10−154−0.02050.01140.072
Docosahexaenoic acid
(DHA, 22:6 ω-3)		rs22362126ELOVL2G65.150.1130.0141.26 × 10−15−0.00850.01100.439
a Summary statistic for GWAS of PUFA from “Genome-wide association study of plasma N6 polyunsaturated fatty acids within the cohorts for heart and aging research in genomic epidemiology consortium” [23], “Genetic loci associated with plasma phospholipid n-3 fatty acids: a meta-analysis of genome-wide association studies from the CHARGE Consortium” [24]. b Summary statistics for GWAS of schizophrenia from “Biological insights from 108 schizophrenia-associated genetic loci” [29]. PUFA—polyunsaturated fatty acid; SNP—single-nucleotide polymorphism; Chr—chromosome; SE—standard error; and OR—odds ratio.
Table
Table 4. The multivariable MR analyses for ω-6 and ω-3 PUFAs with risk for schizophrenia.
Table 4. The multivariable MR analyses for ω-6 and ω-3 PUFAs with risk for schizophrenia.
Number of SNP Instruments (Related to PUFA Group) aNumber of SNPs in AnalysisF StatisticResult in IVW Method a,b
OR (95% CI)p
Arachidonic acid
(AA, 20:4 ω-6)		3885.330.981 (0.976–0.986)0.009
Eicosapentaenoic acid
(EPA, 20:5 ω-3)		6830.671.359 (0.702–2.633)0.090
Linoleic acid
(LA, 18:2 ω-6)		386.461.005 (0.957–1.055)0.805
α-linolenic acid
(ALA, 18:3 ω-3)		686.023.747 (0.019–728.509)0.512
a Summary statistic for GWAS of PUFA from “Genome-wide association study of plasma N6 polyunsaturated fatty acids within the cohorts for heart and aging research in genomic epidemiology consortium” [23], “Genetic loci associated with plasma phospholipid n-3 fatty acids: a meta-analysis of genome-wide association studies from the CHARGE Consortium” [24]. b Summary statistics for GWAS of schizophrenia from “Biological insights from 108 schizophrenia-associated genetic loci” [29]. PUFA–polyunsaturated fatty acid; SNP–single-nucleotide polymorphism; Chr–chromosome; SE–standard error; and OR–odds ratio.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gao, Y.; Hu, X.; Wang, D.; Jiang, J.; Li, M.; Qing, Y.; Yang, X.; Zhang, J.; Zhang, Y.; Wan, C. Association between Arachidonic Acid and the Risk of Schizophrenia: A Cross-National Study and Mendelian Randomization Analysis. Nutrients 2023, 15, 1195. https://doi.org/10.3390/nu15051195

AMA Style

Gao Y, Hu X, Wang D, Jiang J, Li M, Qing Y, Yang X, Zhang J, Zhang Y, Wan C. Association between Arachidonic Acid and the Risk of Schizophrenia: A Cross-National Study and Mendelian Randomization Analysis. Nutrients. 2023; 15(5):1195. https://doi.org/10.3390/nu15051195

Chicago/Turabian Style

Gao, Yan, Xiaowen Hu, Dandan Wang, Jie Jiang, Minghui Li, Ying Qing, Xuhan Yang, Juan Zhang, Yue Zhang, and Chunling Wan. 2023. "Association between Arachidonic Acid and the Risk of Schizophrenia: A Cross-National Study and Mendelian Randomization Analysis" Nutrients 15, no. 5: 1195. https://doi.org/10.3390/nu15051195

APA Style

Gao, Y., Hu, X., Wang, D., Jiang, J., Li, M., Qing, Y., Yang, X., Zhang, J., Zhang, Y., & Wan, C. (2023). Association between Arachidonic Acid and the Risk of Schizophrenia: A Cross-National Study and Mendelian Randomization Analysis. Nutrients, 15(5), 1195. https://doi.org/10.3390/nu15051195

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop