The ctenophore genome and the evolutionary origins of neural systems.

The origins of neural systems remain unresolved. In contrast to other basal metazoans, ctenophores (comb jellies) have both complex nervous and mesoderm-derived muscular systems. These holoplanktonic predators also have sophisticated ciliated locomotion, behaviour and distinct development. Here we present the draft genome of Pleurobrachia bachei, Pacific sea gooseberry, together with ten other ctenophore transcriptomes, and show that they are remarkably distinct from other animal genomes in their content of neurogenic, immune and developmental genes. Our integrative analyses place Ctenophora as the earliest lineage within Metazoa. This hypothesis is supported by comparative analysis of multiple gene families, including the apparent absence of HOX genes, canonical microRNA machinery, and reduced immune complement in ctenophores. Although two distinct nervous systems are well recognized in ctenophores, many bilaterian neuron-specific genes and genes of 'classical' neurotransmitter pathways either are absent or, if present, are not expressed in neurons. Our metabolomic and physiological data are consistent with the hypothesis that ctenophore neural systems, and possibly muscle specification, evolved independently from those in other animals.

The development of large-scale de-identified biomedical databases in the age of genomics-principles and challenges.

Contemporary biomedical databases include a wide range of information types from various observational and instrumental sources. Among the most important features that unite biomedical databases across the field are high volume of information and high potential to cause damage through data corruption, loss of performance, and loss of patient privacy. Thus, issues of data governance and privacy protection are essential for the construction of data depositories for biomedical research and healthcare. In this paper, we discuss various challenges of data governance in the context of population genome projects. The various challenges along with best practices and current research efforts are discussed through the steps of data collection, storage, sharing, analysis, and knowledge dissemination.

The role of alternative Polyadenylation in regulation of rhythmic gene expression.

BACKGROUND: Alternative transcription is common in eukaryotic cells and plays important role in regulation of cellular processes. Alternative polyadenylation results from ambiguous PolyA signals in 3' untranslated region (UTR) of a gene. Such alternative transcripts share the same coding part, but differ by a stretch of UTR that may contain important functional sites. METHODS: The methodoogy of this study is based on mathematical modeling, analytical solution, and subsequent validation by datamining in multiple independent experimental data from previously published studies. RESULTS: In this study we propose a mathematical model that describes the population dynamics of alternatively polyadenylated transcripts in conjunction with rhythmic expression such as transcription oscillation driven by circadian or metabolic oscillators. Analysis of the model shows that alternative transcripts with different turnover rates acquire a phase shift if the transcript decay rate is different. Difference in decay rate is one of the consequences of alternative polyadenylation. Phase shift can reach values equal to half the period of oscillation, which makes alternative transcripts oscillate in abundance in counter-phase to each other. Since counter-phased transcripts share the coding part, the rate of translation becomes constant. We have analyzed a few data sets collected in circadian timeline for the occurrence of transcript behavior that fits the mathematical model. CONCLUSION: Alternative transcripts with different turnover rate create the effect of rectifier. This "molecular diode" moderates or completely eliminates oscillation of individual transcripts and stabilizes overall protein production rate. In our observation this phenomenon is very common in different tissues in plants, mice, and humans. The occurrence of counter-phased alternative transcripts is also tissue-specific and affects functions of multiple biological pathways. Accounting for this mechanism is important for understanding the natural and engineering the synthetic cellular circuits.

Computational workflow for analysis of gain and loss of genes in distantly related genomes.

BACKGROUND: Early evolution of animals led to profound changes in body plan organization, symmetry and the rise of tissue complexity including formation of muscular and nervous systems. This process was associated with massive restructuring of animal genomes as well as deletion, acquisition and rapid differentiation of genes from a common metazoan ancestor. Here, we present a simple but efficient workflow for elucidation of gene gain and gene loss within major branches of the animal kingdom. METHODS: We have designed a pipeline of sequence comparison, clustering and functional annotation using 12 major phyla as illustrative examples. Specifically, for the input we used sets of ab initio predicted gene models from the genomes of six bilaterians, three basal metazoans (Cnidaria, Placozoa, Porifera), two unicellular eukaryotes (Monosiga and Capsospora) and the green plant Arabidopsis as an out-group. Due to the large amounts of data the software required a high-performance Linux cluster. The final results can be imported into standard spreadsheet analysis software and queried for the numbers and specific sets of genes absent in specific genomes, uniquely present or shared among different taxons. RESULTS AND CONCLUSIONS: The developed software is open source and available free of charge on Open Source principles. It allows the user to address a number of specific questions regarding gene gain and gene loss in particular genomes, and user-defined groups of genomes can be formulated in a type of logical expression. For example, our analysis of 12 sequenced genomes indicated that these genomes possess at least 90,000 unique genes and gene families, suggesting enormous diversity of the genome repertoire in the animal kingdom. Approximately 9% of these gene families are shared universally (homologous) among all genomes, 53% are unique to specific taxa, and the rest are shared between two or more distantly related genomes.

A systems biology approach to identify molecular pathways altered by HDAC inhibition in osteosarcoma.

Osteosarcoma (OS) is the most common primary tumor in humans and dogs affecting the skeleton, and spontaneously occurring OS in dogs serves as an extremely useful model. Unacceptable toxicities using current treatment protocols prevent further dose-intensification from being a viable option to improve patient survival and thus, novel treatment strategies must be developed. Histone deacetylase inhibitors (HDACi) have recently emerged as a promising class of therapeutics demonstrating an ability to enhance the anti-tumor activity of traditional chemotherapeutics. To date, gene expression analysis of OS cell lines treated with HDACi has not been reported, and evaluation of the resultant gene expression changes may provide insight into the mechanisms that lead to success of HDACi. Canine OS cells, treated with a clinically relevant concentration of the HDACi valproic acid (VPA), were used for expression analysis on the Affymetrix canine v2.0 genechip. Differentially expressed genes were grouped into pathways based upon functional annotation; pathway analysis was performed with MetaCore and Ingenuity Pathways Analysis software. Validation of microarray results was performed by a combination of qRT-PCR and functional/biochemical assays revealing oxidative phosphorylation, cytoskeleton remodeling, cell cycle, and ubiquitin-proteasome among those pathways most affected by HDACi. The mitomycin C-bioactivating enzyme NQ01 also demonstrated upregulation following VPA treatment, leading to synergistic reductions in cell viability. These results provide a better understanding of the mechanisms by which HDACi exert their effect in OS, and have the potential to identify biomarkers that may serve as novel targets and/or predictors of response to HDACi-containing combination therapies in OS.

Proteome of human subcutaneous adipose tissue stromal vascular fraction cells versus mature adipocytes based on DIGE.

Adipose tissue contains a heterogeneous population of mature adipocytes, endothelial cells, immune cells, pericytes, and preadipocytic stromal/stem cells. To date, a majority of proteomic analyses have focused on intact adipose tissue or isolated adipose stromal/stem cells in vitro. In this study, human subcutaneous adipose tissue from multiple depots (arm and abdomen) obtained from female donors was separated into populations of stromal vascular fraction cells and mature adipocytes. Out of 960 features detected by 2-D gel electrophoresis, a total of 200 features displayed a 2-fold up- or down-regulation relative to each cell population. The protein identity of 136 features was determined. Immunoblot analyses comparing SVF relative to adipocytes confirmed that carbonic anhydrase II was up-regulated in both adipose depots while catalase was up-regulated in the arm only. Bioinformatic analyses of the data set determined that cytoskeletal, glycogenic, glycolytic, lipid metabolic, and oxidative stress related pathways were highly represented as differentially regulated between the mature adipocytes and stromal vascular fraction cells. These findings extend previous reports in the literature with respect to the adipose tissue proteome and the consequences of adipogenesis. The proteins identified may have value as biomarkers for monitoring the physiology and pathology of cell populations within subcutaneous adipose depots.

Rhythms and synchronization patterns in gene expression in the Aedes aegypti mosquito.

BACKGROUND: Aedes aegypti is arguably the most studied of all mosquito species in the laboratory and is the primary vector of both Dengue and Yellow Fever flaviviruses in the field. A large number of transcriptional studies have been made in the species and these usually report transcript quantities observed at a certain age or stage of development. However, circadian oscillation is an important characteristic of gene expression in many animals and plants, modulating both their physiology and behavior. Circadian gene expression in mosquito species has been previously reported but for only a few genes directly involved in the function of the molecular clock. RESULTS: Herein we analyze the transcription profiles of 21,494 messenger RNAs using an Ae. aegypti Agilent® microarray. Transcripts were quantified in adult female heads at 24 hours and then again at 72 hours and eight subsequent time points spaced four hours apart. We document circadian rhythms in multiple molecular pathways essential for growth, development, immune response, detoxification/pesticide resistance. Circadian rhythms were also noted in ribosomal protein genes used for normalization in reverse transcribed PCR (RT-PCR) to determine transcript abundance. We report pervasive oscillations and intricate synchronization patterns relevant to all known biological pathways. CONCLUSION: These results argue strongly that transcriptional analyses either need to be made over time periods rather than confining analyses to a single time point or development stage or exceptional care needs to be made to synchronize all mosquitoes to be analyzed and compared among treatment groups.

Small RNA profiling of Dengue virus-mosquito interactions implicates the PIWI RNA pathway in anti-viral defense.

BACKGROUND: Small RNA (sRNA) regulatory pathways (SRRPs) are important to anti-viral defence in mosquitoes. To identify critical features of the virus infection process in Dengue serotype 2 (DENV2)-infected Ae. aegypti, we deep-sequenced small non-coding RNAs. Triplicate biological replicates were used so that rigorous statistical metrics could be applied. RESULTS: In addition to virus-derived siRNAs (20-23 nts) previously reported for other arbovirus-infected mosquitoes, we show that PIWI pathway sRNAs (piRNAs) (24-30 nts) and unusually small RNAs (usRNAs) (13-19 nts) are produced in DENV-infected mosquitoes. We demonstrate that a major catalytic enzyme of the siRNA pathway, Argonaute 2 (Ago2), co-migrates with a ~1 megadalton complex in adults prior to bloodfeeding. sRNAs were cloned and sequenced from Ago2 immunoprecipitations. Viral sRNA patterns change over the course of infection. Host sRNAs were mapped to the published aedine transcriptome and subjected to analysis using edgeR (Bioconductor). We found that sRNA profiles are altered early in DENV2 infection, and mRNA targets from mitochondrial, transcription/translation, and transport functional categories are affected. Moreover, small non-coding RNAs (ncRNAs), such as tRNAs, spliceosomal U RNAs, and snoRNAs are highly enriched in DENV-infected samples at 2 and 4 dpi. CONCLUSIONS: These data implicate the PIWI pathway in anti-viral defense. Changes to host sRNA profiles indicate that specific cellular processes are affected during DENV infection, such as mitochondrial function and ncRNA levels. Together, these data provide important progress in understanding the DENV2 infection process in Ae. aegypti.

Circadian rhythms in adipose tissue: an update.

PURPOSE OF REVIEW: Over the past decade, evidence has accumulated from basic science, clinical and epidemiological studies linking circadian mechanisms to adipose tissue biology and its related comorbidities, diabetes, metabolic syndrome and obesity. This review highlights recent in-vitro and in-vivo findings from murine, human and model organism studies. RECENT FINDINGS: High-fat diets attenuate circadian mechanisms in murine adipose depots and these effects appear to be due to obesity rather than hyperglycemia. Deletion of circadian regulatory genes such as AMPK1 and nocturnin alter the circadian biology of adipose tissue. Unlike the mouse, circadian gene oscillation in human adipose tissue appears to be independent of BMI and diabetes status, suggesting that circadian mechanistic variation occurs across species. Clues for future directions in this emerging field come from studies of the hibernation and torpor state in mammals and infection models involving the Drosophila metabolic organ or 'fat body'. SUMMARY: There is a growing consensus that circadian rhythms and metabolism are tightly regulated in adipose tissue and peripheral metabolic organs. Although central mechanisms are critical, autonomous clocks exist within the adipocytes themselves. Future circadian advances are likely to result from the studies of adipose tissue-specific gene deletions.

Resveratrol suppresses human colon cancer cell proliferation and induces apoptosis via targeting the pentose phosphate and the talin-FAK signaling pathways-A proteomic approach.

BACKGROUND: We and others have previously reported that resveratrol (RSV) suppresses colon cancer cell proliferation and elevates apoptosis in vitro and/or in vivo, however molecular mechanisms are not fully elucidated. Particularly, little information is available on RSV's effects on metabolic pathways and the cell-extra cellular matrix (ECM) communication that are critical for cancer cell growth. To identify important targets of RSV, we analyzed whole protein fractions from HT-29 advanced human colon cancer cell line treated with solvent control, IGF-1 (10 nM) and RSV (150 µM) using LC/MS/MS-Mud PIT (Multidimensional Protein Identification Technology). RESULTS: Pentose phosphate pathway (PPP), a vital metabolic pathway for cell cycle progression, was elevated and suppressed by IGF-1 and RSV, respectively in the HT-29 cell line. Enzymatic assays confirmed RSV suppression of glucose-6 phosphate dehydrogenase (rate limiting) and transketolase, key enzymes of the PPP. RSV (150 µM) suppressed, whereas IGF-1 (10 nM) elevated focal adhesion complex (FAC) proteins, talin and pFAK, critical for the cell-ECM communication. Western blotting analyses confirmed the suppression or elevation of these proteins in HT-29 cancer cells treated with RSV or IGF-1, respectively. CONCLUSIONS: Proteomic analysis enabled us to establish PPP and the talin-pFAK as targets of RSV which suppress cancer cell proliferation and induce apoptosis in the colon cancer cell line HT-29. RSV (150 µM) suppressed these pathways in the presence and absence of IGF-1, suggesting its role as a chemo-preventive agent even in obese condition.

The American lobster genome reveals insights on longevity, neural, and immune adaptations.

The American lobster, Homarus americanus, is integral to marine ecosystems and supports an important commercial fishery. This iconic species also serves as a valuable model for deciphering neural networks controlling rhythmic motor patterns and olfaction. Here, we report a high-quality draft assembly of the H. americanus genome with 25,284 predicted gene models. Analysis of the neural gene complement revealed extraordinary development of the chemosensory machinery, including a profound diversification of ligand-gated ion channels and secretory molecules. The discovery of a novel class of chimeric receptors coupling pattern recognition and neurotransmitter binding suggests a deep integration between the neural and immune systems. A robust repertoire of genes involved in innate immunity, genome stability, cell survival, chemical defense, and cuticle formation represents a diversity of defense mechanisms essential to thrive in the benthic marine environment. Together, these unique evolutionary adaptations contribute to the longevity and ecological success of this long-lived benthic predator.

Yield and characterization of subcutaneous human adipose-derived stem cells by flow cytometric and adipogenic mRNA analyzes.

BACKGROUND AIMS: Adipose-derived stromal/stem cells (ASC) capable of multipotential differentiation can be isolated with high yields from human subcutaneous lipoaspirates. This study reports our recent experience of isolating and immunophenotypically characterizing ASC from >60 human patients with a mean age of 43.6 and body mass index (BMI) of 27. METHODS: We examined the ASC yield per unit volume of lipoaspirate tissue, the surface antigen profile based on flow cytometry, histochemical differentiation potential along the adipogenic and osteogenic pathways, and expression of adipogenic mRNA by transcriptomic microarray and reverse transcription (RT)-polymerase chain reaction (PCR). RESULTS: The population (n = 64) of predominantly Caucasian (84.3%) female (90.6%) donors had a mean age of 43.6 +/- 11.1 years and a mean BMI of 27.0 +/- 3.8. A yield of 375 +/- 142 x 10(3) ASC was obtained per milliliter of lipoaspirate within a 4.1 +/- 0.7-day culture period (n = 62). The ASC population was uniformly CD29(+) CD34(+) CD44(lo) CD45(lo) CD73(+) CD90(+) CD105(+) and capable of undergoing both adipogenesis and osteogenesis in vitro based on Oil Red O and Alizarin Red staining, respectively. Adipogenic differentiation was associated with a significant induction of multiple mRNA associated with lipid storage and synthesis based on microarray analysis of n = 3 donors. During an adipogenic differentiation time-course, representative mRNA (adiponectin, C/EBPalpha, leptin and LPL) displayed increases of several orders of magnitude. CONCLUSIONS: These findings demonstrate the reproducibility of subcutaneous lipoaspirates as a consistent and abundant source of functional ASC from donors across a spectrum of ages and BMI. These results have relevance for regenerative medical applications exploiting autologous and allogeneic ASC for soft and hard tissue engineering.

Expression profiling in canine osteosarcoma: identification of biomarkers and pathways associated with outcome.

BACKGROUND: Osteosarcoma (OSA) spontaneously arises in the appendicular skeleton of large breed dogs and shares many physiological and molecular biological characteristics with human OSA. The standard treatment for OSA in both species is amputation or limb-sparing surgery, followed by chemotherapy. Unfortunately, OSA is an aggressive cancer with a high metastatic rate. Characterization of OSA with regard to its metastatic potential and chemotherapeutic resistance will improve both prognostic capabilities and treatment modalities. METHODS: We analyzed archived primary OSA tissue from dogs treated with limb amputation followed by doxorubicin or platinum-based drug chemotherapy. Samples were selected from two groups: dogs with disease free intervals (DFI) of less than 100 days (n = 8) and greater than 300 days (n = 7). Gene expression was assessed with Affymetrix Canine 2.0 microarrays and analyzed with a two-tailed t-test. A subset of genes was confirmed using qRT-PCR and used in classification analysis to predict prognosis. Systems-based gene ontology analysis was conducted on genes selected using a standard J5 metric. The genes identified using this approach were converted to their human homologues and assigned to functional pathways using the GeneGo MetaCore platform. RESULTS: Potential biomarkers were identified using gene expression microarray analysis and 11 differentially expressed (p < 0.05) genes were validated with qRT-PCR (n = 10/group). Statistical classification models using the qRT-PCR profiles predicted patient outcomes with 100% accuracy in the training set and up to 90% accuracy upon stratified cross validation. Pathway analysis revealed alterations in pathways associated with oxidative phosphorylation, hedgehog and parathyroid hormone signaling, cAMP/Protein Kinase A (PKA) signaling, immune responses, cytoskeletal remodeling and focal adhesion. CONCLUSIONS: This profiling study has identified potential new biomarkers to predict patient outcome in OSA and new pathways that may be targeted for therapeutic intervention.

Computational analysis of gene expression space associated with metastatic cancer.

BACKGROUND: Prostate carcinoma is among the most common types of cancer affecting hundreds of thousands people every year. Once the metastatic form of prostate carcinoma is documented, the majority of patients die from their tumors as opposed to other causes. The key to successful treatment is in the earliest possible diagnosis, as well as understanding the molecular mechanisms of metastatic progression. A number of recent studies have identified multiple biomarkers for metastatic progression. However, most of the studies consider only direct comparison between metastatic and non-metastatic classes of samples. RESULTS: We propose an alternative concept of analysis that considers the entire multidimensional space of gene expression and identifies the partition of this space in which metastatic development is possible. To apply this concept in cancer gene expression studies we utilize a modification of high-dimension natural taxonomy algorithm FOREL. Our analysis of microarray data containing primary and metastatic cancer samples has revealed not only differentially expressed genes, but also relations between different groups of primary and metastatic cancer. Metastatic samples tend to occupy a distinct partition of gene expression space. Further pathway analysis suggests that this partition is delineated by a specific pattern of gene expression in cytoskeleton remodeling, cell adhesion and apoptosis/cell survival pathways. We compare our findings with both report of original analysis and recent studies in molecular mechanism of metastasis. CONCLUSION: Our analysis indicates a single molecular mechanism of metastasis. The new approach does not contradict previously reported findings, but reveals important details unattainable with traditional methodology.

Persistent fetal infection with bovine viral diarrhea virus differentially affects maternal blood cell signal transduction pathways.

The consequences of viral infection during pregnancy include impact on fetal and maternal immune responses and on fetal development. Transplacental infection in cattle with noncytopathic bovine viral diarrhea virus (ncpBVDV) during early gestation results in persistently infected (PI) fetuses with life-long viremia and susceptibility to infections. Infection of the fetus during the third trimester or after birth leads to a transient infection cleared by a competent immune system. We hypothesized that ncpBVDV infection and presence of an infected fetus would alter immune response and lead to downregulation of proinflammatory processes in pregnant dams. Naïve pregnant heifers were challenged with ncpBVDV2 on day 75 (PI fetus) and day 175 [transiently infected (TI) fetus] or kept uninfected (healthy control fetus). Maternal blood samples were collected up to day 190 of gestation. Genome-wide microarray analysis of gene expression in maternal peripheral white blood cells, performed on days 160 and 190 of gestation, revealed multiple signal transduction pathways affected by ncpBVDV infection. Acute infection and presence of a TI fetus caused upregulation of the type I interferon (IFN) pathway genes, including dsRNA sensors and IFN-stimulated genes. The presence of a PI fetus caused prolonged downregulation of chemokine receptor 4 (CXCR4) and T cell receptor (TCR) signaling in maternal blood cells. We conclude that: 1) infection with ncpBVDV induces a vigorous type I IFN response, and 2) presence of a PI fetus causes downregulation of important signaling pathways in the blood of the dam, which could have deleterious consequences on fetal development and the immune response.

Comprehensive analysis of circadian periodic pattern in plant transcriptome.

BACKGROUND: Circadian rhythm is a crucial factor in orchestration of plant physiology, keeping it in synchrony with the daylight cycle. Previous studies have reported that up to 16% of plant transcriptome are circadially expressed. RESULTS: Our studies of mammalian gene expression revealed circadian baseline oscillation in nearly 100% of genes. Here we present a comprehensive analysis of periodicity in two independent data sets. Application of the advanced algorithms and analytic approached already tested on animal data reveals oscillation in almost every gene of Arabidopsis thaliana. CONCLUSION: This study indicates an even more pervasive role of oscillation in molecular physiology of plants than previously believed. Earlier studies have dramatically underestimated the prevalence of circadian oscillation in plant gene expression.

Systems biology approach to identification of biomarkers for metastatic progression in cancer.

BACKGROUND: Metastases are responsible for the majority of cancer fatalities. The molecular mechanisms governing metastasis are poorly understood, hindering early diagnosis and treatment. Previous studies of gene expression patterns in metastasis have concentrated on selection of a small number of "signature" biomarkers. RESULTS: We propose an alternative approach that puts into focus gene interaction networks and molecular pathways rather than separate genes. We have reanalyzed expression data from a large set of primary solid and metastatic tumors originating from different tissues using the latest available tools for normalization, identification of differentially expressed genes and pathway analysis. Our studies indicate that regardless of the tissue of origin, all metastatic tumors share a number of common features related to changes in basic energy metabolism, cell adhesion/cytoskeleton remodeling, antigen presentation and cell cycle regulation. Analysis of multiple independent datasets indicates significantly reduced oxidative phosphorylation in metastases compared to primary solid tumors. CONCLUSION: Our methods allow identification of robust, although not necessarily highly expressed biomarkers. A systems approach relying on groups of interacting genes rather than single markers is also essential for understanding the cellular processes leading to metastatic progression. We have identified metabolic pathways associated with metastasis that may serve as novel targets for therapeutic intervention.

Digital signal processing reveals circadian baseline oscillation in majority of mammalian genes.

In mammals, circadian periodicity has been described for gene expression in the hypothalamus and multiple peripheral tissues. It is accepted that 10%-15% of all genes oscillate in a daily rhythm, regulated by an intrinsic molecular clock. Statistical analyses of periodicity are limited by the small size of datasets and high levels of stochastic noise. Here, we propose a new approach applying digital signal processing algorithms separately to each group of genes oscillating in the same phase. Combined with the statistical tests for periodicity, this method identifies circadian baseline oscillation in almost 100% of all expressed genes. Consequently, circadian oscillation in gene expression should be evaluated in any study related to biological pathways. Changes in gene expression caused by mutations or regulation of environmental factors (such as photic stimuli or feeding) should be considered in the context of changes in the amplitude and phase of genetic oscillations.

Circadian succession of molecular processes in living tissues.

BACKGROUND: Oscillations of different origin, period and amplitude play an important role in the regulation of cellular processes. Most widely studied is the circadian or approximately daily variation in gene expression activity. Timing of gene expression is controlled by internal molecular clock keeping steady periodic expression. In this study, we shift attention towards a broad range of periodically expressed genes involved in multiple cellular functions which may or may not be under direct control of the intrinsic circadian clock. Are all molecular functions represented in expressed genes at all times? Alternatively, are different molecular functions performed at different times? Is there a pattern of succession for molecular processes and functions throughout their daily activity period? RESULTS: To answer these questions, we re-analyzed a number of mouse circadian gene expression data available from public sources. These data represent the normal function of metabolically active peripheral tissues (white adipose tissue, brown adipose tissue, liver). We applied novel methods for the estimation of confidence in phase assignment to identify groups of synchronous genes peaking at the same time regardless of the amplitude or the absolute intensity of expression. Each synchronous group has been annotated to identify Gene Ontology (GO) terms and molecular pathways. Our analysis identified molecular functions specific to a particular time of the day in different tissues. CONCLUSION: Improved methodology for datamining allowed for the discovery of functions and biological pathways in groups of genes with synchronized peak expression time. In particular, such functions as oxidative phase of energy metabolism, DNA repair, mRNA processing, lipid biosynthesis and others are separated in time. This timewise compartmentalization is important for understanding the cellular circuitry and can be used to optimize the time of intervention with drug or genome medication.

Circadian oscillation of gene expression in murine calvarial bone.

UNLABELLED: The genes encoding the core circadian transcription factors display an oscillating expression profile in murine calvarial bone. More than 26% of the calvarial bone transcriptome exhibits a circadian rhythm, comparable with that observed in brown and white adipose tissues and liver. Thus, circadian mechanisms may directly modulate oxidative phosphorylation and multiple metabolic pathways in bone homeostasis. INTRODUCTION: Although circadian rhythms have been associated historically with central regulatory mechanisms, there is emerging evidence that the circadian transcriptional apparatus exists in peripheral tissues. The aim of this study was to determine the presence and extent of circadian oscillation in the transcriptome of murine calvarial bone. MATERIALS AND METHODS: Cohorts of 8-week-old male AKR/J mice were maintained in a controlled 12-h light:12-h dark cycle on an ad libitum diet for 2 weeks. Groups of three mice were killed every 4 h over a 48-h period. The level of gene expression at successive times-points was determined by quantitative RT-PCR and Affymetrix microarray. Data were analyzed using multiple statistical time series algorithms, including Cosinor, Fisher g-test, and the permutation time test. RESULTS: Both the positive (Bmal1, Npas2) and negative (Cry1, Cry2, Per1, Per2, Per3) elements of the circadian transcriptional apparatus and their immediate downstream targets and mediators (Dbp, Rev-erbalpha, Rev-erbbeta) exhibited oscillatory expression profiles. Consistent with findings in other tissues, the positive and negative elements were in antiphase relative to each other. More than 26% of the genes present on the microarray displayed an oscillatory profile in calvarial bone, comparable with the levels observed in brown and white adipose tissues and liver; however, only a subset of 174 oscillating genes were shared among all four tissues. CONCLUSIONS: Our findings show that the components of the circadian transcriptional apparatus are represented in calvarial bone and display coordinated oscillatory behavior. However, these are not the only genes to display an oscillatory expression profile, which is seen in multiple pathways involving oxidative phosphorylation and lipid, protein, and carbohydrate metabolism.