Circadian signatures in rat liver: from gene expression to pathways

Comparative Analysis of Circadian Transcriptomes Reveals Circadian Characteristics between Arabidopsis and Soybean

The circadian clock, an endogenous timing system, exists in nearly all organisms on Earth. The plant circadian clock has been found to be intricately linked with various essential biological activities. Extensive studies of the plant circadian clock have yielded valuable applications. However, the distinctions of circadian clocks in two important plant species, Arabidopsis thaliana and Glycine max (soybean), remain largely unexplored. This study endeavors to address this gap by conducting a comprehensive comparison of the circadian transcriptome profiles of Arabidopsis and soybean to uncover their distinct circadian characteristics. Utilizing non-linear regression fitting (COS) integrated with weights, we identified circadian rhythmic genes within both organisms. Through an in-depth exploration of circadian parameters, we unveiled notable differences between Arabidopsis and soybean. Furthermore, our analysis of core circadian clock genes shed light on the distinctions in central oscillators between these two species. Additionally, we observed that the homologous genes of Arabidopsis circadian clock genes in soybean exert a significant influence on the regulation of flowering and maturity of soybean. This phenomenon appears to stem from shifts in circadian parameters within soybean genes. These findings highlight contrasting biological activities under circadian regulation in Arabidopsis and soybean. This study not only underscores the distinctive attributes of these species, but also offers valuable insights for further scrutiny into the soybean circadian clock and its potential applications.

Ruminiclostridium 5, Parabacteroides distasonis, and bile acid profile are modulated by prebiotic diet and associate with facilitated sleep/clock realignment after chronic disruption of rhythms

Chronic disruption of rhythms (CDR) impacts sleep and can result in circadian misalignment of physiological systems which, in turn, is associated with increased disease risk. Exposure to repeated or severe stressors also disturbs sleep and diurnal rhythms. Prebiotic nutrients produce favorable changes in gut microbial ecology, the gut metabolome, and reduce several negative impacts of acute severe stressor exposure, including disturbed sleep, core body temperature rhythmicity, and gut microbial dysbiosis. In light of previous compelling evidence that prebiotic diet broadly reduces negative impacts of acute, severe stressors, we hypothesize that prebiotic diet will also effectively mitigate the negative impacts of chronic disruption of circadian rhythms on physiology and sleep/wake behavior. Male, Sprague Dawley rats were fed diets enriched in prebiotic substrates or calorically matched control chow. After 5 weeks on diet, rats were exposed to CDR (12 h light/dark reversal, weekly for 8 weeks) or remained on undisturbed normal light/dark cycles (NLD). Sleep EEG, core body temperature, and locomotor activity were recorded via biotelemetry in freely moving rats. Fecal samples were collected on experimental days -33, 0 (day of onset of CDR), and 42. Taxonomic identification and relative abundances of gut microbes were measured in fecal samples using 16S rRNA gene sequencing and shotgun metagenomics. Fecal primary, bacterially modified secondary, and conjugated bile acids were measured using liquid chromatography with tandem mass spectrometry (LC-MS/MS). Prebiotic diet produced rapid and stable increases in the relative abundances of Parabacteroides distasonis and Ruminiclostridium 5. Shotgun metagenomics analyses confirmed reliable increases in relative abundances of Parabacteroides distasonis and Clostridium leptum, a member of the Ruminiclostridium genus. Prebiotic diet also modified fecal bile acid profiles; and based on correlational and step-wise regression analyses, Parabacteroides distasonis and Ruminiclostridium 5 were positively associated with each other and negatively associated with secondary and conjugated bile acids. Prebiotic diet, but not CDR, impacted beta diversity. Measures of alpha diversity evenness were decreased by CDR and prebiotic diet prevented that effect. Rats exposed to CDR while eating prebiotic, compared to control diet, more quickly realigned NREM sleep and core body temperature (ClockLab) diurnal rhythms to the altered light/dark cycle. Finally, both cholic acid and Ruminiclostridium 5 prior to CDR were associated with time to realign CBT rhythms to the new light/dark cycle after CDR; whereas both Ruminiclostridium 5 and taurocholic acid prior to CDR were associated with NREM sleep recovery after CDR. These results support our hypothesis and suggest that ingestion of prebiotic substrates is an effective strategy to increase the relative abundance of health promoting microbes, alter the fecal bile acid profile, and facilitate the recovery and realignment of sleep and diurnal rhythms after circadian disruption.

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.

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.

Chronopharmacology of glucocorticoids

Glucocorticoids influence a wide array of metabolic, anti-inflammatory, immunosuppressive, and cognitive signaling processes, playing an important role in homeostasis and preservation of normal organ function. Synthesis is regulated by the hypothalamic-pituitary-adrenal (HPA) axis of which cortisol is the primary glucocorticoid in humans. Synthetic glucocorticoids are important pharmacological agents that augment the anti-inflammatory and immunosuppressive properties of endogenous cortisol and are widely used for the treatment of asthma, Crohn's disease, and rheumatoid arthritis, amongst other chronic conditions. The homeostatic activity of cortisol is disrupted by the administration of synthetic glucocorticoids and so there is interest in developing treatment options that minimize HPA axis disturbance while maintaining the pharmacological effects. Studies suggest that optimizing drug administration time can achieve this goal. The present review provides an overview of endogenous glucocorticoid activity and recent advances in treatment options that have further improved patient safety and efficacy with an emphasis on chronopharmacology.

Effects of 4-Hydroxy-2,3,3',4',5-Pentachlorobiphenyl (4-OH-CB107) on Liver Transcriptome in Rats: Implication in the Disruption of Circadian Rhythm and Fatty Acid Metabolism

Polychlorinated biphenyls (PCBs) and their hydroxylated metabolites (OH-PCBs) have been detected in tissues of both wild animals and humans. Several previous studies have suggested adverse effects of OH-PCBs on the endocrine and nervous systems in mammals. However, there have been no studies on transcriptome analysis of the effects of OH-PCBs, and thus, the whole picture and mechanisms underlying the adverse effects induced by OH-PCBs are still poorly understood. We therefore investigated the mRNA expression profile in the liver of adult male Wistar rats treated with 4-hydroxy-2,3,3',4',5-pentachlorobiphenyl (4-OH-CB107) to explore the genes responsive to OH-PCBs and to understand the potential effects of the chemical. Next-generation RNA sequencing analysis revealed changes in the expression of genes involved in the circadian rhythm and fatty acid metabolism, such as nuclear receptor subfamily 1, group D, member 1, aryl hydrocarbon receptor nuclear translocator-like protein 1, cryptochrome circadian clock 1, and enoyl-CoA hydratase and 3-hydroxyacyl-CoA dehydrogenase, in 4-OH-CB107-treated rats. In addition, biochemical analysis of the plasma revealed a dose-dependent increase in the leucine aminopeptidase, indicating the onset of liver damage. These results suggest that OH-PCB exposure may induce liver injury as well as disrupt the circadian rhythm and peroxisome proliferator-activated receptor-related fatty acid metabolism.

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.

The growing role of precision and personalized medicine for cancer treatment

Cancer is a devastating disease that takes the lives of hundreds of thousands of people every year. Due to disease heterogeneity, standard treatments, such as chemotherapy or radiation, are effective in only a subset of the patient population. Tumors can have different underlying genetic causes and may express different proteins in one patient versus another. This inherent variability of cancer lends itself to the growing field of precision and personalized medicine (PPM). There are many ongoing efforts to acquire PPM data in order to characterize molecular differences between tumors. Some PPM products are already available to link these differences to an effective drug. It is clear that PPM cancer treatments can result in immense patient benefits, and companies and regulatory agencies have begun to recognize this. However, broader changes to the healthcare and insurance systems must be addressed if PPM is to become part of standard cancer care.

Understanding Physiology in the Continuum: Integration of Information from Multiple -Omics Levels

In this paper, we discuss approaches for integrating biological information reflecting diverse physiologic levels. In particular, we explore statistical and model-based methods for integrating transcriptomic, proteomic and metabolomics data. Our case studies reflect responses to a systemic inflammatory stimulus and in response to an anti-inflammatory treatment. Our paper serves partly as a review of existing methods and partly as a means to demonstrate, using case studies related to human endotoxemia and response to methylprednisolone (MPL) treatment, how specific questions may require specific methods, thus emphasizing the non-uniqueness of the approaches. Finally, we explore novel ways for integrating -omics information with PKPD models, toward the development of more integrated pharmacology models.

Physiologically-based pharmacokinetic models: approaches for enabling personalized medicine

Personalized medicine strives to deliver the 'right drug at the right dose' by considering inter-person variability, one of the causes for therapeutic failure in specialized populations of patients. Physiologically-based pharmacokinetic (PBPK) modeling is a key tool in the advancement of personalized medicine to evaluate complex clinical scenarios, making use of physiological information as well as physicochemical data to simulate various physiological states to predict the distribution of pharmacokinetic responses. The increased dependency on PBPK models to address regulatory questions is aligned with the ability of PBPK models to minimize ethical and technical difficulties associated with pharmacokinetic and toxicology experiments for special patient populations. Subpopulation modeling can be achieved through an iterative and integrative approach using an adopt, adapt, develop, assess, amend, and deliver methodology. PBPK modeling has two valuable applications in personalized medicine: (1) determining the importance of certain subpopulations within a distribution of pharmacokinetic responses for a given drug formulation and (2) establishing the formulation design space needed to attain a targeted drug plasma concentration profile. This review article focuses on model development for physiological differences associated with sex (male vs. female), age (pediatric vs. young adults vs. elderly), disease state (healthy vs. unhealthy), and temporal variation (influence of biological rhythms), connecting them to drug product formulation development within the quality by design framework. Although PBPK modeling has come a long way, there is still a lengthy road before it can be fully accepted by pharmacologists, clinicians, and the broader industry.

Burn trauma disrupts circadian rhythms in rat liver

Circadian rhythms play an important role in maintaining homeostasis and solid organ function. The purpose of this study is to assess the implications of burn injury in rats on the underlying circadian patterns of gene expression in liver. Circadian-regulated genes and burn-induced genes were identified by applying consensus clustering methodology to temporally differentially expressed probe sets obtained from burn and sham-burn data sets. Of the liver specific genes which we hypothesize that exhibit circadian rhythmicity, 88% are not differentially expressed following the burn injury. Specifically, the vast majority of the circadian regulated-genes representing central carbon and nitrogen metabolism are "up-regulated" after the burn injury, indicating the onset of hypermetabolism. In addition, cell-cell junction and membrane structure related genes showing rhythmic behavior in the control group were not differentially expressed across time in the burn group, which could be an indication of hepatic damage due to the burn. Finally, the suppression of the immune function related genes is observed in the postburn phase, implying the severe "immunosuppression". Our results demonstrated that the short term response (24-h post injury) manifests a loss of circadian variability possibly compromising the host in terms of subsequent challenges.

Network biology in development of monoclonal antibody therapeutics

Monoclonal antibodies (mAbs) are large glycoproteins that recognize and remove/neutralize a specific target. Inflammation and inflammatory diseases are often treated with mAb-based therapeutics. Mathematical modeling is widely used in development of mAbs. Bioinformatics and structural modeling is used for humanization of mAbs and PK/PD modeling is extensively used in preclinical and clinical development. The objective of this commentary is to introduce systems biology-based modeling that can accelerate and improve development of mAbs.

A HIERARCHICAL BAYESIAN MODEL FOR INFERENCE OF COPY NUMBER VARIANTS AND THEIR ASSOCIATION TO GENE EXPRESSION

A number of statistical models have been successfully developed for the analysis of high-throughput data from a single source, but few methods are available for integrating data from different sources. Here we focus on integrating gene expression levels with comparative genomic hybridization (CGH) array measurements collected on the same subjects. We specify a measurement error model that relates the gene expression levels to latent copy number states which, in turn, are related to the observed surrogate CGH measurements via a hidden Markov model. We employ selection priors that exploit the dependencies across adjacent copy number states and investigate MCMC stochastic search techniques for posterior inference. Our approach results in a unified modeling framework for simultaneously inferring copy number variants (CNV) and identifying their significant associations with mRNA transcripts abundance. We show performance on simulated data and illustrate an application to data from a genomic study on human cancer cell lines.

The use of transcriptomics to unveil the role of nutrients in Mammalian liver

Liver is the organ primarily responding to diet, and it is crucial in determining plasma carbohydrate, protein, and lipid levels. In addition, it is mainly responsible for transformation of xenobiotics. For these reasons, it has been a target of transcriptomic analyses. In this review, we have covered the works dealing with the response of mammalian liver to different nutritional stimuli such as fasting/feeding, caloric restriction, dietary carbohydrate, cholesterol, fat, protein, bile acid, salt, vitamin, and oligoelement contents. Quality of fats or proteins has been equally addressed, and has the influence of minor dietary components. Other compounds, not purely nutritional as those represented by alcohol and food additives, have been included due to their relevance in processed food. The influence has been studied not only on mRNA but also on miRNA. The wide scope of the technology clearly reflects that any simple intervention has profound changes in many metabolic parameters and that there is a synergy in response when more compounds are included in the intervention. Standardized arrays to systematically test the same genes in all studies and analyzing data to establish patterns of response are required, particularly for RNA sequencing. Moreover, RNA is a valuable, easy-screening ally but always requires further confirmation.

Effect of fasting on the metabolic response of liver to experimental burn injury

Liver metabolism is altered after systemic injuries such as burns and trauma. These changes have been elucidated in rat models of experimental burn injury where the liver was isolated and perfused ex vivo. Because these studies were performed in fasted animals to deplete glycogen stores, thus simplifying quantification of gluconeogenesis, these observations reflect the combined impact of fasting and injury on liver metabolism. Herein we asked whether the metabolic response to experimental burn injury is different in fed vs. fasted animals. Rats were subjected to a cutaneous burn covering 20% of the total body surface area, or to similar procedures without administering the burn, hence a sham-burn. Half of the animals in the burn and sham-burn groups were fasted starting on postburn day 3, and the others allowed to continue ad libitum. On postburn day 4, livers were isolated and perfused for 1 hour in physiological medium supplemented with 10% hematocrit red blood cells. The uptake/release rates of major carbon and nitrogen sources, oxygen, and carbon dioxide were measured during the perfusion and the data fed into a mass balance model to estimate intracellular fluxes. The data show that in fed animals, injury increased glucose output mainly from glycogen breakdown and minimally impacted amino acid metabolism. In fasted animals, injury did not increase glucose output but increased urea production and the uptake of several amino acids, namely glutamine, arginine, glycine, and methionine. Furthermore, sham-burn animals responded to fasting by triggering gluconeogenesis from lactate; however, in burned animals the preferred gluconeogenic substrate was amino acids. Taken together, these results suggest that the fed state prevents the burn-induced increase in hepatic amino acid utilization for gluconeogenesis. The role of glycogen stores and means to increase and/or maintain internal sources of glucose to prevent increased hepatic amino acid utilization warrant further studies.

Interactions between the circadian clock and metabolism: there are good times and bad times

An endogenous circadian (∼24 h) clock regulates rhythmic processes of physiology, metabolism and behavior in most living organisms. While able to free-run under constant conditions, the circadian clock is coupled to day : night cycles to increase its amplitude and align the phase of circadian rhythms to the right time of the day. Disruptions of the circadian clock are correlated with brain dysfunctions, cardiovascular diseases and metabolic disorders. In this review, we focus on the interactions between the circadian clock and metabolism. We discuss recent findings on circadian clock regulation of feeding behavior and rhythmic expression of metabolic genes, and present evidence of metabolic input to the circadian clock. We emphasize how misalignment of circadian clocks within the body and with environmental cycles or daily schedules leads to the increasing prevalence of metabolic syndromes in modern society.

COMPADRE: an R and web resource for pathway activity analysis by component decompositions

The analysis of biological networks has become essential to study functional genomic data. Compadre is a tool to estimate pathway/gene sets activity indexes using sub-matrix decompositions for biological networks analyses. The Compadre pipeline also includes one of the direct uses of activity indexes to detect altered gene sets. For this, the gene expression sub-matrix of a gene set is decomposed into components, which are used to test differences between groups of samples. This procedure is performed with and without differentially expressed genes to decrease false calls. During this process, Compadre also performs an over-representation test. Compadre already implements four decomposition methods [principal component analysis (PCA), Isomaps, independent component analysis (ICA) and non-negative matrix factorization (NMF)], six statistical tests (t- and f-test, SAM, Kruskal-Wallis, Welch and Brown-Forsythe), several gene sets (KEGG, BioCarta, Reactome, GO and MsigDB) and can be easily expanded. Our simulation results shown in Supplementary Information suggest that Compadre detects more pathways than over-representation tools like David, Babelomics and Webgestalt and less false positives than PLAGE. The output is composed of results from decomposition and over-representation analyses providing a more complete biological picture. Examples provided in Supplementary Information show the utility, versatility and simplicity of Compadre for analyses of biological networks.

AVAILABILITY AND IMPLEMENTATION

Compadre is freely available at http://bioinformatica.mty.itesm.mx:8080/compadre. The R package is also available at https://sourceforge.net/p/compadre.

Dynamics of short-term gene expression profiling in liver following thermal injury

BACKGROUND

Severe trauma, including burns, triggers a systemic response that significantly impacts on the liver, which plays a key role in the metabolic and immune responses aimed at restoring homeostasis. While many of these changes are likely regulated at the gene expression level, there is a need to better understand the dynamics and expression patterns of burn injury-induced genes in order to identify potential regulatory targets in the liver. Herein we characterized the response within the first 24 h in a standard animal model of burn injury using a time series of microarray gene expression data.

METHODS

Rats were subjected to a full thickness dorsal scald burn injury covering 20% of their total body surface area while under general anesthesia. Animals were saline resuscitated and sacrificed at defined time points (0, 2, 4, 8, 16, and 24 h). Liver tissues were explanted and analyzed for their gene expression profiles using microarray technology. Sham controls consisted of animals handled similarly but not burned. After identifying differentially expressed probe sets between sham and burn conditions over time, the concatenated data sets corresponding to these differentially expressed probe sets in burn and sham groups were combined and analyzed using a "consensus clustering" approach.

RESULTS

The clustering method of expression data identified 621 burn-responsive probe sets in four different co-expressed clusters. Functional characterization revealed that these four clusters are mainly associated with pro-inflammatory response, anti-inflammatory response, lipid biosynthesis, and insulin-regulated metabolism. Cluster 1 pro-inflammatory response is rapidly up-regulated (within the first 2 h) following burn injury, while Cluster 2 anti-inflammatory response is activated later on (around 8 h post-burn). Cluster 3 lipid biosynthesis is down-regulated rapidly following burn, possibly indicating a shift in the utilization of energy sources to produce acute phase proteins, which serve the anti-inflammatory response. Cluster 4 insulin-regulated metabolism was down-regulated late in the observation window (around 16 h post-burn), which suggests a potential mechanism to explain the onset of hypermetabolism, a delayed but well-known response that is characteristic of severe burns and trauma with potential adverse outcome.

CONCLUSIONS

Simultaneous analysis and comparison of gene expression profiles for both burn and sham control groups provided a more accurate estimation of the activation time, expression patterns, and characteristics of a certain burn-induced response based on which the cause-effect relationships among responses were revealed.