FinaleMe: Predicting DNA methylation by the fragmentation patterns of plasma cell-free DNA.

Analysis of DNA methylation in cell-free DNA reveals clinically relevant biomarkers but requires specialized protocols such as whole-genome bisulfite sequencing. Meanwhile, millions of cell-free DNA samples are being profiled by whole-genome sequencing. Here, we develop FinaleMe, a non-homogeneous Hidden Markov Model, to predict DNA methylation of cell-free DNA and, therefore, tissues-of-origin, directly from plasma whole-genome sequencing. We validate the performance with 80 pairs of deep and shallow-coverage whole-genome sequencing and whole-genome bisulfite sequencing data.

FinaleMe: Predicting DNA methylation by the fragmentation patterns of plasma cell-free DNA.

Analysis of DNA methylation in cell-free DNA (cfDNA) reveals clinically relevant biomarkers but requires specialized protocols and sufficient input material that limits its applicability. Millions of cfDNA samples have been profiled by genomic sequencing. To maximize the gene regulation information from the existing dataset, we developed FinaleMe, a non-homogeneous Hidden Markov Model (HMM), to predict DNA methylation of cfDNA and, therefore, tissues-of-origin directly from plasma whole-genome sequencing (WGS). We validated the performance with 80 pairs of deep and shallow-coverage WGS and whole-genome bisulfite sequencing (WGBS) data.

NOMe-HiC: joint profiling of genetic variant, DNA methylation, chromatin accessibility, and 3D genome in the same DNA molecule.

Cis-regulatory elements are coordinated to regulate the expression of their targeted genes. However, the joint measurement of cis-regulatory elements' activities and their interactions in spatial proximity is limited by the current sequencing approaches. We describe a method, NOMe-HiC, which simultaneously captures single-nucleotide polymorphisms, DNA methylation, chromatin accessibility (GpC methyltransferase footprints), and chromosome conformation changes from the same DNA molecule, together with the transcriptome, in a single assay. NOMe-HiC shows high concordance with state-of-the-art mono-omic assays across different molecular measurements and reveals coordinated chromatin accessibility at distal genomic segments in spatial proximity and novel types of long-range allele-specific chromatin accessibility.

CRAG: de novo characterization of cell-free DNA fragmentation hotspots in plasma whole-genome sequencing.

The fine-scale cell-free DNA fragmentation patterns in early-stage cancers are poorly understood. We developed a de novo approach to characterize the cell-free DNA fragmentation hotspots from plasma whole-genome sequencing. Hotspots are enriched in open chromatin regions, and, interestingly, 3'end of transposons. Hotspots showed global hypo-fragmentation in early-stage liver cancers and are associated with genes involved in the initiation of hepatocellular carcinoma and associated with cancer stem cells. The hotspots varied across multiple early-stage cancers and demonstrated high performance for the diagnosis and identification of tissue-of-origin in early-stage cancers. We further validated the performance with a small number of independent case-control-matched early-stage cancer samples.

Human fetal cerebellar cell atlas informs medulloblastoma origin and oncogenesis.

Medulloblastoma (MB) is the most common malignant childhood brain tumour1,2, yet the origin of the most aggressive subgroup-3 form remains elusive, impeding development of effective targeted treatments. Previous analyses of mouse cerebella3-5 have not fully defined the compositional heterogeneity of MBs. Here we undertook single-cell profiling of freshly isolated human fetal cerebella to establish a reference map delineating hierarchical cellular states in MBs. We identified a unique transitional cerebellar progenitor connecting neural stem cells to neuronal lineages in developing fetal cerebella. Intersectional analysis revealed that the transitional progenitors were enriched in aggressive MB subgroups, including group 3 and metastatic tumours. Single-cell multi-omics revealed underlying regulatory networks in the transitional progenitor populations, including transcriptional determinants HNRNPH1 and SOX11, which are correlated with clinical prognosis in group 3 MBs. Genomic and Hi-C profiling identified de novo long-range chromatin loops juxtaposing HNRNPH1/SOX11-targeted super-enhancers to cis-regulatory elements of MYC, an oncogenic driver for group 3 MBs. Targeting the transitional progenitor regulators inhibited MYC expression and MYC-driven group 3 MB growth. Our integrated single-cell atlases of human fetal cerebella and MBs show potential cell populations predisposed to transformation and regulatory circuitries underlying tumour cell states and oncogenesis, highlighting hitherto unrecognized transitional progenitor intermediates predictive of disease prognosis and potential therapeutic vulnerabilities.

FinaleDB: a browser and database of cell-free DNA fragmentation patterns.

SUMMARY: Circulating cell-free DNA (cfDNA) is a promising biomarker for the diagnosis and prognosis of many diseases, including cancer. The genome-wide non-random fragmentation patterns of cfDNA are associated with the nucleosomal protection, epigenetic environment and gene expression in the cell types that contributed to cfDNA. However, current progress on the development of computational methods and understanding of molecular mechanisms behind cfDNA fragmentation patterns is significantly limited by the controlled-access of cfDNA whole-genome sequencing (WGS) dataset. Here, we present FinaleDB (FragmentatIoN AnaLysis of cEll-free DNA DataBase), a comprehensive database to host thousands of uniformly processed and curated de-identified cfDNA WGS datasets across different pathological conditions. Furthermore, FinaleDB comes with a fragmentation genome browser, from which users can seamlessly integrate thousands of other omics data in different cell types to experience a comprehensive view of both gene-regulatory landscape and cfDNA fragmentation patterns. AVAILABILITY AND IMPLEMENTATION: FinaleDB service: http://finaledb.research.cchmc.org/. FinaleDB source code: https://github.com/epifluidlab/finaledb_portal, https://github.com/epifluidlab/finaledb_workflow. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Extracellular RNA in a single droplet of human serum reflects physiologic and disease states.

Extracellular RNAs (exRNAs) are present in human serum. It remains unclear to what extent these circulating exRNAs may reflect human physiologic and disease states. Here, we developed SILVER-seq (Small Input Liquid Volume Extracellular RNA Sequencing) to efficiently sequence both integral and fragmented exRNAs from a small droplet (5 µL to 7 µL) of liquid biopsy. We calibrated SILVER-seq in reference to other RNA sequencing methods based on milliliters of input serum and quantified droplet-to-droplet and donor-to-donor variations. We carried out SILVER-seq on more than 150 serum droplets from male and female donors ranging from 18 y to 48 y of age. SILVER-seq detected exRNAs from more than a quarter of the human genes, including small RNAs and fragments of mRNAs and long noncoding RNAs (lncRNAs). The detected exRNAs included those derived from genes with tissue (e.g., brain)-specific expression. The exRNA expression levels separated the male and female samples and were correlated with chronological age. Noncancer and breast cancer donors exhibited pronounced differences, whereas donors with or without cancer recurrence exhibited moderate differences in exRNA expression patterns. Even without using differentially expressed exRNAs as features, nearly all cancer and noncancer samples and a large portion of the recurrence and nonrecurrence samples could be correctly classified by exRNA expression values. These data suggest the potential of using exRNAs in a single droplet of serum for liquid biopsy-based diagnostics.

Gemini surfactant-induced DNA condensation into a beadlike structure.

Using a magnetic tweezers (MT) apparatus and an atomic force microscope (AFM), we studied the condensation of DNA induced by the cationic gemini surfactant hexyl-alpha,omega-bisdodecyldimethylammonium bromide (C12 C6 C12 Br2). Stepwise condensation events were found for forces from 0.5 to 4 pN, with a decrease in DNA extension by approximately 100 nm in each condensation event. Applied stretching forces larger than 6 pN were found to be able to decondense the condensates in a similar intermittent stepwise manner, but with approximately 60 nm of DNA released in each decondensation event. These observations are consistent with AFM images that show beadlike structures on DNA. The results lead to a model in which, during condensation, independent beadlike structures are initially formed on DNA, and as the local density of the beads increases, they compact into higher-order structures while maintaining their independence. The condensation process is different from any previously reported condensation process, but it is somewhat similar to that of chromatin assembly.