Machine learning analysis of bleeding status in venous thromboembolism patients.

Background: Anticoagulation therapy is the mainstay of therapy for patients with venous thromboembolism (VTE). However, continuing or stopping anticoagulants after the first 3 to 6 months is a difficult decision that requires ascertainment of the risk of bleeding and recurrent VTE. Despite the development of several statistical models to predict bleeding, the benefit of machine learning (ML) models has not been investigated in depth.

A deep-learning approach to predict bleeding risk over time in patients on extended anticoagulation therapy.

Background: Thus far, all the clinical models developed to predict major bleeding in patients on extended anticoagulation therapy use the baseline predictors to stratify patients into different risk groups. Therefore, these models do not account for the clinical changes and events that occur after the baseline visit, which can modify risk of bleeding. However, it is difficult to develop predictive models from the routine follow-up clinical interviews, which are irregular sequences of multivariate time series data.

Empirical Bayes single nucleotide variant-calling for next-generation sequencing data.

One of the fundamental computational problems in cancer genomics is the identification of single nucleotide variants (SNVs) from DNA sequencing data. Many statistical models and software implementations for SNV calling have been developed in the literature, yet, they still disagree widely on real datasets. Based on an empirical Bayesian approach, we introduce a local false discovery rate (LFDR) estimator for germline SNV calling.

A Markov constraint to uniquely identify elementary flux mode weights in unimolecular metabolic networks.

Elementary flux modes (EFMs) are minimal, steady state pathways characterizing a flux network. Fundamentally, all steady state fluxes in a network are decomposable into a linear combination of EFMs. While there is typically no unique set of EFM weights that reconstructs these fluxes, several optimization-based methods have been proposed to constrain the solution space by enforcing some notion of parsimony.

Quantification of Muscle Satellite Stem Cell Divisions by High-Content Analysis.

High-content screening is commonly performed on 2D cultured cells, which is high throughput but has low biological relevance. In contrast, single myofiber culture assay preserves the satellite cell niche between the basal lamina and sarcolemma and consequently has high biological relevance but is low throughput. We describe here a high-content screening method that utilizes single myofiber culture that addresses the caveats of both techniques.

Costs of Next-Generation Sequencing Assays in Non-Small Cell Lung Cancer: A Micro-Costing Study.

Next-generation sequencing (NGS) of tumor genomes has changed and improved cancer treatment over the past few decades. It can inform clinicians on the optimal therapeutic approach in many of the solid and hematologic cancers, including non-small lung cancer (NSCLC). Our study aimed to determine the costs of NGS assays for NSCLC diagnostics.

BATL: Bayesian annotations for targeted lipidomics.

Motivation: Bioinformatic tools capable of annotating, rapidly and reproducibly, large, targeted lipidomic datasets are limited. Specifically, few programs enable high-throughput peak assessment of liquid chromatography-electrospray ionization tandem mass spectrometry data acquired in either selected or multiple reaction monitoring modes.

Results: We present here Bayesian Annotations for Targeted Lipidomics, a Gaussian naïve Bayes classifier for targeted lipidomics that annotates peak identities according to eight features related to retention time, intensity, and peak shape.

WACS: improving ChIP-seq peak calling by optimally weighting controls.

Background: Chromatin immunoprecipitation followed by high throughput sequencing (ChIP-seq), initially introduced more than a decade ago, is widely used by the scientific community to detect protein/DNA binding and histone modifications across the genome. Every experiment is prone to noise and bias, and ChIP-seq experiments are no exception. To alleviate bias, the incorporation of control datasets in ChIP-seq analysis is an essential step.

Bayesian correlation is a robust gene similarity measure for single-cell RNA-seq data.

Assessing similarity is highly important for bioinformatics algorithms to determine correlations between biological information. A common problem is that similarity can appear by chance, particularly for low expressed entities. This is especially relevant in single-cell RNA-seq (scRNA-seq) data because read counts are much lower compared to bulk RNA-seq.

Absolute quantification of transcription factors in human erythropoiesis using selected reaction monitoring mass spectrometry.

Quantitative changes in transcription factor (TF) abundance regulate dynamic cellular processes, including cell fate decisions. Protein copy number provides information about the relative stoichiometry of TFs that can be used to determine how quantitative changes in TF abundance influence gene regulatory networks. In this protocol, we describe a targeted selected reaction monitoring (SRM)-based mass-spectrometry method to systematically measure the absolute protein concentration of nuclear TFs as human hematopoietic stem and progenitor cells differentiate along the erythropoietic lineage.

Accuracy and reproducibility of somatic point mutation calling in clinical-type targeted sequencing data.

Background: Treating cancer depends in part on identifying the mutations driving each patient's disease. Many clinical laboratories are adopting high-throughput sequencing for assaying patients' tumours, applying targeted panels to formalin-fixed paraffin-embedded tumour tissues to detect clinically-relevant mutations. While there have been some benchmarking and best practices studies of this scenario, much variant calling work focuses on whole-genome or whole-exome studies, with fresh or fresh-frozen tissue.

Myf6/MRF4 is a myogenic niche regulator required for the maintenance of the muscle stem cell pool.

The function and maintenance of muscle stem cells (MuSCs) are tightly regulated by signals originating from their niche environment. Skeletal myofibers are a principle component of the MuSC niche and are in direct contact with the muscle stem cells. Here, we show that Myf6 establishes a ligand/receptor interaction between muscle stem cells and their associated muscle fibers.

Absolute Quantification of Transcription Factors Reveals Principles of Gene Regulation in Erythropoiesis.

Dynamic cellular processes such as differentiation are driven by changes in the abundances of transcription factors (TFs). However, despite years of studies, our knowledge about the protein copy number of TFs in the nucleus is limited. Here, by determining the absolute abundances of 103 TFs and co-factors during the course of human erythropoiesis, we provide a dynamic and quantitative scale for TFs in the nucleus.

Establishment of macaque trophoblast stem cell lines derived from cynomolgus monkey blastocysts.

The placenta forms a maternal-fetal junction that supports many physiological functions such as the supply of nutrition and exchange of gases and wastes. Establishing an in vitro culture model of human and non-human primate trophoblast stem/progenitor cells is important for investigating the process of early placental development and trophoblast differentiation. In this study, we have established five trophoblast stem cell (TSC) lines from cynomolgus monkey blastocysts, named macTSC #1-5.

Transcriptomically Guided Mesendoderm Induction of Human Pluripotent Stem Cells Using a Systematically Defined Culture Scheme.

Human pluripotent stem cells (hPSCs) are an essential cell source in tissue engineering, studies of development, and disease modeling. Efficient, broadly amenable protocols for rapid lineage induction of hPSCs are of great interest in the stem cell biology field. We describe a simple, robust method for differentiation of hPSCs into mesendoderm in defined conditions utilizing single-cell seeding (SCS) and BMP4 and Activin A (BA) treatment.

Human gene expression variability and its dependence on methylation and aging.

Background: Phenotypic variability of human populations is partly the result of gene polymorphism and differential gene expression. As such, understanding the molecular basis for diversity requires identifying genes with both high and low population expression variance and identifying the mechanisms underlying their expression control. Key issues remain unanswered with respect to expression variability in human populations.

High-resolution genome-wide expression analysis of single myofibers using SMART-Seq.

Skeletal muscle is a heterogeneous tissue. Individual myofibers that make up muscle tissue exhibit variation in their metabolic and contractile properties. Although biochemical and histological assays are available to study myofiber heterogeneity, efficient methods to analyze the whole transcriptome of individual myofibers are lacking.

Preclinical rationale for entinostat in embryonal rhabdomyosarcoma.

Background: Rhabdomyosarcoma (RMS) is the most common soft tissue sarcoma in the pediatric cancer population. Survival among metastatic RMS patients has remained dismal yet unimproved for years. We previously identified the class I-specific histone deacetylase inhibitor, entinostat (ENT), as a pharmacological agent that transcriptionally suppresses the PAX3:FOXO1 tumor-initiating fusion gene found in alveolar rhabdomyosarcoma (aRMS), and we further investigated the mechanism by which ENT suppresses PAX3:FOXO1 oncogene and demonstrated the preclinical efficacy of ENT in RMS orthotopic allograft and patient-derived xenograft (PDX) models.

RECAP reveals the true statistical significance of ChIP-seq peak calls.

Motivation: Chromatin Immunopreciptation (ChIP)-seq is used extensively to identify sites of transcription factor binding or regions of epigenetic modifications to the genome. A key step in ChIP-seq analysis is peak calling, where genomic regions enriched for ChIP versus control reads are identified. Many programs have been designed to solve this task, but nearly all fall into the statistical trap of using the data twice-once to determine candidate enriched regions, and again to assess enrichment by classical statistical hypothesis testing.

The HDAC3-SMARCA4-miR-27a axis promotes expression of the fusion oncogene in rhabdomyosarcoma.

Rhabdomyosarcoma (RMS) is the most common soft tissue sarcoma of childhood with an unmet clinical need for decades. A single oncogenic fusion gene is associated with treatment resistance and a 40 to 45% decrease in overall survival. We previously showed that expression of this fusion oncogene in alveolar RMS (aRMS) mediates tolerance to chemotherapy and radiotherapy and that the class I-specific histone deacetylase (HDAC) inhibitor entinostat reduces PAX3:FOXO1 protein abundance.

Targeting the MTF2-MDM2 Axis Sensitizes Refractory Acute Myeloid Leukemia to Chemotherapy.

Deep sequencing has revealed that epigenetic modifiers are the most mutated genes in acute myeloid leukemia (AML). Thus, elucidating epigenetic dysregulation in AML is crucial to understand disease mechanisms. Here, we demonstrate that metal response element binding transcription factor 2/polycomblike 2 (MTF2/PCL2) plays a fundamental role in the polycomb repressive complex 2 (PRC2) and that its loss elicits an altered epigenetic state underlying refractory AML.

Cis-regulatory determinants of MyoD function.

Muscle-specific transcription factor MyoD orchestrates the myogenic gene expression program by binding to short DNA motifs called E-boxes within myogenic cis-regulatory elements (CREs). Genome-wide analyses of MyoD cistrome by chromatin immnunoprecipitation sequencing shows that MyoD-bound CREs contain multiple E-boxes of various sequences. However, how E-box numbers, sequences and their spatial arrangement within CREs collectively regulate the binding affinity and transcriptional activity of MyoD remain largely unknown.

Mtf2-PRC2 control of canonical Wnt signaling is required for definitive erythropoiesis.

Polycomb repressive complex 2 (PRC2) accessory proteins play substoichiometric, tissue-specific roles to recruit PRC2 to specific genomic loci or increase enzymatic activity, while PRC2 core proteins are required for complex stability and global levels of trimethylation of histone 3 at lysine 27 (H3K27me3). Here, we demonstrate a role for the classical PRC2 accessory protein Mtf2/Pcl2 in the hematopoietic system that is more akin to that of a core PRC2 protein. erythroid progenitors demonstrate markedly decreased core PRC2 protein levels and a global loss of H3K27me3 at promoter-proximal regions.