- Research article
- Open Access
Osteo-Promoter Database (OPD) – Promoter analysis in skeletal cells
BMC Genomicsvolume 6, Article number: 46 (2005)
Increasing our knowledge about the complex expression of genes in skeletal tissue will provide a better understanding of the physiology of skeletal cells. The study summarizes transcriptional regulation factors interacting and cooperating at promoter regions that regulate gene expression. Specifically, we analyzed A/T rich elements along the promoter sequences.
The Osteo-Promoter Database (OPD) is a collection of genes and promoters expressed in skeletal cells. We have compiled a new viewer, OPD, as unique database developed and created as an accessible tool for skeletal promoter sequences. OPD can navigate to identify genes specific to skeletal cDNA databases and promoter analysis sites. OPD offers exclusive access to facilitate a dynamic extraction of promoters' gene-specific analyses in skeletal tissue. The data on promoters included in OPD contains cloned promoters or predicted promoters that were analyzed by bioinformatics tools. OPD offers MAR-analysis, which allocates A/T rich elements along these promoter sequences.
The analysis leads to a better insight of proteins that bind to DNA, regulate DNA, and function in chromatin remodeling. The OPD is a distinctive tool for understanding the complex function of chromatin remodeling and transcriptional regulation of specific gene expression in skeletal tissue.
The Human Genome Project has provided large-scale sequencing and a multitude of genetic databases, based on information from bacteria, plant, Drosophila, vertebrates and mammals [1–5]. These databases represent evolutionary information that is summarized in the public databases NCBI, EMBL, and HUGO. In silico biology explores various tools that enable the translation of raw data into workable models and provides guidance for high-throughout gene analysis to 'make sense' of the genomic data. Computational tools analyzed the genetic information and data mining aided the knowledge interpretation in the functional cellular context. Genes expression from tissue-specific libraries can be developed to new models and add insight into biological information in a context of gene networks, protein pathways ect. Genes in eukaryotic cells are regulated in "active regions," in which chromatin structure is "open" and accessible to DNA-binding proteins and "silent regions" where "packed" chromatin renders the DNA inaccessible. The regulation is performed on several epigenetic levels that include DNA methylation, nucleosome positioning, histone modification, and additional components of higher-order structure. These modifications can activate or repress RNA synthesis that is associated with tissue-specific genes in the absence of transcriptional enhancers. The details of how chromatin structures modified are unresolved, but it is clear that DNA methylation has a direct effect on both histone acetylation and on higher-order chromatin structure that regulates DNA-chromatin structure .
Transcriptional regulation factors interacting and cooperating at promoter regions determine the fate of cells and their function. The promoter is a cis-acting element immediately upstream of the transcription start site (TSS) that controls the rate of the initiation of transcription. The promoter region is important in regulation of various processes in development, morphogenesis, cell differentiation; hormonal communication and cellular stress responses. To date, only a small fraction of promoters is identified and most are not yet cloned. To better understand tissue specificity, there is a need to explore the interactions that occur between chromatin remodeling proteins and the promoters. Such interactions take place next to transcription-factor (TF) binding sites. An important motif presented in the promoter is an A/T-rich region that binds proteins with an A/T-hook motif. Numerous chromatin-remodeling proteins identified by presence of single or multiple A/T hook motifs that form complex of the protein with the DNA . The A/T rich regions are characteristic by MAR frequently observed at the genomic regions and affect the interactions with the A/T hook motif of the proteins .
The knowledge evolve transcription factor binding sites on the promoter will allow the creation of the regulation map of each gene and will enable the recognition of other transcription regulation systems that control gene expression. To learn how promoters operate, it is first necessary to identify them. To date, a relatively small number of promoters have been identified, and, thus, it is possible to predict the location of the promoter upstream to the coding regions of the gene sequence. Gene-structure prediction programs, such as BCM Gene Finder tool, ConPro , Core Promoter , FirstEF , Gene2Promoter , FunSiteP  and MatInspector  enable to classify promoters. Promoters are being analyzed on species levels as well, for example, the Saccharomyces promoter database (SCPD) [14, 15] and C. elegans promoter database (CEPDB) http://cgsigma.cshl.org/CEPDB/, EPD is a Eukaryotic Promoter Database , and LSPD is a liver-specific tissue promoter database http://rulai.cshl.org/software/index1.htm. Transcription factors bind to regulatory elements upstream of transcription start sites and interact with other factors to create an abundant regulatory environment. Each regulatory protein binds to a specific element on the promoter, in response to an appropriate signal. The complexity of the regulatory system can be explored by the use of software that searches for specific transcription elements on the promoter. Genomatix offers S/MARt DB  and S/MARs , aid in the analysis of scaffold/matrix-attached regions and the nuclear matrix proteins.
Construction and content
The OPD homepage contains a list of genes, where each gene has one entry link that enables the user to obtain the specific gene description. Skeletal genes which are cataloged in OPD are obtained from Human Bone libraries such as: Skeletal Gene Database (SGD) http://sgd.nia.nih.gov/ and Human Bone Marrow stroma (lib. 931) http://www.ncbi.nlm.nih.gov/Unigene/library.cgi?ORG=Hs&LID=931/. Bone tissue specific expression was defined using the Unigene database http://www.ncbi.nlm.nihgov/entrez/query.fcgi?CMD=search&DB=unigene. Information about gene products, their cellular function and chromosomal locus imported from the NCBI database. Accession numbers were used to link the gene entries to information in the resource databases. Accession numbers were obtained from NCBI and HUGO. Information about human genetic disorders and their associated genes is derived from the Online Mendelian Inheritance in Man database, OMIM http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=OMIM. All the references were structured with hyperlinked both to the original Web servers and to the NCBI and HUGO database and OMIM. Searching for a gene's promoter was performed using NCBI and PUBMED databases. The accession number of each promoter is linked to the promoter sequence page in the NCBI database. MAR-Analysis in the promoter region is carried out using the MAR-FINDER program http://www.futuresoft.org/modules/MarFinder/. The MAR Finder program algorithm scores A/T nucleotides along promoters using a 200 bp window at 10 bp intervals with MAR potential of 0.1 as minimum significance threshold.
Utility; Osteopromoter Database (OPD) description
The Osteo-Promoter Database (OPD) summarizes the information about genes and promoters were designed based on analysis of resources covering skeletal cells. The functional genes presented in OPD were annotated and catalogued from human bone marrow stroma  and the skeletal gene database (SGD) . The OPD is a novel bone-tissue promoter database that includes hundreds of genes that mediate osteoblast cells proliferation and differentiation [19, 20, 22–25]. The genes and promoters included in OPD from the PUBMED, bibliographic references databases, and cDNA sequences centers and more promoters that were identified using Gene2Promoter. The data is organized in alphabetical order and have cross-references with other databases.
In OPD, we summarized series of genes known to play a role in the regulation of skeletal cells. The information explores the transcription controlled by various factors on skeletal cells [22–25]. OPD presents distribution of transcription-factor sites in promoter regions and the connection with chromatin structure [26, 27]. Promoters were analyzed and enlighten on the connections between transcription factors and chromatin remodeling factors. We choose to focus on A/T-rich region that binds A-T hook proteins and serves as interface between transcription factors and chromatin remodeling proteins. MAR regions that contain A/T-rich elements are embedded along the genome and recognized for their role in the regulation of gene transcription by creating an environment that controls gene expression at the chromatin remodeling level.
Specifically, examples are given to genes related to progenitors proliferate and differentiate to osteoblasts that form the calcified matrix. Table 1 summarizes genes related to cell proliferation, differentiation and matrix deposition in the formed bone. Transcription factors and hormones control the expression of genes responsible for cell proliferation. The cell's differentiation requires sequential activation of regulatory proteins and signaling molecules such as BMPs and hormones. The osteoblast maturation is defined by the biosynthesis and organization of extracellular matrix proteins, such as osteocalcin, bone sialoprotein and osteonectin . We identified from the literature two hundred promoters' that were cloned, and the rest of the promoters included in OPD, were not cloned but identified using bioinformatics tools.
OPD catalogs hundreds of genes and promoters that are organized in alphabetical order, with direct access from the OPD homepage to each specific gene possible via the list. To view a specific gene, one can view the list and choose the gene or alternatively, identify the gene through a search engine. To search for the gene, the user types the full name of the gene/protein into the search engine box and clicks on the "search" button (Fig. 1). Entry into the gene page will provide a description for each gene, including expression analysis and accession numbers structured as specific links. The information for each gene/protein on the promoter includes chromosomal location, catalog international number of sequences, cellular function, mutations (OMIM), and involvement in differentiation or diseases (Fig. 2, Table 2). A/T-rich region analyses are available through the "see picture" link next to the name of the gene. Absence of a link means that promoter of the gene has not been investigated yet. On the left side of the page, there are useful links to the promoter analysis and to the genomic databases sites (Fig. 1). The OPD presents promoter sequences identified and readable from positions within sequences of the HUGO and NCBI databases. Other links to skeletal cDNA databases and promoter analysis sites are presented.
The database is constructed based on Access database and is presented on the Web pages were created dynamically by ASP scripts written in html and VB script language. The A/T-rich graphs designed in Bitmap format.
OPD summarizes genes that are important in the control of osteoblast differentiation. OPD includes hundreds of promoters that were identified by cloning and the study of the promoter function was performed at the molecular level. Other promoters included were not cloned, but are related to genes relevant to skeletal-cell differentiation and function, which we investigated based on each promoter's prediction from the genomic region. The task of OPD focused on resolving A/T-rich sequences in a defined sequence, using the MAR-FINDER algorithm. The analysis by MAR-FINDER is based on the integration of the frequency of the motif in a defined sequence in relation to the probability of their random occurrence in the same sequence. This software scores A/T nucleotides in a given sequence, and statistical significance is associated with the frequency of occurrence of the pattern influenced by A-T content.
Using bioinformatics we predicted promoter structures from genomic contigs that had not been cloned and followed with the analyzed of A/T-rich elements. Herein we present few examples presented at the OPD; TGFβ3 and TGFRβ3 promoters were predicated. TGFβ3 is constructed from 7 exons and the promoter region upstream to TSS in contig NT_026437 identified A/T rich at -1180-1970 bp. TGFRβ3 gene was identified in contig NT_034388 includes 18 exons and the promoter was identified from the same genomic region, includes three A/T-rich elements at -420-1120 bp, -1250-1330 bp, and 1690-1970 bp. An analysis for BMP factors that are part of the TGF-β superfamily know for their role in regulating the growth and differentiation of osteoblasts. Promoters for BMP2, BMP4, and BMP7 were analyzed and identified the A/T-rich element for the BMP2 at -250-820 bp, and for BMP4 promoter three A/T-rich elements were identified at 250–430 bp, 670–930 bp, and 1170–2150 bp, but none were identified for BMP7. The BMP4 promoter includes three A/T-rich elements with high MARs potentials and thus suggests a significant potential to bind protein with an A-T hook motif.
The differences among MAR-potential are reflected by the percentage of A/T-rich sequences, higher MAR regions suggests a higher potential for binding of a protein with an A-T hook domain. MARs play a role in the context of long-distance interactions mediated between locus control regions (LCR) and the target promoters by creating loops of DNA. The DNA loop keeps the flexibility in A/T-rich length for its functionality. The presence of several copies of A/T-rich sites and longer stretches of A/T-rich sequence create a wider surface area for bonding a specific protein [28, 29]. HMG-I/Y protein contains three A-T hooks responsible for multivalent interactions between the A-T hooks and appropriately spaced A/T-rich tracts are required for high-affinity binding. A single A-T hook in a polypeptide or a single A-T tract on the target DNA results in reduced interaction . Furthermore, the C-terminal acidic region and the spacer regions between A-T hooks also affect the protein – DNA binding .
Fibroblast growth factors (FGFs) are regulators of cells division and differentiation  also control alkaline phosphatase; and regulate matrix proteins synthesis such as collagen, osteonectin, and osteocalcin. Four FGF receptors (FGFR 1–4) members of the tyrosine kinase receptors family were identified. FGFR1 and FGFR2 promoters were predicted and analyzed for MAR region results with two A/T-rich sites with high MAR potential for FGFR1 and lower for FGFR2. These different MAR potentials expressed in the promoters of functional genes assumed for the potential of protein to bind at a different affinity in these promoters. The analysis of promoters related to functional genes in osteogenic differentiation results with annotated genes extended for analysis by MAR-FINDER regions (Table 2). The results have shown that 73% of these genes expressing MARs regions suggesting that MAR element are a common regulatory element in transcriptional control. Thus, it indicates for A/T-rich regions a role in control and guidance of activation and repression of osteoblastic differentiation. The accessibility of proteins to A/T-rich sites on promoters of tissue-specific genes along with other factors will promote the comprehensive specificity of the regulatory process of functional genes in bone tissue.
The presence of several A/T-rich sites on promoters and their high MAR potential implies for their involvement in regulation by creating the boundaries for protein interactions. Along the human genome, 100,000 MAR regions contain A/T-rich elements. The MAR regions were recognized along with other regulatory elements at the promoters, such as enhancers and transcriptional-factors, binding sites that form protein-complex regulators for DNA transcription . MAR regions control expression of genes, since they subdivide the chromatin into DNA loops, which mediate between Locus controls regions (LCRs) to distant promoters . Frequency of motifs in MAR elements are influenced by the A/T content  that participates in global transcription regulation by binding RNA polymerase to the nuclear matrix and creating the regulatory structure necessary for transcriptional initiation. Thus, the interaction between A/T-rich sites on the DNA to A/T-hook motif of the protein is an important component of the regulatory complex that controls transcription of genes .
The data presented in OPD will assist in the study of gene regulation by proteins that possess A/T-hook motifs. The OPD provides a broad analysis of the promoters essential for skeletal cell regulation and an accessible resource for the skeletal research community. This analysis will deepen and facilitate our understanding of the mechanism of transcriptional regulation of osteogenic cells. The construction of OPD, a bone-specific promoter database, will lead to the development of a novel skeletal-enhanced promoter-based microarray for studying gene expression in skeletal tissues at different stages of growth, development and the pathological state of malignancies and metabolic bone diseases.
Availability and requirements
Osteo-Promoter Database (OPD), available through web at http://www.OPD.tau.ac.il.
Rombauts S, Dehais P, Van Montagu M, Rouze P: PlantCARE, a plant cis-acting regulatory element database. Nucleic Acids Res. 1999, 27: 295-296. 10.1093/nar/27.1.295.
Higo K, Ugawa Y, Iwamoto M, Korenaga T: Plant cis-acting regulatory DNA elements (PLACE) database: 1999. Nucleic Acids Res. 1999, 27: 297-300. 10.1093/nar/27.1.297.
The FlyBase database of the Drosophila Genome Projects and community literature. The FlyBase Consortium. Nucleic Acids Res. 1999, 27: 85-88. 10.1093/nar/27.1.85.
Blake JA, Eppig JT, Richardson JE, Davisson MT: The Mouse Genome Database (MGD): expanding genetic and genomic resources for the laboratory mouse. The Mouse Genome Database Group. Nucleic Acids Res. 2000, 28: 108-111. 10.1093/nar/28.1.108.
Strausberg RL, Feingold EA, Klausner RD, Collins FS: The mammalian gene collection. Science. 1999, 286: 455-457. 10.1126/science.286.5439.455.
Geiman TM, Robertson KD: Chromatin remodeling, histone modifications, and DNA methylation-how does it all fit together?. J Cell Biochem. 2002, 87: 117-125. 10.1002/jcb.10286.
Aravind L, Landsman D: AT-hook motifs identified in a wide variety of DNA-binding proteins. Nucleic Acids Res. 1998, 26: 4413-4421. 10.1093/nar/26.19.4413.
Liebich I, Bode J, Reuter I, Wingender E: Evaluation of sequence motifs found in scaffold/matrix-attached regions (S/MARs). Nucleic Acids Res. 2002, 30: 3433-3442. 10.1093/nar/gkf446.
Liu R, States DJ: Consensus promoter identification in the human genome utilizing expressed gene markers and gene modeling. Genome Res. 2002, 12: 462-469. 10.1101/gr.198002.
Zhang MQ: Identification of human gene core promoters in silico. Genome Res. 1998, 8: 319-326.
Davuluri RV, Grosse I, Zhang MQ: Computational identification of promoters and first exons in the human genome. Nat Genet. 2001, 29: 412-417. 10.1038/ng780.
Scherf M, Klingenhoff A, Werner T: Highly specific localization of promoter regions in large genomic sequences by PromoterInspector: a novel context analysis approach. J Mol Biol. 2000, 297: 599-606. 10.1006/jmbi.2000.3589.
Hehl R, Wingender E: Database-assisted promoter analysis. Trends Plant Sci. 2001, 6: 251-255. 10.1016/S1360-1385(01)01954-9.
Quandt K, Frech K, Karas H, Wingender E, Werner T: MatInd and MatInspector: new fast and versatile tools for detection of consensus matches in nucleotide sequence data. Nucleic Acids Res. 1995, 23: 4878-4884.
Zhu J, Zhang MQ: SCPD: a promoter database of the yeast Saccharomyces cerevisiae. Bioinformatics. 1999, 15: 607-611. 10.1093/bioinformatics/15.7.607.
Praz V, Perier R, Bonnard C, Bucher P: The Eukaryotic Promoter Database, EPD: new entry types and links to gene expression data. Nucleic Acids Res. 2002, 30: 322-324. 10.1093/nar/30.1.322.
Liebich I, Bode J, Frisch M, Wingender E: S/MARt DB: a database on scaffold/matrix attached regions. Nucleic Acids Res. 2002, 30: 372-374. 10.1093/nar/30.1.372.
Frisch M, Frech K, Klingenhoff A, Cartharius K, Liebich I, Werner T: In silico prediction of scaffold/matrix attachment regions in large genomic sequences. Genome Res. 2002, 12: 349-354. 10.1101/gr.206602. Article published online before print in January 2002.
Jia L, Young MF, Powell J, Yang L, Ho NC, Hotchkiss R, Robey PG, Francomano CA: Gene expression profile of human bone marrow stromal cells: high-throughput expressed sequence tag sequencing analysis. Genomics. 2002, 79: 7-17. 10.1006/geno.2001.6683.
Jia L, Ho NC, Park SS, Powell J, Francomano CA: Comprehensive resource: Skeletal gene database. Am J Med Genet. 2001, 106: 275-281. 10.1002/ajmg.10227.
Singh GB, Kramer JA, Krawetz SA: Mathematical model to predict regions of chromatin attachment to the nuclear matrix. Nucleic Acids Res. 1997, 25: 1419-1425. 10.1093/nar/25.7.1419.
Katagiri T, Takahashi N: Regulatory mechanisms of osteoblast and osteoclast differentiation. Oral Dis. 2002, 8: 147-159. 10.1034/j.1601-0825.2002.01829.x.
Aubin JE: Regulation of osteoblast formation and function. Rev Endocr Metab Disord. 2001, 2: 81-94. 10.1023/A:1010011209064.
McCarthy TL, Ji C, Centrella M: Links among growth factors, hormones, and nuclear factors with essential roles in bone formation. Crit Rev Oral Biol Med. 2000, 11: 409-422.
Karsenty G: Minireview: transcriptional control of osteoblast differentiation. Endocrinology. 2001, 142: 2731-2733. 10.1210/en.142.7.2731.
Bucher P, Trifonov EN: Compilation and analysis of eukaryotic POL II promoter sequences. Nucleic Acids Res. 1986, 14: 10009-10026.
Bucher P: Weight matrix descriptions of four eukaryotic RNA polymerase II promoter elements derived from 502 unrelated promoter sequences. J Mol Biol. 1990, 212: 563-578. 10.1016/0022-2836(90)90223-9.
Maher JF, Nathans D: Multivalent DNA-binding properties of the HMG-1 proteins. Proc Natl Acad Sci U S A. 1996, 93: 6716-6720. 10.1073/pnas.93.13.6716.
Frank O, Schwanbeck R, Wisniewski JR: Protein footprinting reveals specific binding modes of a high mobility group protein I to DNAs of different conformation. J Biol Chem. 1998, 273: 20015-20020. 10.1074/jbc.273.32.20015.
Heng HH, Krawetz SA, Lu W, Bremer S, Liu G, Ye CJ: Re-defining the chromatin loop domain. Cytogenet Cell Genet. 2001, 93: 155-161. 10.1159/000056977.
This study is fulfillment for the requirement of MsC thesis for GI. Grants supported this study Greidinger for cancer research and the Chief Scientist, Ministry of Health to BD.
GI carried out the collections of genes in OPD, made the predication of genes' promoters, constructed the web site and drafted the manuscript. BD participated in the design of the study, conceived of the study, coordination and helped to draft the manuscript. All authors read and approved the final manuscript.