Multi-omic modeling of antidepressant response implicates dynamic immune and inflammatory changes in individuals who respond to treatment

BACKGROUND

Major depressive disorder (MDD) is a leading cause of disability worldwide, and is commonly treated with antidepressant drugs (AD). Although effective, many patients fail to respond to AD treatment, and accordingly identifying factors that can predict AD response would greatly improve treatment outcomes. In this study, we developed a machine learning tool to integrate multi-omic datasets (gene expression, DNA methylation, and genotyping) to identify biomarker profiles associated with AD response in a cohort of individuals with MDD.

MATERIALS AND METHODS

Individuals with MDD (N = 111) were treated for 8 weeks with antidepressants and were separated into responders and non-responders based on the Montgomery-Åsberg Depression Rating Scale (MADRS). Using peripheral blood samples, we performed RNA-sequencing, assessed DNA methylation using the Illumina EPIC array, and performed genotyping using the Illumina PsychArray. To address this rich multi-omic dataset with high dimensional features, we developed integrative Geneset-Embedded non-negative Matrix factorization (iGEM), a non-negative matrix factorization (NMF) based model, supplemented with auxiliary information regarding gene sets and gene-methylation relationships. In particular, we factorize the subjects by features (i.e., gene expression or DNA methylation) into subjects-by-factors and factors-by-features. We define the factors as the meta-phenotypes as they represent integrated composite scores of the molecular measurements for each subject.

RESULTS

Using our model, we identified a number of meta-phenotypes which were related to AD response. By integrating geneset information into the model, we were able to relate these meta-phenotypes to biological processes, including a meta-phenotype related to immune and inflammatory functions as well as other genes related to depression or AD response. The meta-phenotype identified several genes including immune interleukin 1 receptor like 1 (IL1RL1) and interleukin 5 receptor (IL5) subunit alpha (IL5RA), AKT/PIK3 pathway related phosphoinositide-3-kinase regulatory subunit 6 (PIK3R6), and sphingomyelin phosphodiesterase 3 (SMPD3), which has been identified as a target of AD treatment.

CONCLUSIONS

The derived meta-phenotypes and associated biological functions represent both biomarkers to predict response, as well as potential new treatment targets. Our method is applicable to other diseases with multi-omic data, and the software is open source and available on Github (https://github.com/li-lab-mcgill/iGEM).

ERK1/2 activation in human taste bud cells regulates fatty acid signaling and gustatory perception of fat in mice and humans

Obesity is a major public health problem. An in-depth knowledge of the molecular mechanisms of oro-sensory detection of dietary lipids may help fight it. Humans and rodents can detect fatty acids via lipido-receptors, such as CD36 and GPR120. We studied the implication of the MAPK pathways, in particular, ERK1/2, in the gustatory detection of fatty acids. Linoleic acid, a dietary fatty acid, induced via CD36 the phosphorylation of MEK1/2-ERK1/2-ETS-like transcription factor-1 cascade, which requires Fyn-Src kinase and lipid rafts in human taste bud cells (TBCs). ERK1/2 cascade was activated by Ca2+ signaling via opening of the calcium-homeostasis modulator-1 (CALHM1) channel. Furthermore, fatty acid-evoked Ca2+ signaling and ERK1/2 phosphorylation were decreased in both human TBCs after small interfering RNA knockdown of CALHM1 channel and in TBCs from Calhm1-/- mice. Targeted knockdown of ERK1/2 by small interfering RNA or PD0325901 (MEK1/2 inhibitor) in the tongue and genetic ablation of Erk1 or Calhm1 genes impaired preference for dietary fat in mice. Lingual inhibition of ERK1/2 in healthy volunteers also decreased orogustatory sensitivity for linoleic acid. Our data demonstrate that ERK1/2-MAPK cascade is regulated by the opening of CALHM1 Ca2+ channel in TBCs to modulate orogustatory detection of dietary lipids in mice and humans.-Subramaniam, S., Ozdener, M. H., Abdoul-Azize, S., Saito, K., Malik, B., Maquart, G., Hashimoto, T., Marambaud, P., Aribi, M., Tordoff, M. G., Besnard, P., Khan, N. A. ERK1/2 activation in human taste bud cells regulates fatty acid signaling and gustatory perception of fat in mice and humans.