SPlinted Ligation Adapter Tagging (SPLAT), a novel library preparation method for whole genome bisulphite sequencing

Effects of DNA methylation and its application in inflammatory bowel disease (Review)

Inflammatory bowel disease (IBD) is marked by persistent inflammation, and its development and progression are linked to environmental, genetic, immune system and gut microbial factors. DNA methylation (DNAm), as one of the protein modifications, is a crucial epigenetic process used by cells to control gene transcription. DNAm is one of the most common areas that has drawn increasing attention recently, with studies revealing that the interleukin (IL)‑23/IL‑12, wingless‑related integration site, IL‑6‑associated signal transducer and activator of transcription 3, suppressor of cytokine signaling 3 and apoptosis signaling pathways are involved in DNAm and in the pathogenesis of IBD. It has emerged that DNAm‑associated genes are involved in perpetuating the persistent inflammation that characterizes a number of diseases, including IBD, providing a novel therapeutic strategy for exploring their treatment. The present review discusses DNAm‑associated genes in the pathogenesis of IBD and summarizes their application as possible diagnostic, prognostic and therapeutic biomarkers in IBD. This may provide a reference for the particular form of IBD and its related methylation genes, aiding in clinical decision‑making and encouraging therapeutic alternatives.

Uncovering epigenetic and transcriptional regulation of growth in Douglas-fir: identification of differential methylation regions in mega-sized introns

Tree growth performance can be partly explained by genetics, while a large proportion of growth variation is thought to be controlled by environmental factors. However, to what extent DNA methylation, a stable epigenetic modification, contributes to phenotypic plasticity in the growth performance of long-lived trees remains unclear. In this study, a comparative analysis of targeted DNA genotyping, DNA methylation and mRNAseq profiling for needles of 44-year-old Douglas-fir trees (Pseudotsuga menziesii (Mirb.) Franco) having contrasting growth characteristics was performed. In total, we identified 195 differentially expressed genes (DEGs) and 115 differentially methylated loci (DML) that are associated with genes involved in fitness-related processes such as growth, stress management, plant development and energy resources. Interestingly, all four intronic DML were identified in mega-sized (between 100 and 180 kbp in length) and highly expressed genes, suggesting specialized regulation mechanisms of these long intron genes in gymnosperms. DNA repetitive sequences mainly comprising long-terminal repeats of retroelements are involved in growth-associated DNA methylation regulation (both hyper- and hypomethylation) of 99 DML (86.1% of total DML). Furthermore, nearly 14% of the DML was not tagged by single nucleotide polymorphisms, suggesting a unique contribution of the epigenetic variation in tree growth.

Comprehensive characterization of small noncoding RNA profiles in hypoxia-induced pulmonary hypertension (HPH) rat tissues

Hypoxia-induced pulmonary hypertension (HPH) is a fatal cardiovascular disease characterized by an elevation in pulmonary artery pressure, resulting in right ventricular dysfunction and eventual heart failure. Exploring the pathogenesis of HPH is crucial, and small noncoding RNAs (sncRNAs) are gaining recognition as potential regulators of cellular responses to hypoxia. In this study, we conducted a comprehensive analysis of sncRNA profiles in eight tissues of male HPH rats using high-throughput sequencing. Our study unveiled several sncRNAs, with the brain, kidney, and spleen exhibiting the highest abundance of microRNA (miRNA), tRNA-derived small RNA (tDR), and small nucleolar RNA (snoRNA), respectively. Moreover, we identified numerous tissue-specific and hypoxia-responsive sncRNAs, particularly miRNAs and tDRs. Interestingly, we observed arm switching in miRNAs under hypoxic conditions and a significant increase in the abundance of 5' tRNA-halves among the total tDRs during hypoxia. Overall, our study provides a comprehensive characterization of the sncRNA profiles in HPH rats.

Sequencing of N6-methyl-deoxyadenosine at single-base resolution across the mammalian genome

Although DNA N6-methyl-deoxyadenosine (6mA) is abundant in bacteria and protists, its presence and function in mammalian genomes have been less clear. We present Direct-Read 6mA sequencing (DR-6mA-seq), an antibody-independent method, to measure 6mA at base resolution. DR-6mA-seq employs a unique mutation-based strategy to reveal 6mA sites as misincorporation signatures without any chemical or enzymatic modulation of 6mA. We validated DR-6mA-seq through the successful mapping of the well-characterized G(6mA)TC motif in the E. coli DNA. As expected, when applying DR-6mA-seq to mammalian systems, we found that genomic DNA (gDNA) 6mA abundance is generally low in most mammalian tissues and cells; however, we did observe distinct gDNA 6mA sites in mouse testis and glioblastoma cells. DR-6mA-seq provides an enabling tool to detect 6mA at single-base resolution for a comprehensive understanding of DNA 6mA in eukaryotes.

Methylomics and cancer: the current state of methylation profiling and marker development for clinical care

Epigenetic modifications have long been recognized as an essential level in transcriptional regulation linking behavior and environmental conditions or stimuli with biological processes and disease development. Among them, methylation is the most abundant of these reversible epigenetic marks, predominantly occurring on DNA, RNA, and histones. Methylation modification is intimately involved in regulating gene transcription and cell differentiation, while aberrant methylation status has been linked with cancer development in several malignancies. Early detection and precise restoration of dysregulated methylation form the basis for several epigenetics-based therapeutic strategies. In this review, we summarize the current basic understanding of the regulation and mechanisms responsible for methylation modification and cover several cutting-edge research techniques for detecting methylation across the genome and transcriptome. We then explore recent advances in clinical diagnostic applications of methylation markers of various cancers and address the current state and future prospects of methylation modifications in therapies for different diseases, especially comparing pharmacological methylase/demethylase inhibitors with the CRISPRoff/on methylation editing systems. This review thus provides a resource for understanding the emerging role of epigenetic methylation in cancer, the use of methylation-based biomarkers in cancer detection, and novel methylation-targeted drugs.

Genomic, transcriptomic and epigenomic sequencing data of the B-cell leukemia cell line REH

OBJECTIVES

The aim of this data paper is to describe a collection of 33 genomic, transcriptomic and epigenomic sequencing datasets of the B-cell acute lymphoblastic leukemia (ALL) cell line REH. REH is one of the most frequently used cell lines for functional studies of pediatric ALL, and these data provide a multi-faceted characterization of its molecular features. The datasets described herein, generated with short- and long-read sequencing technologies, can both provide insights into the complex aberrant karyotype of REH, and be used as reference datasets for sequencing data quality assessment or for methods development.

DATA DESCRIPTION

This paper describes 33 datasets corresponding to 867 gigabases of raw sequencing data generated from the REH cell line. These datasets include five different approaches for whole genome sequencing (WGS) on four sequencing platforms, two RNA sequencing (RNA-seq) techniques on two different sequencing platforms, DNA methylation sequencing, and single-cell ATAC-sequencing.

The genomic basis and environmental correlates of local adaptation in the Atlantic horse mackerel (Trachurus trachurus)

Understanding how populations adapt to their environment is increasingly important to prevent biodiversity loss due to overexploitation and climate change. Here we studied the population structure and genetic basis of local adaptation of Atlantic horse mackerel, a commercially and ecologically important marine fish that has one of the widest distributions in the eastern Atlantic. We analyzed whole-genome sequencing and environmental data of samples collected from the North Sea to North Africa and the western Mediterranean Sea. Our genomic approach indicated low population structure with a major split between the Mediterranean Sea and the Atlantic Ocean and between locations north and south of mid-Portugal. Populations from the North Sea are the most genetically distinct in the Atlantic. We discovered that most population structure patterns are driven by a few highly differentiated putatively adaptive loci. Seven loci discriminate the North Sea, two the Mediterranean Sea, and a large putative inversion (9.9 Mb) on chromosome 21 underlines the north-south divide and distinguishes North Africa. A genome-environment association analysis indicates that mean seawater temperature and temperature range, or factors correlated to them, are likely the main environmental drivers of local adaptation. Our genomic data broadly support the current stock divisions, but highlight areas of potential mixing, which require further investigation. Moreover, we demonstrate that as few as 17 highly informative SNPs can genetically discriminate the North Sea and North African samples from neighboring populations. Our study highlights the importance of both, life history and climate-related selective pressures in shaping population structure patterns in marine fish. It also supports that chromosomal rearrangements play a key role in local adaptation with gene flow. This study provides the basis for more accurate delineation of the horse mackerel stocks and paves the way for improving stock assessments.

The performance of whole genome bisulfite sequencing on DNBSEQ-Tx platform examined by different library preparation strategies

Background

Whole-genome bisulfite sequencing (WGBS) technology can provide comprehensive DNA methylation at a single-base resolution on a genome-wide scale, and is considered to be the gold standard for the detection of 5-methylcytosine (5 mC). However, the International Human Epigenome Consortium propose a full DNA methylome should have at least 30 fold redundant coverage of the reference genome from a single biological replicate. Therefore, it remains cost prohibitive for large-scale studies. To find a solution, the DNBSEQ-Tx sequencing was developed that can generate up to 6 Tb data in a single run for projects involving large-scale sequencing.

Results

In this study, we provided two WGBS library construction methods DNB_PREBSseq and DNB_SPLATseq optimized for the DNBSEQ-Tx sequencer, and demonstrated the performance of these two methods on the DNBSEQ-Tx platform, using the DNA extracted from four different cell lines. We also compared the sequencing data from these two WGBS library construction methods with HeLa cell line data from ENCODE sequenced on Illumina HiSeq X Ten and WGBS data of two other cell lines sequenced on HiSeq2500. Various quality control (QC) analyses such as the base quality scores, methylation-bias (m-bias), and conversion efficiency indicated that the data sequenced on the DNBSEQ-Tx platform met the WGBS-required quality controls. Meanwhile, our data closely resembled the coverage shown by the data generated by the Illumina platform.

Conclusions

Our study showed that with our optimized methods, DNBSEQ-Tx could generate high-quality WGBS data with relatively good stability for large-scale WGBS sequencing applications. Thus, we conclude that DNBSEQ-Tx can be used for a wide range of WGBS research.

Environmental stress during larval development induces DNA methylation shifts in the migratory painted lady butterfly (Vanessa cardui)

Seasonal environmental fluctuations provide formidable challenges for living organisms, especially small ectotherms such as butterflies. A common strategy to cope with harsh environments is to enter diapause, but some species avoid unsuitable conditions by migrating. Despite a growing understanding of migration in the life cycles of some butterfly species, it remains unknown how individuals register and store environmental cues to determine whether and where to migrate. Here, we explored how competition and host plant availability during larval development affect patterns of DNA methylation in the migratory painted lady (Vanessa cardui) butterfly. We identify a set of potentially functional methylome shifts associated with differences in the environment, indicating that DNA methylation is involved in the response to different conditions during larval development. By analysing the transcriptome for the same samples used for methylation profiling, we also uncovered a non-monotonic relationship between gene body methylation and gene expression. Our results provide a starting point for understanding the interplay between DNA methylation and gene expression in butterflies in general and how differences in environmental conditions during development can trigger unique epigenetic marks that might be important for behavioural decisions in the adult stage.

Post-bisulfite Adaptor Tagging Based on an ssDNA Ligation Technique (tPBAT)

Post-bisulfite adaptor tagging (PBAT) is a concept that enables the preparation of an efficient sequencing library from bisulfite-treated DNA, and it also means the protocol implemented the concept. Although the previous PBAT or rPBAT was sensitive enough for single-cell methylome analysis, the protocol had several drawbacks owing to the repeated random priming reactions. To resolve these problems, we developed a unique single-strand DNA ligation technique, termed TACS ligation, and established a new protocol called tPBAT. With tPBAT, the data quality improved, with a longer insert and higher mapping rate than that obtained with rPBAT. In addition, paired-end sequencing and indexing were supported by the default. In this chapter, the tPBAT protocol is introduced, and a thorough description of its application to small samples is provided.

Environmentally induced DNA methylation is inherited across generations in an aquatic keystone species

Transgenerational inheritance of environmentally induced epigenetic marks can have significant impacts on eco-evolutionary dynamics, but the phenomenon remains controversial in ecological model systems. We used whole-genome bisulfite sequencing of individual water fleas (Daphnia magna) to assess whether environmentally induced DNA methylation is transgenerationally inherited. Genetically identical females were exposed to one of three natural stressors, or a de-methylating drug, and their offspring were propagated clonally for four generations under control conditions. We identified between 70 and 225 differentially methylated CpG positions (DMPs) in F1 individuals whose mothers were exposed to a natural stressor. Roughly half of these environmentally induced DMPs persisted until generation F4. In contrast, treatment with the drug demonstrated that pervasive hypomethylation upon exposure is reset almost completely after one generation. These results suggest that environmentally induced DNA methylation is non-random and stably inherited across generations in Daphnia, making epigenetic inheritance a putative factor in the eco-evolutionary dynamics of freshwater communities.

•

Naturally induced DNA-methylation persists until generation F4 in Daphnia

•

Drug-induced de-methylation is reset after one generation

•

Methylation is enriched in exons suggesting a gene regulatory function

•

Epigenetic inheritance may influence eco-evolutionary dynamics

scSPLAT, a scalable plate-based protocol for single cell WGBS library preparation

DNA methylation is a central epigenetic mark that has diverse roles in gene regulation, development, and maintenance of genome integrity. 5 methyl cytosine (5mC) can be interrogated at base resolution in single cells by using bisulfite sequencing (scWGBS). Several different scWGBS strategies have been described in recent years to study DNA methylation in single cells. However, there remain limitations with respect to cost-efficiency and yield. Herein, we present a new development in the field of scWGBS library preparation; single cell Splinted Ligation Adapter Tagging (scSPLAT). scSPLAT employs a pooling strategy to facilitate sample preparation at a higher scale and throughput than previously possible. We demonstrate the accuracy and robustness of the method by generating data from 225 single K562 cells and from 309 single liver nuclei and compare scSPLAT against other scWGBS methods.

Decoding ceRNA regulatory network in the pulmonary artery of hypoxia-induced pulmonary hypertension (HPH) rat model

Background

Hypoxia-induced pulmonary hypertension (HPH) is a lethal cardiovascular disease with the characteristic of severe remodeling of pulmonary vascular. Although a large number of dysregulated mRNAs, lncRNAs, circRNAs, and miRNAs related to HPH have been identified from extensive studies, the competitive endogenous RNA (ceRNA) regulatory network in the pulmonary artery that responds to hypoxia remains largely unknown.

Results

Transcriptomic profiles in the pulmonary arteries of HPH rats were characterized through high-throughput RNA sequencing in this study. Through relatively strict screening, a set of differentially expressed RNAs (DERNAs) including 19 DEmRNAs, 8 DElncRNAs, 19 DEcircRNAs, and 23 DEmiRNAs were identified between HPH and normal rats. The DEmRNAs were further found to be involved in cell adhesion, axon guidance, PPAR signaling pathway, and calcium signaling pathway, suggesting their crucial role in HPH. Moreover, a hypoxia-induced ceRNA regulatory network in the pulmonary arteries of HPH rats was constructed according to the ceRNA hypothesis. More specifically, the ceRNA network was composed of 10 miRNAs as hub nodes, which might be sponged by 6 circRNAs and 7 lncRNAs, and directed the expression of 18 downstream target genes that might play important role in the progression of HPH. The expression patterns of selected DERNAs in the ceRNA network were then validated to be consistent with sequencing results in another three independent batches of HPH and normal control rats. The diagnostic effectiveness of several hub mRNAs in ceRNA network was further evaluated through investigating their expression profiles in patients with pulmonary artery hypertension (PAH) recorded in the Gene Expression Omnibus (GEO) dataset GSE117261. Dysregulated POSTN, LTBP2, SPP1, and LSAMP were observed in both the pulmonary arteries of HPH rats and lung tissues of PAH patients.

Conclusions

A ceRNA regulatory network in the pulmonary arteries of HPH rats was constructed, 10 hub miRNAs and their corresponding interacting lncRNAs, circRNAs, and mRNAs were identified. The expression patterns of selected DERNAs were further validated to be consistent with the sequencing result. POSTN, LTBP2, SPP1, and LSAMP were suggested to be potential diagnostic biomarkers and therapeutic targets for PAH.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-022-00762-1.

Analysis and Performance Assessment of the Whole Genome Bisulfite Sequencing Data Workflow: Currently Available Tools and a Practical Guide to Advance DNA Methylation Studies

DNA methylation is associated with transcriptional repression, genomic imprinting, stem cell differentiation, embryonic development, and inflammation. Aberrant DNA methylation can indicate disease states, including cancer and neurological disorders. Therefore, the prevalence and location of 5-methylcytosine in the human genome is a topic of interest. Whole-genome bisulfite sequencing (WGBS) is a high-throughput method for analyzing DNA methylation. This technique involves library preparation, alignment, and quality control. Advancements in epigenetic technology have led to an increase in DNA methylation studies. This review compares the detailed experimental methodology of WGBS using accessible and up-to-date analysis tools. Practical codes for WGBS data processing are included as a general guide to assist progress in DNA methylation studies through a comprehensive case study.

Low guanine content and biased nucleotide distribution in vertebrate mtDNA can cause overestimation of non-CpG methylation

Mitochondrial DNA (mtDNA) methylation in vertebrates has been hotly debated for over 40 years. Most contrasting results have been reported following bisulfite sequencing (BS-seq) analyses. We addressed whether BS-seq experimental and analysis conditions influenced the estimation of the levels of methylation in specific mtDNA sequences. We found false positive non-CpG methylation in the CHH context (fpCHH) using unmethylated Sus scrofa mtDNA BS-seq data. fpCHH methylation was detected on the top/plus strand of mtDNA within low guanine content regions. These top/plus strand sequences of fpCHH regions would become extremely AT-rich sequences after BS-conversion, whilst bottom/minus strand sequences remained almost unchanged. These unique sequences caused BS-seq aligners to falsely assign the origin of each strand in fpCHH regions, resulting in false methylation calls. fpCHH methylation detection was enhanced by short sequence reads, short library inserts, skewed top/bottom read ratios and non-directional read mapping modes. We confirmed no detectable CHH methylation in fpCHH regions by BS-amplicon sequencing. The fpCHH peaks were located in the D-loop, ATP6, ND2, ND4L, ND5 and ND6 regions and identified in our S. scrofa ovary and oocyte data and human BS-seq data sets. We conclude that non-CpG methylation could potentially be overestimated in specific sequence regions by BS-seq analysis.

The SEQC2 epigenomics quality control (EpiQC) study

BACKGROUND

Cytosine modifications in DNA such as 5-methylcytosine (5mC) underlie a broad range of developmental processes, maintain cellular lineage specification, and can define or stratify types of cancer and other diseases. However, the wide variety of approaches available to interrogate these modifications has created a need for harmonized materials, methods, and rigorous benchmarking to improve genome-wide methylome sequencing applications in clinical and basic research. Here, we present a multi-platform assessment and cross-validated resource for epigenetics research from the FDA's Epigenomics Quality Control Group.

RESULTS

Each sample is processed in multiple replicates by three whole-genome bisulfite sequencing (WGBS) protocols (TruSeq DNA methylation, Accel-NGS MethylSeq, and SPLAT), oxidative bisulfite sequencing (TrueMethyl), enzymatic deamination method (EMSeq), targeted methylation sequencing (Illumina Methyl Capture EPIC), single-molecule long-read nanopore sequencing from Oxford Nanopore Technologies, and 850k Illumina methylation arrays. After rigorous quality assessment and comparison to Illumina EPIC methylation microarrays and testing on a range of algorithms (Bismark, BitmapperBS, bwa-meth, and BitMapperBS), we find overall high concordance between assays, but also differences in efficiency of read mapping, CpG capture, coverage, and platform performance, and variable performance across 26 microarray normalization algorithms.

CONCLUSIONS

The data provided herein can guide the use of these DNA reference materials in epigenomics research, as well as provide best practices for experimental design in future studies. By leveraging seven human cell lines that are designated as publicly available reference materials, these data can be used as a baseline to advance epigenomics research.

Enzymatic methyl sequencing detects DNA methylation at single-base resolution from picograms of DNA

Bisulfite sequencing detects 5mC and 5hmC at single-base resolution. However, bisulfite treatment damages DNA, which results in fragmentation, DNA loss, and biased sequencing data. To overcome these problems, enzymatic methyl-seq (EM-seq) was developed. This method detects 5mC and 5hmC using two sets of enzymatic reactions. In the first reaction, TET2 and T4-BGT convert 5mC and 5hmC into products that cannot be deaminated by APOBEC3A. In the second reaction, APOBEC3A deaminates unmodified cytosines by converting them to uracils. Therefore, these three enzymes enable the identification of 5mC and 5hmC. EM-seq libraries were compared with bisulfite-converted DNA, and each library type was ligated to Illumina adaptors before conversion. Libraries were made using NA12878 genomic DNA, cell-free DNA, and FFPE DNA over a range of DNA inputs. The 5mC and 5hmC detected in EM-seq libraries were similar to those of bisulfite libraries. However, libraries made using EM-seq outperformed bisulfite-converted libraries in all specific measures examined (coverage, duplication, sensitivity, etc.). EM-seq libraries displayed even GC distribution, better correlations across DNA inputs, increased numbers of CpGs within genomic features, and accuracy of cytosine methylation calls. EM-seq was effective using as little as 100 pg of DNA, and these libraries maintained the described advantages over bisulfite sequencing. EM-seq library construction, using challenging samples and lower DNA inputs, opens new avenues for research and clinical applications.

Ultrasensitive detection of circulating tumour DNA via deep methylation sequencing aided by machine learning

The low abundance of circulating tumour DNA (ctDNA) in plasma samples makes the analysis of ctDNA biomarkers for the detection or monitoring of early-stage cancers challenging. Here we show that deep methylation sequencing aided by a machine-learning classifier of methylation patterns enables the detection of tumour-derived signals at dilution factors as low as 1 in 10,000. For a total of 308 patients with surgery-resectable lung cancer and 261 age- and sex-matched non-cancer control individuals recruited from two hospitals, the assay detected 52-81% of the patients at disease stages IA to III with a specificity of 96% (95% confidence interval (CI) 93-98%). In a subgroup of 115 individuals, the assay identified, at 100% specificity (95% CI 91-100%), nearly twice as many patients with cancer as those identified by ultradeep mutation sequencing analysis. The low amounts of ctDNA permitted by machine-learning-aided deep methylation sequencing could provide advantages in cancer screening and the assessment of treatment efficacy.

Small RNA-Sequencing: Approaches and Considerations for miRNA Analysis

MicroRNAs (miRNAs) are a class of small RNA molecules that have an important regulatory role in multiple physiological and pathological processes. Their disease-specific profiles and presence in biofluids are properties that enable miRNAs to be employed as non-invasive biomarkers. In the past decades, several methods have been developed for miRNA analysis, including small RNA sequencing (RNA-seq). Small RNA-seq enables genome-wide profiling and analysis of known, as well as novel, miRNA variants. Moreover, its high sensitivity allows for profiling of low input samples such as liquid biopsies, which have now found applications in diagnostics and prognostics. Still, due to technical bias and the limited ability to capture the true miRNA representation, its potential remains unfulfilled. The introduction of many new small RNA-seq approaches that tried to minimize this bias, has led to the existence of the many small RNA-seq protocols seen today. Here, we review all current approaches to cDNA library construction used during the small RNA-seq workflow, with particular focus on their implementation in commercially available protocols. We provide an overview of each protocol and discuss their applicability. We also review recent benchmarking studies comparing each protocol's performance and summarize the major conclusions that can be gathered from their usage. The result documents variable performance of the protocols and highlights their different applications in miRNA research. Taken together, our review provides a comprehensive overview of all the current small RNA-seq approaches, summarizes their strengths and weaknesses, and provides guidelines for their applications in miRNA research.

Adaptor Template Oligo-Mediated Sequencing (ATOM-Seq) is a new ultra-sensitive UMI-based NGS library preparation technology for use with cfDNA and cfRNA

Liquid biopsy testing utilising Next Generation Sequencing (NGS) is rapidly moving towards clinical adoption for personalised oncology. However, before NGS can fulfil its potential any novel testing approach must identify ways of reducing errors, allowing separation of true low-frequency mutations from procedural artefacts, and be designed to improve upon current technologies. Popular NGS technologies typically utilise two DNA capture approaches; PCR and ligation, which have known limitations and seem to have reached a development plateau with only small, stepwise improvements being made. To maximise the ultimate utility of liquid biopsy testing we have developed a highly versatile approach to NGS: Adaptor Template Oligo Mediated Sequencing (ATOM-Seq). ATOM-Seq's strengths and versatility avoid the major limitations of both PCR- and ligation-based approaches. This technology is ligation free, simple, efficient, flexible, and streamlined, and it offers novel advantages that make it perfectly suited for use on highly challenging clinical material. Using reference and clinical materials, we demonstrate detection of known SNVs down to allele frequencies of 0.1% using as little as 20–25 ng of cfDNA, as well as the ability to detect fusions from RNA. We illustrate ATOM-Seq’s suitability for clinical testing by showing high concordance rates between paired cfDNA and FFPE clinical samples.

A low-bias and sensitive small RNA library preparation method using randomized splint ligation

Small RNAs are important regulators of gene expression and are involved in human development and disease. Next generation sequencing (NGS) allows for scalable, genome-wide studies of small RNA; however, current methods are challenged by low sensitivity and high bias, limiting their ability to capture an accurate representation of the cellular small RNA population. Several studies have shown that this bias primarily arises during the ligation of single-strand adapters during library preparation, and that this ligation bias is magnified by 2′-O-methyl modifications (2′OMe) on the 3′ terminal nucleotide. In this study, we developed a novel library preparation process using randomized splint ligation with a cleavable adapter, a design which resolves previous challenges associated with this ligation strategy. We show that a randomized splint ligation based workflow can reduce bias and increase the sensitivity of small RNA sequencing for a wide variety of small RNAs, including microRNA (miRNA) and tRNA fragments as well as 2′OMe modified RNA, including Piwi-interacting RNA and plant miRNA. Finally, we demonstrate that this workflow detects more differentially expressed miRNA between tumorous and matched normal tissues. Overall, this library preparation process allows for highly accurate small RNA sequencing and will enable studies of 2′OMe modified RNA with new levels of detail.

Manual and automated preparation of single-stranded DNA libraries for the sequencing of DNA from ancient biological remains and other sources of highly degraded DNA

It has been shown that highly fragmented DNA is most efficiently converted into DNA libraries for sequencing if both strands of the DNA fragments are processed independently. We present an updated protocol for library preparation from single-stranded DNA, which is based on the splinted ligation of an adapter oligonucleotide to the 3' ends of single DNA strands, the synthesis of a complementary strand using a DNA polymerase and the addition of a 5' adapter via blunt-end ligation. The efficiency of library preparation is determined individually for each sample using a spike-in oligonucleotide. The whole workflow, including library preparation, quantification and amplification, requires two work days for up to 16 libraries. Alternatively, we provide documentation and electronic protocols enabling automated library preparation of 96 samples in parallel on a Bravo NGS Workstation (Agilent Technologies). After library preparation, molecules with uninformative short inserts (shorter than ~30-35 base pairs) can be removed by polyacrylamide gel electrophoresis if desired.

A ligation-based single-stranded library preparation method to analyze cell-free DNA and synthetic oligos

Background

Cell-free DNA (cfDNA), present in circulating blood plasma, contains information about prenatal health, organ transplant reception, and cancer presence and progression. Originally developed for the genomic analysis of highly degraded ancient DNA, single-stranded DNA (ssDNA) library preparation methods are gaining popularity in the field of cfDNA analysis due to their efficiency and ability to convert short, fragmented DNA into sequencing libraries without altering DNA ends. However, current ssDNA methods are costly and time-consuming.

Results

Here we present an efficient ligation-based single-stranded library preparation method that is engineered to produce complex libraries in under 2.5 h from as little as 1 nanogram of input DNA without alteration to the native ends of template molecules. Our method, called Single Reaction Single-stranded LibrarY or SRSLY, ligates uniquely designed Next-Generation Sequencing (NGS) adapters in a one-step combined phosphorylation/ligation reaction that foregoes end-polishing. Using synthetic DNA oligos and cfDNA, we demonstrate the efficiency and utility of this approach and compare with existing double-stranded and single-stranded approaches for library generation. Finally, we demonstrate that cfDNA NGS data generated from SRSLY can be used to analyze DNA fragmentation patterns to deduce nucleosome positioning and transcription factor binding.

Conclusions

SRSLY is a versatile tool for converting short and fragmented DNA molecules, like cfDNA fragments, into sequencing libraries while retaining native lengths and ends.

From Genetics to Epigenetics, Roles of Epigenetics in Inflammatory Bowel Disease

Inflammatory bowel disease (IBD) is a destructive, recurrent, and heterogeneous disease. Its detailed pathogenesis is still unclear, although available evidence supports that IBD is caused by a complex interplay between genetic predispositions, environmental factors, and aberrant immune responses. Recent breakthroughs with regard to its genetics have offered valuable insights into the sophisticated genetic basis, but the identified genetic factors only explain a small part of overall disease variance. It is becoming increasingly apparent that epigenetic factors can mediate the interaction between genetics and environment, and play a fundamental role in the pathogenesis of IBD. This review outlines recent genetic and epigenetic discoveries in IBD, with a focus on the roles of epigenetics in disease susceptibility, activity, behavior and colorectal cancer (CRC), and their potential translational applications.

High-throughput methods in aptamer discovery and analysis

Aptamers are small, functional nucleic acids that bind a variety of targets, often with high specificity and affinity. Genomic aptamers constitute the ligand-binding domains of riboswitches, whereas synthetic aptamers find applications as diagnostic and therapeutic tools, and as ligand-binding domains of regulatory RNAs in synthetic biology. Discovery and characterization of aptamers has been limited by a lack of high-throughput approaches that uncover the target-binding domains and the biochemical properties of individual sequences. With the advent of high-throughput sequencing, large-scale analysis of in vitro selected populations of aptamers (and catalytic nucleic acids, such as ribozymes and DNAzmes) became possible. In recent years the development of new experimental approaches and software tools has led to significant streamlining of the selection-pool analysis. This article provides an overview of post-selection data analysis and describes high-throughput methods that facilitate rapid discovery and biochemical characterization of aptamers.

"Same difference": comprehensive evaluation of four DNA methylation measurement platforms

Background

DNA methylation in CpG context is fundamental to the epigenetic regulation of gene expression in higher eukaryotes. Changes in methylation patterns are implicated in many diseases, cellular differentiation, imprinting, and other biological processes. Techniques that enrich for biologically relevant genomic regions with high CpG content are desired, since, depending on the size of an organism’s methylome, the depth of sequencing required to cover all CpGs can be prohibitively expensive. Currently, restriction enzyme-based reduced representation bisulfite sequencing and its modified protocols are widely used to study methylation differences. Recently, Agilent Technologies, Roche NimbleGen, and Illumina have ventured to both reduce sequencing costs and capture CpGs of known biological relevance by marketing in-solution custom-capture hybridization platforms. We aimed to evaluate the similarities and differences of these four methods considering each platform targets approximately 10–13% of the human methylome.

Results

Overall, the regions covered per platform were as expected: targeted capture-based methods covered > 95% of their designed regions, whereas the restriction enzyme-based method covered > 70% of the expected fragments. While the total number of CpG loci shared by all methods was low, ~ 24% of any platform, the methylation levels of CpGs covered by all platforms were concordant. Annotation of CpG loci with genomic features revealed roughly the same proportions of feature annotations across the four platforms. Targeted capture methods comprise similar types and coverage of annotations and, relative to the targeted methods, the restriction enzyme method covers fewer promoters (~ 9%), CpG shores (~ 8%) and unannotated loci (~ 11%).

Conclusions

Although all methods are largely consistent in terms of covered CpG loci, the commercially available capture methods result in covering nearly all CpG sites in their target regions with few off-target loci and covering similar proportions of annotated CpG loci, the restriction-based enrichment results in more off-target and unannotated CpG loci. Quality of DNA is very important for restriction-based enrichment and starting material can be low. Conversely, quality of the starting material is less important for capture methods, and at least twice the amount of starting material is required. Pricing is marginally less for restriction-based enrichment, and the number of samples that can be prepared is not restricted to the number of capture reactions a kit supports. However, the advantage of capture libraries is the ability to custom design areas of interest. The choice of the technique would be decided by the number of samples, the quality and quantity of DNA available and the biological areas of interest since comparable data are obtained from all platforms.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0190-4) contains supplementary material, which is available to authorized users.

Data quality of whole genome bisulfite sequencing on Illumina platforms

The powerful HiSeq X sequencers with their patterned flowcell technology and fast turnaround times are instrumental for many large-scale genomic and epigenomic studies. However, assessment of DNA methylation by sodium bisulfite treatment results in sequencing libraries of low diversity, which may impact data quality and yield. In this report we assess the quality of WGBS data generated on the HiSeq X system in comparison with data generated on the HiSeq 2500 system and the newly released NovaSeq system. We report a systematic issue with low basecall quality scores assigned to guanines in the second read of WGBS when using certain Real Time Analysis (RTA) software versions on the HiSeq X sequencer, reminiscent of an issue that was previously reported with certain HiSeq 2500 software versions. However, with the HD.3.4.0 /RTA 2.7.7 software upgrade for the HiSeq X system, we observed an overall improved quality and yield of the WGBS data generated, which in turn empowers cost-effective and high quality DNA methylation studies.

Capturing the 'ome': the expanding molecular toolbox for RNA and DNA library construction

All sequencing experiments and most functional genomics screens rely on the generation of libraries to comprehensively capture pools of targeted sequences. In the past decade especially, driven by the progress in the field of massively parallel sequencing, numerous studies have comprehensively assessed the impact of particular manipulations on library complexity and quality, and characterized the activities and specificities of several key enzymes used in library construction. Fortunately, careful protocol design and reagent choice can substantially mitigate many of these biases, and enable reliable representation of sequences in libraries. This review aims to guide the reader through the vast expanse of literature on the subject to promote informed library generation, independent of the application.

Comparison of whole-genome bisulfite sequencing library preparation strategies identifies sources of biases affecting DNA methylation data

Background

Whole-genome bisulfite sequencing (WGBS) is becoming an increasingly accessible technique, used widely for both fundamental and disease-oriented research. Library preparation methods benefit from a variety of available kits, polymerases and bisulfite conversion protocols. Although some steps in the procedure, such as PCR amplification, are known to introduce biases, a systematic evaluation of biases in WGBS strategies is missing.

Results

We perform a comparative analysis of several commonly used pre- and post-bisulfite WGBS library preparation protocols for their performance and quality of sequencing outputs. Our results show that bisulfite conversion per se is the main trigger of pronounced sequencing biases, and PCR amplification builds on these underlying artefacts. The majority of standard library preparation methods yield a significantly biased sequence output and overestimate global methylation. Importantly, both absolute and relative methylation levels at specific genomic regions vary substantially between methods, with clear implications for DNA methylation studies.

Conclusions

We show that amplification-free library preparation is the least biased approach for WGBS. In protocols with amplification, the choice of bisulfite conversion protocol or polymerase can significantly minimize artefacts. To aid with the quality assessment of existing WGBS datasets, we have integrated a bias diagnostic tool in the Bismark package and offer several approaches for consideration during the preparation and analysis of WGBS datasets.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1408-2) contains supplementary material, which is available to authorized users.

SALP, a new single-stranded DNA library preparation method especially useful for the high-throughput characterization of chromatin openness states

BACKGROUND

Next-generation sequencing (NGS) is fundamental to the current biological and biomedical research. Construction of sequencing library is a key step of NGS. Therefore, various library construction methods have been explored. However, the current methods are still limited by some shortcomings.

RESULTS

This study developed a new NGS library construction method, Single strand Adaptor Library Preparation (SALP), by using a novel single strand adaptor (SSA). SSA is a double-stranded oligonucleotide with a 3' overhang of 3 random nucleotides, which can be efficiently ligated to the 3' end of single strand DNA by T4 DNA ligase. SALP can be started with any denatured DNA fragments such as those sheared by Tn5 tagmentation, enzyme digestion and sonication. When started with Tn5-tagmented chromatin, SALP can overcome a key limitation of ATAC-seq and become a high-throughput NGS library construction method, SALP-seq, which can be used to comparatively characterize the chromatin openness state of multiple cells unbiasly. In this way, this study successfully characterized the comparative chromatin openness states of four different cell lines, including GM12878, HepG2, HeLa and 293T, with SALP-seq. Similarly, this study also successfully characterized the chromatin openness states of HepG2 cells with SALP-seq by using 10⁵ to 500 cells.

CONCLUSIONS

This study developed a new NGS library construction method, SALP, by using a novel kind of single strand adaptor (SSA), which should has wide applications in the future due to its unique performance.

Epigenetics in pediatric acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common malignancy in children. ALL arises from the malignant transformation of progenitor B- and T-cells in the bone marrow into leukemic cells, but the mechanisms underlying this transformation are not well understood. Recent technical advances and decreasing costs of methods for high-throughput DNA sequencing and SNP genotyping have stimulated systematic studies of the epigenetic changes in leukemic cells from pediatric ALL patients. The results emerging from these studies are increasing our understanding of the epigenetic component of leukemogenesis and have demonstrated the potential of DNA methylation as a biomarker for lineage and subtype classification, prognostication, and disease progression in ALL. In this review, we provide a concise examination of the epigenetic studies in ALL, with a focus on DNA methylation and mutations perturbing genes involved in chromatin modification, and discuss the future role of epigenetic analyses in research and clinical management of ALL.