Inferring Contagion in Regulatory Networks

Several gene regulatory network models containing concepts of directionality at the edges have been proposed. However, only a few reports have an interpretable definition of directionality. Here, differently from the standard causality concept defined by Pearl, we introduce the concept of contagion in order to infer directionality at the edges, i.e., asymmetries in gene expression dependences of regulatory networks. Moreover, we present a bootstrap algorithm in order to test the contagion concept. This technique was applied in simulated data and, also, in an actual large sample of biological data. Literature review has confirmed some genes identified by contagion as actually belonging to the TP53 pathway.

GRASSIUS: A Platform for Comparative Regulatory Genomics across the Grasses

Transcription factors (TFs) are major players in gene regulatory networks and interactions between TFs and their target genes furnish spatiotemporal patterns of gene expression. Establishing the architecture of regulatory networks requires gathering information on TFs, their targets in the genome, and the corresponding binding sites. We have developed GRASSIUS (Grass Regulatory Information Services) as a knowledge-based Web resource that integrates information on TFs and gene promoters across the grasses. In its initial implementation, GRASSIUS consists of two separate, yet linked, databases. GrassTFDB holds information on TFs from maize (Zea mays), sorghum (Sorghum bicolor), sugarcane (Saccharum spp.), and rice (Oryza sativa). TFs are classified into families and phylogenetic relationships begin to uncover orthologous relationships among the participating species. This database also provides a centralized clearinghouse for TF synonyms in the grasses. GrassTFDB is linked to the grass TFome collection, which provides clones in recombination-based vectors corresponding to full-length open reading frames for a growing number of grass TFs. GrassPROMDB contains promoter and cis-regulatory element information for those grass species and genes for which enough data are available. The integration of GrassTFDB and GrassPROMDB will be accomplished through GrassRegNet as a first step in representing the architecture of grass regulatory networks. GRASSIUS can be accessed from www.grassius.org.

Fur controls iron homeostasis and oxidative stress defense in the oligotrophic alpha-proteobacterium Caulobacter crescentus

In most bacteria, the ferric uptake regulator (Fur) is a global regulator that controls iron homeostasis and other cellular processes, such as oxidative stress defense. In this work, we apply a combination of bioinformatics, in vitro and in vivo assays to identify the Caulobacter crescentus Fur regulon. A C. crescentus fur deletion mutant showed a slow growth phenotype, and was hypersensitive to H(2)O(2) and organic peroxide. Using a position weight matrix approach, several predicted Fur-binding sites were detected in the genome of C. crescentus, located in regulatory regions of genes not only involved in iron uptake and usage but also in other functions. Selected Fur-binding sites were validated using electrophoretic mobility shift assay and DNAse I footprinting analysis. Gene expression assays revealed that genes involved in iron uptake were repressed by iron-Fur and induced under conditions of iron limitation, whereas genes encoding iron-using proteins were activated by Fur under conditions of iron sufficiency. Furthermore, several genes that are regulated via small RNAs in other bacteria were found to be directly regulated by Fur in C. crescentus. In conclusion, Fur functions as an activator and as a repressor, integrating iron metabolism and oxidative stress response in C. crescentus.

Functionally conserved cis-regulatory elements of COL18A1 identified through zebrafish transgenesis

Type XVIII collagen is a component of basement membranes, and expressed prominently in the eye, blood vessels, liver, and the central nervous system. Homozygous mutations in COL18A1 lead to Knobloch Syndrome, characterized by ocular defects and occipital encephalocele. However, relatively little has been described on the role of type XVIII collagen in development, and nothing is known about the regulation of its tissue-specific expression pattern. We have used zebrafish transgenesis to identify and characterize cis-regulatory sequences controlling expression of the human gene. Candidate enhancers were selected from non-coding sequence associated with COL18A1 based on sequence conservation among mammals. Although these displayed no overt conservation with orthologous zebrafish sequences, four regions nonetheless acted as tissue-specific transcriptional enhancers in the zebrafish embryo, and together recapitulated the major aspects of col18a1 expression. Additional post-hoc computational analysis on positive enhancer sequences revealed alignments between mammalian and teleost sequences, which we hypothesize predict the corresponding zebrafish enhancers; for one of these, we demonstrate functional overlap with the orthologous human enhancer sequence. Our results provide important insight into the biological function and regulation of COL18A1, and point to additional sequences that may contribute to complex diseases involving COL18A1. More generally, we show that combining functional data with targeted analyses for phylogenetic conservation can reveal conserved cis-regulatory elements in the large number of cases where computational alignment alone falls short. (C) 2009 Elsevier Inc. All rights reserved.

Transcriptional regulation differs in affected facioscapulohumeral muscular dystrophy patients compared to asymptomatic related carriers

Facioscapulohumeral muscular dystrophy (FSHD) is a progressive muscle disorder that has been associated with a contraction of 3.3-kb repeats on chromosome 4q35. FSHD is characterized by a wide clinical inter- and intrafamilial variability, ranging from wheelchair-bound patients to asymptomatic carriers. Our study is unique in comparing the gene expression profiles from related affected, asymptomatic carrier, and control individuals. Our results suggest that the expression of genes on chromosome 4q is altered in affected and asymptomatic individuals. Remarkably, the changes seen in asymptomatic samples are largely in products of genes encoding several chemokines, whereas the changes seen in affected samples are largely in genes governing the synthesis of GPI-linked proteins and histone acetylation. Besides this, the affected patient and related asymptomatic carrier share the 4qA161 haplotype. Thus, these polymorphisms by themselves do not explain the pathogenicity of the contracted allele. Interestingly, our results also suggest that the miRNAs might mediate the regulatory network in FSHD. Together, our results support the previous evidence that FSHD may be caused by transcriptional dysregulation of multiple genes, in cis and in trans, and suggest some factors potentially important for FSHD pathogenesis. The study of the gene expression profiles from asymptomatic carriers and related affected patients is a unique approach to try to enhance our understanding of the missing link between the contraction in D4Z4 repeats and muscle disease, while minimizing the effects of differences resulting from genetic background.

CztR, a LysR-Type Transcriptional Regulator Involved in Zinc Homeostasis and Oxidative Stress Defense in Caulobacter crescentus

Caulobacter crescentus is a free-living alphaproteobacterium that has 11 predicted LysR-type transcriptional regulators (LTTRs). Previously, a C. crescentus mutant strain with a mini-Tn5lacZ transposon inserted into a gene encoding an LTTR was isolated; this mutant was sensitive to cadmium. In this work, a mutant strain with a deletion was obtained, and the role of this LTTR (called CztR here) was evaluated. The transcriptional start site of this gene was determined by primer extension analysis, and its promoter was cloned in front of a lacZ reporter gene. beta-Galactosidase activity assays, performed with the wild-type and mutant strains, indicated that this gene is 2-fold induced when cells enter stationary phase and that it is negatively autoregulated. Moreover, this regulator is essential for the expression of the divergent cztA gene at stationary phase, in minimal medium, and in response to zinc depletion. This gene encodes a hypothetical protein containing 10 predicted transmembrane segments, and its expression pattern suggests that it encodes a putative zinc transporter. The cztR strain was also shown to be sensitive to superoxide (generated by paraquat) and to hydrogen peroxide but not to tert-butyl hydroperoxide. The expression of katG and ahpC, but not that of the superoxide dismutase genes, was increased in the cztR mutant. A model is proposed to explain how CztR binding to the divergent regulatory regions could activate cztA expression and repress its own transcription.

Transcriptional Processing of the pst Operon of Escherichia coli

The pst operon of Escherichia coli is composed of five genes that encode a high-affinity phosphate transport system. pst belongs to the PHO regulon, which is a group of genes and operons that are induced in response to phosphate limitation. The pst operon also has a regulatory role in the repression of PHO genes` transcription under phosphate excess conditions. Transcription of pst is initiated at the promoter located upstream to the first gene, pstS. Immediately after its synthesis, the primary transcript of pst is cleaved into shorter mRNA molecules in a ribonuclease E-dependent manner. Other ribonucleases, such as RNase III and MazF, do not play a role in pst mRNA processing. RNase E is thus at least partially responsible for processing the pst primary transcript.

Intrinsically Multivariate Predictive Genes

Canalizing genes possess such broad regulatory power, and their action sweeps across a such a wide swath of processes that the full set of affected genes are not highly correlated under normal conditions. When not active, the controlling gene will not be predictable to any significant degree by its subject genes, either alone or in groups, since their behavior will be highly varied relative to the inactive controlling gene. When the controlling gene is active, its behavior is not well predicted by any one of its targets, but can be very well predicted by groups of genes under its control. To investigate this question, we introduce in this paper the concept of intrinsically multivariate predictive (IMP) genes, and present a mathematical study of IMP in the context of binary genes with respect to the coefficient of determination (CoD), which measures the predictive power of a set of genes with respect to a target gene. A set of predictor genes is said to be IMP for a target gene if all properly contained subsets of the predictor set are bad predictors of the target but the full predictor set predicts the target with great accuracy. We show that logic of prediction, predictive power, covariance between predictors, and the entropy of the joint probability distribution of the predictors jointly affect the appearance of IMP genes. In particular, we show that high-predictive power, small covariance among predictors, a large entropy of the joint probability distribution of predictors, and certain logics, such as XOR in the 2-predictor case, are factors that favor the appearance of IMP. The IMP concept is applied to characterize the behavior of the gene DUSP1, which exhibits control over a central, process-integrating signaling pathway, thereby providing preliminary evidence that IMP can be used as a criterion for discovery of canalizing genes.

Human regulatory protein Ki-1/57 has characteristics of an intrinsically unstructured protein

The human protein Ki-1/57 was first identified through the cross reactivity of the anti-CD30 monoclonal antibody Ki-1; in Hodgkin lymphoma cells. The expression of Ki-1/57 in diverse cancer cells and its phosphorylation in peripheral blood leukocytes after mitogenic activation suggested its possible role in cell signaling. Ki-1/57 interacts with several other regulatory proteins involved in cellular signaling, transcriptional regulation and RNA metabolism, suggesting it may have pleiotropic functions. In a previous spectroscopic analysis, we observed a low content of secondary structure for Ki-1/57 constructs. Here, Circular dichroism experiments, in vitro RNA binding analysis, and limited proteolysis assays of recombinant Ki-1/57(122-413) and proteolysis assays of endogenous full length protein from human HEK293 cells suggested that Ki-1/57 has characteristics of an intrinsically unstructured protein. Small-angle X-ray scattering (SAXS) experiments were performed with the C-terminal fragment Ki-1/57(122-413). These results indicated an elongated shape and a partially unstructured conformation of the molecule in solution, confirming the characteristics of an intrinsically unstructured protein. Experimental curves together with ab initio modeling approaches revealed an extended and flexible molecule in solution. An elongated shape was also observed by analytical gel filtration. Furthermore, sedimentation velocity analysis suggested that Ki-1/57 is a highly asymmetric protein. These findings may explain the functional plasticity of Ki-1/57, as suggested by the wide array of proteins with which it is capable of interacting in yeast two-hybrid interaction assays.