Potential role of the apelin-APJ pathway in sex-related differential cardiotoxicity induced by doxorubicin in mice

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₁ mice indicated more differentially expressed genes in males than females (224 vs. 19; ≥1.5-fold, False Discovery Rate [FDR] < 0.05) at 1 week after receiving 24 mg/kg total cumulative DOX dose that induced cardiac lesions only in males. Pathway analysis further revealed probable inactivation of cardiac apelin fibroblast signaling pathway (p = 0.00004) only in DOX-treated male mice that showed ≥1.25-fold downregulation in the transcript and protein levels of the apelin receptor, APJ. In hearts of DOX-treated females, the transcript levels of apelin (1.24-fold) and APJ (1.47-fold) were significantly (p < 0.05) increased compared to saline-treated controls. Sex-related differential DOX effect was also observed on molecular targets downstream of the apelin-APJ pathway in cardiac fibroblasts and cardiomyocytes. In cardiac fibroblasts, upregulation of Tgf-β2, Ctgf, Sphk1, Serpine1, and Timp1 (fibrosis; FDR < 0.05) in DOX-treated males and upregulation of only Tgf-β2 and Timp1 (p < 0.05) in females suggested a greater DOX toxicity in hearts of males than females. Additionally, Ryr2 and Serca2 (calcium handling; FDR < 0.05) were downregulated in conjunction with 1.35-fold upregulation of Casp12 (sarcoplasmic reticulum-mediated apoptosis; FDR < 0.05) in DOX-treated male mice. Drug effect on the transcript level of these genes was less severe in female hearts. Collectively, these data suggest a likely role of the apelin-APJ axis in sex-related differential DOX-induced cardiotoxicity in our mouse model.

MicroRNA-34a-5p as a promising early circulating preclinical biomarker of doxorubicin-induced chronic cardiotoxicity

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 one 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 7 miRNAs were 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 (FDR<0.1) elevated at all total cumulative doses compared to 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<0.1) in plasma and heart compared to 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.

Doxorubicin-induced delayed-onset subclinical cardiotoxicity in mice

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 B6C3F₁ 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.

A comprehensive rat transcriptome built from large scale RNA-seq-based annotation

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.

Hepatic Transcript Profiles of Cytochrome P450 Genes Predict Sex Differences in Drug Metabolism

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 cytochrome P450 (P450) enzymes 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 (preweaning) to 104 weeks (elderly). Large differences in the transcript profiles of 29 drug metabolizing enzymes and transporters were found between adult males and females (8-52 weeks). Using the PharmaPendium data base, 41 drugs were found to be metabolized by one or two P450 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.

Candidate early predictive plasma protein markers of doxorubicin-induced chronic cardiotoxicity in B6C3F1 mice

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₁ mice given a weekly intravenous dose of 3 mg/kg DOX or saline (SAL) for 2, 3, 4, 6, or 8 weeks (6, 9, 12, 18, or 24 mg/kg cumulative DOX doses, respectively). Also, a group of mice received the cardio-protectant, dexrazoxane (DXZ; 60 mg/kg; intraperitoneal) 30 min before a weekly DOX or SAL dose. Proteomic analysis in plasma collected a week after the last dose showed a significant ≥1.2-fold change in level of 18 proteins in DOX-treated mice compared to SAL-treated counterparts during 8-week exposure. Of these, neurogenic locus notch homolog protein 1 (NOTCH1), von Willebrand factor (vWF), mitochondrial glutamate carrier 2, Wnt inhibitory factor 1, legumain, and mannan-binding lectin serine protease 1 were increased in plasma at 6 mg/kg cumulative DOX dose, prior to the release of myocardial injury marker, cardiac troponin I at 12 mg/kg and higher cumulative doses. These six proteins also remained significantly elevated following myocardial injury or pathology at 24 mg/kg. Pretreatment of mice with DXZ significantly attenuated DOX-induced elevated levels of only NOTCH1 and vWF with mitigation of cardiotoxicity. This suggests NOTCH1 and vWF as candidate early biomarkers of DOX cardiotoxicity, which may help in addressing a clinically important question of identifying cancer patients at risk for cardiotoxicity.

Transcript profiling in the testes and prostates of postnatal day 30 Sprague-Dawley rats exposed prenatally and lactationally to 2-hydroxy-4-methoxybenzophenone

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.

Transcriptomics analysis of early embryonic stem cell differentiation under osteoblast culture conditions: Applications for detection of developmental toxicity

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 72h, 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.

Sex and age differences in the expression of liver microRNAs during the life span of F344 rats

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) <5% and ≥1.5-fold change) sex- and age-related differences in expression, respectively. Thirty-eight miRNAs showed 2-week-specific expression, of which 31 miRNAs were found to be encoded within the Dlk1-Dio3 cluster located on chromosome 6. This cluster has been associated with tissue proliferation and differentiation, and liver energy homeostasis in postnatal development. Predictive pathway analysis linked sex-biased miRNA expression with sexually dimorphic molecular functions and toxicological functions that may reflect sex differences in hepatic physiology and disease. The expression of miRNAs (miR-18a, miR-99a, and miR-203, miR-451) was also found to associate with specific sexually dimorphic hepatic histopathology. The expression of miRNAs involved in regulating cell death, cell proliferation, and cell cycle was found to change as the rats matured from adult to old age.

Conclusions

Overall, significant age- and sex-related differences in liver miRNA expression were identified and linked to histopathological findings and predicted functional pathways that may underlie susceptibilities to liver toxicity and disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s13293-017-0127-9) contains supplementary material, which is available to authorized users.

Stably Expressed Genes Involved in Basic Cellular Functions

Stably Expressed Genes (SEGs) whose expression varies within a narrow range may be involved in core cellular processes necessary for basic functions. To identify such genes, we re-analyzed existing RNA-Seq gene expression profiles across 11 organs at 4 developmental stages (from immature to old age) in both sexes of F344 rats (n = 4/group; 320 samples). Expression changes (calculated as the maximum expression / minimum expression for each gene) of >19000 genes across organs, ages, and sexes ranged from 2.35 to >10⁹-fold, with a median of 165-fold. The expression of 278 SEGs was found to vary ≤4-fold and these genes were significantly involved in protein catabolism (proteasome and ubiquitination), RNA transport, protein processing, and the spliceosome. Such stability of expression was further validated in human samples where the expression variability of the homologous human SEGs was significantly lower than that of other genes in the human genome. It was also found that the homologous human SEGs were generally less subject to non-synonymous mutation than other genes, as would be expected of stably expressed genes. We also found that knockout of SEG homologs in mouse models was more likely to cause complete preweaning lethality than non-SEG homologs, corroborating the fundamental roles played by SEGs in biological development. Such stably expressed genes and pathways across life-stages suggest that tight control of these processes is important in basic cellular functions and that perturbation by endogenous (e.g., genetics) or exogenous agents (e.g., drugs, environmental factors) may cause serious adverse effects.

Sex-related differential susceptibility to doxorubicin-induced cardiotoxicity in B6C3F1 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₁ 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.

Status of hepatic DNA methylome predetermines and modulates the severity of non-alcoholic fatty liver injury in mice

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.

Early metabolomics changes in heart and plasma during chronic doxorubicin treatment in B6C3F1 mice

The present study aimed to identify molecular markers of early stages of cardiotoxicity induced by a potent chemotherapeutic agent, doxorubicin (DOX). Male B6C3F1 mice were dosed with 3 mg kg(-1) DOX or saline via tail vein weekly for 2, 3, 4, 6 or 8 weeks (cumulative DOX doses of 6, 9, 12, 18 or 24 mg kg(-1) , respectively) and euthanized a week after the last dose. Mass spectrometry-based and nuclear magnetic resonance spectrometry-based metabolic profiling were employed to identify initial biomarkers of cardiotoxicity before myocardial injury and cardiac pathology, which were not noted until after the 18 and 24 mg kg(-1) cumulative doses, respectively. After a cumulative dose of 6 mg kg(-1) , 18 amino acids and four biogenic amines (acetylornithine, kynurenine, putrescine and serotonin) were significantly increased in cardiac tissue; 16 amino acids and two biogenic amines (acetylornithine and hydroxyproline) were significantly altered in plasma. In addition, 16 acylcarnitines were significantly increased in plasma and five were significantly decreased in cardiac tissue compared to saline-treated controls. Plasma lactate and succinate, involved in the Krebs cycle, were significantly altered after a cumulative dose of 6 mg kg(-1) . A few metabolites remained altered at higher cumulative DOX doses, which could partly indicate a transition from injury processes at 2 weeks to repair processes with additional injury happening concurrently before myocardial injury at 8 weeks. These altered metabolic profiles in mouse heart and plasma during the initial stages of injury progression due to DOX treatment may suggest these metabolites as candidate early biomarkers of cardiotoxicity. Published 2016. This article is a U.S. Government work and is in the public domain in the USA.

A rat RNA-Seq transcriptomic BodyMap across 11 organs and 4 developmental stages

The rat has been used extensively as a model for evaluating chemical toxicities and for understanding drug mechanisms. However, its transcriptome across multiple organs, or developmental stages, has not yet been reported. Here we show, as part of the SEQC consortium efforts, a comprehensive rat transcriptomic BodyMap created by performing RNA-Seq on 320 samples from 11 organs of both sexes of juvenile, adolescent, adult and aged Fischer 344 rats. We catalogue the expression profiles of 40,064 genes, 65,167 transcripts, 31,909 alternatively spliced transcript variants and 2,367 non-coding genes/non-coding RNAs (ncRNAs) annotated in AceView. We find that organ-enriched, differentially expressed genes reflect the known organ-specific biological activities. A large number of transcripts show organ-specific, age-dependent or sex-specific differential expression patterns. We create a web-based, open-access rat BodyMap database of expression profiles with crosslinks to other widely used databases, anticipating that it will serve as a primary resource for biomedical research using the rat model.

Gene expression is highly variable between tissues, and changes during development and with age. Here, the authors provide a comprehensive RNA-Seq analysis of the rat transcriptome, spanning eleven organs, four developmental stages and both sexes.

Effects of oral exposure to bisphenol A on gene expression and global genomic DNA methylation in the prostate, female mammary gland, and uterus of NCTR Sprague-Dawley rats

Bisphenol A (BPA), an industrial chemical used in the manufacture of polycarbonate and epoxy resins, binds to the nuclear estrogen receptor with an affinity 4-5 orders of magnitude lower than that of estradiol. We reported previously that "high BPA" [100,000 and 300,000 µg/kg body weight (bw)/day], but not "low BPA" (2.5-2700 µg/kg bw/day), induced clear adverse effects in NCTR Sprague-Dawley rats gavaged daily from gestation day 6 through postnatal day (PND) 90. The "high BPA" effects partially overlapped those of ethinyl estradiol (EE2, 0.5 and 5.0 µg/kg bw/day). To evaluate further the potential of "low BPA" to induce biological effects, here we assessed the global genomic DNA methylation and gene expression in the prostate and female mammary glands, tissues identified previously as potential targets of BPA, and uterus, a sensitive estrogen-responsive tissue. Both doses of EE2 modulated gene expression, including of known estrogen-responsive genes, and PND 4 global gene expression data showed a partial overlap of the "high BPA" effects with those of EE2. The "low BPA" doses modulated the expression of several genes; however, the absence of a dose response reduces the likelihood that these changes were causally linked to the treatment. These results are consistent with the toxicity outcomes.

Persistence of furan-induced epigenetic aberrations in the livers of F344 rats

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.

Age and sex differences in kidney microRNA expression during the life span of F344 rats

Background

Growing evidence suggests that epigenetic mechanisms of gene regulation may play a role in susceptibilities to specific toxicities and adverse drug reactions. MiRNAs in particular have been shown to be important regulators in cancer and other diseases and show promise as predictive biomarkers for diagnosis and prognosis. In this study, we characterized the global kidney miRNA expression profile in untreated male and female F344 rats throughout the life span. These findings were correlated with sex-specific susceptibilities to adverse renal events, such as male-biased renal fibrosis and inflammation in old age.

Methods

Kidney miRNA expression was examined in F344 rats at 2, 5, 6, 8, 15, 21, 78, and 104 weeks of age in both sexes using Agilent miRNA microarrays. Differential expression was determined using filtering criteria of ≥1.5 fold change and ANOVA or pairwise t-test (FDR <5%) to determine significant age and sex effects, respectively. Pathway analysis software was used to investigate the possible roles of these target genes in age- and sex-specific differences.

Results

Three hundred eleven miRNAs were found to be expressed in at least one age and sex. Filtering criteria revealed 174 differentially expressed miRNAs in the kidney; 173 and 34 miRNAs exhibiting age and sex effects, respectively. Principal component analysis revealed age effects predominated over sex effects, with 2-week miRNA expression being much different from other ages. No significant sexually dimorphic miRNA expression was observed from 5 to 8 weeks, while the most differential expression (13 miRNAs) was observed at 21 weeks. Potential target genes of these differentially expressed miRNAs were identified.

Conclusions

The expression of 56% of detected renal miRNAs was found to vary significantly with age and/or sex during the life span of F344 rats. Pathway analysis suggested that 2-week-expressed miRNAs may be related to organ and cellular development and proliferation pathways. Male-biased miRNA expression at older ages correlated with male-biased renal fibrosis and mononuclear cell infiltration. These miRNAs showed high representation in renal inflammation and nephritis pathways, and included miR-214, miR-130b, miR-150, miR-223, miR-142-5p, miR-185, and miR-296*. Analysis of kidney miRNA expression throughout the rat life span will improve the use of current and future renal biomarkers and inform our assessments of kidney injury and disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s13293-014-0019-1) contains supplementary material, which is available to authorized users.

Sexual dimorphism in the expression of mitochondria-related genes in rat heart at different ages

Cardiovascular disease (CVD) is the leading cause of mortality worldwide. Moreover, sex and age are considered major risk factors in the development of CVDs. Mitochondria are vital for normal cardiac function, and regulation of mitochondrial structure and function may impact susceptibility to CVD. To identify potential role of mitochondria in sex-related differences in susceptibility to CVD, we analyzed the basal expression levels of mitochondria-related genes in the hearts of male and female rats. Whole genome expression profiling was performed in the hearts of young (8-week), adult (21-week), and old (78-week) male and female Fischer 344 rats and the expression of 670 unique genes related to various mitochondrial functions was analyzed. A significant (p<0.05) sexual dimorphism in expression levels of 46, 114, and 41 genes was observed in young, adult and old rats, respectively. Gene Ontology analysis revealed the influence of sex on various biological pathways related to cardiac energy metabolism at different ages. The expression of genes involved in fatty acid metabolism was significantly different between the sexes in young and adult rat hearts. Adult male rats also showed higher expression of genes associated with the pyruvate dehydrogenase complex compared to females. In young and adult hearts, sexual dimorphism was not noted in genes encoding oxidative phosphorylation. In old rats, however, a majority of genes involved in oxidative phosphorylation had higher expression in females compared to males. Such basal differences between the sexes in cardiac expression of genes associated with energy metabolism may indicate a likely involvement of mitochondria in susceptibility to CVDs. In addition, female rats showed lower expression levels of apoptotic genes in hearts compared to males at all ages, which may have implications for better preservation of cardiac mass in females than in males.

Simultaneous quantification of t(14;18) and HPRT exon 2/3 deletions in human lymphocytes

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.

Early biomarkers of doxorubicin-induced heart injury in a mouse model

Cardiac troponins, which are used as myocardial injury markers, are released in plasma only after tissue damage has occurred. Therefore, there is a need for identification of biomarkers of earlier events in cardiac injury to limit the extent of damage. To accomplish this, expression profiling of 1179 unique microRNAs (miRNAs) was performed in a chronic cardiotoxicity mouse model developed in our laboratory. Male B6C3F1 mice were injected intravenously with 3mg/kg doxorubicin (DOX; an anti-cancer drug), or saline once a week for 2, 3, 4, 6, and 8weeks, resulting in cumulative DOX doses of 6, 9, 12, 18, and 24mg/kg, respectively. Mice were euthanized a week after the last dose. Cardiac injury was evidenced in mice exposed to 18mg/kg and higher cumulative DOX dose whereas examination of hearts by light microscopy revealed cardiac lesions at 24mg/kg DOX. Also, 24 miRNAs were differentially expressed in mouse hearts, with the expression of 1, 1, 2, 8, and 21 miRNAs altered at 6, 9, 12, 18, and 24mg/kg DOX, respectively. A pro-apoptotic miR-34a was the only miRNA that was up-regulated at all cumulative DOX doses and showed a significant dose-related response. Up-regulation of miR-34a at 6mg/kg DOX may suggest apoptosis as an early molecular change in the hearts of DOX-treated mice. At 12mg/kg DOX, up-regulation of miR-34a was associated with down-regulation of hypertrophy-related miR-150; changes observed before cardiac injury. These findings may lead to the development of biomarkers of earlier events in DOX-induced cardiotoxicity that occur before the release of cardiac troponins.

MicroRNA-155 deficient mice experience heightened kidney toxicity when dosed with cisplatin

The development of nephrotoxicity limits the maximum achievable dosage and treatment intervals for cisplatin chemotherapy. Therefore, identifying mechanisms that regulate this toxicity could offer novel methods to optimize cisplatin delivery. MicroRNAs are capable of regulating many different genes, and can influence diverse cellular processes, including cell death and apoptosis. We previously observed miR-155 to be highly increased following ischemic or toxic injury to the kidneys and, therefore, sought to determine whether mice deficient in miR-155 would respond differently to kidney injury. We treated C57BL/6 and miR-155(-/-) mice with 20 mg/kg of cisplatin and found a significantly higher level of kidney injury in the miR-155(-/-) mice. Genome-wide expression profiling and bioinformatic analysis indicated the activation of a number of canonical signaling pathways relating to apoptosis and oxidative stress over the course of the injury, and identified potential upstream regulators of these effects. One predicted upstream regulator was c-Fos, which has two confirmed miR-155 binding sites in its 3' UTR and, therefore, can be directly regulated by miR-155. We established that the miR-155(-/-) mice had significantly higher levels of c-Fos mRNA and protein than the C57BL/6 mice at 72 h after cisplatin exposure. These data indicate a role for miR-155 in the cisplatin response and suggest that targeting of c-Fos could be investigated to reduce cisplatin-induced nephrotoxicity.

Differential gene expression in human hepatocyte cell lines exposed to the antiretroviral agent zidovudine

Zidovudine (3'-azido-3'-deoxythymidine; AZT) is the most widely used nucleoside reverse transcriptase inhibitor for the treatment of AIDS patients and prevention of mother-to-child transmission of HIV-1. Previously, we demonstrated that AZT had significantly greater growth inhibitory effects upon the human liver carcinoma cell line HepG2 as compared to the immortalized human liver cell line THLE2. We have now used gene expression profiling to determine the molecular pathways associated with toxicity in both cell lines. HepG2 cells were incubated with 0, 2, 20, or 100 μM AZT for 2 weeks; THLE2 cells were treated with 0, 50, 500, or 2,500 μM AZT, concentrations that were equi-toxic to those used in the HepG2 cells. After the treatment, total RNA was isolated and subjected to microarray analysis. Global analysis of gene expression, with a false discovery rate ≤0.01 and a fold change ≥1.5, indicated that 6- to 70-fold more genes were differentially expressed in a significant concentration-dependent manner in HepG2 cells when compared to THLE2 cells. Comparative analysis indicated that 7 % of these genes were common to both cell lines. Among the common differentially expressed genes, 70 % changed in the same direction, most of which were associated with cell death and survival, cell cycle, cell growth and proliferation, and DNA replication, recombination, and repair. As determined by the uptake of [methyl-(3)H]AZT, the intracellular levels of total AZT were approximately twofold higher in THLE2 cells than in HepG2 cells. The expression of thymidine kinase 1 (TK1) and UDP-glucuronosyltransferase 2B7 (UGT2B7) genes that regulate the metabolic activation and deactivation of AZT, respectively, was increased in HepG2 cells but decreased in THLE2 cells after treatment with AZT. This differential response in AZT metabolism was confirmed by real-time PCR, western blotting, and/or enzymatic assays. These data indicate that molecular pathways involved with cell death and survival, cell cycle, cell growth and proliferation, and DNA replication, recombination, and repair are involved in the toxicities associated with AZT in both human cell lines, and that the difference in expression of TK1 and UGT2B7 in response to AZT treatment in HepG2 cells and THLE2 cells might explain why HepG2 cells are more sensitive than THLE2 cells to the toxicity of AZT.

Metzincins and related genes in experimental renal ageing: towards a unifying fibrosis classifier across species

BACKGROUND

We have previously described a transcriptomic classifier consisting of metzincins and related genes (MARGS) discriminating kidneys and other organs with or without fibrosis from human biopsies. We now apply our MARGS-based algorithm to a rat model of age-associated interstitial renal fibrosis.

METHODS

Untreated Fisher 344 rats (n = 76) were sacrificed between 2 to 104 weeks of age. For gene expression studies, we used single colour (Cy3) Agilent Whole Rat Genome 4 × 44k microarrays; 4-5 animals of each sex were profiled at each of the following ages: 2, 5, 6, 8, 15, 21, 78 and 104 weeks. Intensity data were subjected to variance stabilization (www.Partek.com). Data were analysed with ANOVA and other statistical methods.

RESULTS

Sixty MARGS were differentially expressed across age groups. More MARGS were differentially expressed in older males than in older females. Principal component analysis showed gene expression induced segregation of age groups by sex from 6 to 104 weeks of age. The expression level of MMP7 correlated best with fibrosis grade. Severity of fibrosis was determined in 20 animals at 78 and 104 weeks of age. Expression values of 15 of 19 genes of the original classifier present on the Agilent array, in conjunction with linear discriminant analysis, was sufficient to correctly classify these 20 samples into non-fibrosis and fibrosis. Overrepresentation of MMP2 protein and CD44 protein in fibrosis was confirmed by immunofluorescence.

CONCLUSIONS

Based on these results and our previous work, the MARGS classifier represents a cross-organ and cross-species classifier of fibrosis irrespective of aetiology. This finding provides evidence for a common pathway leading to fibrosis and will help to design a PCR-based clinical test.

Transcriptomic responses provide a new mechanistic basis for the chemopreventive effects of folic acid and tributyrin in rat liver carcinogenesis

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 498, 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.

Toxicogenomics and cancer susceptibility: advances with next-generation sequencing

The aim of this review is to comprehensively summarize the recent achievements in the field of toxicogenomics and cancer research regarding genetic-environmental interactions in carcinogenesis and detection of genetic aberrations in cancer genomes by next-generation sequencing technology. Cancer is primarily a genetic disease in which genetic factors and environmental stimuli interact to cause genetic and epigenetic aberrations in human cells. Mutations in the germline act as either high-penetrance alleles that strongly increase the risk of cancer development, or as low-penetrance alleles that mildly change an individual's susceptibility to cancer. Somatic mutations, resulting from either DNA damage induced by exposure to environmental mutagens or from spontaneous errors in DNA replication or repair are involved in the development or progression of the cancer. Induced or spontaneous changes in the epigenome may also drive carcinogenesis. Advances in next-generation sequencing technology provide us opportunities to accurately, economically, and rapidly identify genetic variants, somatic mutations, gene expression profiles, and epigenetic alterations with single-base resolution. Whole genome sequencing, whole exome sequencing, and RNA sequencing of paired cancer and adjacent normal tissue present a comprehensive picture of the cancer genome. These new findings should benefit public health by providing insights in understanding cancer biology, and in improving cancer diagnosis and therapy.

Life cycle analysis of kidney gene expression in male F344 rats

Age is a predisposing condition for susceptibility to chronic kidney disease and progression as well as acute kidney injury that may arise due to the adverse effects of some drugs. Age-related differences in kidney biology, therefore, are a key concern in understanding drug safety and disease progression. We hypothesize that the underlying suite of genes expressed in the kidney at various life cycle stages will impact susceptibility to adverse drug reactions. Therefore, establishing changes in baseline expression data between these life stages is the first and necessary step in evaluating this hypothesis. Untreated male F344 rats were sacrificed at 2, 5, 6, 8, 15, 21, 78, and 104 weeks of age. Kidneys were collected for histology and gene expression analysis. Agilent whole-genome rat microarrays were used to query global expression profiles. An ANOVA (p<0.01) coupled with a fold-change>1.5 in relative mRNA expression, was used to identify 3,724 unique differentially expressed genes (DEGs). Principal component analyses of these DEGs revealed three major divisions in life-cycle renal gene expression. K-means cluster analysis identified several groups of genes that shared age-specific patterns of expression. Pathway analysis of these gene groups revealed age-specific gene networks and functions related to renal function and aging, including extracellular matrix turnover, immune cell response, and renal tubular injury. Large age-related changes in expression were also demonstrated for the genes that code for qualified renal injury biomarkers KIM-1, Clu, and Tff3. These results suggest specific groups of genes that may underlie age-specific susceptibilities to adverse drug reactions and disease. This analysis of the basal gene expression patterns of renal genes throughout the life cycle of the rat will improve the use of current and future renal biomarkers and inform our assessments of kidney injury and disease.

Development of doxorubicin-induced chronic cardiotoxicity in the B6C3F1 mouse model

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.

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

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.

Characterization of whole genome amplified (WGA) DNA for use in genotyping assay development

Background

Genotyping assays often require substantial amounts of DNA. To overcome the problem of limiting amounts of available DNA, Whole Genome Amplification (WGA) methods have been developed. The multiple displacement amplification (MDA) method using Φ29 polymerase has become the preferred choice due to its high processivity and low error rate. However, the uniformity and fidelity of the amplification process across the genome has not been extensively characterized.

Results

To assess amplification uniformity, we used array-based comparative genomic hybridization (aCGH) to evaluate DNA copy number variations (CNVs) in DNAs amplified by two MDA kits: GenomiPhi and REPLI-g. The Agilent Human CGH array containing nearly one million probes was used in this study together with DNAs from a normal subject and 2 cystic fibrosis (CF) patients. Each DNA sample was amplified 4 independent times and compared to its native unamplified DNA. Komogorov distances and Phi correlations showed a high consistency within each sample group. Less than 2% of the probes showed more than 2-fold CNV introduced by the amplification process. The two amplification kits, REPLI-g and GenomiPhi, generate very similar amplified DNA samples despite the differences between the unamplified and amplified DNA samples. The results from aCGH analysis indicated that there were no obvious CNVs in the CFTR gene region due to WGA when compared to unamplified DNA. This was confirmed by quantitative real-time PCR copy number assays at 10 locations within the CFTR gene. DNA sequencing analysis of a 2-kb region within the CFTR gene showed no mutations introduced by WGA.

Conclusion

The relatively high uniformity and consistency of the WGA process, coupled with the low replication error rate, suggests that WGA DNA may be suitable for accurate genotyping. Regions of the genome that were consistently under-amplified were found to contain higher than average GC content. Because of the consistent differences between the WGA DNA and the native unamplified DNA, characterization of the genomic region of interest, as described here, will be necessary to ensure the reliability of genotyping results from WGA DNA.

The liver toxicity biomarker study phase I: markers for the effects of tolcapone or entacapone

The Liver Toxicity Biomarker Study is a systems toxicology approach to discover biomarkers that are indicative of a drug's potential to cause human idiosyncratic drug-induced liver injury. In phase I, the molecular effects in rat liver and blood plasma induced by tolcapone (a "toxic" drug) were compared with the molecular effects in the same tissues by dosing with entacapone (a "clean" drug, similar to tolcapone in chemical structure and primary pharmacological mechanism). Two durations of drug exposure, 3 and 28 days, were employed. Comprehensive molecular analysis of rat liver and plasma samples yielded marker analytes for various drug-vehicle or drug-drug comparisons. An important finding was that the marker analytes associated with tolcapone only partially overlapped with marker analytes associated with entacapone, despite the fact that both drugs have similar chemical structures and the same primary pharmacological mechanism of action. This result indicates that the molecular analyses employed in the study are detecting substantial "off-target" markers for the two drugs. An additional interesting finding was the modest overlap of the marker data sets for 3-day exposure and 28-day exposure, indicating that the molecular changes in liver and plasma caused by short- and long-term drug treatments do not share common characteristics.

Modeling chemical interaction profiles: II. Molecular docking, spectral data-activity relationship, and structure-activity relationship models for potent and weak inhibitors of cytochrome P450 CYP3A4 isozyme

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 ¹³C-NMR and 1D ¹⁵N-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.

Modeling chemical interaction profiles: I. Spectral data-activity relationship and structure-activity relationship models for inhibitors and non-inhibitors of cytochrome P450 CYP3A4 and CYP2D6 isozymes

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 pollutant—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 Mold² 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.

Expression analysis of hepatic mitochondria-related genes in mice exposed to acrylamide and glycidamide

Acrylamide (AA) is an industrial chemical that has been extensively investigated for central nervous system (CNS), reproductive, and genetic toxicity. However, AA effects on the liver, a major organ of drug metabolism, have not been adequately explored. In addition, the role of mitochondria in AA-mediated toxicity is still unclear. Changes in expression levels of genes associated with hepatic mitochondrial function of male transgenic Big Blue (BB) mice administered 500 mg/L AA or an equimolar concentration (600 mg/L) of its reactive metabolite glycidamide (GA) in drinking water for 3 and 4 wk, respectively, were examined. Transcriptional profiling of 542 mitochondria-related genes indicated a significant downregulation of genes associated with the 3-beta-hydroxysteroid dehydrogenase family in AA- and GA-treated mice, suggesting a possible role of both chemicals in altering hepatic steroid metabolism in BB mice. In addition, genes associated with lipid metabolism were altered by both treatments. Interestingly, only the parental compound (AA) significantly induced expression levels of genes associated with oxidative phosphorylation, in particular ATP synthase, which correlated with elevated ATP levels, indicating an increased energy demand in liver during AA exposure. Acrylamide-treated mice also showed significantly higher activity of glutathione S-transferase in association with decreased levels of reduced glutathione (GSH), which may imply an enhanced rate of conjugation of AA with GSH in liver. These results suggest different hepatic mechanisms of action of AA and GA and provide important insights into the involvement of mitochondria during their exposures.

Identification and categorization of liver toxicity markers induced by a related pair of drugs

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 and their predicted target messenger RNAs are deregulated by exposure to a carcinogenic dose of comfrey in rat liver

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.

Age and sex dependent changes in liver gene expression during the life cycle of the rat

Background

Age- and sex-related susceptibility to adverse drug reactions and disease is a key concern in understanding drug safety and disease progression. We hypothesize that the underlying suite of hepatic genes expressed at various life cycle stages will impact susceptibility to adverse drug reactions. Understanding the basal liver gene expression patterns is a necessary first step in addressing this hypothesis and will inform our assessments of adverse drug reactions as the liver plays a central role in drug metabolism and biotransformation. Untreated male and female F344 rats were sacrificed at 2, 5, 6, 8, 15, 21, 52, 78, and 104 weeks of age. Liver tissues were collected for histology and gene expression analysis. Whole-genome rat microarrays were used to query global expression profiles.

Results

An initial list of differentially expressed genes was selected using criteria based upon p-value (p < 0.05) and fold-change (+/- 1.5). Three dimensional principal component analyses revealed differences between males and females beginning at 2 weeks with more divergent profiles beginning at 5 weeks. The greatest sex-differences were observed between 8 and 52 weeks before converging again at 104 weeks. K-means clustering identified groups of genes that displayed age-related patterns of expression. Various adult aging-related clusters represented gene pathways related to xenobiotic metabolism, DNA damage repair, and oxidative stress.

Conclusions

These results suggest an underlying role for genes in specific clusters in potentiating age- and sex-related differences in susceptibility to adverse health effects. Furthermore, such a comprehensive picture of life cycle changes in gene expression deepens our understanding and informs the utility of liver gene expression biomarkers.

Coupling global methylation and gene expression profiles reveal key pathophysiological events in liver injury induced by a methyl-deficient diet

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.

Metabolism, genotoxicity, and carcinogenicity of comfrey

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.

Microarray analysis of virulence gene profiles in Salmonella serovars from food/food animal environment

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. 

The accurate prediction of protein family from amino acid sequence by measuring features of sequence fragments

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.

Application of microarray-based analysis of gene expression in the field of toxicogenomics

The field of toxicogenomics, which is becoming an important sub-discipline of toxicology, resulted from the natural convergence of the field of conventional toxicological research and the emergent field of functional genomics. One technology that has played a significant role in the field of toxicogenomics (in addition to many others) is the gene expression microarray. In this chapter, the authors provide an example of the application of gene expression microarrays to the field of toxicogenomics by detailing the strategy that was used for obtaining, analyzing, and interpreting gene expression data generated from RNA isolated from the liver of toxicant-exposed rats.

Effect of (+)-usnic acid on mitochondrial functions as measured by mitochondria-specific oligonucleotide microarray in liver of B6C3F1 mice

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.

The Liver Toxicity Biomarker Study: Phase I Design and Preliminary Results

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.

Effect of short-term exposure to zidovudine (AZT) on the expression of mitochondria-related genes in skeletal muscle of neonatal mice

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.

Designing Toxicogenomics Studies that use DNA Array Technology

Background

Bioassays are routinely used to evaluate the toxicity of test agents. Experimental designs for bioassays are largely encompassed by fixed effects linear models. In toxicogenomics studies where DNA arrays measure mRNA levels, the tissue samples are typically generated in a bioassay. These measurements introduce additional sources of variation, which must be properly managed to obtain valid tests of treatment effects.

Results

An analysis of covariance model is developed which combines a fixed-effects linear model for the bioassay with important variance components associated with DNA array measurements. These models can accommodate the dominant characteristics of measurements from DNA arrays, and they account for technical variation associated with normalization, spots, dyes, and batches as well as the biological variation associated with the bioassay. An example illustrates how the model is used to identify valid designs and to compare competing designs.

Conclusions

Many toxicogenomics studies are bioassays which measure gene expression using DNA arrays. These studies can be designed and analyzed using standard methods with a few modifications to account for characteristics of array measurements, such as multiple endpoints and normalization. As much as possible, technical variation associated with probes, dyes, and batches are managed by blocking treatments within these sources of variation. An example shows how some practical constraints can be accommodated by this modelling and how it allows one to objectively compare competing designs.

Designing toxicogenomics studies that use DNA array technology

BACKGROUND

Bioassays are routinely used to evaluate the toxicity of test agents. Experimental designs for bioassays are largely encompassed by fixed effects linear models. In toxicogenomics studies where DNA arrays measure mRNA levels, the tissue samples are typically generated in a bioassay. These measurements introduce additional sources of variation, which must be properly managed to obtain valid tests of treatment effects.

RESULTS

An analysis of covariance model is developed which combines a fixed-effects linear model for the bioassay with important variance components associated with DNA array measurements. These models can accommodate the dominant characteristics of measurements from DNA arrays, and they account for technical variation associated with normalization, spots, dyes, and batches as well as the biological variation associated with the bioassay. An example illustrates how the model is used to identify valid designs and to compare competing designs.

CONCLUSIONS

Many toxicogenomics studies are bioassays which measure gene expression using DNA arrays. These studies can be designed and analyzed using standard methods with a few modifications to account for characteristics of array measurements, such as multiple endpoints and normalization. As much as possible, technical variation associated with probes, dyes, and batches are managed by blocking treatments within these sources of variation. An example shows how some practical constraints can be accommodated by this modelling and how it allows one to objectively compare competing designs.

The balance of reproducibility, sensitivity, and specificity of lists of differentially expressed genes in microarray studies

BACKGROUND

Reproducibility is a fundamental requirement in scientific experiments. Some recent publications have claimed that microarrays are unreliable because lists of differentially expressed genes (DEGs) are not reproducible in similar experiments. Meanwhile, new statistical methods for identifying DEGs continue to appear in the scientific literature. The resultant variety of existing and emerging methods exacerbates confusion and continuing debate in the microarray community on the appropriate choice of methods for identifying reliable DEG lists.

RESULTS

Using the data sets generated by the MicroArray Quality Control (MAQC) project, we investigated the impact on the reproducibility of DEG lists of a few widely used gene selection procedures. We present comprehensive results from inter-site comparisons using the same microarray platform, cross-platform comparisons using multiple microarray platforms, and comparisons between microarray results and those from TaqMan - the widely regarded "standard" gene expression platform. Our results demonstrate that (1) previously reported discordance between DEG lists could simply result from ranking and selecting DEGs solely by statistical significance (P) derived from widely used simple t-tests; (2) when fold change (FC) is used as the ranking criterion with a non-stringent P-value cutoff filtering, the DEG lists become much more reproducible, especially when fewer genes are selected as differentially expressed, as is the case in most microarray studies; and (3) the instability of short DEG lists solely based on P-value ranking is an expected mathematical consequence of the high variability of the t-values; the more stringent the P-value threshold, the less reproducible the DEG list is. These observations are also consistent with results from extensive simulation calculations.

CONCLUSION

We recommend the use of FC-ranking plus a non-stringent P cutoff as a straightforward and baseline practice in order to generate more reproducible DEG lists. Specifically, the P-value cutoff should not be stringent (too small) and FC should be as large as possible. Our results provide practical guidance to choose the appropriate FC and P-value cutoffs when selecting a given number of DEGs. The FC criterion enhances reproducibility, whereas the P criterion balances sensitivity and specificity.