Skip to main content

Analysis of drug-induced interstitial lung disease caused by herbal medicine using the Japanese Adverse Drug Event Report database

Abstract

Background

Drug-induced interstitial lung disease (DIILD) is a severe adverse event leading to morbidity and mortality. This study evaluated the adverse event indicators of DIILD and time-to-onset profiles following the daily intake of herbal drugs (Scutellariae radix [“ogon” in Japanese], Bupleuri radix [“saiko” in Japanese], and Pinelliae tuber [“hange” in Japanese]) using the Japanese Adverse Drug Event Report database. DIILD was defined in accordance with the Medical Dictionary for Regulatory Activities.

Methods

The Japanese Adverse Drug Event Report database contained 830,079 reports published between April 2004 and April 2023. The association between herbal medicines and DILLD was evaluated using the pharmacovigilance index as the reporting odds ratio (ROR), logistic regression models, propensity score-matching techniques, and Weibull shape parameters.

Results

The adjusted RORs using multivariate logistic regression models for Scutellariae radix (daily intake), Pinelliae tuber (daily intake), sex (male), age (≥ 60 years), Scutellariae radix (daily intake)*age (≥ 60 years), and Scutellariae radix (daily intake)* Pinelliae tuber (daily intake) were 1.47 (1.36 − 1.59), 1.05 (1.01 − 1.10), 1.45 (1.34 − 1.57), 1.92 (1.74 − 2.11), 3.35 (3.12 − 3.60), and 1.49 (1.46 − 1.53), respectively. DIILD onset profiles were evaluated using the Weibull shape parameter. A logistic plot of daily intake and onset of DIILD was drawn. ROR signals were detected in 32 of 54 herbal medicines, including Scutellariae radix, Bupleuri radix, and Pinelliae tuber. The median duration (days) (interquartile range) to DIILD onset was 36.0 (27.0–63.0) for Saikokaryukotsuboreito, 35.0 (21.0–55.0) for Saireito, and 31.0 (13.5–67.5) for Shosaikoto. The Weibull shape parameter beta (95% confidence interval) values for Saikokaryukotsuboreito, Saireito, and Shosaikoto were 1.36 (1.08–1.67), 1.36 (1.20–1.52), and 1.31 (0.98–1.68), respectively.

Conclusions

DIILD demonstrated a dose-dependent to crude drugs. Clinicians should strive for the early detection of DIILD and avoid the inadvertent administration of herbal medicines.

Peer Review reports

Background

Drug-induced interstitial lung disease (DIILD) is the most common drug-induced lung toxicity that occurs when drug exposure causes inflammation and, eventually, fibrosis of the lung interstitium [1, 2]. DIILD is a severe adverse event (AE) that leads to increased morbidity and mortality rates. In recent years, DIILD has been reported as a risk factor for coronavirus disease 2019 [3]. Therefore, all healthcare professionals should be familiar with potential DIILD as early detection is important to prevent its progression.

Drugs known to cause DIILD include amiodarone, interferons, anticancer drugs (e.g., gefitinib), and Japanese herbal medicines “Kampo” (e.g., Shosaikoto and Saireito) [4]. Kampo is formulated from natural agents, and in Japan, 148 “Kampo extract formulations for prescription” have been approved for ethical use [5, 6]. The percentage of events related to lung injury was 27.8% [7]. Tsukiyama et al. reported a clinical case of DIILD caused by Shosaikoto in 1989 [8], and numerous other clinical cases have been reported. Saikozai is a general term for Japanese herbal medicines mainly containing Scutellariae radix (SR, “ogon” in Japanese) and Bupleuri radix (BR, “saiko” in Japanese) (e.g., Shosaikoto and Saikokeishikankyoto) and has been reported to cause DIILD [9,10,11,12,13].

Although Kampo is a mixture of multiple crude drugs, other crude drugs, such as SR, BR, and Pinelliae tuber (PT, “hange” in Japanese), show positive results in the lymphocyte stimulation test and are thought to be associated with DIILD [8, 11,12,13]. Both SR and BR are thought to cause lung injury [14]. According to a review by Enomoto et al., the most common formulas are Shosaikoto (26%), Saireito (16%), Seishinrenshiin (8%), and Bofutsushosan (8%) [15]. SR occurred in 86% of the formulas, BR in 59%, and PT in 58% [15]. Furthermore, the Ministry of Health, Labour, and Welfare of the Government of Japan has issued “Pharmaceuticals and Medical Devices Safety Information” on DIILD caused by herbal medicines [16, 17]. However, few studies have examined the relationship between the risk of DIILD and each herbal ingredient or its contents. Therefore, we focused on the constituents of herbal medicines, such as SR, BR, and PT.

Although the risk of DIILD has been described in several studies, this information does not reflect the complexities in real-world clinical practice. Spontaneous reporting systems (SRSs), such as the Japanese Adverse Drug Event Report (JADER) database of the Pharmaceuticals and Medical Devices Agency (PMDA), have been used for pharmacovigilance assessments. The JADER database is a valuable tool for postmarketing surveillance that reflects real-world settings.

In this study, we evaluated the relationship between herbal medicines and DIILD using data from the JADER database. We evaluated AE signals using a pharmacovigilance index known as the reporting odds ratio (ROR), which is currently employed by the Japanese PMDA. We assessed the possible relationships between age, sex, and daily SR, BR, and PT intakes using adjusted RORs and propensity score (PS) matching techniques. Furthermore, we assessed the time-to-onset profiles of DIILD associated with herbal medicines.

Methods

Figure 1 summarizes the workflow of downloading JADER data from the PMDA website, excluding reports that were unsuitable for analysis due to missing values, and finally providing the data for analysis (Fig. 1). The JADER database provides information on cases reported by pharmaceutical corporations or medical institutions since 2004, based on the Act on Securing Quality, Efficacy, and Safety of Products Including Pharmaceuticals and Medical Devices (Act No. 145 of 1960), and cases reported by medical institutions since 2013, based on the Immunization Law. Healthcare professionals, marketing approval holders, patients, and consumers voluntarily report adverse events to the PMDA. All data from the JADER database are cleaned and fully anonymized by the regulatory authority of Japan. These data are provided in CSV format, are made publicly available, and can be downloaded from the PMDA website (https://www.pmda.go.jp) [18]. When additional information or corrections are reported in the CSV file, the PMDA updates the appropriate information accordingly. The data in the CSV file are not evaluated by the PMDA with regards to the association between individual drugs and adverse events. The database comprised four tables of data: patient demographic information, drug information, AEs, and primary illness. We integrated a relational database based on the four tables using FileMaker Pro 12 software (FileMaker, Inc., Santa Clara, CA, United States) according to the ASCII Entity Relationship Diagram, which is publicly available on the PMDA website [18]. Drug data included the following role codes assigned to each drug according to its association with AEs: “suspected,” “concomitant,” and “interacting drugs.” We extracted and analyzed the “suspected drug” records.

Fig. 1
figure 1

Flowchart of data analysis

We listed 147 Kampo extract products approved by the Ministry of Health, Labor, and Welfare [19] and selected 54 herbal medicines containing SR, BR, or PT (Fig. 1). Regardless of whether Kampo included SR, BR, or PT, we classified Kampo into 2^3 = 8 groups. We analyzed the signals for seven groups, excluding those that did not include the three types of crude drugs. When a report contained several Kampo, we counted them as long as the Kampo were in the same group (e.g., a combination of Bofutsushosan [SR] and Orengedokuto [SR] was counted as one case; a combination of Bofutsushosan [SR] and Yokukansan [BR] was not counted as one case; and a combination of Bofutsushosan [SR] and Otsujito [SR_BR] was not counted as one case). When analyzing signals for sex and age, we eliminated incomplete reports (e.g., missing, unknown, and aged) on sex and age from the subsets (Fig. 1).

To extract cases of DIILD, we used “interstitial lung disease” as the preferred term (code:10,022,611) following the terminology preferred by the Medical Dictionary for Regulatory Activities (MedDRA) (version 23.0) [20]. According to the Introductory Guide MedDRA Version 23.0, a preferred term is defined as a clear descriptor (a single medical concept) regarding symptoms, signs, disease, diagnosis, therapeutic indication, and surgical or medical procedure [21]. The contributors only reported AEs according to the International Council for Harmonisation of Technical Requirements for Pharmaceuticals for Human Use (ICH) E2B and international safety reporting guidelines, and relied on the definitions provided by MedDRA.

ROR is the ratio of the odds of reporting DIILD versus all other AEs for a given herbal medicine to the reporting odds for all other herbal medicines or drugs present in the database [22, 23]. The ROR was calculated using a two-by-two contingency table, where “a” is the number of AE reports in which patients received the herbal medicine and manifested DIILD; “b” is the number of AE reports in which patients received herbal medicine but did not present DIILD; “c” is the number of AE reports in which patients did not receive herbal medicine of interest and had DIILD; and “d” is the number of AE reports in which patients did not receive herbal medicine of interest and did not present DIILD (Fig. 2). In the case of the worked example, it is the ratio of the odds of DIILD divided by the odds of all other AEs (ROR of “Bofutsushosan,” (103/36745)/(307/793334) = 7.3, Table 1). RORs were expressed as point estimates with 95% confidence intervals (CIs). For signal detection, general qualitative judgments are viable; detection of a signal depends on whether the signal indices exceed predefined thresholds, and crude ROR values < 1 indicate no potential exposure-event associations. Crude ROR is an indicator of AEs when the lower limit of the 95% CI is > 1 and the reported number of AEs is > 2 [22, 23].

Fig. 2
figure 2

Two-by-two used for the calculation of the reporting odds ratio

Table 1 Number of reports and ROR of drug-induced interstitial lung disease for each herbal medicine for analysis of whole data

We used only reports with complete or calculable data on sex, age, and daily intake of each crude drug. The effectiveness of the explanatory variables was evaluated using a stepwise method with a significance level of 0.05 (forward and backward) [24]. We input the reporting year parameter into the following formula and evaluated the influence of the explanatory variable using the likelihood ratio test:

$$\begin{array}{l}\text{log}\left(odds\right)={\mathrm\beta}_{10}+{\mathrm\beta}_{11}Y+{\mathrm\beta}_{12}D1+{\mathrm\beta}_{13}D2+{\mathrm\beta}_{14}S+{\mathrm\beta}_{15}A+{\mathrm\beta}_{16}D1\ast A+{\mathrm\beta}_{17}D1\ast D2\\(\mathrm D1\;=\;\mathrm{daily}\;\mathrm{intake}\;\lbrack\mathrm{SR}\rbrack,\;\mathrm D2\;=\;\mathrm{daily}\;\mathrm{intake}\;\lbrack\mathrm{PT}\rbrack,\;\mathrm S\;=\;\mathrm{sex}\;\lbrack\mathrm{male}\;=\;1,\;\mathrm{female}\;=\;0\rbrack,\;\mathrm A\;=\;\mathrm{age}\;\lbrack\geq60\;\mathrm{years}\;=\;1,\;<60\;\mathrm{years}\;=0\rbrack)\end{array}$$
(1)

PS matching was used as an assessment approach to reduce selection bias. PS matching is a statistical matching technique used to construct matched sets with similar covariate distributions without requiring close or exact matches for all individual variables [25,26,27]. Sex, age, and daily crude drug intake were included in multivariate logistic regression analysis. The presence or absence of DIILD was also assessed. Nearest-neighbor matching was performed based on the calculated PS between the presence and absence of crude drugs. The caliper width of 0.2 of the standard deviation of the PS logit was used. The standard mean difference (SMD) was used as a covariate balance indicator between the presence and absence of crude drugs. An SMD value less than 0.1 was regarded as balanced.

Using a logistic regression model, we estimated the relationship between the daily intake of crude drugs (SR, BR, and PT) and DIILD onset. All of the products used in this analysis are granule “Kampo extract formulations for prescription.” The amount of Kampo extract formulation prescribed per package was specified for each crude drug. For example, two packages (3.75 g) of Shosaikoto granules contained 2.25 g of Shosaikoto extract. Bupleurum Root, Pinellia Tuber, Scutellaria Root, Jujube, Ginseng, Glycyrrhiza, and Ginger were present in the Shosaikoto extract at proportions of 3.5, 2.5, 1.5, 1.5, 1.0, and 0.5, respectively. The daily intake of each crude drug was calculated based on the amount of herbal medicine intake in the JADER reports and information on the package insert of each herbal medicine. The following formula is used:

$$\begin{array}{c}\text{log}\left(\mathit{odds}\right)=\beta_{20}+\beta_{21}D\\(\text{D}\hspace{0.17em}=\hspace{0.17em}\mathrm{dailyintake}\lbrack\text{SR},\mathrm{BR},\mathrm{PT}\rbrack)\end{array}$$
(2)

The receiver operating characteristic (ROC) curve determined the accuracy of model predictions for treatment allocation. For the ROC curve, the horizontal axis shows the false-positive fraction (1-specificity) and the vertical axis shows the true-positive fraction (sensitivity); the cut-off point has the highest sensitivity while maintaining high specificity.

The period from the prescription of specific drugs to the occurrence of a specific AE was evaluated using median time, quartile, and Weibull shape parameters (WSP) [28]. Duplicate prescriptions for the same patient and reports without a complete date of drug administration or AE onset were excluded. Truncation should be considered when estimating the time to onset of AEs using SRS data. We chose an analysis period of 365 days after the start of drug administration to focus on the onset of AEs within one year after drug administration. The occurrence rate of AEs after prescription depends on the causal mechanism, and often varies over time. In contrast, AEs not associated with the drug occurred at a constant rate.

Median, quartile, and WSP tests were used to evaluate the time-to-onset data. Reports that did not include complete AE occurrences or prescription start times were excluded. Time-to-onset data were calculated using the start date of drug administration and the date that the AE occurred. It is necessary to consider truncation when estimating the time-to-onset of adverse events from JADER data. We selected an analysis period of 365 days after the start of drug administration to focus on the onset of AEs within a year after drug administration. We applied a Weibull distribution to the survival analysis. In this study, we did not use a proportional hazard model or a parametric modeling method such as an accelerated failure time model. The WSP test was used for the statistical analysis of the time-to-onset data and can describe the non-constant ratio of the incidence of AEs. The scale parameter α and shape parameter β determine the scale and shape of the distribution function, respectively. The shape parameter β of the Weibull distribution indicates the hazard without a reference population. If β = 1, then the hazard is estimated to be constant over time (random failure type). If β > 1 and the 95% CI of β excludes the value of 1, the hazard is estimated to increase over time (wear-out failure type). Finally, if β < 1 and the 95% CI β excludes value 1, the hazard is estimated to decrease over time (initial failure type) [29, 30].

Data analysis was performed using the JMP version 11.0 software (SAS Institute, Inc., Cary, NC, USA).

Results

The JADER database contains 830,079 reports published between April 2004 and April 2023. The ROR signals, including SR, BR, and PT, were detected in 32 of the 54 herbal medicines. The ROR (95% CI) for SR, BR, PT, and the combination of SR, BR, and PT (SR_BR_PT) were 9.0 (7.9 − 10.3), 2.6 (2.1 − 3.3), 4.1 (3.2 − 5.2), and 13.6 (11.9 − 15.4), respectively (Table 1). After excluding missing reports, the number of reports for sex and age were 6,845 and 6,741, respectively. The reported number of male, female, < 20 years, the 20 s, 30 s, 40 s, 50 s, 60 s, 70 s, and ≥ 80 s were 808, 719, 6, 7, 16, 50, 153, 385, 534, and 354, respectively (Table 2).

Table 2 Number of reports and ROR of drug-induced interstitial lung disease for sex and age

The adjusted RORs for SR (daily intake), PT (daily intake), sex (male), age (≥ 60 years), SR (daily intake)*age (≥ 60 years), and SR (daily intake)*PT (daily intake) were 1.47 (1.36 − 1.59), 1.05 (1.01 − 1.10), 1.45 (1.34 − 1.57), 1.92 (1.74 − 2.11), 3.35 (3.12 − 3.60), and 1.49 (1.46 − 1.53), respectively (Table 3).

Table 3 The statistics of multivariate logistic regression of daily intake and onset of drug-induced interstitial lung disease

The ROC curve of PS determined the accuracy of the model predictions for treatment allocation. The areas under the ROC curves for the SR, BR, and PT were 0.77708, 0.72868, and 0.79195, respectively (data not shown). In the presence or absence of the SR, the SMD decreased for all factors (Table 4). In contrast, the SMD for age increased with the presence or absence of BR, and the SMD for sex increased with the presence or absence of PT (Table 4). Before PS matching, the crude ROR (95% CI) for SR, BR, and PT were 2.79 (2.42–3.21), 1.55 (1.34–1.79), and 1.97 (1.69–2.29), respectively. After PS matching, the RORs for the SR, BR, and PT were 2.84 (2.25–3.58), 0.89 (0.73–1.09), and 1.06 (0.84–1.35), respectively.

Table 4 Comparison of the number of reports for each factor before and after propensity score matching

We analyzed the relationship between the daily intake of crude drugs and the onset of DIILD (Table 5 and Fig. 3). The estimated β21 values (95% CI) for SR, BR, and PT were 0.40 (0.35–0.46) (p < 0.0001), 0.14 (0.11–0.17) (p < 0.0001), and 0.14 (0.10–0.17) (p < 0.0001), respectively (Table 5). The slope of the graph was the largest for SR (Table 521 = 0.40) and Fig. 3). According to the ROC curves for SR (AUC = 0.63572), BR (AUC = 0.56849), and PT (AUC = 0.57047), the cut-off values (g/day) were 0.96, 3.30, and 1.60, respectively (Fig. 4).

Table 5 The statistics of logistic regression of daily intake and onset of drug-induced interstitial lung disease
Fig. 3
figure 3

Logistic plots of daily intake (gram per day) and the onset of drug-induced interstitial lung disease (yes = 1, no = 0)

Fig. 4
figure 4

Receiver operating characteristic curve of SR, BR, and PT

To analyze the time-to-onset of DIILD, we evaluated 21 herbal medicines used in > 10 reported cases. The median duration (days) (interquartile range) to the onset of DIILD was 36.0 (27.0–63.0) days for Saikokaryukotsuboreito (containing SR_BR_PT), 35.0 (21.0–55.0) days for Saireito (containing SR_BR_PT), and 31.0 (13.5–67.5) days for Shosaikoto (containing SR_BR_PT) (Figs. 5 and 6). The Weibull shape parameters β (95% CI) of Saikokaryukotsuboreito (containing SR_BR_PT), Saireito (containing SR_BR_PT), and Shosaikoto (containing SR_BR_PT) were 1.36 (1.08–1.67), 1.36 (1.20–1.52), and 1.31 (0.98–1.68), respectively.

Fig. 5
figure 5

Time-to-onset profiles of drug-induced interstitial lung disease

Fig. 6
figure 6

Kaplan–Meier curve plot for each group (SR, BR, PT, SR_BR, SR_PT, and SR_BR_PT)

The underlying diseases of patients with DIILD are summarized in Table 6); the number of reported cases of hypertension and constipation in the underlying disease of SR was 160 and 81, respectively. The number of reported cases of hypertension and dementia in the underlying diseases of BR were 161 and 107, respectively. The number of reported cases of hypertension and gastroesophageal reflux disease among the underlying diseases of PT were 87 and 28, respectively. The number of reported cases of hemorrhoids and hypertension in the underlying diseases of SR_BR were 60 and 33, respectively. The number of reported cases of hypertension and periarthritis in the underlying diseases of SR_PT were 33 and 21, respectively. The number of reported cases of diabetes and depression in the underlying diseases of SR_BR_PT were 40 and 33, respectively.

Table 6 The underlying diseases of the patients with DIILD

Discussion

AE signals for DIILD were detected in 59% of the herbal medicines containing SR, BR, and PT in the JADER database. We could not evaluate the effect of AEs on the incidence rates according to sex and age based on the crude ROR, as the crude ROR indicated an increased risk of AE reporting, but did not present the risk of AE occurrence in absolute terms. Multivariate logistic regression offers the advantage of controlling for covariates and can be used to analyze the use of interaction terms in more detail [31,32,33]. Our adjusted ROR results suggest an interaction between daily SR intake and age. Male sex and age ≥ 60 years were risk factors for DIILD [2]. Mortality rates in interstitial lung disease are highest in men and older adults [34]. In Japan, approximately 60% of patients using Japanese herbal medicines are over the age of 60 years [35,36,37]. This information is important for healthcare professionals.

Both the ROR (1.50 [1.42–1.58]; Table 5) and ROR (2.84 [2.25–3.58]) of SR after PS matching were higher than those of BR and PT. We observed that increased daily intake of SR, BR, and PT was associated with an increased incidence of DIILD and determined the dose response to crude drugs for DIILD onset. The inadvertent administration of herbal medicines is a risk factor for DIILD. The slope (β21) of SR was larger than that of the other drugs. This result may reflect the differences in the pharmacological actions of SR and other herbal medicines. However, further in vivo studies are required to confirm these results. Based on ROC curve analysis, if the daily intake of SR was greater than 0.96 g/day, this trend was more likely to appear than the following administration at lower doses.

Two mechanisms are involved in the pathogenesis of DIILD: direct dose-dependent toxicity and immune-mediated toxicity [38]. Both SR and BR are thought to cause lung injury [14]. SR mainly contains baicalein and baicalin as chemical ingredients, which may be associated with the onset of DIILD. A previous study suggested that baicalein induces apoptosis in human lung fibroblasts [39], which is essential for the growth and differentiation of type 2 alveolar epithelial cells [40], which are necessary for normal lung repair [41]. Pulmonary fibrosis is an interstitial lung disease caused by fibrosis due to remodeling associated with defective re-epithelization [42, 43]. Baicalein reduced the number of lung fibroblasts, negatively affecting the growth and differentiation of type 2 alveolar epithelial cells. Under these conditions, the lungs are not repaired properly, leading to fibrillation. Based on these findings, SR was considered a crude drug associated with DIILD. A high dose of baicalin (400–1600 mg/kg) induced fibrosis in Sprague–Dawley rats by activating the TGF-β/Smad signaling pathway [44]. According to the package inserts of Japanese herbal medicines, the daily dose of baicalein is estimated to be less than 200 mg/d [45]. The results from the animal studies cannot be extrapolated to patients; however, our cut-off value of 0.96 g/day may be worthy of further investigation.

Several studies of BR and PT have demonstrated dose-dependent toxicity in vitro and in organs such as the liver [46,47,48,49,50,51]. For PT, we observed a few reports on lymphocyte stimulation tests [8, 11,12,13]. Previous studies on BR and PT have reported pulmonary toxicity [11,12,13]. Our results did not demonstrate clear correlations between daily intake, BR, PT, and DIILD onset. The crude components of herbal medicines should be considered when determining their appropriateness for daily consumption.

More than 50% of DIILD cases associated with herbal medicine in all groups, except BR were observed within 50 days in the real-world data used in the time-to-onset analysis. Herbal medicines containing BR showed a longer time to DIILD onset. DIILD can occur after 180 days (six months) of administration of Orengedokuto (SR), Hochuekkito (BR), Yokukansan (BR), Bakumondoto (PT), Otsujito (SR_BR), Hangeshashinto (SR_PT), Saibokuto (SR_BR_PT), and Saireito (SR_BR_PT). This information may enable clinicians to perform early interventions to prevent DIILD progression. Herbal medicines (56% [9/16]) containing SR exhibited wear-off failure effects, and the rate of DIILD increased over time. SR may be cytotoxic because it is often expressed immediately if the expression mechanism is allergic, and delayed if it is cytotoxic [1, 38].

Our study has several limitations. SRS, such as JADER, are subject to overreporting, underreporting, missing data, exclusion of data from healthy individuals, lack of a denominator, and the presence of confounding factors. There was a lack of comparisons between groups. Moreover, it should be emphasized in SRS studies that disproportionality measures (crude ROR) do not allow for risk quantification. RORs offer a rough indication of the signal strength. In absolute terms, the ROR indicates an increased risk of AE reporting, and not the risk of AE occurrence. Comparing the intensity of risk based on measures, such as the crude ROR value, in a disproportionality analysis is not recommended [52].

The ROR is a pharmacovigilance index used by Japanese and Netherlands administrative authorities. In SRS studies, it is common practice to summarize the number of reports and RORs to provide an overview of the data (Table 1). Multiple logistic regression analysis and PS matching were used to avoid inappropriate data interpretation of the crude ROR for conventional signal detection. In this study, the adjusted ROR by multiple logistic regression analysis was described, and the ROR of the three herbal medicines were compared based on the adjusted ROR using the dataset after PS matching (Table 4). Several previous studies have compared RORs after matching backgrounds using PS matching, and we believe that our analysis method is valid. However, the effects of multivariate logistic regression and PS matching techniques are partial and subject to bias due to unmeasured confounders. Therefore, careful attention must be paid when interpreting results from the JADER database. Further epidemiological studies are required to determine the effects of these factors.

Depending on the analysis period of the time-to-onset profile, the hazard trends (increase, decrease, or constant) may differ [28]. Therefore, the follow-up duration must be carefully determined. To our knowledge, there is no gold standard for analyzing SRS research. The parameters were determined according to the dataset and the purpose of the study. In our previous research on DIILD of 110 drugs, including herbal medicines, we chose an analysis period of 730 days to focus on the onset of AEs within 2 years [4]. In the research, the median of “Saireito” and “Shosaikoto” were 35.0 and 33.0 days, respectively. Furthermore, the median time to onset of DIILD by anticancer drugs known to be cytotoxic is 262 days at most [53,54,55]. In this study, we set the analysis period to one year because the longer the observation period, the more unknown factors other than drugs may affect AEs. Furthermore, one limitation of the WSP tool is that several AE mechanisms are not time dependent.

According to the National Health and Nutrition Survey by the Ministry of Health, Labor, and Welfare, smoking is a risk factor for DIILD, and men exhibit a higher rate of smoking than women in Japan [56]. Smoking history was not evaluated in this study because of a lack of data in the JADER database. Patients with DIILD have various underlying diseases such as hypertension; therefore, the evaluation of DIILD risk due to these differences in patient backgrounds should be considered in the future.

Herbal medicines contain several crude drugs, and DIILD is caused by herbal medicines that do not contain SR, BR, or PT [57]. A difference in the amount of baicalin in herbal formulations among pharmaceutical companies, even for the same formula, has been reported, which was not considered in our study [45]. In this study, we selected three herbal medicines that have been widely recognized as risky, based on previous studies.

Conclusions

We determined the possible risk and time-of-onset of DIILD caused by herbal medicines. We observed a correlation between crude drug intake and DIILD onset. Clinicians should exercise caution when prescribing herbal medicines and should carefully consider the timing and dose of administration. Additionally, inadvertent administration of herbal medicines should be avoided.

Availability of data and materials

Data used in this study are available from the following. links: https://www.pmda.go.jp/safety/info-services/drugs/adr-info/suspected-adr/0004.html.

Abbreviations

DIILD:

Drug-induced interstitial lung disease

SR:

Scutellariae Radix

BR:

Bupleuri radix

PT:

Pinelliae tuber

SRS:

Spontaneous reporting system

ROR:

Reporting odds ratio

AE:

Adverse event

JADER:

Japanese Adverse Drug Event Report

PMDA:

Pharmaceutical and Medical Devices Agency

PS:

Propensity score

MedDRA:

Medical Dictionary for Regulatory Activities

WSP:

Weibull shape parameter

SMD:

Standard mean difference

ROC:

Receiver operating characteristic

References

  1. Schwaiblmair M, Behr W, Haeckel T, Märkl B, Foerg W, Berghaus T. Drug induced interstitial lung disease. Open Respir Med J. 2012;6:63–74. https://doi.org/10.2174/1874306401206010063. PMID: 22896776.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Ministry of Health, Labour and Welfare. The manual for handling disorders due to adverse drug reactions, interstitial lung disease. 2019. https://www.mhlw.go.jp/topics/2006/11/tp1122-1b.html. Accessed 15 Jun 2022.

  3. Podolanczuk AJ, Wong AW, Saito S, Lasky JA, Ryerson CJ, Eickelberg O. Update in Interstitial lung disease 2020. Am J Respir Crit Care Med. 2021;203:1343–52. https://doi.org/10.1164/rccm.202103-0559UP. PMID: 33835899.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Matsumoto K, Nakao S, Hasegawa S, Matsui T, Shimada K, Mukai R, et al. Analysis of drug-induced interstitial lung disease using the Japanese Adverse Drug Event Report database. SAGE Open Med. 2020;8:2050312120918264. https://doi.org/10.1177/2050312120918264. PMID: 32528682.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Tsutani K. Kusuri ha risuku : Kampoyaku kara seiyoyaku wo miru. Jpn J Pharmacoepidemiol. 2010;15:31 Japanese.

    Google Scholar 

  6. Kuchta K, Cameron S. Editorial: Kampo medicine in a modern context: ethnopharmacological perspectives. Front Pharmacol. 2022;13:971254. https://doi.org/10.3389/fphar.2022.971254. PMID: 36120328.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Shimada Y, Fujimoto M, Nogami T, Watari H. Adverse events associated with ethical kampo formulations: analysis of the Domestic Adverse-Event Data Reports of the Ministry of Health, Labor, and Welfare in Japan. Evid Based Complement Alternat Med. 2019;2019:1643804. https://doi.org/10.1155/2019/1643804. PMID: 31118950.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Tsukiyama K, Tasaka Y, Nakajima M, Hino J, Nakahama C, Okimoto N, et al. A case of pneumonitis due to Sho-saiko-to. Nihon Kyobu Shikkan Gakkai Zasshi. 1989;27:1556–61 PMID: 2630777.

    CAS  PubMed  Google Scholar 

  9. Itoh T, Fujimoto H, Umekawa K, Rensha K, Minami K, Shoji S. A case of Sai-rei-to-induced pneumonitis. Nihon Kokyuki Gakkai Zasshi. 2006;44:833–7 PMID: 17144582.

    PubMed  Google Scholar 

  10. Miyagawa T, Mochizuki Y, Nakahara Y, Kawamura T, Sasaki S, Tsukamoto H, et al. A case of drug-induced pneumonitis due to Sai-rei-to. Nihon Kokyuki Gakkai Zasshi. 2009;47:47–51 PMID: 19198236.

    PubMed  Google Scholar 

  11. Katou K, Mori K. Autoimmune hepatitis with drug-induced pneumonia due to Sho-saiko-to. Nihon Kokyuki Gakkai Zasshi. 1999;37:641–6 PMID: 10496105.

    CAS  PubMed  Google Scholar 

  12. Okada Y, Watanabe K, Suzuki Y, Suzuki K, Ito G, Muranushi A, et al. A case of hepatitis and intersitial pneumonitis induce by Hange-shasin-to and Sho-saiko-to. Kampo Med. 1999;50:57–65.

    Article  Google Scholar 

  13. Yamamoto T, Tsutsui N, Kazawa T, Sasagawa M, Sato Y. A case of drug-induced pneumonitis due to Ougon and Saiko involved in Sai-rei-to showing strong uptake of Ga scintigraphy. Niigata Igakkai Zasshi. 2016;130:699–705.

    Google Scholar 

  14. Komiya K, Ishii H, Ohama M, Uchida M, Tsubone T, Iwashita T, et al. Sai-rei-to-induced lung injury: a case report and brief review of the literature. Intern Med. 2012;51:3421–5. https://doi.org/10.2169/internalmedicine.51.8705. PMID: 23257532.

    Article  PubMed  Google Scholar 

  15. Enomoto Y, Nakamura Y, Enomoto N, Fujisawa T, Inui N, Suda T. Japanese herbal medicine-induced pneumonitis: a review of 73 patients. Respir Investig. 2017;55:138–44. https://doi.org/10.1016/j.resinv.2016.11.007. PMID: 28274529.

    Article  PubMed  Google Scholar 

  16. Ministry of Health, Labour and Welfare. Pharmaceutical and Food Safety Bureau. Pharmaceuticals and Medical Devices Safety Information No. 146. 1998. https://www.pmda.go.jp/safety/info-services/drugs/calling-attention/safety-info/0147.html. Accessed 18 May 2022.

  17. Ministry of Health, Labour and Welfare. Pharmaceutical and Food Safety Bureau. Pharmaceuticals and Medical Devices Safety Information No. 158. 2000. https://www.pmda.go.jp/safety/info-services/drugs/calling-attention/safety-info/0089.html. Accessed 18 May 2022.

  18. PMDA website. https://www.pmda.go.jp/safety/info-services/drugs/adr-info/suspected-adr/0004.html. 2004. Accessed 18 May 2022.

  19. Pharmaceuticals and Medical Devices Agency. Find Review reports, PI; 2022. https://www.pmda.go.jp/english/search_index.html. Accessed 15 Jun 2022.

  20. Japanese Maintenance Organization. Medical Dictionary for Regulatory Activities MedDRA, version 23.0; 2021. https://www.jmo.pmrj.jp. Accessed 15 Jun 2022.

  21. The international council for harmonisation of technical requirements for pharmaceuticals for human use (ICH), Introductory guide MedDRA Version 23.0. https://admin.meddra.org/sites/default/files/guidance/file/intguide_%2024_0_English.pdf Accessed 28 Feb. 2024.

  22. van Puijenbroek EP, Bate A, Leufkens HG, Lindquist M, Orre R, Egberts AC. A comparison of measures of disproportionality for signal detection in spontaneous reporting systems for adverse drug reactions. Pharmacoepidemiol Drug Saf. 2002;11:3–10. https://doi.org/10.1002/pds.668. PMID: 11998548.

    Article  CAS  PubMed  Google Scholar 

  23. Poluzzi E, Raschi E, Piccinni C, De Ponti F, Karahoca A. Data Mining Techniques in Pharmacovigilance: analysis of the publicly accessible FDA Adverse Event Reporting System (AERS). InTech. 2012. p.265–302. https://doi.org/10.5772/50095.

  24. Nakao S, Hasegawa S, Shimada K, Mukai R, Tanaka M, Matsumoto K, et al. Evaluation of anti-infective-related clostridium difficile-associated colitis using the Japanese Adverse Drug Event Report database. Int J Med Sci. 2020;17:921–30. https://doi.org/10.7150/ijms.43789. PMID: 32308545.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Rosenbaum PR, Rubin DB. The central role of the propensity score in observational studies for causal effects. Biometrika. 1983;70:41–55.

    Article  MathSciNet  Google Scholar 

  26. Akimoto H, Oshima S, Negishi A, Ohara K, Ohshima S, Inoue N, et al. Assessment of the risk of suicide-related events induced by concomitant use of antidepressants in cases of smoking cessation treatment with varenicline and assessment of latent risk by the use of varenicline. PLoS One. 2016;11:e0163583. https://doi.org/10.1371/journal.pone.0163583. PMID: 27657721.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Nakao S, Hasegawa S, Umetsu R, Shimada K, Mukai R, Tanaka M, et al. Pharmacovigilance study of anti-infective-related acute kidney injury using the Japanese Adverse Drug Event Report database. BMC Pharmacol Toxicol. 2021;22:47. https://doi.org/10.1186/s40360-021-00513-x. PMID: 34462002.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Sauzet O, Carvajal A, Escudero A, Molokhia M, Cornelius VR. Illustration of the Weibull shape parameter signal detection tool using electronic healthcare record data. Drug Saf. 2013;36:995–1006. https://doi.org/10.1007/s40264-013-0061-7. PMID: 23673816.

    Article  PubMed  Google Scholar 

  29. Sasaoka S, Matsui T, Hane Y, Abe J, Ueda N, Motooka Y, et al. Time-to-onset analysis of drug-induced long QT syndrome based on a spontaneous reporting system for adverse drug events. PLoS One. 2016;11:e0164309. https://doi.org/10.1371/journal.pone.0164309. PMID: 27723808.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Nakao S, Hatahira H, Sasaoka S, Hasegawa S, Motooka Y, Ueda N, et al. Evaluation of Drug-Induced Photosensitivity Using the Japanese Adverse Drug Event Report (JADER) Database. Biol Pharm Bull. 2017;40:2158–65. https://doi.org/10.1248/bpb.b17-00561. PMID: 29199239.

    Article  CAS  PubMed  Google Scholar 

  31. van Puijenbroek EP, Egberts AC, Meyboom RH, Leufkens HG. Signalling possible drug–drug interactions in a spontaneous reporting system: delay of withdrawal bleeding during concomitant use of oral contraceptives and itraconazole. Br J Clin Pharmacol. 1999;47:689–93. https://doi.org/10.1046/j.1365-2125.1999.00957.x. PMID: 10383548.

    Article  PubMed  PubMed Central  Google Scholar 

  32. van Puijenbroek EP, Egberts AC, Heerdink ER, Leufkens HG. Detecting drug–drug interactions using a database for spontaneous adverse drug reactions: an example with diuretics and non-steroidal anti-inflammatory drugs. Eur J Clin Pharmacol. 2000;56:733–8. https://doi.org/10.1007/s002280000215. PMID: 11214785.

    Article  PubMed  Google Scholar 

  33. Qian Y, Ye X, Du W, Ren J, Sun Y, Wang H, et al. A computerized system for detecting signals due to drug–drug interactions in spontaneous reporting systems. Br J Clin Pharmacol. 2010;69:67–73. https://doi.org/10.1111/j.1365-2125.2009.03557.x. PMID: 20078614.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Montesi SB, Fisher JH, Martinez FJ, Selman M, Pardo A, Johannson KA. Update in Interstitial Lung Disease 2019. Am J Respir Crit Care Med. 2020;202:500–7. https://doi.org/10.1164/rccm.202002-0360UP. PMID: 32412784.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Cho SJ, Stout-Delgado HW. Aging and Lung Disease. Annu Rev Physiol. 2020;82:433–59. https://doi.org/10.1146/annurev-physiol-021119-034610. PMID: 31730381.

    Article  CAS  PubMed  Google Scholar 

  36. Selman M, Pardo A. Revealing the pathogenic and aging-related mechanisms of the enigmatic idiopathic pulmonary fibrosis. An integral model. Am J Respir Crit Care Med. 2014;189:1161–72. https://doi.org/10.1164/rccm.201312-2221PP. PMID: 24641682.

    Article  CAS  PubMed  Google Scholar 

  37. Katayama K, Yoshino T, Munakata K, Yamaguchi R, Imoto S, Miyano S, et al. Prescription of kampo drugs in the Japanese health care insurance program. Evid Based Complement Alternat Med. 2013;2013:576973. https://doi.org/10.1155/2013/576973. PMID: 24550992.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Matsuno O. Drug-induced interstitial lung disease: mechanisms and best diagnostic approaches. Respir Res. 2012;13:39. https://doi.org/10.1186/1465-9921-13-39. PMID: 22651223.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Liu ZL, Tanaka S, Horigome H, Hirano T, Oka K. Induction of apoptosis in human lung fibroblasts and peripheral lymphocytes in vitro by Sho-saiko-to derived phenolic metabolites. Biol Pharm Bull. 2002;25:37–41. https://doi.org/10.1248/bpb.25.37. PMID: 11824553.

    Article  CAS  PubMed  Google Scholar 

  40. Barkauskas CE, Cronce MJ, Rackley CR, Bowie EJ, Keene DR, Stripp BR, et al. Type 2 alveolar cells are stem cells in adult lung. J Clin Invest. 2013;123:3025–36. https://doi.org/10.1172/JCI68782. PMID: 23921127.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Adamson IY, Hedgecock C, Bowden DH. Epithelial cell-fibroblast interactions in lung injury and repair. Am J Pathol. 1990;137:385–92 PMID: 1696785.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. Antoniou KM, Margaritopoulos GA, Tomassetti S, Bonella F, Costabel U, Poletti V. Interstitial lung disease. Eur Respir Rev. 2014;23:40–54. https://doi.org/10.1183/09059180.00009113. PMID: 24591661.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Burgoyne RA, Fisher AJ, Borthwick LA. The role of epithelial damage in the pulmonary immune response. Cells. 2021;10:2763. https://doi.org/10.3390/cells10102763. PMID: 34685744.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Cai Y, Ma W, Xiao Y, Wu B, Li X, Liu F, et al. High doses of baicalin induces kidney injury and fibrosis through regulating TGF-β/Smad signaling pathway. Toxicol Appl Pharmacol. 2017;333:1–9. https://doi.org/10.1016/j.taap.2017.08.003. PMID: 28803990.

    Article  CAS  PubMed  Google Scholar 

  45. Shii T, Kuroda M, Shamoto N, Mimaki Y. An analysis of the ingredients in decoctions and extracts of Kampo medicines: amounts of baicalin and baicalein in Kampo medicines containing Scutellariae radix. Nihon Ronen Igakkai Zasshi. 2020;57:72–80. https://doi.org/10.3143/geriatrics.57.72. PMID: 32074563.

    Article  PubMed  Google Scholar 

  46. Li X, Li X, Huang N, Liu R, Sun R. A comprehensive review and perspectives on pharmacology and toxicology of saikosaponins. Phytomedicine. 2018;50:73–87. https://doi.org/10.1016/j.phymed.2018.09.174. PMID: 30466994.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Wang Q, Zheng XL, Yang L, Shi F, Gao LB, Zhong YJ, et al. Reactive oxygen species-mediated apoptosis contributes to chemosensitization effect of saikosaponins on cisplatin-induced cytotoxicity in cancer cells. J Exp Clin Cancer Res. 2010;29:159. https://doi.org/10.1186/1756-9966-29-159. PMID: 21143894.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Huang W, Sun R, Zhang Z. “Dose-time-toxicity” relationship study on hepatotoxicity caused by multiple dose of total Bupleurum saponin crude extracts to rats. Zhongguo Zhong Yao Za Zhi. 2010;35:3344–7 Chinese PMID: 21438404.

    PubMed  Google Scholar 

  49. Zhang X, Cai Y, Wang L, Liu H, Wang X. Optimization of processing technology of Rhizoma Pinelliae Praeparatum and its anti-tumor effect. Afr Health Sci. 2015;15:101–6. https://doi.org/10.4314/ahs.v15i1.14. PMID: 25834537.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Zhao X, Zhao J, Hu R, Yao Q, Zhang G, Shen H, et al. Ruanjian Sanjie decoction exhibits antitumor activity by inducing cell apoptosis in breast cancer. Oncol Lett. 2017;13:3071–9. https://doi.org/10.3892/ol.2017.5832. PMID: 28529560.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Zhang ZH, Zhao YY, Cheng XL, Dai Z, Zhou C, Bai X, et al. General toxicity of Pinellia ternata (Thunb.) Berit. in rat: a metabonomic method for profiling of serum metabolic changes. J Ethnopharmacol. 2013;149:303–10. https://doi.org/10.1016/j.jep.2013.06.039. PMID: 23831079.

    Article  PubMed  Google Scholar 

  52. Montastruc JL, Sommet A, Bagheri H, Lapeyre-Mestre M. Benefits and strengths of the disproportionality analysis for identification of adverse drug reactions in a pharmacovigilance database. Br J Clin Pharmacol. 2011;72:905–8. https://doi.org/10.1111/j.1365-2125.2011.04037.x. PMID: 21658092.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Komada F, Nakayama Y, Takara K. Analysis of time-to-onset and onset-pattern of interstitial lung disease after the administration of monoclonal antibody agents. Yakugaku Zasshi. 2018;138:1587–94. https://doi.org/10.1248/yakushi.18-00094. PMID: 30504674.

    Article  CAS  PubMed  Google Scholar 

  54. Komada F. Analysis of time-to-onset of interstitial lung disease after the administration of small molecule molecularly-targeted drugs. Yakugaku Zasshi. 2018;138:229–35. https://doi.org/10.1248/yakushi.17-00194. PMID: 29386436.

    Article  CAS  PubMed  Google Scholar 

  55. Yonemori K, Hirakawa A, Kawachi A, Kinoshita F, Okuma H, Nishikawa T, et al. Drug induced interstitial lung disease in oncology phase I trials. Cancer Sci. 2016;107:1830–6. https://doi.org/10.1111/cas.13087. PMID: 27685762.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Ministry of Health, Labour and Welfare. National Health and Nutrition Surve; 2019. https://www.mhlw.go.jp/bunya/kenkou/kenkou_eiyou_chousa.html. Accessed 15 Jun 2022.

  57. Fujita T, Nagakawa H, Izawa T, Okada T, Tanabe N, Takiguchi Y, et al. Case of Shakuyaku-kanzo-to-induced CD4 dominant pneumonitis diagnosed on day eight of the challenge test. Nihon Kokyuki Gakkai Zasshi. 2008;46:717–21 PMID: 18939414.

    PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This research was partially supported by the Japan Society for the Promotion of Science (JSPS) KAKENHI Grant Number, 21K06646 and 20K10408.

Author information

Authors and Affiliations

Authors

Contributions

All authors have contributed to this scientific work and approved the final version of the manuscript. MN, and OK designed this study, performed the data analyses, and wrote the manuscript. SH, MM, JL, MN involved in methodology and software. MT, KM, RS, MI, YY, WW, SH, MI, TS, MM, SN, and MN assisted the data curation and validation. MI, JL, KI and MN supervised the drafting of the manuscript. All authors took responsibility for the integrity of the data and accuracy of the data analysis.

Corresponding author

Correspondence to Mitsuhiro Nakamura.

Ethics declarations

Ethics approval and consent to participate

Not applicable. Our research does not fall within the purview of any of the following laws and guidelines: “Clinical Trials Act (Act No. 16 of April 14, 2017),” “Act on Securing Quality, Efficacy and Safety of Products Including Pharmaceuticals and Medical Devices (Law number: Act No. 145 of 1960, Last Version: Amendment of Act No. 50 of 2015),” “Guideline for good clinical practice E6 (R1), https://www.pmda.go.jp/int-activities/int-harmony/ich/0076.html,” “Ethical guidelines for human genome and gene analysis research, https://www.mhlw.go.jp/general/seido/kousei/i-kenkyu/genome/0504sisin.html,” and “Ethical Guidelines for Medical and Health Research Involving Human Subjects, https://www.mhlw.go.jp/stf/seisakunitsuite/bunya/hokabunya/kenkyujigyou/i-kenkyu/index.html#HID1_mid1.” Therefore, it is not subject to ethical examination. The study was an observational study without any research subjects. No consent to participate was required due to the retrospective nature of this study.

No administrative permissions or licenses were required to access the raw data from the Japanese Adverse Drug Event Report database because all results were obtained from data openly available online from the PMDA website. All data from the JADER database were fully anonymized by the regulatory authority before we accessed them.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Oura, K., Tanaka, M., Matsumoto, K. et al. Analysis of drug-induced interstitial lung disease caused by herbal medicine using the Japanese Adverse Drug Event Report database. BMC Complement Med Ther 24, 121 (2024). https://doi.org/10.1186/s12906-024-04428-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12906-024-04428-y

Keywords