Skip to content

Advertisement

Volume 8 Supplement 7

Proceedings of the Fourth Annual MCBIOS Conference. Computational Frontiers in Biomedicine

  • Proceedings
  • Open Access

Comparison of gene expression profiles altered by comfrey and riddelliine in rat liver

BMC Bioinformatics20078 (Suppl 7) :S22

https://doi.org/10.1186/1471-2105-8-S7-S22

©  Guo et al; licensee BioMed Central Ltd. 2007

  • Published:

Abstract

Background

Comfrey (Symphytum officinale) is a perennial plant and has been consumed by humans as a vegetable, a tea and an herbal medicine for more than 2000 years. It, however, is hepatotoxic and carcinogenic in experimental animals and hepatotoxic in humans. Pyrrolizidine alkaloids (PAs) exist in many plants and many of them cause liver toxicity and/or cancer in humans and experimental animals. In our previous study, we found that the mutagenicity of comfrey was associated with the PAs contained in the plant. Therefore, we suggest that carcinogenicity of comfrey result from those PAs. To confirm our hypothesis, we compared the expression of genes and processes of biological functions that were altered by comfrey (mixture of the plant with PAs) and riddelliine (a prototype of carcinogenic PA) in rat liver for carcinogenesis in this study.

Results

Groups of 6 Big Blue Fisher 344 rats were treated with riddelliine at 1 mg/kg body weight by gavage five times a week for 12 weeks or fed a diet containing 8% comfrey root for 12 weeks. Animals were sacrificed one day after the last treatment and the livers were isolated for gene expression analysis. The gene expressions were investigated using Applied Biosystems Rat Whole Genome Survey Microarrays and the biological functions were analyzed with Ingenuity Analysis Pathway software. Although there were large differences between the significant genes and between the biological processes that were altered by comfrey and riddelliine, there were a number of common genes and function processes that were related to carcinogenesis. There was a strong correlation between the two treatments for fold-change alterations in expression of drug metabolizing and cancer-related genes.

Conclusion

Our results suggest that the carcinogenesis-related gene expression patterns resulting from the treatments of comfrey and riddelliine are very similar, and PAs contained in comfrey are the main active components responsible for carcinogenicity of the plant.

Keywords

  • Ingenuity Pathway Analysis
  • Pyrrolizidine Alkaloid
  • Gene Expression Alteration
  • Main Active Component
  • Gene Encode Heat Shock Protein

Background

Comfrey (Symphytum officinale) is consumed by humans as a vegetable and a tea. It has been used as an herbal medicine for more than 2000 years to treat broken bones, tendon damage, ulcerations in the gastrointestinal tract, lung congestion, and joint inflammation, and to promote wound healing [1]. It, however, has been reported that comfrey is hepatotoxic in livestock and humans, and carcinogenic in experimental animals. Comfrey induced hepatic veno-occlusive lesion (VOD) in humans [2–4] and hepatocellular adenomas and hemangioendothelial sarcomas in rats [5]. Therefore, the regular use of comfrey is a potential health risk for development of liver cancers. In 2001, the US Food and Drug Administration requested voluntary compliance for the removal of products containing comfrey [6].

It is still not clear about the mechanism of tumor induction by comfrey because comfrey is a mixture of many different substances and the active components responsible for the carcinogenesis have not been identified. It has been suggested that the induction of hepatic tumors has been associated with the pyrrolizidine alkaloids (PAs) that are present in comfrey [7–9] since PAs are genotoxic and carcinogenic in liver [10, 11]. Recently, we demonstrated that the PAs in the comfrey plant appear to be responsible for mutation induction in rat liver [12]. Mutations are involved in the etiology of cancer [13]. Research results in molecular cancer genetics have identified inherited and somatic cell mutations associated with cancer in oncogenes, tumor suppressor genes, DNA repair genes and other related genes [14–17]. Therefore, we hypothesize that PAs contained in comfrey are the main active components resulting in tumors in liver.

Riddelliine is one of the tumorigenic PAs and has been studied as a prototype of PA. The toxicity and carcinogenicity of riddelliine have been studied by the National Toxicological Program (NTP) [18–20], and the mechanism of riddelliine-induced tumorigenicity in experimental animals has been studied at the National Center for Toxicological Research (NCTR) [21–27]. Results showed that (1) riddelliine was metabolized to the major metabolites 6,7-dihydro-1-hydroxymethyl-5H-pyrrolizine (DHP) and riddelliine N-oxide; (2) DHP-derived DNA adducts were formed both in vivo and in vitro; (3) G:C → T:A transversions were the major type of mutation induced in the liver of riddelliine-treated rats; and (4) riddelliine caused liver tumors in male mice and both sexes of rats, mononuclear cell leukemia in rats, and lung neoplasms in female mice.

DNA microarray, a key advanced technology, has developed rapidly because of its ability to examine the expression levels of thousands of genes simultaneously. It is used increasingly for identifying biomarkers, elucidating patterns of gene expression, and understanding the mechanism of disease and toxicity. Toxicogenomics applies high through-put genomics tools to the study of toxicology and gene expression microarrays have been used extensively [28]. The interpretation of gene expression data, however, can be complicated by the gene expression alterations caused by environmental factors such as diet [29] and time of day [30, 31]. That is, it can be difficult to discern toxin-specific gene expression changes from those due to environmental effects. This can be complicated further when the test article is a complex food such as comfrey.

To examine the possible PA-induced affects on gene expression caused by comfrey, the gene expression profiles in the livers of comfrey-treated rats were compared to the gene expression profiles from rats treated with the purified PA riddelliine. The correlation of the gene expression and biological functions related to carcinogenesis between the two treatments were explored.

Results and discussion

Previously, we investigated the mutagenicity of comfrey and riddelliine in rats [26, 32] under conditions which have been shown to result in liver tumors [5, 20]. In the original mutation studies, 8% comfrey treatment induced about a 4-fold higher mutant frequency (MF) over the control group [32], while 1 mg/kg riddelliine exposure resulted in about a 3-fold increase in MF compared to controls [26]. In these studies, the rats were exposed to the compounds for 12 weeks and then sacrificed. No tumors were visible at this time in our study, although Hirono et al. reported that liver tumors developed in the rats as early as 7 weeks after initiation of feeding comfrey [5]. The mutation spectra for comfrey- and riddelliine-treated rats were significantly different from the controls, but there was no significant difference between the spectra for comfrey and riddelliine. G:C → T:A transversion (35–42%) was the major type of mutation in both comfrey and riddelliine treated rats, with 8–17% tandem base substitutions [12]. These mutational data from comfrey-treated rats suggests that PAs in the plant are responsible for mutation induction and tumor initiation in rat liver.

To further investigate the effects of comfrey and riddelliine on gene expression changes in liver, microarray analysis was performed on liver samples from these high-dose groups. The details of each experiment and gene expression analysis were reported previously [32, 33]. In the present study, gene expression alterations caused by exposure to comfrey were compared to those caused by exposure to the representative PA riddelliine.

Comparison of gene expression profiles induced by comfrey and riddelliine

The intensities of the whole rat gene data were normalized by quantile normalization and then analyzed by Hierarchical Cluster Analysis (HCA) and Principal Component Analysis (PCA) to visualize clusters of samples corresponding to the different treatments (Figures 1 and 2). The result of the HCA of the gene expression data is shown in Figure 1, and shows that samples were grouped together according to the treatments. There were clear separations between of the control, comfrey-treated, and riddelliine-treated groups. The gene expression pattern due to comfrey exposure was distinct from that due to riddelliine exposure (Figure 2), possibly reflecting the effects of many other substances in comfrey other than PAs [9]. This difference is expected to be the result of different properties of two chemicals including pharmacological and toxicological effects.
Figure 1
Figure 1

Hierarchical Cluster Analysis (HCA) of expression profiles for control, 8% comfrey-fed, and 1 mg/kg riddelliine treated groups. The log2 intensity of the entire gene set was scaled by Z-score transformation, and then these values were hierarchically clustered using Euclidean distance metric and average linkage. Each column represents the results from an individual animal. CTR, control; CFY, comfrey treatment; RDL, riddelliine treatment.

Figure 2
Figure 2

Principal Component Analysis (PCA) of expression profiles for control, 8% comfrey-fed, and 1 mg/kg riddelliine treated groups. The intensity of the entire gene set was used, and no specific cut off was applied.

The differentially expressed genes between the treatment and control groups were identified based on simple t-test. The criteria used to classify a gene as differentially expressed were that there was at least a two-fold change in the gene expression compared to the controls and a P-value less than 0.01 for the fold-change difference. By locus link ID, a total of 1814 genes from the comfrey treatment and 639 genes from the riddelliine treatment satisfied the requirements; both with about equal numbers of up- and down-regulated genes in response to the treatment [32, 33]. The number of genes whose expression was altered by comfrey was much higher than that by riddelliine, possibly due to the fact that comfrey is a complex mixture of biologically active compounds. By locus link ID, there were 302 genes that were common between comfrey and riddelliine; the majority of these common genes (90%, 273 genes) were regulated in the same direction (Figure 3). Since PAs generally produce the same types of DNA adducts [10] and both comfrey and riddelliine have similar mutational spectra [12, 32], it would be expected that the common genes altered by riddelliine and comfrey in this study may contribute to PA-induced toxicity.
Figure 3
Figure 3

Numbers of differentially expressed genes (DGEs) regulated by riddelliine and comfrey treatment. A gene was identified as differentially expressed if the fold-change was greater than 2 (up or down) and the P-value was less than 0.01 in comparison to the control group. The color in green refers to the number of genes whose expressions were significantly altered by both comfrey and riddelliine.

Common drug metabolizing genes induced by comfrey and riddelliine

A striking feature of the biotransformation enzymes is that their activities can be induced following exposure of chemicals or drugs. Therefore, examination of drug metabolizing genes (DMGs) can provide important pharmacological information, with potential clinical and toxicity implications for a given chemical. Liver is the major organ for biotransformation of xenobiotics and drugs, and PAs require metabolic activation to exert their biological and toxicological effects. Study and comparison of these genes altered by comfrey and riddelliine can provide us valuable information about common mechanisms on PAs' biotransformation and metabolisms. Also, in these steps PAs are activated into carcinogen, initiating the carcinogenesis. Therefore, we focused our investigation on the expression changes of DMGs. There were 45 and 87 DMGs up- or down-regulated by the treatment of riddelliine and comfrey, when a cut off of 2-fold change and P < 0.01 were used as the criteria to select genes. The expressions of 22 of these genes were altered by both comfrey and riddelliine, and are detailed in Table 1. The changed expression of each gene was always in the same direction for the 2 compounds. These commonly regulated DMGs were grouped into the three drug metabolism phases (phase I, II and III).
Table 1

Genes involved in drug metabolism altered by comfrey and riddelliine treatments in liver

Gene symbol

Gene description

Locus link ID

Fold changea

   

Comfrey

Riddelliine

Phase I metabolism

   

Ces2

carboxylesterase 2 (intestine, liver)

171118

10.4

2.4

Cyp2c

cytochrome P450, subfamily IIC

29277

-24.3

-270.3

Cyp2c40

cytochrome P450, family 2, subfamily c

25011

8.7

42.3

Cyp2e1

cytochrome P450, family 2, subfamily e

25086

2.2

2.1

Fmo5

flavin containing monooxygenase 5

246248

4.1

4.2

Phase II metabolism

   

Gsta3

glutathione S-transferase, alpha 3

14859

22.5

13.4

Inmt

indolethylamine N-methyltransferase

21743

-29.3

-11.7

Nqo1

NAD(P)H dehydrogenase, quinone 1

24314

5.0

2.6

Sult1c1

sulfotransferase family, cytosolic,1C, member1

65185

-2.7

-3.2

Phase III metabolism

   

Abcb1

ATP-binding cassette, sub-family B

24646

97.0

11.3

Abcb9

ATP-binding cassette, sub-family B, member 9

63886

-5.0

-2.7

Abcc3

ATP-binding cassette, sub-family C, member 3

140668

25.1

10.8

Abcc8

ATP-binding cassette, sub-family C, member 8

25559

-9.7

-2.6

Atp13a5

ATPase type 13A5

268878

-3.6

-4.1

Slc13a5

solute carrier family 13, member 5

266998

-2.9

-2.4

Slc16a4

solute carrier family 16, member 4

229699

-2.8

-4.9

Slc22a6

solute carrier family 22, member 6

29509

-3.7

-2.4

Slc22a8

solute carrier family 22, member 8

83500

-12.2

-833.3

Slc25a21

solute carrier family 25, member 21

171151

4.4

2.7

Slc25a30

solute carrier family 25, member 30

67554

-4.5

-2.7

aThe symbol of minus (-) means down-regulation, and the P-value for all genes is less than 0.01.

There are postulated to be three phase 1 pathways for the metabolism of PA; (1) Oxidation via cytochrome P450 (Cyp), (2) N-oxidation via flavin-containing monooxygenases (Fmo), and (3) hydrolysis via carboxylesterases. Hepatic cytochrome P450, Cyp3a and Cyp2b isoforms, were the major metabolizing enzymes involved in PA metabolism in humans [10]. One gene of the Cyp2b family (Cyp2b15) and three genes of the Cyp3a family (Cyp3a2, 3a9 and 3a18) were contained on the microarray. Although riddelliine elevated the level of Cyp3a9 RNA, no significant changes in the expression of these genes was induced by both riddelliine and comfrey exposure in our study. However, a significant induction of Cyp2c12 (Cyp2c40) was found with the expression increased by 8.7- and 42-fold by the treatments with comfrey and riddelliine, respectively (Table 1). This suggests that Cyp2c may play an important role in metabolizing comfrey and riddelliine in addition to Cyp3a and 2b.

Flavin-containing monooxygenases were reported to be involved in the biotransformation of PAs to the N-oxide metabolites [34], and N-oxidation appears to be one of the crucial pathways in determining the toxicity of PAs. However, it is not well known which particular family member contributes to this metabolism and toxicity. Among five flavin-containing monooxygenase genes (Fmo1-5) contained on microarray, Fmo5 was induced about 4-fold by both comfrey and riddelliine treatments. Hydrolysis of PAs is considered to be a major detoxification pathway in phase I. As shown in Table 1, carboxylesterase 2 (Ces2) was induced by the treatments of comfrey and riddelliine in our study. Our results not only confirm the findings reported in the literature but also provide additional detail that Fmo5 may be the form of flavin-containing monooxygenases family involved in N-oxidation and carboxylesterase 2 may play a role in hydrolysis of these chemicals.

Phase II and phase III enzymes usually function as detoxication. Phase II enzymes involve in conjugation of the polar functional groups of phase I metabolites and phase III enzymes are transporters for these metabolites. In general, exposure to some xenobiotics can trigger cellular "stress" response leading to an increase in the gene expression of phases II and III, which ultimately enhances the elimination and clearance of these xenobiotics. Treatment with comfrey and riddelliine has been shown to alter the expression of a number of phases II and III genes along with phase I genes (Table 1). The expression of the glutathione S-transferase, alpha 3 (Gsta3) gene increased 22- and 13-fold over the control with comfrey and riddelliine treatment, respectively. Glutathione S-transferases are major phase II detoxification enzymes found mainly in the cytosol. It has been reported that PAs can be conjugated with glutathione by these enzymes, reduce glutathione concentrations in liver, and increase the activity of these conjugation enzymes after PAs treatment [35–37]. This result is consistent with those previously studied. In the phase III genes, ATP-binding cassette, sub-family B (Abcb1) was induced by comfrey and riddelliine 97- and 11-fold over the control, respectively. No publications have been found about the relationship between PA treatment and induction of ATP-binding cassette. Future research is encouraged for the role of ATP-binding cassette in PAs' metabolism and carcinogenesis.

Not surprisingly, there were more DMGs altered by comfrey treatment than those by riddelliine treatment. This may be because comfrey is a mixture of many substances in addition to PAs, and these other substances may be involved in the induction of other DMGs. Nearly 50% of DMGs altered by the treatment with riddelliine were also found altered by the treatment with comfrey. It should be note that all of these 22 DMGs regulated by both chemicals changed in the same direction (Table 1). The gene expression similarity between the two treatments was further assessed by calculating the correlation coefficients of the log2 fold changes. The result is displayed as a scatter plot in Figure 4. The correlation coefficient was 0.72, indicating a fairly good agreement between the two treatments and suggesting that the toxicity of comfrey results from activation of PA. Although the process of drug metabolism is complicated, our results at least indicate that there is a common mechanism of drug metabolism involved in these two chemical treatments.
Figure 4
Figure 4

Comparison of fold-change of drug metabolizing genes (DMG) whose expression was altered by both comfrey and riddelliine treatments. Twenty-two DMGs detailed in Table 1 were commonly regulated by comfrey and riddelliine.

Common biological processes and genes associated with carcinogenesis

The Ingenuity Pathway Analysis Knowledge database, which provides a classification of gene products into molecular functions, biological processes, and cellular components, was used to help understand the biological consequences of exposure to comfrey or riddelliine. Biological process was examined for the genes from comfrey and riddelliine individually, and the software placed the processes into functional processes or categories such as cancer and cell death. After a cut off of P < 0.01, we observed that there were 83 and 118 functional processes for riddelliine and comfrey treatments, respectively. Among them, 46 of the function processes were altered by both comfrey and riddelliine, which was about half of the function processes altered by riddelliine and one-third of the function processes altered by comfrey (Figure 5). The top categories for these functions are listed in Table 2, including cancer, cell death, cell morphology, cell-to-cell signaling and interaction, and tissue development. These results suggest that common mechanism(s) may be responsible for the toxicity of both comfrey and riddelliine.
Table 2

The major relevant functions altered by both comfrey and riddelliine treatments in liver

Function category

Number of processes

Cancer

35

Cell death

35

Small molecule biochemistry

17

Cell morphology

14

Cell-to-cell signaling and interaction

14

Lipid metabolism

12

Molecular transport

12

Tissue development

11

Cellular movement

10

Figure 5
Figure 5

Number of regulated function processes significantly altered ( P < 0.01) by comfrey and riddelliine treatments. The color in green refers to the number of function processes that were significantly altered by both riddelliine and comfrey.

The finding that both comfrey and riddelliine treatments trigger gene expression alterations associated with cancer and cell death processes in this study is likely due to the PAs. As discussed above, PAs produce hepatotoxicity and carcinogenesis in experimental animals and humans. Functional analysis with Ingenuity Pathway Analysis showed that comfrey which contains 9 PAs [7–9] altered many more genes involved in cancer related pathways than riddelliine. The expressions of 387 cancer-related genes were significantly altered after comfrey exposure and the expressions of 84 cancer-related genes were significantly changed by riddelliine treatment. These genes are involved in the many functions associated with different stages of carcinogenesis. A large number of genes related to both treatments may indicate the PAs' carcinogenic insults. Table 3 shows the detailed information of 42 genes whose expressions were significantly altered by both riddelliine and comfrey treatments. These common cancer-related genes include genes involved in apoptosis and cell death, invasion, cell growth, cell morphology, and cell cycle. Figure 6 shows the strong correlation between the log2 fold-changes in gene expression caused by the two treatments; only two of the genes were regulated in the opposite directions. These results also indicate that the carcinogenicity of comfrey is generated from PAs.
Table 3

Genes involved in carcinogenesis altered by comfrey and riddelliine treatments in liver

Gene symbol

Gene description

Locus link ID

Fold changea

   

Comfrey

Riddelliine

AKR1C3

aldo-keto reductase family 1, member C3

171516

-4.9

-39.5

AR

androgen receptor (dihydrotestosterone receptor)

24208

-13.6

-19.6

ARHGDIG

Rho GDP dissociation inhibitor (GDI) gamma

14570

-8.9

-12.1

C5ORF13

chromosome 5 open reading frame 13

338475

-7.5

-7.6

CAST

Calpastatin (includes EG:831)

25403

-3.3

-5.4

CD36

CD36 molecule (thrombospondin receptor)

29184

3.3

7.9

CDKN1C

cyclin-dependent kinase inhibitor 1C (p57, Kip2)

246060

5.1

2.7

CHEK1

CHK1 checkpoint homolog (S. pombe)

140583

2.1

2.5

CREM

cAMP responsive element modulator

25620

-2.6

-3.3

DOK1

docking protein 1, 62kDa (downstream of tyrosine kinase 1)

13448

4.5

2.9

ECGF1

endothelial cell growth factor 1 (platelet-derived)

72962

-6.7

-2.8

EGR1

early growth response 1

24330

12.2

2.2

FABP2

fatty acid binding protein 2

25598

3.3

2.7

FAS

Fas (TNF receptor superfamily, member 6)

246097

3.3

3.2

FEZ1

fasciculation and elongation protein zeta 1 (zygin I)

81730

-3.2

2.9

FUT1

fucosyltransferase 1

81919

-4.8

-6.9

FYN

FYN oncogene related to SRC, FGR, YES

25150

4.3

2.8

GDF15

growth differentiation factor 15

29455

4.8

2.3

HGF

hepatocyte growth factor (hepapoietin A; scatter factor)

24446

3.6

2.4

HRASLS3

HRAS-like suppressor 3

24913

3.0

4.0

HSPA1A

heat shock 70kDa protein 1A

3303

4.1

8.1

ID1

inhibitor of DNA binding 1,

25261

3.9

2.2

IGFBP2

insulin-like growth factor binding protein 2, 36kDa

25662

3.2

15.7

ITGA4

integrin, alpha 4 (antigen CD49D)

16401

5.4

2.1

LAMA5

laminin, alpha 5

3911

37.5

2.1

LIMK2

LIM domain kinase 2

29524

4.1

2.0

MGAT5

mannoside acetylglucosaminyltransferase 5

65271

-2.2

2.1

MT1A

metallothionein 1A (functional)

24567

5.8

12.2

NCR1

natural cytotoxicity triggering receptor 1

117547

2.5

61.7

NQO1

NAD(P)H dehydrogenase, quinone 1

24314

5.0

2.6

PPARG

peroxisome proliferative activated receptor, gamma

25664

-2.6

-5.3

PRLR

prolactin receptor

24684

29.8

284.4

PTK6

PTK6 protein tyrosine kinase 6

20459

-16.7

-3.7

RTN4

reticulon 4

83765

6.2

3.4

SERPINA5

serpin peptidase inhibitor, clade A, member 5

65051

-3.1

-2.2

SMOX

spermine oxidase

228608

2.0

2.3

ST6GAL1

ST6 beta-galactosamide alpha-2,6-sialyltranferase 1

20440

-2.7

-5.7

ST8SIA1

ST8 alpha-N-acetyl-neuraminide alpha-2,8-sialyltransferase 1

25280

-6.8

-6.8

TIMP3

TIMP metallopeptidase inhibitor 3

25358

2.1

2.1

TNF

tumor necrosis factor (TNF superfamily, member 2)

24835

-4.2

-24.5

TNFSF10

tumor necrosis factor (ligand) superfamily, member 10

246775

2.4

2.6

TUBG1

tubulin, gamma 1

252921

-3.2

-4.2

aThe symbol of minus (-) means down-regulation, and the P-value for all genes is less than 0.01.

Figure 6
Figure 6

Comparison of fold-change of cancer-related genes whose expression was altered by both comfrey and riddelliine treatments. The genes commonly regulated by comfrey and riddelliine are listed in Table 3. Black symbol means regulated in same direction, and red symbol means regulated in opposite direction.

Most of genes involved in apoptosis (Cdkn1c, Egr1, Fas, Gdf15, Hgf, Hrasls3, Hspa1a, Lama5, Smox, Tnfsf10) were up-regulated and may be explained by the removal of cells damaged by PAs. For example, Egr1 (early growth response 1), whose expression increased 12- and 2-fold over the control with comfrey and riddelliine treatments, respectively, has the ability to function in numerous capacities, including differentiation, growth, growth inhibition, and apoptosis depending on the cell type and the stimulus. In response to stress, Egr1 displays a remarkable functional similarity to p53 and p73 [38] and may be involved in protection of the cell against the toxicity of these compounds. We also observed that comfrey and riddelliine significantly up-regulated Hspa1a, the gene encoding heat shock protein (HSP) 70 family members. High levels of inducible HSP70s prevent stress-induced apoptosis and block caspase activity, mitochondrial damage, and nuclear fragmentation [39]. These genes may all be involved in repair and protection of the cell against the damage caused by the PAs.

Among the common cancer-related genes, the expressions of two genes (Fez1 and Mgat5) were altered in opposite directions (down-regulated by comfrey and up-regulated by riddelliine). Mgat5 encodes N-acetylglucosaminyltransferase V (GlcNAc-TV), the Golgi enzyme required in the N-glycan processing pathway that modifies glycoproteins, including the cytokine receptors. Mgat5 gene expression is up-regulated by Ras pathway activation [40]. In Mgat5-/- mice, tumor latency is longer and metastasis is reduced compared with Mgat5+/+ mice [41]. FEZ1 (Fasciculation and elongation protein zeta-1) is a tumor suppressor gene that maps to chromosome 8p22, a chromosomal region frequently deleted in many human malignancies. The FEZ1 gene is expressed almost ubiquitously in normal tissues, and is prone to inactivation in human cancer [42]. In addition, introduction of FEZ1 into Fez1-negative cancer cells results in suppression of tumorigenicity and reduced cell growth with accumulation of cells at late S-G2/M stage of the cell cycle [43]. These two genes appear to be important in the carcinogenesis process and the differences in their expression in the two treatment groups may be an indication that the livers in the two groups are in different stages of the carcinogenic process. Alternatively, this may reflect the contribution of other components in comfrey.

Conclusion

Previously we found that the mutagenicity of comfrey was associated with the PAs contained in the plant. Therefore, we suggested that the toxicity and carcinogenicity of comfrey resulted from these PAs. In this study, we analyzed gene expression profiles in the liver of rats treated with riddelliine, a widely studied genotoxic PA and a proven rodent mutagen and carcinogen, and comfrey, an herbal tea that contains as many as nine PAs. We used the systematic approaches of HCA, PCA, functional analysis, and biological process analysis. Although there were large differences between the genes whose expression was significantly altered by these two treatments, there are common genes and function processes which may be indicative of effects due to PAs. In support of this is the finding that there were strong correlations between the gene expression fold-change alterations caused by the two treatments, in particular, drug metabolizing genes and cancer-related genes. These results suggest that comfrey and riddelliine may share common mechanism(s) of carcinogenesis, and further support our hypothesis that PAs are the main active components for carcinogenesis of comfrey.

Materials and methods

Treatments with comfrey and riddelliine

The details on the description of the in vivo portion of this study has been described previously [26, 32]. Briefly, Big Blue Fisher 344 transgenic rats were fed with a diet containing 8% comfrey root for 12 weeks or gavaged with 1 mg/kg riddelliine 5 times a week for 12 weeks. One day after the last treatment, 6 rats from each vehicle control, comfrey-treated, and riddelliine-treated groups were sacrificed, and the livers were isolated and stored at -80°C for mutagenesis and microarray analysis. Big Blue 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.

Microarray processing

The detailed description of the microarray processing of this study has been reported previously [32, 33]. Briefly, total RNA was isolated from liver tissues of 6 control, 6 comfrey-fed, and 6 riddelliine-treated rats using an RNeasy system (Qiagen, Chatsworth, CA). All RNA targets were labeled using the Applied Biosystems RT-IVT Labeling Kit Version 2.0. Digoxigenin labeled cRNA targets were hybridized to Applied Biosystems Rat Whole Genome Survey Microarrays using the Applied Biosystems Chemiluminescent Detection Kit. 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

Raw microarray intensity data from the Applied Biosystems' Rat Genome Survey Microarray, which is a one channel microarray with chemiluminescence detection and contains 26,857 probes, were normalized with quantile normalization which is recommended by the manufacturer. The normalized data were then input to ArrayTrack, a software system developed by the NCTR for the management, analysis, visualization and interpretation of microarray data [44]. The identification of differentially expressed genes based on fold-change and t-tests cutoffs, and HCA and PCA were conducted within ArrayTrack. Ingenuity Pathway Analysis (Mountain View, CA) was used for network and function analysis.

Declarations

Acknowledgements

The views presented in this article do not necessarily reflect those of the U.S. Food and Drug Administration.

This article has been published as part of BMC Bioinformatics Volume 8 Supplement 7, 2007: Proceedings of the Fourth Annual MCBIOS Conference. Computational Frontiers in Biomedicine. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/8?issue=S7.

Authors’ Affiliations

(1)
Division of Systems Toxicology, National Center for Toxicological Research, FDA, Jefferson, AR 72079, USA
(2)
Division of Genetic and Reproductive Toxicology, National Center for Toxicological Research, FDA, Jefferson, AR 72079, USA

References

  1. Rode D: Comfrey toxicity revisited. Trends Pharmacol Sci 2002,23(11):497–499. 10.1016/S0165-6147(02)02106-5View ArticlePubMedGoogle Scholar
  2. Yeong ML, Swinburn B, Kennedy M, Nicholson G: Hepatic veno-occlusive disease associated with comfrey ingestion. J Gastroenterol Hepatol 1990,5(2):211–214.View ArticlePubMedGoogle Scholar
  3. 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.View ArticleGoogle Scholar
  4. 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.PubMedGoogle Scholar
  5. Hirono I, Mori H, Haga M: Carcinogenic activity of Symphytum officinale. J Natl Cancer Inst 1978,61(3):865–869.PubMedGoogle Scholar
  6. USFDA: Center for Food Safety and Applied Nutrition: FDA advises dietary supplement manufacturers to remove comfrey products from the market. FDA 2001. [http://vmcfsanfdagov/~dms/dspltr06html]Google Scholar
  7. 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/np0004653View ArticlePubMedGoogle Scholar
  8. Schaneberg BT, Molyneux RJ, Khan IA: Evaporative light scattering detection of pyrrolizidine alkaloids. Phytochem Anal 2004,15(1):36–39. 10.1002/pca.715View ArticlePubMedGoogle Scholar
  9. Stickel F, Seitz HK: The efficacy and safety of comfrey. Public Health Nutr 2000,3(4A):501–508.View ArticlePubMedGoogle Scholar
  10. Fu PP, Xia Q, Lin G, Chou MW: Pyrrolizidine alkaloids – genotoxicity, metabolism enzymes, metabolic activation, and mechanisms. Drug Metab Rev 2004,36(1):1–55. 10.1081/DMR-120028426View ArticlePubMedGoogle Scholar
  11. Prakash AS, Pereira TN, Reilly PE, Seawright AA: Pyrrolizidine alkaloids in human diet. Mutat Res 1999,443(1–2):53–67.View ArticlePubMedGoogle Scholar
  12. 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.6602420PubMed CentralView ArticlePubMedGoogle Scholar
  13. 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.View ArticlePubMedGoogle Scholar
  14. Knudson AG: Antioncogenes and human cancer. Proc Natl Acad Sci USA 1993,90(23):10914–10921. 10.1073/pnas.90.23.10914PubMed CentralView ArticlePubMedGoogle Scholar
  15. Moolgavkar SH, Knudson AG Jr: Mutation and cancer: a model for human carcinogenesis. J Natl Cancer Inst 1981,66(6):1037–1052.PubMedGoogle Scholar
  16. Loeb LA: A mutator phenotype in cancer. Cancer Res 2001,61(8):3230–3239.PubMedGoogle Scholar
  17. Greenman C, Stephens P, Smith R, Dalgliesh GL, Hunter C, Bignell G, Davies H, Teague J, Butler A, Stevens C, et al.: Patterns of somatic mutation in human cancer genomes. Nature 2007,446(7132):153–158. 10.1038/nature05610PubMed CentralView ArticlePubMedGoogle Scholar
  18. Chan PC, Haseman JK, Prejean JD, Nyska A: Toxicity and carcinogenicity of riddelliine in rats and mice. Toxicol Lett 2003,144(3):295–311. 10.1016/S0378-4274(03)00240-6View ArticlePubMedGoogle Scholar
  19. Chan PC, Mahler J, Bucher JR, Travlos GS, Reid JB: Toxicity and carcinogenicity of riddelliine following 13 weeks of treatment to rats and mice. Toxicon 1994,32(8):891–908. 10.1016/0041-0101(94)90368-9View ArticlePubMedGoogle Scholar
  20. National Toxicology Program: Toxicology and carcinogenesis studies of riddelliine (CAS No. 23246–96–0) in F344/N rats and B6C3F1 mice (gavage studies). Natl Toxicol Program Tech Rep Ser 2003, (508):1–280.Google Scholar
  21. Yang Y, Yan J, Churchwell M, Beger R, Chan P, Doerge DR, Fu PP, Chou MW: Development of a (32)P-postlabeling/HPLC method for detection of dehydroretronecine-derived DNA adducts in vivo and in vitro. Chem Res Toxicol 2001,14(1):91–100. 10.1021/tx000149oView ArticlePubMedGoogle Scholar
  22. Yang YC, Yan J, Doerge DR, Chan PC, Fu PP, Chou MW: Metabolic activation of the tumorigenic pyrrolizidine alkaloid, riddelliine, leading to DNA adduct formation in vivo. Chem Res Toxicol 2001,14(1):101–109. 10.1021/tx000150nView ArticlePubMedGoogle Scholar
  23. Chou MW, Jian Y, Williams LD, Xia Q, Churchwell M, Doerge DR, Fu PP: Identification of DNA adducts derived from riddelliine, a carcinogenic pyrrolizidine alkaloid. Chem Res Toxicol 2003,16(9):1130–1137. 10.1021/tx030018yView ArticlePubMedGoogle Scholar
  24. Chou MW, Wang YP, Yan J, Yang YC, Beger RD, Williams LD, Doerge DR, Fu PP: Riddelliine N-oxide is a phytochemical and mammalian metabolite with genotoxic activity that is comparable to the parent pyrrolizidine alkaloid riddelliine. Toxicol Lett 2003,145(3):239–247. 10.1016/S0378-4274(03)00293-5View ArticlePubMedGoogle Scholar
  25. Chou MW, Yan J, Nichols J, Xia Q, Beland FA, Chan PC, Fu PP: Correlation of DNA adduct formation and riddelliine-induced liver tumorigenesis in F344 rats and B6C3F(1) mice. Cancer Lett 2003,193(2):119–125. 10.1016/S0304-3835(03)00045-4View ArticlePubMedGoogle Scholar
  26. 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/tx049955bView ArticlePubMedGoogle Scholar
  27. 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.013View ArticlePubMedGoogle Scholar
  28. Lettieri T: Recent applications of DNA microarray technology to toxicology and ecotoxicology. Environ Health Perspect 2006,114(1):4–9.PubMed CentralPubMedGoogle Scholar
  29. Liu-Stratton Y, Roy S, Sen CK: DNA microarray technology in nutraceutical and food safety. Toxicol Lett 2004,150(1):29–42. 10.1016/j.toxlet.2003.08.009View ArticlePubMedGoogle Scholar
  30. Boorman GA, Blackshear PE, Parker JS, Lobenhofer EK, Malarkey DE, Vallant MK, Gerken DK, Irwin RD: Hepatic gene expression changes throughout the day in the Fischer rat: implications for toxicogenomic experiments. Toxicol Sci 2005,86(1):185–193. 10.1093/toxsci/kfi166View ArticlePubMedGoogle Scholar
  31. Desai VG, Moland CL, Branham WS, Delongchamp RR, Fang H, Duffy PH, Peterson CA, Beggs ML, Fuscoe JC: Changes in expression level of genes as a function of time of day in the liver of rats. Mutat Res 2004,549(1–2):115–129.View ArticlePubMedGoogle Scholar
  32. Mei N, Guo L, Zhang L, Shi L, Sun YA, Fung C, Moland CL, Dial SL, Fuscoe JC, Chen T: 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 2006,7(Suppl 2):S16. 10.1186/1471-2105-7-S2-S16PubMed CentralView ArticlePubMedGoogle Scholar
  33. Mei N, Guo L, Liu R, Fuscoe JC, Chen T: Gene expression changes induced by the tumorigenic pyrrolizidine alkaloid riddelliine in liver of Big Blue rats. BMC Bioinformatics 2007, in press.Google Scholar
  34. Williams DE, Reed RL, Kedzierski B, Ziegler DM, Buhler DR: The role of flavin-containing monooxygenase in the N-oxidation of the pyrrolizidine alkaloid senecionine. Drug Metab Dispos 1989,17(4):380–386.PubMedGoogle Scholar
  35. Lame MW, Morin D, Jones AD, Segall HJ, Wilson DW: Isolation and identification of a pyrrolic glutathione conjugate metabolite of the pyrrolizidine alkaloid monocrotaline. Toxicol Lett 1990,51(3):321–329. 10.1016/0378-4274(90)90075-WView ArticlePubMedGoogle Scholar
  36. Lin G, Cui YY, Hawes EM: Microsomal formation of a pyrrolic alcohol glutathione conjugate of clivorine. Firm evidence for the formation of a pyrrolic metabolite of an otonecine-type pyrrolizidine alkaloid. Drug Metab Dispos 1998,26(2):181–184.PubMedGoogle Scholar
  37. Yan CC, Huxtable RJ: Effects of monocrotaline, a pyrrolizidine alkaloid, on glutathione metabolism in the rat. Biochem Pharmacol 1996,51(3):375–379. 10.1016/0006-2952(95)02189-2View ArticlePubMedGoogle Scholar
  38. Yu J, Baron V, Mercola D, Mustelin T, Adamson ED: A network of p73, p53 and Egr1 is required for efficient apoptosis in tumor cells. Cell Death Differ 2007,14(3):436–446. 10.1038/sj.cdd.4402029View ArticlePubMedGoogle Scholar
  39. Mosser DD, Caron AW, Bourget L, Denis-Larose C, Massie B: Role of the human heat shock protein hsp70 in protection against stress-induced apoptosis. Mol Cell Biol 1997,17(9):5317–5327.PubMed CentralView ArticlePubMedGoogle Scholar
  40. Dennis JW, Granovsky M, Warren CE: Glycoprotein glycosylation and cancer progression. Biochim Biophys Acta 1999,1473(1):21–34.View ArticlePubMedGoogle Scholar
  41. Granovsky M, Fata J, Pawling J, Muller WJ, Khokha R, Dennis JW: Suppression of tumor growth and metastasis in Mgat5-deficient mice. Nat Med 2000,6(3):306–312. 10.1038/73163View ArticlePubMedGoogle Scholar
  42. Ishii H, Baffa R, Numata SI, Murakumo Y, Rattan S, Inoue H, Mori M, Fidanza V, Alder H, Croce CM: The FEZ1 gene at chromosome 8p22 encodes a leucine-zipper protein, and its expression is altered in multiple human tumors. Proc Natl Acad Sci USA 1999,96(7):3928–3933. 10.1073/pnas.96.7.3928PubMed CentralView ArticlePubMedGoogle Scholar
  43. Ishii H, Vecchione A, Murakumo Y, Baldassarre G, Numata S, Trapasso F, Alder H, Baffa R, Croce CM: FEZ1/LZTS1 gene at 8p22 suppresses cancer cell growth and regulates mitosis. Proc Natl Acad Sci USA 2001,98(18):10374–10379. 10.1073/pnas.181222898PubMed CentralView ArticlePubMedGoogle Scholar
  44. Tong W, Cao X, Harris S, Sun H, Fang H, Fuscoe J, Harris A, Hong H, Xie Q, Perkins R, et al.: ArrayTrack – supporting toxicogenomic research at the U.S. Food and Drug Administration National Center for Toxicological Research. Environ Health Perspect 2003,111(15):1819–1826.PubMed CentralView ArticlePubMedGoogle Scholar

Copyright

© Guo et al; licensee BioMed Central Ltd. 2007

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Advertisement

BMC Bioinformatics

ISSN: 1471-2105

Contact us