Siddharth Sukumaran

Circadian rhythms in gene expression: Relationship to physiology, disease, drug disposition and drug action.

Circadian rhythms (24h cycles) are observed in virtually all aspects of mammalian function from expression of genes to complex physiological processes. The master clock is present in the suprachiasmatic nucleus (SCN) in the anterior part of the hypothalamus and controls peripheral clocks present in other parts of the body. Components of this core clock mechanism regulate the circadian rhythms in genome-wide mRNA expression, which in turn regulate various biological processes. Disruption of circadian rhythms can be either the cause or the effect of various disorders including metabolic syndrome, inflammatory diseases and cancer. Furthermore, circadian rhythms in gene expression regulate both the action and disposition of various drugs and affect therapeutic efficacy and toxicity based on dosing time. Understanding the regulation of circadian rhythms in gene expression plays an important role in both optimizing the dosing time for existing drugs and in the development of new therapeutics targeting the molecular clock.

Circadian variations in gene expression in rat abdominal adipose tissue and relationship to physiology

Circadian rhythms occur in all levels of organization from expression of genes to complex physiological processes. Although much is known about the mechanism of the central clock in the suprachiasmatic nucleus, the regulation of clocks present in peripheral tissues as well as the genes regulated by those clocks is still unclear. In this study, the circadian regulation of gene expression was examined in rat adipose tissue. A rich time series involving 54 animals euthanized at 18 time points within the 24-h cycle (12:12 h light-dark) was performed. mRNA expression was examined with Affymetrix gene array chips and quantitative real-time PCR, along with selected physiological measurements. Transcription factors involved in the regulation of central rhythms were examined, and 13 showed circadian oscillations. Mining of microarray data identified 190 probe sets that showed robust circadian oscillations. Circadian regulated probe sets were further parsed into seven distinct temporal clusters, with >70% of the genes showing maximum expression during the active/dark period. These genes were grouped into eight functional categories, which were examined within the context of their temporal expression. Circadian oscillations were also observed in plasma leptin, corticosterone, insulin, glucose, triglycerides, free fatty acids, and LDL cholesterol. Circadian oscillation in these physiological measurements along with the functional categorization of these genes suggests an important role for circadian rhythms in controlling various functions in white adipose tissue including adipogenesis, energy metabolism, and immune regulation.

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.

Light-dark oscillations in the lung transcriptome: implications for lung homeostasis, repair, metabolism, disease, and drug action

Diurnal-nocturnal, or circadian-like, rhythms are 24-h variations in biological processes, evolved for the efficient functioning of living organisms. Such oscillations and their regulation in many peripheral tissues are still unclear. In this study, we used Affymetrix gene chips in a rich time-series experiment involving 54 animals killed at 18 time points within the 24-h cycle to examine light-dark cycle patterns of gene expression in rat lungs. Data mining identified 646 genes (represented by 1,006 probe sets) showing robust oscillations in expression in lung that were parsed into 8 distinct temporal clusters. Surprisingly, more than two-thirds of the probe sets showing cyclic expression peaked during the animals light/inactive period. Six core clock genes and nine clock-related genes showed rhythmic oscillations in their expression in lung. Many of the genes that peaked during the inactive period included genes related to extracellular matrix, cytoskeleton, and protein processing and trafficking, which appear to be mainly involved in the repair and remodeling of the organ. Genes coding for growth factor ligands and their receptors, which play important roles in maintaining normal lung function, also showed rhythmic expression. In addition, genes involved in the metabolism and transport of endogenous compounds, xenobiotics, and therapeutic drugs, along with genes that are biomarkers or potential therapeutic targets for many lung diseases, also exhibited 24-h cyclic oscillations, suggesting an important role for such rhythms in regulating various aspects of the physiology and pathophysiology of lung.

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.

Evidence for a glucocorticoid receptor beta splice variant in the rat and its physiological regulation in liver.

Glucocorticoids are important regulators of metabolism and immune function. Synthetic glucocorticoids are extensively used for immunosuppression/anti-inflammatory therapy. Since the glucocorticoid receptor (GR) is central to most hormone effects; its in vivo regulation will influence hormone/drug action. An alternative splice variant, GRβ, is present in humans and may function as a dominant negative regulator of GR transcriptional activity. Recently, a similar splice variant was reported in mouse, although the mechanism of alternative splicing differs from that in humans. We present evidence that a splice variant of GR with an alternative C-terminus also occurs in the rat by a mechanism of intron inclusion. A highly quantitative qRT-PCR assay for the simultaneous measurement of both splice variants in a single sample was developed in order to accurately measure their regulation. We used this assay to assess the tissue specific expression of both mRNAs, and demonstrate that GRα is predominant in all tissues. In addition, the regulation of both GRα and GRβ mRNA by various physiological factors in rat liver was assessed. GRα showed a robust circadian rhythm, which was entrained with the circadian oscillation of the endogenous hormone. Time series experiments showed that both corticosteroids and LPS but not insulin dosing resulted in the transient down-regulation of GRα mRNA. LPS treatment also resulted in down-regulation of GRβ expression. A modest up-regulation in GRβ expression was observed only in animals having chronically elevated plasma insulin concentrations. However the expression of GRβ was significantly lower than that of GRα in all cases.

Differential muscle gene expression as a function of disease progression in Goto-Kakizaki diabetic rats

The Goto-Kakizaki (GK) rat, a polygenic non-obese model of type 2 diabetes, is a useful surrogate for study of diabetes-related changes independent of obesity. GK rats and appropriate controls were killed at 4, 8, 12, 16 and 20 weeks post-weaning and differential muscle gene expression along with body and muscle weights, plasma hormones and lipids, and blood cell measurements were carried out. Gene expression analysis identified 204 genes showing 2-fold or greater differences between GK and controls in at least 3 ages. Array results suggested increased oxidative capacity in GK muscles, as well as differential gene expression related to insulin resistance, which was also indicated by HOMA-IR measurements. In addition, potential new biomarkers in muscle gene expression were identified that could be either a cause or consequence of T2DM. Furthermore, we demonstrate here the presence of chronic inflammation evident both systemically and in the musculature, despite the absence of obesity.

MECHANISTIC MODELING OF THE EFFECTS OF GLUCOCORTICOIDS AND CIRCADIAN RHYTHMS ON ADIPOKINE EXPRESSION

A mechanism-based model was developed to describe the effects of methylprednisolone (MPL), circadian rhythms, and the glucose/free fatty acid (FFA)/insulin system on leptin and adiponectin expression in white adipose tissue in rats. Fifty-four normal Wistar rats received 50 mg/kg MPL intramuscularly and were sacrificed at various times. An additional set of 54 normal Wistar rats were sacrificed at 18 time points across the 24-h light/dark cycle and served as controls. Measurements included plasma MPL, glucocorticoid receptor (GR) mRNA, leptin mRNA, adiponectin mRNA, plasma leptin, adiponectin, glucose, FFA, and insulin. MPL pharmacokinetics was described by a two-compartment model with two absorption components. All measured plasma markers and mRNA expression exhibited circadian patterns except for adiponectin and were described by Fourier harmonic functions. MPL caused significant down-regulation in GR mRNA with the nadir occurring at 5 h. MPL disrupted the circadian patterns in plasma glucose and FFA by stimulating their production. Plasma glucose and FFA subsequently caused an increase in plasma insulin. Furthermore, MPL disrupted the circadian patterns in leptin mRNA expression by stimulating its production. This rise was closely followed by an increase in plasma leptin. Both leptin mRNA and plasma leptin peaked at 12 h after MPL and eventually returned back to their circadian baselines. MPL and insulin had opposing effects on adiponectin mRNA expression and plasma adiponectin, which resulted in biphasic pharmacodynamic profiles. This small systems model quantitatively describes, integrates, and provides additional insights into various factors controlling adipokine gene expression.

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.

Circadian signatures in rat liver: from gene expression to pathways

BackgroundCircadian 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.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.ResultsWe 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.ConclusionsRather 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.