Importance of replication in analyzing time-series gene expression data: Corticosteroid dynamics and circadian patterns in rat liver

<p>Abstract</p> <p>Background</p> <p>Microarray technology is a powerful and widely accepted experimental technique in molecular biology that allows studying genome wide transcriptional responses. However, experimental data usually contain potential sources of uncertainty and thus many experiments are now designed with repeated measurements to better assess such inherent variability. Many computational methods have been proposed to account for the variability in replicates. As yet, there is no model to output expression profiles accounting for replicate information so that a variety of computational models that take the expression profiles as the input data can explore this information without any modification.</p> <p>Results</p> <p>We propose a methodology which integrates replicate variability into expression profiles, to generate so-called 'true' expression profiles. The study addresses two issues: (i) develop a statistical model that can estimate 'true' expression profiles which are more robust than the average profile, and (ii) extend our previous micro-clustering which was designed specifically for clustering time-series expression data. The model utilizes a previously proposed error model and the concept of 'relative difference'. The clustering effectiveness is demonstrated through synthetic data where several methods are compared. We subsequently analyze <it>in vivo </it>rat data to elucidate circadian transcriptional dynamics as well as liver-specific corticosteroid induced changes in gene expression.</p> <p>Conclusions</p> <p>We have proposed a model which integrates the error information from repeated measurements into the expression profiles. Through numerous synthetic and real time-series data, we demonstrated the ability of the approach to improve the clustering performance and assist in the identification and selection of informative expression motifs.</p

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

<p>Abstract</p> <p>Background</p> <p>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 <it>in vivo </it>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.</p> <p>Results</p> <p>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.</p> <p>Conclusions</p> <p>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 <it>in vivo </it>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.</p

Circadian signatures in rat liver: from gene expression to pathways

<p>Abstract</p> <p>Background</p> <p>Circadian rhythms are 24 hour oscillations in many behavioural, physiological, cellular and molecular processes that are controlled by an endogenous clock which is entrained to environmental factors including light, food and stress. Transcriptional analyses of circadian patterns demonstrate that genes showing circadian rhythms are part of a wide variety of biological pathways.</p> <p>Pathway activity method can identify the significant pattern of the gene expression levels within a pathway. In this method, the overall gene expression levels are translated to a reduced form, pathway activity levels, via singular value decomposition (SVD). A given pathway represented by pathway activity levels can then be as analyzed using the same approaches used for analyzing gene expression levels. We propose to use pathway activity method across time to identify underlying circadian pattern of pathways.</p> <p>Results</p> <p>We used synthetic data to demonstrate that pathway activity analysis can evaluate the underlying circadian pattern within a pathway even when circadian patterns cannot be captured by the individual gene expression levels. In addition, we illustrated that pathway activity formulation should be coupled with a significance analysis to distinguish biologically significant information from random deviations. Next, we performed pathway activity level analysis on a rich time series of transcriptional profiling in rat liver. The over-represented five specific patterns of pathway activity levels, which cannot be explained by random event, exhibited circadian rhythms. The identification of the circadian signatures at the pathway level identified 78 pathways related to energy metabolism, amino acid metabolism, lipid metabolism and DNA replication and protein synthesis, which are biologically relevant in rat liver. Further, we observed tight coordination between cholesterol biosynthesis and bile acid biosynthesis as well as between folate biosynthesis, one carbon pool by folate and purine-pyrimidine metabolism. These coupled pathways are parts of a sequential reaction series where the product of one pathway is the substrate of another pathway.</p> <p>Conclusions</p> <p>Rather than assessing the importance of a single gene beforehand and map these genes onto pathways, we instead examined the orchestrated change within a pathway. Pathway activity level analysis could reveal the underlying circadian dynamics in the microarray data with an unsupervised approach and biologically relevant results were obtained.</p

Identification of Global Transcriptional Dynamics

One of the challenges in exploiting high throughput measurement techniques such as microarrays is the conversion of the vast amounts of data obtained into relevant knowledge. Of particular importance is the identification of the intrinsic response of a transcriptional experiment and the characterization of the underlying dynamics.The proposed algorithm seeks to provide the researcher a summary as to various aspects relating to the dynamic progression of a biological system, rather than that of individual genes. The approach is based on the identification of smaller number of expression motifs that define the transcriptional state of the system which quantifies the deviation of the cellular response from a control state in the presence of an external perturbation. The approach is demonstrated with a number of data sets including a synthetic base case and four animal studies. The synthetic dataset will be used to establish the response of the algorithm on a "null" dataset, whereas the four different experimental datasets represent a spectrum of possible time course experiments in terms of the degree of perturbation associated with the experiment as well as representing a wide range of temporal sampling strategies. This wide range of experimental datasets will thus allow us to explore the performance of the proposed algorithm and determine its ability identify relevant information.In this work, we present a computational approach which operates on high throughput temporal gene expression data to assess the information content of the experiment, identify dynamic markers of important processes associated with the experimental perturbation, and summarize in a concise manner the evolution of the system over time with respect to the experimental perturbation

Modeling of corticosteroid pharmacogenomics in rat liver using gene microarrays.&quot;

ABSTRACT Corticosteroid (CS) pharmacogenomics was studied using gene microarrays in rat liver. Methylprednisolone (MPL) was administered intravenously at 50 mg/kg. Rats were sacrificed and liver excised at 17 time points over 72 h. RNAs from individual livers were used to query Affymetrix GeneChips that contain sequences for 8000 genes. Cluster analysis revealed six temporal patterns consisting of 197 CS-responsive probes representing 143 genes. Based on our fifth-generation model of CS pharmacokinetics/pharmacodynamics (PK/PD), mechanistic models were developed to describe the time pattern for each CS-responsive gene. Two clusters showed increased expression with different effect duration. PK/PD models assuming CS stimulation of mRNA synthesis were applied. Another two clusters showed an initial decline followed by delayed increase, suggesting two mechanisms might be involved jointly. The initial suppression was captured by CS inhibition of mRNA synthesis or stimulation of degradation. CS may also stimulate the production of a biosignal (transcription factors or other hormones), which can cause secondary induction of the target mRNA. One cluster showed a very abrupt increase in message followed by rapid decrease. These genes were lymphocytic in origin and were modeled combining the fast gene induction effect of CS in lymphoid cells and its direct lymphocyte trafficking effect. Another cluster showed reduction persisting for 18 h, which was described by CS inhibition of mRNA synthesis. Our results reveal the marked diversity of genes regulated by CS via a limited array of mechanisms. These PK/PD models provide quantitation of CS pharmacogenomics and new hypotheses regarding understanding of diverse mechanisms of CS receptor-gene mediated action

Modeling Corticosteroid Effects in a Rat Model of Rheumatoid Arthritis I: Mechanistic Disease Progression Model for the Time Course of Collagen-Induced Arthritis in Lewis Rats

ABSTRACT A mechanism-based model was developed to describe the time course of arthritis progression in the rat. Arthritis was induced in male Lewis rats with type II porcine collagen into the base of the tail. Disease progression was monitored by paw swelling, bone mineral density (BMD), body weights, plasma corticosterone (CST) concentrations, and tumor necrosis factor (TNF)-␣, interleukin (IL)-1␤, IL-6, and glucocorticoid receptor (GR) mRNA expression in paw tissue. Bone mineral density was determined by PIXImus II dual energy X-ray densitometry. Plasma CST was assayed by high-performance liquid chromatography. Cytokine and GR mRNA were determined by quantitative real-time polymerase chain reaction. Disease progression models were constructed from transduction and indirect response models and applied using S-ADAPT software. A delay in the onset of increased paw TNF-␣ and IL-6 mRNA concentrations was successfully characterized by simple transduction. This rise was closely followed by an up-regulation of GR mRNA and CST concentrations. Paw swelling and body weight responses peaked approximately 21 days after induction, whereas bone mineral density changes were greatest at 23 days after induction. After peak response, the time course in IL-1␤, IL-6 mRNA, and paw edema slowly declined toward a disease steady state. Model parameters indicate TNF-␣ and IL-1␤ mRNA most significantly induce paw edema, whereas IL-6 mRNA exerted the most influence on BMD. The model for bone mineral density captures rates of turnover of cancellous and cortical bone and the fraction of each in the different regions analyzed. This small systems model integrates and quantitates multiple factors contributing to arthritis in rats

Adipose Tissue Deficiency and Chronic Inflammation in Diabetic Goto-Kakizaki Rats

Type 2 diabetes (T2DM) is a heterogeneous group of diseases that is progressive and involves multiple tissues. Goto-Kakizaki (GK) rats are a polygenic model with elevated blood glucose, peripheral insulin resistance, a non-obese phenotype, and exhibit many degenerative changes observed in human T2DM. As part of a systems analysis of disease progression in this animal model, this study characterized the contribution of adipose tissue to pathophysiology of the disease. We sacrificed subgroups of GK rats and appropriate controls at 4, 8, 12, 16 and 20 weeks of age and carried out a gene array analysis of white adipose tissue. We expanded our physiological analysis of the animals that accompanied our initial gene array study on the livers from these animals. The expanded analysis included adipose tissue weights, HbA1c, additional hormonal profiles, lipid profiles, differential blood cell counts, and food consumption. HbA1c progressively increased in the GK animals. Altered corticosterone, leptin, and adiponectin profiles were also documented in GK animals. Gene array analysis identified 412 genes that were differentially expressed in adipose tissue of GKs relative to controls. The GK animals exhibited an age-specific failure to accumulate body fat despite their relatively higher calorie consumption which was well supported by the altered expression of genes involved in adipogenesis and lipogenesis in the white adipose tissue of these animals, including Fasn, Acly, Kklf9, and Stat3. Systemic inflammation was reflected by chronically elevated white blood cell counts. Furthermore, chronic inflammation in adipose tissue was evident from the differential expression of genes involved in inflammatory responses and activation of natural immunity, including two interferon regulated genes, Ifit and Iipg, as well as MHC class II genes. This study demonstrates an age specific failure to accumulate adipose tissue in the GK rat and the presence of chronic inflammation in adipose tissue from these animals

Effects of early feeding on growth velocity and overweight/obesity in a cohort of HIV unexposed South African infants and children

BACKGROUND: South Africa has the highest prevalence of overweight/obesity in Sub-Saharan Africa. Assessing the effect of modifiable factors such as early infant feeding on growth velocity and overweight/obesity is therefore important. This paper aimed to assess the effect of infant feeding in the transitional period (12 weeks) on 12–24 week growth velocity amongst HIV unexposed children using WHO growth velocity standards and on the age and sex adjusted body mass index (BMI) Z-score distribution at 2 years. METHODS: Data were from 3 sites in South Africa participating in the PROMISE-EBF trial. We calculated growth velocity Z-scores using the WHO growth standards and assessed feeding practices using 24-hour and 7-day recall data. We used quantile regression to study the associations between 12 week infant feeding and 12–24 week weight velocity (WVZ) with BMI-for-age Z-score at 2 years. We included the internal sample quantiles (70th and 90th centiles) that approximated the reference cut-offs of +2 (corresponding to overweight) and +3 (corresponding to obesity) of the 2 year BMI-for-age Z-scores. RESULTS: At the 2-year visit, 641 children were analysed (median age 22 months, IQR: 17–26 months). Thirty percent were overweight while 8.7% were obese. Children not breastfed at 12 weeks had higher 12–24 week mean WVZ and were more overweight and obese at 2 years. In the quantile regression, children not breastfed at 12 weeks had a 0.37 (95% CI 0.07, 0.66) increment in BMI-for-age Z-score at the 50th sample quantile compared to breast-fed children. This difference in BMI-for-age Z-score increased to 0.46 (95% CI 0.18, 0.74) at the 70th quantile and 0.68 (95% CI 0.41, 0.94) at the 90th quantile . The 12–24 week WVZ had a uniform independent effect across the same quantiles. CONCLUSIONS: This study demonstrates that the first 6 months of life is a critical period in the development of childhood overweight and obesity. Interventions targeted at modifiable factors such as early infant feeding practices may reduce the risks of rapid weight gain and subsequent childhood overweight/obesity.Scopu

Supplement: "Localization and broadband follow-up of the gravitational-wave transient GW150914" (2016, ApJL, 826, L13)

This Supplement provides supporting material for Abbott et al. (2016a). We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands