Search

Preclinical and clinical findings suggest sexual dimorphism in cardiotoxicity induced by a chemotherapeutic drug, doxorubicin (DOX). However, molecular alterations leading to sex-related differential vulnerability of heart to DOX toxicity are not fully explored. In the present study, RNA sequencing in hearts of B6C3F

The rat is an important model organism in biomedical research for studying human disease mechanisms and treatments, but its annotated transcriptome is far from complete. We constructed a Rat Transcriptome Re-annotation named RTR using RNA-seq data from 320 samples in 11 different organs generated by the SEQC consortium. Totally, there are 52 807 genes and 114 152 transcripts in RTR. Transcribed regions and exons in RTR account for ∼42% and ∼6.5% of the genome, respectively. Of all 73 074 newly annotated transcripts in RTR, 34 213 were annotated as high confident coding transcripts and 24 728 as high confident long noncoding transcripts. Different tissues rather than different stages have a significant influence on the expression patterns of transcripts. We also found that 11 715 genes and 15 852 transcripts were expressed in all 11 tissues and that 849 house-keeping genes expressed different isoforms among tissues. This comprehensive transcriptome is freely available at http://www.unimd.org/rtr/. Our new rat transcriptome provides essential reference for genetics and gene expression studies in rat disease and toxicity models.

Cardiotoxicity is a serious adverse effect of an anticancer drug, doxorubicin (DOX), which can occur within a year or decades after completion of therapy. The present study was designed to address a knowledge gap concerning a lack of circulating biomarkers capable of predicting the risk of cardiotoxicity induced by DOX. Profiling of 2083 microRNAs (miRNAs) in mouse plasma revealed 81 differentially expressed miRNAs 1 week after 6, 9, 12, 18, or 24 mg/kg total cumulative DOX doses (early-onset model) or saline (SAL). Among these, the expression of seven miRNAs was altered prior to the onset of myocardial injury at 12 mg/kg and higher cumulative doses. The expression of only miR-34a-5p was significantly (false discovery rate [FDR] lt; 0.1) elevated at all total cumulative doses compared with concurrent SAL-treated controls and showed a statistically significant dose-related response. The trend in plasma miR-34a-5p expression levels during DOX exposures also correlated with a significant dose-related increase in cardiac expression of miR-34a-5p in these mice. Administration of a cardioprotective drug, dexrazoxane, to mice before DOX treatment, significantly mitigated miR-34a-5p expression in both plasma and heart in conjunction with attenuation of cardiac pathology. This association between plasma and heart may suggest miR-34a-5p as a potential early circulating marker of early-onset DOX cardiotoxicity. In addition, higher expression of miR-34a-5p (FDR lt; 0.1) in plasma and heart compared with SAL-treated controls 24 weeks after 24 mg/kg total cumulative DOX dose, when cardiac function was altered in our recently established delayed-onset cardiotoxicity model, indicated its potential as an early biomarker of delayed-onset cardiotoxicity.

Subclinical cardiotoxicity at low total cumulative doxorubicin (DOX) doses can manifest into cardiomyopathy in long-term cancer survivors. However, the underlying mechanisms are poorly understood. In male B6C3F1 mice, assessment of cardiac function by echocardiography was performed at 1, 4, 10, 17, and 24 weeks after exposure to 6, 9, 12, and 24 mg/kg total cumulative DOX doses or saline (SAL) to monitor development of delayed-onset cardiotoxicity. The 6- or 9-mg/kg total cumulative doses resulted in a significant time-dependent decline in systolic function (left ventricular ejection fraction (LVEF) and fractional shortening (FS)) during the 24-week recovery although there was not a significant alteration in % LVEF or % FS at any specific time point during the recovery. A significant decline in systolic function was elicited by the cardiotoxic cumulative DOX dose (24 mg/kg) during the 4- to 24-week period after treatment compared to SAL-treated counterparts. At 24 weeks after DOX treatment, a significant dose-related decrease in the expression of genes and proteins involved in sarcoplasmic reticulum (SR) calcium homeostasis (Ryr2 and Serca2) was associated with a dose-related increase in the transcript level of Casp12 (SR-specific apoptosis) in hearts. These mice also showed enhanced apoptotic activity in hearts indicated by a significant dose-related elevation in the number of apoptotic cardiomyocytes compared to SAL-treated counterparts. These findings collectively suggest that a steady decline in SR calcium handling and apoptosis might be involved in the development of subclinical cardiotoxicity that can evolve into irreversible cardiomyopathy later in life.

Specific recurring chromosomal translocations and deletions are found in a variety of cancers. In hematopoietic malignancies, many of these chromosomal aberrations result from mistakes involving V(D)J recombination. V(D)J recombination is required for the formation of functional T-cell receptor genes in T-cells and antibody genes in B-cells. This is an inherently dangerous process, however, because double-strand breaks are introduced into the chromosomes. Molecular evidence indicates that failure of the fidelity of this process results in the activation of proto-oncogenes or the inactivation of tumor-suppressor genes. Here we describe sensitive, quantitative PCR assays for the measurement of such events in human lymphocytes. One assay measures the frequency of t(14;18) translocations that result in the dysfunctional regulation of the anti-apoptotic gene BCL-2. The other assay measures the frequency of a deletion caused by illegitimate V(D)J recombination in the X-linked HPRT gene. The findings and conclusions presented in this article are the author's and do not necessarily reflect those of the Food and Drug Administration.

2-hydroxy-4-methoxybenzophenone (HMB) is an ultraviolet light-absorbing compound that is used in sunscreens, cosmetics and plastics. HMB has been reported to have weak estrogenic activity by in vivo and in vitro studies, making it a chemical with potential reproductive concern. To explore if prenatal and lactational HMB exposure alters gene expression profiles of the developing reproductive organs, we performed microarray analysis using the prostate and testis of postnatal day (PND) 30 male Sprague-Dawley rats offspring exposed to 0, 3000, or 30,000 ppm of HMB from gestational day 6 through PND 21. Gene expression profiles of the prostate and testis were differentially affected by HMB dose with significant alterations observed at the 30,000 ppm HMB group. Tissue-specific gene expression was also identified. These genes, whose expression was altered by HMB exposure, may be considered as candidate biomarker(s) for testicular or prostatic toxicity; however, further studies are necessary to explore this potential.

The mouse embryonic stem cell test (mEST) is a promising in vitro assay for predicting developmental toxicity. In the current study, early differentiation of D3 mouse embryonic stem cells (mESCs) under osteoblast culture conditions and embryotoxicity of cadmium sulfate were examined. D3 mESCs were exposed to cadmium sulfate for 24, 48 or 72 h, and whole genome transcriptional profiles were determined. The results indicate a track of differentiation was identified as mESCs differentiate. Biological processes that were associated with differentiation related genes included embryonic development and, specifically, skeletal system development. Cadmium sulfate inhibited mESC differentiation at all three time points. Functional pathway analysis indicated biological pathways affected included those related to skeletal development, renal and reproductive function. In summary, our results suggest that transcriptional profiles are a sensitive indicator of early mESC differentiation. Transcriptomics may improve the predictivity of the mEST by suggesting possible modes of action for tested chemicals.

Safety assessments of new drug candidates are an important part of the drug development and approval process. Often, possible sex-associated susceptibilities are not adequately addressed, and better assessment tools are needed. We hypothesized that hepatic transcript profiles of drug metabolizing enzymes and transporters (DMETs) can be used to predict sex-associated differences in drug metabolism, and possible adverse events. Comprehensive hepatic transcript profiles were generated for F344 rats of both sexes at nine ages, from 2 weeks (pre-weaning) to 104 weeks (elderly). Large differences in the transcript profiles of 29 DMETs were found between adult males and females (8-52 weeks). Using the PharmaPendium database, 41 drugs were found to be metabolized by one or two cytochrome P450 (Cyp) enzymes encoded by sexually dimorphic mRNAs, and thus were candidates for evaluation of possible sexually dimorphic metabolism and/or toxicities. Suspension cultures of primary hepatocytes from three male and three female adult rats (10-13 weeks old) were used to evaluate the metabolism of 11 drugs predicted to have sexually dimorphic metabolism. The pharmacokinetics of the drug or its metabolite was analyzed by liquid chromatography/tandem mass spectrometry using multiple reaction monitoring. Of those drugs with adequate metabolism, the predicted significant sex-different metabolism was found for six of seven drugs, with half-lives 37%- 400% longer in female hepatocytes than in male hepatocytes. Thus, in this rat model, transcript profiles may allow identification of potential sex-related differences in drug metabolism. SIGNIFICANCE STATEMENT The present study showed that sex-different expression of genes coding for drug metabolizing enzymes, specifically cytochrome P450s, could be used to predict sex-different drug metabolism, and, thus, provide a new tool for protecting susceptible subpopulations from possible adverse drug events.

Sex is a risk factor for development of cardiotoxicity, induced by the anti-cancer drug, doxorubicin (DOX), in humans. To explore potential mechanisms underlying differential susceptibility to DOX between sexes, 8-week old male and female B6C3F1 mice were dosed with 3mg/kg body weight DOX or an equivalent volume of saline via tail vein once a week for 6, 7, 8, and 9 consecutive weeks, resulting in 18, 21, 24, and 27mg/kg cumulative DOX doses, respectively. At necropsy, one week after each consecutive final dose, the extent of myocardial injury was greater in male mice compared to females as indicated by higher plasma concentrations of cardiac troponin T at all cumulative DOX doses with statistically significant differences between sexes at the 21 and 24mg/kg cumulative doses. A greater susceptibility to DOX in male mice was further confirmed by the presence of cytoplasmic vacuolization in cardiomyocytes, with left atrium being more vulnerable to DOX cardiotoxicity. The number of TUNEL-positive cardiomyocytes was mostly higher in DOX-treated male mice compared to female counterparts, showing a statistically significant sex-related difference only in left atrium at 21mg/kg cumulative dose. DOX-treated male mice also had an increased number of γ-H2A.X-positive (measure of DNA double-strand breaks) cardiomyocytes compared to female counterparts with a significant sex effect in the ventricle at 27mg/kg cumulative dose and right atrium at 21 and 27mg/kg cumulative doses. This newly established mouse model provides a means to identify biomarkers and access potential mechanisms underlying sex-related differences in DOX-induced cardiotoxicity.

Cardiotoxicity is a serious adverse effect of doxorubicin (DOX) treatment in cancer patients. Currently, there is a lack of sensitive biomarkers to predict the risk of DOX-induced cardiotoxicity. Using SOMAmer-based proteomic technology, 1129 proteins were profiled to identify potential early biomarkers of cardiotoxicity in plasma from male B6C3F

Background Physiological factors such as age and sex have been shown to be risk factors for adverse effects in the liver, including liver diseases and drug-induced liver injury. Previously, we have reported age- and sex-related significant differences in hepatic basal gene expression in rats during the life span that may be related to susceptibility to such adverse effects. However, the underlying mechanisms of the gene expression changes were not fully understood. In recent years, increasing evidence for epigenetic mechanisms of gene regulation has fueled interest in the role of microRNAs (miRNAs) in toxicogenomics and biomarker discovery. We therefore proposed that significant age and sex differences exist in baseline liver miRNA expression, and that comprehensive profiling of miRNAs will provide insights into the epigenetic regulation of gene expression in rat liver. Methods To address this, liver tissues from male and female F344 rats were examined at 2, 5, 6, 8, 15, 21, 52, 78, and 104 weeks of age for the expression of 677 unique miRNAs. Following data processing, predictive pathway analysis was performed on selected miRNAs that exhibited prominent age and/or sex differences in expression. Results Of the 314 miRNAs found to be expressed, 214 were differentially expressed; 65 and 212 miRNAs showed significant (false discovery rate (FDR)

Furan is a heterocyclic organic compound produced in the chemical manufacturing industry and also found in a broad range of food products, including infant formulas and baby foods. Previous reports have indicated that the adverse biological effects of furan, including its liver tumorigenicity, may be associated with epigenetic abnormalities. In the present study, we investigated the persistence of epigenetic alterations in rat liver. Male F344 rats were treated by gavage 5 days per week with 8 mg furan/kg body weight (bw)/day for 90 days. After the last treatment, rats were divided randomly into 4 groups; 1 group of rats was sacrificed 24 h after the last treatment, whereas other groups were maintained without further furan treatment for an additional 90, 180, or 360 days. Treatment with furan for 90 days resulted in alterations in histone lysine methylation and acetylation, induction of base-excision DNA repair genes, suggesting oxidative damage to DNA, and changes in the gene expression in the livers. A majority of these furan-induced molecular changes was transient and disappeared after the cessation of furan treatment. In contrast, histone H3 lysine 9 and H3 lysine 56 showed a sustained and time-depended decrease in acetylation, which was associated with formation of heterochromatin and altered gene expression. These results indicate that furan-induced adverse effects may be mechanistically related to sustained changes in histone lysine acetylation that compromise the ability of cells to maintain and control properly the expression of genetic information.

The concordance of RNA-sequencing (RNA-seq) with microarrays for genome-wide analysis of differential gene expression has not been rigorously assessed using a range of chemical treatment conditions. Here we use a comprehensive study design to generate Illumina RNA-seq and Affymetrix microarray data from the same liver samples of rats exposed in triplicate to varying degrees of perturbation by 27 chemicals representing multiple modes of action (MOAs). The cross-platform concordance in terms of differentially expressed genes (DEGs) or enriched pathways is linearly correlated with treatment effect size (R-2 approximate to 0.8). Furthermore, the concordance is also affected by transcript abundance and biological complexity of the MOA. RNA-seq outperforms microarray (93% versus 75%) in DEG verification as assessed by quantitative PCR, with the gain mainly due to its improved accuracy for low-abundance transcripts. Nonetheless, classifiers to predict MOAs perform similarly when developed using data from either platform. Therefore, the endpoint studied and its biological complexity, transcript abundance and the genomic application are important factors in transcriptomic research and for clinical and regulatory decision making.

We present primary results from the Sequencing Quality Control (SEQC) project, coordinated by the US Food and Drug Administration. Examining Illumina HiSeq, Life Technologies SOLiD and Roche 454 platforms at multiple laboratory sites using reference RNA samples with built-in controls, we assess RNA sequencing (RNA-seq) performance for junction discovery and differential expression profiling and compare it to microarray and quantitative PCR (qPCR) data using complementary metrics. At all sequencing depths, we discover unannotated exon-exon junctions, with >80% validated by qPCR. We find that measurements of relative expression are accurate and reproducible across sites and platforms if specific-filters are used. In contrast, RNA-seq and microarrays do not provide accurate absolute measurements, and gene-specific biases are observed for all examined platforms, including qPCR. Measurement performance depends on the platform and data analysis pipeline, and variation is large for transcript-level profiling. The complete SEQC data sets, comprising >100 billion reads (10Tb), provide unique resources for evaluating RNA-seq analyses for clinical and regulatory settings.

The steady increase in the incidence and mortality of hepatocellular carcinoma (HCC) signifies a crucial need to understand better its pathogenesis to improve clinical management and prevention of the disease. The aim of this study was to investigate molecular mechanisms for the chemopreventive effects of folic acid and tributyrin alone or in combination on rat hepatocarcinogenesis. Male Wistar rats were subjected to a classic “resistant hepatocyte” model of liver carcinogenesis and treated with folic acid and tributyrin alone or in combination for 5 weeks during promotion stage. Treatment with folic acid and tributyrin alone or in combination strongly inhibited the development of glutathione-S-transferase placental form (GSTP)-positive foci. Microarray analysis showed significant changes in gene expression. A total of 501, 655, and 940 of differentially expressed genes, involved in cell cycle, p53-signaling, angiogenesis, and Wnt pathways, was identified in the livers of rats treated with folic acid, tributyrin or folic acid and tributyrin. A detailed analysis of these differentially expressed genes revealed that treatments inhibited angiogenesis in the preneoplastic livers. This was evidenced by the fact that 30 out of 77 differentially expressed genes common to all three treatments are involved in the regulation of the angiogenesis pathway. The inhibition of angiogenesis was confirmed by reduced levels of CD34 protein. In conclusion, the tumor-suppressing activity of folic acid and tributyrin is associated with inhibition of angiogenesis at early stages of rat liver carcinogenesis. Importantly, the combination of folic acid and tributyrin has stronger chemopreventive effect than each of the compounds alone.

Serum levels of cardiac troponins serve as biomarkers of myocardial injury. However, troponins are released into the serum only after damage to cardiac tissue has occurred. Here, we report development of a mouse model of doxorubicin (DOX)-induced chronic cardiotoxicity to aid in the identification of predictive biomarkers of early events of cardiac tissue injury. Male B6C3F(1) mice were administered intravenous DOX at 3mg/kg body weight, or an equivalent volume of saline, once a week for 4, 6, 8, 10, 12, and 14weeks, resulting in cumulative DOX doses of 12, 18, 24, 30, 36, and 42mg/kg, respectively. Mice were sacrificed a week following the last dose. A significant reduction in body weight gain was observed in mice following exposure to a weekly DOX dose for 1week and longer compared to saline-treated controls. DOX treatment also resulted in declines in red blood cell count, hemoglobin level, and hematocrit compared to saline-treated controls after the 2nd weekly dose until the 8th and 9th doses, followed by a modest recovery. All DOX-treated mice had significant elevations in cardiac troponin T concentrations in plasma compared to saline-treated controls, indicating cardiac tissue injury. Also, a dose-related increase in the severity of cardiac lesions was seen in mice exposed to 24mg/kg DOX and higher cumulative doses. Mice treated with cumulative DOX doses of 30mg/kg and higher showed a significant decline in heart rate, suggesting drug-induced cardiac dysfunction. Altogether, these findings demonstrate the development of DOX-induced chronic cardiotoxicity in B6C3F(1) mice.

Sex is a risk factor for development of cardiotoxicity, induced by the anti-cancer drug, doxorubicin (DOX), in humans. To explore potential mechanisms underlying differential susceptibility to DOX between sexes, 8-week old male and female B6C3F

Background Nonalcoholic fatty liver disease (NAFLD) is a major health problem and a leading cause of chronic liver disease in the United States and Western countries. In humans, genetic factors greatly influence individual susceptibility to NAFLD; nonetheless, the effect of inter-individual differences in the normal liver epigenome with regard to the susceptibility to NAFLD has not been determined. Results In the present study, we investigated the association between the DNA methylation status in the livers of A/J and WSB/EiJ mice and the severity of NAFLD-associated liver injury. We demonstrate that A/J and WSB/EiJ mice, which are characterized by significant differences in the severity of liver injury induced by a choline- and folate-deficient (CFD) diet exhibit substantial differences in cytosine DNA methylation in their normal livers. Furthermore, feeding A/J and WSB/EiJ mice a CFD diet for 12 weeks resulted in different trends and changes in hepatic cytosine DNA methylation. Conclusion Our findings indicate a primary role of hepatic DNA methylation in the pathogenesis of NAFLD and suggest that individual variations in DNA methylation across the genome may be a factor determining and influencing the vulnerability to NAFLD. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-2617-2) contains supplementary material, which is available to authorized users.

Nonalcoholic fatty liver disease (NAFLD) is a major health problem and a leading cause of chronic liver disease in the United States and developed countries. In humans, genetic factors greatly influence individual susceptibility to NAFLD. The goals of this study were to compare the magnitude of interindividual differences in the severity of liver injury induced by methyl-donor deficiency among individual inbred strains of mice and to investigate the underlying mechanisms associated with the variability. Feeding mice a choline- and folate-deficient diet for 12 wk caused liver injury similar to NAFLD. The magnitude of liver injury varied among the strains, with the order of sensitivity being A/J ≈ C57BL/6J ≈ C3H/HeJ < 129S1/SvImJ ≈ CAST/EiJ < PWK/PhJ < WSB/EiJ. The interstrain variability in severity of NAFLD liver damage was associated with dysregulation of genes involved in lipid metabolism, primarily with a down-regulation of the peroxisome proliferator receptor α (PPARα)-regulated lipid catabolic pathway genes. Markers of oxidative stress and oxidative stress-induced DNA damage were also elevated in the livers but were not correlated with severity of liver damage. These findings suggest that the PPARα-regulated metabolism network is one of the key mechanisms determining interstrain susceptibility and severity of NAFLD in mice.—Tryndyak, V., de Conti, A., Kobets, T., Kutanzi, K., Koturbash, I., Han, T., Fuscoe, J. C., Latendresse, J. R., Melnyk, S., Shymonyak, S., Collins, L., Ross, S. A., Rusyn, I., Beland, F. A., Pogribny, I. P. Interstrain differences in the severity of liver injury induced by a choline- and folate-deficient diet in mice are associated with dysregulation of genes involved in lipid metabolism.

An interagency collaboration was established to model chemical interactions that may cause adverse health effects when an exposure to a mixture of chemicals occurs. Many of these chemicals--drugs, pesticides, and environmental pollutants--interact at the level of metabolic biotransformations mediated by cytochrome P450 (CYP) enzymes. In the present work, spectral data-activity relationship (SDAR) and structure-activity relationship (SAR) approaches were used to develop machine-learning classifiers of inhibitors and non-inhibitors of the CYP3A4 and CYP2D6 isozymes. The models were built upon 602 reference pharmaceutical compounds whose interactions have been deduced from clinical data, and 100 additional chemicals that were used to evaluate model performance in an external validation (EV) test. SDAR is an innovative modeling approach that relies on discriminant analysis applied to binned nuclear magnetic resonance (NMR) spectral descriptors. In the present work, both 1D ¹³C and 1D ¹⁵N-NMR spectra were used together in a novel implementation of the SDAR technique. It was found that increasing the binning size of 1D ¹³C-NMR and ¹⁵N-NMR spectra caused an increase in the tenfold cross-validation (CV) performance in terms of both the rate of correct classification and sensitivity. The results of SDAR modeling were verified using SAR. For SAR modeling, a decision forest approach involving from 6 to 17 Mold2 descriptors in a tree was used. Average rates of correct classification of SDAR and SAR models in a hundred CV tests were 60% and 61% for CYP3A4, and 62% and 70% for CYP2D6, respectively. The rates of correct classification of SDAR and SAR models in the EV test were 73% and 86% for CYP3A4, and 76% and 90% for CYP2D6, respectively. Thus, both SDAR and SAR methods demonstrated a comparable performance in modeling a large set of structurally diverse data. Based on unique NMR structural descriptors, the new SDAR modeling method complements the existing SAR techniques, providing an independent estimator that can increase confidence in a structure-activity assessment. When modeling was applied to hazardous environmental chemicals, it was found that up to 20% of them may be substrates and up to 10% of them may be inhibitors of the CYP3A4 and CYP2D6 isoforms. The developed models provide a rare opportunity for the environmental health branch of the public health service to extrapolate to hazardous chemicals directly from human clinical data. Therefore, the pharmacological and environmental health branches are both expected to benefit from these reported models.

Polypharmacy increasingly has become a topic of public health concern, particularly as the U.S. population ages. Drug labels often contain insufficient information to enable the clinician to safely use multiple drugs. Because many of the drugs are bio-transformed by cytochrome P450 (CYP) enzymes, inhibition of CYP activity has long been associated with potentially adverse health effects. In an attempt to reduce the uncertainty pertaining to CYP-mediated drug-drug/chemical interactions, an interagency collaborative group developed a consensus approach to prioritizing information concerning CYP inhibition. The consensus involved computational molecular docking, spectral data-activity relationship (SDAR), and structure-activity relationship (SAR) models that addressed the clinical potency of CYP inhibition. The models were built upon chemicals that were categorized as either potent or weak inhibitors of the CYP3A4 isozyme. The categorization was carried out using information from clinical trials because currently available in vitro high-throughput screening data were not fully representative of the in vivo potency of inhibition. During categorization it was found that compounds, which break the Lipinski rule of five by molecular weight, were about twice more likely to be inhibitors of CYP3A4 compared to those, which obey the rule. Similarly, among inhibitors that break the rule, potent inhibitors were 2–3 times more frequent. The molecular docking classification relied on logistic regression, by which the docking scores from different docking algorithms, CYP3A4 three-dimensional structures, and binding sites on them were combined in a unified probabilistic model. The SDAR models employed a multiple linear regression approach applied to binned 1D 13C-NMR and 1D 15N-NMR spectral descriptors. Structure-based and physical-chemical descriptors were used as the basis for developing SAR models by the decision forest method. Thirty-three potent inhibitors and 88 weak inhibitors of CYP3A4 were used to train the models. Using these models, a synthetic majority rules consensus classifier was implemented, while the confidence of estimation was assigned following the percent agreement strategy. The classifier was applied to a testing set of 120 inhibitors not included in the development of the models. Five compounds of the test set, including known strong inhibitors dalfopristin and tioconazole, were classified as probable potent inhibitors of CYP3A4. Other known strong inhibitors, such as lopinavir, oltipraz, quercetin, raloxifene, and troglitazone, were among 18 compounds classified as plausible potent inhibitors of CYP3A4. The consensus estimation of inhibition potency is expected to aid in the nomination of pharmaceuticals, dietary supplements, environmental pollutants, and occupational and other chemicals for in-depth evaluation of the CYP3A4 inhibitory activity. It may serve also as an estimate of chemical interactions via CYP3A4 metabolic pharmacokinetic pathways occurring through polypharmacy and nutritional and environmental exposures to chemical mixtures.

Drug-induced liver injury (DILI) is the primary adverse event that results in the withdrawal of drugs from the market and a frequent reason for the failure of drug candidates in the pre-clinical or clinical phases of drug development. This paper presents an approach for identifying potential liver toxicity genomic biomarkers from a liver toxicity biomarker study involving the paired compounds entacapone (“non-liver toxic drug”) and tolcapone (“hepatotoxic drug”). Molecular analysis of the rat liver and plasma samples, combined with statistical analysis, revealed many similarities and differences between the in vivo biochemical effects of the two drugs. Six hundred and ninety-five genes and 61 pathways were selected based on the classification scheme. Of the 61 pathways, 5 were specific to treatment with tolcapone. Two of the 12 animals in the tolcapone group were found to have high ALT, AST, or TBIL levels. The gene Vars2 (valyl-tRNA synthetase 2) was identified in both animals and the pathway to which it belongs, the aminoacyl-tRNA biosynthesis pathway, was one of the three most significant tolcapone-specific pathways identified.

MicroRNAs (MiRNAs) are small noncoding RNAs that function as regulators of gene expression to control cell growth and differentiation. In this study, we analyzed miRNA and mRNA expression in the livers of rats treated with a carcinogenic dose of comfrey (Symphytum officinale) for 12 weeks. Groups of six rats were fed a normal diet or a diet containing 8% comfrey root. The animals were sacrificed 1 day after the last treatment and the livers were isolated for miRNA expression analysis using LC Sciences miRNA microarrays and for mRNA expression analysis using Affymetrix rat genome microarrays. MiRNA expression levels were significantly changed by comfrey treatment. The treated samples were separated clearly from the control samples in both principal component analysis (PCA) and hierarchical clustering analysis (HCA). Quantitative measurements of seven miRNAs using TaqMan real-time PCR were consistent with the microarray results in terms of fold-change and the direction of the change in expression. Forty-five miRNAs (P < 0.01) and 1,921 mRNAs (q = 0) were significantly changed by comfrey treatment. Using a target prediction algorithm, 434 differentially expressed genes (DEGs) were predicted to be targeted by the differentially expressed miRNAs (DEMs). The DEM-targeted DEGs were more likely to be involved in carcinogenesis than the DEGs that were not targeted by the DEMs. The nontargeted DEGs were enriched in noncancer-related biological processes. Our data suggest that comfrey may exert its carcinogenic effects by disturbing miRNA expression resulting in altered mRNA levels of the DEM-targeted genes that are functionally associated with carcinogenesis. Environ. Mol. Mutagen., 2011. © Published 2011 Wiley-Liss, Inc.

Comfrey has been consumed by humans as a vegetable and a tea and used as an herbal medicine for more than 2000 years. Comfrey, however, produces hepatotoxicity in livestock and humans and carcinogenicity in experimental animals. Comfrey contains as many as 14 pyrrolizidine alkaloids (PA), including 7-acetylintermedine, 7-acetyllycopsamine, echimidine, intermedine, lasiocarpine, lycopsamine, myoscorpine, symlandine, symphytine, and symviridine. The mechanisms underlying comfrey-induced genotoxicity and carcinogenicity are still not fully understood. The available evidence suggests that the active metabolites of PA in comfrey interact with DNA in liver endothelial cells and hepatocytes, resulting in DNA damage, mutation induction, and cancer development. Genotoxicities attributed to comfrey and riddelliine (a representative genotoxic PA and a proven rodent mutagen and carcinogen) are discussed in this review. Both of these compounds induced similar profiles of 6,7-dihydro-7-hydroxy-1-hydroxymethyl-5H-pyrrolizine (DHP)-derived DNA adducts and similar mutation spectra. Further, the two agents share common mechanisms of drug metabolism and carcinogenesis. Overall, comfrey is mutagenic in liver, and PA contained in comfrey appear to be responsible for comfrey-induced toxicity and tumor induction.

Scope: A methyl-deficient diet induces liver injury similar to human nonalcoholic steatohepatitis, one of the main risk factors for the development of hepatocellular carcinoma. Previous studies have demonstrated that this diet perturbs DNA methylation by causing a profound loss of global cytosine methylation, predominantly at heavily methylated repetitive sequences. However, whether methyl deficiency affects the methylation status of gene promoters has not been explored. Methods and results: Mouse gene expression and CpG island microarrays were used to characterize the gene expression and CpG island methylation profiles in the livers of C57BL/6J mice fed a methyl-deficient diet. We detected 164 genes that were differentially expressed and exhibited an inverse relationship between the gene expression and the extent of CpG island methylation. Furthermore, these genes were associated with altered lipid and glucose metabolism, DNA damage and repair, apoptosis, the development of fibrosis, and liver tissue remodeling. Although there were both increased and decreased levels of CpG island methylation, the number of hypomethylated genes was substantially greater than the number of hypermethylated genes. Conclusion: The results this study demonstrate that pairing methylation profiles with gene expression profiles is a powerful approach to identify dysregulated high-priority fundamental pathophysiological pathways associated with disease development.

Introduction: Rapid, accurate and inexpensive analysis of the disease-causing potential of foodborne pathogens is an important consideration in food safety and biodefense, particularly in developing countries. The objective of this study is to demonstrate the use of a robust and inexpensive microarray platform to assay the virulence gene profiles in Salmonella from food and/or the food animal environment, and then use ArrayTrackTM for data analysis. Methodology: The spotted array consisted of 69 selected Salmonella-specific virulence gene probes (65bp each). These probes were printed on poly-L-lysine-coated slides. Genomic DNA was digested with Sau3AI, labeled with Cy3 dye, hybridized to the gene probes, and the images were captured and analyzed by GenePix 4000B and ArrayTrackTM, a free software developed by Food and Drug Administration (FDA) researchers. Results: Nearly 58% of the virulence-associated genes tested were present in all Salmonella strains tested. In general, genes belonging to inv, pip, prg, sic, sip, spa or ttr families were detected in more than 90% of the isolates, while the iacP, avrA, invH, rhuM, sirA, sopB, sopE or sugR genes were detected in 40 to 80% of the isolates. The gene variability was independent of the Salmonella serotype. Conclusions: This hybridization array presents an accurate and cost-effective method for evaluating the disease-causing potential of Salmonella in outbreak investigations by targeting a selective set of Salmonella-associated virulence genes.

Gene expression data from microarrays are being applied to predict preclinical and clinical endpoints, but the reliability of these predictions has not been established. In the MAQC-II project, 36 independent teams analyzed six microarray data sets to generate predictive models for classifying a sample with respect to one of 13 endpoints indicative of lung or liver toxicity in rodents, or of breast cancer, multiple myeloma or neuroblastoma in humans. In total, >30,000 models were built using many combinations of analytical methods. The teams generated predictive models without knowing the biological meaning of some of the endpoints and, to mimic clinical reality, tested the models on data that had not been used for training. We found that model performance depended largely on the endpoint and team proficiency and that different approaches generated models of similar performance. The conclusions and recommendations from MAQC-II should be useful for regulatory agencies, study committees and independent investigators that evaluate methods for global gene expression analysis.

The rapid advances in proteomic analyses coupled with the completion of multiple genomes have led to an increased demand for determining protein functions. The first step is classification or prediction into families. A method was developed for the prediction of protein family based only on protein sequence using support vector machine (SVM) models. In these models, the amino acids were classified into three categories (apolar, polar, and charged). Consecutive fragments ranging from one to five were annotated by amino acid type to define the protein features of each protein. SVM models were constructed based on the protein features of a training set of proteins and then examined with an independent set of proteins. The approach was tested for 20 protein families from the iProClass database of Protein Information Resources (PIR). For two-class SVM models, an average prediction accuracy of 0.9985 was achieved, while for multi-class SVM models an accuracy of 0.9941 was achieved. This study demonstrates that SVM based methods can accurately recognize and predict the protein family to which a sequence belongs based solely on its primary amino acid sequence.

Usnic acid is a lichen metabolite used as a weight-loss dietary supplement due to its uncoupling action on mitochondria. However, its use has been associated with severe liver disorders in some individuals. Animal studies conducted thus far evaluated the effects of usnic acid on mitochondria primarily by measuring the rate of oxygen consumption and/or ATP generation. To obtain further insight into usnic acid-mediated effects on mitochondria, we examined the expression levels of 542 genes associated with mitochondrial structure and functions in liver of B6C3F(1) female mice using a mitochondria-specific microarray. Beginning at 8 weeks of age, mice received usnic acid at 0, 60, 180, and 600 ppm in ground, irradiated 5LG6 diet for 14 days. Microarray analysis showed a significant effect of usnic acid on the expression of several genes only at the highest dose of 600 ppm. A prominent finding of the study was a significant induction of genes associated with complexes I through IV of the electron transport chain. Moreover, several genes involved in fatty acid oxidation, the Krebs cycle, apoptosis, and membrane transporters were over-expressed. Usnic acid is a lipophilic weak acid that can diffuse through mitochondrial membranes and cause a proton leak (uncoupling). The up-regulation of complexes I-IV may be a compensatory mechanism to maintain the proton gradient across the mitochondrial inner membrane. In addition, induction of fatty acid oxidation and the Krebs cycle may be an adaptive response to uncoupling of mitochondria.

Zidovudine (3'-azido-3'-deoxythymidine; AZT) is the main anti-retroviral drug given to HIV-1-infected pregnant women during pregnancy and to their infants after birth to reduce mother-to-child transmission of the virus. In animal studies, however, a significant mitochondrial morphological damage has been reported in skeletal muscle as a consequence of transplacental or perinatal exposure to AZT. Because proper muscle function is highly dependent on efficient mitochondrial function and information on AZT-induced mitochondrial toxicity during neonatal exposure is limited, we investigated the effect of AZT on the expression of 542 mitochondria-related genes encoded by both nuclear and mitochondrial DNA in the skeletal muscle of infant male and female mice using microarray technology. Animals were treated orally by gavage with AZT at 0, 10, 50, 100, and 200mg/kg body weight/day from postnatal day (PND) 1 through 8 and were sacrificed at 1- and 2-h following the last dose on PND 8. These doses in mice correspond to 0, 1.1, 5.5, 11.0, and 22.0mg/kg AZT in human infants [Center for Drug Evaluation and Research (CDER) 2005. Pharmacology and Toxicology, Guidance for industry. Estimating the maximum safe dose in initial clinical trials for therapeutics in adult healthy volunteers, p. 7. http://www.fda.gov/cder/guidance/index.htm.]. Microarray data were analyzed for effects of time, sex, treatment, and their interactions using a fixed effect linear model. The results showed modest, but significant, dose-related responses in the expression level of genes associated with apoptosis, fatty acid metabolism, mitochondrial DNA maintenance, and various mitochondrial membrane transporters. The transcription levels were not significantly different at both time points and were not sex dependent. The results suggest that changes in expression of mitochondria-related genes in skeletal muscle may be an initial response to short-term AZT exposure in infant mice.

Drug-induced liver injury (DILI) is the primary adverse event that results in withdrawal of drugs from the market and a frequent reason for the failure of drug candidates in development. The Liver Toxicity Biomarker Study (LTBS) is an innovative approach to investigate DILI because it compares molecular events produced in vivo by compound pairs that (a) are similar in structure and mechanism of action, (b) are associated with few or no signs of liver toxicity in preclinical studies, and (c) show marked differences in hepatotoxic potential. The LTBS is a collaborative preclinical research effort in molecular systems toxicology between the National Center for Toxicological Research and BG Medicine, Inc., and is supported by seven pharmaceutical companies and three technology providers. In phase I of the LTBS, entacapone and tolcapone were studied in rats to provide results and information that will form the foundation for the design and implementation of phase II. Molecular analysis of the rat liver and plasma samples combined with statistical analyses of the resulting datasets yielded marker analytes, illustrating the value of the broad-spectrum, molecular systems analysis approach to studying pharmacological or toxicological effects.