DNA methylation profiling of pseudogene-parental gene pairs and two gene families

Divergent DNA Methylation Provides Insights into the Evolution of Duplicate Genes in Zebrafish

The evolutionary mechanism, fate and function of duplicate genes in various taxa have been widely studied; however, the mechanism underlying the maintenance and divergence of duplicate genes in Danio rerio remains largely unexplored. Whether and how the divergence of DNA methylation between duplicate pairs is associated with gene expression and evolutionary time are poorly understood. In this study, by analyzing bisulfite sequencing (BS-seq) and RNA-seq datasets from public data, we demonstrated that DNA methylation played a critical role in duplicate gene evolution in zebrafish. Initially, we found promoter methylation of duplicate genes generally decreased with evolutionary time as measured by synonymous substitution rate between paralogous duplicates (Ks). Importantly, promoter methylation of duplicate genes was negatively correlated with gene expression. Interestingly, for 665 duplicate gene pairs, one gene was consistently promoter methylated, while the other was unmethylated across nine different datasets we studied. Moreover, one motif enriched in promoter methylated duplicate genes tended to be bound by the transcription repression factor FOXD3, whereas a motif enriched in the promoter unmethylated sequences interacted with the transcription activator Sp1, indicating a complex interaction between the genomic environment and epigenome. Besides, body-methylated genes showed longer length than body-unmethylated genes. Overall, our results suggest that DNA methylation is highly important in the differential expression and evolution of duplicate genes in zebrafish.

DNA methylation of the LIN28 pseudogene family

Background

DNA methylation directs the epigenetic silencing of selected regions of DNA, including the regulation of pseudogenes, and is widespread throughout the genome. Pseudogenes are decayed copies of duplicated genes that have spread throughout the genome by transposition. Pseudogenes are transcriptionally silenced by DNA methylation, but little is known about how pseudogenes are targeted for methylation or how methylation levels are maintained in different tissues.

Results

We employed bisulfite next generation sequencing to examine the methylation status of the LIN28 gene and four processed pseudogenes derived from LIN28. The objective was to determine whether LIN28 pseudogenes maintain the same pattern of methylation as the parental gene or acquire a methylation pattern independent of the gene of origin. In this study, we determined that the methylation status of LIN28 pseudogenes does not resemble the pattern evident for the LIN28 gene, but rather these pseudogenes appear to acquire methylation patterns independent of the parental gene. Furthermore, we observed that methylation levels of the examined pseudogenes correlate to the location of insertion within the genome. LIN28 pseudogenes inserted into gene bodies were highly methylated in all tissues examined. In contrast, pseudogenes inserted into genomic regions that are not proximal to genes were differentially methylated in various tissue types.

Conclusions

Our analysis suggests that Lin28 pseudogenes do not aquire patterns of tissue-specific methylation as for the parental gene, but rather are methylated in patterns specific to the local genomic environment into which they were inserted.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1487-3) contains supplementary material, which is available to authorized users.

Discrimination of pseudogene and parental gene DNA methylation using allelic bisulfite sequencing

Determining the methylation status of genes with pseudogenes can be technically challenging due to sequence homology. High sequence homology can result in the amplification of both pseudogene and parental gene alleles, potentially leading to data misinterpretation. Allelic bisulfite sequencing allows for detection of the methylation status of individual alleles at nucleotide resolution and represents the most reliable method for discriminating pseudogene and parental gene sequences. Here, we discuss important points that should be considered when investigating pseudogene and parental gene methylation status and we describe the method of allelic bisulfite sequencing, including assay design.

Direct bisulfite sequencing for examination of DNA methylation with gene and nucleotide resolution from brain tissues

DNA methylation is an epigenetic modification that is essential for the development and mature function of the central nervous system. Due to the relevance of this modification to the transcriptional control of gene expression, it is often necessary to examine changes in DNA methylation patterns with both gene and single-nucleotide resolution. Here, we describe an in-depth basic protocol for direct bisulfite sequencing of DNA isolated from brain tissue, which will permit direct assessment of methylation status at individual genes as well as individual cytosine molecules/nucleotides within a genomic region. This method yields analysis of DNA methylation patterns that is robust, accurate, and reproducible, thereby allowing insights into the role of alterations in DNA methylation in brain tissue.

Early de novo DNA methylation and prolonged demethylation in the muscle lineage

Myogenic cell cultures derived from muscle biopsies are excellent models for human cell differentiation. We report the first comprehensive analysis of myogenesis-specific DNA hyper- and hypo-methylation throughout the genome for human muscle progenitor cells (both myoblasts and myotubes) and skeletal muscle tissue vs. 30 non-muscle samples using reduced representation bisulfite sequencing. We also focused on four genes with extensive hyper- or hypo-methylation in the muscle lineage (PAX3, TBX1, MYH7B/MIR499 and OBSCN) to compare DNA methylation, DNaseI hypersensitivity, histone modification, and CTCF binding profiles. We found that myogenic hypermethylation was strongly associated with homeobox or T-box genes and muscle hypomethylation with contractile fiber genes. Nonetheless, there was no simple relationship between differential gene expression and myogenic differential methylation, rather only for subsets of these genes, such as contractile fiber genes. Skeletal muscle retained ~30% of the hypomethylated sites but only ~3% of hypermethylated sites seen in myogenic progenitor cells. By enzymatic assays, skeletal muscle was 2-fold enriched globally in genomic 5-hydroxymethylcytosine (5-hmC) vs. myoblasts or myotubes and was the only sample type enriched in 5-hmC at tested myogenic hypermethylated sites in PAX3/CCDC140 andTBX1. TET1 and TET2 RNAs, which are involved in generation of 5-hmC and DNA demethylation, were strongly upregulated in myoblasts and myotubes. Our findings implicate de novo methylation predominantly before the myoblast stage and demethylation before and after the myotube stage in control of transcription and co-transcriptional RNA processing. They also suggest that, in muscle, TET1 or TET2 are involved in active demethylation and in formation of stable 5-hmC residues.

The influence of DNA sequence on epigenome-induced pathologies

Clear cause-and-effect relationships are commonly established between genotype and the inherited risk of acquiring human and plant diseases and aberrant phenotypes. By contrast, few such cause-and-effect relationships are established linking a chromatin structure (that is, the epitype) with the transgenerational risk of acquiring a disease or abnormal phenotype. It is not entirely clear how epitypes are inherited from parent to offspring as populations evolve, even though epigenetics is proposed to be fundamental to evolution and the likelihood of acquiring many diseases. This article explores the hypothesis that, for transgenerationally inherited chromatin structures, "genotype predisposes epitype", and that epitype functions as a modifier of gene expression within the classical central dogma of molecular biology. Evidence for the causal contribution of genotype to inherited epitypes and epigenetic risk comes primarily from two different kinds of studies discussed herein. The first and direct method of research proceeds by the examination of the transgenerational inheritance of epitype and the penetrance of phenotype among genetically related individuals. The second approach identifies epitypes that are duplicated (as DNA sequences are duplicated) and evolutionarily conserved among repeated patterns in the DNA sequence. The body of this article summarizes particularly robust examples of these studies from humans, mice, Arabidopsis, and other organisms. The bulk of the data from both areas of research support the hypothesis that genotypes predispose the likelihood of displaying various epitypes, but for only a few classes of epitype. This analysis suggests that renewed efforts are needed in identifying polymorphic DNA sequences that determine variable nucleosome positioning and DNA methylation as the primary cause of inherited epigenome-induced pathologies. By contrast, there is very little evidence that DNA sequence directly determines the inherited positioning of numerous and diverse post-translational modifications of histone side chains within nucleosomes. We discuss the medical and scientific implications of these observations on future research and on the development of solutions to epigenetically induced disorders.

DNA methylation rebalances gene dosage after mammalian gene duplications

Although gene duplication plays a major role in organismal evolution, it may also lead to gene dosage imbalance, thereby having an immediate adverse effect on an organism's fitness. Investigating the evolution of the expression patterns of genes that duplicated after the divergence of rodents and primates, we confirm that adaptive evolution has been involved in dosage rebalance after gene duplication. To understand mechanisms underlying this process, we examined 1) microRNA (miRNA)-mediated gene regulation, 2) cis-regulatory sequence modifications, and 3) DNA methylation. Neither miRNA-mediated regulation nor cis-regulatory changes was found to be associated with expression reduction of duplicate genes. By contrast, duplicate genes, especially lowly expressed copies, were heavily methylated in the upstream region. However, for duplicate genes encoding proteins that are members of macromolecular complexes, heavy methylation in the genic region was not consistently observed. This result held after controlling potential confounding factors, such as enrichment in functional categories. Our results suggest that during mammalian evolution, DNA methylation plays a dominant role in dosage rebalance after gene duplication by inhibiting transcription initiation of duplicate genes.

Sequence overlap between autosomal and sex-linked probes on the Illumina HumanMethylation27 microarray

The Illumina Infinium HumanMethylation27 BeadChip (Illumina 27k) microarray is a high-throughput platform capable of interrogating the human DNA methylome. In a search for autosomal sex-specific DNA methylation using this microarray, we discovered autosomal CpG loci showing significant methylation differences between the sexes. However, we found that the majority of these probes cross-reacted with sequences from sex chromosomes. Moreover, we determined that 6-10% of the microarray probes are non-specific and map to highly homologous genomic sequences. Using probes targeting different CpGs that are exact duplicates of each other, we investigated the precision of these repeat measurements and concluded that the overall precision of this microarray is excellent. In addition, we identified a small number of probes targeting CpGs that include single-nucleotide polymorphisms. Overall, our findings address several technical issues associated with the Illumina 27k microarray that, once considered, will enhance the analysis and interpretation of data generated from this platform.