Skip to main content

Spatial organization of the chicken beta-globin gene domain in erythroid cells of embryonic and adult lineages

Abstract

Background

The β-globin gene domains of vertebrate animals constitute popular models for studying the regulation of eukaryotic gene transcription. It has previously been shown that in the mouse the developmental switching of globin gene expression correlates with the reconfiguration of an active chromatin hub (ACH), a complex of promoters of transcribed genes with distant regulatory elements. Although it is likely that observations made in the mouse β-globin gene domain are also relevant for this locus in other species, the validity of this supposition still lacks direct experimental evidence. Here, we have studied the spatial organization of the chicken β-globin gene domain. This domain is of particular interest because it represents the perfect example of the so-called ‘strong’ tissue-specific gene domain flanked by insulators, which delimit the area of preferential sensitivity to DNase I in erythroid cells.

Results

Using chromosome conformation capture (3C), we have compared the spatial configuration of the β-globin gene domain in chicken red blood cells (RBCs) expressing embryonic (3-day-old RBCs) and adult (9-day-old RBCs) β-globin genes. In contrast to observations made in the mouse model, we found that in the chicken, the early embryonic β-globin gene, Ε, did not interact with the locus control region in RBCs of embryonic lineage (3-day RBCs), where this gene is actively transcribed. In contrast to the mouse model, a strong interaction of the promoter of another embryonic β-globin gene, ρ, with the promoter of the adult β-globin gene, βA, was observed in RBCs from both 3-day and 9-day chicken embryos. Finally, we have demonstrated that insulators flanking the chicken β-globin gene domain from the upstream and from the downstream interact with each other, which places the area characterized by lineage-specific sensitivity to DNase I in a separate chromatin loop.

Conclusions

Taken together, our results strongly support the ACH model but show that within a domain of tissue-specific genes, the active status of a promoter does not necessarily correlate with the recruitment of this promoter to the ACH.

Background

The domain of chicken beta-globin genes is located on chromosome 1 and has a length of approximately 33 kb. It includes a cluster of four beta-globin genes: ρ (HBG1), βH (HBE1), βA (HBG2) and Ε (HBE) and several distant regulatory regions, which are marked with sites of hypersensitivity to DNase I (HS) and are necessary for the regulation of transcription, replication and chromatin status of the domain[1, 2]. The locus control region of the domain (LCR) is located upstream of the embryonic β-globin gene ρ and is composed of three blocks co-localizing with the erythroid cell-specific HSs 1 to 3[3, 4]. A constitutive HS4 located upstream of the LCR marks the position of the well-studied CTCF-dependent (CTCF ‐ СССTC-binding protein factor) insulator[59]. The CTCF-dependent enhancer-blocking element is also located at the downstream end of the domain in the area marked by the so-called 3′ constitutive HS[10, 11]. In addition to the LCR, the domain also possesses an internal enhancer (β/Ε-enhancer), which is located between the genes βA and Ε and is believed to contribute to the control of the activity of both genes[1214]. The partitioning of tasks between the LCR and the β/Ε-enhancer is not quite clear. The results of experiments with transgenic mice carrying the chicken β-globin gene domain[4] suggested that the expression of one of the embryonic genes (ρ) is controlled by the LCR, while the expression of the other embryonic gene (Ε) is controlled by the β/Ε-enhancer. In addition, both regulatory elements appeared to be necessary for the correct developmental expression of the adult β-globin genes[3]. The domain of the chicken β-globin gene represents one of the best-studied genomic areas in higher eukaryotes[15]. This domain possesses all of the characteristic features of the so-called closed or strong chromatin domains[16, 17], which change the pattern of sensitivity to DNase I in a cell lineage-specific fashion. Detailed studies of the distribution of different histone modifications along the domain have contributed much to the present knowledge about the histone code[18, 19]. Nevertheless, the mechanisms of transcriptional switching of chicken β-globin gene expression remain poorly understood[15]. The developmental switching of globin gene transcription is typical for all vertebrate organisms[2]. In chickens, the first molecules of hemaglobin appear in embryonic erythroblasts 32 h after fertilization[20]. The beta-chains of the hemaglobin represent the products of genes ρ and Ε, which are transcribed from the 2nd to the 5th day of the embryonic development. After the 5th day a new population of erythroid cells expressing the genes of fetal and adult globins βН and βA appears. These cells gradually replace the embryonic-type erythroblasts that express genes 1ρ and Ε[21]. In erythroblasts of adult lineage, the transcription of embryonic β-type globin genes is repressed.

Figure 1
figure 1

Analysis of β-globin mRNA levels in chicken 3-day and 9-day embryonic red blood cells. The vertical axis shows the relative amounts of β-globin gene primary transcripts, as determined by RT-qPCR analysis of RNA samples prepared from the blood of 3- and 9-day old chicken embryos. (A) Normalization of data to the level of β-actin mRNA. (B) Normalization of data to the amount of cells used for RNA preparation. The amounts of β-globin gene transcripts detected in different RNA samples were normalized to the number of cells used for preparation of these RNA samples. The amount of the most abundant transcript (ρ-transcript in 3-day RBCs) was set as 100 relative units and the other data were normalized to this value. The error bars represent the S.E.M. for two independent experiments.

Recent studies of the spatial configuration of the mouse and human β-globin gene domain suggest that in erythroid cells, all regulatory elements of the domain become assembled into a common activation complex (the ACH) where the promoters of different globin genes are recruited in a developmental stage-specific fashion[2224].

In the present work, we have analyzed the spatial configuration of the chicken β-globin gene domain in circulating red blood cells (RBCs) both before (3-day RBCs) and after (9-day RBCs) the β-globin switching. Using the quantitative chromosome conformation capture assay (3C-qPCR)[25, 26], we have demonstrated that the switching of chicken β-globin gene expression correlates with prominent changes in the spatial configuration of the domain. However, the observed spectrum of changes does not perfectly fit the simplistic model that postulates that only the transcriptionally active genes are recruited to the ACH.

Results

Measurement of β-globin mRNA levels in red blood cells of 3- and 9-day-old chicken embryos

We first estimated the relative levels of primary transcripts of each of the β-globin genes in RBCs (red blood cells) of 3-and 9-day-old chicken embryos using reverse transcription with subsequent TaqMan™ real-time PCR. To exclude the spliced mRNA from this analysis, we used test amplicons that spanned exon-intron junctions. The results obtained (Figure1) were consistent with the previously reported data[1, 21]. In 3-day RBCs, only embryonic Ε and ρ genes were actively transcribed, with the level of ρ transcripts being approximately 4.5-fold higher than the level of Ε transcripts. In 9-day RBCs, transcripts of all four β-globin genes were detected, and as expected, the transcript of the βA gene was the most abundant. The presence of small amounts of Ε and ρ gene transcripts in 9-day RBCs is again consistent with a previous observation that approximately 10% of chicken embryonic RBCs express Ε- and ρ-globin even on the 10th day of development[1].

Figure 2
figure 2

MboI and NlaIII q3C analysis of the chicken β-globin gene domain. (A) Map showing the positions of genes and known regulatory elements of the chicken β-globin gene domain. The open boxes indicate genes, the blue lines represent exons, and the arrows indicate the direction of transcription. The closed green boxes show HSs. Upstream insulator, downstream enhancer-blocking elements and the β/Ε enhancer are shown by orange ellipses. FOLR1 - gene encoding folate receptor, transcribed in earlier erythroid progenitors; HAS - tissue-specific erythroid DNase I-hypersensitive site; HS1-3 - locus control region (LCR) of β-globin gene domain; OR51M1 - open reading frame for odorant receptor. The map was made using nucleotide sequences NW_001471556 and L17432.1 [GenBank:NW_001471556, GenBank:L17432.1]. (B-G) The relative frequencies of crosslinking of the anchor fragments containing the HS2 (B), ρ-promoter (C), upstream portion of gene βA(D), β/Ε-enhancer (E), upstream portion of gene Ε (F) and Ε-promoter (G). (B-F) show results of MboI 3C analysis and (G) shows the results of NlaIII 3C analysis. The dark grey rectangle in each panel indicates the anchor DNA fragment, and the light grey rectangles indicate the test fragments. The white rectangles indicate the restriction fragments that were not analyzed. The borders between neighboring fragments are indicated by dark grey lines. The tailless arrows represent the primers used for q3C analysis. The positions of test fragments on the restriction map are plotted on the x-axis. The cross-linking frequencies are plotted on the y-axis; the highest cross-linking frequency observed was set as 100 relative units (the cross-linking frequency between the fragment containing the ρ-promoter and the fragment containing the βA gene in 3-day RBCs for MboI 3C analysis, and the cross-linking frequency between the fragment containing HS2 and the fragment containing HS1 in 3-day RBCs for NlaIII 3C analysis (data not shown)). The results of 3C analysis for 3-day RBCs, 9-day RBCs and CEF are shown by the blue, red and green lines, respectively. The error bars represent the S.E.M. for three (B-F) or two (G) independent experiments.

To be able to compare the relative levels of transcription of different globin genes in 3-day and 9-day RBCs, the experimentally detected amounts of globin gene primary transcripts were normalized to the experimentally detected amount of β-actin mRNA (Figure1A). However, it is worth mentioning that the total amount of RNA routinely extracted from the same number of 3-day and 9-day RBCs was much higher in the case of 3-day RBCs (approximately 50 μg and 1.7 μg from 10 million 3-day and 9-day RBCs, respectively). Consequently, the amounts of mRNAs of housekeeping genes differ drastically in samples prepared from the same number of 3-day and 9-day RBCs. For this reason, normalization to the amount of cells (3-day and 9-day RBCs) used for the extraction of RNA (Figure1B) appears to be more appropriate for comparison of the absolute levels of globin gene transcription in these cells. Indeed, using this normalization one can see that total amounts of β-globin gene transcripts are comparable in 3-day and 9-day RBCs while the spectrums of transcribed β-globin genes are different.

Interaction of the locus control region and β/Ε enhancer with promoters of β-globin genes in red blood cells of 3-day and 9-day chicken embryos

The quantitative chromosome conformation capture (3C-qPCR) assay was used to investigate the spatial organization of the chicken β-globin gene domain in RBCs of 3- and 9-day-old embryos. It is well known that the accurate analysis and interpretation of raw 3C data is not a trivial task[27], and it can be difficult to distinguish the relevant chromatin interactions from background and cross-linking artifacts. One problem in the interpretation of the 3C data is that neighboring restriction fragments are always ligated to a certain extent due to cross-links between nucleosomal proteins, even if these fragments are not involved in the assembly of a chromatin hub. To discriminate the functionally relevant interactions from the background cross-linking, we performed the 3C analysis of the β-globin gene domain in chicken embryonic fibroblasts (CEF). These cells do not express globin genes. Thus, there are good reasons to expect that the domain in these cells does not adopt any functionally relevant spatial configuration.

In the main set of 3C experiments, we used MboI restriction enzyme digestion to cleave the DNA before the proximity ligation. This enzyme cuts the β-globin gene domain into a number of relatively small fragments, such that the promoters of β-globin genes and regulatory elements were located in different restriction fragments separated by several cleavage sites. It was thus possible to independently analyze all of these regulatory elements. In some cases, however, it was not possible to place the primers and TaqMan™ probes exactly on the regulatory element under study due to the small size of the corresponding restriction fragment or the presence of repetitive sequences. In these cases, the primers and TaqMan™ probes were placed at the adjacent restriction fragments. For example, for βA and Ε genes, the probes were placed at the fragments located just downstream of the promoters.

We have first fixed the anchor at the restriction fragment containing HS2 (a part of β-globin domain LCR). In the 3-day RBCs, this anchor was shown to interact preferentially with the fragments containing the promoter of the embryonic gene ρ, the adult gene βA and the β/Ε-enhancer (Figure2B, blue curve). In addition, an elevated frequency of interaction of the HS2 with the upstream and the downstream insulators of the β-globin gene domain was observed. Interestingly, no interaction of the HS2 with the embryonic gene Ε was detected. As expected, no interaction between any of the functional elements studied was detected in embryonic fibroblasts (Figure2B, green curve). In 9-day RBCs the frequency of interactions of the HS2 with the βA gene increased. In addition, a strong interaction between the HS2 and the promoter of βH gene was established (Figure2B, red curve). The validity of the observations made with the anchor placed at HS2 was confirmed in experiments with the anchor placed on the promoter of the ρ gene (Figure2C), the upstream portion of the βA gene (Figure2D), the upstream portion of the Ε gene (Figure2F) and the β/Ε-enhancer (Figure2E). The anchor placed on the promoter of the ρ gene demonstrated strong interaction frequencies with the HSs 1‐3 of the LCR in 3-day RBCs (Figure2C, blue curve). The frequency of interaction of this anchor with the HS1 remained about the same in 9-day RBCs while the frequencies of interaction with the HSs 2 and 3 decreased (Figure2C, red curve). Unexpectedly, in 3-day RBCs the promoter of the ρ gene demonstrated a strong interaction with the upstream region of the βA gene. The frequency of this interaction decreased in 9-day RBCs but still remained quite prominent (Figure2C, blue and red curves). The interaction between ρ and βA genes was also observed in embryonic fibroblasts (Figure2C, green curve). It may thus reflect some internal feature of the β-globin gene domain spatial organization that is not related to the transcriptional status of the domain. The interaction between the βA and ρ gene promoters observed in the experiment with the anchor placed on the ρ gene promoter was confirmed in the reciprocal experiment with the anchor placed on the upstream part of the βA gene (Figure2D, all curves). No interaction of the ρ gene promoter with the promoter of the minor adult βH gene was observed in 3-day RBC (Figure2C, blue curve). In contrast, in the same cells (3-day RBC) the upstream region of the βA gene demonstrated strong interactions with both the ρ and βH gene promoters (Figure2D, blue curve). As expected, the upstream region of the βA gene demonstrated a strong interaction with the HSs 1 to 3 of the LCR in 9-day RBCs, but not in 3-day RBCs (Figure2D, red and blue curves). A strong interaction of the βA gene with the β/Ε enhancer was also observed in 9-day but not in 3-day RBCs (Figure2D, red and blue curves). Although in this particular case the anchor (the upstream part of the βA gene) and the test-fragment (β/Ε enhancer) are located close to each other, the much higher interaction frequency observed in 9-day RBCs comparing to 3-day RBCs and CEF argues for the existence of functionally relevant interaction between the βA gene and the β/Ε enhancer in 9-day RBCs. This conclusion is supported by the results of a reciprocal experiment with the anchor placed on the β/Ε enhancer. Again, in 9-day RBCs the frequency of interaction of the β/Ε enhancer with the upstream part of the βA gene was much higher in 9-day RBCs than in 3-day RBCs and CEF (Figure2E). The interaction of the β/Ε enhancer with the three HSs of the LCR and the promoter of ρ gene was observed in both 3-day and 9-day RBCs, but not in CEF. On the other hand, interaction of the β/Ε enhancer with the promoter of gene βH was observed only in 9-day RBCs (Figure2E, red curve). No interaction of the β/Ε enhancer with the gene Ε was observed in either of the cell lines studied (Figure2E). In full accord with this observation, the anchor placed on the upstream area of the gene Ε showed no interaction with any of the test fragments studied in either erythroid or non-erythroid cells (Figure2F). Thus it appears that the Ε gene does not interact with any of the remote regulatory elements of the domain. To verify this interesting observation, we performed 3C experiments using digestion with another restriction enzyme, NlaIII. Distribution of NlaIII recognition sites along the β-globin gene domain allows fixation of the anchor exactly on the promoter of the gene Ε. Using this anchor we performed 3C analysis of the β-globin domain configuration in 3-day RBCs and in CEF. Again, no interaction of the Ε gene promoter with any of the remote regulatory elements was observed (Figure2G). As a matter of fact, the frequencies of interaction of the Ε gene promoter with all tested fragments of the β-globin gene domain were very similar in 3-day RBCs that express the Ε gene and in embryonic fibroblasts that do not express globin genes at all.

Interaction between the upstream and the downstream insulators of the β-globin gene domain

Figure3 demonstrates the results of 3C experiments with the anchors placed on either upstream or downstream insulators of the β-globin gene domain. It is evident that in both 3-day and 9-day RBCs, the insulators strongly interact with each other, thus forming a loop harboring the entire β-globin gene domain. In addition, in 9-day RBCs, both insulators demonstrate an interaction with the β/Ε enhancer, while in 3-day RBCs, the upstream insulator interacts with the βA gene. Finally, in both the 3-day and the 9-day RBCs, both insulators show some interaction with the LCR. In embryonic fibroblasts neither the interaction between the upstream and the downstream insulators of the β-globin gene domain nor an interaction of these insulators with any of the internal fragments of the domain was observed (Figure3A/B, green curves).

Figure 3
figure 3

MboI q3C analysis of spatial interactions of β-globin boundary insulators. The relative frequencies of cross-linking of the anchor fragments containing the 5′-insulator (A) and 3′-insulator (B) of the chicken β-globin gene domain. All designations are as described in Figure2.

Discussion

The current model of the enhancer action suggests that enhancers directly interact with target promoters[23, 2830]. As a result, an assembly of regulatory elements, called the active chromatin hub (ACH), forms. One or several enhancers present in the ACH may control the activity of one or several promoters recruited to the ACH[23, 29]. The transcriptional switching of murine β-globin gene expression correlates with a reconfiguration of the ACH[24]. The aim of the present study was to determine whether the developmental switching of chicken β-globin gene transcription depends on the spatial reconfiguration of the β-globin gene domain. The study of the chicken β-globin gene domain is particularly interesting for two reasons. First, this domain is the best example of the so-called closed chromatin domains[17], as the upstream and the downstream insulators perfectly delimit the area of preferential sensitivity to DNase I in erythroid cells[31], which is not the case in mouse and human β-globin gene domains[32, 33]. Second, the chicken β-globin gene domain contains an internal enhancer element in addition to the upstream locus control region[12]. Once again, this is not the case in the mouse or human. Hence, the regulation of chicken β-globin gene expression may differ from regulation of β-globin gene expression in the mouse and human.

The results of the present study showed that in erythroid cells there is a clear interaction of the upstream and the downstream insulators of the chicken β-globin gene domain. As a result, the domain should be organized in a closed loop. Although we have not studied the mechanism of this loop formation, it is very probable that CTCF is involved, as both insulators possess CTCF-binding sites[9, 10] and it was shown that CTCF might mediate an interaction between distant genomic elements[34, 35]. It is also likely that the β-globin loop ends are fixed at the nuclear matrix, as the 5′ insulator of the chicken β-globin gene domain has been reported to interact with the nuclear matrix[36]. This type of spatial organization perfectly explains the fact that in the chicken β-globin gene domain, the area demonstrating the cell-lineage-specific sensitivity to DNase I perfectly co-localizes with the genomic segment flanked by insulators. In agreement with this model, the interaction between the upstream and the downstream insulators of the chicken β-globin gene domain was not observed in non-erythroid cells where the domain is not preferentially sensitive to DNase I[37]. In should be mentioned, however, that in mouse and human the insulators flanking the β-globin gene domain from the upstream and the downstream also interact with each other[22, 24, 34], although the chromatin loop established by this interaction does not correspond to the region of preferential sensitivity to DNase I[32].

The most interesting feature of the spatial organization of the chicken β-globin gene domain is that the embryonic β-globin gene Ε interacts with neither the LCR nor the β/Ε enhancer. This does not exclude the possibility that the β/Ε enhancer located at a short distance upstream of the Ε gene can activate this gene promoter by some tracking mechanism. However, an activation effect of the LCR on Ε gene transcription seems very improbable. In this respect, it is of interest that in the chicken α-globin gene domain, the embryonic gene π also does not interact with any of the known distant regulatory elements of the domain[38]. It is thus likely that in the chicken, both α-type and β-type early embryonic globin genes possess an autonomous regulatory system.

In regard to the composition of the β-globin domain ACH in 3-day and 9-day RBCs, we noted several unexpected features. First of all, the HS2 of LCR shows a significant frequency of interaction with the promoter of the βA gene already in 3-day RBCs where the level of this gene transcription is very low. Still it is detectable. Furthermore, it was reported that the βA gene is transcribed to some extend even in 2-day RBCs[39]. Hence, the interaction of the βA gene with the LCR in 3-day RBCs may be functionally relevant. On the other hand, we observed a strong interaction between the promoters of the ρ and βA genes, which was especially prominent in 3-day RBCs in which the level of the βA gene transcription is very low. The functional significance of this interaction is not clear at the moment. The fact that it was detected also in non-erythroid cells suggests that it may reflect some inherent feature of the β-globin gene domain spatial organization that is not directly related to the assembly of ACH. It is of interest that despite the above-mentioned interaction between the ρ and βA promoters, the levels of interaction of these promoters with the LCR-HS2 differ significantly in 3-day RBCs (where the promoter of gene ρ interacts preferentially with the LCR-HS2) and 9-day RBCs (where the promoter of gene βA interacts preferentially with the LCR-HS2). Thus, it is likely that the fine structure of the β-globin ACH is quite different in 3-day and 9-day RBCs. We suggest that in 3-day RBCs, the promoter of gene ρ directly interacts with the LCR while the promoter of gene βA is passively pulled to the complex due to the interaction of βA and ρ. In 9-day RBCs the situation is reversed. Thus, in spite of interaction with each other, the genes ρ and βA demonstrate developmental stage-specific association with the LCR which correlates with the expression status of these genes. Even more evident developmental stage-specific recruitment to the ACH shows the promoter of the minor adult β-type globin gene βH. The promoter of this gene demonstrates either no interaction or a strong interaction with the LCR at days 3 and 9 of the embryo development, respectively. Finally, it should be mentioned that in both 3-day and 9-day RBCs, a significant level of interaction between the βA and βH promoters was observed. Taking into consideration the fact that the anchors placed on the LCR-HS2 and ρ gene promoter demonstrated no interaction with βH and a strong interaction with βA in 3-day RBCs, one must assume that at this developmental stage, the promoters of the two adult β-type globin genes can either interact with each other outside the ACH or the βA promoter can be recruited to the ACH via interaction with the ρ gene promoter. The two alternative configurations of the β-globin domain ACH in 3-day RBCs are shown in Figure4.

Figure 4
figure 4

β-globin active chromatin hubs in the nuclei of 3-day and 9-day red blood cells of chicken embryos. The filled disks designate ACH, and the blue ellipses show protein complexes assembled on the boundary insulators.

In 9-day RBCs, the level of mutual interaction of most of the regulatory elements of the β-globin gene domain increases compared to the levels in 3-day RBCs. These observations may reflect more compact packaging of the domain as a whole. Thus, it should be noted that 9-day RBCs differ drastically from the 3-day RBCs in the total amount of RNA, including rRNA, present in the cell (see the first section of the results). Evidently, the 9-day RBCs are less metabolically active compared to 3-day RBCs. In other words, they possess certain features of mature chicken erythrocytes, which have non-active nuclei. Nevertheless, the globin genes are still transcribed in the majority of 9-day RBCs[40]. Thus, the concept that the apparently more compact spatial organization of the β-globin gene domain in 9-day RBCs may simply reflect the presence of a substantial portion of mature erythrocytes in the analyzed cell population would not be correct. Rather, the β-globin gene domain suffers a certain amount of pressure because it is a part of the chromosomal territory that gradually acquires a more compact configuration. Regardless, it remains active and possesses certain features of an active chromatin domain, including a preferential sensitivity to DNase I, even in the compact nuclei of mature chicken erythrocytes[41]. The spatial isolation of this domain from the surrounding chromatin by the formation of a distinct chromatin loop via interaction of the flanking insulators may be of major importance.

Conclusions

In this work, we have analyzed the spatial configuration of the chicken β-globin gene domain in red blood cells of embryonic and adult lineages. The obtained results strongly support the active chromatin hub model, as it has been found that distant regulatory elements located upstream of the β-globin gene cluster form a complex with the β/Ε enhancer located in the downstream portion of the cluster. Promoters of all β-globin genes, with the exception of embryonic gene Ε, are recruited to this complex, demonstrating a specificity of interaction with the LCR in respect to the developmental stage. The embryonic β-globin gene Ε does not interact with any of the distal regulatory elements; thus, it is likely to possess an autonomous regulatory system. In both 3-day and 9-day RBCs the upstream and the downstream insulators of the chicken β-gene domain interact with each other, thus placing the domain in a separate chromatin loop. This may explain the fact that in chickens, the borders of erythroid cell-specific domain with preferential sensitivity to DNase I perfectly co-localize with the borders of β-globin gene domain marked by insulators.

Methods

Real-time reverse transcriptase polymerase chain reaction (RT-PCR)

Chicken red blood cells (RBCs) were collected from 3- and 9-day embryos and directly used for analysis. Chicken embryonic fibroblasts (CEFs) were isolated from 9-day embryos according to the standard protocol[42] and grown in DMEM supplemented with 8% FBS and 2% chicken serum. Total RNA was extracted using the TRIzol™ Reagent (Invitrogen, Life Technologies, Carlsbad, CA, USA) and was treated with DNase I (Fermentas, Vilnius, Lithuania) to remove residual DNA. RNA (0.5 μg) was reverse transcribed in a total volume of 20 μl for 1 h at 42°C using 0.2 μg of random hexamers and 200 U of reverse transcriptase M-MuLV (Fermentas, Vilnius, Lithuania) in the presence of 20 U of ribonuclease inhibitor (Fermentas, Vilnius, Lithuania). The resulting cDNA was analyzed by TaqMan™ real-time PCR. For β-globin genes, the amplicons were designed in such a way that they spanned intron-exon junctions. For β-actin gene the amplicon was within an exon. The PCR primers and TaqMan™ probes are presented in Table1.

Table 1 Primers and TaqMan™ probes used for analysis of β-globin primary transcripts and β-actin mRNA levels

3C analysis

The 3C analysis was performed as described[43] with minor modifications. Briefly, 7.5×107 cells were fixed in DMEM/serum with 2% formaldehyde for 10 min at room temperature. The reaction was stopped by adding glycine to a concentration of 0.125 M and placing the samples on ice. Cross-linked cells were washed with cold PBS and lysed by incubation for 10 min in ice-cold lysis buffer (10 mM Tris (pH 8.0), 10 mM NaCl, 0.2% NP-40, protease inhibitor cocktail R1321 (Fermentas, Vilnius, Lithuania)). The nuclei were harvested, and aliquots containing 3.75×107 nuclei were frozen in liquid nitrogen and stored at ‐70°C.

Approximately 1.25×107 nuclei were suspended in 0.25 ml of an appropriate 1.2x restriction buffer, SDS was added to a final concentration of 0.3%, and the solution was incubated with shaking at 37°C for 1 h. The solution was then diluted 2-fold with the same restriction buffer, and Triton X-100 was added to a final concentration of 1.8% (to sequester the SDS) followed by incubation at the same conditions for another hour. After incubation, 400 U of MboI (NEB, Hitchin, UK) or NlaIII (Fermentas, Vilnius, Lithuania) was added, and the digestion was performed for 20 h at 37°C with continuous shaking. The reaction was terminated by the addition of SDS to a final concentration of 1.3% and incubation at 65°C for 20 min. The solution was diluted by adding 7 ml of 1x ligation buffer (50 mM Tris (pH 7.5), 10 mM MgCl2, 10 mM DTT, 1 mM ATP). Triton X-100 was added to a final concentration of 1%, and the solution was incubated at 37°C for 1 h while shaking. T4 DNA Ligase (100 U, Fermentas, Vilnius, Lithuania) was then added, and DNA was ligated at 14-16°C for 5 h followed by further incubation at room temperature for 30 min with slow agitation. Cross-links were reversed by incubation at 65°C for 16 h in the presence of 500 μg of Proteinase K (Sigma). After cross-link reversion, 300 μg of RNase A was added, and the RNA was digested at 37°C for 45 min. DNA was purified by extraction with phenol (pH 8.0), followed by precipitation with ethanol in the presence of 120 μg of glycogen (Fermentas, Vilnius, Lithuania).

The efficacy of restriction enzyme digestion in 3C reactions was checked using PCR stop analysis and was found to be more than 80% for different sites throughout the locus under study (data not shown).

To generate a random ligation control, 8 μg of a bacterial artificial chromosome containing the chicken β-globin gene domain along with the flanking areas [CHORI BACPAC Resources Center:TAM32-28 J3] was digested with Sau3A or NlaIII and then ligated at a high DNA concentration.

The TaqMan™ real-time PCR technology was employed to analyze the ligation products. Primers and TaqMan™ probes for PCR were designed using the DNA sequence [GenBank: NW_001471556] and Primer Premier 5 computer software (PREMIER Biosoft International, Palo Alto, CA, USA). The sequences of the primers and TaqMan™ probes are present in Table2. A real-time PCR reaction in a total volume of 20 μl included 20‐100 ng of a 3C DNA template or the same amount of digested and re-ligated chicken genomic DNA along with 3, 0.3, 0,003 or 0,0003 ng of the control BAC template, 1x PCR buffer (50 mM Tris (pH 8.6), 50 mM KCl, 1.5 mM MgCl2, 0.1% Tween 20), 10 pmol of each primer, 5 pmol of a TaqMan™ probe (5′-FAM dye, inside BHQ-1quencher), 4 nmol of each dNTP and 0.75 U of Hot start Taq DNA polymerase (Sibenzyme, Novosibirsk, Russia). The reaction was performed as follows: initial denaturation for 5 min at 94°C, 50 cycles of 15 s at 94°C and 60 s at 60°C followed by reading the plate. All real-time PCR component reactions were aliquoted and stored at ‐70°C until required for PCR with each new primer pair.

Table 2 Primers and TaqMan™ probes used for 3C analysis

To account for the qualitative and quantitative differences of independently prepared 3C samples, the qPCR data for each sample were normalized to the results obtained by using primers for the housekeeping gene, ERCC3 (encodes the ATP-dependent DNA-helicase, a component of the complex TFIIH), located in an area that presumably possesses the same structure of chromatin in all types of cells. The research was approved by the institute ethics governing board.

Abbreviations

ACH:

Active chromatin hub

bp:

base pair

CEF:

Chicken embryonic fibroblast

CTCF:

CCCTC-binding factor

HS:

Site of hypersensitivity to DNase I

LCR:

Locus control region

RBC:

Red blood cell

RT-qPCR:

real-time reverse transcriptase polymerase chain reaction

3C-qPCR:

quantitative chromosome conformation capture assay.

References

  1. Reitman M, Felsenfeld G: Developmental regulation of topoisomerase II sites and DNase I-hypersensitive sites in the chicken beta-globin locus. Mol Cell Biol. 1990, 10: 2774-2786.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. Guerrero G, Delgado-Olguin P, Escamilla-Del-Arenal M, Furlan-Magaril M, Rebollar E, De La Rosa-Velazquez IA, Soto-Reyes E, Rincon-Arano H, Valdes-Quezada C, Valadez-Graham V, Recillas-Targa F: Globin genes transcriptional switching, chromatin structure and linked lessons to epigenetics in cancer: a comparative overview. Comp Biochem Physiol A Mol Integr Physiol. 2007, 147: 750-760. 10.1016/j.cbpa.2006.10.037.

    Article  PubMed  Google Scholar 

  3. Reitman M, Lee E, Westphal H: Function of the upstream hypersensitive sites of the chicken beta-globin gene cluster in mice. Nucleic Acids Res. 1995, 23: 1790-1794. 10.1093/nar/23.10.1790.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Mason MM, Lee E, Westphal H, Reitman M: Expression of the chicken b-globin cluster in mice: correct developmental expression and distributed control. Mol Cell Biol. 1995, 15: 407-414.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Chung JH, Whiteley M, Felsenfeld G: A 5′ element of the chicken beta-globin domain serves as an insulator in human erythroid cells and protects against position effect in Drosophila. Cell. 1993, 74: 505-514. 10.1016/0092-8674(93)80052-G.

    Article  CAS  PubMed  Google Scholar 

  6. Pikaart MJ, Recillas-Targa F, Felsenfeld G: Loss of transcriptional activity of a transgene is accompanied by DNA methylation and histone deacetylation and is prevented by insulators. Gen Dev. 1998, 12: 2852-2862. 10.1101/gad.12.18.2852.

    Article  CAS  Google Scholar 

  7. Recillas-Targa F, Pikaart MJ, Burgess-Beusse B, Bell AC, Litt MD, West AG, Gaszner M, Felsenfeld G: Position-effect protection and enhancer blocking by the chicken beta-globin insulator are separable activities. Proc Natl Acad Sci USA. 2002, 99: 6883-6888. 10.1073/pnas.102179399.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Recillas-Targa F, Bell AC, Felsenfeld G: Positional enhancer-blocking activity of the chicken beta-globin insulator in transiently transfected cells. Proc Natl Acad Sci USA. 1999, 96: 14354-14359. 10.1073/pnas.96.25.14354.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Bell AC, West AG, Felsenfeld G: The protein CTCF is required for the enhancer-blocking activity of vertebrate insulators. Cell. 1999, 98: 387-396. 10.1016/S0092-8674(00)81967-4.

    Article  CAS  PubMed  Google Scholar 

  10. Saitoh N, Bell AC, Recillas-Targa F, West AG, Simpson M, Pikaart M, Felsenfeld G: Structural and functional conservation at the boundaries of the chicken beta-globin domain. EMBO J. 2000, 19: 2315-2322. 10.1093/emboj/19.10.2315.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Burgess-Beusse B, Farrell C, Gaszner M, Litt M, Mutskov V, Recillas-Targa F, Simpson M, West A, Felsenfeld G: The insulation of genes from external enhancers and silencing chromatin. Proc Natl Acad Sci USA. 2002, Suppl 4: 16433-16437.

    Article  Google Scholar 

  12. Choi O, Engel JD: A 3′ enhancer is required for temporal and tissue-specific transcriptional activation of the chicken adult β-globin gene. Nature. 1986, 323: 731-734. 10.1038/323731a0.

    Article  CAS  PubMed  Google Scholar 

  13. Emerson BM, Nickol JM, Jackson PD, Felsenfeld G: Analysis of the tissue-specific enhancer at the 3′ end of the chicken adult beta-globin gene. Proc Natl Acad Sci USA. 1987, 84: 4786-4790. 10.1073/pnas.84.14.4786.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Reitman M, Lee E, Westphal H, Felsenfeld G: Site-independent expression of the chicken beta A-globin gene in transgenic mice. Nature. 1990, 348: 749-752. 10.1038/348749a0.

    Article  CAS  PubMed  Google Scholar 

  15. Recillas-Targa F, Razin SV: Chromatin domains and regulation of gene expression: familiar and enigmatic clusters of chicken globin genes. Crit Rev Eukaryot Gene Expr. 2001, 11: 227-242.

    CAS  PubMed  Google Scholar 

  16. Dillon N, Sabbatini P: Functional gene expression domains: defining the functional units of eukaryotic gene regulation. Bio Essays. 2000, 22: 657-665.

    CAS  Google Scholar 

  17. Razin SV, Farrell CM, Recillas-Targa F: Genomic domains and regulatory elements operating at the domain level. Int Rev Cytol. 2003, 226: 63-125.

    Article  CAS  PubMed  Google Scholar 

  18. Litt MD, Simpson M, Recillas-Targa F, Prioleau MN, Felsenfeld G: Transitions in histone acetylation reveal boundaries of three separately regulated neighboring loci. EMBO J. 2001, 20: 2224-2235. 10.1093/emboj/20.9.2224.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Litt MD, Simpson M, Gaszner M, Allis CD, Felsenfeld G: Correlation between histone lysine methylation and developmental changes at the chicken beta-globin locus. Science. 2001, 293: 2453-2455. 10.1126/science.1064413.

    Article  CAS  PubMed  Google Scholar 

  20. Groudine M, Holtzer H, Scherrer K, Therwath A: Lineage-dependent transcription of globin genes. Cell. 1974, 3: 243-247. 10.1016/0092-8674(74)90138-X.

    Article  CAS  PubMed  Google Scholar 

  21. Bruns GA, Ingram VM: Erythropoiesis in the developing chick embryo. Dev Biol. 1973, 30: 455-459. 10.1016/0012-1606(73)90102-4.

    Article  CAS  PubMed  Google Scholar 

  22. Tolhuis B, Palstra RJ, Splinter E, Grosveld F, de Laat W: Looping and interaction between hypersensitive sites in the active beta-globin locus. Mol Cell. 2002, 10: 1453-1465. 10.1016/S1097-2765(02)00781-5.

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  24. Palstra RJ, Tolhuis B, Splinter E, Nijmeijer R, Grosveld F, de Laat W: The beta-globin nuclear compartment in development and erythroid differentiation. Nat Genet. 2003, 35: 190-194. 10.1038/ng1244.

    Article  CAS  PubMed  Google Scholar 

  25. Hagege H, Klous P, Braem C, Splinter E, Dekker J, Cathala G, de Laat W, Forne T: Quantitative analysis of chromosome conformation capture assays (3C-qPCR). Nat Protoc. 2007, 2: 1722-1733. 10.1038/nprot.2007.243.

    Article  CAS  PubMed  Google Scholar 

  26. Gavrilov A, Eivazova E, Priozhkova I, Lipinski M, Razin S, Vassetzky Y: Chromosome conformation capture (from 3C to 5C) and its ChIP-based modification. Methods Mol Biol. 2009, 567: 171-188. 10.1007/978-1-60327-414-2_12.

    Article  CAS  PubMed  Google Scholar 

  27. Dekker J: The three ’ s of chromosome conformation capture: controls, controls, controls. Nat Methods. 2006, 3: 17-21. 10.1038/nmeth823.

    Article  CAS  PubMed  Google Scholar 

  28. Wijgerde M, Grosveld F, Fraser P: Transcription complex stability and chromatin dynamics in vivo. Nature. 1995, 377: 209-213. 10.1038/377209a0.

    Article  CAS  PubMed  Google Scholar 

  29. West AG, Fraser P: Remote control of gene transcription. Hum Mol Genet. 2005, 14 (1): 101-111. 10.1093/hmg/ddi104.

    Article  Google Scholar 

  30. Zhu X, Ling J, Zhang L, Pi W, Wu M, Tuan D: A facilitated tracking and transcription mechanism of long-range enhancer function. Nucleic Acids Res. 2007, 35: 5532-5544. 10.1093/nar/gkm595.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Kukushkin AN, Svetlikova SB, Pospelov VA: A structure of potentially active and inactive genes of chicken erythrocyte chromatin upon decondensation. Nucleic Acids Res. 1988, 16: 8555-8569. 10.1093/nar/16.17.8555.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Bulger M, Schubeler D, Bender MA, Hamilton J, Farrell CM, Hardison RC, Groudine M: A complex chromatin landscape revealed by patterns of nuclease sensitivity and histone modification within the mouse beta-globin locus. Mol Cell Biol. 2003, 23: 5234-5244. 10.1128/MCB.23.15.5234-5244.2003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Schubeler D, Francastel C, Cimbora DM, Reik A, Martin DIK, Groudine M: Nuclear localization and histone acetilation: a pathway for chromatin opening and transcription activation of the human β-globin locus. Gen. Dev. 2000, 14: 940-950.

    CAS  Google Scholar 

  34. Splinter E, Heath H, Kooren J, Palstra RJ, Klous P, Grosveld F, Galjart N, de Laat W: CTCF mediates long-range chromatin looping and local histone modification in the beta-globin locus. Genes Dev. 2006, 20: 2349-2354. 10.1101/gad.399506.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Kurukuti S, Tiwari VK, Tavoosidana G, Pugacheva E, Murrell A, Zhao Z, Lobanenkov V, Reik W, Ohlsson R: CTCF binding at the H19 imprinting control region mediates maternally inherited higher-order chromatin conformation to restrict enhancer access to Igf2. Proc Natl Acad Sci USA. 2006, 103: 10684-10689. 10.1073/pnas.0600326103.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Yusufzai TM, Felsenfeld G: The 5′-HS4 chicken beta-globin insulator is a CTCF-dependent nuclear matrix-associated element. Proc Natl Acad Sci USA. 2004, 101: 8620-8624. 10.1073/pnas.0402938101.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Wood WI, Felsenfeld G: Chromatin structure of the chicken beta-globin gene region. Sensitivity to DNase I, micrococcal nuclease, and DNase II. J Biol Chem. 1982, 257: 7730-7736.

    CAS  PubMed  Google Scholar 

  38. Ioudinkova ES, Ulianov SV, Bunina D, Iarovaia OV, Gavrilov AA, Razin SV: The inactivation of the pi gene in chicken erythroblasts of adult lineage is not mediated by packaging of the embryonic part of the alpha-globin gene domain into a repressive heterochromatin-like structure. Epigenetics. 2011, 6: 1481-1488. 10.4161/epi.6.12.18215.

    Article  CAS  PubMed  Google Scholar 

  39. Alev C, McIntyre BA, Nagai H, Shin M, Shinmyozu K, Jakt LM, Sheng G: BetaA, the major beta globin in definitive red blood cells, is present from the onset of primitive erythropoiesis in chicken. Dev Dyn. 2008, 237: 1193-1197. 10.1002/dvdy.21510.

    Article  CAS  PubMed  Google Scholar 

  40. Iarovaia OV, Borounova VV, Philonenko ES, Kantidze OL, Vassetzky YS, Razin SV: In embryonic chicken erythrocytes actively transcribed alpha globin genes are not associated with the nuclear matrix. J Cell Biochem. 2009, 106: 170-178. 10.1002/jcb.21987.

    Article  CAS  PubMed  Google Scholar 

  41. Goldsmith ME: Release of a globin gene enriched chromatin fraction from chicken erythrocyte nuclei following DNase II digestion. Nucleic Acids Res. 1981, 9: 6471-6485. 10.1093/nar/9.23.6471.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Silim A, El Azhary MA, Roy RS: A simple technique for preparation of chicken-embryo-skin cell cultures. Avian Dis. 1982, 26: 182-185. 10.2307/1590039.

    Article  CAS  PubMed  Google Scholar 

  43. Gavrilov AA, Razin SV: Spatial configuration of the chicken {alpha}-globin gene domain: immature and active chromatin hubs. Nucleic Acids Res. 2008, 36: 4629-4640. 10.1093/nar/gkn429.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by the Ministry of Science and Education of the Russian Federation (contracts 16.740.11.0353, 14.740.12.1344 and 16.740.11.0483), by the Russian Foundation for Support of Basic Research (grants 11-04-00361-а, 11-04-91334-NNIO_а and 12-04-00036-а), by the Presidium of the Russian Academy of Sciences (grants from the Program on Molecular and Cellular Biology), by a grant from the President of the Russian Federation for young scientists (MK-3813.2012.4) and by Dmitri Zimin’s foundation ”Dynasty”. Experiments were performed using the equipment of IGB RAS facilities supported by the Ministry of Science and Education of the Russian Federation (contract 16.552.11.7067). The work of Sergey Ulianov was also supported in part by the Russian Federation (Order 220) research grant.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Sergey V Razin.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

SVU, AAG, and SVR designed the studies; SVU performed the studies; SVU, AAG, and SVR analyzed the data; and SVU, AAG, and SVR prepared the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Ulianov, S.V., Gavrilov, A.A. & Razin, S.V. Spatial organization of the chicken beta-globin gene domain in erythroid cells of embryonic and adult lineages. Epigenetics & Chromatin 5, 16 (2012). https://doi.org/10.1186/1756-8935-5-16

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1756-8935-5-16

Keywords