Impact Factor 5.810 | CiteScore 6.2
More on impact ›


Front. Pharmacol., 01 March 2021 |

Risk of Liver Injury Associated with Intravenous Lipid Emulsions: A Prescription Sequence Symmetry Analysis

www.frontiersin.orgXiao-xiao Li1,2, www.frontiersin.orgYin-chu Cheng1,3, www.frontiersin.orgSuo-di Zhai1,2, www.frontiersin.orgPeng Yao2, www.frontiersin.orgSi-yan Zhan3,4* and www.frontiersin.orgLu-wen Shi2*
  • 1Department of Pharmacy, Peking University Third Hospital, Beijing, China
  • 2Department of Pharmacy Administration and Clinical Pharmacy, School of Pharmaceutical Sciences, Peking University, Beijing, China
  • 3Department of Epidemiology and Biostatistics, School of Public Health, Peking University Health Science Center, Beijing, China
  • 4Research Center of Clinical Epidemiology, Peking University Third Hospital, Beijing, China

Aims: To determine the risk of liver injury associated with the use of different intravenous lipid emulsions (LEs) in large populations in a real-world setting in China.

Methods: A prescription sequence symmetry analysis was performed using data from 2015 Chinese Basic Health Insurance for Urban Employees. Patients newly prescribed both intravenous LEs and hepatic protectors within time windows of 7, 14, 28, 42, and 60 days of each other were included. The washout period was set to one month according to the waiting-time distribution. After adjusting prescribing time trends, we quantify the deviation from symmetry of patients initiating LEs first and those initiating hepatic protectors first, by calculating adjusted sequence ratios (ASRs) and relevant 95% confidence intervals. Analyses were further stratified by age, gender, and different generations of LEs developed.

Results: In total, 416, 997, 1,697, 2,072, and 2,342 patients filled their first prescriptions with both drugs within 7, 14, 28, 42, and 60 days, respectively. Significantly increased risks of liver injury were found across all time windows, and the strongest effect was observed in the first 2 weeks [ASR 6.97 (5.77–8.42) ∼ 7.87 (6.04–10.61)] in overall patients. In subgroup analyses, female gender, age more than 60 years, and soybean oil-based and alternative-LEs showed higher ASRs in almost all time windows. Specially, a lower risk for liver injury was observed in the first 14 days following FO-LEs administration (ASR, 3.42; 95% CI, 0.81–14.47), but the risk started to rise in longer time windows.

Conclusion: A strong association was found between LEs use and liver injury through prescription sequence symmetry analysis in a real-world setting, which aligns with trial evidence and clinical experience. Differences revealed in the risks of liver injury among various LEs need further evaluation.


•Drug safety signal was detected by using prescription sequence symmetry analysis and Chinese Basic Health Insurance database.

•Lipid emulsion use in the first 2 months is associated with 3.6- to 7.9-fold increased risk for liver injury requiring use of hepatic protectors in a real-world setting.

•The association is positive across all time windows and generations.

•It suggests that liver injury after lipid emulsion initiation is common in China.

•It is important to strengthen the appropriate use of LEs at the beginning of administration.


Patients dependent on total or partial parenteral nutrition (PN) are at higher risk of developing a wide range of disruption to liver function, such as cholestasis and steatosis (Beath and Kelly, 2016; Meyerson and Naini, 2019). It occurs as a major consequence of metabolic complications related to PN, besides physiological or anatomical abnormalities, with an incidence ranging from 20% to 80% in both adults and children (Gabe, 2013).

As an integral component of PN, a wide variety of commercial lipid emulsions (LEs) is now available. However, the amount, type, and infusion time of intravenous LEs were reported to affect the risk of inducing liver complications (Badia-Tahull et al., 2015; Kapoor et al., 2019a; b). Exposure to excess intravenous LE >1 g/kg/d in adults may result in hepatic steatosis, which is considered as the first step of liver injury (Javid et al., 2005; Kumar and Teckman, 2015; Beath and Kelly, 2016). Compared to LEs being part of the total nutrient admixture, LEs separately taken are prone to medication errors (MEs), especially when the prescription is not under the supervision of a nutrition support pharmacist (Berlana et al., 2019). To our knowledge, separate use is, however, the most common way of LE prescriptions in China (58.01% of patients received PN in a cross-sectional survey among eight tertiary hospitals during 2011–2014) (The Pharmacy Workgroup of CSPEN, unpublished data), which adds to the risk of liver injury. Besides, owing to excessive polyunsaturated fatty acid and linoleic acid content that might increase lipid peroxidation and inflammatory response, soybean oil-based LEs (S-LEs), though providing enough energy and essential fatty acids, have a tendency to cause cellular damage and liver injury (Carter et al., 2007; Wanten and Calder, 2007; Beath and Kelly, 2016). Current evidences from murine models suggest that new generations of LEs, including alternative-LEs (such as medium chain triglycerides-soybean oil LE (MCT-LEs) and olive-soybean oil LE (O-LE)) and fish oil-based LEs (FO-LEs), may improve biochemical measures of hepatobiliary function (Le et al., 2012; Baker et al., 2019). However, data and evidence are still limited especially from routine clinical practice and high quality clinical studies (Javid et al., 2005; Kumar and Teckman, 2015; Kapoor et al., 2019a; b). Therefore, it is important to study the relationship between LEs and liver injury in a real-world setting. To the best of our knowledge, no large population-based studies have been conducted in adults to quantify the risk of hepatic dysfunction following the initiation of LEs.

Prescription sequence symmetry analysis (PSSA) is a valid method used for rapid signal detection of adverse drug events (ADE) by calculating the sequence ratio between exposure and outcomes (Hallas, 1996; Tsiropoulos et al., 2009), which inherently controls time-constant confounders. For now, only two Chinese studies reported drug-related liver injury using PSSA (Fang et al., 2016; Zhang et al., 2019). Moreover, the health insurance data is an important source used for safety evaluation (Tyree et al., 2006). This study aimed to estimate the safety signal between LEs and hepatic protectors using PSSA in a Chinese health insurance database.

Materials and Methods

Data Source

This study was conducted by analyzing data from the 2015 Chinese Health Insurance Research Association (CHIRA) database, which is a national-level claims database collecting sampled hospital record of patients from the Urban Employee Basic Medical Insurance scheme all over mainland China (Yang et al., 2018). The data was annually resampled with a two-stage sampling design which has been described previously (Xia et al., 2015; Yong et al., 2018). The 2015 CHIRA database contained a total of 4.64 million patients selected from four municipalities, all provincial capitals, one prefecture-level city of each province, and two county-level cities from each province. The sample proportion was 2%, 5%, and 10%, respectively, of patients from municipalities and provincial capitals, prefecture-level cities, and county-level cities. The database maintains detailed hospital record information of both inpatients and outpatients, including demographics, clinical diagnoses, and prescriptions for drugs and procedures. This study was approved by the Ethical Review Board of Peking University Health Science Center, and informed consent of participant was exempted (IRB00001052-15045). The database was de-identified for protection of patients’ privacy.

Study Design

We performed an observational study using PSSA to explore the association between LEs and liver injury. PSSA, first proposed by Hallas (1996), is frequently used as a post-marketing active surveillance tool to detect drug safety signals from large prescription databases, with verified validity (Wahab et al., 2013; Pratt et al., 2014). The principles of PSSA have been fully discussed previously (Lai et al., 2017). In brief, it examines the sequence of a marker drug (used to treat an ADE) initiated before and after an index drug (suspected of inducing an ADE): in the absence of a causal association, the propensity to initiate the index drug before or after the marker drug will theoretically be equal in the patient population, showing a symmetrical pattern; however, if the index drug causes the ADE requiring the marker drug for treatment, there will be an asymmetrical prescribing pattern in which more patients will be observed of initiating the index drug before the marker drug. Since PSSA is based on a self-controlled case-only new user design, it is unaffected by potential between-subject confounding and is robust toward time-invariant confounders such as race, gender, and genetic characteristics (Lai et al., 2014).

Cohort Selection

The index drugs in our study are LEs, consisting of S-LEs, alternative-Les, and FO-LEs (Table 1). The marker drugs are a group of hepatic protectors widely used in China with high specificity as the surrogate of liver injury, chosen by literature search, guideline recommendation, and clinical experience (expert consultation) (Fang et al., 2016; Yu et al., 2017; Zhang et al., 2019), including anti-inflammatory agents, antioxidant agents, antidote agent, choleretics, cell membrane repair agents, and others (Table 1).


TABLE 1. Lipid emulsions and hepatoprotective drugs in the study.

Patients prescribing both an index drug and a marker drug between Jan 1st, 2015, to Dec 31st, 2015, were identified from the 2015 CHIRA database. We applied a washout period based on the “waiting-time distribution” to ensure that patients were new users of the index and marker drugs (Hallas et al., 1997). The observation periods between the initiation of index drugs and marker drugs, namely, time windows, were restricted to 7, 14, 28, 42, and 60 days, respectively, for the sensitivity analysis in consideration of biologic processes leading to steatosis, steatohepatitis, cholestasis, and fibrosis being very likely to occur within 60 days after initiating LEs (Gabe, 2013; Beath and Kelly, 2016) and to reduce within-subject time-variant confounding (Lai et al., 2014). Patients with first prescriptions of the index drug and the marker drug on the same date were excluded because we were not able to determine which drug was prescribed earlier from the data. This helped to lessen patient misclassification induced by prophylactic hepatic protectors. Patients who used more than one generation of LEs (switchers) were also excluded to avoid potential bias from the existing use of the medicine class (Lai et al., 2017). The whole patient selection process and criteria are detailed in Figure 1. For patients included in the final cohort, we extracted their patient unique identification, age, gender, drug prescriptions (generic names with Anatomical Therapeutic Chemical classification codes), and time of prescription.


FIGURE 1. Flowchart of study cohort selection. CHIRA, Chinese Health Insurance Research Association; LEs, lipid emulsions; PSSA, prescription sequence symmetry analysis.

Statistical Analysis

The parameter to measure the association between the index drug and the suspected AE in PSSA is called sequence ratio. A crude sequence ratio (CSR) can be calculated by dividing the number of patients initiating LEs first (causal group) by that of those initiating hepatic protectors first (non-causal group). Drug prescribing trends over time can extraneously affect the treatment sequence, and it may result in a biased effect estimate. To adjust for this, we further calculated adjusted sequence ratio (ASR) using a correction method proposed by Tsiropoulos et al., (2009). First, a null-effect sequence ratio (NESR), which represents the expected sequence ratio given no causal association between the index drug and the marker drug, is derived from the overall average probability (p) that the index drug will be prescribed before the marker drug in the background population, where p is calculated as


Here, m and n indicate the consecutive days of the study period, u is the last day of the study period, d indicates the specified observation time window, LEm indicates the number of patients receiving first LE prescription on date m, and Hn is the number of patients initiating hepatic protectors on date n. Given p, NESR can be generated as p/(1 − P). ASR can be calculated as CSR/NESR. The estimation of confidence interval (CI) is based on binomial distribution, using the Wilson (Score) method (Morris and Gardner, 1988; Newcombe, 1998), calculated as


Here, n is the sample size included in the final analysis and p^ is the ratio of causal group patient number over n. Detailed computational formulas have been clearly described in previous studies (Tsiropoulos et al., 2009; Adimadhyam et al., 2019).

Subgroup analyses were conducted by LE generation (S-LEs, alternative-LEs, and FO-LEs), age (≥60 and <60 years), and gender. We performed sensitivity analyses in different time windows to test the robustness of PSSA results and find out during which time period the adverse effects were more likely to occur. ASRs with the lower limit of 95% confidence interval (CI) bigger than 1 were considered statistically significant. All analyses were performed using SAS (version 9.4, SAS Institute Inc., Cary, NC, USA).


Figure 1 shows the flowchart of study cohort selection. From the 4.64 million participants covered by the 2015 CHIRA database, we identified 2,342 patients who initiated both a LE and a hepatic protector and met the selection criteria. The waiting-time distributions of both index and marker drugs prescribed in the background population showed a rapid decrease in the first month and reached an almost stable plateau thereafter (Figure 2). Hence, we could infer that many of the patients who filled their initial prescriptions during the first month were prevalence users and should be excluded from our study. Accordingly, the washout period used to exclude prevalent users was set to 1 month. Based on the prespecified time windows, the eligible patients were further grouped into those who started both drugs within 7 (n = 416), 14 (n = 997), 28 (n = 1,697), and 42 (n = 2072) days. Of the 2,342 patients with a maximum time window of 60 days, 35.1% were female. The average age was 62.6 ± 15.0 (standard deviation) years. Alternative-LEs were the most commonly used LE in these patients (62.7%), followed by S-LEs (36.5%). FO-LEs were seldom used (0.73%).


FIGURE 2. Waiting-time distributions of lipid emulsions and hepatic protectors. For each consecutive month during January to December 2015, the graph depicts the number of persons who presented their first recorded prescription of index (left Y-axis) and marker (right Y-axis) drug, respectively, during that month.

In the analysis where index and marker drugs were initiated within 60 days, the adjusted prescribing trends resulted in an ASR of 3.60 (95% CI 3.26–3.97), indicating that initiating a LE is associated with a 3.6-fold increase in the rate of liver injury receiving hepatic protectors during a 60-day period after initiating LEs. Moreover, a strong asymmetrical pattern of treatment sequence was revealed in Figure 3, showing much more incident hepatoprotective therapies in the months following LEs initiation than before initiation of LEs. An increased risk for hepatic dysfunction following intravenous infusion with LEs was also observed in other prespecified time windows, which was the highest during the first week (ASR, 7.87; 95% CI, 5.80–10.68), and fell with time (Table 2, Figure 4).


FIGURE 3. Incident use of hepatic protectors before and after initiation of lipid emulsion within a 60-day time window.


TABLE 2. Crude and adjusted sequence ratios, overall and by subgroups.


FIGURE 4. Adjusted sequence ratios. (A) Overall, (B) by type, (C) by gender, and (D) by age.

In subgroup analyses by the type of LEs initiated (Table 2, Figure 4), all the three generations of LEs showed significantly increased risk of inducing liver injury requiring use of hepatic protectors in all time periods, except for FO-LEs during the first 14 days (ASR, 3.42; 95% CI, 0.81–14.47). ASRs of alternative-LEs were slightly higher than that of S-LEs across all prespecified time windows, and the risks of these two generations decreased over time. ASRs of FO-LEs were much lower than that of S-LEs and alternative-LEs in the first 14 days, but the risk started to rise in longer time windows.

The associations of LEs and hepatic dysfunction were also significant in different age and gender groups, but ASRs differed slightly between groups (Table 2, Figure 4). Generally, ASRs were higher in women and elder patients than men and younger patients across all analysis time windows.


In this large, population-based observational study, we found an association between the index drugs (LEs) and the marker drugs (hepatic protectors), with positive and robust ASRs across all time windows and different generations of LEs, suggesting that LEs might possibly induce hepatic dysfunction in Chinese patients. There was a 3.6- to 7.9-fold increase in the risk of hepatotoxicity within the 2 months after initiating LEs. In 1971, Peden et al. was the first to report a case of an infant who had received total PN for 2.5 months before dying from liver failure (Peden et al., 1971). In a prospective cohort conducted in 90 adults on home parenteral nutrition, 58 patients (65%) developed chronic cholestasis after 6 months of follow-up (Cavicchi et al., 2000). Evidences from different studies also reported that one or more factors relating to hepatic dysfunction, including yet not limited to PN (Lal et al., 2018). Generally, increased intestinal permeability combined with administration of PN promotes lipopolysaccharide and toll-like receptor 4 signaling dependent Kupffer cell activation, which is an early event in the pathogenesis of PN-related liver injury (El Kasmi et al., 2012).

Drug-induced liver injury accounts for approximately 20% of inpatients with acute liver injury in China, owing to a wide use of traditional Chinese medicine and anti-tuberculosis drugs (Li et al., 2007; Yu et al., 2017; Fang et al., 2016; Zhang et al., 2019). The marker drugs used in both studies, and also in our study, were recommended by the Chinese Medical Association guidelines for the treatment of DILI (Yu et al., 2017), and they were widely used to treat liver injury in clinical practice (Fang et al., 2016; Zhang et al., 2019). However, as a signal detection method, PSSA is vulnerable to bias and confounding existing in real-world data (Davies et al., 2019). On one hand, a priori knowledge of the drug-event association may affect the validity of PSSA results, mainly in two ways: 1) the high ASRs observed in our study might be overestimated because prescribers with the knowledge of the ADE might avoid prescribing LEs when the patient already had liver injury; 2) when prescribers with the knowledge of the ADE encounter LE users with liver injury, they might also respond by discontinuing LEs rather than prescribing a liver protection agent, which, on the contrary, would attenuate the signal in our study. On the other hand, fat overload syndrome, which could be caused by rapid infusion (Hojsak and Kolacek, 2014) or high dose (Levit et al., 2016) of LEs, has one clinical manifestation of jaundice and thus might be mistaken as liver injury. This could also lead to an overestimated ASR.

Similar increased risks for hepatic dysfunction showed up in the first 2 months following all the three generations of LEs. However, it is worth noting that a lower risk for hepatic dysfunction was observed in the first 2 weeks following FO-LEs administration, and the absence or reduction of plant sterols in FO-LEs might account for it. Being an important constituent in S-LEs and alternative-LEs, stigmasterol had a prominent role in promoting cholestasis, liver injury, and liver macrophage activation. It might be mediated through suppressing canalicular bile transporter expression (Abcd11/BSEP, Abcd2/MRP2) via antagonism of the nuclear receptors (Fxr, Lxr) and failure of upregulation of the hepatic sterol exporters (Abcg5/g8/ABCG5/8)(El Kasmi et al., 2013). Studies in murine models also showed a benefit for liver diseases from FO-LEs (Le et al., 2012; El Kasmi et al., 2013; Baker et al., 2019). As a small sample size in the FO-LEs subgroup, this result was unstable and showed an increased trend after the first 2 weeks. This was in line with the findings of two meta-analyses focusing on the infants, suggesting that FO-LEs, compared with SO-LEs and alternative-LEs, might have an uncertain role in liver complications (Kapoor et al., 2019a; b). A randomized, controlled clinical trial also showed a comparable influence on hepatic function following different generations of LEs in adults with chronic intestinal failure during 1-year follow-up (Hojsak and Kolacek, 2014).

Being opposite to the increasing time trend of risk for liver injury in FO-LEs, the ASRs were highest in the first 14 days after S-LEs or alternative-LEs initiation and decreased thereafter with longer time windows. This suggested that the majority of patients suffered from acute hepatotoxicity after S-LEs or alternative-LEs initiation, possibly resulting from exposure to intravenous LE >1 g/kg/d caused by inappropriate prescription or MEs (rapid infusion speed) (Beath and Kelly, 2016; Hojsak and Kolacek, 2014; Levit et al., 2016). Some patients might discontinue their LEs after first 14-day prescription and were consequently no longer at risk thereafter, also leading to the decreasing trend of ASRs. For most patients, continuous use of LEs might not be necessary, but for some patients, especially those with short bowel syndrome, chronic intestinal pseudo-obstruction, and radiation enteritis, LEs are essential nutrition agents for life support (Matarese et al., 2005; Pironi et al., 2016). Though the increasing trend of ASRs in FO-LEs was in line with the risk from long-term PN reported previously (Javid et al., 2005; Kumar and Teckman, 2015; Beath and Kelly, 2016), the small sample size in the LE subgroup limited its interpretation.

Moreover, the magnitude of risk for liver injury following LEs initiation was higher among female patients and those with elder age. We noted an increased risk of liver injury in the female group, which was also found in a retrospective study conducted in a New York hospital (Huard et al., 2018). The multivariate analysis revealed that female gender was a significant predictor of advanced liver fibrosis among patients receiving PN, who required intestinal transplantation (Huard et al., 2018). Furthermore, altered expression of hepatic β-adrenergic receptors in relation to age-related lipid metabolic dysfunction in liver may explain a slightly higher risk of hepatotoxicity induced by LEs in patients aged more than 60 years (Shi et al., 2018).

Our study has several strengths. To our knowledge, it is the first to evaluate the association of LEs with the potential risk of hepatic dysfunction in a large Chinese population. It provides information on LE-relevant liver injury in a real-world setting that could help us better understand the current situation of safety problems in clinical practice of PN. PSSA has been proved an effective and fast signal detection method in drug safety evaluation, with moderate sensitivity and high specificity (Maclure et al., 2012; Wahab et al., 2013; Pratt et al., 2015). In the rofecoxib and rosiglitazone case, sequence symmetry analysis detected ADE signals earlier than disproportionality analysis (Izyan et al., 2014). As a simple form of a self-controlled design, PSSA is insensitive to between-subject and time-constant confounders (Kubota, 2016). To minimize residual confounding, we also calculated NESR to adjust for the time trends of prescribing. A validation study showed the consistency of PSSA results for detecting ADE regardless of different patterns of medicine utilization and different health settings (Pratt et al., 2015). We further applied a sensitivity analysis with different time windows to test the robustness of the results. Additionally, only drugs used specifically to treat liver injury were selected as marker drugs in our study, to minimize the potential of outcome misclassification.

Our study also has some limitations. Firstly, we only included one-year patient data because the maximum follow-up duration of CHIRA database was one year as a result of the annually resampling data collection strategy. This limited the sample size of study population, especially in the subgroup of FO-LEs, whose effect estimates would thus be more easily affected by random fluctuation. It should also be noted that some delayed events, such as cholestasis, can occur years after PN initiation in adults (Gabe, 2013). However, limiting the study period has the advantage of minimizing the potential of introducing time-dependent confounding factors in PSSA studies (Lai et al., 2017). Moreover, as most hepatic dysfunctions are generally acute events (4–8 weeks) (Palova et al., 2008; Badia-Tahull et al., 2015), one-year follow-up would still be enough to capture related events of interest. Secondly, although a comprehensive set of marker drugs were chosen as a surrogate for hepatotoxicity, we still could not ascertain that it was hepatic dysfunction that provoked the use of these marker drugs in all participants. A cross-validation was not possible because our data did not hold indication of adverse drug reactions, but it would be a valuable topic in future studies. Thirdly, our claims data could not measure over-the-counter use of marker drugs, though most patients indeed got their prescriptions during the hospital stay. Fourth, though the study covers a representative sample from all over China, results should still be generalized with caution as all patient data is from the 2015 CHIRA database, which might be limited to reflect the current prescribing patterns. Fifth, we were not able to explore the dose-response relationship between LEs and hepatotoxicity because of incomplete dosage information.


Our results show that there is a strong association between LEs and hepatotoxicity with an asymmetrically distributed treatment sequence. The findings suggest that hepatic dysfunction after LEs is common in China, and it is important to strengthen the appropriate use of LEs and enhance patient education at the initiation of LEs to reduce the liver injury.

Data Availability Statement

The data analyzed in this study is subject to the following licenses/restrictions: The dataset presented in this article is not readily available because of confidentiality agreement between Chinese Health Insurance Research Association database and Peking University. Requests to access these datasets should be directed to

The code supporting the conclusions of this article will be made available by the authors. Please send a request to the principal investigator of this article. The steering committee of this study will discuss all requests and decide on the basis of the scientific rigor of the proposal whether code sharing is appropriate. All applicants are asked to sign a code access agreement.

Ethics Statement

This study was approved by the Ethical Review Board of Peking University Health Science Center (IRB00001052-15045). Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

All authors took part in the final version for submission and accepted overall accountability for accuracy and integrity of the manuscript. X-xL was responsible for conceptualization, investigation, and writing-original draft; Y-cC for conceptualization, data curation, formal analysis, visualization, funding acquisition, and writing-original draft; S-dZ for project administration; PY for writing-review and editing; S-yZ for supervision, validation, and methodology; L-wS for supervision, validation, and funding acquisition.


This work was supported by the National Natural Science Foundation of China (Grant Number 71904002, 2020) and the National Science and Technology Major Project (Grant Number 2018ZX09721003-001-002, 2018).

Conflict of Interest

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


ADE, adverse drug event; ASR, adjusted sequence ratio; ATC, anatomical therapeutic chemical; B-LE, borage oil lipid emulsion; CHIRA, Chinese Health Insurance Research Association; CI, confidence intervals; CSR, crude sequence ratio; FO-LE, fish oil-based lipid emulsion; LE, lipid emulsion; MCT-LE, medium chain triglycerides-soybean oil lipid emulsion; ME, medication error; NESR, null-effect sequence ratio; O-LE, olive-soybean oil lipid emulsion; PN, parenteral nutrition; PSSA, prescription sequence symmetry analysis; S-LE, soybean oil-based lipid emulsion.


Adimadhyam, S., Schumock, G. T., Calip, G. S., Smith Marsh, D. E., Layden, B. T., and Lee, T. A. (2019). Increased risk of mycotic infections associated with sodium-glucose co-transporter 2 inhibitors: a prescription sequence symmetry analysis. Br. J. Clin. Pharmacol 85 (1), 160–168. doi:10.1111/bcp.13782

PubMed Abstract | CrossRef Full Text | Google Scholar

Badia-Tahull, M. B., Llop-Talaveron, J., and Leiva-Badosa, E. (2015). Impact of intravenous lipid emulsions on liver function tests: contribution of parenteral fish oil. Nutrition 31 (9), 1109–1116. doi:10.1016/j.nut.2015.04.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Baker, M. A., Cho, B. S., Anez-Bustillos, L., Dao, D. T., Pan, A., O'Loughlin, A. A., et al. (2019). Fish oil-based injectable lipid emulsions containing medium-chain triglycerides or added alpha-tocopherol offer anti-inflammatory benefits in a murine model of parenteral nutrition-induced liver injury. Am. J. Clin. Nutr 109 (4), 1038–1050. doi:10.1093/ajcn/nqy370

PubMed Abstract | CrossRef Full Text | Google Scholar

Beath, S. V., and Kelly, D. A. (2016). Total parenteral nutrition-induced cholestasis: prevention and management. Clin. Liver Dis 20 (1), 159–176. doi:10.1016/j.cld.2015.08.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Berlana, D., Almendral, M. A., Abad, M. R., Fernandez, A., Torralba, A., Cervera-Peris, M., et al. (2019). Cost, time, and error assessment during preparation of parenteral nutrition: multichamber bags versus hospital-compounded bags. JPEN - J. Parenter. Enter. Nutr 43 (4), 557–565. doi:10.1002/jpen.1436

CrossRef Full Text | Google Scholar

Carter, B. A., Taylor, O. A., Prendergast, D. R., Zimmerman, T. L., Von Furstenberg, R., Moore, D. D., et al. (2007). Stigmasterol, a soy lipid-derived phytosterol, is an antagonist of the bile acid nuclear receptor FXR. Pediatr. Res 62 (3), 301–306. doi:10.1203/PDR.0b013e3181256492

PubMed Abstract | CrossRef Full Text | Google Scholar

Cavicchi, M., Beau, P., Crenn, P., Degott, C., and Messing, B. (2000). Prevalence of liver disease and contributing factors in patients receiving home parenteral nutrition for permanent intestinal failure. Ann. Intern. Med 132 (7), 525–532. doi:10.7326/0003-4819-132-7-200004040-00003

PubMed Abstract | CrossRef Full Text | Google Scholar

Davies, R., Cock, D., Kearsley-Fleet, L., Southwood, T., Baildam, E., Beresford, M. W., et al. (2019). The risk of uveitis in patients with JIA receiving etanercept: the challenges of analysing real-world data. Rheumatology 14, 33. doi:10.1093/rheumatology/kez449

CrossRef Full Text | Google Scholar

El Kasmi, K. C., Anderson, A. L., Devereaux, M. W., Fillon, S. A., Harris, J. K., Lovell, M. A., et al. (2012). Toll-like receptor 4-dependent Kupffer cell activation and liver injury in a novel mouse model of parenteral nutrition and intestinal injury. Hepatology 55 (5), 1518–1528. doi:10.1002/hep.25500

PubMed Abstract | CrossRef Full Text | Google Scholar

El Kasmi, K. C., Anderson, A. L., Devereaux, M. W., Vue, P. M., Zhang, W., Setchell, K. D., et al. (2013). Phytosterols promote liver injury and Kupffer cell activation in parenteral nutrition-associated liver disease. Sci. Transl. Med 5 (206), 137. doi:10.1126/scitranslmed.3006898

CrossRef Full Text | Google Scholar

Fang, R. F., Li, J. H., Zhang, J., Gao, S. Q., Pan, Y. T., and Zhan, S. Y. (2016). Association between statins use and liver injury based on prescription sequence symmetry analysis. Zhonghua Liuxingbingxue Zazhi 37 (7), 935–939. doi:10.3760/cma.j.issn.0254-6450.2016.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Gabe, S. M. (2013). Lipids and liver dysfunction in patients receiving parenteral nutrition. Curr. Opin. Clin. Nutr. Metab. Care 16 (2), 150–155. doi:10.1097/MCO.0b013e32835dbcf4

PubMed Abstract | CrossRef Full Text | Google Scholar

Hallas, J. (1996). Evidence of depression provoked by cardiovascular medication: a prescription sequence symmetry analysis. Epidemiology 7 (5), 478–484.

PubMed Abstract | CrossRef Full Text | Google Scholar

Hallas, J., Gaist, D., and Bjerrum, L. (1997). The waiting time distribution as a graphical approach to epidemiologic measures of drug utilization. Epidemiology 8 (6), 666–670. doi:10.1097/00001648-199710000-00009

PubMed Abstract | CrossRef Full Text | Google Scholar

Hojsak, I., and Kolacek, S. (2014). Fat overload syndrome after the rapid infusion of SMOFlipid emulsion. JPEN - J. Parenter. Enter. Nutr 38 (1), 119–121. doi:10.1177/0148607113482001

CrossRef Full Text | Google Scholar

Huard, G., Fiel, M. I., Moon, J., Iyer, K., and Schiano, T. D. (2018). Prevalence, evolution, and risk factors for advanced liver fibrosis in adults undergoing intestinal transplantation. JPEN - J. Parenter. Enter. Nutr 42 (7), 1195–1202. doi:10.1002/jpen.1148

CrossRef Full Text | Google Scholar

Izyan, A. W., Pratt, N. L., Kalisch, L. M., and Roughead, E. E. (2014). Comparing time to adverse drug reaction signals in a spontaneous reporting database and a claims database: a case study of rofecoxib-induced myocardial infarction and rosiglitazone-induced heart failure signals in Australia. Drug Saf 37 (1), 53–64. doi:10.1007/s40264-013-0124-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Javid, P. J., Greene, A. K., Garza, J., Gura, K., Alwayn, I. P., Voss, S., et al. (2005). The route of lipid administration affects parenteral nutrition-induced hepatic steatosis in a mouse model. J. Pediatr. Surg 40 (9), 1446–1453. doi:10.1016/j.jpedsurg.2005.05.045

PubMed Abstract | CrossRef Full Text | Google Scholar

Kapoor, V., Malviya, M. N., and Soll, R. (2019a). Lipid emulsions for parenterally fed preterm infants. Cochrane Database Syst. Rev 6, CD013163. doi:10.1002/14651858.CD013163.pub2

PubMed Abstract | CrossRef Full Text | Google Scholar

Kapoor, V., Malviya, M. N., and Soll, R. (2019b). Lipid emulsions for parenterally fed term and late preterm infants. Cochrane Database Syst. Rev 6, CD013171. doi:10.1002/14651858.CD013171.pub2

PubMed Abstract | CrossRef Full Text | Google Scholar

Kubota, K. (2016). Symmetry analysis for monitoring safety of newly marketed drugs. Pharmacoepidemiol. Drug Saf 25 (3), 349–351. doi:10.1002/pds.3886

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, J. A., and Teckman, J. H. (2015). Controversies in the mechanism of total parenteral nutrition induced pathology. Children 2 (3), 358–370. doi:10.3390/children2030358

PubMed Abstract | CrossRef Full Text | Google Scholar

Lai, E. C., Hsieh, C. Y., Kao Yang, Y. H., and Lin, S. J. (2014). Detecting potential adverse reactions of sulpiride in schizophrenic patients by prescription sequence symmetry analysis. PLoS One 9 (2), e89795. doi:10.1371/journal.pone.0089795

PubMed Abstract | CrossRef Full Text | Google Scholar

Lai, E. C., Pratt, N., Hsieh, C. Y., Lin, S. J., Pottegard, A., Roughead, E. E., et al. (2017). Sequence symmetry analysis in pharmacovigilance and pharmacoepidemiologic studies. Eur. J. Epidemiol 32 (7), 567–582. doi:10.1007/s10654-017-0281-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Lal, S., Pironi, L., Wanten, G., Arends, J., Bozzetti, F., Cuerda, C., et al. (2018). Clinical approach to the management of intestinal failure associated liver disease (IFALD) in adults: a position paper from the home artificial nutrition and chronic intestinal failure special interest group of ESPEN. Clin. Nutr 37 (6 Pt A), 1794–1797. doi:10.1016/j.clnu.2018.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Le, H. D., Meisel, J. A., de Meijer, V. E., Fallon, E. M., Gura, K. M., Nose, V., et al. (2012). Docosahexaenoic acid and arachidonic acid prevent essential fatty acid deficiency and hepatic steatosis. JPEN - J. Parenter. Enter. Nutr 36 (4), 431–441. doi:10.1177/0148607111414580

CrossRef Full Text | Google Scholar

Levit, O. L., Calkins, K. L., Gibson, L. C., Kelley-Quon, L., Robinson, D. T., Elashoff, D. A., et al. (2016). Low-dose intravenous soybean oil emulsion for prevention of cholestasis in preterm Neonates. JPEN - J. Parenter. Enter. Nutr 40 (3), 374–382. doi:10.1177/0148607114540005

CrossRef Full Text | Google Scholar

Li, L., Jiang, W., and Wang, J. (2007). Clinical analysis of 275 cases of acute drug-induced liver disease. Front. Med. China 1 (1), 58–61. doi:10.1007/s11684-007-0012-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Maclure, M., Fireman, B., Nelson, J. C., Hua, W., Shoaibi, A., Paredes, A., et al. (2012). When should case-only designs be used for safety monitoring of medical products?. Pharmacoepidemiol. Drug Saf 21 (Suppl 1), 50–61. doi:10.1002/pds.2330

PubMed Abstract | CrossRef Full Text | Google Scholar

Matarese, L. E., O'Keefe, S. J., Kandil, H. M., Bond, G., Costa, G., and Abu-Elmagd, K. (2005). Short bowel syndrome: clinical guidelines for nutrition management. Nutr. Clin. Pract 20 (5), 493–502. doi:10.1177/0115426505020005493

PubMed Abstract | CrossRef Full Text | Google Scholar

Meyerson, C., and Naini, B. V. (2019). Something old, something new: liver injury associated with Total parenteral nutrition therapy and immune checkpoint inhibitors. Hum. Pathol 96, 39–47. doi:10.1016/j.humpath.2019.10.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Morris, J. A., and Gardner, M. J. (1988). Calculating confidence intervals for relative risks (odds ratios) and standardised ratios and rates. Br. Med. J 296 (6632), 1313–1316. doi:10.1136/bmj.296.6632.1313

CrossRef Full Text | Google Scholar

Newcombe, R. G. (1998). Two-sided confidence intervals for the single proportion: comparison of seven methods. Stat. Med 17 (8), 857–872. doi:10.1002/(sici)1097-0258(19980430)17:8<857::aid-sim777>;2-e

PubMed Abstract | CrossRef Full Text | Google Scholar

Palova, S., Charvat, J., and Kvapil, M. (2008). Comparison of soybean oil- and olive oil-based lipid emulsions on hepatobiliary function and serum triacylglycerols level during realimentation. J. Int. Med. Res 36 (3), 587–593. doi:10.1177/147323000803600326

PubMed Abstract | CrossRef Full Text | Google Scholar

Peden, V. H., Witzleben, C. L., and Skelton, M. A. (1971). Total parenteral nutrition. J. Pediatr 78 (1), 180–181. doi:10.1016/s0022-3476(71)80289-5

CrossRef Full Text | Google Scholar

Pironi, L., Arends, J., Bozzetti, F., Cuerda, C., Gillanders, L., Jeppesen, P. B., et al. (2016). ESPEN guidelines on chronic intestinal failure in adults. Clin. Nutr 35 (2), 247–307. doi:10.1016/j.clnu.2016.01.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Pratt, N., Chan, E. W., Choi, N. K., Kimura, M., Kimura, T., Kubota, K., et al. (2015). Prescription sequence symmetry analysis: assessing risk, temporality, and consistency for adverse drug reactions across datasets in five countries. Pharmacoepidemiol. Drug Saf 24 (8), 858–864. doi:10.1002/pds.3780

PubMed Abstract | CrossRef Full Text | Google Scholar

Pratt, N. L., Ilomaki, J., Raymond, C., and Roughead, E. E. (2014). The performance of sequence symmetry analysis as a tool for post-market surveillance of newly marketed medicines: a simulation study. BMC Med. Res. Methodol 14, 66. doi:10.1186/1471-2288-14-66

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, Y., Shu, Z. J., Wang, H., Barnes, J. L., Yeh, C. K., Ghosh, P. M., et al. (2018). Altered expression of hepatic beta-adrenergic receptors in aging rats: implications for age-related metabolic dysfunction in liver. Am. J. Physiol. Regul. Integr. Comp. Physiol 314 (4), R574–R583. doi:10.1152/ajpregu.00372.2017

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsiropoulos, I., Andersen, M., and Hallas, J. (2009). Adverse events with use of antiepileptic drugs: a prescription and event symmetry analysis. Pharmacoepidemiol. Drug Saf 18 (6), 483–491. doi:10.1002/pds.1736

PubMed Abstract | CrossRef Full Text | Google Scholar

Tyree, P. T., Lind, B. K., and Lafferty, W. E. (2006). Challenges of using medical insurance claims data for utilization analysis. Am. J. Med. Qual 21 (4), 269–275. doi:10.1177/1062860606288774

PubMed Abstract | CrossRef Full Text | Google Scholar

Wahab, I. A., Pratt, N. L., Wiese, M. D., Kalisch, L. M., and Roughead, E. E. (2013). The validity of sequence symmetry analysis (SSA) for adverse drug reaction signal detection. Pharmacoepidemiol. Drug Saf 22 (5), 496–502. doi:10.1002/pds.3417

PubMed Abstract | CrossRef Full Text | Google Scholar

Wanten, G. J., and Calder, P. C. (2007). Immune modulation by parenteral lipid emulsions. Am. J. Clin. Nutr 85 (5), 1171–1184. doi:10.1093/ajcn/85.5.1171

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, L., Li, J. H., Zhao, K., and Wu, H. Y. (2015). Incidence and in-hospital mortality of acute aortic dissection in China: analysis of China Health Insurance Research (CHIRA) Data 2011. J. Geriatr. Cardiol 12 (5), 502–506. doi:10.11909/j.issn.1671-5411.2015.05.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., Zhou, X., Gao, S., Lin, H., Xie, Y., Feng, Y., et al. (2018). Evaluation of electronic healthcare databases for post-marketing drug safety surveillance and pharmacoepidemiology in China. Drug Saf 41 (1), 125–137. doi:10.1007/s40264-017-0589-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Yong, M., Xianjun, X., Jinghu, L., and Yunyun, F. (2018). Effect of health insurance on direct hospitalisation costs for in-patients with ischaemic stroke in China. Aust. Health Rev 42 (1), 39–44. doi:10.1071/ah16056

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, Y. C., Mao, Y. M., Chen, C. W., Chen, J. J., Chen, J., Cong, W. M., et al. (2017). CSH guidelines for the diagnosis and treatment of drug-induced liver injury. Hepatol. Int 11 (3), 221–241. doi:10.1007/s12072-017-9793-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Wu, J., Zhang, Z., Qian, H., Wang, Y., Yang, M., et al. (2019). Association of atorvastatin with the risk of hepatotoxicity: a pilot prescription sequence symmetry analysis. Therapeut. Clin. Risk Manag 15, 803–810. doi:10.2147/TCRM.S204860

CrossRef Full Text | Google Scholar

Keywords: lipid emulsion, hepatic protector, prescription sequence symmetry analysis, pharmacoepidemiology, health insurance database, drug safety

Citation: Li X-x, Cheng Y-c, Zhai S-d, Yao P, Zhan S-y and Shi L-w (2021) Risk of Liver Injury Associated with Intravenous Lipid Emulsions: A Prescription Sequence Symmetry Analysis. Front. Pharmacol. 12:589091. doi: 10.3389/fphar.2021.589091

Received: 25 August 2020; Accepted: 12 January 2021;
Published: 01 March 2021.

Edited by:

David Sacerdoti, University of Verona, Italy

Reviewed by:

Jia-bo Wang, Fifth Medical Center of the PLA General Hospital, China
Shaowen Tang, Nanjing Medical University, China

Copyright © 2021 Li, Cheng, Zhai, Yao, Zhan and Shi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Si-yan Zhan,; Lu-wen Shi,

These authors share first authorship