DNA Demethylation Dynamics in the Human Prenatal Germline.

Global DNA demethylation in humans is a fundamental process that occurs in pre-implantation embryos and reversion to naive ground state pluripotent stem cells (PSCs). However, the extent of DNA methylation reprogramming in human germline cells is unknown. Here, we performed whole-genome bisulfite sequencing (WGBS) and RNA-sequencing (RNA-seq) of human prenatal germline cells from 53 to 137 days of development. We discovered that the transcriptome and methylome of human germline is distinct from both human PSCs and the inner cell mass (ICM) of human blastocysts. Using this resource to monitor the outcome of global DNA demethylation with reversion of primed PSCs to the naive ground state, we uncovered hotspots of ultralow methylation at transposons that are protected from demethylation in the germline and ICM. Taken together, the human germline serves as a valuable in vivo tool for monitoring the epigenome of cells that have emerged from a global DNA demethylation event.

Epigenetic factors direct synergistic and antagonistic regulation of transposable elements in Arabidopsis.

Arabidopsis (Arabidopsis thaliana) HISTONE DEACETYLASE 6 (HDA6) and HISTONE DEMETHYLASES LSD-LIKE 1 (LDL1) and LDL2 synergistically regulate the expression of long non-coding RNAs associated with H3Ac and H3K4me2. The underlying mechanisms of such highly coordinated interactions among genetic and epigenetic factors contributing to this collaborative regulation remain largely unclear. We analyzed all transposable elements (TEs) across the Arabidopsis genome and the individual and combined roles of HDA6 and LDL1/LDL2 by dissecting multi-layered epigenomes and their association with transcription. Instead of an individual synergistic effect, we observed dual synergistic and antagonistic effects, which are positively associated with H3Ac and H3K4me2 while maintaining a negative but moderate association with DNA methylation. Specifically, two modes of synergistic regulation were discovered in TEs: 74% are primarily regulated by HDA6, with less dependence on LDL1/LDL2, and the remaining 26% are co-regulated by both. Between the two modes, we showed that HDA6 has a strong effect on TE silencing, whereas LDL1/LDL2 plays a weaker yet crucial role in co-regulation with HDA6. Our results led to a model of epigenomic regulation - the differential de-repression between the two modes of synergistic regulation of TEs was determined by H3Ac and H3K4me2 levels, where TEs are in accessible chromatins free of DNA methylation, and this open chromatin environment precedes transcriptional changes and epigenome patterning. Our results discovered unbalanced effects of genetic factors in synergistic regulation through delicately coordinated multi-layered epigenomes and chromatin accessibility.

Perturbation of periodic spot-generation balance leads to diversified pigmentation patterning of harlequin Phalaenopsis orchids: in silico prediction.

BACKGROUND: A retrotransposon HORT1 in the promoter of the anthocyanin activator gene PeMYB11, microRNA858 (miR858) that targets PeMYB11, and a repressor PeMYBx have been implicated in pigmentation patterning diversity of harlequin Phalaenopsis orchids. However, the interrelationship among them remains to be elucidated. RESULTS: To understand how these factors interact to generate anthocyanin spots in Phalaenopsis, we successfully developed a mathematical model based on the known reaction-diffusion system to simulate their interplay and refined the conceptual biological model. Intriguingly, the expression of both PeMYBx and PeMYB11 were in phase for purple spot formation, even though they showed adverse effects on anthocyanin accumulations. An increase in the self-activation rate of PeMYB11 resulted in the increased size of purple spots, but no effects on spot fusion. Decreased degradation rate of miR858 in the purple regions, led to disruption of the formation of spotted pigmentation patterning and a full-red pigmentation pattern. Significantly, the reduced miR858 level promotes the fusion of large dark purple dots induced by the solo-LTR of HORT1, eventually generating the purple patches. In addition, the spatially heterogeneous insertion of HORT1 caused by the remnant solo-LTR of HORT1 derived from random homologous unequal recombination of HORT1 in individual cells of floral organs could explain the diverse pigmentation patterning of harlequin Phalaenopsis. CONCLUSIONS: This devised model explains how HORT1 and miR858 regulate the formation of the pigmentation patterning and holds great promise for developing efficient and innovative approaches to breeding harlequin Phalaenopsis orchids.

Rice transformation treatments leave specific epigenome changes beyond tissue culture.

During transgenic plant production, tissue culture often carries epigenetic, and genetic changes that underlie somaclonal variations, leading to unpredictable phenotypes. Additionally, specific treatments for rice (Oryza sativa) transformation processes may individually or jointly contribute to somaclonal variations, but their specific impacts on rice epigenomes toward transcriptional variations remain unknown. Here, the impact of individual transformation treatments on genome-wide DNA methylation and the transcriptome were examined. In addition to activating stress-responsive genes, individual transformation components targeted different gene expression modules that were enriched in specific functional categories. The transformation treatments strongly impacted DNA methylation and expression; 75% were independent of tissue culture. Furthermore, our genome-wide analysis showed that the transformation treatments consistently resulted in global hypo-CHH methylation enriched at promoters highly associated with downregulation, particularly when the promoters were colocalized with miniature inverted-repeat transposable elements. Our results clearly highlight the specificity of impacts triggered by individual transformation treatments during rice transformation with the potential association between DNA methylation and gene expression. These changes in gene expression and DNA methylation resulting from rice transformation treatments explain a significant portion of somaclonal variations, that is, way beyond the tissue culture effect.

HISTONE DEACETYLASE 15 and MOS4-associated complex subunits 3A/3B coregulate intron retention of ABA-responsive genes.

Histone deacetylases (HDAs) play an important role in transcriptional regulation of multiple biological processes. In this study, we investigated the function of HDA15 in abscisic acid (ABA) responses. We used immunopurification coupled with mass spectrometry-based proteomics to identify proteins interacting with HDA15 in Arabidopsis (Arabidopsis thaliana). HDA15 interacted with the core subunits of the MOS4-associated complex (MAC), MAC3A and MAC3B, with interaction between HDA15 and MAC3B enhanced by ABA. hda15 and mac3a/mac3b mutants were ABA-insensitive during seed germination and hyposensitive to salinity. RNA sequencing analysis demonstrated that HDA15 and MAC3A/MAC3B co-regulate ABA-responsive intron retention (IR). Furthermore, HDA15 reduced the histone acetylation level of genomic regions near ABA-responsive IR sites and the association of MAC3B with ABA-responsive pre-mRNA was dependent on HDA15. Our results indicate that HDA15 is involved in ABA responses by interacting with MAC3A/MAC3B to mediate splicing of introns.

Mycena genomes resolve the evolution of fungal bioluminescence.

Mushroom-forming fungi in the order Agaricales represent an independent origin of bioluminescence in the tree of life; yet the diversity, evolutionary history, and timing of the origin of fungal luciferases remain elusive. We sequenced the genomes and transcriptomes of five bonnet mushroom species (Mycena spp.), a diverse lineage comprising the majority of bioluminescent fungi. Two species with haploid genome assemblies ∼150 Mb are among the largest in Agaricales, and we found that a variety of repeats between Mycena species were differentially mediated by DNA methylation. We show that bioluminescence evolved in the last common ancestor of mycenoid and the marasmioid clade of Agaricales and was maintained through at least 160 million years of evolution. Analyses of synteny across genomes of bioluminescent species resolved how the luciferase cluster was derived by duplication and translocation, frequently rearranged and lost in most Mycena species, but conserved in the Armillaria lineage. Luciferase cluster members were coexpressed across developmental stages, with the highest expression in fruiting body caps and stipes, suggesting fruiting-related adaptive functions. Our results contribute to understanding a de novo origin of bioluminescence and the corresponding gene cluster in a diverse group of enigmatic fungal species.

SWI3B and HDA6 interact and are required for transposon silencing in Arabidopsis.

Although the interplay of covalent histone acetylation/deacetylation and ATP-dependent chromatin remodelling is crucial for the regulation of chromatin structure and gene expression in eukaryotes, the underlying molecular mechanism in plants remains largely unclear. Here we show a direct interaction between Arabidopsis SWI3B, an essential subunit of the SWI/SNF chromatin-remodelling complex, and the RPD3/HDA1-type histone deacetylase HDA6 both in vitro and in vivo. Furthermore, SWI3B and HDA6 co-repress the transcription of a subset of transposons. Both SWI3B and HDA6 maintain transposon silencing by decreasing histone H3 lysine 9 acetylation, but increasing histone H3 lysine 9 di-methylation, DNA methylation and nucleosome occupancy. Our findings reveal that SWI3B and HDA6 may act in the same co-repressor complex to maintain transposon silencing in Arabidopsis.

Global impacts of chromosomal imbalance on gene expression in Arabidopsis and other taxa.

Changes in dosage of part of the genome (aneuploidy) have long been known to produce much more severe phenotypic consequences than changes in the number of whole genomes (ploidy). To examine the basis of these differences, global gene expression in mature leaf tissue for all five trisomies and in diploids, triploids, and tetraploids of Arabidopsis thaliana was studied. The trisomies displayed a greater spread of expression modulation than the ploidy series. In general, expression of genes on the varied chromosome ranged from compensation to dosage effect, whereas genes from the remainder of the genome ranged from no effect to reduced expression approaching the inverse level of chromosomal imbalance (2/3). Genome-wide DNA methylation was examined in each genotype and found to shift most prominently with trisomy 4 but otherwise exhibited little change, indicating that genetic imbalance is generally mechanistically unrelated to DNA methylation. Independent analysis of gene functional classes demonstrated that ribosomal, proteasomal, and gene body methylated genes were less modulated compared with all classes of genes, whereas transcription factors, signal transduction components, and organelle-targeted protein genes were more tightly inversely affected. Comparing transcription factors and their targets in the trisomies and in expression networks revealed considerable discordance, illustrating that altered regulatory stoichiometry is a major contributor to genetic imbalance. Reanalysis of published data on gene expression in disomic yeast and trisomic mouse cells detected similar stoichiometric effects across broad phylogenetic taxa, and indicated that these effects reflect normal gene regulatory processes.

Epigenomic regulation of OTU5 in Arabidopsis thaliana.

Epigenetic regulation by DNA methylation and histone marks is crucial to plant development. In Arabidopsis, the otu5 mutant exhibited altered root phenotypes resembling those of phosphate-deficient plants. In low phosphate (Pi) conditions, altered H3K4 and H3K27 trimethylation were associated with the expression of Pi homeostasis-related genes. However, the genetic effect of OTU5 on the epigenomes was left unexplored. We assessed genome-wide DNA methylation, gene expression and histone modifications of roots from both Col-0 and otu5 mutants. We found that OTU5 altered DNA methylation profile with a context-specific effect through targeting local genomic regions. Our analysis showed that in otu5 the abundance of H3K4me3 was clearly associated with the changes of DNA methylation, leading to the transcriptional difference from wildtype. We concluded that OTU5 induced cross-talks among epigenomes that altogether impacted the regulation of approximately 7060 genes. Of which 186 genes associated with root development were likely to be epigenetically regulated.

MethGET: web-based bioinformatics software for correlating genome-wide DNA methylation and gene expression.

BACKGROUND: DNA methylation is a major epigenetic modification involved in regulating gene expression. The effects of DNA methylation on gene expression differ by genomic location and vary across kingdoms, species and environmental conditions. To identify the functional regulatory roles of DNA methylation, the correlation between DNA methylation changes and alterations in gene expression is crucial. With the advance of next-generation sequencing, genome-wide methylation and gene expression profiling have become feasible. Current bioinformatics tools for investigating such correlation are designed to the assessment of DNA methylation at CG sites. The correlation of non-CG methylation and gene expression is very limited. Some bioinformatics databases allow correlation analysis, but they are limited to specific genomes such as that of humans and do not allow user-provided data. RESULTS: Here, we developed a bioinformatics web tool, MethGET (Methylation and Gene Expression Teller), that is specialized to analyse the association between genome-wide DNA methylation and gene expression. MethGET is the first web tool to which users can supply their own data from any genome. It is also the tool that correlates gene expression with CG, CHG, and CHH methylation based on whole-genome bisulfite sequencing data. MethGET not only reveals the correlation within an individual sample (single-methylome) but also performs comparisons between two groups of samples (multiple-methylomes). For single-methylome analyses, MethGET provides Pearson correlations and ordinal associations to investigate the relationship between DNA methylation and gene expression. It also groups genes with different gene expression levels to view the methylation distribution at specific genomic regions. Multiple-methylome analyses include comparative analyses and heatmap representations between two groups. These functions enable the detailed investigation of the role of DNA methylation in gene regulation. Additionally, we applied MethGET to rice regeneration data and discovered that CHH methylation in the gene body region may play a role in the tissue culture process, which demonstrates the capability of MethGET for use in epigenomic research. CONCLUSIONS: MethGET is a Python software that correlates DNA methylation and gene expression. Its web interface is publicly available at https://paoyang.ipmb.sinica.edu.tw/Software.html. The stand-alone version and source codes are available on GitHub at https://github.com/Jason-Teng/MethGET.

EpiMOLAS: an intuitive web-based framework for genome-wide DNA methylation analysis.

BACKGROUND: DNA methylation is a crucial epigenomic mechanism in various biological processes. Using whole-genome bisulfite sequencing (WGBS) technology, methylated cytosine sites can be revealed at the single nucleotide level. However, the WGBS data analysis process is usually complicated and challenging. RESULTS: To alleviate the associated difficulties, we integrated the WGBS data processing steps and downstream analysis into a two-phase approach. First, we set up the required tools in Galaxy and developed workflows to calculate the methylation level from raw WGBS data and generate a methylation status summary, the mtable. This computation environment is wrapped into the Docker container image DocMethyl, which allows users to rapidly deploy an executable environment without tedious software installation and library dependency problems. Next, the mtable files were uploaded to the web server EpiMOLAS_web to link with the gene annotation databases that enable rapid data retrieval and analyses. CONCLUSION: To our knowledge, the EpiMOLAS framework, consisting of DocMethyl and EpiMOLAS_web, is the first approach to include containerization technology and a web-based system for WGBS data analysis from raw data processing to downstream analysis. EpiMOLAS will help users cope with their WGBS data and also conduct reproducible analyses of publicly available data, thereby gaining insights into the mechanisms underlying complex biological phenomenon. The Galaxy Docker image DocMethyl is available at https://hub.docker.com/r/lsbnb/docmethyl/. EpiMOLAS_web is publicly accessible at http://symbiosis.iis.sinica.edu.tw/epimolas/.

Similarity between soybean and Arabidopsis seed methylomes and loss of non-CG methylation does not affect seed development.

We profiled soybean and Arabidopsis methylomes from the globular stage through dormancy and germination to understand the role of methylation in seed formation. CHH methylation increases significantly during development throughout the entire seed, targets primarily transposable elements (TEs), is maintained during endoreduplication, and drops precipitously within the germinating seedling. By contrast, no significant global changes in CG- and CHG-context methylation occur during the same developmental period. An Arabidopsis ddcc mutant lacking CHH and CHG methylation does not affect seed development, germination, or major patterns of gene expression, implying that CHH and CHG methylation does not play a significant role in seed development or in regulating seed gene activity. By contrast, over 100 TEs are transcriptionally de-repressed in ddcc seeds, suggesting that the increase in CHH-context methylation may be a failsafe mechanism to reinforce transposon silencing. Many genes encoding important classes of seed proteins, such as storage proteins, oil biosynthesis enzymes, and transcription factors, reside in genomic regions devoid of methylation at any stage of seed development. Many other genes in these classes have similar methylation patterns, whether the genes are active or repressed. Our results suggest that methylation does not play a significant role in regulating large numbers of genes important for programming seed development in both soybean and Arabidopsis. We conclude that understanding the mechanisms controlling seed development will require determining how cis-regulatory elements and their cognate transcription factors are organized in genetic regulatory networks.

Detecting and prioritizing biosynthetic gene clusters for bioactive compounds in bacteria and fungi.

Secondary metabolites (SM) produced by fungi and bacteria have long been of exceptional interest owing to their unique biomedical ramifications. The traditional discovery of new natural products that was mainly driven by bioactivity screening has now experienced a fresh new approach in the form of genome mining. Several bioinformatics tools have been continuously developed to detect potential biosynthetic gene clusters (BGCs) that are responsible for the production of SM. Although the principles underlying the computation of these tools have been discussed, the biological background is left underrated and ambiguous. In this review, we emphasize the biological hypotheses in BGC formation driven from the observations across genomes in bacteria and fungi, and provide a comprehensive list of updated algorithms/tools exclusively for BGC detection. Our review points to a direction that the biological hypotheses should be systematically incorporated into the BGC prediction and assist the prioritization of candidate BGC.

H2B ubiquitylation and the histone chaperone Asf1 cooperatively mediate the formation and maintenance of heterochromatin silencing.

Heterochromatin is a heritable form of gene repression, with critical roles in development and cell identity. Understanding how chromatin factors results in such repression is a fundamental question. Chromatin is assembled and disassembled during transcription, replication and repair by anti-silencing function 1 (Asf1), a highly conserved histone chaperone. Transcription and DNA replication are also affected by histone modifications that modify nucleosome dynamics, such as H2B ubiquitylation (H2Bub). We report here that H2Bub and Asf1 cooperatively promote transcriptional silencing at yeast telomeres and mating loci. Through real time monitoring of HML (Hidden MAT Left) locus silencing, we found that transcriptional repression was slowly initiated and never fully established in mutants lacking both Asf1 and H2Bub. These findings are consistent with impaired HML silencer-binding and spreading of repressor proteins, Sir2 and Sir3. In addition, mutants lacking H2Bub and Asf1 show defects in both nucleosome assembly and higher-order heterochromatin organization at the HML locus. Our findings reveal a novel role for H2Bub and Asf1 in epigenetic silencing at mating loci. Thus, the interplay between H2Hbub and Asf1 may fine-tune nucleosome dynamics and SIR protein recruitment, and represent an ongoing requirement for proper formation and maintenance of heterochromatin.

Polybenzyl Glutamate Biocompatible Scaffold Promotes the Efficiency of Retinal Differentiation toward Retinal Ganglion Cell Lineage from Human-Induced Pluripotent Stem Cells.

Optic neuropathy is one of the leading causes of irreversible blindness caused by retinal ganglion cell (RGC) degeneration. The development of induced pluripotent stem cell (iPSC)-based therapy opens a therapeutic window for RGC degeneration, and tissue engineering may further promote the efficiency of differentiation process of iPSCs. The present study was designed to evaluate the effects of a novel biomimetic polybenzyl glutamate (PBG) scaffold on culturing iPSC-derived RGC progenitors. The iPSC-derived neural spheres cultured on PBG scaffold increased the differentiated retinal neurons and promoted the neurite outgrowth in the RGC progenitor layer. Additionally, iPSCs cultured on PBG scaffold formed the organoid-like structures compared to that of iPSCs cultured on cover glass within the same culture period. With RNA-seq, we found that cells of the PBG group were differentiated toward retinal lineage and may be related to the glutamate signaling pathway. Further ontological analysis and the gene network analysis showed that the differentially expressed genes between cells of the PBG group and the control group were mainly associated with neuronal differentiation, neuronal maturation, and more specifically, retinal differentiation and maturation. The novel electrospinning PBG scaffold is beneficial for culturing iPSC-derived RGC progenitors as well as retinal organoids. Cells cultured on PBG scaffold differentiate effectively and shorten the process of RGC differentiation compared to that of cells cultured on coverslip. The new culture system may be helpful in future disease modeling, pharmacological screening, autologous transplantation, as well as narrowing the gap to clinical application.

BS-Seeker3: ultrafast pipeline for bisulfite sequencing.

BACKGROUND: DNA methylation is an important epigenetic modification critical in regulation and transgenerational inheritance. The methylation level can be estimated at single-nucleotide resolution by whole-genome bisulfite sequencing (BS-seq; WGBS). Current bisulfite aligners provide pipelines for processing the reads by WGBS; however, few are able to analyze the BS-seqs in a reasonable timeframe that meets the needs of the rapid expansion of epigenome sequencing in biomedical research. RESULTS: We introduce BS-Seeker3, an extensively improved and optimized implementation of BS-Seeker2 that leverages the available computational power of a standard bioinformatics lab. BS-Seeker3 adopts all alignment features of BS-Seeker2. It performs ultrafast alignments and achieves both high accuracy and high mappability, more than twice that of the other aligners that we evaluated. Moreover, BS Seeker 3 is well linked with downstream analyzer MethGo for up to 9 types of genomic and epigenomic analyses. CONCLUSIONS: BS-Seeker3 is an accurate, versatile, ultra-fast pipeline for processing bisulfite-converted reads. It also helps the user better visualize the methylation data.

Deubiquitinating Enzyme OTU5 Contributes to DNA Methylation Patterns and Is Critical for Phosphate Nutrition Signals.

Phosphate (Pi) starvation induces a suite of adaptive responses aimed at recalibrating cellular Pi homeostasis. Plants harboring a mutation in OVARIAN TUMOR DOMAIN-CONTAINING DEUBIQUITINATING ENZYME5 (OTU5) showed altered DNA methylation of root hair-related genes and altered Pi-responsive root traits. Unlike the wild type, homozygous otu5 mutants did not respond to Pi starvation by increased lateral root formation and increased root hair length but formed very short root hairs when grown on low-Pi media. Under Pi-replete conditions, otu5 plants developed more root hairs than the wild type due to attenuated primary root growth, a phenotype that resembled that of Pi-deficient plants. Growth of plants on low-Pi media altered both H3K4 and H3K27 trimethylation levels at the transcriptional start site of a subset of genes encoding key players in Pi homeostasis, which was correlated with mRNA abundance changes of these genes. Pi starvation had a minor impact on DNA methylation. Differentially methylated regions were enriched in transposable elements, suggesting that DNA methylation associated with low Pi supply is required for maintaining genome integrity. It is concluded that DNA methylation and histone methylation constitute critical, interdependent regulatory components that orchestrate the activity of a subset of Pi-responsive genes.

TEA: the epigenome platform for Arabidopsis methylome study.

BACKGROUND: Bisulfite sequencing (BS-seq) has become a standard technology to profile genome-wide DNA methylation at single-base resolution. It allows researchers to conduct genome-wise cytosine methylation analyses on issues about genomic imprinting, transcriptional regulation, cellular development and differentiation. One single data from a BS-Seq experiment is resolved into many features according to the sequence contexts, making methylome data analysis and data visualization a complex task. RESULTS: We developed a streamlined platform, TEA, for analyzing and visualizing data from whole-genome BS-Seq (WGBS) experiments conducted in the model plant Arabidopsis thaliana. To capture the essence of the genome methylation level and to meet the efficiency for running online, we introduce a straightforward method for measuring genome methylation in each sequence context by gene. The method is scripted in Java to process BS-Seq mapping results. Through a simple data uploading process, the TEA server deploys a web-based platform for deep analysis by linking data to an updated Arabidopsis annotation database and toolkits. CONCLUSIONS: TEA is an intuitive and efficient online platform for analyzing the Arabidopsis genomic DNA methylation landscape. It provides several ways to help users exploit WGBS data. TEA is freely accessible for academic users at: http://tea.iis.sinica.edu.tw .

Relationship between nucleosome positioning and DNA methylation.

Nucleosomes compact and regulate access to DNA in the nucleus, and are composed of approximately 147 bases of DNA wrapped around a histone octamer. Here we report a genome-wide nucleosome positioning analysis of Arabidopsis thaliana using massively parallel sequencing of mononucleosomes. By combining this data with profiles of DNA methylation at single base resolution, we identified 10-base periodicities in the DNA methylation status of nucleosome-bound DNA and found that nucleosomal DNA was more highly methylated than flanking DNA. These results indicate that nucleosome positioning influences DNA methylation patterning throughout the genome and that DNA methyltransferases preferentially target nucleosome-bound DNA. We also observed similar trends in human nucleosomal DNA, indicating that the relationships between nucleosomes and DNA methyltransferases are conserved. Finally, as has been observed in animals, nucleosomes were highly enriched on exons, and preferentially positioned at intron-exon and exon-intron boundaries. RNA polymerase II (Pol II) was also enriched on exons relative to introns, consistent with the hypothesis that nucleosome positioning regulates Pol II processivity. DNA methylation is also enriched on exons, consistent with the targeting of DNA methylation to nucleosomes, and suggesting a role for DNA methylation in exon definition.