Compilation of transcription regulating proteins

Developmental genome-wide occupancy analysis of bZIP transcription factor NRL uncovers the role of c-Jun in early differentiation of rod photoreceptors in the mammalian retina

The basic motif-leucine zipper (bZIP) transcription factor NRL determines rod photoreceptor cell fate during retinal development, and its loss leads to cone-only retina in mice. NRL works synergistically with homeodomain protein CRX and other regulatory factors to control the transcription of most genes associated with rod morphogenesis and functional maturation, which span over a period of several weeks in the mammalian retina. We predicted that NRL gradually establishes rod cell identity and function by temporal and dynamic regulation of stage-specific transcriptional targets. Therefore, we mapped the genomic occupancy of NRL at four stages of mouse photoreceptor differentiation by CUT&RUN analysis. Dynamics of NRL-binding revealed concordance with the corresponding changes in transcriptome of the developing rods. Notably, we identified c-Jun proto-oncogene as one of the targets of NRL, which could bind to specific cis-elements in the c-Jun promoter and modulate its activity in HEK293 cells. Coimmunoprecipitation studies showed association of NRL with c-Jun, also a bZIP protein, in transfected cells as well as in developing mouse retina. Additionally, shRNA-mediated knockdown of c-Jun in the mouse retina in vivo resulted in altered expression of almost 1000 genes, with reduced expression of phototransduction genes and many direct targets of NRL in rod photoreceptors. We propose that c-Jun-NRL heterodimers prime the NRL-directed transcriptional program in neonatal rod photoreceptors before high NRL expression suppresses c-Jun at later stages. Our study highlights a broader cooperation among cell-type restricted and widely expressed bZIP proteins, such as c-Jun, in specific spatiotemporal contexts during cellular differentiation.

A New Weighted Imputed Neighborhood-Regularized Tri-Factorization One-Class Collaborative Filtering Algorithm: Application to Target Gene Prediction of Transcription Factors

Identifying target genes of transcription factors (TFs) is crucial to understand transcriptional regulation. However, our understanding of genome-wide TF targeting profile is limited due to the cost of large-scale experiments and intrinsic complexity of gene regulation. Thus, computational prediction methods are useful to predict unobserved TF-gene associations. Here, we develop a new Weighted Imputed Neighborhood-regularized Tri-Factorization one-class collaborative filtering algorithm, WINTF. It predicts unobserved target genes for TFs using known but noisy, incomplete, and biased TF-gene associations and protein-protein interaction networks. Our benchmark study shows that WINTF significantly outperforms its counterpart matrix factorization-based algorithms and tri-factorization methods that do not include weight, imputation, and neighbor-regularization, for TF-gene association prediction. When evaluated by independent datasets, accuracy is 37.8 percent on the top 495 predicted associations, an enrichment factor of 4.19 compared with random guess. Furthermore, many predicted novel associations are supported by literature evidence. Although we only use canonical TF-gene interaction data, WINTF can directly be applied to tissue-specific data when available. Thus, WINTF provides a potentially useful framework to integrate multiple omics data for further improvement of TF-gene prediction and applications to other sparse and noisy biological data. The benchmark dataset and source code are freely available at https://github.com/XieResearchGroup/WINTF.

The histone demethylase KDM6B fine-tunes the host response to Streptococcus pneumoniae

Streptococcus pneumoniae is a natural colonizer of the human respiratory tract and an opportunistic pathogen. Although epithelial cells are among the first to encounter pneumococci, the cellular processes and contribution of epithelial cells to the host response are poorly understood. Here, we show that a S. pneumoniae serotype 6B ST90 strain, which does not cause disease in a murine infection model, induces a unique NF-κB signature response distinct from an invasive-disease-causing isolate of serotype 4 (TIGR4). This signature is characterized by activation of p65 and requires a histone demethylase KDM6B. We show, molecularly, that the interaction of the 6B strain with epithelial cells leads to chromatin remodelling within the IL-11 promoter in a KDM6B-dependent manner, where KDM6B specifically demethylates histone H3 lysine 27 dimethyl. Remodelling of the IL-11 locus facilitates p65 access to three NF-κB sites that are otherwise inaccessible when stimulated by IL-1β or TIGR4. Finally, we demonstrate through chemical inhibition of KDM6B with GSK-J4 inhibitor and through exogenous addition of IL-11 that the host responses to the 6B ST90 and TIGR4 strains can be interchanged both in vitro and in a murine model of infection in vivo. Our studies therefore reveal how a chromatin modifier governs cellular responses during infection.

Exercise rejuvenates quiescent skeletal muscle stem cells in old mice through restoration of Cyclin D1

Aging impairs tissue repair. This is pronounced in skeletal muscle, whose regeneration by muscle stem cells (MuSCs) is robust in young adult animals but inefficient in older organisms. Despite this functional decline, old MuSCs are amenable to rejuvenation through strategies that improve the systemic milieu, such as heterochronic parabiosis. One such strategy, exercise, has long been appreciated for its benefits on healthspan, but its effects on aged stem cell function in the context of tissue regeneration are incompletely understood. Here we show that exercise in the form of voluntary wheel running accelerates muscle repair in old animals and improves old MuSC function. Through transcriptional profiling and genetic studies, we discovered that the restoration of old MuSC activation ability hinges on restoration of Cyclin D1, whose expression declines with age in MuSCs. Pharmacologic studies revealed that Cyclin D1 maintains MuSC activation capacity by repressing TGFβ signaling. Taken together, these studies demonstrate that voluntary exercise is a practicable intervention for old MuSC rejuvenation. Furthermore, this work highlights the distinct role of Cyclin D1 in stem cell quiescence.

Expanding Role of Ubiquitin in Translational Control

The eukaryotic proteome has to be precisely regulated at multiple levels of gene expression, from transcription, translation, and degradation of RNA and protein to adjust to several cellular conditions. Particularly at the translational level, regulation is controlled by a variety of RNA binding proteins, translation and associated factors, numerous enzymes, and by post-translational modifications (PTM). Ubiquitination, a prominent PTM discovered as the signal for protein degradation, has newly emerged as a modulator of protein synthesis by controlling several processes in translation. Advances in proteomics and cryo-electron microscopy have identified ubiquitin modifications of several ribosomal proteins and provided numerous insights on how this modification affects ribosome structure and function. The variety of pathways and functions of translation controlled by ubiquitin are determined by the various enzymes involved in ubiquitin conjugation and removal, by the ubiquitin chain type used, by the target sites of ubiquitination, and by the physiologic signals triggering its accumulation. Current research is now elucidating multiple ubiquitin-mediated mechanisms of translational control, including ribosome biogenesis, ribosome degradation, ribosome-associated protein quality control (RQC), and redox control of translation by ubiquitin (RTU). This review discusses the central role of ubiquitin in modulating the dynamism of the cellular proteome and explores the molecular aspects responsible for the expanding puzzle of ubiquitin signals and functions in translation.

Target Gene Prediction of Transcription Factor Using a New Neighborhood-regularized Tri-factorization One-class Collaborative Filtering Algorithm

Identifying the target genes of transcription factors (TFs) is one of the key factors to understand transcriptional regulation. However, our understanding of genome-wide TF targeting profile is limited due to the cost of large scale experiments and intrinsic complexity. Thus, computational prediction methods are useful to predict the unobserved associations. Here, we developed a new one-class collaborative filtering algorithm tREMAP that is based on regularized, weighted nonnegative matrix tri-factorization. The algorithm predicts unobserved target genes for TFs using known gene-TF associations and protein-protein interaction network. Our benchmark study shows that tREMAP significantly outperforms its counterpart REMAP, a bi-factorization-based algorithm, for transcription factor target gene prediction in all four performance metrics AUC, MAP, MPR, and HLU. When evaluated by independent data sets, the prediction accuracy is 37.8% on the top 495 predicted associations, an enrichment factor of 4.19 compared with the random guess. Furthermore, many of the predicted novel associations by tREMAP are supported by evidence from literature. Although we only use canonical TF-target gene interaction data in this study, tREMAP can be directly applied to tissue-specific data sets. tREMAP provides a framework to integrate multiple omics data for the further improvement of TF target gene prediction. Thus, tREMAP is a potentially useful tool in studying gene regulatory networks. The benchmark data set and the source code of tREMAP are freely available at https://github.com/hansaimlim/REMAP/tree/master/TriFacREMAP.

Motif scraper: a cross-platform, open-source tool for identifying degenerate nucleotide motif matches in FASTA files

Summary

Many genomic features are defined not by exact sequence matches, but by degenerate nucleotide motifs that represent multiple compatible matches. While there are databases cataloging genomic features, such as the location of transcription factor motifs, for commonly used model species, identifying the locations of novel motifs, known motifs in non-model genomes, or known motifs in personal whole-genomes is difficult. I designed motif scraper to overcome this limitation, allowing for efficient, multiprocessor motif searches in any FASTA file.

Availability and implementation

The motif scraper package (MIT license) is available via PyPI, and the Python source is available on GitHub at https://github.com/RobersonLab/motif_scraper.

Prediction and Validation of Transcription Factors Binding Sites in the Il9 Locus

Over the past decade, multiple effector T cell subsets have been identified with varying differentiation conditions in the milieu as well as a broad diversity of cytokine expression. Interleukin-9 (IL-9) secreting T helper 9 (Th9) cells are the newest member of this family. T helper cell differentiation including Th9 cells appears to be an epigenetic phenomenon requiring the coordination of a large variety of transcription factors to reshape the chromatin landscape and generate various T helper phenotypes. This chapter details methods for both predicting and validating potential transcription factor binding sites as well as their downstream epigenetic effect using a variety of in silico and in vitro methods in both primary Th9 cells and IL-9-producing T cell lines.

Transcription factors involved in drought tolerance and their possible role in developing drought tolerant cultivars with emphasis on wheat (Triticum aestivum L.)

TFs involved in drought tolerance in plants may be utilized in future for developing drought tolerant cultivars of wheat and some other crops. Plants have developed a fairly complex stress response system to deal with drought and other abiotic stresses. These response systems often make use of transcription factors (TFs); a gene encoding a specific TF together with -its target genes constitute a regulon, and take part in signal transduction to activate/silence genes involved in response to drought. Since, five specific families of TFs (out of >80 known families of TFs) have gained widespread attention on account of their significant role in drought tolerance in plants, TFs and regulons belonging to these five multi-gene families (AP2/EREBP, bZIP, MYB/MYC, NAC and WRKY) have been described and their role in improving drought tolerance discussed in this brief review. These TFs often undergo reversible phosphorylation to perform their function, and are also involved in complex networks. Therefore, some details about reversible phosphorylation of TFs by different protein kinases/phosphatases and the co-regulatory networks, which involve either only TFs or TFs with miRNAs, have also been discussed. Literature on transgenics involving genes encoding TFs and that on QTLs and markers associated with TF genes involved in drought tolerance has also been reviewed. Throughout the review, there is a major emphasis on wheat as an important crop, although examples from the model cereal rice (sometimes maize also), and the model plant Arabidopsis have also been used. This knowledge base may eventually allow the use of TF genes for development of drought tolerant cultivars, particularly in wheat.

CG-containing oligonucleotides and transcription factor-binding motifs are enriched in human pericentric regions

Unsupervised data mining capable of extracting a wide range of information from big sequence data without prior knowledge or particular models is highly desirable in an era of big data accumulation for research on genes, genomes and genetic systems. By handling oligonucleotide compositions in genomic sequences as high-dimensional data, we have previously modified the conventional SOM (self-organizing map) for genome informatics and established BLSOM for oligonucleotide composition, which can analyze more than ten million sequences simultaneously and is thus suitable for big data analyses. Oligonucleotides often represent motif sequences responsible for sequence-specific binding of proteins such as transcription factors. The distribution of such functionally important oligonucleotides is probably biased in genomic sequences, and may differ among genomic regions. When constructing BLSOMs to analyze pentanucleotide composition in 50-kb sequences derived from the human genome in this study, we found that BLSOMs did not classify human sequences according to chromosome but revealed several specific zones, which are enriched for a class of CG-containing pentanucleotides; these zones are composed primarily of sequences derived from pericentric regions. The biological significance of enrichment of these pentanucleotides in pericentric regions is discussed in connection with cell type- and stage-dependent formation of the condensed heterochromatin in the chromocenter, which is formed through association of pericentric regions of multiple chromosomes.

Bioinformatic landscapes for plant transcription factor system research

Diverse bioinformatic resources have been developed for plant transcription factor (TF) research. This review presents the bioinformatic resources and methodologies for the elucidation of plant TF-mediated biological events. Such information is helpful to dissect the transcriptional regulatory systems in the three reference plants Arabidopsis , rice, and maize and translation to other plants. Transcription factors (TFs) orchestrate diverse biological programs by the modulation of spatiotemporal patterns of gene expression via binding cis-regulatory elements. Advanced sequencing platforms accompanied by emerging bioinformatic tools revolutionize the scope and extent of TF research. The system-level integration of bioinformatic resources is beneficial to the decoding of TF-involved networks. Herein, we first briefly introduce general and specialized databases for TF research in three reference plants Arabidopsis, rice, and maize. Then, as proof of concept, we identified and characterized heat shock transcription factor (HSF) members through the TF databases. Finally, we present how the integration of bioinformatic resources at -omics layers can aid the dissection of TF-mediated pathways. We also suggest ways forward to improve the bioinformatic resources of plant TFs. Leveraging these bioinformatic resources and methodologies opens new avenues for the elucidation of transcriptional regulatory systems in the three model systems and translation to other plants.

Identification of a face enhancer reveals direct regulation of LIM homeobox 8 (Lhx8) by wingless-int (WNT)/β-catenin signaling

Development of the mammalian face requires a large number of genes that are expressed with spatio-temporal specificity, and transcriptional regulation mediated by enhancers plays a key role in the precise control of gene expression. Using chromatin immunoprecipitation for a histone marker of active enhancers, we generated a genome-wide map of candidate enhancers from the maxillary arch (primordium for the upper jaw) of mouse embryos. Furthermore, we confirmed multiple novel craniofacial enhancers near the genes implicated in human palate defects through functional assays. We characterized in detail one of the enhancers (Lhx8_enh1) located upstream of Lhx8, a key regulatory gene for craniofacial development. Lhx8_enh1 contained an evolutionarily conserved binding site for lymphoid enhancer factor/T-cell factor family proteins, which mediate the transcriptional regulation by the WNT/β-catenin signaling pathway. We demonstrated in vitro that WNT/β-catenin signaling was indeed essential for the expression of Lhx8 in the maxillary arch cells and that Lhx8_enh1 was a direct target of the WNT/β-catenin pathway. Together, we uncovered a molecular mechanism for the regulation of Lhx8, and we provided valuable resources for further investigation into the gene regulatory network of craniofacial development.

Chromatin and epigenetic determinants of estrogen receptor alpha (ESR1) signaling

The oestrogen receptor alpha (ESR1) is a transcription factor that potentiates the response to diverse stimuli, including oestrogen and growth factors, in various tissue types. Its recruitment to the DNA is directly regulated by the chromatin landscape, inclusive of chromatin compaction and epigenetic modifications. In this review we discuss our current understanding of the interplay between ESR1 signaling and the chromatin landscape. We present how the chromatin landscape primes the lineage-specific response and contributes to stimuli-specific signaling. Finally, we discuss recent efforts to decipher the relationship between genetic and epigenetic as it relates to ESR1 signaling in breast cancer.

Quantifying the information content of homing endonuclease target sites by single base pair profiling

Homing endonucleases (HEs) are native proteins that recognize long DNA sequences with high site specificity in vitro and in vivo. The target site specificity of HEs is high, though not absolute. For example, members of the well-characterized LAGLIDADG family of homing endonucleases (the LHEs) recognize target sites of ~20 base pairs, and can tolerate some target site base pair changes without losing site binding or cleavage activity. This modest degree of target site degeneracy is practically useful once defined and can facilitate the engineering of LHE variants with new DNA recognition specificities. In this chapter, we outline general protocols for systematically profiling HE target site base pair positions in order to define their functional importance in vitro and in vivo, and show how information theory can be used to make sense of the resulting data.

The placental mammal ancestor and the post-K-Pg radiation of placentals

To discover interordinal relationships of living and fossil placental mammals and the time of origin of placentals relative to the Cretaceous-Paleogene (K-Pg) boundary, we scored 4541 phenomic characters de novo for 86 fossil and living species. Combining these data with molecular sequences, we obtained a phylogenetic tree that, when calibrated with fossils, shows that crown clade Placentalia and placental orders originated after the K-Pg boundary. Many nodes discovered using molecular data are upheld, but phenomic signals overturn molecular signals to show Sundatheria (Dermoptera + Scandentia) as the sister taxon of Primates, a close link between Proboscidea (elephants) and Sirenia (sea cows), and the monophyly of echolocating Chiroptera (bats). Our tree suggests that Placentalia first split into Xenarthra and Epitheria; extinct New World species are the oldest members of Afrotheria.

TFClass: an expandable hierarchical classification of human transcription factors

TFClass (http://tfclass.bioinf.med.uni-goettingen.de/) provides a comprehensive classification of human transcription factors based on their DNA-binding domains. Transcription factors constitute a large functional family of proteins directly regulating the activity of genes. Most of them are sequence-specific DNA-binding proteins, thus reading out the information encoded in cis-regulatory DNA elements of promoters, enhancers and other regulatory regions of a genome. TFClass is a database that classifies human transcription factors by a six-level classification schema, four of which are abstractions according to different criteria, while the fifth level represents TF genes and the sixth individual gene products. Altogether, nine superclasses have been identified, comprising 40 classes and 111 families. Counted by genes, 1558 human TFs have been classified so far or >2900 different TFs when including their isoforms generated by alternative splicing or protein processing events. With this classification, we hope to provide a basis for deciphering protein-DNA recognition codes; moreover, it can be used for constructing expanded transcriptional networks by inferring additional TF-target gene relations.

Sequences of acetyl CoA carboxylase promoter for tumour necrosis factor action

Tumour necrosis factor (TNF) inhibits the accumulation of acetyl CoA carboxylase (ACC) mRNA by decreasing the rate of ACC gene transcription. The ACC mRNA species found in 30A5 cells are generated from promoter II and TNF inhibits the accumulation of class 2 type mRNAs. By using 5' deletion mutants of promoter II fused to the bacterial chloramphenicol acetyltransferase (CAT) gene, the DNA mobility shift assay and the DNase I footprinting assay, the authors have identified the 30 bp from −389 to −359 as the TNF responsive element in promoter II. TNF treatment causes a decrease in the binding activity of nuclear protein(s) specific to the TNF responsive element. When the fragment containing the TNF responsive element was incorporated into the thymidine kinase promoter, the chimeric gene exhibited TNF induced inhibition of expression.

Allele-specific transcription of the asthma-associated PHD finger protein 11 gene (PHF11) modulated by octamer-binding transcription factor 1 (Oct-1)

BACKGROUND

Asthma is a common, chronic inflammatory airway disease of major public health importance with multiple genetic determinants. Previously, we found by positional cloning that PHD finger protein 11 (PHF11) on chromosome 13q14 modifies serum immunoglobulin E (IgE) concentrations and asthma susceptibility. No coding variants in PHF11 were identified.

OBJECTIVE

Here we investigate the 3 single nucleotide polymorphisms (SNPs) in this gene most significantly associated with total serum IgE levels--rs3765526, rs9526569, and rs1046295--for a role in transcription factor binding.

METHODS

We used electrophoretic mobility shift assays to examine the effect of the 3 SNPs on transcription factor binding in 3 cell lines relevant to asthma pathogenesis. Relative preferential expression of alleles was investigated by using the allelotyping method.

RESULTS

Electrophoretic mobility shift assays show that rs1046295 modulates allele-specific binding by the octamer-binding transcription factor 1 (Oct-1). Analysis of the relative expression levels of the 2 alleles of this SNP in heterozygous individuals showed a modest, but highly significant (P = 6.5 × 10(-16)), preferential expression of the A allele consistent with a functional role for rs1046295.

CONCLUSION

These results suggest a mechanism by which rs1046295 may act as a regulatory variant modulating transcription at this locus and altering asthma susceptibility.

Up regulation of GW112 Gene by NF kappaB promotes an antiapoptotic property in gastric cancer cells

To clarify the regulatory mechanism of GW112 gene expression, 5'-flanking region of the human GW112 gene was isolated and characterized in the present study. 5'-RACE analysis showed a single transcription start site, which is located 142 nucleotides upstream of the translation initiation site. Transient transfection studies with serial deletion constructs and close examination of the sequences identified a putative NF kappaB binding sequence between -442 and -430, which could be responsible for efficient expression of the GW112 gene. Indeed, GW112 gene was found to be regulated by NF kappaB signals including overexpressed p65 and I kappaB alpha, IKK inhibitor, and proteasome inhibitor. Binding of NF kappaB to its putative site was confirmed by EMSA and ChIP assays. These results suggest that NF kappaB is an essential regulatory factor for GW112 transcription. Based on this finding, we next confirmed that inhibition of GW112 expression could induce apoptosis in the presence of cytotoxic agent in gastric cancer cells. Furthermore, knocking-down or overexpression of GW112 gene in gastric cancer cells demonstrated that GW112 has an antiapoptotic property against the cytotoxic agents-induced apoptosis. Taken together, these results suggest that GW112 could be an important mediator in NF kappaB-dependent tumorigenesis of digestive tract tissues.

Identification of promoter regions and regulatory sites

Promoter sequences are the main regulatory elements of gene expression. Their recognition by computer algorithms is fundamental for understanding gene expression patterns, cell specificity and development. This chapter describes the advanced approaches to identify promoters in animal, plant and bacterial sequences. Also, we discuss an approach to identify statistically significant regulatory motifs in genomic sequences.

Regulation by transcription factors in bacteria: beyond description

Transcription is an essential step in gene expression and its understanding has been one of the major interests in molecular and cellular biology. By precisely tuning gene expression, transcriptional regulation determines the molecular machinery for developmental plasticity, homeostasis and adaptation. In this review, we transmit the main ideas or concepts behind regulation by transcription factors and give just enough examples to sustain these main ideas, thus avoiding a classical ennumeration of facts. We review recent concepts and developments: cis elements and trans regulatory factors, chromosome organization and structure, transcriptional regulatory networks (TRNs) and transcriptomics. We also summarize new important discoveries that will probably affect the direction of research in gene regulation: epigenetics and stochasticity in transcriptional regulation, synthetic circuits and plasticity and evolution of TRNs. Many of the new discoveries in gene regulation are not extensively tested with wetlab approaches. Consequently, we review this broad area in Inference of TRNs and Dynamical Models of TRNs. Finally, we have stepped backwards to trace the origins of these modern concepts, synthesizing their history in a timeline schema.

Prediction of Cis-Regulatory Elements Controlling Genes Differentially Expressed by Retinal and Choroidal Vascular Endothelial Cells

Cultured endothelial cells of the human retina and choroid demonstrate distinct patterns of gene expression. We hypothesized that differential gene expression reflected differences in the interactions of transcription factors and respective cis-regulatory motifs(s) in these two endothelial cell subpopulations, recognizing that motifs often exist as modules. We tested this hypothesis in silico by using TRANSFAC Professional and CisModule to identify cis-regulatory motifs and modules in genes that were differentially expressed by human retinal versus choroidal endothelial cells, as identified by analysis of a microarray data set. Motifs corresponding to eight transcription factors were significantly (p < 0.05) differentially abundant in genes that were relatively highly expressed in retinal (i.e., glucocorticoid receptor, high mobility group AT-hook 1, heat shock transcription factor 1, p53, vitamin D receptor) or choroidal (i.e., transcription factor E2F, Yin Yang 1, zinc finger 5) endothelial cells. Predicted cis-regulatory modules were quite different for these two groups of genes. Our findings raise the possibility of exploiting specific cis-regulatory motifs to target therapy at the ocular endothelial cells subtypes responsible for neovascular age-related macular degeneration or proliferative diabetic retinopathy.

Identification of direct target genes using joint sequence and expression likelihood with application to DAF-16

A major challenge in the post-genome era is to reconstruct regulatory networks from the biological knowledge accumulated up to date. The development of tools for identifying direct target genes of transcription factors (TFs) is critical to this endeavor. Given a set of microarray experiments, a probabilistic model called TRANSMODIS has been developed which can infer the direct targets of a TF by integrating sequence motif, gene expression and ChIP-chip data. The performance of TRANSMODIS was first validated on a set of transcription factor perturbation experiments (TFPEs) involving Pho4p, a well studied TF in Saccharomyces cerevisiae. TRANSMODIS removed elements of arbitrariness in manual target gene selection process and produced results that concur with one's intuition. TRANSMODIS was further validated on a genome-wide scale by comparing it with two other methods in Saccharomyces cerevisiae. The usefulness of TRANSMODIS was then demonstrated by applying it to the identification of direct targets of DAF-16, a critical TF regulating ageing in Caenorhabditis elegans. We found that 189 genes were tightly regulated by DAF-16. In addition, DAF-16 has differential preference for motifs when acting as an activator or repressor, which awaits experimental verification. TRANSMODIS is computationally efficient and robust, making it a useful probabilistic framework for finding immediate targets.

Complete genome and phylogenetic position of bovine papillomavirus type 7

Six bovine papillomavirus (BPV) types and 16 putative BPV types have been reported previously. Here, the complete genome sequence of BAPV6, a novel putative BPV type isolated from cattle in Japan, was determined by using multiple-primed rolling-circle amplification. The genome consisted of 7412 bp (G+C content of 46 mol%) that encoded five early (E1, E2, E4, E6 and E7) and two late (L1 and L2) genes, but did not encode the E5 gene. The E6 protein contained a non-consensus CxxC(x)33CxxC and a consensus CxxC(x)29CxxC zinc-binding domain, and the E7 protein lacked the LxCxE motif. The nucleotide sequence of the L1 open reading frame (ORF) was related most closely (57–58 %) to the L1 ORF of member(s) of the genera Betapapillomavirus, Gammapapillomavirus and Pipapillomavirus. Phylogenetic analysis based on the complete L1 ORF suggests that BAPV6 should be classified in a novel genus in the family Papillomaviridae as BPV-7.

Integrative content-driven concepts for bioinformatics “beyond the cell”

Bioinformatics has delivered great contributions to genome and genomics research, without which the world-wide success of this and other global ('omics') approaches would not have been possible. More recently, it has developed further towards the analysis of different kinds of networks thus laying the foundation for comprehensive description, analysis and manipulation of whole living systems in modern "systems biology". The next step which is necessary for developing a systems biology that deals with systemic phenomena is to expand the existing and develop new methodologies that are appropriate to characterize intercellular processes and interactions without omitting the causal underlying molecular mechanisms. Modelling the processes on the different levels of complexity involved requires a comprehensive integration of information on gene regulatory events, signal transduction pathways, protein interaction and metabolic networks as well as cellular functions in the respective tissues / organs.

In vitro analysis of DNA-protein interactions by proximity ligation

Protein-binding DNA sequence elements encode a variety of regulated functions of genomes. Information about such elements is currently in a state of rapid growth, but improved methods are required to characterize the sequence specificity of DNA-binding proteins. We have established an in vitro method for specific and sensitive solution-phase analysis of interactions between proteins and nucleic acids in nuclear extracts, based on the proximity ligation assay. The reagent consumption is very low, and the excellent sensitivity of the assay enables analysis of as few as 1-10 cells. We show that our results are highly reproducible, quantitative, and in good agreement with both EMSA and predictions obtained by using a motif finding software. This assay can be a valuable tool to characterize in-depth the sequence specificity of DNA-binding proteins and to evaluate effects of polymorphisms in known transcription factor binding sites.

Sequence analysis and organization of the Neodiprion abietis nucleopolyhedrovirus genome

Of 30 baculovirus genomes that have been sequenced to date, the only nonlepidopteran baculoviruses include the dipteran Culex nigripalpus nucleopolyhedrovirus and two hymenopteran nucleopolyhedroviruses that infect the sawflies Neodiprion lecontei (NeleNPV) and Neodiprion sertifer (NeseNPV). This study provides a complete sequence and genome analysis of the nucleopolyhedrovirus that infects the balsam fir sawfly Neodiprion abietis (Hymenoptera, Symphyta, Diprionidae). The N. abietis nucleopolyhedrovirus (NeabNPV) is 84,264 bp in size, with a G+C content of 33.5%, and contains 93 predicted open reading frames (ORFs). Eleven predicted ORFs are unique to this baculovirus, 10 ORFs have a putative sequence homologue in the NeleNPV genome but not the NeseNPV genome, and 1 ORF (neab53) has a putative sequence homologue in the NeseNPV genome but not the NeleNPV genome. Specific repeat sequences are coincident with major genome rearrangements that distinguish NeabNPV and NeleNPV. Genes associated with these repeat regions encode a common amino acid motif, suggesting that they are a family of repeated contiguous gene clusters. Lepidopteran baculoviruses, similarly, have a family of repeated genes called the bro gene family. However, there is no significant sequence similarity between the NeabNPV and bro genes. Homologues of early-expressed genes such as ie-1 and lef-3 were absent in NeabNPV, as they are in the previously sequenced hymenopteran baculoviruses. Analyses of ORF upstream sequences identified potential temporally distinct genes on the basis of putative promoter elements.

Studies on the activity of a bidirectional promoter of Mungbean yellow mosaic India virus by agroinfiltration

The AV promoter expressing AV1 and AV2 genes and AC1 promoter expressing AC1 gene are present in opposite orientation in the intergenic region of Mungbean yellow mosaic India virus (MYMIV). Transient Agrobacterium-mediated delivery of putative promoter constructs into Nicotiana benthamiana and different legumes, followed by reporter gene (beta-d-glucuronidase, GUS) assay, identified the promoter region of both AC1 and AV genes that is necessary for transcriptional initiation. Transcription activator protein-independent activity of AV promoter and differential regulation of AC1 promoter are unique to MYMIV. The AV promoter is a composite core promoter having both TATA box and Initiator elements (TATA(+)Inr(+)). Many transcription factor binding sites were identified in the upstream promoter sequences of both virion and complementary sense genes, which might be used in the transcription regulation studies of the host plant as well as the virus.

TRANSFAC and its module TRANSCompel: transcriptional gene regulation in eukaryotes

The TRANSFAC database on transcription factors, their binding sites, nucleotide distribution matrices and regulated genes as well as the complementing database TRANSCompel on composite elements have been further enhanced on various levels. A new web interface with different search options and integrated versions of Match and Patch provides increased functionality for TRANSFAC. The list of databases which are linked to the common GENE table of TRANSFAC and TRANSCompel has been extended by: Ensembl, UniGene, EntrezGene, HumanPSD and TRANSPRO. Standard gene names from HGNC, MGI and RGD, are included for human, mouse and rat genes, respectively. With the help of InterProScan, Pfam, SMART and PROSITE domains are assigned automatically to the protein sequences of the transcription factors. TRANSCompel contains now, in addition to the COMPEL table, a separate table for detailed information on the experimental EVIDENCE on which the composite elements are based. Finally, for TRANSFAC, in respect of data growth, in particular the gain of Drosophila transcription factor binding sites (by courtesy of the Drosophila DNase I footprint database) and of Arabidopsis factors (by courtesy of DATF, Database of Arabidopsis Transcription Factors) has to be stressed. The here described public releases, TRANSFAC 7.0 and TRANSCompel 7.0, are accessible under http://www.gene-regulation.com/pub/databases.html.

Binding sites for metabolic disease related transcription factors inferred at base pair resolution by chromatin immunoprecipitation and genomic microarrays

We present a detailed in vivo characterization of hepatocyte transcriptional regulation in HepG2 cells, using chromatin immunoprecipitation and detection on PCR fragment-based genomic tiling path arrays covering the encyclopedia of DNA element (ENCODE) regions. Our data suggest that HNF-4alpha and HNF-3beta, which were commonly bound to distal regulatory elements, may cooperate in the regulation of a large fraction of the liver transcriptome and that both HNF-4alpha and USF1 may promote H3 acetylation to many of their targets. Importantly, bioinformatic analysis of the sequences bound by each transcription factor (TF) shows an over-representation of motifs highly similar to the in vitro established consensus sequences. On the basis of these data, we have inferred tentative binding sites at base pair resolution. Some of these sites have been previously found by in vitro analysis and some were verified in vitro in this study. Our data suggests that a similar approach could be used for the in vivo characterization of all predicted/uncharacterized TF and that the analysis could be scaled to the whole genome.

PromoLign: a database for upstream region analysis and SNPs

The study of transcriptional regulation at the genomic level has been hindered by the lack of functional annotation in the putative regulatory regions. Phylogenetic footprinting, in which cross-species sequence alignment among orthologous genes is applied to locate conserved sequence blocks, is an effective strategy to attack this problem. Single nucleotide polymorphisms (SNPs) in transcription factor (TF) binding sites contribute to the heterogeneity of TF binding sites and might disrupt or enhance their regulatory activity. The correlation of SNPs with the TF sites will not only help in functional evaluation of SNPs, but will also help in the study of transcription regulation by focusing attention on specific TF sites. PromoLign (http://polly.wustl.edu/promolign/main.html) is an online database application that presents SNPs and TF binding profiles in the context of human-mouse orthologous sequence alignment with a hyperlinked graphical interface. PromoLign could be applied to a variety of SNPs and transcription related studies, including association genetics, population genetics, and pharmacogenetics.

TRANSFAC: transcriptional regulation, from patterns to profiles

The TRANSFAC database on eukaryotic transcriptional regulation, comprising data on transcription factors, their target genes and regulatory binding sites, has been extended and further developed, both in number of entries and in the scope and structure of the collected data. Structured fields for expression patterns have been introduced for transcription factors from human and mouse, using the CYTOMER database on anatomical structures and developmental stages. The functionality of Match, a tool for matrix-based search of transcription factor binding sites, has been enhanced. For instance, the program now comes along with a number of tissue-(or state-)specific profiles and new profiles can be created and modified with Match Profiler. The GENE table was extended and gained in importance, containing amongst others links to LocusLink, RefSeq and OMIM now. Further, (direct) links between factor and target gene on one hand and between gene and encoded factor on the other hand were introduced. The TRANSFAC public release is available at http://www.gene-regulation.com. For yeast an additional release including the latest data was made available separately as TRANSFAC Saccharomyces Module (TSM) at http://transfac.gbf.de. For CYTOMER free download versions are available at http://www.biobase.de:8080/index.html.

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.

Gene regulation of Wiskott-Aldrich syndrome protein and the human homolog of the Drosophila Su(var)3-9: WASP and SUV39H1, two adjacent genes at Xp11.23

The region Xp11.23 is a gene-rich, light giemsa-staining segment on the short arm of the X chromosome. In this study, we have characterized the transcriptional regulatory elements in this interval for two adjacent genes: SUV39H1, a regulator of chromatin organization, and the Wiskott-Aldrich syndrome protein (WASP). The WASP gene exhibits two alternate promoters, both of which demonstrate transcription factor binding elements specific to blood cell lineages. Reporter gene expression analyses indicate that both WASP promoters show high levels of expression in different hematopoietic cell lines. The human homolog of the Drosophila Su(var)3-9 gene was identified by sequence analysis of the region downstream from WASP. SUV39H1 is ubiquitously expressed, and the promoter sequence consists mostly of general transcription factors. The presence of putative binding sites for GAGA and Adf1 transcription factors may indicate a cross regulatory mechanism with other chromatin regulators.

TRANSFAC: an integrated system for gene expression regulation

TRANSFAC is a database on transcription factors, their genomic binding sites and DNA-binding profiles (http://transfac.gbf.de/TRANSFAC/). Its content has been enhanced, in particular by information about training sequences used for the construction of nucleotide matrices as well as by data on plant sites and factors. Moreover, TRANSFAC has been extended by two new modules: PathoDB provides data on pathologically relevant mutations in regulatory regions and transcription factor genes, whereas S/MARt DB compiles features of scaffold/matrix attached regions (S/MARs) and the proteins binding to them. Additionally, the databases TRANSPATH, about signal transduction, and CYTOMER, about organs and cell types, have been extended and are increasingly integrated with the TRANSFAC data sources.

Interaction of the human NF-kappaB p52 transcription factor with DNA-PNA hybrids mimicking the NF-kappaB binding sites of the human immunodeficiency virus type 1 promoter

We determined whether peptide nucleic acids (PNAs) are able to interact with NF-kappaB p52 transcription factor. The binding of NF-kappaB p52 to DNA-DNA, DNA-PNA, PNA-DNA, and PNA-PNA hybrid molecules carrying the NF-kappaB binding sites of human immunodeficiency type 1 long terminal repeat was studied by (i) biospecific interaction analysis (BIA) using surface plasmon resonance technology, (ii) electrophoretic mobility shift, (iii) DNase I footprinting, and (iv) UV cross-linking assays. Our results demonstrate that NF-kappaB p52 does not efficiently bind to PNA-PNA hybrids. However, a DNA-PNA hybrid molecule was found to be recognized by NF-kappaB p52, although the molecular complexes generated exhibited low stability. From the theoretical point of view, our results suggest that binding of NF-kappaB p52 protein to target DNA motifs is mainly due to contacts with bases; interactions with the DNA backbone are, however, important for stabilization of the protein-DNA complex. From the practical point of view, our results suggest that DNA-PNA hybrid can be recognized by NF-kappaB p52 protein, although with an efficiency lower than DNA-DNA NF-kappaB target molecules; therefore, our results should encourage studies on modified PNAs in order to develop potential agents for the decoy approach in gene therapy.

The identification and characterization of two promoters and the complete genomic sequence for the Wiskott-Aldrich syndrome gene

The Wiskott-Aldrich syndrome (WAS) is an X-linked disorder characterized by immunodeficiency, eczema and thrombocytopenia. The gene responsible for WAS was identified through positional cloning, and the function of the encoded protein (WASP) is still the subject of much speculation. WASP is currently thought to be involved in the regulation of actin polymerization in hematopoietic cells. To study the elements that regulate the WASP gene, we have identified the sites for transcription initiation. We found that two promoters were responsible for controlling WASP expression. Multiple transcription initiation sites were found immediately adjacent to the translation start site, however an alternate exon with a second promoter region was identified 6 kb upstream. Examination of the 5' sequence adjacent to the initiation sites in both promoters failed to reveal a TATA or CCAAT box, but numerous putative transcription factor binding sites including Sp1, Ets, c-Myb and PU.1 were apparent. Reporter constructs generated from each promoter showed functional activity in the Jurkat T-cell and HEL erythro-megakaryocytic cell lines. Although the alternate exon sequence was extremely GC rich and contained several potential binding elements, the primary promoter was stronger than the upstream promoter in the cell lines assayed. The transcription factor binding site profiles within each promoter suggested that they may play different roles in regulating WASP expression depending on the stage of differentiation and development, and the cell lineage. In this study we have also reported the complete nucleotide sequence of the coding and intervening sequences for the WASP gene. A comprehensive knowledge of the genomic structure and the further characterization of WASP gene expression will facilitate the continued investigation of mutations in WAS patients, and the eventual prospect of gene therapy.

Expanding the TRANSFAC database towards an expert system of regulatory molecular mechanisms

TRANSFAC is a database on transcription factors, their genomic binding sites and DNA-binding profiles. In addition to being updated and extended by new features, it has been complemented now by a series of additional database modules. Among them, modules which provide data about signal transduction pathways (TRANSPATH) or about cell types/organs/developmental stages (CYTOMER) are available as well as an updated version of the previously described COMPEL database. The databases are available on the WWW at http://transfac.gbf.de/

Mechanism of angiotensin II-mediated regulation of fibronectin gene in rat vascular smooth muscle cells

This study was performed to investigate a mechanism of angiotensin II (Ang II)-mediated activation of the fibronectin (FN) gene in rat vascular smooth muscle cells. Actinomycin D and CV11974 completely inhibited Ang II-mediated increase in FN mRNA levels. Inhibitors of protein kinase C (PKC), protein-tyrosine kinase (PTK), phosphatidylinositol-specific phospholipase C, Ras, phosphatidylinositol 3-kinase, p70 S6 kinase, and Ca2+/calmodulin kinase also decreased Ang II-induced activation of FN mRNA. In contrast, cycloheximide; PD123319; or inhibitors of Gi, protein kinase A, or mitogen-activated protein kinase kinase did not affect the induction. FN promoter contained a putative AP-1 binding site (rFN/AP-1; -463 to -437), and the results of a transient transfection and electrophoretic mobility shift assay showed that Ang II enhanced rFN/AP-1 activity. CV11974 and inhibitors of PKC or PTK suppressed Ang II-mediated increases in rFN/AP-1 activity, although neither PD123319 nor a protein kinase A inhibitor affected the induction. Furthermore, mutation of rFN/AP-1 that disrupted nuclear binding suppressed Ang II-induced transcription in the native FN promoter (-1908 to +136) context. Thus, Ang II activates transcription of the FN gene through the Ang II type 1 receptor in vascular smooth muscle cells, at least in part, via the activation of AP-1 by a signaling mechanism dependent on PKC and PTK.

The Maastricht strain and England strain of rat cytomegalovirus represent different betaherpesvirus species rather than strains

The major immediate early (MIE) locus of the Maastricht strain of rat cytomegalovirus (RCMV) was found to comprise five exons of which the first is noncoding. The first three exons are spliced to either exon 4, generating IE1, or exon 5, generating IE2. An additional splicing event unique to RCMV (Maastricht) was identified in exon 5, resulting in a 466-bp deletion. IE1 transcripts were detected exclusively during the IE phase of infection in vitro, whereas IE2 transcripts were detected during both the IE and late phase of infection. The similarities between amino acid sequences derived from the MIE gene of RCMV (Maastricht) and murine cytomegalovirus are low (22 and 37% for IE1 and IE2, respectively). Surprisingly, the similarities between the MIE proteins of RCMV (Maastricht) and the England strain of RCMV are also low (23 and 32% for IE1 and IE2, respectively). This suggests that these RCMV strains represent different betaherpesvirus species rather than strains. This is underscored by the difference between both viruses in genome size as well as growth characteristics. The existence of two different RCMV-like species might have important implications for the use of these viruses as models for human cytomegalovirus.

Transcript heterogeneity of the human reduced folate carrier results from the use of multiple promoters and variable splicing of alternative upstream exons

We previously identified three separate cDNAs (KS6, KS32 and KS43) for the human reduced folate carrier (RFC) with unique 5' untranslated regions (5' UTRs) [Wong, Proefke, Bhushan and Matherly (1995) J. Biol. Chem. 270, 17468-17475]. Multiple RFC transcripts were confirmed in CCRF-CEM cells and transport-up-regulated K562.4CF cells by 5' rapid amplification of cDNA ends (5' RACE) and/or primer extension analysis. Two groups of 5' RACE clones were identified, one containing a variable length sequence identical with the KS43 cDNA 5' UTR, and another consisting of variants of the KS32 5' UTR, apparently generated by alternative splicing. The 5' UTR for the KS6 cDNA was not detected. A single band was detected on Southern blots of CCRF-CEM genomic DNA probed with a 326 bp genomic fragment common to all three cDNA species. The unique 5' UTRs for the KS43 and KS32 transcripts were localized to separate non-coding exons (exons 1 and 2 respectively), upstream from a large (approx. 3.42 kb) intron; the KS6 5'UTR also mapped to exon 1. Exons 1 and 2 were contiguous with 996 and 342 bp GC-rich 5' flanking regions (designated Pro43 and Pro32 respectively) that contained multiple SP1 and AP2 but no TATA or CAAT boxes. Both Pro43 and Pro32 exhibited strong promoter activities when cloned in front of a luciferase reporter gene and transfected into HT1080 and K562 cells. By an analysis of promoter deletion mutants we identified two 89 bp tandem repeats that seemed to increase Pro32 activity, and a 240 bp distal sequence that repressed Pro43 activity. Taken together, our results show that multiple human RFC transcripts are encoded by a single gene locus and that the heterogeneous 5' UTRs result from multiple transcriptional starts and variable splicing of alternative non-coding exons transcribed from separate promoters.

Purification and characterization of transcription factor IIIA from higher plants

Transcription factor IIIA (TF IIIA) binds and specifically activates transcription of eukaryotic 5S rRNA genes. It also forms a 7S ribonucleoprotein complex with mature 5S rRNA. Here, we describe the purification and properties of pTF IIIA from higher plants. The purified protein from tulip (Tulipa whittalii) has a molecular mass of about 40 kDa and also binds 5S rRNA and 5S rRNA genes. pTF IIIA also facilitates the transcription of a 5S rRNA gene in a HeLa cell extract.

A novel proximal element mediates the regulation of mouse Ren-1C promoter by retinoblastoma protein in cultured cells

The protein product of the retinoblastoma susceptibility gene, RB, is a nuclear phosphoprotein that modulates transcription of genes involved in growth control via interactions with transcription factors. Renin is a rate-limiting enzyme of the renin-angiotensin system that regulates blood pressure and water-electrolyte balance. Renin gene expression is regulated in a tissue-specific and developmentally linked manner. Similarly, the expression of RB is controlled in a differentiation-linked manner. Thus, to investigate whether RB is involved in the regulation of renin gene expression, we examined the effects of RB on transcriptional activity of the mouse renin (Ren-1C) promoter. The Ren-1C promoter contains two transcriptionally important elements; the RU-1 (-224 to -138) and RP-2 (-75 to -47) elements. RB activated the Ren-1C promoter in human embryonic kidney cells. The promoter element responsible for RB-mediated transcriptional regulation was the RP-2 element. The results of DNA-protein binding experiments showed that RB increased nuclear binding activity to the RP-2 element, and site-directed mutation which disrupted binding of nuclear factors to the RP-2 element markedly reduced RB-mediated activation of Ren-1C promoter in human embryonic kidney cells. These results indicate that the RP-2 element plays an important role in RB-mediated transcriptional regulation of Ren-1C promoter activity in human embryonic kidney cells, thereby suggesting an interesting mechanism by which RB may modulate the renin-angiotensin system.

A statistical model for locating regulatory regions in genomic DNA

In addition to genes, chromosomal DNA contains sequences that serve as signals for turning on and off gene expression. These signals are thought to be distributed as clusters in the regulatory regions of genes. We develop a Bayesian model that views locating regulatory regions in genomic DNA as a change-point problem, with the beginning of regulatory and non-regulatory regions corresponding to the change points. The model is based on a hidden Markov chain. The data consist of nucleotide positions of protein-binding elements in a genomic DNA sequence. These positions are identified using a reference catalogue containing elements that interact with transcription factors implicated in controlling the expression of protein-encoding genes. Among the protein-binding elements in a genomic DNA sequence, the statistical model automatically selects those that tend to predict regulatory regions. We test the model using viral sequences that include known regulatory regions and provide the results obtained for human genomic DNA corresponding to the beta globin locus on chromosome 11.

TRANSFAC, TRRD and COMPEL: towards a federated database system on transcriptional regulation

Three databases that provide data on transcriptional regulation are described. TRANSFAC is a database on transcription factors and their DNA binding sites. TRRD (Transcription Regulatory Region Database) collects information about complete regulatory regions, their regulation properties and architecture. COMPEL comprises specific information on composite regulatory elements. Here, we describe the present status of these databases and the first steps towards their federation.

The EMBL Nucleotide Sequence Database

The EMBL Nucleotide Sequence Database is a comprehensive database of DNA and RNA sequences directly submitted from researchers and genome sequencing groups and collected from the scientific literature and patent applications. In collaboration with DDBJ and GenBank the database is produced, maintained and distributed at the European Bioinformatics Institute (EBI) and constitutes Europe's primary nucleotide sequence resource. Database releases are produced quarterly and are distributed on CD-ROM. EBI's network services allow access to the most up-to-date data collection via Internet and World Wide Web interface, providing database searching and sequence similarity facilities plus access to a large number of additional databases.