Intrinsic noise in gene regulatory networks

Cells are intrinsically noisy biochemical reactors: low reactant numbers can lead to significant statistical fluctuations in molecule numbers and reaction rates. Here we use an analytic model to investigate the emergent noise properties of genetic systems. We find for a single gene that noise is essentially determined at the translational level, and that the mean and variance of protein concentration can be independently controlled. The noise strength immediately following single gene induction is almost twice the final steady-state value. We find that fluctuations in the concentrations of a regulatory protein can propagate through a genetic cascade; translational noise control could explain the inefficient translation rates observed for genes encoding such regulatory proteins. For an autoregulatory protein, we demonstrate that negative feedback efficiently decreases system noise. The model can be used to predict the noise characteristics of networks of arbitrary connectivity. The general procedure is further illustrated for an autocatalytic protein and a bistable genetic switch. The analysis of intrinsic noise reveals biological roles of gene network structures and can lead to a deeper understanding of their evolutionary origin.

A Late Mitotic Regulatory Network Controlling Cyclin Destruction in Saccharomyces cerevisiae

Exit from mitosis requires the inactivation of mitotic cyclin-dependent kinase–cyclin complexes, primarily by ubiquitin-dependent cyclin proteolysis. Cyclin destruction is regulated by a ubiquitin ligase known as the anaphase-promoting complex (APC). In the budding yeast Saccharomyces cerevisiae, members of a large class of late mitotic mutants, including cdc15, cdc5, cdc14, dbf2, and tem1, arrest in anaphase with a phenotype similar to that of cells expressing nondegradable forms of mitotic cyclins. We addressed the possibility that the products of these genes are components of a regulatory network that governs cyclin proteolysis. We identified a complex array of genetic interactions among these mutants and found that the growth defect in most of the mutants is suppressed by overexpression of SPO12, YAK1, and SIC1 and is exacerbated by overproduction of the mitotic cyclin Clb2. When arrested in late mitosis, the mutants exhibit a defect in cyclin-specific APC activity that is accompanied by high Clb2 levels and low levels of the anaphase inhibitor Pds1. Mutant cells arrested in G1 contain normal APC activity. We conclude that Cdc15, Cdc5, Cdc14, Dbf2, and Tem1 cooperate in the activation of the APC in late mitosis but are not required for maintenance of that activity in G1.

Epithelial and endothelial expression of the green fluorescent protein reporter gene under the control of bovine prion protein (PrP) gene regulatory sequences in transgenic mice

The expression of the cellular form of the prion protein (PrPc) gene is required for prion replication and neuroinvasion in transmissible spongiform encephalopathies. The identification of the cell types expressing PrPc is necessary to understanding how the agent replicates and spreads from peripheral sites to the central nervous system. To determine the nature of the cell types expressing PrPc, a green fluorescent protein reporter gene was expressed in transgenic mice under the control of 6.9 kb of the bovine PrP gene regulatory sequences. It was shown that the bovine PrP gene is expressed as two populations of mRNA differing by alternative splicing of one 115-bp 5′ untranslated exon in 17 different bovine tissues. The analysis of transgenic mice showed reporter gene expression in some cells that have been identified as expressing PrP, such as cerebellar Purkinje cells, lymphocytes, and keratinocytes. In addition, expression of green fluorescent protein was observed in the plexus of the enteric nervous system and in a restricted subset of cells not yet clearly identified as expressing PrP: the epithelial cells of the thymic medullary and the endothelial cells of both the mucosal capillaries of the intestine and the renal capillaries. These data provide valuable information on the distribution of PrPc at the cellular level and argue for roles of the epithelial and endothelial cells in the spread of infection from the periphery to the brain. Moreover, the transgenic mice described in this paper provide a model that will allow for the study of the transcriptional activity of the PrP gene promoter in response to scrapie infection.

Scavenger receptor A gene regulatory elements target gene expression to macrophages and to foam cells of atherosclerotic lesions.

Transcription of the macrophage scavenger receptor A gene is markedly upregulated during monocyte to macrophage differentiation. In these studies, we demonstrate that 291 bp of the proximal scavenger receptor promoter, in concert with a 400-bp upstream enhancer element, is sufficient to direct macrophage-specific expression of a human growth hormone reporter in transgenic mice. These regulatory elements, which contain binding sites for PU.1, AP-1, and cooperating ets-domain transcription factors, are also sufficient to mediate regulation of transgene expression during the in vitro differentiation of bone marrow progenitor cells in response to macrophage colony-stimulating factor. Mutation of the PU.1 binding site within the scavenger receptor promoter severely impairs transgene expression, consistent with a crucial role of PU.1 in regulating the expression of the scavenger receptor gene. The ability of the scavenger receptor promoter and enhancer to target gene expression to macrophages in vivo, including foam cells of atherosclerotic lesions, suggests that these regulatory elements will be of general utility in the study of macrophage differentiation and function by permitting specific modifications of macrophage gene expression.

DBTBS: a database of Bacillus subtilis promoters and transcription factors

With the completion of the determination of its entire genome sequence, one of the next major targets of Bacillus subtilis genomics is to clarify the whole gene regulatory network. To this end, the results of systematic experiments should be compared with the rich source of individual experimental results accumulated so far. Thus, we constructed a database of the upstream regulatory information of B.subtilis (DBTBS). The current version was constructed by surveying 291 references and contains information on 90 binding factors and 403 promoters. For each promoter, all of its known cis-elements are listed according to their positions, while these cis-elements are aligned to illustrate their consensus sequence for each transcription factor. All probable transcription factors coded in the genome were classified with the Pfam motifs. Using this database, we compared the character of B.subtilis promoters with that of Escherichia coli promoters. Our database is accessible at http://elmo.ims.u-tokyo.ac.jp/dbtbs/.

Integration of circadian and phototransduction pathways in the network controlling CAB gene transcription in Arabidopsis

The transcription of CAB genes, encoding the chlorophyll a/b-binding proteins, is rapidly induced in dark-grown Arabidopsis seedlings following a light pulse. The transient induction is followed by several cycles of a circadian rhythm. Seedlings transferred to continuous light are known to exhibit a robust circadian rhythm of CAB expression. The precise waveform of CAB expression in light–dark cycles, however, reflects a regulatory network that integrates information from photoreceptors, from the circadian clock and possibly from a developmental program. We have used the luciferase reporter system to investigate CAB expression with high time resolution. We demonstrate that CAB expression in light-grown plants exhibits a transient induction following light onset, similar to the response in dark-grown seedlings. The circadian rhythm modulates the magnitude and the kinetics of the response to light, such that the CAB promoter is not light responsive during the subjective night. A signaling pathway from the circadian oscillator must therefore antagonize the phototransduction pathways controlling the CAB promoter. We have further demonstrated that the phase of maximal CAB expression is delayed in light–dark cycles with long photoperiods, due to the entrainment of the circadian oscillator. Under short photoperiods, this pattern of entrainment ensures that dawn coincides with a phase of high light responsiveness, whereas under long photoperiods, the light response at dawn is reduced.

The TRANSFAC system on gene expression regulation

The TRANSFAC database on transcription factors and their DNA-binding sites and profiles (http://www.gene-regulation.de/) has been quantitatively extended and supplemented by a number of modules. These modules give information about pathologically relevant mutations in regulatory regions and transcription factor genes (PathoDB), scaffold/matrix attached regions (S/MARt DB), signal transduction (TRANSPATH) and gene expression sources (CYTOMER). Altogether, these distinct database modules constitute the TRANSFAC system. They are accompanied by a number of program routines for identifying potential transcription factor binding sites or for localizing individual components in the regulatory network of a cell.

Iron regulatory protein 1 is not required for the modulation of ferritin and transferrin receptor expression by iron in a murine pro-B lymphocyte cell line

Iron regulatory proteins (IRPs) are cytoplasmic RNA binding proteins that are central components of a sensory and regulatory network that modulates vertebrate iron homeostasis. IRPs regulate iron metabolism by binding to iron responsive element(s) (IREs) in the 5′ or 3′ untranslated region of ferritin or transferrin receptor (TfR) mRNAs. Two IRPs, IRP1 and IRP2, have been identified previously. IRP1 exhibits two mutually exclusive functions as an RNA binding protein or as the cytosolic isoform of aconitase. We demonstrate that the Ba/F3 family of murine pro-B lymphocytes represents the first example of a mammalian cell line that fails to express IRP1 protein or mRNA. First, all of the IRE binding activity in Ba/F3-gp55 cells is attributable to IRP2. Second, synthesis of IRP2, but not of IRP1, is detectable in Ba/F3-gp55 cells. Third, the Ba/F3 family of cells express IRP2 mRNA at a level similar to other murine cell lines, but IRP1 mRNA is not detectable. In the Ba/F3 family of cells, alterations in iron status modulated ferritin biosynthesis and TfR mRNA level over as much as a 20- and 14-fold range, respectively. We conclude that IRP1 is not essential for regulation of ferritin or TfR expression by iron and that IRP2 can act as the sole IRE-dependent mediator of cellular iron homeostasis.

Molecular Interaction Map of the Mammalian Cell Cycle Control and DNA Repair Systems

Eventually to understand the integrated function of the cell cycle regulatory network, we must organize the known interactions in the form of a diagram, map, and/or database. A diagram convention was designed capable of unambiguous representation of networks containing multiprotein complexes, protein modifications, and enzymes that are substrates of other enzymes. To facilitate linkage to a database, each molecular species is symbolically represented only once in each diagram. Molecular species can be located on the map by means of indexed grid coordinates. Each interaction is referenced to an annotation list where pertinent information and references can be found. Parts of the network are grouped into functional subsystems. The map shows how multiprotein complexes could assemble and function at gene promoter sites and at sites of DNA damage. It also portrays the richness of connections between the p53-Mdm2 subsystem and other parts of the network.

UV light selectively coinduces supply pathways from primary metabolism and flavonoid secondary product formation in parsley

The UV light-induced synthesis of UV-protective flavonoids diverts substantial amounts of substrates from primary metabolism into secondary product formation and thus causes major perturbations of the cellular homeostasis. Results from this study show that the mRNAs encoding representative enzymes from various supply pathways are coinduced in UV-irradiated parsley cells (Petroselinum crispum) with two mRNAs of flavonoid glycoside biosynthesis, encoding phenylalanine ammonia-lyase and chalcone synthase. Strong induction was observed for mRNAs encoding glucose 6-phosphate dehydrogenase (carbohydrate metabolism, providing substrates for the shikimate pathway), 3-deoxyarabinoheptulosonate 7-phosphate synthase (shikimate pathway, yielding phenylalanine), and acyl-CoA oxidase (fatty acid degradation, yielding acetyl-CoA), and moderate induction for an mRNA encoding S-adenosyl-homocysteine hydrolase (activated methyl cycle, yielding S-adenosyl-methionine for B-ring methylation). Ten arbitrarily selected mRNAs representing various unrelated metabolic activities remained unaffected. Comparative analysis of acyl-CoA oxidase and chalcone synthase with respect to mRNA expression modes and gene promoter structure and function revealed close similarities. These results indicate a fine-tuned regulatory network integrating those functionally related pathways of primary and secondary metabolism that are specifically required for protective adaptation to UV irradiation. Although the response of parsley cells to UV light is considerably broader than previously assumed, it contrasts greatly with the extensive metabolic reprogramming observed previously in elicitor-treated or fungus-infected cells.

Mice with a targeted mutation in the thyroid hormone β receptor gene exhibit impaired growth and resistance to thyroid hormone

Patients with mutations in the thyroid hormone receptor β (TRβ) gene manifest resistance to thyroid hormone (RTH), resulting in a constellation of variable phenotypic abnormalities. To understand the molecular basis underlying the action of mutant TRβ in vivo, we generated mice with a targeted mutation in the TRβ gene (TRβPV; PV, mutant thyroid hormone receptor kindred PV) by using homologous recombination and the Cre/loxP system. Mice expressing a single PVallele showed the typical abnormalities of thyroid function found in heterozygous humans with RTH. Homozygous PV mice exhibit severe dysfunction of the pituitary–thyroid axis, impaired weight gains, and abnormal bone development. This phenotype is distinct from that seen in mice with a null mutation in the TRβ gene. Importantly, we identified abnormal expression patterns of several genes in tissues of TRβPV mice, demonstrating the interference of the mutant TR with the gene regulatory functions of the wild-type TR in vivo. These results show that the actions of mutant and wild-type TRβ in vivo are distinct. This model allows further study of the molecular action of mutant TR in vivo, which could lead to better treatment for RTH patients.

The Ath5 proneural genes function upstream of Brn3 POU domain transcription factor genes to promote retinal ganglion cell development

During retinogenesis, the Xenopus basic helix–loop–helix transcription factor Xath5 has been shown to promote a ganglion cell fate. In the developing mouse and chicken retinas, gene targeting and overexpression studies have demonstrated critical roles for the Brn3 POU domain transcription factor genes in the promotion of ganglion cell differentiation. However, the genetic relationship between Ath5 and Brn3 genes is unknown. To understand the genetic regulatory network(s) that controls retinal ganglion cell development, we analyzed the relationship between Ath5 and Brn3 genes by using a gain-of-function approach in the chicken embryo. We found that during retinogenesis, the chicken Ath5 gene (Cath5) is expressed in retinal progenitors and in differentiating ganglion cells but is absent in terminally differentiated ganglion cells. Forced expression of both Cath5 and the mouse Ath5 gene (Math5) in retinal progenitors activates the expression of cBrn3c following central-to-peripheral and temporal-to-nasal gradients. As a result, similar to the Xath5 protein, both Cath5 and Math5 proteins have the ability to promote the development of ganglion cells. Moreover, we found that forced expression of all three Brn3 genes also can stimulate the expression of cBrn3c. We further found that Ath5 and Brn3 proteins are capable of transactivating a Brn3b promoter. Thus, these data suggest that the expression of cBrn3c in the chicken and Brn3b in the mouse is initially activated by Ath5 factors in newly generated ganglion cells and later maintained by a feedback loop of Brn3 factors in the differentiated ganglion cells.

Multilocus analysis of extracellular putative virulence proteins made by group A Streptococcus: Population genetics, human serologic response, and gene transcription

Species of pathogenic microbes are composed of an array of evolutionarily distinct chromosomal genotypes characterized by diversity in gene content and sequence (allelic variation). The occurrence of substantial genetic diversity has hindered progress in developing a comprehensive understanding of the molecular basis of virulence and new therapeutics such as vaccines. To provide new information that bears on these issues, 11 genes encoding extracellular proteins in the human bacterial pathogen group A Streptococcus identified by analysis of four genomes were studied. Eight of the 11 genes encode proteins with a LPXTG(L) motif that covalently links Gram-positive virulence factors to the bacterial cell surface. Sequence analysis of the 11 genes in 37 geographically and phylogenetically diverse group A Streptococcus strains cultured from patients with different infection types found that recent horizontal gene transfer has contributed substantially to chromosomal diversity. Regions of the inferred proteins likely to interact with the host were identified by molecular population genetic analysis, and Western immunoblot analysis with sera from infected patients confirmed that they were antigenic. Real-time reverse transcriptase–PCR (TaqMan) assays found that transcription of six of the 11 genes was substantially up-regulated in the stationary phase. In addition, transcription of many genes was influenced by the covR and mga trans-acting gene regulatory loci. Multilocus investigation of putative virulence genes by the integrated approach described herein provides an important strategy to aid microbial pathogenesis research and rapidly identify new targets for therapeutics research.

DNA bending by an adenine–thymine tract and its role in gene regulation

To gain insight into the structural basis of DNA bending by adenine–thymine tracts (A-tracts) and their role in DNA recognition by gene-regulatory proteins, we have determined the crystal structure of the high-affinity DNA target of the cancer-associated human papillomavirus E2 protein. The three independent B-DNA molecules of the crystal structure determined at 2.2-Å resolution are examples of A-tract-containing helices where the global direction and magnitude of curvature are in accord with solution data, thereby providing insights, at the base pair level, into the mechanism of DNA bending by such sequence motifs. A comparative analysis of E2–DNA conformations with respect to other structural and biochemical studies demonstrates that (i) the A-tract structure of the core region, which is not contacted by the protein, is critical for the formation of the high-affinity sequence-specific protein–DNA complex, and (ii) differential binding affinity is regulated by the intrinsic structure and deformability encoded in the base sequence of the DNA target.

Yeast histone H3 and H4 N termini function through different GAL1 regulatory elements to repress and activate transcription.

Previous work has shown that N-terminal deletions of yeast histone H3 cause a 2- to 4-fold increase in the induction of GAL1 and a number of other genes involved in galactose metabolism. In contrast, deletions at the H4 N terminus cause a 10- to 20-fold decrease in the induction of these same GAL genes. However, H3 and H4 N-terminal deletions each decrease PHO5 induction only 2- to 4-fold. To define the GAL1 gene regulatory elements through which the histone N termini activate or repress transcription, fusions were made between GAL1 and PHO5 promoter elements attached to a beta-galactosidase reporter gene. We show here that GAL1 hyperactivation caused by the H3 N-terminal deletion delta 4-15 is linked to the upstream activation sequence. Conversely, the relative decrease in GAL1 induction caused by the H4N-terminal deletion delta 4-28 is linked to the downstream promoter which contains the TATA element. These data indicate that the H3 N terminus is required for the repression of the GAL1 upstream element, whereas the H4N terminus is required for the activation of the GAL1 downstream promoter element.