Correlation Between DNase I Hypersensitive Site Distribution and Gene Expression in HeLa S3 Cells

Mapping DNase I hypersensitive sites (DHSs) within nuclear chromatin is a traditional and powerful method of identifying genetic regulatory elements. DHSs have been mapped by capturing the ends of long DNase I-cut fragments (>100,000 bp), or 100-1200 bp DNase I-double cleavage fragments (also called double-hit fragments). But next generation sequencing requires a DNA library containing DNA fragments of 100-500 bp. Therefore, we used short DNA fragments released by DNase I digestion to generate DNA libraries for next generation sequencing. The short segments are 100-300 bp and can be directly cloned and used for high-throughput sequencing. We identified 83,897 DHSs in 2,343,479 tags across the human genome. Our results indicate that the DHSs identified by this DHS assay are consistent with those identified by longer fragments in previous studies. We also found: (1) the distribution of DHSs in promoter and other gene regions of similarly expressed genes differs among different chromosomes; (2) silenced genes had a more open chromatin structure than previously thought; (3) DHSs in 3'untranslated regions (3'UTRs) are negatively correlated with level of gene expression.     

Identification of a conserved sequence in the non-coding regions of many human genes.

We have analyzed a sequence of approximately 70 base pairs (bp) that shows a high degree of similarity to sequences present in the non-coding regions of a number of human and other mammalian genes. The sequence was discovered in a fragment of human genomic DNA adjacent to an integrated hepatitis B virus genome in cells derived from human hepatocellular carcinoma tissue. When one of the viral flanking sequences was compared to nucleotide sequences in GenBank, more than thirty human genes were identified that contained a similar sequence in their non-coding regions. The sequence element was usually found once or twice in a gene, either in an intron or in the 5' or 3' flanking regions. It did not share any similarities with known short interspersed nucleotide elements (SINEs) or presently known gene regulatory elements. This element was highly conserved at the same position within the corresponding human and mouse genes for myoglobin and N-myc, indicating evolutionary conservation and possible functional importance. Preliminary DNase I footprinting data suggested that the element or its adjacent sequences may bind nuclear factors to generate specific DNase I hypersensitive sites. The size, structure, and evolutionary conservation of this sequence indicates that it is distinct from other types of short interspersed repetitive elements. It is possible that the element may have a cis-acting functional role in the genome.     

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.     

Chromatin Accessibility Data Sets Show Bias Due to Sequence Specificity of the DNase I Enzyme

Background

DNase I is an enzyme which cuts duplex DNA at a rate that depends strongly upon its chromatin environment. In combination with high-throughput sequencing (HTS) technology, it can be used to infer genome-wide landscapes of open chromatin regions. Using this technology, systematic identification of hundreds of thousands of DNase I hypersensitive sites (DHS) per cell type has been possible, and this in turn has helped to precisely delineate genomic regulatory compartments. However, to date there has been relatively little investigation into possible biases affecting this data.

Results

We report a significant degree of sequence preference spanning sites cut by DNase I in a number of published data sets. The two major protocols in current use each show a different pattern, but for a given protocol the pattern of sequence specificity seems to be quite consistent. The patterns are substantially different from biases seen in other types of HTS data sets, and in some cases the most constrained position lies outside the sequenced fragment, implying that this constraint must relate to the digestion process rather than events occurring during library preparation or sequencing.

Conclusions

DNase I is a sequence-specific enzyme, with a specificity that may depend on experimental conditions. This sequence specificity is not taken into account by existing pipelines for identifying open chromatin regions. Care must be taken when interpreting DNase I results, especially when looking at the precise locations of the reads. Future studies may be able to improve the sensitivity and precision of chromatin state measurement by compensating for sequence bias.     

Principal component analysis- and tensor decomposition-based unsupervised feature extraction to select more suitable differentially methylated cytosines: Optimization of standard deviation versus state-of-the-art methods

In contrast to RNA-seq analysis, which has various standard methods, no standard methods for identifying differentially methylated cytosines (DMCs) exist. To identify DMCs, we tested principal component analysis and tensor decomposition-based unsupervised feature extraction with optimized standard deviation, which has been shown to be effective for differentially expressed gene (DEG) identification. The proposed method outperformed certain conventional methods, including those that assume beta-binomial distribution for methylation as the proposed method does not require this, especially when applied to methylation profiles measured using high throughput sequencing. DMCs identified by the proposed method also significantly overlapped with various functional sites, including known differentially methylated regions, enhancers, and DNase I hypersensitive sites. The proposed method was applied to data sets retrieved from The Cancer Genome Atlas to identify DMCs using American Joint Committee on Cancer staging system edition labels. This suggests that the proposed method is a promising standard method for identifying DMCs.     

Moloney murine leukemia virus-induced tumors: recombinant proviruses in active chromatin regions.

The DNase I sensitivity of chromosomal DNA regions carrying integrated proviral genomes of Moloney (M-MuLV) and AKR Murine Leukemia Virus (AKR-MuLV), and the cellular homologue of the mos-gene (c-mos) of Moloney Sarcoma Virus (MSV) were studied in tumor tissues of leukemic mice. The genetically transmitted sequences of M-MuLV, AKR-MuLV, and the c-mos gene are all in DNase I resistant chromatin conformations in M-MuLV-induced tumors. Each M-MuLV-induced tumor contained at least one somatically acquired integrated recombinant MuLV genome that displayed two main characteristic features of active chromatin: a) a configuration hypersensitive to DNase I, and b) extensive hypomethylation. DNase I hypersensitive sites were mapped at the junction of cellular sequences and the 5'-viral large terminal repeat (LTR). Expression of a recombinant MuLV seems therefore to be a necessary feature to maintain the transformed state.     

DNase I hypersensitive sites flank the mouse class II major histocompatibility complex during B cell development

The mouse class II major histocompatibility complex (MHC) encodes a polymorphic, multigene family important in the immune response, and is expressed mainly on mature B cells, on certain types of dendritic cells and is also inducible by gamma-interferon on antigen presenting cells. To study the regulatory elements which control this expression pattern, we have examined the chromatin structure flanking the class II MHC region, in particular during B cell differentiation. Using a panel of well-characterised mouse cell lines specific for different stages of B cell development (pre-B, B, plasma cell) as well as non-B cell lines, we have mapped the DNase I hypersensitive (DHS) sites adjacent to the mouse MHC class II region. The results presented show, for the first time that there are specific hypersensitive sites flanking the class II MHC locus during pre B cell, B cell and plasma cell stages of B cell differentiation, irrespective of the status of class II MHC expression. These hypersensitive sites are not found in T cell, fibroblast or uninduced myelomonocytic cell lines. This suggests that these DHS sites define a developmentally stable, chromatin structure, which can be used as a marker of B cell lineage commitment and may indicate that a combination of these hypersensitive sites reflect regulatory proteins involved in the immediate expression of a particular class II MHC gene or possibly control of the entire locus.     

Preferential binding of adriamycin and nogalamycin to DNase-I hypersensitive sites of Sarcoma-180 chromatin

Binding of nogalamycin and adriamycin with Sarcoma-180 ascites tumor cell chromatin was studied by a spectrofluorometric method. There was significant reduction in the number of available drug binding sites per nucleotide when the chromatin was digested with DNase I for a period which releases only 7% of the chromosomal DNA. Results indicate preferential binding of these drugs with DNase I hypersensitive sites of chromatin. The DNase-I hypersensitive sites of chromatin were shown to correlate to the sequences required for gene expression. Further digestion with DNase I increases availability of drug binding sites, probably due to relaxation of the compact chromatin.