Defining Candidate Imprinted loci in Bos taurus.

Using a whole-genome assembly of Bos taurus, I applied my bioinformatics strategy to locate candidate imprinting control regions (ICRs) genome-wide. In mammals, genomic imprinting plays essential roles in embryogenesis. In my strategy, peaks in plots mark the locations of known, inferred, and candidate ICRs. Genes in the vicinity of candidate ICRs correspond to potential imprinted genes. By displaying my datasets on the UCSC genome browser, one could view peak positions with respect to genomic landmarks. I give two examples of candidate ICRs in loci that influence spermatogenesis in bulls: CNNM1 and CNR1. I also give examples of candidate ICRs in loci that influence muscle development: SIX1 and BCL6. By examining the ENCODE data reported for mice, I deduced regulatory clues about cattle. I focused on DNase I hypersensitive sites (DHSs). Such sites reveal accessibility of chromatin to regulators of gene expression. For inspection, I chose DHSs in chromatin from mouse embryonic stem cells (ESCs) ES-E14, mesoderm, brain, heart, and skeletal muscle. The ENCODE data revealed that the SIX1 promoter was accessible to the transcription initiation apparatus in mouse ESCs, mesoderm, and skeletal muscles. The data also revealed accessibility of BCL6 locus to regulatory proteins in mouse ESCs and examined tissues.

Along the Bos taurus genome, uncover candidate imprinting control regions.

BACKGROUND: In mammals, Imprinting Control Regions (ICRs) regulate a subset of genes in a parent-of-origin-specific manner. In both human and mouse, previous studies identified a set of CpG-rich motifs occurring as clusters in ICRs and germline Differentially Methylated Regions (gDMRs). These motifs consist of the ZFP57 binding site (ZFBS) overlapping a subset of MLL binding units known as MLL morphemes. MLL or MLL1 (Mixed Lineage Leukemia 1) is a relatively large multidomain protein that plays a central role in the regulation of transcription. The structures of both MLL1 and MLL2 include a domain (MT) that binds CpG-rich DNA and a conserved domain (SET) that methylates lysine 4 in histone H3 producing H3K4me3 marks in chromatin. RESULTS: Since genomic imprinting impacts many developmental and key physiological processes, we followed a previous bioinformatics strategy to pinpoint ICR positions in the Bos taurus genome. Initial genome-wide analyses involved finding the positions of ZFP57 binding sites, and the CpG-rich motifs (ZFBS-morph overlaps) along cattle chromosomal DNA. By creating plots displaying the density of ZFBS-morph overlaps, we removed background noise and thus improved signal detection. With the density-plots, we could view the positions of peaks locating known and candidate ICRs in cattle DNA. Our evaluations revealed the correspondence of peaks in plots to reported known and inferred ICRs/DMRs in cattle. Beside peaks pinpointing such ICRs, the density-plots also revealed additional peaks. Since evaluations validated the robustness of our approach, we inferred that the additional peaks may correspond to candidate ICRs for imprinted gene expression. CONCLUSION: Our bioinformatics strategy offers the first genome-wide approach for systematically localizing candidate ICRs. Furthermore, we have tailored our datasets for upload onto the UCSC genome browser so that researchers could find known and candidate ICRs with respect to a wide variety of annotations at all scales: from the positions of Single Nucleotide Polymorphisms (SNPs), to positions of genes, transcripts, and repeated DNA elements. Furthermore, the UCSC genome browser offers tools to produce enlarged views: to uncover the genes in the vicinity of candidate ICRs and thus discover potential imprinted genes for experimental validations.

Discovering candidate imprinted genes and imprinting control regions in the human genome.

BACKGROUND: Genomic imprinting is a process thereby a subset of genes is expressed in a parent-of-origin specific manner. This evolutionary novelty is restricted to mammals and controlled by genomic DNA segments known as Imprinting Control Regions (ICRs) and germline Differentially Methylated Regions (gDMRs). Previously, I showed that in the mouse genome, the fully characterized ICRs/gDMRs often includes clusters of 2 or more of a set of composite-DNA-elements known as ZFBS-morph overlaps. RESULTS: Because of the importance of the ICRs to regulating parent-of-origin specific gene expression, I developed a genome-wide strategy for predicting their positions in the human genome. My strategy consists of creating plots to display the density of ZFBS-morph overlaps along the entire chromosomal DNA sequences. In initial evaluations, I found that peaks in these plots pinpointed several of the known ICRs/gDMRs along the DNA in chromosomal bands. I deduced that in density-plots, robust peaks corresponded to actual or candidate ICRs in the DNA. By locating the genes in the vicinity of candidate ICRs, I could discover potential imprinting genes. Additionally, my assessments revealed a connection between several of the potential imprinted genes and human developmental anomalies. Examples include Leber congenital amaurosis 11, Coffin-Siris syndrome, progressive myoclonic epilepsy-10, microcephalic osteodysplastic primordial dwarfism type II, and microphthalmia, cleft lip and palate, and agenesis of the corpus callosum. CONCLUSION: With plots displaying the density of ZFBS-morph overlaps, researchers could locate candidate ICRs and imprinted genes. Since the datafiles are available for download and display at the UCSC genome browser, it is possible to examine the plots in the context of Single nucleotide polymorphisms (SNPs) to design experiments to discover novel ICRs and imprinted genes in the human genome.

Datasets on the genomic positions of the MLL1 morphemes, the ZFP57 binding site, and ZFBS-Morph overlaps in the build mm9 of the mouse genome.

While MLL1 activates gene expression in most tissues, ZFP57 represses transcription. MLL1 selectively interacts with a group of nonmethylated DNA sequences known as the MLL1 morphemes. ZFP57 associates with a methylated hexamer (ZFBS), dispersed in the genomic DNA segments known as Imprinted Control Regions (ICRs) and germline Differentially Methylated Regions (gDMRs), to maintain allele-specific gene repression. We have identified a set of composite DNA elements (ZFBS-Morph overlaps) that provides the sequence context of ZFBS in the canonical ICRs/gDMRs. This report provides tables listing the nucleotide sequences of the MLL1 morphemes and ZFBS-Morph overlaps. The report also offers links to the data repository at Purdue University, for downloading the positions of the MLL1 morphemes, the ZFP57 binding site, and the ZFBS-Morph overlaps in the mouse genome.

Imprinted control regions include composite DNA elements consisting of the ZFP57 binding site overlapping MLL1 morphemes.

Mammalian genomes include DNA segments that are imprinted (CpG-methylated) only on one of the two parental chromosomes, leading to parent-of-origin-specific gene expression. The process is regulated by Imprinting Control Regions (ICRs) and germline Differentially Methylated Regions (gDMRs). Previously, ZFP57 was shown to recognize a methylated hexanucleotide in ICRs to maintain allele-specific gene repression. In Bioinformatics analyses, I found that the hexamer occurred frequently in mouse chromosomal DNA, suggesting that beside the ZFP57 binding site (ZFBS), ICRs contained sequence features with unknown characteristics. To identify such features, I examined chromosomal abundance of motifs in which the length of the hexamer was extended by one or several nucleotides. Results led to the discovery of a group of functionally significant composite DNA elements (ZFBS-Morph overlaps) that may play dual roles in the regulation of allele-specific gene expression. Importantly, results of genome-wide evaluations revealed that nearly 90% of the gDMRs included closely-spaced ZFBS-Morph overlaps.

Impact of the MLL1 morphemes on codon utilization and preservation in CpG islands.

Previous studies have shown that Mixed Lineage Leukemia 1 (MLL1 or MLL) binds a group of CpG-rich motifs known as morphemes. To examine whether occurrences of MLL1 morphemes in genomic DNA may influence codon utilization, we analyzed the frequency of various 9-mers in human cDNAs and in total human genomic DNA. We uncovered preferential utilization of GGC for Gly, GCG for Ala, CCG for Pro, and TCG for Ser, in coding sequences (CDSs) that included MLL1 morphemes. We also examined weighted occurrences of CDS 9-mers in a 30-base window that moved along each human chromosome. In plots, we observed peaks with fluctuating intensities. High intensity peaks appeared within promoter and exons localized in CpG islands, encompassing sequences that included MLL1 morphemes. High intensity peaks included CCG/GGC repeats, whose expansion may cause neurological disorders and congenital malformations. Such repeats are generated from overlap of a morpheme (CGCCG/CGGCG), which depending on reading frame and orientation would produce runs of Ala, Gly, or Pro in proteins. Overall, our results point to a role for morpheme occurrences on synonymous codon utilization in human genomic DNA and indicate that regulatory instructions are dispersed not only in promoters but also in exons of human genes.

Discovery of MLL1 binding units, their localization to CpG Islands, and their potential function in mitotic chromatin.

BACKGROUND: Mixed Lineage Leukemia 1 (MLL1) is a mammalian ortholog of the Drosophila Trithorax. In Drosophila, Trithorax complexes transmit the memory of active genes to daughter cells through interactions with Trithorax Response Elements (TREs). However, despite their functional importance, nothing is known about sequence features that may act as TREs in mammalian genomic DNA. RESULTS: By analyzing results of reported DNA binding assays, we identified several CpG rich motifs as potential MLL1 binding units (defined as morphemes). We find that these morphemes are dispersed within a relatively large collection of human promoter sequences and appear densely packed near transcription start sites of protein-coding genes. Genome wide analyses localized frequent morpheme occurrences to CpG islands. In the human HOX loci, the morphemes are spread across CpG islands and in some cases tail into the surrounding shores and shelves of the islands. By analyzing results of chromatin immunoprecipitation assays, we found a connection between morpheme occurrences, CpG islands, and chromatin segments reported to be associated with MLL1. Furthermore, we found a correspondence of reported MLL1-driven "bookmarked" regions in chromatin to frequent occurrences of MLL1 morphemes in CpG islands. CONCLUSION: Our results implicate the MLL1 morphemes in sequence-features that define the mammalian TREs and provide a novel function for CpG islands. Apparently, our findings offer the first evidence for existence of potential TREs in mammalian genomic DNA and the first evidence for a connection between CpG islands and gene-bookmarking by MLL1 to transmit the memory of highly active genes during mitosis. Our results further suggest a role for overlapping morphemes in producing closely packed and multiple MLL1 binding events in genomic DNA so that MLL1 molecules could interact and reside simultaneously on extended potential transcriptional maintenance elements in human chromosomes to transmit the memory of highly active genes during mitosis.

Gene regulation.

This review concisely highlights the complexity of regulatory events. It provides examples of how interconnectivity of regulatory hubs could maintain transcriptional synergy and orchestrate the proper spatial and temporal patterns of gene expression.

Discovering sequences with potential regulatory characteristics.

We developed a computational model to explore the hypothesis that regulatory instructions are context dependent and conveyed through specific 'codes' in human genomic DNA. We provide examples of correlation of computational predictions to reported mapped DNase I hypersensitive segments in the HOXA locus in human chromosome 7. The examples show that statistically significant 9-mers from promoter regions may occur in sequences near and upstream of transcription initiation sites, in intronic regions, and within intergenic regions. Additionally, a subset of 9-mers from coding sequences appears frequently, as clusters, in regulatory regions dispersed in noncoding regions in genomic DNA. The results suggest that the computational model has the potential of decoding regulatory instructions to discover candidate transcription factor binding sites and to discover candidate epigenetic signals that appear in both coding and regulatory regions of genes.

The genome browser at UCSC for locating genes, and much more!

For beginners in the field, this review highlights the key features of the genome browser at UCSC for data display, and provides nearly step-by-step procedures for creating publication quality maps. The browser offers an engine (Blat) for searching a known genomic DNA for correspondence with protein and DNA sequences specified by the user. The results provide links to graphical displays, known as maps. Users can create "designer maps" by adding Tracks to view various types of data and specific landmarks. The browser offers an extensive list of options. They include the position of annotated genes, the position of reference cDNA sequences (RefSeq from GenBank), the position of alternatively spliced mRNA species, and predictions derived from computational models to identify potential transcription start sites and potential protein binding elements in genomic DNA. Several tracks can be tailored for comparative genomics. The browser also offers tracks for displaying large-scale experimental data including gene expression profiles, exon chips, and single-nucleotide-polymorphisms.

Use of genome browsers to locate your favorite genes.

The completion of whole-genome sequencing projects offers the opportunity of creating high-resolution maps of specific segments in a known genomic DNA sequence. For this purpose, several genome browsers have been created. They include the map-view (http://www.ncbi.nlm.nih.gov/mapview/), the Ensembl genome browser (http://www.ensembl.org/), and the genome browser at UCSC (http://genome.ucsc.edu/). For the beginners in the field, through a specific example, this chapter provides a step-by-step procedure for creating a map using the genome browser at UCSC. The example describes mapping, in the human genome, the promoter region of the NF-IL6 gene. The procedure is applicable to creating maps of the desired regions in genomes of other species available at the genome browser at UCSC.

Identification and mapping of paralogous genes on a known genomic DNA sequence.

The completion of whole genome sequencing projects offers the opportunity to examine the organization of genes and the discovery of evolutionarily related genes in a given species. For the beginners in the field, through a specific example, this chapter provides a step-by-step procedure for identifying paralogous genes, using the genome browser at UCSC (http://genome.ucsc.edu/). The example describes identification and mapping in the human genome, the paralogs of TCF12/HTF4. The example identifies TCF3 and TCF4 as paralogs of the TCF12/HTF4 gene. The example also identifies a related sequence, corresponding to a pseudogene, in one of the introns of the JAK2 gene. The procedure described should be applicable to the discovery and creation of maps of paralogous genes in the genomic DNA sequences that are available at the genome browser at UCSC.

A database of 9-mers from promoter regions of human protein-coding genes.

Discovery of lexical characteristics of specific sequence motifs in human genomic DNA can help with predicting and classifying regulatory cis elements according to the genes they control. In lexical models, some "words" may serve as downstream targets of signaling systems, whereas other "words" may specify sequences that selectively control the expression of a subset of genes to produce the various cell types and tissues. To discover lexical features of potential regulatory "words," we have created a database of 9-mers derived from the promoter regions of a subset of human protein-coding genes. This report describes the procedure for extracting information from that database through the web.

A program toolkit for the analysis of regulatory regions of genes.

A major challenge in systems biology is to discover and reconstruct the cis-regulatory networks through which the expression of genes is controlled. Even though a variety of sequences have been shown to interact with the transcription factors that bind DNA, extensive work is needed to discover and classify regulatory "codes" and to elucidate the role played by the sequence context of genomic DNA in the regulation of genes. Databases of sequence elements extracted from regulatory regions may facilitate this process. This report provides a Toolkit and instructions for creating a database for collecting and analyzing 9-base elements (9-mers) from a large collection of DNA sequences. A reference set consisting of all possible 9-mers is included for extracting potential control elements, irrespective of their orientation and order in DNA.

Exploring transcription factor binding properties of several non-coding DNA sequence elements in the human NF-IL6 gene.

We examined several DNA segments upstream of the transcription start site of the human NF-IL6 gene to evaluate the predictions of two computational models developed to identify potential regulatory elements in the non-coding regions of genes. One model, comparative genomics, is based on the hypothesis that functional regulatory sequences can be localized in alignments of genomic DNA from several species. The other model is based on the hypothesis that protein-binding sites in genomic DNA may include sequence elements that occur frequently in proximal promoters of genes. The segments selected for DNA binding and functional evaluations included: (1) two conserved regions identified in multi-species sequence alignments; (2) a region containing several localized hits with 9-mers that ranked highly in studies of proximal promoters of human genes; and (3) two regions that were either GC-rich and/or contained tracts of G. The assays were done under nearly identical experimental conditions, using a cell line (U937) representing human monocytes/macrophages. The experiments also aimed at evaluating what effect, if any, cellular stimulation could have on the interactions of nuclear proteins with naturally occurring GC-rich elements in a human genomic DNA. In DNA binding assays, several complexes were formed with the conserved regions identified in multi-species sequence alignment. Furthermore, these regions were active in functional assays. The region containing several matches with 9-mers derived from proximal promoters of human genes was not conserved but formed several complexes with nuclear proteins including Sp1, Egr-1, and an unidentified protein. In addition, this region was active in functional assays and responded to cellular stimulations. Overall, the results of the assays suggest an important role for the sequence context of genomic DNA in protein binding and selection.

Exploring the characteristics of sequence elements in proximal promoters of human genes.

Central to reconstruction of cis-regulatory networks is identification and classification of naturally occurring transcription factor-binding sites according to the genes that they control. We have examined salient characteristics of 9-mers that occur in various orders and combinations in the proximal promoters of human genes. In evaluations of a dataset derived with respect to experimentally defined transcription initiation sites, in some cases we observed a clear correspondence of highly ranked 9-mers with protein-binding sites in genomic DNA. Evaluations of the larger dataset, derived with respect to the 5' end of human ESTs, revealed that a subset of the highly ranked 9-mers corresponded to sites for several known transcription factor families (including CREB, ETS, EGR-1, SP1, KLF, MAZ, HIF-1, and STATs) that play important roles in the regulation of vertebrate genes. We identified several highly ranked CpG-containing 9-mers, defining sites for interactions with the CREB and ETS families of proteins, and identified potential target genes for these proteins. The results of the studies imply that the CpG-containing transcription factor-binding sites regulate the expression of genes with important roles in pathways leading to cell-type-specific gene expression and pathways controlled by the complex networks of signaling systems.

Regulation of HIV-1 transcription by NF-IL6 in activated Jurkat T cells.

To examine the mechanism of HIV-1 regulation by NF-IL6 in activated human cells, we selected a Jurkat cell line that did not contain endogenous NF-IL6. In this cellular environment, we evaluated the effect of exogenous NF-IL6 on transcription mediated by native and deleted LTR sequences. In Jurkat cells stimulated with LPS and PMA, LTR-mediated transcription was enhanced by NF-IL6. The results of deletion studies revealed a central role for the basal LTR region and the TATA element in the LTR, in upregulation of reporter gene expression by NF-IL6 in activated cells. In the selected cellular environment, regulation of transcription by NF-IL6 was not evident in studies of promoter regions of other genes. The results implied that the basal region of HIV-1 LTR includes molecular properties that support activation of HIV-1 by NF-IL6 in stimulated cells.

Regulation of HIV-1 transcription in activated monocyte macrophages.

DNA-binding and functional assays examined the role played by NF-IL6 in regulation of HIV-1 transcription in human monocyte/macrophages (U937 cells), stimulated with LPS+PMA. When incubated with nuclear extracts from stimulated cells, a region (-189/-147), containing the major NF-IL6-binding sequence and the USF site, interacted selectively with USF1 and USF2. Anti-C/EBPbeta reacted poorly with the complexes produced with the wild-type probe. In contrast, complex formation with NF-IL6 was clearly evident in experiments analyzing a probe containing an insertion in the USF site. In functional assays, increasing concentrations of a decoy against NF-IL6 reduced gene expression from the LTR of the wild-type HIV-1 variant, supporting a critical role for NF-IL6 in regulation of HIV-1 transcription in stimulated monocyte/macrophages. The decoy also reduced gene expression from a deletion construct lacking NF-IL6-binding sequences. The results implied that in LPS+PMA-stimulated monocyte/macrophages, the endogenous NF-IL6 could act via a site-independent pathway in upregulation of HIV-1 transcription. Analysis of a short DNA segment, containing the -189/-147 region, suggested functional interactions of NF-IL6 and USF. In activated cells exogenous NF-IL6 enhanced dramatically gene expression through a short DNA segment containing the NF-kappaB sites, supporting functional interactions of NF-IL6 and NF-kappaB.

Organization of the promoter region of the human NF-IL6 gene.

In monocyte/macrophages, the human NF-IL6 gene was activated by LPS or PMA. However, a robust response required stimulation of cells with both LPS and PMA. To examine the molecular basis of this response, we isolated human genomic DNA and determined the nucleotide sequence of a segment (6.4 kb) that included the transcription initiation site of the gene. The unique sequences in the 6.4-kb DNA include several potential transcription factor-binding elements that may explain the molecular basis of the activation of the human NF-IL6 gene by signaling molecules that control the immune and inflammatory responses. Deletion analysis localized an LPS+PMA responsive region downstream position -287, with respect to the transcription initiation site of the NF-IL6 gene. The responsive region includes a potential site for interactions with CREB and a region (-287 to -247) that interacts with SP1 and SP3. In functional assays, the potential CREB site responded to cellular stimulation. The region that interacted with SP1 and SP3 augmented the overall level of activity produced in response to LPS+PMA.