A DNA methylation haplotype block landscape in human tissues and preimplantation embryos reveals regulatory elements defined by comethylation patterns.

DNA methylation and associated regulatory elements play a crucial role in gene expression regulation. Previous studies have focused primarily on the distribution of mean methylation levels. Advances in whole-genome bisulfite sequencing (WGBS) have enabled the characterization of DNA methylation haplotypes (MHAPs), representing CpG sites from the same read fragment on a single chromosome, and the subsequent identification of methylation haplotype blocks (MHBs), in which adjacent CpGs on the same fragment are comethylated. Using our expert-curated WGBS data sets, we report comprehensive landscapes of MHBs in 17 representative normal somatic human tissues and during early human embryonic development. Integrative analysis reveals MHBs as a distinctive type of regulatory element characterized by comethylation patterns rather than mean methylation levels. We show the enrichment of MHBs in open chromatin regions, tissue-specific histone marks, and enhancers, including super-enhancers. Moreover, we find that MHBs tend to localize near tissue-specific genes and show an association with differential gene expression that is independent of mean methylation. Similar findings are observed in the context of human embryonic development, highlighting the dynamic nature of MHBs during early development. Collectively, our comprehensive MHB landscapes provide valuable insights into the tissue specificity and developmental dynamics of DNA methylation.

Beyond the base pairs: comparative genome-wide DNA methylation profiling across sequencing technologies.

Deoxyribonucleic acid (DNA) methylation plays a key role in gene regulation and is critical for development and human disease. Techniques such as whole-genome bisulfite sequencing (WGBS) and reduced representation bisulfite sequencing (RRBS) allow DNA methylation analysis at the genome scale, with Illumina NovaSeq 6000 and MGI Tech DNBSEQ-T7 being popular due to their efficiency and affordability. However, detailed comparative studies of their performance are not available. In this study, we constructed 60 WGBS and RRBS libraries for two platforms using different types of clinical samples and generated approximately 2.8 terabases of sequencing data. We systematically compared quality control metrics, genomic coverage, CpG methylation levels, intra- and interplatform correlations, and performance in detecting differentially methylated positions. Our results revealed that the DNBSEQ platform exhibited better raw read quality, although base quality recalibration indicated potential overestimation of base quality. The DNBSEQ platform also showed lower sequencing depth and less coverage uniformity in GC-rich regions than did the NovaSeq platform and tended to enrich methylated regions. Overall, both platforms demonstrated robust intra- and interplatform reproducibility for RRBS and WGBS, with NovaSeq performing better for WGBS, highlighting the importance of considering these factors when selecting a platform for bisulfite sequencing.

mHapBrowser: a comprehensive database for visualization and analysis of DNA methylation haplotypes.

DNA methylation acts as a vital epigenetic regulatory mechanism involved in controlling gene expression. Advances in sequencing technologies have enabled characterization of methylation patterns at single-base resolution using bisulfite sequencing approaches. However, existing methylation databases have primarily focused on mean methylation levels, overlooking phased methylation patterns. The methylation status of CpGs on individual sequencing reads represents discrete DNA methylation haplotypes (mHaps). Here, we present mHapBrowser, a comprehensive database for visualizing and analyzing mHaps. We systematically processed data of diverse tissues in human, mouse and rat from public repositories, generating mHap format files for 6366 samples. mHapBrowser enables users to visualize eight mHap metrics across the genome through an integrated WashU Epigenome Browser. It also provides an online server for comparing mHap patterns across samples. Additionally, mHap files for all samples can be downloaded to facilitate local processing using downstream analysis toolkits. The utilities of mHapBrowser were demonstrated through three case studies: (i) mHap patterns are associated with gene expression; (ii) changes in mHap patterns independent of mean methylation correlate with differential expression between lung cancer subtypes; and (iii) the mHap metric MHL outperforms mean methylation for classifying tumor and normal samples from cell-free DNA. The database is freely accessible at http://mhap.sibcb.ac.cn/.

TRAmHap: accurate prediction of transcriptional activity from DNA methylation haplotypes in bisulfite-sequencing data.

Deoxyribonucleic acid (DNA) methylation (DNAm) is an important epigenetic mechanism that plays a role in chromatin structure and transcriptional regulation. Elucidating the relationship between DNAm and gene expression is of great importance for understanding its role in transcriptional regulation. The conventional approach is to construct machine-learning-based methods to predict gene expression based on mean methylation signals in promoter regions. However, this type of strategy only explains about 25% of gene expression variation, and hence is inadequate in elucidating the relationship between DNAm and transcriptional activity. In addition, using mean methylation as input features neglects the heterogeneity of cell populations that can be reflected by DNAm haplotypes. We here developed TRAmaHap, a novel deep-learning framework that predicts gene expression by utilizing the characteristics of DNAm haplotypes in proximal promoters and distal enhancers. Using benchmark data of human and mouse normal tissues, TRAmHap shows much higher accuracy than existing machine-learning based methods, by explaining 60~80% of gene expression variation across tissue types and disease conditions. Our model demonstrated that gene expression can be accurately predicted by DNAm patterns in promoters and long-range enhancers as far as 25 kb away from transcription start site, especially in the presence of intra-gene chromatin interactions.

Discovery of novel coumarin-based KRAS-G12C inhibitors from virtual screening and Rational structural optimization.

KRAS-G12C inhibitors has been made significant progress in the treatment of KRAS-G12C mutant cancers, but their clinical application is limited due to the adaptive resistance, motivating development of novel structural inhibitors. Herein, series of coumarin derivatives as KRAS-G12C inhibitors were found through virtual screening and rational structural optimization. Especially, K45 exhibited strong antiproliferative potency on NCI-H23 and NCI-H358 cancer cells harboring KRAS-G12C with the IC50 values of 0.77 µM and 1.50 µM, which was 15 and 11 times as potent as positive drug ARS1620, respectively. Furthermore, K45 reduced the phosphorylation of KRAS downstream effectors ERK and AKT by reducing the active form of KRAS (KRAS GTP) in NCI-H23 cells. In addition, K45 induced cell apoptosis by increasing the expression of anti-apoptotic protein BAD and BAX in NCI-H23 cells. Docking studies displayed that the 3-naphthylmethoxy moiety of K45 extended into the cryptic pocket formed by the residues Gln99 and Val9, which enhanced the interaction with the KRAS-G12C protein. These results indicated that K45 was a potent KRAS-G12C inhibitor worthy of further study.

Robust Cooperative Fault-Tolerant Control for Uncertain Multi-Agent Systems Subject to Actuator Faults.

This article investigates the robust cooperative fault-tolerant control problem of multi-agent systems subject to mismatched uncertainties and actuator faults. During the design process of the intermediate variable estimator, there is no need to satisfy fault estimation matching conditions, and this overcomes a crucial constraint of traditional observers and estimators. The feedback term of the designed estimator contains the centralized estimation errors and the distributed estimation errors of the agent, and this further improves the design freedom of the proposed estimator. A novel fault-tolerant control protocol is designed based on the fault estimation information. In this work, the bounds of the fault and its derivatives are unknown, and the considered method is applicable to both directed and undirected multi-agent systems. Furthermore, the parameters of the estimator are determined through the resolution of a linear matrix inequality (LMI), which is decoupled by employing coordinate transformation and Schur decomposition. Lastly, a numerical simulation result is used to demonstrate the effectiveness of the proposed method.

mHapTk: a comprehensive toolkit for the analysis of DNA methylation haplotypes.

SUMMARY: Bisulfite sequencing remains the gold standard technique to detect DNA methylation profiles at single-nucleotide resolution. The DNA methylation status of CpG sites on the same fragment represents a discrete methylation haplotype (mHap). The mHap-level metrics were demonstrated to be promising cancer biomarkers and explain more gene expression variation than average methylation. However, most existing tools focus on average methylation and neglect mHap patterns. Here, we present mHapTk, a comprehensive python toolkit for the analysis of DNA mHap. It calculates eight mHap-level summary statistics in predefined regions or across individual CpG in a genome-wide manner. It identifies methylation haplotype blocks, in which methylations of pairwise CpGs are tightly correlated. Furthermore, mHap patterns can be visualized with the built-in functions in mHapTk or external tools such as IGV and deepTools. AVAILABILITY AND IMPLEMENTATION: https://jiantaoshi.github.io/mhaptk/index.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Progression signature underlies clonal evolution and dissemination of multiple myeloma.

Clonal evolution drives tumor progression, dissemination, and relapse in multiple myeloma (MM), with most patients dying of relapsed disease. This multistage process requires tumor cells to enter the circulation, extravasate, and colonize distant bone marrow (BM) sites. Here, we developed a fluorescent or DNA-barcode clone-tracking system on MM PrEDiCT (progression through evolution and dissemination of clonal tumor cells) xenograft mouse model to study clonal behavior within the BM microenvironment. We showed that only the few clones that successfully adapt to the BM microenvironment can enter the circulation and colonize distant BM sites. RNA sequencing of primary and distant-site MM tumor cells revealed a progression signature sequentially activated along human MM progression and significantly associated with overall survival when evaluated against patient data sets. A total of 28 genes were then computationally predicted to be master regulators (MRs) of MM progression. HMGA1 and PA2G4 were validated in vivo using CRISPR-Cas9 in the PrEDiCT model and were shown to be significantly depleted in distant BM sites, indicating their role in MM progression and dissemination. Loss of HMGA1 and PA2G4 also compromised the proliferation, migration, and adhesion abilities of MM cells in vitro. Overall, our model successfully recapitulates key characteristics of human MM disease progression and identified potential new therapeutic targets for MM.

Design, synthesis, and biological evaluation of quinoxalinone derivatives as potent BRD4 inhibitors.

The bromodomain-containing protein 4 (BRD4) has gained growing interest as an effective drug target for the treatment of hepatocellular carcinoma (HCC). Herein, we designed and synthesized a series of quinoxalinone derivatives as BRD4 inhibitors via scaffold hopping. The representative compound X9 showed potent BRD4 inhibitory activity (with IC50 = 82.3 nM), and preferable antiproliferative activity against HepG2 cells (with IC50 = 1.13 ± 0.07 µM), as well as less toxicity against GES-1 cells (with IC50 = 57.24 ± 5.46 µM). Furthermore, compound X9 dose-dependently inhibited colony formation and blocked the migration of HepG2 cells by down-regulating the expression of Snail and MMP-9 while up-regulating the E-cadherin and Occludin. Besides, compound X9 efficiently down-regulated the expression of c-Myc in HepG2 cells, induced apoptosis, and arrested at G0/G1 phase. In total, quinoxalinone was a potential core as BRD4 inhibitor and compound X9 might be effective for liver cancer therapy.

1H-Indazoles derivatives targeting PI3K/AKT/mTOR pathway: Synthesis, anti-tumor effect and molecular mechanism.

The PI3K/AKT/mTOR signaling pathway is one of the most common abnormal activation pathways in tumor cells, and has associated with multiple functions such as tumor cell growth, proliferation, migration, invasion, and tumor angiogenesis. Here, a series of 3-amino-1H-indazole derivatives were synthesized, and their antiproliferative activities against HT-29, MCF-7, A-549, HepG2 and HGC-27 cells were evaluated. Among them, W24 exhibited the broad-spectrum antiproliferative activity against four cancer cells with IC50 values of 0.43-3.88 µM. Mechanism studies revealed that W24 inhibited proliferation by affecting the DNA synthesis, induced G2/M cell cycle arrest and apoptosis by regulating Cyclin B1, BAD and Bcl-xL, meanwhile induced the change of intracellular ROS and mitochondrial membrane potential in HGC-27 cells. Moreover, W24 inhibited the migration and invasion of HGC-27 cells by decreasing EMT pathway related proteins and reducing the mRNA expression levels of Snail, Slug and HIF-1α. Furthermore, W24 displayed low tissue toxicity profile and good pharmacokinetic properties in vivo. Therefore, 3-amino-1H-indazole derivatives might serve as a new scaffold for the development of PI3K/AKT/mTOR inhibitor and anti-gastric cancer reagent.

Epigenetic restriction of extraembryonic lineages mirrors the somatic transition to cancer.

In mammals, the canonical somatic DNA methylation landscape is established upon specification of the embryo proper and subsequently disrupted within many cancer types. However, the underlying mechanisms that direct this genome-scale transformation remain elusive, with no clear model for its systematic acquisition or potential developmental utility. Here, we analysed global remethylation from the mouse preimplantation embryo into the early epiblast and extraembryonic ectoderm. We show that these two states acquire highly divergent genomic distributions with substantial disruption of bimodal, CpG density-dependent methylation in the placental progenitor. The extraembryonic epigenome includes specific de novo methylation at hundreds of embryonically protected CpG island promoters, particularly those that are associated with key developmental regulators and are orthologously methylated across most human cancer types. Our data suggest that the evolutionary innovation of extraembryonic tissues may have required co-option of DNA methylation-based suppression as an alternative to regulation by Polycomb-group proteins, which coordinate embryonic germ-layer formation in response to extraembryonic cues. Moreover, we establish that this decision is made deterministically, downstream of promiscuously used-and frequently oncogenic-signalling pathways, via a novel combination of epigenetic cofactors. Methylation of developmental gene promoters during tumorigenesis may therefore reflect the misappropriation of an innate trajectory and the spontaneous reacquisition of a latent, developmentally encoded epigenetic landscape.

The DNA methylation haplotype (mHap) format and mHapTools.

SUMMARY: Bisulfite sequencing (BS-seq) is currently the gold standard for measuring genome-wide DNA methylation profiles at single-nucleotide resolution. Most analyses focus on mean CpG methylation and ignore methylation states on the same DNA fragments [DNA methylation haplotypes (mHaps)]. Here, we propose mHap, a simple DNA mHap format for storing DNA BS-seq data. This format reduces the size of a BAM file by 40- to 140-fold while retaining complete read-level CpG methylation information. It is also compatible with the Tabix tool for fast and random access. We implemented a command-line tool, mHapTools, for converting BAM/SAM files from existing platforms to mHap files as well as post-processing DNA methylation data in mHap format. With this tool, we processed all publicly available human reduced representation bisulfite sequencing data and provided these data as a comprehensive mHap database. AVAILABILITY AND IMPLEMENTATION: https://jiantaoshi.github.io/mHap/index.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Design, synthesis and biological evaluation of novel pyrazinone derivatives as PI3K/HDAC dual inhibitors.

PI3Ks and HDACs play essential roles in the occurrence and progression of leukemia. Herein, a series of novel pyrazin-2(1H)-one derivatives were rationally designed and synthesized as novel dual PI3K and HDAC inhibitors based on scaffold replacement and heterozygous strategies. Most of the target compounds showed potent inhibitory potency to PI3Kα and HDAC6. Especially, compound 9q displayed PI3Kα and HDAC6 inhibitory with IC50 values of 372 nM and 4.5 nM, and anti-proliferative activity against MV4-11 cells with IC50 value of 0.093 ± 0.012 µM. Further mechanistic studies revealed that 9q induced apoptosis, arrested the cell cycle in the G2/M phase, promoted the acetylation of α-tubulin, and blocked the PI3K/AKT/mTOR signal way in MV4-11 cells. All the results demonstrated that 9q was a promising lead candidate for further development of novel PI3K/HDAC dual inhibitors for leukemia treatment.

Design, synthesis and biological evaluation of coumarin derivatives as potential BRD4 inhibitors.

The bromodomain and extra-terminal (BET) bromodomains, particularly BRD4, have been identified as promising therapeutic targets in the treatment of many human disorders such as cancer. Coumarin is a highly privileged moiety for the development of novel anticancer drugs which has been identified in clinical trials for the treatment of various cancers. Herein, we modified BRD4i ABBV-075 with a coumarin ring and synthesized a novel series of coumarin derivatives as BRD4 inhibitors. Among them, the representative compound 27d showed excellent BRD4 inhibitory activities with an IC50 value of 99 nM in the TR-FRET assay. Compared with ABBV-075, compound 27d displayed a favorable cell proliferation inhibitory activity in solid tumors, such as MCF-7, HGC-27 and HepG-2. Further mechanism investigation illustrated that 27d-treatment resulted in G0/G1 phase arrest and promoted apoptosis of MCF-7 cells. Compound 27d also blocked colony formation in a concentration-dependent manner in McF-7 cell lines. As the downstream-protein of BRD4, the expression of c-Myc was decreased in a dose-dependent manner after the treatment of compound 27d. Moreover, compound 27d also exhibited good in vivo and in vitro metabolic stability. All the findings meaningfully make it as a promising lead compound for further drug development.

DIFFpop: a stochastic computational approach to simulate differentiation hierarchies with single cell barcoding.

SUMMARY: DIFFpop is an R package designed to simulate cellular differentiation hierarchies using either exponentially-expanding or fixed population sizes. The software includes functionalities to simulate clonal evolution due to the emergence of driver mutations under the infinite-allele assumption as well as options for simulation and analysis of single cell barcoding and labeling data. The software uses the Gillespie Stochastic Simulation Algorithm and a modification of expanding or fixed-size stochastic process models expanded to a large number of cell types and scenarios. AVAILABILITY AND IMPLEMENTATION: DIFFpop is available as an R-package along with vignettes on Github (https://github.com/ferlicjl/diffpop). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Design, synthesis and biological evaluation of novel 4-phenoxypyridine based 3-oxo-3,4-dihydroquinoxaline-2-carboxamide derivatives as potential c-Met kinase inhibitors.

Blocking c-Met kinase activity by small-molecule inhibitors has been identified as a promising approach for the treatment of cancers. Herein, we described the design, synthesis, and biological evaluation of a series of 4-phenoxypyridine-based 3-oxo-3,4-dihydroquinoxaline derivatives as c-Met kinase inhibitors. Inhibitory activitives against c-Met kinase evaluation indicated that most of compounds showed excellent c-Met kinase activity in vitro, and IC50 values of ten compounds (23a, 23e, 23f, 23l, 23r, 23s, 23v, 23w, 23x and 23y) were less than 10.00 nM. Notably, three of them (23v, 23w and 23y) showed remarkable potency with IC50 values of 2.31 nM, 1.91 nM and 2.44 nM, respectively, and thus they were more potent than positive control drug foretinib (c-Met, IC50 = 2.53 nM). Cytotoxic evaluation indicated the most promising compound 23w showed remarkable cytotoxicity against A549, H460 and HT-29 cell lines with IC50 values of 1.57 µM, 0.94 µM and 0.65 µM, respectively. Furthermore, the acridine orange/ethidium bromide (AO/EB) staining, cell apoptosis assays by flow cytometry, wound-healing assays and transwell migration assays on HT-29 and/or A549 cells of 23w were performed. Especially compound 23w, which displayed potent antitumor, apoptosis induction and antimetastatic activity, could be used as a promising lead for further development. Meanwhile, their preliminary structure-activity relationships (SARs) were also discussed.

Prognostic role of circulating exosomal miRNAs in multiple myeloma.

Exosomes, secreted by several cell types, including cancer cells, can be isolated from the peripheral blood and have been shown to be powerful markers of disease progression in cancer. In this study, we examined the prognostic significance of circulating exosomal microRNAs (miRNAs) in multiple myeloma (MM). A cohort of 156 patients with newly diagnosed MM, uniformly treated and followed, was studied. Circulating exosomal miRNAs were isolated and used to perform a small RNA sequencing analysis on 10 samples and a quantitative reverse transcription polymerase chain reaction (qRT-PCR) array on 156 samples. We studied the relationship between miRNA levels and patient outcomes, including progression-free survival (PFS) and overall survival (OS). We identified miRNAs as the most predominant small RNAs present in exosomes isolated from the serum of patients with MM and healthy controls by small RNA sequencing of circulating exosomes. We then analyzed exosomes isolated from serum samples of 156 patients using a qRT-PCR array for 22 miRNAs. Two of these miRNAs, let-7b and miR-18a, were significantly associated with both PFS and OS in the univariate analysis and were still statistically significant after adjusting for the International Staging System and adverse cytogenetics in the multivariate analysis. Our findings support the use of circulating exosomal miRNAs to improve the identification of patients with newly diagnosed MM with poor outcomes. The results require further validation in other independent prospective MM cohorts.

Exome sequencing reveals recurrent germ line variants in patients with familial Waldenström macroglobulinemia.

Familial aggregation of Waldenström macroglobulinemia (WM) cases, and the clustering of B-cell lymphoproliferative disorders among first-degree relatives of WM patients, has been reported. Nevertheless, the possible contribution of inherited susceptibility to familial WM remains unrevealed. We performed whole exome sequencing on germ line DNA obtained from 4 family members in which coinheritance for WM was documented in 3 of them, and screened additional independent 246 cases by using gene-specific mutation sequencing. Among the shared germ line variants, LAPTM5(c403t) and HCLS1(g496a) were the most recurrent, being present in 3/3 affected members of the index family, detected in 8% of the unrelated familial cases, and present in 0.5% of the nonfamilial cases and in <0.05 of a control population. LAPTM5 and HCLS1 appeared as relevant WM candidate genes that characterized familial WM individuals and were also functionally relevant to the tumor clone. These findings highlight potentially novel contributors for the genetic predisposition to familial WM and indicate that LAPTM5(c403t) and HCLS1(g496a) may represent predisposition alleles in patients with familial WM.

Endogenous hydrogen sulphide attenuates NLRP3 inflammasome-mediated neuroinflammation by suppressing the P2X7 receptor after intracerebral haemorrhage in rats.

BACKGROUND: Emerging studies have demonstrated the important physiological and pathophysiological roles of hydrogen sulphide (H2S) as a gasotransmitter for NOD-like receptor family pyrin domain-containing 3 (NLRP3) inflammasome-associated neuroinflammation in the central nervous system. However, the effects of H2S on neuroinflammation after intracerebral haemorrhage (ICH), especially on the NLRP3 inflammasome, remain unknown. METHODS: We employed a Sprague-Dawley rat of collagenase-induced ICH in the present study. The time course of H2S content and the spatial expression of cystathionine-ß-synthase (CBS) after ICH, the effects of endogenous and exogenous H2S after ICH, the effects of endogenous and exogenous H2S on NLRP3 inflammasome activation under P2X7 receptor (P2X7R) overexpression after ICH, and the involvement of the P2X7R in the mechanism by which microglia-derived H2S prevented NLRP3 inflammasome activation were investigated. RESULTS: We found ICH induced significant downregulation of endogenous H2S production in the brain, which may be the result of decreasing in CBS, the predominant cerebral H2S-generating enzyme. Administration of S-adenosyl-L-methionine (SAM), a CBS-specific agonist, or sodium hydrosulfide (NaHS), a classical exogenous H2S donor, not only restored brain and plasma H2S content but also attenuated brain oedema, microglial accumulation and neurological deficits at 1 day post-ICH by inhibiting the P2X7R/NLRP3 inflammasome cascade. Endogenous H2S production, which was derived mainly by microglia and above treatments, was verified by adenovirus-overexpressed P2X7R and in vitro primary microglia studies. CONCLUSIONS: These results indicated endogenous H2S synthesis was impaired after ICH, which plays a pivotal role in the P2X7R/NLRP3 inflammasome-associated neuroinflammatory response in the pathogenesis of secondary brain injury. Maintaining appropriate H2S concentrations in the central nervous system may represent a potential therapeutic strategy for managing post-ICH secondary brain injury and associated neurological deficits.

edgeRun: an R package for sensitive, functionally relevant differential expression discovery using an unconditional exact test.

Next-generation sequencing platforms for measuring digital expression such as RNA-Seq are displacing traditional microarray-based methods in biological experiments. The detection of differentially expressed genes between groups of biological conditions has led to the development of numerous bioinformatics tools, but so far, few exploit the expanded dynamic range afforded by the new technologies. We present edgeRun, an R package that implements an unconditional exact test that is a more powerful version of the exact test in edgeR. This increase in power is especially pronounced for experiments with as few as two replicates per condition, for genes with low total expression and with large biological coefficient of variation. In comparison with a panel of other tools, edgeRun consistently captures functionally similar differentially expressed genes.