A cis-regulatory map of the Drosophila genome.

Systematic annotation of gene regulatory elements is a major challenge in genome science. Direct mapping of chromatin modification marks and transcriptional factor binding sites genome-wide has successfully identified specific subtypes of regulatory elements. In Drosophila several pioneering studies have provided genome-wide identification of Polycomb response elements, chromatin states, transcription factor binding sites, RNA polymerase II regulation and insulator elements; however, comprehensive annotation of the regulatory genome remains a significant challenge. Here we describe results from the modENCODE cis-regulatory annotation project. We produced a map of the Drosophila melanogaster regulatory genome on the basis of more than 300 chromatin immunoprecipitation data sets for eight chromatin features, five histone deacetylases and thirty-eight site-specific transcription factors at different stages of development. Using these data we inferred more than 20,000 candidate regulatory elements and validated a subset of predictions for promoters, enhancers and insulators in vivo. We identified also nearly 2,000 genomic regions of dense transcription factor binding associated with chromatin activity and accessibility. We discovered hundreds of new transcription factor co-binding relationships and defined a transcription factor network with over 800 potential regulatory relationships.

Of elections and cell-death decisions.

A new study in the journal Nature (Spencer et al., 2009) argues that cell-to-cell variation in the decision to undergo apoptosis is not due to genetic, epigenetic, or cell-cycle differences, nor due to random molecular noise, but instead is determined by differences in protein abundances.

A protein turnover signaling motif controls the stimulus-sensitivity of stress response pathways.

Stimulus-induced perturbations from the steady state are a hallmark of signal transduction. In some signaling modules, the steady state is characterized by rapid synthesis and degradation of signaling proteins. Conspicuous among these are the p53 tumor suppressor, its negative regulator Mdm2, and the negative feedback regulator of NFκB, IκBα. We investigated the physiological importance of this turnover, or flux, using a computational method that allows flux to be systematically altered independently of the steady state protein abundances. Applying our method to a prototypical signaling module, we show that flux can precisely control the dynamic response to perturbation. Next, we applied our method to experimentally validated models of p53 and NFκB signaling. We find that high p53 flux is required for oscillations in response to a saturating dose of ionizing radiation (IR). In contrast, high flux of Mdm2 is not required for oscillations but preserves p53 sensitivity to sub-saturating doses of IR. In the NFκB system, degradation of NFκB-bound IκB by the IκB kinase (IKK) is required for activation in response to TNF, while high IKK-independent degradation prevents spurious activation in response to metabolic stress or low doses of TNF. Our work identifies flux pairs with opposing functional effects as a signaling motif that controls the stimulus-sensitivity of the p53 and NFκB stress-response pathways, and may constitute a general design principle in signaling pathways.

Characterizing the relationship between steady state and response using analytical expressions for the steady states of mass action models.

The steady states of cells affect their response to perturbation. Indeed, diagnostic markers for predicting the response to therapeutic perturbation are often based on steady state measurements. In spite of this, no method exists to systematically characterize the relationship between steady state and response. Mathematical models are established tools for studying cellular responses, but characterizing their relationship to the steady state requires that it have a parametric, or analytical, expression. For some models, this expression can be derived by the King-Altman method. However, King-Altman requires that no substrate act as an enzyme, and is therefore not applicable to most models of signal transduction. For this reason we developed py-substitution, a simple but general method for deriving analytical expressions for the steady states of mass action models. Where the King-Altman method is applicable, we show that py-substitution yields an equivalent expression, and at comparable efficiency. We use py-substitution to study the relationship between steady state and sensitivity to the anti-cancer drug candidate, dulanermin (recombinant human TRAIL). First, we use py-substitution to derive an analytical expression for the steady state of a published model of TRAIL-induced apoptosis. Next, we show that the amount of TRAIL required for cell death is sensitive to the steady state concentrations of procaspase 8 and its negative regulator, Bar, but not the other procaspase molecules. This suggests that activation of caspase 8 is a critical point in the death decision process. Finally, we show that changes in the threshold at which TRAIL results in cell death is not always equivalent to changes in the time of death, as is commonly assumed. Our work demonstrates that an analytical expression is a powerful tool for identifying steady state determinants of the cellular response to perturbation. All code is available at http://signalingsystems.ucsd.edu/models-and-code/ or as supplementary material accompanying this paper.

Ex vivo innate immune cytokine signature of enhanced risk of relapsing brucellosis.

BACKGROUND: Brucellosis, a zoonotic infection caused by one of the Gram-negative intracellular bacteria of the Brucella genus, is an ongoing public health problem in Perú. While most patients who receive standard antibiotic treatment recover, 5-40% suffer a brucellosis relapse. In this study, we examined the ex vivo immune cytokine profiles of recovered patients with a history of acute and relapsing brucellosis. METHODOLOGY/PRINCIPAL FINDINGS: Blood was taken from healthy control donors, patients with a history of acute brucellosis, or patients with a history of relapsing brucellosis. Peripheral blood mononuclear cells were isolated and remained in culture without stimulation or were stimulated with a panel of toll-like receptor agonists or heat-killed Brucella melitensis (HKBM) isolates. Innate immune cytokine gene expression and protein secretion were measured by quantitative real-time polymerase chain reaction and a multiplex bead-based immunoassay, respectively. Acute and relapse patients demonstrated consistently elevated cytokine gene expression and secretion levels compared to controls. Notably, these include: basal and stimulus-induced expression of GM-CSF, TNF-α, and IFN-γ in response to LPS and HKBM; basal secretion of IL-6, IL-8, and TNF-α; and HKBM or Rev1-induced secretion of IL-1ß, IL-2, GM-CSF, IFN-Υ, and TNF-α. Although acute and relapse patients were largely indistinguishable by their cytokine gene expression profiles, we identified a robust cytokine secretion signature that accurately discriminates acute from relapse patients. This signature consists of basal IL-6 secretion, IL-1ß, IL-2, and TNF-α secretion in response to LPS and HKBM, and IFN-γ secretion in response to HKBM. CONCLUSIONS/SIGNIFICANCE: This work demonstrates that informative cytokine variations in brucellosis patients can be detected using an ex vivo assay system and used to identify patients with differing infection histories. Targeted diagnosis of this signature may allow for better follow-up care of brucellosis patients through improved identification of patients at risk for relapse.

A framework for modeling the relationship between cellular steady-state and stimulus-responsiveness.

In cell signaling systems, the abundances of signaling molecules are generally thought to determine the response to stimulation. However, the kinetics of molecular processes, for example receptor trafficking and protein turnover, may also play an important role. Few studies have systematically examined this relationship between the resting state and stimulus-responsiveness. Fewer still have investigated the relative contribution of steady-state concentrations and reaction kinetics. Here we describe a mathematical framework for modeling the resting state of signaling systems. Among other things, this framework allows steady-state concentration measurements to be used in parameterizing kinetic models, and enables comprehensive characterization of the relationship between the resting state and the cellular response to stimulation.

Mass spectrometric detection of tissue proteins in plasma.

It has long been thought that blood plasma could serve as a window into the state of one's organs in health and disease because tissue-derived proteins represent a significant fraction of the plasma proteome. Although substantial technical progress has been made toward the goal of comprehensively analyzing the blood plasma proteome, the basic assumption that proteins derived from a variety of tissues could indeed be detectable in plasma using current proteomics technologies has not been rigorously tested. Here we provide evidence that such tissue-derived proteins are both present and detectable in plasma via direct mass spectrometric analysis of captured glycopeptides and thus provide a conceptual basis for plasma protein biomarker discovery and analysis.

UniPep--a database for human N-linked glycosites: a resource for biomarker discovery.

There has been considerable recent interest in proteomic analyses of plasma for the purpose of discovering biomarkers. Profiling N-linked glycopeptides is a particularly promising method because the population of N-linked glycosites represents the proteomes of plasma, the cell surface, and secreted proteins at very low redundancy and provides a compelling link between the tissue and plasma proteomes. Here, we describe UniPep http://www.unipep.org--a database of human N-linked glycosites--as a resource for biomarker discovery.