BACKGROUND
high-throughput assays are being widely employed in various fields of biology.
for example in genomics, dna microarrays are used to simultaneously measure the expression levels of nearly all genes of a genome  <cit> .
recently, a new high-throughput method has been developed for a whole genome mapping of protein-dna interactions that is based on the chromatin immunoprecipitation and next generation sequencing technology   <cit> .
these two high-throughput methods, when combined, are instrumental to study how transcription factors regulate gene expression at a global genomic scale.
as the costs of new-generation sequencing and dna microarrays decrease, such high-throughput assays should be increasingly used.
in addition, new software tools are emerging rapidly, allowing faster and easier analyses of large-scale genomic datasets  <cit> .
however, extracting the significant and biologically relevant information from such massive datasets still represents a great challenge.
in purely experimental studies, the use of negative controls such as the blank or mock conditions is absolutely necessary.
however, genome-wide computer analyses may lack an adequate negative control.
in such case, a randomly selected portion of the total dataset can be used as an in silico negative control, such as for instance a randomly picked sample of all genomic loci.
the size of the control dataset can then be chosen to parallel that of the experimental set, to simplify the statistical analysis  <cit> .
as a study group of proteins, we chose nfi/ctf family of mammalian transcription factors.
nfi/ctf represents a family of transcription-replication factors comprising polypeptides encoded by four paralogous genes located on different chromosomes in mammals   <cit> .
nfi family of proteins displays the unusual property of regulating not only the initiation of transcription but also of mediating dna replication  <cit> .
nfi recognition sequence were found in the promoter sequences of many cellular genes  <cit> , where they may act as activator or repressor of transcription  <cit> .
recently, it has been proposed that nfi may be involved in a long range regulation of gene expression, through the formation of chromatin barrier and by blocking the propagation of a heterochromatic structure  <cit> .
nfi binds as a dimer, and its preferred binding sequence is a palindrome composed of two half sites ttggcannntgccaa.
a position weight matrix for the nfi/ctf was established using a collection of over  <dig>  selex-sage selected sites, allowing the prediction of its binding affinity to any genomic sequence  <cit> .
however, since this prediction matrix is based on nfi binding specificity in vitro, the specificity of this family of proteins may be different from that observed in the cell, where interactions with other transcription factors may take place and dna accessibility may be restrained by chromatin.
here we assessed the in vivo binding preferences of nfi/ctf, its global functional properties regarding the regulation of gene expression and the relationship of nfi binding sites with different histone methylation markers typical of either an open or closed chromatin structure.
RESULTS
nfi preferentially binds upstream of transcription initiation sites in mouse genome
statistical analysis in genomics often relies on the sub-sampling of datasets, which requires random sampling algorithms.
we devised a random sampling algorithm that can be conveniently applied to large genomic datasets.
the random sampling algorithm c++ source code is available as a text file online .
each randomly generated number is used to extract an entry line from the main dataset to generate a subset of the desired size.
simulation experiments indicated that subsampling can be applied to sets of normally distributed values without loosing the statistical robustness of the comparisons, provided that relatively large subsets of data are retained .
we first used this random sampling method to compare data from a chip-seq experiment performed on primary mouse embryonic fibroblasts for the nfi transcription factor family relative to the in silico predictions of its binding sites.
the mouse genome  was found to contain a set of  <dig>  nfi predicted sites that were defined using a previously established position weight matrix  <cit> .
the predicted sites were defined with a matrix score threshold >  <dig>  which corresponds to a medium in vitro binding affinity in the range of scores that extends from a minimum of - <dig> to a maximum of  <dig>  within this set,  <dig>  predicted sites overlapped dna sequences covering the refseq annotated transcription start sites  and  <dig> kb of upstream sequences.
actual in vivo binding site occupancy was estimated using the average chip-seq tag count covering each subset of nfi predicted sites.
when tag counts occurring at tss-proximal binding sites were compared to those that do not occur near known tss, a higher tag counts was noted at the predicted sites, but also in regions extending several kb away from the known binding sites .
since the dna fragments used in this assay did not cover a size range extending over several kb, we reasoned that this might either result from a surprisingly frequent association of the protein to sequences near binding sites, or, alternatively, that it may stem from the smaller size of the dataset comprising the tss-proximal  <dig>  binding sites relative to the non-tss datasets of  <dig>  sequences, resulting in a relatively noisier distribution at non-bound sequences, or from some other artefactual effects associated to the method.
to assess the latter possibility further, we generated control datasets of the same size as the experimental set by randomly selecting  <dig>  sequences from the  <dig>  predicted sites.
comparison of distinct randomly selected data subsets indicated comparable tag counts at predicted binding sites, and similar signal-to-noise ratio.
thus, the sampling method provided a reliable estimation of the binding site occupancy, since several random groups of predicted sites did not differ markedly in their protein occupancy .
as before, the dataset corresponding to tss-proximal sites showed a more prominent tag count around the predicted sites.
however, comparable tag counts were observed within  <dig> to  <dig> bp windows around the predicted binding sites, when comparing the experimental profiles to those of control datasets of the same size.
thus, we conclude that the group of tss-proximal predicted sites displayed higher protein occupancy at predicted sites when compared to the randomly selected binding sites, but that it was not the overall promoter region that was more frequently bound by the protein.
this suggested that binding sites within genomic loci upstream of core promoters may bind nfi with higher apparent affinity.
nfi-bound genes show higher expression levels
next, we assessed whether genes that bind nfi may have some common features regarding their gene expression levels.
we selected  <dig>  in vivo nfi sites genome-wide from the chip-seq data, which corresponds to the genomic loci obtained from the collection of precipitated dna fragments.
out of this number,  <dig>  in vivo sites were located within  <dig> kb upstream regions of refseq annotated genes.
as multiple in vivo nfi sites may occur within the  <dig> kb upstream regions, a total of  <dig> refseq-annotated genes were identified to contain one or more in vivo nfi sites.
control groups were generated to consist of  <dig>  genes randomly selected from the total database of refseq genes.
the expression levels of these genes was assessed using microarray profiling data of mrnas obtained from murine embryo fibroblasts  <cit> .
again, randomly selected groups of genes did not differ significantly in their overall expression levels or in their distribution profiles, with the most prominent peak corresponding to lowly expressed genes .
however, the transcriptional levels of the group of nfi-containing genes were distinctly higher than those of the control subgroups , with a median expression value of  <dig>  for nfi-bound genes versus  <dig> ,  <dig>  and  <dig>  for the randomly selected groups of genes .
thus, we conclude that nfi preferentially occupies promoters or upstream regulatory regions of genes that exhibit high expression levels in the cell.
nfi binding correlates with specific histone methylation patterns
to test whether the preferential binding of nfi to upstream regions of expressed genes may be associated with specific chromatin modifications, we used chip-seq data for  <dig> different histone h <dig> methylations, as previously obtained from mouse embryonic fibroblasts  <cit> .
we created average chip-seq tag profiles surrounding the tss for each of these modifications.
the group of  <dig>  nfi bound genes showed higher levels of h3k4me <dig> and h3k36me <dig>  which are markers of active promoters and transcribed regions, respectively, when compared to the three groups of  <dig>  randomly chosen genes .
interestingly, nfi binding at tss did not correlate with a modification often associated with non-expressed or bivalent chromatin, h3k27me <dig>  when compared with the profiles of the control gene subgroups , indicating that nfi binding does not occur at promoter regions that are generally enriched in all histone modifications.
interestingly, we found h3k9me <dig>  a modification recently associated with nfi chromatin-domain boundary activity at telomeres  <cit> , to be also enriched surrounding the tss of the group of nfi-bound genes .
nfi-c most often acts as an activator of gene expression
the preferential association of nfi with expressed genes prompted us to test whether nfi family members may have the potential to activate gene expression, or whether their binding is rather the consequence of histone modifications and chromatin accessibility at promoters.
for this purpose, we used gene expression microarray data from mouse embryonic fibroblasts in which one of the nfi factors  was knocked out by insertional mutagenesis  <cit> .
we focused on differences in expression levels when comparing wild-type  and knock-out cells , and defined a set of  <dig> genes which expression was most decreased in the absence of nfi-c .
we also defined the set of  <dig> genes whose expression is most increased in the absence of nfi-c .
we also randomly sampled  <dig> genes from all refseq genes as control sets.
differences in the microarray signals of wt and ko cells  were determined for each gene within defined groups and plotted as an absolute value.
the group of genes activated by nfi-c showed the highest median Δe value, implicating this factor as a potent gene activator .
nfi-c suppressed genes also appeared to be significantly regulated as compared to the control groups.
however, nfi-c suppressed genes displayed moderate Δe values, between those of the control groups and of the group of nfi-c activated genes, indicating that nfi-c is globally a less potent inhibitor of gene expression.
when considering the expression levels in wt cells , the group of nfi-c suppressed genes had levels of expression that did not differ significantly from the control subgroups , again indicating that nfi-c acts as a moderate inhibitor of gene expression .
overall, the group of nfi-c activated genes had significantly higher levels of expression than the control groups .
these results suggest that nfi-c member of the nfi family acts globally as an effective gene activator.
discussion and
CONCLUSIONS
in this study, we used a random sampling procedure as a general method to obtain reliable control datasets in the analysis of high-throughput genomic assays.
we find that datasets of more than  <dig> individual values can be used without decreasing the robustness of statistical analysis, and that independently generated random subsets of data have statistically indistinguishable global properties.
thus, subsampling can provide a convenient way to display and compare the noise and signals from experimental and control datasets of the same size.
first, we showed that nfi binds preferentially those predicted sites that are located upstream of the initiation sites of transcription .
several interpretations may be given to the preferential association of nfi to binding sites in the proximity of tss rather than to other locations of the genome.
it is known that nfi occupies the promoters of many genes where it may bind synergistically with some other transcription factors such as hepatocyte nuclear factor  <dig> alpha, estrogen receptor, brg-associated factor  <cit> .
thus, the preferred occupancy of tss proximal sites may at least in part reflect the synergistic association of nfi with other factors.
we also found that nfi occupy promoters or upstream regions of the group of genes that are significantly more expressed than the representative randomly selected control groups.
since correlation does not necessarily imply causal relationship, this observation does not allow the conclusion that nfi-family members actually activate the expression of these genes.
for instance, nfi might bind highly expressed genes to suppress in part their expression, but still leaving relatively high transcription levels.
however, taken together with previous observations that nfi activates the expression of many genes in higher eukaryotes  <cit> , we rather conclude that the observed correlation may originate from a direct up-regulation of gene expression by nfi, at least for a significant proportion of its target genes.
the hypothesis that nfi family members may directly activate genes appears to be true for at least one of the member of the family , as mrna profiling analysis performed on wild-type and nfi-c knock-out cells revealed that nfi-c is a more potent gene activator than a repressor.
the  <dig> genes that are most up-regulated by nfi-c had significantly higher change in their expression levels than the top  <dig> down-regulated genes.
in addition, up-regulated genes showed significantly higher expression levels than representative control gene samples selected from the total gene population, implicating again that this factor is a potent activator of gene expression.
since the selected in vivo nfi binding sites are located up to  <dig> kb from their tss, which is a relatively large distance, nfi might act as well through some of the types of remote regulation, for instance by the establishment of a chromatin domain boundary that would prevent the propagation of a silencing chromatin structure towards the promoter  <cit> .
histone h <dig> methylations such as the h3k4me <dig> and h3k36me <dig> modifications were found to be enriched around the tss of nfi-occupied genes when compared with control gene groups.
this finding is consistent with the model that nfi acts predominantly as an activator of transcription, since h3k4me <dig> and h3k36me <dig>  but not h3k27me <dig>  were proposed as markers of active gene transcription  <cit> .
this indicates that nfi binding to the upstream regions may contribute to the recruitment of the specific enzymes for the h3k4me <dig> and h3k36me <dig> modifications.
a genome-wide correlation of the occurrence of h3k27me <dig> was also observed around tss occurring close to nfi-bound sites, however it was indistinguishable to that of the control group of genes.
this indicates that this correlation results from an enrichment of h3k27me <dig> around at least some of the tss, and that nfi is not involved in the recruitment of enzymes mediating this modification.
thus, the enrichment of h3k27me <dig> modification over the nfi bound genes represents a false positive genome-wide correlation.
interestingly, we also found the h3k9me <dig> modification to be slightly enriched in the group of nfi bound genes.
although h3k9me <dig> has been associated with a closed chromatin structure, this suggests that nfi may be involved in the recruitment of enzymes that mediate this modification.
interestingly, this modification was recently associated with a chromatin domain boundary effect at telomeric regions in human cells  <cit> .
in this study, nfi was shown to prevent the propagation of a silencing chromatin structure from the telomere, and the expressed genes protected from telomeric silencing by nfi were shown to have elevated h3k9me <dig> marks at specific telomeric positions.
thus, we may conclude from these studies that the enrichment in h3k9me <dig> may be a hallmark of gene expression activation by nfi.

 high-throughput assays are being widely employed in various fields of biology. for example in genomics, dna microarrays are used to simultaneously measure the expression levels of nearly all genes of a genome. recently, a new high-throughput method has been developed for a whole genome mapping of protein-dna interactions that is based on the chromatin immunoprecipitation and next generation sequencing technology. these two high-throughput methods, when combined, are instrumental to study how transcription factors regulate gene expression at a global genomic scale. as the costs of new-generation sequencing and dna microarrays decrease, such high-throughput assays should be increasingly used. in addition, new software tools are emerging rapidly, allowing faster and easier analyses of large-scale genomic datasets. however, extracting the significant and biologically relevant information from such massive datasets still represents a great challenge. in purely experimental studies, the use of negative controls such as the blank or mock conditions is absolutely necessary. however, genome-wide computer analyses may lack an adequate negative control. in such case, a randomly selected portion of the total dataset can be used as an in silico negative control, such as for instance a randomly picked sample of all genomic loci. the size of the control dataset can then be chosen to parallel that of the experimental set, to simplify the statistical analysis. as a study group of proteins, we chose nfi/ctf family of mammalian transcription factors. nfi/ctf represents a family of transcription-replication factors comprising polypeptides encoded by four paralogous genes located on different chromosomes in mammals. nfi family of proteins displays the unusual property of regulating not only the initiation of transcription but also of mediating dna replication. nfi recognition sequence were found in the promoter sequences of many cellular genes, where they may act as activator or repressor of transcription. recently, it has been proposed that nfi may be involved in a long range regulation of gene expression, through the formation of chromatin barrier and by blocking the propagation of a heterochromatic structure. nfi binds as a dimer, and its preferred binding sequence is a palindrome composed of two half sites ttggcannntgccaa. a position weight matrix for the nfi/ctf was established using a collection of over selex-sage selected sites, allowing the prediction of its binding affinity to any genomic sequence. however, since this prediction matrix is based on nfi binding specificity in vitro, the specificity of this family of proteins may be different from that observed in the cell, where interactions with other transcription factors may take place and dna accessibility may be restrained by chromatin. here we assessed the in vivo binding preferences of nfi/ctf, its global functional properties regarding the regulation of gene expression and the relationship of nfi binding sites with different histone methylation markers typical of either an open or closed chromatin structure.  nfi preferentially binds upstream of transcription initiation sites in mouse genome statistical analysis in genomics often relies on the sub-sampling of datasets, which requires random sampling algorithms. we devised a random sampling algorithm that can be conveniently applied to large genomic datasets. the random sampling algorithm c++ source code is available as a text file online. each randomly generated number is used to extract an entry line from the main dataset to generate a subset of the desired size. simulation experiments indicated that subsampling can be applied to sets of normally distributed values without loosing the statistical robustness of the comparisons, provided that relatively large subsets of data are retained. we first used this random sampling method to compare data from a chip-seq experiment performed on primary mouse embryonic fibroblasts for the nfi transcription factor family relative to the in silico predictions of its binding sites. the mouse genome was found to contain a set of nfi predicted sites that were defined using a previously established position weight matrix. the predicted sites were defined with a matrix score threshold > which corresponds to a medium in vitro binding affinity in the range of scores that extends from a minimum of - to a maximum of within this set, predicted sites overlapped dna sequences covering the refseq annotated transcription start sites and kb of upstream sequences. actual in vivo binding site occupancy was estimated using the average chip-seq tag count covering each subset of nfi predicted sites. when tag counts occurring at tss-proximal binding sites were compared to those that do not occur near known tss, a higher tag counts was noted at the predicted sites, but also in regions extending several kb away from the known binding sites. since the dna fragments used in this assay did not cover a size range extending over several kb, we reasoned that this might either result from a surprisingly frequent association of the protein to sequences near binding sites, or, alternatively, that it may stem from the smaller size of the dataset comprising the tss-proximal binding sites relative to the non-tss datasets of sequences, resulting in a relatively noisier distribution at non-bound sequences, or from some other artefactual effects associated to the method. to assess the latter possibility further, we generated control datasets of the same size as the experimental set by randomly selecting sequences from the predicted sites. comparison of distinct randomly selected data subsets indicated comparable tag counts at predicted binding sites, and similar signal-to-noise ratio. thus, the sampling method provided a reliable estimation of the binding site occupancy, since several random groups of predicted sites did not differ markedly in their protein occupancy. as before, the dataset corresponding to tss-proximal sites showed a more prominent tag count around the predicted sites. however, comparable tag counts were observed within to bp windows around the predicted binding sites, when comparing the experimental profiles to those of control datasets of the same size. thus, we conclude that the group of tss-proximal predicted sites displayed higher protein occupancy at predicted sites when compared to the randomly selected binding sites, but that it was not the overall promoter region that was more frequently bound by the protein. this suggested that binding sites within genomic loci upstream of core promoters may bind nfi with higher apparent affinity. nfi-bound genes show higher expression levels next, we assessed whether genes that bind nfi may have some common features regarding their gene expression levels. we selected in vivo nfi sites genome-wide from the chip-seq data, which corresponds to the genomic loci obtained from the collection of precipitated dna fragments. out of this number, in vivo sites were located within kb upstream regions of refseq annotated genes. as multiple in vivo nfi sites may occur within the kb upstream regions, a total of refseq-annotated genes were identified to contain one or more in vivo nfi sites. control groups were generated to consist of genes randomly selected from the total database of refseq genes. the expression levels of these genes was assessed using microarray profiling data of mrnas obtained from murine embryo fibroblasts. again, randomly selected groups of genes did not differ significantly in their overall expression levels or in their distribution profiles, with the most prominent peak corresponding to lowly expressed genes. however, the transcriptional levels of the group of nfi-containing genes were distinctly higher than those of the control subgroups, with a median expression value of for nfi-bound genes versus, and for the randomly selected groups of genes. thus, we conclude that nfi preferentially occupies promoters or upstream regulatory regions of genes that exhibit high expression levels in the cell. nfi binding correlates with specific histone methylation patterns to test whether the preferential binding of nfi to upstream regions of expressed genes may be associated with specific chromatin modifications, we used chip-seq data for different histone h methylations, as previously obtained from mouse embryonic fibroblasts. we created average chip-seq tag profiles surrounding the tss for each of these modifications. the group of nfi bound genes showed higher levels of h3k4me and h3k36me which are markers of active promoters and transcribed regions, respectively, when compared to the three groups of randomly chosen genes. interestingly, nfi binding at tss did not correlate with a modification often associated with non-expressed or bivalent chromatin, h3k27me when compared with the profiles of the control gene subgroups, indicating that nfi binding does not occur at promoter regions that are generally enriched in all histone modifications. interestingly, we found h3k9me a modification recently associated with nfi chromatin-domain boundary activity at telomeres, to be also enriched surrounding the tss of the group of nfi-bound genes. nfi-c most often acts as an activator of gene expression the preferential association of nfi with expressed genes prompted us to test whether nfi family members may have the potential to activate gene expression, or whether their binding is rather the consequence of histone modifications and chromatin accessibility at promoters. for this purpose, we used gene expression microarray data from mouse embryonic fibroblasts in which one of the nfi factors was knocked out by insertional mutagenesis. we focused on differences in expression levels when comparing wild-type and knock-out cells, and defined a set of genes which expression was most decreased in the absence of nfi-c. we also defined the set of genes whose expression is most increased in the absence of nfi-c. we also randomly sampled genes from all refseq genes as control sets. differences in the microarray signals of wt and ko cells were determined for each gene within defined groups and plotted as an absolute value. the group of genes activated by nfi-c showed the highest median Δe value, implicating this factor as a potent gene activator. nfi-c suppressed genes also appeared to be significantly regulated as compared to the control groups. however, nfi-c suppressed genes displayed moderate Δe values, between those of the control groups and of the group of nfi-c activated genes, indicating that nfi-c is globally a less potent inhibitor of gene expression. when considering the expression levels in wt cells, the group of nfi-c suppressed genes had levels of expression that did not differ significantly from the control subgroups, again indicating that nfi-c acts as a moderate inhibitor of gene expression. overall, the group of nfi-c activated genes had significantly higher levels of expression than the control groups. these results suggest that nfi-c member of the nfi family acts globally as an effective gene activator. discussion and  in this study, we used a random sampling procedure as a general method to obtain reliable control datasets in the analysis of high-throughput genomic assays. we find that datasets of more than individual values can be used without decreasing the robustness of statistical analysis, and that independently generated random subsets of data have statistically indistinguishable global properties. thus, subsampling can provide a convenient way to display and compare the noise and signals from experimental and control datasets of the same size. first, we showed that nfi binds preferentially those predicted sites that are located upstream of the initiation sites of transcription. several interpretations may be given to the preferential association of nfi to binding sites in the proximity of tss rather than to other locations of the genome. it is known that nfi occupies the promoters of many genes where it may bind synergistically with some other transcription factors such as hepatocyte nuclear factor alpha, estrogen receptor, brg-associated factor. thus, the preferred occupancy of tss proximal sites may at least in part reflect the synergistic association of nfi with other factors. we also found that nfi occupy promoters or upstream regions of the group of genes that are significantly more expressed than the representative randomly selected control groups. since correlation does not necessarily imply causal relationship, this observation does not allow the conclusion that nfi-family members actually activate the expression of these genes. for instance, nfi might bind highly expressed genes to suppress in part their expression, but still leaving relatively high transcription levels. however, taken together with previous observations that nfi activates the expression of many genes in higher eukaryotes, we rather conclude that the observed correlation may originate from a direct up-regulation of gene expression by nfi, at least for a significant proportion of its target genes. the hypothesis that nfi family members may directly activate genes appears to be true for at least one of the member of the family, as mrna profiling analysis performed on wild-type and nfi-c knock-out cells revealed that nfi-c is a more potent gene activator than a repressor. the genes that are most up-regulated by nfi-c had significantly higher change in their expression levels than the top down-regulated genes. in addition, up-regulated genes showed significantly higher expression levels than representative control gene samples selected from the total gene population, implicating again that this factor is a potent activator of gene expression. since the selected in vivo nfi binding sites are located up to kb from their tss, which is a relatively large distance, nfi might act as well through some of the types of remote regulation, for instance by the establishment of a chromatin domain boundary that would prevent the propagation of a silencing chromatin structure towards the promoter. histone h methylations such as the h3k4me and h3k36me modifications were found to be enriched around the tss of nfi-occupied genes when compared with control gene groups. this finding is consistent with the model that nfi acts predominantly as an activator of transcription, since h3k4me and h3k36me but not h3k27me were proposed as markers of active gene transcription. this indicates that nfi binding to the upstream regions may contribute to the recruitment of the specific enzymes for the h3k4me and h3k36me modifications. a genome-wide correlation of the occurrence of h3k27me was also observed around tss occurring close to nfi-bound sites, however it was indistinguishable to that of the control group of genes. this indicates that this correlation results from an enrichment of h3k27me around at least some of the tss, and that nfi is not involved in the recruitment of enzymes mediating this modification. thus, the enrichment of h3k27me modification over the nfi bound genes represents a false positive genome-wide correlation. interestingly, we also found the h3k9me modification to be slightly enriched in the group of nfi bound genes. although h3k9me has been associated with a closed chromatin structure, this suggests that nfi may be involved in the recruitment of enzymes that mediate this modification. interestingly, this modification was recently associated with a chromatin domain boundary effect at telomeric regions in human cells. in this study, nfi was shown to prevent the propagation of a silencing chromatin structure from the telomere, and the expressed genes protected from telomeric silencing by nfi were shown to have elevated h3k9me marks at specific telomeric positions. thus, we may conclude from these studies that the enrichment in h3k9me may be a hallmark of gene expression activation by nfi. 
as a study model, we chose the nfi/ctf family of mammalian transcription factors and we compared the results obtained from a genome-wide study of its binding sites with chromatin structure assays, gene expression microarray data, and in silico binding site predictions. we found that nfi/ctf family members preferentially bind their dna target sites when they are located around transcription start sites when compared to control datasets generated from the random subsampling of the complete set of nfi binding sites. taken together, we illustrate approaches that can be taken to analyze large genomic data, and provide evidence that nfi family members may act in conjunction with specific chromatin modifications to activate gene expression. however, inferring valuable biological information from the comparisons of very large genomic datasets still represents an enormous challenge. this was documented for one member of the family, which revealed as a more potent gene activator than repressor in global gene expression analysis. nfi proteins preferably associate with the upstream regions of genes that are highly expressed and that are enriched in active chromatin modifications such as h3k4me and h3k36me we postulate that this is a causal association and that nfi proteins mainly act as activators of transcription.
 multiplex experimental assays coupled to computational predictions are being increasingly employed for the simultaneous analysis of many specimens at the genome scale, which quickly generates very large amounts of data. however, inferring valuable biological information from the comparisons of very large genomic datasets still represents an enormous challenge.  as a study model, we chose the nfi/ctf family of mammalian transcription factors and we compared the results obtained from a genome-wide study of its binding sites with chromatin structure assays, gene expression microarray data, and in silico binding site predictions. we found that nfi/ctf family members preferentially bind their dna target sites when they are located around transcription start sites when compared to control datasets generated from the random subsampling of the complete set of nfi binding sites. nfi proteins preferably associate with the upstream regions of genes that are highly expressed and that are enriched in active chromatin modifications such as h3k4me and h3k36me we postulate that this is a causal association and that nfi proteins mainly act as activators of transcription. this was documented for one member of the family, which revealed as a more potent gene activator than repressor in global gene expression analysis. interestingly, we also discovered the association of nfi with the tri-methylation of lysine of histone h a chromatin marker previously associated with the protection against silencing of telomeric genes by nfi.  taken together, we illustrate approaches that can be taken to analyze large genomic data, and provide evidence that nfi family members may act in conjunction with specific chromatin modifications to activate gene expression. 