Functional and binding studies of HS3.2 of the beta-globin locus control region

Patterns of regulatory activity across diverse human cell types predict tissue identity, transcription factor binding, and long-range interactions

Regulatory elements recruit transcription factors that modulate gene expression distinctly across cell types, but the relationships among these remains elusive. To address this, we analyzed matched DNase-seq and gene expression data for 112 human samples representing 72 cell types. We first defined more than 1800 clusters of DNase I hypersensitive sites (DHSs) with similar tissue specificity of DNase-seq signal patterns. We then used these to uncover distinct associations between DHSs and promoters, CpG islands, conserved elements, and transcription factor motif enrichment. Motif analysis within clusters identified known and novel motifs in cell-type-specific and ubiquitous regulatory elements and supports a role for AP-1 regulating open chromatin. We developed a classifier that accurately predicts cell-type lineage based on only 43 DHSs and evaluated the tissue of origin for cancer cell types. A similar classifier identified three sex-specific loci on the X chromosome, including the XIST lincRNA locus. By correlating DNase I signal and gene expression, we predicted regulated genes for more than 500K DHSs. Finally, we introduce a web resource to enable researchers to use these results to explore these regulatory patterns and better understand how expression is modulated within and across human cell types.

Sequence composition similarities with the 7SL RNA are highly predictive of functional genomic features

Transposable elements derived from the 7SL RNA gene, such as Alu elements in primates, have had remarkable success in several mammalian lineages. The results presented here show a broad spectrum of functions for genomic segments that display sequence composition similarities with the 7SL RNA gene. Using thoroughly documented loci, we report that DNaseI-hypersensitive sites can be singled out in large genomic sequences by an assessment of sequence composition similarities with the 7SL RNA gene. We apply a root word frequency approach to illustrate a distinctive relationship between the sequence of the 7SL RNA gene and several classes of functional genomic features that are not presumed to be of transposable origin. Transposable elements that show noticeable similarities with the 7SL sequence include Alu sequences, as expected, but also long terminal repeats and the 5′-untranslated regions of long interspersed repetitive elements. In sequences masked for repeated elements, we find, when using the 7SL RNA gene as query sequence, distinctive similarities with promoters, exons and distal gene regulatory regions. The latter being the most notoriously difficult to detect, this approach may be useful for finding genomic segments that have regulatory functions and that may have escaped detection by existing methods.

Primary sequence and epigenetic determinants of in vivo occupancy of genomic DNA by GATA1

DNA sequence motifs and epigenetic modifications contribute to specific binding by a transcription factor, but the extent to which each feature determines occupancy in vivo is poorly understood. We addressed this question in erythroid cells by identifying DNA segments occupied by GATA1 and measuring the level of trimethylation of histone H3 lysine 27 (H3K27me3) and monomethylation of H3 lysine 4 (H3K4me1) along a 66 Mb region of mouse chromosome 7. While 91% of the GATA1-occupied segments contain the consensus binding-site motif WGATAR, only ∼0.7% of DNA segments with such a motif are occupied. Using a discriminative motif enumeration method, we identified additional motifs predictive of occupancy given the presence of WGATAR. The specific motif variant AGATAA and occurrence of multiple WGATAR motifs are both strong discriminators. Combining motifs to pair a WGATAR motif with a binding site motif for GATA1, EKLF or SP1 improves discriminative power. Epigenetic modifications are also strong determinants, with the factor-bound segments highly enriched for H3K4me1 and depleted of H3K27me3. Combining primary sequence and epigenetic determinants captures 52% of the GATA1-occupied DNA segments and substantially increases the specificity, to one out of seven segments with the required motif combination and epigenetic signals being bound.

Evaluation of regulatory potential and conservation scores for detecting cis-regulatory modules in aligned mammalian genome sequences

Techniques of comparative genomics are being used to identify candidate functional DNA sequences, and objective evaluations are needed to assess their effectiveness. Different analytical methods score distinctive features of whole-genome alignments among human, mouse, and rat to predict functional regions. We evaluated three of these methods for their ability to identify the positions of known regulatory regions in the well-studied HBB gene complex. Two methods, multispecies conserved sequences and phastCons, quantify levels of conservation to estimate a likelihood that aligned DNA sequences are under purifying selection. A third function, regulatory potential (RP), measures the similarity of patterns in the alignments to those in known regulatory regions. The methods can correctly identify 50%-60% of noncoding positions in the HBB gene complex as regulatory or nonregulatory, with RP performing better than do other methods. When evaluated by the ability to discriminate genomic intervals, RP reaches a sensitivity of 0.78 and a true discovery rate of approximately 0.6. The performance is better on other reference sets; both phastCons and RP scores can capture almost all regulatory elements in those sets along with approximately 7% of the human genome.

Factors Binding a Non-classical Cis-element Prevent Heterochromatin Effects on Locus Control Region Activity

A locus control region (LCR) is a cis-acting gene-regulatory element capable of transferring the expression characteristics of its gene locus of origin to a linked transgene. Furthermore, it can do this independently of the site of integration in the genome of transgenic mice. Although most LCRs contain subelements with classical transcriptional enhancer function, key aspects of LCR activity are supported by cis-acting sequences devoid of the ability to act as direct transcriptional enhancers. Very few of these "non-enhancer" LCR components have been characterized. Consequently, the sequence requirements and molecular bases for their functions, as well as their roles in LCR activity, are poorly understood. We have investigated these questions using the LCR from the mouse T cell receptor (TCR) alpha/Dad1 gene locus. Here we focus on DNase hypersensitive site (HS) 6 of the TCRalpha LCR. HS6 does not support classical enhancer activity, yet has gene regulatory activity in an in vivo chromatin context. We have identified three in vivo occupied factor-binding sites within HS6, two of which interact with Runx1 and Elf-1 factors. Deletion of these sites from the LCR impairs its activity in vivo. This mutation renders the transgene locus abnormally inaccessible in chromatin, preventing the normal function of other LCR subelements and reducing transgene mRNA levels. These data show these factor-binding sites are required for preventing heterochromatin formation and indicate that they function to maintain an active TCRalpha LCR assembly in vivo.

Effect of fetal hemoglobin-stimulating medicines on the interaction of DNA and protein of important erythroid regulatory elements

β-Thalassemia is the most common single gene disorder in the world, which is caused by the imbalance between α-globin chain and β-globin chain synthesis. Several medicines, such as 5-azacytidine, hydroxyurea, cytarabine, vinblatine, butyrate, and myleran, have been shown to be able to reactivate γ-globin chain synthesis during the adult stage, and some of them (5-azacytidine, hydroxyurea, myleran, and butyrate) have been used clinically to treat thalas semia and sickle cell disease. Much research efforts are focusing on the determination of the underlying mechanisms of medicine action. In this experiment, as an effort to probe the underlying mechanism of medicine action, we used ligation-mediated polymerase chain reaction and in vivo footprinting methods to study the DNA-protein interaction at critical erythroid regulatory elements after hydroxyurea or myleran administration to mice. Our results showed that the patterns of in vivo footprints at both the hypersensitive site 2 of the locus control region and the β-globin gene promoter were changed after medicine treatment. We proposed based on these results that the medicines' administration might result in a change in the interaction between trans-acting factors and cis-acting elements at these regions. These changes might influence the assembly of the transcription complex and, lastly, influence the expression of the β-globin gene.Key words: hydroxyurea, in vivo footprinting, ligation-mediated PCR, LCR, β-globin.

Hemoglobin switching: new insights

During the past year, many interesting advances have been made regarding molecular mechanisms controlling beta-like globin gene switching. Throughout the beta locus, -acting elements exist that are dynamically bound by trans-acting proteins, including transcription factors, coactivators, repressors, and chromatin modifiers. Characterization of transcription factors, their interaction with one another, and an ever-increasing role for chromatin structure in gene expression have enhanced understanding of the mechanism of globin gene switching. The studies reviewed here contribute new insights on the interplay between -acting elements, transcription factors, and chromatin modifiers that underlie globin gene switching during development.

PipTools: a computational toolkit to annotate and analyze pairwise comparisons of genomic sequences

Sequence conservation between species is useful both for locating coding regions of genes and for identifying functional noncoding segments. Hence interspecies alignment of genomic sequences is an important computational technique. However, its utility is limited without extensive annotation. We describe a suite of software tools, PipTools, and related programs that facilitate the annotation of genes and putative regulatory elements in pairwise alignments. The alignment server PipMaker uses the output of these tools to display detailed information needed to interpret alignments. These programs are provided in a portable format for use on common desktop computers and both the toolkit and the PipMaker server can be found at our Web site (http://bio.cse.psu.edu/). We illustrate the utility of the toolkit using annotation of a pairwise comparison of the mouse MHC class II and class III regions with orthologous human sequences and subsequently identify conserved, noncoding sequences that are DNase I hypersensitive sites in chromatin of mouse cells.