- Open Access
Analysis of gene expression changes in relation to toxicity and tumorigenesis in the livers of Big Blue transgenic rats fed comfrey (Symphytum officinale)
BMC Bioinformatics volume 7, Article number: S16 (2006)
Comfrey is consumed by humans as a vegetable and a tea, and has been used as an herbal medicine for more than 2000 years. Comfrey, however, is hepatotoxic in livestock and humans and carcinogenic in experimental animals. Our previous study suggested that comfrey induces liver tumors by a genotoxic mechanism and that the pyrrolizidine alkaloids in the plant are responsible for mutation induction and tumor initiation in rat liver.
In this study, we identified comfrey-induced gene expression profile in the livers of rats. Groups of 6 male transgenic Big Blue rats were fed a basal diet and a diet containing 8% comfrey roots, a dose that resulted in liver tumors in a previous carcinogenicity bioassay. The animals were treated for 12 weeks and sacrificed one day after the final treatment. We used a rat microarray containing 26,857 genes to perform genome-wide gene expression studies. Dietary comfrey resulted in marked changes in liver gene expression, as well as in significant decreases in the body weight and increases in liver mutant frequency. When a two-fold cutoff value and a P-value less than 0.01 were selected, 2,726 genes were identified as differentially expressed in comfrey-fed rats compared to control animals. Among these genes, there were 1,617 genes associated by Ingenuity Pathway Analysis with particular functions, and the differentially expressed genes in comfrey-fed rat livers were involved in metabolism, injury of endothelial cells, and liver injury and abnormalities, including liver fibrosis and cancer development.
The gene expression profile provides us a better understanding of underlying mechanisms for comfrey-induced hepatic toxicity. Integration of gene expression changes with known pathological changes can be used to formulate a mechanistic scheme for comfrey-induced liver toxicity and tumorigenesis.
Comfrey belongs to the family Boraginaceae. Three plant species in the genus Symphytum contribute to the crop known as comfrey, Symphytum officinale L. (wild or common comfrey; the major comfrey species), S. asperum Lepechin (prickly or rough comfrey), and S. x uplandicum Nyman (quaker, Russian, or blue comfrey; a natural hybrid of S. officinale L. and S. asperum Lepechin). Symphytum officinale L. is a tall perennial with large hairy leaves and small purple flowers . Comfrey has been used as an herbal medicine for more than two thousand years for the treatment of broken bones, tendon damage, ulcerations in the gastrointestinal tract, and lung congestion, as well as for wound healing and/or reducing joint inflammation when it is applied externally .
In addition to essential nutrients, comfrey also contains pyrrolizidine alkaloids (PAs). PAs are constituents of over 6000 plants, and many of them are hepatotoxic and carcinogenic in humans and animals . In the liver, PAs are transformed to pyrroles by the mixed-function oxidases. Pyrroles exert their toxic effect by reacting with cellular macromolecules, including proteins and DNA . Therefore, comfrey's therapeutic use might increase the risk of liver toxicity. Many countries including Canada, Germany, and the UK, have restricted its availability. In 2001, the US Food and Drug Administration requested voluntary compliance for the removal of products containing comfrey .
The major hepatotoxic manifestation in humans ingesting comfrey is the hepatic veno-occlusive lesion (VOD) [6, 7], also called sinusoidal obstruction syndrome (SOS) . Comfrey is also carcinogenic in rats, which suggests the potential tumorigenic effects of the plant . Hepatocellular adenomas were induced in rats receiving a diet containing comfrey. Feeding rats comfrey leaves produced a dose-dependent reduction in survival and an increase in liver tumor incidence . Comfrey roots are much more toxic than the leaves.
In a previous study, we developed evidence indicating that the liver tumors induced by feeding rats 2% comfrey root were generated by a genotoxic mechanism and that the PAs in the plant were responsible for mutation induction and tumor initiation in rat liver . Considering the fact that rats tolerate diets containing up to 33% comfrey leaves and 8% comfrey roots for relatively long periods of time (at least 6 months) , in the present study, we evaluated the mutagenicity of 8% comfrey root. Using a toxicogenomic approach, we analyzed the changes in global gene expressions in the liver of rats following comfrey-treatment.
Growth curve of rats fed with 8% comfrey root
Male Big Blue transgenic rats were fed with 8% comfrey root for 12 weeks. The mean body weight of the comfrey-fed rats was less than that of the vehicle controls throughout the study (Figure 1). Rats fed with comfrey weighed 5%, 23%, and 35% less than the control rats after 1, 6, and 12 weeks of the study, respectively, and displayed little weight gain after 6 weeks of feeding with comfrey.
Mutant frequency (MF) in the liver cII gene of comfrey-fed rats
The results of cII MF analyses in the comfrey-fed and control rats are shown in Figure 2. DNA from each liver was packaged 2–4 times either to confirm the MF or to obtain a minimum of 2 × 105 plaque-forming units for mutant detection. The MF for rats fed with 8% comfrey was 139 ± 35 (SD) × 10-6, which was similar to the MF previously detected in 2% comfrey-fed rats  and significantly increased over the control group (30 ± 16 × 10-6, P < 0.001).
Mutation spectrum in the liver cII gene from comfrey-fed rats
Comfrey-induced mutations in the liver cII gene were evaluated by DNA sequence analysis of 106 mutants isolated from 6 rats fed a diet containing 8% comfrey root. Since mutations that were found more than once among the mutants isolated from a single animal were assumed to be siblings and to represent only one independent mutation, a total of 99 independent mutations were identified. Table 1 summarizes the types of cII mutations observed in the livers of rats fed 8% comfrey compared with mutation spectra of control and 2% comfrey-fed rats that we reported previously . The overall pattern of mutations in 8% comfrey-fed rats differed significantly from controls (P < 0.0001), but did not differ from 2% comfrey-fed rats. Among the independent mutations, about 85% from both the comfrey-treated and control rats were base pair substitutions. G:C → T:A transversion (41%) was the major type of mutation in the 8% comfrey-fed rats, whereas G:C → A:T transition was the predominant mutation in the controls. In addition, a 13% frequency of tandem base substitutions was observed among the mutations from the 8% comfrey-fed rats. The 8% comfrey-fed rats had more large deletions than control rats, although the percentages of frameshift mutations were the same.
Gene expression regulated by comfrey-treatment
In order to determine gene expression changes associated with comfrey exposure, we employed the Genome Survey Microarray, which contains 26,857 verified rat genes. We isolated total RNA from the livers of 6 control and 6 comfrey-fed rats. For the purpose of visualizing the data, the intensities of the whole rat gene data were analyzed by Principal Components Analysis within ArrayTrack (Figure 3). It demonstrates that one array of sample 6 from the control group appears to be quite different from the rest of the arrays, which was further identified as an outlier array by the Pearson's correlation coefficient of pair-wise log2 intensity correlation (data not shown). This outlier was excluded from further data analysis. A separation between control and comfrey-treated groups was clearly observed, suggesting that there was a clear comfrey-treatment effect on liver gene expression (Figure 3). To select significant genes, minimum requirements were established for both a two-fold change in the gene expression compared to the controls and a P-value less than 0.01 for the difference. A total of 2,726 genes satisfied the requirements; 1,235 genes were up-regulated and 1,491 genes were down-regulated in response to comfrey treatment (Figure 4). Among the regulated genes, 1,671 were identified by Ingenuity Pathway Analysis. In this study, we focused on genes involved in metabolism, injury of endothelial cells, and liver injury and abnormalities.
Genes associated with drug metabolizing genes
Since the PAs in comfrey require metabolic activation to exert their biological effects, we investigated the gene expression changes of drug metabolizing genes. We observed the up- or down-regulation of many cytochrome P450 genes (e.g., Cyp2c12, Cyp3a18, Cyp4a12, Cyp26), glutathione S-transferases (Gsta3, Gstm3, and Gstp1), ATP-binding cassette transporters (e.g., Abcb9 and Abcc3), and other metabolism-associated genes, including NAD(P)H oxidoreductase (Nqo1) and aldehyde dehydrogenase (Aldh1a1). Table 2 shows the significant changes in gene expression related to phase I, II, and III drug metabolizing genes.
Genes involved in injury of endothelial cells
Considering the involvement of sinusoidal endothelial cells and sometimes hepatic venular endothelial cells in VOD/SOS , we focused on genes expressed in endothelial cells. Table 3 summarizes the alteration of genes related to cell death, apoptosis, and cell growth of endothelial cells. For example, we observed the up-regulation of previously reported tumor necrosis factors (Tnfsf10, Tnfrsf6, and Tnfrsf12a), transforming growth factor (Tgfb1 and Tgfb2), and other genes (such as, Mmp2, PAI-1 or Serpine1, Plaur, Col4a2, and Edn1) [8, 11, 12]. Several novel genes were responsive to comfrey treatment, including activating transcription factor 3 (Atf3), B-cell cLL/lymphoma 2 (Bcl2), caspase 1 (Casp1), mitogen-activated protein kinase 9 (Mapk9), and secreted phosphoprotein 1 (Spp1).
Genes involved in liver injury and abnormalities
Endothelial damage can lead to hepatocellular injury, and fibrotic reaction in the sinusoids is characterized as the later stages of VOD/SOS. Chronic comfrey treatment also resulted in the induction of a number of genes involved in liver injury and abnormalities. Significantly changed genes were divided into subsets based on functionality, and categories included cell death, growth, fibrosis, morphology, and liver cancer development (Table 4). Sixteen genes were functionally categorized in liver cell death. Besides the tumor necrosis factors and transforming growth factor, comfrey treatment resulted in strong up-regulation of cyclin-dependent kinase inhibitor 1A (Cdkn1a) and heme oxyenase 1 (Hmox1), as well as down-regulation of epidermal growth factor (Egf) and CCAAT/enhancer binding protein beta (Cebpb). Of note, genes involved in cell growth were up-regulated along with the genes involved in liver cell death, such as cyclin E1 (Ccne1), hepatocyte growth factor (Hgf), and tissue inhibitor of metalloproteinase 1 (Timp1).
Liver fibrosis is the hallmark of all chronic liver diseases, irrespective of their origin. Ingenuity Pathway Analysis showed that sixteen significantly regulated genes appeared to have a role in liver fibrosis (Table 4). Most of them encode cytokines (Hgf, Igf1, Lif, Tgfb1, and Tnfrsf6), angiogenic molecules (Agt, and Agtr1a), chemokines and their receptors (Ccr2, Cxcl4, and Cxcl12), coagulation cascade genes (PAI-1/Serpine1 and F11), and cytoskeletal molecules (Krt2-8). Moreover, eight genes related to the development of liver cancer were among the most discriminating genes involved in cancer development, such as growth arrest and DNA damage inducible 45 alpha (Gadd45a), B-cell cLL/lymphoma 2 (Bcl2), N-myc downstream regulated 1 (Ndr1), and xeroderma pigmentosum complementation group C (Xpc).
Comfrey was one of the most popular herbal teas in the world, including the United States. Although its popularity has declined due to the understanding of its dangers, it is still available commercially in several forms. The regular use of comfrey is a potential health risk owing to the presence of PAs. Comfrey contains as many as nine PAs, including acetyl intermedine, acetyl lycopsamine, echimidine, intermedine, lasiocarpine, lycopsamine, symlandine, symphytine, and symviridine [13–15]. The PA content of comfrey is less than 1% and depends on the plant part . Higher PA concentrations occur in the roots of comfrey than in the leaves, and commercially available comfrey tablets containing high levels of PAs are likely to be derived from comfrey roots [16, 17]. PAs are the leading plant toxins associated with disease in humans and animals.
Several cases of VOD/SOS associated with comfrey ingestion have been reported in humans [6, 7, 18–20], as well as in experimental animals . Deleve and her colleagues  have developed a reproducible animal model of hepatic VOD, in which the rats are gavaged with a single dose of PA, monocrotaline. The model exhibits the characteristic clinical and histological features of hepatic VOD, including the earliest manifestations (progressive injury to the sinusoidal endothelial cells and central vein endothelium), early VOD (centrilobular coagulative necrosis and severe sinusoidal injury, hemorrhage, and central vein endothelial damage), and late VOD (fibrotic occulusion of the central veins) . Comfrey and the specific alkaloids in comfrey (e.g., symphytine and lasiocarpine) induce hepatoadenoma and hemangioendothelial sarcomas in rats [9, 23, 24]. The mechanisms by which toxicity and carcinogenicity are produced are still not fully understood. Gene expression profiling offers a powerful approach for identifying differentially expressed genes and identifying mechanisms.
Gene expression was markedly affected (Figure 3) in the livers of rats exposed to 8% comfrey root, a dose that resulted in significant decreases in body weight (Figure 1) and increases in liver MF (Figure 2). Out of 26,857 genes evaluated, the expression of 4,132 (15%) and 9,937 (37%) genes were altered more than 2-fold and 1.5-fold, respectively. At P-values of 0.01 and 0.05, 7,518 (28% of expressed genes) and 1,0341 (39%) genes, respectively, displayed a significant effect after comfrey-treatment compared to control group. In this study, differential gene expression was considered significant for genes showing at least a 2-fold up- or down-change, and a P < 0.01. In total, 2,726 genes (10%) satisfied the requirements and about half of them were down-regulated and half up-regulated in response to comfrey exposure (Figure 4). Such a large number of significantly altered genes may partly reflect the therapeutic effects of comfrey exerted through plant components other than PAs . In the present study, we concentrated on the analysis of genes involved in metabolism, injury of endothelial cells, and liver injury and abnormalities.
Liver is the major organ for biotransformation of xenobiotics and drugs. PAs are metabolically activated to toxic, alkylating pyrroles by mixed-function oxidases. The cytochrome P450 (Cyp) superfamily contains 57 genes and plays a critical role in the phase I metabolism of a variety of xenobiotics including drugs, carcinogens, steroids and eicosanoids . In the present study, comfrey exposure resulted in changes in the expression of 19 Cyp genes (Table 2). Among phase I, II, and III drug metabolizing genes, Cyp2c12, Cyp7a1, Cyp26, Gsta3, and Abcc3 were increased 6-21-fold. In contrast, Cyp2c, Cyp39a1, Gstm3, Abcc8, and others were reduced in the comfrey-treated livers. It is known that many herbal/dietary constituents form reactive intermediates capable of irreversibly inhibiting various Cyps (reviewed in ). The resultant metabolites lead to Cyp inactivation by chemical modification of the heme, the apoprotein, or both, as a result of covalent binding of modified heme to the apoprotein. Phase II consists of conjugating enzymes, such as glutathione S-transferases (GSTs), UDP-glucuronosyltransferases (UGTs), and sulfatases. GSTs also participate in oxidative stress release pathways. The altered expression of phase I and phase II enzymes along with altered drug transport proteins (phase III) could contribute to the increased susceptibility of rats to carcinogenic chemicals, such as comfrey.
Sinusoidal endothelial cells are more susceptible than hepatocytes to PAs that cause VOD/SOS . Functional annotation extracted from Ingenuity Pathway Analysis revealed that many of the transcriptional responses were associated with the apoptosis, cell death, adhesion, and cell movement of endothelial cells (Table 3). The genes in these pathways were highly expressed in comfrey-treated livers; these genes included the endothelial cell markers plasminogen activator inhibitor type 1 (PAI-1, also known as Serpine1 and Serpine2) and tissue plasminogen activator (Plat), and the cytokine tumor necrosis factors (Tnfrsf6, Tnfrsf12a, and Tnfsf10), and transforming growth factors (Tgfb1 and Tgfb2). Induction of a number of genes involved in the injury of endothelial cells was also detected, including endothelin 1 (Edn1), urokinase plasminogen activator receptor (Plaur), collagen type IV alpha 2 (Col4a2), matrix metalloproteinase 2 (Mmp2), mitogen-activated protein kinase 9 (Mapk9), and secreted phosphoprotein 1 (Spp1). It has been reported that Edn1 is a mediator of hepatic sinusoidal constriction, and increased activity of matrix metalloproteinases is responsible for changes of sinusoidal endothelial cells . Elevated plasma PAI-1 levels are useful in distinguishing VOD/SOS . Endothelial injury is the initiating event in the cascade of events leading to the hepatic changes and clinical manifestation of VOD/SOS . Our results offer a more comprehensive overview of the molecular responses to comfrey exposure by expression of multiple genes in liver endothelial cells.
PAs in comfrey can reach the hepatocytes via the sinusoidal blood, and their toxic metabolites lead to immediate damage to the hepatocytes . Genes involved in liver cell death and growth were also induced or repressed in response to comfrey treatment (Table 4), including tumor necrosis factors, transforming growth factor β, chemokine receptor CCR2 gene (Ccr2), heme oxygenase 1 (Hmox1), immediate early response 3 (Ier3), cyclin-dependent kinase inhibitor (Cdkn1a), inhibin beta A (Inhba), and tissue inhibitor of metalloproteinase 1 (Timp1). It is well known that Hmox1 induction is a protective mechanism against the oxidative stress associated with liver injury , and that elevated Hgf could protect hepatocytes from injury or promote hepatocellular regeneration . Since Inhba, so-called activin A, is a negative regulator of hepatocyte cell growth , the decreased expression of Inhba observed after comfrey treatment suggests the induction of hepatic growth.
Necrosis of hepatocytes and mesenchymal cells follows comfrey-induced liver cell injury, and functional cells are replaced by fibrotic tissues [27, 31, 32]. In the present study, we observed 16 genes involved in the function of liver fibrosis by Ingenuity Pathway Analysis (Table 4). Liver fibrosis is characterized by cell proliferation and the accumulation of extracellular matrix components and is mediated by cytokines and growth factors, of which TGF-β1 appears to be a key mediator . The up-regulation of cytokines Hgf and Tnfrsf6 and down-regulation of Igf1 and Lif play an important role in the pathogenesis of liver injury and fibrosis. Decreased serum Igf-1 levels provide a useful index of hepatocellular dysfunction and impaired nutritional status, and increased Hgf appears to limit liver fibrosis . Agt and Agtr1a, cytokines with vasoactive properties, also regulate liver fibrogenesis. Chemokines have a much wider biological role including angiogenesis, carcinogenesis, and cell cycle control . Chemokines in the liver (Cxcl4 and Cxcl12) may modulate the progression of liver fibrosis through their actions on hepatic stellate cells.
PA-induced DNA damage in the liver (endothelial cells and hepatocytes), if not repaired prior to DNA synthesis, might produce replication errors and mutations, which eventually could result in the development of neoplasmas in the treated animals. We determined MFs in the liver cII gene of Big Blue transgenic rats. After feeding with 8% comfrey root for 12 weeks, we observed a 4-fold higher MF in the liver cII gene compared to the controls (Figure 2). The induction of mutation was similar to that reported previously for rats fed with 2% comfrey root . These observations suggest that the rats could not tolerate the feeding of roots in concentrations over 2%, in terms of mutation induction. Furthermore, the overall pattern of mutations induced by 8% comfrey in liver was similar to that in the livers of rats fed 2% comfrey root (Table 1), whereas both the 2% and 8% comfrey-induced mutation spectra were significantly different from liver controls. In contrast to the G:C → A:T transition that was the predominant mutation in the controls, the major type of mutation in the 8% comfrey-fed rats was G:C → T:A transversion (41%), a mutation that was also induced by riddelliine, a representative genotoxic PA that is tumorigenic for rat liver . In addition, 13% of mutations from the 8% comfrey-fed rats were tandem base substitution, which has been suggested as a mutational signature for the genetic damage of PAs . G:C → T:A transversion may cause the initiation of tumors in the liver of rats fed with comfrey, because it has been reported that more than half of riddelliine-induced liver hemangiosarcomas have a G → T mutation at K-ras codon 12 . p53 mutation also has been detected at an early stage of riddelliine exposure . Mutations are thought to be involved in carcinogenesis because the transition from a normal somatic cell to a cancer cell is due to mutations in protooncogenes, tumor suppressor genes and/or genes that function in the maintenance of genomic stability [39, 40]. Comfrey-induced neoplasms in the rat were mostly found in the liver. Hepatocellular adenomas were induced in all experimental groups that received diets containing 1–8% comfrey root or 8–33% comfrey leaves. In addition, a few rats bearing hepatocellular adenomas simultaneously had hemangioendothelial sarcoma of the liver . The first hepatocellular carcinoma (HCC) appearance was 7 months after initiating the 8% comfrey diet.
Ingenuity Pathway Analysis found 8 genes (2 down- and 6 up -regulated) involved in liver cancer development that were altered due to comfrey treatment (Table 5). In mammals, the nucleotide excision repair process is the most important repair pathway for elimination of DNA damage caused by exogenous agents, including UV light, DNA-reactive carcinogens, and some endogenously generated oxidative lesions . Xeroderma pigmentosum group C (Xpc) is implicated in the early steps of this repair pathway. A significantly higher incidence of chemically induced liver and lung tumors is observed in Xpc null mice . Gap junction membrane channel protein beta 1 (Gjb1), also called connexin 32 (Cx32), is the main gap junction protein in hepatocytes and plays an important role in the regulation of signal transfer and growth control in the liver. It has been reported that Cx32 expression decreases gradually as liver disease progresses to cirrhosis and HCC , and a low expression of Cx32 mRNAs in HCC tissues is also predictive of the postoperative recurrence of HCCs . Bcl-2 is characterized as an antiapoptotic/oncogenic protein and also functions as an antioxidant. Increased Bcl-2 expression in cirrhotic patients correlates with the development of HCC . Bcl-2 is also expressed in HCC tissues and the increasing Bcl-2 expression associated with HCC progression suggests that the Bcl-2 protein takes part in the formation of HCC . Hgf, identified originally as the most potent mitogen for hepatocytes, is now known to be a cytokine with numerous functions in a wide variety of cells . It is up-regulated in inflammatory liver diseases and stimulates DNA synthesis preferentially in initiated hepatocytes, presumably resulting in tumour promotion . N-myc downstream-regulated gene 1 (Ndr1, or Ndrg1) plays a role in growth arrest and cell differentiation, is induced by several stress conditions, and is overexpressed in many cancers . Expression of Ndrg1 was significantly up-regulated in HCC tissues compared to that of noncancerous and normal liver tissues . Cyclin E (Ccne1), a regulatory subunit of cyclin-dependent kinase 2, is an important regulator for entry into the S phase of the mammalian cell cycle. Overexpression of Ccne1 has been observed in many tumors including primary HCCs ; overexpression results in chromosome instability and thus may contribute to tumorigenesis . Growth arrest and DNA damage 45 alpha (Gadd45a) is a nuclear protein involved in the maintenance of genomic stability, DNA repair, and the suppression of cell growth . Gadd45a protein levels are higher in liver cirrhotic and neoplastic tissues . Tissue inhibitor of metalloproteinase 1 (Timp1) is a contributory factor to fibrosis of a variety of organs including the liver. Timp1 and other extracellular matrix remodeling genes are implicated in the transition from mild to moderate fibrosis in patients with chronic hepatitis C .
The integration of gene expression changes with a mechanistic pathway analysis suggests a scheme for comfrey treatment leading to tumorigenesis via mutation induction (Figure 5). The available evidence suggests that active metabolites of PAs in comfrey interact with endothelial and hepatocyte DNA, causing damage to hepatic endothelial cells and hepatocytes. This may result in liver fibrosis and increases in mutation induction, which may be associated with the development of HCC and hemagiosarcoma. We have identified 2,726 genes in the livers of comfrey-fed rats that were differentially expressed. Some of the gene changes are associated with the metabolism, injury of endothelial cells, and liver injury and abnormalities that are postulated to occur as a result of comfrey exposure. This approach provides further insight into the mechanisms involved in the development of VOD/SOS and tumorigenesis after exposure to comfrey.
Materials and methods
Plant material and animals
Comfrey roots (Symphytum officinale) were purchased from Camas Prairie Products (Trout Lake, WA). Male Big Blue Fisher 344 transgenic rats were obtained from Taconic Laboratories (Germantown, NY) through purchase from Stratagene (La Jolla, CA). All animal procedures followed the recommendations of the NCTR Institutional Animal Care and Use Committee for the handling, maintenance, treatment, and sacrifice of the rats.
Comfrey diet and treatments
The comfrey roots were ground into powder using a Wiley Mill, and the comfrey root powder was stored at room temperature until use. NIH-31 pellets (Purina Mills International, Brentwood, MO) were autoclaved and ground into meal form. The base diet was blended with comfrey root powder in a Hobart Mixer to make 8% comfrey root dosed diet, which was provided in rat feeders. The treatment schedule was based on the protocol used in a carcinogenesis assay . Male, 6-week-old Big Blue rats were fed without (vehicle control group) or with 8% comfrey roots for 12 weeks. Six rats from each treatment group were sacrificed at the end of the treatment. The livers were isolated, frozen quickly in liquid nitrogen, and stored at -80°C.
cII mutation assay
High-molecular-weight genomic DNA was extracted from rat livers using the RecoverEase DNA Isolation Kit (Stratagene) and stored at 4°C until DNA packaging was performed. The packaging of the phage, plating the packaged DNA samples, and determination of MF were carried out following the manufacturer's instructions for the λ Select-cII Mutation Detection System for Big Blue Rodents (Stratagene). The shuttle vector containing the cII target gene was rescued from total genomic DNA with phage packaging extract (Transpack, Stratagene). The plating was performed with the Escherichia coli host strain G1250. To determine the total titer of packaged phages, G1250 bacteria were mixed with 1:3000 dilutions of phage, plated on TB1 plates, and incubated overnight at 37°C (nonselective conditions). For mutant selection, the packaged phages were mixed with G1250, plated on TB1 plates, and incubated at 24°C for about 42 h (conditions for cII-selection). Under these conditions, phages with wild-type cII genes undergo lysogenization and become part of the developing bacterial lawn, whereas phages with mutated cII genes undergo lytic growth and give rise to plaques. When incubated at 37°C, phages with wild-type cII genes also undergo a lytic cycle, resulting in plaque formation. Assays were repeated until a minimum of 2 × 105 plaque-forming units from each sample were examined for mutation. The cII MF is defined as the total number of mutant plaques (determined at 24°C) divided by the total number of plaques screened (determined at 37°C).
Sequence analysis of the cII mutants
The mutants were sequenced according to the method of Mei et al. . The cII mutant plaques were selected at random from different animals and replated at low density to verify the mutant phenotype. Single, well-isolated plaques were selected from these plates and transferred to a microcentrifuge tube containing 100 μl of sterile distilled water. The tube was heated at 100°C for 5 min and centrifuged at 12,000 g for 3 min. The cII target DNA for sequencing was amplified by PCR using primers 5'-AAAAAGGGCATCAAATTAACC-3' (upstream) and 5'-CCGAAGTTGAGTATTTTTGCTG-3' (downstream). For PCR amplification, 10 μl of the supernatant were added to 10 μl of a PCR Master Mix (Promega, Madison, WI) and the primers. The final concentrations of the reagents were: 1× Taq polymerase reaction buffer, 0.2 μM of each primer, 200 μM of each dNTP, 1.5 mM MgCl2, and 0.25 U of Taq DNA polymerase. The PCR reactionwas performed using a PCR System 9700 (Applied Biosystems, Foster City, CA), with the following cycling parameters: a 3 min denaturation at 95°C, followed by 35 cycles of 30 s at 95°C, 1 min at 60°C, and 1 min at 72°C, with a final extension of 10 min at 72°C. The PCR products were isolated using a PCR purification kit (Qiagen, Chatsworth, CA). The cII mutant DNA was sequenced with a CEQ Dye Terminator Cycle Sequencing Kit and a CEQ 8000 Genetic Analysis System (Beckman Coulter, Fullerton, CA). The primer for cII mutation sequencing was the upstream primer used for the PCR.
Statistical analyses for mutagenicity data
Analyses were performed using the SigmaStat 2.03 program (SPSS, Chicago, IL). All of the MF data were expressed as the mean ± standard deviation (SD) from 6 rats per group. Statistical significance was determined by one-way analysis of variance (ANOVA) followed by the Holm-Sidak test. Mutational spectra were compared using the computer program written by Cariello and colleagues  for the Monte Carlo analysis developed by Adams and Skopek .
RNA isolation and quality control
Total RNA was isolated from liver tissues of 6 control and 6 comfrey-fed rats using an RNeasy system (Qiagen). The yield of the extracted RNA was determined spectrophotometrically by measuring the optical density at 260 nm. The purity and quality of extracted RNA were evaluated using the RNA 6000 LabChip and Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA). RNA samples with RNA integrity numbers (RINs) greater than 7.5 were used for microarray experiments performed using Applied Biosystems' Rat Genome Survey Microarray platform, which is a one channel microarray with chemiluminescence detection, and contains 26,857 probes (60-mer) for the interrogation of 27,088 genes and 1592 controls that track system performance through each experiment.
Preparation of digoxigenin labeled in vitro transcribed cRNA
All RNA targets were labeled using the Applied Biosystems RT-IVT Labeling Kit Version 2.0. Briefly, 1.5 μg of total RNA was reverse transcribed via 2 h incubation at 42°C with ArrayScript RT enzyme (Ambion, Austin, TX) and oligo dT-T7 primer. Double stranded cDNA was produced following 2 h incubation with E. coli DNA polymerase and RNase H at 16°C. Double-stranded cDNA was purified according to the RT-IVT kit protocol. In vitro transcription was performed by incubation of the cDNA product with T7 RNA polymerase, 0.75 mM Digoxigenin-11-UTP (Roche Applied Science, Indianapolis, IN) and all other NTPs for 9 h. Labeled cRNA was purified according to the RT-IVT kit protocol and analyzed for quality and quantity using standard UV spectrometry and the Bioanalyzer.
Hybridization of labeled cRNA to microarrays and microarray imaging
Digoxigenin labeled cRNA targets were hybridized to Applied Biosystems Rat Whole Genome Survey Microarrays using the Applied Biosystems Chemiluminescent Detection Kit. Briefly, 15 μg of labeled cRNA targets were fragmented via incubation with fragmentation buffer provided in the kit for 30 min at 60°C. Fragmented targets were hybridized to microarrays during a 16 h incubation at 55°C with buffers and reagents from the Chemiluminescent Detection Kit. Post-hybridization washes and anti-Digoxigenin-Alkaline Phosphatase binding were performed according to the protocol of the kit. Chemiluminescence detection, image acquisition and analysis were performed using Applied Biosystems Chemiluminescence Detection Kit and Applied Biosystems 1700 Chemiluminescent Microarray Analyzer following the manufacturer's protocols. Images were auto-gridded and the chemiluminescent signals were quantified, corrected for background, and finally, spot- and spatially-normalized using the Applied Biosystems 1700 Chemiluminescent Microarray Analyzer software version 1.1.
Microarray data analysis
Gene expression data from the Applied Biosystems' Rat Genome Survey Microarray were input to ArrayTrack, a software system developed by the FDA's National Center for Toxicological Research for the management, analysis, visualization and interpretation of microarray data . Raw microarray intensity data were normalized per chip to the same median intensity value of 500. Chemiluminescent signals from 1529 control probes that track system performance through each experiment were not used in normalization. The identification of differentially expressed genes based on t-tests and fold-change cutoffs, and Principal Component Analysis were conducted within ArrayTrack. Ingenutity Pathway Analysis (Mountain View, CA) was used for pathway and function analysis.
Winship KA: Toxicity of comfrey. Adverse Drug React Toxicol Rev 1991, 10(1):47–59.
Rode D: Comfrey toxicity revisited. Trends Pharmacol Sci 2002, 23(11):497–499. 10.1016/S0165-6147(02)02106-5
Fu PP, Xia Q, Lin G, Chou MW: Pyrrolizidinealkaloids – genotoxicity, metabolism enzymes, metabolic activation, and mechanisms. Drug Metab Rev 2004, 36(1):1–55. 10.1081/DMR-120028426
Prakash AS, Pereira TN, Reilly PE, Seawright AA: Pyrrolizidine alkaloids in human diet. Mutat Res 1999, 443(1–2):53–67.
USFDA, Center for Food Safety and Applied Nutrition: FDA advises dietary supplement manufacturers to remove comfrey products from the market. FDA 2001. [http://vm.cfsan.fda.gov/~dms/dspltr06.html]
Weston CF, Cooper BT, Davies JD, Levine DF: Veno-occlusive disease of the liver secondary to ingestion of comfrey. Br Med J (Clin Res Ed) 1987, 295(6591):183.
Ridker PM, Ohkuma S, McDermott WV, Trey C, Huxtable RJ: Hepatic venocclusive disease associated with the consumption of pyrrolizidine-containing dietary supplements. Gastroenterology 1985, 88(4):1050–1054.
DeLeve LD, Shulman HM, McDonald GB: Toxic injury to hepatic sinusoids: sinusoidal obstruction syndrome (veno-occlusive disease). Semin Liver Dis 2002, 22(1):27–42. 10.1055/s-2002-23204
Hirono I, Mori H, Haga M: Carcinogenic activity of Symphytum officinale. J Natl Cancer Inst 1978, 61(3):865–869.
Mei N, Guo L, Fu PP, Heflich RH, Chen T: Mutagenicity of comfrey (Symphytum Officinale) in rat liver. Br J Cancer 2005, 92(5):873–875. 10.1038/sj.bjc.6602420
DeLeve LD, Wang X, Kanel GC, Ito Y, Bethea NW, McCuskey MK, Tokes ZA, Tsai J, McCuskey RS: Decreased hepatic nitric oxide production contributes to the development of rat sinusoidal obstruction syndrome. Hepatology 2003, 38(4):900–908. 10.1053/jhep.2003.50383
Paradis V, Bieche I, Dargere D, Cazals-Hatem D, Laurendeau I, Saada V, Belghiti J, Bezeaud A, Vidaud M, Bedossa P, et al.: Quantitative gene expression in Budd-Chiari syndrome: a molecular approach to the pathogenesis of the disease. Gut 2005, 54(12):1776–1781. 10.1136/gut.2005.065144
Kim NC, Oberlies NH, Brine DR, Handy RW, Wani MC, Wall ME: Isolation of symlandine from the roots of common comfrey (Symphytum officinale) using countercurrent chromatography. J Nat Prod 2001, 64(2):251–253. 10.1021/np0004653
Schaneberg BT, Molyneux RJ, Khan IA: Evaporative light scattering detection of pyrrolizidine alkaloids. Phytochem Anal 2004, 15(1):36–39. 10.1002/pca.715
Stickel F, Seitz HK: The efficacy and safety of comfrey. Public Health Nutr 2000, 3(4A):501–508.
Couet CE, Crews C, Hanley AB: Analysis, separation, and bioassay of pyrrolizidine alkaloids from comfrey (Symphytum officinale). Nat Toxins 1996, 4(4):163–167.
Huxtable RJ, Luthy J, Zweifel U: Toxicity of comfrey-pepsin preparations. N Engl J Med 1986, 315(17):1095.
Bach N, Thung SN, Schaffner F: Comfrey herb tea-induced hepatic veno-occlusive disease. Am J Med 1989, 87(1):97–99. 10.1016/S0002-9343(89)80492-9
Ridker PM, McDermott WV: Comfrey herb tea and hepatic veno-occlusive disease. Lancet 1989, 1(8639):657–658. 10.1016/S0140-6736(89)92154-5
Yeong ML, Swinburn B, Kennedy M, Nicholson G: Hepatic veno-occlusive disease associated with comfrey ingestion. J Gastroenterol Hepatol 1990, 5(2):211–214.
Yeong ML, Clark SP, Waring JM, Wilson RD, Wakefield SJ: The effects of comfrey derived pyrrolizidine alkaloids on rat liver. Pathology 1991, 23(1):35–38.
DeLeve LD, McCuskey RS, Wang X, Hu L, McCuskey MK, Epstein RB, Kanel GC: Characterization of a reproducible rat model of hepatic veno-occlusive disease. Hepatology 1999, 29(6):1779–1791. 10.1002/hep.510290615
Svoboda DJ, Reddy JK: Malignant tumors in rats given lasiocarpine. Cancer Res 1972, 32(5):908–913.
Hirono I, Haga M, Fujii M, Matsuura S, Matsubara N, Nakayama M, Furuya T, Hikichi M, Takanashi H, Uchida E, et al.: Induction of hepatic tumors in rats by senkirkine and symphytine. J Natl Cancer Inst 1979, 63(2):469–472.
Zhou S, Koh HL, Gao Y, Gong ZY, Lee EJ: Herbal bioactivation: the good, the bad and the ugly. Life Sci 2004, 74(8):935–968. 10.1016/j.lfs.2003.09.035
Salat C, Holler E, Kolb HJ, Reinhardt B, Pihusch R, Wilmanns W, Hiller E: Plasminogen activator inhibitor-1 confirms the diagnosis of hepatic veno-occlusive disease in patients with hyperbilirubinemia after bone marrow transplantation. Blood 1997, 89(6):2184–2188.
Kumar S, DeLeve LD, Kamath PS, Tefferi A: Hepatic veno-occlusive disease (sinusoidal obstruction syndrome) after hematopoietic stem cell transplantation. Mayo Clin Proc 2003, 78(5):589–598.
Bauer M, Bauer I: Heme oxygenase-1: redox regulation and role in the hepatic response to oxidative stress. Antioxid Redox Signal 2002, 4(5):749–758. 10.1089/152308602760598891
Oe H, Kaido T, Mori A, Onodera H, Imamura M: Hepatocyte growth factor as well as vascular endothelial growth factor gene induction effectively promotes liver regeneration after hepatectomy in Solt-Farber rats. Hepatogastroenterology 2005, 52(65):1393–1397.
Gold EJ, Francis RJ, Zimmermann A, Mellor SL, Cranfield M, Risbridger GP, Groome NP, Wheatley AM, Fleming JS: Changes in activin and activin receptor subunit expression in rat liver during the development of CCl4-induced cirrhosis. Mol Cell Endocrinol 2003, 201(1–2):143–153. 10.1016/S0303-7207(02)00417-3
Yeong ML, Wakefield SJ, Ford HC: Hepatocyte membrane injury and bleb formation following low dose comfrey toxicity in rats. Int J Exp Pathol 1993, 74(2):211–217.
Bataller R, Brenner DA: Liver fibrosis. J Clin Invest 2005, 115(2):209–218. 10.1172/JCI200524282
Sanz S, Pucilowska JB, Liu S, Rodriguez-Ortigosa CM, Lund PK, Brenner DA, Fuller CR, Simmons JG, Pardo A, Martinez-Chantar ML, et al.: Expression of insulin-like growth factor I by activated hepatic stellate cells reduces fibrogenesis and enhances regeneration after liver injury. Gut 2005, 54(1):134–141. 10.1136/gut.2003.024505
Marra F: Chemokines in liver inflammation and fibrosis. Front Biosci 2002, 7: d1899–1914.
Mei N, Chou MW, Fu PP, Heflich RH, Chen T: Differential mutagenicity of riddelliine in liver endothelial and parenchymal cells of transgenic big blue rats. Cancer Lett 2004, 215(2):151–158. 10.1016/j.canlet.2004.06.013
Mei N, Heflich RH, Chou MW, Chen T: Mutations induced by the carcinogenic pyrrolizidine alkaloid riddelliine in the liver cII gene of transgenic big blue rats. Chem Res Toxicol 2004, 17(6):814–818. 10.1021/tx049955b
Hong HL, Ton TV, Devereux TR, Moomaw C, Clayton N, Chan P, Dunnick JK, Sills RC: Chemical-specific alterations in ras, p53, and beta-catenin genes in hemangiosarcomas from B6C3F1 mice exposed to o-nitrotoluene or riddelliine for 2 years. Toxicol Appl Pharmacol 2003, 191(3):227–234. 10.1016/S0041-008X(03)00165-0
Nyska A, Moomaw CR, Foley JF, Maronpot RR, Malarkey DE, Cummings CA, Peddada S, Moyer CF, Allen DG, Travlos G, et al.: The hepatic endothelial carcinogen riddelliine induces endothelial apoptosis, mitosis, S phase, and p53 and hepatocytic vascular endothelial growth factor expression after short-term exposure. Toxicol Appl Pharmacol 2002, 184(3):153–164. 10.1006/taap.2002.9485
Hennings H, Glick AB, Greenhalgh DA, Morgan DL, Strickland JE, Tennenbaum T, Yuspa SH: Critical aspects of initiation, promotion, and progression in multistage epidermal carcinogenesis. Proc Soc Exp Biol Med 1993, 202(1):1–8.
Loeb LA: A mutator phenotype in cancer. Cancer Res 2001, 61(8):3230–3239.
Friedberg EC: How nucleotide excision repair protects against cancer. Nat Rev Cancer 2001, 1(1):22–33. 10.1038/35094000
Cheo DL, Burns DK, Meira LB, Houle JF, Friedberg EC: Mutational inactivation of the xeroderma pigmentosum group C gene confers predisposition to 2-acetylaminofluorene-induced liver and lung cancer and to spontaneous testicular cancer in Trp53-/-mice. Cancer Res 1999, 59(4):771–775.
Nakashima Y, Ono T, Yamanoi A, El-Assal ON, Kohno H, Nagasue N: Expression of gap junction protein connexin32 in chronic hepatitis, liver cirrhosis, and hepatocellular carcinoma. J Gastroenterol 2004, 39(8):763–768. 10.1007/s00535-003-1386-2
Sheen IS, Jeng KS, Wang PC, Shih SC, Chang WH, Wang HY, Chen CC, Shyung LR: Are gap junction gene connexins 26, 32 and 43 of prognostic values in hepatocellular carcinoma? A prospective study. World J Gastroenterol 2004, 10(19):2785–2790.
Frommel TO, Yong S, Zarling EJ: Immunohistochemical evaluation of Bcl-2 gene family expression in liver of hepatitis C and cirrhotic patients: a novel mechanism to explain the high incidence of hepatocarcinoma in cirrhotics. Am J Gastroenterol 1999, 94(1):178–182. 10.1111/j.1572-0241.1999.00792.x
Duan XX, Ou JS, Li Y, Su JJ, Ou C, Yang C, Yue HF, Ban KC: Dynamic expression of apoptosis-related genes during development of laboratory hepatocellular carcinoma and its relation to apoptosis. World J Gastroenterol 2005, 11(30):4740–4744.
Kaido T, Imamura M: Hepatocyte growth factor: clinical implications in hepatobiliary pancreatic surgery. J Hepatobiliary Pancreat Surg 2001, 8(1):65–75. 10.1007/s005340170052
Drucker C, Parzefall W, Teufelhofer O, Grusch M, Ellinger A, Schulte-Hermann R, Grasl-Kraupp B: Non-parenchymal liver cells support the growth advantage in the first stages of hepatocarcinogenesis. Carcinogenesis 2006, 27(1):152–161. 10.1093/carcin/bgi202
Cangul H: Hypoxia upregulates the expression of the NDRG1 gene leading to its overexpression in various human cancers. BMC Genet 2004, 5: 27. 10.1186/1471-2156-5-27
He J, Zhou R: [Correlation of NDRG1 gene with liver tissue differentiation and hepatocarcinogenesis]. Beijing Da Xue Xue Bao 2003, 35(5):471–475.
Yasui K, Arii S, Zhao C, Imoto I, Ueda M, Nagai H, Emi M, Inazawa J: TFDP1, CUL4A, and CDC16 identified as targets for amplification at 13q34 in hepatocellular carcinomas. Hepatology 2002, 35(6):1476–1484. 10.1053/jhep.2002.33683
Spruck CH, Won KA, Reed SI: Deregulated cyclin E induces chromosome instability. Nature 1999, 401(6750):297–300. 10.1038/45836
Hollander MC, Sheikh MS, Bulavin DV, Lundgren K, Augeri-Henmueller L, Shehee R, Molinaro TA, Kim KE, Tolosa E, Ashwell JD, et al.: Genomic instability in Gadd45a-deficient mice. Nat Genet 1999, 23(2):176–184. 10.1038/13802
Gramantieri L, Chieco P, Giovannini C, Lacchini M, Trere D, Grazi GL, Venturi A, Bolondi L: GADD45-alpha expression in cirrhosis and hepatocellular carcinoma: relationship with DNA repair and proliferation. Hum Pathol 2005, 36(11):1154–1162. 10.1016/j.humpath.2005.07.017
Asselah T, Bieche I, Laurendeau I, Paradis V, Vidaud D, Degott C, Martinot M, Bedossa P, Valla D, Vidaud M, et al.: Liver gene expression signature of mild fibrosis in patients with chronic hepatitis C. Gastroenterology 2005, 129(6):2064–2075. 10.1053/j.gastro.2005.09.010
Cariello NF: Software for the analysis of mutations at the human hprt gene. Mutat Res 1994, 312(2):173–185.
Adams WT, Skopek TR: Statistical test for the comparison of samples from mutational spectra. J Mol Biol 1987, 194(3):391–396. 10.1016/0022-2836(87)90669-3
Fang H, Tong W, Perkins R, Shi L, Hong H, Cao X, Xie Q, Yim SH, Ward JM, Pitot HC, et al.: Bioinformatics approaches for cross-species liver cancer analysis based on microarray gene expression profiling. BMC Bioinformatics 2005, 6(Suppl 2):S6. [http://] 10.1186/1471-2105/6/S2/S6
We thank Drs. Robert H. Heflich and Tao Han from NCTR for their helpful discussions, comments, and criticisms.
The views presented in this article do not necessarily reflect those of the Food and Drug Administration.
NM performed the animal treatment and mutagenicity analysis, was involved in the analysis of microarray data, and wrote the manuscript. LG, LS, and JCF helped conceive the experiments, analysis of the data, and writing of the manuscript. LZ, YAS, and CF conducted the microarray experiment and generated the raw data. CLM and SLD performed technical support. TC had the original idea for this study, and was involved in designing the experiment and writing the manuscript. All authors approved the final version of manuscript.
About this article
Cite this article
Mei, N., Guo, L., Zhang, L. et al. Analysis of gene expression changes in relation to toxicity and tumorigenesis in the livers of Big Blue transgenic rats fed comfrey (Symphytum officinale). BMC Bioinformatics 7, S16 (2006). https://doi.org/10.1186/1471-2105-7-S2-S16
- Liver Fibrosis
- Ingenuity Pathway Analysis
- Sinusoidal Endothelial Cell
- Liver Cell Death