pENCODE: a plant encyclopedia of DNA elements

ENCODE projects exist for many eukaryotes, including humans, but as of yet no defined project exists for plants. A plant ENCODE would be invaluable to the research community and could be more readily produced than its metazoan equivalents by capitalizing on the preexisting infrastructure provided from similar projects. Collecting and normalizing plant epigenomic data for a range of species will facilitate hypothesis generation, cross-species comparisons, annotation of genomes, and an understanding of epigenomic functions throughout plant evolution. Here, we discuss the need for such a project, outline the challenges it faces, and suggest ways forward to build a plant ENCODE.

Comparative epigenomics reveals evolution of duplicated genes in potato and tomato

The evolution of duplicated genes after polyploidization has been the subject of many evolutionary biology studies. Potato (Solanum tuberosum) and tomato (Solanum lycopersicum) are the first two sequenced genomes of asterids, and share a common polyploidization event. However, the epigenetic role of DNA methylation on the evolution of duplicated genes derived from polyploidization is not fully understood. Here, we explore the role of the DNA methylation in the evolution of duplicated genes in potato and tomato. The overall levels of DNA methylation are different, although patterns of DNA methylation are similar in potato and tomato. Different types of duplicated genes can display different methylation patterns in potato and tomato. In addition, we found that differences in the methylation levels between duplicated genes were associated with gene expression divergence. In particular, for the majority of duplicated gene pairs, one copy is always hyper- or hypo-methylated compared with the other copy across different tomato fruit ripening stages, and these genes are enriched for specific function related to transcription factor (TF) activity. Furthermore, transcription of hundreds of duplicated TFs was shown to be regulated by DNA methylation during fruit ripening stages in tomato, some of which are well-known fruit ripening TFs. Taken together, our results support the notion that DNA methylation may facilitate divergent evolution of duplicated genes and play roles in important biological processes such as tomato fruit ripening.

Comparative Epigenomic Profiling of the DNA Methylome in Mouse and Zebrafish Uncovers High Interspecies Divergence

The DNA methylation landscape is dynamically patterned during development and distinct methylation patterns distinguish healthy from diseased cells. However, whether tissue-specific methylation patterns are conserved across species is not known. We used comparative methylome analysis of base-resolution DNA methylation profiles from the liver and brain of mouse and zebrafish generated by reduced representation bisulfite sequencing to identify the conserved and divergent aspects of the methylome in these commonly used vertebrate model organisms. On average, 24% of CpGs are methylated in mouse livers and the pattern of methylation was highly concordant among four male mice from two different strains. The same level of methylation (24.2%) was identified in mouse brain. In striking contrast, zebrafish had 63 and 70% of CpG methylation in the liver and brain, respectively. This is attributed, in part, to the higher percentage of the zebrafish genome occupied by transposable elements (52% vs. 45% in mice). Thus, the species identity was more significant in determining methylome patterning than was the similarity in organ function. Conserved features of the methylome across tissues and species was the exclusion of methylation from promoters and from CpG islands near transcription start sites, and the clustering of methylated CpGs in gene bodies and intragenic regions. These data suggest that DNA methylation reflects species-specific genome structure, and supports the notion that DNA methylation in non-promoter regions may contribute to genome evolution.

Comparative Epigenomics Reveals that RNA Polymerase II Pausing and Chromatin Domain Organization Control Nematode piRNA Biogenesis

Piwi-interacting RNAs (piRNAs) are important for genome regulation across metazoans, but their biogenesis evolves rapidly. In Caenorhabditis elegans, piRNA loci are clustered within two 3-Mb regions on chromosome IV. Each piRNA locus possesses an upstream motif that recruits RNA polymerase II to produce an ∼28 nt primary transcript. We used comparative epigenomics across nematodes to gain insight into the origin, evolution, and mechanism of nematode piRNA biogenesis. We show that the piRNA upstream motif is derived from core promoter elements controlling snRNA transcription. We describe two alternative modes of piRNA organization in nematodes: in C. elegans and closely related nematodes, piRNAs are clustered within repressive H3K27me3 chromatin, while in other species, typified by Pristionchus pacificus, piRNAs are found within introns of active genes. Additionally, we discover that piRNA production depends on sequence signals associated with RNA polymerase II pausing. We show that pausing signals synergize with chromatin to control piRNA transcription.

•

Nematode piRNA transcription evolved from small nuclear RNA biogenesis

•

Clustered piRNAs are produced from regulated (H3K27me3) chromatin domains

•

Dispersed piRNAs are produced from active (H3K36me3) chromatin domains

•

RNA polymerase II pausing determines the short (∼28 nt) length of piRNA precursors

Genome-wide mapping of histone modifications during axenic growth in two species of Leptosphaeria maculans showing contrasting genomic organization

Leptosphaeria maculans 'brassicae' (Lmb) and Leptosphaeria maculans 'lepidii' (Lml) are closely related phytopathogenic species that exhibit a large macrosynteny but contrasting genome structure. Lmb has more than 30% of repeats clustered in large repeat-rich regions, while the Lml genome has only a small amount of evenly distributed repeats. Repeat-rich regions of Lmb are enriched in effector genes, expressed during plant infection. The distinct genome structures of Lmb and Lml provide an excellent model for comparing the organization of pathogenicity genes in relation to the chromatin landscape in two closely related phytopathogenic fungi. Here, we performed chromatin immunoprecipitation (ChIP) during axenic culture, targeting histone modifications typical for heterochromatin or euchromatin, combined with transcriptomic analysis to analyze the influence of chromatin organization on gene expression. In both species, we found that facultative heterochromatin is enriched with genes lacking functional annotation, including numerous effector and species-specific genes. Notably, orthologous genes located in H3K27me3 domains are enriched with effector genes. Compared to other fungal species, including Lml, Lmb is distinct in having large H3K9me3 domains associated with repeat-rich regions that contain numerous species-specific effector genes. Discovery of these two distinctive heterochromatin landscapes now raises questions about their involvement in the regulation of pathogenicity, the dynamics of these domains during plant infection and the selective advantage to the fungus to host effector genes in H3K9me3 or H3K27me3 domains.

Identifying co-opted transposable elements using comparative epigenomics

The human genome gives rise to different epigenomic landscapes that define each cell type and can be deregulated in disease. Recent efforts by ENCODE, the NIH Roadmap and the International Human Epigenome Consortium (IHEC) have made significant advances towards assembling reference epigenomic maps of various tissues. Notably, these projects have found that approximately 80% of human DNA was biochemically active in at least one epigenomic assay while only approximately 10% of the sequence displayed signs of purifying selection. Given that transposable elements (TEs) make up at least 50% of the human genome and can be actively transcribed or act as regulatory elements either for their own purposes or be co-opted for the benefit of their host; we are interested in exploring their overall contribution to the "functional" genome. Traditional methods used to identify functional DNA have relied on comparative genomics, conservation analysis and low throughput validation assays. To discover co-opted TEs, and distinguish them from noisy genomic elements, we argue that comparative epigenomic methods will also be important.

Predicted Archaic 3D Genome Organization Reveals Genes Related to Head and Spinal Cord Separating Modern from Archaic Humans

High coverage sequences of archaic humans enabled the reconstruction of their DNA methylation patterns. This allowed comparing gene regulation between human groups, and linking such regulatory changes to phenotypic differences. In a previous work, a detailed comparison of DNA methylation in modern humans, archaic humans, and chimpanzees revealed 873 modern human-derived differentially methylated regions (DMRs). To understand the regulatory implications of these DMRs, we defined differentially methylated genes (DMGs) as genes that harbor DMRs in their promoter or gene body. While most of the modern human-derived DMRs could be linked to DMGs, many others remained unassigned. Here, we used information on 3D genome organization to link ~70 out of the remaining 288 unassigned DMRs to genes. Combined with the previously identified DMGs, we reinforce the enrichment of these genes with vocal and facial anatomy, and additionally find significant enrichment with the spinal column, chin, hair, and scalp. These results reveal the importance of 3D genomic organization in understanding gene regulation by DNA methylation.

Toward understanding the evolution of vertebrate gene regulatory networks: comparative genomics and epigenomic approaches

Vertebrates, as most animal phyla, originated >500 million years ago during the Cambrian explosion, and progressively radiated into the extant classes. Inferring the evolutionary history of the group requires understanding the architecture of the developmental programs that constrain the vertebrate anatomy. Here, I review recent comparative genomic and epigenomic studies, based on ChIP-seq and chromatin accessibility, which focus on the identification of functionally equivalent cis-regulatory modules among species. This pioneer work, primarily centered in the mammalian lineage, has set the groundwork for further studies in representative vertebrate and chordate species. Mapping of active regulatory regions across lineages will shed new light on the evolutionary forces stabilizing ancestral developmental programs, as well as allowing their variation to sustain morphological adaptations on the inherited vertebrate body plan.

Comparative Analyses of Sperm DNA Methylomes Among Three Commercial Pig Breeds Reveal Vital Hypomethylated Regions Associated With Spermatogenesis and Embryonic Development

Identifying epigenetic changes is essential for an in-depth understanding of phenotypic diversity and pigs as the human medical model for anatomizing complex diseases. Abnormal sperm DNA methylation can lead to male infertility, fetal development failure, and affect the phenotypic traits of offspring. However, the whole genome epigenome map in pig sperm is lacking to date. In this study, we profiled methylation levels of cytosine in three commercial pig breeds, Landrace, Duroc, and Large White using whole-genome bisulfite sequencing (WGBS). The results showed that the correlation of methylation levels between Landrace and Large White pigs was higher. We found that 1,040-1,666 breed-specific hypomethylated regions (HMRs) were associated with embryonic developmental and economically complex traits for each breed. By integrating reduced representation bisulfite sequencing (RRBS) public data of pig testis, 1743 conservated HMRs between sperm and testis were defined, which may play a role in spermatogenesis. In addition, we found that the DNA methylation patterns of human and pig sperm showed high similarity by integrating public data from WGBS and chromatin immunoprecipitation sequencing (ChIP-seq) in other mammals, such as human and mouse. We identified 2,733 conserved HMRs between human and pig involved in organ development and brain-related traits, such as NLGN1 (neuroligin 1) containing a conserved-HMR between human and pig. Our results revealed the similarities and diversity of sperm methylation patterns among three commercial pig breeds and between human and pig. These findings are beneficial for elucidating the mechanism of male fertility, and the changes in commercial traits that undergo strong selection.

Comparative Epigenomics Reveals Host Diversity of the Trichinella Epigenomes and Their Effects on Differential Parasitism

Comparative epigenomics provides new insights on evolutionary biology in relation with complex interactions between species and their environments. In the present study, we focus on deciphering the conservation and divergence of DNA methylomes during Trichinella evolution. Whole-genome bisulfite sequencing and RNA-seq were performed on the two clades of Trichinella species, in addition to whole-genome sequencing. We demonstrate that methylation patterns of sing-copy orthologous genes (SCOs) of the 12 Trichinella species are host-related and can mirror known phylogenetic relationships. Among these SCOs, we identify a panel of genes exhibiting hyper-/hypo-methylated features in gene-bodies or respective promoters that play pivotal roles in transcriptome regulation. These hyper-/hypo-methylated SCOs are also of functional significance across developmental stages, as they are highly enriched species-specific and stage-specific expressed genes both in Ad and ML stages. We further identify a set of parasitism-related functional genes that exhibit host-related differential methylation and expression among those SCOs, including p53-like transcription factor and Cdc37 that are of functional significance for elucidating differential parasitology between the two clades of Trichinella. This comparative epigenome study can help to decipher the environmental effects on differential adaptation and parasitism of the genus Trichinella.

Stability of DNA methylation and chromatin accessibility in structurally diverse maize genomes

Accessible chromatin and unmethylated DNA are associated with many genes and cis-regulatory elements. Attempts to understand natural variation for accessible chromatin regions (ACRs) and unmethylated regions (UMRs) often rely upon alignments to a single reference genome. This limits the ability to assess regions that are absent in the reference genome assembly and monitor how nearby structural variants influence variation in chromatin state. In this study, de novo genome assemblies for four maize inbreds (B73, Mo17, Oh43, and W22) are utilized to assess chromatin accessibility and DNA methylation patterns in a pan-genome context. A more complete set of UMRs and ACRs can be identified when chromatin data are aligned to the matched genome rather than a single reference genome. While there are UMRs and ACRs present within genomic regions that are not shared between genotypes, these features are 6- to 12-fold enriched within regions between genomes. Characterization of UMRs present within shared genomic regions reveals that most UMRs maintain the unmethylated state in other genotypes with only ∼5% being polymorphic between genotypes. However, the majority (71%) of UMRs that are shared between genotypes only exhibit partial overlaps suggesting that the boundaries between methylated and unmethylated DNA are dynamic. This instability is not solely due to sequence variation as these partially overlapping UMRs are frequently found within genomic regions that lack sequence variation. The ability to compare chromatin properties among individuals with structural variation enables pan-epigenome analyses to study the sources of variation for accessible chromatin and unmethylated DNA.

Examining the Utility of the Mammalian Methylation Array for Pan-Mammalian Analysis of Monozygotic Twinning

BACKGROUND/OBJECTIVES

Human identical twins are born at a rate of 3-4 per 1000 live births. Many other mammals also occasionally produce monozygotic twins, referred to as sporadic polyembryony. The underlying mechanisms are unknown. Through epigenome-wide association studies (EWAS), we identified a robust DNA methylation signature in somatic tissues from human monozygotic (MZ) twins, comprising 834 differentially methylated positions (MZ-DMPs). The results point to a connection between monozygotic twinning and early genome programming and enable new angles to study monozygotic twinning.

METHODS

The mammalian methylation array (MMA) measures 38,608 CpGs focusing on regions that are well-conserved across many mammalian species, allowing for pan-mammalian comparative epigenomic studies. Here, we successfully map human MZ-DMPs to probes of the mammalian methylation array across 157 mammalian genomes.

RESULTS

As expected, based on the modest probe overlap between Illumina 450k/EPIC and mammalian methylation array probes, only a subset of MZ-DMPs reside in conserved regions covered by the mammalian methylation array. These include probes mapping to NPAS3, KLHL35, CASZ1, and ATP2B2. Re-analysis restricting the original EWAS in humans to conserved MMA regions yielded additional MZ-DMPs, suggesting that more loci may be detected by application of the mammalian array to monozygotic twins.

CONCLUSIONS

In conclusion, the mammalian methylation array may prove to be a promising platform to study whether a shared DNA methylation signature of sporadic polyembryony exists across diverse mammalian species. This may potentially point to shared underlying mechanisms.

Comparative analysis of sperm DNA methylation supports evolutionary acquired epigenetic plasticity for organ speciation

Aim: To perform a comparative epigenomic analysis of DNA methylation in spermatozoa from humans, mice, rats and mini-pigs. Materials & methods: Genome-wide DNA methylation analysis was used to compare the methylation profiles of orthologous CpG sites. Transcription profiles of early embryo development were analyzed to provide insight into the association between sperm methylation and gene expression programming. Results: We identified DNA methylation variation near genes related to the central nervous system and signal transduction. Gene expression dynamics at different time points of preimplantation stages were modestly associated with spermatozoal DNA methylation at the nearest promoters. Conclusion: Conserved genomic regions subject to epigenetic variation across different species were associated with specific organ functions, suggesting their potential contribution to organ speciation and long-term adaptation to the environment.