Modularity and evolutionary constraints in a baculovirus gene regulatory network

Abstract Background The structure of regulatory networks remains an open question in our understanding of complex biological systems. Interactions during complete viral life cycles present unique opportunities to understand how host-parasite network take shape and behave. The Anticarsia gemmatalis multiple nucleopolyhedrovirus (AgMNPV) is a large double-stranded DNA virus, whose genome may encode for 152 open reading frames (ORFs). Here we present the analysis of the ordered cascade of the AgMNPV gene expression. Results We observed an earlier onset of the expression than previously reported for other baculoviruses, especially for genes involved in DNA replication. Most ORFs were expressed at higher levels in a more permissive host cell line. Genes with more than one copy in the genome had distinct expression profiles, which could indicate the acquisition of new functionalities. The transcription gene regulatory network (GRN) for 149 ORFs had a modular topology comprising five communities of highly interconnected nodes that separated key genes that are functionally related on different communities, possibly maximizing redundancy and GRN robustness by compartmentalization of important functions. Core conserved functions showed expression synchronicity, distinct GRN features and significantly less genetic diversity, consistent with evolutionary constraints imposed in key elements of biological systems. This reduced genetic diversity also had a positive correlation with the importance of the gene in our estimated GRN, supporting a relationship between phylogenetic data of baculovirus genes and network features inferred from expression data. We also observed that gene arrangement in overlapping transcripts was conserved among related baculoviruses, suggesting a principle of genome organization. Conclusions Albeit with a reduced number of nodes (149), the AgMNPV GRN had a topology and key characteristics similar to those observed in complex cellular organisms, which indicates that modularity may be a general feature of biological gene regulatory networks.

A novel phase variation mechanism in the meningococcus driven by a ligand-responsive repressor and differential spacing of distal promoter elements

Phase variable expression, mediated by high frequency reversible changes in the length of simple sequence repeats, facilitates adaptation of bacterial populations to changing environments and is frequently important in bacterial virulence. Here we elucidate a novel phase variable mechanism for NadA expression, an adhesin and invasin of Neisseria meningitidis. The NadR repressor protein binds to operators flanking the phase variable tract of the nadA promoter gene and contributes to the differential expression levels of phase variant promoters with different numbers of repeats, likely due to different spacing between operators. It is shown that IHF binds between these operators, and may permit looping of the promoter, allowing interaction of NadR at operators located distally or overlapping the promoter. The 4-hydroxyphenylacetic acid, a metabolite of aromatic amino acid catabolism that is secreted in saliva, induces nadA expression by inhibiting the DNA binding activity of the NadR repressor. When induced, only minor differences are evident between NadR-independent transcription levels of promoter phase variants, which are likely due to differential RNA polymerase contacts leading to altered promoter activity. These results suggest that NadA expression is under both stochastic and tight environmental-sensing regulatory control, and both regulations are mediated by the NadR repressor that and may be induced during colonization of the oropharynx where it plays a major role in the successful adhesion and invasion of the mucosa. Hence, simple sequence repeats in promoter regions may be a strategy used by host-adapted bacterial pathogens to randomly switch between expression states that may nonetheless still be induced by appropriate niche-specific signals.

An interferon regulatory factor element controls the major immune modulator of murine cytomegalovirus

Immune modulation by herpesviruses, such as cytomegalovirus, is critical for the establishment of acute and persistent infection confronting a vigorous antiviral immune response of the host. Therefore, the action of immune-modulatory proteins has long been the subject of research, with the final goal to identify new strategies for antiviral therapy.rnIn the case of murine cytomegalovirus (mCMV), the viral m152 protein has been identified to play a major role in targeting components of both the innate and the adaptive immune system in terms of infected host-cell recognition in the effector phase of the antiviral immune response. On the one hand, it inhibits cell surface expression of RAE-1 and thereby prevents ligation of the activating natural killer (NK)-cell receptor NKG2D. On the other hand, it decreases cell surface expression of peptide-loaded MHC class I molecules thereby preventing antigen presentation to CD8 T cells. Ultimately, the outcome of CMV infection is determined by the interplay between viral and cellular factors.rnIn this context, the work presented here has revealed a novel and intriguing connection between viral m152 and cellular interferon (IFN), a key cytokine of the immune system: rnthe m152 promoter region contains an interferon regulatory factor element (IRFE) perfectly matching the consensus sequence of cellular IRFEs.rnThe biological relevance of this regulatory element was first suggested by sequence comparisons revealing its evolutionary conservation among various established laboratory strains of mCMV and more recent low-passage wild-derived virus isolates. Moreover, search of the mCMV genome revealed only three IRFE sites in the complete sequence. Importantly, the functionality of the IRFE in the m152 promoter was confirmed with the use of a mutant virus, representing a functional deletion of the IRFE, and its corresponding revertant virus. In particular, m152 gene expression was found to be inhibited in an IRFE-dependent manner in infected cells. Essentially, this inhibition proved to have a severe impact on the immune-modulatory function of m152, first demonstrated by a restored direct antigen presentation on infected cells for CD8 T-cell activation. Even more importantly, this effect of IRFE-mediated IFN signaling was validated in vivo by showing that the protective antiviral capacity of adoptively-transferred, antigen-specific CD8 T cells is also significantly restored by the IRFE-dependent inhibition of m152. Somewhat curious and surprising, the decrease in m152 protein simultaneously prevented an enhanced activation of NK cells in acute-infected mice, apparently independent of the RAE-1/NKG2D ligand/receptor interaction but rather due to reduced ‘missing-self’ recognition.rnTaken together, this work presents a so far unknown mechanism of IFN signaling to control mCMV immune modulation in acute infection.rnrn

Design and construction of a new Drosophila species, D.synthetica, by synthetic regulatory evolution

Here, I merge the principles of synthetic biology1,2 and regulatory evolution3-11 to create a new species12-15 with a minimal set of known elements. Using preexisting transgenes and recessive mutations of Drosophila melanogaster, a transgenic population arises with small eyes and a different venation pattern that fulfills the criteria of a new species according to Mayr's "Biological Species Concept"7,10. The genetic circuit entails the loss of a non-essential transcription factor and the introduction of cryptic enhancers. Subsequent activation of those enhancers causes hybrid lethality. The transition from "transgenic organisms" towards "synthetic species", such as Drosophila synthetica, constitutes a safety mechanism to avoid hybridization with wild type populations and preserve natural biodiversity16-18. Drosophila synthetica is the first transgenic organism that cannot hybridize with the original wild type population but remains fertile when crossed with other transgenic animals.

Cis-acting elements and developmental regulators govern toxin gene expression in Bacillus anthracis

Expression of the structural genes for the anthrax toxin proteins is coordinately controlled by host-related signals such as elevated CO2 , and the trans-acting positive regulator, AtxA. No specific binding of AtxA to the toxin gene promoters has been demonstrated and no sequence-based similarities are apparent in the promoter regions of toxin genes. We hypothesized that the toxin genes possess common structural features that are required for positive regulation. To test this hypothesis, I performed an extensive characterization of the toxin gene promoters. I determined the minimal sequences required for atxA-mediated toxin gene expression and compared these sequences for structural similarities. In silico modeling and in vitro experiments indicated significant curvature within these regions. Random mutagenesis revealed that point mutations associated with reduced transcriptional activity, mostly mapped to areas of high curvature. This work enabled the identification of two potential cis-acting elements implicated in AtxA-mediated regulation of the toxin genes. In addition to the growth condition requirements and AtxA, toxin gene expression is under growth phase regulation. The transition state regulator AbrB represses atxA expression to influence toxin synthesis. Here I report that toxin gene expression also requires sigH, a gene encoding the RNA polymerase sigma factor associated with development in B. subtilis. In the well-studied B. subtilis system, σH is part of a feedback control pathway that involves AbrB and the major response regulator of sporulation initiation, Spo0A. My data indicate that in B. anthracis, regulatory relationships exist between these developmental regulators and atxA . Interestingly, during growth in toxin-inducing conditions, sigH and abrB expression deviates from that described for B. subtilis, affecting expression of the atxA gene. These findings, combined with previous observations, suggest that the steady state level of atxA expression is critical for optimal toxin gene transcription. I propose a model whereby, under toxin-inducing conditions, control of toxin gene expression is fine-tuned by the independent effects of the developmental regulators on the expression of atxA . The growth condition-dependent changes in expression of these regulators may be crucial for the correct timing and uninterrupted expression of the toxin genes during infection. ^

Control of the master virulence regulatory gene atxA in Bacillus anthracis

Transcription of the Bacillus anthracis structural genes for the anthrax toxin proteins and biosynthetic operon for capsule are positively regulated by AtxA, a transcription regulator with unique properties. Consistent with the role of atxA in virulence factor expression, a B. anthracis atxA-null mutant is avirulent in a murine model for anthrax. In batch culture, multiple signals impact atxA transcript levels, and the timing and steady state level of atxA expression is critical for optimal toxin and capsule synthesis. Despite the apparent complex control of atxA transcription, only one trans-acting protein, the transition state regulator AbrB, has been demonstrated to directly interact with the atxA promoter. The AbrB-binding site has been described, but additional cis-acting control sequences have not been defined. Using transcriptional lacZ fusions, electrophoretic mobility shift assays, and Western blot analysis, the cis-acting elements and trans-acting factors involved in regulation of atxA in B. anthracis strains containing either both virulence plasmids, pXO1 and pXO2, or only one plasmid, pXO1, were studied. This work demonstrates that atxA transcription from the major start site P1 is dependent upon a consensus sequence for the housekeeping sigma factor SigA, and an A+T-rich upstream element (UP-element) for RNA polymerase (RNAP). In addition, the data show that a trans-acting protein(s) other than AbrB negatively impacts atxA transcription when it binds specifically to a 9-bp palindrome within atxA promoter sequences located downstream of P1. Mutation of the palindrome prevents binding of the trans-acting protein(s) and results in a corresponding increase in AtxA and anthrax toxin production in a strain- and culture-dependent manner. The identity of the trans-acting repressor protein(s) remains elusive; however, phenotypes associated with mutation of the repressor binding site have revealed that the trans-acting repressor protein(s) indirectly controls B. anthracis development. Mutation of the repressor binding site results in misregulation and overexpression of AtxA in conditions conducive for development, leading to a marked sporulation defect that is both atxA- and pXO2-61-dependent. pXO2-61 is homologous to the sensor domain of sporulation sensor histidine kinases and is proposed to titrate an activating signal away from the sporulation phosphorelay when overexpressed by AtxA. These results indicate that AtxA is not only a master virulence regulator, but also a modulator of proper B. anthracis development. Also demonstrated in this work is the impact of the developmental regulators AbrB, Spo0A, and SigH on atxA expression and anthrax toxin production in a genetically incomplete (pXO1+, pXO2-) and genetically complete (pXO1+, pXO2+) strain background. AtxA and anthrax toxin production resulting from deletion of the developmental regulators are strain-dependent suggesting that factors on pXO2 are involved in control of atxA. The only developmental deletion mutant that resulted in a prominent and consistent strain-independent increase in AtxA protein levels was an abrB-null mutant. As a result of increased AtxA levels, there is early and increased production of anthrax toxins in an abrB-null mutant. In addition, the abrB-null mutant exhibited an increase in virulence in a murine model for anthrax. In contrast, virulence of the atxA promoter mutant was unaffected in a murine model for anthrax despite the production of 5-fold more AtxA than the abrB-null mutant. These results imply that AtxA is not the only factor impacting pathogenesis in an abrB-null mutant. Overall, this work highlights the complex regulatory network that governs expression of atxA and provides an additional role for AtxA in B. anthracis development.

Sterol regulatory element binding protein-1c is a major mediator of insulin action on the hepatic expression of glucokinase and lipogenesis-related genes

Hepatic glucokinase plays a key role in glucose metabolism as underlined by the anomalies associated with glucokinase mutations and the consequences of tissue-specific knock-out. In the liver, glucokinase transcription is absolutely dependent on the presence of insulin. The cis-elements and trans-acting factors that mediate the insulin effect are presently unknown; this is also the case for most insulin-responsive genes. We have shown previously that the hepatic expression of the transcription factor sterol regulatory element binding protein-1c (SREBP-1c) is activated by insulin. We show here in primary cultures of hepatocytes that the adenovirus-mediated transduction of a dominant negative form of SREBP-1c inhibits the insulin effect on endogenous glucokinase expression. Conversely, in the absence of insulin, the adenovirus-mediated transduction of a dominant positive form of SREBP-1c overcomes the insulin dependency of glucokinase expression. Hepatic fatty acid synthase and Spot-14 are insulin/glucose-dependent genes. For this latter class of genes, the dominant positive form of SREBP-1c obviates the necessity for the presence of insulin, whereas glucose potentiates the effect of SREBP-1c on their expression. In addition, the insulin dependency of lipid accumulation in cultured hepatocytes is overcome by the dominant positive form of SREBP-1c. We propose that SREBP-1c is a major mediator of insulin action on hepatic gene expression and a key regulator of hepatic glucose/lipid metabolism.

Multiple domains of repressor activator protein 1 contribute to facilitated binding of glycolysis regulatory protein 1

The function of repressor activator protein 1 (Rap1p) at glycolytic enzyme gene upstream activating sequence (UAS) elements in Saccharomyces cerevisiae is to facilitate binding of glycolysis regulatory protein 1 (Gcr1p) at adjacent sites. Rap1p has a modular domain structure. In its amino terminus there is an asymmetric DNA-bending domain, which is distinct from its DNA-binding domain, which resides in the middle of the protein. In the carboxyl terminus of Rap1p lie its silencing and putative activation domains. We carried out a molecular dissection of Rap1p to identify domains contributing to its ability to facilitate binding of Gcr1p. We prepared full-length and three truncated versions of Rap1p and tested their ability to facilitate binding of Gcr1p by gel shift assay. The ability to detect ternary complexes containing Rap1p⋅DNA⋅Gcr1p depended on the presence of binding sites for both proteins in the probe DNA. The DNA-binding domain of Rap1p, although competent to bind DNA, was unable to facilitate binding of Gcr1p. Full-length Rap1p and the amino- and carboxyl-truncated versions of Rap1p were each able to facilitate binding of Gcr1p at an appropriately spaced binding site. Under these conditions, Gcr1p displayed an approximately 4-fold greater affinity for Rap1p-bound DNA than for otherwise identical free DNA. When spacing between Rap1p- and Gcr1p-binding sites was altered by insertion of five nucleotides, the ability to form ternary Rap1p⋅DNA⋅Gcr1p complexes was inhibited by all but the DNA-binding domain of Rap1p itself; however, the ability of each individual protein to bind the DNA probe was unaffected.

Functional antioxidant responsive elements

Exposure of human and rodent cells to a wide variety of chemoprotective compounds confers resistance against a broad set of carcinogens. For a subset of the chemoprotective compounds, protection is generated by an increase in the abundance of protective enzymes like glutathione S-transferases (GST). Antioxidant responsive elements (AREs) mediate the transcriptional induction of a battery of genes which comprise much of this chemoprotective response system. Past studies identified a necessary ARE “core” sequence of RTGACnnnGC, but this sequence alone is insufficient to mediate induction. In this study, the additional sequences necessary to define a sufficient, functional ARE are identified through systematic mutational analysis of the murine GST Ya ARE. Introduction of the newly identified necessary nucleotides into the regions flanking a nonresponsive, ARE-like, GST-Mu promoter sequence produced an inducible element. A screen of the GenBank database with the newly identified ARE consensus identified 16 genes which contained the functional ARE consensus sequence in their promoters. Included within this group was an ARE sequence from the murine ferritin-L promoter that mediated induction when tested. In an electrophoretic mobility-shift assay, the ferritin-L ARE was bound by ARE–binding protein 1, a protein previously identified as the likely mediator of the chemoprotective response. A three-level ARE classification system is presented to account for the distinct induction strengths observed in our mutagenesis studies. A model of the ARE as a composite regulatory site, where multiple transcription factors interact, is presented to account for the complex characteristics of ARE-mediated chemoprotective gene expression.

Major recent and independent changes in levels and patterns of expression have occurred at the b gene, a regulatory locus in maize

The b locus encodes a transcription factor that regulates the expression of genes that produce purple anthocyanin pigment. Different b alleles are expressed in distinct tissues, causing tissue-specific anthocyanin production. Understanding how phenotypic diversity is produced and maintained at the b locus should provide models for how other regulatory genes, including those that influence morphological traits and development, evolve. We have investigated how different levels and patterns of pigmentation have evolved by determining the phenotypic and evolutionary relationships between 18 alleles that represent the diversity of b alleles in Zea mays. Although most of these alleles have few phenotypic differences, five alleles have very distinct tissue-specific patterns of pigmentation. Superimposing the phenotypes on the molecular phylogeny reveals that the alleles with strong and distinctive patterns of expression are closely related to alleles with weak expression, implying that the distinctive patterns have arisen recently. We have identified apparent insertions in three of the five phenotypically distinct alleles, and the fourth has unique upstream restriction fragment length polymorphisms relative to closely related alleles. The insertion in B-Peru has been shown to be responsible for its unique expression and, in the other two alleles, the presence of the insertion correlates with the phenotype. These results suggest that major changes in gene expression are probably the result of large-scale changes in DNA sequence and/or structure most likely mediated by transposable elements.

Silencer elements controlling the B29 (Igβ) promoter are neither promoter- nor cell-type-specific

The murine B29 (Igβ) promoter is B cell specific and contains essential SP1, ETS, OCT, and Ikaros motifs. Flanking 5′ DNA sequences inhibit B29 promoter activity, suggesting this region contains silencer elements. Two adjacent 5′ DNA segments repress transcription by the murine B29 promoter in a position- and orientation-independent manner, analogous to known silencers. Both these 5′ segments also inhibit transcription by several heterologous promoters in B cells, including mb-1, c-fos, and human B29. These 5′ segments also inhibit transcription by the c-fos promoter in T cells suggesting they are not B cell-specific elements. DNase I footprint analyses show an approximately 70-bp protected region overlapping the boundary between the two negative regulatory DNA segments and corresponding to binding sites for at least two different DNA-binding proteins. Within this footprint, two unrelated 30-bp cis-acting DNA motifs (designated TOAD and FROG) function as position- and orientation-independent silencers when located directly 5′ of the murine B29 promoter. These two silencer motifs act cooperatively to restrict the transcriptional activity of the B29 promoter. Neither of these motifs resembles any known silencers. Mutagenesis of the TOAD and FROG motifs in their respective 5′ DNA segments eliminates the silencing activity of these upstream regions, indicating these two motifs as the principal B29 silencer elements within these regions.

Restriction-modification gene complexes as selfish gene entities: Roles of a regulatory system in their establishment, maintenance, and apoptotic mutual exclusion

We have reported some type II restriction-modification (RM) gene complexes on plasmids resist displacement by an incompatible plasmid through postsegregational host killing. Such selfish behavior may have contributed to the spread and maintenance of RM systems. Here we analyze the role of regulatory genes (C), often found linked to RM gene complexes, in their interaction with the host and the other RM gene complexes. We identified the C gene of EcoRV as a positive regulator of restriction. A C mutation eliminated postsegregational killing by EcoRV. The C system has been proposed to allow establishment of RM systems in new hosts by delaying the appearance of restriction activity. Consistent with this proposal, bacteria preexpressing ecoRVC were transformed at a reduced efficiency by plasmids carrying the EcoRV RM gene complex. Cells carrying the BamHI RM gene complex were transformed at a reduced efficiency by a plasmid carrying a PvuII RM gene complex, which shares the same C specificity. The reduction most likely was caused by chromosome cleavage at unmodified PvuII sites by prematurely expressed PvuII restriction enzyme. Therefore, association of the C genes of the same specificity with RM gene complexes of different sequence specificities can confer on a resident RM gene complex the capacity to abort establishment of a second, incoming RM gene complex. This phenomenon, termed “apoptotic mutual exclusion,” is reminiscent of suicidal defense against virus infection programmed by other selfish elements. pvuIIC and bamHIC genes define one incompatibility group of exclusion whereas ecoRVC gene defines another.

Activation of iron regulatory protein-1 by oxidative stress in vitro

Iron regulatory protein-1 (IRP-1), a central cytoplasmic regulator of cellular iron metabolism, is rapidly activated by oxidative stress to bind to mRNA iron-responsive elements. We have reconstituted the response of IRP-1 to extracellular H2O2 in a system derived from murine B6 fibroblasts permeabilized with streptolysin-O. This procedure allows separation of the cytosol from the remainder of the cells (cell pellet). IRP-1 in the cytosolic fraction fails to be directly activated by addition of H2O2. IRP-1 activation requires the presence of a nonsoluble, possibly membrane-associated component in the cell pellet. The streptolysin-O-based in vitro system faithfully recapitulates characteristic hallmarks of IRP-1 activation by H2O2 in intact cells. We show that the H2O2-mediated activation of IRP-1 is temperature dependent and sensitive to treatment with calf intestinal alkaline phosphatase (CIAP). Although IRP-1 activation is unaffected by addition of excess ATP or GTP to this in vitro system, it is negatively affected by the nonhydrolyzable nucleotide analogs adenylyl-imidodiphosphate and guanylyl-imidophosphate and completely blocked by ATP-γS and GTP-γS. The in vitro reconstitution of this oxidative stress-induced pathway has opened a different avenue for the biochemical dissection of the regulation of mammalian iron metabolism by oxidative stress. Our data show that H2O2 must be sensed to stimulate a pathway to activate IRP-1.

Novel role of phosphorylation in Fe–S cluster stability revealed by phosphomimetic mutations at Ser-138 of iron regulatory protein 1

Animals regulate iron metabolism largely through the action of the iron regulatory proteins (IRPs). IRPs modulate mRNA utilization by binding to iron-responsive elements (IRE) in the 5′ or 3′ untranslated region of mRNAs encoding proteins involved in iron homeostasis or energy production. IRP1 is also the cytosolic isoform of aconitase. The activities of IRP1 are mutually exclusive and are modulated through the assembly/disassembly of its [4Fe–4S] cluster, reversibly converting it between an IRE-binding protein and cytosolic aconitase. IRP1 is also phosphoregulated by protein kinase C, but the mechanism by which phosphorylation posttranslationally increases IRE binding activity has not been fully defined. To investigate this, Ser-138 (S138), a PKC phosphorylation site, was mutated to phosphomimetic glutamate (S138E), aspartate (S138D), or nonphosphorylatable alanine (S138A). The S138E IRP1 mutant and, to a lesser extent, the S138D IRP1 mutant were impaired in aconitase function in yeast when grown aerobically but not when grown anaerobically. Purified wild-type and mutant IRP1s could be reconstituted to active aconitases anaerobically. However, when exposed to oxygen, the [4Fe–4S] cluster of the S138D and S138E mutants decayed 5-fold and 20-fold faster, respectively, than was observed for wild-type IRP1. Our findings suggest that stability of the Fe–S cluster of IRP1 can be regulated by phosphorylation and reveal a mechanism whereby the balance between the IRE binding and [4Fe–4S] forms of IRP1 can be modulated independently of cellular iron status. Furthermore, our results show that IRP1 can function as an oxygen-modulated posttranscriptional regulator of gene expression.

Polymorphism in the 3′-untranslated region of TNFα mRNA impairs binding of the post-transcriptional regulatory protein HuR to TNFα mRNA

Tumor necrosis factor α (TNFα) acts as a beneficial mediator in the process of host defence. In recent years major interest has focused on the AU-rich elements (AREs) present in the 3′-untranslated region (3′-UTR) of TNFα mRNA as this region plays a pivotal role in post-transcriptional control of TNFα production. Certain stimuli, such as lipopolysaccharides, a component of the Gram-negative bacterial cell wall, have the ability to relinquish the translational suppression of TNFα mRNA imposed by these AREs in macrophages, thereby enabling the efficient production of the TNFα. In this study we show that the polymorphism (GAU trinucleotide insertional mutation) present in the regulatory 3′-UTR of TNFα mRNA of NZW mice results in the hindered binding of RNA-binding proteins, thereby leading to a significantly reduced production of TNFα protein. We also show that the binding of macrophage proteins to the main ARE is also decreased by another trinucleotide (CAU) insertion in the TNFα 3′-UTR. One of the proteins affected by the GAU trinucleotide insertional mutation was identified as HuR, a nucleo-cytoplasmic shuttling protein previously shown to play a prominent role in the stability and translatability of mRNA containing AREs. Since binding of this protein most likely modulates the stability, translational efficiency and transport of TNFα mRNA, these results suggest that mutations in the ARE of TNFα mRNA decrease the production of TNFα protein in macrophages by hindering the binding of HuR to the ARE.