Linking Hematopoietic Differentiation to Co-Expressed Sets of Pluripotency-Associated and Imprinted Genes and to Regulatory microRNA-Transcription Factor Motifs.

Maintenance of cell pluripotency, differentiation, and reprogramming are regulated by complex gene regulatory networks (GRNs) including monoallelically-expressed imprinted genes. Besides transcriptional control, epigenetic modifications and microRNAs contribute to cellular differentiation. As a model system for studying the capacity of cells to preserve their pluripotency state and the onset of differentiation and subsequent specialization, murine hematopoiesis was used and compared to embryonic stem cells (ESCs) as a control. Using published microarray data, the expression profiles of two sets of genes, pluripotent and imprinted, were compared to a third set of known hematopoietic genes. We found that more than half of the pluripotent and imprinted genes are clearly upregulated in ESCs but subsequently repressed during hematopoiesis. The remaining genes were either upregulated in hematopoietic progenitors or in differentiated blood cells. The three gene sets each consist of three similarly behaving gene groups with similar expression profiles in various lineages of the hematopoietic system as well as in ESCs. To explain this co-regulation behavior, we explored the transcriptional and post-transcriptional mechanisms of pluripotent and imprinted genes and their regulator/target miRNAs in six different hematopoietic lineages. Therewith, lineage-specific transcription factor (TF)-miRNA regulatory networks were generated and their topologies and functional impacts during hematopoiesis were analyzed. This led to the identification of TF-miRNA co-regulatory motifs, for which we validated the contribution to the cellular development of the corresponding lineage in terms of statistical significance and relevance to biological evidence. This analysis also identified key miRNAs and TFs/genes that might play important roles in the derived lineage networks. These molecular associations suggest new aspects of the cellular regulation of the onset of cellular differentiation and during hematopoiesis involving, on one hand, pluripotent genes that were previously not discussed in the context of hematopoiesis and, on the other hand, involve genes that are related to genomic imprinting. These are new links between hematopoiesis and cellular differentiation and the important field of epigenetic modifications.

Cellular functions of genetically imprinted genes in human and mouse as annotated in the gene ontology.

By analyzing the cellular functions of genetically imprinted genes as annotated in the Gene Ontology for human and mouse, we found that imprinted genes are often involved in developmental, transport and regulatory processes. In the human, paternally expressed genes are enriched in GO terms related to the development of organs and of anatomical structures. In the mouse, maternally expressed genes regulate cation transport as well as G-protein signaling processes. Furthermore, we investigated if imprinted genes are regulated by common transcription factors. We identified 25 TF families that showed an enrichment of binding sites in the set of imprinted genes in human and 40 TF families in mouse. In general, maternally and paternally expressed genes are not regulated by different transcription factors. The genes Nnat, Klf14, Blcap, Gnas and Ube3a contribute most to the enrichment of TF families. In the mouse, genes that are maternally expressed in placenta are enriched for AP1 binding sites. In the human, we found that these genes possessed binding sites for both, AP1 and SP1.

Computational studies of imprinted genes.

Computational studies on imprinted genes can have very different purposes: one major aim of these studies is the identification of DNA elements that distinguish imprinted genes from biallelically expressed genes. Comparative studies may help to identify imprinting regulatory elements and to understand common mechanisms of imprinted gene regulation in mammalian species. To date, the continuously growing number of genomic and epigenetic data sets makes detailed, genome-wide analyses on imprinted genes feasible. However, imprinted genes are characterized by genomic features that can influence statistics and can make such studies difficult. Hence, comparative computational studies can get very complex and require a tight interaction between bioinformaticians and biologists. Furthermore, analyses of raw data that are generated by micro-array hybridization and high-throughput sequencing technologies require computational approaches that have been designed especially for the epigenetic field. This chapter gives an overview about databases and software that is suitable for analyses of imprinted genes. Furthermore, possible difficulties that are typical for computational and statistical analyses of imprinted genes are described.

Unique patterns of evolutionary conservation of imprinted genes.

During mammalian evolution, complex systems of epigenetic gene regulation have been established: Epigenetic mechanisms control tissue-specific gene expression, X chromosome inactivation in females and genomic imprinting. Studying DNA sequence conservation in imprinted genes, it becomes evident that evolution of gene function and evolution of epigenetic gene regulation are tightly connected. Furthermore, comparative studies allow the identification of DNA sequence features that distinguish imprinted genes from biallelically expressed genes. Among these features are CpG islands, tandem repeats and retrotransposed elements that are known to play major roles in epigenetic gene regulation. Currently, more and more genetic and epigenetic data sets become available. In future, such data sets will provide the basis for more complex investigations on epigenetic variation in human populations. Therein, an exciting topic will be the genetic and epigenetic variability of imprinted genes and its input on human disease.

Imprinted genes show unique patterns of sequence conservation.

BACKGROUND: Genomic imprinting is an evolutionary conserved mechanism of epigenetic gene regulation in placental mammals that results in silencing of one of the parental alleles. In order to decipher interactions between allele-specific DNA methylation of imprinted genes and evolutionary conservation, we performed a genome-wide comparative investigation of genomic sequences and highly conserved elements of imprinted genes in human and mouse. RESULTS: Evolutionarily conserved elements in imprinted regions differ from those associated with autosomal genes in various ways. Whereas for maternally expressed genes strong divergence of protein-encoding sequences is most prominent, paternally expressed genes exhibit substantial conservation of coding and noncoding sequences. Conserved elements in imprinted regions are marked by enrichment of CpG dinucleotides and low (TpG+CpA)/(2·CpG) ratios indicate reduced CpG deamination. Interestingly, paternally and maternally expressed genes can be distinguished by differences in G+C and CpG contents that might be associated with unusual epigenetic features. Especially noncoding conserved elements of paternally expressed genes are exceptionally G+C and CpG rich. In addition, we confirmed a frequent occurrence of intronic CpG islands and observed a decelerated degeneration of ancient LINE-1 repeats. We also found a moderate enrichment of YY1 and CTCF binding sites in imprinted regions and identified several short sequence motifs in highly conserved elements that might act as additional regulatory elements. CONCLUSIONS: We discovered several novel conserved DNA features that might be related to allele-specific DNA methylation. Our results hint at reduced CpG deamination rates in imprinted regions, which affects mostly noncoding conserved elements of paternally expressed genes. Pronounced differences between maternally and paternally expressed genes imply specific modes of evolution as a result of differences in epigenetic features and a special response to selective pressure. In addition, our data support the potential role of intronic CpG islands as epigenetic key regulatory elements and suggest that evolutionary conserved LINE-1 elements fulfill regulatory functions in imprinted regions.

Divergence of imprinted genes during mammalian evolution.

BACKGROUND: In contrast to the majority of mammalian genes, imprinted genes are monoallelically expressed with the choice of the active allele depending on its parental origin. Due to their special inheritance patterns, maternally and paternally expressed genes might be under different evolutionary pressure. Here, we aimed at assessing the evolutionary history of imprinted genes. RESULTS: In this study, we investigated the conservation of imprinted genes in vertebrate genomes and their exposition to natural selection. In a genome-wide comparison, orthologs of imprinted genes show a stronger divergence on cDNA and protein level in mammals. This pattern is most pronounced for maternally expressed genes in rodents in comparison to their non-rodent orthologs. The divergence is not attributable to increased mutation of CpG positions. It is contrasted by strong conservation of paternally expressed genes in mouse and rat. Interestingly, we found that the early divergence of imprinted genes was accompanied by an unusually strict conservation of their paralogs. CONCLUSIONS: The apparent degeneration of maternally expressed genes may reflect a relaxation of selective pressure due to counteracting effects on maternal and embryonic fitness. Functional redundancy provided by the presence of highly conserved (non-imprinted) paralogs may have facilitated the divergence. Moreover, intensification of imprinting in modern rodents seems to have shifted the evolutionary fate of imprinted genes towards strong purifying selection.

The Begain gene marks the centromeric boundary of the imprinted region on mouse chromosome 12.

Although the central portion of the imprinted region on mouse chromosome 12 has been intensively analysed in the past, little is known about the neighbouring centromeric genes. A DNA sequence comparison shows that the region upstream of Dlk1 and Gtl2 is dominated by an expanded cluster of repetitive elements in the mouse. These elements separate the paternally expressed Dlk1 gene from the centromeric Begain gene. Despite the long physical distance to the IG-DMR imprinting centre, Begain is subjected to genomic imprinting. Similar to the ovine Begain gene, the homologous mouse gene encodes two different transcript variants, one of which shows a strong bias toward paternal transcription. Nevertheless, imprinting effects do not spread further centromeric to the Wdr25 gene, which is biallelically expressed as the previously studied neighbouring Wars and Yy1 genes.

Expression profile and transcription factor binding site exploration of imprinted genes in human and mouse.

BACKGROUND: In mammals, imprinted genes are regulated by an epigenetic mechanism that results in parental origin-specific expression. Though allele-specific regulation of imprinted genes has been studied for several individual genes in detail, little is known about their overall tissue-specific expression patterns and interspecies conservation of expression. RESULTS: We performed a computational analysis of microarray expression data of imprinted genes in human and mouse placentae and in a variety of adult tissues. For mouse, early embryonic stages were also included. The analysis reveals that imprinted genes are expressed in a broad spectrum of tissues for both species. Overall, the relative tissue-specific expression levels of orthologous imprinted genes in human and mouse are not highly correlated. However, in both species distinctive expression profiles are found in tissues of the endocrine pathways such as adrenal gland, pituitary, pancreas as well as placenta. In mouse, the placental and embryonic expression patterns of imprinted genes are highly similar. Transcription factor binding site (TFBS) prediction reveals correlation of tissue-specific expression patterns and the presence of distinct TFBS signatures in the upstream region of human imprinted genes. CONCLUSION: Imprinted genes are broadly expressed pre- and postnatally and do not exhibit a distinct overall expression pattern when compared to non-imprinted genes. The relative expression of most orthologous gene pairs varies significantly between human and mouse suggesting rapid species-specific changes in gene regulation. Distinct expression profiles of imprinted genes are confined to certain human and mouse hormone producing tissues, and placentae. In contrast to the overall variability, distinct expression profiles and enriched TFBS signatures are found in human and mouse endocrine tissues and placentae. This points towards an important role played by imprinted gene regulation in these tissues.

Identifying CpG islands by different computational techniques.

CpG islands (CGIs) are generally regarded as important epigenetic regulatory elements due to their association with promoter regions. However, identification of functional CGIs is hampered by repetitive elements and species-specific particularities. Here, we compared the performance of different CGI detection programs on genomic sequences of human and mouse genes. Although mouse CGIs are shorter and G+C poorer than their human counterparts, the different tools tested in our study reliably identify CGIs in promoter regions in both species. Our study confirms that substantially fewer murine than human CGIs coincide with repetitive elements and indicates that such CGIs are subject to accelerated cytosine deamination. In addition, CpG depletion appears to anticorrelate with the epigenetic features of functional regulatory CGIs. Taking into account different deamination rates in unmethylated CGIs versus those in methylated CGIs might support the detection of functional CGIs in other species for which there is little epigenetic information available.

Structural conservation versus functional divergence of maternally expressed microRNAs in the Dlk1/Gtl2 imprinting region.

BACKGROUND: MicroRNAs play an important functional role in post-transcriptional gene regulation. One of the largest known microRNA clusters is located within the imprinted Dlk1/Gtl2 region on human chromosome 14 and mouse chromosome 12. This cluster contains more than 40 microRNA genes that are expressed only from the maternal chromosome in mouse. RESULTS: To shed light on the function of these microRNAs and possible crosstalk between microRNA-based gene regulation and genomic imprinting, we performed extensive in silico analyses of the microRNAs in this imprinted region and their predicted target genes.Bioinformatic analysis reveals that these microRNAs are highly conserved in both human and mouse. Whereas the microRNA precursors at this locus mostly belong to large sequence families, the mature microRNAs sequences are highly divergent. We developed a target gene prediction approach that combines three widely used prediction methods and achieved a sufficiently high prediction accuracy. Target gene sets predicted for individual microRNAs derived from the imprinted region show little overlap and do not differ significantly in their properties from target genes predicted for a group of randomly selected microRNAs. The target genes are enriched with long and GC-rich 3' UTR sequences and are preferentially annotated to development, regulation processes and cell communication. Furthermore, among all analyzed human and mouse genes, the predicted target genes are characterized by consistently higher expression levels in all tissues considered. CONCLUSION: Our results suggest a complex evolutionary history for microRNA genes in this imprinted region, including an amplification of microRNA precursors in a mammalian ancestor, and a rapid subsequent divergence of the mature sequences. This produced a broad spectrum of target genes. Further, our analyses did not uncover a functional relation between imprinted gene regulation of this microRNA-encoding region, expression patterns or functions of predicted target genes. Specifically, our results indicate that these microRNAs do not regulate a particular set of genes. We conclude that these imprinted microRNAs do not regulate a particular set of genes. Rather, they seem to stabilize expression of a variety of genes, thereby being an integral part of the genome-wide microRNA gene regulatory network.

Inter-individual variation of DNA methylation and its implications for large-scale epigenome mapping.

Genomic DNA methylation profiles exhibit substantial variation within the human population, with important functional implications for gene regulation. So far little is known about the characteristics and determinants of DNA methylation variation among healthy individuals. We performed bioinformatic analysis of high-resolution methylation profiles from multiple individuals, uncovering complex patterns of inter-individual variation that are strongly correlated with the local DNA sequence. CpG-rich regions exhibit low and relatively similar levels of DNA methylation in all individuals, but the sequential order of the (few) methylated among the (many) unmethylated CpGs differs randomly across individuals. In contrast, CpG-poor regions exhibit substantially elevated levels of inter-individual variation, but also significant conservation of specific DNA methylation patterns between unrelated individuals. This observation has important implications for experimental analysis of DNA methylation, e.g. in the context of epigenome projects. First, DNA methylation mapping at single-CpG resolution is expected to uncover informative DNA methylation patterns for the CpG-poor bulk of the human genome. Second, for CpG-rich regions it will be sufficient to measure average methylation levels rather than assaying every single CpG. We substantiate these conclusions by an in silico benchmarking study of six widely used methods for DNA methylation mapping. Based on our findings, we propose a cost-optimized two-track strategy for mammalian methylome projects.

CpG island mapping by epigenome prediction.

CpG islands were originally identified by epigenetic and functional properties, namely, absence of DNA methylation and frequent promoter association. However, this concept was quickly replaced by simple DNA sequence criteria, which allowed for genome-wide annotation of CpG islands in the absence of large-scale epigenetic datasets. Although widely used, the current CpG island criteria incur significant disadvantages: (1) reliance on arbitrary threshold parameters that bear little biological justification, (2) failure to account for widespread heterogeneity among CpG islands, and (3) apparent lack of specificity when applied to the human genome. This study is driven by the idea that a quantitative score of "CpG island strength" that incorporates epigenetic and functional aspects can help resolve these issues. We construct an epigenome prediction pipeline that links the DNA sequence of CpG islands to their epigenetic states, including DNA methylation, histone modifications, and chromatin accessibility. By training support vector machines on epigenetic data for CpG islands on human Chromosomes 21 and 22, we identify informative DNA attributes that correlate with open versus compact chromatin structures. These DNA attributes are used to predict the epigenetic states of all CpG islands genome-wide. Combining predictions for multiple epigenetic features, we estimate the inherent CpG island strength for each CpG island in the human genome, i.e., its inherent tendency to exhibit an open and transcriptionally competent chromatin structure. We extensively validate our results on independent datasets, showing that the CpG island strength predictions are applicable and informative across different tissues and cell types, and we derive improved maps of predicted "bona fide" CpG islands. The mapping of CpG islands by epigenome prediction is conceptually superior to identifying CpG islands by widely used sequence criteria since it links CpG island detection to their characteristic epigenetic and functional states. And it is superior to purely experimental epigenome mapping for CpG island detection since it abstracts from specific properties that are limited to a single cell type or tissue. In addition, using computational epigenetics methods we could identify high correlation between the epigenome and characteristics of the DNA sequence, a finding which emphasizes the need for a better understanding of the mechanistic links between genome and epigenome.

Tandem repeats in the CpG islands of imprinted genes.

In contrast to most genes in mammalian genomes, imprinted genes are monoallelically expressed depending on the parental origin of the alleles. Imprinted gene expression is regulated by distinct DNA elements that exhibit allele-specific epigenetic modifications, such as DNA methylation. These so-called differentially methylated regions frequently overlap with CpG islands. Thus, CpG islands of imprinted genes may contain special DNA elements that distinguish them from CpG islands of biallelically expressed genes. Here, we present a detailed study of CpG islands of imprinted genes in mouse and in human. Our study shows that imprinted genes more frequently contain tandem repeat arrays in their CpG islands than randomly selected genes in both species. In addition, mouse imprinted genes more frequently possess intragenic CpG islands that may serve as promoters of allele-specific antisense transcripts. This feature is much less pronounced in human, indicating an interspecies variability in the evolution of imprinting control elements.

CpG island methylation in human lymphocytes is highly correlated with DNA sequence, repeats, and predicted DNA structure.

CpG island methylation plays an important role in epigenetic gene control during mammalian development and is frequently altered in disease situations such as cancer. The majority of CpG islands is normally unmethylated, but a sizeable fraction is prone to become methylated in various cell types and pathological situations. The goal of this study is to show that a computational epigenetics approach can discriminate between CpG islands that are prone to methylation from those that remain unmethylated. We develop a bioinformatics scoring and prediction method on the basis of a set of 1,184 DNA attributes, which refer to sequence, repeats, predicted structure, CpG islands, genes, predicted binding sites, conservation, and single nucleotide polymorphisms. These attributes are scored on 132 CpG islands across the entire human Chromosome 21, whose methylation status was previously established for normal human lymphocytes. Our results show that three groups of DNA attributes, namely certain sequence patterns, specific DNA repeats, and a particular DNA structure, are each highly correlated with CpG island methylation (correlation coefficients of 0.64, 0.66, and 0.49, respectively). We predicted, and subsequently experimentally examined 12 CpG islands from human Chromosome 21 with unknown methylation patterns and found more than 90% of our predictions to be correct. In addition, we applied our prediction method to analyzing Human Epigenome Project methylation data on human Chromosome 6 and again observed high prediction accuracy. In summary, our results suggest that DNA composition of CpG islands (sequence, repeats, and structure) plays a significant role in predisposing CpG islands for DNA methylation. This finding may have a strong impact on our understanding of changes in CpG island methylation in development and disease.

High-resolution map and imprinting analysis of the Gtl2-Dnchc1 domain on mouse chromosome 12.

The imprinted Dlk1-Dio3 region on mouse chromosome 12 contains six imprinted genes and a number of maternally expressed snoRNAs and miRNAs. Here we present a high-resolution sequence analysis of the 1.1-Mb segment telomeric to Gtl2 in mouse and a homology comparison to the human. Ppp2r5c and Dnchc1 at the telomeric end of the analyzed sequence are biallelically expressed, suggesting that the imprinted domain does not extend beyond the paternally expressed Dio3 gene. RT-PCR experiments support the predicted presence of a maternally expressed intergenic transcript(s) encompassing Gtl2, Rian, and Mirg. These maternally expressed genes, and also the intergenic transcript(s), show pronounced expression in the adult mouse brain, whereas the paternally transcribed Dio3 and the nonimprinted Ppp2r5c and Dnchc1 are expressed in different tissues. Hence, tissue-specific coregulation of maternally expressed genes might be an important feature of this domain.

BiQ Analyzer: visualization and quality control for DNA methylation data from bisulfite sequencing.

SUMMARY: Manual processing of DNA methylation data from bisulfite sequencing is a tedious and error-prone task. Here we present an interactive software tool that provides start-to-end support for this process. In an easy-to-use manner, the tool helps the user to import the sequence files from the sequencer, to align them, to exclude or correct critical sequences, to document the experiment, to perform basic statistics and to produce publication-quality diagrams. Emphasis is put on quality control: The program automatically assesses data quality and provides warnings and suggestions for dealing with critical sequences. The BiQ Analyzer program is implemented in the Java programming language and runs on any platform for which a recent Java virtual machine is available. AVAILABILITY: The program is available without charge for non-commercial users and can be downloaded from http://biq-analyzer.bioinf.mpi-inf.mpg.de/

Evolution of the Beckwith-Wiedemann syndrome region in vertebrates.

In the animal kingdom, genomic imprinting appears to be restricted to mammals. It remains an open question how structural features for imprinting evolved in mammalian genomes. The clustering of genes around imprinting control centers (ICs) is regarded as a hallmark for the coordinated imprinted regulation. Hence imprinted clusters might be structurally distinct between mammals and nonimprinted vertebrates. To address this question we compared the organization of the Beckwith Wiedemann syndrome (BWS) gene cluster in mammals, chicken, Fugu (pufferfish), and zebrafish. Our analysis shows that gene synteny is apparently well conserved between mammals and birds, and is detectable but less pronounced in fish. Hence, clustering apparently evolved during vertebrate radiation and involved two major duplication events that took place before the separation of the fish and mammalian lineages. A cross-species analysis of imprinting center regions showed that some structural features can already be recognized in nonimprinted amniotes in one of the imprinting centers (IC2). In contrast, the imprinting center IC1 is absent in chicken. This suggests a progressive and stepwise evolution of imprinting control elements. In line with that, imprinting centers in mammals apparently exhibit a high degree of structural and sequence variation despite conserved epigenetic marking.

Asymmetric regulation of imprinting on the maternal and paternal chromosomes at the Dlk1-Gtl2 imprinted cluster on mouse chromosome 12.

Genomic imprinting causes parental origin-specific gene expression. Cis-acting regulatory elements that control imprinting are not fully understood but involve regions that become differentially methylated on the two parental chromosomes during male and female gametogenesis. Understanding properties of maternally and paternally inherited imprints provides insight into the mechanisms and evolution of genomic imprinting. Previously we identified an intergenic germline-derived differentially methylated region (IG-DMR) that is a candidate control element for an imprinted domain on distal mouse chromosome 12 (ref. 5). The 1-Mb cluster contains the paternally expressed protein-coding genes Dlk1 (refs. 6,7) and Dio3 (ref. 8,9) and several maternally expressed non-coding RNAs, including Gtl2 (refs. 6,7,10) and C/D snoRNAs. A retrotransposon-like gene (Rtl1) is expressed from the paternal chromosome and has an antisense transcript expressed from the maternal chromosome containing two microRNAs with full complementarity to Rtl1 (ref. 12). Here we show that deletion of the IG-DMR from the maternally inherited chromosome causes bidirectional loss of imprinting of all genes in the cluster. When the deletion is transmitted from the father, imprinting is unaltered. These results prove that the IG-DMR is a control element for all imprinted genes on the maternal chromosome only and indicate that the two parental chromosomes control allele-specific gene expression differently.

The potential role of gene duplications in the evolution of imprinting mechanisms.

Using the completed genomic sequences of mouse and human we performed a comparative analyses of imprinted genes and gene clusters. For many imprinted genes we could detect imprinted as well as non-imprinted paralogues. The inter- and intrachromosomal similarities between paralogues and their linkage to imprinting clusters suggests that imprinted genes were dispersed throughout the genome by gene duplications as well as translocation and transposition events. Our findings indicate that imprinting clusters may have been linked together on one (or a few) ancestral pre-imprinted chromosome(s), arguing for a common mechanistic origin of imprinting control. Imprinting may originally have evolved on a simple basis of dosage compensation required for some duplicated genes (chromosomes) followed by selection of sex-biased expression control.

Imprinted microRNA genes transcribed antisense to a reciprocally imprinted retrotransposon-like gene.

MicroRNAs (miRNAs) are an abundant class of RNAs that are approximately 21-25 nucleotides (nt) long, interact with mRNAs and trigger either translation repression or RNA cleavage (RNA interference, RNAi) depending on the degree of complementarity with their targets. Here we show that the imprinted mouse distal chromosome 12 locus encodes two miRNA genes expressed from the maternally inherited chromosome and antisense to a retrotransposon-like gene (Rtl1) expressed only from the paternal allele.