LiBis: an ultrasensitive alignment augmentation for low-input bisulfite sequencing.

The cell-free DNA (cfDNA) methylation profile in liquid biopsy has been utilized to diagnose early-stage disease and estimate therapy response. However, typical clinical procedures are capable of purifying only very small amounts of cfDNA. Whole-genome bisulfite sequencing (WGBS) is the gold standard for measuring DNA methylation; however, WGBS using small amounts of fragmented DNA introduces a critical challenge for data analysis, namely a low-mapping ratio. The resulting low sequencing depth and low coverage of CpG sites genome-wide is a bottleneck for the clinical application of cfDNA-based WGBS assays. We developed LiBis (Low-input Bisulfite Sequencing), a novel method for low-input WGBS data alignment. By dynamically clipping initially unmapped reads and remapping clipped fragments, we judiciously rescued those reads and uniquely aligned them to the genome. By substantially increasing the mapping ratio by up to 88%, LiBis dramatically improved the number of informative CpGs and the precision in quantifying the methylation status of individual CpG sites. LiBis significantly improved the cost efficiency of low-input WGBS experiments by dynamically removing contamination introduced by random priming. The high sensitivity and cost effectiveness afforded by LiBis for low-input samples will allow the discovery of genetic and epigenetic features suitable for downstream analysis and biomarker identification using liquid biopsy.

Increased Pathogenicity and Virulence of Middle East Respiratory Syndrome Coronavirus Clade B In Vitro and In Vivo.

Middle East respiratory syndrome coronavirus (MERS-CoV) causes severe acute respiratory disease in humans. MERS-CoV strains from early epidemic clade A and contemporary epidemic clade B have not been phenotypically characterized to compare their abilities to infect cells and mice. We isolated the clade B MERS-CoV ChinaGD01 strain from a patient infected during the South Korean MERS outbreak in 2015 and compared the phylogenetics and pathogenicity of MERS-CoV EMC/2012 (clade A) and ChinaGD01 (clade B) in vitro and in vivo Genome alignment analysis showed that most clade-specific mutations occurred in the orf1ab gene, including mutations that were predicted to be potential glycosylation sites. Minor differences in viral growth but no significant differences in plaque size or sensitivity to beta interferon (IFN-ß) were detected between these two viruses in vitro ChinaGD01 virus infection induced more weight loss and inflammatory cytokine production in human DPP4-transduced mice. Viral titers were higher in the lungs of ChinaGD01-infected mice than with EMC/2012 infection. Decreased virus-specific CD4+ and CD8+ T cell numbers were detected in the lungs of ChinaGD01-infected mice. In conclusion, MERS-CoV evolution induced changes to reshape its pathogenicity and virulence in vitro and in vivo and to evade adaptive immune response to hinder viral clearance.IMPORTANCE MERS-CoV is an important emerging pathogen and causes severe respiratory infection in humans. MERS-CoV strains from early epidemic clade A and contemporary epidemic clade B have not been phenotypically characterized to compare their abilities to infect cells and mice. In this study, we showed that a clade B virus ChinaGD01 strain caused more severe disease in mice, with delayed viral clearance, increased inflammatory cytokines, and decreased antiviral T cell responses, than the early clade A virus EMC/2012. Given the differences in pathogenicity of different clades of MERS-CoV, periodic assessment of currently circulating MERS-CoV is needed to monitor potential severity of zoonotic disease.

Human pluripotent stem cell-derived chondroprogenitors for cartilage tissue engineering.

The cartilage of joints, such as meniscus and articular cartilage, is normally long lasting (i.e., permanent). However, once damaged, especially in large animals and humans, joint cartilage is not spontaneously repaired. Compensating the lack of repair activity by supplying cartilage-(re)forming cells, such as chondrocytes or mesenchymal stromal cells, or by transplanting a piece of normal cartilage, has been the basis of therapy for biological restoration of damaged joint cartilage. Unfortunately, current biological therapies face problems on a number of fronts. The joint cartilage is generated de novo from a specialized cell type, termed a 'joint progenitor' or 'interzone cell' during embryogenesis. Therefore, embryonic chondroprogenitors that mimic the property of joint progenitors might be the best type of cell for regenerating joint cartilage in the adult. Pluripotent stem cells (PSCs) are expected to differentiate in culture into any somatic cell type through processes that mimic embryogenesis, making human (h)PSCs a promising source of embryonic chondroprogenitors. The major research goals toward the clinical application of PSCs in joint cartilage regeneration are to (1) efficiently generate lineage-specific chondroprogenitors from hPSCs, (2) expand the chondroprogenitors to the number needed for therapy without loss of their chondrogenic activity, and (3) direct the in vivo or in vitro differentiation of the chondroprogenitors to articular or meniscal (i.e., permanent) chondrocytes rather than growth plate (i.e., transient) chondrocytes. This review is aimed at providing the current state of research toward meeting these goals. We also include our recent achievement of successful generation of "permanent-like" cartilage from long-term expandable, hPSC-derived ectomesenchymal chondroprogenitors.

Structural basis for effects of CpA modifications on C/EBPß binding of DNA.

CCAAT/enhancer binding proteins (C/EBPs) regulate gene expression in a variety of cells/tissues/organs, during a range of developmental stages, under both physiological and pathological conditions. C/EBP-related transcription factors have a consensus binding specificity of 5'-TTG-CG-CAA-3', with a central CpG/CpG and two outer CpA/TpG dinucleotides. Methylation of the CpG and CpA sites generates a DNA element with every pyrimidine having a methyl group in the 5-carbon position (thymine or 5-methylcytosine (5mC)). To understand the effects of both CpG and CpA modification on a centrally-important transcription factor, we show that C/EBPß binds the methylated 8-bp element with modestly-increased (2.4-fold) binding affinity relative to the unmodified cognate sequence, while cytosine hydroxymethylation (particularly at the CpA sites) substantially decreased binding affinity (36-fold). The structure of C/EBPß DNA binding domain in complex with methylated DNA revealed that the methyl groups of the 5mCpA/TpG make van der Waals contacts with Val285 in C/EBPß. Arg289 recognizes the central 5mCpG by forming a methyl-Arg-G triad, and its conformation is constrained by Val285 and the 5mCpG methyl group. We substituted Val285 with Ala (V285A) in an Ala-Val dipeptide, to mimic the conserved Ala-Ala in many members of the basic leucine-zipper family of transcription factors, important in gene regulation, cell proliferation and oncogenesis. The V285A variant demonstrated a 90-fold binding preference for methylated DNA (particularly 5mCpA methylation) over the unmodified sequence. The smaller side chain of Ala285 permits Arg289 to adopt two alternative conformations, to interact in a similar fashion with either the central 5mCpG or the TpG of the opposite strand. Significantly, the best-studied cis-regulatory elements in RNA polymerase II promoters and enhancers have variable sequences corresponding to the central CpG or reduced to a single G:C base pair, but retain a conserved outer CpA sequence. Our analyses suggest an important modification-dependent CpA recognition by basic leucine-zipper transcription factors.

GsmPlot: a web server to visualize epigenome data in NCBI.

BACKGROUND: Epigenetic regulation is essential in regulating gene expression across a variety of biological processes. Many high-throughput sequencing technologies have been widely used to generate epigenetic data, such as histone modification, transcription factor binding sites, DNA modifications, chromatin accessibility, and etc. A large scale of epigenetic data is stored in NCBI Gene Expression Omnibus (GEO). However, it is a great challenge to reanalyze these large scale and complex data, especially for researchers who do not specialize in bioinformatics skills or do not have access to expensive computational infrastructure. RESULTS: GsmPlot can simply accept GSM IDs to automatically download NCBI data or can accept user's private bigwig files as input to plot the concerned data on promoters, exons or any other user-defined genome locations and generate UCSC visualization tracks. By linking public data repository and private data, GsmPlot can spark data-driven ideas and hence promote the epigenetic research. CONCLUSIONS: GsmPlot web server allows convenient visualization and efficient exploration of any NCBI epigenetic data in any genomic region without need of any bioinformatics skills or special computing resources. GsmPlot is freely available at https://gsmplot.deqiangsun.org/.

Decoding the dynamic DNA methylation and hydroxymethylation landscapes in endodermal lineage intermediates during pancreatic differentiation of hESC.

Dynamic changes in DNA methylation and demethylation reprogram transcriptional outputs to instruct lineage specification during development. Here, we applied an integrative epigenomic approach to unveil DNA (hydroxy)methylation dynamics representing major endodermal lineage intermediates during pancreatic differentiation of human embryonic stem cells (hESCs). We found that 5-hydroxymethylcytosine (5hmC) marks genomic regions to be demethylated in the descendent lineage, thus reshaping the DNA methylation landscapes during pancreatic lineage progression. DNA hydroxymethylation is positively correlated with enhancer activities and chromatin accessibility, as well as the selective binding of lineage-specific pioneer transcription factors, during pancreatic differentiation. We further discovered enrichment of hydroxymethylated regions (termed '5hmC-rim') at the boundaries of large hypomethylated functional genomic regions, including super-enhancer, DNA methylation canyon and broad-H3K4me3 peaks. We speculate that '5hmC-rim' might safeguard low levels of cytosine methylation at these regions. Our comprehensive analysis highlights the importance of dynamic changes of epigenetic landscapes in driving pancreatic differentiation of hESC.

Smurf2, an E3 ubiquitin ligase, interacts with PDE4B and attenuates liver fibrosis through miR-132 mediated CTGF inhibition.

We previously reported that Smad ubiquitin regulatory factor 2 (Smurf2) activity was decreased in human fibrotic livers. Here, we overexpressed Smurf2 in livers of transgenic mice and observed inhibited collagen deposition and hepatic stellate cell activation in fibrotic model induced by carbon tetrachloride treatment or bile duct ligation. Hepatic Smurf2 overexpression also inhibited the production of connective tissue growth factor (CTGF), a central mediator of liver fibrosis. Using miRNA array and bioinformatics analyses, we identified miR-132 as a mediator of this inhibitory effect. miR-132 directly targets the 3'-untranslated region of CTGF and was transcriptionally upregulated by cAMP-PKA-CREB signaling. In addition, Smurf2 activated cAMP-PKA-CREB pathway by interacting with phosphodiesterase 4B (PDE4B) and facilitating its degradation. Thus, we have demonstrated a previously unrecognized anti-fibrotic pathway controlled by Smurf2.

Engineered Split-TET2 Enzyme for Inducible Epigenetic Remodeling.

The Ten-eleven translocation (TET) family of 5-methylcytosine (5mC) dioxygenases catalyze the conversion of 5mC into 5-hydroxymethylcytosine (5hmC) and further oxidized species to promote active DNA demethylation. Here we engineered a split-TET2 enzyme to enable temporal control of 5mC oxidation and subsequent remodeling of epigenetic states in mammalian cells. We further demonstrate the use of this chemically inducible system to dissect the correlation between DNA hydroxymethylation and chromatin accessibility in the mammalian genome. This chemical-inducible epigenome remodeling tool will find broad use in interrogating cellular systems without altering the genetic code, as well as in probing the epigenotype-phenotype relations in various biological systems.

Oncogenic targets Mmp7, S100a9, Nppb and Aldh1a3 from transcriptome profiling of FAP and Pirc adenomas are downregulated in response to tumor suppression by Clotam.

Intervention strategies in familial adenomatous polyposis (FAP) patients and other high-risk colorectal cancer (CRC) populations have highlighted a critical need for endoscopy combined with safe and effective preventive agents. We performed transcriptome profiling of colorectal adenomas from FAP patients and the polyposis in rat colon (Pirc) preclinical model, and prioritized molecular targets for prevention studies in vivo. At clinically relevant doses in the Pirc model, the drug Clotam (tolfenamic acid, TA) was highly effective at suppressing tumorigenesis both in the colon and in the small intestine, when administered alone or in combination with Sulindac. Cell proliferation in the colonic crypts was reduced significantly by TA, coincident with increased cleaved caspase-3 and decreased Survivin, ß-catenin, cyclin D1 and matrix metalloproteinase 7. From the list of differentially expressed genes prioritized by transcriptome profiling, Mmp7, S100a9, Nppb and Aldh1a3 were defined as key oncogene candidates downregulated in colon tumors after TA treatment. Monthly colonoscopies revealed the rapid onset of tumor suppression by TA in the Pirc model, and the temporal changes in Mmp7, S100a9, Nppb and Aldh1a3, highlighting their value as potential early biomarkers for prevention in the clinical setting. We conclude that TA, an "old drug" repurposed from migraine, offers an exciting new therapeutic avenue in FAP and other high-risk CRC patient populations.

BSeQC: quality control of bisulfite sequencing experiments.

MOTIVATION: Bisulfite sequencing (BS-seq) has emerged as the gold standard to study genome-wide DNA methylation at single-nucleotide resolution. Quality control (QC) is a critical step in the analysis pipeline to ensure that BS-seq data are of high quality and suitable for subsequent analysis. Although several QC tools are available for next-generation sequencing data, most of them were not designed to handle QC issues specific to BS-seq protocols. Therefore, there is a strong need for a dedicated QC tool to evaluate and remove potential technical biases in BS-seq experiments. RESULTS: We developed a package named BSeQC to comprehensively evaluate the quality of BS-seq experiments and automatically trim nucleotides with potential technical biases that may result in inaccurate methylation estimation. BSeQC takes standard SAM/BAM files as input and generates bias-free SAM/BAM files for downstream analysis. Evaluation based on real BS-seq data indicates that the use of the bias-free SAM/BAM file substantially improves the quantification of methylation level. AVAILABILITY AND IMPLEMENTATION: BSeQC is freely available at: http://code.google.com/p/bseqc/.