Emergent genome-wide control in wildtype and genetically mutated lipopolysaccarides-stimulated macrophages

Damage-mediated macrophage polarization in sterile inflammation

Most of the leading causes of death, such as cardiovascular diseases, cancer, dementia, neurodegenerative diseases, and many more, are associated with sterile inflammation, either as a cause or a consequence of these conditions. The ability to control the progression of inflammation toward tissue resolution before it becomes chronic holds significant clinical potential. During sterile inflammation, the initiation of inflammation occurs through damage-associated molecular patterns (DAMPs) in the absence of pathogen-associated molecules. Macrophages, which are primarily localized in the tissue, play a pivotal role in sensing DAMPs. Furthermore, macrophages can also detect and respond to resolution-associated molecular patterns (RAMPs) and specific pro-resolving mediators (SPMs) during sterile inflammation. Macrophages, being highly adaptable cells, are particularly influenced by changes in the microenvironment. In response to the tissue environment, monocytes, pro-inflammatory macrophages, and pro-resolution macrophages can modulate their differentiation state. Ultimately, DAMP and RAMP-primed macrophages, depending on the predominant subpopulation, regulate the balance between inflammatory and resolving processes. While sterile injury and pathogen-induced reactions may have distinct effects on macrophages, most studies have focused on macrophage responses induced by pathogens. In this review, which emphasizes available human data, we illustrate how macrophages sense these mediators by examining the expression of receptors for DAMPs, RAMPs, and SPMs. We also delve into the signaling pathways induced by DAMPs, RAMPs, and SPMs, which primarily contribute to the regulation of macrophage differentiation from a pro-inflammatory to a pro-resolution phenotype. Understanding the regulatory mechanisms behind the transition between macrophage subtypes can offer insights into manipulating the transition from inflammation to resolution in sterile inflammatory diseases.

Is Cancer Reversible? Rethinking Carcinogenesis Models-A New Epistemological Tool

A growing number of studies shows that it is possible to induce a phenotypic transformation of cancer cells from malignant to benign. This process is currently known as "tumor reversion". However, the concept of reversibility hardly fits the current cancer models, according to which gene mutations are considered the primary cause of cancer. Indeed, if gene mutations are causative carcinogenic factors, and if gene mutations are irreversible, how long should cancer be considered as an irreversible process? In fact, there is some evidence that intrinsic plasticity of cancerous cells may be therapeutically exploited to promote a phenotypic reprogramming, both in vitro and in vivo. Not only are studies on tumor reversion highlighting a new, exciting research approach, but they are also pushing science to look for new epistemological tools capable of better modeling cancer.

The sesquiterpene lactone estafiatin exerts anti-inflammatory effects on macrophages and protects mice from sepsis induced by LPS and cecal ligation puncture

BACKGROUND

Previously, we found that the water extract of Artermisia scoparia Waldst. & Kit suppressed the cytokine production of lipopolysaccharide (LPS)-stimulated macrophages and alleviated carrageenan-induced acute inflammation in mice. Artemisia contains various sesquiterpene lactones and most of them exert immunomodulatory activity.

PURPOSE

In the present study, we investigated the immunomodulatory effect of estafiatin (EST), a sesquiterpene lactone derived from A. scoparia, on LPS-induced inflammation in macrophages and mouse sepsis model.

STUDY DESIGN AND METHODS

Murine bone marrow-derived macrophages (BMDMs) and THP-1 cells, a human monocytic leukemia cell line, were pretreated with different doses of EST for 2 h, followed by LPS treatment. The gene and protein expression of pro-inflammatory cytokines interleukin (IL)-6, tumor necrosis factor (TNF)-α, and inducible nitric oxide synthase (iNOS) were measured by quantitative real-time polymerase chain reaction (qPCR) and Western blot analysis. The activation of nuclear factor kappa B (NF-κB) and mitogen-activated protein kinases (MAPKs) was also evaluated at the level of phosphorylation. The effect of EST on inflammatory cytokine production, lung histopathology, and survival rate was assessed in an LPS-induced mice model of septic shock. The effect of EST on the production of cytokines in LPS-stimulated peritoneal macrophages was evaluated by in vitro and ex vivo experiments and protective effect of EST on cecal ligation and puncture (CLP) mice was also assessed.

RESULTS

The LPS-induced expression of IL-6, TNF-α, and iNOS was suppressed at the mRNA and protein levels in BMDMs and THP-1 cells, respectively, by pretreatment with EST. The half-maximal inhibitory concentration (IC50) of EST on IL-6 and TNF-α production were determined as 3.2 μM and 3.1 μM in BMDMs, 3 μM and 3.4 μM in THP1 cells, respectively. In addition, pretreatment with EST significantly reduced the LPS-induced phosphorylation p65, p38, JNK, and ERK in both cell types. In the LPS-induced mice model of septic shock, serum levels of IL-6, TNF-α, IL-1β, CXCL1, and CXCL2 were lower in EST-treated mice than in the control animals. Histopathology analysis revealed that EST treatment ameliorated LPS-induced lung damage. Moreover, while 1 of 7 control mice given lethal dose of LPS survived, 3 of 7 EST-treated (1.25 mg/kg) mice and 5 of 7 EST-treated (2.5 mg/kg) mice were survived. Pretreatment of EST dose-dependently suppressed the LPS-induced production of IL-6, TNF-α and CXCL1 in peritoneal macrophages. In CLP-induced mice sepsis model, while all 6 control mice was dead at 48 h, 1 of 6 EST-treated (1.25 mg/kg) mice and 3 of 6 EST-treated (2.5 mg/kg) mice survived for 96 h.

CONCLUSION

These results demonstrated that EST exerts anti-inflammatory effects on LPS-stimulated macrophages and protects mice from sepsis. Our study suggests that EST could be developed as a new therapeutic agent for sepsis and various inflammatory diseases.

ScatLay: utilizing transcriptome-wide noise for identifying and visualizing differentially expressed genes

Differential expressed (DE) genes analysis is valuable for understanding comparative transcriptomics between cells, conditions or time evolution. However, the predominant way of identifying DE genes is to use arbitrary threshold fold or expression changes as cutoff. Here, we developed a more objective method, Scatter Overlay or ScatLay, to extract and graphically visualize DE genes across any two samples by utilizing their pair-wise scatter or transcriptome-wide noise, while factoring replicate variabilities. We tested ScatLay for 3 cell types: between time points for Escherichia coli aerobiosis and Saccharomyces cerevisiae hypoxia, and between untreated and Etomoxir treated Mus Musculus embryonic stem cell. As a result, we obtain 1194, 2061 and 2932 DE genes, respectively. Next, we compared these data with two widely used current approaches (DESeq2 and NOISeq) with typical twofold expression changes threshold, and show that ScatLay reveals significantly larger number of DE genes. Hence, our method provides a wider coverage of DE genes, and will likely pave way for finding more novel regulatory genes in future works.

Attractor Concepts to Evaluate the Transcriptome-wide Dynamics Guiding Anaerobic to Aerobic State Transition in Escherichia coli

For any dynamical system, like living organisms, an attractor state is a set of variables or mechanisms that converge towards a stable system behavior despite a wide variety of initial conditions. Here, using multi-dimensional statistics, we investigate the global gene expression attractor mechanisms shaping anaerobic to aerobic state transition (AAT) of Escherichia coli in a bioreactor at early times. Out of 3,389 RNA-Seq expression changes over time, we identified 100 sharply changing genes that are key for guiding 1700 genes into the AAT attractor basin. Collectively, these genes were named as attractor genes constituting of 6 dynamic clusters. Apart from the expected anaerobic (glycolysis), aerobic (TCA cycle) and fermentation (succinate pathways) processes, sulphur metabolism, ribosome assembly and amino acid transport mechanisms together with 332 uncharacterised genes are also key for AAT. Overall, our work highlights the importance of multi-dimensional statistical analyses for revealing novel processes shaping AAT.

ABioTrans: A Biostatistical Tool for Transcriptomics Analysis

Here we report a bio-statistical/informatics tool, ABioTrans, developed in R for gene expression analysis. The tool allows the user to directly read RNA-Seq data files deposited in the Gene Expression Omnibus or GEO database. Operated using any web browser application, ABioTrans provides easy options for multiple statistical distribution fitting, Pearson and Spearman rank correlations, PCA, k-means and hierarchical clustering, differential expression (DE) analysis, Shannon entropy and noise (square of coefficient of variation) analyses, as well as Gene ontology classifications.

The GABAA Receptor α2 Subunit Activates a Neuronal TLR4 Signal in the Ventral Tegmental Area that Regulates Alcohol and Nicotine Abuse

Alcoholism initiates with episodes of excessive alcohol drinking, known as binge drinking, which is one form of excessive drinking (NIAAA Newsletter, 2004) that is related to impulsivity and anxiety (Ducci et al., 2007; Edenberg et al., 2004) and is also predictive of smoking status. The predisposition of non-alcohol exposed subjects to initiate binge drinking is controlled by neuroimmune signaling that includes an innately activated neuronal Toll-like receptor 4 (TLR4) signal. This signal also regulates cognitive impulsivity, a heritable trait that defines drug abuse initiation. However, the mechanism of signal activation, its function in dopaminergic (TH+) neurons within the reward circuitry implicated in drug-seeking behavior [viz. the ventral tegmental area (VTA)], and its contribution to nicotine co-abuse are still poorly understood. We report that the γ-aminobutyric acid_A receptor (GABA_AR) α2 subunit activates the TLR4 signal in neurons, culminating in the activation (phosphorylation/nuclear translocation) of cyclic AMP response element binding (CREB) but not NF-kB transcription factors and the upregulation of corticotropin-releasing factor (CRF) and tyrosine hydroxylase (TH). The signal is activated through α2/TLR4 interaction, as evidenced by co-immunoprecipitation, and it is present in the VTA from drug-untreated alcohol-preferring P rats. VTA infusion of neurotropic herpes simplex virus (HSV) vectors for α2 (pHSVsiLA2) or TLR4 (pHSVsiTLR4) but not scrambled (pHSVsiNC) siRNA inhibits signal activation and both binge alcohol drinking and nicotine sensitization, suggesting that the α2-activated TLR4 signal contributes to the regulation of both alcohol and nicotine abuse.

Self-Organizing Global Gene Expression Regulated through Criticality: Mechanism of the Cell-Fate Change

BACKGROUND

A fundamental issue in bioscience is to understand the mechanism that underlies the dynamic control of genome-wide expression through the complex temporal-spatial self-organization of the genome to regulate the change in cell fate. We address this issue by elucidating a physically motivated mechanism of self-organization.

PRINCIPAL FINDINGS

Building upon transcriptome experimental data for seven distinct cell fates, including early embryonic development, we demonstrate that self-organized criticality (SOC) plays an essential role in the dynamic control of global gene expression regulation at both the population and single-cell levels. The novel findings are as follows: i) Mechanism of cell-fate changes: A sandpile-type critical transition self-organizes overall expression into a few transcription response domains (critical states). A cell-fate change occurs by means of a dissipative pulse-like global perturbation in self-organization through the erasure of initial-state critical behaviors (criticality). Most notably, the reprogramming of early embryo cells destroys the zygote SOC control to initiate self-organization in the new embryonal genome, which passes through a stochastic overall expression pattern. ii) Mechanism of perturbation of SOC controls: Global perturbations in self-organization involve the temporal regulation of critical states. Quantitative evaluation of this perturbation in terminal cell fates reveals that dynamic interactions between critical states determine the critical-state coherent regulation. The occurrence of a temporal change in criticality perturbs this between-states interaction, which directly affects the entire genomic system. Surprisingly, a sub-critical state, corresponding to an ensemble of genes that shows only marginal changes in expression and consequently are considered to be devoid of any interest, plays an essential role in generating a global perturbation in self-organization directed toward the cell-fate change.

CONCLUSION AND SIGNIFICANCE

'Whole-genome' regulation of gene expression through self-regulatory SOC control complements gene-by-gene fine tuning and represents a still largely unexplored non-equilibrium statistical mechanism that is responsible for the massive reprogramming of genome expression.

SMT and TOFT: Why and How They are Opposite and Incompatible Paradigms

The Somatic Mutation Theory (SMT) has been challenged on its fundamentals by the Tissue Organization Field Theory of Carcinogenesis (TOFT). However, a recent publication has questioned whether TOFT could be a valid alternative theory of carcinogenesis to that presented by SMT. Herein we critically review arguments supporting the irreducible opposition between the two theoretical approaches by highlighting differences regarding the philosophical, methodological and experimental approaches on which they respectively rely. We conclude that SMT has not explained carcinogenesis due to severe epistemological and empirical shortcomings, while TOFT is gaining momentum. The main issue is actually to submit SMT to rigorous testing. This concern includes the imperatives to seek evidence for disproving one's hypothesis, and to consider the whole, and not just selective evidence.

Emergent Self-Organized Criticality in Gene Expression Dynamics: Temporal Development of Global Phase Transition Revealed in a Cancer Cell Line

BACKGROUND

The underlying mechanism of dynamic control of the genome-wide expression is a fundamental issue in bioscience. We addressed it in terms of phase transition by a systemic approach based on both density analysis and characteristics of temporal fluctuation for the time-course mRNA expression in differentiating MCF-7 breast cancer cells.

METHODOLOGY

In a recent work, we suggested criticality as an essential aspect of dynamic control of genome-wide gene expression. Criticality was evident by a unimodal-bimodal transition through flattened unimodal expression profile. The flatness on the transition suggests the existence of a critical transition at which up- and down-regulated expression is balanced. Mean field (averaging) behavior of mRNAs based on the temporal expression changes reveals a sandpile type of transition in the flattened profile. Furthermore, around the transition, a self-similar unimodal-bimodal transition of the whole expression occurs in the density profile of an ensemble of mRNA expression. These singular and scaling behaviors identify the transition as the expression phase transition driven by self-organized criticality (SOC).

PRINCIPAL FINDINGS

Emergent properties of SOC through a mean field approach are revealed: i) SOC, as a form of genomic phase transition, consolidates distinct critical states of expression, ii) Coupling of coherent stochastic oscillations between critical states on different time-scales gives rise to SOC, and iii) Specific gene clusters (barcode genes) ranging in size from kbp to Mbp reveal similar SOC to genome-wide mRNA expression and ON-OFF synchronization to critical states. This suggests that the cooperative gene regulation of topological genome sub-units is mediated by the coherent phase transitions of megadomain-scaled conformations between compact and swollen chromatin states.

CONCLUSION AND SIGNIFICANCE

In summary, our study provides not only a systemic method to demonstrate SOC in whole-genome expression, but also introduces novel, physically grounded concepts for a breakthrough in the study of biological regulation.

Tracking global gene expression responses in T cell differentiation

Upon receiving antigens from the innate immune cells, CD4(+) T cells differentiate into distinct effector cells. To probe the global responses of distinct effector cells, we analyzed transcriptome-wide expressions of Th1, Th2, Treg and Th17 using Pearson correlation, entropy and principal component analyses, with Th0 as a control. Although the global response of Th0 was quite distinct from Th17, surprisingly, it was highly similar to Th1, Th2 and Treg. Moreover, 8 major temporal groups consisting of 5704 differentially expressed genes were revealed for both Th0 and Th17. Gene functional enrichment analysis showed immune responses and metabolic processes were mainly activated between Th0 and Th17, while genes related to cell cycle and replication were differentially regulated. Moreover, we found the upregulation of several novel genes for Th0 and Th17. Overall, we deduce that Th0 is globally similar to Th1, Th2 and Treg. Our results indicate that Th0 is a differentiated state and, therefore, may not be used as a control cell type.

Transcriptome-wide variability in single embryonic development cells

Molecular heterogeneity of individual molecules within single cells has been recently shown to be crucial for cell fate diversifications. However, on a global scale, the effect of molecular variability for embryonic developmental stages is largely underexplored. Here, to understand the origins of transcriptome-wide variability of oocytes to blastocysts in human and mouse, we examined RNA-Seq datasets. Evaluating Pearson correlation, Shannon entropy and noise patterns (η²vs.μ), our investigations reveal a phase transition from low to saturating levels of diversity and variability of transcriptome-wide expressions through the development stages. To probe the observed behaviour further, we utilised a stochastic transcriptional model to simulate the global gene expressions pattern for each development stage. From the model, we concur that transcriptome-wide regulation initially begins from 2-cell stage, and becomes strikingly variable from 8-cell stage due to amplification and quantal transcriptional activity.

Global genetic response in a cancer cell: self-organized coherent expression dynamics

Understanding the basic mechanism of the spatio-temporal self-control of genome-wide gene expression engaged with the complex epigenetic molecular assembly is one of major challenges in current biological science. In this study, the genome-wide dynamical profile of gene expression was analyzed for MCF-7 breast cancer cells induced by two distinct ErbB receptor ligands: epidermal growth factor (EGF) and heregulin (HRG), which drive cell proliferation and differentiation, respectively. We focused our attention to elucidate how global genetic responses emerge and to decipher what is an underlying principle for dynamic self-control of genome-wide gene expression. The whole mRNA expression was classified into about a hundred groups according to the root mean square fluctuation (rmsf). These expression groups showed characteristic time-dependent correlations, indicating the existence of collective behaviors on the ensemble of genes with respect to mRNA expression and also to temporal changes in expression. All-or-none responses were observed for HRG and EGF (biphasic statistics) at around 10–20 min. The emergence of time-dependent collective behaviors of expression occurred through bifurcation of a coherent expression state (CES). In the ensemble of mRNA expression, the self-organized CESs reveals distinct characteristic expression domains for biphasic statistics, which exhibits notably the presence of criticality in the expression profile as a route for genomic transition. In time-dependent changes in the expression domains, the dynamics of CES reveals that the temporal development of the characteristic domains is characterized as autonomous bistable switch, which exhibits dynamic criticality (the temporal development of criticality) in the genome-wide coherent expression dynamics. It is expected that elucidation of the biophysical origin for such critical behavior sheds light on the underlying mechanism of the control of whole genome.

Network theory inspired analysis of time-resolved expression data reveals key players guiding P. patens stem cell development

Transcription factors (TFs) often trigger developmental decisions, yet, their transcripts are often only moderately regulated and thus not easily detected by conventional statistics on expression data. Here we present a method that allows to determine such genes based on trajectory analysis of time-resolved transcriptome data. As a proof of principle, we have analysed apical stem cells of filamentous moss (P. patens) protonemata that develop from leaflets upon their detachment from the plant. By our novel correlation analysis of the post detachment transcriptome kinetics we predict five out of 1,058 TFs to be involved in the signaling leading to the establishment of pluripotency. Among the predicted regulators is the basic helix loop helix TF PpRSL1, which we show to be involved in the establishment of apical stem cells in P. patens. Our methodology is expected to aid analysis of key players of developmental decisions in complex plant and animal systems.

Enhancing apoptosis in TRAIL-resistant cancer cells using fundamental response rules

The tumor necrosis factor related apoptosis-inducing ligand (TRAIL) induces apoptosis in malignant cells, while leaving other cells mostly unharmed. However, several carcinomas remain resistant to TRAIL. To investigate the resistance mechanisms in TRAIL-stimulated human fibrosarcoma (HT1080) cells, we developed a computational model to analyze the temporal activation profiles of cell survival (IκB, JNK, p38) and apoptotic (caspase-8 and -3) molecules in wildtype and several (FADD, RIP1, TRAF2 and caspase-8) knock-down conditions. Based on perturbation-response approach utilizing the law of information (signaling flux) conservation, we derived response rules for population-level average cell response. From this approach, i) a FADD-independent pathway to activate p38 and JNK, ii) a crosstalk between RIP1 and p38, and iii) a crosstalk between p62 and JNK are predicted. Notably, subsequent simulations suggest that targeting a novel molecule at p62/sequestosome-1 junction will optimize apoptosis through signaling flux redistribution. This study offers a valuable prospective to sensitive TRAIL-based therapy.

Anti-inflammatory modulation of immune response by probiotic Escherichia coli Nissle 1917 in human blood mononuclear cells

Escherichia coli Nissle 1917 (EcN) bears a defect in its LPS biosynthesis leading to truncated variable oligosaccharide-antigen chains and a semi-rough phenotype. It is effectively inactivated by complement factors due to resolved serum resistance and is, therefore, safe as a probiotic strain, i.e. for the treatment of inflammatory gastrointestinal diseases. It is unknown whether the modification of LPS in EcN contributes to its probiotic properties. Purified LPS from EcN and wild-type LPS from uropathogenic E. coli W536 together with raw lysates of both strains were analyzed for their gene expression activity with human PBMCs measured by microarrays. Comparing the two LPS molecules and the two lysate variants with each other, respectively, no differences of transcriptional patterns were observed. However, when comparing LPS with lysate patterns, pro-inflammatory cytokine IL-12p40 was up-regulated by both LPS molecules and anti-inflammatory IL-10 by both lysates. The higher the lysate concentration, the higher IL-10 release from PBMCs, clearly exceeding LPS induced IL-12p40 release. Furthermore, inflammatory chemokine CCL24 (eotaxin) was down-regulated by lysates and quantitative real-time PCR revealed that EcN compared to wild-type LPS was 8 times stronger in down-regulation of CCL24. We conclude that truncated LPS may down-regulate CCL24-mediated inflammation and that EcN lysate contains as yet unidentified factors which preferably induce anti-inflammatory activity. Both effects may contribute to the probiotic properties of EcN.

Collective dynamics of specific gene ensembles crucial for neutrophil differentiation: the existence of genome vehicles revealed

Cell fate decision remarkably generates specific cell differentiation path among the multiple possibilities that can arise through the complex interplay of high-dimensional genome activities. The coordinated action of thousands of genes to switch cell fate decision has indicated the existence of stable attractors guiding the process. However, origins of the intracellular mechanisms that create “cellular attractor” still remain unknown. Here, we examined the collective behavior of genome-wide expressions for neutrophil differentiation through two different stimuli, dimethyl sulfoxide (DMSO) and all-trans-retinoic acid (atRA). To overcome the difficulties of dealing with single gene expression noises, we grouped genes into ensembles and analyzed their expression dynamics in correlation space defined by Pearson correlation and mutual information. The standard deviation of correlation distributions of gene ensembles reduces when the ensemble size is increased following the inverse square root law, for both ensembles chosen randomly from whole genome and ranked according to expression variances across time. Choosing the ensemble size of 200 genes, we show the two probability distributions of correlations of randomly selected genes for atRA and DMSO responses overlapped after 48 hours, defining the neutrophil attractor. Next, tracking the ranked ensembles' trajectories, we noticed that only certain, not all, fall into the attractor in a fractal-like manner. The removal of these genome elements from the whole genomes, for both atRA and DMSO responses, destroys the attractor providing evidence for the existence of specific genome elements (named “genome vehicle”) responsible for the neutrophil attractor. Notably, within the genome vehicles, genes with low or moderate expression changes, which are often considered noisy and insignificant, are essential components for the creation of the neutrophil attractor. Further investigations along with our findings might provide a comprehensive mechanistic view of cell fate decision.

Hematopoietic differentiation: a coordinated dynamical process towards attractor stable states

Background

The differentiation process, proceeding from stem cells towards the different committed cell types, can be considered as a trajectory towards an attractor of a dynamical process. This view, taking into consideration the transcriptome and miRNome dynamics considered as a whole, instead of looking at few 'master genes' driving the system, offers a novel perspective on this phenomenon. We investigated the 'differentiation trajectories' of the hematopoietic system considering a genome-wide scenario.

Results

We developed serum-free liquid suspension unilineage cultures of cord blood (CB) CD34⁺ hematopoietic progenitor cells through erythroid (E), megakaryocytic (MK), granulocytic (G) and monocytic (Mo) pathways. These cultures recapitulate physiological hematopoiesis, allowing the analysis of almost pure unilineage precursors starting from initial differentiation of HPCs until terminal maturation. By analyzing the expression profile of protein coding genes and microRNAs in unilineage CB E, MK, G and Mo cultures, at sequential stages of differentiation and maturation, we observed a coordinated, fully interconnected and scalable character of cell population behaviour in both transcriptome and miRNome spaces reminiscent of an attractor-like dynamics. MiRNome and transcriptome space differed for a still not terminally committed behaviour of microRNAs.

Conclusions

Consistent with their roles, the transcriptome system can be considered as the state space of a cell population, while the continuously evolving miRNA space corresponds to the tuning system necessary to reach the attractor. The behaviour of miRNA machinery could be of great relevance not only for the promise of reversing the differentiated state but even for tumor biology.

Meta-analysis of lineage-specific gene expression signatures in mouse leukocyte populations

In order to address fundamental questions associated with the relationships between mononuclear phagocytes and other myeloid and lymphoid cell populations, we have taken advantage of the growing body of expression data available in the public domain. We collated a large number of published expression studies on mouse haemopoietic cell lineages comprising 304 cell samples from 29 independent experiments performed on a single microarray platform (Affymetrix MOE430-2). The data were subjected to network-based cluster analysis using Biolayout Express(3D). Genes with related function clustered together in distinct regions of the graph reaffirming many known associations between gene expression and role in specific pathways and defining most major cell types of the immune system. Promoters of genes within individual clusters were distinguished by over-representation of regulatory motifs recognised by specific transcription factors. However, these data indicate that commonly used myeloid subpopulation markers, such as CD11c (Itgax), do not correlate with expression of other genes, and further bring into question their use in defining myeloid cell lineage, activation (M1 vs. M2) and antigen-presenting cell function. In particular, there were few mRNA markers that clearly distinguished classical dendritic cells (DC) from macrophages, other than low expression of genes required for phagocytic activity. Bone marrow-derived DC, grown in GM-CSF, were clearly identified as phagocytes and distinguished from isolated lymphoid tissue DC. Thus, through pooling datasets from public data and examining the gene expression clusters within, we can learn a great deal about the transcriptional networks that underpin the differences in functional activities between cell populations of the immune system.

Functional clustering and lineage markers: insights into cellular differentiation and gene function from large-scale microarray studies of purified primary cell populations

Very large microarray datasets showing gene expression across multiple tissues and cell populations provide a window on the transcriptional networks that underpin the differences in functional activity between biological systems. Clusters of co-expressed genes provide lineage markers, candidate regulators of cell function and, by applying the principle of guilt by association, candidate functions for genes of currently unknown function. We have analysed a dataset comprising pure cell populations from hemopoietic and non-hemopoietic cell types (http://biogps.gnf.org). Using a novel network visualisation and clustering approach, we demonstrate that it is possible to identify very tight expression signatures associated specifically with embryonic stem cells, mesenchymal cells and hematopoietic lineages. Selected examples validate the prediction that gene function can be inferred by co-expression. One expression cluster was enriched in phagocytes, which, alongside endosome-lysosome constituents, contains genes that may make up a 'pathway' for phagocyte differentiation. Promoters of these genes are enriched for binding sites for the ETS/PU.1 and MITF families. Another cluster was associated with the production of a specific extracellular matrix, with high levels of gene expression shared by cells of mesenchymal origin (fibroblasts, adipocytes, osteoblasts and myoblasts). We discuss the limitations placed upon such data by the presence of alternative promoters with distinct tissue specificity within many protein-coding genes.

Collective motions and specific effectors: a statistical mechanics perspective on biological regulation

BACKGROUND

The interaction of a multiplicity of scales in both time and space is a fundamental feature of biological systems. The complementation of macroscopic (entire organism) and microscopic (molecular biology) views with a mesoscopic level of analysis able to connect the different planes of investigation is urgently needed. This will allow to both obtain a general frame of reference for rationalizing the burden of data coming from high throughput technologies and to derive effective operational views on biological systems.

RESULTS

The network paradigm in which microscopic level elements (nodes) are each other related by functional links so giving rise to both global (entire network) and local (specific) behavior is a promising metaphor to try and develop a statistical mechanics inspired approach for biological systems. Here we show the application of this paradigm to different systems going from yeast metabolism to murine macrophages response to immune stimulation.

CONCLUSIONS

The need to complement the purely molecular view with mesoscopic approaches is evident in all the studied examples that in turn demonstrate the untenability of the simple ergodic approach dominant in molecular biology in which the data coming from huge ensemble of cells are considered as relative to a single 'average' cell.

Heat-induced and spontaneous expression of Hsp70.1Luciferase transgene copies localized on Xp22 in female bovine cells

Background

Expression of several copies of the heat-inducible Hsp70.1Luciferase (LUC) transgene inserted at a single X chromosome locus of a bull (Bos taurus) was assessed in females after X-chromosome inactivation (XCI). Furthermore, impact of the chromosomal environment on the spontaneous expression of these transgene copies before XCI was studied during early development in embryos obtained after in vitro fertilization (IVF), when the locus was carried by the X chromosome inherited from the bull, and after somatic cell nuclear transfer (SCNT) cloning, when the locus could be carried by the inactive Xi or the active Xa chromosome in a female donor cell, or by the (active) X in a male donor cell.

Findings

Transgene copies were mapped to bovine Xp22. In XX^LUC female fibroblasts, i.e. after random XCI, the proportions of late-replicating inactive and early-replicating active X^LUC chromosomes were not biased and the proportion of cells displaying an increase in the level of immunostained luciferase protein after heat-shock induction was similar to that in male fibroblasts. Spontaneous transgene expression occurred at the 8-16-cell stage both in transgenic (female) embryos obtained after IVF and in male and female embryos obtained after SCNT.

Conclusions

The X^LUC chromosome is normally inactivated but at least part of the inactivated X-linked Hsp70.1Luciferase transgene copies remains heat-inducible after random XCI in somatic cells. Before XCI, the profile of the transgenes' spontaneous expression is independent of the epigenetic origin of the X^LUC chromosome since it is similar in IVF female, SCNT male and SCNT female embryos.