Periodicity of DNA bend sites in human epsilon-globin gene region. Possibility of sequence-directed nucleosome phasing

Gene therapy of hemoglobinopathies: progress and future challenges

Recently, gene therapy clinical trials have been successfully applied to hemoglobinopathies, such as sickle cell disease (SCD) and β-thalassemia. Among the great discoveries that led to the design of genetic approaches to cure these disorders is the discovery of the β-globin locus control region and several associated transcription factors, which determine hemoglobin switching as well as high-level, erythroid-specific expression of genes at the ß-globin locus. Moreover, increasing evidence shows that lentiviral vectors are efficient tools to insert large DNA elements into nondividing hematopoietic stem cells, showing reassuring safe integration profiles. Alternatively, genome editing could restore expression of fetal hemoglobin or target specific mutations to restore expression of the wild-type β-globin gene. The most recent clinical trials for β-thalassemia and SCD are showing promising outcomes: patients were able to discontinue transfusions or had reduced transfusion requirements. However, toxic myeloablation and the high cost of current ex vivo hematopoietic stem cell gene therapy platforms represent a barrier to a widespread application of these approaches. In this review, we summarize these gene therapy strategies and ongoing clinical trials. Finally, we discuss possible strategies to improve outcomes, reduce myeloablative regimens and future challenges to reduce the cost of gene therapy platform.

A conserved regulatory element in the mammalian β-globin promoters

We provide here evidence for a conserved regulatory element for transcription of the β-family globin genes based on a comparative study of 32 genes from 16 mammals. The element is characterized by the appearance of AA or TT dinucleotides in the A + T-rich region located 200-400 bp upstream of the cap sites. G-tracts 3-5 nucleotides long exist between the A + T-rich region and the conserved transcription factor binding sites (GATA-1 site and the CACCC, CCAAT, and ATA boxes) apparently dividing the regions. The average periodicity of AA or TT dinucleotides in the region from a total of 18 β-family globin genes from four species was approximately 10 bp, suggesting that the DNA in these regions shows right-handed superhelicity. The proposed biological function of this element is to adjust the spatial positions for the first interaction of the transcription factor(s) which can recognize specific DNA sequences in the presence of packed chromatin.

MFTZ-1 reduces constitutive and inducible HIF-1α accumulation and VEGF secretion independent of its topoisomerase II inhibition

The macrolide compound MFTZ-1 has been identified as a novel topoisomerase II (Top2) inhibitor with potent in vitro and in vivo anti-tumour activities. In this study, we further examined the effects of MFTZ-1 on hypoxia-inducible factor-1α (HIF-1α) accumulation, vascular endothelial growth factor (VEGF) secretion and angiogenesis. MFTZ-1 reduced HIF-1α accumulation driven by hypoxia or growth factors in human cancer cells. Mechanistic studies revealed that MFTZ-1 did not affect the degradation of HIF-1α protein or the level of HIF-1α mRNA. By contrast, MFTZ-1 apparently inhibited constitutive and inducible activation of both phosphatidylinositol-3-kinase (PI3K)-Akt and p42/p44 mitogen-activated protein kinase (MAPK) pathways. Further studies revealed that MFTZ-1 abrogated the HIF-1α-driven increase in VEGF mRNA and protein secretion. MFTZ-1 also lowered the basal level of VEGF secretion. The results reveal an important feature that MFTZ-1 can reduce constitutive, HIF-1α-independent VEGF secretion and concurrently antagonize inducible, HIF-1α-dependent VEGF secretion. Moreover, MFTZ-1 disrupted tube formation of human umbilical vein endothelial cells (HUVECs) stimulated by hypoxia with low-concentration serum or by serum at normoxia, and inhibited HUVECs migration at normoxia. MFTZ-1 also prevented microvessel outgrowth from rat aortic ring. These data reflect the potent anti-angiogenesis of MFTZ-1 under different conditions. Furthermore, using specific small interfering RNA targeting Top2α or Top2-defective HL60/MX2 cells, we showed that MFTZ-1 affected HIF-1α accumulation and HUVECs tube formation irrelevant to its Top2 inhibition. Taken together, our data collectively reveal that MFTZ-1 reduces constitutive and inducible HIF-1α accumulation and VEGF secretion possibly via PI3K-Akt and MAPK pathways, eliciting anti-angiogenesis independently of its Top2 inhibition.

Analysis of DNA bending by transient electric birefringence

Transient electric birefringence has been used to analyze DNA bending in six restriction fragments containing 171, 174, 207, 263, 289, and 471 bp in three different low ionic strength buffers. The target fragments contain sequences corresponding to the apparent bend centers in pUC19 and Litmus 28, previously identified by the circular permutation assay (Strutz, K.; Stellwagen, N. C. Electrophoresis 1996, 17, 989-995). The target fragments migrate anomalously slowly in polyacrylamide gels and exhibit birefringence relaxation times that are shorter than those of restriction fragments of the same size, taken from nonbent regions of the same plasmids. Apparent bend angles ranging from 30 degrees to 41 degrees were calculated for the target fragments by tau-ratio method. The bend angles of four of the target fragments were independent of temperature from 4 degrees C to 20 degrees C, but decreased when the temperature was increased to 37 degrees C. The bend angles of the other two target fragments were independent of temperature over the entire range examined, 4 degrees -37 degrees C. Hence, the thermal stability of sequence-dependent bends in random-sequence DNA is variable. The bend angles of five of the six target fragments were independent of the presence or absence of Mg2+ ions in the solution, indicating most of the target fragments were stably bent or curved, rather than anisometrically flexible. Restriction fragments containing 219 and 224 bp, with sequences somewhat offset from the sequence of the 207 bp fragment, were also studied. Comparison of the tau-ratios of these overlapping fragments allowed both the bend angle and bend position to be independently determined. These methods should be useful for analyzing sequence-dependent bending in other random-sequence DNAs.

Dinucleosome DNA of human K562 cells: experimental and computational characterizations

Dinucleosome formation is the first step in the organization of the higher order chromatin structure. With the ultimate aim of elucidating the dinucleosome structure, we constructed a library of human dinucleosome DNA. The library consists of PCR-amplifiable DNA fragments obtained by treatment of nuclei of erythroid K562 cells with micrococcal nuclease followed by extraction of DNA and adaptor ligation to the blunt-ended DNA fragments. The library was then cloned using a plasmid vector and the sequences of the clones were determined. The dominating clones containing the Alu elements were removed. A total of 1002 clones, which comprised a dinucleosome database, contained 84 and 918 clones from the clones before and after removing Alu elements, respectively. Approximately 70% of the clones were between 300 and 400 bp in size and they were distributed to various locations of all chromosomes except the Y chromosome. The clones containing A(2)N(8)A(2)N(8)A(2) or T(2)N(8)T(2)N(8)T(2) sequences were classified into three types, Type I (N shape), Type II (V shape) and Type III (M shape) according to DNA curvature plots. The locations of experimentally determined curved DNA segments matched well with the calculated ones though the clones of Types I and III showed additional curved DNA segments as revealed by the curvature plots. The distributions of complementary dinucleotides in the nucleosome DNA, at the ends of the dinucleosome DNA clones, allowed us to predict the positions of the nucleosome dyad axis, and estimate the size of the nucleosome core DNA, 125nt. The distributions of AA and TT dinucleotides, as well as other RR and YY dinucleotides, showed a periodicity with an average period of 10.4 bases, close to the values observed before. Mapping of nucleosome positions in the dinucleosome database based on the observed periodicity revealed that the nucleosomes were separated by a linker of 7.5+ approximately 10 x n nt. This indicates that the nucleosome-nucleosome orientations are, typically, halfway between parallel and antiparallel. Also an important finding is that the distributions of AA/TT and other RR/YY dinucleotides, apparently, reflect both DNA curvature and DNA bendability, cooperatively contributing to the nucleosome formation.

What positions nucleosomes?--A model

Here we propose a new determinant for localization of nucleosomes along genomic DNA, in addition to sequence-dependent features. The new specific class of chromatin scaling signals involves curved DNA. According to the observed positional distribution of DNA curvature, the new synchronizing signal occurs once per four nucleosomes on average. This new factor in nucleosome positioning should substantially influence the efficiency of biological reactions through regulatory factors microscopically and the entire chromatin structure through the 30 nm fiber structure macroscopically. Allocation of the new type of signals is found to be fixed evolutionarily although they could be shifted in accordance with the hierarchy of functional genomic structures.

Ligand-dependent transcriptional enhancement by DNA curvature between two half motifs of the estrogen response element in the human estrogen receptor alpha gene

We previously reported five DNA bend sites (ERB-4 to -1, and ERB+1) in the promoter region of the human estrogen receptor alpha (ERalpha) gene [FEBS Lett. 444 (1999) 117]. One of these sites, ERB-2, was accompanied by two half motifs of the estrogen response element (ERE) and several short poly(dA)(.)poly(dT) tracts including an A(4) tract located next to a half ERE motif. This A(4) tract and the 20 bp immediate flanking sequence containing a half ERE motif (T3B) exhibited DNA curvature. Transcription assays using luciferase as a reporter gene indicated that T3B sequence conferred positive estrogen responsiveness. Mutations introduced in this sequence indicated that both bendability and estrogen responsiveness were synergistically associated with the A(4) tract located next to the half ERE motif. This motif and a mutant sequence, GGTTA, had affinity for ERalpha protein, which seems to account for ERalpha protein binding to the region without an ERE motif. These findings suggest that some DNA curvature acts as a transcriptional modulator by modifying the state of ligand effects.

Mapping of intrinsic bent DNA sites in the upstream region of DNA puff BhC4-1 amplified gene

We have identified bent DNA sites in the distal and proximal DNA puff BhC4-1 amplified gene promoter region of Bradysia hygida. The 2D modeling of the 3D DNA path and the ENDS ratio values calculated in this promoter region resulted in the identification of ten pronounced bent sites named BhC4B - 9 to + 1. The 1847 bp fragment (- 3697 to - 1850) in relation to the transcription start site shows multiple bending sites, BhC4B - 9 to BhC4B - 4, with periodicity approximately 300 bp. The analysis of the other identified bent region, starting at position - 957, reveals that the BhC4B + 1 bent site colocalizes with the putative BhC4-1 minimal promoter. The sequence analysis of bent site BhC4B - 4 shows a distribution of dA*dT at approximately 10 bp intervals between the middle of each tract, but intervals with more than one turn, approximately 20 bp, two helix turns, were detected in the other bent sites described here. The bent sites BhC4B - 6 and BhC4B - 4, contain two consensus sequences, with 60 bp each. The apparent molecular weight of fragments in the BhC4-1 promoter region were estimated in agarose gels and compared with the data obtained in polyacrylamide gels without and with ethidium bromide. The mobility reduction ratios (R-values) were determined, and a high R-value, 1.80, for a 1215 bp fragment in the distal promoter region and a 1.23 significant R-value for a 662 bp fragment in the proximal segment were found. To further analyze the predicted bent DNA sites in these fragments, the 2D trajectories of the 3D DNA path and other parameters, AT percentage, roll angle, ENDS ratio and DeltaG, were determined. The role of these bent sites in the BhC4-1 transcription regulation is discussed.

Intrinsically bent DNA in the promoter regions of the yeast GAAL1-10 and GAL80 genes

Circular permutation analysis has detected fairly strong sites of intrinsic DNA bending on the promoter regions of the yeast GAL1-10 and GAL80 genes. These bends lie in functionally suggestive locations. On the promoter of the GAL1-10 structural genes, strong bends bracket nucleosome B, which lies between the UAS(G) and the GAL1 TATA. These intrinsic bends could help position nucleosome B. Nucleosome B plus two other promoter nucleosomes protect the TATA and start site elements in the inactive state of expression but are completely disrupted (removed) when GAL1-10 expression is induced. The strongest intrinsic bend ( approximately 70 degrees ) lies at the downstream edge of nucleosome B; this places it approximately 30 base pairs upstream of the GAL1 TATA, a position that could allow it to be involved in GAL1 activation in several ways, including the recruitment of a yeast HMG protein that is required for the normally robust level of GAL1 expression in the induced state (Paull, T., Carey, M., and Johnson, R. (1996) Genes Dev. 10, 2769-2781). On the regulatory gene GAL80, the single bend lies in the non-nucleosomal hypersensitive region, between a GAL80-specific far upstream promoter element and the more gene-proximal promoter elements. GAL80 promoter region nucleosomes contain no intrinsically bent DNA.

Localization of curved DNA and its association with nucleosome phasing in the promoter region of the human estrogen receptor alpha gene

We determined DNA bend sites in the promoter region of the human estrogen receptor (ER) gene by the circular permutation assay. A total of five sites (ERB-4 to -1, and ERB+1) mapped in the 3 kb region showed an average distance of 688 bp. Most of the sites were accompanied by short poly(dA) x poly(dT) tracts including the potential bend core sequence A2N8A2N8A2 (A/A/A). Fine mapping of the ERB-2 site indicated that this A/A/A and the 20 bp immediate flanking sequence containing one half of the estrogen response element were the sites of DNA curvature. All of the experimentally mapped bend sites corresponded to the positions of DNA curvature as well as to nucleosomes predicted by computer analysis. In vitro nucleosome mapping at ERB-2 revealed that the bend center was located 10-30 bp from the experimental and predicted nucleosome dyad axes.

Binding of HMG-I(Y) elicits structural changes in a silencer of the human beta-globin gene

Proteins involved in repression of the human beta-globin gene may be useful in the treatment of sickle cell anemia, in conjunction with therapy to reactivate fetal globin genes. If there is a reciprocal elevation of gamma-globin expression upon repression, this approach could be useful in additional hemoglobinopathies. We previously showed that repression of the beta-globin gene appears to be mediated through two DNA sequences, silencers I and II, and identified a protein termed BP1 which binds to both silencer sequences. In this study, we cloned two cDNAs encoding proteins which bind to an oligonucleotide in silencer I containing a BP1 binding site. These cDNAs correspond to HMG-I and HMG-Y, isoforms regarded as architectural proteins. We demonstrate that binding of HMG-I(Y) to this oligonucleotide causes bending/flexure of the DNA. HMG-I(Y) also binds to a second oligonucleotide containing a BP1 binding site located in a negative control region upstream of the delta-globin gene, suggesting a role for HMG-I(Y) in repression of adult globin genes. Expression studies revealed that HMG-I(Y) is ubiquitously expressed in human tissues that do not express beta-globin, being present in 48 of 50 tissues and six hematopoietic cell lines examined. Furthermore, HMG-I(Y) expression is down-regulated during differentiation of primary erythroid cells. We present a model in which HMG-I(Y) alters DNA conformation to allow binding of repressor proteins, and in which the relative amount of HMG-I(Y) helps to determine the repressive state of the beta-globin gene.

Expression-dependent perturbation of nucleosomal phases at HS2 of the human beta-LCR: possible correlation with periodic bent DNA

DNA bend sites appear periodically at average intervals of 680 bp, corresponding to a length of four nucleosomes, in the human epsilon-, beta- and Ggamma-Agamma-psibeta-globin gene regions. We found that the HS2 region flanked by two DNA bend sites accommodated five nucleosomes and they were regularly phased throughout the region with the exception of that located in the middle, which corresponded to the precise location of HS2 and included the binding site for NF-E2. There appeared to be several phases in this region in the reconstituted chromatin and in erythroid K562 cells where the globin genes are expressed, whereas only one phase was adopted in non-erythroid HeLa cells. Meanwhile, almost unique phases were adopted at the flanking bend sites in vitro as well as in vivo. Sequences of 30 bp containing the bend centers cloned into the vector alone showed identical nucleosomal phases to those observed with the in vitro and in vivo experiments and removing the bend sites caused disruption of the phases at the bend sites as well as those in their direct vicinity. Finally, the nucleosome in this HS2 region had an inhibitory effect on NF-E2 binding, although remodeling occurred with the nuclear extract from K562 cells in the presence of ATP. This suggests that HS2 is placed at a region of weak nucleosome phasing activity along with factor binding sites.

Conservation and continuity of periodic bent DNA in genomic rearrangements between the c-myc and immunoglobulin heavy chain mu loci

Periodic bent DNA was mapped in the human c- myc and immunoglobulin heavy chain mu (Ig mu) loci. A total of 12 DNA bend sites in the c- myc gene and 11 sites in the Ig mu locus were aligned at average intervals of 694.2 +/- 281.4 and 654.5 +/- 222.7 bp respectively. Although some of the bend sites retained the distance of 700 bp, their periodicity was disturbed at several locations, including the exons of the c- myc gene and the enhancer element present in the Ig mu locus. Analysis of rearrangements that resulted in tumorigenesis of lymphocytes showed that the continuity of DNA bend sites was conserved in three lymphoma cell lines, Manca, BL22 and Ramos, suggesting that the genomic rearrangements gain stability by retaining their periodicity. This adds further evidence that the periodic bent DNA plays a crucial role in genomic structure.

Curvature and sequence analysis of eukaryotic promoters

We have calculated the curvature of 504 eukaryotic promoters predicted by the bent A-tract model of Bolshoy et al. (Proc. Natl. Acad. Sci. USA, 88(6), pp. 2312-16) and the bent non-A-tract models of Calladine et al. (J. Mol. Biol., 201, pp. 127-37) and Satchwell et al. (J. Mol. Biol., 191, pp. 659-75) and found in each case a correlation between TBP binding sites and DNA curvature. Characterizing the TBP binding sites revealed that in addition to the classical TATA box (TATAAA) five more elements occur significantly often in the promoters, nearly all of them being one point mutations of the classical TATA box element. Separate curvature calculations for promoters with canonical and non-canonical TATA boxes have shown that in both cases the strong curvature of the helix axes in the domain of the binding sites is maintained (classical TBP binding sites: + 64-135%, non-classical TBP binding sites: + 27-49%). These results support the proposition that beside DNA flexibility and DNA-protein interactions intrinsic curvature of DNA is one further important criterion for the recognition of different DNA elements by TBP.

Segmented structure of separate and transposable DNA and RNA elements as suggested by their size distributions

A collection of about 1000 different eukaryotic and prokaryotic DNA mobile and separate elements is compiled from literature-transposons, plasmids, extrachromosomal circular DNA, insertion sequences, as well as viral genomes and separate genome segments. Only small elements are collected, upto 2000 base pairs. Analysis of the sequence length distributions of the elements reveals that certain sizes are clearly preferred, namely those which correspond to multiples of about 345 bp in eukaryotes and multiples of about 210 bp in prokaryotes. This provides additional evidence in support of the theory (1) that segmented structure is characteristic of not only protein-coding sequences (2) but rather of genomes in general. In particular, it confirms the prediction (1) that mobile and separate elements would also be segmented.

Common DNA structural features exhibited by eukaryotic ribosomal gene promoters

Nucleotide sequences of DNA regions containing eukaryotic ribosomal promoters were analysed using strategies designed to reveal sequence-directed structural features. DNA curvature, duplex stability and pattern of twist angle variation were studied by computer modelling. Although ribosomal promoters are known to lack sequence homology (unless very closely related species are considered), investigation of these structural characteristics uncovered striking homologies in all the taxonomic groups examined so far. This wide conservation of DNA structures, while DNA sequence is not conserved, suggests that the determined structures are fundamental for ribosomal promoter function. Moreover, this result agrees well with the recent observations showing that RNA polymerase I transcription factors have not evolved as intensively as previously suspected.

Conservation and periodicity of DNA bend sites in the human beta-globin gene locus

A total of seven DNA bend sites were mapped in the 4.4-kilobase human beta-globin gene region by the circular permutation assay. The periodicity of these sites (except one) was about every 700 (average 685.5 +/- 267.7) base pairs. All of the sites contained the sequence feature of short poly(dA) tracts, which are typical of DNA bending. The relative positions of the sites to the cap site were identical to those in the epsilon-globin gene region, suggesting that the bend sites were conserved during molecular evolution of the two globin genes. To explain this periodicity and conservation of the sites within the evolutionary unstable noncoding regions, we focused upon the appearance of a potential bend core sequence, A2N8A2N8A2 (A/A/A), and its complement, T2N8T2N8T2 (T/T/T). These sequences appeared in or very close to most of the bend sites of the globin gene regions, whereas other A+T-rich sequences or candidates for DNA bending did not. The distances between any two of the core sequences in the entire beta-globin locus showed a strong bias to a length of about 700 base pairs and its multiples, suggesting that the periodicity exists throughout the locus. The data presented here strengthen the idea of sequence-directed nucleosome phasing.