the changes of gene expression level after estrogen stimulation could be the result of hormone  receptor interaction via estrogen response elements ( eres ) , activation of the target gene transcription through interaction with co - activators or over non - genomic pathways ( 1,2 ) .
a number of genes whose transcription is affected by estrogen have been well documented in manually curated estrogen responsive genes database ( ergdb ) ( 3 ) .
although ergdb provides information about genes responsive to estrogen proved in experiments , it does not give any clue about the way these genes are regulated and what functional properties these genes may have .
thus , to facilitate functional analysis of estrogen responsive genes we do need a different framework not available in ergdb .
the new knowledgebase for estrogen responsive genes ( kberg ) provides users with such options to explore some of the genes ' functional properties and complements the experimental - information - rich ergdb.since the set of estrogen responsive genes covers a broad spectrum of functionality , belongs to many different pathways and may be related to many diseases , it is thus a significant challenge to characterize functions of these genes and to provide clues on how these genes are regulated and what potential effects we should expect . as a part of functional analysis ,
transcription regulation of these genes is an important step to unravel parts of the underlying estrogen driven signal pathways and to help develop insights of estrogen effects that can bring us closer to efficient therapies for estrogen - related diseases .
the binding of transcription factors ( tfs ) to specific regulatory elements is a primary mechanism for controlling gene transcription .
promoter region is the control center for such transcription regulation , where multiple transcription factor binding sites ( tfbss ) are clustered forming specific regulatory modules ( 4 ) .
a comparative genome study found that the complexity of the organism is not correlated to the number of genes ( 5 ) but to the complexity of gene regulation systems .
recent studies ( 6,7 ) document that the richness of transcript variants stemming from multiple transcription start sites ( tsss ) in the vast majority of mammalian gene loci are controlled in various ways at the promoter level . in brief
, dynamic gene regulation is achieved through multi - level manipulation in which transcription level control is an essential step .
the sophisticated gene regulatory modes characterized by a set of tfbss determine when , where and how gene will express based on the integration of internal and external signals .
broad gene groups are subjected to different regulatory programs as shown in ( 8) . thus ,
identification of the structure and organization of tfbss in promoters is the primary step towards gaining understanding of the gene 's transcription regulation .
however , detection of functional tfbss is challenging due to features such as short motif length , high degeneration and flexible distance location .
position weight matrix ( pwm ) approach has been widely used in tfbs detection although high false positive rate is the main drawback ( 9 ) .
comparative genomic approach that searches for highly preserved motifs across different species has been applied to improve the accuracy of prediction ( 10 ) .
the method has been corroborated in a genome - wide screen for high - affinity ere in human and mouse ( 11 ) . through analyzing remarkably consistent patterns in orthologous promoters ,
we aim to provide highly selective promoter motifs that could be involved in the control of transcription of genes responsive to estrogen .
in addition to the detection of conserved motifs in promoter regions , mapping these motifs to known tfbs provides hints of the interaction with tfs and the formation of transcription initiation complex .
as typical tfs , estrogen receptors ( ers ) influence their target genes by either directly binding to ere or indirectly interacting with a wide range of co - regulators forming complexes that bind to promoters ( 12 ) .
however , the experimental detection of co - regulators is time and labor consuming . in kberg
, we provided partial solution regarding these issues by matching the highly preserved motifs in promoters to the transfac database .
, we found that sp1-binding site is the most frequent tfbs in promoter regions in our collection .
this prediction result is in good concordance with promoter deletion studies that specifies the essential role of er / sp1 complex for gene 's response to estrogen at transcription level ( 13,14 ) .
the top 10 most frequently found tfbss are given on the kberg website under the section 
statistic. systematic assessment of gene functionality can help us to understand the complexity of living cells .
co - expression is a biological property that reflects functional relationship among a group of genes .
co - expressed genes tend to share certain regulatory mechanism by using similar sets of tfs .
the phenomenon of co - expression can also be a result of the biological dependency among a group of proteins .
it has been observed that proteins working closely in the same interaction network are also likely to be co - expressed under the same biological condition
. sets of nf-b target genes were overexpressed in aggressive breast cancer comparing with non - aggressive breast cancer ( 15 ) .
thus , clustering genes based on their expression profile is one method to potentially deduce gene functional characteristics .
another systematic approach to cluster gene for functional annotation is to group them based on gene ontology ( go ) ( 16 ) categories by using well - defined terminology to describe biological process , molecular function and cellular component where genes exert their activity .
owing to the important impact of estrogen in human physiology and its involvement in numerous diseases , we believe that functionally characterizing estrogen responsive genes can benefit both basic research and clinical applications .
however , up to date , there is no supporting system that facilitates such functional exploration . to complement currently available estrogen - related database and to provide researchers with comprehensive information related to potential functionality of estrogen responsive genes , we developed kberg
. users can obtain integrated information covering the following aspects : orthologous gene information from human , mouse and rat ; ab initio determined dna motifs preserved in promoter regions ; possible relation between motifs and known tfs based on the transfac database ( 17 ) .
the system also helps user to deduce functional association of gene groups annotated by the same go term or possessing the same expression characteristics . links to unigene , entrez gene , homologene , genbank , go ( 16 ) , evoc ( 18 ) and ergdb ( 3 ) make kberg a unique site for exploring estrogen responsive genes . at current stage , kberg contains information of 1362 homolog groups and 1526 estrogen responsive genes from human , mouse and rat .
although kberg is specialized for estrogen responsive gene , the framework can be adapted to other gene groups for their functional exploration and annotation .
kberg is free for academic and non - profit users and can be accessed at .
we have created a comprehensive system that presents functional properties of estrogen responsive genes from multiple viewpoints . all genes in kberg are experimentally proved to respond to estrogen  their expression levels are significantly altered by estrogen .
these genes were derived from a manually curated estrogen - specific database , ergdb ( 3 ) .
each module performs specialized functions that provide users a convenient way to find desired information .
the goal of this module is to identify highly conserved promoter motifs from groups of orthologous genes . to facilitate promoter analysis
, we built an in - house promoter database . using genes derived from ergdb as a primary list ,
orthologous genes from human ( homo sapiens ) , mouse ( mus musculus ) and rat ( rattus norvegicus ) genomes are downloaded from the ncbi homologene site ( ) .
recent genome comparison study ( 7 ) indicates that the availability of precise tsss increases the accuracy of finding conserved motifs by focusing on phylogenetically relevant promoter sequences . to accurately analyze functional motifs from selected promoters , we first identified tsss with solid experimental evidence .
each mouse and human tss is identified through multiple selection steps . supporting evidences include ( i ) full - length cdna from fantom3 database ( ) ; ( ii ) cage tags from riken tag - cluster database ( ) ; and ( iii ) 5-utr sequences based on ucsc human gene sorter database ( ) .
we limited search for conserved motifs to the region [ 1200 , + 500 ] relative to tss as most functional tfbss are located in these region , although some can be found several thousand base pairs away from the tss .
the shorter promoter sequence also reduces the noise in motif identification . dragon motif builder ( 19,20 )
the highly preserved motifs were detected without any presumption of what may be tfs that control activation of the respective genes .
however , sets of motifs within close mutual distance and in preserved orders in orthologous sequences are considered embedding more reliable and more significant information for transcriptional regulation .
we expect some of these predicted motifs will turn out to be functional tfbss , although further experimental validation is needed . in the report page , we provide both graphic overview of the motif arrangements and details of the motif sequence information , which enable users to further investigate these motifs experimentally , or to compare them to the known tfbs through other data sources .
the length of motifs we considered was 9 bp . in each orthologous promoter group , five motifs were selected based on their sequence preservation , motif combination order and maximum allowed distance between motifs ( 50 bp ) .
details of motif family consensus , matrix , positions , strands and rough graphic distribution in the orthologous promoter group are provided in the report page .
we also matched each motif to the transfac professional database , version 9.2 , to suggest possible tfbss that these motifs contain or overlap with .
we allowed one mismatch in 9 bp motif using the patch program ( transfac ) and the similarity score was given in the report page .
although the evaluation of motif prediction is beyond the scope of this paper , we would like to suggest the user to take full advantage of the information in the report page .
the matching score of potential tfbs with transfac , the total information content and position distribution , the degree of sequence preservation and the combination pattern with neighboring motifs are valuable indicators to estimate the accuracy of prediction .
this module enables users to group estrogen responsive genes that share common go category annotation . in this way genes in kberg
can be clustered into functionally related groups . using go terms we associate wide range of biological knowledge to groups of gene . by examining the frequency of the used go terms for the gene group
the  browse gene ontology function in kberg allows user to select estrogen responsive genes that share the same go annotation term and the statistic information about the term .
for example , there are 101 genes under the go term  receptor activity ( go:0004872 ) and only 22 genes belong to  steroid hormone receptor activity ( go:0003707 ) .
users can easily obtain two relevant but not identical lists of genes to study receptors with specific annotation and with broad annotation .
we believe that clustering genes according to go terms is a complementary way to associate genes with their functions and a logical way to assess genes .
allowing user to download the list of classified genes is one of the helpful features of the system .
evoc ( 18 ) associates human gene expression data with controlled ontology terms in sets of hierarchical vocabularies , covering categories such as anatomical system , cell type , pathology , developmental stage , etc .
we use evoc ontology in our system to provide an alternative approach to cluster gene for functional assessment .
the expression data provide extensive information about the gene co - expression under various biological conditions , which , although coarse , becomes essential indicators for gene co - regulation and interaction .
expression results also represent important gene behaviors that should not be missed out in gene functional characterization . in kberg , we integrated evoc ( version 2.7 ) to reveal gene expression in terms of human organ , cell type , disease and age
. thus , a user can easily find out some of the associated information for gene expression ( organ and cell type , whether in normal or diseased tissues and in the age group ) .
the downloadable list of classified genes is provided and could be useful for user 's further exploration . designed for estrogen - related research
, kberg provides a wide range of information about estrogen responsive genes through promoter motifs and functional associations via well - defined vocabularies from go and evoc .
we summarized the main functions of the system in table 1 . in brief , through searching this system , within a short period of time , a user can gain valuable information related to either individual gene or a group of genes that share the same ontology description .
however , there is limitation of the system that user needs to bear in mind .
the promoter motif prediction is a pure computational estimation and needs to be verified by experiments .
also , classifying gene based on ontology categories may not always be the best way to reveal its biological function .
the goal of this module is to identify highly conserved promoter motifs from groups of orthologous genes . to facilitate promoter analysis
, we built an in - house promoter database . using genes derived from ergdb as a primary list ,
orthologous genes from human ( homo sapiens ) , mouse ( mus musculus ) and rat ( rattus norvegicus ) genomes are downloaded from the ncbi homologene site ( ) .
recent genome comparison study ( 7 ) indicates that the availability of precise tsss increases the accuracy of finding conserved motifs by focusing on phylogenetically relevant promoter sequences . to accurately analyze functional motifs from selected promoters , we first identified tsss with solid experimental evidence .
each mouse and human tss is identified through multiple selection steps . supporting evidences include ( i ) full - length cdna from fantom3 database ( ) ; ( ii ) cage tags from riken tag - cluster database ( ) ; and ( iii ) 5-utr sequences based on ucsc human gene sorter database ( ) .
we limited search for conserved motifs to the region [ 1200 , + 500 ] relative to tss as most functional tfbss are located in these region , although some can be found several thousand base pairs away from the tss . the shorter promoter sequence also reduces the noise in motif identification . dragon motif builder ( 19,20 ) was used to identify ab initio dna motifs in each orthologous promoter sequence group .
the highly preserved motifs were detected without any presumption of what may be tfs that control activation of the respective genes .
however , sets of motifs within close mutual distance and in preserved orders in orthologous sequences are considered embedding more reliable and more significant information for transcriptional regulation .
we expect some of these predicted motifs will turn out to be functional tfbss , although further experimental validation is needed . in the report page , we provide both graphic overview of the motif arrangements and details of the motif sequence information , which enable users to further investigate these motifs experimentally , or to compare them to the known tfbs through other data sources .
the length of motifs we considered was 9 bp . in each orthologous promoter group , five motifs were selected based on their sequence preservation , motif combination order and maximum allowed distance between motifs ( 50 bp ) .
details of motif family consensus , matrix , positions , strands and rough graphic distribution in the orthologous promoter group are provided in the report page .
we also matched each motif to the transfac professional database , version 9.2 , to suggest possible tfbss that these motifs contain or overlap with .
we allowed one mismatch in 9 bp motif using the patch program ( transfac ) and the similarity score was given in the report page .
although the evaluation of motif prediction is beyond the scope of this paper , we would like to suggest the user to take full advantage of the information in the report page .
the matching score of potential tfbs with transfac , the total information content and position distribution , the degree of sequence preservation and the combination pattern with neighboring motifs are valuable indicators to estimate the accuracy of prediction .
this module enables users to group estrogen responsive genes that share common go category annotation . in this way genes in kberg
can be clustered into functionally related groups . using go terms we associate wide range of biological knowledge to groups of gene . by examining the frequency of the used go terms for the gene group
the  browse gene ontology function in kberg allows user to select estrogen responsive genes that share the same go annotation term and the statistic information about the term .
for example , there are 101 genes under the go term  receptor activity ( go:0004872 ) and only 22 genes belong to  steroid hormone receptor activity ( go:0003707 ) .
users can easily obtain two relevant but not identical lists of genes to study receptors with specific annotation and with broad annotation .
we believe that clustering genes according to go terms is a complementary way to associate genes with their functions and a logical way to assess genes .
allowing user to download the list of classified genes is one of the helpful features of the system .
evoc ( 18 ) associates human gene expression data with controlled ontology terms in sets of hierarchical vocabularies , covering categories such as anatomical system , cell type , pathology , developmental stage , etc .
we use evoc ontology in our system to provide an alternative approach to cluster gene for functional assessment .
the expression data provide extensive information about the gene co - expression under various biological conditions , which , although coarse , becomes essential indicators for gene co - regulation and interaction .
expression results also represent important gene behaviors that should not be missed out in gene functional characterization . in kberg
, we integrated evoc ( version 2.7 ) to reveal gene expression in terms of human organ , cell type , disease and age .
thus , a user can easily find out some of the associated information for gene expression ( organ and cell type , whether in normal or diseased tissues and in the age group ) .
the downloadable list of classified genes is provided and could be useful for user 's further exploration .
designed for estrogen - related research , kberg provides a wide range of information about estrogen responsive genes through promoter motifs and functional associations via well - defined vocabularies from go and evoc .
we summarized the main functions of the system in table 1 . in brief , through searching this system , within a short period of time , a user can gain valuable information related to either individual gene or a group of genes that share the same ontology description . however , there is limitation of the system that user needs to bear in mind .
the promoter motif prediction is a pure computational estimation and needs to be verified by experiments .
also , classifying gene based on ontology categories may not always be the best way to reveal its biological function .
 estrogen has a profound impact on human physiology affecting transcription of numerous genes . to decipher functional characteristics of estrogen responsive genes , we developed knowledgebase for estrogen responsive genes ( kberg ) . 
 genes in kberg were derived from estrogen responsive gene database ( ergdb ) and were analyzed from multiple aspects . 
 we explored the possible transcription regulation mechanism by capturing highly conserved promoter motifs across orthologous genes , using promoter regions that cover the range of [ 1200 , + 500 ] relative to the transcription start sites . 
 the motif detection is based on ab initio discovery of common cis - elements from the orthologous gene cluster from human , mouse and rat , thus reflecting a degree of promoter sequence preservation during evolution . 
 the identified motifs are linked to transcription factor binding sites based on the transfac database . 
 in addition , kberg uses two established ontology systems , go and evoc , to associate genes with their function . 
 users may assess gene functionality through the description terms in go . 
 alternatively , they can gain gene co - expression information through evidence from human est libraries via evoc . 
 kberg is a user - friendly system that provides links to other relevant resources such as ergdb , unigene , entrez gene , homologene , go , evoc and genbank , and thus offers a platform for functional exploration and potential annotation of genes responsive to estrogen . 
 kberg database can be accessed at . 
INTRODUCTION
DATABASE CONSTRUCTION
Promoter analysis and motif detection module
Gene Ontology grouping module
eVOC module
SUMMARY