Genetic Susceptibility Loci in Rheumatoid Arthritis Establish Transcriptional Regulatory Networks with Other Genes

Linkage studies have identified the human leukocyte antigen (HLA)-DRB1 as a putative rheumatoid arthritis (RA) susceptibility locus (SL). Nevertheless, it was estimated that its contribution was partial, suggesting that other non-HLA genes may play a role in RA susceptibility. To test this hypothesis, we conducted microarray transcription profiling of peripheral blood mononuclear cells in 15 RA patients and analyzed the data, using bioinformatics programs (significance analysis of microarrays method and GeneNetwork), which allowed us to determine the differentially expressed genes and to reconstruct transcriptional networks. The patients were grouped according to disease features or treatment with tumor necrosis factor blocker. Transcriptional networks that were reconstructed allowed us to identify the interactions occurring between RA SL and other genes, for example, HLA-DRB1 interacting with FNDC3A (fibronectin type III domain containing 3A). Given that fibronectin fragments can stimulate mediators of matrix and cartilage destruction in RA, this interaction is of special interest and may contribute to a clearer understanding of the functional role of HLA-DRB1 in RA pathogenesis.

Genes that code for T cell signaling proteins establish transcriptional regulatory networks during thymus ontogeny

Gene expression profiling by cDNA microarrays during murine thymus ontogeny has contributed to dissecting the large-scale molecular genetics of T cell maturation. Gene profiling, although useful for characterizing the thymus developmental phases and identifying the differentially expressed genes, does not permit the determination of possible interactions between genes. In order to reconstruct genetic interactions, on RNA level, within thymocyte differentiation, a pair of microarrays containing a total of 1,576 cDNA sequences derived from the IMAGE MTB library was applied on samples of developing thymuses (14-17 days of gestation). The data were analyzed using the GeneNetwork program. Genes that were previously identified as differentially expressed during thymus ontogeny showed their relationships with several other genes. The present method provided the detection of gene nodes coding for proteins implicated in the calcium signaling pathway, such as Prrg2 and Stxbp3, and in protein transport toward the cell membrane, such as Gosr2. The results demonstrate the feasibility of reconstructing networks based on cDNA microarray gene expression determinations, contributing to a clearer understanding of the complex interactions between genes involved in thymus/thymocyte development.

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.

The impact of measurement errors in the identification of regulatory networks

Background: There are several studies in the literature depicting measurement error in gene expression data and also, several others about regulatory network models. However, only a little fraction describes a combination of measurement error in mathematical regulatory networks and shows how to identify these networks under different rates of noise. Results: This article investigates the effects of measurement error on the estimation of the parameters in regulatory networks. Simulation studies indicate that, in both time series (dependent) and non-time series (independent) data, the measurement error strongly affects the estimated parameters of the regulatory network models, biasing them as predicted by the theory. Moreover, when testing the parameters of the regulatory network models, p-values computed by ignoring the measurement error are not reliable, since the rate of false positives are not controlled under the null hypothesis. In order to overcome these problems, we present an improved version of the Ordinary Least Square estimator in independent (regression models) and dependent (autoregressive models) data when the variables are subject to noises. Moreover, measurement error estimation procedures for microarrays are also described. Simulation results also show that both corrected methods perform better than the standard ones (i.e., ignoring measurement error). The proposed methodologies are illustrated using microarray data from lung cancer patients and mouse liver time series data. Conclusions: Measurement error dangerously affects the identification of regulatory network models, thus, they must be reduced or taken into account in order to avoid erroneous conclusions. This could be one of the reasons for high biological false positive rates identified in actual regulatory network models.

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.

Complex networks: the key to systems biology

Though introduced recently, complex networks research has grown steadily because of its potential to represent, characterize and model a wide range of intricate natural systems and phenomena. Because of the intrinsic complexity and systemic organization of life, complex networks provide a specially promising framework for systems biology investigation. The current article is an up-to-date review of the major developments related to the application of complex networks in biology, with special attention focused on the more recent literature. The main concepts and models of complex networks are presented and illustrated in an accessible fashion. Three main types of networks are covered: transcriptional regulatory networks, protein-protein interaction networks and metabolic networks. The key role of complex networks for systems biology is extensively illustrated by several of the papers reviewed.

Evidence for a network transcriptional control of promiscuous gene expression in medullary thymic epithelial cells

The expression of peripheral tissue antigens (PTAs) in the thymus by medullary thymic epithelial cells (mTECs) is essential for the central self-tolerance in the generation of the T cell repertoire. Due to heterogeneity of autoantigen representation, this phenomenon has been termed promiscuous gene expression (PGE), in which the autoimmune regulator (Aire) gene plays a key role as a transcription factor in part of these genes. Here we used a microarray strategy to access PGE in cultured murine CD80(+) 3.10 mTEC line. Hierarchical clustering of the data allowed observation that PTA genes were differentially expressed being possible to found their respective induced or repressed mRNAs. To further investigate the control of PGE, we tested the hypothesis that genes involved in this phenomenon might also be modulated by transcriptional network. We then reconstructed such network based on the microarray expression data, featuring the guanylate cyclase 2d (Gucy2d) gene as a main node. In such condition, we established 167 positive and negative interactions with downstream PTA genes. Silencing Aire by RNA interference, Gucy2d while down regulated established a larger number (355) of interactions with PTA genes. T- and G-boxes corresponding to AIRE protein binding sites located upstream to ATG codon of Gucy2d supports this effect. These findings provide evidence that Aire plays a role in association with Gucy2d, which is connected to Several PTA genes and establishes a cascade-like transcriptional control of promiscuous gene expression in mTEC cells. (C) 2009 Elsevier Ltd. All rights reserved.

Transcriptional changes in the nuc-2A mutant strain of Neurospora crassa cultivated under conditions of phosphate shortage

The molecular mechanism that controls the response to phosphate shortage in Neurospora crassa involves four regulatory genes - nuc-2, preg, pgov, and nuc-1. Phosphate shortage is sensed by the nuc-2 gene, the product of which inhibits the functioning of the PREG-PGOV complex. This allows the translocation of the transcriptional factor NUC-1 into the nucleus, which activates the transcription of phosphate-repressible phosphatases. The nuc-2A mutant strain of N. crassa carries a loss-of-function mutation in the nuc-2 gene, which encodes an ankyrin-like repeat protein. In this study, we identified transcripts that are downregutated in the nuc-2A mutant strain. Functional grouping of these expressed sequence tags allowed the identification of genes that play essential roles in different cellular processes such as transport, transcriptional regulation, signal transduction, metabolism, protein synthesis, protein fate, and development. These results reveal novel aspects of the phosphorus-sensing network in N. crassa. (C) 2009 Elsevier GmbH. All rights reserved.

Activation of Neural and Pluripotent Stem Cell Signatures Correlates with Increased Malignancy in Human Glioma

The presence of stem cell characteristics in glioma cells raises the possibility that mechanisms promoting the maintenance and self-renewal of tissue specific stem cells have a similar function in tumor cells. Here we characterized human gliomas of various malignancy grades for the expression of stem cell regulatory proteins. We show that cells in high grade glioma co-express an array of markers defining neural stem cells (NSCs) and that these proteins can fulfill similar functions in tumor cells as in NSCs. However, in contrast to NSCs glioma cells co-express neural proteins together with pluripotent stem cell markers, including the transcription factors Oct4, Sox2, Nanog and Klf4. In line with this finding, in high grade gliomas mesodermal-and endodermal-specific transcription factors were detected together with neural proteins, a combination of lineage markers not normally present in the central nervous system. Persistent presence of pluripotent stem cell traits could only be detected in solid tumors, and observations based on in vitro studies and xenograft transplantations in mice imply that this presence is dependent on the combined activity of intrinsic and extrinsic regulatory cues. Together these results demonstrate a general deregulated expression of neural and pluripotent stem cell traits in malignant human gliomas, and indicate that stem cell regulatory factors may provide significant targets for therapeutic strategies.

Complete Genome Viral Phylogenies Suggests the Concerted Evolution of Regulatory Cores and Accessory Satellites

We consider the concerted evolution of viral genomes in four families of DNA viruses. Given the high rate of horizontal gene transfer among viruses and their hosts, it is an open question as to how representative particular genes are of the evolutionary history of the complete genome. To address the concerted evolution of viral genes, we compared genomic evolution across four distinct, extant viral families. For all four viral families we constructed DNA-dependent DNA polymerase-based (DdDp) phylogenies and in addition, whole genome sequence, as quantitative descriptions of inter-genome relationships. We found that the history of the polymerase gene was highly predictive of the history of the genome as a whole, which we explain in terms of repeated, co-divergence events of the core DdDp gene accompanied by a number of satellite, accessory genetic loci. We also found that the rate of gene gain in baculovirus and poxviruses proceeds significantly more quickly than the rate of gene loss and that there is convergent acquisition of satellite functions promoting contextual adaptation when distinct viral families infect related hosts. The congruence of the genome and polymerase trees suggests that a large set of viral genes, including polymerase, derive from a phylogenetically conserved core of genes of host origin, secondarily reinforced by gene acquisition from common hosts or co-infecting viruses within the host. A single viral genome can be thought of as a mutualistic network, with the core genes acting as an effective host and the satellite genes as effective symbionts. Larger virus genomes show a greater departure from linkage equilibrium between core and satellites functions.

Inference of gene regulatory networks from time series by Tsallis entropy

Background: The inference of gene regulatory networks (GRNs) from large-scale expression profiles is one of the most challenging problems of Systems Biology nowadays. Many techniques and models have been proposed for this task. However, it is not generally possible to recover the original topology with great accuracy, mainly due to the short time series data in face of the high complexity of the networks and the intrinsic noise of the expression measurements. In order to improve the accuracy of GRNs inference methods based on entropy (mutual information), a new criterion function is here proposed. Results: In this paper we introduce the use of generalized entropy proposed by Tsallis, for the inference of GRNs from time series expression profiles. The inference process is based on a feature selection approach and the conditional entropy is applied as criterion function. In order to assess the proposed methodology, the algorithm is applied to recover the network topology from temporal expressions generated by an artificial gene network (AGN) model as well as from the DREAM challenge. The adopted AGN is based on theoretical models of complex networks and its gene transference function is obtained from random drawing on the set of possible Boolean functions, thus creating its dynamics. On the other hand, DREAM time series data presents variation of network size and its topologies are based on real networks. The dynamics are generated by continuous differential equations with noise and perturbation. By adopting both data sources, it is possible to estimate the average quality of the inference with respect to different network topologies, transfer functions and network sizes. Conclusions: A remarkable improvement of accuracy was observed in the experimental results by reducing the number of false connections in the inferred topology by the non-Shannon entropy. The obtained best free parameter of the Tsallis entropy was on average in the range 2.5 <= q <= 3.5 (hence, subextensive entropy), which opens new perspectives for GRNs inference methods based on information theory and for investigation of the nonextensivity of such networks. The inference algorithm and criterion function proposed here were implemented and included in the DimReduction software, which is freely available at http://sourceforge.net/projects/dimreduction and http://code.google.com/p/dimreduction/.

Functional and bioinformatics analyses reveal conservation of Cis-regulatory elements between sciaridae and Drosophilidae

The sciarid DNA puff C4 BhC4-1 gene is amplified and transcribed in salivary glands at the end of the larval stage. In transgenic Drosophila, the BhC4-1 promoter drives transcription in prepupal salivary glands and in the ring gland of late embryos. A bioinformatics analysis has identified 162 sequences similar to distinct regions of the BhC4-1 proximal promoter, which are predominantly located either in 5` or 3` regions or introns in the Drosophila melanogaster genome. A significant number of the identified sequences are found in the regulatory regions of Drosophila genes that are expressed in the salivary gland. Functional assays in Drosophila reveal that the BhC4-1 proximal promoter contains both a 129 bp (-186/-58) salivary gland enhancer and a 67 bp (-253/-187) ring gland enhancer that drive tissue specific patterns of developmentally regulated gene expression, irrespective of their orientation.

Comprehensive gene expression profiling in lungs of mice infected with Mycobacterium tuberculosis following DNAhsp65 immunotherapy

Background The continued increase in tuberculosis (TB) rates and the appearance of extremely resistant Mycobacterium tuberculosis strains (XDR-TB) worldwide are some of the great problems of public health. In this context, DNA immunotherapy has been proposed as an effective alternative that could circumvent the limitations of conventional drugs. Nonetheless, the molecular events underlying these therapeutic effects are poorly understood. Methods We characterized the transcriptional signature of lungs from mice infected with M. tuberculosis and treated with heat shock protein 65 as a genetic vaccine (DNAhsp65) combining microarray and real-time polymerase chain reaction analysis. The gene expression data were correlated with the histopathological analysis of lungs. Results The differential modulation of a high number of genes allowed us to distinguish DNAhsp65-treated from nontreated animals (saline and vector-injected mice). Functional analysis of this group of genes suggests that DNAhsp65 therapy could not only boost the T helper (Th)1 immune response, but also could inhibit Th2 cytokines and regulate the intensity of inflammation through fine tuning of gene expression of various genes, including those of interleukin-17, lymphotoxin A, tumour necrosis factor-cl, interleukin-6, transforming growth factor-beta, inducible nitric oxide synthase and Foxp3. In addition, a large number of genes and expressed sequence tags previously unrelated to DNA-therapy were identified. All these findings were well correlated with the histopathological lesions presented in the lungs. Conclusions The effects of DNA therapy are reflected in gene expression modulation; therefore, the genes identified as differentially expressed could be considered as transcriptional biomarkers of DNAhsp65 immunotherapy against TB. The data have important implications for achieving a better understanding of gene-based therapies. Copyright (C) 2008 John Wiley & Sons, Ltd.

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.

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.