TRANSFAC: an integrated system for gene expression regulation

Insulin and glucocorticoids differentially regulate leptin transcription and secretion in brown adipocytes

Leptin, the ob gene product, is produced by adipose tissue and is submitted to a complex hormonal and metabolic regulation. Leptin plays a critical role in the balance of body weight. Here we report on secretion and hormonal regulation of leptin by brown adipocytes. Using the recently established T37i cell line, we show that leptin expression and secretion occurred as a function of cell differentiation. In differentiated T37i cells, insulin induced leptin release ( approximately 0.25 ng/10(6) cells/h) in a concentration-dependent manner (EC50=0.1 nM), and this was totally suppressed by beta3-adrenergic ligand, thiazolidinedione, cycloheximide, or actinomycin D. Insulin induced a strong, rapid (within 2 h) but transient fivefold increase in leptin mRNA levels. This transcriptional control of ob gene expression by insulin involved both phosphatidylinositol 3-kinase- and MAP kinase-dependent pathways. Glucocorticoids inhibited both insulin-stimulated leptin secretion and ob gene expression without affecting leptin mRNA stability (t(1/2)=3h05). Altogether, our results demonstrate that brown adipocytes express and secrete leptin, whose hormonal regulation clearly differs from that described in white adipose tissue. These findings point to tissue-specific molecular mechanisms and suggest that leptin might exert direct effects on energy homeostasis through an autocrine mechanism.

Cloning and expression pattern of a gene encoding an alpha-xylosidase active against xyloglucan oligosaccharides from Arabidopsis

An alpha-xylosidase active against xyloglucan oligosaccharides was purified from cabbage (Brassica oleracea var. capitata) leaves. Two peptide sequences were obtained from this protein, the N-terminal and an internal one, and these were used to identify an Arabidopsis gene coding for an alpha-xylosidase that we propose to call AtXYL1. It has been mapped to a region of chromosome I between markers at 100.44 and 107.48 cM. AtXYL1 comprised three exons and encoded a peptide that was 915 amino acids long, with a potential signal peptide of 22 amino acids and eight possible N-glycosylation sites. The protein encoded by AtXYL1 showed the signature regions of family 31 glycosyl hydrolases, which comprises not only alpha-xylosidases, but also alpha-glucosidases. The alpha-xylosidase activity is present in apoplastic extractions from Arabidopsis seedlings, as suggested by the deduced signal peptide. The first eight leaves from Arabidopsis plants were harvested to analyze alpha-xylosidase activity and AtXYL1 expression levels. Both increased from older to younger leaves, where xyloglucan turnover is expected to be higher. When this gene was introduced in a suitable expression vector and used to transform Saccharomyces cerevisiae, significantly higher alpha-xylosidase activity was detected in the yeast cells. alpha-Glucosidase activity was also increased in the transformed cells, although to a lesser extent. These results show that AtXYL1 encodes for an apoplastic alpha-xylosidase active against xyloglucan oligosaccharides that probably also has activity against p-nitrophenyl-alpha-D-glucoside.

Computer-assisted identification of cell cycle-related genes: new targets for E2F transcription factors

The processes that take place during development and differentiation are directed through coordinated regulation of expression of a large number of genes. One such gene regulatory network provides cell cycle control in eukaryotic organisms. In this work, we have studied the structural features of the 5' regulatory regions of cell cycle-related genes. We developed a new method for identifying composite substructures (modules) in regulatory regions of genes consisting of a binding site for a key transcription factor and additional contextual motifs: potential targets for other transcription factors that may synergistically regulate gene transcription. Applying this method to cell cycle-related promoters, we created a program for context-specific identification of binding sites for transcription factors of the E2F family which are key regulators of the cell cycle. We found that E2F composite modules are found at a high frequency and in close proximity to the start of transcription in cell cycle-related promoters in comparison with other promoters. Using this information, we then searched for E2F sites in genomic sequences with the goal of identifying new genes which play important roles in controlling cell proliferation, differentiation and apoptosis. Using a chromatin immunoprecipitation assay, we then experimentally verified the binding of E2F in vivo to the promoters predicted by the computer-assisted methods. Our identification of new E2F target genes provides new insight into gene regulatory networks and provides a framework for continued analysis of the role of contextual promoter features in transcriptional regulation. The tools described are available at http://compel.bionet.nsc.ru/FunSite/SiteScan.html.

Transcriptional regulation of the leptin gene promoter in rat GH3 pituitary and C6 glioma cells

Leptin was originally believed to be an exclusively adipocyte-derived hormone regulating appetite and energy balance. It has recently become apparent that leptin is actively expressed in a number of other tissues including the CNS and pituitary, as well as brain- and pituitary-derived cell lines. However, the factors controlling leptin expression in cells of neuroectodermal origin are unknown. The mouse leptin gene 5'-flanking DNA contains multiple AP-1 and SRF-1 binding sites as well as a consensus CRE site at -491 to -482 bp. In addition, a number of potential PIT1 and Oct-1 binding sites may contribute to leptin gene transcription in pituitary and brain. We have used leptin promoter-luciferase reporter constructs to examine the regulation of the leptin promoter in 3T3-L1 preadipocytes, C6 glioma cells, and GH3 pituitary cells in response to serum and hormonal stimuli. Cells were transiently transfected with reporter constructs containing either the proximal 500 bp of the leptin promoter (-500-luc) or 6000 bp of the leptin gene 5' flanking region (-6000-luc). Functional analysis indicates that the leptin promoter is constitutively active in all 3 cell lines. Transcriptional activity was significantly higher with a -500 to +9 promoter than with a construct containing -6000 to +9 bp of 5' flanking DNA, indicating the presence of repressor elements which may contribute to the tissue-specific regulation of leptin expression. However, qualitatively similar results were observed with both constructs in response to serum and hormonal manipulation. Leptin promoter activity was significantly stimulated by serum in all cell lines, although to varying extents. In contrast, the response of the leptin promoter to insulin, IGF-1 and dibutyryl cAMP was cell-type specific and dependent on the presence or absence of FBS in the culture medium. Insulin, IGF-1 and dibutyryl cAMP each caused an approximately two-fold stimulation of leptin promoter activity in 3T3-L1 cells under serum-free conditions, but had no significant effect in the presence of 10% FBS. In contrast, dibutyryl cAMP markedly stimulated leptin promoter activity (5-8-fold) in C6 or GH3 cells in the presence or absence of FBS, whereas insulin or IGF-1 had minimal effects. These findings support our previous studies on the regulation of leptin steady state mRNA levels in C6 cells and demonstrate tissue-specific differences in the regulation of leptin gene transcription in adipose vs. neuroectodermal tissues.

The tumor suppressor protein p53 requires a cofactor to activate transcriptionally the human BAX promoter

An important regulator of the proapoptotic BAX is the tumor suppressor protein p53. Unlike the p21 gene, in which p53-dependent transcriptional activation is mediated by a response element containing two consensus p53 half-sites, it previously was reported that activation of the BAX element by p53 requires additional sequences. Here, it is demonstrated that the minimal BAX response element capable of mediating p53-dependent transcriptional activation consists of two p53 half-sites plus an adjacent 6 base pairs (5'-GGGCGT-3'). This GC-rich region constitutes a "GC box" capable both of binding members of the Sp family of transcription factors, including Sp1 in vitro, and of conferring Sp1-dependent transcriptional activation on a minimal promoter in cells. Mutations within this GC box abrogated the ability of p53 to activate transcription without affecting the affinity of p53 for its binding site, demonstrating that these 6 bases are required for p53-dependent activation. In addition, a positive correlation was observed between the ability of p53 to activate transcription in cells and the ability of Sp1 to bind this response element in vitro. Mutations that inhibited Sp1 binding also blocked the ability of p53 to activate transcription through this element. Together, these results suggest a model in which p53 requires the cooperation of Sp1 or a Sp1-like factor to mediate transcriptional activation of the human BAX promoter.

Integrated genomic and proteomic analyses of a systematically perturbed metabolic network

We demonstrate an integrated approach to build, test, and refine a model of a cellular pathway, in which perturbations to critical pathway components are analyzed using DNA microarrays, quantitative proteomics, and databases of known physical interactions. Using this approach, we identify 997 messenger RNAs responding to 20 systematic perturbations of the yeast galactose-utilization pathway, provide evidence that approximately 15 of 289 detected proteins are regulated posttranscriptionally, and identify explicit physical interactions governing the cellular response to each perturbation. We refine the model through further iterations of perturbation and global measurements, suggesting hypotheses about the regulation of galactose utilization and physical interactions between this and a variety of other metabolic pathways.

Therapeutic modulation of transcription factor activity

Aberrant gene expression is a fundamental cause of many disease-associated pathophysiologies. The pharmacological modulation of transcription factor activity therefore represents an attractive therapeutic approach to such disorders. With the exception of nuclear receptors, which are the direct targets of pharmaceuticals, other known classes of transcription factors are largely regulated indirectly by drugs that impact upon those signal transduction cascades that alter transcription factor phosphorylation and dephosphorylation and/or nuclear import. However, recent advances in drug discovery technologies now enable high-throughput screens that can identify molecules that act directly at the level of transcription factor complexes.

A novel regulatory element in the dnmt1 gene that responds to co-activation by Rb and c-Jun

Rb, c-Jun and dnmt1 play critical roles in the process of cellular differentiation. We demonstrate that a regulatory region of murine dnmt1 contains an element which is responsible for transactivation by Rb and c-Jun in P19 embryocarcinoma cells which is not observed in Y1 adrenocarcinoma cells. During differentiation of P19 cells, the induction of Rb and c-Jun coincides with an increase of dnmt1 mRNA. Using linker scanning mutagenesis we identify the element that is responsible for this activation to be a non-canonical AP-1 site. Our data is an example of how a proto-oncogene activates its downstream effectors by recruiting a tumor suppressor. This interaction of Rb and a proto-oncogene might play an important role in differentiation. The responsiveness of dnmt1 to this type of signal is consistent with an important role for regulated expression of dnmt1 during cellular differentiation.

Expression of CYP6B1 and CYP6B3 cytochrome P450 monooxygenases and furanocoumarin metabolism in different tissues of Papilio polyxenes (Lepidoptera: Papilionidae)

The CYP6B1 and CYP6B3 cytochrome P450 monooxygenases in the midgut of the black swallowtail participate in the metabolism of toxic furanocoumarins present in its host plants. In this study, biochemical analyses indicate that the fat body metabolizes significant amounts of the linear furanocoumarins bergapten and xanthotoxin after larvae feed on xanthotoxin. Northern analyses of the combined CYP6B1/3 transcript expression patterns indicate that transcripts in this P450 subfamily are induced in the midgut and fat body by xanthotoxin. Semi-quantitative RT-PCR analyses of individual CYP6B1/CYP6B3 mRNAs indicate that CYP6B1 transcripts are induced by xanthotoxin in all tissues examined and that CYP6B3 transcripts are induced in the fat body only. These results indicate that the fat body participates in the P450-mediated metabolism of excess furanocoumarins unmetabolized by the midgut. Although transcripts of both genes were detected and CYP6B1 transcripts were induced by xanthotoxin in the integument, furanocoumarin metabolism was not detected. Comparison of these P450 promoters with the promoters of alcohol dehydrogenase genes expressed in the fat bodies of several Drosophila species suggest that the xanthotoxin inducibilities of these P450 genes in fat bodies are regulated by elements other than those modulating expression of Adh genes.

The human ITGB4BP gene is constitutively expressed in vitro, but highly modulated in vivo

The ITGB4BP gene encodes for a highly conserved protein, named p27BBP (also known as eIF6), originally identified in mammals as a cytoplasmic interactor of beta4 integrin. In vitro and in vivo studies demonstrated that p27BBP is essential for cell viability and has a primary function in the biogenesis of the 60S ribosomal subunit. Here we report the genomic organization of the human ITGB4BP gene and show that its gene product is expressed with features of a housekeeping element in vitro, but is regulated in a cell specific fashion in vivo. The human gene spans 10 kb and comprises seven exons and six introns. The 5' flanking region shows a TATA-less promoter, canonical CpG islands, and binding sites for serum responsive elements. In cultured cells, p27BBP mRNA and protein are constitutively expressed and stable. A gradual loss of p27BBP mRNA can be observed only after prolonged serum starvation, and heat shock treatment. In contrast, p27BBP mRNA and protein levels in vivo are variable among different organs. More strikingly, immunohistochemical analysis shows that the p27BBP protein is present in a cell specific fashion, even within the same tissue. Taken together, these data show that ITGB4BP gene expression is highly regulated in vivo, possibly by the combination of tissue specific factors and protein synthesis pathways.

Analysis of PBX1 as a candidate gene for type 2 diabetes mellitus in Pima Indians

The human proto-oncogene PBX1 codes for a homeodomain containing protein that modulates expression of several genes, including those contributing to regulation of insulin action and glucose metabolism. PBX1 is located on chromosome 1q22, a region linked with type 2 diabetes in Pima Indians, Caucasians, and an Old Order Amish population. We have investigated the PBX1 genomic sequence to identify polymorphisms that may contribute to diabetes susceptibility in the Pimas. PBX1 is composed of nine exons spanning approx. 117 kb and is located within 300 kb of microsatellite D1S1677, which marks the peak of linkage to diabetes susceptibility in the Pima Indians. We detected 16 single nucleotide polymorphisms in PBX1 including one causing a glycine to serine substitution at residue 21. Comparison of the frequencies of the polymorphisms between affected and unaffected Pima Indians did not detect any significant differences, indicating that mutations in PBX1 do not explain the linkage of 1q with type 2 diabetes in this population. The genomic structure of PBX1 provides a basis for similar systematic examinations of this candidate locus in other populations in relation to both type 2 diabetes and other metabolic disorders.

Type B leukemogenic virus has a T-cell-specific enhancer that binds AML-1

Type B leukemogenic virus (TBLV) induces rapidly appearing T-cell tumors in mice. TBLV is highly related to mouse mammary tumor virus (MMTV) except that TBLV long terminal repeats (LTRs) have a deletion of negative regulatory elements and a triplication of sequences flanking the deletion. To determine if the LTR triplication represents a viral enhancer element, we inserted the triplication upstream and downstream in either orientation relative to the thymidine kinase promoter linked to the luciferase gene. These experiments showed that upregulation of reporter gene activity by the TBLV triplication was relatively orientation independent, consistent with the activity of eukaryotic enhancer elements. TBLV enhancer activity was observed in T-cell lines but not in fibroblasts, B cells, or mammary cells, suggesting that enhancer function is cell type dependent. To analyze the transcription factor binding sites that are important for TBLV enhancer function, we prepared substitution mutations in a reconstituted C3H MMTV LTR that recapitulates the deletion observed in the TBLV LTR. Transient transfections showed that a single mutation (556M) decreased TBLV enhancer activity at least 20-fold in two different T-cell lines. This mutation greatly diminished AML-1 (recently renamed RUNX1) binding in gel shift assays with a mutant oligonucleotide, whereas AML-1 binding to a wild-type TBLV oligomer was specific, as judged by competition and supershift experiments. The 556 mutation also reduced TBLV enhancer binding of two other protein complexes, called NF-A and NF-B, that did not appear to be related to c-Myb or Ets. AML-1 overexpression in a mammary cell line enhanced expression from the TBLV LTR approximately 30-fold. These data suggest that binding of AML-1 to the TBLV enhancer, likely in combination with other factors, is necessary for optimal enhancer function.

Molecular characterization of the 5' regulatory region of rat sodium-dependent multivitamin transporter gene

Previous studies have demonstrated the involvement of a specialized, Na(+)-dependent carrier-mediated system for biotin uptake in mammalian intestine. The molecular identity of the carrier protein, the Na(+)-dependent multivitamin transporter (SMVT), has recently been identified. Upon characterization of transcript expression in the rat intestine, four distinct transcript variants (I-IV) due to heterogeneity at the 5'-untranslated region were found (Chatterjee NS, Kumar CK, Ortiz A, Rubin SA, and Said HM. Am J Physiol Cell Physiol 277: C605-C613, 1999). This finding raised the possibility that multiple promoters may be involved in driving the transcription of the SMVT gene. To test this possibility, we cloned the 5' regulatory region of the SMVT gene by genome walking. A 6.5-kb genomic DNA fragment was identified and sequenced. Three putative promoters (P1, P2, and P3) that were separated by exons of the four previously identified transcript variants were, indeed, found. P1 was found to contain multiple putative regulatory regions like GATA-1, AP-1, AP-2, and C/EBP, including several repeats of purine-rich regions and two TATA-like elements. P2 and P3 were GC rich and also revealed the presence of many putative regulatory elements including several SP-1 consensus sequences. The functional identity of each promoter and the minimal regions required for its function were established by the luciferase assay following transfection of rat-derived cultured intestinal epithelial IEC-6 cells. The highest functional activity of the cloned promoters was found to be in the order of P1 > P2 > P3. These findings represent the first characterization of the 5' regulatory region of any mammalian SMVT gene and should assist in the understanding of transcriptional regulation of this important gene.

Identification of CCAAT displacement protein (CDP/cut) as a locus-specific repressor of major histocompatibility complex gene expression in human tumor cells

Human major histocompatibility (MHC) class I antigen expression is important in controlling the metastatic growth of malignant tumors. Locus-specific down-regulation of MHC class I gene expression is frequently observed in human tumors, leading to decreased susceptibility to cytotoxic T-cell-mediated lysis. The mechanism of this down-regulation is incompletely understood. Here, we describe the identification of human CCAAT displacement protein (CDP/cut) as a locus-specific repressor of HLA-B and C gene expression. Transient and stable transfections in HeLa and K562 cells demonstrated the presence of a repressor element 650 base pairs upstream of the first exon of HLA-B7. A specific binding complex with the HLA-B7 and Cw2 repressor elements was demonstrated by EMSA. Formation of the EMSA complex was inhibited specifically with polyclonal antiserum to human CDP/cut, demonstrating that CDP/cut binds the HLA-B7 repressor element. The corresponding region of the HLA-A2 promoter neither repressed HLA-A2 gene expression nor bound CDP/cut. Overexpression of CDP/cut in cell lines deficient in CDP/cut resulted in a nearly 4-fold repression of reporter constructs containing the HLA-B7 repressor element but not the corresponding region of the HLA-A2 promoter. Repression of HLA-B and C gene expression by CDP/cut does not involve displacement of NF-Y, nor is CDP/cut associated with the histone deacetylase HDAC1 when bound to the HLA-B7 repressor element. To our knowledge, these results identify CDP/cut as the first example of a locus-specific repressor of MHC class I gene transcription in human tumor cells.

Dissection of the promoter region of the inositol 1,4,5-trisphosphate receptor gene, itr-1, in C. elegans: a molecular basis for cell-specific expression of IP3R isoforms

Inositol 1,4,5-trisphosphate receptors in Caenorhabditis elegans are encoded by a single gene, itr-1. This provides a powerful system in which to dissect the mechanisms that control the tissue-specific expression of molecules that determine the specificity of calcium signalling. We first identified the Caenorhabditis briggsae orthologue of itr-1, Cbitr-1. Comparison of the two itr-1 genes revealed that the chromosomal organisation, gene structure and predicted cDNA and protein sequences were all conserved. The conserved gene structure supports the hypothesis that the itr-1 gene has three promoters, each of which gives rise to an alternative mRNA and hence unique protein. To test this and to identify the roles of the three putative promoters (pA, pB and pC) in regulating itr-1 expression we fused each promoter to the green fluorescent protein gene and identified their expression patterns. Introduction of these transgenes into C. elegans identified unique and defined patterns of green fluorescent protein expression directed by each promoter: pA directs expression in the pharyngeal terminal bulb, the rectal epithelial cells and vulva; pB directs expression in the motor neurone PDA, the amphid socket cells and the spermatheca; pC directs expression in the spermathecal valve, uterine sheath cells, pharyngeal isthmus and intestine. Thus tissue-specific expression of itr-1 variants is directed by three promoters and this results in adjacent cells in the same tissue containing different inositol trisphosphate receptor isoforms. Within pA, four short regions (pA-A to pA-D) of sequence conservation between C. elegans and C. briggsae were identified. Deletion analysis demonstrated that the region containing pA-C is required for expression in the terminal bulb and rectal epithelial cells and the region containing pA-D is required for expression in the vulva. pA-C includes sequences similar to the binding sites for transcription factors that have been demonstrated to be important in pharyngeal development and gene expression.

Nerve growth factor, but not epidermal growth factor, increases Fra-2 expression and alters Fra-2/JunD binding to AP-1 and CREB binding elements in pheochromocytoma (PC12) cells

In pheochromocytoma (PC12) cells nerve growth factor (NGF) and epidermal growth factor (EGF) activate similar receptor tyrosine kinase signaling pathways but evoke strikingly different biological outcomes: NGF induces differentiation and EGF acts as a mitogen. A novel approach was developed for identifying transcription factor activities associated with NGF-activated, but not EGF-activated, signaling, using random oligonucleotide clones from a DNA recognition library to isolate specific DNA binding proteins from PC12 nuclear extracts. A protein complex from NGF-treated, but not EGF-treated, cells was identified that exhibits increased mobility and DNA binding activity in gel mobility shift assays. The binding complex was identified in supershift assays as Fra-2/JunD. The clones used as probes contain either AP-1 or cAMP response element binding (CREB) recognition elements. Time course experiments revealed further differences in NGF and EGF signaling in PC12 cells. NGF elicits a more delayed and sustained ERK phosphorylation than EGF, consistent with previous reports. Both growth factors transiently induce c-fos, but NGF evokes a greater response than EGF. NGF specifically increases Fra-1 and Fra-2 levels at 4 and 24 hr. The latter is represented in Western blots by bands in the 40-46 kDa range. NGF, but not EGF, enhances the upper bands, corresponding to phosphorylated Fra-2. These findings suggest that prolonged alterations in Fra-2 and subsequent increases in Fra-2/JunD binding to AP-1 and CREB response elements common among many gene promoters could serve to trigger broadly an NGF-specific program of gene expression.

Centromeric tandem repeat from the chaffinch genome: Isolation and molecular characterization

A new family of avian centromeric satellites is described. The highly repeated sequence, designated FCP (Fringilla coelebs PstI element), was cloned from the 500-bp PstI digest fraction of the chaffinch (Fringilla coelebs L.) genomic DNA, sequenced, and characterized. The FCP repeat was found to have 505–506 bp length of monomer, 57% content of GC, to compose about 0.9% of the chaffinch genome, and to be highly methylated. Results of Southern-blot hybridization of cloned FCP element onto genomic DNA digested with different restriction enzymes, and sequencing directly from total genomic DNA using FCP-specific primers and ThermoFidelase enzyme (Fidelity Systems Inc.) were in agreement with a tandem arrangement of this repeat in the chaffinch genome. Five positions of single-nucleotide polymorphism (SNP) were found in the FCP monomers using direct genomic sequencing. Fluorescence in situ hybridization (FISH) with FCP probe and primed in situ labelling (PRINS) with FCP specific primers showed that the FCP elements occupy pericentric regions of all chaffinch chromosomes. On chromosome spreads, the fluorescent signals were also observed in the intercentromeric connectives between nonhomologous chromosomes. The results suggest that the centromeric FCP repeat is responsible for chromosome ordering during mitosis in chaffinch.Key words: satellite DNA, centromeric heterochromatin, Fringilla coelebs.

Genomic strategies to identify mammalian regulatory sequences

With the continuing accomplishments of the human genome project, high-throughput strategies to identify DNA sequences that are important in mammalian gene regulation are becoming increasingly feasible. In contrast to the historic, labour-intensive, wet-laboratory methods for identifying regulatory sequences, many modern approaches are heavily focused on the computational analysis of large genomic data sets. Data from inter-species genomic sequence comparisons and genome-wide expression profiling, integrated with various computational tools, are poised to contribute to the decoding of genomic sequence and to the identification of those sequences that orchestrate gene regulation. In this review, we highlight several genomic approaches that are being used to identify regulatory sequences in mammalian genomes.

Identification of a cis-acting element and a novel trans-acting factor of the human insulin receptor gene in HepG2 and rat liver cells

The liver is a major target organ of insulin and is important for glucose homeostasis. We analyzed the tissue specific regulation of the insulin receptor gene in the liver by studying the cis-acting element and trans-acting factor of the human insulin receptor gene in human hepatoma cell line, HepG2 cells. In the chloramphenicol acetyl transferase (CAT) assay with chimeric plasmids containing various deletions and insertions of the human insulin receptor promoter/CAT gene, a HepG2 cell specific cis-acting element was identified between nt -592 to -577 of the promoter. In electrophoretic mobility shift assay and UV cross-link analysis, a 35-kDa nuclear protein that bound to 5'-TCCCTCCC-3' (nt -588 to -581) sequence was identified in HepG2 cells as well as in rat hepatocytes. This nuclear protein, designated as hepatocyte-specific transcription factor of the insulin receptor gene (HTFIR), might play an important role in tissue-specific expression of the insulin receptor gene in the liver.

Assessing clusters and motifs from gene expression data

Large-scale gene expression studies and genomic sequencing projects are providing vast amounts of information that can be used to identify or predict cellular regulatory processes. Genes can be clustered on the basis of the similarity of their expression profiles or function and these clusters are likely to contain genes that are regulated by the same transcription factors. Searches for cis-regulatory elements can then be undertaken in the noncoding regions of the clustered genes. However, it is necessary to assess the efficiency of both the gene clustering and the postulated regulatory motifs, as there are many difficulties associated with clustering and determining the functional relevance of matches to sequence motifs. We have developed a method to assess the potential functional significance of clusters and motifs based on the probability of finding a certain number of matches to a motif in all of the gene clusters. To avoid problems with threshold scores for a match, the top matches to a motif are taken in several sample sizes. Genes from a sample are then counted by the cluster in which they appear. The probability of observing these counts by chance is calculated using the hypergeometric distribution. Because of the multiple sample sizes, strong and weak matching motifs can be detected and refined and significant matches to motifs across cluster boundaries are observed as all clusters are considered. By applying this method to many motifs and to a cluster set of yeast genes, we detected a similarity between Swi Five Factor and forkhead proteins and suggest that the currently unidentified Swi Five Factor is one of the yeast forkhead proteins.

The EMBL nucleotide sequence database

The EMBL Nucleotide Sequence Database (http://www.ebi.ac.uk/embl/) is maintained at the European Bioinformatics Institute (EBI) in an international collaboration with the DNA Data Bank of Japan (DDBJ) and GenBank at the NCBI (USA). Data is exchanged amongst the collaborating databases on a daily basis. The major contributors to the EMBL database are individual authors and genome project groups. Webin is the preferred web-based submission system for individual submitters, whilst automatic procedures allow incorporation of sequence data from large-scale genome sequencing centres and from the European Patent Office (EPO). Database releases are produced quarterly. Network services allow free access to the most up-to-date data collection via ftp, email and World Wide Web interfaces. EBI's Sequence Retrieval System (SRS), a network browser for databanks in molecular biology, integrates and links the main nucleotide and protein databases plus many specialized databases. For sequence similarity searching a variety of tools (e.g. Blitz, Fasta, BLAST) are available which allow external users to compare their own sequences against the latest data in the EMBL Nucleotide Sequence Database and SWISS-PROT.

rSNP_Guide, a database system for analysis of transcription factor binding to target sequences: application to SNPs and site-directed mutations

rSNP_Guide is a novel curated database system for analysis of transcription factor (TF) binding to target sequences in regulatory gene regions altered by mutations. It accumulates experimental data on naturally occurring site variants in regulatory gene regions and site-directed mutations. This database system also contains the web tools for SNP analysis, i.e., active applet applying weight matrices to predict the regulatory site candidates altered by a mutation. The current version of the rSNP_Guide is supplemented by six sub-databases: (i) rSNP_DB, on DNA-protein interaction caused by mutation; (ii) SYSTEM, on experimental systems; (iii) rSNP_BIB, on citations to original publications; (iv) SAMPLES, on experimentally identified sequences of known regulatory sites; (v) MATRIX, on weight matrices of known TF sites; (vi) rSNP_Report, on characteristic examples of successful rSNP_Tools implementation. These databases are useful for the analysis of natural SNPs and site-directed mutations. The databases are available through the Web, http://wwwmgs.bionet.nsc.ru/mgs/systems/rsnp/.

Genomic organization of the human phosphodiesterase PDE11A gene. Evolutionary relatedness with other PDEs containing GAF domains

PDE11A is a dual-substrate, cAMP and cGMP, cyclic nucleotide phosphodiesterase (PDE). Presently four unique variants carrying distinct GAF sequences in the N-terminal region have been identified. While human PDE11A3 and PDE11A4 are known to be specifically expressed in testis and prostate, respectively, PDE11A1 was mainly detected in skeletal muscle. The human PDE11A gene was investigated and revealed to span > 300 kb, contain 23 exons and be mapped on chromosome 2q31. The transcription start sites of PDE11A1, PDE11A3 and PDE11A4 were determined, and the promoter sequences were revealed. Although 5' flanking genomic regions of PDE11A1 and PDE11A3 had a consensus TATA motif, that of PDE11A4 was a TATA-less but contained CCAAT box and Sp1-binding sequence. Interestingly, we found that the exon 2 sequence for N-terminal region of PDE11A3 encoded an N-terminal sequence of the cytochrome c pseudogene in an alternate reading frame, and that C-terminal region of the cytochrome c pseudogene in intron 2 was disrupted by the insertion of Alu repetitive sequence. Furthermore, we examined the exon-intron organization of the PDE2A gene and compared the exon organization among GAF-PDE family. The exon organization of the PDE11A catalytic domain was very similar to those of PDE5A and PDE6B. However, other GAF-PDEs, PDE2A and PDE10A, displayed different exon organization from PDE11A although these three PDEs are similar in their amino-acid sequences to each other. The findings suggested that PDE11A has a common ancestral gene with PDE5A and PDE6s, whereas PDE2A and PDE10A are generated separately from these three GAF-PDEs.

ACTIVITY: a database on DNA/RNA sites activity adapted to apply sequence-activity relationships from one system to another

ACTIVITY is a database on DNA/RNA site sequences with known activity magnitudes, measurement systems, sequence-activity relationships under fixed experimental conditions and procedures to adapt these relationships from one measurement system to another. This database deposits information on DNA/RNA affinities to proteins and cell nuclear extracts, cutting efficiencies, gene transcription activity, mRNA translation efficiencies, mutability and other biological activities of natural sites occurring within promoters, mRNA leaders, and other regulatory regions in pro- and eukaryotic genomes, their mutant forms and synthetic analogues. Since activity magnitudes are heavily system-dependent, the current version of ACTIVITY is supplemented by three novel sub-databases: (i) SYSTEM, measurement systems; (ii) KNOWLEDGE, sequence-activity relationships under fixed experimental conditions; and (iii) CROSS_TEST, procedures adapting a relationship from one measurement system to another. These databases are useful in molecular biology, pharmacogenetics, metabolic engineering, drug design and biotechnology. The databases can be queried using SRS and are available through the Web, http://wwwmgs. bionet.nsc.ru/systems/Activity/.

Mendel-GFDb and Mendel-ESTS: databases of plant gene families and ESTs annotated with gene family numbers and gene family names

There is no control over the information provided with sequences when they are deposited in the sequence databases. Consequently mistakes can seed the incorrect annotation of other sequences. Grouping genes into families and applying controlled annotation overcomes the problems of incorrect annotation associated with individual sequences. Two databases (http://www.mendel.ac.uk) were created to apply controlled annotation to plant genes and plant ESTs: Mendel-GFDb is a database of plant protein (gene) families based on gapped-BLAST analysis of all sequences in the SWISS-PROT family of databases. Sequences are aligned (ClustalW) and identical and similar residues shaded. The families are visually curated to ensure that one or more criteria, for example overall relatedness and/or domain similarity relate all sequences within a family. Sequence families are assigned a 'Gene Family Number' and a unified description is developed which best describes the family and its members. If authority exists the gene family is assigned a 'Gene Family Name'. This information is placed in Mendel-GFDb. Mendel-ESTS is primarily a database of plant ESTs, which have been compared to Mendel-GFDb, completely sequenced genomes and domain databases. This approach associated ESTs with individual sequences and the controlled annotation of gene families and protein domains; the information being placed in Mendel-ESTS. The controlled annotation applied to genes and ESTs provides a basis from which a plant transcription database can be developed.

Nerve growth factor, but not epidermal growth factor, increases Fra-2 expression and alters Fra-2/JunD binding to AP-1 and CREB binding elements in pheochromocytoma (PC12) cells

In pheochromocytoma (PC12) cells nerve growth factor (NGF) and epidermal growth factor (EGF) activate similar receptor tyrosine kinase signaling pathways but evoke strikingly different biological outcomes: NGF induces differentiation and EGF acts as a mitogen. A novel approach was developed for identifying transcription factor activities associated with NGF-activated, but not EGF-activated, signaling, using random oligonucleotide clones from a DNA recognition library to isolate specific DNA binding proteins from PC12 nuclear extracts. A protein complex from NGF-treated, but not EGF-treated, cells was identified that exhibits increased mobility and DNA binding activity in gel mobility shift assays. The binding complex was identified in supershift assays as Fra-2/JunD. The clones used as probes contain either AP-1 or cAMP response element binding (CREB) recognition elements. Time course experiments revealed further differences in NGF and EGF signaling in PC12 cells. NGF elicits a more delayed and sustained ERK phosphorylation than EGF, consistent with previous reports. Both growth factors transiently induce c-fos, but NGF evokes a greater response than EGF. NGF specifically increases Fra-1 and Fra-2 levels at 4 and 24 hr. The latter is represented in Western blots by bands in the 40-46 kDa range. NGF, but not EGF, enhances the upper bands, corresponding to phosphorylated Fra-2. These findings suggest that prolonged alterations in Fra-2 and subsequent increases in Fra-2/JunD binding to AP-1 and CREB response elements common among many gene promoters could serve to trigger broadly an NGF-specific program of gene expression.

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.

Isolation of a bi-directional promoter directing expression of the mouse GABPalpha and ATP synthase coupling factor 6 genes

The GA-binding protein (GABP) is a ubiquitous heteromeric transcription factor implicated in the regulation of several genes involved in mitochondrial energy metabolism including subunits of cytochrome c oxidase, ATP synthase, and mitochondrial transcription factor 1 (mtTF1). GABPalpha subunit binds the PEA3/Ets binding sites (EBS), while GABPbeta contains a transcription activation domain and mediates alphabeta dimer and alpha(2)beta(2) tetramer formation essential for activation of transcription. Here we report the cloning of 2449 bp of the mouse (m) GABPalpha promoter region including 201 bp of the 5' end of the published mGABPalpha cDNA sequence. Surprisingly, sequences homologous to the 5'UTR of mouse, rat and human mitochondrial ATP synthase coupling factor 6 (ATPsynCF6) cDNAs were found165-240 bp upstream of the mGABPalpha cDNA. A search of the non-redundant nucleotide database revealed a human genomic sequence derived from chromosome 21 (21q22) bearing significant homology to the mGABPalpha/ATPsynCF6 promoter region and encompassed the entire hGABPalpha and hATPsynCF6 genes. Primer extension analysis revealed multiple transcription start sites for both mGABPalpha and mATPsynCF6 mRNAs that mapped near the published cDNA 5' ends. Sequence analysis identified several binding sites upstream of the GABPalpha cDNA sequence including sites for GABP (-86, -104, -169, -257, and -994), YY1 (-57), Sp1 (-242 and -226), and NRF1 (-5). No 'TATA' motif was identified near either the GABPalpha or ATPsynCF6 transcription start sites. The human and mouse promoters retain significant sequence identity including binding sites for several tissue-specific transcription factors. Transient transfection assays using Luciferase reporter constructs containing the intergenic region and flanking sequences confirmed that this region of DNA promotes transcription in both directions.

Identification and characterization of two CD40-inducible enhancers in the mouse TRAF1 gene locus

We have shown that CD40 engagement induces TRAF1 gene expression in B lymphocytes. Here we report that CD40-dependent TRAF1 gene transcription in murine B cells is controlled by two enhancer regions. One region is located approximately 2 kb upstream of the transcription start site and the other lies in the intron between exons 5 and 6. The upstream enhancer contains a single NF-kappaB site in addition to sites that bind constitutive transcription factors. Mutation of this NF-kappaB site completely abrogates CD40-driven TRAFl transcription. The intronic enhancer contains two sites that strongly bind the CD40-inducible factors NF-kappaB and AP-1. Simultaneous mutation of the AP-1 site and of the NF-kappaB site abolishes transcription driven by this enhancer. When cloned together into reporter constructs, the two TRAF1 enhancers do not synergize, suggesting that each enhancer may separately participate in the induction of TRAF1 transcription in B cells following CD40 activation.

The human phosphodiesterase PDE10A gene genomic organization and evolutionary relatedness with other PDEs containing GAF domains

PDE10A is a cyclic nucleotide phosphodiesterase (PDE) exhibiting properties of a cAMP PDE and a cAMP-inhibited cGMP PDE. The transcripts are specifically expressed in the striatum. The human gene encoding PDE10A was cloned and investigated. The PDE10A gene spanned > 200 kb and contained 24 exons. The exon-intron organization of PDE10A was different from those of PDE5A and PDE6B, although these three PDEs include two GAF domains and have similar amino-acid sequences. The promoter sequence of PDE10A was highly GC-rich and did not contain a TATA motif and a CAAT box, suggesting it is a housekeeping gene. In Caenorhabditis elegans, the C32E12.2 gene encoding a probable PDE that is 48% identical to the human PDE10A protein showed similar exon organization to PDE10A but not PDE5A and PDE6B. This, together with the phylogenic tree analysis, suggested that the ancestral gene for PDE10A existed in a lower organism such as C. elegans.

An intron transcriptional enhancer element regulates IL-4 gene locus accessibility in mast cells

The cell type-specific expression of a gene is dependent on developmentally regulated modifications in chromatin structure that allow accessibility of basal and inducible transcription factors. In this study, we demonstrate that a cis-acting element in the second intron of the murine IL-4 gene has a dual function in regulating transcription in mast cells as well as chromatin accessibility of the IL-4 gene locus through its influence on the methylation state of the gene. Previous studies have shown that mast cell-restricted transcription factors GATA-1/2 and PU.1 associate with the intron element and regulate its activity. In this study, we use DNase I footprinting and mutational analyses to identify two additional sites that contribute to the element's ability to enhance transcription. One of these sites associates preferentially with STAT5a and STAT5b. We also demonstrate that deletion of the element or mutation of the GATA binding site in the context of a stably integrated IL-4 genomic construct prevents maintenance of a demethylated locus in IL-4-producing mast cells. These data indicate that, analogous to Ig and TCR intron regulatory elements, the intron enhancer has an essential role in maintaining developmentally regulated demethylation at the IL-4 gene locus. In addition, they indicate that members of the GATA family of transcription factors likely play an important role in these processes.

Representing and analysing molecular and cellular function using the computer

Determining the biological function of a myriad of genes, and understanding how they interact to yield a living cell, is the major challenge of the post genome-sequencing era. The complexity of biological systems is such that this cannot be envisaged without the help of powerful computer systems capable of representing and analysing the intricate networks of physical and functional interactions between the different cellular components. In this review we try to provide the reader with an appreciation of where we stand in this regard. We discuss some of the inherent problems in describing the different facets of biological function, give an overview of how information on function is currently represented in the major biological databases, and describe different systems for organising and categorising the functions of gene products. In a second part, we present a new general data model, currently under development, which describes information on molecular function and cellular processes in a rigorous manner. The model is capable of representing a large variety of biochemical processes, including metabolic pathways, regulation of gene expression and signal transduction. It also incorporates taxonomies for categorising molecular entities, interactions and processes, and it offers means of viewing the information at different levels of resolution, and dealing with incomplete knowledge. The data model has been implemented in the database on protein function and cellular processes 'aMAZE' (http://www.ebi.ac.uk/research/pfbp/), which presently covers metabolic pathways and their regulation. Several tools for querying, displaying, and performing analyses on such pathways are briefly described in order to illustrate the practical applications enabled by the model.

Genomic organization and promoter analysis of the murine Id3 gene

Overexpression of the helix-loop-helix motif-containing transcription inhibitor Id3 has been shown to repress muscle-specific gene expression. Consistent with its putative negative regulatory role in the myogenic process, Id3 is highly expressed in proliferating myoblasts but down regulated when myoblasts are induced to differentiate. To investigate how Id3 expression may be transcriptionally regulated, we isolated a mouse Id3 genomic DNA fragment and characterized its organization and promoter activity. Comparison of the Id3 gene from human and mouse demonstrated a conserved exon-intron organization in which the first intron interrupts the C-terminal protein coding region and the second intron interrupts the 3' untranslated region at analogous positions in the two species. Sequence analysis of the 5'-flanking region revealed an unexpected mouse strain-specific genetic polymorphism due to a single base substitution. Deletion analysis revealed that as little as 180 base pairs of the mouse Id3 promoter upstream of the transcription start site is sufficient for a high level of gene expression in proliferating C2C12 myoblasts. In particular, the region between the nucleotide position -180 and -34 appeared to be crucial for maximal reporter gene activity and interacted specifically with C2C12 nuclear proteins. Finally, we showed that, despite the creation of a putative transcription factor-binding site, the genetic polymorphism observed did not affect Id3 promoter activity in proliferating C2C12 cells.

Alfresco--a workbench for comparative genomic sequence analysis

Comparative analysis of genomic sequences provides a powerful tool for identifying regions of potential biologic function; by comparing corresponding regions of genomes from suitable species, protein coding or regulatory regions can be identified by their homology. This requires the use of several specific types of computational analysis tools. Many programs exist for these types of analysis; not many exist for overall view/control of the results, which is necessary for large-scale genomic sequence analysis. Using Java, we have developed a new visualization tool that allows effective comparative genome sequence analysis. The program handles a pair of sequences from putatively homologous regions in different species. Results from various different existing external analysis programs, such as database searching, gene prediction, repeat masking, and alignment programs, are visualized and used to find corresponding functional sequence domains in the two sequences. The user interacts with the program through a graphic display of the genome regions, in which an independently scrollable and zoomable symbolic representation of the sequences is shown. As an example, the analysis of two unannotated orthologous genomic sequences from human and mouse containing parts of the UTY locus is presented.