Functional analysis of DNA sequences located within a cluster of DNase U hypersensitive sites colocalizing with a MAR element at the upstream border of the chicken α‐globin gene domain

Matrix association region/scaffold attachment region: the crucial player in defining the positions of chromosome breaks mediated by bile acid-induced apoptosis in nasopharyngeal epithelial cells

Background

It has been found that chronic rhinosinusitis (CRS) increases the risk of developing nasopharyngeal carcinoma (NPC). CRS can be caused by gastro-oesophageal reflux (GOR) that may reach nasopharynx. The major component of refluxate, bile acid (BA) has been found to be carcinogenic and genotoxic. BA-induced apoptosis has been associated with various cancers. We have previously demonstrated that BA induced apoptosis and gene cleavages in nasopharyngeal epithelial cells. Chromosomal cleavage occurs at the early stage of both apoptosis and chromosome rearrangement. It was suggested that chromosome breaks tend to cluster in the region containing matrix association region/scaffold attachment region (MAR/SAR). This study hypothesised that BA may cause chromosome breaks at MAR/SAR leading to chromosome aberrations in NPC. This study targeted the AF9 gene located at 9p22 because 9p22 is a deletion hotspot in NPC.

Methods

Potential MAR/SAR sites were predicted in the AF9 gene by using MAR/SAR prediction tools. Normal nasopharyngeal epithelial cells (NP69) and NPC cells (TWO4) were treated with BA at neutral and acidic pH. Inverse-PCR (IPCR) was used to identify chromosome breaks in SAR region (contains MAR/SAR) and non-SAR region (does not contain MAR/SAR). To map the chromosomal breakpoints within the AF9 SAR and non-SAR regions, DNA sequencing was performed.

Results

In the AF9 SAR region, the gene cleavage frequencies of BA-treated NP69 and TWO4 cells were significantly higher than those of untreated control. As for the AF9 non-SAR region, no significant difference in cleavage frequency was detected between untreated and BA-treated cells. A few breakpoints detected in the SAR region were mapped within the AF9 region that was previously reported to translocate with the mixed lineage leukaemia (MLL) gene in an acute lymphoblastic leukaemia (ALL) patient.

Conclusions

Our findings suggest that MAR/SAR may be involved in defining the positions of chromosomal breakages induced by BA. Our report here, for the first time, unravelled the relation of these BA-induced chromosomal breakages to the AF9 chromatin structure.

Electronic supplementary material

The online version of this article (10.1186/s12920-018-0465-4) contains supplementary material, which is available to authorized users.

Conservation of the PRM1→PRM2→TNP2 Domain

In mouse and human, the genes encoding protamines PRM1, PRM2 and transition protein TNP2 are found clustered together on chromosome 16. In addition, these three genes lie in the same orientation to one another and are coordinately expressed in a haploid-specific manner during spermatogenesis. Previously, we have shown that the human PRM1 --> PRM2 --> TNP2 locus exists as a single chromatin domain bounded by two male germ cell-specific MARs, i.e. Matrix Attachment Regions. A third, somatic-specific MAR element lies immediately 3' of the PRM1 --> PRM2 --> TNP2 domain. This MAR maps to a conserved CpG island 5' of the human SOCS-1 gene. Similarly, two candidate MARs flank the mouse Prm1 --> Prm2 --> Tnp2 domain. Comparative analysis of the mouse and human promoter regions identified several conserved regulatory motifs for each of the genes of this cluster. This further establishes the synteny of this region. Global structural similarities and the functional relevance of the associated candidate regulatory elements are discussed.

Early replication timing of the chicken alpha-globin gene domain correlates with its open chromatin state in cells of different lineages

The vertebrate alpha-globin gene domain is an open chromatin domain overlapping a neighboring house-keeping gene. The tissue-specific cluster of alpha-globin genes and the overlapping housekeeping gene share the same replication origin. We have studied the replication timing of chicken alpha-globin genes in cells of different lineages using the FISH-based approach and found that alpha-globin genes replicate early both in erythroid and in non-erythroid cells, i.e. regardless of their transcriptional activity. Early replication timing of chicken alpha-globin genes in cells of different lineages was in good correlation with the open chromatin configuration of the alpha-globin gene domain in both erythroid and non-erythroid cells. We propose that active transcription of the housekeeping gene overlapping the alpha-globin gene domain enables an access of Origin Recognition Complex (ORC) proteins to the replication origin resulting in early replication of alpha-globin genes even in non-erythroid cells.

Spatial configuration of the chicken alpha-globin gene domain: immature and active chromatin hubs

The spatial configuration of the chicken α-globin gene domain in erythroid and lymphoid cells was studied by using the Chromosome Conformation Capture (3C) approach. Real-time PCR with TaqMan probes was employed to estimate the frequencies of cross-linking of different restriction fragments within the domain. In differentiated cultured erythroblasts and in 10-day chick embryo erythrocytes expressing ‘adult’ α^A and α^D globin genes the following elements of the domain were found to form an ‘active’ chromatin hub: upstream Major Regulatory Element (MRE), −9 kb upstream DNase I hypersensitive site (DHS), −4 kb upstream CpG island, α^D gene promoter and the downstream enhancer. The α^A gene promoter was not present in the ‘active’ chromatin hub although the level of α^A gene transcription exceeded that of the α^D gene. Formation of the ‘active’ chromatin hub was preceded by the assembly of multiple incomplete hubs containing MRE in combination with either −9 kb DHS or other regulatory elements of the domain. These incomplete chromatin hubs were present in proliferating cultured erythroblasts which did not express globin genes. In lymphoid cells only the interaction between the α^D promoter and the CpG island was detected.

YY1 and GATA-1 interaction modulate the chicken 3'-side alpha-globin enhancer activity

Studying the chicken alpha-globin domain as a model system of gene regulation, we have previously identified contiguous silencer-enhancer elements located on the 3'-side of the domain. To better characterize the enhancer we performed a systematic functional analysis to define its expression influence range and the ubiquitous and stage-specific transcriptional regulators interacting with this control element. In contrast to previous reports, we found that, in addition to a core element that includes three GATA-1 binding sites, the enhancer incorporates a 120 base-pair DNA fragment where EKLF, NF-E2 and a fourth GATA-1 factor could interact. Functional experiments demonstrate that the enhancer activity over the adult alpha(D) promoter is differentially regulated. We found that the transcriptional factor Ying Yang 1 (YY1) binds to the 120 base-pair DNA fragment and its effect over the enhancer activity is GATA-1-dependent. In addition, we characterize a novel physical interaction between GATA-1 and YY1 that influences the enhancer function. Experiments using a histone deacetylation inhibitor indicate that, in pre-erythroblasts, the enhancer down-regulation could be influenced by a closed chromatin conformation. Our observations show that the originally defined enhancer possesses a more complex composition than previously assumed. We propose that its activity is modulated through differential nuclear factor interactions and chromatin modifications at distinct erythroid stages.

Genomic domains and regulatory elements operating at the domain level

The sequencing of the complete genomes of several organisms, including humans, has so far not contributed much to our understanding of the mechanisms regulating gene expression in the course of realization of developmental programs. In this so-called "postgenomic" era, we still do not understand how (if at all) the long-range organization of the genome is related to its function. The domain hypothesis of the eukaryotic genome organization postulates that the genome is subdivided into a number of semiindependent functional units (domains) that may include one or several functionally related genes, with these domains having well-defined borders, and operate under the control of special (domain-level) regulatory systems. This hypothesis was extensively discussed in the literature over the past 15 years. Yet it is still unclear whether the hypothesis is valid or not. There is evidence both supporting and questioning this hypothesis. The most conclusive data supporting the domain hypothesis come from studies of avian and mammalian beta-globin domains. In this review we will critically discuss the present state of the studies on these and other genomic domains, paying special attention to the domain-level regulatory systems known as locus control regions (LCRs). Based on this discussion, we will try to reevaluate the domain hypothesis of the organization of the eukaryotic genome.

The 33 kb transcript of the chicken ?-globin gene domain is part of the nuclear matrix

Giant nuclear transcripts, and in particular the RNAs of the globin gene domains which are much larger than their canonical pre-mRNAs, have been an enigma for many years. We show here that in avian erythroblastosis virus (AEV)-transformed chicken erythroleukaemic cells, where globin gene expression is abortive, the whole domain of alpha-globin genes is transcribed for about 33 kb in the globin direction and that this RNA is part of the nuclear matrix. Northern blot hybridisation with strand-specific riboprobes, recognising genes and intergenic sequences, and RT-PCR with downstream primers, show that the continuous full domain transcript (FDT) starts in the vicinity of a putative LCR and includes all the genes as well as known regulatory sites, the replication origin, and the DNA loop anchorage region in the upstream area. Absent in chicken fibroblasts, the globin FDT overlaps the major part of the ggPRX housekeeping gene that is transcribed in the opposite direction. RT-PCR and in situ hybridisation with genic and extra-genic globin probes demonstrated that the globin FDT is a component of the nuclear matrix. We suggest that the globin FDTs keep the domain in an active state, and the globin RNAs on the processing pathway are a component of the nuclear matrix. They may take part in the dynamic nuclear architecture when productively processed, or turn over slowly when globins are not synthesised.

Genomic DNA breakpoints in AML1/RUNX1 and ETO cluster with topoisomerase II DNA cleavage and DNase I hypersensitive sites in t(8;21) leukemia

The translocation t(8;21)(q22;q22) is one of the most frequent chromosome translocations in acute myeloid leukemia (AML). AML1/RUNX1 at 21q22 is involved in t(8;21), t(3;21), and t(16;21) in de novo and therapy-related AML and myelodysplastic syndrome as well as in t(12;21) in childhood B cell acute lymphoblastic leukemia. Although DNA breakpoints in AML1 and ETO (at 8q22) cluster in a few introns, the mechanisms of DNA recombination resulting in t(8;21) are unknown. The correlation of specific chromatin structural elements, i.e., topoisomerase II (topo II) DNA cleavage sites, DNase I hypersensitive sites, and scaffold-associated regions, which have been implicated in chromosome recombination with genomic DNA breakpoints in AML1 and ETO in t(8;21) is unknown. The breakpoints in AML1 and ETO were clustered in the Kasumi 1 cell line and in 31 leukemia patients with t(8;21); all except one had de novo AML. Sequencing of the breakpoint junctions revealed no common DNA motif; however, deletions, duplications, microhomologies, and nontemplate DNA were found. Ten in vivo topo II DNA cleavage sites were mapped in AML1, including three in intron 5 and seven in intron 7a, and two were in intron 1b of ETO. All strong topo II sites colocalized with DNase I hypersensitive sites and thus represent open chromatin regions. These sites correlated with genomic DNA breakpoints in both AML1 and ETO, thus implicating them in the de novo 8;21 translocation.

Extensive methylation of a part of the CpG island located 3.0-4.5 kbp upstream to the chicken alpha-globin gene cluster may contribute to silencing the globin genes in non-erythroid cells

Here, we show that in the chicken genome, the domain of alpha-globin genes is preceded by a CpG island of which the downstream part ( approximately 0.65 kbp) is heavily methylated in lymphoid cells; it is either non-methylated or undermethylated in erythroid cells. Recombinant plasmids were constructed with the corresponding DNA fragment (called "uCpG") placed upstream to a reporter CAT gene expressed from the promoter of the alpha(D) chicken globin gene. Selective methylation of CpG dinucleotides within the uCpG fragment suppressed fivefold the expression of the CAT gene, when neither this gene itself nor the alpha(D) promoter were methylated. Methylation of CpG dinucleotides within the alpha(D) gene promoter did not modify the suppression effect exerted by methylated uCpG. We interpret these results within the frame of the hypothesis postulating, that methylation of the upstream CpG island of the chicken alpha-globin gene domain may play an essential role in silencing the alpha-globin genes in non-erythroid cells.