Thermodynamic State Ensemble Models of cis-Regulation

A major goal in computational biology is to develop models that accurately predict a gene's expression from its surrounding regulatory DNA. Here we present one class of such models, thermodynamic state ensemble models. We describe the biochemical derivation of the thermodynamic framework in simple terms, and lay out the mathematical components that comprise each model. These components include (1) the possible states of a promoter, where a state is defined as a particular arrangement of transcription factors bound to a DNA promoter, (2) the binding constants that describe the affinity of the protein–protein and protein–DNA interactions that occur in each state, and (3) whether each state is capable of transcribing. Using these components, we demonstrate how to compute a cis-regulatory function that encodes the probability of a promoter being active. Our intention is to provide enough detail so that readers with little background in thermodynamics can compose their own cis-regulatory functions. To facilitate this goal, we also describe a matrix form of the model that can be easily coded in any programming language. This formalism has great flexibility, which we show by illustrating how phenomena such as competition between transcription factors and cooperativity are readily incorporated into these models. Using this framework, we also demonstrate that Michaelis-like functions, another class of cis-regulatory models, are a subset of the thermodynamic framework with specific assumptions. By recasting Michaelis-like functions as thermodynamic functions, we emphasize the relationship between these models and delineate the specific circumstances representable by each approach. Application of thermodynamic state ensemble models is likely to be an important tool in unraveling the physical basis of combinatorial cis-regulation and in generating formalisms that accurately predict gene expression from DNA sequence

Genome-Wide Identification of Bcl11b Gene Targets Reveals Role in Brain-Derived Neurotrophic Factor Signaling

B-cell leukemia/lymphoma 11B (Bcl11b) is a transcription factor showing predominant expression in the striatum. To date, there are no known gene targets of Bcl11b in the nervous system. Here, we define targets for Bcl11b in striatal cells by performing chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) in combination with genome-wide expression profiling. Transcriptome-wide analysis revealed that 694 genes were significantly altered in striatal cells over-expressing Bcl11b, including genes showing striatal-enriched expression similar to Bcl11b. ChIP-seq analysis demonstrated that Bcl11b bound a mixture of coding and non-coding sequences that were within 10 kb of the transcription start site of an annotated gene. Integrating all ChIP-seq hits with the microarray expression data, 248 direct targets of Bcl11b were identified. Functional analysis on the integrated gene target list identified several zinc-finger encoding genes as Bcl11b targets, and further revealed a significant association of Bcl11b to brain-derived neurotrophic factor/neurotrophin signaling. Analysis of ChIP-seq binding regions revealed significant consensus DNA binding motifs for Bcl11b. These data implicate Bcl11b as a novel regulator of the BDNF signaling pathway, which is disrupted in many neurological disorders. Specific targeting of the Bcl11b-DNA interaction could represent a novel therapeutic approach to lowering BDNF signaling specifically in striatal cells

Generation of the first BAC-based physical map of the common carp genome

<p>Abstract</p> <p>Background</p> <p>Common carp (<it>Cyprinus carpio</it>), a member of Cyprinidae, is the third most important aquaculture species in the world with an annual global production of 3.4 million metric tons, accounting for nearly 14% of the all freshwater aquaculture production in the world. Apparently genomic resources are needed for this species in order to study its performance and production traits. In spite of much progress, no physical maps have been available for common carp. The objective of this project was to generate a BAC-based physical map using fluorescent restriction fingerprinting.</p> <p>Result</p> <p>The first generation of common carp physical map was constructed using four- color High Information Content Fingerprinting (HICF). A total of 72,158 BAC clones were analyzed that generated 67,493 valid fingerprints (5.5 × genome coverage). These BAC clones were assembled into 3,696 contigs with the average length of 476 kb and a N50 length of 688 kb, representing approximately 1.76 Gb of the common carp genome. The largest contig contained 171 BAC clones with the physical length of 3.12 Mb. There are 761 contigs longer than the N50, and these contigs should be the most useful resource for future integrations with linkage map and whole genome sequence assembly. The common carp physical map is available at <url>http://genomics.cafs.ac.cn/fpc/WebAGCoL/Carp/WebFPC/</url>.</p> <p>Conclusion</p> <p>The reported common carp physical map is the first physical map of the common carp genome. It should be a valuable genome resource facilitating whole genome sequence assembly and characterization of position-based genes important for aquaculture traits.</p

Interrogating domain-domain interactions with parsimony based approaches

<p>Abstract</p> <p>Background</p> <p>The identification and characterization of interacting domain pairs is an important step towards understanding protein interactions. In the last few years, several methods to predict domain interactions have been proposed. Understanding the power and the limitations of these methods is key to the development of improved approaches and better understanding of the nature of these interactions.</p> <p>Results</p> <p>Building on the previously published Parsimonious Explanation method (PE) to predict domain-domain interactions, we introduced a new Generalized Parsimonious Explanation (GPE) method, which (i) adjusts the granularity of the domain definition to the granularity of the input data set and (ii) permits domain interactions to have different costs. This allowed for preferential selection of the so-called "co-occurring domains" as possible mediators of interactions between proteins. The performance of both variants of the parsimony method are competitive to the performance of the top algorithms for this problem even though parsimony methods use less information than some of the other methods. We also examined possible enrichment of co-occurring domains and homo-domains among domain interactions mediating the interaction of proteins in the network. The corresponding study was performed by surveying domain interactions predicted by the GPE method as well as by using a combinatorial counting approach independent of any prediction method. Our findings indicate that, while there is a considerable propensity towards these special domain pairs among predicted domain interactions, this overrepresentation is significantly lower than in the iPfam dataset.</p> <p>Conclusion</p> <p>The Generalized Parsimonious Explanation approach provides a new means to predict and study domain-domain interactions. We showed that, under the assumption that all protein interactions in the network are mediated by domain interactions, there exists a significant deviation of the properties of domain interactions mediating interactions in the network from that of iPfam data.</p

Cinnamomum cassia Bark in Two Herbal Formulas Increases Life Span in Caenorhabditis elegans via Insulin Signaling and Stress Response Pathways

Background: Proving the efficacy and corresponding mode of action of herbal supplements is a difficult challenge for evidence-based herbal therapy. A major hurdle is the complexity of herbal preparations, many of which combine multiple herbs, particularly when the combination is assumed to be vitally important to the effectiveness of the herbal therapy. This issue may be addressed through the use of contemporary methodology and validated animal models. Methods and Principal Findings: In this study, two commonly used traditional herbal formulas, Shi Quan Da Bu Tang (SQDB) and Huo Luo Xiao Ling Dan (HLXL) were evaluated using a survival assay and oxidative stress biomarkers in a well-established C. elegans model of aging. HLXL is an eleven herb formula modified from a top-selling traditional herbal formula for the treatment of arthritic joint pain. SQDB consists of ten herbs often used for fatigue and energy, particularly in the aged. We demonstrate here that SQDB significantly extend life span in a C. elegans model of aging. Among all individual herbs tested, two herbs Cinnamomum cassia bark (Chinese pharmaceutical name: Cinnamomi Cortex, CIN) and Panax ginseng root (Chinese pharmaceutical name: Ginseng Radix, GS) significantly extended life span in C. elegans. CIN in both SQDB and HLXL formula extended life span via modulation of multiple longevity assurance genes, including genes involved in insulin signaling and stress response pathways. All the life-span-extending herbs (SQDB, CIN and GS) also attenuated levels of H2O2 and enhanced small heat shock protein expression. Furthermore, the life spanextending herbs significantly delayed human amyloid beta (Aβ)-induced toxicity in transgenic C. elegans expressing human Aβ. Conclusion/Significance:These results validate an invertebrate model for rapid, systematic evaluation of commonly used Chinese herbal formulations and may provide insight for designing future evidence-based herbal therapy(s). Copyright: © 2010 Yu et al.published_or_final_versio

VILIP-1 Downregulation in Non-Small Cell Lung Carcinomas: Mechanisms and Prediction of Survival

VILIP-1, a member of the neuronal Ca++ sensor protein family, acts as a tumor suppressor gene in an experimental animal model by inhibiting cell proliferation, adhesion and invasiveness of squamous cell carcinoma cells. Western Blot analysis of human tumor cells showed that VILIP-1 expression was undetectable in several types of human tumor cells, including 11 out of 12 non-small cell lung carcinoma (NSCLC) cell lines. The down-regulation of VILIP-1 was due to loss of VILIP-1 mRNA transcripts. Rearrangements, large gene deletions or mutations were not found. Hypermethylation of the VILIP-1 promoter played an important role in gene silencing. In most VILIP-1-silent cells the VILIP-1 promoter was methylated. In vitro methylation of the VILIP-1 promoter reduced its activity in a promoter-reporter assay. Transcriptional activity of endogenous VILIP-1 promoter was recovered by treatment with 5′-aza-2′-deoxycytidine (5′-Aza-dC). Trichostatin A (TSA), a histone deacetylase inhibitor, potently induced VILIP-1 expression, indicating that histone deacetylation is an additional mechanism of VILIP-1 silencing. TSA increased histone H3 and H4 acetylation in the region of the VILIP-1 promoter. Furthermore, statistical analysis of expression and promoter methylation (n = 150 primary NSCLC samples) showed a significant relationship between promoter methylation and protein expression downregulation as well as between survival and decreased or absent VILIP-1 expression in lung cancer tissues (p<0.0001). VILIP-1 expression is silenced by promoter hypermethylation and histone deacetylation in aggressive NSCLC cell lines and primary tumors and its clinical evaluation could have a role as a predictor of short-term survival in lung cancer patients

SARS Coronavirus 3b Accessory Protein Modulates Transcriptional Activity of RUNX1b

BACKGROUND: The causative agent of severe acute respiratory syndrome, SARS coronavirus (SARS-CoV) genome encodes several unique group specific accessory proteins with unknown functions. Among them, accessory protein 3b (also known as ORF4) was lately identified as one of the viral interferon antagonist. Recently our lab uncovered a new role for 3b in upregulation of AP-1 transcriptional activity and its downstream genes. Thus, we believe that 3b might play an important role in SARS-CoV pathogenesis and therefore is of considerable interest. The current study aims at identifying novel host cellular interactors of the 3b protein. METHODOLOGY/PRINCIPAL FINDINGS: In this study, using yeast two-hybrid and co-immunoprecipitation techniques, we have identified a host transcription factor RUNX1b (Runt related transcription factor, isoform b) as a novel interacting partner for SARS-CoV 3b protein. Chromatin immunoprecipitaion (ChIP) and reporter gene assays in 3b expressing jurkat cells showed recruitment of 3b on the RUNX1 binding element that led to an increase in RUNX1b transactivation potential on the IL2 promoter. Kinase assay and pharmacological inhibitor treatment implied that 3b also affect RUNX1b transcriptional activity by regulating its ERK dependent phosphorylation levels. Additionally, mRNA levels of MIP-1α, a RUNX1b target gene upregulated in SARS-CoV infected monocyte-derived dendritic cells, were found to be elevated in 3b expressing U937 monocyte cells. CONCLUSIONS/SIGNIFICANCE: These results unveil a novel interaction of SARS-CoV 3b with the host factor, RUNX1b, and speculate its physiological relevance in upregulating cytokines and chemokine levels in state of SARS virus infection

Peregrine and saker falcon genome sequences provide insights into evolution of a predatory lifestyle

As top predators, falcons possess unique morphological, physiological and behavioral adaptations that allow them to be successful hunters: for example, the peregrine is renowned as the world's fastest animal. To examine the evolutionary basis of predatory adaptations, we sequenced the genomes of both the peregrine (Falco peregrinus) and saker falcon (Falco cherrug), and we present parallel, genome-wide evidence for evolutionary innovation and selection for a predatory lifestyle. The genomes, assembled using Illumina deep sequencing with greater than 100-fold coverage, are both approximately 1.2 Gb in length, with transcriptome-assisted prediction of approximately 16,200 genes for both species. Analysis of 8,424 orthologs in both falcons, chicken, zebra finch and turkey identified consistent evidence for genome-wide rapid evolution in these raptors. SNP-based inference showed contrasting recent demographic trajectories for the two falcons, and gene-based analysis highlighted falcon-specific evolutionary novelties for beak development and olfaction and specifically for homeostasis-related genes in the arid environment–adapted saker