Combinatorial interactions among transcription factors are critical to directing tissue-specific gene expression. To build a global atlas of these combinations, we have screened for physical interactions among the majority of human and mouse DNA-binding transcription factors (TFs). The complete networks contain 762 human and 877 mouse interactions. Analysis of the networks reveals that highly connected TFs are broadly expressed across tissues, and that roughly half of the measured interactions are conserved between mouse and human. The data highlight the importance of TF combinations for determining cell fate, and they lead to the identification of a SMAD3/FLI1 complex expressed during development of immunity. The availability of large TF combinatorial networks in both human and mouse will provide many opportunities to study gene regulation, tissue differentiation, and mammalian evolution.

We present a sound and complete graphical criterion for reading dependencies from the minimal undirected independence map G of a graphoid M that satisfies weak transitivity. Here, complete means that it is able to read all the dependencies in M that can be derived by applying the graphoid properties and weak transitivity to the dependencies used in the construction of G and the independencies obtained from G by vertex separation. We argue that assuming weak transitivity is not too restrictive. As an intermediate step in the derivation of the graphical criterion, we prove that for any undirected graph G there exists a strictly positive discrete probability distribution with the prescribed sample spaces that is faithful to G. We also report an algorithm that implements the graphical criterion and whose running time is considered to be at most O(n(2)(e + n)) for n nodes and e edges. Finally, we illustrate how the graphical criterion can be used within bioinformatics to identify biologically meaningful gene dependencies.

Complex regulatory dynamics is ubiquitous in molecular networks composed of genes and proteins. Recent progress in computational biology and its application to molecular data generate a growing number of complex networks. Yet, it has been difficult to understand the governing principles of these networks beyond graphical analysis or extensive numerical simulations. Here the authors exploit several simplifying biological circumstances which thereby enable to directly detect the underlying dynamical regularities driving periodic oscillations in a dynamical nonlinear computational model of a protein-protein network. System analysis is performed using the cell cycle, a mathematically well-described complex regulatory circuit driven by external signals. By introducing an explicit time delay and using a tearing-and-zooming approach the authors reduce the system to a piecewise linear system with two variables that capture the dynamics of this complex network. A key step in the analysis is the identification of functional subsystems by identifying the relations between state-variables within the model. These functional subsystems are referred to as dynamical modules operating as sensitive switches in the original complex model. By using reduced mathematical representations of the subsystems the authors derive explicit conditions on how the cell cycle dynamics depends on system parameters, and can, for the first time, analyse and prove global conditions for system stability. The approach which includes utilising biological simplifying conditions, identification of dynamical modules and mathematical reduction of the model complexity may be applicable to other well-characterised biological regulatory circuits.

Recently, important insights into static network topology for biological systems have been obtained, but still global dynamical network properties determining stability and system responsiveness have not been accessible for analysis. Herein, we explore a genome-wide gene-to-gene regulatory network based on expression data from the cell cycle in Saccharomyces cerevisae (budding yeast). We recover static properties like hubs (genes having several out-going connections), network motifs and modules, which have previously been derived from multiple data sources such as whole-genome expression measurements, literature mining, protein-protein and transcription factor binding data. Further, our analysis uncovers some novel dynamical design principles; hubs are both repressed and repressors, and the intra-modular dynamics are either strongly activating or repressing whereas inter-modular couplings are weak. Finally, taking advantage of the inferred strength and direction of all interactions, we perform a global dynamical systems analysis of the network. Our inferred dynamics of hubs, motifs and modules produce a more stable network than what is expected given randomised versions. The main contribution of the repressed hubs is to increase system stability, while higher order dynamic effects (e.g. module dynamics) mainly increase system flexibility. Altogether, the presence of hubs, motifs and modules induce few flexible modes, to which the network is extra sensitive to an external signal. We believe that our approach, and the inferred biological mode of strong flexibility and stability, will also apply to other cellular networks and adaptive systems.

Working memory capacity, the maximum number of items that we can transiently store in working memory, is a good predictor of our general cognitive abilities. Neural activity in both dorsolateral prefrontal cortex and posterior parietal cortex has been associated with memory retention during visuospatial working memory tasks. The parietal cortex is thought to store the memories. However, the role of the dorsolateral prefrontal cortex, a top-down control area, during pure information retention is debated, and the mechanisms regulating capacity are unknown. Here, we propose that a major role of the dorsolateral prefrontal cortex in working memory is to boost parietal memory capacity. Furthermore, we formulate the boosting mechanism computationally in a biophysical cortical microcircuit model and derive a simple, explicit mathematical formula relating memory capacity to prefrontal and parietal model parameters. For physiologically realistic parameter values, lateral inhibition in the parietal cortex limits mnemonic capacity to a maximum of 2-7 items. However, at high loads inhibition can be counteracted by excitatory prefrontal input, thus boosting parietal capacity. Predictions from the model were confirmed in an fMRI study. Our results show that although memories are stored in the parietal cortex, interindividual differences in memory capacity are partly determined by the strength of prefrontal top-down control. The model provides a mechanistic framework for understanding top-down control of working memory and specifies two different contributions of prefrontal and parietal cortex to working memory capacity.

Environmental exposures filtered through the genetic make-up of each individual alter the transcriptional repertoire in organs central to metabolic homeostasis, thereby affecting arterial lipid accumulation, inflammation, and the development of coronary artery disease (CAD). The primary aim of the Stockholm Atherosclerosis Gene Expression (STAGE) study was to determine whether there are functionally associated genes (rather than individual genes) important for CAD development. To this end, two-way clustering was used on 278 transcriptional profiles of liver, skeletal muscle, and visceral fat (n=66/tissue) and atherosclerotic and unaffected arterial wall (n=40/tissue) isolated from CAD patients during coronary artery bypass surgery. The first step, across all mRNA signals (n=15,042/12,621 RefSeqs/genes) in each tissue, resulted in a total of 60 tissue clusters (n=3958 genes). In the second step (performed within tissue clusters), one atherosclerotic lesion (n=49/48) and one visceral fat (n=59) cluster segregated the patients into two groups that differed in the extent of coronary stenosis (P=0.008 and P=0.00015). The associations of these clusters with coronary atherosclerosis were validated by analyzing carotid atherosclerosis expression profiles. Remarkably, in one cluster (n=55/54) relating to carotid stenosis (P=0.04), 27 genes in the two clusters relating to coronary stenosis were confirmed (n=16/17, P<10^-27and-30). Genes in the transendothelial migration of leukocytes (TEML) pathway were overrepresented in all three clusters, referred to as the atherosclerosis module (A-module). In a second validation step, using three independent cohorts, the A-module was found to be genetically enriched with CAD risk by 1.8-fold (P<0.004). The transcription co-factor LIM domain binding 2 (LDB2) was identified as a potential high-hierarchy regulator of the A-module, a notion supported by subnetwork analysis, cellular and lesion expression of LDB2, and the expression of 13 TEML genes in Ldb2-deficient arterial wall. Thus, the A-module appears to be important for atherosclerosis development and together with LDB2 merits further attention in CAD research.

Background: Plasmid encoded (CTX)-C-bla-M enzymes represent an important sub-group of class A beta-lactamases causing the ESBL phenotype which is increasingly found in Enterobacteriaceae including Klebsiella spp. Molecular typing of clinical ESBL-isolates has become more and more important for prevention of the dissemination of ESBL-producers among nosocomial environment.

Methods: Multiple displacement amplified DNA derived from 20 K. pneumoniae and 34 K. oxytoca clinical isolates with an ESBL-phenotype was used in a universal CTX-M PCR amplification assay. Identification and differentiation of (CTX)-C-bla-M and (OXY)-O-bla/K1 sequences was obtained by DNA sequencing of M13-sequence-tagged CTX-M PCR-amplicons using a M13-specific sequencing primer.

Results: Nine out of 20 K. pneumoniae clinical isolates had a (CTX)-C-bla-M genotype. Interestingly, we found that the universal degenerated primers also amplified the chromosomally located K1-gene in all 34 K. oxytoca clinical isolates. Molecular identification and differentiation between (CTX)-C-bla-M and (OXY)-O-bla/K1-genes could only been achieved by sequencing of the PCR-amplicons. In silico analysis revealed that the universal degenerated CTX-M primer-pair used here might also amplify the chromosomally located (OXY)-O-bla and K1-genes in Klebsiella spp. and K1-like genes in other Enterobacteriaceae.

Conclusion: The PCR-based molecular typing method described here enables a rapid and reliable molecular identification of (CTX)-C-bla-M, and (OXY)-O-bla/K1-genes. The principles used in this study could also be applied to any situation in which antimicrobial resistance genes would need to be sequenced.

The quest to determine cause from effect is often referred to as reverse engineering in the context of cellular networks. Here we propose and evaluate an algorithm for reverse engineering a gene regulatory network from time-series kind steady-state data. Our algorithmic pipeline, which is rather standard in its parts but not in its integrative composition, combines ordinary differential equations, parameter estimations by least angle regression, and cross-validation procedures for determining the in-degrees and selection of nonlinear transfer functions. The result of the algorithm is a complete directed net-work, in which each edge has been assigned a score front it bootstrap procedure. To evaluate the performance, we submitted the outcome of the algorithm to the reverse engineering assessment competition DREAM2, where we used the data corresponding to the InSillico1 and InSilico2 networks as input. Our algorithm outperformed all other algorithms when inferring one of the directed gene-to-gene networks.

Using deep sequencing (deepCAGE), the FANTOM4 study measured the genome-wide dynamics of transcription-start-site usage in the human monocytic cell line THP-1 throughout a time course of growth arrest and differentiation. Modeling the expression dynamics in terms of predicted cis-regulatory sites, we identified the key transcription regulators, their time-dependent activities and target genes. Systematic siRNA knockdown of 52 transcription factors confirmed the roles of individual factors in the regulatory network. Our results indicate that cellular states are constrained by complex networks involving both positive and negative regulatory interactions among substantial numbers of transcription factors and that no single transcription factor is both necessary and sufficient to drive the differentiation process.

Background: Free fatty acids released from adipose tissue affect the synthesis of apolipoprotein B-containing lipoproteins and glucose metabolism in the liver. Whether there also exists a reciprocal metabolic arm affecting energy metabolism in white adipose tissue is unknown. Methods and Findings: We investigated the effects of apoB-containing lipoproteins on catecholamine-induced lipolysis in adipocytes from subcutaneous fat cells of obese but otherwise healthy men, fat pads from mice with plasma lipoproteins containing high or intermediate levels of apoB100 or no apoB100, primary cultured adipocytes, and 3T3-L1 cells. In subcutaneous fat cells, the rate of lipolysis was inversely related to plasma apoB levels. In human primary adipocytes, LDL inhibited lipolysis in a concentration-dependent fashion. In contrast, VLDL had no effect. Lipolysis was increased in fat pads from mice lacking plasma apoB100, reduced in apoB100-only mice, and intermediate in wild-type mice. Mice lacking apoB100 also had higher oxygen consumption and lipid oxidation. In 3T3-L1 cells, apoB100-containing lipoproteins inhibited lipolysis in a dose-dependent fashion, but lipoproteins containing apoB48 had no effect. ApoB100-LDL mediated inhibition of lipolysis was abolished in fat pads of mice deficient in the LDL receptor (Ldlr-/- Apob100/100). Conclusions: Our results show that the binding of apoB100-LDL to adipocytes via the LDL receptor inhibits intracellular noradrenaline-induced lipolysis in adipocytes. Thus, apoB100-LDL is a novel signaling molecule from the liver to peripheral fat deposits that may be an important link between atherogenic dyslipidemias and facets of the metabolic syndrome. © 2008 Skogsberg et al.