Association between circulating fatty acid metabolites and asthma risk: a two-sample bidirectional Mendelian randomization study

Fatty acids are involved in a wide range of immunological responses in humans. Supplementation of polyunsaturated fatty acids has been reported to help alleviate symptoms and airway inflammation in asthma patients, whereas the effects of fatty acids on the actual risk of asthma remain controversial. This study comprehensively investigated the causal effects of serum fatty acids on asthma risk using two-sample bidirectional Mendelian Randomization (MR) analysis.

Genetic variants strongly associated with 123 circulating fatty acid metabolites were extracted as instrumental variables, and a large GWAS data of asthma was used to test effects of the metabolites on this outcome. The inverse-variance weighted method was used for primary MR analysis. The weighted median, MR-Egger regression, MR-PRESSO, and leave-one-out analyses were utilized to evaluate heterogeneity and pleiotropy. Potential confounders were adjusted by performing multivariable MR analyses. Reverse MR analysis was also conducted to estimate the causal effect of asthma on candidate fatty acid metabolites. Further, we performed colocalization analysis to examine the pleiotropy of variants within the fatty acid desaturase 1 (FADS1) locus between the significant metabolite traits and the risk of asthma. Cis-eQTL-MR and colocalization analysis were also performed to determine the association between RNA expression of FADS1 and asthma.

Genetically instrumented higher average number of methylene groups was causally associated with a lower risk of asthma in primary MR analysis, while inversely, the higher ratio of bis-allylic groups to double bonds and the higher ratio of bis-allylic groups to total fatty acids, were associated with higher probabilities of asthma. Consistent results were obtained in multivariable MR when adjusted for potential confounders. However, these effects were completely eliminated after SNPs correlated with the FADS1 gene were excluded. The reverse MR also found no causal association. The colocalization analysis suggested that the three candidate metabolite traits and asthma likely share causal variants within the FADS1 locus. In addition, the cis-eQTL-MR and colocalization analyses demonstrated a causal association and shared causal variants between FADS1 expression and asthma.

Our study supports a negative association between several PUFA traits and the risk of asthma. However, this association is largely attributed to the influence of FADS1 polymorphisms. The results of this MR study should be carefully interpreted given the pleiotropy of SNPs associated with FADS1.