Skip to main content

Comparative analysis of chromatin landscape in regulatory regions of human housekeeping and tissue specific genes

Abstract

Background

Global regulatory mechanisms involving chromatin assembly and remodelling in the promoter regions of genes is implicated in eukaryotic transcription control especially for genes subjected to spatial and temporal regulation. The potential to utilise global regulatory mechanisms for controlling gene expression might depend upon the architecture of the chromatin in and around the gene. In-silico analysis can yield important insights into this aspect, facilitating comparison of two or more classes of genes comprising of a large number of genes within each group.

Results

In the present study, we carried out a comparative analysis of chromatin characteristics in terms of the scaffold/matrix attachment regions, nucleosome formation potential and the occurrence of repetitive sequences, in the upstream regulatory regions of housekeeping and tissue specific genes. Our data show that putative scaffold/matrix attachment regions are more abundant and nucleosome formation potential is higher in the 5' regions of tissue specific genes as compared to the housekeeping genes.

Conclusion

The differences in the chromatin features between the two groups of genes indicate the involvement of chromatin organisation in the control of gene expression. The presence of global regulatory mechanisms mediated through chromatin organisation can decrease the burden of invoking gene specific regulators for maintenance of the active/silenced state of gene expression. This could partially explain the lower number of genes estimated in the human genome.

Background

Eukaryotic gene transcription is largely known to be orchestrated by protein factors like activators, co-activators and co-repressors [1]. However, nucleosomal organisation, non-passive structural scaffolds and global structure of chromatin are increasingly being recognised as major players in the regulation of gene expression. The ability of sequences to position nucleosomes and to be anchored to the nuclear matrix to provide a spatial context for regulation of expression are measurable parameters that may influence the interactions with transcription machinery [2, 3]. This level of regulation may be distinctly different for genes whose expression is constitutive in comparison to genes that exhibit tissue specific expression. The latter would demand an open chromatin configuration in certain tissues and repressive organisation in others. In this study, we examined whether the potential to utilise global regulatory mechanisms to control gene expression through chromatin organisation varies between housekeeping and tissue specific genes (Hkg and Tsg respectively) by virtue of their organisation. An in-silico comparison of chromatin related organisational differences in the 5' and 3' regulatory regions of housekeeping and tissue specific genes was carried out to shed light in this direction.

Results and discussion

Chromatin landscape of a region plays a major role in determining and modulating the expression status of its neighbouring genes [4]. The role played by chromatin in the 5' regulatory regions of genes in transcriptional regulation has been extensively studied [5, 6]. In the present study, we have taken 2 distinct sets of genes differing predominantly in their spatial expression aspect, namely, housekeeping and tissue specific, to understand the various attributes of the regulatory role played by chromatin organisation in the 5' region.

Analysis of scaffold/matrix associated sequences

Scaffold/matrix attachment regions (S/MARs) are defined as sequences, which can attach themselves to the nuclear matrix and hence help in the formation of independent chromatin loops [7]. Transcriptional regulation of gene expression is known to involve formation of dynamic chromatin loops mediated by S/MAR attachment to the nuclear matrix [3]. The attachment of a DNA sequence to the matrix will place the neighbouring genes in proximity of the transcription factors. The abundance of S/MARs in the 5' cis-regulatory regions of genes further demonstrates their role in transcriptional regulation [8]. We have analysed the predicted S/MAR sites in the 5' and 3' flanking regions of human Hkg and Tsg (Table 1). We used MAR Finder (new version) and ChrClass programs for predicting S/MAR binding sites in the sequences (Table 1). Glazko et al have classified 5' flanking regions up to 1500 bp of human tissue specific genes as an out-group, assuming that these regions have no significant association with S/MAR binding [7]. On the contrary, our study reveals that S/MAR binding sequences are enriched in 5' regulatory regions of Tsg in comparison to the Hkg. The common predictions of both the programs were taken for the analysis. This data indicates a significant enrichment of S/MAR binding sequences in the 5' flanking regions of Tsg and depletion of S/MARs in the 3' Hkg regions as compared to Tsg. Chi-square test was applied for both 5' and 3' region S/MAR predictions of Hkg and Tsg, to ascertain whether the distributions are significantly different. The chi-square value of 11.37 (df = 1) and P-value ≤ 0.001 obtained for the distribution of S/MARs in 5' regions of Hkg and Tsg indicate a significant difference in the distribution of S/MAR elements between the two sets. Similarly, for the distribution of S/MARs in 3' regions of Hkg and Tsg the chi-square value of 5.033 (df = 1) and P-value of ≤ 0.025 show that the Hkg 3' regions are significantly depleted of S/MARs as compared to Tsg.

Table 1 Distribution of putative S/MARs in housekeeping and tissue specific genes.

The observation that the 5' regulatory regions of Hkg are less enriched in S/MARs in comparison with Tsg might be related to the distribution of housekeeping genes in the genome. Housekeeping genes cluster in chromosomes and therefore, they often would be present in distinct chromatin domains along with housekeeping genes that have a co-ordinated expression [9, 10]. The data showing preferential absence of S/MARs in the 3' regions in Hkg further lend support to this hypothesis. On the other hand, tissue specific genes are known to be dispersed in gene dense as well as heterochromatic regions [9, 11]. It may be necessary for them to shield themselves against the effects of positive and negative cis-acting elements of adjacent regions in order to maintain tissue specific expression profile. In this context, the boundary elements or the insulator model has been proposed earlier [11]. S/MARs function as boundary elements and their co-localisation with insulators such as the Drosophila gypsy element is also reported [12, 13]. They also function as boundary elements in in vitro systems by shielding away the position effect [14]. Some earlier reports have suggested a role for S/MARs in maintaining tissue specific gene expression [15]. More recently, the 5'-HS4 chicken-globin insulator is known to have a CTCF protein binding dependent matrix association [16]. Hence, the over representation of S/MARs seen in Tsg set might possibly be associated with a boundary element function.

Our results on the prediction performance of the programs have been quite different from the previous reports [7]. We find that MAR Finder (an under predictor) predicts more number of S/MAR regions in our dataset in comparison to ChrClass program (an over predictor) [7]. This may be attributed to the use of the advanced version of MAR Finder in our study wherein, new parameters/features have been added in the form of the "New MAR Rules" option.

Analysis of nucleosomal organisation

The primary template for local and global changes in the chromatin structure of a chromosome is the nucleosomal unit [4]. Chromatin structure and nucleosomal organisation over the promoter regions play a major role in regulation of expression of downstream gene(s) [6, 17]. The nucleosome distribution would depend upon the occurrence of nucleosome destabilising elements as well as nucleosome forming sequences. We have analysed both these parameters in our study.

Nucleosome destabilising elements

Nucleosome destabilising/excluding elements such as poly (dA.dT) and (CCGNN)n in promoter regions have been implicated in maintaining constitutive gene expression [18–21]. At the functional level, it is known that poly (dA.dT) elements increase the accessibility of promoters of HIS3, URA3 and Ilv1 in yeast to the cognate transcription factor [18]. With the increasing length of poly (dA.dT) repeat, the availability of the sequences to transcription factors improves and similarly, with increasing lengths, the propensity to exclude nucleosomes increases for (CCGNN)n sequence motif as demonstrated in yeast and mammalian systems [19–21]. It has been demonstrated that (CCGNN)n sequences promote meiotic recombination and activated HIS4 expression by generating open chromatin [22].

We hypothesised that the differential distribution of nucleosome exclusion elements might be one of the mechanisms involved in maintaining distinct nucleosomal organisation of the housekeeping and tissue specific genes. The frequency of pure poly (dA.dT) stretches >10 bp and (CCGNN)2–5 in the 2000 bp 5' cis-regulatory regions of human Hkg and Tsg(s) were analysed. A significant enrichment of poly (dA.dT) elements in the upstream regions of Hkg is seen in comparison to Tsg (Table 2). The t-test for the difference in distribution of poly (dA.dT) stretches (>10 bp) between Hkg and Tsg show significant P-values in the different lengths of the stretches examined.

Table 2 Distribution of poly (dA.dT) repeats of various lengths in the 5' upstream regions of housekeeping and tissue specific genes.

In Hkg, 670 repeats of (CCGNN)2–5 were detected as against 430 in Tsg. (CCGNN)2 was the most prevalent repeat unit and uninterrupted repeat units (>5 mers) were not found in the sequence sets. Although shorter repeat units (2–5 mers) have not been studied for nucleosome exclusion, they might play a role in destabilising the histone octamer [20]. Further, many of them form a part of longer interrupted stretches. The t-test for difference in distribution of (CCGNN)2–5 between Hkg and Tsg shows a significant P-value of 1.71E-06.

Nucleosome formation potential scores and expression level of genes

Using Recon, Levitsky et al (2001) have examined the nucleosome formation potential of 3 classes of human genes namely, Hkg, Tsg and widely expressed genes that differ in their spatial expression status [2]. Their report, based on a small sample size of around 200 genes shows the difference in the nucleosome formation potential between these 3 classes of genes in the upstream 50 bp from the transcription start site. In this study, we examined the nucleosome formation potential values in upstream 2000 bp of 5' regions of Hkg and Tsg and their correlation with gene expression levels with the complete set of 1083 genes.

The Tsg and Hkg sequences show a considerable difference in their nucleosome formation potential scores over an extended upstream region of 2000 bp (Figures 1 and 2). The Tsg region is enriched in nucleosome formation potential scores (peak at 1) in all upstream positions analysed (till 2000 bp). For Hkg, the distribution seems to be shifted towards the negative scores at 400 bp region and this shift diminishes gradually as we move further upstream to finally peak at 1 in 2000 bp upstream region (Figure 1). t-test was applied to ascertain the difference in distribution of Recon scores between Hkg and Tsg (Table 3). The resultant P-values in various intervals of relevance (0.8 to 1, 1 to 1.2, -0.8 to -1 and -1 to -1.2) reflect that the scores in the upstream 400 bp from the gene start site show the maximum difference in all the intervals and at 2000 bp, the difference gradually fades away in intervals 0.8 to 1 and 1 to 1.2 (Table 3).

Figure 1
figure 1

Nucleosome formation potential score distributions for 5' regions of housekeeping and tissue specific genes. The 5' sequences of human housekeeping and tissue specific genes were analysed by Recon for distribution of nucleosome formation potential scores. Frequency distribution histograms were plotted for scores in various intervals (range -3.2 to +3.2). (A) and (B) show the distribution of nucleosome formation potential scores at 400 and 2000 bp upstream from the gene start site respectively. Nucleosomal density is significantly lower for housekeeping genes as compared to tissue specific ones, in regions close to the gene start site.

Figure 2
figure 2

Nucleosome formation potential score distributions for 5' regions at different positions from the gene start site in housekeeping and tissue specific genes. The 5' regions of 800, 1200 and 1600 bp from the gene start site of housekeeping and tissue specific genes were taken for the analysis. Frequency distribution histograms were plotted for Recon scores in various intervals (range -3.2 to +3.2). (A), (B) and (C) show the distribution of nucleosome formation potential scores at 800, 1200 and 1600 bp upstream from the gene start site respectively. As we move upstream from the gene start site, the difference in the nucleosome formation potentials between housekeeping and tissue specific genes gradually fades away.

Table 3 t-test P-values for the difference in the distribution of nucleosome formation potential scores between housekeeping and tissue specific genes.

A correlation analysis between nucleosome formation potential and expression levels was carried out considering the Recon scores at upstream 400bp region, where the P-values reflect the largest difference and the log10 values of expression levels were taken as inputs ["see Additional file 1"]. Initially, we analysed the gross dependence of total expression levels on nucleosome potential in the upstream regions of the two sets of genes (Table 4). In all the four intervals, no correlation is seen, indicating that chromatin plays an insignificant role in global modulation of levels of expression in these two sets of genes. These results are similar to that observed in case of Saccharomyces cerevisae whole genome analysis (unpublished results).

Table 4 Correlation coefficients of total expression levels (log10) with nucleosome formation potential scores in housekeeping (Hkg) and tissue specific genes (Tsg).

Further, we refined the analysis to examine the correlation, if any, between nucleosome formation potential in upstream regions and extreme expression levels of genes. The Hkg and Tsg groups were further categorised separately into high and low expression level groups as described under "Methods" section and their correlation with the nucleosome formation potential was analysed (Table 5). The high and low expression genes of Hkg show a low negative correlation with scores in intervals 0.8 to 1.0 and 1 to 1.2 and a low positive correlation with scores in intervals -1.2 to -1 and -1.0 to -0.8. In Tsg, except in one interval, there was no valid correlation seen. This solitary value was not considered since the correlation coefficients in other intervals didn't reflect this trend.

Table 5 Comparison of the level of correlation between nucleosome formation potential scores and contrasting expression levels of genes.

Our data restates that chromatin in 5' region plays a major role in determining the ubiquitous or restricted tissue expression of a gene as shown by Levitsky et al (2001) [2]. The abundance of nucleosome exclusion elements in Hkg 5' regions and the low Recon scores reflect their poor preference for nucleosome assembly. The expression analysis suggests that although chromatin plays a role in bringing about extreme variations of gene expression levels in certain classes of genes such as the housekeeping genes, the relation is not linearly correlated with the total, wider range of expression levels. It is possible that nucleosomes might be involved in fine-tuning of expression levels that may escape our attention, since the difference in the range of expression considered is fairly large. The difference detected in nucleosome formation potential between the two sets might reflect the accessibility to basal transcription factors for Hkg and gene/tissue specific transcription factors for Tsg, considering the difference in spatial and temporal expression patterns of the two groups.

Analysis of repetitive sequences

Repetitive sequences are implicated in chromatin organisation and heterochromatinisation [23–25]. They are differentially enriched in various functional categories of genes and are predicted to play an important role in gene regulation [24, 26]. We analysed the distribution of various repeat classes in the 5' regions of Hkg and Tsg using RepeatMasker software. The total repeat content in Hkg regions is seen to be more than in Tsg regions. As reported earlier, our data shows enrichment of SINES (Alu) in comparison to other classes of repetitive sequences in both the sets [24]. Further, the 5' sequences of Hkg are more enriched in Alu sequences in comparison to those of Tsg regions (Table 6). The difference in the distribution of Alu repeats in the two classes of sequences was determined by applying t-test for the number of repeats and the repeat content in terms of length in base pairs in each sequence set (Table 7). The low total repeat content seen in Tsg upstream regions lends support to the hypothesis that condensed chromatin disfavours transposable element insertions in comparison to open chromatin (Hkg promoters)[27].

Table 6 The distribution of Alu repeats in 5' upstream regions of housekeeping (Hkg) and tissue specific genes (Tsg) is represented in terms of the number of copies and basepairs covered by Alu repeats.
Table 7 t-test P-values for the difference in the distribution of Alu repeats in 5' upstream regions of housekeeping and tissue specific genes.

Genes with high expression levels are clustered in genomic regions known as ridges. These gene rich regions also have high (G+C) content, SINES and genes with short introns [9]. Eisenberg and Levanon [28] have reported the presence of significantly shorter introns and an overall compact gene structure in Hkg as compared to non-Hkg [28]. We have used the gene list provided by Eisenberg and Levanon [28] for our analysis. The enrichment of SINES in the 5' regions of Hkg suggests that Hkg might be localised in the ridge regions of the genome. More recently, it has been suggested that the contrasting attributes of gene compactness, GC content and the length of the intronic and intergenic sequences in Hkg and Tsg might be involved in chromatin mediated regulation for maintaining distinct expression patterns in the gene sets [29]. Recently, Alu elements have been shown to house transcription factor binding sites and the presence of such regulatory elements might influence the chromatin structure and gene expression [30].

The paradigm for regulation of gene expression in human tissues has shifted the focus from involvement of a battery of transcription regulators to global regulatory mechanisms [31]. These mechanisms have also gained significance in the context of the low estimates of gene numbers in the human genome [32]. It is in this framework that we have analysed the chromatin characteristics of two groups of genes, one that needs almost a continuous and ubiquitous expression and another demanding tissue specific regulation. It had been predicted that the nucleosomal density in a chromatin domain and the buffering of supercoiling waves by repetitive DNA will play a major role in establishing coordinated gene regulation in a domain in the context of the relevance of maintenance of repetitive sequences during evolution [[25, 33], and [34]]. A recent report also infers the role of chromatin-mediated mechanisms in the differential gene expression patterns seen in housekeeping and tissue specific genes [29]. Our data and analyses lend support to these hypotheses (Figure 3). Another recent report, which addresses the chromatin architecture of the human genome, provides experimental evidence that open chromatin correlates with high gene density regions but not with gene expression [35]. This data further supports our in-silico observations and strengthens the domain concept for concerted expression of clustered genes. The role of nucleosome formation potential is apparent from the present analysis in both the housekeeping genes as well as tissue specific genes but with an opposing correlation. Housekeeping genes apparently discourage nucleosome formation to match their expression profile in space and time by ensuring accessibility to transcription machinery. In addition, they also show a significant enrichment in poly (dA.dT) stretches, which are known to destabilise nucleosomes. On the other hand, the tissue specific genes show higher scores for nucleosome formation potential through which they perhaps provide selective accessibility to the transcriptional machinery. Further, our analysis suggests that tissue specific genes resort to additional global regulatory features such as matrix association, which would facilitate maintenance of functionally distinct domains to insulate themselves from both silencing and activating regulatory influence of adjacent domains. The differential distribution of repetitive sequences in housekeeping and tissue specific genes might also play an important role in maintaining distinct chromatin landscape over these regions.

Figure 3
figure 3

A model for chromatin landscape in 5' regions of tissue specific and housekeeping genes. (A) depicts the repressive role of chromatin in maintaining tissue specific gene expression profiles in a chromosome. The chromatin organisation in the 5' regions of Tsg1 and Tsg2, two different tissue specific genes dispersed in the chromosome is shown. Nucleosome formation potentials and S/MARs – the boundary elements, are enriched in their upstream regions and might play a major role in facilitating tissue specific expression. This is likely to be a local effect since neighbouring genes might have a different expression pattern. (B) depicts the chromatin organisation in the 5' regions of Hkg1, Hkg2 and Hkg3, three housekeeping genes clustered in the chromosome. The presence of low nucleosome formation potential regions and enrichment of nucleosome destabilising elements ensure an open chromatin configuration in this domain. As Hkg generally cluster together, they are depleted in S/MARs relative to tissue specific genes as shown in the present analysis by the significant absence of predicted S/MARs in both 5' and 3' regions of housekeeping genes as compared to tissue specific genes.

Conclusion

We have demonstrated that the regulatory regions of housekeeping and tissue specific genes have differential chromatin architecture with respect to S/MAR binding, nucleosome positioning potential and repetitive sequences. This has potential implications for regulation of gene expression in eukaryotic genomes.

Methods

In this study, the 5' and 3' flanking regions of genes were analysed for various attributes of chromatin organisation. The list of human housekeeping genes (Hkg) was retrieved from http://www.compugen.co.il/supp_info/Housekeeping_genes.html[28, 36]. 532 genes have been categorised as housekeeping because of their ubiquitous and high expression levels in 47 tissues. The list and expression levels of the human tissue specific genes were obtained from Eli Eisenberg (personal communication). 566 genes expressed in only a single tissue were taken as tissue specific genes (Tsg) and analysed. We could unambiguously retrieve sequences of 525 Hkg and 558 Tsg from human genome build 33 (NCBI). Approximately, 2000 bp of the 5' and 3' regions from each of these genes were taken for analysis.

Scaffold/matrix associated regions (S/MAR) analysis

MAR Finder was used for prediction of S/MAR regions [37, 38]. All the default options and the "New MAR Rules" were selected for predicting S/MARs. ChrClass program was used for S/MAR prediction [39, 40].

Nucleosome organisation and gene expression correlation analysis

The upstream regions (2000 bp) were scanned for nucleosome exclusion elements [18, 20] – poly (dA.dT) pure stretches of >10 bp length and [5' (CCGNN) 3']2–5 using in-house programs. Recon was used for evaluating nucleosome formation potential in the sequences [2, 41]. The score outputs of the 5' regions were categorised in frequency intervals of 0.2 with a range from -3.2 to +3.2. The Recon scores around +1 and -1 imply strong nucleosome formation and exclusion potentials respectively. The scores in the four intervals of relevance (0.8 to 1, 1 to 1.2, -0.8 to -1 and -1 to -1.2) were taken for all the analyses. Since the promoter region information was not retrieved for these genes, the 2000 bp upstream region from the gene start site was split into 400, 800, 1200 & 1600 bp and analysed.

The Recon scores at 400 bp were used to draw correlation between the nucleosome formation potential and expression levels in the two sets of genes. In each sequence set, genes with expression levels <500 and >5000 affymetrix expression units were classified as low and high expression genes respectively. We considered a minimum ten fold difference in the expression levels of genes as a relevant criterion for classifying them as high and low expression genes. In Hkg, this criterion yielded 33 low expression and 35 high expression genes. In Tsg, we categorised 416 low expression genes and 24 high expression genes.

Repetitive sequence analysis

RepeatMasker version: 20040306-web was used to calculate the repeat content in 2000 bp upstream sequences of the two groups of genes [42].

Abbreviations

S/MAR:

scaffold/matrix attachment regions

Hkg:

housekeeping genes

Tsg:

tissue specific genes

df:

degree of freedom

References

  1. Lemon B, Tjian R: Orchestrated response: a symphony of transcription factors for gene control. Genes Dev 2000, 14: 2551–2569. 10.1101/gad.831000

    Article  CAS  PubMed  Google Scholar 

  2. Levitsky VG, Podkolodnaya OA, Kolchanov NA, Podkolodny NL: Nucleosome formation potential of eukaryotic DNA: calculation and promoters analysis. Bioinformatics 2001, 17: 998–1010. 10.1093/bioinformatics/17.11.998

    Article  CAS  PubMed  Google Scholar 

  3. Bode J, Benham C, Knopp A, Mielke C: Transcriptional augmentation: modulation of gene expression by scaffold/matrix-attached regions (S/MAR elements). Crit Rev Eukaryot Gene Expr 2000, 10: 73–90.

    Article  CAS  PubMed  Google Scholar 

  4. Grewal SI, Moazed D: Heterochromatin and epigenetic control of gene expression. Science 2003, 301: 798–802. 10.1126/science.1086887

    Article  CAS  PubMed  Google Scholar 

  5. Boeger H, Griesenbeck J, Strattan JS, Kornberg RD: Nucleosomes unfold completely at a transcriptionally active promoter. Mol Cell 2003, 11: 1587–1598. 10.1016/S1097-2765(03)00231-4

    Article  CAS  PubMed  Google Scholar 

  6. Wolffe AP: Transcriptional activation. Switched-on chromatin. Curr Biol 1994, 4: 525–528. 10.1016/S0960-9822(00)00114-7

    Article  CAS  PubMed  Google Scholar 

  7. Glazko GV, Rogozin IB, Glazkov MV: Comparative study and prediction of DNA fragments associated with various elements of the nuclear matrix. Biochim Biophys Acta 2001, 1517: 351–364.

    Article  CAS  PubMed  Google Scholar 

  8. Glazko GV, Koonin EV, Rogozin IB, Shabalina SA: A significant fraction of conserved noncoding DNA in human and mouse consists of predicted matrix attachment regions. Trends Genet 2003, 19: 119–124. 10.1016/S0168-9525(03)00016-7

    Article  CAS  PubMed  Google Scholar 

  9. 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.1649303

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Lercher MJ, Urrutia AO, Hurst LD: Clustering of housekeeping genes provides a unified model of gene order in the human genome. Nat Genet 2002, 31: 180–183. 10.1038/ng887

    Article  CAS  PubMed  Google Scholar 

  11. de Laat W, Grosveld F: Spatial organisation of gene expression: the active chromatin hub. Chromosome Res 2003, 11: 447–459. 10.1023/A:1024922626726

    Article  CAS  PubMed  Google Scholar 

  12. Byrd K, Corces VG: Visualization of chromatin domains created by the gypsy insulator of Drosophila. J Cell Biol 2003, 162: 565–574. 10.1083/jcb.200305013

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Nabirochkin S, Ossokina M, Heidmann T: A nuclear matrix/scaffold attachment region co-localizes with the gypsy retrotransposon insulator sequence. J Biol Chem 1998, 273: 2473–2479. 10.1074/jbc.273.4.2473

    Article  CAS  PubMed  Google Scholar 

  14. Kim JM, Kim JS, Park DH, Kang HS, Yoon J, Baek K, Yoon Y: Improved recombinant gene expression in CHO cells using matrix attachment regions. J Biotechnol 2004, 107: 95–105. 10.1016/j.jbiotec.2003.09.015

    Article  CAS  PubMed  Google Scholar 

  15. Bonifer C, Yannoutsos N, Kruger G, Grosveld F, Sippel AE: Dissection of the locus control function located on the chicken lysozyme gene domain in transgenic mice. Nucleic Acids Res 1994, 22: 4202–4210.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Yusufzai TM, Felsenfeld G: The 5'-HS4 chicken beta-globin insulator is a CTCF-dependent nuclear matrix-associated element. Proc Natl Acad Sci U S A 2004, 101: 8620–8624. 10.1073/pnas.0402938101

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Khorasanizadeh S: The nucleosome: from genomic organisation to genomic regulation. Cell 2004, 116: 259–272. 10.1016/S0092-8674(04)00044-3

    Article  CAS  PubMed  Google Scholar 

  18. Suter B, Schnappauf G, Thoma F: Poly(dA.dT) sequences exist as rigid DNA structures in nucleosome-free yeast promoters in vivo. Nucleic Acids Res 2000, 28: 4083–4089. 10.1093/nar/28.21.4083

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Koch KA, Thiele DJ: Functional analysis of a homopolymeric (dA-dT) element that provides nucleosomal access to yeast and mammalian transcription factors. J Biol Chem 1999, 274: 23752–23760. 10.1074/jbc.274.34.23752

    Article  CAS  PubMed  Google Scholar 

  20. Wang YH, Griffith JD: The [(G/C)3NN]n motif: a common DNA repeat that excludes nucleosomes. Proc Natl Acad Sci U S A 1996, 93: 8863–8867. 10.1073/pnas.93.17.8863

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Iyer V, Struhl K: Poly(dA:dT), a ubiquitous promoter element that stimulates transcription via its intrinsic DNA structure. Embo J 1995, 14: 2570–2579.

    PubMed Central  CAS  PubMed  Google Scholar 

  22. Kirkpatrick DT, Wang YH, Dominska M, Griffith JD, Petes TD: Control of meiotic recombination and gene expression in yeast by a simple repetitive DNA sequence that excludes nucleosomes. Mol Cell Biol 1999, 19: 7661–7671.

    PubMed Central  CAS  PubMed  Google Scholar 

  23. Grover D, Mukerji M, Bhatnagar P, Kannan K, Brahmachari SK: Alu repeat analysis in the complete human genome: trends and variations with respect to genomic composition. Bioinformatics 2004, 20: 813–817. 10.1093/bioinformatics/bth005

    Article  CAS  PubMed  Google Scholar 

  24. Jordan IK, Rogozin IB, Glazko GV, Koonin EV: Origin of a substantial fraction of human regulatory sequences from transposable elements. Trends Genet 2003, 19: 68–72. 10.1016/S0168-9525(02)00006-9

    Article  CAS  PubMed  Google Scholar 

  25. Brahmachari SK, Meera G, Sarkar PS, Balagurumoorthy P, Tripathi J, Raghavan S, Shaligram U, Pataskar S: Simple repetitive sequences in the genome: structure and functional significance. Electrophoresis 1995, 16: 1705–1714. 10.1002/elps.11501601283

    Article  CAS  PubMed  Google Scholar 

  26. Grover D, Majumder PP, Rao CB, Brahmachari SK, Mukerji M: Nonrandom distribution of alu elements in genes of various functional categories: insight from analysis of human chromosomes 21 and 22. Mol Biol Evol 2003, 20: 1420–1424. 10.1093/molbev/msg153

    Article  CAS  PubMed  Google Scholar 

  27. Vijaya S, Steffen DL, Robinson HL: Acceptor sites for retroviral integrations map near DNase I-hypersensitive sites in chromatin. J Virol 1986, 60: 683–692.

    PubMed Central  CAS  PubMed  Google Scholar 

  28. Eisenberg E, Levanon EY: Human housekeeping genes are compact. Trends Genet 2003, 19: 362–365. 10.1016/S0168-9525(03)00140-9

    Article  CAS  PubMed  Google Scholar 

  29. Vinogradov AE: Compactness of human housekeeping genes: selection for economy or genomic design? Trends Genet 2004, 20: 248–253. 10.1016/j.tig.2004.03.006

    Article  CAS  PubMed  Google Scholar 

  30. Oei SL, Babich VS, Kazakov VI, Usmanova NM, Kropotov AV, Tomilin NV: Clusters of regulatory signals for RNA polymerase II transcription associated with Alu family repeats and CpG islands in human promoters. Genomics 2004, 83: 873–882. 10.1016/j.ygeno.2003.11.001

    Article  CAS  PubMed  Google Scholar 

  31. Nemeth A, Langst G: Chromatin higher order structure: opening up chromatin for transcription. Brief Funct Genomic Proteomic 2004, 2: 334–343.

    Article  CAS  PubMed  Google Scholar 

  32. Pennisi E: Human genome. A low number wins the GeneSweep Pool. Science 2003, 300: 1484. 10.1126/science.300.5625.1484b

    Article  CAS  PubMed  Google Scholar 

  33. Brahmachari SK, Ramesh N, Shouche YS, Mishra RK, Bagga R, Meera G: Unusual DNA Structures: Sequence Requirements and Role in Transcriptional Control. In Structure and Methods. Volume 2. Edited by: Sarma RH, Sarma MH. New York: Adenine Press; 1990:33–49.

    Google Scholar 

  34. Conrad M, Brahmachari SK, Sasisekharan V: DNA structural variability as a factor in gene expression and evolution. Biosystems 1986, 19: 123–126. 10.1016/0303-2647(86)90024-9

    Article  CAS  PubMed  Google Scholar 

  35. Gilbert N, Boyle S, Fiegler H, Woodfine K, Carter NP, Bickmore WA: Chromatin architecture of the human genome: gene-rich domains are enriched in open chromatin fibers. Cell 2004, 118: 555–566. 10.1016/j.cell.2004.08.011

    Article  CAS  PubMed  Google Scholar 

  36. List of housekeeping genes[http://www.compugen.co.il/supp_info/Housekeeping_genes.html]

  37. 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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. MAR Finder[http://futuresoft.org/MAR-Wiz/]

  39. Rogozin IB, Glazko GV, Glazkov MV: Computer prediction of sites associated with various elements of the nuclear matrix. Brief Bioinform 2000, 1: 33–44.

    Article  CAS  PubMed  Google Scholar 

  40. ChrClass[http://ftp.bionet.nsc.ru/pub/biology/chrclass/chrclass.zip]

  41. Recon[http://wwwmgs.bionet.nsc.ru/mgs/programs/recon/]

  42. RepeatMasker[http://www.repeatmasker.org/]

Download references

Acknowledgements

MG acknowledges the financial support provided by Council for Scientific and Industrial Research (CSIR), India. SKB thanks Council for Scientific and Industrial Research, India and VB thanks Indian Council for Medical Research (ICMR), for financial assistance through a grant. The authors wish to acknowledge Dr. Beena Pillai, Dr. Rakesh Sharma, Dr. Neeraj Pandey and Dr. Mitali Mukerji for their valuable discussions and suggestions. We would also like to acknowledge Samira for careful checking and helping with the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Samir K Brahmachari.

Additional information

Authors' contributions

MG, VB and SKB contributed to the study design, analyses and in drafting the manuscript. MG, PS, SKDS, DD, KK and GP were involved in the data retrieval and analyses. All authors read and approved the final manuscript.

Electronic supplementary material

12859_2004_451_MOESM1_ESM.xls

Additional File 1: gene list and expression levels of housekeeping genes, gene list and expression levels of tissue specific genes. 'supplementaryfile1.xls' contains the list of housekeeping and tissue specific genes and their expression levels used for the analysis. (XLS 117 KB)

Authors’ original submitted files for images

Rights and permissions

Reprints and permissions

About this article

Cite this article

Ganapathi, M., Srivastava, P., Sutar, S.K.D. et al. Comparative analysis of chromatin landscape in regulatory regions of human housekeeping and tissue specific genes. BMC Bioinformatics 6, 126 (2005). https://doi.org/10.1186/1471-2105-6-126

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-6-126

Keywords