Introduction
Traditional Chinese medicine compatibility (TCMC) is a basic form of clinical medication of traditional Chinese medicine (TCM) for the treatment and prevention of diseases. Proper TCMC cannot only improve clinical efficacy but also reduce toxicity caused by drug interactions. However, administering incompatible drugs may increase toxicity, create new toxicity, and even negatively affect treatment efficacy [
1]. Traditional Chinese medicine incompatibility (TCMI), which primarily comprises “eighteen incompatible pairs” and “nineteen medicaments of mutual antagonism,” makes great progresses and plays a crucial role in reducing the risk of clinical medication [
2]. However, with the endless number of reports about liver injury related to traditional nontoxic medicines, new challenges emerged, and they should be viewed objectively and managed scientifically. Furthermore, TCMI that affects clinical medication safety should be further studied.
Epimedii Folium (EF) combined with Psoraleae Fructus (PF) is a common TCM preparation used to invigorate the liver and kidney; however, related preparations can lead to liver injury. The Chinese Food and Drug Administration has reported that liver injury can be caused by two Chinese patent medicines, namely, Zhuangguguanjiewan pills (ZGW) and Xianlinggubao capsule (XLG), which both contain EF and PF [
3–
9]. Increasing clinical reports have indicated potential hepatotoxicity in patients who use PF-containing drugs [
10,
11]. One clinical case demonstrated that a PF-induced liver injury can obtain a score of 3 by Roussel Uclaf Causality Assessment Method (RUCAM)-based causality, which may be related to drug metabolism
in vivo and immune responses [
12]. Moreover, with the updated RUCAM, a major step forward has been made to facilitate causality assessment in suspected herb-induced liver injury (HILI) cases; then, the updated RUCAM is improved by providing a better definition of the elements to consider and more accuracy in data and elements to assist the exclusion of alternative causes [
13]. EF could induce liver injury according to clinical research, but no hepatotoxicity was observed in a long-term toxicity experiment [
14]. Additionally, ZGW, which contains both EF and PF, has shown immune idiosyncratic liver injury on the basis of clinical case reports and a nonhepatotoxic dose of lipopolysaccharide (LPS)-mediated rat model of idiosyncratic drug-induced liver injury (IDILI). Nevertheless, the insight into whether EF or PF alone or in combination can induce idiosyncratic liver injury remains to be studied. In this study, we evaluated the liver injury induced by EF or PF alone or in combination and investigated the related mechanism by using the LPS model.
Materials and methods
Chemicals and reagents
Acetonitrile was purchased from Merck (Darmstadt, Germany). Methanol was purchased from Burdick and Jackson (Ulsan, Korea). Formic acid was obtained from Sigma-Aldrich Co. (St. Louis, USA). These three substances were of HPLC grade. Furthermore, double-distilled water was purified using the Millipore water purification system (Millipore, Bedford, MA, USA). LPS and pentobarbital sodium (20151217) were purchased from Sinopharm Chemical Reagent Co., Ltd. LPS was derived from Escherichia coli 055:B5 (Lot #057M4013V). Assay kits to detect plasma alanine aminotransferase (ALT) and aspartate aminotransferase (AST) activities were purchased from Jiancheng Biological Technology, Co., Ltd. (Nanjing, China). All of the other reagents and solvents were of the highest commercially available grade.
Preparation and analysis of test materials
EF and PF were supplied by Sanjiu Medical and Pharmaceutical Co., Ltd. (Batch No. 1512007S; shelf time: 3 years). Epimedium koreanum Nakai and Psoralea corylifolia L. were authenticated by Dr. Xiaohe Xiao, director of Integrative Medical Centre, the Fifth Medical Centre, Chinese PLA General Hospital (Beijing, China). Voucher specimens were deposited in the Herbarium of Sanjiu Medical and Pharmaceutical Co., Ltd., China. In addition, EF and PF were crushed into fine powders, and the quality of these herbs for use in animal studies were tested and verified. Next, EF and PF powders were suspended in 0.5% sodium carboxymethyl cellulose (CMC-Na) for administration.
Animal handling and experimental design
Male Sprague–Dawley rats (180–220 g) were obtained from the Laboratory Animal Center of the Academy of Military Medical Sciences (License No. SCXK2012-004, Beijing, China). The rats were housed in the Laboratory Animal Center of the Fifth Medical Centre, Chinese PLA General Hospital (animal ethics committee approval No. IACUC-2017-003). The room temperature was regulated at 25±2 °C with 50%–60% humidity. A 12 h light/dark cycle was maintained, and the animals were allowed to access food and water ad libitum. All of the animals were acclimatized for 1 week prior to the start of the experiments.
The assessment of idiosyncratic hepatotoxicity was based on our previously reported rat model, which was modified from the literature. On the basis of one of our published articles and the specification of ZGW, we adopted thrice the clinical dose of ZGW (3.78 g/kg) to perform our research [
3]. Then, on the basis of the manufacturing process provided by the manufacturer, 0.22 g/kg of EF and PF were adopted to accomplish our research. The animals were randomly divided into eight separate groups, which were administered with deionized water (control group, Con); EF powder, 0.22 g/kg (EF); PF powder, 0.22 g/kg (PF); EF powder and PF powder, 0.22 g/kg (PE); LPS, 2.8 mg/kg (model group, Mod); EF powder (0.22 g/kg) and LPS (2.8 mg/kg) (EL); PF powder (0.22 g/kg) and LPS (2.8 mg/kg) (PL); or EF powder, PF powder (0.22 g/kg), and LPS (2.8 mg/kg) (PEL). The animals were injected with LPS or normal saline in the tail vein by using standard techniques, followed by the intragastric administration of the aforementioned herbs or an equivalent volume of normal saline 2 h later. Then, 6 h after oral administration, the rats were anesthetized with pentobarbital sodium. Blood samples were collected from the inferior vena cava by using a syringe containing sodium heparin after a midline abdominal incision, and the livers were removed from the rats immediately after sacrifice. The isolated livers were utilized for histopathological examination, and the plasma samples were separated from the collected blood and analyzed for ALT and AST activities. The experimental protocol was approved by the Committee on the Ethics of Animal Experiments of the Fifth Medical Centre, Chinese PLA General Hospital.
Serum biochemistry, histopathological analysis, and immunohistochemistry
After centrifugation (3500 rpm, 10 min, 4 °C), plasma ALT and AST activities were determined according to the microplate assay kit instructions. The left lateral liver lobes obtained were fixed and preserved in 10% neutral-buffered formalin for at least 72 h before being processed for histologic analysis. Paraffin-embedded samples were sectioned at approximately 5 µm and stained with hematoxylin and eosin (HE). TdT-mediated dUTP Nick-End Labeling (TUNEL) staining was performed using the In Situ Cell Death Detection Kit, POD (Roche, Switzerland), following the manufacturer’s instructions. The samples were visualized with a fluorescence microscope (Nikon, Japan) and then imaged.
Serum sample handling
In a polypropylene tube (1.5 mL), 200 µL thawed serum sample and 600 µL methanol were mixed and incubated for 20 min at 4 °C. After the solid debris was removed by centrifugation at 12 000 rpm for 10 min at 4 °C, the supernatants were moved into a polypropylene tube and filtered with a syringe filter (0.22 µm). Finally, the UHPLC–mass spectrometry (MS) system was injected with 4 µL supernatant.
Metabolomics analysis
Metabolic profiling analysis of the biofluids was performed using the Waters Xevo G2-XS QTOF/MS (Waters, Manchester, UK). An Acquity UPLC HSS T3 C18 column analytical column (temperature 30 °C) was injected with 4 µL aliquots of each sample. For positive electrospray ionization source (ESI) analysis, samples were isolated using a 30 min linear gradient of solvent A (water spiked with 0.1% formic acid) and solvent B (acetonitrile spiked with 0.1% formic acid) as mobile phases. The flow rate was fixed at 0.30 mL/min. Then, 10 µL of each sample was drawn as a quality control sample to ensure that the system was stable and the analyses were repeatable. Every 10th sample was injected with the control sample and subsequently analyzed. In analysis, the adopted mobile phases were (A) water and (B) acetonitrile. The same conditions were applied for other analyses. MS was performed with the Waters Xevo G2-XS QTOF/MS instrument (Waters, Manchester, UK) by using an ESI in both positive and negative modes. Progenesis QI (v. 2.0, Waters Technologies, UK) was used to preprocess all of the data. The online software program MetaboAnalyst 3.0 was used to estimate the missing value and to filter and normalize the data. For multivariate statistical analyses, including principal component analysis (PCA) and orthogonal partial least-squares discriminant analysis (OPLS-DA), the resultant data matrices were loaded into SIMCA 13.0 (Umetrics, Umeå, Sweden). All of the variables from the data matrix were mean-centered and scaled to the Pareto variance before PCA. The score plot of the PCA was used to demonstrate the natural interrelationship of the observations. Variables with a high VIP and |p(corr)| value (VIP≥1.0 and |p(corr)|≥0.5) in the OPLS-DA model and with significant differences between groups (P<0.05 and folder change>1.5) were selected as potential biomarkers. All of the selected variables were identified by Agilent MassHunter PCDL Manager software with the KEGG1 database.
Statistical analysis
The data were analyzed using SPSS (v. 17.0, Chicago, IL, USA). Unless otherwise indicated, the data were expressed as the mean±SD. A one-way analysis of variance (ANOVA) was used to statistically analyze the results. The differences were considered significant when P<0.05 and highly significant when P<0.01.
Results
Development of liver injury in LPS-treated rats after EF, PF, or PE administration
The serum biochemistry results are presented in Fig. 1A and 1B. Plasma ALT and AST are well-recognized markers of various types of liver damage and were thereby used for the analysis. As shown in Fig. 1A and 1B, the ALT and AST levels of rats remained unchanged when exposed to low-dose LPS (model group) compared with those of the control group, suggesting that this dose of LPS alone was not associated with liver injury. In addition, no significant differences were found between the AST and ALT levels of rats treated with EF, PF, or PE alone compared with those of the control group. In contrast, the plasma ALT and AST levels of LPS-exposed rats were markedly increased when administered with EF, PF, or PE compared with those of both the control and model groups. Furthermore, the transaminase activity in PEL-treated rats was the highest among all the groups, and the increase in plasma transaminase activity in EL-treated rats was slightly higher than that in PL-treated rats.
Liver tissues were examined by microscopy to obtain visual evidence of the degree of liver damage. As shown in Fig. 1C, the liver sections of the control group showed normal hepatocyte structures. The liver samples from EF-, PF-, and PE-treated rats were almost indistinguishable from those from the normal rats, although the EF-treated rats exhibited minimal steatosis. The samples from the model group displayed slight inflammatory cell infiltration in the portal area but no evident hepatocyte injury. Liver samples from rats co-treated with PF and LPS exhibited hepatocyte focal necrosis, loss of central vein intima, and inflammatory cell infiltration in portal areas. Notably, the liver damage in specimens co-treated with EF and LPS exhibited more severe hepatocyte focal necrosis and extensive inflammatory cell infiltration in portal areas. Moreover, hepatocyte focal necrosis was most severe in the group co-treated with PE and LPS, further showing the highest plasma ALT level among the groups. The drug-induced liver injury observed in the PEL group was more serious than that observed in the other groups, indicating that PF aggravated the EF hepatotoxicity.
The TUNEL assay did not reveal any apparent hepatocyte apoptosis in the control group and model group, with positive TUNEL signals occasionally scattered in the nuclei as small round or circular shapes (Fig. 2). Compared with the model group, the PL group (P<0.05) exhibited both scattered and clustered positive signals, whereas the EL and PEL groups exhibited significant hepatocyte apoptosis (P<0.001). Furthermore, hepatocyte apoptosis was significantly increased in the PEL group compared with that in the EL group, and a large number of positive signals (block distribution) was observed in the PEL group. Thus, compared with treatment with EF or PF alone, treatment with PE led to a significantly increased hepatocyte apoptosis.
EF or PF alone or PE modulated the inflammatory cytokine expression in an LPS model
ELISA was used to determine the TNF-α, IL-1β, IL-6, and IFN-g levels in rat plasma (Fig. 3). Compared with those in the control group, the TNF-α, IL-1β, and IL-6, levels in the model group were significantly increased, thus demonstrating that low-dose LPS injected into the tail vein did not cause liver toxicity but instead stimulated the production of considerably numerous inflammatory cytokines. The PF, EF, and PE groups exhibited no significant difference compared with the control group. Compared with those in the model group, the TNF-α, IL-6, and IFN-g levels were significantly increased in the PL group, whereas plasma TNF-α, IL-1β, IL-6, and IFN-g levels were significantly increased in the EL group. However, the TNF-α, IL-1b, and IFN-g levels increased significantly in the PEL group compared with those in the model group. Thus, EF, PF, or PE enhanced the LPS-stimulated production of inflammatory cytokines, implying that IL-1β is a specific indicator of EF-induced immune idiosyncratic hepatotoxicity.
Metabolomics analysis of EF, PF, or PE in an LPS model
The UHPLC system provides a rapid, effective, and convenient method to analyze the variance in chemical constituents between different rat samples. The base peak intensity chromatograms (BPC) of samples from the control, model, EL, PL, and PEL groups in positive and negative ion modes are presented in Figs. S1 and S2, respectively. Visual inspection of these spectra revealed differences in BPC profiles among the control, model, EL, PL, and PEL groups, indicating that metabolite levels were altered by LPS co-administered with EF, PF, or PE.
As the PEL group exhibited the most severe liver damage, the control, model, EL, PL, and PEL groups were specifically selected for explicit classification. Initially, PCA was used as an unsupervised statistical method to study metabolic differences between the control, model, EL, PL, and PEL groups. The score plots for the PCA analysis derived from ESI+ mode and ESI− mode data are shown in Fig. S3A and S3B, respectively. As shown in Fig. S3, the QC samples clustered closely in both PCA score plots, demonstrating the stability of the LC/MS system throughout the analysis. Furthermore, apparent separation was observed among the control, model, EL, PL, and PEL groups in both PCA models, indicating considerable metabolite differences between the five groups.
OPLS-DA was used for the in-depth investigation of the differences between the control and model groups, the model and PL groups, the model and EL groups, and the model and PEL groups to identify potential biomarkers for discriminating idiosyncratic liver injury and drug types. Fig. S4A–S4J displays the result of an OPLS-DA model derived from ESI+ data analysis. Moreover, data analysis of the ESI− mode was performed, and the results are shown in Fig. S5. Finally, the metabolites from the ESI+ and ESI− mode analyses were combined and subsequently identified by their molecular formulas. All of the biomarkers were tentatively identified with the accurate mass–charge ratio by the online METLIN database. The metabolomics analysis results shown in Figs. 4 and 5 and Table 1 revealed that the levels of biomarkers such as phosphatidic acid (PA), phosphatidylethanolamine, phosphatidylinositol (PI), phosphatidylglycerol (PG), phosphatidylserine (PS), sphingomyelin, and N-acylsphingosine in the PEL, EL, and PL groups were significantly different from those in the control group. The specific trend is shown in Table 1. These biomarkers were mainly enriched in the glycerophospholipid metabolic pathway, the biomarker intensity in the EL group changed, whereas that in the PL group changed more evidently. The specificity of certain chemical compounds in the PEL group was significantly increased, indicating that the LPS treatment combined with EF and PF caused a disordered glycerophospholipid metabolism, which in turn caused liver damage. Additional analysis showed increased levels of PA, sphingomyelin, and N-acylsphingosine, as well as decreased levels of phosphatidylethanolamine, PG, and PS in the PL, EL, and PEL groups, especially in the PL group. Thus, liver damage was potentially caused by the disruption of PA and sphingomyelin metabolism.
Discussion
In recent years, the safety problems of traditional nontoxic Chinese medicine have frequently occurred, and some literature have reported many traditional nontoxic Chinese medicines, such as Dictamni Cortex, Menthae Haplocalycis Herba, PF, Bupleuri Radix, Toosendan Fructus, Puerariae Lobatae Radix, Polygoni Multiflori Radix (PM), and Kava, that can cause HILI [
15–
18]. HILI can be divided into intrinsic or idiosyncratic types. Intrinsic HILI is generally dose or time dependent and predictable with relatively short latency and slight individual differences. However, idiosyncratic HILI is unpredictable, affects only susceptible individuals, is less dose dependent, and has variable latency, presentation, and course [
19]. Therefore, the insight into how to scientifically evaluate idiosyncratic HILI induced by traditional nontoxic Chinese medicines has become one of the most challenging problems in the field of TCM safety research [
20,
21].
RUCAM is the worldwide most commonly used causality assessment method (CAM) for HILI and provides a straightforward application in HILI with scored items. Especially, the updated RUCAM can identify HILI cases early in clinical development, enabling physicians and researchers to take measures to minimize the risk of hepatic reactions [
22–
24]. In recent years, we have established an integrated chain of evidence for the diagnosis of liver injury induced by Chinese herbal medicines; the integrated evidence-chain method (iEC) is then proposed [
25]. On the basis of this method, we demonstrated that the objectivity of liver injury was caused by PM, and the idiosyncratic liver injury of PM may be immune mediated, as proven through an IDILI model [
19]. We also have verified that ZGW, which contains both EF and PF, has a property of immune idiosyncratic hepatotoxicity according to iEC and a nonhepatotoxic dose of LPS model [
3]. PF may be the main substance in ZGW that induced idiosyncratic liver injury [
5]. Furthermore, Teschke used the structured causality assessment to indicate that Ayurvedic herb PF is potentially hepatotoxic [
26]. In another study, oral PF administered at 0.4 g versus 1.6 g (high dose) daily for 12 weeks in rats caused liver injury in the form of steatosis and hepatocyte necrosis [
27]. Meanwhile, Bavachin, which is one of the flavonoids in PF, is a hepatotoxic substance [
28]. Taken together, these results showed that PF can trigger a series of toxic events, including liver injury. However, the reason why minority individuals are more prone to the adverse reactions remains unanswered. In addition, EF in ZGW also has certain hepatotoxicity. Cheng verified that mice fed with EF for 3 days exhibited vomiting, poor appetite, and decreased activity, and hepatic steatosis was observed after 15 days of administration [
29]. In chronic toxicity study, the rats were given 40, 80, and 160 times of human clinical dosage of total flavonoids of EF for 12 weeks, suggesting that the chronic toxicity of this drug was also weak [
14]. To sum up, the insight into whether EF can induce liver injury remains controversial
As mentioned, ZGW, which contains EF and PF, can induce immune idiosyncratic hepatotoxicity. Therefore, we used an LPS model to clarify the reasonable or unreasonable compatibility of EF and PF. This study demonstrated that EF, PF, and PE combined with LPS not only increased plasma ALT and AST levels but also caused liver injury. ALT and AST are generally used as sensitive evaluation indexes of liver injury. In fact, liver injury often indicates degeneration or necrosis of hepatocytes, resulting in the leakage of cytosolic enzymes into the blood stream and the increase of serum enzymes. Moreover, ALT mainly exists in the cytoplasm of hepatocytes, but AST is mainly distributed in the mitochondria and cytoplasm of hepatocytes [
30–
32]. Thus, plasma ALT activity increased significantly in the EL, PL, and PEL groups, but no increase was observed in the remaining groups (Fig. 1). To further prove the result, we examined liver histology and TUNEL staining (Figs. 1 and 2) and found that TUNEL staining (marker for apoptosis) increased in rats treated with PL, EL, and PEL, indicating an increase in cell death caused by these combinations. Similar results were found in the liver histology evaluations. Hence, a minor inflammatory response triggered by a nonhepatotoxic dose of LPS decreased the hepatotoxicity threshold for EF and PF, leading to an increase in inflammatory cell infiltration and liver cell apoptosis. Notably, EF or PF did not cause liver injury in normal rats, whereas EF or PF could induce liver injury in LPS-treated rats. EF combined with PF can induce more severe liver injury in LPS-treated rats than the treatments administered separately. Therefore, EF combined with PF can cause serious adverse reactions during inflammation. Thus, this combination can be categorized as immune idiosyncratic hepatotoxicity.
Our inflammatory cytokine data revealed that TNF-α, IL-6, IL-1β, and IFN-
g positively correlated with EL-, PL-, and PEL-induced liver injury [
33–
35]. LPS interacted with cell surface receptors such as CD14 and TLR4, to induce NF-
kB pathway activation and stimulate monocytes, macrophages, neutrophils, lymphocytes, and other cell types to release considerably numerous cellular factors and inflammatory mediators, such as IL-6 and TNF-α, generating a series of inflammatory reactions and liver injuries [
36–
39]. Our results revealed that TNF-α levels gradually increased, following the trend PL<EL<PEL (Fig. 3A). The most vital consequence of the idiosyncratic hepatotoxicity lethal stress is the promotion and amplification of inflammation through activating the transcription of cytokines and chemokines such as TNF-α, IL-6, and IL-1β in IDILI [
40]. IDILI-associated drugs sensitize hepatocytes to cell death signaling from cytokines such as TNF-α and IFN-
g [
41]. IL-1
b, which is processed and released by NLRP3 inflammasome activation, is the most potent proinflammatory cytokine, and it is implicated as the effector molecule in many NLRP3-driven diseases, including IDILI [
42–
47]. Our study revealed that EF might significantly increase IL-1β secretion and activation compared with PF in LPS-treated rats. Thus, IL-1β is a marker for EF-related immune idiosyncratic hepatotoxicity. The results also showed that EF may induce immune idiosyncratic hepatotoxicity mainly by regulating LPS-mediated NLRP3 inflammasome activation. In summary, the activation of either receptor initiates signaling pathways that lead to NF-
kB or NLRP3 inflammasome activation, resulting in proinflammatory cytokine production. Thus, an episode of inflammation during EF treatment predisposes animals to induce liver injury, suggesting that EF combined with LPS is an important factor for the occurrence of idiosyncratic drug-induced liver injury.
The metabolomics analysis showed that PF affects more metabolites in the glycerophospholipid and sphingolipid metabolic pathways compared with EF in an LPS model. Specifically, PA (20:5(5Z,8Z,11Z,14Z,17Z)/13:0); PG (22:4(7Z,10Z,13Z,16Z)/0:0); PS (20:3(8Z,11Z,14Z)/0:0); Sphingosine-1-phosphate; Cer (t18:0/20:0); PI (18:3(9Z,12Z,15Z)/0:0); and PG (15:0/14:0) were greatly affected on the PL group. However, metabonomics analysis showed that the mechanism of liver injury induced by PF is mainly metabolic dysregulation, but EL and PEL groups also affect metabolism. Furthermore, EL, PL, and PEL groups might be regulated by TNF-α, with a gradually increasing trend (PL group<EL group<PEL group) observed for the sphingolipid pathway. TNF-α affects sphingolipid biosynthesis via sphingomyelin phosphodiesterase 2 in the liver [
48]. Under normal circumstances, upregulated N-acylsphingosine levels increase sphingosine-1-phosphate levels [
49,
50]. However, the levels of sphingosine-1-phosphate and its metabolite phosphatidylethanolamine both decreased in our study, potentially causing inhibition and indicating that the PEL, EL, and PL groups inhibited the metabolic conversion of N-acylsphingosine to sphingosine-1-phosphate. Of the two key enzymes in this process, namely, ceramidase and sphingosine kinase, sphingosine kinase is activated by TNF-α to regulate PA levels. However, both TNF-α and PA levels increased in our study; thus, we speculated that ceramidase but not sphingosine kinase was greatly affected. Ceramidase activity is closely related to hepatocyte apoptosis; although ceramidase mainly hydrolyzes N-acylsphingosine, ceramidase inhibition inevitably leads to the accumulation of upstream N-acylsphingosine and affects the metabolism of sphingomyelin [
51,
52]. Upstream sphingomyelin and PA metabolic disorders caused the abnormal metabolism of downstream glycerophospholipid metabolites such as phosphatidylethanolamine, PC, and PI. PC and phosphatidylethanolamine are sources of arachidonic acid (AA) and can be hydrolyzed by phospholipase to produce AA, contributing to inflammation by releasing factors such as prostaglandins and leukotrienes [
53–
58]. PIs are hydrolyzed by phosphatidylinositol 3-kinase (PI3K), which plays a key regulatory role in the occurrence and development of inflammation [
59–
61]. In fact, PF increases the responsiveness of the liver to LPS or other inflammatory mediators via the modulation of multiple metabolic pathways, thereby leading to hepatotoxicity. Therefore, PF may induce liver injury primarily by causing glycerophospholipid and sphingolipid metabolic dysfunctions under immunological stress conditions.
Overall, the combination of EF and PF can lead to a serious liver injury under immunological stress conditions. Furthermore, the main mechanism of liver injury induced by EF combined with PF may regulate immune inflammation and metabolic dysfunctions (Fig. 6), suggesting that immune and metabolic factors of the body should also be comprehensively considered when assessing the safety of EF combined with PF. Thus, EF combined with PF can be defined as a new incompatibility of TCM, thereby enriching our knowledge of the 18 incompatible pairs, one of the traditional theories of TCMI.
Higher Education Press and Springer-Verlag GmbH, Germany, part of Springer Nature