Hiccups can significantly reduce the quality of life of patients and can occur as a drug side effect. Previous report using the Drug Adverse Event Database revealed sex-specific differences in the incidence of drug-induced hiccups. However, the pathogenesis of drug-induced hiccups remains unknown, and there is limited evidence on its treatment or prevention. This study examined molecular initiating events (MIEs), which are the starting point of adverse events, to investigate the drug-induced pathways of hiccups. We extracted drugs suspected to cause hiccups using FDA Adverse Event Reporting System, a large database on adverse drug reactions. Information on drugs suspected to be associated with hiccups were extracted from the overall population and sex-specific subgroups divided. In each data table, the predicted activity values of nuclear receptors (NRs) and stress response pathways (SPs) for each drug were calculated using the toxicity predictor, a machine-learning model. Univariate analysis extracted NRs and SPs associated with hiccups, and a multivariate analysis was performed to evaluate independent factors. Transforming growth factor-beta and antioxidant response element were considered an independent factor for hiccups in the male and female subgroups, respectively. This report first examined one of the mechanisms of drug-induced hiccups and identified MIEs associated with drug-induced hiccups.

Keywords:

Subject: Medicine and Pharmacology - Pharmacology and Toxicology

1. Introduction

Hiccups are a common symptom experienced by many people and are mainly caused by myoclonus of the diaphragm [1]. Hiccups have been associated with speech, sleep, and swallowing disorders and weight loss, fatigue, and insomnia [2]. Although not life-threatening, hiccups can significantly reduce the quality of life of patients. If hiccups occur as a drug side effect, they can also impact the choice of drug treatment. Therefore, controlling hiccups that are drug side effects is important for drug therapy.

Hiccups are brief, involuntary spasms of the muscles of the diaphragm accompanied by coordinated contractions of the glottis closure muscle group [1]. The afferent and efferent pathways of the hiccup reflex arch involve the glossopharyngeal nerve (9th cranial nerve), vagus nerve (10th cranial nerve), nucleus of the solitary bundle tract, nucleus ambiguous, and peroneal nerve [3]. The exact mechanism of the hiccup reflex arch is unclear. However, neurotransmitters such as gamma-aminobutyric acid (GABA), serotonin, and dopamine have been associated with the onset of hiccups. Reflex arches are mediated by central neurotransmitters (GABA, dopamine, and serotonin) and peripheral neurotransmitters (epinephrine, norepinephrine, acetylcholine, and histamine) [4,5]. There are no reports on sex differences or physical information about common hiccups experienced by healthy individuals. However, male patients have a high incidence of intractable or persistent hiccups [6].

Several reports have described drug-induced hiccups. Lee et al. revealed a male predominance for peripheral hiccups based on a meta-analysis. However, there were no sex difference in hiccups caused by central nervous system disorders [6]. In a study about the risk factors of patient-induced hiccups, low body mass index, nausea and vomiting, and cancer drug therapy were associated with hiccups [7]. In other reports, cisplatin [8,9] and dexamethasone [10] have been suspected to affect the development of hiccups. However, no studies have analyzed the risk factors for drug-induced hiccups using data from large databases. Therefore, in a previous study, drug and patient information related to hiccups were analyzed using data from the Japanese Adverse Drug Event Report (JADER) database and the FDA Adverse Event Reporting System (FAERS) database via data mining [11,12]. Dexamethasone and several anticancer drugs were identified as independent risk factors for hiccups. We also established a visualization method for suspected drugs in the adverse drug reaction database. In a prior study using data from JADER, data on the drugs suspected to affect hiccups were extracted, and results showed that a patient’s height is a risk factor for hiccups [11]. Furthermore, in an analysis using FAERS data stratified according to sex, the drugs suspected to affect hiccups differed significantly between male and female patients. Among them, nicotine was commonly suspected to be associated with hiccups in male and female individuals [12]. Hiccups have male predominance [6]; however, their cause has not been elucidated. No studies have yet validated whether the male predominance of hiccup occurrence is attributed to drug sensitivity or neurological sex differences.

As shows in a previous study, several case reports and database studies have assessed the association between hiccups and drugs. However, to the best of our knowledge, only a few studies have investigated the mechanism of drug-induced hiccups. The small number of cases involving drug-induced hiccups and the difficulty in predicting their onset are some reasons of the causes for the lack of research. To treat and prevent drug-induced hiccups, the causative mechanisms should be identified. Therefore, the molecular chemistry related to the induction of hiccups based on the suspected drug information identified using FAERS. By handling information on suspected drugs obtained from several sources, rare cases of adverse drug reactions, such as hiccups, can be evaluated. These results may provide leads for novel research to elucidate the pathogenic mechanism of hiccups.

Numerous factors are involved in the mechanism of the effects and adverse drug reactions. From the perspective of toxicology, molecular initiating events (MIEs) in the adverse outcome pathway (AOP) represent an important concept when considering the mechanism of drug-induced adverse effects [13]. MIEs denote the first interaction between a molecule and a biomolecule or the biological system in the AOP. Their targets include nuclear receptors (NRs) and stress response pathways (SPs). NRs are intracellular proteins that regulate DNA transcription in the cell nucleus by binding to hormones and other substances. NRs bind directly to DNA and regulate gene expression, and they are associated with the development, homeostasis, and metabolism of the organism. If a ligand binds to a nuclear receptor, a conformational change occurs, which activates the receptor and subsequently regulates gene expression [14]. SPs are biochemical processes that regulate the response to various stress signals in cells. Endocrine-disrupting chemicals disrupt the endocrine system by interacting with NRs and SPs, causing various adverse developmental, reproductive, neural, and immunological effects in humans and wildlife [15]. Moreover, NRs and SPs may be involved in adverse drug events.

Drug-induced hiccups may develop hours or days after drug administration and may persist for an extended period [7,10]. Therefore, these phenomena are not controlled on a millisecond-by-millisecond basis but are thought to develop through several molecular pathways from the time of drug administration. Given these facts, we believe that searching for NRs/SRs as MIEs involved in the induction of drug-induced hiccups will help to elucidate their mechanisms. Therefore, this study aimed to identify NRs and SPs associated with drug-induced hiccups using Toxicity Predictor, a machine learning model.

2. Results

2.1. Prediction of MIEs Associated with Hiccups

There were 169 drugs with >1,000 reports and those considered significant (p ≤ 0.001 based on the Fisher’s exact test). Supplementary File S1 shows a list of 157 drugs, including the ATC names and classifications of the suspected drugs. Univariate analysis using lnROR as the binary classification value for the objective variable revealed that the significant factors were the estrogen receptor alpha with antagonist (ERα), progesterone receptor (PR), androgen receptor with antagonist (ARα), transforming growth factor-beta (TGF-β), and sonic hedgehog (Shh). Figure 1 presents a scatter plot (volcano plot) depicting the results of univariate analysis for all patient data. Table 1 depicts the activity types and the results of both univariate and multivariate analyses for each MIE.

2.2. Prediction of MIEs Associated with Hiccups in Male Patients

There were 122 drugs with >1,000 reports and those considered significant (p ≤ 0.001 based on the Fisher’s exact test) in men. Supplementary File S1 depicts a list of 118 drugs, including the ATC names and classifications of the suspected drugs. Univariate analysis using lnROR as the binary classification value for the objective variable revealed that progesterone receptor (PR), constitutive androstane receptor (CAR), Shh, androgen receptor lbd (ARlbd), and TGF-β were identified as significant factors. Figure 2 presents a volcano plot of the key suspect drugs for hiccups in male patients. Table 2 presents the activity types and the results of both univariate and multivariate analyses for each MIE.

2.3. Prediction of MIEs Associated with Hiccups in Female Patients

There were 45 drugs with >1,000 reports and those considered significant (p ≤ 0.001 based on the Fisher’s exact test) in female participants. Supplementary File S1 shows a list of 42 drugs, including the ATC names and classifications of the suspected drugs. Based on the univariate analysis using lnROR as the binary classification value for the objective variable, the significant factors were antioxidant response element (ARE), CAR, and ARa. Figure 3 presents a volcano plot of the key suspect drugs for hiccups in female patients. Table 3 depicts the activity types and the results of both univariate and multivariate analyses for each MIE.

Figure 4 presents the results of univariate analysis for the onset of hiccups and activity of each MIE during drug treatment. By taking the negative value of the ordinary logarithm of Fisher's exact test on the vertical axis and the natural logarithm of the odds ratio on the horizontal axis to create a scatter plot, the MIEs likely to be strongly associated with the development of hiccups are indicated at the top. We visually confirmed that the associated MIEs differed for men and women.

This volcano plot was created by plotting the negative logarithm of the p value (−log10p) using the Fisher’s exact test on the y-axis and the natural logarithm of the ROR (lnROR) on the x-axis. The red line on the y-axis represents p = 0.001. In female patients, drugs with >1,000 reports and those considered significant (p ≤ 0.001 using the Fisher’s exact test) are included in the volcano plot. The color of the individual points represents differences in the log of the number of reports for each drug. The gray plots indicate drugs with p > 0.001 or with < 1,000 reported cases. In this scatter plot, the signal is larger for the points (drugs) plotted in the upper right corner. The blue-to-red colors represent the number of times an adverse drug reaction was reported.

3. Discussion

This study used data from the FAERS database, an adverse drug reaction database, to identify the key suspected drugs causing hiccups. Furthermore, the association between the suspected drugs and MIEs was investigated. To the best of our knowledge, there is limited information on the mechanism of drug-induced hiccups, and no reports based on adverse events reporting systems have revealed the mechanism of drug-induced hiccups. In addition, this was the first study that aimed to examine the MIEs , especially NRs/SRs, of drug-induced hiccups using Toxicity Predictor and FAERS. The Toxicity Predictor can predict the agonist and antagonist activities of a drug in the MIE, which is the starting point of action of an adverse effect. Previous reports have described methods for evaluating the MIEs of adverse effects using the Toxicity Predictor [13]. We believe that such methods are reliable [16].

To investigate the association between each suspected drug and NRs/SPs, the objective variable was a binary variable of lnROR > 0 (positive signal) or < 0 (negative signal). A method to identify prophylactic and therapeutic agents against side effects by focusing on drugs with positive and negative signals is a method of drug repositioning and has been attracting significant attention in recent years [17,18,19]. This study examined the involvement of NRs/SPs in hiccups by comparing the MIE activities of drugs with positive and negative signals for drug-induced hiccups.

The data used in this study were classified into the overall cohort and the male and female subgroups. Drug-induced hiccups occur more frequently in men [6]. In our previous study using the adverse drug reaction database, hiccups were also more commonly in men, and the distinctions regarding the drugs suspected were noted between male and female patients [11,12]. Based on these findings, our objective was to identify the differences in NRs and SPs involved in the development of hiccups between male and female patients. Multivariate analysis identified the TGF-β agonist as an independent factor associated with hiccups in the overall dataset. In the male subgroup, TGF-β was an independent factor based on the univariate and multivariate analyses. In the female group, ARE was an independent factor in the univariate and multivariate analyses.

To the best of our knowledge, this report first explored the association between drug-induced hiccups and TGF-β. Toxicity Predictor, a toxicity prediction program, centers its predictions on the toxicity results of Tox21. To assess the cytotoxicity of a compound in Tox21, the transduction process by which TGF-β binds to the TGF-β receptor and acts on nuclear gene expression via the Smad protein is being monitored.

TGF-β is an important cytokine that maintains bodily homeostasis. It has been detected in a broad array of tissues and organs and is known for its roles in inhibiting cell proliferation and inducing cell differentiation and apoptosis in various cell types. In addition, TGF-β is involved in processes such as cell differentiation, migration, and adhesion. Moreover, it plays significant roles in diverse mechanisms including ontogeny, tissue remodeling, wound healing, inflammation, immunity, and cancer invasion and metastasis. TGF-β was identified in the cerebral spinal fluid of mice exhibiting exhaustion. Thus, its use as a fatigue marker has gained attention [20].

TGF-β, which is ubiquitous and produced in several tissues and cells, is released in an inactive, latent form (latent TGF-β) that cannot bind to receptors. It becomes activated when in the vicinity of target cells where it transforms into its active form (active TGF-β) that can bind to receptors and exert its effects. Astrocytes are known to be one of the major cells producing and releasing TGF-β in the brain where GABA receptors are located, with drug action and viral infection having been reported as facilitating factors [21,22]. Known external factors that activate TGF-β include acids, alkalis, heat, reactive oxygen species, vitamin A, vitamin D, anti-estrogens, bleomycin, and dexamethasone, among others [23]. The TGF-β family receptors have been categorized into the following three types: type I, type II, and type III. Within the signaling pathway, ligand-bound type II receptors initiate the activation of type I receptors via phosphorylation, followed by autophosphorylation. This process activates and binds Smad2 and Smad3. After phosphorylation, Smad2/3 associates with Smad4, leading to the translocation of this complex into the nucleus where it functions as a transcription factor (Figure 5). The TGF-β/Smad signaling pathway is important for regulating cell development and growth. Disruptions in this pathway have also been significantly associated with tumor development [24].

In this study, TGF-β was identified as a novel factor implicated in the induction of hiccups. Previous research has shown that drug-induced hiccups are more prevalent in men and that there are sex differences in the drugs commonly implicated [12]. Based on these findings, TGF-β was identified as an associated factor in both the overall dataset and specifically in the male subgroup. This discovery could significantly enhance our understanding of sex differences in hiccups. Moreover, hiccups are associated with GABA, an inhibitory system neuron, with the GABA derivative baclofen used for treating hiccups [25]. TGF-β promotes the growth of dopamine nerves and regulates the GABAergic nervous system. Although the exact intracellular mechanism of TGF-β is not completely understood, some studies have revealed that Erk1/2 and GSK3β might increase GABAergic neurotransmission by inhibiting the phosphorylation of gephyrin, a scaffolding protein for the GABA receptor [26]. Hence, TGF-β can be associated with the onset of hiccups via its influence on the GABAergic nervous system. Current evidence showing an association between TGF-β and drug-induced hiccups is limited. Hence, it is considered a possible MIE in the mechanistic understanding of hiccups. Ongoing research about the TGF-β signaling pathway can further clarify its role in neurotransmission.

The TGF-β superfamily interacts with receptors located on the plasma membrane. Upon binding with TGF-β, type 2 receptors undergo phosphorylation and subsequently form a complex with type 1 receptors. Activation of this receptor complex leads to the phosphorylation of Smad, an intracellular signaling molecule, which then forms a complex. It has been hypothesized that the interaction of TGF-β with GABA receptors may inhibit the phosphorylation of gephyrin, a scaffolding protein essential for GABA receptors, thereby enhancing GABAergic neurotransmission. However, the specific details of this interaction and its mechanisms remain significantly unexplored.

In the analysis of the female subgroup, ARE, one of the mechanisms associated with protection against oxidative stress, was found to be a factor associated with hiccups. The odds ratio [OR: 0.08 [0.01–0.70]) indicates a negative signal for drug-induced hiccups, which may be associated with hiccup inhibition. The Tox21 program monitors ARE activation via the Nrf2/antioxidant response signaling pathway to assess compound toxicity.

Nrf2-ARE is an important signaling pathway that regulates the expression of antioxidant enzymes. In the cytosol, Nrf2 binds to Kelch-like ECH-associated protein 1 (KEAP1) in the cytoplasm, and its activity remains low under normal physiological conditions [27]. When cells are exposed to oxidative stress, Nrf2 unbinds KEAP1 and moves into the nucleus to bind with ARE (Figure 6). This promotes the transcription and expression of a series of antioxidant enzymes, including heme oxygenase-1 (HO-1), superoxide dismutase (SOD), and thioredoxin (Trx) [28].

Currently, there are no reports on hiccups and oxidative stress. However, drugs that have been reported as important suspect drugs for hiccups (e.g., nicotine, dexamethasone, and anticancer drugs) increase oxidative stress. Nicotine causes increased oxidative stress, leading to greater neuronal apoptosis, DNA damage, reactive oxygen species, and lipid peroxides. Nicotinic acetylcholine receptors (nAChRs) have been identified in the tissues other than those in the nervous system, and their effects on nicotinic receptors have been associated with acute and chronic effects [29]. Glucocorticoid levels increase if the organism is stressed, which is accompanied by an increase in free radicals. Dexamethasone ingestion mimics the adverse effects of increased corticosterone in vivo [30]. These findings may explain why hiccups are associated with oxidative stress. The fact that the current study showed the involvement of ARE in the inhibitory side of hiccups may also support an association between hiccups and oxidative stress. We believe that ARE was an independent factor only in the female subgroup because of sex differences regarding the suspected drug for hiccups. Hence, future studies on oxidative stress and sex differences in antioxidant mechanisms should be performed. Regarding the independence of factors in this subgroup analysis, the results of multiple logistic regression analysis showed that ARE had a p value of <0.05 based on the likelihood ratio test, but not the Wald test (Table 3). As a possible cause for the uncertainty of the results, the small number of female patients used in this analysis might have affected the test. It is necessary to wait for the accumulation cases involving female patients to corroborate the current results and ensure the validity of the discussion.

The results of the current study indicated that should pathways related to NRs and SPs be associated with the induction of hiccups, the mechanism of their development may differ between men and women (Figure 4). However, given the clear sex differences in drug-induced hiccups, some sex hormone receptor involvement in the mechanism underlying the onset of drug-induced hiccups can be expected. The nuclear receptor activities we used to predict hiccups included agonist and antagonist activities of the androgen and estrogen receptors. However, no significant correlation had been observed between these activities and the induction of hiccups. The limited number of reports on hiccups compared with other side effects may explain why no relationship between sex hormone receptors and hiccups could be detected.

Another study on TGF-β and sex differences observed a difference in gene expression in the lens of male and female rats with cataract formations induced by TGF-β [31]. The mentioned report suggests that sex differences may be associated with the expression and activation of TGF-β. Given that reports on sex differences in TGF-β in humans are limited, future studies need to examine these sex differences in more detail.

Nuclear factor erythroid 2-related factor 2 (Nrf2)–ARE is an important signaling pathway that regulates the expression of antioxidant enzymes. Nrf2 binds to KEAP1 in the cytoplasm and maintains a low activity under normal physiological conditions. If cells are exposed to oxidative stress, Nrf2 unbinds KEAP1 and translocate into the nucleus to bind with ARE. This promotes the transcription and expression of a series of antioxidant enzymes, including heme oxygenase-1 (HO-1), superoxide dismutase (SOD), and thioredoxin (Trx).

4. Materials and Methods

4.1. Database Information

This study used data from the FAERS database. FAERS is a large database of information on adverse drug events collected from all regions globally. In total, 14,836,467 cases reported between January 1, 2004, and March 31, 2022, were downloaded from the FDA website [32]. Figure 7 shows the procedure for creating the data tables.

The FAERS data used in this study comprised five tables: demographic characteristics of the patients and management information (DEMO), drug and biological information on reported adverse events (DRUG), and Medical Dictionary for Regulatory Activities (MedDRA) terminology coding adverse events (REAC). Duplicate data were removed from each table. The DRUG and REAC tables were incorporated into the DEMO table using the primary ID. This table was defined as the data table for analysis. All patient data contained in the DEMO table were used. Using this data table for analysis, the male and female subgroups were created.

4.2. Definitions of Adverse Events and Suspected Drugs

The International Council for Harmonization of Technical Requirements for Pharmaceuticals for Human Use developed MedDRA [33] as a repository of medical terminology for symptoms, signs, and diseases. In this study, hiccups were treated as a common term for adverse drug reactions. According to the adverse drug reaction reported, the drugs listed in FAERS are classified as primary, secondary, concomitant, and interacting drugs. In this study, all drugs were treated as suspected drugs.

4.3. Extraction of Drugs Suspected of Causing Hiccups

The drugs suspected to induce hiccups vary widely according to sex [11,12]. Thus, the overall analysis data table and subgroup tables divided based on sex were used to extract the suspected drugs (Figure 7). In each group, cross-tabulations were performed based on two categories (the presence/absence of hiccups and the presence/absence of the suspect drug) to calculate the reported odds ratio (ROR), and any significant result was evaluated using the Fisher’s exact test (Table 4). Supplementary File S1 shows the p values, ROR, and number of reports for each drug. A scatterplot (volcano plot) was then created by plotting the negative log of the p value from the Fisher’s exact test for each group on the y-axis and the natural log of the ROR (lnROR) on the x-axis. The drugs with a p value of < 0.001 based on the Fisher’s exact test and with am ROR of > 1 were suspected to cause hiccups. Supplementary File S1 shows the list of important suspected drugs and their ATC classification, which were divided according to the overall population and sex-specific subgroups.

R O R = \frac{n_{11} / n_{12}}{n_{21} / n_{22}} = \frac{n_{11} \cdot n_{22}}{n_{12} \cdot n_{21}}

4.4. MIE Activity Prediction Using the Toxicity Predictor

Toxicity Predictor [34] is a web-based application developed and operated by Meiji Pharmaceutical University. In this study, it was used to explore the AOPs of drug-induced hiccups. Further, Toxicity Predictor, which was developed as part of the Drug Discovery Information System Development Project of the Japan Agency for Medical Research and Development, can convert molecular structures obtained from input files into three-dimensional (3D) structures. In addition, it is a quantitative structure–activity relationship (QSAR) model-based toxicity prediction system that can evaluate agonist and antagonist activity against 56 MIEs. QSAR is a machine-learning algorithm that can learn mathematical associations between the chemical structure of a molecule and its biological or chemical activity to create a predictive model. MIE 3D structures and predicted results can be downloaded in the SDF and CSV formats, respectively [16].

To predict MIE activity using the Toxicity Predictor, the simplified molecular-input line-entry system was combined for each drug in each table. Using the Toxicity Predictor, the agonist and antagonist activity (MIE activity) of drugs targeting NRs and SPs for each group was calculated. The cutoff values for the predicted MIE activity were determined using the Youden method. The calculated MIE activities were normalized to establish a cutoff value of 0.5. Thus, the predicted labels of compounds with normalized predicted values of >0.5 were assigned with a value of 1. Meanwhile, those with predicted values of <0.5 were assigned with a value of 0. Supplementary File S2 depicts the predicted values.

4.5. MIE Associated with Hiccups

In each data table for the overall population and for the male and female subgroups, only the drugs suspected to be associated with hiccups were extracted. The selection criteria included the presence of at least 1,000 adverse event reports and a p value of ≤0.05, as determined using the Fisher’s exact test. The objective variable was a binary variable with an lnROR of > 0 or lnROR of <0. Univariate analysis was performed using the MIE activity (1-0) of each drug as the explanatory variable. In addition, to extract independent risk factors for each factor, logistic regression analysis was performed with lnROR as the objective variable and NRs and SPs, which were significant in the univariate analysis as explanatory variables.

4.6. Statistical Analysis

In this study, the objective variable was the binary classification of the ROR (lnROR ≥ 0, lnROR < 0). Hence, the ROR of a drug is significantly higher for a given adverse effect. For example, the adverse event is more likely to occur in other drugs. The explanatory variables were the probability for MIE activity calculated using Toxicity Predictor, normalized to a cutoff value of 0.5, and binary classified values. Univariate analysis was conducted to evaluate the NRs and SPs associated with hiccups. If the result of the Fisher’s exact test was significant in the positive direction, it was considered a factor associated with the induction of hiccups. By contrast, if the result of the Fisher’s exact test was significant in the negative direction, it was a factor related to the inhibition of hiccup induction. Using MIEs that were significant in the univariate analysis, a multivariate analysis was performed to extract independent risk factors for hiccup inversions. The likelihood ratio test and the Wald test were used to evaluate the independence of each factor. The significance level was set at p < 0.05. The pairwise method was used to evaluate internal correlation. An internal correlation was considered to exist if the Spearman's rank-order correlation coefficient [ρ2] was > 0.9. In this study, no internal correlation was found, and each factor was treated as independent. All analyses were performed using the JMP Pro 17 software (SAS Institute Inc., Cary, NC, USA).

4.7. Limitation

This study is entirely database based. Therefore, it has several limitations. The nuclear receptor activity values used in this study were predicted using QSAR and may differ from the activity values obtained experimentally. FAERS, the other database used in this study, which comprises a spontaneous report of adverse drug reactions, is particularly useful for detecting rare and serious adverse drug reactions. Moreover, it is an important source of information for evaluating the safety of drugs. In particular, FAERS has a stronger reporting bias because it includes cases in which patients who have taken the drug are enrolled as reporters [35]. In addition, the presence of drugs that have already been found to be associated with adverse events may affect the occurrence of adverse drug event-associated signals [36].

In this study, the number of reports and the p value obtained using the exact test were examined along with the ROR obtained via univariate analysis when detecting the signal, thereby avoiding a simple comparison of RORs and treating them semiquantitatively. This strategy was based on the hypothesis that a signal detection indicator considered highly significant based on the number of reports and p values would have excellent reliability [12].

5. Conclusion

The FAERS database was used to examine the association between NRs and SPs as well as drug-induced hiccups. We successfully extracted MIEs associated with drug-induced hiccups. Current case reports have reported about hiccups. However, to the best of our knowledge, this report first considered the pathogenesis of hiccups. This method can be applied to other adverse drug reactions, and it may be useful for elucidating the mechanism of adverse drug events.

Supplementary Materials

The following supporting information can be downloaded at: www.mdpi.com/xxx/s1, Table S1: The suspected drug list and ATC classification; Table S2: The MIE activity using QSAR Toxicity Predictor

Author Contributions

Conceptualization, Ryuichiro Hosoya; Data curation, Ryuichiro Hosoya; Funding acquisition, Ryuichiro Hosoya and Yoshihiro Uesawa; Investigation, Ryuichiro Hosoya; Methodology, Ryuichiro Hosoya, Tomoko Terajima and Yoshihiro Uesawa; Project administration, Reiko Ishii-Nozawa and Yoshihiro Uesawa; Resources, Ryuichiro Hosoya; Supervision, Yoshihiro Uesawa; Writing – original draft, Ryuichiro Hosoya; Writing – review & editing, Ryuichiro Hosoya, Reiko Ishii-Nozawa, Tomoko Terajima, Hajime Kagaya and Yoshihiro Uesawa.

Funding

This work was supported by a grant from the Japanese Society for Pharmaceutical Palliative Care and Sciences (grant number: 202102) and JSPS KAKENHI (grant number: 22K06707).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data were contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

Becker, D.E. Nausea, Vomiting, and Hiccups: A Review of Mechanisms and Treatment. Anesth. Prog. 2010, 57. [Google Scholar] [CrossRef] [PubMed]
Kolodzik, P.W.; Filers, M.A. Hiccups (Singultus): Review and Approach to Management. Ann. Emerg. Med. 1991, 20. [Google Scholar] [CrossRef]
Kondo, T.; Toyooka, H.; Arita, H. Hiccup Reflex Is Mediated by Pharyngeal Branch of Glossopharyngeal Nerve in Cats. Neurosci. Res. 2003, 47, 317–321. [Google Scholar] [CrossRef] [PubMed]
Nausheen, F.; Mohsin, H.; Lakhan, S.E. Neurotransmitters in Hiccups. Springerplus 2016, 5. [Google Scholar] [CrossRef]
Oshima, T.; Sakamoto, M.; Tatsuta, H.; Arita, H. GABAergic Inhibition of Hiccup-like Reflex Induced by Electrical Stimulation in Medulla of Cats. Neurosci. Res. 1998, 30, 287–293. [Google Scholar] [CrossRef]
Lee, G.W.; Kim, R.B.; Go, S. Il; Cho, H.S.; Lee, S.J.; Hui, D.; Bruera, E.; Kang, J.H. Gender Differences in Hiccup Patients: Analysis of Published Case Reports and Case-Control Studies. J. Pain Symptom Manage. 2016, 51, 278–283. [Google Scholar] [CrossRef]
Hosoya, R.; Tanaka, I.; Ishii-Nozawa, R.; Amino, T.; Kamata, T.; Hino, S.; Kagaya, H.; Uesawa, Y. Risk Factors for Cancer Chemotherapy-Induced Hiccups (CIH). Pharmacol. Pharm. 2018, 09, 331–343. [Google Scholar] [CrossRef]
Liaw, C.C.; Wang, C.H.; Chang, H.K.; Wang, H.M.; Huang, J.S.; Lin, Y.C.; Chen, J.S. Cisplatin-Related Hiccups: Male Predominance, Induction by Dexamethasone, and Protection against Nausea and Vomiting. J. Pain Symptom Manage. 2005, 30. [Google Scholar] [CrossRef] [PubMed]
Hayashi, M.; Sugimura, H.; Suga, Y.; Kawahara, M.; Aimiya, K.; Miyamoto, K. Study on Risk Factors for Hiccups Induced by Cisplatin-Based Chemotherapy. Iryo Yakugaku (Jpn. J. Pharm. Heal. Care Sci.) 2009, 35. [Google Scholar] [CrossRef]
Kang, J.H.; Hui, D.; Kim, M.J.; Kim, H.G.; Kang, M.H.; Lee, G.W.; Bruera, E. Corticosteroid Rotation to Alleviate Dexamethasone-Induced Hiccup: A Case Series at a Single Institution. J. Pain Symptom Manage. 2012, 43, 625–630. [Google Scholar] [CrossRef]
Hosoya, R.; Uesawa, Y.; Ishii-Nozawa, R.; Kagaya, H. Analysis of Factors Associated with Hiccups Based on the Japanese Adverse Drug Event Report Database. PLoS One 2017, 12. [Google Scholar] [CrossRef]
Hosoya, R.; Ishii-nozawa, R.; Kurosaki, K.; Uesawa, Y. Analysis of Factors Associated with Hiccups Using the FAERS Database. 2022. [Google Scholar]
Okunaka, M.; Kano, D.; Uesawa, Y. Nuclear Receptor and Stress Response Pathways Associated with Antineoplastic Agent-Induced Diarrhea. Int. J. Mol. Sci. 2022, 23. [Google Scholar] [CrossRef] [PubMed]
Olefsky, J.M. Nuclear Receptor Minireview Series. J. Biol. Chem. 2001, 276, 36863–36864. [Google Scholar] [CrossRef] [PubMed]
Diamanti-Kandarakis, E.; Bourguignon, J.P.; Giudice, L.C.; Hauser, R.; Prins, G.S.; Soto, A.M.; Zoeller, R.T.; Gore, A.C. Endocrine-Disrupting Chemicals: An Endocrine Society Scientific Statement. Endocr. Rev. 2009, 30, 293–342. [Google Scholar] [CrossRef] [PubMed]
Kurosaki, K.; Wu, R.; Uesawa, Y. A Toxicity Prediction Tool for Potential Agonist/Antagonist Activities in Molecular Initiating Events Based on Chemical Structures. Int. J. Mol. Sci. 2020, 21, 1–20. [Google Scholar] [CrossRef]
Oh, J.M.; Kim, I. Drug Repositioning Prediction for Psoriasis Using the Adverse Event Reporting Database. Front. Med. 2023, 10(1159453). [Google Scholar] [CrossRef]
Kinoshita, S.; Hosomi, K.; Yokoyama, S.; Takada, M. Inverse Association between Metformin and Amiodarone-Associated Extracardiac Adverse Events. Int. J. Med. Sci. 2020, 17, 302–309. [Google Scholar] [CrossRef] [PubMed]
Takada, M.; Fujimoto, M.; Motomura, H.; Hosomi, K. Inverse Association between Sodium Channel-Blocking Antiepileptic Drug Use and Cancer : Data Mining of Spontaneous Reporting and Claims Databases. 2016, 13. [Google Scholar] [CrossRef] [PubMed]
Inoue, K.; Yamazaki, H.; Manabe, Y.; Fukuda, C.; Hanai, K.; Fushiki, T. Transforming Growth Factor-Beta Activated during Exercise in Brain Depresses Spontaneous Motor Activity of Animals. Relevance to Central Fatigue. Brain Res. 1999, 846, 145–153. [Google Scholar] [CrossRef]
Dhandapani, K.M.; Wade, F.M.; Mahesh, V.B.; Brann, D.W. Astrocyte-Derived Transforming Growth Factor-β Mediates the Neuroprotective Effects of 17β-Estradiol: Involvement of Nonclassical Genomic Signaling Pathways. Endocrinology 2005, 146, 2749–2759. [Google Scholar] [CrossRef]
Wahl, S.M.; Allen, J.B.; McCartney-Francis, N.; Morganti-Kossmann, M.C.; Kossmann, T.; Ellingsworth, L.; Mai, U.E.H.; Mergenhagen, S.E.; Orenstein, J.M. Macrophage- and Astrocyte-Derived Transforming Growth Factor β as a Mediator of Central Nervous System Dysfunction in Acquired Immune Deficiency Syndrome. J. Exp. Med. 1991, 173, 981–991. [Google Scholar] [CrossRef]
KONDO, W.; OKUNO, M.; KOJIMA, S. Molecular Mechanism and Regulation of Latent TGF-ΒActivation. Jpn. J. Thromb. Hemost. 2003, 14, 210–219. [Google Scholar] [CrossRef]
Katsuno, M.; Adachi, H.; Banno, H.; Suzuki, K.; Tanaka, F.; Sobue, G. Transforming Growth Factor-β Signaling in Motor Neuron Diseases. Curr. Mol. Med. 2011, 11. [Google Scholar] [CrossRef]
Chang, F.Y.; Lu, C.L. Hiccup: Mystery, Nature and Treatment. J. Neurogastroenterol. Motil. 2012, 18, 123–130. [Google Scholar] [CrossRef] [PubMed]
Muñoz, M.D.; de la Fuente, N.; Sánchez-capelo, A. TGF-β/Smad3 Signalling Modulates GABA Neurotransmission: Implications in Parkinson’s Disease. Int. J. Mol. Sci. 2020, 21. [Google Scholar] [CrossRef] [PubMed]
Zhou, L.; Ouyang, R.; Luo, H.; Peng, Y.; Chen, P.; Ren, S.; Liu, G. Dysfunction of Nrf2-ARE Signaling Pathway: Potential Pathogenesis in the Development of Neurocognitive Impairment in Patients with Moderate to Severe Obstructive Sleep Apnea-Hypopnea Syndrome. Oxid. Med. Cell. Longev. 2018, 2018. [Google Scholar] [CrossRef]
Zhang, M.; An, C.; Gao, Y.; Leak, R.K.; Chen, J.; Zhang, F. Emerging Roles of Nrf2 and Phase II Antioxidant Enzymes in Neuroprotection. Prog. Neurobiol. 2013, 100, 30–47. [Google Scholar] [CrossRef]
Mishra, A.; Chaturvedi, P.; Datta, S.; Sinukumar, S.; Joshi, P.; Garg, A. Harmful Effects of Nicotine. Indian J. Med. Paediatr. Oncol. 2015, 36, 24–31. [Google Scholar] [CrossRef]
Fathi, M.; Saeedyan, S.; Kaoosi, M. Gamma-Amino Butyric Acid (GABA) Supplementation Alleviates Dexamethasone Treatment-Induced Oxidative Stress and Inflammation Response in Broiler Chickens. Stress 2023, 26. [Google Scholar] [CrossRef]
Sun, J.K.; Iwata, T.; Zigler, J.S.; Carper, D.A. Differential Gene Expression in Male and Female Rat Lenses Undergoing Cataract Induction by Transforming Growth Factor-β (TGF-β). Exp. Eye Res. 2000, 70, 169–181. [Google Scholar] [CrossRef]
FDA FAERS. Available online: https://www.fda.gov/drugs/questions-and-answers-fdas-adverse-event-reporting-system-faers/fda-adverse-event-reporting-system-faers-public-dashboard.
MedDRA. Available online: https://www.meddra.org/26.
Department of Medical Molecular Informatics, M.P.U. Toxicity Predictor. Available online: http://mmi-03.my-pharm.ac.jp/tox1/prediction_groups/new (accessed on 8 May 2023).
Maeda, R. JADER from Pharmacovigilance Point of View. Jpn J Pharmacoepidemiol 2014, 19, 51–56. [Google Scholar] [CrossRef]
Avillach, P.; Salvo, F.; Thiessard, F.; Miremont-Salamé, G.; Fourrier-Reglat, A.; Haramburu, F.; Bégaud, B.; Moore, N.; Pariente, A. Pilot Evaluation of an Automated Method to Decrease False-Positive Signals Induced by Co-Prescriptions in Spontaneous Reporting Databases. Pharmacoepidemiol. Drug Saf. 2014, 23. [Google Scholar] [CrossRef]

Figure 1. Drugs associated with hiccups. This volcano plot was created by plotting the negative logarithm of the p value (−log₁₀p) based on the Fisher’s exact test on the y-axis and the natural logarithm of the ROR (lnROR) on the x-axis. In all patients, drugs with > 1,000 reports and those considered as significant (p ≤ 0.001 based on the Fisher’s exact test) were included in the volcano plot. The color of the individual points represents differences in the log of the number of reports for each drug. The gray plots indicate drugs with a p value of > 0.001 or < 1,000 reported cases. The red line on the y-axis represents a p = 0.001. In this scatter plot, the signal is larger for the points (drugs) plotted in the upper right corner. The blue-to-red colors represent the number of times an adverse drug reaction was reported.

Figure 2. Drugs associated with hiccups in male participants This volcano plot was created by plotting the negative logarithm of the p value (−log₁₀p) based on the Fisher’s exact test on the y-axis and the natural logarithm of the ROR (lnROR) on the x-axis. The red line on the y-axis represents p = 0.001. In male patients, drugs with > 1,000 reports and those that are significant (p ≤ 0.001 based on the Fisher’s exact test) are included in the volcano plot. The color of the individual points represents differences in the log of the number of reports for each drug. The gray plots indicate drugs with p > 0.001 or <1,000 reported cases. In this scatter plot, the signal is larger for the points (drugs) plotted in the upper right corner. The line on the y-axis represents the total average. The blue-to-red colors represent the number of times an adverse drug reaction was reported.

Figure 3. Medicines associated with hiccups in female patients.

Figure 4. MIEs associated with hiccup induction in pharmacotherapy. This volcano plot was created by plotting the negative logarithm of the p value (−log10 p) using the Fisher’s exact test on the y-axis and the natural logarithm of the OR (ln OR) on the x-axis. The red line on the y-axis represents p = 0.05. In this scatter plot, signals associated with hiccup induction (positive signals) are plotted in the upper right, whereas those associated with inhibition (negative signals) are plotted in the upper left. (a) is the volcano plot for the male table, whereas (b) is the volcano plot for the female group. PR, progesterone receptor; TGFb, transforming growth factor-beta; Shh, sonic hedgehog signaling; ARlbd, androgen receptor lbd; ARE, antioxidant response element; CAR, constitutive androstane receptor; ARa, androgen receptor with antagonist.

Figure 5. TGF-β signaling and the GABA neurotransmission interaction.

Figure 6. Oxidative stress and Nrf2–ARE signaling.

Figure 7. Data analysis table.

Table 1. Univariate and multivariate analyses of MIEs associated with drug-induced hiccups.

MIEs	Activity type	Univariate analysis				Multivariate analysis
				95% CI			95% CI
		p value	Odds ratio	Lower	Upper	Odds ratio	Lower	Upper	p value (Likelihood ratio test)	p value (Wald test)
ERaant	agonist	0.010	0.21	0.06	0.69	0.28	0.06	1.26	0.088	0.097
PR	antagonist	0.011	0.43	0.23	0.82	0.61	0.30	1.26	0.181	0.181
ARaant	agonist	0.011	0.33	0.14	0.76	0.65	0.21	1.97	0.445	0.443
TGFb	agonist	0.014	3.37	1.30	8.75	4.59*	1.54	13.67	0.002	0.006
PR	agonist	0.014	0.20	0.05	0.76	0.66	0.12	3.55	0.625	0.628
Shh	agonist	0.026	0.40	0.18	0.89	0.63	0.25	1.60	0.330	0.329

MIEs: molecular initiating events, ERaant: estrogen receptor alfa with antagonist, PR: progesterone receptor, ARant: androgen receptor with antagonist, TGFb: transforming growth factor-β, Shh: sonic hedgehog. *; significant odds ratio on multivariate analysis.

Table 2. Univariate and multivariate analyses of MIEs associated with drug-induced hiccups in the male group.

MIEs	Activity type	Univariate analysis				Multivariate analysis
				95% CI			95% CI
		p value	Odds ratio	Lower	Upper	Odds ratio	Lower	Upper	p value (Likelihood ratio test)	p value (Wald test)
PR	agonist	0.005	0.326	0.152	0.698	0.60	0.22	1.62	0.318	0.318
CAR	antagonist	0.009	0.347	0.160	0.753	0.56	0.23	1.41	0.222	0.220
Sonic hedgehog	agonist	0.021	0.389	0.176	0.861	0.70	0.25	1.97	0.504	0.502
ARlbd	antagonist	0.025	0.401	0.189	0.851	0.62	0.22	1.70	0.351	0.348
TGFb	agonist	0.030	2.857	1.100	7.423	3.67*	1.29	10.45	0.011	0.015

PR: progesterone receptor, CAR: Constitutive androstane receptor, Shh: sonic hedgehog, ARlbd: androgen receptor lbd, TGFb: transforming growth factor-β. *; significant odd ratio on multivariate analysis.

Table 3. The activity type and result of univariate and multivariate of MIE.

MIEs	Activity type	Univariate analysis				Multivariate analysis
				95% CI		95% CI
		p value	Odds ratio	Lower	Upper	Odds ratio	Lower	Upper	p value (Likelihood ratio test)	p value (Wald test)
ARE	Agonist	0.01	0.08	0.01	0.70	0.12*	0.01	1.32	0.048	0.083
CAR	Antagonist	0.04	0.18	0.03	0.94	0.50	0.07	3.62	0.483	0.489
ARant	Agonist	0.05	0.20	0.05	0.88	0.24	0.04	1.20	0.073	0.080

ARE: antioxidant response element, CAR: constitutive androstane receptor, ARant: androgen receptor with antagonist, *: indicates the significant odds ratios on multivariate analysis.

Table 4. Cross-tabulation and formula for reportedodds ratios (RORs) of hiccups.

	Hiccups	Other adverse events
Suspected drugs	n₁₁	n₁₂
Other drugs	n₂₁	n₂₂

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.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer

Association between Molecular Initiating Events (MIEs) in Drug-Induced Hiccups

Abstract