Hepatotoxic pyrrolizidine alkaloids induce DNA damage response in rat liver in a 28-day feeding study

Modeling HepaRG metabolome responses to pyrrolizidine alkaloid exposure for insight into points of departure and modes of action

The new challenges in toxicology demand novel and innovative in vitro approaches for deriving points of departure (PODs) and determining the mode of action (MOA) of chemicals. Therefore, the aim of this original study was to couple in vitro studies with untargeted metabolomics to model the concentration-response of extra- and intracellular metabolome data on human HepaRG cells treated for 48 h with three pyrrolizidine alkaloids (PAs): heliotrine, retrorsine and lasiocarpine. Modeling revealed that the three PAs induced various monotonic and, importantly, biphasic curves of metabolite content. Based on unannotated metabolites, the endometabolome was more sensitive than the exometabolome in terms of metabolomic effects, and benchmark concentrations (BMCs) confirmed that lasiocarpine was the most hepatotoxic PA. Regarding its MOA, impairment of lipid metabolism was highlighted at a very low BMC (first quartile, 0.003 µM). Moreover, results confirmed that lasiocarpine targets bile acids, as well as amino acid and steroid metabolisms. Analysis of the endometabolome, based on coupling concentration-response and PODs, gave encouraging results for ranking toxins according to their hepatotoxic effects. Therefore, this novel approach is a promising tool for next-generation risk assessment, readily applicable to a broad range of compounds and toxic endpoints.

Mass spectrometric analysis strategies for pyrrolizidine alkaloids

Pyrrolizidine alkaloids (PAs) in food and natural preparations have received widespread attention due to their hepatotoxicity, genotoxicity, and embryotoxicity. Mass spectrometry (MS), as a high resolution, high sensitive, and high throughput detection tool, has been the most commonly used technique for the determination of PAs. The continuous advancement of new technologies, methods, and strategies in the field of MS has contributed to the improvement of the analytical efficiency and methodological enhancement of PAs. This paper provides an overview of the structure, toxicity properties and commonly employed analytical methods, focusing on the concepts, advances, and novel techniques and applications of MS-based methods for the analysis of PAs. Additionally, the remaining challenges, future perspectives, and trends for PA detection are discussed. This review provides a reference for toxicological studies of PAs, content monitoring, and the establishment of quality control and safety standards for herbal and food products.

DNA modifications: Biomarkers for the exposome?

The concept of the exposome is the encompassing of all the environmental exposures, both exogenous and endogenous, across the life course. Many, if not all, of these exposures can result in the generation of reactive species, and/or the modulation of cellular processes, that can lead to a breadth of modifications of DNA, the nature of which may be used to infer their origin. Because of their role in cell function, such modifications have been associated with various major human diseases, including cancer, and so their assessment is crucial. Historically, most methods have been able to only measure one or a few DNA modifications at a time, limiting the information available. With the development of DNA adductomics, which aims to determine the totality of DNA modifications, a far more comprehensive picture of the DNA adduct burden can be gained. Importantly, DNA adductomics can facilitate a "top-down" investigative approach whereby patterns of adducts may be used to trace and identify the originating exposure source. This, together with other 'omic approaches, represents a major tool for unraveling the complexities of the exposome and hence allow a better a understanding of the environmental origins of disease.

OCT1-dependent uptake of structurally diverse pyrrolizidine alkaloids in human liver cells is crucial for their genotoxic and cytotoxic effects

Pyrrolizidine alkaloids (PAs) are important plant hepatotoxins, which occur as contaminants in plant-based foods, feeds and phytomedicines. Numerous studies demonstrated that the genotoxicity and cytotoxicity of PAs depend on their chemical structure, allowing for potency ranking and grouping. Organic cation transporter-1 (OCT1) was previously shown to be involved in the cellular uptake of the cyclic PA diesters monocrotaline, retrorsine and senescionine. However, little is known about the structure-dependent transport of PAs. Therefore, we investigated the impact of OCT1 on the uptake and toxicity of three structurally diverse PAs (heliotrine, lasiocarpine and riddelliine) differing in their degree and type of esterification in metabolically competent human liver cell models and hamster fibroblasts. Human HepG2-CYP3A4 liver cells were exposed to the respective PA in the presence or absence of the OCT1-inhibitors D-THP and quinidine, revealing a strongly attenuated cytotoxicity upon OCT1 inhibition. The same experiments were repeated in V79-CYP3A4 hamster fibroblasts, confirming that OCT1 inhibition prevents the cytotoxic effects of all tested PAs. Interestingly, OCT1 protein levels were much lower in V79-CYP3A4 than in HepG2-CYP3A4 cells, which correlated with their lower susceptibility to PA-induced cytotoxicity. The cytoprotective effect of OCT1 inhibiton was also demonstrated in primary human hepatocytes following PA exposure. Our experiments further showed that the genotoxic effects triggered by the three PAs are blocked by OCT1 inhibition as evidenced by strongly reduced γH2AX and p53 levels. Consistently, inhibition of OCT1-mediated uptake suppressed the activation of the DNA damage response (DDR) as revealed by decreased phosphorylation of checkpoint kinases upon PA treatment. In conclusion, we demonstrated that PAs, independent of their degree of esterification, are substrates for OCT1-mediated uptake into human liver cells. We further provided evidence that OCT1 inhibition prevents PA-triggered genotoxicity, DDR activation and subsequent cytotoxicity. These findings highlight the crucial role of OCT1 together with CYP3A4-dependent metabolic activation for PA toxicity.

Hepatotoxicity screening and ranking of structurally different pyrrolizidine alkaloids in zebrafish

Pyrrolizidine alkaloids (PAs) are phytotoxins distributed in ∼6000 plant species. PA-contaminated/containing foodstuffs/herbs/supplements pose a potential threat to human health. Various regulatory authorities established different PA margins of exposure assuming an equal hepatotoxic potency of structurally diverse PAs, although they exhibit different toxic potencies. Therefore, understanding hepatotoxic potencies of different PAs would facilitate a more appropriate risk assessment of PA exposure. In this study, a zebrafish model, which mimics physiological processes of absorption, distribution, metabolism, and excretion, was selected to evaluate acute hepatotoxic potency of different PAs (7 PAs and 2 PA N-oxides) and explore possible physiological pathways involved in PA-induced hepatotoxicity. After 6 h oral administration, PAs caused distinct structure-dependent hepatotoxicity with a series of biochemical and histological changes in zebrafish. Based on the measured toxicological endpoints, the relative toxic potency order of different PAs was derived as lasiocarpine ∼ retrorsine > monocrotaline > riddelliine > clivorine > heliotrine > retrorsine N-oxide ∼ riddelliine N-oxide≫>platyphyline. Further, compared to control group, different upregulation/downregulation of mRNA expression in PA-treated groups indicated that inflammation, apoptosis, and steatosis were involved in PA-induced hepatotoxicity in zebrafish. These findings demonstrate that zebrafish model is useful for screening and ranking hepatotoxicity of PAs with diverse structures, which would facilitate the more accurate risk assessment of PA exposure.

Chromatographic Methods for Detection and Quantification of Pyrrolizidine Alkaloids in Flora, Herbal Medicines, and Food: An Overview

Pyrrolizidine alkaloids (PAs) are natural toxins produced by some plants that gained special interest due to their potential hazardous effects in humans and animals. These substances have been found in wild flora, herbal medicines and food products raising health concerns. Recently, maximum concentration levels of PAs were established for some food products; however, maximum daily intake frequently surpasses the upper limit set by the competent authorities posing a health risk. Given the scarcity or absence of occurrence data on PAs in many products, there is an urgent need to measure their levels and establish safety intake levels. Analytical methods have been reported to detect and quantify PAs in different matrices. The commonly used chromatographic methodologies provides accurate and reliable results. Analytical methods include diverse steps as extraction and sample preparation procedures that are critical for sensitivity and selectivity of the analytical method. Great efforts have been directed toward optimization of extraction procedures, clean up and chromatographic conditions to improve recovery, reduce matrix effects, and achieve low limits of detection and quantification. Therefore, this paper aims to give a general overview about the occurrence of PAs in flora, herbal medicines, and foodstuff; and discuss the different chromatographic methodologies used for PAs analysis, namely extraction and sample preparation procedures and chromatographic conditions.

Comparative Pharmacokinetic Study of Two Pyrrolizidine Alkaloids Lasiocarpine and Heliotrine in Rats

Lasiocarpine (LAS) and heliotrine (HEL) are two different ester types of toxic pyrrolizidine alkaloids (PAs): open-chain diester and monoester. However, the pharmacokinetics of these two types of PAs in rats have not been reported. In the present study, two LC-MS/MS methods for determining LAS and HEL were established and validated. The methods exhibited good linearity, accuracy, and precision and were then applied to a comparative pharmacokinetic study. After intravenous administration to male rats at 1 mg/kg, the AUC_0-t values of LAS and HEL were 336 ± 26 ng/mL × h and 170 ± 5 ng/mL × h. After oral administration at 10 mg/kg, the AUC_0-t of LAS was much lower than that of HEL (18.2 ± 3.8 ng/mL × h vs. 396 ± 18 ng/mL × h), while the C_max of LAS was lower than that of HEL (51.7 ± 22.5 ng/mL × h vs. 320 ± 26 ng/mL × h). The absolute oral bioavailability of LAS was 0.5%, which was significantly lower than that of HEL (23.3%). The results revealed that the pharmacokinetic behaviors of LAS differed from that of HEL.

Potency ranking of pyrrolizidine alkaloids in metabolically competent human liver cancer cells and primary human hepatocytes using a genotoxicity test battery

Pyrrolizidine alkaloids (PAs) occur as contaminants in plant-based foods and herbal medicines. Following metabolic activation by cytochrome P450 (CYP) enzymes, PAs induce DNA damage, hepatotoxicity and can cause liver cancer in rodents. There is ample evidence that the chemical structure of PAs determines their toxicity. However, more quantitative genotoxicity data are required, particularly in primary human hepatocytes (PHH). Here, the genotoxicity of eleven structurally different PAs was investigated in human HepG2 liver cells with CYP3A4 overexpression and PHH using an in vitro test battery. Furthermore, the data were subject to benchmark dose (BMD) modeling to derive the genotoxic potency of individual PAs. The cytotoxicity was initially determined in HepG2-CYP3A4 cells, revealing a clear structure-toxicity relationship for the PAs. Importantly, experiments in PHH confirmed the structure-dependent toxicity and cytotoxic potency ranking of the tested PAs. The genotoxicity markers γH2AX and p53 as well as the alkaline Comet assay consistently demonstrated a structure-dependent genotoxicity of PAs in HepG2-CYP3A4 cells, correlating well with their cytotoxic potency. BMD modeling yielded BMD values in the range of 0.1-10 µM for most cyclic and open diesters, followed by the monoesters. While retrorsine showed the highest genotoxic potency, monocrotaline and lycopsamine displayed the lowest genotoxicity. Finally, experiments in PHH corroborated the genotoxic potency ranking, and revealed genotoxic effects even in the absence of detectable cytotoxicity. In conclusion, our findings strongly support the concept of grouping PAs into potency classes and help to pave the way for a broader acceptance of relative potency factors in risk assessment.

Physiologically-Based Kinetic Modeling Predicts Similar In Vivo Relative Potency of Senecionine N-Oxide for Rat and Human at Realistic Low Exposure Levels

SCOPE

This study aims to determine if previously developed physiologically-based kinetic (PBK) model in rat can be modified for senecionine (SEN) and its N-oxide (SENO), and be used to investigate potential species differences between rat and human in relative potency (REP) of the N-oxide relative to the parent pyrrolizidine alkaloid (PA).

METHODS AND RESULTS

In vitro derived kinetic parameters including the apparent maximum velocities (V_max ) and Michaelis-Menten constants (K_m ) for SENO reduction and SEN clearance are used to define the PBK models. The rat model is validated with published animal data, and the toxicokinetic profiles of SEN from either orally-administered SENO or SEN are simulated. REP values of SENO relative to SEN amount to 0.84 and 0.89 in rat and human, respectively.

CONCLUSION

The REP value can be dose- and species-dependent, with the values for rat and human being comparable at low realistic exposure scenarios. In summary, PBK modeling serves as a valuable New Approach Methodology (NAM) tool for predicting REP values of PA-N-oxides and may actually result in more accurate REP values for human risk assessment than what would be defined using in vivo animal experiments.

Current Trends in Toxicity Assessment of Herbal Medicines: A Narrative Review

Even in modern times, the popularity level of medicinal plants and herbal medicines in therapy is still high. The World Health Organization estimates that 80% of the population in developing countries uses these types of remedies. Even though herbal medicine products are usually perceived as low risk, their potential health risks should be carefully assessed. Several factors can cause the toxicity of herbal medicine products: plant components or metabolites with a toxic potential, adulteration, environmental pollutants (heavy metals, pesticides), or contamination of microorganisms (toxigenic fungi). Their correct evaluation is essential for the patient’s safety. The toxicity assessment of herbal medicine combines in vitro and in vivo methods, but in the past decades, several new techniques emerged besides conventional methods. The use of omics has become a valuable research tool for prediction and toxicity evaluation, while DNA sequencing can be used successfully to detect contaminants and adulteration. The use of invertebrate models (Danio renio or Galleria mellonella) became popular due to the ethical issues associated with vertebrate models. The aim of the present article is to provide an overview of the current trends and methods used to investigate the toxic potential of herbal medicinal products and the challenges in this research field.

Transcriptome of pituitary function changes in rat model of high altitude cerebral edema

High altitude cerebral edema (HACE) is a serious subtype of acute mountain sickness (AMS). Studies have suggested that increased expression of corticotropin releasing hormone receptor 1 (CRFR1) in pituitary is related to the development of HACE, but no study has revealed the molecular landscape of pituitary function changes in this process. Rat model of HACE was established by simulating the high-altitude hypobaric hypoxia environment. Then RNA-sequencing was performed of rat pituitary gland (PG) in HACE and non-HACE groups. The function annotations, enrichment analysis, protein-protein interaction (PPI) network, chromosome location and drug repositioning of differentially expressed genes (DEGs) were explored based on the transcriptomic data. And we found pituitary secretion function was disordered in HACE, which was partly due to activated inflammation and oxidative stress. In addition, we identified potential biomarkers for early recognition of pituitary dysfunction and potential protective drugs for pituitary function in HACE.

Hematotoxicity and nephrotoxicity of long-term administration of Guiera senegalensis (J.F. Gme), Cassia occidentalis (Linn), and Ziziphus mauritiana (Lam) leaves obtained in Birnin Kebbi, Nigeria

Introduction: Previously, we established the phytochemical composition and short-term administration safety of Guiera senegalensis (sabara), Cassia occidentalis (coffee senna), and Ziziphus mauritiana (jujube) leaves, which are common medicinal plants in Northern Nigeria. In the current study, heavy metal contents and long-term administration effects of the plants’ leaf extracts on hematological parameters and the kidneys of albino rats (Rattus norvegicus) were investigated. The heavy metals analyzed were copper, lead, cadmium, nickel, and manganese, while the hematological parameters evaluated were packed cell volume, hemoglobin, red blood cells, white blood cells, lymphocytes, and monocytes. Methods: Twenty-four mixed-sex rats were distributed into four groups of six rats each. Group 1 was made the control, while groups 2, 3, and 4 were administered 1000 mg kg-1 one of the plants extracts for 90 days. Blood and kidney samples were collected across the groups for hematological and histopathological examinations. Results: The heavy metals were present in the extracts within the World Health Organization’s acceptable limits. The treated rats were anemic compared to the control. However, on average, only the C. occidentalis group showed significant differences (P<0.05) in hematological parameters. Unlike the control, the kidneys of the rats fed with Z. mauritiana and G. senegalensis showed vacuolation of cytoplasm and tubular degeneration, while the C. occidentalis-fed rats had inflammation and dilated Bowman’s capsules. Conclusion: These findings reveal that constant administration of high doses of the extracts for a long time may cause health hazards. People are advised to seek an expert’s advice before using the plants.

The chemical structure impairs the intensity of genotoxic effects promoted by 1,2-unsaturated pyrrolizidine alkaloids in vitro

1,2-unsaturated pyrrolizidine alkaloids (PAs) represent a large group of secondary plant metabolites exhibiting hepatotoxic, genotoxic, and carcinogenic properties upon bioactivation. To examine how the degree of esterification affects the genotoxic profile of PA we investigated cytotoxicity, histone H2AX phosphorylation, DNA strand break induction, cell cycle perturbation, micronuclei formation, and aneugenic effects in different cell models. Analysis of cytotoxicity and phosphorylation of histone H2AX was structure- and concentration-dependent: diester-type PAs (except monocrotaline) showed more pronounced effects than monoester-type PAs. Cell cycle analysis identified that diester-type PAs induced a S-phase arrest and a decrease in the occurrence of cells in the G1-phase. The same structure-dependency was observed by flow-cytometric analysis of PA-induced micronuclei in CYP3A4-overexpressing V79 cells. Analysis of centromeres induced by lasiocarpine in the micronuclei by fluorescence in situ hybridization indicated an aneugenic effect in V79_h3A4 cells. Comet assays revealed no significant induction of DNA strand breaks for all investigated PAs. Overall, diester-type PAs induced more pronounced effects than monoester-type PAs. Furthermore, our results indicate aneugenic effects upon exposure towards lasiocarpine in vitro. These data improve our understanding how structural features of PA influence the genotoxic profile. Especially, the monoester-type PAs seem to induce less severe effects than other PAs.

Plasma concentration profiles for hepatotoxic pyrrolizidine alkaloid senkirkine in humans extrapolated from rat data sets using a simplified physiologically based pharmacokinetic model

AIM

The main aim of the current study was to obtain forward dosimetry assessments of pyrrolizidine alkaloid senkirkine plasma and liver concentrations by setting up a human physiologically based pharmacokinetic (PBPK) model based on the limited information available.

BACKGROUND

The risks associated with plant-derived pyrrolizidine alkaloids as natural toxins have been assessed.

OBJECTIVE

The pyrrolizidine alkaloid senkirkine was investigated because it was analyzed in a European transcriptomics study of natural hepatotoxins and in a study of the alkaloidal constituents of traditional Japanese food plants Petasites japonicus. The in silico human plasma and liver concentrations of senkirkine were modeled using doses reported for acute-term toxicity in humans.

METHODS

Using a simplified PBPK model established using rat pharmacokinetic data, forward dosimetry was conducted. Since in vitro rat and human intrinsic hepatic clearances were similar; an allometric scaling approach was applied to rat parameters to create a human PBPK model.

RESULTS

After oral administration of 1.0 mg/kg in rats in vivo, water-soluble senkirkine was absorbed and cleared from plasma to two orders of magnitude below the maximum concentration in 8 h. Human in silico senkirkine plasma concentration curves were generated after virtual daily oral administrations of 3.0 mg/kg senkirkine (the dose involved in an acute fatal hepatotoxicity case). A high concentration of senkirkine in the culture medium caused in vitro hepatotoxicity as evidenced by lactate dehydrogenase leakage from human hepatocyte-like HepaRG cells.

CONCLUSION

Higher virtual concentrations of senkirkine in human liver and plasma than those in rat plasma were estimated using the current rat and human PBPK models. Current simulations suggest that if P. japonicus (a water-soluble pyrrolizidine alkaloid-producing plant) is ingested daily as food, hepatotoxic senkirkine could be continuously present in human plasma and liver.

Fasting augments pyrrolizidine alkaloid-induced hepatotoxicity

Pyrrolizidine alkaloids (PAs) are phytotoxins widely present in various natural products and foodstuffs. The present study aims to investigate the effects of fasting on PA-induced hepatotoxicity and the underlying biochemical mechanisms. The results of hepatotoxic study showed that 15-h overnight fasting significantly exacerbated the hepatotoxicity of retrorsine (RTS, a representative toxic PA) in fasted rats compared to fed rats, as indicated by remarkably elevated plasma ALT and bilirubin levels and obvious liver histological changes. Further toxicokinetic studies revealed that fasting significantly enhanced cytochromes P450 enzymes (CYPs)-mediated metabolic activation of RTS leading to increased formation of pyrrole-protein adducts and thus decreased the in vivo exposure and excretion of both parent RTS and its N-oxide metabolite. Metabolic studies demonstrated that fasting induced enzyme activities of CYP1A2, CYP2B6 and CYP2E1 that participated in catalyzing RTS to its reactive pyrrolic metabolites. Moreover, fasting also dramatically decreased hepatic glutathione (GSH) content, which restricted the detoxification of GSH by neutralizing the reactive pyrrolic metabolite of RTS, further contributing to the enhanced hepatotoxicity. The present findings may have an impact on future PA toxicity tests with different dietary styles and/or risk assessment of metabolite-mediated toxins by considering the profound effects of fasting.

Pyrrolizidine alkaloid-induced transcriptomic changes in rat lungs in a 28-day subacute feeding study

Pyrrolizidine alkaloids (PAs) are secondary plant metabolites synthesized by a wide range of plants as protection against herbivores. These toxins are found worldwide and pose a threat to human health. PAs induce acute effects like hepatic sinusoidal obstruction syndrome and pulmonary arterial hypertension. Moreover, chronic exposure to low doses can induce cancer and liver cirrhosis in laboratory animals. The mechanisms causing hepatotoxicity have been investigated previously. However, toxic effects in the lung are less well understood, and especially data on the correlation effects with individual chemical structures of different PAs are lacking. The present study focuses on the identification of gene expression changes in vivo in rat lungs after exposure to six structurally different PAs (echimidine, heliotrine, lasiocarpine, senecionine, senkirkine, and platyphylline). Rats were treated by gavage with daily doses of 3.3 mg PA/kg bodyweight for 28 days and transcriptional changes in the lung and kidney were investigated by whole-genome microarray analysis. The results were compared with recently published data on gene regulation in the liver. Using bioinformatics data mining, we identified inflammatory responses as a predominant feature in rat lungs. By comparison, in liver, early molecular consequences to PAs were characterized by alterations in cell-cycle regulation and DNA damage response. Our results provide, for the first time, information about early molecular effects in lung tissue after subacute exposure to PAs, and demonstrates tissue-specificity of PA-induced molecular effects.

The Food Contaminants Pyrrolizidine Alkaloids Disturb Bile Acid Homeostasis Structure-Dependently in the Human Hepatoma Cell Line HepaRG

Pyrrolizidine alkaloids (PAs) are a group of secondary plant metabolites being contained in various plant species. The consumption of contaminated food can lead to acute intoxications in humans and exert severe hepatotoxicity. The development of jaundice and elevated bile acid concentrations in blood have been reported in acute human PA intoxication, indicating a connection between PA exposure and the induction of cholestasis. Additionally, it is considered that differences in toxicity of individual PAs is based on their individual chemical structures. Therefore, we aimed to elucidate the structure-dependent disturbance of bile acid homeostasis by PAs in the human hepatoma cell line HepaRG. A set of 14 different PAs, including representatives of all major structural characteristics, namely, the four different necine bases retronecine, heliotridine, otonecine and platynecine and different grades of esterification, was analyzed in regard to the expression of genes involved in bile acid synthesis, metabolism and transport. Additionally, intra- and extracellular bile acid levels were analyzed after PA treatment. In summary, our data show significant structure-dependent effects of PAs on bile acid homeostasis. Especially PAs of diester type caused the strongest dysregulation of expression of genes associated with cholestasis and led to a strong decrease of intra- and extracellular bile acid concentrations.

Metabolism-mediated cytotoxicity and genotoxicity of pyrrolizidine alkaloids

Pyrrolizidine alkaloids (PAs) and PA N-oxides are common phytotoxins produced by over 6000 plant species. Humans are frequently exposed to PAs via ingestion of PA-containing herbal products or PA-contaminated foods. PAs require metabolic activation to form pyrrole-protein adducts and pyrrole-DNA adducts which lead to cytotoxicity and genotoxicity. Individual PAs differ in their metabolic activation patterns, which may cause significant difference in toxic potency of different PAs. This review discusses the current knowledge and recent advances of metabolic pathways of different PAs, especially the metabolic activation and metabolism-mediated cytotoxicity and genotoxicity, and the risk evaluation methods of PA exposure. In addition, this review provides perspectives of precision toxicity assessment strategies and biomarker development for the risk control and translational investigations of human intoxication by PAs.

Transcriptomics analysis of hepatotoxicity induced by the pesticides imazalil, thiacloprid and clothianidin alone or in binary mixtures in a 28-day study in female Wistar rats

Co-occurrence of pesticide residues in food commodities raises a potential safety issue as their mixture effects on human health are largely unknown. In a previous study, we reported the toxicological effects (pathology and histopathology) of imazalil (IMZ), thiacloprid (THI), and clothianidin (CTD) alone and in binary mixtures in a 28-day oral gavage study in female Wistar rats. Five dose levels (up to 350 mg/kg body weight/day) ranging from a typical toxicological reference value to a clear effect dose were applied. In the present study, we undertook a transcriptomics analysis of rat livers by means of total RNA sequencing (RNA-Seq). Bioinformatic data analysis involving Ingenuity Pathway Analysis (IPA) was used to gain mechanistic information on hepatotoxicity-related pathways affected after treatment with the pesticides, alone and in mixtures. Our data show that 2986 genes were differentially regulated by CTD while IMZ and THI had effects on 194 and 225 genes, respectively. All three individual compounds shared a common subset of genes whose network is associated with xenobiotic metabolism and nuclear receptor activation. Similar networks were retrieved for the mixtures. Alterations in the expression of individual genes were in line with the assumption of dose addition. Our results bring new insight into the hepatotoxicity mechanisms of IMZ, THI, and CTD and their mixtures.

Toxic Prediction of Pyrrolizidine Alkaloids and Structure-Dependent Induction of Apoptosis in HepaRG Cells

Pyrrolizidine alkaloids (PAs) are common phytotoxins and could cause liver genotoxicity/carcinogenicity following metabolic activation. However, the toxicity of different structures remains unclear due to the wide variety of PAs. In this study, the absorption, distribution, metabolism, excretion, and toxicity (ADMET) of 40 PAs were analyzed, and their toxicity was predicted by Komputer Assisted Technology (TOPKAT) using Discovery Studio software. The in silico results showed that all PAs except retronecine had good intestinal absorption, and all PAs were predicted to have different toxicity ranges. To verify the predictive results, 4 PAs were selected to investigate cell injury and possible mechanisms of the differentiation in HepaRG cells, including retronecine type of twelve-membered cyclic diester (retrorsine), eleven-membered cyclic diester (monocrotaline), noncyclic diester (retronecine), and platynecine type (platyphylline). After 24 h exposure, retronecine-type PAs exhibited concentration-dependent cytotoxicity. The high-content screening assay showed that cell oxidative stress, mitochondrial damage, endoplasmic reticulum stress, and the concentration of calcium ions increased, and neutral lipid metabolism was changed notably in HepaRG cells. Induced apoptosis by PAs was indicated by cell cycle arrest in the G2/M phase, disrupting the mitochondrial membrane potential. Overall, our study revealed structure-dependent cytotoxicity and apoptosis after PA exposure, suggesting that the prediction results of in silico have certain reference values for compound toxicity. A 1,2-membered cyclic diester seems to be a more potent apoptosis inducer than other PAs.

Lung injury induced by pyrrolizidine alkaloids depends on metabolism by hepatic cytochrome P450s and blood transport of reactive metabolites

Pyrrolizidine alkaloids (PAs) are common phytotoxins with both hepatotoxicity and pneumotoxicity. Hepatic cytochrome P450 enzymes are known to bioactivate PAs into reactive metabolites, which can interact with proteins to form pyrrole-protein adducts and cause intrahepatic cytotoxicity. However, the metabolic and initiation biochemical mechanisms underlying PA-induced pneumotoxicity remain unclear. To investigate the in vivo metabolism basis for PA-induced lung injury, this study used mice with conditional deletion of the cytochrome P450 reductase (Cpr) gene and resultant tissue-selective ablation of microsomal P450 enzyme activities. After oral exposure to monocrotaline (MCT), a pneumotoxic PA widely used to establish animal lung injury models, liver-specific Cpr-null (LCN) mice, but not extrahepatic Cpr-low (xh-CL) mice, had significantly lower level of pyrrole-protein adducts in the serum, liver and lungs compared with wild-type (WT) mice. While MCT-exposed LCN mice had significantly higher blood concentration of intact MCT, compared to MCT-exposed WT or xh-CL mice. Consistent with the MCT in vivo bioactivation data, MCT-induced lung injury, represented by vasculature damage, in WT and xh-CL mice but not LCN mice. Furthermore, reactive metabolites of MCT were confirmed to exist in the blood efflux from the hepatic veins of MCT-exposed rats. Our results provide the first mode-of-action evidence that hepatic P450s are essential for the bioactivation of MCT, and blood circulating reactive metabolites of MCT to the lung causes pneumotoxicity. Collectively, this study presents the scientific basis for the application of MCT in animal lung injury models, and more importantly, warrants public awareness and further investigations of lung diseases associated with exposure to not only MCT but also different PAs.