REEF: searching REgionally Enriched Features in genomes
© Coppe et al; licensee BioMed Central Ltd. 2006
Received: 15 September 2006
Accepted: 16 October 2006
Published: 16 October 2006
In Eukaryotic genomes, different features including genes are not uniformly distributed. The integration of annotation information and genomic position of functional DNA elements in the Eukaryotic genomes opened the way to test novel hypotheses of higher order genome organization and regulation of expression.
REEF is a new tool, aimed at identifying genomic regions enriched in specific features, such as a class or group of genes homogeneous for expression and/or functional characteristics. The method for the calculation of local feature enrichment uses test statistic based on the Hypergeometric Distribution applied genome-wide by using a sliding window approach and adopting the False Discovery Rate for controlling multiplicity. REEF software, source code and documentation are freely available at http://telethon.bio.unipd.it/bioinfo/reef/.
REEF can aid to shed light on the role of organization of specific genomic regions in the determination of their functional role.
In Eukaryotic genomes, different features including genes are not uniformly distributed. The integration of annotation information and genomic position of functional DNA elements in the Eukaryotic genomes opened the way to test novel hypotheses of higher order genome organization and regulation of expression. Different studies attempted at searching genomic regions enriched in specific gene categories and at correlating the localization of genes with their expression characteristics. In this way, clusters of highly expressed genes in the human genome were identified [1, 2], whereas frequent co-localization of co-expressed genes was found in the Drosophila genome  and, in a broad study, in six Eukaryotic genomes . Other studies described chromosomal clustering of muscle expressed genes and of reproductive genes in Caenorhabditis elegans [5, 6] and chromosomal domains of genes highly expressed in liver and others in colon, containing well conserved genes among human, mice and rat . The reconstruction of the human heart transcriptome showed then that most heart genes are clustered in small groups of elements, which show in their promoters similar regulatory elements and the products of which share more functional features than expected by chance . In these studies, diverse experimental schemas were adopted to search for local clusters of special gene categories and different systematic analyses were carried out by custom-tailored scripts applied on large-scale datasets. As far as program availability, to our knowledge dedicated software for this kind of analysis is still unavailable if we exclude some R functions specifically devoted to the analysis of differential gene expression along the human genome and to detect genomic regions involved in tumor-related imbalances [9–11] or yeast specific software . In this paper we present a novel program, called REEF (Regionally Enriched FEatures), allowing the systematic analysis of any genome for regions that show significant local enrichment of specific features and to display results in the context of genome annotation. REEF scans the considered genome using a sliding window approach and adopts the False Discovery Rate to give a genome-wide significance of observed local enrichment. A feature class may be established as a group of sequences or of genomic positions with specific characteristics. For instance a feature class may be the group of occurrences of a sequence motif in a genome or it may be a set of genes showing a specific expression pattern or a group of genes encoding products involved in a given biological process, in a molecular function, as encoded by Gene Ontology terms. Thus REEF can be used to inspect for local clustering of different kind of features, aiding to shed light on the role of organization of specific genomic regions in the determination of their functional role.
REEF aims at identifying regions of a genome enriched in specific features, as compared with a reference landscape of features density, by using a sliding window approach. REEF takes as input a list of reference features (RF; e.g. human genes) mapped to a genomic DNA sequence, a list of features selected (SF) among the RF (e.g. human genes specifically expressed in a given tissue), along with their genomic positions and the number and the length of chromosomes in the considered genome. Once selected the values for the size of the sliding window and for the shift between adjacent windows, the significance of regional enrichment in SF observed in each window, is calculated by using the Hypergeometric Distribution . Given the known distribution of RF, a p-value is associated to each window, corresponding to the probability of observing by chance a number of SF equal or greater than the observed. Let S be the total number of SF over the entire genome, R the total number of RF over the entire genome, and r the number of RF in a given window (with R ≥ r and S ≥ r). The probability of observing by chance at least k SF out of r RF in the window is the pointwise significance of the observed numbers of SF in the window (p-value, p):
The False Discovery Rate (FDR)  is used to circumvent the problem of multiple testing when calculating the genome-wide statistical significance for the observed enrichment in SF in a given region. After sorting windows by p-values over the entire genome, q-values (FDR) are calculated. Q (q-value) for each window is defined as Q = (p*N)/i, where p is the p-value of the window, N the total number of windows considered and i the number of windows with a p-value not higher than p. Given a global threshold for the genome wide FDR (e.g. 5%), the number of windows "significantly enriched in SF" is determined. The position and length of one cluster of SF is defined as the genome coverage of the maximum number of adjacent windows statistically significant.
REEF is an application written in Python , which runs both under Microsoft Windows and GNU/Linux. The wxPython multiplatform GUI toolkit is used to provide a graphical user interface. The software and its manual are freely available for download  [see additional files 1, 2, 3].
Display of results
Since REEF exploits UCSC Genome Browser Custom Annotation Tracks facility, clusters and features pertaining to each cluster are visualized as custom tracks, together with standard tracks from UCSC Genome Browser. A "cluster" track shows the position on the chromosome and the width of each given cluster, whereas, a "cluster elements" track shows the position and the width of the different SF in the cluster, identified by the name/ID given by the user.
Identification of clusters of tissue-specific genes in the human genome
For the sample analysis, we collected a set of 8,959 human genes, mapped on the genomic sequences and associated to good quality expression data in normal tissues. This group of genes was used as a set of RF, depicting the gene density along human chromosomes. Among these 8,959 genes, we selected a group of 896 genes expressed prevalently in single tissue, which constituted the set of SF. REEF was applied to this dataset, to search for clusters of tissue-specific genes in human chromosomes.
We considered the data of the GNFAtlas2 gene expression survey in 79 human tissues . From the "gnfHumanAtlas2Median table" of the UCSC Genome Browser annotation database hg18 (reporting for each tissue or cell type, median values of expression level measurements in different replicates) we selected expression data pertaining to ten differentiated human tissues under normal conditions (whole brain, lymph node, bone marrow, whole blood, skin, adipocyte, pancreas, skeletal muscle, lung and liver). We used only the subset of expression level measurements belonging to 33,148 Affimetrix_HG_133A probesets, for which quality checked annotation data were available at the GeneAnnot database of Weizmann Institute . GeneAnnot provides a revised and improved annotation of Affymetrix probesets, whose assignment to one or more GeneCards IDs is ranked by sensitivity and specificity scores of probeset/gene matches. We selected only probe-sets reported to be associated to an unique human gene, according to the GeneAnnot database (Release 13). Among them, 16,388, showing either sensitivity and specificity score of 1, were selected. By using GeneALaCart, 15,985 of them were associated to 10,539 EntrezGene IDs. For each of the 3,280 EntrezGene IDs (31%) represented by different probesets, we adopted a jackknifing procedure to integrate expression data in a unique expression vector. For each gene g, represented by N probesets, we calculated N vectors of median values, each for the group of probesets obtained by excluding the i probeset (with i ranging from 1 to N). If Spearman correlation between the calculated vector of median values and of the expression vector of the excluded probeset resulted higher than 0.5, the probeset was discarded. The remained vectors were used to calculate a vector of median values for the gene. In this way we obtained 8,959 unique expression vectors, each for a specific EntrezGene IDs. We then associated a genomic position to these genes, by using the start and end coordinates of the corresponding KnownGene at the UCSC database. For the majority of EntrezGene IDs resulting associated to more than one KnownGene ID, and thus to different genomic positions, the gene width was defined as the maximum genomic region covered by the different KnownGene elements pertaining to the given EntrezGene ID.
These 8,959 human genes were then ranked by their tendency of being regulated in a tissue-dependent manner, such as by their tissue-specificity. For each gene expression vector, we calculated the Shannon Entropy index (H) of tissue-specificity  and we selected the 10% of genes apparently most tissue-specific, which constituted the sample of 896 SF for subsequent analyses [see additional file 4]. The entire group of 8,959 genes constituted the RF sample [see additional files 5 and 6]. By running REEF on these data and selecting window size and shift of 500 and 100 Kb respectively, 44 significant (FDR set to 0.05) clusters of at least 3 tissue-specific genes were identified, accounting for 150 (17%) out of 896 genes included in clusters. Figures 1 and 2 show the graphical display of results produced by REEF.
Discussion and conclusion
REEF software aims at the identification of genomic regions containing features locally clustered. It can be useful, for example, to find regions enriched in genes with specific expression characteristics, where transcription regulation may operate at chromosomal domain level. As shown in the sample analysis, we identified by REEF clusters of tissue-specific genes in the human genome.
REEF graphical display facilities give several advantages and allow an intuitive and efficient evaluation of analysis results at a glance. Data are presented at different levels: at genome level, all clusters of locally enriched features are shown along the entire genome; at chromosome level, a plot of the quantity (1 - p-value) for each window in the selected chromosome is displayed. The chromosome plot outlines regions showing statistically significant local enrichment of the selected features (clusters) and may also reveal regions associated to a statistical score just lower than the threshold, but which may be of special interest since connecting strictly significant clusters in larger regions. Moreover, the integration of REEF with the UCSC Genome Browser Custom Annotation Tracks allows exploration of genomic regions pertaining to selected cluster, giving the possibility to reach all information about the region available at the UCSC Genome Browser database, such as DNA sequence and annotation regarding gene expression data, protein domain and structure information, Gene Ontology, homologous genes and interspecies sequence conservation.
Information on clusters position and content are also given in text format, thus allowing post-processing with custom made scripts or other available software. For instance the user can identify by REEF clusters of a special gene category, defined by analysis of gene expression data, but subsequently analyze obtained lists of clustered genes on the basis of additional different criteria such as Gene Ontology-based function annotation, analyzed to find terms enrichment within clusters.
Due to the very short execution time, subsequent analyses with different settings can be carried out quickly and the genome-wide visualization of results make easy to efficiently compare obtained results.
In order to obtain informative results, search parameters must be chosen appropriately. The "window width" parameter, changing the dimension of the window used to scan the genome by the sliding window approach, should be chosen taking into account the average dimension of considered features, the possible number of features per clusters and the density of features along the considered genomes. We have chosen the 500 Kb width, with about 6 genes per windows expected, since the median value of the span of genes included in our dataset is 26 Kb and the median distance between not overlapping genes is 64 Kb. The "window shift" parameter changes the distance between the starts of adjacent windows in the sliding window algorithm. Both the "window shift" and the "minimum number of features in cluster" parameter influence the FDR calculation, by acting directly on the number of windows on which the statistical test is calculated: a shorter "window shift" makes more stringent the statistical analysis, whereas rising the "minimum number of features in cluster" decreases the stringency.
As said in the Rationale, only scenario-specific programs are available to the researchers for the discovery of local enrichment of features/genes in genomes. The differences in aims and input data of these programs do not allow any direct evaluation of the performances of REEF by comparative algorithms experimentation.
Since REEF can be used on different genomes, to analyze different types of features (genes, DNA motifs, DNA sequences with specific characteristics), it represents an innovative and general tool for detecting the localization of genomic regions of clustered features, thus helping to deepen the knowledge on the architecture and function of genomes. REEF can be useful for studies aiming at evaluating the contribution of chromosomal patterning of co-expressed genes in gene expression regulation, during development and differentiation, with respect to the role of classical gene-specific induction/repression acting on proximal regulatory regions.
Availability and requirements
Project name: REEF
Project home page: http://telethon.bio.unipd.it/bioinfo/reef/
Operating system(s): Platform independent
Programming language: Python
Other requirements: python interpreter (version 2.3 or higher), wxPython GUI toolkit
Licence: GNU GPL
Any restrictions to use by non-academics: none
The authors acknowledge the financial support of the University of Padova and of the Italian Ministry of University and Scientific and Technological Development. They thank UCSC Genome Browser and GeneAnnot for providing data and Dr. Chiara Romualdi for helpful discussion on statistical methods. AC is recipient of a PhD fellowship of the PhD school in Biosciences of University of Padova.
- Caron H, van Schaik B, van der Mee M, Baas F, Riggins G, van Sluis P, Hermus MC, van Asperen R, Boon K, Voute PA, Heisterkamp S, van Kampen A, Versteeg R: The human transcriptome map: clustering of highly expressed genes in chromosomal domains. Science 2001, 291: 1289–1292. 10.1126/science.1056794View ArticlePubMedGoogle Scholar
- Versteeg R, van Schaik BD, van Batenburg MF, Roos M, Monajemi R, Caron H, Bussemaker HJ, van Kampen AH: The human transcriptome map reveals extremes in gene density, intron length, GC content, and repeat pattern for domains of highly and weakly expressed genes. Genome Res 2003, 13: 1998–2004. 10.1101/gr.1649303PubMed CentralView ArticlePubMedGoogle Scholar
- Spellman PT, Rubin GM: Evidence for large domains of similarly expressed genes in the Drosophila genome. J Biol 2002, 1: 5. 10.1186/1475-4924-1-5PubMed CentralView ArticlePubMedGoogle Scholar
- Fukuoka Y, Inaoka H, Kohane IS: Inter-species differences of co-expression of neighboring genes in eukaryotic genomes. BMC Genomics 2004, 5(1):4. 10.1186/1471-2164-5-4PubMed CentralView ArticlePubMedGoogle Scholar
- Roy PJ, Stuart JM, Lund J, Kim SK: Chromosomal clustering of muscle-expressed genes in Caenorhabditis elegans. Nature 2002, 418: 975–979.PubMedGoogle Scholar
- Miller MA, Cutter AD, Yamamoto I, Ward S, Greenstein D: Clustered organization of reproductive genes in the C. elegans genome. Curr Biol 2004, 14: 1284–1290. 10.1016/j.cub.2004.07.025View ArticlePubMedGoogle Scholar
- Yamashita T, Honda M, Takatori H, Nishino R, Hoshino N, Kaneko S: Genome-wide transcriptome mapping analysis identifies organ-specific gene expression patterns along human chromosomes. Genomics 2004, 84: 867–875. 10.1016/j.ygeno.2004.08.008View ArticlePubMedGoogle Scholar
- Vogel JH, von Heydebreck A, Purmann A, Sperling S: Chromosomal clustering of a human transcriptome reveals regulatory background. BMC Bioinformatics 2005, 6: 230. 10.1186/1471-2105-6-230PubMed CentralView ArticlePubMedGoogle Scholar
- Husing J, Zeschnigk M, Boes T, Jockel KH: Combining DNA expression with positional information to detect functional silencing of chromosomal regions. Bioinformatics 2003, 19: 2335–2342. 10.1093/bioinformatics/btg314View ArticlePubMedGoogle Scholar
- Zhou Y, Luoh SM, Zhang Y, Watanabe C, Wu TD, Ostland M, Wood WI, Zhang Z: Genome-wide identification of chromosomal regions of increased tumor expression by transcriptome analysis. Cancer Res 2003, 63: 5781–5784.PubMedGoogle Scholar
- Kano M, Nishimura K, Ishikawa S, Tsutsumi S, Hirota K, Hirose M, Aburatani H: Expression imbalance map: a new visualization method for detection of mRNA expression imbalance regions. Physiol Genomics 2003, 13: 31–46.View ArticlePubMedGoogle Scholar
- Chang CF, Wai KM, Patterton HG: Calculating the statistical significance of physical clusters of co-regulated genes in the genome: the role of chromatin in domain-wide gene regulation. Nucleic Acids Res 2004, 32: 1798–1807. 10.1093/nar/gkh507PubMed CentralView ArticlePubMedGoogle Scholar
- Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci USA 2003, 100: 9440–9445. 10.1073/pnas.1530509100PubMed CentralView ArticlePubMedGoogle Scholar
- The Python Programming Language[http://www.python.org]
- REEF website[http://telethon.bio.unipd.it/bioinfo/reef/]
- Su AI, Wiltshire T, Batalov S, Lapp H, Ching KA, Block D, Zhang J, Soden R, Hayakawa M, Kreiman G, Cooke MP, Walker JR, Hogenesch JB: A gene atlas of the mouse and human protein-encoding transcriptomes. Proc Natl Acad Sci USA 2004, 101: 6062–6067. 10.1073/pnas.0400782101PubMed CentralView ArticlePubMedGoogle Scholar
- Chalifa-Caspi V, Yanai I, Ophir R, Rosen N, Shmoish M, Benjamin-Rodrig H, Shklar M, Stein TI, Shmueli O, Safran M, Lancet D: GeneAnnot: comprehensive two-way linking between oligonucleotide array probesets and GeneCards genes. Bioinformatics 2004, 20: 1457–1458. 10.1093/bioinformatics/bth081View ArticlePubMedGoogle Scholar
- Schug J, Schuller WP, Kappen C, Salbaum JM, Bucan M, Stoeckert CJ Jr: Promoter features related to tissue specificity as measured by Shannon entropy. Genome Biol 2005, 6: R33. 10.1186/gb-2005-6-4-r33PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.