Development of mitochondria-specific mouse oligonucleotide microarray and validation of data by real-time PCR

In Vitro Modulation of Redox and Metabolism Interplay at the Brain Vascular Endothelium: Genomic and Proteomic Profiles of Sulforaphane Activity

Sulforaphane (SFN) has been shown to protect the brain vascular system and effectively reduce ischemic injuries and cognitive deficits. Given the robust cerebrovascular protection afforded by SFN, the objective of this study was to profile these effects in vitro using primary mouse brain microvascular endothelial cells and focusing on cellular redox, metabolism and detoxification functions. We used a mouse MitoChip array developed and validated at the FDA National Center for Toxicological Research (NCTR) to profile a host of genes encoded by nuclear and mt-DNA following SFN treatment (0-5 µM). Corresponding protein expression levels were assessed (ad hoc) by qRT-PCR, immunoblots and immunocytochemistry (ICC). Gene ontology clustering revealed that SFN treatment (24 h) significantly up-regulated ~50 key genes (>1.5 fold, adjusted p < 0.0001) and repressed 20 genes (<0.7 fold, adjusted p < 0.0001) belonging to oxidative stress, phase 1 & 2 drug metabolism enzymes (glutathione system), iron transporters, glycolysis, oxidative phosphorylation (OXPHOS), amino acid metabolism, lipid metabolism and mitochondrial biogenesis. Our results show that SFN stimulated the production of ATP by promoting the expression and activity of glucose transporter-1, and glycolysis. In addition, SFN upregulated anti-oxidative stress responses, redox signaling and phase 2 drug metabolism/detoxification functions, thus elucidating further the previously observed neurovascular protective effects of this compound.

Molecular signatures and functional analysis of beige adipocytes induced from in vivo intra-abdominal adipocytes

Beige adipocytes can be induced from white adipocytes and precursors upon stimulation by cold temperatures and act like brown adipocytes to increase energy expenditure. Most in vivo studies examining the mechanisms for the induction of beige adipocytes have focused on subcutaneous white adipose tissue (sWAT; benign fat) in the mouse. How intra-abdominal WAT (aWAT; malignant fat) develops into beige adipocytes remains obscure, largely because there is a lack of a good animal model for the induction of beige adipocytes from aWAT. To better understand the development of beige adipocytes from mammalian WATs, especially aWAT, we induced beige adipocytes from bat aWAT and mouse sWAT by exposure to cold temperatures and analyzed their molecular signatures. RNA sequencing followed by whole genome-wide expression analysis shows that beige adipocytes induced from bat aWAT, rather than sWAT, have molecular signatures resembling those of mouse sWAT-induced beige adipocytes and exhibit dynamic profiles similar to those of classical brown adipocytes. In addition, we identified molecular markers that were highly enriched in beige adipocytes and conserved between bat aWAT and mouse sWAT, a set that included the genes Uqcrc1 and Letm1. Furthermore, knockdown of Uqcrc1 and Letm1 expression shows that they are required not only for beige adipocyte differentiation but also for preadipocyte maturation. This study presents a new model for research into the induction of beige adipocytes from aWAT in vivo, which, when combined with models where beige adipocytes are induced from sWAT, provides insight into therapeutic approaches for combating obesity-related diseases in humans.

Early transcriptional changes in cardiac mitochondria during chronic doxorubicin exposure and mitigation by dexrazoxane in mice

Identification of early biomarkers of cardiotoxicity could help initiate means to ameliorate the cardiotoxic actions of clinically useful drugs such as doxorubicin (DOX). Since DOX has been shown to target mitochondria, transcriptional levels of mitochondria-related genes were evaluated to identify early candidate biomarkers in hearts of male B6C3F1 mice given a weekly intravenous dose of 3mg/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 was pretreated (intraperitoneally) with the cardio-protectant, dexrazoxane (DXZ; 60 mg/kg) 30 min before each weekly dose of DOX or SAL. At necropsy a week after the last dose, increased plasma concentrations of cardiac troponin T (cTnT) were detected at 18 and 24 mg/kg cumulative DOX doses, whereas myocardial alterations were observed only at the 24 mg/kg dose. Of 1019 genes interrogated, 185, 109, 140, 184, and 451 genes were differentially expressed at 6, 9, 12, 18, and 24 mg/kg cumulative DOX doses, respectively, compared to concurrent SAL-treated controls. Of these, expression of 61 genes associated with energy metabolism and apoptosis was significantly altered before and after occurrence of myocardial injury, suggesting these as early genomics markers of cardiotoxicity. Much of these DOX-induced transcriptional changes were attenuated by pretreatment of mice with DXZ. Also, DXZ treatment significantly reduced plasma cTnT concentration and completely ameliorated cardiac alterations induced by 24 mg/kg cumulative DOX. This information on early transcriptional changes during DOX treatment may be useful in designing cardioprotective strategies targeting mitochondria.

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.

Evaluation of Hepatic Mitochondria and Hematological Parameters in Zidovudine-Treated B6C3F(1) Mice

The effects of 12-week exposure to zidovudine (AZT) at 400, 500, and 600 mg/kg/d were examined on expression of 542 mitochondria-related genes and mitochondrial DNA (mtDNA) copy number in the liver of male and female B6C3F₁ mice to understand mitochondrial role in sex-related differences in development of lactic acidosis. Plasma lactate levels and hematologic parameters were also examined. Results indicated increased red blood cell (RBC) count in vehicle-treated controls, whereas a dose-related decline in the RBC count was noted in AZT-treated mice compared to the basal levels before treatments began. These decreases were associated with significant dose-related increases in mean corpuscular volume and mean corpuscular hemoglobin levels. This effect was greater in AZT-treated females compared to males. In both sexes, 12-week AZT or vehicle exposure significantly reduced plasma lactate levels compared to the basal levels. Results also showed modest, but significant, changes in the expression of genes associated with apoptosis and lipid metabolism at 600 mg/kg/d AZT. Neither drug nor sex influenced hepatic mtDNA copy number. Altogether, 12-week AZT exposure as high as 600 mg/kg/d did not impair hepatic mitochondria or induce lactic acidosis in B6C3F₁ mice. However, AZT-mediated hematologic toxicity appeared to be greater in females compared to males.

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.

Rat mitochondrion-neuron focused microarray (rMNChip) and bioinformatics tools for rapid identification of differential pathways in brain tissues

Mitochondrial function is of particular importance in brain because of its high demand for energy (ATP) and efficient removal of reactive oxygen species (ROS). We developed rat mitochondrion-neuron focused microarray (rMNChip) and integrated bioinformatics tools for rapid identification of differential pathways in brain tissues. rMNChip contains 1,500 genes involved in mitochondrial functions, stress response, circadian rhythms and signal transduction. The bioinformatics tool includes an algorithm for computing of differentially expressed genes, and a database for straightforward and intuitive interpretation for microarray results. Our application of these tools to RNA samples derived from rat frontal cortex (FC), hippocampus (HC) and hypothalamus (HT) led to the identification of differentially-expressed signal-transduction-bioenergenesis and neurotransmitter-synthesis pathways with a dominant number of genes (FC/HC = 55/6; FC/HT = 55/4) having significantly (p<0.05, FDR<10.70%) higher (≥1.25 fold) RNA levels in the frontal cortex than the others, strongly suggesting active generation of ATP and neurotransmitters and efficient removal of ROS. Thus, these tools for rapid and efficient identification of differential pathways in brain regions will greatly facilitate our systems-biological study and understanding of molecular mechanisms underlying complex and multifactorial neurodegenerative diseases.

Mitochondrial gene profiling: translational perspectives

The last decade has witnessed the development of multiple microarray platforms designed to study, in a comprehensive fashion, the expression and sequence of both mitochondrial and nuclear genes that encode mitochondrial proteins. Mitochondrial dysfunction has been implicated in a number of severe medical conditions including cancer, metabolic diseases (i.e., cardiovascular, diabetes and obesity) and neurodegenerative disorders and it is responsible for the adverse effects of numerous drugs. Profiling of the genetic and genomic status of mitochondria with focused microarrays offers the promise of rapidly and robustly identifying novel biomarkers for early disease diagnoses and prognoses, predicting of drug safety, liability, and selecting and stratifying of patients in clinical trials.

Gene expression profiling as an initial approach for mechanistic studies of toxicity and tumorigenicity of herbal plants and herbal dietary supplements

Dietary supplements are consumed by more than 300 million people worldwide, and herbal dietary supplements represent the most rapidly growing portion of this industry. Even though adverse health effects of many herbal dietary supplements have been reported, safety assurances are not being addressed adequately. Toxicological data on the identification of genotoxic and tumorigenic ingredients in many raw herbs are also lacking. Currently, more than 30 herbal dietary supplements and active ingredients have been selected by the National Toxicology Program (NTP) for toxicity and tumorigenicity studies. Due to the complexity of the chemical components present in plant extracts, there are no established methodologies for determining the mechanisms of toxicity (particularly tumorigenicity) induced by herbs, such as Gingko biloba leaf extract (GBE) and other herbal plant extracts. Consequently, the understanding of toxicity of herbal dietary supplements remains limited. We have proposed that application of DNA microarrays could be a highly practical initial approach for revealing biological pathways and networks associated with toxicity induced by herbal dietary supplements and the generation of hypotheses to address likely mechanisms. The changes in expression of subsets of genes of interest, such as the modulation of drug metabolizing genes, can be analyzed after treatment with an herbal dietary supplement. Although levels of gene expression do not represent fully the levels of protein activities, we propose that subsequent biochemical and genomic experiments based on these initial observations will enable elucidation of the mechanisms leading to toxicity, including tumorigenicity. This review summarizes the current practices of microarray analysis of gene expressions in animals treated with herbal dietary supplements and discusses perspectives for the proposed strategy.

Transgenic expression of proximal tubule peroxisome proliferator-activated receptor-alpha in mice confers protection during acute kidney injury

Our previous studies suggest that peroxisome proliferator-activated receptor-alpha (PPARalpha) plays a critical role in regulating fatty acid beta-oxidation in kidney tissue and this directly correlated with preservation of kidney morphology and function during acute kidney injury. To further study this, we generated transgenic mice expressing PPARalpha in the proximal tubule under the control of the promoter of KAP2 (kidney androgen-regulated protein 2). Segment-specific upregulation of PPARalpha expression by testosterone treatment of female transgenic mice improved kidney function during cisplatin or ischemia-reperfusion-induced acute kidney injury. Ischemia-reperfusion injury or treatment with cisplatin in wild-type mice caused inhibition of fatty-acid oxidation, reduction of mitochondrial genes of oxidative phosphorylation, mitochondrial DNA, fatty-acid metabolism, and the tricarboxylic acid cycle. Similar injury in testosterone-treated transgenic mice resulted in amelioration of these effects. Similarly, there were increases in the levels of 4-hydroxy-2-hexenal-derived lipid peroxidation products in wild-type mice, which were also reduced in the transgenic mice. Similarly, necrosis of the S3 segment was reduced in the two injury models in transgenic mice compared to wild type. Our results suggest proximal tubule PPARalpha activity serves as a metabolic sensor. Its increased expression without the use of an exogenous PPARalpha ligand in the transgenic mice is sufficient to protect kidney function and morphology, and to prevent abnormalities in lipid metabolism associated with acute kidney injury.

Toxicogenomic biomarkers for liver toxicity

Toxicogenomics (TGx) is a widely used technique in the preclinical stage of drug development to investigate the molecular mechanisms of toxicity. A number of candidate TGx biomarkers have now been identified and are utilized for both assessing and predicting toxicities. Further accumulation of novel TGx biomarkers will lead to more efficient, appropriate and cost effective drug risk assessment, reinforcing the paradigm of the conventional toxicology system with a more profound understanding of the molecular mechanisms of drug-induced toxicity. In this paper, we overview some practical strategies as well as obstacles for identifying and utilizing TGx biomarkers based on microarray analysis. Since clinical hepatotoxicity is one of the major causes of drug development attrition, the liver has been the best documented target organ for TGx studies to date, and we therefore focused on information from liver TGx studies. In this review, we summarize the current resources in the literature in regard to TGx studies of the liver, from which toxicologists could extract potential TGx biomarker gene sets for better hepatotoxicity risk assessment.

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.

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.

Testing for treatment effects on gene ontology

In studies that use DNA arrays to assess changes in gene expression, it is preferable to measure the significance of treatment effects on a group of genes from a pathway or functional category such as gene ontology terms (GO terms, ) because this facilitates the interpretation of effects and may markedly increase significance. A modified meta-analysis method to combine p-values was developed to measure the significance of an overall treatment effect on such functionally-defined groups of genes, taking into account the correlation structure among genes. For hypothesis testing that allows gene expression to change in both directions, p-values are calculated under the null distribution generated by a Monte Carlo method.

As a test of this procedure, we attempted to distinguish altered pathways in microarray studies performed with Mitochips, oligonucleotide microarrays specific to mitochondrial DNA-encoded transcripts. We found that our analytic method improves the specificity of selection for altered pathways, due to incorporation of the inter-gene correlation structure in each pathway. It is thus a practical method to measure treatment effects on GO groups. In many actual applications, microarray experiments measure treatment effects under complicated design structures and with small sample sizes. For such applications to real data of limited statistical power, and also in computer simulations, we demonstrate that our method gives reasonable test results.

Nucleoside reverse transcriptase inhibitors (NRTIs)-induced expression profile of mitochondria-related genes in the mouse liver

Mitochondrial dysfunction has been implicated in the adverse effects of nucleoside reverse transcriptase inhibitors (NRTIs) used to treat HIV-1 infections. To gain insight into the mechanism by which NRTIs alter mitochondrial function, the expression level of 542 genes associated with mitochondrial structure and functions was determined in the livers of p53 haplodeficient (+/-) C3B6F1 female mouse pups using mouse mitochondria-specific oligonucleotide microarray. The pups were transplacentally exposed to zidovudine (AZT) at 240 mg/kg bw/day or a combination of AZT and lamivudine (3TC) at 160 and 100mg/kg bw/day, respectively, from gestation day 12 through 18, followed by continuous treatment by oral administration from postnatal day 1-28. In addition, AZT/3TC effect was investigated in wild-type (+/+) C3B6F1 female mice. The genotype did not significantly affect the gene expression profile induced by AZT/3TC treatment. However, the transcriptional level of several genes associated with oxidative phosphorylation, mitochondrial tRNAs, fatty acid oxidation, steroid biosynthesis, and a few transport proteins were significantly altered in pups treated with AZT and AZT/3TC compared to their vehicle counterparts. Interestingly, AZT/3TC altered the expression level of 153 genes with false discovery rate of less than 0.05, in contrast to only 20 genes by AZT alone. These results suggest that NRTI-related effect on expression level of genes associated with mitochondrial functions was much greater in response to AZT/3TC combination treatment than AZT alone.