Search

The term "transcriptional network" refers to the mechanism(s) that underlies coordinated expression of genes, typically involving transcription factors (TFs) binding to the promoters of multiple genes, and individual genes controlled by multiple TFs. A multitude of studies in the last two decades have aimed to map and characterize transcriptional networks in the yeast Saccharomyces cerevisiae. We review the methodologies and accomplishments of these studies, as well as challenges we now face. For most yeast TFs, data have been collected on their sequence preferences, in vivo promoter occupancy, and gene expression profiles in deletion mutants. These systematic studies have led to the identification of new regulators of numerous cellular functions and shed light on the overall organization of yeast gene regulation. However, many yeast TFs appear to be inactive under standard laboratory growth conditions, and many of the available data were collected using techniques that have since been improved. Perhaps as a consequence, comprehensive and accurate mapping among TF sequence preferences, promoter binding, and gene expression remains an open challenge. We propose that the time is ripe for renewed systematic efforts toward a complete mapping of yeast transcriptional regulatory mechanisms.

The complement system is a potent component of the innate immune response, promoting inflammation and orchestrating defense against pathogens. However, dysregulation of complement is critical to several autoimmune and inflammatory syndromes. Elevated expression of the proinflammatory cytokine IL-1β is often linked to such diseases. In this study, we reveal the mechanistic link between complement and IL-1β secretion using murine dendritic cells. IL-1β secretion occurs following intracellular caspase-1 activation by inflammasomes. We show that complement elicits secretion of both IL-1β and IL-18 in vitro and in vivo via the NLRP3 inflammasome. This effect depends on the inflammasome components NLRP3 and ASC, as well as caspase-1 activity. Interestingly, sublethal complement membrane attack complex formation, but not the anaphylatoxins C3a and C5a, activated the NLRP3 inflammasome in vivo. These findings provide insight into the molecular processes underlying complement-mediated inflammation and highlight the possibility of targeting IL-1β to control complement-induced disease and pathological inflammation.

Gene duplication results in two identical paralogs that diverge through mutation, leading to loss or gain of interactions with other biomolecules. Here, we comprehensively characterize such network rewiring for C. elegans transcription factors (TFs) within and across four newly delineated molecular networks. Remarkably, we find that even highly similar TFs often have different interaction degrees and partners. In addition, we find that most TF families have a member that is highly connected in multiple networks. Further, different TF families have opposing correlations between network connectivity and phylogenetic age, suggesting that they are subject to different evolutionary pressures. Finally, TFs that have similar partners in one network generally do not in another, indicating a lack of pressure to retain cross-network similarity. Our multiparameter analyses provide unique insights into the evolutionary dynamics that shaped TF networks., (Copyright © 2013 Elsevier Inc. All rights reserved.)

RNA-binding proteins are key regulators of gene expression, yet only a small fraction have been functionally characterized. Here we report a systematic analysis of the RNA motifs recognized by RNA-binding proteins, encompassing 205 distinct genes from 24 diverse eukaryotes. The sequence specificities of RNA-binding proteins display deep evolutionary conservation, and the recognition preferences for a large fraction of metazoan RNA-binding proteins can thus be inferred from their RNA-binding domain sequence. The motifs that we identify in vitro correlate well with in vivo RNA-binding data. Moreover, we can associate them with distinct functional roles in diverse types of post-transcriptional regulation, enabling new insights into the functions of RNA-binding proteins both in normal physiology and in human disease. These data provide an unprecedented overview of RNA-binding proteins and their targets, and constitute an invaluable resource for determining post-transcriptional regulatory mechanisms in eukaryotes.

The membrane attack complex of complement (MAC), apart from its classical role of lysing cells, can also trigger a range of non-lethal effects on cells, acting as a drive to inflammation. In the present study, we chose to investigate these non-lethal effects on inflammasome activation. We found that, following sublytic MAC attack, there is increased cytosolic Ca(2+) concentration, at least partly through Ca(2+) release from the endoplasmic reticulum lumen via the inositol 1,4,5-triphosphate receptor (IP3R) and ryanodine receptor (RyR) channels. This increase in intracellular Ca(2+) concentration leads to Ca(2+) accumulation in the mitochondrial matrix via the 'mitochondrial calcium uniporter' (MCU), and loss of mitochondrial transmembrane potential, triggering NLRP3 inflammasome activation and IL-1β release. NLRP3 co-localises with the mitochondria, probably sensing the increase in calcium and the resultant mitochondrial dysfunction, leading to caspase activation and apoptosis. This is the first study that links non-lethal effects of sublytic MAC attack with inflammasome activation and provides a mechanism by which sublytic MAC can drive inflammation and apoptosis.

To evaluate the effect of CG methylation on DNA binding of sequence-specific B-ZIP transcription factors (TFs) in a high-throughput manner, we enzymatically methylated the cytosine in the CG dinucleotide on protein binding microarrays. Two Agilent DNA array designs were used. One contained 40,000 features using de Bruijn sequences where each 8-mer occurs 32 times in various positions in the DNA sequence. The second contained 180,000 features with each CG containing 8-mer occurring three times. The first design was better for identification of binding motifs, while the second was better for quantification. Using this novel technology, we show that CG methylation enhanced binding for CEBPA and CEBPB and inhibited binding for CREB, ATF4, JUN, JUND, CEBPD, and CEBPG. The CEBPB|ATF4 heterodimer bound a novel motif CGAT|GCAA 10-fold better when methylated. The electrophoretic mobility shift assay (EMSA) confirmed these results. CEBPB ChIP-seq data using primary female mouse dermal fibroblasts with 50× methylome coverage for each strand indicate that the methylated sequences well-bound on the arrays are also bound in vivo. CEBPB bound 39% of the methylated canonical 10-mers ATTGC|GCAAT in the mouse genome. After ATF4 protein induction by thapsigargin which results in ER stress, CEBPB binds methylated CGAT|GCAA in vivo, recapitulating what was observed on the arrays. This methodology can be used to identify new methylated DNA sequences preferentially bound by TFs, which may be functional in vivo.

Nucleosomes in all eukaryotes examined to date adopt a characteristic architecture within genes and play fundamental roles in regulating transcription, yet the identity and precise roles of many of the trans-acting factors responsible for the establishment and maintenance of this organization remain to be identified. We profiled a compendium of 50 yeast strains carrying conditional alleles or complete deletions of genes involved in transcriptional regulation, histone biology, and chromatin remodeling, as well as compounds that target transcription and histone deacetylases, to assess their respective roles in nucleosome positioning and transcription. We find that nucleosome patterning in genes is affected by many factors, including the CAF-1 complex, Spt10, and Spt21, in addition to previously reported remodeler ATPases and histone chaperones. Disruption of these factors or reductions in histone levels led genic nucleosomes to assume positions more consistent with their intrinsic sequence preferences, with pronounced and specific shifts of the +1 nucleosome relative to the transcription start site. These shifts of +1 nucleosomes appear to have functional consequences, as several affected genes in Ino80 mutants exhibited altered expression responses. Our parallel expression profiling compendium revealed extensive transcription changes in intergenic and antisense regions, most of which occur in regions with altered nucleosome occupancy and positioning. We show that the nucleosome-excluding transcription factors Reb1, Abf1, Tbf1, and Rsc3 suppress cryptic transcripts at their target promoters, while a combined analysis of nucleosome and expression profiles identified 36 novel transcripts that are normally repressed by Tup1/Cyc8. Our data confirm and extend the roles of chromatin remodelers and chaperones as major determinants of genic nucleosome positioning, and these data provide a valuable resource for future studies., Competing Interests: The authors have declared that no competing interests exist.

Background: Research on complement factor H (fH) in human disease is hampered by lack of an assay suitable for use in large-scale epidemiological studies. We describe the development and validation of a high throughput nephelometric assay for fH., Methods: Reagents from a commercial radial immunodiffusion (RID) assay (The Binding Site) were adapted for use on the Siemens BNII high throughput nephelometric instrument. The assay was calibrated with a highly purified human fH preparation with rigorously determined concentration, and assay performance was comprehensively evaluated using samples from healthy human volunteers, with the commercial RID assay as a comparator. The distribution and determinants of circulating fH concentration in humans were then investigated in a large representative population sample., Results: The nephelometric assay had recovery close to 100%, was reproducible with intra- and inter-assay CV's of 11% and 5-15% respectively, and had a wider operating range than the RID assay. fH values were unaffected after multiple freeze-thaw cycles demonstrating that it is evidently a stable analyte for immunoassay. fH concentration was unaltered by an acute inflammatory stimulus. The population study showed that plasma fH concentration is associated with circulating lipids and indices of body fat., Conclusion: We present the first high throughput assay for circulating fH; the assay is accurate and reliable with reproducible measures from stored samples. It has established the distribution of fH values at a population level and demonstrated important associations with circulating lipids and indices of body fat, thus providing an important reference for future clinical and epidemiological investigations., (Copyright © 2013 Elsevier B.V. All rights reserved.)

Genomic analyses often involve scanning for potential transcription factor (TF) binding sites using models of the sequence specificity of DNA binding proteins. Many approaches have been developed to model and learn a protein's DNA-binding specificity, but these methods have not been systematically compared. Here we applied 26 such approaches to in vitro protein binding microarray data for 66 mouse TFs belonging to various families. For nine TFs, we also scored the resulting motif models on in vivo data, and found that the best in vitro-derived motifs performed similarly to motifs derived from the in vivo data. Our results indicate that simple models based on mononucleotide position weight matrices trained by the best methods perform similarly to more complex models for most TFs examined, but fall short in specific cases (<10% of the TFs examined here). In addition, the best-performing motifs typically have relatively low information content, consistent with widespread degeneracy in eukaryotic TF sequence preferences.

Large-scale genome-wide SNP association studies have identified an association between variants of CR1, the gene encoding complement component receptor 1, and the sporadic form of Alzheimer's disease. The role of CR1 and the complement system in Alzheimer's disease remains far from clear. In rodents the closest ortholog of CR1 is the Crry gene (Cr1-related protein Y). To begin to explore its role in Alzheimer's disease we examined hippocampal lysates from Crry(-/-) mice and age matched controls by immunoblotting. We measured complement factor H, a component of the complement system and biomarker for Alzheimer's disease progression, and tau phosphorylation at the serine 235 site, hyperphosphorylated forms of tau being a defining neuropathological hallmark of the disease. We found that levels of CFH and of tau phosphorylation at serine 235 were strongly and significantly reduced in Crry(-/-) samples. These observations provide a starting point for further attempts to determine the role of CR1 in the neuropathological process driving Alzheimer's disease., (Copyright © 2012 Elsevier Ireland Ltd. All rights reserved.)

Trafficking of lung dendritic cells (DCs) to the draining lymph node (dLN) is a crucial step for the initiation of T cell responses upon pathogen challenge. However, little is known about the factors that regulate lung DC migration to the dLN. In this study, using a model of influenza infection, we demonstrate that complement component C3 is critically required for efficient emigration of DCs from the lung to the dLN. C3 deficiency affect lung DC-mediated viral antigen transport to the dLN, resulting in severely compromised priming of virus-specific T cell responses. Consequently, C3-deficient mice lack effector T cell response in the lungs that affected viral clearance and survival. We further show that direct signaling by C3a and C5a through C3aR and C5aR respectively expressed on lung DCs is required for their efficient trafficking. However, among lung DCs, only CD103(+) DCs make a significant contribution to lung C5a levels and exclusively produce high levels of C3 and C5 during influenza infection. Collectively, our findings show that complement has a profound impact on immune regulation by controlling tissue DC trafficking and highlights a potential utility for complement as an adjuvant in novel vaccine strategies.

Microglial priming predisposes the brain to neurodegeneration and affects disease progression. The signal to switch from the quiescent to the primed state is unknown. We show that deleting the C3 convertase regulator complement receptor 1-related protein y (Crry) induces microglial priming. Mice that were double-knockout for Crry and either C3 or factor B did not show priming, demonstrating dependence on alternative pathway activation. Colocalization of C3b/iC3b and CR3 implicated the CR3/iC3b interaction in priming. Systemic lipopolysaccharide challenge overactivated primed microglia with florid expression of proinflammatory molecules, which were blocked by complement inhibition. Relevance for neurodegenerative disease is exemplified by human multiple sclerosis (MS) and by experimental autoimmune encephalomyelitis (EAE), a model of MS. In human MS, microglial priming was evident in perilesional white matter, in close proximity to C3b/iC3b deposits. EAE was accelerated and exacerbated in Crry-deficient mice, and was dependent on C activation. In summary, C3-dependent microglial priming confers susceptibility to other challenges. Our observations are relevant to progression in MS and other neurological diseases exacerbated by acute insults.

The yeast Saccharomyces cerevisiae is a prevalent system for the analysis of transcriptional networks. As a result, multiple DNA-binding sequence specificities (motifs) have been derived for most yeast transcription factors (TFs). However, motifs from different studies are often inconsistent with each other, making subsequent analyses complicated and confusing. Here, we have created YeTFaSCo (The Yeast Transcription Factor Specificity Compendium, http://yetfasco.ccbr.utoronto.ca/), an extensive collection of S. cerevisiae TF specificities. YeTFaSCo differs from related databases by being more comprehensive (including 1709 motifs for 256 proteins or protein complexes), and by evaluating the motifs using multiple objective quality metrics. The metrics include correlation between motif matches and ChIP-chip data, gene expression patterns, and GO terms, as well as motif agreement between different studies. YeTFaSCo also features an index of 'expert-curated' motifs, each associated with a confidence assessment. In addition, the database website features tools for motif analysis, including a sequence scanning function and precomputed genome-browser tracks of motif occurrences across the entire yeast genome. Users can also search the database for motifs that are similar to a query motif.

Background: Cannabis sativa has been cultivated throughout human history as a source of fiber, oil and food, and for its medicinal and intoxicating properties. Selective breeding has produced cannabis plants for specific uses, including high-potency marijuana strains and hemp cultivars for fiber and seed production. The molecular biology underlying cannabinoid biosynthesis and other traits of interest is largely unexplored., Results: We sequenced genomic DNA and RNA from the marijuana strain Purple Kush using shortread approaches. We report a draft haploid genome sequence of 534 Mb and a transcriptome of 30,000 genes. Comparison of the transcriptome of Purple Kush with that of the hemp cultivar 'Finola' revealed that many genes encoding proteins involved in cannabinoid and precursor pathways are more highly expressed in Purple Kush than in 'Finola'. The exclusive occurrence of Δ9-tetrahydrocannabinolic acid synthase in the Purple Kush transcriptome, and its replacement by cannabidiolic acid synthase in 'Finola', may explain why the psychoactive cannabinoid Δ9-tetrahydrocannabinol (THC) is produced in marijuana but not in hemp. Resequencing the hemp cultivars 'Finola' and 'USO-31' showed little difference in gene copy numbers of cannabinoid pathway enzymes. However, single nucleotide variant analysis uncovered a relatively high level of variation among four cannabis types, and supported a separation of marijuana and hemp., Conclusions: The availability of the Cannabis sativa genome enables the study of a multifunctional plant that occupies a unique role in human culture. Its availability will aid the development of therapeutic marijuana strains with tailored cannabinoid profiles and provide a basis for the breeding of hemp with improved agronomic characteristics.

Atherosclerosis, the leading cause of death in the Western world, is driven by chronic inflammation within the artery wall. Elements of the complement cascade are implicated in the pathogenesis, because complement proteins and their activation products are found in the atherosclerotic plaque. We examined the role of CD55, a membrane inhibitor of the complement component 3 (C3) convertase, which converts C3 into C3a and C3b, in atherosclerosis. CD55-deficient (CD55(-/-)) mice were crossed onto the atherosclerosis-prone apolipoprotein E (apoE)-deficient (apoE(-/-)) background. High fat-fed male apoE(-/-)/CD55(-/-) mice were strongly protected from developing atherosclerosis compared with apoE(-/-) controls. Lipid profiling showed significantly lower levels of triglycerides, nonesterified fatty acids, and cholesterol in apoE(-/-)/CD55(-/-) mice than that in controls after high-fat feeding, whereas body fat in apoE(-/-)/CD55(-/-) mice content was increased. Plasma levels of C3 fell, whereas concentrations of C3adesArg (alias acylation stimulating protein; ASP), produced by serum carboxypeptidase N-mediated desargination of C3a, increased in nonfasted high fat-fed apoE(-/-)/CD55(-/-) mice, indicating complement activation. Thus, complement dysregulation in the absence of CD55 provoked increased C3adesArg production that, in turn, caused altered lipid handling, resulting in atheroprotection and increased adiposity. Interventions that target complement activation in adipose tissue should be explored as lipid-decreasing strategies., (Copyright © 2011 American Society for Investigative Pathology. Published by Elsevier Inc. All rights reserved.)

The cAMP response element-binding protein (CREB) is a highly conserved transcription factor that integrates signaling through the cAMP-dependent protein kinase A (PKA) in many eukaryotes. PKA plays a critical role in Dictyostelium development but no CREB homologue has been identified in this system. Here we show that Dictyostelium utilizes a CREB-like protein, BzpF, to integrate PKA signaling during late development. bzpF(-) mutants produce compromised spores, which are extremely unstable and germination defective. Previously, we have found that BzpF binds the canonical CRE motif in vitro. In this paper, we determined the DNA binding specificity of BzpF using protein binding microarray (PBM) and showed that the motif with the highest specificity is a CRE-like sequence. BzpF is necessary to activate the transcription of at least 15 PKA-regulated, late-developmental target genes whose promoters contain BzpF binding motifs. BzpF is sufficient to activate two of these genes. The comparison of RNA sequencing data between wild type and bzpF(-) mutant revealed that the mutant fails to express 205 genes, many of which encode cellulose-binding and sugar-binding proteins. We propose that BzpF is a CREB-like transcription factor that regulates spore maturation and stability in a PKA-related manner., (Copyright © 2011 Elsevier Inc. All rights reserved.)

Alternative splicing (AS) is a key process underlying the expansion of proteomic diversity and the regulation of gene expression. Here, we identify an evolutionarily conserved embryonic stem cell (ESC)-specific AS event that changes the DNA-binding preference of the forkhead family transcription factor FOXP1. We show that the ESC-specific isoform of FOXP1 stimulates the expression of transcription factor genes required for pluripotency, including OCT4, NANOG, NR5A2, and GDF3, while concomitantly repressing genes required for ESC differentiation. This isoform also promotes the maintenance of ESC pluripotency and contributes to efficient reprogramming of somatic cells into induced pluripotent stem cells. These results reveal a pivotal role for an AS event in the regulation of pluripotency through the control of critical ESC-specific transcriptional programs., (Copyright © 2011 Elsevier Inc. All rights reserved.)

H-NS and Lsr2 are nucleoid-associated proteins from Gram-negative bacteria and Mycobacteria, respectively, that play an important role in the silencing of horizontally acquired foreign DNA that is more AT-rich than the resident genome. Despite the fact that Lsr2 and H-NS proteins are dissimilar in sequence and structure, they serve apparently similar functions and can functionally complement one another. The mechanism by which these xenogeneic silencers selectively target AT-rich DNA has been enigmatic. We performed high-resolution protein binding microarray analysis to simultaneously assess the binding preference of H-NS and Lsr2 for all possible 8-base sequences. Concurrently, we performed a detailed structure-function relationship analysis of their C-terminal DNA binding domains by NMR. Unexpectedly, we found that H-NS and Lsr2 use a common DNA binding mechanism where a short loop containing a "Q/RGR" motif selectively interacts with the DNA minor groove, where the highest affinity is for AT-rich sequences that lack A-tracts. Mutations of the Q/RGR motif abolished DNA binding activity. Netropsin, a DNA minor groove-binding molecule effectively outcompeted H-NS and Lsr2 for binding to AT-rich sequences. These results provide a unified molecular mechanism to explain findings related to xenogeneic silencing proteins, including their lack of apparent sequence specificity but preference for AT-rich sequences. Our findings also suggest that structural information contained within the DNA minor groove is deciphered by xenogeneic silencing proteins to distinguish genetic material that is self from nonself.

C2H2 zinc fingers (C2H2-ZFs) are the most prevalent type of vertebrate DNA-binding domain, and typically appear in tandem arrays (ZFAs), with sequential C2H2-ZFs each contacting three (or more) sequential bases. C2H2-ZFs can be assembled in a modular fashion, providing one explanation for their remarkable evolutionary success. Given a set of modules with defined three-base specificities, modular assembly also presents a way to construct artificial proteins with specific DNA-binding preferences. However, a recent survey of a large number of three-finger ZFAs engineered by modular assembly reported high failure rates (∼70%), casting doubt on the generality of modular assembly. Here, we used protein-binding microarrays to analyze 28 ZFAs that failed in the aforementioned study. Most (17) preferred specific sequences, which in all but one case resembled the intended target sequence. Like natural ZFAs, the engineered ZFAs typically yielded degenerate motifs, binding dozens to hundreds of related individual sequences. Thus, the failure of these proteins in previous assays is not due to lack of sequence-specific DNA-binding activity. Our findings underscore the relevance of individual C2H2-ZF sequence specificities within tandem arrays, and support the general ability of modular assembly to produce ZFAs with sequence-specific DNA-binding activity.

Tiling microarrays have proven to be a valuable tool for gaining insights into the transcriptomes of microbial organisms grown under various nutritional or stress conditions. Here, we describe the use of such an array, constructed at the level of 20 nt resolution for the Escherichia coli MG1655 genome, to observe genome-wide changes in the steady-state RNA levels in mutants defective in either RNase E or RNase III. The array data were validated by comparison to previously published results for a variety of specific transcripts as well as independent northern analysis of additional mRNAs and sRNAs. In the absence of RNase E, 60% of the annotated coding sequences showed either increases or decreases in their steady-state levels. In contrast, only 12% of the coding sequences were affected in the absence of RNase III. Unexpectedly, many coding sequences showed decreased abundance in the RNase E mutant, while more than half of the annotated sRNAs showed changes in abundance. Furthermore, the steady-state levels of many transcripts showed overlapping effects of both ribonucleases. Data are also presented demonstrating how the arrays were used to identify potential new genes, RNase III cleavage sites and the direct or indirect control of specific biological pathways.

This chapter briefly summarizes the topics in this volume.

The RNA-Binding Protein DataBase (RBPDB) is a collection of experimental observations of RNA-binding sites, both in vitro and in vivo, manually curated from primary literature. To build RBPDB, we performed a literature search for experimental binding data for all RNA-binding proteins (RBPs) with known RNA-binding domains in four metazoan species (human, mouse, fly and worm). In total, RPBDB contains binding data on 272 RBPs, including 71 that have motifs in position weight matrix format, and 36 sets of sequences of in vivo-bound transcripts from immunoprecipitation experiments. The database is accessible by a web interface which allows browsing by domain or by organism, searching and export of records, and bulk data downloads. Users can also use RBPDB to scan sequences for RBP-binding sites. RBPDB is freely available, without registration at http://rbpdb.ccbr.utoronto.ca/.

Transcription factors (TFs) play key roles in the regulation of gene expression by binding in a sequence-specific manner to genomic DNA. In eukaryotes, DNA binding is achieved by a wide range of structural forms and motifs. TFs are typically classified by their DNA-binding domain (DBD) type. In this chapter, we catalogue and survey 91 different TF DBD types in metazoa, plants, fungi, and protists. We briefly discuss well-characterized TF families representing the major DBD superclasses. We also examine the species distributions and inferred evolutionary histories of the various families, and the potential roles played by TF family expansion and dimerization.

Classifying proteins into subgroups with similar molecular function on the basis of sequence is an important step in deriving reliable functional annotations computationally. So far, however, available classification procedures have been evaluated against protein subgroups that are defined by experts using mainly qualitative descriptions of molecular function. Recently, in vitro DNA-binding preferences to all possible 8-nt DNA sequences have been measured for 178 mouse homeodomains using protein-binding microarrays, offering the unprecedented opportunity of evaluating the classification methods against quantitative measures of molecular function. To this end, we automatically derive homeodomain subtypes from the DNA-binding data and independently group the same domains using sequence information alone. We test five sequence-based methods, which use different sequence-similarity measures and algorithms to group sequences. Results show that methods that optimize the classification robustness reflect well the detailed functional specificity revealed by the experimental data. In some of these classifications, 73-83% of the subfamilies exactly correspond to, or are completely contained in, the function-based subtypes. Our findings demonstrate that certain sequence-based classifications are capable of yielding very specific molecular function annotations. The availability of quantitative descriptions of molecular function, such as DNA-binding data, will be a key factor in exploiting this potential in the future.

Background: A gene network's capacity to process information, so as to bind past events to future actions, depends on its structure and logic. From previous and new microarray measurements in Saccharomyces cerevisiae following gene deletions and overexpressions, we identify a core gene regulatory network (GRN) of functional interactions between 328 genes and the transfer functions of each gene. Inferred connections are verified by gene enrichment., Results: We find that this core network has a generalized clustering coefficient that is much higher than chance. The inferred Boolean transfer functions have a mean p-bias of 0.41, and thus similar amounts of activation and repression interactions. However, the distribution of p-biases differs significantly from what is expected by chance that, along with the high mean connectivity, is found to cause the core GRN of S. cerevisiae's to have an overall sensitivity similar to critical Boolean networks. In agreement, we find that the amount of information propagated between nodes in finite time series is much higher in the inferred core GRN of S. cerevisiae than what is expected by chance., Conclusions: We suggest that S. cerevisiae is likely to have evolved a core GRN with enhanced information propagation among its genes.

We have investigated the role of the 3'untranslated region (3'UTR) in the expression of decay accelerating factor (DAF), one of the major membrane regulators of Complement activation. We show here that the 3'UTR of DAF contains an adenylate uridine rich element (ARE) AUUUAUUUAUAUUUAUUUA, which belongs to Class II Cluster 4 of the AU-rich element-containing mRNA (ARED) database. Enhanced Green Fluorescent Protein (EGFP) Reporter constructs containing the DAF 3'UTR showed reduced levels of expression when transfected into a variety of cell lines compared to 3'UTR reporter constructs without the ARE sequence. Furthermore, the inhibitor of mRNA transcription Actinomycin D had a much stronger effect on mRNA half-life of the ARE-containing 3'UTR demonstrating that this ARE destabilises the mRNA. Electrophoretic Mobility Shift Assays (EMSA) using biotinylated RNA probes, demonstrated that cytoplasmic Human antigen R (HuR) bound to the DAF ARE. Transfection experiments using HuR specific siRNA increased DAF expression whilst plasmids containing the coding sequence of HuR had the opposite effect, demonstrating that HuR reduced the stability of DAF mRNA and suggesting that it is of importance in regulating the expression of DAF. These data suggest that modulators of HuR could potentially be used to alter DAF expression and therefore increase the susceptibility of malignant cells to immunotherapy., (Copyright © 2010 Elsevier Ltd. All rights reserved.)

Summary: Computational gene function prediction can serve to focus experimental resources on high-priority experimental tasks. FuncBase is a web resource for viewing quantitative machine learning-based gene function annotations. Quantitative annotations of genes, including fungal and mammalian genes, with Gene Ontology terms are accompanied by a community feedback system. Evidence underlying function annotations is shown. For example, a custom Cytoscape viewer shows functional linkage graphs relevant to the gene or function of interest. FuncBase provides links to external resources, and may be accessed directly or via links from species-specific databases., Availability: FuncBase as well as all underlying data and annotations are freely available via http://func.med.harvard.edu/

Members of the large ETS family of transcription factors (TFs) have highly similar DNA-binding domains (DBDs)-yet they have diverse functions and activities in physiology and oncogenesis. Some differences in DNA-binding preferences within this family have been described, but they have not been analysed systematically, and their contributions to targeting remain largely uncharacterized. We report here the DNA-binding profiles for all human and mouse ETS factors, which we generated using two different methods: a high-throughput microwell-based TF DNA-binding specificity assay, and protein-binding microarrays (PBMs). Both approaches reveal that the ETS-binding profiles cluster into four distinct classes, and that all ETS factors linked to cancer, ERG, ETV1, ETV4 and FLI1, fall into just one of these classes. We identify amino-acid residues that are critical for the differences in specificity between all the classes, and confirm the specificities in vivo using chromatin immunoprecipitation followed by sequencing (ChIP-seq) for a member of each class. The results indicate that even relatively small differences in in vitro binding specificity of a TF contribute to site selectivity in vivo.

Metazoan genomes encode hundreds of RNA-binding proteins (RBPs). These proteins regulate post-transcriptional gene expression and have critical roles in numerous cellular processes including mRNA splicing, export, stability and translation. Despite their ubiquity and importance, the binding preferences for most RBPs are not well characterized. In vitro and in vivo studies, using affinity selection-based approaches, have successfully identified RNA sequence associated with specific RBPs; however, it is difficult to infer RBP sequence and structural preferences without specifically designed motif finding methods. In this study, we introduce a new motif-finding method, RNAcontext, designed to elucidate RBP-specific sequence and structural preferences with greater accuracy than existing approaches. We evaluated RNAcontext on recently published in vitro and in vivo RNA affinity selected data and demonstrate that RNAcontext identifies known binding preferences for several control proteins including HuR, PTB, and Vts1p and predicts new RNA structure preferences for SF2/ASF, RBM4, FUSIP1 and SLM2. The predicted preferences for SF2/ASF are consistent with its recently reported in vivo binding sites. RNAcontext is an accurate and efficient motif finding method ideally suited for using large-scale RNA-binding affinity datasets to determine the relative binding preferences of RBPs for a wide range of RNA sequences and structures.

The genetic code-the binding specificity of all transfer-RNAs--defines how protein primary structure is determined by DNA sequence. DNA also dictates when and where proteins are expressed, and this information is encoded in a pattern of specific sequence motifs that are recognized by transcription factors. However, the DNA-binding specificity is only known for a small fraction of the approximately 1400 human transcription factors (TFs). We describe here a high-throughput method for analyzing transcription factor binding specificity that is based on systematic evolution of ligands by exponential enrichment (SELEX) and massively parallel sequencing. The method is optimized for analysis of large numbers of TFs in parallel through the use of affinity-tagged proteins, barcoded selection oligonucleotides, and multiplexed sequencing. Data are analyzed by a new bioinformatic platform that uses the hundreds of thousands of sequencing reads obtained to control the quality of the experiments and to generate binding motifs for the TFs. The described technology allows higher throughput and identification of much longer binding profiles than current microarray-based methods. In addition, as our method is based on proteins expressed in mammalian cells, it can also be used to characterize DNA-binding preferences of full-length proteins or proteins requiring post-translational modifications. We validate the method by determining binding specificities of 14 different classes of TFs and by confirming the specificities for NFATC1 and RFX3 using ChIP-seq. Our results reveal unexpected dimeric modes of binding for several factors that were thought to preferentially bind DNA as monomers.

A series of reports over the last few years have indicated that a much larger portion of the mammalian genome is transcribed than can be accounted for by currently annotated genes, but the quantity and nature of these additional transcripts remains unclear. Here, we have used data from single- and paired-end RNA-Seq and tiling arrays to assess the quantity and composition of transcripts in PolyA+ RNA from human and mouse tissues. Relative to tiling arrays, RNA-Seq identifies many fewer transcribed regions ("seqfrags") outside known exons and ncRNAs. Most nonexonic seqfrags are in introns, raising the possibility that they are fragments of pre-mRNAs. The chromosomal locations of the majority of intergenic seqfrags in RNA-Seq data are near known genes, consistent with alternative cleavage and polyadenylation site usage, promoter- and terminator-associated transcripts, or new alternative exons; indeed, reads that bridge splice sites identified 4,544 new exons, affecting 3,554 genes. Most of the remaining seqfrags correspond to either single reads that display characteristics of random sampling from a low-level background or several thousand small transcripts (median length = 111 bp) present at higher levels, which also tend to display sequence conservation and originate from regions with open chromatin. We conclude that, while there are bona fide new intergenic transcripts, their number and abundance is generally low in comparison to known exons, and the genome is not as pervasively transcribed as previously reported., Competing Interests: The authors have declared that no competing interests exist.

Sequence-specific binding by transcription factors (TFs) interprets regulatory information encoded in the genome. Using recently published universal protein binding microarray (PBM) data on the in vitro DNA binding preferences of these proteins for all possible 8-base-pair sequences, we examined the evolutionary conservation and enrichment within putative regulatory regions of the binding sequences of a diverse library of 104 nonredundant mouse TFs spanning 22 different DNA-binding domain structural classes. We found that not only high affinity binding sites, but also numerous moderate and low affinity binding sites, are under negative selection in the mouse genome. These 8-mers occur preferentially in putative regulatory regions of the mouse genome, including CpG islands and non-exonic ultraconserved elements (UCEs). Of TFs whose PBM "bound" 8-mers are enriched within sets of tissue-specific UCEs, many are expressed in the same tissue(s) as the UCE-driven gene expression. Phylogenetically conserved motif occurrences of various TFs were also enriched in the noncoding sequence surrounding numerous gene sets corresponding to Gene Ontology categories and tissue-specific gene expression clusters, suggesting involvement in transcriptional regulation of those genes. Altogether, our results indicate that many of the sequences bound by these proteins in vitro, including lower affinity DNA sequences, are likely to be functionally important in vivo. This study not only provides an initial analysis of the potential regulatory associations of 104 mouse TFs, but also presents an approach for the functional analysis of TFs from any other metazoan genome as their DNA binding preferences are determined by PBMs or other technologies., (Copyright 2010 Elsevier Inc. All rights reserved.)

Activation of complement occurs during autoimmune retinal and intraocular inflammatory disease as well as neuroretinal degenerative disorders. The cleavage of C5 into fragments C5a and C5b is a critical event during the complement cascade. C5a is a potent proinflammatory anaphylatoxin capable of inducing cell migration, adhesion and cytokine release, while membrane attack complex C5b-9 causes cell lysis. Therapeutic approaches to prevent complement-induced inflammation include the use of blocking monoclonal antibodies (mAb) to prevent C5 cleavage. In these current experiments, the rat anti-mouse C5 mAb (BB5.1) was utilized to investigate the effects of inhibition of C5 cleavage on disease progression and severity in experimental autoimmune uveoretinitis (EAU), a model of organ-specific autoimmunity in the eye characterized by structural retinal damage mediated by infiltrating macrophages. Systemic treatment with BB5.1 results in significantly reduced disease scores compared with control groups, while local administration results in an earlier resolution of disease. In vitro, contemporaneous C5a and interferon-gamma signalling enhanced nitric oxide production, accompanied by down-regulation of the inhibitory myeloid CD200 receptor, contributing to cell activation. These experiments demonstrate that C5 cleavage contributes to the full expression of EAU, and that selective C5 blockade via systemic and local routes of administration can suppress disease. This presents great therapeutic potential to protect against tissue damage during autoimmune responses in the retina or inflammation-induced degenerative disease.

Active eukaryotic regulatory sites are characterized by open chromatin, and yeast promoters and transcription factor binding sites (TFBSs) typically have low intrinsic nucleosome occupancy. Here, we show that in contrast to yeast, DNA at human promoters, enhancers, and TFBSs generally encodes high intrinsic nucleosome occupancy. In most cases we examined, these elements also have high experimentally measured nucleosome occupancy in vivo. These regions typically have high G+C content, which correlates positively with intrinsic nucleosome occupancy, and are depleted for nucleosome-excluding poly-A sequences. We propose that high nucleosome preference is directly encoded at regulatory sequences in the human genome to restrict access to regulatory information that will ultimately be utilized in only a subset of differentiated cells.

Aims: To examine the roles of the membrane attack complex of complement and its sole membrane regulator, CD59, in atherosclerosis., Methods: C6 (C6(-/-)) deficient and CD59a (Cd59a(-/-)) knockout mice were separately crossed onto the apolipoprotein E knockout (apoE(-/-)) background. The double knockout mice were fed high-fat diet in order to study the effects of absence of C6 or CD59a on the progression of atherosclerosis., Results: C6 deficiency significantly reduced plaque area and disease severity. CD59a had the opposite effect in that deficiency was associated with a significant increase in plaque area, correlating with increased membrane attack complex (MAC) deposition in the plaque and increased smooth muscle cell proliferation in early plaques., Conclusions: Our results demonstrate that the MAC contributes to the development of atherosclerosis, C6 deficiency being protective and CD59a deficiency exacerbating disease., ((c) 2009 Elsevier Ltd. All rights reserved.)

Regulatory regions with similar transcriptional output often have little overt sequence similarity, both within and between genomes. Although cis- and trans-regulatory changes can contribute to sequence divergence without dramatically altering gene expression outputs, heterologous DNA often functions similarly in organisms that share little regulatory sequence similarities (e.g. human DNA in fish), indicating that trans-regulatory mechanisms tend to diverge more slowly and can accommodate a variety of cis-regulatory configurations. This capacity to 'tinker' with regulatory DNA probably relates to the complexity, robustness and evolvability of regulatory systems, but cause-and-effect relationships among evolutionary processes and properties of regulatory systems remain a topic of debate. The challenge of understanding the concrete mechanisms underlying cis-regulatory evolution - including the conservation of function without the conservation of sequence - relates to the challenge of understanding the function of regulatory systems in general. Currently, we are largely unable to recognize functionally similar regulatory DNA.

A recent study reveals a surprisingly high degree of change in the occupancy patterns of two transcription factors in the livers of five vertebrates.

We propose definitions and procedures for comparing nucleosome maps and discuss current agreement and disagreement on the effect of histone sequence preferences on nucleosome organization in vivo., (© 2010 BioMed Central Ltd)

Background: The relative preference of nucleosomes to form on individual DNA sequences plays a major role in genome packaging. A wide variety of DNA sequence features are believed to influence nucleosome formation, including periodic dinucleotide signals, poly-A stretches and other short motifs, and sequence properties that influence DNA structure, including base content. It was recently shown by Kaplan et al. that a probabilistic model using composition of all 5-mers within a nucleosome-sized tiling window accurately predicts intrinsic nucleosome occupancy across an entire genome in vitro. However, the model is complicated, and it is not clear which specific DNA sequence properties are most important for intrinsic nucleosome-forming preferences., Results: We find that a simple linear combination of only 14 simple DNA sequence attributes (G+C content, two transformations of dinucleotide composition, and the frequency of eleven 4-bp sequences) explains nucleosome occupancy in vitro and in vivo in a manner comparable to the Kaplan model. G+C content and frequency of AAAA are the most important features. G+C content is dominant, alone explaining approximately 50% of the variation in nucleosome occupancy in vitro., Conclusions: Our findings provide a dramatically simplified means to predict and understand intrinsic nucleosome occupancy. G+C content may dominate because it both reduces frequency of poly-A-like stretches and correlates with many other DNA structural characteristics. Since G+C content is enriched or depleted at many types of features in diverse eukaryotic genomes, our results suggest that variation in nucleotide composition may have a widespread and direct influence on chromatin structure.

Cwc21 (complexed with Cef1 protein 21) is a 135 amino acid yeast protein that shares homology with the N-terminal domain of human SRm300/SRRM2, a large serine/arginine-repeat protein shown previously to associate with the splicing coactivator and 3'-end processing stimulatory factor, SRm160. Proteomic analysis of spliceosomal complexes has suggested a role for Cwc21 and SRm300 at the core of the spliceosome. However, specific functions for these proteins have remained elusive. In this report, we employ quantitative genetic interaction mapping, mass spectrometry of tandem affinity-purified complexes, and microarray profiling to investigate genetic, physical, and functional interactions involving Cwc21. Combined data from these assays support multiple roles for Cwc21 in the formation and function of splicing complexes. Consistent with a role for Cwc21 at the core of the spliceosome, we observe strong genetic, physical, and functional interactions with Isy1, a protein previously implicated in the first catalytic step of splicing and splicing fidelity. Together, the results suggest multiple functions for Cwc21/SRm300 in the splicing process, including an important role in the activation of splicing in association with Isy1.

The last decade has seen an explosion of interest in new classes of non-coding RNA. While some are now firmly established as new categories of legitimate functional RNAs, the purpose and even existence of others remain to be solidified. Here, we discuss the challenges associated with discovery and characterization of non-traditional categories of non-coding RNA.

We describe a fluorescent reporter system that exploits the functional genomic tools available in budding yeast to systematically assess consequences of genetic perturbations on gene expression. We used our Reporter-Synthetic Genetic Array (R-SGA) method to screen for regulators of core histone gene expression. We discovered that the histone chaperone Rtt106 functions in a pathway with two other chaperones, Asf1 and the HIR complex, to create a repressive chromatin structure at core histone promoters. We found that activation of histone (HTA1) gene expression involves both relief of Rtt106-mediated repression by the activity of the histone acetyltransferase Rtt109 and restriction of Rtt106 to the promoter region by the bromodomain-containing protein Yta7. We propose that the maintenance of Asf1/HIR/Rtt106-mediated repressive chromatin domains is the primary mechanism of cell-cycle regulation of histone promoters. Our data suggest that this pathway may represent a chromatin regulatory mechanism that is broadly used across the genome.

Metazoan genomes encode hundreds of RNA-binding proteins (RBPs) but RNA-binding preferences for relatively few RBPs have been well defined. Current techniques for determining RNA targets, including in vitro selection and RNA co-immunoprecipitation, require significant time and labor investment. Here we introduce RNAcompete, a method for the systematic analysis of RNA binding specificities that uses a single binding reaction to determine the relative preferences of RBPs for short RNAs that contain a complete range of k-mers in structured and unstructured RNA contexts. We tested RNAcompete by analyzing nine diverse RBPs (HuR, Vts1, FUSIP1, PTB, U1A, SF2/ASF, SLM2, RBM4 and YB1). RNAcompete identified expected and previously unknown RNA binding preferences. Using in vitro and in vivo binding data, we demonstrate that preferences for individual 7-mers identified by RNAcompete are a more accurate representation of binding activity than are conventional motif models. We anticipate that RNAcompete will be a valuable tool for the study of RNA-protein interactions.

The development and application of genomic reagents and techniques has fuelled progress in our understanding of regulatory networks that control gene expression in eukaryotic cells. However, a full description of the network of regulator-gene interactions that determine global gene expression programs remains elusive and will require systematic genetic as well as biochemical assays. Here, we describe a functional genomics approach that combines reporter technology, genome-wide array-based reagents and high-throughput imaging to discover new regulators controlling gene expression patterns in Saccharomyces cerevisiae. Our strategy utilizes the synthetic genetic array (SGA) method to systematically introduce promoter-GFP (green fluorescent protein) reporter constructs along with a control promoter-RFP (red fluorescent protein) gene into the array of approximately 4500 viable yeast deletion mutants. Fluorescence intensities from each reporter are assayed from individual colonies arrayed on solid agar plates using a scanning fluorimager and the ratio of GFP to RFP intensity reveals deletion mutants that cause differential GFP expression. We are exploiting this screening approach to construct a detailed map describing the interplay of regulators controlling the eukaryotic cell cycle. The method is extensible to any transcription factor or signalling pathway for which an appropriate reporter gene can be devised.

Sequence preferences of DNA binding proteins are a primary mechanism by which cells interpret the genome. Despite the central importance of these proteins in physiology, development, and evolution, comprehensive DNA binding specificities have been determined experimentally for only a few proteins. Here, we used microarrays containing all 10-base pair sequences to examine the binding specificities of 104 distinct mouse DNA binding proteins representing 22 structural classes. Our results reveal a complex landscape of binding, with virtually every protein analyzed possessing unique preferences. Roughly half of the proteins each recognized multiple distinctly different sequence motifs, challenging our molecular understanding of how proteins interact with their DNA binding sites. This complexity in DNA recognition may be important in gene regulation and in the evolution of transcriptional regulatory networks.

Motivation: Recognition of specific DNA sequences is a central mechanism by which transcription factors (TFs) control gene expression. Many TF-binding preferences, however, are unknown or poorly characterized, in part due to the difficulty associated with determining their specificity experimentally, and an incomplete understanding of the mechanisms governing sequence specificity. New techniques that estimate the affinity of TFs to all possible k-mers provide a new opportunity to study DNA-protein interaction mechanisms, and may facilitate inference of binding preferences for members of a given TF family when such information is available for other family members., Results: We employed a new dataset consisting of the relative preferences of mouse homeodomains for all eight-base DNA sequences in order to ask how well we can predict the binding profiles of homeodomains when only their protein sequences are given. We evaluated a panel of standard statistical inference techniques, as well as variations of the protein features considered. Nearest neighbour among functionally important residues emerged among the most effective methods. Our results underscore the complexity of TF-DNA recognition, and suggest a rational approach for future analyses of TF families.

Background: Transfer RNA (tRNA) gene predictions are complicated by challenges such as structural variation, limited sequence conservation and the presence of highly reiterated short interspersed sequences (SINEs) that originally derived from tRNA genes or tRNA-like transcription units. Annotation of "tRNA genes" in sequenced genomes generally have not been accompanied by experimental verification of the expression status of predicted sequences., Results: To address this for mouse tRNA genes, we have employed two programs, tRNAScan-SE and ARAGORN, to predict the tRNA genes in the nuclear genome, resulting in diverse but overlapping predicted gene sets. From these, we removed known SINE repeats and sorted the genes into predicted families and single-copy genes. In particular, four families of intron-containing tRNA genes were predicted for the first time in mouse, with introns in positions and structures similar to the well characterized intron-containing tRNA genes in yeast. We verified the expression of the predicted tRNA genes by microarray analysis. We then confirmed the expression of appropriately sized RNA for the four intron-containing tRNA gene families, as well as the other 30 tRNA gene families creating an index of expression-verified mouse tRNAs., Conclusions: These confirmed tRNA genes represent all anticodons and all known mammalian tRNA structural groups, as well as a variety of predicted "rogue" tRNA genes within families with altered anticodon identities.

The unique species specificity of the bacterial cytolysin intermedilysin is explained by its requirement for the human complement regulator CD59 as the primary receptor. Binding studies using individual domains of intermedilysin mapped the CD59-binding site to domain 4 and swap mutants between human and rabbit (non-intermedilysin-binding) CD59 implicated a short sequence (residues 42-59) in human CD59 in binding intermedilysin. We set out to map more closely the CD59 binding site in intermedilysin. We first looked for regions of homology between domain 4 in intermedilysin and the terminal complement components that bind CD59, C8 and C9. A nine amino acid sequence immediately adjacent the undecapeptide segment in intermedilysin domain 4 matched (5 of 9 identical, 3 of 9 conserved) a sequence in C9. A peptide containing this sequence caused dose-dependent inhibition of intermedilysin-mediated lysis of human erythrocytes and rendered erythrocytes more susceptible to complement lysis. Surface plasmon resonance analysis of intermedilysin binding to immobilized CD59 revealed saturable fast-on, fast-off binding and a calculated affinity of 4.9 nM. Substitution of three residues from the putative binding site caused a 5-fold reduction in lytic potency of intermedilysin and reduced affinity for immobilized CD59 by 2.5-fold. The demonstration that a peptide modeled on the CD59-binding site inhibits intermedilysin-mediated haemolysis leads us to suggest that such peptides might be useful in treating infections caused by intermedilysin-producing bacteria.