Whole-genome long-read TAPS deciphers DNA methylation patterns at base resolution using PacBio SMRT sequencing technology.

Long-read sequencing provides valuable information on difficult-to-map genomic regions, which can complement short-read sequencing to improve genome assembly, yet limited methods are available to accurately detect DNA methylation over long distances at a whole-genome scale. By combining our recently developed TET-assisted pyridine borane sequencing (TAPS) method, which enables direct detection of 5-methylcytosine and 5-hydroxymethylcytosine, with PacBio single-molecule real-time sequencing, we present here whole-genome long-read TAPS (wglrTAPS). To evaluate the performance of wglrTAPS, we applied it to mouse embryonic stem cells as a proof of concept, and an N50 read length of 3.5 kb is achieved. By sequencing wglrTAPS to 8.2× depth, we discovered a significant proportion of CpG sites that were not covered in previous 27.5× short-read TAPS. Our results demonstrate that wglrTAPS facilitates methylation profiling on problematic genomic regions with repetitive elements or structural variations, and also in an allelic manner, all of which are extremely difficult for short-read sequencing methods to resolve. This method therefore enhances applications of third-generation sequencing technologies for DNA epigenetics.

Modular Oxidation of Cytosine Modifications and Their Application in Direct and Quantitative Sequencing of 5-Hydroxymethylcytosine.

Selective, efficient, and controllable oxidation of cytosine modifications is valuable for epigenetic analyses, yet only limited progress has been made. Here, we present two modular chemical oxidation reactions: conversion of 5-hydroxymethylcytosine (5hmC) into 5-formylcytosine (5fC) using 4-acetamido-2,2,6,6-tetramethylpiperidine-1-oxoammonium tetrafluoroborate (ACT+BF4-) and further transformation of 5fC into 5-carboxycytosine (5caC) through Pinnick oxidation. Both reactions are mild and efficient on double-stranded DNA. We integrated these two oxidations with borane reduction to develop chemical-assisted pyridine borane sequencing plus (CAPS+), for direct and quantitative mapping of 5hmC. Compared with CAPS, CAPS+ improved the conversion rate and false-positive rate. We applied CAPS+ to mouse embryonic stem cells, human normal brain, and glioblastoma DNA samples and demonstrated its superior sensitivity in analyzing the hydroxymethylome.

Endonuclease enrichment TAPS for cost-effective genome-wide base-resolution DNA methylation detection.

Whole genome base-resolution methylome sequencing allows for the most comprehensive analysis of DNA methylation, however, the considerable sequencing cost often limits its applications. While reduced representation sequencing can be an affordable alternative, over 80% of CpGs in the genome are not covered. Building on our recently developed TET-assisted pyridine borane sequencing (TAPS) method, we here described endonuclease enrichment TAPS (eeTAPS), which utilizes dihydrouracil (DHU)-cleaving endonuclease digestion of TAPS-converted DNA to enrich methylated CpG sites (mCpGs). eeTAPS can accurately detect 87% of mCpGs in the mouse genome with a sequencing depth equivalent to 4× whole genome sequencing. In comparison, reduced representation TAPS (rrTAPS) detected less than 4% of mCpGs with 2.5× sequencing depth. Our results demonstrate eeTAPS to be a new strategy for cost-effective genome-wide methylation analysis at single-CpG resolution that can fill the gap between whole-genome and reduced representation sequencing.

Plant Regulomics: a data-driven interface for retrieving upstream regulators from plant multi-omics data.

High-throughput technology has become a powerful approach for routine plant research. Interpreting the biological significance of high-throughput data has largely focused on the functional characterization of a large gene list or genomic loci that involves the following two aspects: the functions of the genes or loci and how they are regulated as a whole, i.e. searching for the upstream regulators. Traditional platforms for functional annotation largely help resolving the first issue. Addressing the second issue is essential for a global understanding of the regulatory mechanism, but is more challenging, and requires additional high-throughput experimental evidence and a unified statistical framework for data-mining. The rapid accumulation of 'omics data provides a large amount of experimental data. We here present Plant Regulomics, an interface that integrates 19 925 transcriptomic and epigenomic data sets and diverse sources of functional evidence (58 112 terms and 695 414 protein-protein interactions) from six plant species along with the orthologous genes from 56 whole-genome sequenced plant species. All pair-wise transcriptomic comparisons with biological significance within the same study were performed, and all epigenomic data were processed to genomic loci targeted by various factors. These data were well organized to gene modules and loci lists, which were further implemented into the same statistical framework. For any input gene list or genomic loci, Plant Regulomics retrieves the upstream factors, treatments, and experimental/environmental conditions regulating the input from the integrated 'omics data. Additionally, multiple tools and an interactive visualization are available through a user-friendly web interface. Plant Regulomics is available at http://bioinfo.sibs.ac.cn/plant-regulomics.

Polycomb repressive complex 2 attenuates ABA-induced senescence in Arabidopsis.

The phytohormone abscisic acid (ABA)-induced leaf senescence facilitates nutrient reuse and potentially contributes to enhancing plant stress tolerance. However, excessive senescence causes serious reductions in crop yield, and the mechanism by which senescence is finely tuned at different levels is still insufficiently understood. Here, we found that the double mutant of core enzymes of the polycomb repressive complex 2 (PRC2) is hypersensitive to ABA in Arabidopsis thaliana. To elucidate the interplay between ABA and PRC2 at the genome level, we extensively profiled the transcriptomic and epigenomic changes triggered by ABA. We observed that H3K27me3 preferentially targets ABA-induced senescence-associated genes (SAGs). In the double, but not single, mutant of PRC2 enzymes, these SAGs were derepressed and could be more highly induced by ABA compared with the wild-type, suggesting a redundant role for the PRC2 enzymes in negatively regulating ABA-induced senescence. Contrary to the rapid transcriptomic changes triggered by ABA, the reduction of H3K27me3 at these SAGs falls far behind the induction of their expression, indicating that PRC2-mediated H3K27me3 contributed to long-term damping of ABA-induced senescence to prevent an oversensitive response. The findings of this study may serve as a paradigm for a global understanding of the interplay between the rapid effects of a phytohormone such as ABA and the long-term effects of the epigenetic machinery in regulating plant senescence processes and environmental responses.

Antitumor activity of Z-endoxifen in aromatase inhibitor-sensitive and aromatase inhibitor-resistant estrogen receptor-positive breast cancer.

BACKGROUND: The tamoxifen metabolite, Z-endoxifen, demonstrated promising antitumor activity in endocrine-resistant estrogen receptor-positive (ER+) breast cancer. We compared the antitumor activity of Z-endoxifen with tamoxifen and letrozole in the letrozole-sensitive MCF7 aromatase expressing model (MCF7AC1), as well as with tamoxifen, fulvestrant, exemestane, and exemestane plus everolimus in a letrozole-resistant MCF7 model (MCF7LR). METHODS: MCF7AC1 tumor-bearing mice were randomized to control (no drug), letrozole (10 µg/day), tamoxifen (500 µg/day), or Z-endoxifen (25 and 75 mg/kg). Treatment in the letrozole arm was continued until resistance developed. MCF7LR tumor-bearing mice were then randomized to Z-endoxifen (50 mg/kg) or tamoxifen for 4 weeks and tumors harvested for microarray and immunohistochemistry analysis. The antitumor activity of Z-endoxifen in the MCF7LR tumors was further compared in a second in vivo study with exemestane, exemestane plus everolimus, and fulvestrant. RESULTS: In the MCF7AC1 tumors, both Z-endoxifen doses were significantly superior to control and tamoxifen in reducing tumor volumes at 4 weeks. Additionally, the 75 mg/kg Z-endoxifen dose was additionally superior to letrozole. Prolonged letrozole exposure resulted in resistance at 25 weeks. In MCF7LR tumor-bearing mice, Z-endoxifen significantly reduced tumor volumes compared to tamoxifen, letrozole, and exemestane, with no significant differences compared to exemestane plus everolimus and fulvestrant. Additionally, compared to tamoxifen, Z-endoxifen markedly inhibited ERα target genes, Ki67 and Akt expression in vivo. CONCLUSION: In endocrine-sensitive and letrozole-resistant breast tumors, Z-endoxifen results in robust antitumor and antiestrogenic activity compared to tamoxifen and aromatase inhibitor monotherapy. These data support the ongoing development of Z-endoxifen.

Critical roles of DNA demethylation in the activation of ripening-induced genes and inhibition of ripening-repressed genes in tomato fruit.

DNA methylation is a conserved epigenetic mark important for genome integrity, development, and environmental responses in plants and mammals. Active DNA demethylation in plants is initiated by a family of 5-mC DNA glycosylases/lyases (i.e., DNA demethylases). Recent reports suggested a role of active DNA demethylation in fruit ripening in tomato. In this study, we generated loss-of-function mutant alleles of a tomato gene, SlDML2, which is a close homolog of the Arabidopsis DNA demethylase gene ROS1 In the fruits of the tomato mutants, increased DNA methylation was found in thousands of genes. These genes included not only hundreds of ripening-induced genes but also many ripening-repressed genes. Our results show that SlDML2 is critical for tomato fruit ripening and suggest that active DNA demethylation is required for both the activation of ripening-induced genes and the inhibition of ripening-repressed genes.

Correction: Arabidopsis Flower and Embryo Developmental Genes are Repressed in Seedlings by Different Combinations of Polycomb Group Proteins in Association with Distinct Sets of Cis-regulatory Elements.

[This corrects the article DOI: 10.1371/journal.pgen.1005771.].

The heterochronic gene Oryza sativa LIKE HETEROCHROMATIN PROTEIN 1 modulates miR156b/c/i/e levels.

The juvenile-to-adult transition in plants involves changes in vegetative growth and plant architecture; the timing of this transition has important implications for agriculture. The microRNA miR156 regulates this transition and shoot maturation in plants. In Arabidopsis thaliana, deposition of histone H3 trimethylation on lysine 27 (H3K27me3, a repressive mark) at the MIR156A/C loci is regulated by Polycomb Repressive Complex 1 (PRC1) or PRC2, depending on the developmental stage. The levels of miR156 progressively decline during shoot maturation. The amount of H3K27me3 at MIR156A/C loci affects miR156 levels; however, whether this epigenetic regulation is conserved remains unclear. Here, we found that in rice (Oryza sativa), the putative PRC1 subunit LIKE HETEROCHROMATIN PROTEIN 1 (OsLHP1), with the miR156-SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) module, affects developmental phase transitions. Loss of OsLHP1 function results in ectopic expression of MIR156B/C/I/E, phenocopy of miR156 overexpression, and reduced H3k27me3 levels at MIR156B/C/I/E. This indicates that OsLHP1 has functionally diverged from Arabidopsis LHP1. Genetic and transcriptome analyses of wild-type, miR156b/c-overexpression, and Oslhp1-2 mutant plants suggest that OsLHP1 acts upstream of miR156 and SPL during the juvenile-to-adult transition. Therefore, modifying the OsLHP1-miR156-SPL pathway may enable alteration of the vegetative period and plant architecture.

Arabidopsis Flower and Embryo Developmental Genes are Repressed in Seedlings by Different Combinations of Polycomb Group Proteins in Association with Distinct Sets of Cis-regulatory Elements.

Polycomb repressive complexes (PRCs) play crucial roles in transcriptional repression and developmental regulation in both plants and animals. In plants, depletion of different members of PRCs causes both overlapping and unique phenotypic defects. However, the underlying molecular mechanism determining the target specificity and functional diversity is not sufficiently characterized. Here, we quantitatively compared changes of tri-methylation at H3K27 in Arabidopsis mutants deprived of various key PRC components. We show that CURLY LEAF (CLF), a major catalytic subunit of PRC2, coordinates with different members of PRC1 in suppression of distinct plant developmental programs. We found that expression of flower development genes is repressed in seedlings preferentially via non-redundant role of CLF, which specifically associated with LIKE HETEROCHROMATIN PROTEIN1 (LHP1). In contrast, expression of embryo development genes is repressed by PRC1-catalytic core subunits AtBMI1 and AtRING1 in common with PRC2-catalytic enzymes CLF or SWINGER (SWN). This context-dependent role of CLF corresponds well with the change in H3K27me3 profiles, and is remarkably associated with differential co-occupancy of binding motifs of transcription factors (TFs), including MADS box and ABA-related factors. We propose that different combinations of PRC members distinctively regulate different developmental programs, and their target specificity is modulated by specific TFs.

Auxin-Independent NAC Pathway Acts in Response to Explant-Specific Wounding and Promotes Root Tip Emergence during de Novo Root Organogenesis in Arabidopsis.

Plants have powerful regenerative abilities that allow them to recover from damage and survive in nature. De novo organogenesis is one type of plant regeneration in which adventitious roots and shoots are produced from wounded and detached organs. By studying de novo root organogenesis using leaf explants of Arabidopsis (Arabidopsis thaliana), we previously suggested that wounding is the first event that provides signals to trigger the whole regenerative process. However, our knowledge of the role of wounding in regeneration remains limited. In this study, we show that wounding not only triggers the auxin-mediated fate transition of regeneration-competent cells, but also induces the NAC pathway for root tip emergence. The NAC1 transcription factor gene was specifically expressed in response to wounding in the leaf explant, but not in the wounded leaf residue of the source plant. Inhibition of the NAC1 pathway severely affected the emergence of adventitious root tips. However, the NAC1 pathway functioned independently of auxin-mediated cell fate transition and regulates expression of CEP genes, which encode proteins that might have a role in degradation of extensin proteins in the cell wall. Overall, our results suggest that wounding has multiple roles in de novo root organogenesis and that NAC1 acts as one downstream branch in regulating the cellular environment for organ emergence.

Inhibition of p38 MAPK attenuates renal atrophy and fibrosis in a murine renal artery stenosis model.

Renal artery stenosis (RAS) is an important cause of chronic renal dysfunction. Recent studies have underscored a critical role for CCL2 (MCP-1)-mediated inflammation in the progression of chronic renal damage in RAS and other chronic renal diseases. In vitro studies have implicated p38 MAPK as a critical intermediate for the production of CCL2. However, a potential role of p38 signaling in the development and progression of chronic renal disease in RAS has not been previously defined. We sought to test the hypothesis that inhibition of p38 MAPK ameliorates chronic renal injury in mice with RAS. We established a murine RAS model by placing a cuff on the right renal artery and treated mice with the p38 inhibitor SB203580 or vehicle for 2 wk. In mice treated with vehicle, the cuffed kidney developed interstitial fibrosis, tubular atrophy, and interstitial inflammation. In mice treated with SB203580, the RAS-induced renal atrophy was reduced (70% vs. 39%, P < 0.05). SB203580 also reduced interstitial inflammation and extracellular matrix deposition but had no effect on the development of hypertension. SB203580 partially blocked the induction of CCL2, CCL7 (MCP-3), CC chemokine receptor 2 (CCR2), and collagen 4 mRNA expression in the cuffed kidneys. In vitro, blockade of p38 hindered both TNF-α and TGF-ß-induced CCL2 upregulation. Based on these observations, we conclude that p38 MAPK plays a critical role in the induction of CCL2/CCL7/CCR2 system and the development of interstitial inflammation in RAS.

Genetic deficiency of Smad3 protects the kidneys from atrophy and interstitial fibrosis in 2K1C hypertension.

Although the two-kidney, one-clip (2K1C) model is widely used as a model of human renovascular hypertension, mechanisms leading to the development of fibrosis and atrophy in the cuffed kidney and compensatory hyperplasia in the contralateral kidney have not been defined. Based on the well-established role of the transforming growth factor (TGF)-ß signaling pathway in renal fibrosis, we tested the hypothesis that abrogation of TGF-ß/Smad3 signaling would prevent fibrosis in the cuffed kidney. Renal artery stenosis (RAS) was established in mice with a targeted disruption of exon 2 of the Smad3 gene (Smad3 KO) and wild-type (WT) controls by placement of a polytetrafluoroethylene cuff on the right renal artery. Serial pulse-wave Doppler ultrasound assessments verified that blood flow through the cuffed renal artery was decreased to a similar extent in Smad3 KO and WT mice. Two weeks after surgery, systolic blood pressure and plasma renin activity were significantly elevated in both the Smad3 KO and WT mice. The cuffed kidney of WT mice developed renal atrophy (50% reduction in weight after 6 wk, P < 0.0001), which was associated with the development of interstitial fibrosis, tubular atrophy, and interstitial inflammation. Remarkably, despite a similar reduction of renal blood flow, the cuffed kidney of the Smad3 KO mice showed minimal atrophy (9% reduction in weight, P = not significant), with no significant histopathological alterations (interstitial fibrosis, tubular atrophy, and interstitial inflammation). We conclude that abrogation of TGF-ß/Smad3 signaling confers protection against the development of fibrosis and atrophy in RAS.

n-3 Fatty acids block TNF-α-stimulated MCP-1 expression in rat mesangial cells.

Monocyte chemoattractant protein 1 (MCP-1) is a CC cytokine that fundamentally contributes to the pathogenesis of inflammatory renal disease. MCP-1 is highly expressed in cytokine-stimulated mesangial cells in vitro and following glomerular injury in vivo. Interventions to limit MCP-1 expression are commonly effective in assorted experimental models. Fish oil, an abundant source of n-3 fatty acids, has anti-inflammatory properties, the basis of which remains incompletely defined. We examined potential mechanisms whereby fish oil reduces MCP-1 expression and thereby suppresses inflammatory responses to tissue injury. Cultured mesangial cells were treated with TNF-α in the presence of the n-3 fatty acids docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA); equimolar concentrations of the n-6 fatty acids LA and OA served as controls. MCP-1 mRNA expression was assessed by Northern blotting, and transcriptional activity of the MCP-1 promoter was assessed by transient transfection. The involvement of the ERK and NF-κB pathways was evaluated through transfection analysis and the use of the MEK inhibitor U0126. DHA and EPA decreased TNF-α-stimulated MCP-1 mRNA expression by decreasing transcription of the MCP-1 gene. DHA and EPA decreased p-ERK expression and nuclear translocation of NF-κB, both of which are necessary for TNF-α-stimulated MCP-1 expression. Both NF-κB and AP-1 sites were involved in transcriptional regulation of the MCP-1 gene by DHA and EPA. We conclude that DHA and EPA inhibit TNF-α-stimulated transcription of the MCP-1 gene through interaction of signaling pathways involving ERK and NF-κB. We speculate that such effects may contribute to the salutary effect of fish oil in renal and vascular disease.

Subtraction-free and bisulfite-free specific sequencing of 5-methylcytosine and its oxidized derivatives at base resolution.

Although various methods have been developed for sequencing cytosine modifications, it is still challenging for specific and quantitative sequencing of individual modification at base-resolution. For example, to obtain both true 5-methylcytosine (5mC) and true 5-hydroxymethylcytosine (5hmC) information, the two major epigenetic modifications, it usually requires subtraction of two methods, which increases noise and requires high sequencing depth. Recently, we developed TET-assisted pyridine borane sequencing (TAPS) for bisulfite-free direct sequencing of 5mC and 5hmC. Here we demonstrate that two sister methods, TAPSß and chemical-assisted pyridine borane sequencing (CAPS), can be effectively used for subtraction-free and specific whole-genome sequencing of 5mC and 5hmC, respectively. We also demonstrate pyridine borane sequencing (PS) for whole-genome profiling of 5-formylcytosine and 5-carboxylcytosine, the further oxidized derivatives of 5mC and 5hmC. This work completes the set of versatile borane reduction chemistry-based methods as a comprehensive toolkit for direct and quantitative sequencing of all four cytosine epigenetic modifications.

Homology-mediated inter-chromosomal interactions in hexaploid wheat lead to specific subgenome territories following polyploidization and introgression.

BACKGROUND: Polyploidization and introgression are major events driving plant genome evolution and influencing crop breeding. However, the mechanisms underlying the higher-order chromatin organization of subgenomes and alien chromosomes are largely unknown. RESULTS: We probe the three-dimensional chromatin architecture of Aikang 58 (AK58), a widely cultivated allohexaploid wheat variety in China carrying the 1RS/1BL translocation chromosome. The regions involved in inter-chromosomal interactions, both within and between subgenomes, have highly similar sequences. Subgenome-specific territories tend to be connected by subgenome-dominant homologous transposable elements (TEs). The alien 1RS chromosomal arm, which was introgressed from rye and differs from its wheat counterpart, has relatively few inter-chromosome interactions with wheat chromosomes. An analysis of local chromatin structures reveals topologically associating domain (TAD)-like regions covering 52% of the AK58 genome, the boundaries of which are enriched with active genes, zinc-finger factor-binding motifs, CHH methylation, and 24-nt small RNAs. The chromatin loops are mostly localized around TAD boundaries, and the number of gene loops is positively associated with gene activity. CONCLUSIONS: The present study reveals the impact of the genetic sequence context on the higher-order chromatin structure and subgenome stability in hexaploid wheat. Specifically, we characterized the sequence homology-mediated inter-chromosome interactions and the non-canonical role of subgenome-biased TEs. Our findings may have profound implications for future investigations of the interplay between genetic sequences and higher-order structures and their consequences on polyploid genome evolution and introgression-based breeding of crop plants.

Cell-free DNA TAPS provides multimodal information for early cancer detection.

Multimodal, genome-wide characterization of epigenetic and genetic information in circulating cell-free DNA (cfDNA) could enable more sensitive early cancer detection, but it is technologically challenging. Recently, we developed TET-assisted pyridine borane sequencing (TAPS), which is a mild, bisulfite-free method for base-resolution direct DNA methylation sequencing. Here, we optimized TAPS for cfDNA (cfTAPS) to provide high-quality and high-depth whole-genome cell-free methylomes. We applied cfTAPS to 85 cfDNA samples from patients with hepatocellular carcinoma (HCC) or pancreatic ductal adenocarcinoma (PDAC) and noncancer controls. From only 10 ng of cfDNA (1 to 3 ml of plasma), we generated the most comprehensive cfDNA methylome to date. We demonstrated that cfTAPS provides multimodal information about cfDNA characteristics, including DNA methylation, tissue of origin, and DNA fragmentation. Integrated analysis of these epigenetic and genetic features enables accurate identification of early HCC and PDAC.

Accurate targeted long-read DNA methylation and hydroxymethylation sequencing with TAPS.

We present long-read Tet-assisted pyridine borane sequencing (lrTAPS) for targeted base-resolution sequencing of DNA methylation and hydroxymethylation in regions up to 10 kb from nanogram-level input. Compatible with both Oxford Nanopore and PacBio Single-Molecule Real-Time (SMRT) sequencing, lrTAPS detects methylation with accuracy comparable to short-read Illumina sequencing but with long-range epigenetic phasing. We applied lrTAPS to sequence difficult-to-map regions in mouse embryonic stem cells and to identify distinct methylation events in the integrated hepatitis B virus genome.

Temporal analysis of signaling pathways activated in a murine model of two-kidney, one-clip hypertension.

Unilateral renal artery stenosis (RAS) leads to atrophy of the stenotic kidney and compensatory enlargement of the contralateral kidney. Although the two-kidney, one-clip (2K1C) model has been extensively used to model human RAS, the cellular responses in the stenotic and contralateral kidneys, particularly in the murine model, have received relatively little attention. We studied mice 2, 5, and 11 wk after unilateral RAS. These mice became hypertensive within 1 wk. The contralateral kidney increased in size within 2 wk after surgery. This enlargement was associated with a transient increase in expression of phospho-extracellular signal-regulated kinase (p-ERK), the proliferation markers proliferating cell nuclear antigen and Ki-67, the cell cycle inhibitors p21 and p27, and transforming growth factor-beta, with return to baseline levels by 11 wk. The size of the stenotic kidney was unchanged at 2 wk but progressively decreased between 5 and 11 wk. Unlike the contralateral kidney, which showed minimal histopathological alterations, the stenotic kidney developed progressive interstitial fibrosis, tubular atrophy, and interstitial inflammation. Surprisingly, the stenotic kidney showed a proliferative response, which involved largely tubular epithelial cells. The atrophic kidney had little evidence of apoptosis, despite persistent upregulation of p53; expression of cell cycle regulatory proteins in the stenotic kidney was persistently increased through 11 wk. These studies indicate that in the 2K1C model, the stenotic kidney and contralateral, enlarged kidney exhibit a distinct temporal expression of proteins involved in cell growth, cell survival, apoptosis, inflammation, and fibrosis. Notably, an unexpected proliferative response occurs in the stenotic kidney that undergoes atrophy.