The genomic response of skeletal muscle to methylprednisolone using microarrays: tailoring data mining to the structure of the pharmacogenomic time series

Delay the progression of glucocorticoid-induced osteoporosis: Fraxin targets ferroptosis via the Nrf2/GPX4 pathway

Glucocorticoid-induced osteoporosis (GIOP) commonly accelerates bone loss, increasing the risk of fractures and osteonecrosis more significantly than traditional menopausal osteoporosis. The extracellular environment influenced by glucocorticoids heightens fracture and osteonecrosis risks. Fraxin (Fra), a key component of the traditional Chinese herbal remedy Cortex Fraxini, is known for its wide-ranging pharmacological effects, but its impact on GIOP remains unexplored. This investigation aims to delineate the effects and underlying mechanisms of Fra in combating dexamethasone (Dex)-induced ferroptosis and GIOP. We established a mouse model of GIOP via intraperitoneal injections of Dex and cultured osteoblasts with Dex treatment for in vitro analysis. We evaluated the impact of Fra on Dex-treated osteoblasts through assays such as C11-BODIPY and FerroOrange staining, mitochondrial functionality tests, and protein expression analyses via Western blot and immunofluorescence. The influence of Fra on bone microarchitecture of GIOP in mice was assessed using microcomputerized tomography, hematoxylin and eosin staining, double-labeling with Calcein-Alizarin Red S, and immunohistochemistry at imaging and histological levels. Based on our data, Fra prevented Dex-induced ferroptosis and bone loss. In vitro, glutathione levels increased and malondialdehyde, lipid peroxidation, and mitochondrial reactive oxygen species decreased. Fra treatment also increases nuclear factor erythroid 2-related factor 2 (Nrf2), glutathione peroxidase 4 (GPX4), and COL1A1 expression and promotes bone formation. To delve deeper into the mechanism, the findings revealed that Fra triggered the activation of Nrf2/GPX4 signaling. Moreover, the use of siRNA-Nrf2 blocked the beneficial effect of Fra in osteoblasts cultivated with Dex. Fra effectively combats GIOP by activating the Nrf2/GPX4 signaling pathway to inhibit ferroptosis.

Pathway-level analysis of genome-wide circadian dynamics in diverse tissues in rat and mouse

A computational framework is developed to enable the characterization of genome-wide, multi-tissue circadian dynamics at the level of "functional groupings of genes" defined in the context of signaling, cellular/genetic processing and metabolic pathways in rat and mouse. Our aim is to identify how individual genes come together to generate orchestrated rhythmic patterns and how these may vary within and across tissues. We focus our analysis on four tissues (adipose, liver, lung, and muscle). A genome-wide pathway-centric analysis enables us to develop a comprehensive picture on how the observed circadian variation at the individual gene level, orchestrates functional responses at the pathway level. Such pathway-based "meta-data" analysis enables the rational integration and comparison across platforms and/or experimental designs evaluating emergent dynamics, as opposed to comparisons of individual elements. One of our key findings is that when considering the dynamics at the pathway level, a complex behavior emerges. Our work proposes that tissues tend to coordinate gene's circadian expression in a way that optimizes tissue-specific pathway activity, depending of each tissue's broader role in homeostasis.

Transitioning from Basic toward Systems Pharmacodynamic Models: Lessons from Corticosteroids

Technology in bioanalysis, -omics, and computation have evolved over the past half century to allow for comprehensive assessments of the molecular to whole body pharmacology of diverse corticosteroids. Such studies have advanced pharmacokinetic and pharmacodynamic (PK/PD) concepts and models that often generalize across various classes of drugs. These models encompass the "pillars" of pharmacology, namely PK and target drug exposure, the mass-law interactions of drugs with receptors/targets, and the consequent turnover and homeostatic control of genes, biomarkers, physiologic responses, and disease symptoms. Pharmacokinetic methodology utilizes noncompartmental, compartmental, reversible, physiologic [full physiologically based pharmacokinetic (PBPK) and minimal PBPK], and target-mediated drug disposition models using a growing array of pharmacometric considerations and software. Basic PK/PD models have emerged (simple direct, biophase, slow receptor binding, indirect response, irreversible, turnover with inactivation, and transduction models) that place emphasis on parsimony, are mechanistic in nature, and serve as highly useful "top-down" methods of quantitating the actions of diverse drugs. These are often components of more complex quantitative systems pharmacology (QSP) models that explain the array of responses to various drugs, including corticosteroids. Progressively deeper mechanistic appreciation of PBPK, drug-target interactions, and systems physiology from the molecular (genomic, proteomic, metabolomic) to cellular to whole body levels provides the foundation for enhanced PK/PD to comprehensive QSP models. Our research based on cell, animal, clinical, and theoretical studies with corticosteroids have provided ideas and quantitative methods that have broadly advanced the fields of PK/PD and QSP modeling and illustrates the transition toward a global, systems understanding of actions of diverse drugs. SIGNIFICANCE STATEMENT: Over the past half century, pharmacokinetics (PK) and pharmacokinetics/pharmacodynamics (PK/PD) have evolved to provide an array of mechanism-based models that help quantitate the disposition and actions of most drugs. We describe how many basic PK and PK/PD model components were identified and often applied to the diverse properties of corticosteroids (CS). The CS have complications in disposition and a wide array of simple receptor-to complex gene-mediated actions in multiple organs. Continued assessments of such complexities have offered opportunities to develop models ranging from simple PK to enhanced PK/PD to quantitative systems pharmacology (QSP) that help explain therapeutic and adverse CS effects. Concurrent development of state-of-the-art PK, PK/PD, and QSP models are described alongside experimental studies that revealed diverse CS actions.

Modeling Pathway Dynamics of the Skeletal Muscle Response to Intravenous Methylprednisolone (MPL) Administration in Rats: Dosing and Tissue Effects

A model-based approach for the assessment of pathway dynamics is explored to characterize metabolic and signaling pathway activity changes characteristic of the dosing-dependent differences in response to methylprednisolone in muscle. To consistently compare dosing-induced changes we extend the principles of pharmacokinetics and pharmacodynamics and introduce a novel representation of pathway-level dynamic models of activity regulation. We hypothesize the emergence of dosing-dependent regulatory interactions is critical to understanding the mechanistic implications of MPL dosing in muscle. Our results indicate that key pathways, including amino acid and lipid metabolism, signal transduction, endocrine regulation, regulation of cellular functions including growth, death, motility, transport, protein degradation, and catabolism are dependent on dosing, exhibiting diverse dynamics depending on whether the drug is administered acutely of continuously. Therefore, the dynamics of drug presentation offer the possibility for the emergence of dosing-dependent models of regulation. Finally, we compared acute and chronic MPL response in muscle with liver. The comparison revealed systematic response differences between the two tissues, notably that muscle appears more prone to adapt to MPL.

Pathway-Based Analysis of the Liver Response to Intravenous Methylprednisolone Administration in Rats: Acute Versus Chronic Dosing

Pharmacological time-series data, from comparative dosing studies, are critical to characterizing drug effects. Reconciling the data from multiple studies is inevitably difficult; multiple in vivo high-throughput -omics studies are necessary to capture the global and temporal effects of the drug, but these experiments, though analogous, differ in (microarray or other) platforms, time-scales, and dosing regimens and thus cannot be directly combined or compared. This investigation addresses this reconciliation issue with a meta-analysis technique aimed at assessing the intrinsic activity at the pathway level. The purpose of this is to characterize the dosing effects of methylprednisolone (MPL), a widely used anti-inflammatory and immunosuppressive corticosteroid (CS), within the liver. A multivariate decomposition approach is applied to analyze acute and chronic MPL dosing in male adrenalectomized rats and characterize the dosing-dependent differences in the dynamic response of MPL-responsive signaling and metabolic pathways. We demonstrate how to deconstruct signaling and metabolic pathways into their constituent pathway activities, activities which are scored for intrinsic pathway activity. Dosing-induced changes in the dynamics of pathway activities are compared using a model-based assessment of pathway dynamics, extending the principles of pharmacokinetics/pharmacodynamics (PKPD) to describe pathway activities. The model-based approach enabled us to hypothesize on the likely emergence (or disappearance) of indirect dosing-dependent regulatory interactions, pointing to likely mechanistic implications of dosing of MPL transcriptional regulation. Both acute and chronic MPL administration induced a strong core of activity within pathway families including the following: lipid metabolism, amino acid metabolism, carbohydrate metabolism, metabolism of cofactors and vitamins, regulation of essential organelles, and xenobiotic metabolism pathway families. Pathway activities alter between acute and chronic dosing, indicating that MPL response is dosing dependent. Furthermore, because multiple pathway activities are dominant within a single pathway, we observe that pathways cannot be defined by a single response. Instead, pathways are defined by multiple, complex, and temporally related activities corresponding to different subgroups of genes within each pathway.

Functional proteomic analysis of corticosteroid pharmacodynamics in rat liver: Relationship to hepatic stress, signaling, energy regulation, and drug metabolism

Corticosteroids (CS) are anti-inflammatory agents that cause extensive pharmacogenomic and proteomic changes in multiple tissues. An understanding of the proteome-wide effects of CS in liver and its relationships to altered hepatic and systemic physiology remains incomplete. Here, we report the application of a functional pharmacoproteomic approach to gain integrated insight into the complex nature of CS responses in liver in vivo. An in-depth functional analysis was performed using rich pharmacodynamic (temporal-based) proteomic data measured over 66h in rat liver following a single dose of methylprednisolone (MPL). Data mining identified 451 differentially regulated proteins. These proteins were analyzed on the basis of temporal regulation, cellular localization, and literature-mined functional information. Of the 451 proteins, 378 were clustered into six functional groups based on major clinically-relevant effects of CS in liver. MPL-responsive proteins were highly localized in the mitochondria (20%) and cytosol (24%). Interestingly, several proteins were related to hepatic stress and signaling processes, which appear to be involved in secondary signaling cascades and in protecting the liver from CS-induced oxidative damage. Consistent with known adverse metabolic effects of CS, several rate-controlling enzymes involved in amino acid metabolism, gluconeogenesis, and fatty-acid metabolism were altered by MPL. In addition, proteins involved in the metabolism of endogenous compounds, xenobiotics, and therapeutic drugs including cytochrome P450 and Phase-II enzymes were differentially regulated. Proteins related to the inflammatory acute-phase response were up-regulated in response to MPL. Functionally-similar proteins showed large diversity in their temporal profiles, indicating complex mechanisms of regulation by CS.

SIGNIFICANCE

Clinical use of corticosteroid (CS) therapy is frequent and chronic. However, current knowledge on the proteome-level effects of CS in liver and other tissues is sparse. While transcriptomic regulation following methylprednisolone (MPL) dosing has been temporally examined in rat liver, proteomic assessments are needed to better characterize the tissue-specific functional aspects of MPL actions. This study describes a functional pharmacoproteomic analysis of dynamic changes in MPL-regulated proteins in liver and provides biological insight into how steroid-induced perturbations on a molecular level may relate to both adverse and therapeutic responses presented clinically.

Highly multiplexed and reproducible ion-current-based strategy for large-scale quantitative proteomics and the application to protein expression dynamics induced by methylprednisolone in 60 rats

A proteome-level time-series study of drug effects (i.e., pharmacodynamics) is critical for understanding mechanisms of action and systems pharmacology, but is challenging, because of the requirement of a proteomics method for reliable quantification of many biological samples. Here, we describe a highly reproducible strategy, enabling a global, large-scale investigation of the expression dynamics of corticosteroid-regulated proteins in livers from adrenalectomized rats over 11 time points after drug dosing (0.5–66 h, N = 5/point). The analytical advances include (i) exhaustive tissue extraction with a Polytron/sonication procedure in a detergent cocktail buffer, and a cleanup/digestion procedure providing very consistent protein yields (relative standard deviation (RSD%) of 2.7%–6.4%) and peptide recoveries (4.1–9.0%) across the 60 animals; (ii) an ultrahigh-pressure nano-LC setup with substantially improved temperature stabilization, pump-noise suppression, and programmed interface cleaning, enabling excellent reproducibility for continuous analyses of numerous samples; (iii) separation on a 100-cm-long column (2-μm particles) with high reproducibility for days to enable both in-depth profiling and accurate peptide ion-current match; and (iv) well-controlled ion-current-based quantification. To obtain high-quality quantitative data necessary to describe the 11 time-points protein expression temporal profiles, strict criteria were used to define “quantifiable proteins”. A total of 323 drug-responsive proteins were revealed with confidence, and the time profiles of these proteins provided new insights into the diverse temporal changes of biological cascades associated with hepatic metabolism, response to hormone stimuli, gluconeogenesis, inflammatory responses, and protein translation processes. Most profile changes persisted well after the drug was eliminated. The developed strategy can also be broadly applied in preclinical and clinical research, where the analysis of numerous biological replicates is crucial.

Tissue-specific gene expression and regulation in liver and muscle following chronic corticosteroid administration

Although corticosteroids (CSs) affect gene expression in multiple tissues, the array of genes that are regulated by these catabolic steroids is diverse, highly tissue specific, and depends on their functions in the tissue. Liver has many important functions in performing and regulating diverse metabolic processes. Muscle, in addition to its mechanical role, is critical in maintaining systemic energy homeostasis and accounts for about 80% of insulin-directed glucose disposal. Consequently, a better understanding of CS pharmacogenomic effects in these tissues would provide valuable information regarding the tissue-specificity of transcriptional dynamics, and would provide insights into the underlying molecular mechanisms of action for both beneficial and detrimental effects. We performed an integrated analysis of transcriptional data from liver and muscle in response to methylprednisolone (MPL) infusion, which included clustering and functional annotation of clustered gene groups, promoter extraction and putative transcription factor (TF) identification, and finally, regulatory closeness (RC) identification. This analysis allowed the identification of critical transcriptional responses and CS-responsive functions in liver and muscle during chronic MPL administration, the prediction of putative transcriptional regulators relevant to transcriptional responses of CS-affected genes which are also potential secondary bio-signals altering expression levels of target-genes, and the exploration of the tissue-specificity and biological significance of gene expression patterns, CS-responsive functions, and transcriptional regulation. The analysis provided an integrated description of the genomic and functional effects of chronic MPL infusion in liver and muscle.

Glucocorticoid effects on adiponectin expression

Maintenance of energy metabolism and glucose homeostasis is achieved by the regulatory effects of many hormones and their interactions. Glucocorticoids produced from adrenal cortex and adiponectin produced by adipose tissue play important roles in the production, distribution, storage, and utilization of energy substrates. Glucocorticoids are involved in the activation of a number of catabolic processes by affecting the expression of a plethora of genes, while adiponectin acts primarily as an insulin sensitizer. Both are regulated by a number of physiological and pharmacological factors. Although the effects of glucocorticoids on adiponectin expression have been extensively studied in different in vitro, animal and clinical study settings, no consensus has been reached. This report reviews the primary literature concerning the effects of glucocorticoids on adiponectin expression and identifies potential reasons for the contradictory results between different studies. In addition, methods to gain better insights pertaining to the regulation of adiponectin expression are discussed.

Comparative analysis of acute and chronic corticosteroid pharmacogenomic effects in rat liver: transcriptional dynamics and regulatory structures

Background

Comprehensively understanding corticosteroid pharmacogenomic effects is an essential step towards an insight into the underlying molecular mechanisms for both beneficial and detrimental clinical effects. Nevertheless, even in a single tissue different methods of corticosteroid administration can induce different patterns of expression and regulatory control structures. Therefore, rich in vivo datasets of pharmacological time-series with two dosing regimens sampled from rat liver are examined for temporal patterns of changes in gene expression and their regulatory commonalities.

Results

The study addresses two issues, including (1) identifying significant transcriptional modules coupled with dynamic expression patterns and (2) predicting relevant common transcriptional controls to better understand the underlying mechanisms of corticosteroid adverse effects. Following the orientation of meta-analysis, an extended computational approach that explores the concept of agreement matrix from consensus clustering has been proposed with the aims of identifying gene clusters that share common expression patterns across multiple dosing regimens as well as handling challenges in the analysis of microarray data from heterogeneous sources, e.g. different platforms and time-grids in this study. Six significant transcriptional modules coupled with typical patterns of expression have been identified. Functional analysis reveals that virtually all enriched functions (gene ontologies, pathways) in these modules are shown to be related to metabolic processes, implying the importance of these modules in adverse effects under the administration of corticosteroids. Relevant putative transcriptional regulators (e.g. RXRF, FKHD, SP1F) are also predicted to provide another source of information towards better understanding the complexities of expression patterns and the underlying regulatory mechanisms of those modules.

Conclusions

We have proposed a framework to identify significant coexpressed clusters of genes across multiple conditions experimented from different microarray platforms, time-grids, and also tissues if applicable. Analysis on rich in vivo datasets of corticosteroid time-series yielded significant insights into the pharmacogenomic effects of corticosteroids, especially the relevance to metabolic side-effects. This has been illustrated through enriched metabolic functions in those transcriptional modules and the presence of GRE binding motifs in those enriched pathways, providing significant modules for further analysis on pharmacogenomic corticosteroid effects.

Pharmacodynamic/pharmacogenomic modeling of insulin resistance genes in rat muscle after methylprednisolone treatment: exploring regulatory signaling cascades

Corticosteroids (CS) effects on insulin resistance related genes in rat skeletal muscle were studied. In our acute study, adrenalectomized (ADX) rats were given single doses of 50 mg/kg methylprednisolone (MPL) intravenously. In our chronic study, ADX rats were implanted with Alzet mini-pumps giving zero-order release rates of 0.3 mg/kg/h MPL and sacrificed at various times up to 7 days. Total RNA was extracted from gastrocnemius muscles and hybridized to Affymetrix GeneChips. Data mining and literature searches identified 6 insulin resistance related genes which exhibited complex regulatory pathways. Insulin receptor substrate-1 (IRS-1), uncoupling protein 3 (UCP3), pyruvate dehydrogenase kinase isoenzyme 4 (PDK4), fatty acid translocase (FAT) and glycerol-3-phosphate acyltransferase (GPAT) dynamic profiles were modeled with mutual effects by calculated nuclear drug-receptor complex (DR(N)) and transcription factors. The oscillatory feature of endothelin-1 (ET-1) expression was depicted by a negative feedback loop. These integrated models provide testable quantitative hypotheses for these regulatory cascades.

Mathematical modeling of corticosteroid pharmacogenomics in rat muscle following acute and chronic methylprednisolone dosing

The pharmacogenomic effects of a corticosteroid (CS) were assessed in rat skeletal muscle using microarrays. Adrenalectomized (ADX) rats were treated with methylprednisolone (MPL) by either 50 mg/kg intravenous injection or 7-day 0.3 mg/kg/h infusion through subcutaneously implanted pumps. RNAs extracted from individual rat muscles were hybridized to Affymetrix Rat Genome Genechips. Data mining yielded 653 and 2316 CS-responsive probe sets following MPL bolus and infusion treatments. Of these, 196 genes were controlled by MPL under both dosing conditions. Cluster analysis revealed that 124 probe sets exhibited three typical expression dynamic profiles following acute dosing. Cluster A consisted of up-regulated probe sets which were grouped into five subclusters each exhibiting unique temporal patterns during the infusion. Cluster B comprised down-regulated probe sets which were divided into two subclusters with distinct dynamics during the infusion. Cluster C probe sets exhibited delayed down-regulation under both bolus and infusion conditions. Among those, 104 probe sets were further grouped into subclusters based on their profiles following chronic MPL dosing. Several mathematical models were proposed and adequately captured the temporal patterns for each subcluster. Multiple types of dosing regimens are needed to resolve common determinants of gene regulation as chronic exposure results in unexpected differences in gene expression compared to acute dosing. Pharmacokinetic/pharmacodynamic (PK/PD) modeling provides a quantitative tool for elucidating the complexities of CS pharmacogenomics in skeletal muscle.

Pharmacodynamic modeling of acute and chronic effects of methylprednisolone on hepatic urea cycle genes in rats

Corticosteroids (CS) regulate many enzymes at both mRNA and protein levels. This study used microarrays to broadly assess regulation of various genes related to the greater urea cycle and employs pharmacokinetic/pharmacodynamic (PK/PD) modeling to quantitatively analyze and compare the temporal profiles of these genes during acute and chronic exposure to methylprednisolone (MPL). One group of adrenalectomized male Wistar rats received an intravenous bolus dose (50 mg/kg) of MPL, whereas a second group received MPL by a subcutaneous infusion (Alzet osmotic pumps) at a rate of 0.3 mg/kg/hr for seven days. The rats were sacrificed at various time points over 72 hours (acute) or 168 hours (chronic) and livers were harvested. Total RNA was extracted and Affymetrix gene chips (RG_U34A for acute and RAE 230A for chronic) were used to identify genes regulated by CS. Besides five primary urea cycle enzymes, many other genes related to the urea cycle showed substantial changes in mRNA expression. Some genes that were simply up- or down-regulated after acute MPL showed complex biphasic patterns upon chronic infusion indicating involvement of secondary regulation. For the simplest patterns, indirect response models were used to describe the nuclear steroid-bound receptor mediated increase or decrease in gene transcription (e.g. tyrosine aminotransferase, glucocorticoid receptor). For the biphasic profiles, involvement of a secondary biosignal was assumed (e.g. ornithine decarboxylase, CCAAT/enhancer binding protein) and more complex models were derived. Microarrays were used successfully to explore CS effects on various urea cycle enzyme genes. PD models presented in this report describe testable hypotheses regarding molecular mechanisms and quantitatively characterize the direct or indirect regulation of various genes by CS.

Pharmacodynamic Modeling of Acute and Chronic Effects of Methylprednisolone on Hepatic Urea Cycle Genes in Rats

Corticosteroids (CS) regulate many enzymes at both mRNA and protein levels. This study used microarrays to broadly assess regulation of various genes related to the greater urea cycle and employs pharmacokinetic/pharmacodynamic (PK/PD) modeling to quantitatively analyze and compare the temporal profiles of these genes during acute and chronic exposure to methylprednisolone (MPL). One group of adrenalectomized male Wistar rats received an intravenous bolus dose (50 mg/kg) of MPL, whereas a second group received MPL by a subcutaneous infusion (Alzet osmotic pumps) at a rate of 0.3 mg/kg/hr for seven days. The rats were sacrificed at various time points over 72 hours (acute) or 168 hours (chronic) and livers were harvested. Total RNA was extracted and Affymetrix® gene chips (RGU34A for acute and RAE 230A for chronic) were used to identify genes regulated by CS. Besides five primary urea cycle enzymes, many other genes related to the urea cycle showed substantial changes in mRNA expression. Some genes that were simply up- or down-regulated after acute MPL showed complex biphasic patterns upon chronic infusion indicating involvement of secondary regulation. For the simplest patterns, indirect response models were used to describe the nuclear steroid-bound receptor mediated increase or decrease in gene transcription (e.g. tyrosine aminotransferase, glucocorticoid receptor). For the biphasic profiles, involvement of a secondary biosignal was assumed (e.g. ornithine decarboxylase, CCAAT/enhancer binding protein) and more complex models were derived. Microarrays were used successfully to explore CS effects on various urea cycle enzyme genes. PD models presented in this report describe testable hypotheses regarding molecular mechanisms and quantitatively characterize the direct or indirect regulation of various genes by CS.

Bayesian state space models for inferring and predicting temporal gene expression profiles

Prediction of gene dynamic behavior is a challenging and important problem in genomic research while estimating the temporal correlations and non-stationarity are the keys in this process. Unfortunately, most existing techniques used for the inclusion of the temporal correlations treat the time course as evenly distributed time intervals and use stationary models with time-invariant settings. This is an assumption that is often violated in microarray time course data since the time course expression data are at unequal time points, where the difference in sampling times varies from minutes to days. Furthermore, the unevenly spaced short time courses with sudden changes make the prediction of genetic dynamics difficult. In this paper, we develop two types of Bayesian state space models to tackle this challenge for inferring and predicting the gene expression profiles associated with diseases. In the univariate time-varying Bayesian state space models we treat both the stochastic transition matrix and the observation matrix time-variant with linear setting and point out that this can easily be extended to nonlinear setting. In the multivariate Bayesian state space model we include temporal correlation structures in the covariance matrix estimations. In both models, the unevenly spaced short time courses with unseen time points are treated as hidden state variables. Bayesian approaches with various prior and hyper-prior models with MCMC algorithms are used to estimate the model parameters and hidden variables. We apply our models to multiple tissue polygenetic affymetrix data sets. Results show that the predictions of the genomic dynamic behavior can be well captured by the proposed models.

Analysis of time-series gene expression data: methods, challenges, and opportunities

Monitoring the change in expression patterns over time provides the distinct possibility of unraveling the mechanistic drivers characterizing cellular responses. Gene arrays measuring the level of mRNA expression of thousands of genes simultaneously provide a method of high-throughput data collection necessary for obtaining the scope of data required for understanding the complexities of living organisms. Unraveling the coherent complex structures of transcriptional dynamics is the goal of a large family of computational methods aiming at upgrading the information content of time-course gene expression data. In this review, we summarize the qualitative characteristics of these approaches, discuss the main challenges that this type of complex data present, and, finally, explore the opportunities in the context of developing mechanistic models of cellular response.

Microarray analysis of the temporal response of skeletal muscle to methylprednisolone: comparative analysis of two dosing regimens

The transcriptional response of skeletal muscle to chronic corticosteroid exposure was examined over 168 h and compared with the response profiles observed following a single dose of corticosteroid. Male adrenalectomized Wistar rats were given a constant-rate infusion of 0.3 mg x kg(-1) x h(-1) methylprednisolone for up to 7 days via subcutaneously implanted minipumps. Four control and forty drug-treated animals were killed at ten different time points during infusion. Liver total RNAs were hybridized to 44 individual Affymetrix REA230A gene chips. Previously, we described a filtration approach for identifying genes of interest in microarray data sets developed from tissues of rats treated with methylprednisolone (MPL) following acute dosing. Here, a similar approach involving a series of three filters was applied sequentially to identify genes of interest. These filters were designed to eliminate probe sets that were not expressed in the tissue, not regulated by the drug, or did not meet defined quality control standards. Filtering eliminated 86% of probe sets, leaving a remainder of 2,316 for further consideration. In a previous study, 653 probe sets were identified as MPL regulated following administration of a single (acute) dose of the drug. Comparison of the two data sets yielded 196 genes identified as regulated by MPL in both dosing regimens. Because of receptor downregulation, it was predicted that genes regulated by receptor-glucocorticoid response element interactions would exhibit tolerance in chronic profiles. However, many genes did not exhibit steroid tolerance, indicating that present perspectives on the mechanism of glucocorticoid action cannot entirely explain all temporal profiles.

A microarray analysis of the temporal response of liver to methylprednisolone: a comparative analysis of two dosing regimens

Microarray analyses were performed on livers from adrenalectomized male Wistar rats chronically infused with methylprednisolone (MPL) (0.3 mg/kg.h) using Alzet mini-osmotic pumps for periods ranging from 6 h to 7 d. Four control and 40 drug-treated animals were killed at 10 different times during drug infusion. Total RNA preparations from the livers of these animals were hybridized to 44 individual Affymetrix REA230A gene chips, generating data for 15,967 different probe sets for each chip. A series of three filters were applied sequentially. These filters were designed to eliminate probe sets that were not expressed in the tissue, were not regulated by the drug, or did not meet defined quality control standards. These filters eliminated 13,978 probe sets (87.5%) leaving a remainder of 1989 probe sets for further consideration. We previously described a similar dataset obtained from animals after administration of a single dose of MPL (50 mg/kg given iv). That study involved 16 time points over a 72-h period. A similar filtering schema applied to the single-bolus-dose dataset identified 1519 probe sets as being regulated by MPL. A comparison of datasets from the two different dosing regimens identified 358 genes that were regulated by MPL in response to both dosing regimens. Regulated genes were grouped into 13 categories, mainly on gene product function. The temporal profiles of these common genes were subjected to detailed scrutiny. Examination of temporal profiles demonstrates that current perspectives on the mechanism of glucocorticoid action cannot entirely explain the temporal profiles of these regulated genes.

Application of scaling factors in simultaneous modeling of microarray data from diverse chips

PURPOSE

Microarrays have been utilized in many biological, physiological and pharmacological studies as a high-throughput genomic technique. Several generations of Affymetrix GeneChip microarrays are widely used in gene expression studies. However, differences in intensities of signals for different probe sets that represent the same gene on various types of Affymetrix chips make comparison of datasets complicated.

MATERIALS AND METHODS

A power coefficient scaling factor was applied in the pharmacokinetic/pharmacodynamic (PK/PD) modeling to account for differences in probe set sensitivities (i.e., signal intensities). Microarray data from muscle and liver following methylprednisolone 50 mg/kg i.v. bolus and 0.3 mg/kg/h infusion regimens were taken as an exemplar.

RESULTS

The scaling factor applied to the pharmacodynamic output function was used to solve the problem of intensity differences between probe sets. This approach yielded consistent pharmacodynamic parameters for the applied models.

CONCLUSIONS

Modeling of pharmacodynamic/pharmacogenomic (PD/PG) data from diverse chips should be performed with caution due to differential probe set intensities. In such circumstances, a power scaling factor can be applied in the modeling.

Temporal changes in gene expression in rainbow trout exposed to ethynyl estradiol

We examined changes in the genomic response during continuous exposure to the xenoestrogen ethynyl estradiol. Isogenic rainbow trout Oncorhynchus mykiss were exposed to nominal concentrations of 100 ng/L ethynyl estradiol (EE2) for a period of 3 weeks. At fixed time points within the exposure, fish were euthanized, livers harvested and RNA extracted. Fluorescently labeled cDNA were generated and hybridized against a commercially available Salmonid array (GRASP project, University of Victoria, Canada) spotted with 16,000 cDNAs. The slides were scanned to measure abundance of a given transcript in each sample relative to controls. Data were analyzed via Genespring (Silicon Genetics) to identify a list of up and down regulated genes, and to determine gene clustering patterns that can be used as "expression signatures". Gene ontology was determined using the annotation available from the GRASP website. Our analysis indicates each exposure time period generated specific gene expression profiles. Changes in gene expression were best understood by grouping genes by their gene expression profiles rather than examining fold change at a particular time point. Many of the genes commonly used as biomarkers of exposure to xenoestrogens were not induced initially and did not have gene expression profiles typical of the majority of genes with altered expression.

Bioinformatics analysis of the early inflammatory response in a rat thermal injury model

Background

Thermal injury is among the most severe forms of trauma and its effects are both local and systemic. Response to thermal injury includes cellular protection mechanisms, inflammation, hypermetabolism, prolonged catabolism, organ dysfunction and immuno-suppression. It has been hypothesized that gene expression patterns in the liver will change with severe burns, thus reflecting the role the liver plays in the response to burn injury. Characterizing the molecular fingerprint (i.e., expression profile) of the inflammatory response resulting from burns may help elucidate the activated mechanisms and suggest new therapeutic intervention. In this paper we propose a novel integrated framework for analyzing time-series transcriptional data, with emphasis on the burn-induced response within the context of the rat animal model. Our analysis robustly identifies critical expression motifs, indicative of the dynamic evolution of the inflammatory response and we further propose a putative reconstruction of the associated transcription factor activities.

Results

Implementation of our algorithm on data obtained from an animal (rat) burn injury study identified 281 genes corresponding to 4 unique profiles. Enrichment evaluation upon both gene ontologies and transcription factors, verifies the inflammation-specific character of the selections and the rationalization of the burn-induced inflammatory response. Conducting the transcription network reconstruction and analysis, we have identified transcription factors, including AHR, Octamer Binding Proteins, Kruppel-like Factors, and cell cycle regulators as being highly important to an organism's response to burn response. These transcription factors are notable due to their roles in pathways that play a part in the gross physiological response to burn such as changes in the immune response and inflammation.

Conclusion

Our results indicate that our novel selection/classification algorithm has been successful in selecting out genes with play an important role in thermal injury. Additionally, we have demonstrated the value of an integrative approach in identifying possible points of intervention, namely the activation of certain transcription factors that govern the organism's response.

Corticosteroid-regulated genes in rat kidney: mining time series array data

Kidney is a major target for adverse effects associated with corticosteroids. A microarray dataset was generated to examine changes in gene expression in rat kidney in response to methylprednisolone. Four control and 48 drug-treated animals were killed at 16 times after drug administration. Kidney RNA was used to query 52 individual Affymetrix chips, generating data for 15,967 different probe sets for each chip. Mining techniques applicable to time series data that identify drug-regulated changes in gene expression were applied. Four sequential filters eliminated probe sets that were not expressed in the tissue, not regulated by drug, or did not meet defined quality control standards. These filters eliminated 14,890 probe sets (94%) from further consideration. Application of judiciously chosen filters is an effective tool for data mining of time series datasets. The remaining data can then be further analyzed by clustering and mathematical modeling. Initial analysis of this filtered dataset identified a group of genes whose pattern of regulation was highly correlated with prototype corticosteroid enhanced genes. Twenty genes in this group, as well as selected genes exhibiting either downregulation or no regulation, were analyzed for 5' GRE half-sites conserved across species. In general, the results support the hypothesis that the existence of conserved DNA binding sites can serve as an important adjunct to purely analytic approaches to clustering genes into groups with common mechanisms of regulation. This dataset, as well as similar datasets on liver and muscle, are available online in a format amenable to further analysis by others.

Temporal profiling of the transcriptional basis for the development of corticosteroid-induced insulin resistance in rat muscle

Elevated systemic levels of glucocorticoids are causally related to peripheral insulin resistance. The pharmacological use of synthetic glucocorticoids (corticosteroids) often results in insulin resistance/type II diabetes. Skeletal muscle is responsible for close to 80% of the insulin-induced systemic disposal of glucose and is a major target for glucocorticoid-induced insulin resistance. We used Affymetrix gene chips to profile the dynamic changes in mRNA expression in rat skeletal muscle in response to a single bolus dose of the synthetic glucocorticoid methyl-prednisolone. Temporal expression profiles (analyzed on individual chips) were obtained from tissues of 48 drug-treated animals encompassing 16 time points over 72 h following drug administration along with four vehicle-treated controls. Data mining identified 653 regulated probe sets out of 8799 present on the chip. Of these 653 probe sets we identified 29, which represented 22 gene transcripts, that were associated with the development of insulin resistance. These 29 probe sets were regulated in three fundamental temporal patterns. 16 probe sets coding for 12 different genes had a profile of enhanced expression. 10 probe sets coding for eight different genes showed decreased expression and three probe sets coding for two genes showed biphasic temporal signatures. These transcripts were grouped into four general functional categories: signal transduction, transcription regulation, carbohydrate/fat metabolism, and regulation of blood flow to the muscle. The results demonstrate the polygenic nature of transcriptional changes associated with insulin resistance that can provide a temporal scaffolding for translational and post-translational data as they become available.