Skip to main content

ORIGINAL RESEARCH article

Front. Pharmacol., 20 December 2021
Sec. Pharmacogenetics and Pharmacogenomics
This article is part of the Research Topic Translation and Implementation of Pharmacogenomic Testing in Daily Clinical Practice: Considering Current Challenges and Future Needs View all 6 articles

Neuropsychiatric Adverse Events of Montelukast: An Analysis of Real-World Datasets and drug−gene Interaction Network

Ryogo UmetsuRyogo Umetsu1Mizuki TanakaMizuki Tanaka1Yoko NakayamaYoko Nakayama1Yamato KatoYamato Kato1Natsumi UedaNatsumi Ueda1Yuri NishibataYuri Nishibata1Shiori HasegawaShiori Hasegawa1Kiyoka MatsumotoKiyoka Matsumoto1Noriaki TakeyamaNoriaki Takeyama2Kazuhiro IguchiKazuhiro Iguchi3Hiroyuki TanakaHiroyuki Tanaka4Eiichi Hinoi,Eiichi Hinoi2,4Naoki InagakiNaoki Inagaki4Masatoshi IndenMasatoshi Inden5Yoshinori MutoYoshinori Muto6Mitsuhiro Nakamura
Mitsuhiro Nakamura1*
  • 1Laboratory of Drug Informatics, Gifu Pharmaceutical University, Gifu, Japan
  • 2Laboratory of Pharmacology, Gifu Pharmaceutical University, Gifu, Japan
  • 3Laboratory of Community Pharmacy, Gifu Pharmaceutical University, Gifu, Japan
  • 4United Graduate School of Drug Discovery and Medical Information Sciences, Gifu University, Gifu, Japan
  • 5Laboratory of Medical Therapeutics and Molecular Therapeutics, Gifu Pharmaceutical University, Gifu, Japan
  • 6Department of Functional Bioscience, Gifu University School of Medicine, Gifu, Japan

Montelukast is a selective leukotriene receptor antagonist that is widely used to treat bronchial asthma and nasal allergy. To clarify the association between montelukast and neuropsychiatric adverse events (AEs), we evaluated case reports recorded between January 2004 and December 2018 in the Food and Drug Administration Adverse Event Reporting System (FAERS). Furthermore, we elucidated the potential toxicological mechanisms of montelukast-associated neuropsychiatric AEs through functional enrichment analysis of human genes interacting with montelukast. The reporting odds ratios of suicidal ideation and depression in the system organ class of psychiatric disorders were 21.5 (95% confidence interval (CI): 20.3–22.9) and 8.2 (95% CI: 7.8–8.7), respectively. We explored 1,144 human genes that directly or indirectly interact with montelukast. The molecular complex detection (MCODE) plug-in of Cytoscape detected 14 clusters. Functional analysis indicated that several genes were significantly enriched in the biological processes of “neuroactive ligand–receptor interaction.” “Mood disorders” and “major depressive disorder” were significant disease terms related to montelukast. Our retrospective analysis based on the FAERS demonstrated a significant association between montelukast and neuropsychiatric AEs. Functional enrichment analysis of montelukast-associated genes related to neuropsychiatric symptoms warrant further research on the underlying pharmacological mechanisms.

Introduction

Montelukast is a selective leukotriene receptor antagonist used to treat bronchial asthma and nasal allergy. Although montelukast is generally well tolerated, several clinical trials and post-marketing studies have reported serious neuropsychiatric adverse events (AEs) (Philip et al., 2009a; Philip et al., 2009b; Kelsay, 2009; Calapai et al., 2014; Haarman et al., 2017). The potential association between montelukast and suicidal behavior has previously been demonstrated based on the results of a literature search of MEDLINE, EMBASE, International Pharmaceutical Abstracts, and the Food and Drug Administration (FDA) adverse event reporting system (FAERS) (Schumock et al., 2011). This information has led the FDA to issue multiple warnings concerning an increased risk of neuropsychiatric AEs after taking montelukast and other leukotriene antagonists, including aggressive behavior, anxiety, depression, abnormal dreams, excitement, hallucinations, insomnia, irritability, and potential suicidality (Marchand et al., 2013; Perona et al., 2016; Merck & Co, 2020). On March 4, 2020, the United States FDA issued a safety announcement regarding the necessity of boxed warnings about serious neuropsychiatric AEs for montelukast (Singulair) (Food and Drug Administration 2020).

The FAERS is a spontaneous reporting system (SRS) involving reports of AEs in a real-world setting that are voluntarily submitted by healthcare professionals, pharmaceutical companies, and patients. The FAERS database is publicly available, can be downloaded from the FDA website (http://www.fda.gov), and is used in the post-marketing safety assessments of approved drugs. The objective of this study was to evaluate the association between neuropsychiatric AEs using well-established pharmacovigilance indices such as reporting odds ratio (ROR).

The pharmacological mechanisms causing neuropsychiatric alterations are currently unclear (Khalid et al., 2017). Most drugs act via interactions with several proteins encoded by different genes. An analysis of drug–gene interactions improved our understanding of drug toxicity (Ludovini et al., 2016). In recent years, integrated analysis using FAERS data and drug–gene interaction analysis data has been proposed as a method to expand our knowledge of AEs (Wu et al., 2016; Lin et al., 2017; Tanaka et al., 2021). To better understand the toxicological mechanisms underlying montelukast-associated neuropsychiatric AEs, we extracted a data set of human genes interacting with montelukast from public databases and constructed a drug–gene interaction network. Functional enrichment analysis of these genes was performed to elucidate the potential toxicological mechanisms of montelukast-associated neuropsychiatric AEs.

Methods

Data Source

Data from April 2004 to December 2018 were extracted from the FAERS database on the FDA website. The informatic structure of the FAERS database is based on the international safety reporting guidelines issued by the International Council on Harmonization (ICH), known as ICH E2B guidelines (U. S. Department of Health and Human Services, 2014). We integrated our database from the FAERS dataset using FileMaker Pro Advanced software (FileMaker, Inc., Santa Clara, CA, United States), according to the ASCII Entity Relationship Diagram, which is publicly available from the FDA website (https://www.fda.gov).

Following the FDA’s recommendation, we excluded duplicate reports of the same patient from different reporting sources from the analysis and extracted reports. Drugs in FAERS are classified into four categories: primary suspect drug (PS), secondary suspect drug (SS), concomitant (C), and interacting (I), according to the anticipated degree of involvement for AEs. Only reports with the PS drug code were included in this analysis.

Definition of Adverse Events

AEs were coded with terms found in the Medical Dictionary for Regulatory Activities (MedDRA, https://www.meddra.org), which is the dictionary for terminology used in the FAERS database. This study relied on the definitions provided by MedDRA version 21.0. To evaluate montelukast-associated AEs, we utilized the system organ classes (SOCs) of “psychiatric disorders,” “general disorders and administration site conditions,” “nervous system disorders,” “respiratory, thoracic and mediastinal disorders,” and “gastrointestinal disorders” (Table 1). The preferred terms (PTs) related to each SOC are summarized in Table 1.

TABLE 1
www.frontiersin.org

TABLE 1. Number of reports and reporting odds ratio related to montelukast in the FAERS (January 2004−November 2018).

Signal Detection

We used the ROR to analyze the association between montelukast and AEs. The ROR is the ratio of the odds of reporting an AE relative to all other AEs associated with the drug of interest compared with the reporting odds for all other drugs in the FAERS database (Poluzzi et al., 2012). ROR is calculated based on the two-by-two contingency table. RORs are expressed as point estimates with 95% confidence interval (CI). The signal was considered positive when the lower limit of 95% CI was >1 and the number of reports was ≥2 (Poluzzi et al., 2012).

Drug−gene Interaction Network

The drug−gene interaction network was constructed on the basis of drug−gene and gene−gene interactions. Montelukast-associated genes were retrieved from DGIdb (drug–gene interaction database, https://www.dgidb.org), DSigDB (drug signatures database, http://dsigdb.tanlab.org), and STITCH (https://stitch.embl.de). The indirectly associated genes were retrieved from iRefIndex 15.0 (“9606. mitab,” https://irefindex.vib.be) (Razick et al., 2008). Molecular complex detection (MCODE) is an approach for detecting highly interconnected regions in protein−protein interaction networks (Bader and Hogue, 2003). The clusters likely to be involved in common biological function were investigated using the MCODE plug-in (version 1.5.1) of Cytoscape version 3.7 (http://cytoscape.org). This plug-in was utilized to choose hub modules of the gene−gene interaction network in Cytoscape with a degree cutoff = 2, node score cutoff = 0.2, k-core = 2, and Max. Depth from seed = 100 as the criteria. Next, we used “clusterProfiler (version 1.4.0),” an R package, to perform functional analysis and visualization of functional profiles for genes and gene clusters. We used the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis to explore the biological significance. KEGG enrichment analysis was performed using clusterProfiler with organism = “hsa,” pvalueCutoff = 0.05, pAdjustMethod = “BH,” and qvalueCutoff = 0.1. The thresholds in the KEGG enrichment analysis are pvalueCutoff = 0.05 and qvalueCutoff = 0.1. The default thresholds in the KEGG enrichment analysis are pvalueCutoff = 0.05 and qvalueCutoff = 0.2. Lin et al. applied pvalueCutoff = 0.05 and qvalueCutoff (not listed) (Lin et al., 2017). We could not find a gold standard for the thresholds. Finally, disease enrichment analysis based on DisGeNET was performed using the function enrichDGN in the R package named DOSE: Disease Ontology Semantic and Enrichment analysis (version 3.2).

Results

The FAERS database contains 11,527,470 reports from January 2004 to December 2018. After excluding duplicates according to the FDA recommendations, 9,702,166 were analyzed. The RORs of suicidal ideation, suicide attempts, and depression in the SOC of psychiatric disorders were 21.5 (95% CI: 20.3–22.9), 9.5 (95% CI: 8.5–10.5), and 8.2 (95% CI: 7.8–8.7), respectively.

We primarily searched DGIdb (drug–gene interaction database, https://www.dgidb.org), DSigDB (drug signatures database, https://dsigdb.tanlab.org), and STITCH (https://stitch.embl.de) and retrieved 26 genes (ABCC1, AHR, ALOX5, ATAD5, ATG4B, CCL11, CYP2C8, CYSLTR1, CYSLTR2, IL13, IL4, IL5, KDM4A, LTA4H, LTB4R, LTB4R2, LTC4S, PLA2G1B, POLH, POLI, POLK, PPP1CA, S1PR1, S1PR3, S1PR4, SLCO2B) that interact with montelukast directly. All the genes of “9606. mitab” from iRefIndex 15.0 (https://irefindex.vib.be) were integrated into a network with 20,877 nodes and 429,350 edges. The genes that directly or indirectly interact with the above 26 genes were integrated into a network with 1,144 nodes and 35,384 edges. MCODE plug-in found 14 clusters (Table 2). To translate the network into biological insights, we further performed functional enrichment analysis using the KEGG pathways. The top four clusters stratified by biological process are shown in Figure 1 (Supplemental Figures 1A–D).

TABLE 2
www.frontiersin.org

TABLE 2. Clusters of networks analyzed by MCODE.

FIGURE 1
www.frontiersin.org

FIGURE 1. Top four gene interaction networks based on the MCODE plug-in of Cytoscape.

The genes interacting with montelukast were enriched in a number of gene sets involved in “neuroactive ligand−receptor interaction” and “chemokine signaling pathway” in cluster 1, and “neuroactive ligand−receptor interaction” and “calcium signaling pathway” in cluster 2 (Figure 2).

FIGURE 2
www.frontiersin.org

FIGURE 2. Dot plot of KEGG functional enrichment analysis.

Furthermore, data retrieved from DisGeNET was used to characterize diseases associated with montelukast. We found significant enrichment in genes involved in the following diseases related to montelukast (Figure 3): “pneumonia” (adjusted p-value = 2.85 × 10–15) and “respiratory syncytial virus infections” (adjusted p-value = 1.36 × 10–17) in cluster 1 and “mood disorders” (adjusted p-value = 1.21 × 10–12) and “major depressive disorder” (adjusted p-value = 1.83 × 10–7) in cluster 2.

FIGURE 3
www.frontiersin.org

FIGURE 3. Dot plot of DisGeNET disease enrichment analysis.

Discussion

We elucidated the AE profile of montelukast using the SOC of psychiatric disorders associated with the drug in the FAERS database. The lower limits of the 95% CIs of RORs related to the SOC of “psychiatric disorders” were more than 1, and the signal was significantly detected. Neuropsychiatric AEs have been found to account for the most important costs associated with comorbidity in asthma and have a negative impact on the patients’ quality of life (Schumock et al., 2011; Chen et al., 2016; Khalid et al., 2017), although most symptoms improve upon stopping montelukast therapy.

To better understand the toxicological mechanisms of montelukast-associated neuropsychiatric AEs, we curated drug−gene interactions from public databases. A total of 1,144 human genes interacting with montelukast were investigated. Some of these genes that were highly enriched in DisGeNET were related to “mood disorders” and “major depressive disorder” (Figure 3). HCRT (hypocretin neuropeptide precursor), HTR2A (5-hydroxytryptamine receptor 2A), and KALRN (kalirin RhoGEF kinase) genes were enriched in the modules “mood disorders” and “major depressive disorder” in cluster 2 (Table 2). HCRT encodes hypocretin, a hypothalamic neuropeptide precursor protein that gives rise to two mature neuropeptides, orexin A and orexin B, via proteolytic processing. The hypothalamic-pituitary-adrenal (HPA) axis plays an important role in the network mediated by stress-related neurotransmitters and have been proposed to affect depression (Bao et al., 2012) and suicide (Turecki et al., 2012). Hypocretins produced in the hypothalamus (Hunt et al., 2015) have functional interactions with the HPA axis and regulate sleep, feeding, energy balance, sexual behavior, and stress response, which are affected in depression (Nollet and Leman, 2013). HTR2A encodes 5-HT2A receptors, which are associated with major depressive disorder, schizophrenia, and suicidality (Niculescu et al., 2017). KALRN is a protein-coding gene that has been associated with stroke (Krug et al., 2010), coronary heart disease (Wang et al., 2007; Beręsewicz et al., 2008; Krug et al., 2010), schizophrenia (Hill et al., 2006; Hayashi-Takagi et al., 2010; Bradshaw and Porteous, 2012), and adult attention-deficit/hyperactivity disorder (Lesch et al., 2008). These findings suggest that montelukast could increase the risk of “psychiatric disorders.”

The common AEs caused by montelukast are upper airway infections, anaphylaxis, nausea, vomiting, diarrhea, elevated levels of liver enzymes, agitation, anxiety, depression, sleep disturbance, and eosinophilic granulomatosis with polyangiitis (EGPA), also known as Churg–Strauss syndrome (Calapai et al., 2014; Merck & Co, 2020). The presence of the ROR signal of EGPA in our study indicates the association of EGPA with the use of montelukast. However, the hypothesis that EGPA is not attributed to montelukast but to the reduction in the dose of glucocorticoid used in combination with montelukast has recently been accepted (Bibby et al., 2010). We considered that the value of the ROR related to EGPA was only apparently high.

Many studies supporting an association between leukotriene-modifying agents including montelukast and suicidality are primarily based on reviews of individual safety reports in AE databases which are subject to reporting bias and confounding factors. On the contrary, case-control and cohort studies, and clinical trials do not support an association between the two. Ecological studies have demonstrated a lack of positive association between leukotriene-modifying agents and suicidality at the population level (Khalid et al., 2017). Although our study is based on the FAERS database, it also has some limitations that are worth mentioning. As the FAERS is an SRS, it has several limitations including biases (under-reporting, over-reporting, missing data, and comorbidities), a lack of detailed information about the patients, and the exclusion of healthy individuals as a reference group. Therefore, ROR cannot be used for assessing true risks and ranking AEs. The risk of suicidal behavior increases among patients with respiratory diseases such as asthma (Schumock et al., 2011; Khalid et al., 2017). It has been elucidated how co-morbidities render FAERS data difficult to interpret (compared with controlled study data). Another limitation is that the SOCs in the FAERS data analysis and “mood disorders” and “major depressive disorder” in drug–gene analyses might not exactly represent the same clinical outcomes. Therefore, our results from the FAERS database must be interpreted considering these limitations. Further epidemiological studies using a large number of patients and well-controlled trials are required to confirm the safety risks of montelukast. When prescribing montelukast, clinicians should carefully monitor patients who may be at elevated risk for suicidal ideation or depression, according to the boxed warnings.

Some limitations of our functional enrichment analysis should also be noted. Our results do not offer any hard evidence regarding the potential mechanisms of montelukast-associated neuropsychiatric AEs. For now, the drug–gene interactions investigated have not been validated in any experimental model or in vitro and in vivo experiments because of our currently limited knowledge about disease-associated proteins and their interactions. Therefore, the association between montelukast and genes should be confirmed experimentally. Furthermore, we identified a list of 26 proteins interacting with montelukast. However, not all genes from databases like DGIdb are direct pharmacological targets; many may be indirectly affected by drugs. Our analysis uses a protein–protein interaction network to map these 26 affected proteins to larger networks and demonstrate that the networks are enriched in genes pertaining to mood disorders. We seized on the genes HCRT, KALRN, and HTR2A to substantiate the connection to mood disorders. As these genes are not direct targets of montelukast, this approach must be validated by showing that a given “hotspot” in such a protein–protein interaction network distinguishes drugs that cause “mood disorders” from drugs that do not. The modular assembly of drug safety subnetworks (MADSS) algorithm may be suitable for solving this problem (Lorberbaum et al., 2015).

Conclusion

Our retrospective analysis demonstrated a significant association between montelukast and neuropsychiatric AEs. The genes that were thought to be associated with neuropsychiatric symptoms due to their interaction with montelukast were found to be significantly enriched in functional categories of psychiatric disease, which necessitates future pharmacological research.

Data Availability Statement

All datasets generated for this study are included in the article/Supplementary Material.

Author Contributions

RU, MT, YM, and MN contributed to the overall concept and design of the study. RU, MT, and MN wrote the main manuscript. RU, YK, and MN created the relational FAERS database system used in the work. YN, NU, YN, SH, KM, and NT carried out data extraction and statistical analysis. YM, MT, and MN contributed important intellectual content in the drug–gene interaction analysis. KI, HT, EH, NI, and MI revised the article critically for important intellectual content. All authors have reviewed the manuscript.

Funding

This research was partially supported by JSPS KAKENHI Grant Number, 17K08452, 20K10408, and 21K06646.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2021.764279/full#supplementary-material

References

Bader, G. D., and Hogue, C. W. (2003). An Automated Method for Finding Molecular Complexes in Large Protein Interaction Networks. BMC Bioinformatics 4, 2. doi:10.1186/1471-2105-4-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Bao, A. M., Ruhé, H. G., Gao, S. F., and Swaab, D. F. (2012). Neurotransmitters and Neuropeptides in Depression. Handb. Clin. Neurol. 106, 107–136. doi:10.1016/B978-0-444-52002-9.00008-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Beręsewicz, M., Kowalczyk, J. E., and Zabłocka, B. (2008). Kalirin-7, a Protein Enriched in Postsynaptic Density, Is Involved in Ischemic Signal Transduction. Neurochem. Res. 33, 1789–1794. doi:10.1007/s11064-008-9631-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Bibby, S., Healy, B., Steele, R., Kumareswaran, K., Nelson, H., and Beasley, R. (2010). Association between Leukotriene Receptor Antagonist Therapy and Churg-Strauss Syndrome: An Analysis of the FDA AERS Database. Thorax 65, 132–138. doi:10.1136/thx.2009.120972

PubMed Abstract | CrossRef Full Text | Google Scholar

Bradshaw, N. J., and Porteous, D. J. (2012). DISC1-binding Proteins in Neural Development, Signalling and Schizophrenia. Neuropharmacology 62, 1230–1241. doi:10.1016/j.neuropharm.2010.12.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Calapai, G., Casciaro, M., Miroddi, M., Calapai, F., Navarra, M., and Gangemi, S. (2014). Montelukast-induced Adverse Drug Reactions: A Review of Case Reports in the Literature. Pharmacology 94, 60–70. doi:10.1159/000366164

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, W., Lynd, L. D., FitzGerald, J. M., Marra, C. A., Balshaw, R., To, T., et al. (2016). Excess Medical Costs in Patients with Asthma and the Role of Comorbidity. Eur. Respir. J. 48, 1584–1592. doi:10.1183/13993003.01141-2016

PubMed Abstract | CrossRef Full Text | Google Scholar

Food and Drug Administration (2020). FDA Requires Boxed Warning about Serious Mental Health Side Effects for Asthma and Allergy Drug Montelukast (Singulair); Advises Restricting Use for Allergic Rhinitis: Risks May Include Suicidal Thoughts or Actions. Available at: https://www.fda.gov/drugs/drug-safety-and-availability/fda-requires-boxed-warning-about-serious-mental-health-side-effects-asthma-and-allergy-drug (Accessed August 23, 2021).

Google Scholar

Haarman, M. G., van Hunsel, F., and de Vries, T. W. (2017). Adverse Drug Reactions of Montelukast in Children and Adults. Pharmacol. Res. Perspect. 5, e00341. doi:10.1002/prp2.341

CrossRef Full Text | Google Scholar

Hayashi-Takagi, A., Takaki, M., Graziane, N., Seshadri, S., Murdoch, H., Dunlop, A. J., et al. (2010). Disrupted-in-Schizophrenia 1 (DISC1) Regulates Spines of the Glutamate Synapse via Rac1. Nat. Neurosci. 13, 327–332. doi:10.1038/nn.2487

PubMed Abstract | CrossRef Full Text | Google Scholar

Hill, J. J., Hashimoto, T., and Lewis, D. A. (2006). Molecular Mechanisms Contributing to Dendritic Spine Alterations in the Prefrontal Cortex of Subjects with Schizophrenia. Mol. Psychiatry 11, 557–566. doi:10.1038/sj.mp.4001792

PubMed Abstract | CrossRef Full Text | Google Scholar

Hunt, N. J., Rodriguez, M. L., Waters, K. A., and Machaalani, R. (2015). Changes in Orexin (Hypocretin) Neuronal Expression with normal Aging in the Human Hypothalamus. Neurobiol. Aging 36, 292–300. doi:10.1016/j.neurobiolaging.2014.08.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Kelsay, K. (2009). Assessing Risk: Data from Montelukast Clinical Trials. J. Allergy Clin. Immunol. 124, 697–698. doi:10.1016/j.jaci.2009.09.001

CrossRef Full Text | Google Scholar

Khalid, F., Aftab, A., and Khatri, S. (2017). The Association between Leukotriene-Modifying Agents and Suicidality: A Review of Literature. Psychosomatics 59, 19–27. doi:10.1016/j.psym.2017.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Krug, T., Manso, H., Gouveia, L., Sobral, J., Xavier, J. M., Albergaria, I., et al. (2010). Kalirin: A Novel Genetic Risk Factor for Ischemic Stroke. Hum. Genet. 127, 513–523. doi:10.1007/s00439-010-0790-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Lesch, K. P., Timmesfeld, N., Renner, T. J., Halperin, R., Röser, C., Nguyen, T. T., et al. (2008). Molecular Genetics of Adult ADHD: Converging Evidence from Genome-wide Association and Extended Pedigree Linkage Studies. J. Neural Transm. (Vienna) 115, 1573–1585. doi:10.1007/s00702-008-0119-3

CrossRef Full Text | Google Scholar

Lin, Y., He, S., Feng, R., Xu, Z., Chen, W., Huang, Z., et al. (2017). Digoxin-induced Anemia Among Patients with Atrial Fibrillation and Heart Failure: Clinical Data Analysis and Drug-Gene Interaction Network. Oncotarget 8, 57003–57011. doi:10.18632/oncotarget.18504

PubMed Abstract | CrossRef Full Text | Google Scholar

Lorberbaum, T., Nasir, M., Keiser, M. J., Vilar, S., Hripcsak, G., and Tatonetti, N. P. (2015). Systems Pharmacology Augments Drug Safety Surveillance. Clin. Pharmacol. Ther. 97, 151–158. doi:10.1002/cpt.2

PubMed Abstract | CrossRef Full Text | Google Scholar

Ludovini, V., Bianconi, F., Siggillino, A., Piobbico, D., Vannucci, J., Metro, G., et al. (2016). Gene Identification for Risk of Relapse in Stage I Lung Adenocarcinoma Patients: a Combined Methodology of Gene Expression Profiling and Computational Gene Network Analysis. Oncotarget 7, 30561–30574. doi:10.18632/oncotarget.8723

PubMed Abstract | CrossRef Full Text | Google Scholar

Marchand, M.-S., Jonville-Béra, A.-P., and Autret-Leca, E. (2013). Troubles psychiatriques associés au montélukast : données de la base nationale de pharmacovigilance. Arch. de Pédiatrie 20, 269–273. doi:10.1016/j.arcped.2012.12.006

CrossRef Full Text | Google Scholar

Merck & Co (2020). HIGHLIGHTS of PRESCRIBING INFORMATION: Singulair. Available at: https://www.merck.com/product/usa/pi_circulars/s/singulair/singulair_pi.pdf (Accessed August 23, 2021).

Google Scholar

Niculescu, A. B., Le-Niculescu, H., Levey, D. F., Phalen, P. L., Dainton, H. L., Roseberry, K., et al. (2017). Precision Medicine for Suicidality: From Universality to Subtypes and Personalization. Mol. Psychiatry 22, 1250–1273. doi:10.1038/mp.2017.128

PubMed Abstract | CrossRef Full Text | Google Scholar

Nollet, M., and Leman, S. (2013). Role of Orexin in the Pathophysiology of Depression: Potential for Pharmacological Intervention. CNS Drugs 27, 411–422. doi:10.1007/s40263-013-0064-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Perona, A. A., García-Sáiz, M., and Sanz Álvarez, E. (2016). Psychiatric Disorders and Montelukast in Children: A Disproportionality Analysis of the VigiBase. Drug Saf. 39, 69–78. doi:10.1007/s40264-015-0360-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Philip, G., Hustad, C., Noonan, G., Malice, M. P., Ezekowitz, A., Reiss, T. F., et al. (2009b). Reports of Suicidality in Clinical Trials of Montelukast. J. Allergy Clin. Immunol. 124, 691–e6. doi:10.1016/j.jaci.2009.08.010

CrossRef Full Text | Google Scholar

Philip, G., Hustad, C. M., Malice, M. P., Noonan, G., Ezekowitz, A., Reiss, T. F., et al. (2009a). Analysis of Behavior-Related Adverse Experiences in Clinical Trials of Montelukast. J. Allergy Clin. Immunol. 124, 699–e8. doi:10.1016/j.jaci.2009.08.011

CrossRef Full Text | Google Scholar

Poluzzi, E., Raschi, E., Piccinni, C., and De, F. (2012). Data Mining Techniques in Pharmacovigilance: Analysis of the Publicly Accessible FDA Adverse Event Reporting System (AERS). London: INTECH. 265–302. doi:10.5772/50095

CrossRef Full Text | Google Scholar

Razick, S., Magklaras, G., and Donaldson, I. M. (2008). iRefIndex: a Consolidated Protein Interaction Database with Provenance. BMC Bioinformatics 9, 405. doi:10.1186/1471-2105-9-405

PubMed Abstract | CrossRef Full Text | Google Scholar

Schumock, G. T., Lee, T. A., Joo, M. J., Valuck, R. J., Stayner, L. T., and Gibbons, R. D. (2011). Association between Leukotriene-Modifying Agents and Suicide: what Is the Evidence? Drug Saf. 34, 533–544. doi:10.2165/11587260-000000000-00000

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanaka, M., Matsumoto, K., Satake, R., Yoshida, Y., Inoue, M., Hasegawa, S., et al. (2021). Gentamicin-induced Hearing Loss: a Retrospective Study Using the Food and Drug Administration Adverse Event Reporting System and a Toxicological Study Using Drug−gene Network Analysis. Heliyon 7, e07429. doi:10.1016/j.heliyon.2021.e07429

PubMed Abstract | CrossRef Full Text | Google Scholar

Turecki, G., Ernst, C., Jollant, F., Labonté, B., and Mechawar, N. (2012). The Neurodevelopmental Origins of Suicidal Behavior. Trends Neurosci. 35, 14–23. doi:10.1016/j.tins.2011.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

U. S. Department of Health and Human Services (2014). “Food and Drug Administration, Center for Drug Evaluation and Research (CDER), and Center for Biologics Evaluation and Research (CBER),” in Guidance for Industry, E2B (R3) Electronic Transmission of Individual Case Safety Reports (ICSRs), Implementation Guide – Data Elements and Message Specification. Available at: https://www.fda.gov/media/81904/download (Accessed August 23, 2021).

Google Scholar

Wang, L., Hauser, E. R., Shah, S. H., Pericak-Vance, M. A., Haynes, C., Crosslin, D., et al. (2007). Peakwide Mapping on Chromosome 3q13 Identifies the Kalirin Gene as a Novel Candidate Gene for Coronary Artery Disease. Am. J. Hum. Genet. 80, 650–663. doi:10.1086/512981

CrossRef Full Text | Google Scholar

Wu, M., Yu, Q., and Li, Q. (2016). Differences in Reproductive Toxicology between Alopecia Drugs: an Analysis on Adverse Events Among Female and Male Cases. Oncotarget 7, 82074–82084. doi:10.18632/oncotarget.12617

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: montelukast, neuropsychiatric adverse events, food and drug administration adverse event reporting system, drug-gene intraction, protein-protein interaction, enrichment analysis

Citation: Umetsu R, Tanaka M, Nakayama Y, Kato Y, Ueda N, Nishibata Y, Hasegawa S, Matsumoto K, Takeyama N, Iguchi K, Tanaka H, Hinoi E, Inagaki N, Inden M, Muto Y and Nakamura M (2021) Neuropsychiatric Adverse Events of Montelukast: An Analysis of Real-World Datasets and drug−gene Interaction Network. Front. Pharmacol. 12:764279. doi: 10.3389/fphar.2021.764279

Received: 25 August 2021; Accepted: 22 November 2021;
Published: 20 December 2021.

Edited by:

George P Patrinos, University of Patras, Greece

Reviewed by:

Bo Liu, Sichuan University, China
Lise Aagaard, Independent researcher, Copenhagen, Denmark

Copyright © 2021 Umetsu, Tanaka, Nakayama, Kato, Ueda, Nishibata, Hasegawa, Matsumoto, Takeyama, Iguchi, Tanaka, Hinoi, Inagaki, Inden, Muto and Nakamura. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Mitsuhiro Nakamura, mnakamura@gifu-pu.ac.jp

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.